线型碳链LiC_(2n)Li的结构和电子光谱的密度泛函理论研究

应用密度泛函理论,在B3LYP/6 31G(d)水平上优化得到了线型簇合物LiC2nLi(n=1 ~10,D∞h)的基态平衡几何构型,并计算了它们的谐振动频率.利用含时密度泛函理论,计算了簇合物LiC2nLi的X1∑ +g→11∑ +u跃迁的垂直激发能,以及相应的振子强度.基于计算结果,建立了跃迁能和体系大小n的解析关系式. 同时也计算了体系的第一绝热电离能,讨论了体系的电离能与体系大小n的关系.     

乙硼烷离子和自由基异构体重排与体系势能面的能量子化学计算研究

在B3LYP/6-311G(d,p)和QCISD(T)/6-311++G(3df,2p)(单点)水平下, 对B2H+5阳离子和B2H*5自由基全优化得到9个几何异构体: B2H+5单态体系(D3h, C1), B2H+5三重态(Cs, Cs, C1), B2H*5自由基(C2v, Cs, Cs, Cs). 得到势能面上与体系异构化过程相联系的5种过渡态.     

吡咯-HCN体系在气相及溶液中相互作用的理论研究

用量子化学B3LYP方法在6-311G(d, p)水平上优化了吡咯-HCN氢键复合物,通过振动频率分析确定了两个吡咯-HCN体系稳定构型.为了得到更加精确的氢键作用能,采用相关一致基组aug-cc-pVDZ以及Boys 和Bernardi的CP(counterpoise)校正方法消除基组重叠误差后得到C－H…π和N－H…N型复合物的氢键相互作用能.为了确定B3LYP方法计算的相互作用能的可靠性,在MP2/aug-cc-pVDZ水平计算了复合物的氢键相互作用能,结果分别为－25.10和－19.30 kJ·mol－1.采用自然键轨道(NBO)分析考察了吡咯与HCN分子间轨道相互作用.以自洽场理论(SCRF)中的Onsager模型研究了不同极性溶剂对吡咯-氰化氢体系N－H…N型氢键几何构型,频率位移,电荷分布以及相对能量的影响.研究发现,当溶液的介电常数在1.5～30.0范围时,溶液作用十分显著,而当介电常数超过30.0以后,溶液作用已经达到了极限.     

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

在MP2／6—31G(d)和MP2(FC)／6—311 G(d，p)水平上，对H2CO和HCOOH以及设计的4种构型H2CO…HCOOH复合物等进行几何全优化计算，经振动频率分析，确认它们为势能面上的稳定驻点．然后在MP2／6—311 G(2df，2p)水平上进行单点能计算和基组重叠误差(BSSE)校正以获得相互作用能，并利用自然键轨道理论探讨H2CO和H(X)OH相互作用的本质。     

水和甲醇混合溶剂氢键相互作用及对高分子链溶解能力的理论研究

用密度泛函理论对水和甲醇混合溶剂体系的氢键结构进行了详细研究.通过构象和频率分析发现在水团簇中五聚体和六聚体环状结构最为稳定,同时发现一个全新的特征,即甲醇分子能与水五聚体和六聚体形成双氢键.根据各相互作用的稳定化能,分析了水和甲醇混合溶剂对PNIPAM溶解能力的影响,并对实验现象给予了合理解释.     

木犀草素与胞嘧啶相互作用的理论研究

采用密度泛函理论中的B3LYP方法,在6-31＋G＊基组水平上对木犀草素、胞嘧啶、木犀草素-胞嘧啶复合物进行结构优化和振动频率分析,得到了12种稳定复合物.并应用分子中的原子理论（AIM）分析、自然键轨道（NBO）理论分析得到复合物氢键性质和特征.通过基组重叠误差（BSSE）校正后的相互作用能、成键临界点电荷密度、二阶...     

