Open Access
Issue
A&A
Volume 667, November 2022
Article Number A20
Number of page(s) 15
Section Interstellar and circumstellar matter
DOI https://doi.org/10.1051/0004-6361/202244312
Published online 01 November 2022

© E. Artur de la Villarmois et al. 2022

Licence Creative CommonsOpen Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (https://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

This article is published in open access under the Subscribe-to-Open model. Subscribe to A&A to support open access publication.

1 Introduction

The formation and evolution of protoplanetary disks are fundamental in the process of low-mass star formation, such as the formation of our own Solar System. A typical low-mass star forms when a molecular cloud with angular momentum collapses, and a protostar is formed at the central part with an infalling-rotating envelope whose inner part evolves to a circumstellar disk (Terebey et al. 1984; Shu et al. 1993; Hartmann 1998). Eventually, the star reaches its final mass, the envelope dissipates, and planets form in the disk. As a consequence, the final composition of planets is strongly dependent on the physical and chemical processes within the circumstellar disk. However, as disks first arise in the early stages of young stars (Jørgensen et al. 2009; Harsono et al. 2014; Yen et al. 2015) and the first steps of planet formation may occur when they are still deeply embedded (e.g., Harsono et al. 2018; Tychoniec et al. 2020), the chemical evolution of the material as it is accreted from the infalling envelope may play a key role.

The process of low-mass star formation comprises different stages (Robitaille et al. 2006), and Class I sources link the deeply embedded Class 0 sources (where the envelope is the dominant mass component) with the emergence of Class II disks (Keplerian disks with a negligible envelope). Class I sources are therefore the perfect candidates to study the connection between the envelope and the disk and, additionally, to investigate the dynamics and chemical composition of the young disk.

Theoretical models predict that the material from the envelope falls on the circumstellar disk and produces accretion shocks at the envelope-disk interface (Stahler et al. 1994; Yorke & Bodenheimer 1999; Krasnopolsky & Königl 2002). These accretion shocks have been invoked to explain the observed jump in density and drastic enhancement of SO toward the Class 0 and I sources L1527 and TMC-1A (Sakai et al. 2014, 2016), the asymmetric accretion found toward TMC-1A (Hanawa et al. 2022), and the emission of SO and SO2 at the edge of the disks from two Class I/II sources, DG Tau and HL Tau (Garufi et al. 2022). Accretion shocks in dense (≥108 cm−3) gas induce an increase in the dust temperature, and species that are locked in grain mantles are subsequently released into the gas phase, which affects the chemical content of the early disk (van Gelder et al. 2021). Although the presence of accretion shocks explains the jump in abundances observed for shock-related species and is the most plausible mechanism deduced from numerical simulations (Miura et al. 2017), only a few low-mass protostars show evidence of accretion shocks to date (e.g., Lee et al. 2014; Sakai et al. 2014; Garufi et al. 2022), and their physical parameters are not well constrained observationally. Apart from accretion shocks, contributions from disk winds or outflows would also be important (e.g., Bjerkeli et al. 2016; Alves et al. 2017; Tabone et al. 2017; Harsono et al. 2021). Therefore, observations at disk scales (~ 100 au) need to be performed to confirm the existence of accretion shocks, understand the origin of the observed abundances, and assess the physical parameters associated with this mechanism.

A suitable source for proving the nature of accretion shocks is Oph-IRS 44, a Class I source located in the Ophiuchus molecular cloud at a distance of 139 pc (average value for the L1688 cloud; Cánovas et al. 2019). Artur de la Villarmois et al. (2019) detected strong SO2 emission toward a compact region (≤60 au) in IRS 44, with an angular resolution of 0/4 (~60 au). This particular SO2 transition (184,14–183,15) is associated with an upper-level energy (Eup) of ~200 K and its line profile shows a velocity range of ~20 km s−1. The angular resolution of 0/4 of the data was not high enough to resolve the SO2 emission and provide strong conclusions for the possible origin scenarios: accretions shocks, disk winds, or outflows.

IRS 44 was first identified as YLW 16A by Young et al. (1986) through IRAS observations, and other common names are Oph-emb 13, ISO-Oph 143, LFAM 35, and [GY92] 269, among others. It is associated with a bolometric temperature (7W) of 280 K, a bolometric luminosity (Lbol) of 7.1 L (Evans et al. 2009), and an envelope mass (Menv) of 0.051 M (for a distance of 139 pc; Jørgensen et al. 2009). IRS 44 has been proposed to be a protobinary system with a separation of ~0″.3, based on observations with the Hubble Space Telescope (HST; Allen et al. 2002), the Very Large Telescope (VLT; Duchêne et al. 2007), and the Spitzer Space Telescope (McClure et al. 2010). Nevertheless, there is no evidence of a binary component in the submillimeter regime, through ALMA band 6 and band 7 observations (Sadavoy et al. 2019; Artur de la Villarmois et al. 2019).

In this paper we present high angular resolution 0″.1 (14 au) ALMA observations of multiple SO2 molecular lines toward IRS 44. We discuss their potential to trace accretions shocks, and provide values of the physical parameters for the emitting gas. Section 2 describes the observational procedure, calibration, and the parameters of the observed molecular transitions. The observational results are presented in Sect. 3, while Sect. 4 is dedicated to the analysis of the data, with position-velocity diagrams, radiative-transfer models, estimations of rotational and excitation temperatures, and calculations of molecular column densities. We discuss the structure and kinematics of IRS 44 in Sect. 5, and end with a summary in Sect. 6.

2 Observations

IRS 44 was observed with ALMA during 2021 May 17 and 18 as part of the program 2019.1.00362.S (PI: Elizabeth Artur de la Villarmois). At the time of the observations, 47 and 45 antennas were available, respectively, in the array providing baselines between 15 and 2517 m. The observations targeted nine different spectral windows to observe multiple SO2 lines, the less abundant 34SO2 isotopolog, and SO. The observed molecular transitions and their spectroscopic data are summarized in Table 1.

The calibration and imaging were done in CASA1 version 6.1.1 (McMullin et al. 2007). Gain and bandpass calibrations were performed through the observation of the quasars J1517-2422 and J1700-2610. Imaging was performed using the tclean task in CASA, where the Briggs weighting with a robust parameter of 0.5 was employed. The automasking option was chosen and the channel resolution is 0.21 km s−1. The resulting dataset has a beam size of 0″.13 × 0″.09 (18 × 13 au) with a position angle (PA; measured from north to east) of −81° and a largest angular scale (LAS) of 2″.3. The continuum rms level is 0.08 mJy beam−1 and the rms level of each spectral window is listed in Table 1.

3 Results

3.1 Continuum emission

The continuum emission is shown in Fig. 1, where the horizontal component is slightly more extended than the vertical component, and the emission above 5σ is contained within a radius of 0″.2 (~30 au). Two-dimensional (2D) Gaussians are used to fit emission in the image plane, obtaining an integrated flux of 22.9 ± 1.0 mJy, a peak flux of 16.91 ± 0.48 mJy beam−1, and a deconvolved size of (0″.07 ± 0″.01) × (0″.06 ± 0″.01) with a PA of 119 ± 74° (see the magenta ellipse in Fig. 1). The continuum peak position corresponds to α = 16h27m27s.9858 ± 0s.0002 and δ = −24°39′34″.063 ± 0″.001.

The disk mass at 0.87 mm was calculated from the continuum flux (22.9 ± 1.0 mJy) and using Eq. (2) from Artur de la Villarmois et al. (2018), which assumed optically thin emission, an opacity of 0.0175 cm−2 per gram of gas at 0.87 mm, and a dust temperature of 30 K. A total mass Mgas+ dust of (4.0 ± 0.2) × 10−3 M was obtained, adopting a dust temperature (Tdust) of 15 K, the value proposed by Dunham et al. (2014) for Class I sources. If Tdust = 30 K is assumed, the total mass decreases by a factor of ~3. Given that the dust emission at 0.87 mm could be optically thick toward a Class I source, the calculated Mgas+dust represents a lower limit for the total mass.

Allen et al. (2002) and Duchêne et al. (2007) suggested that IRS 44 is a protobinary system with a separation of ~0/3. However, the continuum emission at 0.87 mm shows no binary detection in our ALMA data and we can only set an upper limit of 7 × 10−5 M for the total mass of a possible binary component (for a value of 5σ and adopting the same parameters as in the previous paragraph).

Table 1

Spectral setup and parameters of the observed molecular transitions.

thumbnail Fig. 1

Continuum emission (0.87 mm) toward IRS 44 above 5σ (σ = 0.08 mJy beam−1). The white contours represent the weakest emission of [5, 10, 15, 20, and 25σ] for clarity. The black star shows the position of the continuum peak and the synthesized beam is indicated by the white filled ellipse. The magenta ellipse represents the deconvolved size from the 2D Gaussian fit.

3.2 Molecular transitions

All the targeted molecular lines listed in Table 1 were detected toward IRS 44, with the exception of the 34SO2 96,4–105,5 line. This nondetection is consistent with the low Einstein A coefficient (Aij = 4 × 10−5 s−1) of the transition and it being a less abundant isotopolog. The six detected SO2 lines have different upper level energies Eup, covering a broad range from 36 to 293 K. The brightest emission toward the continuum peak is from the SO2 184,14–183,15 line with Eup = 197 K, while there is an offset region located at a distance of ~3″.0 (~400 au) from the protostar that shows bright emission of the SO2 line related with the lowest energy: SO2 53,3–42,2 with Eup = 36 K.

Figure 2 presents the spectrum, and moment 0 and 1 maps of the SO2 184,14–183,15 line toward the central region. The spectrum was taken over a circular region with r = 0″.2 and shows a broad-line profile, from −20 to 20 km s−1, and a decrease in the emission around the systemic velocity (Vsys) of 3.7 km s−1, estimated from previous APEX observations (Lindberg et al. 2017). The moment 0 map reveals that the emission is concentrated around the protostar; however, the emission peak is slightly offset from the continuum peak, ~0″.1 (~14 au), and corresponds to the redshifted component. The moment 1 map shows a clear rotational signature from northwest to southeast, with a PA of 157 ± 3°. The PA for the SO2 184,14–183,15 line emission was obtained from a 2D Gaussian fit of the moment 0 map. We note that this PA value is not perpendicular to the outflow direction (PA = 20°), which was estimated by van der Marel et al. (2013) using single-dish observations of CO 3–2. For the other detected lines (five SO2, two 34SO2, and one SO line), the spectra and moment 0 and 1 maps are presented in the appendix, in Figs. A.1, A.2, and A.3, showing that SO2, 34SO2, and SO exhibit a similar nature: broad spectra, emission concentrated around the protostar, and a clear rotational signature. In addition, all the detected transitions show that the peak of emission is offset south from the continuum peak position, at a distance of ~0″.1 (~14 au). On average, the six SO2 transitions show a full width at half maximum (FWHM) value of 12 km s−1 for the blueshifted emission and 14 km s−1 for the redshifted emission. Integrated fluxes of the observed transitions are presented in Table A.1 in the appendix.

Figure 3 shows contour maps of the SO2 184,14–183,15 line for different velocity ranges. Low-velocity contours (between −2 and 2 km s−1) are concentrated around the protostar, but the weakest contours also present emission toward the west. Intermediate velocities, between ±6 and ±10 km s−1, show that the red-shifted emission is more extended than the blueshifted emission, possibly related with a protrusion from a localized streamer. Finally, a clear and symmetric rotating signature around the protostar is seen for high velocities (≥10 km s−1), which will be referred to as a disk-envelope structure.

At larger angular scales, the SO2 53,3–42,2 line (Eup = 36 K) shows bright emission toward an offset region, located at a distance of 2″.8 (~400 au) from the protostar, and its spectra and moment 0 map are presented in Fig. 4. The spectra were taken over the central region (gray) and the offset region (red), revealing that the offset region is associated with low velocities (≤2 km s−1), in contrast with the broad-line profile observed toward the central region, suggesting a different and more quiescent origin. Two other SO2 lines (with Eup value of 197 and 199 K) show weaker emission toward the offset region and their moment 0 maps are presented in Fig. A.4. Given that the SO2 line with the lowest Eup value (36 K) shows the brightest emission, the offset region is associated with colder gas, more consistent with a cloudlet or an SO2 knot with a PA of 125 ± 7°. In this case the PA value was calculated by projecting a line that connects the continuum peak with the brightest pixel of the offset region, and it is consistent with the direction of the weakest contours seen in the first two panels of Fig. 3. The offset region is henceforth referred to as an SO2 knot.

thumbnail Fig. 2

Emission of the SO2 184,14–183,15 line. Left: spectrum rebinned by a factor of 4, integrated over a circular region with r = 0″.2, and centered on the continuum peak position. Center, moment 0 map above 3σ (color scale), integrated over 60 km s−1, and continuum emission (black contours), starting at 20σ and following steps of 40σ. Right: moment 1 map (color scale) and moment 0 map (white contours) above 3σ. The blue and red arrows show the direction of the outflow from van der Marel et al. (2013), the yellow star indicates the continuum peak position, and the synthesized beam is represented by the black filled ellipse in the bottom right corner. The adopted systemic velocity is 3.7 km s−1.

thumbnail Fig. 3

Contour maps of SO2 184,14–183,15 for different velocity ranges, shifted to 0 velocity. The contours start at 5σ and follow steps of 5σ. The panel that includes the systemic velocity (3.7 km s−1) is shown by the green contours in the first panel, while blue- and redshifted emission is represented by the blue and red contours, respectively. The yellow star shows the position of the source and the synthesized beam is indicated by the black filled ellipse in the bottom left corner of the first panel.

thumbnail Fig. 4

Emission of the SO2 53,3–42,2 line. Left: spectra integrated over a circular region with r = 0″.2, centered on the continuum peak position (gray) and centered on the offset region (red). The spectrum taken at the continuum peak position has been rebinned by a factor of 4. Right: moment 0 map above 3σ integrated over 60 km s−1. The gray and red circles represent the regions from which the spectra in the left panel were taken. The blue and red arrows show the direction of the outflow, the yellow star indicates the continuum peak position, and the synthesized beam is represented by the black filled ellipse in the bottom right corner. The adopted systemic velocity is 3.7 km s−1.

thumbnail Fig. 5

Position-velocity diagram for SO2 184,14–183,15. Left: emission above 3σ% employing a PA of 157°. The blue and red dots represent blue- and redshifted emission peaks above ±3 km s−1, respectively. The black line represents an infalling-rotating profile with M* = 1.5 M, an inclination of 70°, and rCB = 0708 (magenta dashed lines). The rCR is shown by the green dashed lies. Center: zoomed-in version for the redshifted emission. Right: zoomed-in version for the blueshifted emission. The systemic velocity is 3.7 km s−1.

4 Analysis

4.1 Position-velocity diagrams

Figure 5 shows a position-velocity (PV) diagram for the SO2 184,44–183,15 line, employing a PA of 157°, with the peak emission of each channel superimposed. The peak emission was obtained through the CASA task imfit and the offset position was calculated by projecting the peak emission onto the disk position angle. The redshifted emission is more extended than the blueshifted emission: the former shows emission up to 0″.35 (~50 au) and the latter up to 0″.21 (~30 au). The central and right panels of Fig. 5 are zoomed-in versions of the red- and blueshifted emission. The high-velocity points are best fitted with an infalling-rotating profile (Vrot+inf), employing the equation Vrot+inf=2GM*sin(i)rCBr,${V_{{\rm{rot + inf}}}} = {{\sqrt {2G{M_*}\sin \left( i \right){r_{{\rm{CB}}}}} } \over r},$(1)

where G is the gravitational constant, M* the protostellar mass, rCB the radius of the centrifugal barrier, i the inclination of the disk, and r the distance from the protostar. Equation (1) is from Oya et al. (2014), and the inclination term has been added explicitly. The rCB is given by the maximum radial velocity (Sakai et al. 2014) and can be estimated from the PV diagram; the maximum radial velocity of 17.5 km s−1 corresponds to rCB = 0″.08 (~11 au). The maximum radial velocity changes depending on the assumption of Vsys; therefore, if Vsys changes by 0.5 km s−1, rCB will change by 0701 (~1.4 au). This leaves us with a degeneracy in the protostellar mass and the inclination, given by the term M*sin(i). A protostellar mass of 1.5 M is obtained if an inclination value of 70° is assumed, following the interpretation of Terebey et al. (1992) that the outflow axis of IRS 44 lies close to the plane of the sky (from VLA observation of water masers). If we use other inclination values, such as 50° and 90°, the points are well fitted with an infalling-rotating profile with a M* of 1.8 M and 1.4 M, respectively. Seifried et al. (2016) proposed that the protostellar mass can be estimated by fitting the maximum velocity offset in the PV diagram (i.e., the borders above 3σ, which correspond to the outer envelope), instead of fitting the peak emission of each channel. Following this procedure, a protostellar mass of 4 M is obtained.

The centrifugal barrier is the radius at which most of the gas kinetic energy contained in infalling motion is converted to rotational motion. The gas motion of the disk-envelope system outside rCB can be regarded as infalling-rotating motion, while that inside can be regarded as Keplerian motion (Sakai et al. 2014; Oya et al. 2018). The rCB is half of the centrifugal radius (rCR) beyond which the gas is falling (Oya et al. 2018). From the SO2 184,14–183,15 line, rCB = 0″.08 (~11 au) and rCR = 0″.16 (~22 au). Beyond rCR the more extended redshifted emission is seen, while no blueshifted counterpart is observed. This is consistent with the redshifted protrusion seen in the contour maps of Fig. 3 at velocities between 6 and 10 km s−1, suggesting that a localized streamer might be infalling toward the system and, when entering the centrifugal radius at 0″.16, an infalling-rotating profile dominates the dynamics. A Keplerian disk is expected inside the centrifugal barrier of 0″.08; however, this is close to the resolution of our data and the presence of a Keplerian disk is not conclusive with the current data. If a Keplerian disk exists toward IRS 44, its radius will be ≤0″.08 (~11 au). Given that no Keplerian motions are observed in our data, the rotational signature seen in the moment 1 map of SO2 (Fig. 2) suggests the presence of a disk-envelope structure and not a rotationally supported disk.

4.2 Column densities, kinetic temperatures, and optical depth

In this section we estimate kinetic temperatures (Tkin), SO2 and SO molecular column densities NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ and NSO), and the optical depth of the lines by employing the non-LTE radiative transfer code RADEX (van der Tak et al. 2007). Later on, rotational temperatures (Trot) and NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ of optically thin lines are estimated from the rotational diagram method, and excitation temperatures (Tex) are assessed from optically thick lines.

