Chinese Journal of Chemical Physics  2017, Vol. 30 Issue (1): 16-24

#### The article information

Zhi-qiang Zhao, Shu Liu, Dong H. Zhang

Differential Cross Sections for the H+D2O→HD+OD Reaction: a Full Dimensional State-to-State Quantum Dynamics Study
H+D2O→HD+OD反应的微分截面:一种全维态-态量子动力学计算
Chinese Journal of Chemical Physics, 2017, 30(1): 16-24

http://dx.doi.org/10.1063/1674-0068/30/cjcp1608163

### Article history

Received on: August 22, 2016
Accepted on: September 17, 2016
Differential Cross Sections for the H+D2O→HD+OD Reaction: a Full Dimensional State-to-State Quantum Dynamics Study
Zhi-qiang Zhaoa,b, Shu Liua, Dong H. Zhanga
Dated: Received on August 22, 2016; Accepted on September 17, 2016
a. State Key Laboratory of Molecular Reaction Dynamics and Center for Theoretical Computational Chemistry, Dalian Institute of Chemical Physics, Chinese Academy of Sciences, Dalian 116023, China;
b. University of Chinese Academy of Sciences, Beijing 100049, China
*Author to whom correspondence should be addressed. E-mail:liushu1985@dicp.ac.cn; E-mail:zhangdh@dicp.ac.cn
Abstract: The time-dependent wave-packet method was employed to calculate the first full-dimensional state-to-state differential cross sections (DCS) for the title reaction with D2O in the ground and the first symmetric (100) and asymmetric stretching (001) excited states. The calculated DCSs for these three initial states are strongly backward peaked at low collision energies. With the increase of collision energy, these DCSs become increasingly broader with the peak position shifting gradually to a smaller angle, consistent with the fact that the title reaction is a direct reaction via an abstraction mechanism. It is found that the (100) and (001) states not only have roughly the same integral cross sections, but also have essentially identical DCS, which are very close to that for the ground state at the same total energy of reaction. The reaction produces a small fraction of OD in the v=1 state, with the population close to the relative reactivity between the ground and vibrationally excited states, therefore confirming the experimental result of Zare et al. and the local mode picture[J. Phys. Chem. 97, 2204 (1993)]. Unexpectedly, the stretching excitation reduces the rotation excitation of product HD at the same total energy.
Key words: State-to-state    Differential cross section    Time-dependent wave-packet method
I. INTRODUCTION

H+H2O and its isotopically substituted reactions are the prototype for tetra-atomic reactions, in much the same way that the H+H2 reaction served as the prototype for triatomics. Theoretically, because three of the four atoms are hydrogens, the system is an ideal candidate for high quality ab initio calculation of a potential energy surface (PES), as well as for accurate quantum reactive scattering calculations. Schatz and Elgersma fitted the first analytic PES, known as the WDSE PES [1]. After that, a number of new PES have been constructed with increasing level of accuracy for the reaction system [2-5]. Extensive quantum reactive scattering calculations have been performed on these PES to study dynamics of this benchmark four-atom system.

The H+H2O→H2+OH is also considered as the simplest system in which different vibrational modes in the reactant can play an important role in the reaction dynamics. In the past decades, substantial efforts have been devoted to understanding the dependence of reactivity on initial vibrational excitation for the title reaction and its isotopic analogies. Crim, Zare, and co-workers carried out a serial of experimental studies and observed strong mode specific reactivity in the H+H2O/D2O reactions with H2O/D2O in overtone or fundamental vibrational excitations [6-8]. Bond selectivity has also been demonstrated in the H+HOD reaction with a preferential cleavage of the HOD vibrationally excited bond [9-11]. Theoretically, besides the quasiclassical trajectory (QCT) studies of Schatz et al.\ [12-14] there have been a number of quantum dynamics studies [15-19], which are in qualitative agreement with the experimental observations. Recently, Fu et al. calculated the exact coupled-channel (CC) integral cross sections (ICS) for H2O and HOD initially in the first and second stretching excited states and simultaneous excitations of both bending and stretching modes on the YZCL2 PES, and observed strong mode specific and bond selected reactivity as discovered in the experiments [20-21]

