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

Astrid Bergeat, Alexandre Faure, Laurent Wiesenfeld, Chloé Miossec, Sébastien B. Morales, Christian Naulin

► To cite this version:

Astrid Bergeat, Alexandre Faure, Laurent Wiesenfeld, Chloé Miossec, Sébastien B. Morales, et al. Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D2O with Normal-H2. Molecules, 2022, 27 (21), pp.7535. 10.3390/molecules27217535. hal-03875820

HAL Id: hal-03875820

https://hal.science/hal-03875820

Submitted on 17 May 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
Article

# Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D₂O with Normal-H₂

Astrid Bergeat<fnref n="1" />,*, Astrid Faure<fnref n="2" />, Laurent Wiesenfeld<fnref n="3" />, Chloé Miossec<fnref n="1" />, Sébastien B. Morales<fnref n="1" />
and Christian Naulin<fnref n="1" />

**Abstract:** We present a combined experimental and theoretical study on the rotationally inelastic scattering of heavy water, D₂O, with normal-H₂. Crossed-molecular beam measurements are performed in the collision energy range between 10 and 100 cm⁻<fnref n="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 D₂O using the REMPI technique. These measurements are complemented by quantum close-coupling scattering calculations based on a high-accuracy D₂O-H₂ 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 (D₂O versus H₂O) are also discussed by comparing the shape and magnitude of the respective cross-sections.

**Keywords:** inelastic collisions; cross-sections; water; hydrogen

<figure><img src="image_1.png" /></figure>

check for
updates

**Citation:** Bergeat, A.; Faure, A.; Wiesenfeld, L.; Miossec, C.; Morales, S.B.; Naulin, C. Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D₂O with Normal-H₂. *Molecules* **2022**, *27*, 7535. https://doi.org/10.3390/molecules27217535

Academic Editor: Niels Engholm Henriksen

Received: 21 September 2022

Accepted: 29 October 2022

Published: 3 November 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_2.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. 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⁻⁵, 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₂O are of special importance because they can help to understand the origin of water in protostars [3].

D₂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₂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₂ 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

<footnote n="1">¹ Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France</footnote>

<footnote n="2">² Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France</footnote>

<footnote n="3">³ Laboratoire Aimé-Cotton, CNRS, Université Paris-Saclay, Bât 505, F-91405 Orsay, France</footnote>

<footnote>* Correspondence: astrid.bergeat@u-bordeaux.fr</footnote>
Article

# Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D₂O with Normal-H₂

Astrid Bergeat<fnref n="1" />,*, Astrid Faure<fnref n="2" />, Laurent Wiesenfeld<fnref n="3" />, Chloé Miossec<fnref n="1" />, Sébastien B. Morales<fnref n="1" />
and Christian Naulin<fnref n="1" />

**Abstract:** We present a combined experimental and theoretical study on the rotationally inelastic scattering of heavy water, D₂O, with normal-H₂. Crossed-molecular beam measurements are performed in the collision energy range between 10 and 100 cm⁻<fnref n="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 D₂O using the REMPI technique. These measurements are complemented by quantum close-coupling scattering calculations based on a high-accuracy D₂O-H₂ 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 (D₂O versus H₂O) are also discussed by comparing the shape and magnitude of the respective cross-sections.

**Keywords:** inelastic collisions; cross-sections; water; hydrogen

<figure><img src="image_2.png" /></figure>

check for
updates

**Citation:** Bergeat, A.; Faure, A.; Wiesenfeld, L.; Miossec, C.; Morales, S.B.; Naulin, C. Near-Threshold and Resonance Effects in Rotationally Inelastic Scattering of D₂O with Normal-H₂. *Molecules* **2022**, *27*, 7535. https://doi.org/10.3390/molecules27217535

Academic Editor: Niels Engholm Henriksen

Received: 21 September 2022

Accepted: 29 October 2022

Published: 3 November 2022

**Publisher's Note:** MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

<figure><img src="image_3.png" /></figure>

