Estimation of the relative stability of chair forms of hexahydropyrimidine

Nonempirical approximations HF/6-31G(d), MP2/6-31G(d)//HF/6-31G(d), PBE/3z, and RI-MP2/λ2 were applied to the calculation of the relative stability and parameters of the conformational equilibrium of three chair forms of the unsubstituted hexahydropyrimidine with different orientation of the N-H bonds. It was established that in the framework of the used methods the most stable are the conformers of diaxial and axial-equatorial chairs. The occupancy of the diequatorial chair form at the room temperature does not exceed 1.7%.     

The photoelectron spectrum and conformation of phenylphosphine and phenylarsine

He(I) and He(II) photoelectron spectra of phenylphosphine and phenylarsine have been investigated and assigned. The rotational barrier of the phosphino group has been investigated at the MP2/6-31G(d,p)//MP2/6-31G(d,p) and HF/6-31G(d,p)//HF/6-31G(d,p) levels of theory, and that of the arsino group at the HF/6-31G(d,p)//6-31G(d,p) levels of theory. The rotational barrier of the two molecules is nearly the same. The energy difference between the two possible conformers of the molecules is low (1.5 kJ/mol at the MP2/6-31G(d,p) level of theory), allowing nearly free rotation about the P-C bond. The photoelectron spectrum cannot be interpreted by considering the most stable rotamer, but all possible conformers should be taken into account. The present interpretation is consistent with the smalln _p - interaction concluded from other investigations. The rotational barrier ofo-phosphinophenol is significantly larger than for phenylphosphine, and the photoelectron spectrum of this compound can be interpreted by considering a single conformer, and no appreciable interaction between the -system of the ring and the phosphorus lone pair.     

Steric and Electronic Structure of C6H5XCF3 Molecules (X = O or S): A Quantum-Chemical Study

The potential functions of internal rotation around the Csp ²-X bond in C₆H₅XCF₃ molecules (X = O or S) were obtained by quantum-chemical calculations in the HF/6-31G(d), MP2(f)/6-31G(d), and B3LYP/6-31G(d) approximations. The calculations were performed in the range of torsion angles (angle between the planes of the benzene ring and Csp ²-X-Csp ³ bonds) from 0° to 90° with a 15° step. The barriers to rotation around the Csp ²-X bonds (kJ mol^{- 1}) were evaluated: for C₆H₅XCF₃, 7.60 (HF), 3.04 (MP2), and 1.04 (B3LYP); for C₆H₅XCF₃, 16.57 (HF), 14.67 (MP2), and 8.73 (B3LYP). The geometries (bond angles and bond lengths), Koopmans ionization potentials, and dipole moments of the molecules were calculated. The hybridization, energy, and population of the lone electron pairs of the heteroatoms, and also the energy of their resonance interaction with antibonding orbitals and the natural atomic charges were evaluated using the NBO approach.     

Conformational and electronic properties of N-methacryloyl-(L)-glutamic acid

Conformers of N-methacryloyl-(L)-glutamic acid (MAGA, C₉H₁₃NO₅) in the gas phase, benzene or methanol environment were searched by the density functional theory (DFT) and time-dependent DFT. Selected conformational isomers based on potential energy surface analysis were studied by considering the structural and electronic properties. B3LYP functional, and HF and MP2 levels were used with the 6-31+G(d,p), 6-311+G(3df,p) and aug-cc-pVDZ basis sets. We explored the effect of medium on the conformational preferences, UV spectra and frontier molecular orbitals of the conformers. The outcomes of this research will be useful for future studies including moieties analogous to MAGA.     

An ab initio study of the structure,dissociation energy,and heat of formation of Na2S

The equilibrium geometries and fundamental frequencies of Na₂S are calculated at HF, MP2(FC, FU), and MP3 with the 6–31G(d) basis set and at HF and MP2(FC, FU) with the 6–31G(d) basis set, respectively. The total energy at MP2(FU)/6–31G(d)-optimized geometry is computed at MP4 with 6–311G(d, p), 6–311 + G(d, p), and 6–311G(2df, p), at QCISD(T)/6–311G(d, p), and at MP2/6–311G(3df, 2p) levels, respectively. The dissociation energy, the atomization energy, and the heat of formation for Na₂S are evaluated using the G1 and G2 models. The calculated results indicated that Na₂S in its ground state was a bent structure (C_2v). Electron correlation corrections on the bending angle are very significant. The equilibrium geometrical parameters are R_e(Na-S) = 2.45 Å and ∠Na-S-Na = 111.13° at the MP2(FU)/6–31G(d) level. The theoretically estimated dissociation energy, total atomization energy, and heat of formation are 67.07, 117.55, and 0.35 kcal mol⁻¹, respectively, at 298.15 K. © 1997 John Wiley & Sons, Inc.     