Another important question of mode specificity is how the initial vibrational excitation influences the product state distribution. The Crim group observed that the reaction of water molecules excited to the |04〉- overtone state predominantly produces OH(v=0) while reaction from the |13〉- state produces mostly OH(v=1) [7]. Zare group showed that for the D2O asymmetric stretch fundamental states (001) and (011), the OD product is produced with little vibrational excitation (OD(v=l):OD(v=0)≈1:20) [8]. At first glance, such an OD(OH) product vibration distribution might seem surprising. In a normal mode picture, both OD bonds in D2O are equally excited in the symmetrically or asymmetrically excited state. Breaking one OD bond by the incoming H atom should leave the other OD in vibrationally excited state. Zare and coworkers postulated a local mode picture in which the D2O asymmetric stretch state is thought as a linear combination of local mode stretches with the H atom reacting preferentially with the excited OD bond, leaving the other OD group in its vibrational ground state. Until recently, we report the first full-dimensional state-to-state study for the H+H2O→H2+OH reaction with H2O in the ground and the first symmetric and asymmetric stretching excited states [22]. It is found that the (100) and (001) initial excited states both lead to the OH product mostly in its ground vibrational state, confirming the local mode picture of H2O vibration, and revealing the spectator nature of the non-reacting OH bond.

In this work, we report the first full-dimensional (6D) quantum differential cross sections (DCS) for the title reaction with initial nonrotating D2O in (000), (100) and (001) vibrational states to investigate the initial vibrational excitation influences on the product state distribution and DCS, by using the reactant-product decoupling (RPD) scheme based time-dependent wave packet (TDWP) method. Through the development of the TDWP method, it is routine now to calculate fully converged integral cross sections (ICS) without any dynamical approximations for some four-atom reactions [20, 21, 23-27], by only using one set of reactant Jacobi coordinates. However to extract S-matrix elements in a TDWP calculation, one needs to transform the wave function initially prepared in reactant Jacobi coordinates to product Jacobi coordinates, making the calculations much difficult. Up to now, there are basically three state-to-state approaches to deal with the coordinate problem, namely, "reactant coordinates based", "reactant-product decoupling", and "product coordinates based" approach. In 2002 and 2006, Zhang calculated the first state-to-state DCS for the H+H2O→H2+OH reaction from the ground initial state, but in five dimensions with the nonreactive OH bond length fixed [28, 29].

Soon after that, Althorpe and coworkers introduced another coordinate transformation method and applied it to calculate state-to-state reaction probability in full dimensions for the total angular momentum J=0 [30-32]. In 2011, the TDWP method was developed to compute full-dimensional DCSs for four-atom reactions and applied to the prototypical HD+OH→ H2O+D and D2+OH→HOD+D reactions [33-35]. Excellent agreements were achieved for the first time for a four-atom reaction between the theoretical DCS and high-resolution crossed-molecular beam experimental results, declaring that theory is now able to calculate complete dynamical information for some four-atom reactions.

In this work, we present the theory of the full-dimensional state-to-state quantum dynamics treatment to atom-triatom reactions and some computational details.

II. Theory

The RPD method, originally introduced by Peng et al. and Huang et al. [36-37], partitions the full time-dependent wavefunction into a sum of reactant component (Ψr) and all product components (Ψp, p+1, 2, 3, ···) that satisfy the following decoupled equations,

 (1)
 (2)

where -iVp is the negative imaginary potential (absorption potential) employed to prevent the wavefunction Ψr(t) from entering the pth product arrangement. Then, Althorpe and co-workers employed an exact partitioned version of the split-operator propagator [30], which is more efficient than partitioning the Schrödinger equation. Ψr(t) can be propagated in the reactant Jacobi coordinates just as for total reaction probability calculations,

 (3)

Every product wave function, Ψp, can be propagated in its own coordinates as in a normal wave packet propagation except for a source term, ξp, prepared by Eq.(3),

 (4)

In the original RPD approach, the source term is saved and transformed from reactant to product Jacobi coordinates at every propagation time step. To reduce the computational cost, we use multi-step RPD [32, 38] by saving and transforming the source term at every M time steps in the actual calculation,

 (5)

when mod((t+△)-t0, M△)=0,

 (6)

