Accurate ab initio predictions of ionization energies and heats of formation for the 2-propyl, phenyl, and benzyl radicals

The ionization energies (IEs) for the 2-propyl (2-C(3)H(7)), phenyl (C(6)H(5)), and benzyl (C(6)H(5)CH(2)) radicals have been calculated by the wave-function-based ab initio CCSD(T)/CBS approach, which involves the approximation to the complete basis set (CBS) limit at the coupled cluster level with single and double excitations plus quasiperturbative triple excitation [CCSD(T)]. The zero-point vibrational energy correction, the core-valence electronic correction, and the scalar relativistic effect correction have been also made in these calculations. Although a precise IE value for the 2-C(3)H(7) radical has not been directly determined before due to the poor Franck-Condon factor for the photoionization transition at the ionization threshold, the experimental value deduced indirectly using other known energetic data is found to be in good accord with the present CCSD(T)/CBS prediction. The comparison between the predicted value through the focal-point analysis and the highly precise experimental value for the IE(C(6)H(5)CH(2)) determined in the previous pulsed field ionization photoelectron (PFI-PE) study shows that the CCSD(T)/CBS method is capable of providing an accurate IE prediction for C(6)H(5)CH(2), achieving an error limit of 35 meV. The benchmarking of the CCSD(T)/CBS IE(C(6)H(5)CH(2)) prediction suggests that the CCSD(T)/CBS IE(C(6)H(5)) prediction obtained here has a similar accuracy of 35 meV. Taking into account this error limit for the CCSD(T)/CBS prediction and the experimental uncertainty, the CCSD(T)/CBS IE(C(6)H(5)) value is also consistent with the IE(C(6)H(5)) reported in the previous HeI photoelectron measurement. Furthermore, the present study provides support for the conclusion that the CCSD(T)/CBS approach with high-level energy corrections can be used to provide reliable IE predictions for C(3)-C(7) hydrocarbon radicals with an uncertainty of +/-35 meV. Employing the atomization scheme, we have also computed the 0 K (298 K) heats of formation in kJ/mol at the CCSD(T)/CBS level for 2-C(3)H(7)/2-C(3)H(7) (+) ,C(6)H(5)/C(6)H(5) (+), and C(6)H(5)CH(2)/C(6)H(5)CH(2) (+) to be 105.2/822.7 (90.0/806.4), 351.4/1148.5 (340.4/1138.8), and 226.2/929.0 (210.3/912.7), respectively. Comparing these values with the available experimental values, we find that the discrepancies for the 0 and 298 K heats of formation values are < or =2.6 kJ/mol for 2-C(3)H(7)/2-C(3)H(7) (+),< or =4.1 kJ/mol for C(6)H(5)/C(6)H(5) (+), and < or =3.2 kJ//mol for C(6)H(5)CH(2)C(6)H(5)CH(2) (+).     

Properties of pyridine and pyrazine from ab initio molecular orbital wave functions

Various molecular properties have been calculated for pyridine and pyrazine from Gaussian lobe ab initio SCF molecular wave functions. Values are compared with available experimental data. In general, agreement is satisfactory with the exception of the rather sensitive asymmetry parameter of the quadrupole coupling tensor. The distributions of total electronic charge, and of selected molecular orbitals have been displayed as plots of the charge density contours in two dimensions.     

Molecular structure of ethynylbenzene from electron diffraction and ab initio molecular orbital calculations

The molecular structure and benzene ring distortions of ethynylbenzene have been investigated by gas-phase electron diffraction and ab initio MO calculations at the HF/6-31G^* and 6-3G^** levels. Least-squares refinement of a model withC _2v, symmetry, with constraints from the MO calculations, yielded the following important bond distances and angles:r _g(C _i -C _o )=1.407±0.003 Å,r _g(C _o -C _m )=1.397±0.003 Å,r _g(C _m -C _p )=1.400±0.003 Å,r _g(Cr _i -CCH)=1.436 ±0.004 Å,r _g(C=C)=1.205±0.005 Å, C _o -C _i -C _o =119.8±0.4°. The deformation of the benzene ring of ethynylbenzene given by the MO calculations, including o-C_i-C_o=119.4°, is insensitive to the basis set used and agrees with that obtained by low-temperature X-ray crystallography for the phenylethynyl fragment, C₆H₅-CC-, in two different crystal environments. The partial substitution structure of ethynylbenzene from microwave spectroscopy is shown to be inaccurate in the ipso region of the benzene ring.     

In cation interactions with some organics: ab initio molecular orbital and density functional theory

