Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “accelerated 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 145 records · Page 8
All sourcesNASA NTRSDOE OSTI
DOE OSTI2024

Atomic Diffusivities of Yttrium, Titanium and Oxygen Calculated by Ab Initio Molecular Dynamics in Molten 316L Oxide-Dispersion-Strengthened Steel Fabricated via Additive Manufacturing

Oxide-dispersion-strengthened (ODS) steels have long been viewed as a prime solution for harsh environments. However, conventional manufacturing of ODS steels limits the final product geometry, is difficult to scale up to large components, and is expensive due to multiple highly involved, solid-state processing steps required. Additive manufacturing (AM) can directly incorporate dispersion elements (e.g., Y, Ti and O) during component fabrication, thus bypassing the need for an ODS steel supply chain, the scale-up challenges of powder processing routes, the buoyancy challenges associated with casting ODS steels, and the joining issues for net-shape component fabrication. In the AM process, the diffusion of the dispersion elements in the molten steel plays a key role in the precipitation of the oxide particles, thereby influencing the microstructure, thermal stability and high-temperature mechanical properties of the resulting ODS steels. In this work, the atomic diffusivities of Y, Ti, and O in molten 316L stainless steel (SS) as functions of temperature are determined by ab initio molecular dynamics simulations. The latest Vienna Ab initio Simulation Package (VASP) package that incorporates an on-the-fly machine learning force field for accelerated computation is used. At a constant temperature, the time-dependent coordinates of the target atoms in the molten 316L SS were analyzed in the form of mean square displacement in order to obtain diffusivity. The values of the diffusivity at multiple temperatures are then fitted to the Arrhenius form to determine the activation energy and the pre-exponential factor. Given the challenges in experimental measurement of atomic diffusivity at such high temperatures and correspondingly the lack of experimental data, this study provides important physical parameters for future modeling of the oxide precipitation kinetics during AM process.

Chemistry↗
DOE OSTI2024

Acceleration of Solvation Free Energy Calculation via Thermodynamic Integration Coupled with Gaussian Process Regression and Improved Gelman–Rubin Convergence Diagnostics

The determination of the solvation free energy of ions and molecules holds profound importance across a spectrum of applications spanning chemistry, biology, energy storage, and the environment. Molecular dynamics simulations are powerful tools for computing this critical parameter. Nevertheless, the accurate and efficient calculation of the solvation free energy becomes a formidable endeavor when dealing with complex systems characterized by potent Coulombic interactions and sluggish ion dynamics and, consequently, slow transition across various metastable states. Here, in the present study, we expose limitations stemming from the conventional calculation of the statistical inefficiency g in the thermodynamic integration method, a factor that can hinder the determination of convergence of the solvation free energy and its associated uncertainty. Instead, we propose a robust scheme based on Gelman–Rubin convergence diagnostics. We leverage this improved estimation of uncertainties to introduce an innovative accelerated thermodynamic integration method based on the Gaussian Process regression. This methodology is applied to the calculation of the solvation free energy of trivalent rare-earth elements immersed in ionic liquids, a scenario in which the aforementioned challenges render standard approaches ineffective. The proposed method proves to be effective in computing solvation free energy in situations where traditional thermodynamic integration methods fall short.

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

Semi-Empirical Shadow Molecular Dynamics: A PyTorch Implementation

Here, extended Lagrangian Born–Oppenheimer molecular dynamics (XL-BOMD) in its most recent shadow potential energy version has been implemented in the semiempirical PyTorch-based software PySeQM. The implementation includes finite electronic temperatures, canonical density matrix perturbation theory, and an adaptive Krylov subspace approximation for the integration of the electronic equations of motion within the XL-BOMB approach (KSA-XL-BOMD). The PyTorch implementation leverages the use of GPU and machine learning hardware accelerators for the simulations. The new XL-BOMD formulation allows studying more challenging chemical systems with charge instabilities and low electronic energy gaps. The current public release of PySeQM continues our development of modular architecture for large-scale simulations employing semi-empirical quantum-mechanical treatment. Applied to molecular dynamics, simulation of 840 carbon atoms, one integration time step executes in 4 s on a single Nvidia RTX A6000 GPU.

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

