吡啶与HCl和CHCl3形成分子间红移氢键和蓝移氢键的理论研究

采用量子化学从头算的MP2方法, 分别在6-31G(d,p), 6-311+G(d,p)和AUG-cc-pVDZ基组下, 研究了复合物C5H5N…HCl(1), C5H5N…HCCl3(2)和C5H5N…HCCl3(3)的分子间氢键. 计算结果表明, 在复合物1中, HCl中Cl—H键伸长, 形成Cl—H…N红移氢键; 在复合物2中, HCCl3中C—H键伸长, 形成C—H…N 红移氢键; 在复合物3中, HCCl3中C—H键收缩, 形成C—H…π蓝移氢键. 自然键轨道(NBO)分析表明, 影响氢键红移和氢键蓝移主要有3个因素: n(Y)→σ*(X—H)超共轭作用、X—H键轨道再杂化和质子供体电子密度重排. 其中, 超共轭作用属于键伸长效应, 电子密度重排和轨道再杂化属于键收缩效应. 在复合物1和2中, 由于键伸长效应处于优势地位导致形成红移氢键; 在复合物3中, 由于键收缩效应处于优势地位导致形成蓝移氢键.     

吡咯与HX (X＝F,Cl, Br)分子间多种氢键的电子密度拓扑研究

采用密度泛函B3LYP/6-311＋＋G(d,p)方法, 对吡咯与HX (X＝F, Cl, Br)形成的经典氢键和π型氢键, 从其几何参数、电子密度的拓扑性质和电子积分等方面进行了研究. 在对π型氢键的讨论中我们将π电子与σ电子分离, 得到了π型氢键体系的π电子的密度等值线和拉普拉斯量等值线图以及各原子的π电子积分, 形象地说明了π型氢键的作用本质.     

HCCO自由基与水分子氢键复合物的从头算研究

在MP2/aug-cc-pvdz水平下,对二体氢键复合物H2O…HCCO(Ⅰ)和HCCO…H2O(Ⅱ)以及三体氢键复合物(H2O)2…HCCO(Ⅲ),H2O…H2O…HCCO(Ⅳ)和H2O…HCCO…H2O(Ⅴ)的几何和相互作用能进行了计算.轨道分析表明Ⅲ中HCCO中的H(1),C(2)通过2个氢键与2个水分子形成了环...     

肾上腺素与胞嘧啶相互作用的密度泛函理论研究

采用密度泛函理论（DFT）B3LYP方法和6-311＋G（d,p）基组对肾上腺素-胞嘧啶复合物进行结构优化和频率计算,得到15种稳定的复合物.研究发现,所有的复合物进行基组重叠误差（BSSE）校正后的相互作用能为-11.43^-48.96kJ/mol,符合氢键能量范围,相互作用能主要由氢键所贡献.结构和振动频率分析显示,氢键的形成使相应O（N）—H键的键长变长,对称伸缩振动频率减小,说明复合物中形成的氢键都是正常的红移型氢键.应用自然键轨道（NBO）理论和分子中的原子（AIM）理论对15种复合物的氢键性质和特征进行分析,发现氢键对于复合物的稳定性起着重要作用,当复合物形成2个或更多的氢键时,氢键的数目、类型及强度共同决定着复合物的稳定性,复合物基本符合三氢键〉二氢键〉单氢键的稳定顺序,三氢键复合物4是最稳定的,复合物3存在单氢键O—H…O,比部分二氢键复合物要稳定.     

邻二氮杂苯-水复合物的氢键结构与性质

用密度泛函理论B3LYP方法和MP2方法对邻二氮杂苯-水复合物基态的氢键结构与相互作用能进行了理论计算,结果表明复合物之间存在较强的氢键N…H－O.在复合物中,水的H－O对称伸缩振动频率明显红移.同时,使用含时密度泛函理论方法计算了邻二氮杂苯单体及复合物的低占据1(n,π*) 和1(π,π*) 态的垂直激发能,计算结果与实验值吻合较好.     

C6H5—H…X分子间氢键的理论计算

利用密度泛函理论B3LYP方法, 在6-311+G(3df,2p)水平上对C6H5—H…X型分子间氢键进行了几何构型优化、氢键相互作用能、电子密度分布等计算. 其中C6H6为质子供体, HCOH、H2O、NH3、CH2NH和HCN为质子受体. 从电荷布居分析、自然键轨道等角度详细地讨论了C6H5—H…X 体系中, 共轭π键、O和N的不同键型结构对氢键形成的影响以及孤电子对与C—H 反键轨道之间的相互作用(n→σ*)等.     

