Ab initio valence-bond calculations of H2O

The first and second bond dissociation energies for H₂O have been calculated in anab initio manner using a multistructure valence-bond scheme. The basis set consisted of a minimal number of non-orthogonal atomic orbitals expressed in terms of gaussian-lobe functions. The valence-bond structures considered properly described the change in the molecular system as the hydrogen atoms were individually removed to infinity. The calculated equilibrium geometry for the H₂O molecule has an O-H bond length of 1.83 Bohrs and an HOH bond angle of 106.5°. With 49 valence-bond structures the energy of H₂O at this geometry was ?76.0202 Hartrees. The calculated equilibrium bond length for the OH radical was 1.86 Bohrs and the energy, using the same basis set, was ?75.3875 Hartrees. After correction for zero point energies the calculated bond dissociation energies are: H₂O → OH + H, D₁=75.38 kcal/mole and OH → O+H, D₂=54.79 kcal/mole.     

AB initio SCF LCAO MO study of the HOH…Cl− hydrogen bond

Ab initio SCF LCAO MO calculations for the [H₂O…Cl]^? complex have been performed. The energy of the linear hydrogen bond has been found to be lower than the energy of the bifurcated one. The difference of the energies is about 3 kcal/mole. The calculated equilibrium distance between the oxygen and chlorine atoms equals 5.75 au. The interaction energy of the chlorine anion and the rigid water molecule amounts to ?19 kcal/mole. The optimization of the OH bond length in the complex (linear hydrogen bond) leads to an interaction energy of ?19.5 kcal/mole (the experimental value equals ?13.1 kcal/mole). As a result of the hydrogen bond formation the OH bond length increases by 0.08 au.     

Ab initio calculation of the lowest singlet and triplet states in CH2, CHF,CF2, and CHCH3

The lowest singlet and triplet states of the radicals CH₂, CHF, CF₂, and CHCH₃ have been investigated both in SCF and IEPA approximation (“independent electron pair approach” to account for electron correlation). The SCF calculations yield triplet ground states for CH₂, CHF, and CHCH₃, and a singlet ground state for CF₂. Electron correlation stabilizes the singlet state by about 14 kcal/mole with respect to the triplet for all four radicals leading to a singlet ground state also for CHF. The final triplet-singlet energy separations are 10, 6, ?11, ?47 kcal/mole for CH₂, CHCH₃, CHF, CF₂, respectively. Values for equilibrium bond angles, ionization potentials and bond energies are also given.     

SCF LCAO MO studies on the hydration of ions: The system F− · 2H2O

The energy surface of the dihydrated fluoride anion (F·2H₂O)^–1 is studied for a number of different geometry points near the equilibrium structure within the SCF LCAO MO framework, using an extended gaussian basis set to approximate the molecular wavefunctions. For the first and second hydration step of the fluoride anion the corresponding hydration energies are calculated to beB ₁ ^scf =24.1 kcal/mole andB ₂ ^SCF =20.8 kcal/mole (experimental measurements: 23.3 kcal/mole and 16.6 kcal/mole, respectively). The hydration energies and equilibrium bond distances obtained for the dihydrated fluoride anion (F·2H₂O)^– are compared with those found for the monohydrate (FHOH)^– and with corresponding results of the dihydrated lithium cation (Li · 2H₂O)⁺. The system (F·2H₂O)^– is taken as a very simple model to discuss some basic features of the hydration process of small ions and to study the influence of a negative ion on an adjacent hydrogen bond.We would like to thank our technical staff for valuable help in carrying out these calculations.     

The dynamics of reactions of O(3P) atoms with allene and methylacetylene

The reactions of O(³P) atoms with allene and methylacetylene: O+CH₂=C=CH₂

CO+C₂H₄,ΔH₁⁰ = ?119.4 kcal/mole, O+CH₃-C

CH

