Theoretical study of the enol imine <--> enaminone tautomeric equilibrium in organic solvents

The tautomeric enol imine <--> enaminone (phenol <--> quinone) equilibrium of the 1-hydroxy-2-naphthaldehyde Schiff base (2-phenyliminomethyl-naphthalen-1-ol) was investigated by density functional theory (B3LYP) and ab initio (MP2) methods in the IEF-PCM polarizable continuum dielectric solvent approximation and by a combined ab initio + FEP/MC study by considering an explicit solvent model. Special emphasis was put on the effect of solvation on this equilibrium by using an apolar (CCl4), polar aprotic (CH3CN), and polar protic (CH3OH) solvent. Compared with experimental tautomerization Gibbs free energies, the IEF-PCM/B3LYP calculations apparently overestimate the stability of the quinone form both when the 6-31G(d,p) and the 6-311++G(d,p) basis sets are applied. IEF-PCM/MP2 studies with the above basis sets predict the preference of the aromatic phenol tautomer, in contrast to the experiment in methanol and acetonitrile solvent. Calculation of the total relative free energy as DeltaG(tot) = DeltaE(int)(IEF-PCM/QCISD(T)/6-31G(d)) + DeltaG(solv, FEP/MC) + DeltaG(thermal) provided agreement with the experimental values up to 0.6 kcal/mol in the three solvents, and the predominant tautomer was always correctly predicted. In-solution relevant atomic charges, derived by a fit to the molecular electrostatic potential generated by the IEF-PCM/B3LYP/6-31G(d,p) wave function, show strong dependence on the fitting procedure (CHELPG or RESP) and are fairly insensitive to the chemical nature of the actual solvent. Use of the CHELPG charges in FEP/MC simulations revealed to be superior in comparison with the use of the RESP charge set.     

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

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

低硝酸盐离子N2O22-和亚硝酸HONO的异构体及其异构化反应的量化研究

采用Gaussian-98程序进行，在HF／6—31G(d)，B3LYP／6—31G(d)和MP2／6—31G(d)水平下优化分子结构并寻找过渡态，对于MP2／6—31G(d)结果在QCISD(T，E4T)，MP4／6—311 G(d，p)，MP4／6—311 G(2df，p)水平下重新计算能级．并用内禀反应坐标(IRC)法研究了N2O2^2-和亚硝酸HONO的异构化反应机理。     

Ab initio Study on Structures and Isomerization of Magnesium Fluorosilylenoid H2SiFMgF

The structures and isomerization of magnesium fluorosilylenoid H2SiFMgF were investigated by ab initio molecular orbital theory. Four equilibrium structures and three isomeric transition states were located and fully optimized at the B3LYP/6-31G(d,p) and G3MP2B3 levels, respectively. Based on the B3LYP/6-31G(d,p) optimized geometries, harmonic frequencies of various structures were obtained and 29Si chemical shifts were calculated. The solvent effects were investigated by means of the polarizable continuum model using THF as a solvent at B3LYP/6-31G(d,p) level. Isomerization paths for isomers were confirmed by intrinsic reaction coordinate calculations. The calculated results show that tetrahedral structure has the lowest energy and is the most stable; tetrahedral, three-membered ring, and p-complex structures are suggested to be the experimentally detectable ones; and σ-complex structure has the highest energy and will not exist.     

A dual-level ab initio and hybrid density functional theory dynamics study on the unimolecular decomposition reaction C2H5O-->CH2O + CH3

We present a direct ab initio and hybrid density functional theory dynamics study of the thermal rate constants of the unimolecular decomposition reaction of C2H5O-->CH2O + CH3 at a high-pressure limit. MPW1K/6-31+G(d,p), MP2/6-31+G(d,p), and MP2(full)/6-31G(d) methods were employed to optimize the geometries of all stationary points and to calculate the minimum energy path (MEP). The energies of all the stationary points were refined at a series of multicoefficient and multilevel methods. Among all methods, the QCISD(T)/aug-cc-pVTZ energies are in good agreement with the available experimental data. The rate constants were evaluated based on the energetics from the QCISD(T)/aug-cc-pVTZ//MPW1K/6-31+G(d,p) level of theory using both microcanonical variational transition state theory (microVT) and RRKM theory with the Eckart tunneling correction in the temperature range of 300-2500 K. The calculated rate constants at the QCISD(T)/aug-cc-pVTZ/MPW1K/6-31+G(d,p) level of theory are in good consistent with experimental data. The fitted three-parameter Arrhenius expression from the microVT/Eckart rate constants in the temperature range 200-2500 K is k = 2.52 x 10(12)T(0.41)e(-8894.0/T) s(-1). The falloff curves of pressure-dependent rate constants are performed using master-equation method within the temperature range of 391-471 K. The calculated results are in good agreement with the available experimental data.     

