Hydrogenation of the silanone groups (≡Si-O)2Si=0. Experimental and quantum-chemical studies

The reactivity of bis(siloxy)silanone groups (Si-0)₂Si=O stabilized on a silica surface with respect to H₂ molecules was studied. The reaction was found to give the (Si-O)₂SiH(OH) groups. The rate constant for this process was determined. Its activation energy in the 300–580 K temperature range is 13.4±0.3 kcal mol^–1, and the enthalpy is 54±5 kcal mol^–1. The activation energy for the reverse reaction,viz., elimination of a hydrogen molecule, is equal to 65 kcal mol^–1. Quantum-chemical calculations of hydrogenation of F₂Si=O and (HO)₂Si=O, which are the simplest molecular models of the silanone groups, were performed. Data on the geometrical and electronic structures of transition states and on the effects of substituents at the silicon atom on the reactivity of the silanone groups in this process were obtained. The optical absorption band of the surface silanone groups was quantitatively characterized. Its maximum is located at 5.65±0.1 eV; the extinction coefficient at the maximum (220 nm) is (3±0.5) · 10^–18 cm² molec.^–1.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 8, pp. 1951–1958, August, 1996.     

Proton migrations in associates of two molecules of formic acid with one molecule of hydrazine or hydrogen peroxide

The mechanisms of the proton transfer in associates of two molecules of formic acid with one molecule of hydrazine or hydrogen peroxide were studied usingab initio (SCFj6-31G^**) method. The mechanism of cooperative (concerted, one-step) four-proton transfer is realized in the associate with the hydrazine molecule. The proton transfer occurs stepwisevia an intermediate in the associate with a hydrogen peroxide molecule. The calculated activation barriers to the proton transfer in the associates investigated are 34.7 kcal mol^–1 and 27.1 kcal mol^–1, respectively.Translated fromlzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 11, pp. 2631–2635, November, 1996.     

Far-infrared spectrum,conformational stability,barriers to internal rotation,normal coordinate calculations,and vibrational assignment for chloroacetaldehyde

The far infrared spectrum [350 to 25 cm^–1] of gaseous chloroacetaldehyde, ClCH₂CHO, has been recorded at a resolution of 0.10 cm^–1. The first excited-state transition of the asymmetric torsion of the more stable near s-cis [chlorine atom s-cis to the aldehyde hydrogen atom] conformer has been observed at 26.9 cm^–1, with seven additional upper state transitions falling to higher frequency. Additionally, the fundamental torsional transition of the s-trans conformer has been observed at 58.9 cm^–1 with two excited states also falling to higher frequency. From these data, the asymmetric torsional potential coefficients have been determined to be:V ₁=414±11;V ₂ = 191±3;V ₃=–203±5;V ₄=44±1 andV ₆=–26±1 cm^–1. The s-cis to s-trans barrier is 500±5 cm^–1 (1.43±0.01 kcal mol^–1) with the s-cis conformer being more stable by 267±19 cm^–1 (0.76±0.05 kcal mol^–1) than the s-trans form. The Raman [4000 to 100 cm^–1] and infrared (4000 to 400 cm^–1] spectra of the gas have been recorded. Additionally, the Raman spectrum of the liquid has been recorded and qualitative depolarization values obtained. Complete vibrational assignments are proposed for both conformers based on band contours, depolarization values, and group frequencies. The assignments are supported by ab initio Hartree-Fock gradient calculations employing the 3–21G^* basis set to obtain the frequencies and the potential energy distributions for the normal vibrations for both rotamers. Additional ab initio calculations at the MP4/6-31G^* level have been carried out to determine the structural parameters for both conformers. The results are discussed and compared with the corresponding quantities obtained for some similar molecules.This contribution taken in part from the thesis of C. L. Tolley which will be submitted to the Department of Chemistry in partial fulfillment of the Ph.D. degree.     

The kinetics of the dissolution of monoclinic pyrrhotite in aqueous acidic solutions

Using the rotating disc method, the rates of dissolution of natural monoclinic pyrrhotite, FeS_1.14, in oxygen-free aqueous solutionsS([H⁺]=0.1, [Na⁺]=0.9, [ClO ₄ ^– ]=1.0 mol kg^–1) were determined. In the temperature range 40–90 °C the dissolution reaction occurs under kinetic control; the activation energy being 14±1 kcal mol^–1 (50±5 kJ mol^–1).

---

Die Kinetik der Auflösung von monoklinem Pyrrhotin in sauren wäßrigen Lösungen

