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 775 records · Page 43
All sourcesNASA NTRSDOE OSTI
DOE OSTI2023

Graph-based quantum response theory and shadow Born–Oppenheimer molecular dynamics

Graph-based linear scaling electronic structure theory for quantum-mechanical molecular dynamics simulations [A. M. N. Niklasson et al., J. Chem. Phys. 144, 234101 (2016)] is adapted to the most recent shadow potential formulations of extended Lagrangian Born–Oppenheimer molecular dynamics, including fractional molecular-orbital occupation numbers [A. M. N. Niklasson, J. Chem. Phys. 152, 104103 (2020) and A. M. N. Niklasson, Eur. Phys. J. B 94, 164 (2021)], which enables stable simulations of sensitive complex chemical systems with unsteady charge solutions. The proposed formulation includes a preconditioned Krylov subspace approximation for the integration of the extended electronic degrees of freedom, which requires quantum response calculations for electronic states with fractional occupation numbers. For the response calculations, we introduce a graph-based canonical quantum perturbation theory that can be performed with the same natural parallelism and linear scaling complexity as the graph-based electronic structure calculations for the unperturbed ground state. Further, the proposed techniques are particularly well-suited for semi-empirical electronic structure theory, and the methods are demonstrated using self-consistent charge density-functional tight-binding theory both for the acceleration of self-consistent field calculations and for quantum-mechanical molecular dynamics simulations. Graph-based techniques combined with the semi-empirical theory enable stable simulations of large, complex chemical systems, including tens-of-thousands of atoms.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2023

Molecular dynamics study of carbon dioxide and nitrogen selectivity through poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) membrane

Here, a molecular dynamics simulation model was developed to comprehensively understand carbon dioxide over nitrogen (CO 2 /N 2 ) selectivity through polyphosphazene-based membrane comprising of poly[bis((methoxyethoxy)ethoxy)phosphazene] (MEEP) selective layer at the molecular level. The effects of temperature, pressure, and initial feed gas composition on the CO 2 transport on a polymer membrane were studied. The computed free energy and density profile of the permeating gas mixture exhibited that CO 2 molecules express higher interactions with the membrane than N 2 molecules, resulting in higher CO 2 /N 2 selectivity. Statistical analysis of gas molecules (CO 2 , water (H 2 O), and N 2 ) transportation suggested that hydro- and CO 2 -philic functional groups in the membrane significantly impact CO 2 permeability and CO 2 /N 2 selectivity. This study suggested that Lewis acid–base and hydrogen bonding combinations contribute to CO 2 permeation and CO 2 /N 2 selectivity. An equal CO 2 /N 2 selectivity was observed with and without water vapor in the feed gas suggesting that water does not hinder CO 2 transport through the membrane.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Curvature-induced hydrophobicity at imogolite–water interfaces

Imogolite, a nanotubular aluminosilicate mineral, is commonly found in volcanic soils, where it exerts a control on carbon dynamics. Synthetic imogolites are used for the removal of contaminants from industrial effluents and are considered for a range of other applications including gas adsorption and functionalised heterogeneous catalysts. In spite of their environmental and industrial relevance, the properties of imogolite–water interfaces remain poorly understood. In this work, an experimental and computational study is presented in which the structure and energetics of water are characterized on the curved external surface of imogolite and the hydrophilicity of this surface is contrasted with that of gibbsite, its planar counterpart. Atomic force spectroscopy experiments show that in spite of their identical surface structure, imogolite has a lower hygroscopicity than gibbsite. Molecular dynamics simulations provide an explanation for this observation: the curvature of imogolite prevents the formation of in-plane H-bonds along the directions of the nanotube circumference, lowering the enthalpy of adsorption of water molecules. The different arrangement of surface H-bonds and the resulting differences in hydration properties also affects the acidity constants of surface hydroxyl groups. This ‘nanotube effect’ may be relevant to other nanotubular systems with high curvatures, potentially impacting their wetting properties, their colloidal stability and their affinity towards hydrophobic organic moieties.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Variable domain mutational analysis to probe the molecular mechanisms of high viscosity of an IgG 1 antibody

