苯酚-酰胺系列氢键复合物从头计算研究

运用G94W量子化学程序包,在HF/6-31G基组水平上对酰胺(DMF,DMA,HCONH2,HCONHCH3andCH3CONH2)与苯酚形成的系列氢键复合物(看作超分子)进行从头计算研究。根据计算结果探讨复合物的稳定性、施体和受体间的电荷转移及几何参数变化等规律。结果表明苯酚与上述一系列酰胺都可形成稳定的氢键复合物,其稳定性次序为CH3CONH2～HCONHCH3＞HCONH2＞DMA＞DMF。结果还表明形成氢键复合物的过程包含着电荷转移,电荷由供体酰胺转移到受体苯酚中,酰胺中C=O键长和苯酚中的O-H键长都明显有规律性地变长。计算结果与实验规律相符。     

一系列氢键诱导液晶复合物的合成及表征

以２－氰基－３－（４－（２－氯－３－甲基－１－丁酰氧基）－苯基丙烯酸（Ａ）为质子给体，Ｎ－（４－吡啶基亚甲基）－４－烷氧基苯胺（ｎＳＳＺ）为质子受体，合成了一系列新的氢键复和物，经红外光谱证实了分子间氢键的存在，通过ＤＳＣ，偏光方法及Ｘ射线衍射方法对其液晶行为进行研究，结果表明复合物呈现近晶相行为。     

取代基对1-甲基尿嘧啶与N-甲基乙酰胺氢键复合物中氢键强度的影响

使用密度泛函理论B3LYP方法和二阶微扰理论MP2方法对由1-甲基尿嘧啶与N-甲基乙酰胺所形成的氢键复合物中的氢键强度进行了理论研究, 探讨了不同取代基取代氢键受体分子1-甲基尿嘧啶中的氢原子对氢键强度的影响和氢键的协同性. 研究表明: 供电子取代基使N－H…O＝C氢键键长r(H…O)缩短, 氢键强度增强; 吸电子取代基使N－H…O＝C氢键键长r(H…O)伸长, 氢键强度减弱. 自然键轨道(NBO)分析表明: 供电子基团使参与形成氢键的氢原子的正电荷增加, 使氧原子的负电荷增加, 使质子供体和受体分子间的电荷转移量增多; 吸电子基团则相反. 供电子基团使N－H…O＝C氢键中氧原子的孤对电子轨道n(O)对N－H的反键轨道σ^*(N－H)的二阶相互作用稳定化能增强, 吸电子基团使这种二阶相互作用稳定化能减弱. 取代基对与其相近的N－H…O＝C氢键影响更大.     

氢键复合物自组装性质的研究

超分子复合物体系因近年来发展较快及应用范围广而颇具吸引力。作为形成超分子复合物的主要手段,氢键自组装对于高分子聚合材料以及生命科学等领域有着重要的意义。本文根据氢键的多重性及不同氢键的缔合方式对氢键复合物进行分类,并对氢键自组装复合物的性质和研究状况作了综述,同时介绍了本研究组在此方面的理论研究工作。     

Fapy-G对碱基氢键复合物影响的理论研究

应用量子化学方法对2,6-二氨基-4-羟基-5-甲酰胺嘧啶(Fapy-G)与正常碱基作用形成的20种氧化碱基对的多种性质进行了理论分析,碱基G的C8位被氧化后N7和N9位分别增加一个H原子,使其由氢键受体变成氢键供体,N7,N9及O6原子所带的电荷变负,同时O6作为氢键受体的能力增强.与碱基单体相比,碱基对中形成氢键的受体原子所带的电荷平均减小0.05 e;供体H原子所带的电荷平均增大0.02 e.Fapy-G分子中六元环上受体N原子参与形成氢键时,环的呼吸振动模式和N与对位C的振动模式的振动频率蓝移;与氢键相关的振动频率红移.所有氧化碱基对中,NH…N比NH…O氢键作用更强,且在NH…N氢键中,在六元环上的供体N原子形成的氢键比在氨基或开环上的供体N原子形成的氢键强.Fapy-G与碱基A作用结合能区域顺序为1>2>4>3,与碱基T(R)作用区域顺序为3=4>1>2;水溶液使Fapy-G与碱基C作用的结合能减弱程度最大,结合能达到41.84~58.58 k J/mol,且使碱基对结合能力次序发生改变.     

