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.

The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms.

Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.

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.

Aims: In this paper, we aim to constrain the properties of dust structures in the central first parsecs of Active Galactic Nuclei (AGN). Our goal is to study the required optical depth and composition of different dusty and ionised structures. Methods: We developed a radiative transfer code, MontAGN, optimised for polarimetric observations in the infrared. With both this code and STOKES, designed to be relevant from the hard X-ray band to near-infrared wavelengths, we investigate the polarisation emerging from a characteristic model of the AGN environment. For that purpose, we compare predictions of our models with previous infrared observations of NGC 1068, and try to reproduce several key polarisation patterns revealed by polarisation mapping. Results: We constrain the required dust structures as well as their densities. More precisely, we find out that the electron density inside the ionisation cone is about $2.0 \times 10^9$ m$^{-3}$ . With structures constituted of spherical grains of constant density, we also highlight that the torus should be thicker than 20 in term of K band optical depth to block direct light from the centre. It should also have a stratification in density, with a lesser dense outer rim with an optical depth at 2.2 {\mu}m typically between 0.8 and 4 for observing the double scattering effect previously proposed. Conclusions: We bring constraints on the dust structures in the inner parsecs of an AGN model supposed to describe NGC 1068. When compared to observations, this leads to optical depth of at least 20 in Ks band for the torus of NGC 1068, corresponding to $\tau_V \approx 170$, which is within the range of current estimation based on observations. In the future, we will improve our study by including non uniform dust structures and aligned elongated grains to constrain other possible interpretations of the observations.

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.

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.