Multiple chemical tracers finally unveil the intricate NGC\,1333 IRAS\,4A outflow system. FAUST XVI

The exploration of outflows in protobinary systems presents a challenging yet crucial endeavour, offering valuable insights into the dynamic interplay between protostars and their evolution. In this study, we examine the morphology and dynamics of jets and outflows within the IRAS\,4A protobinary system. This analysis is based on ALMA observations of SiO(5--4), H$_2$CO(3$_{0,3}$--2$_{0,3}$), and HDCO(4$_{1,4}$--3$_{1,3}$) with a spatial resolution of $\sim$150\,au. Leveraging an astrochemical approach involving the use of diverse tracers beyond traditional ones has enabled the identification of novel features and a comprehensive understanding of the broader outflow dynamics. Our analysis reveals the presence of two jets in the redshifted emission, emanating from IRAS\,4A1 and IRAS\,4A2, respectively. Furthermore, we identify four distinct outflows in the region for the first time, with each protostar, 4A1 and 4A2, contributing to two of them. We characterise the morphology and orientation of each outflow, challenging previous suggestions of bends in their trajectories. The outflow cavities of IRAS\,4A1 exhibit extensions of 10$''$ and 13$''$ with position angles (PA) of 0$^{\circ}$ and -12$^{\circ}$, respectively, while those of IRAS\,4A2 are more extended, spanning 18$''$ and 25$''$ with PAs of 29$^{\circ}$ and 26$^{\circ}$. We propose that the misalignment of the cavities is due to a jet precession in each protostar, a notion supported by the observation that the more extended cavities of the same source exhibit lower velocities, indicating they may stem from older ejection events.


INTRODUCTION
It has long been known that the gravitational collapse, at the basis of the formation of a star and planetary system, is inevitably accompanied by a simultaneous ejection process, mediated by the magnetic field (e.g.Frank et al. 2014;Bally et al. 2014;Lee 2020).These ejections form the so-called "outflow systems", in general, constituted by supersonic jets, bow-shocks and cavities, all more or less delineated and identified by shocked material (see, e.g., Fig. 3 of Bally 2016;Rivera-Ortiz et al. 2023).Several molecular lines are generally used for these analyses.Particularly useful are lines from species that have a very low abundance in molecular clouds because they are trapped in the interstellar grains, both in the refractory cores and volatile mantles.The paradigmatic species is SiO, considered a clear-cut shock hallmark (e.g.Schilke et al. 1997;Caselli et al. 1997;Gusdorf et al. 2008).Formaldehyde is also a good tracer of material ejected from the grain mantles, while deuterated formaldehyde traces the most recent ejected material (e.g.Fontani et al. 2014).
Often, but not always, the outflows have a bipolar structure.Outflows have been the focus of a multitude of studies at the beginning because they were unexpectedly observed towards protostars before any sign of collapse was found (e.g.Kwan & Scoville 1976;Bally & Lada 1983;Frank et al. 2014;Bally 2016).Soon after, it was realised the huge importance of the outflows in getting rid of the angular momentum of the protostellar infalling matter, as well as in shaping and disrupting the protostellar envelope and molecular cloud hosting the protostar.Finally, molecular outflows cause different shocks, which in one way or another increase the gas chemical complexity (e.g.Bachiller 1996) and are therefore very precious chemical laboratories (e.g.Codella et al. 2017).
In addition to all the above reasons, the study of outflows in protobinary systems provides a fascinating glimpse into the complexity of star birth.It sheds light on the genesis and early dynamical interactions of multiple protostars, providing essential insights into their intertwined evolution.This is particularly true for systems where the ★ E-mail: layal.chahine@univ-grenoble-alpes.fr† E-mail: cecilia.ceccarelli@univ-grenoble-alpes.frbinary companions are close enough to possibly interact and/or influence each other.In fact, how the presence of a companion affects the ejection phenomenon is still not elucidated, especially in very young protobinary systems, where the surrounding envelope can mask and likely influence the evolution of the outflow itself.It is not by chance that the best-studied systems, especially when it comes to chemistry, refer to single outflows (e.g.HH 212; Codella et al. 2007;Lee et al. 2015, L1157-mm Bachiller et al. 1995;Arce et al. 2013;Lefloch et al. 2017;Codella et al. 2017 andHH 30 Pety et al. 2006;Louvet et al. 2018 among others).
A specific system that has piqued the interest of many researchers due to its distinctive features and complex dynamics and chemistry is IRAS 4A.IRAS 4A is a binary system, harbouring two protostars IRAS 4A1 and IRAS 4A2 (hereafter 4A1 and 4A2 respectively), located in the star-forming region NGC 1333 within the Perseus molecular cloud, at a distance of 293 ± 22 pc (Zucker et al. 2018).Notably, the system stands out for its remarkable outflow phenomena, which have been the subject of numerous studies.Indeed, the outflows within the IRAS 4A system have been mapped using various molecular tracers, such as CO (e.g.Blake et al. 1995;Girart et al. 1999;Jørgensen et al. 2007;Yıldız et al. 2012;Santangelo et al. 2015;Ching et al. 2016;Podio et al. 2021;Chuang et al. 2021), CS (Di Francesco et al. 2001;Taquet et al. 2020), SiO (e.g.Lefloch et al. 1998;Choi 2001;Santangelo et al. 2015;Koumpia et al. 2017;De Simone et al. 2022), SO and other sulphur-bearing species (e.g.Santangelo et al. 2015;Taquet et al. 2020;Podio et al. 2021;Chuang et al. 2021), as well as HCN (e.g.Choi 2001), and H 2 CO (e.g.Di Francesco et al. 2001;Su et al. 2019).The presence of more complex molecules (CH 3 OH and CH 3 CHO; De Simone et al. 2022) was revealed along the outflows, highlighting the rich chemistry of this region.
The SiO observations at an angular resolution of 2 ′′ (Choi 2005;Choi et al. 2006) suggested the presence of a bipolar northeastsouthwest outflow driven by 4A2, and a tentative south outflow possibly driven by 4A1, with no distinct northern counterpart being evident.To determine the outflow direction, Choi (2005) employed a least-squares fit method based on the positions of compact peaks near the driving source, identifying a trajectory passing near 4A2.Similar results were also found with other tracers.Using a higher angular resolution ∼1 ′′ , Santangelo et al. (2015) successfully distinguished the two outflows in the southern blue-shifted emission and, notably, at extremely high velocities (41 -55 km/s), they could resolve the redshifted counterpart of the 4A1 southern outflow for the first time.The same feature was also observed at similar velocity ranges with other tracers.There is no agreement, however, over the direction of the redshifted counterpart from 4A1.While the SiO emission mapped by Santangelo et al. (2015) suggests that the jet is oriented north-south (N-S), Ching et al. ( 2016) considered a position angle (PA) of -9 • , aligning with the trajectory of the southern counterpart of the jet.Nonetheless, Ching et al. ( 2016) have suggested that the northern outflow of IRAS 4A1 is later bent due to the interaction with the outflow of IRAS 4A2.A recent study using CO and SO at an angular resolution of 0.3 ′′ suggested the presence of two northern outflow lobes, one from each protostar, both bent towards the northeast (Chuang et al. 2021).Another study at similar angular resolution using H 2 CO suggested that both outflows are oriented north-south (Su et al. 2019).Meanwhile, the morphology of the southern blueshifted lobes is similar to that seen in previous studies, with the IRAS 4A1 lobe oriented southeast and the IRAS 4A2 lobe directed southwest.It is also worth noting that an abrupt bend at 20 ′′ from 4A1 was also observed in several studies (e.g.Choi 2005;Yıldız et al. 2012;De Simone et al. 2022).Choi (2005) referred to this as "directional variability" and proposed that this feature may be attributed to external medium perturbations, specifically a deflection resulting from interactions with an external core or cloud, distinct from the one associated with the IRAS 4A system.
Despite numerous studies on this system, the intricate morphology of the outflows has posed challenges in definitively attributing them to the respective protostars.Consequently, a comprehensive understanding of the broader scenario remains elusive.Achieving a more coherent picture demands advancements in several aspects, including higher spatial and spectral resolution, enhanced sensitivity, wider field of view coverage, and the use of diverse tracers.In this paper, we have examined the morphology and kinematics of the jets and outflows towards the IRAS 4A system, via the study of the line emission of three standard shock tracers (see above): SiO, H 2 CO and HDCO.We employ these three tracers because a posteriori and after a fast check of the maps for this source we found that they offer complementary insights into the observed dynamics, with H 2 CO and HDCO playing a crucial role in identifying the cavities and enriching our understanding of the system's dynamics.We aimed to conduct a comprehensive characterisation, pinpointing and associating distinct features with their respective protostars, using observations from the Large Program FAUST (Fifty AU STudy of the chemistry in the disc/envelope system of solar-like protostars; 2018.1.01205.L, PI: S. Yamamoto; Codella et al. 2021).The present article is structured as follows.In Sect. 2 we describe the observations.In Sect. 3 we present the molecular line maps, together with the main results of the analysis.In Sect. 4 we discuss the structure of the outflows in the IRAS 4A protobinary system, as derived by the FAUST observations, and compare it with the interpretations in previous literature.In Sect. 5 we summarise the conclusions of this work.

