Arxiv:1406.3917V1 [Cond-Mat.Mtrl-Sci]

Home , Silver fulminate

arXiv:1406.3917v1 [cond-mat.mtrl-sci] 16 Jun 2014 uto ftogssC n N and CO gases two of duction n ufliae aesrne ntaigpwrthan power initiating stronger have Hg(CNO) fulminates Cu and tn rpte fvrosfliae eeeaie by examined were wohler fulminates and various Martin crackers of Christmas propeties in ating used is it ml uniyi sd nprtcnc,fiewrs and works, fire pyrotechnics, in pistols used) toy is quantity small sa olw:2AgCNO follows: as is fACOi ucett e aiu uptfrom output maximum get (RDX)) to 1,3,5-trinitro-1,3,5-triazine sufficient Hexahydro the is AgCNO of ehlntomn n 0 oeecetfrexploding for efficient tri-nitro-xylene more 30% and methyl-nitroamine nteIainNavy Italian the in o ogtm ihgo eoaigproperties detonating good with time long a for saotfite ie seceta h ecr ful- mercury the as efficient as times (Hg(CNO) minate fifteen about is ie yahtwire hot a by nited neetv eoao n lofudt ea effective as is be which (Pb(N to azide found fulminates, lead also inorganic and as of detonator class effective an azides, the corresponding to the belongs cyanamides with explosives. and electronic secondary cyanates iso for are initiators They as applications find vrg eoainvlct fACOsml fthick- of sample AgCNO ness of velocity detonation average oyopimadtemdnmcgon tt fSle fulm Silver of state ground thermodynamic and Polymorphism nrai umntsaepiayepoie n they and explosives primary are fulminates Inorganic ∼ . mi nteodro 70ms 1700 of order the in is mm 0.5 2 6 oee,Hg(CNO) However, . h eopsto ehns fAgCNO of mechanism decomposition The . etepeerdtemdnmceulbimpaeudrstu under phase equilibrium Hitherto thermodynamic preferred the be ae norttleeg acltosuigDTD method, DFT-D2 using calculations energy total i our methods on correction Based dispersion the semiempirical predicting in by role captured relat Densi significant a on on play based interactions temperature Waals approaches and pseudopotential pressure planewave of using effect the investigated rhrobc( orthorhombic nue oyopi hs rniini enuigDTfunc DFT using seen is transition phase polymorphic induced h acltdeetoi tutr hw that shows structure compressibili Tran-Bla electronic the and calculated within prope state The method optical wave of plane and augmented structure equation linearized Electronic The investigated. were method. DFT-D2 with ASnumbers: PACS pcr.Oeal h rsn td rvdsa al indic early an provides study unstable present of stat formation the polymorph natur of the Overall, The of density light partial respectively. to spectra. sensitivity and eV, photo 4.43 plots and density and transitions charge 3.51 the of values through gap band a ivrfliae(gN)i rmr xlsv,wihexis which explosive, primary a is (AgCNO) fulminate Silver 3 .INTRODUCTION I. ti nw ob estv explosive sensitive a be to known is It . 2 8–10 o xldn 2,4,6-tri-nitro-phenyl- exploding for ) 4 5 gN eoaoshsbe used been has detonators AgCNO . tfid oeapiain (when applications some finds It . 3 Cmcm ) n hyfudta h g Cd Ag, the that found they and , 2 → nispr om(x 0mg 60 (Ex. form pure its in ) dacdCnr fRsac nHg nryMtras(ACRHE Materials Energy High in Research of Centre Advanced A C N + 2CO + 2Ag 1 Cmcm ivrfliae(AgCNO) fulminate Silver . and e al est ucinlcalculations functional density Waals der 2 ahbwi yeaa-5006 nhapaeh India pradesh, Andhra 046, 500 Hyderabad- Gachibowli, R n rgnl( trigonal and ) iedsic os when noise distinct give 7 β 2 hssdntdas denoted phases 3 ¯ oevr h initi- the Moreover, . fr fAgCNO. of -form nvriyo yeaa,Po.C .RoRoad, Rao R. C. Prof. Hyderabad, of University tnsfis hc is which first stands .Ydknau n .Vaitheeswaran G. and Yedukondalu, N. 2 n h pro- the and ; − 1 hnig- when , Dtd uut2,2018) 27, August (Dated: R )frsa min odtos ntepeetsuy ehave we study, present the In conditions. ambient at forms 3) ¯ 2,3 2 the , It . α and α , β β frso gN,rsetvl.As pressure a Also respectively. AgCNO, of -forms umnts esrdfiueo nestvns Fo )to I) Hg(CNO) of (F for insensitiveness the impact of to figure compared Measured insensitive fulminates. and stability thermal better o en ieyue ssadr rmr explosive, Pb(N primary by replaced standard largely as was it used but widely being now omuega oordet ogepsr olgtand light nature to hygroscopic exposure non long possess transform to crystals crystals these due white colour crys- the mouse-gray fresh that been to like exception as has an properties with it similar tals about However, exhibit water ago under years stored in]. 40 crystals there AgCNO that references reported and 2 [Ref 3)adtmeaueo nnto To )frAgCNO for I) of (T ingnition of (170 temperature and (30) Pb(N The toxicity. and stability n ae eonsdb Gay-Lussac by recognised later and a rtpitdotb Pandey by out AgCNO of pointed phases first ploymorphic two was of formation and try oycnrdotohmi tutr aigsaegroup space Ibam having structure orthorhombic centred body Singh 6.04 = c h w oyopi hssaerpre sorthorhombic of as ( structures reported are crystal phases polymorphic redetermined two the the Further, (1965). aaeesa=3.88 = a parameters n ,adtioa ( trigonal and 4, = Z and 9.087 = a nlssstudies analysis mal Cmcm hssaeidrc adgpisltr with insulators gap band indirect are phases gN a rtpeae yBrugnatelli by prepared first was AgCNO o C) hs tblt n hycnb effectively be can they and stability phase r nlzdfo h acltdoptical calculated the from analyzed are s 15 te eecluae sn ulpotential full using calculated were rties ihltiecntnsa=3.88 = a constants lattice with cnrs osadr F functionals. DFT standard to ncontrast hs yBarrick by phase ) A ˚ to oeprmnait oaodthe avoid to experimentalists to ation eotdtecytlsrcueo gN obe to AgCNO of structure crystal the reported s pia nstoy electric-dipole anisotropy, Optical es. yFntoa hoy(F) a der van (DFT). polymorphs Theory the Functional of ty stability phase ive inl n h aewsntobserved not was same the and tionals idpesr n eprtr range. temperature and pressure died < amdfidBceJhsnpotential. Becke-Johnson modified ha , si w oyopi hssnamely phases polymorphic two in ts A ˚ α n .Teacat rsa symme- crystal accuarte The 4. = Z and , Hg(CNO) fceia odn sanalyzed is bonding chemical of e 115.73 = the yo ohplmrhcphases polymorphic both of ty Cmcm ∗ 20 2 A ˚ 2 ◦ 10.752 = b , niaen vdnefrtemper- for evidence no indicate (10) n .Tedffrnilther- differential The 6. = Z and (210 tutr sfudto found is structure R nt tde rmvan from studied inate )paeb Britton by phase 3) ¯ M), tal et o < C) gN (22) AgCNO 18 < 3 16 ) 17) ihunit-cell with (1979), 2 3 Pb(N A ˚ 14 n Britton and ) A ˚ a on ohave to found was 2 11.20 = b , 11 n 5.804 = c and 12) n1959, In (1824). 3 u oispoor its to due ) 2 19 < 2,12 ( > 13 19) a (1991), Pb(N . 300 (1798) tal et A ˚ and o 3 C) A ˚ 17 ) 2 , 2

o o ature induced phase changes between 20 C and 145 C. only on the chemical species, and Rij is an interatomic 1 To the best of our knowledge, there are no studies to distance. fdamp = − − is a damping function 1+e d(Rij /R0 1) explore the relative stability under pressure and temper- which is necessary to avoid divergence for small values ature, electronic structure, and chemical bonding of the of Rij and R0 is the sum of atomic vdW radii. The two polymorphic phases in the literature from both the- DFT-D2 method has been applied to several molecular oretical and experimental point of view for the energetic systems such as secondary explosives31–34 and inorganic solid AgCNO. Therefore, we systematically investigated azides35–37, and these results demostrated that disper- the above mentioned properties for AgCNO using ﬁrst sion corrections are important in determining the ground principles calculations based on Density Functional The- state properties. Hence, we attempted to study the phase ory (DFT). The rest of the article is organized as follows, stability of AgCNO using the DFT-D2 method and the in section II, we brieﬂy describe the methodology of our results are discussed in the following section. calculations. In section III, the polymorphic phase sta- The thermodynamic stabilities were investigated by bility, electronic struture, chemical bonding, and optical comparing the free energy of the polymorphs. In order to properties of AgCNO are discussed. Finally, in section predict the most stable form of AgCNO, we have calcu- IV, we summarize the results, which concludes our paper. lated the Gibbs free energy of both polymorphic phases. The Gibbs free energy for any material at a given temperature is as follows. II. METHODOLOGY OF CALCULATIONS

G(P,T )= F (V,T )+ P V = Fvib + Fperfect + P V (3) The polymorphic phases of AgCNO crystals were studied using periodic DFT calculations performed with two Where F, P, T, and V are the Helmoltz free energy, different codes such as Plane Wave Self Consistent Field pressure, absolute temperature and volume, respectively. (PWSCF)21 and Cambridge series of total energy pack- F(V,T) is summation of vibrational free energy and per- 22 age (CASTEP) , which are based on planewave pseu- fect lattice energy i.e F(V, T) = Fvib + Fperfect. The dopotential (PW-PP) approach. We have used Van- vibrational free energy Fvib is calculated within the har- derbilt type ultrasoft pseudo potentials23 to treat ion- monic approximation and is given by electron interactions, while electron-electron interactions − hω hωi are treated with the generalized gradient approximation i K T F = K T + ln(1 − e B ) (4) 24 vib B X 2K T (GGA) by Perdew-Burke-Ernzerhof (PBE) and Perdew i B and Wang (PW91)25. The Broyden-Fletcher-Goldfarb- 26 Shanno (BFGS) minimization scheme was used in ge- Where ωi is the phonon frequency, h is the Planck con- ometry optimization. The plane wave cutoff energy was stant and KB is the Boltzmann constant. Where Fperfect set to 60 Ry and a k-point grid of 7×7×4 and 5×5×5 =E0 +Eel - TSel,E0 and Eel- TSel are the contributions for Cmcm and R3¯ phases, respectively using Monkhorst- from the lattice and electronic excitations, respectively. Pack grid scheme27. The self-consistent energy conver- If the electronic excitations and pressure changes are ne- − gence accuracy was set to 1.0×10 6 eV/atom and the glected then the Gibbs free energy can be calculated by −4 force to be 1.0×10 eV/A. G=E0 + Fvib. The standard LDA/GGA functionals are unabale to describe the non-covalent interactions in molecular crystals. On the other hand, semiemipirical dispersion cor- III. RESULTS AND DISCUSSION rection methods such as the Ortmann, Bechstedt, and 28 Schmidt (DFT-OBS) correction to PW91, as well as A. Crystal structure and phase stability of 29 the Tkatchenko and Scheffler (DFT-TS) and Grimme polymorphs (DFT-D2)30 corrections to PBE were employed to handle vdW interactions. These methods are found to improve Extreme sensitivity, poor stability and high cost of Sil- the results to a greater extent. Here the total energy ver as a raw material prohibited AgCNO in commercial after inclusion of vdW correction is given by or military priming and detonating devices38,39 but still it finds applications in pyrotechincs, fire works, toy pistols EDF T +D = EDF T + Edisp (1) and in the navy. However, the relative stability and nec- Where EDF T is the self-consistent Kohn-Sham energy, essary conditions for the formation of two polymorphic Edisp is empirical dispersion correction given by phases are still unknown for this material. As discussed in section I, AgCNO exists in two polymorphic phases C E = −S ij f (R ) (2) at ambient conditions, and it is noticed that except the disp 6 X R6 damp ij 17–19 i

