Search NASA⌕ Search

SEARCH · Search NASA

Results for “molecular dynamics (MD)”

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 271 records · Page 15
All sourcesNASA NTRSDOE OSTI
DOE OSTI2021

Multiscale Concurrent Atomistic-Continuum (CAC) modeling of multicomponent alloys

We report strengthening in complex multicomponent systems such as solid solution alloys is controlled primarily by the dynamic interactions between dislocation lines and heterogeneously distributed solute species. Modeling of extended defect length scales in such multicomponent systems becomes prohibitively expensive, motivating the development of reduced order approaches. This work explores the application of the Concurrent Atomistic-Continuum (CAC) method to model dislocation mobility in random alloys at extended length scales. By employing recently developed average-atom interatomic potentials, the average “bulk” material response in coarse-grained regions interacts with true random solute species in the atomistic-scale domain. We demonstrate that spurious stresses in domain resolution transition regions are eliminated entirely due to the CAC formulation. Simultaneously, the key details of local stress fluctuation due to randomness in the dislocation core region are captured, and fluctuating stress smoothly decays to the long-range dislocation stress field response. Dislocation mobility calculations, for line lengths over 400 nm, are computed as a function of alloy composition in the model FeNiCr system and compared to full molecular dynamics (MD). The results capture the composition-dependent trends, while reducing degrees of freedom by nearly 40%. This approach can be readily extended to any system described by an EAM potential and facilitates the study of large-scale defect dynamics in complex solute environments to support computational alloy design.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Extreme Metastability of Diamond and its Transformation to the BC8 Post-Diamond Phase of Carbon

Diamond possesses exceptional physical properties due to its remarkably strong carbon–carbon bonding, leading to significant resilience to structural transformations at very high pressures and temperatures. Despite several experimental attempts, synthesis and recovery of the theoretically predicted post-diamond BC8 phase remains elusive. Through quantum-accurate multimillion atom molecular dynamics (MD) simulations, we have uncovered the extreme metastability of diamond at very high pressures, significantly exceeding its range of thermodynamic stability. We predict the post-diamond BC8 phase to be experimentally accessible only within a narrow high pressure–temperature region of the carbon phase diagram. The diamond to BC8 transformation proceeds through premelting followed by BC8 nucleation and growth in the metastable carbon liquid. Here we propose a double-shock compression pathway for BC8 synthesis, which is currently being explored in experiments at the National Ignition Facility.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2021

Understanding the plasticity contributions during laser-shock loading and spall failure of Cu microstructures at the atomic scales

A hybrid atomic-scale and continuum modeling framework is used to study the microstructural evolution during the laser-induced shock deformation and failure (spallation) of copper microstructures. A continuum two-temperature model (TTM) is used to account for the interaction of Cu atoms with a laser in molecular dynamics (MD) simulations. The MD-TTM simulations study the effect of laser loading conditions (laser fluence) on the microstructure (defects) evolution during various stages of shock wave propagation, reflection, and interaction in single-crystal (sc) Cu systems. In addition, the role of the microstructure is investigated by comparing the defect evolution and spall response of sc-Cu and nanocrystalline (nc) Cu systems. The defect (stacking faults and twin faults) evolution behavior in the metal at various times is further characterized using virtual in situ selected area electron diffraction and x-ray diffraction during various stages of evolution of microstructure. Here, the simulations elucidate the uncertain relation between spall strength and strain-rate and the much stronger relation between the spall strength and the temperatures generated due to laser shock loading for the small Cu sample dimensions considered here.

36 MATERIALS SCIENCE↗
DOE OSTI2023

Molecular-dynamics simulation methods for macromolecular crystallography

It is investigated whether molecular-dynamics (MD) simulations can be used to enhance macromolecular crystallography (MX) studies. Historically, protein crystal structures have been described using a single set of atomic coordinates. Because conformational variation is important for protein function, researchers now often build models that contain multiple structures. Methods for building such models can fail, however, in regions where the crystallographic density is difficult to interpret, for example at the protein–solvent interface. To address this limitation, a set of MD–MX methods that combine MD simulations of protein crystals with conventional modeling and refinement tools have been developed. In an application to a cyclic adenosine monophosphate-dependent protein kinase at room temperature, the procedure improved the interpretation of ambiguous density, yielding an alternative water model and a revised protein model including multiple conformations. The revised model provides mechanistic insights into the catalytic and regulatory interactions of the enzyme. The same methods may be used in other MX studies to seek mechanistic insights.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2023

