Molecular structure of thiourea

The molecular structure of thiourea has been investigated under C(s), C(2), and C(2v) symmetry constraints. At the coupled-cluster level in conjunction with a triple-ζ basis set, only the C(2) conformer has been found to be a real minimum on the potential energy surface. Its equilibrium structure has therefore been accurately evaluated using both theoretical and experimental data. With respect to the former, high-level quantum-chemical calculations at the coupled-cluster level in conjunction with correlation-consistent basis sets ranging in size from triple- to quintuple-zeta have been carried out. Extrapolation to the complete basis-set limit as well as core-correlation effects and inclusion of full treatment of triple excitations in the cluster operator have been considered. On the basis of the vibrational ground-state rotational constants available for five isotopic species and the corresponding computed vibrational corrections, the semiexperimental equilibrium geometry of thiourea has also been determined for the first time.     

Ab initio prediction of the potential energy surface and vibrational-rotational energy levels of dialuminum monoxide, Al2O

The equilibrium structure and potential energy surface of dialuminum monoxide, Al(2)O, have been determined from large-scale ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of triple- through quintuple-zeta quality. The effects of core-electron correlation on the calculated molecular parameters were investigated. The vibrational-rotational energy levels of the Al(2) (16)O and Al(2) (18)O isotopic species were calculated by a variational approach. The predicted energy levels are in remarkably good agreement with the available experimental spectroscopic data (from laser-induced fluorescence), demonstrating that the Al(2)O molecule is linear at equilibrium in its ground electronic state. The reported theoretical data settle controversies between the experimental studies about the equilibrium structure and assignment of vibrational fundamentals of the Al(2)O molecule.     

Ab initio prediction of the potential energy surface and vibrational-rotational energy levels of calcium dihydride, CaH2

The equilibrium structure and potential energy surface of calcium dihydride, CaH(2), have been determined from large-scale ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of quadruple- and quintuple-zeta quality. The CaH(2) molecule was found to be quasilinear. The HCaH bending potential function was predicted to be extraordinarily flat near the minimum, located at the HCaH angle of 164 degrees. The barrier to linearity was calculated to be just 6 cm(-1). The vibrational-rotational energy levels of various isotopomers were predicted using the variational method. The calculated vibrational fundamental frequencies are in good agreement with the results of matrix-isolation studies, and the other predicted spectroscopic constants can assist in the future detection of calcium dihydride in the gas phase.     

The structures,binding energies and vibrational frequencies of Ca3 and Ca4 — An application of the CCSD(T) method

Summary The Ca₃ and Ca₄ metallic clusters have been investigated using state-of-the-artab initio quantum mechanical methods. Large atomic natural orbital basis sets have been used in conjunction with the singles and doubles coupled-cluster (CCSD) method, a coupled-cluster method that includes a perturbational estimate of connected triple excitations, denoted CCSD(T), and the multireference configuration interaction (MRCI) method. The equilibrium geometries, binding energies and harmonic vibrational frequencies have been determined with each of the methods so that the accuracy of the coupled-cluster methods may be assessed. Since the CCSD(T) method reproduces the MRCI results very well, cubic and quartic force fields of Ca₃ and Ca₄ have been determined using this approach and used to evaluate the fundamental vibrational frequencies. The infrared intensities of both thee mode of Ca₃ and thet ₂ mode of Ca₄ are found to be small. The results obtained in this study are compared and contrasted with those from our earlier studies on small Be and Mg clusters.Dedicated to Prof. Klaus Ruedenberg on the occasion of his 70th birthday     

Ab initio studies of CIO, reactions: prediction of the rate constants of CIO + NO for the forward and reverse processes.

