# Gas phase Elemental abundances in Molecular cloudS (GEMS) V. Methanol in Taurus

S. Spezzano, A. Fuente, P. Caselli, A. Vasyunin, D. Navarro-Almaid, M.

Rodríguez-Baras, A. Punanova, C. Vastel, V. Wakelam

► To cite this version:

S. Spezzano, A. Fuente, P. Caselli, A. Vasyunin, D. Navarro-Almaid, et al., Gas phase Elemental abundances in Molecular cloudS (GEMS) V. Methanol in Taurus. Astronomy and Astrophysics - A&A, 2022, 657, 10.1051/0004-6361/202141971. insu-03672069

HAL Id: insu-03672069

https://insu.hal.science/insu-03672069

Submitted on 19 May 2022

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# Gas phase Elemental abundances in Molecular cloudS (GEMS)
## V. Methanol in Taurus*

S. Spezzano<fnref n="1" />, A. Fuente<fnref n="2" />, P. Caselli<fnref n="1" />, A. Vasyunin<fnref n="3" />, D. Navarro-Almaida<fnref n="2" />, M. Rodríguez-Baras<fnref n="2" />, A. Punanova<fnref n="3" />, C. Vastel<fnref n="4" />, and V. Wakelam<fnref n="5" />

Received 5 August 2021 / Accepted 4 October 2021

## ABSTRACT

**Context.** Methanol, one of the simplest complex organic molecules in the interstellar medium, has been shown to be present and extended in cold environments such as starless cores. Studying the physical conditions at which CH₃OH starts its efficient formation is important to understand the development of molecular complexity in star-forming regions.

**Aims.** We aim to study methanol emission across several starless cores and investigate the physical conditions at which methanol starts to be efficiently formed, as well as how the physical structure of the cores and their surrounding environment affect its distribution.

**Methods.** Methanol and C¹⁸O emission lines at 3 mm have been observed with the IRAM 30 m telescope within the large programme Gas phase Elemental abundances in Molecular CloudS towards 66 positions across 12 starless cores in the Taurus Molecular Cloud. A non-LTE (local thermodynamic equilibrium) radiative transfer code was used to compute the column densities in all positions. We then used state-of-the-art chemical models to reproduce our observations.

**Results.** We have computed N(CH₃OH)/N(C¹⁸O) column density ratios for all the observed offsets, and the following two different behaviours can be recognised: the cores where the ratio peaks at the dust peak and the cores where the ratio peaks with a slight offset with respect to the dust peak (~10 000 AU). We suggest that the cause of this behaviour is the irradiation on the cores due to protostars nearby which accelerate energetic particles along their outflows. The chemical models, which do not take irradiation variations into account, can reproduce the overall observed column density of methanol fairly well, but they cannot reproduce the two different radial profiles observed.

**Conclusions.** We confirm the substantial effect of the environment on the distribution of methanol in starless cores. We suggest that the clumpy medium generated by protostellar outflows might cause a more efficient penetration of the interstellar radiation field in the molecular cloud and have an impact on the distribution of methanol in starless cores. Additional experimental and theoretical work is needed to reproduce the distribution of methanol across starless cores.

**Key words.** ISM: clouds – ISM: molecules – radio lines: ISM

## 1. Introduction

Methanol, CH₃OH, and even more complex organic molecules (COMs) have been widely observed in star-forming regions. COMs are defined as organic molecules with ≥6 atoms (van Dishoeck 2009), and their formation is considered to be the first step in the chemical complexity that will eventually be inherited by forming planets in the process of star and planetary system formation. In the past decade several papers have reported on the detection of methanol and other COMs towards starless cores showing that COMs are present and their emission is extended in cold and shielded environments (Bacmann et al. 2012; Bizzocchi et al. 2014; Vastel et al. 2014; Jiménez-Serra et al. 2016; Scibelli & Shirley 2020; Scibelli et al. 2021). Starless cores provide the ingredients that will eventually build-up a planet, both the organic as well as the refractory material. Therefore, it is valuable to understand the chemistry that drives

the build-up of chemical complexity in stellar nurseries such as starless cores.

COMs have been widely observed towards hot cores and hot corinos, which are both characterised by an active chemistry driven by ice sublimation at T ≥ 100 K (e.g. Caselli et al. 1993; Cazaux et al. 2003; Caux et al. 2011; Jørgensen et al. 2012; Belloche et al. 2014). Complex organic molecules are mostly formed on the surface of a dust grain via the followig: (i) hydrogenation of carbon monoxide, producing small COMs such as formaldehyde and methanol (Charnley et al. 1995; Watanabe & Kouchi 2002), and (ii) radical-radical (associative) reactions, producing larger COMs such as methyl formate and dimethyl ether (Garrod & Herbst 2006). In order for radical-radical COM formation to happen, it is crucial to have a central heating source such as a young stellar object (YSO) because the mobility of radicals on the surface appears at T ≥ 30 K. Because of the lack of a central heating source, the detection of COMs towards starless cores challenged our understanding of COMs in the interstellar medium. Large COMs such as acetaldehyde, dimethyl ether, and

<footnote n="1">¹ Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstrasse 1, 85748 Garching, Germany
e-mail: spezzano@mpe.mpg.de</footnote>

<footnote n="2">² Observatorio Astronómico Nacional (OAN), 28014 Madrid, Spain</footnote>

<footnote n="3">³ Ural Federal University, 620002, Mira st. 19, Yekaterinburg, Russia</footnote>

<footnote n="4">⁴ IRAP, Université de Toulouse, CNRS, UPS, CNES, 31400 Toulouse, France</footnote>

<footnote n="5">⁵ Laboratoire d'astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N, allée Geoffroy Saint-Hilaire, 33615 Pessac, France</footnote>

<footnote>* Based on observations carried out with the IRAM NOEMA interferometer. IRAM is supported by INSU/CNRS (France), MPG (Germany), and IGN (Spain).</footnote>
<table><caption>**Table 1.** Spectroscopic parameters of the observed lines.</caption><thead><tr><th>Molecule</th><th>Transition</th><th>Rest frequency<sup>(a)</sup><br>(MHz)</th><th><i>E</i><sub>up</sub><br>(K)</th><th><i>A</i><br>(×10<sup>-5</sup> s<sup>-1</sup>)</th><th><i>n</i><sup>*</sup><sup>(b)</sup><br>cm<sup>-3</sup></th></tr></thead><tbody><tr><td>CH<sub>3</sub>OH</td><td><i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 1<sub>0,1</sub>-0<sub>0,0</sub> (A<sup>+</sup>)</td><td>48 372.460(2)</td><td>2.3<sup>(c)</sup></td><td>3.5 × 10<sup>-7</sup></td><td>7 × 10<sup>3</sup></td></tr><tr><td>CH<sub>3</sub>OH</td><td><i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>)</td><td>96 739.358(2)</td><td>12.53<sup>(c)</sup></td><td>2.6 × 10<sup>-4</sup></td><td>3 × 10<sup>4</sup></td></tr><tr><td>CH<sub>3</sub>OH</td><td><i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 2<sub>0,2</sub>-1<sub>0,1</sub> (A<sup>+</sup>)</td><td>96 741.375(2)</td><td>6.96<sup>(c)</sup></td><td>3.4 × 10<sup>-4</sup></td><td>3 × 10<sup>4</sup></td></tr><tr><td>CH<sub>3</sub>OH</td><td><i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 2<sub>0,2</sub>-1<sub>0,1</sub> (E<sub>1</sub>)</td><td>96 744.545(2)</td><td>20.08<sup>(c)</sup></td><td>3.4 × 10<sup>-4</sup></td><td>3 × 10<sup>4</sup></td></tr><tr><td>CH<sub>3</sub>OH</td><td><i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 2<sub>0,2</sub>-1<sub>0,1</sub> (E<sub>1</sub>-E<sub>2</sub>)</td><td>108 893.945(2)</td><td>13.10<sup>(c)</sup></td><td>1.47 × 10<sup>-5</sup></td><td>5 × 10<sup>5</sup> <sup>(d)</sup></td></tr><tr><td>C<sup>18</sup>O</td><td><i>J</i> = 1–0</td><td>109 782.173(2)</td><td>5.27</td><td>6.3 × 10<sup>-8</sup></td><td>2 × 10<sup>3</sup></td></tr><tr><td>C<sup>18</sup>O</td><td><i>J</i> = 2–1</td><td>219 560.354(2)</td><td>15.81</td><td>6.0 × 10<sup>-7</sup></td><td>2 × 10<sup>4</sup></td></tr></tbody></table>

**Notes.** <sup>(a)</sup>Frequencies and uncertainties from the CDMS (Müller et al. 2005). <sup>(b)</sup><i>n</i>* is the critical density, calculated at 10 K. <sup>(c)</sup>Energy relative to the ground 0<sub>0,0</sub>, A rotational state. <sup>(d)</sup>Calculated with the rate coefficient reported in Bizzocchi et al. (2014).

methyl formate have been observed towards dark cloud cores and starless cores (Bacmann et al. 2012; Vastel et al. 2014; Jiménez-Serra et al. 2016; Nagy et al. 2019; Scibelli & Shirley 2020; Scibelli et al. 2021; Jiménez-Serra et al. 2021). The high-level of molecular complexity revealed by these studies towards cold cores has been a challenge for chemical models (Öberg et al. 2010). Recent laboratory work has shown that COMs such as acetaldehyde, vinyl alcohol, and also glyicine is possible in interstellar ice analogues with non-energetic processes (Chuang et al. 2020; Ioppolo et al. 2021). The formation of COMs on the grains depends on the efficiency for the large radicals to diffuse on the surfaces. Although it is an uncertain parameter, a larger diffusion efficiency will produce more COMs on the surfaces (Ruaud et al. 2015; Walsh et al. 2016). The formation in the ices also depends on the chemical network itself, which is probably not complete for the larger molecules. Ruaud et al. (2015) for instance have proposed new routes to form COMs at low temperature (and low diffusion efficiency) via reactions induced by the formation of complexes with the main ice components. New gas-phase pathways have also been proposed but still require the efficient production and desorption of methanol (Vasyunin & Herbst 2013; Balucani et al. 2015). If the efficiency on the surface is sufficient, remains the problem of desorbing these molecules strongly bound to the surfaces. The most promising non-thermal desorption mechanisms seem to be chemical reactive desorption from CO-rich surfaces (Minissale et al. 2016; Vasyunin et al. 2017), chemical explosions (Rawlings et al. 2013; Ivlev et al. 2015), and grain sputtering induced by cosmic-rays (Wakelam et al. 2021).

Methanol is one of the simplest COMs and the precursor of many of the O-bearing COMs, hence understanding the processes responsible for its release in the gas phase is of pivotal importance to constrain chemical models and reproduce the chemical complexity that we observe in starless cores.

In this paper we present our results on the emission of methanol towards 12 starless cores as well as towards the clouds where they are embedded. We study the dependency, or the lack thereof, between the emission of methanol and physical parameters such as visual extinction, volume density, as well as environment. The paper is structured as follows: Sect. 2 describes the observations, Sect. 3 presents the dataset as well as the TMC-1 molecular cloud and the B213 filament, where the 12 observed starless cores are located. The analysis of the observed data and its results are in Sects. 4 and 5. The results of the chemical modelling are presented in Sect. 6.2. Our conclusions are discussed in Sect. 7. More detailed information on each source is presented in the Appendix.

## 2. Observations

The molecular transitions used in this study are reported in Table 1. The 3 mm methanol lines were observed within the Gas phase Elemental abundances in Molecular CloudS (GEMS) (PI: Asunción Fuente) IRAM 30 m Large Program and the observing procedure and receiver setups were described by Fuente et al. (2019). The HPBW is varying with the frequency as HPBW'' = 2460 /ν where ν is in GHz. The observing mode was frequency switching with a frequency throw of 6 MHz well adapted to remove standing waves between the secondary mirror and the receivers. The Eight MIXer Receivers (EMIR) and the Fast Fourier Transform Spectrometers (FTS) with a spectral resolution of 49 kHz were used for these observations. The intensity scale is T<sub>MB</sub>, which is related with T<sub>A</sub><sup>*</sup> by T<sub>MB</sub> = (F<sub>eff</sub>/B<sub>eff</sub>)T<sub>A</sub><sup>*</sup> (see Table B.1 in Fuente et al. 2019). GEMS observations in TMC 1 were completed with C<sup>18</sup>O 2→1 observations carried out with the IRAM 30 m telescope in a previous project.

