Ammonia is a classic interstellar thermometer. The estimation of the rotational and kinetic temperatures can be affected by the blended Hyperfine Components (HFCs). We developed a new recipe, referred to as the HyperFine Group Ratio (HFGR), which utilizes only direct observables, namely the intensity ratios between the grouped HFCs. As tested on the model spectra, the empirical formulae in HFGR can derive the rotational temperature ($T_{\rm rot}$) from the HFC group ratios in an unambiguous manner. We compared HFGR to the two conventional methods, hyperfine fitting and line intensity ratio, based on both real data and simulated spectra. The HFGR has three major improvements. First, HFGR does not require modeling the HFC or fitting the line profiles, thus is more robust against the effect of HFC blending. Second, the simulation-enabled empirical formulae are much faster than fitting the spectra over the parameter space, so the computer time and human time can be both largely saved. Third, the statistical uncertainty of the temperature $\Delta T_{\rm rot}$ as a function of the signal-to-noise ratio (SNR) are also provided. HFGR can keep an internal error of $\Delta T_{\rm rot}\leq0.5$ K over a broad parameter space of temperature (10 to 70 K), line width (0.3 to 4 \kms), and optical depth (0 to 5). When applied to the noisy spectra, HFGR can keep an uncertainty of $\Delta T_{\rm rot}\leq 1.0$ K (1 $\sigma$) when SNR$>4$.

We present an overview of the CMZoom survey and its first data release. CMZoom is the first blind, high-resolution survey of the Central Molecular Zone (CMZ; the inner 500 pc of the Milky Way) at wavelengths sensitive to the pre-cursors of high-mass stars. CMZoom is a 500-hour Large Program on the Submillimeter Array (SMA) that mapped at 1.3 mm all of the gas and dust in the CMZ above a molecular hydrogen column density of 10^23 cm^-2 at a resolution of ~3" (0.1 pc). In this paper, we focus on the 1.3 mm dust continuum and its data release, but also describe CMZoom spectral line data which will be released in a forthcoming publication. While CMZoom detected many regions with rich and complex substructure, its key result is an overall deficit in compact substructures on 0.1 - 2 pc scales (the compact dense gas fraction: CDGF). In comparison with clouds in the Galactic disk, the CDGF in the CMZ is substantially lower, despite having much higher average column densities. CMZ clouds with high CDGFs are well-known sites of active star formation. The inability of most gas in the CMZ to form compact substructures is likely responsible for the dearth of star formation in the CMZ, surprising considering its high density. The factors responsible for the low CDGF are not yet understood but are plausibly due to the extreme environment of the CMZ, having far-reaching ramifications for our understanding of the star formation process across the cosmos.

We report results of a project to map HCN and HCO+ J = 1-0 emission toward a sample of molecular clouds in the inner Galaxy, all containing dense clumps that are actively engaged in star formation. We compare these two molecular line tracers with millimeter continuum emission and extinction, as inferred from 13CO, as tracers of dense gas in molecular clouds. The fraction of the line luminosity from each tracer that comes from the dense gas, as measured by AV > 8 mag, varies substantially from cloud to cloud. In all cases, a substantial fraction (in most cases, the majority) of the total luminosity arises in gas below the AV > 8 mag threshold and outside the region of strong mm continuum emission. Measurements of the luminosity of HCN toward other galaxies will likely be dominated by such gas at lower surface density. Substantial, even dominant, contributions to the total line luminosity can arise in gas with densities typical of the cloud as a whole (densities about 100 per cubic cm). Defining the dense clump from the HCN or HCO+ emission itself, similarly to previous studies, leads to a wide range of clump properties, with some being considerably larger and less dense than in previous studies. HCN and HCO+ have similar ability to trace dense gas for the clouds in this sample. For the two clouds with low virial parameters, the 13CO is definitely a worse tracer of the dense gas, but for the other four, it is equally good (or bad) at tracing dense gas.

Aims: Our aim is to identify the dominant molecular cooling lines and characteristic emission features in the 1.3 mm window of distinct regions in the northern part of the Orion A molecular cloud. By defining and analysing template regions, we also intend to help with the interpretation of observations from more distant sources which cannot be easily spatially resolved. Methods: We analyse an imaging line survey covering the area of OMC-1 to OMC-3 from 200.2 to 281.8 GHz obtained with the PI230 receiver at the APEX telescope. Masks are used to define regions with distinct properties (e.g. column density or temperature ranges) from which we obtain averaged spectra. Lines of 29 molecular species (55 isotopologues) are fitted for each region to obtain the respective total intensity. Results: We find that strong sources like Orion KL have a clear impact on the emission on larger scales. Although not spatially extended, their line emission contributes substantially to spectra averaged over large regions. Conversely, the emission signatures of dense, cold regions like OMC-2 and OMC-3 (e.g. enhanced N$_{2}$H$^{+}$ emission and low HCN/HNC ratio) seem to be difficult to pick up on larger scales, where they are eclipsed by signatures of stronger sources. In all regions, HCO$^{+}$ appears to contribute between 3% and 6% to the total intensity, the most stable value for all bright species. N$_{2}$H$^{+}$ shows the strongest correlation with column density, but not with typical high-density tracers like HCN, HCO$^{+}$, H$_{2}$CO, or HNC. Common line ratios associated with UV illumination, CN/HNC and CN/HCO$^{+}$, show ambiguous results on larger scales, suggesting that the identification of UV illuminated material may be more challenging. The HCN/HNC ratio may be related to temperature over varying scales.

We present new observations of C-band continuum emission and masers to assess high-mass ($>$8 $M_\odot$) star formation at early evolutionary phases in the inner 200 pc of the Central Molecular Zone (CMZ) of the Galaxy. The continuum observation is complete to free-free emission from stars above 10-11 $M_\odot$ in 91% of the covered area. We identify 104 compact sources in the continuum emission, among which five are confirmed ultracompact H II regions, 12 are candidates of ultra-compact H II regions, and the remaining 87 sources are mostly massive stars in clusters, field stars, evolved stars, pulsars, extragalactic sources, or of unknown nature that is to be investigated. We detect class II CH$_3$OH masers at 23 positions, among which six are new detections. We confirm six known H$_2$CO masers in two high-mass star forming regions, and detect two new H$_2$CO masers toward the Sgr C cloud, making it the ninth region in the Galaxy that contains masers of this type. In spite of these detections, we find that current high-mass star formation in the inner CMZ is only taking place in seven isolated clouds. The results suggest that star formation at early evolutionary phases in the CMZ is about 10 times less efficient than expected by the dense gas star formation relation, which is in line with previous studies that focus on more evolved phases of star formation. This means that if there will be any impending, next burst of star formation in the CMZ, it has not yet begun.

The Event Horizon Telescope (EHT) recently produced the first horizon-scale image of a supermassive black hole. Expanding the array to include a 3-meter space telescope operating at >200 GHz enables mass measurements of many black holes, movies of black hole accretion flows, and new tests of general relativity that are impossible from the ground.

We have investigated the formation and kinematics of sub-mm continuum cores in the Orion A molecular cloud. A comparison between sub-mm continuum and near infrared extinction shows a continuum core detection threshold of $A_V\sim$ 5-10 mag. The threshold is similar to the star formation extinction threshold of $A_V\sim$ 7 mag proposed by recent work, suggesting a universal star formation extinction threshold among clouds within 500 pc to the Sun. A comparison between the Orion A cloud and a massive infrared dark cloud G28.37+0.07 indicates that Orion A produces more dense gas within the extinction range 15 mag $\lesssim A_V \lesssim$ 60 mag. Using data from the CARMA-NRO Orion Survey, we find that dense cores in the integral-shaped filament (ISF) show sub-sonic core-to-envelope velocity dispersion that is significantly less than the local envelope line dispersion, similar to what has been found in nearby clouds. Dynamical analysis indicates that the cores are bound to the ISF. An oscillatory core-to-envelope motion is detected along the ISF. Its origin is to be further explored.

Understanding the collapse of clouds and the formation of protoplanetary disks is essential to understanding the formation of stars and planets. Infall and accretion, the mass-aggregation processes that occur at envelope and disk scales, drive the dynamical evolution of protostars. While the observations of protostars at different stages constrain their evolutionary tracks, the impact of variability due to accretion bursts on dynamical and chemical evolution of the source is largely unknown. The lasting effects on protostellar envelopes and disks are tracked through multi-wavelength and time domain observational campaigns, requiring deep X-ray, infrared, and radio imaging and spectroscopy, at a sufficient level of spatial detail to distinguish contributions from the various substructures (i.e., envelope from disk from star from outflow). Protostellar models derived from these campaigns will illuminate the initial chemical state of protoplanetary disks during the epoch of giant planet formation. Insight from individual star formation in the Milky Way is also necessary to understand star formation rates in extragalactic sources. This cannot be achieved with ground-based observatories and is not covered by currently approved instrumentation. Requirements: High (v < 10 km/s for survey; v < 1 km/s for followup) spectral resolution capabilities with relatively rapid response times in the IR (3-500 um), X-ray (0.1-10 keV), and radio (cm) are critical to follow the course of accretion and outflow during an outburst. Complementary, AU-scale radio observations are needed to probe the disk accretion zone, and 10 AU-scale to probe chemical and kinematic structures of the disk-forming regions, and track changes in the dust, ice, and gas within protostellar envelopes.

The cycling of baryons in and out of galaxies is what ultimately drives galaxy formation and evolution. The circumgalactic medium (CGM) represents the interface between the interstellar medium and the cosmic web, hence its properties are directly shaped by the baryon cycle. Although traditionally the CGM is thought to consist of warm and hot gas, recent breakthroughs are presenting a new scenario according to which an important fraction of its mass may reside in the cold atomic and molecular phase. This would represent fuel that is readily available for star formation, with crucial implications for feeding and feedback processes in galaxies. However, such cold CGM, especially in local galaxies where its projected size on sky is expected to be of several arcminutes, cannot be imaged by ALMA due to interferometric spatial scale filtering of large-scale structures. We show that the only way to probe the multiphase CGM including its coldest component is through a large (e.g. 50-m) single dish (sub-)mm telescope.

The past two decades have seen extensive surveys of the far-infrared to submillimeter continuum emission in the plane of our Galaxy. We line out prospects for the coming decade for corresponding molecular and atomic line surveys which are needed to fully understand the formation of the dense structures that give birth to clusters and stars out of the diffuse interstellar medium. We propose to work towards Galaxy wide surveys in mid-J CO lines to trace shocks from colliding clouds, Galaxy-wide surveys for atomic Carbon lines in order to get a detailed understanding of the relation of atomic and molecular gas in clouds, and to perform extensive surveys of the structure of the dense parts of molecular clouds to understand the importance of filaments/fibers over the full range of Galactic environments and to study how dense cloud cores are formed from the filaments. This work will require a large (50m) Single Dish submillimeter telescope equipped with massively multipixel spectrometer arrays, such as envisaged by the AtLAST project.

