Modeling of adiabatic and diabatic potential energy surfaces of Cl(2P)⋯H2(1∑g+) prereactive complex from ab initio calculations

Infrared-active spin-orbit transitions of halogen atom dopants in solid parahydrogen: the role of trapping site geometry

We present theoretical calculations of the (2)P(1/2) ← (2)P(3/2) spin-orbit transition of Cl dopants embedded as substitutional impurities in solid parahydrogen (pH2) matrices. In the lower-energy (2)P(3/2) spin-orbit level, the Cl atom's electron density distribution is anisotropic, and slightly distorts the geometry of the atom's trapping site. This distortion leads to a blue shift in the spin-orbit transition energy; the blue shift is enhanced when we account for the large-amplitude zero point motions of the pH2 molecules surrounding the Cl dopant. We also show that the intensity of the transition depends on the geometry of the trapping site. In the gas phase, the (2)P(1/2) ← (2)P(3/2) atomic transition is electric dipole forbidden. However, when the Cl atom resides in trapping sites that mimic the hexagonal close packed morphology of pure solid pH2, the transition becomes electric dipole allowed through interaction-induced transition dipole moments. These transition dipole moments originate in the anisotropic electron density distribution of the lower-energy (2)P(3/2) spin-orbit level.

Spin-orbit quenching of Cl(2P(1∕2)) by H2

We report fully-quantum, time-independent, scattering calculations for the spin-orbit quenching of Cl((2)P(1∕2)) by H(2) molecules at low and moderate temperature. Our calculations take into account chemical reaction channels. Cross sections are calculated for total energies up to 5000 cm(-1) which are used to determine, by thermal averaging, state-to-state rate coefficients at temperatures ranging from 50 to 500 K. Spin-orbit relaxation of chlorine is dominated by collisions with H(2) in the rotationally excited states j = 2 and j = 3. In the former case the near-resonant energy transfer is the primary relaxation mechanism. The inclusion of the reactive channel could lead differences compared to pure inelastic calculations. Good agreement is obtained with experimental relaxation measurements at room temperature.

Dynamical regimes on the Cl + H2 collisions: inelastic rainbow scattering

While Cl + H(2) reactive collisions have been a subject of numerous experimental and theoretical studies, inelastic collisions leading to rotational energy transfer and/or vibrational excitation have been largely ignored. In this work, extensive quantum mechanical calculations covering the 0.5-1.5 eV total energy range and various initial rovibrational states have been carried out and used to perform a joint study of inelastic and reactive Cl + H(2) collisions. Quasiclassical trajectories calculations complement the quantum mechanical results. The analysis of the inelastic transition probabilities has revealed the existence of two distinct dynamical regimes that correlate with low and high impact parameters, b, and are neatly separated by glory scattering. It has been found that while high-b collisions are mainly responsible for |Δj| = 2 transitions which dominate the inelastic scattering, they are very inefficient in promoting higher |Δj| transitions. The effectiveness of this type of collision also drops with rotational excitation of H(2). In contrast, reactive scattering, that competes with |Δj| > 2 inelastic transitions, is exclusively caused by low-b collisions, and it is greatly favored when the reactants get rotationally excited. Previous studies focusing on the reactivity of the Cl + H(2) system established that the van der Waals well located in the entrance channel play a key role in determining the mechanism of the collisions. Our results prove this to be also a case for inelastic processes, where the origin of the double dynamical regime can be traced back to the influence exerted by this well that shapes the topology of the entrance channel of the Cl-H(2) system.

Spin-orbit relaxation of Cl(P1∕22) and F(P1∕22) in a gas of H2

The authors present quantum scattering calculations of rate coefficients for the spin-orbit relaxation of F(2P1/2) atoms in a gas of H2 molecules and Cl(2P1/2) atoms in a gas of H2 and D2 molecules. Their calculation of the thermally averaged rate coefficient for the electronic relaxation of chlorine in H2 agrees very well with an experimental measurement at room temperature. It is found that the spin-orbit relaxation of chlorine atoms in collisions with hydrogen molecules in the rotationally excited state j=2 is dominated by the near-resonant electronic-to-rotational energy transfer accompanied by rotational excitation of the molecules. The rate of the spin-orbit relaxation in collisions with D2 molecules increases to a great extent with the rotational excitation of the molecules. They have found that the H2/D2 isotope effect in the relaxation of Cl(2P1/2) is very sensitive to temperature due to the significant role of molecular rotations in the nonadiabatic transitions. Their calculation yields a rate ratio of 10 for the electronic relaxation in H2 and D2 at room temperature, in qualitative agreement with the experimental measurement of the isotope ratio of about 5. The isotope effect becomes less significant at higher temperatures.

Interactions of 2P Atoms with Closed-Shell Diatomic Molecules: Alternative Diabatic Representations for the Electronic Anisotropy

The matrices of electrostatic and spin-orbit Hamiltonians for the system of a 2P atom interacting with a closed shell diatomic molecule in uncoupled, coupled, and complex-valued representations for electronic diabatic basis functions are rederived, and the unitary transformations connecting them are given explicitly. The links to previous derivations are established and existing inconsistencies are identified and eliminated. It is proven that the block-diagonalization of a 6 x 6 matrix of the electronic Hamiltonian is a result of using the basis functions with well-defined properties with respect to time reversal. Consideration of time-reversal symmetry also enforces phase consistency relevant for applications to multisurface reactive scattering and photodetachment spectroscopy calculations, as well as for perspective studies of inelastic effects in cold and ultracold environments. These and further developments are briefly sketched.

Interactions in open-shell clusters: ab initio study of pre-reactive complex O(3P) + HCl

Van der Waals interactions between the ground-state triplet O(3P) atom and the closed-shell HCl molecule are investigated in the pre-reactive region. Three adiabatic (two of A'' symmetry and one of A' symmetry) and four non-relativistic diabatic potential energy surfaces are obtained by combining a restricted open-shell coupled cluster approach with the multireference configuration interaction method. The lower A'' adiabatic potential surface has a single minimum (D(e) = 589 cm(-1)) for a linear O...HCl configuration. The upper A'' potential has a weak (D(e) = 65 cm(-1)) minimum for a linear HCl...O configuration. The A' adiabatic potential has a weak (124 cm(-1)) T-shaped minimum. Adiabatic potentials intersect once in the O...HCl linear configuration and twice in the linear HCl...O geometry. The role of electrostatic interactions in shaping these potentials is discussed. The effects of spin-orbit coupling on this interaction are also investigated assuming a constant value of the SO parameter.