Engineering Papers⌕ Search

SEARCH · Engineering Papers

Results for “molecular explosives”

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

Sensitivity of pore collapse heating to the melting temperature and shear viscosity of HMX

A multiscale modeling strategy is used to quantify factors governing the temperature rise in hot spots formed by pore collapse from supported and unsupported shock waves in the high explosive HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). Two physical aspects are examined in detail, namely the melting temperature and liquid shear viscosity. All-atom molecular dynamics simulations of phase coexistence are used to predict the pressure-dependent melting temperature up to 5 GPa. Equilibrium simulations and the Green–Kubo formalism are used to obtain the temperature- and pressure-dependent liquid shear viscosity. Starting from a simplified continuum-based grain-scale model of HMX, in this study we systematically increase the complexity of treatments for the solid–liquid phase transition and liquid shear viscosity in simulations of pore collapse. Using a realistic pressure-dependent melting temperature completely suppresses melting for supported shocks, which is otherwise predicted when treating it as a constant determined at atmospheric pressure. Alternatively, melt pools form around collapsed pores when the pressure (and melting temperature) are reduced during the release stage of unsupported shocks. Capturing the pressure dependence of the shear viscosity increases the peak temperature of melt pools by hundreds of Kelvin through viscous work. The complicated interplay of the solid-phase plastic work, solid–liquid phase transition, and liquid-phase viscous work identified here motivate taking a systematic approach to building increasingly complex grain-scale models.

36 MATERIALS SCIENCE↗
DOE OSTI2026

Anisotropic Hot Spot Formation at a Grain Boundary in Shock-Compressed TATB High Explosive Crystal

Secondary high explosives (HEs) exhibit rich microstructure that promotes the formation of hot spots responsible for detonation initiation, but the role of microstructural interfaces remains poorly quantified. To this end, we develop extensions for the generalized crystal-cutting method (GCCM) to prepare molecular dynamics (MD) simulation cells containing grain boundaries (GBs) and other crystal–crystal interfaces with prescribed tilt and twist orientations. Using the GCCM, we perform MD simulations of shock interactions with a GB between the (001) and (100) crystal facets in the secondary HE TATB (1,3,5-triamino-2,4,6-trinitrobenzene). Our MD simulations reveal a strong directional dependence to the formation of a hot spot at the GB interface. In particular, transmission of the shock from the (001) grain to the (100) grain yields a hot spot in the (100) grain at the GB interface, whereas no hot spot is produced when an equivalent shock transits the GB in the opposite direction. We trace the origin of this GB anisotropy to three dominant factors: (1) the intrinsic differences in shock-deformation mechanisms and wave structures for the bulk (100) and (001) grains, which leads to distinct geometries and mechanical impedances upon shock arrival to the GB depending on which grains are donor or acceptor for the transmitted shock; (2) the different time intervals separating the initial shock rise and the formation of steady wave structures in the respective donor–acceptor configurations; and (3) the differences in time scales required to re-establish local thermal equilibrium. Interfacial hot spots form when these factors combine to impede development of a steady two-wave structure and instead induce a localized, pseudosingly shocked region that undergoes a higher rate of work production (resulting in a higher temperature) compared to when the steady two-wave structure develops further from the interface. The extensions to the GCCM approach presented here are anticipated to facilitate a wide range of MD studies that focus on understanding the role of crystal–crystal interfaces in molecular materials.

organic↗
DOE OSTI2020

Comprehensive Data Analysis and Analytic Method Development for PBX9501 Material (FY2019 Annual Report of Aging and Lifetimes Program)

