The density structure of the interstellar medium (ISM) determines where stars form and release energy, momentum, and heavy elements, driving galaxy evolution. Density variations are seeded and amplified by gas motion, but the exact nature of this motion is unknown across spatial scale and galactic environment. Although dense star-forming gas likely emerges from a combination of instabilities, convergent flows, and turbulence, establishing the precise origin is challenging because it requires quantifying gas motion over many orders of magnitude in spatial scale. Here we measure the motion of molecular gas in the Milky Way and in nearby galaxy NGC 4321, assembling observations that span an unprecedented spatial dynamic range ($10^{-1}{-}10^3$ pc). We detect ubiquitous velocity fluctuations across all spatial scales and galactic environments. Statistical analysis of these fluctuations indicates how star-forming gas is assembled. We discover oscillatory gas flows with wavelengths ranging from $0.3{-}400$ pc. These flows are coupled to regularly-spaced density enhancements that likely form via gravitational instabilities. We also identify stochastic and scale-free velocity and density fluctuations, consistent with the structure generated in turbulent flows. Our results demonstrate that ISM structure cannot be considered in isolation. Instead, its formation and evolution is controlled by nested, interdependent flows of matter covering many orders of magnitude in spatial scale.

Knowledge about the distribution of CO emission in the Milky Way is essential to understand the impact of Galactic environment on the formation and evolution of structures in the interstellar medium. However, currently our insight about the fraction of CO in spiral arm and interarm regions is still limited by large uncertainties in assumed rotation curve models or distance determination techniques. In this work we use a Bayesian approach to obtain the current best assessment of the distribution of 13CO from the Galactic Ring Survey. We performed two different distance estimates that either included or excluded a model for Galactic features. We also include a prior for the solution of the kinematic distance ambiguity that was determined from a compilation of literature distances and an assumed size-linewidth relationship. We find that the fraction of 13CO emission associated with spiral arm features varies from 76% to 84% between the two distance runs. The vertical distribution of the gas is concentrated around the Galactic midplane showing FWHM values of ~75 pc. We do not find any significant difference between gas emission properties associated with spiral arm and interarm features. In particular the distribution of velocity dispersion values of gas emission in spurs and spiral arms is very similar. We detect a trend of higher velocity dispersion values with increasing heliocentric distance, which we attribute to beam averaging effects caused by differences in spatial resolution. We argue that the true distribution of the gas emission is likely more similar to a combination of our two distance results, and highlight the importance of using complementary distance estimations to safeguard against the pitfalls of any single approach. We conclude that the methodology presented in this work is a promising way to determine distances to gas emission features in Galactic plane surveys.

CO isotopologue transitions are routinely observed in molecular clouds to probe the column density of the gas, the elemental ratios of carbon and oxygen, and to trace the kinematics of the environment. We aim at estimating the abundances, excitation temperatures, velocity field and velocity dispersions of the three main CO isotopologues towards a subset of the Orion B molecular cloud. We use the Cramer Rao Bound (CRB) technique to analyze and estimate the precision of the physical parameters in the framework of local-thermodynamic-equilibrium excitation and radiative transfer with an additive white Gaussian noise. We propose a maximum likelihood estimator to infer the physical conditions from the 1-0 and 2-1 transitions of CO isotopologues. Simulations show that this estimator is unbiased and efficient for a common range of excitation temperatures and column densities (Tex > 6 K, N > 1e14 - 1e15 cm-2). Contrary to the general assumptions, the different CO isotopologues have distinct excitation temperatures, and the line intensity ratios between different isotopologues do not accurately reflect the column density ratios. We find mean fractional abundances that are consistent with previous determinations towards other molecular clouds. However, significant local deviations are inferred, not only in regions exposed to UV radiation field but also in shielded regions. These deviations result from the competition between selective photodissociation, chemical fractionation, and depletion on grain surfaces. We observe that the velocity dispersion of the C18O emission is 10% smaller than that of 13CO. The substantial gain resulting from the simultaneous analysis of two different rotational transitions of the same species is rigorously quantified. The CRB technique is a promising avenue for analyzing the estimation of physical parameters from the fit of spectral lines.

Molecular clouds form from the atomic phase of the interstellar medium. However, characterizing the transition between the atomic and the molecular interstellar medium (ISM) is a difficult observational task. Here we address cloud formation processes by combining HSIA with molecular line data. One scenario proposed by numerical simulations is that the column density probability density functions (N-PDF) evolves from a log-normal shape at early times to a power-law-like shape at later times. In this paper, we study the cold atomic component of the giant molecular filament GMF38a (d=3.4 kpc, length$\sim230$ pc). We identify an extended HISA feature, which is partly correlated with the 13CO emission. The peak velocities of the HISA and 13CO observations agree well on the eastern side of the filament, whereas a velocity offset of approximately 4 km s$^{-1}$ is found on the western side. The sonic Mach number we derive from the linewidth measurements shows that a large fraction of the HISA, which is ascribed to the cold neutral medium (CNM), is at subsonic and transonic velocities. The column density of the CNM is on the order of 10$^{20}$ to 10$^{21}$ cm$^{-2}$. The column density of molecular hydrogen is an order of magnitude higher. The N-PDFs from HISA (CNM), HI emission (WNM+CNM), and 13CO (molecular component) are well described by log-normal functions, which is in agreement with turbulent motions being the main driver of cloud dynamics. The N-PDF of the molecular component also shows a power law in the high column-density region, indicating self-gravity. We suggest that we are witnessing two different evolutionary stages within the filament. The eastern subregion seems to be forming a molecular cloud out of the atomic gas, whereas the western subregion already shows high column density peaks, active star formation and evidence of related feedback processes.

The analysis of large molecular line surveys of the Galactic plane is essential for our understanding of the gas kinematics on Galactic scales, in particular its link with the formation and evolution of dense structures in the interstellar medium. An approximation of the emission peaks with Gaussian functions allows for an efficient and straightforward extraction of useful physical information contained in the shape and Doppler-shifted frequency of the emission lines contained in these enormous data sets. In this work we present an overview and first results of a Gaussian decomposition of the entire Galactic Ring Survey (GRS) 13CO (1-0) data that consists of about 2.3 million spectra. We performed the decomposition with the fully automated GaussPy+ algorithm and fitted about 4.6 million Gaussian components to the GRS spectra. We discuss the statistics of the fit components and relations between the fitted intensities, velocity centroids, and velocity dispersions. We find that the magnitude of the velocity dispersion values increase toward the inner Galaxy and around the Galactic midplane, which we speculate is partly due to the influence of the Galactic bar and regions with higher non-thermal motions located in the midplane, respectively. We also use our decomposition results to infer global properties of the gas emission and find that the number of fit components used per spectrum is indicative for the amount of structure along the line of sight. We find that the emission lines from regions located on the far side of the Galaxy show increased velocity dispersion values, likely due to beam averaging effects. We demonstrate how this trend has the potential to aid in characterising Galactic structure by disentangling emission that is belonging to the nearby Aquila Rift molecular cloud from emission that is more likely associated with the Perseus and Outer spiral arms.

The Vista Variables in the Via Lactea (VVV) survey has performed a multi-epoch near-infrared imaging of the inner Galactic plane. High-fidelity photometric catalogs are needed to utilize the data. We aim at producing a deep, point-spread-function (PSF) photometric catalog for the VVV survey J, H, and Ks band data. Specifically, we aim at taking advantage of all the epochs of the survey to reach high limiting magnitudes. We develop an automatic PSF-fitting pipeline based on the DaoPHOT algorithm and perform photometry on the stacked VVV images in J, H, and Ks bands. We present a PSF photometric catalog in the Vega system that contains about 926 million sources in the J, H, and Ks filters. About 10% of the sources are flagged as possible spurious detections. The 5 sigma limiting magnitudes of the sources with high reliability are about 20.8, 19.5, and 18.7 mag in the J, H, and Ks band, respectively, depending on the local crowding condition. Our photometric catalog reaches on average about one magnitude deeper than the previously released PSF DoPHOT photometric catalog. It also includes less spurious detections. There are significant differences in the brightnesses of faint sources between our catalog and the previously released one. The likely origin of these differences is in the different photometric algorithms that are utilized; it is not straightforward to assess which catalog is more accurate in which situations. Our new catalog is beneficial especially for science goals that require high limiting magnitudes; our catalog reaches such in fields that have a relatively uniform source number density. Overall, the limiting magnitudes and completeness are different in the fields with different crowding conditions.

Using a method based on Faraday rotation measurements, Tahani et al. 2018 find the line-of-sight component of magnetic fields in Orion-A and show that their direction changes from the eastern side of this filamentary structure to its western side. Three possible magnetic field morphologies that can explain this reversal across the Orion-A region are toroidal, helical, and bow-shaped morphologies. In this paper we construct simple models to represent these three morphologies and compare them with the available observational data to find the most probable morphology(ies). To compare the observations with the models, we use probability values and a Monte-Carlo analysis to determine the most likely magnetic field morphology among these three morphologies. We find that the bow morphology has the highest probability values and that our Monte-Carlo analysis suggests that the bow morphology is more likely. We suggest that the bow morphology is the most likely and the most natural of the three morphologies that could explain a magnetic field reversal across the Orion-A filamentary structure (i.e., bow, helical and toroidal morphologies).

Context: OH masers trace diverse physical processes, from the expanding envelopes around evolved stars to star-forming regions or supernovae remnants. Aims: We identify the ground-state OH masers at 18cm wavelength in the area covered by ``The HI/OH/Recombination line survey of the Milky Way (THOR)''. We present a catalogue of all OH maser features and their possible associated environments. Methods: The THOR survey covers longitude and latitude ranges of 14.3<l<66.8 and b<1.25 deg. All OH ground state lines at 1612, 1665, 1667 and 1720MHz have been observed, employing the Very Large Array. The spatial resolution of the data varies between 12.5'' and 19'', the spectral resolution is 1.5km/s, and the rms sensitivity of the data is ~10mJy/beam per channel. Results: We identify 1585 individual maser spots distributed over 807 maser sites. Based on different criteria from spectral profiles to literature comparison, we try to associate the maser sites with astrophysical source types. Approximately 51\% of the sites exhibit the double-horned 1612MHz spectra typically emitted from the expanding shells of evolved stars. The separations of the two main velocity features of the expanding shells typically vary between 22 and 38km/s. In addition to this, at least 20% of the maser sites are associated with star-forming regions. While the largest fraction of 1720MHz maser spots (21 out of 53) is associated with supernova remnants, a significant fraction of the 1720MHz maser spots (17) are also associated with star-forming regions. We present comparisons to the thermal 13CO(1-0) emission as well as to other surveys of class II CH3OH and H2O maser emission. The catalogue attempts to present associations to astrophysical sources where available, and the full catalogue is available in electronic form.

