Next Article in Journal
Photosynthesis 2.0
Previous Article in Journal
Molecular Mechanisms of the Protective Effects of Olive Leaf Polyphenols against Alzheimer’s Disease
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Theoretically Revealing the Response of Intermolecular Vibration Energy Transfer and Decomposition Process of the DNTF System to Electric Fields Using Two-Dimensional Infrared Spectra

1
Xi’an Modern Chemistry Research Institute, Xi’an 710065, China
2
Department of Physics and Engineering Physics, University of Saskatchewan, Saskatoon, SK S7N 5E2, Canada
3
College of Food Science and Engineering, Henan University of Technology, Zhengzhou 450001, China
4
College of Life Science and Biotechnology, Shanghai Jiao Tong University, Shanghai 200240, China
5
National Key Laboratory for Shock Wave and Detonation Physics Research, Institute of Fluid Physics, Chinese Academy of Engineering Physics, Mianyang 621000, China
*
Authors to whom correspondence should be addressed.
Int. J. Mol. Sci. 2023, 24(5), 4352; https://doi.org/10.3390/ijms24054352
Submission received: 10 December 2022 / Revised: 12 February 2023 / Accepted: 20 February 2023 / Published: 22 February 2023
(This article belongs to the Topic Theoretical, Quantum and Computational Chemistry)

Abstract

:
The external electric field (E-field), which is an important stimulus, can change the decomposition mechanism and sensitivity of energetic materials. As a result, understanding the response of energetic materials to external E-fields is critical for their safe use. Motivated by recent experiments and theories, the two-dimensional infrared (2D IR) spectra of 3,4-bis (3-nitrofurazan-4-yl) furoxan (DNTF), which has a high energy, a low melting point, and comprehensive properties, were theoretically investigated. Cross-peaks were observed in 2D IR spectra under different E-fields, which demonstrated an intermolecular vibration energy transfer; the furazan ring vibration was found to play an important role in the analysis of vibration energy distribution and was extended over several DNTF molecules. Measurements of the non-covalent interactions, with the support of the 2D IR spectra, indicated that there were obvious non-covalent interactions among different DNTF molecules, which resulted from the conjugation of the furoxan ring and the furazan ring; the direction of the E-field also had a significant influence on the strength of the weak interactions. Furthermore, the calculation of the Laplacian bond order, which characterized the C-NO2 bonds as trigger bonds, predicted that the E-fields could change the thermal decomposition process of DNTF while the positive E-field facilitates the breakdown of the C-NO2 in DNTFⅣ molecules. Our work provides new insights into the relationship between the E-field and the intermolecular vibration energy transfer and decomposition mechanism of the DNTF system.

Graphical Abstract

1. Introduction

The response of energetic materials to an accidental stimulus such as impact, shock, or heat, which can trigger an unintended detonation, is a persistent concern with regard to the safe usage of energetic materials. The electric field (E-field) is one of the most important external stimuli that has been employed as an electrically controlled approach for the ignition of energetic materials [1]. Several reports have indicated that the introduction of an external E-field into energetic materials can increase the amount of energy and consequently result in an increase in detonation pressure and velocity [2,3,4]. Tasker et al. reported that the application of a parallel E-field (0.02 V/nm) caused a significant reduction in the failure thickness of the explosive PBX-9502 [2]. In similar experiments, Lee studied the effects of the E-field on the initiation and growth of a reaction in an HMX-based, cast-cured explosive with the use of the modified gap test; modest effects were observed in the test [3,4]. In the actual experiment, the external E-field was imposed on energetic materials that were considered hazardous. In addition to the dangers that might result from experimental operations and equipment, such as strong E-field installations, high-energy molecules may be more unstable in the presence of an external E-field. Computational techniques have recently become popular for determining the behaviors of energetic materials in the presence of an external E-field [5,6,7,8,9,10,11,12]. Polizer and co-workers systematically studied the influences of an E-field (−1.5426~1.5426 V/nm) on the “trigger linkage” bonds, including C-NO2, N-NO2, and O-NO2, in nitromethane and dimethylnitramine. The results indicated that these bonds could be strengthened along the E-field direction through the reinforcement of intrinsic molecular polarities [5,6,7]. Ren et al. theoretically established linear correlations between the strength of the E-field (−5.142~5.142 V/nm) and the structural parameters, including bond length and bond stretching frequency, and demonstrated that the E-field could trigger the rupture of the N-O bond at the initial stage of decomposition [8,9,10]. Liu et al. investigated the structural arrangement of nitromethane under an E-field ranging from 0 to 5 V/nm with the use of molecular dynamics (MD). They found that the critical E-field strength was between 2 and 3 V/nm, which could induce the transition of nitromethane molecules from a relatively disordered distribution to a solid-like, ordered, and compact arrangement with a large density [11]. Wood et al. studied the initial reactions and subsequent decomposition of α-HMX, which was excited via an E-field ranging from 0.56 V/nm to 3.56 V/nm at various frequencies through MD with the reactive potential ReaxFF. Their results revealed that the energy increased with an increase in the rate of the energy input and plateaued as the process became athermal for high loading rates [12]. However, there is still a lack of understanding regarding the role of the E-field in the intermolecular interactions of energetic materials.
As a promising high-density material, 3,4-bis (3-nitrofurazan-4-yl) furoxan (DNTF), with its high energy, low melting point, and excellent comprehensive properties, has become the focus of researchers studying energetic materials in recent years [13,14,15,16,17,18,19,20]. It was found that a coordinated oxygen atom exists in the furoxan ring and that the C-NO2 bond is the longest and most likely to be broken down upon heating or impact [13,14]. Nan et al. conducted experiments that determined the thermal decomposition properties of DNTF in the condensed phase. They employed differential scanning calorimetry and fast scanning Fourier transform infrared (FTIR) spectroscopy and found that in the thermal decomposition process, C-NO2 broke down first, followed by the breakdown of the C-C bond connecting the furazan and furoxan rings [15]. The fast thermolysis of DNTF at 0.1~0.4 Mpa was investigated using temperature-jump FTIR spectroscopy, and the results revealed that the two competitive reactions, C-NO2 homolysis (to form NO2) and isomerization (to form NO), could occur during the fast thermolysis of DNTF [16]. Theoretical and experimental works demonstrated that the loss of the nitro group (NO2) was an important step in the shock initiation of DNTF, and that the intermolecular non-covalent interactions contributed to the stability of DNTF.
The rate of energy transfer in the presence and absence of intermolecular non-covalent interactions is of interest to the greater explosive community, and it is important to understand the nature of these interactions in explosive materials [21,22,23,24,25,26,27,28]. Ostrander et al. reported two-dimensional infrared (2D IR) measurements of pentaerythritol tetranitrate (PETN) thin films that resolved vibrational modes and were coherently delocalized over at least 15–30 molecules through the measurement of the transition dipole strength [23]. Shi et al. investigated the structural dynamics and the vibrational energy transfer dynamics of hexahydro-1,3,5-trinitro-1,3,5-triazine (RDX) using 2D IR spectroscopy, and the results indicated that vibrational energy transfer was more efficient in RDX microcrystals than in the solvated form [24]. Using theoretical, quasi-static 2D IR spectra, we recently calculated the vibrational energy transfer rate between 2,4,6-trinitrotoluene (TNT) and 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12-hexaazaisowurtzitane (CL-20) in the TNT/CL-20 cocrystal [25]. Considering the proposed importance of vibrational energy transfer in energy materials, 2D IR was especially well-suited for providing insight into the molecular structure and vibrational dynamics that may impact explosive design [26,27,28].
In this manuscript, we theoretically reveal the vibrational energy transfer process in the DNTF system using 2D IR spectra under different E-fields, with the use of the asymmetric NO2 stretching modes as the vibrational probes. Moreover, the stability of the trigger C-NO2 bonds and the non-covalent interactions were studied, respectively, using the Laplacian bond order (LBO) method [29] and the independent gradient model based on the Hirshfeld partition (IGMH) method [30]. This investigation serves as a theoretical reference for the efficient addition of an E-field to the DNTF system as well as a means to avoid a catastrophic explosion in an external E-field.

