Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D2O with Normal-H2

We present a combined experimental and theoretical study on the rotationally inelastic scattering of heavy water, D2O, with normal-H2. Crossed-molecular beam measurements are performed in the collision energy range between 10 and 100 cm−1, corresponding to the near-threshold regime in which scattering resonances are most pronounced. State-to-state excitation cross-sections are obtained by probing three low-lying rotational levels of D2O using the REMPI technique. These measurements are complemented by quantum close-coupling scattering calculations based on a high-accuracy D2O–H2 interaction potential. The agreement between experiment and theory is within the experimental error bars at 95% confidence intervals, leading to a relative difference of less than 7%: the near-threshold rise and the overall shape of the cross-sections, including small undulations due to resonances, are nicely reproduced by the calculations. Isotopic effects (D2O versus H2O) are also discussed by comparing the shape and magnitude of the respective cross-sections.


Introduction
Water is the third most abundant molecule in the interstellar medium (ISM) and has ubiquitously been observed by ground-and space-based telescopes since its first detection in 1969 in the Orion nebula [1,2]. Water thus is a key molecule for understanding the energy balance and the physical-chemical processes that occur in these environments. Despite the low elemental deuterium abundance in the Galaxy, D/H∼10 −5 , a spectacular deuterium enrichment of various interstellar molecules, has been observed in star forming regions. Among these molecules, the deuterated water isotopologues HDO and D 2 O are of special importance because they can help to understand the origin of water in protostars [3]. D 2 O was identified for the first time in 2007 toward IRAS 16293-2422 [4] in its para form and detected in its ortho form in 2010 in the same cold and diluted protostar envelope [5], which has then allowed to derive the ortho-to-para ratio of D 2 O (< 2.6. For comparison, the thermal equilibrium value is 2).
However, in order to interpret the astrophysical observations [2] in terms of relative abundances or even local physical conditions, radiative transfer modeling is necessary, as the molecular populations are generally far from local thermodynamical equilibrium. At low density, the heavy water population is indeed determined by a competition between the radiative and collisional processes, thus requiring the knowledge of rates for collisional (de)excitation. Its principal collision partner obviously is H 2 because of its high abundance in the cold ISM.
Inelastic collisions between molecules are fundamental processes in which energy is transferred between their relative translational motion and their internal degrees of freedom. Water-hydrogen collisions play an essential role in a large variety of properties Previous theoretical studies of the H2O-H2 interactions are too extensive to be enumerated in detail herein. The reader may refer to the last articles on rotational inelastic collisions and citations herein [8,9] and, therefore, we will limit our review to only a portion of the literature pertinent to our experiments, mainly the various experimental results Previous theoretical studies of the H 2 O-H 2 interactions are too extensive to be enumerated in detail herein. The reader may refer to the last articles on rotational inelastic collisions and citations herein [8,9] and, therefore, we will limit our review to only a portion of the literature pertinent to our experiments, mainly the various experimental results on H 2 O-H 2 and the theoretical calculations and experiments on D 2 O-H 2 .
In the first molecular beam scattering experiment on H 2 O + H 2 elastic differential crosssections [10], well resolved diffractive oscillations were observed, as well as the rainbow maximum. The spherically symmetric model potential parameters (e.g., the well depth and location of zero potential) were deduced assuming a Lennard-Jones (12,6) potential. The H 2 -H 2 O van der Waals complex was then studied by infrared spectroscopy in the gas phase by the group of D. Nesbitt [11,12] and also in matrix [13]. There are also studies on H 2 O pressure broadening by H 2 [14,15] or inelastic collisions [16,17]. The Perugia group has reported determination of the potential parameters for the isotropic component of the D 2 O-D 2 interaction by elastic scattering studies [18]. The group of Nesbitt and coworkers has also studied the infrared spectrum of the H 2 O-D 2 complex in the bend mode frequency region [19].
A large number of studies can thus be found in the literature providing information on H 2 O-H 2 , whereas calculations or experiments on D 2 O-H 2 are scarce. The complex has been studied in solid neon matrix [13] where the infrared absorptions principally arose from complexes involving ortho-H 2 , for which j H2 = 1. Excitation rate constants of the D 2 O low lying rotational transitions by para-H 2 were calculated for temperatures between 1 and 30 K [20] and later completed for temperatures below 100 K [21]. The calculations at the close-coupling level were performed using the highly accurate full 9D (but vibrationally averaged to 5D) PES of Valiron et al. [22]. In the cold or near cold regime (5-30 K, relevant for astrophysical applications), prominent open channel (orbiting) and closed channel (Feshbach) resonances were observed. The first state-to-state differential cross-sections (DCSs) for rotationally inelastic scattering of D 2 O(0 00 or 1 01 ) + H 2 (j H2 = 0, 2 or 1) were measured by Sarma et al. [23] at the collision energy of 584 cm −1 . At collision energies above 500 cm −1 , this crossed-beam study could only probe the PES in the short-range and well regions. The first experimental check at long range was achieved by measuring the state-to-state scattering ICSs of rotationally inelastic excitation of D 2 O in the 2 02 level by collisions with para-H 2 [8], in the very low collision energy range 10-100 cm 1 (near rotational thresholds). The observed scattering resonances were unambiguously identified by quantum mechanical calculations performed on the PES of Valiron et al. [22]. The aim of this study is (i) to present the full set of experimental ICSs for the D 2 O low lying rotational transitions colliding with normal-H 2 in the vicinity of the thresholds and (ii) to provide new stringent tests of theory as the experimental ICS resonance structures (amplitude and position) are very sensitive to the shape and even to tiny details of the PES.
In the results section, our beam characteristics and the experimentally measured crossed beam ICSs will be presented followed by the main theoretical results in relation with our experimental studies. The experimental observations will be compared with the inelastic scattering calculations in the discussion section.

