Energetics and mechanism of the decomposition of trifluoromethanol

The thermal instability of alpha-fluoroalcohols is generally attributed to a unimolecular 1,2-elimination of HF, but the barrier to intramolecular HF elimination from CF3OH is predicted to be 45.1 +/- 2 kcal/mol. The thermochemical parameters of trifluoromethanol were calculated using coupled-cluster theory (CCSD(T)) extrapolated to the complete basis set limit. High barriers of 42.9, 43.1, and 45.0 kcal/mol were predicted for the unimolecular decompositions of CH2FOH, CHF2OH, and CF3OH, respectively. These barriers are lowered substantially if cyclic H-bonded dimers of CF3OH with complexation energies of approximately 5 kcal/mol are involved. A six-membered ring dimer has an energy barrier of 28.7 kcal/mol and an eight-membered dimer has an energy barrier of 32.9 kcal/mol. Complexes of CF3OH with HF lead to strong H-bonded dimers, trimers and tetramers with complexation energies of approximately 6, 11, and 16 kcal/mol, respectively. The dimer, CH3OH:HF, and the trimers, CF3OH:2HF and (CH3OH)2:HF, have decomposition energy barriers of 26.7, 20.3, and 22.8 kcal/mol, respectively. The tetramer (CH3OH:HF)2 gives rise to elimination of two HF molecules with a barrier of 32.5 kcal/mol. Either CF3OH or HF can act as catalysts for HF-elimination via an H-transfer relay. Because HF is one of the decomposition products, the decomposition reactions become autocatalytic. If the energies due to complexation for the CF3OH-HF adducts are not dissipated, the effective barriers to HF elimination are lowered from approximately 20 to approximately 9 kcal/mol, which reconciles the computational results with the experimentally observed stabilities.     

Trimethyl phosphate-acetylene interaction: a matrix-isolation infrared and ab initio study

Trimethyl phosphate (TMP) and acetylene were codeposited in nitrogen and argon matrices and adducts of these species were identified using infrared spectroscopy. Formation of the adducts was evidenced by shifts in the vibrational frequencies of the modes involving the TMP and acetylene submolecules. The structures of these adducts, energies and the vibrational frequencies were computed at the HF/6-31G** level. Both the experimental and computational studies indicated that two types of TMP-acetylene complexes were formed; one in which the hydrogen in acetylene was bonded to the phosphoryl oxygen and another in which the bonding was at the alkoxy oxygen of the phosphate. In addition to the primary hydrogen bonded interaction at the phosphoryl oxygen, this complex, also appeared to be stablilized by a secondary and weaker interaction involving a methyl hydrogen in TMP and the pi cloud in acetylene--a case of a H...pi interaction. The computed vibrational frequencies in the adducts agreed well with the observed frequencies for the modes involving the TMP submolecule, while the agreement was relatively poor for the modes involving the acetylene submolecule. The stabilization energies of these adducts, corrected for both zero-point energies and basis set superposition errors, were approximately 3 kcal/mol for the phosphoryl complex and, approximately 1 kcal/mol for the alkoxy complex.     

Unimolecular elimination of HF and HCl from chemically activated CF3CFClCH2Cl

The unimolecular reactions of CF3CFClCH2Cl molecules formed with 87 kcal mol(-1) of vibrational energy by recombination of CF3CFCl and CH2Cl radicals at room temperature have been characterized by the chemical activation technique. The 2,3-ClH and 2,3-FH elimination reactions, which have rate constants of (2.5 +/- 0.8) x 10(4) and (0.38 +/- 0.11) x 10(4) s(-1), respectively, are the major reactions. The 2,3-FCl interchange reaction was not observed. The trans (or E)-isomers of CF3CFCHCl and CF3CClCHCl are favored over the cis (or Z)-isomers. Density functional theory at the B3PW91/6-31G(d',p') level was used to evaluate thermochemistry and structures of the molecule and transition states. This information was used to calculate statistical rate constants. Matching the calculated to the experimental rate constants for the trans-isomers gave threshold energies of 62 and 63 kcal mol(-1) for HCl and HF elimination, respectively. The threshold energy for FCl interchange must be 3-4 kcal mol(-1) higher than for HF elimination. The results for CF3CFClCH2Cl are compared to those from CF3CFClCH3; the remarkable reduction in rate constants for HCl and HF elimination upon substitution of one Cl atom for one H atom is a consequence of both a lower E and higher threshold energies for CF3CFClCH2Cl.     

