## Abstract

It has long been known that there is a conical intersection (CI) between the ground and first excited electronic state in the H_{3} system. Its associated geometric phase (GP) effect has been theoretically predicted to exist below the CI since a long time. However, the experimental evidence has not been established yet and its dynamical origin is waiting to be elucidated. Here we report a combined crossed molecular beam and quantum reactive scattering dynamics study of the H+HD → H_{2}+D reaction at 2.28 eV, which is well below the CI. The GP effect is clearly identified by the observation of distinct oscillations in the differential cross section around the forward direction. Quantum dynamics theory reveals that the GP effect arises from the phase alteration of a small part of the wave function, which corresponds to an unusual roaming-like abstraction pathway, as revealed by quasi-classical trajectory calculations.

## Introduction

The geometric phase (GP) is an important concept in adiabatic theory. It represents the phase acquired when a quantum state is adiabatically transported around a closed circuit in the parameter space, while remaining in the same Hilbert space. The concept of the GP has given rise to various observable effects, and has now become a central topic across various research fields ranging from chemistry to condensed matter physics^{1,2}.

A classical example of GP in condensed phase physics is the anomalous Hall effect^{3}. In a molecule, the presence of a conical intersection (CI) leads to the appearance of the GP and enables rapid conversion of the excess electronic energy into nuclear motion. The GP induces a sign change of the adiabatic electronic wave functions along a closed path of the nuclear configuration space encircling the CI^{4,5,6,7}. For a molecular system with a CI, the theoretical treatment in the adiabatic picture requires the inclusion of the GP to make sure that at each nuclear geometry, the total wave function can be single-valued, although the nuclear wave function is localized away from the region of the CI. The GP effect in the spectra of molecular vibrational bound states has been well studied with Jahn–Teller molecules^{8,9,10}. Recently, the significant role of the GP in photodissociation process of phenol has been reported^{7,11,12}. In larger molecules, the increased number of degrees of freedom and electronic states leads to the general presence of CIs in various molecular photoinduced chemical dynamics processes involving electronically excited states, where CIs play a very important role^{13}.

H+H_{2} is one of the most important chemical reactions since it can be accurately described with current quantum theoretical methods. It exhibits a typical *D*_{3h} CI between the ground and first excited electronic states at equilateral triangle geometries, corresponding to the total energy of about 2.75 eV. Using the ground adiabatic potential energy surface (PES), product state-resolved integral and differential cross sections (DCS) for the H+H_{2} reaction and its isotopologues without consideration of the GP were first reported in 1988^{14,15}. Later, it was found that the adiabatic quantum dynamics predictions without the GP (i.e. NGP) agree with experiments at some collision energies quite well, suggesting that the GP may not be important in this reaction and its isotopologues^{16,17,18,19,20,21}. Mead and Truhlar found that the GP would only influence observables when the nuclear wavefunction encircled the CI, and the GP effect could be involved by the introduction of a vector potential^{22}.

In the last decades, a series of quantum mechanical theoretical studies on the H+H_{2} reactive scattering were carried out by Kendrick and by Althorpe and coworkers. Their results demonstrated that the GP effect should be negligible at total energies below 1.6 eV^{23,24,25,26}, and only becomes significant when the total energy is above 3.5 eV. Quantum wave-packet results by Althorpe and coworkers also suggested that a shift of the fast oscillatory angular distributions of the sideways scattered products could serve as a clear signature to reveal the effect of GP^{19,27,28}. Recently, Kendrick and his coworkers found that the GP effectively controls the reactivity of the H+H_{2} (*v* = 4, *j* = 0) and OH+O → O_{2}+H reactions in ultracold conditions^{29,30}.

In the past several decades, high-resolution crossed molecular beam (CMB) studies have been carried out on many important elementary reactions using the Rydberg state H-atom tagging method^{31,32}, including the H+HD and H+D_{2} reactions at a set of different collision energies^{16,18,20,33}. No fast angular oscillations in DCS were observed, and this is most probably due to the limited experimental angular resolution. More recently, the PHOTOLOC technique has been employed in the study of the GP effect in the H+H_{2} reaction and its related isotopologues, but no clear evidence of the GP was found^{34,35}.