2. Results and Discussion

The theoretical and experimental geometries of DNTF in the absence of an E-field are presented in Table 1 [13,14], and the theoretical geometries of DNTF under different E-fields are given in Tables S3–S6. It can be seen that the bond lengths of N7O2 are well consistent with the experimental values of 1.207 Å and 1.217 Å, while the bond lengths of N8O2 vary greatly and are slightly smaller than the experimental values. The theoretical bond lengths of the furoxan ring range from 1.301 Å to 1.477 Å, which agree with the experimental scope of 1.302~1.440 Å. Compared with the experimental values, the optimized furoxan rings are within the reasonable scope of 1.295 Å~1.433 Å. Moreover, the theoretical trigger bond lengths of C-N show excellent agreement with the experimental values of 1.459 Å and 1.442 Å. Finally, the average deviation between the values for the theoretical DNTF molecules and the experimental optimization is calculated, and the results indicate that DNTFⅠ has the smallest value at 0.0128 Å, while DNTFⅡ has the largest value at 0.0144 Å.
Due to the fact that NO2 asymmetric stretching modes ranging from 1500 cm−1 to 1700 cm−1 are often used in the vibrational probes of energetic materials [21,22], details of the vibrational assignment of the NO2 asymmetric stretching vibration are presented in Figure 1 and Table S2. The theoretical anharmonic IR peak at 1640.8 cm−1 is consistent with our experimental value of 1639.5 cm−1 and other reports [13,14,17,18], and it is attributed to the NO2 asymmetric and furazan ring stretching vibration in DNTFⅣ; however, there is an obvious difference in the relative strength. The theoretical anharmonic IR peak at 1612.0 cm−1 is ascribed to the NO2 asymmetric and furazan ring stretching vibration in DNTFⅡ, which is in agreement with the experimental value of 1612.5 cm−1 [13,14,17,18]. Moreover, there is a relatively strong experimental peak at 1586.4 cm−1 [13,14,17,18] that is confirmed by the theoretical IR peak at 1579.7 cm−1. This is contributed by the NO2 asymmetric and furazan ring stretching vibrations in DNTFⅡ and the NO2 asymmetric stretching vibrations in DNTFⅣ. However, some theoretical IR peaks have not been reported in experiments. For example, the theoretical peak at 1623.2 cm−1 is attributed to the NO2 asymmetric and furazan ring stretching vibration in DNTFⅢ and the furoxan ring stretching vibration in DNTFⅠ. The theoretical peak at 1640.8 cm−1 is ascribed to the NO2 asymmetric and furazan ring stretching vibrations in DNTFⅣ. Through the above comparisons of the structures and IR peaks, the conclusion is that the optimized DNTF structure can well reproduce the experimental one, despite some differences.
The normalized quasi-static 2D IR spectra of NO2 asymmetric stretching modes under different E-fields are presented in Figure 2. The diagonal peaks of the 2D IR spectra reflect their corresponding infrared spectra, as shown in the top panel and bottom panel of Figure 2; their intensities are scaled with the transition dipole as μ [4,31,32]. As can be seen, a network of diagonal and off-diagonal peaks comprises a pair of positive (in blue) and negative (in yellow) signals, which are identified as the ground-state bleaching and simulated emission (v = 0→v = 1, where v is the vibration quantum number), and as the first excited-state absorption (v = 1→v = 2), respectively [33,34]. Due to relatively narrower absorption peaks, four main pairs of positive and negative signals are observed, and the positive signals usually contain several vibration modes from different molecules based on vibration energy distribution. However, the negative anharmonicity and small energy gap between the two NO2 modes in DNTF lead to the overlap of negative and positive peaks, while the E-field is 1.0284 V/nm (Figure 2d). The anharmonicity of the transition (ω = 1611.9 cm−1), which is attributed to the N7O2 asymmetric and the corresponding furazan ring stretching vibration in DNTFⅠ, is −19.8 cm−1, and the anharmonicity of the transition (ω = 1617.7 cm−1), which mainly comes from the NO2 asymmetric and furazan ring stretching vibration in DNTFⅡ, is 19.9 cm−1. Besides the diagonal peaks, cross-peaks are also seen in the upper left and lower right panels of Figure 2, which indicates that the vibrational energy transfer among the four eigenstates extends over multiple DNTF molecules. For example, when the E-field is 2.0568 V/nm (Figure 2e), there are four cross-peaks between the N8O2 asymmetric and the corresponding furazan ring stretching vibration in DNTFⅣ (ω = 1577.2 cm−1) and between the N7O2 asymmetric and the corresponding furazan ring stretching vibration in DNTFⅡ (ω = 1593.9 cm−1); this demonstrates that although DNTFⅣ and DNTFⅡ are spatially separated by DNTFⅢ, the vibration energy transfer between them still exists. Conversely, when the E-field is −1.0284 V/nm (Figure 2b), the cross-peaks between the N7O2 asymmetric stretching modes in DNTFⅡ (ω = 1626.6 cm−1) and DNTFⅢ (ω = 1656.4 cm−1) disappear, which indicates that the vibrational energy transfer process between them is not possible. The reasons are as follows: For linear spectroscopy in a coupled dimer, it is reasonable to consider the one-exciton Hamiltonian [33,34,35,36]
H = ω i β ij β ij ω j
where β and ω represent the coupling constant and the frequency between chromophore i and j, respectively. The exciton eigenvalues are given as follows
E i , j = ω i + ω j 4 β ij 2 + ω i ω j 2 2
The eigenstates of Ei and Ej are
ϕ i = cos α 10 sin α | 01
ϕ j = sin α 10 cos α | 01
where α is the mixing angle
tan 2 α = 2 β ij ω j ω i
In the weak coupling regime, the coupling constants are smaller than the frequency splitting
β ij ω i ω j
Their mixing angle α is small, and the exciton states will be localized on the individual sites. Furthermore, the two excitonic states closely resemble the uncoupled local modes
E i , j ω i , j 2 β ij 2 ω j ω i
Then, the coupling constant βij can be obtained based on the perturbation theory
β ij Δ ij Δ i ω i ω j
where Δij and Δi represent the off-diagonal anharmonicity and the diagonal anharmonicity, respectively.
Based on the Fermi golden rule, the energy transfer rate κij is proportional to the square of the coupling constant βij, that is
κ ij β ij 2
From Equations (8) and (9), when the off-diagonal anharmonicity Δij is zero, that is, when the cross-peaks vanish, the energy transfer rate is also zero.
In addition, based on the vibration energy distribution and comparison of the 2D IR spectra, it can be concluded that the furazan ring stretching vibrations play an important role in intermolecular vibration energy transfer.
In the above discussions, it was established that the vibration energy transfer usually occurs among the vibration modes, including the furazan ring stretching vibration in DNTF. Considering the fact that intermolecular vibration energy transfer is the bridge of intermolecular interactions, we opted to use the IGMH method [30,37], where the Hirshfeld partitioning of the actual molecular electron density is used to derive the atomic densities in order to be able to analyze intermolecular non-covalent interactions in the DNTF system under different E-fields (Figure 3). To maximally reveal the interactions between the DNTF molecules, the isovalue was set to 0.01 a.u. in the maps. As expected, due to the conjugation of the furoxan and furazan rings, a series of thin and broad isosurfaces appear between adjacent DNTF molecules, which ideally exhibit π-π stacking interactions and can also explain the existence of cross-peaks in the 2D IR spectra. More specifically, there are mainly π-π stacking interactions between the nitro group and the furazan ring in DNTFⅠ and DNTFⅡ; between the furazan ring connecting with the N8O2 in DNTFⅡ and the furazan and furoxan rings in DNTFⅢ; and between the furoxan ring in DNTFⅢ and the furazan ring connecting with the N8O2 and the furoxan ring in DNTFⅣ. Moreover, a van deer Walls (vdW) interaction is observed between the nitro group in DNTFⅠ and that in DNTFⅢ. Interestingly, when the direction of the E-field is opposite to the X-axis (Figure 3a,b), the area of the green isosurface between DNTFⅢ and DNTFⅣ, especially between the furoxan ring in DNTFⅢ and the furazan ring connecting with the N8O2 in DNTFⅣ (Figure 3, top right), is larger than that when the E-field is along the X-axis (Figure 3d,e); this indicates that the direction of the E-field has a significant influence on the strength of weak interactions.
It is believed that a more sensitive energetic molecule generally includes a trigger bond with a smaller bond order calculated using quantum chemistry. Furthermore, the experimental results demonstrate that the thermal decomposition process in DNTF involves the C-NO2 bonds breaking down first [15,16,19,20]. Therefore, the stability of the C-NO2 bond as the trigger bond is investigated using the Laplacian bond order (LBO) [38]. The LBO, which corresponds to the Laplacian of electron density in the fuzzy overlap space, is more comprehensive in investigating the stability of covalent bonds and is defined between atoms A and B as follows [38]
2 ρ ( r ) = 2 ρ ( r ) x 2 + 2 ρ ( r ) y 2 + 2 ρ ( r ) z 2
LBO A , B = 10 × 2 ρ < 0 w A ( r ) w B ( r ) 2 ρ ( r ) d r
where ρ(r) denotes the electron density at point r, and wA(r) stands for an atomic weighting function and corresponds directly to the scope of the atomic space. The prefactor −10 is introduced to make the magnitude of the LBO consistent with one’s chemical intuition for typical covalent bonds. Generally, the value of the LBO is positively related to the stability of the covalent bonds. As the DNTF primitive cell contains four DNTF molecules [13,14], the LBO values of eight C-NO2 bonds are presented in Figure 4. Due to the orientation of the coordinated oxygen on the furoxan ring [13,14], the LBO values of C-N8O2 are usually larger than those of C-N7O2, which reveals that the C-N7O2 bonds tend to induce decomposition. More specifically, regardless of the E-field change, the breakdown of the C-N7O2 bond in DNTFⅠ is the first step in the thermal decomposition process. The LBO value of the C-N8O2 bond in DNTFⅡ is the largest and is almost unchanged, which may be attributed to the fact that N8O2 in DNTFⅡ has strong π-π stacking interactions with the furoxan and furazan rings in other molecules (Figure 3). When the E-field is at 1.0284 V/nm, the LBO value of the C-N7O2 bond in DNTFⅠ is the largest, demonstrating that the DNTF in this E-field is relatively more stable during the decomposition process. However, with an E-field change, the next steps in the decomposition process are competitive. For example, when the E-field is −2.0568 V/nm, the second and third steps involve the cleavage of the C-N7O2 bonds in DNTFⅡ and DNTFⅣ, respectively, which are the same as those in the absence of an E-field. Nevertheless, when the E-field is 2.0568 V/nm, the second and third steps in the decomposition process represent the cleavage of the C-N7O2 bonds in DNTFⅣ and that of the C-N8O2 bonds in DNTFⅢ, respectively.
The highest occupied molecular orbital (HOMO) and the lowest unoccupied molecular orbital (LUMO) play a significant role in chemical reactivity [39]. Particularly, their energy gap can reflect the chemical reactivity and optical polarizability of energy materials [40,41]. The distribution of the HOMO and LUMO of DNTF under E-fields, along with their energy gaps, is presented in Figure 5. It can be seen that the LUMO is mainly distributed on the N7O2 and the corresponding furazan ring in DNTFⅣ under negative E-fields, and that it is mainly distributed on the N8O2 and the corresponding furazan ring in DNTFⅡ in the absence of an E-field and under positive E-fields. The HOMO is contributed by the furoxan and furazan rings connecting with the N7O2 in DNTFⅡ in the absence of an E-field and under negative E-fields, but the HOMO is contributed by the furoxan and furazan rings connecting with the N7O2 in DNTFⅣ under positive fields. In addition, larger HOMO and LUMO energy gaps usually lead to lower chemical reactivity, which allows for molecules to have higher chemical stability [40,41]. It is also revealed that when the E-field is −2.0568 V/nm, the value of the energy gap is the smallest, indicating that the DNTF molecule is the most active; on the other hand, the value of the energy gap is the largest in the absence of an E-field, which suggests that the DNTF molecule has the highest chemical stability. Therefore, E-fields may be considered to have an important impact on the chemical stability and position of the HOMO and LUMO.

