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

在MP2／6—31G(d)和MP2(FC)／6—311 G(d，p)水平上，对H2CO和HCOOH以及设计的4种构型H2CO…HCOOH复合物等进行几何全优化计算，经振动频率分析，确认它们为势能面上的稳定驻点．然后在MP2／6—311 G(2df，2p)水平上进行单点能计算和基组重叠误差(BSSE)校正以获得相互作用能，并利用自然键轨道理论探讨H2CO和H(X)OH相互作用的本质。     

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)下,考察了溶剂化效应.结果表明,所考察的水、二甲亚砜、乙醇和乙醚等四种溶剂均使单体间的相互作用能增大,且溶剂对复合物中的锂键结构及其振动频率具有显著的影响,而对复合物中的氢键的振动频率影响不大.     

反常蓝移单电子锂键Y…Li-CH3[Y=CH3, CH2CH3, CH(CH3)2, C(CH3)3]体系的结构与性质

在密度泛函理论B3LYP/6-311++G(d,p)及MP2/6-311++G(d,p)水平上研究了单电子锂键复合物Y…Li—CH3[Y=CH3, CH2CH3, CH(CH3)2, C(CH3)3]的结构与性质. 结果表明, 三种单电子锂键复合物H3CH2C…Li—CH3(II), (H3C)2HC…Li—CH3(III)和(H3C)3C…Li—CH3(IV)单电子锂键强度依II(-26.7 kJ·mol-1)相似文献   

复合物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分子间相互作用能,这说明三个氢原子间存在着较强的相互作用,并对分子的稳定性起着重要作用.     

CH_3SH+CN微观反应机理的理论研究

采用BMC-CCSD//B3LYP/6-311G(d,p)方法对CH3SH+CN反应机理进行了详细的理论研究.反应中涉及的各稳定点的构型、振动频率和零点能在B3LYP/6-311G(d,p)水平下计算得到,计算结果表明,该反应存在两种反应机理,5条可能的反应通道.SN2机理由于能垒太高,与直接氢抽提机理相比可以忽略.该反应的最可行通道为CN中的C原子进攻SH中的H原子经由一个前期和一个后期分子络合物生成产物CH3S和HCN.计算得到的反应焓变与已有实验值非常吻合.     

高级量子化学从头计算法研究N2和H2O分子间相互作用

在MP2/6-311++G(3d,3p)电子相关校正水平上,对N₂和H₂O分子间可能存在的氢键复合物进行全自由度能量梯度优化,发现了一个接近于直线的弱氢键总能量极小结构(1),进一步在高级电子相关校正的MP4SDTQ和CCSD(T)水平,用6-311++G(3d,3p)基组加上(3s3p2d1f)键函数,用MP4和CCSD(T)计算的结构1的结合能分别为-5.061kJ/mol和-4.715kJ/mol.     

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

三接触弯曲构型的蓝移氢键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伸缩振动频率的蓝移.     

CH_3CN+CH_3→CH_2CN+CH_4氢消除反应的直接动力学研究

采用直接动力学方法,对乙腈与甲基的反应进行了理论研究.在BHandHLYP/6-311G(d,p)和MP2/6-311G(d,p)水平下获得,稳定点的几何结构、振动频率及最小能量路径(MEP),在G3(MP2)和MC-QCISD水平下对能量信息进一步确认.利用正则变分过渡态理论,结合小曲率隧道效应校正(CVT/SCT)方法计算了该反应在220K～2000K的速率常数,与实验值符合得很好.     

三接触弯曲构型的蓝移氢键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伸缩振动频率的蓝移.     

Can cytosine, thymine and uracil be formed in interstellar regions? A theoretical study

This theoretical study investigates possible synthetic routes to cytosine, uracil and thymine in the gas phase from precursor molecules that have been detected in interstellar media. Studies at the CCSD(T)/6-311++G(d,p)//B3LYP/6-311++G(d,p) level of theory suggest that: The reactions between :CCCNH and :CCCO with monosolvated urea may constitute viable interstellar syntheses of cytosine and uracil. No low energy equilibration between cytosine and uracil has been demonstrated. The interaction of :CH(2) with the 5 C-H bond of uracil may form thymine in an energetically favourable reaction, but competing reactions where :CH(2) reacts with double bonds and other CH and NH bonds of uracil, reduce the effectiveness of this synthesis. The reaction between the hydrated propional enolate anion and isocyanic acid may produce thymine, in a reaction sequence where ΔG(reaction)(298 K) is -22 kJ mol(-1) and the maximum energy requirement (barrier to the first transition state) is only 47 kJ mol(-1).     

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