4.2.1 Radiative transfer

The six different SO2 transitions were employed to derive the gas density and temperature by comparing the observed relative intensities with those predicted by RADEX. The observed relative intensities are the quotient between the moment 0 maps, which present emission up to a radius of ~072. RADEX was run for a set of kinetic temperatures from 30 to 300 K, SO2 column densities from 1012 to 1018 cm−2, and H2 number density nH between 103 and 109 cm−3. Collisional rates for SO2 were taken from the Leiden atomic and molecular database (LAMDA; Balança et al. 2016). A value of 5 km s−1 was used for the broadening parameter (b), which corresponds to the line width observed in pixels far from the SO2 peak. The brightest SO2 line, which is associated with an Eup of 197 K, is used as a reference line.

RADEX models with nH between 103 and 107 cm−3 were unable to explain the observed line ratios. The observed relative intensities are shown in Fig. B.1, and they are compared with RADEX results for a H2 number density of 109 and 108 cm−3. The observed values provide a range of possibilities for Tkin and 8×1016NSO28×1017${N_{{\rm{S}}{{\rm{O}}_2}}}$, given nH∙ For nH = 108 cm−3, there are no possible values that satisfy all the observed ranges, implying that the SO2 emitting region is associated with nH > 108 cm−3. On the other hand, for nH = 109 cm−3, Tkin should be higher than 90 K and NSO2$8 \times {10^{16}} \le {N_{{\rm{S}}{{\rm{O}}_2}}} \le 8 \times {10^{17}}$ cm−2. This possible values are shown in Fig. B.2.

