可乐定分子构象异构和互变异构的理论研究

采用HF/6-31G*方法, 对氨基型和亚胺型可乐定中性分子以及质子化的可乐定分子进行了构象分析. 之后采用B3LYP/6-31G**方法对处于势能面上的能量极小构象异构体进行全自由度几何优化和频率分析. Onsager反应场溶剂模型用于水相的计算. 结果表明, 在气相和水相中可乐定中性分子主要以亚胺型存在, 这同已有实验结果一致. 进一步, 寻找构象异构化过渡态和氨基型亚胺型互变异构化过渡态, 探讨质子化和水溶剂化效应对异构化过程的几何结构和能量的影响. 为了考察氯苯的共轭效应对可乐定互变异构体稳定性的影响, 在B3LYP/6-31G**水平上, 研究了2-氨基-2-咪唑啉的互变异构化反应机理.     

聚偏二氯乙烯构象分布的理论研究

在MP2/6-311++G**水平下, 对2,2,4,4-四氯戊烷与2,2,4,4,6,6-六氯庚烷旋转异构体构象进行几何优化和能量计算. 结果表明, 对于2,2,4,4-四氯戊烷, 采用gauche-gauche排列的旋转异构体的能量较低; 2,2,4,4,6,6-六氯庚烷旋转异构体中, 采用trans-gauche-trans-gauche排列的构象能量较低. 反之, 完全采用trans-trans排列的旋转异构体构象能量较高, 不稳定. 通过比较模型分子不同旋转异构体构象的能量差值可以得到一级和二级特征的相互作用能差, 进而计算统计权重参数. 在此基础上, 应用计算得到的模型分子的几何构型与统计权重参数, 分别构建针对—CH2—和—CCl2—中心的聚偏二氯乙烯的6态旋转异构态模型. 通过旋转异构态模型可以计算聚偏二氯乙烯分子中各种构象的分布.     

5-氟尿嘧啶和5-氯尿嘧啶及其互变异构体的理论计算研究

采用HF/3-21G方法, 对6种气相和水相中可能存在的5-氟尿嘧啶(和5-氯尿嘧啶)互变异构体进行了构象分析.采用B3LYP/6-311＋G**方法对处于优势构象时的各互变异构体进行了几何全优化, 并计算出它们的总能量、焓、熵、吉布斯自由能. Onsager反应场溶剂模型用于水相的计算. 计算结果表明, 5-氟尿嘧啶和5-氯尿嘧啶在气相中和水相中主要以双酮形式存在. 5-氟尿嘧啶和5-氯尿嘧啶的熵效应小, 对互变异构平衡没有显著的影响, 而焓变对互变异构产生了主要的影响. 讨论了水溶剂化作用对异构体的能量、电荷分布和偶极矩的影响. 溶剂化自由能与异构体的气相偶极矩存在相关性. 另外, 详细地将5-氟尿嘧啶和5-氯尿嘧啶与尿嘧啶进行了对比, 获得三者最稳定异构体间电子结构异同的有用信息.     

2,4-二氟戊烷的构象与能量的理论研究

在MP2/6-311++G**和G3水平下,对2,4-二氟戊烷异构体的构象进行几何优化和能量计算.氟的强电负效应使构型发生偏转,同时导致了相隔超过3个键的原子间的电子相互作用.     

5-氟胞嘧啶互变异构的密度泛函理论计算

采用BH-HLYP/6-311＋G**方法对10种气相和水相中可能存在的5-氟胞嘧啶互变异构体进行了几何全优化, 并计算出它们的总能量和吉布斯自由能. Onsager反应场溶剂模型用于水相的计算. 计算结果表明, 5-氟胞嘧啶在气相中主要以烯醇式-氨基式形式存在, 在水相中主要以酮式-氨基式形式存在. 溶剂化自由能与异构体的气相偶极矩存在相关性．进一步求得互变异构化以及构象异构化和顺反异构化的过渡态, 探讨异构化过程中的几何结构和能量的变化.     

氟比洛芬构象异构体的理论计算

