Wave Packet Approach to Adiabatic and Nonadiabatic Dynamics of Cold Inelastic Scatterings

Due to the extremely large de Broglie wavelength of cold molecules, cold inelastic scattering is always characterized by the time-independent close-coupling (TICC) method. However, the TICC method is difficult to apply to collisions of large molecular systems. Here, we present a new strategy for characterizing cold inelastic scattering using wave packet (WP) method. In order to deal with the long de Broglie wavelength of cold molecules, the total wave function is divided into interaction, asymptotic and long-range regions (IALR). The three regions use different numbers of ro-vibrational basis functions, especially the long-range region, which uses only one function corresponding to the initial ro-vibrational state. Thus, a very large grid range can be used to characterize long de Broglie wavelengths in scattering coordinates. Due to its better numerical scaling law, the IALR-WP method has great potential in studying the inelastic scatterings of larger collision systems at cold and ultracold regimes.


Introduction
The study of molecular collisions at low temperatures has attracted great interest due to the special properties of cold molecules [1][2][3][4][5]. Cold atoms and molecules can be prepared in a single quantum state, and their motion is controlled by quantum effects. Cold atoms and molecules have extremely long de Broglie wavelengths and the long-range potentials play an important role in collisions. In addition, only one or a few partial waves contribute to collisions, allowing for resonance in cross sections to have better resolution. These properties of cold molecules are of great significance for understanding the nature of molecular collisions, and thus promote the study of molecular collisions at cold and ultracold temperatures [6][7][8][9].
Molecular inelastic collisions are common in the interstellar medium [10], combustion [11] and chemical lasers [12]. The study of energy transfer process in molecular inelastic collisions is very important for understanding many physical phenomena. The distribution of intramolecular states in molecular inelastic scattering is affected by many factors at low temperature, such as long-range van der Waals interaction, shape resonance, Feschbach resonance, initial orientation of molecules, etc. Additionally, only the full quantum dynamics calculation can reveal the influence of these factors on the molecular inelastic scattering. Up to now, most of the theoretical studies on inelastic scattering of cold/ultracold molecules have used the time-independent close coupling (TICC) method [13][14][15][16]. For TICC method, the numerical costs scale as N [3], where N is the number of basis set used in calculations. Therefore, TICC method is currently only applied to collision systems with fewer open channels. At present, full-dimensional TICC calculations can be performed for most of the inelastic collisions of triatomic systems and some tetratomic systems. In order to improve the practicability of TICC method, many approximations have been proposed to reduce the number of channels in TICC calculation. The most widely used approximation is the rigid rotor model [17], which ignores the vibrational dimension of the colliding molecules and

