新型均三氮苯类衍生物构效关系的模式识别研究

用量子化学密度泛函理论(DFT)、分子力学(MM)及模式识别方法,对34个新型均三氮苯类衍生物进行了结构活性关系研究.结果表明,影响此类衍生物分子除草活性的主要因素有C(1)上所带的电荷Qc(1),分子总能量E,分子生成热Ef,溶剂可及面积SGrid,分子体积V和C(3)上所带的电荷Qc(3)等参数.通过在C(1)上加上吸电子基,降低C(1)上的电荷,以及在C(3)上连接单取代氨基可以提高除草活性.所得模型对化合物生物活性有较好的预测效果.     

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.     

Investigation into the molecular structure and energetic stability of endohedral and exohedral metallofullerene derivatives of C24

Density-functional theory was applied to the investigation of the structural and electronic properties of C₂₄ fullerene derivatives. Transition metals (TMs) from groups 11 and 12, in various oxidation and spin-states, are inserted at either endohedral (TM@C₂₄) or exohedral (TM-C₂₄) sites and their subsequent energetic stabilities are assessed. With the exception of Ag@C₂₄, all derivatives are predicted to occupy a minimum on the potential energy surface. The optimized exohedral TM-C₂₄ geometries yield TM-C bond lengths that are consistent with comparable carbon-metal bond lengths, and the overwhelming majority of the derivatives result in a slight deformation of the C₂₄ cage as the bonding carbon takes on more sp³ character. All of the TM@C₂₄ equilibrium structures maintain the integrity of the cage structure with a moderate increase in the diameter. All neutral exohedral and endohedral complexes favor the low spin-state; conversely, all of the charged exohedral complexes prefer the high spin-state, with the exception of Cu-C₂₄¹⁺ molecular ion. The Group 12 charged endohedral derivatives prefer the low spin-state, whereas the Group 11 molecular ions do not necessarily exhibit a definitive trend. Analysis of the energetic data predicts that of the lowest energy endohedral molecular species only four are predicted to be energetically favorable in terms of insertion energy and an advantageous HOMO-LUMO gap: Cu@C₂₄²⁺, Ag@C₂₄¹⁺, Au@C₂₄³⁺, and Zn@C₂₄²⁺.     

Theoretical calculation of nitro-1-(2,4,6-trinitrophenyl)-1H-azoles energetic compounds

In order to study the properties of new energetic compounds formed by introducing nitroazoles into 2,4,6-trinitrobezene, the density, heat of formation and detonation properties of 36 nitro-1-(2,4,6-trinitrobenzene)-1H-azoles energetic compounds are studied by density functional theory, and their stability and melting point are predicted. The results show that most of target compounds have good detonation properties and stability. And it is found that nitro-1-(2,4,6-Trinitrophenyl)-1H-pyrrole compounds and nitro-1-(2,4,6-trinitrop-enyl)-1H-Imidazole compounds have good thermal stability, and their weakest bond is C NO₂ bond, the bond dissociation energy of the weakest bond is 222–238 kJ mol⁻¹ and close to 2,4,6-trinitrotoluene (235 kJ mol⁻¹). The weakest bond of the other compounds may be the C NO₂ bond or the N N bond, and the strength of the N N bond is related to the nitro group on azole ring.     

Molecular dynamics simulation and density functional theory insight into the cocrystal explosive of hexaazaisowurtzitane/nitroguanidine

Theoretical methods involving molecular dynamics (MD) simulation and density functional theory were performed to investigate the different molecular ratios, mechanical Properties, structure, trigger bond, and intermolecular interaction of hexaazaisowurtzitane (CL‐20)/nitroguanidine (NQ) cocrystal explosive. Results of MD simulation show that CL‐20 and NQ packed in ratios of 1:1 present the larger binding energy and better mechanical properties than any other molecular ratios, which indicates 1:1 cocrystal can form the stable crystal structure. Shorter length and larger dissociation energy of trigger bond in composite structure than in isolated CL‐20 component suggests that the cocrystal may exhibit less sensitive than CL‐20. Analyses of atoms in molecules, reduced density gradient, and natural bond orbital confirm that intermolecular interactions are mainly derived from a series of weak hydrogen bond and strong vdW forces, involving of NH···O, CH···O, CH···N, O···N, and O···O. Additionally, composite structures of 2 and 3 bringing us more attractive performance will act as a key role in constructing of CL‐20/NQ cocrystal explosive. © 2015 Wiley Periodicals, Inc.     