Ab initio molecular orbital and density functional theory studies were undertaken to investigate the structural and energetic characteristics of complexes of In⁺ with several different organic molecules for the first time. HF, MP2, QCISD, and CCD levels of theory in ab initio MO as well as B3LYP, B3PW91 hybrid functionals in density functional theory were used. A valence TZ+P basis set with relativistic effective core potentials was used for the In atom while the 6-311++G(3d, 2p) basis set was utilized for all other atoms. Both closed-shell (H₂O, CH₄, CH₃OH, and C₆H₆) and open-shell (CH₃ and C₂H₃) molecules were considered for complexation with In⁺. In⁺ affinities of 21.5, 24.8, 28.6, 18.4, and 23.0 kcal/mol were obtained with the B3PW91 hybrid functional for H₂O, CH₃OH, C₆H₆, CH₃, and C₂H₃, respectively. The large values for the calculated affinities indicate the validity of our recent experimental detection of In⁺ ion attachment to some organic molecules.     

An ab initio molecular orbital study of the deprotonation of amines

The geometries of the amines NH₂X and amido anions NHX^?, where X = H, CH₃, NH₂, OH, F, C₂H, CHO, and CN have been optimized using ab initio molecular orbital calculations with a 4-31G basis set. The profiles to rotation about the N? X bonds in CH₃NH^?, NH₂NH^?, and HONH^? are very similar to those for the isoprotic and isoelectronic neutral compounds CH₃OH, NH₂OH, and HOOH. The amines with unsaturated bonds adjacent to the nitrogen atoms undergo considerable skeletal rearrangement on deprotonation such that most of the negative charge of the anion is on the substituent. The computed order of acidity for the amines NH₂X is X = CN > HCO > F ≈ C₂H > OH > NH₂ > CH₃ > H and for the reaction NHX^? + H⁺ → NH₂X the computed energies vary over the range 373–438 kcal/mol.     

An ab initio molecular orbital study of the protonation of amines

Ab initio molecular orbital calculations using an 8 ^s , 3 ^p ; 3 ^s Gaussian basis set, with contraction, have been used to study a series of primary amines XNH₂, where X = H, CH₃, OH, F, CN, CHO, and NO₂. The geometries of the corresponding ammonium ions have been optimised and the energy differences have been used to estimate relative proton affinities. The 1 ^s orbital energies for both the amines and ammonium ions, when corrected for the effects of charges on the other atoms in the molecule by use of an ESCA equation, give a good correlation with the computed charge on the nitrogen atom.     

An ab initio molecular dynamics and density functional theory study of the formation of phosphate chains from metathiophosphates

The reactions that occur between metathiophosphate (MTP) molecules are identified and examined through ab initio molecular dynamics simulations and static quantum chemical calculations at the density functional level of theory. The simulations show that certain types of MTPs can react to yield phosphate chains, while others only dimerize. These differences are rationalized in terms of reaction energies and the electronic structures of these molecules. In the reaction leading to the formation of phosphate chains, the reactive center, a tri-coordinate phosphorus atom, is continually regenerated. A polymerization mechanism linking MTPs to phosphate chains is developed on the basis of these results. This information sheds light on the underlying processes that may be responsible for the formation of phosphates under high-temperature conditions and may prove useful in the development of protocols for the rational synthesis of complex phosphate structures.     

The combined use of ab initio molecular orbital theory and matrix isolation infrared spectroscopy in the study of molecular interactions

The infrared spectra of a linear (C_s) and a cyclic (C_i) water dimer and a linear (C_s) and a cyclic (C_2h) ammonia dimer have been calculated, at the second order level of MØller-Plesset perturbation theory, using the 6-31G^** basis set. The calculated spectra have been compared with the spectra of water and ammonia isolated in nitrogen and argon matrices. In the case of water, theory predicts the linear to be the more stable isomer, and the observed spectrum can only be interpreted in terms of the linear structure. For ammonia, while the experimentally determined spectrum fits the calculated spectrum of the linear dimer fairly closely, the prediction of which is the more stable structure is still ambiguous.     

The structure of 1-chlorosilatrane: An ab initio molecular orbital and a density functional theory study

