Showing 1–48 of 48 results for author: Ortiz, M

Thermalized and mixed meanfield ADP potentials for magnesium hydrides

Authors: M. Molinos, M. Ortiz, M. P. Ariza

Abstract: We develop meanfield approximation and numerical quadrature schemes for the evaluation of Angular-Dependent interatomic Potentials (ADPs) for magnesium and magnesium hydrides at finite temperature (thermalization) and arbitrary atomic molar fractions (mixing) within a non-equilibrium statistical mechanical framework and derive local equilibrium relations. We numerically verify and experimentally v… ▽ More We develop meanfield approximation and numerical quadrature schemes for the evaluation of Angular-Dependent interatomic Potentials (ADPs) for magnesium and magnesium hydrides at finite temperature (thermalization) and arbitrary atomic molar fractions (mixing) within a non-equilibrium statistical mechanical framework and derive local equilibrium relations. We numerically verify and experimentally validate the accuracy and fidelity of the resulting thermalized/mixed ADPs (TADPs) by means of selected numerical tests including free entropy, heat capacity, thermal expansion, molar volumes, equation of state and elastic constants. We show that the local equilibrium properties predicted by TADPs agree closely with those computed directly from ADP by means of Molecular Dynamics (MD). △ Less

Submitted 13 July, 2024; originally announced July 2024.

Comments: 8 figures

Fracture and size effect in mechanical metamaterials

Authors: J. Ulloa, M. P. Ariza, J. E. Andrade, M. Ortiz

Abstract: We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be char… ▽ More We resort to variational methods to evaluate the asymptotic behavior of fine metamaterials as a function of cell size. To zeroth order, the metamaterial behaves as a micropolar continuum with both displacement and rotation degrees of freedom, but exhibits linear-elastic fracture mechanics scaling and therefore no size effect. To higher order, the overall energetics of the metastructure can be characterized explicitly in terms of the solution of the zeroth-order continuum problem by the method of Γ-expansion. We present explicit expressions of the second-order correction for octet frames. As an application, we evaluate the compliance of double-cantilever octet specimens to second order and use the result to elucidate the dependence of the apparent toughness of the specimen on cell size. The analysis predicts the discreteness of the metamaterial lattice to effectively shield the crack-tip, a mechanism that we term lattice shielding. The theory specifically predicts anti-shielding, i. e., coarser is weaker, in agreement with recent experimental observations. △ Less

Submitted 26 June, 2024; originally announced July 2024.

The energy-stepping Monte Carlo method: an exactly symmetry-preserving, a Hamiltonian Monte Carlo method with a 100% acceptance ratio

Authors: Ignacio Romero, Michael Ortiz

Abstract: We introduce the energy-stepping Monte Carlo (ESMC) method, a Markov chain Monte Carlo (MCMC) algorithm based on the conventional dynamical interpretation of the proposal stage but employing an energy-stepping integrator. The energy-stepping integrator is quasi-explicit, symplectic, energy-conserving, and symmetry-preserving. As a result of the exact energy conservation of energy-stepping integrat… ▽ More We introduce the energy-stepping Monte Carlo (ESMC) method, a Markov chain Monte Carlo (MCMC) algorithm based on the conventional dynamical interpretation of the proposal stage but employing an energy-stepping integrator. The energy-stepping integrator is quasi-explicit, symplectic, energy-conserving, and symmetry-preserving. As a result of the exact energy conservation of energy-stepping integrators, ESMC has a 100\%\ acceptance ratio of the proposal states. Numerical tests provide empirical evidence that ESMC affords a number of additional benefits: the Markov chains it generates have weak autocorrelation, it has the ability to explore distant characteristic sets of the sampled probability distribution and it yields smaller errors than chains sampled with Hamiltonian Monte Carlo (HMC) and similar step sizes. Finally, ESMC benefits from the exact symmetry conservation properties of the energy-stepping integrator when sampling from potentials with built-in symmetries, whether explicitly known or not. △ Less

Submitted 12 December, 2023; originally announced December 2023.

Towards Quantum Computational Mechanics

Authors: Burigede Liu, Michael Ortiz, Fehmi Cirak

Abstract: The advent of quantum computers, operating on entirely different physical principles and abstractions from those of classical digital computers, sets forth a new computing paradigm that can potentially result in game-changing efficiencies and computational performance. Specifically, the ability to simultaneously evolve the state of an entire quantum system leads to quantum parallelism and interfer… ▽ More The advent of quantum computers, operating on entirely different physical principles and abstractions from those of classical digital computers, sets forth a new computing paradigm that can potentially result in game-changing efficiencies and computational performance. Specifically, the ability to simultaneously evolve the state of an entire quantum system leads to quantum parallelism and interference. Despite these prospects, opportunities to bring quantum computing to bear on problems of computational mechanics remain largely unexplored. In this work, we demonstrate how quantum computing can indeed be used to solve representative volume element (RVE) problems in computational homogenisation with polylogarithmic complexity of~$ \mathcal{O}((\log N)^c)$, compared to~$\mathcal{O}(N^c)$ in classical computing. Thus, our quantum RVE solver attains exponential acceleration with respect to classical solvers, bringing concurrent multiscale computing closer to practicality. The proposed quantum RVE solver combines conventional algorithms such as a fixed-point iteration for a homogeneous reference material and the Fast Fourier Transform (FFT). However, the quantum computing reformulation of these algorithms requires a fundamental paradigm shift and a complete rethinking and overhaul of the classical implementation. We employ or develop several techniques, including the Quantum Fourier Transform (QFT), quantum encoding of polynomials, classical piecewise Chebyshev approximation of functions and an auxiliary algorithm for implementing the fixed-point iteration and show that, indeed, an efficient implementation of RVE solvers on quantum computers is possible. We additionally provide theoretical proofs and numerical evidence confirming the anticipated~$ \mathcal{O} \left ((\log N)^c \right) $ complexity of the proposed solver. △ Less

Submitted 7 January, 2024; v1 submitted 6 December, 2023; originally announced December 2023.

A Simple Quantitative Model of Neuromodulation. Part I: Ion Flow Through Neural Ion Channels

Authors: Linda Werneck, Mertcan Han, Erdost Yildiz, Marc-André Keip, Metin Sitti, Michael Ortiz

Abstract: We develop a simple model of ionic current through neuronal membranes as a function of membrane potential and extracellular ion concentration. The model combines a simplified Poisson-Nernst-Planck (PNP) model of ion transport through individual mechanosensitive ion channels with channel activation functions calibrated from ad hoc in-house experimental data. The simplified PNP model is validated ag… ▽ More We develop a simple model of ionic current through neuronal membranes as a function of membrane potential and extracellular ion concentration. The model combines a simplified Poisson-Nernst-Planck (PNP) model of ion transport through individual mechanosensitive ion channels with channel activation functions calibrated from ad hoc in-house experimental data. The simplified PNP model is validated against bacterial Gramicidin A ion channel data. The calibrated model accounts for the transport of calcium, sodium, potassium, and chloride and exhibits remarkable agreement with the experimentally measured current-voltage curves for the differentiated human neural cells. All relevant data and code related to the ion flow models are available at DaRUS. △ Less

Submitted 4 September, 2023; originally announced September 2023.

Comments: 17 pages, 8 figures

Modeling hard-soft block copolymers as a liquid crystalline polymer

Authors: M. Manav, M. Ponga, M. Ortiz

Abstract: We report a new computational approach to model hard-soft block copolymers like polyurea as a liquid crystalline polymer to understand their microstructural evolution due to mechanical loading. The resulting microstructure closely resembles the microstructure observed in polyurea. The stress-strain relations in uniaxial compression and tension loading obtained from the model are also in close quan… ▽ More We report a new computational approach to model hard-soft block copolymers like polyurea as a liquid crystalline polymer to understand their microstructural evolution due to mechanical loading. The resulting microstructure closely resembles the microstructure observed in polyurea. The stress-strain relations in uniaxial compression and tension loading obtained from the model are also in close quantitative agreement with the experimental data for polyurea. We use the model to elucidate the evolution of the hard and the soft domains during loading, which is consistent with the experimental measurements characterizing microstructural evolution in polyurea. △ Less

Submitted 11 May, 2023; originally announced May 2023.

Homogenizing elastic properties of large digital rock images by combining CNN with hierarchical homogenization method

Authors: Rasool Ahmad, Mingliang Liu, Michael Ortiz, Tapan Mukerji, Wei Cai