B炸药主要组分TNT和RDX分子间相互作用的理论研究

在B3LYP/6-31G(d)水平上研究了B炸药的主要成分--2,4,6-三硝基甲苯(TNT)与环三亚甲基三硝胺(RDX)分子间的相互作用, 得到了10种TNT＋RDX的全优化构型. 讨论了稳定构型在几何参数、稳定性、红外光谱和电荷分布上的差异. 借助自然键轨道(NBO)理论揭示了TNT与RDX分子间相互作用的本质, 主要由氢键所贡献. 分子间相互作用能在－3.930～－14.652 kJ•mol^－1之间, 经基组叠加误差(BSSE)校正, 相互作用能顺序为VI＞III＞V＞IV＞X＞I＞IX＞II＞VII＞VIII. 对全优化构型进行了热力学性质的分析, 探讨了由单体分子形成混合体系的热力学性质的变化, 结果发现, 形成分子间氢键是个放热过程. 运用Kamlet-Jacobs方程基于理论密度(ρ)估算了混合体系TNT＋RDX的爆轰性质爆速(D)和爆压(p), 与文献值进行比较表明理论计算方法和结果是可靠的.     

Theoretical investigation on density,detonation properties,and pyrolysis mechanism of nitro derivatives of benzene and aminobenzenes

Nitro derivatives of benzene and aminobenzenes are optimized at the DFT‐B3LYP/6‐31G* level. The heat of formation (ΔH_f) and crystal theoretical density (ρ) are estimated to evaluate the detonation properties using the modified Kamlet–Jacobs equations. Thermal stability and the pyrolysis mechanism of the title compounds are investigated by calculating the bond dissociation energies (BDE) at the unrestricted B3LYP/6‐31G* level. The kinetic parameter and the static electronic structural parameters can be used to predict the stability and the relative magnitude of the impact sensitivity of homologues. According to the quantitative standard of the energy and the stability as an HEDC, the title compounds having more than four nitro groups satisfy this requirement. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Nitrogen-rich salts based on energetic nitroaminodiazido[1,3,5]triazine and guanazine

Highly dense nitrogen-rich ionic compounds are potential high-performance energetic materials for use in military and industrial venues. Guanazinium salts with promising energetic anions and a family of energetic salts based on nitrogen-rich cations and the 6-nitroamino-2,4-diazido[1,3,5]triazine anion (NADAT) were prepared and fully characterized by elemental analysis, IR spectroscopy, (1)H NMR and (13)C NMR spectroscopy, and differential scanning calorimetry (DSC). The crystal structures of neutral NADAT (2) and its biguanidinium salt 5 were determined by single-crystal X-ray diffraction (2: orthorhombic, Pnma; 5: monoclinic, P2(1)). Additionally, the isomerization behavior of 2 in solution was investigated by proton-decoupled (13)C and (15)N NMR spectroscopy. All the new salts exhibit desirable physical properties, such as relatively high densities (1.63-1.78 g cm(-3)) and moderate thermal stabilities (T(d) = 130-196 °C for 3-10 and 209-257 °C for 11-15). Theoretical performance calculations (Gaussian 03 and Cheetah 5.0) gave detonation pressures and velocities for the ionic compounds 3-15 in the range of 21.0-30.3 GPa and 7675-9048 m s(-1), respectively, which makes them competitive energetic materials.     

Molecular design of all nitrogen pentazole‐based high energy density compounds with oxygen balance equal to zero

