Quantum reactive scattering with a transmission-free absorbing potential

Coupled-trajectory surface hopping with sign consistency

The framework of exact factorization (XF) has inspired a series of trajectory-based nonadiabatic dynamics methods by introducing different approximations. Recently, the coupled-trajectory surface hopping (CTSH) method has been proposed to combine the key advantages of the coupled-trajectory mixed quantum-classical method based on XF and the fewest switches surface hopping. We here present a novel variant of CTSH, namely, sign-consistent CTSH (SC-CTSH), which considers proper trajectory clustering to reconstruct the nuclear density distribution and the consistency between wave function and active states to introduce decoherence. Using the exact quantum solutions as references, the high performance of SC-CTSH is benchmarked in the widely studied scattering models and compared with other related XF-based methods. Due to the incorporation of new trajectory clustering and sign consistency algorithms, SC-CTSH obtains more accurate quantum momentum and decoherence during the nonadiabatic dynamics, which makes the combination of XF and surface hopping more consistent and reliable. This study further highlights the significance of internal consistency between wave function and active states, which is important in the further development of mixed quantum-classical dynamics methods.

Detailed Complementary Consistency: Wave Function Tells Particle How to Hop, Particle Tells Wave Function How to Collapse

In mixed quantum-classical dynamics, the quantum subsystem can have both wave function and particle-like descriptions. However, they may yield inconsistent results for the expectation value of the same physical quantity. We here propose a novel detailed complementary consistency (DCC) method based on the principle of detailed internal consistency. Namely, the wave function along each trajectory tells the particle how to hop, while the particle tells the wave function how to collapse based on active states in the trajectory ensemble. As benchmarked in a diverse array of representative models with localized nonadiabatic couplings, DCC not only achieves fully consistent results (i.e., identical populations calculated based on wave functions and active states) but also closely reproduces the exact quantum results. Due to the high performance, our new DCC method has great potential to give a consistent and accurate mixed quantum-classical description of general nonadiabatic dynamics after further development.

Higher-Order Split Operator Schemes for Solving Tetratomic Reactions Using the Time-Dependent Wave Packet Method

In this work, using the time-dependent quantum wave packet method, quite a few typical higher-order split operators (HOSOs) were for the first time applied to calculate the tetratomic reactive scattering processes in the hyperspherical coordinate. It was found that the HOSOs were hardly efficient for a tetratomic reaction calculation, unlike those for a triatomic reactive scattering calculation. We proposed an efficient HOSO with a force gradient (denoted as 2G1 in the main text) for efficiently and accurately calculating a tetratomic reaction using the quantum wave packet method. Several typical tetratomic reactions, such as H2 + OH, HF + OH, and H2 + OH+, are calculated for demonstrating the effectiveness of the proposed 2G1 in terms of (product state-resolved) reaction probability and inelastic probability, by comparing with the performance of the previously reported various HOSOs. We suggest that the 2G1 propagator could be applied to efficiently calculate a general tetratomic reaction.

Consistent Construction of the Density Matrix from Surface Hopping Trajectories

Proper construction of the density matrix based on surface hopping trajectories remains a difficult problem. Due to the well-known overcoherence in traditional surface hopping simulations, the electronic wave function cannot be used directly. In this work, we propose a consistent density matrix construction method, which takes the advantage of occupation of active states to rescale the coherence calculated by wave functions and ensures the intrinsic consistency of the density matrix. This new trajectory analysis method can be used for both Tully's fewest switches surface hopping (FSSH) and our recently proposed branching corrected surface hopping (BCSH). As benchmarked in both one- and two-dimensional standard scattering models, the new approach combined with BCSH trajectories achieves highly accurate time-dependent spatial distributions of adiabatic populations and coherence compared to exact quantum results.

Surface Hopping with Reliable Wave Function by Introducing Auxiliary Wave Packets to Trajectory Branching

