FAUST XI: Enhancement of the complex organic material in the shocked matter surrounding the [BHB2007] 11 protobinary system

Interstellar complex organic molecules (also called iCOMs) are defined in the astronomical community as carbon-bearing species containing at least six atoms (Herbst & van Dishoeck, 2009; Ceccarelli et al., 2017). The detection of iCOMs in solar-type protostars and their precursors in the prestellar core phase is crucial for understanding how molecular complexity emerges and evolves in the different stages that lead to the formation of a solar system such our own. The protostellar phase starts with the young accreting protostar phase with a central object at T $\geq$ 100 K that eventually will become a star, surrounded by a colder envelope from which the future star accretes matter (also called Class 0). In a simple scenario, the molecules trapped in the ices are injected back into the gas phase due to thermal desorption, giving rise to a rich and peculiar chemistry in which iCOMS are found in these regions, called hot corinos (Cazaux et al., 2003; Bottinelli et al., 2004; Ceccarelli, 2004). The name comes from their high-mass analogues, which are called hot cores (Morris et al., 1980). Many iCOMS have been detected in these objects, such as methyl formate (CH${}_{3}$OCHO), dimethyl ether (CH${}_{3}$OCH${}_{3}$), acetaldehyde (CH${}_{3}$CHO), ethanol (CH${}_{3}$CH${}_{2}$OH), and formamide (NH${}_{2}$CHO) (see, e.g., the recent review by Ceccarelli et al., 2023).
In this context, the Atacama Large Millimeter/submillimeter Array (ALMA) Large Program¹¹1http://faust-alma.riken.jp Fifty AU STudy of the chemistry in the disk and envelope system of Solar-like protostars: 2018.1.01205.L (FAUST) is designed to survey the chemical composition of a sample of 13 Class 0/I protostars at a scale of about 50 to 1000 au (all sources have a distance smaller than 250 pc). The description of the project is reported in Codella et al. (2021). One of these sources is the [BHB2007] 11 protobinary system, which is embedded in the darkest parts of the Pipe nebula. This is a quiescent complex of interstellar clouds located at a distance of 163 $\pm$ 5 pc (Dzib et al., 2018) from the Sun. Only a handful of embedded sources are known in the entire complex in the B59 molecular cloud, with [BHB2007] 11 being the youngest member of a protocluster of low-mass young stellar objects (Brooke et al., 2007). At a higher angular resolution ($\sim$ 0.04${}^{\prime\prime}$), Alves et al. (2019) previously uncovered a young binary protostar system embedded in circumstellar disks with radii of 2 to 3 au. The protostars are named [BHB2007] 11 A (hereafter 11A) for the northern source and [BHB2007] 11 B (hereafter 11B) for the southern source. Both cores are surrounded by a complex filamentary structure (within $\sim$ 100${}^{\prime\prime}$) connecting to the larger circumbinary disk. Alves et al. (2019) suggested that these filaments are inflow streamers (also called accretion streamers) from the extended circumbinary disk onto the circumstellar disks of the protobinary system. Their CO observations primarily trace the outflow emission, but the high-velocity components indicate acceleration of infalling gas from the circumbinary disk onto [BHB2007] 11 B. No high-velocity components are detected near [BHB2007] 11 A. Evans et al. (2023) recently presented large-scale formaldehyde emission that seems to be consistent with an asymmetric molecular outflow launched by a wide-angle disk wind. Vastel et al. (2022) reported the emission of a compact hot methanol emission with three velocity components at -2, 2.8, and 9.9 km s${}^{-1}$ and surprisingly found that the methanol components are not associated with the dust emission peaks, which would suggest a hot-corino origin. An estimation of the dust opacities using the 225 GHz dust continuum emission, which could lead to a high attenuation of the methanol emission at the dust peak, had uncertainties that were too large for any conclusive evidence of whether the foreground dust could completely attenuate the methanol emission from the hot corinos. Another possibility would be that methanol emission originates from the shocked positions at the intersection of intertwined filaments, leading to a different origin than the classical hot-corino origin. This could result in the sputtering of the grain mantles in these shocked regions due to the streaming of the gas towards 11A and 11B, which impacts the quiescent gas of the circumbinary envelope or, possibly, the two circumstellar disks (Vastel et al., 2022). Evans et al. (2023) also presented some HDCO and D${}_{2}$CO emission that was only found at the same location of the methanol compact emission, deriving an average D/H ratio of 0.02${}^{+0.02}_{-0.01}$.
In the present work, we investigate the iCOM emission found in the FAUST observations towards [BHB2007] 11 as well as some precursors. Our goal is to understand the origin of the detected iCOM emission by locating and comparing it with the most simple iCOM, methanol, which was previously studied by (Vastel et al., 2022). We present in Section 2 the observations and describe in Section 3 the method we used for the line identification of the iCOMs and precursors. Section 4 presents the spatial distribution of isolated transitions in order to determine the location of their emission. To compute the column densities of the species of interest, we perform a radiative transfer analysis of the emission lines in Section 5. A discussion is finally developed in Section 6 to understand the origin of the emission of the iCOMs in [BHB2007] 11 and compare with other protostellar sources.

