Theoretical studies of the benzene oxide—oxepin valence tautomerism

We have calculated the geometry and energy of the valence tautomers benzene oxide and oxepin using the semiempirical AM1 model and the 6–31G and 6–31G^* basis sets utilizing full geometry optimization. In the oxide the folding angle, the angle between the epoxide ring and the adjacent plane containing four carbon atoms, is about 106°. The carbon skeleton is almost planar, the folding angle, the angle between the two four-carbon atom planes being about 175°. In contrast, oxepin is found to have a marked boat-shaped structure with the corresponding and angles about 137° and 159°, respectively. The AM1, 6–31G, and 6–31G^* calculations give –11.4, –10.8, and –2.9 kcal mol^–1 for the energy change that accompanies the valence tautomerism, oxide-oxepin, compared to an experimental value of about +0.3 kcal mol^–1. Single point calculations of the energies at the 6–31 G^* geometry using Møller-Plesset perturbation theory to second order (MP2/6–31 G^*) and third order (MP3/6–31G^*) give E _T =+3.3 and +0.8 kcal mol^–1. The values for the energy change in the transfer of epoxide oxygen from ethylene oxide to benzene using AM1, 6–31G, and 6–31G^* are in good agreement, viz., +31.1, +34.5, and +33.6 kcal mol^–1, respectively. A large positive energy change is to be expected in view of the loss of benzene aromaticity.     

Double H-bridged and single H-bridged diboryl radicals

The structures of B₂H₅·, B₂H₅CO·, and B₂H₅N₂· radicals are investigated using the 6–31G^* basis set. Both double H-bridged and single H-bridged isomers are found to be local minima on the potential energy surface. The effects of electron correlation are taken into account using single point MP4/6–31G^* calculations and, for the diboryl radicals, complete MP3/6–31G^* optimizations. In all cases the single H-bridged isomers are found to be more stable than the corresponding double H-bridged isomers.The transition state for the double H-bridged to single H-bridged B₂H₅· isomerization reaction is calculated to be 2.54 kcal mol^–1 above the double H-bridged radical at the MP4SDTQ/6-31G^*//UHF/ 6–31 G^* level when corrected for zero point energy. Barrier tunneling increased the reaction rate by a factor of 2.5–3.0, strongly suggesting the system is fluxional at this temperature.The addition of CO and N₂ to the diboryl radicals leads to relocation of the unpaired electron and rehybridization of the C and N atoms adjacent to the boron atoms. The isomers of B₂H₅CO· and B₂H₅N₂· are different and should be distinguishable experimentally. While the CO moiety is bound to the diboryl radicals isomers by over 19 kcal mol^–1, no binding energy is evident for N₂.     

Reexamination of the C2H4F potential surface. the status of the classical 2-fluoroethyl cation: a local minimum or transition structure?

According to ab initio molecular orbital calculations carried out with full geometry optimization at the MP2/6–31G^** level, the classical 2-fluoroethyl cation, FCH₂CH₂⁺, is a transition structure for H-scrambling in CH₃CHF⁺. Single point MP4/6–31G^** calculations at the optimized geometries predict the cyclic ethylene fluoronium ion to lie 24.2 kcal mol⁻¹ above CH₃CHF⁺ and 5.4 kcal mol⁻¹ below the 2-fluoroethyl cation. ΔG‡ for ring opening of the cyclic fluoronium ion at -60° is estimated to be ca 15 kcal mol⁻¹. This barrier is largely attributable to the powerful negative fluorine hyperconjugation in the transition state as described by Hoffmann and coworkers. When electron correlation effects are ignored a qualitatively different potential surface is obtained on which the 2-fluoroethyl cation is calculated to be a local minimum separated from the stable 1-fluoroethyl cation by an H-bridged transition state.     

C36 fullerenes and quasi-fullerenes: computational search through 598 cages

C₃₆ is computed at the SAM1 level and partially also at the HF/4-31G and B3LYP/6-31G^* levels. Altogether 598 cages are generated by a topological Stone–Wales treatment. Three cages contribute by more than 10% to the high-temperature equilibrium isomeric mixture – two conventional fullerenes with D_2d and C_2v symmetries and a C_s quasi-fullerene containing one four-membered ring.     

Isomerism in OBe3F3 + cation: anab initio study