High-mass stars and star clusters commonly form within hub-filament systems. Monoceros R2, harbors one of the closest such systems, making it an excellent target for case studies. We investigate the morphology, stability and dynamical properties of the hub-filament system on basis of 13CO and C18O observations obtained with the IRAM-30m telescope and H2 column density maps derived from Herschel dust emission observations. We identified the filamentary network and characterized the individual filaments as either main (converging into the hub) or secondary (converging to a main filament) filaments. The main filaments have line masses of 30-100 Msun/pc and show signs of fragmentation. The secondary filaments have line masses of 12-60 Msun/pc and show fragmentation only sporadically. In the context of Ostriker's hydrostatic filament model, the main filaments are thermally super-critical. If non-thermal motions are included, most of them are trans-critical. Most of the secondary filaments are roughly trans-critical regardless of whether non-thermal motions are included or not. From the main filaments, we estimate a mass accretion rate of 10(-4)-10(-3) Msun/pc into the hub. The secondary filaments accrete into the main filaments with a rate of 0.1-0.4x10(-4) Msun/pc. The main filaments extend into the hub. Their velocity gradients increase towards the hub, suggesting acceleration of the gas. We estimate that with the observed infall velocity, the mass-doubling time of the hub is ~2.5 Myr, ten times larger than the free-fall time, suggesting a dynamically old region. These timescales are comparable with the chemical age of the HII region. Inside the hub, the main filaments show a ring- or a spiral-like morphology that exhibits rotation and infall motions. One possible explanation for the morphology is that gas is falling into the central cluster following a spiral-like pattern.

Our understanding of the dynamics of the interstellar medium is informed by the study of the detailed velocity structure of emission line observations. One approach to study the velocity structure is to decompose the spectra into individual velocity components; this leads to a description of the dataset that is significantly reduced in complexity. However, this decomposition requires full automation lest it becomes prohibitive for large datasets, such as Galactic plane surveys. We developed GaussPy+, a fully automated Gaussian decomposition package that can be applied to emission line datasets, especially large surveys of HI and isotopologues of CO. We built our package upon the existing GaussPy algorithm and significantly improved its performance for noisy data. New functionalities of GaussPy+ include: i) automated preparatory steps, such as an accurate noise estimation, which can also be used as standalone applications; ii) an improved fitting routine; iii) an automated spatial refitting routine that can add spatial coherence to the decomposition results by refitting spectra based on neighbouring fit solutions. We thoroughly tested the performance of GaussPy+ on synthetic spectra and a test field from the Galactic Ring Survey. We found that GaussPy+ can deal with cases of complex emission and even low to moderate signal-to-noise values.

Young massive clusters (YMCs) are the most compact, high-mass stellar systems still forming at the present day. The precursor clouds to such systems are, however, rare due to their large initial gas mass reservoirs and rapid dispersal timescales due to stellar feedback. Nonetheless, unlike their high-z counterparts, these precursors are resolvable down to the sites of individually forming stars, and hence represent the ideal environments in which to test the current theories of star and cluster formation. Using high angular resolution (1$^{\prime\prime}$ / 0.05pc) and sensitivity ALMA observations of two YMC progenitor clouds in the Galactic Centre, we have identified a suite of molecular line transitions -- e.g. c-C$_{3}$H$_{2} $($7-6$) -- that are believed to be optically thin, and reliably trace the gas structure in the highest density gas on star-forming core scales. We conduct a virial analysis of the identified core and proto-cluster regions, and show that half of the cores (5/10) and both proto-clusters are unstable to gravitational collapse. This is the first kinematic evidence of global gravitational collapse in YMC precursor clouds at such an early evolutionary stage. The implications are that if these clouds are to form YMCs, then they likely do so via the "conveyor-belt" mode, whereby stars continually form within dispersed dense gas cores as the cloud undergoes global gravitational collapse. The concurrent contraction of both the cluster-scale gas and embedded (proto)stars ultimately leads to the high (proto)stellar density in YMCs.

We introduce the histogram of oriented gradients (HOG), a tool developed for machine vision that we propose as a new metric for the systematic characterization of observations of atomic and molecular gas and the study of molecular cloud formation models. In essence, the HOG technique takes as input extended spectral-line observations from two tracers and provides an estimate of their spatial correlation across velocity channels. We characterize HOG using synthetic observations of HI and $^{13}$CO(J=1-0) emission from numerical simulations of MHD turbulence leading to the formation of molecular gas after the collision of two atomic clouds. We find a significant spatial correlation between the two tracers in velocity channels where $v_{HI}\approx v_{^{13}CO}$, independent of the orientation of the collision with respect to the line of sight. We use HOG to investigate the spatial correlation of the HI, from the THOR survey, and the $^{13}$CO(J=1-0) emission, from the GRS, toward the portion of the Galactic plane 33.75$\lt l\lt$35.25$^{o}$ and $|b|\lt$1.25$^{o}$. We find a significant spatial correlation between the tracers in extended portions of the studied region. Although some of the regions with high spatial correlation are associated with HI self-absorption features, suggesting that it is produced by the cold atomic gas, the correlation is not exclusive to this kind of region. The HOG results also indicate significant differences between individual regions: some show spatial correlation in channels around $v_{HI}\approx v_{^{13}CO}$ while others present this correlation in velocity channels separated by a few km/s. We associate these velocity offsets to the effect of feedback and to the presence of physical conditions that are not included in the atomic-cloud-collision simulations, such as more general magnetic field configurations, shear, and global gas infall.

Observations have discovered numerous giant molecular filaments (GMFs) in the Milky Way. However, their role in the Galactic star formation and Galaxy-scale evolution of dense gas is still unknown. We investigate systematically the star-forming content of all currently known GMFs. This allows us to estimate the star formation rates (SFRs) of the GMFs and to establish relationships between the SFRs and the GMF properties. We identify and classify the young stellar object (YSO) population of each GMF using multi-wavelength photometry from near- to far-infrared. We estimate the total SFRs assuming a universal and fully sampled initial mass function and luminosity function. We uniformly estimate the physical properties of 57 GMFs. The GMFs show correlations between the 13CO line width, mass, and size, similar to Larson's relations. We identify 36 394 infrared excess sources in 57 GMFs and obtain SFRs for 46 GMFs. The median SFR surface density (Sigma_SFR) and star formation efficiency (SFE) of GMFs are 0.62 Msun Myr-1 pc-2 and 1%, similar to the nearby star-forming clouds. The star formation rate per free-fall time of GMFs is between 0.002-0.05 with the median value of 0.02. We also find a strong correlation between SFR and dense gas mass that is defined as gas mass above a visual extinction of 7 mag, which suggests that the SFRs of the GMFs scale similarly with dense gas as those of nearby molecular clouds. We also find a strong correlation between the mean SFR per unit length and dense gas mass per unit length. The origin of this scaling remains unknown, calling for further studies that can link the structure of GMFs to their SF activity and explore the differences between GMFs and other molecular clouds

The hydroxyl radical (OH) is present in the diffuse molecular and partially atomic phases of the interstellar medium (ISM), but its abundance relative to hydrogen is not clear. We aim to evaluate the abundance of OH with respect to molecular hydrogen using OH absorption against cm-continuum sources over the first Galactic quadrant. This OH study is part of the HI/OH/Recombination line survey (THOR). THOR is a Karl G. Jansky Very Large Array large program of atomic, molecular and ionized gas in the range 15{\deg}$\leq$l$\leq$67{\deg} and |b|$\leq$1{\deg}. It is the highest-resolution unbiased OH absorption survey to date towards this region. We combine the derived optical depths with literature 13CO(1-0) and HI observations to determine the OH abundance. We detect absorption in the 1665 and 1667 MHz transitions for continuum sources stronger than $F_{\rm cont}\geq$0.1 Jy/beam. OH absorption is found against $\sim$15% of these continuum sources with increasing fractions for stronger sources. Most of the absorption is associated with Galactic HII regions. We find OH and 13CO gas to have similar kinematic properties. The OH abundance decreases with increasing hydrogen column density. The OH abundance with respect to the total hydrogen nuclei column density (atomic and molecular phase) is in agreement with a constant abundance for $A_V$ < 10-20. Towards the lowest column densities, we find sources that exhibit OH absorption but no 13CO emission, indicating that OH is a well suited tracer of the low column density molecular gas. We present spatially resolved OH absorption towards W43. The unbiased nature of the THOR survey opens a new window onto the gas properties of the ISM. The characterization of the OH abundance over a large range of hydrogen gas column densities contributes to the understanding of OH as a molecular gas tracer and provides a starting point for future investigations.

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.

An increasing number of hundred-parsec scale, high line-mass filaments have been detected in the Galaxy. Their evolutionary path, including fragmentation towards star formation, is virtually unknown. We characterize the fragmentation within the Nessie filament, covering size-scales between $\sim$ 0.1-100 pc. We also connect the small-scale fragments to the star-forming potential of the cloud. We combine near-infrared data from the VVV survey with mid-infrared GLIMPSE data to derive a high-resolution dust extinction map and apply a wavelet decomposition technique on it to analyze the fragmentation characteristics of the cloud, which are compared with predictions from fragmentation models. We compare the detected objects to those identified in $\sim$ 10 times coarser resolution from ATLASGAL data. We present a high-resolution extinction map of Nessie. We estimate the mean line-mass of Nessie to be $\sim$ 627 M$_\odot$/pc and the distance to be $\sim$ 3.5 kpc. We find that Nessie shows fragmentation at multiple size scales. The nearest-neighbour separations of the fragments at all scales are within a factor of 2 of the Jeans' length at that scale. However, the relationship between the mean densities of the fragments and their separations is significantly shallower than expected for Jeans' fragmentation. The relationship is similar to the one predicted for a filament that exhibits a Larson-like scaling between size-scale and velocity dispersion; such a scaling may result from turbulent support. Based on the number of YSOs in Nessie, we estimate that the star formation rate is $\sim$ 371 M$_\odot$/Myr; similar values result if using the number of dense cores, or the amount of dense gas, as the proxy of star formation. The star formation efficiency is 0.017. These numbers indicate that Nessie's star-forming content is comparable to the Solar neighborhood giant molecular clouds like Orion A.