Bayesian force fields from active learning for simulation of inter-dimensional transformation of stanene

Abstract We present a way to dramatically accelerate Gaussian process models for interatomic force fields based on many-body kernels by mapping both forces and uncertainties onto functions of low-dimensional features. This allows for automated active learning of models combining near-quantum accuracy, built-in uncertainty, and constant cost of evaluation that is comparable to classical analytical models, capable of simulating millions of atoms. Using this approach, we perform large-scale molecular dynamics simulations of the stability of the stanene monolayer. We discover an unusual phase transformation mechanism of 2D stanene, where ripples lead to nucleation of bilayer defects, densification into a disordered multilayer structure, followed by formation of bulk liquid at high temperature or nucleation and growth of the 3D bcc crystal at low temperature. The presented method opens possibilities for rapid development of fast accurate uncertainty-aware models for simulating long-time large-scale dynamics of complex materials.

Chemistry↗
DOE OSTI2020

Electron correlation effects in attosecond photoionization of CO 2

A technique for measuring photoionization time delays with attosecond precision is combined with calculations of photoionization matrix elements to demonstrate how multielectron dynamics affect photoionization time delays in carbon dioxide. Electron correlation is observed to affect the time delays through two mechanisms: autoionization of molecular Rydberg states and accelerated escape from a continuum shape resonance.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2021

Anomalous fast atomic dynamics in bulk metallic glasses

It is believed that the strong glass formers tend to be dense packing and have rather slow atomic dynamics. Here, we report that the anomalous fast atomic dynamics exist in an Au 55 Cu 20 Ag 5 Si 20 metallic glass (MG) with good glass forming ability (GFA) by combining X-ray photon correlation spectroscopy measurements with ab initio molecular dynamics simulations. The fast dynamics are mainly originated from the increased mobility of Cu atoms by Ag addition. Ag atoms are inclined to locate at Cu-centered <0, 2, 8, 1> motifs, which slow or accelerate the dynamics of shell Cu atoms, depending on whether Ag atoms appear at the first shells or not, leading to the formation of structural heterogeneity at the atomic level. The results enrich the understanding of atomic dynamics in Cue-Ag-containing MGs and will trigger more studies on the dynamics-GFA relationship in MGs.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2025

Photoinduced twist and untwist of moiré superlattices

This paper directly observes femtosecond photoinduced twist–untwist motion of WSe₂/MoSe₂ moiré superlattices using ultrafast electron diffraction, showing that optical excitation can dynamically reconfigure moiré lattice structure and interlayer coupling. It supports the proposal by establishing the team’s ability to probe ultrafast interfacial and moiré dynamics, which is relevant to using moiré TMDs as reconfigurable platforms for controlling spin, charge, and molecular interactions.

72 PHYSICS OF ELEMENTARY PARTICLES AND FIELDS↗
DOE OSTI2022

Radiation effects on the structure and alteration behavior of an SiO 2 –Al 2 O 3 –B 2 O 3 –Na 2 O glass

As borosilicate glasses are used in many countries to immobilize fission products and minor actinides after spent fuel reprocessing before storage in a deep geological repository, assessing that their chemical durability is of paramount importance. Here, pristine and preirradiated (952 MeV, 136Xe) SiO 2 –B2O 3 –Al2O 3 –Na 2 O glasses with the same molar ratios as in the French SON68 and ISG glasses have been subjected to aqueous corrosion in deionized water and in silica-saturated solution to measure the initial and longer term alteration rates. Pristine and preirradiated glasses corrode following the same mechanisms, but the preirradiation has a strong impact on the initial dissolution rate (increase by a factor of 5.6), and on the alteration layer depth in silica-saturated conditions (by two- to threefolds). The later result is related to the formation of a more porous, less passivating gel on the preirradiated glass specimen. Using both experimental spectroscopies (NMR, IR, and SFG) and classical molecular dynamics, the radiation effects on the glass structure and water diffusion have been assessed. After preirradiation, the density and the polymerization degree of the glass decrease, whereas the topological disorder increases. In consequence, water diffusion accelerates. These observations allow to correlate the radiation impact on the alteration behavior to the structural changes.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Polymer Chain Diffusion in All-Polymer Nanocomposites: Confinement vs Chain Acceleration

