1,2,4-三氮杂苯-(H2O)n复合物氢键相互作用的密度泛函理论研究

用密度泛函理论方法在B3LYP/6-31++G**水平上对1,2,4-三氮杂苯-(H₂O)_n (n＝1, 2, 3)氢键复合物的基态进行了结构优化和能量计算, 结果表明复合物之间存在较强的氢键作用, 所有稳定复合物结构中形成一个N…H—O氢键并终止于弱O…H—C氢键的氢键水链的构型最稳定. 同时, 用含时密度泛函理论方法(TD-DFT)在TD-B3LYP/6-31++G**水平上计算了1,2,4-三氮杂苯单体及其氢键复合物的单重态第一¹(n, π^*)垂直激发能.     

乙醇-水分子团簇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*)激发态的垂直激发能.     

含有二维蜂窝层/网格结构的新颖氢键型三维超分子体系

在液相反应体系中设计合成了两个新型超分子配合物 [Cu(H₂dhbd)(3-pyOH)(H₂O)]·3-pyOH·2H₂O (1) 和 [Cu₂(dhbd)(dpa)₂(H₂O)]·6H₂O (2) (其中H₄dhbd = 2,3-二羟基丁二酸，3-pyOH = 3-羟基吡啶， dpa = 2,2′-二吡啶胺)，并通过元素分析、紫外与红外光谱、热分析及X射线单晶衍射等技术对其进行了表征. 晶体结构分析表明，在四对强O···H−O氢键的导向作用下，配合物1的环形双核单元与紧邻的四个单元相连形成了封有游离的3-羟基吡啶的二维蜂窝层. 有趣的是，在二维蜂窝层间发现了作为氢键桥的双核水簇，其可进一步将上述二维蜂窝层连接为氢键型三维结构.配合物2则含有由手性双核单元通过O···H−O和O···H−N 氢键构筑的二维网格, 这些二维网格又通过其他结晶水的氢键连接扩展为具有隧道的三维超分子结构. 磁性研究结果表明，两个配合物均具有弱的反铁磁相互作用.     

N,N-二甲基乙酰胺与醇分子间相互作用的理论研究

采用密度泛函理论在B3LYP/6-311＋G*水平上对乙醇、丙醇、丁醇、戊醇与N,N-二甲基乙酰胺形成的1∶1氢键复合物进行计算研究. 结果表明: 醇与N,N-二甲基乙酰胺形成的复合物存在强的氢键, 表现为羰基氧原子的孤对电子与醇羟基反键σ轨道的相互作用. 振动分析显示, 分子间C＝O…H—O氢键的形成使C＝O和H—O伸缩振动频率明显红移. 溶剂对氢键产生较大的影响, 随着溶剂极性的增加, 复合物氢键有蓝移趋势.     

煤炭燃烧过程中N2O消除反应机理的密度泛函理论研究

用密度泛函理论B3LYP方法对煤炭燃烧过程中N₂O的消除反应进行研究。选用6-311++G**和aug-cc-pVTZ基组,优化了反应通道上反应物、过渡态和产物的几何构型。预测了它们的热力学性质(总能量、焓、熵和吉布斯自由能)及其随温度的变化。预测N₂O+CO反应的活化能为200 kJ·mol^-1,与实验值193±8 kJ·mol^-1较一致。计算了500~1 800 K 温度范围的反应速率常数。在N₂O的分解中,N₂O与H和CN自由基的反应为动力学优先进行的反应,其活化能为50~55 kJ·mol^-1。在B3LYP/aug-cc-pVTZ level水平下,N₂O+CN反应是热力学最有利的自发反应,其吉布斯自由能变化为-407 kJ·mol^-1。     

单电子卤键复合体系H3C···Br-Y(Y=H, CN, NC, CCH, C2H3)的作用机制

