用ab initio方法研究NaN_3的几何构型、电子结构和热力学、动力学性质

用ａｂｉｎｉｔｉｏＭＯ方法，在ＭＰ２（ｆｕｌ）／６３１１Ｇ水平下，全优化计算了叠氮化钠（ＮａＮ３）分子的线状和环状两种稳定构型及其转化过渡态的几何参数、电荷分布、分子总能量和振动频率，并研究了它们的热力学性质及转化速率常数和平衡常数．结果表明，线状比环状构型稍稳定（能量低６．０４ｋＪ／ｍｏｌ）；两者相互转化的能垒分别为１３．１５ｋＪ／ｍｏｌ（线型→环状）和７．１１ｋＪ／ｍｏｌ（环状→线型）．热力学和动力学计算均表明：ＮａＮ３通常主要以线型结构存在（占８５％以上），且随温度升高而增多（在１０００Ｋ大于９１％）．     

Ab initio Study on the Molecular Geometry,Electronic Structure and Thermodynamic Properties of Benzotrifuroxan(BTF)

１ＩＮＴＲＯＤＵＣＴＩＯＮＢｅｎｚｏｔｒｉｆｕｒｏｘａｎ（ＢＴＦ）ｉｓａｐｏｗｅｒｆｕｌ,ｈｙｄｒｏｇｅｎ ｆｒｅａｎｄｒｅｌａｔｉｖｅｌｙｓｅｎｓｉｔｉｖｅｅｘ ｐｌｏｓｉｖｅ．Ｓｉｎｃｅｔｈｉｓｃｏｍｐｏｕｎｄｗａｓｆｉｒｓｔｐｒｅｐａｒｅｄ,ｉ...     

Ab initio Study on the Molecular Geometry,Electronic Structure and Thermodynamic Properties of Benzotrifuroxan(BTF)

１ＩＮＴＲＯＤＵＣＴＩＯＮＢｅｎｚｏｔｒｉｆｕｒｏｘａｎ（ＢＴＦ）ｉｓａｐｏｗｅｒｆｕｌ，ｈｙｄｒｏｇｅｎｆｒｅａｎｄｒｅｌａｔｉｖｅｌｙｓｅｎｓｉｔｉｖｅｅｘｐｌｏｓｉｖｅ．Ｓｉｎｃｅｔｈｉｓｃｏｍｐｏｕｎｄｗａｓｆｉｒｓｔｐｒｅｐａｒｅｄ，ｉ...     

四唑及其衍生物的理论研究Ⅷ: 硝氨四唑衍生物的从头算

运用从头计算法,在HF/6-31G^*水平下,全优化计算了7种硝氨四唑衍生物的分子几何和电子结构。结果表明,标题物的四唑环近似为平面构型;2H式中性分子的芳香性大于相应的1H式异构体。5-硝氨-1H四唑分子内氢键使硝氨基与环共面;其余标题物中硝氨基与环垂直。不同水平下的总能量计算表明,标题物中电子相关效应显著;1H式中性分子较2H式的能量高,5-硝氨四唑负离子在三个负离子中最稳定。根据电荷分布阐明了四唑环质子化位置和标题物与金属的配位方式。此外还计算了各标题物的红外光谱及热力学性质。     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Difluoroamine Complex

Introduction　　ApplicationsofabinitiocalculationstointermolecularinteractionsincludingeitherweakvanderWaalsorstrongerhydrogenbondinghavedrawnmuchattentioninthepastdecadesbecausetheyareimportantinawiderangeofphysical,chemicalandbiologicalfields .1 5Inre centyears ,wehaveappliedtheintermolecularinteractionstoenergeticsystemsandobtainedsomemeaningfulinfor mationthatisvaluableforthestudyofenergeticmateri als .6 14 Thebehaviorofmolecularcomplexesisusuallybe tweentwoextremes :thegasphaseandthecrys…     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Dimethylnitroamine Clusters