All-polymer nanocomposites, in which soft, polymer-based nanoparticles are dispersed in the polymer matrix, have received great interest lately due to their potential use in a range of applications, including drug delivery and self-healing materials. However, the impact of this new class of nanoparticles on the dynamics of a linear polymer matrix in an all-polymer nanocomposite is still largely unknown. In this work, we report that the addition of polystyrene soft nanoparticles accelerates the diffusion of high molecular weight linear PS chains over a range of nanoparticle loadings. Our results show that at nanoparticle loadings below 1%, the diffusion of the linear matrix increases, presumably via a constraint release mechanism. At loadings above 1%, the increase in diffusion is mitigated by confinement effects of the nanoparticles. Thus, the response of these all-polymer nanocomposites is dominated by the balance of entropic confinement of the chain, which slows diffusion and a constraint release mechanism that speeds up the diffusion. However, the diffusion of the linear chain in the all-polymer nanocomposite is faster than that of the same chain in the melt at all loadings, diverging from the behavior of most nanocomposites with hard, impenetrable nanoparticles. Thus, the mechanism that accelerates the chain diffusion dominates in these systems. This behavior is unusual and fundamentally different than what has been reported for nanocomposites with hard inorganic nanoparticles, indicating that new perspectives are needed for these materials where the control of loading can either accelerate or decelerate the dynamics of the matrix in a distinctive manner.

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

Ultra-Fast High-Precision Metallic Nanoparticle Synthesis using Laser-Accelerated Protons

Laser-driven proton acceleration, as produced during the interaction of a high-intensity (I > 1 × 10 18 W/cm 2 ), short pulse (<1 ps) laser with a solid target, is a prosperous field of endeavor for manifold applications in different domains, including astrophysics, biomedicine and materials science. These emerging applications benefit from the unique features of the laser-accelerated particles such as short duration, intense flux and energy versatility, which allow obtaining unprecedented temperature and pressure conditions. In this paper, we show that laser-driven protons are perfectly suited for producing, in a single sub-ns laser pulse, metallic nanocrystals with tunable diameter ranging from tens to hundreds of nm and very high precision. Our method relies on the intense and very quick proton energy deposition, which induces in a bulk material an explosive boiling and produces nanocrystals that aggregate in a plasma plume composed by atoms detached from the proton-irradiated surface. The properties of the obtained particles depend on the deposited proton energy and on the duration of the thermodynamical process. Suitably controlling the irradiated dose allows fabricating nanocrystals of a specific size with low polydispersity that can easily be isolated in order to obtain a monodisperse nanocrystal solution. Molecular Dynamics simulations confirm our experimental results.

77 NANOSCIENCE AND NANOTECHNOLOGY↗
DOE OSTI2020

Molecular simulation of the separation of toluene and p-xylene with the thermally-robust ionic liquid triphenyl-p-phenyl sulfonyl phenyl phosphonium

In this work, we use molecular dynamics simulations to evaluate the behavior of toluene (TOL) and p-xylene (XYL) within three different solvents: n -hexane (HEX), n -heptane (HEP), and a thermally-robust ionic liquid (triphenyl-p-phenyl sulfonyl phenyl phosphonium + bis(trifluoromethylsulfonyl)imide (TPSP + Tf 2 N)). Several different temperatures are explored, while the solvation structures, energetics, and dynamics of TOL and XYL are evaluated. The solvation free energy (ΔG solv ) and transfer free energy (ΔG transfer ) of TOL and XYL are calculated using thermodynamic integration (TI). The transfer free energy of XYL from HEP/HEX to TPSP + Tf 2 N is more favorable than that of TOL over the entire temperature range, while the transfer of both TOL and XYL becomes less favorable as the temperature increases. The diffusion rates of TOL and XYL are lower in the ionic liquid when compared to HEP and HEX, but the thermal stability of TPSP + Tf 2 N may allow for higher operating temperatures and accelerated mass transport rates.

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

