Halogenated PETN derivatives: interplay between physical and chemical factors in explosive sensitivity.

Determining the factors that influence and can help predict energetic material sensitivity has long been a challenge in the explosives community. Decades of literature reports identify a multitude of factors both chemical and physical that influence explosive sensitivity; however no unifying theory has been observed. Recent work by our team has demonstrated that the kinetics of "trigger linkages" (i.e., the weakest bonds in the energetic material) showed strong correlations with experimental drop hammer impact sensitivity. These correlations suggest that the simple kinetics of the first bonds to break are good indicators for the reactivity observed in simple handling sensitivity tests. Herein we report the synthesis of derivatives of the explosive pentaerythritol tetranitrate (PETN) in which one, two or three of the nitrate ester functional groups are substituted with an inert group. Experimental and computational studies show that explosive sensitivity correlates well with Q (heat of explosion), due to the change in the number of trigger linkages removed from the starting material. In addition, this correlation appears more significant than other observed chemical or physical effects imparted on the material by different inert functional groups, such as heat of formation, heat of explosion, heat capacity, oxygen balance, and the crystal structure of the material.

Radiolytic degradation of dodecane substituted with common energetic functional groups.

Explosives exist in and are expected to withstand a variety of harsh environments up to and including ionizing radiation, though little is known about the chemical consequences of exposing explosives to an ionizing radiation field. This study focused on the radiation-induced chemical changes to a variety of common energetic functional groups by utilizing a consistent molecular backbone. Dodecane was substituted with azide, nitro, nitrate ester, and nitramine functional groups and γ-irradiated with 60Co in order to study how the functional group degraded along with what the relative stability to ionizing radiation was. Chemical changes were assessed using a combination of analysis techniques including: nuclear magnetic resonance (NMR) spectroscopy, gas chromatography of both the condensed and gas phases, Raman spectroscopy, and Fourier transform infrared (FTIR) spectroscopy. Results revealed that much of the damage to the molecules was on the energetic functional group and often concentrated on the trigger linkage, also known as the weakest bond in the molecule. The general trend from most to least susceptible to radiolytic damage was found to be D-ONO2 → D-N3 → D-NHNO2 → D-NO2. These results also appear to be in line with the relative stability of these functional groups to things such as photolysis, thermolysis, and explosive insults.

Chemical Descriptors for a Large-Scale Study on Drop-Weight Impact Sensitivity of High Explosives.

The drop-weight impact test is an experiment that has been used for nearly 80 years to evaluate handling sensitivity of high explosives. Although the results of this test are known to have large statistical uncertainties, it is one of the most common tests due to its accessibility and modest material requirements. In this paper, we compile a large data set of drop-weight impact sensitivity test results (mainly performed at Los Alamos National Laboratory), along with a compendium of molecular and chemical descriptors for the explosives under test. These data consist of over 500 unique explosives, over 1000 repeat tests, and over 100 descriptors, for a total of about 1500 observations. We use random forest methods to estimate a model of explosive handling sensitivity as a function of chemical and molecular properties of the explosives under test. Our model predicts well across a wide range of explosive types, spanning a broad range of explosive performance and sensitivity. We find that properties related to explosive performance, such as heat of explosion, oxygen balance, and functional group, are highly predictive of explosive handling sensitivity. Yet, models that omit many of these properties still perform well. Our results suggest that there is not one or even several factors that explain explosive handling sensitivity, but that there are many complex, interrelated effects at play.

Modeling Photolytic Decomposition of Energetically Functionalized Dodecanes.

The photolytic stability of explosives and energetic functional groups is of importance for those who regularly handle or are exposed to explosives in typical environmental conditions. This study models the photolytic degradation of dodecane substituted with various energetic functional groups: azide, nitro, nitrate ester, and nitramine. For the studied molecules, it was found that excitons localize on the energetic functional group, no matter where they were initially formed, and thus, the predominant degradation pathway involves the degradation of the energetic functional group. The relative trends for both 4 and 8 eV excitation energies followed with what is expected from the relative stability of the energetic functional groups to thermal and sub-shock degradation. The one notable exception was the azide functional group; more work should be done to further understand the photolytic effects on the azide functional group.

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.

Effects of Low-Level Gamma Radiation on Common Nitroaromatic, Nitramine, and Nitrate Ester Explosives.