Ab initio SCF and Mφller-Plesset correlation correction methods in combination with counterpose procedure for BSSE correction have been applied to the theroetical studying of dimethylnitroamine and its dimers and trimers.Three optimized stable dimers and two trimers have been obtained.The corrected binding energies of the most stable dimer and trimer were predicted to be -24.68kJ/mol and -47.27kJ/mol,respectively at the MP2/6-31G^*//HF/6-31G^* level.The proportion of correlated interation energies to their total interaction energies for all clusters was at least 29.3 percent,and the BSSE of ΔE(MP2) was at least 10.0kJ/mol.Dispersion and/or electrostatic force were dominant in all clusters.There exist cooperative effects in both the chain and the cyclic trimers.The vibrational frequencies associated with N-O stretches or wags exhibit slight red shifts,but the modes associated with the motion of hydrogen atoms of the methyl group show somewhat blue shifts with respect to those of monomer.Thermodynamic properties of dimethylnitroamine and its clusters at different temperatures have been calculated on the basis of vibrational analyses.The changes of the Gibbs free energies for the aggregation from monomer to the most stable dimer and trimer were predicted to be 14.37kJ/mol and 30.40kJ/mol,respectively,at 1 atm and 298.15K.     

Ab initio Study on the Intermolecular Interaction and Thermodynamic Properties of Methyl Nitrate Dimer

Three stable dimers of methyl nitrate have been obtained and their geometries have been fully optimized at the HF/6‐31G,. level. Binding energies have been calculated with correction for the basis set superposition error (BSSE) and zero point energy (ZPE). The cyclic overlap‐type structure, the binding energy of which is 11.97 kJ/mol at the MP4SDTQ/6‐31G. / HF/6‐31G. level, is the most stable. No intermolecular hydrogen bond was found, and the charge transfer between two subsystems is minute. The thermodynamic properties of methyl nitrate and its dimers have been calculated based on the vibrational analysis and statistical thermodynamics.     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

钒(Ⅴ)巯基苯酚双齿配合物从头算研究

运用Ｇ９４Ｗ量子化学程序包,采用ＬＡＮＬ２ＤＺ基组,对钒与Ｈ２ｍｐ（ｏ－巯基苯酚）形成的双齿配合物结构单元Ｖ（ｍｐ）３ｎ－（ｍｐ＝ＯＣ６Ｈ４Ｓ２－,ｎ＝１,２）进行从头算研究．探讨这些配合物结构单元的稳定性、分子轨道能量、原子净电荷布居规律以及一些前沿分子轨道的组成特征等,为过渡金属双齿配合物的合成、分子组装分析及其在酶和蛋白质中的活性部位研究提供理论参考．     

2,4-二吡啶基方酸衍生物非线性光学性质的从头算研究

采用HF/4-31G方法优化分子几何构型,在此基础上用CPHF研究了系列2,4-二吡啶基方酸衍生物的二阶非线性光学系数βvec,并对影响βvec的因素进行了探讨.结果表明:此类化合物均有很大的βvec,它们的非线性光学性质与其分子结构有着密切的关系.     

戊唑醇的分子结构、理论计算及热力学性质研究

在甲苯溶剂中利用缓慢蒸发法得到1-(4-氯苯基)-4,4-二甲基-3-(1H-1,2,4-三唑基甲基)戊醇-3（戊唑醇）的单晶,通过 X射线单晶结构分析法测定其晶体结构,晶体属单斜晶系,空间群为P2(1)/c,晶胞参数为：a = 1.1645(1) nm,b = 1.6768(2) nm,c = 1.7478(2) nm,β= 92.055(2),Dcalc.= 1.199 g/cm3,Z = 4,F(000)= 264。运用密度泛函理论 (DFT) B3LYP得到其优化几何构型并得到其频率。计算得到的结构参数与相应的实验值十分接近。运用微热量仪对标题物进行比热容测定,在所测温度范围283~353 K内,比热容随温度呈稳定的线性变化,根据测定的比热容方程,计算出戊唑醇以298.15 K为基础在283~353 K温区的的热力学函数：焓、熵和吉布斯自由能。     