Theory
In this section, the IALR-WP method for solving bimolecular inelastic scattering is introduced. In the IALR-WP method, total grid space is decomposed into interaction region (R ∈ [R 1 , R 2 ]), asymptotic region (R ∈ [R 2 , R 3 ]) and long-range region (R ∈ [R 3 , R 4 ]) by the radial distance between colliding particles, as shown in Figure 1.
The total time dependent wave packet in body-fixed (BF) frame can be expanded as Ψ J Mp (R, r, t) = ∑ n,v,j,K where R is radial distance between colliding atoms/molecules, and r represents internal degree of molecule. u n vj (R) is radial basis function, and ψ vj is ro-vibrational wave function of molecule. For A + BC collisions, v and j are vibrational and rotational quantum states of BC molecule, respectively. For AB + CD collisions, v denotes (v 1 , v 2 ) and j denotes (j 1 , j 2 ), where v 1 (v 2 ) and j 1 (j 2 ) are vibrational and rotational quantum states of AB (CD) molecule, respectively.  (1) where R is radial distance between colliding atoms/molecules, and r represents internal degree of molecule. un vj (R) is radial basis function, and ψvj is ro-vibrational wave function of molecule. For A + BC collisions, v and j are vibrational and rotational quantum states of BC molecule, respectively. For AB + CD collisions, v denotes (v1, v2) and j denotes (j1, j2), where v1 (v2) and j1 (j2) are vibrational and rotational quantum states of AB (CD) molecule, respectively.
In the IALR-WP method, different numbers of ψvj function are used for three regions. Three sets of wave functions Ψ I , Ψ II and Ψ III are defined in the calculation, corresponding to the interaction, asymptotic and long-range regions, respectively. The long-range region is mainly used to describe the initial wave packet and the elastic scattered wave packet, so the Ψ III is defined in the range of R ∈ [R1, R4], and only a single ψvj with v = v0, j = j0 is included. The asymptotic region is mainly used to describe the inelastic scatterings, so the Ψ II is defined in the range of R ∈ [R1, R3], and a few ψvj functions with v ∈ [0, vAsy] and j ∈ [0, jAsy] are included. The interaction region is mainly used to describe the reactive scatterings, so the Ψ I is defined in the range of R ∈ [R1, R2], and a large amount of ψvj functions with v ∈ [0, vInt] and j ∈ [0, jInt] are included. It can be found that there are overlapping regions between the three wave functions. During wave packet propagation, three wave functions exchange information with each other through the overlapping regions. The rovibrational state dependent radial basis function is defined as I  I  I   II  II  II   III  III  III sin( ) , with ICSs show that rotational quenching ICSs are much larger that of vibrational quenching. The total vibrational quenching ICS in Figure 2a is multiplied by 10 [5]. Three In the IALR-WP method, different numbers of ψ vj function are used for three regions. Three sets of wave functions Ψ I , Ψ II and Ψ III are defined in the calculation, corresponding to the interaction, asymptotic and long-range regions, respectively. The long-range region is mainly used to describe the initial wave packet and the elastic scattered wave packet, so the Ψ III is defined in the range of R ∈ [R 1 , R 4 ], and only a single ψ vj with v = v 0 , j = j 0 is included. The asymptotic region is mainly used to describe the inelastic scatterings, so the Ψ II is defined in the range of R ∈ [R 1 , R 3 ], and a few ψ vj functions with v ∈ [0, v Asy ] and j ∈ [0, j Asy ] are included. The interaction region is mainly used to describe the reactive scatterings, so the Ψ I is defined in the range of R ∈ [R 1 , R 2 ], and a large amount of ψ vj functions with v ∈ [0, v Int ] and j ∈ [0, j Int ] are included. It can be found that there are overlapping regions between the three wave functions. During wave packet propagation, three wave functions exchange information with each other through the overlapping regions. The ro-vibrational state dependent radial basis function is defined as where N R I , N R II and N R III are grid numbers used in Ψ I , Ψ II and Ψ III , respectively. ∆R is the step size.
The initial wave packet is product of Gaussian wave packet G(R) and initial rovibrational state wave function. The G(R) is given as where, µ R is the reduced mass for R coordinate. R 0 and δ define the position and width of the initial wave packet, respectively, and E 0 is the mean collision energy. The initial wave packet is propagated by second-order split-operator method [34], which is described by whereT andV are the kinetic energy operator and reference potential energy. For nonadiabatic calculations, the potential energy operator is not diagonal in the diabatic representation, the wave packet needs to be first transformed to the adiabatic representation when operating the potential energy terms: where E i and T are eigenvalues and eigenvectors of the reference diabatic potential energy matrix.
In the IALR-WP method, different damping functions are used for asymptotic and long-range regions. The ro-vibrational state dependent damping function for the Ψ II has following form: and the damping function for the Ψ III has following form: where, C R II (C R III ) defines the strength of damping function, with R a II (R a III ) and R b II (R b III ) as the corresponding starting and ending points of damping functions. Since reactive scattering is not involved in this work, the damping function is not used for r coordinate in the calculations.
To collect the elastic and inelastic scattering information form WP calculations, a projection plane is set in the asymptotic region. After the propagation of wave packet, the time-independent scattering wavefunction at projection plane ψ + (E;R ∞ ) can be obtained by a time-energy Fourier transform of time-dependent wave packet at the projection plane. The state-to-state scattering S-matrix is obtained by where k f is the momenta in the exit channel, h l f is outgoing Ricatti-Hankel function, χ v f j f is molecular ro-vibrational wavefunction. The R ∞ is position of the projection plane. The α(E) is amplitude of the initial Gaussian wave packet at a specified collision energy, and given by where, k v 0 j 0 is the momenta in the entrance channel, l 0 is the initial orbital angular momentum quantum number. Normally, in the well-established WP method, the projection plane should be set between the tail of the initial wave packet and the starting point of the damping function.
However, in the IALR-WP method, the initial wave packet is always constructed in the long-range region, behind the damping function in the asymptotic region. This means that we cannot use one projection plane to collect both elastic and inelastic information. Since elastic and inelastic scattered waves are treated separately in the IALR-WP method, it allows us to use two different projection planes to collect elastic and inelastic information separately. One projection plane is placed in the asymptotic region to collect information of inelastic scattering, and the other is placed in the long-range region to collect information of elastic scattering, as shown in Figure 1. Since the calculation of inelastic scattering wavefunction is independent from elastic one, and the inelastic scattering projection plane is just set at the asymptotic region, the needed total propagation time is greatly shortened for inelastic scattering compared with elastic scattering.
After propagation for a series of J and l 0 , the state-to-state inelastic cross section for atom-diatom collision is obtained by Additionally, diatom-diatom inelastic scattering cross section is obtained by

