硼、氮取代二聚苯中碳原子导致层间相互作用本质变化研究

在MP2/aug-cc-pVTZ水平下优化了二聚苯及其碳被硼、氮原子取代的几何构型;计算了这些二聚体的CCSD(T)/CBS相互作用能;并用SAPT2+/aug-cc-pVDZ分析了相互作用能成分.探索了硼、氮取代二聚苯中碳原子所导致层间相互作用本质变化.结果表明:稳定构型而言,从苯(C_6H_6)、1-氮-2-硼杂苯(BNC_4H_6)、1,3-二氮-2,4-二硼杂苯(B_2N_2C_2H_6)二聚体的平行移位(PD)构型向无机苯(B_3N_3H_6)二聚体的夹心(S)构型转变,其中C_6H_6和B_3N_3H_6二聚体稳定构型的堆叠型式分别与之相对应的石墨烯(GE)和六方氮化硼(h-BN)2D层间材料的堆叠型式相一致.硼、氮原子取代二聚苯中的碳原子后使其相互作用能增大,其中BNC_4H_6和B_2N_2C_2H_6二聚体的相互作用能增大较为明显.所研究体系二聚体稳定构型均以色散能为主导、静电能次之、诱导能相对较小.硼、氮取代二聚苯中碳原子后其静电能对总吸引能的贡献明显增大.     

酞菁及其过渡金属配合物与富勒烯非共价相互作用的密度泛函理论方法研究

精确计算酞菁、卟啉及其过渡金属配合物与富勒烯分子间非共价相互作用能,探讨其相互作用本质,需要完成由上百个原子构成的超分子复合物、涉及数千个基函数第一性原理计算,对于当前的理论方法和计算机发展水平都构成严峻的挑战.本文首先选择酞菁…C₆₀、酞菁镍…C₆₀和酞菁锌…C₆₀作为此类非共价复合物的代表.在保留它们分子间相互作用区域结构特征前提下,裁剪成四氮杂卟啉…碗烯、四氮杂卟啉镍…碗烯和四氮杂卟啉锌…碗烯三个模型.通过大规模的CCSD(T)/CBS计算,获得了此三个模型的达到化学精度的分子间相互作用能.以此为基准,从27种密度泛函及其色散和非局域校正方法中,筛选出ωB97M-V,ωB97X-D和BLYP-D3BJ方法,其结果与CCSD(T)/CBS基准值相差均在化学精度范围,特别是ωB97M-V方法,最大误差仅为0.05 kcal/mol,很好地重现了CCSD(T)/CBS的计算结果,说明ωB97M-V方法不仅适用于主族元素分子间非共价相互作用的研究,对于涉及Ni, Zn等闭壳层过渡金属配合物与富勒烯相互作用能计算同样是非常优...     

TATB分子的结构及分子间相互作用

采用Hartree-Fock方法和密度泛函BPW91方法,对TATB双分子系统的几何结构和能量进行了优化和计算,得到双分子共平面的平衡结构.结果表明在TATB分子内部H和O之间存在着明显的氢键.分子之间的相互作用使TATB分子的对称性下降.一般地,没有充分计入分子之间的相互作用的影响时,计算得到的C-N键长比实验测得值大,而计算得到的C-C键长比实验测得值小.计算并讨论了TATB分子之间的相互作用能,表明TATB分子形成晶体结构时,应当以层状晶体结构最为稳定.     

O3与O3分子间相互作用的理论研究

采用从头算方法,在B3LYP水平下对O3分子单体及其二聚体进行了几何结构优化,得到四种O3分子二聚体的稳定构型,并经MP2法、均衡法(Massage)分别校正了电子相关和基组叠加误差(BSSE),通过对分子中原子间的最短距离,分子间质心距离与各构型相互作用能的关系进行了分析、比较,得到了最稳定的构型,其对应的O3分子二聚体相互作用能为:-8.788 kJ/mol.     

苯与PtH分子间电子相互作用的理论研究

采用从头算法在B3LYP水平分别对苯和PtH分子单体及其分子复合物体系(苯-PtH)进行了几何结构优化和振动频率计算,得到了3种苯-PtH分子复合物构型,研究了苯与PtH分子间相互作用对Pt-H键振动频率的影响,并得出了苯与PtH分子之间的相互作用能.     

水二聚体分子间相互作用势

