Structural and vibrational properties of solid nitromethane under high pressure by density functional theory

Analysis of the initial reaction mechanism of TKX-50 based on Raman intensity

CONTEXT

Dihydroxylammonium 5,5'-biotetrazolium-1,1'-diolate (TKX-50) has two important properties of typical azole energy-containing ionic salts, including high energy and safety. Therefore, in today's era where more emphasis is placed on explosive performance and explosive detonation control conditions, TKX-50 is a very important object of research, and its reaction process in the initial stage of detonation is gradually receiving more and more attention from researchers in the field of energy-containing materials research.

METHODS

In this paper, based on first-principles density-functional theory (DFT), the mechanism of chemical bond breakage of TKX-50 under pressure was determined based on the analysis of the strength and stability of chemical bonds inside the TKX-50 molecules using Raman spectroscopy relative intensity analysis. The results show that TKX-50 is dominated by N-H bond breaking and followed by H-O bond breaking in the initial reaction stage. These reactions lead to the reorganization and structural changes within the molecule, which eventually lead to the decomposition of TKX-50.

Triggering the mechanism of the initial reaction of energetic materials under pressure based on Raman intensity analysis

The Raman intensity and other stoichiometric calculations of nitromethane (NM) and 2-nitrimino-5-nitro-hexahydro-1,3,5-triazine (NNHT) have been made by using first-principles density functional theory. We propose a method to judge the initial reaction mechanism of NM and NNHT under pressure based on the Raman intensity. Both the resulting NM and NNHT undergo hydrogen transfer and conventional trigger bond cleavage. And the results obtained from the Raman peak intensities infer a reaction path that is not inferior to the traditional C-NO₂ and N-NO₂ bond cleavage, thus verifying our results.

Performance of Dispersion-Inclusive Density Functional Theory Methods for Energetic Materials

Molecular crystals of energetic materials (EMs) are denser than typical molecular crystals and are characterized by distinct intermolecular interactions between nitrogen-containing moieties. To assess the performance of dispersion-inclusive density functional theory (DFT) methods, we have compiled a data set of experimental sublimation enthalpies of 31 energetic materials. We evaluate the performance of three methods: the semilocal Perdew-Burke-Ernzerhof (PBE) functional coupled with the pairwise Tkatchenko-Scheffler (TS) dispersion correction, PBE with the many-body dispersion (MBD) method, and the PBE-based hybrid functional (PBE0) with MBD. Zero-point energy contributions and thermal effects are described using the quasi-harmonic approximation (QHA), including explicit treatment of thermal expansion, which we find to be non-negligible for EMs. The lattice energies obtained with PBE0+MBD are the closest to experimental sublimation enthalpies with a mean absolute error of 9.89 kJ/mol. However, the state-of-the-art treatment of vibrational and thermal contributions makes the agreement with experiment worse. Pressure-volume curves are also examined for six representative materials. For pressure-volume curves, all three methods provide reasonable agreement with experimental data with mean absolute relative errors of 3% or less. Most of the intermolecular interactions typical of EMs, namely nitro-amine, nitro-nitro, and nitro-hydrogen interactions, are more sensitive to the choice of the dispersion method than to the choice of the exchange-correlation functional. The exception is π-π stacking interactions, which are also very sensitive to the choice of the functional. Overall, we find that PBE+TS, PBE+MBD, and PBE0+MBD do not perform as well for energetic materials as previously reported for other classes of molecular crystals. This highlights the importance of testing dispersion-inclusive DFT methods for diverse classes of materials and the need for further method development.

Nitroethane at high density: an experimental and computational vibrational study

The vibrational spectrum of liquid and solid nitroethane was measured as a function of pressure. Both Raman scattering and absorption IR spectroscopies were applied to samples of nitroethane, statically compressed at ambient temperature to a maximum pressure of 8.0 GPa and 16.9 GPa, respectively. A new amorphous to crystalline transition pressure was found to lie between 1.59-1.63 GPa. Davydov splitting of internal modes into two components suggests two molecules associated with the unit cell, which is consistent with the DFT predictions made in a previous study. For most bands below 1200 cm-1, pressure induced mode progression was consistent with DFT predictions. Conversely, observed mode shifts in the 2950-3100 cm-1 region were generally stiffer than their DFT counterparts. A discontinuity in mode evolution between 3.7-4.3 GPa was observed for a number of modes and shown to coincide with hydrogen bond rearrangement in this pressure region. Preferred orientation and crystallite strain might explain the increased scatter between the various pressure induced mode shift cycles. Time intervals on the order of ∼30 h may be required between spectra, in order to give the crystallites time to equilibrate their strain.

