Ab initio calculations for ground states of CH2Li− and CH2Be

The ground state of CH₂Li^? and CH₂Be molecules has been investigated by an SCF calculation using a contracted Gaussian basis set. Only for the second system a bound state with respect to the ground states of the molecular fragments has been found.     

Simplified methods for ab initio calculations. The valence states of CH2 and CH

Three simplifying methodds are discussed and applied to the four lowest valence states of CH₂(³B₁, ¹A₁, ¹B₁ and ¹Σ(¹A)) and to the two lowest of CH(²A₁ and ²∏_u(²B₁)). These methods concern: (1) the development of polarization functions for Gaussian-lobe basis sets by least-square fitting of numebrical multiconfigurational atomic fuinctions (this approach is tested also on (C₂H₂, (2) the use of intermediate Hamiltonians to calculate avoided crossings between potential hypersurfaces, and (3) thecalculation of correlation energies using an atoms-in-molecule approach. The calculated equilibrium geometries of the CH₂ States are within 0.02 Å and 5° of available experimental data. The calculated term values and ionization potentials, T_e(¹A₁ = 0.35 eV, T_e (¹B₁) = 1.22 eV, T_e (¹Σ(¹A)) 2.48 eV, I.P. (²A₁) = 10.39 eV and I. P. (²∏_u(²B₁)) = 10.58 eV, are in agreement with some recent theoretical studies, and are very close to existing experimental information.     

Nonempirical SCF-LCAO-MO calculations for (CH)+

Nonempirical molecular orbital calculations, using a minimal basis set of Slater-type atomic orbitals, are reported for the ground state of the (CH)⁺ ion. The C-H bond length calculated using optimized exponents is in excellent agreement with the experimental value, whereas the distances predicted using the conventional free atom exponents are too large by up to 17%.     

Multichannel RRKM-TST and CVT rate constant calculations for reactions of CH2OH or CH3O with HO2

The kinetics and mechanism of the gas-phase reactions between hydroxy methyl radical (CH(2)OH) or methoxy radical (CH(3)O) with hydroproxy radical (HO(2)) have been theoretically investigated on their lowest singlet and triplet surfaces. Our investigations indicate the presence of one deep potential well on the singlet surface of each of these systems that play crucial roles on their kinetics. We have shown that the major products of CH(2)OH + HO(2) system are HCOOH, H(2)O, H(2)O(2), and CH(2)O and for CH(3)O + HO(2) system are CH(3)OH and O(2). Multichannel RRKM-TST calculations have been carried out to calculate the individual rate constants for those channels proceed through the formation of activated adducts on the singlet surfaces. The rate constants for direct hydrogen abstraction reactions on the singlet and triplet surfaces were calculated by means of direct-dynamics canonical variational transition-state theory with small curvature approximation for the tunneling.     

Extended group function calculations for H2O,NH3 and CH4

Wave functions expressed as antisymmetrized products of strongly orthogonal geminals have been evaluated for H₂O, NH₃ and CH₄. The geminals have been expressed as linear combinations of 2 × 2 detors constructed with localized SCF -MO 's. Several ground state observables have been computed together with the electric polarizabilities and magnetic susceptibilities. In addition, a configuration interaction calculation limited to all possible double group excitations has been carried out.     

Thermochemistry and accurate quantum reaction rate calculations for H2/HD/D2 + CH3

Accurate quantum-mechanical results for thermodynamic data, cumulative reaction probabilities (for J = 0), thermal rate constants, and kinetic isotope effects for the three isotopic reactions H2 + CH3 --> CH4 + H, HD + CH3 --> CH4 + D, and D2 + CH3 --> CH(3)D + D are presented. The calculations are performed using flux correlation functions and the multiconfigurational time-dependent Hartree (MCTDH) method to propagate wave packets employing a Shephard interpolated potential energy surface based on high-level ab initio calculations. The calculated exothermicity for the H2 + CH3 --> CH4 + H reaction agrees to within 0.2 kcal/mol with experimentally deduced values. For the H2 + CH3 --> CH4 + H and D2 + CH3 --> CH(3)D + D reactions, experimental rate constants from several groups are available. In comparing to these, we typically find agreement to within a factor of 2 or better. The kinetic isotope effect for the rate of the H2 + CH3 --> CH4 + H reaction compared to those for the HD + CH3 --> CH4 + D and D2 + CH3 --> CH(3)D + D reactions agree with experimental results to within 25% for all data points. Transition state theory is found to predict the kinetic isotope effect accurately when the mass of the transferred atom is unchanged. On the other hand, if the mass of the transferred atom differs between the isotopic reactions, transition state theory fails in the low-temperature regime (T < 400 K), due to the neglect of the tunneling effect.     

Gradient optimization of polarization exponents inab initio MO calculations on H2SO → HSOH and CH3SH → CH2SH2

Summary Analytical gradients were used to optimize the polarization function exponents in the 6-31G(d) and 6-31G(d, p) basis sets for the reactants, transition structures and products in the reactions H₂SO HSOH and CH₃SH CH₂SH₂. The optimizedd exponents on the heavy atoms change by ±10% in the course of the reactions and depend on the bonding of the heavy atoms. Thep exponents on the hydrogens change by as much as a factor of 5 and depend on the element to which the hydrogen is bonded and its valency. The effect of exponent optimization on the relative energies is small (±3 kcal/mol). With the 6-31G(d, p) basis set, optimization of the polarization exponents can make some of the bonds significantly more polar, as judged by the Mulliken charges.     

Accurate potential energy surface and quantum reaction rate calculations for the H+CH4-->H2+CH3 reaction

Calculations for the cumulative reaction probability N(E) (for J=0) and the thermal rate constant k(T) of the H+CH(4)-->H(2)+CH(3) reaction are presented. Accurate electronic structure calculations and a converged Shepard-interpolation approach are used to construct a potential energy surface which is specifically designed to allow the precise calculation of k(T) and N(E). Accurate quantum dynamics calculations employing flux correlation functions and multiconfigurational time-dependent Hartree wave packet propagation compute N(E) and k(T) based on this potential energy surface. The present work describes in detail the various convergence test performed to investigate the accuracy of the calculations at each step. These tests demonstrate the predictive power of the present calculations. In addition, approximate approaches for reaction rate calculations are discussed. A quite accurate approximation can be obtained from a potential energy surface which includes only interpolation points on the minimum energy path.     

Vibrational calculations in formaldehyde: the CH stretch system

An alternative procedure for the calculation of highly excited vibrational levels in S0 formaldehyde was developed to apply to larger molecules. It is based on a new set of symmetrized vibrational valence coordinates. The fully symmetrized vibrational kinetic energy operator is derived in these coordinates using the Handy expression [Molec. Phys. 61, 207 (1987)]. The potential energy surface is expressed as a fully symmetrized quartic expansion in the coordinates. We have performed ab initio electronic computations using GAMESS to obtain all force constants of the S0 formaldehyde quartic force field. Our large scale vibrational calculations are based on a fully symmetrized vibrational basis set, in product form. The vibrational levels are calculated one by one using an artificial intelligence search/selection procedure and subsequent Lanczos iteration, providing access to extremely high vibrational energies. In this work special attention has been given to the CH stretch system by calculating the energies up to the fifth CH stretch overtone at ∼16000 cm⁻¹, but the method has also been tested on two highly excited combination levels including other lower frequency modes.      

AB Initio calculations and vibrational spectra of CH3SC≡CH and CH3SC≡CSCH3

Methylthio- (MTE) and bis-methylthioethyne (BMTE) molecules are calculated by the SCF MO method (geometry optimization, basis set 6–31G*/MP2). The calculated internal rotation barriers of methyl groups are 7.12 kJ/lmole for MTE and 12.86 kJ/mole for BMTE (both groups are simultaneously rotated). The s-gosh-orientation of the thiomethyl fragments corresponds to a stable conformation of BMTE. The estimated values of the s-cis- and s-trans-barriers of mutual rotation of SCH₃ groups about the axis of the C≡C bond are 13.61 and 12.54 kJ/mole, respectively. Conformationally sensitive MOs and vibration frequencies are established. An analysis of the experimental IR absorption and Raman spectra and the calculated vibrational spectrum makes it possible to conclude that in the liquid phase the BMTE molecules also have an s-gosh-conformation. Translated fromZhumal Strukturnoi Khimii, Vol. 39, No. 4, pp. 602–609, July–August, 1998.     

Synthesis and Characterization of （CH3CH2CH2CH2NH3）2SnBr6

1 INTRODUCTION Recently, the researches on inorganic-organic hy-brid compounds represent an advanced field in mate-rial science[1]. At the molecular level, the combina-tion of two extremely different components providesan avenue to design new hybrid materials as well asthe ability to modulate properties of one or more ofthe components[2～6]. Some attractive properties, suchas efficient luminescence[2～4], ideal thermal and me-chanical stability, interesting magnetic[5], non-linearoptical[…     

Thermodynamic calculations in the system CH4-H2O and methane hydrate phase equilibria

Using the Gibbs function of reaction, equilibrium pressure, temperature conditions for the formation of methane clathrate hydrate have been calculated from the thermodynamic properties of phases in the system CH4-H2O. The thermodynamic model accurately reproduces the published phase-equilibria data to within +/-2 K of the observed equilibrium boundaries in the range 0.08-117 MPa and 190-307 K. The model also provides an estimate of the third-law entropy of methane hydrate at 273.15 K, 0.1 MPa of 56.2 J mol(-1) K(-1) for 1/nCH4.H2O, where n is the hydrate number. Agreement between the calculated and published phase-equilibria data is optimized when the hydrate composition is fixed and independent of the pressure and temperature for the conditions modeled.     

Secondary kinetics of methanol decomposition: theoretical rate coefficients for 3CH2 + OH, 3CH2 + 3CH2, and 3CH2 + CH3

Direct variable reaction coordinate transition state theory (VRC-TST) rate coefficients are reported for the (3)CH(2) + OH, (3)CH(2) + (3)CH(2), and (3)CH(2) + CH(3) barrierless association reactions. The predicted rate coefficient for the (3)CH(2) + OH reaction (approximately 1.2 x 10(-10) cm(3) molecule(-1) s(-1) for 300-2500 K) is 4-5 times larger than previous estimates, indicating that this reaction may be an important sink for OH in many combustion systems. The predicted rate coefficients for the (3)CH(2) + CH(3) and (3)CH(2) + (3)CH(2) reactions are found to be in good agreement with the range of available experimental measurements. Product branching in the self-reaction of methylene is discussed, and the C(2)H(2) + 2H and C(2)H(2) + H2 products are predicted in a ratio of 4:1. The effect of the present set of rate coefficients on modeling the secondary kinetics of methanol decomposition is briefly considered. Finally, the present set of rate coefficients, along with previous VRC-TST determinations of the rate coefficients for the self-reactions of CH(3) and OH and for the CH(3) + OH reaction, are used to test the geometric mean rule for the CH(3), (3)CH(2), and OH fragments. The geometric mean rule is found to predict the cross-combination rate coefficients for the (3)CH(2) + OH and (3)CH(2) + CH(3) reactions to better than 20%, with a larger (up to 50%) error for the CH(3) + OH reaction.     

Dative sigma- and pi-bonding in boron-nitrogen compounds: molecular structures of (CH3)2NB(CH3)N(CH3)B(CH3)2 and [(CH3)2B]2NN(CH3)2 determined by gas electron diffraction and quantum chemical calculations

The molecular structures of dimethylamino[(dimethylboryl)methylamino]methylborane, Me2NBMeNMeBMe2 (1) and 1,1-bis(dimethylboryl)-2,2-dimethylhydrazine, (Me2B)2NNMe2 (2) have been determined by gas electron diffraction (GED), density functional theory calculations at the B3PW91/6-311++G** level and ab initio calculations at the MP2/6-311++G** level. 1 adopts an open structure similar to that of the isoelectronic hydrocarbon molecule permethylbutadiene; the central B-N bond distance at 148.0/149.3(7) pm (MP2/GED) corresponds to a single covalent N--B bond distance, the two terminal distances, 140.9/140.5(4) pm and 141.8/141.3(4) pm, correspond to the distance between N and B atoms joined by a covalent sigma-bond and a dative pi-bond. A closed form where the establishment of a dative bond between the terminal N and B atoms has led to the formation of a four-membered ring also corresponds to a minimum on the potential energy surface, but the energy is calculated to be 14.3 kJ mol(-1) higher at the MP2 level. This structure is also incompatible with the GED data. 2 adopts a structure in which a dative sigma-bond between the dimethylamino N atom and one of the boron atoms has led to the formation of a three-membered N(2)B ring. The dative sigma-bond distance is 165.5/164.0(13) pm, the two other bond distances in the ring are N-B=150.6/148.9(9) pm corresponding to a covalent sigma-bond and N-N=145.1/145.4(3) pm. The terminal B--N distance 139.6/138.9(9) pm is consistent with a covalent sigma-bond augmented by a dative pi-bond. An open Y-shaped structure also corresponds to a minimum on the potential energy surface, but the energy is 18.7 kJ mol(-1) higher (MP2) and it is incompatible with the GED data.     

Kinetic and mechanistic study of the reactions of atomic chlorine with CH3CH2Br, CH3CH2CH2Br, and CH2BrCH2Br

A laser flash photolysis-resonance fluorescence technique has been employed to investigate the reactions of atomic chlorine with three alkyl bromides (R-Br) that have been identified as short-lived atmospheric constituents with significant ozone depletion potentials (ODPs). Kinetic data are obtained through time-resolved observation of the appearance of atomic bromine that is formed by rapid unimolecular decomposition of radicals generated via abstraction of a β-hydrogen atom. The following Arrhenius expressions are excellent representations of the temperature dependence of rate coefficients measured for the reactions Cl + CH(3)CH(2)Br (eq 1 ) and Cl + CH(3)CH(2)CH(2)Br (eq 2 ) over the temperature range 221-436 K (units are 10(-11) cm(3) molecule(-1) s(-1)): k(1)(T) = 3.73?exp(-378/T) and k(2)(T) = 5.14?exp(+21/T). The accuracy (2σ) of rate coefficients obtained from the above expressions is estimated to be ±15% for k(2)(T) and +15/-25% for k(1)(T) independent of T. For the relatively slow reaction Cl + CH(2)BrCH(2)Br (eq 3 ), a nonlinear ln k(3) vs 1/T dependence is observed and contributions to observed kinetics from impurity reactions cannot be ruled out; the following modified Arrhenius expression represents the temperature dependence (244-569 K) of upper-limit rate coefficients that are consistent with the data: k(3)(T) ≤ 3.2 × 10(-17)T(2)?exp(-184/T) cm(3) molecule(-1) s(-1). Comparison of Br fluorescence signal strengths obtained when Cl removal is dominated by reaction with R-Br with those obtained when Cl removal is dominated by reaction with Br(2) (unit yield calibration) allows branching ratios for β-hydrogen abstraction (k(ia)/k(i), i = 1,2) to be evaluated. The following Arrhenius-type expressions are excellent representations of the observed temperature dependences: k(1a)/k(1) = 0.85?exp(-230/T) and k(2a)/k(2) = 0.40 exp(+181/T). The accuracy (2σ) of branching ratios obtained from the above expressions is estimated to be ±35% for reaction 1 and ±25% for reaction 2 independent of T. It appears likely that reactions 1 and 2 play a significant role in limiting the tropospheric lifetime and, therefore, the ODP of CH(3)CH(2)Br and CH(3)CH(2)CH(2)Br, respectively.     

High-level ab initio calculations on CH+2(2A1) + PO(2II) reactions

We present ab initio calculations carried out in the framework of the G 2 theory on the singlet and triplet potential energy surfaces corresponding to the gas-phase between CH⁺₂ and PO. The global minimum of both potential energy surfaces is a cyclic singlet-state cation. Oxygen attachment of PO to CH⁺₂ in a triplet configuration is accompanied by a P(SINGLEBOND)O bond fission, with the result that the corresponding global minimum is an ion-dipole complex between P⁺(³P) and formaldehyde. This is also consistent with the fact that our results predict the formation of formaldehyde to be highly exothermic, either as a neutral or as radical cation. Both charge-transfer processes yielding CH₂(³B₁) or CH₂(¹A₁) are also exothermic. The formation of other carbon and oxygen containing species are endothermic. © 1996 John Wiley & Sons, Inc.     

The elusive methene ylides CH2CIH,CH2FH and CH2OH2

The possible existence of the ylides CH₂CIH, CH₂₅FH and CH₂OH₂ as stable neutral species in the gas phase has been investigated by the neutralization–reionization (NR) mass spectrometry of their radical cations using a double-focusing mass spectrometer of reversed geometry. The experiments were, for the most part, performed under single-collision conditions with Xe as the neutralization target gas and He and O₂ as reionization agents. For each ylidion a peak was observed in their NR mass spectrum which indicated that the neutral ylide had apparently been produced. However for CH₂FH⁺˙ and CH₂OH₂⁺˙ the m/z 34 and m/z 32 peaks, respectively, were attribut-able to interferences from the natural isotopic abundance of ions of lower mass. For CH₂CIH⁺˙, the NR recovery signal was found to arise from the presence of CH₃CI⁺˙ as an impurity in the ylidion flux. This was proved by examination of the collisional activation mass spectra of the [C, H₃, CI]⁺˙ ions produced in the NR mass spectra of the conventional ions and ylidions, an experiment performed using a triple-sector mass spectrometer.     

Thermal decomposition of di-t-butyl peroxide in the presence of (CH3)2CCH2: Reactions of CH3, (CH3)2ĊCH2CH3, and (CH3)2ĊCH2C(CH3)2CH2CH3 radicals

A study of the reaction initiated by the thermal decomposition of di-t-butyl peroxide (DTBP) in the presence of (CH₃)₂C?CH₂ (B) at 391–444 K has yielded kinetic data on a number of reactions involving CH₃ (M·), (CH₃)₂CCH₂CH₃ (MB·) and (CH₃)₂?CH₂C(CH₃)₂CH₂CH₃ (MBB·) radicals. The cross-combination ratio for M· and MB· radicals, rate constants for the addition to B of M· and MB· radicals relative to those for their recombination reactions, and rate constants for the decomposition of DTBP, have been determined. The values are, respectively, where θ = RT ln 10 and the units are dm^3/2 mol^?1/2 s^?1/2 for k₂/k and k₉/k, s^?1 for k₀, and kJ mol^?1 for E. Various disproportionation-combination ratios involving M·, MB·, and MBB· radicals have been evaluated. The values obtained are: Δ₁(M·, MB·) = 0.79 ± 0.35, Δ₁(MB·, MB·) = 3.0 ± 1.0, Δ₁(MBB·, MB·) = 0.7 ± 0.4, Δ₁(M·, MBB·) = 4.1 ± 1.0, Δ₁(MB·, MBB·) = 6.2 ± 1.4, and Δ₁(MBB·, MBB·) = 3.9 ± 2.3, where Δ₁ refers to H-abstraction from the CH₃ group adjacent to the center of the second radical, yielding a 1-olefin. © 1994 John Wiley & Sons, Inc.     

Ab initio molecular calculations with pseudopotentials: calculations of double-zeta quality on BeH2, BH3, CH4, and C2H6

A new form of pseudopotential for applicaiton in ab initio molecular calculations is described. A method for determining pseudopotential parameters is suggested and pseudopotential parameters of double-zeta quality are presented for the first row atoms of the periodic table. The pseudopotential is especially well suited for incorporation into the floating-spherical-Gaussian-orbital (FSGO ) method, though it is not restricted to any particular method. Applications of the resulting pseudo-FSGO method to BeH₂, BH₃, CH₄, and C₂H₆ are presented.