The majority of stars are part of gravitationally bound stellar systems, such as binaries. Observations of protobinary systems constrain the conditions that lead to stellar multiplicity and subsequent orbital evolution. We report high-angular resolution observations of the circumbinary disk around [BHB2007] 11, a young binary protostar system. The two protostars are embedded in circumstellar disks that have radii of 2 to 3 astronomical units and probably contain a few Jupiter masses. These systems are surrounded by a complex structure of filaments connecting to the larger circumbinary disk. We also observe accretion and radio jets associated with the protobinary system. The accretion is preferentially onto the lower-mass protostar, consistent with theoretical predictions.

The giant molecular cloud Sagittarius B2 (hereafter SgrB2) is the most massive region with ongoing high-mass star formation in the Galaxy. In the southern region of the 40-pc large envelope of SgrB2, we encounter the SgrB2(DS) region which hosts more than 60 high-mass protostellar cores distributed in an arc shape around an extended HII region. We use the Very Large Array in its CnB and D configurations, and in the frequency bands C (4--8 GHz) and X (8--12 GHz) to observe the whole SgrB2 complex. Continuum and radio recombination line maps are obtained. We detect radio continuum emission in SgrB2(DS) in a bubble-shaped structure. From 4 to 12 GHz, we derive a spectral index between -1.2 and -0.4, indicating the presence of non-thermal emission. We decompose the contribution from thermal and non-thermal emission, and find that the thermal component is clumpy and more concentrated, while the non-thermal component is more extended and diffuse. The radio recombination lines in the region are found to be not in local thermodynamic equilibrium (LTE) but stimulated by the non-thermal emission. The thermal free-free emission is likely tracing an HII region ionized by an O7 star, while the non-thermal emission can be generated by relativistic electrons created through first-order Fermi acceleration. We have developed a simple model of the SgrB2(DS) region and found that first-order Fermi acceleration can reproduce the observed flux density and spectral index.

We have used an unbiased, spectral line-survey that covers the frequency range from 211 to 275 GHz and was obtained with ALMA (angular resolution of 0.4 arcsec) to study the small-scale structure of the dense gas in Sagittarius B2 (north). Eight filaments are found converging to the central hub and extending for about 0.1 pc. The spatial structure, together with the presence of the massive central region, suggest that these filaments may be associated with accretion processes. In order to derive the kinematic properties of the gas in a chemically line-rich source like Sgr B2(N), we have developed a new tool that stacks all the detected transition lines of any molecular species. This permits to increase the signal-to-noise ratio of our observations and average out line blending effects, which are a common problem in line-rich regions. We derive velocity gradients along the filaments of about 20-100 km s$^{-1}$ pc$^{-1}$, which are 10-100 times larger than those typically found on larger scales (1 pc) in other star-forming regions. The mass accretion rates of individual filaments are about 0.05 M$_\odot$ yr$^{-1}$, which result in a total accretion rate of 0.16 M$_\odot$ yr$^{-1}$. Some filaments harbor dense cores that are likely forming stars and stellar clusters. The stellar content of these dense cores is on the order of 50% of the total mass. We conclude that the cores may merge in the center when already forming stellar clusters but still containing a significant amount of gas, resulting in a "damp" merger. The high density and mass of the central region, combined with the presence of converging filaments with high mass, high accretion rates and embedded dense cores already forming stars, suggest that Sgr B2(N) may have the potential to evolve into a super stellar cluster.

The determination of the specific angular momentum radial profile, $j(r)$, in the early stages of star formation is crucial to constrain star and circumstellar disk formation theories. The specific angular momentum is directly related to the largest Keplerian disk possible, and it could constrain the angular momentum removal mechanism. We determine $j(r)$ towards two Class 0 objects and a first hydrostatic core candidate in the Perseus cloud, which is consistent across all three sources and well fit with a single power-law relation between 800 and 10,000\,au: $j_{fit}(r)=10^{-3.60\pm0.15}\left(r/\textrm{1,000au}\right)^{1.80\pm0.04}$ km s$^{-1}$ pc. This power-law relation is in between solid body rotation ($\propto r^2$) and pure turbulence ($\propto r^{1.5}$). This strongly suggests that even at 1,000\,au, the influence of the dense core's initial level of turbulence or the connection between core and the molecular cloud is still present. The specific angular momentum at 10,000\,au is $\approx3\times$ higher than previously estimated, while at 1,000\,au it is lower by $2\times$. We do not find a region of conserved specific angular momentum, although it could still be present at a smaller radius. We estimate an upper limit to the largest Keplerian disk radius of 60\,au, which is small but consistent with published upper limits. Finally, these results suggest that more realistic initial conditions for numerical simulations of disk formation are needed. Some possible solutions include: a) use a larger simulation box to include some level of driven turbulence or connection to the parental cloud, or b) incorporate the observed $j(r)$ to setup the dense core kinematics initial conditions.

The evolution of molecular clouds in galactic centres is thought to differ from that in galactic discs due to a significant influence of the external gravitational potential. We present a set of numerical simulations of molecular clouds orbiting on the 100-pc stream of the Central Molecular Zone (the central $\sim500$ pc of the Galaxy) and characterise their morphological and kinematic evolution in response to the background potential and eccentric orbital motion. We find that the clouds are shaped by strong shear and torques, by tidal and geometric deformation, and by their passage through the orbital pericentre. Within our simulations, these mechanisms control cloud sizes, aspect ratios, position angles, filamentary structure, column densities, velocity dispersions, line-of-sight velocity gradients, spin angular momenta, and kinematic complexity. By comparing these predictions to observations of clouds on the Galactic Centre 'dust ridge', we find that our simulations naturally reproduce a broad range of key observed morphological and kinematic features, which can be explained in terms of well-understood physical mechanisms. We argue that the accretion of gas clouds onto the central regions of galaxies, where the rotation curve turns over and the tidal field is fully compressive, is accompanied by transformative dynamical changes to the clouds, leading to collapse and star formation. This can generate an evolutionary progression of cloud collapse with a common starting point, which either marks the time of accretion onto the tidally-compressive region or of the most recent pericentre passage. Together, these processes may naturally produce the synchronised starbursts observed in numerous (extra)galactic nuclei.

