High-mass Star-forming Clumps Research Articles

Digging into the Interior of Hot Cores with ALMA (DIHCA). IV. Fragmentation in High-mass Star-forming Clumps

Fragmentation contributes to the formation and evolution of stars. Observationally, high-mass stars are known to form multiple-star systems, preferentially in cluster environments. Theoretically, Jeans instability has been suggested to determine characteristic fragmentation scales, and thermal or turbulent motion in the parental gas clump mainly contributes to the instability. To search for such a characteristic fragmentation scale, we have analyzed Atacama Large Millimeter/submillimeter Array (ALMA) 1.33 mm continuum observations toward 30 high-mass star-forming clumps taken by the Digging into the Interior of Hot Cores with ALMA survey. We have identified 573 cores using the dendrogram algorithm and measured the separation of cores by using the Minimum Spanning Tree technique. The core separation corrected by projection effects has a distribution peaked around 5800 au. In order to remove biases produced by different distances and sensitivities, we further smooth the images to a common physical scale and perform completeness tests. Our careful analysis finds a characteristic fragmentation scale of ∼7000 au, comparable to the thermal Jeans length of the clumps. We conclude that thermal Jeans fragmentation plays a dominant role in determining the clump fragmentation in high-mass star-forming regions, without the need to invoke turbulent Jeans fragmentation.

ATOMS: ALMA three-millimetre observations of massive star-forming regions – XVI. Neutral versus ion line widths

ABSTRACT It has been suggested that the line width of ions in molecular clouds is narrower than that of the co-existing neutral particles, which has been interpreted as an indication of the decoupling of neutral turbulence from magnetic fields within a partially ionized environment. We calculate the principal component analysis (PCA) correlation coefficients of CCH versus H$^{13}$CO$^{+}$ and H$^{13}$CN versus H$^{13}$CO$^{+}$. We find aside from H$^{13}$CN, CCH could also be strongly spatial correlated with H$^{13}$CO$^{+}$ in high-mass star-forming regions. CCH and H$^{13}$CO$^{+}$ line emissions are strongly spatial correlated with each other in 48 per cent sources with a PCA correlation coefficient over 0.7. So, we investigate the ambipolar diffusion (AD) effect using CCH and H$^{13}$CO$^{+}$ lines as a neutral/ion pair in a sample of 129 high-mass star-forming clumps. We conduct a careful analysis of line widths of the CCH–H$^{13}$CO$^{+}$ pair pixel-by-pixel in 12 sources, which show a strong correlation in CCH–H$^{13}$CO$^+$ emission and no obvious outflows or multiple velocity components. The mean velocity dispersion of CCH is about the same as H$^{13}$CO$^{+}$ in 12 sources. In low-density regions of most sources, CCH shows a broader velocity dispersion than H$^{13}$CO$^{+}$. However, the AD effect is not significant from a statistical point of view.

The ALMA Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES). XI. Statistical Study of Early Fragmentation

Fragmentation during the early stages of high-mass star formation is crucial for understanding the formation of high-mass clusters. We investigated fragmentation within 39 high-mass star-forming clumps as part of the Atacama Large Millimeter/submillimeter Array Survey of 70 μm Dark High-mass Clumps in Early Stages (ASHES) survey. Considering projection effects, we have estimated core separations for 839 cores identified from the continuum emission and found mean values between 0.08 and 0.32 pc within each clump. We find compatibility of the observed core separations and masses with the thermal Jeans length and mass, respectively. We also present subclump structures revealed by the 7 m array continuum emission. Comparison of the Jeans parameters using clump and subclump densities with the separation and masses of gravitationally bound cores suggests that they can be explained by clump fragmentation, implying the simultaneous formation of subclumps and cores within rather than a step-by-step hierarchical fragmentation. The number of cores in each clump positively correlates with the clump surface density and the number expected from the thermal Jeans fragmentation. We also find that the higher the fraction of protostellar cores, the larger the dynamic range of the core mass, implying that the cores are growing in mass as the clump evolves. The ASHES sample exhibits various fragmentation patterns: aligned, scattered, clustered, and subclustered. Using the Q -parameter, which can help distinguish between centrally condensed and subclustered spatial core distributions, we finally find that in the early evolutionary stages of high-mass star formation, cores tend to follow a subclustered distribution.

