Search NASA⌕ Search

SEARCH · Search NASA

Results for “classical molecular dynamic simulations”

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 OSTI2025

Understanding Viscoelasticity of an Entangled Silicone Copolymer via Coarse-Grained Molecular Dynamics Simulations

Entangled dynamics is important for understanding rheological properties of long-chain polymers. For entangled homopolymers, the classic tube-reptation model and its refinements have been successfully applied to quantify properties like diffusion coefficient and zero-rate viscosity. However, the application of such models to copolymers has been limited despite scientific and industrial importance. Here, we study the entangled melt dynamics of poly(dimethyl-co-diphenyl)siloxane random copolymer for a range of mean-composition-ratio ϕ of the diphenyl component via long-term molecular dynamics simulation with a recently developed coarse-grained model. We found that the segmental relaxation is heterogeneous at the monomeric level because of compositional fluctuations. However, at the chain-entanglement level and higher length scales, the viscoelastic response is homogeneous with compositional dependence only through the overall diphenyl fraction ϕ. The relaxation modulus of the entangled copolymer melt conforms to the Likhtman–McLeish model, and the viscosity predicted using our current coarse-grained parameters is in good quantitative agreement with experimental data.

Copolymers↗
DOE OSTI2026

Prediction and Experimental Verification of Electrolyte Solvation Structure from an OMol25-Trained Interatomic Potential

A molecular-level understanding of electrolyte solvation structure and ion–ion correlations is critical to developing next-generation battery chemistries. Atomistic simulation capabilities with sufficient accuracy, speed, and transferability to deliver reliable structural insights while avoiding arduous system-specific reparameterization are thus highly desirable. Machine learning interatomic potentials (MLIPs) trained on large, chemically diverse data sets are revolutionizing computational chemistry, enabling molecular dynamics simulations of battery electrolytes with near-DFT accuracy over 10,000× faster than DFT. While previous MLIP training data sets with suitable elemental coverage for electrolytes have been based on inorganic materials, the Open Molecules 2025 (OMol25) data set provides large-scale molecular DFT MLIP training data with broad elemental coverage and specifically samples tens of millions of electrolyte configurations. Here, we integrate computational modeling with experimental validation to systematically assess the ability of large-scale MLIPs pretrained on materials data or on OMol25 to accurately resolve nanoscale structural organization and ion-solvation characteristics in Na-ion battery electrolytes across diverse physicochemical conditions and compositional regimes. We find that the OMol25-trained Universal Model of Atoms (UMA-OMol) predicts experimentally measured densities and X-ray structure factors in substantially better agreement compared to state-of-the-art models trained only on inorganic materials data. Using UMA-OMol, we further analyze systematic trends in solvation structure as a function of cation identity, anion chemistry, salt concentration, and solvent topology. We observe that increasing system temperature amplifies the heterogeneity within the solvation environment, perturbing cation–solvent interactions and promoting the formation of contact ion pairs (CIPs). Moreover, subtle variations in the solvent topology of glyme-based electrolytes cause pronounced changes in ion correlations and solvation structure. The experimental agreement and microscopic insights shown here position OMol25-trained MLIPs as a practical route to predictive, high-throughput electrolyte simulations beyond the limits of classical force fields and direct DFT molecular dynamics, serving as a powerful tool for accelerating the design of next-generation Na-ion battery electrolytes and beyond.

MLIPs↗
DOE OSTI2024

Impact of classical statistics on thermal conductivity predictions of BAs and diamond using machine learning molecular dynamics

Machine learning interatomic potentials (MLIPs) have greatly enhanced molecular dynamics (MD) simulations, achieving near-first-principles accuracy in thermal conductivity studies. In this work, we reveal that this accuracy, observed in BAs and diamond at sub-Debye temperatures, stems from an accidental error cancelation: classical statistics overestimates specific heat while underestimating phonon lifetimes, balancing out in thermal conductivity predictions. However, this balance is disrupted when isotopes are introduced, leading MLIP-based MD to significantly underpredict thermal conductivity compared to experiments and quantum statistics-based Boltzmann transport equation. This discrepancy arises not from classical statistics affecting phonon–isotope scattering rates but from its impact on the interplay between phonon–isotope and phonon–phonon scattering in the normal scattering-dominated BAs and diamond. In conclusion, this work underscores the limitations of MLIP-based MD for thermal conductivity studies at sub-Debye temperatures.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Cellular Phosphate Sensing and Anion Binding by an Azacrown‐Calixpyrrole Hybrid