利用超分子二阶到四阶多体微扰理论和扩展的相关一致基组(aug-cc-pVTZ)结合有效的中点键函数(3s3p2d1f1g)计算水二聚物的平衡结构和分子间相互作用势,并用平衡方法修正基组重叠误差.在MP2/aug-cc-pVTZ理论级别优化水二聚物几何构型.与正常优化的结果相比,平衡修正优化得到的RO-O和α值分别轻微的增加0.002 nm和0.19°,同时,θ值减小0.013°.在MP2水平利用扩展的相关一致基组结合有效键函数,预言了RO-O和ΔECP值分别为0.0923 nm和-4.86 kcal/mol,计算结果与实验值符合得很好.用exp-4.2势函数拟合分子间相互作用能的离散点,拟合结果与从头算计算的结果一致.     

甲醛与甲酰胺相互作用的从头算研究

在MP2 /6 3 1G(d)和MP2 (FC) /6 3 11+ +G(d ,p)水平上 ,对H2 CO和HCONH2 以及设计的 3种构型H2 CO…HCONH2 复合物等进行几何全优化计算 ,经振动频率分析 ,确认它们为势能超曲面上的稳定驻点 .然后在MP2 /6 3 11+ +G(2df,3p)水平上进行单点能计算和基组重叠误差 (BSSE)校正以获得相互作用能 .并利用自然键轨道理论和分子中的原子理论探讨H2 CO和HCONH2 相互作用的本质 .分子间相互作用的能量分解分析显示 ,静电能在H2 CO…HCONH2 相互作用能中占主导地位     

甲醛与乙腈相互作用的理论研究

在HF/ 6 311G(d ,p)、MP2 / 6 311G(d ,p)和B3LYP/ 6 311G(d ,p)水平上 ,对H2 CO和CH3 CN以及设计的 4种结构H2 CO…CH3 CN复合物等进行几何全优化和振动频率计算 ,排除振动频率为负值的非局域极小点结构 ,并对稳定的环状构型复合物结合能进行基组重叠误差校正和零点振动能校正 .分子间相互作用的能量分解分析显示 ,静电能在H2 CO…CH3 CN相互作用能量中占主导地位 ,电荷转移能居第二位 .     

弱结合分子Kr-HF结构与相关效应

采用二阶、三阶MollerPlesset微扰理论方法(MP2,MP3),组态相互作用方法(QCISD)在不同的基组下对弱结合分子体系KrHF进行了abinitio计算,得到了KrHF体系的两个不同的线型平衡几何结构:KrHF和KrFH.对于KrHF分子,在MP2/6311++G水平上得到Kr与HF分子中心的间距为037787nm,离解能为61480eV,谐振频率分别为ν1(σ)=41879638cm-1,ν2(π)=1622953cm-1,ν3(σ)=418689cm-1.并计算得到了这两个分子构型的热力 关键词：     

HMX与含硼化合物相互作用的理论计算

在PBX体系中加入合适的键合剂通常能改善HMX与高聚物的界面作用.为了预测其界面作用强度及粘合效果,本文采用B3LYP密度泛函方法和6-31G*基组计算HMX与含硼化合物BR1R2R3的相互作用.优化得到了存在B…O相互作用的HMX…BR1R2R3复合物的几何构型.计算结果表明,HMX…BH3和HMX…BH2(CH3)中的B与O原子间相互作用强,二者分子间相互作用能在MP2/6-31 G**//B3LYP/6-31G*水平上分别为-65.3 kJ/mol和-40.2 kJ/mol.由于空间位阻效应HMX中的O原子与BH(CH3)2、B(CH3)3和B(OH)(CH3)(CH2NH2)中的B原子距离大于3.理论计算结果与实验结果基本一致.     

迭代的多组态相互作用方法

本文提出了迭代的组态相互作用方法(IMRCI). IMRCI被用来计算H₂O和CH₂(单重态和三重态)在平衡态和远离平衡时的电子能量. IMRCI、MP2、MP3、MP4、CCSD和CCSD(T)还用来计算H₂O、CH₂(单重态和双重态)和N₂的势能曲线. 这些计算结果表明IMRCI的结果不依赖初始的多参考态组态函数,并能较快地收敛到完全组态相互作用的计算结果. 相比完全组态相互作用的结果,仅需2至4次迭代,IMRCI的误差就能达到10^-5 hartree数量级. 另外,IMRCI还提供了寻找势能面上主要电子组态的一个有效途径. 这将有助于单参考态和多参考态理论模型得到准确的计算结果.     

Ne-HF体系的相互作用势及散射截面的密耦计算

