Pathways of the reaction of nucleophilic addition of ammonia to formaldehyde in the gas phase and in the complex with formic acid:ab initio calculations

The gradient pathways of the reaction of nucleophilic addition of ammonia to formaldehyde were calculated for free molecules and in the NH₃...H₂CO...HC(O)OH complex by theab initio RHF/6-31G^**, MP2(fc)/6-31G^**, and MP2(full)/6-311++G^** methods. Both reactions proceed concertedly. In the first case, the reaction successively passes through two transitional states with an energy barrier exceeding 35 kcal mol⁻¹. In the case of the complex with formic acid, the reaction follows a conventional pathway, although its activation barrier calculated by the RHF/6-31G^** and MP2(fc)/6-31G^** methods decreases to 12.6 and 3.8 kcal mol⁻¹, respectively. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 13–20, January, 1998.     

A study of methanetetraol dehydration to carbonic acid

Two alternative dehydration reactions C(OH)₄ → (HO)₂CO + H₂O and C(OH)₄ + H₂O → (HO)₂CO + 2H₂O are studied by ab initio Becke3LYP/6–311 + G^** and MP2/6–31G^** methods. Calculated energy and geometry characteristics of intermediates and transition states predict a catalytic effect of one water molecule and the exothermism of the transformations. Relevant HF/6–311 + G^**, HF/6–31G^**, HF/6–31G, and HF/3–21G calculations were performed for comparison. © 1997 John Wiley & Sons, Inc.     

Stabilization of nonclassical types of valence bond orientation at the carbon atom in organoboron compounds

The complete topological structure of the potential energy surface (PES) of methane in the inversion region was studied by theab initio CCD(full)/6-311++G^** method. The necessity of taking into account nuclear motions was shown. Penta- and hexacoordination of carbon atoms in boron-containing organic compounds was investigated by theab initio MP2(full)/6-31G^** and MP2(full)/6-311++G^** methods. The CB₄H₄Li₂, CB₆N₂H₂, and CB₆O₂ systems containing hexacoordinated carbon atoms correspond to rather deep minima on the relevant PES and can be the subject of synthetic studies. According to theab initio calculations, pyramidal boron-containing systems with hypercoordinated carbon atoms, which fulfill the “8e rule,” also correspond to rather deep PES minima and can be detected experimentally. Dedicated to Academician V. I. Minkin on the occasion of his 65th birthday. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 5, pp. 786–796, May, 2000     

Intramolecular nucleophilicS N 2 substitution at the tetrahedral carbon atom: Anab initio study

Pentacoordination of carbon atom in bicyclic organic compounds of the pentalene type was studied by theab initio RHF/6-31G^** and MP2(full)/6-31G^** methods. It was shown that intramolecularS _N 2 reactions with energy barriers within the energy scale of NMR spectroscopy can occur in systems in which a linear orientation of the attacking and leaving groups is realized. The barrier to the intramolecular nucleophilic substitution reaction in 2,3-dihydro-3-formylmethylenefuran is 36.9 (RHF) and 27.7 kcal mol⁻¹ (MP2) and decreases to 16.4 and 19.4 kcal mol⁻¹, respectively, in the case of diprotonation at the O atoms in this system. For model pentalene type compounds containing electron-deficient B atoms in the ring, theab initio calculations predict a further decrease in the barrier height (down to less than 10 kcal mol⁻¹). Translated fromIzvestiya Akademii Nauk, Seriya Khimicheskaya, No. 7, pp. 1246–1256, July, 1999.     

Intramolecular Hypervalent X←X Interaction in Mono- and Bicyclic Heteroelement Compounds Containing Se and Te Atoms

Chalcogen-containing heterapentalene and quasimonocyclic compounds having SeÄSeÄSe and TeÄTeÄTe triads or SeÄSe and TeÄTe diads were studied by the ab initio [MP2(full)/6-31G**, MP2(fc)/6-31+G**, and MP2(fc)/LANL2DZ] and DFT methods (B3LYP/6-31G^**, B3LYP/6-31+G^**, and B3LYP/LANL2DZ). Heterapentalene compounds were found to be stable as planar bicyclic structures having a C _2v symmetry. The stability of quasimonocycic -chalcogenovinyl aldehydes increases with increase in the electron-acceptor power of the substituent at the X atom.     

