A methodology to generate crystal-based molecular structures for atomistic simulations.

We propose a systematic method to construct crystal-based molecular structures often needed as input for computational chemistry studies. These structures include crystal 'slabs' with periodic boundary conditions (PBCs) and non-periodic solids such as Wulff structures. We also introduce a method to build crystal slabs with orthogonal PBC vectors. These methods are integrated into our code,Los Alamos Crystal Cut(LCC), which is open source and thus fully available to the community. Examples showing the use of these methods are given throughout the manuscript.

Identifying the Molecular Properties that Drive Explosive Sensitivity in a Series of Nitrate Esters.

Energetic materials undergo hundreds of chemical reactions during exothermic runaway, generally beginning with the breaking of the weakest chemical bond, the "trigger linkage." Herein we report the syntheses of a series of pentaerythritol tetranitrate (PETN) derivatives in which the energetic nitrate ester groups are systematically substituted by hydroxyl groups. Because all the PETN derivatives have the same nitrate ester-based trigger linkages, quantum molecular dynamics (QMD) simulations show very similar Arrhenius kinetics for the first reactions. However, handling sensitivity testing conducted using drop weight impact indicates that sensitivity decreases precipitously as nitrate esters are replaced by hydroxyl groups. These experimental results are supported by QMD simulations that show systematic decreases in the final temperatures of the products and the energy release as the nitrate ester functional groups are removed. To better interpret these results, we derive a simple model based only on the specific enthalpy of explosion and the kinetics of trigger linkage rupture that accounts qualitatively for the decrease in sensitivity as nitrate ester groups are removed.

Reaction Rates in Nitromethane under High Pressure from Density Functional Tight Binding Molecular Dynamics Simulations.

We use density functional tight binding (DFTB) molecular dynamics (MD) simulations to determine the reaction rates of nitromethane CH3NO2 (NM) under high pressure (P = 14-28 GPa), and temperature (T = 1450-1850 K). DFTB-MD simulations performed with the same initial conditions (P0, T0) reveal a stochastic behavior, both in terms of reaction times and chemical paths. By running series of MD simulations, we are able to obtain average reaction times with quantified errors and devise a simple two-step model for NM explosion: ignition/explosion. While our model bypasses the chemical complexity due to the numerous reaction paths and intermediates observed during reactions, the chemistry is accounted for via the accurate parameterization of the DFTB model, and our results suggest a single main reaction pathway for the pressure range considered here, dominated in the earlier stages by the formation of the aci-ion, CH2NOO-. By fitting our data to a Frank-Kamenetskii model, we extract prefactors and pressure-independent activation energies and volumes for the ignition and explosion stages. A two-step model is then built and compared to experimental observations. Single and two-step Arrhenius models are also provided for comparison with literature data. This work presents an efficient way of investigating the reactivity of high explosives by performing electronic structure-based MD simulations and provides reaction rates for simplified models that can be implemented into hydrocodes.

Spatially-varying inversion near grain boundaries in MgAl2O4 spinel.

Complex materials, containing multiple chemical species, often exhibit chemical disorder or inversion. Typically, this disorder is viewed as spatially homogeneous throughout the material. Here, we show, using a simple grain boundary in MgAl2O4 spinel, that this is not the case and that the level of inversion at the grain boundary plane is different than in the bulk. This has ramifications for the energetics of the boundary and how defects interact with it, as exemplified by the relative formation energy of vacancies. Using these results as motivation, we construct a simple model of inversion versus grain size that captures the salient behavior observed in experiments and allows us to extract inversion-relevant properties from those same experiments, suggesting that grain boundaries in the experimentally prepared material are essentially fully inverse. Together, these results highlight the role that microstructure plays on the inversion in the material.

Correction: Spatially-varying inversion near grain boundaries in MgAl2O4 spinel.

[This corrects the article DOI: 10.1039/D0RA00700E.].