Mixed quantum-classical reaction path dynamics of C2H5F --> C2H4 + HF

A mixed quantum-classical method for calculating product energy partitioning based on a reaction path Hamiltonian is presented and applied to HF elimination from fluoroethane. The goal is to describe the effect of the potential energy release on the product energies using a simple model of quantized transverse vibrational modes coupled to a classical reaction path via the path curvature. Calculations of the minimum energy path were done at the B3LYP/6-311++G(2d,2p) and MP2/6-311++G** levels of theory, followed by energy-partitioning dynamics calculations. The results for the final HF vibrational state distribution were found to be in good qualitative agreement with both experimental studies and quasiclassical trajectory simulations.     

Theoretical study of pyridine and 4,4'-bipyridine adsorption on the lewis acid sites of alumina surfaces based on ab initio and density functional cluster calculations

The interactions of pyridine and 4,4'-bipyridine with the Lewis acid sites of alumina surfaces are investigated using ab initio and density functional calculations. Four cluster models of different sizes and shapes are chosen to represent the Lewis acid sites: three hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12) and one non-hydrogenated cluster Al(4)O(6). The Hartree-Fock (HF) and B3LYP approaches with two basis sets 6-31G and 6-31+G are used to calculate the geometries, the electronic structures, the vibrational frequencies, and the adsorption energies of the complexes formed upon interaction of pyridine or 4,4'-bipyridine ligands on the cluster surfaces. Electronic structures are determined by the electrostatic potential (ESP) analysis of charges. Adsorption energies are calculated with corrections made for zero-point energies (ZPE) and basis set superposition error (BSSE). The ESP analysis of atomic charges reveals that the charge-transfer effects are more important in Lewis complexes formed with Al(4)O(6) cluster than in those formed with hydrogenated clusters Al(OH)(3), Al(4)O(9)H(6), and Al(10)O(21)H(12). The significantly larger charge transferred from pyridine or 4,4'-bipyridine ligand to Al(4)O(6) cluster should increase the adsorption energy of these complexes. Consequently, at all levels of calculation, the adsorption energies of pyridine and 4,4'-bipyridine complexed to Al(4)O(6) cluster ( approximately 46 kcal/mol), which compare very well to experiment, are strongly larger than those obtained for both pyridine and 4,4'-bipyridine ligands complexed to Al(OH)(3) (32 kcal/mol), Al(4)O(9)H(6) (24 kcal/mol) and Al(10)O(21)H(12) (25 kcal/mol) clusters. The corrected adsorption energy is found to be insensitive to basis set and electron correlation effects. It essentially depends on the ionic character of the cluster model rather than on its size. For 4,4'-bipyridine complexes, similar results to those obtained for pyridine are found, and the geometry and the amount of charge of the unbound pyridyl ring are unchanged upon complexation. The calculated vibrational frequencies and frequency shifts are little sensitive to the size and shape of the cluster model. The two ring stretching modes 8a and 19b of pyridine and 4,4'-bipyridine observed in the 1400-1600 cm(-1) region are the most affected modes upon adsorption, in good agreement with the available infrared and Raman data.     

CCSD(T) study of CD3-O-CD3 and CH3-O-CD3 far-infrared spectra

From a vibrationally corrected 3D potential energy surface determined with highly correlated ab initio calculations (CCSD(T)), the lowest vibrational energies of two dimethyl-ether isotopologues, (12)CH(3)-(16)O-(12)CD(3) (DME-d(3)) and (12)CD(3)-(16)O-(12)CD(3) (DME-d(6)), are computed variationally. The levels that can be populated at very low temperatures correspond to the COC-bending and the two methyl torsional modes. Molecular symmetry groups are used for the classification of levels and torsional splittings. DME-d(6) belongs to the G(36) group, as the most abundant isotopologue (12)CH(3)-(16)O-(12)CH(3) (DME-h(6)), while DME-d(3) is a G(18) species. Previous assignments of experimental Raman and far-infrared spectra are discussed from an effective Hamiltonian obtained after refining the ab initio parameters. Because a good agreement between calculated and experimental transition frequencies is reached, new assignments are proposed for various combination bands corresponding to the two deuterated isotopologues and for the 020 → 030 transition of DME-d(6). Vibrationally corrected potential energy barriers, structural parameters, and anharmonic spectroscopic parameters are provided. For the 3N - 9 neglected vibrational modes, harmonic and anharmonic fundamental frequencies are obtained using second-order perturbation theory by means of CCSD and MP2 force fields. Fermi resonances between the COC-bending and the torsional modes modify DME-d(3) intensities and the band positions of the torsional overtones.     

