First-Principles Study of the Four Polymorphs of Crystalline Octahydro

Oct 12, 2007 - Institute for Computation in Molecular and Materials Science and Department of Chemistry, Nanjing University of Science and Technology,...
1 downloads 12 Views 271KB Size
J. Phys. Chem. B 2007, 111, 12715-12722

12715

First-Principles Study of the Four Polymorphs of Crystalline Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine Weihua Zhu,*,† Jijuan Xiao,† Guangfu Ji,‡ Feng Zhao,‡ and Heming Xiao*,† Institute for Computation in Molecular and Materials Science and Department of Chemistry, Nanjing UniVersity of Science and Technology, Nanjing 210094, China, and Laboratory for Shock WaVe and Detonation Physics, Institute of Fluid Physics, China Academy of Engineering Physics, Mianyang 621900, China ReceiVed: June 28, 2007; In Final Form: August 15, 2007

The electronic structure and vibrational properties of the four polymorphs of crystalline octahydro-1,3,5,7tetranitro-1,3,5,7-tetrazocine (HMX) have been studied using density functional theory within the local density approximation. The results show that the states of N in the ring make more important contributions to the valence bands than these of C and N of NO2 and so N in the ring acts as an active center. From the low frequency to high-frequency region, the molecular motions of the vibrational frequencies for the four HMX polymorphs present unique features. It is also noted that there is a relationship between the band gap and impact sensitivity for the four HMX polymorphs. From the cell bond order per unit volume, we may infer the variation order of crystal bonding for the four polymorphs and so predict their impact sensitivity order as follows: β-HMX < γ-HMX < R-HMX < δ-HMX, which is in agreement with their experimental order.

1. Introduction Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine, known as HMX, is an important and commonly used energetic ingredient in various high performance explosives and propellant formulations due to its thermal stability and high detonation velocity relative to other explosives.1,2 HMX is known to exist in four crystalline phases, denoted as R, β, δ, and γ.3,4 Although β-HMX is the thermodynamically stable form under ambient conditions, post-mortem analysis5 of samples recovered from safety experiments involving low-velocity projectile impacts on the HMX-based plastic-bonded explosive (PBX)-95016 has shown the formation of δ-HMX in the vicinity of damaged regions within the material. It is also known that a layer of δ-HMX is formed at the solid-melt interphase in deflagrating HMX. These observations present a safety concern since δ-HMX is considerably more sensitive than β-HMX. Therefore, a detailed and comparative study of the four HMX polymorphs offers the possibility of determining microscopic properties and understanding their explosive properties. Materials such as propellants and explosives contain tightly bonded groups of atoms that retain their molecular character until a sufficient stimulus is applied to cause exothermic dissociation. This in turn triggers further dissociation leading to initiation or ignition. The macroscopic behavior is ultimately controlled by microscopic properties such as the electronic structure and interatomic forces. Thus, a desire to probe more fundamental questions relating to the basic properties of HMX as a solid energetic material is generating significant interest in the basic solid-state properties of such energetic systems. Although the detailed decomposition mechanisms by which energetic materials release energy under mechanical shock are * To whom correspondence should be addressed. Fax: +86-2584303919. E-mail: [email protected]; [email protected]. † Nanjing University of Science and Technology. ‡ China Academy of Engineering Physics.

still not well understood, it has been suggested that these decompositions may result from transferring thermal and mechanical energy into the internal degrees of freedom of the molecules in energetic solids.7-9 It is thus very important to understand the electronic and vibrational properties of the four HMX polymorphs. Many experimental studies have been performed to investigate the vibrational spectra of the HMX polymorphs.10-14 The investigation of the microscopic properties of energetic materials, which possess a complex chemical behavior, remains to be a challenging task. Theoretical calculations can play an important role in investigating the physical and chemical properties of complex solids at the atomic level and the establishment of the relationships between their structure and function. Previously, density functional theory (DFT) calculations based on norm-conserving pseudopotentials and a local density approximation (LDA) were performed to examine the energetics of the three pure polymorphic HMX (R, β, δ).15 Afterward, Lewis16 used the same method to investigate the energetics of HONO formation in the three pure polymorphs of condensed-phase HMX, where intermolecular hydrogen transfer occurs. Recently, Ye and Koshi17 have evaluated the energy transfer rates of the three pure HMX polymorphs in terms of the density of vibrational states and the unharmonic vibronphonon coupling term, which were calculated by using a flexible potential containing both intra- and intermolecular terms. As the electronic structure and vibrational properties of the four HMX polymorphs are not systematically investigated and compared, there is a clear need to gain an understanding of those at the ab initio level. In this study, we report a systematic study of the electronic structure and vibrational properties of the four HMX crystalline phases (R, β, δ, and γ) from DFT. Our main purpose here is to examine the differences in the microscopic properties of the polymorphs and to understand their structure-function relationships.

10.1021/jp075056v CCC: $37.00 © 2007 American Chemical Society Published on Web 10/12/2007

12716 J. Phys. Chem. B, Vol. 111, No. 44, 2007

Zhu et al.

Figure 1. Unit cell of R-HMX (a), β-HMX (b), δ-HMX (c), and γ-HMX (d). Gray, blue, red, and white spheres stand for C, N, O, and H atoms, respectively.

