Clump Mass Function Research Articles

Star formation towards the Scutum tangent region and the effects of Galactic environment

By positional matching to the catalogue of Galactic Ring Survey molecular clouds, we have derived distances to 793 Bolocam Galactic Plane Survey (BGPS) sources out of a possible 806 located within the region defined by Galactic longitudes l = 28.5 degr to 31.5 degr and latitudes |b| < 1 degr. This section of the Galactic Plane contains several major features of Galactic structure at different distances, mainly mid-arm sections of the Perseus and Sagittarius spiral arms and the tangent of the Scutum-Centarus arm, which is coincident with the end of the Galactic Long Bar. By utilising the catalogued cloud distances plus new kinematic distance determinations, we are able to separate the dense BGPS clumps into these three main line-of-sight components to look for variations in star-formation properties that might be related to the different Galactic environments. We find no evidence of any difference in either the clump mass function or the average clump formation efficiency (CFE) between these components that might be attributed to environmental effects on scales comparable to Galactic-structure features. Despite having a very high star-formation rate, and containing at least one cloud with a very high CFE, the star formation associated with the Scutum-Centarus tangent does not appear to be in any way abnormal or different to that in the other two spiral-arm sections. Large variations in the CFE are found on the scale of individual clouds, however, which may be due to local triggering agents as opposed to the large-scale Galactic structure.

Search for starless clumps in the ATLASGAL survey

In this study, we present an unbiased sample of the earliest stages of massive star formation across 20 square-degree of the sky. Within the region 10deg < l < 20deg and |b| < 1deg, we search the ATLASGAL survey at 870 micron for dense gas condensations. These clumps are carefully examined for indications of ongoing star formation using YSOs from the GLIMPSE source catalog as well as sources in the 24 micron MIPSGAL images, to search for starless clumps. We calculate the column densities as well as the kinematic distances and masses for sources where the v_lsr is known from spectroscopic observations. Within the given region, we identify 210 starless clumps with peak column densities > 1 x 10e23 cm^(-2). In particular, we identify potential starless clumps on the other side of the Galaxy. The sizes of the clumps range between 0.1 pc and 3 pc with masses between a few tens of solar masses up to several ten thousands of solar masses. Most of them may form massive stars, but in the 20 square-degree we only find 14 regions massive enough to form stars more massive than 20 solar masses and 3 regions with the potential to form stars more massive than 40 40 solar masses. The slope of the high-mass tail of the clump mass function for clumps on the near side of the Galaxy is 2.2 and, therefore, Salpeter-like. We estimate the lifetime of the most massive starless clumps to be 60000 yr. The sample offers a uniform selection of starless clumps. In the large area surveyed, we only find a few potential precursors of stars in the excess of 40 solar masses. It appears that the lifetime of these clumps is somewhat shorter than their free-fall times, although both values agree within the errors. In addition, these are ideal objects for detailed studies and follow-up observations.

The Carina Flare

$^{13}$CO(J=2--1) and C$^{18}$O(J=2--1) observations of the molecular cloud G285.90+4.53 (Cloud~16) in the Carina Flare supershell (GSH287+04-17) with the APEX telescope are presented. With an algorithm DENDROFIND we identify 51 fragments and compute their sizes and masses. We discuss their mass spectrum and interpret it as being the result of the shell fragmentation process described by the pressure assisted gravitational instability - PAGI. We conclude that the explanation of the clump mass function needs a combination of gravity with pressure external to the shell.

FEEDBACK-REGULATED STAR FORMATION: IMPLICATIONS FOR THE KENNICUTT-SCHMIDT LAW

We derive a metallicity dependent relation between the surface density of the star formation rate (Sigma_{SFR}) and the gas surface density (Sigma_{g}) in a feedback regulated model of star formation in galactic disks. In this model, star formation occurs in gravitationally bound protocluster clumps embedded in larger giant molecular clouds with the protocluster clump mass function following a power law function with a slope of -2. Metallicity dependent feedback is generated by the winds of OB stars (M > 5 Msol) that form in the clumps. The quenching of star formation in clumps of decreasing metallicity occurs at later epochs due to weaker wind luminosities, thus resulting in higher final star formation efficiencies (SFE_{exp}). By combining SFE_{exp} with the timescales on which gas expulsion occurs, we derive the metallicity dependent star formation rate per unit time in this model as a function of Sigma_{g}. This is combined with the molecular gas fraction in order to derive the global dependence of Sigma_{SFR} on Sigma_{g}. The model reproduces very well the observed star formation laws extending from low gas surface densities up to the starburst regime. Furthermore, our results show a dependence of $\Sigma_{SFR}$ on metallicity over the entire range of gas surface densities in contrast to other models, and can also explain part of the scatter in the observations. We provide a tabulated form of the star formation laws that can be easily incorporated into numerical simulations or semi-analytical models of galaxy formation and evolution.

