An ab initio study of the geometry,harmonic and anharmonic force fields,and fundamental vibrational frequencies of cis- and trans- thiolformic acid

The geometry, harmonic and anharmonic force fields, and fundamental vibrational frequencies of cis- and trans-thiolformic acid are studied ab initio in the 4-31G basis set. An extensive comparison is made between changes in diagonal and off-diagonal quadratic and cubic force constants and diagonal stretching quartic constants in going from the chain to the ring structure in thiolformic acid and formic acid. The changes in the force constants are indicative of a much weaker interaction in the trans conformer between S? H and O?C, compared with O–H and O?C, in keeping with the weaker hydrogenbonding property of the S? H group in general.     

An ab initio molecular orbital study of the grignard reagents CH3MgCl and [CH3MgCl]2: The Schlenk equilibrium

Ab initio molecular orbital calculations are used to study the modified Schlenk equilibrium: 2RMgCl (RMgCl)² MgR₂ + MgCl₂ Mg(Cl₂)MgR₂ with R=H and CH₃. In the absence of any solvents, calculations indicate that the formation of the various possible bridged dimers (RMgCl)₂ is substantially exothermic. However, using dimethylether as a model solvent, we show that the formation of the dimer (Me₂O)(CH₃)Mg(Cl₂)Mg(CH₃)(OMe₂) is exothermic only when entropic effects are included.     

Binding of the two-valence-electron metal ions Sc(+), Ti(2+), and V(3+) to the second-row hydrides BeH(2), BH(3), NH(3), OH(2), and FH: a computational study

We report results from a computational study of the binding in complexes formed from one of the transition-metal ions Sc(+), Ti(2+), or V(3+), each of which has two valence electrons outside an argon core, and one of the second-row hydrides FH, OH(2), NH(3), BH(3), or BeH(2). The complexes that involve the electron-rich ligands FH, OH(2), and NH(3) have strong ion-dipole components to their binding. There are large stabilization energies for sigma-interactions that transfer charge from occupied lone-pair natural bond orbitals on the F, O, or N atom of the (idealized) Lewis structure into empty non-Lewis orbitals on the metal ions; these interactions effectively increase electron density in the bonding region between the metal ion and liganded atom, and the metal ions in these complexes act in the capacity of Lewis acids. The complexes formed from the electron-poor hydrides BH(3) and BeH(2) consistently incorporate bridging hydrogen atoms to support binding, and there are large stabilization energies for interactions that transfer charge from the Be-H or B-H bonds into the region between the metal ion and liganded atom. The metal ions in Sc(+)-BeH(2), Ti(2+)-BeH(2), Ti(2+)-BH(3), and V(3+)-BH(3) act in the capacity of Lewis acids, whereas the scandium ion in Sc(+)-BH(3) acts as a Lewis base.     

Intramolecular Nonbonded Interactions Between Oxygen and Group VIA Elements: An Ab Initio Molecular Orbital and Density Functional Theory Investigation of the Structures of HX—CH2—COOH (X=S,Se, and Te)

Ab initio molecular orbital and density functional methods have been used to study the potential energy surfaces of the substituted acetic acids HX—CH₂—COOH, where X is one of the Group VIA Chalcophiles S, Se, or Te. The various conformers adopted by these compounds provide information regarding the energetic importance of nonbonded and hydrogen bonding interactions involving oxygen atoms with different hybridizations. Density functional and ab initio molecular orbital methods yield similar structural and energetic trends for these compounds. Calculations show that the structure of the lowest-energy conformer of each of these acids has the X—C—C—O backbone substantially twisted from planarity, similar to that previously observed for the corresponding aldehydes, HX—CH₂—CHO. In the twisted acid structures the shortest distance is within about 0.1 Å of the sum of the X and O van der Waals radii, which reduces overcrowding of the lone pairs of electrons on these atoms. In conformers where the heavy atom backbone is planar, one of the distances is significantly shorter than the sum of the van der Waals radii, and the total molecular energy of these conformers is higher than that of the twisted forms. The variation of X—H vibrational frequencies among conformers reflects the extent of X—H hydrogen bonding, and indicates that formation of this hydrogen bond is not the dominant factor in determining the lowest-energy conformation. When X is oxygen (HO—CH₂—COOH), the lowest-energy conformer is also nonplanar, whereas for the corresponding aldehyde, HO—CH₂—CHO, the lowest-energy conformer is a planar structure with C_S symmetry. The conformational preferences of these simple species provide reference points for inter- and intramolecular interactions in more complex systems of biological interest.     