Probing ultrafast shock-induced chemistry in liquids using broad-band mid-infrared absorption spectroscopy.

We probe shock-induced chemistry in two organic liquids by measuring broadband, midinfrared absorption in the 800-1400 cm-1 frequency range. To test this new method and understand the signatures of chemical reactions in time resolved vibrational spectra, we compared liquid benzene shocked to unreactive conditions (shocked to a pressure of 18 GPa for a duration of 300 ps) to nitromethane under reactive conditions (25 GPa). We see clear signatures of shock-induced chemistry that are distinguishable from the pressure- and temperature-induced changes in vibrational mode shapes. While shocked benzene shows primarily a broadening and shifting of the vibrational modes, the nitromethane vibrational modes vanish once the shock wave enters the liquid and simultaneously, a spectrally broad feature appears that we interpret as the infrared spectrum of the complex mixture of product and intermediate species. To further interpret these measurements, we compare them to reactive quantum molecular dynamics simulations, which gives qualitatively consistent results. This work demonstrates a promising method for time resolving shock-induced chemistry, illustrating that chemical reactions produce distinct changes in the vibrational spectra.

The impact of chemical order on defect transport in mixed pyrochlores.

Using temperature accelerated dynamics, an accelerated molecular dynamics method, we examine the relationship between composition and cation ordering and defect transport in the mixed pyrochlore Gd2(Ti1-xZrx)2O7, using the oxygen vacancy as a representative defect structure. We find that the nature of transport is very sensitive to the cation structure, transitioning, as a function of composition, from three-dimensional migration to two-dimensional to pseudo-one-dimensional to becoming essentially immobile before reverting back to three-dimensional as the Zr content is increased. The rates of migration are also affected by the cation structure in the various compositions. This behavior is driven by the connectivity of Ti polyhedra in the material, with more extensive networks of Ti ions leading to a greater ability of the vacancy to traverse the material. Our results indicate that the nature of transport is dictated by the cation structure of the material and that, conversely, the cation structure could be used to control transport and potentially other functionalities in mixed pyrochlores.

Linear Scaling Pseudo Fermi-Operator Expansion for Fractional Occupation.

Recursive Fermi-operator expansion methods for the calculation of the idempotent density matrix are valid only at zero electronic temperature with integer occupation numbers. We show how such methods can be modified to include fractional occupation numbers of an approximate or pseudo Fermi-Dirac distribution and how the corresponding entropy term of the free energy is calculated. The proposed methodology is demonstrated and evaluated for different electronic structure methods, including density functional tight-binding theory, Kohn-Sham density functional theory using numerical orbitals, and quantum chemistry Hartree-Fock theory using Gaussian basis functions.

Examining the chemical and structural properties that influence the sensitivity of energetic nitrate esters.

The sensitivity of explosives is controlled by factors that span from intrinsic chemical reactivity and chemical intramolecular effects to mesoscale structure and defects, and has been a topic of extensive study for over 50 years. Due to these complex competing chemical and physical elements, a unifying relationship between molecular framework, crystal structure, and sensitivity has yet to be developed. In order to move towards this goal, ideally experimental studies should be performed on systems with small, systematic structural modifications, with modeling utilized to interpret experimental results. Pentaerythritol tetranitrate (PETN) is a common nitrate ester explosive that has been widely studied due to its use in military and commercial explosives. We have synthesized PETN derivatives with modified sensitivity characteristics by substituting one -CCH2ONO2 moiety with other substituents, including -CH, -CNH2, -CNH3X, -CCH3, and -PO. We relate the handling sensitivity properties of each PETN derivative to its structural properties, and discuss the potential roles of thermodynamic properties such as heat capacity and heat of formation, thermal stability, crystal structure, compressibility, and inter- and intramolecular hydrogen bonding on impact sensitivity. Reactive molecular dynamics (MD) simulations of the C/H/N/O-based PETN-derivatives have been performed under cook-off conditions that mimic those accessed in impact tests. These simulations infer how changes in chemistry affect the subsequent decomposition pathways.