Zusammenfassung Die Auflösungsgeschwindigkeit von natürlichem monoklinen Pyrrhotin, FeS_1.14, wurde in sauerstofffreien LösungenS([H⁺]=0.1, [Na⁺]=0.9, [ClO ₄ ^– ]=1.0 mol kg^–1) mit Hilfe der Methode der rotierenden Scheibe bestimmt. Im Temperaturbereich von 40–90° erfolgt die Auflösungsreaktion kinetisch kontrolliert, wobei eine Aktivierungsenergie von 14±1 kcal mol^–1 (59±5 kJ mol^–1) gefunden wurde.

     

Very empirical treatment of solvation and entropy: a force field derived from Log Po/w

A non-covalent interaction force field model derived from the partition coefficient of 1-octanol/water solubility is described. This model, HINT for Hydropathic INTeractions, is shown to include, in very empirical and approximate terms, all components of biomolecular associations, including hydrogen bonding, Coulombic interactions, hydrophobic interactions, entropy and solvation/desolvation. Particular emphasis is placed on: (1) demonstrating the relationship between the total empirical HINT score and free energy of association, G _interaction; (2) showing that the HINT hydrophobic-polar interaction sub-score represents the energy cost of desolvation upon binding for interacting biomolecules; and (3) a new methodology for treating constrained water molecules as discrete independent small ligands. An example calculation is reported for dihydrofolate reductase (DHFR) bound with methotrexate (MTX). In that case the observed very tight binding, G _interaction–13.6 kcal mol^–1, is largely due to ten hydrogen bonds between the ligand and enzyme with estimated strength ranging between –0.4 and –2.3 kcal mol^–1. Four water molecules bridging between DHFR and MTX contribute an additional –1.7 kcal mol^–1 stability to the complex. The HINT estimate of the cost of desolvation is +13.9 kcal mol^–1.     

Association of metal cations with alkanes: Na(CH4)+ versus Cu(CH4)+ as molecular models

Summary Ab initio molecular orbital calculations give small stabilization energies for the various Na(CH₄)⁺ adducts (less than 4 kcal mol^–1), but predict a stronger binding for the copper compounds (about 13 kcal mol^–1). The different behaviour of Na⁺ and Cu⁺, already present at the SCF level, is reinforced by electron correlation. This can be attributed to an important contribution of the dispersion energy to the binding energy of the copper ion: about 40% of the total, including basis set superposition corrections.Dedicated to Mrs A. Pullman     

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.     

Double proton migrations in the associates of formic acid with nitrous,nitric, orthophosphoric,and sulfuric acids

Mechanisms of the proton transfer in dimeric associates of formic acid with nitrous, nitric, orthophosphoric, and sulfuric acids were studied by theab initio (HF/6-31G^**) method. The mechanism of the cooperative (concerted or one-step) proton transfer was shown to occur in all cases. The calculated activation barriers of the proton transfer reactions for the associates investigated are equal to 19.9, 14.2, 13.3, and 10.7 kcal mol^–1, respectively.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 2184–2189, September, 1996.     

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.     

Ab initio MP2 calculations of the products of acetylene addition to HgCl2

Gas-phase reaction of acetylene with HgCl₂ resulting in -chlorovinylmercury derivatives and their interaction with Cl^– and I^– anions and KI molecule was studied by the ab initio MP2 method with the Dunning—Hay double zeta basis set and LanL pseudopotential for Hg, K, and I atoms. The reaction was shown to proceed via a -complex of acetylene and HgCl₂ (the calculated enthalpy of formation is –6.5 kcal mol^–1). According to calculations, the activation energy of formation of cis--chlorovinylmercury chloride from acetylene and HgCl₂ is 31 kcal mol^–1. Chloride and iodide anions and KI molecule are readily added to both cis- and trans-isomer of -chlorovinylmercury chloride to give stable species.     

Thermal isomerization of acetylnitrene: a quantum-chemical study

The electronic structure and pathways of thermal isomerization of formylnitrene and acetylnitrene were studied by the B3LYP/6-311G(d,p) density functional method and ab initio G2(MP2,SVP) computational procedure using the geometries obtained from B3LYP calculations. According to G2 calculations, both nitrenes have singlet ground states while the energies of the corresponding triplet states are 2.8 and 5.7 kcal mol^–1 higher. For acetylnitrene, the activation barrier to the nitrene isocyanate isomerization was estimated at 28.9 kcal mol^–1 (G2). Calculations revealed no pathway for single-step isomerization of nitrene into cyanate in both systems. The formation of methyl cyanate from isocyanate is thermodynamically unfavorable (E = 26.5 kcal mol^–1) and requires a high activation barrier (89.4 kcal mol^–1) should be overcome. Based on the results obtained, the pathways of transformation of nitrene formed in thermal decomposition of acetyl azide (Curtius rearrangement) were analyzed.     

