# Accretion versus core-filament collision: Implications for streamer formation in Per-emb-2

Fumitaka Nakamura, Quang Nguyen-Luong, Kousuke Ishihara, Aoto Yoshino

► To cite this version:

Fumitaka Nakamura, Quang Nguyen-Luong, Kousuke Ishihara, Aoto Yoshino. Accretion versus core-filament collision: Implications for streamer formation in Per-emb-2. Astronomy & Astrophysics - A&A, 2024, 690, 10.1051/0004-6361/202450854. insu-04843510

HAL Id: insu-04843510

https://insu.hal.science/insu-04843510v1

Submitted on 17 Dec 2024

HAL is a multi-disciplinary open access
archive for the deposit and dissemination of sci-
entific research documents, whether they are pub-
lished or not. The documents may come from
teaching and research institutions in France or
abroad, or from public or private research centers.

L'archive ouverte pluridisciplinaire HAL, est
destinée au dépôt et à la diffusion de documents
scientifiques de niveau recherche, publiés ou non,
émanant des établissements d'enseignement et de
recherche français ou étrangers, des laboratoires
publics ou privés.

<figure><img src="image_1.png" /></figure>

Distributed under a Creative Commons Attribution 4.0 International License
LETTER TO THE EDITOR

**Accretion versus core-filament collision**

**Implications for streamer formation in Per-emb-2**

Fumitaka Nakamura<fnref n="1" />,<fnref n="2" />,<fnref n="3" />,*, Quang Nguyen-Luong<fnref n="4" />,<fnref n="5" />,<fnref n="6" />, Kousuke Ishihara<fnref n="1" />,<fnref n="2" />, and Aoto Yoshino<fnref n="1" />,<fnref n="3" />

Received 24 May 2024 / Accepted 3 September 2024

ABSTRACT

*Context.* Recent millimetre and sub-millimetre observations have unveiled elongated and asymmetric structures around protostars. These structures, referred to as streamers, often exhibit coherent velocity gradients, seemingly indicating a directed gas flow towards the protostars. However, their origin and role in star formation remain uncertain.

*Aims.* The protostellar core Per-emb-2, located in Barnard 1, has a relatively large streamer of $10^4$ au that is more prominent in emission from carbon-chain molecules. We aim to unveil the formation mechanism of this streamer.

*Methods.* We conducted mapping observations towards Per-emb-2 using the Nobeyama 45 m telescope. We targeted carbon-chain molecular lines such as CCS, $HC_3N$, and $HC_5N$ at 45 GHz.

*Results.* Using *astrodendro*, we identified one protostellar and four starless cores, including three new detections, on the *Herschel* column density map. The starless and protostellar cores are more or less gravitationally bound. We discovered strong CCS and $HC_3N$ emissions extending from the north to the south, appearing to bridge the gap between the protostellar core and the starless core to its north. This bridge spans $3 \times 10^4$ au with velocities of 6.5–7.0 km s$^{-1}$. The velocity gradient of the bridge is opposite that of the streamer. Thus, the streamer is unlikely to be connected to the bridge, suggesting that the former does not have an accretion origin.

*Conclusions.* We propose that a collision between a spherical core and the filament has shaped the density structure in this region, consequently triggering star formation within the head-tail-shaped core. In this core-filament collision scenario, the collision appears to have fragmented the filament into two structures. The streamer is a bow structure, while the bridge is a remnant of the shock-compressed filament. Thus, we conclude that the Per-emb-2 streamer does not significantly contribute to the mass accumulation towards the protostar.

**Key words.** stars: formation – stars: protostars – ISM: clouds – ISM: kinematics and dynamics – ISM: structure

**1. Introduction**

Stars form from dense cores in molecular clouds. According to the standard scenario of star formation (Shu et al. 1987), a nearly spherical dense core undergoes gravitational contraction, forming a nearly axi-symmetric system consisting of a protostar and a rotationally supported disk. However, recent observations have uncovered significant non-axisymmetric structures surrounding young protostellar systems, commonly referred to as streamers or spirals (Pineda et al. 2020; Sanhueza et al. 2019; Chen et al. 2021; Valdivia-Mena et al. 2022; Hsieh et al. 2023; Olguin et al. 2023). Hereafter, we use the term ‘streamer’ to refer to such structures for the sake of simplicity. These structures often exhibit coherent velocity patterns, interpreted as infall towards the central protostar (Pineda et al. 2020; Sanhueza et al. 2021). Understanding the role of streamers in mass loading towards central stars is crucial for unveiling the complexities of the stellar mass accumulation process.

Numerical studies have proposed various mechanisms for streamer formation. For example, in many numerical simulations of cluster formation, intermediate-density structures connected to dense cores and sink particles, which resemble the streamers, are often observed, and these structures play a role in mass accretion onto the central protostellar system (e.g. Bonnell & Bate 2006; Nakamura & Li 2007; Kuffmeier et al. 2017; Padoan et al. 2020). On the other hand, Tu et al. (2024) demonstrate that streamer-like structures could arise from local turbulent compression within magnetised parent cores, without the need for external accretion of ambient gas. In this case, the protostellar system does not need to have structures connected to the inter-core medium, and streamers can be generated through the accretion process within the natal core. Therefore, the accretion of ambient gas likely plays only a minor role in supplying the gas, and the stellar mass is roughly determined by the parent core mass. Yano et al. (2024) propose that streamers could form through dense core collisions (DCCs; see also Kinoshita & Nakamura 2022). They estimated the frequencies of DCCs in nearby star-forming molecular clouds and found

