Ab initio and density function theory computational studies of the CH4 + H → CH3 + H2 reaction

The activation barrier for the CH₄ + H → CH₃ + H₂ reaction was evaluated with traditional ab initio and Density Functional Theory (DFT) methods. None of the applied ab initio and DFT methods was able to reproduce the experimental activation barrier of 11.0-12.0 kcal/mol. All ab initio methods (HF, MP2, MP3, MP4, QCISD, QCISD(T), G1, G2, and G2MP2) overestimated the activation energy. The best results were obtained with the G2 and G2MP2 ab initio computational approaches. The zero-point corrected energy was 14.4 kcal mol⁻¹. Some of the exchange DFT methods (HFB) computed energies which were similar to the highly accurate ab initio methods, while the B3LYP hybrid DFT methods underestimated the activation barrier by 3 kcal mol⁻¹. Gradient-corrected DFT methods underestimated the barrier even more. The gradient-corrected DFT method that incorporated the PW91 correlational functional even generated a negative reaction barrier. The suitability of some computational methods for accurately predicting the potential energy surface for this hydrogen radical abstraction reaction was discussed.     

A Gaussian-2 ab initio study of [C2H5S] ions: III. H2 and CH4 eliminations from CH3CHSH and CH3SCH2

Gaussian-2 ab initio calculations were performed to examine the six modes of unimolecular dissociation of cis-CH₃CHSH⁺ (1⁺), trans-CH₃CHSH⁺ (2⁺), and CH₃SCH₂⁺ (3⁺): 1⁺→CH₃⁺+trans-HCSH (1); 1⁺→CH₃+trans-HCSH⁺ (2); 1⁺→CH₄+HCS⁺ (3); 1⁺→H₂+c-CH₂CHS⁺ (4); 2⁺→H₂+CH₃CS⁺ (5); and 3⁺→H₂+c-CH₂CHS⁺ (6). Reactions (1) and (2) have endothermicities of 584 and 496 kJ mol⁻¹, respectively. Loss of CH₄ from 1⁺ (reaction (3)) proceeds through proton transfer from the S atom to the methyl group, followed by cleavage of the C–C bond. The reaction pathway has an energy barrier of 292 kJ mol⁻¹ and a transition state with a wide spectrum of nonclassical structures. Reaction (4) has a critical energy of 296 kJ mol⁻¹ and it also proceeds through the same proton transfer step as reaction (3), followed by elimination of H_2. Formation of CH₃CS⁺ from 2⁺ (reaction (5)) by loss of H₂ proceeds through protonation of the methine (CH) group, followed by dissociation of the H₂ moiety. Its energy barrier is 276 kJ mol⁻¹. On both the MP2/6-31G* and QCISD/6-31G* potential-energy surfaces, the H₂ 1,1-elimination from 3⁺ (reaction (6)) proceeds via a nonclassical intermediate resembling c-CH₃SCH₂⁺ and has a critical energy of 269 kJ mol⁻¹.     

Thermochemistry of the higher chlorine oxides ClOx (x=3, 4) and Cl2Ox (x=3–7)†

Heats of formation for ClO₃, ClO₄, Cl₂O₃, Cl₂O₄, Cl₂O₅, Cl₂O₆ and Cl₂O₇ molecules are determined at the B3LYP, B3PW91, mPW1PW91 and B1LYP levels of the density functional theory employing a series of extended basis sets, and using Gaussian-3 model chemistries. Modified Gaussian-3 calculations, which employ accurate B3LYP/6-311+G(3d2f) molecular geometries and vibrational frequencies, were also performed. Heats of formation were calculated from both total atomization energies and isodesmic reaction schemes. The latter method in conjunction with Gaussian-3 models leads to the most reliable results. The best values at 298 K for ClO₃, ClO₄, Cl₂O₃ and Cl₂O₄ as derived from an average of G3//B3LYP and G3//B3LYP/6-311+G(3d2f) calculations are 43.1, 54.8, 31.7 and 37.4 kcal mol⁻¹. From calculations carried out at the G3(MP2)//B3LYP and G3(MP2)//B3LYP/6-311+G(3d2f) levels, heats of formation for Cl₂O₅, Cl₂O₆ and Cl₂O₇ are predicted to be 53.2, 52.2 and 61.5 kcal mol⁻¹. All best values are reproduced within 1 kcal mol⁻¹ by using mPW1PW91/6-311+G(3d2f) isodesmic energies. Enthalpy changes for relevant Cl–O bond fission reactions are reported. Comparisons with previous thermodynamics data are made.     

