Scattering lengths for Li-Cs, Na-Cs, K-Cs, and Rb-Cs ultracold collisions

We present triplet and singlet scattering lengths for X--133Cs collisions (X=6,7Li, 23Na, 39,40,41K, and 85,87Rb). We consider the short-range potential-energy curves based on high-quality ab initio calculations and the long-range potential described by accurate dispersion coefficients available in the literature. Using a Fermi function to smoothly connect the short- and long-range terms we construct the potential for all R and evaluate the scattering length for each mixed species X--Cs. In particular, we predict that the cross section for inelastic processes between 23Na and 133Cs and between 40K and 133Cs will be small as compared with the respective elastic cross sections.     

Low energy collisions of CN(X2Σ+) with He in magnetic fields

A theoretical investigation of the He-CN((2)Σ(+)) complex is presented. We perform ab initio calculations of the interaction potential energy surface and carry out accurate calculations of bound energy levels of the complex including the molecular fine structure. We find the potential has a shallow minimum and supports seven and nine bound levels in complex with (3)He and (4)He, respectively. Based on the potential the quantum scattering calculation is then implemented for elastic and inelastic cross sections of the magnetically trappable low-field-seeking state of CN((2)Σ(+)) in collision with (3)He atom. The cold collision properties and the influence of the external magnetic field as well as the effect of the uncertainty of interaction potential on the collisionally induced Zeeman relaxation are explored and discussed in detail. The ratios of elastic to inelastic cross sections are large over a wide range of collision energy, magnetic field, and scaling factor of the potential, suggesting helium buffer gas loading and cooling of CN in a magnetic trap is a good prospect.     

Scattering resonances in slow NH3-He collisions

We theoretically study slow collisions of NH(3) molecules with He atoms, where we focus in particular on the observation of scattering resonances. We calculate state-to-state integral and differential cross sections for collision energies ranging from 10(-4) cm(-1) to 130 cm(-1), using fully converged quantum close-coupling calculations. To describe the interaction between the NH(3) molecules and the He atoms, we present a four-dimensional potential energy surface, based on an accurate fit of 4180 ab initio points. Prior to collision, we consider the ammonia molecules to be in their antisymmetric umbrella state with angular momentum j = 1 and projection k = 1, which is a suitable state for Stark deceleration. We find pronounced shape and Feshbach resonances, especially for inelastic collisions into the symmetric umbrella state with j = k = 1. We analyze the observed resonant structures in detail by looking at scattering wavefunctions, phase shifts, and lifetimes. Finally, we discuss the prospects for observing the predicted scattering resonances in future crossed molecular beam experiments with a Stark-decelerated NH(3) beam.     

New exchange-Coulomb N2-Ar potential-energy surface and its comparison with other recent N2-Ar potential-energy surfaces

The reliability of five N2-Ar potential-energy surfaces in representing the N2-Ar interaction has been investigated by comparing their abilities to reproduce a variety of experimental results, including interaction second viral coefficients, bulk transport properties, relaxation phenomena, differential scattering cross sections, and the microwave and infrared spectra of the van der Waals complexes. Four of the surfaces are the result of high-level ab initio quantal calculations; one of them utilized fine tuning by fitting to microwave data. To date, these four potential-energy surfaces have only been tested against experimental microwave data. The fifth potential-energy surface, based upon the exchange-Coulomb potential-energy model for the interaction of closed-shell species, is developed herein: it is a combination of a damped dispersion energy series and ab initio calculations of the Heitler-London interaction energy, and has adjustable parameters determined by requiring essentially simultaneous agreement with selected quality interaction second viral coefficient and microwave data. Comparisons are also made with the predictions of three other very good literature potential-energy surfaces, including the precursor of the new exchange-Coulomb potential-energy surface developed here. Based upon an analysis of a large body of information, the new exchange-Coulomb and microwave-tuned ab initio potential-energy surfaces provide the best representations of the N2-Ar interaction; nevertheless, the other potential-energy surfaces examined still have considerable merit with respect to the prediction of specific properties of the N2-Ar van der Waals complex. Of the two recommended surfaces, the new exchange-Coulomb surface is preferred on balance due to its superior predictions of the effective cross sections related to various relaxation phenomena, and to its reliable, and relatively simple, representation of the long-range part of the potential-energy surface. Moreover, the flexibility still inherent in the exchange-Coulomb potential form can be further exploited, if required, in future studies of the N2-Ar interaction.     