Very recently, with a high-resolution time-sliced velocity map imaging (VMI) – CMB apparatus, we detected D products of the H+HD reactions by using the near threshold ionization scheme^{21,36}. The product ro-vibrational state-resolved DCS with fast oscillatory structures in the forward angular distributions have been observed at the collision energies of 1.35 and 2.77 eV^{21,36}. At 1.35 eV, which is far below the CI, the reaction appears to take place via a simple direct abstraction mechanism, where the GP plays a negligible role. At 2.77 eV (about 2.99 eV in total energy relative to the equilibrium energy of the H_{2} molecule), which is 0.24 eV above the CI, the measured fast angular oscillations of the forward scattered products in the product ro-vibrational state-resolved DCS could only be reproduced by quantum mechanical calculations on the lowest adiabatic PES with GP included, clearly demonstrating the existence of the GP effect in the reaction.

In previous work^{36}, quantum dynamics calculations using accurate diabatic coupled PESs had to be carried out to exclude the role of the adiabatic electronically excited state since the studied energy was above the CI. This follows an intriguing question: at an energy well below the CI, where the role of the adiabatic electronically excited state could be simply excluded, will there be observable effects due to the GP?

Here we report a high-resolution CMB study on the H+HD → H_{2}+D reaction at the collision energy of 2.28 eV, corresponding to 2.50 eV in total energy, or 0.25 eV below the CI. In addition, we carry out accurate adiabatic quantum mechanical time-dependent wave packet calculations with and without considering the GP effect in hyperspherical coordinates to investigate the role of GP on this reaction^{37,38,39}. The GP effect is analyzed by calculating the time-independent scattering wave functions at 2.28 eV in hyperspherical coordinates. At the same time, the dynamical mechanism determining the GP effect is revealed using quasi-classical trajectory (QCT) calculations.

## Results

### High-resolution experimental image

The H+HD→H_{2}+D reactive scattering experiment is carried out using a CMB apparatus with a VMI detector at the collision energy of 2.28 eV. A more detailed experimental procedure is presented in the Supplementary Note 1. The experimental image of the D product (Fig. 1) exhibits a series of well-resolved rings coming from different ro-vibrational states of the H_{2} product, which are clearly assignable. Fine oscillatory features in the forward angular distributions are successfully captured, similar to those observed at 1.35 eV and 2.77 eV^{21,36}. There are two very pronounced rings labeled by yellow arrows in Fig. 1. Each of them corresponds to two ro-vibrational states with very small energy differences. The experimental angular distributions in the forward scattering direction are presented in Fig. 2 for the H_{2} product at the (*v*′ = 0, *j*′ = 11 & *v*′ = 1, *j*′ = 7) and (*v*′ = 1, *j*′ = 9 & *v*′ = 2, *j*′ = 3) levels. Here “&” is used to connect the two ro-vibrational states with very small energy differences. For the corresponding rings, the angular distributions are derived by extracting the signals at various scattering angles (refer to Supplementary Note 2).

### GP identification via quantum dynamics calculations

In order to understand the reaction dynamics of the H+H_{2} reaction and its isotopologues, recently we have developed a state-to-state time-dependent quantum reactive scattering theory in hyperspherical coordinates with and without the inclusion of the GP as a vector potential, similar to that in the Jacobi coordinates^{25,36}. In hyperspherical coordinates, it is straightforward to express the vector potential for the CI of *D*_{3h}, which is the function of only the two hyperangles *θ* and *χ*, thus it is convenient to analyze the GP effect in the reaction with wave functions in hyperspherical coordinates.

To justify whether there is a detectable GP effect at this energy for the H+HD reaction, the adiabatic quantum dynamics calculations are first carried out on the accurate adiabatic BKMP2 PES without including the GP effect (NGP) (Fig. 2a, b). The DCS from the NGP calculations clearly deviated from the experimental results. Especially for the plot in Fig. 2b, the slight oscillations in the 10°–30° angular range of the experimental angular distribution for the H_{2} (*v*′ = 1, *j*′ = 9 & *v*′ = 2, *j*′ = 3) product states are nearly completely out-of-phase with the theoretical NGP DCS, with the peaks in the calculated DCS located at the valley positions of the experimental one.