CO+C₂H₄,ΔH₂⁰ = ?117.8 kcal/mole were studied at 293 K with a CO laser resonant absorption and a discharge-flow GC-sampling method. The CO formed in reaction (1) was found to have a vibrational temperature of 5100 ± 100 K, compared with 2400 ± 200 K in (2). The good agreement between the observed CO vibrational distributions and those predicted by simple statistical models indicates that the reaction energies were completely randomized.The present results also showed unambiguously that CH₃CH, instead of C₂H₄, was produced initially in reaction (2).     

Nonadditivity of methyl group in single‐electron hydrogen bond of methyl radical‐water complex

The nonadditivity of methyl group in the single‐electron hydrogen bond of the methyl radical‐water complex has been studied with quantum chemical calculations at the UMP2/6‐311++G(2df,2p) level. The bond lengths and interaction energies have been calculated in the four complexes: CH₃? H₂O, CH₃CH₂? H₂O, (CH₃)₂CH? H₂O, and (CH₃)₃C? H₂O. With regard to the radicals, tert‐butyl radical forms the strongest hydrogen bond, followed by iso‐propyl radical and then ethyl radical; methyl radical forms the weakest hydrogen bond. These properties exhibit an indication of nonadditivity of the methyl group in the single‐electron hydrogen bond. The degree of nonadditivity of the methyl group is generally proportional to the number of methyl group in the radical. The shortening of the C···H distance and increase of the binding energy in the (CH₃)₂CH? H₂O and (CH₃)₃C? H₂O complexes are less two and three times as much as those in the CH₃CH₂? H₂O complex, respectively. The result suggests that the nonadditivity among methyl groups is negative. Natural bond orbital (NBO) and atom in molecules (AIM) analyses also support such conclusions. © 2008 Wiley Periodicals, Inc. Int J Quantum Chem, 2009     

The energy of the N→O bond in 3-nitroisoxazoline N-oxide

The enthalpy of combustion of 3-nitroisoxazoline has been determined as ΔH _c ^298.15 =?414±0.3 kcal/mole and that of 3-nitroisoxazoline N-oxide as ΔH _c ^298.15 =?406.6±0.5 kcal/mole. From the values for the heats of combustion and evaporation, the standard enthalpies of formation have been calculated and the energy of the N→O bond has been evaluated at 64±3 kcal/mole.     

Kinetics of the gas-phase reaction of acetone with iodine: Heat of formation and stabilization energy of the acetonyl radical

The kinetics of the gas-phase reaction CH₃COCH₃ + I₂ ? CH₃COCH₂I + HI have been measured spectrophotometrically in a static system over the temperature range 340–430°. The pressure of CH₃COCH₃ was varied from 15 to 330 torr and of I₂ from 4 to 48 torr, and the initial rate of the reaction was found to be consistent with \documentclass{article}\pagestyle{empty}\begin{document}$ {\rm CH}_3 {\rm COCH}_3 + {\rm I}^{\rm .} \stackrel{1}{\rightarrow}{\rm CH}_{\rm 3} {\rm COCH} + {\rm HI} $\end{document} as the rate-determining step. An Arrhenius plot of the variation of k₁ with temperature showed considerable scatter of the points, depending on the conditioning of the reaction vessel. After allowance for surface catalysis, the best line drawn by inspection yielded the Arrhenius equation, log [k₁/(M^?1 sec^?1)] = (11.2 ± 0.8) – (27.7 θ 2.3)/θ, where θ = 2.303 R T in kcal/mole. This activation energy yields an acetone C? H bond strength of 98 kcal/mole and δH (CH₃CO?H₂) radical = ?5.7 ± 2.6 kcal/mole. As the acetone bond strength is the same as the primary C? H bond strength in isopropyl alcohol, there is no resonance stabilization of the acetonyl radical due to delocalization of the radical site. By contrast, the isoelectronic allyl resonance energy is 10 kcal/mole, and reasons for the difference are discussed in terms of the π-bond energies of acetone and propene.     