Ab initio MP2/6-31G^*//HF/6-31G^*+ZPE(HF/6-31G^*) calculations of the potential energy surface in the vicinity of stationary points and the pathways of intramolecular rearrangements between low-lying structures of the OBe₃F₃ ⁺ cation detected in the mass spectra of μ₄-Be₄O(CF₃COO)₆ were carried out. Ten stable isomers with di- and tricoordinate oxygen atoms were localized. The relative energies of six structures lie in the range 0–8 kcal mol⁻¹ and those of the remaining four structures lie in the range 20–40 kcal mol⁻¹. Two most favorable isomers, aC _2v isomer with a dicoordinate oxygen atom, planar six-membered cycle, and one terminal fluorine atom and a pyramidalC _3v isomer with a tricoordinate oxygen atom and three bridging fluorine atoms, are almost degenerate in energy. The barriers to rearrangements with the breaking of one fluorine bridge are no higher than 4 kcal mol⁻¹, except for the pyramidalC _3v isomer (∼16 kcal mol⁻¹). On the contrary, rearrangements with the breaking of the O−Be bond occur with overcoming of a high energy barrier (∼24 kcal mol⁻¹). A planarD _3h isomer with a tricoordinate oxygen atom and linear O−Be−H fragments was found to be the most favorable for the OBe₃H₃ ⁺ cation, a hydride analog of the OBe₃F₃ ⁺ ion; the energies of the remaining five isomers are more than 25 kcal mol⁻¹ higher. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 3, pp. 420–430, March, 1999.     

An ab initio Molecular Orbital Study of the Conformational Properties of all-( Z )-Cyclododeca-1,4,7,10-tetraene

 Ab initio HF/6-31G^* and MP2/6-31G^*//HF/6-31G^* methods were used to calculate the structure optimization and conformational interconversion pathways for all-(Z )-cyclododeca-1,4,7,10-tetraene. This compound adopts the symmetrical crown (C _4v) conformation. Ring inversion takes place via symmetrical intermediates, such as boat-chair (BC, C _s) and twist (C _2h) conformers and requires about 22.3 kJ · mol⁻¹. The calculated strain energies for BC and twist conformers are 5.9 and 13.5 kJ · mol⁻¹. The results of semiempirical AM1 calculations for structural parameters and relative energies of the important geometries of the title compound are in good agreement with the results of ab initio methods.     

Minimal basis set MINI-1 — powerful tool for calculating of molecular interactions. I. Neutral complexes

The Huzinaga MINI-1 minimal basis set was applied to the following complexes: H₃N...HF, H₂O...HF, (HF)₂, (H₂O)₂, HCN...HF, OCO...HF, HF...HNH₂, (HCl)₂, CH₄...OH₂ and (H₃Si)₂O...H₂O. The optimized geometries and stabilization energies agree well with the 6–31G^* values, while the computing time is reduced considerably. The MINI-1 basis set was further tested for the stacking structure of the ethylene dimer and the cyclic structure of the formamide...formamidine complex, where fair agreement with the 4–31G results was obtained. The normal vibrational frequencies calculated for some complexes are comparable to those evaluated at the 4–31G (6–31G) level.     

FT-IR, FT-Raman vibrational spectra and molecular structure investigation of 2-chloro-4-methylaniline: A combined experimental and theoretical study

In this work, the experimental and theoretical vibrational spectra of 2-chloro-4-methylaniline (2Cl4MA, C₇H₈NCl) were studied. FT-IR and FT-Raman spectra of 2Cl4MA in the liquid phase have been recorded in the region 4000–400 cm⁻¹ and 3500–50 cm⁻¹, respectively. The structural and spectroscopic data of the molecule in the ground state have been calculated by using Hartree-Fock (HF) and density functional method (B3LYP) with the 6-31G(d), 6-31G(d,p), 6-31+G(d,p), 6-31++G(d,p) and 6-311G(d), 6-311G(d,p), 6-311+G(d,p), 6-311++G(d,p) basis sets. The vibrational frequencies have been calculated and scaled values have been compared with experimental FT-IR and FT-Raman spectra. The observed and calculated frequencies are found to be in good agreement. The complete assignments were performed on the basis of the total energy distribution (TED) of the vibrational modes, calculated with scaled quantum mechanics (SQM) method. The DFT-B3LYP/6-311++G(d,p) calculations have been found more reliable than the ab initio HF/6-311++G(d,p) calculations for the vibrational study of 2Cl4MA. The optimized geometric parameters (bond lengths and bond angles) were compared with experimental values of aniline and p-methylaniline molecules.     

Synthesis and vibrational study of platinum(II) and palladium(II) complexes of glyoxilic acid oxime