The aging of high explosives in an ionizing radiation field is not well understood, and little work has been done in the low dose and low dose rate regime. In this study, four explosives were exposed to low-level gamma irradiation from a 137Cs source: PETN, PATO, and PBX 9501 both with and without the Irganox 1010 stabilizer. Post-irradiation analysis included GC-MS of the headspace gas, SEM of the pellets and powder, NMR spectroscopy, DSC analysis, impact sensitivity tests, and ESD sensitivity tests. Overall, no significant change to the materials was seen for the dose and dose rate explored in this study. A small change in the 1H NMR spectrum of PETN was observed and SEM and ESD results suggest a surface energy change in PATO, but these differences are minor and do not appear to have a substantial impact on the handling safety.

Understanding Explosive Sensitivity with Effective Trigger Linkage Kinetics.

We present a simple linear model for ranking the drop weight impact sensitivity of organic explosives that is based explicitly on chemical kinetics. The model is parameterized to specific heats of explosion, Q, and Arrhenius kinetics for the onset of chemical reactions that are obtained from gas-phase Born-Oppenheimer molecular dynamics simulations for a chemically diverse set of 24 molecules. Reactive molecular dynamics simulations sample all possible decomposition pathways of the molecules with the appropriate probabilities to provide an effective reaction barrier. In addition, the calculations of effective trigger linkage kinetics can be accomplished without prior physical intuition of the most likely decomposition pathways. We found that the specific heat of explosion tends to reduce the effective barrier for decomposition in accordance with the Bell-Evans-Polanyi principle, which accounts naturally for the well-known correlations between explosive performance and sensitivity. Our model indicates that sensitive explosives derive their properties from a combination of weak trigger linkages that react at relatively low temperatures and large specific heats of explosion that further reduce the effective activation energy.

Chemical Evaluation and Performance Characterization of Pentaerythritol Tetranitrate (PETN) under Melt Conditions.

Pentaerythritol tetranitrate (PETN) is a nitrate ester explosive commonly used in commercial detonators. Although its degradation properties have been studied extensively, very little information has been collected on its thermal stability in the molten state due to the fact that its melting point is only ∼20 °C below its onset of decomposition. Furthermore, studies that have been performed on PETN thermal degradation often do not fully characterize or quantify the decomposition products. In this study, we heat PETN to melt temperatures and identify thermal decomposition products, morphology changes, and mass loss by ultrahigh-pressure liquid chromatography coupled to quadrupole time of flight mass spectrometry, scanning electron microscopy, nuclear magnetic resonance spectroscopy, and differential scanning calorimetry. For the first time, we quantify several decomposition products using independently prepared standards and establish the resulting melting point depression after the first melt. We also estimate the amount of decomposition relative to sublimation that we measure through gas evolution and evaluate the performance behavior of the molten material in commercial detonator configurations.

Atom Equivalent Energies for the Rapid Estimation of the Heat of Formation of Explosive Molecules from Density Functional Tight Binding Theory.

Atom equivalent energies have been derived from which the gas-phase heat of formation of explosive molecules can be estimated from fast, semiempirical density functional tight binding total energy calculations. The root-mean-square deviation and maximum deviation of the heats of formation from the experimental values for the set of 45 energetic molecules compiled by Byrd and Rice [ J. Phys. Chem. A, 2006, 110, 1005-1013] are 10.4 and 25.5 kcal/mol, respectively, using 4 atom equivalent energies and 7.4 and 15.0 kcal/mol, respectively, using 7 atom equivalent energies. These errors are around a factor of 2-3 larger than those obtained from density functional theory calculations but are smaller than those obtained from other semiempirical electronic structure methods. Heats of formation calculated with density functional tight binding theory using the 4 and 7 atom equivalent energies, the Byrd and Rice scheme, and the atom pair contribution method for a new set of 531 energetic molecules that contain only carbon, hydrogen, nitrogen, and oxygen are provided.

Analysis of Ignition Sites for the Explosives 3,3'-Diamino-4,4'-azoxyfurazan (DAAF) and 1,3,5,7-Tetranitro-1,3,5,7-tetrazoctane (HMX) Using Crush Gun Impact Testing.

