Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Molecular simulations”

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 163 records · Page 9
All sourcesNASA NTRSDOE OSTI
DOE OSTI2020

Development of a grain growth model for U 3 Si 2 using experimental data, phase field simulation and molecular dynamics

The purpose of this work is to develop a model for normal grain growth in U 3 Si 2 . The average grain boundary energy was determined from previously published molecular dynamics simulations. The grain growth kinetics were quantified at various temperatures by annealing nanocrystalline samples. Here, the mobility was determined by comparing phase field grain growth simulations to the experimental data.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Update 1.1 to “pysimm: A python package for simulation of molecular systems”, (PII: S2352711016300395)

This update of the pysimm application programming interface, pysimm 1.1, provides both infrastructural as well as functional updates. Moreover, improvements to the random walk application that allow it to construct polymers with controlled tacticity are highlighted. Additions to the forcefield module include an update to enable working with the family of CHARMM forcefields and automated typing with the CHARMM generalized forcefield (CGenFF). Finally, new detailed examples demonstrating new features are also provided.

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

A review of displacement cascade simulations using molecular dynamics emphasizing interatomic potentials for TPBAR components

This review explores molecular dynamics simulations for studying radiation damage in Tritium Producing Burnable Absorber Rod (TPBAR) materials, emphasizing the role of interatomic potentials in displacement cascades. Recent machine learning potentials (MLPs), trained on quantum data, enhance prediction accuracy over traditional models like EAM. We highlight temperature, PKA energy, and composition effects on damage evolution in TPBAR components, recommending suitable potentials and discussing advancements for materials in extreme radiation environments.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Interfacial Reactivity and Speciation Emerging from Na-Montmorillonite Interactions with Water and Formic Acid at 200 °C: Insights from Reactive Molecular Dynamics Simulations, Infrared Spectroscopy, and X-ray Scattering Measurements

Reactive organic fluid - mineral interactions at elevated temperatures contribute to the evolution of planetary matter. One of the less studied but important transformations in this regard involves the reactions of formic acid with naturally occurring clays such as sodium montmorillonite. To advance a mechanistic understanding of these interactions, we use ReaxFF reactive molecular dynamics simulations in conjunction with infrared (IR) spectroscopy and X-ray scattering experiments to investigate the speciation behavior of water-formic acid mixtures on sodium montmorillonite interfaces at 473 K and 1 atm. Using a newly developed reactive forcefield, we show that the experimental IR spectra of unreacted and reacted mixture can be accurately reproduced by ReaxFF/MD. We further benchmark the simulation predictions of sodium carbonate and bicarbonate formation in the clay interlayers using Small and Wide-Angle X-ray Scattering measurements. Subsequently, leveraging the benchmarked forcefield, we interrogate the pathway of speciation reactions with emphasis on carbonate, formate, and hydroxide groups elucidating the energetics, transition states, intermediates, and preferred products. Further, we also delineate the differences in reactivities and catalytic effects of clay edges, facets, and interlayers owing to their local chemical environments, which have far reaching consequences in their speciation behavior. The experimental and simulation approaches described in this study and the transferable forcefields can be applied translationally to advance the science of clay-fluid interactions for several applications including subsurface fluid storage and recovery and clay-pollutant dynamics.

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

Molecular dynamics simulation studies of 1,3-dimethyl imidazolium nitrate ionic liquid with water

The fundamental understanding of intermolecular interactions of ionic liquids (ILs) with water is essential in predicting IL–water thermodynamic properties. In this study, intermolecular or noncovalent interactions were studied for 1,3-dimethyl imidazolium [DMIM]+ cation and nitrate [NO 3 ] – anion with water, employing quantum mechanics and molecular dynamics simulations. Molecular dynamics simulations were performed using a revised multipolar polarizable force field. The effect of water on ionic liquids was evaluated in terms of thermodynamic and dynamic properties. Thermodynamic properties included liquid densities ρ, excess molar volumes ΔV E , and liquid structures g(r)⁠. Dynamic properties included self-diffusion coefficients D of mixture constituents as a function of water concentration. The density of ionic liquid–water mixtures monotonically decrease with increasing concentration of water. A negative excess volume was obtained for low and high water concentrations, demonstrating strong intermolecular interactions of water with ionic liquid components. Liquid structures of ionic liquid–water mixtures revealed a tendency for anions to interact with cations at shorter intermolecular distances when the water concentration is increased. Diffusion rates were found to increase for all mixture components with increase in the fraction of water. A significant change in the diffusion rate was found at ~0.3 weight fraction of water. However, the water self-diffusion coefficient was dominant at all concentrations. Furthermore, the ratio of water/anion and anion/cation self-diffusion coefficients was found to decrease linearly with increasing concentration of water molecules.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Accelerators for Classical Molecular Dynamics Simulations of Biomolecules