**Copyright:** © 2022 by the authors. Licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution (CC BY) license (https://creativecommons.org/licenses/by/4.0/).

## 1. 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⁻⁵, 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₂O are of special importance because they can help to understand the origin of water in protostars [3].

D₂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₂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₂ 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

<footnote n="1">¹ Univ. Bordeaux, CNRS, Bordeaux INP, ISM, UMR 5255, F-33400 Talence, France</footnote>

<footnote n="2">² Univ. Grenoble Alpes, CNRS, IPAG, F-38000 Grenoble, France</footnote>

<footnote n="3">³ Laboratoire Aimé-Cotton, CNRS, Université Paris-Saclay, Bât 505, F-91405 Orsay, France</footnote>

<footnote>* Correspondence: astrid.bergeat@u-bordeaux.fr</footnote>
or research fields (thermal conductivity, transport properties for water diluted in molecular hydrogen, or the dynamics of molecular hydrogen confined in the cages of clathrate hydrates) and allow to better understand the non-covalent interactions between molecules of great importance for energy applications, combustion chemistry, and astrochemistry. The efficiency of the inelastic process (i.e., whether a collision does or does not promote the molecule in a given state into another state) is governed by the state-to-state integral cross-section (ICS, σ), which is a function of molecular energy state, colliding partner, and relative translational energies. Excitation requires that the collision energy of the impact be sufficient to reach the upper state, which at sufficiently low temperature becomes energetically impossible. It should be noticed that the rate coefficients are calculated by averaging the product σ × v<sub>r</sub> over a thermal distribution of relative velocities v<sub>r</sub>.

If inelastic scattering ICS is zero below the threshold energy (the energy needed to
promote the molecule from its initial state to an upper final state), above this threshold,
classical mechanics predicts that it rises sharply to a maximum and then decreases smoothly
at higher energies [6]. However, quantum mechanics show that at low kinetic energies,
the inelastic scattering ICS of simple molecular species colliding with H<sub>2</sub> do not follow
such a simple threshold law but are instead highly structured, particularly in the vicinity
of the thresholds of the lowest molecular rotational excitations. Such behavior of the ICS
characterizes the interaction potential between the colliding partners and is very sensitive
to the shape of this potential energy surface (PES).

The rotational energy levels of D₂O are labelled by the quantum number *j*<sub>D2O</sub> and the two pseudo-quantum numbers *k*<sub>a</sub> and *k*<sub>c</sub>, which are the projections of *j*<sub>D2O</sub> along the principal inertia axis, *a*- and *c*-axes, respectively. The C<sub>2</sub> symmetry axis corresponds to the *b*-axis, the *a*-axis being in-plane and perpendicular to the *b*-axis. The rotation around the *c*-axis is in the molecular plane. For a given *j*<sub>D2O</sub>, *k*<sub>a</sub> = *j*<sub>D2O</sub> and *k*<sub>c</sub> = *j*<sub>D2O</sub> correspond to rotation around the *a*- and *c*-axes, respectively. *Para*- and *ortho*-D<sub>2</sub>O levels are characterized by *k*<sub>a</sub> + *k*<sub>c</sub> odd and even, respectively, in accordance with the symmetry relations for bosons, due to the nuclear spin of D atoms. As no atom exchange is possible at the studied collision energies and as the *ortho* and *para* levels do not interconvert in inelastic collisions, these are presented separately in Figure 1. The rotational levels of H<sub>2</sub> are labelled by the angular momentum *j*<sub>H2</sub>, with *j*<sub>H2</sub> even for *para*-H<sub>2</sub> and *j*<sub>H2</sub> odd for *ortho*-H<sub>2</sub> due to the fermionic nuclei.

<figure><img src="image_3.png" /><figcaption>**Figure 1.** Rotational energy levels of the water isotopologue, D<sub>2</sub>O, for E < 100 cm<sup>-1</sup> [7]. Levels are labelled by *j*<sub>kakc</sub>. Ortho levels (D<sub>2</sub>O, *k*<sub>a</sub> + *k*<sub>c</sub> even) are depicted in left panel and para levels in right panel.</figcaption></figure>
or research fields (thermal conductivity, transport properties for water diluted in molecular hydrogen, or the dynamics of molecular hydrogen confined in the cages of clathrate hydrates) and allow to better understand the non-covalent interactions between molecules of great importance for energy applications, combustion chemistry, and astrochemistry. The efficiency of the inelastic process (i.e., whether a collision does or does not promote the molecule in a given state into another state) is governed by the state-to-state integral cross-section (ICS, σ), which is a function of molecular energy state, colliding partner, and relative translational energies. Excitation requires that the collision energy of the impact be sufficient to reach the upper state, which at sufficiently low temperature becomes energetically impossible. It should be noticed that the rate coefficients are calculated by averaging the product σ × v<sub>r</sub> over a thermal distribution of relative velocities v<sub>r</sub>.

If inelastic scattering ICS is zero below the threshold energy (the energy needed to
promote the molecule from its initial state to an upper final state), above this threshold,
classical mechanics predicts that it rises sharply to a maximum and then decreases smoothly
at higher energies [6]. However, quantum mechanics show that at low kinetic energies,
the inelastic scattering ICS of simple molecular species colliding with H<sub>2</sub> do not follow
such a simple threshold law but are instead highly structured, particularly in the vicinity
of the thresholds of the lowest molecular rotational excitations. Such behavior of the ICS
characterizes the interaction potential between the colliding partners and is very sensitive
to the shape of this potential energy surface (PES).

The rotational energy levels of D₂O are labelled by the quantum number *j*<sub>D2O</sub> and the two pseudo-quantum numbers *k*<sub>a</sub> and *k*<sub>c</sub>, which are the projections of *j*<sub>D2O</sub> along the principal inertia axis, *a*- and *c*-axes, respectively. The C₂ symmetry axis corresponds to the *b*-axis, the *a*-axis being in-plane and perpendicular to the *b*-axis. The rotation around the *c*-axis is in the molecular plane. For a given *j*<sub>D2O</sub>, *k*<sub>a</sub> = *j*<sub>D2O</sub> and *k*<sub>c</sub> = *j*<sub>D2O</sub> correspond to rotation around the *a*- and *c*-axes, respectively. *Para*- and *ortho*-D₂O levels are characterized by *k*<sub>a</sub> + *k*<sub>c</sub> odd and even, respectively, in accordance with the symmetry relations for bosons, due to the nuclear spin of D atoms. As no atom exchange is possible at the studied collision energies and as the *ortho* and *para* levels do not interconvert in inelastic collisions, these are presented separately in Figure 1. The rotational levels of H₂ are labelled by the angular momentum *j*<sub>H2</sub>, with *j*<sub>H2</sub> even for *para*-H₂ and *j*<sub>H2</sub> odd for *ortho*-H₂ due to the fermionic nuclei.

<figure><img src="image_4.png" /><figcaption>**Figure 1.** Rotational energy levels of the water isotopologue, D₂O, for E < 100 cm⁻¹ [7]. Levels are labelled by *j*<sub>kakc</sub>. Ortho levels (D₂O, *k*<sub>a</sub> + *k*<sub>c</sub> even) are depicted in left panel and para levels in right panel.</figcaption></figure>
Previous theoretical studies of the H₂O–H₂ 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₂O–H₂ and the theoretical calculations and experiments on D₂O–H₂.

In the first molecular beam scattering experiment on H₂O + H₂ elastic differential cross-sections [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₂–H₂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₂O pressure broadening by H₂ [14,15] or inelastic collisions [16,17]. The Perugia group has reported determination of the potential parameters for the isotropic component of the D₂O–D₂ interaction by elastic scattering studies [18]. The group of Nesbitt and coworkers has also studied the infrared spectrum of the H₂O–D₂ complex in the bend mode frequency region [19].

A large number of studies can thus be found in the literature providing information on H₂O–H₂, whereas calculations or experiments on D₂O–H₂ are scarce. The complex has been studied in solid neon matrix [13] where the infrared absorptions principally arose from complexes involving ortho-H₂, for which jH₂ = 1. Excitation rate constants of the D₂O low lying rotational transitions by para-H₂ 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₂O(0₀₀ or 1₀₁) + H₂(jH₂ = 0, 2 or 1) were measured by Sarma et al. [23] at the collision energy of 584 cm⁻¹. At collision energies above 500 cm⁻¹, 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₂O in the 2₀₂ level by collisions with para-H₂ [8], in the very low collision energy range 10–100 cm¹ (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₂O low lying rotational transitions colliding with normal-H₂ 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.

## 2. Results

### 2.1. 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, γ:

$$E_{\text{collision}} = \frac{1}{2}\mu v_r^2 = \frac{1}{2}\mu(v_{\text{D2O}}^2 + v_{\text{H2}}^2 - 2v_{\text{D2O}}v_{\text{H2}}^2 \cos \gamma) \quad (1)$$
Previous theoretical studies of the H₂O–H₂ 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₂O–H₂ and the theoretical calculations and experiments on D₂O–H₂.

In the first molecular beam scattering experiment on H₂O + H₂ elastic differential cross-sections [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₂–H₂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₂O pressure broadening by H₂ [14,15] or inelastic collisions [16,17]. The Perugia group has reported determination of the potential parameters for the isotropic component of the D₂O–D₂ interaction by elastic scattering studies [18]. The group of Nesbitt and coworkers has also studied the infrared spectrum of the H₂O–D₂ complex in the bend mode frequency region [19].

A large number of studies can thus be found in the literature providing information on H₂O–H₂, whereas calculations or experiments on D₂O–H₂ are scarce. The complex has been studied in solid neon matrix [13] where the infrared absorptions principally arose from complexes involving ortho-H₂, for which jH₂ = 1. Excitation rate constants of the D₂O low lying rotational transitions by para-H₂ 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₂O(0₀₀ or 1₀₁) + H₂(jH₂ = 0, 2 or 1) were measured by Sarma et al. [23] at the collision energy of 584 cm⁻¹. At collision energies above 500 cm⁻¹, 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₂O in the 2₀₂ level by collisions with para-H₂ [8], in the very low collision energy range 10–100 cm¹ (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₂O low lying rotational transitions colliding with normal-H₂ 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.

## 2. Results

### 2.1. 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, γ:

$$E_{\text{collision}} = \frac{1}{2}\mu v_r^2 = \frac{1}{2}\mu(v_{\text{D2O}}^2 + v_{\text{H2}}^2 - 2v_{\text{D2O}}v_{\text{H2}}^2 \cos \gamma) \quad (1)$$
where $v_{D2O}$ and $v_{H2}$ are the velocities of the D₂O and H₂ beams, $\mu$ is the reduced mass of the system D₂O-H₂, 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.

<figure><img src="image_4.png" /><figcaption>**Figure 2.** (a) The crossed molecular beam apparatus with variable angle, $\gamma$ from 12.5° to 90°; between the two beam axes generated by pulsed valves (PV); (b) Newton diagrams for D₂O seeded in Ne collisions with normal-H₂ at 13° and 70° and for D₂O seeded in He collisions with normal-H₂ at 35°. The radii of the circles represent the velocities of rotationally excited D₂O in the center-of-mass (c.m.) frame and the green vector, the relative velocities of the colliders.</figcaption></figure>

The D₂O or normal-H₂ 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><caption>**Table 1.** Characteristics of the molecular beams used in the determination of the integral cross-sections.</caption><thead><tr><th>Beam</th><th>Velocity <sup>1</sup><br>/m s<sup>-1</sup></th><th>Velocity Spread <sup>1</sup><br>HWHM/m s<sup>-1</sup></th><th>Pulse Duration <sup>2</sup><br>HWHM/µs</th><th>Angle Spread <sup>3</sup><br>HWHM/°</th></tr></thead><tbody><tr><td>D<sub>2</sub>O in He</td><td>1858 ± 30</td><td>44</td><td>15.0</td><td>0.79</td></tr><tr><td>normal-H<sub>2</sub> at 145 K</td><td>1959 ± 30</td><td>44</td><td>7.9</td><td>1.17</td></tr><tr><td>D<sub>2</sub>O in Ne</td><td>853 ± 24</td><td>14</td><td>19.3</td><td>0.71</td></tr><tr><td>normal-H<sub>2</sub> at 10 K</td><td>947 ± 28</td><td>49</td><td>11.5</td><td>1.17</td></tr></tbody></table>

### 2.1.1. D₂O Beam

Two carrier gases were used for the D₂O beam, namely He and Ne, leading to beams with different characteristics (see Table 1). From the REMPI spectra reported in Figure 3, the D₂O rotational distribution was determined using the spectroscopic data of the C¹B₁, $v' = 0 \leftarrow X¹A₁$, $v = 0$ transitions, from Yang et al. [25] provided in PGopher software [26]. In our previous study on D₂O + para-H₂ [8], a rotational temperature of 12 K was determined for the D₂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

<footnote n="1">¹ Beam velocity peak values and half width at half maximum (HWHM) of the velocity spread deduced from temporal profiles at the crossing point and 393.3 mm downstream. ² Pulse duration (HWHM) at the crossing point. ³ Angular divergence (HWHM).</footnote>
where $v_{D2O}$ and $v_{H2}$ are the velocities of the D₂O and H₂ beams, $\mu$ is the reduced mass of the system D₂O-H₂, 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.

<figure><img src="image_5.png" /><figcaption>**Figure 2.** (a) The crossed molecular beam apparatus with variable angle, $\gamma$ from 12.5° to 90°; between the two beam axes generated by pulsed valves (PV); (b) Newton diagrams for D₂O seeded in Ne collisions with normal-H₂ at 13° and 70° and for D₂O seeded in He collisions with normal-H₂ at 35°. The radii of the circles represent the velocities of rotationally excited D₂O in the center-of-mass (c.m.) frame and the green vector, the relative velocities of the colliders.</figcaption></figure>

The D₂O or normal-H₂ 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><caption>**Table 1.** Characteristics of the molecular beams used in the determination of the integral cross-sections.</caption><thead><tr><th>Beam</th><th>Velocity <sup>1</sup><br>/m s<sup>-1</sup></th><th>Velocity Spread <sup>1</sup><br>HWHM/m s<sup>-1</sup></th><th>Pulse Duration <sup>2</sup><br>HWHM/µs</th><th>Angle Spread <sup>3</sup><br>HWHM/°</th></tr></thead><tbody><tr><td>D<sub>2</sub>O in He</td><td>1858 ± 30</td><td>44</td><td>15.0</td><td>0.79</td></tr><tr><td>normal-H<sub>2</sub> at 145 K</td><td>1959 ± 30</td><td>44</td><td>7.9</td><td>1.17</td></tr><tr><td>D<sub>2</sub>O in Ne</td><td>853 ± 24</td><td>14</td><td>19.3</td><td>0.71</td></tr><tr><td>normal-H<sub>2</sub> at 10 K</td><td>947 ± 28</td><td>49</td><td>11.5</td><td>1.17</td></tr></tbody></table>

### 2.1.1. D₂O Beam

Two carrier gases were used for the D₂O beam, namely He and Ne, leading to beams with different characteristics (see Table 1). From the REMPI spectra reported in Figure 3, the D₂O rotational distribution was determined using the spectroscopic data of the C¹B₁, $v' = 0 \leftarrow X¹A₁$, $v = 0$ transitions, from Yang et al. [25] provided in PGopher software [26]. In our previous study on D₂O + para-H₂ [8], a rotational temperature of 12 K was determined for the D₂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

<footnote n="1">¹ Beam velocity peak values and half width at half maximum (HWHM) of the velocity spread deduced from temporal profiles at the crossing point and 393.3 mm downstream. ² Pulse duration (HWHM) at the crossing point. ³ Angular divergence (HWHM).</footnote>
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.

<figure><img src="image_5.png" /><figcaption>**Figure 3.** A (2 + 1) REMPI spectrum of the C$^1$B$_1$, $v' = 0 \leftarrow X^1$A$_1$, $v = 0$ transition of D$_2$O in the supersonic beam when Ne (top) or He (bottom) is used as a carrier gas. The inset is the enlargement of a part of the spectrum. The labels indicate the rotational levels, $j_{kakc}$, from which the transitions originate. See Yang et al. [25] for more details.</figcaption></figure>

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.
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.

<figure><img src="image_6.png" /><figcaption>**Figure 3.** A (2 + 1) REMPI spectrum of the C$^1$B$_1$, $v' = 0 \leftarrow X^1$A$_1$, $v = 0$ transition of D$_2$O in the supersonic beam when Ne (top) or He (bottom) is used as a carrier gas. The inset is the enlargement of a part of the spectrum. The labels indicate the rotational levels, $j_{kakc}$, from which the transitions originate. See Yang et al. [25] for more details.</figcaption></figure>

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:

$$ \frac{1}{2}mv^2 = \int_{T}^{T_0} C_p dT, \qquad (2) $$

where the temperature of the beam is $T \ll 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 $\sqrt{(20/4)} = 2.2$ for the two carrier gas, which is experimentally observed (see Table 1, where $1858/853 = 2.2$).

### 2.1.2. Normal-H₂ Beam

As the best collisional energy resolution is obtained when the two beams have almost the same velocity (Figure 2), for the H₂ 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₂ beam remains a mixture of *para-* and *ortho*-H₂ with the room temperature ratio, i.e., 3:1. The H₂ beam was probed (Figure 4) when recording ($C^1\Pi_u; v = 2 \leftarrow X^1\Sigma^+_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₂ beam is composed of 25% of $H_2(j_{H2} = 0)$ and 75% of $H_2(j_{H2} = 1)$.

<figure><img src="image_6.png" /><figcaption>**Figure 4.** The (3 + 1) REMPI spectra [29] of normal-H₂ beams. Temperature (145 K or 10 K) is the set point of the cold head. Vertical error bars correspond to statistical uncertainties at 95% of the confidence interval.</figcaption></figure>
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:

$$ \frac{1}{2}mv^2 = \int_{T}^{T_0} C_p dT, \qquad (2) $$

where the temperature of the beam is $T \ll 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 $\sqrt{(20/4)} = 2.2$ for the two carrier gas, which is experimentally observed (see Table 1, where 1858/853 = 2.2).

### 2.1.2. Normal-H₂ Beam

As the best collisional energy resolution is obtained when the two beams have almost the same velocity (Figure 2), for the H₂ 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₂ beam remains a mixture of *para-* and *ortho*-H₂ with the room temperature ratio, i.e., 3:1. The H₂ beam was probed (Figure 4) when recording ($C^1\Pi_u; v = 2 \leftarrow X^1\Sigma^+_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₂ beam is composed of 25% of H₂($j_{H2} = 0$) and 75% of H₂($j_{H2} = 1$).

<figure><img src="image_7.png" /><figcaption>**Figure 4.** The (3 + 1) REMPI spectra [29] of normal-H₂ beams. Temperature (145 K or 10 K) is the set point of the cold head. Vertical error bars correspond to statistical uncertainties at 95% of the confidence interval.</figcaption></figure>
2.1.3. Integral Cross-Sections

Three rotational levels of D₂O were probed: 1₁₁, 2₀₂ for *ortho*-D₂O, and 1₁₀ for *para*-D₂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₂O and H₂ 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⁻¹ for the 0₀₀–1₁₁ and 1₀₁–1₁₀ 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₂O and Ne–D₂O differ by 30 cm⁻¹, and the D₂O + H₂ resonances are observed at the same collision energies for the two water carrier gases. However, compared to the ICSs measured when D₂O is seeded in Ne, they are broadened and slightly less marked in the case of D₂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₁₁ rotational population starts decreasing when H₂ collides with D₂O. Indeed, in the D₂O:He beam, the first two rotational levels of *ortho*-D₂O and *para*-D₂O are populated, leading to the de-excitation of the 1₁₁ excited state until the collision energy is large enough to promote the excitation from the ground state (0₀₀ → 1₁₁). When the 1₁₀ 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₀₁ → 1₁₀), not allowing the observation of the de-excitation. In the case of D₂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₀₀ → 1₁₁ and 1₀₁ → 1₁₀ transitions in Figure 5a,b, respectively.

<figure><img src="image_7.png" /><figcaption>**Figure 5.** Experimental cross-sections for different conditions: in orange and in green, for D₂O seeded in Ne and in He, respectively: (a) for the D₂O(0₀₀) → D₂O(1₁₁) transition; (b) for the D₂O(1₀₁) → D₂O(1₁₀) transition. Negative values for the cross-sections mean that the rotational population probed decreases: loss of the population on the 1₁₀ level by transition to the 1₀₁ or 2₁₂ levels and loss of the 1₁₁ population by transition to the 0₀₀ or 2₀₂ levels. Vertical error bars correspond to statistical uncertainties at 95% of the confidence interval.</figcaption></figure>
### 2.1.3. Integral Cross-Sections

Three rotational levels of D₂O were probed: $1_{11}$, $2_{02}$ for *ortho*-D₂O, and $1_{10}$ for *para*-D₂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₂O and H₂ 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⁻¹ for the 0₀₀–1₁₁ and 1₀₁–1₁₀ 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₂O and Ne–D₂O differ by 30 cm⁻¹, and the D₂O + H₂ resonances are observed at the same collision energies for the two water carrier gases. However, compared to the ICSs measured when D₂O is seeded in Ne, they are broadened and slightly less marked in the case of D₂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₂ collides with D₂O. Indeed, in the D₂O:He beam, the first two rotational levels of *ortho*-D₂O and *para*-D₂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₀₀ → 1₁₁). 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} \rightarrow 1_{10}$), not allowing the observation of the de-excitation. In the case of D₂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} \rightarrow 1_{11}$ and $1_{01} \rightarrow 1_{10}$ transitions in Figure 5a,b, respectively.

<figure><img src="image_8.png" /><figcaption>**Figure 5.** Experimental cross-sections for different conditions: in orange and in green, for D₂O seeded in Ne and in He, respectively: (a) for the D₂O(0₀₀) → D₂O(1₁₁) transition; (b) for the D₂O(1₀₁) → D₂O(1₁₀) transition. Negative values for the cross-sections mean that the rotational population probed decreases: loss of the population on the 1₁₀ level by transition to the 1₀₁ or 2₁₂ levels and loss of the 1₁₁ population by transition to the 0₀₀ or 2₀₂ levels. Vertical error bars correspond to statistical uncertainties at 95% of the confidence interval.</figcaption></figure>
<figure><img src="image_8.png" /><figcaption>**Figure 6.** (a) Interaction time for an isotropic distribution; (b) collision energy spread calculated for the experimental conditions: in orange (green) solid lines, D₂O seeded in Ne(He) for the D₂O(0₀₀) → D₂O(1₁₁) or D₂O(1₀₁) → D₂O(1₀₁) transition, and in purple dashed line, D₂O seeded in He for the D₂O(0₀₀) → D₂O(2₀₂) transition.</figcaption></figure>

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.

## 2.2. *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₂O ground state by collisions with H₂ are plotted in Figures 2 and 3.

Calculations were performed for H₂ in its ground *para* (*j*<sub>H2</sub> = 0) and *ortho* (*j*<sub>H2</sub> = 1) states, and cross-sections for *normal*-H₂ were obtained for an *ortho*-to-*para* ratio fixed at 3. The rotational basis set for D₂O included angular momenta up to *j*<sub>D2O</sub> = 6, corresponding to the lowest 25 *para*- and lowest 25 *ortho*-D₂O levels. The lowest two levels were included in the basis set for *para*-H₂ (*j*<sub>H2</sub> = 0, 2) and *ortho*-H₂ (*j*<sub>H2</sub> = 1, 3). Enlarging the D₂O or H₂ 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*<sub>e</sub> = 238 cm<sup>-1</sup>) and anisotropy of the D₂O–H₂ PES, the inclusion of closed-channels was found to be crucial, even for *ortho*-H₂.

As previously noted for the H₂O–H₂ system [17] and also for other molecule-H₂ systems [30], resonances are more conspicuous for *para*-H₂ (*j*<sub>H2</sub> = 0) than for *ortho*-H₂ (*j*<sub>H2</sub> = 1). This is particularly striking here for the 0₀₀ → 2₀₂ transition (Figure 7) where resonances have magnitudes up to twice that of the background cross-section for *para*-H₂ (*j*<sub>H2</sub> = 0) while for *ortho*-H₂ (*j*<sub>H2</sub> = 1) are much less intense. This can be explained by the non-vanishing quadrupole of *ortho*-H₂ (*j*<sub>H2</sub> = 1) and the corresponding stronger binding with *ortho*-H₂ (*j*<sub>H2</sub> = 1) than with *para*-H₂ (*j*<sub>H2</sub> = 0). This was observed in the H₂O–H₂ complex where the dissociation energy and number of bound-states are indeed significantly larger for *ortho*-H₂ (*j*<sub>H2</sub> = 1) monomers than *para*-H₂ (*j*<sub>H2</sub> = 0) monomers [31]. Overlapping and
<figure><img src="image_9.png" /><figcaption>**Figure 6.** (a) Interaction time for an isotropic distribution; (b) collision energy spread calculated for the experimental conditions: in orange (green) solid lines, D₂O seeded in Ne(He) for the D₂O(0₀₀) → D₂O(1₁₁) or D₂O(1₀₁) → D₂O(1₀₁) transition, and in purple dashed line, D₂O seeded in He for the D₂O(0₀₀) → D₂O(2₀₂) transition.</figcaption></figure>

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.

## 2.2. *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₂O ground state by collisions with H₂ are plotted in Figures 2 and 3.

Calculations were performed for H₂ in its ground *para* (*j*<sub>H2</sub> = 0) and *ortho* (*j*<sub>H2</sub> = 1) states, and cross-sections for *normal*-H₂ were obtained for an *ortho*-to-*para* ratio fixed at 3. The rotational basis set for D₂O included angular momenta up to *j*<sub>D2O</sub> = 6, corresponding to the lowest 25 *para*- and lowest 25 *ortho*-D₂O levels. The lowest two levels were included in the basis set for *para*-H₂ (*j*<sub>H2</sub> = 0, 2) and *ortho*-H₂ (*j*<sub>H2</sub> = 1, 3). Enlarging the D₂O or H₂ 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*<sub>e</sub> = 238 cm⁻¹) and anisotropy of the D₂O–H₂ PES, the inclusion of closed-channels was found to be crucial, even for *ortho*-H₂.

