ABSTRACT
CONTEXT AND RESULTS: The structure, mechanical, electronic, vibration, and hydrogen bonding properties of a novel high-energy and low-sensitivity 5, 5'-dinitroamino-3, 3'-azo-oxadiazole 4, 7-diaminopyridazino [4, 5-c] furoxan salt have been studied by density functional theory. The calculated vibrational properties show that the low-frequency mode is mainly contributed by the vibration of the -NO2 group, and the high-frequency mode is mainly contributed by the vibration of the -NH2 group and the N7-H3 bond which protonates the cation. In addition, it is analyzed that the first bond to break may be the N-NO2 bond. The calculated hydrogen bond properties indicate that the hydrogen bond between water molecules and cations is N7-H3 O5 (1.563 Å), which is the shortest hydrogen bond among all hydrogen bonds. The presence of this exceptionally short hydrogen bond renders the N7-H3 and H6-O5 bonds resistant to disruption at high frequencies, underscoring the pivotal role of hydrogen bonding in stabilizing the structure of energetic materials. Given the absence of experimental and theoretical data on the electronic, mechanical, and vibrational properties of the material thus far, our calculations offer valuable theoretical insights into the ionic salts of high energy and low sensitivity. COMPUTATIONAL METHODS: All calculations have been carried out based on density functional theory (DFT) and implemented in the CASTEP code. The mode-conserving pseudopotential is utilized to describe the plane wave expansion function, while the PBE functional within the generalized gradient approximation (GGA) is employed to characterize the exchange-correlation interaction. Additionally, dispersion correction is applied using Grimme's DFT-D method.
ABSTRACT
CONTEXT: Energy-containing materials such as explosives have attracted considerable interest recently. In the field of high-energy materials, tetrazine and its derivatives can largely meet the requirements of high nitrogen content and oxygen balance. Nitrogen-rich energetic salts are important research subjects. Nitrogen-rich salt of 3,6-dinitramino-1,2,4,5-tetrazine is a high-energy nitrogen-rich material, but there are few related studies. This paper systematically studies the crystal structure and electronic, vibrational, and thermodynamic properties of (NH4)2(DNAT). The lattice parameters of (NH4)2(DNAT) are observed to align well with the experimental values. The properties of electrons are analyzed by band structure and density of states (DOS). The phonon dispersion curves indicate that the compound is dynamically stable. The vibrational modes of bonds and chemical groups are described in detail, and the peaks in the Raman and infrared spectra are assigned to different vibration modes. Based on the vibration characteristics, thermodynamic properties such as enthalpy (H), Helmholtz free energy (F), entropy (S), Gibbs free energy (G), constant volume heat capacity (CV), and Debye temperature (Θ) are analyzed. This article can pave the way for subsequent work or provide data support to other researchers, promoting further research. METHODS: In this study, we utilized the density functional theory (DFT) for our calculations. The exchange-correlation potential and van der Waals interactions were characterized based on the GGA-PBE + G function calculation. We obtained Brillouin zone integrals using Monkhorst-Pack k-point grids, with the k-point of the Brillouin zone set to a 2 × 2 × 2 grid. During the self-consistent field operation, we set the total energy convergence tolerance to 5 × 10-6 eV per atom. The cut-off energy for the calculation was established at 830 eV. Additionally, the states of H (1s1), C (2s2 2p2), N (2s2 2p3), and O (2s2 2p4) were treated as valence electrons in our study.
ABSTRACT
CONTEXT AND RESULTS: The structure, electronic and optical properties of single-layer transition metallic chalcogenides ZrX3 (X = S, Se, Te) have been studied by density functional theory. The electron energy dispersion curve shows that ZrX3 has semiconductor properties, in which the conduction band is mainly contributed by the correlated states of the Zr-d orbital, and the valence band is mainly contributed by the correlated states of the X-p orbital. It is found that b-axis and biaxial strain have great influence on the bandgap and the shift of density of states is also large. At the same time, the peak value of density of states increases greatly when biaxial strain is applied. It is of guiding significance for selecting suitable substrates to prepare two-dimensional ZrX3 materials to study their electronic properties. The calculation of optical constants confirms that ZrX3 has strong optical anisotropy. In the visible range, the light absorption efficiency of ZrX3 in the direction of electric field polarization [100] is higher than that in the direction of [010]. The reflectance spectral results show that ZrS3 and ZrSe3 in the [100] directions have the highest reflectance, and ZrTe3 in the [010] direction has the highest reflectance, even in the long electromagnetic radiation range (up to 10 eV), which is of great significance for the construction of visible optical devices. COMPUTATIONAL METHOD: All computations have been carried out based on density functional theory (DFT) as implemented in the CASTEP code. The pseudo-potential is adopted by the norm conserving, and the exchange correlation functional is adopted by the Perdew-Burke-Ernzerhof in local generalized gradient approximation (GGA).
ABSTRACT
Hydrides of alkaline-earth and rare-earth metals have garnered significant interest in high-temperature superconductor research due to their excellent electron-phonon coupling and high T c upon pressurization. This study explores the electronic structures and electron-phonon coupling of metal hydrides XHn (n = 4,6), where X includes Ca, Mg, Sc, and Y. The involvement of d-orbital electrons alters the Fermi surface, leading to saddle-point nesting and a charge density wave (CDW) phase transition, which opens the superconducting gap. For instance, in YH6, the exchange coupling between Y-4d and H-1s holes in the phonon softening region results in T c values up to 230 K. The study suggests that factors, such as the origin of the CDW order, hydrogen concentration, and d-orbital contributions are crucial to superconductivity. This work proposes a new rule for high T c superconductors, emphasizing the importance of double gaps and electron-phonon interactions at exchange coupling sites, and predicts potential high-quality superconductors among rare-earth hydrides.
ABSTRACT
CONTEXT AND RESULTS : In this paper, the crystal structure, electronic, optical, and mechanical properties of SrVO3 have been systematically studied by first-principles calculation. The results show that the calculated lattice parameters are in good agreement with the experimental values of X-ray diffraction. The density of states is described in detail in this paper. By analyzing the crystal structure and electronic properties of SrVO3, the magnetic properties of SrVO3 are obtained from the one unpaired electrons of V and the exchange interaction between two V ions. At the same time, a detailed analysis of the optical properties of SrVO3 was conducted, and it was found that it is transparent in the visible light range. Finally, the mechanical properties of SrVO3 are calculated, which can provide some references for future research. COMPUTATIONAL METHOD: In this paper, a first-principles method based on density functional theory (DFT) is reported for PBE-GGA analysis using the plane wave-pseudo potential method in a quantum concentrate packet, U value of 7 eV to V-d and a U value of 2 eV to O-p, Grimme correction by DFT-D method. The k points in the Brillouin region are set to 4 × 4 × 4. The energy convergence criterion for self-consistent field calculation is set at 5.0 × 10-6 eV/atom, and the cutoff energy is 1170 eV. In this paper, the force acting on each atom is not more than 0.01 eV/Å, the maximum stress is not more than 0.02GPa, and the maximum atomic displacement is 5 × 10-4 Å.
ABSTRACT
CONTEXT: The key factor in designing heat-resistant energetic materials is their thermal sensitivity. Further research and prediction of thermal sensitivity remains a great challenge for us. This study is based on first-principles calculations and establishes a theoretical model, which comprehensively considers band gap, density of states, and Young's modulus to obtain a empirical parameter Ψ. A quantitative relationship was established between the new parameter and the thermal decomposition temperature. The value of Ψ is calculated for 10 energetic materials and is found to have a strong correlation with the experimental thermal decomposition temperature. This further proves the reliability of our model. Specifically, the larger the value of Ψ, the higher the thermal decomposition temperature, and the more stable the energetic material will be. Therefore, to some extent, we can use the new parameter Ψ calculated by the model to predict thermal sensitivity. METHODS: Based on first-principles, this paper used the Cambridge Serial Total Energy Package (CASTEP) module of Materials Studio (MS) for calculations. The Perdew-Burke-Ernzerhof (PBE) functionals in Generalized Gradient Approximation (GGA) method as well as the Grimme dispersion correction was used in this paper.
ABSTRACT
CONTEXT: The addition of central metal atoms to hydrogen clathrate structures is thought to provide a certain amount of "internal chemical pressure" to offset some of the external physical pressure required for compound stability. The size and valence of the central atoms significantly affect the minimum pressure required for the stabilization of hydrogen-rich compounds and their superconducting transition temperature. In recent years, many studies have calculated the minimum stable pressure and superconducting transition temperature of compounds with H24, H29, and H32 hydrogen clathrates, with centrally occupied metal atoms. In order to investigate the stability and physical properties of compounds with H cages in which the central atoms change in the same third group B, herein, based on first-principles calculations, we systematically investigated the lattice parameters, crystal volume, band structures, density of states, Mulliken analysis, charge density, charge density difference, and electronic localization function in I m 3 ¯ m -MH6 and P63/mmc-MH9 systems with different centered rare earth atoms M (M = Sc, Y, La) under a series of pressures. We find that for MH9, the pressure mainly changes the crystal lattice parameters along the c-axis, and the contributions of the different H atoms in MH9 to the Fermi level are H3 > H1 > H2. The density of states at the Fermi level of MH6 is mainly provided by H 1 s. Moreover, the size of the central atom M is particularly important for the stability of the crystal. By observing a series of properties of the structures with H24 and H29 cages wrapping the same family of central atoms under a series of pressures, our theoretical study is helpful for further understanding the formation mechanism of high-temperature superconductors and provides a reference for future research and design of high-temperature superconductors. METHODS: The first principles based on the density functional theory and density functional perturbation theory were employed to execute all calculations by using the CASTEP code in this work.
ABSTRACT
In order to investigate the impact of an external electric field on the sensitivity of ß-HMX explosives, we employ first-principles calculations to determine the molecular structure, dipole moment, and electronic properties of both ß-HMX crystals and individual ß-HMX molecules under varying electric fields. When the external electric field is increasing along the [100], [010], and [001] crystallographic directions of ß-HMX, the calculation results indicate that an increase in the bond length (N1-N3/N1'-N3') of the triggering bond, an increase in the main Qnitro (N3, N3') value, an increase in the minimum surface electrostatic potential, and a decrease in band gap all contribute to a reduction in its stability. Among these directions, the [010] direction exhibits the highest sensitivity, which can be attributed to the significantly smaller effective mass along the [010] direction compared with the [001] and [100] directions. Moreover, the application of an external electric field along the Y direction of the coordinate system on individual ß-HMX molecules reveals that the strong polarization effect induced by the electric field enhances the decomposition of the N1-N3 bonds. In addition, due to the periodic potential energy of ß-HXM crystal, the polarization effect of ß-HMX crystal caused by an external electric field is much smaller than that of a single ß-HXM molecule.
ABSTRACT
CONTEXT: In order to study the relationship between the sensitivity and pressure of energetic materials, six kinds of energetic materials were selected as the research object. The crystal structure, electronic, and phonon properties under hydrostatic pressure of 0 ~ 45 GPa were calculated by first principles. The calculation results show that the lattice parameters and band gap values of these six energetic materials decrease with the increase of pressure. The peak of the density of states decreases and moves to the low energy direction, and the electrons become more active. Meanwhile, the effect of pressure on the sensitivity of the energetic materials is analyzed based on the multi-phonon up-pumping theory. The number of doorway modes and integral of projected phonon density of states under high pressure is calculated. The results show that both of them increase with the increase of pressure. And the smaller the value of the band gap, the larger the number of doorway modes and integral of projected phonon density of states, and the more sensitive the energetic material is. METHODS: All calculations are performed using the Materials Studio software based on density functional theory. The Perdew-Burke-Ernzerhof (PBE) functional of the generalized gradient approximation (GGA) is used to calculate the exchange correlation function, and the Grimme dispersion correction method is used to deal with the weak intermolecular interaction. The structure of the compound was optimized by BFGS algorithm. The linear response is used to calculate the phonon properties of energetic materials. The plane wave cutoff energy was set to 830 eV. The K-point grids of TATB, FOX-7, TNX, RDX, TNT, and HMX were chosen as 2 × 2 × 2, 2 × 2 × 1, 2 × 1 × 1, 1 × 1 × 1, 1 × 2 × 1, and 2 × 1 × 2.
ABSTRACT
CONTEXT: Based on first principles, the structure, elasticity, mechanics, electronics, and optical properties of cubic K2Pb2O3 were studied. The structural parameters calculated by this method are close to the previous theoretical results. The elastic constant, bulk modulus, shear modulus, Young's modulus, Poisson's ratio, and mechanical stability are studied, and it is shown that cubic K2Pb2O3 is mechanically stable, isotropic, and brittleness. The electrical conductivity and chemical bonding of cubic K2Pb2O3 were analyzed based on the calculated band structure, density of states (DOS), and bond populations. The dispersion of optical functions, including the dielectric function, refractive index, extinction coefficient, reflectivity, absorption coefficient, and loss function, is displayed and analyzed. METHODS: All computations have been carried out based on density functional theory (DFT) as implemented in the CASTEP code. The norm conservation pseudopotential method is used to exchange correlation functionals within the generalized gradient approximation (GGA).
ABSTRACT
Predicting the melting temperature of materials has always been a topic of great concern. This article proposes an alternative model for determining the melting temperature of materials based on the main idea of the Lindemann melting criterion combined with the first-principles calculations of density functional theory. To verify the accuracy of the melting model, this article selected typical ionic crystals of MgO and 10 alkali metal halides as the validation objects. The calculation results indicate that the melting temperature of the MgO crystals and I-VII compounds is in good agreement with the experimental results.
ABSTRACT
INTRODUCTION: The electronic and optical properties of ß-Ga2O3 have been investigated by CASTEP using first principles. It is found that ß-Ga2O3 has an indirect band gap and the conduction band base is located at the Γ point. The stability of ß-Ga2O3 is demonstrated by the calculation of elastic constants, and the ductility of ß-Ga2O3 is demonstrated by the ratio of Poisson's ratio to shear modulus. The optical property analysis shows that ß-Ga2O3 has a high absorption capacity in the ultraviolet region, but a low absorption capacity in visible and infrared light. CONTEXT: The structure, optical, and electronic properties of ß-Ga2O3 are calculated and analyzed based on first-principles calculation. The optimized structures of ß-Ga2O3 are in good agreement with previously studied. In this paper, the elastic, electronic, and optical properties of ß-Ga2O3 are calculated. METHODS: The CASTEP code was employed to execute these calculations in the present work, where the exchange-correlation interactions were treated in the generalized gradient approximation (GGA) using the Perdew-Burke-Ernzerhof (PBE) functional in the geometry optimizations and electronic and elastic properties.
ABSTRACT
METHODS: This study used molecular dynamics (MD) to simulate three materials (HMX, FOX-7, and TATB) under the NVT ensemble. Six temperatures (100 K, 200 K, 300 K, 400 K, 500 K, and 600 K) were simulated. In addition, the trigger bond lengths, energy bands, and density of states of three materials were obtained at different temperatures and compared with the calculated results at 0 K. CONTEXT: The results indicate that the trigger bond lengths of the three materials are very close to the experimental values. Overall, the maximum and average bond lengths of the trigger bonds increase with increasing temperature. The band gap value decreases with increasing temperature. The changes in trigger bond length and band gap value are consistent with the experimental fact that sensitivity increases with increasing temperature. And Eg > 1 eV is consistently found within the temperature range of 0-600 K, indicating that all three materials are non-metallic compounds.
ABSTRACT
The determination of impact sensitivity of energetic materials traditionally relies on expensive and safety-challenged experimental means. This has instigated a shift towards scientific computations to gain insights into and predict the impact response of energetic materials. In this study, we refine the phonon-vibron coupling coefficients ζ in energetic materials subjected to impact loading, building upon the foundation of the phonon up-pumping model. Considering the full range of interactions between high-order phonon overtones and molecular vibrational frequencies, this is a pivotal element for accurately determining phonon-vibron coupling coefficients ζ. This new coupling coefficient ζ relies exclusively on phonon and molecular vibrational frequencies within the range of 0-700 cm-1. Following a regression analysis involving ζ and impact sensitivity (H50) of 45 molecular nitroexplosives, we reassessed the numerical values of damping factors, establishing a = 2.5 and b = 35. This coefficient is found to be a secondary factor in determining sensitivity, secondary to the rate of decomposition propagation and thermodynamic factor (heat of explosion). Furthermore, the relationship between phonon-vibron coupling coefficients ζ and impact sensitivity was studied in 16 energetic crystalline materials and eight nitrogen-rich energetic salts. It was observed that as the phonon-vibron coupling coefficient increases, the tendency for reduced impact sensitivity H50 still exists.
ABSTRACT
CONTEXT: 1,3,5-triamino-2,4,6-trinitrobenzene (TATB) is a typical insensitive energetic material. It can be used in explosive formulations, such as PBX-9502 and LX-17-0. TATB is an intriguing and unusual explosive for another reason: it crystallizes into a wide array of planar hydrogen bonds, forming a graphite-like layered structure. Therefore, TATB is one of the important research objects, and its surface structure needs to be deeply understood. In this research work, the electronic and energetic properties of TATB (001) surface are explored. METHODS: In this paper, the structural, electronic, energetic properties and impact sensitivity of TATB (001) surface structure at 0 and -3 GPa along with x-axis were calculated in this study using the first-principles calculations. The calculations in this paper are performed in the CASTEP code, which is based on the density functional theory with the first-principles calculation method using the plan-wave pseudopotential approach. The exchange-correlation interaction was adopted by the generalized gradient approximation (GGA) with the Perdew-Burke-Ernzerhof (PBE) functional. The DFT-D method with the Grimme correction accurately models van der Waals interactions. To model the surface structures of TATB, the planar slab method was employed. We constructed TATB (001) periodic slabs including three layers with a 15-Å vacuum layer.
ABSTRACT
CONTEXT: To find the potential value of Ca3N2 in the field of optoelectronics, the physical properties of Ca3N2 will be analyzed. It can be concluded from the electronic properties that the Ca-N bonds of α-Ca3N2 are more stable than those of δ-Ca3N2 and ε-Ca3N2. The dielectric function, reflectivity function, and absorption function of three types of Ca3N2 were accurately calculated, and it was concluded that α-Ca3N2, δ-Ca3N2, and ε-Ca3N2 have greater transmittance for visible light and exhibit optical transparency in the near-infrared frequency domain. Combined with the high hardness, strong bonding, high melting point, and wear resistance of Ca3N2, Ca3N2 can be used as a new generation of window heat-resistant materials. The α-Ca3N2, δ-Ca3N2, and ε-Ca3N2 are indirect, direct, and indirect narrow bandgap compounds, respectively, that is, δ-Ca3N2 is more suitable for luminescent materials than α-Ca3N2 and ε-Ca3N2. α-Ca3N2 and δ-Ca3N2 have high reflective properties in the ultraviolet region and can be used as UV protective coatings. All three Ca3N2 materials can be used industrially to synthesize photovoltaic devices that operate in the ultraviolet region. METHODS: Based on the first-principles of density functional theory calculations, the structures, electronic properties, and optical properties of α-Ca3N2, δ-Ca3N2, and ε-Ca3N2 were calculated. The calculation results show that although the α-Ca3N2, δ-Ca3N2, and ε-Ca3N2 have similar electronic structures, some phases have better properties in some aspects.
ABSTRACT
Aluminized explosive has attracted more and more attention in recent years because of its high explosive heat and high power. Al2O and AlO are indispensable aluminum oxides in the explosion process of aluminized explosives. The study of the physical properties of solid Al2O and AlO under pressure may play an important role in the understanding of the explosion mechanism of aluminized explosives. CONTEXT: The structures, cold-pressed lines and electronic properties of cubic Al2O and AlO are calculated and analyzed based on first-principles calculation in this paper. The optimized structures of Al2O and AlO are in good agreement with those previously studied. The cold pressure line shows that the specific volumes of Al2O and AlO decrease with increasing pressure. The peak values and peak positions of density of state of Al2O and AlO change greatly under pressure. METHODS: The CASTEP code was used to execute these calculations throughout the present work, where the plane-wave basis set and norm conserving pseudopotential were employed.
ABSTRACT
The sensitivity of energetic materials along different crystal directions is not the same and is anisotropic. In order to explore the difference in friction sensitivity of different surfaces, we calculated the structure, excess energy, surface energy, electronic structure, and the nitro group along (1 1 1), (1 1 0), (1 0 1), (0 1 1), (0 0 1), (0 1 0), and (1 0 0) surfaces of EDNA based on density functional theory. The analysis results showed that relative to other surfaces, the (0 0 1) surface has the shortest N-N average bond length, largest N-N average bond population, smallest excess energy and surface energy, widest band gap, and the largest nitro group charge value, which indicates that the (0 0 1) surface has the lowest friction sensitivity compared to other surfaces. Furthermore, the conclusions obtained by analyzing the excess energy are consistent with the results of the N-N bond length and bond population, band gap, and nitro charge. Therefore, we conclude that the friction sensitivity of different surfaces of EDNA can be evaluated using excess energy.
ABSTRACT
The metallization of alkaline earth metal hydrides offers a way to achieve near-room temperature superconductivity. In order to explore the metallization mechanism of these hydrides under pressure, a detailed understanding of the property changes of alkaline earth metal hydrides is required. Based on first-principles calculations, we have systematically investigated the dihydrides (XH2, X = Be, Mg, Ca, Sr, Ba) and tetrahydrides (XH4, X = Mg, Ca, Sr, Ba) of alkaline earth metals, respectively. By applying external pressure, we show that the structures of these alkaline earth metal hydrides undergo a series of phase transitions. Moreover, we investigate how the size of the bandgap decreases and eventually closes and reveal the role of electronegativity of metal elements in the critical pressure of hydride metallization. Remarkably, the hydrogen units (H6 or H8) formed in XH4 can accelerate the metallization process. The increase of the energy level difference in hydrogen units promotes the electroacoustic coupling effect, which is conducive to realization of high superconducting transition temperature (Tc). Our theoretical findings identify MgH4-I4/mmm as having potential to be a high-temperature superconductor and provide unusual ideas for the search of unknown high-temperature superconducting materials.
ABSTRACT
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.