Atomistic Molecular Dynamics (MD) simulations provide researchers the ability to model biomolecular structures such as proteins and their interactions with drug-like small molecules with greater spatiotemporal resolution than is otherwise possible using experimental methods. MD simulations are notoriously expensive computational endeavors that have traditionally required massive investment in specialized hardware to access biologically relevant spatiotemporal scales. Our goal is to summarize the fundamental algorithms that are employed in the literature to then highlight the challenges that have affected accelerator implementations in practice. We consider three broad categories of accelerators: Graphics Processing Units (GPUs), Field-Programmable Gate Arrays (FPGAs), and Application Specific Integrated Circuits (ASICs). These categories are comparatively studied to facilitate discussion of their relative trade-offs and to gain context for the current state of the art. We conclude by providing insights into the potential of emerging hardware platforms and algorithms for MD.

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

Near-surface damage and mixing in Si-Cl 2 -Ar atomic layer etching processes: Insights from molecular dynamics simulations

Silicon-chlorine-argon (Si-Cl 2 -Ar) atomic layer etching (ALE) is simulated using classical molecular dynamics (MD). The simulations provide a detailed view into the near-surface region during ALE processing. Bombardment of Ar + ions creates a mixed amorphous region that significantly differs from the picture of ideal ALE. There is also a significant change in the Si etch yield and the etch product distribution as a function of Ar + ion fluence. The Si etch yield is the highest at the beginning of the bombardment step but eventually decays to the physical sputtering yield. Atomic Cl and silicon chlorides are major etch products at the start of an ion bombardment step, but quickly decay. Atomic Si yields remain relatively constant as a function of Ar + ion fluence. As a result, a new schematic of Si-Cl 2 -Ar ALE is presented in order to emphasize the complex behavior observed in MD simulations.

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

The selective heating effect of microwave irradiation on a binary mixture of water and polyethylene oxide: a molecular dynamics simulation approach

In this study, we investigate the molecular mechanisms of a microwave-driven selective heating process by performing molecular dynamics simulations for three different systems including pure water, pure polyethylene oxide (PEO), and water-PEO mixed systems in the presence of a microwave with two different intensities of electric field such as 0.001 V/Å -1 and 0.01 V/Å -1 at a frequency of 100 GHz. First, from performing molecular dynamics simulations of CO and CO 2 in the presence of the microwave, it is confirmed that the molecular dipole moment is responsible for the rotational motion induced by the oscillating electric field. Second, by analyzing the MD simulations of the pure water system, we discover that the dipole moment of water exhibits a time lag with respect to the microwave. During the heating process, however, the temperature, kinetic, and potential energies increase synchronously with the oscillating electric field of the microwave, showing that the heating of the water system is caused by the molecular reaction of water to the microwave. Comparing the water-PEO mixed system to the pure water and pure PEO systems, the water-PEO mixed system has a higher heating rate than the pure PEO system but a lower heating rate than the pure water system. Therefore, we conclude that heating the water-PEO mixed system is driven by water molecules selectively activated by microwave irradiation. Furthermore, we also calculate the diffusion coefficients of water molecules and PEO chains by describing their mean square displacements, demonstrating that the diffusion coefficients are increased in the presence of microwaves for both water and PEO in pure and mixed systems. Lastly, during the microwave heating process, the structures of the water-PEO mixed system are altered as a function of the intensity of electric field, which is mainly driven by the response of water molecules.

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

Minimal Experimental Bias on the Hydrogen Bond Greatly Improves Ab Initio Molecular Dynamics Simulations of Water

Experiment Directed Simulations (EDS) is a method within a class of techniques seeking to improve molecular simulations by minimally biasing the system Hamiltonian to reproduce certain experimental observables. In a previous application of EDS to ab initio molecular dynamics (AIMD) simulation based on electronic density functional theory (DFT), the AIMD simulations of water were biased to reproduce its experimentally derived solvation structure. In particular, by solely biasing the O-O pair correlation functions, other structural and dynamical properties that were not biased were improved. In this work, the hypothesis is tested that directly biasing the OH pair correlation (and hence the H-O∙∙∙H hydrogen bonding), will provide an even better improvement of DFT-based water properties in AIMD simulations. The logic behind this hypothesis is that for most electronic DFT descriptions of water the hydrogen bonding is known to be deficient due to anomalous charge transfer and over polarization in the DFT. Using recent advances to the EDS learning algorithm, we thus train a minimal bias on AIMD water that reproduces the O-H radial distribution function derived from the highly accurate MB-pol model of water. Finally, it is then confirmed that biasing the O-H pair correlation alone can lead to improved AIMD water properties, with structural and dynamical properties in even closer to experiment than the previous EDS-AIMD model.

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

