## The water line emission and ortho-to-para ratio in the Orion Bar photon-dominated region

T. Putaud 1 , X. Michaut 1 , F. Le Petit 1 , E. Roue GLYPH&lt;11&gt; 1 , and D. C. Lis 2 1 ;

- 1 Sorbonne Universit, Observatoire de Paris, Universit PSL, CNRS, LERMA, F-75014, Paris, France
- 2 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, CA 91109, USA

c 2019. All rights reserved. Received date GLYPH&lt;13&gt; / Accepted date

## ABSTRACT

Context. The ortho-to-para ratio (OPR) of water in the interstellar medium (ISM) is often assumed to be related to the formation temperature of water molecules, making it a potentially interesting tracer of the thermal history of interstellar gas.

Aims. Avery low OPR of 0.1-0.5 was previously reported in the Orion Bar photon-dominated region (PDR), based on observations of two optically thin H 18 2 Olines which were analyzed by using a single-slab large velocity gradient (LVG) model. The corresponding spin temperature does not coincide with the kinetic temperature of the molecular gas in this UV-illuminated region. This was interpreted as an indication of water molecules being formed on cold icy grains which were subsequently released by UV photodesorption.

Methods. A more complete set of water observations in the Orion Bar, including seven H 16 2 O lines and one H 18 2 O line, carried out using Herschel / HIFI instrument, was reanalyzed using the Meudon PDR code to derive gas-phase water abundance and the OPR. The model takes into account the steep density and temperature gradients present in the region.

Results. The model line intensities are in good agreement with the observations assuming that water molecules formed with an OPR corresponding to thermal equilibrium conditions at the local kinetic temperature of the gas and when solely considering gas-phase chemistry and water gas-grain exchanges through adsorption and desorption. Gas-phase water is predicted to arise from a region deep into the cloud, corresponding to a visual extinction of A V GLYPH&lt;24&gt; 9, with a H 16 2 O fractional abundance of GLYPH&lt;24&gt; 2 GLYPH&lt;2&gt; 10 = 7 and column density of (1.4 GLYPH&lt;6&gt; 0.8) GLYPH&lt;2&gt; 10 15 cm 2 = for a total cloud depth of A V = 15. A line-of-sight average ortho-to-para ratio of 2.8 GLYPH&lt;6&gt; 0.2 is derived.

Conclusions. The observational data are consistent with a nuclear spin isomer repartition corresponding to the thermal equilibrium at a temperature of (36 GLYPH&lt;6&gt; 2) K, much higher than the spin temperature previously reported for this region and close to the gas kinetic temperature in the water-emitting gas.

Key words. ISM: molecules - ISM: individual objects: Orion Bar - (ISM:) photon-dominated region (PDR) - ISM: lines and bands - Submillimeter: ISM

## 1. Introduction

Photon-dominated regions (PDRs) are surface layers of molecular clouds irradiated by a strong UV photon flux, such as those present in star-forming regions. Studies of the physical and chemical properties of molecules in such regions are of great interest for understanding the feedback of young stars on their surrounding medium.