<footnote n="1">¹ National Astronomical Observatory of Japan, National Institutes of Natural Sciences, 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan</footnote>

<footnote n="2">² Department of Astronomical Science, SOKENDAI (The Graduate University for Advanced Studies), 2-21-1 Osawa, Mitaka, Tokyo 181-8588, Japan</footnote>

<footnote n="3">³ Department of Astronomy, Graduate School of Science, The University of Tokyo, 7-3-1 Hongo, Bunkyo-ku, Tokyo 113-0033, Japan</footnote>

<footnote n="4">⁴ Astrophysical Data Sciences, CSMES, the American University of Paris, PL111, 2 bis, passage Landrieu 75007, Paris, France</footnote>

<footnote n="5">⁵ Université Paris-Saclay, Université Paris Cité, CEA, CNRS, AIM, 91191 Gif-sur-Yvette, France</footnote>

<footnote n="6">⁶ TNU Observatory, Tay Nguyen University, 567 Lé Duan, Ea Tam, Buôn Ma Thuột, Dăk Lăk 630000, Vietnam</footnote>

<footnote>* Corresponding author; fumitaka.nakamura@ao.ac.jp</footnote>
<figure><img src="image_2.png" /><figcaption>**Fig. 1.** Overall distributions of molecular gas. Panel (a): Observed area as seen on the Herschel hydrogen molecule column density map. The black box is the 5' × 5' observation box. The five enclosed red contours indicate the dense cores identified by astrodendro (see Table 1 for the core numbers). The yellow dot located at the centre is the position of the protostar. The blue dots show the positions of the starless cores identified by Pezzuto et al. (2021). Panels (b)–(d): Velocity-integrated maps (from 6.0 to 8.5 km s⁻¹) of CCS (J_N = 4₃ − 3₂) (b), HC₃N (J = 5 − 4) (c), and HC₅N (J = 17 − 16) (d) observed with the NRO 45 m telescope. The red lines follow the large-scale filament in this area. The filament is more prominent in a map that covers a larger area (see Fig. D.1). The blue ellipses in panels (b), (c), and (d) are the map effective beam size.</figcaption></figure>

that the typical dense core experiences at least one collision in
its lifetime, particularly in clustered environments. In this DCC
model, the compressed layer created between two cores evolves
into streamers and spirals, creating a rotating envelope around
the protostellar system. In this case, the streamer supplies addi-
tional mass to the central protostar, but the mass accretion occurs
only temporarily and depends on conditions such as the impact
parameter and collision speed. However, observational evidence
is needed to further elucidate these mechanisms.

To shed light on streamer formation, we conducted obser-
vations towards Per-emb-2, a core with a streamer towards
the central protostar (Pineda et al. 2020). It is located in the
southern part of Barnard 1 in the Perseus molecular cloud
at a distance of 300 pc from the Sun (Zucker et al. 2019).
Gas within the streamer appears to exhibit a velocity gradi-
ent towards the central protostar – from 7.5 km s⁻¹ (far side) to
7.0 km s⁻¹ (protostar) – and the streamer is particularly promi-
nent in emissions from chemically early-phase molecules such
as CCS and HC₃N, as observed with the Northern Extended Mil-
limetre Array (NOEMA) interferometer. Since the abundance
of these molecules decreases steeply with time in high-density
environments, Pineda et al. (2020) claim that the gas of the
streamer is supplied from outside via accretion. Very recently,
Taniguchi et al. (2024) suggested that the ambient gas may be
provided from the northern part of the protostar.

We present detailed observational information in Sect. 2 and
the analysis of the velocity and spatial distributions of molecu-
lar gas obtained from these observations in Sect. 3. In Sect. 4
we propose our core-filament collision (CFC) model, which
explains the structure in the observed region.

**2. Observations and data**

We carried out the on-the-fly mapping observations towards
the protostellar core Per-emb-2 using the Z45 receiver of
the Nobeyama Radio Observatory (NRO) 45m telescope
(Nakamura et al. 2014) in the 2022–2023 season. Details of the
observations are presented in Appendix A. The mapped area was
determined such that the protostar is located at the centre of the
5' × 5' box (Fig. 1a).

The target lines are CCS (J<sub>N</sub> = 4<sub>3</sub> − 3<sub>2</sub>), HC<sub>3</sub>N (J = 5 − 4), and HC<sub>5</sub>N (J = 17 − 16). The parameters are summarised in Table A.1. See Nakamura et al. (2024) for the accurate rest frequency of CCS (J<sub>N</sub> = 4<sub>3</sub> − 3<sub>2</sub>). The telescope has a full width at half maximum beam size of 38'' at 43 GHz. The standard chopper wheel method was used to convert the output signal to the antenna temperatures (T<sub>A</sub>*) and to correct for the atmospheric attenuation. The main beam efficiency of Z45 is η<sub>Z45</sub> ≈ 0.73 at 45 GHz. The main-beam temperature was calculated as T<sub>mb</sub> = T<sub>A</sub>*/η<sub>Z45</sub>. We adopted a spheroidal function as a gridding convolution function to calculate the intensity at each grid point of the final cube data with a spatial grid size of 15'' and a frequency resolution of 3.81 kHz, which corresponds to ~0.025 km s<sup>-1</sup> at 45 GHz. The final effective angular resolution of the map was 49'', corresponding to 15 000 au, or 0.07 pc at the distance of 300 pc (Zucker et al. 2019).

We downloaded the high-resolution H₂ column density and dust temperature maps obtained by the Herschel Gould Belt Survey (André et al. 2010; Pezzuto et al. 2021). The angular resolution of the maps is 18''.

