Computational Design of Novel Energetic Materials: Dinitro-bis

Mar 24, 2015 - Energetic Materials Center, Lawrence Livermore National ... The quantum-chemical computational methods are used to design a new highly ...
1 downloads 0 Views 1MB Size
Subscriber access provided by SETON HALL UNIV

Article

Computational Design of Novel Energetic Materials: Dinitro-Bis-Triazolo-Tetrazine (DNBTT) Roman V. Tsyshevsky, Philip F. Pagoria, and Maija M. Kuklja J. Phys. Chem. C, Just Accepted Manuscript • DOI: 10.1021/acs.jpcc.5b01086 • Publication Date (Web): 24 Mar 2015 Downloaded from http://pubs.acs.org on March 28, 2015

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

The Journal of Physical Chemistry C is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Computational Design of Novel Energetic Materials: Dinitro-Bis-Triazolo-Tetrazine (DNBTT) Roman V. Tsyshevsky1, Philip Pagoria2, and Maija M. Kuklja1,* 1

Materials Science and Engineering Department, University of Maryland College Park, College Park, MD, 20742 2 Energetic Materials Center, Lawrence Livermore National Laboratory, 7000 East Ave., L-282, Livermore, CA 94550 Abstract The quantum-chemical computational methods are used to design a new highly energetic

heterocyclic molecule DNBTT, 2,7-dinitro-4H,9H-bis([1,2,4]triazolo)[1,5-b:1',5'-e][1,2,4,5] tetrazine. We analyze and predict the structure and a range of its properties. DNBTT has high energy content (i.e., high performance) and exhibits a high stability (i.e., low sensitivity) due to two triazole rings connected via a central tetrazine ring. A relatively high activation barrier needs to be overcame to trigger the thermal decomposition of DNBTT, indicating that the sensitivity of DNBTT is on par with (or better than) the benchmark stability of TATB. Additionally, decomposition chemistry can be controlled optically or electronically. Thus, it is expected that DNBTT is an excellent candidate energetic substance with low sensitivity and high performance and is attractive on its own or as a component in composite energetic formulations to improve their properties.

Key words: 2,7-dinitro-4H,9H-bis([1,2,4]triazolo)[1,5-b:1',5'-e][1,2,4,5]tetrazine, heterocyclic high explosives, sensitivity to detonation initiation, decomposition kinetics, density functional theory, BNFF, BNFF-1, and ANFF-1

*

Author to whom correspondence should be addressed, e-mail: [email protected] and [email protected];

phone: 703-292-4940.

ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 2 of 27

I. Introduction Search for new energetic materials that have a stunningly wide range of applications has been long a stubborn challenge for researchers and engineers1. New energetics must satisfy specific requirements, such as an increased performance, reliably high stability to external stimuli, costefficiency and ease of synthesis, be environmentally benign, and be safe for handling and transportation. Among the most important and probably the most challenging prerequisites is low sensitivity of energetic materials to detonation initiation, dictating that explosive decomposition chemistry will be initiated only on demand and not accidently in response to unintended, mild external stimuli. During last decade, the attention of researchers has drifted from widely used nitroether-, nitramine-, and nitroaromatic-based explosives to nitrogen-rich heterocyclic compounds2,3,4,5,6,7. Despite significant efforts and many suggested empirical correlations between properties and structural elements in energetic materials, clear understanding of the relationships between thermal stability of materials, their molecular and crystalline structures, and chemical properties remain unachieved8. The micro-scale theory of decomposition of energetic materials has yet to be developed. A typical approach to design novel high energy density materials is largely of Edisonian nature and involves sophisticated synthesis procedures combined with extensive sensitivity characterization tests. Such empirical explorations are time and effort consuming and often relatively expensive while the successful outcomes are never guaranteed. Thus, many novel energetic materials were attained only to be rejected because they turned out to be so dangerously sensitive that it was problematic to use and, sometimes, even to characterize them9. With tremendous progress in computer architectures and processors’ speed, computational modeling and simulation techniques gradually become a powerful and practical tool to analyze and even predict physicochemical properties of explosive materials including the geometry and electronic structure,10,11,12,13 thermal stability,14,15,16,17,18,19 and their response to mechanical impact20,21,22,23,24,25,26,27. Hence, a close collaboration between synthetic chemists, researchers who perform advanced characterization tests of materials, and those who analyze materials with theoretical and computational quantum-chemical methods prove to be particularly effective1,28. In this article, we report results of such a collaboration to design a potential energetic molecule 2,7-dinitro-4H,9H-bis([1,2,4]triazolo)[1,5-b:1',5'-e][1,2,4,5]tetrazine (DNBTT) with superior properties prior to its synthesis and experimental characterization. The molecular 1 ACS Paragon Plus Environment

Page 3 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

structure and a range of properties (including the formation enthalpy, oxygen balance, decomposition kinetics, optical absorption spectrum, electronic affinity, and ionization potential) are obtained and carefully analyzed by means of quantum-chemical computational methods. Our study concludes that the DNBTT is an outstanding candidate new energetic molecule and the material composed out of these molecules will exhibit high performance and low sensitivity, the properties most desirable for high energy density solids.

II. Computational Details Decomposition reactions in the gas phase were studied within density functional theory (DFT)29,30 with M0631 functionals as implemented in the GAUSSIAN code32. All calculations were performed by using valence split double-zeta 6-31+G(2df,p) basis set. Minimal energy paths were investigated by conducting intrinsic reaction coordinate computations using the Hessian-based Predictor-Corrector integrator algorithm33,34 for each transition state. Reaction rates were calculated by applying transition state theory (TST)35 for reactions with well-defined transition states and variational TST36 for reactions with no barriers as detailed elsewhere15,16,17. It has been recently claimed that M06-based functionals show an outstanding accuracy in predictions of various physical and chemical parameters31. Hence, we carried out a series of test calculations that assess the ability of M06 functional with 6-31+G(2df,p) basis set to predict formation enthalpies (∆fH0298) and C-N bond dissociation enthalpies (BDE) of some C-nitro compounds as well as their ionization potentials (IP) and electronic affinities (EA). These tests served as a quantitative calibration of our approach. Table 1 shows that the calculated ∆fH0298 and BDE values are in good agreement with experiment, with absolute deviations not exceeding ~1-3 kcal/mol, which is within best accuracy attainable for advanced DFT methods. The calculated ionization potentials (IP) and electron affinities (EA) collected in Table 1 are also in respectful agreement with experiment and discrepancies do not exceed a fraction of eV. Table S1 of Supplementary Information shows that M06 predicts well formation enthalpies, dissociation enthalpies of C-NO2 bond of C1-C3 mono nitroalkanes, and activation barriers of the HONO elimination reaction from C2-C3 mononitroalkanes. For example, absolute deviations of calculated activation barriers from experimental data do not exceed 1 kcal/mol (Table S1). Therefore, one should conclude that the M06 functional with 6-31+G(2df,p) basis set can be

2 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 27

