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.

Study on the UV FEL single-shot damage threshold of an Au thin film

The damage threshold of an Au-coated flat mirror, one of the reflective optics installed on the FEL-2 beamline of the Dalian Coherent Light Source, China, upon far-UV free-electron laser irradiation is evaluated. The surface of the coating is characterized by profilometer and optical microscope. A theoretical approach of the phenomenon is also presented, by application of conventional single-pulse damage threshold calculations, a one-dimensional thermal diffusion model, as well as finite-element analysis with ANSYS.

Technologies for investigating single-molecule chemical reactions

Single molecules, the smallest independently stable units in the material world, serve as the fundamental building blocks of matter. Among different branches of single-molecule sciences, single-molecule chemical reactions, by revealing the behavior and properties of individual molecules at the molecular scale, are particularly attractive because they can advance the understanding of chemical reaction mechanisms and help to address key scientific problems in broad fields such as physics, chemistry, biology and materials science. This review provides a timely, comprehensive overview of single-molecule chemical reactions based on various technical platforms such as scanning probe microscopy, single-molecule junction, single-molecule nanostructure, single-molecule fluorescence detection and crossed molecular beam. We present multidimensional analyses of single-molecule chemical reactions, offering new perspectives for research in different areas, such as photocatalysis/electrocatalysis, organic reactions, surface reactions and biological reactions. Finally, we discuss the opportunities and challenges in this thriving field of single-molecule chemical reactions.

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.

Quantum interference in atom-exchange reactions

Chemical reactions, where bonds break and form, are highly dynamic quantum processes. A fundamental question is whether coherence can be preserved in chemical reactions and then harnessed to generate entangled products. Here we investigated this question by studying the 2KRb → K2 + Rb2 reaction at 500 nK, focusing on the nuclear spin degrees of freedom. We prepared the initial nuclear spins in KRb in an entangled state by lowering the magnetic field to where the spin-spin interaction dominates and characterized the preserved coherence in nuclear spin wavefunction after the reaction. We observed an interference pattern that is consistent with full coherence at the end of the reaction, suggesting that entanglement prepared within the reactants could be redistributed through the atom-exchange process.

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.

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.

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.

Multiple Dissociation Pathways in HNCO Decomposition Governed by Potential Energy Surface Topography

The exquisite features of molecular photochemistry are key to any complete understanding of the chemical processes governed by potential energy surfaces (PESs). It is well established that multiple dissociation pathways relate to nonadiabatic transitions between multiple coupled PESs. However, little detail is known about how the single PES determines reaction outcomes. Here we perform detailed experiments on HNCO photodissociation, acquiring the state-specific correlations of the NH (a1Δ) and CO (X1Σ+) products. The experiments reveal a trimodal CO rotational distribution. Dynamics simulations based on a full-dimensional machine-learning-based PES of HNCO unveil three dissociation pathways exclusively occurring on the S1 excited electronic state. One pathway, following the minimum energy path (MEP) via the transition state, contributes to mild rotational excitation in CO, while the other two pathways deviating substantially from the MEP account for relatively cold and hot CO rotational state populations. These peculiar dynamics are unambiguously governed by the S1 state PES topography, i.e., a narrow acceptance cone in the vicinity of the transition state region. The dynamical picture shown in this work will serve as a textbook example illustrating the importance of the PES topography in molecular photochemistry.

Topological Metamaterials

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

Geometric Phase Effect in Attosecond Stimulated X-ray Raman Spectroscopy

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

Control of Raman Scattering Quantum Interference Pathways in Graphene

Graphene is an ideal platform to study the coherence of quantum interference pathways by tuning doping or laser excitation energy. The latter produces a Raman excitation profile that provides direct insight into the lifetimes of intermediate electronic excitations and, therefore, on quantum interference, which has so far remained elusive. Here, we control the Raman scattering pathways by tuning the laser excitation energy in graphene doped up to 1.05 eV. The Raman excitation profile of the G mode indicates its position and full width at half-maximum are linearly dependent on doping. Doping-enhanced electron-electron interactions dominate the lifetimes of Raman scattering pathways and reduce Raman interference. This will provide guidance for engineering quantum pathways for doped graphene, nanotubes, and topological insulators.

Globally Accurate Gaussian Process Potential Energy Surface and Quantum Dynamics Studies on the Li(2S) + Na2 → LiNa + Na Reaction at Low Collision Energies