Time-dependent adiabatic quantum dynamics calculations are then carried out for the H+HD → H_{2}+D reaction on the BKMP2 PES. A vector potential is introduced to include the effect of GP. The calculated DCS with the GP included are shown in Fig. 2c, d. The theoretical GP DCS is in excellent agreement with the experimental results, with the theoretical angular oscillations in the DCS right in phase (and amplitude) with the experimental results, while the NGP results show clear differences with the experimental one. This suggests that, in the adiabatic picture, the GP effect can be unambiguously observed for the benchmark reaction at an energy well below the CI.

## Discussion

At the current relative low collision energy, it is impressive that the GP effect on the H+HD reaction can be identified so clearly around the forward direction in the DCS. The GP only bring in a phase change to the wave functions. The topological argument^{25} therefore infers that the reaction dynamics will not be influenced by the GP if the reaction only take place via a single pathway. As a result, a second reaction pathway, besides the normal abstraction one, is expected for the H+HD → H_{2}+D reaction at 2.28 eV, similar to that for the H+H_{2} reaction. The GP can then alter the interference pattern between these two pathways thus leading to variations in the pattern of the DCS of the products.

In order to investigate the mechanisms accounting for the GP effect in the H+HD → H_{2}+D reaction, we have calculated the time-independent scattering wave function as

with *E*_{c} = 2.28 eV and *E*_{0} as the zero point energy. At *ρ* ~ 2.5 bohr, where the CI has the lowest energy, the GP is expected to play the most significant role. Thus, a closer examination of the wave function at *ρ* ~ 2.5 bohr would be helpful for understanding the mechanisms.

The wave functions at *ρ* ~ 2.5 bohr as a function of *x* = sin(*θ*)cos(2*χ*) and *y* = sin(*θ*)sin(2*χ*) of the NGP and GP calculation are presented in Fig. 3a, b in logarithmic scale (base 10). It is not surprising that the majority part of these two wave functions, which are away from and above the CI in the figure, are almost identical. These regions correspond to the normal abstraction process from the reactant H´ + HD to products H_{2} + D and H´D + H, which experiences only one transition state. We name it as Path 1 (as shown in Fig. 3a).

However, in the region just below the CI, there is a significant difference between these two wave functions, although the amplitudes of the wave functions are quite small. In that region, the GP wave function acquires a node structure that does not appear in the NGP wave function. This part of the wave function should correspond to the reaction path which experiences two transition states^{23}. We name it as Path 2 (as shown in Fig. 3a).

The phases of the wave functions are further calculated and presented in Fig. 3c, d for the NGP and GP calculations, respectively. The line, which starts from the CI and extends up to the right, results from the boundary condition of double values for the nuclear wave function introduced by the vector potential. The phases of the NGP and GP wave functions in most regions are very similar. However, in the area to the right of the node of the wave function below the CI and below the double values boundary condition line, which is mainly located in the (x > 0, y < 0) quadrant in Fig. 3c, d, the phase of the GP wave function differs by about 180° from that of the NGP wave function. This is the physical origin of the difference between the GP and NGP DCS. This additional node in the GP wavefunction in the chemical reaction is very similar to the GP wavefunction in a photodissociation process^{7,12}.

Althorpe and coworkers also extracted the contributions of the two classes of reaction pathways with an approach based on topological arguments^{25,26,27}. In their studies, the nuclear wave functions from Path 1 and Path 2 can be derived by *ψ*_{1} = \((\psi _{NGP} + \psi _{GP})/\sqrt 2\) and *ψ*_{2} =\((\psi _{NGP} - \psi _{GP})/\sqrt 2\), respectively. Subscripts “NGP” and “GP” are used to denote the wave functions without and with the introduction of the GP effect (similarly hereinafter), respectively. The scattering amplitudes from these two pathways can be expressed as,

respectively. The DCSs for the individual paths are given by the square moduli of *f*_{1}(θ) and *f*_{2}(θ), |*f*_{1}(θ)|^{2} and |*f*_{2}(θ)|^{2}. The total DCS for the reaction can be described as,