**3. Results and analysis**

**3.1. Dense core identification and their dynamical states**

In Fig. 1a we present the Herschel column density map of the observed area. We applied astrodendro (Rosolowsky et al. 2008) to this image and identified five cores, four of which are fully within our observation box (see Appendix B.1 for more details). The boundaries of the identified cores are outlined in red. Their physical quantities are summarised in Table 1. In the same area, Pezzuto et al. (2021) identified three cores using get-sources: a protostellar core (HGBS 447) and two starless cores (HGBSs 442 and 449). Our cores 1 and 3 correspond to HGBS 442 and 447, respectively. Cores 2, 4, and 5 are starless and newly detected. We could not identify HGBS 449, which has the lowest column density. The central protostellar core (core 3) has a mass of ~10 M<sub>☉</sub>. According to Pineda et al. (2020), the protostellar core has a systemic velocity of ~7.0 km s⁻¹ and has a streamer with a length of 10⁴ au (Pineda et al. 2020).
<table><caption>**Table 1.** Properties of cores detected in the mapping region by **astrodendro**.</caption><thead><tr><th>Number</th><th>Peak RA<br>(J2000)</th><th>Peak Dec<br>(J2000)</th><th><i>M</i><sub>core</sub><br>[M<sub>☉</sub>]</th><th><i>R</i><sub>core</sub><br>[pc]</th><th>⟨<i>N</i>⟩<br>[10<sup>22</sup> cm<sup>-2</sup>]</th><th>⟨<i>n</i>⟩<br>[10<sup>5</sup> cm<sup>-3</sup>]</th><th>⟨<i>T</i><sub>dust</sub>⟩<br>[K]</th><th>σ<sub>CCS</sub><br>[km s<sup>-1</sup>]</th><th>α<sub>vir</sub>(α<sub>BE</sub>)</th><th>Type/HGBS</th></tr></thead><tbody><tr><td>1</td><td>3:32:19.13</td><td>30:51:44.67</td><td>0.508</td><td>0.0148</td><td>1.28</td><td>5.4</td><td>12.7</td><td>0.37</td><td>1.422 (0.565)</td><td>Starless/449</td></tr><tr><td>2</td><td>3:32:16.36</td><td>30:51:08.40</td><td>0.398</td><td>0.0077</td><td>1.25</td><td>29.6</td><td>12.6</td><td>0.41</td><td>1.014 (0.377)</td><td>Starless/new</td></tr><tr><td>3</td><td>3:32:18.50</td><td>30:49:50.63</td><td>9.723</td><td>0.0884</td><td>1.97</td><td>0.5</td><td>12.1</td><td>0.84</td><td>1.077 (0.176)</td><td>Protostellar/447</td></tr><tr><td>4</td><td>3:32:26.94</td><td>30:48:21.45</td><td>0.200</td><td>0.0130</td><td>0.64</td><td>3.1</td><td>13.8</td><td>–</td><td>– (1.265)</td><td>Starless/new</td></tr><tr><td>5</td><td>3:32:06.46</td><td>30:47:55.46</td><td>1.083</td><td>0.0565</td><td>0.86</td><td>0.2</td><td>12.7</td><td>0.45 <sup>(a)</sup></td><td>2.954 (1.013)</td><td>Starless/new</td></tr><tr><td>HGBS 442</td><td>3:32:12.04</td><td>30:47:59.0</td><td>0.105</td><td>0.0300</td><td>0.25</td><td>0.05</td><td>0.02</td><td>10.4</td><td>(3.9)</td><td>Starless</td></tr></tbody></table>

**Notes.** The core numbers are from the *Herschel* core catalogue of Pezzuto et al. (2021). <sup>(a)</sup>We estimated this velocity dispersion only within the area overlapping with the observation box.

From the CCS cube data, we derived the core’s CCS velocity dispersion (σ<sub>CCS</sub>) and then converted it to the intrinsic velocity dispersion (σ<sub>tot</sub>) of gas with a mean molecular weight of 2.33m<sub>H</sub> to assess its dynamical state (see Appendix B.2 for the details of the analysis). Cores 1, 2, 3, and 5 have α<sub>vir</sub> ~ 1; the virial parameter α<sub>vir</sub> was computed under the assumption of a centrally condensed (×r<sup>-2</sup>) sphere. The protostellar core has a significantly small Bonnor-Ebert ratio of α<sub>BE</sub> ~ 0.18, but its α<sub>vir</sub> is close to unity (the Bonnor-Ebert ratio is the ratio between the Bonnor-Ebert critical mass of the core and the observed mass). Core 4 does not have significant CCS emission, and we could not measure the velocity dispersion, but its Bonnor-Ebert ratio (α<sub>BE</sub> ~ 1.3) suggests it is close to bound. We also calculated the external pressures exerted on the core surfaces. We find that they play a minor role except for core 5, for which we calculated the core’s external pressure (3ρ<sub>amb</sub>σ<sup>2</sup><sub>tot,amb</sub>) from the ambient density (ρ<sub>amb</sub>) and the ambient gas intrinsic velocity dispersion (σ<sub>tot,amb</sub>) measured in the branch just below the corresponding core. The branch is an intermediate structure identified by **astrodendro**. Core 5 is dynamically compressed by the ambient pressure and likely did not form via spontaneous gravitational fragmentation, a process that leads to cores becoming self-gravitating. The other cores are close to being in dynamical equilibrium.

