# ALMA observations of the young protostellar system Barnard 1b: Signatures of an incipient hot corino in B1b-S

N. Marcelino, M. Gerin, J. Cernicharo, A. Fuente, H. Wootten, E. Chapillon, J. Pety, D. Lis, E. Roueff, B. Commerçon, et al.

► To cite this version:

N. Marcelino, M. Gerin, J. Cernicharo, A. Fuente, H. Wootten, et al., ALMA observations of the young protostellar system Barnard 1b: Signatures of an incipient hot corino in B1b-S. Astronomy and Astrophysics - A&A, 2018, 620, pp.A80. 10.1051/0004-6361/201731955. hal-02371327

HAL Id: hal-02371327

https://hal.science/hal-02371327

Submitted on 2 Sep 2021

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

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

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

Distributed under a Creative Commons Attribution 4.0 International License
# ALMA observations of the young protostellar system Barnard 1b: Signatures of an incipient hot corino in B1b-S*

N. Marcelino<fnref n="1" />, M. Gerin<fnref n="2" />, J. Cernicharo<fnref n="1" />, A. Fuente<fnref n="3" />, H. A. Wootten<fnref n="4" />, E. Chapillon<fnref n="5" />,<fnref n="6" />, J. Pety<fnref n="5" />,<fnref n="2" />, D. C. Lis<fnref n="2" />, E. Roueff<fnref n="2" />,<fnref n="7" />, B. Commerçon<fnref n="8" />, and A. Ciardi<fnref n="2" />

Received 14 September 2017 / Accepted 17 September 2018

## ABSTRACT