乙炔与H2O分子之间的氢键结构与性质

采用B3LYP方法,在6-311++G水平上优化得到了H2O…C2H2氢键复合物的σ-n和H-π型两种稳定构型,并进行频率分析,讨论了相关自然键红外振动光谱的红移现象.采用NBO理论对σ-n和H-π氢键复合物形成过程中的电荷转移的类型进行了分析讨论.分子间的氢键相互作用能结果表明,σ-n型比H-π型氢键复合物更稳定.     

准确快速计算含多肽酰胺和核酸碱基的氢键复合物的氢键强度和氢键作用势能曲线

N-H···O=C、C-H···O=C、N-H···N和C-H···N等氢键作用是蛋白质a-螺旋结构、b-折叠结构和DNA双螺旋结构形成的主要因素,在生物分子识别、蛋白质复制以及遗传信息传递等过程中起重要作用。准确快速计算生物体系中存在的N-H···O=C、C-H···O=C、N-H···N和C-H···N等氢键作用强度以及氢键强度随分子几何结构(距离和角度)变化的势能曲线对正确模拟(从而正确认识和理解)蛋白质折叠机制和DNA双螺旋结构形成机制等生物过程意义重大,对设计合成具有特殊功能的生物分子材料有重要指导价值。本文主要介绍了近年来建立的偶极-偶极氢键作用模型及其在快速预测多肽-多肽分子间和核酸碱基-多肽分子间氢键作用强度和氢键作用势能曲线方面的应用。     

HCCO自由基与LiX(X=F,Cl,Br)锂键复合物的从头算研究

在MP2/6-311++G**水平上优化乙烯酮自由基与LiX(X=F,Cl,Br)形成锂键复合物.当卤素的电负性很强(如F元素),使得Li原子处于缺电子状态,此时,电子给体会把电子偏移向锂,形成共价性较强的锂键.而当卤素的电负性减弱时,锂键中主要成分逐渐变为离子键,并且此时锂键性质还要受电子给体影响.另外,由于HCCO为缺电子结构,电负性较弱且体积较大的卤素中的孤对电子会与HCCO之间通过静电相互作用,使得HCCO…Li—X键夹角变小,接近120°.锂键性质对HCCO…LiX(X=F,Cl,Br)复合物中Li—X的伸缩振动频率有直接影响.当锂键表现为共价性时,该频率红移,而当锂键表现为离子性时,该频率蓝移.但是,由于Cl的电负性与O的接近,C的电负性与Br接近所以,在O…Li…Cl和C…Li…Br中容易形成共振结构,导致远大于在其他复合物中的红移.     

核酸碱基与甘氨酸二肽间氢键作用的最佳位点

优化得到了碱基腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶与甘氨酸二肽分子形成的28个氢键复合物的稳定结构并计算了结合能, 探讨了五种碱基与甘氨酸二肽分子间氢键作用的最佳位点. 本文研究发现: 每种碱基均可以通过不同位点与二肽分子形成氢键复合物, 腺嘌呤、胸腺嘧啶、尿嘧啶、鸟嘌呤及胞嘧啶分别最倾向使用A3、T1、U1、G3及C1位点与甘氨酸二肽分子形成氢键复合物; 碱基分子某位点的质子化反应焓变越负所形成的氢键复合物越稳定, 去质子化反应焓变越小所形成的氢键复合物越稳定; 由氢键复合物的结合能计算得到的稳定性次序与由碱基分子质子化和去质子化反应焓变推得的稳定性次序一致.     

槲皮素与腺嘌呤氢键作用的最佳位点

本文优化得到了16个由槲皮素与腺嘌呤形成的氢键复合物的稳定结构,并计算了它们的结合能.研究发现,在气相和水相中,槲皮素均通过qu1位点与腺嘌呤作用形成稳定的氢键复合物.比较了腺嘌呤与槲皮素形成的氢键复合物、腺嘌呤与胸腺嘧啶形成的Watson-Crick碱基对的相对稳定性.在气相条件下Watson-Crick碱基对更稳定,在水相条件下腺嘌呤与槲皮素形成的氢键复合物更稳定,说明水相条件下腺嘌呤与槲皮素之间的相互作用强于与胸腺嘧啶之间的相互作用.基于标准反应Gibbs自由能变的计算结果估算了水相条件下腺嘌呤与槲皮素形成的氢键复合物和Watson-Crick碱基对的相对平衡浓度.     