3. Methods and Materials

The molecule structure and atomic coordination of DNTF were taken from its crystal structure, determined by X-ray diffraction CCDC (Cambridge Crystallographic Data Center) number 665080 [13]. The primitive cell includes four DNTF molecules, wherein both nitro groups were determined to be twisted in relation to the adjacent furazan rings, and both furazan rings were twisted relative to the central furoxan ring (Figure 6 and Table S1). Geometric optimization and vibration frequency calculations of DNTF were performed via density function theory (DFT) at the level of B3LYP (D3), with the 6-311G (d, p) basis set as implemented in the Gaussian 16 program and the in-house BDF program [33,34,35,38,42,43,44,45,46]. The static E-field along with the X-axis (Figure 6, red arrow) was in the range of −2.0568~2.0568 V/nm, which has been extensively investigated [5,6,7,8,9,10,11,12]. The positive values mean that the direction of the E-field is along the X-axis, while the negative values indicate that the direction of the E-field is opposite to the X-axis. The absence of negative frequencies revealed by harmonic vibrational analysis confirms that the structure is a true minima. Then, through the vibration energy distribution analysis (VEDA) code [36], the frequencies were designated as the normal vibration modes. Anharmonic vibrational analysis was carried out on the reduced subspace of the asymmetric NO2 stretching normal modes. The anharmonic treatment of the vibrations was performed using the generalized second-order vibrational perturbation theory (GVPT2) implemented in Gaussian 16 [34,47], including terms up to the third and fourth derivatives of the potential energy with respect to the normal mode coordinates.
The calculation of quasi-static 2D IR spectra has been described in detail elsewhere [25,31]. In summary, the ground-state, two-exciton, and combination band frequencies of the asymmetric NO2 stretching modes were used to construct the two-exciton Hamiltonian, where the anharmonicity was calculated with the use of anharmonic frequencies; the two-exciton Hamiltonian was subsequently diagonalized. The transition dipole moment vectors between the vibrational ground state and the one-exciton state for each of the normal modes were acquired from the DFT results [34]. Then, the harmonic approximation, that is, the corresponding ground state of one-exciton state transition dipole moment vector multiplied by √2 42, was used to calculate the transition dipole moment vector between the one-exciton and two-exciton states. Moreover, the waiting time Tw was set to 3.0 ps. In a given vibrational mode, the fourth power of the transition dipole moments involved is generally proportional to the 2D IR diagonal peak intensity [31,32].