The Barnard 1b core shows signatures of being at the earliest stages of low-mass star formation, with two extremely young and deeply embedded proto-stellar objects. Hence, this core is an ideal target to study the structure and chemistry of the first objects formed in the collapse of prestellar cores. We present ALMA Band 6 spectral line observations at ~0.6" of angular resolution towards Barnard 1b. We have extracted the spectra towards both protostars, and used a local thermodynamic equilibrium (LTE) model to reproduce the observed line profiles. B1b-S shows rich and complex spectra, with emission from high energy transitions of complex molecules, such as CH₃OCOH and CH₃CHO, including vibrational level transitions. We have tentatively detected for the first time in this source emission from NH₂CN, NH₂CHO, CH₃CH₂OH, CH₂OHCHO, CH₃CH₂OCOH and both *aGg'* and *gGg'* conformers of (CH₂OH)₂. This is the first detection of ethyl formate (CH₃CH₂OCOH) towards a low-mass star forming region. On the other hand, the spectra of the FHSC candidate B1b-N are free of COMs emission. In order to fit the observed line profiles in B1b-S, we used a source model with two components: an inner hot and compact component (200 K, 0.35") and an outer and colder one (60 K, 0.6"). The resulting COM abundances in B1b-S range from 10⁻¹³ for NH₂CN and NH₂CHO, up to 10⁻⁹ for CH₃OCOH. Our ALMA Band 6 observations reveal the presence of a compact and hot component in B1b-S, with moderate abundances of complex organics. These results indicate that a hot corino is being formed in this very young Class 0 source.

**Key words.** astrochemistry – ISM: clouds – ISM: individual objects: Barnard 1b – ISM: abundances – stars: formation – stars: low-mass

## 1. Introduction

In the innermost regions of low-mass Class 0 objects, gas and dust can reach temperatures of >100 K, at which water and other grain mantle components can evaporate into the gas phase. In particular, the large contribution of complex organic molecules (COMs; defined as species with six or more atoms) to their spectra, makes them similar to hot cores in massive star forming regions, and they are therefore known as hot corinos (Herbst & van Dishoeck 2009; Caselli & Ceccarelli 2012). These regions have small sizes, of <100 au, and their study has greatly improved in the last few years thanks to interferometers operating at millimetre and sub-millimetre wavelengths such as NOEMA and ALMA. In particular, ALMA observations at unprecedented angular resolution, have allowed us to probe the inner hot regions of low-mass protostars and compare with sophisticated MHD models (see, e.g. Commerçon et al. 2012; Hincelin et al. 2016).

Of particular interest are the earliest stages of low-mass star formation: the initial collapse, the formation of accretion disks

and the launching of protostellar outflows. The Barnard 1 cloud is a moderately active star forming region in the Perseus molecular complex, at a distance of 230 pc (Hirota et al. 2008). It contains several dense cores at different stages of star formation, two of them, B1a and B1c, hosting class 0 sources associated with high velocity outflows (Hatchell et al. 2005, 2007). The proto-stellar core B1b contains three remarkable sources: the infrared source B1b-W, also observed recently at millimetre and sub-millimetre wavelengths (Gerin et al. 2017; Cox et al. 2018), and two extremely young protostellar objects, B1b-N and B1b-S. These two sources have been suggested to be first hydrostatic core (FHSC) candidates, which is the first object to be formed in the collapse of a prestellar core (Larson 1969). Therefore, this core has been a target of multiple studies (see, e.g. Pezzuto et al. 2012; Hirano & Liu 2014). Recent interferometric observations with NOEMA and ALMA have shown the incipient star formation activity in B1b. Gerin et al. (2015) observed young and slow outflows driven by both sources, and obtained outflow properties consistent with B1b-S being a very early Class 0 and B1b-N a FHSC candidate. ALMA continuum emission at 350 GHz at a very high resolution (0.1"), indicates the presence of small compact structures that were interpreted as nascent disks by

<footnote n="1">¹ Instituto de Física Fundamental, CSIC, C/ Serrano 123, 28006 Madrid, Spain
e-mail: nuria.marcelino@csic.es</footnote>

<footnote n="2">² Sorbonne Université, Observatoire de Paris, Université PSL, École normale supérieure, CNRS, LERMA, 75014 Paris, France</footnote>

<footnote n="3">³ Observatorio Astronómico Nacional (OAN, IGN), Apdo 112, 28803 Alcalá de Henares, Spain</footnote>

<footnote n="4">⁴ National Radio Astronomy Observatory, 520 Edgemont Road, Charlottesville, VA 22903, USA</footnote>

<footnote n="5">⁵ Institut de Radioastronomie Millimétrique (IRAM), 300 rue de la Piscine, 38406 Saint Martin d'Hères, France</footnote>

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

<footnote n="7">⁷ Sorbonne Université, Observatoire de Paris, Université PSL, CNRS, LERMA, 92190 Meudon, France</footnote>

<footnote n="8">⁸ Univ. Lyon, ENS de Lyon, Univ. Lyon1, CNRS, Centre de Recherche Astrophysique de Lyon UMR5574, 69007 Lyon, France</footnote>

<footnote>* The full Table B.2 is only available at the CDS via anonymous ftp to cdsarc.u-strasbg.fr (130.79.128.5) or via http://cdsarc.u-strasbg.fr/viz-bin/qcat?J/A+A/620/A80</footnote>
Gerin et al. (2017). Using NOEMA at ~4'', Fuente et al. (2017) imaged the emission of NH<sub>2</sub>D, which is tracing the pseudo-disk, and detected rotation around B1b-S. Here we present spectral line ALMA observations at 0.6'' angular resolution of the two young proto-stars B1b-N and B1b-S.

## 2. Observations

ALMA Cycle 3 observations were performed on 2016 June 21–24, in four different tracks of about an hour each. The ALMA 12 m array was used in the C40-4 configuration, with 36–40 antennas and baselines between 15.1 and 704.1 m. The mean angular resolution achieved (interferometric synthesised beam) at the observed frequencies is ~0.7 × 0.55" (160 × 126 au), and the largest scales to which we are sensitive are ~4". We used the Band 6 receiver centred at 218 and 233 GHz in the lower and upper side band, respectively, connected to the ALMA correlator in FDM mode. In total, 13 spectral windows were used, one at low resolution (~40 km s<sup>-1</sup> and 2 GHz of bandwidth) for the continuum, and 12 windows with spectral resolutions of ~0.16–0.2 km s<sup>-1</sup> and a 58.5 MHz bandwidth.

We performed a seven-point mosaic centred at $\alpha_{J2000} = 03^h33^m21.3^s$ and $\delta_{J2000} = 31^\circ07'34.0''$, covering an area of 40 arcsec<sup>2</sup> containing B1b-N, B1b-S, and the infrared source B1b-W. All the fields were observed in scans of approximately seven minutes duration, with 30 s scans for phase calibration inbetween. Total time on-source is 44 min per track. The phase calibrator used is J0336+3218, located at a distance of ~1.36° of the B1b cores. Two other quasars were observed for bandpass (J0237+2848) and flux (J0238+1636) calibration. Flux density for J0238+1636 was set to ~1.14–1.22 Jy, depending on the day and frequency, based on close flux measurements. Calibration was performed by the ALMA Pipeline in CASA<sup>1</sup> 4.5.3. Images of all spectral windows were also produced in CASA, with a spectral resolution of 0.2 km s<sup>-1</sup> and Briggs weighting. Table B.1 gives the frequency ranges, the final beam and rms of the images. For analysis purposes we used the GILDAS<sup>2</sup> software. We have extracted the spectra, within one synthesised beam size, towards the positions of the B1b-N and B1b-S cores. Depending on the frequency and the beam for each spectral window, we have converted the flux into brightness temperature to analyse the spectra (K Jy<sup>-1</sup> ~ 63).

## 3. Results

The obtained spectra towards the two sources look very different: while B1b-S shows intense and rich spectra, B1b-N is mostly free of line emission. The right-hand panels in Fig. 1 show one of the observed spectral windows, between 215 810 and 215 870 MHz, towards B1b-N (top) and B1b-S (bottom), and Figs. B.1–B.3 present the full spectral coverage. In the northern source, we observe lines from <sup>12</sup>CO, <sup>13</sup>CO, C<sup>18</sup>O, <sup>13</sup>CS, DCO<sup>+</sup>, DCN, N<sub>2</sub>D<sup>+</sup>, and SO. All of them show extended emission and towards the position of the protostar present complex line profiles, which in some cases are only seen in absorption (e.g. <sup>12</sup>CO, and SO). One transition of H<sub>2</sub>C<sup>17</sup>O is marginally detected in emission. All these species are also observed in B1b-S, together with lines from H<sub>2</sub>CO, H<sub>2</sub><sup>13</sup>CO, D<sub>2</sub>CO, CH<sub>3</sub>OH, <sup>13</sup>CH<sub>3</sub>OH, CH<sub>2</sub>DOH, CH<sub>3</sub>OD and several COMs. In total, we have detected 190 lines towards B1b-S. Of those, 113 spectral features have been assigned to 155 transitions from 36 species (including isotopologues and vibrational states). Table B.2 lists

all assigned transitions and their associated species in B1b-S, with line parameters obtained from Gaussian line fits. Line identification was carried out using MADEX<sup>3</sup> (Cernicharo 2012). The other 77 lines for which we could not find reasonable carriers are labeled as unidentified (~40% of the total number of lines). Of those, 16 have also been observed in Orion KL, that is only 21% of the unidentified lines. All transitions listed in Table B.2 are above the 3$\sigma$ level (with $\sigma = \text{rms} \sqrt{Dv\dot{v}}$, $Dv$ the linewidth, and $\dot{v}$ the spectral resolution), while U lines were considered above 4$\sigma$. Several lines, mainly from the most abundant species and low energy transitions, show extended emission when looking at the maps and present complex line profiles towards the position of B1b-S (lines in absorption, inverse P Cygni profiles, line wings, etc.). These lines are indicated in Table B.2 by CP or Ext (for complex profile and extended, respectively), and we do not include them in the analysis below. While we have identified all the lines in the observed spectra (see Figs. B.1–B.3 and Table B.2), the analysis of the extended emission is outside the scope of this paper. Furthermore, it is likely that a significant fraction of the flux is resolved out in the ALMA 12 m-array observations.

In the following we focus on the emission from COMs in the protostar B1b-S. We have detected multiple transitions of acetaldehyde (CH<sub>3</sub>CHO) and methyl formate (CH<sub>3</sub>OCOH), including emission from their deuterated substitutions and vibrational states. We have also observed emission from isocyanic acid (HNCO), cyanamide (NH<sub>2</sub>CN), formamide (NH<sub>2</sub>CHO), vinyl cyanide (CH<sub>2</sub>CHCN), propenal (t-CH<sub>2</sub>CHCHO), ethanol (t-CH<sub>3</sub>CH<sub>2</sub>OH), dimethyl ether (CH<sub>3</sub>OCH<sub>3</sub>), ethyl formate (t-CH<sub>3</sub>CH<sub>2</sub>OCOH), glycolaldehyde (CH<sub>2</sub>OHCHO) and both *a*Gg' and *g*Gg' conformers of ethylene glycol ((CH<sub>2</sub>OH)<sub>2</sub>). However, some of these should be considered tentative detections. In particular, we have detected only one unblended line of HNCO, NH<sub>2</sub>CN, NH<sub>2</sub>CHO, t-CH<sub>2</sub>CHCHO, CH<sub>2</sub>OHCHO, t-CH<sub>3</sub>CH<sub>2</sub>OH, and CH<sub>3</sub>OCH<sub>3</sub>. Of the 16 transitions observed from the ethylene glycol conformers, 8 are blended with transitions from other species. We clearly detect three unblended lines from CH<sub>3</sub>CH<sub>2</sub>OCOH, which has only been previously detected towards Sgr B2, Orion KL, and the hot core in W51 (Belloche et al. 2009; Tercero et al. 2013; Rivilla et al. 2017). We would like to point out the large number of unassigned lines and the high level of blending and complexity in the spectra (see Figs. B.1–B.3 and Table B.2), similar to that observed typically in hot cores and hot corinos. The emission of all these complex species is unresolved in the ALMA images. We have checked that at one beam separation from the position of the B1b-S protostar the spectra are free of COMs emission.

## 4. Column densities and abundances

We have used MADEX (Cernicharo 2012) to reproduce the observed spectra of COMs in B1b-S, and to obtain column densities and kinetic temperatures of the emitting region. We assumed local thermodynamic equilibrium (LTE) conditions, since the densities at the centre of the core are expected to be very high (Gerin et al. 2017). We have focused on CH<sub>3</sub>OCOH and CH<sub>3</sub>CHO, and their isotopic substitutions and vibrationally excited states, for which many transitions are available, to obtain the source model which was then used for the other observed species.

Methyl formate is the species for which the largest number of transitions were identified. In total 24 transitions from both A and E species were observed (18 unblended), with energy

<footnote n="1">¹ https://casa.nrao.edu/</footnote>

<footnote n="2">² http://www.iram.fr/IRAMFR/GILDAS</footnote>

<footnote n="3">³ https://nanocosmos.iff.csic.es/?$page_id=1619</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** *Left panel:* Continuum image at 232 GHz (1.290 mm). The colour scale ranges from 0.001 to 0.2 Jy beam⁻¹. The contour levels are indicated on the colour bar. The ALMA beam is plotted on the lower left corner. *Right panel:* Observed spectra between 215 810 and 215 870 MHz, towards B1b-N (top panel) and B1b-S (bottom panel). Brightness temperature scale, in K, is given on the left. Line features with no labels are unidentified lines. The LTE simulations are over-plotted in red for B1b-S.</figcaption></figure>

levels between 100 and 500 K, allowing to constrain the prop-
erties of the emitting region from the LTE model. We ran a
set of models using one source component of 0.6" in size (the
angular resolution of the ALMA band 6 observations), and we
modified the temperature and column density to fit the observed
spectra. However, it was difficult to fit all the lines with only
one component. The problem became more evident when using
the deuterated substitutions, CH₃OCOD and CH₂DOCOH, and
the vibrationally excited states νt = 1 and νt = 2. We have
detected 12 transitions from CH₃OCOD and four transitions
from CH₂DOCOH (three and two, respectively, are blended with
other species), with upper energy levels between 90 and 170 K.
In both cases the model gave better results when decreasing the
temperature with respect to the best CH₃OCOH model. On the
other hand, the vibrationally excited lines (13 unblended transi-
tions from νt = 1, and one transition from νt = 2), supported
the high temperature model. Therefore, we built a model of the
source with two components. One warm and compact compo-
nent at 200 K and a size of 0.35" (that of the continuum compact
emission obtained by ALMA Band 7 observations, Gerin et al.
2017), and another larger (0.6", approximately the angular res-
olution of the Band 6 observations) and colder component at
60 K. Figure B.4 shows two examples of the best one compo-
nent models obtained using the main CH₃OCOH (200 K), and
the deuterated methyl formate lines (60 K), together with the
final two component model. While some species give reason-
able agreement with observations using either model (mostly
those which observed transitions have similar energy levels),
there are others (besides the ones discussed above) that fail
in one of the models (e.g. HNCO and t-CH₃CH₂OCOH). We
have checked that using other temperatures (between 50 and
250 K) in the one component models does not reproduce sat-
isfactorily the observations. For the vibrationally excited states,
we only used the compact component at 200 K, since the low
temperature one does not contribute much to the line inten-
sities. Indeed, energy levels of the observed transitions range
between 300 and 500 K, supporting this scenario. We have also
observed five transitions of ¹³CH₃OCOH, with energies ≤200
K and intensities of the order of 500 mK or below. Since the
lines are weak, any of the models give reasonable agreement
with the observations. Acetaldehyde is the other species for
which many lines are observed. We have detected 15 transitions