Vibrational energy redistribution in crystalline nitromethane simulated by ab initio molecular dynamics

Ab initio molecular dynamics simulations (AIMD) are systematically performed to study the Vibrational Energy Redistribution (VER) in solid nitromethane (NM) by combining normal mode decomposition and short-time Fourier transform technique. After the selective excitations of all fourteen intramolecular vibrational modes above 400 cm⁻¹, four three-dimensional (3D) excitation and detected vibrational spectra are obtained. The evolution of the kinetic energy proportion of all vibrations are also given and discussed quantitatively. These results show that, as the daughter modes, NO₂ symmetric stretches, CH₃ stretches and bends are usually excited quickly and relatively conspicuously compared with the other vibrations. Interestingly, we found that, although the stretching vibration of the CN bond which is a bridge between the methyl and nitro group can not respond immediately to the selective excitations, it always accumulates the vibrational energy slowly and steadily. Then, the underlying mechanisms are discussed based on the response of vibrational modes in both the time and frequency domain. As a result, we found that anharmonic transfers following symmetry rules which involve the couplings assisted by the overtones and rotations, as well as the transfers among the adjacent modes, play important roles in the VER of solid NM.

Ab initio molecular dynamics simulations (AIMD) are systematically performed to study the Vibrational Energy Redistribution (VER) in solid nitromethane (NM) by combining normal mode decomposition and short-time Fourier transform technique.

Effects of molecular vacancy and ethylenediamine on structural and electronic properties of CH3NO2 surfaces

The structural and electronic properties of (100) surface for nitromethane (NM) are studied using density functional theory (DFT) with the generalized gradient approximation and Perdew-Burke-Ernzerhof functional (GGA-PBE). Molecular vacancy and ethylenediamine (C₂H₈N₂) substitution are considered in this work. We find that ethylenediamine substitution significantly decreases the band gap, while molecular vacancy increases the band gap slightly. It indicates that ethylenediamine substitution has a positive effect on the impact sensitivity of NM. Also, the formation energies are calculated and the reasons for the decrease of band gap for ethylenediamine substitution and the increase of band gap for CH₃NO₂ vacancy are explained.

A systematic study of the surface structures and energetics of CH3NO2 surfaces by first-principles calculations

Density functional theory (DFT) has been employed within the generalized gradient approximation and Perdew-Burke-Ernzerhof functional (GGA-PBE) to study the structural and electronic properties of nitromethane (NM) surface models. Different surfaces, including (100), (001), (101), (110), and (111), are considered in this work. The corresponding properties of bulk crystal for NM were also calculated to form a contrast to the slab models. Results with anisotropic characteristics of different surfaces have been observed in this study. There was an obviously great anisotropy in electronic parameters, especially the band gaps of different surfaces, indicating the anisotropic impact sensitivity along different directions of NM. The band gap value for (111) surface, 2.687 eV, was smaller than that of other surfaces, showing a higher impact sensitivity for NM. The estimated anisotropy has been revealed in surface energies for different surfaces. Graphical Abstract The valence band minimum (VBM) and conduction band maximum (CBM) of the nitromethane (100), (001), (101), (110) and (111) surface models.

Effects of different dopant elements on structures, electronic properties, and sensitivity characteristics of nitromethane

In this study, the doped defects in nitromethane crystals were investigated using first-principles calculations for the first time. We introduce dopant atoms in the interstitial sites of the nitromethane lattice, aiming to study the effects of element-doping on the structural properties, electronic properties, and sensitivity characteristics. The obtained results show that doped defects obviously affect the neighboring nitromethane molecules. The modification of electronic properties shows that the band gaps are significantly influenced by doped defects. Partial density of states and population analysis further reveal the mechanism for sensitivity control of nitromethane. It is shown that the new electronic states were introduced in the forbidden bands and the doped defects resulted in charge redistributions in the systems. Graphical abstract The valence and conduction band edge positions as well as defect levels of pure and X-doped NM.

