Search NASA⌕ Search

SEARCH · Search NASA

Results for “Molecular dynamics computer 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.

At least 127 records · Page 7
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Molecular growth upon ionization of van der Waals clusters containing HCCH and HCN is a pathway to prebiotic molecules

The growth mechanisms of organic molecules in an ionizing environment such as the interstellar medium are not completely understood. Here we examine by means of ab initio molecular dynamics (AIMD) simulations and density functional theory (DFT) computations the possibility of bond formation and molecular growth upon ionization of van der Waals clusters of pure HCN clusters, and mixed clusters of HCN and HCCH, both of which are widespread in the interstellar medium. Ionization of van der Waals clusters can potentially lead to growth in low temperature and low-density environments. Our results show, that upon ionization of the pure HCN clusters, strongly bound stable structures are formed that contain N∙∙∙H bonds, and growth beyond pairwise HCN molecules is seen only in a small percentage of cases. In contrast, mixed clusters, where HCCH is preferentially ionized over HCN, can grow up to 3 or 4 units long with new carbon–carbon and carbon–nitrogen covalent bonds. Moreover, cyclic molecules formed, such as the radical cation of pyridine, which is a prebiotic molecule. The results presented here are significant as they provide a feasible pathway for molecular growth of small organic molecules containing both carbon and nitrogen in cold and relatively denser environments such as in dense molecular clouds but closer to the photo-dissociation regions, and protoplanetary disks. In the mechanism we propose, first, a neutral van der Waals cluster is formed. Once the cluster is formed it can undergo photoionization which leads to chemical reactivity without any reaction barrier.

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

Reinforcement learning based hybrid bond-order coarse-grained interatomic potentials for exploring mesoscale aggregation in liquid–liquid mixtures

Exploring mesoscopic physical phenomena has always been a challenge for brute-force all-atom molecular dynamics simulations. Although recent advances in computing hardware have improved the accessible length scales, reaching mesoscopic timescales is still a significant bottleneck. Coarse-graining of all-atom models allows robust investigation of mesoscale physics with a reduced spatial and temporal resolution but preserves desired structural features of molecules, unlike continuum-based methods. Here, we present a hybrid bond-order coarse-grained forcefield (HyCG) for modeling mesoscale aggregation phenomena in liquid–liquid mixtures. The intuitive hybrid functional form of the potential offers interpretability to our model, unlike many machine learning based interatomic potentials. We parameterize the potential with the continuous action Monte Carlo Tree Search (cMCTS) algorithm, a reinforcement learning (RL) based global optimizing scheme, using training data from all-atom simulations. The resulting RL-HyCG correctly describes mesoscale critical fluctuations in binary liquid–liquid extraction systems. cMCTS, the RL algorithm, accurately captures the mean behavior of various geometrical properties of the molecule of interest, which were excluded from the training set. The developed potential model along with the RL-based training workflow could be applied to explore a variety of other mesoscale physical phenomena that are typically inaccessible to all-atom molecular dynamics simulations.

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

Molecular Dynamics Investigation of the Structural and Mechanical Properties of Off-Stoichiometric Epoxy Resins

Molecular dynamic (MD) simulations were performed to compute the mechanical properties of off-stoichiometric epoxy resins as a function of hardener/epoxy mixture ratio (r). Properties were characterized in relation to their microscopic structures. Such resins have been used recently for adhesive-free bonding of large-scale composite structures using the co-curing-ply method. In this process, two partially precured composite panels with hardener-poor (HP) off-stoichiometric resins are coupled with ply(ies) of complementary hardener-rich (HR) formulations and then cured simultaneously. This bonding process has the potential to produce reliable and certifiable composite joints without the need for additional fasteners, which are often required for many conventional bonding methods because even small amounts of contamination can cause a weak bond. The reflow and mixing of the HP/HR resin in this bonding process result in a joint with no discernable interface that should not be susceptible to surface contamination. However, incomplete mixing of the two offset resins may result in chemical heterogeneity of the cured polymeric joint. Thus, different r values may be obtained across the joint. Classical MD simulations were performed to compute the Young’s modulus of polymers with different r values and correlate their properties to network structures. High stiffness was associated with molecular packing due to chemical crosslinking, leading to a single network structure. Moreover, the networks became denser as the ratio approached the stoichiometric value r = 1. Thus, the r = 1 systems were single clusters, with high stiffness, high molecular weight, and a high degree of crosslinking. Structural properties such as radius of gyration and mean square displacement were determined to investigate the variation in the stiffness with respect to r. This MD simulation study was validated with experimental measurements.