reliably used to explore the decomposition mechanisms of DNBTT molecule and its thermodynamic properties and electronic properties. Table 1. Calculated formation enthalpies (∆fH0,298, kcal/mol), enthalpies of the C-N bond dissociation (BDE, kcal/mol), electronic affinities (EA, eV), ionization potentials (IP, eV), and vertical ionization potentials (IPvert, eV) for a series of nitro molecules. Molecule

∆fH0, 298

BDE (kcal/mol)

EA

IP

IPvert

NO2

5.8 [7.9]37a

--

2.26 [2.27]38,39

9.50 [9.59]40,41

11.49 [11.23]42,43

CH3NO2

-20.9 [-17.8]44,45

61.4 [60.8]46

0.16 [0.17]47,39

11.29 [11.28]48,49

11.70 [11.80]50,49

C2H5NO2

-27.1 [-24.4]51

61.8 [60.7]46

0.15

10.59 [11.92]52,53

11.24 [11.02]54,49

C6H5NO2

16.9 [16.4]55,45

73.9 [70.7]56

1.10 [1.0]57,39

9.78 [9.8]58,43

9.95 [9.93]59,49

TNT

5.1 [5.8]60,45

64.3 [58.961, 58.062]

2.3

10.43 [10.59]63,49

10.70

a

Corresponding experimental data are in parentheses

III. Results and Discussion 3.1 The Structure Our approach to theoretically design and computationally model a novel energetic material with superior properties capitalizes on many years of experience coupled with interdisciplinary expertise in the field. In our recent study, we performed a detailed holistic analysis of the structure, stability, and decomposition mechanisms of oxadiazole-based BNFF28, BNFF-164 and ANFF-165 (Fig. 1), novel improved energetic materials that were recently synthesized and characterized in Lawrence Livermore National Laboratory. High density (~1.8−1.9 g/cm3) together with low sensitivity to mechanical impact and low melting point (85−110 °C) make these materials attractive melt castable explosives and an excellent energetic ingredient of modern explosive composites. It was established that these three relevant materials are comparatively thermally stable, as the dissociation of the critical C-NO2 chemical bond requires overcome a high energy barrier of ~60-70 kcal/mol, similar to the stability benchmark TATB. However, in a series of gas-phase64 and solid-state28 calculations, we also discovered another plausible decomposition channel, common for these materials, the cleavage of the oxadiazole rings, which requires a lower energy, ~45-50 kcal/mol, and potentially can prompt an increase of sensitivity to detonation initiation. 3 ACS Paragon Plus Environment

Page 5 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

In designing a new material, we formulate specific criteria that it has to meet, in particular: (i) the energetic molecule should be comprised of several nitrogen and/or oxygen fused hetero rings to ensure high energy content; (ii) the nitro group has to be attached to a carbon atom to ensure high thermal stability, as the C-NO2 bond has the highest dissociation energy as compared to O-NO2 and N-NO2 bonds; and (iii) the molecular structure should serve to additionally stabilize heterocyclic rings to prevent their dissociation as the primary decomposition step. Our search resulted in the following compound 2,7-dinitro-4H,9H-bis([1,2,4]triazolo)[1,5b:1',5'-e][1,2,4,5]tetrazine depicted in Fig. 2. The DNBTT molecule belongs to C2 symmetry point group and has an inversion as symmetry operation. It has a planar structure with two triazole rings connected via a central tetrazine ring. Two nitro groups are inversely symmetrically attached to the carbon atoms of the triazole rings (Fig. 2). This molecule possesses a zero dipole moment due to symmetry, which will additionally stabilize the molecule. Quantum-chemical calculations demonstrated that the existence of the molecular dipole moment may be a reason for a crucial drop in thermal stability of energetic materials due to polar surface effect66,67. For example, among chemically related C-nitrocompounds, polar DADNE is found to be more sensitive than nonpolar TATB20,21,68,69,70, and among stoichiometrically similar nitramines, polar δ-HMX is more sensitive than nonpolar βHMX16,17,66,67.

Figure 1. Sketched structures of BNFF, BNFF-1 and ANFF-1 molecules

4 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

H O N N N N + N N N N O N H

Page 6 of 27

O + N O

Figure 2. The suggested molecular structure of DNBTT (bonds are in Å).

3.2 Formation Enthalpy and Oxygen Balance Our calculations demonstrate that the DNBTT molecule has a high positive formation enthalpy of 177.1 kcal/mol. The formation enthalpy of a nitro compound is directly related to enthalpy (heat) of combustion according to Hess’ law: ∆H0combustion = Σ∆H0f(P) - ∆H0f(NC) - ∆H0f(O2),

(1)

where Σ∆H0f(P) is a sum of formation enthalpies of products, which usually are H2O, CO2, CO, N2; ∆H0f(NC) is the formation enthalpy of the corresponding nitro compound; ∆H0f(O2) – enthalpy of formation of O2, which is equal to zero. Taking into consideration that Σ∆H0f(P) is in most cases a constant value, the higher the formation enthalpy of a nitro compound, the higher the heat released during its combustion, as it is seen from Eq. 1. In other words, the formation enthalpy is a measure of thermal energy stored in the material (and released in explosion), or the materials performance. The formation enthalpy is also indirectly related to the velocity of detonation of explosive materials71.

5 ACS Paragon Plus Environment

Page 7 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

The calculated formation enthalpy of DNBTT (∆H0f(DNBTT)=177.1 kcal/mol)

is

appreciably higher than formation enthalpy of widely used materials, HMX (∆H0f(HMX)=32.7 kcal/mol)72, TNT (∆H0f(TNT)=5.1 kcal/mol)72, and recently synthesized triazole-based explosives (for more details see, for example, ref.4). Another important parameter related to the molecular formula of a compound is an oxygen balance (OB) which indicates the degree to which an explosive can be oxidized:

OB = [16(Z-2X+Y/2)/Mw] × 100%,

(2)

where X, Y, and Z represent a number of atoms of carbon, hydrogen, and oxygen in the molecule, and Mw stands for the molecular weight. Widely used secondary explosives, such as PETN, HMX, TATB and TNT, are considered oxygen deficient, which means that they have negative oxygen balances of -10.1, -21.6, -55.8, and -74.0%, respectively. Negative OB values show that the oxygen content in all these compounds is insufficient to react with the available carbon and hydrogen. DNBTT is also oxygen deficient material with OB = -31.5%.

3.3 Decomposition Mechanisms Now we will explore fragmentation mechanisms of the DNBTT molecule to establish its stability against thermal excitation. In modeling the thermal decomposition of the DNBTT molecule, we studied five most plausible decomposition channels (Fig. 3, 4): I) a cyclization of the nitro-triazole ring into the oxazete N-oxide isomer (CNON cyclization), II) a homolytic cleavage of the C-NO2 bond, III) an opening of the outer triazole ring (OR cleavage) via a scission of the N-N bond accompanied by the C-NO2 bond homolysis and the detachment of NO2, IV) a nitro-nitrite isomerization (CONO) followed by the NO loss, and V) an opening of the central tetrazine ring (CR cleavage). Structures of some reaction products and intermediates are collected in Fig. 5a-c, and Cartesian coordinates of all intermediates, products and transitions states are presented in Supplementary materials.

6 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

I

H N