Ab-initio molecular geometry and normal coordinate analysis of tetrahydrothiophene molecule.

The molecular geometry of tetrahydrothiophene (THT) was quantum mechanically calculated using the split valence 6-31G** basis set. Electron correlation energy has been computed employing MP2 method. The molecule showed a twist form puckered structure with a twist torsion angle of 13 degrees and has a total energy of -347,877.514 kcal/mol of which a 436.715 kcal/mol electron correlation energy. The envelope form of the molecule showed an inter-plane angle of 22 degrees and has a total energy of -347,874.430 kcal/mol involving -436.558 kcal/mol electron correlation energy. The normal coordinates of the molecule were theoretically analyzed and the fundamental vibrational frequencies were calculated. The IR and laser Raman spectra of THT molecule was measured. All the observed vibrational bands including combination bands and overtones were assigned to normal modes with the aid of the potential energy distribution values obtained from normal coordinate calculations. The molecular force field was determined by refining the initial set of force constants using the least square fit method instead of using the less accurate scaling factor methods. The determined molecular force field has produced simulated frequencies which best match the observed values. The lowest-energy modes of vibration were two molecular out-of-plane deformations, observed at 114 and 166 cm(-1). The barrier of ring twisting estimated from the observed ring out-of-plane vibrational mode at 114 cm(-1) was estimated.     

Unimolecular rate constants, kinetic isotope effects and threshold energies for FH and FD elimination from CF3CHFCH3 and CF3CHFCD3

The combination of CF(3)CHF and CH(3) or CD(3) radicals was used to prepare vibrationally excited CF(3)CHFCH(3) or CF(3)CHFCD(3) molecules with 97 kcal mol(-1) of internal energy. The experimental unimolecular rate constants were 3.7 x 10(6) s(-1) for 2,3-FH elimination from CF(3)CHFCH(3) and 1.3 x 10(6) s(-1) for 2,3-DF elimination from CF(3)CHFCD(3). Unimolecular rate constants for 1,2-FH elimination reaction were approximately 230 and 98 times smaller for CF(3)CHFCH(3) and CF(3)CHFCD(3), respectively, than the corresponding rate constants for 2,3-FH elimination. Density functional theory (DFT) was used to calculate the structures and vibrational frequencies of the molecules and transition states; this information was subsequently employed for calculations of RRKM rate constants. Comparison of the experimental and calculated rate constants gave a threshold energy of 73 +/- 2 kcal mol(-1) for the 1,2-FH elimination process and 60.5 +/- 1.5 kcal mol(-1) for the 2,3-FH elimination reaction from CF(3)CHFCH(3). The calculated kinetic-isotope effects agree with the experimental results. The experimentally derived threshold energies for 1,2-FH and 2,3-FH elimination reactions from several fluoropropanes and fluorochloropropanes are summarized and compared to those from DFT calculations.     

The ground and two lowest-lying singlet excited electronic states of copper hydroxide (CuOH)

Various ab initio methods, including self-consistent field (SCF), configuration interaction, coupled cluster (CC), and complete-active-space SCF (CASSCF), have been employed to study the electronic structure of copper hydroxide (CuOH). Geometries, total energies, dipole moments, harmonic vibrational frequencies, and zero-point vibrational energies are reported for the linear 1Sigma+ and 1Pi stationary points, and for the bent ground-state X 1A', and excited-states 2 1A' and 1 1A". Six different basis sets have been used in the study, Wachters/DZP being the smallest and QZVPP being the largest. The ground- and excited-state bending modes present imaginary frequencies for the linear stationary points, indicating that bent structures are more favorable. The effects of relativity for CuOH are important and have been considered using the Douglas-Kroll approach with cc-pVTZ/cc-pVTZ_DK and cc-pVQZ/cc-pVQZ_DK basis sets. The bent ground and two lowest-lying singlet excited states of the CuOH molecule are indeed energetically more stable than the corresponding linear structures. The optimized geometrical parameters for the X 1A' and 1 1A" states agree fairly well with available experimental values. However, the 2 1A' structure and rotational constants are in poor agreement with experiment, and we suggest that the latter are in error. The predicted adiabatic excitation energies are also inconsistent with the experimental values of 45.5 kcal mol(-1) for the 2 1A' state and 52.6 kcal mol(-1) for the 1 1A" state. The theoretical CC and CASSCF methods show lower adiabatic excitation energies for the 1 1A" state (53.1 kcal mol(-1)) than those for the corresponding 2 1A' state (57.6 kcal mol(-1)), suggesting that the 1 1A" state might be the first singlet excited state while the 2 1A' state might be the second singlet excited state.     