It is well-known that the widely utilized fewest switches surface hopping method suffers from the severe overcoherence problem, and thus adiabatic populations calculated by wave functions are generally inferior to those based on active states. More importantly, to achieve a complete description of nonadiabatic dynamics, the density matrix is essential. In this paper, we present an auxiliary branching corrected surface hopping (A-BCSH) method that introduces auxiliary wave packets (WPs) on the adiabatic potential energy surfaces for trajectory branching. Both rapid and gradual separation of WP components on different surfaces are characterized, and thus the correct decoherence time along each trajectory is captured. As demonstrated in the three standard Tully models, A-BCSH exhibits excellent internal consistency. Namely, close adiabatic populations are obtained based on both wave functions and active states. In particular, A-BCSH successfully obtains a reliable time-dependent spatial distribution of the density matrix, which relies only on electronic wave functions. Due to its high performance, our A-BCSH method provides a new and highly promising perspective on further development of more consistent surface hopping with reliable wave function.

What is Missing in the Mean Field Description of Spatial Distribution of Population? Important Role of Auxiliary Wave Packets in Trajectory Branching

When the traditional Ehrenfest mean field approach is employed to simulate nonadiabatic dynamics, an effective wave packet (WP) on the average potential energy surface (PES) is utilized to describe the nuclear motion. In the fully quantum picture, however, the WP components on different adiabatic PESs gradually separate in space because they evolve under different velocities and forces. Due to trajectory branching of the WP components, proper decoherence needs to be taken into account, and the spatial distribution of population cannot be described by a single effective WP. Here, we propose an auxiliary branching corrected mean field (A-BCMF) method, where trajectories of auxiliary WPs on adiabatic PESs are introduced. As benchmarked in the three standard Tully models, A-BCMF not only gives correct channel populations but also captures an accurate time-dependent spatial distribution of population. Thereby, we reveal the important role of auxiliary WPs in solving intrinsic problems of the widely used mean field description of nonadiabatic dynamics.

Quantum Transport in Large-Scale Patterned Nitrogen-Doped Graphene

It has recently been demonstrated how the nitrogen dopant concentration in graphene can be controlled spatially on the nano-meter scale using a molecular mask. This technique may be used to create ballistic electron optics-like structures of high/low doping regions; for example, to focus electron beams, harnessing the quantum wave nature of the electronic propagation. Here, we employ large-scale Greens function transport calculations based on a tight-binding approach. We first benchmark different tight-binding models of nitrogen in graphene with parameters based on density functional theory (DFT) and the virtual crystal approximation (VCA). Then, we study theoretically how the random distribution within the masked regions and the discreteness of the nitrogen scattering centers impact the transport behavior of sharp n-p and n-n' interfaces formed by different, realistic nitrogen concentrations. We investigate how constrictions for the current can be realized by patterned high/low doping regions with experimentally feasible nitrogen concentrations. The constrictions can guide the electronic current, while the quantized conductance is significantly washed out due to the nitrogen scattering. The implications for device design is that a p-n junction with nitrogen corrugation should still be viable for current focusing. Furthermore, a guiding channel with less nitrogen in the conducting canal preserves more features of quantized conductance and, therefore, its low-noise regime.

Learning Decoherence Time Formulas for Surface Hopping from Quantum Dynamics

Surface hopping simulations have achieved great success in many different fields, but their reliability has long been limited by the overcoherence problem. We here present a general machine learning assisted approach to identify optimal decoherence time formulas for surface hopping using exact quantum dynamics as references. In order to avoid computationally expensive force calculations, we use the nuclear kinetic energy and the adiabatic energy difference to iteratively generate the descriptor space. Through multilayer screening of the candidate descriptors and discrete optimization of the relevant parameters, we obtain new energy-based decoherence time formulas. As benchmarked in thousands of diverse multilevel systems and six standard scattering models, surface hopping with our new decoherence time formulas nearly reproduces the exact quantum dynamics while maintaining high efficiency. Thereby, our approach provides a promising avenue for systematically improving the accuracy of surface hopping simulations in complex systems from quantum dynamics data.