OBSERVATIONS
The IRAS 4A system was observed at 1.2 mm with ALMA during its Cycle 6 operations, between October 2018 and September 2019, as part of the Large Program FAUST.The observations were performed in Band 6 using the 12-m array (C43-3 and C43-6 configurations, including 46 and 50 antennas, respectively) and the 7-m Atacama Compact Array (ACA).The baselines for the 12-m array were between 15.1 m and 3.6 km, probing angular scales from 0. ′′ 03 (∼9 au) to 8. ′′ 2 (∼2400 au).
The quasar J0237+2848 was used for bandpass and flux calibration, while J0336+3218 and J0328+3139 were used for phase calibration.The absolute flux calibration uncertainty is estimated to be <15%.The data calibration was performed using the standard ALMA calibration pipeline with the Common Astronomy Software Applications package 5.6.1-8(CASA 1 , CASA Team et al. 2022).An additional calibration routine2 has been used to correct for the   and for spectral data normalisation.The spectrum of IRAS 4A is particularly rich and the identification of line-free continuum channels had to be done by hand, outside the ALMA pipeline.Line-free channels were then averaged per-spectral window and used to create a continuum model for each ALMA configuration, and both amplitude and phase self-calibration were employed to correct for remaining tropospheric phase fluctuations, position offsets, and amplitude calibration differences between the execution blocks.For each ALMA configuration, two rounds of phase-only selfcal were performed, followed by two rounds of amplitude and phase selfcal once the model was sufficiently complete.The resulting improvement in the continuum dynamic range was typically a factor of 5-7.The continuum model was then subtracted from the line channels to produce continuum-subtracted visibilities.A similar self-calibration technique was then used to align data across configurations.The visibility data from the 12-m arrays were then combined to produce the continuum-subtracted line cubes.The resulting continuumsubtracted cubes were cleaned using Högbom algorithm (Högbom 1974) and natural weighting in the IRAM-GILDAS software package3 .The data analysis was performed using the same package.The resulting beam size and rms are summarised in Table 1.The FWHM Field of View (FoV) is 29 ′′ for SiO, 28.9 ′′ for H 2 CO, and 25.5 ′′  () Frequencies and spectroscopic parameters have been extracted from Müller et al. (2013) for SiO, Bocquet et al. (1996); Müller & Lewen (2017) for H 2 CO, and Bocquet et al. (1999) for HDCO.
for HDCO.The maps shown in the paper are not corrected for the primary beam attenuation, as we are more interested in the structures' morphology (rather than the flux) that was not affected by the primary beam attenuation.By not applying correction for primary beam attenuation, we ensure that features extending beyond the field of view, are not inadvertently excluded in our analysis.This approach allows us to preserve the integrity of the observed structures, facilitating a comprehensive understanding of their spatial distribution and morphology.

RESULTS
We have imaged the emission from SiO, H 2 CO, and HDCO towards the IRAS 4A system.The comprehensive understanding of the observed dynamics originates from the combined insights derived from the analysis of the three tracers.For a concise overview of our findings, readers are directed to Section 3.4.

