NH3BH3的双接触蓝移双氢键的理论研究

运用量子化学从头算方法研究了NH3BH3与含大的分子内超共轭的质子给体CHF3, H2CO, HCOOH, HCOCl和HNO形成的双接触弯曲双氢键B—H2…H—X(X＝C, N)的分子结构、电子密度拓扑性质与频率位移特征. 计算结果表明, 在所有体系中, 由于双氢键的形成, B—H键拉长且伸缩振动频率红移, 而X—H键长减小且伸缩频率蓝移. 弯曲的双接触构型导致分子间超共轭减小、X—H键的大的正重极化与正重杂化以及分子内超共轭减小三个因素导致了X—H键蓝移; B—H键红移的主要原因是B—H键的负的重极化与负的重杂化.     

双接触弯曲蓝移氢键CH_2…Y的理论研究

用量子化学从头算方法在MP2/6-311 G(d,p)与MP2/6-311 G(2df,2p)两个理论级别上研究了双接触弯曲氢键Y…H2CZ(Z=O,S,Se)和Y…H2CZ2(Z=F,Cl,Br)(Y=Cl-,Br-).计算表明,在这些复合物中都存在两个等价的Y…H—C蓝移氢键;相互作用能和蓝移都比较大,每个Y…H—C氢键的能量为15～27kJ/mol,CH键长变化为-0.1～-0.5pm,CH伸缩振动频率位移为30～80cm-1.自然键轨道分析表明,3个因素导致了这些蓝移氢键的形成:(1)存在较大的重杂化;(2)弯曲的氢键构型使得分子间直接超共轭相互作用比较小,而存在相当的分子间间接超共轭相互作用;(3)质子给体的分子内超共轭相互作用较大地减小.电子密度拓扑性质的研究表明,在这些氢键复合物中都存在3个分子间临界点:在接受体原子Y与每个H原子之间存在一个键临界点,也存在相应的键径和原子间界面;在YHCH四边形内部存在一个环临界点.因此这些分子间相互作用可以严格地看作氢键.     

三接触弯曲构型的蓝移氢键CH3…Y的理论研究

运用量子化学从头算方法, 在MP2/6-311＋＋G(d,p), MP2/6-311＋＋G(2df,2p), MP2/6-311＋＋G(3df,3pd)和QCISD/6-311＋＋G(d,p)水平上, 研究了CH₃F, CH₃Cl和CH₃Br作为质子给体与Cl^－, Br^－作为质子接受体形成的氢键CH₃…Y. 计算结果表明: 6种复合物中C—H键收缩, 伸缩振动频率增大, 形成蓝移氢键. 分子中原子(Atoms in Molecules, AIM)分析表明, 这些复合物的电子密度拓扑性质与普通氢键有着本质的不同, 在Y…H之间不存在键临界点, 而在Y与C之间存在键临界点, 因此这些相互作用严格地不能称为氢键. 自然键轨道(Natural bond orbital, NBO)分析表明, 在这些复合物中弯曲的CH…Y的特殊结构使得分子间超共轭n(Y)®σ*(C—H)减小到可以忽略; 质子接受体的电子密度没有转移到σ*(C—H)上, 而是转移到了σ*(C—X) (X＝F, Cl, Br)上; 存在一定程度的重杂化; 分子内超共轭相互作用减小使得σ*(C—H)的电子密度减少. 这些因素共同导致C—H伸缩振动频率的蓝移.     

三接触弯曲构型的蓝移氢键CH3…Y的理论研究

运用量子化学从头算方法, 在MP2/6-311＋＋G(d,p), MP2/6-311＋＋G(2df,2p), MP2/6-311＋＋G(3df,3pd)和QCISD/6-311＋＋G(d,p)水平上, 研究了CH₃F, CH₃Cl和CH₃Br作为质子给体与Cl^－, Br^－作为质子接受体形成的氢键CH₃…Y. 计算结果表明: 6种复合物中C—H键收缩, 伸缩振动频率增大, 形成蓝移氢键. 分子中原子(Atoms in Molecules, AIM)分析表明, 这些复合物的电子密度拓扑性质与普通氢键有着本质的不同, 在Y…H之间不存在键临界点, 而在Y与C之间存在键临界点, 因此这些相互作用严格地不能称为氢键. 自然键轨道(Natural bond orbital, NBO)分析表明, 在这些复合物中弯曲的CH…Y的特殊结构使得分子间超共轭n(Y)®σ*(C—H)减小到可以忽略; 质子接受体的电子密度没有转移到σ*(C—H)上, 而是转移到了σ*(C—X) (X＝F, Cl, Br)上; 存在一定程度的重杂化; 分子内超共轭相互作用减小使得σ*(C—H)的电子密度减少. 这些因素共同导致C—H伸缩振动频率的蓝移.     

