Conformers and intramolecular hydrogen bonding of the oxalic acid monomer and its anions

Density functional theory, B3LYP/6‐31G** and B3LYP/6‐311+G(2d,p), and ab initio MP2/6‐31G** calculations have been carried out to investigate the conformers, transition states, and energy barriers of the conformational processes of oxalic acid and its anions. QCISD/6‐31G** geometrical optimization is also performed in the stable forms. Its calculated energy differences between the two most stable conformers are very near to the related observed value at 7.0 kJ/mol. It is found that the structures and relative energies of oxalic acid conformers predicted by these methods show similar results, and that the conformer L1 (C_2h) with the double‐interfunctional‐groups hydrogen bonds is the most stable conformer. The magnitude of hydrogen bond energies depends on the energy differences of various optimized structures. The hydrogen bond energies will be about 32 kJ/mol for interfunctional groups, 17 kJ/mol for weak interfunctional groups, 24 kJ/mol for intra‐COOH in (COOH)₂, and 60 kJ/mol for interfunctional groups in (COOH)COO⁻¹ ion if calculated using the B3LYP/6‐311+G(2d,p) method. © 2000 John Wiley & Sons, Inc. Int J Quant Chem 76: 541–551, 2000     

Potential energy scans and vibrational assignments of cyclopropanecarboxylic acid and cyclopropanecarboxamide

The structural stability and internal rotations in cyclopropanecarboxylic acid and cyclopropanecarboxamide were investigated by the DFT-B3LYP and the ab initio MP2 calculations using 6-311G** and 6-311+G** basis sets. The computations were extended to the MP4//MP2/6-311G** and CCSD(T)//MP2/6-311G** single-point calculations. From the calculations the molecules were predicted to exist predominantly in the cis (C=O group eclipses the cyclopropane ring) with a cis-trans barrier of about 4-6kcal/mol. The OCOH torsional barrier in the acid was estimated to be about 12-13kcal/mol while the corresponding OCNH torsional barrier in the amide was calculated to be about 20kcal/mol. The equilibrium constant k for the cis<-->trans interconversion in cyclopropanecarboxylic acid was calculated to be 0.1729 at 298.15K that corresponds to an equilibrium mixture of about 85% cis and 15% trans. The vibrational frequencies were computed at the DFT-B3LYP level. Normal coordinate calculations were carried out and potential energy distributions were calculated for the low energy cis conformer of the molecules. Complete vibrational assignments were made on the basis of normal coordinate calculations and comparison with experimental data of the molecules.     

Theoretical study on keto–enol tautomerism and isomerization in pyruvic acid

