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

Theoretical study on the hydrolytic deamination reaction mechanism of guanine and (H2O)n

The hydrolytic deamination reaction mechanism of guanine (G) and (H₂O)n (n = 1–4) has been theoretically investigated, including solvent effects at the B3LYP/6-31G (d,p) and MP2/6-311G (d,p) levels. The results show that the hydrolytic deamination reaction of G involves two steps. In the first step, a tetra-coordinated intermediate is generated by a hydrolysis reaction, and then deamination reaction is carried out. In the second step, the final products form via a hydrogen transfer. For the hydrolytic deamination reaction of guanine and 3H₂O or 4H₂O, only two molecules of water participate in the reaction, and other water molecule make the transition state more stable as a catalyst. The study on the potential energy surface shows that the deamination reaction of G and nH₂O does not take place because of a higher barrier for the opening system, which is in agreement with the experimental result.     

A computational mechanistic study of the deamination reaction of melamine

A detailed computational study of the deamination reaction of melamine by OH^–, n H₂O/OH^–, n H₂O (where n = 1, 2, 3), and protonated melamine with H₂O, has been carried out using density functional theory and ab initio calculations. All structures were optimized at M06/6‐31G(d) level of theory, as well as with the B3LYP functional with each of the basis sets: 6‐31G(d), 6‐31 + G(d), 6‐31G(2df,p), and 6‐311++G(3df,3pd). B3LYP, M06, and ω B97XD calculations with 6‐31 + G(d,p) have also been performed. All structures were optimized at B3LYP/6‐31 + G(d,p) level of theory for deamination simulations in an aqueous medium, using both the polarizable continuum solvation model and the solvation model based on solute electron density. Composite method calculations have been conducted at G4MP2 and CBS‐QB3. Fifteen different mechanistic pathways were explored. Most pathways consisted of two key steps: formation of a tetrahedral intermediate and in the final step, an intermediate that dissociates to products via a 1,3‐proton shift. The lowest overall activation energy, 111 kJ mol^?1 at G4MP2, was obtained for the deamination of melamine with 3H₂O/OH^?.     

苯并氧化呋咱稳定性和异构化的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)计算还支持了邻二亚硝基苯为苯并氧化呋咱“自-自”互变重排反应的中间体。     

Radical-molecule reactions HCO/HOC + C2H2: mechanistic study

A detailed computational study is performed on the unknown radical-molecule reactions between HCO/HOC and acetylene (C2H2) at the CCSD(T)/6-311G(2d,p)//B3LYP/6-311G(d,p)+ZPVE, Gaussian-3//B3LYP/6-31G(d), and Gaussian-3//MP2(full)/6-31G(d) levels. For the HCO + C2H2 reaction, the most favorable pathway is direct C-addition forming the intermediate HC=CHCH=O followed by a 1,3-H-shift leading to H2C=CHC=O, which finally dissociates to the product C2H3 + CO. The overall reaction barrier is 13.8, 10.5, and 11.3 kcal/mol, respectively, at the three levels. The quasi-direct H-donation process to produce C2H3 + CO with barriers of 14.0, 14.1, and 14.1 kcal/mol is less competitive. Thus only at higher temperatures could the HCO + C2H2 reaction play a role. In contrast, the HOC + C2H2 reaction can barrierlessly generate C2H3 + CO via the quasi-direct H-donation mechanism proceeding via a prereactive complex with OH...C2 hydrogen bonding. This is suggestive of the potential importance of the HOC + C2H2 reaction in both combustion and interstellar processes. However, the direct C-addition channel is much less competitive. For both reactions, the possible formation of the intriguing interstellar molecules propadiene and propynal is also discussed. The present theoretical study represents the first attempt to probe the reaction mechanism between HOC and pi-systems. Future laboratory investigations on both reactions (particularly HOC + C2H2) are recommended.     

硝基甲烷异构化反应势能面的ab initio研究

在B3LYP/6-311+ +G(2d,2p)水平上,优化得到硝基甲烷CH3NO2的10种异构体和23个异构化反应过渡态,并用G2MP2方法进行了单点能计算.根据计算得到的G2MP2相对能量,探讨了CH3NO2势能面上异构化反应的微观机理.研究表明,反应初始阶段的CH3NO2异构化过程具有较高的能垒,其中CH3NO2的两个主要异构化反应通道,即CH3NO2→CH3ONO和CH3NO2→CH2N(O)OH的活化能分别为270.3和267.8 kJ/mol,均高于CH3NO2的C-N键离解能.因而,从动力学角度考虑, CH3NO2的异构化反应较为不利.     