SiO emission
Figure 1 illustrates a thorough view of the SiO emission in six distinct panels, depicting the spatial and velocity features of the studied region.The upper left panel shows the SiO(5-4) integrated intensity map (moment 0), the lower left panel shows the SiO velocity map (moment 1) both integrated between -27 and 35 km s −1 , and the remaining four panels show the integrated channel maps, each representing SiO emission at distinct velocity intervals.The SiO moment 0 map reveals different features, including bow-shocks, jets, and outflow, alongside the finger-like structures previously reported by De Simone et al. (2022) in the south.The SiO emission in the figure spans from -27 km s −1 (V= -34 km s −1 ; V lsr = 7 km s −1 ) to 35 km s −1 (V= +28 km s −1 ).The narrow spectral window limits detection beyond 35 km s −1 , however, the emission is expected to be relatively low beyond this velocity threshold.The moment 1 map indicates that features located north of the protostars are characterised by redshifted velocities, while those situated south exhibit blueshifted velocities.The integrated channel maps in the figure correspond to four different velocity intervals, namely the high-velocity 1 (HV1; V= -34 (+28) -±20 km s −1 ), high-velocity 2 (HV2; V= ±12 -±20 km s −1 ), intermediate-velocity (IV; V= ±4 -±12 km s −1 ), and low-velocity intervals (LV; V= ±1 -±4 km s −1 ).

High Velocity (HV1 and HV2)
HV1: In the redshifted HV1 regime, a north-south jet (J1) stemming from source 4A1 is evident.The jet shows bullet-like structures and reaches almost ∼5 ′′ (∼1500 au) from the protostar.At distances of 12 ′′ -18 ′′ (∼3500-5300 au) from 4A2, another jet (J2) is visible, clearly directed towards it, and distinguished by the presence of identifiable bullets along this feature.A zoom-in onto the innermost region showing of J1 from J2 is depicted in Fig. A1.
On the other hand, in the blueshifted emission, we detect a ringlike feature, directed almost N-S, associated with the first blueshifted outflow cavity of 4A1 (C1-blue).HV2: In the panel corresponding to the HV2 interval, the emission starts to overlap at the north.Conversely, in the south, the components attributed to 4A1 can be distinguished from those of 4A2.At ∼13 ′′ (∼3800 au) to the south of 4A2, we identify a prominent bowshock structure associated with an outflow cavity of 4A2 (C2-blue).Additionally, there is observable emission around ∼23 ′′ (∼6700 au) to the south of 4A2, which is likely associated with another cavity (C4-blue; See Sect.3.3).

Intermediate Velocity (IV)
Transitioning to the redshifted intermediate velocity (IV) regime, shown in the bottom left panel of Fig. 1, the SiO emission presents a complex overlap of the structures originating from 4A1 and 4A2, making interpretation difficult.These features are better discerned in HDCO and will be discussed further (See Sect.3.2).However, with the detection of a bow-shock at ∼18 ′′ (∼5300 au) from 4A2 (Bow2), which is also perpendicular to the trajectory of its jet 4A2 (J2), we observe that the extended northeast-southwest outflow cavity wall (C2-w) reaches Bow2.This observation suggests that we are tracing the high-velocity component of the 4A2 outflow (C2-red; further explanation is provided in the subsequent sections).The outflow shows an S-shape bending within 4 ′′ (∼1200 au) from IRAS 4A2 that was also observed in previous studies (Santangelo et al. 2015;Chuang et al. 2021).It is worth noting that Bow2 was present in the CS maps of Taquet et al. (2020), but was not explicitly identified by the authors.Additionally, an arc-like feature appears at approximately 20-25 ′′ (∼5900-7300 au) from 4A2 (C4-w).This feature is coincident with what was interpreted as a sharp bend (also referred to as deflection or directional variability) by Choi (2005), which has also been observed by Yıldız et al. (2012); De Simone et al. (2022), as illustrated in Fig.

2.
Initially, this structure was attributed to a deflection of the 4A2 outflow caused by a dense core (Choi 2005).However, we suggest that it is rather linked with the LV component of the 4A2 outflow (C4red; see Sect.3.3), being the remnant of its western cavity wall.The emission of this outflow is better discerned with H 2 CO.Within the same velocity interval, the blue-shifted emission shows the presence of two bow-shocks directed towards 4A1, one in the N-S direction reaching 10 ′′ (∼3000 au) from the protostar and the other with a PA of ∼-12 • reaching ∼13 ′′ (∼3800) from it.Although less clear in SiO, the bow-shocks are better observed in the other tracers, and they illustrate two outflow cavities driven by 4A1 (C1-blue and C3blue).To the south of 4A2, distinguishing the emission becomes more intricate, and we likely observe two outflows extending ∼13 ′′ (C2-blue) and ∼23 ′′ (C4-blue) from the protostar.The 4A2 outflow resembles that of HH 46/47 (Arce et al. 2013) showing a bifurcation along its eastern wall.5-4) towards the NGC 1333 IRAS 4A system.In the moment 0 and moment 1 panels the SiO emission is integrated between -27 and 35 km s −1 and a 5 threshold cut was applied 5 ( = 0.81 mJy beam −1 km s −1 ).The channel maps correspond to four different velocity intervals: High-Velocity 1 (HV1), High-Velocity 2 (HV2), Intermediate-Velocity (IV) and Low-Velocity (LV).The intervals and their corresponding labels are shown at the upper right corners of each channel maps panel.The contour levels are [5,25,40,60,100]  for the redshifted HV1 regime and at [5,25,50,100]  for the other regimes with  = (1.3mJy beam −1 km s −1 , 1.5 mJy beam −1 km s −1 , 1.7 mJy beam −1 km s −1 , 0.9 mJy beam −1 km s −1 , respectively).The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The jets from each source are marked with arrows.The synthesised beam is depicted in the lower-left corner of the top-left panel.The FWHM of the Field of View is depicted with a dashed grey circle in the lower-left panel.

