CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

(abridged) The far-UV wavelength range (912-2000A) provides access to atomic and molecular transitions of many species the interstellar medium (ISM), circumgalactic medium (CGM), and intergalactic medium, within phases spanning a wide range of ionization, density, temperature, and molecular gas fraction. Far-UV space telescopes have enabled detailed studies of the ISM in the Milky Way thanks to absorption features appearing in the UV spectra of hot stars and yielding fundamental insights into the composition and physical characteristics of all phases of the ISM along with the processes that influence them. However, we have yet to design a spectrometer able to observe the full UV domain at resolving power R>10^5 with a signal-to-noise ratio SNR>500. Such a resolution is necessary to resolve lines from both the cold molecular hydrogen and the warm metal ions with a turbulent velocity of about 1 km s-1, and to differentiate distinct velocity components. Future UV spectroscopic studies of the Milky Way ISM must revolutionize our understanding of the ISM as a dynamical, unstable, and magnetized medium, and rise to the challenge brought forward by current theories. Another interesting prospect is to transpose the same level of details that has been reached for the Milky Way to the ISM in external galaxies, in particular in metal-poor galaxies, where the ISM chemical composition, physical conditions, and topology change dramatically, with significant consequences on the star-formation properties. Finally, we need to be able to perform statistical analyses of background quasar lines of sight intersecting the CGM of galaxies at various redshifts and to comprehend the role of gas exchanges and flows for galaxy evolution.

Filaments are a key step on the path that leads from molecular clouds to star formation. However, their characteristics are heavily debated, and the exact processes that lead to their formation and fragmentation into dense cores still remain to be fully understood. We aim at characterising the mass, kinematics, and stability against gravitational collapse of a statistically significant sample of filaments in the Orion B molecular cloud, which is renown for its very low star formation efficiency. We characterise the gas column densities and kinematics over a field of 1.9 deg$^2$, using C$^{18}$O(J=1-0) data from the IRAM-30m large programme ORION-B. Using two different Hessian-based filters, we extract and compare two filamentary networks, each containing over 100 filaments. Independent of the extraction method, the filaments have widths of 0.12$\pm$0.04 pc, and show a wide range of linear (1 - 100 $M_{\odot}$pc$^{-1}$) and volume densities (2.10$^3$ - 2.10$^5$ cm$^{-3}$). Compared to previous studies, the filament population is dominated by low-density, thermally sub-critical structures, suggesting that most of the identified filaments are not collapsing to form stars. In fact, only ~1% of the Orion B cloud mass covered by our observations can be found in super-critical, star-forming filaments, explaining the low star formation efficiency of the region. The velocity profiles observed across the filaments show quiescence in the centre, and coherency in the plane of the sky, despite being mostly supersonic. The filaments in Orion B apparently belong to a continuum which contains a few elements comparable to already studied star-forming filaments as well as many lower-density, gravitationally unbound structures. This comprehensive study of the Orion B filaments shows that the mass fraction in super-critical filaments is a key factor in determining star formation efficiency.

The radiative and mechanical interaction of stars with their environment drives the evolution of the ISM and of galaxies as a whole. The far-IR emission (lambda ~30 to 350 microns) from atoms and molecules dominates the cooling of the warm gas in the neutral ISM, the material that ultimately forms stars. Far-IR lines are thus the most sensitive probes of stellar feedback processes, and allow us to quantify the deposition and cycling of energy in the ISM. While ALMA (in the (sub)mm) and JWST (in the IR) provide astonishing sub-arcsecond resolution images of point sources and their immediate environment, they cannot access the main interstellar gas coolants, nor are they designed to image entire star-forming regions (SFRs). Herschel far-IR photometric images of the interstellar dust thermal emission revealed the ubiquitous large-scale filamentary structure of SFRs, their mass content, and the location of thousands of prestellar cores and protostars. These images, however, provide a static view of the ISM: not only they dont constrain the cloud dynamics, moreover they cannot reveal the chemical composition and energy transfer within the cloud, thus giving little insight into the regulation process of star formation by stellar feedback. In this white paper we emphasize the need of a space telescope with wide-field spectral-imaging capabilities in the critical far-IR domain.

