Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Solvation models”

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.

At least 91 records · Page 5
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

ezAlign: A Tool for Converting Coarse-Grained Molecular Dynamics Structures to Atomistic Resolution for Multiscale Modeling

Soft condensed matter is challenging to study due to the vast time and length scales that are necessary to accurately represent complex systems and capture their underlying physics. Multiscale simulations are necessary to study processes that have disparate time and/or length scales, which abound throughout biology and other complex systems. Herein we present ezAlign, an open-source software for converting coarse-grained molecular dynamics structures to atomistic representation, allowing multiscale modeling of biomolecular systems. The ezAlign v1.1 software package is publicly available for download at github.com/LLNL/ezAlign. Its underlying methodology is based on a simple alignment of an atomistic template molecule, followed by position-restraint energy minimization, which forces the atomistic molecule to adopt a conformation consistent with the coarse-grained molecule. The molecules are then combined, solvated, minimized, and equilibrated with position restraints. Validation of the process was conducted on a pure POPC membrane and compared with other popular methods to construct atomistic membranes. Additional examples, including surfactant self-assembly, membrane proteins, and more complex bacterial and human plasma membrane models, are also presented. By providing these examples, parameter files, code, and an easy-to-follow recipe to add new molecules, this work will aid future multiscale modeling efforts.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2026

Deciphering the Solvation Structure of Aqueous ZnCl 2 Solutions from X-ray Absorption Spectra Using the Interpretable Graph Neural Network

Machine learning (ML) provides powerful pathways for predicting spectroscopic observables from atomic structures, but its broader impact depends on making model predictions interpretable in terms of physical and chemical principles. Here, we introduce a physics-guided graph neural network (GNN) model that predicts Zn K-edge X-ray spectroscopy (XAS) spectra of aqueous ZnCl 2 solutions. Training data are generated from ab initio XAS calculations on molecular dynamics snapshots obtained using a machine learning interatomic potential. The GNN reproduces experimental spectra across concentrations from dilute (<0.1 m) to highly concentrated (30 m, “water-in-salt”) regimes and scales efficiently to large, disordered liquid systems beyond the reach of conventional ab initio approaches. Gradient-based attribution analysis reveals that the model learns physically meaningful structure-spectrum relationships. Ligand-specific attributions reflect orbital hybridization patterns and the origin of the excitations derived from the density functional theory. Bond-length attributions recover spectral shifts consistent with multiple-scattering theory. Finally, this work bridges data-driven prediction with electronic-structure theory, establishing a general paradigm for interpretable ML that links atomic structure, electronic structure, and spectroscopic observables.

25 ENERGY STORAGE↗
DOE OSTI2025

Microscopic insights into the solvation of polyethylene glycol chains in water: A machine learning potential approach

Polyethylene glycol (PEG) is a structurally simple, nontoxic, and water-soluble polymer widely utilized in medical and pharmaceutical applications. Notably, when a PEG chain is immersed in water, the surrounding water molecules play a key role in driving conformational changes of this macromolecule. In this study, we explore the solvation behavior of PEG under mechanical strain using molecular dynamics simulations, with an interatomic potential obtained from machine learning. Our focus is on the transition from the favored coil-like conformation to an extended one under external force. Through analyses of radial distribution functions, hydrogen bonding, and solvation dynamics, we uncover how mechanical stretching influences the local hydration environment. Furthermore, we disentangle the enthalpic and entropic contributions to the conformational stability of PEG in water. Surprisingly, our neural network potential model identifies dewetting of PEG C-atoms, and not water H-bonding with PEG O-atoms, as the main enthalpic driving force for the coiling of PEG in water.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
DOE OSTI2024

Zn 2+ -mediated catalysis for fast-charging aqueous Zn-ion batteries