Unexpected Kinetic Solvent Effects Enhance Activity and Selectivity in Biphasic Systems

Biphasic dehydration of fructose to 5-hydroxymethylfurfural (HMF) has shown unprecedented increases in productivity, but a mechanistic understanding is lacking. Herein, we couple fast experimental reaction kinetics, multiscale modeling (phase behavior, classical molecular dynamics(MD), and quantum mechanics/molecular mechanics MD), in situ sampling, and IR and 13 C-NMR spectroscopy to elucidate the complex effects of nonpolar extracting organic solvents on the kinetics of fructose dehydration. We show that these organic solvents can reach significant mutual solubility with water at reaction temperatures, enabling the partition of the sugar and catalyst into the extracting phase. In the organic-rich environment, the dehydration of fructose proceeds faster and more selectively than in water due to increased relative abundance of the reactive furanose isomer, enhanced water–catalyst–substrate interactions driven by nanophase separation, and higher product stability stemming from preferential solvation. Furthermore, we demonstrate that these solvent effects impact other critical biphasic reactions in biomass upgrading and provide qualitative principles for solvent selection.

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

Interfacial water asymmetry at ideal electrochemical interfaces

Controlling electrochemical reactivity requires a detailed understanding of the charging behavior and thermodynamics of the electrochemical interface. Experiments can independently probe the overall charge response of the electrochemical double layer by capacitance measurements and the thermodynamics of the inner layer with potential of maximum entropy measurements. Relating these properties by computational modeling of the electrochemical interface has so far been challenging due to the low accuracy of classical molecular dynamics (MD) for capacitance and the limited time and length scales of ab initio MD. Here, we combine large ensembles of long-time-scale classical MD simulations with charge response from electronic density functional theory to predict the potential-dependent capacitance of a family of ideal aqueous electrochemical interfaces with different peak capacitances. We show that while the potential of maximum capacitance varies, this entire family exhibits an electrode charge of maximum capacitance (CMC) between -2.9 and -2.2 μC/cm 2 , regardless of the details in the electronic response. Simulated heating of the same interfaces reveals that the entropy peaks at a charge of maximum entropy (CME) of -5.1 ± 0.6 μC/cm 2 , in agreement with experimental findings for metallic electrodes. The CME and CMC both indicate asymmetric response of interfacial water that is stronger for negatively charged electrodes, while the difference between CME and CMC illustrates the richness in behavior of even the ideal electrochemical interface.

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

Effect of free surfaces on localized plastic deformation in single-crystal nickel containing helium bubbles and radiation-induced self-interstitial atom clusters

This paper presents molecular dynamics (MD) simulations of plastic deformation in single-crystal Ni models containing internal distributions of helium (He) bubbles and self-interstitial atom (SIA) clusters corresponding to experimental characterization of ex-service Inconel X-750 components. Here, we simulate single crystals with and without free surfaces at room temperature and 588K, the latter temperature corresponding to the Inconel X-750 service conditions. Both models are deformed under multiple slip to a maximum deviatoric strain of ~10 %. A dislocation slip band forms in the model without free surfaces, but not in the model with free surfaces. However, the latter exhibits nano-twin nucleation from surfaces. Helium (He) bubbles that overlap with a slip band or a nano-twinned region deform more severely than other bubbles in the model. Nevertheless, the bubbles do not link up and do not form configurations that resemble incipient cracks.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Interface and Thermophysical Properties of R 32 Refrigerant

Driven by the urgent demand for efficient cooling in microelectronics and advanced thermal management systems, difluoromethane (R32/CH 2 F 2 ) has emerged as a promising candidate owing to its favorable thermophysical properties, including high heat transfer efficiency and low viscosity. While bulk properties such as density, viscosity, and thermal conductivity have been widely studied, interfacial properties, including surface tension and interfacial thickness, remain comparatively underexplored, despite their importance in phase-transition dynamics. Here, we perform molecular dynamics (MD) simulations from 180 to 300 K using an optimized transferable force field for fluoropropenes with enhanced electrostatics to assess both bulk and interfacial behavior of R32. Simulations reproduced density within ±2.1%, viscosity within 3.05%, and thermal conductivity within 7.41% of NIST reference data. Heat capacities (C p and C v ) were predicted within 5%. For interfacial properties, surface tension trends were reproduced within 13.58% deviation, and the vapor–liquid coexistence curve closely matched reference data, yielding a critical temperature of 345.7 K (1.6% deviation) and a critical density of 0.397 g/cm 3 (6.4% deviation). Importantly, the vapor–liquid interface exhibited pronounced temperature-dependent broadening across the 180–290 K range. This behavior correlates with increasing molecular kinetic energy, reduction in intermolecular cohesive interactions, and a progressive loss of preferential dipole alignment, which collectively enhance thermal fluctuation amplitudes at elevated temperatures. These validated results provide predictive molecular-level insights, particularly for interfacial properties that remain less characterized. By reducing property prediction errors in key parameters such as critical temperature, this work provides reliable inputs for heat-exchanger and system models. Such correlations can support optimized component sizing, improved performance, and reduced refrigerant charge. Beyond R32, the methodology offers a transferable framework for blended and next-generation low-GWP refrigerants, contributing to sustainable thermal management aligned with the 2027 EU F-Gas regulation and 2030 Kigali Amendment.