Young massive stars regulate the physical conditions, ionization, and fate of their natal molecular cloud. It is important to find tracers that help quantifying the stellar feedback processes that take place at different scales. We present ~85 arcmin^2 velocity-resolved maps of several submm molecular lines toward the closest high-mass star-forming region, OMC-1. The observed rotational lines include probes of warm and dense molecular gas that are difficult to detect from ground-based telescopes: CH+ (1-0), CO (10-9), HCO+ (6-5), and HCN (6-5). These lines trace an extended but thin layer of molecular gas at high thermal pressure, P_th ~ 1e7-1e9 K/cm3, associated with the FUV-irradiated surface of OMC-1. The intense FUV field, emerging from massive stars in the Trapezium cluster, heats, compresses and photoevaporates the cloud edge. It also triggers the formation of reactive molecules such as CH+. The CH+ (1-0) emission spatially correlates with the flux of FUV photons impinging the cloud: G_0 from 1e3 to 1e5. This correlation is supported by isobaric PDR models in the parameter space P_th/G_0 ~ [5e3-8e4] K/cm3 where many PDRs seem to lie. The CH+ (1-0) emission correlates with the extended emission from vibrationally excited H2, and with that of [CII]158um and CO 10-9, all emerging from FUV-irradiated gas. These correlations link the presence of CH+ to the availability of C+ ions and of FUV-pumped H2(v>0) molecules. The parsec-scale CH+ emission and narrow-line (dv ~ 3 km/s) mid-J CO emission arises from extended PDRs and not from fast shocks. PDR line tracers are the smoking gun of the stellar feedback from young massive stars. The PDR component in OMC-1 represents 5 to 10% of the total gas mass, however, it dominates the emitted line luminosity. These results provide insights into the source of submm CH+ and mid-J CO emission from distant star-forming galaxies.

In bright photodissociation regions (PDRs) associated to massive star formation, the presence of dense "clumps" that are immersed in a less dense interclump medium is often proposed to explain the difficulty of models to account for the observed gas emission in high-excitation lines. We aim at presenting a comprehensive view of the modeling of the CO rotational ladder in PDRs, including the high-J lines that trace warm molecular gas at PDR interfaces. We observed the 12CO and 13CO ladders in two prototypical PDRs, the Orion Bar and NGC 7023 NW using the instruments onboard Herschel. We also considered line emission from key species in the gas cooling of PDRs (C+, O, H2) and other tracers of PDR edges such as OH and CH+. All the intensities are collected from Herschel observations, the literature and the Spitzer archive and are analyzed using the Meudon PDR code. A grid of models was run to explore the parameter space of only two parameters: thermal gas pressure and a global scaling factor that corrects for approximations in the assumed geometry. We conclude that the emission in the high-J CO lines, which were observed up to Jup=23 in the Orion Bar (Jup=19 in NGC7023), can only originate from small structures of typical thickness of a few 1e-3 pc and at high thermal pressures (Pth~1e8 K cm-3). Compiling data from the literature, we found that the gas thermal pressure increases with the intensity of the UV radiation field given by G0, following a trend in line with recent simulations of the photoevaporation of illuminated edges of molecular clouds. This relation can help rationalising the analysis of high-J CO emission in massive star formation and provides an observational constraint for models that study stellar feedback on molecular clouds.

We investigate the physical conditions of the CO gas near the young star cluster, Trumpler 14 of the Carina Nebula. The observations presented in this work are taken with the Fourier Transform Spectrometer (FTS) of the Spectral and Photometric Imaging REceiver (SPIRE) onboard the Herschel Space Observatory. Our field of view covers the edge of a cavity carved by Trumpler 14 about $1\,\mathrm{Myr}$ ago and marks the transition from HII regions to photo-dissociation regions. With the state-of-the-art Meudon PDR code, we successfully derive the physical conditions, which include the thermal pressure ($P$) and the scaling factor of radiation fields ($G_{\mathrm{UV}}$), from the observed CO spectral line energy distributions~(SLEDs) in the observed region. The derived $G_{\mathrm{UV}}$ values generally show an excellent agreement with the UV radiation fields created by nearby OB-stars and thus confirm that the main excitation source of the observed CO emission are the UV-photons provided by the massive stars. The derived thermal pressure is between $0.5-3\,\times\,10^{8}\,\mathrm{K\,cm^{-3}}$ with the highest values found along the ionization front in Car I-E region facing Trumpler 14, hinting that the cloud structure is similar to the recent observations of the Orion Bar. Comparing the derived thermal pressure with the radiation fields, we report the first observationally-derived and spatially-resolved $P \sim 2\times10^4\,G_{\mathrm{UV}}$ relationship. As direct comparisons of the modeling results to the observed $^{13}\mathrm{CO}$, [OI] $63\,\mathrm{\mu m}$, and [CII] $158\,\mathrm{\mu m}$ intensities are not straightforward, we urge the readers to be cautious when constraining the physical conditions of PDRs with combinations of $^{12}\mathrm{CO}$, $^{13}\mathrm{CO}$, [CI], [OI] $63\,\mathrm{\mu m}$, and [CII] $158\,\mathrm{\mu m}$ observations.