Rate constants and kinetic isotope effects for unimolecular 1,2-HX or DX (X = F or Cl) elimination from chemically activated CF3CFClCH3-d0, -d1, -d2, and -d3

Chemically activated CF(3)CFClCH(3), CF(3)CFClCD(3), CF(3)CFClCH(2)D, and CF(3)CFClCHD(2) molecules with 94 kcal mol(-1) of internal energy were formed by the combination of CF(3)CFCl radicals with CH(3), CD(3), CH(2)D, and CHD(2) radicals, which were generated from UV photolysis of CF(3)CFClI and CH(3)I, CD(3)I, CH(2)DI, or CHD(2)I. The total (HF + HCl) elimination rate constants for CF(3)CFClCH(3) and CF(3)CFClCD(3) were 5.3 x 10(6) and 1.7 x 10(6) s(-1) with product branching ratios of 8.7 +/- 0.6 in favor of HCl (or DCl). The intermolecular kinetic isotope effects were 3.22 and 3.18 for the HCl and HF channels, respectively. The product branching ratios were 10.3 +/- 1.9 and 11.8 +/- 1.8 (10.8 +/- 3.8 and 11.6 +/- 1.7) for HCl/HF and DCl/DF, respectively, from CF(3)CFClCH(2)D (CF(3)CFClCHD(2)). The intramolecular kinetic-isotope effects (without correction for reaction path degeneracy) for HCl/DCl and HF/DF elimination from CF(3)CFClCH(2)D (CF(3)CFClCHD(2)) were 2.78 +/- 0.16 and 2.98 +/- 0.12 (0.82 +/- 0.04 and 0.91 +/- 0.03), respectively. Density function theory at the B3PW91/6-311+G(2d,p) and B3PW91/6-31G(d',p') levels was investigated, and the latter was chosen to calculate frequencies and moments of inertia for the molecules and transition states. Rate constants, branching ratios and kinetic-isotope effects then were calculated using RRKM theory with torsional motions treated as hindered internal rotations. Threshold energies for HF and HCl elimination from CF(3)CFClCH(3) were assigned as 61.3 +/- 1.5 and 58.5 +/- 1.5 kcal mol(-1), respectively. The threshold energy for Cl-F interchange was estimated as 67 kcal mol(-1). The difference between the transition states for HCl and HF elimination is discussed.     

Molecular mechanics (MM3) calculations on conjugated hydrocarbons

The MM3 molecular mechanics program has been extended to conjugated systems. A VESCF method is applied to the pi-system to calculate bond orders, from which various stretching and torsional parameters are obtained. The procedure gives somewhat better results than the analogous MM2 calculations. It has been applied to a study of 81 compounds of aromatic and other conjugated hydrocarbons, as well as 45 alkenes and unconjugated polyenes. The structures calculated are generally in good agreement with experiment, and the heats of formation of these compounds can be calculated with a rms value of 0.62 kcal/mol, which may be compared with the average experimental error of 0.61 kcal/mol. In addition, vibrational frequencies for five representative conjugated model structures are calculated, with an rms value of 46 cm^?1, and from these, other properties such as entropy can be calculated.     

Quantum and classical studies of the O(3P) + H2(v = 0-3,j = 0) --> OH + H reaction using benchmark potential surfaces