Classical trajectory study of collisions of Ar with alkanethiolate self-assembled monolayers: potential-energy surface effects on dynamics

We have investigated collisions between Ar and alkanethiolate self-assembled monolayers (SAMs) using classical trajectory calculations with several potential-energy surfaces. The legitimacy of the potential-energy surfaces is established through comparison with molecular-beam data and ab initio calculations. Potential-energy surfaces used in previous work overestimate the binding of Ar to the SAM, leading to larger energy transfer than found in the experiments. New calculations, based on empirical force fields that better reproduce ab initio calculations, exhibit improved agreement with the experiments. In particular, polar-angle-dependent average energies calculated with explicit-atom potential-energy surfaces are in excellent agreement with the experiments. Polar- and azimuthal-angle-dependent product translational energies are examined to gain deeper insight into the dynamics of Ar+SAM collisions.     

Efficient method for quantum calculations of molecule-molecule scattering properties in a magnetic field

We show that the cross sections for molecule-molecule collisions in the presence of an external field can be computed efficiently using a total angular momentum basis, defined either in the body-fixed frame or in the space-fixed coordinate system. This method allows for computations with much larger basis sets than previously possible. We present calculations for (15)NH-(15)NH collisions in a magnetic field. Our results support the conclusion of the previous study that the evaporative cooling of rotationally ground (15)NH molecules in a magnetic trap has a prospect of success.     

Ultracold O2 + O2 collisions in a magnetic field: on the role of the potential energy surface

The collision dynamics of (17)O(2)((3)Σ(g)(-)) + (17)O(2)((3)Σ(g)(-)) in the presence of a magnetic field is studied within the close-coupling formalism in the range between 10 nK and 50 mK. A recent global ab initio potential energy surface (PES) is employed and its effect on the dynamics is analyzed and compared with previous calculations where an experimentally derived PES was used [T. V. Tscherbul et al., New J. Phys 134, 055021 (2009)]. Compared to the results using the older PES, magnetic-field dependence of the low-field-seeking state in the ultracold regime is characterized by a very large background scattering length, a(bg), and cross sections exhibit broader and more pronounced Feshbach resonances. The marked resonance structure is somewhat surprising considering the influence of inelastic scattering but it can be explained by resorting to the analytical van der Waals theory, where the short-range amplitude of the entrance channel wavefunction is enhanced by the large a(bg). This strong sensitivity to the short range of the ab initio PES persists up to relatively high energies (10 mK). After this study and despite quantitative predictions are very difficult, it can be concluded that the ratio between elastic and spin relaxation scattering is generally small, except for magnetic fields which are either low or close to an asymmetric Fano-type resonance. Some general trends found here, such as a large density of quasibound states and a propensity toward large scattering lengths, could be also characteristic of other anisotropic molecule-molecule systems.     

Potential-energy surface for the electronic ground state of NH3 up to 20,000 cm-1 above equilibrium

Ab initio coupled cluster calculations with single and double substitutions and a perturbative treatment of connected triple excitations [CCSD(T)] with the augmented correlation-consistent polarized valence triple-zeta aug-cc-pVTZ basis at 51 816 geometries provide a six-dimensional potential-energy surface for the electronic ground state of NH3. At 3814 selected geometries, CBS+ energies are obtained by extrapolating the CCSD(T) results for the aug-cc-pVXZ(X=T,Q,5) basis sets to the complete basis set (CBS) limit and adding corrections for core-valence correlation and relativistic effects. CBS** ab initio energies are generated at 51,816 geometries by an empirical extrapolation of the CCSD(T)/aug-cc-pVTZ results to the CBS+ limit. They cover the energy region up to 20,000 cm-1 above equilibrium. Parametrized analytical functions are fitted through the ab initio points. For these analytical surfaces, vibrational term values and transition moments are calculated by means of a variational program employing a kinetic-energy operator expressed in the Eckart-Sayvetz frame. Comparisons against experiment are used to assess the quality of the generated potential-energy surfaces. A "spectroscopic" potential-energy surface of NH3 is determined by a slight empirical adjustment of the ab initio potential to the experimental vibrational term values. Variational calculations on this refined surface yield rms deviations from experiment of 0.8 cm-1 for 24 inversion splittings and 0.4 (3.0) cm-1 for 34 (51) vibrational term values up to 6100 (10,300) cm-1.     

