CO plays an important role in interstellar molecular clouds, both as a coolant, and as a diagnostic molecule. However, a proper evaluation of the cooling rate due to CO requires a determination of the populations of many levels, the spontaneous and stimulated radiative de-excitation rates between these levels, and the transfer of the emitted multi-line radiation; additionally, this must be done for three isotopologues. It would be useful to have a simple analytic formulation that avoided these complications; this could then be used in situations where CO plays an important role as a coolant, but the details of this role are not the main concern. We derive such a formulation here, by first considering the two asymptotic forms that obtain in the limits of (a) low volume-density and optical depth, and (b) high volume-density and optical depth. These forms are then combined in such a way as to fit the detailed numerical results from Goldsmith \& Langer (1978), which cover low temperatures, and a range of physical conditions where the interplay of thermal and sub-thermal excitation, optical-depth effects, and the contributions from rare isotopologues, are all important. The fit is obtained using the Metropolis-Hastings method, and reproduces the results of GL78 well. It is a purely local and analytic function of state --- specifically a function of the density, $\rho$, isothermal sound speed, $a$, CO abundance, $X_{_{\rm CO}}$, and velocity divergence, $\nabla\cdot{\boldsymbol\upsilon}$. As an application, we consider the cooling layer following a slow steady non-magnetic planar J-shock. We show that, if the post-shock cooling is dominated by CO and its isotopologues, the thickness of the post-shock cooling layer is very small and approximately independent of the pre-shock velocity, $\upsilon_o$, or pre-shock isothermal sound speed, $a_o$.

The multifractal geometry remains an under-exploited approach to describe and quantify the large-scale structure of interstellar clouds. In this paper, the typical tools of multifractal analysis are applied to Herschel far-infrared (70-500 $\mu$m) dust continuum maps, which represent an ideal case of study. This dust component is a relatively optically thin tracer at these wavelengths and the size in pixel of the maps is well suitable for this statistical analysis. We investigate the so-called multifractal spectrum and generalised fractal dimension of six Hi-GAL regions in the third Galactic quadrant. We measure systematic variations of the spectrum at increasing wavelength, which generally correspond to an increasing complexity of the image, and we observe peculiar behaviours of the investigated fields, strictly related to the presence of high-emission regions, which in turn are connected to star formation activity. The same analysis is applied to synthetic column density maps, generated from numerical turbulent molecular cloud models and from fractal Brownian motion (fBm), allowing for the confrontation of the observations with models with well controlled physical parameters. This comparison shows that, in terms of multifractal descriptors of the structure, fBm images exhibit a very different, quite unrealistic behaviour when compared with Hi-GAL observations, whereas the numerical simulations appear in some cases (depending on the specific model) more similar to the observations. Finally, the link between mono-fractal parameters (commonly found in the literature) and multifractal indicators is investigated: the former appear to be only partially connected with the latter, with the multifractal analysis offering a more extended and complete characterization of the cloud structure.

Molecular-line observations of filaments in star-forming regions have revealed the existence of elongated coherent features within the filaments; these features are termed fibres. Here we caution that, since fibres are traced in PPV space, there is no guarantee that they represent coherent features in PPP space. We illustrate this contention using simulations of the growth of a filament from a turbulent medium. Synthetic C$^{18}$O observations of the simulated filaments reveal the existence of fibres very similar to the observed ones, i.e. elongated coherent features in the resulting PPV data-cubes. Analysis of the PPP data-cubes (i.e. 3D density fields) also reveals elongated coherent features, which we term sub-filaments. Unfortunately there is very poor correspondence between the fibres and the sub-filaments in the simulations. Both fibres and sub-filaments derive from inhomogeneities in the turbulent accretion flow onto the main filament. As a consequence, fibres are often affected by line-of-sight confusion. Similarly, sub-filaments are often affected by large velocity gradients, and even velocity discontinuities. These results suggest that extreme care should be taken when using velocity coherent features to constrain the underlying substructure within a filament.

We present simulations of the dynamical evolution of young embedded star clusters. Our initial conditions are directly derived from X-ray, infrared, and radio observations of local systems, and our models evolve both gas and stars simultaneously. Our regions begin with both clustered and extended distributions of stars, and a gas distribution which can include a filamentary structure in addition to gas surrounding the stellar subclusters. We find that the regions become spherical, monolithic, and smooth quite quickly, and that the dynamical evolution is dominated by the gravitational interactions between the stars. In the absence of stellar feedback, the gas moves gently out of the centre of our regions but does not have a significant impact on the motions of the stars at the earliest stages of cluster formation. Our models at later times are consistent with observations of similar regions in the local neighbourhood. We conclude that the evolution of young proto-star clusters is relatively insensitive to reasonable choices of initial conditions. Models with more realism, such as an initial population of binary and multiple stars and ongoing star formation, are the next step needed to confirm these findings.

