Strong geometric-phase effects in the hydrogen-exchange reaction at high collision energies: II. Quasiclassical trajectory analysis

Electronic nonadiabatic effects in the state-to-state dynamics of the H + H2 → H2 + H exchange reaction with a vibrationally excited reagent

Out of the many major breakthroughs that the hydrogen-exchange reaction has led to, electronic nonadiabatic effects that are mainly due to the geometric phase has intrigued many. In this work we investigate such effects in the state-to-state dynamics of the H + H2 (v = 3, 4, j = 0) → H2 (v', j') + H reaction with a vibrationally excited reagent at energies corresponding to thermal conditions. The dynamical calculations are performed by a time-dependent quantum mechanical method both on the lower adiabatic potential energy surface (PES) and also using a two-states coupled diabatic theoretical model to explicitly include all the nonadiabatic couplings present in the 1E' ground electronic manifold of the H3 system. The nonadiabatic couplings are considered here up to the quadratic term; however, the effect of the latter on the reaction dynamics is found to be very small. Adiabatic population analysis showed a minimal participation of the upper adiabatic surface even for the vibrationally excited reagent. A strong nonadiabatic effect appears in the state-to-state reaction probabilities and differential cross sections (DCSs). This effect is manifested as "out-of-phase" oscillations in the DCSs between the results of the uncoupled and coupled surface situations. The oscillations persist as a function of both scattering angle and collision energy in both the backward and forward scattering regions. The origins of these oscillations are examined in detail. The oscillations that appear in the forward direction are found to be different from those due to glory scattering, where the latter showed a negligibly small nonadiabatic effect. The nonadiabatic effects are reduced to a large extent when summed over all product quantum states, in addition to the cancellation due to integration over the scattering angle and partial wave summation.

Effect of Reagent Vibration and Rotation on the State-to-State Dynamics of the Hydrogen Exchange Reaction, H + H2 → H2 + H

State-to-state dynamics of the benchmark hydrogen exchange reaction H + H₂ (v = 0-4, j = 0-3) → H₂ (v', j') + H is investigated with the aid of the real wave packet approach of Gray and Balint-Kurti (J. Chem. Phys. 1998, 108, 950-962) and electronic ground BKMP2 potential energy surface of Boothroyd et al. (J. Chem. Phys. 1996, 104, 7139-7152). Initial state-selected and product state-resolved reaction probabilities, integral cross section, and product diatom vibrational and rotational level populations at a few collision energies are reported to elucidate the energy disposal mechanism. State-specific thermal rate constants are also calculated and compared with the available literature results. Coriolis coupling terms of the nuclear Hamiltonian are included, and calculations are parallelized over the helicity quantum number, Ω'. Attempts are made, in particular, to study the effect of reagent vibrational and rotational excitations on the dynamical attributes. It is found that the calculations become computationally expensive with reagent vibrational and rotational excitation. Reagent vibrational excitation is found to enhance the reactivity and has significant impact on the energy disposal to the vibrational and rotational degrees of freedom of the product. The interplay of reagent translational and vibrational energy on the product vibrational distribution unfolds an important aspect of the energy disposal mechanism. The effect of reagent rotation on the state-to-state dynamics is found not to be very significant, and the weak effect turns out to be specific to v'.

Semiclassical initial value representation: From Møller to Miller

The initial value representation (IVR) was proposed five decades ago by Miller [J. Chem. Phys. 53, 3578 (1970)] in order to improve the feasibility and accuracy of semiclassical (SC) scattering calculations. Møller operators, which play a fundamental role in quantum scattering theory, do not appear in his formulation based on action-angle coordinates. These operators were introduced much later by Garashchuk and Light [J. Chem. Phys. 114, 1060 (2001)] in SC-IVR calculations performed in Cartesian coordinates within the Tannor and Weeks [J. Chem. Phys. 98, 3884 (1993)] formulation of quantum scattering theory. Remarkably, Møller operators were found to boost the numerical efficiency of SC-IVR calculations. The purpose of this work is to show within a simple model of light-induced rotational transitions that, in fact, Møller operators were already underlying Miller's pioneering formulation. In line with the results of Garashchuk and Light [J. Chem. Phys. 114, 1060 (2001)], removing the action of these operators in Miller's theory strongly decreases its numerical efficiency.

When classical trajectories get to quantum accuracy: II. The scattering of rotationally excited H2 on Pd(111)

The classical trajectory method in a quantum spirit assigns statistical weights to classical paths on the basis of two semiclassical corrections: Gaussian binning and the adiabaticity correction.

Nonadiabatic Ultracold Quantum Reactive Scattering of Hydrogen with Vibrationally Excited HD(v = 5-9)

The results from electronically non-adiabatic and adiabatic quantum reactive scattering calculations are presented for the H + HD(v = 5-9) → H + HD(v', j') reaction at ultracold collision energies from 10 nK to 60 K. Several experimentally verifiable signatures of the geometric phase are reported in the total and vibrationally and rotationally resolved rate coefficients. Most notable is the predicted 2 orders of magnitude enhancement of the rotationally resolved ultracold rates of odd symmetry relative to those of even symmetry. Prominent shape resonances appear at higher collision energies (100 mK to 20 K), which could be measured experimentally. Significant geometric phase effects are also reported on the resonance energies and lifetimes. In particular, an enhancement (suppression) of the l = 1 (l = 2) shape resonances for HD(v = 5, 6) is predicted for even symmetry relative to those of odd symmetry.