where t0 is the starting point for performing wave function transformation. At other time steps, we carried out the standard split-operator propagation [39]for Ψr(t) and Ψp(t) without the absorption potential, Vp, related terms. In this way, we can cut the computational time for wave function transformation from reactant coordinates to product coordinates by a factor of M.

 FIG. 1 The reagent Jacobi coordinates (R, r1, r2, θa1, θa2, ψ) for the A+BCD atom-triatom arrangement, and product Jacobi coordinates (R', r'1, r'2, θ'1, θ'2, ψ') for the AB+CD diatom-diatom arrangement. △ is the angle between $\hat{R}$ and $\hat{R}$'.

The coordinate systems and wavefunctions representations for both atom-triatom and diatom-diatom scattering studies are outlined here. For details, please refer to [40, 41]. Figure 1 shows the reagent Jacobi coordinates (R, r1, r2, θa1, θa2, ψ) for the H+D2O atom-triatom arrangement, and product Jacobi coordinates (R', r'1, r'2, θ'1, θ'2, ψ') for the HD+OD diatom-diatom arrangement. As can be seen r2 and r'2 share the same vector.

The system Hamiltonian in the H+D2O arrangement for a given total angular momentum J can be written as

 (7)

where μ is the reduced mass between the center-of-mass of H and D2O, J is the total angular momentum operator of the system, j12 is the total angular momentum operator of D2O, j2 is the rotational angular momentum operator of OD, and j1=j12-j2 is the orbital angular momentum of D2O. The reference Hamiltonian hi(ri) (i=1, 2) is defined as

 (8)

The time-dependent wave function can be expanded in terms of the translational basis of R, the vibrational basis , and the body-fixed (BF) rovibrational eigenfunction as

 (9)

The BF total angular momentum eigenfunctions can be written as,

 (10)

where is the Wigner rotation matrix [42] with three Euler angles , ε is the parity of the system defined as ε=(-1)j1+j2+L with L=J-j12 being the orbital angular momentum of H and D2O, and Yj1j2j12 K is the angular momentum eigenfunction of j12 defined as,

 (11)

where yjm are spherical harmonics. Note that in Eq.(10), the restriction ε(-1)j1+j2+j12+J=1 for K=0 partitions the whole rotational basis set into even and odd parities. Thus a K=0 initial state can only appear in one of these two parity blocks. For K>0, however, there is no such restriction, the basis set is the same for even and odd parities. Hence a K>0 initial state can appear in both parities.

For the the HD+OD arrangement, similarly to the discussion about H+D2O arrangement, we have the same forms of equations, except the quantities with unprimed indices replaced by primed ones. The system Hamiltonian for the HD+OD arrangement, corresponding to Eq.(7),

 (12)

where μ' is the reduced mass between the center of mass of HD and OD, j'1 and j'2 are the rotational angular momentum operators of HD and OD, which are coupled to form j'12.

Corresponding to the Eq.(9), we have

 (13)

It should be pointed out here that the functional form for the product arrangement basis (primed) are different from the one for the reagent arrangement basis (unprimed).

BF total angular momentum basis in Eq.(13) are defined as,

 (14)

where ε is the parity of the system defined as ε=(-1)j'1+j'2+N with N=J-j'12 being the orbital angular momentum of HD and OD, and is the angular momentum eigenfunction of j'12,

 (15)

In our calculations we employ the CXZ PES [5] constructed by using neural network (NN) fitting [43-45] based on 17, 000 ab initio points calculated at UCCSD(T)-F12a/AVTZ level of theory. It is the best PES available for the system, substantially more accurate and smooth than YZCL2 and XXZ PES. We carry out state-to-state calculations for all the total angular momentum J up to 40 to converge DCSs for collision energies up to 1.5 eV for (000) initial state, and for collision energies up to 1.2 eV for (100) and (001) states. The numerical parameters used in reactant atom-triatom coordinates are as the following.