Large scale surveys have brought about a revolution in astronomy. To analyse the resulting wealth of data, we need automated tools to identify, classify, and quantify the important underlying structures. We present here a method for classifying and quantifying a pixelated structure, based on its principal moments of inertia. The method enables us to automatically detect, and objectively compare, centrally condensed cores, elongated filaments and hollow rings. We illustrate the method by applying it to (i) observations of surface-density from Hi-GAL, and (ii) simulations of filament growth in a turbulent medium. We limit the discussion here to 2D data; in a future paper we will extend the method to 3D data.

Hi-GAL is a large-scale survey of the Galactic plane, performed with Herschel in five infrared continuum bands between 70 and 500 $\mu$m. We present a band-merged catalogue of spatially matched sources and their properties derived from fits to the spectral energy distributions (SEDs) and heliocentric distances, based on the photometric catalogs presented in Molinari et al. (2016a), covering the portion of Galactic plane $-71.0^{\circ}< \ell < 67.0^{\circ}$. The band-merged catalogue contains 100922 sources with a regular SED, 24584 of which show a 70 $\mu$m counterpart and are thus considered proto-stellar, while the remainder are considered starless. Thanks to this huge number of sources, we are able to carry out a preliminary analysis of early stages of star formation, identifying the conditions that characterise different evolutionary phases on a statistically significant basis. We calculate surface densities to investigate the gravitational stability of clumps and their potential to form massive stars. We also explore evolutionary status metrics such as the dust temperature, luminosity and bolometric temperature, finding that these are higher in proto-stellar sources compared to pre-stellar ones. The surface density of sources follows an increasing trend as they evolve from pre-stellar to proto-stellar, but then it is found to decrease again in the majority of the most evolved clumps. Finally, we study the physical parameters of sources with respect to Galactic longitude and the association with spiral arms, finding only minor or no differences between the average evolutionary status of sources in the fourth and first Galactic quadrants, or between "on-arm" and "inter-arm" positions.

Many young star clusters appear to be fractal, i.e. they appear to be concentrated in a nested hierarchy of clusters within clusters. We present a new algorithm for statistically analysing the distribution of stars to quantify the level of sub-structure. We suggest that, even at the simplest level, the internal structure of a fractal cluster requires the specification of three parameters. (i) The 3D fractal dimension, $\mathcal{D}$, measures the extent to which the clusters on one level of the nested hierarchy fill the volume of their parent cluster. (ii) The number of levels, $\mathcal{L}$, reflects the finite ratio between the linear size of the large root-cluster at the top of the hierarchy, and the smallest leaf-clusters at the bottom of the hierarchy. (iii) The volume-density scaling exponent, $\mathcal{C}=-\textrm{d}\ln[\delta n]/\textrm{d}\ln[L]$ measures the factor by which the excess density, $\delta n$, in a structure of scale $L$, exceeds that of the background formed by larger structures; it is similar, but not exactly equivalent, to the exponent in Larson's scaling relation between density and size for molecular clouds. We describe an algorithm which can be used to constrain the values of $({\cal D},{\cal L},{\cal C})$ and apply this method to artificial and observed clusters. We show that this algorithm is able to reliably describe the three dimensional structure of an artificial star cluster from the two dimensional projection, and quantify the varied structures observed in real and simulated clusters.

We present the results of SPH simulations in which two clouds, each having mass $M_{_{\rm{o}}}\!=\!500\,{\rm M}_{_\odot}$ and radius $R_{_{\rm{o}}}\!=\!2\,{\rm pc}$, collide head-on at relative velocities of $\Delta v_{_{\rm{o}}} =2.4,\;2.8,\;3.2,\;3.6\;{\rm and}\;4.0\,{\rm km}\,{\rm s}^{-1}$. There is a clear trend with increasing $\Delta v_{_{\rm{o}}}$. At low $\Delta v_{_{\rm{o}}}$, star formation starts later, and the shock-compressed layer breaks up into an array of predominantly radial filaments; stars condense out of these filaments and fall, together with residual gas, towards the centre of the layer, to form a single large-$N$ cluster, which then evolves by competitive accretion, producing one or two very massive protostars and a diaspora of ejected (mainly low-mass) protostars; the pattern of filaments is reminiscent of the hub and spokes systems identified recently by observers. At high $\Delta v_{_{\rm{o}}}$, star formation occurs sooner and the shock-compressed layer breaks up into a network of filaments; the pattern of filaments here is more like a spider's web, with several small-$N$ clusters forming independently of one another, in cores at the intersections of filaments, and since each core only spawns a small number of protostars, there are fewer ejections of protostars. As the relative velocity is increased, the {\it mean} protostellar mass increases, but the {\it maximum} protostellar mass and the width of the mass function both decrease. We use a Minimal Spanning Tree to analyse the spatial distributions of protostars formed at different relative velocities.