For nH = 109 cm−3, the optical depth of the six SO2 lines is analyzed, taking into account the possible values of Tkin and NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$. Figure B.3 shows that, from the six SO2 lines, two are optically thick (SO2 184,14–183,15 and SO2 201,19–192,18), two are optically thin (SO2 167,9–176,12 and SO2 106,4–115,7), and nothing conclusive can be said about the remaining two (SO2 53,3–42,2 and SO2 242,22–233,21).

thumbnail Fig. 6

Temperature structure of IRS 44. Left: rotational diagram at the source position, where only the optically thin SO2 transitions (blue dots) and the 34SO2 lines (red dots) are used for the fit. The abundance ratio 32S/34S = 22 is from Wilson (1999). Optically thick SO2 transitions (green dots) and those lines where the optical depth is not conclusive (open dots) show a significant offset with respect to the other lines, and they were not included in the calculation of the rotational temperature and SO2 column densities. Center, rotational temperature map created from optically thin SO2 and 34SO2 transitions. Contours show specific values of 140, 160, 180, and 200 K. Right: moment 1 map (right panel of Fig. 2) with the same specific values of the rotational temperature as those of the central panel, showing that the highest temperatures coincide with the redshifted protrusion. The synthesized beam is shown by the black ellipse in the upper left corner and the yellow star indicates the position of the source.

4.2.2 Rotational diagram

For optically thin SO2 lines and the less abundant isotopolog 34SO2, the beam-averaged column densities and rotational temperatures can be assessed by the rotational diagram analysis, summarized by Goldsmith & Langer (1999). The gas is assumed to be under local thermodynamic equilibrium (LTE); therefore, all the molecular transitions can be characterized by a single excitation temperature, also called rotational temperature (Trot). In this regime the following equation is valid: lnNugu=1TrotEuk+lnNQ(Trot).$\ln {{{N_{\rm{u}}}} \over {{g_{\rm{u}}}}} = {{ - 1} \over {{T_{{\rm{rot}}}}}}{{{E_{\rm{u}}}} \over k} + \ln {N \over {Q\left( {{T_{{\rm{rot}}}}} \right)}}.$(2)

Here Nu is the column density of the upper level, gu the level degeneracy, Eu/k the energy of the upper level in K, k the Boltzmann constant, N the total column density of the molecule, and Q(Trot) the partition function that depends on the rotational temperature.

Under the optically thin condition, Nu is obtained from Nu=8πkv2Whc2Aul,${N_{\rm{u}}} = {{8\pi k{v^2}W} \over {h{c^2}{A_{{\rm{ul}}}}}},$(3)

where v is the line frequency, W the integrated line intensity, c the speed of light, and Aul the Einstein coefficient for spontaneous emission. Equation (3) can be rewritten as Nu=1943.59(v1 GHz)2(W1 K km s1)(1s1Aul),${N_{\rm{u}}} = 1943.59{\left( {{v \over {1\,{\rm{GHz}}}}} \right)^2}\left( {{W \over {1\,{\rm{K}}\,{\rm{km}}\,{{\rm{s}}^{ - 1}}}}} \right)\left( {{{1\,{{\rm{s}}^{ - 1}}} \over {{A_{{\rm{ul}}}}}}} \right),$(4)

where Nu is obtained in units of cm−2.

Equations (2) and (4) were used to calculate Trot and create the map shown in Fig. 6. For each pixel, only the optically thin SO2 lines and 34SO2 isotopologs were used to fit the rotational temperature. For 34SO2, an abundance ratio 32S/34S = 22 (Wilson 1999) was adopted. The left panel of Fig. 6 shows the example of the fit from the pixel that corresponds to the source position and a clear offset is seen between optically thin (blue and red dots) and optically thick lines (green dots). The detection of optically thin lines and the less abundant isotopolog, 34Sθ2, is crucial for an accurate estimate of the rotational temperature, and consequently for the SO2 column density as well. The Trot map (central and right panels of Fig. 6) shows high temperatures (≥120 K) in the region where the SO2 emission arises. In addition, the warmest region, southeast from the protostar, seems to correlate with the redshifted protrusion. When infalling material reaches the surface layers of the disk-envelope structure, it generates accretion shocks that are predicted to increase the temperature and the density by up to two orders of magnitude (~109 cm−3; van Gelder et al. 2021). If the dust temperature exceeds 60 K, SO2 molecules can efficiently desorb from dust grains.

