Looking for high energy density compounds among polynitraminepurines

A series of purine derivatives with nitramine groups are calculated by using density functional theory (DFT). The molecular theory density, heats of formation, bond dissociation energies and detonation performance are investigated at DFT-B3LYP/6-311G** level. The isodesmic reaction method is employed to calculate the HOFs of the energies obtained from electronic structure calculations. Results show that the position of nitramine groups can influence the values of HOFs. The bond dissociation energies and the impact sensitivity are analyzed to investigate the thermal stability of the purine derivatives. The calculated bond dissociation energies of ring-NHNO₂ and NH-NO₂ bond show that the NH-NO₂ bond should be the trigger bond in pyrolysis processes. The H₅₀ of most compounds are larger than that of CL-20 and RDX.     

Computational study of energetic nitrogen-rich derivatives of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole

The heats of formation (HOFs), electronic structure, energetic properties, and thermal stabilities for a series of 1,4-bis(1-azo-2,4-dinitrobenzene)-iminotetrazole derivatives with different substituents and substitution positions and numbers of nitrogen atoms in the nitrobenzene rings were studied using the DFT-B3LYP method. All the substituted compounds have higher HOFs than their parent compounds. As the number of nitrogen atoms in the nitrobenzene ring increases, the HOFs of the derivatives with the same substituent rise gradually. Replacing carbon atoms in the nitrobenzene with nitrogen atoms to form N–N bonds is very helpful in improving their HOFs. Most of the substituted compounds have higher HOMO–LUMO gaps than the corresponding unsubstituted compounds. Substitution of the –NO₂, –NF₂, or –ONO₂ group and an increase in the number of nitrogen atoms in the nitrobenzene rings are useful for enhancing their detonation performance. The substituents’ substitution is not favorable for improving thermal stability. Considering detonation performance and thermal stability, five compounds may be considered potential candidates for high energy density compounds (HEDCs).     

Molecular design of aminopolynitroazole-based high-energy materials

The density functional theory (DFT) was employed to calculate the energetic properties of several aminopolynitroazoles. The calculations were performed to study the effect of amino and nitro substituents on the heats of formation, densities, detonation performances, thermal stabilities, and sensitivity characteristics of azoles. DFT-B3LYP, DFT-B3PW91, and MP2 methods utilizing the basis sets 6-31 G* and 6-311 G (2df, 3p) were adopted to predict HOFs via designed isodesmic reactions. All of the designed aminopolynitroazoles had heats of formation of >220 kJ mol(-1). The crystal densities of the aminopolynitroazoles were predicted with the cvff force field. All of the energetic azoles had densities of >1.83 g/cm(3). The detonation velocities and pressures were evaluated using the Kamlet-Jacobs equations, utilizing the predicted densities and heats of formation. It was found that aminopolynitroazoles have a detonation velocity of about 9.1 km/s and detonation pressure of 36 GPa. The bond dissociation energies for the C-NO(2) and N-NO(2) bonds were analyzed to investigate the stabilities of the designed molecules. The charge on the nitro group was used to assess impact sensitivity in the present study. The results obtained imply that the designed molecules are stable and are expected to be candidates for high-energy materials (HEMs).     

Theoretical investigation on the heats of formation and detonation performance in polydinitroaminocubanes

A series of polydinitroaminocubanes have been designed computationally. We calculated the heats of formation, the detonation velocity (D) and detonation pressure (P) of the title compounds by density function theory (DFT) with 6-311?G** basis set. The relationship between the heats of formation and the molecular structures is discussed. The result shows that all cubane derivatives have high and positive heats of formation, which increase with increasing number of dinitroamino groups. The detonation performances of the title compound were estimated by Kamlet-Jacobs equation, and the result indicated that most cubane derivatives have good detonation performance over RDX (hexahydro-1,3,5-trinitro-1,3,5-trizine) and HMX (1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane). In addition, we also found that the heat of detonation (Q) is another very important impact in increasing detonation performance except density. The relative stabilities of the title compound are discussed in the terms of the calculated heats of formation, and the energy gaps between the frontier orbitals. The results have not only shown that these compounds may be used as high energy density compounds (HEDCs), but also provide some useful information for further investigation.     