The molecular geometries of the 1-chloro-, 1-fluoro-, 1-methyl-, and 1-hydrogenosilatranes were fully optimized by the restricted Hartree-Fock (HF) method supplemented with 3-21G, 3-21G(d), 6-31G(d), and CEP-31G(d) basis sets; by MP2 calculations using 6-31G(d) and CEP-31G(d) basis sets; and by GGA-DFT calculations using 6-31G(d5) basis set with the aim of locating the positions of the local minima on the energy hypersurface. The HF/6-31G(d) calculations predict long (>254 pm) and the MP2/CEP calculations predicted short (∼225 pm) equilibrium Si(SINGLE BOND)N distances. The present GGA-DFT calculations reproduce the available gas phase experimental Si(SINGLE BOND)N distances correctly. The solid phase experimental results predict that the Si(SINGLE BOND)N distance is shorter in 1-chlorosilatrane than in 1-fluorosilatrane. In this respect the HF results show a strong basis set dependence, the MP2/CEP results contradict the experiment, and the GGA-DFT results in electrolytic medium agree with the experiment. The latter calculations predict that 1-chlorosilatrane is more polarizable than 1-fluorosilatrane and also support a general Si(SINGLE BOND)N distance shortening trend for silatranes during the transition from gas phase to polar liquid or solid phase. The calculations predict that the ethoxy links of the silatrane skeleton are flexible. Consequently, it is difficult to measure experimentally the related bond lengths and bond and torsion angles. This is the probable origin of the surprisingly large differences for the experimental structural parameters. On the basis of experimental analogies, ab initio calculations, and density functional theory (DFT) calculations, a gas phase equilibrium (r_e) geometry is predicted for 1-chlorosilatrane. The semiempirical methods predict a so-called exo minimum (at above 310 pm Si(SINGLE BOND)N distance); however, the ab initio and GGA-DFT calculations suggest that this form is nonexistent. The GGA-DFT geometry optima were characterized by frequency analysis. © 1996 by John Wiley & Sons, Inc.     

Molecular conformations of methyl formate and methyl vinyl ether from ab initio molecular orbital calculations

Ab initio molecular orbital theory with minimal and extended basis sets and a flexible rotor geometric model has been used to investigate the rotational potential surfaces of methyl formate and methyl vinyl ether. For both molecules, the most stable structures (IA and IIA, respectively) are planar cis; additional potential minima are found which correspond to planar trans structures (IB and IIB). The latter lie respectively about 4—8 and 1—2 kcal mol^?1 above the corresponding cis rotational isomers. Methyl rotational barriers have been determined for cis and trans structures of each molecule. For trans methyl formate, there is a slight but unexpected preference for an eclipsed arrangement of the methyl group.     

Approximate molecular orbital theory: the ese mo formalism. II. Direct comparison with ab initio results

The utility of the ESE MO formalism is examined in relation to ab initio results, when a typical small gaussian basis is employed. Good electronic structure predictions result from simplified schemes based on the ESE MO formalism, despite adverse features of the chosen basis set, for the diverse group of molecules FH, OH₂, NH₃, FCN, O₃ and OF₂.     

An ab initio SCF molecular orbital study on the conformation of serotonin and bufotenine

Ab initio SCF molecular orbital computations on the conformation of cationic serotonin and bufotenine indicate a preference for a perpendicular, or nearly so, arrangement of the ethylamine side chain with respect to the ring. The planar extended forms observed, among others, in the crystals of cationic indolealkylamines do not represent intrinsically stable conformations of these molecules. Their occurrence must be attributed to the effect of environmental forces.This work was supported by the A.T.P. N ° A. 655-2303 of the C.N.R.S.     

Conformational analysis of 4h-1,3,2-benzodioxaphosphorin 2-sulfides with ab initio molecular orbital calculations

Among some 2-substituted 4H-1,3,2-benzodioxaphosphorin 2-sulfides, the 2-OMe derivative has the sulfur atom mainly in the quasiequatorial position. On the other hand, the conformation in which the sulfur is quasiaxial is favored in the 2-NHMe and 2-Et derivatives. We have done ab initio molecular orbital calculations for the model compounds, and such a substituent effect was reproduced. Moreover, the effect was discussed in terms of the intramolecular electrostatic interactions and the degree of charge distribution.     

An ab initio molecular orbital study of intramolecular H-bonding: 1,3 propanediol

An ab initio molecular orbital calculation has been carried out for three different conformations of 1,3 propanediol, one of which permits intramolecular H-bond studied by ab initio quantum mechanical methods. The ΔE for H-bonding formation is compated to be 0.9 kcal/mole and the charge redistributions and molecular orbital energy changes are compared to those found in intermolecular H-bonds.     

A further study on closo boron hydrides B16H2−16 (D2) and B16H16Td) using ab initio molecular orbital theory

Ab initio molecular orbital calculations of doubly negative charged B₁₆H^2?₁₆(D₂) and neutral B₁₆H₁₆(Td) have been done at the HF/6-31G level. They are predicted to be chemically and kinetically stable by vibrational analyses on their respective energy hypersurface of the HF/6-31G level. The geometrical structure of the species B₁₆H²₁₆ (D₂) was discussed.     

ab initio molecular orbital and density functional analysis of acetylene + O2 reactions with CHEMKIN evaluation