An ab initio study of some structural features and the rotational barriers in performic acid,formic acid and hydrogen peroxide

Calculations on performic acid at the 4-31G level, with and without bond functions with complete geometry optimization, and at the (9, 5) level, with and without polarization functions and rigid rotation, all give no sign of a well in the potential energy curve for rotation about the O/O bond axis in the region of 50° – 90° ; and all but the unaugmented 4-31G basis set find the cis-cis planar conformer to be the most stable form. Calculations at the (9,5) level with rigid rotation find the energies of the other planar conformers, relative to the cis-cis conformer, to be 0.94, 1.50 and 14.80 kcal mol^?1 for the trans-trans, cis-trans and trans-cis structures respectively. These energies and also that for the barrier separating the cis-cis and cis-trans conformers, 1–2 kcal mol^?1, are discussed in relation to corresponding data for formic acid, hydrogen peroxide and several planar four heavy-atom molecules. Dipole moment calculations using the same basis sets would seem to favor a skew conformation as the most stable for performic acid, but comparisons between calculated and experimental values for formic acid and for hydrogen peroxide cast doubt on the validity of such results.     

An ab initio study of the geometry,energy, and selected force constants for the three planar conformers of carbonic acid,and the bicarbonate ion; and of the energy for the reaction H2O + CO2 → H2CO3

Ab initio calculations using the unscaled 4-31G basis set have been carried out on the cc, tc, and tt conformers of carbonic acid and the bicarbonate ion, with full geometry optimization assuming the structures to be planar. The complete harmonic force field is reported for the (most stable) tt conformer and for the bicarbonate ion, also selected quadratic force constants for the cc and tc conformers. The changes in certain bond lengths and stretching force constants in the cc → tc, tc → tt, and cc → tt conformer conversion reactions are indicative of intramolecular hydrogen bonding, C?O…H? O and H? O…H? O, which is examined in greater detail by partitioning the overall conformer conversion energy into distortion and bonding energy components. The fundamental vibration frequencies for the tt conformer and the bicarbonate ion are calculated from the force constant matrices, and hence, using a scaling factor based on a comparison of calculated and experimental values for the bicarbonate ion and trans-formic acid, a value is predicted for the zero-point energy of the tt conformer. A new estimate of ΔH? for the hydration reaction, H₂O + CO₂ → H₂CO₃, at 298 K in the gas phase; is made from thermochemical data, +20.2 ± 3.4 kJ mol^?1, which, together with estimates of (H₂₉₈? – H₀?) and the zero-point energy for H₂CO₃, gives +8.1 ± 7.0 kJ mol^?1 for ΔE_T(expt). ΔE_T calculated from the 4-31G basis set data is -29.1 kJ mol^?1. Comparison of the experimental value, the Hartree–Fock limit value, and values calculated with a variety of basis sets for the bond separation reaction, CO₂ + CH₄ → 2H₂CO, suggests that the differences, ΔE_T(expt) minus ΔE_T(SCF ), are due mainly to basis set limitations and not substantial correlation energy contributions.     

A molecular orbital study of the rotation about the C-C bond in 1,3-butadiene

