Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular dynamics”

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 217 records · Page 12
All sourcesNASA NTRSDOE OSTI
DOE OSTI2025

Quantum Simulation of Molecular Dynamics Processes─A Benchmark Study Using a Classical Simulator and Present-Day Quantum Hardware

Here, we explore how the fundamental problems in quantum molecular dynamics can be modeled using classical simulators (emulators) of quantum computers and the actual quantum hardware available to us today. The list of problems we tackle includes propagation of a free wave packet, vibration of a harmonic oscillator, and tunneling through a barrier. Each of these problems starts with the initial wave packet setup. Although Qiskit provides a general method for initializing wave functions, in most cases it generates deep quantum circuits. While these circuits perform well on noiseless simulators, they suffer from excessive noise on quantum hardware. To overcome this issue, we designed a shallower quantum circuit for preparing a Gaussian-like initial wave packet, which improves the performance of real hardware. Next, quantum circuits are implemented to apply the kinetic and potential energy operators for the evolution of a wave function over time. The results of our modeling on classical emulators of quantum hardware agree perfectly with the results obtained using the traditional (classical) methods. This serves as a benchmark and demonstrates that the quantum algorithms and Qiskit codes we developed are accurate. However, the results obtained on the actual quantum hardware available today, such as IBM’s superconducting qubits and IonQ’s trapped ions, indicate large discrepancies due to hardware limitations. This work highlights both the potential and challenges of using quantum computers to solve fundamental quantum molecular dynamics problems.

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

Identification of preferred multimodal ligand‐binding regions on IgG1 F C using nuclear magnetic resonance and molecular dynamics simulations

Abstract In this study, the binding of multimodal chromatographic ligands to the IgG1 F C domain were studied using nuclear magnetic resonance and molecular dynamics simulations. Nuclear magnetic resonance experiments carried out with chromatographic ligands and a perdeuterated 15 N‐labeled F C domain indicated that while single‐mode ion exchange ligands interacted very weakly throughout the F C surface, multimodal ligands containing negatively charged and aromatic moieties interacted with specific clusters of residues with relatively high affinity, forming distinct binding regions on the F C . The multimodal ligand‐binding sites on the F C were concentrated in the hinge region and near the interface of the C H 2 and C H 3 domains. Furthermore, the multimodal binding sites were primarily composed of positively charged, polar, and aliphatic residues in these regions, with histidine residues exhibiting some of the strongest binding affinities with the multimodal ligand. Interestingly, comparison of protein surface property data with ligand interaction sites indicated that the patch analysis on F C corroborated molecular‐level binding information obtained from the nuclear magnetic resonance experiments. Finally, molecular dynamics simulation results were shown to be qualitatively consistent with the nuclear magnetic resonance results and to provide further insights into the binding mechanisms. An important contribution to multimodal ligand‐F C binding in these preferred regions was shown to be electrostatic interactions and π–π stacking of surface‐exposed histidines with the ligands. This combined biophysical and simulation approach has provided a deeper molecular‐level understanding of multimodal ligand–F C interactions and sets the stage for future analyses of even more complex biotherapeutics.

Gudhka, Ronak B.↗
DOE OSTI2022

Correlations for the specific heat capacity of ( U x Pu 1 - x ) 1 - y Gd y O 2 - z derived from molecular dynamics

We report UO 2 is the primary conventional fuel used in most nuclear reactors with Gd 2 O 3 commonly added as a burnable absorber to produce a more level power distribution in the reactor core at the beginning of operation. It can also be mixed with other actinide oxides to produce mixed oxide (MOx) fuel. In this study, molecular dynamics simulations were used to predict the specific heat capacity of Gd-doped PuO 2 , UO 2 and (U, Pu)O 2 MOx accommodating Gd 3+ substituted at cation sites via two charge compensation mechanisms - oxygen vacancy formation and the oxidation of U 4+ to U 5+ . The specific heat capacity values for PuO 2 and UO 2 are in good agreement with other studies showing a distinct peak at high temperatures - above 1800 K. As Gd 3+ is added, the peak height reduces for each composition considered. An analytical fit was applied to the data where Gd 3+ was fully charge compensated by either oxygen vacancies or U 5+ . The expression was then validated by predicting the specific heat capacity for three compositions of (Ux Pu 1-x ) 1-y Gd y O 2-z containing both oxygen vacancies and U 5+ , and compared to molecular dynamics data.