Low Velocity (LV)
In the redshifted low-velocity regime, the emission is dominated by the components linked with 4A2.Indeed, prominent outflow cavity walls directed towards 4A2 are distinctly visible towards the west and the east (i.e.C2-w, C2-e and C4-w) together with the Bow2 bow-shock.We suggest that C2-w, C2-e and Bow2 delineate the HV component of the 4A2 outflow (C2-red) and that C4-w is linked with the LV component of the 4A2 outflow cavity wall (C4-red; See Sect.

for further explanation)
Lastly, it is noteworthy that in the moment 0 map of SiO, a bow-like feature is discernible to the east of the 4A1 blueshifted outflow (C1blue).This structure was previously identified by De Simone et al. (2022) at lower angular resolution, where it resembled a finger-like structure similar to the one observed at the west of the A2 outflow.The authors proposed these fingers to be formed by a train of shocks due to an expanding bubble crashing against NGC 1333.The higher resolution of the ALMA data reveals that the finger-like structure at the west of the 4A1 outflow has an arc-like shape at its end.However, since there is no apparent association with any specific outflow emanating from 4A1 or 4A2, we postpone a dedicated investigation of these fingers to future work.

HDCO emission
In Fig. 3 we present the moment 0 map of HDCO integrated between (-3.6 and 20.1 km s −1 ).HDCO emission delineates the outflow cavities towards the IRAS 4A system, helping to better discern different components, particularly to the north of the protostar, despite some overlap in emission.Notably, HDCO emission reveals, for the first time, a bow-like feature at ∼10 ′′ from 4A1 (Bow1) perpendicular to  its N-S jet.Stronger emission contours also show two curved features (C1-w and C2-w), one aligned with 4A1 and another with 4A2.For a detailed inspection, we show the channel map of the HDCO redshifted emission towards the north (V = 8.4-17.8km s −1 ; Fig. A2).We note that C1-w is emanating from A1 and has a higher velocity compared with the other features (including C2-w).This feature (C1-w) extends up to ∼10 ′′ from 4A1 reaching the bowshock Bow1, suggesting that Bow1 and C1-w delineate the 4A1 cavity, oriented in the N-S direction (C1-red).Conversely, C2-w emission extends beyond 10 ′′ (see also Fig. A2), forming a V-shape emerging from 4A2, together with C2-e (detected in SiO).The overlap between the SiO and HDCO emission (Fig. 4), shows that these two features and the Bow2 bow-shock observed in SiO are actually tracing the same outflow cavity, which we named C2-red, where C2-w is the western cavity wall and C2-e is the eastern.Further, emission beyond 10 ′′ from the protostars aligns with the western wall of C2-red, suggesting an association with A2.Notably, while C2-w and C2-e belong to the same cavity, C2-w exhibits a higher velocity by almost 3 km s −1 , indicating possible cavity rotation.Moving southward, an extended cavity from 4A1 with PA∼-12 • east (C3-blue) and part of the eastern wall of the cavity from 4A2 can be observed.Intriguingly, within the 4A1 cavity, another "nested" cavity, C1-blue, directed almost N-S, is noticeable (see Sect. 3.3 for further discussion).It is worth noting that while C1-blue has C1-red as a redshifted counterpart, C3-blue lacks a counterpart, suggesting it could be exiting the cloud or hidden by the emission of the 4A2 cavity.

H 2 CO emission
Figure 5 presents the maps of the IRAS 4A system using H 2 CO as a molecular tracer.The H 2 CO emission is illustrated in six panels, the upper left panel shows the H 2 CO moment 0 map, and the remaining five panels show the integrated channel maps.H 2 CO effectively maps the distribution of the bulk gas within the IRAS 4A system; it also traces the intricate outflow cavities originating from both the 4A1 and 4A2 sources, as well as the jet from 4A1.The emission in the figure spans from -16 km s −1 (V= -23 km s −1 ; V lsr = 7 km s −1 ) to 35 km s −1 (V= +28 km s −1 ).The channel maps velocity intervals correspond to five different velocity intervals, the HV2 and IV are as in Fig. 1, while for the HV1 regime (V= -23 (+28) -±20 km s −1 ).To better illustrate the features of the IV regime, we split it into IV1 and IV2.We note that the SiO map encompasses a slightly broader range of velocities compared to H 2 CO, which is indicative of the latter exhibiting dimmer emissions that remain undetected at the highest velocities.

High Velocity (HV1 and HV2)
HV1: In the HV1 panel, the H 2 CO emission displays the presence of the previously identified north-south jet originating from 4A1 (J1).HV2: In the redshifted HV2 interval, we detected the previously observed curved feature C1-w from 4A1 (i.e. the western wall of the redshifted cavity), consistent with our HDCO observations.Additionally, we identified another feature oriented northeast-southwest, pointing towards 4A2 and aligned with the western wall of C2-red (seen in SiO).Transitioning to the blueshifted regime, we identify a distinct N-S outflow cavity from 4A1 (C1-blue).

Intermediate Velocities (IV1, IV2)
Moving further into the blueshifted IV1 regime, as in SiO, clear differentiation emerges between the distinct outflow cavities, attributed to 4A1 and 4A2.Within the 4A1 cavity across both IV1 and IV2 regimes, akin to HDCO observations, we observe two bow-like features.One is nearly perpendicular to the N-S jet direction (C1-blue), while the other is tilted about 12 • eastward (C3-blue), affirming the presence of dual cavities from this source (see also the sketch in Fig. 6).Additionally, the cavities manifest at different velocities; C1-blue is observable starting from -12.5 km/s, whereas the bowshock of C3-blue is discernible at -6.9 km/s, with the cavity becoming more apparent at lower velocities (See Fig A3).To the south of 4A2, at least one outflow is observed.
The interpretation of the emission in the redshifted regime presents considerable complexity, as in the other tracers.The 4A1 outflow cavity (C1-red) is detected.The western lobe (C1-w) extends up to 10 ′′ from 4A1, while we see only 3 ′′ from the eastern lobe.The features delineating C2-red (i.e.C2-w, Bow2 and C2-e) are also detected.As in HDCO, the western cavity wall C2-w is observed at higher velocities, being 7 km s −1 faster than eastern C2-e (as shown in the detailed channel map A4) suggesting a possible cavity rotation.