Recent Herschel and ALMA observations of Photodissociation Regions (PDRs) have revealed the presence of a high thermal pressure (P ~ 10^7-10^8 K cm-3) thin compressed layer at the PDR surface where warm molecular tracer emission (e.g. CH+, SH+, high-J CO, H2,...) originate. These high pressures (unbalanced by the surrounding environment) and a correlation between pressure and incident FUV field (G0) seem to indicate a dynamical origin with the radiation field playing an important role in driving the dynamics. We investigate whether photoevaporation of the illuminated edge of a molecular cloud could explain these high pressures and pressure-UV field correlation. We developed a 1D hydrodynamical PDR code coupling hydrodynamics, EUV and FUV radiative transfer and time-dependent thermo-chemical evolution. We applied it to a 1D plane-parallel photoevaporation scenario where a UV-illuminated molecular cloud can freely evaporate in a surrounding low-pressure medium. We find that photoevaporation can produce high thermal pressures and the observed P-G0 correlation, almost independently from the initial gas density. In addition, we find that constant-pressure PDR models are a better approximation to the structure of photoevaporating PDRs than constant-density PDR models, although moderate pressure gradients are present. Strong density gradients from the molecular to the neutral atomic region are found, which naturally explain the large density contrasts (1-2 orders of magnitude) derived from observations of different tracers. The photoevaporating PDR is preceded by a low velocity shock (a few km/s) propagating into the molecular cloud. Photoevaporating PDR models offer a promising explanation to the recent observational evidence of dynamical effects in PDRs.

Molecular hydrogen is the most abundant molecule in the universe. It is the first one to form and survive photo-dissociation in tenuous environments. Its formation involves catalytic reactions on the surface of interstellar grains. The micro-physics of the formation process has been investigated intensively in the last 20 years, in parallel of new astrophysical observational and modeling progresses. In the perspectives of the probable revolution brought by the future satellite JWST, this article has been written to present what we think we know about the H$_2$ formation in a variety of interstellar environments.

Previous attempts at segmenting molecular line maps of molecular clouds have focused on using position-position-velocity data cubes of a single line to separate the spatial components of the cloud. In contrast, wide field spectral imaging with large spectral bandwidth in the (sub)mm domain now allows to combine multiple molecular tracers to understand the different physical and chemical phases that constitute giant molecular clouds. We aim at using multiple tracers (sensitive to different physical processes) to segment a molecular cloud into physically/chemically similar regions (rather than spatially connected components). We use a machine learning clustering method (the Meanshift algorithm) to cluster pixels with similar molecular emission, ignoring spatial information. Simple radiative transfer models are used to interpret the astrophysical information uncovered by the clustering. A clustering analysis based only on the J=1-0 lines of 12CO, 13CO and C18O reveals distinct density/column density regimes (nH~100, 500, and >1000 cm-3), closely related to the usual definitions of diffuse, translucent and high-column-density regions. Adding two UV-sensitive tracers, the (1-0) lines of HCO+ and CN, allows us to distinguish two clearly distinct chemical regimes, characteristic of UV-illuminated and UV-shielded gas. The UV-illuminated regime shows overbright HCO+ and CN emission, which we relate to photochemical enrichment. We also find a tail of high CN/HCO+ intensity ratio in UV-illuminated regions. Finer distinctions in density classes (nH~7E3, and 4E4 cm-3) for the densest regions are also identified, likely related to the higher critical density of the CN and HCO+ (1-0) lines. The association of simultaneous multi-line, wide-field mapping and powerful machine learning methods such as the Meanshift algorithm reveals how to decode the complex information available in molecular tracers.