Chang Woon Jang↗
NASA NTRS1981

Computer simulation of surface and film processes

A molecular dynamics technique based upon Lennard-Jones type pair interactions is used to investigate time-dependent as well as equilibrium properties. The case study deals with systems containing Si and O atoms. In this case a more involved potential energy function (PEF) is employed and the system is simulated via a Monte-Carlo procedure. This furnishes the equilibrium properties of the system at its interfaces and surfaces as well as in the bulk.

Tiller, W. A.↗
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 OSTI2025

Breaking the mold: Overcoming the time constraints of molecular dynamics on general-purpose hardware

The evolution of molecular dynamics (MD) simulations has been intimately linked to that of computing hardware. For decades following the creation of MD, simulations have improved with computing power along the three principal dimensions of accuracy, atom count (spatial scale), and duration (temporal scale). Since the mid-2000s, computer platforms have, however, failed to provide strong scaling for MD, as scale-out central processing unit (CPU) and graphics processing unit (GPU) platforms that provide substantial increases to spatial scale do not lead to proportional increases in temporal scale. Important scientific problems therefore remained inaccessible to direct simulation, prompting the development of increasingly sophisticated algorithms that present significant complexity, accuracy, and efficiency challenges. While bespoke MD-only hardware solutions have provided a path to longer timescales for specific physical systems, their impact on the broader community has been mitigated by their limited adaptability to new methods and potentials. In this work, we show that a novel computing architecture, the Cerebras wafer scale engine, completely alters the scaling path by delivering unprecedentedly high simulation rates up to 1.144 M steps/s for 200 000 atoms whose interactions are described by an embedded atom method potential. This enables direct simulations of the evolution of materials using general-purpose programmable hardware over millisecond timescales, dramatically increasing the space of direct MD simulations that can be carried out. In this paper, we provide an overview of advances in MD over the last 60 years and present our recent result in the context of historical MD performance trends.

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

Accelerated Modeling of Lithium Diffusion in Solid State Electrolytes using Artificial Neural Networks

Abstract Previous efforts to understand structure‐function relationships in high ionic conductivity materials for solid state batteries have predominantly relied on density functional theory (DFT‐) based ab initio molecular dynamics (MD). Such simulations, however, are computationally demanding and cannot be reasonably applied to large systems containing more than a hundred atoms. Here, an artificial neural network (ANN) is trained to accelerate the calculation of high accuracy atomic forces and energies used during such MD simulations. After carefully training a robust ANN for four and five element systems, nearly identical lithium ion diffusivities are obtained for Li 10 GeP 2 S 12 (LGPS) when benchmarking the ANN‐MD results with DFT‐MD. Applying the ANN‐MD approach, the effect of chlorine doping on the lithium diffusivity is calculated in an LGPS‐like structure and it is found that a dopant concentration of 1.3% maximizes ionic conductivity. The optimal concentration balances the competing consequences of effective atomic radii and dielectric constants on lithium diffusion and agrees with the experimental composition. Performing simulations at the resolution necessary to model experimentally relevant and optimal concentrations would be infeasible with traditional DFT‐MD. Systems that require a large number of simulated atoms can be studied more efficiently while maintaining high accuracy with the proposed ANN‐MD framework.

Rao, Karun K.↗
DOE OSTI2023

Multifunctional Chiral Chemically‐Powered Micropropellers for Cargo Transport and Manipulation