Results and Discussion
In this section, in order to test the accuracy of the IALR-WP method in characterizing inelastic scattering at cold temperatures, the dynamics calculations are performed for three typical collision systems. The three collision systems selected in this paper represent three different inelastic collision cases: the adiabatic atom-diatom collision, the nonadiabatic atom-diatom collision and the adiabatic diatom-diatom collision. The adiabatic calculations results are compared with those of the TICC method, and the non-adiabatic calculations results are compared with those of the L-shaped WP method.
First, in order to verify that the IALR-WP method can accurately describe the dynamics of inelastic collisions in adiabatic triatomic systems, the ro-vibrational quenching ICSs of Na(3s) + HD(v = 1, j = 2) collision are calculated at low temperatures and compared with the accurate TICC results in Figure 2. The PES constructed by Hack et al. [35] is used in the calculations. The calculations are performed for J = 0~6 partial waves, and ro-vibrational quenching ICSs are obtained in collision energy range of 0.01 to 1 K. The numerical parameters used in the IALR-WP calculation are listed in Table 1. In the TICC calculations, energy cutoff is 3.0 eV, a total of 10 vibrational basis is used for r direction, and the other parameters are also carefully tested. As seen from Figure 2, the IALR-WP calculated ICSs are in good agreement with TICC results, which demonstrates the accuracy of IALR-WP method for characterizing adiabatic atom-diatom inelastic scatterings. The product ro-vibrational state resolved ICSs show that rotational quenching ICSs are much larger that of vibrational quenching. The total vibrational quenching ICS in Figure 2a is multiplied by 10 [5]. Three resonance peaks are found in the ICS curves. To obtain more insights into the formation of these peaks, the partial wave resolved inelastic scattering ICSs of Na(3s) + HD(v = 1, j = 2) → Na(3s) + HD(v = 1, j = 1) are shown in Figure 2b. As seen from Figure 2b, the first peak locates at 0.05 K, and mainly contributed by J = 0 partial wave. The second peak locates at 0.1 K, and J = 1, 2, 3 and 4 partial waves all contribute, but the contribution of J = 2 is the most obvious. The third peak locates at 0.6 K and it mainly contributed by J = 5 partial wave. of these peaks, the partial wave resolved inelastic scattering ICSs of Na(3s) + HD(v = 1, j = 2) → Na(3s) + HD(v′ = 1, j′ = 1) are shown in Figure 2b. As seen from Figure 2b, the first peak locates at 0.05 K, and mainly contributed by J = 0 partial wave. The second peak locates at 0.1 K, and J = 1, 2, 3 and 4 partial waves all contribute, but the contribution of J = 2 is the most obvious. The third peak locates at 0.6 K and it mainly contributed by J = 5 partial wave.  Next, we verify the accuracy of the IALR-WP method for inelastic scattering in nonadiabatic atom-diatom collisions at low temperatures. When the collision contains more than one electronic state, the diabatic representation is used in the WP calculation to avoid the singularity caused by the non-adiabatic couplings [36][37][38]. When the non-adiabatic coupling vanishes at long-range region, the IALR-WP method could be readily extended to non-adiabatic calculations as described in our previous work [32]. In such circumstances, the diabatic PESs approach to adiabatic PESs as R increases, therefore the initial wave packet could be constructed in diabatic representation. Additionally, the elastic and inelastic scattered wave packets can also be separated in diabatic representation by using IALR grid method. By using this strategy, the non-adiabatic IALR-WP calculations are performed for Na(3p) + H 2 (v = 0, 1) collisions in this work.
For the Na(3p) + H 2 collision, the processes of electronic to vibrational and rotational energy transfer (Na(3p) + H 2 → Na(3s) + H 2 (v , j )) and chemical reaction (Na(3p) + H 2 → NaH + H) through an excited state complex have been widely studied by experimentally and theoretically [39][40][41][42][43][44]. The electronic state of collision system changes during the collision for both electronic quenching and reactive collisions. To study the non-adiabatic collisions between Na(3p) and H 2 , Hack et al. [35] reported an analytical diabatic potential energy matrix for the two lowest electronic states of NaH 2 system. The long-range forces between Na(3p) and H 2 are included and non-diagonal matrix elements at long-range region are damped to zero in this non-adiabatic PESs, therefore it is a good candidate for testing the accuracy and ability of IALR-WP method for cold non-adiabatic collisions. For the Na(3p) + H 2 → NaH + H reaction, it is an endothermic reaction and the endothermicity is about 0.47 eV, therefore the reaction channel is not considered in the current calculations. To better describe the non-adiabatic processes in the collision between Na(3p) and H 2 , a 3D plot of diabatic PESs as functions of R and r at C 2v geometry is shown in Figure 3a, and contour plot is shown in Figure 3b. There is an exciplex near the R = 4.2 a 0 and r = 1.4 a 0 , and the r of the exciplex is close to the equilibrium bond length of H 2 molecule. The black dash line in Figure 3b represents the crossing of two diabatic PESs. The minimum potential energy of the crossing line is almost close to the Na(3p) + H 2 asymptotic energy, indicating that the electronic quenching should be no barriers. The values of r on the crossing line are greater than that equilibrium bond length of H 2 molecule, thus H 2 should be stretched before electronic quenching occurs. Therefore, a relatively large r grid range is required in the dynamics calculations of the Na(3p) + H 2 collision.  To test the accuracy of IARL-WP method in characterizing electronic quenching sc tering, the total electronic quenching probabilities of Na(3p) + H2(v = 0, 1) collisions calculated at higher collision energies by IARL-WP method, and the results are compa with those of L-shaped WP method. In Figure 4, J = 0 total electronic quenching proba ities of Na(3p) + H2(v = 0, 1) collisions obtained by IALR-WP method and L-shaped W method are compared. The agre To test the accuracy of IARL-WP method in characterizing electronic quenching scattering, the total electronic quenching probabilities of Na(3p) + H 2 (v = 0, 1) collisions are calculated at higher collision energies by IARL-WP method, and the results are compared with those of L-shaped WP method. In Figure 4, J = 0 total electronic quenching probabilities of Na(3p) + H 2 (v = 0, 1) collisions obtained by IALR-WP method and L-shaped WP method are compared. The agreement between two methods demonstrates the accuracy of the IALR-WP method in characterizing electronic quenching scatterings. There are many resonance peaks in the probability curves, which should be caused by the potential well on the excited state PES. In addition, although the vibrational excitation of H 2 molecule opens up more product channels, it reduces the total electronic quenching probability. This is because there are ro-vibrational quenching channels Na(3p) + H 2 (v = 0, j ) in addition to the electronic quenching channel in the Na(3p) + H 2 (v = 1) collision. ro-vibrational quenching channels Na(3p) + H2(v′ = 0, j′) in addition to the electronic quenching channel in the Na(3p) + H2(v = 1) collision. After validating the IALR-WP method in treating non-adiabatic inelastic scatterings, it is employed to compute electronic quenching ICSs of Na(3p) + H2(v = 0, 1) collisions at low temperatures. Only J = 0~8 partial waves contribute to electronic quenching scattering in the collision energy range of 0.1 to 10 K. Since the electronic quenching of Na(3p) + H2 collision can produce molecules with very high ro-vibrational states, a relatively large amount of ro-vibrational basis is used in both the interaction region and the asymptotic region. The main parameters used in the IALR-WP calculations are listed in Table 2. In Figure 5, the electronic quenching ICSs of Na(3p) + H2 collision with two different initial vibrational states v = 0 and v = 1 are shown. It can be seen that the vibrational excitation of H2 molecule reduces the probability of electronic quenching, and the same result is obtained when the collision energy is high. In addition, the vibrational excitation of H2 molecule also has a great influence on the product state distribution. The products of Na(3p) + H2(v = 0) collision are mainly distributed in the v′ = 3 state, while the products of Na(3p) + H2(v = 1) collision are mainly distributed in the v′ = 5 state. To obtain more insights into the product state distribution, the ro-vibrational state distribution of Na(3p) + H2(v = 0, 1) collisions at collision energy of 0.1 K are shown in Figure 6. For the collision with v = 0, the products are mainly distributed in the rotational ground state and the vibrational excited state, which accords with the characteristics of the diabatic potential energy surface discussed above, and also After validating the IALR-WP method in treating non-adiabatic inelastic scatterings, it is employed to compute electronic quenching ICSs of Na(3p) + H 2 (v = 0, 1) collisions at low temperatures. Only J = 0~8 partial waves contribute to electronic quenching scattering in the collision energy range of 0.1 to 10 K. Since the electronic quenching of Na(3p) + H 2 collision can produce molecules with very high ro-vibrational states, a relatively large amount of ro-vibrational basis is used in both the interaction region and the asymptotic region. The main parameters used in the IALR-WP calculations are listed in Table 2. In Figure 5, the electronic quenching ICSs of Na(3p) + H 2 collision with two different initial vibrational states v = 0 and v = 1 are shown. It can be seen that the vibrational excitation of H 2 molecule reduces the probability of electronic quenching, and the same result is obtained when the collision energy is high. In addition, the vibrational excitation of H 2 molecule also has a great influence on the product state distribution. The products of Na(3p) + H 2 (v = 0) collision are mainly distributed in the v = 3 state, while the products of Na(3p) + H 2 (v = 1) collision are mainly distributed in the v = 5 state. To obtain more insights into the product state distribution, the ro-vibrational state distribution of Na(3p) + H 2 (v = 0, 1) collisions at collision energy of 0.1 K are shown in Figure 6. For the collision with v = 0, the products are mainly distributed in the rotational ground state and the vibrational excited state, which accords with the characteristics of the diabatic potential energy surface discussed above, and also indicates that the collision tends to occur in the T-geometry. When H 2 molecule is excited to v = 1 state, the electronic quenching products have a wide rotational state distribution, but mainly distribute at lower rotational states. Table 2. Numerical parameters used in IALR-WP calculations for the Na(3p) + H 2 (v = 0, 1) → Na(3s) + H 2 collisions at low temperatures.    Finally, we extend the IALR-WP method to diatom-diatom inelastic scattering at low temperatures. Since the H 2 + H 2 is the simplest neutral diatom-diatom system, the H 2 + H 2 and its isotopic variants are treated as a benchmark system to studying the quantum dynamics of diatom-diatom scatterings and testing new numerical calculation methods [45][46][47][48][49]. Here, we use the HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) collision as an example to illustrate the accuracy and ability of IALR-WP method in characterizing cold diatom-diatom inelastic scattering. The strategy is basically the same as in the case of atom-diatom scattering, except that the ro-vibrational quantum states of molecule are slightly more complicated for diatom-diatom collisions. For the HD(v = 0, j = 1) + H 2 (v = 0, j = 0) collision, Sultanov et al. [50] reported the rotational quenching ICS at low temperatures by using rigid-rotor model. To make a comparison with previous results, four-dimensional (one vibrational basis is used for r 1 and r 2 ) IALR-WP and TICC calculations are performed for HD(v = 0, j = 1) + H 2 (v = 0, j = 0) collision, and BMKP [51] PES is used. Main parameter used in the IALR-WP calculation are listed in Table 3. In the TICC method, E cut is set to 1.0 eV and the other parameters are chosen to ensure the convergence. The partial wave J is calculated from 0 to 5 and the converged ICSs are obtained at collision energy ranges from 0.01 to 20 K. A comparison of rotational quenching ICSs of HD(v = 0, j = 1) + H 2 (v = 0, j = 0) collision calculated by 4D TICC method, 4D IALR-WP method and rigid-rotor TICC method are shown in Figure 7a. The agreement of ICSs in Figure 7a demonstrates the accuracy of the results of 4D TICC and 4D IALR-WP in this work. To obtain more insights in to the contribution of partial waves for rotational quenching ICS, the partial wave resolved inelastic scattering ICSs are shown in Figure 7b. There are three peaks found in the total ICS of HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) collision. The first peak is located at 0.6 K, and mainly contributed by J = 1, 3 partial waves. The second peak is located at 6 K, and mainly contributed by J = 2, 4 partial waves. The third peak is located at 10.5 K and it mainly contributed by J = 3, 5 partial waves. The agreement between IALR-WP method and TICC method in partial wave resolved ICSs further demonstrate the accuracy of IALR-WP method and its ability in characterizing cold diatom-diatom inelastic scatterings.
Molecules 2022, 27, x FOR PEER REVIEW 7 of 8 Table 2. Numerical parameters located at 6 K, and mainly contributed by J = 2, 4 partial waves. The third peak is located at 10.5 K and it mainly contributed by J = 3, 5 partial waves. The agreement between IALR-WP method and TICC method in partial wave resolved ICSs further demonstrate the accuracy of IALR-WP method and its ability in characterizing cold diatom-diatom inelastic scatterings.  In some cases, the rigid rotor approximation produces large errors, so full-dimensional quantum calculations are required to describe molecular collisions. In Figure 8, the 4D and 6D calculated rotational quenching ICSs of the HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) collision are shown. In the 6D calculation, four vibrational basis functions and three vibrational functions are used for HD and H 2 molecule, respectively. Only small differences are found in the different dimensional calculated ICSs, indicating that the rotational quenching process of HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) can be described well by the rigid-rotor approximation. In addition, the results obtained by full dimensional IALR-WP method are in good agreement with those calculated by TICC method. Table 3. Numerical parameters used in IALR-WP calculations for the HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) collision at low temperatures. H2(v = 0, j = 0) collision are shown. In the 6D calculation, four vibrational basis functions and three vibrational functions are used for HD and H2 molecule, respectively. Only small differences are found in the different dimensional calculated ICSs, indicating that the rotational quenching process of HD(v = 0, j = 1) + H2(v = 0, j = 0) → HD(v = 0, j = 0) + H2(v = 0, j = 0) can be described well by the rigid-rotor approximation. In addition, the results obtained by full dimensional IALR-WP method are in good agreement with those calculated by TICC method.

