SEARCH · Search NASA
Results for “Quantum molecular dynamics simulation”
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.
Development of neural network force fields for corrosion studies
To fully understand the chemistry and physics of corrosion, novel methods of simulation must be developed. One approach is designing machine learning (ML) algorithms integrated with density functional theory to develop adaptive force fields to gain insight into corrosion behavior namely at the surface of metal oxides. Current methods of modeling corrosion are slow due to the computational cost of resolving both reaction mechanics and mass transport processes. Machine learning methods can be implemented to obtain structure-activity relationships at both the molecular and bulk scale while still retaining the accuracy of density functional theory (DFT) and significantly decreasing the time needed for simulations of complex chemical processes in the various environments of corrosion. Multiscale models are needed for corrosion studies to fully understand its processes not only at the atomic length scale (chemical bonding, energies, and forces), but also at the nano and meso length scales (solid-state physics and material science processes). Current methods of study include DFT, molecular dynamics, and Monte Carlo. The limitation of DFT is that only a small number of atoms or molecules can be simulated at that level of theory. Density functional theory is used to study the electronic structure of atoms and molecules, and calculate the force component of each atom. However, these calculations are limited to about 1000 atoms. Custom periodic boundary conditions (PBC) can be used to describe the various environments and defects that affect the atomic forces to produce a large data set from which a training set can be derived. Machine learning can be utilized to overcome the barrier of modeling macroscopic and multi-scale processes from ab initio calculations through the development of adaptive force fields. Local environments determine the atomic forces of a given system, therefore adaptive force fields must be created to produce reliable quantum mechanical calculations. This can be achieved by developing a learning algorithm that uses the mapped atomic forces or fingerprint as an input to produce energies and magnetic moments as output. A systematic approach was used to begin to build a data set in order to accurately describe the atomic forces in various environments. In Figure 4 below, a simple PBC cell of Fe{sub 2}O{sub 3} was first optimized. A surface optimization was performed next, followed by a hydroxylated surface optimization. Once this calculation has converged, the adsorption of halide species to the hydroxylated surface will be investigated. TensorFlow is an open source platform for machine learning developed by Google. Using a high level application program interface (API) such as Keras allows for building and training ML models easily in a number of different environments and languages. For this project, a neural network was developed within Anaconda in Python. Future Work: Further development of reference data set; Refining neural network and learning algorithm; Fingerprinting atomic environment to enable mapping of atomic force components; Choosing appropriate training set from reference data; Learning from training set and enabling non-linear mapping of training set fingerprints and the atomic forces; Estimation of uncertainty to identify ranges of outside applicability; Testing and analysis of molecular dynamic simulations.
Path sampling of recurrent neural networks by incorporating known physics
Recurrent neural networks have seen widespread use in modeling dynamical systems in varied domains such as weather prediction, text prediction and several others. Often one wishes to supplement the experimentally observed dynamics with prior knowledge or intuition about the system. While the recurrent nature of these networks allows them to model arbitrarily long memories in the time series used in training, it makes it harder to impose prior knowledge or intuition through generic constraints. In this work, we present a path sampling approach based on principle of Maximum Caliber that allows us to include generic thermodynamic or kinetic constraints into recurrent neural networks. We show the method here for a widely used type of recurrent neural network known as long short-term memory network in the context of supplementing time series collected from different application domains. These include classical Molecular Dynamics of a protein and Monte Carlo simulations of an open quantum system continuously losing photons to the environment and displaying Rabi oscillations. Our method can be easily generalized to other generative artificial intelligence models and to generic time series in different areas of physical and social sciences, where one wishes to supplement limited data with intuition or theory based corrections.
Fock-Space Schrieffer–Wolff Transformation: Classically-Assisted Rank-Reduced Quantum Phase Estimation Algorithm
We present an extension of many-body downfolding methods to reduce the resources required in the quantum phase estimation (QPE) algorithm. In this paper, we focus on the Schrieffer–Wolff (SW) transformation of the electronic Hamiltonians for molecular systems that provides significant simplifications of quantum circuits for simulations of quantum dynamics. We demonstrate that by employing Fock-space variants of the SW transformation (or rank-reducing similarity transformations (RRST)) one can significantly increase the locality of the qubit-mapped similarity-transformed Hamiltonians. The practical utilization of the SW-RRST formalism is associated with a series of approximations discussed in the manuscript. In particular, amplitudes that define RRST can be evaluated using conventional computers and then encoded on quantum computers. The SW-RRST QPE quantum algorithms can also be viewed as an extension of the standard state-specific coupled-cluster downfolding methods to provide a robust alternative to the traditional QPE algorithms to identify the ground and excited states for systems with various numbers of electrons using the same Fock-space representations of the downfolded Hamiltonian. The RRST formalism serves as a design principle for developing new classes of approximate schemes that reduce the complexity of quantum circuits.
Diabatization with Electrostatic Embedding for Studying Photophysics in Organic Molecular Crystals
Highly emissive organic molecular crystals find applications in several areas, such as organic electronics, solar cells, and sensors. Understanding the excited-state mechanisms underlying these applications is essential for optimizing and controlling them effectively. Exciton models coupled with nonadiabatic dynamics, particularly quantum dynamics, provide crucial insights into photochemical and photophysical processes in molecular crystals. Nevertheless, there remains a lack of general tools and automated workflows to facilitate such simulations. In this paper, we present a computational strategy to investigate the photoactivated dynamics of organic molecular crystals, bridging methodologies traditionally used for molecular systems and materials science, with a particular focus on the interplay between local excitations and charge transfer (CT) processes. We have implemented an interface between the fromage and Overdia programs, enabling the construction of vibronic Hamiltonians for molecular crystals within an excited-state ONIOM(QM:QM′) framework, incorporating long-range electrostatics through a RESP-based Ewald summation. Fragment-based diabatization provides a route to quantum dynamics simulations in weak-to-intermediate coupling regimes. The method was applied to the photophysics of dibenzo[g,p]chrysene (DBC) crystals using time-dependent DFT. The fromage/ Overdia interface was employed to compute the couplings of local excitations and CT states for 18 unique DBC dimers in the crystal and to quantify the influence of electrostatic embedding, which was found to be modest (10−20%). Simulations on π-stacked dimers reproduced the small red shift observed experimentally from solution to crystal, attributed to electronic interactions among fixed monomers rather than crystal electrostatics. Quantum dynamics simulations revealed ultrafast population transfer from bright local excitations to CT states. This approach establishes a robust framework linking molecular and solid-state excited-state dynamics, with potential applications for studying excitations, defects, and impurities in molecular crystals.
Breaking the Million-Electron and 1 EFLOP/s Barriers: Biomolecular-Scale Ab Initio Molecular Dynamics Using MP2 Potentials
The accurate simulation of complex biochemical phenomena has historically been hampered by the computational requirements of high-fidelity molecular-modeling techniques. Quantum mechanical methods, such as ab initio wave-function (WF) theory, deliver the desired accuracy, but have impractical scaling for modeling biosystems with thousands of atoms. Combining molecular fragmentation with MP2 perturbation theory, this study presents an innovative approach that enables biomolecular-scale ab initio molecular dynamics (AIMD) simulations at WF theory level. Leveraging the resolution-of-the-identity approximation for Hartree-Fock and MP2 gradients, our approach eliminates computationally intensive four-center integrals and their gradients, while achieving near-peak performance on modern GPU architectures. The introduction of asynchronous time steps minimizes time step latency, overlapping computational phases and effectively mitigating load imbalances. Utilizing up to 9,400 nodes of Frontier and achieving 59% (1006.7 PFLOP/s) of its double-precision floating-point peak, our method enables us to break the million-electron and 1EFLOP/s barriers for AIMD simulations with quantum accuracy.
Quantum Simulations of Vibrational Strong Coupling via Path Integrals
A quantum simulation of vibrational strong coupling (VSC) in the collective regime via thermostated ring-polymer molecular dynamics (TRPMD) is reported. For a collection of liquid-phase water molecules resonantly coupled to a single lossless cavity mode, the simulation shows that as compared with a fully classical calculation, the inclusion of nuclear and photonic quantum effects does not lead to a change in the Rabi splitting but does broaden polaritonic line widths roughly by a factor of 2. Moreover, under thermal equilibrium, both quantum and classical simulations predict that the static dielectric constant of liquid water is largely unchanged inside vs outside the cavity. Furthermore, this result disagrees with a recent experiment demonstrating that the static dielectric constant of liquid water can be resonantly enhanced under VSC, suggesting either limitations of our approach or perhaps other experimental factors that have not yet been explored.
Computational methods based on density functional theory for reactions and processes involving electronic spin (Final Technical Report)
This award supports one post-doctoral researcher for 1.5 years. Publications that acknowledge this grant: Refs. 1–14. Refs. 3,10,12 assess current methodology for the evaluation of magnetic exchange couplings in transition metal complexes. In particular, Ref. 10 validates the use of an approximate (non-iterative) Green’s function approach for the calculation of magnetic exchange couplings and will be the foundation for Thrust 2 in this proposal. Refs. 3 and 12 focus on widely used density functional approaches based on the standard energy differences methodology for the particular case of oxo-bridged Fe(III) complexes. Refs. 2,4–7,11 apply current methodologies to problems of practical interest in molecular magnetism. Ref. 13 presents a methodology to explicitly simulate the dynamics of open quantum systems within density functional theory (DFT) calculations based on the Liouville-von Neumann equation of motion for quantum systems driven out-of-equilibrium. Ref. 8 uses non-collinear spin DFT to explain the mechanical behavior of magnetic mono-atomic Pt wires produced in break-junction experiments in the presence of a magnetic field.
Failure of molecular dynamics to provide appropriate structures for quantum mechanical description of the aqueous chloride ion charge-transfer-to-solvent ultraviolet spectrum
The lowest band in the charge-transfer-to-solvent ultraviolet absorption spectrum of aqueous chloride ion is studied by experiment and computation. Interestingly, the experiments indicate that at concentrations up to at least 0.25 M, where calculations indicate ion pairing to be significant, there is no notable effect of ionic strength on the spectrum. The experimental spectra are fitted to aid comparison with computations. Classical molecular dynamic simulations are carried out on dilute aqueous Cl - , Na + , and NaCl, producing radial distribution functions in reasonable agreement with experiment and, for NaCl, clear evidence of ion pairing. Clusters are extracted from the simulations for quantum mechanical excited state calculations. Accurate ab initio coupled-cluster benchmark calculations on a small number of representative clusters are carried out and used to identify and validate an efficient protocol based on time-dependent density functional theory. The latter is used to carry out quantum mechanical calculations on thousands of clusters. The resulting computed spectrum is in excellent agreement with experiment for the peak position, with little influence from ion pairing, but is in qualitative disagreement on the width, being only about half as wide. It is concluded that simulation by classical molecular dynamics fails to provide an adequate variety of structures to explain the experimental CTTS spectrum of aqueous Cl - .
Imaging multiphoton ionization dynamics of CH 3 I at a high repetition rate XUV free-electron laser
XUV multiphoton ionization of molecules is commonly used in free-electron laser experiments to study charge transfer dynamics. However, molecular dissociation and electron dynamics, such as multiple photon absorption, Auger decay, and charge transfer, often happen on competing time scales, and the contributions of individual processes can be difficult to unravel. We experimentally investigate the Coulomb explosion dynamics of methyl iodide upon core–hole ionization of the shallow inner-shell of iodine (4d) and classically simulate the fragmentation by phenomenologically introducing ionization dynamics and charge transfer. Under our experimental conditions with medium fluence and relatively long XUV pulses (~75 fs), we find that fast Auger decay prior to charge transfer significantly contributes to the charging mechanism, leading to a yield enhancement of higher carbon charge states upon molecular dissociation. Furthermore, we argue for the existence of another charging mechanism for the weak fragmentation channels leading to triply charged carbon atoms. This study shows that classical simulations can be a useful tool to guide the quantum mechanical description of the femtosecond dynamics upon multiphoton absorption in molecular systems.
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.
Active learning for robust, high-complexity reactive atomistic simulations
Machine learned reactive force fields based on polynomial expansions have been shown to be highly effective for describing simulations involving reactive materials. Nevertheless, the highly flexible nature of these models can give rise to a large number of candidate parameters for complicated systems. In these cases, reliable parameterization requires a well-formed training set, which can be difficult to achieve through standard iterative fitting methods. In this paper, we present an active learning approach based on cluster analysis and inspired by Shannon information theory to enable semi-automated generation of informative training sets and robust machine learned force fields. The use of this tool is demonstrated for development of a model based on linear combinations of Chebyshev polynomials explicitly describing up to four-body interactions, for a chemically and structurally diverse system of C/O under extreme conditions. We show that this flexible training database management approach enables development of models exhibiting excellent agreement with Kohn–Sham density functional theory in terms of structure, dynamics, and speciation.
Effective Fragment Potential-Based Molecular Dynamics Studies of Diffusion in Acetone and Hexane
To facilitate more reliable descriptions of transport properties in liquids, molecular dynamics (MD) simulations are performed based on the effective fragment potential (EFP) method derived from first-principles quantum mechanics (in contrast to MD based upon empirically fitted potentials). The EFP method describes molecular interactions in terms of Coulomb, polarization/induction, dispersion, exchange-repulsion, and charge-transfer interactions. The EFP MD simulations described in this paper, performed on hexane and acetone, are able to track the mean-square displacement of molecules for sufficient time to reliably extract translational diffusion coefficients. Finally, the results reported here are in reasonable agreement with experiment.
Melting curves of ice polymorphs in the vicinity of the liquid–liquid critical point
The possible existence of a liquid–liquid critical point in deeply supercooled water has been a subject of debate due to the challenges associated with providing definitive experimental evidence. The pioneering work by Mishima and Stanley [Nature 392, 164–168 (1998)] sought to shed light on this problem by studying the melting curves of different ice polymorphs and their metastable continuation in the vicinity of the expected liquid–liquid transition and its associated critical point. Based on the continuous or discontinuous changes in the slope of the melting curves, Mishima [Phys. Rev. Lett. 85, 334 (2000)] suggested that the liquid–liquid critical point lies between the melting curves of ice III and ice V. Herein we explore this conjecture using molecular dynamics simulations with a machine learning model based on ab initio quantum-mechanical calculations. We study the melting curves of ices III, IV, V, VI, and XIII and find that all of them are supercritical and do not intersect the liquid–liquid transition locus. We also find a pronounced, yet continuous, change in the slope of the melting lines upon crossing of the liquid locus of maximum compressibility. Finally, we analyze the literature in light of our findings and conclude that the scenario in which the melting curves are supercritical is favored by the most recent computational and experimental evidence. Although the preponderance of evidence is consistent with the existence of a second critical point in water, the behavior of ice polymorph melting lines does not provide strong evidence in support of this viewpoint, according to our calculations.
Ab initio molecular dynamics on quantum computers
Ab initio molecular dynamics (AIMD) is a valuable technique for studying molecules and materials at finite temperatures where the nuclei evolve on potential energy surfaces obtained from accurate electronic structure calculations. In this work, we present an approach to running AIMD simulations on noisy intermediate-scale quantum (NISQ)-era quantum computers. The electronic energies are calculated on a quantum computer using the variational quantum eigensolver (VQE) method. Algorithms for computation of analytical gradients entirely on a quantum computer require quantum fault-tolerant hardware, which is beyond NISQ-era. Therefore, we compute the energy gradients numerically using finite differences, the Hellmann-Feynman theorem, and a correlated sampling technique. This method only requires additional classical calculations of electron integrals for each degree of freedom without any additional computations on a quantum computer beyond the initial VQE run. As a proof of concept, AIMD simulations are demonstrated for the H-2 molecule on IBM quantum devices. In addition, we demonstrate the validity of the method for larger molecules using full configuration interaction wave functions. As quantum hardware and noise mitigation techniques continue to improve, the method can be utilized for studying larger molecular systems.
Macroscale properties and atomic-scale mechanisms of ash removal in low-temperature hydrothermal carbonization
Biogenic ash is a significant impediment to the utilization of agricultural residues in biofuel production. Such challenge can be addressed by various treatments, as demonstrated in this study on the experimental and computational mechanisms involved in the hydrothermal treatment (HT) of wheat straw. A combination of classical (all-atom) molecular dynamics simulations of cellulose carrying silica and calcium species, along with first principles quantum chemical calculations, indicates the dissociation of inorganics from the cellulose with increased HT temperature. Here, this observation is confirmed by experimental evidence of effective ash removal by HT, showing at least 50% removal of sulfur, chlorine, potassium, and calcium, and 12.5% of silica, leading to a reduced total ash content (from 6.7% to 4.2%).Changes in structural features upon HT, such as surface cellular structure and porosity, were revealed, accompanied by an increased specific surface area (from 1.17 to 6.34 m 2 /g). Our simulations suggest that silica binds tightly to the hydrophobic face of cellulose at room temperature, but HT significantly reduces the binding free energy of association with both hydrophobic and hydrophilic surfaces. Most significantly, ash removal leads to an increased calorific value, rising from approximately 16 MJ/kg to about 19 MJ/kg, along with improved thermal behavior. The improved integration combustion index parameter S indicates that the combustion properties improve with ash removal efficiency. The proposed atomic-level mechanism for the observed removal of inorganics during mild HT underscores the potential of such treatment in producing energy-dense wheat straw, a widely available agricultural residue.
Ground and excited state gradients with end-to-end differentiable semiempirical quantum chemistry
Accurate and efficient gradients of molecular energy with respect to nuclear degrees of freedom are essential for geometry optimization and molecular dynamics, including simulations that go beyond the Born–Oppenheimer regime. A common approach involves deriving analytical formulas for new electronic structure methods, which is often conceptually difficult and requires tedious coding. Here, we implement analytical, semi-numerical, and automatic differentiation (AD)-based gradient pathways for semiempirical Hamiltonian models in the PYSEQM software package, leveraging both graphics processing unit (GPU) and central processing unit (CPU) architectures. We further extend these capabilities to excited states calculated using the configuration interaction singles and time-dependent Hartree–Fock ansätze. We benchmark wall time, peak memory usage, and accuracy across three molecular families of varying chemical complexity, including systems of up to a thousand atoms. For ground-state simulations, analytical and AD gradients achieve near-identical GPU runtimes, while semi-numerical gradients are slower on GPU but remain competitive on CPU. For excited states, both analytical and custom AD approaches using implicit differentiation show similar performance and low memory requirements, whereas gradients with full AD are memory-limited. AD gradients match analytical ones in accuracy across all tested systems, aided by a quaternion-based diatomic frame rotation for two-center quantities that ensures smooth energy surfaces. Overall, automatic differentiation emerges as a practical alternative to analytical gradients in semiempirical quantum chemistry, offering high accuracy while allowing seamless integration in AI-driven workflows and popular packages, such as PyTorch and JAX. Our results provide actionable guidance for selecting optimal gradient strategies in large-scale ground- and excited-state molecular dynamics simulations.
Neutron scattering and neural-network quantum molecular dynamics investigation of the vibrations of ammonia along the solid-to-liquid transition
Abstract Vibrational spectroscopy allows us to understand complex physical and chemical interactions of molecular crystals and liquids such as ammonia, which has recently emerged as a strong hydrogen fuel candidate to support a sustainable society. We report inelastic neutron scattering measurement of vibrational properties of ammonia along the solid-to-liquid phase transition with high enough resolution for direct comparisons to ab-initio simulations. Theoretical analysis reveals the essential role of nuclear quantum effects (NQEs) for correctly describing the intermolecular spectrum as well as high energy intramolecular N-H stretching modes. This is achieved by training neural network models using ab-initio path-integral molecular dynamics (PIMD) simulations, thereby encompassing large spatiotemporal trajectories required to resolve low energy dynamics while retaining NQEs. Our results not only establish the role of NQEs in ammonia but also provide general computational frameworks to study complex molecular systems with NQEs.