Theoretical study of hydrogen‐bonded complexes of benzene with hydrides of astrochemical interest

Post Hartree–Fock and DFT calculations have been performed for studying the possibility for a benzene support to be linked to various hydrides through a quasi Bz···H? A bond. Interaction energy of compounds, including C? H bonds (CH₄, CH₃F, CH₂O, CHN, CHN? O), N? H bonds (NH₃, NH₂F, NHC, NHCO, NH₃O), and O? H bonds (OH₂, OHF, NCOH), were evaluated, taking basis set superposition error (BSSE) and zero point vibrational energy (ZPVE) corrections into account. Numerical convergence of results with respect to the ingredients included at different steps of theory (basis set, DFT functionals, correlation treatments, geometry optimization) was tested mainly on the example of the water adduct and, for comparison, the Bz···H₃O⁺ system containing a cation instead of a neutral molecule. A rather large range of adsorption energies is obtained, from about 1 kcal/mol for methane to more than 6 kcal/mol for cyanic acid, according to the acidic character of the adsorbed species in each family of Bz···H? A bonds. Some consequences for astrophysical problems involving PAHs in the interstellar medium are pointed out. © 2007 Wiley Periodicals, Inc. Int J Quantum Chem, 2008     

Heterocycles. CII. Barriers to rotation in some N′-Heteroaryl N,N-dimethylformamidines

The free energies of activation about the =CH? NMe₂ bond in N′-heteroaryl N,N-dimethylformamidines have been found in the range from 15.6 kcal/mole to 23 kcal/mole.     

Long‐range π‐type hydrogen bond in dimers CH2?HF,CH2O?H2O,and CH2O?NH3