双接触弯曲蓝移氢键CH2…Y的理论研究

用量子化学从头算方法在MP2／6-311＋＋G（d,p）与MP2／6—311＋＋G（2df,2p）两个理论级别上研究了双接触弯曲氢键Y…H2CZ（z=O,S,Se）和Y…H2CZ2（Z=F,Cl,Br）（Y=Clˉ,Br-）．计算表明,在这些复合物中都存在两个等价的Y…H—C蓝移氢键;相互作用能和蓝移都比较大,每个Y…H—C氢键的能量为15～27kJ／mol,CH键长变化为-0．1～-0．5pm,CH伸缩振动频率位移为30—80cm^-1．自然键轨道分析表明,3个因素导致了这些蓝移氢键的形成：（1）存在较大的重杂化;（2）弯曲的氢键构型使得分子间直接超共轭相互作用比较小,而存在相当的分子间间接超共轭相互作用;（3）质子给体的分子内超共轭相互作用较大地减小．电子密度拓扑性质的研究表明,在这些氢键复合物中都存在3个分子间临界点：在接受体原子Y与每个H原子之间存在一个键临界点,也存在相应的键径和原子间界面;在YHCH四边形内部存在一个环临界点．因此这些分子间相互作用可以严格地看作氢键．     

乙二醇的分子间氢键结构动力学的飞秒非线性红外光谱

利用稳态线性红外光谱和飞秒泵浦-探测红外光谱技术, 研究了在乙腈(MeCN)、丙酮(AC)、四氢呋喃(THF)和二甲基亚砜(DMSO)溶剂中乙二醇(EG)的结构和羟基(―OH)伸缩振动动力学. 结果表明, 乙二醇的―OH伸缩振动的频率位置、峰宽以及振动弛豫动力学都表现出强烈的溶剂依赖性. 乙二醇溶液中至少存在两种形式的分子间氢键, 一种是溶质-溶剂团簇的分子间氢键, 另一种是溶质-溶质团簇的分子间氢键. 量子化学计算预测的―OH伸缩振动频率的溶剂依赖性与我们的红外光谱实验观测结果一致. 进一步, 我们发现在乙腈中参与形成溶质-溶剂团簇氢键的乙二醇―OH伸缩振动具有最慢的弛豫动力学, 丙酮和四氢呋喃次之, 而最快的弛豫动力学过程发生在二甲基亚砜中. 在每一溶剂条件下, 乙二醇/乙二醇溶质团簇中―OH伸缩振动弛豫都更快一些. 本文结果有助于认识在溶质-溶质、溶质-溶剂分子团簇共存的体系中不同分子间氢键的结构动力学特性.     

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

采用密度泛函理论（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,比部分二氢键复合物要稳定.     

呋喃与HCl和CHCl3构成的分子间氢键的理论研究

量子化学从头算方法MP2研究了呋喃-HCl体系和呋喃-CHCl₃体系的分子间氢键的本质. 主要研究了这两个体系中新的氢键类型C(Cl)—H…O和C(Cl)—H…π的相互作用. 研究表明, 氯仿与呋喃分子之间的相互作用使氯仿中C—H键长缩短, 振动频率增大(蓝移), 而HCl与呋喃分子之间的相互作用使H—Cl键长增长, 振动频率减小(红移). 利用自然键轨道(NBO)分析表明电荷从质子受体转移到C—H反键轨道和Cl原子的孤对电子轨道上.     

用红外光谱技术研究氢键的键合方式

氢键的形成能够使参与形成氢键的原化学键力常数降低,吸收频率移向低波数方向,同时吸收强度增加。为了介绍怎样利用红外光谱技术研究氢键的键合方式,本文以联酰胺衍生物为例,主要运用变温红外光谱技术分析NH伸缩振动波数、强度以及形成氢键的键长随温度的变化来研究羰基(C O)与氨基(H—N)之间的氢键形式。结果表明,本文举例中联酰胺基团中的C O与H—N以分子间氢键形式存在。     

CH_3SH和CH_3SLi分子间的锂键和氢键共存型相互作用

