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

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

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, π^*)垂直激发能.     

NCBBH…HNa(H2O)n和CNBBH…HNa(H2O)n(n=1~5)结构与性质的理论研究

为了探索缺电子B-H键作为质子供体形成双氢键复合物的溶剂化效应,分别采用DFT-B3LYP/6-311++G**和CCSD(T)/6-311++G**方法对NCBBH…HNa和CNBBH…HNa及其水合物NCBBH…HNa(H2O)n和CNBBH…HNa(H2O)n(n=1~5)进行了结构优化和相互作用能计算,并利用AIM(atom in molecule)方法分析了H…H键特征,借助前线分子轨道理论探讨了水合物中双氢键形成H-H共价键的本质。结果表明:随着H2O分子数的增加,B-H键拉长,H…H距离缩短,双氢键由离子型向共价型过渡;当H2O分子数达到4时,双氢键相互作用能和NCBBH…HNa与水分子间的相互作用能分别达到-374.21和-306.50 kJ.mol-1,形成了H-H共价键;缺电子B-H键作为质子供体形成双氢键复合物的水合物析出H2的能力比FH…HLi(H2O)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*)激发态的垂直激发能.     

乙醛二聚体内蓝移型氢键的理论研究

利用理论方法研究了乙醛二聚体内的氢键. 在MP2/6-31＋G(d), B3LYP/6-31＋G(d), B3LYP/6-311＋＋G(d,p)和B3LYP/6-311＋＋G(3df,2p)水平上, 利用常规方法和均衡校正方法对3种稳定的乙醛二聚体进行了几何优化和振动频率计算. 计算结果表明: 在二聚体A和C中乙醛中C—H键强烈收缩, 存在显著的C—H…O蓝移型氢键. 自然键轨道(NBO)分析表明, 电子供体轨道和电子受体轨道之间相互作用的稳定化能、分子内电子密度重排、轨道再杂化和结构重组是决定氢键红移和蓝移的主要因素. 其中, 轨道间稳定化能属于键伸长效应, 分子内电子密度重排、轨道再杂化和电子受体内部结构重组属于键收缩效应. 在二聚体A和C中, 由于键收缩效应处于优势地位导致C—H…O蓝移氢键存在.     

Al6(OH)18(H2O)6的结构及成键方位的从头算及密度泛函分析Ⅱ

用RHF法6-31G、6-31G*、6-31G**水平的量子化学从头算及RB3LYP/6-31G、B3LYP/6-31G*,B3LYP/6-31G**水平的密度泛函方法,采用SCRF=Dipole溶剂模型,计算了三水铝石有利生长基元最优结构Al6(OH)18(H2O)6的总能量、布居数、原子静电荷等.计算结果表明,桥联OH基团更易成键,Al6(OH)18(H2O)6较为有利的成键方位是桥联OH基团方位.     

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

用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势能参数提供了坚实的基础,此参数对分子动力学模拟天然气水合物是非常有用的.     

喹吖啶酮光物理性质的理论表征

运用密度泛函理论对喹吖啶酮的光物理性质进行了研究. 结果表明, B3LYP/6-31G(d)//TD-B3LYP-PCM/6-311＋G(d)和CIS/6-31G(d)//TD-B3LYP-PCM/6-311＋G(d)方法能分别较好地描述喹吖啶酮的光吸收和光发射性质. 在吸收和发射的电子跃迁过程中, 主要是π-π*跃迁, 电荷在N和O之间的转移, 既包括单分子内的电荷跃迁, 也包括了通过氢键或π-π作用形成的二聚体的电荷跃迁. 溶剂对其光物理性质是有影响的, 这种影响是通过改变其前线轨道能来实现的, 受影响最大的是最高占据分子轨道(HOMO).     

复合物C2H6·(H2O)2中分子间相互作用的理论研究

在B3LYP/6-311 G**水平上得到C2H6.(H2O)2复合物势能面上四种稳定构型。在相同基组下经MP2电子相关能和基组叠加误差(BSSE)进行单点能量校正,求得单体间相互作用能的大小。结果发现:四种稳定构型都通过CH…O氢键而形成,相应σ(CH)键都出现了较小的收缩,导致伸缩振动发生蓝移,在最稳定的复合物Comp lex2和Comp lex3中,H2O(A)分子的一个H原子与C2H6的两个H原子相对距离较短,并且具有最大的总相互作用能和两个单体AC分子间相互作用能,这说明三个氢原子间存在着较强的相互作用,并对分子的稳定性起着重要作用.     

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

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