As previously noted for the H₂O–H₂ system [17] and also for other molecule-H₂ systems [30], resonances are more conspicuous for *para*-H₂ (*j*<sub>H2</sub> = 0) than for *ortho*-H₂ (*j*<sub>H2</sub> = 1). This is particularly striking here for the 0₀₀ → 2₀₂ transition (Figure 7) where resonances have magnitudes up to twice that of the background cross-section for *para*-H₂ (*j*<sub>H2</sub> = 0) while for *ortho*-H₂ (*j*<sub>H2</sub> = 1) are much less intense. This can be explained by the non-vanishing quadrupole of *ortho*-H₂ (*j*<sub>H2</sub> = 1) and the corresponding stronger binding with *ortho*-H₂ (*j*<sub>H2</sub> = 1) than with *para*-H₂ (*j*<sub>H2</sub> = 0). This was observed in the H₂O–H₂ complex where the dissociation energy and number of bound-states are indeed significantly larger for *ortho*-H₂ (*j*<sub>H2</sub> = 1) monomers than *para*-H₂ (*j*<sub>H2</sub> = 0) monomers [31]. Overlapping and
broadened resonances are thus expected to be more numerous in D₂O-*ortho*-H₂ (*j*<sub>H2</sub> = 1) with respect to D₂O-*para*-H₂ (*j*<sub>H2</sub> = 0). In addition, it is observed that the background cross-sections for *para*-H₂ (*j*<sub>H2</sub> = 0) and *ortho*-H₂ (*j*<sub>H2</sub> = 1) have similar magnitudes for the transition *0*<sub>00</sub> → *2*<sub>02</sub> while *ortho*-H₂ (*j*<sub>H2</sub> = 1) cross-sections are much larger than *para*-H₂ (*j*<sub>H2</sub> = 0) cross-sections for other transitions, e.g., *0*<sub>00</sub> → *1*<sub>11</sub> and *1*<sub>11</sub> → *2*<sub>02</sub> (Figure 8). This is clearly due to the dipole–quadrupole interaction, which dominates in D₂O dipolar transitions (i.e., those with Δ*j*<sub>D2O</sub> = 0, ±1, Δ*k*<sub>a</sub> = ±1, Δ*k*<sub>c</sub> = ±1) while D₂O quadrupolar transitions (Δ*j*<sub>D2O</sub> = ±2) are much less sensitive to the quadrupole of H₂.