Abstract: Determining effective elastic properties of rocks from their pore-scale digital images is a key goal of digital rock physics (DRP). Direct numerical simulation (DNS) of elastic behavior, however, incurs high computational cost; and surrogate machine learning (ML) model, particularly convolutional neural network (CNN), show promises to accelerate homogenization process. 3D CNN models, however, are… ▽ More Determining effective elastic properties of rocks from their pore-scale digital images is a key goal of digital rock physics (DRP). Direct numerical simulation (DNS) of elastic behavior, however, incurs high computational cost; and surrogate machine learning (ML) model, particularly convolutional neural network (CNN), show promises to accelerate homogenization process. 3D CNN models, however, are unable to handle large images due to memory issues. To address this challenge, we propose a novel method that combines 3D CNN with hierarchical homogenization method (HHM). The surrogate 3D CNN model homogenizes only small subimages, and a DNS is used to homogenize the intermediate image obtained by assembling small subimages. The 3D CNN model is designed to output the homogenized elastic constants within the Hashin-Shtrikman (HS) bounds of the input images. The 3D CNN model is first trained on data comprising equal proportions of five sandstone (quartz mineralogy) images, and, subsequently, fine-tuned for specific rocks using transfer learning. The proposed method is applied to homogenize the rock images of size 300x300x300 and 600x600x600 voxels, and the predicted homogenized elastic moduli are shown to agree with that obtained from the brute-force DNS. The transferability of the trained 3D CNN model (using transfer learning) is further demonstrated by predicting the homogenized elastic moduli of a limestone rock with calcite mineralogy. The surrogate 3D CNN model in combination with the HHM is thus shown to be a promising tool for the homogenization of large 3D digital rock images and other random media △ Less

Submitted 10 May, 2023; originally announced May 2023.

Comments: 16 pages, 13 figures

An optimal-transport finite-particle method for mass diffusion

Authors: Anna Pandolfi, Laurent Stainier, Michael Ortiz

Abstract: We formulate a class of velocity-free finite-particle methods for mass transport problems based on a time-discrete incremental variational principle that combines entropy and the cost of particle transport, as measured by the Wasserstein metric. The incremental functional is further spatially discretized into finite particles, i.e., particles characterized by a fixed spatial profile of finite widt… ▽ More We formulate a class of velocity-free finite-particle methods for mass transport problems based on a time-discrete incremental variational principle that combines entropy and the cost of particle transport, as measured by the Wasserstein metric. The incremental functional is further spatially discretized into finite particles, i.e., particles characterized by a fixed spatial profile of finite width, each carrying a fixed amount of mass. The motion of the particles is then governed by a competition between the cost of transport, that aims to keep the particles fixed, and entropy maximization, that aims to spread the particles so as to increase the entropy of the system. We show how the optimal width of the particles can be determined variationally by minimization of the governing incremental functional. Using this variational principle, we derive optimal scaling relations between the width of the particles, their number and the size of the domain. We also address matters of implementation including the acceleration of the computation of diffusive forces by exploiting the Gaussian decay of the particle profiles and by instituting fast nearest-neighbor searches. We demonstrate the robustness and versatility of the finite-particle method by means of a test problem concerned with the injection of mass into a sphere. There test results demonstrate the meshless character of the method in any spatial dimension, its ability to redistribute mass particles and follow their evolution in time, its ability to satisfy flux boundary conditions for general domains based solely on a distance function, and its robust convergence characteristics. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Interatomic-Potential-Free, Data-Driven Molecular Dynamics

Authors: J. Bulin, J. Hamaekers, M. P. Ariza, M. Ortiz

Abstract: We present a Data-Driven (DD) paradigm that enables molecular dynamics calculations to be performed directly from sampled force-field data such as obtained, e.g., from ab initio calculations, thereby eschewing the conventional step of modeling the data by empirical interatomic potentials entirely. The data required by the DD solvers consists of local atomic configurations and corresponding atomic… ▽ More We present a Data-Driven (DD) paradigm that enables molecular dynamics calculations to be performed directly from sampled force-field data such as obtained, e.g., from ab initio calculations, thereby eschewing the conventional step of modeling the data by empirical interatomic potentials entirely. The data required by the DD solvers consists of local atomic configurations and corresponding atomic forces and is, therefore, fundamental, i.e., it is not beholden to any particular model. The resulting DD solvers, including a fully explicit DD-Verlet algorithm, are provably convergent and exhibit robust convergence with respect to the data in selected test cases. We present an example of application to C60 buckminsterfullerenes that showcases the feasibility, range and scope of the DD molecular dynamics paradigm. △ Less

Submitted 5 August, 2022; originally announced August 2022.

Comments: 22 pages, 10 figures

Computation of effective elastic moduli of rocks using hierarchical homogenization

Authors: Rasool Ahmad, Mingliang Liu, Michael Ortiz, Tapan Mukerji, Wei Cai

Abstract: This work focuses on computing the homogenized elastic properties of rocks from 3D micro-computed-tomography (micro-CT) scanned images. The accurate computation of homogenized properties of rocks, archetypal random media, requires both resolution of intricate underlying microstructure and large field of view, resulting in huge micro-CT images. Homogenization entails solving the local elasticity pr… ▽ More This work focuses on computing the homogenized elastic properties of rocks from 3D micro-computed-tomography (micro-CT) scanned images. The accurate computation of homogenized properties of rocks, archetypal random media, requires both resolution of intricate underlying microstructure and large field of view, resulting in huge micro-CT images. Homogenization entails solving the local elasticity problem computationally which can be prohibitively expensive for a huge image. To mitigate this problem, we use a renormalization method inspired scheme, the hierarchical homogenization method, where a large image is partitioned into smaller subimages. The individual subimages are separately homogenized using periodic boundary conditions, and then assembled into a much smaller intermediate image. The intermediate image is again homogenized, subject to the periodic boundary condition, to find the final homogenized elastic constant of the original image. An FFT-based elasticity solver is used to solve the associated periodic elasticity problem. The error in the homogenized elastic constant is empirically shown to follow a power law scaling with exponent -1 with respect to the subimage size across all five microstructures of rocks. We further show that the inclusion of surrounding materials during the homogenization of the small subimages reduces error in the final homogenized elastic moduli while still respecting the power law with the exponent of -1. This power law scaling is then exploited to determine a better approximation of the large heterogeneous microstructures based on Richardson extrapolation △ Less

Submitted 3 August, 2022; originally announced August 2022.

Journal ref: Journal of the Mechanics and Physics of Solids, 2023

Model-Free Data-Driven Viscoelasticity in the Frequency Domain

Authors: Hossein Salahshoor, Michael Ortiz

Abstract: We develop a Data-Driven framework for the simulation of wave propagation in viscoelastic solids directly from dynamic testing material data, including data from Dynamic Mechanical Analysis (DMA), nano-indentation, Dynamic Shear Testing (DST) and Magnetic Resonance Elastography (MRE), without the need for regression or material modeling. The problem is formulated in the frequency domain and the me… ▽ More We develop a Data-Driven framework for the simulation of wave propagation in viscoelastic solids directly from dynamic testing material data, including data from Dynamic Mechanical Analysis (DMA), nano-indentation, Dynamic Shear Testing (DST) and Magnetic Resonance Elastography (MRE), without the need for regression or material modeling. The problem is formulated in the frequency domain and the method of solution seeks to minimize a distance between physically admissible histories of stress and strain, in the sense of compatibility and equilibrium, and the material data. We metrize the space of histories by means of the flat-norm of their Fourier transform, which allows consideration of infinite wave trains such as harmonic functions. Another significant advantage of the flat norm is that it allows the response of the system at one frequency to be inferred from data at nearby frequencies. We demonstrate and verify the approach by means of two test cases, a polymeric truss structure characterized by DMA data and a 3D soft gel sample characterized by MRE data. The examples demonstrate the ease of implementation of the Data-Driven scheme within conventional commercial codes and its robust convergence properties, both with respect to the solver and the data. △ Less

Submitted 13 May, 2022; originally announced May 2022.

Hierarchical multiscale quantification of material uncertainty

Authors: Burigede Liu, Xingsheng Sun, Kaushik Bhattacharya, Michael Ortiz