The interaction region was defined by a rectangular box of [1.0, 6.0]a0 in the r1 coordinate, and [1.0, 6.0]a0 for the R coordinate. The number of vibrational basis functions used was 65 for r1 coordinate. For the R coordinate, we used 40 sine discrete variable representation (DVR) points. The asymptotic region was defined from 6.0 a0 to 12.0 a0 with 42 sine DVR points for the R coordinate and 8 vibrational basis functions for the r1 coordinates. The number of vibrational basis functions for the nonreactive OD is 12 for interaction and 7 for asymptotic regions. For the rotational motion, we used j1max=28 and j2max=60, which results in 22475 rotational basis functions for K=0 and even parity. The number of K blocks used in the calculation increases with the total angular momentum, from 1 for J=0, up to 7 for J=40. The initial Gaussian wave packet was centered at R0=10 a.u. For lower total angular momentum, J, we propagated the wave packets for 6000 a.u. of time with a time step of 10 to converge the low energy reaction probability. The absorption potential on the r1 coordinate for the MRPD calculation takes the form:

 (16)

with r10=4.5 a.u., △r11=1.5 a.u., C=0.2, n=1.5. To minimize the computational cost, the coordinate transformation is carried out at every 8 propagation time steps.

Continuous propagation in the product diatom-diatom coordinates only involves a total number of 120 sine functions for the translational coordinate R' in a range of [3.0, 15.0]a0, 5 vibrational basis functions for HD bond, and 3 vibrational basis functions for OD bond. We used j'1max=30, j'2max=18, resulting in a total of 4750 coupled rotational basis functions for K=0 with even parity. The number of K blocks used in the calculation increases from 1 for J=0, up to 10 for J=40. A dividing surface is placed at R'=10.0 a.u. to extract S matrix elements. Since one OD bond is treated as a nonreactive bond in this work, the reaction probabilities and cross sections should be doubled.

III. Results

Figure 2 shows the total reaction probabilities for the (000), (100), and (001) three initial states as a function of collision energy on the CXZ PES, from the initial state selected wave packet (ISSWP) approach only using reactant Jacobi coordinates, in comparison with those from the RPD state-to-state calculations. As can be seen, for both the ground and excited initial states, the agreements between these two kinds of results are excellent in the entire energy region considered here, indicating the RPD partition of the wave function as well as the continuous propagation in diatom-diatom coordinates is of high accuracy. The reaction probability for the ground state increases steadily with increasing collision energy starting from about 0.9 eV, and is very small in the energy region considered here. The overall behavior of reaction probabilities for the first symmetric stretching excited (100) and asymmetric stretching excited (001) states is quite similar to that for the ground state, although the former is much larger and exhibits some oscillatory structures.

 FIG. 2 The total reaction probabilities for the three initial states (000), (100) and (001) of the title reaction for J=0 as a function of collision energy, obtained from the initial state selected wave packet approach only using reactant Jacobi coordinates (in solid lines), in comparison with those from the RPD state-to-state calculations (in circle symbols).

Figure 3(a) shows the ICS for the three initial states as a function of collision energy from the ISSWP approach only using reactant Jacobi coordinates, in comparison with those from the RPD state-to-state calculations. As can be seen, the RPD based state-to-state calculation is also of high accuracy for J>0. The overall behavior of the ICS is fairly similar to that of the total reaction probabilities shown in Fig. 2, except that the oscillatory structures on the (100) and (001) probability curves vanish on the cross section curves.

The ICS for (000) state is quite small with the value of roughly 0.023 a02 at Ec=1.5 eV, which is smaller than the experimental results by a factor of 1.3 to 4.8 (according to the data points in Ref.[46]). In addition, the reaction cross section for H+D2O is two to three times smaller than that for H+H2O, due to less mobility of D atom compared to H atom. We can also see from Fig. 3(a) that one-quantum excitations of stretching modes of D2O substantially enhance the reactivity. The enhancement factors are 37.7 at the collision energy of 1.2 eV. Zare and co-workers have previously measured that one-quantum excitation of the asymmetric stretch in D2O (or the combination of one quantum each of bend and asymmetric stretch) increases the cross section by at least a factor of 5 (and possibly a factor of 20) [8]. In the low collision energy region, the (001) state is slightly more reactive than the (100) state because of the higher excitation energy for the (001) state. Figure 3(b) shows the ICS for the three initial states as a function of total energy, which are measured with respect to the energy of the D2O ground rovibrational state. The ICS for (100) state is only slightly larger than that for (001) state in the high total energy region, indicating the first symmetric and asymmetric stretching excitation of D2O have nearly equal efficacy on promoting the reaction.

 FIG. 3 (a) ICS for the three initial states as a function of collision energy obtained from the initial state selected wave packet approach only using reactant Jacobi coordinates (in solid lines), in comparison with those from the RPD state-to-state calculations (in X symbols); (b) ICS for the three initial states as a function of total energy.