As a first step, we have performed full structural op- posed that Cmcm is stable than R3¯ phase due to the timization by taking experimental data as a input.18,19 presence of linear form of CNO anions in their crys- The lattice parameters and fractional co-ordinates of tal structures.45 The predicted behavior with DFT-D2 both the phases are allowed to relax to get the most method are in good agreement with the experimental stable configurations at ambient pressure. The opti- observations17–19,45 and this clearly shows the success of mized crystal structures of orthorhombic (Cmcm) and DFT-D2 method in predicting the relative phase stability rhombohedral (R3)¯ phases are shown in Fig. 1. The of the ploymorphic phases of AgCNO. crystal structures of two polymorphic phases mainly differ by planar zigzag chains (orthorhombic phase in yz plane) and cyclic hexamers (rhombohedral phase in xy B. Pressure and temperature effects on the 17 plane) (see Fig. 1) and these AgCNO molecules are polymorphs bind through weak vdW forces within the unitcell. We have calculated the structural properties and their pressure dependence within standard DFT functionals using The investigation of energetic materials at extreme conditions is a challenging task because of their sensi- PWSCF and CASTEP codes. The obtained volumes for Cmcm and R3¯ phases are overestimated by 22.5% and tivity, complex chemical behavior and risk of decomposition. Several accidents occured during the preparation 13.5% with PBE-GGA, respectively. This clearly rep- 2 resents that PBE-GGA functional is inadequate to pre- of AgCNO and hence it is very difficult to perform experiments without any prior knowledge about this kind dict the ground state properties of this energetic molec- of energetic materials. Therefore, theoretical modeling ular solid. Therefore, we have used various dispersion corrected (DFT-D) methods to capture vdW interac- and simulations are efficient tools to predict the physical and chemical properties of complex energetic solids tions to reproduce the ground state properties comparable with the experiments.17–19 The obtained volumes are at extreme conditions and also very useful in predicting results a priori to the experimentalists. DFT is a power- overestimated by about 12.6%, 7.8% using DFT-OBS; 3.7%, 4.2% using DFT-TS; and 2.2%, 4.5% using DFT- ful tool in predicting the behavior of complex solid state systems at extreme conditions. Hence, we attempted to D2 methods for Cmcm, R3¯ phases, respectively. Among the three DFT-D methods, DFT-TS and DFT-D2 meth- study the relative phase stability and possible polymorphic structural transition in AgCNO under hydrostatic ods work well for the molecular solid AgCNO, still there are dicrepancies between the theoretical values at 0 K and pressure upto 5 GPa with a step size of 0.5 GPa. As illustrated in Fig. 2, The calculated enthalpy curves as experimental data at 298 K17–19. The order of discrepan- cies about ∼3-4% are previously reported for seconadry a function of pressure using standard DFT functionals show that R3¯ phase undergoes a structural ploymorphic explosive molecular crystals with DFT-D methods at 0 Cmcm K32,40,41 and these errors were further reduced by inclu- phase transition to phase at about 2.5 GPa within GGA using PWSCF and the corresponding transition sion of thermal effects in the calculations40,42 and the obtained results are in excellent agreement with the ex- pressures obtained from CASTEP code is 2.7 GPa (see Fig. 1 of supplementary material61). Overall, we ob- perimental data. In the present study, the equilibrium volumes are overestimated within static DFT-D2, if we served a similar trends using both of the approaches with small deviations in the calculated transition pressures. include the temperature effects there will be an incre- Cmcm ment in the volumes due to thermal expansion which Above the transition pressures the enthalpy of the phase is lower than that of the R3¯ phase, indicating that leads to large dicrepancies over static DFT-D2 results Cmcm structure becomes stable. The Cmcm phase is unless the polymorphs have negative thermal expansion. The calculated ground state lattice parameters, unitcell stable in the pressure range from transition pressure to 5 GPa, which is the highest pressure we accomplished in volume, and density of both polymorphic phases using various dispersion corrected DFT functionals are com- the present work. The enthalpy (total energy) difference at 0 GPa is ∼ 6 KJ/mole which increases with pressure pared with experimental data and are presented in Table ∼ I. The obatined ground state properties with and with- and it is found to be 21 KJ/mole at 5 GPa within DFT-D2 method. In contrast to LDA/GGA function- out dispersion correction methods using PWSCF and CASTEP are in good accord with each other (see Ta- als, the DFT-D2 method using both of the approaches (PWSCF and CASTEP) reveals that there is no pressure ble I). In general, the energy difference between different polymorphic forms are mostly in the order of 0-10 induced ploymorphic phase transition seen in AgCNO and the Cmcm phase is found to be the most stable KJ/mole.43 The DFT-D2 method is successful in predicting the polymorphs of Benzamide (P1 and P3) with form of AgCNO over studied pressure range using both PW-PP approaches (see Fig. 2). Hitherto, we represent an energy difference by 1.9 KJ/mole.44 Hence, we have Cmcm R¯ α β used DFT-D2 method for further calculations. The ob- and 3 phases as and -AgCNO, respectively. Our recent high pressure study on KClO 46 also revealed tained total energy difference for Cmcm phase lowered 3 by ∼6 KJ/mole per molecule than R3¯ phase which im- that DFT-D2 method is good enough for calculating the transition pressures for molecular solids in opposition to plies that Cmcm phase found to be the thermodynamic ground state at ambient pressure. Britton et al17 pro- standard LDA/GGA functionals. As shown in Fig. 3 (see Fig. 2 of supplementary material61), the calculated total 4 energies ploted as a function of pressure demonstrates efficients under pressure for β-AgCNO. This clearly in- that the DFT-D2 functional predict the α-phase to be dicates that the polymorphs behave anisotropically un- the most stable form of AgCNO over studied pressure der the hydrostatic pressure. Further, to understand the range in contrast to usual LDA/GGA functionals. relative compressibilities of both the polymorphic phases In addition, the contribution of lattice vibration needs bondlengths and bond angles as a function of hydrostatic to be considered at elevated temperatures. We have used pressure are plotted within DFT-D2 method. As shown norm-conserving pseudo potentials47 for calculating the in Fig. 7, Ag-O bond is more compressible in both the Gibbs free energy as they are well suited for phonon cal- phases along a-crystallographic direction whereas Ag-Ag culations as implemented in CASTEP code. Fig. 4 shows is the most compressible in β-AgCNO. The bonds such as the Gibbs free energy as a function of temperature in the Ag-C, C-N, and N-O show similar and less compressible range of 0-450 K within DFT-D2 method. The difference nature over studied pressure range due to strong covalent in Gibbs free energy (∆G) is ∼ 7 KJ/mole at low temper- bonding between Ag and C as well as within the CNO ature region and it decreases monotonically with increase anion molecule and this can be clearly understood from in temperature between the two α- and β-polymorphs their electronic structure. Similarly the corresponding of AgCNO. However, α-AgCNO remains with the low- bond angles Ag-C-Ag decrease because of the orienta- est free energy throughout the investigated temperature tion of Ag-C bonds due to high compressibility nature range, despite the contribution of lattice vibration tend- of Ag-Ag bonds (see Fig. 7) whereas Ag-C-N bond an- ing to decrease the ∆G gap at elevated temperatures. We gle increases under compression (see Fig. 8). Overall, have not observed any temperature driven phase transi- we observe that Ag-C, C-N, and N-O bonds are stiffer, tion between two polymorphs over studied temperature while Ag-Ag and Ag-O are more compressible under the range and this is consistent with results of differential application of hydrostatic pressure. Further, this com- thermal analysis20. Gibbs free energies at temperatures pressibility behavior can be clearly understand by ana- beyond 450 K are not plotted because α-AgCNO phase lyzing the nature of chemical bonding in AgCNO and typically decompose at higher temperatures between 393- hence we have investigated the electronic struture and 450 K20. From our present study, we confirmed that α- chemical bonding of both the polymorphic phases in the phase is the most stable polymorph of AgCNO under following section. studied pressure and temperature range.

