从头算研究3-氨基-1-丙醇的电离能

采用MP2/6-311++G**方法研究了3-氨基-1-丙醇势能面上能量最低的14种构象. 运用B3LYP, MP3和MP4计算了构象的相对能量,结果说明含有分子内OH…N氢键的两种构象最稳定. 电子传播子理论的P3/aug-cc-pVTZ近似方法计算了构象的电离能,与光电子能谱的实验数据符合较好. 运用自然键轨道理论解释了构象间最外价壳层轨道电离能的差异. 结合统计力学原理,计算了不同温度下的构象分布,解释了光电子能谱依赖温度的原因.     

从头计算研究CH3C(O)OSSOC(O)CH3的构型和能量

采用从头计算B3LYP,MP3和MP4方法结合Aug-cc-pVDZ基组研究了CH₃C(O)OSSOC(O)CH₃最稳定的五种构象及其阳离子构型.理论计算了五种稳定构象的转动常数和偶极矩;运用电子传播子理论P3近似方法计算稳定构象外价壳层轨道的电离能,计算结果与光电子能谱实验结果符合的较好.根据构象的相对能量以及理论模拟电离能谱和实验光电子能谱之间的比较,说明在气相光电子能谱实验中至少存在两种构象.与中性构型相比,电离后的五种阳离子构型均发生了明显的结构弛豫,尤 关键词： 构象 电离能 相对能量 光电子能谱     

L-苯丙氨酸构象和电离能的理论研究

我们对L-苯丙氨酸进行了全势能面搜索,采用B3LYP方法优化了L-苯丙氨酸的648种可能构象,最终得到了37种稳定存在的构象.分别采用B3LYP、B3PW91、M06-2X、MP2和CCSD(T)计算了L-苯丙氨酸最稳定的10种构象的相对能量,其中M06-2X和MP2方法能够给出较好的结果.对比不同的基组,说明采用aug-cc-pVDZ已经接近达到基组收敛极限.用电子传播子理论P3近似方法计算稳定构象外价壳层轨道的垂直电离能与光电子能谱实验符合的很好;根据构象的相对能量以及理论模拟与实验的光电能谱的比对,说明对气相光电子能谱至少四种构象有贡献.     

L-苯丙氨酸构象和电离能的理论研究

摘要: 我们对L-苯丙氨酸进行了全势能面搜索,采用B3LYP方法优化了L-苯丙氨酸的648种可能构象,最终得到了37种稳定存在的构象。分别采用B3LYP、B3PW91、M06-2X、MP2和CCSD(T)计算了L-苯丙氨酸最稳定的10种构象的相对能量,其中M06-2X和MP2方法能够给出较好的结果。对比不同的基组,说明采用aug-cc-pVDZ已经接近达到基组收敛极限。用电子传播子理论P3近似方法计算稳定构象外价壳层轨道的垂直电离能与光电子能谱实验符合的很好;根据构象的相对能量以及理论模拟与实验的光电能谱的比对,说明对气相光电子能谱至少四种构象有贡献.     

烯丙醇和烯丙硫醇构象的Dyson轨道和电离能

用B3LYP/aug-cc-pVTZ方法研究了烯丙醇和烯丙硫醇的构象. 采用MP3、MP4(SDQ)和CCSD(T)方法计算了它们的相对能量,这两种分子Gg0构象均最稳定. 根据理论计算的电离能模拟的光电子能谱说明在气相实验中烯丙醇和烯丙硫醇至少存在四种构象. 烯丙硫醇最外壳层和次外壳层的Dyson轨道均显示出强烈的n_S和π_C=C混合特征,这是由于S的孤对电子轨道和C=C_π轨道之间的共振效应和诱导效应引起的.     

2-氟乙醇分子外价轨道的电子动量谱学

利用不共面不对称(e,2e)谱仪测量了2-氟乙醇分子外价分子轨道的束缚能谱和电子动量分布,并用外价格林函数方法和密度泛函理论计算了分子的电离能和轨道波函数. 通过理论计算与实验结果的比较,明确标识了实验测量的电离能谱. 考虑了2-氟乙醇分子五种异构体的玻尔兹曼热力学统计权重后,理论计算的电子动量分布能够较好地解释实验测量结果.     

乙醇分子的高分辨电子动量谱学研究:单个构象异构体的轨道电子动量分布