For TMC 1, we use observations of the CH<sub>3</sub>OH <i>J</i><sub>K<sub>a</sub>,K<sub>c</sub></sub> = 1<sub>0,1</sub> → 0<sub>0,0</sub> line carried out with the Yebes 40 m radiotelescope in Q-band (Tercero et al. 2020) during March-April 2018. The 40 m telescope is equipped with HEMT receivers for the 2.2–50 GHz range, and a SIS receiver for the 85–116 GHz range. Single-dish observations in K-band (21–25 GHz) and Q-band (41–50 GHz) were performed simultaneously. The backends consisted of FFTS covering a bandwidth of ~2 GHz in band K and ~9 GHz in band Q, with a spectral resolution of ~38 kHz. Central frequencies were 23 000 MHz and 44 750 MHz for the K and Q band receivers, respectively. The observing procedure was position-switching, and the OFF-positions are RA(J2000) = 04h42m24s.24 Dec(2000):25°41′27″.6 for TMC 1-CP, RA(J2000) = 04h42m29s.52 Dec(2000):25°48′07″.2 for TMC 1-NH<sub>3</sub>, RA(J2000) = 04h42m32s.16 Dec(J2000):25°59′42″.0 for TMC 1-C. These positions were checked to be empty of emission before the observations. The intensity scale is T<sub>MB</sub> with conversion factors of 4.1 Jy/K in band K (T<sub>MB</sub>/T<sub>A</sub><sup>*</sup> = 1.3) and in 5.7 Jy/K in band Q (T<sub>MB</sub>/T<sub>A</sub><sup>*</sup> = 2.1). The HPBW of the telescope is 42″ at 7 mm and 84″ at 1.3 cm. As an example of our dataset, Fig. 1 shows to the left the H<sub>2</sub> column density map of the TMC-1 molecular cloud (upper panel) and the B213 filament (lower panel) computed from Herschel and Planck data (Palmeirim et al. 2013; Rodríguez-Baras et al. 2021), and on the right a zoom-in into the H<sub>2</sub> column density map of B213-C1 where the observed offsets are marked (upper panel) and the spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol are overlaid with the 1–0 transition of C<sup>18</sup>O (lower panel). The
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Left panels: H₂ column density map of the TMC-1 molecular cloud (lower panel) and the B213 filament (upper panel) in Taurus, computed from Herschel and Planck data (Palmeirim et al. 2013; Rodríguez-Baras et al. 2021). The offsets observed towards the starless cores in our dataset are marked with a triangle. Right panel: H₂ column density of B213-C1 (upper panel). The triangles mark the positions observed with the GEMS large project, and the red circles show the beam size of the IRAM 30 m telescope at 97 GHz (~26"). The spectra of the 2<sub>1,2</sub>−1<sub>1,1</sub> (*E*₂) transition of methanol (in red) overlaid with the *J* = 1–0 transition of C<sup>18</sup>O (in black) in all offsets towards B213-C1 are shown in the lower panel.</figcaption></figure>

spectra observed towards the other sources in our sample are
shown in Figs. A.1–A.11.

3. Dataset

We have observed several lines of CH₃OH and two lines of C¹⁸O towards the GEMS sources in Taurus with the aim of understanding how these two chemically related molecules behave towards several starless cores within the same environment. The GEMS dataset allows us to study the abundance of both molecules, as well as their ratios, towards different offsets within and outside the cores, and hence tracing gas with different physical properties such as volume density and temperature. Our source sample is listed in Table 2.

All transitions of methanol studied in this paper, except the
J<sub>K<sub>a</sub>,K<sub>c</sub></sub> = 1<sub>0,1</sub>-0<sub>0,0</sub> (A<sup>+</sup>) transition at 48 GHz, have been observed
towards every offset in our source sample (although they were
not always detected). Both C<sup>18</sup>O transitions have been observed
towards every offsets in the starless cores in TMC-1, while only
the J = 1-0 transition has been observed towards the starless
cores in the B213 filament. Figures A.1-A.11 show the spectra
of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the 1-
0 transition of C<sup>18</sup>O towards all offsets in our dataset, with the
exception of B213-C1, shown in Fig. 1.

3.1. TMC-1 cloud

Taurus Molecular Cloud-1 (TMC-1) is a cold and dense cloud
with a filamentary structure of about 5' × 40' in the centre of the

<table><caption>**Table 2.** Source sample.</caption>
  <thead>
    <tr>
      <th></th>
      <th>Right ascension<br>(J2000)</th>
      <th>Declination<br>(J2000)</th>
    </tr>
  </thead>
  <tbody>
    <tr>
      <td>TMC-1 C</td>
      <td>04:41:38.80</td>
      <td>+25:59:42.0</td>
    </tr>
    <tr>
      <td>TMC-1 CP</td>
      <td>04:41:41.90</td>
      <td>+25:41:27.1</td>
    </tr>
    <tr>
      <td>TMC-1 NH<sub>3</sub></td>
      <td>04:41:21.30</td>
      <td>+25:48:07.0</td>
    </tr>
    <tr>
      <td>B213-1</td>
      <td>04:17:41.80</td>
      <td>+28:08:47.0</td>
    </tr>
    <tr>
      <td>B213-2</td>
      <td>04:17:50.60</td>
      <td>+27:56:01.0</td>
    </tr>
    <tr>
      <td>B213-5</td>
      <td>04:18:03.80</td>
      <td>+28:23:06.0</td>
    </tr>
    <tr>
      <td>B213-6</td>
      <td>04:18:08.40</td>
      <td>+28:05:12.0</td>
    </tr>
    <tr>
      <td>B213-7</td>
      <td>04:18:11.50</td>
      <td>+27:35:15.0</td>
    </tr>
    <tr>
      <td>B213-10</td>
      <td>04:19:37.60</td>
      <td>+27:15:31.0</td>
    </tr>
    <tr>
      <td>B213-12</td>
      <td>04:19:51.70</td>
      <td>+27:11:33.0</td>
    </tr>
    <tr>
      <td>B213-16</td>
      <td>04:21:21.00</td>
      <td>+27:00:09.0</td>
    </tr>
    <tr>
      <td>B213-17</td>
      <td>04:27:54.00</td>
      <td>+26:17:50.0</td>
    </tr>
  </tbody>
</table>

Taurus molecular cloud at a distance of 140.2<sup>+1.3</sup><sub>-1.3</sub> pc (Galli et al. 2019). The chemistry of TMC-1 has been the object of many studies in the literature (e.g. Guelin et al. 1982; Pratap et al. 1997; Saito et al. 2002; Schnee et al. 2007; Suutarinen et al. 2011). The most peculiar feature of TMC-1 is the chemical differentiation between the southern part, the cyanopolyyne peak (TMC-1 CP), where carbon-chain molecules are more abundant, and the northern part, the ammonia peak (TMC-1 NH₃), where NH₃ and N₂H⁺ are more abundant. The origin of this differentiation has
been extensively discussed, and it can be explained by either
a small variation in the chemical timescale (~10^5 yr), and/or
by the different density (with the ammonia peak being denser)
(Hirahara et al. 1992; Suzuki et al. 1992). Also a variation in the
C/O ratio can reproduce an abundance gradient similar to that
observed (Pratap et al. 1997). We present here the observations
of methanol and C^18O towards the cyanopolyyne and ammo-
nia peaks (TMC-1 CP and NH3). We also present observations
towards the pre-stellar core TMC1-C (Schnee et al. 2007, 2010).
Several observational evidences suggest that this core is a later
evolutionary stage than TMC1 (CP), presenting higher deutera-
tion fractions, more similar to those in L1544 (Crapsi et al. 2005;
Navarro-Almaida et al. 2021). The spectra of the 21,2–11,1 (E2)
transition of methanol overlaid with the 1–0 transition of C^18O
towards TMC-1C, TMC-1 CP and TMC-1 NH3 are shown in
Figs. A.1–A.3.

3.2. B213 filament

B213 (L1495) is one of the longest and most prominent filaments in the Taurus molecular cloud, as it extends for over 80', and it is located towards the north-west of the cloud at a distance of 129.9$^{+0.4}_{-0.3}$ pc (Galli et al. 2019). The gas towards B213 presents a very complex velocity structure that has been interpreted as the result of colliding filaments (Duvert et al. 1986). B213 has been observed in CO isotopologues with large scale maps (Onishi et al. 1996; Goldsmith et al. 2008). The denser gas has been studied using ammonia, $H^{13}CO^{+}$, $N_2H^+$ and SO (Benson & Myers 1989; Hacar et al. 2013). Furthermore, some dense cores in B213 contain young stellar objects (Rebull et al. 2010). We present here the observations of methanol and C$^{18}$O towards nine starless cores in B213. The spectra of the 21,2–11,1 (E2) transition of methanol overlaid with the 1–0 transition of C$^{18}$O towards the starless cores we observed in B213 are shown in Figs. 1, A.4–A.11.

4. Analysis

We have observed the three rotational transitions belonging to
the 2–1 transitions of methanol at 96 GHz (symmetry A, E1 and
E2) as well as the 1–0 transition at 108 GHz (symmetry E2–E1),
in all cores of our sample. In addition, we have observed the
$J_{K_o,K_c} = 1_{0,1}-0_{0,0}$ ($A^+$) transition of methanol at 48 GHz towards
the three cores in TMC-1. The intensity ratios among the four
lines at 3 mm in our dataset cannot be reproduced assuming LTE,
implying that we have to use a non-LTE code if we want to derive
a reliable column density of methanol in all positions within
our dataset. Since methanol lines have a relatively high critical
density (between 0.7 and $3 \times 10^4$ cm$^{-3}$, see Table 1), the cor-
respondent column density is very sensitive to variations in the
H$_2$ volume density. The H$_2$ volume density towards all observed
positions within the GEMS catalogue has been calculated by
modelling the emission of CS, C$^{34}$S, and $^{13}$CS with the radiative
transfer code RADEX (van der Tak et al. 2007) in Rodríguez-
Baras et al. (2021). However, the lines of methanol and CS could
come from different velocity components and/or different layers
within the same cloud. The line profiles of C$^{18}$O and CH$_3$OH
are different in many positions. Therefore, in order to be con-
sistent and derive accurate column densities we need to use the
volume density derived from CH$_3$OH data. We hence decided to
use the Markov chain Monte Carlo method (MCMC) together
with the RADEX non-LTE radiative transfer code within the
CASSIS software (Vastel et al. 2015) in order to compute the

methanol column density, the excitation temperature, and the gas
volume density towards all positions in our dataset. To assess
the goodness of the fit, the reduced χ² was computed for each
model. The reduced χ² in our sample ranges from 1.5 to 4.5,
with most of the models having a reduced χ² ~ 2. The results
of the MCMC+RADEX analysis with CASSIS on methanol are
reported in Table B.1. We used the collisional rates for methanol
reported in Rabli & Flower (2010) at 10 K.

We did not use the MCMC method for the C¹⁸O because for most of the cores we observed only one transition (with the exception of the three cores in TMC-1), and the observed lines of C¹⁸O are not very sensitive to variations in the volume densities. Furthermore, C¹⁸O often has more velocity components than CH₃OH, and we are interested in deriving a column density for C¹⁸O from the velocity component that corresponds to the one observed with CH₃OH. We used two approaches to derive the column densities of C¹⁸O in our dataset. In the offsets with a good match between the rest velocity and line width of the lines of methanol and C¹⁸O, we fitted the lines of C¹⁸O and used RADEX to derive the column density, assuming the T<sub>kin</sub> and n<sub>H2</sub> derived for CH₃OH. In these cases, the errors derived on the v<sub>LSR</sub> and δν from the line fit have been propagated to derive the error on the column density of C¹⁸O, see for example the values for TMC-1C and B213-C6 reported in Table B.1. For the offsets with a bad match between the rest velocity and line width of the lines of methanol and C¹⁸O (see for example the spectra of TMC-1 NH₃ in Fig. A.3), we have used the GUI of CASSIS and performed a single-line fit with RADEX of the portion of the spectrum of C¹⁸O that is comparable with methanol in velocity and line-width. In order to do this we have fixed the v<sub>LSR</sub> and adjusted the line-width in order to be similar to the line-width of methanol, and at the same time model well the line shape of C¹⁸O. The error on the column densities derived in this manner, and reported in Table B.1, has been estimated to be 15%. While in the case of CH₃OH our analysis shows that the molecules are not in LTE (as the resulting T<sub>ex</sub> is lower than the kinetic temperature), this is not the case for the C¹⁸O 1–0 transition, where T<sub>ex</sub> ~ T<sub>kin</sub>.

The visual extinction tabulated in Table B.1 has been
calculated from the H₂ column density maps reported in
Rodríguez-Baras et al. (2021), using the formula $A_V$ [mag] =
N(H₂)[cm⁻²]/9.4 × 10²⁰[cm⁻² mag⁻¹] (Bohlin et al. 1978).

We do not report on the upper limits for the column density of methanol towards positions where the line was not detected because it is not possible to derive them in the same fashion as we did for the methanol column densities reported in Table B.1, and we would have to assume a value for the H₂ volume density. The uncertainty induced by this assumption would make the comparison of the upper limits on N(CH₃OH) with the N(CH₃OH) computed with the MCMC+RADEX method not meaningful.

5. Results

Figure 2 shows the column density ratio of methanol and C18O in all observed positions in our sample plotted against the visual extinction and the dust temperature derived from Herschel and Planck data (Palmeirim et al. 2013; Rodríguez-Baras et al. 2021), as well as against the volume density derived from the methanol MCMC+RADEX analysis, presented here in Sect. 4. Figure 2 shows no correlation among N(CH3OH)/N(C18O) and A<sub>V</sub>, T<sub>dust</sub>, or n<sub>H2</sub>. This also holds if we consider B213 and TMC-1 separately or if we consider solely the column density of methanol instead of the column density ratio with respect to C18O. In
<figure><img src="image_3.png" /><figcaption>**Fig. 2.** Methanol and C<sup>18</sup>O column density ratios computed in all positions within our sample plotted against the: visual extinctions derived from Herschel and Planck data (Palmeirim et al. 2013; Rodríguez-Baras et al. 2021) (upper panel); volume density derived from the MCMC+RADEX analysis of the methanol lines (central panel); T<sub>dust</sub> derived from Herschel data in Rodríguez-Baras et al. (2021) (lower panel).</figcaption></figure>

Fig. 2, the points belonging to B213 and TMC-1 are shown as full circles and as crosses, respectively. In the case of multiple velocity components reported for methanol and C<sup>18</sup>O in Table B.1, we assumed that most of the dust emission comes from the fibre with a larger methanol column density (supposedly the densest). As a consequence, only the velocity component with the

largest methanol column density is plotted against $A_V$ and $T_{dust}$ in Fig. 2. In particular, in the case of TMC-1 CP and NH<sub>3</sub>, we plot the velocity component at 5.7 km s<sup>-1</sup>. As we aim at studying methanol in different positions across different cores, comparing solely the column density of methanol would give us a biased picture because of the different densities of the positions involved. On the other hand, comparing methanol abundances with respect to molecular hydrogen is also not an optimal solution given the presence of different velocity components in several spectra, corresponding to different fibres on the line of sight. We decided to use the column density ratio of methanol and carbon monoxide because they are chemically related, as the depletion of carbon monoxide on the surface of dust grains is necessary to the formation of methanol. Using the rare isotopologue C<sup>18</sup>O we can limit, if not exclude, the risk of large optical depth.

A possible reason for the lack of definite trends in Fig. 2 might be that methanol emits from the outer layers of starless cores (e.g. in L1544 as shown in Bizzocchi et al. 2014; Vastel et al. 2014). Furthermore, the emission of methanol is not homogeneous around pre-stellar and starless cores, and its distribution has been shown to depend on the density structure around the core, and hence the illumination onto the core due to the interstellar radiation field (Spezzano et al. 2016, 2020). In Spezzano et al. (2016), for example, the emission maps of methanol and cyclopropenylidene towards the inner 2.5" × 2.5" of the pre-stellar core L1544 show that both molecules have asymmetric distribution around the core, with methanol peaking towards the north-east and cyclopropenylidene towards the south-west. It is important to note that the visual extinction at the methanol and the c-C<sub>3</sub>H<sub>2</sub> emission peak in L1544 is the same, within error-bars. This shows that the H<sub>2</sub> column density of the medium computed on the line-of-sight is not a good enough indicator to discriminate among the regions where the illumination is effective enough to keep more Carbon in its atomic form and hence available to form hydrocarbons like c-C<sub>3</sub>H<sub>2</sub>, and the regions where there is enough shielding to allow carbon to be locked in to carbon monoxide, and eventually form methanol. In addition to the density structure around the core, also larger scale environmental effects have shown to influence the distribution of methanol around starless and pre-stellar cores. In Spezzano et al. (2020), methanol was mapped towards four pre-stellar and two starless cores, and it was shown that the methanol peak was influenced by the presence of nearby massive stars in the two cores mapped in Ophiucus (OphD and HMM-1).

The GEMS dataset offers the opportunity to study the emission of methanol towards twelve starless core, a larger sample with respect to previous studies. Furthermore, the GEMS dataset is not limited to the cores, giving us information on the transition from core to cloud. In order to spot trends that might be related to the density structure surrounding the core, in Figs. A.12–A.23 we plot the column density ratios ordered by offset position across the observed cut within the core for each source, while the offsets are shown on top of the N(H<sub>2</sub>) column density maps of the cores in the upper panel of the Figure to facilitate the comparison. The individual results are discussed in detail in Appendix A. Generally, we can recognise two different behaviours, shown in Fig. 3: the cores where N(CH<sub>3</sub>OH)/N(C<sup>18</sup>O) peaks at the dust peak, and the cores where the N(CH<sub>3</sub>OH)/N(C<sup>18</sup>O) ratios peak has a slight offset with respect to the dust peak (~10 000 AU). This trend is particularly clear in the cores where more than 5 offsets have been observed. B213-C1 and C5 for example show a clear peak of the column density ratio towards the dust peak
<figure><img src="image_4.png" /><figcaption>**Fig. 3.** Upper panel: methanol and C¹⁸O column density ratio for two cores in the north of B213, C1 and C5. Lower panel: methanol and C¹⁸O column density ratio for two cores in the south of B213, C10 and C16. The column density ratios are plotted with respect to the observed offsets, and offset 1 always refers to the dust peak (also shown as a dashed green line). The positions of the offsets in the cores are shown in Figs. A.15, A.17, A.20, and A.22.</figcaption></figure>

(see Figs. 3, A.15, and A.17), while in B213-C10 and C16 the peak is slightly shifted from the dust peak (~10 000 au), (see Figs. 3, A.20 and A.22). One reason for this behaviour might be the external irradiation on the cores. B213-C1 and C5 are in fact located towards the northern part of B213, that is known to be actively forming stars with ~40 Class 0/I/II sources present in the vicinity of B213-C1, C5, C2, C6 and C7. B213-C10 and C16, on the other side, are located in the central and southern part of B213, which only contains ~15 Class 0/I/II sources around B213-C10, C12 and C16, and 5 around B213-C17. In Fig. B.2 the positions of the young stellar objects in B213 are marked as red triangles on the H₂ column density map of the filament. Low-mass protostars can accelerate energetic particles along the outflow, and consequently increase the local cosmic rays flux (Padovani et al. 2016; Fitz Axen et al. 2021). This will in turn increase the destruction of CO, and consequently decrease the production of methanol towards the outer layers of starless cores in the north of B213.

The line-width of the methanol lines in the cores located in the north of B213 are ~25% broader than the lines of methanol observed towards the central and southern cores of B213, hinting at more prominent effects due to stellar activity in the northern part for B213. This effect has already been observed with

ammonia in Seo et al. (2015). An additional indication of the environmental effects is the fact that the dust temperature in the southern cores of B213 is lower with respect to the northern cores, at comparable visual extinctions, see Fig. B.1. In addition to the irradiation, there are also differences in the velocity structure. In B213-C1, C2, and C7, the profiles of C¹⁸O and CH₃OH are very different. In B213-C6, C10 and C16 we have however a better agreement between the profiles of C¹⁸O and CH₃OH.

The N(CH₃OH)/N(C¹⁸O) column density ratios towards the three cores in TMC-1 can be only calculated towards few offsets close to the dust peak, because methanol was not detected at larger radii, in contrast with B213, likely because the cores in B213 are denser. Towards TMC-1 NH₃ the two different velocities show peaks of the N(CH₃OH)/N(C¹⁸O) ratio at different offsets, maybe hinting at the fact that the gas they are tracing does not match in chemical age. Towards TMC-1 CP, if we exclude the offset 4 for the component at 5.6 km s⁻¹ because of the large error bar, both velocity components show the same trend, an increase of the N(CH₃OH)/N(C¹⁸O) ratio towards the dust peak. Towards TMC-1 C instead, the N(CH₃OH)/N(C¹⁸O) ratio is constant within error-bars.

# 6. Chemical modelling

## 6.1. Physical structure for spherical 1D static chemical model

For a proper chemical modelling of the region under study, it is necessary to have a good knowledge of the density and temperature across the region. However, sometimes our knowledge is limited to a few points, and geometrical considerations and simplifications are mandatory. Temperature parameterizations of prestellar cores can be found in, for example, Crapsi et al. (2007), in which they choose the spherically symmetric profile

$$T(r) = T_{\text{out}} - \frac{T_{\text{out}} - T_{\text{in}}}{1 + \left(\frac{r}{r_{\text{TO}}}\right)^{\alpha}}, \quad (1)$$

where Tout and Tin are the temperatures at the outermost and innermost observed positions, respectively. The dependence of the temperature with the radius is controlled by rTO, the flat radius of the temperature, and α, an asymptotic power index. A similar functional dependence with the radius is found in Plummer-like density profiles. These profiles assume spherical symmetry, and they have become popular density profiles to describe dense cores (Tafalla et al. 2002; Priestley et al. 2018):

$$n_{\text{H}}(r) = \frac{n_0}{1 + \left(\frac{r}{r_0}\right)^{\alpha}}, \quad (2)$$

where nH is the total atomic hydrogen number density across the profile, n0 is the central density, r0 is the flat radius, and α is the asymptotic power index. These profiles allow us to obtain density and temperature profiles across the B213 cuts in Table B.1. To have more data to derive physical profiles, we may fit data of different cuts with similar temperatures simultaneously. Therefore, we merge and average the data from the cut B213-C2 with B213-C6 (B213-C2-C6), B213-C10 with B213-C12 (B213-C10-C12), and consider B213-C16 as independent from the rest. The temperature profiles of B213-C2-C6, B213-C10-C12, and B213-C16 can be determined by finding the set of parameters that best fit the experimental data from Table B.1. The results are shown in Table 3.

We may perform a similar computation to derive the parameters of the density profiles. This time, however, we have information about the total hydrogen column density in the measured
<table><caption>**Table 3.** Parameter values for the temperature and density profiles presented in Sect. 6.1.</caption><thead><tr><th>Parameters</th><th>B213-C2-C6</th><th>B213-C10-C12</th><th>B213-C16</th></tr></thead><tbody><tr><td><i>T</i><sub>in</sub> (K)</td><td>10.94 ± 0.05</td><td>10.89 ± 0.04</td><td>10.25 ± 0.07</td></tr><tr><td><i>T</i><sub>out</sub> (K)</td><td>14.40 ± 0.05</td><td>14.78 ± 0.24</td><td>13.78 ± 0.41</td></tr><tr><td><i>r</i><sub>T0</sub> (au)</td><td>(7.03 ± 0.20) × 10<sup>3</sup></td><td>(1.46 ± 0.14) × 10<sup>4</sup></td><td>(1.04 ± 0.27) × 10<sup>4</sup></td></tr><tr><td>α</td><td>2.68 ± 0.16</td><td>1.56 ± 0.13</td><td>1.00 ± 0.16</td></tr><tr><td><i>n</i><sub>0</sub> (cm<sup>-3</sup>)</td><td>(3.95 ± 7.22) × 10<sup>5</sup></td><td>(2.06 ± 3.85) × 10<sup>5</sup></td><td>(1.27 ± 8.47) × 10<sup>6</sup></td></tr><tr><td><i>r</i><sub>0</sub> (au)</td><td>(2.60 ± 4.37) × 10<sup>3</sup></td><td>(5.00 ± 9.08) × 10<sup>3</sup></td><td>(4.20 ± 37.28) × 10<sup>2</sup></td></tr><tr><td>α</td><td>2.40 ± 1.58</td><td>2.52 ± 2.88</td><td>1.24 ± 1.39</td></tr></tbody></table>

visual extinction across B213-C2-C6, B213-C10-C12, and B213-C16 in Table B.1. Given the density profile, we can estimate the total hydrogen column density, assuming spherical symmetry, performing a projection along the line of sight:

$$N_H(r) = 2 \times \sum_i \Delta l_i \frac{n_H(s_i) + n_H(s_{i+1})}{2}, \quad (3)$$

where *r* is the impact parameter, $\Delta l_i = l_{i+1} - l_i$, $s_i = \sqrt{r^2 + l_i^2}$, $l_i$ is a discretisation of the segment along the line of sight $l_{max} > \dots > l_{i+1} > l_i > \dots > 0$, with $l_{max} = \sqrt{r_{max}^2 - r^2}$, and $r_{max}$ the radius of the density profile. The parameters that best describe the density of the B213 cuts would be those that best fit the visual extinction after the projection along the line of sight. The total hydrogen column density $N_H$ is compared to the visual extinction in Table B.1 using the relationship $A_v = N_H/1.88 \times 10^{21}$ (Bohlin et al. 1978), where it is assumed all hydrogen is in molecular form. To perform the fitting of the parameters to the extinction values, we have subtracted the background extinction, that is, the extinction at the outermost positions, to each cut. The results of the fitting process are shown in Table 3.

## 6.2. Chemical modelling approaches

Two approaches were used to reproduce our observations with chemical models. In the first, we used the physical structure of some of the B213 cores in our sample derived in Sect. 6.1, and ran a spherical 1D static chemical model to compute the radial abundance profiles of methanol and CO. In the second approach, applied to the low-density parts of the filament beyond the starless cores, we ran a 0D chemical model assuming a grid of temperature and $A_v$ to compute the abundances of methanol.

To model the emission of methanol in cold starless cores, it is essential to utilise a so-called three-phase astrochemical model. By three phases, we mean the gas phase, surface of icy mantles of interstellar grains, and bulk ice of grain mantles. The usage of a three-phase model, a more advanced approach in comparison to two-phase (gas phase + ice without distinction between surface and bulk), is justified by the importance of the details of chemical processes on interstellar grains for the formation of methanol and its abundance in the gas phase. Since Geppert et al. 2006, the consensus is that observed abundances of interstellar methanol cannot be formed in gas-phase chemical reactions. In contrast, formation of methanol during hydrogenation of CO molecules on interstellar grains is shown to be very efficient (e.g. Watanabe & Kouchi 2002). The delivery of methanol formed on cold (~10 K) grains to the gas phase is most likely due to the process of reactive desorption (Garrod et al. 2007; Minissale

et al. 2016). The efficiency of reactive desorption depends on the details of a particular chemically reacting system and the composition of the environment, that is the underlying ice surface (Minissale et al. 2016; Chuang et al. 2018). Thus, in order to model properly a complicated relation between CO and CH$_3$OH in prestellar clouds, we need to utilise the reasonably detailed approach to grain chemistry and gas-grain interaction. To model the emission of methanol, we utilised the MONACO model previously applied to studies of chemistry in a number of prestellar cores (Vasyunin & Herbst 2013; Vasyunin et al. 2017; Nagy et al. 2019; Lattanzi et al. 2020; Harju et al. 2020; Scibelli et al. 2021; Jiménez-Serra et al. 2021). MONACO is a rate equations-based numerical code capable of simulation chemistry in interstellar medium under a three-phase approach. As described in details in Vasyunin et al. (2017), the code includes treatment of chemical reactions in the gas phase, as well as diffusive chemistry on surfaces of icy mantles on interstellar grains, and in the bulk of icy mantles. The code includes treatment of thermal desorption, cosmic ray-induced desorption (Hasegawa & Herbst 1993), photodesorption with photodesorption yield equal to $10^{-3}$ particles per incident photon (Öberg et al. 2008) and reactive desorption following parametrisation by Minissale et al. (2016). The efficiency of reactive desorption is further adjusted to the fraction of ice surface covered with non-water species (Vasyunin et al. 2017). We assume low diffusion-to-desorption energy ratio $E_{diff}/E_{des} = 0.3$ for thermal diffusion of species on ice surface. In addition, quantum tunneling as a source of mobility for atomic and molecular hydrogen is assumed through a rectangular barrier for diffusion with thickness of 1.2 angstroms (Vasyunin et al. 2017). For diffusion of bulk species, $E_{diff}/E_{des}$ is assumed twice that of for surface, that is 0.6 following Garrod (2013). The results of the spherical 1D static chemical models are presented in Figs. 4 and B.3. More details on the physical profiles used for the cores can be found in Sect. 6.1.

To convert abundances of species per unit volume to directly observed values of column densities, we assumed that modelled starless cores are spherically symmetric with radial density profiles derived as described in Sect. 6.1. Convolution with IRAM 30 m beam has also been taken into account. The procedure is described in more details in Jiménez-Serra et al. (2016).

To explore the abundances of CO and CH$_3$OH in the molecular gas that surrounds starless cores in the filament, we utilised the same chemical model as for the cores. However in contrast to starless cores, the surrounding gas does not exhibit a well-shaped spatial structure. Thus, we utilised simple 0D chemical models with constant values of temperature and density. Since both parameters can vary in the surrounding gas, we calculated grid of 0D models on a parameter space that we believe reasonable covers the ranges of temperature and density in the filament gas. The results of such gridded modelling allows us to explore the behaviour of CO and CH$_3$OH abundances beyond the starless cores in B213 filament. The results of the 0D models are presented in Fig. 5.

## 6.3. Modelling results

Figure 4 compares the N(CH$_3$OH) derived from the spherical 1D static chemical model used with the physical structures derived for B213 C2-C6, C10-12 and C16 at $10^5$, $10^6$ and $10^7$ years, with the column densities derived in this paper for B213-C6, C10 and C16 and reported in Table B.1. The models at $10^5$ yr reproduces fairly well (within a factor of a few) the column densities of methanol. The model using the physical structure derived for
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** Methanol column densities calculated with the spherical 1D static chemical model using the physical structures derived for the cores B213-C2 and C6 (**left panel**), C10 (**central panel**), and C16 (**right panel**). The yellow, green and grey squares represent the results of the models extracted at $10^5$, $10^6$, and $10^7$ yr. The column densities observed towards B213-C6, C10 and C16 are shown as dark blue, red and light blue in the **left**, **central** and **right** panels, respectively.</figcaption></figure>

<figure><img src="image_6.png" /><figcaption>**Fig. 5.** Methanol column densities calculated with the grid 0D chemical model at $10^5$ yr (**left panel**), at $10^6$ yr (**central panel**), and at $10^7$ yr (**right panel**). The volume density used is a function of $A_V$ defined by the density profiles of C2-C6 cores used for the spherical 1D static chemical model.</figcaption></figure>

the cores C2 and C5 (**left panel** in Fig. 4) reproduces very well also the radial profile of observed column densities of methanol towards B213-C5. We decided to plot only the N(CH$_3$OH) computed for C5 and not for C2 because the latter has larger error bars. The radial profiles observed for C10 and C16 (**central and right panels** in Fig. 4) are instead not very well reproduced by the models. While it is important to keep in mind the intrinsic large uncertainties of the column densities computed by chemical models (Vasyunin et al. 2004), we observe that the changes in N(CH$_3$OH) as a function of radii cannot be reproduced by the chemical models for C10 and C16, while it can reproduced fairly well for C6. This is shown even if we look at the radial profiles of methanol and C$^{18}$O in Figs. A.12–A.23, with respect to the one predicted by the models in Fig. B.3. It is important to note that the physical structures used to run the spherical 1D static chemical model have been derived from Herschel observations, which are not sensitive to the central and densest part of the cores because of the large beam (~40′′). The volume density in starless cores could in fact have a steep increase towards the inner 40′′, especially in the more evolved pre-stellar cores (see for example Fig. 5 in Crapsi et al. 2007). Figure B.3 shows the radial distribution of the N(CH$_3$OH)/N(C$^{18}$O) column density ratio, as well as the column densities of CH$_3$OH and C$^{18}$O (calculated from the main isotopologue assuming a $^{16}$O/$^{18}$O ratio of 500, Solar value in Penzias 1981) in order to facilitate the comparison with our observational results shown in Figs. A.12–A.23. Independently from the physical structure used, the model predicts a drop of CH$_3$OH column density of about one order of magnitude

from $10^5$ to $10^6$ yr, and an additional, although less significant, drop from $10^6$ to $10^7$ yr. The column density of C$^{18}$O is underpredicted by the model. However, the model is only considering the core, and not also the cloud where the core is embedded. To check if the difference in C$^{18}$O column densities between model and observation is only due to the lack of the contribution from the cloud, we have calculated the C$^{18}$O column density at offsets outside of the cores, for example offset 7 in B213-C2, and then added it to the column densities predicted by the model, and we saw that the contribution of the cloud can account for the missing CO.

In some of our cores the observed cuts cover over 40 000 AU, hence merging into the cloud, while the physical models of the cores are derived assuming a Bonnor-Ebert sphere with a radius of 10 000 AU. We therefore computed a grid of models in the temperature range of 8–14 K and $A_V$ range of 2–8 mag. Please note in this case the $A_V$ is the local $A_V$, so it needs to be multiplied by a factor of two to be compared with the $A_V$ tabulated in Table B.1. The column density of methanol resulting from the grid models at $10^5$, $10^6$ and $10^7$ yr are shown in Fig. 5. Also in this case, the models cannot reproduce the variety of profiles that we observe. They predict a substantial drop of the methanol in gas-phase between $10^5$ and $10^6$ yr because of freeze-out, as well as a defined peak of the methanol column density, independent from the dust temperature, at $A_V \sim 4$ mag (corresponding to a total visual extinction of ~8 mag along the line of sight). This is in contrast with what we observe, that is the cores where the dust temperature is higher, tend to have the methanol peaking at
larger $A_V$. It is important however to note that protostellar feedback (e.g. higher cosmic-ray ionisation rate and/or shocks) is not included in the models.

## 7. Conclusions

Our observations of methanol with the GEMS dataset show that methanol is present and its emission is extended in cold starless cores and in their surrounding clouds. In our dataset in fact, methanol has been detected at visual extinction ranging from 4 to 27 mag. Furthermore, the different behaviours in our sample of starless cores in Taurus show that local environmental differences are significative and have an impact on the distribution of methanol.

Previous work on emission maps of methanol towards four starless and two pre-stellar cores suggested that the asymmetric distribution of methanol is linked to the large scale density structure around the cores as well as different amount of illumination from the external radiation field (Spezzano et al. 2016, 2020). With this work we can confirm that environmental effects have a strong impact on the methanol distribution. However, we show that the column density and temperature of the medium computed on the line-of-sight are not good enough indicators to identify the effects of the environment on the methanol, and that the level of star formation activity needs to be taken into account. It is clear from our data that low-mass protostars in the surrounding environment have an impact on the distribution of methanol towards starless cores. Low-mass protostars in fact can accelerate energetic particles and increase the local cosmic-ray flux and ionisation rate. This can locally increase the abundance of $He^+$, the main destruction partner of CO in dark clouds, thus reducing its abundance in the gas and solid form (the latter being the first step toward $CH_3OH$ production). Increased fluxes of cosmic rays can also increase the CO desorption rate (e.g. Hasegawa & Herbst 1993), again reducing the amount of available solid CO for its transformation into methanol via successive hydrogenation. The higher CO desorption rate could compensate the gas phase destruction of CO via $He^+$; indeed we do not see any local reduction of the $C^{18}O$ column density around the cores in the north. Dedicated observations targeting molecular ions are needed to map the variation of ionisation rate across low-mass star forming regions in order to test the effect on the chemistry with the chemical models. Another consequence of the presence of low-mass stars nearby is that they contribute to the clumpiness of the medium, escavated by the cones of the outflows, which could lead to a more efficient penetration of the interstellar radiation field in the molecular cloud. In these conditions, only the densest regions (the centres of the starless cores) will be screened enough from UV photons to allow CO to accumulate on the surface of dust grains. The higher dust temperature measured by Herschel in the northern region of B213 could indeed hint to this scenario, but higher resolution mapping of the dust continuum emission and gas tracers are needed to make quantitative conclusions.

The comparison with state-of-the-art chemical models shows that overall the agreement of the modelled column densities with the observed ones is quite good, and that reactive desorption is a very efficient method to release the methanol that formed on the surface of dust grains into the gas phase, reproducing the observations. The discrepancies among models and observations in the reproduction of the radial profiles might be due to incomplete formation, destruction, and desorption mechanisms in the chemical models, to the fact that the protostellar feedback is not included in the models, as well as to the intrinsic limit

of using spherical models to reproduce an asymmetric emission. Additional experimental and theoretical work on, for example, the chemical desorption and cosmic ray sputtering on different ice mixtures would be helpful for a better match between observations and models. Collisional rates for methanol with molecular hydrogen at temperatures lower than 10 K would also be beneficial.

**Acknowledgements.** The authors are grateful to the anonymous referee for insightful comments. A large part of the data analysis described in this paper was performed during the spring of 2020, in the beginning of the COVID pandemic and during a hard lockdown. S.S. wishes to thank the Max Planck Society for the flexibility that was allowed during the pandemic, because it contributed to maintaining a clear and focus mind during the hours that she could dedicate to her work, and overall to keep calm, while waiting for the 'storm' to pass. Based on analysis carried out with the CASSIS software (http://cassis.irap.omp.eu) and CDMS and JPL spectroscopic databases and LAMDA molecular databases. CASSIS has been developed by IRAP-UPS/CNRS. S.S. wishes to thank the Max Planck Society for the Independent Max Planck Research Group funding. A.F., D.N.A. and M.R.B. are funded by Spanish MICINN through PID2010-106235GB-I00 national research project. V.W. acknowledges the CNRS program Physique et Chimie du Milieu Interstellaire (PCMI) co-funded by the Centre National d'Etudes Spatiales (CNES). A.V. and A.P. are the members of the Max Planck Partner Group at the Ural Federal University. A.V. and A.P. acknowledge the support of the Russian Ministry of Science and Education via the State Assignment Contract no. FEUZ-2020-0038.

## References

André, P., Men'shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102

Bacmann, A., Taquet, V., Faure, A., et al. 2012, A&A, 541, L12

Balucani, N., Ceccarelli, C., & Taquet, V. 2015, MNRAS, 449, L16

Belloche, A., Müller, H. S. P., Menten, K. M., et al. 2014, A&A, 561, C1

Benson, P. J., & Myers, P. C. 1989, ApJS, 71, 89

Bernard, J.-P., Paradis, D., Marshall, D. J., et al. 2010, A&A, 518, L88

Bizzocchi, L., Caselli, P., Spezzano, S., et al. 2014, A&A, 569, A27

Bohlin, R. C., Savage, B. D., & Drake, J. F. 1978, ApJ, 224, 132

Caselli, P., Hasegawa, T. I., & Herbst, E. 1993, ApJ, 408, 548

Caselli, P., Walmsley, C. M., Zucconi, A., et al. 2002, ApJ, 565, 331

Caux, E., Kahane, C., Castets, A., et al. 2011, A&A, 532, A23

Cazaux, S., Tielens, A. G. G. M., Ceccarelli, C., et al. 2003, ApJ, 593, L51

Charnley, S. B., Kress, M. E., Tielens, A. G. G. M., et al. 1995, ApJ, 448, 232

Chuang, K.-J., Fedoseev, G., Qasim, D., et al. 2018, A&A, 617, A87

Chuang, K.-J., Fedoseev, G., Qasim, D., et al. 2020, A&A, 635, A199

Crapsi, A., Caselli, P., Walmsley, C. M., et al. 2005, ApJ, 619, 379

Crapsi, A., Caselli, P., Walmsley, M. C., et al. 2007, A&A, 470, 221

Dartois, E., Augé, B., Boduch, P., et al. 2015, A&A, 576, A125

Duvert, G., Cernicharo, J., & Baudry, A. 1986, A&A, 164, 349

Fitz Axen, M., Offner, S. S. S., Gaches, B. A. L., et al. 2021, ApJ, 915, 43

Fuente, A., Navarro, D. G., Caselli, P., et al. 2019, A&A, 624, A105

Galli, P. A. B., Loinard, L., Bouy, H., et al. 2019, A&A, 630, A137

Garrod, R. T. 2013, ApJ, 778, 158

Garrod, R. T., & Herbst, E. 2006, A&A, 457, 927

Garrod, R. T., Wakelam, V., & Herbst, E. 2007, A&A, 467, 1103

Geppert, W. D., Hamberg, M., Thomas, R. D., et al. 2006, Faraday Discuss., 133, 177

Goldsmith, P. F., Heyer, M., Narayanan, G., et al. 2008, ApJ, 680, 428

Guelin, M., Langer, W. D., & Wilson, R. W. 1982, A&A, 107, 107

Hacar, A., Tafalla, M., Kauffmann, J., et al. 2013, A&A, 554, A55

Harju, J., Pineda, J. E., Vasyunin, A. I., et al. 2020, ApJ, 895, 101

Hasegawa, T. I., & Herbst, E. 1993, MNRAS, 261, 83

Hirahara, Y., Suzuki, H., Yamamoto, S., et al. 1992, ApJ, 394, 539

Ioppolo, S., Fedoseev, G., Chuang, K.-J., et al. 2021, Nat. Astron., 5, 197

Ivlev, A. V., Padovani, M., Galli, D., et al. 2015, ApJ, 812, 135

Jiménez-Serra, I., Vasyunin, A. I., Caselli, P., et al. 2016, ApJ, 830, L6

Jiménez-Serra, I., Vasyunin, A. I., Spezzano, S., et al. 2021, ApJ, 917, 44

Jørgensen, J. K., Favre, C., Bisschop, S. E., et al. 2012, ApJ, 757, L4

Keto, E., & Caselli, P. 2010, MNRAS, 402, 1625

Lattanzi, V., Bizzocchi, L., Vasyunin, A. I., et al. 2020, A&A, 633, A118

Minissale, M., Moudens, A., Baouche, S., et al. 2016, MNRAS, 458, 2953

Müller, H. S. P., Schlöder, F., Stutzki, J., et al. 2005, J. Mol. Struc., 742, 215

Navarro-Almaida, D., Le Gal, R., Fuente, A., et al. 2020, A&A, 637, A39

Navarro-Almaida, D., Fuente, A., Majumdar, L., et al. 2021, A&A, 653, A15

Nagy, Z., Spezzano, S., Caselli, P., et al. 2019, A&A, 630, A136

Öberg, K. I., van Dishoeck, E. F., & Linnartz, H. 2008, Organic Matter in Space, 251, 449
Öberg, K. I., Bottinelli, S., Jørgensen, J. K., et al. 2010, ApJ, 716, 825

Onishi, T., Mizuno, A., Kawamura, A., et al. 1996, ApJ, 465, 815

Onishi, T., Mizuno, A., Kawamura, A., et al. 2002, ApJ, 575, 950

Padovani, M., Marcowith, A., Hennebelle, P., et al. 2016, A&A, 590, A8

Palmeirim, P., André, P., Kirk, J., et al. 2013, A&A, 550, A38

Penzias, A. A. 1981, ApJ, 249, 518

Pratap, P., Dickens, J. E., Snell, R. L., et al. 1997, ApJ, 486, 862

Priestley, F. D., Viti, S., & Williams, D. A. 2018, AJ, 156, 51

Rabli, D., & Flower, D. R. 2010, MNRAS, 406, 95

Rawlings, J. M. C., Williams, D. A., Viti, S., et al. 2013, MNRAS, 436, L59

Reboussin, L., Wakelam, V., Guilloteau, S., et al. 2014, MNRAS, 440, 3557

Rebull, L. M., Padgett, D. L., McCabe, C.-E., et al. 2010, ApJS, 186, 259

Rodríguez-Baras, M., Fuente, A., Riviére-Marichalar, P., et al. 2021, A&A, 648, A120

Ruaud, M., Loison, J. C., Hickson, K. M., et al. 2015, MNRAS, 447, 4004

Saito, S., Aikawa, Y., Herbst, E., et al. 2002, ApJ, 569, 836

Schnee, S., Caselli, P., Goodman, A., et al. 2007, ApJ, 671, 1839

Schnee, S., Enoch, M., Noriega-Crespo, A., et al. 2010, ApJ, 708, 127

Scibelli, S., & Shirley, Y. 2020, ApJ, 891, 73

Scibelli, S., Shirley, Y., Vasyunin, A., et al. 2021, MNRAS, 504, 5754

Seo, Y. M., Shirley, Y. L., Goldsmith, P., et al. 2015, ApJ, 805, 185

Spezzano, S., Bizzocchi, L., Caselli, P., et al. 2016, A&A, 592, L11

Spezzano, S., Caselli, P., Pineda, J. E., et al. 2020, A&A, 643, A60

Suutarinen, A., Geppert, W. D., Harju, J., et al. 2011, A&A, 531, A121

Suzuki, H., Yamamoto, S., Ohishi, M., et al. 1992, ApJ, 392, 551

Tafalla, M., Myers, P. C., Caselli, P., et al. 2002, ApJ, 569, 815

Tercero, B., Cernicharo, J., Cuadrado, S., et al. 2020, A&A, 636, L7

van der Tak, F. F. S., Black, J. H., Schöier, F. L., et al. 2007, A&A, 468, 627

van Dishoeck, E. F. 2009, Astrophys. Space Sci. Proc., 10, 187

Vastel, C., Ceccarelli, C., Lefloch, B., et al. 2014, ApJ, 795, L2

Vastel, C., Bottinelli, S., Caux, E., et al. 2015, SF2A-2015: Proceedings of the Annual meeting of the French Society of Astronomy and Astrophysics, 313

Vasyunin, A. I., & Herbst, E. 2013, ApJ, 769, 34

Vasyunin, A. I., Caselli, P., Dulieu, F., et al. 2017, ApJ, 842, 33

Wakelam, V., Dartois, E., Chabot, M., et al. 2021, A&A, 652, A63

Walsh, C., Loomis, R. A., Öberg, K. I., et al. 2016, ApJ, 823, L10

Ward-Thompson, D., Scott, P. F., Hills, R. E., et al. 1994, MNRAS, 268, 276

Watanabe, N., & Kouchi, A. 2002, ApJ, 571, L173
# Appendix A: Additional information on the single cores

## Appendix A.1: TMC-1 C

We have observed CH₃OH and C¹⁸O (both lines) towards four positions in TMC-1 C. The $J_{K_a,K_c} = 2_{1,2-1,1} (E_2)$, $2_{0,2-1,0,1} (A^+)$, $2_{0,2-1,0,1} (E_1)$, and $2_{0,2-1,0,1} (E_1-E_2)$ transitions of methanol have been observed towards all four positions, while the $1_{0,1-0,0,0} (A^+)$ transition has been observed only towards the offsets 3 and 4. All lines have been detected with a signal to noise larger than 3 towards all positions. The spectra of the $2_{1,2-1,1,1} (E_2)$ transition of methanol overlaid with the $J = 1-0$ transition of C¹⁸O are shown in Figure A.1, and they show a very good match both in $v_{LSR}$ and in line-width.

Figure A.12 shows in the upper panel the observed offsets within TMC-1 C on the H₂ column density map derived from *Herschel* and *Planck* data (Rodríguez-Baras et al. 2021) data, and the variation of the CH₃OH and C¹⁸O column density ratio, as well as the single column densities N(CH₃OH) and N(C¹⁸O) in the observed cut across TMC-1 C in the lower panel. Both N(CH₃OH) and N(C¹⁸O) increase when going towards dust peak of TMC-1 C. The column density ratio is constant, within error bars, across the observed cut.

## Appendix A.2: TMC-1 CP

We have observed CH₃OH and C¹⁸O (both lines) towards four positions in TMC-1 CP. The $J_{K_a,K_c} = 2_{1,2-1,1} (E_2)$, $2_{0,2-1,0,1} (A^+)$, $2_{0,2-1,0,1} (E_1)$, and $2_{0,2-1,0,1} (E_1-E_2)$ transitions of methanol have been observed towards all four positions, while the $1_{0,1-0,0,0} (A^+)$ transition has been observed only towards the offset 4. All lines have been detected with a signal to noise larger than 3 towards all positions. The spectra of the $2_{1,2-1,1,1} (E_2)$ transition of methanol overlaid with the $J = 1-0$ transition of C¹⁸O are shown in Figure A.2. TMC-1 CP is known to have multiple velocity components (Fuente 2019 and references therein). From the spectra shown in Figure A.2 we see that methanol is tracing two of the many velocity components seen in C¹⁸O.

Two velocity components are present towards TMC-1 CP, one at 5.6 km/s and one at 6.0 km/s, with the latter having column densities of both CH₃OH and C¹⁸O larger by about a factor of 2 with respect to the lower velocity component. Figure A.13 shows the observed offsets on the H₂ column density map in the upper panel, and the variation of the CH₃OH and C¹⁸O column density ratio, as well as the single column densities N(CH₃OH) and N(C¹⁸O) in the observed cut across the core for both velocity components, in the lower panel. In the lower velocity component the column density of both methanol and C¹⁸O increase towards the dust peak. Excluding the offset 4, that has a larger error in comparison with the other positions, also the column density ratio increases towards the dust peak.

## Appendix A.3: TMC-1 NH₃

We have observed CH₃OH and C¹⁸O (both lines) towards four positions in TMC-1 NH₃. The $J_{K_a,K_c} = 2_{1,2-1,1} (E_2)$, $2_{0,2-1,0,1} (A^+)$, $2_{0,2-1,0,1} (E_1)$, and $2_{0,2-1,0,1} (E_1-E_2)$ transitions of methanol have been observed towards all four positions, while the $1_{0,1-0,0,0} (A^+)$ transition has been observed only towards the offset 4. All lines have been detected with a signal to noise larger than 3 towards all positions. The spectra of the $2_{1,2-1,1,1} (E_2)$ transition of methanol overlaid with the $J = 1-0$ transition of C¹⁸O are shown in Figure A.3. The velocity structure that we can see in the

C¹⁸O lines is quite complex, and the mismatch with the peaks of the methanol lines suggests the presence of more velocity components that we can actually resolve, at least in C¹⁸O.

TMC-1 NH₃ is the ammonia peak in TMC-1 and supposed to be the most evolved core, with respect to TMC-1 CP for example (Hirahara et al. 1992). Two velocity components are present for methanol, while a more complex structure is present in C¹⁸O, see Figure A.3. Figure A.14 shows in the upper panel the observed offsets within TMC-1 NH₃ on the H₂ column density map derived from *Herschel* and *Planck* data (Rodríguez-Baras et al. 2021), and the variation of the CH₃OH and C¹⁸O column density ratio, as well as the single column densities N(CH₃OH) and N(C¹⁸O) in the observed cut across the core for both velocity components, in the lower panel. In the component at 5.7 km s⁻¹, the column densities of both molecules have a peak at the offset 3, and then decreases moving towards the dust peak. The column density ratio shows the same behaviour. In the component at 6.0 km/s instead, while the column density of C¹⁸O decreases towards the dust peak, the methanol column density has a peak at the offset 2, and so does the column density ratio.

## Appendix A.4: B213 C1

We have observed all lines of CH₃OH reported in Table 1, with the exception of the $1_{0,1-0,0,0} (A^+)$ transition, and the 1-0 transition of C¹⁸O towards nine positions in B213-C1. The observed lines have been detected with a signal to noise larger than 3 towards all positions. The spectra of the $2_{1,2-1,1,1} (E_2)$ transition of methanol overlaid with the $J = 1-0$ transition of C¹⁸O are shown in Figure 1. The velocity structure is quite complex, with both molecules showing multiple velocity components. A shift in the $v_{LSR}$ of the main velocity component across the observed cut, is also clearly seen.

The upper panel in Figure A.15 shows the H₂ column density map of B213-C1 as well as the offsets of our observations. The observed stripe is centred at the denser peak, labeled as offset 1, but also passes through a less dense tail towards the east (offsets 6 and 7). In the lower panel of Figure A.15 are shown both the column density ratio of CH₃OH/C¹⁸O and the single column densities across the observed stripe, which all decrease rather sharply towards the east after the offset 3 and peak towards the dust peak.

## Appendix A.5: B213 C2

We have observed all lines of CH₃OH reported in Table 1, with the exception of the $1_{0,1-0,0,0} (A^+)$ transition, and the 1-0 transition of C¹⁸O towards nine positions in B213-C2. The observed lines have been detected with a signal to noise larger than 3 towards positions 1, 2, 8 and 9. The spectra of the $2_{1,2-1,1,1} (E_2)$ transition of methanol overlaid with the $J = 1-0$ transition of C¹⁸O are shown in Figure A.4. The C¹⁸O transition shows two velocity components in all position with the exception of positions 6 and 7. The methanol transition instead shows a brighter velocity component at ~7 km/s.

The H₂ column density map of B213-C2 is shown in the upper panel of Figure A.16 with all the offsets observed within B213-C2 in the GEMS large project, however methanol was only detected towards the densest offsets (1,2,8 and 9). Both the variation of the CH₃OH and C¹⁸O column density ratio, as well as the single column densities N(CH₃OH) and N(C¹⁸O) are shown in the lower panel. The column densities of CH₃OH and C¹⁸O show different behaviours. While C¹⁸O is at its lowest at the dust peak and increases towards the west (offsets 8 and 9), methanol peaks
at the dust peak and decreases towards offsets 8 and 9, although the error bars on the methanol column density at the dust peak are rather large so it might be only a slight increase. The column density ratio also peaks towards the dust peak of B213-C2.

## Appendix A.6: B213 C5

We have observed all lines of CH<sub>3</sub>OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A<sup>+</sup>) transition, and the 1-0 transition of C<sup>18</sup>O towards nine positions in B213-C5. The observed lines have been detected with a signal to noise larger than 3 towards all positions. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O are shown in Figure A.5. The velocity structure is rather complex, with the methanol line tracing a portion of a lower velocity wing in the C<sup>18</sup>O line towards the positions 1, 2, 3, and 8. Two main velocity components are visible in both lines towards the positions 4, 5, and 9. While a single velocity component in both lines is present at the positions 6 and 7.

The upper panel in Figure A.17 shows the H<sub>2</sub> column density map of B213-C5 as well as the offsets of our observations. The observed stripe is centred at a dense peak towards the west of the map, labeled as offset 1, and also passes through an another dense clump towards the east (offsets 6 and 7). The column densities and column density ratio plots in the lower panel of Figure A.17 show a very clear behaviour: the C<sup>18</sup>O column density decreases when moving from east to the west in the observed offsets, and both the methanol column density and the column density ratio have a very well defined peak at the dust peak of B213-C5 (offset 1).

## Appendix A.7: B213 C6

We have observed all lines of CH<sub>3</sub>OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A<sup>+</sup>) transition, and the 1-0 transition of C<sup>18</sup>O towards nine positions in B213-C6. The observed lines have been detected with a signal to noise larger than 3 towards positions 1, 2, 3, 8, and 9. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O are shown in Figure A.6. In this case, the velocity structure is quite simple compared to the others cores in B213, with only one prominent velocity component. Both v<sub>LSR</sub> and line widths match very well for C<sup>18</sup>O and CH<sub>3</sub>OH.