The Orion Molecular Cloud 1 (OMC-1) contains two embedded star-forming regions: Orion BN KL and Orion South, as / well as a group of young massive stars known as the Trapezium Cluster. The Trapezium, and in particular the brightest O6 type star Θ 1 Ori C, has irradiated the surrounding molecular cloud creating an HII region, which is bordered by the Orion Bar on its southeastern side and by the Orion Ridge on its western side (Rodriguez-Franco et al. 1998). Close to Earth (414 pc; Menten et al. 2007) and with a nearly edge-on geometry (Hogerheijde et al. 1995; Jansen et al. 1995) explained by the blister model (see for example Wen &amp; O'Dell 1995 and references therein), the Orion Bar is an excellent place to test PDR models (Tielens &amp; Hollenbach 1985; Sternberg &amp; Dalgarno 1989; Le Petit et al. 2006; Hollenbach et al. 2009; AndreeLabsch et al. 2017).

The Orion Bar has been observed extensively in multiple gas tracers (Fuente et al. 1996; Larsson et al. 2003; Leurini et al. 2006; van der Tak et al. 2012, 2013; Cuadrado et al.

2015, 2017; Nagy et al. 2017; Joblin et al. 2018). Previous studies have led to a UV flux of G 0 = 1 GLYPH&lt;0&gt; 4 GLYPH&lt;2&gt; 10 4 in Habing units (Tielens &amp; Hollenbach 1985; Marconi et al. 1998) impinging upon a molecular cloud with a mean H2 density of about 5 GLYPH&lt;2&gt; 10 4 cm = 3 (Wyrowski et al. 1997) and a kinetic temperature of 85 K (Hogerheijde et al. 1995). A model of the Orion Bar with a clumpy medium, based on high-density clumps ( n H = 10 6 GLYPH&lt;0&gt; 10 7 cm GLYPH&lt;0&gt; 3 ) embedded in an interclump medium ( n H = 5 GLYPH&lt;2&gt; 10 4 GLYPH&lt;0&gt; 10 5 cm GLYPH&lt;0&gt; 3 ), was proposed to explain the following: the excitation of atomic lines, some excited CO lines, and warm H2 observations (Parmar et al. 1991; Tauber et al. 1994; van der Werf et al. 1996). The presence of high-density clumps was also suggested by observations of H2CO, HCO + , and HCN (Hogerheijde et al. 1995; Young Owl et al. 2000; Lis &amp; Schilke 2003). However, the clumpy model is challenged by the very high resolution observations provided by the Atacama Large Millimeter submillimeter / Array (ALMA) radiotelescope (Goicoechea et al. 2016, 2017), which show that, to the first order, the structure of the PDR is a compressed layer at high-pressure where a warm chemistry takes place, leading to the presence of molecules such as SH + and high-J excited CO (Joblin et al. 2018). This structure is also observed in other PDRs, such as Trumpler 14 in the Carina nebula (Wu et al. 2018). High-density structures exist inside the PDR, as seen in OH (Goicoechea et al. 2011; Parikka et al. 2017) and H 13 CN maps (Lis &amp; Schilke 2003), as well as in ALMA maps (Goicoechea et al. 2016). The high-pressure layer can be

T. Putaud et al.: The water line emission and OPR in the Orion Bar PDR

Table 1: Observations used for the analysis of water line emission in the Orion Bar.

| Molecule   | S a   | Transition         |   GLYPH<23> b (GHz) |   E up b (K) |   g up c |   HPBW d (') |   GLYPH<17> mb d | HIFI band   | Observing mode                             | t int e (sec)    | Observation Id.                             |
|------------|-------|--------------------|---------------------|--------------|----------|--------------|------------------|-------------|--------------------------------------------|------------------|---------------------------------------------|
| H 16 2 O   | o     | 1 10 GLYPH<0> 1 01 |             556.936 |        60.96 |        9 |           37 |             0.62 | 1B 1A 1B    | Load Chop Load Chop On-The-Fly             | 103.3 126.3 4.0  | 1342215923 f 1342218527 g 1342215922        |
| H 16 2 O   | p     | 1 11 GLYPH<0> 0 00 |            1113.34  |        53.43 |        3 |           19 |             0.63 | 4B 4B       | Load Chop On-The-Fly                       | 105.8 2.9        | 1342217720 1342215970                       |
| H 16 2 O   | o     | 2 12 GLYPH<0> 1 01 |            1669.9   |       114.38 |       15 |           12 |             0.58 | 6B 6B 6B    | Fast DBS Fast DBS On-The-Fly               | 1274 134.6 4.0   | 1342229840 h 1342218426 i 1342251054        |
| H 16 2 O   | p     | 2 02 GLYPH<0> 1 11 |             987.927 |       100.85 |        5 |           21 |             0.64 | 4A 4A       | Load Chop On-The-Fly                       | 105.8 2.9        | 1342218628 1342218217                       |
| H 16 2 O   | p     | 2 11 GLYPH<0> 2 02 |             752.033 |       136.94 |        5 |           28 |             0.65 | 2B 2B 2B 2B | Load Chop On-The-Fly On-The-Fly On-The-Fly | 94.0 3.9 4.9 2.9 | 1342216378 1342190848 1342190849 1342203221 |
|            | o     | 2 21 GLYPH<0> 2 12 |            1661.01  |       194.09 |       15 |           13 |             0.58 | 2B 6B       | On-The-Fly Fast DBS                        | 3.9 1274         | 1342203222 1342229840                       |
|            | o     | 3 12 GLYPH<0> 3 03 |            1097.37  |       249.44 |       21 |           19 |             0.63 | 4B          | Load Chop                                  | 105.8            | 1342217720                                  |
| H 18 2 O   | o     | 1 10 GLYPH<0> 1 01 |             547.676 |        60.46 |        9 |           38 |             0.62 | 1A 1A       | Freq. Switch On-The-Fly                    | 5278 3.96        | 1342205273 1342216355                       |
|            | p     | 1 11 GLYPH<0> 0 00 |            1101.7   |        52.87 |        3 |           19 |             0.63 | 4B          | Load Chop                                  | 105.8            | 1342217720 j                                |

Notes. ( a ) Nuclear spin state, o stands for ortho and p for para. ( b ) Frequency of the transition (up-down) and rotational energy of the upper level in relation to the para 0 00 level from the HITRAN2016 database (Gordon et al. 2017) ( c ) Statistical weight of the upper level. ( d ) Half-power beam width and main-beam e GLYPH&lt;14&gt; ciency of the Herschel telescope, averaged between H and V polarizations, at the water line frequency (Mueller et al. 2014) . ( e ) Integration time from FITS header. ( f ) Hereafter, this line is labeled H 16 2 O 110 GLYPH&lt;0&gt; 101-a. ( g ) Hereafter, this line is labeled H 16 2 O 110 GLYPH&lt;0&gt; 101-b. ( h ) Hereafter, this line is labeled H 16 2 O212 GLYPH&lt;0&gt; 101-a. ( ) i Hereafter, this line is labeled H 16 2 O212 GLYPH&lt;0&gt; 101-b. ( j ) H 18 2 O111 GLYPH&lt;0&gt; 000 para line is not unambiguously detected.

interpreted by the photo-evaporation e GLYPH&lt;11&gt; ect as explained by Bron et al. (2018) with the PDR hydrodynamical code Hydra.

An important tracer for studying the physical and chemical evolution of the interstellar medium (ISM) is H2O. This molecule (like other hydrides) has two identical protons and exists in two nuclear configurations, called nuclear spin isomers, where the hydrogen nuclear spins are either parallel (ortho) or antiparallel (para). At the thermal equilibrium, the ortho-to-para ratio (OPR) of gas-phase water molecules is determined by the kinetic temperature of the gas. Some non-thermal equilibrium OPRs have been reported in di GLYPH&lt;11&gt; erent regions (Hogerheijde et al. 2011; Lis et al. 2013; Choi et al. 2014). In particular, Choi et al. (2014) derived a very low OPR of 0.1-0.5 in the Orion Bar, based on observations of two H 18 2 O emission lines. With such a low OPR, the so-called spin temperature does not exceed 10 K, which is much lower than the gas kinetic temperature in this region.

sophisticated models, has sometimes modified the reported OPR of hydride molecules in the ISM. Better calibrated observational data and improved computational models have indeed shown that the low OPR value measured in TW Hydrae protoplanetary disk (Hogerheijde et al. 2011) is very model-dependent (Salinas et al. 2016). Furthermore, the observations are consistent, within the observational uncertainties, with the high-temperature equilibrium value.

To explain the discrepancy between kinetic and spin temperatures, it was suggested that the spin temperature might reflect the temperature of the grains at the time the molecules formed. Indeed, a PDR model which considers oxygen grain chemistry has shown that the photodesorption of water molecules formed on grains could be an e GLYPH&lt;14&gt; cient pathway for gas-phase water production (Hollenbach et al. 2009). After desorption from grains, the non-thermal equilibrium OPR could have been preserved due to the very low pressure of the ISM (Cacciani et al. 2012). However, a recent laboratory study of water molecules photodesorbed from cold surfaces using a UV laser produced a gas-phase water with an OPR in agreement with the high thermal equilibrium value of three (Hama et al. 2016).

Further studies of the impact of grains on the abundance and quantum state of gas-phase molecules have yet to be carried out. However, the reanalysis of the observational data, using more

We present a new analysis of the H 16 2 O and H 18 2 O lines observed toward the Orion Bar with Herschel / HIFI, which are reduced using the latest pipeline and with the latest calibration (Mueller et al. 2014). Careful corrections are applied to the spectra to account for the spatial o GLYPH&lt;11&gt; sets and beam coupling effects. The Meudon PDR code (Le Petit et al. 2006) is used to constrain the physical conditions inside the Orion Bar and derive the gas-phase water OPR. The impact of grain processes, such as water adsorption and thermal or photo-induced desorption, is also investigated.

## 2. Observations

Observations of the Orion Bar discussed here were carried out at the CO + peak position (Stoerzer et al. 1995) ( GLYPH&lt;11&gt; 2000 = 05 h 35 m 20 61 s : , GLYPH&lt;14&gt; 2000 = GLYPH&lt;0&gt; 05 25 14 0 GLYPH&lt;14&gt; 0 : 00 ) as part of the HEXOS 1 (Bergin et al. 2010) guaranteed time key program (GTKP). The Heterodyne Instrument for the FarInfrared (HIFI, de Graauw et al. 2010) aboard the Herschel Space Observatory (Pilbratt et al. 2010), operating with the wide-band spectrometer (WBS) as a backend, provided spectra at 1.1 MHz resolution for H and V orthogonal polarizations. The data were reduced with the Herschel Interactive

1 Herschel observations of the EXtra-Ordinary Sources

Processing Environement (HIPE, Ott 2010) version 14.1.0. Within the whole HIFI's frequency range (480-1250 GHz and 1410-1910 GHz), seven rotational emission lines of H 16 2 O and one of H 18 2 O were detected from March 2011 to October 2011, either in load chop, fast dual beam switch (DBS), or frequency switch observing mode. The upper limit of the ground state para H 18 2 O line was also included in the analysis. Observations used in this article are listed in Table 1. Most of these lines were already presented in an extensive spectral survey carried out by Nagy et al. (2017).

In addition to the spectral scan observations, nine on-the-fly (OTF) HIFI maps using the WBS as a backend were also analyzed (see Table 1). Within the HEXOS GTKP (Bergin et al. 2010), Habart (2011) proposal, and calibration programs , ob2 servations were obtained from February 2010 to September 2012. The total field of view is roughly 50 00 GLYPH&lt;2&gt; 100 00 (except for H 16 2 O 110 GLYPH&lt;0&gt; 101 and 212 GLYPH&lt;0&gt; 101 where only a 100 00 perpendicular strip through the Bar was performed) with a position angle perpendicular to the Bar. The maps were regridded using the GILDAS 3 software leading to a spatial resolution typically 5 % lower than the Herschel telescope half-power beam width (HPBW). All water spectra and water maps are shown in the Appendix A.

Further calibration corrections must be applied to the data. First, the antenna temperature on the T A GLYPH&lt;3&gt; scale, given by the pipeline, should be corrected by the ratio of the forward e GLYPH&lt;14&gt; -ciency ( GLYPH&lt;17&gt; l = 0 96, Roelfsema et al. 2012) and the main-beam : e GLYPH&lt;14&gt; ciency ( GLYPH&lt;17&gt; mb ) to produce spectra on the main-beam brightness temperature scale Tmb (Mueller et al. 2014).

$$T _ { m b } & = \frac { \eta _ { l } } { \eta _ { m b } } T _ { A } ^ { * } & ( 1 ) & \frac { \tilde { \text{ion} } \ B \cdot } { \text{correct} \varepsilon \cdot }$$

This correction only considers the main detection lobe of the telescope. The latest values of GLYPH&lt;17&gt; mb and the HPBW, derived from Mars observations and complete optical model of the Herschel telescope (Mueller et al. 2014; Shipman et al. 2017), were used in the analysis. These numbers di GLYPH&lt;11&gt; er up to 20 % from the previous estimates assuming a simplified Gaussian beam shape (Roelfsema et al. 2012). The exact values used for each H and V observations are listed in Table A.1 in the Appendix A.

Furthermore, owing to the variation in the HPBW with frequency, the beam coupling correction has to be carefully considered, as discussed below. To quantify this e GLYPH&lt;11&gt; ect, three high angular resolution CO maps were convolved to the beam widths corresponding to the various water transitions. The beam coupling factor Ω is defined as the ratio between the convolved intensity at the targeted coordinates and the maximum intensity toward the Bar at high resolution. The CO 1 GLYPH&lt;0&gt; 0 map is part of the CARMANROOrion Survey (Kong et al. 2018) and was obtained by combining the CARMA 4 interferometric data with NRO45 5 singledish observations leading to a resolution of 9 00 over a 2 GLYPH&lt;14&gt; GLYPH&lt;2&gt; 2 GLYPH&lt;14&gt; field of view. The 12 CO 6 GLYPH&lt;0&gt; 5 and 13 CO 6 GLYPH&lt;0&gt; 5 maps were obtained using the Caltech Submillimeter Observatory (CSO) (Lis &amp;Schilke 2003) at 11 00 resolution and cover respectively a 3 0 GLYPH&lt;2&gt; 4 0 and 3 0 GLYPH&lt;2&gt; 2.5 0 regions.

Fig. 1: Averaged spectra of H 16 2 O and H 18 O lines 2 toward the Orion Bar on the main-beam brightness temperature scale, not corrected for the beam coupling factor. The gray areas represent the velocity integration range. The blue dash-dotted line and red dashed line are the expected velocity components from Orion Ridge (9.0 km s = 1 ) and Orion Bar (10.6 km s = 1 ). H 18 2 Olines were multiplied by a factor of 5.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-723b59f1-75b8-463d-8582-73d4ae366053.jpg)

## 3. Results

Integrated line intensities of the water spectra shown in Fig. 1 are given in Table 2. On average, the integration was performed from a local standard of rest velocity ( V LSR) of 7.7 km s = 1 to 13.5 km s = 1 . Precise integration ranges for each line are listed in Table A.1. Linear baseline subtraction was carefully applied except for frequency switch observation mode where a fourth order polynomial baseline was used because of residual standing waves. The uncertainties of the integrated line intensities, in units of K km s = 1 , are computed as

$$\underset { \substack { \text{souping} \\ \text{and} \\ \text{ARMA} - \\ \dots \dots } } { \text{souping} }$$

where the summation runs over all the channels N outside the line integration velocity range, Nint is the number of channels

2 calibration pvhifi 37 and calibration pvhifi 85

3 https: // www.iram.fr IRAMFR GILDAS / / /

- 4 Combined Array for the Research in Millimeter Astronomy

5 Nobeyama Radio Observatory 45 m telescope

inside the integration velocity range and ∆ V is the width of a velocity channel.

The average intensities of the H and V polarizations agree with the values previously reported (Nagy et al. 2017; Choi 2015), but are a factor of two higher than the intensities derived from SPIRE observations averaged along the Bar (Habart et al. 2010). The ortho H 18 2 Oline intensity is also very close to the one derived by Choi et al. (2014). However the para H 18 2 O line (a 4.7 GLYPH&lt;27&gt; detection in Choi et al. 2014) is not clearly detected in our data set given the noise level in the spectrum.

## 3.1. Line profiles

Figure 1 shows the line profiles of the eight detected lines, and the spectrum in the para H 18 2 O 111 GLYPH&lt;0&gt; 000 line region. Some of the lines have a profile inconsistent with a single Gaussian, in particular the H 16 2 O ground state lines (111 GLYPH&lt;0&gt; 000, 110 GLYPH&lt;0&gt; 101, and 212 GLYPH&lt;0&gt; 101) exhibit two velocity components around 9.5 km s = 1 and 11.5 km s = 1 . Previous studies of tracers arising close to the CO + peak reported velocity components from the Orion Bar at 10.6 km s = 1 , as well as features from the Orion Ridge at 9 km s = 1 (van der Tak et al. 2013; Goicoechea et al. 2015; Nagy et al. 2017; Cuadrado et al. 2017). Moreover lines between VLSR = 8 km s GLYPH&lt;0&gt; 1 and VLSR = 9 km s GLYPH&lt;0&gt; 1 were attributed to OMC1 as suggested by IRAM 30 m telescope observations of small hydrocarbons (Cuadrado et al. 2015). Strips across the Bar have shown that CH + peak velocity shifts from 9 km s = 1 in front of the Bar to 11 km s = 1 behind (Parikka et al. 2017).

In the case of water emission, the two-component line profiles are more likely explained by self-absorption in the lowenergy lines (Choi 2015), which are expected to be optically thicker than the excited lines. Such an explanation is supported by the fact that excited lines exhibit a single component at an intermediate velocity, close to the expected velocity of the Bar (Nagy et al. 2017). The self-absorption dip also matches this velocity. Furthermore, a full width at half maximum (FWHM) of GLYPH&lt;24&gt; 4.5 km s = 1 is measured for the lines showing a self-absorption dip in good agreement with optically thick lines in the Orion Bar (Hogerheijde et al. 1995; Nagy et al. 2017) and GLYPH&lt;24&gt; 2.5 km s = 1 for the other ones which is close to the typical line width observed in this region (Choi et al. 2014; Cuadrado et al. 2015; Nagy et al. 2017). The FWHM line widths are reported in Table A.1.

Nevertheless, CO lines (Joblin et al. 2018) show a weak asymmetry at low velocities that could represent signal coming from the Orion Ridge or OMC-1. Assuming a peak velocity of 10.6 km s = 1 and considering a symmetric line profile, velocity integration range is derived for each line from the half width measured at high velocity (see Table A.1). Even if a possible contamination of the line by surrounding features could still lead to an overestimate of the intensity, the resulting intensities are our best-estimates, because the optically thick character of the lines makes a fitting method unreliable.

## 3.2. Beam dilution

The observed water lines cover a wide range of frequencies, leading to a significant variation in the Herschel telescope beam width. Extreme beam sizes at the spectral scan coordinates are shown as white circles in the H 16 2 O202 GLYPH&lt;0&gt; 111 map in Fig. 2. In this map, corresponding to an intermediate frequency, the maximum HIFI beam size is larger than the spatial width of the Bar. Thus the coupling of the beam to the water-emitting region has to be

Fig. 2: H 16 2 O202 GLYPH&lt;0&gt; 111 emission map at 988 GHz integrated from 7.7 km s = 1 to 13.5 km s = 1 . Contours represent the emission of 13 CO 6 GLYPH&lt;0&gt; 5 at 95%, 90%, 75% and 50% of the maximum intensity, integrated from 2 km s = 1 to 17 km s = 1 , observed with the CSO at a resolution of 11 00 . Black star marks the CO + peak position, white circles show the HIFI HPBW at the extreme frequencies, black and magenta disks show the CSO and 988 GHz HIFI beam width (see Table 1). FUV field from the Trapezium cluster comes from the upper right corner.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-f743e8cb-7149-461e-9a8e-098c83550daa.jpg)

considered when comparing the observed line intensities at different frequencies with model predictions. In addition, Figure 2 clearly shows that water emission peaks deeper in the molecular cloud with respect to the CO + peak position. Thus the observations were not centered on the water emission peak in the Orion Bar and only a fraction of the emission is included in the beam.