whereas the interference between Path 1 and Path 2 originates from the last two crossing terms. Since the GP alters the relative phases between the two pathways, a difference in the DCS appears. With this approach, we calculate the total integral cross section (ICS) for the reaction via Path 1 and Path 2 at 2.28 eV. The reaction occurring through Path 2 amounts to only about 0.35% of the total product, which is much smaller than that at 2.77 eV (~2.3% at 2.77 eV)^{36}. Thus the GP effect in the DCS at 2.28 eV is weaker and much more difficult to observe experimentally. With increasing collision energies, a large portion of the wave function encircles the CI in the clockwise rotation direction from the reactant H + HD channel to the product H_{2} + D channel, thus, the DCS of the reaction manifests a stronger GP effect.

The DCS for H_{2} (*v*′ *=* 0, *j*′ = 11), (*v*′ *=* 1, *j*′ = 7), (*v*′ *=* 1, *j*′ = 9), and (*v*′ *=* 2, *j*′ = 3) from Path 1 and Path 2 at 2.28 eV are presented in Fig. 4. The results show that there are very different angular distributions for the two reaction paths. Path 1 leads to an angular distribution with the products predominantly sideways scattered compared to the relatively small forward direction, whereas Path 2 for product H_{2} (*v*′ *=* 0, *j*′ = 11) leads to oscillatory scattering evenly in the whole scattering angular range, for products H_{2} (*v*′ *=* 1, *j*′ = 7 & 9) to significant both forward and backward scattering, and for product H_{2} (*v*′ *=* 2, *j*′ = 3) to dominant backward scattering. The DCS from Path 2 is much smaller, especially for the product H_{2} (*v*′ *=* 0, *j*′ = 11), less than 1% of that from Path 1. Consequently, the oscillations around the forward scattering direction observed in Fig. 2 are weaker than those observed for product H_{2} (*v*′ *=* 0, *j*′ = 7) at the collision energy of 2.77 eV, where the DCS from the two paths are of comparable amplitude.

To shed light on the reaction mechanism about Path 2, QCT calculations are carried out using the adiabatic BKMP2 PES. As we know, the QCT theory has long been used and is able to provide an intuitive picture of the reaction mechanism^{40}, although it will break down owing to the existence of quantum effects. Similar to previous work^{28}, we found that the QCT theory is valid for calculating the product ro-vibrational state-resolved ICSs for both Path 1 and Path 2. This observation confirms that the QCT calculation should be credible for investigating the mechanism of Path 2 in the reaction.

Similar to the previous findings with the H + H_{2} reaction^{23}, the QCT results suggest that the majority of the products were generated through Path 1. However, the reaction mechanism for Path 2 is distinct from that in the H + H_{2} reaction, as the snapshots of representative trajectories for product H_{2} (*v*′ *=* 2, *j*′ = 3) shown in the movie sequence of Fig. 5a–e illustrations. In frame 20–26 fs, the hydrogen H approaches the HD molecule and enters the CI slope region^{41}. Then the HD bond stretches and rotates about its center of mass with the D end toward the incoming H atom, as shown in frame 26–35 fs. At the same time, the incoming H atom is slowed by the linear well, generated by the stretch of the HD molecule^{41}. In frame 35–43 fs, the incoming H atom turns around and passes by the D atom of the HD molecule. In the later frame 43–50 fs, the incoming H atom migrates toward the H atom of the HD molecule, whose bond is much stretched. Finally, as shown in frame 50–58 fs, the incoming H atom impacts at the H atom of the HD molecule to form product molecule as the HD bond breaks. The H_{2} and D products are formed, with product H_{2} scattered into the forward direction. The movie for this representative trajectory of Path 2 is provided as Supplementary Movie 1 for a better understanding. The trajectory is also plotted in hyperspherical coordinates in panel f, where it is seen clearly that the reaction by Path 2 experiences two transition states, as that in the H + H_{2} reaction^{23}. One may compare the trajectory in panel f with the corresponding part of the wave function in Fig. 3. It is noted that for the product H_{2} in other ro-vibrational states, the trajectories are very similar.