For the protostellar core, the column density distribution is not symmetric: it exhibits a head-tail shape. The centre of gravity determined from this distribution is significantly offset from the position of the protostar itself (see Fig. D.1). This suggests that the streamer does not point directly towards the gravity centre of the protostellar core. If its curved structure is generated by the angular momentum of the infalling gas, a large rotating envelope (~10<sup>3</sup> au) should form around the protostar and connect to the streamers. However, there is currently no observational evidence of such a large structure (e.g. Tobin et al. 2016; Pineda et al. 2020). Therefore, it is unlikely that the streamer is driven by gravitational infall from outside the core.

## 3.2. CCS/HC<sub>3</sub>N bridge with a length of 3 × 10<sup>4</sup> au (0.2 pc)

In Fig. 1 we present the velocity-integrated maps of CCS, HC<sub>3</sub>N, and HC<sub>5</sub>N alongside the *Herschel* column density map. In the CCS integrated intensity map, we recognise the two prominent CCS elongated features: (1) an elongated feature, stretching from north to south and bridging the gap between the starless cores 1 and 2 and the protostellar core 3, and (2) a compact emission positioned to the west of the protostellar core.

The distribution of HC<sub>3</sub>N emission appears to resemble that of CCS, with strong concentrations at the southern periphery of the starless cores 1 and 2. Another significant peak is evident

immediately above the protostar, and a fainter peak is located to the west, coinciding with a local maximum in CCS intensity. HC<sub>5</sub>N, which exhibits the weakest intensity of the three lines, has a peak at a similar position as the second peak of HC<sub>3</sub>N, or the northern portion of the protostellar core.

Figures 2 and C.1 display the velocity channel maps of CCS, HC<sub>5</sub>N, and HC<sub>3</sub>N. The starless cores 1 and 2 appear to be associated with relatively strong compact emission around a velocity of 6.5 km s<sup>-1</sup>.

The CCS emission within the velocity range 6.5–6.7 km s<sup>-1</sup> exhibits a fragmented filamentary structure, extending along a line from the north-east to the south-west that coincides with the large-scale filament seen in the H<sub>2</sub> column density image (see Fig. D.1). An elongated structure is discernible in the south-east corner of the CCS map, though the peak is faint in the H<sub>2</sub> map (see Fig. 1).

At a velocity of 6.6–7.0 km s<sup>-1</sup>, the elongated structure observed in CCS gradually connects with the central protostellar core, forming a bridge. A similar bridging structure is also observed in the corresponding velocity channel for HC<sub>5</sub>N. Overall, these emissions demonstrate a weak velocity gradient that spans from 6.5 km s<sup>-1</sup> to 7.0 km s<sup>-1</sup> towards the protostar.

## 3.3. Mismatch in velocity structure between the bridge and streamer

From the observations detailed in the previous section, it becomes evident that the protostellar core 3 and the northern starless cores 1 and 2 exhibit a connection characterised by a bridge structure. In Fig. 3 we present a CCS position-velocity plot that traces the 0.2 pc bridge structure.

The protostellar core 3 exhibits a velocity of ~7 km s<sup>-1</sup>, indicating a velocity disparity between the two cores bridged by an intermediate-velocity region. The velocity gradient along the bridge is estimated to be ~0.5 km s<sup>-1</sup>/0.2 pc ~ 2.5 km s<sup>-1</sup> pc<sup>-1</sup>. The existence of the intermediate-velocity gas bridging the two structures with different velocities suggests a dynamical interaction between the cores, a feature similar to larger-scale cloud-cloud collisions (e.g. Nakamura et al. 2012; Nguyen-Lu’o’ng et al. 2013; Fukui et al. 2021; Kinoshita et al. 2021). The enhanced CCS emission observed could potentially stem from the colliding region, where CCS formation may be accelerated due to increased density. Both CCS and HC<sub>3</sub>N exhibit critical densities of around 10<sup>4</sup> cm<sup>-3</sup>, implying that while the shocked layer may have a higher density, the column density might not be sufficiently high to be recognised in the 2D *Herschel* image.

The streamer seemingly follows the ~10<sup>22</sup> cm<sup>-2</sup> contour of the protostellar core and spatially overlaps with the distribution
<figure><img src="image_3.png" /></figure>

**Fig. 2.** Colour-scale image showing channel maps of CCS ($J_N = 4_3 - 3_2$). The red contours are the CCS ($J_N = 8_7 - 7_6$) integrated intensity contours from NOEMA obtained by Pineda et al. (2020). The black contours are Herschel column densities derived by Pezzuto et al. (2021), from $1 \times 10^{22}$ to $10 \times 10^{22} \text{ cm}^{-2}$ in steps of $0.2 \times 10^{22} \text{ cm}^{-2}$. The blue ellipse is the map effective beam size.

<figure><img src="image_4.png" /></figure>

**Fig. 3.** Velocity structures in Per-emb-2. Left: Integrated-intensity map of CCS (from 6.0 to 8.0 km s⁻¹) in colour overlaid with the Herschel column density map in contours (as in Fig. 1) and the position of the velocity cut along the bridge (in white). Right: Position-velocity plots along the bridge structure from the cut along the bridge seen on the left. White and red stars indicate the positions of starless cores and the protostellar core. The yellow cross indicates the position of the CCS peak. The blue ellipse is the map effective beam size.

of the CCS emission (see Fig. 2), but the velocities are very different. The streamer has a much steeper velocity gradient in the opposite direction to that of the bridge (see the right panel of Fig. 3). Therefore, we conclude that the bridge is not a structure connected to the streamer.