The o GLYPH&lt;11&gt; set e GLYPH&lt;11&gt; ect can be easily taken into account using water maps by comparing the intensity at the targeted coordinates with the mean intensity along the Bar. This gives a rough estimate of the beam coupling factor of about 0.8. Nevertheless, applying this method does not consider the fact that the maximum intensity is already impacted by the size of the beam, in particular for low frequency transitions for which the HPBW of the Herschel telescope is the largest.

To retrieve the maximal intensity of each line in the Bar, the e GLYPH&lt;11&gt; ect of the pointing o GLYPH&lt;11&gt; set and the beam width variations is investigated from high resolution maps. The intensity along the Bar at the full resolution is derived and compared to the intensity at the observed coordinates after convolution at the HPBW corresponding to the line frequency. Three CO maps were used to derive the beam coupling factor. These maps have a higher resolution than the most resolved water maps (9 00 for 12 CO1 GLYPH&lt;0&gt; 0, Kong et al. 2018 and 11 00 for 12 CO 6 GLYPH&lt;0&gt; 5 and 13 CO 6 GLYPH&lt;0&gt; 5, Lis &amp; Schilke 2003) and are large enough to be convolved without edge e GLYPH&lt;11&gt; ect. Moreover, they span di GLYPH&lt;11&gt; erent excitation conditions and di GLYPH&lt;11&gt; erent spatial extensions. Indeed, 12 CO emission is slightly wider than the water one, whereas 13 CO peaks slightly deeper into the cloud (4 00 , see H 16 2 O111 GLYPH&lt;0&gt; 000 and 212 GLYPH&lt;0&gt; 101 maps in Fig. A.5). So averaging the beam coupling factors derived from each CO maps should give a good estimate of its value and its uncertainty (Table. 2). We note that the 4 00 o GLYPH&lt;11&gt; set between 13 CO

Table 2: Observed intensities of water lines toward the Orion Bar and their correction for the beam coupling factor.

| Molecule   | Transition         | I Obs a              | I Corr               | Ω b                | Ω b                  | Ω b                  | Ω b                  | 1 = Ω c            |
|------------|--------------------|----------------------|----------------------|--------------------|----------------------|----------------------|----------------------|--------------------|
|            |                    | (Kkms = 1 )          | (Kkms = 1 )          | H 2 O              | 12 CO 1 GLYPH<0> 0 d | 12 CO 6 GLYPH<0> 5 e | 13 CO 6 GLYPH<0> 5 e |                    |
| H 16 2 O   | 1 10 GLYPH<0> 1 01 | 21.13 GLYPH<6> 0.95  | 27.48 GLYPH<6> 1.19  | 0.83 GLYPH<6> 0.07 | 0.80 GLYPH<6> 0.02   | 0.85 GLYPH<6> 0.02   | 0.66 GLYPH<6> 0.04   | 1.32 GLYPH<6> 0.17 |
| H 16 2 O   | 1 11 GLYPH<0> 0 00 | 10.28 GLYPH<6> 0.75  | 12.57 GLYPH<6> 0.89  | 0.61 GLYPH<6> 0.04 | 0.86 GLYPH<6> 0.03   | 0.85 GLYPH<6> 0.03   | 0.73 GLYPH<6> 0.08   | 1.24 GLYPH<6> 0.13 |
| H 16 2 O   | 2 12 GLYPH<0> 1 01 | 8.91 GLYPH<6> 0.60   | 12.46 GLYPH<6> 0.83  | 0.78 GLYPH<6> 0.04 | 0.86 GLYPH<6> 0.03   | 0.79 GLYPH<6> 0.01   | 0.72 GLYPH<6> 0.04   | 1.28 GLYPH<6> 0.11 |
| H 16 2 O   | 2 02 GLYPH<0> 1 11 | 8.06 GLYPH<6> 0.59   | 9.41 GLYPH<6> 0.66   | 0.3 GLYPH<6> 0.2   | 0.85 GLYPH<6> 0.04   | 0.86 GLYPH<6> 0.03   | 0.73 GLYPH<6> 0.07   | 1.27 GLYPH<6> 0.13 |
| H 16 2 O   | 2 11 GLYPH<0> 2 02 | 7.82 GLYPH<6> 0.60   | 9.95 GLYPH<6> 0.75   | 0.77 GLYPH<6> 0.01 | 0.82 GLYPH<6> 0.03   | 0.87 GLYPH<6> 0.03   | 0.69 GLYPH<6> 0.05   | 1.24 GLYPH<6> 0.16 |
| H 16 2 O   | 2 21 GLYPH<0> 2 12 | 2.62 GLYPH<6> 0.31   | 3.13 GLYPH<6> 0.34   |                    | 0.87 GLYPH<6> 0.03   | 0.79 GLYPH<6> 0.01   | 0.74 GLYPH<6> 0.04   | 1.26 GLYPH<6> 0.10 |
| H 16 2 O   | 3 12 GLYPH<0> 3 03 | 1.85 GLYPH<6> 0.18   | 2.24 GLYPH<6> 0.20   |                    | 0.86 GLYPH<6> 0.03   | 0.85 GLYPH<6> 0.03   | 0.73 GLYPH<6> 0.08   | 1.24 GLYPH<6> 0.13 |
| H 18 2 O   | 1 10 GLYPH<0> 1 01 | 0.29 GLYPH<6> 0.02   | 0.36 GLYPH<6> 0.02   | 0.6 GLYPH<6> 0.4   | 0.79 GLYPH<6> 0.03   | 0.85 GLYPH<6> 0.02   | 0.65 GLYPH<6> 0.05   | 1.33 GLYPH<6> 0.18 |
| H 18 2 O   | 1 11 GLYPH<0> 0 00 | < 0.34 GLYPH<6> 0.15 | < 0.42 GLYPH<6> 0.16 |                    | 0.86 GLYPH<6> 0.03   | 0.85 GLYPH<6> 0.03   | 0.73 GLYPH<6> 0.08   | 1.24 GLYPH<6> 0.13 |

Notes. ( a ) Integration range is around 7.7 km s = 1 and 13.5 km s = 1 . Precise integration range for each line is given in Table A.1. ( b ) Hand V averaged beam coupling factors from each tracer. ( c ) Mean global corrective coe GLYPH&lt;14&gt; cients from CO maps. ( d ) Kong et al. (2018). ( e ) Lis &amp; Schilke (2003). Values for 13 CO are obtained with a 4 00 shift toward the ionization front.

and water emission corresponds to a typical pointing error of the CSO telescope at this high frequency. Since only the spatial extension of the emission is considered in the estimation of the beam coupling factor, the 13 CO map was shifted toward the ionization front by 4 00 to match the water emission, in order not to add artificial bias due to the position.

The incoming UV radiation field is based on the interstellar radiation field (ISRF) (Mathis et al. 1983), scaled by a multiplicative factor G 0 . Mathis UV radiation field is about 1.3 factor as high as Habing one (Habing 1968) in the spectral range of 91.2 to 111 nm (Wu et al. 2018).

The beam coupling factors adopted for each observed transitions are given in Table A.1 together with our best-estimates for the intensity of each lines. The beam coupling factors are close to 0.8 whereas the values retrieved for high-J CO lines assuming a 2 00 filament for the Bar range from 0.07 to 0.23 (Joblin et al. 2018). Each H and V spectrum is corrected separately for beam coupling to account for a pointing o GLYPH&lt;11&gt; set between the two polarizations and the resulting spectra are then averaged with root mean square (rms) weighting:

$$T _ { \text{Weight} } & = \frac { w _ { H } T _ { H } + w _ { V } T _ { V } } { w _ { H } + w _ { V } } & ( 3 ) & \text{theory,} \frac { \text{tensity,} } { \text{choric on} }$$

where TH V = is the main-beam brightness temperature for H and V polarizations, corrected for the beam coupling factor, and wH V = the weight for each polarizations determined by 1 =GLYPH&lt;27&gt; 2 , with GLYPH&lt;27&gt; the rms noise measured outside the line window. The averaged spectra without beam coupling corrections are shown in Fig. 1 and the line intensities are listed in Table 2. The final uncertainties for the corrected intensities are computed as the quadratic sum of the spectral rms related to the signal-to-noise ratio, the absolute intensity calibration error (Teyssier et al. 2017), the main-beam coe GLYPH&lt;14&gt; cient uncertainty (Mueller et al. 2014), and the beam coupling factor uncertainty.

## 4. Analysis

To model the water vapor emission in the Orion Bar, we use the Meudon PDR code 6 (Le Petit et al. 2006) to fit the observed line intensities. The Meudon PDR code simulates a stationary planeparallel slab of gas and dust and computes the radiative transfer along a line-of-siglht. At each position within the cloud, thermal and chemical balances are computed, as well as non-local thermodynamic equilibrium (non-LTE) level populations. Level populations and resulting line intensities are computed considering collisional and radiative processes, as well as chemical formation and destruction in various rotational levels. For radiative processes, the code takes into account non-local pumping by the continuum (background and dust emission) and line emission, as explained in Gonzalez Garcia et al. (2008).

6 https: // ism.obspm.fr

In the currently released version of the code (1.5.2), the formation of water molecules only includes gas-phase chemistry. The nuclear spin distribution of the water molecules formed follows the kinetic temperature at each point within the cloud with a non-thermal equilibrium rotational excitation.

## 4.1. Gas-phase chemistry model

Investigation of the input parameters of the Meudon PDR code such as the inclination angle, cloud depth, UV radiation field intensity, and density or thermal pressure (considering either isochoric or isobaric models) is performed by using the Interstellar Medium Database 7 (ISMDB) to compare the observed line intensities with grids of model predictions (precomputed with the Meudon PDR code 1.5.2 and available online). Intensities given by the ISMDB (in cgs units) are compared with observations using the following conversion:

$$\inf _ { \substack { \text{Fig. 1} \\ \text{certain-} \\ \text{qualratic} } } \ I = \frac { 2 k } { \lambda ^ { 3 } } \int T _ { m b } d V$$

with k the Boltzmann constant (erg K = 1 ) and GLYPH&lt;21&gt; the wavelength (cm). Intensities in cgs units are given in Table 4.

As already shown using more than twenty lines emitted by nine species, the line emission at the edge of the Orion Bar is better explained by an isobaric model than by a constant density model (Joblin et al. 2018). This is understandable because a steep density and temperature gradient is expected to exist at the edge of a UV-illuminated molecular cloud (see Figs. 4 and 5). This agrees with other studies of the Orion Bar (Marconi et al. 1998; Allers et al. 2005), as well as in the northwest PDR in NGC 7023 (Joblin et al. 2018) and Carina PDR (Wu et al. 2018). A theoretical explanation, based on the photo-evaporation process, was provided by Bron et al. (2018). Consequently, an isobaric model is adopted.