Magnetic fields pervade in the interstellar medium (ISM) and are believed to be important in the process of star formation, yet probing magnetic fields in star formation regions is challenging. We propose a new method to use Faraday rotation measurements in small scale star forming regions to find the direction and magnitude of the component of magnetic field along the line-of-sight. We test the proposed method in four relatively nearby regions of Orion A, Orion B, Perseus, and California. We use rotation measure data from the literature. We adopt a simple approach based on relative measurements to estimate the rotation measure due to the molecular clouds over the Galactic contribution. We then use a chemical evolution code along with extinction maps of each cloud to find the electron column density of the molecular cloud at the position of each rotation measure data point. Combining the rotation measures produced by the molecular clouds and the electron column density, we calculate the line-of-sight magnetic field strength and direction. In California and Orion A, we find clear evidence that the magnetic fields at one side of these filamentary structures are pointing towards us and are pointing away from us at the other side. Even though the magnetic fields in Perseus might seem to suggest the same behavior, not enough data points are available to draw such conclusions. In Orion B, as well, there are not enough data points available to detect such behavior. This behavior is consistent with a helical magnetic field morphology. In the vicinity of available Zeeman measurements in OMC-1, OMC-B, and the dark cloud Barnard 1, we find magnetic field values of $-23\pm38~\mu$G, $-129\pm28~\mu$G, and $32\pm101~\mu$G, respectively, which are in agreement with the Zeeman Measurements.

The fragmentation of filaments in molecular clouds has attracted a lot of attention as there seems to be a relation between the evolution of filaments and star formation. The study of the fragmentation process has been motivated by simple analytical models. However, only a few comprehensive studies have analysed the evolution of filaments using numerical simulations where the filaments form self-consistently as part of molecular clouds. We address the early evolution of pc-scale filaments that form within individual clouds. We focus on three questions: How do the line masses of filaments evolve? How and when do the filaments fragment? How does the fragmentation relate to the line masses of the filaments? We examine three simulated molecular clouds formed in kpc-scale numerical simulations performed with the FLASH code. We compare the properties of the identified filaments with the predictions of analytic filament stability models. The line masses and mass fraction enclosed in the identified filaments increase continuously after the onset of self-gravity. The first fragments appear early when the line masses lie well below the critical line mass of Ostriker's hydrostatic equilibrium solution. The average line masses of filaments identified in 3D density cubes increases far more quickly than those identified in 2D column density maps. Our results suggest that hydrostatic or dynamic compression from the surrounding cloud has a significant impact on the early dynamical evolution of filaments. A simple model of an isolated, isothermal cylinder may not provide a good approach for fragmentation analysis. Caution must be exercised in interpreting distributions of properties of filaments identified in column density maps, especially in the case of low-mass filaments. Comparing or combining results from studies that use different filament finding techniques is strongly discouraged.

[abridged] We determine the relationship between turbulence energy and gas density variance for 15 molecular clouds in the Solar neighbourhood. We use the linewidths of the CO molecule as the probe of the turbulence energy (sonic Mach number, $\cal{M}_\mathrm{s}$) and three-dimensional models to reconstruct the density probability distribution function ($\rho$-PDF) of the clouds, derived using near-infrared extinction and Herschel dust emission data, as the probe of the density variance ($\sigma_\mathrm{s}$). We find no significant correlation between $\cal{M}_\mathrm{s}$ and $\sigma_\mathrm{s}$ among the studied clouds, however, we also cannot rule out a weak correlation. In the context of turbulence-dominated gas, the range of the $\cal{M}_\mathrm{s}$ and $\sigma_\mathrm{s}$ values corresponds with the model predictions. The data cannot constrain whether or not the turbulence driving parameter, $b$, and/or thermal-to-magnetic pressure ratio, $\beta$, vary among the sample clouds. Most clouds are not in agreement with field strengths stronger than given by $\beta \lesssim 0.05$. A model with $b^2 \beta / (\beta+1) = 0.30 \pm 0.06$ provides an adequate fit to the cloud sample as a whole. When considering the average behaviour of the sample, we can rule out three regimes: (i) strong compression combined with a weak magnetic field ($b \gtrsim 0.7$ and $\beta \gtrsim 3$), (ii) weak compression ($b \lesssim 0.35$), and (iii) strong magnetic field ($\beta \lesssim 0.1$). Including independent magnetic field strength estimates to the analysis, the data rule out solenoidal driving ($b < 0.4$) for the majority of the Solar neighbourhood clouds. However, most clouds have $b$ parameters larger than unity, which indicates a discrepancy with the turbulence-dominated picture; we discuss the possible reasons for this.

The early evolution during massive star cluster formation is still uncertain. Observing embedded clusters at their earliest stages of formation can provide insight into the spatial and temporal distribution of the stars and thus probe different star cluster formation models. We present near-infrared imaging of an 8'*13'(5.4pc*8.7pc) region around the massive infalling clump G286.21+0.17(also known as BYF73). The stellar content across the field is determined and photometry is derived in order to { obtain} stellar parameters for the cluster members. We find evidence for some sub-structure (on scales less than a pc diameter) within the region with apparently at least three different sub-clusters associated with the molecular clump based on differences in extinction and disk fractions. At the center of the clump we identify a deeply embedded sub-cluster. Near-infrared excess is detected for 39-44% in the two sub-clusters associated with molecular material and 27% for the exposed cluster. Using the disk excess as a proxy for age this suggests the clusters are very young. The current total stellar mass is estimated to be at least 200 Msun. The molecular core hosts a rich population of pre-main sequence stars. There is evidence for multiple events of star formation both in terms of the spatial distribution within the star forming region and possibly from the disk frequency.

Density profiles of isolated cores derived from thermal dust continuum emission rely on models of dust properties, such as mass opacity, which are poorly constrained. With complementary measures from near-infrared extinction maps, we can assess the reliability of commonly-used dust models. In this work, we compare Herschel-derived maps of the optical depth with equivalent maps derived from CFHT WIRCAM near-infrared observations for three isolated cores: CB68, L429, and L1552. We assess the dust opacities provided from four models: OH1a, OH5a, Orm1, and Orm4. Although the consistency of the models differs between the three sources, the results suggest that the optical properties of dust in the envelopes of the cores are best described by either silicate and bare graphite grains (e.g., Orm1) or carbonaceous grains with some coagulation and either thin or no ice mantles (e.g., OH5a). None of the models, however, individually produced the most consistent optical depth maps for every source. The results suggest that either the dust in the cores is not well described by any one dust property model, the application of the dust models cannot be extended beyond the very center of the cores, or more complex SED fitting functions are necessary.

We determine the physical properties and turbulence driving mode of molecular clouds formed in numerical simulations of a Milky Way-type disc galaxy with parsec-scale resolution. The clouds form through gravitational fragmentation of the gas, leading to average values for mass, radii and velocity dispersion in good agreement with observations of Milky Way clouds. The driving parameter (b) for the turbulence within each cloud is characterised by the ratio of the density contrast (sigma_rho) to the average Mach number (Mach) within the cloud, b = sigma_rho/Mach. As shown in previous works, b ~ 1/3 indicates solenoidal (divergence-free) driving and b ~ 1 indicates compressive (curl-free) driving. We find that the average b value of all the clouds formed in the simulations has a lower limit of b > 0.2. Importantly, we find that b has a broad distribution, covering values from purely solenoidal to purely compressive driving. Tracking the evolution of individual clouds reveals that the b value for each cloud does not vary significantly over their lifetime. Finally, we perform a resolution study with minimum cell sizes of 8, 4, 2 and 1 pc and find that the average b value increases with increasing resolution. Therefore, we conclude that our measured b values are strictly lower limits and that a resolution better than 1 pc is required for convergence. However, regardless of the resolution, we find that b varies by factors of a few in all cases, which means that the effective driving mode alters significantly from cloud to cloud.

We study the fragmentation of the nearest high line-mass filament, the integral shaped filament (ISF, line-mass $\sim$ 400 M$_\odot$ pc$^{-1}$) in the Orion A molecular cloud. We have observed a 1.6 pc long section of the ISF with the Atacama Large Millimetre/submillimeter Array (ALMA) at 3 mm continuum emission, at a resolution of $\sim$3" (1 200 AU). We identify from the region 43 dense cores with masses about a solar mass. 60% of the ALMA cores are protostellar and 40\% are starless. The nearest neighbour separations of the cores do not show a preferred fragmentation scale; the frequency of short separations increases down to 1 200 AU. We apply a two-point correlation analysis on the dense core separations and show that the ALMA cores are significantly grouped at separations below $\sim$17 000 AU and strongly grouped below $\sim$6 000 AU. The protostellar and starless cores are grouped differently: only the starless cores group strongly below $\sim$6 000 AU. In addition, the spatial distribution of the cores indicates periodic grouping of the cores into groups of $\sim$30 000 AU in size, separated by $\sim$50 000 AU. The groups coincide with dust column density peaks detected by Herschel. These results show hierarchical, two-mode fragmentation in which the maternal filament periodically fragments into groups of dense cores. Critically, our results indicate that the fragmentation models for lower line-mass filaments ($\sim$ 16 M$_\odot$ pc$^{-1}$) fail to capture the observed properties of the ISF. We also find that the protostars identified with Spitzer and Herschel in the ISF are grouped at separations below $\sim$17 000 AU. In contrast, young stars with disks do not show significant grouping. This suggests that the grouping of dense cores is partially retained over the protostar lifetime, but not over the lifetime of stars with disks.