D. Electronic structure and chemical bonding C. Equation of state and compressibilities of the polymorphs Inorganic fulminates are isoelectronic with azides, cynates, and cynamides. Iqbal and his co-workers1 ex- To study the effect of hydrostatic pressure on crystal plained that these fulminates are found to be the most structures of both the polymorphs, we have used vari- sensitive explosives to shock and heat among the ener- able cell optimization technique. The calculated volumes getic materials such as azides, cynates, thiocynates, and decrease monotonically as a function of pressure using cynamides due to asymmetric charge distribution in their standard LDA/GGA functionals as well as with DFT- structures. Also, the heavy metal complex salts are un- D2 method are shown in Fig. 5 (see Fig. 3 of supple- stable than alkali metal salts because of the asymmet- mentary material61) and the obtained equilibrium bulk ric inter ionic distances45. In order to understand the moduli for α and β-AgCNO are found to be 20.0 (19.4) macroscopic energetic behavior of these sensitive explo- GPa and 17.8 (19.1) GPa within DFT-D2 method us- sive materials, it is necessary to understand the elec- ing PWSCF (CASTEP) by fitting pressure-volume data tronic structure, and chemical bonding of these mate- to second order Birch-Murnaghan equation of state.48 In rials at the microscopic level. Iqbal et al1 carried out a order to understand the behavior of unit-cell parameters detailed experimental study on the electronic structure and their relative compressibilities under compression, we and stability of the inorganic fulminates, which reveals have presented the lattice constants as a function of pres- that sodium, potasium, and thallous fulminates to be sure in Fig. 6. The calculated lattice constant ’a’ de- ionic salts whereas silver and mercury salts are covalent creases monotonically with pressure whereas the b and in nature and this will be reflected in the order of stabil- c lattice constants are decreasing non-monotonically in ity of the fulminate salts. However, the electronic struc- the pressure range of 2.5 to 5 GPa. The computed lat- ture and bonding properties of the polymorphs of silver tice constants are very well reproduced with fifth-degree fulminate salts are not reported in the literature. There- polynamials and the fitting pressure coefficients are given fore in the present study we made a detailed analysis in Table II as previosly reported for explosive nitrate es- of electronic structure and chemical bonding of both the ter 1 (NEST-1).49 As illustrated in Fig. 6 and from the polymorphic phases at ambient pressure. calculated first order pressure coefficients, the c-axis and It is well known fact that the standard LDA/GGA b-axis are the most and least compressible, respectively functionals severely underestimate the band gaps for for α-AgCNO while rhombohedral lattice constant (a) semiconductors and insulators. In order to get rela- and angle (α) are decreasing with distinct pressure co- iable energy band gaps several methods have been pro- 5 posed such as LDA+U, LDA+DMFT, hybrid fuctionals, is a charge sharing between C-p and Ag-d states repre- and GW approximation. Unfortunately these methods senting the covalent nature of Ag-C bond. The nature are computationally very expensive. However, Tran and of chemical bond between two atoms can be predicted Blaha modified Becke Johnson (TB-mBJ) potential50 is from their electronegativity difference. When the differ- computatinally less expensive method, which provides ences are greater or equal to 1.7 (ionic) and less than accurate energy band gaps as comparable with more so- 1.7 (covalent), if the difference is greater than 0.5 the phisticated methods within the Kohn-Sham frame work. covalent bond has some degree of polarity. According Recently, Koller et al 51,52 reported the merits and limi- to Pauling scale the electronegetivity values Ag (1.9), C tations of the TB-mBJ functional. So far, several groups (2.5), N (3.0), and O (3.4) shows that there exists a co- used the TB-mBJ potential for the calculation of elec- valent bond between N-O, C-N and Ag-C bonds. The tronic structure and optical properties of different kinds electronegetivity difference in N-O (0.4), C-N (0.5) and of materials53–55 and confirmed that bandgaps obtained Ag-C (0.6) exhibits strong and polar covalent nature of using TB-mBJ potential are improved when compared to N-O, C-N and Ag-C bonds, respectively. The N-O, C-N LDA/GGA functional for a wide range of materials (see and Ag-C bonds show less compressibility behavior with Fig. 1 from Ref. 56). The accurate prediction of band increasing pressure (see Fig. 7), this is due to strong hy- gaps using TB-mBJ functional for diverse materials mo- bridization between Ag-d and s, p-states of C, N, and O tivated us to use this functional for high energetic mate- atoms leads to strong covalent character. Further, this rials. This semi local functional is implemented through can be clearly understood from electronic charge density WIEN2k package57. Since all these energetic materials plots which are used for accurate description of chemical are semiconductors or insulators, we have used the TB- bonds.58 The calculated valence charge density using TB- mBJ functional to calculate energy band gap values for mBJ functional for both the polymorphs are as shown in the two ploymorphs. We have used experimental crystal Fig. 11. It shows anisotropic bonding interactions and structures18,19 to calculate the electronic structure and the charge cloud is distributed within the CNO molecule optical properties of the both polymorphs at ambient indicating covalent character. Overall, the C, N, and O pressure. The calculated band gaps are found to be 3.51 atoms are covalenlty bonded within CNO anion and the and 4.43 eV for α and β-phases of AgCNO, respectively metal atom is also covalently bonded with CNO group and the respective band gaps using LDA functional are through C atom. X-ray electron spectroscopy59 study on 2.00, 2.85 eV. The obtained TB-mBJ band gap of 3.51 inorganic azides reveals that Heavy Metal azides (HMAs) eV for α-AgCNO is closely comparable with the mea- are more covalent than Alkali metal azides (AMAs), im- sured optical energy gap of 4.0 eV1 over LDA functional. plying that HMAs are more sensitive than AMAs and Also, the calculated TB-mBJ band structures show that find applications as intiators for secondary explosives. α and β-AgCNO phases are indirect band gap insulators Therefore, the presence of covalent bonding in AgCNO along S-(Γ-Z) and T-(F-Γ) high symmetry directions of makes it more sensitive than the above mentioned ionic the Brillouin zone as shown in Fig. 9. inorganic fulminates and it can be used as an initiator as HMAs. Further, the nature of chemical bonding in the two polymorphic phases was investigated by examining the total and partial density of states (DOS) as illustrated in Fig. 10. The relative phase stability of the α and β E. Optical properties polymorphs of AgCNO can also be explained on the basis of the DOS. As shown in Fig. 10, when compared to Energetic materials become unstable and they undergo β-AgCNO the peaks of the DOS in the valence bands of photochemical decomposition by the action of light. α-AgCNO has a tendency to shift towards lower energy, Hence it is interesting to study the optical properties of indicating better stability of the α polymorphic phase. these materials to understand the decomposition mecha- Also, it can be seen that total DOS of both phases ex- nisms. Electronic structure calculations could provide an hibit some similar features, hence we have analyzed the information about the nature and location of interband bonding based on PDOS of α-AgCNO. It is clearly vi- transitions in crystals. The complex dielectric function sualized from the DOS that AgCNO polymorphs have ǫ(ω) = ǫ1(ω) + iǫ2(ω) can be used to describe the lin- molecular character, which arises from a strong overlap ear response of the system to electromagnetic radiation contributions of partial DOS of Silver and Fulminate an- which is related to the interaction of photons with elec- ions in the valence and conduction bands and this can trons. The imaginary part of dielectric function ǫ2(ω) is be commonly seen in molecular crystalline solids. The obtained from the momentum matrix elements between conduction band is mainly due to p-states of C, N, O the occupied and unoccupied wavefunctions within selec- and s,p,d-states of metal (Ag) atoms. The lowest ly- tion rules. The optical properties of AgCNO polymorphs ing states about -9 eV are due to hybridized C-s and were calculated using TB-mBJ functional at a denser k- s, p states of N and O atoms. The states from -5 to mesh of 20 × 20 × 12 and 17 × 17 × 17 for α- and -7 eV are due to hybridization of Ag-d and p-sates of β- phases, respectively. The AgCNO polymorphs crys- C, N and O atoms. Also, the top of the valence band tallize in the orthorhombic and hexgonal symmetry, and is mainly dominated by O-p and Ag-d states and there this allow non zero components of the dielectric tensors 6 three for α and two for β phase along [100], [010] and IV. CONCLUSIONS [001] directions. Fig. 12 shows the real ǫ1(ω) and imaginary ǫ2(ω) part of the dielectric function as a function In summary, ab-initio calculations have been per- of photon energy for both the polymorphs. The peaks in formed to investigate the relative phase stability, struc- ǫ2(ω) mainly arises due to electric-dipole transitions be- tural transition and compressbilities of both polymor- tween valence and conduction bands. The major peaks phic orthorhombic (Cmcm) and trigonal (R3)¯ phases in ǫ2(ω) at 4.57 eV (along 100), 5.94 eV (along 010), 5.06 of AgCNO using PW-PP approaches. The standard eV (along 001) in α-phase and 5.45 eV (along 100), 6.13 LDA/GGA functionals are inadequate to predict the rel- eV (along 001) in β-phase arise due to interband tran- ative phase stability and ground state properties of these → sitions from Ag(4d) N(2p)-states. It can be clearly polymorpic phases. From our total energy calculations seen from DOS, the top of valnce band is mainly domi- using DFT-D2 method, Cmcm phase was found to be nated by 4d-states of Silver atom (see Fig. 10) in both preferred thermodynamic equilibrium phase under pres- the polymorphs, which is similar to the case of Silver sure and temperature range and hence we confirmed 60 azide (AgN3) . As discussed in the section I, AgCNO Cmcm and R3¯ phases as α and β-AgCNO, respectively. is iso-electronic and it consists of same chemical species The present study also reveals that there is no structural as AgN3 in addition to C and O atoms, hence one can phase transition observed using DFT-D2 method whereas expect that major optical transitions arise from Ag(4d) we observe a pressure induced polymorphic phase tran- → N(2p)-states in both energetic materials AgCNO and sition from β → α at about 2.5 (2.7) GPa within GGA 60 AgN3 . using PWSCF (CASTEP) package. Overall, the phase stability and ground state properties obtained from both the PW-PP approaches are in good accord with each other for both of the polymorphs. We have also calculated the electronic structure and density of states of both polymorphic phases using TB-mBJ functional at ambient The real part ǫ1(ω) of can be derived from ǫ2(ω) using the Kramer-Kronig relations. The calculated real static pressure. The calculated electronic band structure show dielectric constant along all the crystallographic direc- that α, β phases are indirect band gap insulators with a tions are 2.24, 5.01, 2.78 for α-phase and 2.59, 4.38 for band gap values of 3.51 and 4.43 eV, respectively. The β-phase. The static refractive indices calculated using detailed analysis of electronic charge density distribution and partial density of states reveal covalent bonding is the dielectric function n = ǫ(0) are given by n100 = p predominant in the energetic AgCNO polymorphs. The 1.50, n010 = 2.24, n001 = 1.67 and n100 = 1.61, n001 = 2.09 for α and β-phases of AgCNO, respectively. The major peaks in ǫ2(ω) are due to interband transitions → obtained refractive indices are distinct in all the crys- between Ag(4d) N(2p)-states. The calculated absorp- tallographic directions, which indicates the anisotropy of tion spectra reveal that the polymorphs show considerable anisotropy and decompose under the action of UV the α and β polymorphs. Using ǫ1, ǫ2 one can dervive the optical constants such as absorption and photoconductiv- light. Finally, the present study suggests that the α-form ity of the materials. The calculated absorption spectra is is stable and it can be used in different applications and shown in Fig. 13 and absorption starts after the energy β-form can be avoided by careful control of the pressure 3.51 eV in α-phase, 4.43 eV in β-phase, which is energy or temperature in the manufacturing process. band gap between the valence band maximum and conduction band minimum. The first absorption peaks along V. ACKNOWLEDGMENTS three (100, 010 and 001) directions are at 4.5 eV, 7.7 eV, 7.9 eV and the corresponding absorption coefficients are found to be 1.1 × 107m−1, 1.4 × 108m−1, 3.9 × 108m−1 Authors would like to thank Defence Research for α-phase, while the same are 1.3 × 108m−1 at 6.8 eV and Development Organisation (DRDO) through and 4.0 × 107m−1 at 7.2 eV along two (100 and 001) ACRHEM for the financial support under grant No. directions for β-phase. Photoconductivity is due to the DRDO/02/0201/2011/00060:ACREHM-PHASE-II, and increase in the number of free carriers when photons are the CMSD, University of Hyderabad, for providing absorbed. The calculated photoconductivity shows (see computational facilities. NYK acknowledges Dr. V. Fig. 13) a wide photocurrent response in the absorption Kanchana, Department of Physics, Indian Institute region of 3.5-35 eV and 4.4-35 eV in α and β phases of of Technology Hyderabad for critical reading of the manuscript. AgCNO, respectively. overall, the energetic polymorphs ∗ show a strong anisotropic and wide range of absorption Author for Correspondence, E-mail: (absorption coefficients ∼ 107m−1). This results indicate [email protected] the possible decomposition of AgCNO into Ag, CO and N2 under the action of ultraviolet (UV) light. Therefore AgCNO polymorphs decompose under the action of UV light and they may explode due to photochemical decomposition. 7