4. Conclusions

In this work, with the use of DFT/B3LYP (D3)/6-311G (d, p) calculations, we performed a comprehensive investigation of the effects of E-fields on the explosive DNTF, which included the study of intermolecular vibration energy transfer processes and non-covalent interaction using 2D IR spectra and the IGMH method, respectively. When there was an absence of an E-field, the optimized structural parameters of the DNTF primitive cell containing four DNTF molecules were in good agreement with the available experimental results, and the anharmonic theoretical spectra of the NO2 asymmetric stretching well reproduced the experiment peaks. In combination with the vibration energy distribution, cross-peaks in the 2D IR spectra were observed, which revealed that intermolecular vibrational energy transfer is extended over the DNTF molecules under different E-fields and that the furazan ring vibrations play an important role in the intermolecular vibration transfer process. It could be seen on the isosurface of the weak interaction that obvious non-covalent interactions exist among different DNTF molecules due to the conjugation of the furoxan and furazan rings, which were confirmed by the cross-peaks in the 2D IR spectra. Furthermore, it was found that the direction of the E-field has a significant influence on the strength of the weak interactions. More importantly, the LBO value of the C-NO2 bonds acting as the trigger bonds under different E-fields indicated that the E-fields could change the thermal decomposition process of DNTF, with the non-covalent interactions significantly contributing to the relative sensitivity of energy materials and the positive E-field facilitating the breakdown of C-NO2 in DNTFⅣ molecules. When the E-field was 1.0284 V/nm, DNTF was relatively more stable during the decomposition process. Finally, measurements of the HOMO and LUMO energy gaps indicate that DNTF is activated by E-fields and that the strength of E-fields has a relatively significant influence on the chemical reactivity of DNTF.

