Infrared matrix isolation and ab initio quantum mechanical studies of purine and adenine

Results are presented from ab initio SCF (3-21G) calculations for the geometries of the N(9)H and N(7)H tautomers of purine and adenine and vibrational spectra (wavenumbers and intensities) of the N(9)H forms. All these results are compared with available geometries from crystallographic studies and with reported infrared spectra of the molecules isolated in inert gas low-temperature matrices. The N(9)⇌N(7)H tautomerism of the molecules in question is also briefly discussed.     

Infrared matrix isolation and ab initio studies of 4-oxopyrimidine; separation of the spectra of tautomers based on the phototautomeric effect

The i.r. absorption spectra of two simultaneously existing tautomers]—4-oxopyrimidine and 4- hydroxypyrimidine—isolated in several inert gas matrices are reported, and interpreted with the aid of ab initio 3–21G normal coordinate calculations. The phototautomeric effect, i.e. the change of the tautomeric form oxo → hydroxy upon u.v. irradiation is observed and used to separate the spectra of both tautomers.     

Theoretical study of “protonated pyruvate”: A methylhydroxycarbene—carbon dioxide complex—implications for the decarboxylation of pyruvic acid

Two conformers of protonated pyruvate, CH₃C⁺(OH)COO, with the OH group either trans or cis to the methyl group and the carboxylate group in the C? C? C plane have been studied using the ab initio SCF/3-21G method, as well as by some semiempirical AM1 calculations. Both ab initio SCF and AM1 curves for the potential energy as a function of the C? COO distance exhibit a minimum corresponding to a complex of methylhydroxycarbene, CH₃COH, associated with carbon dioxide, but only the AM1 curves predict an inner minimum corresponding to a covalently bonded protonated pyruvate molecule with a C? COO distance of 1.6–1.7 Å. The two models also disagree on the dissociation pathway for pyruvic acid, with the AM1 calculations predicting formation of acetyl and HOCO radicals while the ab initio method predicts dissociation into methylhydroxycarbene and carbon dioxide following an initial intramolecular proton transfer. The weakly bound complexes of methylhydroxycarbene and carbon dioxide have been studied in some detail using ab initio SCF and MP2 methods in conjunction with 6-311G** basis sets, obtaining equilibrium geometries and vibrational frequencies. In addition, the lactone-type isomer of protonated pyruvate, which contains a C? C? O ring, was also studied. The conclusions of these calculations are consistent with those from earlier work using the smaller 3-21G basis set. The most stable complex is predicted to occur between trans-methylhydroxycarbene and carbon dioxide where substantial stabilization is provided by an OH ? OC hydrogen bond. © 1993 John Wiley & Sons, Inc.     

Conformational stability of trivinylborane from vibrational spectra and ab initio calculations

The conformational stability, barriers to internal rotation and vibrational frequencies of trivinylborane have been determined from the vibrational spectra and ab initio calculations. The ab initio calculations have been carried out utilizing the RHF/3-21G, RHF/6-31G*, and MP2/6-31G* basis sets and support the vibrational data that there are two stable conformations in the fluid phases separated by a relatively small energy difference. One of the conformations is a near-planar form which has the three vinyl groups twisted out of the BC₃ plane and belongs to the C₃ point group. The other conformer has a non-planar structure and belongs to the C₁ point group. These and other calculated results are compared to the corresponding quantities obtained from the experiment.     

Conformation and tautomerism of methoxy‐substituted 4‐phenyl‐4‐thiazoline‐2‐thiones: a combined crystallographic and ab initio investigation

4‐Phenyl‐4‐thiazoline‐2‐thiol is an active pharmaceutical compound, one of whose activities is as a human indolenamine dioxygenase inhibitor. It has been shown recently that in both the solid state and the gas phase, the thiazolinethione tautomer should be preferred. As part of both research on this lead compound and a medicinal chemistry program, a series of substituted arylthiazolinethiones have been synthesized. The molecular conformations and tautomerism of 4‐(2‐methoxyphenyl)‐4‐thiazoline‐2‐thione and 4‐(4‐methoxyphenyl)‐4‐thiazoline‐2‐thione, both C₁₀H₉NOS₂, are reported and compared with the geometry deduced from ab initio calculations [PBE/6‐311G(d,p)]. Both the crystal structure analyses and the calculations establish the thione tautomer for the two substituted arylthiazolinethiones. In the crystal structure of the 2‐methoxyphenyl regioisomer, the thiazolinethione unit was disordered over two conformations. Both isomers exhibit similar hydrogen‐bond patterns [R₂²(8) motif] and form dimers. The crystal packing is further reinforced by short S…S interactions in the 2‐methoxyphenyl isomer. The conformations of the two regioisomers correspond to stable geometries calculated from an ab initio energy‐relaxed scan.     