The LiNa2 reactive system has recently received great attention in the experimental study of ultracold chemical reactions, but the corresponding theoretical calculations have not been carried out. Here, we report the first globally accurate ground-state LiNa2 potential energy surface (PES) using a Gaussian process model based on only 1776 actively selected high-level ab initio training points. The constructed PES had high precision and strong generalization capability. On the new PES, the quantum dynamics calculations on the Li(2S) + Na2(v = 0, j = 0) → LiNa + Na reaction were carried out in the 0.001-0.01 eV collision energy range using an improved time-dependent wave packet method. The calculated results indicate that this reaction is dominated by a complex-forming mechanism at low collision energies. The presented dynamics data provide guidance for experimental research, and the newly constructed PES could be further used for ultracold reaction dynamics calculations on this reactive system.

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.

Ultrafast light-induced dynamics in the microsolvated biomolecular indole chromophore with water

Interactions between proteins and their solvent environment can be studied in a bottom-up approach using hydrogen-bonded chromophore-solvent clusters. The ultrafast dynamics following UV-light-induced electronic excitation of the chromophores, potential radiation damage, and their dependence on solvation are important open questions. The microsolvation effect is challenging to study due to the inherent mix of the produced gas-phase aggregates. We use the electrostatic deflector to spatially separate different molecular species in combination with pump-probe velocity-map-imaging experiments. We demonstrate that this powerful experimental approach reveals intimate details of the UV-induced dynamics in the near-UV-absorbing prototypical biomolecular indole-water system. We determine the time-dependent appearance of the different reaction products and disentangle the occurring ultrafast processes. This approach ensures that the reactants are well-known and that detailed characteristics of the specific reaction products are accessible - paving the way for the complete chemical-reactivity experiment.

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.

Energy, Momentum, and Angular Momentum Transfer between Electrons and Nuclei

The recently developed exact factorization approach condenses all electronic effects on the nuclear subsystem into scalar and vector potentials that appear in an effective time dependent Schrödinger equation. Starting from this equation, we derive subsystem Ehrenfest identities characterizing the energy, momentum, and angular momentum transfer between electrons and nuclei. An effective electromagnetic force operator induced by the electromagnetic field corresponding to the effective scalar and vector potentials appears in all three identities. The effective magnetic field has two components that can be identified with the Berry curvature calculated with (a) different Cartesian coordinates of the same nucleus and (b) arbitrary Cartesian coordinates of two different nuclei. (a) has a classical interpretation as the induced magnetic field felt by the nucleus, while (b) has no classical analog. Subsystem Ehrenfest identities are ideally suited for quantifying energy transfer in electron-phonon systems. With two explicit examples we demonstrate the usefulness of the new identities.

Low-temperature reaction dynamics of paramagnetic species in the gas phase

Radicals are abundant in a range of important gas-phase environments. They are prevalent in the atmosphere, in interstellar space, and in combustion processes. As such, understanding how radicals react is essential for the development of accurate models of the complex chemistry occurring in these gas-phase environments. By controlling the properties of the colliding reactants, we can also gain insights into how radical reactions occur on a fundamental level. Recent years have seen remarkable advances in the breadth of experimental methods successfully applied to the study of reaction dynamics involving paramagnetic species-from improvements to the well-known crossed molecular beams approach to newer techniques involving magnetically guided and decelerated beams. Coupled with ever-improving theoretical methods, quantum features are being observed and interesting insights into reaction dynamics are being uncovered in an increasingly diverse range of systems. In this highlight article, we explore some of the exciting recent developments in the study of chemical dynamics involving paramagnetic species. We focus on low-energy reactive collisions involving neutral radical species, where the reaction parameters are controlled. We conclude by identifying some of the limitations of current methods and exploring possible new directions for the field.

Control of reactive collisions by quantum interference

In this study, we achieved magnetic control of reactive scattering in an ultracold mixture of ²³Na atoms and ²³Na⁶Li molecules. In most molecular collisions, particles react or are lost near short range with unity probability, leading to the so-called universal rate. By contrast, the Na + NaLi system was shown to have only ~4% loss probability in a fully spin-polarized state. By controlling the phase of the scattering wave function via a Feshbach resonance, we modified the loss rate by more than a factor of 100, from far below to far above the universal limit. The results are explained in analogy with an optical Fabry-Perot resonator by interference of reflections at short and long range. Our work demonstrates quantum control of chemistry by magnetic fields with the full dynamic range predicted by our models.