New platinum(II) and palladium(II) complexes of glyoxilic acid oxime (gao) have been prepared and characterised by infrared (4000–150 cm⁻¹) and Raman (4000–200 cm⁻¹) spectra. The gao acts as bidentate ligand bonding through the oxime nitrogen and carboxyl oxygen atoms to form neutral bis-chelate square-planar complexes. The lowest energy conformer of the gao ligand (ectt) was selected among 16 theoretically possible conformers on the basis of ab initio calculations at HF/3-21G*, HF/6-31G* and HF/6-311** levels of the theory from which structural parameters and conformational stabilities have been obtained. A complete vibrational assignment of the gao was performed for the lowest energy ectt conformer on the basis of ab initio optimised parameters and normal coordinate analysis calculations (PED). NCA calculations of the complexes studied were also performed.     

Investigation of the effects of ionic hydrogen bonds (COH–C and N–HOC) in crystalline dl-proline by ab initio and DFT calculated NQR parameters

In this paper, we have calculated the nuclear quadrupole resonance (NQR) parameters of the quadrupole nuclei involved in the hydrogen bonds (CO⁻H–C and ⁺N–H⁻OC) in the monomer and pentameric cluster of dl-proline by HF and B3LYP methods and basis sets of 6-311+G^* and 6-311++G^**. These computations are performed on the basis of X-ray diffraction structural data of dl-proline. The results indicate that the calculations including hydrogen-bonding (HB) interactions (in pentamer) are in better agreement with the experimental data than those in which these interactions are neglected (in monomer). The quantum chemical calculations show that the intermolecular hydrogen-bonding interactions play an important role in determination of the NQR parameters of ¹⁴N, ²H of group and ¹⁷O.     

The aromatic diboracyclopropenyl (diboriranyl) anion; CB2H3 −: An ab initio study

The most stable structure of CB₂H₃ ^–, as established computationally, is the aromatic diboracyclopropenyl (diboriranyl) anion (5), while open-chainC _2v, isomer H₂BCBH (7) is only 3 kcal/mol higher in energy at the QCISD(T)/6-311 +G^**//MP2/6-31+G^*+ZPE (HF/6-31 +G^*). The 47-kcal/mol barrier between cyclic,5, and open-chain,7, structures suggests that both of them may be observed. The aromatic stabilization energy of the diboriranyl anion (18 kcal/mol) is half the value in the isoelectronic cyclopropenium ion, C₃H₃ ⁺. The computed, by IGLO method (5a), and experimental (6a) chemical shifts,(¹³C) and(¹¹B), agree within 4 ppm range. The theoretical vibrational frequencies of the most stable isomers,5 and7, are presented for experimental verification of these species.     

Molecular structure and vibrational assignment of melaminium phthalate by density functional theory (DFT) and ab initio Hartree-Fock (HF) calculations

The molecular geometry, the normal mode frequencies and corresponding vibrational assignment of melaminium phthalate (C₃H₇N₆⁺·C₈H₅O₄⁻) in the ground state were performed by HF and B3LYP levels of theory using the 6-31G(d) basis set. The optimized bond length numbers with bond angles are in good agreement with the X-ray data. The vibrational spectra of melaminium phthalate which is calculated by HF and B3LYP methods, reproduces vibrational wave numbers with an accuracy which allows reliable vibrational assignments. The title compound has been studied in the 4000–100 cm⁻¹ region where the theoretical evaluation and assignment of all observed bands were made.     

Effect of intramolecular hydrogen bonding on the empty level structure of gas-phase salicylaldehyde:electron affinities and anion yield upon electron impact

The gas-phase empty level structures of γ-butyrolactone, α-hydroxy-γ-butyrolactone, benzaldehyde, 3-hydroxybenzaldehyde and 2-hydroxybenzaldehyde (salicylaldehyde) are studied by means of electron transmission spectroscopy (ETS) and ab initio 6-31G^** calculations. The most stable conformer of salicylaldehyde is characterised by intramolecular hydrogen bonding which causes relevant geometrical variations and, in particular, a lengthening of the C=O double bond. The empty π^* molecular orbitals of salicylaldehyde with large contribution from the carbonyl group are significantly stabilised with respect to their counterparts in benzaldehyde and in the 3-hydroxy derivative. Dissociative electron attachment spectroscopy (DAS) did not reveal formation of negative fragments at low energy, that is, the presence of dissociative decay channels for the π^* resonances observed in the ET spectra. The total anion current measured at the walls of the collision chamber (close to the electron beam), as a function of the incident electron energy, displays signals peaking at zero energy and at about 0.8 eV in benzaldehyde and in the two hydroxy derivatives. However, the only signal detected using a quadrupole mass filter (which requires a lifetime ≥10⁻⁶ s) is a zero energy peak in salicylaldehyde, due to the molecular anion.     