Rechargeable aqueous zinc-ion batteries (AZIBs), renowned for their safety, high energy density and rapid charging, are prime choices for grid-scale energy storage. Historically, ion-shuttling models centring on ion-migration behaviour have dominated explanations for charge/discharge processes in aqueous batteries, like classical ion insertion/extraction and pseudocapacitance mechanisms. However, these models struggle to account for the exceptional performance of AZIBs compared to other aqueous metal-ion batteries. Here, in this study, we present a catalysis model elucidating the Zn 2+ anomaly in aqueous batteries, explaining it through the concept of adsorption in catalysis. Such behaviour can serve the charge/discharge role, predominantly dictated by solvated metal cations and cathode materials. First-principles calculations suggest optimal adsorption/desorption behaviour (water dissociation process) with the Zn 2+ -vanadium nitride (VN) combination. Experimentally, AZIBs implementing VN cathodes demonstrate fast-charging kinetics, showing a capacity of 577.1 mAh g -1 at a current density of 300,000 mA g -1 . The grasp of catalysis steps within AZIBs can drive solutions beyond state-of-the-art fast-charging batteries.

25 ENERGY STORAGE↗
DOE OSTI2025

Machine Learning-Driven Solvent Screening for Biobased 2,3-Butanediol Extraction

Biobased 2,3-butanediol (2,3-BDO) is a valuable biomass-derived chemical due to its versatility in being transformed into a wide variety of products. However, the separation and purification of 2,3-BDO from fermentation broth remain a significant challenge owing to its high boiling point and hydrophilic nature. Herein, we developed a machine learning (ML)-based screening workflow that uses molecular calculations as training data and requires only a small number of experimental measurements for validation to identify alternative solvent candidates for the liquid–liquid extraction (LLE) of 2,3-BDO from aqueous solution. In particular, 130 density functional theory (DFT) calculations with the implicit solvation method not only built a correlation between the computational partition coefficient and the experimental distribution coefficient of 2,3-BDO but also parameterized an Extra-Trees ML model to screen the distribution coefficient for a wider range of 6717 organic solvents. The experimental measurements of only 24 solvents were needed to validate the computational results. A list of 50 prioritized solvents was proposed for 2,3-BDO LLE, and seven additional experimental measurements were conducted to further verify our selected solvents. The impact of the extraction temperature and solvent-to-feed ratio was also investigated for selected solvents in experiments. Furthermore, this work suggested alternative solvents for 2,3-BDO LLE and proposed a versatile workflow that requires fewer experiments and can be applied to a broader range of LLE studies.

Extraction↗
NASA NTRS2004

Computational studies of Ras and PI3K

Until recently, experimental techniques in molecular cell biology have been the primary means to investigate biological risk upon space radiation. However, computational modeling provides an alternative theoretical approach, which utilizes various computational tools to simulate proteins, nucleotides, and their interactions. In this study, we are focused on using molecular mechanics (MM) and molecular dynamics (MD) to study the mechanism of protein-protein binding and to estimate the binding free energy between proteins. Ras is a key element in a variety of cell processes, and its activation of phosphoinositide 3-kinase (PI3K) is important for survival of transformed cells. Different computational approaches for this particular study are presented to calculate the solvation energies and binding free energies of H-Ras and PI3K. The goal of this study is to establish computational methods to investigate the roles of different proteins played in the cellular responses to space radiation, including modification of protein function through gene mutation, and to support the studies in molecular cell biology and theoretical kinetics models for our risk assessment project.

Ren, Lei↗
DOE OSTI2025

Unlocking enhanced gas capture via core scrambling of porous-organic cages