Bimolecular Chemistry in the Ultracold Regime

Advances in atomic, molecular, and optical physics techniques allowed the cooling of simple molecules down to the ultracold regime ([Formula: see text]1 mK) and opened opportunities to study chemical reactions with unprecedented levels of control. This review covers recent developments in studying bimolecular chemistry at ultralow temperatures. We begin with a brief overview of methods for producing, manipulating, and detecting ultracold molecules. We then survey experimental works that exploit the controllability of ultracold molecules to probe and modify their long-range interactions. Further combining the use of physical chemistry techniques such as mass spectrometry and ion imaging significantly improved the detection of ultracold reactions and enabled explorations of their dynamics in the short range. We discuss a series of studies on the reaction KRb + KRb → K₂ + Rb₂ initiated below 1 [Formula: see text]K, including the direct observation of a long-lived complex, the demonstration of product rotational state control via conserved nuclear spins, and a test of the statistical model using the complete quantum state distribution of the products. Expected final online publication date for the Annual Review of Physical Chemistry, Volume 73 is April 2022. Please see http://www.annualreviews.org/page/journal/pubdates for revised estimates.

Quantum mechanical double slit for molecular scattering

[Figure: see text].

Nearside-Farside Analysis of the Angular Scattering for the State-to-State H + HD → H2 + D Reaction: Nonzero Helicities

We theoretically analyze the differential cross sections (DCSs) for the state-to-state reaction, H + HD(v_i = 0, j_i = 0, m_i = 0) → H₂(v_f = 0, j_f = 1,2,3, m_f = 1,..,j_f) + D, over the whole range of scattering angles, where v, j, and m are the vibrational, rotational, and helicity quantum numbers for the initial and final states. The analysis extends and complements previous calculations for the same state-to-state reaction, which had j_f = 0,1,2,3 and m_f = 0, as reported by XiahouC.; ConnorJ. N. L.Phys. Chem. Chem. Phys.2021, 23, 13349–1336934096934. Motivation comes from the state-of-the-art experiments and simulations of Yuan et al.Nature Chem.2018, 10, 653–65829686377 who have measured, for the first time, fast oscillations in the small-angle region of the degeneracy-averaged DCSs for j_f = 1 and 3 as well as slow oscillations in the large-angle region. We start with the partial wave series (PWS) for the scattering amplitude expanded in a basis set of reduced rotation matrix elements. Then our main theoretical tools are two variants of Nearside-Farside (NF) theory applied to six transitions: (1) We apply unrestricted, restricted, and restrictedΔ NF decompositions to the PWS including resummations. The restricted and restrictedΔ NF DCSs correctly go to zero in the forward and backward directions when m_f > 0, unlike the unrestricted NF DCSs, which incorrectly go to infinity. We also exploit the Local Angular Momentum theory to provide additional insights into the reaction dynamics. Properties of reduced rotation matrix elements of the second kind play an important role in the NF analysis, together with their caustics. (2) We apply an approximate N theory at intermediate and large angles, namely, the Semiclassical Optical Model of Herschbach. We show there are two different reaction mechanisms. The fast oscillations at small angles (sometimes called Fraunhofer diffraction/oscillations) are an NF interference effect. In contrast, the slow oscillations at intermediate and large angles are an N effect, which arise from a direct scattering, and are a “distorted mirror image” mechanism. We also compare these results with the experimental data.

Atomistic dynamics of elimination and nucleophilic substitution disentangled for the F- + CH3CH2Cl reaction

Chemical reaction dynamics are studied to monitor and understand the concerted motion of several atoms while they rearrange from reactants to products. When the number of atoms involved increases, the number of pathways, transition states and product channels also increases and rapidly presents a challenge to experiment and theory. Here we disentangle the dynamics of the competition between bimolecular nucleophilic substitution (S_N2) and base-induced elimination (E2) in the polyatomic reaction F^- + CH₃CH₂Cl. We find quantitative agreement for the energy- and angle-differential reactive scattering cross-sections between ion-imaging experiments and quasi-classical trajectory simulations on a 21-dimensional potential energy hypersurface. The anti-E2 pathway is most important, but the S_N2 pathway becomes more relevant as the collision energy is increased. In both cases the reaction is dominated by direct dynamics. Our study presents atomic-level dynamics of a major benchmark reaction in physical organic chemistry, thereby pushing the number of atoms for detailed reaction dynamics studies to a size that allows applications in many areas of complex chemical networks and environments.