We designed a new family of pentazole‐based high energy density compounds with oxygen balance equal to zero by introducing −NH₂, −NO₂, −N₃, −CF₂NF₂, and −C[NO₂]₃, and the properties including density, heats of formation, detonation performances, and impact sensitivity were investigated using density functional theory. The results show that half of these new energetic molecules exhibit higher densities than RDX (1.82 g/cm³), in which H5 gives the highest density of 2.09 g/cm³. Among all the 54 designed molecules, 22 compounds have higher D and P than RDX and eleven compounds have higher D and P than HMX, indicating that designing the pentazole‐based derivatives with oxygen balance equal to zero is a very effective way to obtain potential energetic compounds with outstanding detonation properties. Taking both the detonation performance and stability into consideration, nine compounds may be recognized as potential candidates of high energy density compounds. It is expected that our results will contribute to the theoretical design of new‐generation energetic explosives.     

Looking for high energy density compounds applicable for propellant among the derivatives of DPO with -N3, -ONO2, and -NNO2 groups

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.     

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 predictions on pentaerythritol tetranitrate-based high energy density compounds

A novel family of pentaerythritol tetranitrate (PETN) derivatives based parent PETN skeleton were designed by introducing two energetic groups –NF₂ and –NO₂. Their electronic structure, heats of formation, detonation properties, impact sensitivity, and thermal stability were investigated by using density functional theory. The findings reveal that most of the title compounds have good detonation performance. The –NF₂ group played an important role in improving the densities, heats of detonation, and detonation properties of the designed molecules. The values of h₅₀ for almost all the PETN derivatives are higher than that of octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine. An analysis of bond dissociation energy suggests that the N-NO₂ bond tends to be a trigger bond in thermal decomposition. Taking both detonation properties and thermal stabilities into consideration, the three compounds may be selected as potential high-energy-density compounds.     

环三甲撑三硝胺(RDX)结构和性质的DFT研究

环三甲撑三硝胺(RDX)是综合性能好、使用极广的高能炸药。本文用密度泛函理论(DFT)B3LYP方法,取6-31G*基组,求得其分子几何构型、电子结构、IR谱和热力学性质。全优化几何构型和电子结构均具有C3V对称性。在相邻原子之间以NNO2键的Mulliken集居数最小,表明其间电子分布较少,预示其为热解和起爆的引发键。IR谱与实验结果良好相符。计算所得298-1200K的热力学性质,对进一步深入研究RDX的反应和性质有助。     

Computational Evaluation of Me2TCCP as Lewis Acid

Supramolecular adducts between dimethyl-2,2,3,3-tetracyanocyclopropane (Me₂TCCP) with 21 small (polar) molecules and 10 anions were computed with DFT (B3LYP-D3/def2-TZVP). Their optimized geometries were used to obtain interaction energies, and perform energy decomposition and ‘atoms-in-molecules’ analyses. A set of 38 other adducts were also evaluated for comparison purposes. Selected examples were further scrutinized by inspection of the molecular electrostatic potential maps, Noncovalent Interaction index plots, the Laplacian, the orbital interactions, and by estimating the Gibbs free energy of complexation in hexane solution. These calculations divulge the thermodynamic feasibility of Me₂TCCP adducts and show that complexation is typically driven by dispersion with less polarized partners, but by orbital interactions when more polarized or anionic guests are deployed. Most Me₂TCCP adducts are more stable than simple hydrogen bonding with water, but less stable than traditional Lewis adducts involving Me₃B, or a strong halogen bond such as with Br₂. Several bonding analyses showed that the locus of interaction is found near the electron poor sp³-hydridized (NC)₂C−C(CN)₂ carbon atoms. An empty hybrid σ*/π* orbital on Me₂TCCP was identified that can be held responsible for the stability of the most stable adducts due to donor-acceptor interactions.     

Design and properties prediction of modified CL-20 energetic derivatives

We designed a series of energetic compounds based on the CL-20 molecular skeleton, and the properties including molecular geometric structures, electronic structures, density, heat of formation, detonation performances, and impact sensitivity were evaluated using density functional theory (DFT). The results indicate that five molecules have higher density values than that of Octahydro-1,3,5,7-tetranitro-1,3,5,7-tetrazocine (HMX; 1.91 g/cm³) and A4 has a larger density value (2.07 g/cm³) than that of CL-20 (2.04 g/cm³). In addition, most of the molecules have better detonation performances and stability than those of CL-20, with A4 showing much greater detonation velocity (9.93 km/s) and pressure (47.32 GPa) than those of CL-20 with a h₅₀ value of 14.02 cm. Taking both excellent detonation performance and low sensitivity into consideration, all seven compounds except for A3 and A5 are considered as potential energetic compounds. These theoretically calculated results would be conducive to the design and synthesis of novel nitramine energetic compounds.     