The SEDIGISM survey: A search for molecular outflows

Context. The formation processes of massive stars are still unclear, but a picture is emerging involving accretion disks and molecular outflows in what appears to be a scaled-up version of low-mass star formation. A census of outflow activity toward high-mass star-forming clumps in various evolutionary stages has the potential to shed light on high-mass star formation. Aims. We conducted an outflow survey toward ATLASGAL (APEX Telescope Large Area Survey of the Galaxy) clumps using SEDIGISM (structure, Excitation, and Dynamics of the Inner Galactic InterStellar Medium) data and aimed to obtain a large sample of clumps exhibiting outflow activity in different evolutionary stages. Methods. We identify the high-velocity wings of the 13CO lines, which indicate outflow activity, toward ATLASGAL clumps by (1) extracting the simultaneously observed 13CO (2–1) and C18O (2–1) spectra from SEDIGISM, and (2) subtracting Gaussian fits to the scaled C18O (core emission) from the 13CO line after considering opacity broadening. Results. We detected high-velocity gas toward 1192 clumps out of a total sample of 2052, corresponding to an overall detection rate of 58%. Outflow activity has been detected in the earliest (apparently) quiescent clumps (i.e., 70 μm weak) to the most evolved H II region stages (i.e., 8 μm bright with tracers of massive star formation). The detection rate increases as a function of evolution (quiescent = 51%, protostellar = 47%, YSO = 57%, UC H II regions = 76%). Conclusions. Our sample is the largest outflow sample identified so far. The high detection rate from this large sample is consistent with the results of similar studies reported in the literature and supports the scenario that outflows are a ubiquitous feature of high-mass star formation. The lower detection rate in early evolutionary stages may be due to the fact that outflows in the early stages are weak and difficult to detect. We obtain a statistically significant sample of outflow clumps for every evolutionary stage, especially for outflow clumps in the earliest stage (i.e., 70 μm dark). The detections of outflows in the 70 μm dark clumps suggest that the absence of 70 μm emission is not a robust indicator of starless and/or pre-stellar cores.

Core Mass Function of a Single Giant Molecular Cloud Complex with ∼10,000 Cores

Abstract Similarity in shape between the initial mass function (IMF) and the core mass functions (CMFs) in star-forming regions prompts the idea that the IMF originates from the CMF through a self-similar core-to-star mass mapping process. To accurately determine the shape of the CMF, we create a sample of 8431 cores with the dust continuum maps of the Cygnus X giant molecular cloud complex, and design a procedure for deriving the CMF considering the mass uncertainty, binning uncertainty, sample incompleteness, and the statistical errors. The resultant CMF coincides well with the IMF for core masses from a few M ⊙ to the highest masses of 1300 M ⊙ with a power-law of dN / dM ∝ M − 2.30 ± 0.04 , but does not present an obvious flattened turnover in the low-mass range as the IMF does. More detailed inspection reveals that the slope of the CMF steepens with increasing mass. Given the numerous high-mass star-forming activities of Cygnus X, this is in stark contrast with the existing top-heavy CMFs found in high-mass star-forming clumps. We also find that the similarity between the IMF and the mass function of cloud structures is not unique at core scales, but can be seen for cloud structures of up to several parsec scales. Finally, our SMA observations toward a subset of the cores do not present evidence for the self-similar mapping. The latter two results indicate that the shape of the IMF may not be directly inherited from the CMF.

The 3D Structure of CO Depletion in High-mass Prestellar Regions