# 4. Discussion

## 4.1. Accreting streamer or not?

Taniguchi et al. (2024) claim, based on the integrated intensity map, that the streamer connects to the northern region through the bridge. They suggest that the streamer is an accreting flow. However, as mentioned above, the gas in the bridge is likely a separate component since the streamer and bridge have different velocities: the velocity of the bridge ranges from 6.5 (in the north) to 7.0 km s⁻¹ (south) with a velocity gradient of 2.5 km s⁻¹ pc⁻¹ (hereafter, we define this direction as positive), whereas the streamer has velocities of 7.5 (north) to 7.0 km s⁻¹ (south) with a velocity gradient of -0.5 km s⁻¹/0.05 pc = -10 km s⁻¹ pc⁻¹. From the velocity distribution, there is no evidence that the streamer connects to the bridge.

In addition, the centre of gravity of the protostellar core is significantly offset from the protostar, by ~12'' (~7000 au). If the streamer originates from the inflow from the ambient gas driven by the gravity, it likely points towards the core's gravity centre (see Fig. D.1). If the offset were due to rotation, this object would have a large rotating envelope with a size of ~10³ au. Such a structure is not detected towards the protostar (Tobin et al. 2016). This also implies that the streamer is not an origin of accretion, contrasting with the interpretation of Pineda et al. (2020). The northern cores (1 and 2) are close to gravitationally bound and therefore could potentially evolve into protostars. Their merging could trigger star formation. The mass of the protostellar core is significant; a part of it may eventually accrete onto the central protostar before the ambient gas falls in.

Yano et al. (2024) used the DCC scenario to reproduce the density structure in Per-emb-2 in numerical simulations. According to this model, when two cores collide with each other, a shock-compressed layer is formed in between them and connects to the protostar. The gas is inflowing towards the protostar. However, if this is the case, both the bridge and streamer should have velocity gradients in the same direction (see Kinoshita & Nakamura 2022). The velocity and density
<figure><img src="image_5.png" /><figcaption>**Fig. 4.** CFC model of Per-emb-2. (1) A dense core is approaching a filament (grey) with a speed of ~1 km s⁻¹. (2) They collide, protostar formation is triggered at the shock-compressed part, and a streamer (dark blue) of size 3 × 10⁴ au is formed. (3) The dense core shreds the filament into two unequal parts, and cores 1, 2, and 5 are created. The shock-compressed parts of the filament become a bridge (yellow), where the formation of carbon-chain molecules is accelerated by the increased density.</figcaption></figure>

structures appear more complicated than what the DCC model predicts. Therefore, the formation mechanism of the streamer is neither ambient gas inflow nor two-round core collision.

## 4.2. A core-filament collision (CFC) model

Based on the evidence of the existence of the large-scale filament and a positive gradient towards the streamer of the protostar, we propose a CFC model to explain the configuration of Per-emb-2. Figure 4 depicts the schematic representation of this CFC model.

We propose that: (1) A spherical core, with a relative speed of ~1–1.5 km s⁻¹, approached from the west side of the filament and collided with it. (2) The collision initiated star formation on the far side of the eastern extremity of the core, where the initial shock compression occurred, leading to the formation of a protostellar system. The arc-like streamer may be the bow structure formed by the CFC, and therefore the velocity gradient of the streamer is not due to gravitational accretion. If its dynamics are not significantly affected by the stellar gravity, the flow along the streamer should be pointed towards the outside (i.e. it is an outflow and not infalling). The bow typically has an almost linear velocity structure because of its curvature. Such a feature is consistent with the observation. When the stellar gravity is sufficiently strong to attract the gas towards the star, the flow may be pointing towards the protostar. Since it is difficult to judge the direction (outflow or inflow) from the current observational data, numerical experiments are needed to verify this scenario more quantitatively.

The head-tail morphology observed in core 3 agrees well with what would result from such a collision. The CCS arc (streamer) structure (Pineda et al. 2020) appears to trace the edge of the protostellar core (~1–2 × 10²² cm⁻²; see Fig. 2).

Furthermore, the collision led to the division of the filament into two distinct structures (cores 1 and 2 and core 5). The shock-compressed layers in the northern and southern regions attained densities of ~10⁴ cm⁻³, resulting in relatively strong CCS and HC₃N emissions. The northern shock-compressed layer became interconnected with the central protostellar system. There is a possibility that it could evolve into an accreting streamer in the future. However, the total accretion mass via the bridge (future streamer) may be low. The protostellar core envelope still contains enough material (~10 M☉). Therefore, the accretion along the streamers, if it exits, plays only a minor role in the mass accumulation towards the central protostar.

The total velocity difference is about 1 (= 7.5–6.5) km s⁻¹, and the filament width is about 0.1 pc; therefore, it takes

0.1 pc/1 km s⁻¹ ~ 10⁵ yr for the core to cross the filament. This is comparable to the typical lifetime of a Class 0 protostar and is consistent with the triggered star formation proposed in the CFC model scenario. Furthermore, if the bridge is the remnant of the shock-compressed layer, its density must have increased from $n_{amb} \sim 10^3$ cm⁻³ to ~10⁴ cm⁻³ due to the isothermal shock with a Mach number of 1.5 km s⁻¹/sound speed ~6. This condition in the shocked area, or streamer, is similar to the condition assumed by Suzuki et al. (1992), who started the chemical evolution calculation assuming a steady state with a uniform isothermal gas of 10⁴ cm⁻³ and T = 10 K. According to Fig. 13 of Suzuki et al. (1992), CCS and HC₃N become most abundant in ~10⁵ yr in such a density range, and declines steeply after 0.5 Myr (see also Leung et al. 1984). These timescales match the CFC scenario. In our CFC or DCC scenarios, shocks are likely to increase the gas temperature somewhat, but the molecular gas is highly radiative and therefore the gas temperature stays low even after shock compression. In such a situation, the streamers formed by the collision model may have emission of SO or CH₃OH, whose sublimation temperatures are about 50 and 80 K, respectively. These lines could exhibit broader widths as a result of shock interactions. This hypothesis should be validated with chemical evolution calculations using, for example, the Paris-Durham shock code (Godard 2019).