Raman spectra,conformational stability,structural parameters,vibrational assignment,and ab initio calculations for epifluorohydrin

The Raman spectra (3200–100 cm⁻¹) of epifluorohydrin, OCH₂CH(CH₂F), in variable solvents, as well as that of the gas have been recorded and several of the bands due to the two less stable conformers have been identified. The variable solvent studies were inconclusive on the relative conformer stabilities. The conformational energy differences and optimized geometries for all three conformers have been obtained from ab initio calculations with the 3–21G, 4–31G and 6–31G* basis sets. A normal coordinate analysis has also been performed for each conformer with a force field determined from the 3–21G basis set. Assignment of the vibrational fundamentals observed in the Raman spectra of the fluid phases is proposed based on the normal coordinate calculations. In the liquid phase, one of the conformers with a large dipole moment predominates and it appears to be the gauche-I form which is the only one found in the solid. Utilizing the three rotational constants previously reported for each conformer, along with restricted relative distances for several of the structural parameters among the conformers from ab initio calculations, r₀ structural parameters for the heavy atoms have been determined.     

Vibrational and theoretical studies on CF3SNH2: A compound containing an SNH2 group

The structure and conformational behaviour of CF₃SNH₂ have been studied by ab initio calculations (HF/3-21G* (with additional d-functions on nitrogen, α_d: 0.8) and HF/6-31G**) and vibrational spectroscopy. Both theoretical methods predict a trans form (CS bond in trans position with respect to the N lone pair, C_s symmetry) to be more stable than the corresponding cis conformer by 1.5 and 0.7 kcal/mol, respectively. The calculated barriers to inversion are 2.7 (HF/3-21G*) and 1.6 kcal/mol (HF/6-31G**). A depolarized band observed in the Raman spectra of the liquid is in accordance with the existence of a plane of symmetry in the molecule. A normal coordinate analysis was performed for the A′ modes of vibration.     

Ab initio and molecular mechanics (MM2 and MM3) calculations of nonconjugated positively charged nitrogen-containing compounds

High-level ab initio calculations have been performed on N-methyl-N-methyleneammonium and related compounds to obtain accurate rotational barriers, structures, and vibrational frequencies. The 6-31G** basis set has been utilized at the Hartree-Fock level of theory for these calculations because little experimental data are available. The MM2(91) and MM3(94) force fields have been parameterized to include these nonconjugated charged nitrogen-containing compounds. Molecular mechanics geometries and vibrational frequencies compare well with the ab initio results. © 1995 John Wiley & Sons, Inc.     

Vibrational absorption intensities in chemical analysis—VIII. Dimethyl sulphide

The structure and electric dipole moment of dimethyl sulphide have been calculated with a variety of basis sets. It is shown that the force fields computed with the 3-21G, 3-21G and 6-311G** bases are very similar after scaling. The vibrational assignments are discussed. The force field is compared with that of trimethyl phosphine. There is an excellent accord for the CH stretching modes, but the ab initio calculations indicate that the X-C/X-C (X = S,P) interaction terms have opposite signs. This appears to be incompatible with the experimental data. Effective atomic charges and atomic polar tensors (APT) are computed with the above bases and also with bond centred sp functions added to the 3-21G* basis. It is shown that the APT are much more sensitive to the basis set than are the effective atomic charges. The effects of the lone pair electrons on the bond electrical properties are much less than with oxy-ethers, though the overall perturbation of the sulphur atom on the intensities is greater. A number of inadequacies of the SCF calculations are identified.     

Spectra and structure of small ring compounds. LXII. Conformational stability,structural parameters,ab initio calculations and vibrational assignment of cyclopropylmethyl ketone

The conformational stability, barriers to internal rotation, and fundamental vibrational frequencies of cyclopropylmethyl ketone, c-C₃H₅C(CH₃)O, have been obtained from Hartree—Fock ab initio calculations with the RHF/3-21G and RHF/6-31G* basis sets, as well as the 6-31G* basis set with electron correlation at the MP2 level, and the results are compared to those obtained from experiment. The data are consistent with the predominant rotamer having the cis conformation (carbonyl bond cis to the ring). A second form, having a “near” trans structure, is calculated to have a larger total dipole moment than the cis form, which accounts for its increased abundance in the liquid compared to that in the gas. A complete vibrational assignment is proposed based on experimental data and normal coordinate results from the ab initio calculations. The asymmetric torsional barrier has been calculated to be approximately 2000 cm⁻¹ and this result along with others is compared to the corresponding data obtained from both experiment and theory for the cyclopropylcarbonyl halides.     

Structure and tautomerism of 4-substituted 3(5)-aminopyrazoles in solution and in the solid state: NMR study and Ab initio calculations