利用非线性最小二乘法拟合在CCSD(T)/aug-cc-pVQZ理论水平下计算的相互作用能,得到了基态Ne-HF体系相互作用势的解析表达式.基于拟合的CCSD(T)势,通过密耦计算得到了入射能量分别为60,75,100和150meV 时,Ne-HF散射的微分截面和分波截面,详细讨论了散射截面随能量的变化趋势以及态-态激发截面对总非弹性散射截面的影响. 关键词： 相互作用势 散射截面 密耦计算 Ne-HF体系     

The keto-enol equilibrium in substituted acetaldehydes: focal-point analysis and ab initio limit

High-level ab initio electronic structure calculations up to the CCSD(T) theory level, including extrapolations to the complete basis set (CBS) limit, resulted in high precision energetics of the tautomeric equilibrium in 2-substituted acetaldehydes (XH₂C-CHO). The CCSD(T)/CBS relative energies of the tautomers were estimated using CCSD(T)/aug-cc-pVTZ, MP3/aug-cc-pVQZ, and MP2/aug-cc-pV5Z calculations with MP2/aug-cc-pVTZ geometries. The relative enol (XHC?=?CHOH) stabilities (ΔE _{e,CCSD(T)/CBS}) were found to be 5.98?±?0.17, ?1.67?±?0.82, 7.64?±?0.21, 8.39?±?0.31, 2.82?±?0.52, 10.27?±?0.39, 9.12?±?0.18, 5.47?±?0.53, 7.50?±?0.43, 10.12?±?0.51, 8.49?±?0.33, and 6.19?±?0.18?kcal?mol^?1 for X?=?BeH, BH₂, CH₃, Cl, CN, F, H, NC, NH₂, OCH₃, OH, and SH, respectively. Inconsistencies between the results of complex/composite energy computations methods Gn/CBS (G2, G3, CBS-4M, and CBS-QB3) and high-level ab initio methods (CCSD(T)/CBS and MP2/CBS) were found. DFT/aug-cc-pVTZ results with B3LYP, PBE0 (PBE1PBE), TPSS, and BMK density functionals were close to the CCSD(T)/CBS levels (MAD?=?1.04?kcal?mol^?1).     

Ab initio potential energy surface and intermolecular vibrational frequency of C3-He complex

The intermolecular potential energy surface for C₃–He complex has been constructed using supermolecular CCSD(T) and MP4 methods. The potential surfaces have been calculated for 27 values of R ranging from 2.8 to 8.0 Å and 19 values of θ equally spaced between 0° and 180°. Both CCSD(T) and MP4 potentials have similar global behaviors. The global minimum in each of the potentials corresponds to the slightly distorted T-shaped geometry. On the basis of these two potentials, the intermolecular vibrational energies and wavefunctions were calculated. The energy level pattern of the vdW vibrational states was predicted for C₃–He complex. The zero point bending motion of this complex has a range of 180°. The calculated fundamental frequency of vdW bending is 3.16 cm^?1 at CCSD(T) level, and 5.38 cm^?1 at the MP4 level. In addition, we have also constructed the intermolecular potential energy surface with C₃ bending coordinate of 160° by using supermolecular CCSD(T) method. Two local minima including arrow-shaped and Y-shaped configurations were determined. The rotational constants of three C₃–He structures including T-shaped, arrow-shaped and Y-shaped configurations at CCSD(T) level were also reported.     

A high level computational study of the CH4/CF4 dimer: how does it compare with the CH4/CH4 and CF4/CF4 dimers?

The interaction within the methane–methane (CH₄/CH₄), perfluoromethane–perfluoromethane (CF₄/CF₄) methane–perfluoromethane dimers (CH₄/CF₄) was calculated using the Hartree–Fock (HF) method, multiple orders of Møller–Plesset perturbation theory [MP2, MP3, MP4(DQ), MP4(SDQ), MP4(SDTQ)], and coupled cluster theory [CCSD, CCSD(T)], as well as the PW91, B97D, and M06-2X density functional theory (DFT) functionals. The basis sets of Dunning and coworkers (aug-cc-pVxZ, x?=?D, T, Q), Krishnan and coworkers [6-311++G(d,p), 6-311++G(2d,2p)], and Tsuzuki and coworkers [aug(df, pd)-6-311G(d,p)] were used. Basis set superposition error (BSSE) was corrected via the counterpoise method in all cases. Interaction energies obtained with the MP2 method do not fit with the experimental finding that the methane–perfluoromethane system phase separates at 94.5?K. It was not until the CCSD(T) method was considered that the interaction energy of the methane–perfluoromethane dimer (?0.69?kcal?mol^?1) was found to be intermediate between the methane (?0.51?kcal?mol^?1) and perfluoromethane (?0.78?kcal?mol^?1) dimers. This suggests that a perfluoromethane molecule interacts preferentially with another perfluoromethane (by about 0.09?kcal?mol^?1) than with a methane molecule. At temperatures much lower than the CH₄/CF₄ critical solution temperature of 94.5?K, this energy difference becomes significant and leads perfluoromethane molecules to associate with themselves, forming a phase separation. The DFT functionals yielded erratic results for the three dimers. Further development of DFT is needed in order to model dispersion interactions in hydrocarbon/perfluorocarbon systems.     