<figure><img src="image_9.png" /><figcaption>**Figure 7.** Rotational excitation cross-sections for the D₂O 0₀₀–2₀₂ transition in blue *para*-H₂(*j*<sub>H2</sub> = 0), in red *ortho*-H₂(*j*<sub>H2</sub> = 1), and in black *normal*-H₂(*j*<sub>H2</sub> = 0 and 1).</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>**Figure 8.** Rotational excitation cross-sections in blue D₂O + *para*-H₂(*j*<sub>H2</sub> = 0), in red D₂O + *ortho*-H₂(*j*<sub>H2</sub> = 1), and in black D₂O + *normal*-H₂(*j*<sub>H2</sub> = 0 and 1): (a) for the 0₀₀–1₁₁ transition; (b) for the 1₁₁–2₀₂ transition.</figcaption></figure>
broadened resonances are thus expected to be more numerous in D₂O-*ortho*-H₂ (*j*<sub>H2</sub> = 1) with respect to D₂O-*para*-H₂ (*j*<sub>H2</sub> = 0). In addition, it is observed that the background cross-sections for *para*-H₂ (*j*<sub>H2</sub> = 0) and *ortho*-H₂ (*j*<sub>H2</sub> = 1) have similar magnitudes for the transition *0*<sub>00</sub> → *2*<sub>02</sub> while *ortho*-H₂ (*j*<sub>H2</sub> = 1) cross-sections are much larger than *para*-H₂ (*j*<sub>H2</sub> = 0) cross-sections for other transitions, e.g., *0*<sub>00</sub> → *1*<sub>11</sub> and *1*<sub>11</sub> → *2*<sub>02</sub> (Figure 8). This is clearly due to the dipole–quadrupole interaction, which dominates in D₂O dipolar transitions (i.e., those with Δ*j*<sub>D2O</sub> = 0, ±1, Δ*k*<sub>a</sub> = ±1, Δ*k*<sub>c</sub> = ±1) while D₂O quadrupolar transitions (Δ*j*<sub>D2O</sub> = ±2) are much less sensitive to the quadrupole of H₂.

<figure><img src="image_10.png" /><figcaption>**Figure 7.** Rotational excitation cross-sections for the D₂O *0*₀−*2*₀₂ transition in blue *para*-H₂(*j*<sub>H2</sub> = 0), in red *ortho*-H₂(*j*<sub>H2</sub> = 1), and in black *normal*-H₂(*j*<sub>H2</sub> = 0 and 1).</figcaption></figure>

<figure><img src="image_11.png" /><figcaption>**Figure 8.** Rotational excitation cross-sections in blue D₂O + *para*-H₂(*j*<sub>H2</sub> = 0), in red D₂O + *ortho*-H₂(*j*<sub>H2</sub> = 1), and in black D₂O + *normal*-H₂(*j*<sub>H2</sub> = 0 and 1): (a) for the *0*₀−*1*₁₁ transition; (b) for the *1*₁₁−*2*₀₂ transition.</figcaption></figure>
### 3. 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₂O and 3:1 for H₂ (see Table 2). In Figure 9a, only one state-to-state contribution is observed, the D₂O(1₀₁) + normal-H₂ → D₂O(1₁₀) + normal-H₂. This is because only the lowest *para-* and *ortho*-D₂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₂O(0₀₀) to D₂O(1₁₁) 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₀₁ to 2₁₂ transition was added to take into account the overlap of the REMPI peak from the probed 1₁₁ rotational level of D₂O with those from the 2₁₂ level. Then, the REMPI intensities for the D₂O(1₀₁) 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₀₁ to 1₀₁ 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₁₁ to 0₀₀ population ratio was then fitted from the 2₀₂ experimental ICSs (Figure 11), using the overlap value of the REMPI peaks from 2₁₂ and 2₀₂ 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₂O(1₁₁) 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.

<table><caption>**Table 2.** Rotational population distributions in the D₂O beam seeded in He or Ne and contribution to the main rotational level probed by another rotational level.</caption><thead><tr><th>Contributions</th><th>D<sub>2</sub>O in He</th><th>D<sub>2</sub>O in Ne</th></tr></thead><tbody><tr><td>ortho-D<sub>2</sub>O:para-D<sub>2</sub>O <sup>1</sup></td><td>2/3:1/3</td><td>2/3:1/3</td></tr><tr><td>0<sub>00</sub>:1<sub>11</sub> <sup>2</sup></td><td>90:10</td><td>100:0</td></tr><tr><td>1<sub>01</sub>:1<sub>10</sub> <sup>2</sup></td><td>83:17</td><td>100:0</td></tr><tr><td>2<sub>12</sub> in 2<sub>02</sub> <sup>3</sup></td><td>0.20</td><td></td></tr><tr><td>2<sub>12</sub> in 1<sub>11</sub> <sup>3</sup></td><td>0.12</td><td>0.12</td></tr></tbody></table>

The resonances observed for the 0₀₀–2₀₂ 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₂, and D₂O being *j*<sub>D₂O</sub> = 0 and *j*<sub>H₂</sub> = 0. An analysis of the resonances was put forward in the preceding paper [8].

<footnote n="1">¹ The *ortho-to-para* ratio of D₂O remains at the room temperature equilibrium ratio of 2. ² The relative population ratio of the two *ortho-* or *para*-D₂O first rotational levels. ³ The relative contribution of the 2₁₂ rotational level to the signal of the 2₀₂ or 1₁₁ level due to an overlap of the lines in the REMPI spectra (see the text for details).</footnote>
### 3. 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₂O and 3:1 for H₂ (see Table 2). In Figure 9a, only one state-to-state contribution is observed, the D₂O(1₀₁) + normal-H₂ → D₂O(1₁₀) + normal-H₂. This is because only the lowest *para-* and *ortho*-D₂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₂O(0₀₀) to D₂O(1₁₁) 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₀₁ to 2₁₂ transition was added to take into account the overlap of the REMPI peak from the probed 1₁₁ rotational level of D₂O with those from the 2₁₂ level. Then, the REMPI intensities for the D₂O(1₀₁) 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₀₁ to 1₀₁ 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₁₁ to 0₀₀ population ratio was then fitted from the 2₀₂ experimental ICSs (Figure 11), using the overlap value of the REMPI peaks from 2₁₂ and 2₀₂ 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₂O(1₁₁) 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.

<table><caption>**Table 2.** Rotational population distributions in the D₂O beam seeded in He or Ne and contribution to the main rotational level probed by another rotational level.</caption><thead><tr><th>Contributions</th><th>D<sub>2</sub>O in He</th><th>D<sub>2</sub>O in Ne</th></tr></thead><tbody><tr><td>ortho-D<sub>2</sub>O:para-D<sub>2</sub>O <sup>1</sup></td><td>2/3:1/3</td><td>2/3:1/3</td></tr><tr><td>0<sub>00</sub>:1<sub>11</sub> <sup>2</sup></td><td>90:10</td><td>100:0</td></tr><tr><td>1<sub>01</sub>:1<sub>10</sub> <sup>2</sup></td><td>83:17</td><td>100:0</td></tr><tr><td>2<sub>12</sub> in 2<sub>02</sub> <sup>3</sup></td><td>0.20</td><td></td></tr><tr><td>2<sub>12</sub> in 1<sub>11</sub> <sup>3</sup></td><td>0.12</td><td>0.12</td></tr></tbody></table>

The resonances observed for the 0₀₀–2₀₂ 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₂, and D₂O being *j*<sub>D₂O</sub> = 0 and *j*<sub>H₂</sub> = 0. An analysis of the resonances was put forward in the preceding paper [8].

<footnote n="1">¹ The *ortho-to-para* ratio of D₂O remains at the room temperature equilibrium ratio of 2. ² The relative population ratio of the two *ortho-* or *para*-D₂O first rotational levels. ³ The relative contribution of the 2₁₂ rotational level to the signal of the 2₀₂ or 1₁₁ level due to an overlap of the lines in the REMPI spectra (see the text for details).</footnote>
<figure><img src="image_11.png" /><figcaption>**Figure 9.** Experimental and theoretical cross-sections for the D₂O + normal-H₂ → D₂O(l₁₀) + normal-H₂ inelastic collisions. Black squares are experimental ICSs for the l₁₀ rotational level, with error bars corresponding to statistical uncertainties at 95% of the confidence interval, in the conditions of (Table 1): (a) D₂O seeded in Ne and normal-H₂ at 10 K (b) D₂O seeded in He and normal-H₂ at 145 K. Negative values for the cross-sections mean that the l₁₀ rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions in blue, the l₀₁ → l₁₀ transition, in blue dotted line, the loss is due to the l₁₀ → l₀₁ transition, and in gray, the loss is due to the l₁₀ → 2l₂ transition, assuming the initial rotational population given in Table 2.</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>**Figure 10.** Experimental and theoretical cross-sections for the D₂O + normal-H₂ → D₂O(l₁₁) + normal-H₂ inelastic collisions. Black squares are experimental ICSs for the l₁₁ rotational level, with error bars corresponding to statistical uncertainties at 95% of the confidence interval, in the conditions (Table 1) of: (a) D₂O seeded in Ne and normal-H₂ at 10 K (b) D₂O seeded in He and normal-H₂ at 145 K. Negative values for the cross-sections mean that the l₁₁ rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions of the 0l₀ → l₁₁ transition (in orange solid line), the l₀₁ → 2l₂ transition (in magenta dotted line), l₁₀ → 2l₂ transition (in gray), the loss due to the l₁₁ → 2l₂ transition (in olive), and the loss due to the l₁₁ → 2l₂ transition (in orange dotted line), assuming the initial rotational population and a partial contribution from the 2l₂ state, when probing the l₁₁ state, given in Table 2.</figcaption></figure>
<figure><img src="image_12.png" /><figcaption>**Figure 9.** Experimental and theoretical cross-sections for the D₂O + normal-H₂ → D₂O(l₁₀) + normal-H₂ inelastic collisions. Black squares are experimental ICSs for the l₁₀ rotational level, with error bars corresponding to statistical uncertainties at 95% of the confidence interval, in the conditions of (Table 1): (a) D₂O seeded in Ne and normal-H₂ at 10 K (b) D₂O seeded in He and normal-H₂ at 145 K. Negative values for the cross-sections mean that the l₁₀ rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions in blue, the l₀₁ → l₁₀ transition, in blue dotted line, the loss is due to the l₁₀ → l₀₁ transition, and in gray, the loss is due to the l₁₀ → 2l₂ transition, assuming the initial rotational population given in Table 2.</figcaption></figure>

