硝酰胺二聚体静电能和交换能的理论计算

根据对称性匹配微扰理论, 并运用多个微扰和非微扰方法, 计算了硝酰胺二聚体在不同分子间距(R)时静电能和交换能. 这些具有轨道弛豫的静电能不仅含有第4级单、双、四和三重激发态项, 而且含有CCSD的第5级和更高级的能量项. 同时发现: 第4级的三重激发态能量项比第5级和更高级的能量项之和还重要. 求得的含有分子内电子相关效应的交换能达到了CCSD水平. 用于计算交换能关联校正项的单对交换近似在硝酰胺二聚体的范德华最小距离0.42 nm附近区域才较合理. 在R为0.32～1.42 nm范围, 静电能与R的关系有两种: 一是在小于等于0.47 nm时, (R^－7.64; 二是大于0.47 nm时, (R^－3.97. 交换能具有明显的短程作用特点, 其与R间的关系为指数衰减: 21.061exp(－R/0.318). 最后发现: 在硝酰胺二聚体中, 分子内的电子相关效应对和的影响很显著.     

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.     

Ab initio Studies on Intermolecular Interaction of Formamide and Hydroxyacetonitrile Dimers

The structures, the binding energies and the thermodynamic properties of formamide and hydroxyacetonitrile(HAN) dimers have been studied by means of the self-consistent ab initio Hartree-Fock and the second-order Mφller-Plesset correlation energy correction methods. The counterpoise procedure was used to check the basis set superposition error(BSSE) of the binding energies. There exist cyclic structures in a formamide dimer(Ⅰ), a HAN dimer(Ⅱ) and their heterodimer(Ⅲ). The corrected binding energies for dimers Ⅰ, Ⅱ and Ⅲ are respectively -45.53, -45.83 and -43.89 kJ/mol at the MP2/aug-cc-p VDZ//HF/aug-cc-p VDZ level. The change of the Gibbs free energies(ΔG) in the process of Ⅰ Ⅱ→2Ⅲ was predicted to be -2.74 kJ/mol at 298.15 K. Dimer Ⅲ can be spontaneously produced in the mixture of formamide and HAN, which is in agreement with the experimental fact that most cyanohydrins are capable of interacting with dipeptide cyclo-His-Phe(CHP).     

CH3F分子间相互作用的ab initio研究

在MP2/6-311+ +G(3d,3p)电子相关校正水平上,对CH₃F二聚体可能存在的几何构型进行全自由度能量梯度优化和频率验证,发现3种势能面上有极小点的构型.进一步在高级电子相关校正的MP4S-DTQ、CCSD(T)/6-311+ +G(3df,3pd)方法水平上,对其中总能量最小的构型进行精确计算,得到二聚体的结合能为-9.707kJ/mol.研究结果支持了由光谱实验结果推测的构型,解释了CH₃F二聚体的谐振频率的多样性.     

硝酸甲酯分子间相互作用的DFT和ab initio比较

用密度泛函理论(DFT)和从头算(ab initio)方法,分别在B3LYP/6 31G和HF/6 31G水平上求得硝酸甲酯三种二聚体的全优化几何构型和电子结构,并用6 311G和6 311＋＋G基组进行总能量计算.对HF/6 31G计算结果进行MP4SDTQ电子相关校正.在各基组下均进行基组叠加误差(BSSE)和零点能(ZPE)校正求得结合能.对6 31G优化构型作振动分析并基于统计热力学求得200～600 K温度下单体和二聚体的热力学性质.详细比较两种方法的相应计算结果,发现DFT求得的分子间距离较短,分子内键长较长,所得结合能均小于相应ab initio计算值.     

Intermolecular interaction energy of CH4 trimer by symmetry-adapted perturbation theory

The interaction energy for the cyclic CH₄ trimer is studied in terms of symmetry-adapted perturbation theory. The interaction energy around the van der Waals minimum is dominated by attractive dispersion energy, and the repulsive contribution at the smaller angle region is due to the first-order exchange energy. The total interaction energy is approximated by additive two-body components, because of a mutual cancellation between nonadditive three-body ones.     

Ab initio studies of methanol, methanethiol and methylamine dimer

Quantum chemical calculations are used to provide structural, vibrational and energetical information on the dimers of the methanol, methylamine and methanethiol systems. These systems were studied employing the DFT(B3LYP) and MP2 methods together with the 6-31+G^** and 6-311+G^** basis sets. We found two distinct potential minima for methylamine (one of them is a transition structure) and methanethiol, and one for the methanol dimer. The properties of these dimers are compared with those of the dimers (H₂O)₂, (NH₃)₂ and (CH₃SH)₂. The interactions in these dimers were analyzed using electron density properties at the bond critical point.     

Xe-N2O复合物的分子间势能面和振转光谱的理论研究

采用CCSD(T)方法研究了范德华分子体系Xe—N2O复合物的势能面和振转光谱性质,研究表明,该势能面有两个极小点,分别对应T构型和线性Xe—ONN构型,采用离散变量表象和Lanczos算法计算了体系的振转能级,计算结果表明,CCSD(T)势能面支持97个振动束缚态,并对能级进行了指认,计算得到的Xe—N2O转动跃迁频率与实验值吻合得很好。     

He-N2O的从头算势能面及振转能级

采用超分子MP4方法和较大的基组计算得到了He-N₂O体系的分子间势能面,发现该势能面有3个极小值点,分别对应T形构型及两个线性He-ONN和He-NNO构型.同时采用离散变量表象方法预测了体系的振转能级,计算结果表明,MP4势能面支持5个振动束缚态.     

Finite-Temperature Dimer Method for Finding Saddle Points on Free Energy Surfaces

The dimer method and its variants have been shown to be efficient in finding saddle points on potential surfaces. In the dimer method, the most unstable direction is approximately obtained by minimizing the total potential energy of the dimer. Then, the force in this direction is reversed to move the dimer toward saddle points. When the finite-temperature effect is important for a high-dimensional system, one usually needs to describe the dynamics in a low-dimensional space of reaction coordinates. In this case, transition states are collected as saddle points on the free energy surface. The traditional dimer method cannot be directly employed to find saddle points on a free energy surface since the surface is not known a priori. Here, we develop a finite-temperature dimer method for searching saddle points on the free energy surface. In this method, a constrained rotation dynamics of the dimer system is used to sample dimer directions and an efficient average method is used to obtain a good approximation of the most unstable direction. This approximated direction is then used in reversing the force component and evolving the dimer toward saddle points. Our numerical results suggest that the new method is efficient in finding saddle points on free energy surfaces. © 2019 Wiley Periodicals, Inc.     

InSb的Li替位形成能的从头计算

采用基于平面波展开的第一原理赝势法计算了锑化铟在锂替位到铟位置时的各种情况下的形成能与电子结构,讨论了锂替位的体积变化,电荷分布、能带结构及电子态密度等性质.结果表明,对于闪锌矿结构的InSb,锂的各种替位形成能大致在每个锂原子-2.2eV左右.该结果表明,不可能在嵌入初期Li插入到间隙位置之前发生替位反应,与实验结果一致.     

Ab initio conformational analysis of cyclooctane molecule

The potential energy surface (PES) for the cyclooctane molecule was comprehensively investigated at the Hartree–Fock (HF) level of theory employing the 3–21G, 6–31G, and 6–31G* basis sets. Six distinct true minimum energy structures (named B, BB, BC, CROWN, TBC, and TCC₁), characterized through harmonic frequency analysis, were located on the multidimensional PES. Two transition state structures were also located on the PES for the cyclooctane molecule. Electron correlation effects were accounted for using the Møller–Plesset second-order perturbation theory (MP2) approach. The predicted global minimum energy structure on the ab initio PES for the cyclooctane molecule is the BC conformer. A gas phase electron diffraction study at 300 K suggested a conformational mixture while an NMR study in solution at 161.5 K predicted the BC conformer as the predominant form. The equilibrium constants reported in the present study, which were evaluated from the ab initio calculated total Gibbs free energy change values, were in good agreement with both experimental investigations. The ab initio results showed that the low temperature condition significantly favored the BC conformer while above room temperature both BC and CROWN structures can coexist. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 524–534, 1998     

Ab initio and DFT studies on hydrolyses of phosphorus halides

Hydrolyses of phosphorus halides, (RO)(2)POX where R = H or Me and X = F or Cl, in the gas phase and in the reaction field have been investigated theoretically with ab initio and the density functional theory (DFT). The free energy of activation in the reaction field was also estimated using the Onsager method with a correction of entropy change and basis set superposition error (BSSE). The reaction of (MeO)(2)POF proceeds through a path with bifunctional catalysis regardless of the medium, but the reaction of (MeO)(2)POCl proceeds through bifunctional and general base catalysis in the gas phase and in water, respectively. The estimated free energy barrier of 23 kcal/mol for the hydrolysis of (MeO)(2)POF is in good agreement with the experimental values of 24 kcal/mol, and relative barrier of 3 kcal/mol to the (MeO)(2)POCl is also in good agreement with the experimental values of 5 kcal/mol of diisopropyl phosphorus halides ((Pr(i)O)(2)POX, X = F and Cl).     

Calculation of correlation energy by many-body diagrammatic perturbation theory

The many-body diagrammatic perturbation theory is used for calculation of the correlation energy of closed-shell molecular systems. We apply Brueckner's concept of the two-particle renormalized interaction defined by a non-linear diagrammatic expression containing all possible (diagonal and/or non-diagonal) particle-particle, hole-hole and particle-hole intermediate elementary processes. Then, a “second-order” simple diagrammatic expression for the correlation energy can be formed, where the correlation energy is approximated by all the diagrams with biexcited intermediate states. An illustrative numerical application for the LiH molecule is presented. This article is dedicated to the memory of our friends and colleagues Dr. Jarka Surá and Dr. Marta Černayová, who tragically died in July 1976.     

ChemInform Abstract: Ab initio SCF Calculation on LinHm Molecules and Cations with Four or Less Atoms.

The stabilities and structures of all possible molecules and monopositive Cations of Li Clusters and hydrides up to four atoms have been analyzed.     

烟酸二聚体的结构与性质

用密度泛函理论B3LYP方法选取6-311+G(d,p)基组对烟酸-烟酸复合物进行了量子化学计算研究, 通过在相同水平下的频率振动分析发现了势能面上存在7个极小值点, 其最稳定构型1对应一N…H—O型强氢键, 其结合能在消除基组重叠误差后为-48.3 kJ·mol-1. 通过自然键轨道(NBO)分析, 研究了电荷转移及轨道相互作用. 通过自洽反应场(SCRF)理论中的极化连续模型(PCM)在介电常数分别为1.0(真空)、2.247(苯)、10.36(二氯乙烷)、20.7(丙酮)、32.63(甲醇)、78.39(水)的不同溶剂环境下重新优化烟酸复合物势能面上最稳定构型1, 研究了溶剂对烟酸复合物几何构型及结合能的影响. 发现溶剂化作用增大了烟酸复合物分子间的结合能, 导致N…H距离减小. 当溶液介电常数在1.0-32.63范围时, 溶剂效应十分显著, 当介电常数大于32.63后, 溶剂化作用趋于稳定.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

Ab initio quasirelativistic calculations on electronic transitions in ICl by the multireference many‐body perturbation theory

A quasirelativistic perturbative method of ab initio calculations on ground and excited molecular electronic states and transition properties within the relativistic effective core potential approximation is presented and discussed. The method is based on the construction of a state‐selective many‐electron effective Hamiltonian in the model space spanned by an appropriate set of Slater determinants by means of the second‐order many‐body multireference perturbation theory. The neglect of effective spin–orbit interactions outside of the model space allows the exploitation of relatively high nonrelativistic symmetry during the evaluation of perturbative corrections and therefore dramatic reduction of the cost of computations without any contraction of the model‐space functions. One‐electron transition properties are evaluated via the perturbative construction of spin‐free transition density matrices. Illustrative calculations on the X0⁺ ? A1, B0⁺, and (ii)1 transitions in the ICl molecule are reported. © 2002 Wiley Periodicals, Inc. Int J Quantum Chem, 2002