SEARCH · Search NASA

Results for “molecular dynamic 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.

SourceResearch connectionFrom yearThrough year

Quote a phrase for an exact phrase match. Source license links do not imply unrestricted reuse.

At least 361 records · Page 20

All sources NASA NTRS DOE OSTI

DOE OSTI2024

Thermally frustrated phase transition at high pressure in B 2-ordered FeV

X-ray diffraction measurements of equiatomic B2-ordered FeV were performed in a diamond-anvil cell at room temperature at several pressure points up to 80 GPa that showed the cubic phase to be stable with no indication of structural phase transitions. Density functional theory at 0 K predicts Fermi surface nesting, an electronic topological transition, and a phonon dynamical instability within the experimentally investigated pressure range. Nevertheless, the instability is absent in phonon dispersion curves extracted from ab initio molecular dynamics simulations below the critical volume at temperatures as low as 10 K, indicating that thermal atomic displacements can frustrate the phase transition by renormalizing the phonon dispersion curves. Ferrimagnetism is critical for the stability of the cubic phase at low temperature, but thermal atomic displacements are enough to support the structure at and above the Néel temperature.

36 MATERIALS SCIENCE↗

NASA NTRS2004

Unfolding of Proteins: Thermal and Mechanical Unfolding

We have employed a Hamiltonian model based on a self-consistent Gaussian appoximation to examine the unfolding process of proteins in external - both mechanical and thermal - force elds. The motivation was to investigate the unfolding pathways of proteins by including only the essence of the important interactions of the native-state topology. Furthermore, if such a model can indeed correctly predict the physics of protein unfolding, it can complement more computationally expensive simulations and theoretical work. The self-consistent Gaussian approximation by Micheletti et al. has been incorporated in our model to make the model mathematically tractable by signi cantly reducing the computational cost. All thermodynamic properties and pair contact probabilities are calculated by simply evaluating the values of a series of Incomplete Gamma functions in an iterative manner. We have compared our results to previous molecular dynamics simulation and experimental data for the mechanical unfolding of the giant muscle protein Titin (1TIT). Our model, especially in light of its simplicity and excellent agreement with experiment and simulation, demonstrates the basic physical elements necessary to capture the mechanism of protein unfolding in an external force field.

Hur, Joe S.↗

DOE OSTI2024

Proton Collective Quantum Tunneling Induces Anomalous Thermal Conductivity of Ice under Pressure

Proton tunneling is believed to be nonlocal in ice, but its range has been shown to be limited to only a few molecules. Here, we measured the thermal conductivity of ice under pressure up to 50 GPa and found it increases with pressure until 20 GPa but decreases at higher pressures. We attribute this nonmonotonic thermal conductivity to the collective tunneling of protons at high pressures, supported by large-scale quantum molecular dynamics simulations. The collective tunneling loops span several picoseconds in time and are as large as nanometers in space, which match the phonon periods and wavelengths, leading to strong phonon scattering at high pressures. Finally, our results show direct evidence of global quantum motion existing in high-pressure ice and provide a new perspective to understanding the coupling between phonon propagation and atomic tunneling.

36 MATERIALS SCIENCE↗

NASA NTRS2009

The Development of Directional Decohesion Finite Elements for Multiscale Failure Analysis of Metallic Polycrystals

Atomistic simulations of intergranular fracture have indicated that grain-scale crack growth in polycrystalline metals can be direction dependent. At these material length scales, the atomic environment greatly influences the nature of intergranular crack propagation, through either brittle or ductile mechanisms, that are a function of adjacent grain orientation and direction of crack propagation. Methods have been developed to obtain cohesive zone models (CZM) directly from molecular dynamics simulations. These CZMs may be incorporated into decohesion finite element formulations to simulate fracture at larger length scales. A new directional decohesion element is presented that calculates the direction of Mode I opening and incorporates a material criterion for dislocation emission based on the local crystallographic environment to automatically select the CZM that best represents crack growth. The simulation of fracture in 2-D and 3-D aluminum polycrystals is used to illustrate the effect of parameterized CZMs and the effectiveness of directional decohesion finite elements.

