三硝基苯-对位取代苯酚负离子荷移复合物从头算研究

运用G94W量子化学程序,在HF/3-21G基组水平上对三硝基苯-对位取代苯酚(取代基:CH~3O-,CH~3-,Cl-)负离子的电荷转移复合物进行从头计算。把电荷转移复合物看成一个超分子,研究该系列复合物的稳定性、电荷转移及几何构型等规律。计算结果表明,苯酚负离子供电中心O沿一倾角指向三硝基苯中的任意相邻两个硝基之间的C位置上,复合物的稳定性按对位取代苯酚取代基CH~3O-,CH~3-,Cl-的次序减小,与实验规律相一致。计算结果还表明,形成该系列复合物具有明显的电荷转移,其稳定性与电荷转移量有关等。     

低阶煤中氢键作用的色散校正密度泛函理论研究

本研究以苯酚…苯酚、苯酚…苯、苯酚…二苯醚、苯酚…喹啉和苯甲酸…苯甲酸为对象，采用色散校正的密度泛函理论分别研究褐煤中自缔合OH、OH-π、OH-醚O、OH-N和COOH-COOH之间形成的氢键。此外，还研究了氢键供体中取代基(CH₃-、CH₃O-、OH-、NH₂-、COOH-和NO₂-)对氢键的影响。对上述复合物进行了几何优化，并计算了能量、Mulliken电荷分布及振动频率。从优化的结构中可以看出上述复合物之间都存在氢键，所有复合物中O-H键键长都比苯酚中自由羟基的长，这表明这些复合物之间存在相互作用。其中，羧酸…羧酸复合物中O-H键的键长最长。此外，通过Mulliken电荷分布可看出上述复合物之间存在电荷转移。基于振动频率分析，所有的O-H键伸缩振动都发生了红移，尤其是羧酸…羧酸和苯酚…喹啉复合物，这可为煤中羟基振动的红外光谱分析提供依据。根据键能不同氢键强度按以下顺序依次递减：COOH-COOH>OH-N > 自缔合OH≈OH-醚O > OH-π，这与振动频率的分析结果一致。此外，不同取代基对氢键作用的影响不同。     

精氨酸侧链和核酸碱基间离子氢键作用强度分析

采用MP2/6-31+G(d,p)方法优化得到了22个由精氨酸侧链与碱基尿嘧啶、 胸腺嘧啶、 胞嘧啶、 鸟嘌呤及腺嘌呤形成的氢键复合物的气相稳定结构, 使用包含BSSE校正的MP2/aug-cc-pVTZ方法计算得到了复合物的气相结合能, 通过MP2/6-31+G(d,p)方法和PCM模型优化得到了复合物的水相稳定结构, 采用MP2/aug-cc-pVTZ方法和PCM模型计算得到了复合物的水相结合能. 研究发现, 精氨酸侧链与碱基间的离子氢键作用强度与单体间电荷转移量、 氢键临界点电子密度及二阶作用稳定化能密切相关. 与中性氢键相比, 离子氢键作用具有更显著的共价作用成分. 研究还发现, 精氨酸侧链和碱基间形成的氢键复合物的稳定性次序可以通过氢键受体碱基分子上氧原子和氮原子的质子化反应焓变进行预测, 质子化反应焓变越负, 形成的氢键复合物越稳定.     

水溶液中氨基酸侧链对G∶C碱基对间氢键影响的理论研究

利用量子化学方法研究了气相和水溶液下,氨基酸侧链与鸟嘌呤和胞嘧啶间的氢键作用.应用B3LYP/6-31+G(d,p)方法优化复合物几何结构,使用MP2/aug-cc-p VDZ方法进行复合物能量、自然键轨道(NBO)电荷和二阶稳定化能的计算.结果表明,水溶液可使氨基酸侧链与碱基或碱基对之间氢键键能显著减小;带电复合物气相和水溶液氢键键能之差范围为50.63~146.48 k J/mol,中性为0.17~24.94 k J/mol;电荷的转移量与氢键键能成正比,电荷转移量越多,复合物越稳定;二阶稳定化能与氢键键长成反比,与电荷转移量成正比,且气相与水溶液氢键二阶稳定化能之比约为两相的电荷转移量之比.水溶液对该类体系中氢键作用具有明显影响.     

芳香氮化物-CHCl3系列复合物密度泛函法研究