11 NUCLEAR FUEL CYCLE AND FUEL MATERIALS↗
DOE OSTI2025

Combining Reactive Quantum-Mechanical Molecular-Dynamics Simulations with Mutagenesis, Crystallography, and Enzyme Kinetics to Reveal Plausible Steps of Isocyanide Hydratase Catalysis

A complete understanding of enzyme mechanisms requires atomistic details of chemical reactions. Quantum-based molecular dynamics simulations (QMD) are a potential source of this information, but trade-offs between accuracy and computational cost have limited their use. We previously developed extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) methods that leverage a negligible compromise in accuracy to substantially decrease the cost of QMD simulations. Here, we develop a reactive QMD approach using the latest XL-BOMD formulation, which enables efficient simulations of highly reactive systems, and use it to investigate mechanisms of intermediate formation in isocyanide hydratase (ICH) catalysis. In QMD simulations, molecular analogs of ICH active site residues reacted with para-nitrophenyl isocyanide, forming a thioimidate. Analysis of simulated atomic configurational and charge dynamics revealed a pathway where protonation of the isocyanide carbon occurs prior to thioimidate formation and suggested a possible role of Asp17 as a proton donor in the early phase of ICH catalysis. To test whether the pathway seen using the reactive QMD approach might be relevant to ICH catalysis, we performed X-ray crystallography and pre-steady-state enzyme kinetics studies of wild-type and D17N mutant ICH. Both the structure and kinetics are sensitive to the D17N mutation in a manner that is consistent with the order of the reaction steps seen in the simulations. Mobile protons play essential roles in many enzymes, yet they are difficult to observe experimentally, making the ordering of proton-dependent steps ambiguous in many enzyme mechanisms. The ability to directly simulate model reactions for the design of experiments that provide information about enzyme mechanisms involving mobile protons demonstrates the significance of our reactive QMD approach and motivates further biological applications.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Dynamic and relativistic effects on Pt–Pt indirect spin–spin coupling in aqueous solution studied by ab initio molecular dynamics and two- vs four-component density functional NMR calculations

Treating 195Pt nuclear magnetic resonance parameters in solution remains a considerable challenge from a quantum chemistry point of view, requiring a high level of theory that simultaneously takes into account the relativistic effects, the dynamic treatment of the solvent–solute system, and the dynamic electron correlation. A combination of Car-Parrinello molecular dynamics (CPMD) and relativistic calculations based on two-component zeroth order regular approximation spin–orbit Kohn–Sham (2c-ZKS) and four-component Dirac–Kohn–Sham (4c-DKS) Hamiltonians is performed to address the solvent effect (water) on the conformational changes and JPtPt1 coupling. A series of bridged PtIII dinuclear complexes [L1–Pt2(NH3)4(Am)2–L2]n+ (Am = α–pyrrolidonate and pivalamidate; L = H2O, Cl−, and Br−) are studied. The computed Pt–Pt coupling is strongly dependent on the conformational dynamics of the complexes, which, in turn, is correlated with the trans influence among axial ligands and with the angle N–C–O from the bridging ligands. The J-coupling is decomposed in terms of dynamic contributions. The decomposition reveals that the vibrational and explicit solvation contributions reduce JPtPt1 of diaquo complexes (L1 = L2 = H2O) in comparison to the static gas-phase magnitude, whereas the implicit solvation and bulk contributions correspond to an increase in JPtPt1 in dihalo (L1 = L2 = X−) and aquahalo (L1 = H2O; L2 = X−) complexes. Relativistic treatment combined with CPMD shows that the 2c-ZKS Hamiltonian performs as well as 4c-DKS for the JPtPt1 coupling.