Fluids↗
DOE OSTI2025

Binding Free Energy Analysis of Colicin D, E3 and E8 to Their Respective Cognate Immunity Proteins Using Computational Simulations

Colicins are antimicrobial proteins produced by bacteria for the purpose of destroying neighboring bacteria. Colicin activity is neutralized by a specific cognate immunity protein in order to protect the host. This study investigates the structural and binding mechanisms underlying the interaction of colicin-D, -E3 and -E8 to their respective immunity proteins (ImD, Im3 and Im8) using structure prediction, molecular dynamics (MD) simulations and MM-PBSA approach of free energy calculations. High-confidence colicin-immunity (Col-Im) complex structures predicted using AlphaFold2 were subjected to MD simulations of 150 ns with GROMACS and were analyzed for the binding free energy calculation using gmx_MMPBSA. Results showed that the complex of Col_E3-Im3 exhibited the most favorable binding free energy, driven by strong van der Waals and electrostatic interactions. Col_D-ImD and Col_E8-Im8 also showed the favorable binding. Electrostatics and hydrogen bonding emerged as a key factor driving binding and stability, while polar solvation acted as a destabilizing factor across all systems. These outcomes provide an understanding of the molecular mechanisms of Col-Im systems, with potential applications for developing natural antimicrobials for food safety.

Biochemistry & Molecular Biology↗
DOE OSTI2025

Nonequilibrium effects in high-gain inertial confinement fusion

Recent experimental demonstrations of ignition and target gain in inertial confinement fusion (ICF) have stimulated interest in exploring the fundamental physics of violent deuterium-tritium (DT) burn in high-gain ICF targets. A significant DT-burn fraction is a necessary condition for high energy gain and large neutron yields (>100MJ). Using classical molecular-dynamics (MD) simulations and a hybrid fluid-kinetic model, we examine how a large fraction of low-energy 𝛼 particles can kick D and T ions out of equilibrium in high-gain ICF targets. The MD results suggest that (1) temperatures of 𝑇 𝐷 and 𝑇 𝑇 can differ by as much as ∼20% of their mean temperature and (2) the deviation of the DT energy distribution from the Maxwell-Boltzmann function can exceed ∼30%. Some of these MD observations, such as the preferential heating of D ions by low-energy 𝛼 particles and the temperature separation, can be explained by a proposed hybrid fluid-kinetic model. Furthermore, the implication of such nonequilibrium effects on the DT reactivity is also discussed.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2022

Mechanistic study of chemical looping reactions between solid carbon fuels and CuO

Copper (Cu) based chemical-looping combustion (CLC) is a promising process that utilizes solid carbon fuel such as coal and biomass. Understanding the reaction kinetics in this process can facilitate the industrial design of CLC units. Here, in this work, molecular dynamics (MD) simulations were performed to investigate the reaction kinetics of n-butane and two simplified solid carbon fuels (lignite and anthracite) with/without copper oxide (CuO) nanoparticle. In addition, experiments were conducted on the thermal characteristics, flammability, and flame speeds of CuO and solid carbon mixtures. For n-butane, oxidation on the CuO surface is the primary reaction since the activation energy of the surface reaction is much lower than that of oxygen (O 2 ) combustion (9.2 vs. 53.3 kcal/mol). On the other hand, for lignite, there is a smaller difference in the activation energy of O 2 combustion and surface reaction (23.04 vs. 6.34 kcal/mol). We hypothesize that some solid carbon fuels can have different reaction kinetics dependent upon temperature. This is proven by the case of a simplified anthracite coal: an Arrhenius plot shows that this solid carbon fuel has two different reaction kinetic regimes and the critical temperature for the change in kinetics is related to the oxygen uncoupling of the CuO nanoparticle. Like the modeling simulations, a change in activation energy is observed in the experimental results, where desorption of molecular oxygen from CuO becomes important.

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