High-mass star formation at high luminosities: W31 at &gt;106 L⊙

Context: High-mass star formation has been a very active field over the last decade, however, most studies targeted regions of luminosities between 10^4 and 10^5 L_sun. Methods: We selected the W31 star-forming complex with a total luminosity of ~6x10^6 L_sun for a multi-wavelength spectral line and continuum study covering wavelengths from the near- and mid-infrared via (sub)mm wavelength observations to radio data in the cm regime. Results: While the overall structure of the multi-wavelength continuum data resembles each other well, there are several intriguing differences. The 24mum emission stemming largely from small dust grains follows tightly the spatial structure of the cm emission tracing the ionized free-free emission. Hence warm dust resides in regions that are spatially associated with the ionized hot gas (~10^4 K) of the HII regions. Furthermore, we find several evolutionary stages within the same complexes, ranging from infrared-observable clusters, via deeply embedded regions associated with active star formation traced by 24\,$\mu$m and cm emission, to at least one high-mass gas clump devoid of any such signature. The 13CO(2-1) and C18O(2-1) spectral line observations reveal a large kinematic breadth in the entire region with a total velocity range of approximately 90 km/s. While the average virial mass ratio for W31 is close to unity, the line width analysis indicates large-scale evolutionary differences between the southern and northern sub-regions (G10.2-0.3 and G10.3-0.1) of the whole W31 complex. The clump mass function - tracing cluster scales and not scales of individual stars - derived from the 875mum continuum data has a slope of 1.5+-0.3, consistent with previous cloud mass functions.

Flow-driven cloud formation and fragmentation: results from Eulerian and Lagrangian simulations

The fragmentation of shocked flows in a thermally bistable medium provides a natural mechanism to form turbulent cold clouds as precursors to molecular clouds. Yet because of the large density and temperature differences and the range of dynamical scales involved, following this process with numerical simulations is challenging. We compare two-dimensional simulations of flow-driven cloud formation without self-gravity, using the Lagrangian Smoothed Particle Hydrodynamics (SPH) code VINE and the Eulerian grid code Proteus. Results are qualitatively similar for both methods, yet the variable spatial resolution of the SPH method leads to smaller fragments and thinner filaments, rendering the overall morphologies different. Thermal and hydro-dynamical instabilities lead to rapid cooling and fragmentation into cold clumps with temperatures below 300K. For clumps more massive than 1 Msun/pc, the clump mass function has an average slope of -0.8. The internal velocity dispersion of the clumps is nearly an order of magnitude smaller than their relative motion, rendering it subsonic with respect to the internal sound speed of the clumps, but supersonic as seen by an external observer. For the SPH simulations most of the cold gas resides at temperatures below 100K, while the grid-based models show an additional, substantial component between 100 and 300K. Independently of the numerical method our models confirm that converging flows of warm neutral gas fragment rapidly and form high-density, low-temperature clumps as possible seeds for star formation.

EVOLUTION OF A NUCLEAR GAS DISK AND GAS SUPPLY TO THE GALACTIC CENTER. I. MIGRATION OF MASSIVE GAS CLUMPS

We investigate the gas supply process from a nuclear gas disk to the Galactic center by performing two-dimensional hydrodynamic simulations including self-gravity and radiative cooling. We assume that the gas is supplied to the nuclear gas disk by the host galaxy disk. We show that many gas clumps are formed in the nuclear gas disk after the nuclear gas disk becomes massive enough to become gravitationally unstable. Coalescence of the gas clumps leads to the formation of more massive gas clumps. Gravitational encounters between these massive gas clumps induce angular momentum transfer between the gas clumps. Massive gas clumps that lose their angular momentum due to gravitational encounters move to the galactic center. In this way, the gas supply to the galactic center is realized. We also perform a numerical simulation with twice the spatial resolution to check numerical resolution effects. Time evolution of azimuthally averaged surface density distributions in both models agree well with each other. The high-mass part of the clump mass function in the low-resolution model is in good agreement with that of the high-resolution model, and lower mass clumps are well resolved in the high-resolution model. Therefore, we conclude that massive gas clumps play a major role in the angular momentum transfer. The gas supply rate to the center is about half of the gas supply rate from the host galaxy when the whole region of the nuclear gas disk becomes self-gravitationally unstable. More massive gas clumps are formed and migrate to the Galactic center in the models of the higher gas supply rate from the host galaxy.

