丁二酰亚胺的结构、振动频率和热力学性质计算

用密度泛函理论（ＤＦＴ）和从头计算（ａｂ ｉｎｉｔｉｏ）方法,在Ｂ３ＬＹＰ／６－３１Ｇ、Ｂ３ＬＹＰ／６－３１１Ｇ和ＭＰ２／６－３１１Ｇ水平上,全优化计算了了二酰亚胺的分子几何构型和电子结构．进行了简正振动频率分析并用校正后的频率计算了２００－６００Ｋ温度范围的标准热力学函数,对计算结果进行了比较和讨论．     

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进行了同样的计算,计算结果与实验值及文献值较好地吻合.     

DFT study on mechanism of the classical Biginelli reaction

The condensation of benzaldehyde, urea, and ethyl acetoacetate according to the procedure described by Biginelli was investigated at the B3LYP/6-31G（d）, B3LYP/6-31 ＋G（d,p）, and B3LYP/6-311＋G（3df,2p）//B3LYP/6-31＋G（d,p） levels to explore the reaction mechanism. According to the mechanism proposed by Kappe, structures of five intermediates were optimized and four transition states were found. The calculation results proved that the mechanism proposed by Kappe is right.     

一些密度泛函方法预测电子亲和势

选取了杂化泛函B3LYP, B3PW91, O3LYP, PBE0, 以及与之相对应的GGA泛函BLYP, BPW91, OLYP和PBE, 还选取了能更好地兼顾强相互作用和弱相互作用的X3LYP泛函和在预测NMR的化学位移有较好表现的OPBE泛函, 以及两种meta-GGA泛函VSXC和TPSS, 共12种泛函, 详细地考察了这些泛函在预测EA方面的准确性.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

丁烯和丁烷自由基阳离子的分子结构和超精细结构的 理论研究

用密度函数理论B3LYP方法和6-31G(d,p),6-311G(d,p)及6-311+G(d,p)基组，分别对1-C4H^+~8,2-C4H^+~8和C4H^+~10进行了构型优化和频率分析计算，预言1-C4H^+~8具有非平面构型，与以往报道的从头算和密度函数理论计算结果不同。在各自由基阳离子的B3LYP构型上，进行了B3LYP、MP2及MRSDCI方法的超精细偶合常数计算，得到了比以往更好的结果，特别是MP2/B3LYP计算值是至今与实验值符合得最好的理论计算结果。     

Cyclic vinyl p-tolyl sulfilimines as chiral dienophiles: theoretical study on the stereoselectivity, Lewis acid catalysis, and solvent effects in their Diels-Alder reactions

The theoretical study reported in the present work deals with chiral cyclic vinyl sulfilimines and their reactivity as dienophiles in [4 + 2] cycloaddition reactions, using B3LYP/6-31G(d)//AM1 and B3LYP/6-31G(d)//B3LYP/6-31G(d) model chemistries. Consideration of Lewis acid catalysis, illustrated by BF(3), decreases the activation energies of the cycloaddition process while the charge transfer from the diene to the sulfilimine is augmented. The [4 + 2] cycloaddition reactions of sulfilimines with both furan and cyclopentadiene occur in the gas phase with endo stereoselectivity, which is more pronounced with the latter diene. Endo-exo energy differences in the gas phase with the B3LYP/6-31+G(d)//B3LYP/6-31+G(d), B3LYP/6-31G(d)//B3LYP/6-31G(d), and B3LYP/6-31G(d)//AM1 model chemistries are almost the same. Solvent effects are responsible for the inversion of the stereoselectivity in the reactions of sulfilimines with furan because of the great difference in the dipole moments in endo and exo approaches.     

D-苯丙氨酸与Cu＋(1S0, 3d10)气相反应理论研究

用量子化学密度泛函(DFT)方法研究D-苯丙氨酸与一价基态金属阳离子Cu^＋在气相中反应的机理. 在B3LYP/6-31G*水平上, 优化了反应包含的4个反应通道的反应物、中间体、过渡态和产物的几何构型, 并采用B3LYP/DZVP, B3LYP/[6-311＋G**(C,H,O)＋Lanl2dz(Cu)], B3LYP/6-311＋G**, MP2/6-311＋G** 等方法对各驻点进行了单点能计算. 通过对计算结果的分析, 获得了其单重态反应势能面的一般轮廓、各驻点几何构型优化参数, 明确了其反应机理.     

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

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