The remainder of this paper is organized as follows. A brief description of our computational method is given in section 2. The results and discussion are presented in section 3, followed by a summary of our conclusions in section 4. 2. Computational Method The calculations performed in this study were done within the framework of DFT,18 using Vanderbilt-type ultrasoft pseudopotentials19 and a plane-wave expansion of the wave functions. The self-consistent ground state of the system was determined by using a band-by-band conjugate gradient technique to minimize the total energy of the system with respect to the plane-wave coefficients. The electronic wave functions were obtained by a density-mixing scheme20 and the structures were relaxed by using the Broyden, Fletcher, Goldfarb, and Shannon (BFGS) method.21 The LDA functional proposed by Ceperley and Alder22 and parametrized by Perdew and Zunder23 named CA-PZ, was employed. The cutoff energy of plane waves was set to 300.0 eV. Brillouin zone sampling was performed by using the Monkhost-Pack scheme with a k-point grid of 3 × 3 × 3, 4 × 2 × 3, 3 × 3 × 1, and 2 × 3 × 2 for R-, β-, δ-, and γ-HMX, respectively. The values of the kinetic energy cutoff and the k-point grid were determined to ensure the convergence of total energies. The R phase of HMX crystallizes in an orthorhombic Fdd2 space group and contains eight H8C4N8O8 molecules per unit cell.24 The β phase contains two H8C4N8O8 molecules per unit cell in a monoclinic lattice with space group P21/c.25 The δ phase crystallizes in a hexagonal P61 space group with six H8C4N8O8 molecules per unit cell.26 The γ phase contains two 2H8C4N8O8·0.5H2O molecules per unit cell in a monoclinic lattice with space group Pn.4 The γ phase is a hydrate and not a true polymorph. When the explosive decomposition of this HMX hydrate occurs water does not work. Strictly speaking only three HMX phases should be named polymorphs. Figure 1 displays the unit cell of the four HMX crystalline phases, and conformation and atomic numbering of H8C4N8O8 molecule in the β-HMX phase is shown in Figure 2. Starting from the above-mentioned experimental structures, the geometry relaxation was performed to allow the ionic configurations, cell shape, and volume to change. In the geometry relaxation, the total energy of the system was converged less than 2.0 × 10-5 eV, the residual force less than 0.05 eV/Å, the displacement of atoms less than 0.002 Å, and the residual bulk stress less than 0.1 GPa. The Mulliken charges and bond populations were investigated using a projection of the plane wave states onto a linear combination of atomic

Figure 2. Conformation and atomic numbering of H8C4N8O8 molecule in β-HMX.

TABLE 1: Experimental and Relaxed Lattice Constants (Å) for β-HMX this work a b c β

LDA

PBE

PW91

experimental25

6.539 11.030 8.689 123.9°

6.625 11.202 8.824 124.3°

6.617 11.188 8.810 124.2°

6.540 11.050 8.700 124.3°

orbitals (LCAO) basis set,27,28 which is widely used to perform charge transfers and populations analysis. The phonon frequencies at the gamma point have been calculated from the response to small atomic displacements.29 3. Results and Discussion 3.1. Bulk Properties. As a base for studying other HMX phases and as a well-studied benchmark, we apply three different functionals to bulk β-HMX as a test. The LDA (CA-PZ) and generalized gradient approximation (GGA) (PBE30 and PW9131) functionals were selected to fully relax the β phase without any constraint. The calculated lattice parameters are given in Table 1 together with their experimental values.25 It is found that the errors in the LDA (CA-PZ) results are slightly smaller than that in the GGA (PBE and PW91) results in comparison with the experimental values. This shows that the accuracy of LDA is better than that of the GGA functionals. Table 2 presents the bond lengths and bond angles of the β phase along with the corresponding experimental data. It is seen that the bond lengths and bond angles compare well with experimental values. The comparisons confirm that our computational parameters are reasonably satisfactory. We thus used LDA in all subsequent calculations.

Four Polymorphs of Crystalline HMX

J. Phys. Chem. B, Vol. 111, No. 44, 2007 12717

TABLE 2: Experimental and Relaxed Bond Lengths (Å) and Bond Angles (deg) for β-HMX bond length N2-C1 N2-C2′ N3-C1 N3-C2 N2-N1 N3-N4 N1-O1 N1-O2 N4-O3 N4-O4 C1-H1 C1-H2 C2-H3 C2-H4

bond angle

LDA

expt25

1.435 1.461 1.440 1.421 1.358 1.374 1.253 1.245 1.246 1.243 1.106 1.108 1.107 1.106

1.448 1.471 1.455 1.437 1.354 1.373 1.233 1.222 1.210 1.204 1.111 1.091 1.101 1.095

O1-N1-O2 N2-N1-O1 N2-N1-O2 N1-N2-C2′ N1-N2-C1 C1-N2-C2′ N2-C1-N3 N2-C1-H1 N2-C1-H2 N3-C1-H1 N3-C1-H2 H1-C1-H2 O3-N4-O4 N3-N4-O3 N3-N4-O4 N4-N3-C1 N4-N3-C2 C1-N3-C2 N2′-C2-N3 N3-C2-H4 N3-C2-H3 N2′-C2-H4 N2′-C2-H3 H3-C2-H4

LDA

expt25

125.7 118.5 115.8 115.8 118.7 123.7 111.4 106.8 111.8 107.3 111.6 108.0 126.5 116.5 117.0 117.9 118.2 123.5 110.6 108.2 109.9 109.9 109.0 109.1