This representative trajectory shares some features of the roaming phenomenon in a photodissociation process^{42}: The incoming H atom first migrates around the D atom of HD for some time, but finally impacts the H atom to form the final H_{2} molecule product. Whereas there is only one simple barrier and no any effective potential well in the PES of the reaction, the migration process is transient and lasts for only about 40 fs. Therefore, we name this mechanism for Path 2 as “roaming-like” abstraction. This roaming-like reaction accounts for nearly 100% of the trajectories of Path 2, and some of them are more like insertion reaction and lead to backward scattering, in contrast to the direct insertion mechanism previously found to account for more than 90% of Path 2 in the H + H_{2} reaction.

In summary, this work clearly confirms that the GP effect does exist in the quantum dynamics of chemical reactions at energy well below the CI. The dynamical reason for the existence of the GP effect is arising from the roaming-like abstraction reaction pathway, in addition to the usual direct abstraction reaction pathway, and this roaming-like reaction accounts for nearly 100% of the trajectories of Path 2. In order to accurately describe the reaction involving CI in the adiabatic picture, one has to introduce the vector potential or similar theoretical methods to account for the GP, although the energy is well below the CI. The high-resolution VMI technique is powerful enough for observing the angular oscillations in the product state-resolved DCS, and this opens a door for studying chemical reactions within quantum mechanical principles.

## Methods

### Experimental setup

The CMBs apparatus is equipped with three differentially pumped vacuum chambers: a fixed and a rotatable source chamber, together with a scattering chamber. A high-vacuum condition (10^{−7} ~ 10^{−8} torr with two molecular beams on) can be achieved in the detection area. In this H + HD→H_{2} + D reaction, the H atom beam is produced by dissociating the supersonically expanded HI molecules with a 213 nm laser light. The supersonic HD molecular beam is generated by a liquid nitrogen cooled pulsed valve (Even-Lavie valve). In the HD beam, about 97% of the HD molecules are populated in the (*v* = 0, *j* = 0) level. During the scattering experiment, the pulsed H atom and HD molecule beams collide at a crossing angle of 160°. The D atom products are detected by a two-color (1+1′) (vacuum ultraviolet + ultraviolet) threshold ionization method. The ionized D atom fragments are guided and accelerated by the ion optics^{43} and then detected by a velocity map ion imaging device. For more experimental details, please refer to the Supplementary Note 1. For the experimental statistical error analysis, please refer to Supplementary Note 2.

### Theories

In the quantum dynamics calculations of the H + HD (*v* = 0, *j* = 0) reaction, the product quantum state-resolved *S* matrices are calculated by time-dependent wave packet method on the BKMP2 PES. The wave packet propagation is carried out in hyperspherical coordinates with the second-order split operator^{37}. When the GP effects are included, a vector potential can be conveniently added to the Hamiltonian operator, in a way similar to our previous work using the reactant Jacobi coordinate^{36}. In the calculations, numerical parameters (listed in Supplementary Note 3) can give converged DCS with collision energies up to 4.0 eV.

The QCT calculations for the this reaction are implemented on the BKMP2 PES^{44,45,46,47}. In the calculations, a total of 23.5 million trajectories are propagated at 2.28 eV translational energy giving an overall statistical error of 0.07% in the total reactive cross section. All reactive trajectories can be classified into “Path 1” and “Path 2” paths by accounting how many times each trajectory crossed the transition states, and it is conveniently executed in the hyperspherical coordinate system. In total, 1.892 million trajectories are found to be reactive giving H_{2} product. Of these, ~99.65% reacted through one transition state (the usual direct abstraction mechanism) and only ~0.35% (6661 trajectories) through two transition states.

## Data availability

Data supporting the findings of this study are available from the corresponding authors upon request.

## Code availability

The code that supports the theoretical calculations in this paper is available from the corresponding authors upon request.

## References

- 1.
Berry, M. V. Quantal phase factors accompanying adiabatic changes.