The demand for low-cost, low-energy, and highly selective gas capture and separations is an ongoing driver of porous material development. Porous liquids have been identified as a promising gas separation material by creating permanent porosity in inorganic solvents through inclusion of nanoporous materials that sterically exclude solvent from their internal porosity. Among the nanoporous materials that can be used to form porous liquids, porous-organic cages (POCs) have been one of the most popular due to the inherent tunability of POCs. “Scrambled” POCs with varying functionalities on the POC vertices have been developed and incorporated into porous liquid compositions, increasing their gas adsorption capacity. An unexplored avenue to tailor the properties of porous liquids is through scrambling the functionality of the core of the POC. Here, therefore, we have synthesized a new POC, a CC3-OH derivative with scrambled hydroxides on the core and evaluated the impact on the CO 2 uptake capacity in silicon oil-based porous liquids. Core scrambling of the POC resulted in a twofold increase CO 2 adsorption capacity in the porous liquid, an emergent property that is a dramatic increase beyond a linear combination of the gas adsorption capacity of the neat solvent and the POC. Density functional theory modeling of the CC3 POC and its hydroxide-based derivatives identified that free rotation of the linker hydroxide allowed for forced interaction between the CO 2 molecule and the hydroxide in the pore window. Solvation of the POC may release scrambled core hydroxides from intramolecular bonding with a neighboring imine, allowing for increased gas uptake in the porous liquid over the neat POC. These results identify a key structural relationship of POCs that enables emergent properties in porous liquids and can guide future development of liquid phase gas capture and separation materials for environmental and industrial applications.

Gas capture↗
DOE OSTI2026

Operando Depth-Resolved Measurement of Solvation Entropy, Interfacial Transport, and Charge-Transfer Kinetics in Lithium-Ion Batteries

Understanding and improving the performance and longevity of lithium-ion batteries critically depends on insight into the dynamic processes occurring at buried electrode-electrolyte interfaces. However, direct, depth-resolved, and operando diagnosis of these interfaces remains a longstanding challenge due to their inaccessibility beneath bulk materials, the limitations of conventional surface- and bulk-sensitive characterization tools, and the difficulty of maintaining realistic cell environments during measurement. These challenges have made it nearly impossible to uniquely resolve important interfacial properties such as charge transfer resistance, SEI (solid electrolyte interphase) resistance, and solvation entropy at the individual electrode interfaces within a working cell, information that is essential for mechanistic insight and accelerated battery design. Here, in this study, we report the development of multiharmonic electro-thermal spectroscopy (METS), an operando technique that enables depth-resolved measurement of solvation entropy, interfacial transport resistance, charge-transfer resistance, and SEI resistance at individual electrode-electrolyte interfaces within practical lithium-ion batteries. By leveraging frequency-dependent, thermal-wave sensing and interface-specific modeling, METS uniquely attributes interfacial properties to specific electrodes, as validated by comparison with traditional electrochemical impedance spectroscopy (EIS). The ability to spatially and temporally resolve interfacial processes in real time provides new diagnostic capabilities that are crucial for mechanistic studies of battery degradation and for the rapid development of next-generation energy storage systems.