The H<sub>2</sub> column density map of B213-C6 is shown in the upper panel of Figure A.18 with all the offsets observed within B213-C2 in the GEMS large project, however methanol was only detected towards the densest offsets (1, 2, 3, 8 and 9). The column densities and column density ratio plots are shown in the lower panel: the C<sup>18</sup>O column density decreases towards the denser offsets, while both the methanol column density and the column density ratio increase towards the dust peak of B213-C6.

## Appendix A.8: B213 C7

We have observed all lines of CH<sub>3</sub>OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A<sup>+</sup>) transition, and the 1-0 transition of C<sup>18</sup>O towards nine positions in B213-C7. The observed lines have been detected with a signal to noise larger than 3 towards positions 1, 2, 3, 4, 5, 8 and 9. Only C<sup>18</sup>O has been observed towards the position 6, and neither C<sup>18</sup>O nor CH<sub>3</sub>OH have been observed towards position 7. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O are shown in Figure A.7. Different velocity components are present in the different positions. Two

velocity components (at ~6 and 7 km/s) are visible in positions 1, 8 and 9 for C<sup>18</sup>O, with only the brightest being observed also in methanol. Only the component at 7 km/s is visible in both molecules towards the positions 2 and 3. Positions 5 and 6 show only one velocity component at ~6 km/s. In positions 4 two velocity components are also present, but with a smaller separation in velocity with respect to positions 1, 8 and 9.