CSO BOLOCAM 1.1 mm CONTINUUM MAPPING OF THE BRAID NEBULA STAR FORMATION REGION IN CYGNUS OB7

We present a 1.1 mm map of the Braid Nebula star formation region in Cygnus OB7 taken using Bolocam on the Caltech Submillimeter Observatory. Within the 1 deg2 covered by the map, we have detected 55 cold dust clumps all of which are new detections. A number of these clumps are coincident with IRAS point sources although the majority are not. Some of the previously studied optical/near-IR sources are detected at 1.1 mm. We estimate total dust/gas masses for the 55 clumps together with peak visual extinctions. We conclude that over the whole region, approximately 20% of the clumps are associated with IRAS sources suggesting that these are protostellar objects. The remaining 80% are classed as starless clumps. In addition, both FU Orionis (FUor) like objects in the field, the Braid Star and HH 381 IRS, are associated with strong millimeter emission. This implies that FUor eruptions can occur at very early stages of pre-main-sequence life. Finally, we determine that the cumulative clump mass function for the region is very similar to that found in both the Perseus and ρ Ophiuchus star-forming regions.

The fragmentation of expanding shells - III. Oligarchic accretion and the mass spectrum of fragments

We use SPH simulations to investigate the gravitational fragmentation of expanding shells through the linear and non--linear regimes. The results are analysed using spherical harmonic decomposition to capture the initiation of structure during the linear regime; the potential-based method of Smith et al. (2009) to follow the development of clumps in the mildly non-linear regime; and sink particles to capture the properties of the final bound objects during the highly non-linear regime. In the early, mildly non--linear phase of fragmentation, we find that the clump mass function still agrees quite well with the mass function predicted by the analytic model. However, the sink mass function is quite different, in the sense of being skewed towards high-mass objects. This is because, once the growth of a condensation becomes non-linear, it tends to be growing non-competitively from its own essentially separate reservoir; we call this Oligarchic Accretion.

Mass functions in fractal clouds: the role of cloud structure in the stellar initial mass function

The possibility that the stellar initial mass function (IMF) arises mostly from cloud structure is investigated with fractal Brownian motion (fBm) clouds that have power-law power spectra. An fBm cloud with a realistic projected power spectrum slope of $\beta=2.8$ is found to have a mass function for clumps exceeding a threshold density that is a power-law with a slope of $\alpha=2.35$, the same as in the Salpeter IMF. Any hierarchically structured cloud has a clump mass function with about the same slope. This result implies that turbulent interstellar clouds produce dense substructure with the observed pre-stellar core mass function built in from the start. Details of the clump formation processes are not critical. The conversion of clumps into stars involves a second step. A one-to-one correspondence between clump mass and star mass is not necessary to convert the clump mass spectrum into an IMF with the same power-law slope. As long as clumps have an internal stellar IMF from sub-fragmentation, protostellar accretion, coalescence and other processes, and the characteristic mass for this internal IMF scales with the clump mass, then the IMF slope above the minimum characteristic mass will equal the clump mass slope. A detailed review of IMF models illustrates the prominence of cloud structure as a major component in a wide class of theories. Tests are proposed to determine the relative importance of cloud structure and competitive accretion in the IMF.

MEASURING THE CLUMP MASS FUNCTION IN THE AGE OF SCUBA2,HERSCHEL, AND ALMA

We use simulated images of star-forming regions to explore the effects of various image acquisition techniques on the derived clump mass function. In particular, we focus on the effects of finite image angular resolution, the presence of noise, and spatial filtering. We find that, even when the image has been so heavily degraded with added noise and lowered angular resolution that the clumps it contains clearly no longer correspond to pre-stellar cores, still the clump mass function is typically consistent with the stellar initial mass function within their mutual uncertainties. We explain this result by suggesting that noise, source blending, and spatial filtering all randomly perturb the clump masses, biasing the mass function toward a lognormal form whose high-mass end mimics a Salpeter power law. We argue that this is a consequence of the central limit theorem and that it strongly limits our ability to accurately measure the true mass function of the clumps. We support this conclusion by showing that the characteristic mass scale of the clump mass function, represented by the ``break mass'', scales as a simple function of the angular resolution of the image from which the clump mass function is derived. This strongly constrains our ability to use the clump mass function to derive a star formation efficiency. We discuss the potential and limitations of the current and next generation of instruments for measuring the clump mass function.

Analogues of Cores and Stars in Simulated Molecular Clouds