Figure 7 shows the SO2 and SO column densities, and the ratio between them. The region where the six SO2 lines are detected shows NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ values between 1.0 and 1.8 × 1017 cm−2, while NSO presents lower values, between 0.6 and 1.3 × 1017 cm−2. Since there is only one observed (and detected) SO line, RADEX was employed using the same temperature and density parameters as SO2 (i.e., nH = 109 cm−3 and Tkin = 90 K), concluding that this SO line in particular is optically thin. The column density ratio between SO2 and SO is shown in the right panel of Fig. 7 and it is found to be higher than 1 toward the SO2 emitting region.

thumbnail Fig. 7

Column densities. Left: SO2 column density obtained from the rotational diagram method. The synthesized beam is shown by the black ellipse in the upper left corner. Center: SO column density, assuming optically thin emission and employing the rotational temperatures from Fig. 6. Right: column density ratio of SO2 to SO. The yellow star indicates the position of the source.

4.2.3 Optically thick lines

As seen in Sect. 4.2.1 two out of six SO2 lines are optically thick. For optically thick lines, the peak temperature (Tpeak) provides a good measure of (Tex) with Tex=Eupk[ log(EupkTpeak+1) ]1.${T_{{\rm{ex}}}} = {{{E_{{\rm{up}}}}} \over k}{\left[ {\log \left( {{{{E_{{\rm{up}}}}} \over {k{T_{{\rm{peak}}}}}} + 1} \right)} \right]^{ - 1}}.$(5)

Equation (5) is from Goicoechea et al. (2016) and, if nH is much higher than the critical density of the transition (ncrit), the line is close to thermalization and Tex approaches Tgas. From the two optically thick SO2 lines, the line with the lowest ncrit (~107 cm−3) corresponds to SO2 184,14–183,15. This transition line was used to create the temperature map shown in Fig. 8, where Tpeak was obtained from a moment 8 map (which provides the maximum value of the spectrum in each pixel). The southern region presents a more extended and elongated structure in the excitation temperature map, consistent with the redshifted protrusion (see also Figs. 2, 3, and 6), and Tex ≥ 70 K are found for the SO2 emitting region. Given that the τ value of this line lies between 1 and 7 (see first panel of Fig. B.3), it may not be fully thermalized, and therefore the temperature map in Fig. 8 represents a lower limit for Tex. These excitation temperatures are consistent with those found in Sect. 4.2.2 from the rotational diagram method using optically thin transitions.

5 Discussion

5.1 Accretion shocks, disk winds, or outflows?

The molecules SO and SO2 are known as shock tracers and there are three main physical origins for these shocks: outflows (e.g., Tafalla et al. 2010; Persson et al. 2012), disk winds (e.g., Tabone et al. 2017), and accretion shocks (e.g., Sakai et al. 2014; Garufi et al. 2022).

For IRS 44 the outflow scenario can be ruled out from the shape of the PV diagram shown in Fig. 5 and the high densities (≥108 cm−3) found for the SO2 emitting region. PV diagrams related with outflow emission show that the velocity linearly increases as a function of the distance to the protostar (e.g., Lee et al. 2000; Arce et al. 2013) and densities below 108 cm−3 have been found in the inner regions of the outflow cavity associated with young protostars (Kristensen et al. 2013). In addition, the broadness of the SO2 lines rules out the envelope origin, where typical line widths are below 2 km s−1 (e.g., Harsono et al. 2021).

As disk winds are related with gas that is ejected at small radial distances from the central source (e.g., Bjerkeli et al. 2016; Alves et al. 2017), some degree of symmetry is expected on the surface layers of the disk, such as a butterfly shape. Tabone et al. (2017) have proposed that the SO and SO2 emission detected toward the Class 0 source HH212 originates from a disk wind between ~50 and ~150 au. Nevertheless, Panoglou et al. (2012) have shown that species such as SO survive between 10 and 100 au in disk winds toward Class 0 sources, but they get destroyed by photodissociation beyond ~ 1 au in disk winds from more evolved Class I sources. The SO2 emission does not show the expected symmetry for a disk wind and the kinematic analysis indicates that the material follows an infalling-rotating profile without a Keplerian signature. If disk winds are present, we expect them to arise from the disk surface layers, likely inside 0″.08 (11 au).

The high temperatures estimated from optically thin (≥120 K) and optically thick (≥70 K) lines (Figs. 6 and 8), the moderate velocities (between 12 and 14 km s−1), and the high densities (≥108 cm−3) found for IRS 44 are in agreement with the accretion shock scenario. van Gelder et al. (2021) have shown that accretion shocks can efficiently desorb SO2 from dust grains when moderate velocities (≥10 km s−1) and high densities (≥108 cm−3) are present. For densities above 3 × 104 cm−3, the gas and the dust are efficiently coupled, Tdust = Tgas (Evans et al. 2001; Galli et al. 2002), and a dust temperature above 62 K is required in order to sublimate SO2 molecules from dust grains (Penteado et al. 2017; van Gelder et al. 2021). In interstellar ices, SO2 is tentatively detected (Boogert et al. 1997; Zasowski et al. 2009); however, chemical models predict that SO2 is the most abundant species in the gas in the warm-up phase, when the protostar is formed (Woods et al. 2015).

Accretion shocks would also desorb SO molecules form dust grains and the gas-phase abundance of SO2 could increase through the reaction of SO with OH (Charnley 1997; van Gelder et al. 2021). Nevertheless, Karska et al. (2018) did not detect OH toward IRS 44 from Herschel/PACS observations, suggesting that the gas-phase formation of SO2 by oxidation of SO could be ruled out. SO2/SO ≥ 1 also suggests that the radiation field from the protostar is not efficiently photo-dissociating SO2 into SO (e.g., Booth et al. 2021) and that the cosmic ray ionization rate is low (ζ = 1.3 × 10−17 s−1, Woods et al. 2015).

In this section we suggest that SO and SO2 molecules toward IRS 44 sublimate from heated dust grains by the accretion shocks with moderate velocity shocks (≥10 km s−1) and high densities (≥108 cm−3). If there is a chemical reaction that contributes to the SO2 abundance in the gas phase, it should be a different one from the reaction of SO with OH. Future observations of other molecular species, such as OCS, H2S, and H2CO, will confirm the formation path of SO2: direct desorption from dust grains, gas-phase formation, or a combination of both. H2S and H2CO are directly linked to the gas-phase formation of SO and SO2, while OCS presents a similar desorption temperature to SO, but it does not participate in the gas-phase chemistry (Charnley 1997).

thumbnail Fig. 8

Excitation temperature of the optically thick SO2 184,14–183,15 transition using Eq. (5). The black contour represents the specific value of 90 K. The yellow star indicates the position of the source and the synthesized beam is shown by the black ellipse in the upper left corner.

thumbnail Fig. 9

Schematic representation of IRS 44. The quiescent SO2 emission (knot) would be part of a streamer that allows material to fall into the disk-envelope through the redshifted protrusion and generate accretion shocks. These shocks sublimate SO2 molecules from dust grains and enhance the SO2 gaseous abundance, showing an emission peak ~0″.1 to the south of the protostar.

5.2 Morphology of IRS 44

Given that (i) quiescent and colder SO2 emission is present at ~2″.8, (ii) a redshifted protrusion is seen at velocities between 2 and 10 km s−1, (iii) the highest temperatures seem to correlate with the redshifted protrusion, (iv) blueshifted material beyond 0″. 2 (~30 au) is absent, and (v) the SO2 emission peak is observed at a distance of ~0″.1 from the continuum peak (at redshifted velocities), a localized streamer might be accreting material to the envelope-disk system and generating accretion shocks that release SO2 molecules from the dust to the gas phase. Figure 9 shows a schematic representation of IRS 44 and the localized streamer, which would be located between the observer and the disk-envelope component and would feed the system toward the redshifted protrusion. As IRS 44 is classified as a Class I source, meaning that the envelope is still present but largely dissipated (Menv = 0.051 M for a distance of 139 pc; Jørgensen et al. 2008, 2009), it is more likely that the infalling of material occurs through streamers and not in a spherically symmetric way. A similar behavior is seen toward the Class I source TMC-1A, as asymmetric CS and SO emission is explained by a cloudlet capture and subsequent formation of an infalling streamer (Hanawa et al. 2022), and the more evolved Class I/II sources DG Tau and HL Tau (Garufi et al. 2022), where accretion shocks traced by SO and SO2 are located along the late infalling streamers still feeding the system.