1 Z. Iqbal and A. D. Yoffe, Proc. R. Soc. Lond. A, 302, 35 Chem. Phys 140, 014105 (2014). (1967). 35 K. Ramesh Babu, G. Vaitheeswaran, Chem. Phys. Lett., 2 P. H. Collins, K. J. Holloway, Propellants and explosives, 533, 35 (2012). 3, 159 (1978). 36 K. Ramesh Babu, G. Vaitheeswaran, J. Chem. Sci. 124, 3 T. L. Davis, The chemistry of powder and explosives (Wi- 1391 (2012). ley, New York 1943). 37 K. Ramesh Babu, G. Vaitheeswaran, Chem. Phys. Lett. 4 F. P. Bowden and H. T. Williams, Proc. R. Soc. Lond. A, 586, 44 (2013). 208, 176 (1951). 38 G. W. C. Taylor, and J. M. Jenkins, Proc. 3rd Simposium 5 T. Urbanski, Chemistry and Technology of explosives on Chemical problems connected with the stability of Ex- (Pergamon Press, Oxford, 1967), vol. 3, p. 157. plosives (Ystad, Sweden, May 28-30, 33, 1973). 6 J. P. Springer, Ger Pat 953, 415 (1956). 39 C. A. Taylor and E. P. Buxton, Army Ordnance, 6, 118 7 D. A. Jackson and A. P. Dicks, J. Chem. Educ., 89, 1267 (1925). (2012). 40 A. C. Landerville, M. W. Conroy, M. M. Budzevich, Y. Lin, 8 F. Martin, Dber Azide und Fulminate Darmstadt (1913). C. T. White, I. I. Oleynik, Appl. Phys. Lett. 97, 251908 9 L. Wohler and F. Martin, Angezu. Chem. 30, 33 (1917). (2010). 10 L. Wohler and F. Martin, Ber. Dtsch. Chem. Ges. 50, 586 41 F. Shimojo, Z. Wu, R. K. Kalia, and P. Vashishta, J. Chem. (1917). Phys., 132, 094106 (2010). 11 F. D. Miles, A history of the Research in Nobel division of 42 Z. Wu, R. K. Kalia, A. Nakano and P. Vashishta, J. Chem. I.C.I (Nobel division, ICI, Birmingham, 1955). Phys., 134, 204509 (2011). 12 A. M. Peter, J. Am. Chem. Soc., 38, 486 (1916). 43 D. Baraga, and F. Grepioni, Making Crystals by Design: 13 B. T. Fedoroff, and O. E. Sheffield, Encyclopedia of Explo- Methods, techinques and Applications, Willey, 2006. sives and Related Items (PATR 2700, U.S ARDEC, New 44 P. Ectors, and D. Zahn, Phys. Chem. Chem. Phys. 15, Jersey, 1974) vol. 6. 9219 (2013). 14 J. L. Gay -Lussac, Ann. Chim. Phys. 27, 199 (1824). 45 Z. Iqbal, Structure and Bonding, Springer, 10, 25 (1972). 15 K. Singh, Acta Cryst, 12, 1053 (1959). 46 N. Yedukondalu, Vikas D. Ghule and G. Vaitheeswaran, J. 16 S. N. Pandey, Indian J. Phys. 36, 657 (1962). Chem. Phys., 138, 174701 (2013). 17 D. Britton, and J. D. Dunitz, Acta. Cryst. 19, 662 (1965). 47 N. Troullier and J. L. Martins, Phys. Rev. B, 43, 1993 18 J. C. Barrick, D. Canfield, and B. C. Giessen, Acta. Cryst. (1991). B35, 464 (1979). 48 F. D. Murnaghan, Proc. Natl. Acad. Sci. USA 30, 244 19 D. Britton, Acta Cryst. C47, 2646 (1991). (1944). 20 T. Boddington and Z. Iqbal, Trans. Faraday Soc., 65, 509 49 A. C. Landerville, M. W. Conroy, I. I. Oleynik, C. T. (1969). White, J. Phys. Chem. Lett. 1, 346 (2010). 21 P. Giannozzi et al, J. Phys.: Condens. Matter, 21 395502 50 F. Tran, and P. Blaha, Phys. Rev. Lett. 102, 226401 (2009). (2009). 22 M. C. Payne, M. P. Teter, D. C. Allen, T. A. Arias and J. 51 D. Koller, F. Tran and P. Blaha, Phys. Rev. B 85, 155109 D. Joannopoulos, Rev. Mod. Phys. 64, 1045 (1992). (2012). 23 D. Vanderbilt, Phys. Rev. B, 41, 7892 (1990). 52 D. Koller, F. Tran and P. Blaha, Phys. Rev. B 83, 195134 24 J. P. Perdew, S. Burke and M. Ernzerhof, Phys. Rev. Lett. (2011). 77, 3865 (1996). 53 D. J. Singh, Phys. Rev. B, 82, 205102 (2010). 25 J. P. Perdew, and Y. Wang, Phys. Rev. B 45, 13244 (1992). 54 J. A. Camargo-Martńez and R. Baquero, Phys. Rev. B, 26 T. H. Fischer and J. Almlof, J. Phys. Chem. 96, 9768 86, 195106 (2012). (1992). 55 H. Jiang, J. Chem. Phys., 138, 134115 (2013). 27 H. J. Monkhorst and J. D. Pack, Phys. Rev. B, 13, 5188 56 H. Dixit, R. Saniz, S. Cottenier, D. Lamoen, and B. Par- (1976). toens, J. Phys.: Condens. Matter, 24, 205503 (2012). 28 F. Ortmann, F. Bechstedt, and W. G. Schmidt, Phys. Rev. 57 P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and B 73, 205101 (2006). J. Luitz, WIEN2K, an Augmented Plane Wave + Local Or- 29 A. Tkatchenko and M. Scheffler, Phys. Rev. Lett. 102, bitals Program for Calculating Crystal Properties, Techn. 073005 (2009) Universitat: Wien, Austria. ISBN: 3-9501031-1-1-2, 2001. 30 S. Grimme, J. Comp. Chem. 27, 1787 (2006). 58 R. Hoffman, Rev. Mod. Phys. 60, 601 (1988) 31 Z. Zheng, J. Xu and J. Zhao, High Press. Res. 30, 301 59 R. J. Colton, and J. W. Rabalais. J. Chem. Phys. 64, 3481 (2010). (1976). 32 D. C. Sorescu and B. M. Rice, J. Phys. Chem. C 114, 6734 60 P. Jain, J. Sahariya, H. S. Mund, M. Sharma and B. L. (2010). Ahuja, Comput. Mat. Sci. 72, 101 (2013) 33 S. Appalakondaiah, G. Vaitheeswaran, S. Lebégue, J. 61 See supplementary material for the figures 1, 2 and 3 cal- Chem. Phys, 138, 184705 (2013). culated using CASTEP code. 34 S. Appalakondaiah, G. Vaitheeswaran, S. Lebégue, J. 8

◦ TABLE I: Calculated (PWSCF and CASTEP) ground state lattice parameters (a, b, c, in A˚), crystallographic angle (α, in ), volume (V, in A˚3), and density (ρ, in gr/cc) of α- and β- polymorphic phases of AgCNO using diﬀerent dispersion corrected DFT functionals along with experimental data. α-AgCNO β-AgCNO Package Method a b c V ρ a α V ρ PWSCF PBE 4.701 10.765 5.869 297.00 3.35 9.571 115.99 442.94 3.37 DFT-D2 3.411 11.467 6.327 247.50 4.02 9.264 115.84 409.03 3.65 CASTEP PBE 4.773 10.743 5.847 299.80 3.32 9.550 116.06 436.64 3.42 DFT-D2 3.416 11.436 6.340 247.69 4.02 9.264 115.86 407.76 3.66 DFT-TS 3.883 11.009 5.872 251.05 3.97 9.223 115.76 407.04 3.67 DFT-OBS 4.352 10.749 5.826 272.54 3.65 9.379 115.90 421.24 3.54 Expt.a,b 3.880 10.752 5.804 242.13 4.11 9.087 115.73 390.58 3.82 a Ref.18, b Ref.19

TABLE II: Fifth degree polynomial ﬁts for pressure dependence of lattice parameters for α- and β- polymorphic phases of n AgCNO using PWSCF within DFT-D2 method. Each constant Cn(n = 0-5) has units of A˚/(GPa) . α-AgCNO β-AgCNO Parameter a b c a α C0 3.4114 11.467 6.3271 9.2677 115.84 −1 C1 (× 10 ) -0.693 -0.446 -1.161 -2.686 -2.263 −2 C2 (× 10 ) 2.085 -6.131 -3.011 -3.928 -9.700 −2 C3 (× 10 ) -0.697 3.032 1.509 4.810 7.062 −2 C4 (× 10 ) 0.109 -0.553 -0.242 -1.201 -1.645 −3 C5 (× 10 ) -0.065 0.324 0.141 0.970 1.294

FIG. 1: (Color online) Crystal structures of a) α-phase (2x3 super cell) along yz-plane, b) β-phase (2x2 super cell) along xy-plane and c) single molecule of β-phase. 9

-179050 DFT DFT-D2 -179150 -179100

-179200 -179150

2.5 GPa -179200 -179250 Enthalpy (KJ/mole/f.u) -179250 -179300

0 1 2 3 4 5 0 1 2 3 4 5 Pressure (GPa) Pressure (GPa)

FIG. 2: (Color online) Calculated enthalpy as a function of pressure for α- and β- polymorphic phases of AgCNO with (DFT-D2) and without (PBE-GGA) dispersion correction method

-179290

-179295

-179300

-179305

α-AgCNO -179310 β

Total Energy (KJ/mole/f.u) -AgCNO

-179315 0 1 2 3 4 5 Pressure (GPa)

FIG. 3: (Color online) Calculated total energy as a function of pressure for α- and β- polymorphic phases of AgCNO within DFT-D2 method. 10

-463700

α-AgCNO β-AgCNO

-463710

-463720

-463730 Gibbs free energy (KJ/mole/f.u)

-463740 0 100 200 300 400 Temperature (K)

FIG. 4: (Color online) Calculated Gibbs free energy as a function of temperature for α- and β- polymorphic phases of AgCNO within DFT-D2 method.

α-AgCNO β-AgCNO 65 /f.u) 3