*Proc. R. Soc. Lond., Ser. A***392**, 45–57 (1984). - 2.
Bohm, A., Mostafazadeh, A., Koizumi, H., Niu, Q. & Zwanziger, J.

*The geometric phase in quantum systems*. (Springer Verlag, Berlin, 2003). - 3.
Taguchi, Y., Oohara, Y., Yoshizawa, H., Nagaosa, N. & Tokura, Y. Spin chirality, Berry phase, and anomalous Hall effect in a frustrated ferromagnet.

*Science***291**, 2573–2576 (2001). - 4.
Zwanziger, J. W., Koenig, M. & Pines, A. Berry’s phase.

*Annu. Rev. Phys. Chem.***41**, 601–646 (1990). - 5.
Bohm, A., Kendrick, B. & Loewe, M. E. The Berry phase in molecular physics.

*Int. J. Quan. Chem.***41**, 53–75 (1992). - 6.
Mead, C. A. The geometric phase in molecular systems.

*Rev. Mod. Phys.***64**, 51–85 (1992). - 7.
Xie, C. J., Malbon, C. L., Guo, H. & Yarkony, D. R. Up to a sign. the insidious effects of energetically inaccessible conical intersections on unimolecular reactions.

*Acc. Chem. Res.***52**, 501–509 (2019). - 8.
Applegate, B. E., Barckholtz, T. A. & Miller, T. A. Explorations of conical intersections and their ramifications for chemistry through the Jahn-Teller effect.

*Chem. Soc. Rev.***32**, 38–49 (2003). - 9.
Kendrick, B. Geometric phase effects in the vibrational spectrum of Na

_{3}(*X*).*Phys. Rev. Lett.***79**, 2431–2434 (1997). - 10.
von Busch, H. et al. Unambiguous proof for Berry’s phase in the sodium trimer: Analysis of the transition

*A*^{2}*E*″ ←*X*^{2}*E*′.*Phys. Rev. Lett.***81**, 4584–4587 (1998). - 11.
Nix, M. G. D., Devine, A. L., Dixon, R. N. & Ashfold, M. N. R. Observation of geometric phase effect induced photodissociation dynamics in phenol.

*Chem. Phys. Lett.***463**, 305–308 (2008). - 12.
Xie, C. J. et al. Nonadiabatic tunneling in photodissociation of phenol.

*J. Am. Chem. Soc.***138**, 7828–7831 (2016). - 13.
Domcke, W. & Yarkony, D. R. Role of conical intersections in molecular spectroscopy and photoinduced chemical dynamics.

*Annu. Rev. Phys. Chem.***63**, 325–352 (2012). - 14.
Zhang, J. Z. H., Chu, S.-I. & Miller, W. H. Quantum scattering via the

*S*-matrix version of the Kohn variational principle.*J. Chem. Phys.***88**, 6233–6239 (1988). - 15.
Zhang, J. Z. H. & Miller, W. H. Quantum reactive scattering via the

*S*-matrix version of the Kohn variational principle: Differential and integral cross sections for D+H_{2}→HD+H.*J. Chem. Phys.***91**, 1528–1547 (1989). - 16.
Wrede, E. et al. The dynamics of the hydrogen exchange reaction at 2.20 eV collision energy: Comparison of experimental and theoretical differential cross sections.

*J. Chem. Phys.***110**, 9971–9981 (1999). - 17.
Chao, S. D. et al. A fully state- and angle-resolved study of the H+HD→D+H

_{2}reaction: Comparison of a molecular beam experiment to ab initio quantum reaction dynamics.*J. Chem. Phys.***117**, 8341–8344 (2002). - 18.
Harich, S. A. et al. Forward scattering due to slow-down of the intermediate in the H+HD→D+H

_{2}reaction.*Nature***419**, 281–284 (2002). - 19.
Althorpe, S. C. et al. Observation and interpretation of a time-delayed mechanism in the hydrogen exchange reaction.

*Nature***416**, 67–70 (2002). - 20.
Schneider, L., Seekamp-Rahn, K., Wrede, E. & Welge, K. H. Experimental determination of quantum state resolved differential cross sections for the hydrogen exchange reaction H+D