Isomerization and tautomerism of 12 isomers of pyruvic acid including 4 keto and 8 enol forms were studied at the MP2 and B3LYP levels of theory using 6‐311++G(2df,p) basis set, separately. Activation energy (E_a), imaginary frequency (υ), and Gibbs free energy (ΔG^#) of the considered isomerization and tautomerism reactions were calculated. Interconversion of the enol forms proceeds through two paths: (i) proton transfer and (ii) internal rotation. Activation energies for the proton transfer paths were in the range of 125–145 kJ/mol and for the internal rotation paths were in the range of 5–45 kJ/mol. Keto–enol tautomerism of pyruvic acid proceeds only through proton transfer route and their activation energies were in the range of 200–300 kJ/mol. Effect of microhydration on the transition state structures and activation energies was also investigated. It was found that the presence of a water molecule catalyzes the isomerization and tautomerism reactions of pyruvic acid so that the activation energies decrease. © 2013 Wiley Periodicals, Inc.     

Theoretical Study on Intermolecular Interactions and Thermodynamic Properties of Nitroamine Dimers

Ab initio self-consistent field(SCF) and Mφller-Plesset correlation correction methods employing 6-31G^** basis set have been applied to the optimizations of nitroamine dimers.The binding energies have been corrected for the basis set superposition error (BSSE) and the zero-point energy.Theree optimized dimers have been obtained.The BSSE corrected binding energy of the most stable dimer is predicted to be -31.85kJ/mol at the MP4/6-31G^**//MP2/6-31G^** level.The energy barriers of the Walden conversion for -NH2 group are 19.7kJ/mol and 18.3kJ/mol for monomer and the most stable dimer,respectively.The molecular interaction makes the internal rotation around N1-N2 even more difficult.The thermodynamic properties of nitroamine and its dimers at different temperatures have been calculated on the basis of vibrational analyses.The change of the Gibbs free energy for the aggregation from monomer to the most stable dimer at standard pressure and 298.2 K is predicted to be 14.05kJ/mol.     

Mechanistic study of the deamination reaction of guanine: a computational study

The mechanism for the deamination of guanine with H(2)O, OH(-), H(2)O/OH(-) and for GuaH(+) with H(2)O has been investigated using ab initio calculations. Optimized geometries of the reactants, transition states, intermediates, and products were determined at RHF/6-31G(d), MP2/6-31G(d), B3LYP/6-31G(d), and B3LYP/6-31+G(d) levels of theory. Energies were also determined at G3MP2, G3MP2B3, G4MP2, and CBS-QB3 levels of theory. Intrinsic reaction coordinate (IRC) calculations were performed to characterize the transition states on the potential energy surface. Thermodynamic properties (ΔE, ΔH, and ΔG), activation energies, enthalpies, and Gibbs free energies of activation were also calculated for each reaction investigated. All pathways yield an initial tetrahedral intermediate and an intermediate in the last step that dissociates to products via a 1,3-proton shift. At the G3MP2 level of theory, deamination with OH(-) was found to have an activation energy barrier of 155 kJ mol(-1) compared to 187 kJ mol(-1) for the reaction with H(2)O and 243 kJ mol(-1) for GuaH(+) with H(2)O. The lowest overall activation energy, 144 kJ mol(-1) at the G3MP2 level, was obtained for the deamination of guanine with H(2)O/OH(-). Due to a lack of experimental results for guanine deamination, a comparison is made with those of cytosine, whose deamination reaction parallels that of guanine.     

Solvent-assisted catalysis mechanism on Keto–enol tautomerism of cyameluric acid

The keto–enol tautomerism of cyameluric acid, both in gas phase and in water and methanol solution, has been studied at the B3LYP/6-31++g(d,P) level of theory in this paper. The harmonic frequencies of all the structures are calculated. The results show that the transition states of the tautomerism are 4-membered ring conformations in gas phase, whereas 6-membered ring conformations in solution. In the first proton transfer, activation energy ΔE^# is 56.4 and 50.9 kJ/mol for water and methanol solution, respectively, which is much lower than that in gas phase (163.2 kJ/mol). Solvent molecules (water and methanol) produce an important catalytic effect in the tautomerism, especially for methanol-solvated system. NBO analysis shows that there is a strong interaction between cyameluric acid and solvent molecules in transition states. AIM charge analysis indicates that the keto–enol tautomerism shows a certain degree of proton transfer character. From the reaction enthalpy and reaction rate point of view, keto–enol tautomerism in water-solvated and methanol-solvated system is easier than that in gas phase. The keto–enol tautomerisms are endothermic both in gas phase and in solution, so the enol forms are less stable than the keto ones.     

Structure and vibrational assignment of the enol form of 1,1,1-trifluoro-2,4-pentanedione

Molecular structure of 1,1,1-trifluoro-pentane-2,4-dione, known as trifluoro-acetylacetone (TFAA), has been investigated by means of Density Functional Theory (DFT) calculations and the results were compared with those of acetylacetone (AA) and hexafluoro-acetylacetone (HFAA). The harmonic vibrational frequencies of both stable cis-enol forms were calculated at B3LYP level of theory using 6-31G** and 6-311++G** basis sets. We also calculated the anharmonic frequencies at B3LYP/6-31G** level of theory for both stable cis-enol isomers. The calculated frequencies, Raman and IR intensities, and depolarization ratios were compared with the experimental results. The energy difference between the two stable cis-enol forms, calculated at B3LYP/6-311++G**, is only 5.89 kJ/mol. The observed vibrational frequencies and Raman and IR intensities are in excellent agreement with the corresponding values calculated for the most stable conformation, 2TFAA. According to the theoretical calculations, the hydrogen bond strength for the most stable conformer is 57 kJ/mol, about 9.5kJ/mol less than that of AA and about 14.5 kJ/mol more than that of HFAA. These hydrogen bond strengths are consistent with the frequency shifts for OH/OD stretching and OH/OD out-of-plane bending modes upon substitution of CH(3) groups with CF(3) groups. By comparing the vibrational spectra of both theoretical and experimental data, it was concluded that 2TFAA is the dominant isomer.     

Uracil–4-Hydroxyuracil tautomerism revisited

The tautomerism of the uracil-4-hydroxyuracil system is examined with the 3–21G basis set and full geometry optimization using gradient techniques. Single-point calculations at the 6-31G level were also done and the results compared with the 2-pyridone–2-hydroxypyridine system. The best ΔE(taut) obtained in this work is 74.1 kJ/mol.     

苯并氧化呋咱稳定性和异构化的DFT和ab initio研究

运用B3LYP/6-31G(d)密度泛函理论（DFT）方法对苯并氧化呋咱、邻二亚硝基苯及其间的异构化反应进行了计算研究。结果表明,苯并氧化呋咱的分子总能量比邻二亚硝基苯的低;由苯并氧化呋咱异构为邻二亚硝基苯的正向反应活化能（Ea+=51.0kJ/mol）,与文献实测值（58.6kJ/mol）较接近,而其逆向反应活化能（Ea-=4.6kJ/mol）很小,从而揭示了苯并氧化呋咱比邻二亚硝基苯更稳定·此外,进行了HF/3-21G、HF/6-31G(d)和MP2/6-31G(d)//6-31G(d)水平下相应的计算,发现B3LYP-DFT的结果较abinitio为优。谐振动频率的B3LYP/6-31G(d)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

An ab initio study of intermolecular interactions of nitromethane dimer and nitromethane trimer

Different geometries of nitromethane dimer and nitromethane trimer have been fully optimized employing the density functional theory B3LYP method and the 6-31++G** basis set. Three-body interaction energy has been obtained with the ab initio supermolecular approach at the levels of MP2/6-31++G**//B3LYP/6-31++G** and MP2/aug-cc-pVDZ//B3LYP/6-31++G**. The internal rotation of methyl group induced by intermolecular interaction has been observed theoretically. For the optimized structures of nitromethane dimer, the strength of C--H...O--N H-bond ranges from -9.0 to -12.4 kJ mol(-1) at the MP2/aug-cc-pVDZ//B3LYP/6-31++G** level, and the B3LYP method underestimates the interaction strength compared with the MP2 method, while MP2/6-31++G**//B3LYP/6-31++G** calculated DeltaE(C) is within 2.5 kJ mol(-1) of the corresponding value at the MP4(SDTQ)/6-31G**//B3LYP/6-31++G** level. The analytic atom-atom intermolecular potential has been successfully regressed by using the MP2/6-31++G**//B3LYP/6-31++G** calculated interaction energies of nitromethane dimer. For the optimized structures of nitromethane trimer the three-body interaction energies occupy small percentage of corresponding total binding energies, but become important for the compressed nitromethane explosive. In addition, it has been discovered that the three-body interaction energy in the cyclic nitromethane trimer is more and more negative as intermolecular distances decrease from 2.2 to 1.7 A.     

Structure of the boronic acid dimer and the relative stabilities of its conformers

Despite the widespread use of boronic acids in materials science and as pharmaceutical agents, many aspects of their structure and reactivity are not well understood. In this research the boronic acid dimer, [HB(OH)(2)](2), was studied by second-order M?ller-Plesset (MP2) perturbation theory and coupled cluster methodology with single and double excitations (CCSD). Pople split-valence 6-31+G*, 6-311G**, and 6-311++G** and Dunning-Woon correlation-consistent cc-pVDZ, aug-cc-pVDZ, cc-pVTZ, and aug-cc-pVTZ basis sets were employed for the calculations. A doubly hydrogen-bonded conformer (1) of the dimer was consistently found to be lowest in energy; the structure of 1 was planar (C(2h)) at most computational levels employed but was significantly nonplanar (C(2)) at the MP2/6-311++G** and CCSD/6-311++G** levels, the result of an intrinsic problem with Pople-type sp-diffuse basis functions on heavy atoms. The dimerization energy, enthalpy, and free energy for the formation of (1) from the exo-endo conformer of the monomer were -10.8, -9.2, and +1.2 kcal/mol, respectively, at the MP2/aug-cc-pVTZ level. Several other hydrogen-bonded conformers of the dimer were local minima on the potential energy surface (PES) and ranged from 2 to 5 kcal/mol higher in energy than 1. Nine doubly OH-bridged conformers, in which the boron atoms were tetracoordinated, were also local minima on the PES, but they were all greater than 13 kcal/mol higher in energy than 1; doubly H-bridged structures proved to be transition states. MP2 and CCSD results were compared to those from the BLYP, B3LYP, OLYP, O3LYP, PBE1PBE, and TPSS functionals with the 6-311++G** and aug-cc-pVTZ basis sets; the PBE1PBE functional performed best relative to the MP2 and CCSD results. Self-consistent reaction field (SCRF) calculations predict that boronic acid dimerization is less favorable in solution than in vacuo.     

异硫氰酸与甲亚胺环加成反应机理的理论研究

对异硫氰酸与甲亚胺形成四元或六元环产物的环加成反应进行了理论研究,结果表明,一分子异硫氰酸与一分子甲亚胺形成四元环(1:1)产物P1的反应(1)为经过一个两性离子中间体的分步反应,其中第二步为速控步骤,其活化热垒为107.86kJ/mol.此外,反应(1)的中间体还可与另一甲亚胺或异硫氰酸分子继续反应形成两个不同的六元环(1:2或2:1)产物P2或P3;这两个反应均为协同反应,其活化势垒分别为15.88和21.82kJ/mol.这些结果与当异硫氰酸酯与亚胺发生环加成反应时,只有类似于P2和P3的两种六元环产物生成的实验事实一致。     

4-乙酰基-1,3-二氟苯硝化反应机理研究

以4-氨基-1,3-二氟苯、乙酸酐和浓硝酸等为原料,合成了4-乙酰基-5-硝基-1,3-二氟苯,并且采用FTIR和1H NMR对产物进行了表征.采用量子化学方法,探讨了在酸催化条件下4-乙酰基-1,3-二氟苯硝化反应过程中的微观机理,在密度泛函理论的B3LYP/6-311G**水平下搜索得到反应路径,并对各中间体、过渡态和产物进行几何优化、振动频率分析和标准热力学函数计算;反应路径经过IRC验证.在MP2/6-311G**水平上对所有中间体、过渡态和产物进行更精确的能量计算.     

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

Differentiation of α‐COOH from β‐COOH in aspartic acids by N‐phosphorylation

The biomimic reactions of N‐phosphoryl amino acids, which involved intramolecular penta‐coordinate phosphoric‐carboxylic mixed anhydrides, are very important in the study of many biochemical processes. The reactivity difference between the α‐COOH group and β‐COOH in phosphoryl amino acids was studied by experiments and theoretical calculations. It was found that the α‐COOH group, and not β‐COOH, was involved in the ester exchange on phosphorus in experiment. From MNDO calculations, the energy of the penta‐coordinate phosphoric intermediate containing five‐member ring from α‐COOH was 35 kJ/mol lower than that of the six‐member one from β‐COOH. This result was in agreement with that predicted by HF/6‐31G** and B3LYP/6‐31G** calculations. Theoretical three‐dimensional potential energy surface for the intermediates predicted that the transition states 4 and 5 involving α‐COOH or β‐COOH group had energy barriers of ΔE=175.8 kJ?mol^?1 and 210.4 kJ?mol^?1, respectively. So the α‐COOH could be differentiated from β‐COOH intramolecularly in aspartic acids by N‐phosphorylation. © 2001 John Wiley & Sons, Inc. Int J Quant Chem 83: 41–51, 2001     

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)).     