We investigate star formation at very early evolutionary phases in five massive clouds in the inner 500 pc of the Galaxy, the Central Molecular Zone. Using interferometer observations of H$_2$O masers and ultra-compact H II regions, we find evidence of ongoing star formation embedded in cores of 0.2 pc scales and $\gtrsim$10$^5$ cm$^{-3}$ densities. Among the five clouds, Sgr C possesses a high (9%) fraction of gas mass in gravitationally bound and/or protostellar cores, and follows the dense ($\gtrsim$10$^4$ cm$^{-3}$) gas star formation relation that is extrapolated from nearby clouds. The other four clouds have less than 1% of their cloud masses in gravitationally bound and/or protostellar cores, and star formation rates 10 times lower than predicted by the dense gas star formation relation. At the spatial scale of these cores, the star formation efficiency is comparable to that in Galactic disk sources. We suggest that the overall inactive star formation in these Central Molecular Zone clouds could be because there is much less gas confined in gravitationally bound cores, which may be a result of the strong turbulence in this region and/or the very early evolutionary stage of the clouds when collapse has only recently started.

In this paper we provide a comprehensive description of the internal dynamics of G0.253+0.016 (a.k.a. 'the Brick'); one of the most massive and dense molecular clouds in the Galaxy to lack signatures of widespread star formation. As a potential host to a future generation of high-mass stars, understanding largely quiescent molecular clouds like G0.253+0.016 is of critical importance. In this paper, we reanalyse Atacama Large Millimeter Array cycle 0 HNCO $J=4(0,4)-3(0,3)$ data at 3 mm, using two new pieces of software which we make available to the community. First, scousepy, a Python implementation of the spectral line fitting algorithm scouse. Secondly, acorns (Agglomerative Clustering for ORganising Nested Structures), a hierarchical n-dimensional clustering algorithm designed for use with discrete spectroscopic data. Together, these tools provide an unbiased measurement of the line of sight velocity dispersion in this cloud, $\sigma_{v_{los}, {\rm 1D}}=4.4\pm2.1$ kms$^{-1}$, which is somewhat larger than predicted by velocity dispersion-size relations for the Central Molecular Zone (CMZ). The dispersion of centroid velocities in the plane of the sky are comparable, yielding $\sigma_{v_{los}, {\rm 1D}}/\sigma_{v_{pos}, {\rm 1D}}\sim1.2\pm0.3$. This isotropy may indicate that the line-of-sight extent of the cloud is approximately equivalent to that in the plane of the sky. Combining our kinematic decomposition with radiative transfer modelling we conclude that G0.253+0.016 is not a single, coherent, and centrally-condensed molecular cloud; 'the Brick' is not a \emph{brick}. Instead, G0.253+0.016 is a dynamically complex and hierarchically-structured molecular cloud whose morphology is consistent with the influence of the orbital dynamics and shear in the CMZ.

The infrared dark clouds (IRDCs) G11.11$-$0.12 and G28.34$+$0.06 are two of the best-studied IRDCs in our Galaxy. These two clouds host clumps at different stages of evolution, including a massive dense clump in both clouds that is dark even at 70 and 100$\mu$m. Such seemingly quiescent massive dense clumps have been speculated to harbor cores that are precursors of high-mass stars and clusters. We observed these two "prestellar" regions at 1mm with the Submillimeter Array (SMA) with the aim of characterizing the nature of such cores. We show that the clumps fragment into several low- to high-mass cores within the filamentary structure of the enveloping cloud. However, while the overall physical properties of the clump may indicate a starless phase, we find that both regions host multiple outflows. The most massive core though 70 $\mu$m dark in both clumps is clearly associated with compact outflows. Such low-luminosity, massive cores are potentially the earliest stage in the evolution of a massive protostar. We also identify several outflow features distributed in the large environment around the most massive core. We infer that these outflows are being powered by young, low-mass protostars whose core mass is below our detection limit. These findings suggest that low-mass protostars have already formed or are coevally formed at the earliest phase of high-mass star formation.

We present an initial overview of the filamentary structure in the Orion A molecular cloud utilizing a high angular and velocity resolution C$^{18}$O(1-0) emission map that was recently produced as part of the CARMA-NRO Orion Survey. The main goal of this study is to build a credible method to study varying widths of filaments which has previously been linked to star formation in molecular clouds. Due to the diverse star forming activities taking place throughout its $\sim$20 pc length, together with its proximity of 388 pc, the Orion A molecular cloud provides an excellent laboratory for such an experiment to be carried out with high resolution and high sensitivity. Using the widely-known structure identification algorithm, DisPerSE, on a 3-dimensional (PPV) C$^{18}$O cube, we identified 625 relatively short (the longest being 1.74 pc) filaments over the entire cloud. We study the distribution of filament widths using FilChaP, a python package that we have developed and made publicly available. We find that the filaments identified in a 2 square degree PPV cube do not overlap spatially, except for the complex OMC-4 region that shows distinct velocity components along the line of sight. The filament widths vary between 0.02 and 0.3 pc depending on the amount of substructure that a filament possesses. The more substructure a filament has, the larger is its width. We also find that despite this variation, the filament width shows no anticorrelation with the central column density which is in agreement with previous Herschel observations.

Cosmic-rays are the primary initiators of interstellar chemistry, and getting a better understanding of the varying impact they have on the chemistry of interstellar clouds throughout the Milky Way will not only expand our understanding of interstellar medium chemistry in our own galaxy, but also aid in extra-galactic studies. This work uses the ALCHEMIC astrochemical modeling code to perform numerical simulations of chemistry for a range of ionization rates. We study the impact of variations in the cosmic ray ionization rate on molecular abundances under idealized conditions, given by constant temperatures and a fixed density of 1e4 cm-3. As part of this study we examine whether observations of molecular abundances can be used to infer the cosmic ray ionization rate in such a simplified case. We find that intense cosmic-ray ionisation results in molecules, in particular the large and complex ones, being largely dissociated, and the medium becoming increasingly atomic. Individual species have limitations in their use as probes of the cosmic ray ionization rate. At early time (<1 Myrs) ions such as N2H+ and HOC+ make the best probes, while at later times, neutral species such as HNCO and SO stand out, in particular due to their large abundance variations. It is however by combining species into pairs that we find the best probes. Molecular ions such as N2H+ combined with different neutral species can provide probe candidates that outmatch individual species, in particular N2H+/C4H, N2H+/C2H, HOC+/O and HOC+/HNCO. These still have limitations to their functional range, but are more functional as probes than individual species previously used.

The dynamics of massive clumps, the environment where massive stars originate, is still unclear. Many theories predict that these regions are in a state of near-virial equilibrium, or near energy equi-partition, while others predict that clumps are in a sub-virial state. Observationally, the majority of the massive clumps are in a sub-virial state with a clear anti-correlation between the virial parameter $\alpha_{vir}$ and the mass of the clumps $M_{c}$, which suggests that the more massive objects are also the more gravitationally bound. Although this trend is observed at all scales, from massive clouds down to star-forming cores, theories do not predict it. In this work we show how, starting from virialized clumps, an observational bias is introduced in the specific case where the kinetic and the gravitational energies are estimated in different volumes within clumps and how it can contribute to the spurious $\alpha_{vir}-M_{c}$ anti-correlation in these data. As a result, the observed effective virial parameter $\tilde{\alpha}_{eff}<\alpha_{vir}$, and in some circumstances it might not be representative of the virial state of the observed clumps.

Analysing the kinematics of filamentary molecular clouds is a crucial step towards understanding their role in the star formation process. Therefore, we study the kinematics of 283 filament candidates in the inner Galaxy, that were previously identified in the ATLASGAL dust continuum data. The $^{13}$CO(2 - 1) and C$^{18}$O(2 - 1) data of the SEDIGISM survey (Structure, Excitation, and Dynamics of the Inner Galactic Inter Stellar Medium) allows us to analyse the kinematics of these targets and to determine their physical properties at a resolution of 30 arcsec and 0.25 km/s. To do so, we developed an automated algorithm to identify all velocity components along the line-of-sight correlated with the ATLASGAL dust emission, and derive size, mass, and kinematic properties for all velocity components. We find two-third of the filament candidates are coherent structures in position-position-velocity space. The remaining candidates appear to be the result of a superposition of two or three filamentary structures along the line-of-sight. At the resolution of the data, on average the filaments are in agreement with Plummer-like radial density profiles with a power-law exponent of p = 1.5 +- 0.5, indicating that they are typically embedded in a molecular cloud and do not have a well-defined outer radius. Also, we find a correlation between the observed mass per unit length and the velocity dispersion of the filament of $m \sim \sigma_v^2$. We show that this relation can be explained by a virial balance between self-gravity and pressure. Another possible explanation could be radial collapse of the filament, where we can exclude infall motions close to the free-fall velocity.

We present the first results from a new, high resolution, $^{12}$CO(1-0), $^{13}$CO(1-0), and C$^{18}$O(1-0) molecular line survey of the Orion A cloud, hereafter referred to as the CARMA-NRO Orion Survey. CARMA observations have been combined with single-dish data from the Nobeyama 45m telescope to provide extended images at about 0.01 pc resolution, with a dynamic range of approximately 1200 in spatial scale. Here we describe the practical details of the data combination in uv space, including flux scale matching, the conversion of single dish data to visibilities, and joint deconvolution of single dish and interferometric data. A $\Delta$-variance analysis indicates that no artifacts are caused by combining data from the two instruments. Initial analysis of the data cubes, including moment maps, average spectra, channel maps, position-velocity diagrams, excitation temperature, column density, and line ratio maps provides evidence of complex and interesting structures such as filaments, bipolar outflows, shells, bubbles, and photo-eroded pillars. The implications for star formation processes are profound and follow-up scientific studies by the CARMA-NRO Orion team are now underway. We plan to make all the data products described here generally accessible; some are already available at https://dataverse.harvard.edu/dataverse/CARMA-NRO-Orion

A brief overview of recent advances in the study of star formation in the Galactic Center (GC) environment is presented. Particular attention is paid to new insights concerning the suppression of star formation in GC molecular clouds. Another focus is the question whether the GC can be used as a template for the understanding of starburst galaxies in the nearby and distant universe: this must be done with care. Some of the particular conditions in the center of the Milky Way do not necessarily play a role in starburst galaxies.

SDSS-V will be an all-sky, multi-epoch spectroscopic survey of over six million objects. It is designed to decode the history of the Milky Way, trace the emergence of the chemical elements, reveal the inner workings of stars, and investigate the origin of planets. It will also create an integral-field spectroscopic map of the gas in the Galaxy and the Local Group that is 1,000x larger than the current state of the art and at high enough spatial resolution to reveal the self-regulation mechanisms of galactic ecosystems. SDSS-V will pioneer systematic, spectroscopic monitoring across the whole sky, revealing changes on timescales from 20 minutes to 20 years. The survey will thus track the flickers, flares, and radical transformations of the most luminous persistent objects in the universe: massive black holes growing at the centers of galaxies. The scope and flexibility of SDSS-V will be unique among extant and future spectroscopic surveys: it is all-sky, with matched survey infrastructures in both hemispheres; it provides near-IR and optical multi-object fiber spectroscopy that is rapidly reconfigurable to serve high target densities, targets of opportunity, and time-domain monitoring; and it provides optical, ultra-wide-field integral field spectroscopy. SDSS-V, with its programs anticipated to start in 2020, will be well-timed to multiply the scientific output from major space missions (e.g., TESS, Gaia, eROSITA) and ground-based projects. SDSS-V builds on the 25-year heritage of SDSS's advances in data analysis, collaboration infrastructure, and product deliverables. The project is now refining its science scope, optimizing the survey strategies, and developing new hardware that builds on the SDSS-IV infrastructure. We present here an overview of the current state of these developments as we seek to build our worldwide consortium of institutional and individual members.