A new scheme for determining the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of glycine and alanine peptides

In this paper a new scheme was proposed to calculate the intramolecular hydrogen-bonding energies in peptides and was applied to calculate the intramolecular seven-membered ring N-H...O=C hydrogen-bonding energies of the glycine and alanine peptides. The density-functional theory B3LYP6-31G(d) and B3LYP6-311G(d,p) methods and the second-order Moller-Plesset perturbation theory MP26-31G(d) method were used to calculate the optimal geometries and frequencies of glycine and alanine peptides and related structures. MP26-311++G(d,p), MP26-311++G(3df,2p), and MP2/aug-cc-pVTZ methods were then used to evaluate the single-point energies. It was found that the B3LYP6-31G(d), MP26-31G(d), and B3LYP6-311G(d,p) methods yield almost similar structural parameters for the conformers of the glycine and alanine dipeptides. MP2/aug-cc-pVTZ predicts that the intramolecular seven-membered ring N-H...O=C hydrogen-bonding strength has a value of 5.54 kcal/mol in glycine dipeptide and 5.73 and 5.19 kcal/mol in alanine dipeptides, while the steric repulsive interactions of the seven-membered ring conformers are 4.13 kcal/mol in glycine dipeptide and 6.62 and 3.71 kcal/mol in alanine dipeptides. It was also found that MP26-311++G(3df,2p) gives as accurate intramolecular N-H...O=C hydrogen-bonding energies and steric repulsive interactions as the much more costly MP2/aug-cc-pVTZ does.     

Binding affinities for models of biologically available potential Cu(II) ligands relevant to Alzheimer's disease: an ab initio study

A systematic study of the binding affinities of the model biological ligands X: = (CH3)2S, CH3S-, CH3NH2, 4-CH3-imidazole (MeImid), C6H5O-, and CH3CO2- to (NH3)i(H2O)3-iCu(II)-H2O (i = 3, 2, 1, 0) complexes has been carried out using quantum chemical calculations. Geometries have been obtained at the B3LYP/ 6-31G(d) level of theory, and binding energies, Delta, relative to H2O as a ligand, have been calculated at the B3LYP/6-311+G(2df,2p)//B3LYP/6-31G(d) level. Solvation effects have been included using the COSMO model, and the relative binding free energies in aqueous solution (Delta) have been determined at pH 7 for processes that are pH dependent. CH3S- (Delta = -16.0 to -53.5 kJ mol(-1)) and MeImid (Delta = -18.5 to -35.2 kJ mol(-1)) give the largest binding affinities for Cu(II). PhO- and (CH3)2S are poor ligands for Cu(II), Delta = 20.6 to -9.7 and 19.8 to -3.7 kJ mol(-1), respectively. The binding affinities for CH3NH2 range from -0.8 to -15.0 kJ mol(-1). CH3CO2- has Cu(II) binding affinities in the ranges Delta = -13.5 to -32.4 kJ mol(-1) if an adjacent OH bond is available for hydrogen bonding and Delta = 10.1 to -4.6 kJ mol(-1) if this interaction is not present. In the context of copper coordination by the Abeta peptide of Alzheimer's disease, the binding affinities suggest preferential binding of Cu(II) to the three histidine residues plus a lysine or the N-terminus. For a 3N1O Cu(II) ligand arrangement, it is more probable that the oxygen ligand comes from an aspartate/glutamate residue side chain than from the tyrosine at position 10. Methionine appears unlikely to be a Cu(II) ligand in Abeta.     

