In attempting to quantify statistically the density structure of the interstellar medium, astronomers have considered a variety of fractal models. Here we argue that, to properly characterise a fractal model, one needs to define precisely the algorithm used to generate the density field, and to specify -- at least -- three parameters: one parameter constrains the spatial structure of the field; one parameter constrains the density contrast between structures on different scales; and one parameter constrains the dynamic range of spatial scales over which self-similarity is expected (either due to physical considerations, or due to the limitations of the observational or numerical technique generating the input data). A realistic fractal field must also be noisy and non-periodic. We illustrate this with the exponentiated fractional Brownian motion (xfBm) algorithm, which is popular because it delivers an approximately lognormal density field, and for which the three parameters are, respectively, the power spectrum exponent, $\beta$, the exponentiating factor, ${\cal S}$, and the dynamic range, ${\cal R}$. We then explore and compare two approaches that might be used to estimate these parameters: Machine Learning and the established $\Delta$-Variance procedure. We show that for $2\leq\beta \leq 4$ and $0\leq{\cal S}\leq 3$, a suitably trained Convolutional Neural Network is able to estimate objectively both $\beta$ (with root-mean-square error $\epsilon_{_\beta}\sim 0.12$) and ${\cal S}$ (with $\epsilon_{_{\cal S}}\sim 0.29$). $\;\Delta$-variance is also able to estimate $\beta$, albeit with a somewhat larger error ($\epsilon_{_\beta}\sim 0.17$) and with some human intervention, but is not able to estimate ${\cal S}$.

We have analysed Herschel observations of M31, using the PPMAP procedure. The resolution of PPMAP images is sufficient (31 pc on M31) that we can analyse far-IR dust emission on the scale of Giant Molecular Clouds. By comparing PPMAP estimates of the far-IR emission optical depth at 300 microns (tau_300), and the near-IR extinction optical depth at 1.1 microns (tau_1.1) obtained from the reddening of RGB stars, we show that the ratio R_OBS.tau = tau_1.1/tau_300 falls in the range 500 to 1500. Such low values are incompatible with many commonly used theoretical dust models, which predict values of R_MODEL.kappa = kappa_1.1/kappa_300 (where kappa is the dust opacity coefficient) in the range 2500 to 4000. That is, unless a large fraction, at least 60%, of the dust emitting at 300 microns is in such compact sources that they are unlikely to intercept the lines of sight to a distributed population like RGB stars. This is not a new result: variants obtained using different observations and/or different wavelengths have already been reported by other studies. We present two analytic arguments for why it is unlikely that at least 60% of the emitting dust is in sufficiently compact sources. Therefore it may be necessary to explore the possibility that the discrepancy between observed values of R_OBS.tau and theoretical values of R_MODEL.kappa is due to limitations in existing dust models. PPMAP also allows us to derive optical-depth weighted mean values for the emissivity index, beta = - dln(kappa_lambda)/dln(lambda), and the dust temperature, T, denoted betabar and Tbar. We show that, in M31, R_OBS.tau is anti-correlated with betabar according to R_OBS.tau = 2042(+/-24)-557(+/-10)betabar. If confirmed, this provides a challenging constraint on the nature of interstellar dust in M31.

We have analysed the Herschel and SCUBA-2 dust continuum observations of the main filament in the Taurus L1495 star forming region, using the Bayesian fitting procedure PPMAP. (i) If we construct an average profile along the whole length of the filament, it has fwhm $\simeq 0.087\pm 0.003\,{\rm pc};\;$, but the closeness to previous estimates is coincidental. (ii) If we analyse small local sections of the filament, the column-density profile approximates well to the form predicted for hydrostatic equilibrium of an isothermal cylinder. (iii) The ability of PPMAP to distinguish dust emitting at different temperatures, and thereby to discriminate between the warm outer layers of the filament and the cold inner layers near the spine, leads to a significant reduction in the surface-density, $\varSigma$, and hence in the line-density, $\mu$. If we adopt the canonical value for the critical line-density at a gas-kinetic temperature of $10\,{\rm K}$, $\mu_{_{\rm CRIT}}\simeq 16\,{\rm M_{_\odot}\,pc^{-1}}$, the filament is on average trans-critical, with ${\bar\mu}\sim \mu_{_{\rm CRIT}};\;$ local sections where $\mu >\mu_{_{\rm CRIT}}$ tend to lie close to pre-stellar cores. (iv) The ability of PPMAP to distinguish different types of dust, i.e. dust characterised by different values of the emissivity index, $\beta$, reveals that the dust in the filament has a lower emissivity index, $\beta\leq1.5$, than the dust outside the filament, $\beta\geq 1.7$, implying that the physical conditions in the filament have effected a change in the properties of the dust.