Experimental Results
The experimental results were obtained in a crossed-molecular beam (CMB) experiment under single collision conditions (see Figure 2). Our machine allows inelastic processes [24] to be studied at low collision energies. For this purpose, two supersonic molecular beams, one composed of deuterated water seeded in a carrier gas and the other of pure hydrogen, are generated with a high velocity resolution, while most of the initial population of molecular species being in their ground rotational state. The collision energy of the colliding partners is tuned by varying the beam intersection angle, γ: where v D2O and v H2 are the velocities of the D 2 O and H 2 beams, µ is the reduced mass of the system D 2 O-H 2 , and v r is their relative velocity. The state-to-state ICS is measured by resonance-enhanced multiphoton ionization (REMPI) coupled with time-of-flight mass spectrometry (TOF-MS), using tunable UV/VUV lasers (see the Methods section for more details). We describe the characteristics of each beam, which has an impact on the experimental results before presenting the experimental ICSs. lecular beams, one composed of deuterated water seeded in a carrier gas and the other of pure hydrogen, are generated with a high velocity resolution, while most of the initial population of molecular species being in their ground rotational state. The collision energy of the colliding partners is tuned by varying the beam intersection angle, γ: where vD2O and vH2 are the velocities of the D2O and H2 beams, µ is the reduced mass of the system D2O-H2, and vr is their relative velocity. The state-to-state ICS is measured by resonance-enhanced multiphoton ionization (REMPI) coupled with time-of-flight mass spectrometry (TOF-MS), using tunable UV/VUV lasers (see the Methods section for more details). We describe the characteristics of each beam, which has an impact on the experimental results before presenting the experimental ICSs.
(a) (b) The D 2 O or normal-H 2 beams are generated using Even-Lavie valves. The measured characteristics of the gas pulses, velocity, pulse duration, velocity spread, and angle spread (see Section 4) are reported in Table 1.  Table 1). From the REMPI spectra reported in Figure 3, the D 2 O rotational distribution was determined using the spectroscopic data of the C 1 B 1 , v = 0 ← X 1 A 1 , v = 0 transitions, from Yang et al. [25] provided in PGopher software [26]. In our previous study on D 2 O + para-H 2 [8], a rotational temperature of 12 K was determined for the D 2 O beam seeded in He, with a spin temperature of 320 K, corresponding to the reservoir temperature. The rotational populations were thus 0.55, 0.19, 0.15, 0.05, and 0.04, for the 0 00 , 1 01 , 1 11 , 1 10 , and 2 02 levels, respectively. Further investigation of the contributions of the different rotational levels on the ICSs for the transitions to the D 2 O 2 02 rotational level allowed us to refine the rotational population ratio of D 2 O (0 00 :1 11 ) and (1 01 :1 10 ) to 0.81:0.19 and 0.73:0.27, respectively. New experimental inelastic collision studies were performed with this beam of D 2 O diluted in He and supplemented with D 2 O seeded in Ne. For the latter beam, the REMPI spectrum is given in Figure 3. Obviously, only the lowest rotational levels of ortho-D 2 O and para-D 2 O are populated: the inset is the enlargement of the spectrum range where the first excited rotational levels may be seen and were estimated to have a population of less than 2%. As these signals are almost within the noise (see the inset of the D 2 O spectrum) even if the transitions are favored for the excited states, only the ground rotational levels will significantly contribute to the inelastic collisions.  Indeed, rotational cooling is more efficient when water molecules are seeded in Ne rather than in He [25]. However, even if the cooling is rapid, hence avoiding condensation, dimers may be formed, i.e., a weekly bound van der Waals structure, which presents shallow minima of ~34 and 7 cm −1 for D2O-Ne [27] and D2O-He clusters [28], respectively. Usually, a large cluster binding energy enhances the formation of clusters, which release their condensation energy into the beam, thence limiting translational cooling. However, the characteristics of the beams are equivalent with or without D2O. Moreover, no signal at higher m/z than those of D2O + was detected. In conclusion, there is no evidence of significant amounts of water clusters in the beam.
As carrier gas mass changes, the beam velocity changes, too (see Table 1). The terminal velocity v of a supersonic beam generated by an adiabatic and isentropic expansion of an ideal gas through a nozzle, from a reservoir at a temperature T0, is defined by:  Indeed, rotational cooling is more efficient when water molecules are seeded in Ne rather than in He [25]. However, even if the cooling is rapid, hence avoiding condensation, dimers may be formed, i.e., a weekly bound van der Waals structure, which presents shallow minima of~34 and 7 cm −1 for D 2 O-Ne [27] and D 2 O-He clusters [28], respectively. Usually, a large cluster binding energy enhances the formation of clusters, which release their condensation energy into the beam, thence limiting translational cooling. However, the characteristics of the beams are equivalent with or without D 2 O. Moreover, no signal at higher m/z than those of D 2 O + was detected. In conclusion, there is no evidence of significant amounts of water clusters in the beam. As carrier gas mass changes, the beam velocity changes, too (see Table 1). The terminal velocity v of a supersonic beam generated by an adiabatic and isentropic expansion of an ideal gas through a nozzle, from a reservoir at a temperature T 0 , is defined by: where the temperature of the beam is T << T 0 (substantial cooling in the expansion), and C p is the average molar heat capacity (for an atom, C p = 5/2 R, where R stands for the perfect gas constant). This equation implies a velocity ratio of √ (20/4) = 2.2 for the two carrier gas, which is experimentally observed (see Table 1, where 1858/853 = 2.2).