对单电子溴键复合物H₃C···Br—Y(Y=H, CCH, CN, NC, C2H3)的结构与性质进行了理论研究. 在B3LYP/6-311++G**水平上计算了稳定构型并做了频率分析. BSSE矫正的相互作用能(E^BSSE)和NBO及AIM分析输入的波函数在MP2/6-311++G**水平下完成. 复合物H₃C···Br—Y中, CH₃(供电子体)自由基均提供一未成对电子与Br—Y中Br(受电子体)形成了单电子溴键, 此单电子溴键也具有“三电子”键的特征. 单电子溴键的形成导致甲基H的背向Y弯曲和Br—Y键的拉长及红移单电子溴键复合物的产生. 考察了电子受体中不同取代基, C(spⁿ)-Br杂化及溶剂的存在对复合物作用的影响, 将单电子氢键, 单电子卤键和单电子锂键的作用强度做了对比, 进一步对Popelier提出的氢键体系中的前三个重要拓扑指标在单电子溴键体系中的重现性进行了探讨.     

气相中O3与HSO自由基间的氢键复合物

气相中O3与HSO自由基之间的相互作用及其反应在大气化学中非常重要. 在DFT-B3LYP/6-311++G**和MP2/6-311++G**水平上求得O3+HSO复合物势能面上的稳定构型, B3LYP方法得到了三种构型(复合物I, II和III), 而MP2方法只能得到一种构型(复合物II). 在复合物I和III中, HSO单元中的1H原子作为质子供体, 与O3分子中的端基O原子作为质子受体相互作用, 形成红移氢键复合物; 而在复合物II中, 虽与复合物I和III中具有相同的质子供体和质子受体, 却形成了蓝移氢键复合物. B3LYP/6-311++G**水平上计算的单体间相互作用能的计算考虑了基组重叠误差(BSSE)和零点振动能(ZPVE)校正, 其值在-3.37到-4.55 kJ·mol-1之间. 采用自然键轨道理论(NBO)对单体间相互作用的本质进行了考查, 并通过分子中原子理论(AIM)分析了三种复合物中氢键的电子密度拓扑性质.     

N-H…O红移氢键和蓝移氢键的理论研究

用理论方法研究了复合物HCHO…HNO, HCOOH…HNO, HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F分子间氢键. 在MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/ 6-31+G(d,p)和B3LYP/6-311++G(d,p)水平上, 利用标准方法和均衡校正方法对6种复合物进行了几何优化和振动频率计算. 计算结果表明: 在复合物HCHO…HNO和HCOOH…HNO中, HNO中N—H键强烈收缩, 存在显著的N—H…O蓝移氢键. 在复合物HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F中, NH3和NH₂F中N—H键伸长, 存在N—H…O红移氢键. 利用自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 电子密度重排、轨道再杂化和结构重组属于键收缩效应. 在复合物HCHO…HNO和HCOOH…HNO中, 由于键收缩效应处于优势地位导致N—H…O蓝移氢键存在. 在复合物HCHO…NH₃, HCOOH…NH₃, HCHO…NH₂F和HCOOH…NH₂F中, 键伸长效应居于主导地位,N—H…O红移氢键出现.     

溶剂对三(苯并咪唑-2-甲基)胺和间苯二胺四乙酸的铜配合物晶体结构的影响

本文利用三(苯并咪唑-2-甲基)胺和间苯二胺四乙酸为配体与硝酸铜在CH₃COCH₃/CH₃OH/H₂O混合溶液中反应得到配合物[Cu(ntb)(H₂mpda)]·0.5CH₃COCH₃·2H₂O (1), 在DMF/CH₃OH/H₂O混合溶液中反应得到配合物[Cu(ntb)(H₂mpda)]·DMF·CH₃OH·2H₂O (2)(ntb=三(苯并咪唑-2-甲基)胺, H₄mpda=间苯二胺四乙酸)。2个配合物的中心的铜离子分别与1个ntb配体的4个氮原子和1个H₂mpda的氧原子配位形成三角双锥的配位构型。受溶剂的影响, 配合物中配体的相对位置和构象有较大的区别, 配合物2的配位构型更加扭曲。两个配合物均通过氢键连接形成不同的复杂三维网络。     

DFT Study of Hydrogen-Bonded 1,3,5-Triazine-Water Complexes