在CH3SLi+CH3SH势能面上求得锂键和氢键共存型复合物的两种稳定构型.频率分析表明,与单体相比复合物中S(5)—Li(6)键伸缩振动频率发生红移,而C(8)—H(10)键伸缩振动频率发生蓝移.经B3LYP/6-311++G**,MP2/6-311++G**及MP2/AUG-CC-PVDZ水平计算的含基组重叠误差(BSSE)校正的复合物Ⅰ中相互作用能分别为-58.99,-57.87和-62.89kJ·mol-1.采用自然键轨道(NBO)理论,分析了复合物中单体轨道间的电荷转移,电子密度重排及其与相关键键长变化的本质等.采用分子中的原子(AIM)理论分析了复合物中氢键和锂键的电子密度拓扑性质.在极化连续模型(PCM)下,考察了溶剂化效应.结果表明,所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大,且溶剂对复合物中的锂键结构及其振动频率具有显著的影响,而对复合物中的氢键的振动频率影响不大.     

On the correlation between bond-length change and vibrational frequency shift in halogen-bonded complexes

The C-Hal (Hal = Cl, Br, or I) bond-length change and the corresponding vibrational frequency shift of the C-Hal stretch upon the C-Hal···Y (Y is the electron donor) halogen bond formation have been determined by using density functional theory computations. Plots of the C-Hal bond-length change versus the corresponding vibrational frequency shift of the C-Hal stretch all give straight lines. The coefficients of determination range from 0.94366 to 0.99219, showing that the correlation between the C-Hal bond-length change and the corresponding frequency shift is very good in the halogen-bonded complexes. The possible effects of vibrational coupling, computational method, and anharmonicity on the bond-length change-frequency shift correlation are discussed in detail.     

Quantitative relationships between bond lengths,stretching vibrational frequencies,bond force constants,and bond orders in the hydrogen-bonded complexes involving hydrogen halides

To uncover the correlation between the bond length change and the corresponding stretching frequency shift of the proton donor D–H upon hydrogen bond formation, a series of hydrogen-bonded complexes involving HF and HCl which exhibit the characteristics of red-shifted hydrogen bond were investigated at the MP2/aug-cc-pVTZ, M062X/aug-cc-pVTZ, and B3LYP/aug-cc-pVTZ(GD3) levels of theory with CP optimizations. A statistical analysis of these complexes leads to the quantitative illustrations of the relations between bond length and stretching vibrational frequency, between bond length and bond force constant, between stretching vibrational frequency and bond force constant, between bond length and bond order for hydrohalides in a mathematical way, which would provide valuable insights into the explanation of the geometrical and spectroscopic behaviors during hydrogen bond formation.     

Electronic basis of improper hydrogen bonding: a subtle balance of hyperconjugation and rehybridization

The X[bond]H bond length in X[bond]H...Y hydrogen bonded complexes is controlled by a balance of two main factors acting in opposite directions. "X[bond]H bond lengthening" due to n(Y)-->sigma(H[bond]X) hyperconjugative interaction is balanced by "X[bond]H bond shortening" due to increase in the s-character and polarization of the X[bond]H. When hyperconjugation dominates, X[bond]H bond elongation is reflected in a concomitant red shift of the corresponding IR stretching frequency. When the hyperconjugative interaction is weak and the X-hybrid orbital in the X[bond]H is able to undergo a sufficient change in hybridization and polarization, rehybridization dominates leading to a shortening of the X[bond]H and a blue shift in the X[bond]H stretching frequency.     

Theoretical analysis based on X-H bonding strength and electronic properties in red- and blue-shifting hydrogen-bonded X-H···π complexes