Normal-H 2 Beam
As the best collisional energy resolution is obtained when the two beams have almost the same velocity (Figure 2), for the H 2 beam, a cryogenically cooled valve was used with He closed-cycle cryostats (Equation (2)). The temperature of the cold head was reduced at its minimum, which corresponds to a setpoint of 10 K to reach the lowest hydrogen beam velocity, which is almost equal to the neon beam velocity (Table 1). In the case of water diluted in He, a setpoint at 145 K was used. However, even if the temperature in the reservoir is decreased to a setpoint of 145 K or 10 K, there is no ortho-to-para conversion of the hydrogen. Normal-H 2 beam remains a mixture of paraand ortho-H 2 with the room temperature ratio, i.e., 3:1. The H 2 beam was probed ( Figure 4) when recording (C 1 Π u ; v = 2 ← X 1 Σ + g ; v = 0) R(0) and R(1) transitions near 96.47 nm, with a (3 + 1) REMPI scheme at 289.4 nm [29]. No evidence of R(2) signal was found for the beam with the highest reservoir temperature: a normal-H 2 beam is composed of 25% of H 2 (j H2 = 0) and 75% of H 2 (j H2 = 1). He closed-cycle cryostats (Equation (2)). The temperature of the cold head was reduced at its minimum, which corresponds to a setpoint of 10 K to reach the lowest hydrogen beam velocity, which is almost equal to the neon beam velocity (Table 1). In the case of water diluted in He, a setpoint at 145 K was used. However, even if the temperature in the reservoir is decreased to a setpoint of 145 K or 10 K, there is no ortho-to-para conversion of the hydrogen. Normal-H2 beam remains a mixture of para-and ortho-H2 with the room temperature ratio, i.e., 3:1. The H2 beam was probed ( Figure 4) when recording (C 1 Πu; v = 2 ← X 1 Σ + g; v = 0) R(0) and R(1) transitions near 96.47 nm, with a (3 + 1) REMPI scheme at 289.4 nm [29]. No evidence of R(2) signal was found for the beam with the highest reservoir temperature: a normal-H2 beam is composed of 25% of H2(jH2 = 0) and 75% of H2(jH2 = 1).

Integral Cross-Sections
Three rotational levels of D 2 O were probed: 1 11 , 2 02 for ortho-D 2 O, and 1 10 for para-D 2 O. The REMPI intensities were recorded versus the angle relative to the beam axes. Collision energies were derived from Equation (1). The cross-section was obtained from the REMPI signal intensity, divided by the relative velocity of D 2 O and H 2 and the interaction time of the two beams (see the Methods section for more details).
The beam characteristics are given in Table 1. In Figure 5, whatever the seeded gas was for the water beam and thus the temperature of the hydrogen beam, the cross-sections are well reproduced: the thresholds are at the correct energies (20.3 and 10.6 cm −1 for the 0 00 -1 11 and 1 01 -1 10 excitation, respectively), and we can observe some resonances at the same energies. There is thus no evidence of any clusters in the beams, as the binding energies of He-D 2 O and Ne-D 2 O differ by 30 cm −1 , and the D 2 O + H 2 resonances are observed at the same collision energies for the two water carrier gases. However, compared to the ICSs measured when D 2 O is seeded in Ne, they are broadened and slightly less marked in the case of D 2 O seeded in He. In fact, the collision energy spread is larger ( Figure 6), resulting in a spread-out of the resonance features over a more extended collision energy range. Moreover, the signal starts from a negative value in Figure 5a: the 1 11 rotational population starts decreasing when H 2 collides with D 2 O. Indeed, in the D 2 O:He beam, the first two rotational levels of ortho-D 2 O and para-D 2 O are populated, leading to the de-excitation of the 1 11 excited state until the collision energy is large enough to promote the excitation from the ground state (0 00 → 1 11 ). When the 1 10 excited rotational level is probed (Figure 5b), the lowest collision energy, which corresponds to our lowest angle between the two beams, is above the threshold (1 01 → 1 10 ), not allowing the observation of the de-excitation. In the case of D 2 O:Ne, only the ground rotational levels are populated (see Figure 3). Thus, in Figure 5, there is only one contribution from the ground rotational level to the probed excited one: the 0 00 → 1 11 and 1 01 → 1 10 transitions in Figure 5a,b, respectively. ICSs measured when D2O is seeded in Ne, they are broadened and slightly less marked in the case of D2O seeded in He. In fact, the collision energy spread is larger ( Figure 6), resulting in a spread-out of the resonance features over a more extended collision energy range. Moreover, the signal starts from a negative value in Figure 5a: the 111 rotational population starts decreasing when H2 collides with D2O. Indeed, in the D2O:He beam, the first two rotational levels of ortho-D2O and para-D2O are populated, leading to the de-excitation of the 111 excited state until the collision energy is large enough to promote the excitation from the ground state (000 → 111). When the 110 excited rotational level is probed (Figure 5b), the lowest collision energy, which corresponds to our lowest angle between the two beams, is above the threshold (101 → 110), not allowing the observation of the deexcitation. In the case of D2O:Ne, only the ground rotational levels are populated (see Figure 3). Thus, in Figure 5, there is only one contribution from the ground rotational level to the probed excited one: the 000 → 111 and 101 → 110 transitions in Figure 5a,b, respectively. Each experimental ICS datum corresponds to an average over a Gaussian distribution of collision energies: the half width at 1/e height, called collision energy spread, is given in Figure 6b. More details are given in the Methods section. This collision energy spread originates from the velocities and angular spread (Table 1). It will also affect the average relative velocity and the interaction time, which is given in Figure 6a. As the deconvolution of the experimental points is not possible, the theoretical cross-sections are convoluted with the collision energy spread to allow for comparison with the experimental data. Moreover, as previously explained, the rotational population distribution is