A number of researchers have indicated that a direct reaction of acetylene with oxygen needs to be included in detailed reaction mechanisms in order to model observed flame speeds and induction times. Four pathways for the initiation of acetylene oxidation to chain propagation are considered and the rate constants are compared with values used in the mechanisms:

• 1 ³O₂ + HCCH to triplet adduct and reaction on the triplet surface
• 2 ³O₂ + HCCH to triplet adduct, conversion of triplet adduct to singlet adduct via collision in the reaction environment, with further reaction of the singlet adduct
• 3 ¹O₂ + HCCH to singlet adduct
• 4 Isomerization of HCCH to vinylidene and then vinylidene insertion reaction with ³O₂

Elementary reaction pathways for oxidation of acetylene by addition reaction of O₂(³Σ) on the triplet surface are analyzed. ab initio molecular orbital and density functional calculations are employed to estimate the thermodynamic properties of the reactants, transition states, and products in this system. Acetylene oxidation reaction over the triplet surface is initiated by addition of molecular oxygen, O₂(³Σ), to a carbon atom, forming a triplet peroxy‐ethylene biradical. The reaction path to major products, either two formyl radicals or glyoxal radical plus hydrogen atom, involves reaction through three transition states: O₂(³Σ) addition to acetylene (TS1), peroxy radical addition at the ipso‐carbon to form a dioxirane (TS2), and cleavage of O O bond in a three‐member ring (TS3). Single‐point QCISD(T) and B3LYP calculations with large basis sets were performed to try to verify barrier heights on important transition states. A second pathway to product formation is through spin conversion of the triplet peroxy‐ethylene biradical to the singlet by collision with bath gas. Rapid ring closure of the singlet peroxy‐ethylene biradical to form a four‐member ring is followed by breaking of the peroxy bond to form glyoxal, which further dissociates to either two formyl radicals or a glyoxal radical plus hydrogen atom. The overall forward rate constant through this pathway is estimated to be k_f = 2.21 × 10⁷ T^1.46e^{−33.1(kcal/mol)/RT}. Two additional pathways from the literature, HCCH + O₂(¹Δ) and pressure‐dependent isomerization of acetylene to vinylidene and then vinylidene reaction with O₂(³Σ), are also evaluated for completeness. CHEMKIN modeling on each of the four proposed pathways is performed and concentration profiles from these reactions are evaluated at 0.013 atm and 1 atm over 35 milliseconds. Through reaction on the triplet surface is evaluated to be not important. Formation of the triplet adduct with conversion (via collision) to a singlet and the vinylidene paths show similar and lower rates than those used in mechanisms, respectively. Our implementation of the HCCH + O₂(¹Δ) pathway of Benson suggests the need to include: (i) reverse reaction, (ii) barriers to further reaction of the initial adduct plus (iii) further evaluation of the O₂(¹Δ) addition barrier. The pathways from triplet adduct with conversion to singlet and from vinylidene are both recommended for initiation of acetylene oxidation. © 2000 John Wiley & Sons, Inc. Int J Chem Kinet 32: 623–641, 2000     

Ab initio molecular orbital study of the insertion of H2 into POSS compounds

The insertion of one and two H₂ molecules into polyhedral oligomeric silsesquioxanes (POSS) was investigated as a function of the size of the cage, using both Hartree-Fock (HF) and second order perturbation theory (MP2) methods. Also investigated was the same reaction into the heavier groups 4 and 14 metal-substituted POSS (metallasilsesquioxanes) such as Ge-POSS, Si/Ge-mixed POSS, and Ti- and Zr-POSS. The properties of these species in comparison with those of POSS are discussed.     

An ab initio molecular orbital study of the C2F2 species: The difluorovinylidene → difluoroacetylene rearrangement

An ab initio MO study, using medium size Gaussian basis sets has been made of vinylidene carbene, acetylene and the isomeric cyclic intermediate presumed to represent the transition state of their interconversion reaction, along with their perfluorinated analogs. The total energies of the acetylenes are lower than the vinylidenes and the estimated activation energy of their interconversion is considerably higher for the fluorinated molecules. The results are in line with experimental observations.     

The geometric and electronic structures of oxocarbons. An ab initio molecular orbital study

Ab initio molecular orbital calculations with the STO-3G and 4-31G basis sets have been carried out for the neutral oxocarbons C_nO_n (n = 3, 4, 5, 6 and 7), the dianions C_nO_n^2- (n = 3, 4, 5, 6 and 7), the monoanions C_nO_nH^? (n = 3 and 4) and the related acids C_nO_nH₂ (n = 3 and 4). Fully optimised geometries have been obtained for all species. The geometries, stabilities and acidities are discussed.