水合氢离子团簇从头算和ABEEM/MM的理论研究

应用高水平的从头计算方法和ABEEM/MM模型, 研究了水合氢离子团簇H3O＋(H2O)n (n＝1～6), 优化得到了低能构象, 探讨了其结合能和稳定性, 显示出H3O＋(H2O)3局域结构的优势存在．对H3O＋(H2O)6VIa团簇的ABEEM电荷分布进行分析, 表明第一水合层水分子与水合氢离子之间的氢键相互作用要明显强于与第二水层水分子的氢键相互作用. 研究结果表明, ABEEM/MM方法计算的结果和从头算得到的结果存在很好的一致性.     

Comparison of ab initio and density functional methods for vibrational analysis of TeCl4

Vibrational analysis of tellurium tetrachloride, TeCl₄, was performed with Hartree–Fock (HF), MP2, and generalized gradient approximation density functional theory (DFT) methods supplemented with polarized double-zeta split valence (DZVP) basis sets and relativistic effective core potentials (RECP) of Hay and Wadt. The molecular geometry is best reproduced at the HF and MP2/RECP+DZVP [polarized Hay and Wadt RECP for Te and 6–31G(d) basis set for Cl] levels of theory. The DFT methods gave rise to poorer results, especially those using Becke's 1988 exchange functional. Generally, the vibrational frequencies calculated by the MP2 and B3-type DFT methods with the all electron and RECP+DZVP basis sets as well as at the HF/RECP level were in satisfactory accord with the experimental data. The agreement was good enough to assist the assignment of the measured vibrational spectra. The best agreement with the experimental vibrational frequencies was achieved with the scaled HF/RECP force field. Consistent results were obtained for the unobserved A₂ (ν₄) fundamental, where the results of the best methods were within 4 cm⁻¹. The best force fields were obtained with the following methods: Becke3–Lee–Yang–Parr and Becke3–Perdew/all electron basis, MP2 and Becke3-Perdew/RECP+DZVP, and HF/RECP. The methods using RECPs are advantageous for large-scale computations. The RECP basis set effectively compensates the errors of the HF method for TeCl₄; however, it provides poor results with correlated methods. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 308–318, 1998     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Nitroamine Dimers

Ab initio self-consistent field(SCF) and Mφller-Plesset correlation correction methods employing 6-31G^** basis set have been applied to the optimizations of nitroamine dimers.The binding energies have been corrected for the basis set superposition error (BSSE) and the zero-point energy.Theree optimized dimers have been obtained.The BSSE corrected binding energy of the most stable dimer is predicted to be -31.85kJ/mol at the MP4/6-31G^**//MP2/6-31G^** level.The energy barriers of the Walden conversion for -NH2 group are 19.7kJ/mol and 18.3kJ/mol for monomer and the most stable dimer,respectively.The molecular interaction makes the internal rotation around N1-N2 even more difficult.The thermodynamic properties of nitroamine and its dimers at different temperatures have been calculated on the basis of vibrational analyses.The change of the Gibbs free energy for the aggregation from monomer to the most stable dimer at standard pressure and 298.2 K is predicted to be 14.05kJ/mol.     

Theoretical study on tetrazole and its derivatives (3) MP2 and thermodynamic calculations of amino derivatives of tetrazole

Fully optimized geometries and electronic structures of amino derivatives of tetrazole are obtained at MP2/6-31G* level. The tetrazole rings are planar and aromatic for all the amino derivatives of tetrazole. The amino group is not co-planar with the ring and its conformation is mainly determined by the lone pair electronic repulsion between the substituent and the ring. N(4) atom is more negatively charged and is the most probable coordination site. The energy gaps between LUMOs and HOMOs of 2H-aminotetrazoles and C-aminotetrazole neutrals are smaller than those of the corresponding 1H-isomers and N-aminotetrazole neutrals respectively. The IR frequencies, thermodynamic properties and temperature-dependent functions for heat capacities in the form (a bT cT2) in the 300-1000K range are reported.     