To reproduce the nearly edge-on geometry of the Bar (Wen &amp; O'Dell 1995; Hogerheijde et al. 1995; Jansen et al. 1995; Walmsley et al. 2000), the maximum viewing angle of 60 GLYPH&lt;14&gt; accepted by the code is adopted. This angle is defined as the angle between the normal to the ionization front and the line-ofsight, meaning that a viewing angle of 90 GLYPH&lt;14&gt; corresponds to an

7 http: // ismdb.obspm.fr

Table 3: Best-fit input parameters of the Meudon PDR code for the corrected intensities.

| Parameters                                                                      | ISMDB model a              | Ads / Des model b          |
|---------------------------------------------------------------------------------|----------------------------|----------------------------|
| GLYPH<18> ( GLYPH<14> )                                                         | 60                         | 60 20                      |
| A tot V                                                                         | 20                         |                            |
| G 0 (Mathis Unit)                                                               | 5 GLYPH<2> 10 4            | 3.1 GLYPH<2> 10 4          |
| P th (Kcm = 3 )                                                                 | 1 GLYPH<2> 10 8            | 2.8 GLYPH<2> 10 8          |
| Cosmic rays (s = 1 per H 2 ) Rv                                                 | 1 GLYPH<2> 10 = 16 3.1     | 5 GLYPH<2> 10 = 16 5.62    |
| N H / E(B-V) (cm = 2 )                                                          | 5.8 GLYPH<2> 10 21         | 1.05 GLYPH<2> 10 22        |
| Mass grains / Mass gas                                                          | 0.01                       | 0.01                       |
| Grain size distribution                                                         | / GLYPH<11> GLYPH<0> 3 : 5 | / GLYPH<11> GLYPH<0> 3 : 5 |
| min grain radius (cm)                                                           | 1 GLYPH<2> 10 = 7          | 3 GLYPH<2> 10 = 7          |
| max grain radius (cm)                                                           | 3 GLYPH<2> 10 = 5          | 3 GLYPH<2> 10 = 5          |
| H 2 O binding energy (K) H 2 O photodesorption yield c H 2 internal energy used | No                         | 5600 2 GLYPH<2> 10 = 3 Yes |

Notes. ( a ) Pure gas-phase chemistry model with standard ISMDB inputs. ( b ) Adsorption and desorption optimized model with Joblin et al. (2018) inputs and H2 internal energy. No additional scaling factor was applied to correct for a more edge-on geometry. ( c ) The water photodesorption yield is in units of molecules per incident photon.

Table 4: Intensities of water lines toward the Orion Bar corrected for the beam coupling factor, in erg cm = 2 s = 1 sr = 1 , and preticted intensities from the Meudon PDR code using input parameters in Table 3 for the pure gas-phase chemsitry ISMDB model and the adsorption and desorption optimized model with H2 internal energy.

| Mol.     | Trans.                                                                                                                               | Corrected intensities a                                                                                                                                              | ISMDB model                                                                                                                                     | Ads / Des model                                                                                                                                 |
|----------|--------------------------------------------------------------------------------------------------------------------------------------|----------------------------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------|-------------------------------------------------------------------------------------------------------------------------------------------------|
| H 16 2 O | 1 10 GLYPH<0> 1 01 1 11 GLYPH<0> 0 00 2 12 GLYPH<0> 1 01 2 02 GLYPH<0> 1 11 2 11 GLYPH<0> 2 02 2 21 GLYPH<0> 2 12 3 12 GLYPH<0> 3 03 | 4.9(2) GLYPH<1> 10 = 6 1.8(1) GLYPH<1> 10 = 5 6.0(4) GLYPH<1> 10 = 5 9.3(7) GLYPH<1> 10 = 6 4.3(3) GLYPH<1> 10 = 6 1.5(2) GLYPH<1> 10 = 5 3.0(3) GLYPH<1> 10 = 6 = 8 | 4.9 GLYPH<1> 10 = 6 1.8 GLYPH<1> 10 = 5 3.4 GLYPH<1> 10 = 5 1.2 GLYPH<1> 10 = 5 5.9 GLYPH<1> 10 = 6 1.6 GLYPH<1> 10 = 5 2.9 GLYPH<1> 10 = 6 = 7 | 3.5 GLYPH<1> 10 = 6 1.5 GLYPH<1> 10 = 5 3.1 GLYPH<1> 10 = 5 1.4 GLYPH<1> 10 = 5 3.2 GLYPH<1> 10 = 6 1.3 GLYPH<1> 10 = 5 1.9 GLYPH<1> 10 = 6 = 7 |
| H 18 2 O | 1 10 GLYPH<0> 1 01                                                                                                                   | 6.1(3) GLYPH<1> 10                                                                                                                                                   | 1.4 GLYPH<1> 10                                                                                                                                 | 1.0 GLYPH<1> 10                                                                                                                                 |
| H 18 2 O | 1 11 GLYPH<0> 0 00                                                                                                                   | < 5.8(22) GLYPH<1> 10 = 7                                                                                                                                            | 6.6 GLYPH<1> 10 = 7                                                                                                                             | 3.6 GLYPH<1> 10 = 7                                                                                                                             |
| H 16 2 O |                                                                                                                                      |                                                                                                                                                                      |                                                                                                                                                 |                                                                                                                                                 |

Notes. ( a ) Observed intensity uncertainties are given in parentheses.

edge-on PDR. Using an isobaric model, the best-fit to our bestestimate H 16 2 O intensities is obtained for a cloud depth corresponding to a visual extinction of A tot V = 20, a UV field of G 0 = 3 5 : GLYPH&lt;2&gt; 10 4 GLYPH&lt;0&gt; 8 GLYPH&lt;2&gt; 10 4 in Mathis units and a thermal pressure Pth = 6 GLYPH&lt;2&gt; 10 7 GLYPH&lt;0&gt; 2 GLYPH&lt;2&gt; 10 8 Kcm 3 . The 211 GLYPH&lt;0&gt; GLYPH&lt;0&gt; 202 and 202 GLYPH&lt;0&gt; 111 line intensities require a UV field 40% lower ( G 0 = 3 GLYPH&lt;2&gt; 10 4 in Mathis units) to match the observations, whereas the H 16 2 O 212 GLYPH&lt;0&gt; 101 line requires a higher pressure by a factor of five. It is worth noting that the intensity derived for this line is 20% higher than the one reported by Nagy et al. (2017) measured from a short integration-time observation . H 8 18 2 O lines are not used in the fitting procedure because the online ISMDB does not include H 18 2 Opredictions. Input parameters for the ISMDB model of the Meudon PDR code are listed in Table 3.

The best-fit parameters obtained for the observed H 16 2 O intensities are close to those derived from high-J CO lines (Joblin et al. 2018) or small hydrocarbons and complex organic

8 H 16 2 O 212 GLYPH&lt;0&gt; 101-b, Obs. Id: 1342218426

Fig. 3: Comparison between observed intensities and the ones predicted by the pure gas-phase chemistry ISMDB model and the adsorption and desorption optimized model with H2 internal energy (Table 3). Bottom panel gives the relative di GLYPH&lt;11&gt; erence between predictions and observations.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-92ffa022-c532-491e-9f8f-37fbe818c924.jpg)

molecules (COM) emission (Cuadrado et al. 2015, 2017). Small deviations from the high-J CO line fit could be explained by the fact that water emits GLYPH&lt;24&gt; 5 00 deeper in the cloud than CO (Parikka et al. 2017). However these deviations may also result from the data reduction method applied. Main-beam e GLYPH&lt;14&gt; ciencies reported by Roelfsema et al. (2012) were used by Joblin et al. (2018) along with a mean value between the antenna temperature and the main-beam brightness temperature to deal with a spatial emission in-between point-source and extended source (Ossenkopf et al. 2013). This method could lead to an underestimate of the main-beam brightness temperature by a factor close to 1.3. Moreover specific parameters of the code have been optimized for high-J CO fitting. In particular, the grain size distribution was adapted to model the extinction curve toward Θ 1 Ori C (Fitzpatrick &amp; Massa 1990; Marconi et al. 1998; Joblin et al. 2018), and a scaling factor was applied to correct the bias for the assumed geometry of the Bar, such as the beam coupling factor and the viewing angle, and get closer to an edge-on model.

The stand-alone version of the Meudon PDR code (1.5.2) was then used with the best-fit parameters derived from the water line observations to compute the intensities of water transitions and in particular those of H 18 2 O. This is performed by extending the gas-phase H 16 2 O chemistry to H 18 2 O, introducing the possible

Fig. 4: Evolution of the the water ortho-to-para ratio (left axis) in agreement with the kinetic temperature (right axis) as a function of the visual extinction inside the molecular cloud for the pure gas-phase chemistry ISMDB model and the adsorption and desorption optimized model with H2 internal energy (Table 3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-0c390f0e-9d22-4515-bce9-bc2f2c4827d5.jpg)

18 Ofractionation reactions (Loison et al. 2019). The radiative excitation processes involving H 18 2 O have been implemented from the HITRAN2016 9 database (Gordon et al. 2017). The H 18 2 Odeexcitation collisional rates are assumed to be identical to those of H 16 2 O. The upward collisional excitation rates of H 18 2 O are computed by introducing the appopriate energy defect for H 18 2 O and considering the di GLYPH&lt;11&gt; erent energy order involving highly-J excited rotational levels. Collisional rates for rotational excitation of H 16 2 Oby He are taken from Green et al. (1993) for which analytical temperature dependence was derived by Gonzalez Garcia et al. (2008). The collisional rates due to H2 are computed from these values by introducing the relative reduced mass factor. The computed intensities are given in Table 4 and compared with the corrected observed line intensities in Fig. 3.

The observed intensities are quite well reproduced by the Meudon PDR code with these parameters. Lines outside the observed uncertainty range are those mentioned above (212 GLYPH&lt;0&gt; 101, 202 GLYPH&lt;0&gt; 111, and 211 GLYPH&lt;0&gt; 202) and the H 18 2 O ortho ground state line, which is clearly overestimated by the model. The value obtained for the H 18 2 O 111 GLYPH&lt;0&gt; 000 para line agrees with the upper limit observed.

In the Meudon PDR code, water vapor is produced with nuclear spin populations in agreement with the local kinetic temperature. The computed kinetic temperature varies from 2000 K close to the ionization front down to 30 K in the deepest part of the cloud (see Fig. 4). The ortho-to-para ratio of water as a function of the depth into the cloud is shown in Fig. 4. To obtain the mean OPR of the water molecules along the line-of-sight, the OPR is averaged along the maximum density region of water vapor. Figure 5 shows that the density of gas-phase water peaks between A V = 5 and A V = 10. Thus the gas-phase water-emitting

9 https: // hitran.org /

Fig. 5: Comparison of the spatial density of gas-phase water (left axis) and H density (right axis) as a function of the visual extinction inside the molecular cloud for the pure gas-phase chemistry ISMDB model and the adsorption and desorption optimized model with H2 internal energy (Table 3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-fe20dbdd-79cb-43fa-8acb-7351dc528cf7.jpg)

region is characterized by a temperature of (34.0 GLYPH&lt;6&gt; 0.5) K, implying an OPR of 2.7 GLYPH&lt;6&gt; 0.1. Computing the OPR with the ratio of ortho and para column densities of gas-phase water in the cloud up to A V = 15 leads to a value in good agreement within the uncertainty range. However this value is clearly higher than the value of 0.1 to 0.5 derived by Choi et al. (2014) from H 18 2 O observations.

The OPR uncertainty derived above only considers the predicted OPR variations through the cloud. To take into account the uncertainty due to discrepancies between the observations and the predicted intensities, models at the boundary UV field and thermal pressure regions giving a good fit lead to an OPR uncertainty of 0.1.

## 4.2. Adsorption and desorption model

We stress that gas-phase water population and spin distribution in the currently released version of the Meudon PDR code is only governed by gas-phase chemistry and thermal equilibrium. The presence of grains in the model only a GLYPH&lt;11&gt; ects the thermal balance through photoelectric heating, while attenuating the UV flux (Le Petit et al. 2006) and catalysing the formation of H2 via surface chemistry (Le Petit et al. 2009; Le Bourlot et al. 2012). The need of surface chemistry was less important in previous studies, for which tracers that appear in a warmer medium (Fig. 6) were used (Cuadrado et al. 2015, 2017; Joblin et al. 2018). Up to A V = 3, Esplugues et al. (2016) have shown, with a PDR model considering adsorption and desorption processes and surface chemistry, that the gas-phase abundances of several species, such as H2O, H2CO, and CH3OH, are independant of the dust surface chemistry. Indeed, close to the ionization front, the high kinetic temperature of the gas and grains (see Fig. 7) is not suitable for long time residence of these molecules on

Fig. 6: Evolution of the gas-phase rotational populations of H 16 2 O (left axis), CO and 13 CO (right axis) as a function of the visual extinction inside the molecular cloud for the adsorption and desorption optimized model with H2 internal energy (Table 3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-b2a3d98d-7ee9-46b2-9381-7269b7701b34.jpg)

