甲烷水合物分子间势能的量子化学研究

用Hartree-Fock SCF和密度泛函(BLYP,B3LYP,MPW1PW91)方法对以结构-Ⅰ为单元的甲烷水合物进行了分子间势能的理论研究.该结构单元为正十二面体,其中包括20个水分子,甲烷分子在其中心.采用从头算HF/6-31G(d,p)对甲烷分子进行几何优化,采用ST2模型对水分子作几何优化.水-水间氢键势能Ehb(l)和水-甲烷间范德华势能Evdw(l)作为边长l的函数进行计算,计算时固定水和甲烷分子的几何形状.所有计算中均使用6-31G(d,p)基组.基组重叠误差(BSSE)经校正其上限和下限为水-水氢键能加以确定.由B3LYP经基组重叠误差(BSSE)校正得到的O—O距离为RO—O=0.280 nm,C—O距离RC—O=0.392 nm,比其他方法更接近实验值的0.282和0.395 nm.结果表明,在天然气水合物结构-Ⅰ中水-水分子对的氢键能(30~36 kJ/mol)大于水的二聚体(H2O)2氢键能(-22.6±2.9)kJ/mol,亦大于六角形冰的(-21.7±0.5)kJ/mol,十二面体结构为一稳定单元.以上分子间相互作用势能的结果为得出Lennard-Jones和Kihara势能参数提供了坚实的基础,此参数对分子动力学模拟天然气水合物是非常有用的.     

脂肪族氨基酸二肽与水团簇的理论研究

用ABEEMσπ/MM模型和MP2/6-31+G(d)//B3LYP/6-31G(d)方法研究水合效应对脂肪族氨基酸二肽的影响.从结构和能量两方面说明Leu残基在蛋白质中起成旋作用,Val和Ⅱe残基在蛋白质中起解旋作用.同时得出:水分子严重影响了二肽分子的骨架二面角;对于结合相同数目水分子的团簇倾向于形成含有环状氢键的结构,并且含有环状氢键团簇的结合能大于含有链状氢键团簇的结合能.     

用GGA密度泛函及其长程、色散校正方法计算各类氢键的结合能

应用广义梯度近似(GGA) (PW91和PBE)、含动能密度的广义梯度近似(meta-GGA) (M06-L)、杂化泛函(hyper-GGA)(M06-2X、X3LYP和B3LYP)及其长程校正泛函LC-DFT(CAM-B3LYP、LC-ωPBE和ωB97X)和色散校正密度泛函(DFT-D)(ωB97X-D和B97-D),用多种基函数对15种不同强度的传统氢键和非传统氢键体系的结合能进行了系统的计算与分析.并与高精度的CCSD(T)/aug-cc-pVQZ结果比较发现:在上述各类泛函中,对于氢键结合能的计算M06-2X和ωB97X-D泛函较为精确与可靠,且没有必要使用过大的基函数,6-311++G(2d,2p)或aug-cc-pVDZ水平的基组就已足够,各类泛函所计算结合能的基组重叠误差(BSSE)均较小,除ωB97X和ωB97X-D外,其它9种泛函不经BSSE校正也能得到同样甚至更准确的结果.     

乙醇-水分子团簇C2H5OH(H2O)n (n=1-9)稳定结构的量子化学研究

采用密度泛函理论B3LYP方法, 在B3LYP/6-311++G(2d,2p)//B3LYP/6-311++G(d,p)基组水平上对乙醇-水分子团簇(C₂H₅OH(H₂O)_n (n=1-9))的各种性质进行研究, 如: 优化的几何构型、结构参数、氢键、结合能、平均氢键强度、自然键轨道(NBO)电荷分布、团簇的生长规律等. 结果表明, 从二维(2-D)环状结构到三维(3-D)笼状结构的过渡出现在n=5的乙醇-水分子团簇中. 此外, 利用团簇结合能的二阶差分、形成能、能隙等性质, 发现在n=6时乙醇-水分子团簇的最低能量结构稳定性较好, 可能为幻数结构. 最后, 为了进一步探讨氢键本质, 将C₂H₅OH(H₂O)_n (n=2-9)最低能量结构的各种性质与纯水分子团簇(H₂O)_n (n=3-10)比较, 结果表明前者与后者中的水分子之间氢键相似.     

s-四嗪-水簇复合物的理论研究