_{2}→HD+D.*J. Chem. Phys.***107**, 6175–6195 (1997). - 21.
Yuan, D. et al. Direct observation of forward -scattering oscillations in the H+HD→H

_{2}+D reaction.*Nat. Chem.***10**, 653–658 (2018). - 22.
Mead, C. A. & Truhlar, D. G. On the determination of Born-Oppenheimer nuclear motion wave functions including complications due to conical intersections and identical nuclei.

*J. Chem. Phys.***70**, 2284–2296 (1979). - 23.
Juanes-Marcos, J. C., Althorpe, S. C. & Wrede, E. Theoretical study of geometric phase effects in the hydrogen-exchange reaction.

*Science***309**, 1227–1230 (2005). - 24.
Kendrick, B. Geometric phase effects in chemical reaction dynamics and molecular spectra.

*J. Phys. Chem. A***107**, 6739–6756 (2003). - 25.
Juanes-Marcos, J. C. & Althorpe, S. C. Geometric phase effects in the H+H

_{2}reaction: Quantum wave-packet calculations of integral and differential cross sections.*J. Chem. Phys.***122**, 204324 (2005). - 26.
Bouakline, F., Althorpe, S. C. & Ruiz, D. P. Strong geometric-phase effects in the hydrogen-exchange reaction at high collision energies.

*J. Chem. Phys.***128**, 124322 (2008). - 27.
Althorpe, S. C. General explanation of geometric phase effects in reactive systems: unwinding the nuclear wave function using simple topology.

*J. Chem. Phys.***124**, 084105 (2006). - 28.
Juanes-Marcos, J. C., Althorpe, S. C. & Wrede, E. Effect of the geometric phase on the dynamics of the hydrogen-exchange reaction.

*J. Chem. Phys.***126**, 044317 (2007). - 29.
Kendrick, B. K., Hazra, J. & Balakrishnan, N. Geometric phase appears in the ultracold hydrogen exchange reaction.

*Phys. Rev. Lett.***115**, 153201 (2015). - 30.
Kendrick, B. K., Hazra, J. & Balakrishnan, N. The geometric phase controls ultracold chemistry.

*Nat. Comm.***6**, 7918 (2015). - 31.
Strazisar, B. R., Lin, C. & Davis, H. F. Mode-specific energy disposal in the four-atom reaction OH+D

_{2}→HOD+D.*Science***290**, 958–961 (2000). - 32.
Qiu, M. et al. Observation of Feshbach resonances in the F+H

_{2}→HF+H reaction.*Science***311**, 1440–1443 (2006). - 33.
Dai, D. et al. Interference of quantized transition-state pathways in the H+D

_{2}→D+HD chemical reaction.*Science***300**, 1730–1734 (2003). - 34.
Jankunas, J., Sneha, M., Zare, R. N., Bouakline, F. & Althorpe, S. C. Hunt for geometric phase effects in H+HD→HD (

*v*',*j*') + H.*J. Chem. Phys.***139**, 144316 (2013). - 35.
Gao, H., Sneha, M., Bouakline, F., Althorpe, S. C. & Zare, R. N. Differential cross sections for the H+D

_{2}→HD (*v*’ = 3,*j*’ = 4 ~ 10) + D reaction above the conical intersection.*J. Phys. Chem. A***119**, 12036–12042 (2015). - 36.
Yuan, D. F. et al. Direct observation of the geometric phase effect in the H+HD→H

_{2}+D reaction.*Science***362**, 1289–1293 (2018). - 37.
Zhao, H. L., Hu, X. X., Xie, D. Q. & Sun, Z. Quantum wavepacket method for state-to-state reactive cross sections in hyperspherical coordinates.

*J. Chem. Phys.***149**, 174103 (2018). - 38.
Zhao, H. L., Hu, X. X., Xie, D. Q. & Sun, Z. An interaction-asymptotic region decomposition method for general state-to-state reactive scatterings.

*J. Chem. Phys.***150**, 134105 (2019). - 39.
Pack, R. T. & Parker, G. A. Quantum reactive scattering in three dimensions using hyperspherical (APH) coordinates. Theory.