An ab initio and density functional theory study of the dimethylzinc-hydrogen selenide adduct: (CH3)2Zn:SeH2

The molecular structure (equilibrium geometry) and binding energy of the dimethylzinc (DMZn)-hydrogen selenide (H₂Se) adduct, (CH₃)₂Zn:SeH₂, have been computed with ab initio molecular orbital and density functional theory (DFT) methods and, where possible, compared with experimental results. The structure of the precursors DMZn and H₂Se are perturbed to only a small extent upon adduct formation. (CH₃)₂Zn:SeH₂ was found to be 3 kcal mol⁻¹ less stable than the precursors at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation, indicating that the (CH₃)₂Zn:SeH₂ adduct is unlikely to be a stable gas-phase species under chemical vapour deposition conditions. Further calculations at the B3LYP/6-311 + G(2d,p)//B3LYP/6-311 + G(2d,p) level of computation suggest that the 1:2 adduct species, (CH₃)₂Zn:(SeH₂)₂, is much less stable than the 1:1 adduct and consequently the precursors by 19 kcal mol⁻¹.     

Ab initio study of the two competitive channels HO2 + H → H2 + O2 and HO2 + H → H2O + O

Saddle point geometries and barrier heights have been calculated for the H abstraction reaction HO₂(²A″)+H(²S) → H₂(¹Σ⁺_g)+O₂(³Σ⁻_g) and the concerted H approach-O removing reaction HO₂ (²A″)+H(²S) → H₂O(¹A₁)+O(³P) by using SDCI wavefunctions with a valence double-zeta plus polarization basis set. The saddle points are found to be of C_s symmetry and the barrier heights are respectively 5.3 and 19.8 kcal by including size consistent correction. Moreoever kinetic parameters have been evaluated within the framework of the TST theory. So activation energies and the rate constants are estimated to be respectively 2.3 kcal and 0.4×10⁹ ℓ mol⁻¹ s⁻¹ for the first reaction, 20.0 kcal and 5.4.10⁻⁵ ℓ mol⁻¹ s⁻¹ for the second. Comparison of these results with experimental determinations shows that hydrogen abstraction on HO₂ is an efficient mechanism for the formation of H₂ + O₂, while the concerted mechanism envisaged for the formation of H₂O + O is highly unlikely.     

Ab initio MO study on the global minimum structure of C3H3 anion: propynl-1-yl versus allenyl anions

The geometrical structures of the C₃H⁻₃ anion are surveyed at the coupled-cluster doubles (CCD) level of theory with the aug-cc-pVDZ basis set. To clarify the CCD geometries, the stable two isomers -- propynl-l-yl 1 and allenyl 2 anions -- are further optimized at the coupled-cluster singles, doubles (triples) (CCSD(T)) level of theory both with the aug-cc-pVDZ and aug-cc-pVTZ basis sets. The final energies are calculated at the CCSD(T) and the complete active space self-consistent field (CASSCF) multi-reference internally contracted CI (MRCI) levels of theory with the aug-cc-pVTZ basis set. At the MRCI level of theory including both the corrections due to the cluster energies (MRCI+Q) and the zero-point vibrational energies, the allenyl anion 2 is about 1.3 kcal mol⁻¹ lower in energy than the propynl-l-yl anion 1. These results contrast with the previous theoretical estimates, where the propynl-l-yl anion 1 is 2-3 kcal mol⁻¹ lower in energy than the allenyl anion 2. The activation energies of the intramolecular hydrogen transfer in the 1 → 2 conversion reactions are 63.5 kcal mol⁻¹ at the MRCI+Q level of theory with the aug-cc-pVTZ basis set including the zero-point energy corrections. The adiabatic electron affinity of the planer propargyl (H₂CCCH) radical, which is the global minimum of the C₃H₃ radical, is calculated to be 0.976 eV (after correction for the zero-point energy changes) at the CCSD(T) level of theory with the aug-cc-pVTZ basis set. The present electron affinity is in fairly good agreement with the experimental one (0.893 eV) observed by Oakes and Ellison.     