Chemistry↗
DOE OSTI2022

Investigation of the thermal conductivity of molten LiF-NaF-KF with experiments, theory, and equilibrium molecular dynamics

Molten salts are being proposed for numerous advanced energy applications, including advanced nuclear reactors, concentrating solar power plants, thermal energy storage, and fusion reactors. Accurate knowledge of the thermophysical properties of molten salts directly impact the performance of these energy systems and are essential for design and safety analyses. Thermal conductivity data for fluoride molten salts and mixtures are especially lacking. In this work, experimental measurements of thermal conductivity using the steady-state variable gap technique were performed on eutectic LiF-NaF-KF from 834 to 1195 K. The experiment accounts for radiative, convective, and conductive heat losses. In addition, theoretical and molecular dynamics models are used, from 750 K up to 1300 K, to estimate the thermal conductivity for comparison with the experimental results. The results of experiments show a weak negative deviation of thermal conductivity with temperature, unlike previous experimental results in the literature. The measured thermal conductivity magnitudes agree with the theoretical and molecular dynamics predictions, aside from the data above 1100 K, where heat losses and radiative errors are the most significant, having a 16% maximum deviation from theory. These experimental results provide new thermal conductivity data for the LiF-NaF-KF system and further validation of the predictive models. The theoretical model was used to map the composition and temperature dependent thermal conductivity of LiF-NaF-KF and the mapping’s deviation from a linear additivity estimation of thermal conductivity. Additionally, this mapping showed the highest deviations from linearity for KF-LiF rich mixtures and increasing deviation with temperature. Notably, the deviation from linearity near the LiF-NaF-KF eutectic composition was around 25%.

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

Nanomechanics of Carbon and CxByNz Nanotubes: Via a Quantum Molecular Dynamics Method

Nanomechanics of single-wall C, BN and BC$_3$ and B doped C nanotubes under axial compression and tension are investigated through a generalized tight-binding molecular dynamics (GTBMD) and {\it ab-initio} electronic structure methods. The dynamic strength of BN, BC$_3$ and B doped C nanotubes for small axial strain are comparable to each other. The main difference is in the critical strain at which structural collapse occurs. For example, even a shallow doping with B lowers the value of critical strain for C nanotubes. The critical strain for BN nanotube is found to be more than that for the similar C nanotube. Once the structural collapse starts to occur we find that carbon nanotubes irreversibly go into plastic deformation regime via the formation of tetrahedral (four-fold coordinated) bonds at the location of sharp pinches or kinks. This finding is considerably different from the classical MD (molecular dynamics) simulation results known so far. The energetics and electronic densities of states of the collapsed structures, investigated with {\it ab-initio) methods, will also be discussed.

Srivastava, Deepak↗
DOE OSTI2025

Sticking Coefficients of Fusion Reactor Impurities from Molecular Dynamics Simulations for the Design of Cryopumps

A cryopump can be utilized as an impurity removal component of a direct internal recirculation (DIR) system for the fusion fuel cycle. The DIR facilitates a low fuel inventory by continuously pumping unburnt fuel while removing impurities from the fusion exhaust stream. A cryopump can target multiple impurity species by maintaining a temperature lower than the gas triple-point temperature that promotes desublimation. The desublimation/condensation of gases in cryopumps can be characterized by the sticking coefficient, which is defined as the probability for a gas particle to stick to a (cryo-)surface upon collision. The sticking coefficient is one of the important design/operation parameters for cryopumps, and it depends on a variety of surface and gas properties. Here, in this study, molecular dynamics simulations were utilized to estimate the sticking coefficients of typical fusion gas impurity species N 2 , CO 2 , and CH 4 over a Cu surface for a range of gas temperatures and surface coverages. The molecular dynamics study showed that the sticking coefficients for gases decrease with an increase in gas temperature. The presence of a single full monolayer of condensate on the metallic surface showed an adverse effect on the sticking of gases; however the sticking improved with two full monolayers of condensate on the surface. The sticking of gases over the mixed condensate on a surface was more favorable than the condensate of the same species for N 2 and CH 4 , with an exception for CO 2 , which showed a decrease in sticking over the mixed condensate.

