Ab initio prediction of the vibrational spectrum of the hydrogen bonded complex O2N...HONO2

The changes in the structural parameters and vibrational characteristics (vibrational frequencies, infrared intensities and Raman activities) arising from the hydrogen bonding between NO(2) and HONO(2) have been studied employing ab initio 6-31G(d, p)/UHF and 6-31+G(d, p)/UHF, and B3LYP/6-31G(d, p) calculations. The charge rearrangement upon hydrogen bonding have been, estimated using the Mulliken population analyses. It was established that the complexation between NO(2) and HONO(2) leads to changes in the structural parameters and the vibrational characteristics of the monomers. The most sensitive to the hydrogen bond formation are the vibrational characteristics of the normal modes of the monomer bonds participating in the hydrogen bonding. The predicted shifts in the vibrational frequencies by ab initio and B3LYP/6-31G(d, p) calculations are in very good agreement with the experimentally observed, which is an evidence for the reliance of the studied structure.     

3-羟基丙烯醛单键旋转异构光反应机理的从头算研究

应用从头算分子轨道法分别在RHF/6-31G**和UHF/6-31G**水平上对3腔┤┓肿拥幕?S0)和三重激发态(T1)单键旋转异构反应机理进行了研究,优化出反应物和产物在S0态和T1态的4种平衡态和过渡态的几何构型,通过振动分析得到的虚频和计算的内禀反应坐标对过渡态进行了确认,并得到了零点能,根据基态和激发态反应途径分析了光反应机理。计算结果表明,基态和激发态的3腔┤┓肿拥ゼ旃狗从ξ焕萁系停嘏湃菀追⑸永砺凵涎橹ち斯夥从κ笛橹屑钢忠旃固宓拇嬖凇?     

氰基苯阴离子与CO2间的内球电子转移

用半经验AM1方法以及从头算方法在3-21G和6－31G*基组水平上研究了从氰基苯阴离子到CO2的电子转移过程．结果表明，对于先驱物(precursor),三种计算方法得出的给体、受体间的距离分别为0.2728nm(AM1)、0.2479nm（UHF/3-21G）和0.2769nm(UHF/6-31G*)．在这样短的距离内给体的HOMO与受体的LUMO轨道具有相当程度的重叠，应产生较强的相互作用，说明此反应是内球电子转移反应，从而解释了此类体系的电子转移反应不符合Marcus理论的原因．计算给出先驱物的束缚能为0.19eV(AM1)和0.26eV(6-31G*)     

Structures and energies of D-galactose and galabiose conformers as calculated by ab initio and semiempirical methods

Optimized geometries and total energies of some conformers of alpha- and beta-D-galactose have been calculated using the RHF/6-31G* ab initio method. Vibrational frequencies were computed at the 6-31G* level for the conformers that favor internal hydrogen bonding, in order to evaluate their enthalpies, entropies, Gibbs free energies, and then their structural stabilities. The semiempirical AM1, PM3, MNDO methods have also been performed on the conformers GG, GT, and TG of alpha- and beta-D-galactose. In order to test the reliability of each semiempirical method, the obtained structures and energies from the AM1, PM3, and MNDO methods have been compared to those achieved using the RHF/6-31G* ab initio method. The MNDO method has not been investigated further, because of the large deviation in the structural parameters compared with those obtained by the ab initio method for the galactose. The semiempirical method that has yielded the best results is AM1, and it has been chosen to perform structural and energy calculations on the galabiose molecule (the disaccharides constituted by two galactose units alpha 1,4 linked). The goal of such calculations is to draw the energy surface maps for this disaccharide. To realize each map, 144 different possible conformations resulting from the rotations of the two torsional angles psi and phi of the glycosidic linkage are considered. In each calculation, at each increment of psi and phi, using a step of 30 degrees from 0 to 330 degrees, the energy optimization is employed. In this article, we report also calculations concerning the galabiose molecule using different ab initio levels such as RHF/6-31G*, RHF/6-31G**, and B3Lyp/6-31G*.     

离子簇合物HeN_2~+的电子光谱的从头算研究