利用理论方法研究了乙醛二聚体内的氢键. 在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蓝移氢键存在.     

Theoretical investigation on stability and isomerizations of CH3SO isomers

The stability and isomerizations of CH3SO isomers have been investigated at B3LYP/6-311G(d,p), MP2/6-311G(d,p), QCISD/6-311G(d,p), and CCSD(T)/6-311G(d,p) levels. Geometries of isomers and transition states (TS) have been optimized at the B3LYP/6-311G(d, p) level. Vibration analysis and the intrinsic reaction coordinate (IRC) calculated at the same level have been applied to validate the connection of the stationary points. The four different methods give similar results: 11 isomers and 9 isomerization channels were found. CH3SO and CH2(S)OH are the most stable species among the 11 isomers. Furthermore, the breakage and formation of the chemical bonds in isomerization reactions have been discussed by the topological analysis method of electronic density. The "energy transition state (ETS)" and the "structure transition state (STS)" of all the isomerizations have been found. The topological analysis shows that the relative positions of ETS and STS are determined by reaction energy. The nonplanar four-member ring structure transition state (STS), which was first found in this paper, extended the concept of ring STS.     

C6H5—H…X分子间氢键的理论计算

利用密度泛函理论B3LYP方法, 在6-311+G(3df,2p)水平上对C6H5—H…X型分子间氢键进行了几何构型优化、氢键相互作用能、电子密度分布等计算. 其中C6H6为质子供体, HCOH、H2O、NH3、CH2NH和HCN为质子受体. 从电荷布居分析、自然键轨道等角度详细地讨论了C6H5—H…X 体系中, 共轭π键、O和N的不同键型结构对氢键形成的影响以及孤电子对与C—H 反键轨道之间的相互作用(n→σ*)等.     

Hydride affinity scale of various substituted arylcarbeniums in acetonitrile

Combined with the integral equation formalism polarized continuum model (IEFPCM), the hydride affinities of 96 various acylcarbenium ions in the gas phase and CH(3)CN were estimated by using the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-311++G(2df,2p)//B3LYP/6-31+G(d), and BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) methods for the first time. The results show that the combination of the BLYP/6-311++G(2df,2p)//B3LYP/6-31+G(d) method and IEFPCM could successfully predict the hydride affinities of arylcarbeniums in MeCN with a precision of about 3 kcal/mol. On the basis of the calculated results from the BLYP method, it can be found that the hydride affinity scale of the 96 arylcarbeniums in MeCN ranges from -130.76 kcal/mol for NO(2)-PhCH(+)-CN to -63.02 kcal/mol for p-(Me)(2)N-PhCH(+)-N(Me)(2), suggesting most of the arylcarbeniums are good hydride acceptors. Examination of the effect of the number of phenyl rings attached to the carbeniums on the hydride affinities shows that the increase of the hydride affinities takes place linearly with increasing number of benzene rings in the arylcarbeniums. Analyzing the effect of the substituents on the hydride affinities of arylcarbeniums indicates that electron-donating groups decrease the hydride affinities and electron-withdrawing groups show the opposite effect. The hydride affinities of arylcarbeniums are linearly dependent on the sum of the Hammett substituent parameters σ(p)(+). Inspection of the correlation of the solution-phase hydride affinities with gas-phase hydride affinities and aqueous-phase pK(R)(+) values reveals a remarkably good correspondence of ΔG(H(-)A)(R(+)) with both the gas-phase relative hydride affinities only if the α substituents X have no large electron-donating or -withdrawing properties and the pK(R)(+) values even though the media are dramatically different. The solution-phase hydride affinities also have a linear relationship with the electrophilicity parameter E, and this dependence can certainly serve as one of the most effective ways to estimate the new E values from ΔG(H(-)A)(R(+)) or vice versa. Combining the hydride affinities and the reduction potentials of the arylcarbeniums, we obtained the bond homolytic dissociation Gibbs free energy changes of the C-H bonds in the corresponding hydride adducts in acetonitrile, ΔG(HD)(RH), and found that the effects of the substituent on ΔG(HD)(RH) are very small. Simple thermodynamic analytic platforms for the three C-H cleavage modes were constructed. It is evident that the present work would be helpful in understanding the nature of the stabilities of the carbeniums and mechanisms of the hydride transfers between carbeniums and other hydride donors.     