cryopump↗
DOE OSTI2025

Molecular Dynamics Simulations of a Liquid Metal Wetting Zirconia Suboxides

Molecular dynamics simulations were performed of a liquid metal (Ag) interacting with various zirconia suboxides to calculate the work of adhesion between the liquid metal and the suboxides. These simulations were then used to assess the impact of oxygen content on wetting. The results of this work show a wetting transition as oxygen content is reduced.

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

Molecular Dynamics Study of Creep Deformation in Nickel-based Superalloy

The present study employs molecular dynamics simulations of Ni-based superalloy to investigate the creep behavior under uniaxial compression test. Dislocation dynamics is analyzed for the nickel-based single crystal superalloy with the presence of void and with varying the distribution of gamma-prime phase The results show that multiple-void systems are more prone to yield than single-void systems and single-void systems are more prone to yield than the system withoutvoid. From the simulations, it has been determined that the creep mechanism in Ni/Ni3Al is subject to change on the applied stress depending on the distribution of gamma-prime phases change. Dislocation behavior is also observed in addition to the time evolution of dislocation nucleation and growth to quantify the propensity and mechanism of creep in nickel-based single crystal superalloy. The studies confirmed that dislocation length prior to plastic deformation decreased as the number of voids increased. In addition, the deformation behavior for finegrained nickel-based polycrystalline superalloy upon grain size variation is also investigated. Our findings show that both the movement mechanisms of the grains and grain boundaries and the relationship between grain size and yield strength are highly affected by the grain size of the materials. All the databases have then been classified to accurately predict the target class vi machine learning using the data mining tool WEKA.

creep, nickel-based superalloy, molecular dynamics↗
DOE OSTI2020

Exploring non-adiabaticity to CO reduction reaction through ab initio molecular dynamics simulation

Non-adiabatic chemical reaction refers to the electronic excitation during reactions. This effect cannot be modeled by the ground-state Born-Oppenheimer molecular dynamics (BO-MD), where the electronic structure is at the ground state for every step of ions' movement. Although the non-adiabatic effect has been explored extensively in gas phase reactions, its role in electrochemical reactions, such as water splitting and CO 2 reduction, in electrolyte has been rarely explored. On the other hand, electrochemical reactions usually involve electron transport; thus, a non-adiabatic process can naturally play a significant role. In this work, using one-step CO 2 reduction as an example, we investigated the role of the non-adiabatic effect in the reaction. The reaction barriers were computed by adiabatic BO-MD and non-adiabatic real-time time dependent density functional theory (rt-TDDFT). We found that by including the non-adiabatic effect, rt-TDDFT could increase the reaction barrier up to 6% compared to the BO-MD calculated barrier when the solvent model is used to represent water. Simulations were carried out using explicit water molecules around the reaction site under different overpotentials, and similar non-adiabatic effects were found.

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

Stress Distribution During Deformation of Polycrystalline Aluminum by Molecular-Dynamics and Finite-Element Modeling

In this paper, a multiscale modelling strategy is used to study the effect of grain-boundary sliding on stress localization in a polycrystalline microstructure with an uneven distribution of grain size. The development of the molecular dynamics (MD) analysis used to interrogate idealized grain microstructures with various types of grain boundaries and the multiscale modelling strategies for modelling large systems of grains is discussed. Both molecular-dynamics and finite-element (FE) simulations for idealized polycrystalline models of identical geometry are presented with the purpose of demonstrating the effectiveness of the adapted finite-element method using cohesive zone models to reproduce grain-boundary sliding and its effect on the stress distribution in a polycrystalline metal. The yield properties of the grain-boundary interface, used in the FE simulations, are extracted from a MD simulation on a bicrystal. The models allow for the study of the load transfer between adjacent grains of very different size through grain-boundary sliding during deformation. A large-scale FE simulation of 100 grains of a typical microstructure is then presented to reveal that the stress distribution due to grain-boundary sliding during uniform tensile strain can lead to stress localization of two to three times the background stress, thus suggesting a significant effect on the failure properties of the metal.