*Acknowledgements.* This work was financially supported by JSPS KAKENHI Grant Numbers JP23H01218 (F.N.). Part of this work was supported by the NAOJ Visiting Research Grant (Q.N.L.). We thank the NRO staff for both operating the 45 m and helping us with the data reduction. We thank the anonymous referee for valuable comments that helped improve the paper.

## References

André, P., Men'shchikov, A., Bontemps, S., et al. 2010, A&A, 518, L102

Bertoldi, F., & McKee, C. F. 1992, ApJ, 395, 140

Bonnell, I. A., & Bate, M. R. 2006, MNRAS, 370, 488

Chen, X., Ren, Z.-Y., Li, D.-L., et al. 2021, ApJ, 923, L20

Fukui, Y., Habe, A., Inoue, T., Enokiya, R., & Tachihara, K. 2021, PASJ, 73, S1

Godard, B., Pineau des Forêts, G., Lesaffre, P., et al. 2019, A&A, 622, A100

Hsieh, T. H., Segura-Cox, D. M., Pineda, J. E., et al. 2023, A&A, 669, A137

Kinoshita, S. W., & Nakamura, F. 2022, ApJ, 937, 69

Kinoshita, S. W., Nakamura, F., Nguyen-Luong, Q., et al. 2021, PASJ, 73, S300

Kuffmeier, M., Haugbølle, T., & Nordlund, Å. 2017, ApJ, 846, 7

Leung, C. M., Herbst, E., & Huebner, W. F. 1984, ApJS, 56, 231

Maruta, H., Nakamura, F., Nishi, R., Ikeda, N., & Kitamura, Y. 2010, ApJ, 714, 680

Nakamura, F., & Li, Z.-Y. 2007, ApJ, 662, 395

Nakamura, F., Miura, T., Kitamura, Y., et al. 2012, ApJ, 746, 25

Nakamura, F., Sugitani, K., Tanaka, T., et al. 2014, ApJ, 791, L23

Nakamura, F., Chiong, C.-C., Taniguchi, K., et al. 2024, PASJ, 76, 563

Nguyen-Lu'o'ng, Q., Motte, F., Carlhoff, P., et al. 2013, ApJ, 775, 88

Olguin, F. A., Sanhueza, P., Chen, H.-R. V., et al. 2023, ApJ, 959, L31

Padoan, P., Pan, L., Juvela, M., Haugbølle, T., & Nordlund, Å. 2020, ApJ, 900, 82

Pezzuto, S., Benedettini, M., Di Francesco, J., et al. 2021, A&A, 645, A55

Pineda, J. E., Segura-Cox, D., Caselli, P., et al. 2020, Nat. Astron., 4, 1158

Rosolowsky, E. W., Pineda, J. E., Kauffmann, J., & Goodman, A. A. 2008, ApJ, 679, 1338

Sanhueza, P., Contreras, Y., Wu, B., et al. 2019, ApJ, 886, 102

Sanhueza, P., Girart, J. M., Padovani, M., et al. 2021, ApJ, 915, L10

Shu, F. H., Adams, F. C., & Lizano, S. 1987, ARA&A, 25, 23

Suzuki, H., Yamamoto, S., Ohishi, M., et al. 1992, ApJ, 392, 551

Taniguchi, K., Pineda, J. E., Caselli, P., et al. 2024, ApJ, 965, 162

Tobin, J. J., Looney, L. W., Li, Z.-Y., et al. 2016, ApJ, 818, 73

Tu, Y., Li, Z.-Y., Lam, K. H., Tomida, K., & Hsu, C.-Y. 2024, MNRAS, 527, 10131

Valdivia-Mena, M. T., Pineda, J. E., Segura-Cox, D. M., et al. 2022, A&A, 667, A12

Yano, Y., Nakamura, F., & Kinoshita, S. W. 2024, ApJ, 964, 119

Zucker, C., Speagle, J. S., Schlafly, E. F., et al. 2019, ApJ, 879, 125
# Appendix A: Observations and target lines

The on-the-fly observations were carried out in October 2022 (2 nights) and in March 2023 (1 daytime). Z45 is dual linear polarisation receiver and we summed up the horizontal and vertical polarisation components to improve the signal-to-noise ratios. The scan interval of the on-the-fly observations was set to 8'', about a fifth of the beam size. The pointing observations were made every one hour with SiO maser lines of NML-Tau (KL Tau). The pointing errors were within 3''–5''. We used SAM45 digital spectrometer as a backend. SAM45 is a highly flexible FX-type digital spectrometer with a finest frequency resolution of 3.81 kHz. The bandwidth of SAM45 with 3.81 kHz was 15.625 MHz. It took about 30 min to obtain a single map. We summed 6 maps to obtain a final map. The total observation time was about 6 hours including the overhead of pointing observations. The on-source time was 3 hours. During the observations, the typical system noise temperatures were about 150–200 K.