60 Volume (Å 55

50 0 1 2 3 4 5 Pressure (GPa)

FIG. 5: (Color online) Calculated volume as a function of pressure for α- and β- polymorphic phases of AgCNO within DFT-D2 method. 11

1 (a) (b) 1

0.98 0.98

0.96 0.96 a α 0.94 0.94 a b 0.92 c 0.92

Normalized lattice parameters 0.9 0 1 2 3 4 5 0 1 2 3 4 5 Pressure (GPa) Pressure (GPa)

FIG. 6: (Color online) Calculated normalized lattice parameters of a) α- and b) β- polymorphic phases of AgCNO as a function of pressure within DFT-D2 method.

1.02 (a) (b)

0.99 0.99

0.96 Ag-C Ag-C Ag-C 0.96 Ag-O 0.93 Ag-O C-N C-N N-O N-O Ag-Ag Ag-Ag Normalized bond lengths

0 1 2 3 4 5 0 1 2 3 4 5 Pressure (GPa) Pressure (GPa)

FIG. 7: (Color online) Calculated normalized bond lengths of a) α- and b) β- polymorphic phases of AgCNO as a function of pressure within DFT-D2 method. 12

1.05 (a) (b) 1.02 1.02

0.99 0.99 Ag-C-Ag 0.96 Ag-C-N C-N-O C-Ag-C 0.96 0.93 Ag-C-Ag Ag-C-N Ag-C-N 0.93 0.9 C-N-O Normalized bond angles C-Ag-C

0 1 2 3 4 5 0 1 2 3 4 5 Pressure (GPa) Pressure (GPa)

FIG. 8: (Color online) Calculated normalized bond angles of a) α- and b) β- polymorphic phases of AgCNO as a function of pressure within DFT-D2 method.

12 12

9 9

6 6

3 3

0 0 Energy(eV) Energy(eV)

-3 -3

-6 -6

-9 -9

Γ Z T Y S X UR Γ TFL Γ

FIG. 9: (Color online) Calculated electronic band structures of α- (left) and β- (right) phases of AgCNO using TB-mBJ functional at the experimental crystal structures.18,19 13

9 30 Total 6 20 Total 3 10 0 0 0.9 0.9 O-s 0.6 0.6 O-s O-p O-p 0.3 0.3 0 0 0.9 0.9 N-s 0.6 0.6 N-s N-p N-p 0.3 0.3 0 0 0.2 C-s 0.2 C-s C-p C-p PDOS (States/eV) 0.1 PDOS (States/eV) 0.1 0 0 1.6 Ag-s 1.6 Ag-s Ag-p Ag-p 0.8 0.8 Ad-d Ag-d 0 0 -9 -6 -3 0 3 6 9 12 -9 -6 -3 0 3 6 9 12 Energy (eV) Energy (eV)

FIG. 10: (Color online) Calculated total and PDOS of α- (left) and β- (right) phases of AgCNO using TB-mBJ functional at the experimental crystal structures.18,19 14

FIG. 11: (Color online) Calculated electronic charge densities for α-phase (top) along (100), (010), (001) and for β- phase (bottom) along (100), (001) planes. 15

9 100 6 010 001 (ω)

1 3 ε

6 (ω) 2 ε 3

0 0 9 18 27 36 0 9 18 27 36 Photon energy (eV)

FIG. 12: (Color online) Calculated real (top) and imaginary (bottom) parts of dielectric function of α- (left) and β- (right) phases of AgCNO using TB-mBJ functional at the experimental crystal structures.18,19 16

150 ) -1 100 m 6 )(10 50 α(ω

0 8 100 010

) 001 6 -1

)(fs 4 ω ( σ 2

0 0 9 18 27 36 0 9 18 27 36 Photon energy (eV)

FIG. 13: (Color online) Calculated absorption (top) and photoconductivity (bottom) of α- (left) and β- (right) phases of AgCNO using TB-mBJ functional at the experimental crystal structures.18,19 arXiv:1406.3917v1 [cond-mat.mtrl-sci] 16 Jun 2014 uto ftogssC n N and CO gases two of duction n ufliae aesrne ntaigpwrthan power initiating stronger have fulminates Cu and tn rpte fvrosfliae eeeaie by examined were wohler fulminates and various Martin crackers of Christmas propeties in ating used is it ml uniyi sd nprtcnc,fiewrs and works, fire pyrotechnics, in pistols used) toy is quantity small sa olw:2AgCNO follows: as is fACOi ucett e aiu uptfrom output maximum get (RDX)) to 1,3,5-trinitro-1,3,5-triazine sufficient Hexahydro the is AgCNO of eh-iraiead3%mr ffiin o exploding for efficient more tri-nitro-xylene 30% and methl-nitroamine nteIainNavy Italian the in o ogtm ihgo eoaigproperties detonating good with time long a for saotfite ie seceta h ecr ful- mercury the as efficient as times (Hg(CNO) minate fifteen about is ie yahtwire hot a by nited neetv eoao n lofudt ea effective as is be which (Pb(N to azide found fulminates, lead also inorganic and as of detonator class effective an azides, the corresponding to the belongs cyanamides with explosives. and electronic secondary cyanates iso for are initiators They as applications find vrg eoainvlct fACOsml fthick- of sample AgCNO ness of velocity detonation average oyopimadtemdnmcgon tt fSle fulm Silver of state ground thermodynamic and Polymorphism nrai umntsaepiayepoie n they and explosives primary are fulminates Inorganic ∼ . mi nteodro 70ms 1700 of order the in is mm 0.5 6 h eopsto ehns fAgCNO of mechanism decomposition The . h rfre hroyai qiiru hs ne press under Hitherto observation. phase experimental equilibrium usin with calculations thermodynamic energy preferred m total the our (DFT+D) upon correction Based dispersion functionals. the semiempirical predicting by in role captured relat Densi significant a on on play based temperature interactions Waals approaches and pseudopotential pressure planewave of using effect the investigated rhrobc( orthorhombic ASnumbers: PACS AgCNO. ihnteTa-lh oie ek-ono oeta.Th wer potential. that Becke-Johnson potentia phases modified full Tran-Blaha polymorphic using the calculated within both were of properties and optical compressibility pressure and studied and the u state over seen method of is DFT+D transition with phase observed not polymorphic a Also respectively. ih fteplmrh r nlzdfo h acltdopt calculated a to the experimentalists from to indication analyzed early electric-d are an th anisotropy, provides polymorphs analyzed study Optical the is of bonding states. light chemical of of density nature partial The respectively. ivrfliae(gN)i rmr xlsv,wihexis which explosive, primary a is (AgCNO) fulminate Silver 3 .INTRODUCTION I. ti nw ob estv explosive sensitive a be to known is It . α 2 , 8–10 o xldn 2,4,6-tri-nitro-phenyl- exploding for ) 4 5 gN eoaoshsbe used been has detonators AgCNO . β tfid oeapiain (when applications some finds It . 3 ) n hyfudta h g Cd Ag, the that found they and , hssaeidrc adgpisltr ihabn a value gap band a with insulators gap band indirect are phases 2 → nispr om(x 0mg 60 (Ex. form pure its in ) dacdCnr fRsac nHg nryMtras(ACRHE Materials Energy High in Research of Centre Advanced A C N + 2CO + 2Ag 1 Cmcm ivrfliae(AgCNO) fulminate Silver . e al est ucinlcalculations functional density Waals der 2 ahbwi yeaa-5006 nhapaeh India pradesh, Andhra 046, 500 Hyderabad- Gachibowli, n rgnl( trigonal and ) iedsic os when noise distinct give 7 oevr h initi- the Moreover, . nvriyo yeaa,Po.C .RoRoad, Rao R. C. Prof. Hyderabad, of University .Ydknau n .Vaitheeswaran G. and Yedukondalu, N. 2 n h pro- the and ; − 1 hnig- when , Dtd uut2,2018) 27, August (Dated: R )frsa min odtos ntepeetsuy ehave we study, present the In conditions. ambient at forms 3) ¯ 2,3 Cmcm 2 the , It . and umnts esrdfiueo nestvns Fo )to I) Hg(CNO) of (F for insensitiveness the impact of to figure compared Measured insensitive fulminates. and stability thermal better a rtpitdotb Pandey by out AgCNO of pointed phases first ploymorphic two was of formation and try oycnrdotohmi tutr aigsaegroup space Ibam having structure orthorhombic centred body Singh Hg(CNO) 6.04 = c o en ieyue ssadr rmr explosive, Pb(N primary by replaced standard largely as was it used but widely being now tblt n oiiy h Pb(N The toxicity. and stability 3)adtmeaueo nnto To )frAgCNO for I) of (T ingnition of (170 temperature and (30) omuega oordet ogepsr olgtand light nature to hygroscopic exposure non long possess transform to crystals crystals these due AgCNO colour with white mouse-gray crystals the to fresh that like exception as years an properties 40 about similar water reported exhibit under been ago has stored it crystals However, AgCNO in]. that there references and 2 h w oyopi hssaerpre sorthorhombic of as ( structures reported are crystal phases polymorphic redetermined two the the Further, (1965). n twsrcgie yGay-Lussac by recognised was it and aaeesa=3.88 = a parameters n ,adtioa ( trigonal and 4, = Z and R Cmcm hssdntdas denoted phases 3 ¯ gN a rtpeae yBrugnatelli by prepared first was AgCNO o ierzdagetdpaewv method wave plane augmented linearized l pl rniin n ht estvt to sensitivity photo and transitions ipole hs tblt n hycnb accurately be can they and stability phase C) 15 F+ method, DFT+D g ihltiecntnsa=3.88 = a constants lattice with odtefraino unstable of formation the void igDTfntoasadtesm was same the and functionals DFT sing hs yBarrick by phase ) < eotdtecytlsrcueo gN obe to AgCNO of structure crystal the reported acltdeetoi tutr shows structure electronic calculated e yFntoa hoy(F) a der van (DFT). polymorphs Theory the Functional of ty stability phase ive r n eprtr en consistent being temperature and ure si w oyopi hsssc as such phases polymorphic two in ts 2 A ˚ netgtd lcrncstructure Electronic investigated. e og h hredniyposand plots density charge the rough tosicnrs osadr DFT standard to incontrast ethods Hg(CNO) oee,Hg(CNO) However, . clseta vrl,tepresent the Overall, spectra. ical n .Teacat rsa symme- crystal accuarte The 4. = Z and , eprtr ag.Teequation The range. temperature ∗ α 2 2 A ˚ and (210 f35 n .3eV, 4.43 and 3.51 of s 10.752 = b , (10) Cmcm R nt tde rmvan from studied inate β )paeb Britton by phase 3) ¯ o frso AgCNO, of -forms M), C) tal et < < sfudt be to found is gN (22) AgCNO 18 2 Pb(N 3 16 a h rtwihis which first the was ) 17) ihunit-cell with (1979), β 2 3 fr of -form A ˚ n Britton and ) A ˚ 14 a on ohave to found was 2 3 11.20 = b , 11 n 5.804 = c and ) 12) n1959, In (1824). 2 u oispoor its to due ( > 300 19 < 2,12 13 19) a (1991), Pb(N o . )[Ref C) (1798) tal et A ˚ and 3 A ˚ 17 ) 2 , 2