Abstract Disentangling the different stages of the star formation process, in particular in the high-mass regime, is a challenge in astrophysics. Chemical clocks could help alleviate this problem, but their evolution strongly depends on many parameters, leading to degeneracy in the interpretation of the observational data. One of these uncertainties is the degree of CO depletion. We present here the first self-consistent magneto-hydrodynamic simulations of high-mass, star-forming regions at different scales, fully coupled with a nonequilibrium chemical network, which includes C–N–O bearing molecules. Depletion and desorption processes are treated time dependently. The results show that full CO depletion (i.e., all gas-phase CO frozen-out on the surface of dust grains) can be reached very quickly, in one-third or even smaller fractions of the freefall time, whether the collapse proceeds on slow or fast timescales. This leads to a high level of deuteration in a short time, both for typical tracers like N2H+, as well as for the main ion H 3 + , the latter being in general larger and more extended. N2 depletion is slightly less efficient, and no direct effects on N-bearing molecules and deuterium fractionation are observed. We show that CO depletion is not the only driver of deuteration, and that there is a strong impact on D frac when changing the grain size. We finally apply a two-dimensional Gaussian point-spread function to our results to mimic observations with single-dish and interferometers. Our findings suggest that the low-values observed in high-mass star-forming clumps are in reality masking a full-depletion stage in the inner 0.1 pc region.

ATLASGAL – molecular fingerprints of a sample of massive star-forming clumps★

We have conducted a 3-mm molecular-line survey towards 570 high-mass star-forming clumps, using the Mopra telescope. The sample is selected from the 10,000 clumps identified by the ATLASGAL survey and includes all of the most important embedded evolutionary stages associated with massive star formation, classified into five distinct categories (quiescent, protostellar, young stellar objects, \hii\ regions and photo-dominated regions). The observations were performed in broadband mode with frequency coverage of 85.2 to 93.4\,GHz and a velocity resolution of $\sim$0.9\,\kms, detecting emission from 26 different transitions. We find significant evolutionary trends in the detection rates, integrated line intensities, and abundances of many of the transitions and also identify a couple of molecules that appear to be invariant to changes in the dust temperature and evolutionary stage (N$_2$H$^+$\,(1-0) and HN$^{13}$C\,(1-0)). We use the K-ladders for CH$_3$C$_2$H\,(5-4) and CH$_3$CH\,(5-4) to calculate the rotation temperatures and find $\sim$1/3 of the quiescent clumps have rotation temperatures that suggest the presence of an internal heating source. These sources may constitute a population of very young protostellar objects that are still dark at 70\,\mum\ and suggest that the fraction of truly quiescent clumps may only be a few per cent. We also identify a number of line ratios that show a strong correlation with the evolutionary stage of the embedded objects and discuss their utility as diagnostic probes of evolution.

Chemical Evolution of N2H+ in Six Massive Star-forming Regions

Abstract To investigate how the abundance of N2H+ varies as massive clumps evolve, here we present a multiwavelength study toward six molecular clouds. All of these clouds contain several massive clumps in different evolutionary stages of star formation. Using archival data of the Herschel infrared Galactic Plane Survey (Hi-GAL), we made H2 column density and dust temperature maps of these regions by the spectral energy distribution method. We found that all of the six clouds show distinct dust temperature gradients, ranging from ∼20 to ∼30 K. This makes them good candidates to study chemical evolution of molecules (such as N2H+) in different evolutionary stages of star formation. Our molecular line data comes from the Millimeter Astronomy Legacy Team Survey at 90 GHz (MALT90). We made column density and then abundance maps of N2H+. We found that when the dust temperature is above 27 K, the abundance of N2H+ begins to decrease or reaches a plateau. We regard that this is because in the photodissociation regions around classical H ii regions, N2H+ is heavily destroyed by free electrons. However, when the dust temperature is below 27 K, the abundance of N2H+ increases with the dust temperature. This seems to be inconsistent with previous chemical models made in low-mass star-forming regions. In order to investigate whether this inconsistency is caused by a different chemistry in high-mass star-forming clumps, higher angular resolution observations are necessary.

Protonated CO2 in massive star-forming clumps

ABSTRACT Interstellar CO2 is an important reservoir of carbon and oxygen, and one of the major constituents of the icy mantles of dust grains, but it is not observable directly in cold gas because it has no permanent dipole moment. Its protonated form, HOCO+, is believed to be a good proxy for gaseous CO2. However, it has been detected in only a few star-forming regions to date, so its interstellar chemistry is not well understood. We present new detections of HOCO+ lines in 11 high-mass star-forming clumps. Our observations more than treble the number of detections in star-forming regions to date. We derived beam-averaged abundances relative to H2 of between 0.3 and 3.8 × 10−11. We compared these values with the abundances of H13CO+, a possible gas-phase precursor of HOCO+, and of CH3OH, a product of surface chemistry. We found a positive correlation with H13CO+, but no correlation with CH3OH. We suggest that the gas-phase formation route starting from HCO+ plays an important role in the formation of HOCO+, perhaps more important than protonation of CO2 (upon evaporation of CO2 from icy dust mantles).

