2,3,7,8-四氯苯并二英分子结构与热力学性质的理论研究

用密度泛函理论(DFT)和从头算(ab initio)方法,在B3LYP/6-31G、 B3LYP/6-31G*、 B3LYP/6-311G*和MP2/6-31G*水平上全优化计算了2,3,7,8-四氯苯并二英(2,3,7,8-TCDD)的几何构型、电子结构和振动频率,并用校正后的频率计算了298～1500 K的标准热力学函数,同时用半经验的PM3 SCF-MO进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

A DFT study for the regioselective 1,3-dipolar cycloadditions of nitrile N-oxides toward alkynylboronates

The mechanism for the 1,3-dipolar cycloaddition of benzonitrile oxide toward ethynyl and propynylboronate has been studied by using density functional theory (DFT) at B3LYP/6-31G* level. These cycloadditions are concerted [3+2] processes. The presence of the two oxygens on the boronic ester precludes the participation of the boron atom on [3+3] processes. The two regioisomeric channels associated to the formation of the isoxazoles bearing the boronic ester unit on the 4- or 5-positions have been characterized. The B3LYP/6-31G* activation parameters are in acceptable agreement with the experiments, allowing to explain the factors controlling these regioselective cycloadditions.     

Synthesis,Characterization and Theoretical Studies on a Thiocarbamide Derivative Containing Schiff Base Group

A thiocarbamide derivative containing Schiff base groups,1,5-bis[4-(dimethylamino)benzylidene]thiocarbonohydrazide,has been synthesized and characterized by elemental analysis,IR,1H NMR,UV and X-ray single-crystal diffraction.Density function theory(DFT) calculations at the B3LYP/6-31G* and PBE0/6-31G* levels for optimized geometries and electronic transition spectra have been performed.Comparative studies show that both B3LYP/6-31G* and PBE0/6-31G* methods can well reproduce the molecular structure,and the latter is more reliable than the former to simulate electronic spectra.NPA calculational results at the B3LYP/6-31G* level indicate the title compound to be a potential multidentate ligand to the metallic ions.Based on the vibrational analysis,thermodynamic properties at different temperatures have been obtained.     

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键扭转振动基频的预测值。     

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.     

NMR spectra, GIAO and charge density calculations of five-membered aromatic heterocycles

The B3LYP/6-31+G(d) molecular geometry optimized structures of 17 five-membered heterocycles were employed together with the gauge including atomic orbitals (GIAO) density functional theory (DFT) method at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p) and B3LYP/6-311+G(2d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants. The method of geometry optimization for pyrrole (1), N-methylpyrrole (2) and thiophene (7) using the larger 6-311++G(d,p) basis sets at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,p) and B3LYP/cc-pVTZ levels of theory gave little difference between calculated and experimental values of coupling constants. In general, the (1)H and 13C chemical shifts for all compounds are in good agreement with theoretical calculations using the smaller 6-31 basis set. The values of nJHH(n=3, 4, 5) and rmnJ(CH)(n=1, 2, 3, 4) were predicted well using the larger 6-31+G(d,p) and 6-311++G(d,p) basis sets and at the B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p), B3LYP/6-31+G(2d,2p) levels of theory. The computed atomic charges [Mülliken; Natural Bond Orbital Analysis (NBO); Merz-Kollman (MK); CHELP and CHELPG] for the B3LYP/6-311++G(d,p) geometry optimized structures of 1-17 were used to explore correlations with the experimental proton and carbon chemical shifts.     

Conformational study of the structure of 12-crown-4-alkali metal cation complexes

A conformational search was performed for the 12-crown-4 (12c4)-alkali metal cation complexes using two different methods, one of them is the CONFLEX method, whereby eight conformations were predicted. Computations were performed for the eight predicted conformations at the HF/6-31+G*, MP2/6-31+G*//HF/6-31+G*, B3LYP/6-31+G*, MP2/6-31+G*//B3LYP/6-31+G*, and MP2/6-31+G* levels. The calculated energies predict a C4 conformation for the 12c4-Na+, -K+, -Rb+, and -Cs+ complexes and a C(s) conformation for the 12c4-Li+ complex to be the lowest energy conformations. For most of the conformations considered, the relative energies, with respect to the C4 conformation, at the MP2/6-31+G*//B3LYP/6-31+G* are overestimated, compared to those at the MP2/6-31+G* level, the highest level of theory considerd in this report, by 0.2 kcal/mol. Larger relative energy differences are attributed to larger differences between the B3LYP and MP2 optimized geomtries. Binding enthalpies (BEs) were calculated at the above-mentioned levels for the eight conformations. The agreement between the calculated and experimental BEs is discussed.     

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

用多种密度泛函理论(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)对此分子的振动基频进行了理论归属。     

Rearrangement of azirine intermediates to nitriles: Theoretical study of cleavage of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to cyanoketene–formaldimine complex