= 9.087A˚, α = 115.73◦ and Z = 6. The differential ther- coefficient for the pair of ith and jth atoms that depends 20 mal analysis studies indicate no evidence for temper- only on the chemical species, and Rij is an interatomic o o 1 ature induced phase changes between 20 C and 145 C. distance. fdamp = − − is a damping function 1+e d(Rij /R0 1) To the best of our knowledge, there are no studies to which is necessary to avoid divergence for small values explore the relative stability under pressure and temper- of Rij and R0 is the sum of atomic vdW radii. The ature, electronic structure, and chemical bonding of the DFT+D method has been applied to several molecular two polymorphic phases in the literature from both the- systems such as secondary explosives30–33 and inorganic oretical and experimental point of view for the energetic azides34–36, and these results demostrates that disper- solid AgCNO. Therefore, we systematically investigated sion corrections are important in determining the ground the above mentioned properties for AgCNO using first state properties. Hence, we attempted to study the phase principles calculations based on Density Functional The- stability of AgCNO using the DFT+D method and the ory (DFT). The rest of the article is organized as follows, results are discussed in the following section. in section II, we briefly describe the methodology of our The thermodynamic stabilities were investigated by calculations. In section III, the polymorphic phase sta- comparing the free energy of the polymorphs. In order to bility, electronic struture, chemical bonding, and optical predict the most stable form of AgCNO, we have calcu- properties of AgCNO are discussed. Finally, in section lated the Gibbs free energy of both polymorphic phases. IV, we summarize the results, which concludes our paper. The Gibbs free energy for any material at given temperature is given by

II. METHODOLOGY OF CALCULATIONS G(P,T )= F (V,T )+ P V = Fvib + Fperfect + P V (3)

The polymorphic phases of AgCNO crystals were stud- Where F, P, T, and V are the Helmoltz free energy, ied using periodic DFT calculations performed with two pressure, absolute temperature and volume respectively. different codes such as Plane Wave Self Consistent Field F(V,T) is summation of vibrational free energy and per- (PWSCF)21 and Cambridge series of total energy pack- fect lattice energy i.e F(V, T) = Fvib + Fperfect. The age (CASTEP)22, which are based on planewave pseu- vibrational free energy Fvib is calculated within the har- dopotential (PW-PP) approach. We have used Van- monic approximation and is given by 23 − derbilt type ultrasoft pseudo potentials to treat ion- hωi hωi − K T electron interactions, while electron-electron interactions Fvib = KBT + ln(1 e B ) (4) X 2KBT are treated with the local density approximation (LDA) i developed by Ceperley and Alder24 parameterized by Where ωi is the phonon frequency, T is the absolute Perdew and Zunger25 and generalized gradient approx- temperature, h is the Planck constant and KB is the imation (GGA) by Perdew-Burke-Ernzerhof (PBE)26. Boltzmann constant. Where Fperfect = E0 + Eel - TSel, The Broyden-Fletcher-Goldfarb-Shanno (BFGS) mini- 27 E0 and Eel- TSel are the contributions from the lattice mization scheme was used in geometry optimiza- and electronic excitations, respectively. If the electronic tion. The plane wave cutoff energy was set to 60 Ry excitations and pressure changes are neglected then the and a k-point grid of 7×7×4 and 5×5×5 for Cmcm Gibbs free energy can be calculated by G = E0 + Fvib. and R3¯ phases, respectively using Monkhorst-Pack grid scheme28. The self-consistent energy convergence ac- −6 curacy was set to 1.0×10 eV/atom and the force to III. RESULTS AND DISCUSSION 1.0×10−4 eV/A. The standard LDA/GGA functionals are unabale to A. Crystal structure and phase stability of describe the non-covalent interactions in molecular crys- polymorphs tals. On the other hand, semiemipirical dispersion correction method (DFT+D) has been proposed by Extreme sensitivity, poor stability and high cost of Sil- Grimme29 implemented through standard PBE-GGA ver as a raw material prohibited AgCNO in commercial functional to treat van der Waals (vdW) interactions and or military priming and detonating devices37,38 but still it are found to improve the results to a greater extent. Here finds applications in pyrotechincs, fire works, toy pistals the total energy is given by and in the navy. However, the relative stability and nec- EDF T +D = EDF T + Edisp (1) essary conditions for the formation of two polymorphic phases are still unknown for this material. As discussed Where E is the self-consistent Kohn-Sham energy, DF T in section I, AgCNO exists in two polymorphic phases E is empirical dispersion correction given by disp at ambient conditions, and it is noticed that except the C crystal structures17–19 of both phases most of the physi- E = −S ij f (R ) (2) disp 6 6 damp ij cal properties are not well understood from experimental X Rij i

tivity, complex chemical behavior and risk of decompo- B. Pressure and temperature eﬀects on the sition. Therefore, theoretical modeling and simulations polymorphs are eﬃcient tools to predict the physical and chemical properties of complex energetic solids at extreme condi- The investigation of energetic materials at extreme tions and also very useful in predicting results a priori conditions is a challenging task because of their sensi- to the experimentalists. DFT is a powerful tool in pre- 4

-179050 reveals that there is no ploymorphic phase transition ob- DFT DFT-D2 -179150served in AgCNO and the Cmcm phase is found to be the -179100 most stable form of AgCNO over studied pressure range using both PW-PP approaches (see Fig. 2). Hitherto, ¯ -179200we represent Cmcm and R3 phases as α and β-AgCNO, 45 -179150 respectively. Our recent high pressure study on KClO3 also reveals that PBE+D method is good enough for cal- 2.5 GPa -179250culating the transition pressures for molecular solids in -179200 opposition to standard LDA/GGA functionals. As shown in Fig. 3 (see Fig. 2 of supplementary material59), the Enthalpy (KJ/mole/f.u) -179250 -179300calculated total energies ploted as a function of pressure demonstrates that the PBE+D functional predict the α- 0 1 2 3 4 5 phase0 to be the1 most2 stable form 3 of AgCNO4 over 5 studied Pressure (GPa) pressure range inPressure contrast to (GPa) usual LDA/GGA functionals. FIG. 2: (Color online) Calculated enthalpy as a function of In addition, the contribution of lattice vibration needs pressure for α and β polymorphic phases of AgCNO. to be considered at elevated temperatures. We have used norm-conserving pseudo potentials46 for calculating the -179290 Gibbs free energy as they are well suited for phonon calculations as implemented in CASTEP code. Fig. 4 shows -179295 the Gibbs free energy as a function of temperature in the 0-450 K range within PBE+D method and it illus-

-179300 trates that the difference in Gibbs free energy (∆G) decreases monotonically with increase in temperature between the two α- and β-polymorphs of AgCNO. How- -179305 ever, α-AgCNO remains the form with the lowest free energy throughout the investigated temperature range, α-AgCNO -179310 despite the contribution of lattice vibration tending to β Total Energy (KJ/mole/f.u) decrease the ∆G gap at elevated-AgCNO temperatures. We have not observed any temperature driven phase transition be- -179315 0 1 2tween two 3 polymorphs over4 studied temperature 5 range Pressureand this(GPa) is consistent with results of differential thermal analysis20. Gibbs free energies at temperatures beyond FIG. 3: (Color online) Calculated total energy as a function 450 K are not plotted because α-AgCNO phase typically of pressure for α and β polymorphic phases of AgCNO. decompose at higher temperatures between 393-450 K20. From our present study, we confirmed that α-phase is the most stable polymorph of AgCNO under studied pressure dicting the behavior of complex solid state systems at ex- and temperature range. treme conditions. Hence, we attempted to study the relative phase stability and possible polymorphic structural transition in AgCNO under hydrostatic pressure upto 50 Kbar with a step size of 5 Kbar. As illustrated in Fig. C. Equation of state and compressibilities of the 2, The calculated enthalpy curves as a function of pres- polymorphs sure using standard DFT functionals show that R3¯ phase undergoes a structural ploymorphic phase transition to To study the effect of hydrostatic pressure on crystal Cmcm phase at about 21.85 Kbar and 24.56 Kbar respec- structures of both the polymorphs, we have used vari- tively, using LDA and GGA with PWSCF and the cor- able cell optimization technique. The calculated vol- responding transition pressures obtained from CASTEP umes decrease monotonically as a function of pressure code are 23.34 and 27.17 Kbar (see Fig. 1 of supplemen- using standard LDA and GGA functionals as well as tary material59). Overall, we observed a similar trends with PBE+D method are shown in Fig. 5 (see Fig. 3 using both of the approaches with small deviations in of supplementary material59) and the obtained equilib- the calculated transition pressures. Above the transition rium bulk moduli for α and β-AgCNO are tabulated in pressures the enthalpy of the Cmcm phase is lower than Table I by fitting pressure-volume data to second order that of the R3¯ phase, indicating that Cmcm structure be- Birch-Murnaghan equation of state.47 Further, to inves- comes stable. The Cmcm phase is stable in the pressure tigate the relative compressibilities of both the polymor- range from transition pressure to 50 Kbar, which is the phic phases, we have calculated the lattice parameters, highest pressure we accomplished in the present work. In bondlengths and bond angles as a function of hydrostatic contrast to LDA/GGA functionals, the PBE+D method pressure using PBE+D method within PWSCF. As illus- using both of the approaches (PWSCF and CASTEP) trated in Fig. 6, the lattice is most compressible along 5

-463700 Ag-O are more compressible under the application of hy- α-AgCNOdrostatic pressure. Further, this compressibility behav- β-AgCNOior can be clearly understand by analyzing the nature of chemical bonding in AgCNO and hence we have inves- -463710 tigated the electronic struture of both the polymorphic phases.

-463720 1 (a)

0.98 -463730 Gibbs free energy (KJ/mole/f.u) 0.96

-463740 0 100 200 300 400 Temperature (K) 0.94 a b FIG. 4: (Color online) Calculated Gibbs free energy as a c function of temperature for α and β polymorphic phases of 0.92 AgCNO. Normalized lattice parameters