Evidence for percolation diffusion of cations and reordering in disordered pyrochlore from accelerated molecular dynamics.

Diffusion in complex oxides is critical to ionic transport, radiation damage evolution, sintering, and aging. In complex oxides such as pyrochlores, anionic diffusion is dramatically affected by cation disorder. However, little is known about how disorder influences cation transport. Here, we report results from classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd2Ti2O7 pyrochlore, on the microsecond timescale. We find that diffusion is slow at low levels of disorder, while higher disorder allows for fast diffusion, which is then accompanied by antisite annihilation and reordering, and thus a slowing of cation transport. Cation diffusivity is therefore not constant, but decreases as the material reorders. We also show that fast cation diffusion is triggered by the formation of a percolation network of antisites. This is in contrast with observations from other complex oxides and disordered media models, suggesting a fundamentally different relation between disorder and mass transport.Diffusion plays an important role in sintering, damage tolerance and transport. Here authors perform classical and accelerated molecular dynamics simulations of vacancy-mediated cation diffusion in Gd2Ti2O7 pyrochlore and report non-monotonic evolution of cation diffusivity.

Disorder-induced transition from grain boundary to bulk dominated ionic diffusion in pyrochlores.

We use molecular dynamics simulations to investigate the role of grain boundaries (GBs) on ionic diffusion in pyrochlores, as a function of the GB type, chemistry of the compound, and level of cation disorder. We observe that the presence of GBs promotes oxygen transport in ordered and low-disordered systems, as the GBs are found to have a higher concentration of mobile carriers with higher mobilities than in the bulk. Thus, in ordered samples, the ionic diffusion is 2D, localized along the grain boundary. When cation disorder is introduced, bulk carriers begin to contribute to the overall diffusion, while the GB contribution is only slightly enhanced. In highly disordered samples, the diffusive behavior at the GBs is bulk-like, and the two contributions (bulk vs. GB) can no longer be distinguished. There is thus a transition from 2D/GB dominated oxygen diffusivity to 3D/bulk dominated diffusivity versus disorder in pyrochlores. These results provide new insights into the possibility of using internal interfaces to enhance ionic conductivity in nanostructured complex oxides.

Determination of Krypton Diffusion Coefficients in Uranium Dioxide Using Atomic Scale Calculations.

We present a study of the diffusion of krypton in UO2 using atomic scale calculations combined with diffusion models adapted to the system studied. The migration barriers of the elementary mechanisms for interstitial or vacancy assisted migration are calculated in the DFT+U framework using the nudged elastic band method. The attempt frequencies are obtained from the phonon modes of the defect at the initial and saddle points using empirical potential methods. The diffusion coefficients of Kr in UO2 are then calculated by combining this data with diffusion models accounting for the concentration of vacancies and the interaction of vacancies with Kr atoms. We determined the preferred mechanism for Kr migration and the corresponding diffusion coefficient as a function of the oxygen chemical potential µO or nonstoichiometry. For very hypostoichiometric (or U-rich) conditions, the most favorable mechanism is interstitial migration. For hypostoichiometric UO2, migration is assisted by the bound Schottky defect and the charged uranium vacancy, VU4-. Around stoichiometry, migration assisted by the charged uranium-oxygen divacancy (VUO2-) and VU4- is the favored mechanism. Finally, for hyperstoichiometric or O-rich conditions, the migration assisted by two VU4- dominates. Kr migration is enhanced at higher µO, and in this regime, the activation energy will be between 4.09 and 0.73 eV depending on nonstoichiometry. The experimental values available are in the latter interval. Since it is very probable that these values were obtained for at least slightly hyperstoichiometric samples, our activation energies are consistent with the experimental data, even if further experiments with precisely controlled stoichiometry are needed to confirm these results. The mechanisms and trends with nonstoichiometry established for Kr are similar to those found in previous studies of Xe.