DFT-D studies on structural variation and absorption properties of crystalline benzotrifuroxan under high pressure

The structural, electronic, and absorption properties of crystalline benzotrifuroxan (BTF) under hydrostatic compression of 0–150 GPa have been studied by using density functional theory with dispersion correction. The crystal structure of BTF was relaxed using three types of van der Waals corrections such as the PBE-G06, PBE-TS, and PW91-OBS functionals at ambient conditions. The results indicate that PBE-G06 is the best functional for studying BTF. The a and c directions may be more sensitive to van der Waals interactions than the b direction. The structure is the stiffest in the b direction in the whole pressure range. At 107 GPa, BTF decomposes by the breaking of the N2−O2 and O1−N6 bonds. An analysis of the band gap and density of states indicates that BTF becomes more and more sensitive under compression. The absorption spectra show that BTF has relatively high optical activity with increasing pressure; moreover, the absorption region is broadened in the energy range of below 1.5 and above 16.5 eV, respectively. This work may provide useful information in understanding how BTF behaves under high pressure.

Self-enhanced catalytic activities of functionalized graphene sheets in the combustion of nitromethane: molecular dynamic simulations by molecular reactive force field

Functionalized graphene sheet (FGS) is a promising additive that enhances fuel/propellant combustion, and the determination of its mechanism has attracted much interest. In the present study, a series of molecular dynamic simulations based on a reactive force field (ReaxFF) are performed to explore the catalytic activity (CA) of FGS in the thermal decay of nitromethane (NM, CH3NO2). FGSs and pristine graphene sheets (GSs) are oxidized in hot NM liquid to increase their functionalities and subsequently show self-enhanced CAs during the decay. The CAs result from the interatomic exchanges between the functional groups on the sheets and the NM liquid, i.e., mainly between H and O atoms. CA is dependent on the density of NM, functionalities of sheets, and temperature. The GSs and FGSs that originally exhibit different functionalities tend to possess similar functionalities and consequently similar CAs as temperature increases. Other carbon materials and their oxides can accelerate combustion of other fuels/propellants similar to NM, provided that they can be dispersed and their key reaction steps in combustion are similar to NM.

DFT study of pressure effects in molecular crystal 4,10-dinitro-2,6,8,12-tetraoxa-4,10-diazatetracyclo-[5.5.0.05,903,11]-dodecane

Density functional theory was used to study the structural, electronic, and optical properties of crystalline 4,10-dinitro-2,6,8,12-tetraoxa-4,10-diazatetracyclo-[5.5.0.05,903,11]-dodecane (TEX) under hydrostatic pressure. The results indicate that there is a displacive transition in TEX under compression that has never been found in experiments. As the pressure increases, the band gap gradually decreases but presents an abnormal increase at 61 GPa, called the structural transition; moreover, the gap reduction is more pronounced in the low-pressure range compared with the high-pressure range. An analysis of density of states shows that the electronic delocalization in TEX is enhanced gradually under the influence of pressure. The peaks of the imaginary parts of the dielectric functions, energy-loss function, and reflectivity may come mainly from the electron transitions between the oxygen 2p and nitrogen 2p states. The electron energy-loss function presents a blue shift under compression. TEX has relatively higher optical activity at high pressure than at ambient conditions.

A DFT study on structural, vibrational properties, and quasiparticle band structure of solid nitromethane