The mechanisms for ClO+NO and its reverse reactions were investigated by means of ab initio molecular orbital and statistical theory calculations. The species involved were optimized at the B3LYP/6-311 +G(3df) level, and their energies were refined at the CCSD(T)/6-311+ G(3df)//B3LYP/6-311 + G(3df) level. Five isomers and the transition states among them were located. The relative stability of these isomers is ClNO2 > cis-ClONO > trans-ClONO > cis-OClNO>trans-OClNO. The heats of formation of the three most-stable isomers were predicted using isodesmic reactions by different methods. The predicted bimolecular reaction rate constant shows that, below 100 atm, the formation of Cl+NO2 is dominant and pressure-independent. The total rate constant can be expressed as: k(ClO+NO)= 1.43 x 10(-9)T(-083)exp(92/ T) cm3 molecule(-1)s(-1) in the temperature range of 200-1000 K, in close agreement with experimental data. For the reverse reaction, Cl+NO2-->ClNO2 and ClONO (cis and trans isomers), the sum of the predicted rate constants for the formation of the three isomers and their relative yields also reproduce the experimental data well. The predicted total third-order rate constants in the temperature range of 200-1000 K can be represented by: k0(He) = 4.89 x 10(-6)T(-5.85) exp(-796/T) cm6 molecule(-1)s(-1) and k0(N2) =5.72 x 10(-15)T(-5.80) exp(-814/T) cm6 molecule(-1)s(-1). The predicted high- and low-pressure limit decomposition rates of CINO2 in Ar in the temperature range 400-1500 K can be expressed, respectively, by: k-(ClNO2) = 7.25 x 10(19)T(-1.89) exp(-16875/T) s(-1) and kd(ClNO2) = 2.51 x 10(38)T(-6.8) exp(-18409/T) cm3 molecule(-1) s(-1). The value of k0(ClNO2) is also in reasonable agreement with available experimental data.     

Accurate ab initio structure, dissociation energy, and vibrational spectroscopy of the F(-)-CH4 anion complex

Accurate equilibrium structure, dissociation energy, global potential energy surface (PES), dipole moment surface (DMS), and the infrared vibrational spectrum in the 0-3000 cm(-1) range of the F(-)-CH4 anion complex have been obtained. The equilibrium electronic structure calculations employed second-order M?ller-Plesset perturbation theory (MP2) and coupled-cluster (CC) method up to single, double, triple, and perturbative quadruple excitations using the aug-cc-p(C)VXZ [X = 2(D), 3(T), 4(Q), and 5] correlation-consistent basis sets. The best equilibrium geometry has been obtained at the all-electron CCSD(T)/aug-cc-pCVQZ level of theory. The dissociation energy has been determined based on basis set extrapolation techniques within the focal-point analysis (FPA) approach considering (a) electron correlation beyond the all-electron CCSD(T) level, (b) relativistic effects, (c) diagonal Born-Oppenheimer corrections (DBOC), and (d) variationally computed zero-point vibrational energies. The final D(e) and D0 values are 2398 +/- 12 and 2280 +/- 20 cm(-1), respectively. The global PES and DMS have been computed at the frozen-core CCSD(T)/aug-cc-pVTZ and MP2/aug-cc-pVTZ levels of theory, respectively. Variational vibrational calculations have been performed for CH4 and F(-)-CH4 employing the vibrational configuration interaction (VCI) method as implemented in Multimode.     

Ab initio chemical kinetics for the ClOO + NO reaction: effects of temperature and pressure on product branching formation

The kinetics and mechanism for the reaction of ClOO with NO have been investigated by ab initio molecular orbital theory calculations based on the CCSD(T)/6-311+G(3df)//PW91PW91∕6-311+G(3df) method, employed to evaluate the energetics for the construction of potential energy surfaces and prediction of reaction rate constants. The results show that the reaction can produce two key low energy products ClNO + (3)O(2) via the direct triplet abstraction path and ClO + NO(2) via the association and decomposition mechanism through long-lived singlet pc-ClOONO and ClONO(2) intermediates. The yield of ClNO + O(2) ((1)△) from any of the singlet intermediates was found to be negligible because of their high barriers and tight transition states. As both key reactions initially occur barrierlessly, their rate constants were evaluated with a canonical variational approach in our transition state theory and Rice-Ramspergen-Kassel-Marcus/master equation calculations. The rate constants for ClNO + (3)O(2) and ClO + NO(2) production from ClOO + NO can be given by 2.66 × 10(-16) T(1.91) exp(341/T) (200-700 K) and 1.48 × 10(-24) T(3.99) exp(1711/T) (200-600 K), respectively, independent of pressure below atmospheric pressure. The predicted total rate constant and the yields of ClNO and NO(2) in the temperature range of 200-700 K at 10-760 Torr pressure are in close agreement with available experimental results.     