Cooperative atomic migrations in the associates of formic acid molecule with H2O...X systems (X=CH3OH,NH2OH,H2O2, HOF,and H2O)

The pathways and activation barriers of cooperative biproton migrations in the associates of the formic acid molecule with H₂O and X molecules (X=CH₃OH, NH₂OH, H₂O₂, FOH, and H₂O) are calculated by an ab initio method (3-21G and 6-31G^** basis sets). A cooperative triproton transfer occurs in the system with X=H₂O. The activation barriers of this transfer calculated in the 3-21G and 6-31G^** basis sets are 6.94 and 27.29 (through the structure of C₂ symmetry) or 7.99 and 26.08 kcal/mole (through the structure of C_s symmetry), respectively. In the systems with X=H₃COH, HOOH, and FOH, the biproton transfer is accompanied by synchronous shifts of two hydroxyl groups and overcomes high activation barriers (>40 kcal/mole), which is accounted for by poor stereochemical similarity for the low-barrier cooperative processes in the given molecular associates. Scientific Research Institute of Physical and Organic Chemistry, Rostov State University. Translated fromZhurnal Strukturnoi Khimii, Vol. 37, No. 5, pp. 845–858, September–October, 1996. Translated by I. Izvekova     

Hydration of the carboxylate group: Anab initio molecular orbital study of acetate-water complexes

The hydration of the carboxylate group in the acetate anion has been investigated by performingab initio molecular orbital calculations on selected conformers of complexes with the form CH₃CO₂ ^– ·nH₂O·mH₂O, wheren andm denote the number of water molecules in the first and second hydration spheres around the carboxylate group, andn + m 7. The results of RHF/6–31G^* optimizations for all the complexes and MP2/6–31+G^** optimizations for several one-water complexes are reported. The primary consequence of hydration on the structure of the acetate anion is a decrease in the length of the C-C bond. Enthalpy and free energy changes calculated at the MP2/6–31+G^** and MP2/6–311+ +G^** levels are reported for the reactions CH₃CO₂ ^– + [H₂O] _P CH₃CO₂ ^– ·nH₂O ·mH₂O where [H₂O] _P is a water cluster containingp water molecules andp=n+m 7. The calculations show that conformers with the lowest enthalpy change on complex formation are often not those with the lowest free energy change, due to a greater entropic loss in complexes with tighter and more favorable enthalpic interactions. Hydrogen bonding of six water molecules directly to the carboxylate group in CH₃CO₂ ^– is found to account for approximately 40% of the enthalpy change and 37% of the free energy change associated with bulk solvation.     

Mechanisms of inversion of bond configuration at the tetrahedral boron atom in five-membered chelate cycles

Mechanisms of inversion of the bond configuration at the tetrahedral boron center in five-membered chelate cycles of the 1,3,2-oxazaborolidine and 1,3,2-oxazaborolidene molecules were studied by theab initio MP2(full)/6-31G^** method. It was shown that enantiotopomerization occurs by a dissociative mechanism with the cleavage of the B←N bond and the formation of acyclic intermediates with tricoordinate planar boron atom. The calculated energy barriers to inversion of tetrahedral bond configurations at boron centers in the two chelate complexes are equal to 13.1 and 15.4 kcal mol⁻¹, respectively. In contrast to 1,3,2-oxazaborolidine, internal rotation about the B−O bond in its unsaturated analog makes an appreciable contribution to the reaction coordinate. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 2, pp. 250–255 February, 1999.     

Anion-molecular interaction through CH3 groups: Model ab initio studies

3-21G, 6–31G and 6–31 +G calculations have been performed on Cr…H₃CCl and LiCl…H₃CCl complexes with two different configurations each. Optimized geometries, stabilization energies, CH force constants and harmonic vibrational frequencies for CH₃Cl and its complexes are reported. Comparison of the calculated frequency shifts of CH-stretching bands of CH₃C1, on complexation, with experimental results of related systems indicate that the interaction of Cl^- with CH₃ group takes place in a linear manner with the CH bond. A lower frequency shift for LiCl…H₃CX with reference to the C1^-…H₃CX complexes is explained on the basis of the reduction of the basicity of Cl^- ion in the presence of counter ion in the former complex considered.     

