Theoretical study of the rovibrational spectrum of H2O-H2

In this paper we report transition frequencies and line strengths computed for H(2)O-H(2) and compare with the experimental observations of [M. J. Weida and D. J. Nesbitt, J. Chem. Phys. 110, 156 (1999)]. To compute the spectra we use a symmetry adapted Lanczos algorithm and an uncoupled product basis set. Our results corroborate the assignments of Weida and Nesbitt and there is good agreement between calculated and observed transitions. Possible candidates for lines that Weida and Nesbitt were not able to assign are presented. Several other bands that may be observable are also discovered. Although all the observed bands are associated with states localized near the global potential minimum, at which H(2)O acts as proton acceptor, a state with significant amplitude near the T-shape secondary potential minimum at which H(2)O acts as proton donor is identified by examining many different probability density plots.     

Calculation of argon trimer rovibrational spectrum

Rovibrational spectra of Ar3 are computed for total angular momenta up to J=6 using row-orthonormal hyperspherical coordinates and an expansion of the wave function on hyperspherical harmonics. The sensitivity of the spectra to the two-body potential and to the three-body corrections is analyzed. First, the best available semiempirical pair potential (HFDID1) is compared with our recent ab initio two-body potential. The ab initio vibrational energies are typically 1-2 cm-1 higher than the semiempirical ones, which is related to the slightly larger dissociation energy of the semiempirical potential. Then, the Axilrod-Teller asymptotic expansion of the three-body correction is compared with our newly developed ab initio three-body potential. The difference is found smaller than 0.3 cm-1. In addition, we define approximate quantum numbers to describe the vibration and rotation of the system. The vibration is represented by a hyper-radial mode and a two-degree-of-freedom hyperangular mode, including a vibrational angular momentum defined in an Eckart frame. The rotation is described by the total angular momentum quantum number, its projection on the axis perpendicular to the molecular plane, and a hyperangular internal momentum quantum number, related to the vibrational angular momentum by a transformation between Eckart and principal-axes-of-inertia frames. These quantum numbers provide a qualitative understanding of the spectra and, in particular, of the impact of the nuclear permutational symmetry of the system (bosonic with zero nuclear spin). Rotational constants are extracted from the spectra and are shown to be accurate only for the ground hyperangular mode.     

Theoretical and experimental studies of the infrared rovibrational spectrum of He2-N2O

Rovibrational spectra of the He(2)-N(2)O complex in the nu(1) fundamental band of N(2)O (2224 cm(-1)) have been observed using a tunable infrared laser to probe a pulsed supersonic jet expansion, and calculated using five coordinates that specify the positions of the He atoms with respect to the NNO molecule, a product basis, and a Lanczos eigensolver. Vibrational dynamics of the complex are dominated by the torsional motion of the two He atoms on a ring encircling the N(2)O molecule. The resulting torsional states could be readily identified, and they are relatively uncoupled to other He motions up to at least upsilon(t) = 7. Good agreement between experiment and theory was obtained with only one adjustable parameter, the band origin. The calculated results were crucial in assigning many weaker observed transitions because the effective rotational constants depend strongly on the torsional state. The observed spectra had effective temperatures around 0.7 K and involved transitions with J < or =3, with upsilon(t) = 0 and 1, and (with one possible exception) with Deltaupsilon(t)=0. Mixing of the torsion-rotation states is small but significant: some transitions with Deltaupsilon(t) not equal 0 were predicted to have appreciable intensity even assuming that the dipole transition moment coincides perfectly with the NNO axis. One such transition was tentatively assigned in the observed spectra, but confirmation will require further work.     

Anharmonic rovibrational analysis for disilacyclopropenylidene (Si2CH2)