Abstract: The macroscopic behavior of many materials is complex and the end result of mechanisms that operate across a broad range of disparate scales. An imperfect knowledge of material behavior across scales is a source of epistemic uncertainty of the overall material behavior. However, assessing this uncertainty is difficult due to the complex nature of material response and the prohibitive computational… ▽ More The macroscopic behavior of many materials is complex and the end result of mechanisms that operate across a broad range of disparate scales. An imperfect knowledge of material behavior across scales is a source of epistemic uncertainty of the overall material behavior. However, assessing this uncertainty is difficult due to the complex nature of material response and the prohibitive computational cost of integral calculations. In this paper, we exploit the multiscale and hierarchical nature of material response to develop an approach to quantify the overall uncertainty of material response without the need for integral calculations. Specifically, we bound the uncertainty at each scale and then combine the partial uncertainties in a way that provides a bound on the overall or integral uncertainty. The bound provides a conservative estimate on the uncertainty. Importantly, this approach does not require integral calculations that are prohibitively expensive. We demonstrate the framework on the problem of ballistic impact of a polycrystalline magnesium plate. Magnesium and its alloys are of current interest as promising light-weight structural and protective materials. Finally, we remark that the approach can also be used to study the sensitivity of the overall response to particular mechanisms at lower scales in a materials-by-design approach. △ Less

Submitted 4 February, 2021; originally announced February 2021.

Nonequilibrium thermomechanics of Gaussian phase packet crystals: application to the quasistatic quasicontinuum method

Authors: Prateek Gupta, Michael Ortiz, Dennis M. Kochmann

Abstract: The quasicontinuum method was originally introduced to bridge across length scales -- from atomistics to significantly larger continuum scales -- thus overcoming a key limitation of classical atomic-scale simulation techniques while solely relying on atomic-scale input (in the form of interatomic potentials). An associated challenge lies in bridging across time scales to overcome the time scale li… ▽ More The quasicontinuum method was originally introduced to bridge across length scales -- from atomistics to significantly larger continuum scales -- thus overcoming a key limitation of classical atomic-scale simulation techniques while solely relying on atomic-scale input (in the form of interatomic potentials). An associated challenge lies in bridging across time scales to overcome the time scale limitations of atomistics. To address the biggest challenge, bridging across both length and time scales, only a few techniques exist, and most of those are limited to conditions of constant temperature. Here, we present a new strategy for the space-time coarsening of an atomistic ensemble, which introduces thermomechanical coupling. We investigate the quasistatics and dynamics of a crystalline solid described as a lattice of lumped correlated Gaussian phase packets occupying atomic lattice sites. By definition, phase packets account for the dynamics of crystalline lattices at finite temperature through the statistical variances of atomic momenta and positions. We show that momentum-space correlation allows for an exchange between potential and kinetic contributions to the crystal's Hamiltonian. Consequently, local adiabatic heating due to atomic site motion is captured. Moreover, within the quasistatic approximation the governing equations reduce to the minimization of thermodynamic potentials such as Helmholtz free energy (depending on the fixed variables), and they yield the local equation of state. We further discuss opportunities for describing atomic-level thermal transport using the correlated Gaussian phase packet formulation and the importance of interatomic correlations. Such a formulation offers a promising avenue for a finite-temperature non-equilibrium quasicontinuum method that may be combined with thermal transport models. △ Less

Submitted 14 April, 2021; v1 submitted 25 January, 2021; originally announced January 2021.

Comments: Content of the new submission has changed significantly along with the author list

Robust triboelectric charging of identical balloons of different radii

Authors: Francisco Vera, Rodrigo Rivera, Manuel Ortiz, Francisco Antonio Horta-Rangel

Abstract: Electrification by rubbing different materials is a well known phenomenon with a history that begun more than five centuries BC. However, simple experiments can lead to contradictory or inconsistent results and the history of this phenomena is plagued with non-intuitive results. For example, triboelectric charging by rubbing identical materials is possible. In this work we want to highlight some h… ▽ More Electrification by rubbing different materials is a well known phenomenon with a history that begun more than five centuries BC. However, simple experiments can lead to contradictory or inconsistent results and the history of this phenomena is plagued with non-intuitive results. For example, triboelectric charging by rubbing identical materials is possible. In this work we want to highlight some historical aspects of triboelectricity that could enrich the discussion of electrostatics in an undergraduate physics course. We will focus on the effect of strain on the triboelectric properties of a sample, which we think is not well known to physics teachers. We will show that it is possible to obtain robust polarities by rubbing identical rubber balloons of different radii and we will also show that this charging method can be very useful in introductory physics courses. △ Less

Submitted 18 November, 2020; originally announced November 2020.

The heterodyne sensing system for the ALPS II search for sub-eV weakly interacting particles

Authors: Ayman Hallal, Giuseppe Messineo, Mauricio Diaz Ortiz, Joseph Gleason, Harold Hollis, D. B. Tanner, Guido Mueller, Aaron Spector

Abstract: ALPS II, the Any Light Particle Search, is a second-generation Light Shining through a Wall experiment that hunts for axion-like particles. The experiment is currently transitioning from the design and construction phase to the commissioning phase, with science runs expected to start in 2021. ALPS II plans to use two different sensing schemes to confirm the potential detection of axion-like partic… ▽ More ALPS II, the Any Light Particle Search, is a second-generation Light Shining through a Wall experiment that hunts for axion-like particles. The experiment is currently transitioning from the design and construction phase to the commissioning phase, with science runs expected to start in 2021. ALPS II plans to use two different sensing schemes to confirm the potential detection of axion-like particles or to verify an upper limit on their coupling strength to two photons of $g_{aγγ}\leq2\times10^{-11}\text{GeV}^{-1}$. This paper discusses a heterodyne sensing scheme (HET) which will be the first scheme deployed to detect the regenerated light. It presents critical details of the optical layout, the length and alignment sensing scheme, design features to minimize spurious signals from stray light, as well as several control and veto channels specific to HET which are needed to commission and operate the instrument and to calibrate the detector sensitivity. △ Less

Submitted 7 February, 2021; v1 submitted 1 October, 2020; originally announced October 2020.

Comments: 12 pages, 5 figures

Report number: DESY 20-164

Design of the ALPS II Optical System

Authors: M. Diaz Ortiz, J. Gleason, H. Grote, A. Hallal, M. T. Hartman, H. Hollis, K. S. Isleif, A. James, K. Karan, T. Kozlowski, A. Lindner, G. Messineo, G. Mueller, J. H. Poeld, R. C. G. Smith, A. D. Spector, D. B. Tanner, L. -W. Wei, B. Willke

Abstract: The Any Light Particle Search II (ALPS II) is an experiment currently being built at DESY in Hamburg, Germany, that will use a light-shining-through-a-wall (LSW) approach to search for axion-like particles. ALPS II represents a significant step forward for these types of experiments as it will use 24 superconducting dipole magnets, along with dual, high-finesse, 122 m long optical cavities. This p… ▽ More The Any Light Particle Search II (ALPS II) is an experiment currently being built at DESY in Hamburg, Germany, that will use a light-shining-through-a-wall (LSW) approach to search for axion-like particles. ALPS II represents a significant step forward for these types of experiments as it will use 24 superconducting dipole magnets, along with dual, high-finesse, 122 m long optical cavities. This paper gives the first comprehensive recipe for the realization of the idea, proposed over 30 years ago, to use optical cavities before and after the wall to increase the power of the regenerated photon signal. The experiment is designed to achieve a sensitivity to the coupling between axion-like particles and photons down to g=2e-11 1/GeV for masses below 0.1 meV, more than three orders of magnitude beyond the sensitivity of previous laboratory experiments. The layout and main components that define ALPS II are discussed along with plans for reaching design sensitivity. An accompanying paper (Hallal, et al [1]) offers a more in-depth description of the heterodyne detection scheme, the first of two independent detection systems that will be implemented in ALPS II. △ Less

Submitted 21 December, 2021; v1 submitted 29 September, 2020; originally announced September 2020.

Comments: 26 pages, 5 figures

Report number: DESY 20-158

Ice penetration by a bluff-body melting probe

Authors: Kerstin Weinberg, Michael Ortiz

Abstract: We analyze the operation of melting probes as a Stefan problem for the liquid/solid interface surrounding the probe. We assume that the liquid layer is thin and, therefore, amenable to analysis by lubrication theory. The resulting Stefan problem is solvable in closed form. The solution determines the dependence of the penetration speed on the temperature differential between the probe and the surr… ▽ More We analyze the operation of melting probes as a Stefan problem for the liquid/solid interface surrounding the probe. We assume that the liquid layer is thin and, therefore, amenable to analysis by lubrication theory. The resulting Stefan problem is solvable in closed form. The solution determines the dependence of the penetration speed on the temperature differential between the probe and the surrounding ice, the size, shape and weight of the probe, the viscosity of liquid water and the thermal properties of solid ice. △ Less

Submitted 6 February, 2020; originally announced February 2020.

A Dynamical Model of Oncotripsy by Mechanical Cell Fatigue: Selective Cancer Cell Ablation by Low-Intensity Pulsed Ultrasound (LIPUS)