<figure><img src="image_13.png" /><figcaption>**Figure 10.** Experimental and theoretical cross-sections for the D₂O + normal-H₂ → D₂O(l₁₁) + normal-H₂ inelastic collisions. Black squares are experimental ICSs for the l₁₁ rotational level, with error bars corresponding to statistical uncertainties at 95% of the confidence interval, in the conditions (Table 1) of: (a) D₂O seeded in Ne and normal-H₂ at 10 K (b) D₂O seeded in He and normal-H₂ at 145 K. Negative values for the cross-sections mean that the l₁₁ rotational population probed decreases. In red solid line, the total theory contributions are convoluted to experimental resolution, with partial contributions of the 0l₀ → l₁₁ transition (in orange solid line), the l₀₁ → 2l₂ transition (in magenta dotted line), l₁₀ → 2l₂ transition (in gray), the loss due to the l₁₁ → 2l₂ transition (in olive), and the loss due to the l₁₁ → 2l₂ transition (in orange dotted line), assuming the initial rotational population and a partial contribution from the 2l₂ state, when probing the l₁₁ state, given in Table 2.</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>**Figure 11.** Experimental and theoretical cross-sections for the D₂O + *normal*-H₂ → D₂O(2₀₂) + *normal*-H₂ inelastic collisions in the conditions (Table 1) of D₂O seeded in He and *normal*-H₂ at 145 K. Black squares are experimental ICSs for the 2₀₂ 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₀₀ → 2₀₂ transition, in olive, the 1₁₁ → 2₀₂ transition, in magenta dotted line, the 1₀₁ → 2₁₂ transition, and in gray 1₁₀ → 2₁₂ transition, assuming the initial rotational population and a partial contribution from the 2₁₂ state, when probing the 2₀₂ state, given in Table 2.</figcaption></figure>

The experiments conducted with *normal*-H₂ 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₂, 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₀₀ to 2₀₂ rotational excitation of D₂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⁻¹.

While the reduced mass of D₂O + H₂ (1.8313 u) is only slightly different from that of H₂O + H₂ (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₂O and H₂O by collisions with *normal*-H₂. To ease the comparison, the rotational energy diagram of H₂O and D₂O was added. Note that *ortho* and *para* systems are inverted due to the nuclear spins (deuteron versus proton).
<figure><img src="image_14.png" /><figcaption>**Figure 11.** Experimental and theoretical cross-sections for the D₂O + *normal*-H₂ → D₂O(2₀₂) + *normal*-H₂ inelastic collisions in the conditions (Table 1) of D₂O seeded in He and *normal*-H₂ at 145 K. Black squares are experimental ICSs for the 2₀₂ 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₀₀ → 2₀₂ transition, in olive, the 1₁₁ → 2₀₂ transition, in magenta dotted line, the 1₀₁ → 2₁₂ transition, and in gray 1₁₀ → 2₁₂ transition, assuming the initial rotational population and a partial contribution from the 2₁₂ state, when probing the 2₀₂ state, given in Table 2.</figcaption></figure>

The experiments conducted with *normal*-H₂ 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₂, 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₀₀ to 2₀₂ rotational excitation of D₂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⁻¹.

While the reduced mass of D₂O + H₂ (1.8313 u) is only slightly different from that of H₂O + H₂ (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₂O and H₂O by collisions with *normal*-H₂. To ease the comparison, the rotational energy diagram of H₂O and D₂O was added. Note that *ortho* and *para* systems are inverted due to the nuclear spins (deuteron versus proton).
<figure><img src="image_14.png" /><figcaption>**Figure 12.** D₂O or H₂O + normal-H₂ rotational inelastic collisions: (a) Theoretical cross-sections of the D₂O 0₀₀ → 1₁₁ transition in orange solid line, the D₂O 0₀₀ → 2₀₂ transition in orange dashed line, the D₂O 1₀₁ → 1₁₀ transition in blue line, the H₂O 0₀₀ → 1₁₁ transition in red solid line, the H₂O 0₀₀ → 2₀₂ transition in red dashed line, and the 1₀₁ → 1₁₀ transition in olive line. (b) Rotational energy diagram up to *j*₁ = 2 with the four transitions shown in panel (a) [7].</figcaption></figure>

We can observe in Figure 12 that the first rotational transitions 0₀₀ → 1₁₁ and 1₀₁ → 1₁₀ for H₂O + H₂ and D₂O + H₂ inelastic collisions. Just after the threshold rise and before the decrease, the excitation ICSs for H₂O tend to be smaller than those for D₂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₂v symmetry) are equivalent for these *b*-type transitions (Δ*k*<sub>a</sub> = ±1, Δ*k*<sub>c</sub> = ±1). At larger collision energies, the differences between D₂O and H₂O tend to be smaller. Thus, for these two dipolar transitions, the differences between D₂O and H₂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₀₀ → 2₀₂ where the cross-section peak for D₂O exceeds the cross-section plateau for H₂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.

**4. Methods**

**4.1. Theoretical Calculations**

As in Bergeat et al. [8], inelastic scattering calculations were performed at the close-coupling level with the MOLSCAT program [32] using a rigid-rotor D₂O–H₂ potential energy surface (PES) adapted from the full-dimensional H₂O–H₂ 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₂O–H₂ to D₂O–H₂ 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₂O–H₂ PES was used in [8] and also herein.

D₂O is an asymmetric top, and the rotational constants were taken as A = 15.41998, B = 7.27299, and C = 4.84529 cm⁻¹. Centrifugal correction terms were unnecessary for the low-lying levels examined herein. For H₂, the rotational constant is B₀ = 59.322 cm⁻¹, as in Bergeat et al. [8]. The reduced mass of the D₂O–H₂ system is 1.8313 u. The coupled
<figure><img src="image_15.png" /><figcaption>**Figure 12.** D₂O or H₂O + normal-H₂ rotational inelastic collisions: (a) Theoretical cross-sections of the D₂O 0₀₀ → 1₁₁ transition in orange solid line, the D₂O 0₀₀ → 2₀₂ transition in orange dashed line, the D₂O 1₀₁ → 1₁₀ transition in blue line, the H₂O 0₀₀ → 1₁₁ transition in red solid line, the H₂O 0₀₀ → 2₀₂ transition in red dashed line, and the 1₀₁ → 1₁₀ transition in olive line. (b) Rotational energy diagram up to j₁ = 2 with the four transitions shown in panel (a) [7].</figcaption></figure>

We can observe in Figure 12 that the first rotational transitions 0₀₀ → 1₁₁ and 1₀₁ → 1₁₀ for H₂O + H₂ and D₂O + H₂ inelastic collisions. Just after the threshold rise and before the decrease, the excitation ICSs for H₂O tend to be smaller than those for D₂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₂v symmetry) are equivalent for these *b*-type transitions (Δkₐ = ±1, Δkₑ = ±1). At larger collision energies, the differences between D₂O and H₂O tend to be smaller. Thus, for these two dipolar transitions, the differences between D₂O and H₂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₀₀ → 2₀₂ where the cross-section peak for D₂O exceeds the cross-section plateau for H₂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.

# 4. Methods

## 4.1. *Theoretical Calculations*

As in Bergeat et al. [8], inelastic scattering calculations were performed at the close-coupling level with the MOLSCAT program [32] using a rigid-rotor D₂O-H₂ potential energy surface (PES) adapted from the full-dimensional H₂O-H₂ 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₂O-H₂ to D₂O-H₂ 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₂O-H₂ PES was used in [8] and also herein.

D₂O is an asymmetric top, and the rotational constants were taken as A = 15.41998, B = 7.27299, and C = 4.84529 cm⁻¹. Centrifugal correction terms were unnecessary for the low-lying levels examined herein. For H₂, the rotational constant is B₀ = 59.322 cm⁻¹, as in Bergeat et al. [8]. The reduced mass of the D₂O-H₂ 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⁻¹.

Convergence was carefully checked with respect to the propagator step size. Finally, we used a fine energy grid of 0.25 cm⁻¹, and the cross-sections were interpolated with a cubic spline in order to convolve them with the experimental collision energy spread (≈ 1 cm⁻¹).

## 4.2. Crossed Molecular Beams

### 4.2.1. Experimental Procedure

We performed our experiments with a crossed-beam apparatus (Figure 2), with a fixed water beam and a rotatable H₂ 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₂O vapor pressure [33] and 0.07% D₂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⁻⁶ mbar under full pulsed valve operation. The resulting D₂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⁻⁶ mbar), where the second supersonic beam crosses, in perfect synchronization.

Normal-H₂ 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⁻⁶ mbar by a turbomolecular pump and is separated to the main chamber by a 2-mm skimmer 33 mm after the nozzle. This H₂ 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₂ beam is triggered at 5 Hz, with signals averaged in alternating pulse mode and recorded as two intensities: signal and background.

Detection of D₂O rotational levels in its (X¹A₁, v = 0) ground vibronic state is achieved in the beam-crossing region by a laser beam perpendicular to the scattering plane. D₂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¹B₁, 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₂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₂O mass and recorded to a computer.