The global minimum on the Si(2)CH(2) electronic singlet potential energy surface has been theoretically predicted to be a peculiar hydrogen bridged (Si···H···Si) disilacyclopropenylidene structure (Si(2)CH(2)). An accurate quartic force field for Si(2)CH(2) has been determined employing ab initio coupled-cluster theory with single and double excitations and a perturbative treatment for triple excitations [CCSD(T)], in combination with the correlation consistent core-valence quadruple zeta (cc-pCVQZ) basis set. The vibration-rotation coupling constants, equilibrium and zero-point vibration corrected rotational constants, centrifugal distortion constants, and harmonic and fundamental vibrational frequencies for six isotopologues of Si(2)CH(2) are predicted using vibrational second-order perturbation theory (VPT2). The anharmonic corrections for the vibrational motions involving the H bridged bonds are found to be more than 5% with respect to the corresponding harmonic vibrational frequencies. In this light, an experimental detection and characterization of disilacyclopropenylidene (Si(2)CH(2)) is highly desired.     

Isotope effects in the CO dimer: millimeter wave spectrum and rovibrational calculations of (12C18O)2

The millimeter wave spectrum of the isotopically substituted CO dimer, (12C18O)2, was studied with the Orotron jet spectrometer, confirming and extending a previous infrared study [A. R. W. McKellar, J. Mol. Spectrosc. 226, 190 (2004)]. A very dilute gas mixture of CO in Ne was used, which resulted in small consumption of 12C18O sample gas and produced cold and simple spectra. Using the technique of combination differences together with the data from the infrared work, six transitions in the 84-127 GHz region have been assigned. They belong to two branches, which connect four low levels of A+ symmetry to three previously unknown levels of A- symmetry. The discovery of the lowest state of A- symmetry, which corresponds to the projection K=0 of the total angular momentum J onto the intermolecular axis, identifies the geared bending mode of the 12C18O dimer at 3.607 cm(-1). Accompanying rovibrational calculations using a recently developed hybrid potential from ab initio coupled cluster [CCSD(T)] and symmetry-adapted perturbation theory calculations [G. W. M. Vissers et al., J. Chem. Phys. 122, 054306 (2005)] gave very good agreement with experiment. The isotopic dependence of the A+/A- energy splitting, the intermolecular separation R, and the energy difference of two ground state isomers, which change significantly when 18O or 13C are substituted into the normal (12C16O)2 isotopolog [L. A. Surin et al., J. Mol. Spectrosc. 223, 132 (2004)], was explained by these calculations. It turns out that the change in anisotropy of the intermolecular potential with respect to the shifted monomer centers of mass is particularly significant.     

