Outflows play a key role in the star and planet formation processes. Some outflows show discrete clumps of cold molecular gas moving at extremely high velocities (EHVs) of similar to 100 km s(-1), known as 'molecular bullets', that are likely closely associated with their primary driving agent. Here we present ALMA CO (J = 3-2) observations of a bright EHV molecular bullet that reveal its morphology in detail down to scales of 30 au and its kinematic structure across the entire intermediate velocity range (similar to 30-100 km s(-1)). These provide important insights into how outflows transfer mass and momentum to the surrounding medium. The observed channel maps display several sequences of ring-like features whose velocity increases and size decreases with projected distance from the driving source, each sequence tracing a thin, bow-shaped shell culminating on axis in a bright EHV head. The shape, kinematics and mass of each shell all agree remarkably well with the simplest textbook models of momentum-conserving bowshocks produced by a time-variable EHV jet. The dynamical timescale between consecutive shells is of a few decades, with the latest ejection event coinciding with the protostar optical/infrared outburst observed in similar to 1990. The very strong evidence for bowshock-driven entrainment induced by jet variability revealed by this work suggests that accretion bursts, and therefore variations in the disk snowlines, should occur on decade timescales, which could substantially impact grain growth and planet formation.

e report high angular resolution observations, made with the Atacama Large Millimeter Array in band 6, of high excitation molecular lines of CH3CN and SO2 and of the H29 alpha radio recombination line toward the G345.0061+01.794 B HC H ii region in order to investigate the physical and kinematical characteristics of its surroundings. Emission was detected in all observed components of the J = 14 -> 13 rotational ladder of CH3CN and in the 30(4,26)-30(3,27), and 32(4,28)-32(3,29) lines of SO2. The peak of the velocity-integrated molecular emission is located similar to 04 northwest of the peak of the continuum emission. The first-order moment images and channel maps show a velocity gradient of 1.1 km s(-1) arcsec(-1) across the source and a distinctive spot of blueshifted emission toward the peak of the zero-order moment. The rotational temperature is found to decrease from 252 +/- 24 K at the peak position to 166 +/- 16 K at its edge, indicating that our molecular observations are probing a hot molecular core that is internally excited. The emission in the H29 alpha line arises from a region of 065 in size, where its peak coincides with that of the dust continuum. We model the kinematical characteristics of the "central blue spot" feature as due to infalling motions, suggesting a central mass of 172.8 +/- 8.8 M-circle dot. Our observations indicate that this HC H ii region is surrounded by a compact structure of hot molecular gas, which is rotating and infalling toward a central mass, that is most likely confining the ionized region. The observed scenario is reminiscent of a "butterfly pattern" with an approximately edge-on torus and ionized gas roughly parallel to its rotation axis.

ABSTRACT A number of works infer radial temperature profiles of envelopes surrounding young stellar objects using several rotational transitions in a pixel-by-pixel or azimuthally averaged basis. However, in many cases the assumption that the rotational temperature is constant along the line of sight is made, while this is not the case when a partially resolved envelope, assumed to be spherically symmetric, is used to obtain values of temperature for different projected radii. This kind of analysis (homogeneous analysis) is intrinsically inconsistent. By using a spherical envelope model to interpret NH3 (1,1) and (2,2) observations, we tested how robust it is to infer radial temperature profiles of an envelope. The temperature and density of the model envelope are power laws of radius, but the density can be flat for an inner central part. The homogeneous analysis was applied to obtain radial temperature profiles, and resulted that for small projected radii, where the optical depth of the lines is high, the homogeneous temperature can be much higher than the actual envelope temperature. In general, for larger projected radii, both the temperature and the temperature power-law index can be underestimated by as much as 40 per cent, and 0.15, respectively. We applied this study to the infrared dark cloud G14.225–0.506 for which the radial temperature profile was previously derived from the dust emission at submillimetre wavelengths and the spectral energy distribution. As expected, the homogeneous analysis underestimated both the temperature and the temperature power-law index.

