Ar3H＋中多体势能拆分的从头计算分析

使用从头算的方法计算了Ar3H＋的稳定构型及红外振动基频，并与Ar2H＋的稳定结构及红外振动频率作了对比，讨论了二者之间的联系与差别.对Ar3H＋一个特殊构型下的势能线作了扫描，在此势能线下将Ar3H＋的四体势能拆解为两体，半三体及三体势能的加和，初步探讨了高阶势能分解为低阶势能的最佳途径.     

β-三氰基乙烯取代杯[4]吡咯与阴离子相互作用的密度泛函理论研究

采用密度泛函B3LYP方法对β-三氰基乙烯取代杯[4]吡咯主体分子及其与阴离子(F-、Cl-、CH3COO-、H2PO4-)形成的复合物进行研究。结果表明,β-三氰基乙烯取代杯[4]吡咯可与阴离子以分子间氢键相互作用形成复合物,并且其与阴离子结合能力大小的顺序为F-> CH3COO-> H2PO4-> Cl-,与实验结果基本一致;通过与杯[4]吡咯对比可见,强吸电子取代基的引入增强了主体分子对阴离子的结合能力。本文还从几何构型、振动光谱、NBO分析及前线轨道等方面来阐述β-三氰基乙烯取代杯[4]吡咯与不同阴离子氢键相互作用的本质以及吸电子取代基的引入对杯[4]吡咯与F-和Cl-之间主-客体相互作用的影响。     

C12簇及其碳管烷C12H12的理论研究

C60和C70等碳笼烯的发现及其在新物质、新材料等研究领域中的重要性，激励人们去探索合成更多新型多面体碳簇化合物与碳元素类似物[1，2]．最近，两类新型多面体碳管笼烯的设想提出来了[3，4]应用HMO和MNDO方法对其稳定性变化规律、结构和成键特征进行了讨论.关于多苯的vanderWaals簇实验上已有广泛的调查[5，6]，理论上对其二聚物（C6H6）。的不同几何构型与稳定性进行了深入的研究[7]．本文采用abinitio计算有效势（effectivecorePotential）方法，对C12碳管元素簇和C12H12，碳管烷的平衡结构、稳定性和价键特征进行了理论预测．1…     

硼氢及客体二十面体簇合物的结构和稳定性

在(n=0、－1、－2、－3、－4)簇合物几何构型及稳定性研究的基础上,进一步对它的各种内含式和外 接式二十面体簇合物(X@B12H122－和XB12H122－,X=H0/＋、Li0/＋、He、Ne、Be0/2＋、Na＋、Mg2＋)进行了优化和 计 算.发现在内含式结构X@ B12H122－中,当X=Li＋、Be2＋、Mg2＋时,构型较稳定;在外接式结构中, XB12H122－(C3v)结构比XB12H122－(C2v)的结构稳定.通过IRC计算,确定XB12H122－(C2v)是X与B12H122－作用生成产物XB12H122－(C3v)的一种过渡态.     

氯代甲苯双电荷离子的单分子解离反应研究

研究了在70 eV电子轰击电离条件下，氯代甲苯及氯化苄产生的双电荷离子[C_7H_7Cl]~(2＋)、[C_7H_6Cl]~(2＋·)和[C_7H_5Cl]~(2＋)为母体的两种类型单分子解离反应．主要讨论了亚稳双电荷离子的异构化反应、失H解高的“偶电子规则”以及单分子电荷分离过渡态的结构．     

N_2H_4-CH_3OH氢键团簇体系的从头计算

用从头计算法研究了 (N2 H4-CH3OH)氢键团簇体系。分别在HF/6 31G 和HF/6 31G 水平上对它们的中性和离子团簇进行几何全优化 ,得到了 3种中性混合团簇稳定构型和离子混合团簇稳定构型 ,并对其能量和稳定性进行了比较。讨论了 3种不同构型离子团簇可能的解离通道。给出了质子化混合团簇的稳定构型 ,并对其可能的解离通道进行了讨论。文中最后计算出N2 H4,CH3OH ,(N2 H4-CH3OH)团簇的质子亲和能 (PA) ,分别为 :2 0 6.7kcal/mol,1 78.3kcal/mol,2 2 7.5kcal/mol,其中质子亲和能PAcalc[N2 H4]与实验值PAexp[N2 H4]=2 0 4 .8kcal/mol符合得很好。     