Practical applications of synthetic self-propelled nano and microparticles for microrobotics, targeted drug delivery, and manipulation at the nanoscale are rapidly expanding. However, fabrication limitations often hinder progress, resulting in relatively simple shapes and limited functionality. Here, taking advantage of 3D nanoscale printing, chiral micropropellers powered by the hydrogen peroxide reduction reaction are fabricated. Due to their chirality, the propellers exhibit multifunctional behavior controlled by an applied magnetic field: spinning in place (loitering), directed migration in the prescribed direction, capture, and transport of polymer cargo particles. Design parameters of the propellers are optimized by computation modeling based on mesoscale molecular dynamics. It is predicted by computer simulations, and confirmed experimentally, that clockwise rotating propellers attract each other and counterclockwise repel. These results shed light on how chirality and shape optimization enhance the functionality of synthetic autonomous micromachines.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Critical mineral substitutions in IN617: A combined computational and experimental approach to performance evaluation and feasibility

Addressing the escalating demand for critical minerals (CMs) driven by global climate change initiatives, this study explores compositional modifications to Inconel 617 (IN617) by substituting cobalt (Co) with manganese (Mn) across various atomic percentages. Here, we conducted a computational feasibility study employing Molecular Dynamics (MD) simulations to provide strategic guidance for experimental validation. The simulations analyzed tensile strength and corrosion resistance for five modified compositions (M1 to M5) to identify optimal properties. Tensile tests on cubic simulation cells were performed to generate stress-strain curves, revealing the impact of Co replacement with Mn on tensile strength—a metric correlated with hardness. Oxygen penetration simulations were conducted to evaluate corrosion resistance, indicating that reduced oxygen penetration depth corresponds to enhanced resistance. Promising compositions underwent phase diagram calculations for assessing phase stability. The optimal composition (M1), characterized by high tensile strength and minimal oxygen penetration, was chosen for experimental validation using induction melting and friction stir consolidation techniques. The materials were further characterized using SEM-EDS, XRD, and Vickers hardness testing. Our findings suggest that Mn substitution in IN617 can yield mechanical performance at par with high-Co alloys in energy-critical applications.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Effect of Hydrophilicity and Interfacial Water Structure on Particle Attachment

Oriented attachment (OA) of nanoparticles is an important and common crystallization pathway. Efforts to isolate and quantify the interparticle forces that underlie OA have only recently been initiated and the role played by the interfacial solvent structure in the confined region between nanocrystal surfaces in enabling and controlling OA remains an open question. In this work, classical molecular dynamics (MD) simulations were performed to compute the free energy landscape of nanoparticles approaching mineral substrates in both vacuum and pure water. Periclase (MgO), muscovite (KAl 3 Si 3 O 10 (OH) 2 ), and pyrophyllite (Al 2 Si 4 O 10 (OH) 2 ) nanoparticles exposing the (100), (001), and (001) surfaces, respectively, were brought to contact with substrates exposing the same surfaces. These mineral surfaces were selected to investigate a range of hydrophilicities (periclase > muscovite > pyrophyllite). The MD simulations revealed that the interfacial water structure created free energy minima for nanoparticle approach that corresponded to integer numbers of intervening water layers. The depth of the free energy minima and the free energy barriers between minima correlated positively with the hydrophilicity of the mineral surface. While the range of solvent-mediated interactions was predicted to be short (a few intervening water layers at most for MgO), the free energy barriers were calculated to be large (fraction of an eV per nm 2 ) and the effect of misalignment diminished rapidly with the number of intervening water layers. The results help explain the formation of mesocrystals and the observed ability of nanoparticles to rotate before attachment while in a long-lived solvent-separated state.

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

Lessons from combined experimental and theoretical examination of the FTIR and 2D-IR spectroelectrochemistry of the amide I region of cytochrome c</em

Amide I difference spectroscopy is widely used to investigate protein function and structure changes. In this article, we show that the common approach of assigning features in amide I difference signals to distinct secondary structure elements in many cases may not be justified. Evidence comes from Fourier transform infrared (FTIR) and 2D-IR spectroelectrochemistry of the protein cytochrome c in the amide I range, in combination with computational spectroscopy based on molecular dynamics (MD) simulations. This combination reveals that each secondary structure unit, such as an alpha-helix or a beta-sheet, exhibits broad overlapping contributions, usually spanning a large part of the amide I region, which in the case of difference absorption experiments (such as in FTIR spectroelectrochemistry) may lead to intensity-compensating and even sign-changing contributions. We use cytochrome c as the test case, as this small electron-transferring redox-active protein contains different kinds of secondary structure units. Upon switching its redox-state, the protein exhibits a different charge distribution while largely retaining its structural scaffold. Our theoretical analysis suggests that the change in charge distribution contributes to the spectral changes and that structural changes are small. However, in order to confidently interpret FTIR amide I difference signals in cytochrome c and proteins in general, MD simulations in combination with additional experimental approaches such as isotope labeling, the insertion of infrared labels to selectively probe local structural elements will be required. In case these data are not available, a critical assessment of previous interpretations of protein amide I 1D- and 2D-IR difference spectroscopy data is warranted.

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