Molecular Dynamics Simulation of Pore-Size Effects on Gas Adsorption Kinetics in Zeolites

Strong gas-mineral interactions or slow adsorption kinetics require a molecular-level understanding of both adsorption and diffusion for these interactions to be properly described in transport models. In this combined molecular simulation and experimental study, noble gas adsorption and mobility is investigated in two naturally abundant zeolites whose pores are similar in size (clinoptilolite) and greater than (mordenite) the gas diameters. Simulated adsorption isotherms obtained from grand canonical Monte Carlo simulations indicate that both zeolites can accommodate even the largest gas (Rn). However, gas mobility in clinoptilolite is significantly hindered at pore-limiting window sites, as seen from molecular dynamics simulations in both bulk and slab zeolite models. Further, experimental gas adsorption isotherms for clinoptilolite confirm the presence of a kinetic barrier to Xe uptake, resulting in the unusual property of reverse Kr/Xe selectivity. Finally, a kinetic model is used to fit the simulated gas loading profiles, allowing a comparison of trends in gas diffusivity in the zeolite pores.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Quantum Simulation of Molecular Electronic States with a Transcorrelated Hamiltonian: Higher Accuracy with Fewer Qubits

Simulation of electronic structure is one of the most promising applications on noisy intermediate-scale quantum (NISQ) era devices. However, NISQ devices suffer from a number of challenges like limited qubit connectivity, short coherence times, and sizable gate error rates. Thus, desired quantum algorithms should require shallow circuit depths and low qubit counts to take advantage of these devices. Here, we attempt to reduce quantum resource requirements for molecular simulations on a quantum computer while maintaining the desired accuracy with the help of classical quantum chemical theories of canonical transformation and explicit correlation. In this work, compact ab initio Hamiltonians are generated classically, in the second quantized form, through an approximate similarity transformation of the Hamiltonian with (a) an explicitly correlated two-body unitary operator with generalized pair excitations that remove the Coulombic electron–electron singularities from the Hamiltonian and (b) a unitary one-body operator to efficiently capture the orbital relaxation effects required for accurate description of the excited states. The resulting transcorrelated Hamiltonians are able to describe both the ground and the excited states of molecular systems in a balanced manner. Using the variational quantum eigensolver (VQE) method based on the unitary coupled cluster with singles and doubles (UCCSD) ansatz and only a minimal basis set (ANO-RCC-MB), we demonstrate that the transcorrelated Hamiltonians can produce ground state energies comparable to the reference CCSD energies with the much larger cc-pVTZ basis set. This leads to a reduction in the number of required CNOT gates by more than 3 orders of magnitude for the chemical species studied in this work. Furthermore, using the quantum equation of motion (qEOM) formalism in conjunction with the transcorrelated Hamiltonian, we are able to reduce the deviations in the excitation energies from the reference EOM-CCSD/cc-pVTZ values by an order of magnitude. In conclusion, the transcorrelated Hamiltonians developed here are Hermitian and contain only one- and two-body interaction terms and thus can be easily combined with any quantum algorithm for accurate electronic structure simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗
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 OSTI2021

Crystallization of the P 3 Sn 4 Phase upon Cooling P 2 Sn 5 Liquid by Molecular Dynamics Simulation Using a Machine Learning Interatomic Potential

We performed molecular dynamics simulations to study the crystallization of the P3Sn4 phase from P 2 Sn 5 liquid using a machine learning (ML) interatomic potential with desirable efficiency and accuracy. Our results capture the liquid properties of P 2 Sn 5 at 1300 K, which is well above the melting temperature. The phase separation and crystallization are observed when P 2 Sn 5 liquid is cooled down below 832 and 505 K, respectively. The simulation results are in good agreement with the experimentally observed phase transformation behaviors and provide useful insights into the complex nucleation and crystallization process at the details of atomistic scale. Our work also demonstrated that ML interatomic potentials based on neural network deep learning are robust and capable of accurately describing the energetics and kinetics of complex materials through molecular dynamics simulations.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