用Hypercbem软件中构象搜寻模块,对氟比洛芬进行构象搜索,寻找分子低能构象．用B3LYP／6-31G（d）法优化计算22个低能构象,PCM溶剂模型用于水相计算,获得低能构象的优化几何结构、分子总能量和标准吉布斯自由能．结果表明,在气相和水相时,Z-型构象异构体稳定性最高．     

2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷构象稳定性的理论研究

采用密度泛函方法在B3LYP/6-31＋G**水平上研究了2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷分子(D4G)的构象. 分别研究在气相中的孤立分子和一水合物异构体的相对稳定性和异构体之间的相互转变过程, 分析了水分子的参与对D4G异构体的相对稳定性和几何结构参数以及自然电荷的影响. 结果表明, 孤立的D4G分子在气相中存在8种稳定构象, 其中构象d4g-2是所有构象中最稳定的, 气相中D4G主要以d4g-2存在. 气相中各构象的相对稳定性为: d4g-2＞d4g-1＞d4g-5＞d4g-3＞d4g-6＞d4g-4＞d4g-8＞d4g-7. 计算得到的各构象键长和键角数据与实验值接近. 一个水分子的加入对D4G分子的构型参数有所影响, 基本不改变D4G分子各构象的稳定性顺序, 但构象转变的能垒有所提高. 氢键在分子构象中发挥了重要作用.     

2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷构象稳定性的理论研究

采用密度泛函方法在B3LYP/6-31＋G**水平上研究了2',3'-二脱氧-2',3'-二去氢鸟嘌呤核苷分子(D4G)的构象. 分别研究在气相中的孤立分子和一水合物异构体的相对稳定性和异构体之间的相互转变过程, 分析了水分子的参与对D4G异构体的相对稳定性和几何结构参数以及自然电荷的影响. 结果表明, 孤立的D4G分子在气相中存在8种稳定构象, 其中构象d4g-2是所有构象中最稳定的, 气相中D4G主要以d4g-2存在. 气相中各构象的相对稳定性为: d4g-2＞d4g-1＞d4g-5＞d4g-3＞d4g-6＞d4g-4＞d4g-8＞d4g-7. 计算得到的各构象键长和键角数据与实验值接近. 一个水分子的加入对D4G分子的构型参数有所影响, 基本不改变D4G分子各构象的稳定性顺序, 但构象转变的能垒有所提高. 氢键在分子构象中发挥了重要作用.     

环型硫分子的从头算及其应用

通过对环型硫分子S6～S20同素异构体的实验结构分析,用从头算(ab initio)RHF/6-311G*(包括BLYP/6-311G* 和MP4/6-311G*)方法进行基态几何结构优化,获得了与实验结构相吻合的新的理论稳定构型.其中S18的两种变体S18(α)和S18(β)的能量相近,都是相互稳定的分子构型.通过计算,从理论上推测S20有D4点群的物相.还就对称性、偶极矩和红外振动光谱与分子构象的相互关系,以及环分子结构与化学活性之间的关系进行了讨论,获得了具有实验意义的结论.     

6-硫代黄嘌呤互变异构体的密度泛函理论计算

在密度泛函B3LYP/6-311G**水平下,对14种气相和水相中可能存在的6-硫代黄嘌呤异构体进行了几何构型的全自由度优化,并计算出它们的总能量、焓、熵、吉布斯自由能。Onsager反应场溶剂模型用于水相的计算.计算结果表明,6-硫代黄嘌呤在气相中和水相中主要以硫酮的形式存在.在气相和水相中,硫酮-N7(H)均比硫酮-N9(H)更稳定.计算结果同已有实验结果一致.6-硫代黄嘌呤异构化的熵效应小,对互变异构平衡几乎没有显著的影响,而焓变对互变异构产生了主要的影响.较详细地讨论了水溶剂化作用对异构体的能量、几何结构、电荷分布和偶极矩的影响.     

电场作用下分子导线的理论研究