Inelastic Scattering of H Atoms from Surfaces

We have developed an instrument that uses photolysis of hydrogen halides to produce nearly monoenergetic hydrogen atom beams and Rydberg atom tagging to obtain accurate angle-resolved time-of-flight distributions of atoms scattered from surfaces. The surfaces are prepared under strict ultrahigh vacuum conditions. Data from these experiments can provide excellent benchmarks for theory, from which it is possible to obtain an atomic scale understanding of the underlying dynamical processes governing H atom adsorption. In this way, the mechanism of adsorption on metals is revealed, showing a penetration-resurfacing mechanism that relies on electronic excitation of the metal by the H atom to succeed. Contrasting this, when H atoms collide at graphene surfaces, the dynamics of bond formation involving at least four carbon atoms govern adsorption. Future perspectives of H atom scattering from surfaces are also outlined.

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.

Mechanistic Insights into Ultracold Chemical Reactions under the Control of the Geometric Phase

Ultracold chemical reactions involve collision temperatures approaching absolute zero, and for molecular systems that exhibit a barrierless and exoergic reaction path significant reactivity can occur. In addition, many molecules contain a conical intersection, and the associated geometric phase has been shown to significantly alter the outcome of ultracold reactions. Here we report a quantum dynamics study for the ultracold O + OH → H + O₂ reaction. An analysis of the scattering wave functions reveals explicitly the nature of the quantum interference between the direct and looping reaction pathways around the conical intersection and thus illustrates how the reaction proceeds under the control of the geometric phase for the first time. The wave function analysis should generalize to other ultracold reactions that contain a conical intersection. Our findings indicate that quantum control techniques such as an optical lattice trap or the initial state orientation may be effective in controlling the reactivity.

Interaction-Asymptotic Region Decomposition Method for an Insertion Reaction: Application to the S(1D) + H2 Reaction

With adjusting principal axes hyperspherical (APH) coordinate in the interaction region, and the Jacobi coordinates in the asymptotic regions, an efficient multidomain interaction-asymptotic region decomposition (IARD) method has been developed to solve the "coordinate problem" in a product-state-resolved reactive scattering calculation using the quantum wave packet method. Although the APH coordinate treats with all three channels equally, and is efficient for describing the interaction region for some direct reactions, it is inefficient for describing the insertion-type reaction due to the singularity problem, such as the S(¹D) + H₂ reaction. To deal with this issue, in this work, the channel-dependent Delves hyperspherical (DH) coordinate is proposed to describe the interaction region using the IARD method. The proposed DH-IARD method was applied to calculate the product-state-resolved reaction probabilities of the H + HD reaction, and the differential and integral cross sections of the typical insertion reaction S(¹D) + H₂. It is found that the new DH-IARD method is much more efficient than the previous APH-IARD method for dealing with insertion reactions. The partial wave resonance structures were observed in the integral cross section. It is found that at a low collision energy, the position of the initial wave packet has to be put far away. Otherwise, the partial wave resonance structures could not be correctly reproduced due to the reef well arising with a large total angular momentum J.

Controlling the Spin-Orbit Branching Fraction in Molecular Collisions

The collision geometry, that is, the relative orientation of reactants before interaction, can have a large effect on how a collision or reaction proceeds. Certain geometries may prevent access to a given product channel, while others might enhance it. In this Letter, we demonstrate how the initial orientation of NO molecules relative to approaching Ar atoms determines the branching between the spin-orbit changing and the spin-orbit conserving rotational product channels. We use a recently developed quantum treatment to calculate differential and integral branching fractions, at any arbitrary orientation, from theoretical and experimental data points. Our results show that a substantial degree of control over the final spin-orbit state of the scattering products can be achieved by tuning the initial collision geometry.

Electric field-catalyzed single-molecule Diels-Alder reaction dynamics

Precise time trajectories and detailed reaction pathways of the Diels-Alder reaction were directly observed using accurate single-molecule detection on an in situ label-free single-molecule electrical detection platform. This study demonstrates the well-accepted concerted mechanism and clarifies the role of charge transfer complexes with endo or exo configurations on the reaction path. An unprecedented stepwise pathway was verified at high temperatures in a high-voltage electric field. Experiments and theoretical results revealed an electric field-catalyzed mechanism that shows the presence of a zwitterionic intermediate with one bond formation and variation of concerted and stepwise reactions by the strength of the electric field, thus establishing a previously unidentified approach for mechanistic control by electric field catalysis.