A comparative study of DFT and traditional ab initio methodologies on the OsO4 molecule

The performance of different conventional ab initio methodologies and density functional procedures is compared through its application to the theoretical calculation of the bond distance and harmonic vibrational frequencies of the OsO₄ molecule. The problem of the basis set is first considered, with up to nine different basis sets being tested in calculations using the hybrid Becke3LYP density functional, and the most appropriate basis set is used in the comparison of Hartree–Fock, post‐Hartree–Fock, and density functional methods. The post‐Hartree–Fock methods analyzed are MP2, CISD, and CCSD(T), and the density functionals tested are SVWN, BLYP, BPW91, and Becke3LYP. The results show that for this particular system density functional methods perform better than do HF‐based methods with the exception of CCSD(T), which gives the best overall results. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 77: 544–551, 2000     

Thermodynamic and ab initio analysis of the controversial enthalpy of formation of formaldehyde.

There are two values, -26.0 and -27.7 kcal mol(-1), that are routinely reported in literature evaluations for the standard enthalpy of formation, Delta(f) H(o)(298), of formaldehyde (CH(2)=O), where error limits are less than the difference in values. In this study, we summarize the reported literature for formaldehyde enthalpy values based on evaluated measurements and on computational studies. Using experimental reaction enthalpies for a series of reactions involving formaldehyde, in conjunction with known enthalpies of formation, its enthalpy is determined to be -26.05+/-0.42 kcal mol(-1), which we believe is the most accurate enthalpy currently available. For the same reaction series, the reaction enthalpies are evaluated using six computational methods: CBS-Q, CBS-Q//B3, CBS-APNO, G2, G3, and G3B3 yield Delta(f) H(o)(298)=-25.90+/-1.17 kcal mol(-1), which is in good agreement to our experimentally derived result. Furthermore, the computational chemistry methods G3, G3MP2B3, CCSD/6-311+G(2df,p)//B3LYP/6-31G(d), CCSD(T)/6-311+G(2df,p)//B3LYP/6-31G(d), and CBS-APNO in conjunction with isodesmic and homodesmic reactions are used to determine Delta(f) H(o)(298). Results from a series of five work reactions at the higher levels of calculation are -26.30+/-0.39 kcal mol(-1) with G3, -26.45+/-0.38 kcal mol(-1) with G3MP2B3, -26.09+/-0.37 kcal mol(-1) with CBS-APNO, -26.19+/-0.48 kcal mol(-1) with CCSD, and -26.16+/-0.58 kcal mol(-1) with CCSD(T). Results from heat of atomization calculations using seven accurate ab initio methods yields an enthalpy value of -26.82+/-0.99 kcal mol(-1). The results using isodesmic reactions are found to give enthalpies more accurate than both other computational approaches and are of similar accuracy to atomization enthalpy calculations derived from computationally intensive W1 and CBS-APNO methods. Overall, our most accurate calculations provide an enthalpy of formation in the range of -26.2 to -26.7 kcal mol(-1), which is within computational error of the suggested experimental value. The relative merits of each of the three computational methods are discussed and depend upon the accuracy of experimental enthalpies of formation required in the calculations and the importance of systematic computational errors in the work reaction. Our results also calculate Delta(f) H(o)(298) for the formyl anion (HCO(-)) as 1.28+/-0.43 kcal mol(-1).     

Spectra and structure of organophosphorus compounds: XLIX. Microwave,infrared, and Raman spectra,electric dipole moment,molecular structure,and ab initio calculations of dimethylphosphonothioic fluoride