Subcutaneous injection is the preferred route of administration for many antibody therapeutics for reasons that include its speed and convenience. However, the small volume limit (typically ≤2 mL) for subcutaneous delivery often necessitates antibody formulations at high concentrations (commonly ≥100 mg/mL), which may lead to physicochemical problems. For example, antibodies with large hydrophobic or charged patches can be prone to self-interaction giving rise to high viscosity. Here, we combined X-ray crystallography with computational modeling to predict regions of an anti-glucagon receptor (GCGR) IgG 1 antibody prone to self-interaction. An extensive mutational analysis was undertaken of the complementarity-determining region residues residing in hydrophobic surface patches predicted by spatial aggregation propensity, in conjunction with residue-level solvent accessibility, averaged over conformational ensembles from molecular dynamics simulations. Dynamic light scattering (DLS) was used as a medium throughput screen for self-interaction of ~200 anti-GCGR IgG 1 variants. A negative correlation was found between the viscosity determined at high concentration (180 mg/mL) and the DLS interaction parameter measured at low concentration (2–10 mg/mL). Additionally, anti-GCGR variants were readily identified with reduced viscosity and antigen-binding affinity within a few fold of the parent antibody, with no identified impact on overall developability. The methods described here may be useful in the optimization of other antibodies to facilitate their therapeutic administration at high concentration.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2025

Thermal disorder and phonon softening in the ferroelectric phase transition of lead titanate

We report a molecular dynamics study of ab initio quality of the ferroelectric phase transition in crystalline PbTi⁢O3. We model anharmonicity accurately in terms of potential energy and polarization surfaces trained on density functional theory data with modern machine learning techniques. Our simulations demonstrate that the transition has a strong order-disorder character, in agreement with diffraction experiments, and provide fresh insight into the approach to equilibrium across the phase transition. We find that the emergence and disappearance of the macroscopic polarization is driven by dipolar switching at the nanometer scale. We also computed the infrared optical absorption spectra in both the ferroelectric and the paraelectric phases, finding good agreement with the experimental Raman frequencies. Often, the almost ideal displacive character of the soft mode detected by Raman scattering in the paraelectric phase has been contrasted with the order-disorder character of the transition suggested by diffraction experiments. We settle this issue by showing that the soft mode coexists with a strong Debye relaxation associated with thermal disordering of the dipoles. The Debye relaxation feature is centered at zero frequency and appears near the transition temperature in both the ferroelectric and the paraelectric phases.

75 CONDENSED MATTER PHYSICS, SUPERCONDUCTIVITY AND↗
DOE OSTI2021

Effects of Al:Si and (Al+Na):Si Ratios on the Properties of the International Simple Glass, Part II: Structure

High-alumina containing high-level waste (HLW) will be vitrified at the Waste Treatment Plant at the Hanford Site. The resulting glasses, high in alumina, will have distinct composition-structure-property (C-S-P) relationships compared to previously studied HLW glasses. These C-S-P relationships determine the processability and product durability of glasses and therefore must be understood. The main purpose of this study is to understand the detailed structural changes caused by Al:Si and (Al+Na):Si substitutions in a simplified nuclear waste model glass (ISG, International Simple Glass) by combining experimental structural characterizations and molecular dynamics (MD) simulations. The structures of these two series of glasses were characterized by neutron total scattering and 27Al, 23Na, 29Si, and 11B solid-state nuclear magnetic resonance (NMR) spectroscopy. Additionally, MD simulations were used to generate atomistic structural models of the borosilicate glasses and simulation results were validated by the experimental structural data. Short-range (e.g., bond distance, coordination number, etc.) and medium-range (e.g., oxygen speciation, network connectivity, polyhedral linkages) structural features of the borosilicate glasses were systematically investigated as a function of the degree of substitution. The results show that bond distance and coordination number of the cation-oxygen pairs are relatively insensitive to Al:Si and (Al+Na):Si substitutions with the exception of the B-O pair. Additionally, the Al:Si substitution results in an increase of tri-bridging oxygen species, while (Al+Na):Si substitution creates non-bridging oxygen species. Charge compensator preferences were found for Si-[NBO] (Na+), [3]B-[NBO] (Na+), [4]B (mostly Ca2+), [4]Al (nearly equally split Na+ and Ca2+), and [6]Zr (mostly Ca2+). The network former-BO-network former linkages preferences were also tabulated; Si-O-Al and Al-O-Al were preferred at the expense of lower Si-O-[3]B and [3]B-O-[3]B linkages. These results provide insights on the structural origins of property changes such as glass transition temperature caused by the substitutions, providing a basis for future improvements of theoretical and computer simulation models.

Lu, Xiaonan↗
DOE OSTI2021

Pressure, temperature, and orientation dependent thermal conductivity of pentaerythritol tetranitrate (PETN)