用量子化学B3LYP方法和6-31＋＋G**基函数研究了s-四嗪-水簇复合物基态分子间相互作用, 并进行了构型优化和频率计算, 分别得到无虚频稳定的s-四嗪-(水)₂复合物、s-四嗪-(水)₃复合物和s-四嗪-(水)₄复合物6个、9个和12个. 复合物存在较强的氢键作用, 复合物结构中形成一个N…H—O氢键并终止于O…H—C氢键的氢键水链构型最稳定. 经基组重叠误差和零点振动能校正后, 最稳定的1∶2, 1∶3和1∶4(摩尔比)复合物的结合能分别是41.35, 70.9和 94.61 kJ/mol. 振动分析显示氢键的形成使复合物中水分子H—O键对称伸缩振动频率减小(红移). 研究表明N…H键越短, N…H—O键角越接近直线, 稳定化能越大, 氢键作用越强. 同时, 用含时密度泛函理论方法在TD-B3LYP/6-31＋＋G**水平计算了s-四嗪单体及其氢键复合物的第一¹(n, p*)激发态的垂直激发能.     

二水合甘氨酸两性离子复合体的结构和性能的理论研究

采用理论计算方法B3LYP, 在6-31++G**基组水平研究使甘氨酸质子化所需的最少水分子数目, 然后讨论水合两性离子复合体的结构和性能, 进而计算了二水合甘氨酸中性分子复合体(2W-GN)到二水合甘氨酸两性离子复合体(2W-GZ)的过渡态, 得到如下结论: (1)两个水分子可以使甘氨酸质子化, 能够形成稳定的二水合两性离子复合体. (2)甘氨酸与水分子之间通过氢键相互作用, 结合能较大, 复合体稳定; 在二水合甘氨酸复合体中, 水合甘氨酸中性分子比水合甘氨酸两性离子稳定. (3)由2W-GN到2W-GZ过程的反应活化能和氢键键能相近.     

水合肼在高岭石层间插层行为的量子化学研究

水合肼以其碱性及吸附性受到越来越多的关注,同时它在粘土中的污染问题也越来越受到重视.本工作构建了高岭石团簇模型为Al6Si6O42H42并在B3LYP/6-31G(d,p),MP2/6-31G(d,p)//B3LYP/6-31G(d,p)和MP2/6-31++G(d,p)//B3LYP/6-31G(d,p)水平下对一水合肼以及二水合肼在高岭石层间的插层性质(如:优化构型、结构参数、结合能、电荷分布、振动光谱、静电势等)进行探究.计算表明,当一水合肼进入层间后,水分子和肼分子之间的相互作用发生了改变.即水与肼分子分别以氢键的形式插层于高岭石层间,且肼与高岭石之间的相互作用要强于肼与水之间的相互作用,同时插层位点多位于高岭石四面体层和八面体层的重叠区域内,这些都是水合肼易进入高岭石层间而难以脱去的重要因素.当二水合肼进入层间后,随着层间距的不断扩大,肼分子与高岭石铝氧层之间的相互作用仍强于肼分子与水分子间的作用.但当层间距超过1.05 nm时,水分子与肼分子之间的作用则强于肼分子与高岭石的作用,这也印证了若要将肼脱附,需将层间距增大以减弱肼分子与高岭石的作用,再用溶剂将其脱附的可行性.     

6-烷基鸟嘌呤与DNA碱基间氢键作用的理论研究

用密度泛函B3LYP方法在6-311+G**基组水平上对鸟嘌呤及顺(cis-)、反式(anti-)-6-烷基鸟嘌呤(O6-AlkylG)与DNA碱基(胸腺嘧啶T、胞嘧啶C、腺嘌呤A、鸟嘌呤G)的氢键二聚体结构进行了优化.在MP2/cc-pVXZ(X=D,T)//B3LYP/6-311+G**水平上,采用完全基组外推方法校正了氢键二聚体的相互作用能,并用完全均衡校正法(CP)校正了基组重叠误差(BSSE).在B3LYP/6-311+G**水平上计算了各氢键碱基对的全电子波函数,并用分子中的原子理论(AIM)分析了碱基间的弱相互作用.计算结果显示,鸟嘌呤6-O烷基化改变了碱基间的氢键作用模式,使碱基对发生了明显的螺旋桨式扭转和不同程度的位移,碱基间的电子密度分布和氢键作用能明显减小.O6-AlkylG对DNA碱基间的氢键作用是去稳定化的,去稳定化影响的顺序为GCGGGA≈GT.计算结果与文献给出的实验结论基本一致.     

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

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

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

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

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

Global search for minimum energy (H2O)n clusters, n = 3-5