*J. Chem. Phys.***87**, 3888–3921 (1987). - 40.
Jambrina, P. G. et al. Quantum interference between H+D

_{2}quasiclassical reaction mechanisms.*Nat. Chem.***7**, 661–667 (2015). - 41.
Greaves, S. J., Murdock, D., Wrede, E., Althorpe, S. C. J. New, unexpected, and dominant mechanisms in the hydrogen exchange reaction.

*Chem. Phys.***128**, 164306 (2008). - 42.
Bowman, J. M. & Houston, P. L. Theories and simulations of roaming.

*Chem. Soc. Rev.***46**, 7615–7624 (2017). - 43.
Lin, J. J., Zhou, J. G., Shiu, W. C. & Liu, K. P. Application of time-sliced ion velocity imaging to crossed molecular beam experiments.

*Rev. Sci. Instrum.***74**, 2495–2500 (2003). - 44.
Truhlar, D. G. & Muckerman, J. T. in

*Atom-Molecule Collision Theory*, edited by R. B. Bernstein (Plenum, New York, 1979). - 45.
Aoiz, F. J., Herrero, V. J. & Sáez Rábanos, V. Effects of translational, rotational, and vibrational energy on the dynamics of the D+H

_{2}exchange reaction. A classical trajectory study.*J. Chem. Phys.***94**, 7991–8007 (1991). - 46.
Bañares, L., Aoiz, F. J., Honvault, P., Bussery-Honvault, B. & Launay, J.-M. Quantum mechanical and quasi-classical trajectory study of the C (

^{1}*D*)+H_{2}reaction dynamics.*J. Chem. Phys.***118**, 565–568 (2003). - 47.
Hase, W. L. et al. VENUS96: a general chemical dynamics computer program.

*QCPE Bull.***16**, 43 (1996).

## Acknowledgements

This work was supported by the National Key R&D Program of China (grant no. 2016YFF0200500), the National Natural Science Foundation of China (grant no. 21688102, and 21825303), the Strategic Priority Research Program of Chinese Academy of Sciences (grant no. XDB17000000). We are gratefully indebted to Prof. Piergiorgio Casavecchia for helpful discussions.

## Author information

### Affiliations

### Contributions

X.Y., Z.S., and X.W. conceived and supervised the research. The experiments were carried out by D.Y., W.C., S.Y., C.L., Y.T., S.W., and X.W. Data analysis and interpretation were performed by D.Y., W.C., S.Y., C.L., Y.T., S.W., X.W., and X.Y. Theoretical calculations were performed by Y.H., H.Z., and Z.S. The manuscript was written by X.Y., Z.S., and X.W., with contributions from all authors. All authors contributed to discussions about the content of the paper.

### Corresponding authors

## Ethics declarations

### Competing interests

The authors declare no competing interests.

## Additional information

**Peer review information** *Nature Communications* thanks Robert Field, Richard Zare and the other, anonymous, reviewer for their contribution to the peer review of this work.

**Publisher’s note** Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.

## Supplementary information

## Rights and permissions

**Open Access** This article is licensed under a Creative Commons Attribution 4.0 International License, which permits use, sharing, adaptation, distribution and reproduction in any medium or format, as long as you give appropriate credit to the original author(s) and the source, provide a link to the Creative Commons license, and indicate if changes were made. The images or other third party material in this article are included in the article’s Creative Commons license, unless indicated otherwise in a credit line to the material. If material is not included in the article’s Creative Commons license and your intended use is not permitted by statutory regulation or exceeds the permitted use, you will need to obtain permission directly from the copyright holder. To view a copy of this license, visit http://creativecommons.org/licenses/by/4.0/.

## About this article

### Cite this article

Yuan, D., Huang, Y., Chen, W. *et al.* Observation of the geometric phase effect in the H+HD→H_{2}+D reaction below the conical intersection.
*Nat Commun* **11, **3640 (2020). https://doi.org/10.1038/s41467-020-17381-4

Received:

Accepted:

Published:

## Comments

By submitting a comment you agree to abide by our Terms and Community Guidelines. If you find something abusive or that does not comply with our terms or guidelines please flag it as inappropriate.