Density functional study on the derivatives of purine

The derivatives of purine are designed through substituting the hydrogen atoms on it for nitro and amino functional groups. Geometries and frequency are analyzed at the B3LYP/6-31 G** level of density functional theory(DFT). Heats of formation (HOF), bond dissociation energy(BDE) and detonation parameters (detonation velocity and detonation pressure) are obtained in detail at the same level. It is found that the BDE values of all derivatives are over 120KJ·mol(-1), and have high positive heats of formation. These derivatives possess excellent detonation properties, for B1, B2, and C, the detonation velocity are 9.58, 9.57,and 9.90 km·s(-1), and the detonation pressure are 43.40,46.05, and 46.37 Gpa, respectively, the detonation performances are better than cyclotrimethylenetrinitramine (RDX)and cyclotetramethylenetetranitramine (HMX). Hence, the derivations of purine may be promising well-behaved high energy density materials.     

Density functional theory calculations on the thermodynamic properties of polynitrosoprismanes

A series of polynitrosoprismanes, C(6)H(6 - n )(NO)( n ) (n?=?1-6), considered as high energy density compounds (HEDCs), have been designed computationally. We calculated the electronic structures, the heats of formation, the specific enthalpies of combustion, the bond dissociation energies, and the strain energies of the title compounds using density functional theory (DFT) with the 6-311G** basis set. It was found that the ΔE (LUMO-HOMO) values of the title compounds decrease as the number of nitroso groups increase, and the energy gaps of the prismane derivatives are much lower than that of TATB. Their high positive heats of formation indicate that polynitrosoprismanes can store a great deal of energy. Furthermore, the HOFs for the nitrosoprismane series were observed to decrease until three nitroso groups were connected to the prismane skeleton. For the polynitrosoprismanes, the trigger bond was confirmed to be the C-C bond in the skeleton. According to our calculations, all nitrosoprismanes appear to have large strain energies, and these calculations can provide basic information that may prove useful for the molecular design of novel high energy density materials.     

Dinitroamino benzene derivatives: a class new potential high energy density compounds

Dinitroamino benzene derivatives are designed and studied in detail with quantum chemistry method. The molecular theory density, heats of formation, bond dissociation energies, impact sensitive and detonation performance are investigated at DFT-B3LYP/6-311G** level. The results of detonation performance indicated most of the compounds have better detonation velocity and pressure than RDX and HMX. The N-N bond can be regard as the trigger bond in explosive reaction, and the bond dissociation energies of trigger bond are almost not affected by the position and number of substituent group. The impact sensitive are calculated by two different theory methods. It is found that the compounds, which can become candidates of high energy materials, have smaller H₅₀ values than RDX and HMX. It is hoped that this work can provide some basis information for further theory and experiment studies of benzene derivatives.     

Looking for high energy density compounds among polynitraminecubanes

Based on fully optimized geometric structures at DFT-B3LYP/6-311G** level, we calculated electronic structures, heats of formation, strain energies, bond dissociation energies and detonation performance (detonation velocity and detonation pressure) for a series of polynitraminecubanes. Our results have shown that energy gaps of cubane derivatives are much higher than that of triaminotrinitrobenzene (TATB), which means that cubane derivatives may be more sensitive than TATB. Polynitraminecubanes have high and positive heats of formation, and a good linear relationship between heats of formation and nitramine group numbers was presented. As the number of nitramine groups in the molecule increases, the enthalpies of combustion values are increasingly negative, but the specific enthalpy of combustion values decreases. It is found that all cubane derivatives have high strain energies, which are affected by the number and position of nitramine group. The calculated bond dissociation energies of C-NHNO₂ and C-C bond show that the C-C bond should be the trigger bond in the pyrolysis process. It is found that detonation velocity (D), detonation pressure (P) and molecule density (ρ) have good linear relationship with substituented group numbers. Heptanitraminecubane and octanitraminecubane have good detonation performance over 1,3,5,7-tetranitro-1,3,5,7-tetraazacyclooctane (HMX), and they can be regarded as potential candidates of high energy density compounds (HEDCs). The results have not only shown that these compounds may be used as HEDCs, but also provide some useful information for further investigation.     

Density functional calculations for a high energy density compound of formula C6H6?n (NO2) n

A series of polynitroprismanes, C(6)H(6-n )(NO(2))(n) (n?=?1-6) intended for use as high energy density compounds (HEDCs) were designed computationally. Their electronic structures, heats of formation, interactions between nitro groups, specific enthalpies of combustion, bond dissociation energies, and explosive performances (detonation velocities and detonation pressures) were calculated using density functional theory (DFT) with the 6-311 G** basis set. The results showed that all of the polynitroprismanes had high positive heats of formation that increased with the number of substitutions for the prismane derivatives, while the specific enthalpy of combustion decreased as the number of nitro groups increased. In addition, the range of enthalpy of combustion reducing is getting smaller. Interactions between ortho (vicinal) groups deviate from the group additivity rule and decrease as the number of nitro groups increases. In terms of thermodynamic stability, all of the polynitroprismanes had higher bond dissociation energies (BDEs) than RDX and HMX. Detonation velocities and detonation pressures were estimated using modified Kamlet-Jacobs equations based on the heat of detonation (Q) and the theoretical density of the molecule (ρ). It was found that ρ, D, and P are strongly linearly related to the number of nitro groups. Taking both their energetic properties and thermal stabilities into account, pentanitroprismane and hexanitroprismane are potential candidate HEDCs.     

Theoretical design of energetic nitrogen-rich derivatives of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane

The heats of formation (HOFs), energetic properties, and thermal stability of a series of 1,7-diamino-1,7-dinitrimino-2,4,6-trinitro-2,4,6-triazaheptane derivatives with different substituents, different numbers of substituents, and different original chains are found by using the DFT-B3LYP method. The results show that -NO₂ or -NH₂ is an effective substituent for increasing the gas-phase HOFs of the title compounds, especially -NO₂ group. As the numbers of substitutents increase, their HOFs enhance obviously. Increasing the length of original chain is helpful for improving their HOFs. The substitution of -NO₂ is useful for enhancing their detonation performances and the effects of the length of original chains on detonation properties are coupled with those of the substituents. An analysis of the BDE of the weakest bonds indicates that the substitution of the -NH₂ groups and replacing the -NO₂ groups of N-NO₂ by the -NH₂ groups are favorable for improving their thermal stability, while the substitution of -NO₂ and increasing the length of original chain decrease their thermal stability. Considering the detonation performance and thermal stability, seven compounds may be considered as the potential candidates of high energy density compounds.     

Variations of the tautomeric preferences and π-electron delocalization for the neutral and redox forms of purine when proceeding from the gas phase (DFT) to water (PCM)

Quantum-chemical calculations were performed for all possible nine neutral tautomers of purine and their oxidized and reduced forms in water {PCM//DFT(B3LYP)/6?311+G(d,p)} and compared to those in the gas phase {DFT(B3LYP)/6?311+G(d,p)}. PCM hydration influences geometries, π-electron delocalization, and relative energies of purine tautomers in different ways. Generally, the harmonic oscillator model of electron delocalization (HOMED) indices increase when proceeding from the gas phase to aequeous solution for the neutral and redox forms of purine. Their changes for the neutral and oxidized tautomers are almost parallel to the relative energies showing that aromaticity plays an important role in the tautomeric preferences. Tautomeric stabilities and tautomeric preferences vary when proceeding from the gas phase to water indicating additionally that intra- and intermolecular interactions affect tautomeric equilibria. The tautomeric mixture of neutral purine in the gas phase consists mainly of the N9H tautomer, whereas two tautomers (N9H and N7H) dominate in water. For oxidized purine, N9H is favored in the gas phase, whereas N1H in water. A gain of one electron dramatically changes the relative stabilities of the CH and NH tautomers that C6H and C8H dominate in the tautomeric mixture in the gas phase, whereas N3H in water. These variations show exceptional sensitivity of the tautomeric purine system on environment in the electron-transfer reactions.     

Theoretical studies on 2-diazo-4,6-dinitrophenol derivatives aimed at finding superior propellants

In an attempt to find superior propellants, 2-diazo-4,6-dinitrophenol (DDNP) and its –NO₂, –NH₂, –CN, –NC, –ONO₂, and –NF₂ 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₂ bond appears to be the trigger bond during the thermolysis process for these compounds, except for the –ONO₂ and –NF₂ 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₂ derivatives were higher than that of DDNP, and also meet the requirements of propellants.     

H-bonded complexes between acetylacetone and two molecules of methanol: HF and DFT level study

Five stable H-bonded complexes (supersystems) between acetylacetone and two methanol molecules were investigated at the B3LYP and HF levels of theory using the 6-311G** and 6-11++G** basis sets. The most stable complex was found as the one with the highest relative bonding and interaction energies. All vibrational frequencies resulting from calculations with the 6-311++G** basis set were compared with the recorded IR spectrum of acetylacetone/methanol mixture in a molar ratio 1:2.     

A B3LYP and MP2 theoretical investigation into host-guest interaction between calix[4]arene and Li+ or Na+

The DFT-B3LYP and MP2 methods with 6-311G** and 6-311++G** basis sets have been applied to study the complexation energies of the host-guest complexes between the cone calix[4]arene and Li⁺ or Na⁺ on the B3LYP optimized geometries. A comparison of the complexation energies obtained from the MP2(full) with those from MP2(fc) method is also carried out. The result shows that it is essential to introduce the diffuse basis set into the geometry optimizations and complexation energy calculations of the alkali-metal cation-π interaction complexes of calix[4]arene, and the D _e values show a maximum of 21.13 kJ mol⁻¹ (14.45% of relative error) between the MP2(full)/6-311++G** and MP2(fc)/6-311++G** method. For Li⁺ cation, the complexation is mainly energetically stabilized by the lower rim/cation (namely O–Li⁺) interaction. However, binding energies and NBO analyses confirm that Na⁺ cation prefers to enter the calix[4]arene cavity and the cation-π interaction is predominant, which contradicts the previous low-level theoretical studies. Furthermore, the complexation with Li⁺ is preferred over that with Na⁺ by at least 12.70 kJ mol⁻¹ at MP2(full)/6-311++G**//B3LYP/6-311++G** level.      

Theoretical potential scans and vibrational spectra of vinyl selenonyl halides CH2=CH–SeO2X (X is F, Cl and Br)

The structure and conformational stability of vinyl selenonyl fluoride, chloride and bromide CH₂=CH–SeO₂X (X is F, Cl and Br) were investigated using density functional B3LYP/6-311+G** and ab initio MP2/6-311+G** calculations. From the calculations the molecules were predicted to exist only in the non-planar gauche conformation with the vinyl C=C group almost eclipsing one of the selenonyl Se=O bonds as a result of conjugation between the two moieties. Single-minimum potential scans were calculated at the DFT level for the molecules. The vibrational frequencies were computed using B3LYP/6-311+G**. Normal coordinate calculations were then carried out and potential energy distributions were calculated for the three molecules in the gauche conformation.Figure Potential function for the asymmetric torsion in vinyl selenonyl fluoride (dotted line), chloride (dashed line) and bromide (solid line) as determined at the DFT-B3LYP/6-311+G** level     

Design and selection of nitrogen-rich bridged di-1,3,5-triazine derivatives with high energy and reduced sensitivity

The heats of formation (HOFs), electronic structures, energetic properties, and thermal stabilities of a series of energetic bridged di-1,3,5-triazine derivatives with different substituents and linkages were studied using density functional theory. It was found that the groups -N(3) and -N=N- are effective structural units for improving the HOF values of the di-1,3,5-triazine derivatives. The effects of the substituents on the HOMO-LUMO gap combine with those of the bridge groups. The calculated detonation velocities and detonation pressures indicate that substituting the -ONO(2), -NF(2), or -N=N- group is very useful for enhancing the detonation performance of these derivatives. Analysis of the bond dissociation energies for several relatively weak bonds suggests that most of the derivatives have good thermal stability. On the whole, the -NH(2), -N(3), -NH-, and -CH=CH- groups are effective structural units for increasing the thermal stabilities of the derivatives. Based on detonation performance and thermal stability, nine of the compounds can be considered potential candidates for high energy density materials with reduced sensitivity.     

Structural, IR, and EPR studies of the bis(methoxyacetato)diaquo-copper(II) complex

The structure, stability, and the IR, and EPR spectroscopic properties of bis(methoxyacetato)diaquo-copper(II) were studied both experimentally using FT-IR and theoretically using B3LYP/6-31G**, B3LYP/6-311G, BWP91/6-31G** methods. The same approaches were used to calculate the harmonic frequencies and to compare them to the experimental solid state values. The g-tensors are calculated using the NMR/GIAO computational method.     

Is 1-nitro-1-triazene a high energy density material?

An azo bridge (–N?=?N–) can not only desensitize explosives but also dramatically increase their heats of formation and explosive properties. Amino and nitro are two important high energy density functional groups. Here, we present calculations on 1-nitro-1-triazene (NH₂–N?=?N–NO₂). Thermal stability and detonation parameters were predicted theoretically at CCSD(T)/6-311G* level, based on the geometries optimized at MP2/6-311G* level. It was found that the p?→?π conjugation interaction and the intramolecular hydrogen bonding that exist in the system together increase the thermal stability of the molecule. Moreover, the detonation parameters were evaluated to be better than those of the famous HMX and RDX. Finally, the compound was demonstrated to be a high energy density material.     

Quantum mechanical and NMR spectroscopy studies on the conformations of the hydroxymethyl and methoxymethyl groups in aldohexosides

The potential energy surfaces of the hydroxymethyl and methoxymethyl groups in methyl hexopyranosides have been extensively studied, employing quantum mechanical calculations and high resolution NMR data. The structure and energy of the C-5-C-6 rotamers were calculated at the B3LYP level of the density functional theory (DFT). For all, geometry optimizations were carried out for 264 conformers of 16 methyl D-gluco- and methyl D-galactopyranoside derivatives 1-16 at the B3LYP/6-31G** level. For all calculated minima, single-point calculations were performed at the B3LYP/6-311++G** level. Solvent effects were considered using a self-consistent reaction field method. Values of the vicinal coupling constants 3J(H-5-H-6R), 3J(H-5-H-6S), 3J(C-4-H-6R), and 3J(C-4-H-6S) for methyl D-glucopyranosides, methyl D-galactopyranosides and their 6-O-methyl derivatives 9-16 were measured in two solvents, methanol and water. The calculated gg, gt, and tg rotamer populations of the hydroxymethyl and methoxymethyl groups in 9-16 agreed well with experimental data. The results clearly showed that the population of gg, gt, and tg rotamers is sensitive to solvent effects. It was concluded that the preference of rotamers in 1-16 is due to the hydrogen bonding and solvent effects.     

Computational Study of the Sulfonylated Amino Acid Hydroxamates Binding to the Zinc Ion within the Active Site of Carbonic Anhydrase

Abstract

The B3LYP/6–311+G(d,p) method and three ONIOM extrapolation methods ONI-OM (B3LYP/6–311+G(d,p): AM1); ONIOM(B3LYP/6–311+G(d,p): MNDO); ONIOM (B3LYP/6–311+G(d,p): HF/3-21G(d)) were used to characterize the complexes of Zn²⁺ cation with anionic sulfonylated amino acid hydroxamates (RSO₂NH-AA-CON(-)OH), possessing an unsubstituted RSO₂NH—amino acyl moiety. According to the R moiety we distinguish between pentafluorophenyl and 4-methoxyphenyl derivates. The amino acid hydroxamates included in the study were the Gly, Ala, and Leu derivates. Of the inhibitors investigated, the weakest zinc affinity exhibits the pentafluorophenyl derivate with Gly amino acid and the strongest affinity the 4-methoxyphenyl derivate with Leu amino acid. The inhibitors form bidentate coordination bonds with the zinc cation by means of the sulfonyl oxygen and the ionized hydroxamate nitrogen atoms, respectively. The zinc affinities computed using the B3LYP/6–311 +G(d,p)//HF/6–31 +G(d,p) method are in very good agreement with the full density functional theory (DFT) B3LYP/6–311+G(d,p)//B3LYP/6- 311+G(d,p) method and this method can be adopted to model larger complexes of inhibitors with the active site of carbonic anhydrase.