We report a detailed theoretical study of the structural and vibrational properties of solid nitromethane using first principles density functional calculations. The ground state properties were calculated using a plane wave pseudopotential code with either the local density approximation, the generalized gradient approximation, or with a correction to include van der Waals interactions. Our calculated equilibrium lattice parameters and volume using a dispersion correction are found to be in reasonable agreement with the experimental results. Also, our calculations reproduce the experimental trends in the structural properties at high pressure. We found a discontinuity in the bond length, bond angles, and also a weakening of hydrogen bond strength in the pressure range from 10 to 12 GPa, picturing the structural transition from phase I to phase II. Moreover, we predict the elastic constants of solid nitromethane and find that the corresponding bulk modulus is in good agreement with experiments. The calculated elastic constants show an order of C11> C22 > C33, indicating that the material is more compressible along the c-axis. We also calculated the zone center vibrational frequencies and discuss the internal and external modes of this material under pressure. From this, we found the softening of lattice modes around 8-11 GPa. We have also attempted the quasiparticle band structure of solid nitromethane with the G0W0 approximation and found that nitromethane is an indirect band gap insulator with a value of the band gap of about 7.8 eV with G0W0 approximation. Finally, the optical properties of this material, namely the absorptive and dispersive part of the dielectric function, and the refractive index and absorption spectra are calculated and the contribution of different transition peaks of the absorption spectra are analyzed. The static dielectric constant and refractive indices along the three inequivalent crystallographic directions indicate that this material has a considerable optical anisotropy.

First-principles study of the structural transformation, electronic structure, and optical properties of crystalline 2,6-diamino-3,5-dinitropyrazine-1-oxide under high pressure

Periodic first-principles calculations have been performed to study the effect of high pressure on the geometric, electronic, and absorption properties of 2,6-diamino-3,5-dinitropyrazine-1-oxide (LLM-105) under hydrostatic pressures of 0-50 GPa. Obvious irregular changes in lattice constants, unit-cell angles, bond lengths, bond angles, and band gaps showed that crystalline LLM-105 undergoes four structural transformations at 8, 17, 25, and 42 GPa, respectively. The intramolecular H-bonds were strong at pressures of 0-41 GPa but weakened in the range 42-50 GPa. The lengths of the intermolecular H-bonds (<1.47 Å) indicated that these H-bonds have covalent character and tend to induce the formation of a new twelve-membered ring. Analysis of the DOS showed that the interactions between electrons, especially the valence electrons, strengthen under the influence of pressure. The p states play a very important role in chemical reactions of LLM-105. The absorption spectrum of LLM-105 displayed more bands--as well as stronger bands--in the fundamental absorption region when the pressure was high rather than low. A new absorption peak due to O-H stretching appeared at 18.3 eV above 40 GPa, indicating that covalent O-H bonds and a new twelve-membered ring are present in LLM-105.

Periodic density functional theory study of the high-pressure behavior of energetic crystalline 1,4-dinitrofurazano[3, 4-b]piperazine

A detailed study of the structural, electronic, and optical absorption properties of crystalline 1,4-dinitrofurazano[3,4-b]piperazine (DNFP) under hydrostatic pressures of 0-100 GPa was performed using periodic density functional theory. As the pressure increases, the lattice constants and cell volumes calculated by LDA gradually approach those obtained by GGA-PW91. It was found that the structure of DNFP is much stiffer in the b direction than along the a and c axes, indicating that the compressibility of the crystal is anisotropic. As the pressure increases, the band gap gradually decreases, and this decrease is more pronounced in the low-pressure range than in the high-pressure region. An analysis of the density of states showed that the electronic delocalization in DNFP gradually increases under the influence of pressure. DNFP exhibits relatively high optical activity at high pressure. As the pressure increases, the bands in the fundamental absorption region of the absorption spectrum of DNFP become more numerous and intense.

Reactive molecular dynamics simulation of solid nitromethane impact on (010) surfaces induced and nonimpact thermal decomposition

Which is the first step in the decomposition process of nitromethane is a controversial issue, proton dissociation or C-N bond scission. We applied reactive force field (ReaxFF) molecular dynamics to probe the initial decomposition mechanisms of nitromethane. By comparing the impact on (010) surfaces and without impact (only heating) for nitromethane simulations, we found that proton dissociation is the first step of the pyrolysis of nitromethane, and the C-N bond decomposes in the same time scale as in impact simulations, but in the nonimpact simulation, C-N bond dissociation takes place at a later time. At the end of these simulations, a large number of clusters are formed. By analyzing the trajectories, we discussed the role of the hydrogen bond in the initial process of nitromethane decompositions, the intermediates observed in the early time of the simulations, and the formation of clusters that consisted of C-N-C-N chain/ring structures.