grains (Hollenbach et al. 2009). However, for water molecules, surface chemistry could have an important role to play. Indeed the oxygen surface chemistry, introduced in the Hollenbach et al. (2009) PDR model, appears to dominate the gas-phase water production. Figures 4 and 5 show that the gas-phase water maximumdensity region from A V = 5 to 10 corresponds to a medium with a gas temperature below 40 K. In addition, the water vapor abundance remains relatively high for A V &gt; 10, where a strong freeze-out of water molecules on grains is expected (Hollenbach et al. 2009; Esplugues et al. 2016). Thus the line emission of water molecules could be strongly modified by their interactions with cold grains.

The implementation of state-of-the-art surface chemistry processes in the Meudon PDR code is in progress. For this work, we have added adsorption and desorption processes of water. The adsorption is governed by a sticking coe GLYPH&lt;14&gt; cient proportional to the inverse of the square root of the temperature for a kinetic temperature higher than 10 K and equals to 1 otherwise. Two desorption processes were considered. The thermal desorption is dependent on a binding energy of 5600 K (Garrod et al. 2009; Wakelam et al. 2017) and the photodesorption is related to a yield of 5 GLYPH&lt;2&gt; 10 = 4 molecules per incident UV photon ( ¨ berg O et al. 2009).

Using our ISMDB best-fit parameters and adsorption and desorption processes, the overall water density was computed with the Meudon PDR code. Depletion of gas-phase water is clearly observed for A V &gt; 5 when the adsorption process is included, in good agreement with Hollenbach et al. (2009). As expected, the depletion of water molecules decreases the computed intensity of water lines with an accentuated e GLYPH&lt;11&gt; ect for the excited lines. The depletion is su GLYPH&lt;14&gt; ciently significant that very high UV flux and thermal pressure are needed to match the observational data using regular parameters in the Meudon PDR code.

Fig. 7: Evolution of the fractional abundance of H 16 2 O gas and solid-phase (left axis) as a function of the visual extinction inside the molecular cloud for the adsorption and desorption optimized model with H2 internal energy (Table 3). Dashed lines show the gas-phase temperature and the smallest (SG GLYPH&lt;24&gt; 3 nm) and biggest (BG GLYPH&lt;24&gt; 300 nm) grain temperature (right axis).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-804c49b1-690a-4883-985b-30e6718993c1.jpg)

Gas-phase water is mainly formed by two processes (van Dishoeck et al. 2013 and references therein) represented by each water density maximum in Fig. 5. The first peak around A V = 1 is the result of the balance between water formation via OH and photodissociation.

$$\colon r \, \varphi \, \text{or} \quad \text{OH} + H _ { 2 } \to H _ { 2 } O + H \\ \ a \, \text{strong} \quad H _ { 2 } O + \text{photon} \to H + \text{OH}$$

Then the gas-phase OH abundance decreases and the photodissociation reduces the gas-phase water density. From A V = 3 to A V = 5, the increase in water vapor abundance is caused by the recombination of H3O + with electrons, still softened by the photodissociation.

$$\mathfrak { f } \text{ water} \colon \ H _ { 3 } O ^ { + } + e ^ { - } \to H _ { 2 } O + H$$

Adsorption and desorption processes compete with H3O + recombination leading to an increase in the water ice abundance up to A V = 10 as represented in Fig. 7. Eventually, the electronic recombination is weakened and the gas-phase water abundance is reduced.

Without water surface chemistry, grains act as water trap which explains the discrepancy between observed and computed intensities. To reproduce the observed intensities, water should be formed or released from the grains more e GLYPH&lt;14&gt; ciently.

Close to the H H2 / transition, an activation energy of 3240 K has to be overcome to initiate the following reaction and increase the water precursor reservoir (van Dishoeck et al. 2013 and references therein).

$$\overset { \, \rightharpoonup } { \, } \, \text{excte} \, \quad O + H _ { 2 } \to O H + H$$

In the standard 1.5.2 version of the Meudon PDR code, energy for chemical reactions is only provided by the kinetic energy of

the reactants due to thermal motion. However, close to the H H2 / transition, H2 is highly ro-vibrationally excited (Hollenbach &amp; Tielens 1997). This energy can be used to overcome the activation energy (Ag´ undez et al. 2010), as was already studied for the C + + H2 reaction (Zanchet et al. 2013b; Herr´ aez-Aguilar et al. 2014) and the S + + H2 reaction (Zanchet et al. 2013a). To implement this e GLYPH&lt;11&gt; ect in the model, the activation energy is thus taken as the di GLYPH&lt;11&gt; erence between the activation energy and the rovibrational energy of H2. Reactions with vibrationally excited H2 were already included in the Meudon PDR code to reproduce the OH emission in the Orion Bar (Goicoechea et al. 2011). Considering this e GLYPH&lt;11&gt; ect enhances the intensity predicted for the excited lines up to a factor of 2.5 due to the increase in the excited populations between A V = 6 and A V = 8.

Deeper into the cloud, where the gas is colder, gas-phase water is formed by the recombination of H3O + with electrons. This reaction is initiated by cosmic ray ionization producing H + 3 and leading to H3O + through ion-neutral reactions (Gerin et al. 2010). Moreover, in the cold region, where water is depleted by adsorption, desorption is mainly achieved by photodesorption triggered by the UV secondary photons induced by cosmic rays (Prasad &amp; Tarafdar 1983). Following these two processes, the cosmic ray ionization rate should have a main impact on line intensities coming from the deepest part of the cloud. From our H column density ( N H) estimate in the water region, a range of cosmic ray ionization rate can be estimated. For A V &gt; 5, N H is higher than 10 22 cm = 2 and the cosmic ray ionization rate could be increased by a factor of ten in relation to the value used by Joblin et al. (2018), up to GLYPH&lt;16&gt; H2 = 5 GLYPH&lt;2&gt; 10 GLYPH&lt;0&gt; 16 s GLYPH&lt;0&gt; 1 per H2 (Padovani et al. 2018). This leads to an increase in the line intensities by a factor of the order of three, except for H 16 2 O 211 GLYPH&lt;0&gt; 202 for which a factor of seven is obtained. Such an enhancement is mainly produced by the increase in the H3O + reservoir.

Finally, the water photodesorption yield used in the PDR Code could be tuned up. Recent study of Cruz-Diaz et al. (2018) reported photodesorption yield up to 2 GLYPH&lt;2&gt; 10 = 3 molecule per photon at a temperature around 60 K, using a microwave discharged hydrogen flow lamp with a strong LyGLYPH&lt;11&gt; emission component. Assuming this value increases the water line intensities by a factor of the order of two.

Considering the H2 internal energy, taking upper limit for the cosmic ray ionization rate and the water photodesorption yield, and using input parameters derived by Joblin et al. (2018) otherwise, the resulting intensities for adsorption and desorption model are in good agreement with the observed intensities (see Fig. 3 and Table 4). Input parameters of the adsorption and desorption optimized model are listed in Table 3. We note that the scaling factor used by Joblin et al. (2018) to correct bias for the assumed geometry of the Bar was not applied in this study.

With this model, the depletion of water by adsorption on grains is compensated in the maximum gas-phase water density region as seen in Fig. 5. Figure 7 shows that the water ice abundance reaches a plateau at a fractional abundance ( f H2O = n H2O = n H) of 2 GLYPH&lt;2&gt; 10 = 4 for A V &gt; 10 leading to a decrease in the gas-phase water fractional abundance which peaks at 2 GLYPH&lt;2&gt; 10 = 7 between A V = 7 and A V = 10. This gas-grains balance in relation to the depth agrees with what was predicted by Hollenbach et al. (2009) with oxygen surface chemistry. However, the fractional abundance derived from our best model is one order of magnitude higher than those reported in this former work for a UV field of G 0 = 1 GLYPH&lt;2&gt; 10 3 and a density of 1 GLYPH&lt;2&gt; 10 4 cm = 3 . Water ice abundance increases from A V = 3 where the gas temperature is below 100 K and the grain temperature below 40 K. The

Fig. 8: Evolution of the line center opacity of H 16 2 O and CO lines as a function of the visual extinction inside the molecular cloud for the adsorption and desorption optimized model with H2 internal energy (Table 3).

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-4f46cc43-d87f-45fe-8dd9-48d9a7c9108f.jpg)

grain temperature evolution with depth is consistent with dust observations (Arab et al. 2012).

Figure 6 shows the populations of gas-phase water rotational levels, which peaks between A V = 5 to A V = 10 without stratification. The averaged ortho-to-para ratio in this region is equal to 2.8 GLYPH&lt;6&gt; 0.2 corresponding to a spin temperature of (36 GLYPH&lt;6&gt; 2) K (see Fig. 4) for a UV field of G 0 = 3 1 : GLYPH&lt;2&gt; 10 4 in Mathis units and a thermal pressure Pth = 2 8 : GLYPH&lt;2&gt; 10 8 Kcm 3 . The spin temperature GLYPH&lt;0&gt; derived from this model agrees with the gas-phase kinetic temperature expected by Hollenbach et al. (2009) for G 0 = 1 GLYPH&lt;2&gt; 10 4 and n H = 1 GLYPH&lt;2&gt; 10 5 cm GLYPH&lt;0&gt; 3 .

## 5. Discussion

## 5.1. Self-absorption

The analysis of water line intensities with the Meudon PDR code brings a justification for the assumption made to derive the observed intensities. The presence of a central dip in line profiles is explained by the self-absorption by foreground gas (Choi 2015). In this scenario, the flux coming from layers deep inside the cloud is absorbed by the foreground layers when escaping the Bar.

The Meudon PDR code prediction of the line center opacities, defined as e GLYPH&lt;0&gt; GLYPH&lt;28&gt; with GLYPH&lt;28&gt; the line center optical depth, are given in Fig. 8. This figure shows that the 110 GLYPH&lt;0&gt; 101, 111 GLYPH&lt;0&gt; 000, and 212 GLYPH&lt;0&gt; 101 transitions, which exhibit the deepest self-absorption, are by far optically thicker than the other lines. We have verified that the excitation temperature of the ground state water lines increases rapidly with A V near the location of the GLYPH&lt;28&gt; = 1 surface ( A V GLYPH&lt;24&gt; 1). The opacity is the highest at the line center and decreases going into the wings. An absorption dip at the central velocity is therefore expected. Thus line profiles predicted by the Meudon PDR code readily reproduce such a dip for the

three most intense and broadest lines, while predicting single peak profiles for the other lines.

## 5.2. H 18 2 O intensities

This work has led to the first prediction of H 18 2 O line intensities by the Meudon PDR code. The comparison with observed intensities in Fig. 3 indicates that the computed intensities are overestimated. Furthermore, our estimated gas-phase column density for H 18 2 Oortho ground state line is two times higher than the one reported by Choi et al. (2014).

It appears that the gas-phase water 16 O 18 O isotopic ratio is / lower than the typical value of 560 reported by Wilson &amp; Rood (1994) for the local ISM and based on H2CO surveys (Gardner &amp; Whiteoak 1981). Moreover, the decreasing of the ratio is accentuated in the water emission region down to 350. This could be explained by the fact that only half of the 52 oxygen bearing molecules considered by the Meudon PDR code have their oxygen isomers taken into account. Thus a larger 18 O reactants reservoir, in relation to 16 O chemistry, is available to form H 18 2 O. Considering the 16 O 18 O isotopic ratio / of 560, and assuming that the intensity is impacted by the difference between the isotopic ratios in the main water emission region, it appears that the computed H 18 O 2 line intensities are overestimated by a factor of 1.6. Applying this correction, the prediction of the H 18 2 O ortho 110 GLYPH&lt;0&gt; 101 intensity goes down to 6.3 GLYPH&lt;2&gt; 10 = 8 erg cm = 2 s = 1 sr = 1 in good agreement with the observed intensity. The para 111 GLYPH&lt;0&gt; 110 intensity predicted by the model would thus be around 2.3 GLYPH&lt;2&gt; 10 = 7 erg cm = 2 s = 1 sr = 1 . Removing this spurious e GLYPH&lt;11&gt; ect requires experimental data for the rare isotope molecules.