The degree of fractal substructure in molecular clouds can be quantified by comparing them with Fractional Brownian Motion (FBM) surfaces or volumes. These fields are self-similar over all length scales and characterised by a drift exponent $H$, which describes the structural roughness. Given that the structure of molecular clouds and the initial structure of star clusters are almost certainly linked, it would be advantageous to also apply this analysis to clusters. Currently, the structure of star clusters is often quantified by applying $\mathcal{Q}$ analysis. $\mathcal{Q}$ values from observed targets are interpreted by comparing them with those from artificial clusters. These are typically generated using a Box-Fractal (BF) or Radial Density Profile (RDP) model. We present a single cluster model, based on FBM, as an alternative to these models. Here, the structure is parameterised by $H$, and the standard deviation of the log-surface/volume density $\sigma$. The FBM model is able to reproduce both centrally concentrated and substructured clusters, and is able to provide a much better match to observations than the BF model. We show that $\mathcal{Q}$ analysis is unable to estimate FBM parameters. Therefore, we develop and train a machine learning algorithm which can estimate values of $H$ and $\sigma$, with uncertainties. This provides us with a powerful method for quantifying the structure of star clusters in terms which relate to the structure of molecular clouds. We use the algorithm to estimate the $H$ and $\sigma$ for several young star clusters, some of which have no measurable BF or RDP analogue.

We suggest that Bipolar H{\sc ii} Regions may be the aftermath of collisions between clouds. Such a collision will produce a shock-compressed layer, and a star cluster can then condense out of the dense gas near the centre of the layer. If the clouds are sufficiently massive, the star cluster is likely to contain at least one massive star, which emits ionising radiation, and excites an H{\sc ii} region, which then expands, sweeping up the surrounding neutral gas. Once most of the matter in the clouds has accreted onto the layer, expansion of the H{\sc ii} Region meets little resistance in directions perpendicular to the mid-plane of the layer, and so it expands rapidly to produce two lobes of ionised gas, one on each side of the layer. Conversely, in directions parallel to the mid-plane of the layer, expansion of the H{\sc ii} Region stalls due to the ram-pressure of the gas that continues to fall towards the star cluster from the outer parts of the layer; a ring of dense neutral gas builds up around the waist of the Bipolar H{\sc ii} Region, and may spawn a second generation of star formation. In this paper we present a dimensionless model for the flow of ionised gas in a Bipolar H{\sc ii} Region created according to the above scenario, and predict the characteristics of the resulting freefree continuum and recombination-line emission. This dimensionless model can be scaled to the physical parameters of any particular system. Our intention is that these predictions will be useful in testing the scenario outlined above, and thereby providing indirect support for the role of cloud/cloud collisions in triggering star formation.

We present new Herschel-derived images of warm dust in the Andromeda Galaxy, M31, with unprecedented spatial resolution (~ 30 pc), column density accuracy, and constraints on the three-dimensional distributions of dust temperature and dust opacity index (hence grain size and composition), based on the new PPMAP Bayesian analysis procedure. We confirm the overall radial variation of dust opacity index reported by other recent studies, including the central decrease within ~ 3 kpc of the nucleus. We also investigate the detailed distribution of dust in the nuclear region, a prominent feature of which is a ~ 500 pc bar-like structure seen previously in H{\alpha}. The nature of this feature has been the subject of some debate. Our maps show it to be the site of the warmest dust, with a mean line-of-sight temperature ~ 30 K. A comparison with the stellar distribution, based on 2MASS data, provides strong evidence that it is a gravitationally induced bar. A comparison with radial velocity maps suggests the presence of an inflow towards the nucleus from opposite directions along this bar, fed presumably by the nuclear spiral with which it appears to connect. Such behaviour is common in large-scale bars in spiral galaxies, as is the phenomenon of nested bars whereby a subkiloparsec nuclear bar exists within a large-scale primary bar. We suggest that M31 represents an example of such nesting.