尿酸质子转移反应的理论研究I——分子内质子转移研究

使用量子化学中的Hartree-Fock方法和密度泛函理论中的B3LYP方法，分别在3－21G^*和6－31G（d)水平上，计算了尿酸分子从三羰基异构体向三羟基异构体的转化。结果表明，转化过程经历了单羟基和双羟基异构体2种中间物和3种过渡态时的分子内质子转移（IPT），转移中的H原邻近的N，O和C原子形成了具有四元环结构的过渡态。随着IPT的进行，N－H键逐渐被削弱和断裂，O－H键则逐渐生成。3个反应的活化能分别为190.3kJ/mol,181.4kJ/mol和249.9kJ/mol(B3LYP/6-31G(d))。较高的活化能表明在室温下，无催化剂的IPT难以进行。     

Molecular structure of the trans and cis isomers of metal-free phthalocyanine studied by gas-phase electron diffraction and high-level quantum chemical calculations: NH tautomerization and calculated vibrational frequencies

The molecular structure of the trans isomer of metal-free phthalocyanine (H2Pc) is determined using the gas electron diffraction (GED) method and high-level quantum chemical calculations. B3LYP calculations employing the basis sets 6-31G**, 6-311++G**, and cc-pVTZ give two tautomeric isomers for the inner H atoms, a trans isomer having D2h symmetry and a cis isomer having C2v symmetry. The trans isomer is calculated to be 41.6 (B3LYP/6-311++G**, zero-point corrected) and 37.3 kJ/mol (B3LYP/cc-pVTZ, not zero-point corrected) more stable than the cis isomer. However, Hartree-Fock (HF) calculations using different basis sets predict that cis is preferred and that trans does not exist as a stable form of the molecule. The equilibrium composition in the gas phase at 471 degrees C (the temperature of the GED experiment) calculated at the B3LYP/6-311++G** level is 99.8% trans and 0.2% cis. This is in very good agreement with the GED data, which indicate that the mole fraction of the cis isomer is close to zero. The transition states for two mechanisms of the NH tautomerization have been characterized. A concerted mechanism where the two H atoms move simultaneously yields a transition state of D2h symmetry and an energy barrier of 95.8 kJ/mol. A two-step mechanism where a trans isomer is converted to a cis isomer, which is converted into another trans isomer, proceeds via two transition states of C(s) symmetry and an energy barrier of 64.2 kJ/mol according to the B3LYP/6-311++G** calculation. The molecular geometry determined from GED is in very good agreement with the geometry obtained from the quantum chemical calculations. Vibrational frequencies, IR, and Raman intensities have been calculated using B3LYP/6-311++G**. These calculations indicate that the molecule is rather flexible with six vibrational frequencies in the range of 20-84 cm(-1) for the trans isomer. The cis isomer might be detected by infrared matrix spectroscopy since the N-H stretching frequencies are very different for the two isomers.     

Electronic structure of the peptide linkage. I. Equilibrium geometry and electronic properties of formhydroxamic acid

The equilibrium geometry of formhydroxamic acid has been calculated within the framework of the INDO –MO formalism. Various structural factors are analyzed and discussed in terms of the calculated force constants and charge distribution. The possibility of internal rotation around the C? N bond of formhydroxamic acid has been examined. The potential energy surface for the amide-imide tautomerism is explored by calculating the geometries and characterizing saddle points on that surface. The cyclic and open dimers of formhydroxamic acid are examined and the hydrogenbond energy and length are calculated.