The microwave spectra of (CH₃)₂PSF, (CH₃)(CD₃)PSF, (CD₃)₂PSF, and (CH₃)₂P³⁴SF have been investigated from 20.0 to 40.0 GHz. Botha-type R branch andc-type Q branch transitions have been measured in the ground states of each isotopic species. From a least-square adjustment to fit 12 rotational constants, the following structural parameters were obtained:r(P–F)=1.582 ± 0.003 Å;r(P=S)=1.902 ± 0.001 Å;r(P-C)=1.800 ± 0.001 Å;r(C-H)=1.088 ± 0.002 Å; HCP=109.28 ± 0.12°; SPF=114.50 ± 0.13°; and SPC=116.33 ± 0.06°. From Stark effect measurements, the dipole moment components have been determined to be ¦ _a ¦ =3.556 ± 0.005; ¦ _c ¦=2.026 ± 0.009; and ¦ _t ¦=4.093 ± 0.009 (D). The Raman spectra (3200 to 100 cm^–1) of each isotopic species have been measured for the solid, and liquid and qualitative depolarization values obtained. Additionally, the mid-infrared spectra (3200 to 500 cm^–1) of the solids have been recorded. Proposed assignments of the normal modes have been made on the basis of Raman depolarization values and group frequencies which are supported by normal coordinate analysis utilizing an ab initio force field. Optimized structural parameters have been obtained with both the 3-21G^* and 6-31G^* basis sets. These results are compared to the corresponding quantities for several similar molecules.For part XLVIII, seeJ. Raman Spectrosc.1922,23, 107.     

Ionic dissociations of chlorosulfonic acid in microsolvated clusters: A density functional theory and ab initio MO study

Ionic dissociation of chlorosulfonic acid (HSO3Cl) in the molecular clusters HSO3Cl-(H2O)n (n = 1-4) and HSO3Cl-NH3-(H2O)n (n = 0-3) was investigated by density functional theory and ab initio molecular orbital theory. The equilibrium structures, binding energies, and thermodynamic properties, such as relative enthalpy and relative Gibbs free energy, and were calculated using the hybrid density func- tional (B3LYP) method and the second order M?ller-Plesset approximation (MP2) method with the 6-311 G** basis set. Chlorosulfonic acid was found to require a minimum of three water molecules for ionization to occur and at least one water molecule to protonate ammonia. The corresponding clusters with fewer water molecules were found to be strongly hydrogen-bonded. The related properties and acid strength of chlorosulfonic acid were discussed and compared to the acid strengths of perchloric acid and sulfuric acid in the context of clusters with ammonia and water. The relative stabilities of these clusters were also investigated.     

The mechanism of formamide hydrolysis in water from ab initio calculations and simulations

The neutral hydrolysis of formamide in water is a suitable reference to quantify the efficiency of proteolytic enzymes. However, experimental data for this reaction has only very recently been obtained and the kinetic constant determined experimentally is significantly higher than that predicted by previous theoretical estimations. In this work, we have investigated in detail the possible mechanisms of this reaction. Several solvent models have been considered that represent a considerable improvement on those used in previous studies. Density functional and ab initio calculations have been carried out on a system which explicitly includes the first solvation shell of the formamide molecule. Its interaction with the bulk has been treated with the aid of a dielectric continuum model. Molecular dynamics simulations at the combined density functional/molecular mechanics level have been carried out in parallel to better understand the structure of the reaction intermediates in aqueous solution. Overall, the most favored mechanism predicted by our study involves two reaction steps. In the first step, the carbonyl group of the formamide molecule is hydrated to form a diol intermediate. The corresponding transition structure involves two water molecules. From this intermediate, a water-assisted proton transfer occurs from one of the hydroxy groups to the amino group. This reaction step may lead either to the formation of a new reaction intermediate with a marked zwitterionic character or to dissociation of the system into ammonia and formic acid. The zwitterionic intermediate dissociates quite easily but its lifetime is not negligible and it could play a role in the hydrolysis of substituted amides or peptides. The predicted pseudo-first-order kinetic constant for the rate-limiting step (the first step) of the hydrolysis reaction at 25 degrees C (3.9x10(-10) s(-1)) is in excellent agreement with experimental data (1.1x10(-10) s(-1)).