A detailed investigation of the reaction path for the thermal rearrangement of 3,4-dihydro-1aH-azirine[2,3-c]pyrrol-2-one to yield a cyanoketene–formaldimine complex is carried out at the MP2/6-31G* and B3LYP/6-31G* levels of theory. The ring opening of the five-membered pyrrolinone ring and the formation of the nitrile group takes place in a concerted manner, presenting a significant strain energy release and allowing for an electronic stabilization by coarctate conjugation of the transition structure (TS). These two factors make possible a moderate energy barrier. Although the structural features B3LYP/6-31G* theoretical levels, it is found that the MP2 energy barrier (28.8) CCSD(T)/6-31G*//MP2/6-31G* value (17.1 kcal/mol). The complex electronic rearrangement can be rationalized using the theory of coarctate transition structures developed by Herges as the evolution of an azirine structure without referring to a hypothetical vinyl nitrene intermediate. © 1998 John Wiley & Sons, Inc. J Comput Chem 19: 912–922, 1998     

2-丁基-四氢噻吩亚砜13C-NMR的理论研究

在RHF/6-31G和B3LYP/6-31G水平上对顺式(Cis-)与反式(Trans-)2-丁基-四氢噻吩亚砜(BTHTO)进行几何优化,应用规范不变原子轨道法(GIAO)在6-31G、6-31+G、6-31++G和6-31+G(2d,p)水平上计算了Cis-和Trans-BTHTO的13C-NMR,对13C-NMR谱进行了归属。结果表明,BTHTO噻吩五元环的稳定构象呈半椅式,Cis-和Trans-BTHTO中与硫原子直接碳原子13C-NMR的显著差异主要是由于空间构型不同引起分子的静电势场对相应碳原子的屏蔽作用不同所致。     

线型碳链PC2nP的结构和电子光谱的密度泛函理论研究

应用密度泛函理论,在B3LYP/6-31G**和B3LYP/6-311G**水平上优化得到了线型簇合物PC2nP(n=1-10)的基态平衡几何构型,计算了它们的谐振动频率.在基态平衡构型下,利用含时密度泛函理论,计算得到了簇合物PC2nP(n=1-10)的垂直激发能和相应的振子强度,导出了激发能与体系大小n的解析关系式.     

硝酸甲与硝酸乙酯的结构、振动频率和热力学性质的密度函数 理论(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)与实验结果较吻合。     

Iminopropadienones R–NCCCO and carbon suboxide, C3O2. Theoretical and experimental C NMR spectra

The ¹³C NMR data of five iminopropadienones R–NCCCO as well as carbon suboxide, C₃O₂, have been examined theoretically and experimentally. The best theoretical results were obtained using the GIAO/B3LYP/6-31+G**//MP2/6-31G* level of theory, which reproduces the chemical shifts of the iminopropadienone substituents extremely well while underestimating those of the cumulenic carbons by 5–10 ppm. The computationally faster GIAO/HF/6-31+G**//B3LYP/6-31G* level is also adequate.     

Theoretical study of light-metal tetrahydroborates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroborates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the structural parameters, the energies of different decomposition pathways, barriers to internal rotation of BH₄ groups, normal mode frequencies, magnetic shielding constants, and spin density distribution (in radicals) have been analyzed in various related series of these compounds. The results obtained by ab initio methods and at the DFT level are compared. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** adequately reproduces the results obtained at the higher level of theory CCSD(T)/6-311++G*s*//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.     

Electronic, optical, and charge transfer properties of donor–bridge–acceptor hydrazone sensitizers

The ground state geometries have been computed by using density functional theory (DFT) at B3LYP/6-31G*, B3LYP/6-31G**, and PCM-B3LYP/6-31G* level of theories. The highest occupied molecular orbitals (HOMOs) are delocalized on whole of the molecule and the lowest unoccupied molecular orbitals (LUMOs) are localized on the tricarbonitrile. The lowest HOMO and LUMO energies have been observed for Dye1 while highest for Dye4. The LUMO energies of Dye1–Dye4 are above the conduction band of TiO₂ and HOMOs are below the redox couple. The absorption spectra have been computed in solvent (methanol) and without solvent by using time-dependant DFT at TD-B3LYP/6-31G*, TD-B3LYP/6-31G**, and PCM-TD-B3LYP/6-31G* level of theories. The calculated maximum absorption wavelengths of the spectra in methanol are in good agreement with experimental evidences. The maximum absorption wavelengths of new designed sensitizers are red shifted compared to parent molecule. The electronic coupling constant and electron injection have been computed by first principle investigations. The improved electronic coupling constant and electron injection revealed that new modeled systems would be efficient sensitizers.     

羰基硼化合物(BCO)n(n==1～12)的理论研究