Page 8 of 27

N

N N+ N+ O N O H H N

II H N N H

N N N

DNBTT

N

N N N H

O N O

III

H N N H

IV1

H N

C

+

NO2

N C N

+

N

N

N N N H

NO2

N O IV 2 O

H N

N

N N N H

O

+

NO

Figure 3. Schematic representation of decomposition mechanisms of DNBTT.

7 ACS Paragon Plus Environment

Page 9 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3.2.1 Isomerization into oxazete N-oxide isomer (path I, Fig. 3) The decomposition mechanism through a formation of the 1,2-oxazete 2-oxide structure, initially proposed for high-temperature pyrolysis of β-nitrostyrenes,73 has been recently reported for nitroalkenes74, DADNE75, and BNFF28and claimed to be the primary initiation reaction75. The calculated activation barriers of this isomerization channel, although are fairly low in energy, differ significantly; and they are 48.1 kcal/mol for nitroethylene74, 27.1 kcal/mol for DADNE75, and 57.4 kcal/mol for BNFF28. The similar isomerization of DNBTT requires significantly higher energy (~90 kcal/mol, Table 2), which suggests to rule out such a reorganization of the molecule as a primary or the most favorable initiation reaction. The structure of oxazete N-oxide isomer of DNBTT is depicted in Figure 5a.

3.2.2 Homolysis of C-NO2 bond (path II, Fig. 3) The C-NO2 bond dissociation is considered a canonical initiation reaction of nitro compounds. The calculated enthalpy of the C-N bond cleavage in the DNBTT molecule (69.6 kcal/mol, Table 2) is in good agreement with theoretical and experimental estimates previously reported for C-nitro derivatives of triazole (~67 kcal/mol)76, dinitro isomers of diazole (~65-70 kcal/mol)77, some nitroaromatic compounds, including nitrobenzene (~70 kcal/mol)78, TATB (60-70 kcal/mol)61,78,79,80 and DADNE (67.0-72 kcal/mol)20,68,81. We note that the cleavage of CNO2 in DNBTT requires a higher energy (by ~6-8 kcal/mol) than the activation of this channel in the oxadiazole-based molecules BNFF-1 (63.1 kcal/mol)64, ANFF-1 (61.5 kcal/mol)64, and BNFF (62.3 kcal/mol)28, alluding to a higher stability of DNBTT.

3.2.3 Cleavage of outer triazole ring (path III, Fig. 3) The cleavage of oxadiazole rings was found to require an appreciably lower energy (~50 kcal/mol) than the dissociation of C-NO2 bonds (~60-65 kcal/mol) and therefore it was suggested as a dominating decomposition mechanism of the gas-phase decomposition of BNFF-164, ANFF164 and BNFF28. Another theoretical study19 proposed that the cleavage of the tetrazole ring and a subsequent elimination of N2, with the activation energy of 37-45 kcal/mol, will be the dominant process in the thermal decomposition of bistetrazole. According to our calculations, an opening of the outer triazole ring in DNBTT proceeds via a rupture of the N-N bond (Fig. 3) and subsequent homolysis of the C-NO2 bond. The calculated 8 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 10 of 27

activation barrier of the outer triazole ring opening, 63.6 kcal/mol (Table 2), is noticeably higher than that in previously reported data for the decomposition of BNFF (49.5 kcal/mol)28, its derivatives (BNFF-1 (50.7 kcal/mol)64 and ANFF-1 (51.8 kcal/mol)64), and bistetrazole (45.1 kcal/mol)19.

O O N

V2

N N

N

N N

O

+

N O

N

NH

HN H N

-

O O N N N + + N N N N N O O N H DNBTT

O O N

V3

N N

N

N N

+

O N O

N

HN

HN

V1 V4

O O N N HN

N N

N N NH

N

O

N N

N

N N NH

N O

N

HN

N O

O

+

NO2

O O N

V5

N N

N N

N

+

NO2

N

NH HN O N O

V6

N N

N N

N

NH

N

HN

O N O

V7

... +

NO

... +

NO

O O N

V8

N

N N

N O V9

N N O

NH N HN Figure 4. Schematic representation of decomposition mechanisms of DNBTT proceeding via opening of central tetrazine ring. 9 ACS Paragon Plus Environment

Page 11 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 5. The structures of a) oxazete N-oxide, b) CONO isomers of DNBTT, and c) intermediate, formed in the course of central tetrazine ring opening.

3.2.4 CONO isomerization and subsequent NO loss (path IV, Fig. 3) The CONO isomerization is often suggested as a competing or a slow background initiation reaction in nitro compounds. The nitro-nitrite isomerization proceeds via a pseudo rotation of a nitro group accompanied by breaking the C-NO2 bond and creating the C-ONO bond (path IV1, Fig. 3). The formed nitrite isomers of DNBTT are shown in Fig. 5b. The CONO-isomerization in DNBTT requires 55.6 kcal/mol, which is consistent with the barrier of the CONO formation in 10 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 12 of 27

C-nitro isomers of triazole (60–65 kcal/mol)76, BNFF (54.2 kcal/mol)28, BNFF-1 (53.2 kcal/mol)64, ANFF-1 (60.3 kcal/mol)64, DADNE (66.4 kcal/mol81, 59 kcal/mol82), nitroethylene (57.9 kcal/mol)83, and TNT (54.9 kcal/mol)62. The CONO isomerization in DNBTT is an exothermic reaction, releasing 7.4 kcal/mol of heat (Table 2). The following NO loss from the CONO isomer of DNBTT (path IV2, Fig. 3) requires only 16.7 kcal/mol (Table 2). 3.2.5 Fragmentation of DNBTT via opening of central tetrazine ring (path V, Fig. 4) The final mechanism of the unimolecular DNBTT decomposition, we consider here, involves opening of the central tetrazine ring during the first step, followed by subsequent fragmentation mechanisms. The opening of the central tetrazine ring (path V1, Fig. 4) proceeds via a cleavage of one N-N bond, accompanied by a rotation of the tetrazine ring about the second N-N bond. The resulting intermediate with the opened tetrazine ring (Fig. 5c) rests 31.0 kcal/mol above the DNBTT equilibrium state. The tetrazine ring opening has the lowest activation barrier (53.9 kcal/mol, Table 2) among all studied mechanisms. Further, we simulated several secondary fragmentation channels, including a cleavage of the central C-N (path V2, Fig. 4) and N-N (path V3, Fig. 4) bonds, C-NO2 bond homolysis (paths V4, V5, Fig. 4), and CONO isomerization mechanism (paths V6, V8, Fig. 4) together with subsequent NO loss reactions (paths V7, V9, Fig. 4). The secondary fragmentation channels were explored to reveal whether or not the primary tetrazine ring opening reaction can be considered as the rate limiting step. An analysis of the secondary reactions, illustrated in Fig. 4 and Table 2, shows that the dissociation of the central C-N bond (path V2) and homolytic loss of NO2 (paths V4, V5) require noticeably higher energies (~63-73 kcal/mol, Table 2) than the primary reaction of the central ring opening (53.9 kcal/mol). The cleavage of the central N-N bond (path V3, Fig. 4) has the lowest activation barrier, 34.7 kcal/mol (Table 2). The CONO isomerization of the DNBTT isomer with the opened tetrazine ring (paths V6, V8, Fig. 4) requires much lower energies (45.9 and 37.5 kcal/mol, Table 2, Fig. 6), as compared to the primary CONO isomerization (55.7 kcal/mol, Table 2, path IV1, Fig. 6). The subsequent NO loss (path V9, Fig. 4) from the CONO isomer requires 16.0 kcal/mol (Table 2)84.

