We present an Atomic Cluster Expansion (ACE) machine learned potential developed for high-fidelity atomistic simulations of hydrocarbons, targeting pressures and temperatures near and above supercritical fluid regimes for molecular fluids. A diverse set of stoichiometries were covered in training, including 1:0 (pure carbon), 1:4 (methane), and 1:1 (benzene), and rich bonding environments sampled at supercritical temperatures, hydrogen rich, reactive mixtures where metastable stoichiometries arise, including 1:2 (ethylene) and 1:3 (ethane). A high-fidelity training database was constructed by performing large-scale quantum molecular dynamic simulations [density functional theory (DFT) MD] of diamond, graphite, methane, and benzene.
View Article and Find Full Text PDFThis study conducted integrated experiments and computational modeling to investigate the speeds of a developing shock within granular salt and analyzed the effect of various impact velocities up to 245 m/s. Experiments were conducted on table salt utilizing a novel setup with a considerable bore length for the sample, enabling visualization of a moving shock wave. Experimental analysis using particle image velocimetry enabled the characterization of shock velocity and particle velocity histories.
View Article and Find Full Text PDFThere is a general consensus that the doctor-patient interview should be as productive and efficient as possible. This is becoming increasingly difficult in a health care insurance system that demands shorter appointment times. Clinicians must therefore find ways to condense the clinical encounter without sacrificing quality.
View Article and Find Full Text PDFThe high cost of density functional theory (DFT) has hitherto limited the ab initio prediction of the equation of state (EOS). In this article, we employ a combination of large scale computing, advanced simulation techniques, and smart data science strategies to provide an unprecedented ab initio performance analysis of the high explosive pentaerythritol tetranitrate (PETN). Comparison to both experiment and thermochemical predictions reveals important quantitative limitations of DFT for EOS prediction and thus the assessment of high explosives.
View Article and Find Full Text PDFComplex plasma mixtures with three or more components are often encountered in astrophysics or in inertial confinement fusion (ICF) experiments. For mixtures containing species with large differences in atomic number Z, the modeling needs to consider at the same time the kinetic theory for low-Z elements combined with the theory of strongly coupled plasma for high-Z elements, as well as all the intermediate situations that can appear in multicomponent systems. For such cases, we study the pair distribution functions, self-diffusions, mutual diffusion, and viscosity for ternary mixtures at extreme conditions.
View Article and Find Full Text PDFWe present the formulation, simulations, and results for multicomponent mutual diffusion coefficients in the warm, dense matter regime. While binary mixtures have received considerable attention for mass transport, far fewer studies have addressed ternary and more complex systems. We therefore explicitly examine ternary systems utilizing the Maxwell-Stefan formulation that relates diffusion to gradients in the chemical potential.
View Article and Find Full Text PDFWe performed a series of plate impact experiments on NH gas initially at room temperature and at a pressure of ∼100 psi. Shocked states were determined by optical velocimetry and the temperatures by optical pyrometry, yielding compression ratios of ∼5-10 and second shock temperatures in excess of 7500 K. A first-principles statistical mechanical (thermochemical) approach that included chemical dissociation yielded reasonable agreement with experimental results on the principal Hugoniot, even with interparticle interactions neglected.
View Article and Find Full Text PDFTransport properties of mixtures of elements in the dense plasma regime play an important role in natural astrophysical and experimental systems, e.g., inertial confinement fusion.
View Article and Find Full Text PDFWe study how concentration changes ionic transport properties along isobars-isotherms for a mixture of hydrogen and silver, representative of turbulent layers relevant to inertial confinement fusion and astrophysics. Hydrogen will typically be fully ionized while silver will be only partially ionized but can have a large effective charge. This will lead to very different physical conditions for the H and Ag.
View Article and Find Full Text PDFWe present an alternative method for interpreting the velocity autocorrelation function (VACF) of a fluid with application to extracting the entropy in a manner similar to the methods developed by Lin et al. [J. Chem.
View Article and Find Full Text PDFOrbital-free molecular dynamics simulations are used to benchmark two popular models for hot dense plasmas: the one component plasma (OCP) and the Yukawa model. A unified concept emerges where an effective OCP (EOCP) is constructed from the short-range structure of the plasma. An unambiguous ionization and the screening length can be defined and used for a Yukawa system, which reproduces the long-range structure with finite compressibility.
View Article and Find Full Text PDFWe have performed finite-temperature density functional theory molecular dynamics simulations on dense methane, ammonia, and water mixtures (CH4:NH3:H2O) for various compositions and temperatures (2000 K ≤ T ≤ 10,000 K) that span a set of possible conditions in the interiors of ice-giant exoplanets. The equation-of-state, pair distribution functions, and bond autocorrelation functions (BACF) were used to probe the structure and dynamics of these complex fluids. In particular, an improvement to the choice of the cutoff in the BACF was developed that allowed analysis refinements for density and temperature effects.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
August 2015
We present simulations of a four-component mixture of HCNO with orbital free molecular dynamics (OFMD). These simulations were conducted for 5-200 eV with densities ranging between 0.184 and 36.
View Article and Find Full Text PDFPurpose: The paradox of primary care is the observation that primary care is associated with apparently low levels of evidence-based care for individual diseases, but systems based on primary care have healthier populations, use fewer resources, and have less health inequality. The purpose of this article is to explore, from a complex systems perspective, mechanisms that might account for the effects of primary care beyond disease-specific care.
Methods: In an 8-session, participatory group model-building process, patient, caregiver, and primary care clinician community stakeholders worked with academic investigators to develop and refine an agent-based computer simulation model to test hypotheses about mechanisms by which features of primary care could affect health and health equity.
Phys Rev E Stat Nonlin Soft Matter Phys
January 2015
A recent and unexpected discrepancy between ab initio simulations and the interpretation of a laser shock experiment on aluminum, probed by x-ray Thomson scattering (XRTS), is addressed. The ion-ion structure factor deduced from the XRTS elastic peak (ion feature) is only compatible with a strongly coupled out-of-equilibrium state. Orbital free molecular dynamics simulations with ions colder than the electrons are employed to interpret the experiment.
View Article and Find Full Text PDFWe present a theory for the number fluctuations of a quasi-two-dimensional (quasi-2D) dipolar Bose-Einstein condensate measured with finite resolution cells. We show that when the dipoles are tilted to have a component parallel to the plane of the trap, the number fluctuations become anisotropic, i.e.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
October 2014
In the warm dense matter (WDM) regime, material properties like diffusion and viscosity can be obtained from lengthy quantum molecular dynamics simulations, where the quantum behavior of the electrons is represented using either Kohn-Sham or orbital-free density functional theory. To reduce the simulation duration, we fit the time dependence of the autocorrelation functions (ACFs) and then use the fit to find values of the diffusion and viscosity. This fitting procedure avoids noise in the long time behavior of the ACFs.
View Article and Find Full Text PDFWe have performed nonequilibrium classical and quantum-mechanical molecular dynamics simulations that follow the interpenetration of deuterium-tritium (DT) and carbon (C) components through an interface initially in hydrostatic and thermal equilibrium. We concentrate on the warm, dense matter regime with initial densities of 2.5-5.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
December 2013
Extending the well-known Thomas-Fermi Z-scaling laws to the Coulomb coupling parameter, we investigate the stabilization of the ionic coupling in isochoric heating [Clérouin et al., Phys. Rev.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
February 2013
We have performed a systematic study of lithium hydride (LiH), using orbital-free molecular dynamics, with a focus on mass transport properties such as diffusion and viscosity by extending our previous studies at the lower end of the warm, dense matter regime to cover a span of densities from ambient to 10-fold compressed and temperatures from 10 eV to 10 keV. We determine analytic formulas for self- and mutual-diffusion coefficients, and viscosity, which are in excellent agreement with our molecular dynamics results, and interpolate smoothly between liquid and dense plasma regimes. In addition, we find the orbital-free calculations begin to agree with the Brinzinskii-Landau formula above about 250 eV at which point the medium becomes fully ionized.
View Article and Find Full Text PDFWe show that the correlation dynamics in coherently excited doubly excited resonances of helium can be followed in real time by two-photon interferometry. This approach promises to map the evolution of the two-electron wave packet onto experimentally easily accessible noncoincident single-electron spectra. We analyze the interferometric signal in terms of a semianalytical model which is validated by a numerical solution of the time-dependent two-electron Schrödinger equation in its full dimensionality.
View Article and Find Full Text PDFWe study the superfluid character of a dipolar Bose-Einstein condensate (DBEC) in a quasi-two dimensional geometry. We consider the dipole polarization to have some nonzero projection into the plane of the condensate so that the effective interaction is anisotropic in this plane, yielding an anisotropic dispersion relation. By performing direct numerical simulations of a probe moving through the DBEC, we observe the sudden onset of drag or creation of vortex-antivortex pairs at critical velocities that depend strongly on the direction of the probe's motion.
View Article and Find Full Text PDFWe use universality in two-body dipolar physics to study three-body recombination. We present results for the universal structure of weakly bound two-dipole states that depend only on the s-wave scattering length (a). We study threshold three-body recombination rates into weakly bound dimer states as a function of the scattering length.
View Article and Find Full Text PDFPhys Rev E Stat Nonlin Soft Matter Phys
July 2009
We describe a method for evolving the projected Gross-Pitaevskii equation (PGPE) for an interacting Bose gas in a harmonic-oscillator potential, with the inclusion of a long-range dipolar interaction. The central difficulty in solving this equation is the requirement that the field is restricted to a small set of prescribed modes that constitute the low-energy c -field region of the system. We present a scheme, using a Hermite-polynomial-based spectral representation, which precisely implements this mode restriction and allows an efficient and accurate solution of the dipolar PGPE.
View Article and Find Full Text PDFWe explore the impact of the short-range interaction on the scattering of ground state polar molecules and study the transition from a weak to strong dipolar scattering over an experimentally reasonable range of energies and electric field values. In the strong dipolar limit, the scattering scales with respect to a dimensionless quantity defined by mass, induced dipole moment, and collision energy. The scaling has implications for all quantum mechanical dipolar scattering.
View Article and Find Full Text PDF