A hybrid receptor-sensor for anions originating from the merging of positively charged ammonium moieties for electrostatic attraction/stronger binding of azacrowns with directionality of calixpyrrole hydrogen bond donors for selectivity is investigated. As demonstrated this hybrid receptor-sensor shows a remarkable selectivity for orthophosphate even in the presence of other phosphates and anions found in cellular materials (K assoc H 2 PO 4 − >H 2 P 2 O 7 2− >AMP − ≫ADP 2− or ATP 3− over halides, nitrate, or hydrogen sulfate; all Na + salts in water) but also cellular polyphosphate or phospholipids. This selectivity is harnessed in a real-time monitoring of cell lysis by lysozyme, which releases orthophosphate and other phosphates and anions from the cells. This sensitive (LOD 0.4 μM) fluorescence-based microscale method compares favorably with the state-of-the-art techniques but can easily be practiced in a high-throughput screening (HTS) manner. The anion binding and selectivity in aqueous solutions were investigated by NMR and put in context with phosphate binding of the parent calix[4]pyrrole. The microscopic understanding of anion binding by the hybrid receptor was then obtained from a combination of density functional theory (DFT), classical molecular dynamics (MD) with explicit water solvation, and ab initio MD (AIMD) simulations. Correlating the NMR and fluorescence binding data with studies of solvation of the receptor, phosphate anion, and the resulting complex confirms the binding is largely driven by entropic component (TΔS) associated with receptor and anion desolvation.

Anions↗
NASA NTRS1973

Equilibrium properties of a one-dimensional kinetic system.

One-dimensional systems of N = 500 and 250 particles in equilibrium are numerically simulated utilizing the method of molecular dynamics. Periodic boundary conditions are imposed. The classical two-body interaction potential is short range, repulsive and has a corresponding finite force. The equations of state are determined for densities both less and greater than one. Corresponding theoretical isochores are determined from models based on nearest-neighbor interactions and on a truncated virial expansion, and a comparison is made with the experimental isochores. Time independent radial distributions are constructed numerically and discussed. A change of state from a solidlike state to a fluid-gas state based on the penetrability of the particles is predicted. The transition temperatures are estimated from the radial distribution functions and the nearest-neighbor model. Self-diffusion is observed and the corresponding constants are determined from the velocity autocorrelation functions.

Williams, J. H.↗
DOE OSTI2026

Discerning Influences from Enthalpy and Entropy at Aqueous Interfaces Involved in Biomass Conversions in Porous Catalysts

Project Summary: The goal in this proposal is to learn how solvent influences the enthalpies and entropies of catalytic species in zeolite pores. Specifically, the physical, chemical, and structural features of solvent, catalytic species, and zeolite pores that determine enthalpies and entropies of solvation in solution phase biomass conversions will be interrogated using multiscale simulations and machine learning. Multiscale simulations are based off of existing strategies and employ quantum mechanics and classical molecular dynamics, providing an excellent balance between chemical accuracy and computational expense. They are capable of calculating enthalpies and entropies of solvation separately and have been validated in prior work to achieve high accuracy compared to their parent methods.

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

Building and Breaking Carbon Composites with REACTER

Carbon-based composites have become indispensable materials in aerospace and other high-performance applications, yet obtaining a detailed, nanoscale understanding of their morphology and failure mechanisms using only experimental methods remains a difficult challenge. REACTER is a versatile computational modeling tool for atomistic molecular dynamics simulations designed to model chemical reactions at the speed and length scales of classical force fields. In this work, several recent features of REACTER were applied to the creation and subsequent mechanical testing of two classes of carbon composites: carbon nanotube (CNT) composites and carbon fiber reinforced polymers (CFRP). A network of CNTs was grown dynamically using the new ‘create atoms’ feature of REACTER. The CNT filler was embedded into a polyarylacetylene (PAA) matrix by simulated in situ polymerization to obtain the final composite model. To generate the second class of carbon composite, fully carbonized (graphitic) carbon fiber morphologies were created by the method of Desai et al. [1], but using the advanced reaction constraints framework of REACTER. Two fiber models were created, representing a circular carbon fiber core and a flat surface, and similarly infiltrated with resin to obtain the final CFRP structure. Failure mechanisms were elucidated by simulating mechanically induced bond breaking, as characterized by third order DFT-based tight-binding simulations, via a reaction constraint on the total potential energy of the involved atoms.