Glories, hidden rainbows and nearside-farside interference effects in the angular scattering of the state-to-state H + HD → H2 + D reaction

Yuan et al. [Nat. Chem., 2018, 10, 653] have reported state-of-the-art measurements of differential cross sections (DCSs) for the H + HD → H2 + D reaction, measuring for the first time fast oscillations in the small-angle forward region of the DCSs. We theoretically analyse the angular scattering dynamics in order to quantitatively understand the physical content of structure in the DCSs. We study the H + HD(vi = 0, ji = 0, mi = 0) → H2(vf = 0, jf = 0,1,2,3, mf = 0) + D reaction for the whole range of scattering angles from θR = 0° to θR = 180°, where v, j, m are the vibrational, rotational and helicity quantum numbers respectively for the initial and final states. The restriction to mf = 0 arises because states with mf ≠ 0 have DCSs that are identically zero in the forward (θR = 0°) and backward (θR = 180°) directions. We use accurate quantum scattering matrix elements computed by Yuan et al. at a translational energy of 1.35 eV for the BKMP2 potential energy surface. The following theoretical techniques are employed to analyse the DCSs: (a) full and nearside-farside (NF) partial wave series (PWS) and local angular momentum theory, including resummations of the full PWS up to third order. We also use window representations of the scattering matrix, which give rise to truncated PWS, (b) six asymptotic (semiclassical) small-angle glory theories and four N rainbow theories, (c) we introduce "CoroGlo" tests, which let us distinguish between glory and corona scattering at small angles for Legendre PWS, (d) the semiclassical optical model (SOM) of Herschbach is employed to understand structure in the DCSs at intermediate and large angles. Our conclusions are: (a) the small-angle peaks in the DCSs arise mainly from glory scattering. For the 000 → 020 transition, there is also a contribution from a broad, or hidden, N rainbow, (b) at larger angles, the fast oscillations in the DCSs arise from NF interference, (c) the N scattering in the fast oscillation region contains a hidden rainbow for the 000, 020, 030 cases. For the 000 → 020 transition, the rainbow extends up to θR ≈ 60°; for the 000 and 030 cases, the angular ranges containing a N rainbow are smaller, (d) at intermediate and backward angles, the slowly varying DCSs, which merge into slow oscillations, are explained by the SOM. Physically it shows this structure in a DCS arises from direct scattering and is a distorted mirror image of the corresponding probability versus total angular momentum quantum number plot.

A fundamental invariant-neural network representation of quasi-diabatic Hamiltonians for the two lowest states of H3

The FI-NN approach is capable of representing highly accurate diabatic PESs with particular and complicated symmetry problems.

Time-dependent quantum mechanical wave packet dynamics

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

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

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

Advances and New Challenges to Bimolecular Reaction Dynamics Theory

Dynamics of bimolecular reactions in the gas phase are of foundational importance in combustion, atmospheric chemistry, interstellar chemistry, and plasma chemistry. These collision-induced chemical transformations are a sensitive probe of the underlying potential energy surface(s). Despite tremendous progress in past decades, our understanding is still not complete. In this Perspective, we survey the recent advances in theoretical characterization of bimolecular reaction dynamics, stimulated by new experimental observations, and identify key new challenges.

Stereodynamics of rotationally inelastic scattering in cold He + HD collisions

Stereodynamics of cold collisions has become a fertile ground for sensitive probe of molecular collisions and control of the collision outcome. A benchmark system for stereodynamic control of rotational transition is He + HD. This system was recently probed experimentally by Perreault et al. by examining quenching from j = 2 to j' = 0 state in the v = 1 vibrational manifold of HD. Here, through explicit quantum scattering calculations on a highly accurate ab initio interaction potential for He + H₂, we reveal how a combination of two shape resonances arising from l = 1 and l = 2 partial waves controls the stereodynamic outcome rather than a single l = 2 partial wave attributed in the experiment. Furthermore, for collision energies below 0.5 cm^-1, it is shown that stereodynamic preference for the integral cross section follows a simple universal trend.

How interference reveals geometric phase

Quantum phase effects are probed at energies below the H + HD reaction conical intersection