11 ACS Paragon Plus Environment

Page 13 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

3.4 Decomposition Kinetics The homolytic NO2 loss (path II, Fig. 3) is to be considered the predominant fragmentation channel (Fig. 7a) at the temperature above ~650 K (Fig. 7a) due to the highest value of preexponential factor log(A, s-1) = 17.9 (Table 2). The CNNO cyclization mechanism (path I, Fig. 3) exhibits the highest barrier among all studied mechanisms (Table 2, Fig. 6) and, hence, would proceed at the lowest rate (Fig. 7a). Although the cleavage of the outer triazole ring (path III, Fig. 3) has a noticeably lower activation barrier (63.6 kcal/mol, Table 2, Fig. 6) than the energy of the C-NO2 bond homolysis (69.6 kcal/mol, Table 2, Fig. 6), this fragmentation channel will unlikely be able to compete with the NO2 loss due to its lower pre-exponential factor (log(A, s-1) = 15.0, Table 2). The CONO isomerization (path IV1, Fig. 3) has even lower pre-exponential factor (log(A, s-1) = 13.5, Table 2) than the outer triazole ring cleavage. Nevertheless, the CONO isomerization will be favored over the C-NO2 bond cleavage at low temperatures (Fig. 7a) due to the relatively low activation barrier height (55.6 kcal/mol, Table 2, Fig. 6). The cleavage of the central tetrazine ring (path V1, Fig. 4) has the lowest activation barrier (53.9 kcal/mol, Table 2, Fig. 6). As a result, this mechanism has the highest rate (Fig. 7a) in the temperature range below ~650 K in spite of quite moderate value of pre-exponential factor log(A, s-1) = 13.7 (Table 2). However, Fig. 6 shows that all studied secondary reactions have higher overall activation barriers than the primary cleavage of the central tetrazine ring. For example, the overall energy of the central C-N bond cleavage (path V2, Fig. 4) and homolytic NO2 loss reaction (paths V4, V5, Fig. 4) exceeds 90 kcal/mol (Fig. 6, Table 2). Such high activation barriers dictate to eliminate these mechanisms from further consideration. The lowest overall activation barriers of 65.7 and 68.5 kcal/mol (Fig. 6, Table 2) were obtained for the secondary steps that involve the central N-N bond rupture (paths V3, Fig. 4) and CONO isomerization (paths V8, Fig. 4), respectively. Figure 7b reveals that the latter two mechanisms have lower rates than the primary cleavage of tetrazine rings and hence they should be considered as the rate limited steps for the stepwise channels V1-V3 and V1-V8- V9 (Fig. 4). Summing up, the primary C-NO2 homolysis (paths II, Fig. 3) and NO loss via CONO isomerization (paths IV1, Fig. 3) will dominate other possible decomposition channels with the NO2 loss being most favorable at higher temperatures (Fig. 7b).

12 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 27

Figure 6. Energy diagram showing main decomposition channels of DNBTT.

Table 2. Calculated activation barriers (E, kcal/mol), reaction energies (given in parentheses in kcal/mol) zero-point energy corrected barriers (EZPE, kcal/mol), activation enthalpies (∆H≠298, kcal/mol), and pre-exponential factors (log(A, s−1)) of DNBTT decomposition reactions. Reaction I CNNO cyclization II NO2 loss III OR cleavage 1 CONO IV 2 NO loss 1 CR cleavage 2 C-N bond cleavage 3 N-N bond cleavage 4 NO2 loss 5 NO2 loss V 6 CONO 7 NO loss 8 CONO 9 NO loss

E 92.5 (87.1) 72.7 (72.7) 67.0 (32.4) 57.9 (-6.4) 19.0 (19.0) 57.2 (32.9) 76.5 (76.5) 36.8 (24.2) 69.6 (69.6) 66.3 (66.3) 48.0 (-10.8) 39.3 (-18.4) 18.6 (18.6)

EZPE 90.1 (85.2) 68.6 (68.6) 62.8 (27.1) 55.6 (-7.8) 15.9 (15.9) 53.9 (30.4) 72.8 (72.8) 34.6 (20.8) 65.2 (65.2) 62.3 (62.3) 45.9 (-12.3) 38.2 (-19.8) 15.3 (15.3)

∆H≠298 89.8 (85.2) 69.6 (69.6) 63.6 (28.5) 55.6 (-7.4) 16.7 (16.7) 53.9 (31.0) 73.1 (73.1) 34.7 (21.3) 65.9 (65.9) 62.9 (62.9) 45.9 (-11.9) 37.5 (-19.3) 16.0 (16.0)

log A 12.7 17.9 15.0 13.5 13.7 15.6 14.9 14.2 -

13 ACS Paragon Plus Environment

Page 15 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

Figure 7. Reaction rates of the simulated decomposition channels of DNBTT.

3.5 Electronic Properties and Optical Absorption Electronic properties serve as a measure of electronic stability of the material against photo stimulated decomposition85,86,87,88,89,90 or decomposition induced by electromagnetic field91. The calculated HOMO-LUMO gap of the DNBTT molecule is 4.97 (Table 3) eV, which is close to that of BNFF (4.87 eV)28 and TATB (5.0 eV)28. The energy gap represents a typical value for organic molecular high explosive materials and indicates that DNBTT crystals will likely to be wide-gap dielectrics. The highest occupied molecular orbital (HOMO) orbital is mainly formed from 2pz atomic functions of tetrazine nitrogen atoms and an additional small 14 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 16 of 27

