Issue |
A&A
Volume 638, June 2020
|
|
---|---|---|
Article Number | A7 | |
Number of page(s) | 8 | |
Section | Interstellar and circumstellar matter | |
DOI | https://doi.org/10.1051/0004-6361/202037815 | |
Published online | 29 May 2020 |
The role of Galactic H II regions in the formation of filaments
High-resolution submilimeter imaging of RCW 120 with ArTéMiS★
1
Aix Marseille Univ, CNRS, CNES, LAM,
Marseille,
France
e-mail: annie.zavagno@lam.fr
2
Institut Universitaire de France,
Paris,
France
3
Laboratoire d’Astrophysique (AIM), CEA/DRF, CNRS, Université Paris-Saclay, Université Paris Diderot,
Sorbonne Paris-Cité,
91191
Gif-sur-Yvette, France
4
Leibniz-Institut für Astrophysik Potsdam (AIP),
An der Sternwarte 16,
14482
Potsdam,
Germany
5
School of Physics & Astronomy, Cardiff University,
Cardiff
CF24 3AA,
UK
6
Department of Physics and Astronomy, Graduate School of Science and Engineering, Kagoshima University,
1-21-35 Korimoto, Kagoshima,
Kagoshima
890-0065,
Japan
7
Instituto de Astrofísica e Ciências do Espaço, Universidade do Porto, CAUP,
Rua das Estrelas,
PT4150-762
Porto,
Portugal
8
Laboratoire d’Astrophysique de Bordeaux, Univ. Bordeaux, CNRS, B18N,
allée Geoffroy Saint-Hilaire,
33615
Pessac, France
9
National Centre for Nuclear Research,
Ul. Pasteura 7,
02-093
Warsaw, Poland
10
Jodrell Bank Centre for Astrophysics, Department of Physics and Astronomy, The University of Manchester,
Oxford Road,
Manchester
M13 9PL,
UK
11
Jeremiah Horrocks Institute, University of Central Lancashire,
Preston PR1 2HE, UK
12
Institut d’Astrophysique de Paris, Sorbonne Université, CNRS (UMR7095),
75014
Paris, France
13
I. Physik. Institut, University of Cologne,
50937
Cologne,
Germany
Received:
25
February
2020
Accepted:
6
April
2020
Context. Massive stars and their associated ionized (H II) regions could play a key role in the formation and evolution of filaments that host star formation. However, the properties of filaments that interact with H II regions are still poorly known.
Aims. To investigate the impact of H II regions on the formation of filaments, we imaged the Galactic H II region RCW 120 and its surroundings where active star formation takes place and where the role of ionization feedback on the star formation process has already been studied.
Methods. We used the large-format bolometer camera ArTéMiS on the APEX telescope and combined the high-resolution ArTéMiS data at 350 and 450 μm with Herschel-SPIRE/HOBYS data at 350 and 500 μm to ensure good sensitivity to a broad range of spatial scales. This allowed us to study the dense gas distribution around RCW 120 with a resolution of 8′′ or 0.05 pc at a distance of 1.34 kpc.
Results. Our study allows us to trace the median radial intensity profile of the dense shell of RCW 120. This profile is asymmetric, indicating a clear compression from the H II region on the inner part of the shell. The profile is observed to be similarly asymmetric on both lateral sides of the shell, indicating a homogeneous compression over the surface. On the contrary, the profile analysis of a radial filament associated with the shell, but located outside of it, reveals a symmetric profile, suggesting that the compression from the ionized region is limited to the dense shell. The mean intensity profile of the internal part of the shell is well fitted by a Plummer-like profile with a deconvolved Gaussian full width at half maximum of 0.09 pc, as observed for filaments in low-mass star-forming regions.
Conclusions. Using ArTéMiS data combined with Herschel-SPIRE data, we found evidence for compression from the inner part of the RCW 120 ionized region on the surrounding dense shell. This compression is accompanied with a significant (factor 5) increase of the local column density. This study suggests that compression exerted by H II regions may play a key role in the formation of filaments and may further act on their hosted star formation. ArTéMiS data also suggest that RCW 120 might be a 3D ring, rather than a spherical structure.
Key words: stars: formation / H II regions / ISM: individual objects: RCW 120 / submillimeter: ISM
© A. Zavagno et al. 2020
Open Access article, published by EDP Sciences, under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/4.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
1 Introduction
Although the existence of filaments has been known for a long time (Schneider & Elmegreen 1979; Ungerechts & Thaddeus 1987), the tight relation between filaments and star formation has only been clearly revealed by Herschel observations, which introduce a new paradigm for Galactic star formation (André et al. 2014; Molinari et al. 2014). Filaments are ubiquitous in the Galactic interstellar medium (ISM) and their properties have been extensively studied and discussed (Arzoumanian et al. 2011; Palmeirim et al. 2013; Schisano et al. 2014, 2020; Cox et al. 2016; Li et al. 2016). One of their most striking properties is their typical inner width of 0.1 pc that seems to be quasi-universal (André et al. 2016; Arzoumanian et al. 2019) and which could result from their formation process (Federrath 2016; André 2017). However, this point is still debated and high-resolution observations are needed to solve the question of filament width universality (Henshaw et al. 2017; Ossenkopf-Okada & Stepanov 2019; Xu et al. 2019). The study of filaments has attracted much attention being the original matter and host zone of star formation. In particular, the mass per unit length of filaments, which could determine the mass of the cores that will be formed (André et al. 2019) together with the filaments’ ability to fuel the surrounding material toward dense cores allowing for a continuous accretion (Gómez et al. 2018), can change our view of the star formation process. Understanding filaments as a whole (formation, evolution) has become a new quest of star formation studies (e.g., Zhang et al. 2019).
Another debated question is the role of expanding H I shell and ionized (H II) regions in the formation of filaments and their influence on star formation. Many observational studies have shown that high-mass star formation is observed at the edges of Galactic (Deharveng et al. 2005, 2010; Thompson et al. 2012) and extragalactic (Bernard et al. 2016) H II regions. This phenomenon is important because at least 30% of the Galactic high-mass star formation is observed at the edges of H II regions (Palmeirim et al. 2017), suggesting that these regions could createthe needed conditions for the formation of high-mass stars and might even increase its efficiency (Motte et al. 2018). Models show the complexity of H II regions’ impact on the surrounding medium (Geen et al. 2017, 2020). However, observational studies demonstrate the importance of ionization feedback in assembling denser material in their photodissociation region (PDR), leading to favorable conditions there for the formation of a new generation of (high-mass) stars (Tremblin et al. 2014; Orkisz et al. 2017). H II regions have also been proposed as a key agent of the “bubble-dominated” scenario for the formation of filamentary molecular clouds proposed by Inutsuka et al. (2015): Previous generations of massive stars generate expanding shells and bubbles; additionally, dense molecular filaments result from multiple large-scale supersonic compressions at the surface of such expanding bubbles and/or at the interface between two bubbles. Three-dimensional magnetohydrodynamic simulations of interacting supershells show filamentary structures at the surface of the shells (Ntormousi et al. 2017; Dawson et al. 2015). Recent modeling and observations of PDR in high-mass star-forming regions (Bron et al. 2018; Goicoechea et al. 2019a,b) also show the importance of understanding the role of ionization feedback on star formation and the role of H II regions in the formation of filaments. However, the link between ionization feedback, the possible filamentary structure of PDRs, and star formation observed at the edges of H II regions is still poorly known. Addressing this question requires studying the dense PDRs that surround H II regions and their role in star formation.
RCW 120 is a galactic H II region located above the galactic plane at l = 348.24°, b = 00.46° (RA (J2000) = 258.10°, Dec (J2000) = − 38.46°). The region is ionized by a single O6-8V/III star (Martins et al. 2010). Due to its simple ovoid morphology and nearby location, 1.34 kpc (Deharveng et al. 2009), it has been extensively studied in the recent past, especially in order to discuss the properties of star formation observed at its edges (Deharveng et al. 2009; Zavagno et al. 2010; Figueira et al. 2017, 2018) and also to propose a scenario for its original formation by a cloud–cloud collision (Torii et al. 2015) as well as its kinematics (Anderson et al. 2015; Sánchez-Cruces et al. 2018; Kohno et al. 2018; Kirsanova et al. 2019). The wealth of existing data on this region makes it ideal to study the impact of ionization on the structure of its dense PDR and its relation with the observed young star formation.
With the aim of characterizing the properties of dense PDRs around H II regions and their relation with the star formation observed at their edges, we present the results of dust continuum mapping observations of the RCW 120 region with the ArTéMiS bolometer camera on the APEX 12 m telescope. The 8′′ resolution of ArTéMiS at 350 μm, corresponding to 0.05 pc at the distance of RCW 120, allowed us to reveal the structure of the dense PDR, traced here by a filament, that surrounds the ionized region and to resolve, for the first time, its width. Section 2 presents the observations and data reduction. Section 3 presents the mapping results, which are then discussed in Sect. 4. Conclusions are summarized in Sect. 5.
2 ArTéMiS observations and data reduction
Observations of RCW 120 were obtained at 350 and 450 μm on August 27 and October 10 in 2017 and on May 6, 21, 22, and June 2 in 2018 with the ArTéMiS camera on the Atacama Pathfinder Experiment (APEX; APEX Service Mode run 0100.F-9305(A)). The total on-source observing time is 12.5 h. ArTéMiS1 is a large-format bolometer array camera, which was built by CEA/Saclay and installed in the Cassegrain cabin of APEX, which has a total of 4608 pixels observing at 350 and 450 μm simultaneously. A complete description of the instrument is given in André et al. (2016). We obtained the 350 and 450 μm maps by coadding the 31 individual maps, obtained at each wavelength, on RCW 120 using the total-power on-the-fly scanning mode with a scanning speed of 30′′ /s. The absolute calibration uncertainty is about 30%. The data reduction procedure is the same as the one described in André et al. (2016). The root mean square (rms) noise is σ ≃ 0.25 Jy/8′′-beam at 350 μm and 0.21 Jy/10′′-beam at 450 μm.
3 Results
3.1 Combination of ArTéMiS and Herschel data
The ArTéMiS raw data are affected by a high level of sky noise, which is strongly correlated over the detectors of the focal plane. Subtracting this correlated sky noise during data reduction makes the resulting image insensitive to angular scales larger than the instantaneous field of view of the camera (≃ 2′× 4′). To restore the missing large-scale information, we combined the ArTéMiS data with the SPIRE 350 and 500 μm data from the Herschel HOBYS keyproject (Motte et al. 2010), which employed a similar technique as the one used in combining millimeter interferometer observationswith single-dish data. The task “immerge” in the Miriad software package is used (Sault et al. 1995). A complete description of the method is given in André et al. (2016). Figure 1 shows the distribution of the ArTéMiS-only 350 μm (left) and 450 μm (right) emissions observed with ArTéMiS. These emissions show the main bright condensations identified previously using continuum millimeter observations (Zavagno et al. 2007; Deharveng et al. 2009) that host star formation, which were studied in detail with Herschel (Zavagno et al. 2010; Figueira et al. 2017) and ALMA (Figueira et al. 2018). The northern part of the dense PDR layer surrounding the H II region is seen with fainter emission that delineates the ovoid shape of the region.
![]() |
Fig. 1 ArTéMiS-only maps of RCW 120 at 350 μm (left) and 450 μm (right) emission observed with ArtéMiS. The white circle shows the location of the ionizing star (Right ascension (J2000) = 258.085833°, Declination (J2000) = −38.490555°). North is up,and east is left. |
3.2 Filaments extraction and distribution
We applied the DisPerSe (Sousbie 2011) and getfilaments (Men’shchikov 2013) algorithms to determine the location of the crest of filaments in RCW 120. Both algorithms were applied to the combined Herschel + ArTéMiS 350 μm image. The main crests tracing the dense PDR are shown as the white solid curves in Fig. 2 together with the young sources identified with Herschel and discussed in Figueira et al. (2017). The sources are distributed along the crests with a higher density of sources in the densest part of the shell, which are observed to the south. Sources are also observed to the north in tight relation with filaments’ crests. In this paper, we discuss the filaments’ profile over the dense shell. The properties of the filaments and their relation with the observed star formation will be discussed in a forthcoming paper.
3.3 Radial intensity profile of the shell
Here, we analyze the profile on the dense shell that surrounds RCW 120’s ionized region and that corresponds to the densest part of RCW 120’s PDR. By taking perpendicular cuts at each pixel along the crest, we constructed a median radial intensity profile. Figure 3 shows the shape of the radial intensity profile toward the internal and external part of the dense shell.
We fit all of the observed radial profiles presented in this paper, I(r) as a function of radius r with both a simple Gaussian model,
(1)
and a Plummer-like model function of the form,
(2)
where I0 is the central peak intensity, FWHM is the physical full width at half maximum of the Gaussian model, RFlat represents the characteristic radius of the flat inner part of the Plummer model profile, p > 1 is the power-law index of the underlying density profile, and IBg is a constant background level. These forms were convolved with the approximately Gaussian beam of the ArTéMiS data (FWHM ≃ 8′′) before we compared them with the observed profile. We assume a dust temperature of Td = 15 K (median dust temperature over the region) and a dust opacity law of κλ = 0.1 × (λ∕300 μm)−β cm2 per g (of gas and dust) with an emissivity index of β = 2 (Hildebrand 1983), which is the same opacity law as the one adopted in the HGBS and HOBYS papers (Roy et al. 2014). Figure 4 shows the linear profile of a median cut through the dense shell; the positive offset is toward the shell center. This profile shows clear asymmetry, suggestive of a compression from the internal part of the shell, which is likely of a radiative origin. Indeed, Martins et al. (2010) show that the ionizing star of RCW 120 has a low stellar wind activity. This means that the ionized bubble has a radiative origin and is not created by stellar winds. For the internal part of the shell (compressed part), the mean intensity profile is well fitted by a classic Plummer-like profile (Casali 1986; Whitworth & Ward-Thompson 2001; Arzoumanian et al. 2011; Palmeirim et al. 2013) with p ≃ 2 and DFlat ≃ 0.08 ± 0.02 pc, which corresponds to a Gaussian FWHM of ≃ 0.09 ± 0.02 pc. Toward the external part of the shell, the Gaussian FWHM is ≃ 0.1−0.2 pc, but the profile is much flatter. The Plummer profile is better fitted with p ≃ 1.
3.4 Radial profile of the eastern and western lateral sides of the shell
In order to study the variation compression over the dense shell, we studied the intensity profiles of the two lateral sides of the dense shell, shown as thick red lines in Fig. 2. These eastern and western parts of the resulting median radial intensity profile are displayed in linear format in Fig. 5. These profiles show that the two lateral sides of the shell are also asymmetric, as observed for the dense shell. This result is important because it indicates a homogeneous compression over the dense shell.
![]() |
Fig. 2 Combined ArTéMiS+Herschel-SPIRE 350 μm map of RCW 120, with DisPerSe skeleton overlaid in white and Herschel young sources (solid black circles). The white circle shows the location of the ionizing star. The red parts underline the eastern and western lateral side of the dense shell together with the purple radial filament, which have also been analyzed as separate parts (see text). |
![]() |
Fig. 3 Median radial intensity profile (black curve) of the RCW 120 ’s shell measured perpendicular to the crest, toward the shell’s center (left) and toward the external part of the shell (right). The green error bars show the standard deviation of the mean intensity profile (same color coding for the other figures). The blue solid curve shows the effective beam profile of the ArTéMiS 350 μm data as measured on Mars. The blue dotted curve shows the best-fit Gaussian model to the observed profiles. The red dashed curve shows the best-fit Plummer model convolved with the beam. |
3.5 Profile of part of a radial filament
In order to study the properties of a radial filament that are connected to the shell but situated outside of it, we analyzed the profile of part of a radial filament of RCW 120. We selected this part of the filament, which is located in the southwest and shown in purple in Fig. 2, because it is bright enough to allow for this analysis. By taking perpendicular cuts at each pixel, we constructed two median radial intensity profiles. Figure 6 shows the shape of the radial intensity profile toward the western and eastern parts of this filament. The profiles are noisier than the ones derived for the dense shell. The secondary emission peaks observed along the radial profile correspond to emissions observed along the cut, which are out of the filament. Figure 7 shows the linear profile of a mean cut through the radial filament; the positive offset is to the east. This profile clearly shows a symmetric shape, which is suggestive of no compression from the surroundings, contrary to what is observed for the shell (see Fig. 4).
![]() |
Fig. 4 Observed radial column density (and flux density) profile (black curve) of the RCW 120 shell. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. The yellow error bars show the ± 1σ dispersion of the distribution of radial intensity profile observed along the filament crest (same color coding for the other figures). |
4 Discussion
We are interested in the role of ionized regions in the formation of filaments. If H II regions impact the formation of filaments, they also probably modify the properties of the sources the filaments host, that is, the cores and the stars that form inside. Therefore, it is important to study the relation between H II regions and the filaments they impact.
4.1 Compression by the H II region
Using the ArTéMiS data presented here combined with Herschel data, we revealed the structure of the dense shell observed toward RCW 120. The filament tracing the dense shell shows a clear asymmetric profile, which is suggestive of a radiative compression from the internal part of the ionized region (see Fig. 4). The compressionappears to be homogeneous over the dense shell, as shown in Fig. 5, suggesting that the observed shape of the dense shell that surrounds RCW 120 results from the expansion of the ionized region (see Sect. 4.3). The radiative impact of H II regions on their surroundings has been studied observationally and theoretically. Using Herschel data, Tremblin et al. (2014) studied the impact of the ionization compression on dense gas distribution and star formation in four Galactic ionized regions, including RCW 120. In fitting the probability distribution function (PDF) of the column density of cold dust located around the ionized gas, they show that the PDF of all cloudsis fitted with two lognormal distributions and a power law tail for high column density parts. For RCW 120, they show that the condensations that host star formation and are located at the edge of the ionized gas have a steep compressed radial profile, which is sometimes recognizable in the flattening of the power-law tail (Schneider et al. 2015). This feature has been suggested to be used as an unambiguous criterion to disentangle triggered from pre-existing star formation. Their analysis of RCW 120 shows that the most massive condensation observed at the southwestern edge (condensation 1, see Zavagno et al. 2010; Figueira et al. 2017) may indicate the role of ionization compression for its formation and its collapse. The ArTéMiS data presented here and the analysis of the dense shell’s profile show that compression occurs at the edges of RCW 120. Asymmetry in the filaments’ profile in the Pipe Nebula has been reported by Peretto et al. (2012) who used Herschel images to suggest the formation of filaments by large-scale compression, which likely originate from the winds of the Sco OB2 association. They revealed a bow-like edge for the filamentary structures, suggesting the influence of the compression flow. In RCW 120, we also note the convexity of filaments, with the curvature oriented toward the ionized radiation, which is clearly seen in Fig. 2 and particularly marked in the southwest where high-mass star formation occurs (Figueira et al. 2018). An asymmetry in the filaments’ profile has also been reported in Orion B (Schneider et al. 2013). They show that external compression broadens the column density probability distribution functions, as observed by Tremblin et al. (2014) for RCW 120. On the modeling side, Bron et al. (2018) investigated whether photoevaporation of the illuminated edge of a molecular cloud could explain the high pressures and pressure versus UV field correlation observed at the PDR surface where warm molecular tracers are observed. Using the 1D Hydra PDR code, they show that photoevaporation can produce high thermal pressures (P∕kB ~ 107−108 K cm−3) and the observed P−G0 correlation, almost independently from the initial gas density. They also show that the photoevaporating PDR is preceded by a low velocity shock (a few km s−1) propagating into the molecular cloud. The high pressure resulting from photoevaporation could play a key role in the compression of the molecular material accompagnied by an increase in the local density. The compression is clearly seen in RCW 120 and is also accompagnied by an increase in density, as was also observed by Marsh & Whitworth (2019). The differential column density radial distribution that they present also shows a sharp profile on the interior side of the bubble (see their Fig. 7), which is suggestive of an internal compression responsible for the sharp density enhancement. Deharveng et al. (2009) presented an unsharp mask image, free of extended emission, which was obtained from the Spitzer-MIPSGAL frame at 24 μm near condensation 1 (see their Fig. 10). This image shows distortions of the ionized-dense PDR interface under the influence of the ionizing stars’ radiation observed perpendicular to it and with the same convexity reported here for the filaments.
In RCW 120, the mean intensity profiles of the shell, which were measured toward the shell’s center (internal part) and external part, are presented in Fig. 3 and are both well fitted by a classical Plummer profile. The FWHM of the internal part is 0.09 pc, which is similar to the width found for filaments in low-mass star-forming regions (Arzoumanian et al. 2011, 2019). This similar width observed in different environments suggests a similar origin for the filaments’ formation. Federrath (2016) studied the universality of filaments using magnetohydrodynamical numerical simulations including gravity, turbulence, magnetic fields, as well as jet and outflow feedback. They found that these simulations reproduce realistic filament properties and they also found a universal filament width of 0.1 ±0.02 pc. They propose that this characteristic filament width is based on the sonic scale (λsonic) of the molecular cloud turbulence. This supports the idea that the filament width is determined by the sonic scale and gives a possible explanation for its universality. Using ALMA data, Figueira et al. (2018) have recently shown that fragmentation is limited toward the most massive condensation (condensation 1) in RCW 120 and they suggest that this can be due to the presence of a higher level of turbulence in this region. If so, the compression and turbulence may also have played a key role in the formation offilaments that host the high-mass star formation observed there. Evidence for compression revealed by the asymmetric profile of RCW 120’s shell allows one to possibly link, for the first time, filaments’ properties and star formation at the edges of an H II region.
![]() |
Fig. 5 Mean radial intensity profile (black curve) of the eastern (left) and western (right) parts of the RCW 120 dense shell. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. |
![]() |
Fig. 6 RCW 120: mean radial intensity profile (black curve) of part of the radial filament, which is outlined in purple in Fig. 2 and was measured perpendicular to the filament’s crest in the western part (left) and in the eastern part (right). The green error bars show the standard deviation of the mean intensity profile. The blue solid curve shows the effective beam profile of the ArTéMiS 350 μm data as measured on Mars. The blue dotted curve shows the best-fit Gaussian model to the observed profiles. The red dashed curve shows the best-fit Plummer model convolved with the beam. |
4.2 Symmetric profile of a radial filament
As shown in Fig. 7 and contrary to what is observed for the dense PDR shell, the analysis of part of the southwest radial filament’s profile (see Fig. 2) shows a clear symmetry, which is indicative of no compression for this filament. This suggests that the radiative compression from the H II region is limited to the dense PDR shell and does not extend in the surrounding medium. However, the radial filaments observed around RCW 120 may result from the leaking of radiation coming from the ionized region, as suggested by Deharveng et al. (2009, see their Fig. 16). The leaking radiation may shape the surrounding material in this case but, because the filament is parallel to the incoming radiation, no compressionoccurs, contrary to what happens in the dense shell where the material is perpendicular to the incoming radiation.
4.3 Cloud-cloud collision
Torii et al. (2015) present 12CO, 13CO, and C18 O data of RCW 120 and suggest that the formation of the RCW 120’s ionizing star could result from a cloud-cloud collision. Following this hypothesis, the authors suggest that the ionized region does not result from an expanding shell. According to the proposed scenario (see their Fig. 12), the observed PDR layer located inside the cavity at the interface between the H II region and the molecular cloud can be interpreted as the cavity created in the larger cloud by the collision of a smaller cloud. The inner surface of the interaction zone is illuminated by the strong ultraviolet radiation after the birth of the ionizing star. In this scenario, the compression over the shell is expected to be higher in the interaction zone between the small and large cloud, where the ionizing star has formed. However, we observe a similar asymmetric profile over the dense shell and on its lateral eastern and western parts (thick red lines in Fig. 2), suggesting an isotropic compression over the shell due to the expansion of the ionized region. Anderson et al. (2015) used MOPRA CO data to study the morphology of the RCW 120 ionized region and found no evidence for expansion of the molecular material associated with RCW 120. However, they point out that the lack of detected expansion is roughly in agreement with models for the time-evolution of an H II region like RCW 120, and this is consistent with an expansion speed of ≲ 1.5 km s−1. Forbidden lines, such as [O I] and [C II], are well suited to reveal the complex dynamics of H II regions (Schneider et al. 2018) and will help to shed light on the internal dynamics of RCW 120 (Luisi et al., in prep.).
![]() |
Fig. 7 RCW 120: transverse column density profile of the radial filament outlined in purple in Fig. 2 (black curve). The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. Positive offsets are to the east. Emissions occuring out of the filament are also seen along the radial profile as secondary peaks. |
4.4 3D morphology of RCW 120
The 3D morphology of RCW 120 is still debated. Beaumont & Williams (2010) present CO (3−2) maps of 43 Galactic H II regions and conclude that the shock fronts driven by massive stars tend to create rings instead of spherical shells. They propose that this morphology arises because the host molecular clouds are oblate. They also suggest that such a morphology implies that expanding shock fronts are poorly bound by molecular gas and that cloud compression by these shocks may be limited. However, the ArTéMiS data of RCW 120 clearly reveal an asymmetric profile of the dense shell.
Figure 8 presents the observed median radial column density profile of the dense shell. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown, together with the profile expected for a spherical shell, assuming a classic Plummer profile fitted with RFlat ≃ 0.04 pc and p = 2 (see Sect. 3.3). As observed in Fig. 8, the model of a spherical shell produces too much emission toward the center of the region, compared with the observed profile. The model produces emission that is too low compared to the observed emission of the shell’s external part, but this could be due to extra emission from radial filaments observed in the emission map but which are not included in the model. To reproduce the observed emission, the model should have a very low RFlat and/or much higherp values, which is in disagreement with the observations and results presented in Sect. 3.3. It is also important to note that the intensity level toward the shell’s interior is slightly higher than that observed toward the exterior, indicating that material is present toward the center of the shell. This was also pointed by Deharveng et al. (2010, see their Fig. 12 and Sect. 6.1). Moreover, Iwasaki et al. (2011) present similar density profiles as the one observed here in their 3D simulations of gavitational fragmentation of expanding shells, pointing out the asymmetry of the density profile. The above arguments tend to favor an annular 3D morphology for RCW 120. This morphology is also favored by Kirsanova et al. (2019) who modeled RCW 120 with a ring-like face-on structure to reproduce the observed 13CO (2−1) and C18O (2−1) lines (see also Pavlyuchenkov et al. 2013).
![]() |
Fig. 8 RCW 120: observed mean radial column density profile (black curve) of RCW 120 dense shell, with positive offsets going away from the shell center. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. The green full line shows the profile expected for a spherical shell (see text). |
5 Conclusions
We have presented ArTéMiS 350 and 450 μm emission maps toward the Galactic H II region RCW 120, combined with Herschel data. We reveal the structure of the dense PDR associated with RCW 120 and evidence for the compression of the internal H II region from the asymmetric profile of the mean radial intensity profile of the dense shell. This compression is accompanied by a sharp increase in the column density and ongoing star formation. The width of the filament (0.09 pc) is similar to the one observed in low mass stars forming regions, pointing toward a universal formation mechanism. The analysis of the intensity profiles over the dense shell suggests that RCW 120 could have an annular 3D morphology and that the homogeneous compression of the H II region over the shell does not favor the cloud-cloud collision mechanism proposed for the formation of RCW 120’s ionizing star.
Acknowledgements
A.Z. thanks the support of the Institut Universitaire de France. D.A. and P.P. acknowledge support from FCT through the research grants UIDB/04434/2020 and UIDP/04434/2020. P.P. receives support from fellowship SFRH/BPD/110176/2015 funded by FCT (Portugal) and POPH/FSE (EC).N.S. acknowledge support by the ANR/DFG grant GENESIS (ANR-16-CE92-0035-01/DFG1591/2-1). We thanks ESO, Onsala and the APEX staff in Chile for their support of the ArTéMiS project. We acknowledge the financial support from the French National Research Agency (Grants ANR-05-BLAN-0215 & ANR-11-BS56-0010, and LabEx FOCUS ANR-11-LABX-0013), as well as supports from the French national programs on stellar and ISM physics (PNPS and PCMI). The research leading to these results has received funding from the European Union’s Horizon 2020 research and innovation programme under grant agreement No 730562 [RadioNet]. Part of this work was also supported by the European Research Council under the European Union’s Seventh Framework Programme (ERC Advanced Grant Agreement No. 291294 - “ORISTARS”). This research has made use of data from the Herschel HOBYS project (http://hobys-herschel.cea.fr). HOBYS is a Herschel Key Project jointly carried out by SPIRE Specialist Astronomy Group 3 (SAG3), scientists of the LAM laboratory in Marseille, and scientists of the Herschel Science Center (HSC).
References
- Anderson, L. D., Deharveng, L., Zavagno, A., et al. 2015, ApJ, 800, 101 [NASA ADS] [CrossRef] [Google Scholar]
- André, P. 2017, CR Geosci., 349, 187 [Google Scholar]
- André, P., Di Francesco, J., Ward-Thompson, D., et al. 2014, Protostars and Planets VI (Tucson, AZ: University of Arizona Press), 27 [Google Scholar]
- André, P., Revéret, V., Könyves, V., et al. 2016, A&A, 592, A54 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- André, P., Arzoumanian, D., Könyves, V., Shimajiri, Y., & Palmeirim, P. 2019, A&A, 629, L4 [Google Scholar]
- Arzoumanian, D., André, P., Didelon, P., et al. 2011, A&A, 529, L6 [Google Scholar]
- Arzoumanian, D., André, P., Könyves, V., et al. 2019, A&A, 621, A42 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Beaumont, C. N., & Williams, J. P. 2010, ApJ, 709, 791 [NASA ADS] [CrossRef] [Google Scholar]
- Bernard, A., Neichel, B., Samal, M. R., et al. 2016, A&A, 592, A77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Bron, E., Agúndez, M., Goicoechea, J. R., & Cernicharo, J. 2018, ArXiv e-prints, unpublished [arXiv:1801.01547] [Google Scholar]
- Casali, M. M. 1986, MNRAS, 223, 341 [NASA ADS] [Google Scholar]
- Cox, N. L. J., Arzoumanian, D., André, P., et al. 2016, A&A, 590, A110 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Dawson, J. R., Ntormousi, E., Fukui, Y., Hayakawa, T., & Fierlinger, K. 2015, ApJ, 799, 64 [NASA ADS] [CrossRef] [Google Scholar]
- Deharveng, L., Zavagno, A., & Caplan, J. 2005, A&A, 433, 565 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Deharveng, L., Zavagno, A., Schuller, F., et al. 2009, A&A, 496, 177 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Deharveng, L., Schuller, F., Anderson, L. D., et al. 2010, A&A, 523, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Federrath, C. 2016, MNRAS, 457, 375 [NASA ADS] [CrossRef] [Google Scholar]
- Figueira, M., Zavagno, A., Deharveng, L., et al. 2017, A&A, 600, A93 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Figueira, M., Bronfman, L., Zavagno, A., et al. 2018, A&A, 616, L10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Geen, S., Soler, J. D., & Hennebelle, P. 2017, MNRAS, 471, 4844 [NASA ADS] [CrossRef] [Google Scholar]
- Geen, S., Pellegrini, E., Bieri, R., & Klessen, R. 2020, MNRAS, 492, 915 [NASA ADS] [CrossRef] [Google Scholar]
- Goicoechea, J. R., Gerin, M., & Bron, E. 2019a, BAAS, 51 99 [Google Scholar]
- Goicoechea, J. R., Santa-Maria, M. G., Bron, E., et al. 2019b, A&A, 622, A91 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Gómez, G. C., Vázquez-Semadeni, E., & Zamora-Avilés, M. 2018, MNRAS, 480, 2939 [Google Scholar]
- Henshaw, J. D., Jiménez-Serra, I., Longmore, S. N., et al. 2017, MNRAS, 464, L31 [NASA ADS] [CrossRef] [Google Scholar]
- Hildebrand, R. H. 1983, QJRAS, 24, 267 [NASA ADS] [Google Scholar]
- Inutsuka, S.-i., Inoue, T., Iwasaki, K., & Hosokawa, T. 2015, A&A, 580, A49 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Iwasaki, K., Inutsuka, S.-i., & Tsuribe, T. 2011, ApJ, 733, 17 [Google Scholar]
- Kirsanova, M. S., Pavlyuchenkov, Y. N., Wiebe, D. S., et al. 2019, MNRAS, 488, 5641 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Kohno, M., Tachihara, K., Fujita, S., et al. 2018, PASJ, 126, 1 [Google Scholar]
- Li, G.-X., Urquhart, J. S., Leurini, S., et al. 2016, A&A, 591, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Marsh, K. A., & Whitworth, A. P. 2019, MNRAS, 483, 352 [NASA ADS] [CrossRef] [Google Scholar]
- Martins, F., Pomarès, M., Deharveng, L., Zavagno, A., & Bouret, J. C. 2010, A&A, 510, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Men’shchikov, A. 2013, A&A, 560, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Molinari, S., Bally, J., Glover, S., et al. 2014, Protostars and Planets VI (Tucson, AZ: University of Arizona Press), 125 [Google Scholar]
- Motte, F., Zavagno, A., Bontemps, S., et al. 2010, A&A, 518, L77 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Motte, F., Nony, T., Louvet, F., et al. 2018, Nat. Astron., 2, 478 [NASA ADS] [CrossRef] [Google Scholar]
- Ntormousi, E., Dawson, J. R., Hennebelle, P., & Fierlinger, K. 2017, A&A, 599, A94 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Orkisz, J. H., Pety, J., Gerin, M., et al. 2017, A&A, 599, A99 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ossenkopf-Okada, V., & Stepanov, R. 2019, A&A, 621, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Palmeirim, P., André, P., Kirk, J., et al. 2013, A&A, 550, A38 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Palmeirim, P., Zavagno, A., Elia, D., et al. 2017, A&A, 605, A35 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Pavlyuchenkov, Y. N., Kirsanova, M. S., & Wiebe, D. S. 2013, Astron. Rep., 57, 573 [NASA ADS] [CrossRef] [Google Scholar]
- Peretto, N., André, P., Könyves, V., et al. 2012, A&A, 541, A63 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Roy, A., André, P., Palmeirim, P., et al. 2014, A&A, 562, A138 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Sánchez-Cruces, M., Castellanos-Ramírez, A., Rosado, M., RodríguezGonzález, A., & Reyes-Iturbide, J. 2018, Rev. Mex. Astron. Astrofis., 54, 375 [NASA ADS] [Google Scholar]
- Sault, R. J., Teuben, P. J., & Wright, M. C. H. 1995, ASP Conf. Ser., 77, 433 [Google Scholar]
- Schisano, E., Rygl, K. L. J., Molinari, S., et al. 2014, ApJ, 791, 27 [NASA ADS] [CrossRef] [Google Scholar]
- Schisano, E., Molinari, S., Elia, D., et al. 2020, MNRAS, 492, 5420 [NASA ADS] [CrossRef] [Google Scholar]
- Schneider, S., & Elmegreen, B. G. 1979, ApJS, 41, 87 [NASA ADS] [CrossRef] [Google Scholar]
- Schneider, N., André, P., Könyves, V., et al. 2013, ApJ, 766, L17 [NASA ADS] [CrossRef] [Google Scholar]
- Schneider, N., Bontemps, S., Girichidis, P., et al. 2015, MNRAS, 453, L41 [CrossRef] [Google Scholar]
- Schneider, N., Röllig, M., Simon, R., et al. 2018, A&A, 617, A45 [Google Scholar]
- Sousbie, T. 2011, MNRAS, 414, 350 [NASA ADS] [CrossRef] [Google Scholar]
- Thompson, M. A., Urquhart, J. S., Moore, T. J. T., & Morgan, L. K. 2012, MNRAS, 421, 408 [NASA ADS] [Google Scholar]
- Torii, K., Hasegawa, K., Hattori, Y., et al. 2015, ApJ, 806, 7 [NASA ADS] [CrossRef] [Google Scholar]
- Tremblin, P., Schneider, N., Minier, V., et al. 2014, A&A, 564, A106 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Ungerechts, H., & Thaddeus, P. 1987, ApJS, 63, 645 [NASA ADS] [CrossRef] [Google Scholar]
- Whitworth, A. P., & Ward-Thompson, D. 2001, ApJ, 547, 317 [NASA ADS] [CrossRef] [Google Scholar]
- Xu, S., Ji, S., & Lazarian, A. 2019, ApJ, 878, 157 [NASA ADS] [CrossRef] [Google Scholar]
- Zavagno, A., Pomarès, M., Deharveng, L., et al. 2007, A&A, 472, 835 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zavagno, A., Russeil, D., Motte, F., et al. 2010, A&A, 518, L81 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
- Zhang, M., Kainulainen, J., Mattern, M., Fang, M., & Henning, T. 2019, A&A, 622, A52 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
See http://www.apex-telescope.org/instruments/pi/artemis/ ArTéMiS stands for “ARchitectures de bolomètres pour des TElescopesà grand champ de vue dans le domaine sub-MIllimétrique au Sol” in French.
All Figures
![]() |
Fig. 1 ArTéMiS-only maps of RCW 120 at 350 μm (left) and 450 μm (right) emission observed with ArtéMiS. The white circle shows the location of the ionizing star (Right ascension (J2000) = 258.085833°, Declination (J2000) = −38.490555°). North is up,and east is left. |
In the text |
![]() |
Fig. 2 Combined ArTéMiS+Herschel-SPIRE 350 μm map of RCW 120, with DisPerSe skeleton overlaid in white and Herschel young sources (solid black circles). The white circle shows the location of the ionizing star. The red parts underline the eastern and western lateral side of the dense shell together with the purple radial filament, which have also been analyzed as separate parts (see text). |
In the text |
![]() |
Fig. 3 Median radial intensity profile (black curve) of the RCW 120 ’s shell measured perpendicular to the crest, toward the shell’s center (left) and toward the external part of the shell (right). The green error bars show the standard deviation of the mean intensity profile (same color coding for the other figures). The blue solid curve shows the effective beam profile of the ArTéMiS 350 μm data as measured on Mars. The blue dotted curve shows the best-fit Gaussian model to the observed profiles. The red dashed curve shows the best-fit Plummer model convolved with the beam. |
In the text |
![]() |
Fig. 4 Observed radial column density (and flux density) profile (black curve) of the RCW 120 shell. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. The yellow error bars show the ± 1σ dispersion of the distribution of radial intensity profile observed along the filament crest (same color coding for the other figures). |
In the text |
![]() |
Fig. 5 Mean radial intensity profile (black curve) of the eastern (left) and western (right) parts of the RCW 120 dense shell. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. |
In the text |
![]() |
Fig. 6 RCW 120: mean radial intensity profile (black curve) of part of the radial filament, which is outlined in purple in Fig. 2 and was measured perpendicular to the filament’s crest in the western part (left) and in the eastern part (right). The green error bars show the standard deviation of the mean intensity profile. The blue solid curve shows the effective beam profile of the ArTéMiS 350 μm data as measured on Mars. The blue dotted curve shows the best-fit Gaussian model to the observed profiles. The red dashed curve shows the best-fit Plummer model convolved with the beam. |
In the text |
![]() |
Fig. 7 RCW 120: transverse column density profile of the radial filament outlined in purple in Fig. 2 (black curve). The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. Positive offsets are to the east. Emissions occuring out of the filament are also seen along the radial profile as secondary peaks. |
In the text |
![]() |
Fig. 8 RCW 120: observed mean radial column density profile (black curve) of RCW 120 dense shell, with positive offsets going away from the shell center. The Gaussian fit (blue dotted curve) and the Plummer fit (red dashed curve) are shown. The green full line shows the profile expected for a spherical shell (see text). |
In the text |
Current usage metrics show cumulative count of Article Views (full-text article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 48-96 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.