气相中开壳型(CH3)2S(O)…HOO红移氢键复合物的结构与性质

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…HOO和(CH3)2O…HOO开壳型氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOO自由基相比, 复合物中H10-O11键伸缩振动频率发生显著的红移, 红移值分别为424.21和374.22 cm-1. 在MP2/6-311++G**水平计算得到, 含基组重叠误差(BSSE)校正和零点振动能(ZPVE)校正的相互作用能分别为-24.68和-31.01 kJ·mol-1. 自然键轨道(NBO)理论分析表明, 在(CH3)2S…HOO复合物中, 引起H10-O11键变长的因素包括两种电荷转移: (1) LP(S1)1→σ*(H10-O11); (2) LP(S1)2→σ*(H10-O11), 其中LP(S1)2→σ*(H10-O11)占主要作用, 总的结果是使σ*(H10-O11)的自然布居数增加了37.27 me; 在(CH3)2O…HOO中也有相似的电荷转移的超共轭作用. AIM理论分析表明, S1…H10间和O1…H10间都存在键鞍点, ▽2ρ(r)分别为0.06196和0.03745, 说明这种相互作用介于共价键和离子键之间, 偏于静电作用.     

碱基腺嘌呤与多肽酰胺间氢键作用的最佳位点

正确理解核酸碱基和蛋白质多肽间的作用机制有助于人们利用这些生物分子有效地进行分子设计,进而制备具有特殊纳米结构和功能的生物分子材料.本文优化得到了碱基腺嘌呤与N-甲基乙酰胺、甘氨酸二肽、丙氨酸二肽形成的20个氢键复合物的结构并计算了结合能,探讨了腺嘌呤与多肽酰胺间氢键作用的最佳位点.研究发现:腺嘌呤可以使用两个不同位点(A1位点和A2位点)与N-甲基乙酰胺形成N―H…N型或者N―H…O=C型氢键复合物,腺嘌呤使用A1位点与N-甲基乙酰胺形成的N―H…N型氢键复合物更稳定;二肽分子可以使用主链上两个不同位点(丙氨酸的Ala7位点和Ala5位点或者甘氨酸的Gly7位点和Gly5位点)与腺嘌呤形成含有N―H…N和N―H…O=C两条氢键的复合物,二肽分子使用Ala7或Gly7位点与腺嘌呤形成的氢键复合物更稳定;腺嘌呤与多肽间的氢键作用强于其与N-甲基乙酰胺的作用.基于分子中的原子理论与自然键轨道计算结果分析了氢键作用的本质.     

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*)激发态的垂直激发能.     

(CH3)2S与HOCI分子间的卤键和氢键相互作用

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型.频率分析表明,与单体HOC1相比,在两种复合物中,10C1-11O和12H-11O键伸缩振动频率发生显著的红移.经MP2/6-311++G**水平计算的含基组重叠误差(BSSE)校正的气相中相互作用能分别为-11.69和-24.16 kJ·mol-1.自然键轨道理论(NBO)分析表明,在(CH3)2s…ClOH卤键复合物中,引起10Cl-11O键变长的因素包括两种电荷转移:(i)孤对电子LP(1S)1→σ*(10C1-11O);(ⅱ)孤对电子LP(1s)2→σ*(10Cl-11O),其中孤对电子LP(lS)2→σ*(10Cl-11O)转移占主要作用,总的结果是使σ(10Cl-11O)的自然布居数增加0.14035e,同时11O原子的再杂化使其与10Cl成键时s成分增加,即具有与电荷转移作用同样的"拉长效应";在(CH3)2s…HOCl氢键复合物中也存在类似的电荷转移,但是11O原子的再杂化不同于前者.自然键共振理论(NRT)进行键序分析表明,在卤键复合物和氢键复合物中,10Cl-11O和12H-11O键的键序都减小.通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.     

New calculations on the ion-molecule processes C2H2(+) + H2 --> C2H3(+) + H and C2H2(+) + H2 --> C2H4+

New high-level quantum chemical calculations have been undertaken to understand the rates and mechanisms of the reactive and associative channels for the reactants C2H2(+) + H2. The reactive channel, which produces C2H3(+) + H, has been shown to be slightly endothermic, confirming earlier calculations at a somewhat lower level and in agreement with some recent experimental work. The associative channel, leading to C2H4+, has been shown to proceed via a transition state with negative energy relative to the reactants, so that association is predicted to be efficient. This result is in conflict with an earlier theoretical study but in agreement with low-temperature experimental measurements.     

A pseudopotential study of the hydrogen bond in H2O·H2S,H2S·H2S and H2O·H2Se systems