An ab initio Molecular Orbital Study of the Conformational Properties of all-(Z)-Cyclododeca-1,4,7,10-tetraene

Summary.  Ab initio HF/6-31G^* and MP2/6-31G^*//HF/6-31G^* methods were used to calculate the structure optimization and conformational interconversion pathways for all-(Z )-cyclododeca-1,4,7,10-tetraene. This compound adopts the symmetrical crown (C _4v) conformation. Ring inversion takes place via symmetrical intermediates, such as boat-chair (BC, C _s) and twist (C _2h) conformers and requires about 22.3 kJ · mol⁻¹. The calculated strain energies for BC and twist conformers are 5.9 and 13.5 kJ · mol⁻¹. The results of semiempirical AM1 calculations for structural parameters and relative energies of the important geometries of the title compound are in good agreement with the results of ab initio methods. Received July 9, 2001. Accepted September 26, 2001     

A molecular orbital study of nitrogen inversion in aniline with extensive geometry optimization

The geometry and energy of aniline have been calculated using the 6-31G and 6-31G^*(5D) basis sets for the planar structure and various pyramidal structures, assuming that the ring and the N-atom bonded to it lie in the same plane, but otherwise with full geometry optimization. With the 6-31G basis set the planar structure is predicated to be the most stable, whereas the inclusion of polarization functions in the 6-31G^*(5D) basis set finds a pyramidal structure with the out-of-plane angle =42.3° to be most stable and the energy barrier to inversion via the planar transition state to be 1.59±0.02 kcal mol^–1, in close agreement with experiment. Completing the optimization, allowing the N-atom and the C- and H-atoms of the ring to take up equilibrium out-of-plane positions increases the calculated energy carrier to inversion by less than 0.1 kcal mol^–1 to 1.66 kcal mol^–1. The ring adopts a very shallow inverted boat-type conformation with N₇-C₁C₄ = 2.0°.     

Quantum chemical study of the coordination of glycolic acid conformers and their conjugate bases to [Ca(OH2)n] (n=0–4) ions

A detailed exploration of the configurational and conformational space of glycolic acid and their conjugate bases has been carried out with the aid of first principles quantum chemical techniques at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory. The most stable configuration among the eight possible glycolic acid conformers corresponds to the E-s-cis, s-trans configuration, while the highest energy E-s-trans, s-cis conformer was found at 10.88 and 12.17 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. Upon dissociation of glycolic acid the s-cis(syn), and s-trans(anti) configurations of the glycolate anion can be formed. The anti conformer was found to be less stable than the syn one by 14.20 and 16.87 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p)) levels of theory, respectively. The computed B3LYP/6-311+G(d,p) proton affinity of the syn conformer for the protonation process affording the more stable E-s-cis, s-trans conformer, in vacuum was found to be 325.35 kcal mol⁻¹ (ΔG⁰ value). From a methodological point of view, our results confirm the reliability of the integrated computational tool formed by the B3LYP density functional model. This model has subsequently been used to investigate the interaction of Ca²⁺ ions with the glycolic acid conformers and their conjugate bases in vacuum and in the presence of extra water ligands. For the complexes of glycolic acid conformers the η²–O,O–(COOH) coordination, that is the structure that arises from the coordination of the Ca²⁺ to the carboxylic group, is the global minimum of the PES, while the η²–O(OH),O–(COOH) coordination is a local minimum found at only 1.0 and 1.3 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. Moreover, the two isomers exhibit nearly the same binding affinities, which are predicted to be 89 and 85 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. The same holds also true for the complexes of the glycolate anion. The η²–O,O–(COO⁻) coordination involving the syn conformer of the glycolato ligand, is the global minimum, while the η²–O(OH),O–(COO⁻) one lies at 1.5 and 5.6 kcal mol⁻¹ higher in energy at the B3LYP/6-311+G(d,p) and CCSD(T)/6-31G(d,p) levels of theory, respectively. The other conformer with an η²–O,O–(COO⁻) coordination involving the anti conformer of the glycolato ligand, is less stable by only 0.2 kcal mol⁻¹ at both levels of theory. Noteworthy is the trend seen for the incremental binding energy due to the successive addition of water molecules to [HOCH₂C(O)O]Ca²⁺ species; the computed values are 30.4, 26.8, 22.9 and 16.2 kcal mol⁻¹ at the B3LYP/6-311+G(d,p) level of theory for the mono-, di-, tri- and tetraaqua complexes, respectively. This trend arising from the repulsion of the dipoles between the water ligands and from unfavorable many body interactions is in accordance with those anticipated from electrostatic considerations. The Ca(II)-water interaction weakens with increasing coordination of the metal. Obviously, it is the electrostatic nature of the Ca(II)-water interactions that accounts well for the computed coordination geometries of the cationic (aqua)(glycolato)calcium complexes. Calculated structures, relative stability and bonding properties of the conformers and their complexes with [Ca(OH₂)_n]²⁺ (n=0–4) ions are discussed with respect to computed electronic and spectroscopic properties, such as charge density distribution, harmonic vibrational frequencies and NMR chemical shifts.     