Supplementary Materials

The supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms24054352/s1.

Author Contributions

Conceptualization, H.R., X.W. and G.J.; methodology, H.R., L.J. and X.J.; software, D.W.; validation, R.L. and J.T.; formal analysis, H.R. and L.J.; data curation, X.W. and G.J.; writing-original draft preparation, H.R.; writing-review and editing, L.J. and D.W.; visualization, X.J. and R.L.; supervision, J.T.; project administration, X.W. and G.J.; funding acquisition, J.T., X.W. and G.J. All authors have read and agreed to the published version of the manuscript.

Funding

This research was supported by National Natural Science Foundation of China (11802281) and Military Commission Science and Technology Commission (2022-JCJQ-ZD-128-12).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The data presented in this study are available upon request from the corresponding authors.

Acknowledgments

We are grateful to Hao Wang for the experimental infrared spectra of DNTF and to HZWTECH for providing the computation facilities.

Conflicts of Interest

H.R.; X.J.; J.T.; R.L. and X.W. are employees of the Xi’an Modern Chemistry Research Institute. Xi’an Modern Chemistry Research Institute had no role in the design of the study, in the collection, analysis, or interpretation of data, in the writing of the manuscript, or in the decision to publish the results.

References

  1. Chung, K.; Rozumov, E.; Kaminsky, D.; Anderson, P.; Cook, P.; Sawka, W.; McPherson, M.; Buescher, T. Development of Electrically controlled energetic materials. ESC Trans. 2013, 50, 59–66. [Google Scholar] [CrossRef]
  2. Tasker, D.G.; Whitley, V.H.; Mace, J.L.; Pemberton, S.J.; Lee, R.J. Electromagnetic effects on explosive reaction and plasma. In Proceedings of the 14th International Detonation Symposium, Coeur d’Alene, Idaho, 11–16 April 2010. [Google Scholar]
  3. Lee, R.J.; Forbes, J.W.; Tasker, D.G.; Orme, R.S. Changes in Blow-off-velocity observed in two explosives at the threshold for sustained ignition using the modified gap test. In Shock Compression of Condensed Matter 2009, Proceedings of the American Physical Society Topical Group on Shock Compression of Condensed Matter, Nashville, TN, USA, 28 June–3 July 2009; Furnish, M.D., Ed.; American Institute of Physics: Melville, NY, USA, 2009. [Google Scholar]
  4. Lee, R.J.; Forbes, J.W.; Palermo, E.; Wilson, W.H. Effects of electric fields on sensitivity of an HMX based explosive. In Shock Compression of Condensed Matter 2007, Proceedings of the Conference of the American Physical Society, Waikoloa, HI, USA, 24–29 June 2007; American Institute of Physics: Melville, NY, USA, 2007; pp. 963–966. [Google Scholar]
  5. Politzer, P.; Murray, J.S.; Lane, P. Computational determination of effects of electric fields upon “trigger linkages” of prototypical energetic molecules. Int. J. Quantum Chem. 2009, 109, 534–539. [Google Scholar] [CrossRef]
  6. Politzer, P.; Murray, J.S. Computed effects of electric fields upon the C-NO2 and N-NO2 bonds of nitromethane and dimethylnitramine. Int. J. Quantum Chem. 2009, 109, 3–7. [Google Scholar] [CrossRef]
  7. Politzer, P.; Murray, J.S.; Concha, M.C.; Lane, P. Effects of Electric Fields Upon Energetic Molecules: Nitromethane and Dimethylnitramine. Cent. Eur. J. Energet. Mater. 2007, 4, 3–21. [Google Scholar]
  8. Ren, F.D.; Cao, D.L.; Shi, W.J.; You, M. A dynamic prediction of stability for nitromethane in external electric field. RSC Adv. 2017, 7, 47063–47072. [Google Scholar] [CrossRef] [Green Version]
  9. Ren, F.D.; Cao, D.L.; Shi, W.J. A dynamics prediction of nitromethane → methyl nitrite isomerization in external electric field. J. Mol. Model. 2016, 22, 96. [Google Scholar] [CrossRef]
  10. Ren, F.D.; Cao, D.L.; Shi, W.J.; You, M.; Li, M. A theoretical prediction of the possible trigger linkage of CH3NO2 and NH2NO2 in an external electric field. J. Mol. Model. 2015, 21, 145. [Google Scholar] [CrossRef]
  11. Liu, Y.; Ma, Y.; Yu, T.; Lai, W.; Guo, W.; Ge, Z.; Ma, Z. Structural Rearrangement of Energetic Materials under an External Electric Field: A Case Study of Nitromethane. J. Phys. Chem. A 2018, 122, 2129–2134. [Google Scholar] [CrossRef]
  12. Wood, M.A.; Van Duin, A.C.; Strachan, A. Coupled thermal and electromagnetic induced decomposition in the molecular explosive αHMX; a reactive molecular dynamics study. J. Phys. Chem. A 2014, 118, 885–895. [Google Scholar] [CrossRef]
  13. Zhou, Y.; Zhang, Z.; Li, J.; Guan, X.; Huang, X.; Zhou, C. Crystal Structure of 3,4-dinitrofurazanofuroxan. Chin. J. Explos. Propell. 2005, 28, 43–52. [Google Scholar]
  14. Sheremerev, A.V.; Ivanova, E.A.; Spiridonova, N.P.; Melnikova, S.F.; Antipin, M.Y. Desilylative nitration of C,N-disilylated 3-amino-4-methylfurazan. J. Heterocycl. Chem. 2005, 42, 1237. [Google Scholar] [CrossRef]
  15. Nan, H.; Pan, Q.; Jiang, F.; Wang, X.; Jin, P.; Zhang, K. In situ infrared thermal decomposition kinetics and mechanism of DNTF. Chin. J. Energet. Mater. 2022, 30, 604–610. [Google Scholar]
  16. Ren, X.; Liu, Z.; Wang, X.; Zhao, X.F.; Xie, M.; Heng, S. Investigated on the flash thermolysis of 3, 4-dinitrofuruzan-furoxan by T-jump/FTIR Spectroscopy. Acta Phys.-Chim. Sin. 2010, 26, 547–551. [Google Scholar]
  17. Ma, H.-X.; Song, J.-R.; Xiao, H.-M.; Hu, R.-Z.; Zhao, F.-Q. Density Functional Theoretical Investigation on 3, 4-Dinitrofurazanfuroxan (DNTF). Chin. J. Explos. Propell. 2006, 3, 43. [Google Scholar]
  18. Li, Y.; Yuan, J.; Zhao, W.; Qu, Y.; Xing, X.; Meng, J.; Liu, Y. Application and Development of 3, 4-Bis (3-nitrofurazan-4-yl) furoxan (DNTF). Russ. J. Gen. Chem. 2021, 91, 445–455. [Google Scholar] [CrossRef]
  19. Nan, H.; Jia, X.; Wang, X.; Liu, H.; Jiang, F.; Zhang, P. Mechanical and Electronic Properties of DNTF Crystals under Different Pressure. Crystals 2021, 11, 1180. [Google Scholar] [CrossRef]
  20. Zhang, G.; Jin, S.; Li, L.; Li, Y.; Li, Z.; Wang, D.; Zhang, B.; Jing, B.; Shu, Q. Thermal stability assessment of 3, 4-bis (3-nitrofurazan-4-yl) furoxan (DNTF) by accelerating rate calorimeter. J. Therm. Anal. Calorim. 2016, 126, 1185–1190. [Google Scholar] [CrossRef]
  21. Cole-Filipiak, N.C.; Marquez, M.; Knepper, R.; Harmon, R.; Wiese-Smith, D.; Schrader, P.; Wood, M.; Ramasesha, K. Ultrafast spectroscopic studies of vibrational energy transfer in energetic materials. In AIP Conference Proceedings; AIP Publishing LLC: Melville, NY, USA, 2020; Volume 2272, p. 060006. [Google Scholar]
  22. Cole-Filipiak, N.C.; Knepper, R.; Wood, M.; Ramasesha, K. Sub-picosecond to Sub-nanosecond Vibrational Energy Transfer Dynamics in Pentaerythritol Tetranitrate. J. Phys. Chem. Lett. 2020, 11, 6664–6669. [Google Scholar] [CrossRef]
  23. Ostrande, J.S.; Knepper, R.; Tappan, A.S.; Kay, J.J.; Zanni, M.T.; Farrow, D.A. Energy Transfer Between Coherently Delocalized States in Thin Films of the Explosive Pentaerythritol Tetranitrate (PETN) Revealed by Two-Dimensional Infrared Spectroscopy. J. Phys. Chem. B 2017, 121, 1352–1361. [Google Scholar] [CrossRef]
  24. Shi, L.; Yu, P.; Zhao, J.; Wang, J. Ultrafast Intermolecular Vibrational Energy Transfer in Hexahydro-1,3,5-Trinitro-1,3,5-Triazine in Molecular Crystal by 2D IR Spectroscopy. J. Phys. Chem. C 2020, 124, 2388–2398. [Google Scholar] [CrossRef]
  25. Ren, H.C.; Ji, L.X.; Jia, X.Z.; Liu, R.P.; Zhang, X.Q.; Wei, D.Q.; Wang, X.F.; Ji, G.F. Intermolecular vibration energy transfer process in two CL-20-based cocrystals theoretically revealed by two-dimensional infrared spectra. Molecules 2022, 27, 2153. [Google Scholar] [CrossRef]
  26. Bernstein, J. Ab initio study of energy transfer rates and impact sensitivities of crystalline explosives. J. Chem. Phys. 2018, 148, 084502. [Google Scholar] [CrossRef]
  27. Ramasesha, K.; Wood, M.; Cole-Filipiak, N.C.; Knepper, R. Experimental and Theoretical Studies of Ultrafast Vibrational Energy Transfer Dynamics in Energetic Materials; Sandia National Lab: Livermorel, CA, USA, 2020. [Google Scholar]
  28. Joshi, K.; Losada, M.; Chaudhuri, S. Intermolecular energy transfer dynamics at a hot-spot interface in RDX crystals. J. Phys. Chem. A 2016, 120, 477–489. [Google Scholar] [CrossRef]
  29. Lu, T.; Chen, F. Bond order analysis based on the Laplacian of electron density in Fuzzy overlap space. J. Phys. Chem. A 2013, 117, 3100–3108. [Google Scholar] [CrossRef]
  30. Lu, T.; Chen, Q. Independent gradient model based on Hirshfeld partition: A new method for visual study of interactions in chemical system. J. Comput. Chem. 2022, 43, 539–555. [Google Scholar] [CrossRef]
  31. Hamm, P.; Zanni, M.T. Concepts and Methods of 2D Infrared Spectroscopy; Cambridge University Press: Cambridge, UK, 2011. [Google Scholar]
  32. Kraack, J.P. Ultrafast Structural Molecular Dynamics Investigated with 2D Infrared Spectroscopy Methods. In Multidimensional Time-Resolved Spectroscopy; Springer: Berlin/Heidelberg, Germany, 2019; pp. 113–205. [Google Scholar]
  33. Becke, A.D. Density-Functional Exchange Approximation with Correct Asymptotic Behavior. Phys. Rev. A. 1988, 38, 3098. [Google Scholar] [CrossRef]
  34. Frisch, M.J.; Trucks, G.W.; Schlegel, H.B.; Scuserria, G.E.; Robb, M.A.; Cheeseman, J.R.; Schlmani, G.; Barone, V.; Mennucci, B.; Petersson, G.A.; et al. Gaussian 16, Revision A. 03; Gaussian, Inc.: Wallingford, CT, USA, 2016. [Google Scholar]
  35. Lee, C.; Yang, W.; Parr, R.G. Development of the Colle-salvetti Correlation-energy Formula into a Functional of the Electron Density. Phys. Rev. B 1988, 37, 785–789. [Google Scholar] [CrossRef] [Green Version]
  36. Jamróz, M.H. Vibrational Energy Distribution Analysis (VEDA): Scopes and limitations. Spectrochim. Acta A 2013, 114, 220–230. [Google Scholar] [CrossRef]
  37. Humphrey, W.; Dalke, A.; Schulten, K. VMD: Visual Molecular Dynamics. J. Mol. Graph. Model. 1996, 14, 33–38. [Google Scholar] [CrossRef]
  38. Goerigk, L.; Grimme, S. A Thorough Benchmark of Density Functional Methods for General Main Group Thermochemistry, Kinetics, and Noncovalent Interactions. Phys. Chem. Chem. Phys. 2011, 13, 6670–6688. [Google Scholar] [CrossRef]
  39. He, P.; Zhang, J.-G.; Wang, K.; Yin, X.; Jin, X.; Zhang, T.-L. Extensive theoretical studies on two new members of the FOX-7 family: 5-(dinitromethylene)-1,4-dinitramino-tetrazole and 1,1′-dinitro-4,4′-diamino-5,5′-bitetrazole as energetic compounds. Phys. Chem. Chem. Phys. 2015, 17, 5840–5848. [Google Scholar] [CrossRef]
  40. Zhai, D.; Wang, J.; Hao, L.; Ma, C.; Ma, P.; Pan, Y.; Jiang, J. Molecular design and properties of bridged energetic pyridines derivatives. RSC Adv. 2019, 9, 37747–37758. [Google Scholar] [CrossRef] [Green Version]
  41. Stobiecka, A.; Sikora, M.; Bonikowski, R.; Kula, J. An exploratory study on the peroxyl-radical-scavenging activity of 2,6-dimethyl-5-hepten-2-ol and its heterocyclic analogues. J. Mol. Struct. 2016, 1107, 82–90. [Google Scholar] [CrossRef]
  42. Zhang, Y.; Suo, B.; Wang, Z.; Zhang, N.; Li, Z.; Lei, Y.; Zou, W.; Gao, J.; Peng, D.; Pu, Z.; et al. BDF: A relativistic electronic structure program package. J. Chem. Phys. 2020, 152, 064113. [Google Scholar] [CrossRef]
  43. Fast, P.L.; Corchado, J.; Sanchez, M.L.; Truhlar, D.G. Optimized Parameters for Scaling Correlation Energy. J. Phys. Chem. A 1999, 103, 3139–3143. [Google Scholar] [CrossRef]
  44. Grimme, S. Density Functional Theory with London Dispersion Corrections. Wiley Comput. Mol. Sci. 2011, 1, 211–228. [Google Scholar] [CrossRef]
  45. Caldeweyher, E.; Mewes, J.; Ehlert, S.; Grimme, S. Extension and evaluation of the D4 London-dispersion model for periodic systems. Phys. Chem. Chem. Phys. 2020, 22, 8499–8512. [Google Scholar] [CrossRef]
  46. Caldweyher, E.; Ehlert, S.; Hansen, A.; Neugebauer, H.; Soicher, S.; Spicher, S.; Bannwarth, C.; Grimme, S. A generally applicable atomic-charge dependent London dispersion correction. J. Chem. Phys. 2019, 150, 154122. [Google Scholar] [CrossRef]
  47. Kesharwani, M.K.; Brauer, B.; Martin, J.M.L. Frequency and Zero-Point Vibrational Energy Scale Factors for Double-Hybrid Densith Functionals (and Other Selected Methods): Can Anharmonic Force Fields Be Avoided? J. Phys. Chem. A 2015, 119, 1701–1714. [Google Scholar] [CrossRef]
