Hyperfine-to-rotational energy transfer in ultracold atom-molecule collisions of Rb and KRb

Energy transfer between different mechanical degrees of freedom in atom-molecule collisions has been studied and largely understood. However, systems involving spins remain less explored. In this study, we directly observed energy transfer from atomic hyperfine to molecular rotation in the 87Rb ( ∣ F a , M F a = ∣ 2 , 2 ) + 40K87Rb (X1Σ+, rotational state N = 0) ⟶ Rb ( ∣ 1 , 1 ) + KRb (N = 0, 1, 2) collision with state-to-state precision. We also performed quantum scattering calculations that rigorously included the coupling between spin and rotational degrees of freedom at short range under the assumption of rigid-rotor KRb monomers moving along a single potential energy surface. The calculated product rotational state distribution deviates from the observations even after extensive tuning of the atom-molecule potential energy surface. In addition, our ab initio calculations indicate that spin-rotation coupling is enhanced close to a conical intersection that is energetically accessible at short range. This, together with the deviation, suggests that vibrational degrees of freedom and conical intersections play an important part in the coupling. Our observations confirm that spin is coupled to mechanical rotation at short range and establish a benchmark for future theoretical studies.

Nonadiabatic Effects in the H + LiD(ν = 0, j = 0) → Li(2s) + HD Reaction Near Cold Collisions

Nonadiabatic dynamic study of the H + LiD(ν = 0, j = 0) → Li(2s) + HD reaction is carried out using the time-dependent wave packet method in a collision energy range of 1-80 cm-1. The total integral cross section exhibits a partial wave resonance near 2 cm-1, corresponding to the opening of the J = 5 partial wave. The nonadiabatic coupling effects inhibit the reactivity, especially for the low-vibrational states. The rotational excitation of products is affected by nonadiabatic coupling effects. The maximum accessible rotational state of the products is higher when nonadiabatic effects are included than when they are omitted. At low collision energies, the product angular distributions are influenced by the resonances. Nonadiabatic results reveal a more pronounced backward scattering of the products than adiabatic results. As collision energy increases, the stripping mechanism gradually becomes dominant, and both adiabatic and nonadiabatic results exhibit significant forward-scattering characteristics.

Reaction Path-Resolved Quantum Transition State Framework Using Hyperspherical (APH) Coordinates: The Geometric Phase Effects in the H + H2 Reaction

The quantum transition state framework was developed to calculate the reaction path-resolved scattering matrix for atom-diatom reactions in hyperspherical (APH) coordinates. This approach allows for simply and directly calculating the reaction path-resolved scattering matrix, especially when the encircling reaction path is negligible. It could be used to determine the reactivities of specific pathways in a chemical reaction, providing insights into phenomena such as geometric phase effects. To validate our method, detailed calculations for the H + H2 reaction at J = 0, 1, 2 were carried out, and the results were compared with those from our previous theoretical models including the geometric phase effects. Analysis of the path-resolved reaction probabilities reveals that encircling paths play a very minor role in this reaction, even at energies exceeding the conical intersection minimum. Furthermore, the analysis suggests that the geometric phase effects arising from the interferences between different reactive pathways could become more significant for highly rotationally excited reactants and products.

Unraveling the effect of reagent vibrational excitation on the scattering mechanism of the benchmark H + H2 → H2 + H hydrogen exchange reaction in the coupled 12E' ground electronic manifold

The hydrogen exchange reaction, H + H2 → H2 + H, along with its isotopic variants, has been the cornerstone for the development of new and novel dynamical mechanisms of gas-phase bimolecular reactions since the 1930s. The dynamics of this reaction are theoretically investigated in this work to elucidate the effect of reagent vibrational excitation on differential cross sections (DCSs) in a nonadiabatic situation. The dynamical calculations are carried out using a time-dependent quantum mechanical method, both on the lower adiabatic potential energy surface and employing a two-state coupled diabatic theoretical model to explicitly include all the nonadiabatic couplings present in the 12E' ground electronic manifold of the H3 system. Towards this effort, the Boothroyd-Keogh-Martin-Peterson (BKMP2) surface of the lower adiabatic component is coupled with the double many-body expansion (DMBE) surface of the upper one. The smooth variation of energy along the D3h seam of the conical intersections is explicitly confirmed. The coupled two-state calculations are performed only for H2 (v = 3-4, j = 0), as the minimum of the conical intersections becomes energetically accessible for these vibrational levels in the considered energy range. Initial state-selected total and state-to-state DCSs are calculated to elucidate various mechanistic aspects of reagent vibrational excitation. The latter enhances the forward scattering and makes the backward scattering less prominent. Important roles of collision energy in the vibrational energy disposal of both forward- and backward-scattered products are examined. Analysis of the state-to-state DCSs of the vibrationally excited reagents reveals an important correlation among scattering angle, and the product rotational angular momentum and its helicity state. Such an analysis establishes a novel mechanism for the forward scattering of the reaction.