The carriers of the diffuse interstellar bands (DIBs) are largely unidentified molecules ubiquitously present in the interstellar medium (ISM). After decades of study, two strong and possibly three weak near-infrared DIBs have recently been attributed to the C60+ fullerene based on observational and laboratory measurements. There is great promise for the identification of the over 400 other known DIBs, as this result could provide chemical hints towards other possible carriers. In an effort to systematically study the properties of the DIB carriers, we have initiated a new large-scale observational survey: the ESO Diffuse Interstellar Bands Large Exploration Survey (EDIBLES). The main objective is to build on and extend existing DIB surveys to make a major step forward in characterising the physical and chemical conditions for a statistically significant sample of interstellar lines-of-sight, with the goal to reverse-engineer key molecular properties of the DIB carriers. EDIBLES is a filler Large Programme using the Ultraviolet and Visual Echelle Spectrograph at the Very Large Telescope at Paranal, Chile. It is designed to provide an observationally unbiased view of the presence and behaviour of the DIBs towards early-spectral-type stars whose lines-of-sight probe the diffuse-to-translucent ISM. Such a complete dataset will provide a deep census of the atomic and molecular content, physical conditions, chemical abundances and elemental depletion levels for each sightline. Achieving these goals requires a homogeneous set of high-quality data in terms of resolution (R ~ 70000 -- 100000), sensitivity (S/N up to 1000 per resolution element), and spectral coverage (305--1042 nm), as well as a large sample size (100+ sightlines). In this first paper the goals, objectives and methodology of the EDIBLES programme are described and an initial assessment of the data is provided.