125.9 118.0 116.1 115.2 118.2 122.4 113.5 108.5 109.8 106.6 111.8 106.2 126.7 115.9 117.4 117.4 118.2 123.8 110.2 107.3 110.5 110.1 109.2 109.6

3.2. Electronic Structure. The calculated total density of states (DOS) and partial DOS (PDOS) for the four HMX phases are displayed in Figures 3-6, respectively. The DOS of the four HMX phases are finite at the Fermi energy level. This is because the DOS contain some form of broadening effect. In the upper valence band, all of the four HMX phases have a sharp peak near the Fermi level, which shows that the top valence bands of their band structures are flat. The top of the DOS valence band shows three main peaks for R- and δ-HMX and two main peaks for β- and γ-HMX. These peaks are predominately from the p states. After that, several main peaks in the upper valence band are superimposed by the s and p states. The conduction band of each phase is dominated by the p states. This indicates that the p states for each HMX phase play a very important role in its chemical reaction. The atom-resolved DOS and PDOS of the four HMX phases are also shown in Figures 3-6, respectively. The main features can be summarized as follows. (i) In the upper valence band, the PDOS of the states of N in the ring are far larger than that of the states of C and N of NO2. It is expected that the states of N in the ring make more important contributions to the valence bands than these of C and N of NO2. This shows that N in the ring acts as an active center. (ii) Some strong peaks occur at the same energy in the PDOS of a particular C atom and a particular N atom in the ring. It can be inferred that the two atoms are strongly bonded. Similarly, a particular N atom of NO2 and a particular N atom in the ring are strongly bonded. (iii) There are some differences in the PDOS of the C atoms for the four HMX phases. This is due to the differences in their local molecular packing. The same is true of the N atoms in the ring and the N atoms of NO2. (iv) In the conduction band region of DOS, the peaks are dominated by the N-p states of ring and the N-p states of NO2. (v) In the energy range from -5.0 to 0 eV, the DOS of R-, β-, and γ-HMX are superimposed by the states of C and N in the ring, whereas that of δ-HMX mainly arises from the states of N in the ring. It is inferred that the C-N bond fission in the ring may be favorable in the decomposition of crystalline R-, β-, and γ-HMX. This is consistent with the previous conclusion32 that in a bulk

Figure 3. Total and partial density of states (DOS) of C-states, N in ring-states, N of NO2-states, and R-HMX. The Fermi energy is shown as a dashed vertical line.

Figure 4. Total and partial density of states (DOS) of C-states, N in ring-states, N of NO2-states, and β-HMX.

Figure 5. Total and partial density of states (DOS) of C-states, N in ring-states, N of NO2-states, and δ-HMX.

condensed phase the C-N bond scission reaction of the ring for R-HMX may be energetically favorable because of the steric constrains that will disfavor N-NO2 bond dissociation. Of

12718 J. Phys. Chem. B, Vol. 111, No. 44, 2007

Figure 6. Total and partial density of states (DOS) of C-states, N in ring-states, N of NO2-states, and γ-HMX.

course, near areas with large surfaces or voids there is no confining environmental and N-N bond dissociation would likely be more favorable. This is in agreement with the previous reports33,34 that the N-N bond scission reaction of hexahydro1,3,5-trinitro-1,3,5-triazine (RDX) in the gas-phase would be favorable. For the N-NO2 and C-N bond fission reactions of δ-HMX, it is difficult to judge which one is more favorable since the contribution of the states of C and of NO2 to the total DOS from -5.0 to 0 eV is very small. Recently, Feng and Li35 have investigated two polymorphs of flufenamic acid and elucidated the difference in the solid-state reaction of the two polymorphs based on their electronic structures. 3.3. Vibrational Properties. Here we investigate the vibrational properties of the four HMX phases. β-HMX contains a ring conformation that the NO2 groups adopt a chairlike arrangement. This gives the entire molecule a center of symmetry.24,25 The ring conformation of R-, δ-, and γ-HMX is a “boat” form that the NO2 groups are positioned on one side of the molecule.24,26 It is interesting to note that the R-, δ-, and γ-HMX phases have similar vibrational spectra, whereas the β phase has a significantly different vibrational spectrum.10-14 Therefore, the calculated vibrational frequencies for β-HMX are presented in Table 3, whereas these for R-, δ-, and γ-HMX are presented in Table 4. The corresponding experimental values are also listed for comparison. The lattice mode spectra of β-, δ-, and γ-HMX displays unique features in the region 27-80 cm-1, whereas there is only one librational vibration mode at 31.2 cm-1 in the R-HMX spectrum. β-HMX has five frequencies in this region, which are in agreement with the experimental reports.13 The lattice mode spectra of the other three phases were not reported in this region by the experiments because their lattice region spectra could not be routinely obtained. In the region of 80-500 cm-1, the molecular motions of the first four frequencies are the twisting of NO2 about N-N bonds, denoted γNN(NO2), and the motions corresponding to the following frequencies is concentrated in the distortions of C4N4 rings. In the region of 500-730 cm-1, the molecular motions of the frequencies are the bending of N-N-O angles, denoted b(NNO), and the stretching of N-N bonds, denoted V(NN). The 637.8 cm-1 mode corresponding to the bending of N-N-O and N-N-C angles and the stretching of C-N bonds appears in γ-HMX but not in any of the other three polymorphs. The modes in the range of 730-760 cm-1 are the wagging of N

