We present a uniform catalog of accurate distances to local molecular clouds informed by the Gaia DR2 data release. Our methodology builds on that of Schlafly et al. First, we infer the distance and ...extinction to stars along sightlines toward the clouds using optical and near-infrared photometry. When available, we incorporate knowledge of the stellar distances obtained from Gaia DR2 parallax measurements. We model these per-star distance-extinction estimates as being caused by a dust screen with a 2D morphology derived from Planck at an unknown distance, which we then fit for using a nested sampling algorithm. We provide updated distances to the Schlafly et al. sightlines toward the Dame et al. and Magnani et al. clouds, finding good agreement with the earlier work. For a subset of 27 clouds, we construct interactive pixelated distance maps to further study detailed cloud structure, and find several clouds which display clear distance gradients and/or are comprised of multiple components. We use these maps to determine robust average distances to these clouds. The characteristic combined uncertainty on our distances is 5%-6%, though this can be higher for clouds at greater distances, due to the limitations of our single-cloud model.
ABSTRACT The all-Galaxy CO survey of Dame et al. is by far the most uniform, large-scale Galactic CO survey. Using a dendrogram-based decomposition of this survey, we present a catalog of 1064 ...massive molecular clouds throughout the Galactic plane. This catalog contains 2.5 × 108 solar masses, or 25 − 5.8 + 10.7 % of the Milky Way's estimated H2 mass. We track clouds in some spiral arms through multiple quadrants. The power index of Larson's first law, the size-linewidth relation, is consistent with 0.5 in all regions-possibly due to an observational bias-but clouds in the inner Galaxy systematically have significantly (∼30%) higher linewidths at a given size, indicating that their linewidths are set in part by the Galactic environment. The mass functions of clouds in the inner Galaxy versus the outer Galaxy are both qualitatively and quantitatively distinct. The inner Galaxy mass spectrum is best described by a truncated power law with a power index of γ = −1.6 0.1 and an upper truncation mass of M0 = (1.0 0.2) × 107 M , while the outer Galaxy mass spectrum is better described by a non-truncating power law with γ = −2.2 0.1 and an upper mass of M0 = (1.5 0.5) × 106 M , indicating that the inner Galaxy is able to form and host substantially more massive GMCs than the outer Galaxy. Additionally, we have simulated how the Milky Way would appear in CO from extragalactic perspectives, for comparison with CO maps of other galaxies.
We present a new technique to determine distances to major star-forming regions across the Perseus Molecular Cloud, using a combination of stellar photometry, astrometric data, and 12CO spectral-line ...maps. Incorporating the Gaia DR2 parallax measurements when available, we start by inferring the distance and reddening to stars from their Pan-STARRS1 and Two Micron All Sky Survey photometry, based on a technique presented by Green et al. and implemented in their 3D "Bayestar" dust map of three-quarters of the sky. We then refine their technique by using the velocity slices of a CO spectral cube as dust templates and modeling the cumulative distribution of dust along the line of sight toward these stars as a linear combination of the emission in the slices. Using a nested sampling algorithm, we fit these per-star distance-reddening measurements to find the distances to the CO velocity slices toward each star-forming region. This results in distance estimates explicitly tied to the velocity structure of the molecular gas. We determine distances to the B5, IC 348, B1, NGC 1333, L1448, and L1451 star-forming regions and find that individual clouds are located between 275 and 300 pc, with typical combined uncertainties of 5%. We find that the velocity gradient across Perseus corresponds to a distance gradient of about 25 pc, with the eastern portion of the cloud farther away than the western portion. We determine an average distance to the complex of 294 17 pc, about 60 pc further than the distance derived to the western portion of the cloud using parallax measurements of water masers associated with young stellar objects. The method we present is not limited to the Perseus Complex, but may be applied anywhere on the sky with adequate CO data in the pursuit of more accurate 3D maps of molecular clouds in the solar neighborhood and beyond.
We investigate the effect of line-of-sight temperature variations and noise on two commonly used methods to determine dust properties from dust-continuum observations of dense cores. One method ...employs a direct fit to a modified blackbody spectral energy distribution (SED); the other involves a comparison of flux ratios to an analytical prediction. Fitting fluxes near the SED peak produces inaccurate temperature and dust spectral index estimates due to the line-of-sight temperature (and density) variations. Longer wavelength fluxes in the Rayleigh-Jeans part of the spectrum ( 600 mm for typical cores) may more accurately recover the spectral index, but both methods are very sensitive to noise. The temperature estimate approaches the density-weighted temperature, or 'column temperature,' of the source as short wavelength fluxes are excluded. An inverse temperature-spectral index correlation naturally results from SED fitting, due to the inaccurate isothermal assumption, as well as noise uncertainties. We show that above some 'threshold' temperature, the temperatures estimated through the flux ratio method can be highly inaccurate. In general, observations with widely separated wavelengths, and including shorter wavelengths, result in higher threshold temperatures; such observations thus allow for more accurate temperature estimates of sources with temperatures less than the threshold temperature. When only three fluxes are available, a constrained fit, where the spectral index is fixed, produces less scatter in the temperature estimate when compared to the estimate from the flux ratio method.
ABSTRACT
We study the formation, evolution, and collapse of dense cores by tracking structures in a magnetohydrodynamic simulation of a star-forming cloud. We identify cores using the dendrogram ...algorithm and utilize machine learning techniques, including Neural Gas prototype learning and Fuzzy c-means clustering to analyse the density and velocity dispersion profiles of cores together with six bulk properties. We produce a 2-d visualization using a Uniform Manifold Approximation and Projection (UMAP), which facilitates the connection between physical properties and three partially-overlapping phases: i) unbound turbulent structures (Phase I), ii) coherent cores that have low turbulence (Phase II), and iii) bound cores, many of which become protostellar (Phase III). Within Phase II, we identify a population of long-lived coherent cores that reach a quasi-equilibrium state. Most prestellar cores form in Phase II and become protostellar after evolving into Phase III. Due to the turbulent cloud environment, the initial core properties do not uniquely predict the eventual evolution, i.e. core evolution is stochastic, and cores follow no one evolutionary path. The phase lifetimes are 1.0 ± 0.1 × 105 yr, 1.3 ± 0.2 × 105 yr, and 1.8 ± 0.3 × 105 yr for Phase I, II, and III, respectively. We compare our results to NH3 observations of dense cores. Known coherent cores predominantly map into Phase II, while most turbulent pressure-confined cores map to Phase I or III. We predict that a significant fraction of observed starless cores have unresolved coherent regions and that ≳20 per cent of observed starless cores will not form stars. Measurements of core radial profiles in addition to the usual bulk properties will enable more accurate predictions of core evolution.
We use data gathered by the COMPLETE survey of star-forming regions to find new calibrations of the 'X-factor' and super(13)CO abundance within the Perseus molecular cloud. We divide Perseus into six ...subregions, using groupings in a dust temperature vs. LSR velocity plot. The standard X- factor, image, is derived both for the whole Perseus complex and for each of the six subregions with values consistent with previous estimates. However, the X-factor is heavily affected by the saturation of the emission above image mag, and variations are also found between regions. Linear fits to relate image and image using only points below 4 mag of extinction yield a better estimate of the image than the X-factor. Linear relations of image, and image with image are derived. The extinction thresholds above which super(13)CO(1-0) and C super(18)O(1-0) are detected are about 1 mag larger than previous estimates, so that a more efficient shielding is needed for the formation of CO than previously thought. The super(12)CO and super(13)CO lines saturate above 4 and 5 mag, respectively, whereas C super(18)O(1-0) never saturates in the whole image range probed by our study (up to 10 mag). Approximately 60% of the positions with super(12)CO(1-0) emission have subthermally excited lines, and almost all positions have excitation temperatures below the dust temperature. PDR models, using the Meudon code, can explain the super(12)CO(1-0) and super(13)CO(1- 0) emission with densities ranging between 10 super(3) and 10 super(4) cm super(-3). In general, local variations in the volume density and nonthermal motions (linked to different star formation activity) can explain the observations. Higher densities are needed to reproduce CO data toward active star-forming sites, such as NGC 1333, where the larger internal motions driven by the young protostars allow more photons from the embedded high-density cores to escape the cloud. In the most quiescent region, B5, the super(12)CO and super(13)CO emission appears to arise from an almost uniform thin layer of molecular material at densities around 10 super(4) cm super(-3), and in this region the integrated intensities of the two CO isotopologues are the lowest in the whole complex.
ABSTRACT
We study the physical drivers of slow molecular cloud mergers within a simulation of a Milky Way-like galaxy in the moving-mesh code arepo, and determine the influence of these mergers on ...the mass distribution and star formation efficiency of the galactic cloud population. We find that 83 per cent of these mergers occur at a relative velocity below 5 km s−1, and are associated with large-scale atomic gas flows, driven primarily by expanding bubbles of hot, ionized gas caused by supernova explosions and galactic rotation. The major effect of these mergers is to aggregate molecular mass into higher-mass clouds: mergers account for over 50 per cent of the molecular mass contained in clouds of mass M > 2 × 106 M⊙. These high-mass clouds have higher densities, internal velocity dispersions and instantaneous star formation efficiencies than their unmerged, lower mass precursors. As such, the mean instantaneous star formation efficiency in our simulated galaxy, with its merger rate of just 1 per cent of clouds per Myr, is 25 per cent higher than in a similar population of clouds containing no mergers.
We present an overview of the first data release (DR1) and first-look science from the Green Bank Ammonia Survey (GAS). GAS is a Large Program at the Green Bank Telescope to map all Gould Belt ...star-forming regions with mag visible from the northern hemisphere in emission from NH3 and other key molecular tracers. This first release includes the data for four regions in the Gould Belt clouds: B18 in Taurus, NGC 1333 in Perseus, L1688 in Ophiuchus, and Orion A North in Orion. We compare the NH3 emission to dust continuum emission from Herschel and find that the two tracers correspond closely. We find that NH3 is present in over 60% of the lines of sight with mag in three of the four DR1 regions, in agreement with expectations from previous observations. The sole exception is B18, where NH3 is detected toward ∼40% of the lines of sight with mag. Moreover, we find that the NH3 emission is generally extended beyond the typical 0.1 pc length scales of dense cores. We produce maps of the gas kinematics, temperature, and NH3 column densities through forward modeling of the hyperfine structure of the NH3 (1, 1) and (2, 2) lines. We show that the NH3 velocity dispersion, , and gas kinetic temperature, TK, vary systematically between the regions included in this release, with an increase in both the mean value and the spread of and TK with increasing star formation activity. The data presented in this paper are publicly available (https://dataverse.harvard.edu/dataverse/GAS_DR1).