Low Velocity (LV)
The redshifted LV regime is primarily characterised by the outflows originating from 4A2, notably the extended northeast-southwest cavity (C2-red) detected in SiO.Alongside this cavity, denoted as C2red, we identify another structure labelled "C4-w," also observed in SiO, situated to the west of 4A2.This feature, stemming from IRAS 4A2, extends approximately ∼25 ′′ from the protostar, coinciding with the distinct directional variability previously documented by Choi (2005); Yıldız et al. (2012); De Simone et al. (2022).With H 2 CO, C4-w is detected between ∼8 and 12 km s −1 (see the detailed channel map in Fig. A4).Within this velocity range, we observe emergence of emission towards the northeast of Bow 2 and the east of C2-e (Fig. A4).The overlap of the H 2 CO LV contours with the HDCO moment 0 shows that the emission to the east of C2-e aligns  , and [5, 50, 100]  for the others, with  = 0.9 mJy beam −1 km s −1 , 1.1 mJy beam −1 km s −1 , 1.2 mJy beam −1 km s −1 , 1.0 mJy beam −1 km s −1 , 1.5 mJy beam −1 km s −1 , respectively).The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The solid ellipses in the last two panels represent the C2-red and C4-red cavities.The synthesised beam is depicted in the lower-left corner of the top-left panel.
with the HDCO feature C4-e, extending approximately ∼15 ′′ from 4A2 (see Fig. 4), suggesting that we are tracing the same feature.We propose that C4-w and C4-e together form an elliptical shape with a PA of ∼26 • , directed towards 4A2, implying they trace a second outflow system driven by 4A2 (C4-red), where C4-w represents the western wall and C4-e is the eastern.Hence, the structure previously thought to be a sharp bend (directional variability; Choi 2005;Yıldız et al. 2012;De Simone et al. 2022) is actually part of the western cavity wall of C4-red.The cavity is not entirely detected, likely due to filtering.In Fig. 4, we depict the "potential" blue counterpart of C4-red with a purple dashed ellipse.This ellipse aligns with the SiO emission observed at 23 ′′ from 4A2, indicating that it may represent a portion of the counterpart of C4-red.However, it is challenging to differentiate it clearly from C2-blue, and its western wall is showing very weak emission (see moment 0 map of H 2 CO).

Overview of the outflows and associated molecules
To simplify and present a concise overview, we have summarised the results in the sketch of Fig. 6 and listed the found components and their major properties in Table 2.A comprehensive figure is also shown in the appendix (Fig. A5, wherein all the features identified in this work are overlaid on the molecular integrated velocity maps of SiO, H 2 CO and HDCO, with corresponding labels.Overall, our observations reveal the presence of two outflows emanating from 4A1 and another two outflows emanating from 4A2.As summarised in the introduction, this is the first time that four outflows are identified in the region: each protostar, 4A1 and 4A2, is responsible for two of them.This conclusion is built on the various components detected in our maps in different lines at different velocities, as briefly listed here to guide the reader through their detailed description in the following sections.

Jets J1 and J2:
The high-velocity maps of SiO and H 2 CO reveal the presence of two red-shifted jets, J1 and J2, north of 4A1 and 4A2, emanating from 4A1 and 4A2, respectively.

Bow-shocks Bow1 and Bow2:
The lower velocity maps of H 2 CO and HDCO show a red-shifted bow-shock, Bow1, associated with J1, while those of SiO and H 2 CO show another bow-shock, Bow2, associated with J2.
A1 Cavities C1-red; C1-blue; C3-blue: Despite the clear detection of only one jet, J1, emanating from 4A1, our SiO, H 2 CO and HDCO maps show the presence of one red-shifted cavity in the north, C1-red, and two blue-shifted cavities in the south, C1-blue and C3-blue.C1red and C1-blue are aligned almost North-South, indicating that they are part of the same bipolar outflow system.Likewise, the alignment of the jet J1, the bow-shock Bow1 and the cavity C1-red suggest that they all belong to the same ejection event that created them, which reaches a distance from 4A1 of about 10 ′′ .
A second cavity, labelled C3-blue, extending to the south of 4A1, is distinctly visible, particularly in the HDCO emission.This cavity exhibits an eastward tilt in comparison to C1-blue.Notably, it lacks a corresponding red-shifted counterpart to the north.We suggest that the absence of a redshifted counterpart might be influenced by the overlap with the redshifted cavity of 4A2.This overlap could potentially obscure the cavity associated with 4A1, or even interact with it, dispersing its emission.As for the prospects of detection, it is plausible that the redshifted counterpart would be observable at larger distances, beyond the influence of the 4A2 cavity's obscuration, however, more observations would be required to investigate this point.

A2 Cavities C2-red, C2-blue, C4-red and C4-blue:
We identify two red-shifted cavities emanating from 4A2 in the north, C2-red, C4red, and we associate them with two blue-shifted cavities in the south, C2-blue and C4-blue, although the separation is less clear.The C2red is visible in the intermediate-velocity maps of SiO, H 2 CO and HDCO, while the C4-red appears in the low-velocity maps of SiO and H 2 CO.Both cavities, C2-red and C4-red, have an aligned blueshifted counterpart in the south, C2-blue and C4-blue, respectively.The C2-red is aligned with the jet emanating from 4A2, J2, with a position angle (PA) of ∼29 • .While C4-red is slightly tilted with a PA of ∼26 • .Finally, we notice that, in our interpretation, the emission observed in the SiO north of 4A2, previously interpreted as part of a deflected jet from 4A2 by (Choi 2005), is, in fact, part of the cavity C4-red.
Overall comments: The structures mentioned earlier are the outcome of assembling evidence that occasionally manifests in distinct tracers and at different velocities.Notably, components of the same cavity may not manifest simultaneously, particularly the eastern segments.

4A1 outflow systems:
Our observations disclose the existence of two discernible outflow systems (C1-red & C1-blue, and C3-blue) originating from 4A1.The first (C1-red & C1-blue) is directed N-S, spanning a distance up to 10 ′′ both to the north and south of the protostar.However, the detection of the red-shifted eastern wall is somewhat limited, with only 3 ′′ from it being observed.The cavity walls observed in this system are delineated by SiO, H 2 CO, and HDCO, exhibiting a velocity of approximately ∼18 km s −1 in red.
The second outflow system (C3-blue), driven by 4A1, is more extended, reaching nearly 15 ′′ to the south of 4A1.Intriguingly, the counterpart of the 4A1 southeast cavity is absent in the northern region, implying its possible escape from the cloud, as in the case of HH 46/47 (Arce et al. 2013).The blue-shifted lobe moves at a velocity of ∼15 km s −1 .This outflow system is tilted by ∼12 • with respect to the N-S outflow system of 4A1, suggesting a possible jet precession.

4A2 outflow systems:
Our observations unravel the intricate details of the 4A2 outflow systems, where two distinct outflow systems (C2-red & C2-blue, and C4-red & C4-blue) emerge, as detailed in the previous sections.The first system (C2-red & C2-blue) extends approximately 18 ′′ both to the north and south of 4A2, though the eastern lobe of the southern outflow appears less distinct.This outflow system has a PA of ∼29 • .The cavity walls move at a velocity of ∼16 km s −1 in the red, and they are observed with the three molecular tracers: SiO, H 2 CO, and HDCO.Notably, we observe a velocity gradient in the red-shifted cavity walls, with the west wall exhibiting a velocity of approximately 3-7 km s −1 red-shifted compared to the east wall.This intriguing observation hints at a potential rotation of the cavity, adding a layer of complexity to the dynamics of the 4A2 outflow system.Rotating outflow cavities have been detected in other objects (e.g.HH 212 Codella et al. 2014, VLA1623 Ohashi et al. 2022).We expect that if the jet is launched by a magnetic centrifugal mechanism it should rotate consistently with the disk (e.g.HH 212 Lee et al. 2017) and transfer the rotational velocity also to the cavities.
The second outflow system stemming from 4A2 (C4-red & C4blue) exhibits a greater extent, reaching almost 25 ′′ both to the north and 23 ′′ to the south, though with a less distinct boundary in the southern direction.This system has a PA of ∼26 • .The cavity walls move at a velocity of approximately 10 km s −1 in SiO and H 2 CO, while HDCO does not exhibit these features.The absence of features in HDCO could stem from a sensitivity issue (not connected though to the excitation of the line), particularly considering that C4-red is the oldest cavity.It is worth noting that HDCO might also be frozen back onto the ice, further complicating its detection.A comprehensive astrochemical modelling could provide better insights into this matter.Consequently, the two cavities of 4A2 show a slight misalignment of about 3 • , indicating a possible precession or rotation of the ejection.Notably, the cavity of the second, more extended outflow system displays lower velocities, aligning with the concept of an older ejection event that has decelerated in the current epoch.
As mentioned in Sect.3.4, we associated the SiO emission labelled "C4-w" in Fig. 6 with the cavity wall of the second outflow system rather than a bend induced by the interaction of the jet with an obstructing dense core (Choi 2005).The SiO(1-0) emission observed by Choi (2005) are similar to that of SiO (2-1) observed more recently by De Simone et al. (2022), and to compare with our results we use the latter for convenience.
In Fig. 2, we illustrate the SiO map observed by De Simone et al. (2020a) alongside those from our study, overlaid with two solid ellipses the two ellipses representing the cavities of the aforementioned two outflow systems.While the second ellipse encompasses "C4-w", it does not reach the farthest SiO blob observed by Choi (2005); De Simone et al. (2022) at ∼35 ′′ from 4A2, prompting a discussion on its nature.
One can speculate that a third event gave rise to a third outflow system, whose bow-shock coincides with the farthest SiO blob observed at ∼35 ′′ from 4A2.We represent it with a dashed ellipse in Fig. 2.This ellipse would be inclined by 31 • , slightly deviating from the orientation of the other two ellipses (also see "Hypothetical cavity" in Fig. 6).The sequence from the most distant ellipse (corresponding to the oldest ejection event) to the closest (and youngest event) reads as 31 • , 26 • , and 29 • , suggesting once more a precession or rotation of the ejections.This hypothesis is further supported by the morphology of the blue outflow lobe of 4A2, which bears a notable resemblance to the red lobe of the HH 46/47 outflow, displaying a V-shape and a subsequent bifurcation in the cavity lobe.As proposed by Arce et al. (2013), this feature could result from multiple mass ejections that were not detected in our observations.
It is worth noting that the SiO blob at ∼35 ′′ from 4A2 observed by Choi (2005); De Simone et al. (2022) is ∼6 ′′ out of our FoV and was not detected in our observations due to filtering.

The comparison between 4A1 and 4A2 ejections
Our data shows that each of the protostars had at least two different episodic events.However, it is worth noting that the outflow cavities of 4A2, particularly C4-red, shows a much larger extension and is detected at lower velocities.This would suggest that 4A2 is older than 4A1, as previously suggested (e.g.Choi et al. 2007;Choi et al. 2011a;Santangelo et al. 2015;De Simone et al. 2020a).2020), (g) Chuang et al. (2021).* These values represent the location of the feature from the corresponding protostar, rather than the extension.