Zhu et al. atoms out of the NO2 plane, denoted σ(NO2). We note that R, δ, and γ-HMX share similarities in band separation, whereas β has a different spectral pattern. This is because the crystal structures of R, δ, and γ-HMX differ from one another primarily through small differences in the relative positioning of the NO2 groups, while β-HMX has a different ring conformation from the other phases. In the region of 760-1260 cm-1, the molecular motions of the frequencies are mainly the stretching of C4N4 rings including the symmetric stretch of C-N-C and N-NC2 bonds, the asymmetric stretch of C-N-C, N-C-N, and N-N-C2 bonds, and the rocking in the CH2 plane. The modes of 1220.4 cm-1 in R, 1234.2 cm-1 in δ, and 1228.2 cm-1 in γ correspond to the stretching of only one N-O bond in NO2 moieties, denoted V(NO), plus the stretching of N-N bonds and the bending of C-N-C angles, which do not appear in β-HMX. The R-HMX phase has the smaller motion modes of the frequencies in this region than any of the other three polymorphs. In the region of 1260-1350 cm-1, the molecular motions of the frequencies involve the symmetric stretching of O-N-O bonds, denoted Vs(NO2), and the twisting about bisector of H-C-H angles, denoted γ(CH2). The modes of the four frequencies in the range of 1350-1500 cm-1 involve the wagging of H atoms out of the CH2 plane, denoted ω(CH2), and the bending of H-C-H angles, denoted b(HCH), which does not appear in R-HMX but appears in the other three polymorphs. In the region of 1500-1750 cm-1, the molecular motions of the frequencies are the asymmetric stretching of O-N-O bonds, denoted Vas(NO2). The calculated frequencies for the β-HMX phase in this region are much larger than the corresponding experimental values. The motions of these frequencies is concentrated in the NO2 moieties; thus, this discrepancies may be due to intermolecular hydrogen bonding interactions present in the crystal lattice, which are not well described by DFT. In the region of 2900-3100 cm-1, the molecular motions of the frequencies are divided into two separate groups. The motions of the first four frequencies involve the symmetric stretching of H-C-H bonds, denoted Vs(CH2), whereas these of the next four frequencies are concentrated in the asymmetric stretching of H-C-H bonds, denoted Vas(CH2). These modes do not appear in R-HMX but appear in the other three polymorphs. Recently, Stevens et al.36 has reported the single-crystal polarized Raman spectra for β-HMX. It can be seen from Table 3 that our calculated frequencies are in agreement with their experimental ones. They have also presented complete symmetry assignments for the Raman-active modes of β-HMX. This is useful for understanding the proposed vibrational mechanisms for the initiation of detonation. In a word, the β-HMX phase has clearly very different vibrational properties from R-, δ-, and γ-HMX. This is supported by previous experimental studies.10,12,14 From the low frequency to high-frequency region, the molecular motions of the frequencies for the four HMX polymorphs present unique features, which could be used to distinguish the polymorphs easily from one another. At lower (below 120 cm-1) or higher (above 1150 cm-1) frequencies, the motions seem very localized in certain atom groups, whereas elsewhere the motions are diffusely distributed among different atom groups. Since the torsional motions of the molecule’s functional groups are usually supposed to be highly coupled to the other moiety of the molecule,37 it is possible that the torsion motions of the NO2 groups act as doorways through which kinetic energy can flow into the molecule from its surroundings. These suggestions show

Four Polymorphs of Crystalline HMX

J. Phys. Chem. B, Vol. 111, No. 44, 2007 12719

TABLE 3: Vibrational Frequencies (cm-1) for β-HMX experimental infrared (IR) mode

assignmenta

1

librational vibration translational vibration rotational lattice vibration librational vibration translational vibration γNN(NO2) γNN(NO2) γNN(NO2) γNN(NO2) σ(CNC) b(NNC) σ(CNC) b(NNC), b(CNC) b(NNC), b(NCN) b(NNC), b(NNO) b(NNC) ν(CN), ν(NN), b(CNC) ν(NN), b(CNC), b(NNC) b(NNO), b(NNC) b(NNO), b(NNC) b(CNC), b(NNC) b(CNC), b(NNO), b(NNC) b(NNO) b(NNO) b(NNO), ν(NN) b(NNO), ν(NN) b(NNO), ν(NN) b(NNO), ν(NN) b(ONO), ν(CN) σ(ONO) σ(ONO) σ(ONO) σ(ONO) b(NCN) νs(NC2) νs(NC2) νs(NNC2) νs(NNC2) νas(NNC2)

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

experimental

Raman (R)

this work ref 14 ref 10 ref 13 ref 12 ref 13

950

953

41

νas(CNN), F(CH2)

930.9

965

966

42

νas(CNN), F(CH2)

938.4

65

43

79

81

1090

1088

90 110 126 140 152

97

1168 1190

1170 1192

1248 1268

1251 1270

1312

1310 1318

1350 1368

1351 1369

1418 1438

1420 1438

1460

1461

1532

1535

1558

1560

2992 3028

2994 3028

3037

3038

63.2 65.9

224

369.4 398.8 404.0 424.8 426.1 593.3 595.5 629.0 630.4 641.0 658.6 722.4 750.2 750.9 751.0 752.7 774.4 837.9 838.2 878.8 879.6 920.7

419

385 420

438

441

600

605