Saether, Erik↗

DOE OSTI2026

Computational simulation of asymmetric phase transformation in cracked Li₇La₃Zr₂O₁₂: Variant selection and chemo-mechanical implications

Coupling between microcracks and phase transformation in ion-conducting ceramics can jointly affect mechanical responses and ion transport. In this work we investigate the cubic-to-tetragonal phase transformation in Li₇La₃Zr₂O₁₂ in the presence of a microcrack under hydrostatic loading and quantify its implications for crack-tip stress concentration and Li-ion transport using phase-field and molecular dynamics simulations. The phase transformation exhibits a strong asymmetry between hydrostatic tension and compression. Under tension, the crack edge nucleates one tetragonal variant that amplifies the crack-tip stress intensity and promotes crack opening. Under compression, the crack tip nucleates a different tetragonal variant that enhances the stress-induced crack-closure tendency. Effective Li diffusivity analysis shows faster transport degradation under compression due to accelerated transformation kinetics, exposing a trade-off between mechanical stability and ionic conductivity. These results highlight the intertwined nature of cracking, phase transformation, and ionic transport in ion-conducting oxides and provide mechanistic insights into chemo-mechanical degradation of solid electrolytes.

36 MATERIALS SCIENCE↗

DOE OSTI2025

Defluorination Mechanisms and Real-Time Dynamics of Per- and Polyfluoroalkyl Substances on Electrified Surfaces

Per- and polyfluoroalkyl substances (PFAS) are persistent environmental contaminants found in groundwater sources and a wide variety of consumer products. In recent years, electrochemical approaches for the degradation of these harmful contaminants have garnered a significant amount of attention due to their efficiency and chemical-free modular nature. However, these electrochemical processes occur in open, highly non-equilibrium systems, and a detailed understanding of PFAS degradation mechanisms in these promising technologies is still in its infancy. To shed mechanistic insight into these complex processes, we present the first constant-electrode potential (CEP) quantum calculations of PFAS degradation on electrified surfaces. These advanced CEP calculations provide new mechanistic details about the intricate electronic processes that occur during PFAS degradation in the presence of an electrochemical bias, which cannot be gleaned from conventional density functional theory calculations. We complement our CEP calculations with large-scale ab initio molecular dynamics simulations in the presence of an electrochemical bias to provide time scales for PFAS degradation on electrified surfaces. Taken together, our CEP-based quantum calculations provide critical reaction mechanisms for PFAS degradation in open electrochemical systems, which can be used to prescreen candidate material surfaces and optimal electrochemical conditions for remediating PFAS and other environmental contaminants.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Data-driven modeling of dislocation mobility from atomistics using physics-informed machine learning

Dislocation mobility, which dictates the response of dislocations to an applied stress, is a fundamental property of crystalline materials that governs the evolution of plastic deformation. Traditional approaches for deriving mobility laws rely on phenomenological models of the underlying physics, whose free parameters are in turn fitted to a small number of intuition-driven atomic scale simulations under varying conditions of temperature and stress. This tedious and time-consuming approach becomes particularly cumbersome for materials with complex dependencies on stress, temperature, and local environment, such as body-centered cubic crystals (BCC) metals and alloys. In this paper, we present a novel, uncertainty quantification-driven active learning paradigm for learning dislocation mobility laws from automated high-throughput large-scale molecular dynamics simulations, using Graph Neural Networks (GNN) with a physics-informed architecture. We demonstrate that this Physics-informed Graph Neural Network (PI-GNN) framework captures the underlying physics more accurately compared to existing phenomenological mobility laws in BCC metals.

36 MATERIALS SCIENCE↗

DOE OSTI2026

Plasma hydrodynamics from mean force kinetic theory