The observations were centred at $\alpha$ (2000) = 17h11m23.125s, $\delta$ (2000) = -27${}^{\circ}$24${}^{\prime}$32.87${}^{\prime\prime}$. They consisted of three frequency setups and 13 spectral windows (spw), as explained in previous FAUST papers on [BHB2007] 11 (Vastel et al., 2022; Evans et al., 2023). Table 1 presents the spw for the FAUST [BHB2007] 11 observations, specifically, those in which detections are presented in this study (in bold), as well as those used to compute the upper limits. All the $spw$ were used to search for transitions of the species of interest. We used the 12m array (TM 1 and TM 2) due to the compact emission of the iCOMs transitions. The bandwidth is 59 MHz for the 12m array, except for the continuum window, with a 1.875 GHz bandwidth. Spectral line imaging was performed with the Common Astronomy Software Applications package (CASA) 5.6.1-8 ²²2https://casa.nrao.edu/ with a modified version of the ALMA calibration pipeline and additional in-house calibration routine³³3https://faust-imaging.readthedocs.io/en/latest/ (Chandler in prep.). The systemic velocity of [BHB2007] 11 is 3.6 km s${}^{-1}$ with respect to the local standard of rest (Onishi et al., 1999). For the following line identification, we extracted two spectra from the spw continuum bands to cover the broad velocity emission, which extends within 0.82${}^{\prime\prime}$ $\times$ 0.77${}^{\prime\prime}$ centred at $\alpha$ (2000) = 17h11m23.107s, $\delta$ (2000) = -27${}^{\circ}$24${}^{\prime}$33.005${}^{\prime\prime}$. The two spectra are presented in Fig 1, and the resulting beams and spectral resolutions are quoted in Table 4. Based on the high spatial resolution observations of the dust emission at 1.3 mm with ALMA by Alves et al. (2019), the coordinates for the two cores are $\alpha$ (2000) = 17h11m23.1058s, $\delta$ (2000) = -27${}^{\circ}$24${}^{\prime}$32.828${}^{\prime\prime}$ for 11A and $\alpha$ (2000) = 17h11m23.1015s, $\delta$ (2000) = -27${}^{\circ}$24${}^{\prime}$32.987${}^{\prime\prime}$ for 11B. There was a small error in the coordinates, as quoted in the Introduction of Vastel et al. (2022).

In this section, we describe the line identification method for the exploration of the molecular complexity in [BHB2007] 11. The line identification was performed using the CASSIS⁵⁵5http://cassis.irap.omp.eu software, which connects to the JPL (Pickett et al., 1998) and CDMS (Müller et al., 2005) databases with direct access to VAMDC⁶⁶6http://www.vamdc.org/ as well to verify all the available spectroscopic databases.