气相中CrO2^＋主活化甲烷C—H键的理论研究

用密度泛函UB3LYP／6-311＋＋G＊＊方法计算研究了气相中CrO2^＋（2^A1/4^A＂）活化甲烷C—H键的微观机理，找到了四条反应通道．对其中涉及的两态反应（TSR）进行了分析，并对影响反应机理和反应速率的势能面交叉现象（potential energy surfaces crossing）进行了详细讨论，进而运用Hammond假设和Yoshizawa等的内禀坐标单点垂直激发计算的方法找出了一系列势能面交叉点[crossing points（CPs）]，并作了相应的讨论．进一步用碎片分子轨道理论[fragment molecular orbital（FMO）]对TS1中的轨道相互作用进行了分析，解释了CrO2^＋活化甲烷C—H键的机理．     

Sc3C2@C80与[12]CPP纳米环之间非共价相互作用的理论研究

金属富勒烯嵌套于纳米环内形成主客体系, 二者产生的主客体作用可诱导内部金属团簇的取向, 影响分子的电子结构等性质. 本文基于密度泛函理论(DFT)计算, 对碳纳米环[12]CPP(CPP=环苯撑, 主体分子)与金属富勒烯Sc₃C₂@C₈₀(客体分子)形成的主客体配合物的结构和性质进行了研究. 计算结果表明, 在最稳定构型中, [12]CPP呈现椭圆形, Sc₃C₂@C₈₀与[12]CPP的质心不再重合. Sc₃C₂@C₈₀在[12]CPP内旋转对构型总体能量影响仅为13.51 kJ/mol. [12]CPP向Sc₃C₂@C₈₀转移了0.03 e, 主客体分子之间存在弱相互作用. 对二者相互作用的分析结果表明, 色散作用在弱相互作用中占主导地位.     

合成文拉法星的2-氮杂-1,3-丁二烯杂Diels-Alder反应与其 单分子环合竞争反应的研究

研究了不同的Lewis酸催化剂、温度、微波等条件下, 4-(4-甲氧苯基)-1-苯基-3-三甲基硅氧基-2-氮杂-1,3-丁二烯(2-ABDE)和亲二烯体环己酮发生杂Diels-Alder反应生成[4＋2]的六元环合产物噁嗪酮衍生物, 伴随的2-ABDE的 [2＋2]单分子环合, 生成单环β-内酰胺衍生物. 结果表明: 在低温条件下如(－78 ℃)[4＋2]反应占主导; 而在高温条件下(如135 ℃)仅进行[2＋2]反应. 微波加热方式可显著提高[2＋2]反应的速率和产率. 不同的Lewis酸催化剂对[2＋2]反应和[4＋2]反应的催化效率不同. Lewis酸的酸性强弱、软硬对2-ABDE的[2＋2]反应的催化能力起决定性作用.     

半胱氨酸阴离子与咪唑阳离子间相互作用的理论研究

采用密度泛函理论在B3LYP/6-311＋G(d,p)水平上对1-乙基-3-甲基咪唑阳离子和半胱氨酸阴离子形成的气态阴阳离子对([Emim][Cys])进行理论研究. 通过几何结构优化以及频率分析得到势能面上7个稳定的离子对构型. 计算结果表明[Emim]＋和[Cys]－之间存在较强的氢键相互作用, 其稳定化能主要来源于[Cys]－中羰基O的孤对电子lp(O) 和[Emim]＋中C—H反键轨道 s*(C—H) 之间的相互作用, lp(O)®s*(C—H). [Emim][Cys]_S1是最稳定的离子对构型, 考虑BSSE的相互作用能为－387.66 kJ/mol. 从NPA和NBO分析以及AIM (Atoms in Molecules)计算等方面阐述了半胱氨酸阴离子与咪唑阳离子之间氢键相互作用的本质, 并初步探讨了阴阳离子对相互作用对氨基酸离子液体性质的影响.     

Energetics of proton transfer in liquid water. I. Ab initio study for origin of many-body interaction and potential energy surfaces

The energetics of proton transfer in liquid water investigated by using ab initio calculation. The molecular electronic interaction of hydrated proton clusters in classified into many-body interaction elements by a new energy decomposition method. It is found that up to three-body molecular interaction is essential to describe the potential energy surface. The three-body effect mainly arises from the (non-classical) charge transfer and strongly depends on their configuration. Higher than three-body effects are small enough to be neglected. To simulate the liquid state reactions, two cluster models including all water molecules up to the second shell in the proton transfer reactions are employed. It is shown that these proton transfer reactions only involve small potential energy barriers of a few kcal/mol or less when structural rearrangement of the solvent is induced along the proton movement.     