Context: The past decade has witnessed a large number of Galactic plane surveys at angular resolutions below 20". However, no comparable high-resolution survey exists at long radio wavelengths around 21cm in line and continuum emission. Methods: Employing the Very Large Array (VLA) in the C-array configuration and a large program, we observe the HI 21cm line, four OH lines, nineteen Halpha radio recombination lines as well as the continuum emission from 1 to 2GHz in full polarization over a large part of the first Galactic quadrant. Results: Covering Galactic longitudes from 14.5 to 67.4deg and latitudes between +-1.25deg, we image all of these lines and the continuum at ~20" resolution. These data allow us to study the various components of the interstellar medium (ISM): from the atomic phase, traced by the HI line, to the molecular phase, observed by the OH transitions, to the ionized medium, revealed by the cm continuum and the Halpha radio recombination lines. Furthermore, the polarized continuum emission enables magnetic field studies. In this overview paper, we discuss the survey outline and present the first data release as well as early results from the different datasets. We now release the first half of the survey; the second half will follow later after the ongoing data processing has been completed. The data in fits format (continuum images and line data cubes) can be accessed through the project web-page http://www.mpia.de/thor. Conclusions: The HI/OH/Recombination line survey of the Milky Way (THOR) opens a new window to the different parts of the ISM. It enables detailed studies of molecular cloud formation, conversion of atomic to molecular gas, and feedback from HII regions as well as the magnetic field in the Milky Way. It is highly complementary to other surveys of our Galaxy, and comparing different datasets allows us to address many open questions.

We study the probability distribution function (PDF) of mass surface densities, $\Sigma$, of infrared dark cloud (IRDC) G028.37+00.07 and its surrounding giant molecular cloud. This PDF constrains the physical processes, such as turbulence, magnetic fields and self-gravity, that are expected to be controlling cloud structure and star formation activity. The chosen IRDC is of particular interest since it has almost 100,000 solar masses within a radius of 8 parsecs, making it one of the most massive, dense molecular structures known and is thus a potential site for the formation of a "super star cluster." We study $\Sigma$ in two ways. First, we use a combination of NIR and MIR extinction maps that are able to probe the bulk of the cloud structure up to $\Sigma\sim1\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq200$ mag). Second, we study the FIR and sub-mm dust continuum emission from the cloud utilizing Herschel PACS and SPIRE images and paying careful attention to the effects of foreground and background contamination. We find that the PDFs from both methods, applied over a $\sim20^\prime$(30 pc)-sized region that contains $\simeq1.5\times10^5\:M_\odot$ and encloses a minimum closed contour with $\Sigma\simeq0.013\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq3$ mag), shows a log-normal shape with the peak measured at $\Sigma\simeq0.021\:{\rm{g\:cm}^{-2}}\:$($A_V\simeq4.7$ mag). There is tentative evidence for the presence of a high-$\Sigma$ power law tail that contains from $\sim3\%$ to 8\% of the mass of the cloud material. We discuss the implications of these results for the physical processes occurring in this cloud.

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.

Context: In the last years, there have been many studies on the omnipresence and structures of filaments in star-forming regions, as well as their role in the process of star formation. Those filaments are normally identified as elongated fibres across the plane of the sky. But how would we detect filaments that are inclined? Aims: We aim to learn more about whether, and how, total column density or dust temperature change with respect to the line of sight. Such variations would enable observers to use dust observations to identify and study filaments at any inclination and gain more insight on the distribution and orientations of filaments within the Galactic plane. Methods: As a first step, we perform numerical calculations on simple cylindrical models to evaluate the influence of filament geometry on the average flux density. After that, we apply our three-dimensional Monte Carlo dust radiative transfer code on two models of star-forming regions and derive maps of effective total column density and dust temperature at different viewing angles. Results: We see only slight changes of average flux density for all cylinders we study. For our more complex models, we find that the effective dust temperature is not sensitive to viewing angle, while the total column density is strongly influenced, with differences exceeding an order of magnitude. The variations are not injective with the viewing angle and depend on the structure of the object. Conclusions: We conclude that there is no single quantity in our analysis that can uniquely trace the inclination and three-dimensional structure of a filament based on dust observations alone. However, observing wide variations in total column density at a given effective dust temperature is indicative of inclined filaments.

Filamentary structures are common morphological features of the cold, molecular interstellar medium (ISM). Recent studies have discovered massive, hundred-parsec-scale filaments that may be connected to the large-scale, Galactic spiral arm structure. Addressing the nature of these Giant Molecular Filaments (GMFs) requires a census of their occurrence and properties. We perform a systematic search of GMFs in the fourth Galactic quadrant and determine their basic physical properties. We identify GMFs based on their dust extinction signatures in near- and mid-infrared and velocity structure probed by ^{13}CO line emission. We use the ^{13}CO line emission and ATLASGAL dust emission data to estimate the total and dense gas masses of the GMFs. We combine our sample with an earlier sample from literature and study the Galactic environment of the GMFs. We identify nine GMFs in the fourth Galactic quadrant; six are located in the Centaurus spiral arm and three in inter-arm regions. Combining this sample with an earlier study using the same identification criteria in the first Galactic quadrant results in 16 GMFs, nine of which are located within spiral arms. The GMFs have sizes of 80-160 pc and ^{13}CO-derived masses between 5-90 x 10^{4} Msun. Their dense gas mass fractions are between 1.5-37%, being higher in the GMFs connected to spiral arms. We also compare the different GMF-identification methods and find that emission and extinction based techniques overlap only partially, highlighting the need to use both to achieve a complete census.

The dust extinction curve is a critical component of many observational programs and an important diagnostic of the physics of the interstellar medium. Here we present new measurements of the dust extinction curve and its variation towards tens of thousands of stars, a hundred-fold larger sample than in existing detailed studies. We use data from the APOGEE spectroscopic survey in combination with ten-band photometry from Pan-STARRS1, 2MASS, and WISE. We find that the extinction curve in the optical through infrared is well characterized by a one-parameter family of curves described by R(V). The extinction curve is more uniform than suggested in past works, with sigma(R(V)) = 0.18, and with less than one percent of sight lines having R(V) > 4. Our data and analysis have revealed two new aspects of Galactic extinction: first, we find significant, wide-area variations in R(V) throughout the Galactic plane. These variations are on scales much larger than individual molecular clouds, indicating that R(V) variations must trace much more than just grain growth in dense molecular environments. Indeed, we find no correlation between R(V) and dust column density up to E(B-V) ~ 2. Second, we discover a strong relationship between R(V) and the far-infrared dust emissivity.

Orion A hosts the nearest massive star factory, thus offering a unique opportunity to resolve the processes connected with the formation of both low- and high-mass stars. Here we present the most detailed and sensitive near-infrared (NIR) observations of the entire molecular cloud to date. With the unique combination of high image quality, survey coverage, and sensitivity, our NIR survey of Orion A aims at establishing a solid empirical foundation for further studies of this important cloud. In this first paper we present the observations, data reduction, and source catalog generation. To demonstrate the data quality, we present a first application of our catalog to estimate the number of stars currently forming inside Orion A and to verify the existence of a more evolved young foreground population. We used the European Southern Observatory's (ESO) Visible and Infrared Survey Telescope for Astronomy (VISTA) to survey the entire Orion A molecular cloud in the NIR $J, H$, and $K_S$ bands, covering a total of $\sim$18.3 deg$^2$. We implemented all data reduction recipes independently of the ESO pipeline. Estimates of the young populations toward Orion A are derived via the $K_S$-band luminosity function. Our catalog (799995 sources) increases the source counts compared to the Two Micron All Sky Survey by about an order of magnitude. The 90% completeness limits are 20.4, 19.9, and 19.0 mag in $J, H$, and $K_S$, respectively. The reduced images have 20% better resolution on average compared to pipeline products. We find between 2300 and 3000 embedded objects in Orion A and confirm that there is an extended foreground population above the Galactic field, in agreement with previous work. The Orion A VISTA catalog represents the most detailed NIR view of the nearest massive star-forming region and provides a fundamental basis for future studies of star formation processes toward Orion.

Filaments play a central role in the molecular clouds' evolution, but their internal dynamical properties remain poorly characterized. To further explore the physical state of these structures, we have investigated the kinematic properties of the Musca cloud. We have sampled the main axis of this filamentary cloud in $^{13}$CO and C$^{18}$O (2--1) lines using APEX observations. The different line profiles in Musca shows that this cloud presents a continuous and quiescent velocity field along its $\sim$6.5 pc of length. With an internal gas kinematics dominated by thermal motions (i.e., $\sigma_{NT}/c_s\lesssim1$) and large-scale velocity gradients, these results reveal Musca as the longest velocity-coherent, sonic-like object identified so far in the ISM. The transonic properties of Musca present a clear departure from the predicted supersonic velocity dispersions expected in the Larson's velocity dispersion-size relationship, and constitute the first observational evidence of a filament fully decoupled from the turbulent regime over multi-parsec scales.

Aims: We resolve the length-scales for filament formation and fragmentation (res. <=0.1pc), in particular the Jeans length and cylinder fragmentation scale. Methods: We observed the prototypical high-mass star-forming filament IRDC18223 with the Plateau de Bure Interferometer (PdBI) in the 3.2mm continuum and N2H+(1-0) line emission in a ten field mosaic at a spatial resolution of ~4'' (~14000AU). Results: The dust continuum emission resolves the filament into a chain of at least 12 relatively regularly spaced cores. The mean separation between cores is ~0.40(+-0.18)pc. While this is approximately consistent with the fragmentation of an infinite, isothermal, gravitationally bound gas cylinder, a high mass-to-length ratio of M/l~1000M_sun/pc requires additional turbulent and/or magnetic support against radial collapse of the filament. The N2H+(1-0) data reveal a velocity gradient perpendicular to the main filament. Although rotation of the filament cannot be excluded, the data are also consistent with the main filament being comprised of several velocity-coherent sub-filaments. Furthermore, this velocity gradient perpendicular to the filament resembles recent results toward Serpens south that are interpreted as signatures of filament formation within magnetized and turbulent sheet-like structures. Lower-density gas tracers ([CI] and C18O) reveal a similar red/blueshifted velocity structure on scales around 60'' east and west of the IRDC18223 filament. This may tentatively be interpreted as a signature of the large-scale cloud and the smaller-scale filament being kinematically coupled. We do not identify a velocity gradient along the axis of the filament. This may either be due to no significant gas flows along the filamentary axis, but it may partly also be caused by a low inclination angle of the filament with respect to the plane of the sky that could minimize such signature.