Yamakov, V.↗
DOE OSTI2021

Molecular dynamics of the early stages of high-temperature corrosion

In this work, we performed molecular dynamics simulations of the high-temperature oxidation of metal alloys composed of Al, Cr, and Fe and compared their behavior with that of pure Fe. The metal alloy elements (Al and Cr) segregated to the surface during oxidation, producing a lower stress gradient at the metal/environment interface compared to pure Fe. We have found that the lowered stress gradients produced in the alloy material appear to play a key role in the development of corrosion. Interfaces with lower stress gradients have reduced rates of H 2 O adsorption, especially for the ferritic (bcc) alloys. The diffusivity of oxygen and hydrogen drops more rapidly for the interfaces with reduced stress gradients. The stress gradient is also diminished when the gas pressure is increased, indicating that the Fe-Cr-Al alloy system is more resistant to oxidation than pure Fe at higher pressures. Therefore, we conclude that the lower stress gradients at the alloy/environment interface reduce the stress concentration and can slow down the rate of the initial oxide scale growth. We also compared bcc and fcc alloys with pure Fe based on our three evaluation criteria (peak stress, stress gradient, and summation of stress in the oxide scale). We found that the alloys have lower values under the three criteria compared to pure Fe. The bcc alloy has the best score under a water rich environment and the fcc alloy is proven to be better for peak stress and summation of stress in the oxide scale under an oxygen rich environment. For surface segregation to occur, we find that a minimum content of Al or Cr in the near-surface region must be achieved. We also learned that the role of Al is more important than that of Cr in terms of corrosion resistant behavior at relatively higher temperatures for the Fe-Cr-Al ternary alloys.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Exploring far-from-equilibrium ultrafast polarization control in ferroelectric oxides with excited-state neural network quantum molecular dynamics

Ferroelectric materials exhibit a rich range of complex polar topologies, but their study under far-from-equilibrium optical excitation has been largely unexplored because of the difficulty in modeling the multiple spatiotemporal scales involved quantum-mechanically. To study optical excitation at spatiotemporal scales where these topologies emerge, we have performed multiscale excited-state neural network quantum molecular dynamics simulations that integrate quantum-mechanical description of electronic excitation and billion-atom machine learning molecular dynamics to describe ultrafast polarization control in an archetypal ferroelectric oxide, lead titanate. Far-from-equilibrium quantum simulations reveal a marked photo-induced change in the electronic energy landscape and resulting cross-over from ferroelectric to octahedral tilting topological dynamics within picoseconds. The coupling and frustration of these dynamics, in turn, create topological defects in the form of polar strings. The demonstrated nexus of multiscale quantum simulation and machine learning will boost not only the emerging field of ferroelectric topotronics but also broader optoelectronic applications.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Path integral molecular dynamics: A high-fidelity approach to quantum dynamics of electrons

We investigate electron transport in the uniform electron gas using ring-polymer molecular dynamics (RPMD). Working in the weakly coupled, non-degenerate regime, we use RPMD to probe how the onset of quantum diffraction effects at high temperature reshapes electron–electron collisions and leads to a classical-to-quantum crossover in macroscopic transport properties. Static thermodynamics obtained with RPMD are consistent with the weak-coupling equation of state, confirming correct quantum Boltzmann sampling. Real-time transport extracted from mean square displacements exhibits the expected ballistic-to-diffusive transition and a systematic reduction of the electronic self-diffusivity as quantum effects strengthen, due to quantum diffraction modifying electron–electron collisions. Direct ring-polymer scattering simulations reveal diffractive “softening” of binary deflections, providing a micro-to-macro link between collision physics and diffusion. The present study establishes RPMD as a quantitative, trajectory-based tool for electron transport across the classical–quantum crossover and furnishes benchmarks for improving Coulomb-log interpolation models. We outline extensions to multi-component plasmas and a path to incorporate Fermi–Dirac statistics within path-integral dynamics.

Electronic transport↗