70 0 1 2 3 4 5 0 1 α-AgCNO Pressureβ-AgCNO (GPa) Pressure (GPa) 65 FIG. 6: (Color online) Calculated normalized lattice param- /f.u) 3 eters of a) α and b) β polymorphic phases of AgCNO as a 60 function of pressure Volume (Å 55 1.02 (a) (b) 50 0 1 2 3 4 5 Pressure (GPa) 0.99 FIG. 5: (Color online) Calculated volume as a function of pressure for α and β polymorphic phases of AgCNO.

0.96 c-axis and least compressible along b-axis for α-AgCNO Ag-C whereas lattice constant (a) and rhombohedral angle (α) Ag-C Ag-C are decreasing with pressure for β-AgCNO. As shown Ag-O Ag-O in Fig. 7, Ag-O bond is more compressible in both the 0.93 C-N phases along the a-crystallographic direction whereas Ag- C-N N-O N-O Ag is the most compressible in β-AgCNO. The bonds Ag-Ag Ag-Ag such as Ag-C, C-N, and N-O show similar and less com- Normalized bond lengths pressible over studied pressure range due to strong covalent bonding between Ag and C as well as within the 0 1 2 3 4 5 0 1 CNO anion molecule and this can be clearly understood Pressure (GPa) Pressure (GPa) from their electronic structure. Similarly the corresponding bond angles Ag-C-Ag decrease because of the orien- FIG. 7: (Color online) Calculated normalized bond lengths of tation of Ag-C bonds due to high compressibility nature a) α and b) β polymorphic phases of AgCNO as a function of Ag-Ag bonds (see Fig. 7) whereas Ag-C-N bond angle of pressure increases under compression. Overall, we observe that Ag-C, C-N, and N-O bonds are stiﬀer, while Ag-Ag and 6

1.05 12 (a) (b)

1.02 9 1.02

0.99 6 0.99 Ag-C-Ag Ag-C-N 0.96 3 C-N-O C-Ag-C 0.96 0.93 Ag-C-Ag Ag-C-N 0

Ag-C-NEnergy(eV) 0.93 0.9 C-N-O Normalized bond angles C-Ag-C -3

0 1 2 3 4 5 0 1 -62 3 4 5 Pressure (GPa) Pressure (GPa)

-9 FIG. 8: (Color online) Calculated normalized bond angles of a) α and b) β polymorphic phases of AgCNO as a function Γ Z T Y S X UR of pressure 12

D. Electronic structure and chemical bonding 9

Inorganic fulminates are isoelectronic with azides, cy- 6 nates, and cynamides. Iqbal and his co-workers1 explained these fulminates to be the most sensitive explosives to shock and heat among the energetic materials 3 such as azides, cynates, thiocynates, and cynamides due to asymmetric charge distribution in their structures. Also, the heavy metal complex salts are unstable than 0 alkali metal salts because of the asymmetric inter ionic Energy(eV) distances44. In order to understand the macroscopic en- -3 ergetic behavior of these sensitive explosive materials, it is necessary to understand the electronic structure, and 1 chemical bonding of these materials. Iqbal et al car- -6 ried out a detailed experimental study on the electronic structure and stability of the inorganic fulminates, which reveals that sodium, potasium, and thallous fulminates -9 to be ionic salts whereas silver and mercury salts are covalent in nature and this will be reflected in the order of Γ TFL Γ stability of the fulminate salts. However, the electronic structure and bonding properties of the polymorphs of FIG. 9: (Color online) Calculated electronic band structures silver fulminate salts are not reported in the literature. of α (left) and β (right) phases of AgCNO using TB-mBJ Therefore in the present study we have made a detailed functional at the experimental crystal structure analysis of electronic structure and chemical bonding of both the polymorphic phases at ambient pressure. Blaha modified Becke Johnson (TB-mBJ) potential48 is It is well known fact that the standard LDA/GGA computatinally less expensive method, which provides functionals severely underestimate the band gaps for accurate energy band gaps as comparable with more so- semiconductors and insulators. In order to get rela- phisticated methods within the Kohn-Sham frame work. iable energy band gaps several methods have been pro- Recently, Koller et al 49,50 reported the merits and limi- posed such as LDA+U, LDA+DMFT, hybrid fuctionals, tations of the TB-mBJ functional. So far, several groups and GW approximation. Unfortunately these methods used the TB-mBJ potential for the calculation of elec- are computationally very expensive. However, Tran and tronic structure and optical properties of different kinds 7 of materials51–53 and confirmed that bandgaps obtained covalent character. Further, this can be clearly under- using TB-mBJ potential are improved when compared stood from electronic charge density plots which are used to LDA/GGA for a wide range of materials (see Fig. 1 for accurate description of chemical bonds.56 The calcu- from Ref. 54). The accurate prediction of band gaps us- lated valence charge density using TB-mBJ functional for ing TB-mBJ functional for diverse materials motivated us both the polymorphs are as shown in Fig. 11. It shows to use this functional for high energetic materials. This anisotropic bonding interactions and the charge cloud is semi local functional is implemented through WIEN2k distributed within the CNO molecule indicating covalent package55. Since all these energetic materials are semi- character. Overall, the C, N, and O atoms are covalenlty conductors or insulators, we have used the TB-mBJ func- bonded within CNO anion and the metal atom is also co- tional to calculate accurate energy band gap values for valently bonded with CNO group through C atom. X-ray the two ploymorphs. The calculated band gaps are found electron spectroscopy57 study on inorganic azides reveals to be 3.51 and 4.43 eV for α and β-phases of AgCNO, re- that Heavy Metal azides (HMAs) are more covalent than spectively and the respective band gaps using LDA func- Alkali metal azides (AMAs), implying that HMAs are tional are 2.00, 2.85 eV. The obtained TB-mBJ band more sensitive than AMAs and find applications as in- gap of 3.51 eV for α-AgCNO is closely comparable with tiators for secondary explosives. Therefore, the presence the measured optical energy gap of 4.0 eV1 over LDA of covalent bonding in AgCNO makes it more sensitive functional. Also, the calculated TB-mBJ band structures than the above mentioned ionic inorganic fulminates and show that α and β-AgCNO phases are indirect band gap it can be used as an initiator as HMAs. insulators along S-(Γ-Z) and T-(F-Γ) high symmetry directions of the Brillouin zone as shown in Fig. 9. The nature of chemical bonding in the two polymorphic phases was investigated by examining the total and E. Optical properties partial density of states (DOS) as illustrated in Fig. 10. The relative phase stability of the α and β polymorphs of Energetic materials become unstable and they undergo AgCNO can also be explained on the basis of the DOS. photochemical decomposition by the action of light. As shown in Fig. 10, when compared to β-AgCNO the Hence it is interesting to study the optical properties of peaks of the DOS in the valence bands of α-AgCNO has these materials to understand the decomposition mecha- a tendency to shift towards lower energy, indicating bet- nisms. Electronic structure calculations could provide an ter stability of the α polymorphic phase. Also, it can be information about the nature and location of interband seen that total DOS of both phases exhibit some simi- transitions in crystals. The complex dielectric function lar features, hence we have analyzed the bonding based ǫ(ω) = ǫ1(ω) + iǫ2(ω) can be used to describe the lin- on PDOS of α-AgCNO. It is clearly visualized from the ear response of the system to electromagnetic radiation DOS that AgCNO polymorphs have molecular character, which is related to the interaction of photons with elec- which arises from a strong overlap contributions of partial trons. The imaginary part of dielectric function ǫ2(ω) is DOS of Silver and Fulminate anions in the valence and obtained from the momentum matrix elements between conduction bands. The conduction band is mainly due to the occupied and unoccupied wavefunctions within selec- p-states of C, N, O and s,p,d-states of metal (Ag) atoms. tion rules. The optical properties of AgCNO polymorphs The lowest lying states about -9 eV are due to hybridized were calculated using TB-mBJ functional at a denser k- C-s and s, p states of N and O atoms. The states from mesh of 20 × 20 × 12 and 17 × 17 × 17 for α- and -5 to -7 eV are due to hybridization of Ag-d and p-sates β- phases, respectively. The AgCNO polymorphs crys- of C, N and O atoms. Also, the top of the valence band tallize in the orthorhombic and hexgonal symmetry, and is mainly dominated by O-p and Ag-d states and there this allow non zero components of the dielectric tensors is a charge sharing between C-p and Ag-d states repre- three for α and two for β phase along [100], [010] and senting the covalent nature of Ag-C bond. The nature of [001] directions. Fig. 12 shows the real ǫ1(ω) and imag- chemical bond between two atoms can be predicted from inary ǫ2(ω) part of the dielectric function as a function their electronegativity difference. When the differences of photon energy for both polymorphs. The peaks in are greater or equal to 1.7 (ionic) and less than 1.7 (co- ǫ2(ω) mainly arises due to electric-dipole transitions be- valent), if the difference is greater than 0.5 the covalent tween valence and conduction bands. The major peaks in bond has some degree of polarity. According to Paul- ǫ2(ω) at 4.57 eV (along 100), 5.94 eV (along 010), 5.06 eV ing scale the electronegetivity values Ag (1.9), C (2.5), (along 001) in α-phase and 5.45 eV (along 100), 6.13 eV N (3.0), and O (3.4) shows that there is a covalent bond (along 001) in β-phase arise due to interband transitions between N-O, C-N and Ag-C bonds. The electronege- from Ag(4d) → N(2p)-states. It can be clearly seen from tivity difference in N-O (0.4), C-N (0.5) and Ag-C (0.6) DOS, the top of valnce band is mainly dominated by 4d- exhibits strong and polar covalent nature of N-O, C-N states of Silver atom (see Fig. 10) in both polymorphs, 58 and Ag-C bonds, respectively. The N-O, C-N and Ag-C which is similar to the case of Silver azide (AgN3) . As bonds show less compressibility with increasing pressure discussed in the section I, AgCNO is iso-electronic and it (see Fig. 7), this is due to strong hybridization between consists of same chemical species as AgN3 in addition to Ag-d and s, p-states of C, N, and O atoms leads to strong C and O atoms, hence one can expect that major opti- 8

9 Total 6 3 0 0.9 0.6 O-s O-p 0.3 0 0.9 0.6 N-s N-p 0.3 0 0.2 C-s C-p PDOS (States/eV) 0.1 0 1.6 Ag-s Ag-p 0.8 Ad-d 0 -9 -6 -3 0 3 6 9 12 Energy (eV) 30 20 FIG. 11: (Color online) Calculated electronic charge densities Total of a) α and b) β phases along (100) and (001), (110) planes, 10 respectively. 0 0.9 0.6 O-s O-p 0.3 0 1.50, n010 = 2.24, n001 = 1.67 and n100 = 1.61, n001 = 0.9 2.09 for α and β-phases of AgCNO, respectively. The 0.6 N-s N-p obtained refractive indices are distinct in all the crys- 0.3 tallographic directions, which indicates the anisotropy of 0 the α and β polymorphs. Using ǫ1, ǫ2 one can dervive the 0.2 C-s optical constants such as absorption and photoconductiv- C-p ity of the materials. The calculated absorption spectra is PDOS (States/eV) 0.1 shown in Fig. 13 and absorption starts after the energy 0 3.51 eV in α-phase, 4.43 eV in β-phase, which is energy 1.6 Ag-s band gap between the valence band maximum and con- Ag-p 0.8 duction band minimum. The first absorption peaks along Ag-d three (100, 010 and 001) directions are at 4.5 eV, 7.7 eV, 0 -9 -6 -3 0 3 6 9 12 7.9 eV and the corresponding absorption coefficients are − − − Energy (eV) found to be 1.1 × 107m 1, 1.4 × 108m 1, 3.9 × 108m 1 for α-phase, while the same are 1.3 × 108m−1 at 6.8 eV and 4.0 × 107m−1 at 7.2 eV along two (100 and 001) FIG. 10: (Color online) Calculated total and PDOS of α (left) and β (right) phases of AgCNO using TB-mBJ functional at directions for β-phase. Photoconductivity is due to the the experimental crystal structure. increase in the number of free carriers when photons are absorbed. The calculated photoconductivity shows (see Fig. 13) a wide photocurrent response in the absorption region of 3.5-35 eV and 4.4-35 eV in α and β phases of → cal transitions arise from Ag(4d) N(2p)-states in both AgCNO, respectively. overall, the energetic polymorphs 58 energetic materials AgCNO and AgN3 . show a strong anisotropic and wide range of absorption 7 −1 The real part ǫ1(ω) of can be derived from ǫ2(ω) using (absorption coefficients ∼ 10 m ). This results indicate the Kramer-Kronig relations. The calculated real static the possible decomposition of AgCNO into Ag, CO and dielectric constant along all the crystallographic direc- N2 under the action of ultraviolet (UV) light. Therefore tions are 2.24, 5.01, 2.78 for α-phase and 2.59, 4.38 for AgCNO polymorphs decompose under the action of UV β-phase. The static refractive indices calculated using light and they explode due to photochemical decomposi- the dielectric function n = ǫ(0) are given by n = tion. p 100 9