We use reverse non-equilibrium molecular dynamics to determine the thermal conductivity tensor, κ, of tetragonal pentaerythritol tetranitrate (PETN). The most stable form under standard ambient conditions (PETN I, with space group P4¯2 1 c) is considered in the temperature and pressure intervals 200–500 K and 0–8 GPa, which covers approximately the stability range for this polymorph. We compute κ along the [100] and [001] directions, which is enough to construct the full thermal conductivity tensor for a system with tetragonal symmetry. In addition, we also determine κ along the [010] direction and confirm that κ 100 ~ κ 010 , with less than 4% average absolute error between the two quantities. We observe an anisotropic response for κ, with κ 100 > κ 001 across the whole (T,P) interval, and 37% difference at 300 K, 0 GPa. Furthermore, we provide analytical functions to interpolate κ(P,T) within the fitting interval and serve as input for continuum-scale simulations.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials: Part II Resveratrol Exemplar

This SAND report summarizes work supported by an Engineering Sciences Research Foundation (ESRF) Lab Directed Research and Development (LDRD) project entitled “Beyond Magic Barrels: Digital manufacturing for crystallization, process development and optimization of explosive materials.” This SAND report is written in two parts with Part 1 discusses recrystallization of our explosive exemplar and Part 2 summarizing our work with recrystallization of resveratrol. We have studied resveratrol recrystallization with a multiscale approach combining experiments, modeling and simulation. At the single crystal scale, microscopy experiments illuminate crystal time-dependent growth rates using advanced image analysis. Bench scale experiments were carried out to look at growth of multiple particles in a small reactor creating thousands of particles and analyzing the results with microscopy and μCT. For the modeling we combine kinetic Monte Carlo (kMC) models with subscale information from density functional theory (DFT) or molecular dynamics. This work is discussed in Part 1 and can also be found in a paper from the project discussing a coarse-grained kMC model specifically developed for resveratrol. For well-mixed systems, we have population balance equations (PBE) linked with species mass conservation forming a set of ordinary differential equations that can be solved quickly. For more complicated geometries, such as the vat crystallization used throughout the complex, a coupled computational fluid dynamic (CFD)/PBE method was developed to account for gradients in temperature and concentration and differences in crystallization rates throughout the domain. These simulations are more complex and require high performance computing. We present results for two cases: 5% seed fast cool with parameters fit to the well-mixed case and 5% seed slow cool using the same parameters. We show reasonable agreement with experiments though are particles are significantly larger than the experiments.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Coordinate transformation methodology for simulating quasistatic elastoplastic solids

Molecular dynamics simulations frequently employ periodic boundary conditions where the positions of the periodic images are manipulated in order to apply deformation to the material sample. For example, Lees-Edwards conditions use moving periodic images to apply simple shear. Here, we examine the problem of precisely comparing this type of simulation to continuum solid mechanics. We employ a hypo-elastoplastic mechanical model, and develop a projection method to enforce quasi-static equilibrium. We introduce a simulation framework that uses a fixed Cartesian computational grid on a reference domain, and imposes deformation via a time-dependent coordinate transformation to the physical domain. As a test case for our method, we consider the evolution of shear bands in a bulk metallic glass using the shear transformation zone theory of amorphous plasticity. We examine the growth of shear bands in simple shear and pure shear conditions as a function of the initial preparation of the bulk metallic glass.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Solid-State Calcium-Ion Diffusion in Ca 1.5 Ba 0.5 Si 5 O 3 N 6

Rechargeable batteries based on multivalent working ions are promising candidates for next-generation high-energy-density batteries. Development of these technologies, however, is largely limited by the low diffusion rate of multivalent ions in solid-state materials, thereby necessitating a better understanding of the design principles that control multivalent-ion mobility. We report Ca 1.5 Ba 0.5 Si 5 O 3 N 6 as a potential calcium solid-state conductor and investigate its Ca migration mechanism by means of ab initio computations and neutron diffraction. This compound contains partially occupied Ca sites in close proximity to each other, providing a unique mechanism for Ca migration. Nuclear density maps obtained with the maximum entropy method from neutron powder diffraction data provide strong evidence for low-energy percolating one-dimensional pathways for Ca-ion migration. Ab initio molecular dynamics simulations further support a low Ca-ion migration barrier of ~400 meV when Ca vacancies are present and reveal a unique "vacancy-adjacent"concerted ion migration mechanism. This work provides a new understanding of solid-state Ca-ion diffusion and insights into the future design of novel cation configurations that utilize the interactions between mobile ions to enable fast multivalent-ion conduction in solid-state materials.

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

Simulation of Inelastic Neutron Scattering Spectra Directly from Molecular Dynamics Trajectories