A Mixed Deterministic-Stochastic Algorithm of the Branching Corrected Mean Field Method for Nonadiabatic Dynamics

We present a new algorithm of the branching corrected mean field (BCMF) method for nonadiabatic dynamics [J. Xu and L. Wang, J. Phys. Chem. Lett. 11, 8283 (2020)], which combines the key advantages of the two existed algorithms, i.e., the deterministic BCMF algorithm based on weights of trajectory branches (BCMF-w) and the stochastic BCMF algorithm with random collapse of the electronic wavefunction (BCMF-s). The resulting mixed deterministic-stochastic BCMF algorithm (BCMF-ws) is benchmarked in a series of standard scattering problems with potential wells on the excited-state surfaces, which are common in realistic systems. In all investigated cases, BCMF-ws holds the same high accuracy while the computational time is reduced about two orders of magnitude compared to the original BCMF-w and BCMF-s algorithms, thus promising for nonadiabatic dynamics simulations of general systems.

Time-dependentab initioapproaches for high-harmonic generation spectroscopy

High-harmonic generation (HHG) is a nonlinear physical process used for the production of ultrashort pulses in XUV region, which are then used for investigating ultrafast phenomena in time-resolved spectroscopies. Moreover, HHG signal itself encodes information on electronic structure and dynamics of the target, possibly coupled to the nuclear degrees of freedom. Investigating HHG signal leads to HHG spectroscopy, which is applied to atoms, molecules, solids and recently also to liquids. Analysing the number of generated harmonics, their intensity and shape gives a detailed insight of, e.g., ionisation and recombination channels occurring in the strong-field dynamics. A number of valuable theoretical models has been developed over the years to explain and interpret HHG features, with the three-step model being the most known one. Originally, these models neglect the complexity of the propagating electronic wavefunction, by only using an approximated formulation of ground and continuum states. Many effects unravelled by HHG spectroscopy are instead due to electron correlation effects, quantum interference, and Rydberg-state contributions, which are all properly captured by anab initioelectronic-structure approach. In this review we have collected recent advances in modelling HHG by means ofab initiotime-dependent approaches relying on the propagation of the time-dependent Schrödinger equation (or derived equations) in presence of a very intense electromagnetic field. We limit ourselves to gas-phase atomic and molecular targets, and to solids. We focus on the various levels of theory employed for describing the electronic structure of the target, coupled with strong-field dynamics and ionisation approaches, and on the basis used to represent electronic states. Selected applications and perspectives for future developments are also given.

Branching and phase corrected surface hopping: A benchmark of nonadiabatic dynamics in multilevel systems

Since the seminal work of Tully [J. Chem. Phys. 93, 1061 (1990)], two-level scattering models have been extensively adopted as the standard benchmark systems to assess the performance of different trajectory surface hopping methods for nonadiabatic dynamics simulations. Here, we extend the branching and phase corrections to multilevel systems and combine them with both the traditional fewest switches surface hopping (FSSH) and its variant global flux surface hopping (GFSH) algorithms. To get a comprehensive evaluation of the proposed methods, we construct a series of more challenging and diverse three-level and four-level scattering models and use exact quantum solutions as references. Encouragingly, both FSSH and GFSH with the branching and phase corrections produce excellent and nearly identical results in all investigated systems, indicating that the new surface hopping methods are robust to describe multilevel problems and the reliability is insensitive to the definition of self-consistent hopping probabilities in the adiabatic representation. Furthermore, the branching correction is found to be especially important when dealing with strongly repulsive potential energy surfaces, which are common in realistic systems, thus promising for general applications.

A systematic construction of Gaussian basis sets for the description of laser field ionization and high-harmonic generation