AbstractWe use images derived from collapsing, turbulent molecular cloud simulations without sinks to explore the effects of finite image angular resolution and noise on the derived clump mass function. These effects randomly perturb the clump masses, producing a lognormal clump mass function with a Salpeter-like high mass end. We show that the characteristic break mass of the simulated clump mass functions changes with the angular resolution of the images in a way that is entirely consistent with the observations. We also present some cautionary tales regarding sink particles and highlight the need to ensure that sinks actually correspond to distinct collapsing objects. We test several popular numerical sink criteria in the literature and compare to converged, non-sink results.

A Guide to Comparisons of Star Formation Simulations with Observations

AbstractWe review an approach to observation-theory comparisons we call “Taste-Testing.” In this approach, synthetic observations are made of numerical simulations, and then both real and synthetic observations are “tasted” (compared) using a variety of statistical tests. We first lay out arguments for bringing theory to observational space rather than observations to theory space. Next, we explain that generating synthetic observations is only a step along the way to the quantitative, statistical, taste tests that offer the most insight. We offer a set of examples focused on polarimetry, scattering and emission by dust, and spectral-line mapping in star-forming regions. We conclude with a discussion of the connection between statistical tests used to date and the physics we seek to understand. In particular, we suggest that the “lognormal” nature of molecular clouds can be created by the interaction of many random processes, as can the lognormal nature of the IMF, so that the fact that both the “Clump Mass Function” (CMF) and IMF appear lognormal does not necessarily imply a direct relationship between them.

Fragmentation in molecular clouds and its connection to the IMF

We present an analysis of star-forming gas cores in an SPH simulation of a Giant Molecular Cloud. We identify cores using their deep potential wells. This yields a smoother distribution with clearer boundaries than density. Additionally, this gives an indication of future collapse, as bound potential cores (p-cores) represent the earliest stages of fragmentation in molecular clouds. We find that the mass function of the p-cores resembles the stellar IMF and the observed clump mass function, although p-core masses (~0.7 Msol) are smaller than typical density clumps. The bound p-cores are generally subsonic, have internal substructure, and are only quasi-spherical. We see no evidence of massive bound cores supported by turbulence. We trace the evolution of the p-cores forward in time, and investigate the connection between the original p-core mass and the stellar mass that formed from it. We find that there is a poor correlation, with considerable scatter suggesting accretion onto the core is dependent on more factors than just the initial core mass. During the accretion process the p-cores accrete from beyond the region first bound, highlighting the importance of the core environment to its subsequent evolution.

The structure of molecular clouds and the universality of the clump mass function

Using an SPH simulation of a star-forming region in a molecular cloud, we show that the emergence of a clump mass function (CMF) resembling the stellar initial mass function (IMF) is a ubiquitous feature of molecular cloud structure, but caution against its over-interpretation. We employ three different techniques to extract the clumps in this study. In the first two, we interpolate the SPH particle data to 2 and 3 dimensional grids before performing the clump-find, using position-position (PP) and position-position-velocity (PPV) information respectively. In the last technique, the clump-finding is performed on the SPH data directly, making use of the full 3 dimensional position information. Although the CMF is typically similar to that observed in regions of nearby star formation, the individual clumps and their masses are found to be unreliable since they depend strongly on the parameters and the method of the clump-finding. In particular we find that the resolution and orientation of the data make a significant difference to the resulting properties of the identified clumps in the PP and PPV cases. We conclude that making comparisons between a CMF and the stellar IMF should be done with caution, since the definition of a clump boundary, and hence the number of clumps and their properties, are arbitrary in the extraction method. This is especially true if molecular clouds are truly scale free.

Clump lifetimes and the initial mass function

Recent studies of dense clumps/cores in a number of regions of low-mass star formation have shown that the mass distribution of these clumps closely resembles the initial mass function (IMF) of field stars. One possible interpretation of these observations is that we are witnessing the fragmentation of the clouds into the IMF, and the observed clumps are bound pre-stellar cores. In this paper, we highlight a potential difficulty in this interpretation, namely that clumps of varying mass are likely to have systematically varying lifetimes. This ‘time-scale’ problem can effectively destroy the similarity between the clump and stellar mass functions, such that a stellar-like clump mass function results in a much steeper stellar IMF. We also discuss some ways in which this problem may be avoided.

The star-forming content of the W3 giant molecular cloud