Hydrogen bonding in acetylacetaldehyde: Theoretical insights from the theory of atoms in molecules

All the possible conformations of tautomeric structures (keto and enol) of acetylacetaldehyde (AAD) were fully optimized at HF, B3LYP, and MP2 levels with 6‐31G(d,p) and 6‐311++G(d,p) basis sets to determine the conformational equilibrium. Theoretical results show that two chelated enol forms have extra stability with respect to the other conformers, but identification of global minimum is very difficult. The high level ab initio calculations G2(MP2) and CBS‐QB3) also support the HF conclusion. It seems that the chelated enol forms have equal stability, and the energy gap between them is probably lies in the computational error range. Finally, the analysis of hydrogen bond in these molecules by quantum theory of atoms in molecules (AIM) and natural bond orbital (NBO) methods fairly support the ab initio results. © 2009 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

Transition state structure and energetics of the N(2)O + X (X = Cl,Br) reactions

The structural and vibrational properties of the transition state of the N(2)O + X (X = Cl,Br) reactions have been characterized by ab initio methods using density functional theory. We have employed Becke's hybrid functional (B3LYP), and transition state optimizations were performed with 6-31G(d), 6-311G(2d,2p), 6-311+G(3d,2p), and 6-311+G(3df,2p) basis sets. For the chlorine atom reaction the coupled-cluster method (CCSD(T)) with 6-31G(d) basis set was also used. All calculations resulted in transition state structures with a planar cis arrangement of atoms for both reactions. The geometrical parameters of transition states at B3LYP are very similar, and the reaction coordinates involve mainly the breaking of the N-O bond. At CCSD(T)/6-31G(d) level a contribution of the O-Cl forming bond is also observed in the reaction coordinate. In addition, several highly accurate ab initio composite methods of Gaussian-n (G1, G2, G3), their variations (G2(MP2), G3//B3LYP), and complete basis set (CBS-Q, CBS-Q//B3LYP) series of models were applied to compute reaction energetics. All model chemistries predict exothermic reactions. The G3 and G2 methods result in the smallest deviations from experiment, 1.8 and 0 kcal mol(-1), for the enthalpies of reaction for N(2)O reaction with chlorine and bromine, respectively. The G3//B3LYP and G1 methods perform best among the composite methods in predicting energies of the transition state, with a deviation of 1.9 and 3.0 kcal mol(-1), respectively, in the activation energies for the above processes. However, the B3LYP/6-311+G(3df,2p) method gives smaller deviations of 0.4 and -1.0 kcal mol(-1), respectively. The performance of the methodologies applied in predicting transition state energies was analyzed.     

Oxo-hydroxy tautomerism of 5-fluorouracil: water-assisted proton transfer

Post-Hartree-Fock ab initio quantum chemical calculations were performed for 5-fluorouracil in the gas phase and in a three-water cluster. Full geometry optimizations of the 5-fluorouracil-water complexes were carried out at the MP2/6-31+G(d,p) level of theory. MP4/6-31+G(d,p)//MP2/6-31+G(d,p) and MP4/6-31++G(d,p)//MP2/6-31+G(d,p) single-point calculations were performed to obtain more accurate energies. In water solution, 5-fluorouracil exists mainly in the 2,4-dioxo form (A). We propose that the populations of the 2-hydroxy-4-oxo (B) and 4-hydroxy-2-oxo (D) tautomers are 1 x 10(-4)% and 3.9 x 10(-8)%, respectively, on the basis of the relative stabilities of the tautomers calculated at the MP4/6-31++G(d,p)//MP2/6-31+G(d,p) level of theory. A profound difference between isolated and hydrated 5-fluorouracil is noted for the height of the tautomerization barrier. In the absence of water, the process of proton transfer is very slow. The addition of water molecules decreases the barrier by 2.3 times, making the process much faster. The minimum energy path (MP2/6-31+G(d,p)) for water-assisted proton transfer in trihydrated 5-fluorouracil was followed. CNDO/S-CI calculations predict singlet pi-pi(*) electron transitions at 312 nm for B and at 318 nm for D. The fluorescence spectrum of 5-fluorouracil in water confirms the presence of the hydroxy tautomer.     