The rest frequencies and transitions of the target lines are summarised in Table A.1 (see Nakamura et al. 2024 for the rest frequency of CCS).

<table><caption>**Table A.1.** Target lines.</caption><thead><tr><th>Species</th><th>Transition</th><th>Frequency (GHz)</th></tr></thead><tbody><tr><td>CCS</td><td>J<sub>N</sub> = 4<sub>3</sub> – 3<sub>2</sub></td><td>45.379033</td></tr><tr><td>HC<sub>3</sub>N</td><td>J = 5 – 4</td><td>45.490316</td></tr><tr><td>HC<sub>5</sub>N</td><td>J = 17 – 16</td><td>45.264720</td></tr></tbody></table>

# Appendix B: Core identification and physical quantities

## B.1. Core identification

We applied *astrodendro* (Rosolowsky et al. 2008) to the high-resolution *Herschel* column density map with an angular resolution of 18''. The *astrodendro* has three parameters for the hierarchical structure analysis: *min_value* which is the threshold value, *min_delta* which is the minimum step for structure identification, and *min_npix* which is the minimum pixel number. We adopt *min_value* = 3σ, *min_delta* = 1.5σ, and *min_npix* = beam area/pixel size, σ (= 6.81 × 10²⁰ cm⁻²) is the rms noise level. The *astrodendro* identifies three structures: a leaf, branch, and trunk. The minimum structure of ‘a leaf’ is defined as a dense core. To estimate the ambient gas pressure around the cores, we adopted the physical quantities averaged in a branch immediately below the corresponding leaf. We also derive the core mass by subtracting the background component which is defined as the trunk.

## B.2. Dynamical states of the cores

We assess the dynamical states of the cores applying the virial theorem. The virial equation for a spherical core is given by

$$ \frac{1}{2} \frac{\partial^2 I}{\partial t^2} = U + W + S, \qquad (B.1) $$

where *I*, *U*, *W*, and *S* are the moment of inertia, internal kinetic energy, gravitational energy, and surface pressure term, respectively, and the magnetic field terms are neglected (Maruta et al.

2010). These individual terms are given as

$$ U = 3M_{\text{core}}\sigma_{\text{tot}}^2 \qquad (B.2) $$

$$ W = -\frac{3}{5}a \frac{GM_{\text{core}}^2}{R_{\text{core}}} \qquad (B.3) $$

$$ S = -4\pi R_{\text{core}}^3 P_{\text{ex}}, \qquad (B.4) $$

where *G*, *M*<sub>core</sub>, *R*<sub>core</sub>, and *P*<sub>ex</sub> are the gravitational constant, core mass, core radius, and the external pressure, respectively. *a* is a dimensionless parameter of order unity which measures the effects of a non-uniform or nonspherical mass distribution (Bertoldi & McKee 1992). For a uniform sphere and a centrally condensed sphere with ρ ∝ r⁻², *a* = 1 and 5/3, respectively. Here, we adopt *a* = 5/3 since the cores tend to be more or less centrally-condensed. σ<sub>tot</sub> is the 1D intrinsic velocity dispersion of the molecule of mean mass,

$$ \sigma_{\text{tot}} = \sqrt{\sigma_{\text{CCS}}^2 + k_B T \left( \frac{1}{\mu m_H} - \frac{1}{m_{\text{CCS}}} \right)}, \qquad (B.5) $$

and *G* is the gravitational constant, *k*<sub>B</sub> is the Boltzmann constant, σ<sub>CCS</sub> is the velocity dispersion measured by CCS, μ (=2.33) is the mean molecular weight, *m*<sub>H</sub> is the mass of a hydrogen atom, and *m*<sub>CCS</sub> is the mass of a CCS molecule. The CCS velocity dispersions towards the cores are computed by the Gaussian fitting of the averaged spectra (see Fig. B.1). It is worth noting that CCS emission may preferentially trace the outer envelopes of cores, given its critical density of 10⁴ – 10⁵ cm⁻³. In contrast, the core densities calculated from the *Herschel* image range from 3 × 10⁶ to 2 × 10⁴ cm⁻³ as shown in Table 1. Despite this, in the core analysis presented here, we adopt the CCS line widths as representative of the cores’ properties.

The external pressure is calculated as

$$ P_{\text{ex}} = 3\rho_{\text{amb}}\sigma_{\text{tot,amb}}^2, \qquad (B.6) $$

where ρ<sub>amb</sub> and σ<sub>tot,amb</sub> are the density and 1D intrinsic velocity dispersion of ambient gas, respectively, and they are derived from the values averaged in the branch immediately below the corresponding leaf. The ambient density was computed under the assumption that the branch is spherical symmetry. The core virial parameter is calculated as

$$ \alpha_{\text{vir}} = \frac{U}{W} = \frac{3\sigma_{\text{tot}}^2 R_{\text{core}}}{GM_{\text{core}}}. \qquad (B.7) $$

Table C.1 summarises the velocity dispersions of the cores and their energy ratios. Figure B.2 shows the energy ratios *W*/*U* versus *S*/*U*, where *W*, *U*, and *S* are the gravity term, internal kinetic and thermal energy term, and the surface term, respectively. In this plot, the cores below the black dashed line are satisfied with the condition *S* + *W* + *U* < 0, for which the cores collapse. The core’s self-gravity is important for the area above the dotted line (|*S*/*U*| < |*W*/*U*|), whereas the surface pressure is dominant in the area below the dotted line. The cores 1, 2, and 3 are located in the gravity-dominant area. They are close to the dynamical equilibrium and self-gravitating. For the core 5, the external pressure plays a crucial role in the core’s dynamical stability, and thus the core is dynamically unstable due to the external pressure. Such a core is difficult to create from the gravitational fragmentation of the filament. The formation of unstable non-self-gravitating core (core 5) appears to agree with the collision scenario.
<figure><img src="image_6.png" /><figcaption>**Fig. B.1.** CCS spectra averaged over the entire mapping field (**top left**), at core 1 (**top middle**), at core 2 (**top right**), at core 3 (**bottom left**), at core 4 (**bottom middle**), and at core 5 (**bottom right**).</figcaption></figure>

