Quantum Wave Packet Study of the H+ + D2 Reaction on Diabatic Potential Energy Surfaces

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.

The impact of non-adiabatic effects on reaction dynamics: a study based on the adiabatic and non-adiabatic potential energy surfaces of CaH2. Phys Chem Chem Phys 2023;25:22744-22754. [PMID: 37605513 DOI: 10.1039/d3cp02416d] [Citation(s) in RCA: 0] [Impact Index Per Article: 0]

The two-state non-adiabatic potential energy matrices of the CaH2+ system are calculated via a diabatization approach by using a neural network model. Subsequently, the adiabatic and non-adiabatic potential energy surfaces (PESs) are constructed based on these non-adiabatic potential energy matrices. Furthermore, based on the adiabatic and non-adiabatic PESs, the Ca+(4s2S) + H2(X1Σ+g) → H(2S) + CaH+(X1Σ+) reaction is studied using the time-dependent wave packet method. Comparative analysis of the experimental and theoretical integral reaction cross-sections (ICSs) indicates that the maximum deviation between the results obtained from the adiabatic PES and the corresponding experimental value is 12.7 bohr2; in contrast, the maximum discrepancy between the theoretical result derived from the non-adiabatic PES and the experimental value is merely 0.42 bohr2. The potential well along the reaction path acts as a 'filter', selectively guiding intermediates with longer lifetimes in the potential well back to the reactant channel. This phenomenon indicates that the non-adiabatic effects significantly influence the reaction dynamics of the CaH2+ system.

Charge Transfer Processes for H + H2+ Reaction Employing Coupled 3D Wavepacket Approach on Beyond Born-Oppenheimer Based Ab Initio Constructed Diabatic Potential Energy Surfaces