New detections of (sub)millimeter hydrogen radio recombination lines towards high-mass star-forming clumps

Aims. Previous radio recombination line (RRL) observations of dust clumps identified in the APEX Telescope Large Area Survey of the Galaxy (ATLASGAL) have led to the detection of a large number of RRLs in the 3 mm range. Here, we aim to study their excitation with shorter wavelength (sub)millimeter radio recombination line (submm-RRL) observations. Methods. We made observations of submm-RRLs with low principal quantum numbers (n ≤ 30) using the APEX 12 m telescope, toward 104 H II regions associated with massive dust clumps from ATLASGAL. The observations covered the H25α, H28α, and H35β transitions. Toward a small subsample the H26α, H27α, H29α, and H30α lines were observed to avoid contamination by molecular lines at adjacent frequencies. Results. We have detected submm-RRLs (signal-to-noise (S∕N)≥ 3 σ) from compact H II regions embedded within 93 clumps. The submm-RRLs are approximately a factor of two brighter than the mm-RRLs and consistent with optically thin emission in local thermodynamic equilibrium (LTE). The average ratio (0.31) of the measured H35β/H28α fluxes is close to the LTE value of 0.28. No indication of RRL maser emission has been found. The Lyman photon flux, bolometric, and submm-RRL luminosities toward the submm-RRL detected sources present significant correlations. The trends of dust temperature and the ratio of bolometric luminosity to clump mass, Lbol ∕Mclump, indicate that the H II regions are related to the most massive and luminous clumps. By estimating the production rate of ionizing photons, Q, from the submm-RRL flux, we find that the Q(H28α) measurements provide estimates of the Lyman continuum photon flux consistent with those determined from 5 GHz radio continuum emission. Six RRL sources show line profiles that are a combination of a narrow and a broad Gaussian feature. The broad features are likely associated with high-velocity ionized flows. Conclusions. We have detected submm-RRLs toward 93 ATLASGAL clumps. Six RRL sources have high-velocity RRL components likely driven by high-velocity ionized flows. Their observed properties are consistent with thermal emission that correlates well with the Lyman continuum flux of the H II regions. The sample of H II regions with mm/submm-RRL detections probes, in our Galaxy, luminous clumps (Lbol > 104 L⊙) with high Lbol∕Mclump. We also provide suitable candidates for further studies of the morphology and kinematics of embedded, compact H II regions with the Atacama Large Millimeter/submillimeter Array (ALMA).

Characterizing the properties of cluster precursors in the MALT90 survey

In the Milky Way there are thousands of stellar clusters each harboring from a hundred to a million stars. Although clusters are common, the initial conditions of cluster formation are still not well understood. To determine the processes involved in the formation and evolution of clusters it is key to determine the global properties of cluster-forming clumps in their earliest stages of evolution. Here, we present the physical properties of 1,244 clumps identified from the MALT90 survey. Using the dust temperature of the clumps as a proxy for evolution we determined how the clump properties change at different evolutionary stages. We find that less-evolved clumps exhibiting dust temperatures lower than 20 K have higher densities and are more gravitationally bound than more-evolved clumps with higher dust temperatures. We also identified a sample of clumps in a very early stage of evolution, thus potential candidates for high-mass star-forming clumps. Only one clump in our sample has physical properties consistent with a young massive cluster progenitor, reinforcing the fact that massive proto-clusters are very rare in the Galaxy.

LINKING DENSE GAS FROM THE MILKY WAY TO EXTERNAL GALAXIES