The 1,3,5-triazine-water hydrogen bonding interactions have been investigated using the density functional theory B3LYP method and 6-31 ＋＋G^＊＊ basis, obtaining one, two and seven energy minima of the ground states for the 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively. The fully optimized geometries and binding energies were reported for the various stationary points. The global minima of 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes have a hydrogen bond N…H-O and a chain of water molecules, terminated by a hydrogen bond O…H-C. The binding energies are 13.38, 39.52 and 67.79 kJ/mol for the most stable 1,3,5-triazine-water, 1,3,5-triazine-（water）2 and 1,3,5-triazine-（water）3 complexes respectively, after the basis set superposition error and zero point energy corrections. The H-O symmetric stretching modes of water in the complexes are red-shifted relative to those of the monomer water. In addition, the NBO analysis indicates that inter-molecule charge transfer is 0.02145 e, 0.02501 e and 0.02777 e for the most stable 1 ： 1, 1 ： 2 and 1 ： 3 complexes between 1,3,5-triazine and water, respectively.     

Blue‐shifted and red‐shifted hydrogen bonds: Theoretical study of the CH3CHO· · ·HNO complexes

The blue‐shifted and red‐shifted H‐bonds have been studied in complexes CH₃CHO…HNO. At the MP2/6‐31G(d), MP2/6‐31+G(d,p) MP2/6‐311++G(d,p), B3LYP/6‐31G(d), B3LYP/6‐31+G(d,p) and B3LYP/6‐311++G(d,p) levels, the geometric structures and vibrational frequencies of complexes CH₃CHO…HNO are calculated by both standard and CP‐corrected methods, respectively. Complex A exhibits simultaneously red‐shifted C? H…O and blue‐shifted N? H…O H‐bonds. Complex B possesses simultaneously two blue‐shifted H‐bonds: C? H…O and N? H…O. From NBO analysis, it becomes evident that the red‐shifted C? H…O H‐bond can be explained on the basis of the two opposite effects: hyperconjugation and rehybridization. The blue‐shifted C? H…O H‐bond is a result of conjunct C? H bond strengthening effects of the hyperconjugation and the rehybridization due to existence of the significant electron density redistribution effect. For the blue‐shifted N? H…O H‐bonds, the hyperconjugation is inhibited due to existence of the electron density redistribution effect. The large blue shift of the N? H stretching frequency is observed because the rehybridization dominates the hyperconjugation. © 2006 Wiley Periodicals, Inc. Int J Quantum Chem, 2006     

Hydrogen Bonding Interaction of Formic Acid-, Formaldehyde-, Formylfluoride-Nitrosyl Hydride： Theoretical Study on the Geometries, Interaction Energies and Blue- or Red-Shifted Hydrogen Bonds

The hydrogen bonding interaction of formic acid-, formaldehyde-, formylfluoride-nitrosyl hydride complexes was investigated by the density functional theory （DFT） and ab inito method in conjunction with 6-311＋＋G（2d,2p） basis set. The geometries, vibrational frequencies and interaction energies of the complexes were calculated by both standard and CP-corrected methods respectively. Moreover, G3B3 method was employed to estimate the interaction energies. There are C--H…O, N--H…O, N--H…F blue-shifted H-bonds and red-shifted O----H…O H-bond in the complexes. Electron density redistribution and rehybridization contribute to the N--H and C--H blue shifts. All geometric reorganizations contribute to the N--H blue shifts and partial geometric reorganizations contribute to the C--H blue shifts. The geometric reorganizations of the complex C except ZH（5）-O（4）-C（1） contribute to the O----H red shift. For the N--H blue shifts, the effect of r（N--O） variation on the N--H blue shifts is larger than that of ZH-N-O variation. Rehybridization plays a dominant role in the degree of N--H blue shifts, whereas the electron density redistribution contributes more to the degree of C--H blue shifts than the other effects do.     

Density Functional Theory Study of Hydrogen Bonds of Bipyridine with 1,3,5-Benzenetricarboxylic Acid

The hydrogen‐bonded dimer and trimer formed between 1,3,5‐benzenetricarboxylic acid and bipyridine have been investigated using a density functional theory (DFT) method and 6‐31++G** basis set. The interaction energies are ?45.783 and ?89.998 kJ·mol^?1 for the most stable dimer and trimer, respectively, after the basis set superposition error and zero‐point corrections. The formation of O–H...N hydrogen bonds makes O–H symmetric stretching modes in the dimer and trimer red‐shifted relative to those of the 1,3,5‐benzenetricarboxylic acid monomer. The natural bond orbit analysis shows that the inter‐molecular charge transfers are 0.60475e and 1.20225e for the dimer and trimer, respectively. Thermodynamic analysis indicates that the formation of trimer is an exothermic and spontaneous process at low and room temperature. A supramolecule can be constructed through the strong N···H–O intermolecular hydrogen bonds between bipyridine and 1,3,5‐benzenetricarboxylic acid, which is in good agreement with the experimental results.     