Geometric phase effects on photodissociation dynamics of diatomics

We investigated the effect of the geometric phase (GP) on photodissociation dynamics at a light-induced conical intersection (LICI) through exact quantum dynamical calculations. By taking the one-photon photodissociation of H 2 + ionic molecules as an example, we explored the conditions wherein the LICI associated GP affects dissociation dynamics. We found that GP leads to a phase shift between the angular distributions of GP included and GP excluded photofragments. This effect is more pronounced when the energy of the initial vibrational level is above the energy of the LICI point.

Observation of the geometric phase effect in the H + HD → H2+ D reaction

Theory has established the importance of geometric phase (GP) effects in the adiabatic dynamics of molecular systems with a conical intersection connecting the ground- and excited-state potential energy surfaces, but direct observation of their manifestation in chemical reactions remains a major challenge. Here, we report a high-resolution crossed molecular beams study of the H + HD → H2+ D reaction at a collision energy slightly above the conical intersection. Velocity map ion imaging revealed fast angular oscillations in product quantum state–resolved differential cross sections in the forward scattering direction for H2products at specific rovibrational levels. The experimental results agree with adiabatic quantum dynamical calculations only when the GP effect is included.

Unambiguous Signature of the Berry Phase in Intense Laser Dissociation of Diatomic Molecules

We report strong evidence of Berry phase effects in intense laser dissociation of D₂⁺ molecules, manifested as Aharonov-Bohm-like oscillations in the photofragment angular distribution (PAD). Our calculations show that this interference pattern strongly depends on the parity of the diatom initial rotational state, (-1) ^j. Indeed, the PAD local maxima (minima) observed in one case ( j odd) correspond to local minima (maxima) in the other case ( j even). Using simple topological arguments, we clearly show that such interference conversion is a direct signature of the Berry phase. The sole effect of the latter on the rovibrational wave function is a sign change of the relative phase between two interfering components, which wind in opposite senses around a light-induced conical intersection (LICI). Therefore, encirclement of the LICI leads to constructive ( j odd) or destructive ( j even) self-interference of the initial nuclear wavepacket in the dissociative limit. To corroborate our theoretical findings, we suggest an experiment of strong-field indirect dissociation of D₂⁺ molecules, comparing the PAD of the ortho and para molecular species in directions nearly perpendicular to the laser polarization axis.

Geometric Phase Appears in the Ultracold Hydrogen Exchange Reaction

Quantum reactive scattering calculations for the hydrogen exchange reaction H+H_{2}(v=4,j=0)→H+H_{2}(v^{'}, j^{'}) and its isotopic analogues are reported for ultracold collision energies. Because of the unique properties associated with ultracold collisions, it is shown that the geometric phase effectively controls the reactivity. The rotationally resolved rate coefficients computed with and without the geometric phase are shown to differ by up to 4 orders of magnitude. The effect is also significant in the vibrationally resolved and total rate coefficients. The dynamical origin of the effect is discussed and the large geometric phase effect reported here might be exploited to control the reactivity through the application of external fields or by the selection of a particular nuclear spin state.

The geometric phase controls ultracold chemistry

The geometric phase is shown to control the outcome of an ultracold chemical reaction. The control is a direct consequence of the sign change on the interference term between two scattering pathways (direct and looping), which contribute to the reactive collision process in the presence of a conical intersection (point of degeneracy between two Born-Oppenheimer electronic potential energy surfaces). The unique properties of the ultracold energy regime lead to an effective quantization of the scattering phase shift enabling maximum constructive or destructive interference between the two pathways. By taking the O+OH→H+O2 reaction as an illustrative example, it is shown that inclusion of the geometric phase modifies ultracold reaction rates by nearly two orders of magnitude. Interesting experimental control possibilities include the application of external electric and magnetic fields that might be used to exploit the geometric phase effect reported here and experimentally switch on or off the reactivity.

Correlated electron-nuclear dynamics: exact factorization of the molecular wavefunction

It was recently shown [A. Abedi, N. T. Maitra, and E. K. U. Gross, Phys. Rev. Lett. 105, 123002 (2010)] that the complete wavefunction for a system of electrons and nuclei evolving in a time-dependent external potential can be exactly factorized into an electronic wavefunction and a nuclear wavefunction. The concepts of an exact time-dependent potential energy surface (TDPES) and exact time-dependent vector potential emerge naturally from the formalism. Here, we present a detailed description of the formalism, including a full derivation of the equations that the electronic and nuclear wavefunctions satisfy. We demonstrate the relationship of this exact factorization to the traditional Born-Oppenheimer expansion. A one-dimensional model of the H(2)(+) molecule in a laser field shows the usefulness of the exact TDPES in interpreting coupled electron-nuclear dynamics: we show how features of its structure indicate the mechanism of dissociation. We compare the exact TDPES with potential energy surfaces from the time-dependent Hartree-approach, and also compare traditional Ehrenfest dynamics with Ehrenfest dynamics on the exact TDPES.