Authors: E. F. Schibber, D. R. Mittelstein, M. Gharib, M. G. Shapiro, P. P. Lee, M. Ortiz

Abstract: The method of oncotripsy, first proposed in [S. Heyden and M. Ortiz (2016). Oncotripsy: Targeting cancer cells selectively via resonant harmonic excitation. Journal of the Mechanics and Physics of Solids, 92:164-175], exploits aberrations in the material properties and morphology of cancerous cells in order to ablate them selectively by means of tuned low-intensity pulsed ultrasound (LIPUS). We pr… ▽ More The method of oncotripsy, first proposed in [S. Heyden and M. Ortiz (2016). Oncotripsy: Targeting cancer cells selectively via resonant harmonic excitation. Journal of the Mechanics and Physics of Solids, 92:164-175], exploits aberrations in the material properties and morphology of cancerous cells in order to ablate them selectively by means of tuned low-intensity pulsed ultrasound (LIPUS). We propose a dynamical model of oncotripsy that follows as an application of cell dynamics, statistical mechanical theory of network elasticity and 'birth-death' kinetics to describe processes of damage and repair of the cytoskeleton. We also develop a reduced dynamical model that approximates the three-dimensional dynamics of the cell and facilitates parametric studies, including sensitivity analysis and process optimization. We show that the dynamical model predicts---and provides a conceptual basis for understanding---the oncotripsy effect and other trends in the data of [D. R. Mittelstein, J. Ye, E. F. Schibber, A. Roychoudhury, L. T. Martinez, M. H. Fekrazad, M. Ortiz, P. P. Lee, M. G. Shapiro, M. Gharib (2019). Selective Ablation of Cancer Cells with Low Intensity Pulsed Ultrasound. BioRxiv] for cells in suspension, including the dependence of cell-death curves on cell and process parameters. △ Less

Submitted 27 November, 2019; originally announced November 2019.

Nonequilibrium Casimir pressures in liquids under shear

Authors: J. M. Ortiz de Zárate, T. R. Kirkpatrick, J. V. Sengers

Abstract: In stationary nonequilibrium states coupling between hydrodynamic modes causes thermal fluctuations to become long ranged inducing nonequilibrium Casimir pressures. Here we consider nonequilibrium Casimir pressures induced in liquids by a velocity gradient. Specifically, we have obtained explicit expressions for the magnitude of the shear-induced pressure enhancements in a liquid layer between two… ▽ More In stationary nonequilibrium states coupling between hydrodynamic modes causes thermal fluctuations to become long ranged inducing nonequilibrium Casimir pressures. Here we consider nonequilibrium Casimir pressures induced in liquids by a velocity gradient. Specifically, we have obtained explicit expressions for the magnitude of the shear-induced pressure enhancements in a liquid layer between two horizontal plates that complete and correct results previously presented in the literature. In contrast to nonequilibrium Casimir pressures induced by a temperature or concentration gradient, we find that in shear nonequilibrium contributions from short-range fluctuations are no longer negligible. In addition, it is noted that currently available computer simulations of model fluids in shear observe effects from molecular correlations at nanoscales that have a different physical origin and do not probe shear-induced pressures resulting from coupling of long-wavelength hydrodynamic modes. Even more importantly, we find that in actual experimental conditions, shear-induced pressure enhancements are caused by viscous heating and not by thermal velocity fluctuations. Hence, isothermal computer simulations are irrelevant for the interpretation of experimental shear-induced pressure enhancements. △ Less

Submitted 27 July, 2019; originally announced July 2019.

Comments: 20 pages, 2 figures, 2 tables. arXiv admin note: text overlap with arXiv:1804.06125

Large scale ab-initio simulations of dislocations

Authors: Mauricio Ponga, Kaushik Bhattacharya, Michael Ortiz

Abstract: We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to i… ▽ More We present a novel methodology to compute relaxed dislocations core configurations, and their energies in crystalline metallic materials using large-scale \emph{ab-intio} simulations. The approach is based on MacroDFT, a coarse-grained density functional theory method that accurately computes the electronic structure but with sub-linear scaling resulting in a tremendous reduction in cost. Due to its implementation in \emph{real-space}, MacroDFT has the ability to harness petascale resources to study materials and alloys through accurate \emph{ab-initio} calculations. Thus, the proposed methodology can be used to investigate dislocation cores and other defects where long range elastic defects play an important role, such as in dislocation cores, grain boundaries and near precipitates in crystalline materials. We demonstrate the method by computing the relaxed dislocation cores in prismatic dislocation loops and dislocation segments in magnesium (Mg). We also study the interaction energy with a line of Aluminum (Al) solutes. Our simulations elucidate the essential coupling between the quantum mechanical aspects of the dislocation core and the long range elastic fields that they generate. In particular, our quantum mechanical simulations are able to describe the logarithmic divergence of the energy in the far field as is known from classical elastic theory. In order to reach such scaling, the number of atoms in the simulation cell has to be exceedingly large, and cannot be achieved with the state-of-the-art density functional theory implementations. △ Less

Submitted 4 July, 2019; originally announced July 2019.

Hypernetwork Science via High-Order Hypergraph Walks

Authors: Sinan G. Aksoy, Cliff Joslyn, Carlos Ortiz Marrero, Brenda Praggastis, Emilie Purvine

Abstract: We propose high-order hypergraph walks as a framework to generalize graph-based network science techniques to hypergraphs. Edge incidence in hypergraphs is quantitative, yielding hypergraph walks with both length and width. Graph methods which then generalize to hypergraphs include connected component analyses, graph distance-based metrics such as closeness centrality, and motif-based measures suc… ▽ More We propose high-order hypergraph walks as a framework to generalize graph-based network science techniques to hypergraphs. Edge incidence in hypergraphs is quantitative, yielding hypergraph walks with both length and width. Graph methods which then generalize to hypergraphs include connected component analyses, graph distance-based metrics such as closeness centrality, and motif-based measures such as clustering coefficients. We apply high-order analogs of these methods to real world hypernetworks, and show they reveal nuanced and interpretable structure that cannot be detected by graph-based methods. Lastly, we apply three generative models to the data and find that basic hypergraph properties, such as density and degree distributions, do not necessarily control these new structural measurements. Our work demonstrates how analyses of hypergraph-structured data are richer when utilizing tools tailored to capture hypergraph-native phenomena, and suggests one possible avenue towards that end. △ Less

Submitted 8 June, 2020; v1 submitted 26 June, 2019; originally announced June 2019.

Comments: Updated to address referee comments, to appear in EPJ Data Science

Optical Supersymmetry in the Time Domain

Authors: Carlos García-Meca, Andrés Macho Ortiz, Roberto Llorente Sáez

Abstract: Originally emerged within the context of string and quantum field theory, and later fruitfully extrapolated to photonics, the algebraic transformations of quantum-mechanical supersymmetry were conceived in the space realm. Here, we introduce a paradigm shift, demonstrating that Maxwell's equations also possess an underlying supersymmetry in the time domain. As a result, we obtain a simple analytic… ▽ More Originally emerged within the context of string and quantum field theory, and later fruitfully extrapolated to photonics, the algebraic transformations of quantum-mechanical supersymmetry were conceived in the space realm. Here, we introduce a paradigm shift, demonstrating that Maxwell's equations also possess an underlying supersymmetry in the time domain. As a result, we obtain a simple analytic relation between the scattering coefficients of a large variety of time-varying optical systems and uncover a wide new class of reflectionless, three dimensional, all-dielectric, isotropic, omnidirectional, polarization-independent, non-complex media. Temporal supersymmetry is also shown to arise in dispersive media supporting temporal bound states, which allows engineering their momentum spectra and dispersive properties. These unprecedented features define a promising design platform for free-space and integrated photonics, enabling the creation of a number of novel reconfigurable reflectionless devices, such as frequency-selective, polarization-independent and omnidirectional invisible materials, compact frequency-independent phase shifters, broadband isolators, and versatile pulse-shape transformers. △ Less

Submitted 24 March, 2019; originally announced March 2019.

Journal ref: Nature Communications 11, 813 (2020)

Definition of frame-invariant Soret coefficients for ternary mixtures

Authors: José M. Ortiz de Zárate