Low-temperature heat capacity and thermodynamic properties of crystalline [Re2(Ala)4(H2O)8](ClO4)6 (Re = Eu, Er; Ala = alanine)

[Re₂(Ala)₄(H₂O)₈](ClO₄)₆ (Re=Eu, Er; Ala=alanine) were synthesized, and the low-temperature heat capacities of the two complexes were measured with a high-precision adiabatic calorimeter over the temperature range from 80 to 370 K. For [Eu₂(Ala)₄(H₂O)₈](ClO₄)₆, two solid–solid phase transitions were found, one in the temperature range from 234.403 to 249.960 K, with peak temperature 243.050 K, the other in the range from 249.960 to 278.881 K, with peak temperature 270.155 K. For [Er₂(Ala)₄(H₂O)₈](ClO₄)₆, one solid–solid phase transition was observed in the range from 270.696 to 282.156 K, with peak temperature 278.970 K. The molar enthalpy increments, ΔH_m, and entropy increments,ΔS_m, of these phase transitions, were determined to be 455.6 J mol⁻¹, 1.87 J K⁻¹ mol⁻¹ at 243.050 K; 2277 J mol⁻¹, 8.43 J K⁻¹ mol⁻¹ at 270.155 K for [Eu₂(Ala)₄(H₂O)₈](ClO₄)₆; and 4442 J mol⁻¹, 15.92 J K⁻¹ mol⁻¹ at 278.970 K for [Er₂(Ala)₄(H₂O)₈](ClO₄)₆. Thermal decompositions of the two complexes were investigated by use of the thermogravimetric (TG) analysis. A possible mechanism for the thermal decomposition is suggested.     

Theoretical study on the structure and stability of hydrogen-ion clusters Hn and Hn (n = 3, 5, 7, 9, 11, 13)

Ab initio SCF and CI calculations have been carried out for H_n⁺ and H_n⁻ (n = 3, 5, 7, 9, 11, 13) clusters with a double-zeta plus polarization basis set. The stabilization energy of negative ion clusters H_n⁻ is very small (less than 1 kcal) and their existence is critical. The structural difference between positive- and negative-ion clusters are discussed in terms of the bonding ability involved. While H_n⁺ is a charge-transfer complex, the stability of H_n⁻ comes mainly from the ion-induced-dipole attractions. The electron correlation effect on the structure and stability of these ion clusters is also discussed.     

Crystal structure of new hydroxide fluorides with isolated F anions: [H3N(CH2)6NH3]2M(F, OH)6(F,OH)·H2O (M = Al, In)

The synthesis and structural characterization of new hydroxo-fluorometallates [H₃N(CH₂)₆NH₃]₂M(F,OH)₇·H₂O (M=Al, In) are presented. Their preparation is achieved in solvothermal conditions by microwave or classical heating. The isotopic structures, determined by single crystal X-ray diffraction, are triclinic with the space group P–1. The structural arrangement can be described from isolated MX₆ (X=F, OH) octahedra connected by diprotonated diaminohexane via a complex network of hydrogen bonds. X⁻ anions and water molecules are found between the organic chains. A study by ¹⁹F NMR of the Al compound confirms a statistical occupancy of fluorine sites by hydroxyl groups and the occurrence of isolated F⁻ anions.     

MNDO study of the proton affinity of fluorinated formaldehydes and acetones