Chalise, Divya [University of California, Berkeley↗
DOE OSTI2025

Environmental photosensitizers can exhibit enhanced actinic absorption in microhydrated clusters compared to solution

Brown carbon chromophores at environmental air-water interfaces often act as photosensitizers that absorb sunlight and subsequently transfer energy to nearby molecules, initiating a wide variety of chemical reactions. Despite their importance to understanding daytime chemistry at these airwater interfaces, little is known about the role of the solvation environment on the photophysical properties of these photosensitizers. In this work, we present a joint experimental-theoretical study of the vibrational and photophysical properties of microhydrated protonated and deprotonated 4-benzoylbenzoic acid (4-BBA), a key model system for environmental photosensitizers. We find that for protonated 4-BBAH + ∙(H 2 O) 0-1 , representing photosensitizers in very acidic conditions, a single bright state dominates the UV-Vis spectrum between 280 and 400 nm. Comparing the experimental UV-Vis spectra and quantum chemistry-predicted spectra of 4-BBA + ∙(H 2 O) 0-2 , we find that the degree of microhydration has little effect on the UV-Vis spectra or the orbitals of the dominant feature. For deprotonated 4-BBA ‒ , representing photosensitizers in basic conditions, quantum chemistry calculations predict that the UV-Vis spectra are ~3x weaker in intensity than the brightest 4-BBAH + ∙(H 2 O) 0-1 features and were not observed experimentally. Quantum chemistry calculations predict one dominant UV-Vis feature is present in the spectra of 4-BBA ‒ ∙(H 2 O) 0-2 , which exhibit minor shifts with degree of microhydration. While 4-BBA in bulk solution over a range of pH values has relatively weak absorption within the solar actinic region, we show that microhydrated 4-BBA has bright transitions within the actinic region. This indicates that the complex structure of environmental air-water interfaces can shift the absorption maximum of photosensitizers into the actinic region for enhanced absorption of sunlight and subsequent enhancement of photosensitizer-driven reactions.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
NASA NTRS2011

Ion-Ice Astrochemistry: Barrierless Low-Energy Deposition Pathways to HCOOH, CH3OH, and CO2 on Icy Grain Mantles from Precursor Cations

A new family of very favorable reaction pathways is explored involving the deposition of ions on icy grain mantles with very low energies. Quantum chemical cluster calculations at the MP2/6-31+G** level in 4H2O clusters and at the B3LYP/6-31+G** level in 17H2O clusters indicate that HCO+ and CH3 + are able to react spontaneously with one of the water molecules in the cluster to form protonated formic acid (HCOOH2 +) and protonated methanol (CH3OH2 +), respectively. It is furthermore found that these initial adducts spontaneously transfer their excess protons to the cluster to form neutral formic acid and methanol, plus solvated hydronium, H3O+. In the final case, if a CO molecule is bound to the surface of the cluster, OH+ may react with it to form protonated carbon dioxide (HCO2 +), which then loses its proton to yield CO2 and H3O+. In the present model, all of these processes were found to occur with no barriers. Discussion includes the analogous gas phase processes, which have been considered in previous studies, as well as the competitive abstraction pathway for HCO(+) + H2O.

Woon, David E.↗
NASA NTRS2002

Protein Solvation in Membranes and at Water-Membrane Interfaces

Different salvation properties of water and membranes mediate a host of biologically important processes, such as folding, insertion into a lipid bilayer, associations and functions of membrane proteins. These processes will be discussed in several examples involving synthetic and natural peptides. In particular, a mechanism by which a helical peptide becomes inserted into a model membrane will be described. Further, the molecular mechanism of recognition and association of protein helical segments in membranes will be discussed. These processes are crucial for proper functioning of a cell. A membrane-spanning domain of glycophorin A, which exists as a helical dimer, serves as the model system. For this system, the free energy of dissociation of the helices is being determined for both the wild type and a mutant, in which dimerization is disrupted.

Pohorille, Andrew↗
DOE OSTI2024

DRIFT: Diffusivity Regulation of Isotopes by nanoconFinementT

Spontaneous isotope fractionation has been reported under nanoconfinement conditions in naturally occurring systems, but the origin of this phenomena is currently unknown. Two existing hypotheses have been proposed, one based on changes in the solvation environment of the isotopes that reduces the non-mass dependent hydrodynamics contribution to diffusion. The other is that isotopes have mass-dependent surface adsorption, varying their total diffusion through nanoconfined channels. To investigate these hypotheses, benchtop experiments, nuclear magnetic resonance (NMR) spectroscopy, and molecule scale modeling were applied. Classical molecular dynamics simulations identified that the Na + and Cl - hydration shells across the three different salt solutions ( 22 Na 35 Cl, 23 Na 35 Cl, 24 Na 35 Cl) did not vary as a function of the Na + isotope, but that there was a significant pore size effect, with larger hydration shells at larger pore sizes. Additionally, while total adsorption times did not vary as a function of the Na + isotope or pore size, the free ion concentration, or those adsorbed on the surface for <5% of the simulation time did exhibit isotope dependence. Experimentally, challenges occurred developing a repeatable experiment, but NMR characterization of water diffusion rates through ordered alumina membranes was able to identify the existence of two distinct water environments associated with water inside and outside the pore. Further NMR studies could be used to confirm variation in hydration shells and diffusion rates of dissolved ions in water. Ultimately, mass-dependence adsorption is a primary driver of variations in isotope diffusion rates, rather than variation in hydration shells that occur under nanoconfinement.

73 NUCLEAR PHYSICS AND RADIATION PHYSICS↗
NASA NTRS1994

Interaction of a Model Peptide with a Water--Bilayer System

We discuss a molecular dynamics study of the alanine dipeptide at the interface between water and a glycerol-1-monooleate (GMO) bilayer. The dipeptide is interfacially active and incorporates into the bilayer without disrupting its structure. The interfacial region that is readily penetrated by the dipeptide spans the entire head group portion of the bilayer. The polar groups of the alanine dipeptide mostly remain well solvated by water and the oxygen atoms of GMO, and conformations of the dipeptide are characterized by (phi, psi) angles typical of alpha-helix and beta-sheet structures. When the molecule is deeper in the bilayer, the C(sub 7eq) state also becomes stable. The barrier to the isomerization reaction at the interface is lower than in bulk phases. After 7 ns of trajectories, the system is not fully equilibrated, due to slow collective motions involving GMO head groups. These result in decreased mobility and lower rates of isomerization of the dipeptide at the interface.

Pohorille, A.↗
DOE OSTI2026

Selectivity in gas–liquid interactions: Molecular beam scattering of CD 4 and ND 3 from an aqueous flat liquid jet

The dynamics of polar and nonpolar molecules colliding with an aqueous surface are characterized by scattering molecular beams of deuterated methane and ammonia, CD 4 and ND 3 (E i = 28.9 and 30.3 kJ mol −1 , respectively), from a flat liquid jet of cold salty water (8 m LiBr, 230 K). Translational energy distributions of scattered species collected as a function of collision geometry probe both impulsive scattering (IS) and thermal desorption (TD) mechanisms. Here, we find that CD 4 scattering is dominated by IS and exhibits a super-specular angular distribution. The fraction of TD scattering events is notably smaller for cold salty water than for dodecane, consistent with a higher free energy of solvation for CD 4 in the water jet. In contrast, no scattering signal is seen for ND 3 from the water jet, a result attributed to the high solubility and efficient protonation of ND 3 in liquid water. The IS channel for CD 4 was analyzed using a soft-sphere model, yielding a higher internal energy (E int ) and lower effective surface mass (m eff ) than was seen for Ne/water; the higher value of E int is attributed to rotational excitation of the scattered CD 4 . These findings demonstrate that the outcomes of a gas–liquid collision—scattering trajectory, surface adherence, and energy transfer—are directed at the molecular level by both the gaseous scatterer and liquid surface.

Foreman, Madison M. [University of California, Ber↗
DOE OSTI2024

Computing the solubility of argon and xenon in molten sodium chloride and potassium chloride salts

Molten salt reactors (MSRs) offer significant advancements in nuclear reactor safety and efficiency by operating at higher temperatures and lower pressures compared to traditional reactors. A critical aspect of MSR operation involves understanding the solubility of fission byproducts, particularly noble gases, in the molten salts used. This study employs molecular dynamics (MD) simulations to compute Henry’s law constants and enthalpies of solvation for argon and xenon in molten sodium chloride (NaCl) and potassium chloride (KCl). We developed a new pairwise potential for the noble gas and salt interactions based on first principles calculations. We then used this potential to calculate Henry’s law constants of the two gases in the molten salts, which were modeled using both a rigid ion model (RIM) and a polarizable ion model (PIM). The solubility calculations, performed using the Widom insertion method, show qualitative agreement with limited experimental data, highlighting the temperature dependence and greater solubility of both gases in KCl compared to NaCl. Additionally, free volume analysis elucidated the role of available space within the molten salts in governing solubility trends. Our findings suggest that PIM trajectories provide more reliable predictions for noble gas solubility than RIM due to their accurate density representation. Furthermore, these results enhance understanding of gas solubility in MSR environments, and the methods can be readily extended to other systems.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
NASA NTRS2013

The First X-ray Diffraction Patterns of Clay Minerals from Gale Crater

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 021 band consistent with a trioctahedral phyllosilicate. A broad peak at approx 10A with a slight inflexion at approx 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60degC in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 degC might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

Bristow, Thomas↗
NASA NTRS2014

Paleoenvironmental Implications of Clay Minerals at Yellowknife Bay, Gale Crater, Mars

The Mars Science Laboratory (MSL) Rover, Curiosity spent approx 150 sols at Yellowknife Bay (YKB) studying a section of fluvio-lacustrine sedimentary rocks (with potential indications of volcanic influence), informally known as the Yellowknife Bay formation. YKB lies in a distal region of the Peace Vallis alluvial fan, which extends from the northern rim of Gale Crater toward the dune field at the base of Mt Sharp. Sedimentological and stratigraphic observations are consistent with the Yellowknife Bay formation being part of a distal fan deposit, which could be as young as middle Hesperian to even early Amazonian in age (approx. 3.5 to 2.5 Ga). The Yellowknife Bay formation hosts a unit of mudstone called the Sheepbed member. Curiosity obtained powdered rock samples from two drill holes in the Sheepbed Member, named John Klein and Cumberland, and delivered them to instruments in Curiosity. Data from CheMin, a combined X-ray diffraction (XRD)/X-ray fluorescence instrument (XRF), has allowed detailed mineralogical analysis of mudstone powders revealing a clay mineral component of approx. 20 wt.% in each sample. The clay minerals are important indicators of paleoenvironmental conditions and sensitive recorders of post-depositional alteration processes. The XRD pattern of John Klein reveals a 02l band consistent with a trioctahedral phyllosilicate. A broad peak at approx. 10A with a slight inflexion at approx. 12A indicates the presence of 2:1 type clay minerals in the John Klein sample. The trioctahedral nature of the clay minerals, breadth of the basal reflection, and presence of a minor component with larger basal spacing suggests that John Klein contains a trioctahedral smectite (probably saponite), whose interlayer is largely collapsed because of the low-humidity conditions. The XRD patterns show no evidence of corrensite (mixed-layer chlorite/smectite) or chlorite, which are typical diagenetic products of trioctahedral smectites when subjected to burial and heating >60 C in the presence of water. Given estimated geothermal gradients on Mars temperatures <60 C might still be consistent with (but do not require) moderate burial. However, our ability to identify interstratified minerals is greatly limited by the lack of access to traditional treatments methods used in the lab (e.g., ethylene glycol solvation). Our preferred explanation for the origin of trioctahedral smectites in Sheepbed mudstone is in situ production via reaction of olivine, water and Si-bearing amorphous material, an important mudstone component detected by XRD. Elevated levels of magnetite in the Sheepbed and the trioctahedral monomineralic nature of the clay minerals support this model. These observations, combined with previous studies of olivine stability, support the persistence of circum-neutral hydrous conditions for thousands of years at YKB.

MSL↗
DOE OSTI2026

MixPI: Mixed-time slicing path integral software for quantized molecular dynamics simulations

We introduce the MixPI software to implement path integral molecular dynamics (PIMD) simulations for the study of condensed phase systems where nuclear quantum effects (NQEs) are important. In contrast to existing PIMD simulation software, MixPI enables the implementation of mixed quantum–classical path integral simulations where only a subset of system degrees of freedom (dofs) are treated quantum mechanically in an extended phase space while the remaining dofs are described classically. We expect this software to be particularly useful for simulations of electron and proton transfer in condensed phase systems, as well as for the study of biological and material systems where only a handful of dofs contribute significantly to the observed NQEs. We demonstrate the use of MixPI in two different systems. The first is a simple water model where we implement a set of mixed quantum–classical simulations to compute average energy and radial distribution functions. We use these simulations to benchmark the effectiveness of MixPI and to demonstrate how it enables systematic investigation into the origin of observed NQEs. We then compute radial distribution functions for a system where MixPI is essential: a solvated metal (M 2+ ) cation described using an explicit quantized electron localized on an M 3+ ion in water.

chemical physics↗