Four-body interaction energy for compressed solid krypton from quantum theory

The importance of the four-body contribution in compressed solid krypton was first evaluated using the many-body expansion method and the coupled cluster theory with full single and double excitations plus perturbative treatment of triples. All different four-atom clusters existing in the first- and second-nearest neighbor shells of face-centered cubic krypton were considered, and both self-consistent-field Hartree-Fock and correlation parts of the four-body interaction were accurately determined from the ambient conditions up to eightfold volume compression. We find that the four-body interaction energy is negative at compression ratio lower than 2, where the dispersive forces play a dominant role. With increasing the compression, the four-body contribution becomes repulsive and significantly cancels the over-softening effects of the three-body potential. The obtained equation of state (EOS) was compared with the experiments and the density-functional theory calculations. It shows that combination of the four-body effects with two- and three-body interactions leads to an excellent agreement with EOS measurements throughout the whole experimental range 0-130 GPa, and extends the prediction to 300 GPa.     

Quantum vibrational analysis and infrared spectra of microhydrated sodium ions using an ab initio potential

We present a full-dimensional potential energy surface and a dipole moment surface (DMS) for hydrated sodium ion. These surfaces are based on an n-body expansion for both the potential energy and the dipole moment, truncated at the two-body level for the H(2)O-Na(+) interaction and also for the DMS. The water-water interaction is truncated at the three-body level. The new full-dimensional two-body H(2)O-Na(+) potential is a fit to roughly 20,000 coupled-cluster single double (triple)/aug-cc-pVTZ energies. Properties of this two-body potential and the potential describing (H(2)O)(n)Na(+) clusters, with n up to 4 are given. We then report anharmonic, coupled vibrational calculations with the "local-monomer model" to obtain infrared spectra and also 0 K radial distribution functions for these clusters. Some comparisons are made with the recent infrared predissociation spectroscopy experiments of Miller and Lisy [J. Am. Chem. Soc. 130, 15381 (2008).].     

Ab initio copper–water interaction potential for the simulation of aqueous solutions

A new ab initio effective two-body potential that aims at mimicking the average copper–water interaction energy of the first solvation shell was developed. This new potential, together with the MCY water–water potential and a three-body ion–water–water induction potential, is tested in simulations of gas-phase clusters [Cu²⁺? (H₂O)₂₀] and diluted solutions [Cu²⁺? (H₂O)₂₀₀] at T = 298 K. The results of simulations with conventional ab initio pair potentials, with and without three-body induction corrections, are also presented. The different types of copper–water interaction potentials are evaluated comparatively and the efficiency of the newly proposed effective pair potential is discussed. © 1993 John Wiley & Sons, Inc.     

Sensitivity of the thermal and acoustic virial coefficients of argon to the argon interaction potential

Second, third, and fourth thermal and acoustic virial coefficients between 100 and 1000 K are computed for different argon interaction models derived from combinations of accurate two- and three-body potentials. Differences between the various interaction models tested mirror the presumed order in the accuracy of these models, but are not well captured at the level of the lowest-order contributions in the virial expansion: While the second- and third-order virial coefficients are found to be rather insensitive to small variations in the two- and three-body potentials, more pronounced differences in higher-order coefficients are currently of limited use in assessing the accuracy of the interaction potential due to difficulties in the unambiguous experimental determination of these higher-order coefficients. In contrast, pressure-volume and speed-of-sound data - both of which are experimentally known to highest accuracies - are found to be insensitive to small variations in the interaction model. All but the least accurate models reproduce experimental pressure-volume and speed-of-sound data near-quantitatively in regions where the (fourth-order) virial expansions apply. All quantities considered are found to be completely unaffected by a non-vanishing quadruple-dipole four-body potential.     

Calculation of argon trimer rovibrational spectrum