Energetics of the O-H bond and of intramolecular hydrogen bonding in HOC6H4C(O)Y (Y = H, CH3, CH2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds

The energetics of the phenolic O-H bond in a series of 2- and 4-HOC 6H 4C(O)Y (Y = H, CH3, CH 2CH=CH2, C[triple bond]CH, CH2F, NH2, NHCH 3, NO2, OH, OCH3, OCN, CN, F, Cl, SH, and SCH3) compounds and of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y, was investigated by using a combination of experimental and theoretical methods. The standard molar enthalpies of formation of 2-hydroxybenzaldehyde (2HBA), 4-hydroxybenzaldehyde (4HBA), 2'-hydroxyacetophenone (2HAP), 2-hydroxybenzamide (2HBM), and 4-hydroxybenzamide (4HBM), at 298.15 K, were determined by micro- or macrocombustion calorimetry. The corresponding enthalpies of vaporization or sublimation were also measured by Calvet drop-calorimetry and Knudsen effusion measurements. The combination of the obtained experimental data led to Delta f H m (o)(2HBA, g) = -238.3 +/- 2.5 kJ.mol (-1), DeltafHm(o)(4HBA, g) = -220.3 +/- 2.0 kJ.mol(-1), Delta f H m (o)(2HAP, g) = -291.8 +/- 2.1 kJ.mol(-1), DeltafHm(o)(2HBM, g) = -304.8 +/- 1.5 kJ.mol (-1), and DeltafHm(o) (4HBM, g) = -278.4 +/- 2.4 kJ.mol (-1). These values, were used to assess the predictions of the B3LYP/6-31G(d,p), B3LYP/6-311+G(d,p), B3LYP/aug-cc-pVDZ, B3P86/6-31G(d,p), B3P86/6-311+G(d,p), B3P86/aug-cc-pVDZ, and CBS-QB3 methods, for the enthalpies of a series of isodesmic gas phase reactions. In general, the CBS-QB3 method was able to reproduce the experimental enthalpies of reaction within their uncertainties. The B3LYP/6-311+G(d,p) method, with a slightly poorer accuracy than the CBS-QB3 approach, achieved the best performance of the tested DFT models. It was further used to analyze the trends of the intramolecular O...H hydrogen bond in 2-HOC 6H 4C(O)Y evaluated by the ortho-para method and to compare the energetics of the phenolic O-H bond in 2- and 4-HOC 6H 4C(O)Y compounds. It was concluded that the O-H bond "strength" is systematically larger for 2-hydroxybenzoyl than for the corresponding 4-hydroxybenzoyl isomers mainly due to the presence of the intramolecular O...H hydrogen bond in the 2-isomers. The observed differences are, however, significantly dependent on the nature of the substituent Y, in particular, when an intramolecular H-bond can be present in the radical obtained upon cleavage of the O-H bond.     

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.     

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.     

Mechanisms of formation of 8-oxoguanine due to reactions of one and two OH* radicals and the H2O2 molecule with guanine: A quantum computational study

Mechanisms of formation of the mutagenic product 8-oxoguanine (8OG) due to reactions of guanine with two separate OH* radicals and with H2O2 were investigated at the B3LYP/6-31G, B3LYP/6-311++G, and B3LYP/AUG-cc-pVDZ levels of theory. Single point energy calculations were carried out with the MP2/AUG-cc-pVDZ method employing the optimized geometries at the B3LYP/AUG-cc-pVDZ level. Solvent effect was treated using the PCM and IEF-PCM models. Reactions of two separate OH* radicals and H2O2 with the C2 position of 5-methylimidazole (5MI) were investigated taking 5MI as a model to study reactions at the C8 position of guanine. The addition reaction of an OH* radical at the C8 position of guanine is found to be nearly barrierless while the corresponding adduct is quite stable. The reaction of a second OH* radical at the C8 position of guanine leading to the formation of 8OG complexed with a water molecule can take place according to two different mechanisms, involving two steps each. According to one mechanism, at the first step, 8-hydroxyguanine (8OHG) complexed with a water molecule is formed ,while at the second step, 8OHG is tautomerized to 8OG. In the other mechanism, at the first step, an intermediate complexed (IC) with a water molecule is formed, the five-membered ring of which is open, while at the second step, the five-membered ring is closed and a hydrogen bonded complex of 8OG with a water molecule is formed. The reaction of H2O2 with guanine leading to the formation of 8OG complexed with a water molecule can also take place in accordance with two different mechanisms having two steps each. At the first step of one mechanism, H2O2 is dissociated into two OH* groups that react with guanine to form the same IC as that formed in the reaction with two separate OH* radicals, and the subsequent step of this mechanism is also the same as that of the reaction of guanine with two separate OH* radicals. At the first step of the other mechanism of the reaction of guanine with H2O2, the latter molecule is dissociated into a hydrogen atom and an OOH* group which become bonded to the N7 and C8 atoms of guanine, respectively. At the second step of this mechanism, the OOH* group is dissociated into an oxygen atom and an OH* group, the former becomes bonded to the C8 atom of guanine while the latter abstracts the H8 atom bonded to C8, thus producing 8OG complexed with a water molecule. Solvent effects of the aqueous medium on certain reaction barriers and released energies are appreciable. 5MI works as a satisfactory model for a qualitative study of the reactions of two separate OH* radicals or H2O2 occurring at the C8 position of guanine.     

A theoretical study of the fluoride exchange between UO2F+(aq) and UO2 2+(aq)

Experimental data on the thermodynamics and reaction mechanism of the inner-sphere fluoride exchange reaction U17O2(2+) + UO2F+ <==> U17O2F+ + UO2(2+) have been compared with different intimate reaction mechanisms using quantum chemical methods. Two models have been tested that start from the outer sphere complexes, (H2O)[U(A)O2F(OH2)4+]...[U(B)O2(OH2)5(2+)] and [U(A)O2F(OH2)4+]...[U(B)O2(OH2)5(2+)]; the geometry and energies of the intermediates and transition states along possible reaction pathways have been calculated using different ab initio methods, SCF, B3LYP and MP2. Both the experimental data and the theoretical results suggest that the fluoride exchange takes place via the formation and breaking of a U-F-U bridge that is the rate determining step. The calculated activation enthalpy DeltaH( not equal) = 30.9 kJ mol(-1) is virtually identical to the experimental value 31 kJ mol(-1); however this agreement may be a coincidence as we do not expect a larger accuracy than 10 kJ mol(-1) with the methods used. The calculations show that the fluoride bridge is formed as an insertion of U(A)O2)F(OH2)4+ into U(B)O2(OH2)5(2+) followed by a subsequent transfer of water from the first to the second coordination sphere of U(B).     