The H<sub>2</sub> column density map of B213-C7 is shown in the upper panel of Figure A.19 with all the offsets observed within B213-C2 in the GEMS large project, however methanol was only detected towards the densest offsets (1, 2, 3, 8 and 9). The column densities and column density ratio plots are shown in the lower panel of Figure A.19. The column density of C<sup>18</sup>O has a peak at the offset 2 and then decreases towards the west. The methanol column density, as well as the column density ratio have an oscillating pattern.

## Appendix A.9: B213 C10

We have observed all lines of CH<sub>3</sub>OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A<sup>+</sup>) transition, and the 1-0 transition of C<sup>18</sup>O towards nine positions in B213-C10. The observed lines have been detected with a signal to noise larger than 3 towards all positions except positions 5 and 6. Only C<sup>18</sup>O has been observed towards the position 5. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O are shown in Figure A.8. Two velocity components at ~5.5 and ~7 km/s are present for C<sup>18</sup>O towards the positions 1, 2, 3, 7, 8, and 9. Only the brightest velocity component is observed in methanol.

The upper panel in Figure A.20 shows the H<sub>2</sub> column density map of B213-C10 as well as the offsets of our observations. Towards B213-C10, methanol was observed in all offsets with the exception of offsets 5 and 6. The column densities and column density ratio plots are shown in the lower panel of Figure A.20. The C<sup>18</sup>O column density shows two peaks, one around offsets 2 and 3, and another at offset 8. The methanol column density, as well as the column density ratio instead have a clear peak only at offsets 2 and 3, towards a tail of in the H<sub>2</sub> column density map.