Using four basis bets, (6‐311G(d,p), 6‐31+G(d,p), 6‐31++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for the dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. In contrast with the above three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD (T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Heterocycles. CLXIX. Barriers to rotation in some N,N-Dimethyl-N′-heteroaryl,N,N-Diethyl-N′-heteroaryl-substituted formamidines and N,N-Dimethyl-N′-heteroaryl-substituted acetamidines

The free energies of the rotational barriers, ΔG*, about ?CH? NMe₂ bond in N′-heteroaryl N,N-dimethylformamidines (A), about ?CH? NEt₂ bond in N-heteroaryl N,N-diethylformamidines (B), and about ?C(Me)? NMe₂ bond in N′-heteroaryl N,N-dimethylacetamidines (C) have been found to be in the range 17.5–20.1 kcal/mole for type A, 18.8–21.6 kcal/mole for type B and 13–14 kcal/mole or below for type C of compounds, respectively. The compounds of the types A and B exist in the forms IIa, IIIa, IV, V, and VI, while the compounds of the type C exist in the forms IIb and IIIb.     

Correlated one‐electron wave functions

Using four basis sets, 6‐311G(d,p), 6‐31+G(d,p), 6‐311++G(2d,2p), and 6‐311++G(3df,3pd), the optimized structures with all real frequencies were obtained at the MP2 level for dimers CH₂O? HF, CH₂O? H₂O, CH₂O? NH₃, and CH₂O? CH₄. The structures of CH₂O? HF, CH₂O? H₂O, and CH₂O? NH₃ are cycle‐shaped, which result from the larger bend of σ‐type hydrogen bonds. The bend of σ‐type H‐bond O…H? Y (Y?F, O, N) is illustrated and interpreted by an attractive interaction of a chemically intuitive π‐type hydrogen bond. The π‐type hydrogen bond is the interaction between one of the acidic H atoms of CH₂O and lone pair(s) on the F atom in HF, the O atom in H₂O, or the N atom in NH₃. By contrast with above the three dimers, for CH₂O? CH₄, because there is not a π‐type hydrogen‐bond to bend its linear hydrogen bond, the structure of CH₂O? CH₄ is a noncyclic shaped. The interaction energy of hydrogen bonds and the π‐type H‐bond are calculated and discussed at the CCSD(T)/6‐311++G(3df,3pd) level. © 2005 Wiley Periodicals, Inc. Int J Quantum Chem, 2005     

Molecular orbital theory of the hydrogen bond. VIII. Hydrogen bonding in H2OH2CO in relaxed singlet and triplet n → π* states

Ab initio SCF CI calculations with a minimal STO-3G basis set have been performed on the hydrogen bonded dimers in which H₂O is the proton donor to H₂CO in its relaxed singlet and triplet n→π^* states. Two dimers which are easily interconverted are found in the singet n→π^* state with hydrogen bond energies of 1.82 and 1.71 kcal/mole. The equilibrium dimer in the triplet state has a hydrogen bond energy of 2.97 kcal/mole. In both states, hydrogen bond formation occurs at the carbon atom. The structures of the dimers and the nature of the intermolecular surfaces in the regions of hydrogen bond formation are examined. Electron densities and distributions are also discussed.     

Conformational analysis of intramolecular bonded amino alcohols : The conformational free energies of some intramolecular OH ⋯ N hydrogen bonds

Equilibrium positions between intramolecular OH ? N hydrogen bonded and free OH forms of some 3-piperidinols, decahydroisoquinolinols, a decahydroquinolinol, lupinine and N-methyl-3-piperidinemethanol have been determined from dilute solution IR spectral data at 33°. Conformational free energies of the H-bonds (ΔG°_OH?N, attractive) have been calculated. The results suggest a linear relationship between the apparent value of ΔG°_OH?N, as defined by the method of calculation, and the strength of the OH ? N bond expressed as Δν, within the limits of 0·5 ± 0·2kcal/mole per 100 cm^?1, from Δν 90 to 350 cm^?1. For cis-decahydroisoquinoline (N-Me or N-H) systems, a 0·4 kcal/mole difference has been calculated between the two possible ring-fused conformations, in favor of the so-called steroid form. For the corresponding cis-decahydroqumoline equilibrium, a 0·8 kcal/mole difference has been calculated, in favor of the nonsteroid form.     

Atomic and molecular forms of oxygen on Ag(331). Theoretical analysis using the DFT method

Oxygen adsorption on Ag(331) is analyzed in a cluster approximation using the density functional theory (DFT) method. Adsorption centers (AC) for the bridge (S2) and three-center (S3) coordinations of oxygen are identified on the stepwise face Ag(331) and the Ag-O bond energies at these centers are calculated. For atomic adsorption, the Ag-O bond strength varies from 50 to 65 kcal/mole, depending on AC. The heat of molecular adsorption DH = 5 kcal/mole for S2(L1-L2) type AC. The molecule is oriented parallel to Ag(110) between the terraces with R(O-O) = 1.34 å Calculations showed that the ground state of the O₂Ag20(331) system is a triplet, but a part of spin density is delocalized on silver atoms, so that the spin density on oxygen ρ_s(O) = 0.46 (ρ_s = 1.0 for the free O₂ molecule). The energy of the singlet state is 9 kcal/mole greater than that of the ground state.     

Unconventional H‐bonds: SH···N interaction

Quantum calculations at the MP2/aug‐cc‐pVDZ level are used to analyze the SH···N H‐bond in complexes pairing H₂S and SH radical with NH₃, N(CH₃)₃, NH₂NH₂, and NH₂N(CH₃)₂. Complexes form nearly linear H‐bonds in which the S? H covalent bond elongates and shifts its stretching frequency to the red. Binding energies vary from 14 kJ/mol for acceptor NH₃ to a maximum of 22 kJ/mol for N(CH₃)₃ and N(CH₃)₂NH₂. Analysis of geometric, vibrational, and electronic data indicate that the SH···N interaction involving SH is slightly stronger than that in which the closed‐shell H₂S serves as donor. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

The remarkably invariant interaction energies of lithium first-row compounds with water and with ammonia

Solvation energies of lithium first-row compounds LiX (X ? H, Li, BeH, BH₂, CH₃, NH₂, OH, F) and of the lithium cation with the model solvents, water and ammonia, have been calculated ab inito (MP2/6-31 + G^*//6-31G^* with zero-point vibrational energy corrections at 3-21G//3-21G). The solvation energies are found to be remarkably constant: ?18.0 ± 1.2 and ?21.5 ± 1.3 kcal/mol for the hydrates and ammonia solvates, respectively. This independence on the nature of X is due largely to the ionic character of the LiX compounds (dipole moments 4.7–6.6 debye). The unexpectedly high solvation energies of the lithium molecule (?14.3 and ?17.8 kcal/mol, respectively) are due to the polarizability of Li₂. At the same level, the lithium cation has interaction energies with H₂O and NH₃ of ?34.1 and ?39.7 kcal/mol, respectively. For the hydrates of LiOH and LiF cyclic structures with hydrogen bonds and somewhat increased solvation energies also are described.     

A kinetic study of the thermal elimination of hydrogen fluoride from 1,2-difluoroethane. Determination of the bond dissociation energies D(CH2F?CH2F) and D(CH2F?H).

The kinetics of the thermal elimination of HF from 1,2-difluoroethane have been studied in a static system over the temperature range 734–820°K. The reaction was shown to be first order and homogeneous, with a rate constant of where θ = 2.303RT in kcal/mole. The A-factor falls within the normal range for such reactions and is in line with transition state theory; the activation energy is similarly consistent with an estimate based on data for the analogous reactions of ethyl fluoride and other alkyl halides. The above activation energy has been compared with values of the critical energy calculated from data on the decomposition of chemically activated 1,2-difluoroethane by the RRKM theory and the bond dissociation energy, D(CH₂F? CH₂F) = 88 ± 2 kcal/mole, derived. It follows from thermochemistry that ΔH_f⁰(CH₂F) = -7.8 and D(CH₂F? H) = 101 ± 2 kcal/mole. Bond dissociation energies in fluoromethanes and fluoroethanes are discussed.     

A study of aromatic three membered rings

The resonance energies (REs) of neutral three membered ring analogs of the cyclopropenyl cation, computed using block localized wave function (BLW) methods, reveal considerable variations. The RE's of cyclopropenes substituted with exocyclic double bonded groups C?X, (X = O, NH, CH₂, S, PH, SiH₂) increase with the electronegativity of X in the same row (SiH₂ < PH < S and CH₂ < NH < O). The extra cyclic resonance energies (ECREs) (an energetic measure of aromaticity based on comparisons with the RE's of acyclic models) of these derivatives range from +10.5 kcal/mol for cyclopropenone (X = O) (somewhat aromatic; the benzene ECRE is 29.3 kcal/mol) to ?2.4 kcal/mol (slightly antiaromatic) for X = SiH₂. Additional disubstitution of the C?C double bond by X′ groups (X′ = CH₃, NH₂, OH, SiH₃, PH₂, SH) increases the REs considerably, but has only small effects on the ECREs. Even the ECRE of deltic acid (X = O, X′ = OH) is only increased to +13.3 kcal/mol. The conclusion based on ECRE's, that all 12 of the three membered rings are only marginally aromatic/antiaromatic, is supported by the satisfactorily plot (R² = 0.92) of ECRE against values of NICS(0)_πzz (a superior nucleus chemical independent shift magnetic index of aromaticity), which range only from ?6.1 ppm (diatropic) for deltic acid (cf., ?35.5 ppm for benzene and ?14.2 ppm for the parent cyclopropenium ion) to +8.9 ppm (paratropic) for the silicon derivative, X = SiH₂, X′ = SiH₃. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011