Inelastic neutron scattering (INS) is a widely used technique to study atomic and molecular vibrations. With the increasing complexity of materials and thus the INS spectra, being able to simulate the spectra from various atomistic models becomes an essential step and also a major bottleneck for INS data analysis. The conventional approach using density functional theory and lattice dynamics often falls short when the materials of interest are complex (e.g., defective, disordered, heterogeneous, amorphous, large-scale), for which molecular dynamics driven by an interatomic force field is a more common approach. In this paper, we demonstrate a method to directly convert molecular dynamics trajectories into simulated INS spectra, including not only fundamental but also higher order excitations. In this work, the results are compared with data collected on various representative samples from different neutron spectrometers. This development will open great opportunities by providing the key tool to perform in-depth analysis of INS data and to validate and optimize computer models.

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

Machine learning nonequilibrium electron forces for spin dynamics of itinerant magnets

Abstract We present a generalized potential theory for conservative as well as nonconservative forces for the Landau-Lifshitz magnetization dynamics. Importantly, this formulation makes possible an elegant generalization of the Behler-Parrinello machine learning (ML) approach, which is a cornerstone of ML-based quantum molecular dynamics methods, to the modeling of force fields in adiabatic spin dynamics of out-of-equilibrium itinerant magnetic systems. We demonstrate our approach by developing a deep-learning neural network that successfully learns the electron-mediated exchange fields in a driven s-d model computed from the nonequilibrium Green’s function method. We show that dynamical simulations with forces predicted from the neural network accurately reproduce the voltage-driven domain-wall propagation. Our work also lays the foundation for ML modeling of spin transfer torques and opens a avenue for ML-based multi-scale modeling of nonequilibrium dynamical phenomena in itinerant magnets and spintronics.

36 MATERIALS SCIENCE↗
NASA NTRSDate not supplied

ICME for NASA Aerospace Applications: Batteries for Electric Aviation

NASA’s approach to computational materials modeling is detailed in the NASA Vision 2040 Roadmap for Multiscale Modeling and Simulation of Materials and Systems. This report is in the spirit of national initiatives such as the Material Genome Initiative (MGI), Integrated Computational Materials Engineering (ICME), and others. We utilize a combination of fundamental modeling, computational high-throughput screening, and data science methods, e.g., machine learning, are used to find innovative solutions to NASA or national technology challenges. Applications of interest are wide ranging from advanced alloys to batteries to coatings, among others. In this talk, we present three examples for recent work related to NASA applications. First, doping advanced sulfur battery cathodes with selenium boosts electrical conductivity important for electric aircraft applications. First principles calculations will be discussed that result in compositional design maps for these materials. Second, development of icephobic coatings is important to mitigate safety hazards associated with icing for aircraft. Molecular dynamics simulations are reported for ice-surface interfaces to understand adhesion mechanisms and help screen optimal ice-phobic coatings. Third, shape memory alloys have numerous applications as actuators, superelastic materials, etc. for aerospace. We report machine learning models that predict martensitic transition temperatures across a broad swath of compositional space.

John Lawson↗
DOE OSTI2022

Investigation of γ - ( U , Z r ) structural properties and its interfacial properties with liquid sodium using ab initio molecular dynamics

In this study, the elastic properties, structural parameters, sound velocity, and Debye temperature of γ–(U,Zr) were computed using ab initio molecular dynamics (AIMD) at temperatures between 1000 K and 1400 K and for Zr content between 0 at. % and 100 at. %. UZr is used as a metallic fuel for Sodium Fast Reactors (SFRs). The study of the mechanical and thermal behavior of these alloys leads to a better data-informed fuel design. The bulk modulus, shear modulus, Young's modulus, and Poisson's ratio were calculated from the elastic constants and their dependence on Zr content and temperature was investigated, comparing the results with previous computational work and the available experimental data in the literature. Interfacial properties between UZr (up to 32 at. % which typically exists in nuclear fuel) and liquid sodium are also of interest due to the presence of a sodium bond between the fuel and the cladding in metallic nuclear fuel. The interfacial energy between γ–(U,Zr) and liquid sodium, the surface tension of liquid sodium, and the work of adhesion were computed at different temperatures and Zr concentrations. It was demonstrated that γ–(U,Zr) is completely wetted by liquid sodium at all the investigated temperatures and Zr concentrations. Finally, this work provides the basis for the determination of interfacial resistances in SFRs and their implementation into heat transfer fuel performance simulations, which will be the subject of future work.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Mechanisms and stability of Li dynamics in amorphous Li-Ti-P-S-based mixed ionic–electronic conductors: A machine learning molecular dynamics study