Quantum-Chemical Study of Electronic and Steric Structure of Vinyltetrazoles: I. 2-Substituted 5-Vinyltetrazoles

The total Mulliken charges on the carbon atoms of the vinyl group, populations of S-trans-(N¹)conformers, and internal rotation energies were calculated ab initio (HF/6-31G**, MP2/6-31G**, and MP2/6-31G**//AM1) for a series of 2R-5-vinyltetrazoles (R = CH₃, C₂H₅, i-C₃H₇, t-C₄H₉, C₆H₅). The calculation results were compared to the available experimental data.     

Theoretical study of the conformational energy hypersurface of cyclotrisarcosyl

The multidimensional Conformational Potential Energy Hypersurface (PEHS) of cyclotrisarcosyl was comprehensively investigated at the DFT (B3LYP/6-31G(d), B3LYP/6-31G(d,p) and B3LYP/6-311++G(d,p)), levels of theory. The equilibrium structures, their relative stability, and the Transition State (TS) structures involved in the conformational interconversion pathways were analyzed. Aug-cc-pVTZ//B3LYP/6-311++G(d,p) and MP2/6-31G(d)//B3LYP/6-311++G(d,p) single point calculations predict a symmetric cis-cis-cis crown conformation as the energetically preferred form for this compound, which is in agreement with the experimental data. The conformational interconversion between the global minimum and the twist form requires 20.88 kcal mol-1 at the MP2/6-31G(d)//B3LYP/6-311++G(d,p) level of theory. Our results allow us to form a concise idea about the internal intricacies of the PEHSs of this cyclic tripeptide, describing the conformations as well as the conformational interconversion processes in this hypersurface. In addition, a comparative analysis between the conformational behaviors of cyclotrisarcosyl with that previously reported for cyclotriglycine was carried out     

On the additivity of basis set effects in some simple fluorine containing systems

The reactions F + H₂ → HF + H, HF → H + F, F → F⁺ + e^? and F + e^? → F^? were used as simple test cases to assess the additivity of basis set effects on reaction energetics computed at the MP4 level. The 6-31G and 6-311G basis sets were augmented with 1, 2, and 3 sets of polarization functions, higher angular momentum polarization functions, and diffuse functions (27 basis sets from 6-31Gd, p) to 6-31 ++ G(3df, 3pd) and likewise for the 6-311G series). For both series substantial nonadditivity was found between diffuse functions on the heavy atom and multiple polarization functions (e.g., 6-31 + G(3d, 3p) vs. 6-31 + G(d, p) and 6-31G(3d, 3p)). For the 6-311G series there is an extra nonadditivity between d functions on hydrogen and multiple polarization functions. Provided that these interactions are taken into account, the remaining basis set effects are additive to within ±0.5 kcal/mol for the reactions considered. Large basis set MP4 calculations can also be estimated to within ±0.5 kcal/mol using MP2 calculations, est. E_MP4(6-31 ++ G(3df, 3pd)) ≈ E_MP4(6-31G(d, p)) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31G(d, p)) or E_MP4(6-31 + G(d, p) + E_MP2(6-31 ++ G(3df, 3pd)) – E_MP2(6-31 + G(d, p)) and likewise for the 6-311G series.     

An Ab Initio Molecular Orbital Theory and Density Functional Theory (DFT) Study of Conformers and Rotamers of 4-Substituted (Methyl, Hydroxymethyl, Methanoyl, Ethanoyl, Cyano, Fluoro, Chloro, Bromo, Acetoxy) Tetrahydro-2H-thiopyran-1,1-dioxides

The structures and energies of axial and equatorial conformers and rotamers of 4-substituted tetrahydro-2H-thiopyran-1,1-dioxides (tetrahydrothiopyran-1,1-dioxides, thiacyclohexane-1,1-dioxides, thiane-1,1-dioxides, and 1,1-dioxothianes; CH₃, CH₂OH, CHO, COCH₃, CN, F, Cl, Br, and OCOCH₃) were calculated using the hybrid density functionals B3LYP, B3P86, and B3PW91, as well as MP2 and the 6-31G(d), 6-31G(2d), 6-31G(3d), 6-31G(d,p), and 6-31+G(d) basis sets. MP2/6-31+G(d)/ /HF/6-31+G(d) [–G° = 1.73 kcal/mol], B3P86/6-31G(d) [–G° = 1.75 kcal/mol], and B3PW91/6-31G(d) [–G° = 1.85 kcal/mol] gave conformational free energy (G°) values at 180 K for 4-methyltetrahydro-2H-thiopyran-1,1-dioxide which were similar to the reported experimental values for methylcyclohexane (–G° = 1.80 kcal/mol), 4-methyltetrahydro-2H-thiopyran (–G° = 1.80 kcal/mol), and other 4-methyl-substituted heterocycles. All levels of theory showed that the conformational preferences of the 4-methanoyl (4-formyl), 4-ethanoyl (4-acetyl), and 4-cyano substituents were small. The HF calculations gave conformational free energy (G°) values for 4-chlorotetrahydro-2H-thiopyran-1,1dioxide which were closer to the experimental value than the MP2 and density functional methods. The best agreement with available experimental data for 4-bromotetrahydro-2H-thiopyran-1,1-dioxide was obtained from the HF/6-31G(2d), HF/6-31G(3d), and B3LYP/6-31G(2d) calculations, and, for 4-acetoxytetrahydro-2H-thiopyran-1,1-dioxide, from the HF/6–31G(3d) calculations. The conformational free energies (G°) and relative energies (E) of the conformers and rotamers have been compared with the correspondingly substituted cyclohexanes and tetrahydro-2H-thiopyrans and are discussed in terms of dipole–dipole (electrostatic) interactions and repulsive nonbonded interactions (steric) in the most stable axial and equatorial conformers. The axial S=O bond lengths are shorter than the equatorial S=O bond lengths and the C2–C3 bond lengths in the substituents with carbon-bonded to the ring are shorter than the C3–C4 and C4–C-5 bond lengths. In contrast, the C2–C3 bond lengths in the 4-halogen and 4-acetoxy substituents are longer than the C3–C4 and C4–C-5 bond lengths.     

Iminopropadienones R–NCCCO and carbon suboxide, C3O2. Theoretical and experimental C NMR spectra

The ¹³C NMR data of five iminopropadienones R–NCCCO as well as carbon suboxide, C₃O₂, have been examined theoretically and experimentally. The best theoretical results were obtained using the GIAO/B3LYP/6-31+G**//MP2/6-31G* level of theory, which reproduces the chemical shifts of the iminopropadienone substituents extremely well while underestimating those of the cumulenic carbons by 5–10 ppm. The computationally faster GIAO/HF/6-31+G**//B3LYP/6-31G* level is also adequate.     

Ab initio HF-SCF study of naphthazarin: Geometries,isomerism, hydrogen bonding,and vibrational spectrum

The structural properties and intramolecular hydrogen bonding of a series of structures of naphthazarin molecule were investigated by ab initio HF-SCF methods. The geometries of theC _2v ,C _2h ,D _2h , andC _s symmetry structures were optimized using split-valence basis sets. MP2/6-31G^*// HF/6-31G single-point energy calculations indicate that theC _2v isomer (5,8-dihydroxy-1,4-naphthoquinone) is the lowest energy structure of the molecule and that theC _2h symmetry one (4,8-dihydroxy-1,5-naphthoquinone), lying 37 kJ/mol above theC _2v form, is the other stable isomer of naphthazarin. At the HF/6-31G level, the intramolecular proton exchange between two equivalentC _2v structures is a two-step process where each proton can be independently transferred through an unsymmetrical potential having a 1,5-quinone intermediate, theC _2h symmetry structure, and two equivalent transition states ofC _s symmetry, with a barrier height equal to 38 kJ/ mol (MP2/6-31G^*//HF/6-31G). The study of naphthazarin molecule is flanked by a theoretical investigation on theC _2v andC _2h isomers of the parent naphthoquinone and dihydroxynaphthalene molecules. The SCF vibrational spectrum of the ground state of naphthazarin, harmonic frequencies, and infrared and Raman band intensities were computed at the HF/6-31G level. The results of the calculations are compared with the matrix isolation FT-IR spectroscopy measurements and with the infrared and Raman spectra of the crystal molecule.     

Quantum-Chemical Study of the Electronic and Steric Structure of Vinyltetrazoles: II. 2-Vinyltetrazoles

The total Mulliken charges on the C and N atoms, populations of the S-trans-(N¹) conformers, and rotation barriers in the molecules of 2-vinyl-5-R-tetrazoles (R = H, CH₃, CH = CH₂, C₆H₅, CH₂Cl, CF₃) were calculated ab initio (HF/6-31G**, MP2/6-31G**). The results were compared with the ¹H and ¹³C NMR data for these compounds.     

Comparative conformational analysis of 1,3-dioxane and 1,3-dithiane

It was shown by ab initio quantum-chemical approximations HF/6-31G(d) and MP2/6-31G(d)//HF/6-31G(d) that the conformational isomerization of 1,3-dioxane and 1,3-dithiane proceeded along common routes. The potential energy surface of both compounds contains six minima including the chair invertomers and enantiomeric flexible forms. They are separated by several potential barriers. It was established by molecular dynamics method that the flexible conformers at heating and keeping at 295–300 K transformed into each other and in the chair conformer.     

Ab initio models for multiple-hydrogen exchange: Comparison of cyclic four- and six-center systems

High-level ab initio calculations {QCISD(T)/6-311 +G**//MP2(fu)/6-31 +G**, with corrections for higher polarization [evaluated at MP2/6-311 +G(3df,2p)] and ΔZPE//MP2(fu)/6-31 +G**, i.e., comparable to Gaussian-2 theory} indicate concerted mechanisms for double- and triple-hydrogen exchange reactions in HF and HCl dimers and trimers, in mixed dimers and trimers containing one NH₃, and in mixed dimers of HF, HCl, and NH₃ with formic acid. All these reactions proceed via cyclic four- or six-center transition structures, the latter being generally more favorable. Calculated activation barriers (ΔH^d? at 0 K, kcal/mol) are 42.3 for (HF)₂, 20.3 for (HF)₃, 41.2 for (HCl)₂, 25.6 for (HCl)₃, 36.0 for NH₃-HF, 10.6 for NH₃(HF)₂, 19.9 for NH₃-HCl, 2.3 for NH₃(HCl)₂, 9.7 for HCO₂H-HF, 7.0 for HCO₂H-HCl, and 11.3, for HCO₂H-NH₃. The barriers are lower for the more ionic systems and when more ion pair character is present. © John Wiley & Sons, Inc.     

Theoretical studies on the conformations of selenamides

Ab initio HF/6-31+G*, MP2/6-31+G*, B3LYP/6-31+G* level calculations have been performed on HSe-NH₂ to estimate the Se-N rotational barriers and N-inversion barriers. Two conformers have been found withsyn andanti arrangement of the NH₂ hydrogens with respect to Se-H bond. The N inversion barriers in selenamide are 1.65, 2.47, 1.93 kcal/mol and the Se-N rotational barriers are 6.58, 6.56 and 6.12 kcal/mol respectively at HF/6-31+G*, MP2/6-31+G* and B3LYP/6-31+G* levels respectively. The n_N →Σ *_Se-H negative hyperconjugation is found to be responsible for the higher rotational barriers.     

AB INITIO MOLECULAR ORBITAL STUDY OF 1,2-DITHIANE AND 1,2,4,5-TETRATHIANE

Ab initio calculations at HF/6-31+G? level of theory for geometry optimization, and MP2/6-31+G?//HF/6-31+G? and B3LYP/6-31+G?//HF/6-31+G? levels for a single-point total energy calculation, are reported for the chair and twist conformations of 1,2-dithiane (1), 3,3,6,6-tetramethyl-1,2-dithiane (2), 1,2,4,5-tetrathiane (3), and 3,3,6,6-tetramethyl-1,2,4,5-tetrathiane (4). The C₂ symmetric chair conformations of 1 and 2 are calculated to be 21.9 and 8.6 kJ mol^?1 more stable than the corresponding twist forms. The calculated energy barriers for chair-to-twist processes in 1 and 2 are 56.3 and 72.8 kJ mol^?1, respectively. The C_2h symmetric chair conformation of 3 is 10.7 kJ mol^?1 more stable than the twist form. Interconversion of these forms takes place via a C₂ symmetric transition state, which is 67.5 kJ mol^?1 less stable than 3-Chair. The D₂ symmetric twist-boat conformation of 4 is calculated to be 4.0 kJ mol^?1 more stable than the C_2h symmetric chair form. The calculated strain energy for twist to chair process is 61.1 kJ mol^?1.     

Theoretical Investigations on Fluorine-Substituted Ethylene Dications C2HnF4-n 2+(n = 0–4)

The optimized geometries and energies of fluorine-substituted ethylene dications C₂H_nF_4-n ²⁺ (n = 0–4) have been investigated by means of ab initio methods. At the MP3/6-31G**//6-31G* + zero-point energy level of theory, the results predict that C₂F₄²⁺ and C₂HF₃²⁺ are planar, while C₂H₄²⁺, C₂H₃F²⁺ and 1,1—C₂H₂F₂²⁺ prefer a perpendicular geometry. For 1,2—C₂H₂F₂²⁺ an energy difference of only 0.3 kcal/mol is found between the (trans) planar and perpendicular structure. The stabilizations attributed to hyperconjugation, fluorine lone-pair donation, and (C? F) double-bond conjugation are discussed. A comparison is made for the C? C and C? F stretching frequencies determined at 6-31G*//6-31G* between the neutral and dicationic species. The theoretically determined ionization energies for the vertical process N⁺ → N²⁺ at the MP3/6-31G*//3-21G level are compared with experimental Q_min values.     

Structures and stability of N9, N9 − and N9 + clusters

 Ab initio molecular orbital calculations for N₉, N⁻ ₉ and N⁺ ₉ isomers were carried out at the HF/ 6-31G*, B3PW91/6-31G*, B3LYP/6-31G* and MP2/ 6-31G* levels of theory. Stable equilibrium geometric structures were determined by harmonic vibrational frequency analyses at the HF/6-31G*, B3PW91/6-31G* and B3LYP/6-31G* levels of theory. The most stable free-radical N₉ cluster is structure 1 with C _{2 v} symmetry and that of anion N⁻ ₉ is structure 3 with C _s symmetry. Only one stable structure of the N⁺ ₉ cation with C _{2 v} symmetry was predicted. Their potential application as high-energy-density materials has been examined. Received: 15 June 1999 / Accepted: 11 October 1999 / Published online: 14 March 2000     

Is tetrahedral H42+ a minimum? Anomalous behavior of popular basis sets with the standard p exponents on hydrogen

The nature of the tetrahedral H₄²⁺ stationary point (minimum or triply degenerate saddle) depends remarkably upon the theoretical level employed. Harmonic vibrational analyses with, e.g., the 6-31G** (and 6-31 + +G**) and Dunning's [4s2p1d;2s1p] [D95(d,p)] basis sets using the standard p exponent suggest (erroneously) that the T_d geometry is a minimum at both the HF and MP2 levels. This is not the case at definitive higher levels. The C₃H₄²⁺ structure with an apical H is another example of the failure of the calculations with the 6-31G**, 6-311G**, and D95(d,p) basis sets. Even at MP2/6-31G** and MP2/ cc-pVDZ levels, the C_3v structure has no negative eigenvalues of the Hessian. Actually, this form is a second-order saddle point as shown by the MP2/6-31G** calculation with the optimized exponent. The D_4h methane dication structure is also an example of the misleading performance of the 6-31G** basis set. In all these cases, energy-optimized hydrogen p exponents give the correct results, i.e., those found with more extended treatments. Optimized values of the hydrogen polarization function exponents eliminate these defects in 6-31G** calculations. Species with higher coordinate hydrogens may also be calculated reliably by using more than one set of p functions on hydrogen [e.g., the 6-31G(d,2p) basis set]. Not all cases are critical. A survey of examples, also including some boron compounds, provides calibration. © 1993 John Wiley & Sons, Inc.     

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.