In the polymer bonded explosive (PBX) 9501, a binder consisting of 2.5 wt% nitroplasticizer (NP) and 2.5 wt% Estane R 5703 (Estane) is combined with 94.9 wt% of the high-explosive HMX and 0.1 wt% stabilizer. Because of its flexibility and tensile strength, this binder lowers the sensitivity and improves the manufacturability of PBX 9501. However, like many plasticizers comprised of low molecular weight components, NP has a tendency to diffuse out of the PBX 9501 matrix and can decompose at moderate temperatures into reactive byproducts, such as NO, NO 2 , H 2 O, and HNO x . Through oxidation and hydrolysis, these molecules can further degrade NP and Estane, ultimately degrading the properties of PBX 9501. While gaseous molecules can readily diffuse out of the PBX 9501 charges, the intermediates with lower volatility are more likely trapped in the condensed phase inside the charges and, in a closed system, will co-exist with Estane for an extended period of time creating an ongoing reactive environment. Understanding the rates of formation for these volatiles and intermediates during the NP degradation is therefore a critical prerequisite for understanding the long-term stability of polymeric binders used in munition systems. To work toward this ultimate goal, in the past year, we have conducted comprehensive studies on the physical properties of NP and finished three-year long aging experiment to understand the aging behavior of NP under thermal treatment. Important results are summarized in this annual report. Although NP is widely used in the DOE complex, their physical properties are rather scattered and inconsistent in the open literature. For example, there are at least two widely different values, 14.5°C and -15°C, for the melting point of NP. Although it is known that NP is a eutectic 50:50 mixture of BDNPA and BDNPF, their eutectic phase diagram is not well documented. Furthermore, the effect of temperature on the miscibility between water and eutectic BDNPA/F mixture is rarely reported. To fill these knowledge gaps, a large set of BDNPA/F mixtures with BDNPA concentration ranging from 0 to 100 wt% was analyzed using DSC techniques. In addition to determining the eutectic melt point as -25°C, a phase diagram of the BDNPA/F system was constructed from -30°C to 45°C. With this phase diagram, the phase transition temperatures and composition can be readily found for the BDNPA/F mixtures with various mass ratios.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Optical Microsystems for Molecular Sensing (CRADA Final Report)

As part of the Cyclotron Road program, NuSight Photonics sought to implement and investigate an optical Raman spectroscopy instrument for material and molecular analysis. Raman spectroscopy is widely used for quick and non-invasive material measurement when accurate results are needed – e.g., for monitoring manufacturing process in pharmaceutical plants, or identifying illicit substances and explosives by Law Enforcement. This technique is also growing rapidly in popularity as it is applied successfully to a broad range of applications, from pathology and non-invasive early cancer detection, to food and drug safety. The biggest challenge for Raman spectroscopy to achieve its full potential is the hardware cost, ranging from $\$ $25k to $\$ $250k, with the lower cost solutions significantly limited in sensitivity—or detection limit. NuSight Photonics aimed to tackle this obstacle through a novel hardware approach: Swept-Source Raman Spectroscopy.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Trace explosive residue detection of HMX and RDX in post-detonation dust from an open-air environment

Explosives are often used in industry, geology, mining, and other applications, but it is not always clear what remains after a detonation or the fate and transport of any residual material. The goal of this study was to determine to what extent intact molecules of high explosive (HE) compounds are detectable and quantifiable from post-detonation dust and particulates in a field experiment with varied topography. We focused on HMX (1,3,5,7-Tetranitro-1,3,5,7-tetrazocane), which is less studied in field detonation literature, as the primary explosive material and RDX (1,3,5-Trinitroperhydro-1,3,5-triazine) as the secondary material. The experiment was conducted at Site 300, Lawrence Livermore National Laboratory’s Experimental Test Site, in California, USA. Two 20.4 kg and one 40.8 kg above ground explosions (primarily comprised of LX-14, an HMX-based polymer-bonded high explosive) were detonated on an open-air firing area on separate days. The complex terrain of the firing area (e.g., buildings, berm, low-height obstacles) was advantageous to study HE deposition in relation to plume dynamics. Three types of samples were collected up to 100 m away from each shot: surface swipes of aluminum plates, surface swipes of fixed objects, and filters from air samples. We used atmospheric flow tube-mass spectrometry (AFT-MS) to quantify picogram levels of molecular residue of HE material in the post-detonation dust. An aliquot of sample extract in methanol (e.g., 1 µL of 0.5 mL) was placed onto a resistive material and then thermally desorbed into the AFT-MS. We successfully detected and quantified both HMX and RDX in many of the samples. Based on mass (pg) detected and solution dilution, we back-calculated the mass collected on the swipe or filter (ng per sample). The aerial distribution of molecular residue was consistent with the path of the plume, which was strongly determined by wind speed and direction at the time of each shot. The quantity of material detected appeared to correlate more with distance from the shot and the wind conditions than with shot size. This study demonstrates that the picogram detection levels of AFT-MS are well-suited for quantification of analytes (e.g., HMX and RDX) in environmental samples.

atmospheric flow tube-mass spectrometry (AFT-MS), ↗
DOE OSTI2025

Visualizing the strong field–induced molecular breakup of C 60 via x-ray diffraction