Theoretical studies on high energetic density nitramine explosives containing pyridine

The insensitive property of explosives containing pyridine is combined with the high energy of nitramine explosives,and the concept of new nitramine explosives containing pyridine is proposed,into which nitramine group with N N bonds is introduced as much as possible.Based on molecular structures of nitramine compounds containing pyridine,density functional theory(DFT) calculation method was applied to study designed molecules at B3LYP/6-31+G(d) level.The geometric and electronic structures,density,heats of formation(HOF),detonation performance and bond dissociation energies(BDE) were investigated and comparable to 1,3,5-trinitro-1,3,5-triazinane(RDX) and 1,3,5,7-tetranitro-1,3,5,7-tetrazocane(HMX).The simulation results reveal that molecules B and D perform similarly to traditionally used RDX.Molecule E outperform RDX,with performance that approach that of HMX and may be considered as potential candidate of high energy density compound(HEDC).These results provide basic information for molecular design of novel high energetic density compounds.     

Computational study of the reactivity of cytosine derivatives

The aim of the present study is to provide computational insight using dispersion‐corrected density‐functional calculations into the reactivity properties of modified cytosine in the gas phase and in aqueous solution, whereby special emphasis is put on systems that are obtained through demethylation and methylation. Since this field is relatively incipient, our goal is to identify relationships between reactivity and stability for the modified compounds to understand their biological functionalities. Our results show that addition of a methyl, hydroxylmethyl, formyl, or carboxyl group reduces the length of the nearest hydrogen bond between the cytosine–guanine (CG) base pair and increases the length of the longest hydrogen bond of the DNA base pair. © 2017 Wiley Periodicals, Inc.     

s-四嗪-水簇复合物的理论研究

用量子化学B3LYP方法和6-31＋＋G**基函数研究了s-四嗪-水簇复合物基态分子间相互作用, 并进行了构型优化和频率计算, 分别得到无虚频稳定的s-四嗪-(水)₂复合物、s-四嗪-(水)₃复合物和s-四嗪-(水)₄复合物6个、9个和12个. 复合物存在较强的氢键作用, 复合物结构中形成一个N…H—O氢键并终止于O…H—C氢键的氢键水链构型最稳定. 经基组重叠误差和零点振动能校正后, 最稳定的1∶2, 1∶3和1∶4(摩尔比)复合物的结合能分别是41.35, 70.9和 94.61 kJ/mol. 振动分析显示氢键的形成使复合物中水分子H—O键对称伸缩振动频率减小(红移). 研究表明N…H键越短, N…H—O键角越接近直线, 稳定化能越大, 氢键作用越强. 同时, 用含时密度泛函理论方法在TD-B3LYP/6-31＋＋G**水平计算了s-四嗪单体及其氢键复合物的第一¹(n, p*)激发态的垂直激发能.     

Theoretical studies on heats of formation,detonation properties,and bond dissociation energies of monofurazan derivatives

The heats of formation (HOFs) for a series of monofurazan derivatives were calculated by using density functional theory. It is found that the ? CN or ? N₃ group plays a very important role in increasing the HOF values of the furazan derivatives. The detonation velocities and detonation pressures of the furazan derivatives are evaluated at two different levels. The results show that the ? NF₂ group is very helpful for enhancing the detonation performance for the furazan derivatives, but the case is quite the contrary for the ? CH₃ group. An analysis of the bond dissociation energies and bond orders for the weakest bonds indicate that the substitutions of ? CN group are favorable and enhances the thermal stability of the furazan derivatives, but the ? NO₂ groups produce opposite effects. These results provide basic information for the molecular design of novel high‐energy density materials. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2010     

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.