Search

The intramolecular hydrogen-bonding interactions and properties of a series of nitroamino[1,3,5]triazine-based guanidinium salts were studied by using the dispersion-corrected density functional theory method (DFT-D). Results show that there are evident LP(N or O; LP = lone pair)→σ*(N-H) orbital interactions related to O⋅⋅⋅H-N or N⋅⋅⋅H-N hydrogen bonds. Quantum theory of atoms in molecules (QTAIM) was applied to characterize the intramolecular hydrogen bonds. For the guanidinium salts studied, the intramolecular hydrogen bonds are associated with a seven- or eight-membered pseudo-ring. The guanylurea cation is more helpful for improving the thermal stabilities of the ionic salts than other guanidinium cations. The contributions of different substituents on the triazine ring to the thermal stability increase in the order of -NO(2)

Periodic density functional theory calculations are performed to study the hydrostatic compression effects on the structure, electronic, and thermodynamic properties of the energetic polyazide 4,4',6,6'-tetra(azido)hydrazo-1,3,5-triazine (TAHT) in the range of 0-100 GPa. At the ambient pressure, the local density approximation/Ceperley-Alder exchange-correlation potential parameterized by Perdew and Zunger relaxed crystal structure compares well with the experimental results. The predicted heat of sublimation is 38.68 kcal/mol, and the evaluated condensed phase of formation (414.04 kcal/mol) approximates to the experimental value. The detonation velocity and detonation pressure for the solid TAHT are calculated to be 7.44 km/s and 23.71 GPa, respectively. When the pressure is exerted less than 35 GPa, the crystal structure and geometric parameters change slightly. However, at 36 GPa, the molecular structure, band structure, and density of states change abnormally because of the azide-tetrazole transformation that has not been observed in gas phase or polar solvents. The azido group cyclizes to form a five-membered tetrazole ring that is coplanar with the riazine ring and contributes to a larger conjunction system. As the pressure augments further to 80 GPa, the hydrogen transfer is found and a new covalent bond H2-N9 is formed. In the studied pressure range, the band gap decreases generally except for some breaks due to the molecular transformation and drops to nearly zero at 100 GPa, which means the electronic character of the crystal changes toward a metallic system. An analysis of the electronic structure shows that an applied pressure increases the impact sensitivity of TAHT.

Density function theory has been employed to study a series of compounds containing pyridine ring at the B3LYP/6-31G* level. Detonation performance was evaluated by using the Kamlet–Jacobs equations based on the calculated densities and heats of formation. Some compounds have high densities (ca. 1.9 g cm−3) and good performance (detonation velocities over 9 km s−1, detonation pressures about 39 GPa) and may be the potential candidates of high energy density materials. The thermal stability and the pyrolysis mechanism of the title compounds were investigated by the bond dissociation energies and the impact sensitivity predicted. Solvent effect has been investigated and it makes the title compounds more stable in solutions.

A new polynitro cage compound 2, 4, 6, 8, 10, 12, 13, 14, 15-nonanitro-2, 4, 6, 8, 10, 12, 13, 14, 15-nonaazaheptcyclo [5.5.1.1(3,11).1(5,9)] pentadecane (NNNAHP) was designed in the present work. Its molecular structure was optimized at the B3LYP/6-31 G(d,p) level of density functional theory (DFT) and crystal structure was predicted using the Compass and Dreiding force fields and refined by DFT GGA-RPBE method. The obtained crystal structure of NNNAHP belongs to the P-1 space group and the lattice parameters are a = 9.99 Å, b = 10.78 Å, c = 9.99 Å, α = 90.01°, β = 120.01°, γ = 90.00°, and Z = 2, respectively. Based on the optimized crystal structure, the band gap, density of state, thermodynamic properties, infrared spectrum, strain energy, detonation characteristics, and thermal stability were predicted. Calculation results show that NNNAHP has detonation properties close to those of CL-20 and is a high energy density compound with moderate stability.

Ab initio calculations have been performed using the complete basis set model (CBS-QB3) to study the reaction mechanism of butane radical (C(4)H(9)•) with oxygen (O(2)). On the calculated potential energy surface, the addition of O(2) to C(4)H(9)• forms three intermediates barrierlessly, which can undergo subsequent isomerization or decomposition reaction leading to various products: HOO• + C(4)H(8), C(2)H(5)• + CH(2)CHOOH, OH• + C(3)H(7)CHO, OH• + cycle-C(4)H(8)O, CH(3)• + CH(3)CHCHOOH, CH(2)OOH• + C(3)H(6). Five pathways are supposed in this study. After taking into account the reaction barrier and enthalpy, the most possible reaction pathway is C(4)H(9)• + O(2) → IM1 → TS5 → IM3 → TS6 → IM4 → TS7 → OH• + cycle-C(4)H(8)O.