627

754

230

231

281 312

282 318

752

755

758 769

762 773

832 872

827 871

833 872

945

967

965

966

ν(NN), F(CH2)

1071.2 1088

1090

1089

44

νas(NNC2)

1086.1

45 46 47 48 49 50 51 52 53 54 55 56

νas(CNN), F(CH2) νas(CNN), F(CH2) νas(NC2) νas(NC2) νas(NC2) νas(NC2) νs(NO2) νs(NO2) νs(NO2) γ(CH2) νs(NO2) γ(CH2)

1146.1 1150.0 1196.9 1213.9 1228.0 1228.4 1275.3 1280.6 1293.4 1309.0 1329.2 1332.0

1325

1320

1325

1348

1348

1146 1204 1239 1279 1296

358

364

57

γ(CH2)

1344.9 1349

412

417

432

437

58 59 60 61

γ(CH2) ω(CH2) ω(CH2) ω(CH2)

1346.4 1367.8 1390.8 1385 1394.5 1395

597

600

638

640

62 63 64 65 66 67 68 69 70 71 72 73 74 75 76 77 78

ω(CH2) b(HCH) b(HCH) b(HCH) b(HCH) νas(NO2) νas(NO2) νas(NO2) νas(NO2) νs(CH2) νs(CH2) νs(CH2) νs(CH2) νas(CH2) νas(CH2) νas(CH2) νas(CH2)

1423.8 1426.8 1442.8 1446.5 1459.5 1700.4 1704.0 1706.9 1715.2 2955.4 2955.5 2955.6 2955.7 3029.9 3030.3 3030.8 3031.1

658

760 772

947

149 155

630 660

ref 14 ref 10 ref 13 ref 12 ref 13

923.2

59

84.0 96.1 111.9 112.9 147.6 156.2 217.0 217.6 256.1 295.6 332.0 337.5

assignmenta

Raman (R)

ν(NN), F(CH2)

36

78.8

mode

this work

40

38.8 57.3

infrared (IR)

662 721 759

663 722 762 763

834

834

881

884

948

1146 1205 1240 1280

1397

1145 1203

1281 1298

1395 1401

1433

1433

1434

1462

1465

1463

1534

1540

1538

1563

1570

1565

2985 2992 3027 3037

2978 2983 3027 3037

a

γNN(XY2): twist of XY2 about NN bond; σ(XY2): wag of X atom out of XY2 plane; b(XYZ): bend of X-Y-Z angle; νs(XY2): symmetric stretch of Y-X-Y bonds; νs(XXY2): symmetric stretch of X-X-Y2 bonds; νas(XY2): asymmetric stretch of Y-X-Y bonds; νas(XXY2): asymmetric stretch of X-X-Y2 bonds; F(XY2): rocking in XY2 plane; γ(XY2): twist about bisector of Y-X-Y angle; ω(XY2): wag of Y atoms out of XY2 plane.

that the NO2 groups play a very important role in the decomposition and detonation of HMX. 3.4. Electronic Structure and Impact Sensitivity. In this section, an attempt is made to correlate the impact sensitivity of the four HMX phases with their electronic structure. Band gap is an important parameter to characterize the electronic structure of solids. Table 5 presents the energy gaps between valence and conduction bands for R-, β-, δ-, and γ-HMX. This shows that the band gap decreases in the sequence of β-, γ-, R-, δ-HMX, whereas their experimental impact energy (shown in Table 5) increases in the following order: β-HMX < γ-HMX ∼ R-HMX < δ-HMX.38 Therefore, there is a relationship between the band gap and impact sensitivity for the four HMX polymorphs. This indicates that the smaller the band gap is,