Filamentary structures are common in molecular clouds. Explaining how they fragment to dense cores is a missing step in understanding their role in star formation. We perform a case study of whether low-mass filaments are close-to hydrostatic prior to their fragmentation, and whether their fragmentation agrees with gravitational fragmentation models. For this, we study the 6.5 pc long Musca molecular cloud that is an ideal candidate for a filament at an early stage of fragmentation. We employ dust extinction mapping in conjunction with near-infrared data from the NEWFIRM instrument, and 870 um dust continuum emission data from the LABOCA instrument, to estimate column densities. We use the data to identify fragments from the cloud and to determine the radial density distribution of its filamentary part. We compare the cloud's morphology with 13CO and C18O line emission observed with the APEX/SHeFI instrument. The Musca cloud is pronouncedly fragmented at its ends, but harbours a remarkably well-defined, 1.6 pc long filament in its Center region. The line mass of the filament is 21-31 Ms pc^-1 and FWHM 0.07 pc. Its radial profile can be fitted with a Plummer profile that has the power-index of 2.6 \pm 11%, flatter than that of an infinite hydrostatic filament. The profile can also be fitted with a hydrostatic cylinder truncated by external pressure. These models imply a central density of 5-10 x 10^4 cm^-3. The fragments in the cloud have a mean separation of 0.4 pc, in agreement with gravitational fragmentation. These properties, together with the subsonic and velocity-coherent nature of the cloud, suggest a scenario in which an initially hydrostatic cloud is currently gravitationally fragmenting. The fragmentation has started a few tenths of a Myr ago from the cloud ends, leaving its center yet relatively non-fragmented, possibly because of gravitational focusing in a finite geometry.

We present the first study of the relationship between the column density distribution of molecular clouds within nearby Galactic spiral arms and their evolutionary status as measured from their stellar content. We analyze a sample of 195 molecular clouds located at distances below 5.5 kpc, identified from the ATLASGAL 870 micron data. We define three evolutionary classes within this sample: starless clumps, star-forming clouds with associated young stellar objects, and clouds associated with HII regions. We find that the N(H2) probability density functions (N-PDFs) of these three classes of objects are clearly different: the N-PDFs of starless clumps are narrowest and close to log-normal in shape, while star-forming clouds and HII regions exhibit a power-law shape over a wide range of column densities and log-normal-like components only at low column densities. We use the N-PDFs to estimate the evolutionary time-scales of the three classes of objects based on a simple analytic model from literature. Finally, we show that the integral of the N-PDFs, the dense gas mass fraction, depends on the total mass of the regions as measured by ATLASGAL: more massive clouds contain greater relative amounts of dense gas across all evolutionary classes.

We compare the structure of star-forming molecular clouds in different regions of Orion A to determine how the column density probability distribution function (N-PDF) varies with environmental conditions such as the fraction of young protostars. A correlation between the N-PDF slope and Class 0 protostar fraction has been previously observed in a low-mass star-formation region (Perseus) by Sadavoy; here we test if a similar correlation is observed in a high-mass star-forming region. We use Herschel data to derive a column density map of Orion A. We use the Herschel Orion Protostar Survey catalog for accurate identification and classification of the Orion A young stellar object (YSO) content, including the short-lived Class 0 protostars (with a $\sim$ 0.14 Myr lifetime). We divide Orion A into eight independent 13.5 pc$^2$ regions; in each region we fit the N-PDF distribution with a power-law, and we measure the fraction of Class 0 protostars. We use a maximum likelihood method to measure the N-PDF power-law index without binning. We find that the Class 0 fraction is higher in regions with flatter column density distributions. We test the effects of incompleteness, YSO misclassification, resolution, and pixel-scale. We show that these effects cannot account for the observed trend. Our observations demonstrate an association between the slope of the power-law N-PDF and the Class 0 fractions within Orion A. Various interpretations are discussed including timescales based on the Class 0 protostar fraction assuming a constant star-formation rate. The observed relation suggests that the N-PDF can be related to an "evolutionary state" of the gas. If universal, such a relation permits an evaluation of the evolutionary state from the N-PDF power-law index at much greater distances than those accesible with protostar counts. (abridged)

Throughout the Milky Way, molecular clouds typically appear filamentary, and mounting evidence indicates that this morphology plays an important role in star formation. What is not known is to what extent the dense filaments most closely associated with star formation are connected to the surrounding diffuse clouds up to arbitrarily large scales. How are these cradles of star formation linked to the Milky Way's spiral structure? Using archival Galactic plane survey data, we have used multiple datasets in search of large-scale, velocity-coherent filaments in the Galactic plane. In this paper, we present our methods employed to identify coherent filamentary structures first in extinction and confirmed using Galactic Ring Survey data. We present a sample of seven Giant Molecular Filaments (GMFs) that have lengths of order $\sim$100 pc, total masses of 10$^4$ - 10$^5$ M$_{\odot}$, and exhibit velocity coherence over their full length. The GMFs we study appear to be inter-arm clouds and may be the Milky Way analogues to spurs observed in nearby spiral galaxies. We find that between 2 and 12% of the total mass (above $\sim$10$^{20}$ cm$^{-2}$) is "dense" (above 10$^{22}$ cm$^{-2}$), where filaments near spiral arms in the Galactic midplane tend to have higher dense gas mass fractions than those further from the arms.

The formation of stars shapes the structure and evolution of entire galaxies. The rate and efficiency of this process are affected substantially by the density structure of the individual molecular clouds in which stars form. The most fundamental measure of this structure is the probability density function of volume densities (rho-PDF), which determines the star formation rates predicted with analytical models. This function has remained unconstrained by observations. We have developed an approach to quantify rho-PDFs and establish their relation to star formation. The rho-PDFs instigate a density threshold of star formation and allow us to quantify the star formation efficiency above it. The rho-PDFs provide new constraints for star formation theories and correctly predict several key properties of the star-forming interstellar medium.

We study how the estimation of the sonic Mach number ($M_s$) from $^{13}$CO linewidths relates to the actual 3D sonic Mach number. For this purpose we analyze MHD simulations which include post-processing to take radiative transfer effects into account. As expected, we find very good agreement between the linewidth estimated sonic Mach number and the actual sonic Mach number of the simulations for optically thin tracers. However, we find that opacity broadening causes $M_s$ to be overestimated by a factor of ~ 1.16 - 1.3 when calculated from optically thick $^{13}$CO lines. We also find that there is a dependency on the magnetic field: super-Alfv\'enic turbulence shows increased line broadening as compared with sub-Alfv\'enic turbulence for all values of optical depth for supersonic turbulence. Our results have implications for the observationally derived sonic Mach number--density standard deviation ($\sigma_{\rho/<\rho>}$) relationship, $\sigma^2_{\rho/<\rho>}=b^2M_s^2$, and the related column density standard deviation ($\sigma_{N/<N>}$) sonic Mach number relationship. In particular, we find that the parameter b, as an indicator of solenoidal vs. compressive driving, will be underestimated as a result of opacity broadening. We compare the $\sigma_{N/<N>}-M_s$ relation derived from synthetic dust extinction maps and $^{13}$CO linewidths with recent observational studies and find that solenoidally driven MHD turbulence simulations have values of $\sigma_{N/<N>}$ which are lower than real molecular clouds. This may be due to the influence of self-gravity which should be included in simulations of molecular cloud dynamics.

We use deep $8\:\mu m$ Spitzer-IRAC imaging of a massive Infrared Dark Cloud (IRDC) G028.37+00.07 to construct a Mid-Infrared (MIR) extinction map that probes mass surface densities up to $\Sigma\:\sim 1\:\rm{g~cm^{-2}}$ ($A_V\sim200\:$mag), amongst the highest values yet probed by extinction mapping. Merging with a NIR extinction map of the region, creates a high dynamic range map that reveals structures down to $A_V\sim1\:$mag. We utilize the map to: (1) Measure a cloud mass $\sim7\times10^4\:M_\odot$ within a radius of $\sim8\:$pc. $^{13}$CO kinematics indicate that the cloud is gravitationally bound. It thus has the potential to form one of the most massive young star clusters known in the Galaxy. (2) Characterize the structures of 16 massive cores within the IRDC, finding they can be fit by singular polytropic spheres with $\rho\propto{r}^{-k_\rho}$ and $k_\rho=1.3\pm0.3$. They have $\overline{\Sigma}\simeq0.1-0.4\:\rm{g~cm^{-2}}$ --- relatively low values that, along with their measured cold temperatures, suggest magnetic fields, rather than accretion-powered radiative heating, are important for controlling fragmentation of these cores. (3) Determine the $\Sigma$ (equivalently column density or $A_V$) probability distribution function (PDF) for a region that is near complete for $A_V>3\:$mag. The PDF is well fit by a single log-normal with mean $\overline{A}_V\simeq9\:$mag, high compared to other known clouds. It does not exhibit a separate high-end power law tail, which has been claimed to indicate the importance of self-gravity. However, we suggest that the PDF does result from a self-similar, self-gravitating hierarchy of structure being present over a wide range of scales in the cloud.

Some theories of dense molecular cloud formation involve dynamical environments driven by converging atomic flows or collisions between preexisting molecular clouds. The determination of the dynamics and physical conditions of the gas in clouds at the early stages of their evolution is essential to establish the dynamical imprints of such collisions, and to infer the processes involved in their formation. We present multi-transition 13CO and C18O maps toward the IRDC G035.39-00.33, believed to be at the earliest stages of evolution. The 13CO and C18O gas is distributed in three filaments (Filaments 1, 2 and 3), where the most massive cores are preferentially found at the intersecting regions between them. The filaments have a similar kinematic structure with smooth velocity gradients of ~0.4-0.8 km s-1 pc-1. Several scenarios are proposed to explain these gradients, including cloud rotation, gas accretion along the filaments, global gravitational collapse, and unresolved sub-filament structures. These results are complemented by HCO+, HNC, H13CO+ and HN13C single-pointing data to search for gas infall signatures. The 13CO and C18O gas motions are supersonic across G035.39-00.33, with the emission showing broader linewidths toward the edges of the IRDC. This could be due to energy dissipation at the densest regions in the cloud. The average H2 densities are ~5000-7000 cm-3, with Filaments 2 and 3 being denser and more massive than Filament 1. The C18O data unveils three regions with high CO depletion factors (f_D~5-12), similar to those found in massive starless cores.

