1, 2-二硒方酸气相酸性和芳香性的理论研究

在RHF/6-311G^*^*水平优化得到1,2-二硒方酸(3,4-二羟基-3-环丁烯-1,2-二硒酮)三种平面构象异构体的平衡几何构型。进一步用MP2(full)/6-311G^*//RHF/6-311G^*^*方法计算三种异构体的单点能量,发现ZZ型异构体是能量最低构象,且ZZ和ZE型能量非常接近。用优化的最稳定构象ZZ型异构体在RHF/6-311G^*^*//RHF/6-311G^*^*,RHF/6-311+G^*^*//RHF/6-311+G^*^*,MP2(full)/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其气相酸性[ΔGⅲ~(~2~9~8~K~)]和同键反应芳香性稳定化能(HASE)。用基团加和法(groupincrementapproach)在RHF/6-311+G^*^*//RHF/6-311+G^*^*和B3LYP/6-311+G^*^*//B3LYP/6-311+G^*^*水平计算其磁化率增量(Λ)。计算结果指出标题化合物的键长发生了平均化,同键反应芳香性稳定化能和磁化率增量均为负值,表明它具有芳香性,实现了标题化合物芳香性的几何、能量和磁性的判定。     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(DFT)研究

用密度函数理论(DFT)的BLYP和B3LYP方法,取6-31G,6-31G^*,6-31G^*^*,6-311G,6-311G^*和6-311G^*^*六种基组,对硝酸甲酯和硝酸乙酯的几何构型和红外振动频率进行了计算研究.结果表明,B3LYP方法在采用极化基组(6-31G^*,6-31G^*^*,6-311G^*和6-311G^*^*)时计算得到的结果均较好,适用于硝酸酯类化合物的研究.而BLYP方法无论采用何种基组均不适用;运用校正后的B3LYP/6-31G^*频率(校正因子0.975)计算得到的热力学性质(C^o~p,H^o和S^o)与实验结果较吻合。     

吲哚分子振动光谱的密度泛函理论研究

用多种密度泛函理论(DFT)方法(BLYP/6-31G^*^*,B3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)对吲哚分子的平衡几何构型进行了优化。在优化构型的基础上计算了吲哚分子的谐力场、振动基频和红外光谱强度。计算得到的振动频率与实验值比较平均偏差对四种计算方法(BLYP/6-31G^*^*,P3LYP/6-31G^*^*,B3PW91/6-31G^*^*和SVWN/6-31G^*^*)分别为16.3,40.5,45.1和26.4cm^-^1。BLYP/6-31G^*^*理论力场被用于吲哚分子的简正坐标分析计算中。根据振动率的势能分布(PEDs)对此分子的振动基频进行了理论归属。     

吡啶-BH~3相互作用复合物的理论研究

对吡啶-BH~3复合物分别用MP2/6-31+G^*和B3LYP/6-31+G^*进行理论计算以预测该复合物的构型及解离能，得到四种构型，在MP2优化构型基础上作CCSD/6-31+G^*单点能量计算以验证MP2与B3LYP结果的可靠性，然后用B3LYP作振动频率分析，计算了各构型的垂直电离势，最后用更大基组作单点能量计算和自然键轨道(NBO)分析。结果表明，N-B直接相连的构型最稳定，其解离能为141.50kJ/mol,MP2和B3LYP对N-H接近的构型结果相关较大，另外两种构型稳定性介于二者之间，解离能分别为15.18kJ/mol,14.06kJ/mol(MP2/6-31+G^*)。     

HeLi^n^+(n=0, 1)簇合物的稳定性和垂直激发态光谱

本文用ab initio研究了簇合物HeLi^n^+(n=0, 1)的几何构型和成键性质。在MP2(FULL)/6-31G**, 水平优化所得LeLi^+的平衡键长为0.2062nm, 与实验值0.205nm十分吻合。比较了HeLi^+(X^1∑^+和a^3∑^+), HeLi(X^2∑^+和a^4II)以及HLi(X^1∑^+)的稳定性, 计算了HeLi^+基态的相关能, 势能曲线和垂直激发态光谱。计算采用了6-31G**, 6-311G**,6-311G(2df, 2pd), 6-311G(3df, 2pd)和6-311+G(3df, 2pd)基组; 采用的方法包括MP2(FULL), MP4, MCSCF, MRSDCI, CCD和ST4CCD。计算表明, 同价HeLi^n^+中激发态的离解能均远比基态的大, 其中HeLi^+(a^3∑^+)的离解能最高(60.49kj/mol),说明激发态是稳定束缚态。HeLi^+基态比等电子体HLi分子基态的稳定性小得多。HeLi^+由A^1∑^+到B^1II的垂直跃迁(3σ→1π)振子强度较大而垂直跃迁能较小。     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。     