Conclusions
The quantum dynamics of molecular collisions at low temperatures are always characterized by the TICC method. Since the TICC method is difficult to extend to complex collision systems due to its steep numerical scaling law, the quantum dynamics calculations of full dimensional molecular collisions have been limited to three or four atom systems. The recent studies show that the IALR-WP method has great potential in the study of quantum dynamics of complex systems at cold and ultracold regimes. In this work, the IALR-WP method is extended to study the molecular inelastic scattering at low temperatures, and the accuracy of the method is verified in three different inelastic scattering cases.
Firstly, to demonstrate the accuracy of IALR-WP method in characterizing the adiabatic atom-diatom inelastic scattering, the ro-vibrational quenching ICSs of Na(3s) + HD(v = 1, j = 2) collision at cold temperatures is calculated, and compared with the results of TICC method. The calculation results show that the vibrational quenching ICS is much lower than that of rotational quenching ICS, and HD molecule is mainly distributed in the v = 1, j = 1 state. In addition, three resonance peaks are found in the ro-vibrational quenching ICSs, and the contribution of J partial waves to the resonances are analyzed.
Secondly, to demonstrate the accuracy of IALR-WP method in characterizing the nonadiabatic atom-diatom inelastic scatterings, the J = 0 electronic quenching probabilities of Na(3p) + H 2 (v = 0, 1) collisions are calculated by IALR-WP method, and the accuracy of the method is verified by comparison with the results of L-shaped WP method. The electronic quenching ICSs at low temperatures show that the vibrational excitation of H 2 reduces the electronic quenching probability and has a great influence on the product state distribution.
Finally, the IARL-WP method is applied to HD(v = 0, j = 1) + H 2 (v = 0, j = 0) → HD(v = 0, j = 0) + H 2 (v = 0, j = 0) inelastic scattering to verify the accuracy of the IALR-WP method in characterizing diatom-diatom inelastic scatterings at low temperatures. Three resonance peaks are found on the rotational quenching ICS and the contribution of different J partial waves to those resonances are analyzed. In addition, only small differences are found in the 4D and 6D calculated ICSs, indicating that the HD(v = 0, j = 1) + H 2 (v = 0, j = 0) collision can be described well by the rigid-rotor approximation.
In recent years, with the development of experimental techniques, collision of highly excited state molecules has attracted extensive attention. Excited molecules, including rotational, vibrational and electronic excitation, can be prepared using several pumping techniques, even under ultracold conditions. The collision of excited state molecules, especially electronic excited state molecules, will produce product molecules with very high ro-vibrational states, which makes the TICC method very difficult. However, the IALR-WP method proposed in this work can describe such collision process effectively. Due to the better computational scaling in the IALR-WP method, this method is expected to be of great use in studying inelastic scattering involving heavy atoms and multi-electronic states, such as inelastic scattering of KRb + Rb, RbCs + Cs at low temperatures.