A theoretical study based on the X-H bond strength of the proton donor fragment and its concomitant classical red-shifting or improper blue-shifting of the pure stretching frequency, in weakly hydrogen-bonded X-H···π complexes, is presented. In this sense, the dissociation energy differences, defined as, ΔD(e) = D(e)(X-H)[complex] - D(e)(X-H) [isolated], showed to be linearly connected with the change in stretching frequencies, Δν = ν(X-H)[complex] - ν(X-H)[isolated], of red- and blue-shifting H-bonds. This relationship allows us to define a threshold for the type of the stretching shift of the X-H bond: ΔD(e)(X-H) > 50.3 kcal mol(-1) leads to blue-shifting whereas ΔD(e)(X-H) < 50.3 kcal mol(-1) leads to red-shifting behavior. Complementarily, natural bond orbital analysis along the X-H stretching coordinate and electric dipole polarizability was performed to investigate the factors involved in red- or blue-shifting hydrogen-bonded complexes. It has been found that a high tendency to deplete the electronic population on the H atom upon X-H stretching is exhibited in blue-shifting H-bonded complexes. On the other hand, these types of complexes present a compact electronic redistribution in agreement with polarizability values. This study has been carried out taking as models the following systems: chloroform-benzene (Cl(3)C-H···C(6)H(6)), fluoroform-benzene (F(3)C-H···C(6)H(6)), chloroform-fluorobenzene, as blue-shifting hydrogen-bonded complexes and cyanide acid-benzene (NC-H···C(6)H(6)), bromide and chloride acids-benzene ((Br)Cl-H···C(6)H(6)) and acetylene-benzene (C(2)H(2)···C(6)H(6)) as red-shifting complexes.     

吡啶与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中, 由于键收缩效应处于优势地位导致形成蓝移氢键.     

Strongly blue-shifted C-H stretches: interaction of formaldehyde with hydrogen fluoride clusters

The equilibrium structures, binding energies, and vibrational spectra of the cyclic, hydrogen-bonded complexes formed between formaldehyde, H(2)CO, and hydrogen fluoride clusters, (HF)(1< or =n < or =4), are investigated by means of large-scale second-order M?ller-Plesset calculations with extended basis sets. All studied complexes exhibit marked blue shifts of the C-H stretching frequencies, exceeding 100 cm(-1) for n = 2-4. It is shown that these blue shifts are, however, only to a minor part caused by blue-shifting hydrogen bonding via C-H...F contacts. The major part arises due to the structural relaxation of the H(2)CO molecule under the formation of a strong C=O...H-F hydrogen bond which strengthens as n increases. The close correlation between the different structural parameters in the studied series of complexes is demonstrated, and the consequences for the frequency shifts in the complexes are pointed out, corroborating thus the suggestion of the primary role of the C=O...H-F hydrogen bonding for the C-H stretching frequency shifts. This particular behavior, that the appearance of an increasingly stronger blue shift of the C-H stretching frequencies is mainly induced by the formation of a progressively stronger C=O...H-F hydrogen bond in the series of H(2)CO...(HF)(1< or =n < or =4), complexes and only to a lesser degree by the formation of the so-called blue-shifting C-H...F hydrogen bond, is rationalized with the aid of selected sections of the intramolecular H(2)CO potential energy surface and by performing a variety of structural optimizations of the H(2)CO molecule embedded in external, differently oriented dipole electric fields, and also by invoking a simple analytical force-field model.     

On the correlation between bond-length change and vibrational frequency shift in hydrogen-bonded complexes: a computational study of Y...HCl dimers (Y = N2, CO, BF)

The H-Cl bond-length change and the harmonic vibrational frequency shift of the H-Cl stretch on formation of the linear isoelectronic Y...H-Cl complexes (Y = N(2), CO, BF) have been determined by ab initio computations at different levels of theory. These shifts are in agreement with predictions from a model based on perturbation theory and involving the first and second derivatives of the interaction energy with respect to displacement of the H-Cl bond length from its equilibrium value in the isolated monomer. At the highest level of theory, blue shifts were obtained for BF...HCl and CO...HCl, while red shifts were obtained for FB...HCl, OC...HCl, and N(2)...HCl. These vibrational characteristics are rationalized by considering the balance between the interaction energy derivatives obtained from the perturbative model. The widely believed correlation between the bond-length change and the sign of the frequency shift obtained on complexation is discussed and found to be unreliable.     

Ab initio study of the influence of trimer formation on one- and two-bond spin-spin coupling constants across an X-H-Y hydrogen bond: AH:XH:YH3 complexes for A, X = 19F, 35Cl and Y = 15N, 31P