The maps presented in the following were obtained using the IRAM/GILDAS⁷⁷7http://www.iram.fr/IRAMFR/GILDAS package (version dec22a). Figure 2 shows the moment 0 maps for the two strong CH${}_{3}$OCH${}_{3}$ transitions. The spectra from both transitions were extracted from a $\sim$ 0.4${}^{\prime\prime}$ beam size, corresponding to the beam for both transitions (see Table 4), which includes the 11A and 11B cores. This corresponds to $\sim$ 30$\sigma$ for the 258.5 GHz transition (upper energy level of 93 K) and also corresponds to $\sim$ 4$\sigma$ for the 257.9 GHz emission (upper energy level at 191 K). The line profiles are presented in Fig. 3. The 257.911 GHz line profile is more widespread in frequency and does not present a simple three-component feature because of some blending with CH${}_{3}$OCHO transitions at 257.9 GHz (see Table 4, $\rm E_{up}$ = 306.01 K). The 258.549 GHz transition of dimethyl ether is the strongest line of the detected iCOMs presented in this study that is not affected by blending from other transitions. The profile is a simple three-fit Gaussian, and we used it to try to determine the exact location, assuming that the emission is compact. In Fig. 2 (bottom), the image component size deconvolved from the beam for the central component (2.8 km s${}^{-1}$) is more extended than the other, and we cannot yet deduce a clear location (see more details in Sec. 6). Furthermore, we used the unblended 258.549 GHz transition of dimethyl ether to compute the velocity field (moment 1) around [BHB2007] 11. The left side of Fig. 4 shows a clear northwest-southeast gradient centred in between 11A and 11B at the systemic velocity of the source. The same velocity field has been obtained for the methanol transition (Vastel et al., 2022).

In order to better understand the left side of Fig. 4, we created a position-velocity (PV) diagram for the same transition. The PV plot shown on the right side of the figure was produced from a cut oriented parallel to the velocity gradient observed in the moment 1 map of the dimethyl ether emission at 258.549 GHz, whose PA is 140${}^{\circ}$ east of north. The cut width is 0.5${}^{\prime\prime}$, comparable to the angular resolution of the FAUST data, and it encompasses cores 11A and 11B, which are separated by 0.2${}^{\prime\prime}$. This figure clearly shows the three velocity structures that were identified with methanol (Vastel et al., 2022) and does not show a smooth Keplerian profile over the velocity range from the circumbinary envelope, as was shown using formaldehyde by Alves et al. (2019).

Considering the angular resolution of the present observations, we cannot determine the exact location of the iCOMs emission within the binary system. However, we applied the same method as in Vastel et al. (2022) and present on the left side of Fig. 5 the red- and blueshifted emission for CH${}_{3}$OCH${}_{3}$ in the high-velocity channels where the S/N is higher than five (the channel at -3.69 km s${}^{-1}$ in the blue and the channel at 11.18 km s${}^{-1}$ in the red). Both components tend to be aligned with source 11B, like the high-velocity components of the CO emission (Alves et al., 2019) and the methanol emission (Vastel et al., 2022). The emission in the 2.8 km s${}^{-1}$ channel (yellow) peaks near 11A in the same way as CO and CH${}_{3}$OH. The three CH${}_{3}$OCH${}_{3}$ components appear in a compact configuration that is compatible with the compact emission seen in the moment 0 map (see Fig. 2).

Since the iCOMs emission in each velocity channel is compact and centrally peaked, we performed 2D Gaussian fits on the CH${}_{3}$OCH${}_{3}$ transition at 258.549 GHz (the strongest unblended transition with a high S/N) in each channel to obtain the position of their intensity peak in the three components [-5, 1] km s${}^{-1}$, [1, 6] km s${}^{-1}$, and [6, 13] km s${}^{-1}$. The right side of Fig. 5 shows this 2D analysis. The accuracy on the relative position is given by $\theta/(2\times S/N)$, where $\theta$ is the synthesised beam size, and S/N is the signal-to-noise ratio of the peak with respect to the map noise (Condon, 1997; Condon et al., 1998). Consequently, the uncertainty on the position will be lower for higher S/N line intensities. In other words, a centroid fitting to the compact molecular emission yields astrometry with a precision better than that of the spatial resolution in the case of a high S/N. For this calculation, we only considered intensity peaks brighter than five times its Gaussian fit error and a S/N higher than ten, corresponding to a position accuracy better than 0.02${}^{\prime\prime}$. This 2D analysis is similar to our analysis of methanol, except for the middle component at $\sim$ 3 km s${}^{-1}$, which is more localised on the B source instead of in between the two cores. We note, however, the limited number of channels satisfying the requirements (S/N $>$ 10, intensity/error $>$ 5) compared to methanol. Fig. 5 (right) shows no clear acceleration of the iCOMs material towards [BHB2007] 11 B (circles with red edges for velocities greater than 6 km s${}^{-1}$ and circles with blue edges for velocities lower than 1 km s${}^{-1}$), in contrast to what is seen in the high-velocity component of CO (circles with green edges) with a darker colour inside the circles nearby the source.