摘要利用从头计算法分别在HF/6-31G, HF/6-31G*, HF/6-31G**, HF/6-31+G, HF/6-31++G, HF/6-31+G*, HF/6-31+G**, HF/D95+*, B3LYP/6-31G*和B3LYP/6-31+G*水平上计算了5个单体的聚乙炔分子导线, 从几何构型、 SCF能量和分子轨道能级三个方面讨论了外电场对分子导线的影响, 给出了聚乙炔分子导线性质与外电场变化的定量关系.     

Mycoepoxydiene核磁共振谱的理论研究

用密度泛涵理论(Density Functional Theory,DFT)的B3PW 91方法在6-311G**基组下优化了标题化合物mycoepoxyd iene的几何构型,同时计算了红外光谱.并用规范不变原子轨道(gauge independent atom ic orb ital)G IAO/B3PW 91以及G IAO/HF方法分别在6-311G**,6-311++G**等基组进行了核磁共振谱研究,计算结果与实验结果吻合很好.     

BrONO2光解反应的理论研究

采用密度泛函方法,研究了大气臭氧层主要破坏物BrONO2的光解反应机理,在UB3LYP/6-311++G**水平上优化了反应物、产物、中间体和过渡态的几何构型,并在UQCISD(T)/6-311++G**水平上计算了单点能量,为了确证过渡态的真实性,在UB3LYP/6-311++G**水平上进行了内禀反应坐标(IRC)计算和频率分析.研究结果表明,BrONO2的光解反应有两条反应通道,其中生成BrO+NO2的反应活化能较小(14.89 kJ·mol-1),较易发生.     

Structure and conformational stability of CH2CHCH2X (X=F, Cl and Br) molecules—post Hartree–Fock and density functional theory methods

The molecular structure and conformational stability of CH₂CHCH₂X (X=F, Cl and Br) molecules were studied using ab initio and density functional theory (DFT) methods. The molecular geometries of 3-fluoropropene were optimized employing BLYP and B3LYP levels of theory of DFT method implementing 6-311+G(d,p) basis set. The MP2/6-31G^*, BLYP and B3LYP levels of theory of ab initio and DFT methods were used to optimize the 3-chloropropene and 3-bromopropene molecules. The structural and physical parameters of the molecules are discussed with the available experimental values. The rotational potential energy surface of the above molecules were obtained at MP2/6-31G^* and B3LYP/6-311+G(d,p) levels of theory. The Fourier decomposition of the rotational potentials were analyzed. The HF/6-31G^* and MP2/6-31G^* levels of theory have predicted the cis conformer as the minimum energy structure for 3-fluoropropene, which is in agreement with the experimental values, whereas the BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory reverses the order of conformation. The ΔE values calculated for 3-chloropropene at MP2/6-31G^*, BLYP/6-311+G(d,p) and B3LYP/6-311+G(d,p) levels of theory show that the gauche form is more stable than the cis form, which is in agreement with the experimental value. The same levels of theory have also predicted that the gauche form is stable than cis for 3-bromopropene molecule. The maximum hardness principle has been able to predict the stable conformer of 3-fluoropropene at HF/6-31G^* level of theory, but the same level of theory reverses the conformational stability of 3-chloropropene and 3-bromopropene molecules and MP2/6-31G^* level of theory predicted the stable conformer correctly.     

CH_3自由基和O（~3P）反应机理的量子化学研究

用分子轨道从头计算MP2（full）方法和密度泛函理论（DFT）中的B3LYP方法 研究了CH_3自由基和三线态O原子反应的微观机理，优化得到了反应途径上的反应 物、过渡态、中间体和产物的几何构型，通过振动分析对过渡态和中间体构型进行 了确认，在G3不平上计算了能量，同时用经典过渡态理论对该反应的绝对速率常数 进行了理论计算。研究结果表明：CH_3自由基与O（~3P）反应有四条不同的放热反 应通道，主反应通道为IM1→TS1→CH_2O + H，同时反应可彻底裂解生成CO, H_2 及H。     

1,4-二氢-5H-(二硝基亚甲基)-四唑(DNMT)的合成、晶体结构和热行为研究

