Numerical convergence of the Sinc discrete variable representation for solving molecular vibrational states with a conical intersection in adiabatic representation

Reactant-Product Decoupling Technique Using the Intermediate Coordinate Method

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

Numerical Equivalence of Diabatic and Adiabatic Representations in Diatomic Molecules

The (time-independent) Schrödinger equation for atomistic systems is solved by using the adiabatic potential energy curves (PECs) and the associated adiabatic approximation. In cases where interactions between electronic states become important, the associated nonadiabatic effects are taken into account via derivative couplings (DDRs), also known as nonadiabatic couplings (NACs). For diatomic molecules, the corresponding PECs in the adiabatic representation are characterized by avoided crossings. The alternative to the adiabatic approach is the diabatic representation obtained via a unitary transformation of the adiabatic states by minimizing the DDRs. For diatomics, the diabatic representation has zero DDR and nondiagonal diabatic couplings ensue. The two representations are fully equivalent and so should be the rovibronic energies and wave functions, which result from the solution of the corresponding Schrödinger equations. We demonstrate (for the first time) the numerical equivalence between the adiabatic and diabatic rovibronic calculations of diatomic molecules using the ab initio curves of yttrium oxide (YO) and carbon monohydride (CH) as examples of two-state systems, where YO is characterized by a strong NAC, while CH has a strong diabatic coupling. Rovibronic energies and wave functions are computed using a new diabatic module implemented in the variational rovibronic code Duo. We show that it is important to include both the diagonal Born-Oppenheimer correction and nondiagonal DDRs. We also show that the convergence of the vibronic energy calculations can strongly depend on the representation of nuclear motion used and that no one representation is best in all cases.

Dynamics studies of the H + HBr reaction: Based on a new potential energy surface

The initial state specific quantum wave packet dynamics studies of the H + HBr (v₀ = 0, j₀ = 0-2) reaction were performed using a new global potential energy surface (PES) of the ground state of the BrH₂ system for the collision energy ranging from 0.01 to 2.0 eV. The PES was constructed using the permutation invariant polynomial neural network method based on approximately 63 000 ab initio points, which were calculated by the multireference configuration interaction method with AVTZ and AVQZ basis sets. To improve the accuracy of the PES, Davidson's correction and spin-orbit coupling effects were considered in the ab initio calculation and the basis set was extrapolated to complete basis set limit. The new PES was compared with the previous ones and also the available experimental data, which suggests that the new PES is more accurate. The state-to-state quantum wave packet dynamics was carried out using the reactant-coordinate based approach. The reaction probabilities, integral and differential cross sections, rovibrational state distributions of product and rate constants, etc., were compared with the available theoretical and experimental studies. In general, the present work is in better agreement with the available experimental data. The quantum dynamics studies suggest that the rotational excitation of HBr has little effect on the reaction.