The handling safety characteristics of energetic materials must be measured in order to ensure the safe transport and use of explosives. Drop-weight impact sensitivity measurements are one of the first standardized tests performed for energetics. They utilize a small amount of the explosive sample and a standard weight, which is dropped on the material from various heights to determine its sensitivity. While multiple laboratories have used the impact sensitivity test as an initial screening tool for explosive sensitivity for the past 60 years, variability exists due to the use of different instruments, different methods to determine the initiation, and the scatter commonly associated with less-sensitive explosives. For example, standard explosives such as 1,3,5,7-tetranitro-1,3,5,7-tetrazoctane (HMX) initiate reliably and consistently on the drop-weight impact test, whereas insensitive explosives such as 3,3'-diamino-4,4'-azoxyfurazan (DAAF) exhibit variability in sound levels and the expended material. Herein we investigate the impact sensitivity of DAAF and HMX along with a more detailed investigation of ignition sites using a novel "crush gun" apparatus: a pneumatically powered drop-weight tower with advanced diagnostics, including high-speed visual and infrared cameras. Using this crush gun assembly, the ignition sites in HMX and DAAF were analyzed with respect to the effects of particle size and the presence of a source of grit. The formation of ignition sites was observed in both explosives; however, only HMX showed ignition sites that propagated to a deflagration at lower firing speeds. Finally, the presence of grit particles was shown to increase the occurrence of ignition sites in DAAF at lower firing speeds, though propagation to a full reaction was not observed on the time scale of the test. These results enable a better understanding of how ignition and propagation occurs during the impact testing of DAAF.

Synthesis of Erythritol Tetranitrate Derivatives: Functional Group Tuning of Explosive Sensitivity.

Understanding the factors that affect explosive sensitivity is paramount to the safe handling and development of new explosives molecules. Erythritol tetranitrate (ETN) is an explosive that recently has attracted significant attention in the explosives community because of its ease of synthesis and physical properties. Herein, we report the synthesis of ETN derivatives using azide, nitramine, and nitrate ester functional groups. Impact, spark, and friction sensitivity measurements, computationally calculated explosive properties, and the crystal structure analysis of the ETN derivatives are reported. Mixing explosive functional groups led to changes in the explosive sensitivity, explosive performance as well as physical properties including melting point and physical state at room temperature. Overall, we have demonstrated that combining functional groups can enable the tuning of explosive and physical properties of a molecule. This tunability can potentially aid in the development of new explosives in which characteristics are varied to meet certain specifications.

Increased handling sensitivity of molten erythritol tetranitrate (ETN).

Erythritol tetranitrate (ETN) is a well-studied homemade explosive (HME), which is known for its ability to be melt-cast at a fairly low temperature. We have observed dramatically increased handling sensitivity of ETN in the molten state, using temperature controlled drop-weight impact sensitivity measurements. Impact testing was performed using ERL Type 12 drop hammer equipment using a 2.5 kg weight, a 0.8 kg striker, an anvil and sound detection equipment. Most experiments were performed in the absence of standard grit paper, due to the elevated temperature measurements with a liquid. At room temperature, ETN exhibited an impact sensitivity of 14.7 ± 3.4 cm, which changed to 1.0 ± 0.6 cm in the liquid state at 65 °C. The change in sensitivity in the liquid state was found to be reversible upon solidification, and did not appear to correlate with temperature. Control experiments were performed in the same setup using standard explosives pentaerythritol tetranitrate (PETN) and triacetone triperoxide (TATP). This is the most sensitive material that we have been able to measure using our instrumentation, and indicates that ETN be handled with extreme caution during the melt-casting process.

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.

In Situ Imaging during Compression of Plastic Bonded Explosives for Damage Modeling.

The microstructure of plastic bonded explosives (PBXs) is known to influence behavior during mechanical deformation, but characterizing the microstructure can be challenging. For example, the explosive crystals and binder in formulations such as PBX 9501 do not have sufficient X-ray contrast to obtain three-dimensional data by in situ, absorption contrast imaging. To address this difficulty, we have formulated a series of PBXs using octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX) crystals and low-density binder systems. The binders were hydroxyl-terminated polybutadiene (HTPB) or glycidyl azide polymer (GAP) cured with a commercial blend of acrylic monomers/oligomers. The binder density is approximately half of the HMX, allowing for excellent contrast using in situ X-ray computed tomography (CT) imaging. The samples were imaged during unaxial compression using micro-scale CT in an interrupted in situ modality. The rigidity of the binder was observed to significantly influence fracture, crystal-binder delamination, and flow. Additionally, 2D slices from the segmented 3D images were meshed for finite element simulation of the mesoscale response. At low stiffness, the binder and crystal do not delaminate and the crystals move with the material flow; at high stiffness, marked delamination is noted between the crystals and the binder, leading to very different mechanical properties. Initial model results exhibit qualitatively similar delamination.

High-pressure far-infrared spectroscopic studies of hydrogen bonding in formic acid.