Ab initio study of the structure,cooperativity, and vibrational properties in the mixed hydrogen‐bonded trimers of hydrogen isocyanide and water

The five trimers of H₂O···HNC···H₂O, H₂O···H₂O···HNC, HNC···H₂O···H₂O, H₂O···HNC···HNC, and HNC···HNC···H₂O have been studied with quantum chemical calculations. Their structures, harmonic vibrational frequencies and interaction energies have been calculated at the B3LYP and MP2 levels with the aug‐cc‐pVDZ and aug‐cc‐pVTZ basis sets. The cooperative effect on these properties has also been studied quantitatively. For HNC:(H₂O)₂ systems, the cyclic H₂O···H₂O···HNC trimer is most stable with an interaction energy of ?16.01 kcal/mol and a large cooperative energy of ?3.25 kcal/mol at the MP2/aug‐cc‐pVTZ level. For H₂O:(HNC)₂ systems, the interaction energy and cooperative energy in the H₂O···HNC···HNC trimer are larger than those in the HNC···HNC···H₂O trimer. The NH stretch frequency has a blue shift for the terminal HNC molecule in the HNC···H₂O···H₂O and HNC···HNC···H₂O trimers and a red shift in other cases. A many‐body analysis has also been performed to understand the interaction energies in these hydrogen‐bonded clusters. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Triangular halogen trimers. A DFT study of the structure, cooperativity, and vibrational properties

Triangular halogen trimers (RX)3, where X = Br, I and R represents H, H3C, H2FC, HF2C, F3C, CH2=CH, CH[triple bond]C, and Ph, have been investigated using the density functional theory in the Perdew, Burke, and Ernzerhof method. We report herein the optimized geometries of the stable structures, their vibrational frequencies, and binding energies with the two- and three-body terms. All trimer structures possess a cyclic array of halogen atoms in the type II approach by virtue of the nonspherical atomic charge distribution around the halogens. The Br...Br interactions in trimers are very weak, whereas the I...I interactions in trimers are relatively stronger. While all bromine trimers and most of iodine trimers are predicted to be noncooperative, three of iodine trimers show weak cooperativity. The analysis of vibration modes reveals that all halogen trimers exhibit no especially remarkable frequency shifts. It is also shown that the electrostatic contribution plays a major role in the halogen...halogen interactions in halogen trimers. In contrast to bromine trimers, the relative contribution of charge-transfer component to the halogen...halogen interactions becomes more important for iodine trimers.     

Ab initio study of the vibrational spectra of the hydrogen-bonded nitric acid dimer.

The changes in the vibrational characteristics characterizing the dimerization of nitric acid have been investigated by ab initio calculations at the MP2 level, with 6-31G(d,p) and 6-31 + G(d,p) basis sets, and B3LYP/6-31G(d,p) calculations. The most consistent agreement between the computed values of the frequency shifts for the planar fully symmetric structure (2A) and those experimentally observed suggests that this structure is preferred. It was established that the most sensitive to the complexation is the stretching O-H vibration. The values of the frequency shift (-306 cm(-1)) is indicative for the formation of the relatively strong hydrogen bonds. The calculations predict an increase of the infrared intensity of the stretching O-H vibration in the nitric acid dimer more than 26 times.     

Ab initio SCF MO study of the hydrogen-bonded dimers of cyanogen with HF and HCl

Large-scale gaussian orbital SCF MO calculations are presented for the hydrogen-bonded complexes NCCN…HF and NCCN…HCl. Calculated equilibrium geometries, hydrogen-bond dissociation energies and selected one-electron properties are given to supplement available experimental data. Changes of electron distribution on complex formation are discussed in terms of Mulliken population indices.     

Ab initio study of L-orientational defect in the hydrogen-bonded pattern of liquid water

Current status of the theory of orientational defects in H-bonded pattern of liquid water is briefly reviewed. Ab initio calculated water clusters from dimer to heptamer are thoroughly analyzed in terms of H-bonded pattern. New water heptamer structure of norbornane type is found via ab initio HF 6-311G** calculation. Its normal vibrations are properly assigned. Two ab initio water hexamers that refer to orientational defects are revealed. This is the first evidence of ab initio orientational defect in H-bonded patterns of water clusters. Some properties including normal vibrations of these defects, are studied.     

Ab initio study on the structures, energies, and vibrational frequencies of acetone complexes with metal monocations and dications