本文报道了离子簇合物HeN2＋的电子光谱的从头算结果。在MCSCF／6－31G（d，p）水平上，对其基态进行了几何优化，用该构型．在SDC1／6－31（d，p）水平上计算了基态的总能量。用SDCI方法计算得到HeN2＋从基态到选择激发态的垂直跃迁能、振子强度、跃迁频率、辐射寿命以及Einstain系数。该结果可以较好的验证maier的实验。     

On the variation of magnetic susceptibility of a molecular crystal with temperature: the 2,4,6-triphenylverdazyl system

Magnetic susceptibilities of spin-1/2 systems of orthorhombic and higher crystal symmetries have been numerically investigated while taking possible anisotropy in the coupling constants along different crystal axes into account. The work relies on the magnon-based theory of ferromagnetic (FM) and antiferromagnetic (AFM) crystal systems of types FFF, AFF, AAF, and AAA [J. Chem. Phys. 111, 9009 (1999)]. The AAF crystal, in particular, shows interesting changes in the temperature dependence of magnetic susceptibility when the ferromagnetic exchange coupling constant is varied. We especially show that the susceptibility anomalies of molecular crystals fit naturally within the framework of the extended magnon-theoretical formalism, and do not necessarily imply a FM --> AFM or a reverse phase transition. A real system, molecular crystal of 2,4,6-triphenylverdazyl (2,4,6-TPV), has been investigated here. It was previously interpreted as an AAF system from observed susceptibility data [Tomiyoshi et al., Phys. Rev. B 49, 16031 (1994)]. The trend of the temperature dependence of magnetic susceptibility studied in the present work also indicates that the crystal belongs to the AAF category with a less prominent FM exchange coupling constant. To reinforce our conclusions, we have adopted a two-pronged strategy. First, the geometry of the 2,4,6-TPV monomer has been optimized here by ab initio unrestricted Hartree-Fock (UHF) calculations using the STO-3G basis set. The optimized geometry is almost planar. A subsequent calculation has been carried out with the phenyl rings twisted out of the plane of the nitrogen atoms. The STO-3G optimized geometry, and the same geometry except for the twisted phenyl rings, have been used to perform ab initio coupled-cluster (UCCSD-T) calculations with the same basis, and UHF as well as density-functional (UB3LYP) calculations using the 6-31G basis set. The calculated data can easily rationalize the twists while the species remains in crystal. The magnetic category of the crystal has been unambiguously confirmed as AFA from ab initio UHF and UB3LYP calculations of the total energy in different spin states of dimers and trimers along the crystal axes. The computed energy values, however, fail to yield accurate estimates of the exchange coupling constants Ja, Jb, and Jc, because the latter are on the order of 1kBK corresponding to energy differences on the order of 10(-6) hartree between different spin states. In the second approach, the observed features of the susceptibility minimum and maximum have been used to determine the best values of the exchange coupling constants from the theoretical formulas for an anisotropic AFA crystal. The AFM (Ja and Jc) and FM (Jb) exchange coupling constants and the Neel temperature (TN) found from this analysis correspond to Ja + Jc = -1.05 kBK, Jb = 1.35 kBK, and TN = 1.75 K. The calculated J values significantly differ from those estimated from a linear Heisenberg chain model, but generate a susceptibility versus temperature graph that mimics the experimental plot.     

自由基CH2CO+的平衡构型和超精细耦合常数的从头计算

本文用UHF从头计算法处理π自由基CH_2CO~+, 确定基态并研究其平衡构型, 电子结构和自旋性质。在自旋性质的计算中, 本文使用了由Amos和Snyder~(1)提出的方法。通过从UHF波函数中投影去掉最主要的污染自旋成分, 得到近似的自旋纯态波函数。由这种波函数得到的关于自旋性质的计算结果比UHF结果大为改善, 计算的hfc常数与ESR实验结果符合得相当好。     

HCO＋HCN→HCHO＋CN反应机理研究