Many transitions of methyl formate (CH${}_{3}$OCHO) were detected in the continuum band (Fig. 1), and Fig. 6 shows the moment 0 maps for three isolated transitions covering a wide upper energy range between 115 and 333 K. The smaller beam of the transition at 233 GHz shows that the emission peaks south of [BHB2007] 11 B appear to clearly show a lack of emission towards [BHB2007] 11 A. There is a good correlation between this higher intensity in the southern region with the kinematic structure in the form of a clear velocity gradient around source [BHB2007] 11 B (see Fig. 5, right). Figure 7 shows the spectra of some isolated transitions of CH${}_{3}$OCHO in the continuum spw.

The left side of Fig. 8 shows the moment 0 map for the only transition of ketene, again, with an enhancement in intensity in the southern region of 11B, following the filament traced with the dust continuum emission presented by Alves et al. (2019), as shown in Fig. 5 (right).

We first performed a local thermodynamic equilibrium (LTE) analysis to estimate the CH${}_{3}$OCH${}_{3}$/CH${}_{3}$OH and CH${}_{3}$OCHO/CH${}_{3}$OH abundance ratios. The methanol column density was computed from non-LTE modelling using its ${}^{13}$C isotopologue (Vastel et al., 2022). For methyl formate and dimethyl ether, we extracted the spectra from the same region as for the CH${}_{3}$OCH${}_{3}$  transition at 257.911 GHz, corresponding to $\sim$ 0.4${}^{\prime\prime}$. We then used the parameters (column density, source size, and kinetic temperatures) found from the non-LTE modelling of the methanol transitions (see Vastel et al., 2022) to compute the column densities of CH${}_{3}$OCH${}_{3}$ and CH${}_{3}$OCHO using the CASSIS software. We fixed the source sizes and let the best-fit model converge around the $V_{LSR}$ line width for the Gaussian fits and the kinetic temperatures of methanol for the excitation temperatures of the three components. We then let the source size vary. Table 2 lists the results from the best fit obtained using all transitions in both continuum windows and in setup 2, spw 47 and 43, taking into account the blended lines of CH${}_{3}$OCH${}_{3}$  and CH${}_{3}$OCHO. Figures 3 and 7 show the best-fit LTE model (in red) obtained by simultaneously fitting all lines of CH${}_{3}$OCH${}_{3}$  and CH${}_{3}$OCHO. The corresponding $V_{LSR}$ for the three components are -1.9, 3.0, and 9.2 km s${}^{-1}$, and the FWHM values are 2, 5, and 4 km s${}^{-1}$ respectively.
With methanol as a reference, the CH${}_{3}$OCH${}_{3}$/CH${}_{3}$OH ratio is [5, 8] $\times$ 10${}^{-3}$ and the CH${}_{3}$OCHO/CH${}_{3}$OH ratio is [2.5, 4.5] $\times$ 10${}^{-2}$ for the three components.

We detected a line at 246.79 GHz (see Fig. 11) that may correspond to the blending of four 13C${}_{2}$ methyl formate (CH${}_{3}$O${}^{13}$CHO) transitions. The spectroscopic parameters are listed in Table 6. Many lines have previously been detected in the spectrum of the Orion star-forming region (Carvajal et al., 2009) at the IRc2 position. Combining our CH${}_{3}$OCHO column density with its isotopologue, we used a ${}^{12}$C/${}^{13}$C ratio of 68 (Milam et al., 2005) to compute a CH${}_{3}$O${}^{13}$CHO ratio of [1.5, 5.0] $\times$ 10${}^{14}$, [3.2, 15] $\times$ 10${}^{14}$ and [4.4, 19] $\times$ 10${}^{14}$ cm${}^{-2}$. With the same source size and excitation temperatures as for CH${}_{3}$OCHO and CH${}_{3}$OCH${}_{3}$, the modelled transitions using an LTE radiative transfer analysis in the two continuum bands are compatible with the observations within the noise. However, the LTE spectrum (Fig. 11, at the frequency of the middle group of 246788.406 MHz) presents a shift in velocity for the three components that might be explained by an unidentified transition at this position that is not reported in the CDMS and JPL databases, or by an incorrect frequency in the CH${}_{3}$O${}^{13}$CHO database. The quoted uncertainty of the database frequency is only 2 kHz, which cannot explain this shift. The 246.788 GHz line also appears in the Orion spectra (Fig. 2 of Carvajal et al., 2009), and it is well underestimated by the model with a clear shift as well. However, the authors did not explain this particular shift. Considering the noise and the spectral resolution of our observations for this transition, we cannot conclude on a firm detection of the methyl formate isotopologue.