A theoretical study on CH<Subscript>2</Subscript>N<Subscript>2</Subscript> isomers: structure and energetics

Five CH₂N₂ isomers, namely cyanamide, carbodiimide, diazomethane, isocyanamide and nitrilimine, have been investigated at a high level of accuracy. The singles and doubles coupled-cluster method including a perturbational correction for connected triple excitations, CCSD(T), in conjunction with correlation-consistent basis sets ranging in size from triple to quintuple zeta have been employed. Extrapolation to the complete basis set limit has been used with treatments of core-valence correlation effects in order to accurately predict structures, relative energies as well as N–H and C–H bond dissociation energies. The latter required to also investigate the HNNC radical with the same methodology used for CH₂N₂ isomers, while HCNN and HNCN data are available in the literature by the same authors (Puzzarini and Gambi in J Chem Phys 122:064316, 2005). For all the species studied, harmonic vibrational frequencies have also been evaluated at the CCSD(T) level in order to obtain zero-point corrections to total energies.     

Theoretical studies of UO2(H2O)n(2+), NpO2(H2O)n(+), and PuO2(H2O)n(2+) complexes (n=4-6) in aqueous solution and gas phase

Extensive ab initio calculations both in gas phase and solution have been carried out to study the equilibrium structure, vibrational frequencies, and bonding characteristics of various actinyl (UO2(2+), NpO2(+), and PuO2(2+)) and their hydrated forms, AnO2(H2O)n(z+) (n=4, 5, and 6). Bulk solvent effects were studied using a continuum method. The geometries were fully optimized at the coupled-cluster singles + doubles (CCSD), density-functional theory (DFT), and M?ller-Plesset (MP2) level of theories. In addition vibrational frequencies have been obtained at the CCSD as well as MP2/DFT levels. The results show that both the short-range and long-range solvent effects are important. The combined discrete-continuum model, in which the ionic solute and the solvent molecules in the first and second solvation shells are treated quantum mechanically while the solvent is simulated by a continuum model, can predict accurately the bonding characteristics. Moreover, our values of solvation free energies suggest that five- and six-coordinations are equally preferred for UO2(2+), and five-coordinated species are preferred for NpO2(+) and PuO2(2+). On the basis of combined quantum-chemical and continuum treatments of the hydrated complexes, we are able to determine the optimal cavity radii for the solvation models. The coupled-cluster computations with large basis sets were employed for the vibrational spectra and equilibrium geometries both of which compare quite favorably with experiment. Our most accurate computations reveal that both five- and six-coordination complexes are important for these species.     

The vibration-rotation emission spectrum of hot BeF2

The high-resolution infrared emission spectrum of BeF2 vapor at 1000 degrees C was rotationally analyzed with the assistance of large-scale ab initio calculations using the coupled-cluster method including single and double excitations and perturbative inclusion of triple excitations, in conjunction with correlation-consistent basis sets up to quintuple-zeta quality. The nu3 fundamental band, the nu1+nu2, nu1+nu3, and 2nu2+nu3 combination bands, and 18 hot bands were assigned. The symmetric stretching (nu1), bending (nu2), and antisymmetric stretching (nu3) mode frequencies were determined to be 769.0943(2), 342.6145(3), and 1555.0480(1) cm-1, respectively, from the band origins of the nu3, nu1+nu3, and nu1+nu2 bands. The observed vibrational term values and B rotational constants were fitted simultaneously to an effective Hamiltonian model with Fermi resonance taken into account, and deperturbed equilibrium vibrational and rotational constants were obtained for BeF2. The equilibrium rotational constant (Be) was determined to be 0.235 354(41) cm-1, and the associated equilibrium bond distance (re) is 1.3730(1) A. The results of our ab initio calculations are in remarkably good agreement with those of our experiment, and the calculated value was 1.374 A for the equilibrium bond distance (re). As in the isoelectronic CO2 molecule, the Fermi resonance in BeF2 is very strong, and the interaction constant k122 was found to be 90.20(4) cm-1.     