Abstract: The definition of the Soret coefficient of a binary mixture includes a concentration prefactor, $x(1-x)$ when mol fraction $x$ is used, or $w(1-w)$ when mass fraction $w$ is used. In this paper the physical reasons behind this choice are reviewed, emphasizing that the use of these prefactors makes the Soret coefficient invariant upon change in the reference frame, either mass or molar. Then, it is… ▽ More The definition of the Soret coefficient of a binary mixture includes a concentration prefactor, $x(1-x)$ when mol fraction $x$ is used, or $w(1-w)$ when mass fraction $w$ is used. In this paper the physical reasons behind this choice are reviewed, emphasizing that the use of these prefactors makes the Soret coefficient invariant upon change in the reference frame, either mass or molar. Then, it is shown how this invariance property can be extended to ternary mixtures by using an appropriate concentration prefactor in matrix form. The work is completed with some considerations of general non-isothermal diffusion fluxes, binary limits of the concentration triangle, selection of the dependent concentration in a ternary mixture, use of molar concentrations and, finally, extension to multi-component mixtures. △ Less

Submitted 7 January, 2019; v1 submitted 25 October, 2018; originally announced October 2018.

Comments: 5 pages, no figures

Model-Free Data-Driven Inelasticity

Authors: Robert Eggersmann, Trenton Kirchdoerfer, Stefanie Reese, Laurent Stainier, Michael Ortiz

Abstract: We extend the Data-Driven formulation of problems in elasticity of Kirchdoerfer and Ortiz (2016) to inelasticity. This extension differs fundamentally from Data-Driven problems in elasticity in that the material data set evolves in time as a consequence of the history dependence of the material. We investigate three representational paradigms for the evolving material data sets: i) materials with… ▽ More We extend the Data-Driven formulation of problems in elasticity of Kirchdoerfer and Ortiz (2016) to inelasticity. This extension differs fundamentally from Data-Driven problems in elasticity in that the material data set evolves in time as a consequence of the history dependence of the material. We investigate three representational paradigms for the evolving material data sets: i) materials with memory, i.e., conditioning the material data set to the past history of deformation; ii) differential materials, i.e., conditioning the material data set to short histories of stress and strain; and iii) history variables, i.e., conditioning the material data set to ad hoc variables encoding partial information about the history of stress and strain. We also consider combinations of the three paradigms thereof and investigate their ability to represent the evolving data sets of different classes of inelastic materials, including viscoelasticity, viscoplasticity and plasticity. We present selected numerical examples that demonstrate the range and scope of Data-Driven inelasticity and the numerical performance of implementations thereof. △ Less

Submitted 10 September, 2018; v1 submitted 31 August, 2018; originally announced August 2018.

Comments: Minor revisions: affiliations, acknowledgments

ACM Class: J.2

Learning Parameters and Constitutive Relationships with Physics Informed Deep Neural Networks

Authors: Alexandre M. Tartakovsky, Carlos Ortiz Marrero, Paris Perdikaris, Guzel D. Tartakovsky, David Barajas-Solano

Abstract: We present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partial… ▽ More We present a physics informed deep neural network (DNN) method for estimating parameters and unknown physics (constitutive relationships) in partial differential equation (PDE) models. We use PDEs in addition to measurements to train DNNs to approximate unknown parameters and constitutive relationships as well as states. The proposed approach increases the accuracy of DNN approximations of partially known functions when a limited number of measurements is available and allows for training DNNs when no direct measurements of the functions of interest are available. We employ physics informed DNNs to estimate the unknown space-dependent diffusion coefficient in a linear diffusion equation and an unknown constitutive relationship in a non-linear diffusion equation. For the parameter estimation problem, we assume that partial measurements of the coefficient and states are available and demonstrate that under these conditions, the proposed method is more accurate than state-of-the-art methods. For the non-linear diffusion PDE model with a fully unknown constitutive relationship (i.e., no measurements of constitutive relationship are available), the physics informed DNN method can accurately estimate the non-linear constitutive relationship based on state measurements only. Finally, we demonstrate that the proposed method remains accurate in the presence of measurement noise. △ Less

Submitted 17 August, 2018; v1 submitted 9 August, 2018; originally announced August 2018.

Ultra-short Pulse Propagation in Nonlinear Optics

Authors: Andrés Macho Ortiz, Roberto Llorente Sáez

Abstract: In this work, we perform a detailed review about the theoretical modeling of the optical propagation of ultra-short pulses in different scenarios of nonlinear optics. In particular, we focus our efforts on optical fibers and uniform media, with special attention on the white light continuum generation by using gases. In this work, we perform a detailed review about the theoretical modeling of the optical propagation of ultra-short pulses in different scenarios of nonlinear optics. In particular, we focus our efforts on optical fibers and uniform media, with special attention on the white light continuum generation by using gases. △ Less

Submitted 27 June, 2018; originally announced June 2018.

Comments: 18 pages, 2 figures

A Line-free Method of Monopoles for 3D Dislocation Dynamics

Authors: A. Deffo, M. P. Ariza, M. Ortiz

Abstract: We develop an approximation scheme for three-dimensional dislocation dynamics in which the dislocation line density is concentrated at points, or monopoles. Every monopole carries a Burgers vector and an element of line. The monopoles move according to mobility kinetics driven by elastic and applied forces. The divergence constraint, expressing the requirement that the monopoles approximate a boun… ▽ More We develop an approximation scheme for three-dimensional dislocation dynamics in which the dislocation line density is concentrated at points, or monopoles. Every monopole carries a Burgers vector and an element of line. The monopoles move according to mobility kinetics driven by elastic and applied forces. The divergence constraint, expressing the requirement that the monopoles approximate a boundary, is enforced weakly. The fundamental difference with traditional approximation schemes based on segments is that in the present approach an explicit linear connectivity, or 'sequence', between the monopoles need not be defined. Instead, the monopoles move as an unstructured point set subject to the weak divergence constraint. In this sense, the new paradigm is 'line-free', i. e., it sidesteps the need to track dislocation lines. This attribute offers significant computational advantages in terms of simplicity, robustness and efficiency, as demonstrated by means of selected numerical examples. △ Less

Submitted 8 June, 2018; originally announced June 2018.

Shear-induced Casimir forces in liquid layers

Authors: J. M. Ortiz de Zárate, T. R. Kirkpatrick, J. V. Sengers

Abstract: In stationary nonequilibrium states a coupling between hydrodynamic modes causes thermal fluctuations to become long ranged inducing nonequilibrium Casimir forces or pressures. Here we consider nonequilibrium Casimir pressures induced in liquids by a velocity gradient. Specifically, we have obtained explicit expressions for the magnitude of the shear-induced pressure enhancement in a liquid layer… ▽ More In stationary nonequilibrium states a coupling between hydrodynamic modes causes thermal fluctuations to become long ranged inducing nonequilibrium Casimir forces or pressures. Here we consider nonequilibrium Casimir pressures induced in liquids by a velocity gradient. Specifically, we have obtained explicit expressions for the magnitude of the shear-induced pressure enhancement in a liquid layer between two horizontal plates that complete and correct results previously presented in the literature. In contrast to nonequiibrium Casimir pressures induced by a temperature gradient, kinetic theory shows that nonequilibrium contributions from short-range fluctuations are no longer negligible. In addition, it is noted that computer simulations of model fluids in shear observe effects from molecular correlations at nanoscales that have a different physical origin. The idea that such computer simulations probe shear-induced pressures resulting from coupling of long-wavelength hydrodynamic modes is erroneous. △ Less

Submitted 24 January, 2019; v1 submitted 17 April, 2018; originally announced April 2018.

Comments: 8 pages. Two tables. No figures

Functional optimality of the sulcus pattern of the human brain

Authors: Stefanie Heyden, Michael Ortiz

Abstract: We develop a mathematical model of information transmission across the biological neural network of the human brain. The overall function of the brain consists of the emergent processes resulting from the spread of information through the neural network. The capacity of the brain is therefore related to the rate at which it can transmit information through the neural network. The particular transm… ▽ More We develop a mathematical model of information transmission across the biological neural network of the human brain. The overall function of the brain consists of the emergent processes resulting from the spread of information through the neural network. The capacity of the brain is therefore related to the rate at which it can transmit information through the neural network. The particular transmission model under consideration allows for information to be transmitted along multiple paths between points of the cortex. The resulting transmission rates are governed by potential theory. According to this theory, the brain has preferred and quantized transmission modes that correspond to eigenfunctions of the classical Steklov eigenvalue problem, with the reciprocal eigenvalues quantifying the corresponding transmission rates. We take the model as a basis for testing the hypothesis that the sulcus pattern of the human brain has evolved to maximize the rate of transmission of information between points in the cerebral cortex. We show that the introduction of sulci, or cuts, in an otherwise smooth domain indeed increases the overall transmission rate. We demonstrate this result by means of numerical experiments concerned with a spherical domain with a varying number of slits on its surface. △ Less

Submitted 18 January, 2018; v1 submitted 17 January, 2018; originally announced January 2018.

Data-Driven Computing in Dynamics

Authors: Trenton Kirchdoerfer, Michael Ortiz