the easier the electron transfers from the valence band to the conduction band, and the more the HMX phase becomes decomposed and exploded. In the previous studies,39,40 a “principle of the easiest transition” (PET) was put forward to investigate the relationship between the band gap and impact sensitivity for the metal azides. Although these calculations are performed at the semiempirical discrete variational XR (DVXR) and extended Hu¨ckel-crystal orbital (EH-CO) levels, the results have shown that their band gap could be correlated with their explosive characters. Our recent study on the heavy-metal azides within the framework of periodic DFT41 has also confirmed the relationship between the band gap and impact sensitivity. In several papers,42-45 Gilman has emphasized that the role of HOMO-LUMO (highest occupied molecular orbital-

12720 J. Phys. Chem. B, Vol. 111, No. 44, 2007

Zhu et al.

TABLE 4: Vibrational Frequencies (cm-1) for r-HMX, δ-HMX, and γ-HMX R-HMX

δ-HMX

experimental mode

assignmenta

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 61 62 63 64 65 66 67 68 69 70 71

librational vibration translational vibration rotational lattice vibration librational vibration translational vibration γNN(NO2) γNN(NO2) γNN(NO2) γNN(NO2) σ(CNC) b(NNC) σ(CNC) b(CNC) b(NNC), b(CNC) b(NNC) ν(CN), ν(NN), b(NCN) ν(NN), ν(CN), b(CNC) ν(NN), b(NNO) b(NNC), b(NNO), F(CH2) b(NNO), b(NNC), F(CH2) b(NNO), b(NNC), F(CH2) b(CNC), F(CH2) b(NNO) b(NNO) b(NNO), ν(NN) ν(NN), b(ONO) ν(NN), b(ONO) b(NNO), b(NNC), ν(CN) ν(NN), b(ONO) σ(ONO) σ(ONO) σ(CNC) σ(ONO) σ(ONO) νs(NC2) νs(NC2) b(ONO), νs(NNC2) b(ONO), νs(NNC2) b(ONO), νs(NNC2) νs(NNC2) ν(NN), b(ONO) ν(CN), F(CH2) νs(NNC2) νs(NNC2), F(CH2) νs(NNC2), F(CH2) νas(NC2) νas(NC2) ν(NO), ν(NN), b(CNC) νas(NC2) νas(NC2) νs(NO2) νs(NO2) νs(NO2) γ(CH2) νs(NO2) γ(CH2) ω(CH2) γ(CH2) γ(CH2) ω(CH2) ω(CH2) b(HCH) ω(CH2) b(HCH) b(HCH) b(HCH) νas(NO2) νas(NO2) νas(NO2) νas(NO2) νs(CH2)

this work

IR14

IR10

R12

31.2

214.0 248.2 332.5 344.1

451

450

603 622

603 625

594 620

646

648

646

691.5 732.4 736.9 741.9 747.7 759.4

714 741 751

714 742 753

712 740 750 758

766 770

768

856.9 867.2 868.8

847 864 878 915

848 862

944

846 878

920 945

928 945

1032 1090 1110

1030 1085

1220

1215

936.3 1075.7 1095.8 1155.8 1208.0 1220.4 1247.9 1250.1 1268.0 1273.0 1310.6 1314.1 1330.5 1335.6

1030 1089 1109 1148 1215 1259 1268 1280 1319 1364 1370 1386 1393 1414 1432 1451 1441

1559.6 1570.3 1574.4

this work 32.6 46.5 55.6 56.2 70.9 92.5 95.4 105.8 105.8 124.9 186.9 187.8 216.6 241.8 339.2 359.7 374.3 395.8 395.9 396.0

90.1 92.0 104.9 108.8 128.3 175.1

406.8 409.9 450.2 454.9 587.0 591.1 611.5 619.1 649.5

1541 1561

1258 1280

1280

1320

1318

1370

1368 1385 1391

1394 1416 1432

1412 1422 1448

IR14

IR10

400

696.6 727.8 733.2 742.5 742.6 751.4 824.1 835.8 866.6 869.1 888.3 906.3 911.6 912.1 989.4 1036.6 1089.5 1154.1 1206.5 1234.2 1254.9 1255.3 1281.5 1293.0 1304.3 1305.2 1344.2 1360.5 1385.4 1390.8 1400.3 1400.3 1400.9 1400.9

599 621

605 625

646

655

654

711 733 750

713 740 750

713 735 751 763

763 767 830 841 867

765 848 866

846 870

913 930 941

910 925 940

930 943

1016 1087 1109 1147 1204 1222

1040 1088 1110

1092 1111

1249 1257 1274 1290

1225

1270

1256 1275 1291

1325 1340

1317 1332

1370 1383 1395

1370

1371 1382 1393

1419

1420

1545 1562 1611.8

1215

1320 1341

1451 1539 1561

experimental R12

392 446 473 590 601 622

470 585.8 589.9 610.4 637.6 639.1

1468 1550

γ-HMX

experimental

1393

1415 1422 1422

this work 27.2 43.7 51.9 53.3 68.6 89.2 93.3 101.4 105.8 129.2 191.6 191.7 216.6 255.4 329.7 347.1 386.9 392.7 404.2 408.9 450.2 457.9 579.7 592.8 610.9 632.2 632.3 637.8 712.6 730.3 730.8 744.7 745.6 759.0 843.8 844.1 863.1 864.6 885.4 905.6 910.8 920.0 1013.9 1081.2 1097.1 1153.3 1208.5 1228.2 1249.4 1249.4 1274.9 1275.9 1296.7 1321.5 1330.2 1336.7 1359.8 1377.2 1382.1 1383.5 1403.0 1420.8 1423.4 1442.5 1442.9

IR10

R12

205 224 328 365 402 458

658

592 618 636 648

712

710

622

735 751 768 845

753 840 847 878

918 941 1016 1090 1112 1149 1221 1255

928 940 995 1090 1112 1166 1190 1225 1259 1271

1280 1320

1375 1398

1294 1319 1330

1422

1375 1385 1392 1411 1419

1454

1440 1452

1545.3

1550

1557

1586.2 2923.7

1570

1563 1573 2918

1455 1560

1536 1562

Four Polymorphs of Crystalline HMX

J. Phys. Chem. B, Vol. 111, No. 44, 2007 12721

TABLE 4: Continued R-HMX

δ-HMX

experimental mode 72 73 74 75 76 77 78

assignmenta

this work

νs(CH2) νs(CH2) νs(CH2) νas(CH2) νas(CH2) νas(CH2) νas(CH2)

IR14

IR10

2970

γ-HMX

experimental R12

2975 3047

3049 3056

this work

2963.3

IR14

IR10

2970

experimental R12

2974

3058 3058 3053

3058

this work 2923.7 2935.5 2935.5 3061.1 3061.2 3070.8 3071.0

IR10

R12 2928 2982 3032 3045

a

γNN(XY2): twist of XY2 about NN bond; σ(XY2): wag of X atom out of XY2 plane; b(XYZ): bend of X-Y-Z angle; νs(XY2): symmetric stretch of Y-X-Y bonds; νs(XXY2): symmetric stretch of X-X-Y2 bonds; νas(XY2): asymmetric stretch of Y-X-Y bonds; νas(XXY2): asymmetric stretch of X-X-Y2 bonds; F(XY2): rocking in XY2 plane; γ(XY2): twist about bisector of Y-X-Y angle; ω(XY2): wag of Y atoms out of XY2 plane.

TABLE 5: Experimental Impact Energy and Calculated Band Gap and Cell Bond Order per Unit Volume for r-HMX, β-HMX, δ-HMX, and γ-HMX impact energy (Kg/cm2)38 band gap (eV) cell bond order per unit volume

R-HMX

β-HMX

δ-HMX

γ-HMX

0.20 2.42 0.073

0.75 3.62 0.080

0.10 0.021 0.057

0. 20 3.38 0.078

lowest unoccupied molecular orbital) gap closure in molecules suffering shear strain. Further reports34,46 on the excitonic mechanism of detonation initiation show that the pressure inside the impact wave front reduces the band gap between valence and conducting bands and promotes the HOMO-LUMO transition within a molecule. Although these studies have suggested that the HOMO-LUMO gap in molecules suffering shear strain, impact wave, or distortion relates directly to the sensitivity, they further support our conclusion here that there is the relationship between the band gap and sensitivity. It has been suggested that initiation of combustion or detonation in energetic solids is associated with small regions of the solid termed “hot spots”.47-49 In general, the initiation of explosion in HMX is concerned with the decomposition of a small volume of material in the region of the “hot spot” at a high temperature and in a short time. The thermal decomposition of the energetic material is prior to its detonation; thus, the information of the thermal decomposition seems to provide an important basis for understanding the explosive characteristics. For the crystal as a whole, a more meaningful quantity is the total cell bond order that is the sum of all bond order values in the crystal.50 Since the four HMX polymorphs contain different number of atoms and bonds within the unit cell of each phase, a much more useful concept is the cell bond order per unit volume. These numbers are listed in Table 5. It turns out that β-HMX has the highest value of 0.080 for the cell bond order per unit volume. The R-, δ-, and γ-HMX phases have the values of 0.073, 0.057, and 0.078, respectively. Therefore, we may conclude that the β-HMX crystal has the strongest crystal bonding among the four polymorphs. This is consistent with the fact that β-HMX is the thermodynamically most stable form under ambient conditions among the HMX polymorphs. The crystal bonding for the HMX phases weakens in the following order: β-HMX > γ-HMX > R-HMX > δ-HMX, so the activation order for thermal decomposition of the four polymorphs is as follows: β-HMX < γ-HMX < R-HMX < δ-HMX. From these discussions, we may infer that the impact sensitivity for the HMX phases increases in the following sequence: β-HMX < γ-HMX < R-HMX < δ-HMX, which is in agreement with their experimental sensitivity order.38

