Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Energetic materials”

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

Semiclassical model for calculating exciton and polaron pair energetics at interfaces

Exciton and polaron pair dissociation is a key functional aspect of photovoltaic devices. To improve upon the current state of interfacial transport models, we augment the existing classical models of dielectric interfaces by incorporating results from ab initio calculations, allowing us to calculate exciton and polaron binding energies more accurately. We demonstrate the predictive capabilities of this new model using two interfaces: (i) the boron subphthalocyanine chloride (SubPc) and C60 interface, which is an archetype for many organic photovoltaic devices; and (ii) pentacene and silicon (1 0 0), which represents a hybrid between organic and inorganic semiconductors. Overall, our calculations predict that the insertion of molecular dipoles at interfaces can be used for improving polaron pair dissociation and that sharp transitions in dielectric permittivity can have a stronger effect on the polaron pair dissociation than even the electron-hole Coulomb interaction.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Oxidative Deboronation of Boronic Acids by Hydrogen Peroxide in Planta Generates Borate for Cross-Linking of Rhamnogalacturonan II

Vascular plants require boron to cross-link the rhamnogalacturonan-II (RG-II) domain of pectin to form functional cell walls. Boronic acids, which form reversible esters with cis-diols like borate, have been proposed to influence RG-II cross-linking, though the mechanism remains unclear. We used suspension-cultured rose cells adapted to grow without boron to investigate the effect of boronic acids on RG-II dimerization. When grown with phenylboronic acid (PBA) as the sole boron source, nearly all RG-II was crosslinked, whereas methylboronic acid (MBA) only partially restored cross-linking. In contrast, in vitro assays showed that homogeneous RG-II monomers did not dimerize with alkyl or aryl boronic acids unless supplemented with hydrogen peroxide (H2O2), which oxidatively converts boronic acids to boric acid. Real-time NMR spectroscopy and density functional theory calculations provided insight into the reaction mechanism and energetics of oxidation respectively. Together, our data show that exogenous boronic acids are a source of boric acid for plants, and that the deboronation reaction generates aryl or alkyl alcohol byproducts that can undergo further chemical modification in planta. The fate and potential roles of these byproducts in planta remain to be determined.

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

Influence of nanoscale interfaces on the dynamic deformation and spall failure of Cu–Fe alloy microstructures

Additively manufactured immiscible Cu–Fe alloys can exhibit a distribution of nanoscale interfaces due to the distribution of nanoscale clusters in the equilibrium and metastable phases. Molecular dynamics simulations investigate the role of such interfaces on the phase stability and transformation behavior during shock compression, as well as the mechanisms of damage nucleation during spall failure. The model multiphase Cu–Fe systems studied here comprise a distribution of Fe clusters in an FCC Cu matrix, as well as Cu clusters in a BCC Fe matrix. The length scales of the nanoscale clusters determine the energetics of the interfaces that can result in FCC/BCC (equilibrium) or FCC/FCC (metastable) interfaces in the Cu system matrix, and BCC/FCC (equilibrium) or BCC/BCC (metastable) interfaces in the Fe matrix. The MD simulations demonstrate that nanoscale metastable interface microstructures can induce plastic deformation and also stabilize phases or suppress phase transformations in metastable phase clusters. In contrast, equilibrium interfaces can influence phase transformation thresholds and serve as additional void nucleation sites during failure.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Investigation of the distribution of laser-induced damage precursors in etched grating sidewalls fabricated using inductively coupled plasma reactive-ion etching

Inductively coupled plasma reactive-ion etching for grating fabrication is explored with emphasis on the type and distribution of process-induced defects that can act as precursors to laser-induced damage initiation. To provide flexibility in the ability to perform a variety of material characterization methods, grating-like features with large pitch were fabricated in silica monolayers. The results indicate that the etch process introduces surface fluorine and subsurface carbon-containing contamination, which reduces the laser-induced damage threshold. The local relative concentrations of the laser-damage precursors introduced during the etch step of the manufacturing process were investigated using angle-resolved damage testing, taking advantage of the shift of the location of the peak electric-field intensity incident on the grating sidewall. The results suggest an increased level of etch-related defects near the base of the sidewall, with directly fluorinated silicon identified as the primary defect.