报道对芳香氮化物吡咯(C4H5N)、苯胺(C6H5NH2)、对位氯代吡啶(C5H4NCl)、吡啶(C5H5N)及吡咯负离子(C4H4N-)与三氯甲烷(CHCl3)形成的系列弱作用复合物的理论计算结果.把复合物看作是一个超分子,在密度泛函B3LYP/6-311G(d,p)的水平上进行计算,探讨该系列复合物的稳定性规律、电荷转移及主要几何参数变化等规律.计算结果表明:此类芳香氮化物与三氯甲烷(CHCl3)可形成一系列弱化学作用复合物,特别是带负电荷的芳香氮化物可与CHCl3形成较稳定的复合物.形成复合物的过程包含着电荷转移,该系列复合物的稳定性与广义H键距离的变化规律相一致,与电荷转移量的变化规律大体一致.     

电子给受体复合物中电荷转移吸收光谱和溶剂效应

用MP2/6-31G^**方法研究了二氯甲烷溶剂分子与电子给体、受体以及电子给受体复合物间的相互作用,结果表明,二氯甲烷与电子受体和电子给受体复合物间有弱氢键相互作用.利用CIS/6-31++G^**方法研究了溶剂与溶质分子间形成氢键对激发态的影响.自然键电荷分析表明,电子给受体复合物的S₀→S₁跃迁导致一个电子从电子给体转移到受体.结合非平衡溶剂化处理和自洽反应场方法研究了溶剂分子与复合物间形成氢键时的电荷转移吸收光谱.计算表明氢键作用导致复合物的电荷转移吸收光谱蓝移.     

1,2,3-三氮杂苯-水簇复合物氢键结构与性质的理论研究

在B3LYP/6-31 G**水平上对1,2,3-三氮杂苯和水形成1∶1,1∶2和1∶3复合物的基态氢键结构进行几何优化和性质计算,结果表明复合物之间存在较强的氢键作用。氢键的形成使水分子中H-O键振动频率减小(红移)。NBO分析表明,最稳定的1∶1,1∶2和1∶3复合物发生分子间电荷转移总量分别为0.0222e,0.0261e和0.0273e。     

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

在DFT-B3LYP/6-311++G**水平上分别求得(CH3)2S…ClOH卤键复合物和(CH3)2S…HOCl氢键复合物势能面上的稳定构型. 频率分析表明, 与单体HOCl相比, 在两种复合物中, 10Cl—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→σ*(10Cl—11O); (ii) 孤对电子LP(1S)2→σ*(10Cl—11O), 其中孤对电子LP(1S)2→σ*(10Cl—11O)转移占主要作用, 总的结果是使σ*(10Cl—11O)的自然布居数增加0.14035e, 同时11O原子的再杂化使其与10Cl成键时s成分增加, 即具有与电荷转移作用同样的“拉长效应”; 在(CH3)2S…HOCl氢键复合物中也存在类似的电荷转移, 但是11O原子的再杂化不同于前者. 自然键共振理论(NRT)进行键序分析表明, 在卤键复合物和氢键复合物中, 10Cl—11O和12H—11O键的键序都减小. 通过分子中原子理论(AIM)分析了复合物中卤键和氢键的电子密度拓扑性质.     

胺(或氨) - CHCl3系CT复合物密度泛函法研究

报道对氨(NH3母体)、甲胺(CH3NH2)、二甲胺[(CH3)2NH]、三甲胺[(CH3)3N]与三氯甲烷(CHCl3)形成的系列复合物的理论研究结果。把复合物看作一个超分子,在B3LYP/6-311G(d,p)的水平上进行密度泛函法计算,并运用完全均衡校正法进行基组超位误差校正(BSSE)。探讨该系列复合物的电子结构与相关性质,如稳定性、电荷转移及主要几何参数的变化规律等,结果表明,形成复合物的稳定性次序为NH3－CHCl3(Ⅰ)>CH3NH2-CHCl3(Ⅱ)>(CH3)2NH-CHCl3(Ⅲ)>(CH3)3N-CHCl3(Ⅳ)。形成复合物的过程包含着电荷转移。该系列复合物的稳定性与电荷转移量、前线轨道能量差△εL-H及广义H键距离良好的线性关系,与结合点的电荷布居也有密切的关系,复合物的稳定性是分子间共价作用和静电作用两方面因素综合的结果。计算结果能较好地解释有关的实验现象和规律。     

胺(或氨) - CHCl3系CT复合物密度泛函法研究