Figure 1. The experimental and theoretical IR spectra of DNTF molecules; the theoretical anharmonic IR spectra only include NO2 asymmetric stretching modes.
Figure 1. The experimental and theoretical IR spectra of DNTF molecules; the theoretical anharmonic IR spectra only include NO2 asymmetric stretching modes.
Ijms 24 04352 g001
Figure 2. Normalized and simulated quasi-static 2D IR spectra and the corresponding IR spectra (top panel) of NO2 asymmetric stretching vibrations at different E-fields (V/nm). E-fields: (a) −2.0568; (b) −1.0284; (c) 0; (d) 1.0284; (e) 2.0568. Roman numbers represent the vibration modes in the corresponding DNTF molecules.
Figure 2. Normalized and simulated quasi-static 2D IR spectra and the corresponding IR spectra (top panel) of NO2 asymmetric stretching vibrations at different E-fields (V/nm). E-fields: (a) −2.0568; (b) −1.0284; (c) 0; (d) 1.0284; (e) 2.0568. Roman numbers represent the vibration modes in the corresponding DNTF molecules.
Ijms 24 04352 g002aIjms 24 04352 g002b
Figure 3. Intermolecular non-covalent interactions among DNTF molecules under different E-fields. The green isosurface represents the 0.01 a.u. region.
Figure 3. Intermolecular non-covalent interactions among DNTF molecules under different E-fields. The green isosurface represents the 0.01 a.u. region.
Ijms 24 04352 g003aIjms 24 04352 g003b
Figure 4. The LBO value of C-NO2 bonds in the DNTF system under different E-fields.
Figure 4. The LBO value of C-NO2 bonds in the DNTF system under different E-fields.
Ijms 24 04352 g004
Figure 5. Variation trends for the ΔEHOMO-LUMO of DNTF under different E-fields.
Figure 5. Variation trends for the ΔEHOMO-LUMO of DNTF under different E-fields.
Ijms 24 04352 g005
Figure 6. The optimized structure of the DNTF primitive cell includes four molecules, which are numbered using Roman numbers. The carbon (C), oxygen (O), and nitrogen (N) atoms are represented by blue-green, red, and blue-colored balls, respectively. Atoms in DNTFⅣ are labeled with elements and numbers, and the same rule is followed in other DNTF molecules.
Figure 6. The optimized structure of the DNTF primitive cell includes four molecules, which are numbered using Roman numbers. The carbon (C), oxygen (O), and nitrogen (N) atoms are represented by blue-green, red, and blue-colored balls, respectively. Atoms in DNTFⅣ are labeled with elements and numbers, and the same rule is followed in other DNTF molecules.
Ijms 24 04352 g006
Table 1. Theoretical and experimental bond lengths (Å) of DNTF in the absence of an E-field.
Table 1. Theoretical and experimental bond lengths (Å) of DNTF in the absence of an E-field.
BondDNTFIDNTFIIDNTFIIIDNTFIVExperiment
N7-O61.2131.2111.2181.2111.207
N7-O51.2171.2211.2181.2211.217
C1-N71.4671.4641.4591.4641.459
C1-C21.4331.4291.4261.4291.421
C1-N11.2941.2971.3001.2971.289
N1-O11.3601.3571.3501.3561.372
N2-O11.3641.3671.3671.3681.375
N2-C21.3111.3111.3081.3111.303
C2-C31.4411.4411.4441.4421.445
C3-N31.3371.3371.3291.3361.336
N3-O41.2011.1981.2071.2011.212
C3-C41.4281.4241.4221.4211.405
N8-O71.2341.2161.2241.2241.255
N8-O81.2041.2201.2101.2111.227
C6-N81.4541.4471.4631.4611.442
C5-C61.4261.4241.4251.4241.415
C6-N61.2971.2951.2971.2971.296
O5-N61.3531.3511.3571.3561.366
O5-N51.3741.3751.3741.3741.379
N5-C51.3051.3051.3051.3041.304
C4-C51.464 1.465 1.463 1.465 1.472
C4-N41.303 1.304 1.302 1.301 1.302
O2-N41.349 1.348 1.354 1.353 1.372
O2-N31.465 1.477 1.469 1.470 1.440
Ave. Dev.0.01280.01440.01340.01290
Note: Ave . Dev . = i The i Exp i 2 N .
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Ren, H.; Ji, L.; Jia, X.; Tao, J.; Liu, R.; Wei, D.; Wang, X.; Ji, G. Theoretically Revealing the Response of Intermolecular Vibration Energy Transfer and Decomposition Process of the DNTF System to Electric Fields Using Two-Dimensional Infrared Spectra. Int. J. Mol. Sci. 2023, 24, 4352. https://doi.org/10.3390/ijms24054352