## 5.3. OPR

Water line intensities are well reproduced using the Meudon PDR code, considering either pure gas-phase chemistry model, or by adding adsorption and desorption of water from grains. In both cases, an ortho-to-para ratio close to 2.8 is derived, whereas a previous estimate of the H 18 2 O OPR in the Orion Bar leads to a very low value between 0.1 and 0.5 (Choi et al. 2014). This low ratio was estimated from the column densities obtained with the RADEX low velocity gradient (LVG) code (van der Tak et al. 2007), assuming a single-slab geometry with a homogeneous density and temperature. As already discussed, the gasphase H 18 2 O column densities, derived from the Meudon PDR code using our adsorption and desorption optimized model, are probably overestimated due to the departure of the isotopic ratio from the typical local ISM value. Considering the ISM value of Wilson &amp; Rood (1994), the gas-phase H 18 O column densi2 ties up to A V = 15 are GLYPH&lt;24&gt; 8.2 GLYPH&lt;2&gt; 10 9 cm = 2 for the para 111 level and GLYPH&lt;24&gt; 2.5 GLYPH&lt;2&gt; 10 10 cm = 2 for the ortho 110 level. The para value is one order of magnitude lower than the one previously reported whereas the ortho value agrees with it (Choi et al. 2014). Using these two values, the H 18 2 O OPR is estimated to be equal to 3.0.

Part of the discrepancy between the low OPR reported and our high thermal value derived could be assigned to the improvement of the data reduction pipeline. Indeed, with the latest HIPE version (14.1.0), the H 18 2 O para ground state line, which has a really low signal-to-noise ratio, is barely detected and our upper limit for its intensity is roughly two times lower than the Choi et al. (2014) value 10 . It is worth noting that Nagy et al. (2017)

10 HIPE version 10.0 was used

have not reported the detection of this line, unlike the H 18 2 Oortho ground state transition.

However a rough estimate of the OPR in LTE, using the H 18 2 O ortho line intensity and the para upper limit still gives a low OPR of 0.7. Moreover, considering the two H 18 2 O intensities predicted by the Meudon PDR code, the OPR in LTE would be equal to 1.3. Thus the model used for the analysis has a major impact on the results. RADEX LVG code (van der Tak et al. 2007), used to derive the H 18 2 O OPR, assumes an isothermal homogeneous medium and solves the radiative transfer in non-LTE using radiative and collisional transition rates. Conversely, the Meudon PDR code derives the population considering chemistry, thermal balance and infrared pumping. With this more complete model and assuming an isobaric model, the kinetic temperature, and thus the density, is predicted to be strongly dependent on the depth into the cloud as shown in Figs. 4 and 5. Alongside with an expanded dataset, the analysis of water lines intensity in the Orion Bar is consistent with the thermal equilibrium at 36 K.

To derive the water OPR, we computed an average over the main gas-phase water reservoir between A V = 5 and A V = 10. This assumption is justified by the spatial distribution of the water line emission in the Orion Bar. Indeed OH and H2O emission seems to be decorrelated whereas OH emission matches high-J COones (Goicoechea et al. 2011; Parikka et al. 2017). As high-J CO lines arise from a layer between A V = 1 and A V = 3 (Joblin et al. 2018), the local water density maximum at A V = 1 GLYPH&lt;0&gt; 2 should not significantly contribute to the observed water intensities.

This OPR agrees with OPRs at high thermal values from other tracers reported for the Orion Bar. Indeed an OPR of 2.8 GLYPH&lt;6&gt; 0.6 was derived for c GLYPH&lt;0&gt; C3H2 (Cuadrado et al. 2015) and OPRs of the order of three were inferred for H2CO, H2CS and H2CCO (Hogerheijde et al. 1995; Cuadrado et al. 2017).

This work follows the trend of water OPR reanalysis initiated by better reduction pipeline for archival observations along with the development of more sophisticated models. Thus fewer very low OPRs are confirmed, as for example in the TW Hydrae protoplanetary disk, where the very low OPR reported by Hogerheijde et al. (2011) has been shown to be modeldependent and is consistent within the uncertainties with the high-temperature limit (Salinas et al. 2016).

## 5.4. Grain surface chemistry

The introduction of simple adsorption and desorption processes of water has emphasized the major role that grains could have for tracers of the deepest part of the cloud. To overcome the depletion produced by the adsorption of water molecules on grains, several parameters of the code had to be tuned up to their upper acceptable limits. However a careful study of other parameters could mitigate the extreme values adopted for the cosmic ray ionization rate and the water photodesorption yield.

First, considering the presence of X-ray photons in the Orion Nebula (Getman et al. 2005; Preibisch et al. 2005) could lead to a reduction of these parameters. The chemistry initiated by cosmic ray ionization, such as gas-phase water production via H3O + , could also be activated by X-ray ionization (Gupta et al. 2010; Cuadrado et al. 2015). Furthermore, water X-ray photodesorption appears to be e GLYPH&lt;14&gt; cient and comparable to far-UV photodesorption (Dupuy et al. 2018).

On the other hand, the 60 GLYPH&lt;14&gt; limitation of the Meudon PDR code could lower the predicted intensity in relation to a more