Chiral discrimination in hydrogen‐bonded complexes of 2‐methylol oxirane with hydrogen peroxide

A systematic quantum chemical study reveals the effects of chirality on the intermolecular interactions between two chiral molecules bound by hydrogen bonds. The methods used are second‐order Møller–Plesset perturbation theory (MP2) with the 6‐311++g(d,p) basis set. Complexes via the O? H···O hydrogen bond formed between the chiral 2‐methylol oxirane (S) and chiral HOOH (P and M) molecules have been investigated, which lead to four diastereomeric complexes. The nomenclature of the complexes used in this article is enantiomeric configuration sign corresponding to English letters. Such as: sm, sp. The relative positions of the methylol group and the hydrogen peroxide are designated as syn (same side) and anti (opposite side). The largest chirodiastaltic energy was ΔE_chir = ?1.329 kcal mol^?1 [9% of the counterpoise correct average binding energy D_e(corr)] between the sm‐syn and sp‐anti in favor of sm‐syn. The largest diastereofacial energy was ?1.428 kcal mol^?1 between sm‐syn and sm‐anti in favor of sm‐syn. To take into account solvents effect, the polarizable continuum model (PCM) method has been used to evaluate the chirodiastaltic energies, and diastereofacial energies of the 2‐methylol oxirane···HOOH complexes. The chiral 2,3‐dimethylol oxirane (S, S) is C₂ symmetry which offers two identical faces. Hence, the chirodiastaltic energy is identical to the diastereomeric energy, and is ΔE_chir = 0.563 kcal mol^?1 or 5.3% of the D_e(corr) in favor of s,s‐p. The optimized structures, interaction energies, and chirodiastaltic energies for various isomers were estimated. The harmonic frequencies, IR intensities, rotational constants, and dipole moments were also reported. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Cooperative influence of water binding to peptides by N?H···OH2 and CO···HOH hydrogen bonds: Study by Ab Initio calculations

In this article, the geometry structures of hydrogen bond chains of formamide and N‐methylacetamide and their hydrogen‐bonded complexes with water were optimized at the MP2/6‐31G* level. Then, we performed Møller–Plesset perturbation method with 6‐311++g**, aug‐cc‐pvtz basis sets to study the cooperative influence to the total hydrogen bond energy by the N? H ··· OH₂ and C?O ··· HOH hydrogen bonds. On the basis of our results, we found that the cooperativity of the hydrogen‐bonded complexes become weaker as N? H ··· OH₂ and C?O ··· HOH hydrogen bonds replacing N? H ··· O?C hydrogen bonds in protein and peptide. It means that the N? H and C?O bonds in peptide prefer to form N? H ··· O?C hydrogen bond rather than to form C?O ··· HOH and N? H ··· OH₂. It is significant for understanding the structures and properties of the helical or sheet structures of protein and peptide in biological systems. © 2011 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

Theoretical Study on the Hydrogen Bonding Interactions in Complexes of 5‐Hydroxytryptamine with Water

The energies, geometries and harmonic vibrational frequencies of 1:1 5‐hydroxytryptamine‐water (5‐HT‐H₂O) complexes are studied at the MP2/6‐311++G(d,p) level. Natural bond orbital (NBO), quantum theory of atoms in molecules (QTAIM) analyses and the localized molecular orbital energy decomposition analysis (LMO‐EDA) were performed to explore the nature of the hydrogen‐bonding interactions in these complexes. Various types of hydrogen bonds (H‐bonds) are formed in these 5‐HT‐H₂O complexes. The intermolecular C4H5^5‐HT···O^w H‐bond in HTW3 is strengthened due to the cooperativity, whereas no such cooperativity is found in the other 5‐HT‐H₂O complexes. H‐bond in which nitrogen atom of amino in 5‐HT acted as proton donors was stronger than other H‐bonds. Our researches show that the hydrogen bonding interaction plays a vital role on the relative stabilities of 5‐HT‐H₂O complexes.     