5.3 Outflow direction versus disk-envelope direction

As seen in Fig. 2, the disk-envelope direction (157°) is not perpendicular to that of the single-dish outflow (20°). The latter is seen at large scales (~30″) and the outflow direction may vary due to the surrounding gas. If the misalignment is real, this could be due to the presence of a binary component, a cloudlet capture, or physical processes during the formation history, such as misalignment between the cloud rotation axis and the initial B-field direction, formation from a turbulent core, or non-ideal magnetohydrodynamics (MHD) effects.

Terebey et al. (2001) proposed that IRS 44 is a protobinary system with a separation of 0727 and PA = 81°, based on HST observations. The primary component has been detected at 1.60, 1.87, and 2.05 µm, while the secondary component is only visible at the longest wavelength, at 2.05 µm. Nevertheless, there is no sign of binarity toward IRS 44 in the submillimeter regime, following this work with an angular resolution of ~0″.1 and the ALMA data presented in Artur de la Villarmois et al. (2019) and Sadavoy et al. (2019), which report an angular resolution of -074 and ~0″.25, respectively. The HST emission in 2.05 µm could therefore be associated with scattered light and not a binary component or the binary could be very faint at submillimeter wavelengths (below our sensitivity).

In the case of a cloudlet capture, each cloudlet should have a different angular momentum vector and the capture process can potentially change the rotation axis of the disk (e.g., Dullemond et al. 2019; Kuffmeier et al. 2020; Hanawa et al. 2022). The presence of a localized streamer toward IRS 44 could be affecting the rotation axis of the disk-envelope, and the result would depend on the mass and the angular momentum vector of the infalling material.

The misalignment between the cloud rotation axis and the initial B-field direction can create a warped disk structure during the protostellar core collapse (Hirano & Machida 2019), and B-fields in protostellar cores appear to be randomly aligned with their respective outflows (Hull et al. 2014; Lee et al. 2017). In the absence of a binary component, this initial misalignment could explain the change in direction observed toward IRS 44. A similar situation was proposed for the Class I disk L1489 (Sai et al. 2020), where the observation of a warped disk is explained by the initial misaligment between the initial B-field direction and the angular momentum vector.

Different velocity gradients between the direction of the rotationally supported disk and the direction of the envelope rotation were seen for a handful of Class 0/I sources (Brinch et al. 2007; Harsono et al. 2014). This misalignment might be due to formation from turbulent cores or non-ideal MHD effects, such as the Hall effect (e.g., Li et al. 2011; Braiding & Wardle 2012).

5.4 Nondetection of C17O and absence of warm CH3OH toward IRS 44

Previous observations of IRS 44 do not detect C17O (3–2) and warm CH3OH (Eup =65 K) at an angular resolution of 0″.4 (~60 au); Artur de la Villarmois et al. 2019). C17O is commonly associated with Keplerian disks in Class I sources, and its nondetection might be related with the absence of a Keplerian disk, at least outside 11 au. CH3OH, on the other hand, is hardly detected in Class I sources (Artur de la Villarmois et al. 2019); however, its gas-phase abundance is enhanced in shocked regions and a correlation between SO2 and CH3OH is expected. That SO2 shows strong emission and CH3OH is not detected toward IRS 44 might be related with one of the following possibilities: (i) CH3OH is being desorbed form dust grains, but later on it is destroyed by the moderate velocities of the shocks (≥ 10 km s−1; Suutarinen et al. 2014); (ii) the formation of CH3OH on the grain surfaces, from H2CO, is not efficient; (iii) the presence of a disk results in colder gas (Lindberg et al. 2014; van Gelder et al. 2022); or (iv) optically thick dust can hide the emission of COMs (Nazari et al. 2022). Future observations of H2CO could clarify the CH3OH nondetection, and clearly there is some uncertainty regarding the origin of the SO2 emission. The origin of this may be related to the uncertain carrier of elemental sulfur in proto-stellar envelopes. This carrier must be subject to destruction in shocks and clearly carry both S and O.

6 Summary

This work presents high angular resolution (~0″.1, 14 au) ALMA observations of the Class I source IRS 44. The continuum emission at 0.87 mm is analyzed, together with molecular species such as SO, SO2, and 34SO2. The main results are summarized below:

  • The continuum emission is contained within a radius of ~0″.2 (30 au) and a total mass (gas + dust) of 4.0 × 10−3 M is calculated for IRS 44. Given that no binary component is detected with our sensitivity, an upper limit of 7 × 10−5 M is estimated for its total dust mass.

  • One SO, six SO2, and two out of three 34SO2 lines are detected; all of the detections show two components in their spectra, a blueshifted one and a redshifted one, both with broad linewidths (between −20 and 20 km s−1). At small scales (<0″.3) the brightest SO2 line is associated with high Eup = 197 K, while the SO2 line with low Eup = 36 K presents the brightest emission at larger angular scales (between 2 and 3″), shows narrow lines below 2 km s−1, and has been associated with a shocked region.

  • Around the protostar, SO2 shows that the redshifted component is more extended than the blueshifted one, likely related with a redshifted protrusion, and the velocity profile is better fitted with an infalling-rotating profile with M* = 1.5 M and rCB = 0″.08. The quiescent shocked region and the red-shifted protrusion seem to be part of a localized streamer, allowing material to fall to the disk-envelope and generate accretion shocks. No evidence of Keplerian motions are found; however, a Keplerian disk is expected inside rCB.

  • The comparison between observed relative intensities of the various lines and RADEX results indicates that the SO2 emission around the protostar arises from a dense region (nH ≥ 108 cm−3) with kinetic temperatures above 90 K. In addition, two SO2 lines are clearly optically thin lines and two others are optically thick lines.

  • The rotational diagram provides kinetic temperatures between 120 and 250 K for the SO2 emitting region, where the warmest regions coincide with the location of the red-shifted protrusion, and SO2 column densities lie between 0.4 and 1.8 × 1017 cm−2. SO column densities are a little lower, between 0.4 and 1.2 × 1017 cm−2, and as a consequence the N(SO2)/N(SO2) ratio lies between 1.0 and 2.0.

  • Optically thick SO2 lines provide Tex values between 70 and 110 K (regarded as lower limits) and a temperature structure consistent with warmer material arising from the south.

  • The high temperatures, compact emission, high nH densities, and moderate velocities agree with the accretion shock scenario, where molecules are being efficiently sublimated from dust grains. We can conclude, therefore, that accretion shocks toward IRS 44 are associated with nH ≥ 108 cm−3, Tkin ≥ 90 K, Trot between 120 and 240 K, Tex ≥ 70 K, SO2 column densities between 0.4 and 1.8 × 1017 cm−2, and velocities between 12 and 14 km s−1.

  • Finally, high-energy SO2 lines (Eup ~ 200 K) seem to be the best tracers of accretion shocks.

Accretion shocks might have important consequences for the chemical content of the disk and the release of neutral species, such as H2O and COMs. It is therefore an important physical process that should be studied in more detail, and high angular resolution observations are essential for this purpose. Future observations of other Class I sources that show bright SO2 emission will be necessary to increase the statistics and achieve a more complete picture of accretion shocks. Other molecular species such as CS, OCS, H2S, and H2CO could provide key additional information. CS is the most abundant sulfur-bearing species in young disks and OCS desorbs from dust grains at a similar temperature than SO, but it does not participate in the gas-phase chemistry below 300 K. H2CO and H2S are key species in the gas-phase formation of SO and SO2. In addition, H2CO is a good tracer of the gas temperature and it has similar desorption temperature to SO2. Finally, a kinematic study of CO isotopologs, in special C18O, could provide information about the existence of a Keplerian disk.

Acknowledgements