A precise understanding of mechanisms governing the dynamics of electrons in atoms and molecules subjected to intense laser fields has a key importance for the description of attosecond processes such as the high-harmonic generation and ionization. From the theoretical point of view, this is still a challenging task, as new approaches to solve the time-dependent Schrödinger equation with both good accuracy and efficiency are still emerging. Until recently, the purely numerical methods of real-time propagation of the wavefunction using finite grids have been frequently and successfully used to capture the electron dynamics in small one- or two-electron systems. However, as the main focus of attoscience shifts toward many-electron systems, such techniques are no longer effective and need to be replaced by more approximate but computationally efficient ones. In this paper, we explore the increasingly popular method of expanding the wavefunction of the examined system into a linear combination of atomic orbitals and present a novel systematic scheme for constructing an optimal Gaussian basis set suitable for the description of excited and continuum atomic or molecular states. We analyze the performance of the proposed basis sets by carrying out a series of time-dependent configuration interaction calculations for the hydrogen atom in fields of intensity varying from 5 × 10¹³ W/cm² to 5 × 10¹⁴ W/cm². We also compare the results with the data obtained using Gaussian basis sets proposed previously by other authors.

Quantum dynamics of H2 in a carbon nanotube: Separation of time scales and resonance enhanced tunneling

Quantum confinement effects are known to affect the behavior of molecules adsorbed in nanostructured materials. In order to study these effects on the transport of a single molecule through a nanotube, we present a quantum dynamics study on the diffusion of H₂ in a narrow (8,0) carbon nanotube in the low pressure limit. Transmission coefficients for the elementary step of the transport process are calculated using the flux correlation function approach and diffusion rates are obtained using the single hopping model. The different time scales associated with the motion in the confined coordinates and the motion along the nanotube's axis are utilized to develop an efficient and numerically exact approach, in which a diabatic basis describing the fast motion in the confined coordinate is employed. Furthermore, an adiabatic approximation separating the dynamics of confined and unbound coordinates is studied. The results obtained within the adiabatic approximation agree almost perfectly with the numerically exact ones. The approaches allow us to accurately study the system's dynamics on the picosecond time scale and resolve resonance structures present in the transmission coefficients. Resonance enhanced tunneling is found to be the dominant transport mechanism at low energies. Comparison with results obtained using transition state theory shows that tunneling significantly increases the diffusion rate at T < 120 K.

Perpendicular Emission, Dichroism, and Energy Dependence in Angle-Resolved Photoemission: The Importance of The Final State

Angle-resolved photoemission spectroscopy has been developed to a very high accuracy. However, effects that depend sensitively on the state of the emitted photoelectron were so far hard to compute for real molecules. We here show that the real-time propagation approach to time-dependent density functional theory allows us to obtain final-state effects consistently from first principles and with an accuracy that allows for the interpretation of experimental data. In a combined theoretical and experimental study we demonstrate that the approach captures three hallmark effects that are beyond the final-state plane-wave approximation: emission perpendicular to the light polarization, circular dichroism in the photoelectron angular distribution, and a pronounced energy dependence of the photoemission intensity.

Theoretical scheme for simultaneously observing forward-backward photoelectron holography

Photoelectron angular momentum distribution of He⁺ driven by a few-cycle laser is investigated numerically. We simultaneously observe two dominant interference patterns with one shot of lasers by solving the 3D time-dependent Schrodinger equation. Analysis of a semiclassical model identifies these two interference patterns as two types of photoelectron holography. The interference pattern with P_z>0 is a type of forward rescattering holography, which comes from the interference between direct (reference) and rescattered (signal) forward electrons ionized in the same quarter-cycle. The interference pattern with P_z<0 is a type of backward rescattering holography, which comes from the interference between a direct electron ionized in the third quarter-cycle and rescattered backward electron ionized in the first quarter-cycle. Moreover, we propose a method to distinguish this backward rescattering holography and intracycle interference patterns of direct electrons.