Simple molecules such as HCOOH, or formic acid, are suggested to have played important roles in planetary physics due to their possibility for high pressure and temperature chemistry under impact conditions. In this study, we have investigated the effect of pressure (up to 50 GPa) on H-bonding and reactivity of formic acid using synchrotron far infrared spectroscopy. Based on the pressure-induced changes to H-bond ν(O-H···O) stretching and γ(O-H···O) deformations, we observe significant reorganization of H-bonding network beginning at ~20 GPa. This is in good agreement with reports of symmetrization of H-bonds reported at 16-21 GPa from X-ray diffraction and Raman spectroscopy studies as well as molecular dynamics simulations. With further increase in pressure, beyond 35 GPa, formic acid undergoes a polymerization process that is complete beyond 45 GPa. Remarkably, upon decompression, the polymeric phase reverts to the crystalline high-pressure phase at 8 GPa.

Role of solvent-oxygen ion pairs in photooxidation of CdSe nanocrystal quantum dots.

Understanding the mechanisms for photodegradation of nanocrystal quantum dots is an important step toward their application in real-world technologies. A usual assumption is that photochemical modifications in nanocrystals, such as their photooxidation, are triggered by absorption of a photon in the dot itself. Here, we demonstrate that, contrary to this commonly accepted picture, nanocrystal oxidation can be initiated by photoexcitation of solvent-oxygen ion pairs that relax to produce singlet oxygen, which then reacts with the nanocrystals. We make this conclusion on the basis of photolysis studies of solutions of CdSe nanocrystals. Our measurements indicate a sharp spectral onset for photooxidation, which depends on solvent identity and is 4.8 eV for hexane and 3.4 eV for toluene. Importantly, the photooxidation onset correlates with the position of a new optical absorption feature, which develops in a neat solvent upon its exposure to oxygen. This provides direct evidence that nanocrystal photooxidation is mediated by excitation of solvent-oxygen pairs and suggests that the stability of the nanocrystals is defined by not only the properties of their surfaces (as has been commonly believed) but also the properties of their environment, that is, of the surrounding solvent or matrix.

Effect of organic passivation on photoinduced electron transfer across the quantum dot/TiO2 interface.

We report a study of the internal quantum efficiency (IQE) of CdSe quantum-dot (QD)-sensitized solar cells prepared by direct adsorption of pre-synthesized QDs, passivated with either tri-n-octylphosphine oxide (TOPO) or n-butylamine (BA), onto a nanocrystalline TiO(2) film.

CdSe quantum-dot-sensitized solar cell with ∼100% internal quantum efficiency.

We have constructed and studied photoelectrochemical solar cells (PECs) consisting of a photoanode prepared by direct deposition of independently synthesized CdSe nanocrystal quantum dots (NQDs) onto a nanocrystalline TiO(2) film (NQD/TiO(2)), aqueous Na(2)S or Li(2)S electrolyte, and a Pt counter electrode. We show that light harvesting efficiency (LHE) of the NQD/TiO(2) photoanode is significantly enhanced when the NQD surface passivation is changed from tri-n-octylphosphine oxide (TOPO) to 4-butylamine (BA). In the PEC the use of NQDs with a shorter passivating ligand, BA, leads to a significant enhancement in both the electron injection efficiency at the NQD/TiO(2) interface and charge collection efficiency at the NQD/electrolyte interface, with the latter attributed mostly to a more efficient diffusion of the electrolyte through the pores of the photoanode. We show that by utilizing BA-capped NQDs and aqueous Li(2)S as an electrolyte, it is possible to achieve ∼100% internal quantum efficiency of photon-to-electron conversion, matching the performance of dye-sensitized solar cells.

Concerted proton-electron transfer in a ruthenium terpyridyl-benzoate system with a large separation between the redox and basic sites.