Comparison with previous studies
Outflowing motion: As highlighted in the Introduction, extensive studies have focused on the outflows in the IRAS 4A system.While the blue-shifted emission from the southern region has been more clearly discernible when it comes to associating the observed feature with one of the protostars, the interpretation of the emission to the north of the protostars has proven to be challenging.This difficulty arises from the complex overlap between components originating from 4A1 and 4A2, hindering a clear understanding of the underlying processes in this region.
In the context of the jets within the IRAS 4A system, our observations have unveiled the previously identified northern jet originating from 4A1 (Santangelo et al. 2015;Taquet et al. 2020).The enhanced angular resolution of our data has provided a clearer view of this jet, allowing us to resolve intricate bullet structures within it.Notably, our observations have also introduced a new element to the understanding of the system: a SiO component associated with the primary jet of 4A2.This component also exhibits bullet-like structures.
Regarding the outflowing motion to the north of the protostars, previous studies presented different hypotheses regarding its origin.Some proposed it solely from the 4A2 outflow (Choi 2005), while others attributed it to the bent outflow of 4A1 (Ching et al. 2016).Alternatively, certain studies argued for the presence of two extended bent cavities-one from 4A1 and another from 4A2 (Chuang et al. 2021).Our observations suggest a scenario distinct from these interpretations, as summarised in Table 2.We contend that neither of the cavities is bent.In fact, the non-detection of the bow-shock Bow1 in the previous studies has led to uncertainties regarding the extent and orientation of the cavity associated with 4A1, where it was thought that it is extended up to ∼17 ′′ and bent to the east.Our observations propose that the cavity is oriented N-S and does not extend farther than 10 ′′ from 4A1.It is also somewhat encompassed within that of 4A2, at least in projection.This N-S outflow cavity of 4A1 has a blue-shifted counterpart to the south of the protostar, where a second cavity (C3-blue) is also detected.The latter is misaligned with the former by ∼12 • .The blueshifted cavities are both present in the SO maps of Chuang et al. (2021), but only C3-blue was identified.
In 4A2, previous studies suggested the existence of an extended outflow with a PA of approximately 18.9 • , that changes direction around 20 ′′ (i.e., directional variability) attributed to an interaction with an obstructing dense core (Choi 2005).Recent higher-resolution observations have confirmed the presence of an outflow cavity with a PA ranging from 19 • to 22 • (e.g.Ching et al. 2016;Chuang et al. 2021).However, Marvel et al. (2008) noted that the position angle of the water maser outflow in IRAS 4A differs from the position angles associated with the larger-scale outflows originating from the two sources within the system.Our current data reveal the existence of a high-velocity outflow cavity originating from 4A2, aligning closely with the jet observed in our dataset with a PA of ∼29 • .Furthermore, we distinguish between LV and HV outflow cavities in 4A2, clarifying that what was previously thought to be "directional variability" of the 4A2 outflow (C4-w) is simply an extension of the LV component of the western cavity wall.Jet precession: The misalignment between the outflow systems in 4A1 and those in 4A2, suggests that the jets of the two sources are actually precessing.To verify that the misalignment is attributed to a precession rather than the potential shifting position of protostars over time, we examined the presence of proper motions of this source.We found no significant proper motion, with the coordinates showing only minor fluctuations within the uncertainties since their initial determination (Looney et al. 2000;Reipurth et al. 2002), as compared to recent derivations (De Simone et al. 2020b).While the possibility of precession in the 4A1 jet has not been proposed before, previous studies hinted at the precession scenario in 4A2 to explain specific features observed in the IRAS 4A system (e.g.Choi et al. 2006;Santangelo et al. 2015;Chuang et al. 2021).Choi et al. (2006) referred to the mechanism as "drifting" rather than precession due to a lack of evidence for periodic change.In another work, Choi et al. (2011b) have suggested that the 4A2 jet is rotating within a cone of open-ing half-angle ∼3.1 • .The cause of the 4A2 jet's precession remains unclear.Several studies suggested that the precession might result from a close unresolved binary system, with a separation of 30-80 au (Choi et al. 2006), 32 au (Marvel et al. 2008) or larger than ∼12 au (Chuang et al. 2021).Alternatively, Chuang et al. (2021) proposed that a misalignment between the magnetic field and the cloud's angular momentum vector may induce precession.This scenario may also apply to IRAS 4A1.Particularly noteworthy is the resemblance of the blueshifted outflow cavities to the misaligned outflows observed in HOPS 373, where one cavity is "nested" within the other (Lee et al. 2024).The authors propose that the observed misaligned outflows may originate from either a single protostar, following a reorientation of its angular momentum axis, or from a close protostellar binary, with each outflow linked to a distinct binary component.They further suggest that, in the first scenario, the reorientation could result from a change in the angular momentum axis of the accreting gas or from a misalignment between the rotation axis of the core and the orientation of the global magnetic field (e.g.Machida et al. 2020;Okoda et al. 2021).With our observations, it is challenging to determine the origin of the precession of 4A1 and 4A2' jets, as the identification of a potential binary companion for each protostar would require an angular resolution of about 12 au.A reorientation of the system could be plausible given the misalignment of the magnetic field in the IRAS 4A envelope with the rotation axis of the circumbinary envelope (Chuang et al. 2021 and references therein).Further observations with higher angular resolution are needed to resolve the discs and investigate the possibility of close binaries for each protostar.