Observations of protostars are often compared with synthetic observations of models in order to infer the underlying physical properties of the protostars. The majority of these models have a single protostar, attended by a disc and an envelope. However, observational and numerical evidence suggests that a large fraction of protostars form as multiple systems. This means that fitting models of single protostars to observations may be inappropriate. We produce synthetic observations of protostellar multiple systems undergoing realistic, non-continuous accretion. These systems consist of multiple protostars with episodic luminosities, embedded self-consistently in discs and envelopes. We model the gas dynamics of these systems using smoothed particle hydrodynamics and we generate synthetic observations by post-processing the snapshots using the \textsc{spamcart} Monte Carlo radiative transfer code. We present simulation results of three model protostellar multiple systems. For each of these, we generate $4\times10^4$ synthetic spectra at different points in time and from different viewing angles. We propose a Bayesian method, using similar calculations to those presented here, but in greater numbers, to infer the physical properties of protostellar multiple systems from observations.

We describe new Hi-GAL based maps of the entire Galactic Plane, obtained using continuum data in the wavelength range 70-500 $\mu$m. These maps are derived with the PPMAP procedure, and therefore represent a significant improvement over those obtained with standard analysis techniques. Specifically they have greatly improved resolution (12 arcsec) and, in addition to more accurate integrated column densities and mean dust temperatures, they give temperature-differential column densities, i.e., separate column density maps in twelve distinct dust temperature intervals, along with the corresponding uncertainty maps. The complete set of maps is available on-line. We briefly describe PPMAP and present some illustrative examples of the results. These include (a) multi-temperature maps of the Galactic HII region W5-E, (b) the temperature decomposition of molecular cloud column-density probability distribution functions, and (c) the global variation of mean dust temperature as a function of Galactocentric distance. Amongst our findings are: (i) a strong localised temperature gradient in W5-E in a direction orthogonal to that towards the ionising star, suggesting an alternative heating source and providing possible guidance for models of the formation of the bubble complex, and (ii) the overall radial profile of dust temperature in the Galaxy shows a monotonic decrease, broadly consistent both with models of the interstellar radiation field and with previous estimates at lower resolution. However, we also find a central temperature plateau within ~ 6 kpc of the Galactic centre, outside of which is a pronounced steepening of the radial profile. This behaviour may reflect the greater proportion of molecular (as opposed to atomic) gas in the central region of the Galaxy.

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 a code for generating synthetic SEDs and intensity maps from Smoothed Particle Hydrodynamics simulation snapshots. The code is based on the Lucy (1999) Monte Carlo Radiative Transfer method, i.e. it follows discrete luminosity packets, emitted from external and/or embedded sources, as they propagate through a density field, and then uses their trajectories to compute the radiative equilibrium temperature of the ambient dust. The density is not mapped onto a grid, and therefore the calculation is performed at exactly the same resolution as the hydrodynamics. We present two example calculations using this method. First, we demonstrate that the code strictly adheres to Kirchhoff's law of radiation. Second, we present synthetic intensity maps and spectra of an embedded protostellar multiple system. The algorithm uses data structures that are already constructed for other purposes in modern particle codes. It is therefore relatively simple to implement.