The geometry and energy of 1,3-butadiene have been calculated using the 6-311G** basis set as a function of the CCCC dihedral angle-0 ° (trans), 30 °, 60 °, 75 °, 90 °, 120 °, 135 °, 150 °, 165 ° and 180 ° (cis)-assuming that the vinyl groups remain planar. Potential minima are located at 0 ° and 141.4 °, with the trans structure more stable than the gauche by 13.2 kJ mol^–1. Potential maxima are located at 76.7 °, giving a barrier height of 25.4 kJ mol^–1 relative to the trans structure, and at 180 ° giving a barrier height of 3.0 kJ mol^–1 relative to the 141.4 °-gauche structure. Using the 6-31G* basis set the inclusion of electron correlation, accounting for about 52% of the correlation energy, was found to produce no significant change in the shape of the potential energy curve. The magnitude of the expectation energy differences is such that both barriers with respect to the 14l.4 °-gauche maximum structure can be categorized unequivocally as attractive-dominant, whereas the values for the energy barrier with respect to the trans structure, although characteristic of a repulsive-dominant barrier at the 6–311G** level, are sufficiently small that higher level calculations might give the opposite result. Analysis of V _nn for the conversion reactions cis 150 °-gauche, trans 60 °-gauche, and trans 90 °-gauche in terms of the individual contributions from the various internuclear interactions shows that nonbonded interactions are important, not only in initiating the destabilization of the crowded cis structure, but also through-out the entire range of CCCC dihedral angles, 0 ° to 180 °.     

Rotation about the C?N bond in 2-aza-1,3-butadiene and the N?N bond in 2,3-diaza-1,3 butadiene: A molecular orbital study

The geometry and energy of 2-aza-1,3-butadiene and 2,3-diaza-1,3-butadiene have been calculated using the 6-31G* basis set as a function of the CNCC and CNNC dihedral angles, respectively. With the 2-aza derivative potential minima are located at 0° (trans) and at about 130° for a gauche structure approximately 9.5 kJ mol^?1 less stable than the trans. Potential maxima are at about 75° giving a gauche barrier height of approximately 19 kJ mol^?1 relative to the trans structure, and at 180° (cis) giving a barrier height of approximately 14.5 kJ mol^?1 relative to the 130° gauche structure. With the 2,3-diaza derivative the gauche barrier has disappeared and there are a series of gauche structures in the region 70°–100° of almost equal energy 12.5-15 kJ mol^?1 less stable than the trans. In addition the cis barrier is much greater, nearly 70 kJ mol^?1 relative to the trans structure. Inclusion of electron correlation, accounting for about 50% of the correlation energy, produces no significant changes in the shape of the potential energy curves. There are systematic and progressive changes in almost all the geometrical parameters as the ?CH? groups in butadiene are replaced by ?N? . The outward tilt and compression within the methylene groups show adverse steric interactions to be operative in the cis structures. The values of V_nn indicate that gauche structures of both the 2-aza and the 2,3-diaza derivatives near the cis structure are more compact (as with butadiene), and gauche structures of the 2-aza derivative near the trans structure are less compact (as with butadiene). Originating in the changes in bond lengths and bond angles, rotation-independent nuclear–nuclear interactions again play an important role.     

Intramolecular interaction effects and the structure of N2OS and N2O2: An Ab initio study

An ab initio study of O?N? N?S with full geometry optimization has been carried out to corroborate the presence of an interaction between the terminal atoms in this type of structure, which, in O?N? N?O, apparently stabilizes the cis conformer. Using the unscaled 4–31G basis set with a full set of d functions on the sulfur, there is a potential minimun at the trans but not the cis geometry. A gauche conformer with a torsional angle of 77.2° is the most stable. With N₂O₂ this basis set gives potential minima at both the cis and trans geometries, but the trans conformer is slightly more stable, contrary to experiment and the results of (7,3) basis-set calculations reported in the literature in which Gaussian lobe functions were employed. Using a (9,5) basis set there is no longer a potential minimum at the cis geometry, and a gauche structure is more stable than the cis conformer as in the case of N₂OS with the less-extended basis set. Force constants (harmonic and anharmonic), compliance constants, relaxed force constants, and interaction-displacement coordinates for both molecules are compared for key structural elements.     

The use of 90°-1,3-butadiene as a reference structure for the evaluation of stabilization energies for benzene and other conjugated cyclic hydrocarbons

Reactions are described that employ 90°-1,3-butadiene as a reference structure for the evaluation of the stabilization energyof the benzenoid and other cyclic conjugated hydrocarbons. The unique benefits of this rotamer of butadiene as a reference molecule within the homodesmotic conceptual framework are discussed. Experimental stabilization energies are presented for a number of cyclic hydrocarbons.