用ａｂｉｎｉｔｉｏ方法在ＵＨＦ／６－３１Ｇ基组上研究了ＨＣＯ　＋ＨＣＮ→ＨＣＨＯ＋ＣＮ　反应的反应机理。经寻找反应的过渡态和计算内禀反应坐标（ＩＲＣ）表明，该反应经一过渡态而形成产物，具有较宽的位垒，并属于后位垒反应。用ＵＭＰ２∥６－３１Ｇ方法计算的活化位垒为１４９．６０２ｋＪ／ｍｏｌ，与实验值（１４３．００１ｋＪ／ｍｏｌ）相当一致。     

Triple hydrogen bonding in a circular arrangement: ab initio, DFT and first-principles MD studies of tris-hydroxyaryl enamines

First-principles Car-Parrinello molecular dynamics, ab initio (MP2) and density functional schemes have been used to explore the tautomeric equilibrium in three tris(amino(R)methylidene)cyclohexane-1,3,5-triones (R?=?hydrogen, methyl or phenyl group). The dynamic nature of the cyclic hydrogen bonding has been studied by the first-principles MD method. The comparison of the results obtained by aforesaid methods has been accomplished on the basis of calculations of structural and spectroscopic characteristics of the compounds. The conformational analysis of the studied compounds has been carried out at the MP2/6-31+G(d,p) and B3LYP/6-31+G(d,p) levels of theory. The influence of steric and electronic effects on the cyclic hydrogen bonding has been analysed. The extent of the proton delocalization has been modified by the substituents according to the sequence: hydrogen?<?phenyl?<?methyl. This fact is verified by the spectroscopic and structural data as well as the energy potential curve. A prevalence of the keto-enamine tautomeric form has been observed in the static ab initio and DFT models, and confirmed by the first-principles MD.     

气态烯丙醇光异构化反应的理论研究

本文采用自洽场分子轨道UHF/6-31G从头计算法, 辅以能量梯度法研究气相烯丙醇光异构化反应的机理。全部优化了T1态势能面上反应物、过渡态、中间体和产物的几何构型。基于Fukui提出的内禀反应坐标理论(IRC)计算这一体系的反应途径, 并针对各驻点进行MP2/6-31G的相关能校正, 得到该反应在激发态进行并为一经历双自由基中间体的分步反应的结论。支持了实验工作者提出的机理。     

硝酸甲酯分子间相互作用的DFT和ab initio比较

用密度泛函理论(DFT)和从头算(ab initio)方法,分别在B3LYP/6 31G和HF/6 31G水平上求得硝酸甲酯三种二聚体的全优化几何构型和电子结构,并用6 311G和6 311＋＋G基组进行总能量计算.对HF/6 31G计算结果进行MP4SDTQ电子相关校正.在各基组下均进行基组叠加误差(BSSE)和零点能(ZPE)校正求得结合能.对6 31G优化构型作振动分析并基于统计热力学求得200～600 K温度下单体和二聚体的热力学性质.详细比较两种方法的相应计算结果,发现DFT求得的分子间距离较短,分子内键长较长,所得结合能均小于相应ab initio计算值.     

双氮桥联1,10-菲咯啉大环化合物分子构型的从头算研究

应用规范不变原子轨道法(GIAO)在RHF/6-31G**和B3LYP/6-31G**水平上计算了质子化双氮桥联1,10-菲咯啉大环化合物(H4HAPP2+)C2h和C2h构型的1HNMR,并用TDDFT法计算了H4HAPP2+电子光谱.结果表明,B3LYP/6-31G*优化的C2h构型为较优构型,经谐振频率验证无虚频,C2h构型是H4HAPP2+合理的对称性构型.     

Vibrational coupling between the nitro groups in 1,4-dinitrobenzene and 1,3-dinitrobenzene molecules: ab initio and normal coordinate analyses treatments