Rovibrational spectra of Ar3 are computed for total angular momenta up to J=6 using row-orthonormal hyperspherical coordinates and an expansion of the wave function on hyperspherical harmonics. The sensitivity of the spectra to the two-body potential and to the three-body corrections is analyzed. First, the best available semiempirical pair potential (HFDID1) is compared with our recent ab initio two-body potential. The ab initio vibrational energies are typically 1-2 cm-1 higher than the semiempirical ones, which is related to the slightly larger dissociation energy of the semiempirical potential. Then, the Axilrod-Teller asymptotic expansion of the three-body correction is compared with our newly developed ab initio three-body potential. The difference is found smaller than 0.3 cm-1. In addition, we define approximate quantum numbers to describe the vibration and rotation of the system. The vibration is represented by a hyper-radial mode and a two-degree-of-freedom hyperangular mode, including a vibrational angular momentum defined in an Eckart frame. The rotation is described by the total angular momentum quantum number, its projection on the axis perpendicular to the molecular plane, and a hyperangular internal momentum quantum number, related to the vibrational angular momentum by a transformation between Eckart and principal-axes-of-inertia frames. These quantum numbers provide a qualitative understanding of the spectra and, in particular, of the impact of the nuclear permutational symmetry of the system (bosonic with zero nuclear spin). Rotational constants are extracted from the spectra and are shown to be accurate only for the ground hyperangular mode.     

醇及脱氧核糖与水分子间三体作用强度的快速准确计算

快速准确预测醇及脱氧核糖分子与水形成的氢键复合物的三体作用强度, 对准确模拟水环境下蛋白质和DNA的结构和功能至关重要. 基于对多体极化作用的理解, 在可极化偶极-偶极作用模型(PBFF)基础上, 将体系中的极性化学键视为化学键偶极, 通过模拟键偶极的极化计算了醇及脱氧核糖与水分子形成的氢键复合物的三体作用能. 通过拟合甲醇与水氢键复合物的三体作用能随分子间距离变化的能量曲线确定了所需的参数. 将模型和所确定的参数应用于计算更多的甲醇、 乙醇及脱氧核糖与水氢键复合物的三体作用能, 检验了模型的准确性和参数的可转移性. 计算结果表明, 可极化偶极-偶极作用模型及所确定的参数能够较好地预测具有不同结构的氢键复合物的三体作用强度, 其精度可与MP2方法的计算精度相当.     

An empirical potential energy surface for the He-Br2 (B3pi(u)) van der Waals complex

An empirical potential energy surface is proposed for the He-Br2 (B3pi(u)) complex. The intermolecular potential is modeled as a sum of pairwise He-Br Morse interactions plus a three-body interaction term. The parameters of the potential are fitted in order to reproduce the spectral blue-shifts and vibrational predissociation line widths measured for He-79Br2 (B, v') in the range v' = 8-48 of Br2 vibrational excitations. The calculated blue-shifts and line widths are in very good agreement with the measurements (typically within experimental error or close to its limits) along the whole range of v' levels studied. It is particularly remarkable to note the accuracy provided by the interaction surface in the region of high v' excitations (v' > 35), where three-body effects become important. The behavior of the potential surface with the Br-Br separation is analyzed and correlated with the experimental findings.     

水团簇(H2O)n中多体作用强度与水分子间距的关系

使用高精度从头算方法(含基组重叠误差校正)计算了水团簇(H₂O)_n (n = 8, 10, 16, 20, 22, 24)中的所有二体、三体和四体作用能,分析了水团簇中的多体效应.研究表明,二体作用对体系总作用能的贡献高达70%以上,三体作用对总作用能的贡献可高达25%,四体作用在总作用能中所占比例不超过3%,五体及以上多体作用能在总作用能中所占比例更小,不超过0.5%.本文研究还表明,两个水分子间距小于0.68 nm的二体作用、三个和四个水分子中最近的两个水分子间距小于0.31 nm的三体和四体作用对体系总作用能的贡献高达99.4%.因此,以生物体系为对象的分子模拟方法应该具备准确地模拟两个水分子间距小于0.68 nm的二体作用、三个和四个水分子中最近的两个分子间距小于0.31 nm的三体和四体作用的能力.     

Diffusion Monte Carlo study of the equation of state of solid ortho-D2

The authors present results of diffusion Monte Carlo calculations for a system of solid ortho-D2 at different densities, for pressure ranging from 0 up to 350 MPa. They compare the equation of state obtained using two of the most used effective intermolecular potentials, i.e., the Silvera-Goldman and the Buck potentials, with experimental data, in order to assess the validity of the model interactions. The Silvera-Goldman potential has been found to provide a satisfactory agreement with experimental results, showing that, as opposed to what recently found for p-H2, three-body forces can be efficiently accounted for by an effective two-body term.