<figure><img src="image_7.png" /><figcaption>**Fig. B.2.** Dynamical stability of the cores based on the CCS velocity dispersion. The black dashed line represents $S + W + U = 0$.</figcaption></figure>

## Appendix D: Centre of gravity of the protostellar core with respect to the large filament

The circle in Fig. D.1 indicates the position of the protostar. The black square indicates the column density weighted position of the protostellar core. If the gas were distributed only on the plane of sky, the position would coincide with the centre of gravity. In other words, the protostar was not formed at the core centre but near the core envelope. This implies that the core formation did not occur spontaneously.

If the curved shape of the streamer identified by Pineda et al. (2020) is due to the angular momentum of infalling gas, the rotational motion must be significant on the scale of $10^3$ au, leading to the formation of a large (~ $10^3$ au) rotating envelope or disk. However, such a large circumstellar structure was not detected in the observations by Tobin et al. (2016). Thus, at this scale, the rotational motion would still be minor and unlikely to create the elongated, curved structure seen in Fig. 2. For example, Yano et al. (2024)'s core collision simulation showed that the streamer formed by the infalling gas from the second core is smoothly connected to a central, rotating disk-like envelope. Such a structure is a plausible outcome of the infalling gas scenario. Since the rotating disk is small for this protostar, the curved structure, if formed by the infalling gas, must almost point directly to the protostar. In Fig. D.1 we have also indicated the position of the filament shown in Fig. 4.

## Appendix C: Channel maps

Figure C.1 shows the channel maps of HC₃N ($J = 5-4$) and HC₅N ($J = 17-16$). The HC₃N has a stronger emission but having several overlapped hyperfine components. However, adopting the rest frequency of $J = 5-4$, $F = 5-4$ as a centre frequency, we simply integrated all the emission to construct the channel maps. The bridge is prominent for HC₃N ($J = 5-4$), and structure similar to CCS was detected. For HC₅N ($J = 17-16$), the strong emission is seen just at the southern edge of cores 1 and 2.
<figure><img src="image_8.png" /><figcaption>**Fig. C.1.** Colour-scale image showing channel maps of HC₃N (J = 5 – 4) (upper) and HC₅N (J = 17 – 16) (lower). The red contour shows the CCS (Jₙ = 8₇ – 7₆) integrated intensity contours from NOEMA obtained by Pineda et al. (2020). The black contours show Herschel column densities derived by Pezzuto et al. (2021), from 1×10²² to 10×10²² cm⁻² in steps of 0.2×10²² cm⁻². Note that the HC₃N has overlapping hyperfine components and that we made the cube data by adopting the rest frequency indicated in Table A.1. The blue ellipse is the map effective beam size.</figcaption></figure>

<table><caption>**Table C.1.** Dynamical properties of cores detected in the mapping area by astrodendro.</caption><thead><tr><th>Number</th><th>σ<sub>CCS</sub><br>(km s<sup>-1</sup>)</th><th>σ<sub>CCS,amb</sub><br>(km s<sup>-1</sup>)</th><th>σ<sub>tot</sub><br>(km s<sup>-1</sup>)</th><th>ρ<sub>amb</sub>/(μm<sub>H</sub>)<br>(cm<sup>-3</sup>)</th><th>σ<sub>tot,amb</sub><br>(km s<sup>-1</sup>)</th><th>W/U</th><th>S/U</th><th>α<sub>vir</sub></th></tr></thead><tbody><tr><td>1</td><td>0.37</td><td>0.41</td><td>0.26</td><td>4.53 × 10<sup>4</sup></td><td>0.28</td><td>-0.703</td><td>-0.111</td><td>1.422</td></tr><tr><td>2</td><td>0.41</td><td>0.41</td><td>0.27</td><td>4.53 × 10<sup>4</sup></td><td>0.27</td><td>-0.986</td><td>-0.019</td><td>1.014</td></tr><tr><td>3</td><td>0.84</td><td>0.61</td><td>0.41</td><td>3.02 × 10<sup>4</sup></td><td>0.33</td><td>-0.928</td><td>-0.734</td><td>1.077</td></tr><tr><td>5</td><td>0.45</td><td>0.61</td><td>0.28</td><td>2.49 × 10<sup>4</sup></td><td>0.34</td><td>-0.339</td><td>-3.025</td><td>2.954</td></tr></tbody></table>
<figure><img src="image_9.png" /><figcaption>**Fig. D.1.** *Herschel* H₂ column density image. The white rectangle indicates the filament that collided with the dense core in our CFC model. The white box with the coordinates is the observation box. The circle at the bottom-left corner represents the beam size of the *Herschel* map. The small circle near the centre is the position of the protostar. The colour represents the H₂ column density in cm⁻². The green contours are drawn at the H₂ column density levels of 5 × 10²¹, 2 × 10²², 3 × 10²², 5 × 10²², and 1 × 10²³ cm⁻³.</figcaption></figure>