Abstract: We formulate extensions to Data Driven Computing for both distance minimizing and entropy maximizing schemes to incorporate time integration. Previous works focused on formulating both types of solvers in the presence of static equilibrium constraints. Here formulations assign data points a variable relevance depending on distance to the solution and on maximum-entropy weighting, with distance min… ▽ More We formulate extensions to Data Driven Computing for both distance minimizing and entropy maximizing schemes to incorporate time integration. Previous works focused on formulating both types of solvers in the presence of static equilibrium constraints. Here formulations assign data points a variable relevance depending on distance to the solution and on maximum-entropy weighting, with distance minimizing schemes discussed as a special case. The resulting schemes consist of the minimization of a suitably-defined free energy over phase space subject to compatibility and a time-discretized momentum conservation constraint. The present selected numerical tests that establish the convergence properties of both types of Data Driven solvers and solutions. △ Less

Submitted 9 June, 2017; originally announced June 2017.

Comments: arXiv admin note: substantial text overlap with arXiv:1702.01574

ACM Class: G.1.8

Data Driven Computing with Noisy Material Data Sets

Authors: Trenton Kirchdoerfer, Michael Ortiz

Abstract: We formulate a Data Driven Computing paradigm, termed max-ent Data Driven Computing, that generalizes distance-minimizing Data Driven Computing and is robust with respect to outliers. Robustness is achieved by means of clustering analysis. Specifically, we assign data points a variable relevance depending on distance to the solution and on maximum-entropy estimation. The resulting scheme consists… ▽ More We formulate a Data Driven Computing paradigm, termed max-ent Data Driven Computing, that generalizes distance-minimizing Data Driven Computing and is robust with respect to outliers. Robustness is achieved by means of clustering analysis. Specifically, we assign data points a variable relevance depending on distance to the solution and on maximum-entropy estimation. The resulting scheme consists of the minimization of a suitably-defined free energy over phase space subject to compatibility and equilibrium constraints. Distance-minimizing Data Driven schemes are recovered in the limit of zero temperature. We present selected numerical tests that establish the convergence properties of the max-ent Data Driven solvers and solutions. △ Less

Submitted 16 February, 2017; v1 submitted 6 February, 2017; originally announced February 2017.

Spectrum-splitting approach for Fermi-operator expansion in all-electron Kohn-Sham DFT calculations

Authors: Phani Motamarri, Vikram Gavini, Kaushik Bhattacharya, Michael Ortiz

Abstract: We present a spectrum-splitting approach to conduct all-electron Kohn-Sham density functional theory (DFT) calculations by employing Fermi-operator expansion of the Kohn-Sham Hamiltonian. The proposed approach splits the subspace containing the occupied eigenspace into a core-subspace, spanned by the core eigenfunctions, and its complement, the valence-subspace, and thereby enables an efficient co… ▽ More We present a spectrum-splitting approach to conduct all-electron Kohn-Sham density functional theory (DFT) calculations by employing Fermi-operator expansion of the Kohn-Sham Hamiltonian. The proposed approach splits the subspace containing the occupied eigenspace into a core-subspace, spanned by the core eigenfunctions, and its complement, the valence-subspace, and thereby enables an efficient computation of the Fermi-operator expansion by reducing the expansion to the valence-subspace projected Kohn-Sham Hamiltonian. The key ideas used in our approach are: (i) employ Chebyshev filtering to compute a subspace containing the occupied states followed by a localization procedure to generate non-orthogonal localized functions spanning the Chebyshev-filtered subspace; (ii) compute the Kohn-Sham Hamiltonian projected onto the valence-subspace; (iii) employ Fermi-operator expansion in terms of the valence-subspace projected Hamiltonian to compute the density matrix, electron-density and band energy. We demonstrate the accuracy and performance of the method on benchmark materials systems involving silicon nano-clusters up to 1330 electrons, a single gold atom and a six-atom gold nano-cluster. The benchmark studies on silicon nano-clusters revealed a staggering five-fold reduction in the Fermi-operator expansion polynomial degree by using the spectrum-splitting approach for accuracies in the ground-state energies of $\sim 10^{-4} Ha/atom$ with respect to reference calculations. Further, numerical investigations on gold suggest that spectrum splitting is indispensable to achieve meaningful accuracies, while employing Fermi-operator expansion. △ Less

Submitted 28 August, 2016; originally announced August 2016.

Comments: 12 pages, 5 figures

A multiscale model of distributed fracture and permeability in solids in all-round compression

Authors: Maria Laura De Bellis, Gabriele Della Vecchia, Michael Ortiz, Anna Pandolfi

Abstract: We present a microstructural model of permeability in fractured solids, where the fractures are described in terms of recursive families of parallel, equidistant cohesive faults. Faults originate upon the attainment of a tensile or shear resistance in the undamaged material. Secondary faults may form in a hierarchical orga- nization, creating a complex network of connected fractures that modify th… ▽ More We present a microstructural model of permeability in fractured solids, where the fractures are described in terms of recursive families of parallel, equidistant cohesive faults. Faults originate upon the attainment of a tensile or shear resistance in the undamaged material. Secondary faults may form in a hierarchical orga- nization, creating a complex network of connected fractures that modify the permeability of the solid. The undamaged solid may possess initial porosity and permeability. The particular geometry of the superposed micro-faults lends itself to an explicit analytical quantification of the porosity and permeability of the dam- aged material. The approach is particularly appealing as a means of modeling low permeability oil and gas reservoirs stimulated by hydraulic fracturing. △ Less

Submitted 17 October, 2016; v1 submitted 11 February, 2016; originally announced February 2016.

Nonequilibrium fluctuation-induced Casimir pressures in liquid mixtures

Authors: T. R. Kirkpatrick, J. M. Ortiz de Zárate, J. V. Sengers

Abstract: In this article we derive expressions for Casimir-like pressures induced by nonequilibrium concentration fluctuations in liquid mixtures. The results are then applied to liquid mixtures in which the concentration gradient results from a temperature gradient through the Soret effect. A comparison is made between the pressures induced by nonequilibrium concentration fluctuations in liquid mixtures a… ▽ More In this article we derive expressions for Casimir-like pressures induced by nonequilibrium concentration fluctuations in liquid mixtures. The results are then applied to liquid mixtures in which the concentration gradient results from a temperature gradient through the Soret effect. A comparison is made between the pressures induced by nonequilibrium concentration fluctuations in liquid mixtures and those induced by nonequilibrium temperature fluctuations in one-component fluids. Some suggestions for experimental verification procedures are also presented. △ Less

Submitted 14 January, 2016; originally announced January 2016.

Comments: 9 pages, 3 figures

Journal ref: Phys. Rev. E 93, 032117 (2016)

Oncotripsy: Targeting cancer cells selectively via resonant harmonic excitation

Authors: Stefanie Heyden, Michael Ortiz

Abstract: We investigate a method of selectively targeting cancer cells by means of ultrasound harmonic excitation at their resonance frequency, which we refer to as oncotripsy. The geometric model of the cells takes into account the cytoplasm, nucleus and nucleolus, as well as the plasma membrane and nuclear envelope. Material properties are varied within a pathophysiologically-relevant range. A first moda… ▽ More We investigate a method of selectively targeting cancer cells by means of ultrasound harmonic excitation at their resonance frequency, which we refer to as oncotripsy. The geometric model of the cells takes into account the cytoplasm, nucleus and nucleolus, as well as the plasma membrane and nuclear envelope. Material properties are varied within a pathophysiologically-relevant range. A first modal analysis reveals the existence of a spectral gap between the natural frequencies and, most importantly, resonant growth rates of healthy and cancerous cells. The results of the modal analysis are verified by simulating the fully-nonlinear transient response of healthy and cancerous cells at resonance. The fully nonlinear analysis confirms that cancerous cells can be selectively taken to lysis by the application of carefully tuned ultrasound harmonic excitation while simultaneously leaving healthy cells intact. △ Less

Submitted 20 July, 2020; v1 submitted 7 December, 2015; originally announced December 2015.

Comments: 22 pages

Physical origin of nonequilibrium fluctuation-induced forces in fluids

Authors: T. R. Kirkpatrick, J. M. Ortiz de Zárate, J. V. Sengers