Ab initio quantum chemical computations have been performed for 1,4-dinitrobenzene and 1,3-dinitrobenzene and their (15)N isotope-labeled derivatives at the 3-21G, 6-31G and 6-31G(*) basis set levels. The results have been compared with the experiment: the best coincidence has been obtained by the 6-31G method for 1,4-dinitrobenzene and by the 6-31G(*) method for 1,3-dinitrobenzene. The isotopic frequency shifts of the nitro group bands, induced by (14)N-(15)N labeling and the vibrational coupling between the nitro groups, calculated with GAMESS software, have been in accordance with the measured ones. A normal coordinate analysis has been carried out for the 1,4-dinitro- benzene molecule using a slightly modified force field of 1,3-dinitrobenzene, reasonable fit between the observed and calculated frequencies has been obtained. The results support the conclusion that the electronic contribution to the vibrational coupling between the two nitro groups in the studied dinitrobenzenes is negligible: thus, it may be suggested that the dynamical coupling will predominate.     

Conformational Composition and Spin Properties of β-Phosphone Nitroxyl Radicals in Gaseous Phase and Solutions

Abstract

As part of our studies of the extremely high sensitivity to substitutients A,B,R, temperature and solvents of phosphorus hyperfine splitting in ESR spectra of β-phosphone nitroxyl radicals ABP(O)CR₂N(Î)Bu¹, we have studied them theoretically. The conformational composition and spin delocalization in H₂(O)CH₂N(Î)Í, (ÌåÌ)₂P(O)CÌå₂N(Î)Ìå and (ÍÎ)₂P(O)CÌ?₂N(Î)Bu' has been studied by UHF/6-31G& and MNDO (GAMESS). Solvent effects on conformational equilibrium, rotational and inversion barriers, charge distributions and spin properties of β-phosphone nitroxyl were calculated by SCRF UHF/6-31G* method.     

(CH2)2N和(CH3)2NH+的密度泛函理论计算

用密度泛函理论UB3LYP／6-31G（d,p）方法研究了二甲胺自由基（CH3）2N及其质子化离子（CH3）2NH 的构象和超精细结构．在由构象研究得到的两种自由基的最稳定结构上,用密度泛函的UB3LYP和UB3PW91方法及从头计算的UHF、UMP2（FULL）方法计算了α-质子、β-质子和N核上的超精细偶合常数A（Hα）、A（Hβ）和A（N）结果表明：两种自由基中甲基内旋转的位垒均很低,分别为0．46kJ·mol-1（（CH3）2NH ）和2．6kJ·mol-1（（CH3）2N）．UB3LYP／6-31G（d,p）和UB3PW91／6-31G（d,p）计算的A（Hα）、A（Hβ）和A（N）与ESR实验结果符合得很好,UMP2／6-31G（d,p）方法的计算值与实验值符合得也较好．     

Theoretical study of the energetics of the reactions of triplet dioxygen with hydroquinone, semiquinone, and their protonated forms: relation to the mechanism of superoxide generation in the respiratory chain