(10 unblended) from the main species, with energies between
100 and 300 K. We have also detected one transition arising from
the deuterated substitution CH3CDO, and four transitions (two
of them blended) from the vibrational state v1 = 1, with energy
levels of ~300-400 K. We used the same model as methyl for-
mate, except for CH3CHO v1 = 1, for which only the compact
and hot component was considered. The model allows to repro-
duce all transitions reasonably well, except for a few lines that
are overestimated: two CH3OCOH transitions (Fig. B.1, second
panel), three CH3OCOH v1 = 1 lines (top panel in Fig. B.1, and
second panel in Fig. B.3), and two lines from CH3CHO (see
Fig. B.3, second and third panels). These are transitions with
the lowest energies (~100 K) and the highest line strengths and
Einstein coefficients, and may consequently have high opacities.
We estimate the error in the fit to be ~50%.

This simple two-component model simulates the tempera-
ture gradient that should exist in the source. It is possible that a
more sophisticated model will reproduce the lines better, how-
ever the limited angular resolution and signal-to-noise ratio of
the current observations do not allow a well-constrained model.
In particular, the angular resolution is not sufficient to separate
the emission from the compact and hot component at 200 K from
that of the lukewarm gas at 60 K. For example, it is possible that
the size of the hot component could be smaller, and therefore we
underestimate the column densities.

For the other COMs, we have used the same two-component
model as for methyl formate and acetaldehyde for consistency.
There are two exceptions: the observed lines from ¹³CH₃OH
and HNCO arise from very excited energy levels (~600 and
700 K, respectively). Hence, we have assumed that only the hot
and compact component contributes to the detected emission.
As mentioned above, we have excluded from the analysis all
the species with lines showing complex profiles and extended
emission. We have attempted to run the LTE models for CH₃OH,
H₂CO and their isotopologues, which do not show absorption or
very complex profiles. All the observed lines of H₂CO, H₂¹³CO,
H₂C¹⁷O, and D₂CO (only one transition each) have low E_up
values and show extended emission. Therefore, these lines will
suffer from flux loss due to filtering from the interferometer, and
the obtained column densities are strongly underestimated. For
CH₃OH and CH₂DOH we observed more than one transition
(two and six, respectively), but they span a considerable E_up
<table><caption>**Table 1.** Results from LTE model.</caption><thead><tr><th>Species</th><th>N (cm<sup>-2</sup>)</th><th>X/H<sub>2</sub><sup>a</sup></th></tr></thead><tbody><tr><td>NH<sub>2</sub>CN</td><td>1.0; 1.0 &times; 10<sup>12</sup></td><td>7.1; 9.1 &times; 10<sup>-14</sup></td></tr><tr><td>HNCO</td><td>3.0; &minus; &times; 10<sup>14</sup></td><td>2.1; &minus; &times; 10<sup>-11</sup></td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OH</td><td>5.0; &minus; &times; 10<sup>15</sup></td><td>3.6; &minus; &times; 10<sup>-10</sup></td></tr><tr><td>CH<sub>3</sub>CHO</td><td>8.0; 1.6 &times; 10<sup>14</sup></td><td>5.7; 1.4 &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>3</sub>CHO v<sub>t</sub> = 1</td><td>4.0; &minus; &times; 10<sup>14</sup></td><td>2.8; &minus; &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>3</sub>CDO</td><td>4.0; 4.0 &times; 10<sup>13</sup></td><td>2.8; 3.6 &times; 10<sup>-12</sup></td></tr><tr><td>NH<sub>2</sub>CHO</td><td>4.0; 4.0 &times; 10<sup>12</sup></td><td>2.8; 3.6 &times; 10<sup>-13</sup></td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>1.0; 1.0 &times; 10<sup>16</sup></td><td>7.1; 9.1 &times; 10<sup>-10</sup></td></tr><tr><td>t-CH<sub>3</sub>CH<sub>2</sub>OH</td><td>2.0; 2.0 &times; 10<sup>14</sup></td><td>1.4; 1.8 &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>2</sub>CHCN</td><td>4.0; 2.0 &times; 10<sup>13</sup></td><td>2.8; 1.8 &times; 10<sup>-12</sup></td></tr><tr><td>t-CH<sub>2</sub>CHCHO</td><td>3.0; 2.0 &times; 10<sup>13</sup></td><td>2.1; 1.8 &times; 10<sup>-12</sup></td></tr><tr><td>CH<sub>2</sub>OCHCHO</td><td>5.0; 3.0 &times; 10<sup>14</sup></td><td>3.6; 2.7 &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>3</sub>OCOH</td><td>1.5; 0.5 &times; 10<sup>16</sup></td><td>1.1; 0.4 &times; 10<sup>-9</sup></td></tr><tr><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1</td><td>2.0; &minus; &times; 10<sup>15</sup></td><td>1.4; &minus; &times; 10<sup>-10</sup></td></tr><tr><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 2</td><td>8.0; &minus; &times; 10<sup>14</sup></td><td>5.7; &minus; &times; 10<sup>-11</sup></td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OCOH</td><td>3.0; 3.0 &times; 10<sup>14</sup></td><td>2.1; 2.7 &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>3</sub>OCOD</td><td>3.0; 5.0 &times; 10<sup>14</sup></td><td>2.1; 4.5 &times; 10<sup>-11</sup></td></tr><tr><td>CH<sub>2</sub>DOCOH</td><td>0.4; 1.0 &times; 10<sup>15</sup></td><td>2.8; 9.1 &times; 10<sup>-11</sup></td></tr><tr><td>aGg<sup>-</sup>-(CH<sub>2</sub>OH)<sub>2</sub></td><td>2.0; 0.8 &times; 10<sup>14</sup></td><td>1.4; 0.7 &times; 10<sup>-11</sup></td></tr><tr><td>gGg<sup>-</sup>-(CH<sub>2</sub>OH)<sub>2</sub></td><td>9.0; 6.0 &times; 10<sup>13</sup></td><td>6.4; 5.4 &times; 10<sup>-12</sup></td></tr><tr><td>t-CH<sub>3</sub>CH<sub>2</sub>OCOH</td><td>4.6; 0.8 &times; 10<sup>14</sup></td><td>3.3; 0.7 &times; 10<sup>-11</sup></td></tr></tbody></table>