Molecular Self‐Recognition: Rotational Spectra of the Dimeric 2‐Fluoroethanol Conformers

Fluoroalcohols show competitive formation of intra‐ and intermolecular hydrogen bonds, a property that may be crucial for the protein‐altering process in a fluoroalcohol/water solution. In this study, we examine the intra‐ and intermolecular interactions of 2‐fluoroethanol (FE) in its dimeric conformers by using rotational spectroscopy and ab initio calculations. Three pairs of homo‐ and heterochiral dimeric FE conformers are predicted to be local minima at the MP2/6‐311++G(d,p) level of theory. They are solely made of the slightly distorted most stable G+g?/G?g+ FE monomer units. Jet‐cooled rotational spectra of four out of the six predicted dimeric conformers were observed and unambiguously assigned for the first time. All four observed dimeric conformers have compact geometries in which the fluoromethyl group of the acceptor tilts towards the donor and ensures a large contact area. Experimentally, the insertion of the O? H group of one FE subunit into the intramolecular O? H???F bond of the other was found to lead to a higher stabilisation than the pure association through an intermolecular O? H???O? H link. The hetero‐ and homochiral combinations were observed to be preferred in the inserted and the associated dimeric conformers, respectively. The experimental rotational constants and the stability ordering are compared with the ab initio calculations at the MP2 level with the 6‐311++G(d,p) and aug‐cc‐pVTZ basis sets. The effects of fluorination and the competing inter‐ and intramolecular hydrogen bonds on the stability of the dimeric FE conformers are discussed.     

Time‐dependent density functional theory study on the electronic excited‐state hydrogen‐bonding dynamics of 4‐aminophthalimide (4AP) in aqueous solution: 4AP and 4AP–(H2O)1,2 clusters

The time‐dependent density functional theory (TDDFT) method has been carried out to investigate the excited‐state hydrogen‐bonding dynamics of 4‐aminophthalimide (4AP) in hydrogen‐donating water solvent. The infrared spectra of the hydrogen‐bonded solute?solvent complexes in electronically excited state have been calculated using the TDDFT method. We have demonstrated that the intermolecular hydrogen bond C? O···H? O and N? H···O? H in the hydrogen‐bonded 4AP?(H₂O)₂ trimer are significantly strengthened in the electronically excited state by theoretically monitoring the changes of the bond lengths of hydrogen bonds and hydrogen‐bonding groups in different electronic states. The hydrogen bonds strengthening in the electronically excited state are confirmed because the calculated stretching vibrational modes of the hydrogen bonding C?O, amino N? H, and H? O groups are markedly red‐shifted upon photoexcitation. The calculated results are consistent with the mechanism of the hydrogen bond strengthening in the electronically excited state, while contrast with mechanism of hydrogen bond cleavage. Furthermore, we believe that the transient hydrogen bond strengthening behavior in electroniclly excited state of chromophores in hydrogen‐donating solvents exists in many other systems in solution. © 2010 Wiley Periodicals, Inc. J Comput Chem, 2010     

Characterization of the N?H···ON and N?H···NO H‐bonds in nitrosamine dimers

Ab initio molecular orbital and DFT calculations have been carried out for three most stable dimers of parent nitrosamine (NA) in order to elucidate the structures and energetics of the dimers. The structures were optimized using HF, B3LYP, and MP2 methods with 6‐311+G(d,p) and 6‐311++G(2d,2p) basis sets. At the optimized geometries obtained at MP2/6‐311++G(2d,2p) level of theory, the energies were evaluated at QCISD/aug‐cc‐pVDZ and CCSD/aug‐cc‐pVDZ levels. The most stable dimer has two N? H···O?N hydrogen bonds and the least stable dimer has two N? H···N?O hydrogen bonds. The natural bond orbital analysis showed that the lpO(N) → BD*(N? N) and lpO(N) → BD*(N? H_b) interactions play a decisive role in the stabilization of the NH···O(N) hydrogen bonds in dimers. The atoms in molecules results reveal that the intermolecular N? H···O(N) H‐bonds in dimers have electrostatic character. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008