2-羟基吡啶质子转移过程的理论研究

采用量子化学中的密度泛函理论,在B3LYP/6-31G(d)基组水平上,计算并考察了2-羟基吡啶分子醇式结构和酮式结构进行结构互变的质子转移过程中的4种可能途径:(a)分子内质子转移,(b)水助催化质子转移,(c)同种二聚体双质子转移和(d)异种二聚体间双质子转移.计算结果表明,途经c所需要的活化能最小(2.6 kJ•mol－1,逆反应则为27.1 kJ•mol－1),而过程a所需要的活化能最大(137.2 kJ•mol－1),途径b和d的活化能居中间(分别为38.7和17.3 kJ•mol－1).研究还表明,氢键在降低反应活化能方面起着重要的作用.     

Ab initio- and density-functional studies of conformational behaviour of N-formylmethionine in gaseous phase

The current work is a study of the conformational space of the non-ionic N-formylmethionine molecule around its seven structurally significant internal backbone torsional angles at B3LYP/6-31++G(d,p) levels of theory in the gaseous phase. The potential energy surface exploration reveals that a total of 432 different conformers would result if all the possible combinations of the internal rotations were to be considered. A set of twelve conformers of the N-formylmethionine molecule are then further analysed in terms of their relative stabilities, theoretically predicted harmonic vibrational frequencies, HOMO-LUMO energy gaps, ESP charges, rotational constants and dipole moments calculated using MP2/6-31++G(d,p) and B3LYP/6-311++G(d,p) levels. The calculated relative energy-range of the conformers at the MP2 level is 11.08 kcal mol^?1 (1 kcal = 4.1868 kJ), whereas the same obtained at the B3LYP level is 10.02 kcal mol^?1. The results of this study provide a good account of the role of four types of intramolecular H-bonds, namely O…H—O, O…H—N, O…H—C and N…H—C, in influencing the energies of the conformers as well as their conformational and vibrational spectroscopic aspects. The relative stability order of the conformers appears to depend on the level of theory used while the vibrational frequencies calculated at the B3LYP level are in better agreement with the experimental values.     

Hydration of DNA Bases and Compounds Containing Small Rings—A Model for Interactions of the Ricin Toxin A-Chain. A Theoretical Ab Initio Study

Ab initio quantum chemical calculations were performed for four neutral gas phase adenine and four pterin tautomers along with guanine and formycin bases. The water complexes of the lowest energy tautomers of these bases have been studied to mimic their interaction with the ricin toxin A chain (RTA). The water molecules create a full first hydration shell around the bases. Full geometry optimizations without any constraints on the planarity of these hydrated complexes were carried out at the HF/6-31G(d,p) level. Single point calculations were also performed at the correlated MP2/6-31G(d,p)//HF/6-31G(d,p) level of theory. Hydration energies were corrected for the basis set superposition error. Hydration energies of adenine and formycin are predicted to be lower (in magnitude) than those for the pterin and guanine. Due to these properties, two pterin tautomers can be considered as potentially useful inhibitors of RTA.     

Electronic transitions of guanine tautomers, their stacked dimers, trimers and sodium complexes

Planar and nonplanar geometries of the keto-N9H and keto-N7H tautomers of the guanine base of DNA as well as the hydrogen bonded complexes of these species with three water molecules each were optimized using the density functional theory at the B3LYP/6-31G** level. Geometries of the isolated bases were also optimized using the ab initio approach at the MP2/6-31G** level. The isolated keto-N9H and keto-N7H tautomers as well as their hydrogen bonded complexes with three water molecules each were solvated in bulk water employing the polarized continuum model (PCM) of the self-consistent reaction field theory (SCRF). Stacked dimers and trimers of both the tautomers of guanine were generated by placing the planar forms of the species at a fixed distance of 3.5 A from the neighboring one and rotating one molecule with respect to the other by 110 degrees for the keto-N9H form and 90 degrees for the keto-N7H form which corresponded to total energy minima at the B3LYP/6-31G** level. Geometry optimization for the cation of the monomer of guanine was performed at the same level of theory, and its solvation in bulk water was treated using the PCM model of the SCRF theory. The geometries of complexes of the two tautomers of guanine with a Na+ ion each were optimized at the B3LYP/6-31G** level, and the Na+ ion is predicted to bind with the keto-N9H tautomer preferentially. While the complex of the keto-N7H form of guanine with three water molecules in gas phase is slightly more stable than the corresponding complex of the keto-N9H form of guanine, the reverse is true in bulk water. Stacking interactions enhance the relative stability of the keto-N9H tautomer over that of the keto-N7H tautomer, suggesting that in bulk solutions, the former would be dominant. Electronic spectra of the isolated tautomers of guanine, those of their complexes with three water molecules each, the (keto-n9h and keto-n7h) cation of guanine, the complexes of the tautomers with a Na+ ion each, the stacked dimers and trimers of the two tautomers were calculated using configuration interaction involving single electron excitations (CIS). The relative absorption intensities of the two tautomers of guanine near 275 and 248 nm in the monomer, dimer, and trimer are predicated to be in the opposite order. Thus the absorption intensity oscillation observed using a guanine aqueous solution can be explained in terms of oscillation of relative populations of the two tautomers of the molecule. The 248 nm absorption peak would be appreciably red-shifted on formation of the cation of guanine. Binding of the Na+ ion with the two tautomers of guanine reduces intensities of their transitions appreciably and also it causes large red-shifts in the same.     

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     