The dynamics of the H + H₂⁺ reaction has been analyzed from the electronically first excited state of diabatic potential energy surfaces constructed by employing the Beyond Born-Oppenheimer theory [J. Chem. Phys. 2014, 141, 204306]. We have employed the coupled 3D time-dependent wavepacket formalism in hyperspherical coordinates for multisurface reactive scattering problems. To be specific, the charge transfer processes have been investigated extensively by calculating state-to-state as well as total reaction probabilities and integral cross sections, when the reaction process is initiated from the first excited electronic state (2¹A'). We have depicted the convergence profiles of reaction probabilities for the competing charge transfer processes, namely, reactive charge transfer (RCT) and nonreactive charge transfer (NRCT) processes for different total energies with respect to total angular momentum, J. Total and state-to-state integral cross sections are calculated as a function of total energy for the initial rovibrational state, namely, v = 0, j = 0 level of H₂⁺ (²Σ_g⁺) molecule and are compared with previous theoretical calculations. Finally, we have calculated temperature-dependent rate constants using our presently evaluated cross sections and compared their average with the experimentally measured one.

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.

Angular momentum-scattering angle quantum correlation: a generalized deflection function

A natural generalization of the classical deflection function, the functional dependence of the deflection angle on the angular momentum (or the impact parameter), is the joint probability density function of these two quantities, revealing the correlation between them. It provides, at a glance, detailed information about the reaction mechanisms and how changes in the impact parameter affect the product angular distribution. It is also useful to predict the presence of quantum phenomena such as interference. However, the classical angular momentum-scattering angle correlation function has a limited use whenever quantum effects become important. Rigorously speaking, there is not a quantum equivalent of the classical joint distribution, as the differential cross section depends on the coherences between the different values of J caused by the cross terms in the expansion of partial waves. In this article, we present a simple method to calculate a quantum analog of this correlation, a generalized deflection function that can shed light onto the reaction mechanism using just quantum mechanical results. Our results show that there is a very good agreement between the quantum and classical correlation functions as long as quantum effects are not all relevant. When this is not the case, it will also be shown that the quantum correlation function is most useful to observe the extent of quantum effects such as interference among different reaction mechanisms.

QUASICLASSICAL TRAJECTORY STUDY OF THE STEREODYNAMICS FOR THE REACTION D+ + H2 (ν=0, j = 0) → HD + H+

The vector correlations between products and reactants for the reactive noncharge transfer ion–molecule collisions D + + H 2 (ν=0, j = 0)→ HD + H + have been determined by means of the quasiclassical trajectory method on the ground state in the KBNN potential energy surface (J Chem Phys116:654, 2002) at collision energies of 0.224, 0.524, 0.824, and 1.024 eV. The calculated differential cross section (DCS) results indicate that the lifetime of the complex [Formula: see text] in the deep well on the ground PES becomes shorter as collision energy increases. The existence of long-lived complex leads to a weak product rotational polarization at a low collision energy of 0.224 eV. However, the product rotational angular momentum j′ aligns and orients preferentially along the positive direction of the y-axis at a high collision energy of 1.024 eV. The distribution of P(θr, ϕr) indicates that the product molecules are preferentially polarized perpendicular to the scattering plane and that the reaction is dominated by an in-plane mechanism.

THEORETICAL STUDIES OF STEREODYNAMICS FOR THE H+ + H2 (ν = 0–3, j = 0) → H2 + H+ REACTION

The stereodynamics calculation was carried out for the title reaction by quasiclassical trajectory method on the ground surface of KBNN potential energy surface. The vector correlations are determined at initial ground state and vibrational excitation of the reagent H 2. The results show that the rotational polarization is affected lightly by collision energy and strongly by reagent excitation for title reaction. The rotational alignments are almost isotropic at several collision energies on initial ground state of the reagent H 2, which means that the product rotational angular momentum is weakly polarized (or no polarized). Nevertheless, the polarization of product rotational angular momentum is enhanced remarkably at the vibrational excitations of the reagent H 2 in collision energy of 0.524 eV.

Influence of nuclear exchange on nonadiabatic electron processes in H(+)+H2 collisions

H(+)+H(2) collisions are studied by means of a semiclassical approach that explicitly accounts for nuclear rearrangement channels in nonadiabatic electron processes. A set of classical trajectories is used to describe the nuclear motion, while the electronic degrees of freedom are treated quantum mechanically in terms of a three-state expansion of the collision wavefunction. We describe electron capture and vibrational excitation, which can also involve nuclear exchange and dissociation, in the E = 2-1000 eV impact energy range. We compare dynamical results obtained with two parametrizations of the potential energy surface of H(3)(+) ground electronic state. Total cross sections for E > 10 eV agree with previous results using a vibronic close-coupling expansion, and with experimental data for E < 10 eV. Additionally, some prototypical features of both nuclear and electron dynamics at low E are discussed.

Coriolis coupling and nonadiabaticity in chemical reaction dynamics

The nonadiabatic quantum dynamics and Coriolis coupling effect in chemical reaction have been reviewed, with emphasis on recent progress in using the time-dependent wave packet approach to study the Coriolis coupling and nonadiabatic effects, which was done by K. L. Han and his group. Several typical chemical reactions, for example, H+D(2), F+H(2)/D(2)/HD, D(+)+H(2), O+H(2), and He+H(2)(+), have been discussed. One can find that there is a significant role of Coriolis coupling in reaction dynamics for the ion-molecule collisions of D(+)+H(2), Ne+H(2)(+), and He+H(2)(+) in both adiabatic and nonadiabatic context.

Non-Born-Oppenheimer dynamics calculations using the coherent switching with decay of mixing method

A theoretical investigation of the nonadiabatic processes that are involved in the full three-dimensional D+H(2), H+D(2), D(+)+H(2), and H(+)+D(2) reaction systems has been performed using the method of coherence switching with decay of mixing (CSDM) developed by Truhlar and coworkers. The electronic density matrix for each trajectory is fully coherent in the CSDM method, and a switching algorithm different to the original "decay of mixing" method is used to determine the pure state toward which the decoherent force drives the system. By solving the equations for the evolution of the electronic state populations along the semiclassical trajectory, the ensemble can present effective physical insight into nonadiabatic dynamics. The calculations for the D+H(2) and H+D(2) systems is based on the double many body expansion potential energy surface. The potential energy surface constructed by Kamisaka et al. is employed in the calculation of D(+)+H(2) and H(+)+D(2) reactions. The cross sections and the reaction probabilities for the total angular momentum J = 0 are calculated for all of these systems. The calculated results from the CSDM method are in good agreement with exact quantum mechanical calculations and experimental measurements.

Cumulative reaction probabilities and transition state properties: a study of the H+ + H2 and H+ + D2 proton exchange reactions

Cumulative reaction probabilities (CRPs) have been calculated by accurate (converged, close coupling) quantum mechanical (QM), quasiclassical trajectory (QCT), and statistical QCT (SQCT) methods for the H(+) + H(2) and H(+) + D(2) reactions at collision energies up to 1.2 eV and total angular momentum J = 0-4. A marked resonance structure is found in the QM CRP, most especially for the H(3)(+) system and J = 0. When the CRPs are resolved in their ortho and para contributions, a clear steplike structure is found associated with the opening of internal states of reactants and products. The comparison of the QCT results with those of the other methods evinces the occurrence of two transition states, one at the entrance and one at the exit. At low J values, except for the quantal resonance structure and the lack of quantization in the product channel, the agreement between QM and QCT is very good. The SQCT model, that reflects the steplike structure associated with the opening of initial and final states accurately, clearly tends to overestimate the value of the CRP as the collision energy increases. This effect seems more marked for the H(+) + D(2) isotopic variant. For sufficiently high J values, the growth of the centrifugal barrier leads to an increase in the threshold of the CRP. At these high J values the discrepancy between SQCT and QCT becomes larger and is magnified with growing collision energy. The total CRPs calculated with the QCT and SQCT methods allowed the determination of the rate constant for the H(+) + D(2) reaction. It was found that the rate, in agreement with experiment, decreases with temperature as expected for an endothermic reaction. In the range of temperatures between 200 and 500 K the differences between SQCT and QCT rate results are relatively minor. Although exact QM calculations are formidable for an exact determination of the k(T), it can be reliably expected that their value will lie between those given by the dynamical and statistical trajectory methods.

Attosecond-resolution quantum dynamics calculations for atoms and molecules in strong laser fields

A parallel quantum electron and nuclei wave packet computer code, LZH-DICP, has been developed to study laser-atom-molecule interaction in the nonperturbative regime with attosecond resolution. The nonlinear phenomena occurring in that regime can be studied with the code in a rigorous way by numerically solving the time-dependent Schrödinger equation of electrons and nuclei. Time propagation of the wave functions is performed using a split-operator approach, and based on a sine discrete variable representation. Photoelectron spectra for hydrogen and kinetic-energy spectra for molecular hydrogen ion in linearly polarized laser fields are calculated using a flux operator scheme, which testifies to the validity and the high efficiency of LZH-DICP.

On the dynamics of the H+ +D2(v=0,j=0)-->HD+D + reaction: a comparison between theory and experiment

The H+ +D2(v=0,j=0)-->HD+D + reaction has been theoretically investigated by means of a time independent exact quantum mechanical approach, a quantum wave packet calculation within an adiabatic centrifugal sudden approximation, a statistical quantum model, and a quasiclassical trajectory calculation. Besides reaction probabilities as a function of collision energy at different values of the total angular momentum, J, special emphasis has been made at two specific collision energies, 0.1 and 0.524 eV. The occurrence of distinctive dynamical behavior at these two energies is analyzed in some detail. An extensive comparison with previous experimental measurements on the Rydberg H atom with D2 molecules has been carried out at the higher collision energy. In particular, the present theoretical results have been employed to perform simulations of the experimental kinetic energy spectra.

Quantum approaches for the insertion dynamics of the H+ + D2 and D+ + H2 reactive collisions

The H(+)+D(2) and D(+)+H(2) reactive collisions are studied using a recently proposed adiabatic potential energy surface of spectroscopic accuracy. The dynamics is studied using an exact wave packet method on the adiabatic surface at energies below the curve crossing occurring at approximately 1.5 eV above the threshold. It is found that the reaction is very well described by a statistical quantum method for a zero total angular momentum (J) as compared with the exact ones, while for higher J some discrepancies are found. For J >0 different centrifugal sudden approximations are proposed and compared with the exact and statistical quantum treatments. The usual centrifugal sudden approach fails by considering too high reaction barriers and too low reaction probabilities. A new statistically modified centrifugal sudden approach is considered which corrects these two failures to a rather good extent. It is also found that an adiabatic approximation for the helicities provides results in very good agreement with the statistical method, placing the reaction barrier properly. However, both statistical and adiabatic centrifugal treatments overestimate the reaction probabilities. The reaction cross sections thus obtained with the new approaches are in rather good agreement with the exact results. In spite of these deficiencies, the quantum statistical method is well adapted for describing the insertion dynamics, and it is then used to evaluate the differential cross sections.

Spin-orbit effect in the energy pooling reaction O2(aΔ1)+O2(aΔ1)→O2(bΣ1)+O2(XΣ3)

Five-dimensional nonadiabatic quantum dynamics studies have been carried out on two new potential energy surfaces of S(2)((1)A(')) and T(7)((3)A(")) states for the title oxygen molecules collision with coplanar configurations, along with the spin-orbit coupling between them. The ab initio calculations are based on complete active state second-order perturbation theory with the 6-31+G(d) basis set. The calculated spin-orbit induced transition probability as a function of collision energy is found to be very small for this energy pooling reaction. The rate constant obtained from a uniform J-shifting approach is compared with the existing theoretical and experimental data, and the spin-orbit effect is also discussed in this electronic energy-transfer process.

A detailed quantum mechanical and quasiclassical trajectory study on the dynamics of the H++H2→H2+H+ exchange reaction

The H+ + H2 exchange reaction has been studied theoretically by means of a different variety of methods as an exact time independent quantum mechanical, approximate quantum wave packet, statistical quantum, and quasiclassical trajectory approaches. Total and state-to-state reaction probabilities in terms of the collision energy for different values of the total angular momentum obtained with these methods are compared. The dynamics of the reaction is extensively studied at the collision energy of E(coll)=0.44 eV. Integral and differential cross sections and opacity functions at this collision energy have been calculated. In particular, the fairly good description of the exact quantum results provided by the statistical quantum method suggests that the dynamics of the process is governed by an insertion mechanism with the formation of a long-lived collision complex.

Time-Dependent Wave Packet Studies on the Cl + HCl Hydrogen Exchange Reaction

The initiation of the hydrogen exchange reaction Cl((2)P)+HCl --> ClH+Cl((2)P) by excitation of the HCl molecular stretch to v=2 is studied for total angular momentum quantum number J=(1)/(2) and both even and odd parity. The calculations were performed using a time-dependent propagation from an initial quasi-bound state and employed all three relevant potential energy surfaces and the nonadiabatic couplings between them. Coriolis and spin-orbit coupling were also taken into account. The electronic and HCl rotational distributions of the products in both dissociation channels are analyzed, and the results are interpreted using features of the potential energy surfaces.