We used the CH${}_{3}$OCHO and CH${}_{3}$OCH${}_{3}$ LTE modelling to obtain a good estimate of the formic acid column density. The 246.106 and 247.514 GHz lines are blended with methyl formate and are underestimated by the LTE model when formic acid is not taken into account. Taking the three components with the same beam size, excitation temperature, and $V_{LSR}$ as derived from the methyl formate and dimethyl ether into account, we obtained the column densities quoted in Table 2. The 246.106 and 247.514 GHz lines are now better reproduced by the LTE model. When methanol is again used as a reference, we obtain a t-HCOOH/CH${}_{3}$OH  ratio of [0.01, 0.5] $\%$.

We present on the right side of Fig. 8 the only transition for ketene, as discussed in Sect. 3. The column densities can be constrained from the LTE modelling, again using the excitation temperature, filling factor, and FWHM as for CH${}_{3}$OCH${}_{3}$ and CH${}_{3}$OCHO. The results are presented in Table 2 with a H${}_{2}$CCO/CH${}_{3}$OH ratio of [0.01, 0.1] $\%$, which is very similar to the t-HCOOH/CH${}_{3}$OH ratio.

An upper limit for acetaldehyde column density can be estimated at about 1/70th compared to the methanol column density. In Fig. 7, the transition at 234.486 GHz (see Table 4) is not well reproduced. However, it might be contaminated by an acetaldehyde transition at 234.486 GHz (see Table 6), which is within the upper limit value calculated above (0.15 K for this transition for a 0.115 K $rms$). We also evaluated the upper limit of the C${}_{2}$H${}_{5}$CN and C${}_{2}$H${}_{3}$CN column densities for each velocity component with the LTE analysis. The derivation of the upper limits is described in the appendix. We obtained C${}_{2}$H${}_{5}$CN/CH${}_{3}$OH and C${}_{2}$H${}_{3}$CN/CH${}_{3}$OH to be $<$0.066 % and $<$0.018 %, respectively.

We presented the ALMA observations towards the protobinary system [BHB2007] 11 within the framework of the FAUST Large Program. We reported the detection of iCOMs and iCOMs precursors, such as CH${}_{3}$OCH${}_{3}$, CH${}_{3}$OCHO, t-HCOOH, and H${}_{2}$CCO. The spatial distribution showed that the emission is confined within the ALMA beam, and the velocity distribution shows three different spatial regions. They are comparable with the methanol emission from previous observations of Vastel et al. (2022). The detected iCOMs are associated with three hot spots, two of which lie close to 11B, and the third hot spot lies near 11A, with a source size of $\sim$ 0.15${}^{\prime\prime}$, which is well below the $\sim$ 0.4${}^{\prime\prime}$ beam size. These molecular species are likely the result of the shocked material from the incoming filaments that streams towards [BHB2007] 11 A and B, shown in our integrated intensity maps compared to the dust emission revealing the filamentary structures. Radiative transfer modelling was used to reproduce the many detected transitions and compute the abundance ratios with respect to methanol for methyl formate, dimethyl ether, ketene, formic acid, and an upper limit for acetaldehyde. We then compared the ratios to Class 0 and Class I protostars as well as with the L1157-B1 shocked region. The CH${}_{3}$OCHO/CH${}_{3}$OH as well as CH${}_{3}$CHO/CH${}_{3}$OH ratios are similar in all sources, but H${}_{2}$CCO/CH${}_{3}$OH and t-HCOOH/CH${}_{3}$OH are much lower. Our CH${}_{3}$OCH${}_{3}$/CH${}_{3}$OCHO ratio of [0.1, 0.3] is slightly lower (factor of 3 to 10) than the correlation found in low- and high-mass protostars. Our observations suggest that the shocks caused by the streaming gas crashing onto quiescent material surrounding [BHB2007] 11 can enrich the gas in iCOMs.