报道对氨(NH3母体)、甲胺(CH3NH2)、二甲胺[(CH3)2NH]、三甲胺[(CH3)3N]与三氯甲烷(CHCl3)形成的系列复合物的理论研究结果。把复合物看作一个超分子,在B3LYP/6-311G(d,p)的水平上进行密度泛函法计算,并运用完全均衡校正法进行基组超位误差校正(BSSE)。探讨该系列复合物的电子结构与相关性质,如稳定性、电荷转移及主要几何参数的变化规律等,结果表明,形成复合物的稳定性次序为NH3－CHCl3(Ⅰ)>CH3NH2-CHCl3(Ⅱ)>(CH3)2NH-CHCl3(Ⅲ)>(CH3)3N-CHCl3(Ⅳ)。形成复合物的过程包含着电荷转移。该系列复合物的稳定性与电荷转移量、前线轨道能量差△εL-H及广义H键距离良好的线性关系,与结合点的电荷布居也有密切的关系,复合物的稳定性是分子间共价作用和静电作用两方面因素综合的结果。计算结果能较好地解释有关的实验现象和规律。     

An ab initio study toward understanding the mechanistic photochemistry of acetamide

The potential energy surfaces for CH(3)CONH(2) dissociation into CH(3) + CONH(2), CH(3)CO + NH(2), CH(3)CN + H(2)O, and CH(3)NH(2) + CO in the ground and lowest triplet states have been mapped with DFT, MP2, and CASSCF methods with the cc-pVDZ and cc-pVTZ basis sets, while the S(1) potential energy surfaces for these reactions were determined by the CASSCF/cc-pVDZ optimizations followed by CASSCF/MRSDCI single-point calculations. The reaction pathways leading to different photoproducts are characterized on the basis of the computed potential energy surfaces and surface crossing points. A comparison of the reactivity among HCONH(2), CH(3)CONH(2), and CH(3)CONHCH(3) has been made, which provides some new insights into the mechanism of the ultraviolet photodissociation of small amides.     

On the cooperativity of cation-pi and hydrogen bonding interactions

Quantum chemical calculations are performed to gauge the effect of cation-pi and hydrogen bonding interactions on each other. M-phenol-acceptor (M = Li (+) and Mg (2+); acceptor = H(2)O, HCOOH, HCN, CH(3)OH, HCONH(2) and NH(3)) is taken as a model ternary system that exhibits the cation-pi and hydrogen bonding interactions. Cooperativity is quantified and the computed positive cooperativity between cation-pi and hydrogen bonding interactions is rationalized through reduced variational space (RVS) and charge analyses.     

甲醛与甲酰胺相互作用的从头算研究

The optimizations geometries and vibrational frequencies of H2CO,HCONH2 and acquired 3 complexes between H2CO?HCONH2 have been calculated by using the ab initio method at the MP2/6-31G（ d）and MP2 （FC）/6-311++G（d,p）level. The non-minimum structures with negative vibrational frequencies are excluded. The lowest energy conformer of these complexes is a cyclic structure with N - H?O and C - H?O hydrogen bonds in a common plane. No significant changes are observed in the geometries of the monomers in their complexed state. The most characteristic geometrical properties of the complex are the lengthening of the contacting N-H bonds by 0.4-1.1 pm,and the general shortening of the contacting C-H bonds by 0.3-0.4 pm with respect to the monomers. The interaction energies of complexes have been corrected by the basis set superposition error （BSSE）using the full Boys-Bernardi counterpoise correction scheme. The corrected complex interaction energies of 3 structures at MP2/6-311++G（2df,3p）/ / MP2（FC）/6-311++G（d,p）level are -29.94, -16.10 and -18.45 kJ/mol,respectively. The interaction energies indicate that C - H?O is a weak hydrogen bond. The results of natural bond orbital population analysis reveals that there is only a small charge-transfer in the process of forming the complexes. The results of natural bond orbital analysis and atom in the molecule scheme appear quite significant in view of their importance for understanding the mechanisms of intermolecular interaction leading to hydrogen bonding. The results of molecular interaction energy decomposition analysis show that the electrostatic interaction plays an essential role in stabilizing the H2CO?HCONH2 complexes.     

Models of fragmentations induced by electron attachment to protonated peptides

