Abstract
Dynamical studies of dense structures within molecular clouds often conclude that the most massive clumps contain too little kinetic energy for virial equilibrium, unless they are magnetized to an unexpected degree. This raises questions about how such a state might arise, and how it might persist long enough to represent the population of massive clumps. In an effort to reexamine the origins of this conclusion, we use ammonia line data from the Green Bank Ammonia Survey and Planck-calibrated dust emission data from Herschel to estimate the masses and kinetic and gravitational energies for dense clumps in the Gould Belt clouds. We show that several types of systematic error can enhance the appearance of low kinetic-to-gravitational energy ratios: insufficient removal of foreground and background material; ignoring the kinetic energy associated with velocity differences across a resolved cloud; and overcorrecting for stratification when evaluating the gravitational energy. Using an analysis designed to avoid these errors, we find that the most massive Gould Belt clumps harbor virial motions, rather than subvirial ones. As a by-product, we present a catalog of masses, energies, and virial energy ratios for 85 Gould Belt clumps.
Export citation and abstract BibTeX RIS
1. Introduction
A key parameter in the study of molecular clouds and their substructures is the virial ratio:
![Equation (1)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn1.gif)
Here, is the object's kinetic energy in its center-of-mass frame, and
is its self-gravitational energy. The value of α as a diagnostic tool arises from the fact that it compares prominent opposing terms in the virial theorem (Chandrasekhar & Fermi 1953; McKee & Zweibel 1992)—that is, in the competition of forces (expressed as energies) that cause inward or outward accelerations of an object's radius (expressed as the trace of its moment-of-inertia tensor). Because any significant imbalance leads to rapid change, there is good reason to expect that some chosen collection of dense interstellar structures is close to a state of equilibrium, at least in a statistical sense, especially if the structures live for at least a single crossing time.
17
Insofar as this is true, one can read the value of α as an indication of the other, less easily observed, forces or energetic terms. A state in which α > 1 suggests the importance of the kinetic surface term that represents confinement by external thermal or turbulent pressure, or by the ram pressure due rapid inflow or outflow (see Goldbaum et al. 2011), or due to colliding flows. For example, α = 1.683 in the critical state of an isothermal, unmagnetized, pressure-bounded sphere (Ebert 1955; Bonnor 1956). If α ≫ 1, then gravity is negligible in comparison to external pressure: an equilibrium object in this state is pressure confined.
In contrast, α < 1 indicates the importance of an additional positive term, corresponding to an outward force that opposes the combination of self-gravity and external pressure. Magnetic fields supply one such force, although for molecular clumps and cores, both Zeeman measurements (Crutcher 2012) and estimates based on the Davis–Chandraskhar–Fermi method (Myers & Basu 2021) indicate median mass-to-flux ratios about twice the critical value. This implies that the quasi-static portion of the magnetic force is rarely sufficient to fully offset gravity. The fluctuating portion of the magnetic energy is also limited in its impact, as it tends to be in equipartition with the turbulent portion of (McKee & Zweibel 1992; Federrath 2016). Another, often-overlooked outward force is the momentum injected by protostellar outflows or photoionized regions when star formation is especially active. These must be included as an inner surface term when their kinetic energies are not included in
. The importance of this term is evident in the fact that the kinetic energy in protostellar outflows can be comparable to
(e.g., Graves et al. 2010). However, their effect cannot overwhelm
in the virial theorem, for the simple reason that star-driven flows go on to stir turbulence within the medium (Matzner 2002, 2007; Nakamura & Li 2007). Similar arguments apply to stellar radiation forces (McKee & Zweibel 1992), which in any case are only significant in the presence of vigorous massive star formation (Krumholz & Matzner 2009; Murray et al. 2010; Raskutti et al. 2016; Jumper & Matzner 2018). Yet another often-overlooked term comes from the gravity of matter outside the clump boundary (Ballesteros-Paredes 2006); however, this should be small compared to
for clumps that are overdense and not bounded by tidal forces.
For these reasons it is difficult to envision a scenario in which any collection of interstellar structures would be strongly subvirial—that is, characterized by α ≪ 1.
It is very puzzling, therefore, that observational studies of dense substructures within molecular clouds often find that α is well below unity for the most massive of these objects (see, for instance, Kauffmann et al. 2013; Urquhart et al. 2014; Traficante et al. 2018a, 2018c). Although the selection of objects varies from one study to another (as does the specific correlation between estimates of α and mass), the substructures in question are all molecular clumps: objects intermediate in scale between molecular clouds and the compact cores from which individual star systems are born. If massive molecular clumps are truly subvirial, this has important implications for the initial conditions for star cluster formation.
Could the strongly subvirial appearance of massive molecular clumps in fact be an artifact of the way that observations are taken or interpreted? Traficante et al. (2018b) advance one reason that it might be. They point out that the data used to determine and
tend to weight different regions of a clump, due to the influence of a critical density on molecular line excitation, and that this may lead to a systematic offset in α.
Here, we explore another possibility: that choices involved in the method used to estimate α may themselves introduce systematic errors. Because it depends on an object's three-dimensional density, temperature, and velocity fields, α cannot be determined from projected data. One must construct a proxy, such as the virial parameter introduced by Bertoldi & McKee (1992, hereafter BM92) or the estimate presented by Singh et al. (2019, hereafter SMJ19). We find that at several steps of this process, common analysis choices have the cumulative effect of suppressing the derived value of α, especially for high-mass clumps.
For our exploration we employ NH3 data from the Green Bank Ammonia Survey (GAS; Friesen et al. 2017) and column densities derived from a new analysis of dust optical depth in the Herschel data (A. Singh & P.G. Martin, 2021, in preparation). These provide sensitive, uniform, well-calibrated, and well-resolved information for an investigation such as ours. As a by-product, we present a catalog of properties for 85 clumps in the Gould Belt, using these high-quality data. Our catalog overlaps previous virial analyses, based on a subset of the same data, conducted by Kirk et al. (2017a), Redaelli et al. (2017), Keown et al. (2017), Chen et al. (2019), and Kerr et al. (2019). However, we use a somewhat different algorithm to define clump boundaries; this allows us to focus on the influence of analysis choices on estimates of α.
In Section 2, we review methods for estimating α. We introduce the data for our study in Section 3, and present the details of our technique in Section 4. As a case study, we highlight the analysis of a single clump in Section 4.1. In Section 5, we present results for our full sample of Gould Belt clumps. Finally, in Section 6, we draw conclusions about the impact of biases on the apparent physical state of massive clumps.
2. Methods
The most widely used method for estimating α involves the virial parameter introduced by BM92:
![Equation (2)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn2.gif)
Here, Mcl and σcl are the mass and one-dimensional velocity dispersion, respectively, of the object under study (a clump, in our case), and Rcl is its effective radius, usually defined so the projected area of the clump is . We adopt that definition as well.
The value of αBM92 derives from the fact that the gravitational energy of an interstellar object is usually similar to that of a uniform sphere with the same mass and effective radius; hence, the correction factor a, defined by , is of order unity. BM92 consider the class of spheroidal clumps with power-law density profiles ρ ∝ r−k
in spheroidal radius coordinate r as an example. For these, a can be decomposed (a = a1
a2) into a stratification factor a1 = (1 − k/3)/(1 − k/2.5), and a geometric factor a2 whose angle average is close to unity except in the case of very prolate spheroids (see Figure 2 of BM92).
Along with the definition of a, Equations (1) and (2) imply
![Equation (3)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn3.gif)
showing that either αBM92, or the refined quantity αBM92/a, can be used to estimate α. An important detail is that σcl must be defined so that is a valid estimate for
we return to this point below.
As an alternative to αBM92 we will consider the method proposed by SMJ19, who reexamine the process of estimating α and suggest a procedure that is robust for non-spheroidal structures despite the effects of projection.
The SMJ19 method begins with how a clump is identified and extracted from projected data. One must start by defining a projected clump boundary, which encompasses a peak in the column density Σ. The projected clump boundary should also contain reliable molecular line data from which to obtain the line-of-sight radial velocity vz , the one-dimensional thermal velocity dispersion σth, and the total line-of-sight velocity dispersion σz , from which the nonthermal portion σNT,z can be derived.
The next step is to extract the cloud or clump column density Σcl from any external material projected within the cloud boundary. Foreground and background removal is especially important when Σ is derived from submillimeter dust emission; note that spatially filtered observations (such as chopped or interferometric ones) accomplish this in an approximate way. SMJ19 demonstrate that an approximation based on the Abel transform is more successful than either keeping all material within the cloud boundary, or using simple interpolation to clip a background level. Abel reconstruction uses the fact that there is an exact relationship between any axisymmetric density distribution and its projection. Applied to the column density, it provides an estimate for the component interior to a three-dimensional surface whose projection is the two-dimensional clump boundary. The removed component, Σenv = Σ − Σcl, is always lower toward the core of a clump than toward its edge, thanks to a projection effect that is analogous to limb brightening.
As only one component of the cloud's kinetic energy is visible in projection, the quantity
![Equation (4)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn4.gif)
where x and y are sky coordinates at the cloud distance, provides an estimate for that is unbiased, in the sense that its average over viewing angles (denoted
) is exact:
![Equation (5)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn5.gif)
We note that BM92 adopt the same quantity, defining in their Appendix C. In Equation (4),
![](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2ueqn1.gif)
is the square of the effective one-dimensional velocity dispersion along each line of sight, and
![Equation (6)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn6.gif)
is the contribution from resolved variations in the line-of-sight velocity vz across the cloud, which we will refer to as bulk kinetic energy. Note that this may include a portion due to rotation. The center-of-mass velocity vCM,z is computed from Mcl vCM,z = ∫Σcl(x, y)vz (x, y)dx dy.
It is worth noting that will be resolution-dependent, in practice, because σcl,z
(x, y) can only be defined at the resolution of a given experiment. In the limit that a cloud is too small to be resolved, it would be described by the single centroid velocity vz
= vCM,z
, and thus have
. However, its observed kinetic energy should still be captured well by Equation (4) because its line width σcl,z
will include these unresolved velocity gradients.
For the denominator of α, SMJ19 define a quantity , derived from Σcl, with the desirable property that
![Equation (7)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn7.gif)
SMJ19 show that can be computed by collapsing the clump mass profile to a sheet in the plane of the sky, obtaining the sheet's gravitational self-energy, and correcting the result by a factor 2/π. By using information from the resolved column density map, this procedure avoids the need to choose Rcl and estimate a. It therefore provides a means to calibrate the characteristic value of a for an ensemble of clouds.
With these definitions for and
, SMJ19's quantity
![Equation (8)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn8.gif)
is a valid estimate for α, in the sense that both the numerator and denominator are exact when averaged over viewing angles.
It is important to note that αSMJ19 and αBM92/a are equivalent if evaluated under the following conditions: (1) they are derived from the same model of the clump column density profile Σcl (and furthermore, to give valid estimates of α, this must have been been cleaned of foreground and background contamination); (2) the quantity σcl is defined so that , thus incorporating
and (3) the BM92 correction factor a is evaluated in a way that properly reflects the clump profile.
As we shall see, other choices for Mcl, σcl, and a tend to introduce biases when αBM92/a (sometimes called Mvir/Mcl) is used as an estimate for α.
3. Data
3.1. GAS Molecular Line Data
All of the thermal and kinetic information we use in the calculation of derives from ammonia line data from GAS (Friesen et al. 2017). GAS observations of the NH3 (1,1) and (2,2) inversion transitions, carried out with the K-Band Focal Plane Array of the Green Bank Telescope, achieve sufficient sensitivity (∼0.1 K median noise), spatial resolution (32'', or 0.047 pc at 300 pc), and frequency resolution (5.7 kHz/23.7 GHz, or 0.07 km s−1) to map and resolve gas temperature, centroid velocity, and velocity dispersion across a wide range of dense core, clump, and filamentary structures. These parameters are obtained from a single-component fit, a point we return to in Section 5.1.
Our source data (GAS Data Release 2: J. Pineda et al., 2021, in preparation) incorporate the latest improvements to the GAS analysis pipeline relative to Data Release 1, including improved sensitivity arising from changes to the multicomponent fitting pipeline, as well as an expanded data set of star-forming regions. The regions examined are listed in Table 1.
Table 1. Cloud Regions and Adopted Distances
Cloud | Herschel Name | Adopted Distance |
---|---|---|
[pc] | ||
B1 | Perseus | 301 a |
L1448 | Perseus | 288 a |
L1451 | Perseus | 279 a |
L1455 | Perseus | 235 b |
NGC 1333 | Perseus | 293 c |
Perseus | Perseus | 235 b |
IC 348 | Perseus | 320 c |
B18 | Taurus | 126.6 d |
HC 2 | Taurus | 138.6 d |
IC 5146 | IC 5146 | 831 e |
Cepheus L1228 | Cepheus | 346 f |
Cepheus L1251 | Cepheus | 346 f |
CrA west | Corona Australis | 154 e |
CrA east | Corona Australis | 154 e |
L1688 | Ophiuchus | 138.4 g |
L1689 | Ophiuchus | 144.2 g |
Serpens MWC 279 | Serpens | 437 h |
Orion A | Orion A | 388 i |
Orion A S | Orion A | 428 i |
Orion B NGC 2023-2024 | Orion B | 420 i |
Orion B NGC 2068-2071 | Orion B | 388 i |
Notes.
a Zucker et al. (2018). b Hirota et al. (2008). c Ortiz-León et al. (2018a). d Galli et al. (2018). e Dzib et al. (2018). f Yan et al. (2019). g Ortiz-León et al. (2018b). h Ortiz-León et al. (2017). i Kounkel et al. (2017).Download table as: ASCIITypeset image
3.2. Herschel-derived Column Densities
To estimate the mass column density we employ dust optical depth maps derived by fitting a spectral energy distribution (SED) to continuum data from the Hershel Space Observatory at 160, 250, 350, and 500 μm. We use the results of an improved analysis to be described in an upcoming work (A. Singh & P.G. Martin, 2021, in preparation). A zero-point correction was applied to the Herschel intensity maps at each wavelength by correlating them with intensity models created from Planck dust models (Planck Collaboration XI 2014). Intensity maps were then fitted by a modified blackbody to estimate the dust temperature and optical depth using the dust emissivity index β determined in each pixel from the Planck dust models. In the Singh & Martin, (2021, in preparation) SED fitting pipeline, various cross comparisons are used to optimize the determination of data and model uncertainties, thereby improving the robustness of the final maps. We adopt a constant dust-plus-gas opacity of κν,0 = 0.1 cm2 g−1 at 1 THz (Hildebrand 1983, with a gas-to-dust mass ratio of 100) to determine Σ.
4. Implementation
We identify objects for our analysis as identifiable peaks in NH3 emission and total column density. We draw projected clump and region boundaries around these according to the following procedure.
To isolate structures for which we have complete and well-resolved data, we start by creating a version of the NH3 column density map that we convolve with a bounded parabolic (Epanechnikov) kernel of 88ʺ radius. Contours of this smoothed map are candidates for the clump boundary, from which we choose the largest for which two criteria are met: First, we must have complete data coverage to determine Σ, vz
, and σz
, and nearly complete coverage for Tk
, within the boundary. This requires that NH3 (1,1) and (2,2) transitions are both well detected. (We fill any gaps in the temperature data using linear interpolation, which adds a negligible uncertainty to .) Second, we require that the inferred mass density
implies nH
> 3.6 × 103 cm−3, the critical density of the NH3 (1,1) transition, which exceeds that of the (2,2) transition. Because NH3 has a reasonably consistent abundance at these densities (e.g., Redaelli et al. 2017) and is not affected by freeze-out or significant line optical depth, this choice ensures that our kinematic information can be used to calculate
in a reasonably unbiased way. Nevertheless, abundance variations (e.g., Crapsi et al. 2007) do inevitably inject some uncertainty.
Once a trial boundary has been identified, it is modified if necessary to exclude regions in which vz falls more than 4 standard deviations outside its overall distribution for the clump. Velocity outliers probably represent errors in the line fitting process, or possibly NH3 emission from unrelated structures along the line of sight, and so it is important to exclude the spurious additional kinetic energy they would imply. We find that this approach affects only a small fraction of the clump area and mass (a few pixels at most). Figure 1 provides one example, for clump L1251-1 in Cepheus.
Figure 1. NH3 line-of-sight velocity data for Clump 1 in Cepheus L1251, displaying a contaminating velocity outlier and our method for removing it. The left panel shows vz
relative to the center-of-mass velocity within the original clump boundary. A single pixel of this map lies outside 4 standard deviations from the clump velocity distribution, as seen on the right panel. We revise the cloud boundary to excise this pixel, and recalculate quantities like vCM,z
before deriving . Variations of vz
across the map contribute a bulk kinetic energy
that can be important: in this case, αSMJ19 = 2.05, but evaluates to 0.65 if
is ignored.
Download figure:
Standard image High-resolution imageAs the figure shows, it is possible for our procedure to produce a clump boundary that includes a hole, indicating where NH3 data is lacking. L1251-1 has by far the largest such hole,
18
covering 2% of its area. A number of other clumps have holes that span a few pixels; we indicate these cases with asterisks in Table 3. Because the lack of NH3 data tends to coincide with a region of low to zero Σcl, it makes no practical difference to or
whether we fill these holes with interpolated data. In computing radii we use the outer area of the mask; excluding the holes would only slightly lessen Rcl.
For every clump, we must also identify the boundary of an enclosing region for the purposes of removing dust emission from foreground and background material, thereby creating a map with compact support for Abel-transform analysis. While the exact choice of region boundary is not significant, it should be separated well enough from the clump boundary to contain the envelope physically associated with the clump, yet also close enough to the clump to sample a similar column density of background and foreground material (or at least, the component of this material on long spatial scales).
To accomplish this, we generate a trial map and choose one of its contours to be the region boundary. We create a version of the NH3 column density smoothed with a Gaussian kernel of width Rcl/8. We then multiply this by a vignetting function that is equal to the clump mask (unity within the clump boundary, zero outside), convolved with the function , where rk
is the distance from the center of the kernel as measured in the plane of the sky. From this trial map, we choose the lowest contour that comes within Rcl/4 of the clump boundary. An example region boundary is displayed in Figure 2.
Figure 2. An example of the clump extraction process, applied to Clump 3 of the Perseus B1 cloud discussed in Section 4.1. Panels (a)–(c) show the molecular data for line-of-sight velocity, ammonia column density, and kinetic temperature, respectively, while (d)–(g), in turn, illustrate the total column Σ, removal of a foreground/background column Σbg (by interpolation), and the clump envelope contribution Σenv (by Abel reconstruction), leaving Σcl. Note that the region boundary and clump boundary are identified in panels (a) and (b), respectively, as described in Section 4. Panels (f) and (g) are shown at higher magnification; the bar at the bottom of each panel indicates 1 pc for scale.
Download figure:
Standard image High-resolution imageGiven the clump and region contours, we use the SMJ19 prescription to obtain Σcl in two steps. First we subtract from Σ a bi-cubic interpolation of its value from the region boundary; this is meant to extract an estimate of the emission from foreground/background dust. We then apply the Abel transform, in the manner discussed by SMJ19, to subtract the emission from the clump's envelope, leaving Σcl. The lower panels of Figure 2 provides an illustration of the decomposition, for the case of Clump 3 in Perseus B1.
From the clump column density profile Σcl we derive an estimate for self-gravitational energy according to Equation (4) of SMJ19; we ensure that
is accurate to within a fraction of 1% using numerical refinement and a correction for discretization errors,which we describe in the Appendix.
In addition to our fiducial procedure, we consider several alternatives to identify what affects the trend of α with Mcl. One such choice involves employing αBM92/a, but assuming a value for a (such as unity). Others involve evaluating Σcl differently. For instance, one might replace the Abel transform with a simple interpolation, leaving a different clump column profile, or one might adopt Σcl = Σ, forgoing any foreground and background removal. We shall also consider the effect of neglecting when evaluating σcl and
, to demonstrate that the velocity dispersion derived from individual beams does not capture the entire kinetic energy along the line of sight.
Note that our definition of Rcl, which derives from the projected area of NH3 emission, differs from definitions involving fitted profiles or moments of the column density distribution. This should be kept in mind when comparing radii among catalogs. Moreover, because radius enters into the evaluation of αBM92, the method used to determine Rcl can introduce a new and potentially significant source of systematic error, which we do not attempt to evaluate here.
4.1. Case Studies: Cepheus L1251-1 and Perseus B1-3
As our primary case study we consider the clump shown in Figure 1, Clump 1 within L1251 in Cepheus. Applying our fiducial analysis we obtain erg and
erg, of which
is 71%. These values imply αSMJ19 = 2.05.
In Table 2 we provide additional parameters for this particular clump, comparing our fiducial (SMJ19) evaluation of α against the BM92 analysis, and against alternative reconstructions of Σcl based on removing an interpolated column density level from the clump boundary ("Edge Interp."), or making no correction for foreground and background contamination. We also consider the effect of neglecting .
Table 2. Impact of Analysis Choices for Cepheus L1251-1
Quantity | Abel | Edge Interp. | No Bg./Fg. Corr. |
---|---|---|---|
Mcl [M⊙] | 56.2 | 43.2 | 82.6 |
Rcl [pc] | 0.32 | 0.32 | 0.32 |
σcl [km s−1] | 0.59 | 0.6 | 0.58 |
![]() | 5.68 × 1044 | 3.62 × 1044 | 1.08 × 1045 |
![]() | 5.83 × 1044 | 4.59 × 1044 | 8.2 × 1044 |
![]() | 3.99 × 1044 | 3.16 × 1044 | 5.55 × 1044 |
αSMJ19 | 2.05 | 2.54 | 1.52 |
αSMJ19 (No ![]() | 0.65 | 0.79 | 0.49 |
αBM92 | 2.3 | 3.05 | 1.49 |
αBM92 (No ![]() | 0.73 | 0.95 | 0.48 |
Note. Columns indicate method used to remove foreground and background emission from Σcl: our fiducial Abel-transform reconstruction, interpolation from the clump boundary, or no correction at all. Rows marked "No " indicate that resolved kinetic energy is omitted from
and σcl.
Download table as: ASCIITypeset image
Several aspects of the analysis affect the result. The method used to extract the clump's column density profile has a significant effect on its derived mass: simple interpolation from the clump boundary yields a 25% lower estimate, whereas making no correction for foreground and background emission attributes 49% more mass to the clump. Of the observed kinetic energy, constitutes 68% in this case. Equating
and
(or equivalently, αSMJ19 and αBM92/a) requires a = 1.12; notably, this is quite close to unity. Taken together, we see that overestimating the mass, excluding
, and adopting a value for a significantly above unity (e.g., a = a1 = 5/3, corresponding to k = 2, as chosen by Roman-Duval et al. 2010 and Csengeri et al. 2017) can lead one to underestimate α for this object by as much as a factor of 7.3.
How representative is Cepheus L1251-1? We note that Perseus B1-3, the clump depicted in Figure 2, yields an identical estimate for a (i.e., a = 1.12) but contains only 22% of its kinetic energy in bulk form. Our estimate for its virial ratio is αSMJ19 = 0.58, the lowest of our entire sample. In its case, neglecting to correct for foreground or background material, omitting , and adopting a = 5/3 would lead one to underestimate α by the more moderate factor of 2.7.
In the next section, we place these examples into context, and show that Cepheus L1251-1 and Perseus B1-3 are not outliers, but do bracket the range of massive clump properties.
5. Results: Are Massive Dense Clumps Truly Subvirial?
We list our findings for all 85 of our Gould Belt NH3 clumps in Table 3, and plot them for various choices of the analysis method in Figures 3 and 4.
Figure 3. Dark and light blue squares αBM92 with and without the contribution of , respectively; red and orange diamonds represent αSMJ19 with and without
, respectively. For αBM92 we set Rcl equal to the geometric mean of semimajor and semiminor axes of an ellipse fit to the clump boundary. Pink lines connect data for the same clump. The overplotted curves are quadratic fits to
for the points that share their line color, and are shown only to clarify trends among the four point families.
Download figure:
Standard image High-resolution imageTable 3. Clump Properties and Virial Ratios Derived from Abel Reconstruction
Clump | R.A. | Decl. | Mass | Radius | αSMJ19 | αSMJ19 | αBM92 | αBM92 |
---|---|---|---|---|---|---|---|---|
[M⊙] | [pc] | (No ![]() | (No ![]() | |||||
B1-1 | 03h33m36s | +31d18m49s | 2.97 | 0.11 | 8.49 | 8.23 | 8.82 | 8.55 |
B1-2 | 03h33m18s | +31d18m31s | 2.89 | 0.13 | 3.34 | 3.08 | 4.32 | 3.99 |
B1-3 | 03h33m18s | +31d04m44s | 74.93 | 0.31 | 0.58 | 0.45 | 0.66 | 0.51 |
B1-4 | 03h32m43s | +30d58m41s | 12.32 | 0.21 | 2.03 | 1.50 | 2.12 | 1.57 |
B1-5 | 03h32m24s | +30d48m01s | 5.48 | 0.11 | 1.62 | 1.44 | 2.25 | 2.00 |
B1-6 | 03h31m26s | +30d44m00s | 3.99 | 0.10 | 3.47 | 2.52 | 3.89 | 2.83 |
L1448-1 | 03h25m41s | +30d42m49s | 58.84 | 0.28 | 1.15 | 0.54 | 1.30 | 0.61 |
L1451-1 | 03h25m35s | +30d20m17s | 2.56 | 0.11 | 2.23 | 1.99 | 2.55 | 2.28 |
L1451-2* | 03h24m33s | +30d22m09s | 0.16 | 0.03 | 11.35 | 11.32 | 10.74 | 10.71 |
L1455-1 | 03h27m50s | +30d10m00s | 18.41 | 0.25 | 2.25 | 1.42 | 2.56 | 1.63 |
L1455-2 | 03h28m07s | +30d04m55s | 3.97 | 0.14 | 2.40 | 1.99 | 2.44 | 2.02 |
L1455-3 | 03h27m36s | +29d57m10s | 0.55 | 0.06 | 5.64 | 5.34 | 6.17 | 5.85 |
NGC 1333-1 | 03h30m00s | +31d37m39s | 2.76 | 0.12 | 1.99 | 1.79 | 2.52 | 2.26 |
NGC 1333-2 | 03h29m29s | +31d34m32s | 0.82 | 0.07 | 4.12 | 3.97 | 4.78 | 4.60 |
NGC 1333-3 | 03h29m29s | +31d31m51s | 1.77 | 0.09 | 2.31 | 2.20 | 3.11 | 2.97 |
NGC 1333-4 | 03h29m28s | +31d26m37s | 7.44 | 0.20 | 2.27 | 1.82 | 2.37 | 1.91 |
NGC 1333-5 | 03h29m05s | +31d18m36s | 20.97 | 0.19 | 1.76 | 1.35 | 2.11 | 1.61 |
NGC 1333-6 | 03h29m10s | +31d12m55s | 80.46 | 0.29 | 1.29 | 0.73 | 1.63 | 0.93 |
NGC 1333-7 | 03h28m42s | +31d13m29s | 15.24 | 0.24 | 3.83 | 1.56 | 4.58 | 1.87 |
NGC 1333-8 | 03h28m44s | +31d04m06s | 8.00 | 0.17 | 3.08 | 1.64 | 3.50 | 1.86 |
Perseus-1 | 03h30m40s | +30d25m04s | 0.95 | 0.06 | 3.62 | 3.46 | 4.38 | 4.18 |
Perseus-2 | 03h30m22s | +30d21m52s | 1.28 | 0.10 | 3.60 | 3.37 | 4.36 | 4.08 |
IC 348-1 | 03h45m23s | +32d03m24s | 0.92 | 0.09 | 4.68 | 4.40 | 5.67 | 5.33 |
IC 348-2 | 03h45m07s | +31d59m09s | 1.02 | 0.09 | 5.08 | 4.95 | 5.37 | 5.23 |
IC 348-3 | 03h44m57s | +31d59m03s | 1.14 | 0.08 | 5.26 | 5.01 | 5.84 | 5.57 |
IC 348-4 | 03h44m43s | +31d56m54s | 0.95 | 0.07 | 8.74 | 6.63 | 9.67 | 7.34 |
IC 348-5 | 03h44m26s | +31d57m31s | 5.68 | 0.17 | 2.96 | 2.33 | 3.25 | 2.56 |
IC 348-6 | 03h44m03s | +32d00m42s | 36.46 | 0.28 | 1.26 | 0.87 | 1.42 | 0.98 |
B18-1 | 04h35m46s | +24d07m49s | 3.62 | 0.08 | 1.81 | 1.29 | 2.17 | 1.55 |
B18-2 | 04h32m53s | +24d22m50s | 4.53 | 0.09 | 1.32 | 1.13 | 1.45 | 1.25 |
B18-3 | 04h32m04s | +24d30m20s | 2.33 | 0.08 | 2.01 | 1.52 | 2.57 | 1.94 |
B18-4 | 04h30m12s | +24d24m10s | 0.36 | 0.04 | 4.82 | 4.69 | 5.12 | 4.98 |
B18-5 | 04h29m30s | +24d33m22s | 2.68 | 0.06 | 1.28 | 1.21 | 1.41 | 1.33 |
B18-6 | 04h27m06s | +24d39m52s | 0.31 | 0.05 | 7.73 | 7.62 | 7.76 | 7.64 |
HC 2-1 | 04h41m39s | +25d59m46s | 4.77 | 0.10 | 1.00 | 0.89 | 1.05 | 0.93 |
HC 2-2 | 04h41m35s | +25d44m15s | 10.01 | 0.16 | 2.04 | 1.28 | 1.67 | 1.05 |
HC 2-3 | 04h40m38s | +25d28m03s | 0.08 | 0.04 | 36.92 | 35.29 | 36.74 | 35.12 |
HC 2-4 | 04h39m36s | +26d24m56s | 1.70 | 0.07 | 1.83 | 1.75 | 2.01 | 1.92 |
HC 2-5 | 04h39m24s | +25d50m26s | 0.51 | 0.05 | 8.53 | 7.61 | 9.76 | 8.71 |
HC 2-6* | 04h39m45s | +25d39m39s | 1.18 | 0.08 | 3.81 | 3.15 | 4.03 | 3.34 |
IC 5146-1 | 21h47m31s | +47d31m17s | 65.38 | 0.43 | 1.51 | 1.11 | 2.00 | 1.47 |
IC 5146-2 | 21h47m15s | +47d31m37s | 7.19 | 0.17 | 3.87 | 3.61 | 4.28 | 4.00 |
IC 5146-3 | 21h46m08s | +47d34m15s | 17.70 | 0.27 | 1.71 | 1.60 | 2.06 | 1.93 |
IC 5146-4 | 21h45m14s | +47d32m02s | 23.69 | 0.34 | 1.37 | 1.20 | 1.48 | 1.30 |
IC 5146-5 | 21h45m06s | +47d38m21s | 38.07 | 0.36 | 1.39 | 1.12 | 1.76 | 1.42 |
Cepheus L1228-1 | 20h58m07s | +77d42m03s | 3.07 | 0.10 | 1.76 | 1.66 | 2.00 | 1.88 |
Cepheus L1228-2 | 20h57m09s | +77d39m51s | 0.57 | 0.06 | 5.39 | 5.36 | 5.90 | 5.86 |
Cepheus L1228-3 | 20h57m44s | +77d34m26s | 11.11 | 0.16 | 1.88 | 1.20 | 2.34 | 1.50 |
Cepheus L1251-1* | 22h39m09s | +75d09m56s | 56.32 | 0.32 | 2.05 | 0.65 | 2.30 | 0.73 |
Cepheus L1251-2* | 22h36m15s | +75d17m16s | 4.25 | 0.14 | 1.77 | 1.48 | 2.26 | 1.88 |
Cepheus L1251-3 | 22h30m55s | +75d12m09s | 26.25 | 0.26 | 1.01 | 0.53 | 1.00 | 0.53 |
Cepheus L1251-4 | 22h28m25s | +75d12m29s | 4.66 | 0.16 | 1.67 | 1.49 | 1.91 | 1.70 |
CrA west-1 | 19h01m55s | −36d58m00s | 7.64 | 0.09 | 3.01 | 2.22 | 3.30 | 2.43 |
CrA east-1 | 19h10m26s | −37d09m45s | 1.84 | 0.05 | 1.81 | 1.69 | 2.11 | 1.97 |
L1688-1 | 16h29m05s | −24d22m08s | 0.85 | 0.05 | 3.11 | 3.00 | 3.61 | 3.48 |
L1688-2 | 16h28m38s | −24d19m58s | 0.94 | 0.06 | 3.23 | 2.98 | 3.70 | 3.42 |
L1688-3 | 16h28m29s | −24d37m48s | 0.45 | 0.05 | 5.21 | 5.03 | 6.71 | 6.48 |
L1688-4 | 16h28m06s | −24d35m03s | 1.04 | 0.06 | 3.08 | 2.91 | 4.02 | 3.79 |
L1688-5 | 16h27m42s | −24d43m57s | 2.48 | 0.09 | 8.29 | 5.39 | 8.11 | 5.28 |
L1688-6 | 16h27m26s | −24d29m39s | 13.28 | 0.12 | 2.47 | 1.89 | 2.64 | 2.02 |
L1688-7 | 16h27m02s | −24d34m14s | 6.89 | 0.12 | 1.74 | 1.56 | 2.07 | 1.85 |
L1688-8 | 16h26m35s | −24d25m00s | 5.83 | 0.07 | 1.75 | 1.54 | 2.54 | 2.25 |
L1689-1 | 16h32m34s | −24d30m10s | 11.67 | 0.09 | 0.83 | 0.53 | 1.12 | 0.72 |
L1689-2 | 16h31m48s | −24d51m46s | 5.23 | 0.08 | 1.45 | 1.24 | 1.66 | 1.43 |
L1689-3* | 16h32m03s | −24d58m59s | 4.68 | 0.11 | 3.12 | 1.99 | 3.38 | 2.15 |
Serpens MWC 297-1 | 18h28m15s | −03d48m49s | 16.44 | 0.19 | 1.03 | 0.84 | 1.24 | 1.00 |
Orion A-1 | 05h35m11s | −04d57m12s | 8.95 | 0.16 | 2.19 | 2.00 | 2.42 | 2.22 |
Orion A-2 | 05h35m11s | −05d37m50s | 46.92 | 0.28 | 4.15 | 1.56 | 4.59 | 1.72 |
Orion A-3 | 05h35m14s | −05d53m21s | 7.09 | 0.16 | 5.22 | 2.63 | 6.76 | 3.41 |
Orion A-4 | 05h36m17s | −06d12m15s | 8.60 | 0.14 | 2.82 | 2.31 | 3.97 | 3.25 |
Orion A-5 | 05h35m17s | −06d15m18s | 13.53 | 0.15 | 6.31 | 3.86 | 7.54 | 4.62 |
Orion A S-1 | 05h39m12s | −07d13m22s | 10.76 | 0.22 | 4.83 | 1.60 | 5.52 | 1.83 |
Orion A S-2 | 05h39m37s | −07d24m32s | 5.58 | 0.23 | 4.26 | 1.95 | 6.46 | 2.97 |
Orion A S-3 | 05h39m19s | −07d24m15s | 36.26 | 0.33 | 1.09 | 0.84 | 1.28 | 0.98 |
Orion A S-4 | 05h40m02s | −07d28m15s | 19.81 | 0.19 | 2.03 | 0.90 | 2.33 | 1.04 |
Orion A S-5 | 05h40m27s | −07d37m43s | 20.12 | 0.23 | 1.32 | 0.76 | 1.43 | 0.82 |
Orion A S-6 | 05h40m30s | −07d44m18s | 14.89 | 0.17 | 1.02 | 0.89 | 1.26 | 1.10 |
Orion B NGC 2023-2024-1 | 05h41m27s | −01d43m38s | 0.72 | 0.07 | 14.87 | 14.05 | 15.51 | 14.65 |
Orion B NGC 2023-2024-2 | 05h41m17s | −01d47m55s | 0.76 | 0.08 | 11.70 | 11.07 | 12.74 | 12.06 |
Orion B NGC 2023-2024-3 | 05h41m49s | −01d55m38s | 30.73 | 0.10 | 1.22 | 1.16 | 1.65 | 1.56 |
Orion B NGC 2023-2024-4 | 05h41m51s | −01d58m00s | 73.42 | 0.11 | 0.79 | 0.76 | 1.45 | 1.40 |
Orion B NGC 2023-2024-5 | 05h41m56s | −02d01m10s | 1.22 | 0.09 | 8.51 | 7.57 | 9.49 | 8.44 |
Orion B NGC 2023-2024-6 | 05h41m38s | −02d19m08s | 50.29 | 0.32 | 1.97 | 1.30 | 2.00 | 1.32 |
Orion B NGC 2023-2024-7 | 05h41m37s | −02d25m49s | 4.42 | 0.15 | 7.77 | 6.82 | 7.89 | 6.94 |
Orion B NGC 2068-2071-1 | 05h46m10s | −00d16m31s | 8.35 | 0.18 | 2.86 | 2.67 | 3.73 | 3.49 |
Note. R.A. and decl. values quoted refer to the geometric center of each clump. Asterisks denote clumps whose maps include regions missing NH3 data, most of which have very small covering fraction.
In Figure 3, we adopt our fiducial method (Abel reconstruction) for removing foreground and background matter, and compare the outcomes of other choices in the estimation of α. Two trends are immediately clear. First, αBM92 and αSMJ19 are very similar across our entire sample. Because we use the same data to derive both quantities, the comparison provides a statistical calibration of BM92's self-gravity parameter a, which our analysis shows is only slightly above unity (geometric mean value 1.125) in our Gould Belt clumps. For more detail, examining the ratio of the fits plotted in Figure 3, we see that a grows from unity to a maximum of 1.15 across the range of Mcl.
Second, excluding suppresses α by an amount that changes with Mcl. The difference is insignificant for Mcl < M⊙, where clumps harbor subsonic motions and are not highly resolved. However, the influence of
grows with mass and amounts to 0.26 dex (a factor of 1.81) for Mcl ∼ 100 M⊙. This is important, as most analyses adopt for α a typical value of the total beam-wise line width σeff(x, y), which has the effect of neglecting
in the estimate of
. (There do exist counterexamples, however: for instance, Roman-Duval et al. 2010's Equation (7) includes all of
.)
Methods for extracting Σcl from the column density map Σ also influence the virial parameter. In Figure 4, we compare several methods: removing emission from around the clump using our default method of Abel reconstruction, removing an interpolated value from the clump boundary, or making no correction (in which case Σcl = Σ within the boundary). The trends in α(Mcl) are shown for these different choices.
Figure 4. A demonstration that the method used to extract the clump column Σcl from the column density map Σ affects the derived virial ratio-mass relation. Solid red diamonds represent our fiducial Abel reconstruction, as in Figure 3; open cyan diamonds represent the removal of a bi-cubic interpolation from the cloud edge; and open orange symbols indicate no correction (all dust emission projected within the clump boundary is attributed to Σcl). Pink solid and purple dashed lines connect values for the same clump. The analysis otherwise follows SMJ19. Overplotted curves are quadratic fits in log space to the point families of the same color.
Download figure:
Standard image High-resolution imageSimple interpolation removes the most mass, while Abel reconstruction, by accounting for projection effects, attributes less material to the envelope. These differences are greatest for the least massive clumps, whose column density contrasts the least with their surroundings. The choice of extraction method affects the inferred virial parameter, mostly through the relation that applies when Σcl is multiplied by a constant. However, the effect on α diminishes for higher clump masses. The effect is also roughly in the direction of the trend line α(Mcl), that is, it is partly degenerate with the original trend. These facts diminish the influence of the extraction method on the trend of α with clump mass.
In summary, although α drops with Mcl within our collection of clumps, it is consistently below unity at high masses only when one ignores spatially resolved component of the bulk kinetic energy, when one overestimates the gravitational self-energy by adopting a > 1.1, or when one does not correct at all for foreground and background matter projected within the clump boundary. Using our most complete analysis, which addresses these shortcomings, we find no evidence that the massive clumps in our sample are consistently subvirial.
Note that the clumps in our sample typically have αSMJ19 ∼5 for Mcl = 1 M⊙, while studies of dense cores, such as that of Johnstone et al. (2000), find that cores approach the critical state of a Bonnor–Ebert sphere (in which α = 1.683) at around a solar mass. The discrepancy is not surprising, considering that our algorithm defines clumps that include as much NH3 emission as possible. Any compact cores in our sample are therefore enclosed within larger clumps.
5.1. Sources of Error
We pause here to discuss several sources of error, approximately in descending order of their expected importance for our conclusions.
- 1.Despite being well calibrated, our estimated column densities are subject to systematic errors arising from the assumed dust emissivity and dust-to-gas ratio. Our assumptions that β is constant along each line of sight, and of a universal ratio between dust optical depth and column density, both introduce systematic uncertainties into our determination of Σcl. Variations of the dust opacity have primarily been observed in densest regions of compact cores and filaments (e.g., Chacón-Tanarro et al. 2019; Howard et al. 2019), implying that the impact of line-of-sight variations on our clump column densities is minor but not negligible.We note that the range of inferred β values within each of our clumps is extremely limited, in that the median clump only spans a range Δβ = 0.005 and the largest variation within a clump is Δβ = 0.051. A caveat, considering the limited resolution of Planck, is that these values could be significantly underestimated. Overall, however, we consider submillimeter dust emission a more reliable tracer than NH3 for determining column density, so long as the clump profile can be separated from other emission along the line of sight.
- 2.Our analysis, which is optimized for well-resolved clumps, does not allow an explicit correction for finite resolution of the type introduced by Rosolowsky & Leroy (2006). In this regard, it is useful to note that the clump angular diameter 2Rcl/distance correlates with clump mass, ranging from 44''–8' across our sample, with median values of 2.8' and 3.8' for those clumps with Mcl < 10 M⊙ and Mcl > 10 M⊙, respectively. Comparing to the GAS resolution of 32'', we infer that finite resolution is likely to affect clump selection and to bias our estimates of α, but predominantly for the low-mass clumps that are not our primary focus.
- 3.We work only with projected data and line-of-sight velocities. As was discussed in SMJ19, this leads to errors in both the numerator and denominator of αSMJ19 that are random insofar as the viewing angle is random (but become systematic if clumps are aligned in a larger structure, as Alves et al. 2020 find). These errors are greater for anisotropic structures and velocity fields, and so they tend to be greater for more massive clumps. Projection-dependent errors are surely responsible for some of the scatter we see in Figure 3, which amounts to 0.21 dex standard deviation around the fitted curves. Notably, this is comparable to variation with viewing angle in the simulation examined by SMJ19. One might be able to calibrate this component of the scatter using a statistical sample of simulated clumps, although that would be beyond the scope of the current work.
- 4.We estimate
using parameters obtained from a single-component fit to the NH3 emission, which could miss some kinetic energy from other components that represent fluid velocities within the clump. Although some differences are indeed seen when single-component and multiple-component fits are compared (Choudhury et al. 2020, 2021), these amount to small corrections for the clump kinematics. Moreover, at GAS sensitivity, most pixels are well fit with single components (e.g., Sokolov et al. 2020).
- 5.To remove foreground and background dust emission from Σcl requires deprojection, which necessarily implies some error (Beaumont et al. 2013). SMJ19 found that our fiducial method based on the Abel transform performs better than either using simple interpolation or making no correction. However, it probably implies both random and systematic errors. Because the choice of method has a greater impact at lower clump masses, we infer that these errors declines with Mcl and are quite minor for the massive clumps of greatest interest here.
- 6.Although we require that our clumps exceed the critical densities of the NH3 (1,1) and (2,2) transitions, excitation variations like those seen by Crapsi et al. (2007) are likely to add biases to the data we use to estimate
and hence α, as argued by Traficante et al. (2018b). We also note that NH3 data could be contaminated by emission projected within the clump boundary, which would cause α to be slightly overestimated (e.g., Choudhury et al. 2020).
- 7.We omit any component not traced by dust and NH3 emission, such as embedded stars and high-velocity outflows. In the case of protostellar outflows, this means that our definition of
applies only to matter at velocities within a few σcl,z of the systemic velocity, so that outflows must be treated as surface terms within the virial theorem (as we discussed in Section 1). An alternative definition of
would explicitly include the outflow kinetic energy. In the case of protostars, our definition means that
lacks the contribution from stellar gravity. As the two effects both correlate with star formation, we expect them to be most important for the most massive, lowest-α clumps in our sample. If included in α, we expect they would affect it in opposite ways, with the positive effect of outflows being more significant than the negative effect of stellar gravity. In the active region NGC 1333, for example, ∼ 20% of the total mass is in protostars (Gutermuth et al. 2008, see Matzner & Jumper 2015), a value typical of embedded protoclusters (Lada & Lada 2003), while the protostellar outflow energy is comparable to
(Plunkett et al. 2013).
- 8.Our method for evaluating gravitational energy converges when Σcl is well resolved (as discussed in the Appendix), but underestimates the magnitude of
when there exists unresolved structure. However, the fact that a is very close to its value for a uniform sphere is strong indication that the gravitational energy in unresolved structure is very minor.
Of these, the systematic errors will affect our conclusions regarding α(Mcl), while random errors will tend to average out.
It is important to note, however, that the conclusions we draw when comparing analysis choices are independent of measurement error, because we use identical data to make these comparisons.
6. Discussion
As first stressed by Kauffmann et al. (2013), many studies have found massive clumps to be moderately to strongly subvirial (Roman-Duval et al. 2010; Pillai et al. 2011; Wienen et al. 2012; Ragan et al. 2012; Li et al. 2013; Tan et al. 2013; Urquhart et al. 2014; Friesen et al. 2016; Svoboda et al. 2016; Kirk et al. 2017b; Keown et al. 2017; Redaelli et al. 2017; Csengeri et al. 2017; Traficante et al. 2018a; Keown et al. 2019; Chen et al. 2019; Kerr et al. 2019; Billington et al. 2020; Traficante et al. 2020). In our fiducial analysis, which follows SMJ19, we do not. We trace the difference to several key features.
First, we account for the component of kinetic energy associated with spatial variations of the line-of-sight velocity in the plane of the sky, which we call . BM92 include this component (see their Appendix C), but many later studies exclude it by adopting a typical value of σeff(x, y) for σcl in the evaluation of αBM92. There do exist exceptions, such as Roman-Duval et al. (2010), who include bulk energy in their Equation (7), and works like Colombo et al. (2019) in which σcl,z
is defined in terms of the intensity-weighted variance of velocity (as adopted, for instance, by Solomon et al. 1987). Including
matters most for massive clumps: these contain supersonic internal motions, but their dense substructures tend to have narrow internal line widths with subsonic nonthermal components. When the clump is well resolved, these substructures will tend to set the line width within each beam. The effect therefore varies with resolution. This has been observed, for instance, in B5 by Pineda et al. (2010, 2015) and in Orion by Friesen et al. (2017) and Monsch et al. (2018).
Second, we compute the gravitational self-energies of our clumps using the robust method of SMJ19. Applied to an ensemble of clumps, this provides a sensitive calibration of the BM92 parameter a. An average over our clump sample indicates a ≃ 1.13, while fits to the mass dependence indicate that a grows from unity for low clump masses to ∼1.15 for massive clumps. Because our clumps have generally modest projected aspect ratios (with the exception of HC 2-2), these represent low values of central concentration: the low-mass result a ≃ 1 is consistent with these cores being essentially uniform-density, pressure-confined objects, and our entire sample is less centrally concentrated than a critical Bonnor–Ebert sphere (for which a = 1.22). Other analyses often ignore a, effectively taking a = 1 (e.g., Billington et al. 2020); or they adopt a specific value, such as a = 1.25 (corresponding to k = 3/2: e.g., Keown et al. 2019) or a = 1.67 (corresponding to k = 2: e.g., Roman-Duval et al. 2010; Csengeri et al. 2017). We note that using the definition of Mvir in Rohlfs & Wilson (2004) to estimate α = Mvir/Mcl (e.g., Wienen et al. 2012) amounts to adopting a = 5/6. We also note that our sample does not support the suggestion by Ballesteros-Paredes et al. (2018) that cases in which clumps are inferred to be unbound (α > 2) might result from their values of a being significantly underestimated.
Third, while we use dust emission to trace mass, we are careful to account for foreground and background emission because overestimating Mcl leads one to underestimate α, roughly according to the scaling . For this we prefer the Abel reconstruction method advanced in SMJ19.
On this point, we note that prior studies determine clump masses and profiles using a wide range of data sources and analysis techniques (see Kauffmann et al. 2013), not all of which exclude foreground and background contamination. Common approaches include subtracting an interpolated background, matching and extracting a template profile, or making no correction. Dendrogram analysis (Rosolowsky et al. 2008) can employ any of these, although by default it makes no correction. Approaches tailored to extract compact cores tend to fit templates (e.g., Johnstone et al. 2000) or subtract an estimated background, as in Men'shchikov et al. (2012). As we mentioned in Section 2, the spatial filtering afforded by interferometric or chopped observations accomplishes background subtraction in an approximate way.
In Figure 5, we review the effect of each of the three key features of our analysis mentioned above. We note that while foreground/background subtraction can dominate the error for an individual clump, this effect is partly degenerate with the trend in α(Mcl), as shown in Figure 4. Focusing on the mass scale 100 M⊙, we see that neglecting depresses one's estimate for αSMJ19 by about a factor of 1.8 (depending slightly on the clump extraction method). Including foreground and background emission suppresses it by a further factor of 1.4. Adopting the extreme choice a = 1.67 over the more accurate value of 1.15 reduces α further by a factor of 1.45. Compounding these leads to an underestimate of α by a factor of 3.6.
Figure 5. Estimates for α compared. We present only quadratic fits, suppressing individual clump data for clarity; each fit is plotted only over the appropriate range of Mcl. In each line pair, is included on the top line and excluded on the bottom line. Analysis choices can lead α being underestimated by a factor of up to 3.5 at Mcl = 100 M⊙, relative to αSMJ19 with Abel reconstruction and including
(our recommended best-practices evaluation).
Download figure:
Standard image High-resolution imageThis analysis indicates that, while all three factors have comparable impact, bulk kinetic energy is the most important effect at the reference mass of 100 M⊙, and the only one that is growing more important with increasing clump mass. Although its impact is clearly resolution-dependent, this behavior indicates that bulk kinetic energy might dominate at masses higher than those included in our sample.
A further and potentially significant difference is that we associate each clump with a definite boundary on the sky, which we derive from the emission of a reliable dense gas tracer (NH3) and which we use to infer quantities like the clump radius and gravitational energy. Compared to approaches in which radii are defined by moments of the emission, or by fitting to predetermined model templates, ours has the advantage that the derived quantities reflect those of the three-dimensional clump in question in a mathematically rigorous way (Equations (5) and (7)).
In summary, we conclude that it is easy to underestimate α by a significant margin as a result of choices in the analysis. Addressing these points carefully, as we do here, alleviates the physical puzzle that a finding of consistently subvirial clumps would pose—namely, that while a mildly subvirial population may be close to equilibrium (McKee & Zweibel 1995; Tan et al. 2013) for the observed degree of magnetization (Crutcher 2012), this is not true of the strongly subvirial state. And, although one could prepare initial conditions far from equilibrium in which α ≪ 1, this state would be ephemeral, persisting only for a fraction of a freefall time, as noted, for instance, by Kauffmann et al. (2013). Moreover, it is not clear how a dense clump could be assembled in this state without concomitantly developing velocities at the virial scale.
Within our Gould Belt sample, we find that clump motions are, in fact, virial in magnitude, with α ≃ 1.34 and little dependence on mass at the upper end of our sample. This impacts any conclusions that depend on the virial ratio, such as the expected mode of protostellar accretion. It remains to be seen whether conclusions drawn from this sample will be relevant to samples that include more massive clumps.
Our result that the virial parameter levels off at a value of order unity, rather than continuing to decline with increasing clump mass, points to a change in physical or evolutionary state for the most massive clumps in our sample. Could this reflect a greater role for dynamical feedback from star formation? This would be consistent with the fact that protostellar outflows are active in a number of our massive clumps (Kun et al. 2008; Devine et al. 2009; Arce et al. 2010; Curtis et al. 2010; Feddersen et al. 2020). It is also consistent with the fact that class 0 protostars are highly segregated toward these massive clumps and the structures that host them, as is evident in Orion A (e.g., Stutz & Kainulainen 2015) and Perseus (e.g., Mercimek et al. 2017).
It is already well established that the local star formation rate has a threshold-like (Onishi et al. 1998) or power-law (Heiderman et al. 2010) dependence on Σ, and that the recent history of star formation correlates with the slope of the local probability density function of Σ (Sadavoy 2013; Stutz & Kainulainen 2015). We find that the clump-averaged values of Σ and Σcl both correlate with Mcl and anticorrelate with α, so the former rule, at least, is qualitatively consistent with and anticorrelation between α and star formation intensity. We consider this a promising avenue for future work.
We thank our referee and Christopher McKee for useful feedback and important suggestions. A.S. thanks Tina Peters for advice on data analysis techniques. The research of A.S. and C.D.M. is supported by an NSERC Discovery Grant. A.P. acknowledges the support from the Russian Ministry of Science and Higher Education via the State Assignment Project FEUZ-2020-0038. A.P. is a member of the Max Planck Partner Group at the Ural Federal University. We acknowledge that this work was conducted on the traditional land of the Huron-Wendat, the Seneca, and the Mississaugas of the Credit; we are grateful for the opportunity to work here.
Software: Astropy (Astropy Collaboration et al. 2013, 2018), Scikit-learn (Pedregosa et al. 2011).
Appendix
We describe here our procedure for managing discretization errors in the computation of .
SMJ19 suggest computing from Σcl(x, y) via
![Equation (A1)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn9.gif)
where
![Equation (A2)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn10.gif)
This procedure is perfectly valid in the continuous limit, but requires a small modification to account for the fact that Σcl is sampled on a discrete grid at some resolution δ
x. The most important consequence of discreteness is that, because the Green's function diverges at the origin, its central element must be replaced with some finite value − ε
G/δ
x, where ε parameterizes the self-energy of a single pixel. We adopt
![Equation (A3)](https://content.cld.iop.org/journals/0004-637X/922/1/87/revision1/apjac20d2eqn11.gif)
the value that correctly yields the self-energy of a square of width δ
x in the case that Σcl(x, y) is uniform. We test this choice using discretized projections of structures of known self-gravitational energy, such as uniform-density spheres, and find that it reduces the error in by a factor of about 2.4 relative to the outcome of choosing ε = 0. Significantly smaller improvements can be obtained by adjusting the noncentral elements of the Green's function to account for the finite sizes of pixels. However, we do not implement these because the error can further be reduced by interpolating Σcl to a finer grid before evaluating
. Finding that the relative error when using ε is ∼ δ
x/2Rcl, we adopt a strategy in which each clump's column density is resampled so that Rcl > 100δ
x, implying errors well below 1%.
Footnotes
- 17
Ephemeral nonequilibrium states are also possible. For instance, α approaches 2 from below in asymptotic, nonrotating, unmagnetized freefall (Ballesteros-Paredes et al. 2018), whereas α ≫ 2 for explosive motions in excess of the escape velocity.
- 18
We exclude a clump in Barnard 59 from our sample, on the grounds that it lacks NH3 data over an extensive region including the column density peak.