The orientations of filamentary molecular clouds in the Gould Belt and their local ICM (inter-cloud media) magnetic fields are studied using near-infrared dust extinction maps and optical stellar polarimetry data. These filamentary clouds are a few-to-ten parsecs in length, and we find that their orientations tend to be either parallel or perpendicular to the mean field directions of the local ICM. This bimodal distribution is not found in cloud simulations with super-Alfvenic turbulence, in which the cloud orientations should be random. ICM magnetic fields that are dynamically important compared to inertial-range turbulence and self-gravity can readily explain both field-filament configurations. Previous studies commonly recognize that strong magnetic fields can guide gravitational contraction and result in filaments perpendicular to them, but few discuss the fact that magnetic fields can also channel sub-Alfvenic turbulence to form filaments aligned with them. This strong-field scenario of cloud formation is also consistent with the constant field strength observed from ICM to clouds (Crutcher et al. 2010) and is possible to explain the "hub-filament" cloud structure (Myers 2009) and the density threshold of cloud gravitational contraction (Kainulainen et al. 2009).

Star formation in molecular clouds is intimately linked to their internal mass distribution. We present an unprecedentedly detailed analysis of the column density structure of a high-mass, filamentary molecular cloud, namely IRDC G11.11-0.12 (G11). We use two novel column density mapping techniques: high-resolution (FWHM=2", or ~0.035 pc) dust extinction mapping in near- and mid-infrared, and dust emission mapping with the Herschel satellite. These two completely independent techniques yield a strikingly good agreement, highlighting their complementarity and robustness. We first analyze the dense gas mass fraction and linear mass density of G11. We show that G11 has a top heavy mass distribution and has a linear mass density (M_l ~ 600 Msun pc^{-1}) that greatly exceeds the critical value of a self-gravitating, non-turbulent cylinder. These properties make G11 analogous to the Orion A cloud, despite its low star-forming activity. This suggests that the amount of dense gas in molecular clouds is more closely connected to environmental parameters or global processes than to the star-forming efficiency of the cloud. We then examine hierarchical fragmentation in G11 over a wide range of size-scales and densities. We show that at scales 0.5 pc > l > 8 pc, the fragmentation of G11 is in agreement with that of a self-gravitating cylinder. At scales smaller than l < 0.5 pc, the results agree better with spherical Jeans' fragmentation. One possible explanation for the change in fragmentation characteristics is the size-scale-dependent collapse time-scale that results from the finite size of real molecular clouds: at scales l < 0.5 pc, fragmentation becomes sufficiently rapid to be unaffected by global instabilities.

We examine the physical parameters that affect the accumulation of gas in molecular clouds to high column densities where the formation of stars takes place. In particular, we analyze the dense gas mass fraction (DGMF) in a set of self-gravitating, isothermal, magnetohydrodynamic turbulence simulations including sink particles to model star formation. We find that the simulations predict close to exponential DGMFs over the column density range N(H2) = 3-25 x 10^{21} cm^{-2} that can be easily probed via, e.g., dust extinction measurements. The exponential slopes correlate with the type of turbulence driving and also with the star formation efficiency. They are almost uncorrelated with the sonic Mach number and magnetic-field strength. The slopes at early stages of cloud evolution are steeper than at the later stages. A comparison of these predictions with observations shows that only simulations with relatively non-compressive driving (b ~< 0.4) agree with the DGMFs of nearby molecular clouds. Massive infrared dark clouds can show DGMFs that are in agreement with more compressive driving. The DGMFs of molecular clouds can be significantly affected by how compressive the turbulence is on average. Variations in the level of compression can cause scatter to the DGMF slopes, and some variation is indeed necessary to explain the spread of the observed DGMF slopes. The observed DGMF slopes can also be affected by the clouds' star formation activities and statistical cloud-to-cloud variations.

We present a spectral investigation of X-ray binaries in NGC 5128 (Cen A), using six 100 ks Chandra observations taken over two months in 2007. We divide our sample into thermally and non-thermally dominated states based on the behavior of the fitted absorption column, and present the spectral parameters of sources with L >2x10^37 erg/s. The majority of sources are consistent with being neutron star low mass X-ray binaries (NS LMXBs) and we identify three transient black hole (BH) LMXB candidates coincident with the dust lane, which is the remnant of a small late-type galaxy. Our results also provide tentative support for the apparent `gap' in the mass distribution of compact objects between ~2-5 Msol. We propose that BH LMXBs are preferentially found in the dust lane, and suggest this is because of the younger stellar population. The majority (~70-80%) of potential Roche-lobe filling donors in the Cen A halo are >12 Gyr old, while BH LMXBs require donors >1 Msol to produce the observed peak luminosities. This requirement for more massive donors may also explain recent results that claim a steepening of the X-ray luminosity function with age at Lx >= 5x10^38 erg/s for the XB population of early-type galaxies; for older stellar populations, there are fewer stars >1 Msol, which are required to form the more luminous sources.

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.

Measuring the mass distribution of infrared dark clouds (IRDCs) over the wide dynamic range of their column densities is a fundamental obstacle in determining the initial conditions of high-mass star formation and star cluster formation. We present a new technique to derive high-dynamic-range, arcsecond-scale resolution column density data for IRDCs and demonstrate the potential of such data in measuring the density variance - sonic Mach number relation in molecular clouds. We combine near-infrared data from the UKIDSS/Galactic Plane Survey with mid-infrared data from the Spitzer/GLIMPSE survey to derive dust extinction maps for a sample of ten IRDCs. We then examine the linewidths of the IRDCs using 13CO line emission data from the FCRAO/Galactic Ring Survey and derive a column density - sonic Mach number relation for them. For comparison, we also examine the relation in a sample of nearby molecular clouds. The presented column density mapping technique provides a very capable, temperature independent tool for mapping IRDCs over the column density range equivalent to A_V=1-100 mag at a resolution of 2". Using the data provided by the technique, we present the first direct measurement of the relationship between the column density dispersion, \sigma_{N/<N>}, and sonic Mach number, M_s, in molecular clouds. We detect correlation between the variables with about 3-sigma confidence. We derive the relation \sigma_{N/<N>} = (0.047 \pm 0.016) Ms, which is suggestive of the correlation coefficient between the volume density and sonic Mach number, \sigma_{\rho/<\rho>} = (0.20^{+0.37}_{-0.22}) Ms, in which the quoted uncertainties indicate the 3-sigma range. When coupled with the results of recent numerical works, the existence of the correlation supports the picture of weak correlation between the magnetic field strength and density in molecular clouds (i.e., B ~ \rho^{0.5}).

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

The initial conditions of massive star and star cluster formation are expected to be cold, dense and high column density regions of the interstellar medium, which can reveal themselves via near, mid and even far-infrared absorption as Infrared Dark Clouds (IRDCs). Elucidating the dynamical state of IRDCs thus constrains theoretical models of these complex processes. In particular, it is important to assess whether IRDCs have reached virial equilibrium, where the internal pressure balances that due to the self-gravitating weight of the cloud plus the pressure of the external environmental. We study this question for the filamentary IRDC G035.39-00.33 by deriving mass from combined NIR & MIR extinction maps and velocity dispersion from C18O (1-0) & (2-1) line emission. In contrast to our previous moderately super-virial results based on 13CO emission and MIR-only extinction mapping, with improved mass measurements we now find that the filament is consistent with being in virial equilibrium, at least in its central parsec-wide region where ~1000 M_Sun snakes along several parsecs. This equilibrium state does not require large-scale net support or confinement by magnetic fields.

We report the discovery of a bright X-ray transient, CXOU J132527.6-430023, in the nearby early-type galaxy NGC 5128. The source was first detected over the course of five Chandra observations in 2007, reaching an unabsorbed outburst luminosity of 1-2*10^38 erg/s in the 0.5-7.0 keV band before returning to quiescence. Such luminosities are possible for both stellar-mass black hole and neutron star X-ray binary transients. Here, we attempt to characterize the nature of the compact object. No counterpart has been detected in the optical or radio sky, but the proximity of the source to the dust lanes allows for the possibility of an obscured companion. The brightness of the source after a >100 fold increase in X-ray flux makes it either the first confirmed transient non-ULX black hole system in outburst to be subject to detailed spectral modeling outside the Local Group, or a bright (>10^38 erg/s) transient neutron star X-ray binary, which are very rare. Such a large increase in flux would appear to lend weight to the view that this is a black hole transient. X-ray spectral fitting of an absorbed power law yielded unphysical photon indices, while the parameters of the best-fit absorbed disc blackbody model are typical of an accreting ~10 Msol black hole in the thermally dominant state.

Context: Infrared Dark Clouds (IRDCs) harbor progenitors of high-mass stars. Little is known of the parental molecular clouds of the IRDCs. Aims: We demonstrate the feasibility of the near-infrared (NIR) dust extinction mapping in tracing the parental molecular clouds of IRDCs at the distances of D = 2.5 - 8 kpc. Methods: We derive NIR extinction maps for 10 prominent IRDC complexes using a color-excess mapping technique and NIR data from the UKIDSS/Galactic Plane Survey. We compare the resulting maps to the 13CO emission line data, to the 8 \mu m dust opacity data, and to the millimeter dust emission data. We derive distances for the clouds by comparing the observed NIR source densities to the Besancon stellar distribution model and compare them to the kinematic distance estimates. Results: The NIR extinction maps provide a view to the IRDC complexes over the dynamical range of Av = 2 - 40 mag, in spatial resolution of 30". The NIR extinction data correlate well with the 13CO data and probe a similar gas component, but also extend to higher column densities. The NIR data reveal a wealth of extended structures surrounding the dense gas traced by the 8 \mu m shadowing features and sub-mm dust emission, showing that the clouds contain typically > 10 times more mass than traced by those tracers. The IRDC complexes of our sample contain relatively high amount of high-column density material, and their cumulative column density distributions resemble active nearby star-forming clouds like Orion rather than less active clouds like California. Conclusions: NIR dust extinction data provide a new powerful tool to probe the mass distribution of the parental molecular clouds of IRDCs up to the distances of D = 8 kpc. This encourages for deeper NIR observations of IRDCs, because the sensitivity and resolution of the data can be directly enhanced with dedicated observations.

Context: The large Coalsack dark cloud is one of the rare molecular clouds without clear signs of star formation. Aims: We investigate the dynamical properties of the gas within the Coalsack. Methods: The two highest extinction regions were mapped with the APEX telescope in 13CO(2-1) comprising a region of ~1 square degree. Results: In addition to the well-known, nearby gas component around -4km/s, we identified additional molecular gas components - in particular a second extended molecular cloud at a velocity of ~-30km/s and an estimated distance of ~3.1kpc - that dominate the column density and visual extinction distributions in the northeastern part of the Coalsack. Although comprising ~2600Msun, the mass of this distant cloud is distributed over an extent of ~73pc, much larger than typical high-mass infrared dark clouds. Its filamentary structure is consistent with a compressible gaseous self-gravitating cylinder, and its low mass per length indicates that it may be stable against gravitational collapse. We find barely any mid-infrared emission in archival MSX data, which is indicative of almost no star-formation activity in the near and far cloud complexes. The nearby clouds have narrow, almost thermal velocity dispersions with median values between 0.2 and 0.4km/s, which is also consistent with low star-formation activity. Only Tapia's Globule 2 exhibits a velocity dispersion increase toward the extinction peak and peak-velocity gradients over the core, which is indicative of a state of elevated dynamical properties. Conclusions: The Coalsack is not one single coherent structure, but consists of several cloud complexes nearby as well as at several kpc distance. All studied clouds appear as starless low-turbulence regions that may not even collapse in the future. Only one globule exhibits more dynamical signatures and is a good candidate for present/future star formation.

We present new optical spectroscopy with FLAMES/VLT, near-IR imaging with HAWK-I/VLT, and 870 micron mapping with APEX/LABOCA of the Coronet cluster. The optical data allow to estimate spectral types, extinction and the presence of accretion in 6 more M-type members, in addition to the 12 that we had previously studied. The submillimeter maps and near-IR data reveal the presence of nebular structures and high extinction regions, which are in some cases associated to known IR, optical, and X-ray sources. Most star formation is associated to two elongated structures crossing in the central part of the cluster. Placing all the 18 objects with known spectral types and extinction in the HR diagram suggests that the cluster is younger than previously thought (<2 Myr, and probably ~0.5-1 Myr). The new age estimate is in agreement with the evolutionary status of the various protostars in the region and with its compactness (<1.3 pc across), but results in a conflict with the low disk and accretion fraction (only 50-65% of low-mass stars appear to have protoplanetary disks, and most transitional and homologously depleted disks are consistent with no accretion) and with the evolutionary features observed in the mid-IR spectra and spectral energy distributions of the disks.

Context. The formation and early evolution of high- and intermediate-mass stars towards the main sequence involves the interplay of stars in a clustered and highly complex environment. To obtain a full census of this interaction, the Formation and Early evolution of Massive Stars (FEMS) collaboration studies a well-selected sample of 10 high-mass star-forming regions. Aims. In this study we examine the stellar content of the high-mass star-forming region centered on IRAS 06084-0611 in the Monoceros R2 cloud. Methods. Using the near-infrared H- and K-band spectra from the VLT/SINFONI instrument on the ESO Very Large Telescope (VLT)and photometric near-infrared NTT/SOFI, 2MASS and Spitzer/IRAC data, we were able to determine the spectral types for the most luminous stars in the cluster. Results. Two very young and reddened massive stars have been detected by SINFONI: a massive Young Stellar Object (YSO) con- sistent with an early-B spectral type and a Herbig Be star. Furthermore, stars of spectral type G and K are detected while still in the Pre-Main Sequence (PMS) phase. We derive additional properties such as temperatures, extinctions, radii and masses. We present a Hertzsprung-Russell diagram and find most objects having intermediate masses between \sim1.5-2.5 M\odot. For these stars we derive a median cluster age of \sim4 Myr. Conclusions. Using Spitzer/IRAC data we confirm earlier studies that the younger class 0/I objects are centrally located while the class II objects are spread out over a larger area, with rough scale size radii of \sim0.5 pc and \sim1.25 pc respectively. Moreover, the presence of a massive YSO, an ultracompact H ii region and highly reddened objects in the center of the cluster suggest a much younger age of < 1 Myr. A possible scenario for this observation would be sequential star formation along the line of sight; from a cluster of intermediate-mass to high-mass stars.

We present an analysis of the large-scale molecular cloud structure and of the stability of clumpy structures in nearby molecular clouds. In our recent work, we identified a structural transition in molecular clouds by studying the probability distributions of gas column densities in them. In this paper, we further examine the nature of this transition. The transition takes place at the visual extinction of A_V^tail = 2-4 mag, or equivalently, at \Sigma^tail = 40-80 Ms pc^{-2}. The clumps identified above this limit have wide ranges of masses and sizes, but a remarkably constant mean volume density of n = 10^3 cm^{-3}. This is 5-10 times larger than the density of the medium surrounding the clumps. By examining the stability of the clumps, we show that they are gravitationally unbound entities, and that the external pressure from the parental molecular cloud is a significant source of confining pressure for them. Then, the structural transition at A_V^tail may be linked to a transition between this population and the surrounding medium. The star formation rates in the clouds correlate strongly with the total mass in the clumps, i.e, with the mass above A_V^tail, dropping abruptly below that threshold. These results imply that the formation of pressure confined clumps introduces a prerequisite for star formation. Furthermore, they give a physically motivated explanation for the recently reported relation between the star formation rates and the amount of dense material in molecular clouds. Likewise, they give rise to a natural threshold for star formation at A_V^tail.

We present a study of dense structures in the L1495 filament in the Taurus Molecular Cloud and examine its star-forming properties. In particular we construct a dust extinction map of the filament using deep near-infrared observations, exposing its small-scale structure in unprecedented detail. The filament shows highly fragmented substructures and a high mass-per-length value of M_line=17 Msun/pc, reflecting star-forming potential in all parts of it. However, a part of the filament, namely B211, is remarkably devoid of young stellar objects. We argue that in this region the initial filament collapse and fragmentation is still taking place and star formation is yet to occur. In the star-forming part of the filament, we identify 39 cores with masses from 0.4...10 Msun and preferred separations in agreement with the local Jeans length. Most of these cores exceed the Bonnor-Ebert critical mass, and are therefore likely to collapse and form stars. The Dense Core Mass Function follows a power law with exponent Gamma=1.2, a form commonly observed in star-forming regions.

We present Herschel observations of the isolated, low-mass star-forming Bok globule CB244. It contains two cold sources, a low-mass Class 0 protostar and a starless core, which is likely to be prestellar in nature, separated by 90 arcsec (~ 18000 AU). The Herschel data sample the peak of the Planck spectrum for these sources, and are therefore ideal for dust-temperature and column density modeling. With these data and a near-IR extinction map, the MIPS 70 micron mosaic, the SCUBA 850 micron map, and the IRAM 1.3 mm map, we model the dust-temperature and column density of CB244 and present the first measured dust-temperature map of an entire star-forming molecular cloud. We find that the column-averaged dust-temperature near the protostar is ~ 17.7 K, while for the starless core it is ~ 10.6K, and that the effect of external heating causes the cloud dust-temperature to rise to ~ 17 K where the hydrogen column density drops below 10^21 cm^-2. The total hydrogen mass of CB244 (assuming a distance of 200 pc) is 15 +/- 5 M_sun. The mass of the protostellar core is 1.6 +/- 0.1 M_sun and the mass of the starless core is 5 +/- 2 M_sun, indicating that ~ 45% of the mass in the globule is participating in the star-formation process.

Aims: We derive the probability density functions (PDFs) of column density for a complete sample of prominent molecular cloud complexes closer than 200 pc. Methods: We derive near-infrared dust extinction maps for 23 molecular cloud complexes, using the "nicest" colour excess mapping technique and data from the 2MASS archive. The extinction maps are then used to examine the column density PDFs in the clouds. Results: The column density PDFs of most molecular clouds are well-fitted by log-normal functions at low column densities (0.5 mag < A_v < 3-5 mag). However, at higher column densities prominent, power-law-like wings are common. In particular, we identify a trend among the PDFs: active star-forming clouds always have prominent non-log-normal wings. In contrast, clouds without active star formation resemble log-normals over the whole observed column density range, or show only low excess of higher column densities. This trend is also reflected in the cumulative PDFs, showing that the fraction of high column density material is significantly larger in star-forming clouds. These observations are in agreement with an evolutionary trend where turbulent motions are the main cloud-shaping mechanism for quiescent clouds, but the density enhancements induced by them quickly become dominated by gravity (and other mechanisms) which is strongly reflected by the shape of the column density PDFs. The dominant role of the turbulence is restricted to the very early stages of molecular cloud evolution, comparable to the onset of active star formation in the clouds.

We reveal the stellar light emerging from the kiloparsec-scale, ring-like structure of the NGC5128 (Centaurus A) galaxy in unprecedented detail. We use arcsecond-scale resolution near infrared images to create a "dust-free" view of the central region of the galaxy, which we then use to quantify the shape of the revealed structure. At the resolution of the data, the structure contains several hundreds of discreet, point-like or slightly elongated sources. Typical extinction corrected surface brightness of the structure is K_S = 16.5 mag/arcsec^2, and we estimate the total near infrared luminosity of the structure to be M = -21 mag. We use diffraction limited (FWHM resolution of ~ 0.1", or 1.6 pc) near infrared data taken with the NACO instrument on VLT to show that the structure decomposes into thousands of separate, mostly point-like sources. According to the tentative photometry, the most luminous sources have M_K = -12 mag, naming them red supergiants or relatively low-mass star clusters. We also discuss the large-scale geometry implied by the reddening signatures of dust in our near infrared images.

We reveal the stellar light emerging from the kiloparsec-scale, ring-like structure of the NGC5128 (Centaurus A) galaxy in unprecedented detail. We use arcsecond-scale resolution near infrared images to create a "dust-free" view of the central region of the galaxy, which we then use to quantify the shape of the revealed structure. At the resolution of the data, the structure contains several hundreds of discreet, point-like or slightly elongated sources. Typical extinction corrected surface brightness of the structure is K_S = 16.5 mag/arcsec^2, and we estimate the total near infrared luminosity of the structure to be M = -21 mag. We use diffraction limited (FWHM resolution of ~ 0.1", or 1.6 pc) near infrared data taken with the NACO instrument on VLT to show that the structure decomposes into thousands of separate, mostly point-like sources. According to the tentative photometry, the most luminous sources have M_K = -12 mag, naming them red supergiants or relatively low-mass star clusters. We also discuss the large-scale geometry implied by the reddening signatures of dust in our near infrared images.

In this paper we derive an improved core mass function (CMF) for the Pipe Nebula from a detailed comparison between measurements of visual extinction and molecular-line emission. We have compiled a refined sample of 201 dense cores toward the Pipe Nebula using a 2-dimensional threshold identification algorithm informed by recent simulations of dense core populations. Measurements of radial velocities using complimentary C18O (1-0) observations enable us to cull out from this sample those 43 extinction peaks that are either not associated with dense gas or are not physically associated with the Pipe Nebula. Moreover, we use the derived C18O, central velocities to differentiate between single cores with internal structure and blends of two or more physically distinct cores, superposed along the same line-of-sight. We then are able to produce a more robust dense core sample for future follow-up studies and a more reliable CMF than was possible previously. We confirm earlier indications that the CMF for the Pipe Nebula departs from a single power-law like form with a break or knee at M ~ 2.7 +/- 1.3 Msun. Moreover, we also confirm that the CMF exhibits a similar shape to the stellar IMF, but is scaled to higher masses by a factor of ~4.5. We interpret this difference in scaling to be a measure of the star formation efficiency (22 +/- 8%). This supports earlier suggestions that the stellar IMF may originate more or less directly from the CMF.

The aims of this study are to determine the properties and spatial distribution of dense cores in Ori B9, and to estimate their ages and dynamical timescales. The cloud was mapped in the 870 micron continuum with APEX/LABOCA, and selected positions were observed in the lines of N2H+ and N2D+ using IRAM-30m. These were used together with our previous H2D+ observations. Moreover, archival FIR Spitzer/MIPS maps were combined with the LABOCA map to distinguish between pre- and protostellar cores, and to estimate the evolutionary stages of protostars. Twelve dense cores were detected at 870 micron in the Ori B9 cloud. The submm cores constitute ~4% of the total mass of the Ori B9 region. There is an equal number of pre- and protostellar cores. Two of the submm sources, which we call SMM 3 and SMM 4, are previously unknown Class 0 candidates. We found a moderate degree of deuteration in N2H+ (0.03-0.04). There is, furthermore, evidence for N2H+ depletion in the core SMM 4. We derive a relatively high degree of ionization (~10^-7) in the clump associated with IRAS 05405-0117. The ambipolar diffusion timescales for two of the cores are ~70-100 times longer than the free-fall time. The distribution and masses of dense cores in Ori B9 are similar to those observed in more active regions in Orion, where the statistical core properties have been explained by turbulent fragmentation. The 50/50 proportions of pre- and protostellar cores suggest that duration of the prestellar phase is comparable to the free-fall time. This timescale can be questioned, however, on the basis of chemical data on the IRAS 05405-0117 region. A possible explanation is that this survey samples only the densest, i.e., dynamically most advanced cores.

Aims: We examine the recoverability and completeness limits of the dense core mass functions (CMFs) derived for a molecular cloud using extinction data and a core identification scheme based on two-dimensional thresholding. Methods: We performed simulations where a population of artificial cores was embedded into the variable background extinction field of the Pipe nebula. We extracted the cores from the simulated extinction maps, constructed the CMFs, and compared them to the input CMFs. The simulations were repeated using a variety of extraction parameters and several core populations with differing input mass functions and differing degrees of crowding. Results: The fidelity of the observed CMF depends on the parameters selected for the core extraction algorithm for our background. More importantly, it depends on how crowded the core population is. We find that the observed CMF recovers the true CMF reliably when the mean separation of cores is larger than their mean diameter (f>1). If this condition holds, the derived CMF is accurate and complete above M > 0.8-1.5 Msun, depending on the parameters used for the core extraction. In the simulations, the best fidelity was achieved with the detection threshold of 1 or 2 times the rms-noise of the extinction data, and with the contour level spacings of 3 times the rms-noise. Choosing larger threshold and wider level spacings increases the limiting mass. The simulations show that when f>1.5, the masses of individual cores are recovered with a typical uncertainty of 25-30 %. When f=1 the uncertainty is ~60 %. In very crowded cases where f<1 the core identification algorithm is unable to recover the masses of the cores adequately. For the cores of the Pipe nebula f~2.0 and therefore the use of the method in that region is justified.

The present work aims at performing a comprehensive census and characterisation of the pre-main sequence (PMS) population in the cometary cloud L1615/L1616, in order to assess the significance of the triggered star formation scenario and investigate the impact of massive stars on its star formation history and mass spectrum. Our study is based on UBVRcIc and JHKs photometry, as well as optical multi-object spectroscopy. We performed a physical parametrisation of the young stellar population in L1615/L1616. We identified 25 new T Tauri stars mainly projected on the dense head of the cometary cloud, almost doubling the current number of known members. We studied the spatial distribution of the cloud members as a function of the age and H$\alpha$ emission. The star formation efficiency in the cloud is about 7-8 %, as expected for molecular clouds in the vicinity of OB associations. The slope of the initial mass function (IMF), in the mass range 0.1<M<5.5 $M_{\odot}$, is consistent with that of other T and OB associations, providing further support of an universal IMF down to the hydrogen burning limit, regardless of environmental conditions. The cometary appearance, as well as the high star formation efficiency, can be explained in terms of triggered star formation induced by the strong UV radiation from OB stars or supernovae shockwaves. The age spread as well as both the spatial and age distribution of the PMS objects provide strong evidence of sequential, multiple events and possibly still ongoing star formation activity in the cloud.

[Abridged] We describe the near-infrared reddening signature of giant molecular clouds (GMCs) in external galaxies. In particular, we examine the E(J-H) and E(H-K) color-excesses, and the effective extinction law observed in discrete GMC regions. We also study the effect of the relative scale height of the GMC distribution to the color-excesses, and to the observed mass function of GMCs. We perform Monte Carlo radiative transfer simulations with 3D models of stellar radiation and clumpy dust distributions, resembling a face-on geometry. The scattered light is included in the models, and near-infrared color maps are calculated from the simulated data. The effective near-infrared reddening law, i.e. the ratio E(J-H)/E(H-K), has a value close to unity in GMC regions. The ratio depends on the relative scale height of GMCs, xi, and for xi values 0.1...0.75 we find the typical ratios of 0.6...1.1. The effective extinction law turns out to be very flat in GMC regions. We find the ratios of apparent extinctions of A(H)/A(K)=1.35...1.55 and A(J)/A(H)=1.15. The effect of the scattered flux on the effective reddening law, as well as on the effective extinction law, is significant. Regarding the GMC mass function, we find no correlation between the input and observed slopes of the mass functions. Rather, the observed slope reflects the parameter $\xi$ and the dynamical range of the mass function. We estimate that only a fraction of 10...20 % of the total mass of GMCs is recovered, if the observed color-excess values are transformed to masses using the Galactic reddening law. In the case of individual clouds the fraction can vary between ~0...50 %.

[Abridged] The giant molecular clouds (GMCs) of external galaxies can be mapped with sub-arcsecond resolution using multiband observations in the near-infrared. However, the interpretation of the observed features, and their transformation into physical quantities, is greatly hampered by the effects arising from the unknown geometry and the scattering of light by dust particles. We examine the relation between the observed NIR reddening and the column density of the dust clouds. In this paper we particularly assess the feasibility of deriving the mass function of GMCs from NIR color-excess data. We perform Monte Carlo radiative transfer simulations with 3D models of stellar radiation and clumpy dust distributions. We include the scattered light in the models and calculate NIR color maps from the simulated data. We extract clumps from the color maps and compare the observed mass function to the true mass function. For the physical configuration chosen in this study, essentially a face-on geometry, the observed mass function is a non-trivial function of the true mass function. The dynamical range of the observed mass function is confined to ~10^3.5... 10^5.5 M_\odot regardless of the dynamical range of the true mass function. The color maps are more sensitive in detecting the high-mass end of the mass function, and on average the masses of clouds are underestimated by a factor of ~10. The simulations show that the cloud mass function derived from JHK color excess data using simple foreground screening geometry cannot be regarded as a one-to-one tracer of the underlying mass function.

We present a study of radial density structure of the star forming globule, DCld303.8-14.2 (DC303), and a non-star forming globule, Thumbprint Nebula (TPN), using near-infrared data taken with the ISAAC instrument on the Very Large Telescope. We derive the extinction through the globules using the color excess technique and examine the radial density distribution using Bonnor-Ebert and power-law models. The two globules have significantly different density structures. The extinction profile of DC303 is best fitted with a single power-law with an exponent p = 2.29 \pm 0.08. An unstable Bonnor-Ebert model with a dimensionless parameter \xi_max = 23 \pm 3 provides equally good fit to data. The extinction profile of TPN flattens at small radii, making the profile significantly different from the profile of DC303. We are unable to fit the Bonnor-Ebert model for TPN in a robust manner, but derive the lower limit \xi_max > 8 for the dimensionless outer edge. The density profile derived for TPN is typical compared to recently observed pre-protostellar globules, with high \xi_max value which could be interpreted as the presence of significant additional support or very slow contraction.

We investigate the relationship between the C18O column density and the visual extinction in Chamaeleon I and in a part of the Chamaeleon III molecular cloud. The C18O column densities, N(C18O), are calculated from J=1-0 rotational line data observed with the SEST telescope. The visual extinctions, Av, are derived using JHK photometry from the 2MASS survey and the NICER color excess technique. In contrast with the previous results of Hayakawa et al. (2001), we find that the average N(C18O) / Av ratios are similar in Cha I and Cha III, and lie close to values derived for other clouds, i.e. N(C18O) ~ 2 x 10^14 cm^-2 (Av - 2). We find, however, clear deviations from this average relationship towards individual clumps. Larger than average N(C18O) / Av ratios can be found in clumps associated with the active star forming region in the northern part of Ch a I. On the other hand, some regions in the relatively quiescent southern part of Cha I show smaller than average N(C18O) / Av ratios and also very shallow proportionality between N(C18O) and Av. The shallow proportionality suggests that C18O is heavily depleted in these regions. As the degree of depletion is proportional to the gas density, these regions probably contain very dense, cold cores, which do not stand out in CO mappings. A comparison with the dust temperature map derived from the ISO data shows that the most prominent of the potentially depleted cores indeed coincides with a dust temperature minimum. It seems therefore feasible to use N(C18O) and Av data together for identifying old, dense cores in large scale mappings.