4. Conclusions In this study, we have performed a detailed density functional theory study of the electronic structure and vibrational properties of the four HMX polymorphs in the local density approximation. Then an attempt has been made to correlate the impact sensitivity of the four phases with their electronic structure. From the density of states for the four HMX phases, it is found that the states of N in the ring make more important contributions to the valence bands than these of C and N of NO2 and so N in the ring acts as an active center. The β-HMX phase has clearly very different vibrational properties from R-, δ-, and γ-HMX. From the low frequency to high-frequency region, the molecular motions of the frequencies for the four HMX polymorphs present unique features, which could be used to distinguish the polymorphs easily from one another. It is also noted that there is a relationship between the band gap and impact sensitivity for the four HMX polymorphs. From the cell bond order per unit volume, we may infer the variation order of crystal bonding for the four polymorphs and so predict their impact sensitivity order. Acknowledgment. We thank the funding from Key Laboratory for Shock Wave and Detonation Physics. This work was partly supported by National Natural Science Foundation of China (Grant No. 10576016) and National “973” Project. References and Notes (1) Cooper, P. W.; Kurowski, S. R. Introduction to the Technology of ExplosiVes; Wiley: New York, 1996. (2) Akhaven, J. The Chemistry of ExplosiVes; Royal Society of Chemistry: Cambridge, U.K., 1998. (3) Cady, H. H.; Smith, L. C. Los Alamos Scientific Laboratory Report LAMS-2652 TID-4500; Los Alamos National Laboratory: Los Alamos, NM, 1961. (4) Main, P.; Cobbledick, R. E.; Small, R. W. H. Acta Crystallogr., Sect. C 1985, 41, 1351. (5) Skidmore, C. B.; Phillips, D. S.; Idar, D. J.; Son, S. F. In Conference Proceedings: Europyro 99 Vol. I; Association Francaise de Pyrotechnie: Brest, 1999; p 2. (6) Idar, D. J.; Lucht, R. A.; Straight, J. W.; Scammon, R. J.; Browning, R. V.; Middleditch, J.; Dienes, J. K.; Skidmore, C. B.; Buntain, G. A. In Proceedings of the EleVenth International Detonation Symposium; Snowmass Village, CO, Aug 31-Sept 4, 1999; Naval Surface Warfare Center: Indian Head, 1999; p 335. (7) Dlott, D. D.; Fayer, M. D. J. Chem. Phys. 1990, 92, 3798. (8) Tokmanoff, A.; Fayer, M. D.; Dlott, D. D. J. Phys. Chem. 1993, 97, 1901. (9) Tarver, C. M. J. Phys. Chem. A 1997, 101, 4845. (10) Goetz, F.; Brill, T. B.; Ferraro, J. R. J. Phys. Chem. 1978, 82, 1912. (11) Yoo, C.-S.; Cynn, H. J. Chem. Phys. 1999, 111, 10229. (12) Goetz, F.; Brill, T. B. J. Phys. Chem. 1979, 83, 340. (13) Iqbal, Z.; Bulusu, S.; Autera, J. R. J. Chem. Phys. 1974, 60, 221.