尿酸质子转移反应的理论研究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难以进行。     

Ab initio calculations on the barrier height for the hydrogen addition to ethylene and formaldehyde. The importance of spin projection

Heats of reaction and barrier heights have been computed for H + CH₂CH₂ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH using unrestricted Hartree-Fock and Møller–Plesset perturbation theory up to fourth order (with and without spin annihilation), using single-reference configuration interaction, and using multiconfiguration self-consistent field methods with 3-21G, 6-31G(d), 6-31G(d,p), and 6-311G(d,p) basis sets. The barrier height in all three reactions appears to be relatively insensitive to the basis sets, but the heats of reaction are affected by p-type polarization functions on hydrogen. Computation of the harmonic vibrational frequencies and infrared intensities with two sets of polarization functions on heavy atoms [6-31G(2d)] improves the agreement with experiment. The experimental barrier height for H + C₂H₄ (2.04 ± 0.08 kcal/mol) is overestimated by 7?9 kcal/mol at the MP2, MP3, and MP4 levels. MCSCF and CISD calculations lower the barrier height by approximately 4 kcal/mol relative to the MP4 calculations but are still almost 4 kcal/mol too high compared to experiment. Annihilation of the largest spin contaminant lowers the MP4SDTQ computed barrier height by 8?9 kcal/mol. For the hydrogen addition to formaldehyde, the same trends are observed. The overestimation of the barrier height with Møller-Plesset perdicted barrier heights for H + C₂H₄ → C₂H₅, H + CH₂O → CH₃O, and H + CH₂O → CH₂OH at the MP4SDTQ /6-31G(d) after spin annihilation are respectively 1.8, 4.6, and 10.5 kcal/mol.     

New insights into an old reaction. High-resolution x-ray powder diffraction of Wiberg's aminoalane intermediate