Quantum theory of molecular collisions in a magnetic field: efficient calculations based on the total angular momentum representation

An efficient method is presented for rigorous quantum calculations of atom-molecule and molecule-molecule collisions in a magnetic field. The method is based on the expansion of the wave function of the collision complex in basis functions with well-defined total angular momentum in the body-fixed coordinate frame. We outline the general theory of the method for collisions of diatomic molecules in the (2)Σ and (3)Σ electronic states with structureless atoms and with unlike (2)Σ and (3)Σ molecules. The cross sections for elastic scattering and Zeeman relaxation in low-temperature collisions of CaH((2)Σ(+)) and NH((3)Σ(-)) molecules with (3)He atoms converge quickly with respect to the number of total angular momentum states included in the basis set, leading to a dramatic (>10-fold) enhancement in computational efficiency compared to the previously used methods [A. Volpi and J. L. Bohn, Phys. Rev. A 65, 052712 (2002); R. V. Krems and A. Dalgarno, J. Chem. Phys. 120, 2296 (2004)]. Our approach is thus well suited for theoretical studies of strongly anisotropic molecular collisions in the presence of external electromagnetic fields.     

Calculations of fine-structure resolved collisional rate coefficients for the NH(X3Σ-)-He system

We present fine-structure-resolved collisional rate coefficients for the NH(X(3)Σ(-))-He van der Waals complex. The calculations are based on the state-of-the-art potential energy surface [Cybulski et al., J. Chem. Phys. 122, 094307 (2005)]. Close-coupling calculations of the collisional excitation cross sections of the fine-structure levels of NH by He are calculated for total energies up to 3500 cm(-1), which yield, after thermal average, rate coefficients up to 350 K. The fine-structure splitting of rotational levels is taken into account rigorously. The propensity rules between fine-structure levels are reported, and it is found that F-conserving cross sections are much larger than F-changing cross sections, as expected from theoretical considerations. The calculated rate coefficients are compared with available experimental measurements at room temperature and a fairly good agreement is found between experimental and theoretical data. The agreement confirms the relatively good quality of the scattering calculations and also the accuracy of the potential energy surface used in this work. The new set of thermal rate coefficients for this system may be used for improvements in astrophysical and atmospherical modeling.     

Resonances in rotationally inelastic scattering of OH(X2Π) with helium and neon

We present detailed calculations on resonances in rotationally and spin-orbit inelastic scattering of OH (X(2)Π, j = 3/2, F(1), f) radicals with He and Ne atoms. We calculate new ab initio potential energy surfaces for OH-He, and the cross sections derived from these surfaces compare well with the recent crossed beam scattering experiment of Kirste et al. [Phys. Rev. A 82, 042717 (2010)]. We identify both shape and Feshbach resonances in the integral and differential state-to-state scattering cross sections, and we discuss the prospects for experimentally observing scattering resonances using Stark decelerated beams of OH radicals.     

Theoretical investigation of rotationally inelastic collisions of CH2(X) with helium

Following our earlier work on collisions of He with the methylene radical in its excited ?(1)A(1) state [L. Ma, M. H. Alexander, and P. J. Dagdigian, J. Chem. Phys. 134, 154307 (2011)], we investigate here the analogous relaxation of CH(2) in its ground X(3)B(1) electronic state. The molecule is treated as semi-rigid, with fixed bond lengths but a varying bond angle. We use an ab initio potential energy surface (PES) which is averaged over the CH(2) bending angle weighted by the square of the bending wave function. The PES for the interaction of He with CH(2) in the X state is considerably less anisotropic than for interaction with the ? state since the two 2p electrons on the C atom are evenly distributed among the bonding and non-bonding molecular orbitals. We report quantum scattering calculations of state-to-state and total removal cross sections as well as total removal rate constants at room temperature. Because of the less pronounced anisotropy, these cross sections and rate constants are considerably smaller than for collisions of CH(2)(?) with He. Finally, we investigate the dependence of rotational inelasticity on the bending vibrational quantum number.     