ABSTRACT In a survey of 65 galaxies, Gao & Solomon found a tight linear relation between the infrared luminosity (L IR, a proxy for the star formation rate) and the HCN(1–0) luminosity ( ). Wu et al. found that this relation extends from these galaxies to the much less luminous Galactic molecular high-mass star-forming clumps (∼1 pc scales), and posited that there exists a characteristic ratio L IR/ for high-mass star-forming clumps. The Gao–Solomon relation for galaxies could then be explained as a summation of large numbers of high-mass star-forming clumps, resulting in the same L IR/ ratio for galaxies. We test this explanation and other possible origins of the Gao–Solomon relation using high-density tracers (including HCN(1–0), N2H+(1–0), HCO+(1–0), HNC(1–0), HC3N(10–9), and C2H(1–0)) for ∼300 Galactic clumps from the Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey. The MALT90 data show that the Gao–Solomon relation in galaxies cannot be satisfactorily explained by the blending of large numbers of high-mass clumps in the telescope beam. Not only do the clumps have a large scatter in the L IR/ ratio, but also far too many high-mass clumps are required to account for the Galactic IR and HCN luminosities. We suggest that the scatter in the L IR/ ratio converges to the scatter of the Gao–Solomon relation at some size-scale ≳1 kpc. We suggest that the Gao–Solomon relation could instead result from of a universal large-scale star formation efficiency, initial mass function, core mass function, and clump mass function.

The global chemical properties of high-mass star forming clumps at different evolutionary stages

A total of 197 relatively isolated high-mass star-forming clumps were selected from the Millimeter Astronomy Legacy Team 90 GHz (MALT90) survey data and their global chemical evolution investigated using four molecular lines, \(\mbox{N}_{2}\mbox{H}^{+}\ \mbox{(1--0)}\), \(\mbox{HCO}^{+}\ \mbox{(1--0)}\), HCN (1–0), and HNC (1–0). The results suggest that the global averaged integrated intensity ratios \(I(\mbox{HCO}^{+})/I(\mbox{HNC})\), \(I(\mbox{HCN})/I(\mbox{HNC})\), \(I(\mbox{N}_{2}\mbox{H}^{+})\mbox{/}I(\mbox{HCO}^{+})\), and \(I(\mbox{N}_{2}\mbox{H}^{+})\mbox{/} I(\mbox{HCN})\) are promising tracers for evolution of high-mass star-forming clumps. The global averaged column densities and abundances of \(\mbox{N}_{2}\mbox{H}^{+}\), \(\mbox{HCO}^{+}\), \(\mbox{HCN}\), and \(\mbox{HNC}\) increase as clumps evolve. The global averaged abundance ratios \(X(\mbox{HCN})\mbox{/}X(\mbox{HNC})\) could be used to trace evolution of high-mass star forming clumps, \(X(\mbox{HCO}^{+})\mbox{/}X(\mbox{HNC})\) is more suitable for distinguishing high-mass star-forming clumps in prestellar (stage A) from those in protostellar (stage B) and HII/PDR region (stage C). These results suggest that the global averaged integrated intensity ratios between HCN (1–0), HNC (1–0), \(\mbox{HCO}^{+}\ \mbox{(1--0)}\) and \(\mbox{N}_{2}\mbox{H}^{+}\ \mbox{(1--0)}\) are more suitable for tracing the evolution of high-mass star forming clumps. We also studied the chemical properties of the target high-mass star-forming clumps in each spiral arm of the Galaxy, and got results very different from those above. This is probably due to the relatively small sample in each spiral arm. For high-mass star-forming clumps in Sagittarius arm and Norma-Outer arm, comparing two groups located on one arm with different Galactocentric distances, the clumps near the Galactic Center appear to be younger than those far from the Galactic center, which may be due to more dense gas concentrated near the Galactic Center, and hence more massive stars being formed there.

Infall through the evolution of high-mass star-forming clumps