Comprehensive study of the interaction between hydrogen halides and methanol derivatives

The R? CH₂? HO…H? X (R = SCl, Cl, SH, NO₂, OMe, CHO, CN, C₂H₅, CH₃, H; X = F, Cl, Br) complexes are considered here as the interest sample for the consideration of different measures of H‐bond strength. The intermolecular interaction energies are predicted by using MP2/6‐31++G(d,p) and B3LYP/6‐31++G(d,p) methods with basis set superposition error and zero‐point energy corrections. The results showed that intermolecular hydrogen bonds for complexes with HF are stronger than such interactions in complexes with HCl and HBr. Quantum theory of “Atoms in Molecules” and natural bond orbitals method were applied to analyzed H‐bond interactions. The gas phase thermodynamic properties of complexes were predicted using quantum mechanical computations. The obtained results showed a strong influence of the R and X substituents on the thermodynamic properties of complexes. Numerous correlations between topological, geometrical, thermodynamic properties and energetic parameters were also found. © 2011 Wiley Periodicals, Inc.     

Transitions to high tautomeric states can be induced in adenine by interactions with carboxylate and sodium ions: DFT calculation data

Twelve triple complexes of nine adenine tautomers with carboxylate ion of acetic acid CH₃COO⁻ and sodium ion Na⁺ of the CH₃COO⁻:Ade:Na⁺ type were studied by the B3LYP/6-311++G(d,p)//B3LYP/6-31G(d,p) quantum chemical method. It was established that three rare adenine tautomers generate more stable complexes than the ground state one. The B3LYP/6-311++G(d,p) full optimization was applied to four most stable triple complexes in order to refine on their stability order. To evaluate contributions of each ligands for stabilization of these complexes, calculations of energies and dipole moments of respective binary complexes of the CH₃COO⁻:Ade and Ade:Na⁺ types were performed at the B3LYP/6-31++G(d,p) level of theory. At the same level, there were computed energies and dipole moments of isolated tautomers. Two of them, first calculated in this work, appeared to be zwitter-ions. Highly cooperative interplay of interactions of both ligands with adenine tautomers was observed in triple complexes. Preferable sites of coordination of Na⁺ with tautomers were ascertained. Biological importance of the results obtained is discussed.     

An analytic potential energy function for the amide–amide and amide–water intermolecular hydrogen bonds in peptides

An analytic potential energy function is proposed and applied to evaluate the amide–amide and amide–water hydrogen‐bonding interaction energies in peptides. The parameters in the analytic function are derived from fitting to the potential energy curves of 10 hydrogen‐bonded training dimers. The analytic potential energy function is then employed to calculate the N? H…O?C, C? H…O?C, N? H…OH₂, and C?O…HOH hydrogen‐bonding interaction energies in amide–amide and amide–water dimers containing N‐methylacetamide, acetamide, glycine dipeptide, alanine dipeptide, N‐methylformamide, N‐methylpropanamide, N‐ethylacetamide and/or water molecules. The potential energy curves of these systems are therefore obtained, including the equilibrium hydrogen bond distances R(O…H) and the hydrogen‐bonding energies. The function is also applied to calculate the binding energies in models of β‐sheets. The calculation results show that the potential energy curves obtained from the analytic function are in good agreement with those obtained from MP2/6‐31+G** calculations by including the BSSE correction, which demonstrate that the analytic function proposed in this work can be used to predict the hydrogen‐bonding interaction energies in peptides quickly and accurately. © 2009 Wiley Periodicals, Inc. J Comput Chem, 2009