To measure velocity distribution, the density temporal profiles of D₂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*₀ and *t*₁ and widths *γ*₀ and *γ*₁ yield the peak velocities *v* = *d*/*(t*₁ − *t*₀)*, and the velocity spreads
differential scattering equations were solved using the hybrid modified log-derivative Airy propagator for total energies up to 300 cm⁻¹.

Convergence was carefully checked with respect to the propagator step size. Finally, we used a fine energy grid of 0.25 cm⁻¹, and the cross-sections were interpolated with a cubic spline in order to convolve them with the experimental collision energy spread (≈ 1 cm⁻¹).

## 4.2. Crossed Molecular Beams

### 4.2.1. Experimental Procedure

We performed our experiments with a crossed-beam apparatus (Figure 2), with a fixed water beam and a rotatable H₂ 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₂O vapor pressure [33] and 0.07% D₂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⁻⁶ mbar under full pulsed valve operation. The resulting D₂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⁻⁶ mbar), where the second supersonic beam crosses, in perfect synchronization.

Normal-H₂ 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⁻⁶ mbar by a turbomolecular pump and is separated to the main chamber by a 2-mm skimmer 33 mm after the nozzle. This H₂ 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₂ beam is triggered at 5 Hz, with signals averaged in alternating pulse mode and recorded as two intensities: signal and background.

Detection of D₂O rotational levels in its (X¹A₁, *v* = 0) ground vibronic state is achieved in the beam-crossing region by a laser beam perpendicular to the scattering plane. D₂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¹B₁, *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₂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₂O mass and recorded to a computer.

To measure velocity distribution, the density temporal profiles of D₂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*₀ and *t*₁ and widths *γ*₀ and *γ*₁ yield the peak velocities *v* = *d*/*(t*₁ − *t*₀)*, 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₂O(jₖkₚc). 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><caption>**Table 3.** Each experimental ICSs of a specific rotationally excited water molecules, D₂O(jₖkₚc) corresponds to a number of laser shots, *n*, per angle, scanning the beam intersection angle in the range γinitial–γfinal with −1° or −0.5° decrement.</caption><thead><tr><th>D<sub>2</sub>O(j<sub>kakc</sub>)</th><th>Wavelength</th><th>n Laser Shots/Angle</th><th>Intersection Angles</th></tr></thead><tbody><tr><td>2<sub>02</sub></td><td>247.340 nm</td><td>2275</td><td>35°–13° (–0.5°)</td></tr><tr><td>1<sub>11</sub> with D<sub>2</sub>O in Ne</td><td>247.670 nm</td><td>3585</td><td>70°–23° (–1°)</td></tr><tr><td>1<sub>11</sub> with D<sub>2</sub>O in He</td><td>247.670 nm</td><td>3691</td><td>27°–12.5° (–0.5°)</td></tr><tr><td>1<sub>10</sub> with D<sub>2</sub>O in Ne</td><td>247.684 nm</td><td>2482</td><td>60°–13° (–1°)</td></tr><tr><td>1<sub>10</sub> with D<sub>2</sub>O in He</td><td>247.684 nm</td><td>2684</td><td>27°–12.5° (–0.5°)</td></tr></tbody></table>