Potential energy surface and quantum dynamics study of rovibrational states for HO(3) (X (2)A')

An analytic potential energy surface has been constructed by fitting to about 28 thousand energy points for the electronic ground-state (X (2)A') of HO(3). The energy points are calculated using a hybrid density functional HCTH and a large basis set aug-cc-pVTZ, i.e., a HCTH/aug-cc-pVTZ density functional theory (DFT) method. The DFT calculations show that the trans-HO(3) isomer is the global minimum with a potential well depth of 9.94 kcal mol(-1) with respect to the OH + O(2) asymptote. The equilibrium geometry of the cis-HO(3) conformer is located 1.08 kcal mol(-1) above that of the trans-HO(3) one with an isomerization barrier of 2.41 kcal mol(-1) from trans- to cis-HO(3). By using this surface, a rigorous quantum dynamics (QD) study has been carried out for computing the rovibrational energy levels of HO(3). The calculated results determine a dissociation energy of 6.15 kcal mol(-1), which is in excellent agreement with the experimental value of Lester et al. [J. Phys. Chem. A, 2007, 111, 4727.].     

Infrared spectra of seeded hydrogen clusters: (paraH2)N - OCS, (orthoH2)N - OCS, and (HD)N - OCS, N = 2 - 7

Infrared spectra of hydrogen-carbonyl sulfide clusters containing paraH2, orthoH2, or HD have been studied in the 2060 cm(-1) region of the C-O stretching vibration. The clusters were formed in pulsed supersonic jet expansions and probed using a tunable infrared diode laser spectrometer. Simple symmetric rotor type spectra were observed and assigned for clusters containing up to N = 7 hydrogen molecules. There was no resolved K structure, and Q-branch features were present for orthoH2 and HD but absent for paraH2. These characteristics can be rationalized in terms of near symmetric rotor structures, very low effective rotational temperatures (0.15 to 0.6 K), and nuclear spin statistics. The observed vibrational shifts were compared with those from recent observations on the same clusters embedded in helium nanodroplets. The observed rotational constants for the paraH2 clusters are in good agreement with a recent quantum Monte Carlo simulation. Some mixed clusters were also observed, such as HD-HD-He-OCS and paraH2 - orthoH2 - OCS.     

Slice imaging of the quantum state-to-state cross section for photodissociation of state-selected rovibrational bending states of OCS (v2=0,1,2/JlM)+hnu-->CO(J)+S(1D2)

Using hexapole quantum state-selection of OCS (v(2)=0,1,2/JlM) and high-resolution slice imaging of quantum state-selected CO(J), the state-to-state cross section OCS (v(2)=0,1,2/JlM)+hnu-->CO(J)+S((1)D(2)) was measured for bending states up to v(2)=2. The population density of the state-selected OCS (v(2)=0,1,2 /JlM) in the molecular beam was obtained by resonantly enhanced multiphoton ionization of OCS and comparison with room temperature bulk gas. A strong increase of the cross section with increasing bending state is observed for CO(J) in the high J region, J=60-67. Integrating over all J states the authors find sigma(v(2)=0):sigma(v(2)=1):sigma(v(2)=2)=1.0:7.0:15.0. A quantitative comparison is made with the dependence of the transition dipole moment function on the bending angle.     

A CASPT2 study of the electronic spectrum of hexacyanoosmate(III)

CASPT2 calculations reveal that the ligand field splitting parameter Delta(o) of [Os(CN)6]3- is much higher than previously proposed values of +/-38,000 cm(-1). In line with the expected increase down a transition-metal group, Delta(o) is found to be +/-55,000 cm(-1), excluding the possible appearance of ligand field transitions in the UV-vis spectrum. Instead, the calculations confirm that the observed spectrum arises from the three lowest symmetry-allowed ligand-to-metal charge-transfer (LMCT) excitations. Spin-orbit coupling in the ground state is found to be about 4350 cm(-1), leading to a spin-orbit coupling constant zeta of +/-2900 cm(-1). Spin-orbit coupling in the 2T(1u) LMCT states is found not to be negligible, contrary to previous belief.     

The Al+ -H(2) cation complex: rotationally resolved infrared spectrum, potential energy surface, and rovibrational calculations

The infrared spectrum of the Al(+)-H(2) complex is recorded in the H-H stretch region (4075-4110 cm(-1)) by monitoring Al(+) photofragments. The H-H stretch band is centered at 4095.2 cm(-1), a shift of -66.0 cm(-1) from the Q(1)(0) transition of the free H(2) molecule. Altogether, 47 rovibrational transitions belonging to the parallel K(a)=0-0 and 1-1 subbands were identified and fitted using a Watson A-reduced Hamiltonian, yielding effective spectroscopic constants. The results suggest that Al(+)-H(2) has a T-shaped equilibrium configuration with the Al(+) ion attached to a slightly perturbed H(2) molecule, but that large-amplitude intermolecular vibrational motions significantly influence the rotational constants derived from an asymmetric rotor analysis. The vibrationally averaged intermolecular separation in the ground vibrational state is estimated as 3.03 A, decreasing by 0.03 A when the H(2) subunit is vibrationally excited. A three-dimensional potential energy surface for Al(+)-H(2) is calculated ab initio using the coupled cluster CCSD(T) method and employed for variational calculations of the rovibrational energy levels and wave functions. Effective dissociation energies for Al(+)-H(2)(para) and Al(+)-H(2)(ortho) are predicted, respectively, to be 469.4 and 506.4 cm(-1), in good agreement with previous measurements. The calculations reproduce the experimental H-H stretch frequency to within 3.75 cm(-1), and the calculated B and C rotational constants to within approximately 2%. Agreement between experiment and theory supports both the accuracy of the ab initio potential energy surface and the interpretation of the measured spectrum.     

Fully quantum rovibrational calculation of the He(H2) bound and resonance states

In a recent paper [J. Chem. Phys. 2005, 122, 124318], a full-dimensional quantum method, designed to efficiently compute the rovibrational states of triatomic systems with long-range interactions, was applied to the benchmark Li-(H2) ion-molecule system. The method incorporates several key features in order to accurately represent the rovibrational Hamiltonian using only modestly sized basis sets: (1) exact analytical treatment of Coriolis coupling; (2) a single bend-angle basis for all rotational states; (3) phase space optimization of the vibrational basis; (4) G(4) symmetry adaptation of the rovibrational basis. In this paper, the same methodology is applied for the first time to a van der Waals complex system, He(H2). As in the Li-(H2) study, all of the rovibrational bound states, and a number of resonance states, are computed to very high accuracy (1/10,000 of a wavenumber or better). Three different isotopologues are considered, all of which are found to have a single bound state with a very low binding energy. Several extremely long-lived Feshbach resonances are also reported.     

Ab initio potential energy surface and rovibrational spectrum of Ar-HCCCN

We report an ab initio intermolecular potential energy surface of the Ar-HCCCN complex using a supermolecular method. The calculations were performed using the fourth-order M?ller-Plesset theory with the full counterpoise correction for the basis set superposition error and a large basis set including bond functions. The complex was found to have a planar T-shaped structure minimum and a linear minimum with the Ar atom facing the H atom. The T-shaped minimum is the global minimum with the well depth of 236.81 cm(-1). A potential barrier separating the two minima is located at R=5.57 A and theta=20.39 degrees with the height of 151.59 cm(-1). The two-dimensional discrete variable representation was employed to calculate the rovibrational energy levels for Ar-HCCCN. The rovibrational spectra including intensities for the ground state and the first excited intermolecular vibrational state are also presented. The results show that the spectra are mostly b-type (Delta K(a)=+/-1) transitions with weak a-type (Delta K(a)=0) transitions in structure, which are in good agreement with the recent experimental results [A. Huckauf, W. Jager, P. Botschwina, and R. Oswald, J. Chem. Phys. 119, 7749 (2003)].     

Accurate quantum calculation of the bound and resonant rovibrational states of Li-(H2)

In a recent paper [B. Poirier, Chem. Phys. 308, 305 (2005)] a full-dimensional quantum method for computing the rovibrational dynamics of triatomic systems was presented, incorporating three key features: (1) exact analytical treatment of Coriolis coupling, (2) three-body "effective potential," and (3) a single bend angle basis for all rotational states. In this paper, these ideas are applied to the Li-(H2) electrostatic complex, to compute all of the rovibrational bound state energies, and a number of resonance energies and widths, to very high accuracy (thousandths of a wave number). This application is very challenging, owing to the long-range nature of the interaction and to narrow level spacings near dissociation. Nevertheless, by combining the present method with a G4 symmetry-adapted phase-space-optimized representation, only modest basis sizes are required for which the matrices are amenable to direct diagonalization. Several new bound levels are reported, as compared with a previous calculation [D. T. Chang, G. Surratt, G. Ristroff, and G. I. Gellene, J. Chem. Phys. 116, 9188 (2002)]. The resonances exhibit a clear-cut separation into shape and Feshbach varieties, with the latter characterized by extremely long lifetimes (microseconds or longer).     

Computational study of ethanol adsorption and reaction over rutile TiO(2) (110) surfaces

Studies of the modes of adsorption and the associated changes in electronic structures of renewable organic compounds are needed in order to understand the fundamentals behind surface reactions of catalysts for future energies. Using planewave density functional theory (DFT) calculations, the adsorption of ethanol on perfect and O-defected TiO(2) rutile (110) surfaces was examined. On both surfaces the dissociative adsorption mode on five-fold coordinated Ti cations (Ti(4+)(5c)) was found to be more favourable than the molecular adsorption mode. On the stoichiometric surface E(ads) was found to be equal to 0.85 eV for the ethoxide mode and equal to 0.76 eV for the molecular mode. These energies slightly increased when adsorption occurred on the Ti(4+)(5c) closest to the O-defected site. However, both considerably increased when adsorption occurred at the removed bridging surface O; interacting with Ti(3+) cations. In this case the dissociative adsorption becomes strongly favoured (E(ads) = 1.28 eV for molecular adsorption and 2.27 eV for dissociative adsorption). Geometry and electronic structures of adsorbed ethanol were analysed in detail on the stoichiometric surface. Ethanol does not undergo major changes in its structure upon adsorption with its C-O bond rotating nearly freely on the surface. Bonding to surface Ti atoms is a σ type transfer from the O2p of the ethanol-ethoxide species. Both ethanol and ethoxide present potential hole traps on O lone pairs. Charge density and work function analyses also suggest charge transfer from the adsorbate to the surface, in which the dissociative adsorptions show a larger charge transfer than the molecular adsorption mode.     

A theoretical study on the ionization of OCS with an analysis of vibrational structures of the photoelectron spectrum

 Ab initio calculations have been performed to study the molecular structures and vibrational levels of the four low-lying ionic states (1, 2²Π, and 1, 2²Σ⁺) of carbonyl sulfide. The global regions of the potential-energy surfaces have been obtained by multireference single and double excitation configuration interaction calculations. Vibrational calculations using explicit vibrational Hamiltonians have been used for vibrational analysis. The equilibrium molecular structures and a vibrational analysis of the four ionic states are presented. The theoretical ionization intensity curves including the vibrational structures of the ionic states are also presented and are compared with the photoelectron spectrum. Received: 20 January 2001 / Accepted: 22 August 2001 / Published online: 30 October 2001     

Computational study of the thermochemistry of organophosphorus(III) compounds

The enthalpies of formation of organophosphorus(III) compounds have been calculated at the G3X, G3X(MP2), and B3LYP/6-311+G(3df,2p)//B3LYP/6-31G(d,p) levels of theory using the atomization energy procedure and the method of isodesmic reactions. The Delta f H298 degree values for 50 relatively large molecules with up to 10 non-hydrogen atoms, such as P(CH3)3, P(C2H5)3, P(OCH3)3, n-C4H9OPCl2, [(CH3)2N]2PCl, (C2H5)2NPCl2, and [(CH3)2N]2PCN, have been calculated directly from the G3X atomization energies. A good agreement between the known experimental values and G3X results for 14 compounds provides support to our predictions for remaining species whose experimental enthalpies of formation are unknown or known with relatively large uncertainties. On the basis of our calculations and sometimes conflicting experimental data a set of internally consistent enthalpies of formation has been recommended for organophosphorus(III) compounds. Our computational results call into question the experimental enthalpies of formation of P(C2H5)3 and P(n-C4H9)3. From comparison with most reliable experimental data, the accuracy of the theoretical enthalpies of formation is estimated as ranging from 5 to 10 kJ/mol. The recommended Delta f H298 degree values were used to derive the group additivity values (GAVs) for 45 groups involving the phosphorus(III) atom. These GAVs significantly extend the applicability of Benson's group additivity method and may be used to estimate the enthalpies of formation of larger organophosphorus(III) compounds, where high level quantum chemical calculations are impracticable.     

Computational study on the existence of organic peroxy radical-water complexes (RO2.H2O)

The existence of a series of organic peroxy radical-water complexes [CH3O2.H2O (methyl peroxy); CH3CH2O2.H2O (ethyl peroxy); CH3C(O)O2.H2O (acetyl peroxy); CH3C(O)CH2O2.H2O (acetonyl peroxy); CH2(OH)O2.H2O (hydroxyl methyl peroxy); CH2(OH)CH2O2.H2O (2-hydroxy ethyl peroxy); CH2(F)O2.H2O (fluoro methyl peroxy); CH2(F)CH2O2.H2O (2-fluoro ethyl peroxy)] is evaluated using high level ab initio calculations. A wide range of binding energies is predicted for these complexes, in which the difference in binding energies can be explained by examination of the composition of the R group attached to the peroxy moiety. The general trend in binding energies has been determined to be as follows: fluorine approximately alkyl < carbonyl < alcohol. The weakest bound complex, CH3O2.H2O, is calculated to be bound by 2.3 kcal mol-1, and the strongest, the CH2(OH)O2.H2O complex, is bound by 5.1 kcal mol-1. The binding energy of the peroxy radical-water complexes which contain carbonyl and alcohol groups indicates that these complexes may perturb the kinetics and product branching ratios of reactions involving these complexes.     

Two Modifications of the Mixed-Ligand Complex ZnPhen(i-C3H7OCS2)2 with Monodentate and Bidentate Ligands i-C3H7OCS 2 -

Two crystal modifications (m-1 and m-2) of the mixed-ligand complex ZnPhen(i-PrOCS₂)₂ have been isolated. Their crystal structures were solved using X-ray diffraction data (CAD-4 diffractometer, MoK radiation, 3141 and 3532 F_hkl, R = 0.0363 and 0.0304). For both modifications, the crystals are monoclinic with unit cell parameters a = 10.543(2), b = 13.494(3), c = 16.875(3) , = 102.08(3)°, V = 2347.6(8) ³, Z = 4, space group P2₁/n (m-1) and a = 10.931(2), b = 12.996(3), c = 16.288(3) , = 92.69(3)°, V = 2311.3(8) ³, Z = 4, space group P2₁/n (m-2). The structures basically consist of discrete monomer molecules in which the Zn atom has a tetragonal pyramidal (m-1: ZnN₂S₃, c.n. 5, bidentate and monodentate i-PrOCS ₂ ^- ligands) or distorted octahedral (m-2: ZnN₂S₄, c.n. 6, bidentate i-PrOCS ₂ ^- ligands) environment. Molecular packings and their interactions in the structures are discussed.     

Electronic spectrum of the AlC(2) radical

An electronic transition of the AlC2 radical (C2v structure) has been observed using laser-induced fluorescence spectroscopy. The molecule was prepared in a supersonic expansion by ablation of an aluminum rod in the presence of acetylene gas. A spectrum was recorded in the 451-453 nm region and assigned to the C 2B2-X 2A1 system (T0 = 22,102.7 cm(-1)) based on a rotational analysis and agreement with calculated molecular parameters and excitation energies. Ab initio results obtained using couple cluster methods are in accord with previous theoretical work which concludes that ground-state AlC2 possesses a T-shaped C2v 2A1 geometry, with the linear 2Sigma+ AlCC isomer 0.70 eV higher in energy. A fit of the experimental spectrum yields rotational constants in the ground and electronically excited states that are in reasonable agreement with the calculated values: A' = 1.7093(107), B' = 0.4052(50), C' = 0.3228(49) cm(-1) for the X 2A1 state, and A' = 1.5621(137), B' = 0.4028(46), C' = 0.3201(54) cm(-1) for C 2B2. Variation in individual fluorescence lifetimes suggests that the emitting C 2B2 state undergoes rovibronic mixing with lower lying electronic states.     

Computational study of the aminolysis of 2-benzoxazolinone

Three possible mechanisms (zwitterionic, neutral stepwise, and neutral concerted) of the ring-opening reaction of 2-benzoxazolinone (BO) upon aminolysis with methylamine were studied at the B3LYP/6-31G* level. In the gas phase, the neutral concerted mechanism is shown to be most favorable, which proceeds via a rate-determining barrier of 28-29 kcal/mol. The transition state, CTS, associated with this barrier is a four-centered one, where 1,2-addition of the N[bond]H of methylamine to the C[bond]O of BO ring occurs. The rate-determining barrier of the neutral stepwise pathway is found to be ca. 42 kcal/mol. The inclusion of solvent effects by a polarizable continuum model (PCM) does not change the conclusions based on the gas-phase study; the barrier at CTS is reduced to 20, 20, and 22 kcal/mol in water, ethanol, and acetonitrile, respectively.