Structural peculiarities and the possibility of the existence of small water cluster anions (H2O) n − withn≤4

Structures of (H₂O) _n ⁻ anions withn≤4 were optimized at the UHF/4-31++G^** level and their stability was estimated at the MP2/4-31++G^** level. The trimer anion has a chain-like structure while the tetramer anion can exist either in a chain-like or a cyclic configuration. In the dimer anion and in the chain-like anions, the excess electron density is localized on the terminal water molecule, an acceptor of the H-bond proton. In the cyclic anion, it is uniformly distributed over the free hydrogen atoms. All considered anions have energy values higher than those of the corresponding neutral oligomers. The detachment of an electron from the anions should proceed with the liberation of energy. However, trimer and larger anions are stable against dissociation into individual water molecules and a free electron. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 1, pp. 41–46, January, 1997.     

The formation of formic acid upon low-temperature condensation of CO2−H2 and CO−H2O gas mixtures dissociated in electric discharge

The formation of formic acid in the low-temperature condensation of CO₂−H₂ and CO−H₂O gas mixtures dissociated in electric discharge was investigated. The gas-phase concentrations of H^., O^., OH^., and O₂ were measured downstream a microwave discharge in a CO₂−H₂ mixture. Low-temperature (77 K) condensates formed from CO₂−H₂ and CO−H₂O mixtures were studied by ESR. The formation of formic acid in the CO₂−H₂ and CO−H₂O systems was found to be due to the reactions of H^., CO, O^., and O₂ on the condensate surface. A single mechanism of the formation of formic acid in the CO₂−H₂ and CO−H₂O systems was proposed.     

Theoretical Study of the Interactions between Isolated <Emphasis Type="Italic">DNA</Emphasis> Bases and Various Groups IA and IIA Metal Ions by <Emphasis Type="Italic">Ab Initio</Emphasis> Calculations

Summary. Interactions of the DNA bases adenine (A), guanine (G), cytosine (C), and thymine (T) with various metal ions (M) of groups IA and IIA of the periodic table of the elements were studied at the HF, MP2, and DFT levels of theory. The structures and thermodynamic stabilities of these species were studied at the gas phase. The calculations uphold that there exist two active sites in G and one in A, C, and T. The calculations also show that the O² atom in T is a more active site for metal ion bindings than that in C. The stability energies for G … M complexes are larger than those for A … M complexes and the stability energies for T … M complexes are larger than those for C … M complexes. As z/r ratio for the metal ion increases, the interaction energy for the complex increases systematically. Thermodynamic quantities such as ΔH, ΔG, ΔS, and ln K were determined for each complexation reaction, [Base+M ⁿ⁺ →(Base … M) ⁿ⁺]. A, G, and C complexation reactions except for C … Rb⁺ are exothermic. The situation is quite different for T complexation reactions and all except for T … Be²⁺ and T … Mg²⁺ are endothermic.     

Pentacoordination of boron, carbon, aluminum, and silicon atoms in organic compounds: Anab initio study

Pentacoordination of boron, carbon, aluminum, and silicon atoms in bicyclic organic compounds of the pentalene type was studied using theab initio RHF/6-31G^** and MP2(full)/6-31G^** methods. It was shown that the ability of the atom to form pentacoordinate structures increases on going from B to Al and from C to Si atom,i.e. as the number of the element of Groups IIIA and IVA of the periodic system increases. At the same time, the reverse tendencies are observed in the 2nd and 3rd periods of the periodic system,viz., the ability of the atom to form pentacoordinate structures increases on going from C to B and from Al to Si atom. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 7, pp. 1238–1245, July, 1999.     

Ab initio study of molecular structure and internal rotation in methyldicyanophosphine and methyldiisocyanophosphine. The Raman spectrum of methyldicyanophosphine and its interpretation with the use of scaling ofab initio force fields