We present results of time-dependent quantum mechanics (TDQM) and quasiclassical trajectory (QCT) studies of the excitation function for O(3P) + H2(v = 0-3,j = 0) --> OH + H from threshold to 30 kcal/mol collision energy using benchmark potential energy surfaces [Rogers et al., J. Phys. Chem. A 104, 2308 (2000)]. For H2(v = 0) there is excellent agreement between quantum and classical results. The TDQM results show that the reactive threshold drops from 10 kcal/mol for v = 0 to 6 for v = 1, 5 for v = 2 and 4 for v = 3, suggesting a much slower increase in rate constant with vibrational excitation above v = 1 than below. For H2(v > 0), the classical results are larger than the quantum results by a factor approximately 2 near threshold, but the agreement monotonically improves until they are within approximately 10% near 30 kcal/mol collision energy. We believe these differences arise from stronger vibrational adiabaticity in the quantum dynamics, an effect examined before for this system at lower energies. We have also computed QCT OH(v',j') state-resolved cross sections and angular distributions. The QCT state-resolved OH(v') cross sections peak at the same vibrational quantum number as the H2 reagent. The OH rotational distributions are also quite hot and tend to cluster around high rotational quantum numbers. However, the dynamics seem to dictate a cutoff in the energy going into OH rotation indicating an angular momentum constraint. The state-resolved OH distributions were fit to probability functions based on conventional information theory extended to include an energy gap law for product vibrations.     

Quantum Monte Carlo calculations of the dissociation energy of the water dimer

We report diffusion quantum Monte Carlo (DMC) calculations of the equilibrium dissociation energy D(e) of the water dimer. The dissociation energy measured experimentally, D(0), can be estimated from D(e) by adding a correction for vibrational effects. Using the measured dissociation energy and the modern value of the vibrational energy Mas et al., [J. Chem. Phys. 113, 6687 (2000)] leads to D(e)=5.00+/-0.7 kcal mol(-1), although the result Curtiss et al., [J. Chem. Phys. 71, 2703 (1979)] D(e)=5.44+/-0.7 kcal mol(-1), which uses an earlier estimate of the vibrational energy, has been widely quoted. High-level coupled cluster calculations Klopper et al., [Phys. Chem. Chem. Phys. 2, 2227 (2000)] have yielded D(e)=5.02+/-0.05 kcal mol(-1). In an attempt to shed new light on this old problem, we have performed all-electron DMC calculations on the water monomer and dimer using Slater-Jastrow wave functions with both Hartree-Fock approximation (HF) and B3LYP density functional theory single-particle orbitals. We obtain equilibrium dissociation energies for the dimer of 5.02+/-0.18 kcal mol(-1) (HF orbitals) and 5.21+/-0.18 kcal mol(-1) (B3LYP orbitals), in good agreement with the coupled cluster results.     

A theoretical study of the rotational isomers of nitrosomethanol by semiempirical (AM1) and ab initio methods

The conformational potential energy surface as a function of the two internal torsion angles in C-nitrosomethanol has been obtained using the semiempirical AM1 method. Optimized geometries are reported for the local minima on this surface and also for the corresponding points on the HF/6-31G, 6-31G*, and 6-31G** surfaces. All methods predict cis and trans minima which occur in degenerate pairs, each pair being connected by a transition state of C_s symmetry. The AM1 structures are found to compare well with the corresponding ab initio structures. Ab initio HF/6-31G and HF/6-31G* harmonic vibrational frequencies are reported for the cis and trans forms of nitrosomethanol. When scaled appropriately the calculated frequencies are found to compare well with experimental frequencies. The ab initio calculations predict the energy barrier for cis → trans isomerization to be between 5.8 and 6.5 kcal/mol with the trans → cis isomerization barrier lying between 2.3 and 6.5 kcal/mol. The corresponding AM1 energy barriers are around 1 kcal/mol lower in energy. The ab initio calculations predict the barrier to conversion between the two cis rotamers to be very small with the AM1 value being around 1 kcal/mol. Both AM1 and ab initio calculations predict interconversion between trans rotamers to require between 1.2 and 1.4 kcal/mol.     

Vibrational investigation on FT-IR and FT-Raman spectra, IR intensity, Raman activity, peak resemblance, ideal estimation, standard deviation of computed frequencies analyses and electronic structure on 3-methyl-1,2-butadiene using HF and DFT (LSDA/B3LYP/B3PW91) calculations

FT-IR and FT-Raman (4000–100 cm⁻¹) spectral measurements of 3-methyl-1,2-butadiene (3M12B) have been attempted in the present work. Ab-initio HF and DFT (LSDA/B3LYP/B3PW91) calculations have been performed giving energies, optimized structures, harmonic vibrational frequencies, IR intensities and Raman activities. Complete vibrational assignments on the observed spectra are made with vibrational frequencies obtained by HF and DFT (LSDA/B3LYP/B3PW91) at 6-31G(d,p) and 6-311G(d,p) basis sets. The results of the calculations have been used to simulate IR and Raman spectra for the molecule that showed good agreement with the observed spectra. The potential energy distribution (PED) corresponding to each of the observed frequencies are calculated which confirms the reliability and precision of the assignment and analysis of the vibrational fundamentals modes. The oscillation of vibrational frequencies of butadiene due to the couple of methyl group is also discussed. A study on the electronic properties such as HOMO and LUMO energies, were performed by time-dependent DFT (TD-DFT) approach. The calculated HOMO and LUMO energies show that charge transfer occurs within the molecule. The thermodynamic properties of the title compound at different temperatures reveal the correlations between standard heat capacities (C) standard entropies (S), and standard enthalpy changes (H).     

Unimolecular rate constants for HX or DX elimination (X = F, Cl) from chemically activated CF3CH2CH2Cl, C2H5CH2Cl, and C2D5CH2Cl: threshold energies for HF and HCl elimination

Vibrationally activated CF(3)CH(2)CH(2)Cl molecules were prepared with 94 kcal mol(-1) of vibrational energy by the combination of CF(3)CH(2) and CH(2)Cl radicals and with 101 kcal mol(-1) of energy by the combination of CF(3) and CH(2)CH(2)Cl radicals at room temperature. The unimolecular rate constants for elimination of HCl from CF(3)CH(2)CH(2)Cl were 1.2 x 10(7) and 0.24 x 10(7) s(-1) with 101 and 94 kcal mol(-1), respectively. The product branching ratio, k(HCl)/k(HF), was 80 +/- 25. Activated CH(3)CH(2)CH(2)Cl and CD(3)CD(2)CH(2)Cl molecules with 90 kcal mol(-1) of energy were prepared by recombination of C(2)H(5) (or C(2)D(5)) radicals with CH(2)Cl radicals. The unimolecular rate constant for HCl elimination was 8.7 x 10(7) s(-1), and the kinetic isotope effect was 4.0. Unified transition-state models obtained from density-functional theory calculations, with treatment of torsions as hindered internal rotors for the molecules and the transition states, were employed in the calculation of the RRKM rate constants for CF(3)CH(2)CH(2)Cl and CH(3)CH(2)CH(2)Cl. Fitting the calculated rate constants from RRKM theory to the experimental values provided threshold energies, E(0), of 58 and 71 kcal mol(-1) for the elimination of HCl or HF, respectively, from CF(3)CH(2)CH(2)Cl and 54 kcal mol(-1) for HCl elimination from CH(3)CH(2)CH(2)Cl. Using the hindered-rotor model, threshold energies for HF elimination also were reassigned from previously published chemical activation data for CF(3)CH(2)CH(3,) CF(3)CH(2)CF(3), CH(3)CH(2)CH(2)F, CH(3)CHFCH(3), and CH(3)CF(2)CH(3). In an appendix, the method used to assign threshold energies was tested and verified using the combined thermal and chemical activation data for C(2)H(5)Cl, C(2)H(5)F, and CH(3)CF(3).     

Si~2Cl~6分子的内旋转能垒和振动基频的理论研究

用从头算方法HF/6-31G^*^*和密度函方法B3LYP/6-31G^*^*,对Si~2Cl~6分子的平衡几何构型进行优化,优化的结果与实验结果吻合得较好.并用上述两种不同的方法计算Si~2Cl~6分子的内旋转能垒,结果分别为8.786和6.694kJ/mol,其中DFT方法的计算结果与实验结果4.18kJ/mol吻合得较好.对Si~2Cl~6分子的振动基频进行计算.用HF/6-31G^*^*SQM力场所计算的频率理论值与实验值的平均误差为7.3cm^-^1,用B3LYP/6-31G^*^*未标度的力场所计算的频率理论值与实验值的平均误差为6.0cm^-^1.该密度泛函方法(B3LYP/~6-31G^*^*)的理论计算值比用HF/6-31G^*^*标度后的SQM力场计算的频率与实验值(除Si--Si键扭转振动基频之外的11条振动基频)吻合得更好.并给出了Si--Si键扭转振动基频的预测值。     

An ab initio based full-dimensional global potential energy surface for FH(2)O(X(2)A(')) and dynamics for the F + H(2)O → HF + HO reaction

A global potential energy surface (PES) for the ground electronic state of FH(2)O is constructed based on more than 30 000 ab initio points at the multi-reference configuration interaction level. The PES features a pre-reaction van der Waals well and two post-reaction hydrogen-bonded complexes, as well as a "reactant-like" transition state with a classical barrier of 3.8 kcal∕mol. The adiabatic F + H(2)O → HF + OH reaction dynamics on this PES was investigated using a standard quasi-classical trajectory method. In agreement with experiment, the HF product contains significant vibrational excitation with limited rotational excitation, while the OH product is internally cold, reflecting its spectator role in the reaction. The products are primarily scattered in the backward direction, consistent with a direct abstraction mechanism.     

The diazocarbene (CNN) molecule: characterization of the X 3Sigma- and A 3Pi electronic states

The ground (X (3)Sigma(-)) and first excited triplet (A (3)Pi) electronic states of diazocarbene (CNN) have been investigated systematically starting from the self-consistent-field theory and proceeding to the coupled cluster with single, double, and full triple excitations (CCSDT) method with a wide range of basis sets. While the linear X (3)Sigma(-) ground state of CNN has a real degenerate bending vibrational frequency, the A (3)Pi state of CNN is subject to the Renner-Teller effect and presents two distinct real vibrational frequencies along the bending coordinate. The bending vibrational frequencies of the A (3)Pi state were evaluated via the equation-of-motion coupled cluster (EOM-CC) techniques. The significant sensitivity to level of theory in predicting the ground-state geometry, harmonic vibrational frequencies, and associated infrared intensities has been attributed to the fact that the reference wave function is strongly perturbed by the excitations of 1pi-->3pi followed by a spin flip. At the highest level of theory with the largest basis set, correlation-consistent polarized valence quadruple zeta (cc-pVQZ) CCSDT, the classical X-A splitting (T(e) value) was predicted to be 68.5 kcal/mol (2.97 eV, 24 000 cm(-1)) and the quantum mechanical splitting (T(0) value) to be 69.7 kcal/mol (3.02 eV, 24 400 cm(-1)), which are in excellent agreement with the experimental T(0) values, 67.5-68.2 kcal/mol (2.93-2.96 eV, 23 600-23 900 cm(-1)). With the EOM-CCSD method the Renner parameter (epsilon) and averaged bending vibrational frequency (omega(2)) for the A (3)Pi state were evaluated to be epsilon=-0.118 and omega(2)=615 cm(-1), respectively. They are in fair agreement with the experimental values of epsilon=-0.07 and nu(2)=525 cm(-1).     

A theoretical investigation of atmospheric sulfur chemistry. 1. The HSO/HOS energy separation and the heat of formation of HSO, HOS, and HS2

The energy separation between the ground-state structures of HSO and HOS has been determined by using two independent ab initio methods. In the first method, the optimized geometry of all species was obtained at the HF/6-31G(d) level, as were harmonic vibrational frequencies for zero-point energy corrections. The energies were calculated by using fourth-order Moller-Plesset perturbation theory and a 6-31G(d,p) basis set. After corrections for extrapolation of the Moller-Plesset series to infinite order and extension of the basis set to include diffuse sp-, extra d-, and f-type Gaussian functions, the predicted energy separation, including zero-point vibrational effects, is 2.5 kcal/mol. HOS is the more stable isomer. The second method uses a double-zeta basis augmented with an extra set of p functions and two sets of d functions on the sulfur and oxygen atoms and a double-zeta + p basis on hydrogen. With this basis, equilibrium structures of HSO and HOS were obtained from MCSCF calculations; the energy separation between these structures was corrected by using large scale configuration interaction. In good agreement with the first method, HOS is the more stable isomer by 3.1 kcal/mol. Through calculation of the energy change in the reaction HO2 + XY --> O2 + HXY, the first method predicts the heats of formation of HXY = HSO, HOS, and HS2 to be -0.4, -2.9, and 26.7 kcal/mol, respectively.