The star formation rate in the Central Molecular Zone (CMZ) is an order of magnitude lower than predicted according to star formation relations that have been calibrated in the disc of our own and nearby galaxies. Understanding how and why star formation appears to be different in this region is crucial if we are to understand the environmental dependence of the star formation process. Here, we present the detection of a sample of high-mass cores in the CMZ's "dust ridge" that have been discovered with the Submillimeter Array as part of the CMZoom survey. These cores range in mass from ~ 50 - 2150 Msun within radii of 0.1 - 0.25 pc. All appear to be young (pre-UCHII), meaning that they are prime candidates for representing the initial conditions of high-mass stars and sub-clusters. We report that at least two of these cores ('c1' and 'e1') contain young, high-mass protostars. We compare all of the detected cores with high-mass cores in the Galactic disc and find that they are broadly similar in terms of their masses and sizes, despite being subjected to external pressures that are several orders of magnitude greater - ~ 10^8 K/cm^3, as opposed to ~ 10^5 K/cm^3. The fact that > 80% of these cores do not show any signs of star-forming activity in such a high-pressure environment leads us to conclude that this is further evidence for an increased critical density threshold for star formation in the CMZ due to turbulence.

The Survey of Water and Ammonia in the Galactic Center (SWAG) covers the Central Molecular Zone (CMZ) of the Milky Way at frequencies between 21.2 and 25.4 GHz obtained at the Australia Telescope Compact Array at $\sim 0.9$ pc spatial and $\sim 2.0$ km s$^{-1}$ spectral resolution. In this paper, we present data on the inner $\sim 250$ pc ($1.4^\circ$) between Sgr C and Sgr B2. We focus on the hyperfine structure of the metastable ammonia inversion lines (J,K) = (1,1) - (6,6) to derive column density, kinematics, opacity and kinetic gas temperature. In the CMZ molecular clouds, we find typical line widths of $8-16$ km s$^{-1}$ and extended regions of optically thick ($\tau > 1$) emission. Two components in kinetic temperature are detected at $25-50$ K and $60-100$ K, both being significantly hotter than dust temperatures throughout the CMZ. We discuss the physical state of the CMZ gas as traced by ammonia in the context of the orbital model by Kruijssen et al. (2015) that interprets the observed distribution as a stream of molecular clouds following an open eccentric orbit. This allows us to statistically investigate the time dependencies of gas temperature, column density and line width. We find heating rates between $\sim 50$ and $\sim 100$ K Myr$^{-1}$ along the stream orbit. No strong signs of time dependence are found for column density or line width. These quantities are likely dominated by cloud-to-cloud variations. Our results qualitatively match the predictions of the current model of tidal triggering of cloud collapse, orbital kinematics and the observation of an evolutionary sequence of increasing star formation activity with orbital phase.

(Abridged) Context: The assumption of a gas-to-dust mass ratio (\gamma) is a common approach to estimate the basic properties of molecular clouds, such as total mass and column density of molecular hydrogen, from (sub)mm continuum observations of the dust. In the Milky Way a single value is used at all galactocentric radii, independently of the observed metallicity gradients. Both models and extragalactic observations suggest that this quantity increases for decreasing metallicity Z, typical of the outer regions in disks, where fewer heavy elements are available to form dust grains. Aims: We aim to investigate the variation of the gas-to-dust ratio as a function of galactocentric radius and metallicity, to allow a more accurate characterisation of the quantity of molecular gas across the galactic disk, as derived from observations of the dust. Methods: Observations of the optically thin C\$^{18}\$O (2-1) transition were obtained with the APEX telescope for a sample of 23 massive and dense star-forming regions in the far outer Galaxy (galactocentric distance greater than 14 kpc). From the modelling of this line and of the spectral energy distribution of the selected clumps we computed the gas-to-dust ratio and compared it to that of well-studied sources from the ATLASGAL TOP100 sample in the inner galactic disk. Results: The gradient in gas-to-dust ratio is found to be 0.087 dex/kpc (or equivalently \gamma\ proportional to Z\$^{-1.4}\$). The dust-to-metal ratio, decreases with galactocentric radius, which is the most common situation also for external late-type galaxies. This suggests that grain growth dominates over destruction. The predicted gas-to-dust ratio is in excellent agreement with the estimates in Magellanic clouds, for the appropriate value of Z.

Emission from high-dipole moment molecules such as HCN allows determination of the density in molecular clouds, and is often considered to trace the "dense" gas available for star formation. We assess the importance of electron excitation in various environments. The ratio of the rate coefficients for electrons and H$_2$ molecules, $\simeq$10$^5$ for HCN, yields the requirements for electron excitation to be of practical importance if $n({\rm H}_2) \leq\ 10^{5.5} ~ \rm cm^{-3}$ and $X({\rm e}^-) \geq\ 10^{-5}$, where the numerical factors reflect critical values $n_{\rm{}c}({\rm H_2})$ and $X^*({\rm{}e}^-)$. This indicates that in regions where a large fraction of carbon is ionized, $X({\rm e}^-)$ will be large enough to make electron excitation significant. The situation is in general similar for other "high density tracers", including HCO$^+$, CN, and CS. But there are significant differences in the critical electron fractional abundance, $X^*({\rm e}^-)$, defined by the value required for equal effect from collisions with H$_2$ and e$^-$. Electron excitation is, for example, unimportant for CO and C$^+$. Electron excitation may be responsible for the surprisingly large spatial extent of the emission from dense gas tracers in some molecular clouds (Pety et al. 2017; Kauffmann, Goldsmith et al. 2017). The enhanced estimates for HCN abundances and HCN/CO and HCN/HCO$^+$ ratios observed in the nuclear regions of luminous galaxies may be in part a result of electron excitation of high dipole moment tracers. The importance of electron excitation will depend on detailed models of the chemistry, which may well be non-steady state and non-static.

Trends observed in galaxies, such as the Gao \& Solomon relation, suggest a linear relation between the star formation rate and the mass of dense gas available for star formation. Validation of such relations requires the establishment of reliable methods to trace the dense gas in galaxies. One frequent assumption is that the HCN ($J=1$--0) transition is unambiguously associated with gas at $\rm{}H_2$ densities $\gg{}10^4~\rm{}cm^{-3}$. If so, the mass of gas at densities $\gg{}10^4~\rm{}cm^{-3}$ could be inferred from the luminosity of this emission line, $L_{\rm{}HCN\,(1\text{--}0)}$. Here we use observations of the Orion~A molecular cloud to show that the HCN ($J=1$--0) line traces much lower densities $\sim{}10^3~\rm{}cm^{-3}$ in cold sections of this molecular cloud, corresponding to visual extinctions $A_V\approx{}6~\rm{}mag$. We also find that cold and dense gas in a cloud like Orion produces too little HCN emission to explain $L_{\rm{}HCN\,(1\text{--}0)}$ in star--forming galaxies, suggesting that galaxies might contain a hitherto unknown source of HCN emission. In our sample of molecules observed at frequencies near 100~GHz (also including $\rm{}^{12}CO$, $\rm{}^{13}CO$, $\rm{}C^{18}O$, CN, and CCH), $\rm{}N_2H^+$ is the only species clearly associated with rather dense gas.

We present the first systematic study of the density structure of clouds found in a complete sample covering all major molecular clouds in the Central Molecular Zone (CMZ; inner $\sim{}200~\rm{}pc$) of the Milky Way. This is made possible by using data from the Galactic Center Molecular Cloud Survey (GCMS), the first study resolving all major molecular clouds in the CMZ at interferometer angular resolution. We find that many CMZ molecular clouds have unusually shallow density gradients compared to regions elsewhere in the Milky Way. This is possibly a consequence of weak gravitational binding of the clouds. The resulting relative absence of dense gas on spatial scales $\sim{}0.1~\rm{}pc$ is probably one of the reasons why star formation (SF) in dense gas of the CMZ is suppressed by a factor $\sim{}10$, compared to solar neighborhood clouds. Another factor suppressing star formation are the high SF density thresholds that likely result from the observed gas kinematics. Further, it is possible but not certain that the star formation activity and the cloud density structure evolve systematically as clouds orbit the CMZ.

We recently reported a population of protostellar candidates in the 20 km s$^{-1}$ cloud in the Central Molecular Zone of the Milky Way, traced by H$_2$O masers in gravitationally bound dense cores. In this paper, we report high-angular-resolution ($\sim$3'') molecular line studies of the environment of star formation in this cloud. Maps of various molecular line transitions as well as the continuum at 1.3 mm are obtained using the Submillimeter Array. Five NH$_3$ inversion lines and the 1.3 cm continuum are observed with the Karl G. Jansky Very Large Array. The interferometric observations are complemented with single-dish data. We find that the CH$_3$OH, SO, and HNCO lines, which are usually shock tracers, are better correlated spatially with the compact dust emission from dense cores among the detected lines. These lines also show enhancement in intensities with respect to SiO intensities toward the compact dust emission, suggesting the presence of slow shocks or hot cores in these regions. We find gas temperatures of $\gtrsim$100 K at 0.1-pc scales based on RADEX modelling of the H$_2$CO and NH$_3$ lines. Although no strong correlations between temperatures and linewidths/H$_2$O maser luminosities are found, in high-angular-resolution maps we notice several candidate shock heated regions offset from any dense cores, as well as signatures of localized heating by protostars in several dense cores. Our findings suggest that at 0.1-pc scales in this cloud star formation and strong turbulence may together affect the chemistry and temperature of the molecular gas.

The origin and life-cycle of molecular clouds are still poorly constrained, despite their importance for understanding the evolution of the interstellar medium. We have carried out a systematic, homogeneous, spectroscopic survey of the inner Galactic plane, in order to complement the many continuum Galactic surveys available with crucial distance and gas-kinematic information. Our aim is to combine this data set with recent infrared to sub-millimetre surveys at similar angular resolutions. The SEDIGISM survey covers 78 deg^2 of the inner Galaxy (-60 deg < l < +18 deg, |b| < 0.5 deg) in the J=2-1 rotational transition of 13CO. This isotopologue of CO is less abundant than 12CO by factors up to 100. Therefore, its emission has low to moderate optical depths, and higher critical density, making it an ideal tracer of the cold, dense interstellar medium. The data have been observed with the SHFI single-pixel instrument at APEX. The observational setup covers the 13CO(2-1) and C18O(2-1) lines, plus several transitions from other molecules. The observations have been completed. Data reduction is in progress, and the final data products will be made available in the near future. Here we give a detailed description of the survey and the dedicated data reduction pipeline. Preliminary results based on a science demonstration field covering -20 deg < l < -18.5 deg are presented. Analysis of the 13CO(2-1) data in this field reveals compact clumps, diffuse clouds, and filamentary structures at a range of heliocentric distances. By combining our data with data in the (1-0) transition of CO isotopologues from the ThrUMMS survey, we are able to compute a 3D realization of the excitation temperature and optical depth in the interstellar medium. Ultimately, this survey will provide a detailed, global view of the inner Galactic interstellar medium at an unprecedented angular resolution of ~30".