Theoretical study of the N-H…O red-shifted and blue-shifted hydrogen bonds

Theoretical calculations are performed to study the nature of the hydrogen bonds in complexes HCHO…HNO, HCOOH…HNO, HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F. The geometric structures and vibrational frequencies of these six complexes at the MP2/6-31+G(d,p), MP2/6-311++G(d,p), B3LYP/6-31+G(d,p) and B3LYP/6-311++G(d,p) levels are calculated by standard and counterpoise-corrected methods, respectively. The results indicate that in complexes HCHO…HNO and HCOOH…HNO the N-H bond is strongly contracted and N-H…O blue-shifted hydrogen bonds are observed. While in complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F, the N-H bond is elongated and N-H…O red-shifted hydrogen bonds are found. From the natural bond orbital analysis it can be seen that the X-H bond length in the X-H…Y hydrogen bond is controlled by a balance of four main factors in the opposite directions: hyperconjugation, electron density redistribution, rehybridization and structural reorganization. Among them hyperconjugation has the effect of elongating the X-H bond, and the other three factors belong to the bond shortening effects. In complexes HCHO…HNO and HCOOH…HNO, the shortening effects dominate which lead to the blue shift of the N-H stretching frequencies. In complexes HCHO…NH3, HCOOH…NH3, HCHO…NH2F and HCOOH…NH2F where elongating effects are dominant, the N-H…O hydrogen bonds are red-shifted.     

H2CO与双卤分子间卤键的电子密度分析

运用B3LYP和MP2方法在6-311++G(d,p)基组水平上, 对H2CO-XY(XY=F2、Cl2、Br2、ClF、BrF、BrCl)卤键体系进行构型全优化, 得到了O…X—Y型卤键复合物. 结果表明, MP2/6-311++G(d,p)计算结果与实验值较吻合. 并在MP2水平下计算了分子间的相互作用能, 用完全均衡校正CP(counterpoise procedure)方法对基函数重叠误差(BSSE)进行了校正. 利用电子密度拓扑分析方法对卤键复合物的电子密度拓扑性质进行了分析研究.     

CH3SS与XO (X＝F,Cl, Br)反应过程的电子密度拓扑研究

利用MP2/6-311＋＋G(d,p)//B3LYP/6-311＋＋G(d,p)对CH3SS与XO (X＝F, Cl, Br)的反应机理进行了研究. 着重从电子密度拓扑分析角度讨论了化学键的生成和断裂. 计算结果表明单线态反应为主要反应通道, 且由于该通道的反应能垒低、放热明显, 说明CH3SS与XO在大气中比较容易进行. 电子密度拓扑分析表明, 在单线态抽氢反应通道中存在着四元环状过渡结构, 随着反应进行, 此四元环状过渡结构通过一个T-型结构变为三元环状过渡结构, 最后环状结构消失得到产物.     

Inverse Halogen Bonds Interactions Involving Br Atom in the Electronic Deficiency Systems of CH3⊕ … Br-Y (Y=H,CCH, CN,NC)

Inverse halogen bonds interactions involving Br in the electronic deficiency sys-tems of CH₃^⊕…Br-Y (Y=H, CCH, CN, NC) have been investigated by B3LYP/6-311++G(d, p) and MP2/6-311++G(d, p) methods. The calculated interaction energies with basis set super-position error correction of the four IXBs com-plexes are 218.87, 219.48, 159.18, and 143.05 kJ/mol (MP2/6-311++G(d, p)), re-spectively. The relative stabilities of the four complexes increased in the order:CH₃^⊕…BrCN3^⊕…BrNC3^⊕…BrH≈CH₃^⊕…BrCCH. Natural bond orbital theory analysis and the chemical shifts calculation of the related atoms revealed that the charges flow from Br-Y to CH₃^⊕. Here, the Br of Br-Y acts as both a halogen bond donor and an electron donor. Therefore, compared with conventional halogen bonds, the IXBs complexes formed between Br-Y and CH₃^⊕. Atoms-in-molecules theory has been used to investigate the topological properties of the critical points of the four IXBs structures which have more covalent content.