**Notes.** First values correspond to the compact and hot component (200 K, 0.35"), second values to the extended and cold one (60 K, 0.60"). For vibrationally excited transitions, HNCO and $^{13}$CH$_{3}$OH, only the hot component was considered. (a) Using N(H$_{2}$) = 1.4 × 10$^{25}$ cm$^{-2}$ at 0.35", and N(H$_{2}$) = 1.1 × 10$^{25}$ cm$^{-2}$ at 0.60" (see text).

range, which makes a simultaneous fit difficult. In particular, lines with low energies and extended emission are significantly overestimated. We have observed only one transition from CH$_{3}$OD (also extended) and $^{13}$CH$_{3}$OH. We have kept only the LTE result for the latter which does not present any extended emission. Lines from DCO$^{+}$, DCN, $^{18}$O, SO, $^{13}$CO, CO, $^{13}$CS, and N$_{2}$D$^{+}$, are seen in absorption or show very complex profiles, and we have not attempted to fit them with the LTE model. For all these species, we show the line identification in Figs. 1 and B.1–B.3.

Table 1 shows, for each source component, the obtained column densities and molecular abundances with respect to H$_{2}$. In order to compute the abundances we have estimated the molecular hydrogen column densities for both components. Based on previous observations of dust continuum emission with ALMA at 850 µm and JVLA at 8–10 cm (those with the best angular resolution, see Appendix A) and assuming dust temperatures from 12 to 50 K, we have computed mean column densities of N(H$_{2}$) = 1.4 × 10$^{25}$ cm$^{-2}$ and 1.1 × 10$^{25}$ cm$^{-2}$ for source sizes of 0.35" and 0.6", respectively. While the dust emission is optically thick at 350 GHz on size scales smaller than 0.3", the 8 cm emission is optically thin. Then it is possible the H$_{2}$ column density estimates could be affected by the opacity, in particular at 0.35" more than at 0.6". Another important source of error comes from the unknown dust properties, since the standard opacity law may be different if dust grains have grown. All these factors could result in uncertainties in the column density estimates up to a factor of two.

## 5. Discussion

Previous spectral studies with the IRAM 30 m telescope already pointed out the presence of complex organics in B1b.