Annular tautomerism of 3(5)-aminopyrazoles containing a cyano, thiocyanato, or aryl substituent in the 4-position has been studied by ¹H and ¹³C NMR in solution, cross-polarization and magic-angle spinning ¹³C NMR in the solid state, and ab initio quantum chemical calculations (B3LYP/6-31G**). The title compounds in the solid state exist as 3-amino tautomers. A rare case of slow (on the NMR time scale) annular prototropic tautomerism has been observed in DMSO-d ₆: signals of particular tautomers (3- and 5-aminopyrazoles) have been detected in the NMR spectra. 4-Cyano and 4-thiocyanato derivatives exist preferentially as 5-amino tautomers, whereas 4-methoxy analog is represented mainly by the 3-amino tautomers. Ab initio calculations (B3LYP/6-31G**) for the gas phase and DMSO solution (in terms of the polarizable continuum model) have shown increase of the relative stability of more polar 5-amino tautomer in going to DMSO.     

Comparison of the use of the MNDO and MINDO/3 methods with ab initio methods to study tautomerism in histamine, 2- and 4-methylhistamine

The semiempirical MNDO and MINDO/3 methods are used to study the various tautomeric forms of histamine, 2-methylhistamine, and 4-methylhistamine. Comparisons of the optimized structures and tautomerization energies are made with values obtained from ab initio Hartree-Fock calculations using the 3-21G and STO-3G basis sets. Based on these results and previous comparisons of STO-3G results with x-ray structures, the present results indicate that while there are some differences in the values of the structural parameters, the changes in structure upon tautomerization and/or protonation are very similar. Further analysis of the MNDO and MINDO/3 structures by means of their utilization in 3-21G and STO-3G calculations indicates that either of these semiempirical methods provides reliable values for the structural parameters. Both methods give good qualitative agreement with the ab initio calculations for the relative energies of the various tautomers in the three compounds. In these studies the MNDO method appears to give better quantitative agreement with the 3-21G and STO-3G results than the MINDO/3 method.     

An investigation of the relative stabilities of the isomers of CF2N2: Comparison of ab initio and MNDO calculations

Ab initio SCF calculations at the HF/3-21G level and semi-empirical MNDO calculations have been used to locate the stationary points on the CF₂N₂ energy surface. Perfluorodiazomethane is predicted to be most stable isomer, but perfluorodiazirine is predicted to lie only ca 41 kJ higher in energy at the SCF level. There are significant differences between the ab initio and MNDO results for the ordering of some of the isomers. Frequency calculations give results in good agreement with the limited experimental data on these molecules.     

An SCF and CI Study of the 1,3 Shift in the HX?CHY⇌XCH?YH Isoelectronic Series: X,YCH2, NH,and O

Ab initio molecular orbital calculations at SCF level with the 3-21G, 6-31G, and 6-31G** basis sets and CI level with the 6-31G basis set have been carried out for an isoelectronic series HX? CH?Y and X?CH? YH, where X, Y can be CH₂, NH, and O. Optimized structures (3-21G and 6-31G**) for both tautomers and the 1,3 hydrogen shift transition states are reported. The relative stabilities of the isomers and the barriers of the 1,3 shift are discussed in terms of proton affinities and bond orders. It is shown that both the relative stabilities of the tautomers and the relative barrier heights can be explained qualitatively using simple proton affinity arguments and that the barrier heights are quantitatively related to bond orders.     

Conformational studies of n-propylamine by combined ab initio MO calculations and Raman spectroscopy

The results of ab initio SCF-MO calculations performed with a 3-21G(N*) basis set, for fully optimized geometries of five conformations of n-propylamine, are presented. The calculated relative order of total energies for these conformers is TT≈GG′>TG>GT>GG. At 300 K, the Boltzmann distribution of populations is 18, 37, 20, 19 and 7%, respectively.Raman spectra of n-propylamine and n-propylamine-N-d₂ in the liquid phase exhibit a number of bands whose temperature-dependent intensities clearly suggest the occurrence of different conformers in simultaneous equilibria. Deuteration of the amine group originates pairs of Raman bands at 428 and 440 cm⁻¹ and at 863 and 885 cm⁻¹. The bands at 428 and 885 cm⁻¹ are favoured by reduction of temperature. Normal coordinate calculations permit the assignment of the Raman and i.r. spectra in good agreement with experimental evidence. Among the five possible conformers of n-propylamine, it is possible to detect the presence of at least three conformations in the liquid phase, corresponding to the skeletal trans (TT and GT) and at least one of the skeletal gauche (TG, GG or GG′) forms. In the solid phase, only the bands ascribed to the TT form were observed.The ab initio results for the isolated molecule show that the all-trans conformation, TT, and the conformation GG′ have the smallest energies. On the other hand, the vibrational results for the liquid and solid phases indicate that the all-trans conformation, TT, is the more populated form. In addition, this conformer presents the highest calculated dipole moment, in good agreement with the liquid phase Raman spectroscopic results which point out that this conformation is favoured by polar solvents. Intermolecular interactions operating in the liquid n-propylamine, possibly of the hydrogen bonding type, are responsible for altering the relative order of conformational stability as predicted by the ab initio SCF-MO results for the isolated molecule.     