Mixed ionic–electronic conductors (MIECs) exhibit both high ionic and electronic conductivity to improve the battery performance. In this work, we investigate the mechanism and stability of transport channels in our recently developed MIEC material, amorphous Ti-doped lithium phosphorus sulfide (LPS), using molecular dynamics (MD) simulations with a 99% accurate machine-learning force field (MLFF) trained on ab initio MD data. The achieved MLFF helps efficient large-scale MD simulations on LPS with three Ti concentrations (10%, 20%, and 30%) and six temperatures (25°C to 225°C) to calculate ionic conductivity, activation energy, Li-ion transport mechanism, and configurational entropy. Results show that ionic conductivities and activation energies are consistent with our recent experimental values. Moreover, Li-ion transport occurs via free-volume diffusion facilitated by the formation of disordered Li-S polyhedra. The enhanced stability of transport channels at 10% and 20% Ti doping, compared to 0% and 30%, is observed by analyzing the vibrational and configurational entropy of these disordered Li-S polyhedra. Overall, this study highlights the utility of MLFF-based large-scale MD simulations in explaining the transport mechanism and the stability of Li-ion in Ti-doped LPS electrolyte with significant computational efficiency.

And configuration entropy↗
DOE OSTI2025

MBX V1.2: Accelerating Data-Driven Many-Body Molecular Dynamics Simulations

The MBX software provides an advanced platform for molecular dynamics simulations, leveraging state-of-the-art MB-pol and MB-nrg data-driven many-body potential energy functions. Developed over the past decade, these potential energy functions integrate physics-based and machine-learned many-body terms trained on electronic structure data calculated at the "gold standard" coupled-cluster level of theory. Recent advancements in MBX have focused on optimizing its performance, resulting in the release of MBX v1.2. While the inherently many-body nature of MB-pol and MB-nrg ensures high accuracy, it poses computational challenges. MBX v1.2 addresses these challenges with significant performance improvements, including enhanced parallelism that fully harnesses the power of modern multicore CPUs. In conclusion, these advancements enable simulations on nanosecond time scales for condensed-phase systems, significantly expanding the scope of high-accuracy, predictive simulations of complex molecular systems powered by data-driven many-body potential energy functions.

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

Acetylcholinesterase: Structure, dynamics, and interactions with organophosphorus compounds

Acetylcholinesterase (AChE) is an enzyme that hydrolyzes the neurotransmitter acetylcholine (ACh), removing it from the synaptic cleft after the transmission of an electrical signal, making it an essential component of chemical neurotransmission. AChE is a serine hydrolase, containing a catalytic triad of Ser/His/Glu. AChE is a prime target for pharmaceuticals treating a variety of neurological disorders. It is also the target of synthetic organophosphorus (OP) compounds that have been used as pesticides and chemical warfare agents. OP compounds contain a potent leaving group, such as fluorine, and act by forming a covalent adduct with the catalytic serine of the AChE active site. A wealth of structural information is available for AChE, including over 300 structures, including a subset of structures in complex with drugs as well as OP compounds. This review will highlight the interactions between OP compounds and AChE from a structural and computational perspective, with a discussion of access to the active site, as well as side reactions that lead to dealkylation of the OP-catalytic serine adduct, a process known as aging. We conclude that while the majority of the conformational changes needed to accommodate the OP compounds are localized to the acyl loop in the crystal structures, molecular dynamics simulations highlight the potential for a far more dynamic enzyme.

59 BASIC BIOLOGICAL SCIENCES↗
DOE OSTI2021

Bridging the gap between simulated and experimental ionic conductivities in lithium superionic conductors

Lithium superionic conductors (LSCs) are of major importance as solid electrolytes for next-generation all-solid-state lithium-ion batteries. While ab initio molecular dynamics have been extensively applied to study these materials, there are often large discrepancies between predicted and experimentally measured ionic conductivities and activation energies due to the high temperatures and short time scales of such simulations. Here, we present a strategy to bridge this gap using moment tensor potentials (MTPs). We show that MTPs trained on energies and forces computed using the van der Waals optB88 functional yield much more accurate lattice parameters, which in turn leads to accurate prediction of ionic conductivities and activation energies for the Li 0·33 La 0·56 TiO 3 , Li 3 YCl 6 and Li 7 P 3 S 11 LSCs. NPT MD simulations using the optB88 MTPs also reveal that all three LSCs undergo a transition between two quasi-linear Arrhenius regimes at relatively low temperatures. This transition can be traced to an increase in the number and diversity of diffusion pathways, in some cases with a change in the dimensionality of diffusion. Furthermore, this work presents not only an approach to develop high accuracy MTPs, but also outlines the diffusion characteristics for LSCs which is otherwise inaccessible through ab initio computation.

25 ENERGY STORAGE↗