电子转移反应O2+O2·-→O2·-+O2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

电子转移反应O~2+O~2^.^-→O~2^.^-+O~2的从头算研究

利用abinitio方法,在UHF,UMP2及不同基组3-21G,6-31G^*,6-311+G^*和UMP2(full)/6-311+G^*水平上,研究了O~2/O~2^.^-自交换电子转移反应。优化了电子转移前后反应物和产物的结构,研究了体系能量的变化,计算了自交换电子转移反应的内重组能。对UHF方法和UMP2方法的计算结果进行了比较,并与实验结果进行了对照。结果表明UHF方法由于没有考虑组态相互作用,计算结果存在较大偏差,UMP2(full)/6-311+G^*水平上的计算结果与实验值吻合较好。在UMP2(full)/6-311+G^*水平上计算了气相自交换电子转移反应速率常数。在优化了电子转移复合物结构的基础上考虑了溶剂效应的影响,计算了水溶液中的溶剂重组能。研究结果表明O~2/O~2^.^-体系电子转移反应的活化能主要来源于溶剂重组能的贡献。最后计算了该反应在水溶液中的反应速率常数。理论计算结果与实验值吻合得很好。     

二氨基二硝基乙烯结构和性质的理论研究

对三种二氨基二硝基乙烯同分异构体进行了HF/6-31G^*^*水平、DFT-B3LYP/6-31G^*^*水平的几何全优化以及MP2/6-31G^*^*//HF-6-31G^*^*水平的总能量计算。结果表明,1,1-二氨基-2,2-二硝基乙烯(Ⅰ)总能量比顺式(Ⅱ)和反式(Ⅲ)1,2-二氨基-1,2-二硝基乙烯的总能量低,即热力学稳定性次序为Ⅰ>Ⅲ>Ⅱ。分子的共轭性和分子内氢键的强度次序为Ⅰ≈Ⅲ>Ⅱ,前沿轨道能级差次序为Ⅰ>Ⅱ>Ⅲ,也均表明(Ⅰ)最稳定。此外还计算研究了标题物的红外光谱;化合物Ⅰ的理论计算与实验值良好相符。在此基础上计算研究了标题物的热力学性质。     

A computational study of intramolecular proton transfer in gaseous protonated glycine

The minimum energy paths for intramolecular proton transfer between the amino nitrogen and carbonyl oxygen atoms in gaseous protonated glycine were estimated at the Hartree-Fock (HF) and second-order M?ller-Plesset Perturbation (MP2) levels of theory. Potential energy profiles and their associated reactant, transition state, and product species calculated at the MP2/6-31G* level were shown to differ significantly from those obtained at the HF/6-31G* level. Effects of electron correlation and basis functions on the calculated geometries and energies of relevant species were examined at the HF, MP2, MP4, CCSD, and B3LYP levels using the 6-31G*, 6-31G**, 6-31+G**, 6-311+G**, 6-31+G(2d,2p), 6-311+G(3df,2p), cc-pVDZ, aug-cc-pVDZ, and cc-pVTZ basis sets. The HF and MP2 optimized levels with the 6-31G*, 6-31G**, 6-31+G**, and 6-311+G** bases were used to calculate the thermodynamic and kinetic properties of the proton transfer reaction at 298.15 K and 1 atm, which include enthalpy, entropy, Gibbs free energy, equilibrium constant, potential energy barriers, tunneling transmission coefficients, and rate constants. Results indicate that the proton in a carbonyl O-protonated glycine undergoes a rapid migration to the amino nitrogen atom, while the reverse process is extremely unfavorable. The objective of this work is to develop practical theoretical procedures for studying proton transfer reactions in amino acids and peptides and to assemble physical data from these model calculations for future references.     

An Ab Initio Study of the Complexes Formed by the Positive Acetylene Ion with the Hydrogen and Nitrogen Molecules and the Argon Atom

