Geometric phase effects in H+O2 scattering. I. Surface function solutions in the presence of a conical intersection

Impact of Geometric Phase on Dynamics of Complex-Forming Reactions: H + O2 → OH + O

Reaction dynamics on the ground electronic state might be significantly influenced by conical intersections (CIs) via the geometric phase (GP), as demonstrated for activated reactions (i.e., the H + H2 exchange reaction). However, there have been few investigations of GP effects in complex-forming reactions. Here, we report a full quantum dynamical study of an important reaction in combustion (H + O2 → OH + O), which serves as a proving ground for studying GP effects therein. The results reveal significant differences in reaction probabilities and differential cross sections (DCSs) obtained with and without GP, underscoring its strong impact. However, the GP effects are less pronounced for the reaction integral cross sections, apparently due to the integral of the DCS over the scattering angle. Further analysis indicated that the cross section has roughly the same contributions from the two topologically distinct paths around the CI, namely, the direct and looping paths.

Exploring Exact-Factorization-Based Trajectories for Low-Energy Dynamics near a Conical Intersection

We study low-energy dynamics generated by a two-dimensional two-state Jahn-Teller Hamiltonian in the vicinity of a conical intersection using quantum wave packet and trajectory dynamics. Recently, these dynamics were studied by comparing the adiabatic representation and the exact factorization, with the purpose to highlight the different nature of topological-phase and geometric-phase effects arising in the two theoretical representations of the same problem. Here, we employ the exact factorization to understand how to accurately model low-energy dynamics in the vicinity of a conical intersection using an approximate description of the nuclear motion that uses trajectories. We find that since nonadiabatic effects are weak but non-negligible, the trajectory-based description that invokes the classical approximation struggles to capture the correct behavior.

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.

Non-adiabatic quantum interference in the ultracold Li + LiNa → Li2 + Na reaction

Electronically non-adiabatic effects play an important role in many chemical reactions. However, how these effects manifest in cold and ultracold chemistry remains largely unexplored. Here for the first time we present from first principles the non-adiabatic quantum dynamics of the reactive scattering between ultracold alkali-metal LiNa molecules and Li atoms. We show that non-adiabatic dynamics induces quantum interference effects that dramatically alter the ultracold rotationally resolved reaction rate coefficients. The interference effect arises from the conical intersection between the ground and an excited electronic state that is energetically accessible even for ultracold collisions. These unique interference effects might be exploited for quantum control applications such as a quantum molecular switch. The non-adiabatic dynamics are based on full-dimensional ab initio potential energy surfaces for the two electronic states that includes the non-adiabatic couplings and an accurate treatment of the long-range interactions. A statistical analysis of rotational populations of the Li2 product reveals a Poisson distribution implying the underlying classical dynamics are chaotic. The Poisson distribution is robust and amenable to experimental verification and appears to be a universal property of ultracold reactions involving alkali metal dimers.

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 in Ultracold Chemical Reactions

The role of the geometric phase effect in chemical reaction dynamics has long been a topic of active experimental and theoretical investigations. The topic has received renewed interest in recent years in cold and ultracold chemistry where it was shown to play a decisive role in state-to-state chemical dynamics. We provide a brief review of these developments focusing on recent studies of O + OH and hydrogen exchange in the H + H 2 and D + HD reactions at cold and ultracold temperatures. Non-adiabatic effects in ultracold chemical dynamics arising from the conical intersection between two electronic potential energy surfaces are also briefly discussed. By taking the hydrogen exchange reaction as an illustrative example it is shown that the inclusion of the geometric phase effect captures the essential features of non-adiabatic dynamics at collision energies below the conical intersection.

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.

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.

Properties of Feshbach and “shape”-resonances in ozone and their role in recombination reactions and anomalous isotope effects

Three reaction pathways for formation of symmetric and asymmetric isotopologues of ozone.

Symmetry and the geometric phase in ultracold hydrogen-exchange reactions

Quantum reactive scattering calculations are reported for the ultracold hydrogen-exchange reaction and its non-reactive atom-exchange isotopic counterparts, proceeding from excited rotational states. It is shown that while the geometric phase (GP) does not necessarily control the reaction to all final states, one can always find final states where it does. For the isotopic counterpart reactions, these states can be used to make a measurement of the GP effect by separately measuring the even and odd symmetry contributions, which experimentally requires nuclear-spin final-state resolution. This follows from symmetry considerations that make the even and odd identical-particle exchange symmetry wavefunctions which include the GP locally equivalent to the opposite symmetry wavefunctions which do not. It is shown how this equivalence can be used to define a constant which quantifies the GP effect and can be obtained solely from experimentally observable rates. This equivalence reflects the important role that discrete symmetries play in ultracold chemistry and highlights the key role that ultracold reactions can play in understanding fundamental aspects of chemical reactivity more generally.

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.

Quantum dynamics study on the CHIPR potential energy surface for the hydroperoxyl radical: the reactions O + OH⇋O2 + H

Quantum scattering calculations of the O((3)P)+OH((2)Π)⇌O2((3)Σg (-))+H((2)S) reactions are presented using the combined-hyperbolic-inverse-power-representation potential energy surface [A. J. C. Varandas, J. Chem. Phys. 138, 134117 (2013)], which employs a realistic, ab initio-based, description of both the valence and long-range interactions. The calculations have been performed with the ABC time-independent quantum reactive scattering computer program based on hyperspherical coordinates. The reactivity of both arrangements has been investigated, with particular attention paid to the effects of vibrational excitation. By using the J-shifting approximation, rate constants are also reported for both the title reactions.

A complex guided spectral transform Lanczos method for studying quantum resonance states