AMA Style

Ren H, Ji L, Jia X, Tao J, Liu R, Wei D, Wang X, Ji G. Theoretically Revealing the Response of Intermolecular Vibration Energy Transfer and Decomposition Process of the DNTF System to Electric Fields Using Two-Dimensional Infrared Spectra. International Journal of Molecular Sciences. 2023; 24(5):4352. https://doi.org/10.3390/ijms24054352

Chicago/Turabian Style

Ren, Haichao, Linxiang Ji, Xianzhen Jia, Jun Tao, Ruipeng Liu, Dongqing Wei, Xiaofeng Wang, and Guangfu Ji. 2023. "Theoretically Revealing the Response of Intermolecular Vibration Energy Transfer and Decomposition Process of the DNTF System to Electric Fields Using Two-Dimensional Infrared Spectra" International Journal of Molecular Sciences 24, no. 5: 4352. https://doi.org/10.3390/ijms24054352

APA Style

Ren, H., Ji, L., Jia, X., Tao, J., Liu, R., Wei, D., Wang, X., & Ji, G. (2023). Theoretically Revealing the Response of Intermolecular Vibration Energy Transfer and Decomposition Process of the DNTF System to Electric Fields Using Two-Dimensional Infrared Spectra. International Journal of Molecular Sciences, 24(5), 4352. https://doi.org/10.3390/ijms24054352

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop