A theoretical investigation on the structures,densities, detonation properties,and pyrolysis mechanism of the nitro derivatives of phenols

The nitro derivatives of phenols 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 hydroxy groups. Thermal stability and pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDEs) at the unrestricted B3LYP/6‐31G* level. The activation energies of H‐transfer reaction is smaller than the BDEs of all bonds and this illustrates that the pyrolysis of the title compounds may be started from breaking O? H bond followed by the isomerization reaction of H transfer. Moreover, the C? NO₂ bond with the smaller bond overlap population and the smaller BDE will also overlap may be before homolysis. According to the quantitative standard of energetics and stability as a high‐energy density compound, pentanitrophenol essentially satisfies this requirement. In addition, we have discussed the effect of the nitro and hydroxy groups on the static electronic structural parameters and the kinetic parameter. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

Theoretical Studies on the Structures,Thermodynamic Properties,Detonation Performance,and Pyrolysis Mechanisms for Six Dinitrate Esters

Density functional theory (DFT) has been employed to study the geometric and electronic structures of six dinitrate esters including ethylene glycol dinitrate (EGDN), diethylene glycol dinitrate (Di‐EGDN), triethylene glycol dinitrate (Tri‐EGDN), tetraethylene glycol dinitrate (Tetra‐EGDN), pentaethylene glycol dinitrate (Penta‐EGDN) and hexaethylene glycol dinitrate (Hexa‐EGDN) at the B3LYP/6‐31G* level. Their IR spectra were obtained and assigned by vibrational analysis. Based on the frequencies scaled by 0.96 and the principle of statistic thermodynamics, the thermodynamic properties were evaluated, which were linearly related with the number of CH₂CH₂O groups as well as the temperature, obviously showing good group additivity. Detonation performances were evaluated by the Kamlet‐Jacobs equations based on the calculated densities and heats of formation. It was found that density, detonation velocity, detonation pressure decreased with the increase of the number of CH₂CH₂O groups. Thermal stability and the pyrolysis mechanism of the title compounds were investigated by calculating the bond dissociation energies (BDE) at the B3LYP/6‐31G* level. For the nitrate esters, the O‐NO₂ bond is a trigger bond during a thermolysis initiation process.     

Theoretical Study of Heterocyclic Nitro Compounds for Use as Novel Explosives and Propellants

The heats of formation (HOFs) of heterocyclic nitro compounds were obtained by using a density functional theory B3LYP method with 6‐31G* and 6‐311+G** basis sets. The isodesmic reactions designed for the evaluation of HOFs keep most of the basic ring structures of the title compounds and thus ensure the credibility of the results. The values of HOFs are 567.90, 874.29 and 975.83 kJ/mol at the B3LYP/6‐31G* level for hexanitrohexazaadamantane ( A ), nonanitrononaza‐tetracyclo[7.3.1.1^3,7.1^5,11] pentadecane ( B ) and tetranitrotetrazacubane ( C ) respectively. The predicted detonation velocities of the title compounds are larger than, and detonation pressures are much larger than that of the widely used 1,3,5,7‐tetranitro‐1,3,5,7‐tetraazacyclooctane (HMX). The dissociation energy for the weakest C‐N bonds in the cage skeleton of the title compounds are 137‐144 kJ/mol at the B3LYP/6‐31G* level.     

Theoretical investigation on the structures,densities, and detonation properties of polynitrotetraazaoctahydroanthracenes

The polynitrotetraazaoctahydroanthracenes were optimized to obtain their molecular geometries and electronic structures at density functional theory–B3LYP/6‐31+G(d) level. Detonation velocities (D) and detonation pressures (P) were estimated for this nitramine compounds using Kamlet‐Jacobs equations, based on the theoretical densities (ρ) and heats of formation. It is found that there are good linear relationships between volume, density, detonation velocity, detonation pressure and the number of nitro group. Thermal stability of the compounds was investigated by calculating the bond dissociation energies and energy gap (ΔE_LUMO–HOMO). The simulation results reveal that molecule H performs similarly to famous explosive RDX. These results provide basic information for molecular design of novel high energetic density compounds. © 2011 Wiley Periodicals, Inc.     

Comparative studies on the structures,infrared spectrum,and thermodynamic properties of phthalocyanine using ab initio Hartree–Fock and density functional theory methods

The molecular structure, vibrational spectrum, standard thermodynamic functions, and enthalpy of formation of free base phthalocyanine (Pc) have been studied using the density functional theory B3LYP procedure, as well as the ab initio Hartree–Fock method. Various basis sets 3‐21G, 6‐31G*, and LANL2DZ have been employed. The results obtained at various levels are discussed and compared with each other and with the available experimental data. It is shown that calculations performed at the Hartree–Fock level cannot produce a reliable geometry and related properties such as the dipole moment of Pc and similar porphyrin‐based systems. Electron correlation must be included in the calculations. The basis set has comparatively less effect on the calculated results. The results derived at the B3LYP level using the smaller 3‐21G and LANL2DZ basis sets are very close to those produced using the medium 6‐31G* basis set. The geometry of Pc obtained at the B3LYP level has D_2h symmetry and the diameter of the central macrocycle is about 4 Å. The enthalpy of formation of Pc in the gas phase has been predicted to be 1518.50 kJ/mol at the B3LYP/6‐311G(2d,2p)//B3LYP/6‐31G* level via an isodesmic reaction. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Kinetic Study and NBO Analysis of the Dehydrogenation Mechanism of Five‐membered Ring Heterocyclic 2,5‐Dihydro‐[furan,thiophene, selenophene

The theoretical study of the dehydrogenation of 2,5‐dihydro‐[furan ( 1 ), thiophene ( 2 ), and selenophene ( 3 )] was carried out using ab initio molecular orbital (MO) and density functional theory (DFT) methods at the B3LYP/6‐311G**//B3LYP/6‐311G** and MP2/6‐311G**//B3LYP/6‐311G** levels of theory. Among the used methods in this study, the obtained results show that B3LYP/6‐311G** method is in good agreement with the available experimental values. Based on the optimized ground state geometries using B3LYP/6‐311G** method, the natural bond orbital (NBO) analysis of donor‐acceptor (bond‐antibond) interactions revealed that the stabilization energies associated with the electronic delocalization from non‐bonding lone‐pair orbitals [LP(e)_X3] to δ*_{C(1)  H(2)} antibonding orbital, decrease from compounds 1 to 3 . The LP(e)_X3→δ*_{C(1)  H(2)} resonance energies for compounds 1 – 3 are 23.37, 16.05 and 12.46 kJ/mol, respectively. Also, the LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the decrease of occupancies of LP(e)_X3 non‐bonding orbitals in ring of compounds 1 – 3 ( 3 > 2 > 1 ). The electronic delocalization from LP(e)_X3 non‐bonding orbitals to δ*_{C(1)  H(2)} antibonding orbital increases the ground state structure stability, Therefore, the decrease of LP(e)_X3→δ*_{C(1)  H(2)} delocalizations could fairly explain the kinetic of the dehydrogenation reactions of compounds 1 – 3 (k ₁ >k ₂ >k ₃ ). Also, the donor‐acceptor interactions, as obtained from NBO analysis, revealed that the (_C(4)C(7)→δ*_{C(1)  H(2)} resonance energies decrease from compounds 1 to 3 . Further, the results showed that the energy gaps between (_C(4)C(7) bonding and δ*_{C(1)  H(2)} antibonding orbitals decrease from compounds 1 to 3 . The results suggest also that in compounds 1 – 3 , the hydrogen eliminations are controlled by LP(e)→δ* resonance energies. Analysis of bond order, natural bond orbital charges, bond indexes, synchronicity parameters, and IRC calculations indicate that these reactions are occurring through a concerted and synchronous six‐membered cyclic transition state type of mechanism.     

Chemical Reactivity of Oxo‐ and Aza‐γ‐lactam Rings

Different mechanisms for the alkaline hydrolysis of oxo and aza‐γ‐lactam rings have been studied by ab initio calculations at the MP2/6‐31+G*//MP2/6‐31+G* and B3LYP/6‐31+G*//B3LYP/6‐31+G* levels. The tetrahedral intermediate can undergo two different reactions, the cleavage of the C₂−N₂ bond (the classical mechanism) and the cleavage of the C₂−X₆ bond (X=O, N). Both compounds present similar energy barriers for the classical fragmentation, and show considerably lower barriers for the alternative mechanism. Because of this reactivity, the compounds studied are expected to be β‐lactamase inhibitors.     

Theoretical studies on heats of formation for cubylnitrates using density functional theory B3LYP method and semiempirical MO methods

The heats of formation (HOF) have been calculated for all the 21 cubylnitrate compounds using the semiemprical molecular orbital (MO) methods (MINDO/3, MNDO, AM1, and PM3) and for 8 of 21 cubylnitrates containing 1–4 ? ONO₂ groups using the density functional theory (DFT) method at the B3LYP/6‐31G* level by means of designed isodesmic reactions. The cubane cage skeletons in cubylnitrate molecules have been kept in setting up isodesmic reactions to produce more accurate and reliable results. It is found that there are good linear relationships between the HOFs of the 8 cubylnitrates calculated using B3LYP/6‐31G* and two semiempirical MO (PM3 and AM1) methods, and the linear correlation coefficients of PM3 and AM1 methods are 0.9901 and 0.9826, respectively. Subsequently, the accurate HOFs at B3LYP/6‐31G* level of other 13 cubylnitrates containing 4–8 ? ONO₂ groups are obtained by systematically correcting their PM3‐calculated HOFs. Compared with noncaged nitrates, all the 21 cubylnitrates have high heats of formation implying that they may be very powerful energetic materials and have highly exploitable value. The relationship between the HOFs and the molecular structures of cubylnitrates has been discussed. © 2001 John Wiley & Sons, Inc. Int J Quantum Chem, 2001     

Thermochemical properties for n‐propyl,iso‐propyl,and tert‐butyl nitroalkanes,alkyl nitrites,and their carbon‐centered radicals

Density functional theory (DFT) based calculations are performed on a series of alkyl nitrites and nitroalkanes representing large‐scale primary, secondary, and tertiary nitro compounds and their radicals resulting from the loss of their skeletal hydrogen atoms. Geometries, vibration frequencies, and thermochemical properties [S°(T) and C°_p(T) (10 K ? T ? 5000 K)] are calculated at the B3LYP/6‐31G(d,p) DFT level. Δ_fH°₂₉₈ values are from B3LYP/6‐31G(d,p), B3LYP/6‐31+G(2d,2p), and the composite CBS‐QB3 levels. Potential energy barriers for the internal rotations have been computed at the B3LYP/6‐31G(d,p) level of theory, and the lower barrier contributions are incorporated into entropy and heat capacity data. The standard enthalpies of formation at 298 K are evaluated using isodesmic reaction schemes with several work reactions for each species. Recommended values derived from the most stable conformers of respective nitro‐ and nitrite isomers include ?30.57 and ?28.44 kcal mol^?1 for n‐propane‐, ?33.89 and ?32.32 kcal mol^?1 for iso‐propane‐, ?42.78 and ?41.36 kcal mol^?1 for tert‐butane‐nitro compounds and nitrites, respectively. Entropy and heat capacity values are also reported for the lower homologues: nitromethane, nitroethane, and corresponding nitrites. © 2010 Wiley Periodicals, Inc. Int J Chem Kinet 42: 181–199, 2010     

Theoretical Studies on Energetic Property and Stability of Pyrazine and Pyridine Derivatives

Density functional calculations at the B3LYP level with 6‐311G** and aug‐cc‐pVDZ basis sets were performed to predict the heats of formation (HOFs) for two pyrazine derivatives and eight pyridine derivatives. In the isodesmic reactions designed for the computation of heats of formation (HOFs), pyrazine and pyridine were chosen as reference compounds. The N‐oxidations for the ring nitrogen of pyrazine and pyridine derivatives decrease the HOF values when N‐oxide oxygen is neighboring with amino groups, but increase when it neighbors with nitro groups. Thermal stability was evaluated via bond dissociation energies (BDE) at the UB3LYP/6‐311G** level. As a whole, the homolysis of C–NO₂ bonds is the main step for bond dissociation of the title compounds. The BDE values of title compounds are influenced by intramolecular hydrogen bonds. The hydrogen bond effects associated with the length of the H···O bonds were analyzed by the electron density at the critical points and natural bond orbital.     

Looking for high‐energy density compounds among hexaazaadamantane derivatives with ?CN, ?NC,and ?ONO2 groups

Density functional theory (DFT) was employed to evaluate the heats of formation (HOFs) for hexaazaadamantane (HAA) derivatives with ? CN, ? NC, and ? ONO₂ groups, respectively. This was done by designing isodesmic reactions at the B3LYP/6‐31G* level of theory, where the HAA cage skeletons were kept unbroken to produce more accurate results, and all HOFs for the required reference compounds, NH₂CN, NH₂NC, NH₂ONO₂, and (CH₂NH)₃, were derived from the G₃ theory calculation based on the atomization energies. The calculation results show that the HOFs of HAA derivatives are mainly affected by the number and the position of substituent groups, all the obtained HOFs are positive, and the ? NC derivatives have the most HOFs among the three types of derivatives with the same number of substituent groups. The detonation velocity (D) and detonation pressure (P) were obtained from the empirical Kamlet–Jacobs equations. All the ? NC and ? CN derivatives of HAA have lower densities (ρ), heats of explosion (Q), D, and P. However, these properties of ? ONO₂ derivatives are rather high and vary with the number of ? ONO₂ groups. Considering the easiness for synthesis and relative stability, 2,4,6,8‐hexaazaadamantanenitrate is finally recommended as a potential candidate of a high‐energy density compound (HEDC). © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Theoretical Study on Gas—phase Pyrolytic Reactions of N—Ethyl,Isopropyl and N—t—Butyl Substituted 2—Aminopyrazine

Density functional theory (DFT) and ab initio methods were used to study gas‐phase pyrolytic reaction mechanisms of iV‐ethyl, N‐isopropyl and N‐t‐butyl substituted 2‐aminopyrazine at B3LYP/6–31G* and MP2/6–31G*, respectively. Single‐point energies of all optimized molecular geometries were calculated at B3LYP/6–311 + G(2d,p) level. Results show that the pyrolytic reactions were carried out through a unimolecular first‐order mechanism which were caused by the migration of atom H(17) via a six‐member ring transition state. The activation energies which were verified by vibrational analysis and correlated with zero‐point energies along the reaction channel at B3LYP/6–311 + G(2d,p) level were 252.02 kJ. mo^?1 (N‐ethyl substituted), 235.92 kJ‐mol^?1 (N‐t‐isopropyl substituted) and 234.27 kJ‐mol^?1 (N‐t‐butyl substituted), respectively. The results were in good agreement with available experimental data.     

Computational studies on carbohydrates: I. Density functional ab initio geometry optimization on maltose conformations

Ab initio geometry optimization was carried out on 10 selected conformations of maltose and two 2‐methoxytetrahydropyran conformations using the density functional denoted B3LYP combined with two basis sets. The 6‐31G* and 6‐311++G** basis sets make up the B3LYP/6‐31G* and B3LYP/6‐311++G** procedures. Internal coordinates were fully relaxed, and structures were gradient optimized at both levels of theory. Ten conformations were studied at the B3LYP/6‐31G* level, and five of these were continued with full gradient optimization at the B3LYP/6‐311++G** level of theory. The details of the ab initio optimized geometries are presented here, with particular attention given to the positions of the atoms around the anomeric center and the effect of the particular anomer and hydrogen bonding pattern on the maltose ring structures and relative conformational energies. The size and complexity of the hydrogen‐bonding network prevented a rigorous search of conformational space by ab initio calculations. However, using empirical force fields, low‐energy conformers of maltose were found that were subsequently gradient optimized at the two ab initio levels of theory. Three classes of conformations were studied, as defined by the clockwise or counterclockwise direction of the hydroxyl groups, or a flipped conformer in which the ψ‐dihedral is rotated by ∼180°. Different combinations of ω side‐chain rotations gave energy differences of more than 6 kcal/mol above the lowest energy structure found. The lowest energy structures bear remarkably close resemblance to the neutron and X‐ray diffraction crystal structures. © 2000 John Wiley & Sons, Inc. * J Comput Chem 21: 1204–1219, 2000     

Struktur und Eigenschaften des Hydridomethyltetrafluorphosphatanions, [CH3PF4H]—

Structure and Properties of the Hydridomethyltetrafluorophosphate Anion, [CH₃PF₄H]^— Methyltrifluorphosphorane reacts with strong fluoride donors like CsF, (CH₃)₄NF and (CH₃)₄PF with formation of the corresponding hydridomethyltetrafluorophosphates. The salts are characterized by NMR, IR and Raman spectroscopy. The experimental results show that only the trans isomer is formed. For both isomers theoretical calculations (B3LYP/6‐31+G* and RHF/6‐31+G*) were carried out. The difference for the Gibbs free energy between the isomers was calculated to be 35.4 kJ/mol (B3LYP/6‐31+G*). The RHF/6‐31+G* calculation yields, for the almost octahedral trans isomer, bond distances of r(PF) = 167 pm, r(PC) = 184.5 pm and r(PH) = 138.1 pm.     

DFT Studies on Detonation Properties,Pyrolysis Mechanism,and Stability of the Nitro and Hydroxyl Derivatives of Benzene

Ninety‐one nitro and hydroxyl derivatives of benzene were studied at the B3LYP/6‐31G?? level of density functional theory. Detonation properties were calculated using the Kamlet‐Jacobs equation. Three candidates (pentanitrophenol, pentanitrobenzene, and hexanitrobenzene) were recommended as potential high energy density compounds for their perfect detonation performances and reasonable stability. The pyrolysis mechanism was studied by analyzing the bond dissociation energy (BDE) and the activation energy (E_a) of hydrogen transfer (H–T) reaction for those with adjacent nitro and hydroxyl groups. The results show that E_a is much lower than BDEs of all bonds, so when there are adjacent nitro and hydroxyl groups in a molecule, the stability of the compound will decrease and the pyrolysis will be initiated by the H–T process. Otherwise, the pyrolysis will start from the breaking of the weakest C–NO₂ bond, and only under such condition, the Mulliken population or BDE of the C–NO₂ bond can be used to assess the relative stability of the compound.     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体热解机理与稳定性的理论研究

运用密度泛函理论和半经验分子轨道方法,对一系列高能杂环硝胺—反式-1,4,5,8-四硝基-1,4,5,8-四氮杂萘烷异构体的热解机理和稳定性进行了系统地计算研究。在B3LYP/6-31G**和PM3水平上,分别计算了标题物的化学键离解能(BDE)和热解反应活化能(Ea),并根据BDE和Ea数值考察了硝胺取代基对化合物稳定性和热解机理的影响;同时,还详细考察了BDE与Ea、化学键重叠布居数、前线轨道能级以及能隙之间的相关性。结果表明,由BDE、Ea和静态电子结构参数推断的标题物热稳定性和热解机理的结论基本是一致的,N-NO2键均裂是标题物的热解引发步骤,间位取代异构体较对位取代异构体稳定,而邻位取代的异构体稳定性最差。     

二氯亚甲基锗烯与乙烯环加成反应机理的理论研究

The mechanism of a cycloaddition reaction between singlet dichloromethylene germylene and ethylene has been investigated with B3LYP/6-31G＊ method, including geometry optimization and vibrational analysis for the involved stationary points on the potential energy surface. Energies for the involved conformations were calculated by CCSD（T）//B3LYP/6-31G＊ method. On the basis of the surface energy profile obtained with CCSD（T）// B3LYP/6-31G＊ method for the cycloaddition reaction between singlet dichloromethylene germylene and ethylene, it can be predicted that the dominant reaction pathway is that an intermediate INT1 is firstly formed between the two reactants through a barrier-free exothermic reaction of 61.7 kJ/mol, and the intermediate INT1 then isomerizes to an active four-membered ring product P2.1 via a transition state TS2, an intermediate INT2 and a transition state TS2.1, in which energy barriers are 57.7 and 42.2 kJ/mol, respectively.     

Theoretical investigations on structure,density, detonation properties,and sensitivity of the derivatives of PYX

The ? NH₂, ? NO₂, ? N₃, ? NHNO₂, and ? ONO₂ substitution derivatives of PYX (2,6‐bis(picrylamino)‐3,5‐dinitropyridine) were studied at the B3LYP/6‐31G** level of density functional theory. The sublimation enthalpies and heats of formation (HOFs) in gas phase and solid state of these compounds were calculated. The theoretical predicted density (ρ), detonation pressure (P), and detonation velocity (D) showed that these derivatives have better detonation performance than PYX. The effects of substituent groups on HOF, ρ, P, and D were discussed. The order of contribution of various groups to P and D was ? ONO₂ > ? NO₂ > ? NHNO₂ > ? N₃ > ? NH₂. Sensitivity was evaluated using the frontier orbital energies, bond orders, bond dissociation enthalpies (BDEs), and characteristic heights (h₅₀). The trigger bonds in the pyrolysis process for these PYX derivatives may be Ring‐NO₂, NH? NO₂, or O? NO₂ varying with the substituents. The h₅₀ of most compounds are larger than that of CL‐20, and those of ? NH₂, ? NO₂, and most ? ONO₂ derivatives are larger than that of RDX. The BDEs of the trigger bonds of all but the ? ONO₂ derivatives are sufficiently large. Taking both detonation performance and sensitivity into consideration, some derivatives of PYX may be good candidates of explosives. © 2012 Wiley Periodicals, Inc.     

Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.