Figure 4 shows the DCS for the three initial states, each at four collision energies, in terms of surface plots for the product translational energy and angle distributions for the reaction. It can be seen at the first glance that the first symmetric (100) and asymmetric (001) stretching excited states have essentially identical DCS. Just like the ground initial state, the DCSs for these two states are strongly backward peaked at the threshold energy. As the collision energy increases, the angular distribution becomes increasingly broader, and the peak position starts to shift gradually to a smaller angle, consistent with the fact that the title reaction is a direct reaction via an abstraction mechanism. At the same collision energy, the DCS for the ground state is very different from those for the vibrationally excited state. However, one may note that the DCS at each column shown in Fig. 4 resemble each other somewhat in shape, indicating that the angular distributions for these three states are similar at the same total energy (because the excitation energy for (100) and (001) states are about 0.331 and 0.346 eV, close to the collision energy difference between the ground and vibrationally excited state shown in the figure).

 FIG. 4 Surface plots for the product translational energy and angle distributions for the three initial states (000), (100), and (001) at several collision energies.

In order to examine the effect of the stretching excitation on the angular distribution in a quantitative manner, we present the total DCS for the three states at collision energies of 0.9, 1.06, and 1.2 eV in Fig. 5(a). The behavior of angle distributions for the excited states is quite different from the ground state at the same collision energy, and the former are much broader and have obvious forward scattering components. The DCS for the (100) and (001) states are very close to each other, indicating that one quantum in the symmetric or antisymmetric stretching excitation of D2O has nearly identical effects to DCS.

 FIG. 5 (a) Differential cross sections for the three initial states at collision energies of 0.9, 1.06, and 1.2 eV. (b) Differential cross sections for the three initial states at the total energies of 0.9, 1.2, and 1.5 eV.

Figure 5 (b) shows the total DCS for the three states at total energies of 0.9, 1.2, and 1.5 eV. Overall, the angle distributions for (000), (100), and (001) vibrational states are quite similar at the same total energy of reaction, although these two stretch vibrational excitations enhance the cross section by at least 9 times relative to that for the ground state. Therefore, the energy initially deposited in stretching vibrations is much more efficient than the translational energy at promoting the reaction, but has rather similar effects on product angular distribution as the translational energy. In addition, the backward scattering components from the initial excited states decrease more slowly than that from the ground state as the total energy increases, and exhibit some oscillatory structures.

Figure 6 shows the product HD(v=1) and OD(v=1) vibrational state populations for the three initial states as a function of total energy. It should be noted that the population for (000) state in Fig. 6(b) was multiplied by a factor of 20 before plotting. When D2O is in the vibrational ground state, both HD and OD products are vibrationally very cold in the energy region considered here. There is a small fraction of HD(v=1) produced only for collision energy higher than 1.25 eV with a population of 4.3% at 1.5 eV, and almost no OD(v=1) produced. The initial vibrational excitations increase the branching ratios of HD(v=1) and OD(v=1) channels at the same total energy. The population of HD(v=1) reaches 10.6% at Etot=1.55 eV. However, the OD vibration excited populations remain small for the initial (100) and (001) states.

 FIG. 6 Populations of product HD(v=1) (a) and OD(v=1) (b) vibrational states for the three initial states as a function of total energy, together with the relative reactivity between the OD(v=0) and OD(v=1) state (violet dashed line in (b)). Note the population of OD(v=1) for (000) state shown in (b) was multiplied by a factor of 20 before plotting.