One-electron reduction of the dioxygen molecule by the reduced form of mitochondrial ubiquinones (Q) of the NADH dehydrogenase (complex I) and mitochondrial cytochrome bc1 (complex III) is believed to be the main source of the superoxide anion radical O2*- and the hydroperoxide radical OOH*. In this work, we modeled the energetics of four possible reactions of the triplet ((3)Sigma(g)) dioxygen-molecule reduction by fully reduced and protonated ubiquinone (QH2; reaction 1), its deprotonated form (QH-; reaction 2), the semiquinone radical (QH*; reaction 3), and the semiquinone anion radical (Q*-; reaction 4), by means of ab initio calculations with the 6-31G(d) and 6-31+G(d) basis set in the restricted open-shell Hartree-Fock (ROHF), unrestricted Hartree-Fock (UHF), and complete active space self-consistent field (CASSCF) with dynamic correlation [at the second-order M?ller-Plesset (MP2) or multiple reference M?ller-Plesset (MRMP), respectively] schemes and the basis set superposition error (BSSE) correction included, as well as semiempirical AM1 and PM3 calculations in the UHF and ROHF schemes. 2-Butene-1,4-dione and p-benzoquinone were selected as model compounds. For the reduced forms of both compounds, reaction 1 turned out to be energetically unfavorable at all levels of theory, this agreeing with the experimentally observed diminished reductive properties of hydroquinone derivatives at low pH. For 2-butene-1,4-dione treated at the most advanced MRMP/CASSCF/6-31+G(d) level, the energies of reactions 1-4 are 4.7, -34.3, -15.0, and -4.1 kcal/mol, respectively. This finding suggests that reactions 2 and 3 are the most likely mechanisms of electron transfer to molecular oxygen in aprotic environments and that proton transfer is involved in this process. Nearly the same energies of reactions 2 and 3 were calculated at the MRMP/CASSCF/6-31+G(d) level for reduced forms of p-benzoquinone. Inclusion of diffuse functions in the basis set and dynamic correlation at the CASSCF level appears essential. Because deprotonated ubiquinol is unlikely to exist in physiological environments, reaction 3 appears to be the most likely mechanism of one-electron reduction of oxygen; however, if oxygen can penetrate cytochrome bc1 as far as the Q(o) center where ubiquinol can be deprotonated, reaction 2 can also come into play. The energies of reactions 2 and 3 calculated at the MRMP/CASSCF/6-31+G(d) level are most closely reproduced in the ab initio and semiempirical UHF PM3 calculations. Additional semiempirical calculations on more realistic models of ubiquinone, 2,3-dimethoxy-6-methyl-p-benzoquinone and 2,3-dimethoxy-5-isoprenyl-6-methyl-p-benzoquinone, gave qualitatively the same relations between the energies of reactions 2 and 3 as those carried out for p-benzoquinone species, thereby suggesting that this method could be used in studying electron-transfer reactions from reduced quinone derivatives to molecular oxygen in more complex systems, such as a model of the Q(o) site of cytochrome bc1, where applying ab initio methods is unfeasible.     

取代基对氮自由基稳定性影响的ab initio研究

用UHF／4-31G基组，全构型优化，研究了九个氮自由基NHR（R＝CH3，CF3，CCl3，CN，NH2，CHO，OH，COOH，F）的构型和稳定化能ΔE，它们的ΔE分别为25．24，－38．53，－20．59，21．46，19．96，58．82，73．69，31．75，63．85kJ•mol-1，表明，除CF3，CCl3以外，其余七个取代基对氮自由基起稳定化作用．     

环已二烯与丙烯加成反应的正则速率常数

采用从头算ＵＨＦ方法在６－３１Ｇ基组上研究了１，３－环己二烯与丙烯的加成反应所有可能的反应通道，优化得到了反应途径的所有驻点的几何构型。该反应可按分步和协同过程进行，分步过程的通道是能量最为有利的反应途径。对于分步过程形成ｅｘｏｘ和ｅｎｄｏ产物各存在两类分步加成机理，即“由丙烯双键端位碳先进攻环己二烯”机理和“由丙烯中心碳原子先进攻环己二烯”机理，共四条反应通道，它们均为途径一双自由基中间体的分步     

Theoretical elucidation of kinetic and thermodynamic control of radical addition regioselectivity

The cyclizations of two structurally similar 2-oxo-5-hexenyl-type radicals have been investigated by ab initio and density functional (UB3LYP/6-31+G**//UHF/6-31G* and UB3LYP/6-31G*//UB3LYP/6-31G*) calculations. The origin of apparently contradictory reports of 6-endo and 5-exo cyclizations is determined. Kinetic control favors 6-endo cyclization, while thermodynamic control gives 5-exo cyclization, and the observation of different products from different research groups arises from the difference in experimental conditions used by the two groups. The outcome of a new cyclization reaction was predicted by using these theoretical techniques. Kinetic control is predicted to yield exclusively the products of 6-endo cyclization, while thermodynamic control would lead to an approximately equal mixture of one 6-endo and one 5-exo cyclized product. Experimental studies revealed that the reaction yields only the products of 6-endo cyclization through kinetic control.     

有机双自由基体系的磁性偶合规律

采用量子化学从头算UHF方法对系列有机双自由基体系的基态自旋耦合规律进行研究, 进一步证实了自由基之间在共轭体系中出现铁磁性耦合的拓扑规则,统一了关于自由基耦合规律的几种解释, 为有机磁性材料的实验合成提供了理论指导。