Low-energy electron scattering with the purine bases of DNA/RNA using the R-matrix method

R-matrix calculations on electron collisions with the purine bases found in DNA and RNA (i.e., adenine and guanine) are presented. Resonant anion states of these systems are identified by employing different approximation levels of ab initio theoretical methods, such as the static exchange, the static exchange plus polarization, and the close-coupling methods. The results are compared with other available calculations and experiments. All of these ab initio approximations, which we refer to as a scattering "model," give four shape resonances of (2)A' (π) symmetry within the energy range of 10 eV for both molecules. For adenine, the most sophisticated method, the close-coupling model, gives two very narrow (2)A' (σ) symmetry Feshbach-type resonances at energies above 5 eV. Quantitative results for the total elastic and electronic excitation cross sections are also presented.     

Analytical ab initio potential-energy surfaces for the ground and the first singlet excited states of HeH2

Analytical potential-energy surfaces have been constructed for the ground and the first excited states of HeH₂. The functions fit ab initio MRD CI calculations with standard deviations of 0.05 and 0.13 eV for the ground and the excited surface respectively. Classical trajectory calculations for collisions of ⁴Hc with HD(B ¹Σ_u⁺, υ = 3, J = 2) at the temperature T = 297 K yields the electronic quenching cross section σ_Q = 6.5 A² and the vibrational cross section σ_3→2 = 3.8 A². The results are in qualitative agreement with the experimental values of Fink, Akins and Moore.     

The role of rotation in the vibrational relaxation of water by hydrogen molecules

Vibrational relaxation cross sections of the H(2)O(upsilon(2) = 1) bending mode by H(2) molecules are calculated on a recent high-accuracy ab initio potential-energy surface using quasiclassical trajectory calculations. The role of molecular rotation is investigated at a collisional energy of 3500 cm(-1) and it is shown that initial rotational excitation significantly enhances the total (rotationally summed) vibrational relaxation cross sections. A strong and complex dependence on the orientation of the water angular momentum is also observed, suggesting the key role played by the asymmetry of water. Despite the intrinsic limitations of classical mechanics, these exploratory results suggest that quantum approximations based on a complete decoupling of rotation and vibration, such as the widely used vibrational close-coupling (rotational) infinite-order-sudden method, would significantly underestimate rovibrationally inelastic cross sections. We also present some rationale for the absence of dynamical chaos in the scattering process.     

He-ThO(1Σ+) interactions at low temperatures: elastic and inelastic collisions, transport properties, and complex formation in cold 4He gas

We present an ab initio study of cold (4)He + ThO((1)Σ(+)) collisions based on an accurate potential energy surface (PES) evaluated by the coupled cluster method with single, double, and noniterative triple excitations using an extended basis set augmented by bond functions. Variational calculations of rovibrational energy levels show that the (4)He-ThO van der Waals complex has a binding energy of 10.9 cm(-1) in its ground J = 0 rotational state. The calculated energy levels are used to obtain the temperature dependence of the chemical equilibrium constant for the formation of the He-ThO complex. We find that complex formation is thermodynamically favored at temperatures below 1 K and predict the maximum abundance of free ground-state ThO(v = 0, j = 0) molecules between 2 and 3 K. The calculated cross sections for momentum transfer in elastic He + ThO collisions display a rich resonance structure below 5 cm(-1) and decline monotonically above this collision energy. The cross sections for rotational relaxation accompanied by momentum transfer decline abruptly to zero at low collision energies (<0.1 cm(-1)). We find that Stark relaxation in He + ThO collisions can be enhanced by applying an external dc electric field of less than 100 kV∕cm. Finally, we present calculations of thermally averaged diffusion cross sections for ThO in He gas, and find these to be insensitive to small variations of the PES at temperatures above 1 K.     