A complex guided spectral transform Lanczos (cGSTL) algorithm is proposed to compute both bound and resonance states including energies, widths, and wavefunctions. The algorithm comprises of two layers of complex-symmetric Lanczos iterations. A short inner layer iteration produces a set of complex formally orthogonal Lanczos polynomials. They are used to span the guided spectral transform function determined by a retarded Green operator. An outer layer iteration is then carried out with the transform function to compute the eigen-pairs of the system. The guided spectral transform function is designed to have the same wavefunctions as the eigenstates of the original Hamiltonian in the spectral range of interest. Therefore, the energies and/or widths of bound or resonance states can be easily computed with their wavefunctions or by using a root-searching method from the guided spectral transform surface. The new cGSTL algorithm is applied to bound and resonance states of HO2, and compared to previous calculations.

Analysis of geometric phase effects in the quantum-classical Liouville formalism

We analyze two approaches to the quantum-classical Liouville (QCL) formalism that differ in the order of two operations: Wigner transformation and projection onto adiabatic electronic states. The analysis is carried out on a two-dimensional linear vibronic model where geometric phase (GP) effects arising from a conical intersection profoundly affect nuclear dynamics. We find that the Wigner-then-Adiabatic (WA) QCL approach captures GP effects, whereas the Adiabatic-then-Wigner (AW) QCL approach does not. Moreover, the Wigner transform in AW-QCL leads to an ill-defined Fourier transform of double-valued functions. The double-valued character of these functions stems from the nontrivial GP of adiabatic electronic states in the presence of a conical intersection. In contrast, WA-QCL avoids this issue by starting with the Wigner transform of single-valued quantities of the full problem. As a consequence, GP effects in WA-QCL can be associated with a dynamical term in the corresponding equation of motion. Since the WA-QCL approach uses solely the adiabatic potentials and non-adiabatic derivative couplings as an input, our results indicate that WA-QCL can capture GP effects in two-state crossing problems using first-principles electronic structure calculations without prior diabatization or introduction of explicit phase factors.

Geometric phase effects in dynamics near conical intersections: symmetry breaking and spatial localization

We show that finite systems with conical intersections can exhibit spontaneous symmetry breaking which manifests itself in spatial localization of eigenstates. This localization has a geometric phase origin and is robust against variation of model parameters. The transition between localized and delocalized eigenstate regimes resembles a continuous phase transition. The localization slows down the low-energy quantum nuclear dynamics at low temperatures.

Geometric phase for collinear conical intersections. I. Geometric phase angle and vector potentials

We present a method for properly treating collinear conical intersections in triatomic systems. The general vector potential (gauge theory) approach for including the geometric phase effects associated with collinear conical intersections in hyperspherical coordinates is presented. The current study develops an introductory method in the treatment of collinear conical intersections by using the phase angle method. The geometric phase angle, η, in terms of purely internal coordinates is derived using the example of a spin-aligned quartet lithium triatomic system. A numerical fit and thus an analytical form for the associated vector potentials are explicitly derived for this triatomic A(3) system. The application of this methodology to AB(2) and ABC systems is also discussed.

Curl condition for a four-state Born-Oppenheimer system employing the Mathieu equation

When a group of four states forms a subspace of the Hilbert space, i.e., appears to be strongly coupled with each other but very weakly interacts with all other states of the entire space, it is possible to express the nonadiabatic coupling (NAC) elements either in terms of s or in terms of electronic basis function angles, namely, mixing angles presumably representing the same sub-Hilbert space. We demonstrate that those explicit forms of the NAC terms satisfy the curl conditions--the necessary requirements to ensure the adiabatic-diabatic transformation in order to remove the NAC terms (could be often singular also at specific point(s) or along a seam in the configuration space) in the adiabatic representation of nuclear SE and to obtain the diabatic one with smooth functional form of coupling terms among the electronic states. In order to formulate extended Born-Oppenheimer (EBO) equations [J. Chem. Phys. 2006, 124, 074101] for a group of four states, we show that there should exist a coordinate independent ratio of the gradients for each pair of ADT/mixing angles leading to zero curls and, thereafter, provide a brief discussion on its analytical validity. As a numerical justification, we consider the first four eigenfunctions of the Mathieu equation to demonstrate the interesting features of nonadiabatic coupling (NAC) elements, namely, the validity of curl conditions and the nature of curl equations around CIs.

Global analytical potential energy surfaces for HO2(X2A") based on high-level ab initio calculations

Two global analytical potential energy surfaces for the HO2(X2A") system have been developed by fitting approximately 15,000 ab initio points at the icMRCI+Qaug-cc-pVQZ level of theory, using the reproducing kernel Hilbert space method. One analytical potential is designed to give a very accurate representation of the low energy range that determines the vibrational spectrum, while the other attempts to provide a fast and uniformly accurate potential function for reaction dynamics. The quality of the fitted potential functions is confirmed by good agreement of the (J=0) HO2 vibrational spectrum and (J=0) quantum reaction probability for the H+O2(ji=0,nui=0) reaction with those obtained using the spline fitted potential. Quasiclassical trajectory calculations carried out on the new potential energy surface provided the reaction probability with a zero impact parameter (b=0), which is in reasonably good agreement with the J=0 quantum results.

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.

Conical intersection between the lowest spin-aligned Li3(A′4) potential-energy surfaces

We have calculated new potential-energy surfaces for the lowest two spin-aligned (4)A(') states of the Li(3) trimer. This calculation shows a seam of conical intersections between these states resulting from the extra symmetry of the system when the atoms are in a collinear arrangement. This seam is especially important because of its proximity to the three-body dissociation limit of the system; ultracold scattering calculations and the bound-state energies of the system will be affected by the presence of this conical intersection. In this paper we discuss the calculation of the potential-energy surface and the location of the conical intersection seam.

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.