We thank the anonymous referee for a number of good suggestions that helped us to improve this work. This paper makes use of the following ALMA data: ADS/JAO.ALMA#2019.1.00362.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. The National Radio Astronomy Observatory is a facility of the National Science Foundation operated under cooperative agreement by Associated Universities, Inc. E.A.dlV. acknowledges financial support provided by FONDECYT grant 3200797. V.G. acknowledges support from FONDECYT Iniciación 11180904, ANID project Basal AFB-170002, and ANID – Millennium Science Initiative Program – NCN19_171. J.K.J. acknowledges support from the Independent Research Fund Denmark (grant No. DFF0135-00123B). D.H. is supported by Centre for Informatics and Computation in Astronomy (CICA) and grant number 110J0353I9 from the Ministry of Education of Taiwan. D.H. acknowledges support from the Ministry of Science of Technology of Taiwan through grant number 111B3005191. N.S. is supported by JSPS KAKENHI grant 20H05845 and pioneering project in RIKEN (Evolution of Matter in the Universe). E.v.D. is supported by the European Research Council (ERC) under the European Union’s Horizon 2020 research and innovation program (grant agreement No. 101019751 MOLDISK).

Appendix A SO2, 34SO2, and SO detections

The integrated fluxes of the detected transitions, and upper limits for nondetections, are presented in Table A.1, where a region with r = 0″. 2 and centered on the continuum peak position was chosen. The spectra, moment 0, and moment 1 maps of the detected transitions at small scales are shown in Figs. A.1, A.2, and A.3, respectively. Figure A.4 presents the large-scale emission, where only three SO2 lines were detected.

Table A.1

Integrated fluxes over a velocity range of 60 km s−1 and taking a circular region with r = 0″.2 centered on the continuum peak position.

thumbnail Fig. A.1

Spectra of SO2 (top), 34SO2 (bottom left), and SO (bottom right), integrated over a circular region with r = 0″.2 and centered on the continuum peak position. The systemic velocity corresponds to 3.7 km s−1.

thumbnail Fig. A.2

Small-scale emission. Top: Moment 0 maps of SO2 above 3σ (color scale) and continuum emission (dashed contours). The moment 0 maps were integrated over 60 km s−1 and the continuum contours start at 20σ and follow steps of 80σ. The Eup value of each transition is indicated in the top left corner of each panel and the synthesized beam is shown by the black filled ellipse in the bottom right corner of the first panel. The color scale is the same for the six panels. Bottom: Moment 1 maps of SO2 above 3σ (color scale) and selected values of their respective moment 0 maps (dashed contours). The contours start at 3σ and follow steps of 7σ, with the exception of the last two panels, which follow steps of 3σ. The Ay value of each transition is indicated in the top left corner of each panel; the color scale is the same for the six panels. The adopted systemic velocity is 3.7 km s−1.

thumbnail Fig. A.3

Same as Fig. A.2, but for 34SO2 and SO. The dashed black contours in the moment 1 maps follow steps of 1σ and 3σ for 34SO2 and SO, respectively. The synthesized beam is shown by the black filled ellipse in the bottom right corner of the first panel.

thumbnail Fig. A.4

Large-scale emission. Moment 0 maps of SO2 integrated over a velocity range of 2 km s−1 and above 1σ. The white star shows the position of the source and the synthesized beam is indicated by the black filled ellipse in the bottom right corner of the first panel. The Eup value of each transition is indicated in the top left corner of each panel and the color scale is the same for all three panels.

Appendix B Radiative transfer

The six different SO2 transitions were employed to compare the observed relative intensities with values obtained from RADEX. The observed relative intensities are shown in the upper row of Fig. B.1. The middle and bottom rows show the RADEX results for a H2 number density of 109 and 108 cm−3, respectively, and the contour levels represent the observed values (shown in the top row). An error value of 1σ (employing error propagation) was added to the observed ratios and is represented by the black dashed contours in the RADEX results.

For nH = 108 cm−3 (bottom row of Fig. B.1), the third panel (197/293) and the fifth panel (197/139) do not present an overlapping region; therefore, this density does not reproduce the observed values and nH should be higher than 108 cm−3. On the other hand, for nH = 109 cm−3 (middle row of Fig. B.1), the possible ranges are presented in Fig. B.2 and the possible values consist of Tkin ≥ 90 K and NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ between 8 × 1016 and 8 × 1017 cm−2.

Figure B.3 shows the optical depth of the six SO2 transitions, obtained with RADEX with nH = 109 cm−3, and the possible values discussed above are shown in gray dashed contours. The brightest transitions, those with Eup values of 197 and 199 K, are optically thick lines, while the two weakest ones (Eup of 245 and 139 K) are optically thin lines. Nothing conclusive can be said for those transitions with Eup values of 36 and 293 K.

thumbnail Fig. B.1

Comparison between observed intensity ratios and RADEX models. Top: Observed intensity ratios between SO2 184.14–183.15 (Eup = 197 K) and the other five transitions, above a 3σ level. The synthesized beam is shown by the black filled ellipse in the bottom left corner of the first panel and the yellow star indicates the position of the source. Center: Intensity ratios from RADEX for the same transitions and employing a H2 number density of 109 cm−3. Bottom: Same analysis from RADEX, but employing a H2 number density of 108 cm−3. The black dashed contours indicate error values oiler and the limits of the possible ranges. All values are possible for the panels without black dashed contours.

thumbnail Fig. B.2

Range of possible values for the SO2 column density and the kinetic temperature (white contours) from the overlap of the observed ranges in Fig. B.1 (for a H2 number density of 109 cm−3), indicated in different colors.

thumbnail Fig. B.3

Optical depth for the six SO2 transitions obtained with RADEX with nH = 109 cm−3. The cyan contour represents τ =1 and the dashed gray contours indicate the range of possible values for NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ and Tkin shown in Fig. B.2.