Theoretical Results
Inelastic scattering calculations were performed in a similar way to Bergeat et al. [8], and a summary can be found in the Methods section. The theoretical cross-sections for three particular rotational excitations from the ortho-D2O ground state by collisions with H2 are plotted in Figures 2 and 3.
Calculations were performed for H2 in its ground para (jH2 = 0) and ortho (jH2 = 1) states, and cross-sections for normal-H2 were obtained for an ortho-to-para ratio fixed at 3. The rotational basis set for D2O included angular momenta up to jD2O = 6, corresponding to the lowest 25 para-and lowest 25 ortho-D2O levels. The lowest two levels were included in the basis set for para-H2 (jH2 = 0, 2) and ortho-H2 (jH2 = 1, 3). Enlarging the D2O or H2 basis sets was found to be negligible at the investigated energies, with cross-sections changing by less than one percent. We note, however, that because of the potential well (De = 238 cm −1 ) and anisotropy of the D2O-H2 PES, the inclusion of closed-channels was found to be crucial, even for ortho-H2.
As previously noted for the H2O-H2 system [17] and also for other molecule-H2 systems [30], resonances are more conspicuous for para-H2 (jH2 = 0) than for ortho-H2 (jH2 = 1). This is particularly striking here for the 000 → 202 transition (Figure 7) where resonances have magnitudes up to twice that of the background cross-section for para-H2 (jH2 = 0) while for ortho-H2 (jH2 = 1) are much less intense. This can be explained by the non-vanishing quadrupole of ortho-H2 (jH2 = 1) and the corresponding stronger binding with ortho-H2 (jH2 = 1) than with para-H2 (jH2 = 0). This was observed in the H2O-H2 complex where the dissociation energy and number of bound-states are indeed significantly larger for ortho-H2 (jH2 = 1) monomers than para-H2 (jH2 = 0) monomers [31]. Overlapping and broadened resonances are thus expected to be more numerous in D2O-ortho-H2 (jH2 = 1) with respect to D2O-para-H2 (jH2 = 0). In addition, it is observed that the background cross-sections for para-H2 (jH2 = 0) and ortho-H2 (jH2 = 1) have similar magnitudes for the transition 000 → 202 while ortho-H2 (jH2 = 1) cross-sections are much larger than para-H2 (jH2 = 0) cross-sections for other transitions, e.g., 000 → 111 and 111 → 202 ( Figure 8). This is clearly due to the dipolequadrupole interaction, which dominates in D2O dipolar transitions (i.e., those with ΔjD2O = 0, ±1, ∆ka = ±1, ∆kc = ±1) while D2O quadrupolar transitions (∆jD2O = ±2) are much less sensitive to the quadrupole of H2. Each experimental ICS datum corresponds to an average over a Gaussian distribution of collision energies: the half width at 1/e height, called collision energy spread, is given in Figure 6b. More details are given in the Methods section. This collision energy spread originates from the velocities and angular spread (Table 1). It will also affect the average relative velocity and the interaction time, which is given in Figure 6a. As the deconvolution of the experimental points is not possible, the theoretical cross-sections are convoluted with the collision energy spread to allow for comparison with the experimental data. Moreover, as previously explained, the rotational population distribution is also taken into account.

Theoretical Results
Inelastic scattering calculations were performed in a similar way to Bergeat et al. [8], and a summary can be found in the Methods section. The theoretical cross-sections for three particular rotational excitations from the ortho-D 2 O ground state by collisions with H 2 are plotted in Figures 2 and 3.
Calculations were performed for H 2 in its ground para (j H2 = 0) and ortho (j H2 = 1) states, and cross-sections for normal-H 2 were obtained for an ortho-to-para ratio fixed at 3. The rotational basis set for D 2 O included angular momenta up to j D2O = 6, corresponding to the lowest 25 paraand lowest 25 ortho-D 2 O levels. The lowest two levels were included in the basis set for para-H 2 (j H2 = 0, 2) and ortho-H 2 (j H2 = 1, 3). Enlarging the D 2 O or H 2 basis sets was found to be negligible at the investigated energies, with cross-sections changing by less than one percent. We note, however, that because of the potential well (D e = 238 cm −1 ) and anisotropy of the D 2 O-H 2 PES, the inclusion of closed-channels was found to be crucial, even for ortho-H 2 .

Discussion
To compare the experimental results with the theoretical cross-sections, the

Discussion
To compare the experimental results with the theoretical cross-sections, the calculated ICSs were convoluted with the experimental collision energy spread, and the different contributions of the rotational populations were added or subtracted with the