The Future of a Myriad of Accelerated Biodiscoveries Lies in AI‐Powered Mass Spectrometry and Multiomics Integration

The intersection of modern artificial intelligence (AI) and mass spectrometry (MS) is set to transform the MS‐based “omics” research fields, particularly proteomics, metabolomics, lipidomics, and glycomics, enabling advancements across a wide range of domains, from health to environment and industrial biotechnology. Beginning with an overview of key challenges inherent in MS software pipelines, this personal perspective explores how AI‐driven solutions can address them to enhance data processing, integration and interpretation. It proposes a paradigm shift in molecular identification and quantitation algorithms, leveraging AI to enable holistic interpretation of MS‐based multiomics data. While centered on MS‐based omics, this holistic AI‐driven paradigm is also critical for connecting dynamic biochemical changes to genomics and transcriptomics contexts, reinforcing the integrative value of MS in multiomics research. Ultimately, this AI‐driven approach could enhance efficiency, accuracy, and molecular breadth of coverage, deepening our systems‐level understanding of biological processes and accelerating a myriad of biodiscoveries.

47 OTHER INSTRUMENTATION↗
DOE OSTI2022

Study of solid molecular deuterium D 2 growth under gas pressure

The injection of high-speed cryogenic pellets made of frozen hydrogen-isotopes, represents to date the most effective method to fuel magnetically confined thermonuclear fusion plasmas. Additionally, the injection of very large pellets composed of cryogenic solid of some suitable impurity (typically a noble-gas such as H 2 , Ne, or H 2 /Ne, D 2 /Ne mixtures), shattered in relatively small fragments just before entering the plasma, seems to be the most promising method to reduce the damage risks for the plasma-facing components in case of a plasma disruption. This technology, known as "Shattered Pellet Injection" (SPI), allows to spread out the plasma energy and mitigate possible damage to the in-vessel components, as well as to densify the plasma to suppress the formation of runaway electrons, and/or dissipate their energy. Several techniques to produce and launch cryogenic pellets have been investigated in the past decades. "Pipe gun" injectors are reliable and relatively simple devices are still commonly used today. They make use of single- or two-stage pneumatic light-gas guns to accelerate the pellet at high speeds. In these injectors, the cryogenic pellets are formed “ in situ ” (i.e., inside the launching barrel), by de-sublimating them directly from the gas phase, i.e., at temperatures and pressures below those of the triple point. The simplest case is pure deuterium pellets ($\mathcal{T}$ < 18.7 K, $\mathcal{P}$ < 171.3 hPa). The production of good quality solid deuterium, capable of withstanding the mechanical stress during the acceleration of the pellets, is a key issue. To this end the phase transition of deuterium from gas to solid (and vice versa) is modeled with extensive molecular-dynamics (MD) simulations. Furthermore, the solid growth from the gas phase is simulated in an ample range of temperatures and pressures, to find the best compromise between growth velocity and mechanical properties of the resulting solid system.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Ultrafast Events in Electrocyclic Ring-Opening Reactions

Electrocyclic reactions are characterized by the concerted formation and cleavage of multiple σ and π bonds in a molecular system and have been extensively studied since they were introduced by Robert Burns Woodward and Roald Hoffmann in 1965. Recent advances and the integration of time-resolved experiments and nonadiabatic quantum molecular dynamics simulations have transformed the traditional understanding of electrocyclic reactions beyond the Woodward–Hoffmann rules. In this review, we focus on recent studies of 1,3-cyclohexadiene and two of its derivatives, α-phellandrene and α-terpinene, to shed light on the underlying mechanisms of electrocyclic photochemical reactions. We highlight recent progress in ultrafast electron diffraction techniques and the simulation approach of ab initio multiple spawning. Together, these approaches can elucidate molecular structure dynamics from femtosecond to picosecond timescales as well as nuclear and electronic responses at conical intersections.

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