Cernicharo et al. (2012) reported the detection of several COMs, including CH$_{3}$CHO, CH$_{3}$OCOH, and CH$_{3}$OCH$_{3}$, also observed with ALMA. These detections are based on low energy transitions, <25 K, and because of the large 30 m beam at 3 mm (30–21") the emission should be arising from the cold envelope. Table B.3 shows the abundances derived from the ALMA data, together with those obtained in the envelope from the IRAM 30 m survey at 3 mm. Except for CH$_{2}$CHCN and HNCO, there is a general trend of COMs being enhanced in the inner core of B1b-S with respect to the colder protostellar envelope, with methyl formate showing the highest increase. This may be related to the evaporation of grain mantles due to the increased dust temperature caused by the heating from the embedded protostar. A similar increase in abundance is seen in other hot corinos with respect to the values observed in their envelopes. In IRAS 16293-2422, Jaber et al. (2014) obtain higher abundances in the inner warm envelope with respect to the outer cold envelope, and these values are lower than those observed at smaller scales using interferometers, in other words, closer to the proto-star (Jørgensen et al. 2012; Coutens et al. 2018). The range of temperatures obtained for B1b-S cores in the LTE model, from 60 K up to at least 200 K in a size of 0.35", and the narrow line-widths are consistent with thermal evaporation. Other mechanisms could also be at work, such as accretion shocks in the disk or launching of the outflow. However, we can not conclude about the origin of COMs since we lack angular resolution. Indeed, it is possible that the size of the compact component in the model (0.35") is smaller, since we do not resolve the emission from the highest energy transitions.

In the ALMA data we have also detected the singly deuterated isotopologues of acetaldehyde and methyl formate, whose abundances decrease at the inner hot component with respect to the component at 60 K, with D/H a factor of five higher in the latter (see Table B.3). These deuterated species are not detected with the IRAM 30 m, possibly because the main species are weak and less abundant in the envelope. However, high deuteration levels are observed for other species, for example, strong HDCO lines were detected by Marcelino et al. (2005) with a D to H ratio of 0.145. This value is similar to those derived in the external B1b-S core for CH$_{3}$OCOH and CH$_{2}$DOCOH (D/H of 0.1 and 0.2, respectively), and smaller than CH$_{3}$CDO/CH$_{3}$CHO = 0.25. High deuterium fractionation is also observed in hot corinos, while hot cores show lower deuterium abundances (see Caselli & Ceccarelli 2012, and references within).

The comparison with other hot corinos is not easy since we have a limited number of lines, in contrast to well-known sources that have been more extensively observed with ALMA and NOEMA, such as IRAS 16293-2422 (Jørgensen et al. 2016), NGC1333 IRAS 2A and 4A (Taquet et al. 2015; López-Sepulcre et al. 2017). In general, we find a similar inventory of COMs but with lower abundances in B1b-S. We are going to discuss only some trends and ratios between them. For example, one similarity with hot corinos are the higher abundances of O-bearing COMs with respect to N-bearing COMs. We observe a similar trend in B1b-S, where the highest abundances are obtained for CH$_{3}$OCOH, CH$_{3}$OCH$_{3}$, and CH$_{3}$CHO, while NH$_{2}$CN and NH$_{2}$CHO show the lowest abundances. Cyanamide has been recently detected towards the low-mass protostars IRAS 16293-2422 and NGC 1333 IRAS2A (Coutens et al. 2018), with observed NH$_{2}$CN to NH$_{2}$CHO ratios of 0.2 and 0.02, respectively. These values are in the range of those observed towards the molecular clouds in the Galactic centre but lower than in Orion KL (Coutens et al. 2018). We obtain in B1b-S an abundance ratio of 0.25, similar to IRAS 16293-2422. Of the three
possible isomers of C₂H₄O₂, we have detected CH₃OCOH and CH₂OHCHO. The observed CH₃OCOH to CH₂OHCHO ratio in B1b-S, ~20, is similar to that observed in the low-mass protostars in NGC 1333 and IRAS 16293-2422 (Taquet et al. 2015; Jørgensen et al. 2012). Acetic acid (CH₃COOH) is the most stable but the least abundant of the three isomers (Lattelais et al. 2010). It has been observed in IRAS 16293-2422 with a ratio with respect to glycolaldehyde of ~11 (Jørgensen et al. 2016), consistent with the 5–15 upper limits in B1b-S. Glycolaldehyde and its corresponding alcohol, ethylene glycol, show similar abundances in B1b-S, slightly higher for CH₂OHCHO. This is in contrast to other hot cores and hot corinos (Fuente et al. 2014; Jørgensen et al. 2016; Favre et al. 2017). However, since detections of these species are based on just a few lines, in particular for glycolaldehyde, column densities may be not well constrained and it is difficult to make definite conclusions. We have detected CH₃CH₂OCOH for the first time in a low-mass protostar, with a relatively high abundance of 10⁻¹¹. This detection shows that the molecular complexity is high in young hot corinos and other species with a larger number of atoms could be present. It is noteworthy that this species has not been previously observed in the well-known hot corinos IRAS 16293-2422 and NGC 133IRAS4A. We have checked in the ALMA spectral line survey of IRAS 16293-2422 (PILS, Jørgensen et al. 2016), and no clear features are seen, although the observed frequencies are different. It is possible that the high level of line confusion and wide linewidths prevent the detection of weak lines in these more evolved sources.

Observed line widths in B1b-S are typically narrow, ~1 km s⁻¹ (see Table B.2), in contrast to the wide line profiles in hot cores and hot corinos (~3–5 km s⁻¹). Besides all lines from COMs are found in emission, with no P Cygni or absorption profiles. Therefore, there is no signature of rotation or infall in the disk or core from COMs. The orientation of the disk should not prevent the detection of rotation, following the observations of the outflows in B1b-S (Gerin et al. 2015). Actually, using NOEMA data, Fuente et al. (2017) detect rotation using NH₂D. The absence of detected rotation could be due to the too low angular resolution since the COM emission is unresolved. In addition, although the sensitivity is excellent, these observations are not deep enough to detect weak features at the 0.3–0.6 K level (see rms in Table B.1). Hence, weak line wings with rotation signatures may have been missed.

## 6. Conclusions

We present new ALMA Band 6 data towards the protostellar system B1b in Perseus. We have detected a rich and complex spectra in B1b-S, with 190 spectral features. Of those, 99 lines are assigned to transitions from complex molecules, while 77 remain unidentified. Previous NOEMA observations do not show any emission from COMs, with spatial resolutions between 2'' and 5'' (Gerin et al. 2015; Fuente et al. 2017). These data suggest that embedded in the B1b-S core a very young and compact object is already warming up its most immediate surroundings, where molecules from grain mantles are evaporating. On the other hand, B1b-N is almost free of line emission. The non-detection of complex species supports the scenario in which B1b-N is at an earlier stage of evolution than B1b-S. Another possibility is that the more active outflow in B1b-S, with higher velocities and more complex profiles (with evidences of bullets), could contribute to a more efficient release of COM in the gas phase in B1b-S, as compared to B1b-N. Luminosity outbursts due to episodic accretion in B1b-S could also increase the dust

temperature and produce the evaporation of COMs from grain mantles or trigger the formation of daughter molecules in the gas phase (Taquet et al. 2016). The origin of COMs in B1b-S remains unclear since the emission is not resolved.

While B1b-S presents a rich spectra similar to other hot corinos, observed line widths are typically narrower and do not show any signature of rotation or infall. This could be due to the limited angular resolution and signal-to-noise ratio. These observations confirm that hot corinos are a common feature in class 0 protostars, but require sub-arcsec angular resolution and high sensitivity reaching a few mJy beam⁻¹ to be detected. Higher angular resolution observations of B1b-S with better sensitivity are desirable to resolve the emission from COMs, and better constrain molecular abundances, the size and the structure of the emitting region.

*Acknowledgements.* We thank the referee for a careful reading of the manuscript and useful comments which have improved the paper significantly. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2015.1.00025.S. ALMA is a partnership of ESO (representing its member states), NSF (USA) and NINS (Japan), together with NRC (Canada), MOST and ASIAA (Taiwan), and KASI (Republic of Korea), in cooperation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI/NRAO and NAOJ. We acknowledge funding support from the European Research Council (ERC Grant 610256: NANOCOSMOS), the ANR project IMOLABS, and the Spanish MINECO under projects AYA2016-75066-C2-2-P and AYA2012-32032. This work was partly supported by the Programme National "Physique et Chimie du Milieu Interstellaire" (PCMI) of CNRS/INSU with INC/INP, co-funded by CEA and CNES. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc.

## References

Belloche, A., Garrod, R. T., Müller, H. S. P., et al. 2009, A&A, 499, 215

Caselli, P., & Ceccarelli, C. 2012, A&ARv, 20, 56

Cernicharo, J. 2012, in ECLA-2011: Proc. of the European Conference on Laboratory Astrophysics, EAS Pub. Ser., eds. C. Stehlé, C. Joblin, & L. d'Hendecourt (Cambridge: Cambridge Univ. Press)

Cernicharo, J., Marcelino, N., Roueff, E., et al. 2012, ApJ, 759, L43

Commerçon, B., Levrier, F., Maury, A. J., Henning, Th., & Launhardt, R. 2012, A&A, 548, A39

Coutens, A., Willis, E. R., Garrod, R. T., et al. 2018, A&A, 612, A107

Cox, E. G., Harris, R. J., Looney, L. W., et al. 2018, ApJ, 855, 92

Daniel, F., Gerin, M., Roueff, E., et al. 2013, A&A, 560, A3

Evans, N. J., II, Dunham, M. M., Jorgensen, J. K., et al. 2009, ApJS, 181, 321

Favre, C., Pagani, L., Goldsmith, P., et al. 2017, A&A, 604, L2

Fuente, A., Cernicharo, J., Caselli, P., et al. 2014, A&A, 568, A65

Fuente, A., Gerin, M., Pety, J., et al. 2017, A&A, 606, L3

Gerin, M., Pety, J., Fuente, A., et al. 2015, A&A, 577, L2

Gerin, M., Pety, J., Commerçon, B., et al. 2017, A&A, 606, A35

Hatchell, J., Richer, J. S., Fuller, G. A., et al. 2005, A&A, 440, 151

Hatchell, J., Fuller, G. A., & Richer, J. S. 2007, A&A, 472, 187

Herbst, E., & van Dishoeck, E. F. 2009, ARA&A, 47, 427

Hincelin, U., Commerçon, B., Wakelam, V., et al. 2016, ApJ, 822, 12

Hirano, N., & Liu, F. -C. 2014, ApJ, 789, 50

Hirota, T., Bushimata, T., Choi, Y. K., et al. 2008, PASJ, 60, 37

Jaber, A. A., Ceccarelli, C., Kahane, C., & Caux, E. 2014, ApJ, 791, 29

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

Jørgensen, J. K., van der Wiel, M. H. D., Coutens, A., et al. 2016, A&A, 595, A117

Larson, R. B. 1969, MNRAS, 145, 271

Lattelais, M., Pauzat, F., Ellinger, Y., & Ceccarelli, C. 2010, A&A, 519, A30

López-Sepulcre, A., Sakai, N., Neri, R., et al. 2017, A&A, 606, A121

Marcelino, N., Cernicharo, J., Roueff, E., et al. 2005, ApJ, 620, 308

Marcelino, N., Cernicharo, J., Tercero, B., & Roueff, E. 2009, ApJ, 690, L27

Pezzuto, S., Elia, D., Schisano, E., et al. 2012, A&A, 547, A54

Pokhrel, R., Myers, P. C., Dunham, M. M., et al. 2018, ApJ, 853, 5

Rivilla, V. M., Beltrán, M. T., Martín-Pintado, J., et al. 2017, A&A, 599, A26

Taquet, V., López-Sepulcre, A., Ceccarelli, C., et al. 2015, ApJ, 804, 81

Taquet, V., Wirst्रöm, E. S., & Charnley, S. B. 2016, ApJ, 821, 46

Tercero, B., Kleiner, I., Cernicharo, J., et al. 2013, ApJ, 770, L13

Tobin, J. J., Looney, L. W., Li, Z.-Y., et al. 2016, ApJ, 818, 73
# Appendix A: Continuum results

The frequency setup included a broad spectral window, in order to detect the continuum emission of the B1b sources, and to facilitate the bandpass and phase calibration on the reference quasars. We present in the left panel in Fig. 1 the resulting continuum image. The three Barnard 1b protostars are detected, B1b-N, B1b-S and even the weak source B1b-W. Table A.1 provides a compilation of the spectral energy distribution for these three sources. Given the extended nature of the sources the listed values depend on the angular resolution of the observations. It is therefore important to compare data measured at the same angular resolution.

<table><caption>**Table A.1.** Observed fluxes at different wavelengths for the B1b sources.</caption><thead><tr><th>λ<br>(μm)</th><th>B1b-N<br>(Jy)</th><th>B1b-S<br>(Jy)</th><th>B1b-W<br>(Jy)</th><th>Noise<br>(Jy)</th><th>Beam</th><th>Refs.</th></tr></thead><tbody><tr><td>24</td><td>&lt;2 × 10<sup>-4</sup></td><td>&lt;2 × 10<sup>-4</sup></td><td>0.17</td><td>0.02</td><td>12</td><td>1, 2</td></tr><tr><td>70</td><td>&lt;0.05</td><td>0.22</td><td>0.36</td><td>0.15</td><td>7</td><td>1, 2, 3</td></tr><tr><td>100</td><td>0.6</td><td>2.29</td><td>1.0</td><td>0.3</td><td>8.5</td><td>2</td></tr><tr><td>160</td><td>3.2</td><td>9.1</td><td>&lt;3</td><td>1</td><td>8.9</td><td>2</td></tr><tr><td>250</td><td>9.5</td><td>14.4</td><td>&lt;5</td><td>1</td><td>14.3</td><td>2</td></tr><tr><td>350</td><td>12.7</td><td>16.9</td><td>&lt;10</td><td>3</td><td>18.4</td><td>2</td></tr><tr><td>350</td><td>6.0</td><td>7.2</td><td>&lt;10</td><td>0.3</td><td>20</td><td>3</td></tr><tr><td>450</td><td>13.5</td><td>19.1</td><td>&lt;6</td><td>4</td><td>11</td><td>2</td></tr><tr><td>500</td><td>14.6</td><td>15.8</td><td>&lt;6</td><td>7</td><td>34.3</td><td>2</td></tr><tr><td>850</td><td>3.3</td><td>3.8</td><td>&lt;3</td><td>1.3</td><td>15</td><td>2</td></tr><tr><td>850</td><td>1.03</td><td>1.24</td><td>&lt;1</td><td>0.03</td><td>15</td><td>3</td></tr><tr><td>850</td><td>0.24</td><td>0.48</td><td>&lt;1</td><td>0.01</td><td>0.3</td><td>4</td></tr><tr><td>850</td><td>0.73</td><td>0.89</td><td>&lt;1</td><td>0.01</td><td>1</td><td>4</td></tr><tr><td>870</td><td>—</td><td>—</td><td>0.0194</td><td>1.4 × 10<sup>-4</sup></td><td>0.3</td><td>5</td></tr><tr><td>1057</td><td>0.31</td><td>0.47</td><td>&lt;0.3</td><td>0.015</td><td>3</td><td>3</td></tr><tr><td>1290</td><td>0.12</td><td>0.23</td><td>6.4 × 10<sup>-3</sup></td><td>1 × 10<sup>-3</sup></td><td>0.6</td><td>6</td></tr><tr><td>1300</td><td>0.19</td><td>0.34</td><td>&lt;0.2</td><td>0.01</td><td>5</td><td>3</td></tr><tr><td>1300</td><td>0.249</td><td>0.375</td><td>—</td><td>0.011</td><td>4</td><td>7</td></tr><tr><td>2070</td><td>0.077</td><td>0.096</td><td>&lt;0.04</td><td>4 × 10<sup>-3</sup></td><td>2.2</td><td>8</td></tr><tr><td>3300</td><td>0.026</td><td>0.031</td><td>&lt;0.02</td><td>3.5 × 10<sup>-3</sup></td><td>3</td><td>3</td></tr><tr><td>3660</td><td>0.018</td><td>0.016</td><td>&lt;1 × 10<sup>-3</sup></td><td>1 × 10<sup>-3</sup></td><td>3.9</td><td>8</td></tr><tr><td>7000</td><td>1.9 × 10<sup>-3</sup></td><td>1.7 × 10<sup>-3</sup></td><td>&lt;2 × 10<sup>-3</sup></td><td>2.8 × 10<sup>-4</sup></td><td>0.5</td><td>3</td></tr><tr><td>8000</td><td>1.25 × 10<sup>-3</sup></td><td>0.91 × 10<sup>-3</sup></td><td>1 × 10<sup>-4</sup></td><td>4 × 10<sup>-5</sup></td><td>0.2</td><td>9</td></tr><tr><td>10000</td><td>6.8 × 10<sup>-4</sup></td><td>5.5 × 10<sup>-4</sup></td><td>8 × 10<sup>-5</sup></td><td>3 × 10<sup>-5</sup></td><td>0.2</td><td>9</td></tr></tbody></table>

**References.** (1) Evans et al. (2009); (2) Pezzuto et al. (2012); (3) Hirano & Liu (2014); (4) Gerin et al. (2017); (5) Cox et al. (2018); (6) This work; (7) Pokhrel et al. (2018); (8) Gerin et al. (2015); (9) Tobin et al. (2016).

# Appendix B: Additional figures and tables

<table><caption>**Table B.1.** Summary of the observed spectral line windows.</caption><thead><tr><th>Frequency range<br>(MHz)</th><th>Beam<br>(")</th><th>PA<sup>a</sup><br>(°)</th><th>σF<sup>b</sup><br>(mJy beam<sup>-1</sup>)</th><th>rms<sup>c</sup><br>(mK)</th></tr></thead><tbody><tr><td>215811–215869</td><td>0.74 × 0.57</td><td>23</td><td>3.29</td><td>128</td></tr><tr><td>216083–216142</td><td>0.74 × 0.57</td><td>23</td><td>3.64</td><td>136</td></tr><tr><td>217076–217134</td><td>0.73 × 0.57</td><td>22</td><td>3.52</td><td>184</td></tr><tr><td>217209–217268</td><td>0.73 × 0.58</td><td>16</td><td>3.78</td><td>153</td></tr><tr><td>218432–218490</td><td>0.74 × 0.56</td><td>26</td><td>3.71</td><td>216</td></tr><tr><td>219531–219590</td><td>0.70 × 0.57</td><td>22</td><td>3.86</td><td>196</td></tr><tr><td>219899–219957</td><td>0.69 × 0.57</td><td>19</td><td>4.14</td><td>169</td></tr><tr><td>220367–220426</td><td>0.69 × 0.57</td><td>18</td><td>4.76</td><td>250</td></tr><tr><td>230519–230577</td><td>0.67 × 0.54</td><td>17</td><td>4.46</td><td>193</td></tr><tr><td>231192–231250</td><td>0.67 × 0.53</td><td>15</td><td>4.57</td><td>194</td></tr><tr><td>231293–231351</td><td>0.59 × 0.44</td><td>16</td><td>6.17</td><td>414</td></tr><tr><td>231381–231439</td><td>0.68 × 0.53</td><td>16</td><td>4.04</td><td>184</td></tr></tbody></table>

**Notes.** (a) Beam position angle. (b) rms noise level in the images. (c) rms noise level measured in the B1b-S spectra.
<figure><img src="image_3.png" /><figcaption>**Fig. B.1.** Observed spectra towards B1b-S (black) and B1b-N (blue), between 215810 and 217268 MHz. Line features with no labels are unidentified lines. The LTE simulations for B1b-S are over-plotted in red.</figcaption></figure>
<figure><img src="image_4.png" /><figcaption>**Fig. B.2.** Observed spectra towards B1b-S (black) and B1b-N (blue), between 218432 and 220426 MHz. Line features with no labels are unidentified lines. The LTE simulations for B1b-S are over-plotted in red.</figcaption></figure>
<figure><img src="image_5.png" /><figcaption>**Fig. B.3.** Observed spectra towards B1b-S (black) and B1b-N (blue), between 230519 and 231440 MHz. Line features with no labels are unidentified lines. The LTE simulations for B1b-S are over-plotted in red.</figcaption></figure>
<figure><img src="image_6.png" /><figcaption>**Fig. B.4.** Best LTE models using only one source component for two spectral windows. *Middle panels* for each frequency range show the results for the 200 K model in red, while *bottom panels* show in blue the 60 K model. *Top panels* in purple show the final two-component model.</figcaption></figure>
<table><caption>**Table B.2.** Observed transitions in B1b-S.</caption><thead><tr><th>Obs. Freq. (MHz)</th><th>Molecule and transition</th><th>Rest freq. (MHz)</th><th>Eup (K)</th><th>S<sub>ij</sub></th><th>A<sub>ij</sub> (s<sup>-1</sup>)</th><th>∫ T<sub>Baν</sub> (K km s<sup>-1</sup>)</th><th>V<sub>LSR</sub> (km s<sup>-1</sup>)</th><th>Δ<sub>V</sub> (km s<sup>-1</sup>)</th><th>T<sub>B</sub> (K)</th><th>Notes</th></tr></thead><tbody><tr><td>215814.58</td><td>t-CH<sub>2</sub>CHCHO 25<sub>1,25</sub> – 24<sub>1,24</sub></td><td>215814.794(4)</td><td>137.1</td><td>24.95</td><td>5.33 × 10<sup>-4</sup></td><td>0.442(94)</td><td>6.916(88)</td><td>0.841(207)</td><td>0.494(98)</td><td></td></tr><tr><td>215822.35</td><td><sup>13</sup>CH<sub>3</sub>OCOH 18<sub>14,5</sub> – 17<sub>14,4</sub> E</td><td>215822.298(13)</td><td>227.8</td><td>7.08</td><td>5.95 × 10<sup>-5</sup></td><td>0.149(45)</td><td>6.517(77)</td><td>0.477(151)</td><td>0.294(83)</td><td></td></tr><tr><td>215829.57</td><td>Unidentified</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.286(52)</td><td>6.626(55)</td><td>0.533(103)</td><td>0.505(110)</td><td></td></tr><tr><td>215832.97</td><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 2, 19<sub>1,18</sub> – 18<sub>1,17</sub> A</td><td>215833.102(3)</td><td>441.4</td><td>52.20</td><td>1.57 × 10<sup>-4</sup></td><td>0.696(21)</td><td>6.789(9)</td><td>0.627(22)</td><td>1.044(86)</td><td></td></tr><tr><td>215835.13</td><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1, 20<sub>0,20</sub> – 19<sub>1,19</sub> A</td><td>215835.178(100)</td><td>299.3</td><td>18.42</td><td>2.43 × 10<sup>-5</sup></td><td>0.769(19)</td><td>6.675(10)</td><td>0.859(27)</td><td>0.841(107)</td><td></td></tr><tr><td>215837.52</td><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1, 20<sub>1,20</sub> – 19<sub>1,19</sub> A</td><td>215837.587(100)</td><td>299.3</td><td>19.84</td><td>1.50 × 10<sup>-4</sup></td><td>1.470(192)</td><td>6.683(53)</td><td>0.847(136)</td><td>1.631(61)</td><td></td></tr><tr><td>215839.38</td><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1, 20<sub>0,20</sub> – 19<sub>0,19</sub> A</td><td>215839.544(100)</td><td>299.3</td><td>19.84</td><td>1.50 × 10<sup>-4</sup></td><td>1.088(241)</td><td>6.847(74)</td><td>0.675(169)</td><td>1.513(49)</td><td></td></tr><tr><td>215841.90</td><td>CH<sub>3</sub>OCOD 18<sub>-10,8</sub> – 17<sub>-10,7</sub> E</td><td>215841.936(6)</td><td>158.1</td><td>12.45</td><td>1.06 × 10<sup>-4</sup></td><td>...</td><td>...</td><td>...</td><td>...</td><td>Blen</td></tr><tr><td>215843.47</td><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1, 20<sub>1,20</sub> – 19<sub>0,19</sub> A</td><td>215842.029(100)</td><td>299.3</td><td>18.42</td><td>2.43 × 10<sup>-5</sup></td><td>...</td><td>...</td><td>...</td><td>...</td><td>Blen</td></tr><tr><td>215848.10</td><td>CH<sub>2</sub>DOH 11<sub>6,6</sub> – 12<sub>5,8</sub> e<sub>1</sub></td><td>215843.384(9)</td><td>293.2</td><td>1.17</td><td>3.83 × 10<sup>-6</sup></td><td>1.143(209)</td><td>6.413(83)</td><td>0.963(219)</td><td>1.114(110)</td><td></td></tr><tr><td>215849.09</td><td>CH<sub>2</sub>DOH 11<sub>6,5</sub> – 12<sub>5,7</sub> e<sub>1</sub></td><td>215843.397(9)</td><td>293.2</td><td>1.17</td><td>3.83 × 10<sup>-6</sup></td><td>...</td><td>...</td><td>...</td><td>...</td><td></td></tr><tr><td>215850.33</td><td>Unidentified</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.733(90)</td><td>6.602(63)</td><td>1.061(153)</td><td>0.649(129)</td><td>Ori</td></tr><tr><td>215853.70</td><td>Unidentified</td><td>...</td><td>...</td><td>...</td><td>...</td><td>0.259(105)</td><td>6.598(125)</td><td>0.605(278)</td><td>0.403(123)</td><td>Ori</td></tr><tr><td>215858.22</td><td>CH<sub>3</sub>OCOD 18<sub>10,9</sub> – 17<sub>10,8</sub> A</td><td>215853.832(6)</td><td>158.2</td><td>12.46</td><td>1.06 × 10<sup>-4</sup></td><td>0.558(173)</td><td>6.699(111)</td><td>0.708(263)</td><td>0.740(105)</td><td></td></tr><tr><td>215861.71</td><td>CH<sub>3</sub>OCOD v<sub>t</sub> = 1, 11<sub>3,9</sub> – 11<sub>2,10</sub> A</td><td>215858.298(27)</td><td>285.7</td><td>3.89</td><td>3.15 × 10<sup>-5</sup></td><td>0.863(63)</td><td>6.738(28)</td><td>0.785(63)</td><td>1.034(108)</td><td></td></tr><tr><td>215863.29</td><td>Unidentified</td><td>...</td><td>158.1</td><td>12.45</td><td>1.06 × 10<sup>-4</sup></td><td>0.316(201)</td><td>6.685(195)</td><td>0.618(486)</td><td>0.480(153)</td><td></td></tr><tr><td>215864.38</td><td>CH<sub>2</sub>OHCHO 33<sub>6,28</sub> – 33<sub>5,29</sub></td><td>215864.490(10)</td><td>335.8</td><td>14.66</td><td>1.39 × 10<sup>-4</sup></td><td>0.408(167)</td><td>6.600(149)</td><td>0.803(410)</td><td>0.478(140)</td><td></td></tr><tr><td>215865.68</td><td>CH<sub>3</sub>CHO v<sub>t</sub> = 1, 11<sub>1,10</sub> – 10<sub>1,9</sub> A</td><td>215865.742(7)</td><td>271.0</td><td>11.85</td><td>3.54 × 10<sup>-4</sup></td><td>1.886(210)</td><td>6.653(50)</td><td>1.180(173)</td><td>1.501(170)</td><td></td></tr></tbody></table>

**Notes.** Line parameters are obtained from Gaussian fits. Rest frequencies and spectroscopic data are taken from the MADEX catalogue (Cernicharo 2012). Abs: line observed in absorption. Blen: line blended with others. CP: line with complex line profile. Edge: line at the edge of the spectral window. Ext: line with extended emission. Ori: unidentified line also detected in Orion-KL. The full Table is available in electronic form at the CDS.
<table><caption>**Table B.3.** Obtained abundances with respect to H₂ in B1b-S (this work), and towards the B1b envelope observed with the IRAM 30 m. The D/H obtained for acetaldehyde and methyl formate for the hot and cold components in B1b-S are also shown.</caption><thead><tr><th rowspan="2">Species</th><th colspan="3"><i>X</i>/<i>H</i><sub>2</sub><sup>a</sup></th><th colspan="2"><i>D</i>/<i>H</i></th></tr><tr><th>Hot (200 K)</th><th>Cold (60 K)</th><th>Env. (10 K)</th><th>Hot (200 K)</th><th>Cold (60 K)</th></tr></thead><tbody><tr><td>NH<sub>2</sub>CN</td><td>7.1 × 10<sup>-14</sup></td><td>9.1 × 10<sup>-14</sup></td><td></td><td></td><td></td></tr><tr><td>HNCO</td><td>2.1 × 10<sup>-11</sup></td><td>—</td><td>1.2 × 10<sup>-10</sup></td><td></td><td></td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OH</td><td>3.6 × 10<sup>-10</sup></td><td>—</td><td>8.6 × 10<sup>-11</sup></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>CHO</td><td>5.7 × 10<sup>-11</sup></td><td>1.4 × 10<sup>-11</sup></td><td>2.0 × 10<sup>-11</sup></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>CHO v<sub>t</sub> = 1</td><td>2.8 × 10<sup>-11</sup></td><td>—</td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>CDO</td><td>2.8 × 10<sup>-12</sup></td><td>3.6 × 10<sup>-12</sup></td><td></td><td>0.05</td><td>0.25</td></tr><tr><td>NH<sub>2</sub>CHO</td><td>2.8 × 10<sup>-13</sup></td><td>3.6 × 10<sup>-13</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>OCH<sub>3</sub></td><td>7.1 × 10<sup>-10</sup></td><td>9.1 × 10<sup>-10</sup></td><td>3.9 × 10<sup>-11</sup></td><td></td><td></td></tr><tr><td><i>t</i>-CH<sub>3</sub>CH<sub>2</sub>OH</td><td>1.4 × 10<sup>-11</sup></td><td>1.8 × 10<sup>-11</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>2</sub>CHCN</td><td>2.8 × 10<sup>-12</sup></td><td>1.8 × 10<sup>-12</sup></td><td>8.2 × 10<sup>-12</sup></td><td></td><td></td></tr><tr><td><i>t</i>-CH<sub>2</sub>CHCHO</td><td>2.1 × 10<sup>-12</sup></td><td>1.8 × 10<sup>-12</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>2</sub>OHCHO</td><td>3.6 × 10<sup>-11</sup></td><td>2.7 × 10<sup>-11</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>OCOH</td><td>1.1 × 10<sup>-9</sup></td><td>0.4 × 10<sup>-9</sup></td><td>3.9 × 10<sup>-11</sup></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 1</td><td>1.4 × 10<sup>-10</sup></td><td>—</td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>OCOH v<sub>t</sub> = 2</td><td>5.7 × 10<sup>-11</sup></td><td>—</td><td></td><td></td><td></td></tr><tr><td><sup>13</sup>CH<sub>3</sub>OCOH</td><td>2.1 × 10<sup>-11</sup></td><td>2.7 × 10<sup>-11</sup></td><td></td><td></td><td></td></tr><tr><td>CH<sub>3</sub>OCOD</td><td>2.1 × 10<sup>-11</sup></td><td>4.5 × 10<sup>-11</sup></td><td></td><td>0.02</td><td>0.10</td></tr><tr><td>CH<sub>2</sub>DOCOH</td><td>2.8 × 10<sup>-11</sup></td><td>9.1 × 10<sup>-11</sup></td><td></td><td>0.03</td><td>0.20</td></tr><tr><td>aGg'-(-CH<sub>2</sub>OH)<sub>2</sub></td><td>1.4 × 10<sup>-11</sup></td><td>0.7 × 10<sup>-11</sup></td><td></td><td></td><td></td></tr><tr><td>gGg'-(-CH<sub>2</sub>OH)<sub>2</sub></td><td>6.4 × 10<sup>-12</sup></td><td>5.4 × 10<sup>-12</sup></td><td></td><td></td><td></td></tr><tr><td><i>t</i>-CH<sub>3</sub>CH<sub>2</sub>OCOH</td><td>3.3 × 10<sup>-11</sup></td><td>0.7 × 10<sup>-11</sup></td><td></td><td></td><td></td></tr></tbody></table>

**Notes.** The second and fifth columns correspond to the compact and hot component (200 K, 0.35"), and the third and sixth columns to the extended and cold one (60 K, 0.60"). For vibrational excited transitions and HNCO, only the hot component was considered. Abundances for the envelope (fourth column) are from Marcelino et al. (2009), Cernicharo et al. (2012), and Marcelino et al. (in prep.). Blank spaces mean the molecule was not detected in the IRAM 30 m observations. <sup>(a)</sup>Using N(H₂) = 1.4 × 10²⁵ cm⁻² at 0.35", N(H₂) = 1.1 × 10²⁵ cm⁻² at 0.60" (see text), and N(H₂) = 7.6 × 10²² cm⁻² for the envelope (Daniel et al. 2013).