With the GREAT receiver at the Stratospheric Observatory for Infrared Astronomy (SOFIA), nine massive molecular clumps have been observed in the ammonia $3_{2+}- 2_{2-}$ line at 1.8~THz in a search for signatures of infall. The sources were selected from the ATLASGAL submillimeter dust continuum survey of our Galaxy. Clumps with high masses covering a range of evolutionary stages based on their infrared properties were chosen. The ammonia line was detected in all sources, leading to five new detections and one confirmation of a previous detection of redshifted absorption in front of their strong THz continuum as a probe of infall in the clumps. These detections include two clumps embedded in infrared dark clouds. The measured velocity shifts of the absorptions compared to optically thin \CSEO\ (3--2) emission are 0.3--2.8~km/s, corresponding to fractions of 3\%\ to 30\% of the free-fall velocities of the clumps. The ammonia infall signature is compared with complementary data of different transitions of HCN, HNC, CS, and HCO$^+$, which are often used to probe infall via their blue-skewed line profiles. The best agreement with the ammonia results is found for the HCO$^+$ (4--3) transitions, but the latter is still strongly blended with emission from associated outflows. This outflow signature is far less prominent in the THz ammonia lines, which confirms it as a powerful probe of infall in molecular clumps. Infall rates in the range from 0.3 to 16~$10^{-3}\,M_\odot/$yr were derived with a tentative correlation with the virial parameters of the clumps. The new observations show that infall on clump scales is ubiquitous through a wide range of evolutionary stages, from $L/M$ covering about ten to several hundreds.

INTERFEROMETRIC OBSERVATIONS OF HIGH-MASS STAR-FORMING CLUMPS WITH UNUSUAL N2H+/HCO+LINE RATIOS

The Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey has detected high-mass star-forming clumps with anomalous N$_2$H$^+$/HCO$^+$(1-0) integrated intensity ratios that are either unusually high ("N$_2$H$^+$ rich") or unusually low ("N$_2$H$^+$ poor"). With 3 mm observations from the Australia Telescope Compact Array (ATCA), we imaged two N$_2$H$^+$ rich clumps, G333.234-00.061 and G345.144-00.216, and two N$_2$H$^+$ poor clumps, G351.409+00.567 and G353.229+00.672. In these clumps, the N$_2$H$^+$ rich anomalies arise from extreme self-absorption of the HCO$^+$ line. G333.234-00.061 contains two of the most massive protostellar cores known with diameters of less than 0.1 pc, separated by a projected distance of only 0.12 pc. Unexpectedly, the higher mass core appears to be at an earlier evolutionary stage than the lower mass core, which may suggest that two different epochs of high-mass star formation can occur in close proximity. Through careful analysis of the ATCA observations and MALT90 clumps (including the G333, NGC 6334, and NGC 6357 star formation regions), we find that N$_2$H$^+$ poor anomalies arise at clump-scales and are caused by lower relative abundances of N$_2$H$^+$ due to the distinct chemistry of H II regions or photodissociation regions.

MALT90: tracing the chemistry and kinematics of molecular clumps within the central molecular zone

AbstractThe MALT90 survey targets more than 2000 high-mass star-forming clumps in the Galactic plane obtaining small maps around each of them, in 16 molecular lines at 90 GHz. By observing several thousand high-mass star forming clumps MALT90 aims to characterize their global chemical and physical evolution. Here we summarize the survey parameters and show examples of the MALT90 data toward three clumps in the central molecular zone.

MALT90: The Millimetre Astronomy Legacy Team 90 GHz Survey

AbstractThe Millimetre Astronomy Legacy Team 90 GHz (MALT90) survey aims to characterise the physical and chemical evolution of high-mass star-forming clumps. Exploiting the unique broad frequency range and on-the-fly mapping capabilities of the Australia Telescope National Facility Mopra 22 m single-dish telescope1, MALT90 has obtained 3′ × 3′ maps towards ~2 000 dense molecular clumps identified in the ATLASGAL 870 μm Galactic plane survey. The clumps were selected to host the early stages of high-mass star formation and to span the complete range in their evolutionary states (from prestellar, to protostellar, and on to $\mathrm{H\,{\scriptstyle {II}}}$ regions and photodissociation regions). Because MALT90 mapped 16 lines simultaneously with excellent spatial (38 arcsec) and spectral (0.11 km s−1) resolution, the data reveal a wealth of information about the clumps’ morphologies, chemistry, and kinematics. In this paper we outline the survey strategy, observing mode, data reduction procedure, and highlight some early science results. All MALT90 raw and processed data products are available to the community. With its unprecedented large sample of clumps, MALT90 is the largest survey of its type ever conducted and an excellent resource for identifying interesting candidates for high-resolution studies with ALMA.