The MNDO molecular orbital method is employed to calculate the proton affinities of fluorinated formaldehydes and acetones. Agreement with experimentally reported proton affinities is good. In the acetone series a decrease in proton affinity is calculated for each successive fluorine substituent. The calculated strength of the intramolecular hydrogen bond in the protonated fluoro-formaldehydes and acetones is 0.6—2.7 kcal mol⁻¹, in good agreement with the experimental value of 2—3 kcal mol⁻¹ in the protonated fluoroacetones. Examination of the calculated charge distribution shows that the trends in proton affinity can be understood qualitatively both in terms of initial-state and final-state effects caused by the fluorine substituents. Protonation at the fluorine atom is less stable by about 25 kcal mol⁻¹ than protonation at the oxygen atom for monofluoroacetone.     

Ab initio molecular orbital study on three feasible mechanisms for substitution of the vinylic carbon in F2C=C(OMs)BMe3−

A substitution on 2,2-difluorovinylic carbon was investigated by using ab initio molecular orbital calculations. Three feasible mechanisms, which are the S_N1-like, the S_N2-type and the addition-elimination mechanisms, were ex- amined for a model borate, 2,2-difluoro-1-mesyloxyvinyl(trimethyl)borate. Four TSs were obtained depending on the position of Li⁺ around the vinylborate although activation energies in the gas phase are rather high (ca. 30–40 kcal mol⁻¹) in comparison with that expected from the experimental conditions. It was confirmed at the SCRF-IPCM calculations that the solvent effect reduces the acti- vation energy of one S_N2-type mechanism very much (4. l kcal mol⁻¹ at the B3LYP/6-31+G^*//RHF/6-31+G^/s* level of theory) while those for the other mechanisms do not change very much. Therefore, the S_N2-type mechanism is applicable to the substitution reaction observed for the vinylborate.     

Thermoynamic properties (Hm, CP,m, Vm, κT) of binary mixtures {x1,3-Dioxane + (1−x)Cyclohexane} at 298.15 K

Molar excess enthalpies H_m^E, isobaric heat capacities C_P,m^E, volumes V_m^E and isothermal compressibilities κT^E for the 1,3-dioxane(3DX) + cyclohexane mixture were measured at 298.15 K, in order to compare to those of the 1,4-dioxane(4DX) + cyclohexane mixture. H_m^E is endothermic and the maximum value about 1.5 kJ mol⁻¹ at x ≈ 0.45, and lower than that of the 4DX mixture by about 80 J mol⁻¹. V_m^E is positive over the whole concentration and the maximum value is about 0.85 cm³ mol⁻¹ at x ≈ 0.45, and lower than that of the 4DX mixture. The above results suggest the energetic unstabilization, resulting in the volume expansion in the mixture. C_P,m^E shows the characteristic W-shaped concentration dependence, which has maximum at x ≈ 0.45 and two minima at x ≈ 0.1 and 0.9. The maximum C_P,m^E value for 3DX mixture shifts toward the positive side, compared to that of 4DX mixture. κT^E were estimated from speeds of sound, densities, thermal expansion coefficients and isobaric heat capacities of the pure component liquids and the mixtures. The κT^E result shows the positive concentration dependence over the whole composition range. The 3DX mixture has the similar thermodynamic properties to the 4DX mixture, despite that 4DX is the nonpolar solvent and 3DX is the dipolar liquid. this means that there exists the local dipolar interaction between 4DX molecules, and the prevalence of “microheterogeneity” in the both mixtures.     

Structure and bonding studies of ClX-C2H4 (X = H or Br) van der Waals complexes

High-quality ab initio calculations at the MPn (n = 2, 3 or 4) levels for the π-donor complex formed by HCl and ethene, and the recently characterized ClBr-ethene complex are presented. Interaction energies were calculated with the inclusion of both basis-set superposition error and zero-point energy corrections, resulting in values of about 1 kcal mol⁻¹ for both complexes. The total charge densities for both complexes yielded molecular graphs indicative of the weak binding in these molecules, and correspondingly, calculated charges showed that the XCl moiety has an overall slightly negative charge, while the ethene moiety is slightly positive. Analysis of the Laplacian of the charge density showed that the geometry of the ClBr- C₂H₄ complex may be understood in terms of the VSEPR model as an example of an AXYE₃ molecule.     

Global minima of (C60)nCa, (C60)nF and (C60)nI clusters

Likely candidates for the lowest potential energy minima of (C₆₀)_nCa²⁺, (C₆₀)_nF⁻ and (C₆₀)_nI⁻ clusters are located using basin-hopping global optimisation. In each case, the potential energy surface is constructed using the Girifalco form for the C₆₀ intermolecular interaction, an averaged Lennard–Jones C₆₀–ion interaction, and a polarisation potential, which depends on the first few non-vanishing C₆₀ multipole polarisabilities. We find that the ions generally occupy the interstitial sites of a (C₆₀)_n cluster, the coordination shell being tetrahedral for Ca²⁺ and F⁻. The I⁻ ion has an octahedral coordination shell in the global minimum for (C₆₀)₆I⁻, however for 12  n  8 the preferred coordination geometry is trigonal prismatic.     

Infrared spectra of (CH3)2O and (CH3)2O + H2O at low temperature

Fourier transform infrared reflection spectroscopy (incidence angle of 5°) was used to characterize thin films of dimethyl ether (DME) and of mixtures containing water and DME between 10 and 160 K under a pressure of 10⁻⁷ mbar. Solid DME has two solid phases: an amorphous phase which is obtained below 65 K and a crystalline phase >65 K. From 90 K, DME begins to sublimate with surface binding energy of 20±2 kJ mol⁻¹. Vibrational spectrum of DME trapped in water ice remains nearly unchanged from 30 to 120 K. Between 120 and 130 K, a large part of DME is released and strong changes in the frequencies and the profile of the absorptions of DME are observed. This behavior suggests the formation of clathrate hydrate. Below 120 K, the trapped DME is hydrogen-bonded to water molecules.     

Heat capacity and thermodynamic properties of N-(2-cyanoethyl) aniline (C9H10N2)

The low temperature heat capacities of N-(2-cyanoethyl)aniline were measured with an automated adiabatic calorimeter over the temperature range from 83 to 353 K. The temperature corresponding to the maximum value of the apparent heat capacity in the fusion interval, molar enthalpy and entropy of fusion of this compound were determined to be 323.33 ± 0.13 K, 19.4 ± 0.1 kJ mol⁻¹ and 60.1 ± 0.1 J K⁻¹ mol⁻¹, respectively. Using the fractional melting technique, the purity of the sample was determined to be 99.0 mol% and the melting temperature for the tested sample and the absolutely pure compound were determined to be 323.50 and 323.99 K, respectively. A solid-to-solid phase transition occurred at 310.63 ± 0.15 K. The molar enthalpy and molar entropy of the transition were determined to be 980 ± 5 J mol⁻¹ and 3.16 ± 0.02 J K⁻¹ mol⁻¹, respectively. The thermodynamic functions of the compound [H_T − H_298.15] and [S_T − S_298.15] were calculated based on the heat capacity measurements in the temperature range of 83–353 K with an interval of 5 K.     

Ion chromatography of organic-rich natural waters from peatlands: I. Cl, NO2, Br, NO3, HPO4, SO4 and oxalate

Organic-rich natural waters from peat bogs in continental (Switzerland) and maritime (Shetland Islands, Scotland) areas were analysed for Cl⁻, NO₂⁻, Br⁻, NO₃⁻, HPO₄²⁻, SO₄²⁻ and oxalate using ion chromatography. These anions can be determined simultaneously in the surface and pore water samples from the continental bogs using a 250-μl injection loop. Using this loop, the detection limits were ca. 5 ng/g for the monovalent anions and SO₄²⁻ and 10 ng/g for HPO₄²⁻ and oxalate. An organics-removal cartridge (Dionex OnGuard P) was used to remove humic materials. These cartridges did not significantly affect the measured concentrations of anions in blind standards. Analyses of deionized water treated with these cartridges are not significantly different from those for untreated deionized water. For the maritime bogs, the relatively high concentrations of Cl⁻ (more than 100μ/g in many samples) and SO₄²⁻ (up to 50 μg/g) require two separate determinations for complete analyses. A 10-μl injection loop was used to determine Cl⁻, Br⁻ and SO₄²⁻. A 250-μl injection loop was used to measure NO₂⁻, NO₃⁻, HPO ₄²⁻ and oxalate. In each instance a Dionex OnGuard P cartridge was used to remove humic materials. In addition, a chloride-removal cartridge (Dionex OnGuard AG) was used to remove Cl⁻ when the larger injection loop was used. This cartridge has no significant effect on the measurement of HPO_4-²⁻ at concentrations of 20 ng/g. In each of the bog water chromatograms there were usually a number of unknown peaks. These are probably due mainly to organic anions.     

Structures and energetics of C3H6S isomers: a Gaussian-3 ab initio study

Twenty-two isomers/conformers of C₃H₆S^+√ radical cations have been identified and their heats of formation (ΔH_f) at 0 and 298 K have been calculated using the Gaussian-3 (G3) method. Seven of these isomers are known and their ΔH_f data are available in the literature for comparison. The least energy isomer is found to be the thioacetone radical cation (4⁺) with C_2v symmetry. In contrast, the least energy C₃H₆O^+√ isomer is the 1-propen-2-ol radical cation. The G3 ΔH_f298 of 4⁺ is calculated to be 859.4 kJ mol⁻¹, ca. 38 kJ mol⁻¹ higher than the literature value, ≤821 kJ mol⁻¹. For allyl mercaptan radical cation (7⁺), the G3 ΔH_f298 is calculated to be 927.8 kJ mol⁻¹, also not in good agreement with the experimental estimate, 956 kJ mol⁻¹. Upon examining the experimental data and carrying out further calculations, it is shown that the G3 ΔH_f298 values for 4⁺ and 7⁺ should be more reliable than the compiled values. For the five remaining cations with available experimental thermal data, the agreement between the experimental and G3 results ranges from fair to excellent.

Cation CH₃CHSCH₂^+√ (10⁺) has the least energy among the eleven distonic radical cations identified. Their ΔH_f298 values range from 918 to 1151 kJ mol⁻¹. Nevertheless, only one of them, CH₂=SCH₂CH₂^+√ (12⁺), has been observed. Its G3 ΔH_f298 value is 980.9 kJ mol⁻¹, in fair agreement with the experimental result, 990 kJ mol⁻¹.

A couple of reactions involving C₃H₆S^+√ isomers CH₂=SCH₂CH₂^+√ (12⁺) and trimethylene sulfide radical cation (13⁺) have also been studied with the G3 method and the results are consistent with experimental findings.     

Reaction of H with H2O2 as observed by optical absorption of perhydroxyl radicals or aliphatic alcohol radicals and of OH with H2O2. A pulse radiolysis study

Two new procedures were employed for studying the reaction of hydrogen atoms with hydrogen peroxide. The absorption in the UV-range was observed either for an acidic aqueous solution containing only hydrogen peroxide or for a similar solution but also containing an aliphatic alcohol. From the increase in absorption of various alcohol radicals, a rate constant of 3.5×10⁷ dm³ mol⁻¹ s⁻¹ was determined. In addition, the rate constant for the reaction of hydroxyl radicals with hydrogen peroxide was determined to be 3.0×10⁷ dm³ mol⁻¹ s⁻¹.     

Determination of the free activation energy for the breaking of the P---Ag bond in [(η-p-cymene)Ru(μ-pz)3Ag(PPh3)]

The Arrhenius equation corresponding to the process P---Ag+P*---Ag*→---P---Ag*+P*---Ag has been determined for [(η⁶-p-cymene)Ru(μ-pz)₃Ag(PPh₃)] (1) by complete line-shape analysis of the ³¹P NMR spectra between −40°C and +30°C. It has the form K = 10^11.8± e^{(−46±5 kJ mol−1/RT)}. The preexponential term, log A = 11.8 corresponds to a small activation entropy, whereas the activation energy, 46 kJ mol⁻¹ is comparable to those determined for other phosphorus—metal compounds.