Discussion
To compare the experimental results with the theoretical cross-sections, the calculated ICSs were convoluted with the experimental collision energy spread, and the different contributions of the rotational populations were added or subtracted with the appropriated weights (Table 2), as explained in the Methods section. To begin, the para to ortho ratios were kept constant at their values at 320 K as the ortho-para conversion rates are extremely slow without any magnetic catalyst: 2:1 for D 2 O and 3:1 for H 2 (see Table 2). In Figure 9a, only one state-to-state contribution is observed, the D 2 O(1 01 ) + normal-H 2 → D 2 O(1 10 ) + normal-H 2 . This is because only the lowest paraand ortho-D 2 O rotational levels are significantly populated when seeded in Ne. The overall excitation function, i.e., the evolution of the ICS versus the collision energy, is in good agreement with the convoluted theoretical calculations. In the case of the excitation transition from D 2 O(0 00 ) to D 2 O(1 11 ) shown in Figure 10a, small undulations due to resonances are observed in perfect agreement with the quantum mechanical calculations. However, we should underline that a small contribution of 0.12 on the 1 01 to 2 12 transition was added to take into account the overlap of the REMPI peak from the probed 1 11 rotational level of D 2 O with those from the 2 12 level. Then, the REMPI intensities for the D 2 O(1 10 ) were recorded when heavy water was seeded in He ( Figure 9b). Due to higher beam velocities, the minimum collision energy probed is above threshold. Furthermore, the signal is small as some population in the excited state is lost by de-excitation to the ground state and excitation to the next upper state. No real structure was observed. The 1 10 to 1 01 population ratio of the water beam in He was deduced from this excitation function and used for the other fits, despite the uncertainty. The 1 11 to 0 00 population ratio was then fitted from the 2 02 experimental ICSs (Figure 11), using the overlap value of the REMPI peaks from 2 12 and 2 02 determined previously [8]. All these population ratios and REMPI peak overlaps summarized in Table 2 were also used to analyze the ICSs of the D 2 O(1 11 ) transitions when water is seeded in He. The very good agreement between theory and experiment in Figure 10b, confirmed the accuracy of our experimental treatment and the theoretical calculations. The resonances observed for the 0 00-2 02 rotational inelastic cross-sections occur because of quasi-bound states become transitorily accessible during the collision. They were characterized by analyzing partial cross-sections differing by their total angular momentum J [8]. The total cross-section (summed over all partial waves J) was found to show, on resonance, large contributions from specific orbital angular momentum l, the rotational states of H 2 , and D 2 O being j D2O = 0 and j H2 = 0. An analysis of the resonances was put forward in the preceding paper [8].
certainty. The 111 to 000 population ratio was then fitted from the 202 experimental ICSs (Figure 11), using the overlap value of the REMPI peaks from 212 and 202 determined previously [8]. All these population ratios and REMPI peak overlaps summarized in Table 2 were also used to analyze the ICSs of the D2O(111) transitions when water is seeded in He. The very good agreement between theory and experiment in Figure 10b, confirmed the accuracy of our experimental treatment and the theoretical calculations. Negative values for the cross-sections mean that the 110 rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions in blue, the 101 → 110 transition, in blue dotted line, the loss is due to the 110 → 101 transition, and in gray, the loss is due to the 110 → 212 transition, assuming the initial rotational population given in Table 2. Negative values for the cross-sections mean that the 1 10 rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions in blue, the 1 01 → 1 10 transition, in blue dotted line, the loss is due to the 1 10 → 1 01 transition, and in gray, the loss is due to the 1 10 → 2 12 transition, assuming the initial rotational population given in Table 2.  Table 2.  Negative values for the cross-sections mean that the 1 11 rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions of the 0 00 → 1 11 transition (in orange solid line), the 1 01 → 2 12 transition (in magenta dotted line), 1 10 → 2 12 transition (in gray), the loss due to the 1 11 → 2 02 transition (in olive), and the loss due to the 1 11 → 2 12 transition (in orange dotted line), assuming the initial rotational population and a partial contribution from the 2 12 state, when probing the 1 11 state, given in Table 2.
Molecules 2022, 27, x FOR PEER REVIEW 13 of Figure 11. Experimental and theoretical cross-sections for the D2O + normal-H2 → D2O(202) + norm H2 inelastic collisions in the conditions (Table 1) of D2O seeded in He and normal-H2 at 145 K. Bla squares are experimental ICSs for the 202 rotational level, with error bars corresponding to statisti uncertainties at 95% of the confidence interval. In red solid line, the total theory contributions a convoluted to experimental resolution, with partial contributions in violet solid line, the 000 → transition, in olive, the 111 → 202 transition, in magenta dotted line, the 101 → 212 transition, and gray 110 → 212 transition, assuming the initial rotational population and a partial contribution fro the 212 state, when probing the 202 state, given in Table 2.
The resonances observed for the 000-202 rotational inelastic cross-sections occur b cause of quasi-bound states become transitorily accessible during the collision. They we characterized by analyzing partial cross-sections differing by their total angular mome tum J [8]. The total cross-section (summed over all partial waves J) was found to show, resonance, large contributions from specific orbital angular momentum l, the rotation states of H2, and D2O being jD2O = 0 and jH2 = 0. An analysis of the resonances was p forward in the preceding paper [8].
The experiments conducted with normal-H2 at a given collision energy can only o serve the summation of all partial wave contributions weighted by the (2J + 1) degenera factor in the ICSs. The contributions of many overlapping partial waves corresponding the energetically allowed J-values tend to average the individual resonance amplitud making them interfere and average out. Unfortunately, for the processes with para-H2, explained previously (see Figures 2 and 3), the theoretical calculations predict too sm transition probabilities to be experimentally probed at the collision energies close to t thresholds, except for the 000 to 202 rotational excitation of D2O [8].
In summary, we were able to observe some resonances in the experimental excitati functions, and the agreement between experiment and theory is excellent in all cases p sented herein, for the collision energy range 10-100 cm −1 .
While the reduced mass of D2O + H2 (1.8313 u) is only slightly different from that H2O + H2 (1.8128 u), the water deuteration can significantly affect inelastic collisional rat the H/D isotopic substitution results in a difference in the reduced mass, rotational co stants, zero-point energies and energies of the quasi-bound states supported by the p tential. Additionally, the center of mass is slightly moved. Consequently, the intensit and positions of the resonant peaks will vary. In Figure 12, we have reported the IC calculated for the first rotational excitation of D2O and H2O by collisions with normal-H  (Table 1) of D 2 O seeded in He and normal-H 2 at 145 K. Black squares are experimental ICSs for the 2 02 rotational level, with error bars corresponding to statistical uncertainties at 95% of the confidence interval. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions in violet solid line, the 0 00 → 2 02 transition, in olive, the 1 11 → 2 02 transition, in magenta dotted line, the 1 01 → 2 12 transition, and in gray 1 10 → 2 12 transition, assuming the initial rotational population and a partial contribution from the 2 12 state, when probing the 2 02 state, given in Table 2.
The experiments conducted with normal-H 2 at a given collision energy can only observe the summation of all partial wave contributions weighted by the (2J + 1) degeneracy factor in the ICSs. The contributions of many overlapping partial waves corresponding to the energetically allowed J-values tend to average the individual resonance amplitudes, making them interfere and average out. Unfortunately, for the processes with para-H 2 , as explained previously (see Figures 2 and 3), the theoretical calculations predict too small transition probabilities to be experimentally probed at the collision energies close to the thresholds, except for the 0 00 to 2 02 rotational excitation of D 2 O [8].
In summary, we were able to observe some resonances in the experimental excitation functions, and the agreement between experiment and theory is excellent in all cases presented herein, for the collision energy range 10-100 cm −1 .
While the reduced mass of D 2 O + H 2 (1.8313 u) is only slightly different from that of H 2 O + H 2 (1.8128 u), the water deuteration can significantly affect inelastic collisional rates: the H/D isotopic substitution results in a difference in the reduced mass, rotational constants, zero-point energies and energies of the quasi-bound states supported by the potential. Additionally, the center of mass is slightly moved. Consequently, the intensities and positions of the resonant peaks will vary. In Figure 12, we have reported the ICSs calculated for the first rotational excitation of D 2 O and H 2 O by collisions with normal-H 2 .
To ease the comparison, the rotational energy diagram of H 2 O and D 2 O was added. Note that ortho and para systems are inverted due to the nuclear spins (deuteron versus proton). We can observe in Figure 12 that the first rotational transitions 000 → 111 and 101 → 110 for H2O + H2 and D2O + H2 inelastic collisions. Just after the threshold rise and before the decrease, the excitation ICSs for H2O tend to be smaller than those for D2O by typically 30% (in qualitative accordance with the energy gap law) while the propensity rules (discussed in the theoretical results section in terms of the dipolar transitions in the C2v symmetry) are equivalent for these b-type transitions (∆ka = ±1, ∆kc = ±1). At larger collision energies, the differences between D2O and H2O tend to be smaller. Thus, for these two dipolar transitions, the differences between D2O and H2O are not substantial and mainly reflect the kinematics effects (mass and velocities) due to the change in rotational constants and the corresponding different rotational thresholds. For dipole-forbidden transitions, however, much larger differences are observed, as observed in Figure 12a for the transition 000 → 202 where the cross-section peak for D2O exceeds the cross-section plateau for H2O by a factor of ~ 2.5. As explained by Scribano et al. [20], this is due to the nondipolar interaction terms in the PES, which are more affected by the isotopic substitution than the dipolar terms. As a result, the H/D substitution results both in kinematics and PES effect whose relative importance depends on the type of transition.