The electron affinity of Al13H cluster: high level ab initio study

Al₁₃H clusters have been considered candidates for cluster assembled materials. Here we have carried out benchmark calculations for the Al₁₃H cluster, both neutral and anionic, with the aim of verifying the nature of stationary points on the potential energy surface, studying dynamics of H atom and determining an adiabatic electron affinity. A range of correlated methods applied include second-order perturbation theory (MP2), spin-component-scaled MP2, coupled electron pair (CEPA) and coupled cluster singles and doubles with perturbative triple corrections (CCSD(T)). These methods are used in combination with the correlation consistent basis sets through aug-cc-pVTZ including extrapolation to the complete basis set (CBS) limit. Performance of several different flavours of density functional theory (DFT) such as generalised gradient approximation (GGA), hybrid GGA, meta-GGA and hybrid-meta-GGA is assessed with respect to the ab initio correlated reference data. The harmonic force constant analysis is systematically performed with the MP2 and DFT methods. The MP2 results show that for neutral Al₁₃H only the hollow structure is a potential energy minimum, with the bridged structure being a transition state for the H shift from the hollow site to the adjacent hollow site. The CCSD(T)/aug-cc-pVTZ (CCSD(T)/CBS) estimate of the energy barrier to this H shift is 2.6 (2.9) kcal/mol, implying that the H atom movement over the Al₁₃H cluster surface is facile. By contrast, the DFT force constant analysis results suggest additional terminal and bridged minima structures. For the anion Al₁₃H^?, exhibiting ‘stiffer’ potential energy surface compared to the neutral, the existence of the hollow and terminal isomers is consistent with the earlier photoelectron spectroscopy assignment. The adiabatic electron affinity of Al₁₃H is determined to be 2.00 and 1.95 eV (the latter including the ΔZPE correction) based on the CCSD(T) energies extrapolated to the CBS limit, whereas the respective CCSD(T)/CBS thermodynamic EA values are 2.79 and 2.80 eV.     

碱金属Li原子与CO分子体系相互作用势从头计算研究

采用三重激发项校正CCSD的CCSD(T)方法和AUG-CC-PVnZ(n=2,3,4)基组,优化了CO基态分子结构,并计算了碱金属Li原子与CO分子的相互作用势,共1010个构型势能点得到体系的势能面.结果表明:同一方法下,不同基组得到的CO基态分子的键长、能量等均与实验符合很好.Li-CO势能面体现较小的各向异性势,存在两个势阱,且都为非严格T型结构;基组重叠误差(BSSE)对相互作用势的影响比较明显,采用CP方法(Counterposie method)消除基组重叠误差,不同基组计算的相互作用势显示较好一致.体系各项异性势势阱值远高于CO基态分子转动常数,使得碰撞产生强烈的非弹性碰撞,这将阻止协同冷却制备超冷分子.     

Enthalpy difference between conformations of normal alkanes: effects of basis set and chain length on intramolecular basis set superposition error

The quantum chemistry of conformation equilibrium is a field where great accuracy (better than 100?cal?mol^?1) is needed because the energy difference between molecular conformers rarely exceeds 1000–3000?cal?mol^?1. The conformation equilibrium of straight-chain (normal) alkanes is of particular interest and importance for modern chemistry. In this paper, an extra error source for high-quality ab initio (first principles) and DFT calculations of the conformation equilibrium of normal alkanes, namely the intramolecular basis set superposition error (BSSE), is discussed. In contrast to out-of-plane vibrations in benzene molecules, diffuse functions on carbon and hydrogen atoms were found to greatly reduce the relative BSSE of n-alkanes. The corrections due to the intramolecular BSSE were found to be almost identical for the MP2, MP4, and CCSD(T) levels of theory. Their cancelation is expected when CCSD(T)/CBS (CBS, complete basis set) energies are evaluated by addition schemes. For larger normal alkanes (N?>?12), the magnitude of the BSSE correction was found to be up to three times larger than the relative stability of the conformer; in this case, the basis set superposition error led to a two orders of magnitude difference in conformer abundance. No error cancelation due to the basis set superposition was found. A comparison with amino acid, peptide, and protein data was provided.