Abstract
Observational studies of star formation reveal spatial distributions of young stellar objects (YSOs) that are ‘snapshots’ of an ongoing star formation process. Using methods from spatial ...statistics it is possible to test the likelihood that a given distribution process could produce the observed patterns of YSOs. The aim of this paper is to determine the usefulness of the spatial statistics tests Diggle’s G function (G), the ‘free-space’ function (F), Ripley’s K, and O-ring for application to astrophysical data. The spatial statistics tests were applied to simulated data containing 2D Gaussian clusters projected on random distributions of stars. The number of stars within the Gaussian cluster and number of background stars were varied to determine the tests’ ability to reject complete spatial randomness (CSR) with changing signal-to-noise. The best performing test was O-ring optimized with overlapping logarithmic bins, closely followed by Ripley’s K. The O-ring test is equivalent to the two-point correlation function. Both F and G (and the minimum spanning tree, of which G is a subset) performed significantly less well, requiring a cluster with a factor of two higher signal-to-noise in order to reject CSR consistently. We demonstrate the tests on example astrophysical datasets drawn from the Spitzer Gould's Belt Survey catalogue.
ABSTRACT
Radiative feedback can influence subsequent star formation. We quantify the heating from OB stars in the local star-forming regions in the James Clerk Maxwell Telescope (JCMT) Gould Belt ...Survey. Dust temperatures are calculated from 450/850 $\mu$m flux ratios from SCUBA-2 observations at the JCMT assuming a fixed dust opacity spectral index β = 1.8. Mean dust temperatures are calculated for each submillimetre clump along with projected distances from the main OB star in the region. Temperature versus distance is fitted with a simple model of dust heating by the OB star radiation plus the interstellar radiation field and dust cooling through optically thin radiation. Classifying the heating sources by spectral type, O-type stars produce the greatest clump average temperature rises and largest heating extent, with temperatures of over 40 K and significant heating out to at least 2.4 pc. Early-type B stars (B4 and above) produce temperatures of over 20 K and significant heating over 0.4 pc. Late-type B stars show a marginal heating effect within 0.2 pc. For a given projected distance, there is a significant scatter in clump temperatures that is due to local heating by other luminous stars in the region, projection effects, or shadowing effects. Even in these local, ‘low-mass’ star-forming regions, radiative feedback is having an effect on parsec scales, with 24 per cent of the clumps heated to at least 3 K above the 15 K base temperature expected from heating by only the interstellar radiation field, and a mean dust temperature for heated clumps of 24 K.
ABSTRACT
Star formation is known to occur more readily where more raw materials are available. This is often expressed by a ‘Kennicutt–Schmidt’ relation where the surface density of young stellar ...objects (YSOs) is proportional to column density to some power, μ. The aim of this work was to determine if column density alone is sufficient to explain the locations of Class 0/I YSOs within Serpens South, Serpens Core, Ophiuchus, NGC 1333, and IC348, or if there is clumping or avoidance that would point to additional influences on the star formation. Using the O-ring test as a summary statistic, 95 per cent confidence envelopes were produced for different values of μ from probability models made using the Herschel column density maps. The YSOs were tested against four distribution models: the best estimate of μ for the region, μ = 0 above a minimum column density threshold, and zero probability elsewhere, μ = 1, and the power law that best represents the five regions as a collective, μ = 2.05 ± 0.20. Results showed that μ = 2.05 model was consistent with the majority of regions and, for those regions, the spatial distribution of YSOs at a given column density is consistent with being random. Serpens South and NGC 1333 rejected the μ = 2.05 model on small scales of ${\sim}0.15\, \mathrm{pc}$ which implies that small-scale interactions may be necessary to improve the model. On scales above 0.15 pc, the positions of YSOs in all five regions can be well described using column density alone.
Performing ground-based submillimetre observations is a difficult task as the measurements are subject to absorption and emission from water vapour in the Earth's atmosphere and time variation in ...weather and instrument stability. Removing these features and other artefacts from the data is a vital process which affects the characteristics of the recovered astronomical structure we seek to study. In this paper, we explore two data reduction methods for data taken with the Submillimetre Common-User Bolometer Array-2 (SCUBA-2) at the James Clerk Maxwell Telescope (JCMT). The JCMT Legacy Reduction 1 (JCMT LR1) and The Gould Belt Legacy Survey Legacy Release 1 (GBS LR1) reduction both use the same software (starlink) but differ in their choice of data reduction parameters. We find that the JCMT LR1 reduction is suitable for determining whether or not compact emission is present in a given region and the GBS LR1 reduction is tuned in a robust way to uncover more extended emission, which better serves more in-depth physical analyses of star-forming regions. Using the GBS LR1 method, we find that compact sources are recovered well, even at a peak brightness of only three times the noise, whereas the reconstruction of larger objects requires much care when drawing boundaries around the expected astronomical signal in the data reduction process. Incorrect boundaries can lead to false structure identification or it can cause structure to be missed. In the JCMT LR1 reduction, the extent of the true structure of objects larger than a point source is never fully recovered.
We develop and apply a Hessian-based filament detection algorithm to submillimetre continuum observations of Orion A North. The resultant filament radial density profiles are fitted with ...beam-convolved line-of-sight Plummer-profiles using Markov chain Monte Carlo techniques. The posterior distribution of the radial decay parameter demonstrates that the majority of filaments exhibit p = 1.5–3, with a mode at p = 2.2, suggesting deviation from the Ostriker p = 4 isothermal, equilibrium, self-gravitating cylinder. The spatial distribution of young stellar objects relative to the high column density filaments is investigated, yielding a lower limit on the star-forming age of the integral-shaped filament ∼1.4 Myr. Additionally, inferred lifetimes of filaments are examined which suggest long-term filament accretion, varying rates of star formation, or both. Theoretical filament stability measures are determined with the aid of HARP C18O J = 3–2 observations and indicate that the majority of filaments are gravitationally subcritical, despite the presence of young protostars. The results from this investigation are consistent with the one-dimensional accretion flow filament model recently observed in numerical simulations.
ABSTRACT The Orion A molecular cloud is one of the most well-studied nearby star-forming regions, and includes regions of both highly clustered and more dispersed star formation across its full ...extent. Here, we analyze dense, star-forming cores identified in the 850 and 450 m SCUBA-2 maps from the JCMT Gould Belt Legacy Survey. We identify dense cores in a uniform manner across the Orion A cloud and analyze their clustering properties. Using two independent lines of analysis, we find evidence that clusters of dense cores tend to be mass segregated, suggesting that stellar clusters may have some amount of primordial mass segregation already imprinted in them at an early stage. We also demonstrate that the dense core clusters have a tendency to be elongated, perhaps indicating a formation mechanism linked to the filamentary structure within molecular clouds.
In this paper, we present the first observations of the Ophiuchus molecular cloud performed as part of the James Clerk Maxwell Telescope (JCMT) Gould Belt Survey (GBS) with the SCUBA-2 instrument. We ...demonstrate methods for combining these data with previous HARP CO, Herschel, and IRAM N2H+ observations in order to accurately quantify the properties of the SCUBA-2 sources in Ophiuchus. We produce a catalogue of all of the sources found by SCUBA-2. We separate these into protostars and starless cores. We list all of the starless cores and perform a full virial analysis, including external pressure. This is the first time that external pressure has been included in this level of detail. We find that the majority of our cores are either bound or virialized. Gravitational energy and external pressure are on average of a similar order of magnitude, but with some variation from region to region. We find that cores in the Oph A region are gravitationally bound prestellar cores, while cores in the Oph C and E regions are pressure-confined. We determine that N2H+ is a good tracer of the bound material of prestellar cores, although we find some evidence for N2H+ freeze-out at the very highest core densities. We find that non-thermal linewidths decrease substantially between the gas traced by C18O and that traced by N2H+, indicating the dissipation of turbulence at higher densities. We find that the critical Bonnor–Ebert stability criterion is not a good indicator of the boundedness of our cores. We detect the pre-brown dwarf candidate Oph B-11 and find a flux density and mass consistent with previous work. We discuss regional variations in the nature of the cores and find further support for our previous hypothesis of a global evolutionary gradient across the cloud from south-west to north-east, indicating sequential star formation across the region.
We present a global study of low-mass, young stellar object (YSO) surface densities (Σ) in nearby (<500 pc) star-forming regions based on a comprehensive collection of Spitzer Space Telescope ...surveys. We show that the distribution of YSO surface densities in the solar neighbourhood is a smooth distribution, being adequately described by a lognormal function from a few to 103 YSOs pc−2, with a peak at ∼22 stars pc−2 and a dispersion of . We do not find evidence for multiple discrete modes of star formation (e.g. clustered and distributed). Comparing the observed surface density distribution to previously reported surface density threshold definitions of clusters, we find that the fraction of stars in clusters is crucially dependent on the adopted definitions, ranging from 40 to 90 per cent. However, we find that only a low fraction (<26 per cent) of stars are formed in dense environments where their formation/evolution (along with their circumstellar discs and/or planets) may be affected by the close proximity of their low-mass neighbours.
ABSTRACT The James Clerk Maxwell Telescope Gould Belt Legacy Survey obtained SCUBA-2 observations of dense cores within three sub-regions of Orion B: LDN 1622, NGC 2023/2024, and NGC 2068/2071, all ...of which contain clusters of cores. We present an analysis of the clustering properties of these cores, including the two-point correlation function and Cartwright's Q parameter. We identify individual clusters of dense cores across all three regions using a minimal spanning tree technique, and find that in each cluster, the most massive cores tend to be centrally located. We also apply the independent M- technique and find a strong correlation between core mass and the local surface density of cores. These two lines of evidence jointly suggest that some amount of mass segregation in clusters has happened already at the dense core stage.