Theoretical Calculations
As in Bergeat et al. [8], inelastic scattering calculations were performed at the closecoupling level with the MOLSCAT program [32] using a rigid-rotor D2O-H2 potential energy surface (PES) adapted from the full-dimensional H2O-H2 PES of Valiron et al. [22]. The shift of center of mass and change of internal geometries (zero-point vibrational effects) were both taken into account. Full details about the geometrical transformation from H2O-H2 to D2O-H2 can be found in [20], but it should be noted that we found a small error in the implementation of Scribano et al. [20] so that in practice a new rigid-rotor D2O-H2 PES was used in [8] and also herein.
D2O is an asymmetric top, and the rotational constants were taken as A = 15.41998, B = 7.27299, and C = 4.84529 cm −1 . Centrifugal correction terms were unnecessary for the low-lying levels examined herein. For H2, the rotational constant is B0 = 59.322 cm −1 , as in We can observe in Figure 12 that the first rotational transitions 0 00 → 1 11 and 1 01 → 1 10 for H 2 O + H 2 and D 2 O + H 2 inelastic collisions. Just after the threshold rise and before the decrease, the excitation ICSs for H 2 O tend to be smaller than those for D 2 O by typically 30% (in qualitative accordance with the energy gap law) while the propensity rules (discussed in the theoretical results section in terms of the dipolar transitions in the C 2v symmetry) are equivalent for these b-type transitions (∆k a = ±1, ∆k c = ±1). At larger collision energies, the differences between D 2 O and H 2 O tend to be smaller. Thus, for these two dipolar transitions, the differences between D 2 O and H 2 O are not substantial and mainly reflect the kinematics effects (mass and velocities) due to the change in rotational constants and the corresponding different rotational thresholds. For dipole-forbidden transitions, however, much larger differences are observed, as observed in Figure 12a for the transition 0 00 → 2 02 where the cross-section peak for D 2 O exceeds the cross-section plateau for H 2 O by a factor of~2.5. As explained by Scribano et al. [20], this is due to the non-dipolar interaction terms in the PES, which are more affected by the isotopic substitution than the dipolar terms. As a result, the H/D substitution results both in kinematics and PES effect whose relative importance depends on the type of transition.

Theoretical Calculations
As in Bergeat et al. [8], inelastic scattering calculations were performed at the closecoupling level with the MOLSCAT program [32] using a rigid-rotor D 2 O-H 2 potential energy surface (PES) adapted from the full-dimensional H 2 O-H 2 PES of Valiron et al. [22]. The shift of center of mass and change of internal geometries (zero-point vibrational effects) were both taken into account. Full details about the geometrical transformation from H 2 O-H 2 to D 2 O-H 2 can be found in [20], but it should be noted that we found a small error in the implementation of Scribano et al. [20] so that in practice a new rigid-rotor D 2 O-H 2 PES was used in [8] and also herein. D 2 O is an asymmetric top, and the rotational constants were taken as A = 15.41998, B = 7.27299, and C = 4.84529 cm −1 . Centrifugal correction terms were unnecessary for the low-lying levels examined herein. For H 2 , the rotational constant is B 0 = 59.322 cm −1 , as in Bergeat et al. [8]. The reduced mass of the D 2 O-H 2 system is 1.8313 u. The coupled differential scattering equations were solved using the hybrid modified log-derivative Airy propagator for total energies up to 300 cm −1 .
Convergence was carefully checked with respect to the propagator step size. Finally, we used a fine energy grid of 0.25 cm −1 , and the cross-sections were interpolated with a cubic spline in order to convolve them with the experimental collision energy spread ( 1 cm −1 ).