Far-UV photons (E<13.6 eV) from hot massive stars regulate, or at least influence, the heating, ionization, and chemistry of most of the neutral interstellar medium. Investigating the interaction between FUV radiation and interstellar matter thus plays an important role in astrochemistry. We have used ALMA to mosaic a small field of the Orion Bar where the critical transition from atomic to molecular gas takes place. These observations provide an unprecedented sharp view of this transition layer (~1" resolution or ~414 AU). The resulting images (so far in the rotational emission of CO, HCO+, H13CO+, SO+, SO, and reactive ions SH+ and HOC+) show the small-scale structure in gas density and temperature, and the steep abundance gradients. The images reveal a pattern of high-density substructures, photo-ablative gas flows and instabilities at the edge of the molecular cloud. We have also used the IRAM 30m telescope to carry out a line-survey of the illuminated edge of the Bar in the mm domain. Our observations reveal the presence of complex organic molecules that were not expected in such a harsh environment. In particular, we have reported the first detection of the unstable cis conformer of formic acid (HCOOH) in the ISM. The energy barrier to internal rotation (the conversion from trans to cis) is approximately 4827 cm-1. Hence, this detection is surprising. The low inferred trans-to-cis abundance ratio of 2.8+/-1.0 supports a photoswitching mechanism: a given conformer absorbs a FUV stellar photon that radiatively excites the molecule to electronic states above the interconversion barrier. Subsequent fluorescent decay leaves the molecule in a different conformer form. This mechanism, which we have specifically studied with ab initio quantum calculations, was not considered so far in astrochemistry although it can affect the structure of a variety of molecules in PDRs. (abridged)

Observations towards L1630 in the Orion B molecular cloud, comprising the iconic Horsehead Nebula, allow us to study the interplay between stellar radiation and a molecular cloud under relatively benign conditions, that is, intermediate densities and an intermediate UV radiation field. Contrary to the well-studied Orion Molecular Cloud 1 (OMC1), which hosts much harsher conditions, L1630 has little star formation. We aim to relate the [CII] fine-structure line emission to the physical conditions predominant in L1630 and compare it to studies of OMC1. The [CII] $158\,\mu\mathrm{m}$ emission from an area of $12' \times 17'$ in L1630 was observed using the upGREAT instrument onboard SOFIA. Of the [CII] emission from the mapped area 95%, $13\,L_{\odot}$, originates from the molecular cloud; the adjacent HII region contributes only 5%, that is, $1\,L_{\odot}$. From comparison with other data (CO (1-0)-line emission, far-infrared (FIR) continuum studies, emission from polycyclic aromatic hydrocarbons (PAHs)), we infer a gas density of the molecular cloud of $n_{\mathrm{H}}\sim 3\cdot 10^3\,\mathrm{cm^{-3}}$, with surface layers, including the Horsehead Nebula, having a density of up to $n_{\mathrm{H}}\sim 4\cdot 10^4\,\mathrm{cm^{-3}}$. The temperature of the surface gas is $T\sim 100\,\mathrm{K}$. The average [CII] cooling efficiency within the molecular cloud is $1.3\cdot 10^{-2}$. The fraction of the mass of the molecular cloud within the studied area that is traced by [CII] is only $8\%$. Our PDR models are able to reproduce the FIR-[CII] correlations and also the CO (1-0)-[CII] correlations. Finally, we compare our results on the heating efficiency of the gas with theoretical studies of photoelectric heating by PAHs, clusters of PAHs, and very small grains, and find the heating efficiency to be lower than theoretically predicted, a continuation of the trend set by other observations.

We report high angular resolution (4.9" x 3.0") images of reactive ions SH+, HOC+, and SO+ toward the Orion Bar photodissociation region (PDR). We used ALMA-ACA to map several rotational lines at 0.8 mm, complemented with multi-line observations obtained with the IRAM 30m telescope. The SH+ and HOC+ emission is restricted to a narrow layer of 2"- to 10"-width (~800 to 4000 AU depending on the assumed PDR geometry) that follows the vibrationally excited H2^* emission. Both ions efficiently form very close to the H/H2 transition zone, at a depth of A_V < 1 mag into the neutral cloud, where abundant C+, S+, and H2^* coexist. SO+ peaks slightly deeper into the cloud. The observed ions have low rotational temperatures (T_rot~10-30 K << T_k) and narrow line-widths (~2-3 km/s), a factor of ~2 narrower that those of the lighter reactive ion CH+. This is consistent with the higher reactivity and faster radiative pumping rates of CH+ compared to the heavier ions, which are driven relatively faster toward smaller velocity dispersion by elastic collisions and toward lower T_rot by inelastic collisions. We estimate column densities and average physical conditions from a non-LTE excitation model (n(H2)~10^5-10^6 cm^-3, n(e^-)~10 cm^-3, and T_k~200 K). Regardless of the excitation details, SH+ and HOC+ clearly trace the most exposed layers of the UV-irradiated molecular cloud surface, whereas SO+ arises from slightly more shielded layers.

Context. The combination of wideband receivers and spectrometers currently available in (sub-)millimeter observatories deliver wide- field hyperspectral imaging of the interstellar medium. Tens of spectral lines can be observed over degree wide fields in about fifty hours. This wealth of data calls for restating the physical questions about the interstellar medium in statistical terms. Aims. We aim at gaining information on the physical structure of the interstellar medium from a statistical analysis of many lines from different species over a large field of view, without requiring detailed radiative transfer or astrochemical modeling. Methods. We coupled a nonlinear rescaling of the data with one of the simplest multivariate analysis methods, namely the Principal Component Analysis, to decompose the observed signal into components that we interpret first qualitatively and then quantitatively based on our deep knowledge of the observed region and of the astrochemistry at play. Results. We identify 3 principal components, linear compositions of line brightness temperatures, that are correlated at various levels with the column density, the volume density and the UV radiation field. Conclusions. When sampling a sufficiently diverse mixture of physical parameters, it is possible to decompose the molecular emission in order to gain physical insight on the observed interstellar medium. This opens a new avenue for future studies of the interstellar medium.

The nature of turbulence in molecular clouds is one of the key parameters that control star formation efficiency: compressive motions, as opposed to solenoidal motions, can trigger the collapse of cores, or mark the expansion of Hii regions. We try to observationally derive the fractions of momentum density ($\rho v$) contained in the solenoidal and compressive modes of turbulence in the Orion B molecular cloud and relate these fractions to the star formation efficiency in the cloud. The implementation of a statistical method developed by Brunt & Federrath (2014), applied to a $^{13}$CO(J=1-0) datacube obtained with the IRAM-30m telescope, allows us to retrieve 3-dimensional quantities from the projected quantities provided by the observations, yielding an estimate of the compressive versus solenoidal ratio in various regions of the cloud. Despite the Orion B molecular cloud being highly supersonic (mean Mach number $\sim$ 6), the fractions of motion in each mode diverge significantly from equipartition. The cloud's motions are on average mostly solenoidal (excess > 8 % with respect to equipartition), which is consistent with its low star formation rate. On the other hand, the motions around the main star-forming regions (NGC 2023 and NGC 2024) prove to be strongly compressive. We have successfully applied to observational data a method that was so far only tested on simulations, and have shown that there can be a strong intra-cloud variability of the compressive and solenoidal fractions, these fractions being in turn related to the star formation efficiency. This opens a new possibility for star-formation diagnostics in galactic molecular clouds.

We aim to develop the Orion B Giant Molecular Cloud (GMC) as a local template for interpreting extra-galactic molecular line observations. We use the wide-band receiver at the IRAM-30m to spatially and spectrally resolve the Orion B GMC. The observations cover almost 1 square degree at 26" resolution with a bandwidth of 32 GHz from 84 to 116 GHz in only two tunings. Among the mapped spectral lines are the 12CO, 13CO, C18O, C17O, HCN, HNC, 12CN, CCH, HCO+, N2H+ (1-0), and 12CS, 32SO, SiO, c-C3H2, CH3OH (2-1) transitions. We introduce the molecular anatomy of the Orion B GMC, including relations between line intensities and gas column density or far-UV radiation fields, and correlations between selected line and line ratios. We also obtain a dust-traced gas mass that is less than about one third the CO-traced mass, using the standard Xco conversion factor. The presence of overluminous CO can be traced back to the dependence of the CO intensity on UV illumination. In fact, while most lines show some dependence on the UV radiation field, CN and CCH are the most sensitive. Moreover dense cloud cores are almost exclusively traced by N2H+. Other traditional high density tracers, such as HCN (1-0), are also easily detected in extended translucent regions at a typical density of about 500 H2 cm-3. In general, we find no straightforward relation between line critical density and the fraction of the line luminosity coming from dense gas regions. Our initial findings demonstrate that the relations between line (ratio) intensities and environment in GMCs are more complicated than often assumed. Sensitivity (i.e., the molecular column density), excitation, and above all chemistry contribute to the observed line intensity distributions. They must be considered together when developing the next generation of extra-galactic molecular line diagnostics of mass, density, temperature and radiation field.

We present Herschel SPIRE Fourier Transform Spectrometer (FTS) observations of N159W, an active star-forming region in the Large Magellanic Cloud (LMC). In our observations, a number of far-infrared cooling lines including CO(4-3) to CO(12-11), [CI] 609 and 370 micron, and [NII] 205 micron are clearly detected. With an aim of investigating the physical conditions and excitation processes of molecular gas, we first construct CO spectral line energy distributions (SLEDs) on 10 pc scales by combining the FTS CO transitions with ground-based low-J CO data and analyze the observed CO SLEDs using non-LTE radiative transfer models. We find that the CO-traced molecular gas in N159W is warm (kinetic temperature of 153-754 K) and moderately dense (H2 number density of (1.1-4.5)e3 cm-3). To assess the impact of the energetic processes in the interstellar medium on the physical conditions of the CO-emitting gas, we then compare the observed CO line intensities with the models of photodissociation regions (PDRs) and shocks. We first constrain the properties of PDRs by modelling Herschel observations of [OI] 145, [CII] 158, and [CI] 370 micron fine-structure lines and find that the constrained PDR components emit very weak CO emission. X-rays and cosmic-rays are also found to provide a negligible contribution to the CO emission, essentially ruling out ionizing sources (ultraviolet photons, X-rays, and cosmic-rays) as the dominant heating source for CO in N159W. On the other hand, mechanical heating by low-velocity C-type shocks with ~10 km/s appears sufficient enough to reproduce the observed warm CO.

Context: Fast surface conversion between ortho- and para-H2 has been observed in laboratory studies, and this mechanism has been proposed to play a role in the control of the ortho-para ratio in the interstellar medium. Observations of rotational lines of H2 in Photo-Dissociation Regions (PDRs) have indeed found significantly lower ortho-para ratios than expected at equilibrium. The mechanisms controlling the balance of the ortho-para ratio in the interstellar medium thus remain incompletely understood, while this ratio can affect the thermodynamical properties of the gas (equation of state, cooling function). Aims: We aim to build an accurate model of ortho-para conversion on dust surfaces based on the most recent experimental and theoretical results, and to validate it by comparison to observations of H2 rotational lines in PDRs. Methods: We propose a statistical model of ortho-para conversion on dust grains with fluctuating dust temperatures, based on a master equation approach. This computation is then coupled to full PDR models and compared to PDR observations. Results: We show that the observations of rotational H2 lines indicate a high conversion efficiency on dust grains, and that this high efficiency can be accounted for if taking dust temperature fluctuations into account with our statistical model of surface conversion. Simpler models neglecting the dust temperature fluctuations do not reach the high efficiency deduced from the observations. Moreover, this high efficiency induced by dust temperature fluctuations is quite insensitive to the values of microphysical parameters of the model. Conclusions: Ortho-para conversion on grains is thus an efficient mechanism in most astrophysical conditions that can play a significant role in controlling the ortho-para ratio.

The H3+ molecule has been detected in many lines of sight within the central molecular zone (CMZ) with exceptionally large column densities and unusual excitation properties compared to diffuse local clouds. The detection of the (3,3) metastable level has been suggested to be the signature of warm and diffuse gas in the CMZ. We use the Meudon PDR code to re-examine the relationship between the column density of H3+ and the cosmic-ray ionization rate, $\zeta$, up to large values of $\zeta$. We study the impact of the various mechanisms that can excite H3+ in its metastable state. We produce grids of PDR models exploring different parameters ($\zeta$, size of clouds, metallicity) and infer the physical conditions that best match the observations toward ten lines of sight in the CMZ. For one of them, Herschel observations of HF, OH+, H2O+, and H3O+ can be used as additional constraints. We check that the results found for H3+ also account for the observations of these molecules. We find that the linear relationship between N(H3+) and $\zeta$ only holds up to a certain value of the cosmic-ray ionization rate, which depends on the proton density. A value $\zeta \sim 1 - 11 \times 10^{-14}$ s$^{-1}$ explains both the large observed H3+ column density and its excitation in the metastable level (3,3) in the CMZ. It also reproduces N(OH+), N(H2O+) and N(H3O+) detected toward Sgr B2(N). We confirm that the CMZ probed by H3+ is diffuse, nH $\lesssim$ 100 cm-3 and warm, T $\sim$ 212-505 K. This warm medium is due to cosmic-ray heating. We also find that the diffuse component probed by H3+ must fill a large fraction of the CMZ. Finally, we suggest the warm gas in the CMZ enables efficient H2 formation via chemisorption sites as in PDRs. This contributes to enhance the abundance of H3+ in this high cosmic-ray flux environment.

The $\mathrm{H}_2$ formation on grains is known to be sensitive to dust temperature, which is also known to fluctuate for small grain sizes due to photon absorption. We aim at exploring the consequences of simultaneous fluctuations of the dust temperature and the adsorbed H-atom population on the $\mathrm{H}_2$ formation rate under the full range of astrophysically relevant UV intensities and gas conditions. The master equation approach is generalized to coupled fluctuations in both the grain's temperature and its surface population and solved numerically. The resolution can be simplified in the case of the Eley-Rideal mechanism, allowing a fast computation. For the Langmuir-Hinshelwood mechanism, it remains computationally expensive, and accurate approximations are constructed. We find the Langmuir-Hinshelwood mechanism to become an efficient formation mechanism in unshielded photon dominated region (PDR) edge conditions when taking those fluctuations into account, despite hot average dust temperatures. It reaches an importance comparable to the Eley-Rideal mechanism. However, we show that a simpler rate equation treatment gives qualitatively correct observable results in full cloud simulations under most astrophysically relevant conditions. Typical differences are a factor of 2-3 on the intensities of the $\mathrm{H}_2$ $v=0$ lines. We also find that rare fluctuations in cloud cores are sufficient to significantly reduce the formation efficiency. Our detailed analysis confirms that the usual approximations used in numerical models are adequate when interpreting observations, but a more sophisticated statistical analysis is required if one is interested in the details of surface processes.