Laser-driven dynamics in polyatomic molecules poses a complex many-body problem. Understanding intense light-matter interaction is crucial for steering intramolecular quantum dynamical processes. Here, we record time-resolved x-ray diffraction images of C 60 molecules during and after their interaction with intense near-infrared fields, giving direct access to structural changes of the molecules and their fragmentation in real time. Tuning the intensity of the excitation pulses, we uncover a transition from a weak-field regime of excited but stable molecules to a high-field regime dominated by Coulomb explosion. In the transition region, the molecules expand by up to 50% of their initial size within just 140 fs, with major fragmentation only setting in afterward. This work demonstrates that x-ray diffractive imaging is capable of retrieving time-resolved structural information of large molecules reshaped by intense laser fields. Laser-driven fragmentation is a first step toward observing molecular processes modified by laser fields of increasing intensity.

Schnorr, Kirsten [Paul Scherrer Inst. (PSI), Villi↗
DOE OSTI2021

All-Atom Simulation of 3D Hot Spot Formation in Shocked TATB Explosive

TATB is an insensitive high explosive (IHE) critical to the stockpile that is challenging to model at the continuum scale. Advanced detonation models in the Cheetah high explosive chemistry code require validation though subscale simulations. High explosive initiation is determined by micron-scale physics of hot spots formed a shock-collapsed pores. Pore sizes between 100 nm and 1 μm are believed to be the most important for determining the shock sensitivity of TATB. This range of pore sizes is difficult to access at the atomic scale through allatom molecular dynamics (MD) simulations, even with Sierra-class computers. Quasi-2D simulations are widely used and allow much larger pore sizes (up to 400 nm) to be studied, but the applicability of 2D simulations to the actual 3D pore response is not understood. Resolving these uncertainties through “full physics” MD modeling is key for generalizing, parameterizing, and validating the kinds of continuum models used to inform design, safety, and performance. This work was a continuation of FY20 efforts pushing simulations to full 3D with the largest-ever all-atom simulations of an explosive. These were the first all-atom full-3D simulations of large hot spots thought to govern explosive detonation and required over a billion atoms. Simulations were performed using LAMMPS, an open SNL science code. MD explosive models present unique challenges, even for established codes such as LAMMPS. Their model forms are more complex than typical models for metals, while simulating high temperature-pressure conditions is demanding and increases computational cost. Scaling problems in GPU-enabled MD algorithms initially limited simulations to <100 million atoms but were resolved through collaboration with SNL. An overall 24x speedup was obtained relative to CPU machines. Specialized analysis of these simulations required a bottom-up refactoring and algorithm parallelization of in-house codes and application of computer vision algorithms to extract meaningful information.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Modeling Framework to Predict Melting Dynamics at Microstructural Defects in TNT-HMX High Explosive Composites

Many high explosive (HE) formulations are composite materials whose microstructure is understood to impact functional characteristics. Interfaces are known to mediate the formation of hot spots that control their safety and initiation. Here, to study such processes at molecular scales, we developed all-atom force fields (FFs) for Octol, a prototypical HE formulation comprised of TNT (2,4,6-trinitrotoluene) and HMX (octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine). We extended a FF for TNT and recasted it in a form that can be readily combined with a well-established FF for HMX. The resulting FF was extensively validated against experimental results and density functional theory calculations. We applied the new combined TNT-HMX FF to predict and rank surface and interface energies, which indicate that there is an energetic driver for coarsening of microstructural grains in TNT-HMX composites. Finally, we assess the impact of several microstructural environments on the dynamic melting of TNT crystal under ultrafast thermal loading. We find that both free surfaces and planar material interfaces are effective nucleation points for TNT melting. However, MD simulations show that TNT crystal is prone to superheating by at least 50 K on subnanosecond time scales and that the degree of superheating is inversely correlated with surface and interface energy. The modeling framework presented here will enable future studies on hot spot formation processes in accident scenarios that are governed by strong coupling between microstructural interfaces, material mechanics, momentum and energy transport, phase transitions, and chemistry.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Phosgene Neutralization Chemistry for the Explosive Destruction System (EDS)

This report summarizes activities at Sandia National Laboratories as part of the Explosive Destruction System (EDS) Phase 3 (P3) System design. An exploration of chemical neutralization strategies of phosgene was conducted for safe disposal of recovered mortars and M79 1000 lb. bombs filled with carbonyl dichloride "phosgene" or "CG agent" (molecular formula = COCl 2 ). The incumbent strategy utilized aqueous sodium hydroxide was found to be the worst-case scenario, producing enough CO 2 gas that would cause an unacceptable pressure and temperature spike. Several chemical neutralization strategies were evaluated based on criteria set by the operating envelope of the P3 design. In the end, it was determined that a pure solution of N-methyl ethanolamine (MeEA) or 90% monoethanolamine (MEA)(aqueous) provided the best balance reaction profile, cost, and safety.

36 MATERIALS SCIENCE↗
DOE OSTI2025

Scale Invariance of Hot Spot Formation in TATB High Explosives

Shock-induced detonation of insensitive high explosives based on 1,3,5-triamino-2,4,6-trinitrobenzene starts with formation of hot spots at microstructural defects but has eluded atomistic modeling treatment at micron length scales. To this end, we performed multimicron scale all-atom molecular dynamics (MD) simulations of hot spots that form during the collapse of cylindrical pores with diameters between 10 and 300 nm. Our MD simulations show that hot spots formed at pores larger than 20 nm exhibit temperature fields with scale-invariant features for sizes up to at least 300 nm. Through a continuum-based grain-scale modeling framework, we span and extend beyond the size scales currently accessible to MD and find that hot spot scale invariance is a general feature that arises when the mechanical strength is insensitive to strain rate. Finally, our results demonstrate the applicability of all-atom MD to simulate the complicated dynamical evolution of micron-sized systems and bolster confidence in insights from MD simulations of materials that exhibit strength with negligible rate dependence over the relevant intervals.

36 MATERIALS SCIENCE↗
DOE OSTI2021

Predicted Reaction Mechanisms, Product Speciation, Kinetics, and Detonation Properties of the Insensitive Explosive 2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105)