Polymer↗
NASA NTRSDate not supplied

REACTER: A Heuristic Method for Reactive Molecular Dynamics

REACTER (www.reacter.org) is a heuristic protocol that enables the simulation of complex reactions using atomistic molecular dynamics (MD) with a fixed-valence force field. Incorporating reactions into classical MD with this approach allows modeling of reactive systems over greatly-increased time scales, enabling systems to be modeled with MD that would not otherwise be feasible. One or more competing multi-step reactions or series of reactions can be invoked simultaneously. Special treatment can be applied to neighboring atoms to relax high energy configurations while the simulation progresses. The original version of REACTER, which was implemented in the open-source LAMMPS simulation package as fix bond/react, was only available for serial simulations. This work describes the expansion of fix bond/react for use in parallel simulations, as well as the addition of various new options, including deletion of reaction by-products, reversible reactions, and custom reaction constraints. These new capabilities are demonstrated through large-scale simulations (200,000+ atoms) of the polymerization of polystyrene and nylon 6,6. The morphologies of both polymers are analyzed after reaching >99% extent of polymerization. Finally, the newly-added reversible reactions feature is demonstrated by rupturing these highly-entangled systems under uniaxial strain by defining a chain scission reaction.

polymer simulations, molecular dynamics↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Machine Learning Interatomic Potentials for Modeling Framework Flexibility and Water Uptake in NbOFFIVE-1-Ni Metal–Organic Framework

Metal–organic frameworks (MOFs), with their distinctive porous structures and tunable chemical properties, have shown immense promise in the separation and storage of gases. Currently, the accurate simulation of their adsorptive properties remains challenging, especially for systems where the molecules fit very tightly into the pores. Traditional simulation methods often approximate the frameworks as rigid and do not account for the framework flexibility seen in materials such as NbOFFIVE-1-Ni. First-principles molecular dynamics (FPMD) simulations offer the desired accuracy in modeling this flexibility but are limited by their extensive computational demands, rendering them impractical for long simulations. Conversely, classical force field-based simulations offer computational efficiency but lack the necessary accuracy. Here, to break this accuracy-efficiency trade-off, we have developed machine learning interatomic potentials trained on energies and forces from FPMD to model the framework flexibility of NbOFFIVE-1-Ni in the presence of water over nanosecond time scales. Furthermore, by integrating MLIP-driven molecular dynamics (MLIP-MD) with grand canonical Monte Carlo (GCMC) simulations, we further incorporated framework flexibility into adsorption predictions, yielding water adsorption isotherms that better align with experimental data compared to those of conventional GCMC simulations. These advances offer new opportunities for the design and optimization of MOFs in gas storage and separation applications.

adsorption↗
DOE OSTI2024

Anomalous Interlayer Exciton Diffusion in WS 2 /WSe 2 Moiré Heterostructure

Stacking van der Waals crystals allows for the on-demand creation of a periodic potential landscape to tailor the transport of quasiparticle excitations. We investigate the diffusion of photoexcited electron-hole pairs, or excitons, at the interface of WS 2 /WSe 2 van der Waals heterostructure over a wide range of temperatures. We observe the appearance of distinct interlayer excitons for parallel and antiparallel stacking and track their diffusion through spatially and temporally resolved photoluminescence spectroscopy from 30 to 250 K. While the measured exciton diffusivity decreases with temperature, it surprisingly plateaus below 90 K. Our observations cannot be explained by classical models like hopping in the moiré potential. A combination of ab initio theory and molecular dynamics simulations suggests that low-energy phonons arising from the mismatched lattices of moiré heterostructures, also known as phasons, play a key role in describing and understanding this anomalous behavior of exciton diffusion. Our observations indicate that the moiré potential landscape is dynamic down to very low temperatures and that the phason modes can enable efficient transport of energy in the form of excitons.

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