Mean force kinetic theory is used to evaluate the electrical conductivity, thermal conductivity, electrothermal coefficient, thermoelectric coefficient, and shear viscosity of a two-component plasma. The results are compared with molecular dynamics simulations. These simulations are made possible by assuming a repulsive Coulomb force for all interactions. Good agreement is found for all coefficients up to a Coulomb coupling strength of Γ ≈ 20. This is over 100-times larger than the coupling strength at which traditional theories break down. It is concluded that mean force kinetic theory is a promising framework for extending hydrodynamics to dense plasmas.

Electrical conductivity↗

NASA NTRS1998

Electronic Properties of Carbon Nanotubes and Junctions

Metallic and semiconducting Single Wall Carbon Nanotubes (CNT) have recently been characterized using scanning tunneling microscopy (STM) and the manipulation of individual CNT has been demonstrated. These developments make the prospect of using CNT as molecular wires and possibly as electronic devices an even more interesting one. We have been modeling various electronic properties such as the density of states and the transmission coefficient of CNT wires and junctions. These studies involve first calculating the stability of junctions using molecular dynamics simulations and then calculating the electronic properties using a pi-electron tight binding Hamiltonian. We have developed the expertise to calculate the electronic properties of both finite-sized CNT and CNT systems with semi-infinite boundary conditions. In this poster, we will present an overview of some of our results. The electronic application of CNT that is most promising at this time is their use as molecular wires. The conductance can however be greatly reduced because of reflection due to defects and contacts. We have modeled the transmission through CNT in the presence of two types of defects: weak uniform disorder and strong isolated scatterers. We find that the conductance is affected in significantly different manners due to these defects Junctions of CNT have also been imaged using STM. This makes it essential to derive rules for the formation of junctions between tubes of different chirality, study their relative energies and electronic properties. We have generalized the rules for connecting two different CNT and have calculated the transmission and density of states through CNT junctions. Metallic and semiconducting CNT can be joined to form a stable junction and their current versus voltage characteristics are asymmetric. CNT are deformed by the application of external forces including interactions with a substrate or other CNT. In many experiments, these deformation are expected to occur naturally. We will present some preliminary results of our calculations of the modification of CNT electronic properties as a result of deformations.

Anantram, M. P.↗

DOE OSTI2025

Ferroelectric phase transition in group-IV monochalcogenides from an equivariant machine learned force field

Group-IV monochalcogenides are a class of layered ferroelectric semiconductors that have demonstrated spontaneous intrinsic polarization above room temperature. Here, in this study, we use the multi-atomic cluster expansion (MACE) machine learning architecture to train and test a force field capable of modeling the structural properties and second-order ferroelectric-to-paraelectric phase transition in a Group-IV monochalcogenide, GeSe. The model captures the double-well potential energy surface associated with the onset of macroscopic polarization in bulk GeSe within 12.5 meV/atom, as well as near-equilibrium properties like the phonon dispersion. The development of this quantitatively accurate force field enables long-time molecular dynamics simulations, which predict the critical temperature of the ferroelectric-to-paraelectric phase transition in bulk GeSe to be T c = 600 K. This study demonstrates the capabilities of equivariant force-fields to accurately describe phenomena associated with structural symmetry breaking.

ferroelectricity↗

NASA NTRS2006

Multiscale Modeling of Grain-Boundary Fracture: Cohesive Zone Models Parameterized From Atomistic Simulations

A multiscale modeling strategy is developed to study grain boundary fracture in polycrystalline aluminum. Atomistic simulation is used to model fundamental nanoscale deformation and fracture mechanisms and to develop a constitutive relationship for separation along a grain boundary interface. The nanoscale constitutive relationship is then parameterized within a cohesive zone model to represent variations in grain boundary properties. These variations arise from the presence of vacancies, intersticies, and other defects in addition to deviations in grain boundary angle from the baseline configuration considered in the molecular dynamics simulation. The parameterized cohesive zone models are then used to model grain boundaries within finite element analyses of aluminum polycrystals.