The triplet state of cytosine and its derivatives: electron impact and quantum chemical study

The excitation of the lowest electronic states and vibrational excitation of cytosine (C) have been studied using electron energy loss spectroscopy (EELS, 0-100 eV) with angular analysis. The singlet states have been found to be in good agreement with UV-VIS absorption results on sublimed films, slightly blueshifted by about 0.1 eV. The EEL spectra recorded at residual energy below 2 eV show clear shoulders at energy losses of 3.50 and 4.25 eV (+/-0.1 eV). They are assigned to the lowest triplet electronic states of cytosine. Energies and molecular structures of the lowest-lying triplet state of C and its methylated and halogenated 5-X-C, 6-X-C, and 5-X, 6-X-C substituted derivatives (X=CH3, F, Cl, and Br) have been studied using quantum chemical calculations with both molecular orbital and density functional methods, in conjunction with the 6-311++G(d,p), 6-311++G(3df,2p), and aug-cc-pVTZ basis sets. The triplet-singlet energy gap obtained using coupled-cluster theory [CCSD(T)] and density functional theory (DFT) methods agrees well with those derived from EELS study. The first C's vertical triplet state is located at 3.6 eV, in good agreement with experiment. The weak band observed at 4.25 eV is tentatively assigned to the second C's vertical triplet excitation. For the substituted cytosines considered, the vertical triplet state is consistently centered at 3.0-3.2 eV above the corresponding singlet ground state but about 1.0 eV below the first excited singlet state. Geometrical relaxation involving out-of-plane distortions of hydrogen atoms leads to a stabilization of 0.6-1.0 eV in favor of the equilibrium triplet. The lowest-lying adiabatic triplet states are located at 2.3-3.0 eV. Halogen substitution at both C(5) and C(6) positions tends to reduce the triplet-singlet separations whereas methylation tends to enlarge it. The vibrational modes of triplet cytosine and the ionization energies of substituted derivatives were also evaluated.     

Ab initio prediction of the potential energy surface and vibration-rotation energy levels of CaCl2

The equilibrium structure and potential energy surface of calcium dichloride (CaCl2) have been determined from accurate ab initio calculations using the coupled-cluster method, CCSD(T), in conjunction with basis sets of quadruple- and quintuple-zeta quality. The CaCl2 molecule was found to be linear at equilibrium. The vibration-rotation energy levels of various CaCl2 isotopomers were predicted by the variational method. The calculated spectroscopic constants could be used to guide future high-resolution spectroscopic experiments on calcium dichloride.     

Ab initio calculation of the vibrational and electronic spectra of trans- and cis-azobenzene

We report accurate geometries and harmonic force fields for trans- and cis-azobenzene determined by second-order M?ller-Plesset perturbation theory. For the trans isomer, the planar structure with C(2h) symmetry, found in a recent gas electron diffraction experiment, is verified. The calculated vibrational spectra are compared with experimental data and density functional calculations. Important vibrational frequencies are localized and discussed. For both isomers, we report UV spectra calculated using the second-order approximate coupled-cluster singles-and-doubles model CC2 with accurate basis sets. Vertical excitation energies and oscillator strengths have been determined for the lowest singlet n(pi)* and (pi)(pi)* transitions. The results are compared with the available experimental data and second-order polarization propagator (SOPPA) and density functional (DFT) calculations. For both isomers, the CC2 results for the excitation energies into the S(1) and S(2) states agree within 0.1 eV with experimental gas-phase measurements.     