Multidimensional tunneling dynamics of proton transfer in malonaldehyde molecule and its isotopomers

The twenty-one-dimensional Hamiltonian of malonaldehyde molecule and a number of its isotopomers (H/D, ¹³C/¹²C) was reconstructed in the low-energy region (<3000 cm^–1). Parameters of the Hamiltonian were obtained from quantum-chemical calculations of the energies, equilibrium geometries, and eigenvectors and eigenfrequencies of normal vibrations at the stationary points corresponding to the ground state and transition state. Despite substantial variation of the barrier height calculated using different quantum-chemical methods (from 2.8 to 10.3 kcal mol^–1), the corresponding potential energy surfaces can be matched with high accuracy by scaling only one parameter (the semiclassical parameter , which defines the scales of potential, energy, and action). Scaling invariance allows optimization of the Hamiltonian in such a way that the calculated ground-state tunneling splitting coincides with the experimental value. The corresponding potential barrier height is estimated at 4.34±0.4 kcal mol^–1. The quantum dynamics problem was solved using the perturbative instanton approach without reducing the number of degrees of freedom. The role of all transverse vibrations in proton tunneling is characterized. Vibration-tunneling spectrum is calculated for the ground state and low-lying excited states and mode-specific isotope effects are predicted.     

Quantum-chemical calculations of the dissociation energy of the C-H bond in hydrocarbons,alcohols, and ethers

The dissociation energy of the C-H bonds in hydrocarbons, alcohols, and ethers were calculated by semiempirical MNDO, AM1, and PM3 methods. The average error of calculations of theD(C-H) values by using various quantum-chemical methods is 1.3 kcal mol^–1.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2858–2861, December, 1996.     

Far infrared spectrum,conformational stability,barriers to internal rotation,vibrational assignment,and ab initio calculations of 2-chloro-3-fluoropropene

The far infrared spectrum (375 to 30 cm^–1) of gaseous 2-chloro-3-fluoropropene, CH₂=C(CH₂F)CI, has been recorded at a resolution of 0.10 cm^–1. The fundamental asymmetric torsional mode is observed at 117.5 cm^–1 with ten excited states falling to low frequency for thes-cis (fluorine atom eclipsing the double bond) conformer. For the higher energy gauche conformer, the asymmetric torsion is estimated to be at 94 cm^–1. From these data the asymmetric torsional potential function has been calculated. The potential function coefficients are calculated to be in cm^–1):V ₁=803±21,V ₂=–94±21,V ₃= 1025±10,V ₄=95±10, andV ₆=2±1, with an enthalpy difference between the more stables-cis and gauche conformera of 550±100 cm^–1 (1.57±0.29 kcal/mol). This function gives values of 1227±50cm^–1(3.51±0.14kcal/mol), 1266±200 cm^–1 (3.62±0.57 kcal/mol), and 665±100 cm^–1 (1.90±0.29 kcal/mol), for thes-cis to gauche, gauche to gauche, and gauche tos-cis barriers, respectively. From the relative intensities of the Raman lines of the gas at 652 cm^–1 (gauche) and 731 cm^–1 (s-cis) as a function temperature, the enthalpy difference is found to be 565±96 cm^–1 (1.62±0.27 kcal/mol). However, the more polar gauche conformer remains in the crystalline solid. The Raman spectrum of the gas has been recorded from 3500 to 70 cm^–1 and, utilizing these data and the previously reported infrared data, a complete vibrational analysis is proposed for both conformers. The conformational stability, barriers to internal rotation, fundamental vibrational frequencies, and structural parameters that have been determined experimentally are compared to those obtained from ab initio Hartree-Fock gradient calculations employing both the 3–21 G^* and 6–31G^* basis sets and to the corresponding quantities for some similar molecules.     

A theoretical study of the relative stability of the isomeric forms of N2O3

The electronic structure, geometrical parameters and relative stability of the isomeric forms of N₂O₃ are analysed by means of ab initio calculations. Total energies of the different isomers are given. The energy difference between the most stable conformers of the symmetric N₂O₃ is 4.31 Kcal mol^–1 as provided by 6–31G basis set. The height of the rotational barrier determined by the ab initio technique is 7.12 kcal mol^–1.Member of the Carrera del Investigador CICPBA, R. Argentina.Member of the Carrera del Investigador CONICET, R. Argentina.Predoctoral fellow of CONICET, R. Argentina.     

Stereochemical nonrigidity of (O—Ge)-chelate bis[(2-oxo-1-hexahydroazepinyl)methyl]dichlorogermane: substantiation of the dissociative mechanism of chelating ligand exchange