2,6-Diamino-3,5-dinitropyrazine-1-oxide (LLM-105) is a relatively new and promising insensitive high-explosive (IHE) material that remains only partially characterized. IHEs are of interest for a range of applications and from a fundamental science standpoint, as the root causes behind insensitivity are poorly understood. In this work, we adopt a multitheory approach based on reactive molecular dynamic simulations performed with density functional theory, density functional tight-binding, and reactive force fields to characterize the reaction pathways, product speciation, reaction kinetics, and detonation performance of LLM-105. We compare and contrast these predictions to 1,3,5-triamino-2,4,6-trinitrobenzene (TATB), a prototypical IHE, and 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX), a more sensitive and higher performance material. The combination of different predictive models allows access to processes operative on progressively longer timescales while providing benchmarks for assessing uncertainties in the predictions. We find that the early reaction pathways of LLM-105 decomposition are extremely similar to TATB; they involve intra- and intermolecular hydrogen transfer. Additionally, the detonation performance of LLM-105 falls between that of TATB and HMX. We find agreement between predictive models for first-step reaction pathways but significant differences in final product formations. Predictions of detonation performance result in a wide range of values, and one-step kinetic parameters show the similar reaction rates at high temperatures for three out of four models considered.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2026

An Updated Picture of Presolar History from Short-lived Radioactive Isotopes and Inferences on the Birth of the Sun

We examine the origin of the short-lived radionuclides (SLRs; defined as having half-lives between 0.1 and 100 Ma) present in the early solar system (ESS) by investigating how predictions of their abundances in the interstellar medium (ISM) from steady-state equilibrium relate to their ESS values. For this, we take into account the non-negligible time t iso elapsed between the isolation of the presolar molecular cloud and the formation of the ESS, during which the SLRs decayed freely. We also consider the alternative scenario in which the presolar molecular cloud remained partially mixed with the ISM, with a mixing timescale t mix . We find that the ESS abundances of 107 Pd and 182 Hf produced by slow neutron captures (s-process) and of 53 Mn and 60 Fe produced by explosive nucleosynthesis can be consistently explained within these scenarios. Their required tiso is 9–12 Ma and their required tmix is 11–14 Ma (with one potential exception of t mix = 38 Ma), depending on galactic uncertainties, such as the galactic star formation history and efficiency and the star-to-gas mass ratio. Another s -process SLR, 205 Pb, has a more uncertain ESS value and falls within only some of these time values. The same applies to the SLRs produced by the p -process ( 92 Nb and 146 Sm), depending on the latter’s half-life. In agreement with previous studies, we find that the ESS abundances of the rapid neutron-capture isotopes ( 129 I, 244 Pu, and 247 Cm) and of the most short-lived radionuclides ( 26 Al, 36 Cl, and 41 Ca) cannot be explained by assuming steady-state equilibrium in the ISM.

79 ASTRONOMY AND ASTROPHYSICS↗
DOE OSTI2023

Predicted viscosity of liquid HMX up to 40 GPa