Experimental and theoretical UV characterizations of acetylacetone and its isomers

Cryogenic matrix isolation experiments have allowed the measurement of the UV absorption spectra of the high-energy non-chelated isomers of acetylacetone, these isomers being produced by UV irradiation of the stable chelated form. Their identification has been done by coupling selective UV-induced isomerization, infrared spectroscopy, and harmonic vibrational frequency calculations using density functional theory. The relative energies of the chelated and non-chelated forms of acetylacetone in the S0 state have been obtained using density functional theory and coupled-cluster methods. For each isomer of acetylacetone, we have calculated the UV transition energies and dipole oscillator strengths using the excited-state coupled-cluster methods, including EOMCCSD (equation-of-motion coupled-cluster method with singles and doubles) and CR-EOMCCSD(T) (the completely renormalized EOMCC approach with singles, doubles, and non-iterative triples). For dipole-allowed transition energies, there is a very good agreement between experiment and theory. In particular, the CR-EOMCCSD(T) approach explains the blue shift in the electronic spectrum due to the formation of the non-chelated species after the UV irradiation of the chelated form of acetylacetone. Both experiment and CR-EOMCCSD(T) theory identify two among the seven non-chelated forms to be characterized by red-shifted UV transitions relative to the remaining five non-chelated isomers.     

High-accuracy theoretical thermochemistry of atmospherically important nitrogen oxide derivatives

High-accuracy quantum chemical calculations were performed for several atmospherically important nitrogen oxide derivatives, such as HOONO, HOONO(2), NH(2)NO(2), FNO, FNO(2), FONO, FONO(2), ClNO, ClONO, ClONO(2), and ClOONO. The stable conformers of the molecules were identified, and the corresponding heats of formation (Δ(f)H(0)° and Δ(f)H(298)°) and entropy values (S(298)°) were computed. On the basis of the thermodynamic functions, equilibrium constants were also calculated for a couple of reactions with importance in the chemistry of the atmosphere. In a number of cases this study provides more reliable estimates for the investigated thermodynamic properties than those can be collected from previous reports.     

ClONO2异构化反应和分解反应的理论研究

The potential energy surface for the decomposition and isomerization of Chlorine Nitrate（ClONO2）is calculated using the G3 theory. Geometries of related species are optimized at the MP2（ full）/6-31+G（d）and B3LYP/6-31+G（d）level. Vibrational frequencies and IR intensities of ClONO2 have been calculated at the same level. Obtained geometries and vibrational frequencies as well as IR intensities for ClONO2 are in good agreement with experimental values. A new stable stereoisomer is verified by CCSD（T）and QCISD（T）methods at 6-311G （d）basis set. Calculated geometries using above methods for this stereoisomer are in accord with each other. The calculated reaction heat of ClONO2 are also in good agreement with the available data in the literature. Among these reactions of ClONO2,the stereoisomer reaction is the hardest one. The barrier height for this reaction is 481.52 kJ/mol relative to ClONO2 and the reaction is endothermic by 299.85 kJ/mol. On the other hand,among these unimolecular reactions of ClONO2,the decomposition reaction of NO2+ClO is the easiest one. These results indicate that ClONO2 is very stable.     

High level ab initio study of the structure and vibrational spectra of HO(2)NO(2)

A high-level ab initio study has been performed on the conformational structure and vibrational spectra of HO(2)NO(2). Calculations carried out with coupled-cluster methods using a series of Pople and Dunning basis sets reveal that there is a significant basis set dependence on the predicted ab initio structure. Higher angular momentum basis sets are shown to be necessary in order to bring the calculated structure into agreement with experimental rotational constants. Harmonic vibrational frequencies of HO(2)NO(2) are computed at the CCSD(T)/aug-cc-pVTZ level of theory while the corresponding vibrational anharmonicities are calculated at the MP2/cc-pVTZ level. In addition, the absorption cross sections of OH stretching overtones in HO(2)NO(2) are calculated using a dipole function computed at the QCISD level of theory and found to be in good agreement with the available experimental data.     