Improvement of the hydrogen bonding correction to MNDO for calculations of biochemical interest

Recent suggestions for correcting H? Acceptor interactions within MNDO, together with results of crystallographic analysis, were used to modify this SCF semiempirical calculation for multiple hydrogen bonded associations. Thermodynamic profiles for model systems of biochemical interest such as H₂O? H₂O, hydration of neutral and charged molecules, dimerizations and proton transfers show the advantages of this method. Its treatment of charges, bonding, geometries, energetics and vibrational frequencies is shown to be comparable to ab initio calculations with various basis sets. However, basic MNDO deficiencies and criteria applied for H-bonding result in some too high barriers for proton transfers.     

Ab initio study of vibrational frequencies,force field and normal coordinate analysis of the cis- and trans-isomets of nitrosomethanol

In this paper, the equilibrium geometries of two isomers of the newly found compound — nitrosomethanol—have been optimized by ab initio SCF MO method with 3–21G basic set by gradient technique. And the second derivatives of potential energy (i. e. the force constant matrix elements) have been calculated analytically. Hence the entire force fields of the two isomers of nitrosomethanol have been obtained theoretically. The theoretical vibrational frequencies and the corresponding normal modes were obtained and compared with the experimental values, and the structures of two isomeric forms of nitrosomethanol are established.     

Study of hydrogen bonding interactions relevant to biomolecular structure and function

Ab initio molecular orbital calculations were used to study hydrogen bonding interactions and interatomic distances of a number of hydrogen bonded complexes that are germane to biomolecular structure and function. The calculations were carried out at the STO-3G, 3-21G, 6-31G*, and MP2/6-31G* levels (geometries were fully optimized at each level). For anionic species, 6-31 + G* and MP2/6-31 + G* were also used. In some cases, more sophisticated calculations were also carried out. Whenever possible, the corresponding enthalpy, entropy, and free energy of complexation were calculated. The agreement with the limited quantity of experimental data is good. For comparison, we also carried out semiempirical molecular orbital calculations. In general, AM1 and PM3 give lower interaction enthalpies than the best ab initio results. With regard to structural results, AM1 tends to favor bifurcated structures for O? H-O and N? HO types of hydrogen bonds, but not for hydrogen bonds involving O-H? S and S-H? O, where the usual hydrogen bond patterns are observed. Overall, AM1 geometries are in general in poor agreement with ab initio structural results. On the other hand, PM3 gives geometries similar to the ab initio ones. Hence, from the structural point of view PM3 does show some improvement over AM1. Finally, insights into the formation of cyclic or open formate–water hydrogen bonded complexes are presented. © 1992 by John Wiley & Sons, Inc.     

Experimental matrix isolation and theoretical ab initio HF/6-31G(d,p) studies of infrared spectra of purine,adenine and 2-chloroadenine

The results of ab initio calculations at the Hartree—Fock level with the 6-31G(d, p) basis set are reported for the harmonic vibrational infrared spectra of purine, adenine and 2-chloroadenine. These were compared with the experimental IR spectra of the compounds isolated in low-temperature matrices. The experimental spectra of 2-chloroadenine are reported for the first time. In the case of purine this comparison resulted in the complete assignment of the experimental IR spectrum. The proposed assignments of both the adenines studied are still incomplete.     

Comparative study of the molecular electrostatic potential obtained from different wavefunctions. Reliability of the semiempirical MNDO wavefunction

A systematic analysis of the molecular electrostatic potential (MEP) is presented. This study has been performed with a twofold purpose: first, to study the MEP dependence with regard to the quality of the basis set used to compute the ab initio SCF wavefunction and second, to develop and to assess a new strategy for computing isoelectrostatic potential maps using the semiempirical MNDO wavefunction. The only differences between this procedure and the ab initio SCF MEP computation lie in the freezing of the inner electrons and in the origin of the first-order density matrix. The statistical analysis of MEPs computed for a large number of molecules from MNDO wavefunction and ab initio SCF wavefunctions obtained using STO-3G, 4-31G, 6-31G, 4-31G*, 6-31G*, and 6-31G** basis sets points out the ability of any wavefunction to reproduce the general topological characteristics of the MEP surfaces. Nevertheless, split-valence basis sets including polarization functions are necessary to obtain accurate MEP minimum energy values. MNDO wavefunction tends to overestimate the MEP minima depth by a constant factor and shows an excellent ability to reflect the relative variation of MEP minima energies derived from a rather sophisticated (6-31G*) basis set, lacking of the shortcomings detected in the semiempirical CNDO approximation.