Context. Accretion and ejection dictate the outcomes of star and planet formation processes. The midinfrared MIR wavelength range offers key tracers of processes that have been difficult to detect and spatially resolve in protostars until now. BR  Aims We aim to characterize the interplay between accretion and ejection in the lowmass Class I protobinary system TMC1 comprising two young stellar objects TMC1W and TMC1E at a 85 au separation. BR  Methods Using the James Webb Space Telescope JWST MidInfrared Instrument MIRI observations in 528 m range we measured the intensities of emission lines of HSUB2SUB atoms and ions for instance the Fe II and Ne II and HI recombination lines. We analyzed the spatial distribution of the different species using the MIRI Medium Resolution Spectrometer MRS capabilities to spatially resolve emission at 0.20.7 scales. we compared these results with the corresponding Atacama Large Millimetersubmillimeter Array ALMA maps tracing cold gas and dust. BR  Results We detected HSUB2SUB outflow coming from TMC1E with no significant HSUB2SUB emission from TMC1W. The HSUB2SUB emission from TMC1E outflow appears narrow and extends to wider opening angles with decreasing ESUBupSUB from S8 to S1 rotational transitions indicating the disk wind as its origin. The outflow from TMC1E protostar shows spatially extended emission lines of Ne II Ne III Ar II and Ar III with their line ratios consistent with UV radiation as a source of ionization. With ALMA we detected an accretion streamer infalling from gt 1000 au scales onto the TMC1E component. The TMC1W protostar powers a collimated jet detected with Fe II and Ni II making it consistent with energetic flow. A much weaker ionized jet is observed from TMC1E and both jets appear strikingly parallel to each other indicating that the disks are coplanar. TMC1W is associated with strong emission from hydrogen recombination lines tracing the accretion onto the young star. BR  Conclusions MIRIMRS observations provide an unprecedented view of protostellar accretion and ejection processes on 20 au scales. Observations of a binary Class I protostellar system show that the two processes are clearly intertwined with accretion from the envelope onto the disk influencing a wideangle wind ejected on disk scales. Finally the accretion from the protostellar disk onto the protostar is associated with the source launching a collimated highvelocity jet within the innermost regions of the disk.
['accretion', 'accretino disks', 'techniques: spectroscopic', 'stars: low-mass', 'stars: protostars', 'ISM: jets and outflows', 'infrared: ISM', 'Astrophysics - Solar and Stellar Astrophysics', 'Astrophysics - Earth and Planetary Astrophysics', 'Astrophysics - Astrophysics of Galaxies']
JWST Observations of Young protoStars JOYS. Linked accretion and ejection in a Class I protobinary system
{'Linked accretion and ejection in a Class I protobinary system': "Łukasz Tychoniec 1 , 2 , Martijn L. van Gelder 1 , Ewine F. van Dishoeck 1 , 3 , Logan Francis 1 , Will R. M. Rocha 1 , 4 , Alessio Caratti o Garatti 5 , Henrik Beuther 6 , Caroline Gieser 3 , Kay Justtanont 7 , Harold Linnartz ⋆ 1 , 4 , Valentin J. M. Le Gouellec 8 , 9 , Giulia Perotti 6 , R. Devaraj 10 , Benoît Tabone 11 , Thomas P. Ray 10 , Nashanty G. C. Brunken 1 , Yuan Chen 1 , Patrick J. Kavanagh 12 , Pamela Klaassen 13 , Katerina Slavicinska 1 , 4 , Manuel Güdel 6 , 14 , 15 , Goran Östlin 16 .  \n- 1 Leiden Observatory, Leiden University, PO Box 9513, 2300RA, Leiden, The Netherlands e-mail: tychoniec@strw.leidenuniv.nl\n- 2 European Southern Observatory, Karl-Schwarzschild-Strasse 2, 85748 Garching bei München, Germany\n- 3 Max-Planck-Institut für Extraterrestrische Physik, Giessenbachstrasse 1, D-85748 Garching, Germany\n- 4 Laboratory for Astrophysics, Leiden Observatory, Leiden University, PO Box 9513, NL 2300 RA Leiden, The Netherlands\n- 5 INAF-Osservatorio Astronomico di Capodimonte, Salita Moiariello 16, I-80131 Napoli, Italy\n- 6 Max Planck Institute for Astronomy, Königstuhl 17, 69117 Heidelberg, Germany\n- 7 Department of Space, Earth and Environment, Chalmers University of Technology, Onsala Space Observatory, 439 92 Onsala, Sweden\n- 8 NASA Ames Research Center, Space Science and Astrobiology Division M.S.245-6 Mo ff ett Field, CA94035, USA\n- 9 NASA Postdoctoral Program Fellow\n- 10 School of Cosmic Physics, Dublin Institute for Advanced Studies, 31 Fitzwilliam Place, Dublin 2, Ireland\n- 11 Université Paris-Saclay, CNRS, Institut d'Astrophysique Spatiale, 91405, Orsay,France\n- 12 Department of Experimental Physics, Maynooth University, Maynooth, Co Kildare, Ireland\n- 13 UK Astronomy Technology Centre, Royal Observatory Edinburgh, Blackford Hill, Edinburgh EH9 3HJ, UK\n- 14 Dept. of Astrophysics, University of Vienna, Türkenschanzstr. 17, A-1180 Vienna, Austria\n- 15 ETH Zürich, Institute for Particle Physics and Astrophysics, Wolfgang-Pauli-Str. 27, 8093 Zürich, Switzerland\n- 16 Department of Astronomy, Oskar Klein Centre, Stockholm University, 106 91 Stockholm, Sweden  \nReceived: December 11, 2023", 'ABSTRACT': "Context. Accretion and ejection dictate the outcomes of star and planet formation processes. The mid-infrared (MIR) wavelength range o ff ers key tracers of processes that have been di ffi cult to detect and spatially resolve in protostars until now.  \nAims. We aim to characterize the interplay between accretion and ejection in the low-mass Class I protobinary system TMC1, comprising two young stellar objects: TMC1-W and TMC1-E at a 85 au separation.  \nMethods. Using the James Webb Space Telescope (JWST) Mid-Infrared Instrument (MIRI) observations in 5 - 28 µ m range, we measured the intensities of emission lines of H2, atoms, and ions, for instance, the [Fe II] and [Ne II], and HI recombination lines. We analyzed the spatial distribution of the di ff erent species using the MIRI Medium Resolution Spectrometer (MRS) capabilities to spatially resolve emission at 0 '' . 2-0 '' . 7 scales. We compared these results with the corresponding Atacama Large Millimeter / submillimeter Array (ALMA) maps tracing cold gas and dust.  \nResults. We detected H2 outflow coming from TMC1-E, with no significant H2 emission from TMC1-W. The H2 emission from TMC1-E outflow appears narrow and extends to wider opening angles with decreasing E up from S(8) to S(1) rotational transitions, indicating the disk wind as its origin. The outflow from TMC1-E protostar shows spatially extended emission lines of [Ne II], [Ne III], [Ar II], and [Ar III], with their line ratios consistent with UV radiation as a source of ionization. With ALMA, we detected an accretion streamer infalling from > 1000 au scaling onto the TMC1-E component. The TMC1-W protostar powers a collimated jet, detected with [Fe II] and [Ni II], making it consistent with energetic flow. A much weaker ionized jet is observed from TMC1-E, and both jets appear strikingly parallel to each other, indicating that the disks are co-planar. TMC1-W is associated with strong emission from hydrogen recombination lines, tracing the accretion onto the young star.  \nConclusions. MIRI-MRS observations provide an unprecedented view of protostellar accretion and ejection processes on 20 au scales. Observations of a binary Class I protostellar system show that the two processes are clearly intertwined, with accretion from the envelope onto the disk influencing a wide-angle wind ejected on disk scales. Finally, the accretion from the protostellar disk onto the protostar is associated with the source launching a collimated high-velocity jet within the innermost regions of the disk.", '1. Introduction': "The interplay between the accretion and ejection of mass sets the outcome of the star formation process (Hartmann et al. 2016). However, the link between those two key processes is di ffi cult  \nto ascertain. It is expected that the excess of angular momentum is lost during the accretion from the envelope to the disk or from the disk to the protostar and is removed by the outflow of material either by the jet (collimated, high-velocity stream of gas) or wide-angle disk wind (with a broader opening angle and lower velocities) For more details we refer to Frank et al.  \n(2014); Pascucci et al. (2022) and references therein. The emerging paradigm of the protostellar outflow launching is that the outflow is magnetohydrodynamically (MHD) driven (Pascucci et al. 2022), with the outflow launching point and radius still a matter of debate. Observational evidence shows wide-angle winds launched from the disk (Bjerkeli et al. 2016; Louvet et al. 2018) and collimated jets launched near or within the dust sublimation radius (Tabone et al. 2017; Lee et al. 2019). However, the balance and interaction between the jets and disk winds have not been well characterized yet.  \nBinary systems are unique laboratories for studying protostellar accretion and ejection since it is possible to compare the properties of the two protostars directly. Moreover, interactions between binary members may have strong e ff ects on the truncation of the disk and triggering of accretion outbursts. Since many low-mass protostars are born in binary systems, capturing their formation and evolution is important for understanding the star formation process. (Raghavan et al. 2010; Tobin et al. 2016, 2020; O ff ner et al. 2023). However, such studies in the mid-infrared (MIR) regime are scarce, limited by the spatial resolution of instruments o ff ered at these wavelengths, while many essential tracers of outflow activity (e.g., H2, [Fe II]) are found in this wavelength range. This regime is crucial, especially for the studies of the youngest protostars, where extinction limits the possibility of studying those sources in optical and near-IR.  \nInfrared emission lines have been studied in protostars and their outflows, especially with Spitzer Infrared Spectrograph (IRS) spectroscopy with a typical aperture of 5 '' (e.g., Maret et al. 2009; Lahuis et al. 2010; Tappe et al. 2012). Spitzer -IRS observations already showed hints that binary systems might present very di ff erent outflow properties such as H2 lines dominating emission from one source and ion lines from the other (Dionatos et al. 2014). However, it has been challenging to map the protostellar outflows and jets in the MIR at sub-arcsecond scales without the integrated field unit (IFU) capabilities to retrieve spectral information from each spatial element on sub-arcsecond scales.  \nAccretion onto the protostar and infall from envelope to disk are typically probed in very di ff erent energy regimes, with stellar accretion primarily studied with hydrogen recombination lines (Alcalá et al. 2017), while infall is probed by kinematic features in the cool gas in the sub-millimeter regime (Jørgensen et al. 2002a; Aso et al. 2015). Therefore, the combination of the two, augmented by MIR diagnostics, can give us a complete understanding of the transfer of the material all the way from the largescale envelope through the young disk onto the protostar.  \nTMC1, also known as IRAS 04381 + 2540, is a nearby Class I protobinary star system located in the Taurus-L1527 region (141.8 ± 1.4 pc; Krolikowski et al. 2021). The system contains two protostars associated with near-IR emission tracing scattered light from the outflow cavity walls and the near-IR jet. (Apai et al. 2005; Terebey et al. 2006). When observed with Atacama Large Millimeter / submillimeter Array (ALMA) at subarcsecond resolution at 1.3 mm, the binary stars were resolved with the projected separation of 85 au (van 't Ho ff et al. 2020).  \nIn this paper, we present observations of the TMC1 protobinary system with Medium Resolution Spectroscopy mode (MRS; Wells et al. 2015; Argyriou et al. 2023) of the MidInfraRed Instrument (MIRI; Wright et al. 2015; Rieke et al. 2015; Wright et al. 2023) on board the James Webb Space Telescope (JWST; Rigby et al. 2023). Complementary ALMA observations already presented in Tychoniec et al. (2021) are used to compare the MIR emission with cold gas tracers.  \nThis paper presents stunning di ff erences in the appearance of outflow in the TMC1 binary system. This work is organized  \nas follows. In Section 2, we introduce JWST / MIRI MRS and ALMA observations; Section 3 describes the results, including spectral and spatial analysis of the system; Section 4 discusses the origin of di ff erences in this system and Section 5 provides our conclusions.", '2.1. JWST-MIRI': "JWST MIRI-MRS observations were obtained on March 4, 2023, as part of the JWST Observations of Young ProtoStars (JOYS) Guaranteed Time Observations (GTO) PID 1290 (PI: E. F. van Dishoeck). Three sub-bands were integrated with 36 groups in the FASTR1 integration mode with two dithers in a negative direction optimized for the extended source. The integration time per sub-band was 200s, so 600s to obtain full spectrum from 4.9 µ m to 28.6 µ m with a spectral resolving power of R 1500-3500. The target position was 4 h 41 m 12 s .7 + 25 · 46 ' 34 '' . 8 (J2000). Integration on the source was preceded by the o ff -source background observation centered on coordinates 4 h 39 m 32. s 3 + 26 · 05 ' 14 '' . 5 (J2000). The background observation was obtained with 72 groups in the FASTR1 integration mode, with no dithering, resulting in the same total integration time as the science target.  \nThe data were reduced with the JWST pipeline v.1.11.2 (Bushouse et al. 2023) using the calibration reference data system (CRDS) context file jwst1100.pmap . The raw telescope data were processed with the Detector 1 pipeline with default settings. Detector calibrated files are constructed in the Spec2 step, including fringe flat for extended sources (Mueller et al. in prep.) and residual fringe correction (Kavanagh et al. in prep.). In this step, the background image is subtracted from the detector image using dedicated background rate files. At this step, a bad-pixel routine is applied to calibrated files using the Vortex Image Processing (VIP) package (Christiaens et al. 2023). Final data cubes for each channel and sub-band are created at Spec3 pipeline step, with the outlier rejection step turned o ff since this is handled by a custom-made bad-pixel routine. The absolute flux calibration uncertainty of MIRI-MRS is estimated from flight performance at 5.6 ± 0.7 % (Argyriou et al. 2023).", '2.2. ALMA': "For the analysis, we also used archival ALMA data. CO J = 21 230.538 GHz data and SO J = 56-45 219.9494 GHz with a synthesized beam (resolution element) of 0 '' . 41 × 0 '' . 29 are from project 2017.1.01350.S (PI: Tychoniec). CO 2-1 and SO 56-45 has 0.63 km s -1 and 0.33 km s -1 velocity resolution data, respectively. CN J = 2-1 225.632 GHz data with a synthesized beam of 0 '' . 57 × 0 '' . 38 and velocity resolution of 0.08 km s -1 are from 2017.1.01413.S (PI: van 't Ho ff ).  \nDue to the higher sensitivity of the latter dataset, we used it to construct a continuum image at 1.3 mm from line-free spectral regions for further analysis. We used pipeline-calibrated data images obtained from the ALMA archive for both projects.", '3. Results': 'The MIRI-MRS cubes of the TMC1 protostellar system reveal two continuum point sources. Across the MIR spectrum, we also detect several spatially extended features associated with emission lines. An overview of the system is presented in Fig. 1.  \nFig. 1. Continuum emission at 5.4 µ m obtained from line-free spectral channels in Channel 1 (left). Center: Integrated emission of H2 0-0 S(7) line at 5.51 µ m in colorscale. Right: Integrated emission of [Fe II] line at 5.34 µ m in color scale. Stars mark the MIR peak continuum position of the protostars. White circles show regions from which the spectra at the outflow positions were obtained, with diameters corresponding to two times the MIRI-MRS point-spread function (PSF) at 5.5 µ m. In the bottom-left corners, MIRI-MRS empirical FWHM of PSF (Law et al. 2023) is indicated as a white circle.  \n<!-- image -->  \n20  \nFig. 2. Spectra extracted from the MIR continuum peaks of the two protostars. Several absorption features associated with ice and dust are seen in both sources and emission lines coming from the disk and central protostars, although some contamination with outflows and winds is also possible. Top-left: Zoom-in on the two hydrogen recombination lines detected in emission towards protostars. Bottom-right: Spectra revealing molecular gas emission and corresponding slab models of CO (green) and H2O (orange).  \n<!-- image -->  \n<!-- image -->  \nFig. 3. MIRI-MRS spectra from the outflow positions. Spectra were continuum subtracted using a spline fit to line-free regions. The H2 transitions are shown in purple, the hydrogen recombination line in yellow, and fine structure atomic lines in green. Top: TMC1-W jet (red) and TMC1-E wind (blue). Bottom: TMC1-E wind spectra. OH catalog emission line positions are shown as black markers (Yousefi et al. 2018; Tabone et al. 2021). Insets: [Co II] 10.521 µ m, HCO + 12.069 µ m, and CO2 14.98 µ mlines identified in the spectra. Flux in the insets converted to mJy.  \n<!-- image -->', '3.1. MIRI-MRS Spectra': "This section analyzes spectral information extracted from positions indicated in Fig. 1. The extracted regions were selected to cover both continuum peaks associated with two protostars of the TMC1 system and the representative outflow positions. These positions allow for optimal disentanglement of the outflow components. The extraction aperture at the point source was set to ra = 0 . 033 × λ + 0 '' . 106, which at 5 µ mis 0 '' . 271 and at 20 µ mis 0 '' . 766. This is twice the size of the FWHM of MIRI-MRS PSF as characterized by Law et al. (2023). At the outflow positions with prevalent extended emission, we used a circular aperture of 0 '' . 25 radius to ensure that equal regions are extracted at di ff erent wavelengths.  \nThe location of extracted spectra for continuum peaks and outflow positions are indicated by star symbol and white circles in Fig.1, respectively. After spectral extraction, additional  \nfringe correction is applied to the spectra (Gasman et al. 2023, Kavanagh et al. in prep.). A Gaussian profile was fitted to each identified emission line in the spectra and the flux value of the line was estimated by integrating the region under the Gaussian. The line intensity was then obtained by dividing the final flux value with the extraction aperture area.  \nTo obtain the true line intensity, an extinction correction should be applied, accounting for attenuation of the flux by the envelope in which the forming star is embedded. For TMC1 protostar, a visual extinction of AV = 49.2 mag is reported in Connelley & Greene (2010) based on spectral energy distribution shape in the 0.8 to 2.4 µ m range. Fiorellino et al. (2023) found a value of 36.77-35.70 mag for birth line and at 1 Myr, respectively. However, both those methods assume that the flux comes from a single star; therefore, these estimates are unreliable. In Apai et al. (2005) near-IR colors of two components are evaluated separately, and TMC1-W and TMC1-E are found to have  \nAV of 43 and 37 mag, respectively. However, the colors of protostellar sources can be significantly a ff ected by the presence of the disks that both protostars have, as evident in submillimeter continuum images (Harsono et al. 2014; van 't Ho ff et al. 2018), and the models of color correction developed for T Tauri stars can underestimate the extinction when used in the case of younger protostars (Doppmann et al. 2005).  \nGiven the absence of reliable extinction estimation for this source, we use MIRI-MRS data to estimate AV with di ff erent methods for on-source and o ff -source positions. On protostars, we measure silicate feature depth after the continuum subtraction using a third-order polynomial function. This process is presented in the Appendix A.1. We find τ 9 . 7 = 1 . 30 and τ 9 . 7 = 1 . 95 for TMC1-W and TMC1-E, respectively.  \nConversion of τ 9 . 7 to AV bears large uncertainties, with factors of A V = 18 . 5 × τ 9 . 7 used for ISM and A V = 9 × τ 9 . 7 used for dense clouds (McClure 2009; Weingartner & Draine 2001). This is due to the total-to-selective extinction ( RV ) in the two environments (Chiar et al. 2007). McClure (2009) shows that extinction values toward star-forming regions are in between those provided for the galactic center and di ff use ISM. Therefore, we use those two limits as a range for the extinction estimation. This results in AV = 11 . 7-24 . 0 mag and AV = 17 . 6-36 . 1 mag for TMC1Wand TMC1-E, respectively.  \nAt the o ff -source position, where no continuum is present to measure the depth of the silicate feature, we use H2 emission lines. We fit simultaneously for the rotational temperature and extinction value (Barsony et al. 2010; Federman et al. 2023). Details are provided in the Appendix A.1. We only use v = 0-0 transitions because other non-thermal e ff ects can contribute to the ratio of fluxes between v = 0-0 and v = 1-1 transitions (Timmermann et al. 1996). The value measured at TMC1-E blueshifted wind position AV = 20 . 8 ± 3 . 1 mag and at TMC1-W blueshifted jet AV = 23 . 1 ± 2 . 8 mag. We use these values to correct for line fluxes measured at both positions.  \nWe note that considering the ranges of the extinction values found for the on-source and outflow positions, there is a possibility that the extinction values in the outflow would be higher than on-source. This is unexpected as much larger extinctions are typically measured close to the protostar given strongly decreasing gas density with increasing envelope radius (Jørgensen et al. 2002b). However, it should be considered that dust grain properties change across the envelope, with rapid dust size increase close to the protostar (Cacciapuoti et al. 2023), making the two extinction values challenging to compare as the total-to-selective extinction depends strongly on dust size (Chiar et al. 2007).", '3.1.1. Spectra at source positions': "Figure 2 presents spectra extracted from the central position of the protostars on the MIRI-MRS detector. From Channel 3 short onward ( > 11.55 µ m), the regions of spectral extractions significantly overlap; however, at shorter wavelengths, it is possible to disentangle the emission from the two sources and discuss their di ff erences. Both sources have spectra characteristic of embedded protostars, with a spectral energy distribution (SED) rising with wavelength. The spectra show strong emission lines of molecular H2, atomic and ionized components, along with deep silicate and ice absorption features, such as H2O (6 µ m), CH3OH (7 µ m), silicates (9 µ m), and CO2 (15 µ m) (Ehrenfreund et al. 1997; Keane et al. 2001; Zasowski et al. 2009; Bottinelli et al. 2010; Yang et al. 2022; Rocha et al. 2022).  \nThe spectra also reveal weaker emission lines towards both sources, such as emission from CO molecular gas in the 5 µ m  \nspectral region and H2O emission lines detected in 5-8 µ mrange (see bottom-right inset in Fig 2). In the same inset, we compare spectra towards both targets with CO and H2O slab gas spectra at 1000 K under the assumption of LTE. In the 9-20 µ mspectral range, OH emission is prominently detected towards TMC1-E, with only a hint towards TMC1-W (which could also be the result of emission from TMC1-E included in the aperture). CO2 gas-phase emission lines at 14.98 µ m are detected towards both sources on top of the deep ice absorption feature.  \nMolecular H2, atomic, and ionized emission lines are detected towards both targets. However, verifying whether the emission comes from the inner disk, protostar, or outflow is difficult. The continuum-subtracted maps of the H2 and [Fe II] line in Fig. 1 clearly indicate that these lines originate mainly o ff -source. Therefore, while the contribution of the central region cannot be discarded, we assume that most of H2 and ionized emissions are coming from the winds or jets of both protostars.  \nThe exceptions are atomic hydrogen recombination lines, which are brightest at the on-source apertures; see top-left inset in Fig 2. In Table B.1, we report extinction-corrected line intensities of the atomic hydrogen transitions detected towards the two binary components of the TMC1 system. We identify six transitions of atomic hydrogen toward both protostars, with line intensities significantly higher towards TMC1-W. These lines are expected to trace high temperature and density environments and are usually associated with stellar accretion (Rigliaco et al. 2015; Alcalá et al. 2017). However, not all MIR lines have derived relations between the accretion rate and line fluxes. The best one to date is the HI 7-6 line for which the relation is provided in Rigliaco et al. (2015).  \nIn Section 3.1, we provided ranges for the extinction values for TMC1-W and TMC1-E based on the di ff erent correlations found for di ff erent grain properties. It is safe to assume that the dust grains should have similar properties toward both protostars. Therefore, the relative error on extinction should not be as large as the range of values provided. In the subsequent analysis, we assume the mean value from both ranges. This is AV = 27 mag and AV = 18 mag for TMC1-E and TMC1-W, respectively.  \nWe obtain the extinction corrected line luminosity for this line to be 3 . 8 × 10 -7 L ⊙ and 1 . 8 × 10 -7 L ⊙ , for TMC1-W and TMC1-E, respectively. Applying the relation between HI 7-6 lines and accretion luminosity derived from H α provided in Eq. 1 in Rigliaco et al. (2015), and accounting for the uncertainty on this relation, this corresponds to the accretion luminosity of 2.3 + 9 . 1 -1 . 8 × 10 -4 and 0 . 5 + 3 . 5 -0 . 48 × 10 -4 L ⊙ , respectively. Next, using Equation 11.5 from Stahler & Palla (2005) (see also Gullbring et al. 1998), we get:  \n˙ M = L acc GM ∗ R -1 ∗ , (1)  \nand taking the range of stellar properties of TMC1, namely, M ⋆ = 0.15-0.20 M ⊙ and R ⋆ = 3.0-2.48 R ⊙ for birthline and 1 Myr, respectively, from self-consistent modeling in Fiorellino et al. (2023), we find the accretion rate of TMC1-E to be (0 . 02 -0 . 43) × 10 -11 M ⊙ yr -1 and of TMC1-W to be (0 . 16 -14 . 8) × 10 -11 M ⊙ yr -1 . This range includes uncertainties on the accretion luminosity.  \nThere are several caveats to this analysis. First, the atomic hydrogen emission lines do not have to trace exclusively accretion but can also be detected in the protostellar jet and wind, which is evident especially in recent spatially resolved JWST observations of protostars (Federman et al. 2023; Harsono et al. 2023). While we do not detect significant hydrogen emission  \nFig. 4. H2 rotational diagram obtained with line intensities measured at the TMC1-E wind position. All detected H2 transitions are presented, but only v = 0-0 (red) transitions are used to fit the temperature.  \n<!-- image -->  \nlines in the jet, they could still contribute to the emission close to the source.  \nSecond, stellar properties in Fiorellino et al. (2023) are obtained under the assumption of a single star dominating the flux and, therefore, more accurate for the brighter TMC1-W, with TMC1-E properties remaining uncertain. Additionally, the flux obtained for the HI (7-6) line towards TMC1-E might be contaminated by contribution from TMC1-W, given the large PSF at 12.37 µ m(FWHM ∼ 0 '' . 5). Given that extraction region size at this wavelength is 0 '' . 5, and the distance between the sources is 0 '' . 65, there can be significant contribution from TMC1-W. Comparing ratios of HI lines between the sources, we find that the HI (9-6) line, at 7.5 µ m, is 17 times brighter for TMC1-W compared with TMC1-E, while the HI (7-6), at 12.37 µ m is only 2 times brighter. This further suggests increasing the contribution of TMC1-W flux to TMC1-E aperture with wavelength. The considerable uncertainty on the extinction correction also adds to the di ffi culty when comparing the two sources. In summary, while highly uncertain, these measurements indicate the stellar accretion conditions for both targets, tentatively pointing to TMC1-W being a stronger accretor in the system.", '3.1.2. Spectra at the outflow positions': "In Fig. 3 (top panel), we show 5-7.6 µ m spectra extracted from outflow positions for TMC1-E and TMC1-W outflows. At longer wavelengths, the spectra are a ff ected by contribution from the central source PSF. Hence, we selected positions with a significant distance from the protostar for spectral extraction but close enough to be within the MRS footprint in all channels. The spectra plotted in Fig. 3 were extracted with an aperture size of 4 × PSF FWHM to minimize the e ff ect of undersampling of the spectral information in the MRS detector, while for the flux measurements, we used a fixed circular aperture size of 0 '' . 25 radius.  \nH2 pure rotational transitions are prominently detected in our data, primarily towards the TMC1-E wind position. The detection of H2 from TMC1-W is likely the result of contamination with the western side of the TMC1-E wind (see Fig. 1). The line intensities of H2 lines towards TMC1-E outflow are reported in Table B.2.  \nTo analyze the physical conditions associated with H2 emission, we construct a rotational diagram using pdrtpy package  \n(Pound & Wolfire 2008). Figure 4 shows a rotational diagram produced using the measured fluxes corrected for the visual extinction (see Section 3.1). We only used v = 0-0 transitions to quantify the rotational temperature. By fitting simultaneously two linear functions to the data, the negative inverse of the slope provides the rotational temperature of two separate components: a warm component at 661 ± 28 K and a hot component at 1849 ± 167 K. Taking the intercept of the linear functions multiplied by the partition function (Herbst et al. 1996), we find the total column density of the two components: 2.5 ± 0.3 × 10 19 cm -2 for the warm and 1.0 ± 0.4 × 10 18 cm -2 for the hot component, with hot gas contributing only 1% to the total column density. We do not provide a similar analysis for the TMC1-W jet region since any detected H2 emission at this position is likely coming from the TMC1-E outflow.  \nWhile the molecular CO and H2O have strong emissions on the source positions, they are not detected in the outflow. Other molecules, namely, HCO + , OH, and CO2, appear to be emitted from the outflow (see bottom panel of Fig. 3); however, the contamination from the source position cannot be excluded. The lines are not bright enough to map them with a su ffi ciently high signal-to-noise ratio (S / N). We detected several atomic and ionized emission lines at the outflow positions, namely, [Fe II], [Ne II], [Ne III], [Ni II], [Ni I], [Ar III], [Ar II], [S I], [Cl II], [Cl I]. We also tentatively detected [Co II] (Fig. 3). We found only one instance of detection of this species in the protostellar jet in Giannini et al. (2015). Due to the lack of other detected transitions, this identification is tentative. However, this is a relatively lowenergy transition to the ground level with E up of 1368 K and a high Aij coe ffi cient of 2 . 2 × 10 -2 s -1 (Nussbaumer & Storey 1988), so this is a plausible association that requires more investigation. Table B.3 reports the extinction-corrected line intensities for both outflows. Table B.4 reports the relevant line ratios for both outflow positions.", '3.2. Emission maps: MIR characterization of the TMC1 system': 'This section presents the emission maps of lines detected toward the TMC1 binary system. Figure 5 shows H2 0-0 S(1), S(3), and S(7) lines, while Fig. 7 shows [Fe II], [Ni II], [Ar II], [Ne II], and [Ne III]. To produce the continuum-free emission line maps, we subtracted the line-free part of the spectral cube close to the line from the integrated emission map centered on the line wavelength. We only consider lines with a sufficiently high S / N to separate the extended emission from the continuum and the noise. Therefore, not all the detected lines are present on the maps. The spectral resolution of the MIRI-MRS (R ∼ 1500 -3500 Labiano et al. 2021) allows for the production of channel maps of the emission lines for the high-velocity gas in protostellar jets. The maps are presented in A.4.', '3.2.1. Molecular hydrogen emission lines': "Molecular hydrogen (H2) is an excellent tracer of the bulk of the outflowing gas from protostars; however, due to its lack of permanent electric dipole moment, it can only be detected through ∆ J = 2 rotational transitions with relatively high excitation energies ( E up > 500 K), thereby tracing the heated gas (Neufeld et al. 2009). We see that H2 emission originates from TMC1-E, with the apex of the parabolic shape of the emission pointing to this source. Figure 5 shows the rotational S-branch transitions of v = 0 -0 state of H2. While we cannot completely exclude  \nFig. 5. Three pure rotational transitions of H2 from left to right: S(7) at 5.51 µ m, S(3) at 9.66 µ m, and S(1) at 17.04 µ m. Integrated emission maps are created from the cubes after subtracting the local continuum. The positions of the protostars measured from the line-free region at 5 µ m are shown with white stars. Maps are centered on the position of the TMC1-W source. The horizontal dotted line in the right panel indicates the cut for which di ff erent H2 0-0 S bands are shown in Fig. 6. The physical components of the system are illustrated in the middle panel. The disk wind is indicated based on H2 S(7), CO outflow is based on the ALMA observations, and scattered light on the HST 1.6 µ mimage in Fig. 8. Bottom-left corners: MIRI-MRS empirical FWHM of PSF (Law et al. 2023) indicated as a white circle.  \n<!-- image -->  \nthe H2 coming also from the TMC1-W (see fainter broad component in H2 S(7) line in Fig. 5), we do not see any further extent of this wind at longer wavelengths, especially at S(1), where the wide-angle wind is completely dominated by the TMC1-E wind. Therefore, a significant contribution to H2 emission from TMC1-W is unlikely.  \nH2 emission is present inside the cavity and shows progressively larger opening angles with decreasing upper energy levels. Additionally, two components can be observed: a narrow V-shaped emission and a much broader parabolic shape (annotated on H2 S(3) map). We interpret those as tracing the less collimated wind and the shocks filling the whole outflow cavity wall, respectively. It can be seen that the emission from the wind component remains narrow and compact across all H2 transitions, while the outflow cavity component becomes more prominent with decreasing E up. Maps reveal outflow emission on the south side of the source as well. However, this emission is much fainter, likely because this is the redshifted side of the flow, hidden on the other side of the protostellar envelope of the TMC1 system.  \nTo illustrate the increasing opening angle of the H2 emission, in Fig. 6, we show radial cuts across the blueshifted outflow of TMC1-E. For the high E up lines, there is a clear two-peaked feature that can be associated with the disk wind, and the whole extent of the wind does not exceed 1 '' , while at lower E up broader component is also detectable extending up to 2 '' . 5.", '3.2.2. Ionized and atomic species': 'Atomic and ionized species in MIR can trace various physical components of the systems. Spatially unresolved Spitzer -IRS observations identify the origin of lines in jets, disks, and photodissociation-like environments in the outflow cavity walls (Hollenbach & Gorti 2009; Lahuis et al. 2010). The MIRI spa-  \ntial resolution is invaluable in shedding new light on the origin of those lines.  \nIn Fig. 7, we present key atomic and ionized tracers detected toward the TMC1 system. Based on their appearance, they can be divided into two categories: refractory and noble gas ions. The refractory ions [Fe II] and [Ni II] show a collimated jet (radially unresolved) appearance. They are much brighter in the jet from TMC1-W where H2 lines are undetected, while only weakly detected towards source TMC1-E, where the H2 emission is prominent. On the other hand, noble gas ions [Ar II], [Ne II], and [Ne III] are present in both sources, but they are distinctly brighter in the TMC1-E outflow. They also appear moderately extended in the radial spatial direction, contrasting the collimated refractory ions. This is particularly clear in the channel maps where [Ni II] and [Ar II] are compared (Fig. A.4).  \nChannel maps reveal velocity information contained in the data obtained with MIRI. In the case of the TMC1-W jet both [Ni II] and [Ar II] are detected at velocities up to 200 km s -1 , while gas at lower velocities is detected in the case of TMC1-E, especially the radially extended component of the TMC1-E wind is not detected at velocities exceeding 100 km s -1 .', '3.3. Molecular emission at submillimeter wavelengths and scatter light near-IR emission': "For young stellar objects (YSOs), putting the MIR observations in the context of colder gas tracers detected in the (sub)millimeter range is particularly insightful. This is because those objects are still embedded in their natal envelope and many are still accreting material from the cloud onto the disk. Additionally, colder molecular gas will be released in entrained outflow. In particular, the lowJ rotational transitions of CO are excellent tracers of the low-velocity outflows (e.g., Arce et al. 2010; Plunkett et al. 2013; van der Marel et al. 2013). The (sub-)millimeter range is also less a ff ected by extinction. The velocity resolution  \nFig. 6. Radial cuts across the TMC1-E outflow at 1 '' . 5 distance from the protostar for each of the H2 (0-0) transitions, with upper energy levels increasing from top to bottom. The reference is set to the protostar position. The position of the cut is shown in Fig. 5. The black bar shows the size of the FWHM of instrumental PSF, as described in Argyriou et al. (2023).  \n<!-- image -->  \nof ALMA below 1 km s -1 could provide valuable insights into interpreting velocity unresolved information from MIRI-MRS.  \nNear-IR or optical images of embedded systems at shorter wavelengths show scattered light coming from the protostellar envelopes or outflow cavity walls. These observations are helpful as they inform about the geometry of the system and the presence of small grains.  \nWe show an archival Hubble Space Telescope NICMOS 1.6 µ m image of the TMC1 system in Fig. 8 (left panel), presented in Apai et al. (2005) and Terebey et al. (2006). The extended emission on the northern side of the protostar is attributed to the scattered light at the outflow cavity walls. The scattered light presents a parabolic shape with an apex at the TMC1-W protostar. A collimated jet is also emerging from the TMC1-W source, with [Fe II] line at 1.64 µ m, likely to be dominating the jet emission in the F160W imaging filter (Terebey et al. 2006). Scattered light most likely shows the boundaries of the outflow cavity walls. Since they appear broader than the CO and H2 outflow currently observed towards TMC1-E, this suggests that in the past, both TMC1-W and TMC1-E contributed to shaping the outflow cavity walls; however, as currently observed, there is no visible wide-angle wind from TMC1-W.  \nWe also note that no scattered light is detected with MIRI, in contrast with Class 0 sources also observed at MIR (Federman et al. 2023). This highlights the more evolved nature of the TMC1 system compared with Class 0 protostars.  \nIn Fig. 8, we present ALMA maps of CO 2-1 (central panel), SO 56-45 (right panel), and CN 2-1 (contours in the central and right panel) from Tychoniec et al. (2021). CO emission traces the  \nentrained outflow material piling up at the outflow cavity walls of TMC1-E. Notably, CO emission has a narrower opening angle compared with the scattered light image, and the center of the symmetry of CO outflow appears to be aligned with the TMC1E. Terebey et al. (2006) found the best fit for the system's inclination of 50 · based on the scattered light modeling and already identified the opening angle of the scattered light as being larger than CO outflow.  \nThe CN emission is narrow and not co-spatial with the outflow cavity walls, as evident from the comparison with the CO map. The velocity gradient is consistent with infall motion and inconsistent with outflow, and velocity dispersion is small < 0.2 km s -1 (Fig. A.3). Also, SO highlights the gas likely associated with weak shocks, which can inform where the infalling material interacts with the envelope and disk (Sakai et al. 2014). It may also be tracing the young, embedded disk not associated with accretion, which can explain its presence towards both TMC1W and TMC1-E (Harsono et al. 2014). It is worth noting that even though the SiO 5-4 transition is covered in the ALMA data, it was not detected towards TMC1. This further highlights the more evolved nature of the protostars. High-velocity SiO jets are commonly observed towards the youngest, Class 0 protostars (Hirano et al. 2010; Tychoniec et al. 2019; Podio et al. 2021), where molecules can form e ffi ciently in the jets due to high densities (Tabone et al. 2020). As the outflow evolves, the mass loss rate decreases, velocities increase, and the jet becomes progressively more ionized (Nisini et al. 2015).  \nAll observed features are consistent with an accretion streamer feeding the system (Le Gouellec et al. 2019; Pineda et al. 2020; Valdivia-Mena et al. 2022). We show CN data as an example of the accretion streamer because (for those lines) the highest spectral resolution data are available. However, closer to the source, there is significant contamination from the envelope and disk material. Moreover, SO emission aids in pinpointing currently shocked gas, from which we can see that the streamer appears infalling mostly onto the disk around the TMC1-E protostar.", '4. Discussion': 'MIRI-MRS observations have revealed for the first time a clear di ff erence between the outflows launched by the two binary protostars. In this section, we investigate the origin of the di ff erent outflows. The key features we discuss are as follows.  \nH2 emission originates from the TMC1-E source, while no significant H2 emission from TMC1-W is detected. The H2 emission from TMC1-E has two components: a wind with a narrow opening angle and broader outflow with a wider opening angle, contained within the molecular CO outflow and scattered light emission. Prominent emission from [Fe II] and [Ni II] toward TMC1-W is detected, while both sources show emission from [Ar II] and [Ne II], with the TMC1-E outflow showing bright and radially resolved emission in those tracers.  \nOur proposed scenario that explains the observed features is that di ff erent modes of accretion influence the origin of two outflows: the wide-angle H2 wind towards TMC1-E is influenced by the accretion from the large-scale envelope onto the disk, which provides additional material, while the disk-to-star accretion is related to the launch of the powerful, collimated jet from TMC1W. Figure 9 shows an illustrative schematic diagram of the observed features. We present below a detailed discussion of the results and the proposed scenarios.  \nFig. 7. Integrated emission maps of selected emission lines of ionized species. Integrated emission maps are created from the cubes after subtracting the local continuum. Positions of the protostars measured from the line-free region at 5 µ m are shown with white stars. Maps are centered on the position of the TMC1-W source. Bottom-left corners: MIRI-MRS empirical FWHM of PSF (Law et al. 2023) indicated as a white circle.  \n<!-- image -->', '4.1. TMC1-E: Disk wind possibly powered by infall from the envelope to the disk': "The H2 emission originates in the warm and hot molecular gas in the outflows; however, without spatially resolving the emission, it is di ffi cult to pinpoint the specific location of the hot gas. With MIRI-MRS, we see that the highest E up transitions of H2 are more collimated than the low E up ones. We interpret this as H2 tracing the disk wind, with a flow radius that is larger than the jet, but clearly distinct from the broader entrained envelope material at the outflow cavity walls 5. This is further corroborated by a comparison with the ALMA CO 2-1 map and HST NICMOS scattered light image (Fig. 8). Both delineate the extent of the outflow cavity walls, showing that high E up H2 transitions are within this broader component traced by CO, scattered light and lower energy H2 lines. Velocity stratification of the wind and interplay between the jet and the surrounding gas can be responsible for the observed morphology of the gas. In theoretical works, both jet-driven bow shocks (e.g., Rabenanahary et al. 2022), and wind-driven cavities can explain those structures (e.g., Shang et al. 2023).  \nWith ALMA, we also detected a prominent accretion streamer infalling onto the binary. Tentatively, we associate the TMC1-E with this streamer's infall point, as SO and CN maps suggest. That infall could trigger the prominent disk wind from this source by providing fresh material to replenish the disk. The ionized tracers [Ar II], [Ne II], and [Ne III] are all bright toward this source. Since those noble gas ions have an ionization potential above 13.6 eV, they originate most likely from UV-irradiated material in the disk wind. We base this on the shape of the emission, clearly extended around the base of the wind and not collimated, as seen in the case of [Fe II] or [Ni II] lines. While the disk wind and related UV-excitation is one possibility, some models predict decollimation of the jet in binary systems, leading to a similarly broad appearance of the jet (Lynch et al. 2020).  \nThe line ratios observed toward the TMC1-E wind position suggest UV-irradiated environments. For example, [Ne III] / [Ne II] = 0.04 is consistent with soft EUV irradiation dominating the excitation and ionization of Ne (Hollenbach & Gorti 2009; Szulágyi et al. 2012; Espaillat et al. 2013). We detected OH lines towards TMC1-E (Fig. 3). Those lines are an indicator of UV photodissociation of H2O (e.g., Tabone et al. 2021). It is interesting to note that H2O appears more prominent in the gas phase  \nFig. 8. HST NICMOS 1.6 µ mimage of the TMC1 system (far-left). ALMA CO (2-1) blueshifted and redshifted integrated emission in colorscale, ALMA 1.3 mm continuum in white contours, and CN peak emission map (moment 8) in cyan contours (center-left). ALMA SO (6-5) integrated emission map in colorscale (center-right). NICMOS 1.6 µ m image in colorscale with ALMA CO(2-1) overlaid in red and blue contours, and H2 S(7) moment 0 map in black contours (far-right).  \n<!-- image -->  \nin TMC1-W compared to TMC1-E (Fig. 2); this could be a result of significant water destruction in the TMC1-E source. Ionization and the presence of charged particles can also indicate strong magnetic fields towards this source, which could translate to more prominent disk wind.  \nWhile UV-irradiation is a plausible explanation for the observed ionized species close to the source, further out in the jet, strong, high-velocity shocks ( > 100 km s -1 ) can also explain this emission (Hollenbach & McKee 1989). This is especially valid in the case of the TMC1-W jet, where velocities are approaching 200 km s -1 , stronger emission coming from, for example, [Ar II] is observed further out in the jet in high-velocity channels (see Fig. A.4).  \nInterestingly, in the best-studied Class I disk wind source, DG Tau B, the ejection rate in the wind is significantly higher than that of the collimated jet and several times higher than the mass accretion rate on the protostar (Podio et al. 2011; de Valon et al. 2020). Our results show that a potential source of this high mass-loss rate in the wind could be the influx of the envelope material. Streamers of gas into the protostar system could therefore have an essential impact on the evolution of the protostellar system. Both observations and simulations suggest they could bring a significant fraction of gas and dust to the protostar and disk - even at later stages of evolution (Valdivia-Mena et al. 2022; Garufi et al. 2022; Ku ff meier et al. 2023).  \nFrom a theoretical point of view, simulations of black hole binaries suggest that if the binaries are formed by disk fragmentation, the accretion from the envelope occurs preferentially onto the secondary component, eventually leading to the two components growing similar in mass (Farris et al. 2014). Tokovinin & Moe (2020) used those simulation results to successfully reproduce protostellar binaries, explaining the scarcity of brown dwarfs in binary systems. This would be consistent with our ALMA observations indicating the infall of the streamer onto the secondary component. This has been observed in another protostellar system BHB2007 (Alves et al. 2019). On the other hand, in Murillo et al. (2022), the large-scale infall occurs on the more massive component of the IRAS16293 protobinary system. However, in the case of that study, a separation between protostars is 720 au, therefore, it is not likely to have been formed as a result of disk fragmentation (O ff ner et al. 2010).  \nWhat is particularly striking in this system is that the protostellar jets are launched almost perfectly in parallel. From this, we can assume that the disks from which those jets are launched have co-planar geometry, supporting the protostars' common  \norigin. Interestingly, no interaction or disturbance of the jet from TMC1-W by the TMC1-E wind is observed, in contrast to expectations of the models, where jets from binary systems are expected to cause precession or even to become completely attenuated or engulfed in joint flow (Murphy et al. 2008; Lynch et al. 2020). Most observations of young binary systems also show precession and deflection of jets (Kwon et al. 2015; Ching et al. 2016; Fridlund & Liseau 1998).Therefore, the TMC1 system presents the most striking example of parallel jets to date. Standing at odds with other dynamically interacting protostars, this indicates that the disks of the two protostars are co-planar.", '4.2. TMC1-W: Jet possibly powered by the accretion onto the protostar': 'We detected six transitions of hydrogen recombination lines toward the TMC1 system. Those lines are expected to trace high temperature and density environments and are usually associated with stellar accretion (Rigliaco et al. 2015; Alcalá et al. 2017). Detections of brighter HI lines toward TMC1-W compared with TMC1-E, even after accounting for extinction can indicate that the TMC1-W source is undergoing more vigorous accretion. A natural consequence of this should be a release of the excess angular momentum close to the accretion streamlines in the inner disk. According to the MHD jet launching paradigm (Matt & Pudritz 2005; Romanova et al. 2005) smaller launching radii should correspond to faster ejection and smaller opening angles, leading to the launch of the highly collimated jet, as observed toward TMC1-W.  \nThe accretion rate measured from the HI (7-6) line is (0 . 16 -14 . 8) × 10 -11 M ⊙ yr -1 . This is a low value compared to the other Class I systems, more typical for Class II systems. (Fiorellino et al. 2022). For the same HI line in a more massive and younger system, values can be orders of magnitude higher (Beuther et al. 2023).  \nIt is possible to have a rough estimate of the ejection rate from the luminosity of the [Fe II] 25.96 µ mline and using Equation 2 in Watson et al. (2016). This latter allows us to derive a result of [Fe II] 25.96 µ m ˙ Mwind relation empirically scaling the relation between [O I] 63.2 µ m and Mwind provided by Hollenbach & McKee (1989). Using the extinction-corrected flux measured towards the total extent of the blueshifted jet of TMC1W, we obtained a line luminosity of 4.8 × 10 -6 L ⊙ . We derived a mass-loss rate of 6.7 × 10 -9 M ⊙ yr -1 .  \nCompared with the accretion rate measured on the source, the ejection rate is 45-4000. It is therefore di ffi cult to reconcile with the models of magnetohydrodynamically-driven winds, where the value ranges from 0.01 to more than 0.5 depending on the details of the model (Nisini et al. 2018). It is possible that the accretion rate is highly variable and was much higher in the past or that the flux from the [Fe II] line is coming from jets of both protostars.  \nOther methods can be also used to extract mass-loss rate from the observed lines, such as measure of the total number of emitting molecules (Dionatos et al. 2009). However, the depletion of iron onto the dust grains would contribute to the uncertainty of the measurement in the context of this work.  \nIt is noteworthy that primarily refractory species such as iron, nickel, and chlorine are present in the jet. All those species are expected to be incorporated into the dust grains in the interstellar medium. Therefore, their presence in the gas phase is directly linked to the destruction of the grains. Shock models reveal that only up to 10% of the refractory material is typically released into the gas phase due to grain destruction in shocks (Gusdorf et al. 2008a,b). Observational work to date also confirms that refractory abundances in the jet are depleted with respect to the solar values, indicating that the bulk refractory material is still in the solid state (Podio et al. 2006, 2011; Dionatos et al. 2009; Giannini et al. 2015).  \nWe interpret the presence of refractories as a sign of jet launching within the dust sublimation radius in the protoplanetary disk, where the iron and nickel are released into the gas phase and further launched with the jet (e.g., Podio et al. 2011) This further supports the idea that the TMC1-W source is dominated by the highly collimated jet launched from the inner regions of the disk triggered by strong accretion onto the protostar. The di ff erence between TMC1-E and TMC1-W can also stem from the two jets having di ff erent dust content, namely, the more dusty TMC1-E jet has a lower fraction of the refractory material released in the gas phase.', '4.3. Alternative scenarios': 'We provide a scenario where the accretion activity in the system is mainly responsible for the striking di ff erence in the outflow behavior of the two protostars. Here, we consider alternative scenarios. It is possible that the two protostars have followed di ff erent evolution histories and that, at present, TMC1-W has ended up to be more evolved than TMC1-E. Its lower extinction and stronger ionized jet, without any significant molecular wind component, would support this.  \nAlternatively, the outflow behavior could be governed by the properties of protoplanetary disks surrounding both protostars. ALMA images show compact dust continuum emission from both targets. If the TMC1-W disk is truncated due to e ff ects such as magnetic braking, this could also make the production of the disk wind di ffi cult. Magnetization of the disk could also play an essential role in the launching of the wind. Therefore, di ff erent magnetic field properties in the components could contribute to these di ff erences. It is known that binarity plays an important role in the accretion process, disk truncation, and outbursting activity (Zagaria et al. 2021) Lastly, even if the system followed a similar evolution process, it is possible that inflow of the material would cause apparent rejuvenation of the protostar by increasing the amount of material available for accretion (Ku ff meier et al. 2023).', '5. Conclusions': "JWST MIRI-MRS observations of Class I protobinary system TMC1have revealed an unprecedented view of the accretion and ejection processes in protostars. The results from our study are summarized as follows:  \n- -Wedetect two distinct outflows from both protostars, TMC1W and TMC1-E. The outflow from TMC1-W shows collimated emission lines of [Fe II] and [Ni II] tracing the energetic high-velocity jet. At the same time, TMC1-E powers a prominent wide-angle traced with H2 and is also associated with more collimated emission from [Ne II], [Ne III], and [Ar II]. Jets show a rare parallel geometry that has not been seen before in binary systems, indicating that disks are co-planar.\n- -Six hydrogen recombination lines are detected toward both targets. The lines are much brighter toward TMC1-W than towards TMC1-E, which translates to a higher accretion luminosity and, thus, to higher accretion rates of the TMC1-W source, even when accounting for a higher extinction towards the TMC1-E source. The accretion activity of TMC1-W is likely related to the stronger ejection of material via a collimated jet from this source.\n- -ALMA CO 2-1 map and scattered light emission from HST at 1.6 µ m show that hot H2 traced with MIRI come from within the outflow cavity walls, most likely tracing the disk wind from the TMC1-E source. At the same time, with ALMA, we detect an infalling streamer of gas with CN and SO, pointing to TMC1-E as the landing point of the infall. This could influence the disk growth and transport of the angular momentum from the envelope to the disk and is possibly related to the powerful disk wind launched from this source.  \nIn this work, we observe a remarkable di ff erence in the appearance of the outflow from this binary system in di ff erent atomic and molecular tracers. This comparison allows us to achieve a better understanding of what physical components are probed by di ff erent MIR tracers. In this context, JWST-MIRI has opened up a new avenue of study. Additionally, by linking the outflow properties with stellar accretion measured with our MIR data and infall tracers, as seen by ALMA, we o ff er a scenario that allows us to explain the di ff erences between outflows that originate in di ff erent modes of accretion and infall governing the ejection towards the protostars. This provides a new insight into the connection between accretion and ejection that the models of protostellar collapse and star and disk formation should further investigate.  \nAcknowledgements. Authors would like to thank the referee for helpful suggestions that improved the clarity of the paper. Ł.T. would like to thank Rajika Kuruwita, Aashish Gupta, and Ugo Lebreuilly for helpful discussions and Marta Tychoniec for contributing to Fig. 9. This work is based on observations made with the NASA / ESA / CSA James Webb Space Telescope. The data were obtained from the Mikulski Archive for Space Telescopes at the Space Telescope Science Institute, which is operated by the Association of Universities for Research in Astronomy, Inc., under NASA contract NAS 5-03127 for JWST. These observations are associated with program #1290. The following National and International Funding Agencies funded and supported the MIRI development: NASA; ESA; Belgian Science Policy O ffi ce (BELSPO); Centre Nationale d'Etudes Spatiales (CNES); Danish National Space Centre; Deutsches Zentrum fur Luft- und Raumfahrt (DLR); Enterprise Ireland; Ministerio De Economiá y Competividad; Netherlands Research School for Astronomy (NOVA); Netherlands Organisation for Scientific Research (NWO); Science and Technology Facilities Council; Swiss Space O ffi ce; Swedish National Space Agency; and UK Space Agency. H.B. acknowledges support from the Deutsche Forschungsgemeinschaft in the Collaborative Research Center (SFB 881) 'The Milky Way System' (subproject B1). EvD, MvG, LF, KS, WR and  \nFig. 9. Schematic of physical components and associated tracers for the TMC1 protobinary system.  \n<!-- image -->  \nHL acknowledge support from ERC Advanced grant 101019751 MOLDISK, TOP-1 grant 614.001.751 from the Dutch Research Council (NWO), the Netherlands Research School for Astronomy (NOVA), the Danish National Research Foundation through the Center of Excellence 'InterCat' (DNRF150), and DFGgrant 325594231, FOR 2634 / 2. P.J.K. acknowledges financial support from the Science Foundation Ireland / Irish Research Council Pathway programme under Grant Number 21 / PATH-S / 9360. A.C.G. has been supported by PRIN-INAF MAIN-STREAM 2017 'Protoplanetary disks seen through the eyes of new generation instruments' and from PRIN-INAF 2019 'Spectroscopically tracing the disk dispersal evolution (STRADE)'. K.J. acknowledges the support from the Swedish National Space Agency (SNSA). G.P. gratefully acknowledges support from the Max Planck Society. This paper makes use of the following ALMA data: ADS / JAO.ALMA# 2017.1.01350.S, ADS / JAO.ALMA# 2017.1.01413.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 co-operation with the Republic of Chile. The Joint ALMA Observatory is operated by ESO, AUI / NRAO and NAOJ. Astrochemistry in Leiden is supported by the Netherlands Research School for Astronomy (NOVA). This research has made use of NASA's Astrophysics Data System Bibliographic Services. This research made use of NumPy (Harris et al. 2020); Astropy, a community-developed core Python package for Astronomy (Astropy Collaboration et al. 2013, 2018); Matplotlib (Hunter 2007); pdrtpy (Kaufman et al. 2006; Pound & Wolfire 2008, 2011, 2023)", 'References': 'Alcalá, J. M., Manara, C. F., Natta, A., et al. 2017, A&A, 600, A20  \nAlves, F. O., Caselli, P., Girart, J. M., et al. 2019, Science, 366, 90 Apai, D., Tóth, L. V., Henning, T., et al. 2005, A&A, 433, L33  \n- Arce, H. G., Borkin, M. A., Goodman, A. A., Pineda, J. E., & Halle, M. W. 2010, ApJ, 715, 1170\n- Argyriou, I., Glasse, A., Law, D. R., et al. 2023, A&A, 675, A111\n- Aso, Y., Ohashi, N., Saigo, K., et al. 2015, ApJ, 812, 27\n- Astropy Collaboration, Price-Whelan, A. M., Sip"ocz, B. M., et al. 2018, AJ, 156, 123\n- Astropy Collaboration, Robitaille, T. P., Tollerud, E. J., et al. 2013, A&A, 558, A33\n- Barsony, M., Wolf-Chase, G. A., Ciardi, D. R., & O\'Linger, J. 2010, ApJ, 720, 64\n- Beuther, H., van Dishoeck, E. F., Tychoniec, L., et al. 2023, A&A, 673, A121\n- Bjerkeli, P., van der Wiel, M. H. D., Harsono, D., Ramsey, J. P., & Jørgensen,\n- J. K. 2016, Nature, 540, 406 Boogert, A. C. A., Tielens, A. G. G. M., Ceccarelli, C., et al. 2000, A&A, 360, 683\n- Bottinelli, S., Boogert, A. C. A., Bouwman, J., et al. 2010, ApJ, 718, 1100 Bushouse, H., Eisenhamer, J., Dencheva, N., et al. 2023, JWST Calibration Pipeline, Zenodo  \nArticle number, page 12 of 17  \n- Cacciapuoti, L., Macias, E., Maury, A. J., et al. 2023, A&A, 676, A4\n- Chiar, J. E., Ennico, K., Pendleton, Y. J., et al. 2007, ApJ, 666, L73\n- Ching, T.-C., Lai, S.-P., Zhang, Q., et al. 2016, ApJ, 819, 159\n- Christiaens, V., Gonzalez, C., Farkas, R., et al. 2023, The Journal of Open Source Software, 8, 4774\n- Connelley, M. S. & Greene, T. P. 2010, AJ, 140, 1214\n- de Valon, A., Dougados, C., Cabrit, S., et al. 2020, A&A, 634, L12\n- Dionatos, O., Jørgensen, J. K., Teixeira, P. S., Güdel, M., & Bergin, E. 2014, A&A, 563, A28  \nDionatos, O., Nisini, B., Garcia Lopez, R., et al. 2009, ApJ, 692, 1  \n- Doppmann, G. W., Greene, T. P., Covey, K. R., & Lada, C. J. 2005, The Astronomical Journal, 130, 1145\n- Ehrenfreund, P., Boogert, A. C. A., Gerakines, P. A., Tielens, A. G. G. M., & van Dishoeck, E. F. 1997, A&A, 328, 649\n- Espaillat, C., Ingleby, L., Furlan, E., et al. 2013, ApJ, 762, 62\n- Farris, B. D., Du ff ell, P., MacFadyen, A. I., & Haiman, Z. 2014, ApJ, 783, 134  \nFederman, S., Megeath, S. T., Rubinstein, A. E., et al. 2023, arXiv e-prints,  \narXiv:2310.03803  \n- Fiorellino, E., Tychoniec, Ł., Cruz-Sáenz de Miera, F., et al. 2023, ApJ, 944, 135 Fiorellino, E., Tychoniec, Ł., Manara, C. F., et al. 2022, ApJ, 937, L9\n- Frank, A., Ray, T. P., Cabrit, S., et al. 2014, in Protostars and Planets VI, ed.\n- H. Beuther, R. S. Klessen, C. P. Dullemond, & T. Henning (Tucson, AZ: Univ. Arizona Press), 451-474\n- Fridlund, C. V. M. & Liseau, R. 1998, ApJ, 499, L75\n- Garufi, A., Podio, L., Codella, C., et al. 2022, A&A, 658, A104\n- Gasman, D., Argyriou, I., Sloan, G. C., et al. 2023, A&A, 673, A102\n- Giannini, T., Antoniucci, S., Nisini, B., Bacciotti, F., & Podio, L. 2015, ApJ, 814, 52\n- Gullbring, E., Hartmann, L., Briceño, C., & Calvet, N. 1998, ApJ, 492, 323 Gusdorf, A., Cabrit, S., Flower, D. R., & Pineau des Forêts, G. 2008a, A&A, 482, 809\n- Gusdorf, A., Pineau des Forêts, G., Cabrit, S., & Flower, D. R. 2008b, A&A, 490, 695\n- Harris, C. R., Millman, K. J., van der Walt, S. J., et al. 2020, Nature, 585, 357\n- Harsono, D., Bjerkeli, P., Ramsey, J. P., et al. 2023, ApJ, 951, L32\n- Harsono, D., Jørgensen, J. K., van Dishoeck, E. F., et al. 2014, A&A, 562, A77 Hartmann, L., Herczeg, G., & Calvet, N. 2016, ARA&A, 54, 135\n- Herbst, T. M., Beckwith, S. V. W., Glindemann, A., et al. 1996, AJ, 111, 2403\n- Hirano, N., Ho, P. P. T., Liu, S.-Y., et al. 2010, ApJ, 717, 58\n- Hollenbach, D. & Gorti, U. 2009, ApJ, 703, 1203\n- Hollenbach, D. & McKee, C. F. 1989, ApJ, 342, 306\n- Hunter, J. D. 2007, Computer Science and Engineering, 9, 90\n- Jones, O. C., Álvarez-Márquez, J., Sloan, G. C., et al. 2023, MNRAS, 523, 2519\n- Jørgensen, J. K., Schöier, F. L., & van Dishoeck, E. F. 2002a, A&A, 389, 908  \nJørgensen, J. K., Schöier, F. L., & van Dishoeck, E. F. 2002b, A&A, 389, 908  \n- Kaufman, M. J., Wolfire, M. G., & Hollenbach, D. J. 2006, ApJ, 644, 283\n- Keane, J. V., Tielens, A. G. G. M., Boogert, A. C. A., Schutte, W. A., & Whittet, D. C. B. 2001, A&A, 376, 254  \nKrolikowski, D. M., Kraus, A. L., & Rizzuto, A. C. 2021, AJ, 162, 110  \nKu ff meier, M., Jensen, S. S., & Haugbølle, T. 2023, European Physical Journal Plus, 138, 272  \nKwon, W., Fernández-López, M., Stephens, I. W., & Looney, L. W. 2015, ApJ,  \n814, 43 Labiano, A., Argyriou, I., Álvarez-Márquez, J., et al. 2021, A&A, 656, A57 Lahuis, F., van Dishoeck, E. F., Jørgensen, J. K., Blake, G. A., & Evans, N. J. 2010, A&A, 519, A3  \nLaw, D. D., Morrison, J. E., Argyriou, I., et al. 2023, AJ, 166, 45 Le Gouellec, V. J. M., Hull, C. L. H., Maury, A. J., et al. 2019, ApJ, 885, 106 Lee, C.-F., Codella, C., Li, Z.-Y., & Liu, S.-Y. 2019, ApJ, 876, 63 Louvet, F., Dougados, C., Cabrit, S., et al. 2018, A&A, 618, A120 Lynch, C. J. R., Smith, M. D., & Glover, S. C. O. 2020, MNRAS, 491, 3082 Maret, S., Bergin, E. A., Neufeld, D. A., et al. 2009, ApJ, 698, 1244 Matt, S. & Pudritz, R. E. 2005, ApJ, 632, L135 McClure, M. 2009, ApJ, 693, L81  \nMurillo, N. M., van Dishoeck, E. F., Hacar, A., Harsono, D., & Jørgensen, J. K. 2022, A&A, 658, A53  \nMurphy, G. C., Lery, T., O\'Sull§ivan, S., et al. 2008, A&A, 478, 453 Neufeld, D. A., Nisini, B., Giannini, T., et al. 2009, ApJ, 706, 170 Nisini, B., Antoniucci, S., Alcalá, J. M., et al. 2018, A&A, 609, A87 Nisini, B., Santangelo, G., Giannini, T., et al. 2015, ApJ, 801, 121 Nussbaumer, H. & Storey, P. J. 1988, A&A, 200, L25  \nO ff ner, S. S. R., Kratter, K. M., Matzner, C. D., Krumholz, M. R., & Klein, R. I. 2010, ApJ, 725, 1485  \nO  \nff  \nner, S. S. R., Moe, M., Kratter, K. M., et al. 2023, in Astronomical Society  \nof the Pacific Conference Series, Vol. 534, Protostars and Planets VII, ed.  \nS. Inutsuka, Y. Aikawa, T. Muto, K. Tomida, & M. Tamura, 275  \nPascucci,  \nI.,  \nCabrit,  \nS.,  \nEdwards,  \nS.,  \net  \nal.  \n2022,  \narXiv  \ne-prints,  \narXiv:2203.10068 Pineda, J. E., Segura-Cox, D., Caselli, P., et al. 2020, Nature Astronomy, 4, 1158 Plunkett, A. L., Arce, H. G., Corder, S. A., et al. 2013, ApJ, 774, 22 Podio, L., Bacciotti, F., Nisini, B., et al. 2006, A&A, 456, 189  \nPodio, L., Eislö ff el, J., Melnikov, S., Hodapp, K. W., & Bacciotti, F. 2011, A&A, 527, A13  \nPodio, L., Tabone, B., Codella, C., et al. 2021, A&A, 648, A45  \nPound, M. W. & Wolfire, M. G. 2008, in Astronomical Society of the Pacific  \nConference Series, Vol. 394, Astronomical Data Analysis Software and Sys-  \ntems XVII, ed. R. W. Argyle, P. S. Bunclark, & J. R. Lewis, 654  \nPound, M. W. & Wolfire, M. G. 2011, PDRT: Photo Dissociation Region Toolbox  \nPound, M. W. & Wolfire, M. G. 2023, AJ, 165, 25  \nRabenanahary, M., Cabrit, S., Meliani, Z., & Pineau des Forêts, G. 2022, A&A,  \n664, A118 Raghavan, D., McAlister, H. A., Henry, T. J., et al. 2010, ApJS, 190, 1 Rieke, G. H., Ressler, M. E., Morrison, J. E., et al. 2015, PASP, 127, 665 Rigby, J., Perrin, M., McElwain, M., et al. 2023, PASP, 135, 048001 Rigliaco, E., Pascucci, I., Duchene, G., et al. 2015, ApJ, 801, 31 Rocha, W. R. M., Rachid, M. G., Olsthoorn, B., et al. 2022, A&A, 668, A63 Romanova, M. M., Ustyugova, G. V., Koldoba, A. V., & Lovelace, R. V. E. 2005,  \nSakai, N., Sakai, T., Hirota, T., et al. 2014, Nature, 507, 78 Shang, H., Liu, C.-F., Krasnopolsky, R., & Wang, L.-Y. 2023, ApJ, 944, 230 Stahler, S. W. & Palla, F. 2005, The Formation of Stars, 865 Szulágyi, J., Pascucci, I., Ábrahám, P., et al. 2012, ApJ, 759, 47 Tabone, B., Cabrit, S., Bianchi, E., et al. 2017, A&A, 607, L6 Tabone, B., Godard, B., Pineau des Forêts, G., Cabrit, S., & van Dishoeck, E. F.  \nApJ, 635, L165 2020, A&A, 636, A60  \nTabone, B., van Hemert, M. C., van Dishoeck, E. F., & Black, J. H. 2021, A&A, 650, A192  \nTappe, A., Forbrich, J., Martín, S., Yuan, Y., & Lada, C. J. 2012, ApJ, 751, 9 Terebey, S., Van Buren, D., Brundage, M., & Hancock, T. 2006, ApJ, 637, 811 Timmermann, R., Bertoldi, F., Wright, C. M., et al. 1996, A&A, 315, L281 Tobin, J. J., Looney, L. W., Li, Z.-Y., et al. 2016, ApJ, 818, 73 Tobin, J. J., Sheehan, P. D., Megeath, S. T., et al. 2020, ApJ, 890, 130 Tokovinin, A. & Moe, M. 2020, MNRAS, 491, 5158 Tychoniec, Ł., Hull, C. L. H., Kristensen, L. E., et al. 2019, A&A, 632, A101 Tychoniec, Ł., van Dishoeck, E. F., van\'t Ho ff , M. L. R., et al. 2021, A&A, 655,  \nA65  \nValdivia-Mena, M. T., Pineda, J. E., Segura-Cox, D. M., et al. 2022, A&A, 667, A12  \nvan der Marel, N., van Dishoeck, E. F., Bruderer, S., et al. 2013, Science, 340, 1199  \nvan \'t Ho  \nff  \n, M. L. R., Harsono, D., Tobin, J. J., et al. 2020, ApJ, 901, 166  \nvan \'t Ho  \nff  \n, M. L. R., Tobin, J. J., Harsono, D., & van Dishoeck, E. F. 2018,  \nA&A, 615, A83 Watson, D. M., Calvet, N. P., Fischer, W. J., et al. 2016, ApJ, 828, 52 Weingartner, J. C. & Draine, B. T. 2001, ApJ, 548, 296 Wells, M., Pel, J. W., Glasse, A., et al. 2015, PASP, 127, 646 Wright, G. S., Rieke, G. H., Glasse, A., et al. 2023, PASP, 135, 048003 Wright, G. S., Wright, D., Goodson, G. B., et al. 2015, PASP, 127, 595 Yang, Y.-L., Green, J. D., Pontoppidan, K. M., et al. 2022, ApJ, 941, L13 Yousefi, M., Bernath, P. F., Hodges, J., & Masseron, T. 2018,  \nJ. Quant. Spectr. Rad. Transf., 217, 416 Zagaria, F., Rosotti, G. P., & Lodato, G. 2021, MNRAS, 504, 2235 Zasowski, G., Kemper, F., Watson, D. M., et al. 2009, ApJ, 694, 459', 'Appendix B: Tables': 'Appendix A.1: Extinction measurement  \nExtinction at the continuum position is measured with the optical depth of the silicate feature at 9.7 µ m. We show the fitting results in Fig. A.1. Polynomial is fit for the selected wavelength ranges: 4.9-5.3 µ m, 23.8-23.9 µ m, and 25-28 µ m. This allows us to avoid major broad absorption features such as water and silicates.  \nAt the outflow positions, the extinction correction is provided by simultaneously fitting rotational temperature and extinction value. Fig. A.2 shows the results of the fitting. We fit the value of K-band extinction A K . The relation between A λ and A K is provided following the extinction curve in McClure (2009) and convert to visual extinction with A V / A K = 7.75 for R V = 5.  \nFigure A.3 presents the intensity-weighted channel map of CN, as observed with ALMA. The velocity gradient in the image of increasing velocity with distance to the source is consistent with the infalling streamer.  \nThe MRS spectral resolving power varies from R = λ/ ∆ λ ∼ 1500 -3500 (Jones et al. 2023). In Fig. A.4, we present the emission from [Ni II], [Ar II], and H2. The emission from the ionized species appears to be velocity-resolved, with clear redand blueshifted components on the south and north sides of the source, respectively. On the other hand, H2 remains unresolved in the velocity regime, with only faint emission detected above 80 km s -1 .  \nFig. A.1. Spectra extracted from 2 x PSF FWHM on the continuum position (top) of TMC1-W (top-left) and TMC1-E (top-right). The dashed blue line represents a polynomial fit to the continuum. optical depth as a function of wavelength resulting from dividing the polynomial fit with the observed continuum (bottom). The orange dashed line represents τ = 0. The range of extincition is calculated using optical depth value at the silicate feature following the equations: AV = τ 9 . 7 × 18 . 5, and AV = τ 9 . 7 × 9, representing the correlation found in di ff use ISM and galactic centre dense clouds, respectively (McClure 2009; Weingartner & Draine 2001; Chiar et al. 2007; Boogert et al. 2000).  \n<!-- image -->  \n<!-- image -->  \nFig. A.2. Rotational diagrams for two positions: in the TMC1-W blueshifted jet positions (left) and TMC1-E blueshifted wind positions (right). Regions are indicated in Fig. 1.  \n<!-- image -->  \nFig. A.3. Intensity weighted velocity (moment 1) map of CN (2-1). The map reveals the gradient towards the protostar, consistent with infalling gas.  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nFig. A.4. Channel maps for [Ni II],[Ar II] and H2 S(3) showing blueshifted emission toward the north and weaker redshifted emission to the south.  \n<!-- image -->  \nTable B.1. Line intensities of hydrogen recombination lines at the positions of the protostars.Table B.2. Line intensities of H2 pure rotational transitions integrated at the outflow position of TMC1-E.  \n<!-- image -->  \n<!-- image -->  \n<!-- image -->  \nTable B.3. Line intensities at the outflow positionsTable B.4. Extinction-corrected line intensity ratios for selected lines'}