第一性原理研究HC(NH2)2PbI3钙钛矿(三方相VS六方相)的结构及光电特性

采用第一性原理,研究三方及六方相FAPbI₃（FA=HC（NH₂）₂⁺）钙钛矿的结构及光电特性.结果表明,FAPbI₃钙钛矿由三方到六方的形变能够扭转PbI₃骨架,改变Pb-I键合特性,进而改变其禁带宽度值.两种晶体均属于直接带隙半导体,三方相FAPbI₃的直接带隙点位于布里渊区Z（0,0,0.5）对称点,具有较为理想的约1.50 eV的禁带宽度;六方相FAPbI₃的直接带隙点位于Γ（0,0,0）对称点,具有约2.50 eV的禁带宽度.FA离子不直接参与电子跃迁过程,仅仅充当电荷供体为PbI₃骨架提供超过0.7 e的电荷.相比于六方相晶体,三方相FAPbI₃具有更小的载流子有效质量,吸收光谱发生了显著的红移,并且其吸收特性优于六方相FAPbI₃和四方相MAPbI₃（MA=CH₃NH₃⁺）钙钛矿.六方相FAPbI₃比三方相晶系更稳定,FA和PbI₃骨架之间的结合作用力强于MA和PbI₃骨架之间的作用力.     

H_2O与NO,CN,OH自由基及负离子相互作用的CP研究

采用量子化学从头算方法在均衡 (Counterpoise)校正和非校正势能面上研究了H2 O与NO ,NO-,CN ,CN-,OH ,OH-之间的相互作用 .比较了CP 梯度优化和非校正梯度优化以及基组的选择对超分子结构和能量的影响 .研究表明 ,6 311++g 基组对于这些体系的研究有很高的效率 .这些自由基和负离子均能与H2 O形成强弱不同的氢键 .按相互作用的强弱次序依次为OH-,CN-,NO-,OH ,CN ,NO ,由相互作用能ΔECP,成键临界点电荷密度 ρ ,二阶稳定化相互作用能E( 2 )分析均得到同样结果 .CP 梯度优化和非校正梯度优化所得能量及BSSE相差不大     

Li/Li(100)面自扩散体系的理论研究

本文应用对势方法构造了Li/Li(100)表面缺陷体系的吸附扩散相互作用模型势,考察了各种台阶、扭结、空位等表面缺陷对锂原子表面吸附扩散行为的影响。根据缺陷表面体系结合能和表面扩散活化势垒的结果,提出和分析了捕获能力和捕获网链的概念。     

MRCI及其近似方法下C2分子X1∑+g, B1△g和B'1∑+g态的势能曲线

运用MRCI及其近似方法计算得到C2分子X1∑+g,B1△+g,和B'1∑+g 3个电子态的势能曲线,并以FCI方法为基准,考察比较各种方法的优缺点.MRCI计算与FCI计算几乎不存在差别,在FCI计算不现实的情况下,MRCI计算足以作为基准考验其他计算方法.多种近似方法均能达到化学精度,并提出一个有效的选择参考组态的方法,提高MRCI计算效率的同时,保证计算的精确性.     

Structural and Spectroscopic Properties of Linear Carbon Chains NC2nN and HC2n+1N (n = 1～10)

Using density functional theory, geometries and vibrational frequencies of linear chains NC2nN and HC2n+1N (n = 1～ 10) have been investigated.Time-dependent density functional theory (TD-DFT) has been used to calculate the vertical transition energies and oscillator strengths for the X 1∑+k → 11∑+u transition in NC2nN (n = 1 ～10) and X1∑+ → 11∑+ transition in HC2n+1N (n =1 ～7).On the basis of present calculations, the explicit expressions for the size dependence of the excitation energy and the first adiabatic ionization energy in both carbon chains have been suggested.     

Theoretical study of interaction of 2,6-dithiopurine with Li+, Na+, K+, Be2+, Mg2+, and Ca2+ metal cations