A series of azidamines were studied at the B3LYP/6-311+G(2df) level of density functional theory. Thermodynamic properties were calculated and increased quantitatively with the increasing temperature as well as the number of azido groups. The azidamines are highly energetic with large enthalpies of formation. The detonation performances of the azidamines were evaluated and their performances are comparable to those of hexahydro-1,3,5-trinitro-1,3,5-trizine and 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane. However, they are sensitive to impact as well as halogen azides according to the small potential energy barriers.

In an attempt to find superior propellants, 2-diazo-4,6-dinitrophenol (DDNP) and its -NO(2), -NH(2), -CN, -NC, -ONO(2), and -NF(2) derivatives were studied at the B3LYP/6-311++G level of density functional theory (DFT). Sensitivity was evaluated using bond dissociation enthalpies (BDEs) and molecular surface electrostatic potentials. The C-NO(2) bond appears to be the trigger bond during the thermolysis process for these compounds, except for the -ONO(2) and -NF(2) derivatives. Electrostatic potential results show that electron-withdrawing substituents make the charge imbalance more anomalous, which may change the strength of the bond, especially the weakest trigger bond. Most of the DDNP derivatives have the impact sensitivities that are higher than that of DDNP, making them favorable for use as solid propellants in micro-rockets. The theoretical densities (ρ), heats of formation (HOFs), detonation energies (Q), detonation pressures (P), and detonation velocities (D) of the compounds were estimated. The effects of various substituent groups on ρ, HOF, Q, D, and P were investigated. Some derivatives exhibit perfect detonation properties. The calculated relative specific impulses (I (r,sp)) of all compounds except for -NH(2) derivatives were higher than that of DDNP, and also meet the requirements of propellants.

Thirty nitro aromatic compounds are studied using the density functional theory B3LYP method with six basis sets (3-21G, 6-31G, 6-31G∗, 6-31G∗∗, 6-311G∗ and 6-31+G∗∗) and semiempirical PM3 method. Based on the geometries optimized at various theoretical levels, the molecular volumes and densities are calculated. Compared with the experimental results, the densities estimated by the PM3 and B3LYP/3-21G methods are all systematically larger, and that obtained with the other five basis sets are better and quite accurate. Considering that a larger basis set demands more computer resource, B3LYP/6-31G or 6-31G∗ is recommended for rapid and reliable prediction of crystalline densities. The effect of various groups on the densities is discussed. Based on the theoretical densities, it is possible to estimate the detonation properties, which will be helpful to rapidly and effectively design and screen promising candidates of high energy density materials (HEDMs).

The nitro derivatives of toluenes are optimized to obtain their molecular geometries and electronic structures at the DFT-B3LYP/6-31G* level. Detonation properties are evaluated using the modified Kamlet–Jacobs equations based on the calculated densities and heats of formation. It is found that there are good linear relationships between density, detonation velocity, detonation pressure and the number of nitro and methyl groups. Thermal stability and the pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies at the unrestricted B3LYP/6-31G* level. The activation energies of H-transfer reaction are smaller than the BDEs of all bonds and this illustrates that the pyrolysis of the title compounds may be started from the isomerization reaction of H transfer. According to the quantitative standard of energetics and stability as an HEDC (high energy density compound), pentanitrotoluene essentially satisfies this requirement. In addition, we have discussed the effect of the nitro and methyl groups on the static electronic structural parameters and the kinetic parameter.

Studies have suggested that octanitrocubane (ONC) is one of the most powerful non-nuclear high energy density material (HEDM) currently known. 2,4,6,8-Tetranitro-1,3,5,7-tetraazacubane (TNTAC) studied in this work may also be a novel HEDM due to its high nitrogen content and crystal density. Density functional theory and molecular mechanics methods have been employed to study the crystal structure, IR spectrum, electronic structure, thermodynamic properties, gas-phase and condensed-phase heat of formation, detonation performance, and pyrolysis mechanism of TNTAC. The TNTAC has a predicted density of about 2.12 g/cm(3), and its detonation velocity (10.42 km/s) and detonation pressure (52.82 GPa) are higher than that of ONC. The crystalline packing is P2(1)2(1)2(1), and the corresponding cell parameters are Z = 4, a = 8.87 Å, b = 8.87 Å, and c = 11.47 Å. Both the density of states of the predicted crystal and the bond dissociation energy of the molecule in gas phase show that the cage C-N bond is the trigger bond during thermolysis. The activation energy of the pyrolysis initiation reaction obtained from the B3LYP/6-311++G(2df,2p) level is 125.98 kJ/mol, which indicates that TNTAC meets the thermal stability request as an exploitable HEDM.