## Appendix A.10: B213 C12

We have observed all lines of CH<sub>3</sub>OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A<sup>+</sup>) transition, and the 1-0 transition of C<sup>18</sup>O towards nine positions in B213-C12. The observed lines have been detected with a signal to noise larger than 3 towards positions 1, 2, 3, 4, 8 and 9. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (E<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O are shown in Figure A.9. The C<sup>18</sup>O lines show a complex velocity structure, with multiple velocity components unresolved towards the positions 1, 2, 3, and 4, and two main velocity components towards positions 8 and 9. The methanol line shows only one velocity component at ~7 km/s.

The H<sub>2</sub> column density map of B213-C12 is shown in the upper panel of Figure A.21 with all the offsets observed within B213-C12 in the GEMS large project, however methanol was only detected towards the densest offsets (1, 2, 3 and 8). Both the variation of the CH<sub>3</sub>OH and C<sup>18</sup>O column density ratio, as well as the single column densities N(CH<sub>3</sub>OH) and N(C<sup>18</sup>O) are shown in the lower panel. The column densities of CH<sub>3</sub>OH and C<sup>18</sup>O show different behaviours: C<sup>18</sup>O increases from offsets 3 to offset 8, passing through the dust peak at offset 1, while methanol shows the opposite behaviour. The column density shows the same trend as the methanol column density.
## Appendix A.11: B213 C16

We have observed all lines of CH₃OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A⁺) transition, and the 1-0 transition of C¹⁸O towards nine positions in B213-C16. The observed lines have been detected with a signal to noise larger than 3 towards all positions except positions 5 and 6, where only the line of C¹⁸O have been observed. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E₂*) transition of methanol overlaid with the *J* = 1-0 transition of C¹⁸O are shown in Figure A.10. The C¹⁸O lines show wings in most positions, hinting at the presence of more than one velocity component on the line of sight. The methanol lines do not show wings, with the exception of the position 4.

