Numerical study of the geometric phase in the H+H2 reaction

Strong Angular Oscillation of Rotationally Resolved Differential Cross Section in the H + HD → H2 + D Reaction at the Collision Energy of 2.07 eV: Evidence of Geometric Phase Effects

Geometric phase (GP) effects in chemical reactions are subtle quantum phenomena that are challenging to identify. In this work, we report a joint experimental and theoretical study of the H + HD → H2 + D reaction at a collision energy of 2.07 eV, which is far below the energy of the conical intersection of 2.53 eV. The rotationally state-resolved differential cross sections were measured by a crossed-beam experiment with the scheme of D-atom Rydberg tagging time-of-flight detection. Experimental angular distributions of three rotational states of H2 products exhibit notable variation near the backward scattering direction. Time-dependent quantum mechanics calculations (TDQMs) were carried out at the same collision energy, with and without the inclusion of GP. The experimental angular distributions are in good agreement with TDQM results with the inclusion of GP but do not agree with TDQM results without the inclusion of GP. This work demonstrates the existence of GP effects at energy far below the conical intersection.

Observation of geometric phase effect through backward angular oscillations in the H + HD → H2 + D reaction

Quantum interference between reaction pathways around a conical intersection (CI) is an ultrasensitive probe of detailed chemical reaction dynamics. Yet, for the hydrogen exchange reaction, the difference between contributions of the two reaction pathways increases substantially as the energy decreases, making the experimental observation of interference features at low energy exceedingly challenging. We report in this paper a combined experimental and theoretical study on the H + HD → H2 + D reaction at the collision energy of 1.72 eV. Although the roaming insertion pathway constitutes only a small fraction (0.088%) of the overall contribution, angular oscillatory patterns arising from the interference of reaction pathways were clearly observed in the backward scattering direction, providing direct evidence of the geometric phase effect at an energy of 0.81 eV below the CI. Furthermore, theoretical analysis reveals that the backward interference patterns are mainly contributed by two distinct groups of partial waves (J ~ 10 and J ~ 19). The well-separated partial waves and the geometric phase collectively influence the quantum reaction dynamics.

Geometric Phase Effects in the H + H2+ Reaction on the Lowest Triplet Ground State

To our knowledge, this is the first time geometric phase (GP) effects in the H + H2+ reaction on its lowest triplet ground state with collision energies lower than 1.85 eV have been studied using the quantum wave packet method and vector potential approach. We obtained the total reaction probabilities including and not including GP (NGP) effects for J ≤ 4. Visible GP effects could be seen at the lower energy regime but are tiny at the higher one. Moreover, they are more obvious in the product rovibrational state-resolved reaction probabilities, and the relative resonance magnitudes between GP and NGP results change with product rotational state values alternatively. The main reasons are the interferences between the one- and two-transition-state (1-TS and 2-TS) reaction paths, in that at the lower energy regime, the reaction probabilities from the 2-TS pathway show peaks of comparable probabilities compared with that of the 1-TS pathway. In addition, the "out-of-phase" trend observed in the H + H2 reaction does not exist rigorously in this system. Importantly, the visible GP effects exist in this H3+ system, which makes it a very useful candidate reaction for nonadiabatic investigations in both theory and experiment.

Direct observation of geometric-phase interference in dynamics around a conical intersection

Conical intersections are ubiquitous in chemistry and physics, often governing processes such as light harvesting, vision, photocatalysis and chemical reactivity. They act as funnels between electronic states of molecules, allowing rapid and efficient relaxation during chemical dynamics. In addition, when a reaction path encircles a conical intersection, the molecular wavefunction experiences a geometric phase, which can affect the outcome of the reaction through quantum-mechanical interference. Past experiments have measured indirect signatures of geometric phases in scattering patterns and spectroscopic observables, but there has been no direct observation of the underlying wavepacket interference. Here we experimentally observe geometric-phase interference in the dynamics of a wavepacket travelling around an engineered conical intersection in a programmable trapped-ion quantum simulator. To achieve this, we develop a technique to reconstruct the two-dimensional wavepacket densities of a trapped ion. Experiments agree with the theoretical model, demonstrating the ability of analogue quantum simulators-such as those realized using trapped ions-to accurately describe nuclear quantum effects.

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.

Topological Metamaterials

The topological properties of an object, associated with an integer called the topological invariant, are global features that cannot change continuously but only through abrupt variations, hence granting them intrinsic robustness. Engineered metamaterials (MMs) can be tailored to support highly nontrivial topological properties of their band structure, relative to their electronic, electromagnetic, acoustic and mechanical response, representing one of the major breakthroughs in physics over the past decade. Here, we review the foundations and the latest advances of topological photonic and phononic MMs, whose nontrivial wave interactions have become of great interest to a broad range of science disciplines, such as classical and quantum chemistry. We first introduce the basic concepts, including the notion of topological charge and geometric phase. We then discuss the topology of natural electronic materials, before reviewing their photonic/phononic topological MM analogues, including 2D topological MMs with and without time-reversal symmetry, Floquet topological insulators, 3D, higher-order, non-Hermitian and nonlinear topological MMs. We also discuss the topological aspects of scattering anomalies, chemical reactions and polaritons. This work aims at connecting the recent advances of topological concepts throughout a broad range of scientific areas and it highlights opportunities offered by topological MMs for the chemistry community and beyond.