We report the discovery of the first 14NH3 (2,2) maser, seen in the Sgr B2 Main star forming region near the center of the Milky Way, using data from the Very Large Array radio telescope. The maser is seen in both lower resolution (3" or ~0.1 pc) data from 2012 and higher resolution (0''.1 or ~1000 AU) data from 2018. In the higher resolution data ammonia (2,2) maser emission is detected toward 5 independent spots. The maser spots are not spatially or kinematically coincident with any other masers in this region, or with the peaks of the radio continuum emission from the numerous ultracompact and hypercompact \hii\, regions in this area. While the (2,2) maser spots are spatially unresolved in our highest resolution observations, they have unusually broad linewidths of several kilometers per second, which suggests that each of these spots consists of multiple masers tracing unresolved velocity structure. No other ammonia lines observed in Sgr B2 Main are seen to be masers, which continues to challenge theories of ammonia, maser emission that predict simultaneous maser emission in multiple ammonia transitions.

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

We aim to investigate the chemistry of internal photon-dominated regions surrounding deeply embedded hypercompact and ultracompact HII regions. We search for specific tracers of this evolutionary stage of massive star formation that can be detected with current astronomical facilities. We modeled hot cores with embedded HC/UCHII regions, by coupling the astrochemical code Saptarsy to a radiative transfer framework obtaining the spatio-temporal evolution of abundances as well as time-dependent synthetic spectra. In these models where we focused on the internal PDR surrounding the HI region, the gas temperature is set to the dust temperature and we do not include dynamics thus the density structure is fixed. We compared this to hot molecular core models and studied the effect on the chemistry of the radiation field which is included in the HII region models only during the computation of abundances. In addition, we investigated the chemical evolution of the gas surrounding HII regions with models of different densities at the ionization front, different sizes of the ionized cavity and different initial abundances. We obtain the time evolution of synthetic spectra for a dozen of selected species as well as ratios of their integrated intensities. We find that some molecules such as C, N2H+, CN, and HCO do not trace the inner core and so are not good tracers to distinguish the HII/PDR regions to the HMCs phase. On the contrary, C+ and O trace the internal PDRs, in the two models starting with different initial abundances, but are unfortunately currently unobservable with the current achievable spatial resolution because of the very thin internal PDR (r < 100 AU). In addition, we find that the abundance profiles are highly affected by the choice of the initial abundances, hence the importance to properly define them.

We model the emission of methyl cyanide (CH3CN) lines towards the massive hot molecular core SgrB2(M). We aim at reconstructing the CH3CN abundance field and investigating the gas temperature distribution as well as the velocity field. SgrB2(M) was observed with the ALMA in a spectral line survey from 211 to 275 GHz. This frequency range includes several transitions of CH3CN (including isotopologues and vibrationally excited states). We employ the three-dimensional radiative transfer toolbox Pandora in order to retrieve the velocity and abundance field by modeling different CH3CN lines. For this purpose, we base our model on the results of a previous study that determined the physical structure of SgrB2(M), i.e.\ the distribution of dust dense cores, ionized regions and heating sources. The morphology of the CH3CN emission can be reproduced by a molecular density field that consists of a superposition of cores with modified Plummer-like density profiles. The averaged relative abundance of CH3CN with respect to H2 ranges from 4x10^{-11} to 2x10^{-8} in the northern part of SgrB2(M) and from 2x10^{-10} to 5x10^{-7} in the southern part. In general, we find that the relative abundance of CH3CN is lower at the center of the very dense and hot cores, causing the general morphology of the CH3CN emission to be shifted with respect to the dust continuum emission. The dust temperature calculated by the radiative transfer simulation based on the available luminosity reaches values up to 900 K. However, in some regions vibrationally excited transitions of CH3CN are underestimated by the model, indicating that the predicted gas temperature, which is assumed to be equal to the dust temperature, is partly underestimated. The determination of the velocity component along the line of sight reveals that a velocity gradient from the north to the south exists in SgrB2(M).

We report ALMA observations with resolution $\approx0.5$" at 3 mm of the extended Sgr B2 cloud in the Central Molecular Zone (CMZ). We detect 271 compact sources, most of which are smaller than 5000 AU. By ruling out alternative possibilities, we conclude that these sources consist of a mix of hypercompact HII regions and young stellar objects (YSOs). Most of the newly-detected sources are YSOs with gas envelopes which, based on their luminosities, must contain objects with stellar masses $M_*\gtrsim8$ M$_\odot$. Their spatial distribution spread over a $\sim12\times3$ pc region demonstrates that Sgr B2 is experiencing an extended star formation event, not just an isolated `starburst' within the protocluster regions. Using this new sample, we examine star formation thresholds and surface density relations in Sgr B2. While all of the YSOs reside in regions of high column density ($N(H_2)\gtrsim2\times10^{23}$ cm$^{-2}$), not all regions of high column density contain YSOs. The observed column density threshold for star formation is substantially higher than that in solar vicinity clouds, implying either that high-mass star formation requires a higher column density or that any star formation threshold in the CMZ must be higher than in nearby clouds. The relation between the surface density of gas and stars is incompatible with extrapolations from local clouds, and instead stellar densities in Sgr B2 follow a linear $\Sigma_*-\Sigma_{gas}$ relation, shallower than that observed in local clouds. Together, these points suggest that a higher volume density threshold is required to explain star formation in CMZ clouds.

STATCONT is a python-based tool designed to determine the continuum emission level in spectral data, in particular for sources with a line-rich spectrum. The tool inspects the intensity distribution of a given spectrum and automatically determines the continuum level by using different statistical approaches. The different methods included in STATCONT are tested against synthetic data. We conclude that the sigma-clipping algorithm provides the most accurate continuum level determination, together with information on the uncertainty in its determination. This uncertainty can be used to correct the final continuum emission level, resulting in the here called `corrected sigma-clipping method' or c-SCM. The c-SCM has been tested against more than 750 different synthetic spectra reproducing typical conditions found towards astronomical sources. The continuum level is determined with a discrepancy of less than 1% in 50% of the cases, and less than 5% in 90% of the cases, provided at least 10% of the channels are line free. The main products of STATCONT are the continuum emission level, together with a conservative value of its uncertainty, and datacubes containing only spectral line emission, i.e., continuum-subtracted datacubes. STATCONT also includes the option to estimate the spectral index, when different files covering different frequency ranges are provided.

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.

The high-mass star forming sites SgrB2(M) and SgrB2(N) have been the target of numerous studies, revealing e.g. a rich chemistry. We want to characterize their physical and chemical structure using ALMA high-angular resolution observations at mm wavelengths, reaching spatial scales of about 4000 au, and covering the whole band 6 (from 211 to 275 GHz). In order to determine the continuum emission in line-rich sources, we use a new statistical method: STATCONT. We detect 27 continuum sources in SgrB2(M) and 20 in SgrB2(N). We study the continuum emission across the ALMA band 6, and compare it with previous SMA 345 GHz and VLA 40 GHz observations, to study the nature of the sources detected. The brightest sources are dominated by (partially optically thick) dust emission, while there is an important degree of contamination from ionized gas free-free emission in weaker sources. While the total mass in SgrB2(M) is distributed in many fragments, most of the mass in SgrB2(N) arises from a single object, with filamentary-like structures converging towards the center. There seems to be a lack of low-mass dense cores in both regions. We determine H2 volume densities for the cores of about 10^5-10^7 Msun pc^-3, one to two orders of magnitude higher than the stellar densities of super star clusters. In general, SgrB2(N) is chemically richer than SgrB2(M). There seems to be a correlation between the chemical richness and the mass of the fragments, with more massive clumps being more chemically rich. Both SgrB2(N) and SgrB2(M) harbour a cluster of hot molecular cores. We compare the continuum images with predictions from a detailed 3D radiative transfer model that reproduces the structure of SgrB2 from 45 pc down to 100 au. This dataset, together with ongoing projects in the range 5 to 200 GHz, better constrain the 3D structure of SgrB2, and allow us to understand its physical and chemical structure.

We report the results of a sensitive search for the 443.952902 GHz $J=1-0$ transition of the LiH molecule toward two interstellar clouds in the Milky Way, W49N and Sgr B2 (Main), that has been carried out using the Atacama Pathfinder Experiment (APEX) telescope. The results obtained toward W49N place an upper limit of $1.9 \times 10^{-11}\, (3\sigma)$ on the LiH abundance, $N({\rm LiH})/N({\rm H}_2)$, in a foreground, diffuse molecular cloud along the sight-line to W49N, corresponding to 0.5% of the solar system lithium abundance. Those obtained toward Sgr B2 (Main) place an abundance limit $N({\rm LiH})/N({\rm H}_2) < 3.6 \times 10^{-13} \,(3\sigma)$ in the dense gas within the Sgr B2 cloud itself. These limits are considerably smaller that those implied by the tentative detection of LiH reported previously for the $z=0.685$ absorber toward B0218+357.

Constraining the temperature and density structure of dense molecular cloud cores is fundamental for understanding the initial conditions of star formation. We use Herschel observations of the thermal FIR dust emission from nearby isolated molecular cloud cores and combine them with ground-based submillimeter continuum data to derive observational constraints on their temperature and density structure. The aim of this study is to verify the validity of a ray-tracing inversion technique developed to derive the dust temperature and density structure of isolated starless cores directly from the dust emission maps and to test if the resulting temperature and density profiles are consistent with physical models. Using this ray-tracing inversion technique, we derive the dust temperature and density structure of six isolated starless cloud cores. We employ self-consistent radiative transfer modeling to the derived density profiles, treating the ISRF as the only heating source. The best-fit values of local strength of the ISRF and the extinction by the outer envelope are derived by comparing the self-consistently calculated temperature profiles with those derived by the ray-tracing method. We find that all starless cores are significantly colder inside than outside, with the core temperatures showing a strong negative correlation with peak column density. This suggests that their thermal structure is dominated by external heating from the ISRF and shielding by dusty envelopes. The temperature profiles derived with the ray-tracing inversion method can be well-reproduced with self-consistent radiative transfer models.

Using spectral-line observations of HNCO, N2H+, and HNC, we investigate the kinematics of dense gas in the central ~250 pc of the Galaxy. We present SCOUSE (Semi-automated multi-COmponent Universal Spectral-line fitting Engine), a line fitting algorithm designed to analyse large volumes of spectral-line data efficiently and systematically. Unlike techniques which do not account for complex line profiles, SCOUSE accurately describes the {l, b, v_LSR} distribution of CMZ gas, which is asymmetric about Sgr A* in both position and velocity. Velocity dispersions range from 2.6 km/s<\sigma<53.1 km/s. A median dispersion of 9.8 km/s, translates to a Mach number, M_3D>28. The gas is distributed throughout several "streams", with projected lengths ~100-250 pc. We link the streams to individual clouds and sub-regions, including Sgr C, the 20 and 50 km/s clouds, the dust ridge, and Sgr B2. Shell-like emission features can be explained by the projection of independent molecular clouds in Sgr C and the newly identified conical profile of Sgr B2 in {l ,b, v_LSR} space. These features have previously invoked supernova-driven shells and cloud-cloud collisions as explanations. We instead caution against structure identification in velocity-integrated emission maps. Three geometries describing the 3-D structure of the CMZ are investigated: i) two spiral arms; ii) a closed elliptical orbit; iii) an open stream. While two spiral arms and an open stream qualitatively reproduce the gas distribution, the most recent parameterisation of the closed elliptical orbit does not. Finally, we discuss how proper motion measurements of masers can distinguish between these geometries, and suggest that this effort should be focused on the 20 km/s and 50 km/s clouds and Sgr C.

We model the dust and free-free continuum emission in the high-mass star-forming region Sagittarius B2 in order to reconstruct the three-dimensional density and dust temperature distribution, as a crucial input to follow-up studies of the gas velocity field and molecular abundances. We employ the three-dimensional radiative transfer program RADMC-3D to calculate the dust temperature self-consistently, provided a given initial density distribution. This density distribution of the entire cloud complex is then recursively reconstructed based on available continuum maps, including both single-dish and high-resolution interferometric maps covering a wide frequency range (40 GHz - 4 THz). The model covers spatial scales from 45 pc down to 100 au, i.e. a spatial dynamic range of 10^5. We find that the density distribution of Sagittarius B2 can be reasonably well fitted by applying a superposition of spherical cores with Plummer-like density profiles. In order to reproduce the spectral energy distribution, we position Sgr B2(N) along the line of sight behind the plane containing Sgr B2(M). We find that the entire cloud complex comprises a total gas mass of 8.0 x 10^6 Msun within a diameter of 45 pc, corresponding to an averaged gas density of 170 Msun/pc^3. We estimate stellar masses of 2400 Msun and 20700 Msun and luminosities of 1.8 x 10^6 Lsun and 1.2 x 10^7 Lsun for Sgr B2(N) and Sgr B2(M), respectively. We report H_2 column densities of 2.9 x 10^24 cm^-2 for Sgr B2(N) and 2.5 x 10^24 cm^-2 for Sgr B2(M) in a 40" beam. For Sgr B2(S), we derive a stellar mass of 1100 Msun, a luminosity of 6.6 x 10^5 Lsun and a H_2 column density of 2.2 x 10^24 cm^-2 in a 40" beam. We calculate a star formation efficiency of 5% for Sgr B2(N) and 50% for Sgr B2(M), indicating that most of the gas content in Sgr B2(M) has already been converted to stars or dispersed.

Hot molecular cores (HMCs) are intermediate stages of high-mass star formation and are also known for their rich emission line spectra at (sub-)mm wavebands. The observed spectral feature of HMCs such as total number of emission lines and associated line intensities are also found to vary with evolutionary stages. We developed various 3D models for HMCs guided by the evolutionary scenarios proposed by recent empirical and modeling studies. We then investigated the spatio-temporal variation of temperature and molecular abundances in HMCs by consistently coupling gas-grain chemical evolution with radiative transfer calculations. We explored the effects of varying physical conditions on molecular abundances including density distribution and luminosity evolution of the central protostar(s). The time-dependent temperature structure of the hot core models provides a realistic framework for investigating the spatial variation of ice mantle evaporation as a function of evolutionary timescales. With increasing protostellar luminosity, the water ice evaporation font ($\sim$100K) expands and the spatial distribution of gas phase abundances of these COMs also spreads out. We simulated the synthetic spectra for these models at different evolutionary timescales to compare with observations. A qualitative comparison of the simulated and observed spectra suggests that these self-consistent hot core models can reproduce the notable trends in hot core spectral variation within the typical hot core timescales of 10$^{5}$ year. These models predict that the spatial distribution of various emission line maps will also expand with evolutionary time. The model predictions can be compared with high resolution observation that can probe scales of a few thousand AU in high-mass star forming regions such as from ALMA.[Abridged]

The massive infrared dark cloud G0.253+0.016 projected 45pc from the Galactic centre contains ~10^5Msun of dense gas whilst being mostly devoid of observed star-formation tracers. To scrutinise the physical properties, dynamics and structure of this cloud with reference to its star-forming potential, we have carried out a concerted SMA and IRAM 30m study of this cloud in dust continuum, CO isotopologues, shock tracing molecules, as well as H$_2$CO to trace the gas temperature. We detect and characterise the dust cores within G0.253+0.016 at ~1.3 mm and find that the kinetic temperature of the gas is >320K on size-scales of ~0.15 pc. Analysis of the position-velocity diagrams of our observed lines show broad linewidths and strong shock emission in the south of the cloud, indicating that G0.253+0.016 is colliding with another cloud at v(LSR)~70 km/s. We confirm via an analysis of the observed dynamics in the CMZ that it is an elongated structure, orientated with Sgr B2 closer to the Sun than Sgr A*, however our results suggest that the actual geometry may be more complex than an elliptical ring. We find that the column density PDF of G0.253+0.016 is log-normal with no discernible power-law tail, consistent with little star formation, and that its width can be explained in the framework of theory predicting the density structure of clouds created by supersonic, magnetised turbulence. We also present the delta-variance spectrum of this region, and show it is consistent with that expected for clouds with no star formation. Using G0.253+0.016 as a test-bed of the conditions required for star formation in a different physical environment to that of nearby clouds, we also conclude that there is not one column density threshold for star formation, but instead this value is dependant on the local physical conditions. [Abbrv.]

Context: The initial conditions for the gravitational collapse of molecular cloud cores and the subsequent birth of stars are still not well constrained. The characteristic cold temperatures (about 10 K) in such regions require observations at sub-millimetre and longer wavelengths. The Herschel Space Observatory and complementary ground-based observations presented in this paper have the unprecedented potential to reveal the structure and kinematics of a prototypical core region at the onset of stellar birth. Aims: This paper aims to determine the density, temperature, and velocity structure of the star-forming Bok globule CB 17. This isolated region is known to host (at least) two sources at different evolutionary stages: a dense core, SMM1, and a Class I protostar, IRS. Methods: We modeled the cold dust emission maps from 100 micron to 1.2 mm with both a modified blackbody technique to determine the optical depth-weighted line-of-sight temperature and column density and a ray-tracing technique to determine the core temperature and volume density structure. Furthermore, we analysed the kinematics of CB17 using the high-density gas tracer N2H+. Results: From the ray-tracing analysis, we find a temperature in the centre of SMM1 of 10.6 K, a flat density profile with radius 9500 au, and a central volume density of n(H) = 2.3x10^5 cm-3. The velocity structure of the N2H+ observations reveal global rotation with a velocity gradient of 4.3 km/s/pc. Superposed on this rotation signature we find a more complex velocity field, which may be indicative of differential motions within the dense core. Conclusions: SMM is a core in an early evolutionary stage at the verge of being bound, but the question of whether it is a starless or a protostellar core remains unanswered.

The mid- and far-infrared view on high-mass star formation, in particular with the results from the Herschel space observatory, has shed light on many aspects of massive star formation. However, these continuum studies lack kinematic information. We study the kinematics of the molecular gas in high-mass star-forming regions. We complemented the PACS and SPIRE far-infrared data of 16 high-mass star-forming regions from the Herschel key project EPoS with N2H+ molecular line data from the MOPRA and Nobeyama 45m telescope. Using the full N2H+ hyperfine structure, we produced column density, velocity, and linewidth maps. These were correlated with PACS 70micron images and PACS point sources. In addition, we searched for velocity gradients. For several regions, the data suggest that the linewidth on the scale of clumps is dominated by outflows or unresolved velocity gradients. IRDC18454 and G11.11 show two velocity components along several lines of sight. We find that all regions with a diameter larger than 1pc show either velocity gradients or fragment into independent structures with distinct velocities. The velocity profiles of three regions with a smooth gradient are consistent with gas flows along the filament, suggesting accretion flows onto the densest regions. We show that the kinematics of several regions have a significant and complex velocity structure. For three filaments, we suggest that gas flows toward the more massive clumps are present.

Spitzer Space Telescope and Herschel Space Observatory imaging of M31 is used, with a physical dust model, to construct maps of dust surface density, dust-to-gas ratio, starlight heating intensity, and PAH abundance, out to R=25kpc. The global dust mass is M_d=5.4x10^7Msol, the global dust/H mass ratio is M_d/M_H=0.0081, and the global PAH abundance is <q_PAH>=0.039. The dust surface density has an inner ring at R=5.6kpc, a maximum at R=11.2kpc, and an outer ring at R=15.1kpc. The dust/gas ratio varies from M_d/M_H=0.026 at the center to ~0.0027 at R=25kpc. From the dust/gas ratio, we estimate the ISM metallicity to vary by a factor ~10, from Z/Zsol=3 at R=0 to ~0.3 at R=25kpc. The dust heating rate parameter <U> peaks at the center, with <U> approx 35, declining to <U> approx 0.25 at R=20kpc. Within the central kpc, the starlight heating intensity inferred from the dust modeling is close to what is estimated from the stars in the bulge. The PAH abundance reaches a peak q_PAH=0.045 at R=11.2kpc. When allowance is made for the different spectrum of the bulge stars, q_PAH for the dust in the central kpc is similar to the overall value of q_PAH in the disk. The silicate-graphite-PAH dust model used here is generally able to reproduce the observed dust spectral energy distribution across M31, but overpredicts 500um emission at R=2-6kpc, suggesting that at R=2-6kpc, the dust opacity varies more steeply with frequency (with beta approx 2.3 between 200 and 600um) than in the model

We present a survey of molecules in a sample of Galactic center molecular clouds using the Karl G. Jansky Very large Array, which includes M0.25+0.01, the clouds near Sgr A, and Sgr B2. The molecules detected are primarily NH3 and HC3N; in Sgr B2-N we also detect nonmetastable NH3, vibrationally-excited HC3N, torsionally-excited CH3OH, and numerous isotopologues of these species. 36 GHz Class I CH3OH masers are ubiquitous in these fields, and in several cases are associated with new NH3 (3,3) maser candidates. We also find that NH3 and HC3N are depleted or absent toward several of the highest dust column density peaks identified in submillimeter observations, which are associated with water masers and are thus likely in the early stages of star formation.

Seven isolated, nearby low-mass starless molecular cloud cores have been observed as part of the Herschel key program Earliest Phases of Star formation (EPoS). By applying a ray-tracing technique to the obtained continuum emission and complementary (sub)mm emission maps, we derive the physical structure (density, dust temperature) of these cloud cores. We present observations of the 12CO, 13CO, and C18O (2-1) and N2H+ (1-0) transitions towards the same cores. Based on the density and temperature profiles, we apply time-dependent chemical and line-radiative transfer modeling and compare the modeled to the observed molecular emission profiles. CO is frozen onto the grains in the center of all cores in our sample. The level of CO depletion increases with hydrogen density and ranges from 46% up to more than 95% in the core centers in the core centers in the three cores with the highest hydrogen density. The average hydrogen density at which 50% of CO is frozen onto the grains is 1.1+-0.4 10^5 cm^-3. At about this density, the cores typically have the highest relative abundance of N2H+. The cores with higher central densities show depletion of N2H+ at levels of 13% to 55%. The chemical ages for the individual species are on average 2+-1 10^5 yr for 13CO, 6+-3 10^4 yr for C18O, and 9+-2 10^4 yr for N2H+. Chemical modeling indirectly suggests that the gas and dust temperatures decouple in the envelopes and that the dust grains are not yet significantly coagulated. We observationally confirm chemical models of CO-freezeout and nitrogen chemistry. We find clear correlations between the hydrogen density and CO depletion and the emergence of N2H+. The chemical ages indicate a core lifetime of less than 1 Myr.

We present kiloparsec (kpc) spatial resolution maps of the CO-to-H2 conversion factor (alpha_co) and dust-to-gas ratio (DGR) in 26 nearby, star-forming galaxies. We have simultaneously solved for alpha_co and DGR by assuming that the DGR is approximately constant on kpc scales. With this assumption, we can combine maps of dust mass surface density, CO integrated intensity and HI column density to solve for both alpha_co and DGR with no assumptions about their value or dependence on metallicity or other parameters. Such a study has just become possible with the availability of high resolution far-IR maps from the Herschel key program KINGFISH, 12CO J=(2-1) maps from the IRAM 30m large program HERACLES and HI 21-cm line maps from THINGS. We use a fixed ratio between the (2-1) and (1-0) lines to present our alpha_co results on the more typically used 12CO J=(1-0) scale and show using literature measurements that variations in the line ratio do not effect our results. In total, we derive 782 individual solutions for alpha_co and DGR. On average, alpha_co = 3.1 Msun pc^-2 (K km s^-1)^-1 for our sample with a standard deviation of 0.3 dex. Within galaxies we observe a generally flat profile of alpha_co as a function of galactocentric radius. However, most galaxies exhibit a lower alpha_co in the central kpc---a factor of ~2 below the galaxy mean, on average. In some cases, the central alpha_co value can be factors of 5 to 10 below the standard Milky Way (MW) value of alpha_co,MW =4.4 Msun pc^-2 (K km s^-1)^-1. While for alpha_co we find only weak correlations with metallicity, DGR is well-correlated with metallicity, with an approximately linear slope. Finally, we present several recommendations for choosing an appropriate alpha_co for studies of nearby galaxies.

We use mid-infrared Spitzer spectroscopy and far-infrared Herschel photometry for a sample of twenty main sequence O9--B2 stars in the Small Magellanic Cloud (SMC) with strong 24 micron excesses to investigate the origin of the mid-IR emission. Either debris disks around the stars or illuminated patches of dense interstellar medium (ISM) can cause such mid-IR emission. In a companion paper, Paper I, we use optical spectroscopy to show that it is unlikely for any of these sources to be classical Be stars or Herbig Ae/Be stars. We focus our analysis on debris disks and cirrus hot spots. We find three out of twenty stars to be significantly extended in the mid-IR, establishing them as cirrus hot spots. We then fit the IR spectral energy distributions to determine dust temperatures and masses. We find the dust masses in the SMC stars to be larger than for any known debris disks, although this evidence against the debris disk hypothesis is circumstantial. Finally, we created a local comparison sample of bright mid-IR OB stars in the Milky Way (MW) by cross-matching the WISE and Hipparcos catalogs. All such local stars in the appropriate luminosity range that can be unambiguously classified are young stars with optical emission lines or are spatially resolved by WISE with sizes too large to be plausible debris disk candidates. We conclude that the very strong mid-IR flux excesses are most likely explained as cirrus hot spots, although we cannot rigorously rule out that a small fraction of the sample is made up of debris disks or transition disks. We present suggestive evidence that bow-shock heating around runaway stars may be a contributing mechanism to the interstellar emission. These sources, interpreted as cirrus hot spots, offer a new localised probe of diffuse interstellar dust in a low metallicity environment. (Abridged)

Aims: We study the fragmentation and dynamical properties of a massive starless gas clump at the onset of high-mass star formation. Methods: Based on Herschel continuum data we identify a massive gas clump that remains far-infrared dark up to 100mum wavelengths. The fragmentation and dynamical properties are investigated by means of Plateau de Bure Interferometer and Nobeyama 45m single-dish spectral line and continuum observations. Results: The massive gas reservoir fragments at spatial scales of ~18000AU in four cores. Comparing the spatial extent of this high-mass region with intermediate- to low-mass starless cores from the literature, we find that linear sizes do not vary significantly over the whole mass regime. However, the high-mass regions squeeze much more gas into these similar volumes and hence have orders of magnitude larger densities. The fragmentation properties of the presented low-to high-mass regions are consistent with gravitational instable Jeans fragmentation. Furthermore, we find multiple velocity components associated with the resolved cores. Recent radiative transfer hydrodynamic simulations of the dynamic collapse of massive gas clumps also result in multiple velocity components along the line of sight because of the clumpy structure of the regions. This result is supported by a ratio between viral and total gas mass for the whole region <1. Conclusions: This apparently still starless high-mass gas clump exhibits clear signatures of early fragmentation and dynamic collapse prior to the formation of an embedded heating source. A comparison with regions of lower mass reveals that the linear size of star-forming regions does not necessarily have to vary much for different masses, however, the mass reservoirs and gas densities are orders of magnitude enhanced for high-mass regions compared to their lower-mass siblings.

The temperature and density structure of molecular cloud cores are the most important physical quantities that determine the course of the protostellar collapse and the properties of the stars they form. Nevertheless, density profiles often rely either on the simplifying assumption of isothermality or on observationally poorly constrained model temperature profiles. With the aim of better constraining the initial physical conditions in molecular cloud cores at the onset of protostellar collapse, we initiated the Guaranteed Time Key Project (GTKP) "The Earliest Phases of Star Formation" (EPoS) with the Herschel satellite. This paper gives an overview of the low-mass sources in the EPoS project, including all observations, the analysis method, and the initial results of the survey. We study the thermal dust emission of 12 previously well-characterized, isolated, nearby globules using FIR and submm continuum maps at up to eight wavelengths between 100 micron and 1.2 mm. Our sample contains both globules with starless cores and embedded protostars at different early evolutionary stages. The dust emission maps are used to extract spatially resolved SEDs, which are then fit independently with modified blackbody curves to obtain line-of-sight-averaged dust temperature and column density maps. We find that the thermal structure of all globules is dominated by external heating from the interstellar radiation field and moderate shielding by thin extended halos. All globules have warm outer envelopes (14-20 K) and colder dense interiors (8-12 K). The protostars embedded in some of the globules raise the local temperature of the dense cores only within radii out to about 5000 AU, but do not significantly affect the overall thermal balance of the globules.

We present continuum observations of the infrared dark cloud (IRDC) G48.66-0.22 (G48) obtained with Herschel, Spitzer, and APEX, in addition to several molecular line observations. The Herschel maps are used to derive temperature and column density maps of G48 using a model based on a modified blackbody. We find that G48 has a relatively simple structure and is relatively isolated; thus this IRDC provides an excellent target to study the collapse and fragmentation of a filamentary structure in the absence of complicating factors such as strong external feedback. The derived temperature structure of G48 is clearly non-isothermal from cloud to core scale. The column density peaks are spatially coincident with the lowest temperatures (~ 17.5 K) in G48. A total cloud mass of ~390Msun is derived from the column density maps. By comparing the luminosity-to-mass ratio of 13 point sources detected in the Herschel/PACS bands to evolutionary models, we find that two cores are likely to evolve into high-mass stars (M>8 Msun). The derived mean projected separation of point sources is smaller than in other IRDCs but in good agreement with theoretical predications for cylindrical collapse. We detect several molecular species such as CO, HCO+, HCN, HNC and N2H+. CO is depleted by a factor of ~3.5 compared to the expected interstellar abundance, from which we conclude that CO freezes out in the central region. Furthermore, the molecular clumps, associated with the sub-millimeter peaks in G48, appear to be gravitationally unbound or just pressure confined. The analysis of critical line masses in G48 show that the entire filament is collapsing, overcoming any internal support.

G18.93-0.03 is a prominent dust complex within an 0.8deg long filament, with the molecular clump G18.93/m being IR dark from near IR wavelength up to 160mu. Spitzer composite images show an IR bubble spatially associated with G18.93. We use GRS 13CO and IRAM 30m H13CO+ data to disentangle the spatial structure of the region. From ATLASGAL submm data we calculate the gas mass, while we use the H13CO+ line width to estimate its virial mass. Using HERSCHEL data we produce temperature maps from fitting the SED. With the MAGPIS 20cm and SuperCOSMOS Halpha data we trace the ionized gas, and the VGPS HI survey provides information on the atomic hydrogen gas. We show that the bubble is spatially associated with G18.93, located at a kinematic near distance of 3.6kpc. With 280Msun, the most massive clump within G18.93 is G18.93/m. The virial analysis shows that it may be gravitationally bound and has neither Spitzer young stellar objects nor mid-IR point sources within. Fitting the SED reveals a temperature distribution that decreases towards its center, but heating from the ionizing source puts it above the general ISM temperature. We find that the bubble is filled by HII gas, ionized by an O8.5 star. Between the ionizing source and the IR dark clump G18.93/m we find a layered structure, from ionized to atomic to molecular hydrogen, revealing a PDR. Furthermore, we identify an additional velocity component within the bubble's 8mu emission rim at the edge of the infrared dark cloud and speculate that it might be shock induced by the expanding HII region. While the elevated temperature allows for the build-up of larger fragments, and the shock induced velocity component may lead to additional turbulent support, we do not find conclusive evidence that the massive clump G18.93/m is prone to collapse because of the expanding HII region.

(Abriged) In the framework of the Herschel GTKP "The earliest phases of star formation", we have imaged B68 between 100 and 500 um. Ancillary (sub)mm data, spectral line maps of the 12/13CO(2-1) transitions as well as a NIR extinction map were added to the analysis. We employed a ray-tracing algorithm to derive the 2D mid-plane dust temperature and volume density distribution without suffering from LoS averaging effects of simple SED fitting procedures. Additional 3D radiative transfer calculations were employed to investigate the connection between the external irradiation and the peculiar crescent shaped morphology found in the FIR maps. For the first time, we spatially resolve the dust temperature and density distribution of B68. We find T_dust dropping from 16.7 K at the edge to 8.2 K in the centre, which is about 4 K lower than the result of the simple SED fitting approach. N_H peaks at 4.3x10^22 cm^-2 and n_H at 3.4x10^5 cm^-3 in the centre. B68 has a mass of 3.1 M_sun of material with A_K > 0.2 mag for an assumed distance of 150 pc. We detect a compact source in the southeastern trunk, which is also seen in extinction and CO. We find the radial density distribution from the edge of the inner plateau outward to be n_H ~ r^-3.5. Such a steep profile can arise from either or both of the following: external irradiation with a significant UV contribution or the fragmentation of filamentary structures. Our 3D radiative transfer model of an externally irradiated core by an anisotropic ISRF reproduces the crescent morphology. Our CO observations show that B68 is part of a chain of globules in both space and velocity, which may indicate that it was once part of a filament which dispersed. We also resolve a new compact source in the SE trunk and find that it is slightly shifted in centroid velocity from B68, lending qualitative support to core collision scenarios.

We aim at deriving the molecular abundances and temperatures of the hot molecular cores in the high-mass star-forming region NGC 6334I and consequently deriving their physical and astrochemical conditions. In the framework of the Herschel guaranteed time key program CHESS, NGC 6334I is investigated by using HIFI aboard the Herschel Space Observatory. A spectral line survey is carried out in the frequency range 480-1907 GHz, and auxiliary interferometric data from the SMA in the 230 GHz band provide spatial information for disentangling the different physical components contributing to the HIFI spectrum. The spectral lines are identified with the aid of former surveys and spectral line catalogs. The observed spectrum is then compared to a simulated synthetic spectrum with XCLASS, assuming local thermal equilibrium, and best fit parameters are derived using the model optimization package MAGIX. A total of 46 molecules are identified, with 31 isotopologues, resulting in about 4300 emission and absorption lines. High- energy levels of the dominant emitter methanol and vibrationally excited HCN are detected. The number of unidentified lines remains low with 75, or less than 2 percent of the lines detected. The modeling suggests that several spectral features need two or more components to be fitted properly. Other components could be assigned to cold foreground clouds or to outflows, most visible in the SiO emission. A chemical variation between the two embedded hot cores is found, with more N-bearing molecules identified in SMA1 and O-bearing molecules in SMA2. Spectral line surveys give powerful insights into the study of the interstellar medium. Different molecules trace different physical conditions like the inner hot core, the envelope, the outflows or the cold foreground clouds. The derived molecular abundances provide further constraints for astrochemical models.

(Abridged) We present an overview of the sample of high-mass star and cluster forming regions observed as part of the Earliest Phases of Star Formation (EPoS) Herschel Guaranteed Time Key Program. A sample of 45 infrared-dark clouds (IRDCs) were mapped at PACS 70, 100, and 160 micron and SPIRE 250, 350, and 500 micron. In this paper, we characterize a population of cores which appear in the PACS bands and place them into context with their host cloud and investigate their evolutionary stage. We construct spectral energy distributions (SEDs) of 496 cores which appear in all PACS bands, 34% of which lack counterparts at 24 micron. From single-temperature modified blackbody fits of the SEDs, we derive the temperature, luminosity, and mass of each core. These properties predominantly reflect the conditions in the cold, outer regions. Taking into account optical depth effects and performing simple radiative transfer models, we explore the origin of emission at PACS wavelengths. The core population has a median temperature of 20K and has masses and luminosities that span four to five orders of magnitude. Cores with a counterpart at 24 micron are warmer and bluer on average than cores without a 24 micron counterpart. We conclude that cores bright at 24 micron are on average more advanced in their evolution, where a central protostar(s) have heated the outer bulk of the core, than 24 micron-dark cores. The 24 micron emission itself can arise in instances where our line of sight aligns with an exposed part of the warm inner core. About 10% of the total cloud mass is found in a given cloud's core population. We uncover over 300 further candidate cores which are dark until 100 micron. These are candidate starless objects, and further observations will help us determine the nature of these very cold cores.

We use new Herschel multi-band imaging of the Andromeda galaxy to analyze how dust heating occurs in the central regions of galaxy spheroids that are essentially devoid of young stars. We construct a dust temperature map of M31 through fitting modified blackbody SEDs to the Herschel data, and find that the temperature within 2 kpc rises strongly from the mean value in the disk of 17 pm 1K to \sim35K at the centre. UV to near-IR imaging of the central few kpc shows directly the absence of young stellar populations, delineates the radial profile of the stellar density, and demonstrates that even the near-UV dust extinction is optically thin in M31's bulge. This allows the direct calculation of the stellar radiation heating in the bulge, U\ast(r), as a function of radius. The increasing temperature profile in the centre matches that expected from the stellar heating, i.e. that the dust heating and cooling rates track each other over nearly two orders of magnitude in U\ast. The modelled dust heating is in excess of the observed dust temperatures, suggesting that it is more than sufficient to explain the observed IR emission. Together with the wavelength dependent absorption cross section of the dust, this demonstrates directly that it is the optical, not UV, radiation that sets the heating rate. This analysis shows that neither young stellar populations nor stellar near-UV radiation are necessary to heat dust to warm temperatures in galaxy spheroids. Rather, it is the high densities of Gyr-old stellar populations that provide a sufficiently strong diffuse radiation field to heat the dust. To the extent which these results pertain to the tenuous dust found in the centres of early-type galaxies remains yet to be explored.

Context: The earliest stages of high-mass star formation are still poorly characterized. It is also unknown whether the initial conditions vary with environment. Aims: We want to investigate the youngest massive gas clumps in the environment of extremely active star formation. Methods: We selected the IRDC18454 complex associated with the W43 Galactic mini-starburst, and observed it in the continuum emission between 70mu and 1.2mm, and in spectral line emission of N2H+ and 13CO with single-dish instruments and an interferometer. Results: The multi-wavelength continuum study allows us to identify clumps that are infrared dark even at 70mu and hence the best candidates to be genuine high-mass starless gas clumps. The spectral energy distributions reveal elevated temperatures and luminosities compared to more quiescent environments. Furthermore, we identify a temperature gradient from the W43 mini-starburst toward the starless clumps. We discuss whether the radiation impact of the nearby mini-starburst changes the fragmentation properties of the gas clumps and by that maybe favors more high-mass star formation in such an environment. The spectral line data reveal two different velocity components of the gas at 100 and 50km/s. While chance projection is a possibility to explain these components, the projected associations of the emission sources as well as the prominent location at the Galactic bar -- spiral arm interface also allow the possibility that these two components may be spatially associated and even interacting. Conclusions: High-mass starless gas clumps can exist in the close environment of very active star formation without being destroyed. The impact of the active star formation sites may even allow for more high-mass stars to form in these 2nd generation gas clumps.