The upper panel in Figure A.22 shows the H₂ column density map of B213-C16 as well as the offsets of our observations. Towards B213-C16, methanol was observed in all offsets with the exception of offsets 5 and 6. The column densities and column density ratio plots are shown in the lower panel of Figure A.22. The C¹⁸O column density shows a very shallow peak around offsets 2 and 3. Both the methanol column density, as well as the column density ratio show a steep increase moving from offset 2 towards offset 4 (opposite with respect to the direction of the dust peak), and have a shallower at offsets 1 and 7.

## Appendix A.12: B213 C17

We have observed all lines of CH₃OH reported in Table 1, with the exception of the 1<sub>0.1</sub>-0<sub>0.0</sub> (A⁺) transition, and the 1-0 transition of C¹⁸O towards nine positions in B213-C17. The observed lines have been detected with a signal to noise larger than 3 towards positions 1, 2, 3, 7, 8, and 9 for C¹⁸O, and positions 1, 7, 8, and 9 for methanol. The spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E₂*) transition of methanol overlaid with the *J* = 1-0 transition of C¹⁸O are shown in Figure A.11. The lines show a single velocity component at the positions 1, 7, 8, and 9, while the C¹⁸O line shows a double peak structure at the positions 2 and 3.

The H₂ column density map of B213-C17 is shown in the upper panel of Figure A.23 with all the offsets observed within B213-C17 in the GEMS large project, however methanol was only towards few offsets (1, 7, 8 and 9). Both the variation of the CH₃OH and C¹⁸O column density ratio, as well as the single column densities N(CH₃OH) and N(C¹⁸O) are shown in the lower panel and show a shallow peak towards the denser offsets observed (7 and 8).
<figure><img src="image_7.png" /><figcaption>Fig. A.1. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*₂) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in TMC-1 C.</figcaption></figure>

<figure><img src="image_8.png" /><figcaption>Fig. A.2. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*₂) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in TMC-1 CP.</figcaption></figure>
<figure><img src="image_9.png" /><figcaption>Fig. A.3. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in TMC-1 NH<sub>3</sub>.</figcaption></figure>

<figure><img src="image_10.png" /><figcaption>Fig. A.4. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C2</figcaption></figure>
<figure><img src="image_11.png" /><figcaption>Fig. A.5. Spectra of the $2_{1,2}-1_{1,1}$ ($E_2$) transition of methanol overlaid with the $J = 1-0$ transition of C$^{18}$O in B213-C5</figcaption></figure>

<figure><img src="image_12.png" /><figcaption>Fig. A.6. Spectra of the $2_{1,2}-1_{1,1}$ ($E_2$) transition of methanol overlaid with the $J = 1-0$ transition of C$^{18}$O in B213-C6</figcaption></figure>
<figure><img src="image_13.png" /><figcaption>Fig. A.7. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C7</figcaption></figure>

<figure><img src="image_14.png" /><figcaption>Fig. A.8. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C10</figcaption></figure>
<figure><img src="image_15.png" /><figcaption>Fig. A.9. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C12</figcaption></figure>

<figure><img src="image_16.png" /><figcaption>Fig. A.10. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C16</figcaption></figure>
<figure><img src="image_17.png" /><figcaption>Fig. A.11. Spectra of the 2<sub>1,2</sub>-1<sub>1,1</sub> (*E*<sub>2</sub>) transition of methanol overlaid with the *J* = 1-0 transition of C<sup>18</sup>O in B213-C17</figcaption></figure>
<figure><img src="image_18.png" /><figcaption>**Fig. A.12.** Upper panel: H₂ column density of TMC-1 C derived from Herschel and Planck data (Rodríguez-Baras et al. 2021). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in TMC-1 C. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_19.png" /><figcaption>**Fig. A.13.** Upper panel: H₂ column density of TMC-1 CP derived from Herschel and Planck data (Rodríguez-Baras et al. 2021). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in TMC-1 CP in the two velocity components observed for methanol. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_20.png" /><figcaption>**Fig. A.14.** Upper panel: H₂ column density of TMC-1 NH₃ derived from Herschel and Planck data (Rodríguez-Baras et al. 2021). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in TMC-1 NH₃ in the two velocity components observed for methanol. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_21.png" /><figcaption>**Fig. A.15.** Upper panel: H₂ column density of B213-C1 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangle shows the position of a Class II/III protostellar core (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C1. In the plot of the CH₃OH and C¹⁸O column density ratio, the A_V in each offset is marked in red.</figcaption></figure>
<figure><img src="image_22.png" /><figcaption>**Fig. A.16.** Upper panel: H₂ column density of B213-C2 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C2. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_23.png" /><figcaption>**Fig. A.17.** Upper panel: H₂ column density of B213-C5 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangles show the positions of Class I/flat (full triangle) or Class II/III (empty triangle) protostellar cores (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C5. In the plot of the CH₃OH and C¹⁸O column density ratio, the A_V in each offset is marked in red.</figcaption></figure>
B213-C6

N(H₂) (cm⁻²)