Theoretical investigation on the GaCl3-catalyzed ring-closing metathesis reaction of N-2,3-butadienyl-2-propynyl-1-amine: three-membered ring versus four-membered ring mechanism

The gallium chloride (GaCl(3))-catalyzed ring-closing metathesis reaction mechanism of N-2,3-butadienyl-2-propynyl-1-amine has been studied at the Becke three-parameter hybrid functional combined with Lee-Yang-Parr correlation functional (B3LYP)/6-31G(d), B3LYP/6-31+G(d,p), B3LYP/6-311++G(d,p)//B3LYP/ 6-31G(d) and the second-order M?ller-Plesset perturbation (MP2)/6-311++G(d,p)//B3LYP/6-31+G(d,p) levels. It was found that the final metathesis product can be yielded via a three-membered or four-membered ring mechanism. The three-membered ring pathway is favorable due to its low energy barrier at the rate determining step. The whole reaction is stepwise and strongly exothermic.     

On the flexibility of beta-peptides

The full conformational space was explored for an achiral and two chiral beta-peptide models: namely For-beta-Ala-NH2, For-beta-Abu-NH2, and For-beta-Aib-NH2. Stability and conformational properties of all three model systems were computed at different levels of theory: RHF/3-21G, B3LYP/6-311++G(d,p)//RHF/3-21G, B3LYP/6-311++G(d,p), MP2//B3LYP/6-311++G(d,p), CCSD//B3LYP/6-311++G(d,p), and CCSD(T)//B3LYP/6-311++G(d,p). In addition, ab initio E = E(phi, micro, psi) potential energy hypersurfaces of all three models were determined, and their topologies were analyzed to determine the inherent flexibility properties of these beta-peptide models. Fewer points were found and assigned than expected on the basis of Multidimensional Conformational Analysis (MDCA). Furthermore, it has been demonstrated, that the four-dimensional surface, E = E(phi, mu, psi), can be reduced into a three-dimensional one: E = E[phi, f(phi), psi]. This reduction of dimensionality of freedom of motion suggests that beta-peptides are less flexible than one would have thought. This agrees with experimental data published on the conformational properties of peptides composed of beta-amino acid residues.     

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.     

2-溴丙酸气相热消除反应的机理

用密度泛函（DFT）方法,在B3LYP/6-31G**水平上对2-溴丙酸气相消除反应机理进行了研究.计算表明,反应主要是通过半极化五元环结构过渡态进行的,羧基上的氢原子协助溴原子离去,羧基氧原子帮助稳定过渡态.在B3LYP/6-311＋＋G(3df,3pd)水平上对B3LYP/6-31G**优化的几何构型进行了单点能计算,计算所得反应的速度控制步骤的活化能为189.461 kJ•mol－1,偏离实验值((180.3±3.4) kJ•mol－1)5.08％.     

DFT Calculations of the Ionization Potential and Electron Affinity of Alaninamide

Adiabatic and vertical ionization potentials (IPs) and valence electron affinities (EAs) of alaninamide in gas phase have been determined using density functional theory (BLYP,B3LYP,B3P86) methods with 6-311++G(d,p) basis set,respectively. IPs and EAs of alaninamide in solutions have been calculated at the B3LYP/6-311++G(d,p) level. Five possible conformers of alaninamide and their charged states have been optimized employing density functional theory B3LYP method with 6-311++(d,p) basis set,respectively.     

Structural elucidation of nitro-substituted five-membered aromatic heterocycles utilizing GIAO DFT calculations

The GIAO (Gauge Including Atomic Orbitals) DFT (Density Functional Theory) method is applied at the B3LYP/6-31+G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311++G(d,p)//B3LYP/6-31+G(d), B3LYP/6-311+G (2d,p)//B3LYP/6-31+G(d) and B3LYP/6-311++G(d,p)//B3LYP/6-311++G(d,p) levels of theory for the calculation of proton and carbon chemicals shifts and coupling constants for 25 nitro-substituted five-membered heterocycles. Difference (1D NOE) spectra in combination with long-range gHMBC experiments were used as tools for the structural elucidation of nitro-substituted five-membered heterocycles. The assigned NMR data (chemical shifts and coupling constants) for all compounds were found to be in good agreement with theoretical calculations using the GIAO DFT method. The magnitudes of one-bond (1JCH) and long-range (nJCH, n>1) coupling constants were utilized for unambiguous differentiation between regioisomers of nitro-substituted five-membered heterocycles.     

