Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “Classical 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 703 records · Page 39
All sourcesNASA NTRSDOE OSTI
DOE OSTI2022

UV-induced dissociation of CH 2 BrI probed by intense femtosecond XUV pulses

The ultraviolet (UV)-induced dissociation and photofragmentation of gas-phase CH 2 BrI molecules induced by intense femtosecond extreme ultraviolet (XUV) pulses at three different photon energies are studied by multi-mass ion imaging. Using a UV-pump–XUV-probe scheme, charge transfer between highly charged iodine ions and neutral CH 2 Br radicals produced by C–I bond cleavage is investigated. In earlier charge-transfer studies, the center of mass of the molecules was located along the axis of the bond cleaved by the pump pulse. In the present case of CH 2 BrI, this is not the case, thus inducing a rotation of the fragment. We discuss the influence of the rotation on the charge transfer process using a classical over-the-barrier model. Our modeling suggests that, despite the fact that the dissociation is slower due to the rotational excitation, the critical interatomic distance for charge transfer is reached faster. Furthermore, we suggest that charge transfer during molecular fragmentation may be modulated in a complex way.

74 ATOMIC AND MOLECULAR PHYSICS↗
DOE OSTI2022

Persistence of negative vacancy and self-interstitial formation energies in atomistic models of amorphous silicon

We use classical potential atomistic modeling to demonstrate that negative vacancy and selfinterstitial formation energies are inherently present in amorphous silicon (a-Si). Our conclusions are based on comprehensive calculations of vacancy and self-interstitial formation energies in nine different models, ranging in size from 1000 to 2744 atoms and formed by quenching liquid Si at rates between 10 10 and 10 12 K/s. We account for a range of possible reference energies corresponding to different choices of reservoirs for atoms inserted into or removed from our models. Moreover, we show that successive addition of either vacancies or self-interstitials does not deplete the number of negative formation energy sites. Instead, it further lowers defect formation energies. Furthermore, these findings lead us to conclude that a-Si is not a structure at or near a local enthalpy minimum nor can it be brought to a local enthalpy minimum through the introduction of constitutional defects.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Modeling the Electronic Absorption Spectra of the Indocarbocyanine Cy3

Accurate modeling of optical spectra requires careful treatment of the molecular structures and vibronic, environmental, and thermal contributions. The accuracy of the computational methods used to simulate absorption spectra is limited by their ability to account for all the factors that affect the spectral shapes and energetics. The ensemble-based approaches are widely used to model the absorption spectra of molecules in the condensed-phase, and their performance is system dependent. The Franck–Condon approach is suitable for simulating high resolution spectra of rigid systems, and its accuracy is limited mainly by the harmonic approximation. In this work, the absorption spectrum of the widely used cyanine Cy3 is simulated using the ensemble approach via classical and quantum sampling, as well as, the Franck–Condon approach. The factors limiting the ensemble approaches, including the sampling and force field effects, are tested, while the vertical and adiabatic harmonic approximations of the Franck–Condon approach are also systematically examined. Our results show that all the vertical methods, including the ensemble approach, are not suitable to model the absorption spectrum of Cy3, and recommend the adiabatic methods as suitable approaches for the modeling of spectra with strong vibronic contributions. We find that the thermal effects, the low frequency modes, and the simultaneous vibrational excitations have prominent contributions to the Cy3 spectrum. The inclusion of the solvent stabilizes the energetics significantly, while its negligible effect on the spectral shapes aligns well with the experimental observations.

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