Intermolecular potential energy curves for the hydrogen bonded systems H₂O·H₂S, H₂O·H₂Se and H₂S·H₂S were calculated with nonempirical pseudopotentials using optimized-in-molecules basis sets augmented by polarization functions. The H₂O·H₂O interaction energy curve has been also considered as a test case. The present results for H₂O·H₂S and H₂S·H₂S indicate much weaker intermolecular interactions than those found in previous ab initio calculations. The H₂O·H₂Se interaction was found to be quite similar to H₂O·H₂S.This work was partly supported by the Polish Academy of Sciences within the Project PAN-09, 7.1.1.1On leave from Quantum Chemistry Laboratory, Dept. of Chemistry, University of Warsaw, Pasteura 1, 02-093. Warsaw, Poland     

An Ab Initio Study of C2H2.H2, C2H2.N2, and C2H2.Ar Complexes

Ab initio calculations at the post Hartree–Fock level were performed on complexes of acetylene with hydrogen, nitrogen, and argon. Total energies, optimum geometries, and binding energies were calculated, using the 6-311G** and the 6-31+G(2df,2pd) basis sets. Calculations showed the complexes to be more stable than the separate entities, with the exception of the acetylene–hydrogen complex.     

Deuterium exchange in the systems of H2O+/H2O and H3O+/H2O

The reactions of H₂O⁺, H₃O⁺, D₂O⁺, and D₃O⁺ with neutral H₂O and D₂O were studied by tandem mass spectrometry. The H₂O⁺ and D₂O⁺ ion reactions exhibited multiple channels, including charge transfer, proton transfer (or hydrogen atom abstraction), and isotopic exchange. The H₃O⁺ and D₃O⁺ ion reactions exhibited only isotope exchange. The variation in the abundances of all ions involved in the reactions was measured over a neutral pressure range from 0 to 2 × 10⁻⁵ Torr. A reaction scheme was chosen, which consisted of a sequence of charge transfer, proton transfer, and isotopic exchange reactions. Exact solutions to two groups of simultaneous differential equations were determined; one group started with the reaction of ionized water, and the other group started with the reactions of protonated water. A nonlinear least-squares regression technique was used to determine the rate coefficients of the individual reactions in the schemes from the ion abundance data. Branching ratios and relative rate coefficients were also determined in this manner.A delta chi-squared analysis of the results of the model fitted to the experimental data indicated that the kinetic information about the primary isotopic exchange processes is statistically the most significant. The errors in the derived values of the kinetic information of subsequent channels increased rapidly. Data from previously published selected ion flow tube (SIFT) study were analyzed in the same manner. Rigorous statistical analysis showed that the statistical isotope scrambling model was unable to explain either the SIFT or the tandem mass spectrometry data. This study shows that statistical analysis can be utilized to assess the validity of possible models in explaining experimentally observed kinetic behaviors.     

Rates of reaction of atomic oxygen with C2H3F,C2H3Cl,C2H3Br, 1,1-C2H2F2, and 1,2-C2H2F2

Rate constants for the reactions of atomic oxygen (O³P) with C₂H₃F, C₂H₃Cl, C₂H₃Br, 1,1-C₂H₂F₂, and 1,2-C₂H₂F₂ have been measured at 307°K using a discharge-flow system coupled to a mass spectrometer. The rate constants for these reactions are (in units of 10¹¹ cm³ mole^?1 s^?1) 2.63 ± 0.38, 5.22 ± 0.24, 4.90 ± 0.34, 2.19 ± 0.18, and 2.70 ± 0.34, respectively. For some of these reactions, the product carbonyl halides were identified.     

Proton magnetic shielding in H2O and (H2O)2

The proton magnetic shielding constants in the water molecule and its linear perpendicular dimer are computed from SCF-MO-LCGO wave functions by using the uncoupled Hartree-Fock variation-perturbation procedure due to Karplus and Kolker. The convergence of the calculated shielding constants as well as their gauge dependence is studied. The final results for 17-term polynomial variation function indicate that the best choice for the gauge origin corresponds to the molecular electronic centroid.The calculated proton magnetic shielding constant in the water molecule is in remarkable agreement with experimental data and favourably compares with the best coupled Hartree-Fock results. It follows from the calculations for the water dimer that the H-bond NMR-shift amounts in this case —1.0 ppm and qualitatively agrees with the experimental data for the liquid water.     

Pseudopotential calculations on Si2H6 and Si2H4

Local pseudopotential calculations have been performed for the ground state of disilane as well as for the lowest singlet and triplet states of disilene and silylsilylene. Comparison with all-electron calculations shows good agreement for geometries and relative stabilities.