Abstract We have surveyed a ∼0.9 square degree area of the W3 giant molecular cloud (GMC) and star-forming region in the 850-μm continuum, using the Submillimetre Common-User Bolometer Array on the James Clerk Maxwell Telescope. A complete sample of 316 dense clumps were detected with a mass range from around 13 to 2500 M⊙. Part of the W3 GMC is subject to an interaction with the H ii region and fast stellar winds generated by the nearby W4 OB association. We find that the fraction of total gas mass in dense, 850-μm traced structures is significantly altered by this interaction, being around 5–13 per cent in the undisturbed cloud but ∼25–37 per cent in the feedback-affected region. The mass distribution in the detected clump sample depends somewhat on assumptions of dust temperature and is not a simple, single power law but contains significant structure at intermediate masses. This structure is likely to be due to crowding of sources near or below the spatial resolution of the observations. There is little evidence of any difference between the index of the high-mass end of the clump mass function in the compressed region and in the unaffected cloud. The consequences of these results are discussed in terms of current models of triggered star formation.

Massive Clumps in the NGC 6334 Star‐forming Region

We report observations of dust continuum emission at 1.2 mm toward the star forming region NGC 6334 made with the SEST SIMBA bolometer array. The observations cover an area of $\sim 2$ square degrees with approximately uniform noise. We detected 181 clumps spanning almost three orders of magnitude in mass (3\Msun$-6\times10^3$ \Msun) and with sizes in the range 0.1--1.0 pc. We find that the clump mass function $dN/d\log M$ is well fit with a power law of the mass with exponent -0.6 (or equivalently $dN/dM \propto M^{-1.6}$). The derived exponent is similar to those obtained from molecular line emission surveys and is significantly different from that of the stellar initial mass function. We investigated changes in the mass spectrum by changing the assumptions on the temperature distribution of the clumps and on the contribution of free-free emission to the 1.2 mm emission, and found little changes on the exponent. The Cumulative Mass Distribution Function is also analyzed giving consistent results in a mass range excluding the high-mass end where a power-law fit is no longer valid. The masses and sizes of the clumps observed in NGC 6334 indicate that they are not direct progenitors of stars and that the process of fragmentation determines the distribution of masses later on or occurs at smaller spatial scales. The spatial distribution of the clumps in NGC 6334 reveals clustering which is strikingly similar to that exhibited by young stars in other star forming regions. A power law fit to the surface density of companions gives $\Sigma\propto \theta^{-0.62}$.

Large‐Area Mapping at 850 μm. V. Analysis of the Clump Distribution in the Orion A South Molecular Cloud

We present results from a 2300 arcmin^2 survey of the Orion A molecular cloud at 450 and 850 micron using the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope. The region mapped lies directly south of the OMC1 cloud core and includes OMC4, OMC5, HH1/2, HH34, and L1641N. We identify 71 independent clumps in the 850 micron map and compute size, flux, and degree of central concentration in each. Comparison with isothermal, pressure-confined, self-gravitating Bonnor-Ebert spheres implies that the clumps have internal temperatures T_d ~ 22 +/- K and surface pressures log (k^-1 P cm^-3 K) = 6.0 +/- 0.2. The clump masses span the range 0.3 - 22 Msun assuming a dust temperature T_d ~ 20 K and a dust emissivity kappa_850 = 0.02 cm^2 g^-1. The distribution of clump masses is well characterized by a power-law N(M) propto M^-alpha with alpha = 2.0 +/- 0.5 for M > 3.0 Msun, indicating a clump mass function steeper than the stellar Initial Mass Function. Significant incompleteness makes determination of the slope at lower masses difficult. A comparison of the submillimeter emission map with an H_2 2.122 micron survey of the same region is performed. Several new Class 0 sources are revealed and a correlation is found between both the column density and degree of concentration of the submillimeter sources and the likelihood of coincident H_2 shock emission.

High‐Mass Star Formation. III. The Functional Form of the Submillimeter Clump Mass Function

We investigate the mass function of cold, dusty clumps in 11 low- and high-mass star-forming regions. Using a homogeneous fitting technique, we analyze the shape of each region's clump mass function and examine the commonalities among them. We find that the submillimeter continuum clump mass function in low-mass star-forming regions is typically best fitted by a lognormal distribution, while that in high-mass star-forming regions is better fitted by a double power law. A single power-law clump mass distribution is ruled out in all cases. Fitting all of the regions with a double power law, we find that the mean power-law exponent at the high-mass end of each mass function is αhigh = -2.4 ± 0.1, consistent with the Salpeter result of α = -2.35. We find no region-to-region trend in αhigh with the mass scale of the clumps in a given region, as characterized by their median mass. Similarly, nonparametric tests show that the shape of the clump mass function does not change much from region to region, despite the obvious changes in the intrinsic mass scale. This result is consistent with the hypothesis that the clump mass distribution is determined by a highly stochastic process, such as turbulent fragmentation. It may also suggest that the data reduction and analysis techniques strongly affect the shape of the derived mass function.