The geometries of the complexes of Li⁺, Na⁺, K⁺, Be²⁺, Mg²⁺, and Ca²⁺ metal cations with different possible 2,6-dithiopurine anions (DTP) were studied. The complexes were optimized at the B3LYP level and the 6-311++G(d, p) basis set. The interactions of the metal cations at different nucleophilic sites of various possible 2,6-dithiopurine anions were considered. It was revealed that metal cations would interact with 2,6-dithiopurine anions in a bicoordinate manner. In the gas phase, the most preferred position for the interaction of Li⁺, Na⁺, and K⁺ cations is between the N₃ and S₂ sites, while all divalent cations Be²⁺, Mg²⁺, and Ca²⁺ prefer binding between the N₇ and S₆ sites of the corresponding 2,6-dithiopurine. The influence of aqueous solvent on the relative stability of different complexes has been examined using the Tomasi’s polarized continuum model. The basis set superposition error (BSSE) corrected interaction energy was also computed for complexes. The AIM theory has been applied to analyze the properties of the bond critical points (electron densities and their Laplacians) involved in the coordination between 2,6-dithiopurine anions and the metal cations. It was revealed that aqueous solution would have significant effect on the relative stability of complexes obtained by the interaction of 2,6-dithiopurine anions with Mg²⁺ and Ca²⁺ cations. The effect of metal cations on different NH and CS stretching vibrational modes of 2,6-dithiopurine has also been discussed.     

Theoretical study on structures and stability of HC2P isomers

 The structures and isomerization pathways of various HC₂P isomers in both singlet and triplet states are investigated at the B3LYP/6-311G(d,p), QCISD/6-311G(d,p) (for isomers only) and single-point CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) levels. At the CCSD(T)/6-311G(d,p)//B3LYP/6-311G(d,p) level, the lowest-lying isomer is a linear HCCP structure ³ 1 in the ³∑⁻ state. The second low-lying isomer has a CPC ring with exocyclic CH bonding ¹ 5 in a singlet state at 10.5 kcal/mol. The following third and fourth low-lying isomers are a singlet bent HCCP structure ¹ 1 at 20.9 kcal/mol and a bent singlet HPCC structure ¹ 3 at 35.8 kcal/mol, respectively. Investigation of the HC₂P potential-energy surface indicates that in addition to the experimentally known isomer ³ 1, the other isomers ¹ 1, ¹ 3 and ¹ 5 also have considerable kinetic stability and may thus be observable. However, the singlet and triplet bent isomers HCPC ¹ 2 and ³ 2 as well as the triplet bent isomer HPCC ³ 3 are not only high-lying but are also kinetically unstable, in sharp contrast to the situation of the analogous HCNC and HNCC species that are both kinetically stable and that have been observed experimentally. Furthermore, the reactivity of various HC₂P isomers towards oxygen atoms is briefly discussed. The results presented here may be useful for future identification of the completely unknown yet kinetically stable HC₂P isomers ¹ 1, ¹ 3 and ¹ 5 either in the laboratory or in interstellar space. Received: 5 November 2000 / Accepted: 25 November 2001 / Published online: 8 April 2002     

Theoretical investigations on analytical potential energy function and spectroscopic parameters for the state b3Πu of dimer 7Li2

The SAC‐CI (symmetry‐adapted‐cluster configuration‐interaction) method presented in Gaussian 03 program package is applied to investigate the adiabatic potential energy curves (PECs) of ⁷Li₂(b³Π_u). These calculations are performed at numbers of basis sets, such as 6‐311++G(3df,3pd), 6‐311++G(2df,2pd), 6‐311++G(df,pd), D95V++, D95(3df,3pd), D95(d,p), cc‐PVTZ, 6‐311++G and 6‐311++G(d,p). All the ab initio calculated points are fitted to the analytic Murrell‐Sorbie functions and then used to compute the spectroscopic parameters. The analytic potential energy function (APEF) for this b³Π_u state is reported. By comparison, the spectroscopic parameters reproduced by the APEF attained at 6‐311++G(2df,2pd) are found to be very close to the latest experimental findings. With the APEF obtained at the SAC‐CI/6‐311++G(2df,2pd) level of theory, a total of 62 vibrational states is found when J = 0. The complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants for these vibrational states are also reported. The reasonable dissociation limit for this state is deduced using the calculated results at present. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2007     