Viscous flow serves as a significant heating mechanism during the formation of hot spots, but the shear viscosity which determines this response is poorly characterized for most high explosives. Most recently, a model was proposed for the shear viscosity of liquid HMX (1,3,5,7-tetranitro-1,3,5,7-tetrazocane) that was fit to pressures reaching 5 GPa, but this work uncovered uncertainties in the viscosity at 0 GPa and remains untested at higher pressures. We use molecular dynamics (MD) simulations and the Green-Kubo formalism to predict the temperature- and pressure-dependent shear viscosity of HMX over the pressure interval 0 GPa≤P≤40 GPa. Reassessment of the viscosity at 0 GPa rules out several potential explanations for discrepancies between earlier reports; we attribute these differences to details of MD trajectory integration impacting molecular flexibility. The shear viscosity of HMX exhibits an Arrhenius temperature dependence at each pressure considered, with exponential prefactor and activation energy terms that are also strong functions of pressure. An analytic form for the viscosity is developed based on an extension of the well-known Andrade equation that simultaneously captures the temperature and pressure dependencies in the MD data up to 40 GPa. Comparison against a recently developed model for the viscosity of liquid RDX (1,3,5-trinitro-1,3,5-triazinane) shows that both materials exhibit similar functional dependencies with the viscosity of HMX being higher by roughly an order of magnitude at a given temperature-pressure state.

36 MATERIALS SCIENCE↗
DOE OSTI2020

Computational study of TATB

We have performed a computational study of TATB. The study is composed of two parts, one where we perform static T = 0 K calculations with several different DFT functionals, to investigate structural properties and obtain a cold curve. Even though the functionals used in this study give poor results at ambient pressure, they perform much better at compression, when the dominance of the van der Waals’ forces in the binding is replaced by more normal interactions. In particular AM05 gives good results, mainly because it doesn’t include any “faulty” van der Waals’. Indeed the van der Waals’ corrected functional (PBE+D2) results shows the worst performance on structural properties. For the cold curve the PBE+D2 functional gives the best results compared to experiments (room temperature) but since the structural properties were bad with this functional we cannot for sure say that this is in any way a validation. It is shown that scaling the AM05 results with a common factor reproduces the experimental data. We leave the question on why to a follow up study but speculate that maybe the experiment and our calculations do not represent the same system set-up. From the structural properties we identify a site in the TATB lattice that contains a high concentration of oxygen and hydrogen. The consequences of this finding, if any, are left to a future project. In the second part of the study we perform extensive high-quality DFT-MD (molecular dynamics with DFT forces) calculations to obtain points on the Hugoniot. This part of the study did not get finished within this project and a more extensive summary of this part will be given in a follow up report at a future date. However, we show that obtaining calculated data is essential for equations of state development since the existing experimental Hugoniot data barely deviates from room temperature isotherm data and thus give no information of the high temperature dependency of materials properties. Our conclusion so far in the project is that being able to calculate materials properties is essential for high explosives since experiments are unable to give information in the high temperature part of phase space. While blind use of DFT cannot (yet) give unambiguous answers, we can still explore qualitatively and use such information to guide us in EOS development. However, the main conclusion is that we need better functionals in DFT to easily extract needed properties from calculations.

36 MATERIALS SCIENCE↗
DOE OSTI2022

Single-Step Mechanism for Regioselective Nitration of 9,10-BN-Naphthalene with Acetyl Nitrate in the Gas Phase

Here, the energetics of the regioselective mononitration of 9,10-BN-naphthalene with acetyl nitrate (H 3 C 2 NO 4 ) were modeled with ab initio simulations in the gas phase and an acetonitrile solvent. The single-electron-transfer (SET) nitration mechanism leading to a σ-complex and a single-step nitration mechanism were modeled. The energy barrier for the single-step mechanism was lower than that for the SET mechanism in the gas phase. However, the two are much more energetically competitive in the solvent. The σ-complex was found to be unstable in the gas phase owing to the interaction with the counterion. Using the single-step mechanism, the carbon site 1 nearest boron had the lowest activation energy for nitration of 22.6 kcal/mol, while site 3 had the second lowest barrier of 24.6 kcal/mol. Finally, details on the molecular structures at intermediate and transition states as well as charges in different configurations are discussed.

71 CLASSICAL AND QUANTUM MECHANICS, GENERAL PHYSIC↗
DOE OSTI2025

Shock wave bending around a dusty plasma void

