High-pressure phases of liquid silane studied using first-principles molecular dynamics simulations
Not Available
Badjo, Y. H. (ORCID:0009000739381225)↗
SEARCH · Engineering Papers
Results for “liquid dynamics”
Search indexed NASA NTRS and DOE OSTI research on propulsion, heat transfer, battery materials and energy systems. Follow report and document links to the original sources.
Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.
Phonon dynamics in the site-disordered Kitaev spin liquid
The Kitaev honeycomb model provides a paradigmatic example of an exactly solvable quantum spin liquid (QSL), where spins fractionalize into itinerant Majorana fermions coupled to a static background of ℤ 2 gauge fluxes. This model has attracted significant interest due to its potential experimental realization in spin-orbit Mott insulators like 𝛼−RuCl 3 . Among various experimental techniques, ultrasound measurements of sound attenuation have emerged as a promising approach to detect spin fractionalization in these materials. However, deviations from the ideal Kitaev model, often due to disorder, introduce localized modes that dominate the low-energy physics. To investigate these effects, we calculate the sound attenuation coefficient in the site-disordered Kitaev honeycomb model under an applied magnetic field that breaks time-reversal symmetry. Our analysis reveals that quasilocalized modes from quasivacancies impact the phonon self-energy, yet the sound attenuation coefficient maintains its sixfold symmetry and linear temperature dependence at low temperatures, even with disorder. This robustness underscores sound attenuation's reliability as a probe for spin fractionalization. As a result, we show that this linear behavior persists under an external field and in the random-flux sector, highlighting the persistent influence of fractionalized excitations, despite disorder.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Trends in Ionic Liquid Diffusion in Graphite Slit Pores from Molecular Dynamics Simulation and Machine Learning
Explore the source record for details and available documents.
Greathouse, Jeffery A.↗
Dynamical Anyon Generation in Kitaev Honeycomb Non-Abelian Spin Liquids
Relativistic Mott insulators known as “Kitaev materials” potentially realize spin liquids hosting non-Abelian anyons. Motivated by fault-tolerant quantum-computing applications in this setting, we introduce a dynamical anyon-generation protocol that exploits universal edge physics. The setup features holes in the spin liquid, which define energetically cheap locations for non-Abelian anyons, connected by a narrow bridge that can be tuned between spin liquid and topologically trivial phases. We show that modulating the bridge from trivial to spin liquid over intermediate time scales—quantified by analytics and extensive simulations—deposits non-Abelian anyons into the holes with O(1) probability. The required bridge manipulations can be implemented by integrating the Kitaev material into magnetic tunnel junction arrays that engender locally tunable exchange fields. Combined with existing readout strategies, our protocol reveals a path to topological qubit experiments in Kitaev materials at zero applied magnetic field.
75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
Solid Electrolyte Interphase Formation at the Ionic Liquid Electrolyte—Lithium-Metal Interface Using an Ab Initio Molecular Dynamics Approach
An interfacial study is performed using ab initio molecular dynamics (AIMD) simulations to elucidate the electrochemical interfacial stability of an ionic liquid (IL) as component of an electrolyte in contact with a Li-metal anode when an additive, vinylene carbonate (VC), is added to the electrolyte. The IL electrolyte is composed of 1-butyl-1-methylpyrrolidinium bis(fluoro-sulfonyl)-imide (PYR14 + FSI − ) with lithium bis(trifluoromethane-sulfonyl)imide salt (Li + TFSI − ) in a 8:2 concentration ratio. Reactions observed at the interface electrolyte/Li-metal during the AIMD simulations are verified by electronic structure calculations using density functional theory calculations considering a continuum solvent model. We found that the addition of VC to the IL prevents further reduction of the counter-ion TFSI − and retards reduction of FSI − in comparison to the case when the additive VC is absent. The PYR14 + reduction is never observed with or without the additive. The 1st, 2nd, and 3rd electron affinities of the PYR14 + are −0.08, −0.36 and −3.65 eV, indicating little to large repulsive behavior against electrons. The opposite is observed for the two negative ions, FSI − and TFSI − , which yield 1st electron affinities of 1.73 and 1.83 eV, respectively. Reactions obtained under open circuit conditions. Dissociation mechanisms are presented for all IL components reduced at the IL/Li-metal interface.
Electrochemistry↗
Ultrafast vibrational dynamics of a trigonal planar anionic probe in ionic liquids (ILs): A two-dimensional infrared (2DIR) spectroscopic investigation
A major impediment limiting the widespread application of ionic liquids (ILs) is their high shear viscosity. Incorporation of a tricyanomethanide (TCM - ) anion in ILs leads to low shear viscosity and improvement of several characteristics suitable for large scale applications. However, properties including interactions of TCM - with the local environment and dynamics of TCM - have not been thoroughly investigated. Herein, we have studied the ultrafast dynamics of TCM - in several imidazolium ILs using linear IR and two-dimensional infrared spectroscopy techniques. The spectral diffusion dynamics of the CN stretching modes of TCM - in all ILs exhibit a nonexponential behavior with a short time component of ~2 ps and a long time component spanning ~9 ps to 14 ps. The TCM - vibrational probe reports a significantly faster relaxation of ILs compared to those observed previously using linear vibrational probes, such as thiocyanate and selenocyanate. Our results indicate a rapid relaxation of the local ion-cage structure embedding the vibrational probe in the ILs. The faster relaxation suggests that the lifetime of the local ion-cage structure decreases in the presence of TCM - in the ILs. Linear IR spectroscopic results show that the hydrogen-bonding interaction between TCM - and imidazolium cations in ILs is much weaker. Shorter ion-cage lifetimes together with weaker hydrogen-bonding interactions account for the low shear viscosity of TCM - based ILs compared to commonly used ILs. In addition, this work demonstrates that TCM - can be used as a potential vibrational reporter to study the structure and dynamics of ILs and other molecular systems.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Atomistic simulations of the thermal conductivity of liquids
We present a method based on sinusoidal approach to equilibrium molecular dynamics (SAEMD) to compute the thermal conductivity of liquids Similar to nonequilibrium molecular dynamics, and unlike equilibrium simulations based on the Green-Kubo formalism, the method only requires the calculation of forces and total energies. The evaluation of heat fluxes and energy densities is not necessary, thus offering the promise of efficiently implementing first principles simulations based on density functional theory or deep molecular dynamics. Our approach is a generalization of SAEMD for solids, where the thermal conductivity is computed in the steady state, instead of a transient regime, thus properly taking into account diffusive terms in the heat equation. We present results for liquid water at ambient conditions and under pressure and discuss simulation requirements to obtain converged values of the thermal conductivity as a function of size and simulation time.
36 MATERIALS SCIENCE↗
Conductivity and Transference Numbers in Lithium Salt-Doped Block Copolymeric Ionic Liquid Electrolytes
We present multiscale molecular dynamics simulation results comparing the conductivity and transference numbers in lithium salt-doped polymeric ionic liquids (PolyILs) and the lamellar phases of block copolymeric ionic liquids (coPolyILs). In both systems, the anion mobilities decreased with salt loading. Lithium ions exhibited negative mobilities in both systems, but the magnitudes decreased with an increase in salt concentrations. More interestingly, the anion mobilities were lower in the lamellar systems compared to homopolymers in magnitude, but the lithium ion mobilities and transference numbers were less negative in such systems. We examine the anion–cation and lithium–anion interactions in terms of radial distribution functions, coordination characteristics, and ion-pair relaxation timescales. Based on such analyses, we rationalize the salt concentration dependencies as a result of the interfacial interactions in lamellar systems and the competition between anion–cation and lithium–anion interactions in both PolyILs and coPolyILs. Altogether, the findings presented in this study demonstrate that the modified anion–cation and lithium–anion interactions in the microphase-separated coPolyILs may provide a strategy for realizing higher lithium ion transference numbers relative to the homopolymeric counterparts.
36 MATERIALS SCIENCE↗
Polar liquids at charged interfaces: A dipolar shell theory
The structure of polar liquids and electrolytic solutions, such as water and aqueous electrolytes, at interfaces underlies numerous phenomena in physics, chemistry, biology, and engineering. In this work, we develop a continuum theory that captures the essential features of dielectric screening by polar liquids at charged interfaces, including decaying spatial oscillations in charge and mass, starting from the molecular properties of the solvent. The theory predicts an anisotropic dielectric tensor of interfacial polar liquids previously studied in molecular dynamics simulations. We explore the effect of the interfacial polar liquid properties on the capacitance of the electrode/electrolyte interface and on hydration forces between two plane-parallel polarized surfaces. In the linear response approximation, we obtain simple formulas for the characteristic decay lengths of molecular and ionic profiles at the interface.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Intermolecular correlations of liquid and glassy CS 2 studied by synchrotron radiation x-ray diffraction
How is the orientation of molecular liquids ordered on cooling? What are the basic structures of molecular glasses, e.g., close to the crystalline structure or some special structures such as icosahedral cluster? These are long-standing questions in liquid and glass physics. Here, we have constructed a novel cryostat to prepare simple molecular glasses by vapor deposition and performed in situ synchrotron radiation x-ray diffraction experiments. The glassy state of a simple molecule CS 2 , which cannot be vitrified by normal liquid quenching, was successfully prepared with this instrument, and its diffraction data were collected in a wide Q-range of 0.16–25.7 Å -1 with a high-energy diffractometer at BL04B2, SPring-8. The diffraction data of liquid CS 2 were also recorded in a wide temperature range of 160–300 K. These diffraction data were analyzed with molecular dynamics simulations and reverse Monte Carlo modelings to investigate orientational correlation. From the obtained 3D structure models, the orientational correlation between neighboring CS 2 molecules was investigated quantitatively as a function of temperature. At room temperature, the parallel and T-shaped arrangements are preferred for the nearest neighbor correlation. On cooling, these arrangements are developed gradually, and its rate became prominent below the melting temperature (162 K). In the glassy state, the slipped-parallel arrangement is dominant as well as the T-shaped arrangement. Both arrangements appear in the CS 2 crystal, indicating that the structure of glassy CS 2 is close to that of crystalline CS 2 .
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Highly tail-asymmetric lipids interdigitate and cause bidirectional ordering
Phospholipids form structurally and compositionally diverse membranes. A less studied type of compositional diversity involves phospholipid tail variety. Some phospholipids contain two acyl tails which differ in length. These tail-asymmetric lipids are shown to contribute to temperature sensitivity, oxygen adaptability, and membrane fluidity. Membranes of a highly virulent intracellular bacterium, Francisella tularensis, contain highly tail-asymmetric 1-lignoceroyl-2-decanoyl-sn-glycero-3-phosphatidylethanolamine (XJPE) lipids which were previously shown to inhibit inflammatory responses in host cells. XJPE tails have unusually high asymmetry, and how they contribute to membrane properties on a molecular level is unknown. Here, we use small angle X-ray scattering and molecular dynamics simulations to investigate how varying XJPE ratios alters properties of simple membranes. Our results demonstrate that at high concentration they promote liquid-to-gel transition in otherwise liquid membranes, while at low concentration they are tolerated well, minimally altering membrane properties. In liquid membranes, XJPE lipids dynamically adopt two main conformations; with the long tail extended into the opposing leaflet or bent-back residing in its own leaflet. When added to both leaflets XJPE primarily adopts an extended confirmation, while asymmetric addition results in more bent-back orientations. The former increases tail ordering and the latter decreases it. XJPE tails adopt different conformations that induce composition- and leaflet-dependent bidirectional effect on membrane fluidity and this suggests that Francisella tularensis could use tail asymmetry to facilitate vesicle fusion and destabilize host cells. The effect of tail-asymmetric lipids on complex membranes should be further investigated to reveal the regulatory roles of high tail asymmetry.
59 BASIC BIOLOGICAL SCIENCES↗
Tuning the melting point of selected ionic liquids through adjustment of the cation’s dipole moment.
In prior research with thermally robust salts [Cassity et al., PCCP, 2017, 19, 31560] it was noted that an increase in the dipole moment of the cation generally led to a decrease in the melting point. Molecular dynamics simulations of the liquid state revealed that an increased dipole moment reduces cation-cation repulsions through dipole-dipole alignment. This was believed to reduce the liquid phase enthalpy, which would tend to lower the melting point of the IL. Here, we also test this principle by replacing hydrogen atoms with fluorine atoms at selected positions within the cation. This allows us to alter the electrostatics of the cation without substantially affecting the sterics. Furthermore, the strength of the dipole moment can be controlled by choosing different positions within the cation for replacement. We studied variants of four different parent cations paired with bistriflimide and determined their melting points, and enthalpies and entropies of fusion through DSC experiments. The decreases in the melting point were determined to be enthalpically driven. We found that the dipole moment of the cation, as determined by quantum chemical calculations, is inversely correlated with the melting point of the given compound. Molecular dynamics simulations of the crystalline and solid states of two isomers showed differences in their enthalpies of fusion that closely matched those seen experimentally. Moreover, this reduction in the enthalpy of fusion was determined to be caused by an increase in the enthalpy of the crystalline state. We provide evidence that dipole-dipole interactions between cations leads to the formation of cationic domains in the crystalline state. These cationic associations partially block favourable cation–anion interactions, which are recovered upon melting. If, however, the dipole- dipole interactions between cations is too strong they have a tendency to form glasses. This study provides a design rule for lowering the melting point of structurally similar ILs by altering their dipole moment.
32 ENERGY CONSERVATION, CONSUMPTION, AND UTILIZATI↗
Atomic dynamics in fluids: Normal mode analysis revisited
Developing microscopic understanding of the thermal properties of liquids is challenging due to their strong dynamic disorder, which prevents characterization of the atomic degrees of freedom. There have been significant research interests in the past few decades to extend the normal mode analysis for solids to instantaneous structures of liquids. However, the nature of normal modes that arise from these unstable structures is still elusive. Here, in this paper, we explore the instantaneous eigenmodes of dynamical matrices of various Lennard-Jones argon liquid and gas systems at high temperatures and show that the normal modes can be interpreted as an interpolation of T →∞ (gas) and T=0 (solid) mode descriptions. We find that normal modes become increasingly collisional and translational, recovering atomistic gaslike behavior rather than vibrational with increase in temperature, suggesting that normal modes in liquids may be described by both solidlike and gaslike modes.
74 ATOMIC AND MOLECULAR PHYSICS↗
Activated relaxation in supercooled monodisperse atomic and polymeric WCA fluids: Simulation and ECNLE theory
Here, we combine simulation and Elastically Collective Nonlinear Langevin Equation (ECNLE) theory to study the activated relaxation in monodisperse atomic and polymeric Weeks–Chandler–Andersen (WCA) liquids over a wide range of temperatures and densities in the supercooled regime under isochoric conditions. By employing novel crystal-avoiding simulations, metastable equilibrium dynamics is probed in the absence of complications associated with size polydispersity. Based on a highly accurate structural input from integral equation theory, ECNLE theory is found to describe well the simulated density and temperature dependences of the alpha relaxation time of atomic fluids using a single system-specific parameter, a c , that reflects the nonuniversal relative importance of local cage and collective elastic barriers. For polymer fluids, the explicit dynamical effect of local chain connectivity is modeled at the fundamental dynamic free energy trajectory level based on a different parameter, N c , that quantifies the degree of intramolecular correlation of bonded segment activated barrier hopping. For the flexible chain model studied, a physically intuitive value of N c ≈ 2 results in good agreement between simulation and theory. A direct comparison between atomic and polymeric systems reveals that chain connectivity can speed up activated segmental relaxation due to weakening of equilibrium packing correlations but can slow down relaxation due to local bonding constraints. The empirical thermodynamic scaling idea for the alpha time is found to work well at high densities or temperatures but fails when both density and temperature are low. The rich and subtle behaviors revealed from simulation for atomic and polymeric WCA fluids are all well captured by ECNLE theory.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Fluid Flow Observations of The Spray Near-Field Using High-Speed X-Ray Imaging
Sprays are important to may industrial processes, but the near-field region where spray formation takes place is difficult to visualize and characterize through visible light imaging because this region is typically optically dense. X-ray imaging is one technique that can penetrate this optically dense region and provide unique observations of the spray formation process. This study uses high-speed white beam X-ray imaging conducted at the Advanced Photon Source to capture spray formation dynamics from a canonical airblast atomizer. In this work, a range of momentum flux ratios are visualized by varying the gas flow rate while the liquid flow rate is held constant. The use of X-rays show overlapping flow structures including bags, ligaments, mushrooms, and webs. Air bubbles are captured within the liquid regions, and a very dynamic crown forms at the liquid needle exit at higher momentum flux ratios. Finally, at the highest momentum flux ratios, air bubbles and liquid droplets are observed migrating upstream in the center of the crown, and this motion is attributed to a toroidal recirculation cell at the nozzle exit.
47 OTHER INSTRUMENTATION↗
Computing the viscous effect in early-time drop impact dynamics
The impact of a liquid drop on a solid surface involves many intertwined physical effects, and is influenced by drop velocity, surface tension, ambient pressure and liquid viscosity, among others. Experiments by Kolinski et al. ( Phys. Rev. Lett. , vol. 112, no. 13, 2014 b , p. 134501) show that the liquid–air interface begins to deviate away from the solid surface even before contact. They found that the lift-off of the interface starts at a critical time that scales with the square root of the kinematic viscosity of the liquid. To understand this, we study the approach of a liquid drop towards a solid surface in the presence of an intervening gas layer. We take a numerical approach to solve the Navier–Stokes equations for the liquid, coupled to the compressible lubrication equations for the gas, in two dimensions. With this approach, we recover the experimentally captured early time effect of liquid viscosity on the drop impact, but our results show that lift-off time and liquid kinematic viscosity have a more complex dependence than the square-root scaling relationship. We also predict the effect of interfacial tension at the liquid–gas interface on the drop impact, showing that it mediates the lift-off behaviour.
42 ENGINEERING↗
Thermodynamics–Structure–Dynamics Correlations and Nonuniversal Effects in the Elastically Collective Activated Hopping Theory of Glass-Forming Liquids
We employ the microscopic Elastically Collective Nonlinear Langevin Equation (ECNLE) theory of activated dynamics in combination with crystal-avoiding simulations to study four inter-related questions for metastable monodisperse hard sphere fluids. The first is how significantly improved integral equation theory structural input (Modified-Verlet (MV) closure) changes the dynamical predictions of ECNLE theory. The main consequence is a modest enhancement of the importance of the collective elastic barrier relative to its local cage contribution, which increases the alpha relaxation time and fragility relative to prior results based on the Percus–Yevick closure. Second, ECNLE-MV theory predictions for the alpha time and self-diffusion constant in the metastable regime are quantitatively compared to our new simulations. The small adjustment of a numerical prefactor that enters the collective elastic barrier leads to quantitative agreement over three decades. Third, using the more accurate MV structural input, ECNLE theory is shown to predict thermodynamics–structure–dynamics “correlations” based on various long and short wavelength scalar properties all related to static two-point collective density fluctuations. The logarithm of the alpha relaxation time scales as a power law with these scalar metrics with an exponent that is significantly lower in the less dense noncooperative activated regime compared to the very dense highly cooperative regime. However, the discovered correlation of activated relaxation with a thermodynamic property (dimensionless compressibility) is not causal in ECNLE theory, but rather reflects a strong connection between the local structural quantities that quantify kinetic constraints in the theory with the amplitude of long wavelength density fluctuations. Fourth, the consequences of chemically specific nonuniversalities associated with the onset condition and relative importance of collective elasticity are studied. In conclusion, the predicted thermodynamics–structure–dynamics correlations are found to be robust, albeit with nontrivial shifts of the onset condition.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
Structure, dynamics, and electrochemistry of choline chloride/ethylene glycol eutectic solvents at an electrode surface explored by molecular dynamics simulations
Choline chloride and ethylene glycol mixtures with 1:2, 1:4, and 1:6 molar ratios on the surfaces of graphite and gold electrodes were studied using classical molecular dynamics simulations. Both neutral and charged electrodes were considered. The liquid composition, solvation structure, molecular orientation, and dynamics at the electrode surface are significantly different from those of the bulk liquid. These properties strongly depend on the electrode material and charge density, whereas they are less sensitive to the overall solvent composition. The effect of the electrode on the composition, structure, and orientation of the liquid fades beyond ∼10 Å from the surface of the electrode. This distance corresponds to about two layers of the structured solvent, despite the fact that the layered structure extends to at least five layers or about 25 Å. However, the electrode influences solvent dynamics over a longer distance. The electrochemical properties of the eutectic solvent at both electrode surfaces were also studied. The simulations captured the experimental differential capacitance shapes for both electrode systems, although the magnitudes and exact shapes differ. The simulations further revealed that the solvent in the first solvation layer can both overscreen and underscreen the electrode charges depending on the electrode material and electrode potential.
37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