Research on Galactic Center star formation is making great advances, in particular due to new data from interferometers spatially resolving molecular clouds in this environment. These new results are discussed in the context of established knowledge about the Galactic Center. Particular attention is paid to suppressed star formation in the Galactic Center and how it might result from shallow density gradients in molecular clouds.

The Central Molecular Zone (CMZ; inner $\sim{}200~\rm{}pc$) of the Milky Way is a star formation (SF) environment with very extreme physical properties. Exploration of SF in this region is important because (i) this region allows us to test models of star formation under exceptional conditions, and (ii) the CMZ clouds might be suitable to serve as templates to understand the physics of starburst galaxies in the nearby and the distant universe. For this reason we launched the Galactic Center Molecular Cloud Survey (GCMS), the first systematic study that resolves all major CMZ clouds at interferometer angular resolution (i.e., a few arc seconds). Here we present initial results based on observations with the Submillimeter Array (SMA) and the Atacama Pathfinder Experiment (APEX). Our study is complemented by dust emission data from the Herschel Space Telescope and a comprehensive literature survey of CMZ star formation activity. Our research reveals (i) an unusually steep linewidth-size relation, $\sigma(v)\propto{}r_{\rm{}eff}^{0.66\pm{}0.18}$, down to velocity dispersions $\sim{}0.6~\rm{}km\,s^{-1}$ at 0.1 pc scale. This scaling law potentially results from the decay of gas motions to transonic velocities in strong shocks. The data also show that, relative to dense gas in the solar neighborhood, (ii) star formation is suppressed by factors $\gtrsim{}10$ in individual CMZ clouds. This observation encourages exploration of processes that can suppress SF inside dense clouds for a significant period of time.

The ATLASGAL survey provides an ideal basis for detailed studies of large numbers of massive star forming clumps covering the whole range of evolutionary stages. The ATLASGAL Top100 is a sample of clumps selected from their infrared and radio properties to be representative for the whole range of evolutionary stages. The ATLASGAL Top100 sources are the focus of a number of detailed follow-up studies that will be presented in a series of papers. In the present work we use the dust continuum emission to constrain the physical properties of this sample and identify trends as a function of source evolution. We determine flux densities from mid-infrared to submm wavelength (8-870 micron) images and use these values to fit their spectral energy distributions (SEDs) and determine their dust temperature and flux. Combining these with recent distances from the literature including maser parallax measurements we determine clump masses, luminosities and column densities. We find trends for increasing temperature, luminosity and column density with the proposed evolution sequence, confirming that this sample is representative of different evolutionary stages of massive star formation. We show that most of the sample has the ability to form massive stars (including the most massive O-type stars) and that the majority is gravitationally unstable and hence likely to be collapsing. The highest column density ATLASGAL sources presented cover the whole range of evolutionary stages from the youngest to the most evolved high-mass star forming clumps. Their study provides a unique starting point for more in-depth research on massive star formation in four distinct evolutionary stages whose well defined physical parameters afford more detailed studies. As most of the sample is closer than 5 kpc, these sources are also ideal for follow-up observations with high spatial resolution.

The inner few hundred parsecs of the Milky Way, the Central Molecular Zone (CMZ), is our closest laboratory for understanding star formation in the extreme environments (hot, dense, turbulent gas) that once dominated the universe. We present an update on the first large-area survey to expose the sites of star formation across the CMZ at high-resolution in submillimeter wavelengths: the CMZoom survey with the Submillimeter Array (SMA). We identify the locations of dense cores and search for signatures of embedded star formation. CMZoom is a three-year survey in its final year and is mapping out the highest column density regions of the CMZ in dust continuum and a variety of spectral lines around 1.3 mm. CMZoom combines SMA compact and subcompact configurations with single-dish data from BGPS and the APEX telescope, achieving an angular resolution of about 4" (0.2 pc) and good image fidelity up to large spatial scales.

The Central Molecular Zone (CMZ), usually referring to the inner 500 pc of the Galaxy, contains a dozen of massive ($\sim10^5$ $M_\odot$) molecular clouds. Are these clouds going to actively form stars like Sgr B2? How are they affected by the extreme physical conditions in the CMZ, such as strong turbulence? Here we present a first step towards answering these questions. Using high-sensitivity, high angular resolution radio and (sub)millimeter observations, we studied deeply embedded star formation in six massive clouds in the CMZ, including the 20 and 50 km s$^{-1}$ clouds, Sgr B1 off (as known as dust ridge clouds e/f), Sgr C, Sgr D, and G0.253-0.016. The VLA water maser observations suggest a population of deeply embedded protostellar candidates, many of which are new detections. The SMA 1.3 mm continuum observations reveal peaks in dust emission associated with the masers, suggesting the existence of dense cores. While our findings confirm that clouds such as G0.253-0.016 lack internal compact substructures and are quiescent in terms of star formation, two clouds (the 20 km s$^{-1}$ cloud and Sgr C) stand out with clusters of water masers with associated dense cores which may suggest a population of deeply embedded protostars at early evolutionary phases. Follow-up observations with VLA and ALMA are necessary to confirm their protostellar nature.

The Central Molecular Zone (CMZ) at the center of our Galaxy is the best template to study star formation processes under extreme conditions, similar to those in high-redshift galaxies. We observed on-the-fly maps of para-H$_{2}$CO transitions at 218 GHz and 291 GHz towards seven Galactic Center clouds. From the temperature-sensitive integrated intensity line ratios of H$_{2}$CO(3$_{2,1}-$2$_{2,0}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$) and H$_{2}$CO(4$_{2,2}-$3$_{2,1}$)/H$_{2}$CO(4$_{0,4}-$3$_{0,3}$) in combination with radiative transfer models, we produce gas temperature maps of our targets. These transitions are sensitive to gas with densities of $\sim$10$^{5}$ cm$^{-3}$ and temperatures <150 K. The measured gas temperatures in our sources are all higher (>40 K) than their dust temperatures ($\sim$25 K). Our targets have a complex velocity structure that requires a careful disentanglement of the different components. We produce temperature maps for each of the velocity components and show that the temperatures of the components differ, revealing temperature gradients in the clouds. Combining the temperature measurements with the integrated intensity line ratio of H$_{2}$CO(4$_{0,4}-$3$_{0,3}$)/H$_{2}$CO(3$_{0,3}-$2$_{0,2}$), we constrain the density of this warm gas to 10$^{4}-$10$^{6}$ cm$^{-3}$. We find a positive correlation of the line width of the main H$_{2}$CO lines with the temperature of the gas, direct evidence for gas heating via turbulence. Our data is consistent with a turbulence heating model with a density of n = 10$^5$ cm$^{-3}$.

We present a 3 mm spectral line and continuum survey of L1451 in the Perseus Molecular Cloud. These observations are from the CARMA Large Area Star Formation Survey (CLASSy), which also imaged Barnard 1, NGC 1333, Serpens Main and Serpens South. L1451 is the survey region with the lowest level of star formation activity---it contains no confirmed protostars. HCO+, HCN, and N2H+ (J=1-0) are all detected throughout the region, with HCO+ the most spatially widespread, and molecular emission seen toward 90% of the area above N(H_2) column densities of 1.9x10^21 cm^-2. HCO+ has the broadest velocity dispersion, near 0.3 km/s on average, compared to ~0.15 km/s for the other molecules, thus representing a range from supersonic to subsonic gas motions. Our non-binary dendrogram analysis reveals that the dense gas traced by each molecule has similar hierarchical structure, and that gas surrounding the candidate first hydrostatic core (FHSC), L1451-mm, and other previously detected single-dish continuum clumps have similar hierarchical structure; this suggests that different sub-regions of L1451 are fragmenting on the pathway to forming young stars. We determined the three-dimensional morphology of the largest detectable dense gas structures to be relatively ellipsoidal compared to other CLASSy regions, which appeared more flattened at largest scales. A virial analysis shows the most centrally condensed dust structures are likely unstable against collapse. Additionally, we identify a new spherical, centrally condensed N2H+ feature that could be a new FHSC candidate. The overall results suggest L1451 is a young region starting to form its generation of stars within turbulent, hierarchical structures.

We report on the young massive clump (G35.20w) in W48 that previous molecular line and dust observations have revealed to be in the very early stages of star formation. Based on virial analysis, we find that a strong field of ~1640 microG is required to keep the clump in pressure equilibrium. We performed a deep Zeeman effect measurement of the 113 GHz CN (1-0) line towards this clump with the IRAM 30 m telescope. We combine simultaneous fitting of all CN hyperfines with Monte Carlo simulations for a large range in realization of the magnetic field to obtain a constraint on the line-of-sight field strength of -687 +/- 420 microG. We also analyze archival dust polarization observations towards G35.20w. A strong magnetic field is implied by the remarkably ordered field orientation that is perpendicular to the longest axis of the clump. Based on this, we also estimate the plane-of-sky component of the magnetic field to be ~740 microG. This allows for a unique comparison of the two orthogonal measurements of magnetic field strength of the same region and at similar spatial scales. The expected total field strength shows no significant conflict between the observed field and that required for pressure equilibrium. By producing a probability distribution for a large range in field geometries, we show that plane-of-sky projections are much closer to the true field strengths than line-of-sight projections. This can present a significant challenge for Zeeman measurements of magnetized structures, even with ALMA. We also show that CN molecule does not suffer from depletion on the observed scales in the predominantly cold and highly deuterated core in an early stage of high-mass star formation and is thus a good tracer of the dense gas.

We report the results of an HCO+ (3-2) and N2D+ (3-2) molecular line survey performed toward 91 dense cores in the Perseus molecular cloud using the James Clerk Maxwell Telescope, to identify the fraction of starless and protostellar cores with systematic radial motions. We quantify the HCO+ asymmetry using a dimensionless asymmetry parameter $\delta_v$, and identify 20 cores with significant blue or red line asymmetries in optically-thick emission indicative of collapsing or expanding motions, respectively. We separately fit the HCO+ profiles with an analytic collapse model and determine contraction (expansion) speeds toward 22 cores. Comparing the $\delta_v$ and collapse model results, we find that $\delta_v$ is a good tracer of core contraction if the optically-thin emission is aligned with the model-derived systemic velocity. The contraction speeds range from subsonic (0.03 km/s) to supersonic (0.4 km/s), where the supersonic contraction speeds may trace global rather than local core contraction. Most cores have contraction speeds significantly less than their free-fall speeds. Only 7 of 28 starless cores have spectra well-fit by the collapse model, which more than doubles (15 of 28) for protostellar cores. Starless cores with masses greater than the Jeans mass (M/M$_J$ > 1) are somewhat more likely to show contraction motions. We find no trend of optically-thin non-thermal line width with M/M$_J$, suggesting that any undetected contraction motions are small and subsonic. Most starless cores in Perseus are either not in a state of collapse or expansion, or are in a very early stage of collapse.