<figure><img src="image_24.png" /><figcaption>**Fig. A.18.** Upper panel: H₂ column density of B213-C6 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C6. In the plot of the CH₃OH and C¹⁸O column density ratio, the A_V in each offset is marked in red.</figcaption></figure>
<figure><img src="image_25.png" /><figcaption>**Fig. A.19.** Upper panel: H₂ column density of B213-C7 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C7. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_26.png" /><figcaption>**Fig. A.20.** Upper panel: H₂ column density of B213-C10 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangle shows the position of a Class I/flat protostellar core (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C10. In the plot of the CH₃OH and C¹⁸O column density ratio, the A_V in each offset is marked in red.</figcaption></figure>
<figure><img src="image_27.png" /><figcaption>**Fig. A.21.** Upper panel: H₂ column density of B213-C12 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangles show the positions of Class I/flat protostellar cores (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C12. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.</figcaption></figure>
<figure><img src="image_28.png" /></figure>

**Fig. A.22.** Upper panel: H₂ column density of B213-C16 derived from Herschel and Planck data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangles show the positions of Class I/flat (full triangle) or Class II/III (empty triangle) protostellar cores (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C16. In the plot of the CH₃OH and C¹⁸O column density ratio, the A<sub>V</sub> in each offset is marked in red.
<figure><img src="image_29.png" /><figcaption>**Fig. A.23.** Upper panel: H₂ column density of B213-C17 derived from *Herschel* and *Planck* data (Palmeirim et al. 2013). The triangles mark the positions observed with the GEMS large project. The red triangle shows the position of a Class I/flat protostellar core (Rebull et al. 2010). Lower panel: CH₃OH and C¹⁸O column densities, and column density ratios computed in the different offsets in B213-C17. In the plot of the CH₃OH and C¹⁸O column density ratio, the A_V in each offset is marked in red.</figcaption></figure>
## Appendix B: Results - additional material

<table><caption>Table B.1. Results from the MCMC and RADEX analysis.</caption><thead><tr><th rowspan="2"></th><th rowspan="2">A<sub>V</sub><sup>a</sup><br>mag</th><th colspan="5">CH<sub>3</sub>OH</th><th colspan="3">C<sup>18</sup>O</th></tr><tr><th>N<sub>tot</sub><br>[10<sup>13</sup>cm<sup>-2</sup>]</th><th>T<sub>kin</sub><br>[K]</th><th>n<sub>H<sub>2</sub></sub><br>[10<sup>4</sup>cm<sup>-3</sup>]</th><th>v<sub>LSR</sub><br>[km s<sup>-1</sup>]</th><th>δv<br>[km s<sup>-1</sup>]</th><th>N<sub>tot</sub><br>[10<sup>15</sup>cm<sup>-2</sup>]</th><th>v<sub>LSR</sub><br>[km s<sup>-1</sup>]</th><th>δv<br>[km s<sup>-1</sup>]</th></tr></thead><tbody><tr><td colspan="10"><strong>TMC-1 C</strong></td></tr><tr><td>1</td><td>19.9</td><td>2.7<sup>+0.4</sup><sub>-0.30</sub></td><td>10.4<sup>+0.2</sup><sub>-0.3</sub></td><td>3.8<sup>+1.0</sup><sub>-0.7</sub></td><td>5.18</td><td>0.34</td><td>3.62(1)</td><td>5.21</td><td>0.46</td></tr><tr><td>2</td><td>18.5</td><td>2.7<sup>+0.3</sup><sub>-0.3</sub></td><td>10.3<sup>+0.5</sup><sub>-0.12</sub></td><td>2.5<sup>+0.5</sup><sub>-0.5</sub></td><td>5.14</td><td>0.35</td><td>3.3(1)</td><td>5.19</td><td>0.46</td></tr><tr><td>3</td><td>13.3</td><td>2.1<sup>+0.5</sup><sub>-0.2</sub></td><td>12.5<sup>+0.3</sup><sub>-0.3</sub></td><td>1.01<sup>+0.2</sup><sub>-0.3</sub></td><td>5.12</td><td>0.32</td><td>2.67(2)</td><td>5.16</td><td>0.46</td></tr><tr><td>4</td><td>4.8</td><td>1.4<sup>+0.3</sup><sub>-0.5</sub></td><td>11.3<sup>+1.0</sup><sub>-0.6</sub></td><td>0.15<sup>+0.08</sup><sub>-0.02</sub></td><td>5.09</td><td>0.39</td><td>1.54(1)</td><td>5.13</td><td>0.51</td></tr><tr><td colspan="10"><strong>TMC-1 CP</strong></td></tr><tr><td>1</td><td>18.2</td><td>2.0<sup>+0.6</sup><sub>-0.3</sub></td><td>10.5<sup>+0.9</sup><sub>-0.4</sub></td><td>1.5<sup>+0.3</sup><sub>-0.5</sub></td><td>5.63</td><td>0.33</td><td>2.3±0.34<sup>b</sup></td><td></td><td>0.28<sup>c</sup></td></tr><tr><td>1</td><td>18.2</td><td>1.2<sup>+0.1</sup><sub>-0.1</sub></td><td>10.7<sup>+0.2</sup><sub>-0.4</sub></td><td>3.87<sup>+0.30</sup><sub>-0.51</sub></td><td>6.07</td><td>0.32</td><td>0.9±0.13<sup>b</sup></td><td></td><td>0.3<sup>c</sup></td></tr><tr><td>2</td><td>16.7</td><td>1.3<sup>+0.3</sup><sub>-0.2</sub></td><td>12.1<sup>+0.3</sup><sub>-0.5</sub></td><td>1.6<sup>+0.2</sup><sub>-0.5</sub></td><td>5.66</td><td>0.38</td><td>2.2±0.33<sup>b</sup></td><td></td><td>0.28<sup>c</sup></td></tr><tr><td>2</td><td>16.7</td><td>1.1<sup>+0.2</sup><sub>-0.1</sub></td><td>10.3<sup>+0.2</sup><sub>-0.2</sub></td><td>3.5<sup>+0.3</sup><sub>-0.8</sub></td><td>6.06</td><td>0.37</td><td>1.0±0.15<sup>b</sup></td><td></td><td>0.28<sup>c</sup></td></tr><tr><td>3</td><td>13.7</td><td>0.6<sup>+0.1</sup><sub>-0.1</sub></td><td>10.7<sup>+0.2</sup><sub>-0.3</sub></td><td>4.2<sup>+0.9</sup><sub>-0.6</sub></td><td>5.68</td><td>0.36</td><td>1.8±0.27<sup>b</sup></td><td></td><td>0.28<sup>c</sup></td></tr><tr><td>3</td><td>13.7</td><td>0.7<sup>+0.1</sup><sub>-0.1</sub></td><td>10.6<sup>+0.2</sup><sub>-0.5</sub></td><td>4.2<sup>+0.3</sup><sub>-0.4</sub></td><td>6.01</td><td>0.36</td><td>1.2±0.18<sup>b</sup></td><td></td><td>0.3<sup>c</sup></td></tr><tr><td>4</td><td>7.3</td><td>2.3<sup>+1.9</sup><sub>-1.0</sub></td><td>10.3<sup>+0.4</sup><sub>-0.2</sub></td><td>0.16<sup>+0.05</sup><sub>-0.04</sub></td><td>5.67</td><td>0.33</td><td>1.2±0.18<sup>b</sup></td><td></td><td>0.3<sup>c</sup></td></tr><tr><td>4</td><td>7.3</td><td>0.3<sup>+0.1</sup><sub>-0.1</sub></td><td>12.2<sup>+0.3</sup><sub>-0.7</sub></td><td>0.34<sup>+0.07</sup><sub>-0.05</sub></td><td>6.02</td><td>0.33</td><td>0.7±0.10<sup>b</sup></td><td></td><td>0.3<sup>c</sup></td></tr><tr><td colspan="10"><strong>TMC-1 NH<sub>3</sub></strong></td></tr><tr><td>1</td><td>17.0</td><td>1.7<sup>+1.5</sup><sub>-0.3</sub></td><td>11.3<sup>+0.4</sup><sub>-0.4</sub></td><td>1.5<sup>+0.8</sup><sub>-1.0</sub></td><td>5.59</td><td>0.43</td><td>0.7±0.10<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>1</td><td>17.0</td><td>1.0<sup>+0.2</sup><sub>-0.2</sub></td><td>10.6<sup>+0.2</sup><sub>-0.4</sub></td><td>2.8<sup>+0.7</sup><sub>-0.8</sub></td><td>6.13</td><td>0.38</td><td>0.45±0.07<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>2</td><td>15.6</td><td>3.0<sup>+0.7</sup><sub>-0.6</sub></td><td>11.0<sup>+0.4</sup><sub>-0.7</sub></td><td>0.7<sup>+0.2</sup><sub>-0.1</sub></td><td>5.60</td><td>0.37</td><td>0.8±0.12<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>2</td><td>15.6</td><td>3.8<sup>+1.4</sup><sub>-0.7</sub></td><td>10.9<sup>+0.9</sup><sub>-0.5</sub></td><td>0.4<sup>+0.1</sup><sub>-0.1</sub></td><td>6.01</td><td>0.39</td><td>0.5±0.07<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>3</td><td>12.9</td><td>8.8<sup>+0.6</sup><sub>-0.7</sub></td><td>13.8<sup>+0.5</sup><sub>-1.5</sub></td><td>0.35<sup>+0.05</sup><sub>-0.04</sub></td><td>5.75</td><td>0.44</td><td>1.1±0.16<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>3</td><td>12.9</td><td>1.7<sup>+0.5</sup><sub>-0.8</sub></td><td>10.2<sup>+0.2</sup><sub>-0.1</sub></td><td>0.29<sup>+0.19</sup><sub>-0.07</sub></td><td>6.01</td><td>0.35</td><td>0.7±0.10<sup>b</sup></td><td></td><td>0.4<sup>c</sup></td></tr><tr><td>4</td><td>10.0</td><td>3.0<sup>+0.7</sup><sub>-0.4</sub></td><td>10.6<sup>+1.5</sup><sub>-0.4</sub></td><td>0.26<sup>+0.08</sup><sub>-0.06</sub></td><td>5.84</td><td>0.44</td><td>0.9±0.13<sup>b</sup></td><td></td><td>0.45<sup>c</sup></td></tr><tr><td colspan="10"><strong>B213-1</strong></td></tr><tr><td>1</td><td>26.9</td><td>6.9<sup>+0.5</sup><sub>-0.5</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.1<sup>+0.1</sup><sub>-0.1</sub></td><td>5.97</td><td>0.50</td><td>1.3±0.19<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr><tr><td>2</td><td>13.6</td><td>5.8<sup>+0.4</sup><sub>-0.6</sub></td><td>10.4<sup>+0.2</sup><sub>-0.2</sub></td><td>1.1<sup>+0.2</sup><sub>-0.1</sub></td><td>6.19</td><td>0.50</td><td>1.5±0.22<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr><tr><td>3</td><td>13.4</td><td>4.6<sup>+1.3</sup><sub>-1.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.0<sup>+0.6</sup><sub>-0.3</sub></td><td>6.21</td><td>0.48</td><td>1.2±0.18<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr><tr><td>3</td><td>13.4</td><td>3.2<sup>+1.7</sup><sub>-2.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.7<sup>+0.7</sup><sub>-0.7</sub></td><td>6.83</td><td>0.43</td><td>1.2±0.18<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr><tr><td>4</td><td>6.0</td><td>1.3<sup>+0.5</sup><sub>-0.3</sub></td><td>10.2<sup>+0.2</sup><sub>-0.1</sub></td><td>1.1<sup>+0.7</sup><sub>-0.4</sub></td><td>6.62</td><td>0.66</td><td>1.3±0.19<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr><tr><td>5</td><td>4.2</td><td>0.4<sup>+0.3</sup><sub>-0.2</sub></td><td>10.3<sup>+0.1</sup><sub>-0.1</sub></td><td>1.0<sup>+0.8</sup><sub>-0.4</sub></td><td>6.53</td><td>0.47</td><td>0.9±0.13<sup>b</sup></td><td></td><td>0.6<sup>c</sup></td></tr></tbody></table>
<table><tbody><tr><td>6</td><td>5.2</td><td>0.8<sup>+0.4</sup><sub>-0.2</sub></td><td>10.3<sup>+0.4</sup><sub>-0.3</sub></td><td>1.1<sup>+0.4</sup><sub>-0.4</sub></td><td>6.65</td><td>0.43</td><td>0.71(1)</td><td>6.63</td><td>0.47</td></tr><tr><td>7</td><td>10.5</td><td>0.7<sup>+0.2</sup><sub>-0.2</sub></td><td>10.3<sup>+0.15</sup><sub>-0.16</sub></td><td>4.4<sup>+1.5</sup><sub>-2.1</sub></td><td>6.80</td><td>0.48</td><td>1.17(1)</td><td>6.84</td><td>0.64</td></tr><tr><td>8</td><td>10.3</td><td>3.1<sup>+0.3</sup><sub>-0.2</sub></td><td>10.1<sup>+0.12</sup><sub>-0.08</sub></td><td>5.6<sup>+1.5</sup><sub>-1.2</sub></td><td>5.97</td><td>0.49</td><td>1.24(3)</td><td>6.06</td><td>0.61</td></tr><tr><td>9</td><td>5.5</td><td>2.8<sup>+0.2</sup><sub>-0.2</sub></td><td>10.7<sup>+0.16</sup><sub>-0.19</sub></td><td>6.5<sup>+0.8</sup><sub>-1.5</sub></td><td>5.98</td><td>0.49</td><td>1.26(3)</td><td>6.02</td><td>0.66</td></tr><tr><td>B213-2</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>20.9</td><td>3.5<sup>+2.2</sup><sub>-1.7</sub></td><td>10.5<sup>+0.7</sup><sub>-0.4</sub></td><td>2.4<sup>+1.0</sup><sub>-1.0</sub></td><td>7.29</td><td>0.55</td><td>0.5±0.07</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>2</td><td>10.6</td><td>0.7<sup>+0.2</sup><sub>-0.2</sub></td><td>10.3<sup>+0.1</sup><sub>-0.2</sub></td><td>4.4<sup>+0.4</sup><sub>-1.3</sub></td><td>7.02</td><td>0.46</td><td>0.6±0.09</td><td><sup>b</sup></td><td>0.6<sup>c</sup></td></tr><tr><td>8</td><td>16.5</td><td>1.5<sup>+1.1</sup><sub>-0.4</sub></td><td>10.9<sup>+0.2</sup><sub>-0.3</sub></td><td>2.4<sup>+1.1</sup><sub>-1.0</sub></td><td>6.97</td><td>0.40</td><td>0.8±0.12</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>9</td><td>9.8</td><td>1.4<sup>+0.2</sup><sub>-0.2</sub></td><td>10.2<sup>+0.2</sup><sub>-0.1</sub></td><td>6.2<sup>+0.7</sup><sub>-2.6</sub></td><td>6.85</td><td>0.33</td><td>0.9±0.05</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>B213-5</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>23.6</td><td>7.8<sup>+0.48</sup><sub>-0.9</sub></td><td>10.2<sup>+0.1</sup><sub>-0.1</sub></td><td>1.1<sup>+0.2</sup><sub>-0.1</sub></td><td>6.26</td><td>0.50</td><td>1.5±0.22</td><td><sup>b</sup></td><td>0.7<sup>c</sup></td></tr><tr><td>2</td><td>13.0</td><td>3.0<sup>+0.5</sup><sub>-0.5</sub></td><td>10.2<sup>+0.3</sup><sub>-0.1</sub></td><td>1.2<sup>+0.3</sup><sub>-0.1</sub></td><td>6.44</td><td>0.70</td><td>1.0±0.15</td><td><sup>b</sup></td><td>0.7<sup>c</sup></td></tr><tr><td>3</td><td>10.9</td><td>1.0<sup>+0.3</sup><sub>-0.1</sub></td><td>10.1<sup>+0.06</sup><sub>-0.07</sub></td><td>2.7<sup>+0.5</sup><sub>-1.5</sub></td><td>6.73</td><td>0.68</td><td>1.0±0.15</td><td><sup>b</sup></td><td>0.7<sup>c</sup></td></tr><tr><td>4</td><td>13.3</td><td>0.3<sup>+0.1</sup><sub>-0.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.0<sup>+0.1</sup><sub>-0.1</sub></td><td>6.37</td><td>0.44</td><td>0.57(1)</td><td>6.3</td><td>0.37</td></tr><tr><td>4</td><td>13.3</td><td>0.8<sup>+0.2</sup><sub>-0.2</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>0.96<sup>+0.2</sup><sub>-0.1</sub></td><td>7.07</td><td>0.37</td><td>4.72(1)</td><td>7.11</td><td>0.55</td></tr><tr><td>5</td><td>14.0</td><td>0.2<sup>+0.1</sup><sub>-0.1</sub></td><td>10.0<sup>+0.1</sup><sub>-0.1</sub></td><td>0.6<sup>+0.6</sup><sub>-0.1</sub></td><td>6.28</td><td>0.42</td><td>0.63(1)</td><td>6.33</td><td>0.35</td></tr><tr><td>5</td><td>14.0</td><td>0.4<sup>+0.1</sup><sub>-0.1</sub></td><td>10.4<sup>+0.1</sup><sub>-0.1</sub></td><td>1.9<sup>+0.1</sup><sub>-0.2</sub></td><td>7.09</td><td>0.37</td><td>4.88(2)</td><td>7.13</td><td>0.55</td></tr><tr><td>6</td><td>14.5</td><td>1.2<sup>+0.7</sup><sub>-0.4</sub></td><td>10.4<sup>+0.17</sup><sub>-0.15</sub></td><td>4.8<sup>+2.0</sup><sub>-1.2</sub></td><td>7.25</td><td>0.55</td><td>6.14(2)</td><td>7.15</td><td>0.60</td></tr><tr><td>7</td><td>18.9</td><td>1.8<sup>+1.7</sup><sub>-1.1</sub></td><td>10.7<sup>+0.7</sup><sub>-0.3</sub></td><td>5.5<sup>+0.8</sup><sub>-0.9</sub></td><td>7.3</td><td>0.51</td><td>7.75(2)</td><td>7.21</td><td>0.59</td></tr><tr><td>8</td><td>12.9</td><td>2.9<sup>+0.3</sup><sub>-0.4</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.2<sup>+0.2</sup><sub>-0.1</sub></td><td>6.36</td><td>0.50</td><td>1.2±0.18</td><td><sup>b</sup></td><td>0.6<sup>c</sup></td></tr><tr><td>9</td><td>10.9</td><td>0.8<sup>+0.3</sup><sub>-0.1</sub></td><td>10.3<sup>+0.1</sup><sub>-0.1</sub></td><td>1.3<sup>+0.4</sup><sub>-0.6</sub></td><td>6.06</td><td>0.48</td><td>0.77(1)</td><td>6.12</td><td>0.52</td></tr><tr><td>9</td><td>10.9</td><td>0.3<sup>+0.2</sup><sub>-0.1</sub></td><td>10.2<sup>+0.1</sup><sub>-0.1</sub></td><td>1.0<sup>+0.6</sup><sub>-0.4</sub></td><td>7.00</td><td>0.37</td><td>3.06(1)</td><td>7.11</td><td>0.55</td></tr><tr><td>B213-6</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>22.2</td><td>3.5<sup>+0.5</sup><sub>-0.4</sub></td><td>10.5<sup>+0.3</sup><sub>-0.2</sub></td><td>2.5<sup>+1.3</sup><sub>-0.6</sub></td><td>6.85</td><td>0.53</td><td>0.69(1)</td><td>6.77</td><td>0.68</td></tr><tr><td>2</td><td>15</td><td>2.4<sup>+0.3</sup><sub>-0.3</sub></td><td>10.6<sup>+0.1</sup><sub>-0.1</sub></td><td>4.8<sup>+0.6</sup><sub>-1.0</sub></td><td>6.77</td><td>0.40</td><td>0.83(1)</td><td>6.79</td><td>0.46</td></tr><tr><td>3</td><td>5.6</td><td>0.7<sup>+0.5</sup><sub>-0.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>3.3<sup>+1.4</sup><sub>-2.3</sub></td><td>6.80</td><td>0.39</td><td>1.11(1)</td><td>6.84</td><td>0.44</td></tr><tr><td>8</td><td>18.2</td><td>4.0<sup>+0.9</sup><sub>-0.5</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>1.7<sup>+0.6</sup><sub>-0.6</sub></td><td>6.94</td><td>0.37</td><td>0.73(1)</td><td>6.92</td><td>0.49</td></tr><tr><td>9</td><td>11.3</td><td>2.5<sup>+2.5</sup><sub>-2.1</sub></td><td>10.7<sup>+0.6</sup><sub>-0.6</sub></td><td>5.2<sup>+3.3</sup><sub>-3.0</sub></td><td>7.06</td><td>0.40</td><td>0.76(2)</td><td>6.99</td><td>0.44</td></tr><tr><td>B213-7</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>20.2</td><td>9.5<sup>+2.7</sup><sub>-1.4</sub></td><td>10.2<sup>+0.1</sup><sub>-0.1</sub></td><td>0.48<sup>+0.15</sup><sub>-0.13</sub></td><td>6.91</td><td>0.44</td><td>2.08(1)</td><td>6.87</td><td>0.41</td></tr><tr><td>2</td><td>14.6</td><td>2.7<sup>+0.3</sup><sub>-0.3</sub></td><td>10.4<sup>+0.1</sup><sub>-0.2</sub></td><td>7.4<sup>+0.5</sup><sub>-3.0</sub></td><td>6.92</td><td>0.41</td><td>2.38(2)</td><td>6.91</td><td>0.47</td></tr><tr><td>3</td><td>7.2</td><td>6<sup>+1</sup><sub>-1</sub></td><td>10.1<sup>+0.2</sup><sub>-0.1</sub></td><td>4.5<sup>+1.7</sup><sub>-1.2</sub></td><td>6.96</td><td>0.43</td><td>1.87(2)</td><td>6.96</td><td>0.55</td></tr></tbody></table>
<table><tbody><tr><td>8</td><td>16.6</td><td>2.5<sup>+0.3</sup><sub>-0.2</sub></td><td>10.3<sup>+0.3</sup><sub>-0.2</sub></td><td>8.9<sup>+0.7</sup><sub>-1.5</sub></td><td>6.91</td><td>0.47</td><td>1.82(2)</td><td>6.84</td><td>0.43</td></tr><tr><td>9</td><td>8.7</td><td>7<sup>+1</sup><sub>-4</sub></td><td>11.0<sup>+0.1</sup><sub>-0.8</sub></td><td>0.31<sup>+0.67</sup><sub>-0.08</sub></td><td>6.80</td><td>0.34</td><td>1.51(1)</td><td>6.88</td><td>0.52</td></tr><tr><td>B213-10</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>20.7</td><td>1.3<sup>+0.2</sup><sub>-0.3</sub></td><td>10.2<sup>+0.2</sup><sub>-0.1</sub></td><td>1.2<sup>+0.7</sup><sub>-0.3</sub></td><td>6.76</td><td>0.33</td><td>1.17(2)</td><td>6.77</td><td>0.45</td></tr><tr><td>2</td><td>17.5</td><td>5.1<sup>+1.4</sup><sub>-1.2</sub></td><td>10.5<sup>+0.2</sup><sub>-0.2</sub></td><td>0.27<sup>+0.13</sup><sub>-0.08</sub></td><td>6.69</td><td>0.27</td><td>1.34(2)</td><td>6.68</td><td>0.47</td></tr><tr><td>3</td><td>11.8</td><td>4.9<sup>+0.8</sup><sub>-1.0</sub></td><td>10.3<sup>+0.2</sup><sub>-0.1</sub></td><td>0.27<sup>+0.11</sup><sub>-0.05</sub></td><td>6.60</td><td>0.31</td><td>1.36(2)</td><td>6.57</td><td>0.47</td></tr><tr><td>4</td><td>6.1</td><td>1.2<sup>+1.0</sup><sub>-0.4</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>0.72<sup>+0.45</sup><sub>-0.39</sub></td><td>6.58</td><td>0.32</td><td>0.90(1)</td><td>6.59</td><td>0.45</td></tr><tr><td>7</td><td>17.6</td><td>2<sup>+1</sup><sub>-1</sub></td><td>10.6<sup>+0.2</sup><sub>-0.4</sub></td><td>0.58<sup>+0.62</sup><sub>-0.21</sub></td><td>6.82</td><td>0.32</td><td>1.11(1)</td><td>6.83</td><td>0.42</td></tr><tr><td>8</td><td>14.5</td><td>0.8<sup>+0.2</sup><sub>-0.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>5.5<sup>+0.6</sup><sub>-2.0</sub></td><td>6.89</td><td>0.31</td><td>1.31(1)</td><td>6.86</td><td>0.40</td></tr><tr><td>9</td><td>9.1</td><td>0.6<sup>+0.1</sup><sub>-0.1</sub></td><td>10.1<sup>+0.1</sup><sub>-0.8</sub></td><td>4.9<sup>+0.4</sup><sub>-0.7</sub></td><td>6.73</td><td>0.29</td><td>1.07(2)</td><td>6.74</td><td>0.42</td></tr><tr><td>B213-12</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>22.1</td><td>1.4<sup>+0.4</sup><sub>-0.3</sub></td><td>10.2<sup>+0.1</sup><sub>-0.1</sub></td><td>0.83<sup>+0.3</sup><sub>-0.2</sub></td><td>6.65</td><td>0.39</td><td>0.5±0.07</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>2</td><td>17.6</td><td>2<sup>+1</sup><sub>-1</sub></td><td>10.8<sup>+0.2</sup><sub>-0.2</sub></td><td>0.3<sup>+0.2</sup><sub>-0.1</sub></td><td>6.69</td><td>0.39</td><td>0.45±0.07</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>3</td><td>11.5</td><td>6<sup>+1</sup><sub>-4</sub></td><td>10.1<sup>+0.1</sup><sub>-0.1</sub></td><td>0.13<sup>+0.78</sup><sub>-0.24</sub></td><td>6.69</td><td>0.37</td><td>0.45±0.07</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>8</td><td>8.7</td><td>2.0<sup>+2.0</sup><sub>-0.8</sub></td><td>10.5<sup>+0.21</sup><sub>-0.12</sub></td><td>0.41<sup>+0.29</sup><sub>-0.19</sub></td><td>6.51</td><td>0.32</td><td>0.7±0.10</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>B213-16</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>24.8</td><td>1.8<sup>+0.1</sup><sub>-0.11</sub></td><td>10.1<sup>+0.09</sup><sub>-0.07</sub></td><td>5.1<sup>+0.6</sup><sub>-1.0</sub></td><td>6.68</td><td>0.34</td><td>1.3±0.19</td><td><sup>b</sup></td><td>0.4<sup>c</sup></td></tr><tr><td>2</td><td>13.3</td><td>1.2<sup>+0.1</sup><sub>-0.1</sub></td><td>10.3<sup>+0.11</sup><sub>-0.10</sub></td><td>4.3<sup>+0.7</sup><sub>-1.1</sub></td><td>6.61</td><td>0.44</td><td>1.5±0.22</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>3</td><td>10.5</td><td>3.3<sup>+0.1</sup><sub>-0.1</sub></td><td>10.2<sup>+0.14</sup><sub>-0.15</sub></td><td>0.37<sup>+0.1</sup><sub>-0.1</sub></td><td>6.48</td><td>0.43</td><td>1.5±0.22</td><td><sup>b</sup></td><td>0.5<sup>c</sup></td></tr><tr><td>4</td><td>6.9</td><td>3.3<sup>+0.2</sup><sub>-0.1</sub></td><td>10.2<sup>+0.07</sup><sub>-0.07</sub></td><td>0.16<sup>+0.1</sup><sub>-0.1</sub></td><td>6.35</td><td>0.40</td><td>1.2±0.18</td><td><sup>b</sup></td><td>0.4<sup>c</sup></td></tr><tr><td>7</td><td>22.6</td><td>1.8<sup>+0.1</sup><sub>-0.1</sub></td><td>10.34<sup>+0.12</sup><sub>-0.17</sub></td><td>4.5<sup>+0.7</sup><sub>-0.9</sub></td><td>6.71</td><td>0.33</td><td>1.2±0.18</td><td><sup>b</sup></td><td>0.35<sup>c</sup></td></tr><tr><td>8</td><td>9.6</td><td>1.3<sup>+0.3</sup><sub>-0.2</sub></td><td>10.5<sup>+0.17</sup><sub>-0.20</sub></td><td>2.4<sup>+0.8</sup><sub>-1.0</sub></td><td>6.73</td><td>0.31</td><td>1.2±0.18</td><td><sup>b</sup></td><td>0.35<sup>c</sup></td></tr><tr><td>9</td><td>5.2</td><td>1.5<sup>+0.9</sup><sub>-0.6</sub></td><td>10.3<sup>+0.46</sup><sub>-0.22</sub></td><td>0.32<sup>+0.39</sup><sub>-0.15</sub></td><td>6.56</td><td>0.41</td><td>1.0±0.15</td><td><sup>b</sup></td><td>0.35<sup>c</sup></td></tr><tr><td>B213-17</td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td><td></td></tr><tr><td>1</td><td>6.4</td><td>0.4<sup>+0.3</sup><sub>-0.2</sub></td><td>10.4<sup>+0.3</sup><sub>-0.3</sub></td><td>1.4<sup>+1.2</sup><sub>-0.7</sub></td><td>7.24</td><td>0.42</td><td>0.8±0.12</td><td><sup>b</sup></td><td>0.4<sup>c</sup></td></tr><tr><td>7</td><td>12.0</td><td>4.1<sup>+2.8</sup><sub>-3.7</sub></td><td>10.3<sup>+0.2</sup><sub>-0.2</sub></td><td>7.4<sup>+0.9</sup><sub>-2.1</sub></td><td>7.39</td><td>0.42</td><td>0.87(9)</td><td>7.22</td><td>0.52</td></tr><tr><td>8</td><td>19.7</td><td>5.6<sup>+3.0</sup><sub>-2.2</sub></td><td>10.5<sup>+0.4</sup><sub>-0.3</sub></td><td>4.5<sup>+1.5</sup><sub>-2.4</sub></td><td>7.31</td><td>0.42</td><td>0.83(1)</td><td>7.21</td><td>0.48</td></tr><tr><td>9</td><td>10.2</td><td>0.8<sup>+0.6</sup><sub>-0.3</sub></td><td>10.2<sup>+0.2</sup><sub>-0.1</sub></td><td>3.7<sup>+1.1</sup><sub>-3.0</sub></td><td>7.07</td><td>0.41</td><td>0.66(1)</td><td>7.13</td><td>0.51</td></tr></tbody></table>

Notes. (a) Derived from Herschel and Planck data (Palmeirim et al. 2013; Rodríguez-Baras et al. 2021). (b) Fixed to the correspondent $v_{LRS}$ of CH$_3$OH. (c) Fixed to a value that matches best the $\delta v$ of CH$_3$OH and the line shape of C$^{18}$O. $v_{LSR}$ and $\delta v$ are reported without uncertainty because the statistical uncertainty is much lower than the spectral resolution.
<figure><img src="image_30.png" /><figcaption>**Fig. B.1.** A_V and T_dust in the B213 cores within our sample are plotted against each other to show the effect of the larger stellar activity in the northern part of B213. In this Figure, the cores in the north of B213 are plotted in different shades of red that go from darker to lighter when moving from north to south. The cores in the central and southern part of B213 are plotted in different shades of blue, with the darkest being the core located more to the south.</figcaption></figure>

<figure><img src="image_31.png" /><figcaption>**Fig. B.2.** H₂ column density map of the B213 filament derived from Herschel and Planck data (Palmeirim et al. 2013). The young stellar objects reported in Rebull et al. (2010) are shown as triangles.</figcaption></figure>
<figure><img src="image_32.png" /><figcaption>**Fig. B.3.** Results of the spherical 1D static chemical model using the physical structure derived for B213 C2 and C6 (left panel), C10 and C12 (central vertical panel) and C16 (right panel). The column density profiles of CH₃OH and C¹⁸O (red and blue points) as well as their ratio (black points) have been computed at 10⁵ years (upper panel), 10⁶ years (lower panel).</figcaption></figure>