The Gaussian-3 (G3) model chemistry method has been used to calculate the relative deltaG(o) values for all possible conformers of neutral clusters of water, (H2O)n, where n = 3-5. A complete 12-fold conformational search around each hydrogen bond produced 144, 1728, and 20,736 initial starting structures of the water trimer, tetramer, and pentamer. These structures were optimized with PM3, followed by HF/6-31G* optimization, and then with the G3 model chemistry. Only two trimers are present on the G3 potential energy hypersurface. We identified 5 tetramers and 10 pentamers on the potential energy and free-energy hypersurfaces at 298 K. None of these 17 structures were linear; all linear starting models folded into cyclic or three-dimensional structures. The cyclic pentamer is the most stable isomer at 298 K. On the basis of this and previous studies, we expect the cyclic tetramers and pentamers to be the most significant cyclic water clusters in the atmosphere.     

Efficiency of numerical basis sets for predicting the binding energies of hydrogen bonded complexes: evidence of small basis set superposition error compared to Gaussian basis sets

Binding energies of selected hydrogen bonded complexes have been calculated within the framework of density functional theory (DFT) method to discuss the efficiency of numerical basis sets implemented in the DFT code DMol3 in comparison with Gaussian basis sets. The corrections of basis set superposition error (BSSE) are evaluated by means of counterpoise method. Two kinds of different numerical basis sets in size are examined; the size of the one is comparable to Gaussian double zeta plus polarization function basis set (DNP), and that of the other is comparable to triple zeta plus double polarization functions basis set (TNDP). We have confirmed that the magnitudes of BSSE in these numerical basis sets are comparative to or smaller than those in Gaussian basis sets whose sizes are much larger than the corresponding numerical basis sets; the BSSE corrections in DNP are less than those in the Gaussian 6-311+G(3df,2pd) basis set, and those in TNDP are comparable to those in the substantially large scale Gaussian basis set aug-cc-pVTZ. The differences in counterpoise corrected binding energies between calculated using DNP and calculated using aug-cc-pVTZ are less than 9 kJ/mol for all of the complexes studied in the present work. The present results have shown that the cost effectiveness in the numerical basis sets in DMol3 is superior to that in Gaussian basis sets in terms of accuracy per computational cost.     

A theoretical study of nonlinearity in heterocyclic hydrogen-bonded complexes

Ab initio calculations are performed at the MP2/6-311++G(d,p) and DFT/B3LYP/6-311++G(d,p) theoretical levels to obtain geometries, H-bond energies and harmonic infrared vibrational properties for the Cs symmetry structures of heterocyclic hydrogen-bonded complexes, CnHmY-HX. The H-bond lengths in DFT/B3LYP calculation level are in better agreement with the experimental values than the MP2 results. The geometry optimization are interpreted in terms of hydrogen bond nonlinearity represented by theta; and phi angles, once the hydrogen bond is formed among n-electrons pairs of the heteroatom in heterocyclic and the hydrogen atom in HX. The hydrogen bond energy after of the zero-point vibrational energy (ZPE) and basis set superposition error (BSSE) corrections are overestimated at DFT/B3LYP, whereas the MP2 BSSE corrections are very large than corresponding DFT/B3LYP. For example, the BSSE corrections for the C2H4S-HNC complex are 7.60 and 0.09 kJ mol(-1) in MP2 and DFT/B3LYP calculations levels, respectively. The new vibrational modes in infrared harmonic spectrum arising from complexation show several interesting features, especially the intermolecular stretching mode.     

2-羟基吡啶与水氢键作用的理论研究

本文采用量子化学的Hatree-Fock方法和密度泛函理论(DFT)的B3LYP方法,在6-31G(d)水平上,研究了2-羟基吡啶分子(Hy)及其酮式互变异构体2(1H)-吡啶酮(Py)与水的相互作用。考察它们之间在形成Hy…H2O,Py…H2O,Hy…Hy,Py…Py和Hy…Py等复合物前后的能量变化和分子结构参数变化特点。计算结果表明,在这些复合物中都形成了较强的氢键作用,在水合物中,Py与水形成复合物时能量降低较多,与实验结果一致。经过零点振动能(ZPVE)和基组叠加误差(BSSE)校正后的复合物离解能分别为38.3,40.8,73.0,82.7和71.1 kJ/mol(B3LYP/6-31G(d)),水合物的离解能远小于二聚体复合物,而酮式结构的二聚体的离解能最大。     

Theoretical study on the hydrogen bonding of five-membered heteroaromatics with water

The hydrogen-bonding ability of five-membered heteroaromatic molecules containing one chalcogen and two heteroatoms with nitrogen in addition to chalcogen, respectively, have been analyzed using density functional and molecular orbital methods through adduct formation with water. The stabilization energies for all the adducts are established at B3LYP/6-31+G* and MP2/6-31+G* levels after correcting for the basis set superposition error by using the counterpoise method and also corrected for zero-point vibrational energies. A natural bond orbital analysis at B3LYP/6-31+G* level and natural energy decomposition analysis at HF/6-31+G* using MP2/6-31+G* geometries have been carried out to understand the nature of hydrogen-bonding interaction in monohydrated heterocyclic adducts. Nucleus-independent chemical shift have been evaluated to understand the correlation between hydrogen bond formation and aromaticity.     