These small populations of OD(v=1) clearly reveal that D2O fundamental stretching excitations produce vibrationally "cold" OD as the ground initial state, in agreement with the experimental observation of Zare and coworkers [8]. In the normal mode picture, the (100) and (001) state correspond, respectively, to the symmetric and antisymmetric vibrational excitations of D2O with both OD bonds vibrationally excited. Breaking one OD bond by an incoming H atom will leave another OD bond vibrationally hot. However, in the local mode picture, the (100) and (001) states correspond to |10〉±=2-1/2(|10〉±|01〉), i.e. each OD bond is a superposition of the ground and first excited quantum states. When the incident H comes in, it will react with each quantum state independently without any effect from the non-reacting OD group if the non-reacting OD is a spectator, producing OD in the corresponding quantum state. Because the reactivity of the OD excited state is much higher than the ground state, the incident H atom tends to react with the excited state, producing OD mainly in the corresponding ground state. On the other hand, the H atom also reacts with the ground OD state, producing a small fraction of the OD(v=1) state. The fraction of OD(v=1) should be approximately equal to the relative reactivity between the OD(v=0) and OD(v=1) states. As can be seen from Fig. 6(b), the fraction of the OD(v=1) population is very close to the relative reactivity. Therefore, the local mode picture works very well here, and the non-reacting OD bond does act as a spectator.

Figure 7 shows the normalized rotational distributions for HD (a) and OD (b) product for the three initial states at some total energy. For the ground state reaction, the general trend in the change of HD or OD rotational distributions with the increase of collision energy is quite similar: the distribution becomes broader, and the jHD=0 or jOD=0 population declines. The jHD=0 population decreases from 29% at Ec=0.9 eV to 2% at Ec=1.5 eV, and the jOD=0 population decreases from 26% to 7%. The OD rotation distribution has a maximum at jOD=2 for all the energies shown, except for Ec=0.9 eV, for which jOD=0 is marginally larger than jOD=2 because of the oscillatory behavior. In contrast, the most probable rotational state for HD becomes higher with the increase of collision energy, changing from jHD=1 at Ec=0.9 eV to jHD=7 at Ec=1.5 eV. We can also see from Fig. 7(b) that the OD products from the (100) and (001) reactions have essentially the same rotational distribution as the ground state reaction, indicating one-quantum excitations of stretching modes of D2O has little effect on the OD rotational distribution.

 FIG. 7 Rotational state distributions for (a) HD and (b) OD product for the three initial states at the total energies of 0.9, 1.2, and 1.5 eV.

This is consistent with the experimental observation of Zare and coworkers [8], and further supports the picture of the non-reacting OD bond as a spectator. Unexpectedly, the HD rotation distribution from vibrationally excited reaction becomes less excited at the same total energy, which is different from the H+H2O reaction [22].

The differences in the vibrational and rotational excitation of HD and OD can also be found from the energy partition information. Figure 8 shows the fraction of total available energy in the product channel going into the rotation and vibration excitation of HD and OD as a function of total energy. For the ground state reaction, less than 32% of the available energy goes into the internal motion of the products in the energy region considered here, and only a very small fraction appears as vibrational energy. The HD rotation excitation carries away an increasingly large portion of the internal energy with the increase of collision energy, while the fraction of energy going into the OD rotation excitation essentially does not change with the collision energy. For the excited state reactions, the OD rotational distributions are close to that for (000) at the same total energy, leading to a small difference in the fraction of energy. The initial vibrational excitations reduce the fraction of HD rotation energy, but increase the fraction of HD vibration energy by a factor of 2.1 than that from the ground state D2O at Etot=1.5 eV.

 FIG. 8 The fraction of the total available energy in the product channel going into the rotation (a) and vibration (b) excitation of HD and OD as a function of the total energy. The corresponding total available energy in the product channel is also shown on the top of each panel. Note the fraction for OD vibration for the ground state was multiplied by a factor of 10 before plotting.
IV. Conclusions