Dent, R. [Univ. of Rochester, NY (United States). ↗
DOE OSTI2026

Unraveling interphase-driven failure pathways in LiMn0.6Fe0.4PO4/graphite pouch cells

LiMnxFe1−xPO4 (LMFP) is a promising high-voltage, thermally stable, and earth-abundant cathode material, yet its practical application is limited by interphase instability and Mn dissolution. In this work, we systematically evaluate LiMn0.6Fe0.4PO4/graphite pouch cells using three electrolyte formulations including control carbonate electrolyte, control + 2 wt% vinylene carbonate (VC), and control + 2 wt% VC + 1 wt% 1,3,2-dioxathiolane 2,2-dioxide (DTD), to establish how electrolyte composition governs interphase chemistry and long-term degradation. Electrochemical testing shows that both additives are preferentially reduced prior to ethylene carbonate (EC) during cell formation, generating robust cathode-electrolyte interphase (CEI) and solid-electrolyte interphase (SEI) layers that suppress gas evolution and raise the first-cycle coulombic efficiency to 89.3%. Additionally, the dual-additive electrolyte delivers the most stable performance, retaining over 85% capacity after 600 cycles while minimizing impedance growth under long-term cycling at C/3 and 40 °C. Soft X-ray absorption spectroscopy confirms that VC + DTD effectively suppresses electrolyte oxidation at the cathode surface, and micro-X-ray fluorescence shows substantially reduced Mn dissolution and deposition on the graphite anode. Density functional theory simulations further provided insights into the structural and energetic influences of alkoxide species on the cathode surface, proposing a Mn2+ extraction mechanism. The combined experimental and computational findings establish a mechanistic link between electrolyte composition and interphase evolution, highlighting the effectiveness of electrolyte engineering for extending the operational lifetime of LMFP-based lithium-ion batteries.

Chak, Chanmonirath Michael↗
DOE OSTI2022

Neutron vibrational spectroscopic evidence for short H∙∙∙H contacts in the RNiInH 1.4; 1.6 (R = Ce, La) metal hydride

Intermetallic metal hydrides are critical materials for hydrogen storage applications, however, metal hydrides with greater storage capacities are still needed. Within metal hydrides, the volumetric storage capacities are limited by the number of hydrogen-accommodating interstitial sites which can be simultaneously occupied given a minimum hydride nearest-neighbor distance of ~2.1 A, according to the Switendick-Westlake criterion. To date, violations of this criterion are rare. Perhaps the most well studied compounds violating this criterion are the RNiInHx compounds (R = Ce, La, Nd). Previous neutron diffraction studies on the deuterated species revealed the presence of Ni-D∙∙∙D-Ni-D∙∙∙D-Ni chains with anomalously close D∙∙∙D contacts of ~1.6 A. Yet there are no neutron vibrational spectroscopic investigations reported for these atypical hydrides. Here we use neutron vibrational spectroscopy (NVS) measurements to probe the hydrogen dynamics in LaNiInHx (x = 0.67, 1.6) and CeNiInH 1.4 . For x > 0.67, the presence of close H...H contacts yields two related features in the vibrational spectrum centered near ~90 meV corresponding to the oscillations of paired H atoms simultaneously occupying neighboring R 3 Ni tetrahedra. Notably, these features are energetically distinct from comparable vibrational motions for "unpaired" H atoms when x = 0.67. To compare, we also present powder neutron diffraction and NVS measurements for the newly characterized, chemically similar Sn compounds CeNiSnH, CeNiSnH 2 , and CeNiSnD 2 . These compounds also contain R 3 Ni tetrahedra, however, the H-occupied tetrahedra are well separated from each other with the closest H∙∙∙H distances exceeding 2.1 A, and the Switendick-Westlake criterion is not violated. Consequently, the spectral signature of the close H∙∙∙H contacts is absent in these hydrides.

08 HYDROGEN↗
DOE OSTI2022

Interfacial Unit-Dependent Catalytic Activity for CO Oxidation over Cerium Oxysulfate Cluster Assemblies

Atomically precise cerium oxo clusters offer a platform to investigate structure–property relationships that are much more complex in the ill-defined bulk material cerium dioxide. We investigated the activity of the MCe 70 torus family (M = Cd, Ce, Co, Cu, Fe, Ni, and Zn), a family of discrete oxysulfate-based Ce 70 rings linked by monomeric cation units, for CO oxidation. CuCe 70 emerged as the best performing MCe 70 catalyst among those tested, prompting our exploration of the role of the interfacial unit on catalytic activity. Temperature-programmed reduction (TPR) studies of the catalysts indicated a lower temperature reduction in CuCe 70 as compared to CeCe 70 . In situ diffuse reflectance infrared Fourier transform spectroscopy (DRIFTS) indicated that CuCe 70 exhibited a faster formation of Ce 3+ and contained CO bridging sites absent in CeCe 70 . Isothermal CO adsorption measurements demonstrated a greater uptake of CO by CuCe 70 as compared to CeCe 70 . The calculated energies for the formation of a single oxygen defect in the structure significantly decreased with the presence of Cu at the linkage site as opposed to Ce. Furthermore, this study revealed that atomic-level changes in the interfacial unit can change the reducibility, CO binding/uptake, and oxygen vacancy defect formation energetics in the MCe 70 family to thus tune their catalytic activity.

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

Defect modeling in semiconductors: the role of first principles simulations and machine learning

Abstract Point defects in semiconductors dictate their electronic and optical properties. Vacancies, interstitials, substitutional defects, and defect complexes can form in the semiconductor lattice and significantly impact its performance in applications such as solar absorption, light emission, electronics, and catalysis. Understanding the nature and energetics of point defects is essential for the design and optimization of next-generation semiconductor technologies. Here, we provide a comprehensive overview of the current state of research on point defects in semiconductors, focusing on the application of density functional theory (DFT) and machine learning (ML) in accelerating the prediction and understanding of defect properties. DFT has been instrumental in accurately calculating defect formation energies, charge transition levels, and other defect-related properties such as carrier recombination rates and lifetimes, and ion migration barriers. ML techniques, particularly neural networks, have emerged as powerful tools for enabling rapid prediction of defect properties at DFT-accuracy in order to overcome the expense of using large supercells and advanced functionals. We begin this article with a discussion of different types of point defects and complexes, their impact on semiconductor properties, and the experimental and DFT approaches typically used for their characterization. Through multiple case studies, we explore how DFT has been successfully applied to understand defect behavior across a variety of semiconductors, and how ML approaches integrated with DFT can efficiently predict defect properties and facilitate the discovery of new materials with tailored defect behavior. Overall, the advent of ‘DFT+ML’ promises to drive advancements in semiconductor technology, catalysis, and renewable energy applications, paving the way for the development of high-performance semiconductors which are defect-tolerant or have desirable dopability.

Rahman, Md Habibur (ORCID:000000027705984X)↗
DOE OSTI2026

A cell-centered AMR-ALE framework for 3D multi-material hydrodynamics. Part I: Lagrangian and indirect Euler AMR algorithms

Many applications of physics and engineering involve wide ranges of time and spatial scales. The numerical simulation of localized small scales such as shock waves and material interfaces requires a large number of computational cells in these regions. For these applications, Lagrangian and Arbitrary-Lagrangian-Eulerian (ALE) related methods are engaging since the moving mesh feature naturally brings mesh cells on shock discontinuities and material interfaces are carefully captured. In addition, Adaptive-Mesh-Refinement (AMR) strategies aim to optimize computational resources by concentrating finer mesh cells only in areas of interest while using coarser cells elsewhere. A key but challenging AMR requirement consists in efficiently distributing the computational effort to achieve high accuracy without the prohibitive computational costs associated with uniformly fine grids. Here, in this document, the coupling of the p4est AMR library with a cell-centered Lagrangian scheme is presented with the goal to perform reliable 3D Lagrangian-AMR and indirect Euler-AMR multi-material simulations. In particular, it is shown that starting from a 3D indirect ALE code, the memory management and load balancing requirements can be delegated to an external library (here the p4est library) to unlock ALE-AMR capabilities. First, we present a strategy to transcribe the octant-based connectivity of the 3D AMR framework with that of an unstructured mesh of polygonal cells used in Lagrangian hydrodynamics. Then, we show how refinement and coarsening operations must be adapted to the particular Lagrangian framework to ensure the conservation of volume during those steps. Finally, several numerical test cases are presented that demonstrate the capabilities of the Lagrangian-AMR and indirect Euler-AMR algorithms.

3D cell-centered Lagrangian numerical scheme↗
DOE OSTI2025

Magnetic Phase Diagram of 2D Magnet Fe 5– x GeTe 2 Revealed by Cryogenic Lorentz 4D-STEM

A two-dimensional van der Waals ferromagnet Fe 5–x GeTe 2 (F5GT) exhibits room-temperature magnetic transition and magnetic anisotropy. Studies have identified diverse magnetic states, including stripe domains, skyrmionic type-I and topologically trivial type-II bubbles, depending on the magnetic and thermal history. Yet, the underlying micromagnetic energetics driving these states remains unclear. Here, we establish the magnetic phase diagram of bulk F5GT using cryogenic Lorentz four-dimensional scanning transmission electron microscopy. We reveal that stripe domains spontaneously form upon zero-field cooling, while metastable magnetic bubbles are stabilized by external fields. Transitions from type-I to type-II bubbles are driven by an oblique external field. Micromagnetic simulations confirm that these transitions arise from the interplay of uniaxial anisotropy, dipolar interactions, and external fields, without requiring significant Dzyaloshinskii–Moriya Interaction (DMI). In conclusion, these findings clarify the micromagnetic origin of spin textures of F5GT and establish it as a tunable platform for nanoscale topological magnetism.

2D Magnets↗
DOE OSTI2022

First-Principles Elucidation of Initial Dehydrogenation Pathways in Mg(BH 4 ) 2

Complex borohydrides such as Mg(BH 4 ) 2 offer one of highest capacities to chemically store hydrogen for onboard applications; however, it suffers greatly from kinetic constraints that prevent realization of full capacity and reversibility. Understanding these kinetic limitations solely from experiments is extremely challenging due to the unusual complexity of various competing elemental reaction steps involved during the de/rehydrogenation reaction. This work aims to map out the energetics associated with initial dehydrogenation of Mg(BH 4 ) 2 from first-principles simulations and to identify the preferred reaction pathways. Our calculations suggest the rate-limiting step during BH 4 – –B 3 H 8 – conversion is the formation of the B 2 H 7 – intermediate. We further emphasize and clarify that the B 3 H 8 – and H – intermediates, formed during initial Mg(BH 4 ) 2 decomposition, appear as molecular species that are embedded in the Mg–BH 4 –Mg matrix as evidenced in the nuclear magnetic resonance measurements and not as bulk MgH 2 and Mg(B 3 H 8 ) 2 as previously assumed in theoretical predictions of the thermodynamics.

08 HYDROGEN↗
DOE OSTI2022

Data-driven design of novel halide perovskite alloys

The great tunability of the properties of halide perovskites presents new opportunities for optoelectronic applications as well as significant challenges associated with exploring combinatorial chemical spaces. Here, in this work, we develop a framework powered by high-throughput computations and machine learning for the design and prediction of mixed cation halide perovskite alloys. In a chemical space of ABX 3 perovskites with a selected set of options for A, B, and X species, pseudo-cubic structures with B-site mixing are simulated using density functional theory (DFT) and several properties are computed, including stability, lattice constant, band gap, vacancy formation energy, refractive index, and optical absorption spectrum, using both semi-local and hybrid functionals. Neural networks (NN) are used to train predictive models for every property using tabulated elemental properties of A, B, and X site species as descriptors. Starting from a DFT dataset of 229 points, we use the trained NN models to predict the structural, energetic, electronic, and optical properties of an extensive dataset of 17955 compounds, and perform high-throughput screening in terms of stability, band gap, and defect tolerance, to obtain 392 promising compounds that are ranked as potential absorbers according to their photovoltaic figure of merit. Compositional trends in the screened set of attractive mixed cation halide perovskites are revealed and additional computations are performed on selected compounds. The data-driven design framework developed here is promising for designing novel mixed compositions and can be extended to a wider perovskite chemical space in terms of A, B, and X species, different kinds of mixing at the A, B, or X sites, non-cubic phases, and other properties of interest.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Defects in epitaxial Ru(0001) on Al2O3(0001): Dislocations, stacking faults, and deformation twins

In this study, defects in epitaxial Ru(0001) films on c-plane sapphire, with nominal thicknesses of 10–80 nm, deposited at 350 °C and step-annealed to 950 °C, were characterized using transmission electron microscopy. The variation of Ru and sapphire lattice parameters with temperature is such that the misfit strain for the observed 30° rotated-honeycomb epitaxial relationship is essentially constant with temperature at 1.5%, resulting in a biaxial stress of 10.0 GPa and an energy density of 150 MJ m –3 in unrelaxed films. Stress relaxation occurs by the formation of defects. For the 20–80 nm thick films, the defects are a- and c-type dislocations and stacking faults, argued to be of I 2 type. In addition, the films show the surprising presence of $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ deformation twins. The 10 nm-thick films were found to be defect free. The critical thickness for misfit strain relaxation via the formation of threading and misfit dislocations is computed as 7±2 nm, depending on the choice of the dislocation core radius. Energetic analysis of twin formation, using both the infinite-matrix and the finite-matrix (Mori–Tanaka) approaches, provides values of the twin aspect ratios, assumed to be ellipsoidal, and shows that the latter but not the former approach can qualitatively explain the formation of the observed twins. In addition to providing the maximum strain relief compared to other potential twin types, $\{11\bar{2}1\}1/3\langle11\bar{2}\bar{6}\rangle$ twins do not require lattice shuffles and have a boundary that is a special boundary, namely, a 35° tilt boundary with a-type dislocations every other {0002} plane, that may also favor their formation.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Evidence for anti-synergism between ion-assisted etching and in-plasma photoassisted etching of silicon in a high-density chlorine plasma

Etching of p-Si in 60 mTorr 10%Cl2/90%Ar Faraday-shielded inductively coupled high density plasmas was investigated under both ionassisted etching (IAE) and photoassisted etching (PAE) conditions. Real-time etching rates and after-etching Si surface chemical compositions were obtained by laser interferometry and vacuum-transfer x-ray photoelectron spectroscopy (XPS), respectively. Precisely controlled ion energy distributions (IEDs) were generated by applying pulsed negative DC bias on the conductive sample stage. Above a 36 eV threshold at a total flow rate of 250 SCCM, the IAE rate increased with the square root of the ion energy. In contrast to the DC bias, etching under RF bias did not exhibit a threshold ion energy because of the wide IED. XPS spectra revealed that the surface layer under PAE conditions had a significantly lower chlorine content, composed of only SiCl. Under IAE conditions, however, silicon dangling bonds (Si•), SiCl 2 , and SiCl 3 were found on the surface, in addition to SiCl, with a relative abundance of SiCl > SiCl 2 , and SiCl 3 . The absence of higher chlorides and Si• under PAE conditions suggested that vacuum ultraviolet photons and above threshold-energy ions interact with the surface very differently. By varying the duty cycle of the pulsed DC bias, it was found that the IAE rate scaled with the energetic ion dose, but only for low duty cycles. For higher duty cycles, the apparent IAE yield fell off with an increasing Cl coverage on the surface, as the duty cycle went up, which pointed to a negative synergy (antisynergism) between PAE and IAE as the explanation. This antisynergism was further supported by the observed decrease of the total etching rate with an increasing period of the pulsed DC bias. A plausible mechanism is that increasing the pulsing period causes more near-surface damage, creating more recombination centers that lead to a higher loss rate of electron-hole pairs through recombination, thereby reducing the PAE rate.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2026

First-principles DFT modeling of nitrobenzene adsorption on the Ag(111) surface at varying monolayer coverages

Here, the adsorption behavior of nitrobenzene on the Ag(111) surface as a function of coverage is investigated using density functional theory. Adsorption energies and optimized geometries are analyzed together with isolated intermolecular interaction calculations and electronic structure analysis, including Bader charge partitioning and projected density of states, to disentangle the roles of adsorbate–surface bonding and through-space adsorbate–adsorbate interactions. At low and intermediate coverages (θ = 1/9 and θ = 2/9), bidentate adsorption at top sites is favored due to strong adsorbate–surface interactions, with additional stabilization at θ = 2/9 arising from favorable intermolecular separations. At higher coverage (θ = 1/3), bidentate adsorption is destabilized by strong intermolecular repulsion, and monodentate adsorption becomes energetically preferred as rotational freedom allows more favorable intermolecular spacing. Charge density differences, Bader charge, and density of states analyses show that charge transfer from the Ag surface is localized primarily on the nitro group and increases with coverage and adsorption denticity, although this increase does not directly correlate with adsorption strength at high coverage due to competing intermolecular interactions. These results demonstrate that surface coverage can induce a transition in preferred adsorption denticity driven by intermolecular interactions, highlighting the importance of adsorbate packing in organic molecule adsorption at metal interfaces.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Density Functional Theory Study of Iron–Oxygen Divacancies in Magnetite (Fe 3 O 4 ) and Hematite (Fe 2 O 3 )

Density functional theory (DFT) calculations are employed to investigate the formation energies, charge redistribution, and binding energies of iron–oxygen divacancies in magnetite (Fe 3 O 4 ) and hematite (Fe 2 O 3 ). For magnetite, we focus on the low-temperature phase to explore variations with local environments. Building on previous DFT calculations of the variations in formation energies for oxygen vacancies with local charge and spin order in magnetite, we extend this analysis to include octahedral iron vacancies before analyzing the iron–oxygen divacancies. We also assessed the relative stability of iron–oxygen divacancies by comparing their formation energies with those of individual vacancies. Our findings reveal a significant energetic driving force for the formation of divacancy clusters, particularly in magnetite, where divacancies in the +1 charge state exhibit formation energies comparable to those of neutral iron vacancies under oxidizing conditions. In hematite, the results indicate a strong tendency for oxygen vacancies to bind to iron vacancies. These results highlight the significance of iron–oxygen vacancy complexes in the transport properties of iron oxides, with particular relevance to diffusion mechanisms under irradiation conditions.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Designing a quantum-accurate machine-learning potential to enable large-scale simulations of deuterium under shock

Large-scale molecular dynamics of deuterium under shock can elucidate kinetic processes vital to the target design in inertial confinement fusion and high-energy-density experiments. However, modeling the complex evolution of this material from an insulating molecular state at ambient pressure to an ionized, atomic fluid under strong shock is beyond the capability of simple pair and even bond order potentials. We thus train a quantum-accurate and broadly transferable machine-learning interatomic potential for deuterium using the Chebyshev Interaction Model for Efficient Simulations framework. We show that due to an improved description of the molecular-to-atomic transition, our model is able to better reproduce the ab initio equation of state, radial distribution functions, and principal Hugoniot than bond order potentials. This represents an important step toward large-scale quantum-accurate and nonequilibrium simulations of complicated systems under dynamic changes including phase transitions.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