edge-on geometry. Indeed, previous studies have estimated the ionization front to be tilted from the line-of-sight by an angle between 3 GLYPH&lt;14&gt; and 20 GLYPH&lt;14&gt; (Wen &amp; O'Dell 1995; Hogerheijde et al. 1995; Jansen et al. 1995; Walmsley et al. 2000) implying a nearly edgeon PDR. Joblin et al. (2018) have estimated that a scaling factor of 1.3 has to be applied to correct for the assumed geometries of the Bar from a model at a viewing angle of 60 GLYPH&lt;14&gt; and to fit highJ CO lines. A more edge-on PDR would be characterized by a scaling factor higher than 1.

Moreover, the grain size distribution adopted by Joblin et al. (2018) is optimized to reproduce the extinction curve toward Θ 1 Ori C (Fitzpatrick &amp; Massa 1990; Marconi et al. 1998). This distribution is kept constant through the cloud, whereas deeper in grains could be larger. Thus in the water-emitting region, the total surface of the grains for water molecules to stick onto would be reduced preventing gas-phase water density from dropping down dramatically.

Furthermore, adding only water adsorption and desorption processes on grains could introduce a bias on the chemistry balance and the radiative transfer computation. Indeed, if a full adsorption and desorption description is adopted, modifying the precursors reservoir by trapping them onto grains could change the chemical state of the cloud. This would also have an impact on the radiative transfer by limiting the UV shielding of the molecules that appear deeper into the cloud, and on the desorption processes that should be enhanced by the increase in the number of UV photons inside the cloud. Moreover, increasing the number of molecules on grains will emphasize the need of considering a full surface chemistry balance. Some PDR models have introduced surface chemistry calculations (Hollenbach et al. 2009; Esplugues et al. 2016) and raised the importance of surface chemistry to explain observed gas-phase abundances. In these models, the main ice water formation pathway begins with oxygen adsorption on grains followed by two reactions with hydrogen atoms. Hollenbach et al. (2009) have estimated that in the highest gas-phase water abundance region, 97 % of gas or solid-phase water production is achieved by surface chemistry. However, owing to the high FUV field, the grain temperature is high enough to prevent O or OH to stick e GLYPH&lt;14&gt; ciently (Hollenbach et al. 2009; Melnick et al. 2012), as seen in Fig. 7. Thus the formation of water ice through this pathway should be reduced.

Implementing these processes in the Meudon PDR code is in progress. As already mentioned for isotopic computations, it is only achievable by gathering qualitative and quantitative processes from experimental data. For example adsorption and desorption considerations are mainly limited by the lack of photodesorption yield for astrophysical relevant molecules. In the past few years, several laboratory experiments have been conducted to describe photodesorption processes ( ¨ berg et al. 2009; O Mu˜ noz Caro et al. 2016; Dupuy et al. 2017). Simulations of water surface chemistry have shown the importance of surface processes for the interstellar gas-phase (Cazaux et al. 2010). The development of experimental set-up to address surface chemistry (Hama &amp; Watanabe 2013 and references therein) unveiled the various processes that could occur on interstellar grains and a GLYPH&lt;11&gt; ect the physical and chemical state of interstellar objects.

Getting experimental data to deal with nuclear spin e GLYPH&lt;11&gt; ects is crucial to simulate and understand nuclear spin population distributions. The nuclear spin conversion mechanism has been investigated to estimate the lifetime conservation of a potential thermal disequilibrium (Fillion et al. 2012; Cacciani et al. 2012; Turgeon et al. 2017), as well as the e GLYPH&lt;11&gt; ect of adsorption and desorption on the OPR (Hama et al. 2016). Furthermore selective reaction or formation of water nuclear spin isomers could also be taken into account by the models (Kilaj et al. 2018). Considering these processes are the steps to examine the assumptions of water produced at the thermal equilibrium in the Meudon PDR code, used to fit the water observations of the Orion Bar.

## 6. Conclusions

An analysis of a set of seven H 16 2 Olines and one H 18 2 Oline measured toward the Orion Bar with Herschel / HIFI is performed using the Meudon PDR code. The para H 18 O ground state line is 2 not unambiguously detected. Considering an isobaric model and assuming that gas-phase water is formed with an ortho-to-para ratio in agreement with the local thermal equilibrium, an OPR of 2.7 GLYPH&lt;6&gt; 0.1 is derived.

Investigation of grain e GLYPH&lt;11&gt; ects is made by adding adsorption and desorption processes of water on grains. After tuning the cosmic ray ionization rate and water photodesorption yield up to their upper acceptable limits, a good agreement with the observed intensity is obtained for an OPR of 2.8 GLYPH&lt;6&gt; 0.2. This orthoto-para ratio corresponds to a temperature of (36 GLYPH&lt;6&gt; 2) K, which is much higher than the spin temperature previously reported for H 18 2 O in the Orion Bar and close to the high-temperature limit.

Acknowledgements. This work was supported by the Programme National de Physique et Chimie du Milieu Interstellaire (PCMI) of the CNRS INSU , the / INC INP co-funded by the CEA and the CNES. Financial support from the / LabEx MiChem, par of the French state funds managed by the ANR within the investissements d'avenir program under reference ANR-11-10EX-000402 is acknowledged. Fundings from the Ile-de-France region DIM ACAV + is acknowledged. Part of this research was carried out at the Jet Propulsion Laboratory, California Institute of Technology, under a contract with the National Aeronautics and Space Administration.

## References

Ag´ undez, M., Goicoechea, J. R., Cernicharo, J., Faure, A., &amp; Roue GLYPH&lt;11&gt; , E. 2010, ApJ, 713, 662

- Allers, K. N., Ja GLYPH&lt;11&gt; e, D. T., Lacy, J. H., Draine, B. T., &amp; Richter, M. J. 2005, ApJ, 630, 368
- Andree-Labsch, S., Ossenkopf-Okada, V., &amp; R¨ ollig, M. 2017, A&amp;A, 598, A2
- Arab, H., Abergel, A., Habart, E., et al. 2012, A&amp;A, 541, A19
- Bergin, E. A., Phillips, T. G., Comito, C., et al. 2010, A&amp;A, 521, L20
- Bron, E., Ag´ undez, M., Goicoechea, J. R., &amp; Cernicharo, J. 2018, arXiv e-prints, arXiv:1801.01547
- Cacciani, P., Cosl´ eou, J., &amp; Khelkhal, M. 2012, Phys. Rev. A, 85, 012521 Cazaux, S., Cobut, V., Marseille, M., Spaans, M., &amp; Caselli, P. 2010, A&amp;A, 522, A74

Choi, Y. 2015, PhD thesis, Groningen

- Choi, Y., van der Tak, F. F. S., Bergin, E. A., &amp; Plume, R. 2014, A&amp;A, 572, L10 Cruz-Diaz, G. A., Mart´ ın-Dom´ enech, R., Moreno, E., Mu˜ noz Caro, G. M., &amp; Chen, Y.-J. 2018, MNRAS, 474, 3080
- Cuadrado, S., Goicoechea, J. R., Cernicharo, J., et al. 2017, A&amp;A, 603, A124

Cuadrado, S., Goicoechea, J. R., Pilleri, P., et al. 2015, A&amp;A, 575, A82

- de Graauw, T., Helmich, F. P., Phillips, T. G., et al. 2010, A&amp;A, 518, L6

Dupuy, R., Bertin, M., F´ eraud, G., et al. 2018, Nat. Astron., 2, 796

- Dupuy, R., Bertin, M., F´ eraud, G., et al. 2017, A&amp;A, 603, A61
- Esplugues, G. B., Cazaux, S., Meijerink, R., Spaans, M., &amp; Caselli, P. 2016, A&amp;A, 591, A52
- Fillion, J.-H., Bertin, M., Lekic, A., et al. 2012, EAS Pub. Ser., 58, 307 Fitzpatrick, E. L. &amp; Massa, D. 1990, ApJS, 72, 163
- Fuente, A., Rodriguez-Franco, A., &amp; Martin-Pintado, J. 1996, A&amp;A, 312, 599 Gardner, F. F. &amp; Whiteoak, J. B. 1981, MNRAS, 194, 37P
- Garrod, R. T., Vasyunin, A. I., Semenov, D. A., Wiebe, D. S., &amp; Henning, T. 2009, ApJ, 700, L43
- Gerin, M., De Luca, M., Black, J., et al. 2010, A&amp;A, 518, L110
- Getman, K. V., Feigelson, E. D., Grosso, N., et al. 2005, ApJS, 160, 353
- Goicoechea, J. R., Cuadrado, S., Pety, J., et al. 2017, A&amp;A, 601, L9
- Goicoechea, J. R., Joblin, C., Contursi, A., et al. 2011, A&amp;A, 530, L16
- Goicoechea, J. R., Pety, J., Cuadrado, S., et al. 2016, Nature, 537, 207
- , E. 2008, A&amp;A, 485,
- Goicoechea, J. R., Teyssier, D., Etxaluze, M., et al. 2015, ApJ, 812, 75 Gonzalez Garcia, M., Le Bourlot, J., Le Petit, F., &amp; Roue GLYPH&lt;11&gt; 127

Gordon, I. E., Rothman, L. S., Hill, C., et al. 2017, J. Quant. Spec. Radiat. Transf., 203, 3

Green, S., Maluendes, S., &amp; McLean, A. D. 1993, ApJS, 85, 181 Gupta, H., Rimmer, P., Pearson, J. C., et al. 2010, A&amp;A, 521, L47 Habart, E. 2011, OT2 ehabart 4: Unveiling the origin and excitation mechanisms of the warm CO, OH and CH + , Herschel Space Observatory Proposal, id.2039

Habart, E., Dartois, E., Abergel, A., et al. 2010, A&amp;A, 518, L116 Habing, H. J. 1968, Bull. Astron. Inst. Netherlands, 19, 421 Hama, T., Kouchi, A., &amp; Watanabe, N. 2016, Science, 351, 65 Hama, T. &amp; Watanabe, N. 2013, Chem. Rev., 113, 8783 Herr´ aez-Aguilar, D., Jambrina, P. G., Men´ endez, M., et al. 2014,

Phys. Chem. Chem. Phys., 16, 24800 Hogerheijde, M. R., Bergin, E. A., Brinch, C., et al. 2011, Science, 334, 338 Hogerheijde, M. R., Jansen, D. J., &amp; van Dishoeck, E. F. 1995, A&amp;A, 294, 792 Hollenbach, D. J., Kaufman, M. J., Bergin, E. A., &amp; Melnick, G. J. 2009, ApJ, 690, 1497

Hollenbach, D. J. &amp; Tielens, A. G. G. M. 1997, ARA&amp;A, 35, 179 Jansen, D. J., Spaans, M., Hogerheijde, M. R., &amp; van Dishoeck, E. F. 1995, A&amp;A, 303, 541

Joblin, C., Bron, E., Pinto, C., et al. 2018, A&amp;A, 615, A129 Kilaj, A., Gao, H., Rsch, D., et al. 2018, Nat. Commun., 9 Kong, S., Arce, H. G., Feddersen, J. R., et al. 2018, ApJS, 236, 25 Larsson, B., Liseau, R., Bergman, P., et al. 2003, A&amp;A, 402, L69 Le Bourlot, J., Le Petit, F., Pinto, C., Roue GLYPH&lt;11&gt; , E., &amp; Roy, F. 2012, A&amp;A, 541, A76 Le Petit, F., Barzel, B., Biham, O., Roue GLYPH&lt;11&gt; , E., &amp; Le Bourlot, J. 2009, A&amp;A, 505, 1153

Le Petit, F., Nehm´ e, C., Le Bourlot, J., &amp; Roue

GLYPH&lt;11&gt;

, E. 2006, ApJS, 164, 506

Leurini, S., Rol

GLYPH&lt;11&gt;

s, R., Thorwirth, S., et al. 2006, A&amp;A, 454, L47

Lis,

D.

C.,

Bergin,

E.

A.,

Schilke,

P.,

&amp; van

Dishoeck,

J. Phys. Chem. A, 117, 9661

Lis, D. C. &amp; Schilke, P. 2003, ApJ, 597, L145

Loison, J.-C., Wakelam, V., Gratier, P., et al. 2019, MNRAS, 485, 5777

Marconi, A., Testi, L., Natta, A., &amp; Walmsley, C. M. 1998, A&amp;A, 330, 696

Mathis, J. S., Mezger, P. G., &amp; Panagia, N. 1983, A&amp;A, 128, 212

Melnick, G. J., Tolls, V., Goldsmith, P. F., et al. 2012, ApJ, 752, 26

Menten, K. M., Reid, M. J., Forbrich, J., &amp; Brunthaler, A. 2007, A&amp;A, 474, 515

Mu˜ noz Caro, G. M., Chen, Y.-J., Aparicio, S., et al. 2016, A&amp;A, 589, A19

Mueller, M., Jellema, W., Olberg, M., Moreno, R., &amp; Teyssier, D. 2014, The

HIFI Beam: Release 1 - Release Note for Astronomers, Tech. Rep. HIFI-

ICC-RP-2014-001, SRON Gronigen

Nagy, Z., Choi, Y., Ossenkopf-Okada, V., et al. 2017, A&amp;A, 599, A22 ¨ berg, K. I., Linnartz, H., Visser, R., &amp; van Dishoeck, E. F. 2009, ApJ, 693, O 1209

Ossenkopf, V., R¨ ollig, M., Neufeld, D. A., et al. 2013, A&amp;A, 550, A57 Ott, S. 2010, in Astronomical Data Analysis Software and Systems XIX, Vol.

- 434, ASP Conf. Ser., ed. Y. Mizumoto, K.-I. Morita, &amp; M. Ohishi, 139 Padovani, M., Galli, D., Ivlev, A. V., Caselli, P., &amp; Ferrara, A. 2018, A&amp;A, 619, A144

Parikka, A., Habart, E., Bernard-Salas, J., et al. 2017, A&amp;A, 599, A20 Parmar, P. S., Lacy, J. H., &amp; Achtermann, J. M. 1991, ApJ, 372, L25 Pilbratt, G. L., Riedinger, J. R., Passvogel, T., et al. 2010, A&amp;A, 518, L1 Prasad, S. S. &amp; Tarafdar, S. P. 1983, ApJ, 267, 603

Preibisch, T., McCaughrean, M. J., Grosso, N., et al. 2005, ApJS, 160, 582 Rodriguez-Franco, A., Martin-Pintado, J., &amp; Fuente, A. 1998, A&amp;A, 329, 1097 Roelfsema, P. R., Helmich, F. P., Teyssier, D., et al. 2012, A&amp;A, 537, A17 Salinas, V. N., Hogerheijde, M. R., Bergin, E. A., et al. 2016, A&amp;A, 591, A122 Shipman, R. F., Beaulieu, S. F., Teyssier, D., et al. 2017, A&amp;A, 608, A49 Sternberg, A. &amp; Dalgarno, A. 1989, ApJ, 338, 197

Stoerzer, H., Stutzki, J., &amp; Sternberg, A. 1995, A&amp;A, 296, L9

Tauber, J. A., Tielens, A. G. G. M., Meixner, M., &amp; Goldsmith, P. F. 1994, ApJ, 422, 136

Teyssier, D., Avruch, I., Beaulieu, S., et al. 2017, HIFI Handbook, HERSCHELHSC-DOC-2097

Tielens, A. G. G. M. &amp; Hollenbach, D. 1985, ApJ, 291, 747

- Turgeon, P.-A., Vermette, J., Alexandrowicz, G., et al. 2017, J. Phys. Chem. A, 121, 1571, pMID: 28157310

van der Tak, F. F. S., Black, J. H., Sch¨ier, F. L., Jansen, D. J., &amp; van Dishoeck, o E. F. 2007, A&amp;A, 468, 627

van der Tak, F. F. S., Nagy, Z., Ossenkopf, V., et al. 2013, A&amp;A, 560, A95 van der Tak, F. F. S., Ossenkopf, V., Nagy, Z., et al. 2012, A&amp;A, 537, L10 van der Werf, P. P., Stutzki, J., Sternberg, A., &amp; Krabbe, A. 1996, A&amp;A, 313, 633

van Dishoeck, E. F., Herbst, E., &amp; Neufeld, D. A. 2013, Chem. Rev., 113, 9043 Wakelam, V., Loison, J.-C., Mereau, R., &amp; Ruaud, M. 2017, Mol. Astrophys., 6, 22

Walmsley, C. M., Natta, A., Oliva, E., &amp; Testi, L. 2000, A&amp;A, 364, 301 Wen, Z. &amp; O'Dell, C. R. 1995, ApJ, 438, 784

E.

F.

2013,

Wilson, T. L. &amp; Rood, R. 1994, ARA&amp;A, 32, 191 Wu, R., Bron, E., Onaka, T., et al. 2018, A&amp;A, 618, A53 Wyrowski, F., Schilke, P., Hofner, P., &amp; Walmsley, C. M. 1997, ApJ, 487, L171 Young Owl, R. C., Meixner, M. M., Wolfire, M., Tielens, A. G. G. M., &amp; Tauber,

J. 2000, ApJ, 540, 886

- Zanchet, A., Ag´ undez, M., Herrero, V. J., Aguado, A., &amp; Roncero, O. 2013a, AJ, 146, 125

Zanchet, A., Godard, B., Bulut, N., et al. 2013b, ApJ, 766, 80

## Appendix A: Observed lines and calibrations

Table A.1: Line intensities integrated with a linear baseline fit outside the line, except for H 18 2 O110 GLYPH&lt;0&gt; 101 where the frequency switch observing mode leads to optical standing waves that had to be corrected by a fouth order polynom. See Table 1 for the definition of a and b labeled transitions.

| Transition            | Pol.   |   HPBW a (') | GLYPH<17> mb a      |   RA-DEC b (') |   RA-DEC b (') | ∆ V int c (kms = 1 )   | R T mb d V (Kkms = 1 )   | ∆ V Line d (kms = 1 )   | Ω e                |   Cal. f (%) |
|-----------------------|--------|--------------|---------------------|----------------|----------------|------------------------|--------------------------|-------------------------|--------------------|--------------|
| GLYPH<0>              | H      |        37.15 | 0.622 GLYPH<6> 0.07 |            1.6 |            4.8 | [6.7-14.5]             | 21.02 GLYPH<6> 0.22      | 4.3 GLYPH<6> 0.3        | 0.74 GLYPH<6> 0.11 |          5.8 |
| 1 10 1 01 -a          | V      |        37.49 | 0.616 GLYPH<6> 0.07 |           -0.2 |           -1.6 | [6.7-14.5]             | 22.58 GLYPH<6> 0.23      | 4.6 GLYPH<6> 0.3        | 0.79 GLYPH<6> 0.10 |          5.8 |
| 1 10 GLYPH<0> 1 01 -b | H      |        37.15 | 0.622 GLYPH<6> 0.07 |            0.9 |            3.4 | [6.7-14.5]             | 20.37 GLYPH<6> 0.24      | 4.6 GLYPH<6> 0.5        | 0.75 GLYPH<6> 0.11 |          5.3 |
|                       | V      |        37.49 | 0.616 GLYPH<6> 0.07 |           -2.7 |           -2.2 | [6.7-14.5]             | 20.56 GLYPH<6> 0.25      | 4.3 GLYPH<6> 0.8        | 0.80 GLYPH<6> 0.09 |          5.3 |
| 1 GLYPH<0> 0          | H      |        18.88 | 0.632 GLYPH<6> 0.08 |           -1   |            1.5 | [6.7-14.5]             | 9.63 GLYPH<6> 0.44       | 4.4 GLYPH<6> 0.8        | 0.78 GLYPH<6> 0.09 |          5.4 |
| 11 00                 | V      |        18.71 | 0.636 GLYPH<6> 0.08 |            2.2 |            0   | [6.7-14.5]             | 10.93 GLYPH<6> 0.45      | 4.3 GLYPH<6> 0.5        | 0.85 GLYPH<6> 0.05 |          5.8 |
|                       | H      |        12.58 | 0.593 GLYPH<6> 0.07 |            0.2 |            0.8 | [7.2-14.0]             | 10.41 GLYPH<6> 0.46      | 5.2 GLYPH<6> 0.4        | 0.78 GLYPH<6> 0.07 |          7.1 |
| 2 12 GLYPH<0> 1 01 -a | V      |        12.41 | 0.572 GLYPH<6> 0.07 |            0.4 |            0.1 | [7.2-14.0]             | 10.10 GLYPH<6> 0.47      | 4.9 GLYPH<6> 1.1        | 0.82 GLYPH<6> 0.07 |          6.9 |
| O                     | H      |        12.58 | 0.593 GLYPH<6> 0.07 |           -1.9 |           -0.3 | [7.2-14.0]             | 6.76 GLYPH<6> 1.33       | 4.9 GLYPH<6> 1.4        | 0.80 GLYPH<6> 0.06 |          5.1 |
| 2 12 GLYPH<0> 1 01 -b | V      |        12.41 | 0.572 GLYPH<6> 0.07 |           -2   |            0.5 | [7.2-14.0]             | 8.38 GLYPH<6> 1.77       | 4.0 GLYPH<6> 0.3        | 0.76 GLYPH<6> 0.08 |          4.8 |
|                       | H      |        21.28 | 0.634 GLYPH<6> 0.08 |           -1.2 |            2   | [7.7-13.5]             | 6.47 GLYPH<6> 0.28       | 2.9 GLYPH<6> 0.8        | 0.78 GLYPH<6> 0.09 |          5.4 |
| 2 02 GLYPH<0> 1 11    | V      |        21.09 | 0.637 GLYPH<6> 0.08 |            1.8 |            0.1 | [7.7-13.5]             | 9.64 GLYPH<6> 0.38       | 3.0 GLYPH<6> 0.8        | 0.84 GLYPH<6> 0.06 |          5.7 |
| 2 11 GLYPH<0> 2 02    | H      |        28    | 0.640 GLYPH<6> 0.09 |           -0.4 |            3   | [8.2-13.0]             | 6.51 GLYPH<6> 0.26       | 2.6 GLYPH<6> 0.6        | 0.77 GLYPH<6> 0.11 |          5.5 |
|                       | V      |        27.91 | 0.659 GLYPH<6> 0.09 |            0.9 |           -1.3 | [8.2-13.0]             | 9.13 GLYPH<6> 0.22       | 2.4 GLYPH<6> 0.4        | 0.82 GLYPH<6> 0.08 |          5.7 |
| 2 GLYPH<0> 2          | H      |        12.65 | 0.579 GLYPH<6> 0.07 |            0.2 |           10.8 | [8.2-13.0]             | 2.04 GLYPH<6> 0.45       | 1.4 GLYPH<6> 1.0        | 0.78 GLYPH<6> 0.07 |          7.1 |
| 21 12                 | V      |        12.48 | 0.572 GLYPH<6> 0.07 |            0.4 |            0.1 | [8.2-13.0]             | 3.20 GLYPH<6> 0.55       | 2.4 GLYPH<6> 0.5        | 0.82 GLYPH<6> 0.07 |          6.9 |
| 3 12 GLYPH<0> 3 03    | H      |        19.16 | 0.632 GLYPH<6> 0.08 |           -1   |            1.5 | [9.2-12.0]             | 1.43 GLYPH<6> 0.24       | 2.1 GLYPH<6> 0.6        | 0.78 GLYPH<6> 0.09 |          5.4 |
|                       | V      |        18.98 | 0.636 GLYPH<6> 0.08 |            2.2 |            0   | [9.2-12.0]             | 2.27 GLYPH<6> 0.23       | 1.9 GLYPH<6> 0.6        | 0.85 GLYPH<6> 0.05 |          5.8 |
|                       | H      |        37.77 | 0.622 GLYPH<6> 0.07 |           -0.7 |           -1.5 | [8.2-13.0]             | 0.31 GLYPH<6> 0.02       | 1.9 GLYPH<6> 0.5        | 0.79 GLYPH<6> 0.10 |          5.2 |
| 1 10 GLYPH<0> 1 01    | V      |        38.12 | 0.616 GLYPH<6> 0.07 |            1.9 |            4.6 | [8.2-13.0]             | 0.26 GLYPH<6> 0.01       | 1.9 GLYPH<6> 0.3        | 0.74 GLYPH<6> 0.11 |          5.3 |
| O                     | H      |        19.08 | 0.632 GLYPH<6> 0.08 |           -1   |            1.5 | [8.2-13.0]             | 0.43 GLYPH<6> 0.30       |                         | 0.78 GLYPH<6> 0.09 |          5.4 |
| 1 11 GLYPH<0> 0 00    | V      |        18.91 | 0.636 GLYPH<6> 0.08 |            2.2 |            0   | [8.2-13.0]             | 0.25 GLYPH<6> 0.30       |                         | 0.85 GLYPH<6> 0.05 |          5.8 |

Notes. ( a ) Half-power beam width and main-beam e GLYPH&lt;14&gt; ciency of the Herschel telescope at the water line frequency (Mueller et al. 2014). ( b ) Right ascension and declination o GLYPH&lt;11&gt; sets from the CO + peak ( GLYPH&lt;11&gt; 2000 = 05 h 35 m 20 61 s : ; GLYPH&lt;14&gt; 2000 = GLYPH&lt;0&gt; 05 25 14 0 GLYPH&lt;14&gt; 0 : 00 ). ( c ) Integration velocity range. ( d ) Observed full width at half maximum of the lines. ( e ) Beam coupling factor. ( f ) Absolute calibration uncertainty.

Fig. A.1: Line profile of H 16 2 O for H (left panel) and V (right panel) polarization. The blue solid line is the spectrum corrected for HIFI calibration, the black solid line is the spectrum corrected for the beam coupling factor, the blue dash-dotted line and red dashed line represent Orion Bar and Orion Ridge velocity features and the gray area is the velocity integration range. See Tables 1 and A.1 for the definition of a and b labeled transitions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-18e22b80-e89f-49f7-a314-05032d2019a8.jpg)

Fig. A.2: Line profile of H 16 2 O for H (left panel) and V (right panel) polarization. The blue solid line is the spectrum corrected for HIFI calibration, the black solid line is the spectrum corrected for the beam coupling factor, the blue dash-dotted line and red dashed line lines represent Orion Bar and Orion Ridge velocity features and the gray area is the velocity integration range. See Tables 1 and A.1 for the definition of a and b labeled transitions.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-0bdb2e18-82c2-4e86-a7d5-b9e0b4f1cc66.jpg)

