Rate constants for the H+ + H2 reaction from 5 K to 3000 K with a statistical quantum method

Ortho-Para Conversion for H+ + H2 Collision in Low Temperature: A Fully Close-Coupled Time-Dependent Wave Packet Study

The collision-induced rate coefficients of ortho-para conversion for the H+ + H2 reaction provide accurate information to probe the lifetime of cold environments in interstellar media. Rotationally resolved reaction probabilities are calculated at the low collision energy regime (0 < Ecol ≤ 0.3 eV) by employing the coupled three-dimensional (3D) time-dependent wave packet (TDWP) formalism in hyperspherical coordinates on a recently constructed ab initio ground adiabatic potential energy surface of H3+ [J. Chem. Phys. 2014, 141, 204306] for the process H+ + H2 (v = 0, j = 0-5) → H+ + H2 (v' = 0, j'). Cross-sections are then computed from the converged reaction probabilities as a function of total angular momentum (J) over the same energy regime and subsequently employed to obtain the rate constants for the ortho-to-para (O-P) and para-to-ortho (P-O) conversions and their ratio. The ratio of ortho-para conversion shows (a) appropriate convergence with the inclusion of a higher number of initial rotational states as well as a reasonable agreement with the results from a quantum statistical method and (b) a peak at lower temperature that could be due to the available collision energy for transitions involving lower rotational states (j = 0 → j' = 1).

Simulating chemical reaction dynamics on quantum computer

The electronic energies of molecules have been successfully evaluated on quantum computers. However, more attention is paid to the dynamics simulation of molecules in practical applications. Based on the variational quantum eigensolver (VQE) algorithm, Fedorov et al. proposed a correlated sampling (CS) method and demonstrated the vibrational dynamics of H2 molecules [J. Chem. Phys. 154, 164103 (2021)]. In this study, we have developed a quantum approach by extending the CS method based on the VQE algorithm (labeled eCS-VQE) for simulating chemical reaction dynamics. First, the CS method is extended to the three-dimensional cases for calculation of first-order energy gradients, and then, it is further generalized to calculate the second-order gradients of energies. By calculating atomic forces and vibrational frequencies for H2, LiH, H+ + H2, and Cl- + CH3Cl systems, we have seen that the approach has achieved the CCSD level of accuracy. Thus, we have simulated dynamics processes for two typical chemical reactions, hydrogen exchange and chlorine substitution, and obtained high-precision reaction dynamics trajectories consistent with the classical methods. Our eCS-VQE approach, as measurement expectations and ground-state wave functions can be reused, is less demanding in quantum computing resources and is, therefore, a feasible means for the dynamics simulation of chemical reactions on the current noisy intermediate-scale quantum-era quantum devices.

Benchmarking an improved statistical adiabatic channel model for competing inelastic and reactive processes

Inelastic collisions and elementary chemical reactions proceeding through the formation and subsequent decay of an intermediate collision complex, with an associated deep well on the potential energy surface, pose a challenge for accurate fully quantum mechanical approaches, such as the close-coupling method. In this study, we report on the theoretical prediction of temperature-dependent state-to-state rate coefficients for these complex-mode processes, using a statistical quantum method. This statistical adiabatic channel model is benchmarked by a direct comparison using accurate rate coefficients from the literature for a number of systems (H₂ + H⁺, HD + H⁺, SH⁺ + H, and CH⁺ + H) of interest in astrochemistry and astrophysics. For all of the systems considered, an error of less than factor 2 was found, at least for the dominant transitions and at low temperatures, which is sufficiently accurate for applications in the above mentioned disciplines.