We employ the TDWP method to calculate the first full-dimensional state-to-state DCS for the H+D2O→HD+OD reaction with D2O in the ground and the first symmetric and asymmetric stretching excited states on the CXZ PES. The calculated DCS for these three initial state are strongly backward peaked at low collision energies, become increasingly broader with the increase of collision energy, with the peak position shifting gradually to a smaller angle, consistent with the fact that the title reaction is a direct reaction via an abstraction mechanism. Our calculations revealed that the stretching excitation is much more efficient at promoting the reactivity than the translational energy, but has a rather similar effect on the angular distribution as the translational energy. Moreover, the DCS for the (100) and (001) states are very close to each other, indicating that one quantum in the symmetric or antisymmetric stretching excitation of D2O has nearly identical effects to DCS as found in the H+H2O reaction \cite{liushu2016HH2O}. Product vibrational and rotational state distributions and energy partitioning information are also presented. The small populations of OD(v=1) from (100) and (001) initial states can be explained clearly by the local mode picture of D2O vibration, and reveals that the non-reacting OD does act as a spectator in the reaction. Unexpectedly, the initial vibrational excitation reduces the rotation excitation of HD at the same total energy. The present full-dimensional state-to-state quantum dynamics study for the title reaction is of considerable importance to our understanding of mode-specific chemistry.

V. Acknowledgements

This work was supported by the National Natural Science Foundation of China (No.21403223, No.21433009, and No.91221301), the Ministry of Science and Technology of China (No.2013CB834601), and the Chinese Academy of Sciences.