Stereochemical nonrigidity of the hexacoordinated (O—Ge)-chelate bis(2-oxo-1-hexahydroazepinylmethyl)dichlorogermane in CDCl₃ was studied by dynamic NMR. The activation parameters of the intramolecular rearrangement at the coordination center are G ^# ₂₉₈ = 12.3±0.2 kcal mol^–1, H ^# = 16.9±0.2 kcal mol^–1, and S ^# = 15.3±0.7 cal mol^–1 K^–1. The dissociative mechanism of ligand exchange involving the cleavage of the OGe coordination bond is discussed based on the positive entropy of activation.     

Modeling the Proton Transport in Orthoperiodic and Orthotelluric Acids and Their Salts

Orthoperiodic and orthotelluric acids, their salts MIO₆H₄ (M = Li, Rb, Cs) and CsH₅TeO₆, and dimers of the salt · acid type are calculated within density functional theory B3LYP and basis set LanL2DZ complemented by the polarizationd,p-functions. According to calculations, the salt · acid dimerization is energetically favorable for compounds MIO₆H₄ · H₅IO₆ (M = Rb, Cs) and CsIO₆H₄ · H₆TeO₆. The dimerization energy is equal to 138–146 kJ mol^–1. With relatively small activation energies equal to 4 kJ mol^–1 (M = Li) and 11 kJ mol^–1 (M = Rb, Cs), possible is rotation of octahedron IO₆ relative to the M atom in monomers of salt molecules. The proton transfer along an octahedron occurs with activation energies of 63–84 kJ mol^–1. The activation energy for the proton transfer between neighboring octahedrons of the type salt · acid acid · salt equals 8–17 kJ mol^–1. Quantum-chemical calculations nicely conform to x-ray diffraction and electrochemical data.     

Conformational mobility of the six-membered ring in 5-oxo-1,3-cyclohexadiene and its derivatives

Molecular mechanics and MNDO calculations showed that the six-membered ring in the molecule of 5-oxo-1,3-cyclohexadiene possesses high conformational mobility. The transition from a planar equilibrium conformation to a distorted sofa conformation in which the C(sp²)-C(=O)-C(sp³)-C(sp²) torsion angle is equal to ±30° increases the energy of the molecule by less than 1 kcal mol^–1. The influence of steric (R = Me, Et, Prⁱ, Bu^t) and electronic (R = NH₂, NO₂) effects of substituents R on the equilibrium conformation and mobility of the carbocycle has been analyzed. Both types of substituents at unsaturated C atoms do not change the equlibrium conformation or flexibility of the six-membered ring. Substituents at saturated C atoms cause the transition of the carbocycle to the distorted sofa conformation and significantly restrict its mobility. The electronic structures of 5-oxo-1,3-cyclohexadiene and its amino and nitro derivatives have been analyzed.Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 849–854, May, 1995.     

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.     

Conformational stability,barriers to internal rotation,normal coordinate analysis,and vibrational assignment of chloroacetyl bromide

The infrared spectra (3200 to 30 cm^–1) of gaseous and solid chloroacetyl bromide, CH₂ClC(O)Br, and the Raman spectra (3200 to 10 cm^–1) of the gas, liquid (with depolarization data), and solid have been recorded. From the observed asymmetric torsional transitions, the potential function governing internal rotation of the CH₂Cl moiety has been determined with the following coefficients:V ₁=336±11,V ₂=73±10,V ₃=757+7,V ₄=103±3, andV ₆=5±2 cm^–1. This potential function is consistent with s-trans to gauche and gauche to gauche barriers of 963±11 and 709±12cm^–1, respectively, and enthalpy difference of 373 ± 24 cm^–1 with the dihedral angle of the gauche rotamer being 115°. The enthalpy difference has been determined experimentally from the studies of the Raman spectra at different temperatures to be 359±68 cm^–1 (1.03±0.19 kcal mol^–1) and 507±24 cm^–1 (1.45±0.07 kcal mol^–1) for the gas and liquid, respectively, with the s-trans conformer being the more stable conformer in the gas and liquid and the only one present in the annealed solid. A complete assignment of the vibrational fundamentals is proposed from spectral data obtained for the gas, liquid, and solid. The assignment is supported by a normal coordinate calculation utilizing a modified valence force field to obtain the frequencies for the normal vibrations and the potential energy distribution. The results are discussed and compared to the corresponding quantities for some similar molecules.Taken in part from the thesis of H. V. Phan, which will be submitted to the Department of Chemistry in partial fulfillment of the Ph.D. degree.