Ab initio characterization of the Ne-I2 van der Waals complex: intermolecular potentials and vibrational bound states

A theoretical study of the potential energy surface and bound states is performed for the ground state of the NeI(2) van der Waals (vdW) complex. The three-dimensional interaction energies are obtained from ab initio coupled-cluster, coupled-cluster single double (triple)/complete basis set, calculations using large basis sets, of quadruple- through quintuple-zeta quality, in conjunction with relativistic effective core potentials for the heavy iodine atoms. For the analytical representation of the surface two different schemes, based on fitting and interpolation surface generation techniques, are employed. The surface shows a double-minimum topology for linear and T-shaped configurations. Full variational quantum mechanical calculations are carried out using the model surfaces, and the vibrationally averaged structures and energetics for the NeI(2) isomers are determined. The accuracy of the potential energy surfaces is validated by a comparison between the present results and the corresponding experimental data available. In lieu of more experimental measurements, we also report our results/predictions on higher bound vibrational vdW levels, and the influence of the employed surface on them is discussed.     

Toward spectroscopic accuracy for open-shell systems: molecular structure and hyperfine coupling constants of H2CN, H2CP, NH2, and PH2 as test cases

In the present paper, we investigate the molecular structure and hyperfine couplings of a series of σ radicals containing first- and second-row atoms (H(2)CN,H(2)CP,NH(2),PH(2)) for which accurate gas-phase microwave results are available. The presence of α- and, especially, β-hydrogen atoms makes the evaluation of magnetic properties of these radicals particularly challenging. Geometrical parameters have been computed by the coupled-cluster ansatz in conjunction with hierarchical series of basis sets, thus accounting for extrapolation to the complete basis-set limit. Core correlation as well as higher excitations in the electronic-correlation treatment have also been taken into account. An analogous approach has been employed for evaluating hyperfine coupling constants with particular emphasis given to basis-set, correlation, and geometrical effects. The corresponding vibrational corrections, required for a meaningful comparison to experimental data, have also been investigated. The remarkable agreement with experiment confirms the reliability of the present computational approach, already validated for π radicals, thus establishing the way for setting up a benchmark database for magnetic properties.     

Nitryl chloride (ClNO2): UV/vis absorption spectrum between 210 and 296 K and O(3P) quantum yield at 193 and 248 nm

Recent studies have shown that the UV/vis photolysis of nitryl chloride (ClNO2) can be a major source of reactive chlorine in the troposphere. The present work reports measurements of the ClNO2 absorption spectrum and its temperature dependence between 210 and 296 K over the wavelength range 200–475 nm using diode array spectroscopy. The room temperature spectrum obtained in this work was found to be in good agreement with the results from Ganske et al. (J. Geophys. Res. 1992, 97, 7651) over the wavelength range common to both studies (200–370 nm) but differs systematically from the currently recommended spectrum for use in atmospheric models. The present results lead to a decrease in the calculated atmospheric ClNO2 photolysis rate by 30%. Including the temperature dependence of the ClNO2 spectrum decreases the calculated atmospheric photolysis rate at lower temperatures (higher altitudes) even further. A parametrization of the wavelength and temperature dependence of the ClNO2 spectrum is presented. O(3P) quantum yields, Φ(ClNO2)(O), in the photolysis of ClNO2 at 193 and 248 nm were measured at 296 K using pulsed laser photolysis combined with atomic resonance fluorescence detection of O(3P) atoms. Φ(ClNO2)(O)(λ) was found to be 0.67 ± 0.12 and 0.15 ± 0.03 (2σ error limits, including estimated systematic errors) at 193 and 248 nm, respectively, indicating that multiple dissociation channels are active in the photolysis of ClNO2 at these wavelengths. The Φ(ClNO2)(O)(λ) values obtained in this work are discussed in light of previous ClNO2 photodissociation studies and the differences are discussed.