RASPA3

RASPA3, a molecular simulation code for computing adsorption and diffusion in nanoporous materials and thermodynamic and transport properties of fluids. It implements force field based classical Monte Carlo/molecular dynamics in various ensembles. RASPA3 is rewritten from the ground up in C++23 with speed and code readability in mind. Transition-matrix Monte Carlo is added to compute the density of states and free energies. The Monte Carlo code for rigid molecules is based on quaternions, and the atomic positions needed in the energy evaluation are recreated from the center of mass position and quaternion orientation. The expanded ensemble methodology for fractional molecules, with a scaling parameter λ between 0 and 1, now also keeps track of analytic expressions of dU/dλ, allowing independent verification of the chemical potential using thermodynamic integration. The source code is freely available under the MIT license on GitHub.

Dubbeldam, David↗
NASA NTRS1984

Computer simulation of surface and film processes

All the investigations which were performed employed in one way or another a computer simulation technique based on atomistic level considerations. In general, three types of simulation methods were used for modeling systems with discrete particles that interact via well defined potential functions: molecular dynamics (a general method for solving the classical equations of motion of a model system); Monte Carlo (the use of Markov chain ensemble averaging technique to model equilibrium properties of a system); and molecular statics (provides properties of a system at T = 0 K). The effects of three-body forces on the vibrational frequencies of triatomic cluster were investigated. The multilayer relaxation phenomena for low index planes of an fcc crystal was analyzed also as a function of the three-body interactions. Various surface properties for Si and SiC system were calculated. Results obtained from static simulation calculations for slip formation were presented. The more elaborate molecular dynamics calculations on the propagation of cracks in two-dimensional systems were outlined.

Tiller, W. A.↗
NASA NTRS2001

Better, Cheaper, Faster Molecular Dynamics