Two new nitramine compounds containing pyridine, 1,3,5,7-tetranitro-8-(nitromethyl) -4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine and its N-oxide 1,3,5,7-tetranitro-8- (nitromethyl)-4-imidazolino[4,5-b]4-imidazolino-[4,5-e]pyridine-4-ol were proposed. Density functional theory (DFT) has been employed to study the molecular geometries, electronic structures, infrared spectra, and thermodynamic properties at the B3LYP/6-31G* level. Their detonation performances evaluated using the Kamlet-Jacobs equations with the calculated densities and heats of formation are superior to those of HMX. The predicted densities of them were ca. 2 g cm(-3), detonation velocities were over 9 km s(-1), and detonation pressures were about 40 GPa, showing that they may be potential candidates of high energy density materials (HEDMs). The natural bond orbital analysis indicated that N-NO(2) bond is the trigger bond during thermolysis process. The stability of the title compounds is slightly lower than that of 2,4,6,8,10,12-hexanitro-2,4,6,8,10,12- hexaazaisowurtzitane (CL-20). The results of this study may provide basic information for the molecular design of new HEDMs.

A new polynitro cage compound with the framework of HNIW and a tetrazole unit, i.e., 10-(1-nitro-1, 2, 3, 4-tetraazol-5-yl)) methyl-2, 4, 6, 8, 12-hexanitrohexaazaisowurtzitane (NTz-HNIW) has been proposed and studied by density functional theory (DFT) and molecular mechanics methods. Properties such as IR spectrum, heat of formation, thermodynamic properties, and crystal structure were predicted. The compound belongs to the Pbca space group, with the lattice parameters a = 15.07 Å, b = 12.56 Å, c = 18.34 Å, Z = 8, and ρ = 1.990 g·cm(-3). The stability of the compound was evaluated by the bond dissociation energies and results showed that the first step of pyrolysis is the rupture of the N-NO(2) bond in the side chain. The detonation properties were estimated by the Kamlet-Jacobs equations based on the calculated crystal density and heat of formation, and the results were 9.240 km·s(-1) for detonation velocity and 40.136 GPa for detonation pressure. The designed compound has high thermal stability and good detonation properties and is probably a promising high energy density compound (HEDC).

Aliphatic polynitrates are studied using the density functional theory B3LYP method with basis set 6-31G*. The assigned infrared spectrum is obtained and is used to compute the thermodynamic properties based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics. On comparison of the theoretical densities with the experimental ones, the reliability of this theoretical method is tested. Detonation properties are evaluated using the modified Kamlet-Jacobs equations based on the calculated densities and heats of formation. According to the largest exothermic principle, the relative specific impulse (Is) is investigated by calculating the enthalpy of combustion (ΔH comb ) and the total heat capacity (C p,gases ). It is found that the introduction of methylene nitrate group could decrease the specific impulses on whole. Moreover, in combination with the energetic properties, xylitol pentanitrate, mannitol hexanitrate, volemitol heptanitrate, and 1,2,3,4,5,6,7,8-octanitrate n-octane are potential candidates for high energy density compounds.

The derivatives of DPO (2,5-dipicryl-1,3,4-oxadiazole) are optimized to obtain their molecular geometries and electronic structures at the DFT-B3LYP/6-31G* level. The bond length is focused to primarily predict thermal stability and the pyrolysis mechanism of the title compounds. Detonation properties are evaluated using the modified Kamlet-Jacobs equations based on the calculated densities and heats of formation. It is found that there are good linear relationships between density, detonation velocity, detonation pressure, and the number of azido, nitrate, and nitramine groups. According to the largest exothermic principle, the relative specific impulse is investigated by calculating the enthalpy of combustion (ΔH(comb)) and the total heat capacity (C(p,gases)). It is found that the introduction of -N(3), -ONO(2), and -NNO(2) groups could increase the specific impulses and II-4, II-5, and III-5 are potential candidates for High Energy Density Materials (HEDMs). The effect of the azido, nitrate, and nitramine groups on the structure and the properties is discussed.