12722 J. Phys. Chem. B, Vol. 111, No. 44, 2007 (14) Brand, H. V.; Rabie, R. L.; Funk, D. J.; Diaz-Acosta, I.; Pulay, P.; Lippert, T. K. J. Phys. Chem. B 2002, 106, 10594. (15) Lewis, J. P.; Sewell, T. D.; Evans, R. B.; Voth, G. A. J. Phys. Chem. B 2000, 104, 1009. (16) Lewis, J. P. Chem. Phys. Lett. 2003, 371, 588. (17) Ye, S.; Koshi, M. J. Phys. Chem. B 2006, 110, 18515. (18) Payne, M. C.; Teter, M. P.; Allan, D. C.; Arias, T. A.; Joannopoulos, J. D. ReV. Mod. Phys. 1992, 64, 1045. (19) Vanderbilt, D. Phys. ReV. B 1990, 41, 7892. (20) Kresse, G.; Furthmu¨ller, J. Phys. ReV. B 1996, 54, 11169. (21) Fischer, T. H.; Almlof, J. J. Phys. Chem. 1992, 96, 9768. (22) Ceperley, D. M.; Alder, B. J. Phys. ReV. Lett. 1980, 45, 566. (23) Perdew, J. P.; Zunger, A. Phys. ReV. B 1981, 23, 5048. (24) Cady, H. H.; Larson, A. C.; Cromer, D. T. Acta Crystallogr. 1963, 16, 617. (25) Choi, C. S.; Boutin, H. P. Acta Crystallogr., Sect. B 1970, 26, 1235. (26) Cobbledick, R. E.; Small, R. W. H. Acta Crystallogr., Sect. B 1974, 30, 1918. (27) Sanchez-Portal, D.; Artacho, E.; Soler, J. M. Solid State Commun. 1995, 95, 685. (28) Segall, M. D.; Shah, R.; Pickard, C. J.; Payne, M. C. Phys. ReV. B 1996, 54, 16317. (29) Gonze, X. Phys. ReV. B 1997, 55, 10337. (30) Perdew, J. P.; Burke, K.; Ernzerhof, M. Phys. ReV. Lett. 1996, 77, 3865. (31) Perdew, J. P.; Chevary, J. A.; Vosko, S. H.; Jackson, K. A.; Pederson, M. R.; Singh, D. J.; Fiolhais, C. Phys. ReV. B 1992, 46, 6671. (32) Lewis, J. P.; Glaesemann, K. R.; VanOpdorp, K.; Voth, G. A. J. Phys. Chem. A 2000, 104, 11384.

Zhu et al. (33) Chakraborty, D.; Muller, R. P.; Dasgupta, S.; Goddard, W. A., III J. Phys. Chem. A 2000, 104, 2261. (34) Luty, T.; Ordon, P.; Eckhardt, C. J. J. Chem. Phys. 2002, 117, 1775. (35) Feng, S.; Li, T. J. Phys. Chem. A 2005, 109, 7258. (36) Stevens, L. L.; Haycraft, J. J.; Eckhardt, C. J. Cryst. Growth Des. 2005, 5, 2060. (37) Kirin, D.; Volovsek, V. J. Chem. Phys. 1997, 106, 9505. (38) McCrone, W. C. In Physics and Chemistry of the Organic Solid State; Fox, D., Labes, M. M., Wessberger, A., Eds.; Wiley: New York, 1965; Vol. II, p 726. (39) Xiao, H.-M.; Li, Y.-F. Sci. Chin. B 1995, 38, 538. (40) Xiao, H.-M.; Li, Y.-F. Banding and Electronic Structures of Metal Azides; Science Press: Beijing, 1996; p 88; in Chinese. (41) Zhu, W.; Xiao, H. J. Comput. Chem. 2007, in press. (42) Gilman, J. J. J. Appl. Phys. 1979, 50, 4059. (43) Gilman, J. J. Philos. Mag. Lett. 1998, 77, 79. (44) Gilman, J. J. Philos. Mag. B 1993, 67, 207. (45) Gilman, J. J. Mech. Mater. 1994, 17, 83. (46) Kuklja, M. M.; Stefanovich, E. V.; Kunz, A. B. J. Chem. Phys. 2000, 112, 3417. (47) Bowden, F. P.; Yoffe, A. D. Fast Reactions in Solids; Butterworths Scientific Publications: London, 1958. (48) Bowden, F. P.; Yoffe, A. D. Initiation and Growth of Explosion in Liquids and Solids; Cambridge University Press: Cambridge, 1952. (49) Bardo, R. D. Shock WaVes in Condensed Matter; Gupta, Y. M., Ed.; Plenum: New York, 1985; p 843. (50) Ching, W. Y.; Ouyang, L.; Yao, H.; Xu, Y. N. Phys. ReV. B 2004, 70, 085105.