Calculated structures and heats of formation of some predicted C,N, O,F molecules

Twelve molecules containing only C, N, O, and F, which are of interest in the context of energetic materials, were found to have true local energy minima at the HF/6-31G^* computational level. Both HF/6-31G^* and MP2/6-31G^* geometries are presented. A density functional procedure was used to find gas phase heats of formation. These are strongly positive for most of the molecules. The nitroso group significantly increases H _o ^f.25C compared to the nitro analogue.     

<Emphasis Type="Italic">Ab initio</Emphasis> Study of Structural and Conformational Properties of Five- to Nine-Membered Cyclic 1,2,3-Trienes

Summary.  The structures and relative energies of fundamental conformations of cyclopenta-1,2,3-triene, cyclohexa-1,2,3-triene, cylohepta-1,2,3-triene, cycloocta-1,2,3-triene, and cyclonona-1,2,3-triene were calculated by the HF/6-31G^* and MP2/6-31G^*//HF/6-31G^* methods. Only a C _2v symmetric planar conformation is available to cyclopenta-1,2,3-triene and cyclohexa-1,2,3-triene. The calculated energy barrier for ring inversion of the C _S symmetric puckerd conformation of cyclohepta-1,2,3-triene via the planar geometry is 62.2 kJ·mol⁻¹. The C ₂ symmetric twist conformation of cycloocta-1,2,3-triene was calculated to be the most stable one. Conformational racemization of the twist form takes place via the C _S symmetric half-chair geometry, which is by 60.8 kJ·mol⁻¹ less stable than the twist conformer. The C _S symmetric chair and unsymmetrical twist-boat conformations of cyclonona-1,2,3-triene were calculated to have similar energies; their interconversion takes place via an unsymmetrical low-energy (18.4 kJ·mol⁻¹) transition state. The twist (C ₂) and boat (C _S) geometries of cyclonona-1,2,3-triene are higher in energy by 13.2 and 33.9 kJ·mol⁻¹, respectively. Ring inversion in chair and twist-boat conformations takes place via a twist form as intermediate and requires 33.6 kJ·mol⁻¹. Corresponding author. E-mail: isayavar@yahoo.com Received March 25, 2002; accepted April 4, 2002     

Molecular Structures and Thermodynamic Stabilities of Cycloalkadienes and Their Methoxy Derivatives: An Ab Initio and DFT Study

The geometry-optimized molecular structures and total energies of 4- to 6-membered cycloalkadienes, and of a number of their monoand dimethoxy derivatives, have been calculated by ab initio (HF/6-31G^*, MP2/6-31G^*//HF/6-31G^*) and DFT (B3LYP/6-31G^*) methods. By comparison with available experimental data, the reliability of these computational methods for an estimation of the relative stabilities (enthalpies) of the isomeric forms of the title compounds was tested. The experimental enthalpies of isomerization proved to agree best with the respective theoretical data based on the mean of the HF/6-31G^* and B3LYP/6-31G^* energies. The theoretical calculations were then extended to several isomeric methoxy-substituted cycloalkadienes, for which no previous thermodynamic data exist. Some structural features of the title dienes were also discussed.     

Alkylhydrazides and their protonated forms: Physical methods andab initio calculations

The basicity of hydrazides of the highest aliphatic carboxylic acids RC(O)NHNH₂ (R = C_nH_2n+1,n = 5-12) has been studied by potentiometric titration, and IR and¹H NMR spectroscopy.Ab initio Hartree-Fork calculations using the 6–31G^* basis set with full optimization of geometry were carried out on the simplest acy1hydrazines and their possible protonated forms. Based on these calculations, and the 1R and¹H NMR spectra, the tautomerism of alkylhydrazides and the structures of their protonated forms are discussed.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 2645–2649, November, 1996.