Abstract: Long-range thermal fluctuations appear in fluids in nonequilibrium states leading to fluctuation-induced Casimir-like forces. Two distinct mechanisms have been identified for the origin of the long-range nonequilibrium fluctuations in fluids subjected to a temperature or concentration gradient. One is a coupling between the heat or mass-diffusion mode with a viscous mode in fluids subjected to a t… ▽ More Long-range thermal fluctuations appear in fluids in nonequilibrium states leading to fluctuation-induced Casimir-like forces. Two distinct mechanisms have been identified for the origin of the long-range nonequilibrium fluctuations in fluids subjected to a temperature or concentration gradient. One is a coupling between the heat or mass-diffusion mode with a viscous mode in fluids subjected to a temperature or concentration gradient. Another one is the spatial inhomogeneity of thermal noise in the presence of a gradient. We show that in fluids fluctuation-induced forces arising from mode coupling are several orders of magnitude larger than those from inhomogeneous noise. △ Less

Submitted 13 January, 2016; v1 submitted 12 November, 2015; originally announced November 2015.

Comments: 8 pages

Journal ref: Phys. Rev. E 93, 012148 (2016)

Data-driven computational mechanics

Authors: Trenton Kirchdoerfer, Michael Ortiz

Abstract: We develop a new computing paradigm, which we refer to as data-driven computing, according to which calculations are carried out directly from experimental material data and pertinent constraints and conservation laws, such as compatibility and equilibrium, thus bypassing the empirical material modeling step of conventional computing altogether. Data-driven solvers seek to assign to each material… ▽ More We develop a new computing paradigm, which we refer to as data-driven computing, according to which calculations are carried out directly from experimental material data and pertinent constraints and conservation laws, such as compatibility and equilibrium, thus bypassing the empirical material modeling step of conventional computing altogether. Data-driven solvers seek to assign to each material point the state from a prespecified data set that is closest to satisfying the conservation laws. Equivalently, data-driven solvers aim to find the state satisfying the conservation laws that is closest to the data set. The resulting data-driven problem thus consists of the minimization of a distance function to the data set in phase space subject to constraints introduced by the conservation laws. We motivate the data-driven paradigm and investigate the performance of data-driven solvers by means of two examples of application, namely, the static equilibrium of nonlinear three-dimensional trusses and linear elasticity. In these tests, the data-driven solvers exhibit good convergence properties both with respect to the number of data points and with regard to local data assignment. The variational structure of the data-driven problem also renders it amenable to analysis. We show that, as the data set approximates increasingly closely a classical material law in phase space, the data-driven solutions converge to the classical solution. We also illustrate the robustness of data-driven solvers with respect to spatial discretization. In particular, we show that the data-driven solutions of finite-element discretizations of linear elasticity converge jointly with respect to mesh size and approximation by the data set. △ Less

Submitted 14 October, 2015; originally announced October 2015.

Non-equilibrium concentration fluctuations in binary liquids with realistic boundary conditions

Authors: J. M. Ortiz de Zarate, T. R. Kirkpatrick, J. V. Sengers

Abstract: Because of the spatially long-ranged nature of spontaneous fluctuations in thermal non-equilibrium systems, they are affected by boundary conditions for the fluctuating hydrodynamic variables. In this paper we consider a liquid mixture between two rigid and impervious plates with a stationary concentration gradient resulting from a temperature gradient through the Soret effect. For liquid mixtures… ▽ More Because of the spatially long-ranged nature of spontaneous fluctuations in thermal non-equilibrium systems, they are affected by boundary conditions for the fluctuating hydrodynamic variables. In this paper we consider a liquid mixture between two rigid and impervious plates with a stationary concentration gradient resulting from a temperature gradient through the Soret effect. For liquid mixtures with a large Lewis number, we are able to obtain explicit analytical expressions for the intensity of the non-equilibrium concentration fluctuations as a function of the frequency $ω$ and the wave number $q$ of the fluctuations. In addition we elucidate the spatial dependence of the intensity of the non-equilibrium fluctuations responsible for a non-equilibrium Casimir effect. △ Less

Submitted 18 September, 2015; v1 submitted 6 May, 2015; originally announced May 2015.

Comments: 9 pages, 2 figures

Journal ref: European Physical Journal E, Vol. 38 (2015), article# 99

Non-equilibrium Casimir-like Forces in Liquid Mixtures

Authors: T. R. Kirkpatrick, J. M. Ortiz de Zárate, J. V. Sengers

Abstract: In this Letter we consider a liquid mixture confined between two thermally conducting walls subjected to a stationary temperature gradient. While in a one-component liquid non-equilibrium fluctuation forces appear inside the liquid layer only, non-equilibrium fluctuations in a mixture induce a Casimir-like force on the walls. The physical reason is that the temperature gradient induces large conce… ▽ More In this Letter we consider a liquid mixture confined between two thermally conducting walls subjected to a stationary temperature gradient. While in a one-component liquid non-equilibrium fluctuation forces appear inside the liquid layer only, non-equilibrium fluctuations in a mixture induce a Casimir-like force on the walls. The physical reason is that the temperature gradient induces large concentration fluctuations through the Soret effect. Unlike temperature fluctuations, non-equilibrium concentration fluctuations are also present near a perfectly thermally conducting wall. The magnitude of the fluctuation-induced Casimir force is proportional to the square of the Soret coefficient and is related to the concentration dependence of the heat and volume of mixing. △ Less

Submitted 19 June, 2015; v1 submitted 12 March, 2015; originally announced March 2015.

Comments: 5 pages, 1 figure

Journal ref: Physical Review Letters, 115 (2015) 035901

Spatial correlations in nonequilibrium reaction-diffusion problems by the Gillespie algorithm

Authors: Jorge Luis Hita, José María Ortiz de Zárate

Abstract: We present a study of the spatial correlation functions of a one-dimensional reaction-diffusion system in both equilibrium and out of equilibrium. For the numerical simulations we have employed the Gillespie algorithm dividing the system in cells to treat diffusion as a chemical process between adjacent cells. We find that the spatial correlations are spatially short ranged in equilibrium but beco… ▽ More We present a study of the spatial correlation functions of a one-dimensional reaction-diffusion system in both equilibrium and out of equilibrium. For the numerical simulations we have employed the Gillespie algorithm dividing the system in cells to treat diffusion as a chemical process between adjacent cells. We find that the spatial correlations are spatially short ranged in equilibrium but become long ranged in nonequilibrium. These results are in good agreement with theoretical predictions from fluctuating hydrodynamics for a one-dimensional system and periodic boundary conditions. △ Less

Submitted 14 May, 2013; originally announced May 2013.

Comments: 11 pages, 4 figures

Journal ref: Phys. Rev. E, 87 (2013), 052802

Coarse-graining Kohn-Sham Density Functional Theory

Authors: Phanish Suryanarayana, Kaushik Bhattacharya, Michael Ortiz

Abstract: We present a real-space formulation for coarse-graining Kohn-Sham Density Functional Theory that significantly speeds up the analysis of material defects without appreciable loss of accuracy. The approximation scheme consists of two steps. First, we develop a linear-scaling method that enables the direct evaluation of the electron density without the need to evaluate individual orbitals. We achiev… ▽ More We present a real-space formulation for coarse-graining Kohn-Sham Density Functional Theory that significantly speeds up the analysis of material defects without appreciable loss of accuracy. The approximation scheme consists of two steps. First, we develop a linear-scaling method that enables the direct evaluation of the electron density without the need to evaluate individual orbitals. We achieve this by performing Gauss quadrature over the spectrum of the linearized Hamiltonian operator appearing in each iteration of the self-consistent field method. Building on the linear-scaling method, we introduce a spatial approximation scheme resulting in a coarse-grained Density Functional Theory. The spatial approximation is adapted so as to furnish fine resolution where necessary and to coarsen elsewhere. This coarse-graining step enables the analysis of defects at a fraction of the original computational cost, without any significant loss of accuracy. Furthermore, we show that the coarse-grained solutions are convergent with respect to the spatial approximation. We illustrate the scope, versatility, efficiency and accuracy of the scheme by means of selected examples. △ Less

Submitted 8 August, 2012; originally announced August 2012.

Hydrodynamic Fluctuations in Laminar Fluid Flow. II. Fluctuating Squire Equation

Authors: Jose M. Ortiz de Zarate, Jan V. Sengers

Abstract: We use fluctuating hydrodynamics to evaluate the enhancement of thermally excited fluctuations in laminar fluid flow using plane Couette flow as a representative example. In a previous publication [J. Stat. Phys. vol. 144 (2011) p. 774] we derived the energy amplification arising from thermally excited wall-normal fluctuations by solving a fluctuating Orr-Sommerfeld equation. In the present paper… ▽ More We use fluctuating hydrodynamics to evaluate the enhancement of thermally excited fluctuations in laminar fluid flow using plane Couette flow as a representative example. In a previous publication [J. Stat. Phys. vol. 144 (2011) p. 774] we derived the energy amplification arising from thermally excited wall-normal fluctuations by solving a fluctuating Orr-Sommerfeld equation. In the present paper we derive the energy amplification arising from wall-normal vorticity fluctuation by solving a fluctuating Squire equation. The thermally excited wall-normal vorticity fluctuations turn out to yield the dominant contribution to the energy amplification. In addition, we show that thermally excited streaks, even in the absence of any externally imposed perturbations, are present in laminar fluid flow. △ Less

