Context. Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. Aims. We perform a case ...study of whether low-mass filaments are close to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. To accomplish this, we study the ~6.5 pc long Musca molecular cloud, which is an ideal candidate for a filament at an early stage of fragmentation. Methods. We employ dust extinction mapping, in conjunction with near-infrared JHKS-band data from the CTIO/NEWFIRM instrument, and 870 μm dust continuum emission data from the APEX/LABOCA instrument to estimate column densities in Musca. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud’s morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. Results. The Musca cloud is pronouncedly fragmented at its ends, but harbors a remarkably well-defined, ~1.6 pc long filament in its center region. The line mass of the filament is 21–31 M⊙ pc-1 and the full width at half maximum (FWHM) 0.07 pc. The radial profile of the filament can be fitted with a Plummer profile, which has the power-index of 2.6 ± 11% and is flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of ~5–10 × 104 cm-3. The fragments in the cloud have a mean separation of ~0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation started a few tenths of a Myr ago from the ends of the cloud, leaving its center still relatively nonfragmented, possibly because of gravitational focusing in a finite geometry.
Context. An increasing number of hundred-parsec-scale, high line-mass filaments are being detected in the Galaxy. Their evolutionary path, including fragmentation towards star formation, is virtually ...unknown. Aims. We characterize the fragmentation within the hundred-parsec-scale, high line-mass Nessie filament, covering size-scales in the range ~0.1–100 pc. We also connect the small-scale fragments to the star-forming potential of the cloud. Methods. We combine near-infrared data from the VISTA Variables in the Via Lactea (VVV) survey with mid-infrared Spitzer/GLIMPSE data to derive a high-resolution dust extinction map for Nessie. We then apply a wavelet decomposition technique on the map to analyze the fragmentation characteristics of the cloud. The characteristics are then compared with predictions from gravitational fragmentation models. We compare the detected objects to those identified at a resolution approximately ten times lower from ATLASGAL 870 μm dust emission data. Results. We present a high-resolution extinction map of Nessie (2″ full-width-half-max, FWHM, corresponding to 0.03 pc). We estimate the mean line mass of Nessie to be ~627 M⊙ pc−1 and the distance to be ~3.5 kpc. We find that Nessie shows fragmentation at multiple size scales. The median nearest-neighbor separations of the fragments at all scales are within a factor of two of the Jeans’ length at that scale. However, the relationship between the mean densities of the fragments and their separations is significantly shallower than expected for Jeans’ fragmentation. The relationship is similar to the one predicted for a filament that exhibits a Larson-like scaling between size-scale and velocity dispersion; such a scaling may result from turbulent support. Based on the number of young stellar objects (YSOs) in the cloud, we estimate that the star formation rate (SFR) of Nessie is ~371 M⊙ Myr−1; similar values result if using the number of dense cores, or the amount of dense gas, as the proxy of star formation. The star formation efficiency is 0.017. These numbers indicate that by its star-forming content, Nessie is comparable to the Solar neighborhood giant molecular clouds like Orion A.
Filaments play a central role in the molecular clouds’ evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we ...have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in 13CO and C18O (2–1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its ~6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e. σNT/cs ≲ 1) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the interstellar medium. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson’s velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.
Aims. Revealing the fragmentation, infall, and outflow processes in the immediate environment around massive young stellar objects is crucial for understanding the formation of the most massive ...stars. Methods. With this goal in mind we present the so far highest spatial-resolution thermal submm line and continuum observations toward the young high-mass protostar NGC 7538 IRS1. Using the Plateau de Bure Interferometer in its most extended configuration at 843 μm wavelength, we achieved a spatial resolution of 0.2″ × 0.17″, corresponding to ~500 AU at a distance of 2.7 kpc. Results. For the first time, we have observed the fragmentation of the dense inner core of this region with at least three subsources within the inner 3000 AU. The outflow exhibits blue- and red-shifted emission on both sides of the central source, indicating that the current orientation has to be close to the line-of-sight, which differs from other recent models. We observe rotational signatures in northeast-southwest direction; however, even on scales of 500 AU, we do not identify any Keplerian rotation signatures. This implies that during the early evolutionary stages any stable Keplerian inner disk has to be very small (≤500) AU). The high-energy line HCN(4−3)v2 = 1 (E − u/k = 1050 K) is detected over an extent of approximately 3000 AU. In addition to this, the detection of red-shifted absorption from this line toward the central dust continuum peak position allows us to estimate infall rates of ~1.8 × 10-3 M⊙ yr-1 on the smallest spatial scales. Although all that gas will not necessarily be accreted onto the central protostar, nevertheless, such inner core infall rates are among the best proxies of the actual accretion rates one can derive during the early embedded star formation phase. These data are consistent with collapse simulations and the observed high multiplicity of massive stars.
Understanding the chemical evolution of young (high-mass) star-forming regions is a central topic in star formation research. Chemistry is employed as a unique tool 1) to investigate the underlying ...physical processes and 2) to characterize the evolution of the chemical composition. With these aims in mind, we observed a sample of 59 high-mass star-forming regions at different evolutionary stages varying from the early starless phase of infrared dark clouds to high-mass protostellar objects to hot molecular cores and, finally, ultra-compact Hii regions at 1 mm and 3 mm with the IRAM 30 m telescope. We determined their large-scale chemical abundances and found that the chemical composition evolves along with the evolutionary stages. On average, the molecular abundances increase with time. We modeled the chemical evolution, using a 1D physical model where density and temperature vary from stage to stage coupled with an advanced gas-grain chemical model and derived the best-fit χ2 values of all relevant parameters. A satisfying overall agreement between observed and modeled column densities for most of the molecules was obtained. With the best-fit model we also derived a chemical age for each stage, which gives the timescales for the transformation between two consecutive stages. The best-fit chemical ages are ~10 000 years for the IRDC stage, ~60 000 years for the HMPO stage, ~40 000 years for the HMC stage, and ~10 000 years for the UCHii stage. Thus, the total chemical timescale for the entire evolutionary sequence of the high-mass star formation process is on the order of 105 years, which is consistent with theoretical estimates. Furthermore, based on the approach of a multiple-line survey of unresolved data, we were able to constrain an intuitive and reasonable physical and chemical model. The results of this study can be used as chemical templates for the different evolutionary stages in high-mass star formation.
Context. The formation processes and the evolutionary stages of high-mass stars are poorly understood compared to low-mass stars. Large-scale surveys are needed to provide an unbiased census of high ...column density sites that can potentially host precursors to high-mass stars. Aims. The ATLASGAL survey covers 420 sq. degree of the Galactic plane, between −80° < ℓ < +60° at 870 μm. Here we identify the population of embedded sources throughout the inner Galaxy. With this catalog we first investigate the general statistical properties of dust condensations in terms of their observed parameters, such as flux density and angular size. Then using mid-infrared surveys we aim to investigate their star formation activity and the Galactic distribution of star-forming and quiescent clumps. Our ultimate goal is to determine the statistical properties of quiescent and star-forming clumps within the Galaxy and to constrain the star formation processes. Methods. We optimized the source extraction method, referred to as MRE-GCL, for the ATLASGAL maps in order to generate a catalog of compact sources. This technique is based on multiscale filtering to remove extended emission from clouds to better determine the parameters corresponding to the embedded compact sources. In a second step we extracted the sources by fitting 2D Gaussians with the Gaussclumps algorithm. Results. We have identified in total 10861 compact submillimeter sources with fluxes above 5σ. Completeness tests show that this catalog is 97% complete above 5σ and >99% complete above 7σ. Correlating this sample of clumps with mid-infrared point source catalogs (MSX at 21.3 μm and WISE at 22 μm), we have determined a lower limit of 33% that is associated with embedded protostellar objects. We note that the proportion of clumps associated with mid-infrared sources increases with increasing flux density, achieving a rather constant fraction of ~75% of all clumps with fluxes over 5 Jy/beam being associated with star formation. Examining the source counts as a function of Galactic longitude, we are able to identify the most prominent star-forming regions in the Galaxy. Conclusions. We present here the compact source catalog of the full ATLASGAL survey and investigate their characteristic properties. From the fraction of the likely massive quiescent clumps (~25%), we estimate a formation time scale of ~ 7.5 ± 2.5 × 104 yr for the deeply embedded phase before the emergence of luminous young stellar objects. Such a short duration for the formation of high-mass stars in massive clumps clearly proves that the earliest phases have to be dynamic with supersonic motions.
Context.
The Galactic plane has been observed extensively by a large number of Galactic plane surveys from infrared to radio wavelengths at an angular resolution below 40′′. However, a 21 cm line and ...continuum survey with comparable spatial resolution is lacking.
Aims.
The first half of THOR data (
l
= 14.0°−37.9°, and
l
= 47.1°−51.2°, |
b
|≤ 1.25°) has been published in our data release 1 paper. With this data release 2 paper, we publish all the remaining spectral line data and Stokes I continuum data with high angular resolution (10′′–40′′), including a new H
I
dataset for the whole THOR survey region (
l
= 14.0−67.4° and |
b
|≤ 1.25°). As we published the results of OH lines and continuum emission elsewhere, we concentrate on the H
I
analysis in this paper.
Methods.
With the
Karl G. Jansky
Very Large Array (VLA) in C-configuration, we observed a large portion of the first Galactic quadrant, achieving an angular resolution of ≤40′′. At
L
Band, the WIDAR correlator at the VLA was set to cover the 21 cm H
I
line, four OH transitions, a series of H
nα
radio recombination lines (RRLs;
n
= 151 to 186), and eight 128 MHz-wide continuum spectral windows, simultaneously.
Results.
We publish all OH and RRL data from the C-configuration observations, and a new H
I
dataset combining VLA C+D+GBT (VLA D-configuration and GBT data are from the VLA Galactic Plane Survey) for the whole survey. The H
I
emission shows clear filamentary substructures at negative velocities with low velocity crowding. The emission at positive velocities is more smeared-out, likely due to higher spatial and velocity crowding of structures at the positive velocities. Compared to the spiral arm model of the Milky Way, the atomic gas follows the Sagittarius and Perseus Arm well, but with significant material in the inter-arm regions. With the C-configuration-only H
I
+continuum data, we produce an H
I
optical depth map of the THOR areal coverage from 228 absorption spectra with the nearest-neighbor method. With this
τ
map, we corrected the H
I
emission for optical depth, and the derived column density is 38% higher than the column density with optically thin assumption. The total H
I
mass with optical depth correction in the survey region is 4.7 × 10
8
M
⊙
, 31% more than the mass derived assuming the emission is optically thin. If we applied this 31% correction to the whole Milky Way, the total atomic gas mass would be 9.4–10.5 × 10
9
M
⊙
. Comparing the H
I
with existing CO data, we find a significant increase in the atomic-to-molecular gas ratio from the spiral arms to the inter-arm regions.
Conclusions.
The high-sensitivity and resolution THOR H
I
dataset provides an important new window on the physical and kinematic properties of gas in the inner Galaxy. Although the optical depth we derive is a lower limit, our study shows that the optical depth correction issignificant for H
I
column density and mass estimation. Together with the OH, RRL and continuum emission from the THOR survey, these new H
I
data provide the basis for high-angular-resolution studies of the interstellar medium in different phases.
ABSTRACT
We present a study of the cold atomic hydrogen (H i) content of molecular clouds simulated within the SILCC-Zoom project for solar neighbourhood conditions. We produce synthetic observations ...of H i at 21 cm, including H i self-absorption (HISA) and observational effects. We find that H i column densities, $N_{\rm H\, \small {\rm I}}$, of ≳1022 cm−2 are frequently reached in molecular clouds with H i temperatures as low as ∼10 K. Hence, HISA observations assuming a fixed H i temperature tend to underestimate the amount of cold H i in molecular clouds by a factor of 3–10 and produce an artificial upper limit of $N_{\rm H\, \small {\rm I}}$ around 1021 cm−2. We thus argue that the cold H i mass in molecular clouds could be a factor of a few higher than previously estimated. Also, $N_{\rm H\, \small {\rm I}}$ PDFs obtained from HISA observations might be subject to observational biases and should be considered with caution. The underestimation of cold H i in HISA observations is due to both the large H i temperature variations and the effect of noise in regions of high optical depth. We find optical depths of cold H i around 1–10, making optical depth corrections essential. We show that the high H i column densities (≳1022 cm−2) can in parts be attributed to the occurrence of up to 10 individual H i–H2 transitions along the line of sight. This is also reflected in the spectra, necessitating Gaussian decomposition algorithms for their in-depth analysis. However, also for a single H i–H2 transition, $N_{\rm H\, \small {\rm I}}$ frequently exceeds 1021 cm−2, challenging one-dimensional, semi-analytical models. This is due to non-equilibrium chemistry effects and the fact that H i–H2 transition regions usually do not possess a one-dimensional geometry. Finally, we show that the H i gas is moderately supersonic with Mach numbers of a few. The corresponding non-thermal velocity dispersion can be determined via HISA observations within a factor of ∼2.
We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass ~400 M⊙ pc-1) in the Orion A molecular cloud. We have observed a 1.6 pc long section of ...the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of ~3″ (1200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below ~17 000 AU and strongly grouped below ~6000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below ~6000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of ~30 000 AU in size, separated by ~50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments (~16 M⊙ pc-1) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below ~17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks. This is in agreement with a scenario where protostars are ejected from the maternal filament by the slingshot mechanism, a model recently proposed for the ISF. The separation distributions of the dense cores and protostars may also provide an evolutionary tracer of filament fragmentation.
Context. Filamentary structures in the interstellar medium are crucial ingredients of the star formation process. They fragment to form individual star-forming cores, and at the same time they may ...also funnel gas toward the central gas cores, providing an additional gas reservoir. Aims. We want to resolve the length scales for filament formation and fragmentation (resolution ≤0.1 pc), in particular the Jeans length and cylinder fragmentation scale. Methods. We have observed the prototypical high-mass star-forming filament IRDC 18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2 mm continuum and N2H+(1–0) line emission in a ten-field mosaic at a spatial resolution of ~ 4′′ (~14 000 au). Results. The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(± 0.18) pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, and gravitationally bound gas cylinder, a high mass-to-length ratio of M/l ≈ 1000 M⊙ pc-1 requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1−0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent subfilaments. Furthermore, this velocity gradient perpendicular to the filament resembles results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers (CI and C18O) reveal a similar red- and blueshifted velocity structure on scales around 60′′ east and west of the filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may be due to no significant gas flows along the filamentary axis, but it may also be partly caused by a low inclination angle of the filament with respect to the plane of the sky minimizing such a signature. Conclusions. The IRDC 18223 3.2 mm continuum data are consistent with thermal fragmentation of a gravitationally bound and compressible gas cylinder. However, the high mass-to-length ratio requires additional support – most likely turbulence and/or magnetic fields – against collapse. The N2H+ spectral line data indicate a kinematic origin of the filament, but we cannot conclusively differentiate whether it has formed out of (pre-existing) velocity-coherent subfilaments, whether magnetized converging gas flows, a larger-scale collapsing cloud, or even whether rotation played a significant role during filament formation.