We discuss the factors influencing the formation and gravitational fragmentation of protostellar discs. We start with a review of how observations of prestellar cores can be analysed statistically to yield plausible initial conditions for simulations of their subsequent collapse. Simulations based on these initial conditions show that, despite the low levels of turbulence in prestellar cores, they deliver primary protostars and associated discs which are routinely subject to stochastic impulsive perturbations; consequently misalignment of the spins and orbits of protostars are common. Also, the simulations produce protostars that collectively have a mass function and binary statistics matching those observed in nearby star formation regions, but only if a significant fraction of the turbulent energy in the core is solenoidal, and accretion onto the primary protostar is episodic with a duty cycle > 3000 yr. Under this circumstance a core typically spawns between 4 and 5 protostars, with high efficiency, and the lower-mass protostars are mainly formed by disc fragmentation. The requirement that a proto-fragment in a disc lose thermal energy on a dynamical timescale dictates that {\bf there is a sweet spot for} disc fragmentation at radii 70 AU < R < 100 AU and temperatures 10 K < T < 20 K, and this might explain the Brown Dwarf Desert.

We use Smoothed Particle Hydrodynamics to explore the circumstances under which an isolated very-low-mass prestellar core can be formed by colliding turbulent flows and collapse to form a brown-dwarf. Our simulations suggest that the flows need not be very fast, but do need to be very strongly convergent, i.e. the gas must flow in at comparable speeds from all sides, which seems rather unlikely. We therefore revisit the object Oph-B11, which Andre, Ward-Thompson and Greaves (2012) have identified as a prestellar core with mass between $\sim 0.020\,\mathrm{M_\odot}$ and $\sim 0.030\,\mathrm{M_\odot}$. We reanalyse the observations using a Markov-chain Monte Carlo method that allows us (i) to include the uncertainties on the distance, temperature and dust mass opacity, and (ii) to consider different Bayesian prior distributions of the mass. We estimate that the posterior probability that Oph-B11 has a mass below the hydrogen burning limit at $\sim 0.075\,\mathrm{M_\odot}$, is between 0.66 and 0.86. We conclude that, if Oph-B11 is destined to collapse, it probably will form a brown dwarf. However, the flows required to trigger this appear to be so contrived that it is difficult to envisage this being the only way, or even a major way, of forming isolated brown dwarfs. Moreover, Oph-B11 could easily be a transient, bouncing, prolate core, seen end-on; there could, indeed should, be many such objects masquerading as very low-mass prestellar cores.

Disc fragmentation provides an important mechanism for producing low mass stars in prestellar cores. Here, we describe Smoothed Particle Hydrodynamics simulations which show how populations of prestellar cores evolve into stars. We find the observed masses and multiplicities of stars can be recovered under certain conditions. First, protostellar feedback from a star must be episodic. The continuous accretion of disc material on to a central protostar results in local temperatures which are too high for disc fragmentation. If, however, the accretion occurs in intense outbursts, separated by a downtime of $\sim10^4\,\mathrm{years}$, gravitational instabilities can develop and the disc can fragment. Second, a significant amount of the cores' internal kinetic energy should be in solenoidal turbulent modes. Cores with less than a third of their kinetic energy in solenoidal modes have insufficient angular momentum to form fragmenting discs. In the absence of discs, cores can fragment but results in a top heavy distribution of masses with very few low mass objects.

We present PPMAP, a Bayesian procedure that uses images of dust continuum emission at multiple wavelengths to produce resolution-enhanced image cubes of differential column-density as a function of dust temperature and position. PPMAP is based on the generic 'point process' formalism, whereby the system of interest (in this case, a dusty astrophysical structure such as a filament or prestellar core) is represented by a collection of points in a suitably defined state space. It can be applied to a variety of observational data, such as Herschel images, provided only that the image intensity is delivered by optically thin dust in thermal equilibrium. PPMAP takes full account of the instrumental point spread functions and does not require all images to be degraded to the same resolution. We present the results of testing using simulated data for a prestellar core and a fractal turbulent cloud, and demonstrate its performance with real data from the Hi-GAL survey. Specifically, we analyse observations of a large filamentary structure in the CMa OB1 giant molecular cloud. Histograms of differential column-density indicate that the warm material (T > 13 K) is distributed log-normally, consistent with turbulence, but the column-densities of the cooler material are distributed as a high density tail, consistent with the effects of self-gravity. The results illustrate the potential of PPMAP to aid in distinguishing between different physical components along the line of sight in star-forming clouds, and aid the interpretation of the associated PDFs of column density.