From molecular to macroscopic: predicting liquid–liquid phase equilibria and small-angle scattering of mixtures of organic liquids from atomistic simulation using Kirkwood–Buff theory

Macroscopic phase equilibria between solutions define the functionality of many biological and industrial processes, yet they are challenging to predict due to the inherent complexity of liquids containing large molecules. This work introduces an approach for the purely predictive calculation of such phase equilibria in temperature-composition space from molecular dynamics (MD) simulations at one temperature in the single-phase region. We use an approach developed previously to obtain the entropic and enthalpic contributions to the free energy of mixing from the atomic-scale information given by MD simulations via Kirkwood–Buff theory. This allows us to accurately estimate the free energy of mixing as a function of temperature, and thus obtain liquid–liquid phase equilibria, including liquid–liquid critical points, associated binodal and spinodal lines, and composition fluctuations across a region of temperature and composition. Results for binary malonamide–alkane systems are validated by comparison to a direct experimental probe of the fluctuations: the small angle X-ray scattering intensity near zero wavenumber. The MDKB → Phase method demonstrated here provides a significant improvement in predicting liquid–liquid equilibria and free energy as a function of temperature for our systems of interest compared to conventional thermodynamic models. The accurate performance of this purely predictive approach lies in its preservation of atomistic details when determining thermodynamic properties. Furthermore, its inherent extensibility to multi-component systems will likely make the MDKB → Phase approach a valuable general tool for connecting molecular interactions to macroscopic phase equilibria and for the computational screening of materials for targeted thermodynamic behavior.

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

Understanding the structural mechanics of ligated DNA crystals via molecular dynamics simulation

DNA self-assembly is a highly programmable method to construct arbitrary architectures based on sequence complementarity. Among various constructs, DNA crystals are macroscopic crystalline materials formed by assembling motifs via sticky end association. Due to their high structural integrity and size ranging from tens to hundreds of micrometers, DNA crystals offer unique opportunities to study the structural properties and deformation behaviors of DNA assemblies. For example, enzymatic ligation of sticky ends can selectively seal nicks resulting in more robust structures with enhanced mechanical properties. However, the research efforts have been mostly on experiments involving different motif designs, structural optimization, or new synthesis methods, while their mechanics are not yet fully understood. The complex properties of DNA crystals are difficult to study via experiments alone, and numerical simulation can complement and aid the experiments. The coarse-grained molecular dynamics (MD) simulation is a powerful tool that can probe the mechanics of DNA assemblies. Here, we investigate DNA crystals made of four different motif lengths with various ligation patterns (full ligation, major directions, connectors, and in-plane) using oxDNA, an open-source, coarse-grained MD platform. We found that several distinct deformation stages emerge in response to mechanical loading and that the number and the location of ligated nucleotides can significantly modulate structural behaviors. These findings should be useful for predicting crystal properties and thus improving the design.

DNA crystal↗
DOE OSTI2020

Permeability of the HIV-1 capsid to metabolites modulates viral DNA synthesis

Reverse transcription, an essential event in the HIV-1 life cycle, requires deoxynucleotide triphosphates (dNTPs) to fuel DNA synthesis, thus requiring penetration of dNTPs into the viral capsid. The central cavity of the capsid protein (CA) hexamer reveals itself as a plausible channel that allows the passage of dNTPs into assembled capsids. Nevertheless, the molecular mechanism of nucleotide import into the capsid remains unknown. Employing all-atom molecular dynamics (MD) simulations, we established that cooperative binding between nucleotides inside a CA hexamer cavity results in energetically favorable conditions for passive translocation of dNTPs into the HIV-1 capsid. Furthermore, binding of the host cell metabolite inositol hexakisphosphate (IP 6 ) enhances dNTP import, while binding of synthesized molecules like benzenehexacarboxylic acid (BHC) inhibits it. The enhancing effect on reverse transcription by IP 6 and the consequences of interactions between CA and nucleotides were corroborated using atomic force microscopy, transmission electron microscopy, and virological assays. Collectively, our results provide an atomistic description of the permeability of the HIV-1 capsid to small molecules and reveal a novel mechanism for the involvement of metabolites in HIV-1 capsid stabilization, nucleotide import, and reverse transcription.