利用高分辨电子动量谱仪测量了乙醇分子外价轨道的电离能谱,通过对一系列角度关联的电离能谱进行解谱,获得了各个电离能峰对应的分子轨道电子动量分布.利用密度泛函理论方法计算了乙醇分子两种构象异构体的轨道电子动量分布,通过与实验结果进行比较,发现实验测量的电离能为14.5和15.2 eV能峰对应的电子动量分布分别与理论计算的单个构象异构体trans 8a''和gauche 9a轨道电子动量分布符合较好.     

1-碘丙烷分子的电子动量谱学研究：自旋-轨道耦合效应与分子内轨道相互作用

利用不对称不共面电子动量谱仪,在2.5 keV碰撞能量下,采用高精度的SAC-CI方法计算了1-碘丙烷分子束缚能谱,同时采用Hartree-Fock、B3LYP/aug-cc-pVTZ(C,H)6-311G^**(I)方法计算其电子动量分布. 并对电离能峰进行了标示. 结合非相对论与相对论计算方法以及自然键轨道分析,对最外层两个轨道(碘的5p孤对)的自旋-轨道耦合效应与分子内轨道相互作用进行了比较. 两种相互作用对电子动量分布的不同影响是可观的. 实验结果与相对论计算的结果一致,表明1-碘丙烷分子内自旋-轨道耦合效应占主导.     

丁酮分子内价轨道1a"的电子动量谱学研究

用高分辨率电子动量谱仪进行丁酮分子的结合能谱和内价轨道1a"电子动量谱的实验工作,以及用HartreeFock和密度泛函理论方法对1a"轨道电子动量谱的理论研究.得到了各价轨道的电离能值以及理论计算的总能、偶极矩和1a"轨道的二维密度图.并比较了内价轨道1a"的电子动量谱的实验和理论计算结果,实验结果与理论计算符合较好.     

N9H9链状异构体的结构与稳定性的理论研究

采用密度泛函理论(DFT)B3LYP方法,在6-311++G**基组水平上对N9H9可能存在的链状异构体进行了几何优化,得到46种稳定链状异构体.应用自然键轨道理论NBO和分子中的原子理论AIM分析了这些链状异构体的成键特征和相对稳定性,G3MP2方法计算了各异构体的精确能量及在298K时的生成热△fH°(298 K),并计算了由Peter Politzer等人所介绍的相对比冲量.研究蛄果表明:各异构体中N原子孤对电子与N=N形成了P→π共轭作用是影响双键相邻的N-N键长变化的主要原因,并且对异构体的稳定性起着重要作用.所有异构体中N=N位于链端的稳定性较兰,其中B9最稳定,B6稳定性最差;C5是所有异构体中生成热最大的,也是相对比冲量最大的.     

计算分子二阶超极化率的基组选择

设计和使用了扩展的6-31G,D95, ANO,aug-cc-pVDZ基组以及Hypol基组,用各种相关能级别的方法SCF,MP2,MP3,MP4D,MP4DQ,MP4SDQ,MP4, CCSD,CCSD(T),计算了H2O,CO,HF,H2,N2等分子的二阶超极化率.分别以4-aug-cc-pVTZ基组在不同方法下的二阶超极化率结果,和CCSD(T)/4-aug-cc-pVTZ水平上的结果作为标准,求得各水平计算相对应的平均相对偏差.结果发现:MP2/aug-cc- pVDZ+是计算二阶超极化率最优的水平.综合讨论认为:可在MP2/aug-cc-pVDZ+水平上计算分子的二阶超极化率,其结果是可靠的.而且如果同时要计算分子极化率和一阶超极化率,其结果也是可靠的.     

6-硫代鸟嘌呤磺化异构体热力学稳定性的计算研究

采用量子化学密度泛函理论B3LYP方法在6-311++G(d,p)基组水平上,对6-硫代鸟嘌呤硫醇式与硫酮式两类共14种磺化异构体的热力学稳定性进行了计算研究,全自由度优化了各物种构型,进行了振动分析,获得了各磺化异构体的零点能、焓、热力学能以及吉布斯自由能等热力学数据.分析了各异构体磺化前后稳定性的变化情况,探讨了磺化异构体的稳定性问题.结果表明,不论是硫酮式还是硫醇式磺化异构体,当H原子连接到N7上时磺化异构体稳定性均比连接到N9上时较高一些.在6-硫代鸟嘌呤磺化异构体中,STG(3,7,10,10)是最稳定的异构体,STG(1,7,10,10)、STG(1,9,10,10)和STG(1,3,7,10)等异构体也可共存,其他异构体存在的可能性均较小.     

硼、氮取代二聚苯中碳原子导致层间相互作用本质变化研究