Unraveling disadhesion mechanism of epoxy/CSH interface under aggressive conditions

Almost all FRP-reinforced concrete exposed to the marine environment are susceptible to durability deficiencies due to failure of the FRP/epoxy/concrete interface, especially the vulnerable epoxy/concrete interface, which is frequently attacked by water and ions. Herein, molecular dynamics simulations of epoxy/CSH interface under different environmental exposures are utilized to specifically elucidate the degradation principles of water and ions on the bonding properties and mechanical responses of the epoxy/concrete interface. The environmental sequences associated with the interfacial bonding degradation are arranged as Na{sub 2}SO{sub 4} + NaCl>NaCl>Na{sub 2}SO{sub 4} > Water>Dry. Remarkably, water molecules weaken the interaction between the epoxy and CSH by robbing interfacial H- and Ca-O bonds, thereby reducing the energy barrier to interfacial peel and shear failure. The weakening effects are exacerbated by the presence of aggressive ions, as Na{sup +} ions attracted to the CSH surface accumulate more water molecules in the interfacial region by forming hydrated ions clusters, which ultimately accelerates the epoxy/CSH interfacial degradation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computationally aided, entropy-driven synthesis of highly efficient and durable multi-elemental alloy catalysts

Multi-elemental alloy nanoparticles (MEA-NPs) hold great promise for catalyst discovery in a virtually unlimited compositional space. However, rational and controllable synthesize of these intrinsically complex structures remains a challenge. Here, we report the computationally aided, entropy-driven design and synthesis of highly efficient and durable catalyst MEA-NPs. The computational strategy includes prescreening of millions of compositions, prediction of alloy formation by density functional theory calculations, and examination of structural stability by a hybrid Monte Carlo and molecular dynamics method. Selected compositions can be efficiently and rapidly synthesized at high temperature (e.g., 1500 K, 0.5 s) with excellent thermal stability. We applied these MEA-NPs for catalytic NH 3 decomposition and observed outstanding performance due to the synergistic effect of multi-elemental mixing, their small size, and the alloy phase. We anticipate that the computationally aided rational design and rapid synthesis of MEA-NPs are broadly applicable for various catalytic reactions and will accelerate material discovery.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Pb dimerization greatly accelerates charge losses in MAPbI 3 : Time-domain ab initio analysis

Metal halide perovskites constitute a new type of semiconducting materials with long charge carrier lifetimes and efficient light-harvesting. The performance of perovskite solar cells and related devices is limited by nonradiative charge and energy losses, facilitated by defects. Combining nonadiabatic molecular dynamics and time-domain density functional theory, we demonstrate that charge losses depend strongly on the defect chemical state. By considering an extra Pb atom in CH 3 NH 3 PbI 3 , which is a common defect in lead halide perovskites, we investigate its influence on charge trapping and recombination. In a chemically inert form as a Pb interstitial, the extra Pb atom has only a mild influence on charge recombination. However, if the extra Pb atom binds to a native Pb atom to form a dimer, the charge trapping and recombination are greatly accelerated because the Pb-dimer creates a localized midgap trap state that couples strongly to the perovskite valence band edge. Holes disappear from the valence band two orders of magnitude faster than in the pristine perovskite and recombine with conduction band electrons one order of magnitude faster. The simulations identify the phonon modes involved in the nonradiative charge trapping and recombination and highlight the importance of rapid decoherence within the electronic subsystem for long carrier lifetimes. The detailed atomistic analysis of the charge trapping and recombination mechanisms enriches the understanding of defect properties and provides theoretical guidance for improving perovskite performance.

37 INORGANIC, ORGANIC, PHYSICAL, AND ANALYTICAL CH↗