Reference
 [1] C. Schatz and H. Elgersman G., Chem. Phys. Lett. 73 , 21 (1980). DOI:10.1016/0009-2614(80)85193-1 [2] O. de Aspuru and D. C. Clary G., J. Phys. Chem A102 , 9631 (1998). [3] G. Wu, G. C. Schatz, G. Lendvay, D. C. Fang, and L. B. Harding, J. Chem. Phys. 113 , 3150 (2000). DOI:10.1063/1.1287329 [4] M. Yang, D. H. Zhang, M. A. Collins, and S. Y. Lee, J. Chem. Phys. 115 , 174 (2001). DOI:10.1063/1.1372335 [5] J. Chen, X. Xu, and D. H. Zhang, J. Chem. Phys. 138 , 154301 (2013). DOI:10.1063/1.4801658 [6] A. Sinha, J. Phys. Chem. 94 , 4391 (1990). DOI:10.1021/j100374a005 [7] A. Sinha, M. C. Hsiao, and F. F. Crim, J. Chem. Phys. 94 , 4928 (1991). DOI:10.1063/1.460578 [8] M. J. Bronikowski, W. R. Simpson, and R. N. Zare, J. Phys. Chem. 97 , 2204 (1993). DOI:10.1021/j100112a022 [9] A. Sinha, M. C. Hsiao, and F. F. Crim, J. Chem. Phys. 92 , 6333 (1990). DOI:10.1063/1.458312 [10] M. J. Bronikowski, W. R. Simpson, B. Girard, and R. N. Zare, J. Chem. Phys. 95 , 8647 (1991). DOI:10.1063/1.461243 [11] M. J. Bronikowski, W. R. Simpson, and R. N. Zare, J. Phys. Chem. 97 , 2194 (1993). DOI:10.1021/j100112a021 [12] Kudla and G. C. Schatz K., Chem. Phys. Lett. 193 , 507 (1992). DOI:10.1016/0009-2614(92)85840-7 [13] Kudla and G. C. Schatz K., J. Chem. Phys. 98 , 4644 (1993). DOI:10.1063/1.464992 [14] D. Troya, González M., and G. C. Schatz, J. Chem. Phys. 114 , 8397 (2001). DOI:10.1063/1.1366334 [15] D. C. Clary, Chem. Phys. Lett. 192 , 34 (1992). DOI:10.1016/0009-2614(92)85423-8 [16] M. Bowman and D. Wang J., J. Chem. Phys. 96 , 7852 (1992). DOI:10.1063/1.462382 [17] Wang and J. M. Bowman D., J. Chem. Phys. 98 , 6235 (1993). DOI:10.1063/1.464817 [18] H. Zhang and J. C. Light D., and Chem. Soc. J., Faraday Trans. 93 , 691 (1997). DOI:10.1039/a605888d [19] B. Jiang, D. Xie, and H. Guo, J. Chem. Phys. 135 , 084112 (2011). DOI:10.1063/1.3626525 [20] Fu and D. H. Zhang B., J. Chem. Phys. 138 , 184308 (2013). DOI:10.1063/1.4803695 [21] Fu and D. H. Zhang B., J. Chem. Phys. 142 , 064314 (2015). DOI:10.1063/1.4907918 [22] Liu and D. H. Zhang S., Chem. Sci. 7 , 261 (2016). DOI:10.1039/C5SC03472H [23] H. Zhang and S. Y. Lee D., J. Chem. Phys. 110 , 4435 (1999). DOI:10.1063/1.478327 [24] M. Goldfield and S. K. Gray E., J. Chem. Phys. 117 , 1604 (2002). DOI:10.1063/1.1487824 [25] D. H. Zhang, M. Yang, and S. Y. Lee, J. Chem. Phys. 117 , 10067 (2002). DOI:10.1063/1.1519009 [26] B. Fu, Y. Zhou, and D. H. Zhang, Chem. Sci. 3 , 270 (2011). [27] Jiang and H. Guo B., J. Am. Chem. Soc. 135 , 15251 (2013). DOI:10.1021/ja408422y [28] D. H. Zhang, D. Xie, M. Yang, and S. Y. Lee, Phys. Rev. Lett. 89 , 283203 (2002). DOI:10.1103/PhysRevLett.89.283203 [29] D. H. Zhang, J. Chem. Phys. 125 , 133102 (2006). DOI:10.1063/1.2217439 [30] T. CvitaŠ and S. C. Althorpe M., J. Phys. Chem A113 , 4557 (2009). [31] T. CvitaŠ and S. C. Althorpe M., Phys. Scr. 80 , 048115 (2009). DOI:10.1088/0031-8949/80/04/048115 [32] T. CvitaŠ and S. C. Althorpe M., J. Chem. Phys. 134 , 024309 (2011). DOI:10.1063/1.3525541 [33] C. Xiao, X. Xu, S. Liu, T. Wang, W. Dong, T. Yang, Z. Sun, D. Dai, D. H. Zhang, and X. Yang, Science 333 , 440 (2011). DOI:10.1126/science.1205770 [34] S. Liu, X. Xu, and D. H. Zhang, J. Chem. Phys. 136 , 144302 (2012). DOI:10.1063/1.3701266 [35] S. Liu, C. Xiao, T. Wang, J. Chen, X. Xu, D. H. Zhang, and X. Yang, Faraday Discuss. 157 , 101 (2012). DOI:10.1039/c2fd20018j [36] Peng and J. Z. H. Zhang T., J. Chem. Phys. 105 , 6072 (1996). DOI:10.1063/1.472444 [37] Y. Huang, W. Zhu, D. J. Kouri, and D. K. Hoffman, J. Phys. Chem. 98 , 1868 (1994). DOI:10.1021/j100058a025 [38] S. C. Althorpe, D. J. Kouri, and D. K. Hoffman, J. Chem. Phys. 106 , 7629 (1997). DOI:10.1063/1.473766 [39] D. Feit and J. A. Fleck M., J. Chem. Phys. 78 , 301 (1983). DOI:10.1063/1.444501 [40] H. Zhang and J. Z. H. Zhang D., J. Chem. Phys. 101 , 1146 (1994). DOI:10.1063/1.467808 [41] H. Zhang and J. C. Light D., J. Chem. Phys. 104 , 4544 (1996). DOI:10.1063/1.471203 [42] M. E. Rose, Elementary Theory of Angular Momentum. New York: John Wiley (1995). [43] T. Hagan and M. B. Menhaj M., and Networks Neural, IEEE Trans. 5 , 989 (1994). [44] W. S. Sarle, in Proceedings of the 27th Symposium on the Interface, (1995). [45] L. M. Raff, M. Malshe, M. Hagan, D. I. Doughan, M. G. Rockley, and R. Komanduri, J. Chem. Phys. 122 , 084104 (2005). DOI:10.1063/1.1850458 [46] S. Koppe, T. Laurent, P. D. Naik, H. Volpp, and J. Wolfrum, Can. J. Chem. 72 , 615 (1994). DOI:10.1139/v94-086 [47] D. H. Zhang, M. A. Collins, and S. Y. Lee, Science 290 , 961 (2000). DOI:10.1126/science.290.5493.961
H+D2O→HD+OD反应的微分截面:一种全维态-态量子动力学计算

a. 中国科学院大连化学物理研究所, 分子反应动力学国家重点实验室和理论与计算化学中心, 大连 116023;
b. 中国科学院大学, 北京 100049