在MP2/aug-cc-pVTZ水平下优化了二聚苯及其碳被硼、氮原子取代的几何构型;计算了这些二聚体的CCSD(T)/CBS相互作用能;并用SAPT2+/aug-cc-pVDZ分析了相互作用能成分.探索了硼、氮取代二聚苯中碳原子所导致层间相互作用本质变化.结果表明:稳定构型而言,从苯(C_6H_6)、1-氮-2-硼杂苯(BNC_4H_6)、1,3-二氮-2,4-二硼杂苯(B_2N_2C_2H_6)二聚体的平行移位(PD)构型向无机苯(B_3N_3H_6)二聚体的夹心(S)构型转变,其中C_6H_6和B_3N_3H_6二聚体稳定构型的堆叠型式分别与之相对应的石墨烯(GE)和六方氮化硼(h-BN)2D层间材料的堆叠型式相一致.硼、氮原子取代二聚苯中的碳原子后使其相互作用能增大,其中BNC_4H_6和B_2N_2C_2H_6二聚体的相互作用能增大较为明显.所研究体系二聚体稳定构型均以色散能为主导、静电能次之、诱导能相对较小.硼、氮取代二聚苯中碳原子后其静电能对总吸引能的贡献明显增大.     

密度泛函方法研究气相胞嘧啶的互变异构化

分别采用7种基组、3种理论方法对胞嘧啶异构体Cyt1的结构进行优化,通过与Cyt1的实验结果进行比较,选取了适合研究胞嘧啶分子的B3LYP/6311+G方法.用该方法对胞嘧啶分子的8种异构体构型进行了充分优化,研究了其中能量较低的6种胞嘧啶异构体的互变异构化过程.对于得到的所有优化构型都进行了频率分析.对于基态构型,所有的频率都是正的;对于过渡态构型,只有一个虚频.同时,做了详尽的内禀反应坐标计算,以保证所得到的过渡态连接相应的始末异构体.所有给出的能量都已做了零点能校正.理论研究结果可以对已有的实验结果给予合理解释.     

Theoretical study on the molecular tautomerism of the 3-hydroxy-pyridin-4-one system

3-hydroxy-pyridin-4-one is a parent molecule for the family of hydroxypyridinones that are known in coordination chemistry as efficient metal ions chelators. In this work, relative stabilities of some possible tautomers were investigated using several quantum chemical methods: CBS (complete basis set methods), Gn, DFT (density functional theory), Hartree–Fock and MP2. Performed calculations show that the system under consideration exists as a mixture of two tautomers with comparable energies. Among them, the hydroxypyridinone structure of the studied molecular system seems to be a bit more stable than the o-dihydroxypyridine one, by a few kJ/mol only. Aromaticity and intra-molecular hydrogen bonding are the main effects influencing the stability of the studied tautomeric structures. Consequently, aromatic effects were calculated using several indices of aromaticity: HOMA (harmonic oscillator model of aromaticity), NICS (nucleus independent chemical shift), H, PDI (para delocalisation index), MCI (multi-centre index) and ASE (aromatic stabilisation energy). The strength of possible intra-molecular hydrogen bonds (H-bonds) was determined by means of the AIM (atoms-in-molecules) method and by calculating enthalpies for theoretical reactions that do or do not involve H-bonds. The AIM method was employed to understand how variations in atomic energies influence the stability of different tautomeric structures.     

Mass spectrometric and ab initio study of the interaction between 9-methylguanine and amino acid amide group