To understand how the separation between the electron and proton-accepting sites affects proton-coupled electron transfer (PCET) reactivity, we have prepared ruthenium complexes with 4'-(4-carboxyphenyl)terpyridine ligands, and studied reactivity with hydrogen atom donors (H-X). Ru(II)(pydic)(tpy-PhCOOH) (Ru(II)PhCOOH), was synthesized in one pot from [(p-cymene)RuCl(2)](2), sodium 4'-(4-carboxyphenyl)-2,2':6',2''-terpyridine ([Na(+)]tpy-PhCOO(-)), and disodium pyridine-2,6-dicarboxylate (Na(2)pydic). Ru(II)PhCOOH plus (n)Bu(4)NOH in DMF yields the deprotonated Ru(II) complex, (n)Bu(4)N[Ru(II)(pydic)(tpy-PhCOO)] (Ru(II)PhCOO(-)). The Ru(III) complex (Ru(III)PhCOO) has been isolated by one-electron oxidation of Ru(II)PhCOO(-) with triarylaminium radical cations (NAr(3)(*+)). The bond dissociation free energy (BDFE) of the O-H bond in Ru(II)PhCOOH is calculated from pK(a) and E(1/2) measurements as 87 kcal mol(-1), making Ru(III)PhCOO a strong hydrogen atom acceptor. There are 10 bonds and ca. 11.2 A separating the metal from the carboxylate basic site in Ru(III)PhCOO. Even with this separation, Ru(III)PhCOO oxidizes the hydrogen atom donor TEMPOH (BDFE = 66.5 kcal mol(-1), DeltaG(o)(rxn) = -21 kcal mol(-1)) by removal of an electron and a proton to form Ru(II)PhCOOH and TEMPO radical in a concerted proton-electron transfer (CPET) process. The second order rate constant for this reaction is (1.1 +/- 0.1) x 10(5) M(-1) s(-1) with k(H)/k(D) = 2.1 +/- 0.2, similar to the observed kinetics in an analogous system without the phenyl spacer, Ru(III)(pydic)(tpy-COO(-)) (Ru(III)COO). In contrast, hydrogen transfer from 2,6-di-tert-butyl-p-methoxyphenol [(t)Bu(2)(OMe)ArOH] to Ru(III)PhCOO is several orders of magnitude slower than the analogous reaction with Ru(III)COO.

Trends in ground-state entropies for transition metal based hydrogen atom transfer reactions.

Reported herein are thermochemical studies of hydrogen atom transfer (HAT) reactions involving transition metal H-atom donors M(II)LH and oxyl radicals. [Fe(II)(H(2)bip)(3)](2+), [Fe(II)(H(2)bim)(3)](2+), [Co(II)(H(2)bim)(3)](2+), and Ru(II)(acac)(2)(py-imH) [H(2)bip = 2,2'-bi-1,4,5,6-tetrahydropyrimidine, H(2)bim = 2,2'-bi-imidazoline, acac = 2,4-pentandionato, py-imH = 2-(2'-pyridyl)imidazole)] each react with TEMPO (2,2,6,6-tetramethyl-1-piperidinoxyl) or (t)Bu(3)PhO(*) (2,4,6-tri-tert-butylphenoxyl) to give the deprotonated, oxidized metal complex M(III)L and TEMPOH or (t)Bu(3)PhOH. Solution equilibrium measurements for the reaction of [Co(II)(H(2)bim)(3)](2+) with TEMPO show a large, negative ground-state entropy for hydrogen atom transfer, -41 +/- 2 cal mol(-1) K(-1). This is even more negative than the DeltaS(o)(HAT) = -30 +/- 2 cal mol(-1) K(-1) for the two iron complexes and the DeltaS(o)(HAT) for Ru(II)(acac)(2)(py-imH) + TEMPO, 4.9 +/- 1.1 cal mol(-1) K(-1), as reported earlier. Calorimetric measurements quantitatively confirm the enthalpy of reaction for [Fe(II)(H(2)bip)(3)](2+) + TEMPO, thus also confirming DeltaS(o)(HAT). Calorimetry on TEMPOH + (t)Bu(3)PhO(*) gives DeltaH(o)(HAT) = -11.2 +/- 0.5 kcal mol(-1) which matches the enthalpy predicted from the difference in literature solution BDEs. A brief evaluation of the literature thermochemistry of TEMPOH and (t)Bu(3)PhOH supports the common assumption that DeltaS(o)(HAT) approximately 0 for HAT reactions of organic and small gas-phase molecules. However, this assumption does not hold for transition metal based HAT reactions. The trend in magnitude of |DeltaS(o)(HAT)| for reactions with TEMPO, Ru(II)(acac)(2)(py-imH) << [Fe(II)(H(2)bip)(3)](2+) = [Fe(II)(H(2)bim)(3)](2+) < [Co(II)(H(2)bim)(3)](2+), is surprisingly well predicted by the trends for electron transfer half-reaction entropies, DeltaS(o)(ET), in aprotic solvents. This is because both DeltaS(o)(ET) and DeltaS(o)(HAT) have substantial contributions from vibrational entropy, which varies significantly with the metal center involved. The close connection between DeltaS(o)(HAT) and DeltaS(o)(ET) provides an important link between these two fields and provides a starting point from which to predict which HAT systems will have important ground-state entropy effects.