contribution from 2pz atomic functions of nitrogen and carbon atoms of triazole rings and oxygen atoms (Fig. 8a). The lowest unoccupied molecular orbital (LUMO) orbital is mainly formed from 2pz atomic functions of oxygen and nitrogen atoms of nitro groups and carbon atoms of triazole rings with an additional relatively moderate contribution from 2pz atomic functions of triazole nitrogen atoms (Fig. 8a). The energies of optical transitions of DNBTT were calculated using time dependent TDM06/6-31+G(2df,p) approximation to predict the absorption spectrum. The energy of the vertical spin-forbidden transition from the ground state to the lowest triplet state (S0→T1) requires 3.15 eV. The energy of the lowest singlet-singlet (S0→S1) transition is 3.85 eV (Table 3, Fig. 8b). The large difference (1.8 eV) between the HOMO-LUMO gap and S0→T1 energy indicates the existence of tightly bound excitons in DNBTT. Our calculations also show that the UV-VIS absorption spectrum of DNBTT (Fig. 8b) has two absorption bands corresponding to a series of symmetry- and spin- allowed singlet-singlet transitions in the range of 3-7 eV with distinct maximums at 4.25 and 6.60 eV. These features are consistent with the optical spectra of many nitro-compounds11,13. In a practical sense, the obtained absorption spectrum means that once the DNBTT molecule is irradiated by photons with the predicted energies, the molecule will likely to suffer rapid bond dissociation. On the other hand, the chemical bonds can be stimulated selectively, for example, by laser irradiation, implying highly controllable initiation reactions. Such a precision control is not achievable in thermally stimulated decomposition. The ionization potential of DNBTT (9.19 eV, Table 3) is comparable to that of BNFF, RDX, HMX, PETN, TNT (~ 10 eV), and TATB (8.7 eV)28. The high IP indicates that the removal of an electron from those energetic molecules is an energetically costly process. Another useful characteristic relevant to molecule and materials stability is an electronic affinity, which tells whether the molecule is electron-withdrawing or electron-donating. The calculated EA of DNBTT is positive and equal to 2.03 eV (Table 3), which is close to EA of BNFF (2.3 eV) and TNT (2.0 eV) and somewhat higher than EA of PETN (1.7 eV), HMX (1.4 eV), and TATB (1.3)28. The high EA of DNBTT suggests that in an electron-rich environment, the DNBTT will easily trap an electron, gain energy by doing so, and the ion radicals may exhibit different activation barriers for dissociating chemical bonds. For example, a dramatic effect of charged states on decomposition chemistry, illustrated on DADNE69, HMX8,66,67 and PETN92, was found to be manifested in altered dominating chemistry with pronounced 15 ACS Paragon Plus Environment

Page 17 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

exothermicity as well as in significantly reduced activation barriers of bond dissociation reactions. The observation indicates that DNBTT decomposition can be efficiently controlled by an electric current, donor impurity dopants, or electromagnetic fields. More so, DNBTT is expected to serve as an attractive additive ingredient of modern energetic formulations to improve their stability and achieve highly controllable decomposition chemistry.

Table 3. Electronic properties of DNBTT energetic molecules, optical energy gap (HOMOLUMO), optical transitions, electronic affinity, and ionization potential (eV). Property

Energy (eV)

HOMO-LUMO

4.97

S0→T1

3.15

S0→S1

3.85

EAad

2.03

IPad

9.19

16 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 18 of 27

Figure 8. a) Isosurfaces of HOMO and LUMO orbitals of DNBTT and b) schematic of the predicted UV-Vis absorption spectrum of DNBTT.

IV. Summary and Conclusions

The molecular structure of the theoretically designed novel energetic heterocyclic molecule 2,7-dinitro-4H,9H-bis([1,2,4]triazolo)[1,5-b:1',5'-e][1,2,4,5] tetrazine, DNBTT, is proposed. This molecule was constructed based on a detailed analysis of available data on structures and thermal stability of several classes of energetic heterocycles. Potential decomposition mechanisms and kinetics were computationally simulated by using density functional theory combined with transition state theory.

17 ACS Paragon Plus Environment

Page 19 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

DNBTT is a fused heterocyclic molecule, which contains two triazole rings connected via a central tetrazine ring. The fused structure provides an additional stabilization of the hetero rings, which improves the overall thermal stability of the molecule, as evidenced by the modeled decomposition mechanisms. The DNBTT material has yet to be synthesized, and, thus, no data on the crystalline structure are available. Speculating on intermolecular interactions or cohesive binding in DNBTT crystals we suggest that the potential material will be organic molecular crystal in which molecules are held together by van der Waals forces and the electronic density is largely localized on individual molecules. The thermal decomposition of DNBTT is predicted to proceed via two possible channels, the fast NO2 loss (at high T) and slow NO loss via CONO isomerization (at low T). In contrast to the behavior observed on BNFF28, BNFF-164, ANFF-165 and bistetrazole19 based energetic materials, the cleavage of the hetero rings can barely compete with the two dominating dissociation mechanisms in DNBTT. DFT calculations suggest that the activation energy of thermal decomposition of DNBTT is relatively high, comparable to stability benchmark TATB (57-70 kcal/mol) and noticeably higher than the activation energy (~45-50 kcal/mol) of relevant materials of BNFF28, BNFF-164, ANFF-165, recently reported in literature. A high positive electronic affinity and optical absorption energies suggest specific ways to control decomposition chemistry of DNBTT and to increase sensitivity of energetic composite formulations. Overall, this study predicts that DNBTT is an excellent candidate low sensitivity energetic material individually and as an ingredient in composite formulations. This study provides a solid background for further design of new energetic materials and targeted improvements of existing composites.

Supporting Information Supporting Information Available: Details of the structure of a DBNTT molecule, the structures of considered molecules in their transition states during decomposition, and atomic coordinates of the molecule in its ground state and in the transition states for each decomposition reaction. This material is available free of charge via the Internet at http://pubs.acs.org.

18 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 20 of 27

ACKNOWLEDGMENT Work performed under the auspices of the U.S. Department of Energy by the Lawrence Livermore National Laboratory under Contract DE-AC52-07NA27344. This research is also supported in part by ONR (Grant N00014-12-1-0529) and NSF. We used NSF XSEDE resources (Grant DMR-130077) and DOE NERSC resources (Contract DE-AC02-05CH11231). MMK is grateful to the Office of the Director of National Science Foundation for support under the IRD program. Any appearance of findings, conclusions, or recommendations expressed in this material are those of the authors and do not necessarily reflect the views of NSF.

19 ACS Paragon Plus Environment

Page 21 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

References: 1

Fried, L. E.; Manaa, M. R.; Pagoria, P. F.; Simpson, R. L. Design and Synthesis of Energetic Materials 1. Annu. Rev. Mater. Res. 2001, 31, 291-321. 2 Gao, H., Shreeve, J. M. Azole-Based Energetic Salts. Chem. Rev. 2011, 111, 7377-7436. 3 Göbel, M.; Karaghiosoff, K.; Klapötke, T. M.; Piercey, D. G.; Stierstorfer, J. Nitrotetrazolate-2 N-oxides and the Strategy of N-oxide Introduction. J. Am. Chem. Soc. 2010, 132, 1721617226. 4 Dippold, A. A.; Klapötke, T. M. A Study of Dinitro-bis-1, 2, 4-triazole-1, 1′-diol and Derivatives: Design of High-Performance Insensitive Energetic Materials by the Introduction of N-Oxides. J. Am. Chem. Soc. 2013, 135, 9931-9938. 5 Zhang, Y.; Parrish, D. A.; Shreeve, J. M. Derivatives of 5-nitro-1, 2, 3-2H-Triazole–High Performance Energetic Materials. J. Mater. Chem. A 2013, 1, 585-593. 6 Fischer, N.; Fischer, D.; Klapötke, T. M.; Piercey, D. G.; Stierstorfer, J. Pushing the Limits of Energetic Materials–the Synthesis and Characterization of Dihydroxylammonium 5, 5′bistetrazole-1, 1′-diolate. J. Mater. Chem. 2012, 22, 20418-20422. 7 Klapötke, T. M.; Carles, M. S., Bistetrazoles: Nitrogen-rich, High-performing, Insensitive Energetic Compounds, Chem. of Mater. 2008, 20.11, 3629-3637. 8 Kuklja, M. M. Quantum-Chemical Modeling of Energetic Materials: Chemical Reactions Triggered by Defects, Deformations, and Electronic Excitations. Adv. Quant. Chem. 2014, 68, 71-146. 9 Klapötke, T. M.; Krumm, B.; Ilg, R.; Troegel, D.; Tacke, R. The Sila-Explosives Si (CH2N3)4 and Si (CH2ONO2)4: Silicon Analogues of the Common Explosives Pentaerythrityl Tetraazide, C (CH2N3)4, and Pentaerythritol Tetranitrate, C (CH2ONO2)4. J. Am. Chem. Soc. 2007, 129, 6908-6915. 10 Kuklja, M.M.; Kunz, A.B., Electronic Structure of Molecular Crystals Containing Edge Dislocations. J. Appl. Phys. 2001, 89, 4962-4970. 11 Tsyshevsky, R.; Sharia, O.; Kuklja, M. M. Energies of Electronic Transitions of PETN Molecules and Crystals. J. Phys. Chem. C 2014, 118, 9324–9335. 12 Kuklja M.M.; Kunz, A.B. Compression-Induced Effect on the Electronic Structure of Cyclotrimethylene Trinitramine Containing an Edge Dislocation. J. Appl. Phys. 2000, 87, 2215-2218. 13 Tsyshevsky, R.; Sharia, O.; Kuklja, M. M. Optical Absorption Energies of Molecular Defects in Pentaerythritol Tetranitrate Crystals: Quantum Chemical Modeling. J. Phys. Chem. C 2014, 118, 26530–26542. 14 Kuklja, M. M. Thermal Decomposition of Solid Cyclotrimethylene Trinitramine. J. Phys. Chem. B 2001, 105, 10159-10162. 15 Sharia, O.; Kuklja, M. M. Ab Initio Kinetics of Gas Phase Decomposition Reactions. J. Phys. Chem. A 2010, 114, 12656–12661. 16 Sharia, O.; Kuklja, M. M. Modeling Thermal Decomposition Mechanisms in Gaseous and Crystalline Molecular Materials: Application to β-HMX. J. Phys. Chem. B 2011, 115, 1267712686. 17 Sharia, O.; Kuklja, M. M. Surface-Enhanced Decomposition Kinetics of Molecular Materials Illustrated with Cyclotetramethylene-tetranitramine. J. Phys. Chem. C 2012, 116, 1107711081.

20 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 22 of 27

18

Tsyshevsky, R.; Sharia, O.; Kuklja, M. M. Thermal Decomposition Mechanisms of Nitroesters: Ab Initio Modeling of Pentaerythritol Tetranitrate. J. Phys. Chem. C 2013, 117, 18144−18153. 19 An, Qi.; Liu, W.-G.; Goddard III, W. A.; Cheng, T.; Zybin, S. V.; Xiao, H. Initial Steps of Thermal Decomposition of Dihydroxylammonium 5,5′-bistetrazole-1,1′-diolate Crystals from Quantum Mechanics. J. Phys. Chem. C 2014, 118, 27175−27181. 20 Kuklja, M.M.; Rashkeev, S. N. Shear-strain-induced Chemical Reactivity of Layered Molecular Crystals. Appl. Phys. Lett. 2007, 90, 151913-1-3. 21 Kuklja, M. M.; Rashkeev, S. N. Self-Accelerated Mechanochemistry in Nitroarenes, J. Phys. Chem. Lett. 2010, 1, 363–367 22 Brennan, J. K.; Lísal, M.; Moore, J. D.; Izvekov, S.; Schweigert, I. V.; Larentzos, J. P. CoarseGrain Model Simulations of Nonequilibrium Dynamics in Heterogeneous Materials. J. Phys. Chem. Lett. 2014, 5, 2144–2149. 23 Landerville, A. C.; Oleynik, I. I.; White, C. T. Reactive Molecular Dynamics of Hypervelocity Collisions of PETN Molecules. J. Phys. Chem. A 2009, 113, 12094–12104. 24 Conroy, M. W.; Oleynik, I. I.; Zybin, S. V.; White, C. T. White. First-Principles Anisotropic Constitutive Relationships in β-Cyclotetramethylene Tetranitramine (β-HMX). J. Appl. Phys. 2008, 104, 053506–1–6. 25 Zhou, T. T.; Zybin, S. V.; Liu, Y.; Huang, F. L.; Goddard III, W. A. Anisotropic Shock Sensitivity for β-Octahydro-1,3,5,7-Tetranitro-1,3,5,7-Tetrazocine Energetic Material Under Compressive-Shear Loading from ReaxFF-Lg Reactive Dynamics Simulations. J. Appl. Phys. 2012, 111, 124904–1–11. 26 An, Q.; Liu, Y.; Zybin, S. V.; Kim, H.; Goddard III, W. A. Anisotropic Shock Sensitivity of Cyclotrimethylene Trinitramine (RDX) from Compress-and-Shear Reactive Dynamics. J. Phys. Chem. C 2012, 116, 10198−10206 27 Liu, L. C.; Liu, Y.; Zybin, S. V.; Goddard III, W. A. ReaxFF-lg: Correction of the ReaxFF Reactive Force Field for London Dispersion, with Applications to the Equations of State for Energetic Materials. J. Phys. Chem. A 2011, 115, 11016–11022. 28 Tsyshevsky, R.; Pagoria, P.; Zhang, M.; Racoveanu, A.; DeHope, A.; Parrish, D. A.; Kuklja, M. M. Searching for Low Sensitivity Cast-melt High Energy Density Materials: Synthesis, Characterization, and Decomposition Kinetics of BNFF. J. Phys. Chem. C 2015, 119, 3509−3521. 29 Hohenberg, P.; Kohn, W. Inhomogeneous Electron Gas. Phys. Rev. 1964, 136, B864-B71. 30 Kohn, W; Sham, L. J. Phys. Rev. A Self-Consistent Equations Including Exchange and Correlation Effects. 1965, 140, A1133-A38. 31 Zhao, Y.; Truhlar, D.G. The M06 Suite of Density Functionals for Main Group Thermochemistry, Thermochemical Kinetics, Noncovalent Interactions, Excited States, and Transition elements: Two New Functionals and Systematic Testing of Four M06-Class Functionals and 12 Other Functionals. Theor. Chem. Acc. 2008, 120, 215–241. 32 Gaussian 09, Revision B.01, Frisch, M. J.; Trucks, G. W.; Schlegel, H. B.; Scuseria, G. E.; Robb, M. A.; Cheeseman, J. R.; Scalmani, G.; Barone, V.; Mennucci, B.; Petersson, G. A. et al Gaussian, Inc., Wallingford CT, 2009. 33 Hratchian, H. P.; Schlegel, H. B. Accurate Reaction Paths Using a Hessian Based Predictor– Corrector Integrator. J. Chem. Phys. 2004, 120, 9918-9924.

21 ACS Paragon Plus Environment

Page 23 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

34