运用HF/3-21G方法和密度泛函理论(DFT)的B3LYP/6-31G*方法, 对羰基硼化合物(BCO)n(n=1~12)的各种可能结构进行了优化, 对在B3LYP/6-31G*水平上得到的几何构型、电子态、结合能、振动频率、核独立化学位移(NICS)、能量二次差分和热力学性质进行了理论研究, 得到了(BCO)n(n=1~12)结构的稳定性信息. 十二种基态结构都是端配位(μ1-CO), (BCO)n(n=1~3, 5, 6) 的基态是线型或平面结构, (BCO)n(n=4, 7~12)的基态是笼状结构; B—C平均键能呈现奇偶交替现象, 偶数的结构比奇数稳定; 能量二次差分得到同样的结论；羰基的振动频率与实验值非常吻合; 热力学性质的研究对实验具有重要的指导意义.     

Neutral homoaromaticity in some heterocyclic systems

[reaction: see text] Neutral homoaromaticity has been evaluated in heterocyclic systems related to the bicyclo[3.2.1]octane skeleton with replacement of CH(2) at C-2 in bicyclo[3.2.1]octa-3,6-diene with X = BH, AlH, Be, Mg, O, S, PH, NH (12); replacement of CH at C-3 in bicyclo[3.2.1]octa-3,6-dien-2-yl anion with PH, S, NH, O (13); and replacement at C-2 and C-3 with N and O (14). Stabilization energies (SE) are evaluated using density functional theory and homodesmotic equations at the B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d) level for series 12-14. Stabilization energies are compared with diamagnetic susceptibility exaltations, Lambda, CSGT-B3LYP/6-31G(d)//B3LYP/6-31G(d), and nucleus-independent chemical shifts (NICS), GIAO-B3LYP/6-311+G(2d,p)//B3LYP/6-31G(d). Analysis of frontier orbitals and geometries, B3LYP/6-31G(d)//B3LYP/6-31G(d), and proton affinities of 2-azabicyclo[3.2.1]octa-3,6-diene, pyrrole, and related model systems, B3LYP/6-311+G(2d,2p)//B3LYP/6-31G(d), provide complementary evidence supporting the division of the substrates evaluated into antihomoaromatic (12, X = BH, AlH, and Be), nonhomoaromatic (12, X = O, S, NH, PH), and homoaromatic (13, X = S, PH, NH, O and 14 X = ON), with 12 (X = Mg) appearing as transitional between anti- and nonhomoaromatic.     

DFT法研究3-羟基丙烯醛的双键旋转异构反应机理

利用密度泛函理论(DFT)分别在B3LYP/6-31G**和B3LYP/6-311＋＋G**的计算水平上优化了基态3-羟基丙烯醛分子在双键旋转异构反应过程中的平衡态以及过渡态的几何构型，分析了反应过程中键参数的变化，计算了该反应的内禀反应坐标(IRC)，发现在重排反应途径上存在一个四元环骨架的中间体.通过振动分析对平衡态和过渡态进行了确认，并得到了零点能.计算结果表明，基态3-羟基丙烯醛分子的双键旋转异构反应经过两步完成，第一步反应位垒稍高，第二步反应位垒较低，存在着发生重排反应的可能性.     

Diketoacid HIV-1 integrase inhibitors: An ab initio study

The stable tautomeric forms of two representative arene-substituted diketoacid HIV-1 integrase inhibitors, 5-ClTEP and L-731,988, were investigated by B3LYP with 6-31G*, 6-31G(d,p), and 6-31+G(d,p) basis sets. Optimization with MP2/6-31G* was also performed for 5-ClTEP. The solvation effect was considered using a conductor-like screening model. With the density functional theory method, the trans diketo conformations are more stable than the cis conformers. The difference is 14 kJ mol(-1) for 5-ClTEP and 33 kJ mol(-1) for L-731,988. Two trans diketo structures were obtained. The difference between these two trans diketo structures is less than 4 kJ mol(-1) calculated at the B3LYP/6-311+G(3df,2p) level. Two enol forms prevail over the diketo tautomers and are calculated to have the same free energy. Because there is no barrier observed between these two enol forms, they can interchange easily such that a delocalized transition state is suggested to be the observed form. Contradictory to the results of the MP2 method that predicts a preference for the trans diketo forms, the B3LYP method predicts a preference for the enol tautomers, which is in agreement with the experimental results.     

Theoretical study of light-metal tetrahydroaluminates

Ab initio calculations of the structural, energetic, vibrational, and magnetic characteristics of the lowest-lying structures for isolated molecules and ions of light-metal tetrahydroaluminates (Li, Na, Be, Mg, and Al) have been performed by the perturbation theory (MP2), quadratic configuration interaction (QCISD(T)), coupled cluster (CCSD(T)), and density functional theory (B3LYP) methods using the 6-31G*, 6-31G**, 6-311+G**, and 6-311++G** basis sets. The trends in the behavior of the molecular characteristics have been analyzed in various related series of these compounds. The results are compared with the data on analogous light-metal tetrahydroborates calculated at the same levels of theory. The differences in structure and stability between analogous hydroborate and alanate complexes are examined. The economical approximation B3LYP/6-311++G**//B3LYP/6-311+G** has been shown to adequately reproduce the results obtained at the higher level of theory CCSD(T)/6-311++G**//MP2/6-31G* even though it requires considerably shorter CPU times and smaller amounts of memory.