The experimental ICSs were obtained from the averaged REMPI signal intensities *I*, the relative velocity *v*<sub>r</sub> of the D₂O and H₂ beams, and the mean interaction time Δ*t* between the two beam pulses, as σ = *I*/*v*<sub>r</sub> Δ*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*<sub>H₂</sub>(*x*, *y*, *t*) × *n*<sub>D₂O</sub>(*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₂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₂O, in the center-of-mass frame (Figure 2b). Under all our experimental conditions, the scattered D₂O product recoils with a velocity in the cm frame, which remains rather low (in the most unfavorable case: 82 m s⁻¹ for a *v*<sub>cm</sub> = 809 m s⁻¹); 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₂O and *para*-H₂ [8], the differential cross-sections were also calculated, and it was shown that the collisions were
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₂O(jₖkₚc). 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><caption>**Table 3.** Each experimental ICSs of a specific rotationally excited water molecules, D₂O(jₖkₚc) corresponds to a number of laser shots, *n*, per angle, scanning the beam intersection angle in the range γinitial–γfinal with −1° or −0.5° decrement.</caption><thead><tr><th>D<sub>2</sub>O(j<sub>kakc</sub>)</th><th>Wavelength</th><th>n Laser Shots/Angle</th><th>Intersection Angles</th></tr></thead><tbody><tr><td>2<sub>02</sub></td><td>247.340 nm</td><td>2275</td><td>35°–13° (–0.5°)</td></tr><tr><td>1<sub>11</sub> with D<sub>2</sub>O in Ne</td><td>247.670 nm</td><td>3585</td><td>70°–23° (–1°)</td></tr><tr><td>1<sub>11</sub> with D<sub>2</sub>O in He</td><td>247.670 nm</td><td>3691</td><td>27°–12.5° (–0.5°)</td></tr><tr><td>1<sub>10</sub> with D<sub>2</sub>O in Ne</td><td>247.684 nm</td><td>2482</td><td>60°–13° (–1°)</td></tr><tr><td>1<sub>10</sub> with D<sub>2</sub>O in He</td><td>247.684 nm</td><td>2684</td><td>27°–12.5° (–0.5°)</td></tr></tbody></table>

The experimental ICSs were obtained from the averaged REMPI signal intensities *I*, the relative velocity *v*<sub>r</sub> of the D<sub>2</sub>O and H<sub>2</sub> beams, and the mean interaction time Δ*t* between the two beam pulses, as σ = *I*/*v*<sub>r</sub> Δ*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*<sub>H₂</sub>(x, y, t) × *n*<sub>D₂O</sub>(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₂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₂O, in the center-of-mass frame (Figure 2b). Under all our experimental conditions, the scattered D₂O product recoils with a velocity in the cm frame, which remains rather low (in the most unfavorable case: 82 m s⁻¹ for a *v*<sub>cm</sub> = 809 m s⁻¹); 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₂O and *para*-H₂ [8], the differential cross-sections were also calculated, and it was shown that the collisions were
backward, forward, backward-forward, or isotropic: 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₂O products.

<figure><img src="image_15.png" /><figcaption>**Figure 13.** The interaction volumes taken into account in the density-to-flux transformation for different crossing angles relatively to the water beam: in black 13°, in violet 21°, in blue 31°, in orange 42°, and in green 70°. The red dashed line is the D₂O:Ne beam axis. The black and green dashed lines are the H₂ beam axis, and the black and green solid lines are the center-of-mass velocity vectors, for a crossing angle of 13° and 70°, respectively.</figcaption></figure>

### 4.2.2. Data Treatment for Comparisons

The total inelastic collision cross-section $\sigma_{i \to f}$ is related to the overall probability that an inelastic transition will take place from the rotational state $|i\rangle$ of D₂O to the state $|f\rangle$ by collisions with normal-H₂. The bimolecular rate will induce a transition $|i\rangle \to |f\rangle$ in D₂O is:

$$ \frac{dn_f}{dt} = n_{H2}n_i k_{i \to f}(v_r) = n_{H2}n_{D2O}p(i)v_r\sigma_{i \to f}(v_r) \quad (3) $$

where $k_{i \to f}(v_r)$ is the rate constant at relative velocity $v_r$, $n_{H2}$ the number density of the H₂ collision partner, $n_{D2O}$ the D₂O number density, and $p(i)$ the population of D₂O in the rotational state $|i\rangle$.
backward, forward, backward-forward, or isotropic: 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₂O products.

<figure><img src="image_16.png" /><figcaption>**Figure 13.** The interaction volumes taken into account in the density-to-flux transformation for different crossing angles relatively to the water beam: in black 13°, in violet 21°, in blue 31°, in orange 42°, and in green 70°. The red dashed line is the D₂O:Ne beam axis. The black and green dashed lines are the H₂ beam axis, and the black and green solid lines are the center-of-mass velocity vectors, for a crossing angle of 13° and 70°, respectively.</figcaption></figure>

### 4.2.2. Data Treatment for Comparisons

The total inelastic collision cross-section $\sigma_{i \to f}$ is related to the overall probability that an inelastic transition will take place from the rotational state $|i\rangle$ of D₂O to the state $|f\rangle$ by collisions with normal-H₂. The bimolecular rate will induce a transition $|i\rangle \to |f\rangle$ in D₂O is:

$$ \frac{dn_f}{dt} = n_{H2}n_i k_{i \to f}(v_r) = n_{H2}n_{D2O}p(i)v_r\sigma_{i \to f}(v_r) \quad (3) $$

where $k_{i \to f}(v_r)$ is the rate constant at relative velocity $v_r$, $n_{H2}$ the number density of the H₂ collision partner, $n_{D2O}$ the D₂O number density, and $p(i)$ the population of D₂O in the rotational state $|i\rangle$.
The rate law for the inelastic collisions between D₂O and H₂ of a particular probed rotational level |f> density ($n_f$) is thus:

$$ \frac{dn_f}{dt} = n_{H2}n_{D2O} \sum_{i \neq f} p(i)k_{i \to f}(v_r) - n_{H2}n_{D2O}p(f) \sum_{l \neq f} k_{f \to l}(v_r) \\ = n_{H2}n_{D2O}v_r \left( \sum_{i \neq f} p(i)\sigma_{i \to f}(v_r) - p(f) \sum_{l \neq f} \sigma_{f \to l}(v_r) \right) \quad (4) $$

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. $\sigma_{i \to f}$ are the cross-sections from all the populated levels |i> to the |f> state when the relative velocity between D₂O, and H₂ is equal to $v_r$, $p(i)$ is the population ratio in the state |i>, and $\sigma_{f \to 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₂O density of the f rotational level probed, $n_f$. However, as we do not know the exact densities of D₂O and H₂ 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₀₂ or 1₁₁ rotational level of D₂O, there is a small contribution from the 2₁₂ level. Thus, the ICSs of the 1₀₁ → 2₁₂ and 1₁₀ → 2₁₂ were considered with weights taking into account the relative rotational populations and the relative strength of the REMPI transitions from the 2₀₂ or 1₁₁ and 2₁₂ levels (see Table 2).

Each theoretical cross-section is calculated for one collision energy or relative velocity of D₂O and H₂ (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.

**Author Contributions:** Conceptualization, A.B. and C.N.; laser, S.B.M. and A.B.; experimental investigation, A.B. and C.M. (Master internship); calculations, L.W. and A.F.; funding acquisition, A.B. and A.F.; writing, A.B., A.F. and L.W. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the Programme National « Physique et Chimie du Milieu Interstellaire » (PCMI) of the CNRS/INSU with INC/INP co-funded by CEA and CNES. The experimental research was also funded by Agence Nationale de la Recherche, grant number ANR-20-CE31-0011.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is contained within the article and may be obtained in another format on request from the authors.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

**Sample Availability:** Samples of the compounds are available from the authors.

## References

1. Cheung, A.C.; Thornton, D.D.; Rank, D.M.; Welch, W.J.; Townes, C.H. Detection of Water in Interstellar Regions by its Microwave Radiation. *Nature* **1969**, 221, 626–628. [CrossRef]

2. Van Dishoeck, E.F.; Kristensen, L.E.; Mottram, J.C.; Benz, A.O.; Bergin, E.A.; Caselli, P.; Herpin, F.; Hogerheijde, M.R.; Johnstone, D.; Liseau, R.; et al. Water in star-forming regions: Physics and chemistry from clouds to disks as probed by Herschel spectroscopy. *A&A* **2021**, 648, A24–A57. [CrossRef]

3. Jensen, S.S.; Jorgensen, J.K.; Furuya, K.; Haugbolle, T.; Aikawa, Y. Modeling chemistry during star formation: Water deuteration in dynamic star-forming regions. *A&A* **2021**, 649, A21–A66. [CrossRef]
The rate law for the inelastic collisions between D₂O and H₂ of a particular probed rotational level |f> density ($n_f$) is thus:

$$ \frac{dn_f}{dt} = n_{H2}n_{D2O} \sum_{i \neq f} p(i)k_{i \to f}(v_r) - n_{H2}n_{D2O}p(f) \sum_{l \neq f} k_{f \to l}(v_r) \\ = n_{H2}n_{D2O}v_r \left( \sum_{i \neq f} p(i)\sigma_{i \to f}(v_r) - p(f) \sum_{l \neq f} \sigma_{f \to l}(v_r) \right) \quad (4) $$

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. $\sigma_{i \to f}$ are the cross-sections from all the populated levels |i> to the |f> state when the relative velocity between D₂O, and H₂ is equal to $v_r$, $p(i)$ is the population ratio in the state |i>, and $\sigma_{f \to 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₂O density of the f rotational level probed, $n_f$. However, as we do not know the exact densities of D₂O and H₂ 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₀₂ or 1₁₁ rotational level of D₂O, there is a small contribution from the 2₁₂ level. Thus, the ICSs of the 1₀₁ → 2₁₂ and 1₁₀ → 2₁₂ were considered with weights taking into account the relative rotational populations and the relative strength of the REMPI transitions from the 2₀₂ or 1₁₁ and 2₁₂ levels (see Table 2).

Each theoretical cross-section is calculated for one collision energy or relative velocity of D₂O and H₂ (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.

**Author Contributions:** Conceptualization, A.B. and C.N.; laser, S.B.M. and A.B.; experimental investigation, A.B. and C.M. (Master internship); calculations, L.W. and A.F.; funding acquisition, A.B. and A.F.; writing, A.B., A.F. and L.W. All authors have read and agreed to the published version of the manuscript.

**Funding:** This research was supported by the Programme National « Physique et Chimie du Milieu Interstellaire » (PCMI) of the CNRS/INSU with INC/INP co-funded by CEA and CNES. The experimental research was also funded by Agence Nationale de la Recherche, grant number ANR-20-CE31-0011.

**Institutional Review Board Statement:** Not applicable.

**Informed Consent Statement:** Not applicable.

**Data Availability Statement:** Data is contained within the article and may be obtained in another format on request from the authors.

**Conflicts of Interest:** The authors declare no conflict of interest. The funders had no role in the design of the study; in the collection, analyses, or interpretation of data; in the writing of the manuscript; or in the decision to publish the results.

**Sample Availability:** Samples of the compounds are available from the authors.

## References

1. Cheung, A.C.; Thornton, D.D.; Rank, D.M.; Welch, W.J.; Townes, C.H. Detection of Water in Interstellar Regions by its Microwave Radiation. *Nature* **1969**, 221, 626–628. [CrossRef]

2. Van Dishoeck, E.F.; Kristensen, L.E.; Mottram, J.C.; Benz, A.O.; Bergin, E.A.; Caselli, P.; Herpin, F.; Hogerheijde, M.R.; Johnstone, D.; Liseau, R.; et al. Water in star-forming regions: Physics and chemistry from clouds to disks as probed by Herschel spectroscopy. *A&A* **2021**, 648, A24–A57. [CrossRef]

3. Jensen, S.S.; Jorgensen, J.K.; Furuya, K.; Haugbolle, T.; Aikawa, Y. Modeling chemistry during star formation: Water deuteration in dynamic star-forming regions. *A&A* **2021**, 649, A21–A66. [CrossRef]
4. Butner, H.M.; Charnley, S.B.; Ceccarelli, C.; Rodgers, S.D.; Pardo, J.R.; Parise, B.; Cernicharo, J.; Davis, G.R. Discovery of interstellar heavy water. *Astrophys. J.* **2007**, *659*, L137–L140. [CrossRef]

5. Vastel, C.; Ceccarelli, C.; Caux, E.; Coutens, A.; Cernicharo, J.; Bottinelli, S.; Demyk, K.; Faure, A.; Wiesenfeld, L.; Scribano, Y.; et al. Ortho-to-para ratio of interstellar heavy water. *A&A* **2010**, *521*, L31–L35. [CrossRef]

6. Levine, R.D.; Bernstein, R.B. *Molecular Reaction Dynamics and Chemical Reactivity*; Oxford University Press: New York, NY, USA, 1987.

7. Gordon, I.E.; Rothman, L.S.; Hargreaves, R.J.; Hashemi, R.; Karlovets, E.V.; Skinner, F.M.; Conway, E.K.; Hill, C.; Kochanov, R.V.; Tan, Y.; et al. The HITRAN2020 molecular spectroscopic database. *J. Quant. Spectrosc. Radiat. Transf.* **2022**, *277*, 107949. [CrossRef]

8. Bergeat, A.; Morales, S.B.; Naulin, C.; Wiesenfeld, L.; Faure, A. Probing Low-Energy Resonances in Water-Hydrogen Inelastic Collisions. *Phys. Rev. Lett.* **2020**, *125*, 143402. [CrossRef]

9. Zoltowski, M.; Lique, F.; Karska, A.; Zuchowski, P.S. Rotational excitation of highly excited H₂O by H₂. *MNRAS* **2021**, *502*, 5356–5361. [CrossRef]

10. Bickes, R.W., Jr.; Duquette, G.; van den Meijdenberg, C.J.N.; Rulis, A.M.; Scoles, G.; Smith, K.M. Molecular beam scattering experiments with polar molecules: Measurement of differential collision cross sections for H₂O+H₂, He, Ne, Ar, H₂O and NH₃+H₂, He, NH₃. *J. Phys. B Atom. Mol. Phys.* **1975**, *8*, 3034–3043. [CrossRef]

11. Weida, M.J.; Nesbitt, D.J. High-resolution diode laser study of H₂–H₂O van der Waals complexes: H₂O as proton acceptor and the role of large amplitude motion. *J. Chem. Phys.* **1999**, *110*, 156–167. [CrossRef]

12. Ziemkiewicz, M.P.; Pluetzer, C.; Nesbitt, D.J.; Scribano, Y.; Faure, A.; van der Avoird, A. Overtone vibrational spectroscopy in H₂-H₂O complexes: A combined high level theoretical ab initio, dynamical and experimental study. *J. Chem. Phys.* **2012**, *137*, 084301. [CrossRef] [PubMed]

13. Forney, D.; Jacox, M.E.; Thompson, W.E. Infrared absorptions of the H₂O-H₂ complex trapped in solid neon. *J. Chem. Phys.* **2004**, *121*, 5977–5984. [CrossRef] [PubMed]

14. Dutta, J.M.; Jones, C.R.; Goyette, T.M.; de Lucia, F.C. The Hydrogen and Helium Pressure Broadening at Planetary Temperatures of the 183 and 380 GHz Transitions of Water Vapor. *Icarus* **1993**, *102*, 232–239. [CrossRef]

15. Drouin, B.; Wiesenfeld, L. Low-temperature water–hydrogen-molecule collisions probed by pressure broadening and line shift. *Phys. Rev. A* **2012**, *86*, 022705. [CrossRef]

16. Yang, C.-H.; Sarma, G.; Parker, D.H.; ter Meulen, J.J.; Wiesenfeld, L. State-to-state differential and relative integral cross sections for rotationally inelastic scattering of H₂O by hydrogen. *J. Chem. Phys.* **2011**, *134*, 204308. [CrossRef]

17. Bergeat, A.; Faure, A.; Morales, S.B.; Moudens, A.; Naulin, C. Low-Energy Water–Hydrogen Inelastic Collisions. *J. Phys. Chem. A* **2020**, *124*, 259–264. [CrossRef]

18. Belpassi, L.; Reca, M.L.; Tarantelli, F.; Roncaratti, L.F.; Pirani, F.; Cappelletti, D.; Faure, A.; Scribano, Y. Charge-Transfer Energy in the Water-Hydrogen Molecular Aggregate Revealed by Molecular-Beam Scattering Experiments, Charge Displacement Analysis, and ab Initio Calculations. *J. Am. Chem. Soc.* **2010**, *132*, 13046–13058. [CrossRef]

19. van der Avoird, A.; Scribano, Y.; Faure, A.; Weida, M.J.; Fair, J.R.; Nesbitt, D.J. Intermolecular potential and rovibrational states of the H₂O–D₂ complex. *Chem. Phys.* **2012**, *399*, 28–38. [CrossRef]

20. Scribano, Y.; Faure, A.; Wiesenfeld, L. Communication: Rotational excitation of interstellar heavy water by hydrogen molecules. *J. Chem. Phys.* **2010**, *133*, 231105. [CrossRef]

21. Faure, A.; Wiesenfeld, L.; Scribano, Y.; Ceccarelli, C. Rotational excitation of mono- and doubly-deuterated water by hydrogen molecules. *MNRAS* **2012**, *420*, 699–704. [CrossRef]

22. Valiron, P.; Wernli, M.; Faure, A.; Wiesenfeld, L.; Rist, C.; Kedzuch, S.; Noga, J. R12-calibrated H₂O-H₂ interaction: Full dimensional and vibrationally averaged potential energy surfaces. *J. Chem. Phys.* **2008**, *129*, 134306–14. [CrossRef] [PubMed]

23. Sarma, G.; Yang, C.-H.; Saha, A.K.; Parker, D.H.; Wiesenfeld, L. Rotational excitation of HDO and D₂O by H₂: Experimental and theoretical differential cross-sections. *J. Chem. Phys.* **2013**, *138*, 024314. [CrossRef] [PubMed]

24. Bergeat, A.; Morales, S.B.; Naulin, C.; Klos, J.; Lique, F. Quantum Behavior of Spin-Orbit Inelastic Scattering of C-Atoms by D₂ at Low Energy. *Front. Chem.* **2019**, *7*, 164. [CrossRef] [PubMed]

25. Yang, C.-H.; Sarma, G.; ter Meulen, J.J.; Parker, D.H.; Western, C.M. REMPI spectroscopy and predissociation of the C¹B₁(ν = 0) rotational levels of H₂O, HOD and D₂O. *Phys. Chem. Chem. Phys.* **2010**, *12*, 13983–13991. [CrossRef]

26. Western, C.M. PGOPHER, A Program for Simulating Rotational, Vibrational and Electronic Spectra. *J. Quant. Spectrosc. Radiat. Trans.* **2017**, *186*, 221–242. [CrossRef]

27. Ziemkiewicz, M.P.; Pluetzer, C.; Wojcik, M.; Loreau, J.; van der Avoird, A.; Nesbitt, D.J. Near infrared overtone (ν<sub>OH</sub> = 2 ← 0) spectroscopy of Ne-H₂O clusters. *J. Chem. Phys.* **2017**, *146*, 104204. [CrossRef]

28. Hou, D.; Ma, Y.-T.; Zhang, X.-L.; Li, H. A full-dimension intra- and inter-molecular *ab initio* potential energy surface and predicted infrared spectra for H₂O-He. *J. Mol. Spectrosc.* **2016**, *330*, 217–227. [CrossRef]

29. Pratt, S.T.; Dehmer, P.M.; Dehmer, J.L. Photoionization of excited molecular states. H₂ C¹Π<sub>u</sub>. *Chem. Phys. Lett.* **1984**, *105*, 28. [CrossRef]

30. Ma, Q.; van der Avoird, A.; Loreau, J.; Alexander, M.H.; van de Meerakker, S.Y.T.; Dagdigian, P.J. Resonances in rotationally inelastic scattering of NH₃ and ND₃ with H₂. *J. Chem. Phys.* **2015**, *143*, 044312. [CrossRef]

31. van der Avoird, A.; Nesbitt, D.J. Rovibrational states of the H₂O–H₂ complex: An ab initio calculation. *J. Chem. Phys.* **2011**, *134*, 044314. [CrossRef]
4. Butner, H.M.; Charnley, S.B.; Ceccarelli, C.; Rodgers, S.D.; Pardo, J.R.; Parise, B.; Cernicharo, J.; Davis, G.R. Discovery of interstellar heavy water. *Astrophys. J.* **2007**, *659*, L137–L140. [CrossRef]

5. Vastel, C.; Ceccarelli, C.; Caux, E.; Coutens, A.; Cernicharo, J.; Bottinelli, S.; Demyk, K.; Faure, A.; Wiesenfeld, L.; Scribano, Y.; et al. Ortho-to-para ratio of interstellar heavy water. *A&A* **2010**, *521*, L31–L35. [CrossRef]

6. Levine, R.D.; Bernstein, R.B. *Molecular Reaction Dynamics and Chemical Reactivity*; Oxford University Press: New York, NY, USA, 1987.

7. Gordon, I.E.; Rothman, L.S.; Hargreaves, R.J.; Hashemi, R.; Karlovets, E.V.; Skinner, F.M.; Conway, E.K.; Hill, C.; Kochanov, R.V.; Tan, Y.; et al. The HITRAN2020 molecular spectroscopic database. *J. Quant. Spectrosc. Radiat. Transf.* **2022**, *277*, 107949. [CrossRef]

8. Bergeat, A.; Morales, S.B.; Naulin, C.; Wiesenfeld, L.; Faure, A. Probing Low-Energy Resonances in Water-Hydrogen Inelastic Collisions. *Phys. Rev. Lett.* **2020**, *125*, 143402. [CrossRef]

9. Zoltowski, M.; Lique, F.; Karska, A.; Zuchowski, P.S. Rotational excitation of highly excited H₂O by H₂. *MNRAS* **2021**, *502*, 5356–5361. [CrossRef]

10. Bickes, R.W., Jr.; Duquette, G.; van den Meijdenberg, C.J.N.; Rulis, A.M.; Scoles, G.; Smith, K.M. Molecular beam scattering experiments with polar molecules: Measurement of differential collision cross sections for H₂O+H₂, He, Ne, Ar, H₂O and NH₃+H₂, He, NH₃. *J. Phys. B Atom. Mol. Phys.* **1975**, *8*, 3034–3043. [CrossRef]

11. Weida, M.J.; Nesbitt, D.J. High-resolution diode laser study of H₂–H₂O van der Waals complexes: H₂O as proton acceptor and the role of large amplitude motion. *J. Chem. Phys.* **1999**, *110*, 156–167. [CrossRef]

12. Ziemkiewicz, M.P.; Pluetzer, C.; Nesbitt, D.J.; Scribano, Y.; Faure, A.; van der Avoird, A. Overtone vibrational spectroscopy in H₂-H₂O complexes: A combined high level theoretical ab initio, dynamical and experimental study. *J. Chem. Phys.* **2012**, *137*, 084301. [CrossRef] [PubMed]

13. Forney, D.; Jacox, M.E.; Thompson, W.E. Infrared absorptions of the H₂O-H₂ complex trapped in solid neon. *J. Chem. Phys.* **2004**, *121*, 5977–5984. [CrossRef] [PubMed]

14. Dutta, J.M.; Jones, C.R.; Goyette, T.M.; de Lucia, F.C. The Hydrogen and Helium Pressure Broadening at Planetary Temperatures of the 183 and 380 GHz Transitions of Water Vapor. *Icarus* **1993**, *102*, 232–239. [CrossRef]

15. Drouin, B.; Wiesenfeld, L. Low-temperature water–hydrogen-molecule collisions probed by pressure broadening and line shift. *Phys. Rev. A* **2012**, *86*, 022705. [CrossRef]

16. Yang, C.-H.; Sarma, G.; Parker, D.H.; ter Meulen, J.J.; Wiesenfeld, L. State-to-state differential and relative integral cross sections for rotationally inelastic scattering of H₂O by hydrogen. *J. Chem. Phys.* **2011**, *134*, 204308. [CrossRef]

17. Bergeat, A.; Faure, A.; Morales, S.B.; Moudens, A.; Naulin, C. Low-Energy Water–Hydrogen Inelastic Collisions. *J. Phys. Chem. A* **2020**, *124*, 259–264. [CrossRef]

18. Belpassi, L.; Reca, M.L.; Tarantelli, F.; Roncaratti, L.F.; Pirani, F.; Cappelletti, D.; Faure, A.; Scribano, Y. Charge-Transfer Energy in the Water-Hydrogen Molecular Aggregate Revealed by Molecular-Beam Scattering Experiments, Charge Displacement Analysis, and ab Initio Calculations. *J. Am. Chem. Soc.* **2010**, *132*, 13046–13058. [CrossRef]

19. van der Avoird, A.; Scribano, Y.; Faure, A.; Weida, M.J.; Fair, J.R.; Nesbitt, D.J. Intermolecular potential and rovibrational states of the H₂O–D₂ complex. *Chem. Phys.* **2012**, *399*, 28–38. [CrossRef]

20. Scribano, Y.; Faure, A.; Wiesenfeld, L. Communication: Rotational excitation of interstellar heavy water by hydrogen molecules. *J. Chem. Phys.* **2010**, *133*, 231105. [CrossRef]

21. Faure, A.; Wiesenfeld, L.; Scribano, Y.; Ceccarelli, C. Rotational excitation of mono- and doubly-deuterated water by hydrogen molecules. *MNRAS* **2012**, *420*, 699–704. [CrossRef]

22. Valiron, P.; Wernli, M.; Faure, A.; Wiesenfeld, L.; Rist, C.; Kedzuch, S.; Noga, J. R12-calibrated H₂O-H₂ interaction: Full dimensional and vibrationally averaged potential energy surfaces. *J. Chem. Phys.* **2008**, *129*, 134306–14. [CrossRef] [PubMed]

23. Sarma, G.; Yang, C.-H.; Saha, A.K.; Parker, D.H.; Wiesenfeld, L. Rotational excitation of HDO and D₂O by H₂: Experimental and theoretical differential cross-sections. *J. Chem. Phys.* **2013**, *138*, 024314. [CrossRef] [PubMed]

24. Bergeat, A.; Morales, S.B.; Naulin, C.; Klos, J.; Lique, F. Quantum Behavior of Spin-Orbit Inelastic Scattering of C-Atoms by D₂ at Low Energy. *Front. Chem.* **2019**, *7*, 164. [CrossRef] [PubMed]

25. Yang, C.-H.; Sarma, G.; ter Meulen, J.J.; Parker, D.H.; Western, C.M. REMPI spectroscopy and predissociation of the C¹B₁(ν = 0) rotational levels of H₂O, HOD and D₂O. *Phys. Chem. Chem. Phys.* **2010**, *12*, 13983–13991. [CrossRef]

26. Western, C.M. PGOPHER, A Program for Simulating Rotational, Vibrational and Electronic Spectra. *J. Quant. Spectrosc. Radiat. Trans.* **2017**, *186*, 221–242. [CrossRef]

27. Ziemkiewicz, M.P.; Pluetzer, C.; Wojcik, M.; Loreau, J.; van der Avoird, A.; Nesbitt, D.J. Near infrared overtone (ν<sub>OH</sub> = 2 ← 0) spectroscopy of Ne-H₂O clusters. *J. Chem. Phys.* **2017**, *146*, 104204. [CrossRef]

28. Hou, D.; Ma, Y.-T.; Zhang, X.-L.; Li, H. A full-dimension intra- and inter-molecular *ab initio* potential energy surface and predicted infrared spectra for H₂O-He. *J. Mol. Spectrosc.* **2016**, *330*, 217–227. [CrossRef]

29. Pratt, S.T.; Dehmer, P.M.; Dehmer, J.L. Photoionization of excited molecular states. H₂ C¹Π<sub>u</sub>. *Chem. Phys. Lett.* **1984**, *105*, 28. [CrossRef]

30. Ma, Q.; van der Avoird, A.; Loreau, J.; Alexander, M.H.; van de Meerakker, S.Y.T.; Dagdigian, P.J. Resonances in rotationally inelastic scattering of NH₃ and ND₃ with H₂. *J. Chem. Phys.* **2015**, *143*, 044312. [CrossRef]

31. van der Avoird, A.; Nesbitt, D.J. Rovibrational states of the H₂O–H₂ complex: An ab initio calculation. *J. Chem. Phys.* **2011**, *134*, 044314. [CrossRef]
32. Hutson, J.M.; Green, S. *MOLSCAT: MOLecular SCATtering*, Version 14; Astrophysics Source Code Library: Houghton, MI, USA, 2012; Available online: https://ui.adsabs.harvard.edu/abs/2012ascl.soft06004H/abstract (accessed on 20 September 2022).

33. Harvey, A.H.; Lemmon, E.W. Correlation for the Vapor Pressure of Heavy Water from the Triple Point to the Critical Point. *J. Phys. Chem. Ref. Data* **2002**, *31*, 173–181. [CrossRef]

34. Naulin, C.; Bergeat, A. Low-energy Scattering in Crossed Molecular Beams. In *Book Cold Chemistry: Molecular Scattering and Reactivity Near Absolute Zero*; Osterwalder, A., Dulieu, O., Eds.; The Royal Society of Chemistry: London, UK, 2018; pp. 92–149.
32. Hutson, J.M.; Green, S. *MOLSCAT: MOLecular SCATtering*, Version 14; Astrophysics Source Code Library: Houghton, MI, USA, 2012; Available online: https://ui.adsabs.harvard.edu/abs/2012ascl.soft06004H/abstract (accessed on 20 September 2022).

33. Harvey, A.H.; Lemmon, E.W. Correlation for the Vapor Pressure of Heavy Water from the Triple Point to the Critical Point. *J. Phys. Chem. Ref. Data* **2002**, *31*, 173–181. [CrossRef]

34. Naulin, C.; Bergeat, A. Low-energy Scattering in Crossed Molecular Beams. In *Book Cold Chemistry: Molecular Scattering and Reactivity Near Absolute Zero*; Osterwalder, A., Dulieu, O., Eds.; The Royal Society of Chemistry: London, UK, 2018; pp. 92–149.