Vibrational State-Resolved Differential Cross Sections of the F + CH4 → HF + CH3 Reaction at the Collision Energies of 3.1-13.8 kcal/mol

We report a high-resolution crossed molecular beam experiment investigating the reaction dynamics of the F + CH4 → HF + CH3 reaction across a broad range of collision energies (3.1-13.8 kcal/mol). By using time-sliced velocity map imaging within the crossed molecular beam apparatus, we obtain correlated angular distributions and branching ratios for various product pairs (CH3(ν), HF(ν')). The resolved reactive rainbow-like features that display a distinct bulge in the angular distribution in different channels reveal the formation of vibrationally ground and excited states of the methyl radical accompanied by different rovibrationally excited HF products. These findings suggest distinct reaction dynamics for channels leading to vibrationally excited methyl radicals compared to those forming ground-state methyl radicals.

Profiling a pulsed molecular beam with cavity-enhanced absorption spectroscopy

The molecular beam plays an important role in chemical dynamics experiments. The density in the beam is one of the critical factors influencing the reaction rate in these studies. Here we present a method based on laser-locked cavity-enhanced absorption spectroscopy to measure the molecular density in the beam. The P(1) transition in the second overtone band of CO was measured in the molecular beam, demonstrating a determination of the number density of molecules in a specific quantum state from the absorption spectrum. This non-destructive spectroscopic method allows the measurement of state-resolved properties of a molecular beam, which could be applied to various studies such as molecular collision dynamics.

Extracting the Geometric Phase from the Ensemble of Trajectories

Traditionally, methods designed to investigate the effects of the geometric phase in reaction dynamics, such as including a vector potential in the nuclear Hamiltonian, necessitate the explicit manipulation of geometric phase-related terms in the adiabatic representation. In contrast, the diabatic representation provides an alternative approach that implicitly addresses the geometric phase and nonadiabatic issues. In this study, we present a method to directly extract the phase information on the geometric phase from the ensemble of interdependent trajectories utilizing the diabatic representation. This approach presents a direct means of quantitatively examining the geometric phase effects in dynamics and has the potential to yield observables suitable for experimental measurement.

Simulations of photoinduced processes with the exact factorization: state of the art and perspectives

This perspective offers an overview of the applications of the exact factorization of the electron-nuclear wavefunction to the domain of theoretical photochemistry, where the aim is to gain insights into the ultrafast dynamics of molecular systems via simulations of their excited-state dynamics beyond the Born-Oppenheimer approximation. The exact factorization offers an alternative viewpoint to the Born-Huang representation for the interpretation of dynamical processes involving the electronic ground and excited states as well as their coupling through the nuclear motion. Therefore, the formalism has been used to derive algorithms for quantum molecular-dynamics simulations where the nuclear motion is treated using trajectories and the electrons are treated quantum mechanically. These algorithms have the characteristic features of being based on coupled and on auxiliary trajectories, and have shown excellent performance in describing a variety of excited-state processes, as this perspective illustrates. We conclude with a discussion on the authors' point of view on the future of the exact factorization.

Quasi-Diabatization Based on Minimizing Derivative Couplings in a Limited Configuration Space: Elimination of Boundary Condition Dependence

Due to the cuspidal ridges of adiabatic potential energy surfaces (PESs) and singularities of nonadiabatic couplings (NACs), obtaining an analytical expression for the adiabatic Hamiltonian is difficult. Thereby, nonadiabatic dynamics simulations are often carried out on-the-fly, which is time-consuming. This motivates us to construct quasi-diabatic representations, which have smooth PESs and diabatic couplings. In this study, we propose a new quasi-diabatization method based on minimizing derivative couplings (MDC) in a limited configuration space. The boundary conditions are first considered and finally released to obtain the adiabatic-to-diabatic rotation angles and transformation matrices. As demonstrated in representative one- and two-dimensional models and the widely studied linear H3 molecule, MDC performs significantly better than the direct integration quasi-diabatization approach. In particular, accurate diabatic potential energy matrices have been successfully obtained even when the NACs of all configurations in the considered space are nonnegligible.

The Interaction-Asymptotic Region Decomposition Method in Jacobi Coordinates: Triatomic Reactive Scatterings

The interaction-asymptotic region decomposition (IARD) technique has been proven to be a good solution to the long-standing coordinate problem in reactive scattering calculations. In this work, the IARD technique was further developed using Jacobi coordinates for the interaction region, instead of the previously used hyperspherical coordinates. Although the Jacobi coordinate may not be as optimal as the hyperspherical coordinates for describing the interaction region in reactive scattering processes, it has simpler kinetic operators and provides a more physically intuitive picture. By developing an intermediate interpolation method, which could efficiently transform the overlapped wave functions from the asymptotic regions to the interaction region, the new implementation of the IARD technique for triatomic reactive scatterings is similarly efficient and accurate. The differential cross sections of the H+H2, and product state-resolved reaction probabilities of the F+HD and16O+36O2 reactions, which involve products of extremely low translational energy and are challenging for a single coordinate-based method, were calculated as numerical examples to show the ability of the new method.

Frustrated charge transfer in vibrationally inelastic Ar++N2 collisions via hard collision glory scattering

Vibrational energy transfer in collisions between ions and neutrals is a fundamental process in interstellar media, planetary atmospheres, and plasmas. The conventional wisdom is that glancing collisions with large impact parameters are forward-scattered with low vibrational excitation, while hard collisions with small impact parameters are sideway- or backward-scattered with relatively high vibrational excitation. Here, we report experimental observations with a three-dimensional velocity-map imaging crossed-beam apparatus in the inelastic scattering process Ar++N2(v'' = 0, J'')→Ar++N2(v', J'), where all the vibrationally excited N2 products are dominated by forward scattering, contradicting the textbook model. Trajectory surface hopping calculations not only reproduced the experimental observation qualitatively, but also revealed that the vibrational excitation mainly occurs through a transient charge-transfer process. The hard collision glory mechanism, which has so far only been observed in inelastic rotational energy transfer between neutrals, is shown to play a major role for vibrational excitation in the inelastic Ar++N2 collision, via the frustrated charge transfer process.

Spiers Memorial Lecture: New directions in molecular scattering

The field of molecular scattering is reviewed as it pertains to gas-gas as well as gas-surface chemical reaction dynamics. We emphasize the importance of collaboration of experiment and theory, from which new directions of research are being pursued on increasingly complex problems. We review both experimental and theoretical advances that provide the modern toolbox available to molecular-scattering studies. We distinguish between two classes of work. The first involves simple systems and uses experiment to validate theory so that from the validated theory, one may learn far more than could ever be measured in the laboratory. The second class involves problems of great complexity that would be difficult or impossible to understand without a partnership of experiment and theory. Key topics covered in this review include crossed-beams reactive scattering and scattering at extremely low energies, where quantum effects dominate. They also include scattering from surfaces, reactive scattering and kinetics at surfaces, and scattering work done at liquid surfaces. The review closes with thoughts on future promising directions of research.

Making Peace with Random Phases: Ab Initio Conical Intersection Quantum Dynamics in Random Gauges

Ab initio modeling of conical intersection wave packet dynamics is crucial for various photochemical, photophysical, and biological processes. However, adiabatic electronic states obtained from electronic structure computations involve random phases, or more generally, random gauge fixings, which cannot be directly used in the modeling of nonadiabatic wave packet simulations. Here we develop a random-gauge local diabatic representation that allows an exact modeling of conical intersection dynamics directly using the adiabatic electronic states with phases randomly assigned during the electronic structure computations. Its utility is demonstrated by an exact ab initio modeling of the two-dimensional Shin-Metiu model with and without an external magnetic field. Our results provide a simple approach to integrating the electronic structure computations into nonadiabatic quantum dynamics, thus paving the way for ab initio modeling of conical intersection dynamics.

Accurate and Efficient Schemes for Mapping Out Two Isolated Seams of Conical Intersections with Neural Networks Solely Based on Adiabatic Energies: A Case Study of H+ + NO(X2Π) → H + NO+(X1Σ+)

A new scheme in the neural network (NN) diabatization approach that solely utilizes the adiabatic energies for constructing the global diabatic potential energy matrices (PEMs) of the molecular systems with two isolated seams of conical intersections (CIs) is proposed. Taking a prototype charge transfer reaction H+ + NO(X2Π) → H + NO+(X1Σ+), where two seams of CIs are located at the different linear geometries N-O-H and O-N-H, for example, the diabatization with the new scheme including a diabatic state constraint is shown to map out the topographies of both two linear CIs with 100% of the success rate in 10 different trainings, while the diabatization without such constraint hardly represents CIs, in which the avoided crossings appear instead. Simultaneously, we propose a scheme to separate the whole reactive space into three different regions and define the minimal Euclidean distances for each region to efficiently sample the energy points for the NN trainings. Through adjusting the minimal Euclidean distances, the number of the adiabatic energy data needed in the construction of diabatic PEM can be heavily decreased, lowered from ∼2000 to ∼280 energy points, which is much less than the number (>22 000) of ab initio energy points used in earlier spline diabatic PEM. Further quantum dynamic calculations show that the reaction probabilities, vibrational state distributions, and vibrational state resolved differential cross sections are well reproduced on the new NN diabatic PEM, validating these schemes for constructing the diabatic PEM.

Quantum state tracking and control of a single molecular ion in a thermal environment

Understanding molecular state evolution is central to many disciplines, including molecular dynamics, precision measurement, and molecule-based quantum technology. Details of this evolution are obscured when observing a statistical ensemble of molecules. Here, we report real-time observations of thermal radiation-driven transitions between individual states ("jumps") of a single molecule. We reversed these jumps through microwave-driven transitions, which resulted in a 20-fold improvement in the time the molecule dwells in a chosen state. The measured transition rates showed anisotropy in the thermal environment, pointing to the possibility of using single molecules as in situ probes for the strengths of ambient fields. Our approaches for state detection and manipulation could apply to a wide range of species, facilitating their uses in fields including quantum science, molecular physics, and ion-neutral chemistry.

Quantum Dynamics Studies of the Li + Na2 (V = 0, j = 0) → Na + NaLi Reaction on a New Neural Network Potential Energy Surface

The ultracold reaction offers a unique opportunity to elucidate the intricate microscopic mechanism of chemical reactions, and the Na2Li system serves as a pivotal reaction system in the investigation of ultracold reactions. In this work, a high-precision potential energy surface (PES) of the Na2Li system is constructed based on high-level ab initio energy points and the neural network (NN) method, and a proper asymptotic functional form is adopted for the long-range interaction, which is suitable for the study of cold or ultracold collisions. Based on the new NN PES, the dynamics of the Li + Na2 (v = 0, j = 0) → Na + NaLi reaction are studied in the collision energy range of 10-7 to 80 cm-1. In the high collision energy range of 8 to 80 cm-1, the dynamics of the reaction is studied using the time-dependent wave packet method and the statistical quantum mechanical (SQM) method. Comparing the results of the two methods, it is found that the SQM method provides a rough description of the product ro-vibrational state distribution but overestimates the integral cross-section values. With the decrease of collision energy, the reaction differential cross section gradually changes from forward-backward symmetric scattering to predominant forward scattering. In the low collision energy range from 10-7 to 8 cm-1, the SQM method is used to study the reaction dynamics, and the rate constant in the Wigner threshold region is estimated to be 2.87 × 10-10 cm3/s.

High-Resolution Crossed-Beam Dynamics Studies of the D + Para-H2 → HD + H Reaction at 1.21 eV

Understanding kinetic isotope effects is important in the study of the reaction dynamics of elementary chemical reactions, particularly those involving hydrogen atoms and molecules. As one of the isotopic variants of the hydrogen exchange reaction, the D + para-H2 reaction has attracted much attention. However, experimental studies of this reaction have been limited primarily due to its strong experimental background noise. In this study, by using the velocity map ion imaging method and the near-threshold ionization technique, together with improvements on the vacuum condition in the vicinity of the collision zone, background noise was reduced significantly, and quantum state-resolved differential cross sections (DCSs) for the D + para-H2 reaction at a collision energy of 1.21 eV were acquired in a crossed molecular beams experiment. Interestingly, clear rotational state-dependent angular distributions were noticed in the quantum state-resolved DCSs. The most intense peak's positions for HD (v', j') products shift to different scattering directions as the product's ro-vibrational quantum number increases. Two different microscopic reaction mechanisms are found to be involved in this reaction for HD products in different vibrational states. The results show a direct correlation between the scattering angle and the product's rotational quantum number, revealing that the contributions of impact parameters are strongly influenced by the corresponding centrifugal barrier.

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.

Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method

Although the reactant-product decoupling (RPD) technique was proposed over two decades ago, it remains an efficient approach for calculating product state-resolved information on some simple direct reactions using the quantum wave packet method. In the past, usually the RPD technique employed the collocation method to transform the wave function between reactant and product arrangements, which requires quite large computational efforts. In this work, the intermediate coordinate (IC) method is employed to realize the RPD technique. Numerical examples demonstrate that this new IC RPD (IRPD) technique has superior computational efficiency compared with the original method employing the collocation method. Especially, the new IRPD technique significantly saves disk space and computer memory. To illustrate the features of our new method, the total reaction probabilities of the H + H2, H + Br2, and F + H2 reactions with J = 0 and the differential cross sections of the H + H2 and F + H2 reactions at a series of collision energy are calculated and presented. With this efficient and effective new RPD technique, the Li + HF reaction, which involves sharp resonances with long-range wave functions in the van der Waals wells in both the reactant and product arrangements, is also calculated with several J at the product state-resolved level to reveal the ability of the RPD technique for describing resonance wave functions. With these numerical examples, it is found that, for the reaction with resonances, the RPD approach should be applied carefully. Otherwise, it is very possible that the resonances could disappear with the application of the RPD technique.

Seeking a quantum advantage with trapped-ion quantum simulations of condensed-phase chemical dynamics

Simulating the quantum dynamics of molecules in the condensed phase represents a longstanding challenge in chemistry. Trapped-ion quantum systems may serve as a platform for the analog-quantum simulation of chemical dynamics that is beyond the reach of current classical-digital simulation. To identify a 'quantum advantage' for these simulations, performance analysis of both analog-quantum simulation on noisy hardware and classical-digital algorithms is needed. In this Review, we make a comparison between a noisy analog trapped-ion simulator and a few choice classical-digital methods on simulating the dynamics of a model molecular Hamiltonian with linear vibronic coupling. We describe several simple Hamiltonians that are commonly used to model molecular systems, which can be simulated with existing or emerging trapped-ion hardware. These Hamiltonians may serve as stepping stones towards the use of trapped-ion simulators for systems beyond the reach of classical-digital methods. Finally, we identify dynamical regimes in which classical-digital simulations seem to have the weakest performance with respect to analog-quantum simulations. These regimes may provide the lowest hanging fruit to make the most of potential quantum advantages.

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.

Accurate diabatization based on combined-hyperbolic-inverse-power-representation: 1,2 2A' states of BeH2

A diabatic potential energy matrix (DPEM) for the two lowest states of BeH2+ has been constructed using the combined-hyperbolic-inverse-power-representation (CHIPR) method. By imposing symmetry constraints on the coefficients of polynomials, the complete nuclear permutation inversion symmetry is correctly preserved in the CHIPR functional form. The symmetrized CHIPR functional form is then used in the diabatization by ansatz procedure. The ab initio energies are reproduced with satisfactory accuracy. In addition, the CHIPR-based DPEM also reproduces the local topology of a conical intersection. Future work will focus on a complete four-state diabatic representation with emphasis on the long-range interactions and spin-orbit couplings, which will enable accurate quantum scattering calculations for the Be+(2P) + H2 → BeH+(X1Σ+) + H(2S) reaction.

Theoretical Development of the Interaction-Asymptotic Region Decomposition Method for Tetratomic Reactive Scattering

An accurate and efficient time-dependent wave packet method is proposed for solving the product state-resolved reaction probabilities of the tetratomic reactive system. In this method, the entire scattering process is divided into the interaction region and multiple asymptotic regions, sharing the same spirit as the interaction-asymptotic region decomposition (IARD) approach in a triatomic reactive scattering process. The hyperspherical coordinate is adopted in the interaction region, while the corresponding Jacobi coordinate is employed in each asymptotic region. Therefore, in this IARD method, the "coordinate problem", the difficulty of expressing the wave function in the entire region using a single coordinate system, can be effectively avoided, and only a very small number of the grid points (or the basis functions) are required. For the numerical illustration, the typical tetratomic reaction H2 + OH with zero total angular momentum is calculated, and compared with other quantum wave packet methods. Our proposed IARD method for the tetratomic reactive system is much more efficient and accurate.

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.

Unveiling Quantum Interference in the D+ + H2 Nonadiabatic Reaction Dynamics at Low Collision Energies

Fully converged nonadiabatic dynamics calculations of the D+ + H2 → H+ + HD reaction are performed at low temperatures using the time-dependent wave packet approach based on a set of precise 3 × 3 diabatic potential energy surfaces (PESs) ( Phys. Chem. Chem. Phys., 2021, 23, 7735-7747, DOI: 10.1039/D0CP04100A). The D+ + H2 reaction is mediated by a dense manifold of resonances associated with the deep potential well on the ground-state PES. The calculated results show that the nonadiabatic coupling can affect the resonance positions, deviating from the expectation based solely on adiabatic considerations. Furthermore, significant forward-backward asymmetry in total differential cross sections (DCSs) is revealed, which is markedly influenced by nonadiabatic effects. The nonadiabatic effects not only affect the contribution of partial waves in the reaction but also make the interference patterns in the DCSs change significantly.

On the Nature of Geometric and Topological Phases in the Presence of Conical Intersections

The observable nature of topological phases related to conical intersections in molecules is studied. Topological phases should be ubiquitous in molecular processes, but their elusive character has often made them a topic of discussion. To shed some light on this issue, we simulate the dynamics governed by a Jahn-Teller Hamiltonian and analyze it employing two theoretical representations of the molecular wave function: the adiabatic and the exact factorization. We find fundamental differences between effects related to topological phases arising exclusively in the adiabatic representation, and thus not related to any physical observable, and geometric phases within the exact factorization that can be connected to an observable quantity. We stress that while the topological phase of the adiabatic representation is an intrinsic property of the Hamiltonian, the geometric phase of the exact factorization depends on the dynamics that the system undergoes and is connected to the circulation of the nuclear momentum field.

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.

Competing quantum effects in heavy-atom tunnelling through conical intersections

Thermally activated chemical reactions are typically understood in terms of overcoming potential-energy barriers. However, standard rate theories break down in the presence of a conical intersection (CI) because these processes are inherently nonadiabatic, invalidating the Born-Oppenheimer approximation. Moreover, CIs give rise to intricate nuclear quantum effects such as tunnelling and the geometric phase, which are neglected by standard trajectory-based simulations and remain largely unexplored in complex molecular systems. We present new semiclassical transition-state theories based on an extension of golden-rule instanton theory to describe nonadiabatic tunnelling through CIs and thus provide an intuitive picture for the reaction mechanism. We apply the method in conjunction with first-principles electronic-structure calculations to the electron transfer in the bis(methylene)-adamantyl cation. Our study reveals a strong competition between heavy-atom tunnelling and geometric-phase effects.

A Discrete-Variable Local Diabatic Representation of Conical Intersection Dynamics

Conical intersections (CIs) are ubiquitous in polyatomic molecules and are responsible for a wide range of phenomena in photochemistry and photophysics. Modeling the conical intersection dynamics with adiabatic electronic states is hindered by the divergence of the first- and second-order derivative couplings at CIs due to electronic degeneracy. We introduce and implement a novel diabatic representation for exact correlated electron-nuclear wave packet dynamics through conical intersections. It directly employs the adiabatic electronic states but avoids the singular first- and second-order derivative couplings and is robust to different gauge choices of the electronic wave function phases. The reference nuclear geometries defining the adiabatic electronic states are determined by a discrete-variable representation of the nuclear coordinates. The nonadiabatic effects are accounted for by the electronic overlap matrix instead of derivative couplings as in the adiabatic representation. Illustrated by a two-mode conical intersection model, this representation captures all nonadiabatic effects, including electronic transitions, electronic coherence, and geometric phases. Thus, this representation provides a singularity-free framework for modeling ab initio conical intersection wave packet dynamics.

Imaging of the charge-transfer reaction of spin-orbit state-selected Ar+(2P3/2) with N2 reveals vibrational-state-specific mechanisms

Charge-transfer reactions are ubiquitous and play important roles in various gaseous environments, but, despite a long history of extensive research, our understanding of their dynamics at the quantum state-to-state level is still lacking. Here we report quantum-state-resolved experiments for the paradigmatic charge-transfer reaction Ar+ + N2 → Ar + N2+ using a three-dimensional velocity-map imaging crossed-beam apparatus with the Ar+ beam prepared exclusively in the spin-orbit state 2P3/2. High-resolution scattering images show strong dependence of rotational and angular distributions on the vibrational quantum number of the N2+ product. Trajectory surface-hopping calculations, which semi-quantitatively reproduce the experimental observations, support the existence of two distinct charge-transfer mechanisms. One of these, in the dominant N2+(v' = 1) channel, is the well-known long-distance harpooning mechanism. However, the highly rotationally excited products in the forward direction are attributed to a hard-collision glory scattering mechanism, which occurs on account of the strong attraction between the collisional partners counterbalanced by the short-range repulsive interaction.

Exploring the vacuum ultraviolet photochemistry of astrochemically important triatomic molecules

The recently constructed vacuum ultraviolet (VUV) free electron laser (FEL) at the Dalian Coherent Light Source (DCLS) is yielding a wealth of new and exquisitely detailed information about the photofragmentation dynamics of many small gas-phase molecules. This Review focuses particular attention on five triatomic molecules-H2O, H2S, CO2, OCS and CS2. Each shows excitation wavelength-dependent dissociation dynamics, yielding photofragments that populate a range of electronic and (in the case of diatomic fragments) vibrational and rotational quantum states, which can be characterized by different translational spectroscopy methods. The photodissociation of an isolated molecule from a well-defined initial quantum state provides a lens through which one can investigate how and why chemical reactions occur, and provides numerous opportunities for fruitful, synergistic collaborations with high-level ab initio quantum chemists. The chosen molecules, their photofragments and the subsequent chemical reaction networks to which they can contribute are all crucial in planetary atmospheres and in interstellar and circumstellar environments. The aims of this Review are 3-fold: to highlight new photochemical insights enabled by the VUV-FEL at the DCLS, notably the recently recognized central atom elimination process that is shown to contribute in all of these triatomic molecules; to highlight some of the potential implications of this rich photochemistry to our understanding of interstellar chemistry and molecular evolution within the universe; and to highlight other and future research directions in areas related to chemical reaction dynamics and astrochemistry that will be enabled by increased access to VUV-FEL sources.

On the geometric phases during radio frequency pulses with sine and cosine amplitude and frequency modulation

In this work, we describe the formation of geometric phases during nonadiabatic frequency swept (FS) radio frequency (RF) pulses with sine amplitude modulation and cosine frequency modulation functions. The geometric phases during the FS pulse were analyzed using a Schrödinger equation formalism, and the unified analytical expression for the geometric phase was derived. We present the solutions for sub-geometric phase components incorporated in spinor wavefunctions for the RF Hamiltonian of spin ½ nuclei. We demonstrate that the geometric phases during sine/cosine RF pulses are opposite in signs for different initial conditions of the spinor and that geometric phases can accumulate in correspondence to different magnetization trajectories. The derived formalism could be extended for the evaluation of the geometric phases during a wide class of amplitude- and frequency-modulated pulses used in MRI and in high-resolution NMR.

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.

Constructing Diabatic Potential Energy Matrices with Neural Networks Based on Adiabatic Energies and Physical Considerations: Toward Quantum Dynamic Accuracy

A permutation invariant polynomial-neural network (PIP-NN) approach for constructing the global diabatic potential energy matrices (PEMs) of the coupled states of molecules is proposed. Specifically, the diabatization scheme is based merely on the adiabatic energy data of the system, which is ideally a most convenient way due to not requiring additional ab initio calculations for the data of the derivative coupling or any other physical properties of the molecule. Considering the permutation and coupling characteristics of the system, particularly in the presence of conical intersections, some vital treatments for the off-diagonal terms in diabatic PEM are essentially needed. Taking the photodissociation of H₂O(X~/B~)/NH₃(X~/A~) and nonadiabatic reaction Na(3p) + H₂ → NaH(Σ⁺) + H for example, this PIP-NN method is shown to build up the global diabatic PEMs effectively and accurately. The root-mean-square errors of the adiabatic potential energies in the fitting for three different systems are all small (<10 meV). Further quantum dynamic calculations show that the absorption spectra and product branching ratios in both H₂O(X~/B~) and NH₃(X~/A~) nonadiabatic photodissociation are well reproduced on the new diabatic PEMs, and the nonadiabatic reaction probability of Na(3p) + H₂ → NaH(Σ⁺) + H obtained on the new diabatic PEMs of the 1²A₁ and 1²B₂ states is in reasonably good agreement with previous theoretical result as well, validating this new PIP-NN method.

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.

Ultrafast Imaging of the Jahn-Teller Topography in Carbon Tetrachloride

We report the direct time-domain observation of ultrafast dynamics driven by the Jahn-Teller effect. Using time-resolved photoelectron spectroscopy with a vacuum-ultraviolet femtosecond source to prepare high-lying Rydberg states of carbon tetrachloride, our measurements reveal the local topography of a Jahn-Teller conical intersection. The pump pulse prepares a configurationally mixed superposition of the degenerate ¹T₂ 4p-Rydberg states, which then distorts through spontaneous symmetry breaking that we identify to follow the t₂ bending motion. Photoionization of these states to three cationic states ²T₁, ²T₂, and ²E reveals a shift in the center-of-mass of the photoelectron peaks associated with the ²T_n states which reveals the local topography of the Jahn-Teller conical intersection region prepared by the pump pulse. Time-dependent density functional theory calculations confirm that the dominant nuclear motion observed in the spectrum is the CCl₄ t₂ bending mode. The large density of states in the VUV spectral region at 9.33 eV of carbon tetrachloride and strong vibronic coupling result in ultrafast decay of the excited-state signal with a time constant of 75(4) fs.

Signature of quantum interference effect in inter-layer Coulomb drag in graphene-based electronic double-layer systems

The distinguishing feature of a quantum system is interference arising from the wave mechanical nature of particles which is clearly central to macroscopic electronic properties. Here, we report the signature of quantum interference effect in inter-layer transport process. Via systematic magneto-drag experiments on graphene-based electronic double-layer systems, we observe low-field correction to the Coulomb-scattering-dominated inter-layer drag resistance in a wide range of temperature and carrier density, with its characteristics sensitive to the band topology of graphene layers. These observations can be attributed to a new type of quantum interference between drag processes, with the interference pathway comprising different carrier diffusion paths in the two constituent conductors. The emergence of such effect relies on the formation of superimposing planar diffusion paths, among which the impurity potentials from intermediate insulating spacer play an essential role. Our findings establish an ideal platform where the interplay between quantum interference and many-body interaction is essential.

Stereodynamical control of the H + HD → H2 + D reaction through HD reagent alignment

Prealigning nonpolar reacting molecules leads to large stereodynamical effects because of their weak steering interaction en route to the reaction barrier. However, experimental limitations in preparing aligned molecules efficiently have hindered the investigation of steric effects in bimolecular reactions involving hydrogen. Here, we report a high-resolution crossed-beam study of the reaction H + HD(v = 1, j = 2) → H₂(v', j') + D at collision energies of 0.50, 1.20, and 2.07 electron volts in which the vibrationally excited hydrogen deuteride (HD) molecules were prepared in two collision configurations, with their bond preferentially aligned parallel and perpendicular to the relative velocity of collision partners. Notable stereodynamical effects in differential cross sections were observed. Quantum dynamics calculations revealed that strong constructive interference in the perpendicular configuration plays an important role in the stereodynamical effects observed.

Electronically Nonadiabatic Effects on the Quantum Dynamics of the Ha + BeHb+ → Be+ + HaHb; Hb + BeHa+ Reactions

Nonadiabatic effects are ubiquitous and play an important role in many chemical processes. Here, the adiabatic and nonadiabatic quantum scattering calculations of the H + BeH⁺ reaction are performed using the time-dependent wave packet method based on an accurate diabatic potential energy matrix that includes the lowest two electronic states and their couplings. The resulting integral cross sections reveal that the nonadiabatic effect significantly inhibits the reactivity of the BeH⁺-depletion channel but enhances that of the H-exchange channel. The vibrational excitation is suppressed, but the translational excitation is promoted for the H₂ product in the BeH⁺-depletion channel when the nonadiabatic coupling is included. However, the nonadiabatic coupling has a mild effect on the H-exchange product-state distribution. When the nonadiabatic effect is considered, the differential cross sections of the H₂ product become less polarized because of the formation of an excited-state complex, whereas the corresponding results of the H-exchange channel only present an increase in the magnitude at the backward region.

Crossed Molecular Beam Study of the H + HD → H2 + D Reaction at 0.60 and 1.26 eV Using the Near-Threshold Ionization Velocity Map Ion Imaging

By using the 1 + 1' near-threshold ionization velocity map ion imaging technique, state-to-state reactive differential cross sections have been measured for the H + HD → H₂ + D reaction. High-resolution images of the D products, with the rotational states of the H₂ co-products clearly resolved, were acquired at the collision energies of 0.60 and 1.26 eV, respectively. It is found that the angular distribution is predominantly backward-scattering at the collision energy of 0.60 eV. However, at 1.26 eV, where the collision energy is higher, the angular distribution becomes forward-backward-scattering. Notably, at both collision energies, the main peaks of backward-scattered products gradually shift from backward toward sideways direction as the rotational quantum number of H₂ increases. Moreover, in the forward direction, fast angular oscillations, which are induced by specific partial waves have also been observed at 1.26 eV. These features show a strong correlation between the product states and angular distributions and also indicate the unique role of partial waves in quantum reactive scattering.

Chemistry of a Light Impurity in a Bose-Einstein Condensate

Similar to an electron in a solid, an impurity in an atomic Bose-Einstein condensate (BEC) is dressed by excitations from the medium, forming a polaron quasiparticle with modified properties. This impurity can also undergo chemical recombination with atoms from the BEC, a process resonantly enhanced when universal three-body Efimov bound states cross the continuum. To study the interplay between these phenomena, we use a Gaussian state variational method able to describe both Efimov physics and arbitrarily many excitations of the BEC. We show that the polaron cloud contributes to bound state formation, leading to a shift of the Efimov resonance to smaller interaction strengths. This shifted scattering resonance marks the onset of a polaronic instability towards the decay into large Efimov clusters and fast recombination, offering a remarkable example of chemistry in a quantum medium.

Stereodynamics-Controlled Product Branching in the Nonadiabatic H + NaD → Na(3s, 3p) + HD Reaction at Low Temperatures

Nonadiabatic processes play an important role at energies near or higher than conical intersection of adiabatic potential energy surfaces in chemical reactions. In this work, dynamics of the nonadiabatic H + NaD reaction at low temperatures are studied by using the quantum wave packet method based on an improved L-shaped grid. The nonadiabatic H + NaD reaction has two exothermic reaction channels: Na(3s) + HD and Na(3p) + HD; the latter can only occur via nonadiabatic transition. The dynamics results show that the product branching of the H + NaD reaction at collision energies ranging from 20 to 80 cm^-1 is controlled by stereodynamics. The Na(3s) and Na(3p) reaction channels occur through collinear collision and side-on collision, respectively. When the collision energy is lower than 20 cm^-1, the resonance-mediated reaction mechanism is dominant in both the Na(3s) and Na(3p) reaction channels.

Mechanism and Kinetics of Methane Combustion. Part II: Potential Energy Surface for Hydrogen-Abstraction Reaction of CH4 + O(3P)

Methane combustion plays an important role in various fields such as combustion chemistry and atmospheric chemistry of the stratosphere. Highly accurate study of its initial reaction remains a key challenge. Here, through extensive studies with a state-of-the-art ab initio and neural network method, we present a potential energy surface of the O(³P) + CH₄ → OH + CH₃ reaction on the ground state 1³A and the first excited state 2³A. In this work, the energies of 10 167 points covering all important regions are obtained with state-averaged complete active space self-consistent field calculations and then fitted using the Levenberg-Marquardt algorithm with a root-mean-square error of 0.391 and 0.442 kcal/mol for the 1³A and 2³A states, respectively. This study explores the characteristics of the radical van der Waals (VdW) complex and reveals a detailed mechanism of the methane combustion initial reaction. Within the scope of this mechanism, this surface gives a fairly accurate description of the regions around the saddle point, conical intersection, and vdW wells in the entrance for efficient computational simulations. As a theoretical study on a prototypical polyatomic reaction, it is hopeful that this work will modify our understanding of the primary process in hydrocarbon combustion.