The equilibrium geometric parameters and structures of transition states of internal rotation for the molecules of methyldicyanophospine MeP(CN)₂ and its isocyano analog MeP(NC)₂ were calculated by the RHF and MP2 methods with the 6–31G^* and 6–31G^** basis sets. At the MP2 level, the total energy of cyanide is −35 kcal mol⁻¹ lower than that of isocyanide and the barriers to internal rotation of methyl group for MeP(CN)₂ and MeP(NC)₂ are 2.2 and 2.7 kcal mol⁻¹, respectively. For both molecules, the one-dimensionalab initio potential functions of internal rotation approximated by a truncated Fourier series were used to determine the frequencies of torsional transitions by solving direct vibrational problems for a non-rigid model. The Raman spectrum of crystalline MeP(CN)₂ was recorded in the range 3500–50 cm⁻¹. The vibrational spectra of this compound were interpreted by scalingab initio force fields calculated by the RHF and MP2 method. The vibrational spectrum of methyldiisocyanophosphine was predicted with the use of the obtained scale factors. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 9, pp. 1703–1714, September, 1998.     

Structure and energy of the positively ionized water clusters

Geometry optimization of small (H₂O)_n⁺ clusters (n ≤ 4) at the UHF/4–31 + + G^** level indicates that the cations consist of two fragments: the OH radical and the H_2n−1 O⁺_n−1 ion. The latter can be considered as a thermodynamically stable combination of a distorted H₃O⁺ ion and (n−2) H₂O molecules. The H bond between the fragments becomes weaker with increasing cluster size. Extrapolation of the adiabatic ionization potentials calculated for the (H₂O)_n oligomers (n ≤ 4) at the MP2 level to an infinite cluster size provides the value of approximately 8.7 eV, which can be presumably necessary for the ionization of liquid water in a vacuum. © 1997 John Wiley & Sons, Inc.     

Molecular structures of acetylene derivatives of tin

The geometry and force fields of the bis(trimethylstannyl)acetylene molecule (a conformer withD _3d symmetry corresponding to a minimum of the total energy of the molecule) were calculated by the RHF and MP2(fc) methods. The effective core potential in SBK form with the optimized 31G^* valence basis set was employed in the case of Sn atoms. The 6–31G^** and 6–311G^** basis sets were used for carbon and hydrogen atoms. Vibrational spectra of the light and perdeuterated isotopomers of bis(trimethylstannyl)acetylene were interpreted using the procedure of scaling the quantum-chemical force fields. For Part 5, see Ref. 1. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 4, pp. 616–626, April, 2000.     

Ions with a strong symmetrical H bond in the HCl−BuiOH system by IR spectroscopy data

Ion-molecular interactions in the HCl−BuⁱOH system with different compositions (from neat isobutyl alcohol to 37 mol.% HCl) were studied by Multiple Attenuated Total Reflectance (MATR) IR spectroscopy at 30 °C. Proton disolvates (Buⁱ(H)O…H…O(H)Buⁱ)⁺ with strong symmetrical H bonds are formed upon the addition of HCl to BuⁱOH. At high concentrations of HCl (C ⁰ _HCl>33 mol.%), (Cl…H…Cl)⁻ ions are formed along with (BuⁱOH)₂H⁺. The spectra of positively and negatively charged proton disolvates were compared to those of similar ions in the HCl−PrⁱOH and HCl−MeOH systems. Translated fromIzvestiya Akademii Nauk. Seriya Khimicheskaya, No. 12, pp. 2496–2500, December, 1998.     

Conformational Analysis,Infrared, and Fluorescence Spectra of 1-Phenyl-1,2-propandione 1-oxime and Related Tautomers: Experimental and Theoretical Study