In accordance with the procedure described by E. Wiberg, Me(3)Al-NH(3) was heated as a bulk material in inert atmosphere to give a colorless liquid which slowly loses methane. Close to the end of this elimination reaction, the melt crystallized to give a microcrystalline powder of (Me(2)AlNH(2))(x)(). The structure of this intermediate has been solved by the method of high-resolution X-ray powder diffraction. The compound crystallizes in the monoclinic space group C2/c with the cell parameters of a = 15.0047(6) A, b = 8.7500(2) A, c = 24.4702(8) A, and beta = 107.290(2) degrees, with eight trimers (Me(2)AlNH(2))(3) per unit cell. These trimers crystallize in a boat conformation in contrast to the known trimers of the same composition where a twist-boat conformation had been found by single crystal determination. Different conformers of (Me(2)AlNH(2))(3) have been investigated by theoretical methods (HF/6-31G(d), B3LYP/6-31G(d), B3LYP/6-311G(d,p), MP2(fc)/6-31G(d), and MP2(fc)/6-311G(d,p)). The twist-boat and the chair conformer correspond to minima at the potential energy surface, whereas the boat conformer corresponds to a first-order transition state (relative energies of 0.45-2.56 kJ/mol (boat) and 6.66-11.91 kJ/mol (chair)). Relaxed scans of the potential energy surface at the HF/6-31G(d) and B3LYP/6-31G(d) levels have shown that the boat conformer (C(s)() symmetry) connects two enantiomers of the twist-boat form (C(2) symmetry).     

Conformational stability, vibrational assignmenents, barriers to internal rotations and ab initio calculations of 2-aminophenol (d 0 and d3)

The Raman (3700-100 cm(-1)) and infrared (4000-400 cm(-1)) spectra of solid 2-aminophenol (2AP) have been recorded. The internal rotation of both OH and NH2 moieties produce ten conformers with either Cs or C1 symmetry. However, the calculated energies as well as the imaginary vibrational frequencies reduce rotational isomerism to five isomers. The molecular geometry has been optimized without any constraints using RHF, MP2 and B3LYP levels of theory at 6-31G(d), 6-311+G(d) and 6-31++G(d,p) basis sets. All calculations predict 1 (cis; OH is directed towards NH2) to be the most stable conformation except RHF/6-31++G(d,p) basis set. The 1 (cis) isomer is found to be more stable than 8 (trans; OH is away from the NH2 moiety and the NH bonds are out-of-plane) by 1.7 kcal/mol (598 cm(-1)) as obtained from MP2/6-31G(d) calculations. Aided by experimental and theoretical vibrational spectra, cis and trans 2AP are coexist in solution but cis isomer is more likely present in the crystalline state. Aided by MP2 and B3LYP frequency calculations, molecular force fields, simulated vibrational spectra utilizing 6-31G(d) basis set as well as normal coordinate analysis, complete vibrational assignments for HOC6H4NH2 and DOC6H4ND2 have been proposed. Furthermore, we carried out potential surface scan, to determine the barriers to internal rotations of NH2 and OH groups. All results are reported herein and compared with similar molecules when appropriate.     

Can quantum-mechanical calculations yield reasonable estimates of hydrogen-bonding acceptor strength? The case of hydrogen-bonded complexes of methanol

The thermodynamics and some vibrational properties of hydrogen-bonded complexes of methanol with 23 hydrogen-bond acceptors (HBAs) have been determined in CCl(4) by FTIR spectrometry. The experimental sample contains carbon, nitrogen, oxygen, sulfur, fluorine, and chlorine organic bases and covers an energetic range of 13 kJ mol(-1) in the basicity scale (-ΔG), 22 kJ mol(-1) in the affinity scale (-ΔH), and 400 cm(-1) in the spectroscopic scale (Δν((OH))) (from benzene to trimethylphosphane oxide and amines). The experimental results in CCl(4) are compared to those computed in the gas phase at various levels of theory. Ninety five percent of the variance of the red shift and 89% of the variance of the intensification of the OH stretching upon hydrogen bonding are explained by gas-phase B3LYP/6-31+G(d,p) calculations. However, this level does not satisfactorily explain the thermodynamic properties. Only 68% of the variance of the methanol affinity (-ΔH) is taken into account. MP2/aug-cc-pVTZ//B3LYP/6-31+G(d,p) affinity calculations raise the explanation to 77% for all HBAs and to 93% when three outliers (Me(2)SO, Me(3)PO, and tetrahydrothiophene) are excluded. Discrepancies are analyzed in terms of experimental errors, calculation approximations, and solvation.