Hratchian, H. P.; Schlegel, H. B. Using Hessian Updating to Increase the Efficiency of a Hessian Based Predictor-Corrector Reaction Path Following Method. J. Chem. Theory and Comput. 2005, 1, 61-69. 35 Hanggi, P; Talkner, P; Borkovec, M. Reaction-Rate Theory: Fifty Years After Kramers. Rev. Mod. Phys. 1990, 62, 251-341. 36 Garrett, B. C.; Truhlar, D. G. Generalized Transition State Theory. Classical Mechanical Theory and Applications to Collinear Reactions of Hydrogen Molecules. J. Phys. Chem.1979, 83, 1052-1079. 37 Chase, M.W., Jr. NIST-JANAF Themochemical Tables, Fourth Edition. J. Phys. Chem. Ref. Data. 1998, 9, 1-1951. 38 Ervin, K. M.; Ho, J.; Lineberger, W.C. Ultraviolet Photoelectron Spectrum of NO2-. J. Phys. Chem. 1988, 92, 5405 39 Bartmess, J.E., "Negative Ion Energetics Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 40 Haber, K.S.; Zwanziger, F.X.; Campos, R.T.; Wiedmann, R.T.; Grant, E.R. Direct Determination of the Adiabatic Ionization Potential of NO2 by Multiresonant Optical Absorption. Chem. Phys. Lett. 1988, 144, 58. 41 Lias, S.G.; Liebman, J.F., "Ion Energetics Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 42 Katsumata, S.; Shiromaru, H.; Mitani, K.; Iwata, S.; Kimura, K. Photoelectron Angular Distribution and Assignments of Photoelectron Spectra of Nitrogen Dioxide, Nitromethane and Nitrobenzene. Chem. Phys. 1982, 69, 423. 43 Lias, S.G.; Bartmess, J.E.; Liebman, J.F.; Holmes, J.L.; Levin, R.D.; Mallard, W.G., "Ion Energetics Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 44 Knobel, Y.K.; Miroshnichenko, E.A.; Lebedev, Y.A. Heats of Combustion of Nitromethane and Dinitromethane: Enthalpies of Formation of Nitromethyl Radicals and Energies of Dissociation of Bonds in Nitro Derivatives of Methane. Bull. Acad. Sci. USSR, Div. Chem. Sci. 1971, 20, 425-428. 45 Afeefy, H.Y.; Liebman, J.F.; Stein, S.E., "Neutral Thermochemical Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 46 Khrapkovskii, G.M.; Tsyshevsky, R.V.; Chachkov, D.V.; Egorov, D.L.; Shamov, A.G. J. Mol. Struct.: THEOCHEM. 2010, 958, 1–6. 47 Adams, C.L.; Schneider, H.; Ervin, K.M.; Weber, J.M. Low-energy Photoelectron Imaging Spectroscopy of Nitromethane Anions: Electron Affinity, Vibrational Features, Anisotropies, and the Dipole-Bound State. J. Chem. Phys. 2009, 130, 074307. 48 Allam, S.H.; Migahed, M.D.; El Khodary, A. Electron Impact Study of Nitrobenzene and Nitromethane. Int. J. Mass Spectrom. Ion Phys. 1981, 39, 117. 22 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 27

49

Lias, S.G.; Levin, R.D.; Kafafi, S.A., "Ion Energetics Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 50 Rao, C.N.R. Lone-Pair Ionization Bands of Chromophores in the Photoelectron Spectra of Organic Molecules. Indian J. Chem. 1975, 13, 950. 51 Cox, J. D.; Pilcher, G. Thermochemistry of Organic and Organometallic Compounds; Academic Press: London and New York, 1970; pp. 643. 52 Dewar, M.J.S.; Shanshal, M.; Worley, S.D. Calculated and Observed Ionization Potentials of Nitroalkanes and of Nitrous and Nitric acids and Esters. Extension of the MINDO Method to Nitrogen-Oxygen Compounds. J. Am. Chem. Soc. 1969, 91, 3590. 53 Rosenstock, H.M.; Draxl, K.; Steiner, B.W.; Herron, J.T., "Ion Energetics Data" in NIST Chemistry WebBook, NIST Standard Reference Database Number 69, Eds. P.J. Linstrom and W.G. Mallard, National Institute of Standards and Technology, Gaithersburg MD, 20899, http://webbook.nist.gov, (retrieved January 29, 2015). 54 Kimura, K.; Katsumata, S.; Achiba, Y.; Yamazaki, T.; Iwata, S. Handbook of He Photoelectron Spectra of Fundamental Organic Molecules. Ionization Energies, Ab Initio Assignments, and Valence Electronic Structure for 200 Molecules; Japan Scientific Societies Press; Tokyo, New York, 1981; pp. 268. 55 Lebedeva, N.D.; Katin, Y.A.; Akhmedova, G.Y. Standard Enthalpy of Formation of Nitrobenzene. Russ. J. Phys. Chem. (Engl. Transl.) 1971, 45, 1192-1193. 56 Luo, Y.R. Handbook of Bond Dissociation Energies in Organic Compounds; CRC Press: Boca Raton, 2003; pp. 392. 57 Desfrancois, C.; Periquet, V.; Lyapustina, S.A.; Lippa, T.P.; Robinson, D.W.; Bowen, K.H.; Nonaka, H.; Compton, R.N. Electron Binding to Valence and Multipole States of Molecules: Nitrobenzene, Para- and Meta-Dinitrobenzenes. J. Chem. Phys. 1999, 111, 4569. 58 Katsumata, S.; Shiromaru, H.; Mitani, K.; Iwata, S.; Kimura, K. Photoelectron Angular Distribution and Assignments of Photoelectron Spectra of Nitrogen Sioxide, Nitromethane and Nitrobenzene. Chem. Phys. 1982, 69, 423. 59 Kobayashi, T., A Simple General Tendency in Photoelectron Angular Distributions of Some Monosubstituted Benzenes. Phys. Lett. 1978, 69, 105. 60 Lenchitz, C.; Velicky, R.W.; Silvestro, G.; Schlosberg, L.P. Thermodynamic Properties of Several Nitrotoluenes. J. Chem. Thermodyn. 1971, 3, 689-692. 61 Rice, B.M.; Sahu, S.; Owens, F.J. Density Functional Calculations of Bond Dissociation Energies for NO2 Scission in Some Nitroaromatic Molecules. J. Mol. Struct.: THEOCHEM. 2002, 583, 69–72. 62 Cohen, R.; Zeiri, Y.; Wurzberg, E.; Kosloff, R. Mechanism of Thermal Unimolecular Decomposition of TNT (2,4,6-Trinitrotoluene): A DFT Study. J. Phys. Chem. A 2007, 111, 11074–11083. 63 Potapov, V.K.; Kardash, I.E.; Sorokin, V.V.; Sokolov, S.A.; Evlasheva, T.I. Photoionization of Heteroaromatic Compounds. Khim. Vys. Energ. 1972, 6, 392. 64 Tsyshevsky, R. V.; Kuklja, M. M. Decomposition Mechanisms and Kinetics of Novel Energetic Molecules BNFF-1 and ANFF-1: Quantum-chemical Modeling. Molecules 2013, 18, 8500–8517.