Submitted 19 June, 2012; originally announced June 2012.

Comments: 21 pages, 4 figures, Journal of Statistical Physics, accepted

Journal ref: Journal of Statistical Physics, vol. 150 (2013), pp. 540-558

Mesoscopic approach to granular crystal dynamics

Authors: Marcial Gonzalez, Jinkyu Yang, Chiara Daraio, Michael Ortiz

Abstract: We present a mesoscopic approach to granular crystal dynamics, which comprises a three-dimensional finite-element model and a one-dimensional regularized contact model. The approach investigates the role of vibrational-energy trapping effects in the dynamic behavior of one-dimensional chains of particles in contact (i.e., granular crystals), under small to moderate impact velocities. The only inpu… ▽ More We present a mesoscopic approach to granular crystal dynamics, which comprises a three-dimensional finite-element model and a one-dimensional regularized contact model. The approach investigates the role of vibrational-energy trapping effects in the dynamic behavior of one-dimensional chains of particles in contact (i.e., granular crystals), under small to moderate impact velocities. The only inputs of the models are the geometry and the elastic material properties of the individual particles that form the system. We present detailed verification results and validate the model comparing its predictions with experimental data. This approach provides a physically sound, first-principle description of dissipative losses in granular systems. △ Less

Submitted 16 November, 2011; originally announced November 2011.

Journal ref: Phys. Rev. E 85, 016604, 2012

The Pierre Auger Observatory V: Enhancements

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, A. Aminaei, L. Anchordoqui, S. Andringa, T. Antičić, A. Anzalone, C. Aramo, E. Arganda, F. Arqueros, H. Asorey, P. Assis, J. Aublin, M. Ave, M. Avenier , et al. (471 additional authors not shown)

Abstract: Ongoing and planned enhancements of the Pierre Auger Observatory Ongoing and planned enhancements of the Pierre Auger Observatory △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: Contributions to 32nd International Cosmic Ray Conference, Beijing, China, August 2011

The Pierre Auger Observatory IV: Operation and Monitoring

Authors: The Pierre Auger Collaboration, P. Abreu, M. Aglietta, E. J. Ahn, I. F. M. Albuquerque, D. Allard, I. Allekotte, J. Allen, P. Allison, J. Alvarez Castillo, J. Alvarez-Muñiz, M. Ambrosio, A. Aminaei, L. Anchordoqui, S. Andringa, T. Antičić, A. Anzalone, C. Aramo, E. Arganda, F. Arqueros, H. Asorey, P. Assis, J. Aublin, M. Ave, M. Avenier , et al. (471 additional authors not shown)

Abstract: Technical reports on operations and monitoring of the Pierre Auger Observatory Technical reports on operations and monitoring of the Pierre Auger Observatory △ Less

Submitted 24 July, 2011; originally announced July 2011.

Comments: Constributions to 32nd International Cosmic Ray Conference, Beijing, China, August 2011

Optimal Uncertainty Quantification

Authors: Houman Owhadi, Clint Scovel, Timothy John Sullivan, Mike McKerns, Michael Ortiz

Abstract: We propose a rigorous framework for Uncertainty Quantification (UQ) in which the UQ objectives and the assumptions/information set are brought to the forefront. This framework, which we call \emph{Optimal Uncertainty Quantification} (OUQ), is based on the observation that, given a set of assumptions and information about the problem, there exist optimal bounds on uncertainties: these are obtained… ▽ More We propose a rigorous framework for Uncertainty Quantification (UQ) in which the UQ objectives and the assumptions/information set are brought to the forefront. This framework, which we call \emph{Optimal Uncertainty Quantification} (OUQ), is based on the observation that, given a set of assumptions and information about the problem, there exist optimal bounds on uncertainties: these are obtained as values of well-defined optimization problems corresponding to extremizing probabilities of failure, or of deviations, subject to the constraints imposed by the scenarios compatible with the assumptions and information. In particular, this framework does not implicitly impose inappropriate assumptions, nor does it repudiate relevant information. Although OUQ optimization problems are extremely large, we show that under general conditions they have finite-dimensional reductions. As an application, we develop \emph{Optimal Concentration Inequalities} (OCI) of Hoeffding and McDiarmid type. Surprisingly, these results show that uncertainties in input parameters, which propagate to output uncertainties in the classical sensitivity analysis paradigm, may fail to do so if the transfer functions (or probability distributions) are imperfectly known. We show how, for hierarchical structures, this phenomenon may lead to the non-propagation of uncertainties or information across scales. In addition, a general algorithmic framework is developed for OUQ and is tested on the Caltech surrogate model for hypervelocity impact and on the seismic safety assessment of truss structures, suggesting the feasibility of the framework for important complex systems. The introduction of this paper provides both an overview of the paper and a self-contained mini-tutorial about basic concepts and issues of UQ. △ Less

Submitted 23 May, 2012; v1 submitted 2 September, 2010; originally announced September 2010.

Comments: 90 pages. Accepted for publication in SIAM Review (Expository Research Papers). See SIAM Review for higher quality figures

Journal ref: SIAM Rev. 55(2):271--345, 2013

Upgrade of the Control System of the IFUNAM's Pelletron Accelerator

Authors: R. Macias, E. Chavez, M. E. Ortiz, K. Lopez, A. Huerta, M. C. Verde

Abstract: In 1995 a 9SDH-2 Pelletron from NEC was installed at IFUNAM (Instituto de Fisica, Universidad Nacional Autonoma de Mexico). Two beam lines have been operational since then and two new lines have been built. In order to perform the planned projects in this grown facility, an upgrading of the original manual control system is required. The proposed new control system takes advantage of the existin… ▽ More In 1995 a 9SDH-2 Pelletron from NEC was installed at IFUNAM (Instituto de Fisica, Universidad Nacional Autonoma de Mexico). Two beam lines have been operational since then and two new lines have been built. In order to perform the planned projects in this grown facility, an upgrading of the original manual control system is required. The proposed new control system takes advantage of the existing devices and incorporates the electronics needed for the newer beam lines. The control software from NEC, has been modified to accommodate the larger requirements. It runs on the same dedicated computer but receives commands from a new installed host. Both computers communicate through a local network sharing the accelerator database. The new host computer also handles all parameters related to the new lines. In the future, the old computer will be replaced in order to expand the possibilities of the system and use a friendlier graphical interface. In this work we present the changes made to the control software, the new acquisition and control devices installed and the integration of the old and new systems. Finally, the beam control procedure based on the analysis of the beam profile that is under development is also shown. △ Less

Submitted 8 November, 2001; originally announced November 2001.

Comments: 3 pages, 6 figures, Paper for ICALEPCS 2001 (Poster), WEAP032, PDF

Journal ref: eConf C011127 (2001) WEAP032

Light Scattering from Nonequilibrium Concentration Fluctuations in a Polymer solution

Authors: W. B. Li, K. J. Zhang, J. V. Sengers, R. W. Gammon, J. M. Ortiz de Zarate

Abstract: We have performed light-scattering measurements in dilute and semidilute polymer solutions of polystyrene in toluene when subjected to stationary temperature gradients. Five solutions with concentrations below and one solution with a concentration above the overlap concentration were investigated. The experiments confirm the presence of long-range nonequilibrium concentration fluctuations which… ▽ More We have performed light-scattering measurements in dilute and semidilute polymer solutions of polystyrene in toluene when subjected to stationary temperature gradients. Five solutions with concentrations below and one solution with a concentration above the overlap concentration were investigated. The experiments confirm the presence of long-range nonequilibrium concentration fluctuations which are proportional to $(\nabla T)^2/k^4$, where $\nabla T$ is the applied temperature gradient and $k$ is the wave number of the fluctuations. In addition, we demonstrate that the strength of the nonequilibrium concentration fluctuations, observed in the dilute and semidilute solution regime, agrees with theoretical values calculated from fluctuating hydrodynamics. Further theoretical and experimental work will be needed to understand nonequilibrium fluctuations in polymer solutions at higher concentrations. △ Less

Submitted 24 February, 2000; originally announced February 2000.

Comments: revtex, 16 pages, 7 figures. J. Chem. Phys., to appear

Journal ref: J. Chem. Phys. 112, 9139 (2000)