Proton affinity of proline and modified prolines using the kinetic method: role of the conformation investigated by ab initio calculations

The proton affinities of proline, cis-3-methylproline and cis-3-ethylproline have been measured by the kinetic method using an ion trap instrument; the values obtained are 936, 940.5, and 943 kJ mol(-1), respectively. The experimental values are consistent with those obtained by high-level ab initio calculations (B3LYP/6-31+G*//B3LYP/6-31G* and B3P86/6-31+G*//B3LYP/6-31G*). Several conformations of neutral and protonated proline were considered, in particular the endo and exo ring structure and the position of the carboxyl group. These results show the importance of the position of the hydrogen atom of the carboxyl group in determining the most stable protonated proline structure.     

反式-双(二甲基苯膦)铂配合物二阶非线性光学性质的密度泛函理论研究

采用密度泛函理论B3LYP方法,在6-31G(d)的水平上对两种反式-双(二甲基苯膦)铂配合物的几何构型进行优化,在获得稳定构型基础上,利用TD-B3LYP方法得到体系的UV-Vis吸收光谱,并用有限场(FF/B3LYP)方法探讨体系的二阶非线性光学性质(NLO).结果表明,此类铂配合物具有较大的二阶极化率,以及在可见光范围内透明等优点,是具有很好应用前景的非线性光学材料.     

具有D~6~h对称性的C~3~6及其衍生物电子结构的密度泛函研究

运用密度泛函方法,比较不同水平的基组(HF/6-311+G^*;B3LYP/6-31G^*;B3LYP/6-311+G^*)对具有D~6~h对称性的C~3~6分子进行构型优化的结果,并分析其几何结构、电子结构、稳定性等性质;采用基组B3LYP/6-31G^*对H@C~3~6,Li@C~3~6,Na@C~3~6,K@C~3~6分子进行构型全优化,分析了不同内嵌原子对其几何结构、电子结构、稳定性等性质的影响;首次在B3LYP/6-311+G^*水平上,对C~3~6H~6,C~3~6H~1~2,X@C~3~6(X=H,Li,Na,K)几何构型及电子结构进行研究并得到其稳定性规律。     

Theoretical Study on the Structures and Isomerization of Germylenoid GeH2LiCl

Using geometrical optimization and DFT method at the B3LYP/6-311++G (3df,3pd) level, four equilibrium geometries and one transition state of GeH2LiCl were identified, and the structures at the MP2/6-311++G(3df,3pd) level were calculated simultaneously. We also studied the solvent effects on the structures of Germylenoid GeH2LiCl at the B3LYP/6-311++G (3df,3pd) level. The two more stable forms are suggested to be the p-complex and three-membered ring. The vibrational frequencies and infrared intensities were computed at the B3LYP/6-311++G (3df,3pd) level.     

A comparison of the behavior of functional/basis set combinations for hydrogen-bonding in the water dimer with emphasis on basis set superposition error

We evaluate the performance of ten functionals (B3LYP, M05, M05-2X, M06, M06-2X, B2PLYP, B2PLYPD, X3LYP, B97D, and MPWB1K) in combination with 16 basis sets ranging in complexity from 6-31G(d) to aug-cc-pV5Z for the calculation of the H-bonded water dimer with the goal of defining which combinations of functionals and basis sets provide a combination of economy and accuracy for H-bonded systems. We have compared the results to the best non-density functional theory (non-DFT) molecular orbital (MO) calculations and to experimental results. Several of the smaller basis sets lead to qualitatively incorrect geometries when optimized on a normal potential energy surface (PES). This problem disappears when the optimization is performed on a counterpoise (CP) corrected PES. The calculated interaction energies (ΔEs) with the largest basis sets vary from -4.42 (B97D) to -5.19 (B2PLYPD) kcal/mol for the different functionals. Small basis sets generally predict stronger interactions than the large ones. We found that, because of error compensation, the smaller basis sets gave the best results (in comparison to experimental and high-level non-DFT MO calculations) when combined with a functional that predicts a weak interaction with the largest basis set. As many applications are complex systems and require economical calculations, we suggest the following functional/basis set combinations in order of increasing complexity and cost: (1) D95(d,p) with B3LYP, B97D, M06, or MPWB1k; (2) 6-311G(d,p) with B3LYP; (3) D95++(d,p) with B3LYP, B97D, or MPWB1K; (4) 6-311++G(d,p) with B3LYP or B97D; and (5) aug-cc-pVDZ with M05-2X, M06-2X, or X3LYP.