References

  1. Allen, L. E., Myers, P. C., Di Francesco, J., et al. 2002, ApJ, 566, 993 [NASA ADS] [CrossRef] [Google Scholar]
  2. Alves, F. O., Girart, J. M., Caselli, P., et al. 2017, A&A, 603, L3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  3. Arce, H. G., Mardones, D., Corder, S. A., et al. 2013, ApJ, 774, 39 [CrossRef] [Google Scholar]
  4. Artur de la Villarmois, E., Kristensen, L. E., Jørgensen, J. K., et al. 2018, A&A, 614, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  5. Artur de la Villarmois, E., Jørgensen, J. K., Kristensen, L. E., et al. 2019, A&A, 626, A71 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  6. Balança, C., Spielfiedel, A., & Feautrier, N. 2016, MNRAS, 460, 3766 [CrossRef] [Google Scholar]
  7. Bjerkeli, P., Jørgensen, J. K., & Brinch, C. 2016, A&A, 587, A145 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  8. Boogert, A. C. A., Schutte, W. A., Helmich, F. P., Tielens, A. G. G. M., & Wooden, D. H. 1997, A&A, 317, 929 [NASA ADS] [Google Scholar]
  9. Booth, A. S., van der Marel, N., Leemker, M., van Dishoeck, E. F., & Ohashi, S. 2021, A&A, 651, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  10. Braiding, C. R., & Wardle, M. 2012, MNRAS, 422, 261 [NASA ADS] [CrossRef] [Google Scholar]
  11. Brinch, C., Crapsi, A., Jørgensen, J. K., Hogerheijde, M. R., & Hill, T. 2007, A&A, 475, 915 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  12. Cánovas, H., Cantero, C., Cieza, L., et al. 2019, A&A, 626, A80 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  13. Charnley, S. B. 1997, ApJ, 481, 396 [NASA ADS] [CrossRef] [Google Scholar]
  14. Duchêne, G., Bontemps, S., Bouvier, J., et al. 2007, A&A, 476, 229 [Google Scholar]
  15. Dullemond, C. P., Küffmeier, M., Goicovic, F., et al. 2019, A&A, 628, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  16. Dunham, M. M., Vorobyov, E. I., & Arce, H. G. 2014, MNRAS, 444, 887 [NASA ADS] [CrossRef] [Google Scholar]
  17. Evans, Neal J.I., Rawlings, J. M. C., Shirley, Y. L., & Mundy, L. G. 2001, ApJ, 557, 193 [NASA ADS] [CrossRef] [Google Scholar]
  18. Evans, N. J.II, Dunham, M. M., Jørgensen, J. K., et al. 2009, ApJS, 181, 321 [NASA ADS] [CrossRef] [Google Scholar]
  19. Galli, D., Walmsley, M., & Gonçalves, J. 2002, A&A, 394, 275 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  20. Garufi, A., Podio, L., Codella, C., et al. 2022, A&A, 658, A104 [CrossRef] [EDP Sciences] [Google Scholar]
  21. Goicoechea, J. R., Pety, J., Cuadrado, S., et al. 2016, Nature, 537, 207 [Google Scholar]
  22. Goldsmith, P. F., & Langer, W. D. 1999, ApJ, 517, 209 [Google Scholar]
  23. Hanawa, T., Sakai, N., & Yamamoto, S. 2022, ApJ, 932, 122 [NASA ADS] [CrossRef] [Google Scholar]
  24. Harsono, D., Jørgensen, J. K., van Dishoeck, E. F., et al. 2014, A&A, 562, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  25. Harsono, D., Bjerkeli, P., van der Wiel, M. H. D., et al. 2018, Nat. Astron., 2, 646 [Google Scholar]
  26. Harsono, D., van der Wiel, M. H. D., Bjerkeli, P., et al. 2021, A&A, 646, A72 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  27. Hartmann, L. 1998, Cambridge Astrophysics Series (Canada: Abebooks), 32 [Google Scholar]
  28. Hirano, S., & Machida, M. N. 2019, MNRAS, 485, 4667 [Google Scholar]
  29. Hull, C. L. H., Plambeck, R. L., Kwon, W., et al. 2014, ApJS, 213, 13 [NASA ADS] [CrossRef] [Google Scholar]
  30. Jørgensen, J. K., Johnstone, D., Kirk, H., et al. 2008, ApJ, 683, 822 [CrossRef] [Google Scholar]
  31. Jørgensen, J. K., van Dishoeck, E. F., Visser, R., et al. 2009, A&A, 507, 861 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  32. Karska, A., Kaufman, M. J., Kristensen, L. E., et al. 2018, ApJS, 235, 30 [NASA ADS] [CrossRef] [Google Scholar]
  33. Krasnopolsky, R., & Königl, A. 2002, ApJ, 580, 987 [NASA ADS] [CrossRef] [Google Scholar]
  34. Kristensen, L. E., van Dishoeck, E. F., Benz, A. O., et al. 2013, A&A, 557, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  35. Kuffmeier, M., Goicovic, F. G., & Dullemond, C. P. 2020, A&A, 633, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  36. Lee, C.-F., Mundy, L. G., Reipurth, B., Ostriker, E. C., & Stone, J. M. 2000, ApJ, 542, 925 [NASA ADS] [CrossRef] [Google Scholar]
  37. Lee, C.-F., Hirano, N., Zhang, Q., et al. 2014, ApJ, 786, 114 [NASA ADS] [CrossRef] [Google Scholar]
  38. Lee, J. W. Y., Hull, C. L. H., & Offner, S. S. R. 2017, ApJ, 834, 201 [NASA ADS] [CrossRef] [Google Scholar]
  39. Li, Z.-Y., Krasnopolsky, R., & Shang, H. 2011, ApJ, 738, 180 [Google Scholar]
  40. Lindberg, J. E., Jørgensen, J. K., Brinch, C., et al. 2014, A&A, 566, A74 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  41. Lindberg, J. E., Charnley, S. B., Jørgensen, J. K., Cordiner, M. A., & Bjerkeli, P. 2017, ApJ, 835, 3 [NASA ADS] [CrossRef] [Google Scholar]
  42. McClure, M. K., Furlan, E., Manoj, P., et al. 2010, ApJS, 188, 75 [NASA ADS] [CrossRef] [Google Scholar]
  43. McMullin, J. P., Waters, B., Schiebel, D., Young, W., & Golap, K. 2007 Astronomical Data Analysis Software and Systems XVI, eds. R.A. Shaw, F. Hill, & D.J. Bell, Astronomical Society of the Pacific Conference Series, 376, 127 [NASA ADS] [Google Scholar]
  44. Miura, H., Yamamoto, T., Nomura, H., et al. 2017, ApJ, 839, 47 [NASA ADS] [CrossRef] [Google Scholar]
  45. Nazari, P., Tabone, B., Rosotti, G. P., et al. 2022, A&A, 663, A58 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  46. Oya, Y., Sakai, N., Sakai, T., et al. 2014, ApJ, 795, 152 [NASA ADS] [CrossRef] [Google Scholar]
  47. Oya, Y., Sakai, N., Watanabe, Y., et al. 2018, ApJ, 863, 72 [NASA ADS] [CrossRef] [Google Scholar]
  48. Panoglou, D., Cabrit, S., Pineau Des Forêts, G., et al. 2012, A&A, 538, A2 [CrossRef] [EDP Sciences] [Google Scholar]
  49. Penteado, E. M., Walsh, C., & Cuppen, H. M. 2017, ApJ, 844, 71 [Google Scholar]
  50. Persson, M. V., Jørgensen, J. K., & van Dishoeck, E. F. 2012, A&A, 541, A39 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  51. Robitaille, T. P., Whitney, B. A., Indebetouw, R., Wood, K., & Denzmore, P. 2006, ApJS, 167, 256 [Google Scholar]
  52. Sadavoy, S. I., Stephens, I. W., Myers, P. C., et al. 2019, ApJS, 245, 2 [Google Scholar]
  53. Sai, J., Ohashi, N., Saigo, K., et al. 2020, ApJ, 893, 51 [NASA ADS] [CrossRef] [Google Scholar]
  54. Sakai, N., Sakai, T., Hirota, T., et al. 2014, Nature, 507, 78 [Google Scholar]
  55. Sakai, N., Oya, Y., López-Sepulcre, A., et al. 2016, ApJ, 820, L34 [NASA ADS] [CrossRef] [Google Scholar]
  56. Seifried, D., Sánchez-Monge, Á., Walch, S., & Banerjee, R. 2016, MNRAS, 459, 1892 [NASA ADS] [CrossRef] [Google Scholar]
  57. Shu, F., Najita, J., Galli, D., Ostriker, E., & Lizano, S. 1993, in Protostars and Planets III, eds. E.H. Levy, & J.I. Lunine (Tucson: University of Arizona Press), 3 [Google Scholar]
  58. Stahler, S. W., Korycansky, D. G., Brothers, M. J., & Touma, J. 1994, ApJ, 431, 341 [NASA ADS] [CrossRef] [Google Scholar]
  59. Suutarinen, A. N., Kristensen, L. E., Mottram, J. C., Fraser, H. J., & van Dishoeck, E. F. 2014, MNRAS, 440, 1844 [Google Scholar]
  60. Tabone, B., Cabrit, S., Bianchi, E., et al. 2017, A&A, 607, L6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  61. Tafalla, M., Santiago-García, J., Hacar, A., & Bachiller, R. 2010, A&A, 522, A91 [CrossRef] [EDP Sciences] [Google Scholar]
  62. Terebey, S., Shu, F. H., & Cassen, P. 1984, ApJ, 286, 529 [Google Scholar]
  63. Terebey, S., Vogel, S. N., & Myers, P. C. 1992, ApJ, 390, 181 [NASA ADS] [CrossRef] [Google Scholar]
  64. Terebey, S., van Buren, D., Hancock, T., et al. 2001, ASP Conf. Ser., 243, 243 [NASA ADS] [Google Scholar]
  65. Tychoniec, Ł., Manara, C. F., Rosotti, G. P., et al. 2020, A&A, 640, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  66. van der Marel, N., Kristensen, L. E., Visser, R., et al. 2013, A&A, 556, A76 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  67. van der Tak, F. F. S., Black, J. H., Schöier, F. L., Jansen, D. J., & van Dishoeck, E. F. 2007, A&A, 468, 627 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  68. van Gelder, M. L., Tabone, B., van Dishoeck, E. F., & Godard, B. 2021, A&A, 653, A159 [CrossRef] [EDP Sciences] [Google Scholar]
  69. van Gelder, M. L., Nazari, P., Tabone, B., et al. 2022, A&A, 662, A67 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
  70. Wilson, T. L. 1999, Rep. Prog. Phys., 62, 143 [Google Scholar]
  71. Woods, P. M., Occhiogrosso, A., Viti, S., et al. 2015, MNRAS, 450, 1256 [Google Scholar]
  72. Yen, H.-W., Koch, P. M., Takakuwa, S., et al. 2015, ApJ, 799, 193 [NASA ADS] [CrossRef] [Google Scholar]
  73. Yorke, H. W., & Bodenheimer, P. 1999, ApJ, 525, 330 [NASA ADS] [CrossRef] [Google Scholar]
  74. Young, E. T., Lada, C. J., & Wilking, B. A. 1986, ApJ, 304, L45 [NASA ADS] [CrossRef] [Google Scholar]
  75. Zasowski, G., Kemper, F., Watson, D. M., et al. 2009, ApJ, 694, 459 [CrossRef] [Google Scholar]