Geometric Phase Effect in Attosecond Stimulated X-ray Raman Spectroscopy

Conical intersections (CIs) are diabolical points in the potential energy surfaces generally caused by point-wise degeneracy of different electronic states, and give rise to the geometric phases (GPs) of molecular wave functions. Here we theoretically propose and demonstrate that the transient redistribution of ultrafast electronic coherence in attosecond Raman signal (TRUECARS) spectroscopy is capable of detecting the GP effect in excited state molecules by applying two probe pulses including an attosecond and a femtosecond X-ray pulse. The mechanism is based on a set of symmetry selection rules in the presence of nontrivial GPs. The model of this work can be realized for probing the geometric phase effect in the excited state dynamics of complex molecules with appropriate symmetries, using attosecond light sources such as free-electron X-ray lasers.

Charge Transfer Processes for H + H2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces

The dynamics of the H + H₂⁺ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (2¹A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H₂⁺ (²Σ_g⁺) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.

Time-dependent quantum mechanical wave packet dynamics

Starting from a model study of the collinear (H, H₂) exchange reaction in 1959, the time-dependent quantum mechanical wave packet (TDQMWP) method has come a long way in dealing with systems as large as Cl + CH₄. The fast Fourier transform method for evaluating the second order spatial derivative of the wave function and split-operator method or Chebyshev polynomial expansion for determining the time evolution of the wave function for the system have made the approach highly accurate from a practical point of view. The TDQMWP methodology has been able to predict state-to-state differential and integral reaction cross sections accurately, in agreement with available experimental results for three dimensional (H, H₂) collisions, and identify reactive scattering resonances too. It has become a practical computational tool in predicting the observables for many A + BC exchange reactions in three dimensions and a number of larger systems. It is equally amenable to determining the bound and quasi-bound states for a variety of molecular systems. Just as it is able to deal with dissociative processes (without involving basis set expansion), it is able to deal with multi-mode nonadiabatic dynamics in multiple electronic states with equal ease. We present an overview of the method and its strength and limitations, citing examples largely from our own research groups.

Quantum interference in H + HD → H2 + D between direct abstraction and roaming insertion pathways

Understanding quantum interferences is essential to the study of chemical reaction dynamics. Here, we provide an interesting case of quantum interference between two topologically distinct pathways in the H + HD → H2 + D reaction in the collision energy range between 1.94 and 2.21 eV, manifested as oscillations in the energy dependence of the differential cross section for the H2 (v′ = 2, j′ = 3) product (where v′ is the vibrational quantum number and j′ is the rotational quantum number) in the backward scattering direction. The notable oscillation patterns observed are attributed to the strong quantum interference between the direct abstraction pathway and an unusual roaming insertion pathway. More interestingly, the observed interference pattern also provides a sensitive probe of the geometric phase effect at an energy far below the conical intersection in this reaction, which resembles the Aharonov–Bohm effect in physics, clearly demonstrating the quantum nature of chemical reactivity.

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.

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.

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.

Direct evaluation of the temperature dependence of the rate constant based on the quantum instanton approximation

A general method for the direct evaluation of the temperature dependence of the quantum-mechanical reaction rate constant in many-dimensional systems is described. The method is based on the quantum instanton approximation for the rate constant, thermodynamic integration with respect to the inverse temperature, and the path integral Monte Carlo evaluation. It can describe deviations from the Arrhenius law due to the coupling of rotations and vibrations, zero-point energy, tunneling, corner-cutting, and other nuclear quantum effects. The method is tested on the Eckart barrier and the full-dimensional H+H(2)-->H(2)+H reaction. In the temperature range from 300 to 1500 K, the error of the present method remains within 13% despite the very large deviations from the Arrhenius law. The direct approach makes the calculations much more efficient, and the efficiency is increased even further (by up to two orders of magnitude in the studied reactions) by using optimal estimators for reactant and transition state thermal energies. Which of the estimators is optimal, however, depends on the system and the strength of constraint in a constrained simulation.

Theoretical studies on bimolecular reaction dynamics

This perspective discusses progress in the theory of bimolecular reaction dynamics in the gas phase. The examples selected show that definitive quantum dynamical computations are providing insights into the detailed mechanisms of chemical reactions.

Effect of the geometric phase on the dynamics of the hydrogen-exchange reaction

A recent puzzle in nonadiabatic quantum dynamics is that geometric phase (GP) effects are present in the state-to-state opacity functions of the hydrogen-exchange reaction, but cancel out in the state-to-state integral cross sections (ICSs). Here the authors explain this result by using topology to separate the scattering amplitudes into contributions from Feynman paths that loop in opposite senses around the conical intersection. The clockwise-looping paths pass over one transition state (1-TS) and scatter into positive deflection angles; the counterclockwise-looping paths pass over two transition states (2-TS) and scatter into negative deflection angles. The interference between the 1-TS and 2-TS paths thus integrates to a very small value, which cancels the GP effects in the ICS. Quasiclassical trajectory (QCT) calculations reproduce the scattering of the 1-TS and 2-TS paths into positive and negative deflection angles and show that the 2-TS paths describe a direct insertion mechanism. The inserting atom follows a highly constrained "S-bend" path, which allows it to avoid both the other atoms and the conical intersection and forces the product diatom to scatter into high rotational states. By contrast, the quantum 2-TS paths scatter into a mainly statistical distribution of rotational states, so that the quantum 2-TS total ICS is roughly twice the QCT ICS at 2.3 eV total energy. This suggests that the S-bend constraint is relaxed by tunneling in the quantum system. These findings on H+H(2) suggest that similar cancellations or reductions in GP effects are likely in many other reactions.

The photoelectron spectrum of elusive cyclic-N3 and characterization of the potential energy surface and vibrational states of the ion

A potential energy surface is constructed for the ground X (1)A(1) electronic state of cyclic-N(3) (+) based on three-dimensional spline interpolation of ab initio points. The vibrational states of this molecular ion are calculated in the range up to 14 500 cm(-1) using hyperspherical coordinates and the coupled-channel (sector-adiabatic) approach. All the vibrational states are analyzed and assigned. The Franck-Condon overlaps of these states with the vibrational states of the neutral are calculated to predict the photoelectron spectrum of cyclic-N(3). Peak intensities are governed by the nodal structure of the vibrational wave functions and reflect the large geometric phase effect predicted for cyclic-N(3). Experimental validation may shed light on the existence of this elusive molecule and confirm the magnitude of the geometric phase effect.

Incorporating the Geometric Phase Effect in Triatomic and Tetraatomic Hyperspherical Harmonics

Hyperspherical harmonics in the democratic row-orthonormal hyperspherical coordinates are very appropriate basis sets for performing reactive scattering calculations for triatomic and tetraatomic systems. The mathematical conditions for incorporating the geometric phase effect in these harmonics are given. These conditions are implemented for triatomic systems, and their explicit analytical expressions in terms of Jacobi polynomials, in both the absence and presence of the geometric phase effect, are given.

Theoretical study of geometric phase effects in the hydrogen-exchange reaction

The crossing of two electronic potential surfaces (a conical intersection) should result in geometric phase effects even for molecular processes confined to the lower surface. However, recent quantum simulations of the hydrogen exchange reaction (H + H2 --> H2 + H) have predicted a cancellation in such effects when product distributions are integrated over all scattering angles. We used a simple topological argument to extract reaction paths with different senses from a nuclear wave function that encircles a conical intersection. In the hydrogen-exchange reaction, these senses correspond to paths that cross one or two transition states. These two sets of paths scatter their products into different regions of space, which causes the cancellation in geometric phase effects. The analysis should generalize to other direct reactions.

Geometric phase effects in the H+H2 reaction: Quantum wave-packet calculations of integral and differential cross sections

We report quantum wave-packet calculations on the H+H(2) reaction, aimed at resolving the controversy over whether geometric phase (GP) effects can be observed in this reaction. Two sets of calculations are reported of the state-to-state reaction probabilities, and integral and differential cross sections (ICSs and DCSs). One set includes the GP using the vector potential approach of Mead and Truhlar; the other set neglects the phase. We obtain unequivocal agreement with recent results of Kendrick [J. Phys. Chem. A 107, 6739 (2003)], predicting GP effects in the state-to-state reaction probabilities, which cancel exactly on summing the partial waves to yield the ICS. Our results therefore contradict those of Kuppermann and Wu [Chem. Phys. Lett. 349 537 (2001)], which predicted pronounced GP effects in the cross sections. We also agree with Kendrick in predicting that there are no significant GP effects in the full DCS at energies below 1.8 eV, and in the partial (0<or=J<or=10) DCS at energies above this. However, we find that in the full DCS above 1.8 eV (which was not reported by Kendrick), there are GP effects, which may be experimentally measurable.

Cyclic-N3. II. Significant geometric phase effects in the vibrational spectra

An accurate theoretical prediction of the vibrational spectra for a pure nitrogen ring (cyclic-N(3)) molecule is obtained up to the energy of the (2)A(2)/(2)B(1) conical intersection. A coupled-channel approach using the hyperspherical coordinates and the recently published ab initio potential energy surface [D. Babikov, P. Zhang, and K. Morokuma, J. Chem. Phys. 121, 6743 (2004)] is employed. Two independent sets of calculations are reported: In the first set, the standard Born-Oppenheimer approximation is used and the geometric phase effects are totally neglected. In the second set, the generalized Born-Oppenhimer approximation is used and the geometric phase effects due to the D(3h) conical intersection are accurately treated. All vibrational states are analyzed and assigned in terms of the normal vibration mode quantum numbers. The magnitude of the geometric phase effect is determined for each state. One important finding is an unusually large magnitude of the geometric phase effects in the cyclic-N(3): it is approximately 100 cm(-1) for the low-lying vibrational states and exceeds 600 cm(-1) for several upper states. On average, this is almost two orders of magnitude larger than in the previously reported studies. This unique example suggests a favorable path to experimental validation.