Invoking a number of theoretical levels ranging from HF/STO-3G to CCSD(T)/aug-cc-pVQZ, we have made a detailed survey of six potential energy surfaces (NH4+, NH4*, [CH3CONHCH3]H+, [CH3CONHCH3]H*, [HCONHCH2CONH2]H+, [HCONHCH2CONH2]H*). In conjunction with this, ab inito direct dynamics calculations have been conducted, tracing out several hundred reaction trajectories to reveal details of the electron-capture dissociation mechanism. The model calculations suggest the possibility of a bimodal pattern where some of the radicals, formed upon recombination, dissociate almost directly within one picosecond, the remaining radicals being subject to complete energy redistribution with a subsequent dissociation occurring at the microsecond timescale. Both processes give rise to c and z backbone fragments, resulting from cleavage of N-Calpha bonds of the peptide chain.     

CH3SH…HOO开壳型氢键复合物的结构及其电子密度拓扑性质

在DFT-B3LYP/6-311++G**水平下求得CH3SH…HOO复合物势能面上的稳定构型. 计算结果表明, 在HOO以其O8—H7作为质子供体与CH3SH分子中的S5原子为质子受体形成的氢键复合物1和2中, O8—H7明显被“拉长”, 且其伸缩振动频率发生显著的红移, 红移值分别为330.1和320.4 cm-1; 在CH3SH分子以其S5—H6作为质子供体与HOO的端基O9原子为质子受体形成的氢键复合物3和4中, 也存在类似的情况, 但S5—H6伸缩振动频率红移不大. 经MP2/6-311++G**水平计算的4种复合物含BSSE校正的相互作用能分别为-20.81, -20.10, -4.46和-4.52 kJ/mol. 自然键轨道理论(NBO)分析表明, 在CH3SH…HOO复合物1和2中, 引起H7—O8键长增加的因素包括两种电荷转移, 即孤对电子n1(S5)→σ*(H7—O8)和孤对电子n2(S5)→σ*(H7—O8), 其中后者为主要作用. 在复合物3和4中也有相似的电荷转移情况, 但轨道间的相互作用要弱一些. AIM理论分析结果表明, 4个复合物中的S5…H7间和O9…H6间都存在键鞍点, 且其Laplacian量▽2ρ(r)都是很小的正值, 说明这种相互作用介于共价键和离子键之间, 偏静电作用为主.     

Proton affinity correlations between hydrogen and dihydrogen bond acceptors

Several series of hydrogen- and dihydrogen-bonded complexes with HCN, C2H2, HF, H2O, CH3CONH2, and CH3COOH as donors and H2O, MeOH, EtOH, MeOMe, NH3, NH2Me, NHMe2, NMe3, NEtMe2, and BH3-NMe3 as acceptors were investigated using the MP2/6-311++G(d,p) level of theory. The total lowering of the X-H stretching frequencies in the hydrogen-bonded complexes were linearly correlated with the proton affinities of the accepting bases. From comparison of hydrogen- and dihydrogen-bonded complexes, a scaling factor to estimate the exact proton affinity of a dihydrogen bond acceptor was developed. Further, the scaling factor involving linear donors (1.204) is marginally higher than that involving nonlinear donor molecules (1.162). Finally, it was found that, given identical conditions, a hydrogen bond will be about 16-20% stronger than a corresponding dihydrogen bond.     

Vacuum ultraviolet (VUV) photodecomposition of urea isolated in cryogenic matrix: first detection of isourea

Vacuum ultraviolet (VUV) irradiation at wavelengths of lambda > 160 nm of urea-h4 (NH2CONH2) and urea-d4 (ND2COND2) has been monitored by Fourier transform infrared spectroscopy in argon and xenon matrixes. Several primary photoproducts, such as HNCO:NH3 (isocyanic acid:ammonia), CO:N2H4 (carbon monoxide:hydrazine) molecular complexes, and isourea (H2N(OH)C=NH), which is reported for the first time, were characterized. The assignment of complexes was achieved by co-depositing the pairs of respective species, whereas the isourea identification was based on the comparison between the experimental and theoretical (B3LYP) infrared spectra. Isourea is found in the argon matrix in its most stable (s-Z)-(E) configuration. It is an intermediate in the VUV decomposition process; its dehydration leads to the NH2CN:H2O complex. In the xenon matrix, the photochemistry of urea yields the HNCO:NH3 complex as a major product, whereas the CO:N2H4 complex is observed in trace amounts. The observed differences between the argon and xenon matrixes suggest the crossing between S1 and T1 potential surfaces of urea to be responsible for the formation of the HNCO:NH3 complex. A comparison is also performed with other carboxamides, such as formamide (HCONH2) or acetamide (CH3CONH2).     

气相中开壳型(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, 说明这种相互作用介于共价键和离子键之间, 偏于静电作用.