Stars form predominantly in groups usually denoted as clusters or associations. The observed stellar groups display a broad spectrum of masses, sizes and other properties, so it is often assumed that there is no underlying structure in this diversity. Here we show that the assumption of an unstructured multitude of cluster or association types might be misleading. Current data compilations of clusters show correlations between cluster mass, size, age, maximum stellar mass etc. In this first paper we take a closer look at the correlation of cluster mass and radius. We use literature data to explore relations in cluster and molecular core properties in the solar neighborhood. We show that for embedded clusters in the solar neighborhood there exists a clear correlation between cluster mass and half-mass radius of the form $M_c = C R_c^{\gamma}$ with gamma = 1.7 +/-0.2. This correlation holds for infra red K band data as well as X-ray sources and for clusters containing a hundred stars up to those consisting of a few tens of thousands of stars. The correlation is difficult to verify for clusters containing <30 stars due to low-number statistics. Dense clumps of gas are the progenitors of the embedded clusters. We find a similar slope for the mass-size relation of dense, massive clumps as for the embedded star clusters. This might point at a direct translation from gas to stellar mass: however, it is difficult to relate size measurements for clusters (stars) to those for gas profiles. Taking into account multiple paths for clump mass into cluster mass, we obtain an average star-formation efficiency of 18%{+9.3}{-5.7} for the embedded clusters in the solar neighborhood. The derived mass-radius relation gives constraints for the theory of clustered star formation. Analytical models and simulations of clustered star formation have to reproduce this relation in order to be realistic (abridged)

We have discovered a new H$_2$CO (formaldehyde) $1_{1,0}-1_{1,1}$ 4.82966 GHz maser in Galactic Center Cloud C, G0.38+0.04. At the time of acceptance, this is the eighth region containing an H$_2$CO maser detected in the Galaxy. Cloud C is one of only two sites of confirmed high-mass star formation along the Galactic Center Ridge, affirming that H$_2$CO masers are exclusively associated with high-mass star formation. This discovery led us to search for other masers, among which we found new SiO vibrationally excited masers, making this the fourth star-forming region in the Galaxy to exhibit SiO maser emission. Cloud C is also a known source of CH$_3$OH Class-II and OH maser emission. There are now two known SiO and H$_2$CO maser containing regions in the CMZ, compared to two and six respectively in the Galactic disk, while there is a relative dearth of H$_2$O and CH$_3$OH Class-II masers in the CMZ. SiO and H$_2$CO masers may be preferentially excited in the CMZ, perhaps due to higher gas-phase abundances from grain destruction and heating, or alternatively H$_2$O and CH$_3$OH maser formation may be suppressed in the CMZ. In any case, Cloud C is a new testing ground for understanding maser excitation conditions.

We report the discovery of a population of deeply embedded protostellar candidates in the 20 km s$^{-1}$ cloud, one of the massive molecular clouds in the Central Molecular Zone (CMZ) of the Milky Way, using interferometric submillimeter continuum and H$_2$O maser observations. The submillimeter continuum emission shows five 1-pc scale clumps, each of which further fragments into several 0.1-pc scale cores. We identify 17 dense cores, among which 12 are gravitationally bound. Among the 18 H$_2$O masers detected, 13 coincide with the cores and probably trace outflows emanating from the protostars. There are also 5 gravitationally bound dense cores without H$_2$O maser detection. In total the 13 masers and 5 cores may represent 18 protostars with spectral types later than B1 or potential growing more massive stars at earlier evolutionary stage, given the non-detection in the centimeter radio continuum. In combination with previous studies of CH$_3$OH masers, we conclude that the star formation in this cloud is at an early evolutionary phase, before the presence of any significant ionizing or heating sources. Our findings indicate that star formation in this cloud may be triggered by a tidal compression as it approaches pericenter, similar to the case of G0.253+0.016 but with a higher star formation rate, and demonstrate that high angular resolution, high sensitivity maser and submillimeter observations are a promising technique to unveil deeply embedded star formation in the CMZ.

The Galactic center is the closest region in which we can study star formation under extreme physical conditions like those in high-redshift galaxies. We measure the temperature of the dense gas in the central molecular zone (CMZ) and examine what drives it. We mapped the inner 300 pc of the CMZ in the temperature-sensitive J = 3-2 para-formaldehyde (p-H$_2$CO) transitions. We used the $3_{2,1} - 2_{2,0} / 3_{0,3} - 2_{0,2}$ line ratio to determine the gas temperature in $n \sim 10^4 - 10^5 $cm$^{-3}$ gas. We have produced temperature maps and cubes with 30" and 1 km/s resolution and published all data in FITS form. Dense gas temperatures in the Galactic center range from ~60 K to > 100 K in selected regions. The highest gas temperatures T_G > 100 K are observed around the Sgr B2 cores, in the extended Sgr B2 cloud, the 20 km/s and 50 km/s clouds, and in "The Brick" (G0.253+0.016). We infer an upper limit on the cosmic ray ionization rate ${\zeta}_{CR} < 10^{-14}$ 1/s. The dense molecular gas temperature of the region around our Galactic center is similar to values found in the central regions of other galaxies, in particular starburst systems. The gas temperature is uniformly higher than the dust temperature, confirming that dust is a coolant in the dense gas. Turbulent heating can readily explain the observed temperatures given the observed line widths. Cosmic rays cannot explain the observed variation in gas temperatures, so CMZ dense gas temperatures are not dominated by cosmic ray heating. The gas temperatures previously observed to be high in the inner ~75 pc are confirmed to be high in the entire CMZ.

We present an analysis of probability distribution functions (pdfs) of column density in different zones of the star-forming region Perseus and its diffuse environment based on the map of dust opacity at 353 GHz available from the Planck archive. The pdf shape can be fitted by a combination of a lognormal function and an extended power-law tail at high densities, in zones centred at the molecular cloud Perseus. A linear combination of several lognormals fits very well the pdf in rings surrounding the cloud or in zones of its diffuse neighbourhood. The slope of the mean density scaling law $\langle\rho\rangle_L \propto L^\alpha$ is steep ($\alpha=-1.93$) in the former case and rather shallow ($\alpha=-0.77\pm0.11$) in the rings delineated around the cloud. We interpret these findings as signatures of two distinct physical regimes: i) a gravoturbulent one which is characterized by nearly linear scaling of mass and practical lack of velocity scaling; and ii) a predominantly turbulent one which is best described by steep velocity scaling and by invariant for compressible turbulence $\langle\rho\rangle_L u_L^3/L$, describing a scale-independent flux of the kinetic energy per unit volume through turbulent cascade. The gravoturbulent spatial domain can be identified with the molecular cloud Perseus while a relatively sharp transition to predominantly turbulent regime occurs in its vicinity.

High-mass Stars are cosmic engines known to dominate the energetics in the Milky Way and other galaxies. However, their formation is still not well understood. Massive, cold, dense clouds, often appearing as Infrared Dark Clouds (IRDCs), are the nurseries of massive stars. No measurements of magnetic fields in IRDCs in a state prior to the onset of high-mass star formation (HMSF) have previously been available, and prevailing HMSF theories do not consider strong magnetic fields. Here, we report observations of magnetic fields in two of the most massive IRDCs in the Milky Way. We show that IRDCs G11.11-0.12 and G0.253+0.016 are strongly magnetized and that the strong magnetic field is as important as turbulence and gravity for HMSF. The main dense filament in G11.11-0.12 is perpendicular to the magnetic field, while the lower density filament merging onto the main filament is parallel to the magnetic field. The implied magnetic field is strong enough to suppress fragmentation sufficiently to allow HMSF. Other mechanisms reducing fragmentation, such as the entrapment of heating from young stars via high mass surface densities, are not required to facilitate HMSF.

We present observations of N2H+(1-0), HCO+(1-0), and HCN(1-0) toward the Serpens Main molecular cloud from the CARMA Large Area Star Formation Survey (CLASSy). We mapped 150 square arcminutes of Serpens Main with an angular resolution of 7 arcsecs. The gas emission is concentrated in two subclusters (the NW and SE subclusters). The SE subcluster has more prominent filamentary structures and more complicated kinematics compared to the NW subcluster. The majority of gas in the two subclusters has subsonic to sonic velocity dispersions. We applied a dendrogram technique with N2H+(1-0) to study the gas structures; the SE subcluster has a higher degree of hierarchy than the NW subcluster. Combining the dendrogram and line fitting analyses reveals two distinct relations: a flat relation between nonthermal velocity dispersion and size, and a positive correlation between variation in velocity centroids and size. The two relations imply a characteristic depth of 0.15 pc for the cloud. Furthermore, we have identified six filaments in the SE subcluster. These filaments have lengths of 0.2 pc and widths of 0.03 pc, which is smaller than a characteristic width of 0.1 pc suggested by Herschel observations. The filaments can be classified into two types based on their properties. The first type, located in the northeast of the SE subcluster, has larger velocity gradients, smaller masses, and nearly critical mass-per-unit-length ratios. The other type, located in the southwest of the SE subcluster, has the opposite properties. Several YSOs are formed along two filaments which have supercritical mass per unit length ratios, while filaments with nearly critical mass-per-unit-length ratios are not associated with YSOs, suggesting that stars are formed on gravitationally unstable filaments.

We present details of the CARMA Large Area Star Formation Survey (CLASSy), while focusing on observations of Barnard 1. CLASSy is a CARMA Key Project that spectrally imaged N2H+, HCO+, and HCN (J=1-0 transitions) across over 800 square arcminutes of the Perseus and Serpens Molecular Clouds. The observations have angular resolution near 7" and spectral resolution near 0.16 km/s. We imaged ~150 square arcminutes of Barnard 1, focusing on the main core, and the B1 Ridge and clumps to its southwest. N2H+ shows the strongest emission, with morphology similar to cool dust in the region, while HCO+ and HCN trace several molecular outflows from a collection of protostars in the main core. We identify a range of kinematic complexity, with N2H+ velocity dispersions ranging from ~0.05-0.50 km/s across the field. Simultaneous continuum mapping at 3 mm reveals six compact object detections, three of which are new detections. A new non-binary dendrogram algorithm is used to analyze dense gas structures in the N2H+ position-position-velocity (PPV) cube. The projected sizes of dendrogram-identified structures range from about 0.01-0.34 pc. Size-linewidth relations using those structures show that non-thermal line-of-sight velocity dispersion varies weakly with projected size, while rms variation in the centroid velocity rises steeply with projected size. Comparing these relations, we propose that all dense gas structures in Barnard 1 have comparable depths into the sky, around 0.1-0.2 pc; this suggests that over-dense, parsec-scale regions within molecular clouds are better described as flattened structures rather than spherical collections of gas. Science-ready PPV cubes for Barnard 1 molecular emission are available for download.