MOLECULAR DYNAMIC SIMULATIONS OF SOLID NITROMETHANE UNDER HIGH PRESSURES

We report ReaxFF molecular dynamic simulations of structure change of crystalline nitromethane and the formation of hydrogen bond under high pressure. Under high pressure, the angles between C–N bonds and X, Y and Z axes have changed. Through the calculation of g(r) of O and H atoms, we found a new peak near 1.6 Å, which indicates the formation of the hydrogen bond between O and H atoms. We calculated the distribution of the angles of the C–N bonds orientations, the distribution of the dihedral angle of CNOO , and the charge distribution of nitromethane molecules under various pressures, and made a comparison between low and high pressures. The effects of hydrogen bonding in high explosive materials are discussed.

Terahertz normal mode relaxation in pentaerythritol tetranitrate

Normal vibrational modes for a three-dimensional defect-free crystal of the high explosive pentaerythritol tetranitrate were obtained in the framework of classical mechanics using a previously published unreactive potential-energy surface [J. Phys. Chem. B 112, 734 (2008)]. Using these results the vibrational density of states was obtained for the entire vibrational frequency range. Relaxation of selectively excited terahertz-active modes was studied using isochoric-isoergic (NVE) molecular dynamics simulations for energy and density conditions corresponding to room temperature and atmospheric pressure. Dependence of the relaxation time on the initial modal excitation was considered for five excitation energies between 10 and 500 kT and shown to be relatively weak. The terahertz absorption spectrum was constructed directly using linewidths obtained from the relaxation times of the excited modes for the case of 10 kT excitation. The spectrum shows reasonably good agreement with experimental results. Dynamics of redistribution of the excited mode energy among the other normal modes was also studied. The results indicate that, for the four terahertz-active initially excited modes considered, there is a small subset of zero wave vector (k = 0) modes that preferentially absorb the energy on a few-picosecond time scale. The majority of the excitation energy, however, is transferred nonspecifically to the bath modes of the system.

Thermal decomposition of the solid phase of nitromethane: ab initio molecular dynamics simulations

The Car-Parrinello molecular dynamics simulations were employed to investigate thermal decomposition of the solid nitromethane. It is found that it undergoes chemical decomposition at about 2200 K under ambient pressure. The initiation of reactions involves both proton transfer and commonly known C-N bond cleavage. About 75 species and 100 elementary reactions were observed with the final products being H2O, CO2, N2, and CNCNC. It represents the first complete simulation of solid-phase explosive reactions reported to date, which is of far-reaching implication for design and development of new energetic materials.

Structure of liquid nitromethane: Comparison of simulation and diffraction studies

Simulation (molecular dynamics and Car-Parrinello [Phys. Rev. Lett. 55, 2471 (1985)]) and diffraction (x-ray and neutron) studies on nitromethane are compared aiming at the determination of the liquid structure. Beyond that, the capabilities of the methods to describe liquid structure are discussed. For the studied liquid, the diffraction methods are performing very well in the determination of intramolecular structure, but they do not give detailed structural information on the intermolecular structure. The good agreement between the diffraction experiments and the results of molecular dynamics simulations justifies the use of simulations for the more detailed description of the liquid structure using partial radial distribution functions and orientational correlation functions. Liquid nitromethane is described as a molecular liquid without strong intermolecular interactions such as hydrogen bonding, but with detectable orientational correlations resulting in preferential antiparallel order of the neighboring molecules.

Ab initio studies of crystalline nitromethane under high pressure

We have studied the mechanical compressibility and band structure of solid nitromethane both in equilibrium and compressed states using Hartree-Fock and density functional theory (DFT) with atom-centered all-electron linear combination of atomic orbitals basis sets. Hartree-Fock calculations with a 6-21G basis set, uncorrected for basis set superposition error, gave the best agreement with experimental compression studies. These results may be due to the cancellation of basis set superposition error with dispersion force errors. The equilibrium DFT band gap is comparable to the lowest-energy feature in electron-impact spectroscopy of nitromethane but underpredicts the optical absorption gap; we interpret these features in terms of the presence of tightly bound excitons. Only minor changes in the gap are observed under hydrostatic compression.