Intrinsic bulk viscosity of the one-component plasma

The intrinsic bulk viscosity of the one-component plasma (OCP) is computed and analyzed using equilibrium molecular dynamics simulations and the Green-Kubo formalism. It is found that bulk viscosity exhibits a maximum at Γ ≈ 1, corresponding to the condition that the average kinetic energy of particles equals the potential energy at the average interparticle spacing. The weakly coupled and strongly coupled limits are analyzed and used to construct a model that captures the full range of coupling strengths simulated: Γ ≈ 10 –2 –10 2 . Simulations are also run of the Yukawa one-component plasma (YOCP) in order to understand the impact of electron screening. It is found that electron screening leads to a smaller bulk viscosity due to a reduction in the excess heat capacity of the system. Bulk viscosity is shown to be at least an order of magnitude smaller than shear viscosity in both the OCP and YOCP. The generalized frequency-dependent bulk viscosity coefficient is also analyzed. This is found to exhibit a peak near twice the plasma frequency in strongly coupled conditions, which is associated with the oscillatory decay observed in the bulk viscosity autocorrelation function. As a result, the generalized shear and bulk viscosity coefficients are found to have a similar magnitude for ω ≳ 2⁢ω p at strongly coupled conditions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Stress Intensity Thresholds for Development of Reliable Brittle Materials

Brittle material failure in high consequence systems can appear random and unpredictable at subcritical stresses. Gaps in our understanding of how structural flaws and environmental factors (humidity, temperature) impact fracture propagation need to be addressed to circumvent this issue. A combined experimental and computational approach composed of molecular dynamics (MD) simulations, numerical modeling, and atomic force microscopy (AFM) has been undertaken to identify mechanisms of slow crack growth in silicate glasses. AFM characterization of crack growth as slow as 10 -13 m/s was observed, with some stepwise crack growth. MD simulations have identified the critical role of inelastic relaxation in crack propagation, including evolution of the structure during relaxation. A numerical model for the existence of a stress intensity threshold, a stress intensity below which a fracture will not propagate, was developed. This transferrable model for predicting slow crack growth is being incorporated into mission-based programs.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Enabling Efficient Water Splitting with Advanced Materials Designed for High pH Membrane Interface

This project was focused on developing the durable, high-performance materials and interfaces for advanced water splitting, enabling a clear pathway for achieving <$2/Kg H2 (on scale) with efficiency of 43 kWh/kg H 2 via anion exchange membrane (AEM)-based electrolysis. We aimed to advance these final goals via an improved fundamental understanding of both hydrogen and oxygen evolution reactions (HER/OER) leading to novel platinum group metal (PGM)-free catalyst materials in conjunction with critical improvements in membrane and ionomers and gas evolution electrodes with corresponding characterization and testing. Northeastern University (NU) lead this effort focusing on catalyst development and characterization (both in situ and ex situ) while project partners lead improvements in ionomer and membrane materials and will aid in the development of specialized electrode and membrane electrode assemblies. In addition, close collaboration occured with the HydroGEN Energy Materials Network (EMN) National Laboratory consortium including efforts related to use of advanced ionomers, durability protocols and validation of electrolyzer materials (e.g. NREL), multiscale modeling and computation (e.g. LBNL), and molecular dynamics (MD) simulations of the membrane catalyst interface (e.g. SNL). The interactions with HydroGEN included exchange of data and materials as needed to facilitate project success.

08 HYDROGEN↗
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↗
DOE OSTI2021

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.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