ORNL_AISD_NiPt_108atoms

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for crystal structures with 108 atoms with 1,900 configurations. The data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system 'MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001' [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in a compressed zipped folders atoms108.zip. The zipped folder contains the data structured in the following way: - Ni_ground_state.cfg --> atomic configuration for the pure nickel - Pt_ground_state.cfg --> atomic configuration for the pure platinum - Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder - Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms108.zip is 66,132. This dataset is an extension to the dataset ORNL_AISD_NiPt [5] that has been previously released with crystal structures of 256 atoms, 864 atoms, and 2,048 atoms, with the same methodology for data collection. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D [5] M. Karabin, M. Lupo Pasini, and M. Eisenbach. ORNL_AISD_NiPt. United States: N. p., 2023. Web. doi:10.13139/OLCF/1958172.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Intermediate time sub-diffusion and stress relaxation in ring polymer melts

The slow dynamics of non-concatenated ring melts remains a frontier problem in polymer science with implications for many soft material environments including cellular biophysics. Here, in this work, we report large-scale simulations of model ring melts that analyze the monomer and center-of-mass (CM) mean square displacements (MSD) and stress relaxation function on intermediate time and length scales. The degree of dynamical slowing down is characterized by the maximally sub-diffusive fractional time scaling exponents. The data span an exceptionally wide range of ring degrees of polymerization and stiffnesses and are not successfully organized based on the classic measure linear chain entanglement, N/N e . Rather, we find that the crossover degree of polymerization, N D , based on ring macromolecular caging that successfully allows master curves to be constructed for the long-time CM self-diffusion constant also collapses these temporal dynamic scaling exponents. Different properties display different exponents and exhibit one or two regimes of linear variation with the logarithm of N D / N . A distinct crossover of the CM-MSD and stress relaxation exponents emerges at sufficiently large N or stiffness that is not found for the monomer MSD, indicating a novel form of dynamic decoupling. This crossover aligns with the predicted critical degree of polymerization for transitioning from a weak to strong caging regime, indicative of activated transport. The latter may reflect the emergence of an intermolecular collective contribution to stress in analogy with dense soft colloidal matter. Suggestions are made for future theoretical work to address the rich patterns of behavior discovered.

Anomalous diffusion↗
DOE OSTI2023

ORNL_AISD_NiPt

This dataset describes the nickel-platinum (NiPt) solid solution binary alloy, where the two constituent elements nickel (Ni) and platinum (Pt) are randomly placed on the face centered cubic (FCC) crystal structure, with the lattice constant of 3.840 angstroms. The dataset comprises data for three different sizes of the crystal structure: 256 atoms, 864 atoms, and 2,048 atoms, each of which contains 1900 configurations. For each size of the crystal structure, the data set was generated for concentrations ranging from 0at% of Pt to 100at% of Pt in the NiPt binary system, with increasing the concentration of Pt in the system every 5at%. For each one of the chemical compositions, 100 random configurations were generated, each with a different random seed. Each of the output files contains the mass, type, atomic coordinates, energy per atom, and forces in x, y, and z directions respectively. For each atomic configuration, the output was collected every 150 steps during the minimization stage and every 1000 steps during the replica exchange stage. Large-scale Atomic/Molecular Massively Parallel Simulator (LAMMPS) [1], which is a molecular dynamics code, was used to generate data for NiPt alloy. The simulation used the interatomic potential for NiPt binary system MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [3] from the OpenKIM library (Open Knowledgebase of Interatomic Models) [2]. This potential was developed based on the second nearest-neighbor modified embedded-atom method (2NN MEAM). The simulation process begins with the generation of the random NiPt structure and follows with the short minimization and replica exchange simulation. The minimization procedure adjusts atomic coordinates and performs energy minimization, which typically leads to a local potential energy minimum. The method used for the minimization was the conjugate gradient algorithm. A short replica exchange (parallel tempering) simulation involves four replicas (ensembles) of a system and follows the minimization stage. Multiple snapshots of the configuration were collected during the minimization and replica exchange stages. NiPt alloy is interesting due to its magnetic and charge transfer properties [4]. The data is provided in three compressed zipped folders: atoms256.zip, atoms864.zip, atoms2048.zip Each zipped folder contains the data that describes crystals of size 256 atoms, 864 atoms, and 2,048 atoms respectively. Each one of the three zipped folders contains the data structured in the following way: -Ni_ground_state.cfg --> atomic configuration for the pure nickel -Pt_ground_state.cfg --> atomic configuration for the pure platinum -Pt#_filtered --> folders containing atomic configurations for #at% concentration of platinum. The folder contains 100 atomic configurations, each saved in a subfolder. Each subfolder named config* is associated with a specific atomic configuration. Each of these subfolders contains files with .cfg format, corresponding to outputs for each atomic configuration The total number of atomic configurations contained in atoms256.zip is 65,046. The total number of atomic configurations contained in atoms864.zip is 63,936. The total number of atomic configurations contained in atoms2048.zip is 61,997. The total number of atomic configurations spanned by the entire dataset is 190,979. References [1] https://www.lammps.org/ [2] https://openkim.org/ [3] https://openkim.org/id/MEAM_LAMMPS_KimSeolJi_2017_PtNi__MO_020840179467_001 [4] El-Gendy, Ahmed A. and Hampel, Silke and Büccchner, Bernd and Klingeler, Rüdiger, Tuneable magnetic properties of carbon-shielded NiPt-nanoalloys, RSC Adv., volume 6, issue 57, pages 52427-52433, 2016, The Royal Society of Chemistry, doi:10.1039/C6RA05910D

36 MATERIALS SCIENCE↗
DOE OSTI2025

Predicting the viscoplastic response of a crystallizing fluoropolymer using transient network theory

We employ a molecular theory of dynamic polymer networks to describe the viscoplastic response of rubbery FK-800, a thermoplastic copolymer of chlorotrifluoroethylene and vinylidene fluoride, over a broad range of thermal histories. The kinetics of crystallization at different annealing temperatures was modeled using a modified Avrami equation, whose parameters were found to evolve through simple relationships over the full temperature range of the rubbery state. By fitting experimental compression data, we discovered predictable trends for the physical parameters in our mechanical model over its full range of crystallinities (up to ≈20%) and provided insights based on molecular-level physics to justify them. Using this, an end-to-end model was developed to predict the yielding and post-yield behavior of rubbery FK-800 for arbitrary thermal histories. The model successfully predicted the highly nonlinear evolution of characteristic mechanical signatures (stiffness, yield point, post-yield drop) throughout the crystallization process. A statistical analysis of variance test was employed to determine that the measured variations in the mechanical behavior of rubbery FK-800 are primarily dictated by its fractional crystallinity, regardless of its exact thermal history.

36 MATERIALS SCIENCE↗
DOE OSTI2024

Onset of kinetic effects on Rayleigh–Taylor instability: Advective–diffusive asymmetry

In nature and engineering applications, the Rayleigh–Taylor instability (RTI) occurs over a wide range of Atwood, Reynolds, Mach, and Knudsen numbers. At low Atwood, Mach, and Knudsen numbers, the classic advective instability causes quasi-symmetric bubble and spike growth on the two sides of the interface. However, recent findings suggest that at high degrees of rarefaction, advective effects are suppressed and molecular diffusion leads to planar growth of the density fronts on either side of the interface. This study aims to investigate the flow physics of the transition from advective to diffusive behavior, focusing on the onset of the kinetic effects. Using the gas kinetic methodology, RTI is simulated over a range of Knudsen and Mach numbers in the transition regime. The simulation results reveal the various stages of transformation from advective instability to diffusive transport. For the first time, the study demonstrates the existence of a Knudsen–Mach parameter regime where the bubble side exhibits advective instability, while the other side shows a planar density front due to molecular diffusion, rather than the canonical advective spike shape. The dominance of different mechanisms on the two sides of the interface leads to the advective–diffusive asymmetry. In conclusion, the findings of this study can lead to a more comprehensive understanding of RTI over a wide range of Mach and Knudsen numbers.

42 ENGINEERING↗
DOE OSTI2024

Diagonalizing the Born–Oppenheimer Hamiltonian via Moyal perturbation theory, nonadiabatic corrections, and translational degrees of freedom

This article describes a method for calculating higher order or nonadiabatic corrections in Born–Oppenheimer theory and its interaction with the translational degrees of freedom. The method uses the Wigner–Weyl correspondence to map nuclear operators into functions on the classical phase space and the Moyal star product to represent operator multiplication on those functions. These are explained in the body of the paper. The result is a power series in κ 2 , where κ = (m/M)1/4 is the usual Born–Oppenheimer parameter. The lowest order term is the usual Born–Oppenheimer approximation, while higher order terms are nonadiabatic corrections. These are needed in calculations of electronic currents, momenta, and densities. The separation of nuclear and electronic degrees of freedom takes place in the context of the exact symmetries (for an isolated molecule) of translations and rotations, and these, especially translations, are explicitly incorporated into our discussion. This article presents an independent derivation of the Moyal expansion in molecular Born–Oppenheimer theory. Here, we show how electronic currents and momenta can be calculated within the framework of Moyal perturbation theory; we derive the transformation laws of the electronic Hamiltonian, the electronic eigenstates, and the derivative couplings under translations; we discuss in detail the rectilinear motion of the molecular center of mass in the Born–Oppenheimer representation; and we show how the elimination of the translational components of the derivative couplings leads to a unitary transformation that has the effect of exactly separating the translational degrees of freedom.

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

Recovery of release cloud from laser shock-loaded graphite and hydrocarbon targets: in search of diamonds

Abstract This work presents first insights into the dynamics of free-surface release clouds from dynamically compressed polystyrene and pyrolytic graphite at pressures up to 200 GPa, where they transform into diamond or lonsdaleite, respectively. These ejecta clouds are released into either vacuum or various types of catcher systems, and are monitored with high-speed recordings (frame rates up to 10 MHz). Molecular dynamics simulations are used to give insights to the rate of diamond preservation throughout the free expansion and the catcher impact process, highlighting the challenges of diamond retrieval. Raman spectroscopy data show graphitic signatures on a catcher plate confirming that the shock-compressed PS is transformed. First electron microscopy analyses of solid catcher plates yield an outstanding number of different spherical-like objects in the size range between ten(s) up to hundreds of nanometres, which are one type of two potential diamond candidates identified. The origin of some objects can unambiguously be assigned, while the history of others remains speculative.

70 PLASMA PHYSICS AND FUSION TECHNOLOGY↗
DOE OSTI2025

Generalizable machine learning potentials for quantum-accurate predictions of non-equilibrium behavior in 2D materials

Machine learning interatomic potentials (ML-IAPs) are emerging as transformative tools in materials modeling, promising quantum-level accuracy at a fraction of the computational cost. However, their ability to generalize beyond equilibrium configurations and to reliably capture defect- and temperature-driven behavior remains underexplored. Here, we develop and benchmark two state-of-the-art ML-IAPs, Spectral Neighbor Analysis Potential (SNAP) and Allegro, on a comprehensive dataset for monolayer MoSe₂. Using density functional theory (DFT) as the reference, we evaluate their performance in capturing stress–strain behavior, phase transition energetics, defect evolution, edge stability, and fracture toughness. Allegro, a deep equivariant neural network potential, surpasses both SNAP and the classical Tersoff potential in accuracy, efficiency, and transferability. Importantly, both ML potentials accurately reproduce experimental fracture measurements and ab initio predictions of inversion domain formation—phenomena well beyond their training sets. Our findings establish ML-IAPs as viable replacements for traditional force fields in the study of non-equilibrium mechanical phenomena, enabling large-scale, high-fidelity simulations in 2D materials and beyond. In conclusion, this work provides a broadly applicable framework for data-driven modeling of structural and functional transformations under extreme conditions.

2D materials↗
DOE OSTI2020

NESSi : The N on- E quilibrium S ystems S imulation package

The nonequilibrium dynamics of correlated many-particle systems is of interest in connection with pump–probe experiments on molecular systems and solids, as well as theoretical investigations of transport properties and relaxation processes. Nonequilibrium Green’s functions are a powerful tool to study interaction effects in quantum many-particle systems out of equilibrium, and to extract physically relevant information for the interpretation of experiments. Here, we present the open-source software package NESSi (The Non-Equilibrium Systems Simulation package) which allows to perform many-body dynamics simulations based on Green’s functions on the L-shaped Kadanoff–Baym contour. NESSi contains the library libcntr which implements tools for basic operations on these nonequilibrium Green’s functions, for constructing Feynman diagrams, and for the solution of integral and integro-differential equations involving contour Green’s functions. The library employs a discretization of the Kadanoff–Baym contour into time points and a high-order implementation of integration routines. The total integrated error scales up to $\mathcal{O}(N^{-7})$, which is important since the numerical effort increases at least cubically with the simulation time. A distributed-memory parallelization over reciprocal space allows large-scale simulations of lattice systems. We provide a collection of example programs ranging from dynamics in simple two-level systems to problems relevant in contemporary condensed matter physics, including Hubbard clusters and Hubbard or Holstein lattice models. The libcntr library is the basis of a follow-up software package for nonequilibrium dynamical mean-field theory calculations based on strong-coupling perturbative impurity solvers.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2021

Multimode two-dimensional vibronic spectroscopy. II. Simulating and extracting vibronic coupling parameters from polarization-selective spectra

We report experimental demonstrations of polarization-selection two-dimensional Vibrational-Electronic (2D VE) and 2D Electronic-Vibrational (2D EV) spectroscopies aim to map the magnitudes and spatial orientations of coupled electronic and vibrational coordinates in complex systems. The realization of that goal depends on our ability to connect spectroscopic observables with molecular structural parameters. In this paper, we use a model Hamiltonian consisting of two anharmonically coupled vibrational modes in electronic ground and excited states with linear and bilinear vibronic coupling terms to simulate polarization-selective 2D EV and 2D VE spectra. We discuss the relationships between the linear vibronic coupling and two-dimensional Huang–Rhys parameters and between the bilinear vibronic coupling term and Duschinsky mixing. We develop a description of the vibronic transition dipoles and explore how the Hamiltonian parameters and non-Condon effects impact their amplitudes and orientations. Using simulated polarization-selective 2D EV and 2D VE spectra, we show how 2D peak positions, amplitudes, and anisotropy can be used to measure parameters of the vibronic Hamiltonian and non-Condon effects. This paper, along with the first in the series, provides the reader with a detailed description of reading, simulating, and analyzing multimode, polarization-selective 2D EV and 2D VE spectra with an emphasis on extracting vibronic coupling parameters from complex spectra.

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

Deep glassy state dynamic data challenge glass models: Configurational entropy models

This work was the result of the student's (Dongjie Chen) qualifying exam response. In it the experimental observation from our labs that the relaxation times of ultrastable glasses (a 20 million year old amber and a vapor deposited amorphous teflon) do not follow the expected Vogel-Fulcher or WLF types of divergence as the temperatures fall below the glass transition temperature. In the work we evaluated several theories by considering them in different ways. We used them as generally used and show that they do not capture the non-diverging behavior. Because the models were all entropy-based models, we also considered the possibility that the temperature dependence of the entropy used in the models is incorrect and chose two ways of estimating the entropy that do not follow the classical theoretical models. First we used the actually measured calorimetric response for the two ultrastable glasses, and we used the Milchev model, which does not go to zero above 0 Kelvin. Both of these approaches led to improved agreement between the theories and the experimental data. Finally, we also asked what the temperature dependence of the entropy would have to be so that the models agreed with the observed data. In some cases the entropy so-calculated was reasonable, while in others it was unphysical. In sum, the important point of the work is that the observation of non-diverging time-scales that suggest that there is no ideal glass transition can be reasonably captured by current theories if they take more realistic models of the entropy than is usually the case.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Site Density Functional Theory and Structural Bioinformatics Analysis of the SARS-CoV Spike Protein and hACE2 Complex

The entry of the SARS-CoV-2, a causative agent of COVID-19, into human host cells is mediated by the SARS-CoV-2 spike (S) glycoprotein, which critically depends on the formation of complexes involving the spike protein receptor-binding domain (RBD) and the human cellular membrane receptor angiotensin-converting enzyme 2 (hACE2). Using classical site density functional theory (SDFT) and structural bioinformatics methods, we investigate binding and conformational properties of these complexes and study the overlooked role of water-mediated interactions. Analysis of the three-dimensional reference interaction site model (3DRISM) of SDFT indicates that water mediated interactions in the form of additional water bridges strongly increases the binding between SARS-CoV-2 spike protein and hACE2 compared to SARS-CoV-1-hACE2 complex. By analyzing structures of SARS-CoV-2 and SARS-CoV-1, we find that the homotrimer SARS-CoV-2 S receptor-binding domain (RBD) has expanded in size, indicating large conformational change relative to SARS-CoV-1 S protein. Protomer with the up-conformational form of RBD, which binds with hACE2, exhibits stronger intermolecular interactions at the RBD-ACE2 interface, with differential distributions and the inclusion of specific H-bonds in the CoV-2 complex. Further interface analysis has shown that interfacial water promotes and stabilizes the formation of CoV-2/hACE2 complex. This interaction causes a significant structural rigidification of the spike protein, favoring proteolytic processing of the S protein for the fusion of the viral and cellular membrane. Moreover, conformational dynamics simulations of RBD motions in SARS-CoV-2 and SARS-CoV-1 point to the role in modification of the RBD dynamics and their impact on infectivity.

3DRISM↗
DOE OSTI2021

High-pressure phase diagram of beryllium from ab initio free-energy calculations

In this report we use first-principles molecular dynamics simulations coupled with the thermodynamic integration method to study the hexagonal close-packed (hcp) to body-centered cubic (bcc) transition and melting of beryllium up to a pressure of 1600 GPa. We derive the melting line by equating solid and liquid Gibbs free energies and represent it by a Simon-Glatzel fit T m = 1564 K [1 + P/(15.6032 GPa)] 0.383 , which is in good agreement with previous two-phase simulations <6000 K. We also derive the hcp-bcc solid-solid phase boundary and show that the quasiharmonic approximation underestimates the stability of the hcp structure, predicting lower transition pressures between hcp and bcc phases. Our results are consistent with the stability regime predicted by the phonon quasiparticle method. We also predict that the hcp-bcc-liquid triple point is located at 164.7 GPa and 4314 K. In addition, we compute the shock Hugoniot curve and show that it is in good agreement with experiments, intersecting our derived melting curve at ~235 GPa and 4900 K. Finally, we make predictions for future ramp compression experiments. Starting with an isentropic compression of the liquid, we predict the path to intersect the melting line at low pressure and temperature, then to continue along the melting line over a large temperature interval of 7000 K as the sample remains in the mixed solid-liquid state before it enters the solid phase.

36 MATERIALS SCIENCE↗
DOE OSTI2020

ChIMES: A Machine-Learned Interatomic Model Targeting Improved Description of Condensed Phase Chemistry in Energetic Materials

In this report we detail completion of a Physics and Engineering Model Level Two Milestone targeting improved reactive interatomic potentials (IAPs) for energetic materials (EM) through machine learning. The specific goals of this milestone were to develop, validate, and document a new reactive molecular dynamics method for EM, based on machine learning by (1) generating databases of first-principles-derived forces, stresses, and energies for HN3 and 3,4-bis(3-nitrofurazan- 4-yl)furoxan (DNTF) (2) generate atomistic force fields from these databases via ML, and (3) benchmark model performance against first principles calculations. These goals were achieved by (1) further developing a machine learned reactive IAP and generation approach (i.e. the Chebyshev Interaction Model for Efficient Simulation or “ChIMES”), for which resulting IAPs can approach the predictive power of quantum-mechanical approaches at a fraction of the computational expense, and (2) applying the ChIMES framework to develop models for HN3 and DNTF. We find that for simple energetic materials like HN3, high accuracy ChIMES models can be obtained through application of a fitting approach that does not use active machine learning. We demonstrate the suitability of ChIMES models for simulations involving EM by using the HN3 model in multiscale shock technique simulations to predict the HN3 Chapman-Jouguet detonation state and investigate chemical evolution out to 1 ns following shock compression. This model is then used in larger direct shock (DS) simulations for a preliminary investigation of how bubbles (i.e. voids) influence material response under shock compression. We find that more complex EM (i.e. DNTF) necessitate a more sophisticated fitting approach, and develop a new active learning method and python tool to meet this challenge. We demonstrate that this fitting approach yields ChIMES models that out-perform commonly used standard reactive IAPs as well as semi-empirical quantum methods, and discuss the systematic improvability of these actively learned ChIMES models. We also describe challenges related to model development for EM such as DNTF, for which few experimental or previous simulation data are available (e.g. which could otherwise inform generation of training data). To overcome this issue, we establish a semi-empirical quantum ChIMES capability which can be used to efficiently map out relevant thermodynamic and configurational space, and generate ChIMES-IAP training data in a multiscale manner. We also show that these semi-empirical quantum ChIMES models can be used to generate predictions for the shock Hugoniot (the Hugoniot is the locus of thermodynamic states found in a shocked material) equation of state, investigate related thermochemistry, and explore carbon condensation following shock compression. This work represents a substantial advance in our atomistic modeling capability for EM that will provide much needed information on the chemistry of detonation for continued development of continuum models based on the Cheetah thermochemical code.

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

Time-resolved X-ray diffraction diagnostic development for the National Ignition Facility

Here we present the development of an experimental platform that can collect four frames of x-ray diffraction data along a single line of sight during laser-driven, dynamic-compression experiments at the National Ignition Facility. The platform is comprised of a diagnostic imager built around ultrafast sensors with a 2-ns integration time, a custom target assembly that serves also to shield the imager, and a 10-ns duration, quasi-monochromatic x-ray source produced by laser-generated plasma. We demonstrate the performance with diffraction data for Pb ramp compressed to 150 GPa and illuminated by a Ge x-ray source that produces ∼7 × 10 11 , 10.25-keV photons/ns at the 400 μm diameter sample.

36 MATERIALS SCIENCE↗