Effect of stellar multiplicity:
The influence of stellar multiplicity on outflow activity during the protostellar stage remains inadequately characterised.One notable effect is the precession of the disc, which can alter the orientation and morphology of the outflows.Precession refers to the rotation of the jet axis within a cone, leading to changes in the direction of the outflow.This process was observed towards several sources and was linked to the presence of a binary companion (e.g.HOPS-370 Tobin et al. 2019b;Chahine et al. 2022), or a close unresolved binary (e.g.L1157 Podio et al. 2016, CepE-A de A. Schutzer et al. 2022, HOPS-317 in the HH 24 complex Reipurth et al. 2023, HOPS-373SW Lee et al. 2024) Additionally, binary and multiple systems may exhibit jet wiggling, where the orientation of the jet remains constant, but the spatial position of the launching point varies as the stars orbit each other (e.g.HH 211 Lee et al. 2010, HH 212 Lee et al. 2017, CepE-A de A. Schutzer et al. 2022, VLA15 in OMC-2 FIR 4 Chahine et al. 2022, SSV63 and Jet C in the HH 24 complex Reipurth et al. 2023).Furthermore, the disturbance of the disc due to the presence of a companion can lead to increased accretion onto the protostars, potentially driving strong outbursts of outflow activity (e.g.HH 175 Reipurth & Friberg 2021).While precession has been observed in the IRAS 4A system, potentially linked to stellar multiplicity, further investigations are warranted, particularly at the scale of the circumstellar discs.This additional scrutiny is necessary to discern any signs of disc disturbance or misaligned accretion onto the protostars.

CONCLUSIONS
In this study, we delved into the morphology and kinematics of the jets and outflows towards the IRAS 4A system using three distinct tracers: SiO, H 2 CO, and HDCO, at a spatial resolution of ∼150 au.The use of these three tracers proved to be indispensable in revealing new features, piecing together the intricate puzzle, and gaining a comprehensive understanding of the larger framework, as each tracer contributed a crucial component to the analysis.Our main conclusions are summarised below.
(i) We have successfully identified two distinct jets in the redshifted emission, originating from 4A1 and 4A2.The jet from 4A1 is oriented in a N-S direction, while that from 4A2, observed for the first time, exhibits a NE-SW orientation with a PA of 29 • .
(ii) Each protostar has two outflow systems, with 4A1 showcasing three outflow cavities and 4A2 presenting four.Contrary to prior studies, we argue that none of the outflows is bent.
(iii) We observe an outflow cavity to the north of IRAS 4A1 and two cavities to its south.One of these southern cavities aligns with a N-S direction, while the other exhibits a PA of -12 • without a corresponding counterpart in the north.The misalignment of these cavities suggests the possibility of a precession occurring within the 4A1 jet.
(iv) The analysis of the outflows originating from IRAS 4A2 reveals the presence of two cavities in both the northern and southern directions, albeit with less distinct emission in the south.The cavities are slightly misaligned, and the more extended is moving at lower velocities, suggesting it is potentially stemming from an older ejection event.In addition, the cavity moving at higher velocities shows an indication of rotation (v) Our study challenges previous interpretations of the directional variability observed in the outflow of IRAS 4A2, suggesting instead that it represents the remnants of the western wall of the lowvelocity outflow cavity.We speculate on the occurrence of a third event that was not detected in our data.
In summary, our study underscores the importance of utilising diverse tracers beyond traditional ones, emphasising the significance of comprehensive studies in unravelling the complexities of outflow systems.