Summary. Conformational analysis and frequency calculation were achieved for 1-phenyl-1,2-propandione 1-oxime and its four tautomers: 1-nitroso-1-phenyl-1-propen-2-ol, 1-nitroso-1-phenyl-2-propanone, 2-hydroxy-1-phenyl-propenone oxime, and 3-nitroso-3-phenyl-propen-2-ol. Calculations were carried out at the Hartree–Fock (HF), Density Functional Theory (B3LYP), and the second-order M?ller–Plesset perturbation (MP2) levels of theory using 6-31G^* and 6-311G^** basis sets. Five conformers with no imaginary vibrational frequency were obtained by free rotations around three single bonds of 1-phenyl-1,2-propandione-1-oxime: Ph–C(NOH)C(O)CH₃, PhC(NOH)–C(O)CH₃, and PhC(N–OH)C(O)CH₃. Similarly, eight structures with no imaginary vibrational frequency were encountered upon rotations around three single bonds of 1-nitroso-1-phenyl-1-propen-2-ol: Ph–C(NO)C(OH)CH₃, PhC(N–O)C(OH)CH₃, and PhC(NO)C(–OH)CH₃. In the same manner, six minima were found through rotations around three single bonds of 1-nitroso-1-phenyl-2-propanone: Ph–CH(NO)C(O)CH₃, PhCH(–NO)C(O)CH₃, and PhCH(NO)–C(O)CH₃. Also, two minima were found through rotations around four single bonds of 2-hydroxy-1-phenyl-propenone oxime: Ph–C(NOH)C(OH)CH₂, PhC(N–OH)C(OH)CH₂, PhC(NOH)–C(OH)CH₂, and Ph-C(NOH)C(–OH)CH₂. Finally, two minima were found through rotations around four single bonds of 3-nitroso-3-phenyl-propen-2-ol: Ph–CH(NO)C(OH)CH₂, PhCH(–NO)C(OH)CH₂, PhCH(NO)–C(OH)CH₂, and PhCH(NO)C(–OH)CH₂. Interconversions within the above sets of conformers were probed through scanning (one and/or two dimensional), and/or QST3 techniques. The order of the stability of global minima encountered was: 1,2-propandione-1-oxime > 1-nitroso-1-phenyl-2-propanone > 1-nitroso-1-phenyl-1-propen-2-ol > 2-hydroxy-1-phenyl-propenone oxime > 3-nitroso-3-phenyl-propen-2-ol. Hydrogen bonding appears significant in tautomers of 1-nitroso-1-phenyl-1-propen-2-ol and 2-hydroxy-1-phenyl-propenone oxime. The CIS simulated λ_max for the first excited singlet state (S₁) of 1-phenyl-1,2-propandione 1-oxime is 300.4 nm, which was comparable to its experimental λ_max of 312.0 nm. The calculated IR spectra of 1-phenyl-1,2-propandione 1-oxime and its tautomers were compared to the experimental spectra.     

Tautomerism of xanthine: The second-order MØller-Plesset study

Four 9H and four 7H tautomers of DNA base xanthine were studied by the ab initio LCAO-MO method at the MP2/6-311G^**//HF/6-31G^** and MP2/6-31G^**//HF/6-31G^** approximations. All calculated structures are minima at the HF/6-31G^** potential energy surface with the dioxo 7H tautomer (A1) being the global minimum. The second most stable tautomer, dioxo-9H (B1) is by 9 kcal/mol less stable. For the A1 B1 transition the calculated MP2 energy gap corresponds to the equilibrium constant of 2 × 10^–7. Therefore, only the major tautomeric form A1 is predicted to be detectable in the gas phase. The 7H and 9H groups of tautomers are discussed separately. Within both groups, the dioxo form (A1-7H, B1-9H) is the most stable one and is succeeded by the 2-dihydroxy (A2, B2) form. However, while the energy difference between A1 and A2 is 10 kcal/mol, the energy difference between B1 a B2 is only 2 kcal/mol. The effect of polar environment was estimated by the SCRF method, using a spherical cavity, at the HF/6-31G^** level. These calculations did not change the gas phase stability order of the tautomers. However, the energy difference between A1 and B1 decreased from 9 kcal/mol at the HF/6-31G^** level to 4 kcal/mol at the SCRF HF/6-31G^** level.     

Theoretical refinements of theendo and theexo structures of tetraborane(8) carbonyl

The molecular structures of theendo (1a) andexo (1b) isomers of B₄H₈CO have been optimized at the ab initio MP2(Full)/6-31G^* level of theory. The agreement of the computed geometrical parameters with the recently published electron-diffraction (GED) data is very good, even though a number of geometrical constraints were applied in the experimental determination. The IGLO (individual gauge for localized orbitals)¹¹B NMR chemical shifts, calculated at the II//MP2/6-31G^* level, are also in accord with experiment. The formation of1a and1b by association of B₄H₈ and CO is computed to be exothermic by 22.8 and 22.2 kcal/mol, respectively, at the MP2(Full)/6-31G^*//MP2(Full)/6-31G^* + ZPE(6-31G^*) level of theory. The Lewis acid strength of B₄H₈ toward CO is comparable to that of BH₃.