We present the results of an ensemble of SPH simulations that follow the evolution of prestellar cores for $0.2\,{\rm Myr}$. All the cores have the same mass, and start with the same radius, density profile, thermal and turbulent energy. Our purpose is to explore the consequences of varying the fraction of turbulent energy, $\delta_\mathrm{sol}$, that is solenoidal, as opposed to compressive; specifically we consider $\delta_\mathrm{sol}=1,\,2/3,\,1/3,\,1/9\;{\rm and}\;0$. For each value of $\delta_\mathrm{sol}$, we follow ten different realisations of the turbulent velocity field, in order also to have a measure of the stochastic variance blurring any systematic trends. With low $\delta_\mathrm{sol}(<\!1/3)$ filament fragmentation dominates and delivers relatively high mass stars. Conversely, with high values of $\delta_\mathrm{sol}(>\!1/3)$ disc fragmentation dominates and delivers relatively low mass stars. There are no discernible systematic trends in the multiplicity statistics obtained with different $\delta_\mathrm{sol}$.

It has recently been suggested that, in the field, $\sim\!\!56\%$ of Sun-like stars ($0.8\,{\rm M}_{_\odot}\lesssim M_\star\lesssim 1.2\,{\rm M}_{_\odot}$) are single. We argue here that this suggestion may be incorrect, since it appears to be based on the multiplicity frequency of systems with Sun-like primaries, and therefore takes no account of Sun-like stars that are secondary (or higher-order) components in multiple systems. When these components are included in the reckoning, it seems likely that only $\sim\!46\%$ of Sun-like stars are single. This estimate is based on a model in which the system mass function has the form proposed by Chabrier, with a power-law Salpeter extension to high masses; there is a flat distribution of mass ratios; and the probability that a system of mass $M$ is a binary is $\,0.50 + 0.46\log_{_{10}}\!\left(M/{\rm M}_{_\odot}\right)\,$ for $\,0.08\,{\rm M}_{_\odot}\leq M\leq 12.5\,{\rm M}_{_\odot}$, $\,0\,$ for $\,M<0.08\,{\rm M}_{_\odot}$, and $\,1\,$ for $\,M>12.5\,{\rm M}_{_\odot}$. The constants in this last relation are chosen so that the model also reproduces the observed variation of multiplicity frequency with primary mass. However, the more qualitative conclusion, that a minority of Sun-like stars are single, holds up for virtually all reasonable values of the model parameters. Parenthetically, it is still likely that the majority of {\it all} stars in the field are single, but that is because most M Dwarfs probably are single.

Lomax et al. have constructed an ensemble of 60 prestellar cores having masses, sizes, projected shapes, temperatures and non-thermal radial velocity dispersions that match, statistically, the cores in Ophiuchus; and have simulated the evolution of these cores using SPH. Each core has been evolved once with no radiative feedback from stars, once with continuous radiative feedback, and once with episodic radiative feedback. Here we analyse the multiplicity statistics from these simulations. With episodic radiative feedback, (i) the multiplicity frequency is ~60% higher than in the field; (ii) the multiplicity frequency and the mean semi-major axis both increase with primary mass; (iii) one third of multiple systems are hierarchical systems with more than two components; (iv) in these hierarchical systems the inner pairings typically have separations of a few au and mass ratios concentrated towards unity, whereas the outer pairings have separations of order 100 au and a flatter distribution of mass ratios. The binary statistics are compatible with observations of young embedded populations, and -- if wider orbits are disrupted preferentially by external perturbations -- with observations of mature field populations. With no radiative feedback, the results are similar to those from simulations with episodic feedback. With continuous radiative feedback, brown dwarfs are under-produced, the number of multiple systems is too low, and the statistical properties of multiple systems are at variance with observation. This suggests that star formation in Ophiuchus may only be representative of global star formation if accretion onto protostars, and hence radiative feedback, is episodic.