New ab initio potential energy surface and quantum dynamics of the reaction H(2S) + NH(X3Σ-) → N(4S) + H2

A new global potential energy surface is reported for the ground state ((4)A(")) of the reaction H((2)S) + NH(X(3)Σ(-)) → N((4)S) + H(2) from a set of accurate ab initio data, which were computed using the multi-reference configuration interaction with a basis set of aug-cc-pV5Z. The many-body expansion and neural network methods have been used to construct the new potential energy surface. The topographical features of the new global potential energy surface are presented. The predicted barrier height is lower than previous theoretical estimates and the heat of reaction with zero-point energy is closer to experimental results. The quantum reactive scattering dynamics calculation was carried out over a range of collision energies (0-1.0 eV) on the new potential energy surface. The reaction probabilities, integral cross-section, and rate constants for the title reaction were calculated. The calculated rate constants are in excellent agreement with the available experimental results.     

A time-dependent wave packet quantum scattering study of the reaction HD+ (v = 0 - 3;j0 = 1) + He --> HeH+(HeD+) + D(H)

Time-dependent wave packet quantum scattering (TWQS) calculations are presented for HD(+) (v = 0 - 3;j(0)=1) + He collisions in the center-of-mass collision energy (E(T)) range of 0.0-2.0 eV. The present TWQS approach accounts for Coriolis coupling and uses the ab initio potential energy surface of Palmieri et al. [Mol. Phys. 98, 1839 (2000)]. For a fixed total angular momentum J, the energy dependence of reaction probabilities exhibits quantum resonance structure. The resonances are more pronounced for low J values and for the HeH(+) + D channel than for the HeD(+) + H channel and are particularly prominent near threshold. The quantum effects are no longer discernable in the integral cross sections, which compare closely to quasiclassical trajectory calculations conducted on the same potential energy surface. The integral cross sections also compare well to recent state-selected experimental values over the same reactant and translational energy range. Classical impulsive dynamics and steric arguments can account for the significant isotope effect in favor of the deuteron transfer channel observed for HD(+)(v<3) and low translational energies. At higher reactant energies, angular momentum constraints favor the proton-transfer channel, and isotopic differences in the integral cross sections are no longer significant. The integral cross sections as well as the J dependence of partial cross sections exhibit a significant alignment effect in favor of collisions with the HD(+) rotational angular momentum vector perpendicular to the Jacobi R coordinate. This effect is most pronounced for the proton-transfer channel at low vibrational and translational energies.     

Fully quantum state-resolved inelastic scattering between He and NO(X (2)Pi)

Quantum mechanical close-coupling calculations have been used to obtain fully quantum state-resolved differential cross sections and opacity functions for the rotationally inelastic collisions of NO(X (2)Pi) with He at collision energies of 63 and 147 meV using the most recent ab initio potential energy surfaces of K?os et al. [J. Chem. Phys. 112, 2195 (2000)]. Double peaks observed in the Lambda-doublet resolved differential cross sections are shown to be related to the presence of analogous peaks in the corresponding opacity functions. These structures can be linked directly to a specific expansion term in the potential, and reflect the fact that NO is not quite homonuclear.     

Experimental and theoretical differential cross sections for the N(2D) + H2 reaction

In this paper, we report a combined experimental and theoretical study on the dynamics of the N(2D) + H2 insertion reaction at a collision energy of 15.9 kJ mol(-1). Product angular and velocity distributions have been obtained in crossed beam experiments and simulated by using the results of quantum mechanical (QM) scattering calculations on the accurate ab initio potential energy surface (PES) of Pederson et al. (J. Chem. Phys. 1999, 110, 9091). Since the QM calculations indicate that there is a significant coupling between the product angular and translational energy distributions, such a coupling has been explicitly included in the simulation of the experimental results. The very good agreement between experiment and QM calculations sustains the accuracy of the NH2 ab initio ground state PES. We also take the opportunity to compare the accurate QM differential cross sections with those obtained by two approximate methods, namely, the widely used quasiclassical trajectory calculations and a rigorous statistical method based on the coupled-channel theory.