The infrared dark cloud (IRDC) G14.225-0.506 is part of the extended and massive molecular cloud located to the southwest of the region M17. The cloud is associated with a network of filaments, which result in two different dense hubs, as well as with several signposts of star formation activity and a rich population of protostars and young stellar objects (YSOs). The aim of this work is to study the centimeter continuum emission in order to characterize the stellar population in both regions, as well as to study the evolutionary sequence across the IRDC G14.225-0.506. We performed deep ($ mu radio continuum observations at 6 and 3.6 cm toward the IRDC G14.225-0.506 using the Karl G.\ Jansky Very Large Array (VLA) in its most extended A configuration (sim prime $). Data at both C and X bands were imaged using the same (u,v) range in order to derive spectral indices. We have also made use of observations taken during different days to study the presence of variability at short timescales toward the detected sources. We detected a total of 66 sources, 32 in the northern region G14.2-N and 34 in the southern region G14.2-S. Ten of the sources are found to be variable, with three located in G14.2-N and seven in G14.2-S. Based on their spectral index, the emission in G14.2-N is mainly dominated by nonthermal sources while G14.2-S contains more thermal emitters. Approximately 75 of the sources present a counterpart at other wavelengths. When considering the inner 0.4 pc region around the center of each hub, the number of infrared (IR) sources in G14.2-N is larger than in G14.2-S by a factor of 4. We also studied the relation between the radio luminosity and the bolometric luminosity, finding that the thermal emission of the studied sources is compatible with thermal radio jets. For our sources with X-ray counterparts, the nonthermal emitters follow a Güdel-Benz relation with $ = 0.03$, as previously suggested for other similar regions. We found similar levels of fragmentation between G14.2-N and G14.2-S, suggesting that both regions are most likely twin hubs. The nonthermal emission found in the less evolved objects, mainly coming from G14.2-N, suggests that G14.2-N may be composed of more massive YSOs as well as being in a more advanced evolutionary stage, consistent with the "filament-halo" gradient in age and mass from previous works. Overall, our results confirm a wider evolutionary sequence from the southwest to northeast starting in G14.2-S as the youngest part, followed by G14.2-N, and ending with the most evolved region M17.

A number of works infer radial temperature profiles of envelopes surrounding young stellar objects using several rotational transitions in a pixel-by-pixel or azimuthally averaged basis. However, in many cases the assumption that the rotational temperature is constant along the line of sight is made, while this is not the case when a partially resolved envelope, assumed to be spherically symmetric, is used to obtain values of temperature for different projected radii. This kind of analysis (homogeneous analysis) is intrinsically inconsistent. By using a spherical envelope model to interpret NH3 (1,1) and (2,2) observations, we tested how robust it is to infer radial temperature profiles of an envelope. The temperature and density of the model envelope are power laws of radius, but the density can be flat for an inner central part. The homogeneous analysis was applied to obtain radial temperature profiles, and resulted that for small projected radii, where the optical depth of the lines is high, the homogeneous temperature can be much higher than the actual envelope temperature. In general, for larger projected radii, both the temperature and the temperature power-law index can be underestimated by as much as 40 per cent, and 0.15, respectively. We applied this study to the infrared dark cloud G14.225-0.506 for which the radial temperature profile was previously derived from the dust emission at submillimetre wavelengths and the spectral energy distribution. As expected, the homogeneous analysis underestimated both the temperature and the temperature power-law index.

The IRDC G14.225-0.506 is associated with a network of filaments, which result in two different dense hubs, as well as with several signposts of star formation activity. The aim of this work is to study the cm continuum emission to characterize the stellar population in G14.2. We performed deep ( 1.5-3 microJy) radio continuum observations at 6 and 3.6 cm using the VLA in the A configuration ( 0.3”). We have also made use of observations taken during different days to study the presence of variability at short timescales. We detected a total of 66 sources, 32 in the northern region G14.2-N and 34 in the southern region G14.2-S. Ten of the sources are found to be variable. Based on their spectral index, the emission in G14.2-N is mainly dominated by non-thermal sources while G14.2-S contains more thermal emitters. Approximately 75 0.4 pc region around the center of each hub, the number of IR sources in G14.2-N is larger than in G14.2-S by a factor of 4. We also studied the relation between the radio luminosity and the bolometric luminosity, finding that the thermal emission of the studied sources is compatible with thermal radio jets. For our sources with X-ray counterparts, the non-thermal emitters follow a Güdel-Benz relation with k = 0.03. We found similar levels of fragmentation between G14.2-N and G14.2-S, suggesting that both regions are most likely twin hubs. The non-thermal emission found in the less evolved objects suggests that G14.2-N may be composed of more massive YSOs as well as being in a more advanced evolutionary stage, consistent with the filament-halo gradient in age and mass from previous works. Our results confirm a wider evolutionary sequence starting in G14.2-S as the youngest part, followed by G14.2-N, and ending with the most evolved region M17.

We present new Atacama Large Millimeter/submillimeter Array 0.98 mm observations of the continuum emission and several molecular lines toward the high-mass protostellar system GGD 27-MM1, driving the HH 80-81 radio jet. The detailed analysis of the continuum and the CH _3 CN molecular emission allows us to separate the contributions from the dust content of the disk (extending up to 190 au), the molecular content of the disk (extending from 140–360 au), and the content of the envelope, revealing the presence of several possible accretion streamers (also seen in other molecular tracers, such as CH _3 OH). We analyze the physical properties of the system, producing temperature and column density maps, and radial profiles for the disk and the envelope. We qualitatively reproduce the trajectories and line-of-sight velocities of the possible streamers using a theoretical model approach. An ad hoc model of a flared disk comprising a hot dust disk embedded in cold gas fits the H _2 S emission, which revealed the molecular disk as a crescent shape with a prominent central absorption. Another fit to the central absorption spectrum suggests that the absorption is probably caused by different external cold layers from the envelope or the accretion streamers. Finally, the analysis of the rotation pattern of the different molecular transitions in the molecular disk suggests that there is an inner zone devoid of molecular content.

Abridged. The bright-rimmed cloud IC1396N hosts CO, H$_2$, and Herbig-Haro outflows powered by millimetre compact sources. We aim to characterise the kinematics and physical conditions of the H$_2$ emission features spread over IC1396N, which appear as chains of knots with a jet-like morphology, tracing different H$_2$ outflows, and to obtain further information about (and an identification of) the driving sources. Low-resolution, long-slit near-infrared spectra were acquired with NICS at the TNG, using grisms KB (R~1200), HK and JH (R~500). Several slit pointings and PA were used throughout IC1396N to sample a number of the H$_2$ knots previously detected in deep H$_2$ 2.12 $\mu$m images. The knots exhibit rich ro-vibrational spectra of H$_2$, consistent with shock-excited excitation, from which radial velocities and relevant physical conditions of the IC1396N H$_2$ outflows were derived. These also allowed estimating extinction ranges towards several features. [FeII] emission was only detected towards a few knots, which also display unusually large H$_2$ 1-0 S(3)/S(1) flux ratios. The obtained radial velocities confirm that most of the outflows are close to the plane of the sky. Nearby knots in the same chain often display different radial velocities, both blue-shifted and red-shifted, which we interpret as due to ubiquitous jet precession in the driving sources or the development of oblique shocks. One of the chains (strand A) appears as a set of features trailing a leading bow-shock structure consistent with the results of 3-D magneto-hydrodynamical models. Either side of the leading bow-shock (A15) exhibits a different radial velocity, whose possible explanations are discussed in the paper. Our data cannot confirm whether strands A and B have both been originated by the intermediate mass young stellar object BIMA 2.

Context. The bright-rimmed cloud IC 1396N is believed to host one of the few known cases where two bipolar CO outflows driven by young stellar objects collide. The CO outflows are traced by chains of knots of H2 emission, with enhanced emission at the position of the possible collision. Aims. The aim of this work is to use the proper motions of the H2 knots to confirm the collision scenario. Methods. A second-epoch H2 image was obtained, and the proper motions of the knots were determined with a time baseline of ~11 yr. We also performed differential photometry on the images to check the flux variability of the knots. Results. For each outflow (N and S), we classified the knots as pre-collision or post-collision. The axes of the pre-collision knots, the position of the possible collision point, and the axes of the post-collision knots were estimated. The difference between the proper motion direction of the post-collision knots and the position angle from the collision point was also calculated. For some of the knots, we obtained the 3D velocity using the radial velocity derived from H2 spectra. Conclusions. The velocity pattern of the H2 knots in the area of interaction (post-collision knots) shows a deviation from that of the pre-collision knots, consistent with being a consequence of the interaction between the two outflows. This favours the interpretation of the IC 1396N outflows as a true collision between two protostellar jets instead of a projection effect.

We present Very Large Array (VLA) and Atacama Large Millimeter/submillimeter Array (ALMA) observations of the close (0.″3 = 90 au separation) protobinary system SVS 13. We detect two small circumstellar disks (radii ∼12 and ∼9 au in dust, and ∼30 au in gas) with masses of ∼0.004–0.009 M ☉ for VLA 4A (the western component) and ∼0.009–0.030 M ☉ for VLA 4B (the eastern component). A circumbinary disk with prominent spiral arms extending ∼500 au and a mass of ∼0.052 M ☉ appears to be in the earliest stages of formation. The dust emission is more compact and with a very high optical depth toward VLA 4B, while toward VLA 4A the dust column density is lower, allowing the detection of stronger molecular transitions. We infer rotational temperatures of ∼140 K, on scales of ∼30 au, across the whole source, and a rich chemistry. Molecular transitions typical of hot corinos are detected toward both protostars, being stronger toward VLA 4A, with several ethylene glycol transitions detected only toward this source. There are clear velocity gradients, which we interpret in terms of infall plus rotation of the circumbinary disk, and pure rotation of the circumstellar disk of VLA 4A. We measured orbital proper motions and determined a total stellar mass of 1 M ☉ . From the molecular kinematics, we infer the geometry and orientation of the system, and stellar masses of ∼0.26 M ☉ for VLA 4A and ∼0.60 M ☉ for VLA 4B.

The early stages of the formation of binary stellar systems are still rather poorly understood observationally, which contrasts with some significant recent improvements in numerical simulations. We present a comprehensive study with the VLA and ALMA of the close (separation = 90 au) proto-binary system SVS 13. Our very high sensitivity and spatial resolution observations trace the dust as well as the ionized and molecular gas in this system, reaching scales as small as ~10 au. We infer the orbital motion and masses of the two protostars. We image two circumstellar disks and a still-forming circumbinary disk with prominent spiral arms extending ~500 au. We study the 3D kinematics of the system and measure the physical properties of the disks. We also find evidence for variation of chemical properties on scales of a few tens of au. Finally, we will discuss on how the properties of the SVS 13 system compare with those of a few other protobinary systems that have been observed with a similar degree of detail, and with the predictions of numerical simulations. This kind of information provides some clues on the final configuration of planetary systems in binary systems. ALMA observation of dust in the disks around SVS 13 (credit: A.K. Diaz-Rodriguez, G. Anglada 2022). Cartoon model of the system. The red-blue colours indicate the motion of the gas. Red – away from us, blue – towards us. The peculiar yin-yang shape results from the combination of infalling and rotation motions (credit: A. K. Diaz-Rodriguez et al. 2022).

Aims. The infall of material onto a protostar, in the case of optically thick line emission, produces an asymmetry in the blue- and red-wing line emissions. For an angularly resolved emission, this translates in a blue central spot in the first-order moment (intensity weighted velocity) map. Methods. An analytical expression for the first-order moment intensity as a function of the projected distance was derived, for the cases of infinite and finite infall radius. The effect of a finite angular resolution, which requires the numerical convolution with the beam, was also studied. Results. This method was applied to existing data of several star-forming regions, namely G31.41+0.31 HMC, B335, and LDN 1287, obtaining good fits to the first-order moment intensity maps, and deriving values of the central masses onto which the infall is taking place (G31.41+0.31 HMC: 70–120 M ⊙ ; B335: 0.1 M ⊙ ; Guitar Core of LDN 1287: 4.8 M ⊙ ). The central-blue-spot infall hallmark appears to be a robust and reliable indicator of infall.

Jets and outflows are associated with young stellar objects across the stellar mass spectrum, from brown dwarf protostars to massive, Ae/Be stars. Frequently, the jet morphology is spatially discontinuous because of the temporal variability of the ejection from the driving source. Images covering a wide field of view around the jet driving-source are useful to map the large-scale jet emission and to explore the mass ejection history. The aim of this work was to search for large-scale optical Herbig-Haro (HH) objects lying in a wide field around a sample of IRAS sources, candidates to trace young stellar objects. Deep, narrow-band images through the H$\alpha$ and [SII] emission lines, and through an off-line continuum filter, covering a wide ($\sim15'$) field around the IRAS targets were acquired. The images in the three filters were analyzed to identify shock-excited line emission (i.e., HH) in contrast to scattered line emission. New images of a sample of fifteen IRAS sources, obtained in an homogeneous way are presented. HH emission was detected in six fields, and the astrometry of the knot features is given. The nature of the extended emission as scattered emission around nine of the IRAS targets is confirmed. For seven IRAS sources, with unclear optical counterpart, a more plausible counterpart is proposed. A refined value of the source distance is reported for seven targets. An update of the main data available for each of the sampled fields, including images from public data archives, is also presented.

Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here, we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. to build a sample of 18 massive dense cores where both fragmentation and magnetic field properties are studied in a uniform way. We measured the fragmentation level, N mm, within the field of view common to all regions of ∼0.15 pc, with a mass sensitivity of ∼0.5 M ☉, and a spatial resolution of ∼1000 au. In order to obtain the magnetic field strength using the Davis–Chandrasekhar–Fermi method, we estimated the dispersion of the polarization position angles, the velocity dispersion of the H13CO+(4–3) gas, and the density of each core, all averaged within 0.15 pc. A strong correlation is found between N mm and the average density of the parental core, although with significant scatter. When large-scale systematic motions are separated from the velocity dispersion and only the small-scale (turbulent) contribution is taken into account, a tentative correlation is found between N mm and the mass-to-flux ratio, as suggested by numerical and theoretical works.

Recent high angular resolution (≃40 mas) ALMA observations at 1.14 mm resolve a compact (R ≃ 200 au), flattened dust structure perpendicular to the HH 80–81 jet emanating from the GGD 27-MM1 high-mass protostar, making it a robust candidate for a true accretion disk. The jet–disk system (HH 80–81/GGD 27-MM1) resembles those found in association with low- and intermediate-mass protostars. We present radiative transfer models that fit the 1.14 mm ALMA dust image of this disk, which allow us to obtain its physical parameters and predict its density and temperature structure. Our results indicate that this accretion disk is compact (Rdisk ≃ 170 au) and massive (≃5 M⊙), at about 20% of the stellar mass of ≃20 M⊙. We estimate the total dynamical mass of the star–disk system from the molecular line emission, finding a range between 21 and 30 M⊙, which is consistent with our model. We fit the density and temperature structures found by our model with power-law functions. These results suggest that accretion disks around massive stars are more massive and hotter than their low-mass siblings, but they still are quite stable. We also compare the temperature distribution in the GGD 27–MM1 disk with that found in low- and intermediate-mass stars and discuss possible implications for the water snow line. We have also carried out a study of the distance based on Gaia DR2 data and the population of young stellar objects in this region and from the extinction maps. We conclude that the source distance is within 1.2 and 1.4 kpc, closer than what was derived in previous studies (1.7 kpc).

We present sensitive and high angular resolution ($\sim$0.2-0.3$''$) (sub)millimeter (230 and 345 GHz) continuum and CO(2$-$1)/CO(3$-$2) line archive observations of the disk star system in UX Tauri carried out with ALMA (The Atacama Large Millimeter/Submillimeter Array). These observations reveal the gas and dusty disk surrounding the young star UX Tauri A with a large signal-to-noise ratio ($>$400 in the continuum and $>$50 in the line), and for the first time is detected the molecular gas emission associated with the disk of UX Tauri C (with a size for the disk of $<$56 au). No (sub)millimeter continuum emission is detected at 5$\sigma$-level (0.2 mJy at 0.85 mm) associated with UX Tauri C. For the component UX Tauri C, we estimate a dust disk mass of $\leq$ 0.05 M$_\oplus$. Additionally, we report a strong tidal disk interaction between both disks UX Tauri A/C, separated 360 au in projected distance. The CO line observations reveal marked spiral arms in the disk of UX Tauri A and an extended redshifted stream of gas associated with the UX Tauri C disk. No spiral arms are observed in the dust continuum emission of UX Tauri A. Assuming a Keplerian rotation we estimate the enclosed masses (disk$+$star) from their radial velocities in 1.4 $\pm$ 0.6 M$_\odot$ for UX Tauri A, and 70 $\pm$ 30 / $\sin i$ Jupiter masses for UX Tauri C (the latter coincides with the mass upper limit value for a brown dwarf). The observational evidence presented here lead us to propose that UX Tauri C is having a close approach of a possible wide, evolving and eccentric orbit around the disk of UX Tauri A causing the formation of spiral arms and the stream of molecular gas falling towards UX Tauri C.

Aims. In this paper, we study the dense gas of the molecular cloud LDN 1287 (L1287), which harbors a double FU Ori system, an energetic molecular outflow, and a still-forming cluster of deeply embedded low-mass young stellar objects that show a high level of fragmentation. Methods. We present optical H alpha and [SII], and VLA NH3 (1, 1) and (2, 2) observations with an angular resolution of similar to 3 '.5. The observed NH3 spectra have been analyzed with the Hyperfine Structure tool, fitting simultaneously three different velocity components. Results. The NH3 emission from L1287 comes from four different structures: a core associated with RNO 1, a guitar-shaped core (the Guitar) and two interlaced filaments (the blue and red filaments) roughly centered toward the binary FU Ori system RNO 1B/1C and its associated cluster. Regarding the Guitar core, there are clear signatures of gas infall onto a central mass that has been estimated to be similar to 2.1M(circle dot). Regarding the two filaments, they have radii of similar to 0.03 pc, masses per unit length of similar to 50M(circle dot) pc(-1), and are in near isothermal equilibrium. A central cavity is identified, probably related with the outflow and also revealed by the H alpha and [SII] emission, with several young stellar objects near its inner walls. Both filaments show clear signs of perturbation by the high-velocity gas of the outflows driven by one or several young stellar objects of the cluster. The blue and red filaments are coherent in velocity and have nearly subsonic gas motions, except at the position of the embedded sources. Velocity gradients across the blue filament can be interpreted either as infalling material onto the filament or rotation. Velocity gradients along the filaments are interpreted as infall motions toward a gravitational well at the intersection of the two filaments.

Theoretical and numerical works indicate that a strong magnetic field should suppress fragmentation in dense cores. However, this has never been tested observationally in a relatively large sample of fragmenting massive dense cores. Here we use the polarization data obtained in the Submillimeter Array Legacy Survey of Zhang et al. (2014) to build a sample of 17 massive dense cores where both fragmentation and magnetic field properties are studied in a uniform way. We measured the fragmentation level, Nmm, within the field of view common to all regions, of ~0.15 pc, with a mass sensitivity of ~0.5 Msun, and a spatial resolution of ~1000 AU. In order to obtain the magnetic field strength using the Davis-Chandrasekhar-Fermi method, we estimated the dispersion of the polarization position angles, the velocity dispersion of the H13CO+(4-3) gas, and the density of each core, all averaged within 0.15 pc. The magnetic field strength was also inferred using the Angular Dispersion Function method. Surprisingly, we found no apparent relation between the fragmentation level and the magnetic field strength, while a possible trend of Nmm with the average density of the parental core was found. In addition, the average masses of the fragments are comparable to the thermal Jeans mass. Therefore, our results suggest that thermal fragmentation and gravity dominate the fragmentation process in this sample. However, when only cores with similar densities are considered, there are hints of a possible anticorrelation between the fragmentation level and the magnetic field strength.

Aims. In this paper, we study the dense gas of the molecular cloud LDN 1287 (L1287), which harbors a double FU Ori system, an energetic molecular outflow, and a still-forming cluster of deeply embedded low-mass young stellar objects that show a high level of fragmentation. Methods. We present optical H α and [S II ], and VLA NH 3 (1, 1) and (2, 2) observations with an angular resolution of ~3′′.5. The observed NH 3 spectra have been analyzed with the Hyperfine Structure tool, fitting simultaneously three different velocity components. Results. The NH 3 emission from L1287 comes from four different structures: a core associated with RNO 1, a guitar-shaped core (the Guitar) and two interlaced filaments (the blue and red filaments) roughly centered toward the binary FU Ori system RNO 1B/1C and its associated cluster. Regarding the Guitar core, there are clear signatures of gas infall onto a central mass that has been estimated to be ~2.1 M ⊙ . Regarding the two filaments, they have radii of ~0.03 pc, masses per unit length of ~50 M ⊙ pc −1 , and are in near isothermal equilibrium. A central cavity is identified, probably related with the outflow and also revealed by the H α and [S II ] emission, with several young stellar objects near its inner walls. Both filaments show clear signs of perturbation by the high-velocity gas of the outflows driven by one or several young stellar objects of the cluster. The blue and red filaments are coherent in velocity and have nearly subsonic gas motions, except at the position of the embedded sources. Velocity gradients across the blue filament can be interpreted either as infalling material onto the filament or rotation. Velocity gradients along the filaments are interpreted as infall motions toward a gravitational well at the intersection of the two filaments.