The role of surfaces, chemical interfaces, and disorder on plutonium incorporation in pyrochlores.

Pyrochlores, a class of complex oxides with formula A2B2O7, are one of the candidates for nuclear waste encapsulation, due to the natural occurrence of actinide-bearing pyrochlore minerals and laboratory observations of high radiation tolerance. In this work, we use atomistic simulations to determine the role of surfaces, chemical interfaces, and cation disorder on the plutonium immobilization properties of pyrochlores as a function of pyrochlore chemistry. We find that both Pu(3+) and Pu(4+) segregate to the surface for the four low-index pyrochlore surfaces considered, and that the segregation energy varies with the chemistry of the compound. We also find that pyrochlore/pyrochlore bicrystals A2B2O7/A2'B2'O7 can be used to immobilize Pu(3+) and Pu(4+) either in the same or separate phases of the compound, depending on the chemistry of the material. Finally, we find that Pu(4+) segregates to the disordered phase of an order/disorder bicrystal, driven by the occurrence of local oxygen-rich environments. However, Pu(3+) is weakly sensitive to the oxygen environment, and therefore only slightly favors the disordered phase. This behavior suggests that, at some concentration, Pu incorporation can destabilize the pyrochlore structure. Together, these results provide new insight into the ability of pyrochlore compounds to encapsulate Pu and suggest new considerations in the development of waste forms based on pyrochlores. In particular, the phase structure of a multi-phase pyrochlore composite can be used to independently getter decay products based on their valence and size.

The effects of cation-anion clustering on defect migration in MgAl2O4.

Magnesium aluminate spinel (MgAl2O4), like many other ceramic materials, offers a range of technological applications, from nuclear reactor materials to military body armor. For many of these applications, it is critical to understand both the formation and evolution of lattice defects throughout the lifetime of the material. We use the Speculatively Parallel Temperature Accelerated Dynamics (SpecTAD) method to investigate the effects of di-vacancy and di-interstitial formation on the mobility of the component defects. From long-time trajectories of the state-to-state dynamics, we characterize the migration pathways of defect clusters, and calculate their self-diffusion constants across a range of temperatures. We find that the clustering of Al and O vacancies drastically reduces the mobility of both defects, while the clustering of Mg and O vacancies completely immobilizes them. For interstitials, we find that the clustering of Mg and O defects greatly reduces O interstitial mobility, but has only a weak effect on Mg. These findings illuminate important new details regarding defect kinetics relevant to the application of MgAl2O4 in extreme environments.

Insights into dynamic processes of cations in pyrochlores and other complex oxides.

Complex oxides are critical components of many key technologies, from solid oxide fuel cells and superionics to inert matrix fuels and nuclear waste forms. In many cases, understanding mass transport is important for predicting performance and, thus, extensive effort has been devoted to understanding mass transport in these materials. However, most work has focused on the behavior of oxygen while cation transport has received relatively little attention, even though cation diffusion is responsible for many phenomena, including sintering, radiation damage evolution, and deformation processes. Here, we use accelerated molecular dynamics simulations to examine the kinetics of cation defects in one class of complex oxides, A2B2O7 pyrochlore. We find that, in some pyrochlore chemistries, B cation defects are kinetically unstable, transforming to A cation defects and antisites at rates faster than they can diffuse. When this occurs, transport of B cations occurs through defect processes on the A sublattice. Further, these A cation defects, either interstitials or vacancies, can interact with antisite disorder, reordering the material locally, though this process is much more efficient for interstitials than vacancies. Whether this behavior occurs in a given pyrochlore depends on the A and B chemistry. Pyrochlores with a smaller ratio of cation radii exhibit this complex behavior, while those with larger ratios exhibit direct migration of B interstitials. Similar behavior has been reported in other complex oxides such as spinels and perovskites, suggesting that this coupling of transport between the A and B cation sublattices, while not universal, occurs in many complex oxides.