Fig. A.3: Line profile of H 16 2 O for H (left panel) and V (right panel) polarization. The blue solid line is the spectrum corrected for HIFI calibration, the black solid line is the spectrum corrected for the beam coupling factor, the blue dash-dotted line0gray area is the velocity integration range.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-011009f0-085c-4db7-883c-5b209f0d4582.jpg)

Fig. A.4: Line profile of H 18 2 O for H (left panel) and V (right panel) polarization. The blue solid line is the spectrum corrected for HIFI calibration, the black solid line is the spectrum corrected for the beam coupling factor, the blue dash-dotted line and red dashed line represent Orion Bar and Orion Ridge velocity features and the gray area is the velocity integration range.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-894fff75-20c4-4d64-8ae6-3a00a993028b.jpg)

Fig. A.5: Water emission maps averaged for H and V polarizations. Integration ranges are given in Table A.1. Black solid contours represent the emission of 13 CO 6 GLYPH&lt;0&gt; 5 at 95%, 90%, 75% and 50% of the maximum intensity integrated from 2 km s = 1 to 17 km s = 1 observed with the CSO telescope and convolved at the averaged water line H and V beam width (Table A.1). Upward and downward triangles mark the targeted observations for H and V polarizations (blue ones are for the b labeled transitions, see Tables 1 and A.1). Black dashed circle is the HIFI HPBW at the water frequency centered on the CO + peak.

![image](https://ac.marcodsn.me/arxiv-markdown-images/1908.00268v2-figure-bd2d206e-e93f-4088-af61-653420c45923.jpg)
T. Putaud X. Michaut F. Le Petit E. Roueff D. C. Lis
The water line emission and ortho-to-para ratio in the Orion Bar photon-dominated region
A very low ortho-to-para ratio (OPR) of 0.1-0.5 was previously reported in
the Orion Bar photon-dominated region (PDR), based on observations of two
optically thin $\mathrm{H_2^{18}O}$ lines which were analyzed by using a
single-slab large velocity gradient model. The corresponding spin temperature
does not coincide with the kinetic temperature of the molecular gas in this
UV-illuminated region. This was interpreted as an indication of water molecules
being formed on cold icy grains which were subsequently released by UV
photodesorption. A more complete set of water observations in the Orion Bar,
including seven $\mathrm{H_2^{16}O}$ lines and one $\mathrm{H_2^{18}O}$ line,
carried out using Herschel/HIFI instrument, was reanalyzed using the Meudon PDR
code to derive gas-phase water abundance and the OPR, taking into account the
steep density and temperature gradients present in the region. The model line
intensities are in good agreement with the observations assuming that water
molecules formed with an OPR corresponding to thermal equilibrium conditions at
the local kinetic temperature of the gas and when solely considering gas-phase
chemistry and water gas-grain exchanges through adsorption and desorption.
Gas-phase water is predicted to arise from a region deep into the cloud,
corresponding to a visual extinction of $A_{\mathrm{V}} \sim 9$, with a
$\mathrm{H_2^{16}O}$ fractional abundance of $\sim 2\times 10^{-7}$ and column
density of $(1.4 \pm 0.8) \times 10^{15}$ cm$^{-2}$ for a total cloud depth of
$A_{\mathrm{V}}=15$. A line-of-sight average ortho-to-para ratio of $2.8 \pm
0.2$ is derived. The observational data are consistent with a nuclear spin
isomer repartition corresponding to the thermal equilibrium at a temperature of
$36 \pm 2$ K, much higher than the spin temperature previously reported for
this region and close to the gas kinetic temperature in the water-emitting gas.