The very long and thin infrared dark cloud "Nessie" is even longer than had been previously claimed, and an analysis of its Galactic location suggests that it lies directly in the Milky Way's mid-plane, tracing out a highly elongated bone-like feature within the prominent Scutum-Centaurus spiral arm. Re-analysis of mid-infrared imagery from the Spitzer Space Telescope shows that this IRDC is at least 2, and possibly as many as 5 times longer than had originally been claimed by Nessie's discoverers (Jackson et al. 2010); its aspect ratio is therefore at least 300:1, and possibly as large as 800:1. A careful accounting for both the Sun's offset from the Galactic plane ($\sim 25$ pc) and the Galactic center's offset from the $(l^{II},b^{II})=(0,0)$ position shows that the latitude of the true Galactic mid-plane at the 3.1 kpc distance to the Scutum-Centaurus Arm is not $b=0$, but instead closer to $b=-0.4$, which is the latitude of Nessie to within a few pc. An analysis of the radial velocities of low-density (CO) and high-density (${\rm NH}_3$) gas associated with the Nessie dust feature suggests that Nessie runs along the Scutum-Centaurus Arm in position-position-velocity space, which means it likely forms a dense `spine' of the arm in real space as well. The Scutum-Centaurus arm is the closest major spiral arm to the Sun toward the inner Galaxy, and, at the longitude of Nessie, it is almost perpendicular to our line of sight, making Nessie the easiest feature to see as a shadow elongated along the Galactic Plane from our location. Future high-resolution dust mapping and molecular line observations of the harder-to-find Galactic "bones" should allow us to exploit the Sun's position above the plane to gain a (very foreshortened) view "from above" of the Milky Way's structure.

We present the N2H+(J=1-0) map of the Serpens South molecular cloud obtained as part of the CARMA Large Area Star Formation Survey (CLASSy). The observations cover 250 square arcminutes and fully sample structures from 3000 AU to 3 pc with a velocity resolution of 0.16 km/s, and they can be used to constrain the origin and evolution of molecular cloud filaments. The spatial distribution of the N2H+ emission is characterized by long filaments that resemble those observed in the dust continuum emission by Herschel. However, the gas filaments are typically narrower such that, in some cases, two or three quasi-parallel N2H+ filaments comprise a single observed dust continuum filament. The difference between the dust and gas filament widths casts doubt on Herschel ability to resolve the Serpens South filaments. Some molecular filaments show velocity gradients along their major axis, and two are characterized by a steep velocity gradient in the direction perpendicular to the filament axis. The observed velocity gradient along one of these filaments was previously postulated as evidence for mass infall toward the central cluster, but these kind of gradients can be interpreted as projection of large-scale turbulence.

The "Airships: A New Horizon for Science" study at the Keck Institute for Space Studies investigated the potential of a variety of airships currently operable or under development to serve as observatories and science instrumentation platforms for a range of space, atmospheric, and Earth science. The participants represent a diverse cross-section of the aerospace sector, NASA, and academia. Over the last two decades, there has been wide interest in developing a high altitude, stratospheric lighter-than-air (LTA) airship that could maneuver and remain in a desired geographic position (i.e., "station-keeping") for weeks, months or even years. Our study found considerable scientific value in both low altitude (< 40 kft) and high altitude (> 60 kft) airships across a wide spectrum of space, atmospheric, and Earth science programs. Over the course of the study period, we identified stratospheric tethered aerostats as a viable alternative to airships where station-keeping was valued over maneuverability. By opening up the sky and Earth's stratospheric horizon in affordable ways with long-term flexibility, airships allow us to push technology and science forward in a project-rich environment that complements existing space observatories as well as aircraft and high-altitude balloon missions.

Whether or not molecular clouds and embedded cloud fragments are stable against collapse is of utmost importance for the study of the star formation process. Only "supercritical" cloud fragments are able to collapse and form stars. The virial parameter, alpha=M_vir/M, which compares the virial to the actual mass, provides one way to gauge stability against collapse. Supercritical cloud fragments are characterized by alpha<2, as indicated by a comprehensive stability analysis considering perturbations in pressure and density gradients. Past research has suggested that virial parameters alpha>2 prevail in clouds. This would suggest that collapse towards star formation is a gradual and relatively slow process, and that magnetic fields are not needed to explain the observed cloud structure. Here, we review a range of very recent observational studies that derive virial parameters <<2 and compile a catalogue of 1325 virial parameter estimates. Low values of alpha are in particular observed for regions of high mass star formation (HMSF). These observations may argue for a more rapid and violent evolution during collapse. This would enable "competitive accretion" in HMSF, constrain some models of "monolithic collapse", and might explain the absence of high--mass starless cores. Alternatively, the data could point at the presence of significant magnetic fields ~1 mG at high gas densities. We examine to what extent the derived observational properties might be biased by observational or theoretical uncertainties. For a wide range of reasonable parameters, our conclusions appear to be robust with respect to such biases.

We present combined VLA and Green Bank Telescope images of \ammonia\ inversion transitions (1,1) and (2,2) toward OMC2 and OMC3. We focus on the relatively quiescent Orion cores, which are away from the Trapezium cluster and have no sign of massive protostars nor evolved star formation, such as IRAS source, water maser, and methanol maser. The 5\arcsec\ angular resolution and $0.6 \rm{}km\,s^{-1}$ velocity resolution of these data enable us to study the thermal and dynamic state of these cores at $\sim{}0.02 \rm{}pc$ scales, comparable to or smaller than those of the current dust continuum surveys. We measure temperatures for a total of 30 cores, with average masses of $11\,\Ms$, radii of $0.039 \rm{}pc$, virial mass ratio $\bar{R_{vir}}$ = 3.9, and critical mass ratio $\bar{R_{C}}$ = 1.5. Twelve sources contain \textit{Spitzer} protostars. The thus defined starless and protostellar subsamples have similar temperature, line width, but different masses, with an average of $7.3\,\Ms$ for the former and $16\,\Ms$ for the latter. Compared to others Gould Belt dense cores, mores Orion cores have a high gravitational-to kinetic energy ratio and more cores have a larger thant unity critical mass ratio. Orion dense cores have velocity dispersion similar to those of cores in low-mass star-forming regions but larger masses for fiven size. Some cores appear to have truly supercritical gravitational-to-kinetic energy ratios, even when considering significant observational uncertainties: thermal and non-thermal gas mothins alone cannot prevent collapse.

(Abridged) Context. Core condensation is a critical step in the star-formation process, but is still poorly characterized observationally. Aims. We have studied the 10 pc-long L1495/B213 complex in Taurus to investigate how dense cores have condensed out of the lower-density cloud material. Results. From the N$_2$H$^+$ emission, we identify 19 dense cores, some starless and some protostellar. They are not distributed uniformly, but tend to cluster with relative separations on the order of 0.25 pc. From the C$^{18}$O emission, we identify multiple velocity components in the gas. We have characterized them by fitting gaussians to the spectra, and by studying the distribution of the fits in position-position-velocity space. In this space, the C$^{18}$O components appear as velocity-coherent structures, and we have identified them automatically using a dedicated algorithm (FIVe: Friends In Velocity). Using this algorithm, we have identified 35 filamentary components with typical lengths of 0.5 pc, sonic internal velocity dispersions, and mass-per-unit-length close to the stability threshold of isothermal cylinders at 10 K. Core formation seems to have occurred inside the filamentary components via fragmentation, with a small number of fertile components with larger mass-per-unit-length being responsible for most cores in the cloud. At large scales, the filamentary components appear grouped into families, which we refer to as bundles. Conclusions. Core formation in L1495/B213 has proceeded by hierarchical fragmentation. The cloud fragmented first into several pc-scale regions. Each of these regions later fragmented into velocity-coherent filaments of about 0.5 pc in length. Finally, a small number of these filaments fragmented quasi-statically and produced the individual dense cores we see today.

We present the first interferometric molecular line and dust emission maps for the Galactic Center (GC) cloud G0.253+0.016, observed using the Combined Array for Research in Millimeter--wave Astronomy (CARMA) and the Submillimeter Array (SMA). This cloud is very dense, and concentrates a mass exceeding the Orion Molecular Cloud Complex (2x10^5 M_sun) into a radius of only 3pc, but it is essentially starless. G0.253+0.016 therefore violates "star formation laws" presently used to explain trends in galactic and extragalactic star formation by a factor ~45. Our observations show a lack of dense cores of significant mass and density, thus explaining the low star formation activity. Instead, cores with low densities and line widths 1km/s---probably the narrowest lines reported for the GC region to date---are found. Evolution over several 10^5 yr is needed before more massive cores, and possibly an Arches--like stellar cluster, could form. Given the disruptive dynamics of the GC region, and the potentially unbound nature of G0.253+0.016, it is not clear that this evolution will happen.

We present a hierarchical Bayesian method for fitting infrared spectral energy distributions (SEDs) of dust emission to observed fluxes. Under the standard assumption of optically thin single temperature (T) sources the dust SED as represented by a power--law modified black body is subject to a strong degeneracy between T and the spectral index beta. The traditional non-hierarchical approaches, typically based on chi-square minimization, are severely limited by this degeneracy, as it produces an artificial anti-correlation between T and beta even with modest levels of observational noise. The hierarchical Bayesian method rigorously and self-consistently treats measurement uncertainties, including calibration and noise, resulting in more precise SED fits. As a result, the Bayesian fits do not produce any spurious anti-correlations between the SED parameters due to measurement uncertainty. We demonstrate that the Bayesian method is substantially more accurate than the chi-square fit in recovering the SED parameters, as well as the correlations between them. As an illustration, we apply our method to Herschel and sub millimeter ground-based observations of the star-forming Bok globule CB244. This source is a small, nearby molecular cloud containing a single low-mass protostar and a starless core. We find that T and beta are weakly positively correlated -- in contradiction with the chi-square fits, which indicate a T-beta anti-correlation from the same data-set. Additionally, in comparison to the chi-square fits the Bayesian SED parameter estimates exhibit a reduced range in values.

Despite over 30 years of study, the mass-area relationship within and among clouds is still poorly understood both observationally and theoretically. Modern extinction datasets should have sufficient resolution and dynamic range to characterize this relationship for nearby molecular clouds, although recent papers using extinction data seem to yield different interpretations regarding the nature and universality of this aspect of cloud structure. In this paper we try to unify these various results and interpretations by accounting for the different ways cloud properties are measured and analyzed. We interpret the mass-area relationship in terms of the column density distribution function and its possible variation within and among clouds. We quantitatively characterize regional variations in the column density PDF. We show that structures both within and among clouds possess the same degree of "universality", in that their PDF means do not systematically scale with structure size. Because of this, mass scales linearly with area.

H2D+ is a primary ion which dominates the gas-phase chemistry of cold dense gas. Therefore it is hailed as a unique tool in probing the earliest, prestellar phase of star formation. Observationally, its abundance and distribution is however just beginning to be understood in low-mass prestellar and cluster-forming cores. In high mass star forming regions, H2D+ has been detected only in two cores, and its spatial distribution remains unknown. Here we present the first map of the 372 GHz ortho-H2D+ and N2H+ 4-3 transition in the DR21 filament of Cygnus-X with the JCMT, and N2D+ 3--2 and dust continuum with the SMA. We have discovered five very extended (<= 34000 AU diameter) weak structures in H2D+ in the vicinity of, but distinctly offset from embedded protostars. More surprisingly, the H2D+ peak is not associated with either a dust continuum or N2D+ peak. We have therefore uncovered extended massive cold dense gas that was undetected with previous molecular line and dust continuum surveys of the region. This work also shows that our picture of the structure of cores is too simplistic for cluster forming cores and needs to be refined: neither dust continuum with existing capabilities, nor emission in tracers like N2D+ can provide a complete census of the total prestellar gas in such regions. Sensitive H2D+ mapping of the entire DR21 filament is likely to discover more of such cold quiescent gas reservoirs in an otherwise active high mass star-forming region.

We present the detection of a dust continuum source at 3-mm (CARMA) and 1.3-mm (SMA), and 12CO(2-1) emission (SMA) towards the L1451-mm dense core. These detections suggest a compact object and an outflow where no point source at mid-infrared wavelengths is detected using Spitzer. An upper limit for the dense core bolometric luminosity of 0.05 Lsun is obtained. By modeling the broadband SED and the continuum interferometric visibilities simultaneously, we confirm that a central source of heating is needed to explain the observations. This modeling also shows that the data can be well fitted by a dense core with a YSO and disk, or by a dense core with a central First Hydrostatic Core (FHSC). Unfortunately, we are not able to decide between these two models, which produce similar fits. We also detect 12CO(2-1) emission with red- and blue-shifted emission suggesting the presence of a slow and poorly collimated outflow, in opposition to what is usually found towards young stellar objects but in agreement with prediction from simulations of a FHSC. This presents the best candidate, so far, for a FHSC, an object that has been identified in simulations of collapsing dense cores. Whatever the true nature of the central object in L1451-mm, this core presents an excellent laboratory to study the earliest phases of low-mass star formation.

We report the detection of a compact (of order 5 arcsec; about 1800 AU projected size) CO outflow from L1148-IRS. This confirms that this Spitzer source is physically associated with the nearby (about 325 pc) L1148 dense core. Radiative transfer modeling suggests an internal luminosity of 0.08 to 0.13 L_sun. This validates L1148-IRS as a Very Low Luminosity Object (VeLLO; L < 0.1 L_sun). The L1148 dense core has unusually low densities and column densities for a star-forming core. It is difficult to understand how L1148-IRS might have formed under these conditions. Independent of the exact final mass of this VeLLO (which is likely < 0.24 M_sun), L1148-IRS and similar VeLLOs might hold some clues about the isolated formation of brown dwarfs.

We aim at understanding the massive star formation (MSF) limit $m(r) = 870 M_{\odot} (r/pc)^{1.33}$ in the mass-size space of molecular structures recently proposed by Kauffmann & Pillai (2010). As a first step, we build on the hypothesis of a volume density threshold for overall star formation and the model of Parmentier (2011) to establish the mass-radius relations of molecular clumps containing given masses of star-forming gas. Specifically, we relate the mass $m_{clump}$, radius $r_{clump}$ and density profile slope $-p$ of molecular clumps which contain a mass $m_{th}$ of gas denser than a volume density threshold $\rho_{th}$. In a second step, we use the relation between the mass of embedded-clusters and the mass of their most-massive star to estimate the minimum mass of star-forming gas needed to form a $10\,M_{\odot}$ star. Assuming a star formation efficiency of $SFE \simeq 0.30$, this gives $m_{th,crit} \simeq 150 M_{\odot}$. In a third step, we demonstrate that, for sensible choices of the clump density index ($p \simeq 1.7$) and of the cluster formation density threshold ($n_{th} \simeq 10^4\,cm^{-3}$), the line of constant $m_{th,crit} \simeq 150 M_{\odot}$ in the mass-radius space of molecular structures equates with the MSF limit for spatial scales larger than 0.3\,pc. Hence, the observationally inferred MSF limit of Kauffmann & Pillai is consistent with a threshold in star-forming gas mass beyond which the star-forming gas reservoir is large enough to allow the formation of massive stars. For radii smaller than 0.3\,pc, the MSF limit is shown to be consistent with the formation of a $10\,M_{\odot}$ star out of its individual pre-stellar core of density threshold $n_{th} \simeq 10^5\,cm^{-3}$. The inferred density thresholds for the formation of star clusters and individual stars within star clusters match those previously suggested in the literature.

We present a new high-resolution study of pre-protocluster regions in tracers exclusively probing the coldest and dense gas (NH_2D). The data are used to constrain the chemical, thermal, kinematic, and physical conditions (i.e., densities) in G29.96e and G35.20w. NH_3, NH_2D, and continuum emission were mapped using the VLA, and PdBI. In particular, NH_2D is a unique tracer of cold, precluster gas at high densities, while NH_3 traces both the cold and warm gas of modest-to-high densities. In G29.96e, Spitzer images reveal two massive filaments, one of them in extinction (infrared dark cloud). We observe very low line widths in NH_3 (FWHM <1km/s). These multi-wavelength, high-resolution observations of high-mass pre-protocluster regions show that the target regions are characterized by (i) turbulent Jeans fragmentation of massive clumps into cores (from a Jeans analysis); (ii) cores and clumps that are "over-bound/subvirial", i.e. turbulence is too weak to support them against collapse, meaning that (iii) some models of monolithic cloud collapse are quantitatively inconsistent with data; (iv) accretion from the core onto a massive star, which can (for observed core sizes and velocities) be sustained by accretion of envelope material onto the core, suggesting that (similar to competitive accretion scenarios) the mass reservoir for star formation is not necessarily limited to the natal core; (v) high deuteration ratios ([NH_2D/NH_3]>6%), which make the above discoveries possible; (vi) and the destruction of NH_2D toward embedded stars. [abridged]

We present a new assessment of the ability of Infrared Dark Clouds (IRDCs) to form massive stars and clusters. This is done by comparison with an empirical mass-size threshold for massive star formation (MSF). We establish m(r)>870M_sun(r/pc)^1.33 as a novel approximate MSF limit, based on clouds with and without MSF. Many IRDCs, if not most, fall short of this threshold. Without significant evolution, such clouds are unlikely MSF candidates. This provides a first quantitative assessment of the small number of IRDCs evolving towards MSF. IRDCs below this limit might still form stars and clusters of up to intermediate mass, though (like, e.g., the Ophiuchus and Perseus Molecular Clouds). Nevertheless, a major fraction of the mass contained in IRDCs might reside in few 10^2 clouds sustaining MSF.

We measure the mass and size of cloud fragments in several molecular clouds continuously over a wide range of spatial scales (0.05 < r / pc < 3). Based on the recently developed "dendrogram-technique", this characterizes dense cores as well as the enveloping clouds. "Larson's 3rd Law" of constant column density, m(r) = C*r^2, is not well suited to describe the derived mass-size data. Solar neighborhood clouds not forming massive stars (< 10 M_sun; Pipe Nebula, Taurus, Perseus, and Ophiuchus) obey m(r) < 870 M_sun (r / pc)^1.33 . In contrast to this, clouds forming massive stars (Orion A, G10.15$-$0.34, G11.11$-$0.12) do exceed the aforementioned relation. Thus, this limiting mass-size relation may approximate a threshold for the formation of massive stars. Across all clouds, cluster-forming cloud fragments are found to be---at given radius---more massive than fragments devoid of clusters. The cluster-bearing fragments are found to roughly obey a mass-size law m = C*r^1.27 (where the exponent is highly uncertain in any given cloud, but is certainly smaller than 1.5).

Power law mass-size and linewidth-size correlations, two of "Larson's laws," are often studied to assess the dynamical state of clumps within molecular clouds. Using the result of a hydrodynamic simulation of a molecular cloud, we investigate how geometric projection may affect the derived Larson relationships. We find that large scale structures in the column density map have similar masses and sizes to those in the 3D simulation (PPP). Smaller scale clumps in the column density map are measured to be more massive than the PPP clumps, due to the projection of all emitting gas along lines of sight. Further, due to projection effects, structures in a synthetic spectral observation (PPV) may not necessarily correlate with physical structures in the simulation. In considering the turbulent velocities only, the linewidth-size relationship in the PPV cube is appreciably different from that measured from the simulation. Including thermal pressure in the simulated linewidths imposes a minimum linewidth, which results in a better agreement in the slopes of the linewidth-size relationships, though there are still discrepancies in the offsets, as well as considerable scatter. Employing commonly used assumptions in a virial analysis, we find similarities in the computed virial parameters of the structures in the PPV and PPP cubes. However, due to the discrepancies in the linewidth- and mass- size relationships in the PPP and PPV cubes, we caution that applying a virial analysis to observed clouds may be misleading due to geometric projection effects. We speculate that consideration of physical processes beyond kinetic and gravitational pressure would be required for accurately assessing whether complex clouds, such as those with highly filamentary structure, are bound.

We use a new contour-based map analysis technique to measure the mass and size of molecular cloud fragments continuously over a wide range of spatial scales (0.05 < r / pc < 10), i.e., from the scale of dense cores to those of entire clouds. The present paper presents the method via a detailed exploration of the Perseus Molecular Cloud. Dust extinction and emission data are combined to yield reliable scale-dependent measurements of mass. This scale-independent analysis approach is useful for several reasons. First, it provides a more comprehensive characterization of a map (i.e., not biased towards a particular spatial scale). Such a lack of bias is extremely useful for the joint analysis of many data sets taken with different spatial resolution. This includes comparisons between different cloud complexes. Second, the multi-scale mass-size data constitutes a unique resource to derive slopes of mass-size laws (via power-law fits). Such slopes provide singular constraints on large-scale density gradients in clouds.

In this paper we present a dust emission map of the starless core TMC-1C taken at 2100 microns. Along with maps at 160, 450, 850 and 1200 microns, we study the dust emissivity spectral index from the (sub)millimeter spectral energy distribution, and find that it is close to the typically assumed value of beta = 2. We also map the dust temperature and column density in TMC-1C, and find that at the position of the dust peak (A_V ~ 50), the line-of-sight-averaged temperature is ~7 K. Employing simple Monte Carlo modeling, we show that the data are consistent with a constant value for the emissivity spectral index over the whole map of TMC-1C.

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 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 micron 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.

We investigate how uncertainties in flux measurements affect the results from modified blackbody SED fits. We show that an inverse correlation between the dust temperature T and spectral index (beta) naturally arises from least squares fits due to the uncertainties, even for sources with a single T and beta. Fitting SEDs to noisy fluxes solely in the Rayleigh-Jeans regime produces unreliable T and beta estimates. Thus, for long wavelength observations (lambda >~ 200 micron), or for warm sources (T >~ 60 K), it becomes difficult to distinguish sources with different temperatures. We assess the role of noise in recent observational results that indicate an inverse and continuously varying T - beta relation. Though an inverse and continuous T - beta correlation may be a physical property of dust in the ISM, we find that the observed inverse correlation may be primarily due to noise.

We present a search for all embedded protostars with internal luminosities < 1 solar luminosity in the sample of nearby, low-mass star-forming regions surveyed by the Spitzer Space Telescope c2d Legacy Project. The internal luminosity (Lint) of a source is the luminosity of the central source and excludes luminosity arising from external heating. On average, the c2d data are sensitive to embedded protostars with Lint > 4E-3 (d/140 pc)^2 solar luminosities, a factor of 25 better than the sensitivity of IRAS to such objects. We present selection criteria used to identify candidates from the Spitzer data and examine complementary data to decide whether each candidate is truly an embedded protostar. We find a tight correlation between the 70 micron flux and internal luminosity of a protostar, an empirical result based on observations and two-dimensional radiative transfer models of protostars. We identify 50 embedded protostars with Lint < 1 solar luminosities; 15 have Lint < 0.1 solar luminosities. The intrinsic distribution of source luminosities increases to lower luminosities. While we find sources down to the above sensitivity limit, indicating that the distribution may extend to luminosities lower than probed by these observations, we are able to rule out a continued rise in the distribution below 0.1 solar luminosities. Between 75-85% of cores classified as starless prior to being observed by Spitzer remain starless to our luminosity sensitivity; the remaining 15-25% harbor low-luminosity, embedded protostars. We compile complete Spectral Energy Distributions for all 50 objects and calculate standard evolutionary signatures, and argue that these objects are inconsistent with the simplest picture of star formation wherein mass accretes from the core onto the protostar at a constant rate.

AIMS: To study the structure of nearby (< 500 pc) dense starless and star-forming cores with the particular goal to identify and understand evolutionary trends in core properties, and to explore the nature of Very Low Luminosity Objects (< 0.1 L_sun; VeLLOs). METHODS: Using the MAMBO bolometer array, we create maps unusually sensitive to faint (few mJy per beam) extended (approx. 5 arcmin) thermal dust continuum emission at 1.2 mm wavelength. Complementary information on embedded stars is obtained from Spitzer, IRAS, and 2MASS. RESULTS: Our maps are very rich in structure, and we characterize extended emission features (``subcores'') and compact intensity peaks in our data separately to pay attention to this complexity. We derive, e.g., sizes, masses, and aspect ratios for the subcores, as well as column densities and related properties for the peaks. Combination with archival infrared data then enables the derivation of bolometric luminosities and temperatures, as well as envelope masses, for the young embedded stars. CONCLUSIONS: (abridged) Starless and star-forming cores occupy the same parameter space in many core properties; a picture of dense core evolution in which any dense core begins to actively form stars once it exceeds some fixed limit in, e.g., mass, density, or both, is inconsistent with our data. Comparison of various evolutionary indicators for young stellar objects in our sample (e.g., bolometric temperatures) reveals inconsistencies between some of them, possibly suggesting a revision of some of these indicators.

We demonstrate the utility of dendrograms at representing the essential features of the hierarchical structure of the isosurfaces for molecular line data cubes. The dendrogram of a data cube is an abstraction of the changing topology of the isosurfaces as a function of contour level. The ability to track hierarchical structure over a range of scales makes this analysis philosophically different from local segmentation algorithms like CLUMPFIND. Points in the dendrogram structure correspond to specific volumes in data cubes defined by their bounding isosurfaces. We further refine the technique by measuring the properties associated with each isosurface in the analysis allowing for a multiscale calculation of molecular gas properties. Using COMPLETE 13CO(1-0) data from the L1448 region in Perseus and mock observations of a simulated data cube, we identify regions that have a significant contribution by self-gravity to their energetics on a range of scales. We find evidence for self-gravitation on all spatial scales in L1448 though not in all regions. In the simulated observations, nearly all of the emission is found in objects that would be self-gravitating if gravity were included in the simulation. We reconstruct the size-line width relationship within the data cube using the dendrogram-derived properties and find it follows the standard relation: s_v ~ R^0.58. Finally, we show that constructing the dendrogram of CO J=1-0 emission from the Orion-Monoceros region allows for the identification of giant molecular clouds in a blended molecular line data set using only a physically motivated definition (self-gravitating clouds with masses 5x10^4 Msun.

We present ammonia observations of 193 dense cores and core candidates in the Perseus molecular cloud made using the Robert F. Byrd Green Bank Telescope. We simultaneously observed the NH3(1,1), NH3(2,2), CCS (2_1 -> 1_0) and CC34S (2_1 -> 1_0) transitions near 23 GHz for each of the targets with a spectral resolution of dv ~ 0.024 km/s. We find ammonia emission associated with nearly all of the (sub)millimeter sources as well as at several positions with no associated continuum emission. For each detection, we have measured physical properties by fitting a simple model to every spectral line simultaneously. Where appropriate, we have refined the model by accounting for low optical depths, multiple components along the line of sight and imperfect coupling to the GBT beam. For the cores in Perseus, we find a typical kinetic temperature of T=11 K, a typical column density of N(NH3)~ 10^14.5 /cm^2 and velocity dispersions ranging from sigma_v = 0.07 km/s to 0.7 km/s. However, many cores with velocity dispersions > 0.2 km/s show evidence for multiple velocity components along the line of sight.

We have mapped the central 10'x10' of the dense core TMC-1C at 450, 850 and 1200 microns using SCUBA on the James Clerk Maxwell Telescope and MAMBO on the IRAM 30m telescope. We show that although one can, in principle, use images at these wavelengths to map the emissivity spectral index, temperature and column density independently, noise and calibration errors would have to be less than 2% to accurately derive these three quantities from a set of three emission maps. Because our data are not this free of errors, we use our emission maps to fit the dust temperature and column density assuming a constant value of the emissivity spectral index and explore the effects of noise on the derived physical parameters. We find that the derived extinction values for TMC-1C are large for a starless core (80 mag Av), and the derived temperatures are low (6 K) in the densest regions of the core, using our derived value of beta = 1.8.

We present Spitzer Space Telescope observations of the "evolved starless core" L1521F which reveal the presence of a very low luminosity object (L < 0.07 Lsun). The object, L1521F-IRS, is directly detected at mid-infrared wavelengths (>5 micron) but only in scattered light at shorter infrared wavelengths, showing a bipolar nebula oriented east-west which is probably tracing an outflow cavity. The nebula strongly suggests that L1521F-IRS is embedded in the L1521F core. Thus L1521F-IRS is similar to the recently discovered L1014-IRS and the previously known IRAM 04191 in its substellar luminosity and dense core environment. However these objects differ significantly in their core density, core chemistry, and outflow properties, and some may be destined to be brown dwarfs rather than stars.

The Spitzer Space Telescope allows for the first time to search systematically for very low luminosity (= 0.1 L_sun) objects (VeLLOs) associated with dense molecular cores. They may be the first candidate Class 0 sources with sub-stellar masses. We describe such a source in the dense molecular core L1148. VeLLO natal cores show properties that are unusual for star-forming cores. The low luminosity and in some cases the lack of prominent outflow could be the result of the small gas supply near the VeLLO.

We present maps of over 1.5 square degrees in Chamaeleon (Cha) II at 24, 70, and 160 micron observed with the Spitzer Space Telescope Multiband Imaging Photometer for Spitzer (MIPS) and a 1.2 square degree millimeter map from SIMBA on the Swedish-ESO Submillimetre Telescope (SEST). The c2d Spitzer Legacy Team's data reduction pipeline is described in detail. Over 1500 24 micron sources and 41 70 micron sources were detected by MIPS with fluxes greater than 10-sigma. More than 40 potential YSOs are identified with a MIPS and 2MASS color-color diagram and by their spectral indices, including two previously unknown sources with 24 micron excesses. Our new SIMBA millimeter map of Cha II shows that only a small fraction of the gas is in compact structures with high column densities. The extended emission seen by MIPS is compared with previous CO observations. Some selected interesting sources, including two detected at 1 mm, associated with Cha II are discussed in detail and their SEDs presented. The classification of these sources using MIPS data is found to be consistent with previous studies.

Spitzer Space Telescope observations of a point-like source, L1014-IRS, close to the dust peak of the low-mass dense core L1014 have questioned its starless nature. The presence of an object with colors expected for an embedded protostar makes L1014-IRS the lowest luminosity isolated protostar known, and an ideal target with which to test star formation theories at the low mass end. In order to study its molecular content and to search for the presence of a molecular outflow, we mapped L1014 in at least one transition of 12CO, N2H+, HCO+, CS and of their isotopologues 13CO, C18O, C17O, N2D+ and H13CO+, using the FCRAO, the IRAM 30 meter and the CSO. The data show physical and chemical properties in L1014 typical of the less evolved starless cores: i.e. H2 central density of a few 10^5 molecules cm^-3, estimated mass of ~2M_sun, CO integrated depletion factor less than 10, N(N2H+)~6*10^12 cm^-2, N(N2D+)/N(N2H+) equal to 10% and relatively broad N2H+(1--0) lines (0.35 km/s). Infall signatures and significant velocity shifts between optically thick and optically thin tracers are not observed in the line profiles. No classical signatures of molecular outflow are found in the 12CO and 13CO observations. In particular, no high velocity wings are found, and no well-defined blue-red lobes of 12CO emission are seen in the channel maps. If sensitive, higher resolution observations confirm the absence of an outflow on a smaller scale than probed by our observations, L1014-IRS would be the only protostellar object known to be formed without driving an outflow.

We present observations of L1014, a dense core in the Cygnus region previously thought to be starless, but data from the Spitzer Space Telescope shows the presence of an embedded source. We propose a model for this source that includes a cold core, heated by the interstellar radiation field, and a low-luminosity internal source. The low luminosity of the internal source suggests a substellar object. If L1014 is representative, other "starless" cores may turn out to harbor central sources.

Following an approach initially outlined by McKee & Holliman, we investigate the structure and stability of dense, starless molecular cloud cores. We model those as spherical clouds in hydrostatic equilibrium and supported against gravity by thermal, turbulent, and magnetic pressure. We determine the gas pressure by solving for thermal equilibrium between heating and cooling, while the turbulent and magnetic pressures are assumed to obey polytropic equations of state. In comparing the models to observed cloud cores we find that the observed peak column densities often exceed the limit for stable equilibria supported by thermal pressure alone, suggesting significant non-thermal pressure if the cores are to be stable. Non-thermal support is also needed to stabilize cores embedded in molecular clouds with high average pressures. Since the observed molecular linewidths of cores suggest that the turbulent pressure is lower than the thermal pressure, magnetic field are likely a dominant pressure component in many such cores.