We have simulated star formation in prestellar cores, using SPH and initial conditions informed by observations of the cores in Ophiuchus. Because the observations are limited to two spatial dimensions plus radial velocity, we cannot infer initial conditions for the collapse of a particular core. However, with a minimum of assumptions (isotropic turbulence with a power-law spectrum, a thermal mix of compressive and solenoidal modes, a critical Bonnor-Ebert density profile) we can generate initial conditions that match, in a statistical sense, the distributions of mass, projected size and aspect ratio, thermal and non-thermal one-dimensional velocity dispersion, observed in Ophiuchus. The time between core-core collisions in Ophiuchus is sufficiently long, that we can simulate single cores evolving is isolation, and therefore we are able to resolve masses well below the opacity limit. We generate an ensemble of 100 cores, and evolve them with no radiative feedback from the stars formed, then with continuous radiative feedback, and finally with episodic radiative feedback. With no feedback the simulations produce too many brown dwarfs, and with continuous feedback too few. With episodic radiative feedback, both the peak of the protostellar mass function (at ~ 0.2 M_sun) and the ratio of H-burning stars to brown dwarfs are consistent with observations. The mass of a star is not strongly related to the mass of the core in which it forms. Low-mass cores (M ~ 0.1 M_sun) tend to collapse into single objects, whereas high-mass cores (1 > M_sun) usually fragment into several objects with a broad mass range.

Using observations of cores to infer their intrinsic properties requires the solution of several poorly constrained inverse problems. Here we address one of these problems, namely to deduce from the projected aspect ratios of the cores in Ophiuchus their intrinsic three-dimensional shapes. Four models are proposed, all based on the standard assumption that cores are randomly orientated ellipsoids, and on the further assumption that a core's shape is not correlated with its absolute size. The first and simplest model, M1, has a single free parameter, and assumes that the relative axes of a core are drawn randomly from a log-normal distribution with zero mean and standard deviation \sigma o. The second model, M2a, has two free parameters, and assumes that the log-normal distribution (with standard deviation \sigma o) has a finite mean, \mu o, defined so that \mu o<0 means elongated (prolate) cores are favoured, whereas \mu o>0 means flattened (oblate) cores are favoured. Details of the third model (M2b, two free parameters) and the fourth model (M4, four free parameters) are given in the text. Markov chain Monte Carlo sampling and Bayesian analysis are used to map out the posterior probability density functions of the model parameters, and the relative merits of the models are compared using Bayes factors. We show that M1 provides an acceptable fit to the Ophiuchus data with \sigma o ~ 0.57+/-0.06; and that, although the other models sometimes provide an improved fit, there is no strong justification for the introduction of their additional parameters.

We perform N-body simulations of young triple systems consisting of two low-mass companions orbiting around a significantly more massive primary. We find that, when the orbits of the companions are coplanar and not too widely separated, the chance of a collision between the two companions can be as high as 20 per cent. Collisions between one of the companions (always the less massive) and the primary can also occur in systems with unequal-mass companions. The chance of collisions is a few per cent in systems with more realistic initial conditions, such as with slightly non-coplanar orbits and unequal-mass companions. If the companions start widely separated then collision are very rare except in some cases when the total mass of the companions is large. We suggest that collisions between members of young multiple systems may explain some unusual young multiple systems such as apparently non-coeval companions.

The extent to which the projected distribution of stars in a cluster is due to a large-scale radial gradient, and the extent to which it is due to fractal sub-structure, can be quantified -- statistically -- using the measure ${\cal Q} = \bar{m}/\bar{s}$. Here $\bar{m}$ is the normalized mean edge length of its minimum spanning tree (i.e. the shortest network of edges connecting all stars in the cluster) and $\bar{s}$ is the correlation length (i.e. the normalized mean separation between all pairs of stars). We show how ${\cal Q}$ can be indirectly applied to grey-scale images by decomposing the image into a distribution of points from which $\bar{m}$ and $\bar{s}$ can be calculated. This provides a powerful technique for comparing the distribution of dense gas in a molecular cloud with the distribution of the stars that condense out of it. We illustrate the application of this technique by comparing ${\cal Q}$ values from simulated clouds and star clusters.