自由基HO2与C2H2反应势能面的理论研究

应用量子化学从头计算和密度泛函理论(DFT)对HO2+C2H2反应体系的反应机理进行了研究.在B3LYP/6-311G**和CCSD(T)/6-311G**水平上计算了HO2+ C2H2反应的二重态反应势能面.计算结果表明,主要反应方式为自由基HO2的H原子和C2H2分子中的C原子结合,经过一系列异构化,最后分解得到主要产物P1 (CH2O+ HCO).此反应是放热反应,化学反应热为-321.99 kJ·mol-1.次要产物为P2 (CO2 +CH3),也是放热反应.     

三重态SiCP2异构体的结构和稳定性的理论研究

采用DFT,QCISD和CCSD(T)等理论计算方法对三重态SiCP2异构体的结构和稳定性进行了理论研究.在B3LYP/6-311G(d)水平下,共计算得到由17个过渡态相连接的15个异构体.在CCSD(T)/6-311 +G(2df)//QCISD/6-311G(d)水平下,考虑重点振动能相对能量最低的三元环状异构体P-cCSiP 8(0.0 kJ/mol)及四元环状结构的cPCSiP 4具有相当大的动力学稳定性,在一定的实验室和星际条件下可能被检测到.另外,对它们的成键性质也进行了分析.     

Theoretical Study on the Potential Energy Curve and Vibration-rotation Spectra of BeF Radical

The potential energy curve (PEC) of BeF(X2Σ+) radical is investigated by using the complete active space self-consistent field (CASSCF) method followed by the highly accurate valence internally contracted multireference configuration interaction (MRCI) approach over the internuclear separation range from 0.0522 to 2.0472 nm. The PEC is fitted to the analytic Murrell-Sorbie function, which is employed to accurately determine the spectroscopic parameters. The present D0, De, Re, ωe, ωeχe, αe and Be are 6.14 eV, 6.22 eV, 0.1372 nm, 1236.12 cm-1, 9.11 cm-1, 0.0175 cm-1 and 1.4651 cm-1, respectively. These parameters have been compared with those of previous investigations reported in the literature. With PEC determined at the present level of theory, a total of 75 vibrational states have been predicted for the first time by numerically solving the radial Schrdinger equation of nuclear motion using the Numerov method. For each vibrational state, the complete vibrational levels, classical turning points, inertial rotation and centrifugal distortion constants are determined for the first time. Comparing with the available experiments and other theories, we find that the present spectroscopic parameter and molecular constant results are more accurate and complete than the previous theoretical investigations.     

La2Cn(n=-1,0,+1)团簇的理论研究

用B3LYP/Lan1 2dz方法优化了La2 Cn(n =- 1 ,0 ,+1 )分子的结构 ,计算了La2 C的电子亲和势和离化能 ,并对计算结果进行了讨论。     

Theoretical study of ScCO2+

The structure, binding energy, and vibrational frequencies have been determined for ScCO₂⁺. The inserted OSc⁺CO structure in the ¹A′ state is the most stable isomer and lies 43.2 kcal/mol below the ground-state Sc⁺+ CO₂ asymptote. The linear η¹-O Sc⁺(SINGLE BOND)OCO ³Δ state is bound by a charge-quadrupole interaction and has a binding energy of 13.9 kcal/mol. © 1997 John Wiley & Sons, Inc. Int J Quant Chem 63: 523–528, 1997