We report on experimental observations of the bending of a dust acoustic shock wave around a dust void region. This phenomenon occurs as a planar shock wavefront encounters a compressible obstacle in the form of a void whose size is larger than the wavelength of the wave. As they collide, the central portion of the wavefront, that is the first to touch the void, is blocked while the rest of the front continues to propagate, resulting in an inward bending of the shock wave. The bent shock wave eventually collapses, leading to the transient trapping of dust particles in the void. Subsequently, a Coulomb explosion of the trapped particles generates a bow shock. The experiments have been carried out in a Direct Current (DC) glow discharge plasma, where the shock wave and the void are simultaneously created as self-excited modes of a three-dimensional dust cloud. The salient features of this phenomenon are reproduced in molecular dynamics simulations, which provide valuable insights into the underlying dynamics of this interaction.

Dust acoustic waves↗
DOE OSTI2025

Layers Can Be Deceiving: A Hopping Model for Small Molecule Diffusion in TATB Crystal

Sorption of small molecule gases in materials can play a significant role in their long‐term stability and compatibility within multi‐material assemblies. While many material properties of the insensitive high explosive TATB (1,3,5‐triamino‐2,4,6‐trinitrobenzene) are well understood, very little is known regarding its permeability to gases. TATB crystal exhibits a graphitic‐like layered packing structure that evokes a mental schema in which the layers form nanoscopic channels, but it is unclear whether this structure promotes gas transport. Here, we use molecular dynamics (MD) simulations to predict transport of small molecules through TATB single crystal. An approach to fit classical force fields is developed to model TATB interactions with H 2 O, He, Ne, and Ar, which is then combined with steered MD to probe gas transport along selected directions in the crystal. We find that small molecule transport occurs via a hopping mechanism that exhibits distinct jumps between interstitial sites and is substantially faster normal to the layers as compared to through them. This result stems from the finding that intralayer junctions between adjacent TATB molecules are the most stable interstitial sites and that energetic barriers are lower for hopping between adjacent layers. An empirical model for diffusion rate based on the MD data shows that the rate decays exponentially with increasing molecular radius and is negligibly small for all molecules larger than He, including common atmospheric gases. These findings have implications for the interpretation of experiments that measure surface area, material response to extreme conditions, and are expected to help constrain models for material aging.

36 MATERIALS SCIENCE↗
DOE OSTI2026

From atomistic models to machine learning: Predictive design of nanocarbons under extreme conditions

The formation of technologically valuable nanocarbon structures under extreme conditions, such as those produced during high-explosive detonations, remains poorly understood but holds significant potential for the development of controlled synthesis pathways. While detonation shockwaves provide the high-pressure, high-temperature environment required for nanodiamond formation, subsequent cooling and decompression dictate whether the diamond phase is preserved or transformed into other nanocarbon structures. Here, in this study, we employ GPU-accelerated reactive molecular dynamics (ReaxFF) simulations to investigate the graphitization and structural remodeling of detonation nanodiamond under nonlinear quench and pressure-release trajectories. We further investigate how the initial nanodiamond morphology; cuboctahedral, octahedral, or hexagonal prism influences the resulting transformation products. Evolution of nanostructure, allotrope (via simulated x-ray diffraction), carbon hybridization, and ring statistics are tracked during a two-stage quench from 5000 K to 60 GPa. Rapid cooling combined with slow decompression optimizes cubic diamond retention, whereas slow cooling with rapid pressure release promotes surface-to-core graphitization, producing concentric sp 2 -hybridized layers and hollowed inner shells. Octahedral nanodiamonds evolve into carbon nano-onions, initially forming bucky diamonds that progressively transform into fully sp 2 -hybridized structures, while hexagonal prisms preferentially form parallel-stacked graphite layers resembling carbon dots. Transient hexagonal diamond (lonsdaleite) emerges as an interfacial phase, suggesting potential reversibility in the shock-induced graphite-to-diamond transformation pathway transformation route. To extend predictive capabilities, we trained machine learning (ML) regressors on over 10 5 node-hours of molecular dynamics (MD) trajectories. A multilayer perceptron (MLP) model reliably predicts the number of graphitized layers from temperature–pressure trajectories with a coefficient of determination (R 2 ) exceeding 0.90. This high predictive fidelity enables efficient, high-throughput mapping of the synthesis parameter space for optimized graphitization outcomes. Collectively, morphological control combined with optimized quench–decompression conditions promote the selective synthesis of nanocarbon allotropes. This work establishes a data-driven framework for the rational, a priori design of carbon nanomaterials for applications in energy storage, sensing, and biomedicine.

Detonation nanodiamond remodeling↗