Recent, revolutionary progress in genomics and structural, molecular and cellular biology has created new opportunities for molecular-level computer simulations of biological systems by providing vast amounts of data that require interpretation. These opportunities are further enhanced by the increasing availability of massively parallel computers. For many problems, the method of choice is classical molecular dynamics (iterative solving of Newton's equations of motion). It focuses on two main objectives. One is to calculate the relative stability of different states of the system. A typical problem that has' such an objective is computer-aided drug design. Another common objective is to describe evolution of the system towards a low energy (possibly the global minimum energy), "native" state. Perhaps the best example of such a problem is protein folding. Both types of problems share the same difficulty. Often, different states of the system are separated by high energy barriers, which implies that transitions between these states are rare events. This, in turn, can greatly impede exploration of phase space. In some instances this can lead to "quasi non-ergodicity", whereby a part of phase space is inaccessible on time scales of the simulation. To overcome this difficulty and to extend molecular dynamics to "biological" time scales (millisecond or longer) new physical formulations and new algorithmic developments are required. To be efficient they should account for natural limitations of multi-processor computer architecture. I will present work along these lines done in my group. In particular, I will focus on a new approach to calculating the free energies (stability) of different states and to overcoming "the curse of rare events". I will also discuss algorithmic improvements to multiple time step methods and to the treatment of slowly decaying, log-ranged, electrostatic effects.

Pohorille, Andrew↗
DOE OSTI2025

Accelerating charge estimation in molecular dynamics simulations using physics-informed neural networks: corrosion applications

Molecular Dynamics (MD) simulations are used to understand the effects of corrosion on metallic materials in salt brine. Reactive force fields in classical MD enable accurate modeling of bond formation and breakage in the aqueous medium and at the metal-electrolyte interface, while also facilitating dynamic partial charge equilibration. However, MD simulations are computationally intensive and unsuitable for modeling the long time scales characteristic of corrosive phenomena. To address this, we develop reduced-order machine learning models that provide accurate and efficient predictions of charge density in corrosive environments. Specifically, we use Long Short-Term Memory (LSTM) networks to forecast charge density evolution based on atomic environments represented by Smooth Overlap of Atomic Positions (SOAP) descriptors. A physics-informed loss function enforces charge neutrality and electronegativity equivalence. The atomic charges predicted by the deep learning model trained on this work were obtained two orders of magnitude faster than those from molecular dynamics (MD) simulations, with an error of less than 3% compared to the MD-obtained charges, even in extrapolative scenarios, while adhering to physical constraints. This demonstrates the excellent accuracy, computational efficiency, and validity of the developed model. Lastly, even though developed for corrosion, these protocols are formulated in a phenomenon-agnostic manner, allowing application to various variable-charge interatomic potentials and related fields.

Atomistic models↗
DOE OSTI2024

Experimental and Computational Studies of Crystal Nucleation in Composition Gradients (Final Report)

The major goals of this project were to develop and validate a predictive nucleation model that incorporates composition gradients and is applicable to a large range of metallic systems which form both stochiometric and non-stochiometric compounds. Computationally we planned to expand the extent of thermodynamics-based theories and lay the groundwork to improve the general understanding and predictive capabilities of the role that gradients play in phase formation, glass formability and stability, and nucleation and growth events. To address these goals, we used Molecular Dynamics (MD) simulations and both isothermal and isochronal nanocalorimetric experiments on amorphous phases with a controlled composition gradient with the hopes of validating and improving the classical nucleation model and its use in solid solutions. While the computations were successful and identified ways to improve the classical nucleation model, the in situ nanocalorimetry studies proved very challenging due to unexpected difficulties in fabricating effective calorimeters and samples. Thus, we could not experimentally validate our predicted influence of composition gradients on nucleation. Nonetheless, important insights were gained and modifications to the classical nucleation theory were suggested.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

Quantum Hardware-Enabled Molecular Dynamics via Transfer Learning

The ability to perform ab initio molecular dynamics simulations using potential energy surfaces provided by quantum computers would open the door to virtually exact dynamics for a variety of chemical and biochemical systems, with impacts on catalysis and biophysics. Nonetheless, performing molecular dynamics on surfaces produced by quantum hardware has been hampered by the noisy energies typically produced by quantum computers and challenges associated with computing gradients and scaling to large systems interest. A recent set of advances in machine learning, known as transfer learning, provides a new path forward for molecular dynamics simulations on quantum hardware. Transfer learning offers a workaround, where one first trains models on larger, less accurate classical datasets and then refines them on smaller, more accurate quantum datasets. We explore this approach by training machine learning models to predict a molecule's potential energy based on its geometric structure using Behler-Parrinello neural networks. When successfully trained, the model enables energy gradient predictions necessary for dynamic simulations. To reduce the quantum resources needed, the model is initially trained with data derived from classical density functional theory and subsequently refined with a smaller dataset obtained from a variational quantum eigensolver optimization of the unitary coupled cluster ansatz. We show that this approach significantly reduces the size of the needed quantum training dataset while capturing the high accuracies needed within quantum chemistry simulations. The success of this two-step training method opens more opportunities to apply machine learning models on quantum data, a significant stride towards efficient quantum-classical hybrid computational models.

quantum computing↗
DOE OSTI2024

Effects of Nonequilibrium Atomic Structure on Ionic Diffusivity in LLZO: A Classical and Machine Learning Molecular Dynamics Study

To improve the performance of electrochemical devices, it is essential to understand the effects of nonequilibrium motifs in solids, such as grain boundaries, amorphous phases, and highly strained regions, on atomic-scale transport and stability. Molecular dynamics simulations are used to explore the combined effect of far-from-equilibrium atomic structures and the choice of interatomic potential on ionic diffusivity predictions for Li 7 La 3 Zr 2 O 12 (LLZO), a promising solid electrolyte for all-solid-state batteries. Amorphization and high strain are considered using both classical Buckingham interatomic potentials and machine learning force fields. Here we find that both crystalline expansion and amorphization tend to slow diffusion, although the different physical encodings in the two potentials impact the properties in different ways. We trace these variations to a combination of structural and transport factors, the contributions of which are deconvoluted computationally. Graph-based analysis reveals that the variations for amorphous LLZO arise from the connectivity of diffusion pathways within the predicted structures, which generally correlates with diffusivity and is notably higher for structures generated by the machine learning force fields. Our study provides additional insight into the relationship between atomic structure and diffusivity in LLZO, while also highlighting the need for care in choosing and validating potentials to simulate far from equilibrium structures.

25 ENERGY STORAGE↗