Ab initio equation-of-motion coupled-cluster singles and doubles (EOM-CCSD) calculations have been carried out to investigate the effect of a third polar near-neighbor on one-bond ((1)J(X)(-)(H) and (1h)J(H)(-)(Y)) and two-bond ((2h)J(X)(-)(Y)) spin-spin coupling constants in AH:XH:YH(3) complexes, where A and X are (19)F and (35)Cl and Y is either (15)N or (31)P. The changes in both one- and two-bond spin-spin coupling constants upon trimer formation indicate that the presence of a third molecule promotes proton transfer across the X-H-Y hydrogen bond. The proton-shared character of the X-H-Y hydrogen bond increases in the order XH:YH(3) < ClH:XH:YH(3) < FH:XH:YH(3). This order is also the order of decreasing shielding of the hydrogen-bonded proton and decreasing X-Y distance, and is consistent with the greater hydrogen-bonding ability of HF compared to HCl as the third molecule. For all complexes, the reduced X-H and X-Y spin-spin coupling constants ((1)K(X)(-)(H) and (2h)K(X)(-)(Y)) are positive, consistent with previous studies of complexes in which X and Y are second-period elements in hydrogen-bonded dimers. (1h)K(H)(-)(Y) is, as expected, negative in these complexes which have traditional hydrogen bonds, except for ClH:FH:NH(3) and FH:FH:NH(3). In these two complexes, the F-H-N hydrogen bond has sufficient proton-shared character to induce a change of sign in (1h)K(H)(-)(Y). The effects of trimer formation on spin-spin coupling constants are markedly greater in complexes in which NH(3) rather than PH(3) is the proton acceptor.     

Halogen-hydride interaction between Z-X (Z = CN, NC; X = F, Cl, Br) and H-Mg-Y (Y = H, F, Cl, Br, CH3)

Halogen-hydride interactions between Z-X (Z = CN, NC and X = F, Cl, Br) as halogen donor and H-Mg-Y (Y = H, F, Cl, Br, CH(3)) as electron donor have been investigated through the use of Becke three-parameter hybrid exchange with Lee-Yang-Parr correlation (B3LYP), second-order M?ller-Plesset perturbation theory (MP2), and coupled-cluster single and double excitation (with triple excitations) [CCSD(T)] approaches. Geometry changes during the halogen-hydride interaction are accompanied by a mutual polarization of both partners with some charge transfer occurring from the electron donor subunit. Interaction energies computed at MP2 level vary from -1.23 to -2.99 kJ/mol for Z-F···H-Mg-Y complexes, indicating that the fluorine interactions are relatively very weak but not negligible. Instead, for chlorine- and bromine-containing complexes the interaction energies span from -5.78 to a maximum of -26.42 kJ/mol, which intimate that the interactions are comparable to conventional hydrogen bonding. Moreover, the calculated interaction energy was found to increase in magnitude with increasing positive electrostatic potential on the extension of Z-X bond. Analysis of geometric, vibrational frequency shift and the interaction energies indicates that, depending on the halogen, CN-X···H interactions are about 1.3-2.0 times stronger than NC-X···H interactions in which the halogen bonds to carbon. We also identified a clear dependence of the halogen-hydride bond strength on the electron-donating or -withdrawing effect of the substituent in the H-Mg-Y subunits. Furthermore, the electronic and structural properties of the resulting complexes have been unveiled by means of the atoms in molecules (AIM) and natural bond orbital (NBO) analyses. Finally, several correlative relationships between interaction energies and various properties such as binding distance, frequency shift, molecular electrostatic potential, and intermolecular density at bond critical point have been checked for all studied systems.     

Halogen as halogen-bonding donor and hydrogen-bonding acceptor simultaneously in ring-shaped H3N·X(Y)·HF (X = Cl, Br and Y = F, Cl, Br) complexes

A series of ring-shaped molecular complexes formed by H(3)N, HF and XY (X = Cl, Br and Y = F, Cl, Br) have been investigated at the MP2/aug-cc-pVTZ level of theory. Their optimized geometry, stretching mode, and interaction energy have been obtained. We found that each complex possesses two red-shifted hydrogen bonds and one red-shifted halogen bond, and the two hydrogen bonds exhibit strong cooperative effects on the halogen bond. The cooperativity among the NH(3)···FH, FH···XY and H(3)N···XY interactions leads to the formations of these complexes. The AIM analysis has been performed at the CCSD(T)/aug-cc-pVQZ level of theory to examine the topological characteristics at the bond critical point and at the ring critical point, confirming the coexistence of the two hydrogen bonds and one halogen bond for each complex. The NBO analysis carried out at the B3LYP/aug-cc-pVTZ level of theory demonstrates the effects of hyperconjugation, hybridization, and polarization coming into play during the hydrogen and halogen bonding formations processes, based on which a clockwise loop of charge transfer was discovered. The molecular electrostatic potential has been employed to explore the formation mechanisms of these molecular complexes.