The complexes formed by the positive acetylene ion with the hydrogen molecule, the nitrogen molecule, and the argon atom are investigated with ab initio calculations using the 6-311G** and the 6-31+G(2df,2pd) basis sets. MP2/6-311G** energies and optimum geometries are obtained, as well as single-point MP3, MP4, and QCISD(T) energies with the MP2/6-311G** optimized geometries. Single-point calculations are performed with the 6-31+G(2df,2pd) basis set at MP2/6-311G** optimized geometries.     

O~3＋NH→HNO＋O~2反应机理的量子化学研究

用密度泛函(DFT)的B3LYP方法(6－31++G^*^*)研究了臭氧与NH自由基反应的微观机理,优化得到反应途径上的反应物,过渡态,中间体和产物的构型,通过振动分析对过渡态和中间体进行了确认。对单点用QCISD(T)/6－31++G^*^*方法计算能量,同时进行零点能校正。研究结果表明：NH与O~3反应有两条不同的反应通道,且均表现为亲电反应特征,两条不同的反应均为强放热反应。     

Absolute configurations, predominant conformations, and tautomeric structures of enantiomeric tert-butylphenylphosphinothioic acid.

Vibrational absorption and circular dichroism spectra of dextrorotatory, levorotatory, and racemic mixture of tert-butylphenylphosphinothioic acid have been measured in CCl(4) solutions in the 2000-900 cm(-1) region. The conformations for both tautomeric structures of (S)-tert-butylphenylphosphinothioic acid are investigated using the B3LYP functional with the 6-31G* basis set. For the most stable conformation, the absorption and VCD spectra are predicted ab initio using the B3LYP functional with 6-31G*, 6-311G(2d, 2p), 6-31+G, and 6-311G(3df, 3pd) basis sets. A different functional, B3PW91, was also used with the 6-31G* basis set. The predicted spectra are compared to the experimental spectra. The comparison indicates that (-)-tert-butylphenylphosphinothioic acid is of the (S)-configuration and exists in only one tautomeric structure and one conformation in CCl(4) solution.     

Conformational analysis of carbonic anhydrase inhibitors using ab initio molecular orbital methods. 1. Rotational isomerism in methane sulfonamide anion, CH3-SO2-NH−

Ab initio calculations have been performed on methane sulfonamide anion. Geometries have been optimized using Hartree-Fock basis sets up to 6-31+G*, and single-point calculations employing those Hartree-Fock geometries have been performed at levels up to MP2/6-311++G**. In addition, geometry optimizations for the 0°, 90°, 150°, and 180° conformers have been carried out at the MP2/6-31G*, MP2/6-31+G*, and MP2/6-311++G** levels. Vibrational frequencies have been calculated using the HF/4-31G*, MP2/6-31G*, and MP2/6-31+G* geometries. All calculations at or above the 4-31G* level agree that H—N—S—C ˜90° is the global minimum. The H—N—S—C = 180° conformer is clearly higher in energy although the relative energy of this conformer varies from 0.36 to 1.03 kcal/mol for the post-HF calculations depending on basis set. The H—N—S—C = 180° conformer appears to be a very shallow local minimum. However, the potential energy surface is quite flat in this region, and the highest-level calculations, including MP2 optimizations and vibrational frequency analysis, are ambiguous on this point. The conformer with an H—N—S—C torsion of 0° is a transition state with a relative energy ˜8 kcal/mol. Received: 3 December 1996 / Accepted: 2 January 1997     

TDDFT investigation on nucleic acid bases: comparison with experiments and standard approach

A comprehensive theoretical study of electronic transitions of canonical nucleic acid bases, namely guanine, adenine, cytosine, uracil, and thymine, was performed. Ground state geometries were optimized at the MP2/6-311G(d,p) level. The nature of respective potential energy surfaces was determined using the harmonic vibrational frequency analysis. The MP2 optimized geometries were used to compute electronic vertical singlet transition energies at the time-dependent density functional theory (TDDFT) level using the B3LYP functional. The 6-311++G(d,p), 6-311(2+,2+)G(d,p), 6-311(3+,3+)G(df,pd), and 6-311(5+,5+)G(df,pd) basis sets were used for the transition energy calculations. Computed transition energies were found in good agreement with the corresponding experimental data. However, in higher transitions, the Rydberg contaminations were also obtained. The existence of pisigma* type Rydberg transition was found near the lowest singlet pipi* state of all bases, which may be responsible for the ultrafast deactivation process in nucleic acid bases.