23 ACS Paragon Plus Environment

Page 25 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

65

Pagoria, P.; Zhang, M.; Racoveanu, A.; DeHope, A.; Tsyshevsky, R.; Kuklja, M. M. 3-(4Amino-1,2,5-oxadiazol-3-yl)-4-(4-nitro-1,2,5-oxadiazol-3-yl)-1,2,5-oxadiazole, Molbank 2014, M824, 1-4. 66 Sharia, O.; Tsyshevsky, R.; Kuklja, M. M. Surface-Accelerated Decomposition of δ-HMX. J. Phys. Chem. Lett. 2013, 4, 730−734. 67 Kuklja, M. M.; Tsyshevsky, R. V.; Sharia, O. Effect of Polar Surfaces on Decomposition of Molecular Materials. J. Am. Chem. Soc. 2014, 136, 13289–13302. 68 Kimmel, A.V.; Sushko, P.V.; Shluger, A.L.; Kuklja, M.M. The Effect of Molecular and Lattice Structure on Hydrogen Transfer in Molecular Crystals Diamino-Dinitroethylene and Triamino-Trinitrobenzene, J. Phys. Chem. A 2008, 112, 4496-4450. 69 Kimmel, A.V.; Sushko, P.V.; Shluger, A.L.; Kuklja, M.M. Effect of Charged and Excited States on the Decomposition of 1,1-Diamino-2,2-Dinitroethylene Molecules. J. Chem. Phys. 2007, 126, 234711–234721. 70 Yuan, B.; Yu, Z.; Bernstein, E. R. Initial Decomposition Mechanism for the Energy Release from Electronically Excited Energetic Materials: FOX-7 (1,1-Diamino-2,2Dinitroethene,C2H4N4O4). J. Chem. Phys. 2014, 140, 074708-1–9. 71 Sikder, A.K.; Maddala, G.; Agrawal, J.P.; Singh, H. Important Aspects of Behaviour of Organic Energetic Compounds: A Review. J. Hazard. Mater. 2001, A84, 1–26. 72 Formation enthalpies were calculated in this work using atomization approach at M06/631+G(2df,p) level of theory. 73 Kinstle, T. H.; Stam, J. G. Vacuum Pyrolysis of Nitrostyrenes. J. Org. Chem. 1970, 35, 1771– 1774. 74 Khrapkovskii, G.M.; Shamov, A.G.; Nikolaeva, E.V.; Chachkov, D.V. Mechanisms of The Gas-Phase Decomposition of C-Nitro Compounds Inferred from Quantum Chemical Calculation. Russ. Chem. Rev. 2009, 78, 903-943. 75 Kiselev, V. G.; Gritsan, N. P. Unexpected Primary Reactions for Thermolysis of 1,1-Diamino2,2-Dinitroethylene (FOX-7) Revealed by Ab Initio Calculations. J. Phys. Chem. A 2014, 118, 8002–8008. 76 Korolev, V.L.; Petukhova, T.V.; Pivina, T.S.; Porollo, A.A.; Sheremetev, A.B.; Suponitskii, K. Yu.; Ivshin, V.P. Thermal Decomposition Mechanisms of Nitro-1,2,4-Triazoles: A Theoretical Study. Russ. Chem. Bull. 2006, 55, 1388–1410. 77 Yu, Z.; Bernstein, E. R. On the Decomposition Mechanisms of New Imidazole-Based Energetic Materials. J. Phys. Chem. A 2013, 117, 1756–1764. 78 Brill, T.B.; James, K.J. Kinetics and Mechanisms of Thermal Decomposition of Nitroaromatic Explosives. Chem. Rev. 1993, 93, 2667–2692. 79 Brill, T.B.; James, K.J. Thermal Decomposition of Energetic Materials. 61. Perfidy in the Amino-2,4,6-Trinitrobenzene Series of Explosives. J. Phys. Chem. 1993, 97, 8752-8758. 80 Shao, J.; Cheng, X.; Yang, X. The C-NO2 Bond Dissociation Energies of Some Nitroaromatic Compounds: DFT Study. Struct. Chem. 2006, 17, 547–550. 81 Politzer, P.; Concha, M.C.; Grice, M.E.; Murray, J.S.; Lane, P.; Habibollazadeh, D. Computational investigation of the structures and relative stabilities of amino/nitro derivatives of ethylene. J. Mol. Struct.: THEOCHEM 1998, 452, 75–83. 82 Gindulyté, A.; Masaa, L.; Huang, L.; Karle, J. Proposed mechanism of 1,1-DiaminoDinitroethylene Decomposition: A Density Functional Theory Study. J. Phys. Chem. A 1999, 103, 11045–11051. 24 ACS Paragon Plus Environment

The Journal of Physical Chemistry

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 27

83

Gindulyté, A.; Masaa, L.; Huang, L.; Karle, J. Ab Initio Study of Unimolecular Decomposition of Nitroethylene. J. Phys. Chem. A 1999, 103, 11040–11044. 84 We were unable to obtain the reaction energy for the NO loss path V8 (Fig. 4), as discussed in supplementary materials. The computational issue is of the technical nature and hence it does not qualitatively affect conclusions regarding thermal stability of DNBTT. 85 Hall, R. B.; Williams, F. Photodecomposition and Electronic Structure of Lead Azide. J. Chem. Phys. 1973, 58, 1036–1042. 86 Aluker, E. D.; Krechetov, A. G.; Mitrofanov, A. Y.; Zverev, A. S.; Kuklja, M. M. Topography of Photo-Initiation in Molecular Materials, Molecules 2013, 18, 14148-14160. 87 Aluker, E. D.; Krechetov, A. G.; Mitrofanov, A. Y.; Zverev, A. S.; Kuklja, M. M. Applicability and Limitations of the Thermal Mechanism of Laser Initiation of Energetic Materials, J. Phys. Chem. C 2012, 116, 24482−24486. 88 Im, H.-S.; Bernstein, E. R. On the Initial Steps in the Decomposition of Energetic Materials From Excited Electronic States. J.Chem. Phys. 2000, 113, 7911−7918. 89 Bhattacharya, A.; Guo, Y. Q.; Bernstein, E. R. Nonadiabatic Reaction of Energetic Molecules. Acc. Chem. Res. 2010, 43, 1476–1485. 90 Bernstein, E. R. On the Release of Stored Energy from Energetic Materials. Adv. Quantum Chem. 2014, 69, 31-71. 91 Kuklja, M. M. On the Initiation of Chemical Reactions by Electronic Excitations in Molecular Solids. Appl. Phys. A 2003, 76, 359–366 92 Tsyshevsky, R. V., Rashkeev, S. N., & Kuklja, M. M. (2015). Defect States at Organic– inorganic Interfaces: Insight from First Principles Calculations for Pentaerythritol Tetranitrate on MgO Surface. Surface Science, 2015, 637, 19-28.

25 ACS Paragon Plus Environment

Page 27 of 27

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The Journal of Physical Chemistry

TOC image

1 ACS Paragon Plus Environment