All Tables

Table 1

Spectral setup and parameters of the observed molecular transitions.

Table A.1

Integrated fluxes over a velocity range of 60 km s−1 and taking a circular region with r = 0″.2 centered on the continuum peak position.

All Figures

thumbnail Fig. 1

Continuum emission (0.87 mm) toward IRS 44 above 5σ (σ = 0.08 mJy beam−1). The white contours represent the weakest emission of [5, 10, 15, 20, and 25σ] for clarity. The black star shows the position of the continuum peak and the synthesized beam is indicated by the white filled ellipse. The magenta ellipse represents the deconvolved size from the 2D Gaussian fit.

In the text
thumbnail Fig. 2

Emission of the SO2 184,14–183,15 line. Left: spectrum rebinned by a factor of 4, integrated over a circular region with r = 0″.2, and centered on the continuum peak position. Center, moment 0 map above 3σ (color scale), integrated over 60 km s−1, and continuum emission (black contours), starting at 20σ and following steps of 40σ. Right: moment 1 map (color scale) and moment 0 map (white contours) above 3σ. The blue and red arrows show the direction of the outflow from van der Marel et al. (2013), the yellow star indicates the continuum peak position, and the synthesized beam is represented by the black filled ellipse in the bottom right corner. The adopted systemic velocity is 3.7 km s−1.

In the text
thumbnail Fig. 3

Contour maps of SO2 184,14–183,15 for different velocity ranges, shifted to 0 velocity. The contours start at 5σ and follow steps of 5σ. The panel that includes the systemic velocity (3.7 km s−1) is shown by the green contours in the first panel, while blue- and redshifted emission is represented by the blue and red contours, respectively. The yellow star shows the position of the source and the synthesized beam is indicated by the black filled ellipse in the bottom left corner of the first panel.

In the text
thumbnail Fig. 4

Emission of the SO2 53,3–42,2 line. Left: spectra integrated over a circular region with r = 0″.2, centered on the continuum peak position (gray) and centered on the offset region (red). The spectrum taken at the continuum peak position has been rebinned by a factor of 4. Right: moment 0 map above 3σ integrated over 60 km s−1. The gray and red circles represent the regions from which the spectra in the left panel were taken. The blue and red arrows show the direction of the outflow, the yellow star indicates the continuum peak position, and the synthesized beam is represented by the black filled ellipse in the bottom right corner. The adopted systemic velocity is 3.7 km s−1.

In the text
thumbnail Fig. 5

Position-velocity diagram for SO2 184,14–183,15. Left: emission above 3σ% employing a PA of 157°. The blue and red dots represent blue- and redshifted emission peaks above ±3 km s−1, respectively. The black line represents an infalling-rotating profile with M* = 1.5 M, an inclination of 70°, and rCB = 0708 (magenta dashed lines). The rCR is shown by the green dashed lies. Center: zoomed-in version for the redshifted emission. Right: zoomed-in version for the blueshifted emission. The systemic velocity is 3.7 km s−1.

In the text
thumbnail Fig. 6

Temperature structure of IRS 44. Left: rotational diagram at the source position, where only the optically thin SO2 transitions (blue dots) and the 34SO2 lines (red dots) are used for the fit. The abundance ratio 32S/34S = 22 is from Wilson (1999). Optically thick SO2 transitions (green dots) and those lines where the optical depth is not conclusive (open dots) show a significant offset with respect to the other lines, and they were not included in the calculation of the rotational temperature and SO2 column densities. Center, rotational temperature map created from optically thin SO2 and 34SO2 transitions. Contours show specific values of 140, 160, 180, and 200 K. Right: moment 1 map (right panel of Fig. 2) with the same specific values of the rotational temperature as those of the central panel, showing that the highest temperatures coincide with the redshifted protrusion. The synthesized beam is shown by the black ellipse in the upper left corner and the yellow star indicates the position of the source.

In the text
thumbnail Fig. 7

Column densities. Left: SO2 column density obtained from the rotational diagram method. The synthesized beam is shown by the black ellipse in the upper left corner. Center: SO column density, assuming optically thin emission and employing the rotational temperatures from Fig. 6. Right: column density ratio of SO2 to SO. The yellow star indicates the position of the source.

In the text
thumbnail Fig. 8

Excitation temperature of the optically thick SO2 184,14–183,15 transition using Eq. (5). The black contour represents the specific value of 90 K. The yellow star indicates the position of the source and the synthesized beam is shown by the black ellipse in the upper left corner.

In the text
thumbnail Fig. 9

Schematic representation of IRS 44. The quiescent SO2 emission (knot) would be part of a streamer that allows material to fall into the disk-envelope through the redshifted protrusion and generate accretion shocks. These shocks sublimate SO2 molecules from dust grains and enhance the SO2 gaseous abundance, showing an emission peak ~0″.1 to the south of the protostar.

In the text
thumbnail Fig. A.1

Spectra of SO2 (top), 34SO2 (bottom left), and SO (bottom right), integrated over a circular region with r = 0″.2 and centered on the continuum peak position. The systemic velocity corresponds to 3.7 km s−1.

In the text
thumbnail Fig. A.2

Small-scale emission. Top: Moment 0 maps of SO2 above 3σ (color scale) and continuum emission (dashed contours). The moment 0 maps were integrated over 60 km s−1 and the continuum contours start at 20σ and follow steps of 80σ. The Eup value of each transition is indicated in the top left corner of each panel and the synthesized beam is shown by the black filled ellipse in the bottom right corner of the first panel. The color scale is the same for the six panels. Bottom: Moment 1 maps of SO2 above 3σ (color scale) and selected values of their respective moment 0 maps (dashed contours). The contours start at 3σ and follow steps of 7σ, with the exception of the last two panels, which follow steps of 3σ. The Ay value of each transition is indicated in the top left corner of each panel; the color scale is the same for the six panels. The adopted systemic velocity is 3.7 km s−1.

In the text
thumbnail Fig. A.3

Same as Fig. A.2, but for 34SO2 and SO. The dashed black contours in the moment 1 maps follow steps of 1σ and 3σ for 34SO2 and SO, respectively. The synthesized beam is shown by the black filled ellipse in the bottom right corner of the first panel.

In the text
thumbnail Fig. A.4

Large-scale emission. Moment 0 maps of SO2 integrated over a velocity range of 2 km s−1 and above 1σ. The white star shows the position of the source and the synthesized beam is indicated by the black filled ellipse in the bottom right corner of the first panel. The Eup value of each transition is indicated in the top left corner of each panel and the color scale is the same for all three panels.

In the text
thumbnail Fig. B.1

Comparison between observed intensity ratios and RADEX models. Top: Observed intensity ratios between SO2 184.14–183.15 (Eup = 197 K) and the other five transitions, above a 3σ level. The synthesized beam is shown by the black filled ellipse in the bottom left corner of the first panel and the yellow star indicates the position of the source. Center: Intensity ratios from RADEX for the same transitions and employing a H2 number density of 109 cm−3. Bottom: Same analysis from RADEX, but employing a H2 number density of 108 cm−3. The black dashed contours indicate error values oiler and the limits of the possible ranges. All values are possible for the panels without black dashed contours.

In the text
thumbnail Fig. B.2

Range of possible values for the SO2 column density and the kinetic temperature (white contours) from the overlap of the observed ranges in Fig. B.1 (for a H2 number density of 109 cm−3), indicated in different colors.

In the text
thumbnail Fig. B.3

Optical depth for the six SO2 transitions obtained with RADEX with nH = 109 cm−3. The cyan contour represents τ =1 and the dashed gray contours indicate the range of possible values for NSO2${N_{{\rm{S}}{{\rm{O}}_2}}}$ and Tkin shown in Fig. B.2.

In the text

Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.

Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.

Initial download of the metrics may take a while.