The structures, interaction energies, and vibrational frequencies of acetone and its complexes with proton and various metal monocations/dications such as H⁺, Li⁺, Na⁺, K⁺, Rb⁺, Cs⁺, Be²⁺, Mg²⁺, Ca²⁺, and Zn²⁺ have been investigated from the ab initio calculations. The linear or bent structure for the cation-acetone complexes has been found to be deeply influenced by the amount of charge transfer. The amount of red-shift of CO stretching frequencies of acetone is almost equivalent regardless of the kind of alkali metal monocations. This behavior can be attributed to the electrostatic nature of the interactions rather than orbital interactions. This phenomenon has been supported by the highest occupied molecular orbitals of the acetone complexes investigated, charge transfer based on natural bond orbital (NBO) atomic charges, and the inversely proportional behavior of the interaction energies to the interatomic distances r(M⁺O).     

Ab initio study of the structure, force fields, and vibrational spectra of alkali metal tellurites

The equilibrium geometrical parameters, force fields, vibration frequencies, and intensities in the IR molecular spectra of M₂TeO₃ (M-Li, Na, K) are found by the Hartree — Fock ab initio method in the bases complemented with polarization and diffusion functions using relativistic effective core potentials. The relative energies of alternative molecular configurations and the energies of the dissociation M₂TeO₃→ M2O + TeO2 are refined using Möller — Plessett second-order perturbation theory. It is found that the chemical bond in these molecules can be approximated by the scheme (M⁺)₂[TeO₃]^2- The equilibrium nuclear configuration has C_s symmetry and corresponds to the bis-bidentate coordination of M₊ cations by the pyramidal anion [TeO₃]^2-. The mono-bidentate structures of C_s symmetry correspond to the saddle points on the potential energy surface of the molecules. The calculated IR molecular spectrum of K₂TeO₃ agrees with the experimental spectrum obtained by the matrix isolation method by J. Ogden et al. [J. Chem. Soc. Dalton Trans., 1957 (1997)]. Based on the results of the previous ab initio studies of alkali metal sulfites and selenites, the barrier of the intramolecular rearrangement (bb)→(mb)→(bb)’ was found to lower in the series Li₂XO₃→Na₂XO₂→K₂XO₃ (X = S, Se, Te) and M₂TeO₃→M₂SeO₃→M₂SO₃. The degree of deformation of the XO fragment in the M₂XO₃ molecules decreases when the M atom is replaced by its heavier analog (Li→Na→K).     

Ab initio calculations of the vibrational force field and IR intensities of the ammonia dimers

Ab initio SCF calculations using the 4-31G basis set have been carried out to determine the equilibrium geometry, force constants and dipole moment derivatives of the linear (C_s and cylic (C_2h) ammonia dimers. The results are compared with monomer calculations and experimental data.     

Ab initio study of the structure,force field,and vibrational spectrum of the LiClO3 molecule

The geometrical structure and the vibrational spectrum of the LiClO₃ molecule are studied by the Hartree-Fock-Roothaan (HF) and configurational interaction (CI) methods taking into account single and double excitations and Davidson's correction for quadruple excitations. Double-zeta basis sets of Huzinaga-Dunning and McLean-Chandler complemented with polarization and diffuse functions are used. Potential surface sections corresponding to migration of the Li⁺ cation around the ClO ₃ ⁻ fragment are investigated. It was found that the LiClO₃ molecule has a single stable configuration of C_s symmetry with the bidentate coordination of the Li⁺ cation by the ClO ₃ ⁻ anion. The cyclic fragment is nonplanar (the dihedral angle θ(LiO₂Cl)=173°). The tridentate configuration of C_3v symmetry lies higher than the equilibrium configuration by 24.5 (HF) or 18.3 (CI) kJ/mole and is not an isomer. The ab initio force field of the molecule was refined by the scaling method. Some assignments of the IR bands of the matrix-isolated molecular forms existing in vapor over lithium chlorate are corrected. The vibration frequencies (cm⁻¹) and IR intensities (km/mole; in parentheses) are calculated with the refined force field: A′ type 1099(236), 856(81), 630(73), 557(119), 481(87), 156(66); A″ type 887(229), 459(35), 367(23). Ivanovo State Chemical Technological Academy. Translated fromZhurnal Strukturnoi Khimii, Vol. 37, No. 3, pp. 440–449, May–June, 1996.