APPENDIX A: MAPS
In Fig. A1 we show a zoomed view of the SiO emission within the HV1 velocity interval.In Figs.A2 and A4 we report the channels maps of HDCO and H 2 CO in the redshifted LV and IV regimes.In Fig. A3 we present the detailed blueshifted channel map of H 2 CO.In Fig. A5 we show the integrated intensity maps of the different tracers used to probe the emission towards IRAS 4A, with superimposed ellipses and arrows representing all the features identified in this work.
This paper has been typeset from a T E X/L A T E X file prepared by the author.

Figure 1 .
Figure1.Integrated intensity (moment 0; top left) map, velocity map (moment 1; bottom left) and channel maps (right panels) of SiO(5-4) towards the NGC 1333 IRAS 4A system.In the moment 0 and moment 1 panels the SiO emission is integrated between -27 and 35 km s −1 and a 5 threshold cut was applied 5 ( = 0.81 mJy beam −1 km s −1 ).The channel maps correspond to four different velocity intervals: High-Velocity 1 (HV1), High-Velocity 2 (HV2), Intermediate-Velocity (IV) and Low-Velocity (LV).The intervals and their corresponding labels are shown at the upper right corners of each channel maps panel.The contour levels are[5, 25, 40, 60, 100]  for the redshifted HV1 regime and at[5, 25, 50, 100]   for the other regimes with  = (1.3mJy beam −1 km s −1 , 1.5 mJy beam −1 km s −1 , 1.7 mJy beam −1 km s −1 , 0.9 mJy beam −1 km s −1 , respectively).The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The jets from each source are marked with arrows.The synthesised beam is depicted in the lower-left corner of the top-left panel.The FWHM of the Field of View is depicted with a dashed grey circle in the lower-left panel.

Figure 2 .
Figure 2. SiO(2-1) integrated intensity contours from De Simone et al. (2022) with superimposed SiO(5-4) integrated intensity map from this work in colours.The outflow cavities emanating from IRAS 4A2 and identified in this work are represented with solid ellipses.The possible third event ("Hypothetical cavity") emanating from IRAS 4A2 is represented with a dashed ellipse.The FWHM of the Field of View in our observations is depicted with a dotted magenta circle.The synthesised beams are depicted in the lower-left corner.

Figure 3 .
Figure 3. Moment 0 map of HDCO integrated between -3.6 and 20.1 km s −1 .A 3 threshold cut was applied ( = 0.40 mJy beam −1 km s −1 ).The contour levels are at 30  with  = 0.40 mJy beam −1 km s −1 .The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The features mentioned in the text are labelled in grey.The solid ellipses represent the IRAS 4A1 outflow cavities observed in HDCO.The synthesised beam is depicted in the lower-right corner.

Figure 4 .
Figure 4. Left: Moment 0 map of HDCO in colour scale with superimposed SiO contours at 5 ( = 0.81 mJy beam −1 km s −1 ).The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The HV component of the 4A2 outflow is represented by a navy blue ellipse.The synthesised beams are depicted in the lower-right corner.Right: Moment 0 map of HDCO in colour scale with superimposed H 2 CO 5 and 30 redshifted LV contours in coral ( = 1.5 mJy beam −1 km s −1 ).The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The LV component of the 4A2 outflow is represented by a purple ellipse.The synthesised beams are depicted in the lower-right corner.

Figure 6 .
Figure 6.Sketch illustrating the comprehensive perspective of the jets and outflows towards the IRAS 4A system as proposed in this study.The left diagram incorporates the annotations of the jets and various cavities proposed in this work.The right diagram highlights the annotations of the distinct features of the redshifted cavities, including the cavity walls (east "e" and west "w") and the bow shocks ("Bow"), as discussed and referenced in the text.

"
Chemical Origins" (PI: L. Podio).GS acknowledges the INAF-Minigrant 2023 TRIESTE ("TRacing the chemIcal hEritage of our originS: from proTostars to planEts"; PI: G. Sabatini).ClCo and LP also acknowledge the financial support under the National Recovery and Resilience Plan (NRRP), Mission 4, Component 2, Investment 1.1, Call for tender No. 104 published on 2.2.2022 by the Italian Ministry of University and Research (MUR), funded by the European Union -NextGenerationEU-Project Title 2022JC2Y93 ChemicalOrigins: linking the fossil composition of the Solar System with the chemistry of protoplanetary disks -CUP J53D23001600006 -Grant Assignment Decree No. 962 adopted on 30.06.2023 by the Italian Ministry of Ministry of University and Research (MUR).EB acknowledges support from the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany´s Excellence Strategy -EXC 2094 -390783311.M.B. acknowledges support from the European Research Council (ERC) Advanced Grant MOPPEX 833460.SBC was supported by the NASA Planetary Science Division Internal Scientist Funding Program through the Fundamental Laboratory Research work package (FLaRe).

J2 J1 Figure A1 .
Figure A1.Zoomed view of the SiO emission within the HV1 velocity interval.The contours are as in the HV1 panel of Fig. 1.The positions of IRAS 4A1 and IRAS 4A2 are marked with stars.The jets from each source are marked with arrows.

Table 1 .
Spectral parameters of the SiO, H 2 CO, and HDCO lines observed towards IRAS 4A.

Table 2 .
(Su et al. 2019components of the outflows emanating from IRAS 4A1 and IRAS 4A2 and detected in the SiO, H 2 CO and HDCO maps presented in this work.The systemic velocity of the system is V  = 7 km s −1(Su et al. 2019).