Glaessgen, Edward H.↗

DOE OSTI2024

Structural Characterization and Dynamics of AdhE Ultrastructures from Clostridium thermocellum Show a Containment Strategy for Toxic Intermediates

Clostridium thermocellum, a cellulolytic thermophilic anaerobe, is considered by many to be a prime candidate for the realization of consolidated bioprocessing (CBP) and is known as an industry standard for biofuel production. C. thermocellum is among the best biomass degraders identified to date in nature and produces ethanol as one of its main products. Many studies have helped increase ethanol titers in this microbe; however, ethanol production using C. thermocellum is still not economically viable. Therefore, a better understanding of its ethanol synthesis pathway is required. The main pathway for ethanol production in C. thermocellum involves the bifunctional aldehyde-alcohol dehydrogenase (AdhE). To better understand the function of the C. thermocellum AdhE, we used cryo-electron microscopy (cryo-EM) to obtain a 3.28 A structure of the AdhE complex. This high-resolution structure, in combination with molecular dynamics simulations, provides insight into the substrate channeling of the toxic intermediate acetaldehyde, indicates the potential role of C. thermocellum AdhE to regulate activity and cofactor pools, and establishes a basis for future engineering studies. The containment strategy found in this enzyme offers a template that could be replicated in other systems where toxic intermediates need to be sequestered to increase the production of valuable biochemicals.

09 BIOMASS FUELS↗

DOE OSTI2024

Demystify radiation-enhanced hydrogen isotope diffusion in Fe-Ni-Cr austenitic stainless steels

Understanding and containing hydrogen isotope diffusion is crucial for many nuclear applications. In situ experiments have consistently shown that radiation significantly enhances isotope diffusion in austenitic stainless steels. Despite extensive research, the mechanism behind this phenomenon remains elusive, as most radiation-induced defects (e.g., vacancies, dislocations, and grain boundaries) typically trap hydrogen, thereby slowing diffusion. While grain boundaries may increase in-plane diffusivity and interstitials may enhance diffusion due to material swelling, these effects are relatively minor. Utilizing an Fe-Ni-Cr-H interatomic potential for stainless steels, we conducted extensive molecular dynamics simulations to investigate the origins of radiation-enhanced diffusion. Here, our findings reveal that when a system is resolidified, mimicking defects created by radiation displacements, the resulting structure contains a mixture of phases, boundaries, and dislocation networks. This defective structure significantly increases hydrogen diffusivity, enhancing it by approximately 1.7 times at 900 K. These results suggest that the complex defect structures formed during radiation displacements are the primary drivers of the observed diffusion enhancement, providing valuable insights into the mechanisms underlying radiation-enhanced diffusion in nuclear materials.

36 MATERIALS SCIENCE↗

DOE OSTI2024

Enhanced Activity from Coordinatively Unsaturated and Dynamic Zeolite-Bound Organoscandium Species

Scandium borohydride grafted into the micropores of the faujasite zeolite HY30 catalyzes the C–H borylation of benzene, whereas silica-grafted species are inactive. This catalytic activity may originate from grafting at a Brønsted acid site leading to a more electron-deficient rare earth center. Herein, we apply multinuclear double-resonance nuclear magnetic resonance (NMR) experiments to probe the structure and dynamics of zeolite- and silica-bound scandium borohydride complexes. The experiments reveal that scandium centers located within the zeolite micropores, in proximity to Al-created Brønsted sites, are more dynamic than rigid scandium sites grafted on silanols. Through a combination of NMR and molecular dynamics simulations, we show that the coordination of the scandium in the zeolite is labile, with the metal exchanging between two binding sites. As a result, the weak electron donation from the support that enables the movement of the Sc center leads to the formation of an undercoordinated metal center that cannot exist on silica, ultimately leading to the new catalytic activity of the species.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗

DOE OSTI2024

Parametrized ion-distribution model for extended x-ray absorption fine-structure analysis at high-energy-density conditions

Experiments today can compress solids near isentropically to pressures approaching 100 × 106 atmospheres; however, determining the temperature of such matter remains a major challenge. Extended x-ray absorption fine-structure (EXAFS) spectroscopy is one of the few techniques sensitive to the bulk temperature of highly compressed solid matter, and the validity of this temperature measurement relies on constraining the local ion structure from the EXAFS spectrum. At high-energy-density (HED) conditions, the local ion structure often becomes distorted, which must be accounted for during the EXAFS analysis. Described here is a technique, using a parametrized ion-distribution model to directly analyze EXAFS spectra that provides a better constraint on the local structure than traditional second- or third-order cumulant expansion techniques at HED conditions. The parametrized ion-distribution model is benchmarked by analyzing EXAFS spectra from nickel molecular-dynamics simulations at ∼100 GPa and shown to provide a 10%–20% improvement in constraining the cumulants of the true ion distribution.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗

DOE OSTI2025

Photoinduced electronic excitations drive polymerization of carbon monoxide: A first-principles study

Under pressure, carbon monoxide (CO) transforms into a polymer that can be recovered to ambient conditions. While this transformation can occur without additional stimuli, experimental observations have shown that laser irradiation can induce a similar transformation at reduced pressure. The resulting polymeric phase, which is metastable under ambient conditions, releases energy through decomposition into more stable configurations. Using time-dependent density functional theory and Born–Oppenheimer molecular dynamics simulations, we investigate the mechanism by which electronic excitation facilitates CO polymerization. Our calculations reveal that electronic excitation enhances carbon–carbon bonding, enabling polymerization at pressures significantly lower than those required by conventional compression methods. In conclusion, these findings suggest that a photo-assisted approach could be employed to synthesize novel, potentially energetic materials under less demanding pressure conditions.

Born-Oppenheimer molecular dynamics↗

DOE OSTI2024

Molecular dynamics study of hydrogen Cottrell atmosphere in aluminum: Influence of solute-solute interactions in the dislocation core

Cottrell atmospheres form when solute atoms segregate to dislocations, and is one of the most basic processes by which solutes alter mechanical properties of materials. In the case of hydrogen, Cottrell atmospheres are believed to contribute to hydrogen embrittlement. Here, in this study, we use direct molecular dynamics simulations to study the formation of hydrogen Cottrell atmospheres around edge dislocations in aluminum. Using coarse-graining techniques, we resolve the temporal and spatial evolution of the hydrogen concentration field. By comparing our results with theories of Cottrell atmosphere formation and thermodynamics, we show that a repulsive, concentration-dependent H-H interaction reduces the peak atmosphere concentration at the dislocation core. Furthermore, this interaction differs significantly from the H-H interaction that is experienced in the bulk away from a dislocation.

Aluminum↗

NASA NTRS2024

Multiscale Modeling of Thermoplastics Using Atomistic-informed Micromechanics

A multiscale repeating unit cell model of a single spherulite containing four disparate length scales was developed to predict the thermoelastic behavior of semicrystalline thermoplastic materials for composite aerospace applications. The continuum level scales were fully coupled and modeled using the generalized method of cells and the high fidelity generalized method of cells micromechanics theories. Data from molecular dynamics simulations were used as inputs for the amorphous and crystalline constituents in the multiscale continuum models. Effective Young’s modulus, shear modulus, Poisson’s ratio, coefficient of thermal expansion, and thermal conductivity were predicted for polyether ether ketone and polyether ketone ketone, showing good agreement with the available experimental data from the open literature. Moreover, it is shown that predicted properties are fairly insensitive to the fidelity of the micromechanics model used at the highest continuum scale or the assumed shape of the spherulite.

thermoplastics↗