The temperature dependent field ionization mass spectrometry method combined with ab initio calculations was used to determine the interaction energies and the structures of 9-methylguanine-acrylamide dimers. Acrylamide mimics the side chain amide group of the natural amino acids asparagine and glutamine. The experimental enthalpy of the dimer formation derived from the van't Hoff plot is ?59.5 ± 3.8 kJ mol^?1. The value is higher than interaction energies between acrylamide and other nucleic acid bases which were determined to be ?57.0 for 1-methylcytosine, ?52.0 for 9-methyladenine, and ?40.6 kJ mol^?1 for 1-methyl-uracil. In total, eight hydrogen bonded dimers formed by the three lowest energy 9-methylguanine tautomers and acrylamide were found in the quantum chemical calculations performed at the DFT/B3LYP/6-31++G?? and MP2/6-31++G?? levels of theory. The relative stability and the interaction energies of the dimers were calculated accounting for the basis set superposition error and the zero-point vibrational energy correction. The lowest energy dimer found in the calculations is formed by acrylamide (Ac) with the keto tautomer of 9-methylguanine (Gk). It is stabilized by two intermolecular H bonds, C6=O(Gk) · · · H—N(Ac) and Nl—H(Gk) · · ·O(Ac), and it is more stable than the second lowest energy dimer by ≈ 25 kJ mol^?1. The calculated interaction energies of the lowest energy 9-methylguanine-acrylamide dimer are ?65.0 kJ mol^?1 and ?67.7 kJ mol^?1 at the MP2 and DFT levels of theory, respectively. The experimental enthalpy of the dimer formation is in good agreement with both the calculated interaction energies of the GkAc dimer and much higher than the interaction energies calculated for all other 9-methylguanine-acrylamide dimers. This proved that only one dimer was present in the experimental samples. To verify whether acrylamide is a good model of the amino acid-amide group, we performed direct calculations of the 9-methylguanine-glutamine dimers at the same levels of theory as used for the complexes involving acrylamide. The interaction energies found for the lowest energy 9-methylguanine-glutamine dimer are ?65.1 kJ mon^?1 (MP2/6-31++G??) and ?66.2 kJ mol^?1 (DFT/B3LYP/6-31++G??) and these values are very close (within 0.5 kJ mol^?1) to the interaction energies obtained for the 9-methylguanine-acrylamide dimers.     

Theoretical kinetics study of thymine tautomerism and interaction of Na+ with its tautomers

The current work is a theoretical study of the tautomerism of thymine in the gas phase. Eighteen structures were found in the isomerisation reaction of thymine, some of which are reported for the first time. Thirty hydrogen transfer reactions were carried out. In 24 of the reactions, the hydrogen abstractions N―H→O, N―H→C and C―H→O were considered. The potential energy surface for all trajectories was determined for 18 tautomers and 40 transition states. The RRKM-TST model was used to calculate the rate constants of the reactions to examine their kinetics. Nonlinear least-squares fitting was used to calculate the rate constants expressions. The interaction of sodium ion and tautomers in the gas phase was also investigated. Three types of interaction of metal cations with thymine were found. In the first, metal cations interact with a lone pair of nitrogen or oxygen tautomers. The second type is the interaction of metal cations with two nitrogen and oxygen of tautomers. The last type is the interaction of metal cations and the electron density of the π-system of thymine in which the metal ion is perpendicular to the ring of tautomers. The stability ranking of the thymine tautomers and their complexes was also determined.     

Enhancement effect of lithium bonding on the strength of pnicogen bonds: XH2P···NCLi···NCY as a working model (X = F,Cl; Y = H,F, Cl,CN)

MP2 calculations with aug-cc-pVDZ basis set were used to analyse intermolecular interactions in XH₂P···NCLi···NCY triads (X = F, Cl; Y = H, F, Cl, CN) which are connected via pnicogen bond and lithium bond. To understand the properties of the systems better, the corresponding dyads are also studied. Molecular geometries and interaction energies of dyads, and triads are investigated at the MP2/aug-cc-pVDZ computational level. Particular attention is paid to parameters such as cooperative energies and many-body interaction energies. All studied complexes, with the simultaneous presence of a lithium bond and a pnicogen bond, show cooperativity with energy values ranging between ?4.73 and ?8.88 kJ mol^?1. A linear correlation was found between the interaction energies and magnitude of the product of most positive and negative electrostatic potentials. According to energy decomposition analysis, it is revealed that the electrostatic interactions are the major source of the attraction in the title complexes.     

Cooperative and diminutive interplay between the sodium bonding with hydrogen and dihydrogen bondings in ternary complexes of NaC3N with HMgH and HCN (HNC)

Ternary complexes of NaC₃N with HMgH and HCN (HNC) are connected by sodium, hydrogen and dihydrogen bonds. Molecular geometries and interaction energies of dyads and triads are investigated at the MøllerPlesset perturbation theory of the second order/aug-cc-pVDZ computational level. Particular attention is paid to parameters, such as cooperative energies and many-body interaction energies. Triads with the HMgH molecule located at the end of the chain show an energetic cooperativity ranging between ?2.13 and ?10.53 kJ mol^?1. When the HMgH molecule is located in the middle, the obtained cluster is diminutive with an energetic effect with values 4.39 and 6.77 kJ mol^?1. The electronic properties of the complexes are analysed using parameters derived from the atoms in molecules methodology. Based on the energy decomposition analysis, it can be seen that the stabilities of the complexes are predicted to be attributable mainly to electrostatic effects.