Time-dependent density functional theory for open systems with a positivity-preserving decomposition scheme for environment spectral functions

Understanding electronic dynamics on material surfaces is fundamentally important for applications including nanoelectronics, inhomogeneous catalysis, and photovoltaics. Practical approaches based on time-dependent density functional theory for open systems have been developed to characterize the dissipative dynamics of electrons in bulk materials. The accuracy and reliability of such approaches depend critically on how the electronic structure and memory effects of surrounding material environment are accounted for. In this work, we develop a novel squared-Lorentzian decomposition scheme, which preserves the positive semi-definiteness of the environment spectral matrix. The resulting electronic dynamics is guaranteed to be both accurate and convergent even in the long-time limit. The long-time stability of electronic dynamics simulation is thus greatly improved within the current decomposition scheme. The validity and usefulness of our new approach are exemplified via two prototypical model systems: quasi-one-dimensional atomic chains and two-dimensional bilayer graphene.

Spin-orbit and rotational couplings in radiative association of C(3P) and N(4S) atoms

The role of spin-orbit and rotational couplings in radiative association of C((3)P) and N((4)S) atoms is investigated. Couplings among doublet electronic states of the CN radical are considered, giving rise to a 6-state model of the process. The solution of the dynamical problem is based on the L(2) method, where a complex absorbing potential is added to the Hamiltonian operator in order to treat continuum and bound levels in the same manner. Comparison of the energy-dependent rate coefficients calculated with and without spin-orbit and rotational couplings shows that the couplings have a strong effect on the resonance structure and low-energy baseline of the rate coefficient.

Study of the H+O2 reaction by means of quantum mechanical and statistical approaches: the dynamics on two different potential energy surfaces

The possible existence of a complex-forming pathway for the H+O(2) reaction has been investigated by means of both quantum mechanical and statistical techniques. Reaction probabilities, integral cross sections, and differential cross sections have been obtained with a statistical quantum method and the mean potential phase space theory. The statistical predictions are compared to exact results calculated by means of time dependent wave packet methods and a previously reported time independent exact quantum mechanical approach using the double many-body expansion (DMBE IV) potential energy surface (PES) [Pastrana et al., J. Phys. Chem. 94, 8073 (1990)] and the recently developed surface (denoted XXZLG) by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. The statistical approaches are found to reproduce only some of the exact total reaction probabilities for low total angular momenta obtained with the DMBE IV PES and some of the cross sections calculated at energy values close to the reaction threshold for the XXZLG surface. Serious discrepancies with the exact integral cross sections at higher energy put into question the possible statistical nature of the title reaction. However, at a collision energy of 1.6 eV, statistical rotationally resolved cross sections managed to reproduce the experimental cross sections for the H+O(2)(v=0,j=1)-->OH(v(')=1,j('))+O process reasonably well.

State-to-state reaction probabilities for the H+O2(v,j)→O+OH(v′,j′) reaction on three potential energy surfaces

We report state-to-state and total reaction probabilities for J=0 and total reaction probabilities for J=2 and 4 for the title reaction, both for ground-state and initially rovibrationally excited reactants. The results for three different potential energy surfaces are compared and contrasted. The potential energy surfaces employed are the DMBE IV surface by Pastrana et al. [J. Phys. Chem. 94, 8073 (1990)], the surface by Troe and Ushakov (TU) [J. Chem. Phys. 115, 3621 (2001)], and the new XXZLG ab initio surface by Xu et al. [J. Chem. Phys. 122, 244305 (2005)]. Our results show that the total reaction probabilities from both the TU and XXZLG surfaces are much smaller in magnitude for collision energies above 1.2 eV compared to the DMBE IV surface. The three surfaces also show different behavior with regards to the effect of initial state excitation. The reactivity is increased on the XXZLG and the TU surfaces and decreased on the DMBE IV surface. Vibrational and rotational product state distributions for the XXZLG and the DMBE IV surface show different behaviors for both types of distributions. Our results show that for energies above 1.25 eV the dynamics on the DMBE IV surface are not statistical. However, there is also evidence that the dynamics on the XXZLG surface are not purely statistical for energies above the onset of the first excited product vibrational state v'=1. The magnitude of the total reaction probability is decreased for J>0 for the DMBE IV and the XXZLG surfaces for ground-state reactants. However, for initially rovibrationally excited reactants, the total reaction probability does not decrease as expected for both surfaces. As a result the total cross section averaged over all Boltzmann accessible rotational states may well be larger than the cross section reported in the literature for j=1.