Experimental Procedure
We performed our experiments with a crossed-beam apparatus (Figure 2), with a fixed water beam and a rotatable H 2 beam.
The gas manifold comprises a 0.5 L steel tank of high purity liquid deuterated water into which pure neon or helium is added at 12 bar. The tank is held at 40 • C, resulting in a 0.8578 kPa D 2 O vapor pressure [33] and 0.07% D 2 O/Ne or He gas mixture. This gas mixture is then introduced into the stagnation region of a pulsed Even-Lavie valve (100-µm aperture, 150-µm conical nozzle,~10 µs pulse-time, 10 Hz repetition rate) maintained at 320 K to avoid any condensation of water before the supersonic expansion. A turbomolecular pump keeps an average chamber pressure of ∼10 −6 mbar under full pulsed valve operation. The resulting D 2 O beam seeded in Ne or He travels ∼95 mm through the source chamber, separated to the main chamber by a 3-mm inner diameter skimmer. Then, the supersonic molecular beam flies~22 cm to the interaction region (with a background pressure < 10 −6 mbar), where the second supersonic beam crosses, in perfect synchronization.
Normal-H 2 beam pulses were obtained using cryogenically cooled Even-Lavie valve (pulse-time of 7.3 and 5.8 µs for a valve setpoint at 145 or 10 K, respectively) with 12.5-13.5 bar of net hydrogen in the reservoir. The pulsed valve chamber is maintained at ∼10 −6 mbar by a turbomolecular pump and is separated to the main chamber by a 2-mm skimmer 33 mm after the nozzle. This H 2 source chamber may rotate from 90 • to 13 • allowing the beam intersection angle to be varied continuously. The distance from the skimmer to the crossing point is~6 cm and the time delay between the two pulsed beams was adjusted to ensure a perfect overlap of both beams at the crossing point [24]. To minimize background contributions from any rotationally excited molecules in the water beam or the main chamber, the H 2 beam is triggered at 5 Hz, with signals averaged in alternating pulse mode and recorded as two intensities: signal and background.
Detection of D 2 O rotational levels in its (X 1 A 1 , v = 0) ground vibronic state is achieved in the beam-crossing region by a laser beam perpendicular to the scattering plane. D 2 O is detected by REMPI with a 2 + 1 scheme in the UV: two-photon resonant and one-photon ionization at~247 nm via the (C 1 B 1 , v' = 0) Rydberg state. The laser system consisted of a dye laser (Coumarin 500) pumped by the third harmonic at 355 nm of an Nd:YAG laser. UV laser pulse energies of 2.5-4.5 mJ per pulse at repetition rate of 10 Hz were generated by doubling the visible laser output in a beta barium borate (BBO) crystal. A small part of the visible dye laser was sent to a wavemeter to control the D 2 O rotational state probed, and a photodiode was monitoring the UV laser to ensure that the energy remained constant during scanning. A two-stage time-of-flight mass spectrometer (TOF MS), which has specially cut plates to allow the beams to pass through, is positioned in the beam scattering plane and tilted by 135 • to the water supersonic beam. A multichannel plate detector is used for ion detection. The signal is pre-amplified, gated, and integrated using a boxcar over the time-of-flight corresponding to the D 2 O mass and recorded to a computer.
To measure velocity distribution, the density temporal profiles of D 2 O by REMPI-MS detection at the crossing point and of Ne or He by a fast ionization gauge (FIG) at d = 393.3 mm downstream were recorded and fitted to Gaussian profiles. The detailed procedure was used to extract the density pulse profile, and the velocity distribution was detailed previously [24,34]. Briefly, the Gaussian functions with peak positions t 0 and t 1 and widths γ 0 and γ 1 yield the peak velocities v = d/(t 1 − t 0 ), and the velocity spreads from the pulse broadening, taking into account the response time of the fast-ionization gauge. The results are reported in Table 1. To measure the angular divergence of the beams, the TOF-MS was replaced by a second FIG inserted perpendicular to the beam axis. The two density profiles obtained by the two FIGs separated by d = 393.3 mm were measured at different positions of the FIGS perpendicularly to the beam axis: Gaussian fits gave the angular spread reported in Table 1 and the exact beam axis locations. It should be noticed that the beam velocity characterizations obtained by the two FIGs or by the association of the REMPI signal with the far FIG gave the same results.
ICSs were obtained from the REMPI signal intensities of the rotationally excited water molecules, D 2 O(j kakc ). The intensities, signal, and background were acquired for each individual laser shot and intersection angle. The points associated with too small UV laser intensities were rejected to stay in the linear regime whereby the REMPI intensity only depends on the water density. After subtraction of the background mean for each angle, each signal intensity was normalized by the average value over the whole angle scan, allowing for accumulation of several short-duration scans measured in 2 or 3 days. The conditions for each result set presented are summarized in Table 3. All the plotted vertical error bars on the experimental ICSs represent the combination of statistical fluctuations at a 95% confidence interval of the signal and background data. The plotted error bars on energy are estimated from velocity and crossing angle uncertainties: it represents the uncertainty of the peak collision energy value, not the collision energy spread that is reported in Figure 6b. Table 3. Each experimental ICSs of a specific rotationally excited water molecules, D 2 O(j kakc ) corresponds to a number of laser shots, n, per angle, scanning the beam intersection angle in the range The experimental ICSs were obtained from the averaged REMPI signal intensities I, the relative velocity v r of the D 2 O and H 2 beams, and the mean interaction time ∆t between the two beam pulses, as σ = I/(v r ∆t). This interaction time, given in Figure 6a, takes full account of the density-to-flux transformation under our working conditions, when modelling pulsed beam densities with spatial and temporal Gaussian distribution functions.
Indeed, the beam crossing region where collisions contribute to the signal is delimited by the waist of the laser beam steered perpendicularly to the beam scattering plane and also by the beam overlap region even far upstream from the crossing point (see Figure 13). In fact, we consider the locus of points (x, y) for which the product of pulse beam densities n H2 (x, y, t) x n D2O (x, y, t) = 1/e. As at low collision energies (or angles), which correspond to near-threshold collisions, the excited water molecules have almost no recoil energies (or velocities); the scattered D 2 O are therefore travelling in phase with the center-of-mass (cm) velocity vector. For each collision energy, there is only one definite value for the velocity of the recoiling product, the rotationally excited D 2 O, in the center-of-mass frame (Figure 2b). Under all our experimental conditions, the scattered D 2 O product recoils with a velocity in the cm frame, which remains rather low (in the most unfavorable case: 82 m s −1 for a v cm = 809 m s −1 ); its velocity in the laboratory frame thus remains very close to the cm velocity. This induces not only an angular spread and a collision energy spread (Figure 6b) but also a highest value of the mean collision energy at a specific angle relative to the beam axes, γ due to a shift of the mean value of crossing angle [34]. In the ∆t calculations, the angular distribution of scattered species was also taken into account. However, in our previous study on inelastic collisions between D 2 O and para-H 2 [8], the differential cross-sections were also calculated, and it was shown that the collisions were backward, forward, backward-forward, or istropic: the dynamics were varying with the collision energy. As the angular distribution (forward or isotropic) was found to have an almost negligible impact on the time-interaction values, we have simplified the analysis, considering an isotropic distribution of the products at all collision energies for all the rotationally excited D 2 O products. ∆t calculations, the angular distribution of scattered species was also taken into However, in our previous study on inelastic collisions between D2O and para-H differential cross-sections were also calculated, and it was shown that the collis backward, forward, backward-forward, or istropic: the dynamics were varying collision energy. As the angular distribution (forward or isotropic) was found t almost negligible impact on the time-interaction values, we have simplified the considering an isotropic distribution of the products at all collision energies f rotationally excited D2O products. Figure 13. The interaction volumes taken into account in the density-to-flux transformat ferent crossing angles relatively to the water beam: in black 13°, in violet 21°, in blue 31° 42°, and in green 70°. The red dashed line is the D2O:Ne beam axis. The black and gre lines are the H2 beam axis, and the black and green solid lines are the center-of-mass veloc for a crossing angle of 13° and 70°, respectively.

Data Treatment for Comparisons
The total inelastic collision cross-section σi→f is related to the overall proba an inelastic transition will take place from the rotational state |i> of D2O to the st collisions with normal-H2. The bimolecular rate will induce a transition |i> → | is: where ki→ f (vr) is the rate constant at relative velocity vr, nH2 the number density collision partner, nD2O the D2O number density, and p(i) the population of D2O in tional state |i>. The rate law for the inelastic collisions between D2O and H2 of a particula rotational level |f> density (nf) is thus:

Data Treatment for Comparisons
The total inelastic collision cross-section σ i→f is related to the overall probability that an inelastic transition will take place from the rotational state |i> of D 2 O to the state |f > by collisions with normal-H 2 . The bimolecular rate will induce a transition |i> → |f > in D 2 O is: dn f dt = n H2 n i k i→ f (v r ) = n H2 n D2O p(i)v r σ i→ f (v r ) where k i→f (v r ) is the rate constant at relative velocity v r, n H2 the number density of the H 2 collision partner, n D2O the D 2 O number density, and p(i) the population of D 2 O in the rotational state |i>.
The rate law for the inelastic collisions between D 2 O and H 2 of a particular probed rotational level |f > density (n f ) is thus: where the first term corresponds to filling up the state |f >, and the second term is related to the loss of population out of this initial level f level. σ i→f are the cross-sections from all the populated levels |i> to the |f > state when the relative velocity between D 2 O, and H 2 is equal to v r , p(i) is the population ratio in the state |i>, and σ f→l the cross-sections of the |f > → |l> transitions if the probed level |f > is initially populated (p(f)). The REMPI intensity is proportional to the D 2 O density of the f rotational level probed, n f . However, as we do not know the exact densities of D 2 O and H 2 in our beams, the experimental ICSs are given in arbitrary units. To compare the experimental ICSs with the theoretical ones, we have thus added or subtracted the calculated cross-sections weighted by the appropriate initial rotational population (see Table 2).
Moreover, at the wavelength of 247.34 or 247.67 nm chosen to probe the 2 02 or 1 11 rotational level of D 2 O, there is a small contribution from the 2 12 level. Thus, the ICSs of the 1 01 → 2 12 and 1 10 → 2 12 were considered with weights taking into account the relative rotational populations and the relative strength of the REMPI transitions from the 2 02 or 1 11 and 2 12 levels (see Table 2).
Each theoretical cross-section is calculated for one collision energy or relative velocity of D 2 O and H 2 (Equation (1)), but the experimental values are resulting from the integration over collision energy spread. To compare theory and experiment, the calculated excitation functions were convoluted with our experimental energy spread (Figure 2), and the experimental data were multiplied by a unique scaling value.