60 APPLIED LIFE SCIENCES↗
DOE OSTI2022

The role of deviatoric stress and dislocations on the α to ω phase transformation in Ti

Under extreme conditions, α-Ti becomes unstable and transforms either into β-Ti at high temperature or into ω-Ti at high pressure. In what concerns the α to ω phase transformation (PT), there has been a wide range of experimentally reported transition pressures from approximately 2 to 15 GPa at room temperature. Deviatoric stresses and internal defects are often assumed to be the root cause of this variation. Here, in this study, these postulates are revisited using both continuum mechanics and molecular dynamics (MD) simulations. First, a simple continuum model, assuming linear elasticity and isotropic plasticity, is developed to describe the effects of applied stress and dislocations on the stability of an ω nucleus in an infinite α domain. Second, a new MD simulation method is developed to generate an ω nucleus in the α domain utilizing the displacement field identified from the topological analysis. Results from MD simulations show that despite the fact that phase diagrams typically delineate the limits between two phases in terms of only P and T, deviatoric stress promotes the α to ω phase transformation by reducing the critical radius above which an ω nucleus is stable. Furthermore, the required deviatoric stress to nucleate and stabilize a nanoscale ω nucleus is likely emanating from the internal stress of defects such as dislocations. The MD-informed micromechanics models are used to identify favorable configurations where dislocations help favor the α to ω transformation. These configurations show that the interaction with a basal or prismatic dislocation reduces the critical radius of a ω nucleus by about 10 or 16 %, respectively. In addition, prismatic edge dislocations are found to promote the growth of ω nucleus when interacting with the ($\bar{1}$$\bar{1}$20) α //(0001) ω interface. Importantly, a simple model of the arrival of dislocation at an ω nucleus suggests that PT does not necessarily require a pile-up to be present but could alternatively be mediated by a constant rapid flow of dislocations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Role of the Mineral in the Self-Healing of Cracks in Human Enamel

Human enamel is an incredibly resilient biological material, withstanding repeated daily stresses for decades. The mechanisms behind this resilience remain an open question, with recent studies demonstrating a crack-deflection mechanism contributing to enamel toughness and other studies detailing the roles of the organic matrix and remineralization. Here, we focus on the mineral and hypothesize that self-healing of cracks in enamel nanocrystals may be an additional mechanism acting to prevent catastrophic failure. To test this hypothesis, we used a molecular dynamics (MD) approach to compare the fracture behavior of hydroxyapatite (HAP) and calcite, the main minerals in human enamel and sea urchin teeth, respectively. In this work, we find that cracks heal under pressures typical of mastication by fusion of crystals in HAP but not in calcite, which is consistent with the resilience of HAP enamel that calcite teeth lack. Scanning transmission electron microscopy (STEM) images of structurally intact ("sound") human enamel show dashed-line nanocracks that resemble and therefore might be the cracks healed by fusion of crystals produced in silico . The fast, self-healing mechanism shown here is common in soft materials and ceramics but has not been observed in single crystalline materials at room temperature. The crack self-healing in sound enamel nanocrystals, therefore, is unique in the human body and unique in materials science, with potential applications in designing bioinspired materials.

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

Role of α → ε → α phase transformation on the spall behavior of iron at atomic scales

Shock compression of iron microstructures above a threshold stress results in a α (BCC) → ε (HCP) transformation and the propagation of the release wave brings the metal back to the α phase following the ε→α transformation. Predicting failure behavior under shock loading conditions (spallation) relies on understanding the evolution of defects in the microstructure as it undergoes the α→ε→α phase transformation. This study uses molecular dynamics (MD) simulations to investigate the role of defect evolution during α→ε→α phase transformation on the spall strength values of single-crystal (sc) Fe microstructures. The MD simulations aim to characterize the ε phase fraction formed during shock compression and the defects during shock release for variations in loading orientations and shock stresses. The simulations are carried out for loading along the [100], [110], [111], and [112] orientations and for impact velocities ranging from 600 m/s to 1 km/s. The ε phase fractions during compression and defects (dislocations, twins) characterized during spall failure show an orientation dependence that affects the spall strength values. Here, the lowest value for spall strength is observed for the <110> loading orientation that shows a high density of twinning at the spall plane, whereas the highest value is observed for the <100> orientation and is associated with a α (BCC) → γ (FCC) transformation at the spall plane. The correlations of the spall strength values with the strain rates and with the ε phase fractions are discussed.

36 MATERIALS SCIENCE↗
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↗