Multichannel scattering calculations using absorbing potentials and mapped grids

The authors investigate the use of absorbing potentials and discrete variable representation grid methods in multichannel time-independent scattering calculations. An exactly solvable, coupled-two-channel problem involving square-well potentials is used to assess the quality of numerical results. Special emphasis is given to the description of scattering resonances and near-threshold regions. Numerical treatment of close vicinities of thresholds requires the introduction of nonequidistant grids through a mapping procedure of the scattering coordinate.

Time dependent wave packet and statistical calculations on the H + O(2) reaction

The H + O(2)--> OH + O reaction has been theoretically investigated by means of an exact time dependent wave packet method and two statistical approaches: a recently developed statistical quantum model and phase-space theory. The exhaustive analysis of reaction probabilities at a zero total angular momentum would, in principle, reveal the existence of a complex-forming mechanism at low collision energies (E(c) = 1.15 eV), whereas deviations from a statistical behaviour at higher energies may be interpreted as the onset of a direct abstraction pathway which favours the production of highly excited rotational states of the OH fragment in its ground vibrational state. The good description by statistical means of previously measured product rotational distributions and excitation functions seems to support such an interpretation. However the statistical predictions clearly overestimate both existing and present exact quantum mechanical reaction probabilities and total cross sections, thereby precluding to conclude definitely the statistical nature of the collision. The exact time dependent method yields values of the integral cross sections in agreement with results by Goldfield and Meijer, and below the experimental findings.

The role of asymptotic vibrational states in H3+

Calculations are discussed which characterize all the vibrational bound states of the H3+ and D2H+ molecular ions using a realistic ab initio potential energy surface. Graphical analysis and calculation of rotational constants show that both ions support a series of atom-diatom-like long-range states: asymptotic vibrational states. The role of these states in the H3+ system and other molecules is discussed. The vibrational calculations are extended above dissociation where the resulting (Feshbach) resonances are shown to be too short-lived to be of importance for the H3+ photodissociation spectrum.

Quantum study of Eley-Rideal reaction and collision induced desorption of hydrogen atoms on a graphite surface. I. H-chemisorbed case

Collision induced (CI) processes involving hydrogen atoms on a graphite surface are studied quantum mechanically within the rigid, flat surface approximation, using a time-dependent wave packet method. The Eley-Rideal (ER) reaction and collision induced desorption (CID) cross sections are obtained with the help of two propagations which use different sets of coordinates, a "product" and a "reagent" set. Several adsorbate-substrate initial states of the target H atom in the chemisorption well are considered, and CI processes are studied over a wide range of projectile energy. Results show that (i) the Eley-Rideal reaction is the major reactive outcome and (ii) CID cross sections do not exceed 4 A2 and present dynamic thresholds for low values of the target vibrational quantum number. ER cross sections show oscillations at high energies which cannot be reproduced by classical and quasiclassical trajectory calculations. They are related to the vibrational excitation of the reaction products, which is a rather steep decreasing function of the collision energy. This behavior causes a selective population of the low-lying vibrational states and allows the quantization of the product molecular states to manifest itself in a collisional observable. A peak structure in the CID cross section is also observed and is assigned to the selective population of metastable states of the transient molecular hydrogen.