利用1-氨基-1-肼基-2,2-二硝基乙烯(AHDNE)和亚硝酸钾在酸性水溶液中合成出了高能富氮化合物1,4-二氢- 5H-(二硝基亚甲基)-四唑(DNMT), 并在水溶液中培养出DNMT•2H2O单晶. 该晶体属正交晶系, 空间群为Pnma, 晶胞参数为: a＝10.392(2) Å, b＝15.809(4) Å, c＝5.0640(11) Å, V＝832.0(3) Å3, Dc＝1.629 g•cm－3, μ＝0.163 mm－1, F(000)＝432, Z＝4, R1＝0.0311, wR2＝0.0885. 运用Gaussian 03程序, 在6-311＋＋G**基组水平上, 用HF和B3LYP两种方法对DNMT分子进行了几何全优化和相应的电荷、轨道能量分析. 理论计算和热分析结果表明DNMT呈现较差的热稳定性.     

Conformational stability, the spectroscopic (FT-IR and UV), first order hyperpolarizability, NBO analysis, HOMO and LUMO analysis of 6,8-diphenylimidazo[1,2-α]pyrazine molecule by ab initio HF and density functional methods

The conformational analysis of 6,8-diphenylimidazo[1,2-α]pyrazine molecule (abbreviated as 68DIP) was performed by using B3LYP/6-31G(d) level of theory to find the most stable form. Two staggered stable conformers were observed on the torsional potential energy surface. The equilibrium geometry, bonding features and vibrational frequencies of 68DIP have been investigated by using the DFT (B3LYP) and HF methods for the lowest energy conformer. The first order hyperpolarizability (β(total)) of this molecular system and related properties (β, μ, <α> and Δα) are calculated using HF/6-311++G(d,p) and B3LYP/6-311++G(d,p) methods based on the finite-field approach. Stability of the molecule arising from hyperconjugative interactions, charge delocalization and C-H?N intramolecular hydrogen-bond-like weak interaction has been analyzed using natural bond orbital (NBO) analysis by using B3LYP/6-311++G(d,p) method. The results show that electron density (ED) in the σ* and π* antibonding orbitals and second order delocalization energies E((2)) confirm the occurrence of intramolecular charge transfer (ICT) within the molecule. UV-vis spectrum of the compound was recorded and electronic properties, such as HOMO, LUMO energies, excitation energies and wavelength were performed by TD-DFT/B3LYP, CIS and TD-HF methods by using 6-311++G(d,p) basis set. Finally, the calculation results were applied to simulated infrared spectra of the title compound which show good agreement with observed spectra.     

An Ab Initio Molecular Orbital Study of Valence Bond Isomers of Silabenzene and Phosphabenzene

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single point total energy calculation are reported for silabenzene ( 7 ), phosphabenzene ( 8 ) and 16 valence bond isomers of silabenzene and phosphabenzene ( 9-24 ). The calculated energy difference (19.78 kcal mol m 1 ) between silabenzene and the most stable valence bond isomer of silabenzene (1-silabenzvalene, 9 ) is much smaller than the difference (73.60 kcal mol m 1 ) between benzene and benzvalene ( 2 ). The energy difference between phosphabenzene and the most stable valence bond isomer of phosphabenzene (1-phosphabenzvalene, 17 ) is calculated to be 43.29 kcal mol m 1 .     

An Ab Initio Molecular Orbital Study of Isophosphinines

Ab initio molecular orbital calculation at HF/6-31G*, HF/6-31G**, HF/6-311G**, HF/6-311++G**, RMP2-FC/6-31G*, and B3LYP/6-31G* levels of theory for geometry optimization and MP4(SDQ)/6-31G* for a single-point total energy calculation are reported for phosphinine and 13 isophosphinines 7-19 . Isomers 7-11 with an allenic system are calculated to be 8-18 kcal mol m 1 more stable than structures 12-17 with an acetylenic moiety. The calculated energy difference (66.19 kcal mol m 1 ) between phosphinine and the most stable isophosphinine (1-phospha-1,2,4-cyclohexatriene, 10 ) is smaller than the difference (78.96 kcal mol m 1 ) between benzene and the most stable isobenzene (cyclohexa-1,2,4-triene, 2 ). The isophosphinines 18 and 19 , with a butatriene moiety, are calculated to be the least stable isomers.