A density functional theory study of the structure and vibrational spectra of beta-carotene, capsanthin, and capsorubin

A theoretical study of the structure and the vibrational spectra of the beta-carotene molecule and its derivatives capsanthin and capsorubin is carried out. We first investigate systematically the theoretical method which provides the best results for beta-carotene by performing ab initio calculations at the HF/6-31G(d), SVWN/6-31G(d), PBE0/6-31G(d), BLYP/6-31G(d), B3LYP/6-31G(d), B3LYP/6-31G(d,p), B3LYP/6-311G(d), and B3LYP/6-311G(d,p) levels and by using previous theoretical results available in the literature obtained at the AM1 and BPW91/6-31G(d) levels. The influence of both the level of calculation and the size of the basis set used in the geometry optimization and in the determination of the IR and Raman spectra of this molecule is thus analyzed. It is confirmed that the hybrid functional B3LYP with the basis 6-31G(d) is the method that gives the best results as a whole. By use of this level of calculation, we next optimize the molecular geometries of related molecules of capsanthin and capsorubin, which to the best of our knowledge have only been studied at the semiempirical AM1 level. In addition we calculate the IR and Raman spectra of these molecules at the B3LYP/6-31G(d) level of theory. The results obtained for capsanthin show on the one hand that the double bond of the beta-ionone ring is outside the polyene chain plane, due to the repulsion between the hydrogen atoms of the ring methyl groups and the hydrogen atoms of the polyene chain, and on the other hand that the carbonyl double bond in the other headgroup is very close to planarity with the polyene chain, since in this case such a repulsion does not exist. For the molecule of capsorubin the two carbonyl groups also take the same coplanar orientation relative to the polyene chain. The IR and Raman spectra theoretically computed for these two molecules are finally compared with their experimental spectra and the vibrational normal modes of the main signals are interpreted.     

Tert-butyl N-(2-bromocyclohex-2-enyl)-N-(2-furylmethyl)carbamate的结构和振动光谱

The molecular structure, conformafional stability, and vibrational frequencies of ten-butyl N-（2- bromocyclohex-2-enyl）-N-（2-furylmethyl）carbamate （TBBFC） were investigated by utilizing the Hartree-Fock （HF） and density functional theory （DFT） ab initio calculations with 6-31G ^＊ and 6-31G^＊ ＊ basis sets. The optimized bond length and angle values obtained by HF method showed the best agreement with the experimental values. Comparison of the observed and calculated fundamental vibrational frequencies indicated that B3LYP was superior to the scaled HF approach for molecular problems. Optimal uniform scaling factors calculated for the title compound are 0.899/0.904, 0.958/0.961, and 0.988/0.989 for HF, B3LYP, and BLYP （6-31G ^＊/6-31G ^＊ ＊）, respectively.     

Evaluation of the individual hydrogen bonding energies in N-methylacetamide chains

The individual hydrogen bonding energies in N-methylacetamide chains were evaluated at the MP2/6-31+G** level including BSSE correction and at the B3LYP/6-311++G(3df,2pd) level including BSSE and van der Waals correction. The calculation results indicate that compared with MP2 results, B3LYP calculations without van der Waals correction underestimate the individual hydrogen bonding energies about 5.4 kJ mol^?1 for both the terminal and central hydrogen bonds, whereas B3LYP calculations with van der Waals correction produce almost the same individual hydrogen bonding energies as MP2 does for those terminal hydrogen bonds, but still underestimate the individual hydrogen bonding energies about 2.5 kJ mol^?1 for the hydrogen bonds near the center. Our calculation results show that the individual hydrogen bonding energy becomes more negative (more attractive) as the chain becomes longer and that the hydrogen bonds close to the interior of the chain are stronger than those near the ends. The weakest individual hydrogen bonding energy is about ?29.0 kJ mol^?1 found in the dimer, whereas with the growth of the N-methylacetamide chain the individual hydrogen bonding energy was estimated to be as large as ?62.5 kJ mol^?1 found in the N-methylacetamide decamer, showing that there is a significant hydrogen bond cooperative effect in N-methylacetamide chains. The natural bond orbital analysis indicates that a stronger hydrogen bond corresponds to a larger positive charge for the H atom and a larger negative charge for the O atom in the N-H?O=C bond, corresponds to a stronger second-order stabilization energy between the oxygen lone pair and the N-H antibonding orbital, and corresponds to more charge transfer between the hydrogen bonded donor and acceptor molecules.