Spectra and structure of silicon containing compounds. XXVII. Raman and infrared spectra, conformational stability, vibrational assignment and ab initio calculations of vinyldichlorosilane

The infrared (3200-30 cm(-1) spectra of gaseous and solid and the Raman spectra of liquid (3200-30 cm(-1), with quantitative depolarization values, and solid vinyldichlorosilane, CH2=CHSiHCl2, have been recorded. Both the gauche and the cis conformers have been identified in the fluid phases. Variable temperature (105-150 degrees C) studies of the infrared spectra of the sample dissolved in liquid krypton have been carried out. From these data the enthalpy difference has been determined to be 20 +/- 5 cm(-1) (235 +/- 59 J mol(-1) with the gauche conformer the more stable rotamer. It was not possible to obtain a single conformer in the solid even with repeated annealing of the sample. The experimental enthalpy difference is in agreement with the prediction from MP2/6-311 + G(2d,2p) ab initio calculations with full electron correlation. However, when smaller basis sets, i.e. 6-31G(d) and 6-311 + G(d,p) were utilized the cis conformer was predicted to be the more stable form. Complete vibrational assignments are proposed for both conformers based on infrared contours, relative infrared and Raman intensities, depolarization values and group frequencies, which are supported by normal coordinate calculations utilizing the force constants from ab initio MP2/6-31G(d) calculations. From the frequencies of the Si-H stretches, the Si-H bond distance of 1.474 A has been determined for both the gauche and the cis conformers. Complete equilibrium geometries have been determined for both rotamers by ab initio calculations employing the 6-31G(d), 6-311 + G(d,p) and 6-311 + (2d,2p) basis sets at level of Hartree-Fock (RHF) and/or Moller Plesset to the second order (MP2) with full electron correlation. The potential energy terms for the conformer interconversion have been obtained from the MP2/6-31G(d) calculations. The results are discussed and compared with those obtained for some similar molecules.     

Ab initio方法研究CH~3+OClO反应的可能通道

利用abinitio方法研究了CH~3+OClO反应的三个可能通道,首次应用UMP2(full)/6-31G(d,p)方法得到各反应物、产物、中间物及过渡态的优化构型和谐振频率;然后采用G2MP2理论计算各通道反应焓变和势垒高度。理论计算表明产物通道CH~2O+HOCl是最可能发生的途径,反应放热为443.80kJ·mol^-^1。可能的反应过程为：CH~3和OClO自由基先经无垒过程生成了一个富能中间物,继而通过较低的势垒解离成HOCl+H~2CO。     

Computational study of the deamination reaction of cytosine with H2O and OH-

The mechanism for the deamination reaction of cytosine with H(2)O and OH(-) to produce uracil was investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), and B3LYP/6-31G(d) levels and for anions at the B3LYP/6-31+G(d) level. Single-point energies were also determined at B3LYP/6-31+G(d), MP2/GTMP2Large, and G3MP2 levels of theory. Thermodynamic properties (DeltaE, DeltaH, and DeltaG), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway that was investigated. Intrinsic reaction coordinate analysis was performed to characterize the transition states on the potential energy surface. Two pathways for deamination with H(2)O were found, a five-step mechanism (pathway A) and a two-step mechanism (pathway B). The activation energy for the rate-determining steps, the formation of the tetrahedral intermediate for pathway A and the formation of the uracil tautomer for pathway B, are 221.3 and 260.3 kJ/mol, respectively, at the G3MP2 level of theory. The deamination reaction by either pathway is therefore unlikely because of the high barriers that are involved. Two pathways for deamination with OH(-) were also found, and both of them are five-step mechanisms. Pathways C and D produce an initial tetrahedral intermediate by adding H(2)O to deprotonated cytosine which then undergoes three conformational changes. The final intermediate dissociates to product via a 1-3 proton shift. Deamination with OH(-), through pathway C, resulted in the lowest activation energy, 148.0 kJ/mol, at the G3MP2 level of theory.     

Mechanisms for the deamination reaction of cytosine with H2O/OH(-) and 2H2O/OH(-): a computational study

Mechanisms for the deamination reaction of cytosine with H 2O/OH (-) and 2H 2O/OH (-) to produce uracil were investigated using ab initio calculations. Optimized geometries of reactants, transition states, intermediates, and products were determined at MP2 and B3LYP using the 6-31G(d) basis set and at B3LYP/6-31+G(d) levels of theory. Single point energies were also determined at MP2/G3MP2Large and G3MP2 levels of theory. Thermodynamic properties (Delta E, Delta H, and Delta G), activation energies, enthalpies, and free energies of activation were calculated for each reaction pathway investigated. Intrinsic reaction coordinate (IRC) analysis was performed to characterize the transition states on the potential energy surface. Seven pathways for the deamination reaction were found. All pathways produce an initial tetrahedral intermediate followed by several conformational changes. The final intermediate for all pathways dissociates to product via a 1-3 proton shift. The activation energy for the rate-determining step, the formation of the tetrahedral intermediate for pathway D, the only pathway that can lead to uracil, is 115.3 kJ mol (-1) at the G3MP2 level of theory, in excellent agreement with the experimental value (117 +/- 4 kJ mol (-1)).     

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.     

Semiempirical PM5 molecular orbital study on chlorophylls and bacteriochlorophylls: comparison of semiempirical, ab initio, and density functional results

The semiempirical PM5 method has been used to calculate fully optimized structures of magnesium-bacteriochlorin, magnesium-chlorin, magnesium-porphin, mesochlorophyll a, chlorophylls a, b, c(1), c(2), c(3), and d, and bacteriochlorophylls a, b, c, d, e, f, g, and h with all homologous structures. Hartree-Fock/6-31G* ab initio and density functional B3LYP/6-31G* methods were used to optimize structures of methyl chlorophyllide a, chlorophyll c(1), and methyl bacteriochlorophyllides a and c for comparison. Spectroscopic transition energies of the chromophores and their 1:1 or 1:2 solvent complexes were calculated with the Zindo/S CIS method. The self-consistent reaction field model was used to estimate solvent shifts. The PM5 calculations predict planar structure of the porphyrin ring and central position of the four coordinated magnesium atoms in all pigments studied, in accord with the experimental, ab initio, and density functional results, a significant improvement as compared to the older semiempirical PM3 approach. Only small differences in PM5 and B3LYP/6-31G* or Hartree-Fock/6-31G* minimum energy geometries of the reference molecules were observed. Calculations show that in 1:1 solvent complexes, where the magnesium atom is five coordinated, the magnesium atom is shifted out of the plane of the porphyrin ring towards the solvent molecule, while the hexa coordinated 1:2 complexes are again planar. The PM5 method gives atomic charges that are comparable with those obtained from the Hartree-Fock/6-31G* and B3LYP/6-31G* calculations. The single point ZINDO/S CIS calculations with PM5 minimum energy structure gave excellent correlations between calculated and experimental transition energies of the chlorophylls and bacteriochlorophylls studied. Such correlations may be used for prediction of transition energies of the chromophores in protein binding sites. Calculations also predict existence of dark electronic states below the main Soret absorption band in all chromophores studied. The results suggest that the semiempirical PM5 method is a fairly reliable and computationally efficient method in predicting molecular parameters of porphyrin-like molecules.     

Theoretical study of the gas-phase thermolysis reaction of alkyl (ethyl,isopropyl, and tert-butyl) N,N-dimethylcarbamates and N,N-diethylcarbamates

Theoretical studies on the thermolysis in the gas phase of alkyl N,N-dialkylcarbamates were carried out using ab initio theoretical methods, at the MP2/6-31G(d), MP2/6-31++G(d,p) and MP2/6-311++G(2d,p)//MP2/6-31G(d) levels. The reactions have two steps: the first one corresponds to the formation of an alkene and a neutral dialkylcarbamic acid intermediate via a six-membered cyclic transition state; the second one is the decarboxylation of this intermediate via a four-membered cyclic transition state, leading to carbon dioxide and the corresponding dialkylamine. The progress of the reactions was followed by means of the Wiberg bond indices. The results indicate that the transition states have character intermediate between reactants and products, and the calculated synchronicities show that the reactions are slightly asynchronous. The bond-breaking processes are more advanced than the bond-forming ones, indicating a bond deficiency in the transition states. The rate constants calculated for all the reactions agree very well with the available experimental data.From the Proceedings of the 28th Congreso de Químicos Teóricos de Expresión Latina (QUITEL 2002)     

Direct dynamic study on the hydrogen abstraction reaction CH3CN + OH-->CH2CN + H2O

The reaction of acetonitrile with hydroxyl has been studied using the direct ab initio dynamics methods. The geometries, vibrational frequencies of the stationary points, as well as the minimum energy paths were computed at the BHandHLYP and MP2 levels of theory with the 6-311G(d, p) basis set. The energies were further refined at the PMP4/6-311+G(2df, 2pd) and QCISD(T)/6-311+G(2df, 2pd) levels of theory based on the structures optimized at BHandHLYP/6-311G(d, p) and MP2/6-311G(d, p) levels of theory. The Polyrate 8.2 program was employed to predict the thermal rate constants using the canonical variational transition state theory incorporating a small-curvature tunneling correction. The computed rate constants are in good agreement with the available experimental data.     

Thermochemical properties of molecules and cations of MgNnHm (n=1,2 and m=1–6): enthalpies of formation, proton affinities and ionization energies

Ab initio molecular orbital calculations are reported for small neutral molecules and cations containing magnesium, nitrogen and hydrogen. Structures have been optimized using gradient techniques at B3LYP/6-31+G(d) and at MP2(full)/6-311++G(d,p). Single-point calculations are reported at QCISD(T)(full)/6-311++G(2df,p) and at CCSD(T)(full)/6-311++G(2df,p) levels using geometries optimized at MP2(full)/6-311++G(d,p). Standard enthalpies of formation at 298 K have been calculated at these two higher levels of theory. Other thermochemical properties calculated include ionization energies and proton affinities. The binding enthalpies of ammonia to Mg⁺, MgNH₂⁺ and MgNH₃⁺ are also reported.     

The determination of the equilibrium structures of oxygen,ozone, and hydrogen peroxide using the ab initio and density functional theory methods

The geometries and energies of small oxygen containing molecules are studied by both the ab initio and density functional theory (DFT) methods. The RHF, MP2, and QCISD(T) ab initio methods, BHandH, BHandHLYP, BeckeSLYP, Becke3P86 DFT hybrid methods, BLYP, and the BP86 non-local DFT methods with the 3-21G¹, 6-31G(d,p), 6-311 + G(2d,2p) and 6-311 + + G(3df,3pd) basis sets were used for the computational study. The obtained results from the different methods were compared to the experimental values. The suitability of the DFT methods for reproducing experimental data were discussed.     

优化几何构型对高级别从头算能量的影响

对《ＣＲＣ物理与化学手册》（第７７版）中第三周期以前的无机双原子分子,当其理论优化构型的相对误差大于２％时,分别在实验构型和最大偏差的理论构型下,计算了ＱＣＩＳＤ（Ｔ）／６－３１１＋Ｇ（３ｄｆ,２ｐ）能量并作了比较。结果表明,大多数能量的差别在４．２ｋＪ．ｍｏｌ＾－１以内。由此说明,目前一般采用的构型优化理论方法,多数情况下不至于明显影响单点高级别从头算的计算精度。同时还发现,Ｇ２（ＱＣＩ）的高级     

Vitamin E models. Conformational analysis and stereochemistry of tetralin, chroman, thiochroman and selenochroman

Tetralin, chroman as well as its' S and Se containing congeners were subjected to ab initio (RHF/3-21G and RHF/6-31G(d)) and DFT (B3LYP/6-31G(d)) computation. Molecular geometries and the activation energies for ring inversions were determined with full geometry optimizations.