9 150 100 )

6 -1 010 100 m 001 6 (ω)

1 3 )(10 ε 50 0 α(ω

0 8 9

) 6 -1 6 (ω) )(fs 2 4 ω ε (

3 σ 2

0 0 0 9 18 27 36 0 90 18 9 27 18 36 27 36 09 Photon energy (eV) Photon energy (eV)

FIG. 12: (Color online) Calculated real (top) and imaginary FIG. 13: (Color online) Calculated absorption (top) and pho- (bottom) parts of dielectric function of α (left) and β (right) toconductivity (bottom) of α (left) and β (right) phases of phases of AgCNO using TB-mBJ functional at the experi- AgCNO using TB-mBJ functional at the experimental crys- mental crystal structure tal structure

IV. CONCLUSIONS

In summary, ab-initio calculations have been performed to investigate the relative phase stability, structural transition and compressbilities of both polymor- cord with each other for both polymorphs. We also have phic orthorhombic (Cmcm) and trigonal (R3)¯ phases calculated the electronic structure and density of states of AgCNO using PW-PP approachs (PWSCF and both polymorphic phases using TB-mBJ functional at CASTEP). The standard LDA/GGA functionals are in- ambient pressure. The calculated electronic band struc- adequate to predict the relative phase stability and ture show that α, β phases are indirect band gap insula- ground state properties of these polymorpic phases. tors with a band gap values of 3.51 and 4.43 eV, respec- From our total energy calculations using DFT+D tively. The detailed analysis of electronic charge density method, Cmcm phase was found to be preferred thermo- distribution and partial density of states reveal covalent dynamic equilibrium phase under pressure and tempera- bonding is predominant in the energetic AgCNO poly- ture range and hence we conﬁrmed Cmcm and R3¯ phases morphs. The major peaks in ǫ2(ω) are due to interband as α and β-AgCNO, respectively. The present study also transitions between Ag(4d) → N(2p)-states. The calcu- reveals that there is no structural phase transition ob- lated absorption spectra reveal that the polymorphs show served using DFT+D method whereas we observe a pres- considerable anisotropy and decompose under the action sure induced polymorphic phase transition from β → α at of UV light. Finally, the present study suggests that the about 21.85 (23.34) and 24.56 (27.17) Kbar within LDA α-form is stable and it can be used in diﬀerent applica- and GGA using PWSCF (CASTEP), respectively. Over- tions (when small quantity is used) and β-form can be all, the phase stability and ground state properties ob- avoided by careful control of the pressure or temperature tained from both PWSCF and CASTEP are in good ac- in the manufacturing process. 10

V. ACKNOWLEDGMENTS

Authors would like to thank Defence Research and Development Organisation (DRDO) through ACRHEM for the ﬁnancial support under grant No. DRDO/02/0201/2011/00060:ACREHM-PHASE-II, and the CMSD, University of Hyderabad, for providing computational facilities. Its pleasure to acknowledge Dr. V. Kanchana, Department of Physics, Indian Institute of Technology Hyderabad for critical reading of the manuscript. ∗Author for Correspondence, E-mail: [email protected] 11

1 Z. Iqbal and A. D. Yoffe, Proc. R. Soc. Lond. A, 302, 35 33 S. Appalakondaiah, G. Vaitheeswaran, S. Lebégue, J. (1967). Chem. Phys 140, 014105 (2014). 2 P. H. Collins, K. J. Holloway, Propellants and explosives, 34 K. Ramesh Babu, G. Vaitheeswaran, Chem. Phys. Lett., 3, 159 (1978). 533, 35 (2012). 3 T. L. Davis, The chemistry of powder and explosives (Wi- 35 K. Ramesh Babu, G. Vaitheeswaran, J. Chem. Sci. 124, ley, New York 1943). 1391 (2012). 4 F. P. Bowden and H. T. Williams, Proc. R. Soc. Lond. A, 36 K. Ramesh Babu, G. Vaitheeswaran, Chem. Phys. Lett. 208, 176 (1951). 586, 44 (2013). 5 T. Urbanski, Chemistry and Technology of Explosives, 3, 37 G. W. C. Taylor, and J. M. Jenkins, Proc. 3rd Simposium 157 (1967). on Chemical problems connected with the stability of Ex- 6 J. P. Springer, Ger Pat 953, 415 (1956). plosives (Ystad, Sweden, May 28-30, 33, 1973). 7 D. A. Jackson and A. P. Dicks, J. Chem. Educ., 89, 1267 38 C. A. Taylor and E. P. Buxton, Army Ordnance, 6, 118 (2012). (1925). 8 F. Martin, Dber Azide und Fulminate Darmstadt (1913). 39 E. F. C. Byrd, and B. M. Rice, J. Phys. Chem. C, 111, 9 L. Wohler and F. Martin, Angezu. Chem. 30, 33 (1917). 2787 (2007). 10 L. Wohler and F. Martin, Ber. Dtsch. Chem. Ges. 50, 586 40 Z. Wu, R. K. Kalia, A. Nakano and P. Vashishta, J. Chem. (1917). Phys., 134, 204509 (2011). 11 F. D. Miles, A history of research in the Nobel Division of 41 J. J. Zhao, and H. Liu, Comput. Mater. Sci. 42, 698 (2008). ICI (1955). 42 H. Liu, J. J. Zhao, J. G. Du, Z. Z. Gong, G. F. Ji, and D. 12 A. M. Peter, J. Am. Chem. Soc., 38, 486 (1916). Q. Wei, Phys. Lett. A, 367, 383 (2007). 13 Encyclopedia of Explosives and Related Items (PATR 2700, 43 L. Qiu, H. M. Xiao, W. H. Zhu, J. J. Xiao, and W. Zhu, Volume 6, Dover, New Jersey, USA, Picatinny Arsenal J. Phys. Chem. B, 110, 10651 (2006). 1974). 44 Z. Iqbal, Structure and Bonding, Springer, 10, 25 (1972). 14 J. L. Gay -Lussac, Ann. Chim. Phys. 27, 199 (1824). 45 N. Yedukondalu, Vikas D. Ghule and G. Vaitheeswaran, J. 15 K. Singh, Acta Cryst, 12, 1053 (1959). Chem. Phys., 138, 174701 (2013). 16 S. N. Pandey, Indian J. Phys. 36, 657 (1962). 46 N. Troullier and J. L. Martins, Phys. Rev. B, 43, 1993 17 D. Britton, and J. D. Dunitz, Acta. Cryst. 19, 662 (1965). (1991). 18 J. C. Barrick, D. Canfield, and B. C. Giessen, Acta. Cryst. 47 F. D. Murnaghan, Proc. Natl. Acad. Sci. USA 30, 244 B35, 464 (1979). (1944). 19 D. Britton, Acta Cryst. C47, 2646 (1991). 48 F. Tran, and P. Blaha, Phys. Rev. Lett. 102, 226401 20 T. Boddington and Z. Iqbal, Trans. Faraday Soc., 65, 509 (2009). (1969). 49 D. Koller, F. Tran and P. Blaha, Phys. Rev. B 85, 155109 21 P. Giannozzi et al, J. Phys.: Condens. Matter, 21 395502 (2012). (2009). 50 D. Koller, F. Tran and P. Blaha, Phys. Rev. B 83, 195134 22 M. C. Payne, M. P. Teter, D. C. Allen, T. A. Arias and J. (2011). D. Joannopoulos, Rev. Mod. Phys. 64, 1045 (1992). 51 D. J. Singh, Phys. Rev. B, 82, 205102 (2010). 23 D. Vanderbilt, Phys. Rev. B, 41, 7892 (1990). 52 J. A. Camargo-Martńez and R. Baquero, Phys. Rev. B, 24 D. M. Ceperley and B. J. Alder, Phys. Rev. Lett. 45, 566 86, 195106 (2012). (1980). 53 H. Jiang, J. Chem. Phys., 138, 134115 (2013). 25 J. P. Perdew and A. Zunger, Phys. Rev. B, 23, 5048 (1981). 54 H. Dixit, R. Saniz, S. Cottenier, D. Lamoen, and B. Par- 26 J. P. Perdew, S. Burke and M. Ernzerhof, Phys. Rev. Lett. toens, J. Phys.: Condens. Matter, 24, 205503 (2012). 77, 3865 (1996). 55 P. Blaha, K. Schwarz, G. K. H. Madsen, D. Kvasnicka, and 27 T. H. Fischer and J. Almlof, J. Phys. Chem. 96, 9768 J. Luitz, WIEN2K, an Augmented Plane Wave + Local Or- (1992). bitals Program for Calculating Crystal Properties, Techn. 28 H. J. Monkhorst and J. D. Pack, Phys. Rev. B, 13, 5188 Universitat: Wien, Austria. ISBN: 3-9501031-1-1-2, 2001. (1976). 56 R. Hoffman, Rev. Mod. Phys. 60, 601 (1988) 29 S. Grimme, J. Comp. Chem. 27, 1787 (2006). 57 R. J. Colton, and J. W. Rabalais. J. Chem. Phys. 64, 3481 30 Z. Zheng, J. Xu and J. Zhao, High Press. Res. 30, 301 (1976). (2010). 58 P. Jain, J. Sahariya, H. S. Mund, M. Sharma and B. L. 31 D. C. Sorescu and B. M. Rice, J. Phys. Chem. C 114, 6734 Ahuja, Comput. Mat. Sci. 72, 101 (2013) (2010). 59 See supplementary material for the figures 1, 2 and 3 cal- 32 S. Appalakondaiah, G. Vaitheeswaran, S. Lebégue, J. culated using CASTEP code. Chem. Phys, 138, 184705 (2013).