We present new observations of the active star-formation region NGC 1333 in the Perseus molecular cloud complex from the James Clerk Maxwell Telescope (JCMT) B-Fields In Star-forming Region Observations (BISTRO) survey with the POL-2 instrument. The BISTRO data cover the entire NGC 1333 complex (~ 1.5 pc x 2 pc) at 0.02 pc resolution and spatially resolve the polarized emission from individual filamentary structures for the first time. The inferred magnetic field structure is complex as a whole, with each individual filament aligned at different position angles relative to the local field orientation. We combine the BISTRO data with low- and high- resolution data derived from Planck and interferometers to study the multi-scale magnetic field structure in this region. The magnetic field morphology drastically changes below a scale of ~ 1 pc and remains continuous from the scales of filaments (~ 0.1 pc) to that of protostellar envelopes (~ 0.005 pc or ~ 1000 au). Finally, we construct simple models in which we assume that the magnetic field is always perpendicular to the long axis of the filaments. We demonstrate that the observed variation of the relative orientation between the filament axes and the magnetic field angles are well reproduced by this model, taking into account the projection effects of the magnetic field and filaments relative to the plane of the sky. These projection effects may explain the apparent complexity of the magnetic field structure observed at the resolution of BISTRO data towards the filament network.

Recent analyses of ${\it Herschel}$ observations suggest that in nearby disc galaxies the dust mass opacity at $500 \, {\rm \mu m}$, $\kappa_{500}$, decreases with increasing gas surface density, $\Sigma_{\rm ISM}$ (Clark et al. 2019). This apparent anti-correlation between $\kappa_{500}$ and $\Sigma_{\rm ISM}$ is opposite to the behaviour expected from theoretical dust evolution models; in such models, dust in denser, cooler regions (i.e. regions of increased $\Sigma_{\rm ISM}$) tends to grow and therefore to have increased $\kappa_{500}$. We show, using a toy model, that the presence of a range of dust temperatures along the line of sight can lead to spuriously low estimated values of $\kappa_{500}$. If in regions of higher $\Sigma_{\rm ISM}$ the range of dust temperatures extends to lower values (as seems likely), the magnitude of this effect may be sufficient to explain the apparent anti-correlation between $\kappa_{500}$ and $\Sigma_{\rm ISM}$. Therefore there may not be any need for spatial variation in the intrinsic dust properties that run counter to theoretical expectations.

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}$.

The recent data collected by {\it Herschel} have confirmed that interstellar structures with filamentary shape are ubiquitously present in the Milky Way. Filaments are thought to be formed by several physical mechanisms acting from the large Galactic scales down to the sub-pc fractions of molecular clouds, and they might represent a possible link between star formation and the large-scale structure of the Galaxy. In order to study this potential link, a statistically significant sample of filaments spread throughout the Galaxy is required. In this work we present the first catalogue of $32,059$ candidate filaments automatically identified in the Hi-GAL survey of the entire Galactic Plane. For these objects we determined morphological (length, $l^{a}$, and geometrical shape) and physical (average column density, $N_{\rm H_{2}}$, and average temperature, $T$) properties. We identified filaments with a wide range of properties: 2$'$\,$\leq l^{a}\leq$\, 100$'$, $10^{20} \leq N_{\rm H_{2}} \leq 10^{23}$\,cm$^{-2}$ and $10 \leq T\leq$ 35\,K. We discuss their association with the Hi-GAL compact sources, finding that the most tenuous (and stable) structures do not host any major condensation and we also assign a distance to $\sim 18,400$ filaments for which we determine mass, physical size, stability conditions and Galactic distribution. When compared to the spiral arms structure, we find no significant difference between the physical properties of on-arm and inter-arm filaments. We compared our sample with previous studies, finding that our Hi-GAL filament catalogue represents a significant extension in terms of Galactic coverage and sensitivity. This catalogue represents an unique and important tool for future studies devoted to understanding the filament life-cycle.

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 present the POL-2 850 $\mu$m linear polarization map of the Barnard 1 clump in the Perseus molecular cloud complex from the B-fields In STar-forming Region Observations (BISTRO) survey at the James Clerk Maxwell Telescope. We find a trend of decreasing polarization fraction as a function of total intensity, which we link to depolarization effects towards higher density regions of the cloud. We then use the polarization data at 850 $\mu$m to infer the plane-of-sky orientation of the large-scale magnetic field in Barnard 1. This magnetic field runs North-South across most of the cloud, with the exception of B1-c where it turns more East-West. From the dispersion of polarization angles, we calculate a turbulence correlation length of $5.0 \pm 2.5$ arcsec ($1500$ au), and a turbulent-to-total magnetic energy ratio of $0.5 \pm 0.3$ inside the cloud. We combine this turbulent-to-total magnetic energy ratio with observations of NH$_3$ molecular lines from the Green Bank Ammonia Survey (GAS) to estimate the strength of the plane-of-sky component of the magnetic field through the Davis-Chandrasekhar-Fermi method. With a plane-of-sky amplitude of $120 \pm 60$ $\mu$G and a criticality criterion $\lambda_c = 3.0 \pm 1.5$, we find that Barnard 1 is a supercritical molecular cloud with a magnetic field nearly dominated by its turbulent component.

We present the 850 $\mu$m polarization observations toward the IC5146 filamentary cloud taken using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) and its associated polarimeter (POL-2), mounted on the James Clerk Maxwell Telescope (JCMT), as part of the B-fields In STar forming Regions Observations (BISTRO). This work is aimed at revealing the magnetic field morphology within a core-scale ($\lesssim 1.0$ pc) hub-filament structure (HFS) located at the end of a parsec-scale filament. To investigate whether or not the observed polarization traces the magnetic field in the HFS, we analyze the dependence between the observed polarization fraction and total intensity using a Bayesian approach with the polarization fraction described by the Rice likelihood function, which can correctly describe the probability density function (PDF) of the observed polarization fraction for low signal-to-noise ratio (SNR) data. We find a power-law dependence between the polarization fraction and total intensity with an index of 0.56 in $A_V\sim$ 20--300 mag regions, suggesting that the dust grains in these dense regions can still be aligned with magnetic fields in the IC5146 regions. Our polarization maps reveal a curved magnetic field, possibly dragged by the contraction along the parsec-scale filament. We further obtain a magnetic field strength of 0.5$\pm$0.2 mG toward the central hub using the Davis-Chandrasekhar-Fermi method, corresponding to a mass-to-flux criticality of $\sim$ $1.3\pm0.4$ and an Alfv\'{e}nic Mach number of $<$0.6. These results suggest that gravity and magnetic field is currently of comparable importance in the HFS, and turbulence is less important.

We present the first release of the data and compact-source catalogue for the JCMT Large Program SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE). SCOPE consists of 850-um continuum observations of 1235 Planck Galactic Cold Clumps (PGCCs) made with the Submillimetre Common-User Bolometer Array 2 on the James Clerk Maxwell Telescope. These data are at an angular resolution of 14.4 arcsec, significantly improving upon the 353-GHz resolution of Planck at 5 arcmin, and allowing for a catalogue of 3528 compact sources in 558 PGCCs. We find that the detected PGCCs have significant sub-structure, with 61 per cent of detected PGCCs having 3 or more compact sources, with filamentary structure also prevalent within the sample. A detection rate of 45 per cent is found across the survey, which is 95 per cent complete to Planck column densities of $N_{H_{2}}$ $>$ 5 $\times$ 10$^{21}$ cm$^{-2}$. By positionally associating the SCOPE compact sources with YSOs, the star formation efficiency, as measured by the ratio of luminosity to mass, in nearby clouds is found to be similar to that in the more distant Galactic Plane, with the column density distributions also indistinguishable from each other.

We report 850~$\mu$m dust polarization observations of a low-mass ($\sim$12 $M_{\odot}$) starless core in the $\rho$ Ophiuchus cloud, Ophiuchus C, made with the POL-2 instrument on the James Clerk Maxwell Telescope (JCMT) as part of the JCMT B-fields In STar-forming Region Observations (BISTRO) survey. We detect an ordered magnetic field projected on the plane of sky in the starless core. The magnetic field across the $\sim$0.1~pc core shows a predominant northeast-southwest orientation centering between $\sim$40$^\circ$ to $\sim$100$^\circ$, indicating that the field in the core is well aligned with the magnetic field in lower-density regions of the cloud probed by near-infrared observations and also the cloud-scale magnetic field traced by Planck observations. The polarization percentage ($P$) decreases with an increasing total intensity ($I$) with a power-law index of $-$1.03 $\pm$ 0.05. We estimate the plane-of-sky field strength ($B_{\mathrm{pos}}$) using modified Davis-Chandrasekhar-Fermi (DCF) methods based on structure function (SF), auto-correlation (ACF), and unsharp masking (UM) analyses. We find that the estimates from the SF, ACF, and UM methods yield strengths of 103 $\pm$ 46 $\mu$G, 136 $\pm$ 69 $\mu$G, and 213 $\pm$ 115 $\mu$G, respectively. Our calculations suggest that the Ophiuchus C core is near magnetically critical or slightly magnetically supercritical (i.e. unstable to collapse). The total magnetic energy calculated from the SF method is comparable to the turbulent energy in Ophiuchus C, while the ACF method and the UM method only set upper limits for the total magnetic energy because of large uncertainties.

Young low-mass protostars undergo short phases of high accretion and outburst activity leading to lumpy outflows. Recent observations have shown that the position-velocity and mass-velocity diagrams of such outflows exhibit individual bullet-like features; some of these bullets subscribe to a `Hubble Law' velocity relation, and others are manifest as `Hubble wedges'. In order to explore the origin of these features, we have developed a new episodic outflow model for the SPH code {\sc gandalf}, which mimics the accretion and ejection behaviour of FU Ori type stars. We apply this model to simulations of star formation, invoking two types of initial conditions: spherically symmetric cores in solid-body rotation with $\rho\propto r^{-2}$, and spherically symmetric turbulent cores with density proportional to the density of a Bonnor-Ebert sphere. For a wide range of model parameters, we find that episodic outflows lead to self-regulation of the ejected mass and momentum, and we achieve acceptable results, even with relatively low resolution. Using this model, we find that recently ejected outflow bullets produce a `Hubble wedge' in the position-velocity relation. However, once such a bullet hits the leading shock front, it decelerates and aligns with older bullets to form a `Hubble-law'. Bullets can be identified as bumps in the mass-velocity relation, which can be fit with a power-law, $dM/d\upsilon_{_{\rm RAD}}\propto\upsilon_{_{\rm RAD}}^{-1.5}$.

We present resolution-enhanced images of warm dust at multiple temperatures and opacity index values in the star-forming bubble/HII region, RCW 120. The image set, representing a 4D hypercube of differential column density, was obtained using our Bayesian procedure, ppmap. The cool peripheral material ($\sim16$-22 K) exhibits ragged clumpy structure as noted previously by others. However, at higher temperatures ($\stackrel{>}{_\sim}26$ K) the geometry changes dramatically, showing a bubble boundary which is accurately circular in projection, except for the previously-reported opening in the north. Comparison with Spitzer 8 $\mu$m data suggests that the $\sim26$-30 K dust seen by Herschel resides in the photodissociation region (PDR) surrounding the HII region. Its projected radial profile is consistent with that of a spherical shell, thus arguing against previous suggestions of cylindrical or planar geometry. The inferred geometry is, in fact, consistent with previous interpretations of the HII region as a classical Str\"omgren sphere, except for the fact that the ionising star (CD -38.11636; O8V) is displaced by more than half a radius from its geometric centre. None of the previously published models has satisfactorily accounted for that displacement. It could, however, be explained by proper motion of the O star at $\sim2$-4 km s$^{-1}$ since its formation, possibly due to a cloud-cloud collision. We suggest that the current spherical bubble constitutes the fossilised remnant of the initial expansion of the HII region following the formation of the star, which now continues to flee its formation site.

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.

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

We rehearse the physical and theoretical considerations that define the nature of Brown Dwarfs, in particular the maximum mass for a Brown Dwarf (set by the Hydrogen-Burning Limit) and the minimum mass for a Brown Dwarf (set by the Opacity Limit). We then review the range of mechanisms that have been invoked to explain the formation of Brown Dwarfs and their statistical properties. These include turbulent fragmentation, fragmentation of filaments and discs, dynamical ejection of stellar embryos, and photoerosion. The primary contenders would seem to be turbulent fragmentation and disc fragmentation, and the observations needed to evaluate their relative importance may soon be available.

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.

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 the results of dust emission polarization measurements of Ophiuchus-B (Oph-B) carried out using the Submillimetre Common-User Bolometer Array 2 (SCUBA-2) camera with its associated polarimeter (POL-2) on the James Clerk Maxwell Telescope (JCMT) in Hawaii. This work is part of the B-fields In Star-forming Region Observations (BISTRO) survey initiated to understand the role of magnetic fields in star formation for nearby star-forming molecular clouds. We present a first look at the geometry and strength of magnetic fields in Oph-B. The field geometry is traced over $\sim$0.2 pc, with clear detection of both of the sub-clumps of Oph-B. The field pattern appears significantly disordered in sub-clump Oph-B1. The field geometry in Oph-B2 is more ordered, with a tendency to be along the major axis of the clump, parallel to the filamentary structure within which it lies. The degree of polarization decreases systematically towards the dense core material in the two sub-clumps. The field lines in the lower density material along the periphery are smoothly joined to the large scale magnetic fields probed by NIR polarization observations. We estimated a magnetic field strength of 630$\pm$410 $\mu$G in the Oph-B2 sub-clump using a Davis-Chandeasekhar-Fermi analysis. With this magnetic field strength, we find a mass-to-flux ratio $\lambda$= 1.6$\pm$1.1, which suggests that the Oph-B2 clump is slightly magnetically supercritical.

We present 850 $\mu$m imaging polarimetry data of the $\rho$ Oph-A core taken with the Submillimeter Common-User Bolometer Array-2 (SCUBA-2) and its polarimeter (POL-2), as part of our ongoing survey project, BISTRO (B-fields In STar forming RegiOns). The polarization vectors are used to identify the orientation of the magnetic field projected on the plane of the sky at a resolution of 0.01 pc. We identify 10 subregions with distinct polarization fractions and angles in the 0.2 pc $\rho$ Oph A core; some of them can be part of a coherent magnetic field structure in the $\rho$ Oph region. The results are consistent with previous observations of the brightest regions of $\rho$ Oph-A, where the degrees of polarization are at a level of a few percents, but our data reveal for the first time the magnetic field structures in the fainter regions surrounding the core where the degree of polarization is much higher ($> 5 \%$). A comparison with previous near-infrared polarimetric data shows that there are several magnetic field components which are consistent at near-infrared and submillimeter wavelengths. Using the Davis-Chandrasekhar-Fermi method, we also derive magnetic field strengths in several sub-core regions, which range from approximately 0.2 to 5 mG. We also find a correlation between the magnetic field orientations projected on the sky with the core centroid velocity components.

Understanding the processes that determine the stellar Initial Mass Function (IMF) is a critical unsolved problem, with profound implications for many areas of astrophysics. In molecular clouds, stars are formed in cores, gas condensations which are sufficiently dense that gravitational collapse converts a large fraction of their mass into a star or small clutch of stars. In nearby star-formation regions, the core mass function (CMF) is strikingly similar to the IMF, suggesting that the shape of the IMF may simply be inherited from the CMF. Here we present 1.3 mm observations, obtained with ALMA, the world's largest interferometer, of the active star-formation region W43-MM1, which may be more representative of the Galactic-disk regions where most stars form. The unprecedented resolution of these observations reveals, for the first time, a statistically robust CMF at high masses, with a slope that is markedly shallower than the IMF. This seriously challenges our understanding of the origin of the IMF.

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.

We describe an OctTree algorithm for the MPI-parallel, adaptive mesh-refinement code {\sc FLASH}, which can be used to calculate the gas self-gravity, and also the angle-averaged local optical depth, for treating ambient diffuse radiation. The algorithm communicates to the different processors only those parts of the tree that are needed to perform the tree walk locally. The advantage of this approach is a relatively low memory requirement, important in particular for the optical depth calculation, which needs to process information from many different directions. This feature also enables a general tree-based radiation transport algorithm that will be described in a subsequent paper, and delivers excellent scaling up to at least 1500 cores. Boundary conditions for gravity can be either isolated or periodic, and they can be specified in each direction independently, using a newly developed generalisation of the Ewald method. The gravity calculation can be accelerated with the {\em adaptive block update} technique by partially re-using the solution from the previous time-step. Comparison with the {\sc Flash} internal multi-grid gravity solver shows that tree based methods provide a competitive alternative, particularly for problems with isolated or mixed boundary conditions. We evaluate several multipole acceptance criteria (MACs) and identify a relatively simple APE MAC which provides high accuracy at low computational cost. The optical depth estimates are found to agree very well with those of the {\sc RADMC-3D} radiation transport code, with the tree solver being much faster. Our algorithm is available in the standard release of the {\sc FLASH} code in version 4.0 and later.

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.

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.

We present a comprehensive statistical analysis of star-forming objects located in the vicinities of 1 360 bubble structures throughout the Galactic Plane and their local environments. The compilation of ~70 000 star-forming sources, found in the proximity of the ionized (Hii) regions and detected in both Hi-GAL and GLIMPSE surveys, provided a broad overview of the different evolutionary stages of star-formation in bubbles, from prestellar objects to more evolved young stellar objects (YSOs). Surface density maps of star-forming objects clearly reveal an evolutionary trend where more evolved star-forming objects are found spatially located near the center, while younger star-forming objects are found at the edge of the bubbles. We derived dynamic ages for a subsample of 182 Hii regions for which kinematic distances and radio continuum flux measurements were available. We detect ~80% more star-forming sources per unit area in the direction of bubbles than in the surrounding fields. We estimate ~10% clump formation efficiency (CFE) of Hi-GAL clumps in bubbles, twice the CFE in fields not affected by feedback. We find higher CFE of protostellar clumps in younger bubbles, whose density of the bubble shells is higher. We argue that the formation rate from prestellar to protostellar phase is probably higher during the early stages of the bubble expansion. Evaluation of the fragmentation time inside the shell of bubbles advocates the preexistence of clumps in the medium before the bubble, as supported by numerical simulations. Approximately 23% of the Hi-GAL clumps are found located in the direction of a bubble, with 15% for prestellar clumps and 41% for protostellar clumps. We argue that the high fraction of protostellar clumps may be due to the acceleration of the star-formation process cause by the feedback of the (Hii) bubbles.

We present the first results from the B-fields In STar-forming Region Observations (BISTRO) survey, using the Sub-millimetre Common-User Bolometer Array 2 (SCUBA-2) camera, with its associated polarimeter (POL-2), on the James Clerk Maxwell Telescope (JCMT) in Hawaii. We discuss the survey's aims and objectives. We describe the rationale behind the survey, and the questions which the survey will aim to answer. The most important of these is the role of magnetic fields in the star formation process on the scale of individual filaments and cores in dense regions. We describe the data acquisition and reduction processes for POL-2, demonstrating both repeatability and consistency with previous data. We present a first-look analysis of the first results from the BISTRO survey in the OMC 1 region. We see that the magnetic field lies approximately perpendicular to the famous 'integral filament' in the densest regions of that filament. Furthermore, we see an 'hour-glass' magnetic field morphology extending beyond the densest region of the integral filament into the less-dense surrounding material, and discuss possible causes for this. We also discuss the more complex morphology seen along the Orion Bar region. We examine the morphology of the field along the lower-density north-eastern filament. We find consistency with previous theoretical models that predict magnetic fields lying parallel to low-density, non-self-gravitating filaments, and perpendicular to higher-density, self-gravitating filaments.

We present the results of smoothed particle hydrodynamic simulations investigating the evolution and fragmentation of filaments that are accreting from a turbulent medium. We show that the presence of turbulence, and the resulting inhomogeneities in the accretion flow, play a significant role in the fragmentation process. Filaments which experience a weakly turbulent accretion flow fragment in a two-tier hierarchical fashion, similar to the fragmentation pattern seen in the Orion Integral Shaped Filament. Increasing the energy in the turbulent velocity field results in more sub-structure within the filaments, and one sees a shift from gravity-dominated fragmentation to turbulence-dominated fragmentation. The sub-structure formed in the filaments is elongated and roughly parallel to the longitudinal axis of the filament, similar to the fibres seen in observations of Taurus, and suggests that the fray and fragment scenario is a possible mechanism for the production of fibres. We show that the formation of these fibre-like structures is linked to the vorticity of the velocity field inside the filament and the filament's accretion from an inhomogeneous medium. Moreover, we find that accretion is able to drive and sustain roughly sonic levels of turbulence inside the filaments, but is not able to prevent radial collapse once the filaments become supercritical. However, the supercritical filaments which contain fibre-like structures do not collapse radially, suggesting that fibrous filaments may not necessarily become radially unstable once they reach the critical line-density.

We investigate, using 3D hydrodynamic simulations, the fragmentation of pressure-confined, vertically stratified, self-gravitating gaseous layers. The confining pressure is either thermal pressure acting on both surfaces, or thermal pressure acting on one surface and ram-pressure on the other. In the linear regime of fragmentation, the dispersion relation we obtain agrees well with that derived by Elmegreen & Elmegreen (1978), and consequently deviates from the dispersion relations based on the thin shell approximation (Vishniac 1983) or pressure assisted gravitational instability (W\"unsch et al. 2010). In the non-linear regime, the relative importance of the confining pressure to the self-gravity is a crucial parameter controlling the qualitative course of fragmentation. When confinement of the layer is dominated by external pressure, self- gravitating condensations are delivered by a two-stage process: first the layer fragments into gravitationally bound but stable clumps, and then these clumps coalesce until they assemble enough mass to collapse. In contrast, when external pressure makes a small contribution to confinement of the layer, the layer fragments monolithically into gravitationally unstable clumps and there is no coalescence. This dichotomy persists whether the external pressure is thermal or ram. We apply these results to fragments forming in a shell swept up by an expanding H II region, and find that, unless the swept up gas is quite hot or the surrounding medium has low density, the fragments have low-mass ( ~< 3 M_Sun ), and therefore they are unlikely to spawn stars that are sufficiently massive to promote sequential self-propagating star formation.

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.

In turbulent fragmentation, star formation occurs in condensations created by converging flows. The condensations must be sufficiently massive, dense and cool to be gravitationally unstable, so that they start to contract; {\it and} they must then radiate away thermal energy fast enough for self-gravity to remain dominant, so that they continue to contract. For the metallicities and temperatures in local star forming clouds, this second requirement is only met robustly when the gas couples thermally to the dust, because this delivers the capacity to radiate across the full bandwidth of the continuum, rather than just in a few discrete spectral lines. This translates into a threshold for vigorous star formation, which can be written as a minimum ram-pressure Pcrit ~ 4 10^-11 dyne. Pcrit is independent of temperature, and corresponds to flows with molecular hydrogen number-density nH2 and velocity v satisfying nH2 v^2 > 800 cm^-3 (km/s)^2. This in turn corresponds to a minimum molecular hydrogen column-density for vigorous star formation, NH2crit ~ 4 10^21 cm^-2 (SIGMAcrit ~ 100 MSun pc^-2), and a minimum visual extinction AVcrit ~ 9. The characteristic diameter and line-density for a star-forming filament when this threshold is just exceeded -- a sweet spot for local star formation regions -- are 2Rfil ~ 0.1 pc and mufil ~ 13 MSun pc^-2. The characteristic diameter and mass for a prestellar core condensing out of such a filament are 2Rcore ~ 0.1 pc, and Mcore ~ MSun. We also show that fragmentation of a shock-compressed layer is likely to commence while the convergent flows creating the layer are still ongoing, and we stress that, under this circumstance, the phenomenology and characteristic scales for fragmentation of the layer are fundamentally different from those derived traditionally for pre-existing layers.

Image cubes of differential column density as a function of dust temperature are constructed for Galactic Centre molecular cloud G0.253+0.016 ("The Brick") using the recently described PPMAP procedure. The input data consist of continuum images from the Herschel Space Telescope in the wavelength range 70-500 $\mu$m, supplemented by previously published interferometric data at 1.3 mm wavelength. While the bulk of the dust in the molecular cloud is consistent with being heated externally by the local interstellar radiation field, our image cube shows the presence, near one edge of the cloud, of a filamentary structure whose temperature profile suggests internal heating. The structure appears as a cool ($\sim 14$ K) tadpole-like feature, $\sim 6$ pc in length, in which is embedded a thin spine of much hotter ($\sim$ 40-50 K) material. We interpret these findings in terms of a cool filament whose hot central region is undergoing gravitational collapse and fragmentation to form a line of protostars. If confirmed, this would represent the first evidence of widespread star formation having started within this cloud.

We use smoothed particle hydrodynamic simulations to investigate the growth of perturbations in infinitely long, initially sub-critical but accreting filaments. The growth of these perturbations leads to filament fragmentation and the formation of cores. Most previous work on this subject has been confined to the growth and fragmentation of equilibrium filaments and has found that there exists a preferential fragmentation length scale which is roughly 4 times the filament's diameter. Our results show a more complicated dispersion relation with a series of peaks linking perturbation wavelength and growth rate. These are due to gravo-acoustic oscillations along the longitudinal axis during the sub-critical phase of growth. The positions of the peaks in growth rate have a strong dependence on both the mass accretion rate onto the filament and the temperature of the gas. When seeded with a multi-wavelength density power spectrum there exists a clear preferred core separation equal to the largest peak in the dispersion relation. Our results allow one to estimate a minimum age for a filament which is breaking up into regularly spaced fragments, as well as a maximum accretion rate. We apply the model to observations of filaments in Taurus by Tafalla & Hacar (2015) and find accretion rates consistent with those estimated by Palmeirim et al. (2013).

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

Massive clumps within the swept-up shells of bubbles, like that surrounding the galactic HII region RCW 120, have been interpreted in terms of the Collect and Collapse (C&C) mechanism for triggered star formation. The cold, dusty clumps surrounding RCW 120 are arranged in an almost spherical shell and harbour many young stellar objects. By performing high-resolution, three-dimensional SPH simulations of HII regions expanding into fractal molecular clouds, we investigate whether the formation of massive clumps in dense, swept-up shells necessarily requires the C&C mechanism. In a second step, we use RADMC-3D to compute the synthetic dust continuum emission from our simulations, in order to compare them with observations of RCW 120 made with APEX-LABOCA at 870 micron. We show that a distribution of clumps similar to the one seen in RCW 120 can readily be explained by a non-uniform initial molecular cloud structure. Hence, a shell-like configuration of massive clumps does not imply that the C&C mechanism is at work. Rather, we find a hybrid form of triggering, which combines elements of C&C and Radiatively Driven Implosion (RDI). In addition, we investigate the reliability of deriving clump masses from their 870 micron emission. We find that for clumps with more than 100 M_sun the observational estimates are accurate to within a factor of two and that, even at these long wavelengths, it is important to account for the radiative heating from triggered, embedded protostars.

We use Smoothed Particle Hydrodynamic simulations of cold, uniform density, self-gravitating filaments, to investigate their longitudinal collapse timescales; these timescales are important because they determine the time available for a filament to fragment into cores. A filament is initially characterised by its line-mass, $\mu$, its radius, $R$ (or equivalently its density $\rho\!=\!\mu/\pi R^2$), and its aspect ratio, $A\;\,(\equiv Z/R$, where $Z$ is its half-length). The gas is only allowed to contract longitudinally, i.e. parallel to the symmetry axis of the filament (the $z$-axis). Pon et al. (2012) have considered the global dynamics of such filaments analytically. They conclude that short filaments ($A\! < \!5$) collapse along the $z$-axis more-or-less homologously, on a time-scale $t_{_{\rm HOM}} \sim 0.44\,A\,(G\rho)^{-1/2}$; in contrast, longer filaments ($A\! > \!5$) undergo end-dominated collapse, i.e. two dense clumps form at the ends of the filament and converge on the centre sweeping up mass as they go, on a time-scale $t_{_{\rm END}} \sim 0.98\,A^{1/2}\,(G\rho)^{-1/2}$. Our simulations do not corroborate these predictions. First, for all $A\! > \!2$, the collapse time satisfies a single equation \[t_{_{\rm COL}}\;\sim\;(0.49+0.26A)(G\rho)^{-1/2}\,,\] which for large $A$ is much longer than the Pon et al. prediction. Second, for all $A\! > \!2$, the collapse is end-dominated. Third, before being swept up, the gas immediately ahead of an end-clump is actually accelerated outwards by the gravitational attraction of the approaching clump, resulting in a significant ram pressure. For high aspect ratio filaments the end-clumps approach an asymptotic inward speed, due to the fact that they are doing work both accelerating and compressing the gas they sweep up. Pon et al. appear to have neglected the outward acceleration and its consequences.

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 review current theories for the origin of the Stellar Initial Mass Function (IMF) with particular focus on the extent to which the IMF can be considered universal across various environments. To place the issue in an observational context, we summarize the techniques used to determine the IMF for different stellar populations, the uncertainties affecting the results, and the evidence for systematic departures from universality under extreme circumstances. We next consider theories for the formation of prestellar cores by turbulent fragmentation and the possible impact of various thermal, hydrodynamic and magneto-hydrodynamic instabilities. We address the conversion of prestellar cores into stars and evaluate the roles played by different processes: competitive accretion, dynamical fragmentation, ejection and starvation, filament fragmentation and filamentary accretion flows, disk formation and fragmentation, critical scales imposed by thermodynamics, and magnetic braking. We present explanations for the characteristic shapes of the Present-Day Prestellar Core Mass Function and the IMF and consider what significance can be attached to their apparent similarity. Substantial computational advances have occurred in recent years, and we review the numerical simulations that have been performed to predict the IMF directly and discuss the influence of dynamics, time-dependent phenomena, and initial conditions.

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.

The time evolution of the probability density function (PDF) of the mass density is formulated and solved for systems in free-fall using a simple appoximate function for the collapse of a sphere. We demonstrate that a pressure-free collapse results in a power-law tail on the high-density side of the PDF. The slope quickly asymptotes to the functional form $\mathrm{P}_v(\rho)\propto\rho^{-1.54}$ for the (volume-weighted) PDF and $\mathrm{P}_m(\rho)\propto\rho^{-0.54}$ for the corresponding mass-weighted distribution. From the simple approximation of the PDF we derive analytic descriptions for mass accretion, finding that dynamically quiet systems with narrow density PDFs lead to retarded star formation and low star formation rates. Conversely, strong turbulent motions that broaden the PDF accelerate the collapse causing a bursting mode of star formation. Finally, we compare our theoretical work with observations. The measured star formation rates are consistent with our model during the early phases of the collapse. Comparison of observed column density PDFs with those derived from our model suggests that observed star-forming cores are roughly in free-fall.

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.

Infrared shells and bubbles are ubiquitous in the Galaxy and can generally be associated with HII regions formed around young, massive stars. In this paper, we use high-resolution 3D SPH simulations to explore the effect of a single O7 star emitting photons at 10^49 1/s and located at the centre of a molecular cloud with mass 10^4 M_sun and radius 6.4 pc; the internal structure of the cloud is characterised by its fractal dimension, D (with 2.0 <= D <= 2.8), and the variance of its (log-normal) density distribution, sigma_0^2 (with 0.36 <= sigma_0^2 <= 1.42). Our study focuses on the morphology of the swept-up cold gas and the distribution and statistics of the resulting star formation. If the fractal dimension is low, the border of the HII region is dominated by extended shell-like structures, and these break up into a small number of massive high-density clumps which then spawn star clusters; star formation occurs relatively quickly, and delivers somewhat higher stellar masses. Conversely, if the fractal dimension is high, the border of the HII region is dominated by a large number of pillars and cometary globules, which contain compact dense clumps and tend to spawn single stars or individual multiple systems; star formation occurs later, the stellar masses are somewhat lower, and the stars are more widely distributed.

Observations indicate that the central portions of the Present-Day Prestellar Core Mass Function (CMF) and the Stellar Initial Mass Function (IMF) both have approximately log-normal shapes, but that the CMF is displaced to higher mass than the IMF by a factor F = 4+/-1. This has lead to suggestions that the shape of the IMF is directly inherited from the shape of the CMF - and therefore, by implication, that there is a self-similar mapping from the CMF onto the IMF. If we assume a self-similar mapping, it follows (i) that F = N0/eta, where eta is the mean fraction of a core's mass that ends up in stars, and N0 is the mean number of stars spawned by a single core; and (ii) that the stars spawned by a single core must have an approximately log-normal distribution of relative masses, with universal standard deviation sigma0. Observations can be expected to deliver ever more accurate estimates of F, but this still leaves a degeneracy between eta and N0; and sigma0 is also unconstrained by observation. Here we show that these parameters can be estimated by invoking binary statistics. Specifically, if (a) each core spawns one long-lived binary system, and (b) the probability that a star of mass M is part of this long-lived binary is proportional to M^alpha, current observations of the binary frequency as a function of primary mass and the distribution of mass ratios strongly favour eta = 1.0 +/- 0.3, N0 = 4.3 +/- 0.4, sigma0 = 0.3 +/- 0.03 and alpha = 0.9 +/- 0.6. The mapping from CMF to IMF is not necessarily self-similar - there are many possible motivations for a non self-similar mapping - but if it is not, then the shape of the IMF cannot be inherited from the CMF. Given the limited observational constraints currently available and the ability of a self-similar mapping to satisfy them, the possibility that the shape of the IMF is inherited from the CMF cannot be ruled out at this juncture.

We use three-dimensional smoothed particle hydrodynamics simulations to study the non-linear thin shell instability (NTSI) in supersonic colliding flows. We show that for flows with monochromatic perturbations and for flows with white-noise perturbations, growth speeds approximate quite well to the analytic predictions of Vishniac (1994). For flows with subsonic turbulence, growth speeds match Vishniac's predictions only at short wavelengths where the turbulence is weaker. We find that supersonic turbulence, of a lower Mach number than the colliding flows, completely suppresses the NTSI. Our results provide a diagnostic for identifying the presence of the NTSI in colliding flows with turbulence.

Numerical simulations of star formation frequently rely on the implementation of sink particles, (a) to avoid expending computational resource on the detailed internal physics of individual collapsing protostars, (b) to derive mass functions, binary statistics and clustering kinematics (and hence to make comparisons with observation), and (c) to model radiative and mechanical feedback; sink particles are also used in other contexts, for example to represent accreting black holes in galactic nuclei. We present a new algorithm for creating and evolving sink particles in SPH simulations, which appears to represent a significant improvement over existing algorithms {\refrpt -- particularly in situations where sinks are introduced after the gas has become optically thick to its own cooling radiation and started to heat up by adiabatic compression}. (i) It avoids spurious creation of sinks. (ii) It regulates the accretion of matter onto a sink so as to mitigate non-physical perturbations in the vicinity of the sink. (iii) Sinks accrete matter, but the associated angular momentum is transferred back to the surrounding medium. With the new algorithm -- and modulo the need to invoke sufficient resolution to capture the physics preceding sink formation -- the properties of sinks formed in simulations are essentially independent of the user-defined parameters of sink creation, or the number of SPH particles used.

Protostars grow in mass by accreting material through their discs, and this accretion is initially their main source of luminosity. The resulting radiative feedback heats the environments of young protostars, and may thereby suppress further fragmentation and star formation. There is growing evidence that the accretion of material onto protostars is episodic rather than continuous; most of it happens in short bursts that last up to a few hundred years, whereas the intervals between these outbursts of accretion could be thousands of years. We have developed a model to include the effects of episodic accretion in simulations of star formation. Episodic accretion results in episodic radiative feedback, which heats and temporarily stabilises the disc, suppressing the growth of gravitational instabilities. However, once an outburst has been terminated, the luminosity of the protostar is low, and the disc cools rapidly. Provided that there is enough time between successive outbursts, the disc may become gravitationally unstable and fragment. The model suggests that episodic accretion may allow disc fragmentation if (i) the time between successive outbursts is longer than the dynamical timescale for the growth of gravitational instabilities (a few kyr), and (ii) the quiescent accretion rate onto the protostar is sufficiently low (at most a few times 1e-7 Msun/yr). We also find that after a few protostars form in the disc, their own episodic accretion events shorten the intervals between successive outbursts, and sup- press further fragmentation, thus limiting the number of objects forming in the disc. We conclude that episodic accretion moderates the effect of radiative feedback from young protostars on their environments, and, under certain conditions, allows the formation of low-mass stars, brown dwarfs, and planetary-mass objects by fragmentation of protostellar discs.

The role of feedback from massive stars is believed to be a key element in the evolution of molecular clouds. We use high-resolution 3D SPH simulations to explore the dynamical effects of a single O7 star located at the centre of a molecular cloud with mass 10^4M_sun and radius 6.4pc. The initial internal structure of the cloud is characterised by its fractal dimension, D=2.0 - 2.8, and its log-normal density PDF. (i) As regards star formation, in the short term ionising feedback is positive, in the sense that star formation occurs much more quickly in gas that is compressed by the high pressure of the ionised gas. However, in the long term ionising feedback is negative, in the sense that most of the cloud is dispersed with an outflow rate of up to ~0.01M_sun/yr, on a timescale comparable with the sound-crossing time for the ionised gas (~1-2Myr), and triggered star formation is therefore limited to a few percent of the cloud's mass. (ii) As regards the morphology of the ionisation fronts (IFs) bounding the HII region and the systematics of outflowing gas, we distinguish two regimes. For low D<=2.2, the initial cloud is dominated by large-scale structures, so the neutral gas tends to be swept up into a few extended coherent shells, and the ionised gas blows out through a few large holes between these shells; we term these HII regions "shell-dominated". Conversely, for high D>=2.6, the initial cloud is dominated by small-scale structures, and these are quickly overrun by the advancing IF, thereby producing neutral pillars whilst the ionised gas blows out through a large number of small holes between the pillars; we term these HII regions "pillar-dominated". (iii) As regards the injection of bulk kinetic energy, by ~1Myr, the expansion of the HII region has delivered a rms velocity of ~6km/s; this represents less than 0.1% of the total energy radiated by the O7 star.

Several mechanisms have been proposed for the formation of brown dwarfs, but there is as yet no consensus as to which -- if any -- are operative in nature. Any theory of brown dwarf formation must explain the observed statistics of brown dwarfs. These statistics are limited by selection effects, but they are becoming increasingly discriminating. In particular, it appears (a) that brown dwarfs that are secondaries to Sun-like stars tend to be on wide orbits, $a\ga 100\,{\rm AU}$ (the Brown Dwarf Desert), and (b) that these brown dwarfs have a significantly higher chance of being in a close ($a\la 10\,{\rm AU}$) binary system with another brown dwarf than do brown dwarfs in the field. This then raises the issue of whether these brown dwarfs have formed {\it in situ}, i.e. by fragmentation of a circumstellar disc; or have formed elsewhere and subsequently been captured. We present numerical simulations of the purely gravitational interaction between a close brown-dwarf binary and a Sun-like star. These simulations demonstrate that such interactions have a negligible chance ($<0.001$) of leading to the close brown-dwarf binary being captured by the Sun-like star. Making the interactions dissipative by invoking the hydrodynamic effects of attendant discs might alter this conclusion. However, in order to explain the above statistics, this dissipation would have to favour the capture of brown-dwarf binaries over single brown-dwarfs, and we present arguments why this is unlikely. The simplest inference is that most brown-dwarf binaries -- and therefore possibly also most single brown dwarfs -- form by fragmentation of circumstellar discs around Sun-like protostars, with some of them subsequently being ejected into the field.

In a previous paper we identified cores within infrared dark clouds (IRDCs). We regarded those without embedded sources as the least evolved, and labelled them starless. Here we identify the most isolated starless cores and model them using a three-dimensional, multi-wavelength, Monte Carlo, radiative transfer code. We derive the cores' physical parameters and discuss the relation between the mass, temperature, density, size and the surrounding interstellar radiation field (ISRF) for the cores. The masses of the cores were found not to correlate with their radial size or central density. The temperature at the surface of a core was seen to depend almost entirely on the level of the ISRF surrounding the core. No correlation was found between the temperature at the centre of a core and its local ISRF. This was seen to depend, instead, on the density and mass of the core.

We have used data taken as part of the Herschel infrared Galactic Plane survey (Hi-GAL) to study 3171 infrared-dark cloud (IRDC) candidates that were identified in the mid-infrared (8 {\mu}m) by Spitzer (we refer to these as 'Spitzer-dark' regions). They all lie in the range l=300 - 330 \circ and |b| 6 1 \circ. Of these, only 1205 were seen in emission in the far-infrared (250-500 {\mu}m) by Herschel (we call these 'Herschel-bright' clouds). It is predicted that a dense cloud will not only be seen in absorption in the mid-infrared, but will also be seen in emission in the far-infrared at the longest Herschel wavebands (250-500 {\mu}m). If a region is dark at all wavelengths throughout the mid-infrared and far-infrared, then it is most likely to be simply a region of lower background infrared emission (a 'hole in the sky'). Hence, it appears that previous surveys, based on Spitzer and other mid-infrared data alone, may have over-estimated the total IRDC population by a factor of 2. This has implications for estimates of the star formation rate in IRDCs in the Galaxy.We studied the 1205 Herschel-bright IRDCs at 250 {\mu}m, and found that 972 of them had at least one clearly defined 250-{\mu}m peak, indicating that they contained one or more dense cores. Of these, 653 (67 per cent) contained an 8-{\mu}m point source somewhere within the cloud, 149 (15 per cent) contained a 24-{\mu}m point source but no 8-{\mu}m source, and 170 (18 per cent) contained no 24-{\mu}m or 8-{\mu}m point sources. We use these statistics to make inferences about the lifetimes of the various evolutionary stages of IRDCs.

Recent discoveries of strongly misaligned transiting exoplanets pose a challenge to the established planet formation theory which assumes planetary systems to form and evolve in isolation. However, the fact that the majority of stars actually do form in star clusters raises the question how isolated forming planetary systems really are. Besides radiative and tidal forces the presence of dense gas aggregates in star-forming regions are potential sources for perturbations to protoplanetary discs or systems. Here we show that subsequent capture of gas from large extended accretion envelopes onto a passing star with a typical circumstellar disc can tilt the disc plane to retrograde orientation, naturally explaining the formation of strongly inclined planetary systems. Furthermore, the inner disc regions may become denser, and thus more prone to speedy coagulation and planet formation. Pre-existing planetary systems are compressed by gas inflows leading to a natural occurrence of close-in misaligned hot Jupiters and short-period eccentric planets. The likelihood of such events mainly depends on the gas content of the cluster and is thus expected to be highest in the youngest star clusters.

It is speculated that the accretion of material onto young protostars is episodic. We present a computational method to include the effects of episodic accretion in radiation hydrodynamic simulations of star formation. We find that during accretion events protostars are "switched on", heating and stabilising the discs around them. However, these events typically last only a few hundred years, whereas the intervals in between them may last for a few thousand years. During these intervals the protostars are effectively "switched off", allowing gravitational instabilities to develop in their discs and induce fragmentation. Thus, episodic accretion promotes disc frag- mentation, enabling the formation of low-mass stars, brown dwarfs and planetary-mass objects. The frequency and the duration of episodic accretion events may be responsible for the low-mass end of the IMF, i.e. for more than 60% of all stars.

The collapse of weakly turbulent prestellar cores is a critical stage in the process of star formation. Being highly non-linear and stochastic, the outcome of collapse can only be explored theoretically by performing large ensembles of numerical simulations. Standard practice is to quantify the initial turbulent velocity field in a core in terms of the amount of turbulent energy (or some equivalent) and the exponent in the power spectrum (n \equiv -d log Pk /d log k). In this paper, we present a numerical study of the influence of the details of the turbulent velocity field on the collapse of an isolated, weakly turbulent, low-mass prestellar core. We show that, as long as n > 3 (as is usually assumed), a more critical parameter than n is the maximum wavelength in the turbulent velocity field, {\lambda}_MAX. This is because {\lambda}_MAX carries most of the turbulent energy, and thereby influences both the amount and the spatial coherence of the angular momentum in the core. We show that the formation of dense filaments during collapse depends critically on {\lambda}_MAX, and we explain this finding using a force balance analysis. We also show that the core only has a high probability of fragmenting if {\lambda}_MAX > 0.5 R_CORE (where R_CORE is the core radius); that the dominant mode of fragmentation involves the formation and break-up of filaments; and that, although small protostellar disks (with radius R_DISK <= 20 AU) form routinely, more extended disks are rare. In turbulent, low-mass cores of the type we simulate here, the formation of large, fragmenting protostellar disks is suppressed by early fragmentation in the filaments.

We propose an evolutionary path for prestellar cores on the radius-mass diagram, which is analogous to stellar evolutionary paths on the Hertzsprung-Russell Diagram. Using James Clerk Maxwell Telescope (JCMT) observations of L1688 in the Ophiuchus star-forming complex, we analyse the HCO+ (J=4\rightarrow3) spectral line profiles of prestellar cores. We find that of the 58 cores observed, 14 show signs of infall in the form of a blue-asymmetric double-peaked line profile. These 14 cores all lie beyond the Jeans mass line for the region on a radius-mass plot. Furthermore another 10 cores showing tentative signs of infall, in their spectral line profile shapes, appear on or just over the Jeans mass line. We therefore propose the manner in which a prestellar core evolves across this diagram. We hypothesise that a core is formed in the low-mass, low-radius region of the plot. It then accretes quasistatically, increasing in both mass and radius. When it crosses the limit of gravitational instability it begins to collapse, decreasing in radius, towards the region of the diagram where protostellar cores are seen.

We present simulations of initially stable isothermal clouds exposed to ionising radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN and an HEALPix-based photo-ionisation algorithm to simulate the propagation of the ionising radiation and the resulting dynamical evolution of the cloud. We find that the incident ionising flux, $\Phi_{_{\rm LyC}}$, is the critical parameter determining the cloud evolution. At moderate fluxes, a large fraction of the cloud mass is converted into stars. As the flux is increased, the fraction of the cloud mass that is converted into stars and the mean masses of the individual stars both decrease. Very high fluxes simply disperse the cloud. Newly-formed stars tend to be concentrated along the central axis of the cloud (i.e. the axis pointing in the direction of the incident flux). For given cloud parameters, the time, $t_{_\star}$, at which star formation starts is proportional to $\Phi_{_{\rm LyC}}^{-1/3}$. The pattern of star formation found in the simulations is similar to that observed in bright-rimmed clouds.

We study the influence of gas metallicity, turbulence, and non-equilibrium chemistry on the evolution of the two-phase interstellar medium (warm and cold atomic phases), and thereby constrain the initial conditions for star formation prevailing in turbulent gas. We perform high-resolution simulations in three dimensions, including a realistic non-equilibrium treatment of the ionization state of the gas, and examine both driven and decaying turbulence. This allows us to explore variations in the metallicity Z. In this paper, we study solar metallicity, Z=Z_sun, and low metallicity, Z=1d-3 Z_sun, gas. For driven, large-scale turbulence, we find that the influence of the metallicity on the amount of mass in the cold gas component is small. However, in decaying turbulent conditions this picture is much changed. While cold regions survive in the case of solar metallicity, they are quickly heated and dispersed in low-metallicity gas. This result suggests that star formation can be suppressed in environments of low metallicity, unless a strong turbulent driver is acting on time scales shorter than a few turbulent crossing times. Inter alia this finding could explain the overall inefficient star formation as well as the burst-like mode of star formation found in metal-poor, gas-rich systems like dwarf galaxies.

A star acquires much of its mass by accreting material from a disc. Accretion is probably not continuous but episodic. We have developed a method to include the effects of episodic accretion in simulations of star formation. Episodic accretion results in bursts of radiative feedback, during which a protostar is very luminous, and its surrounding disc is heated and stabilised. These bursts typically last only a few hundred years. In contrast, the lulls between bursts may last a few thousand years; during these lulls the luminosity of the protostar is very low, and its disc cools and fragments. Thus, episodic accretion enables the formation of low-mass stars, brown dwarfs and planetary-mass objects by disc fragmentation. If episodic accretion is a common phenomenon among young protostars, then the frequency and duration of accretion bursts may be critical in determining the low-mass end of the stellar initial mass function.

We present SEREN, a new hybrid Smoothed Particle Hydrodynamics and N-body code designed to simulate astrophysical processes such as star and planet formation. It is written in Fortran 95/2003 and has been parallelised using OpenMP. SEREN is designed in a flexible, modular style, thereby allowing a large number of options to be selected or disabled easily and without compromising performance. SEREN uses the conservative `grad-h' formulation of SPH, but can easily be configured to use traditional SPH or Godunov SPH. Thermal physics is treated either with a barotropic equation of state, or by solving the energy equation and modelling the transport of cooling radiation. A Barnes-Hut tree is used to obtain neighbour lists and compute gravitational accelerations efficiently, and an hierarchical time-stepping scheme is used to reduce the number of computations per timestep. Dense gravitationally bound objects are replaced by sink particles, to allow the simulation to be evolved longer, and to facilitate the identification of protostars and the compilation of stellar and binary properties. At the termination of a hydrodynamical simulation, SEREN has the option of switching to a pure N-body simulation, using a 4th-order Hermite integrator, and following the ballistic evolution of the sink particles (e.g. to determine the final binary statistics once a star cluster has relaxed). We describe in detail all the algorithms implemented in SEREN and we present the results of a suite of tests designed to demonstrate the fidelity of SEREN and its performance and scalability. Further information and additional tests of SEREN can be found at the web-page http://www.astro.group.shef.ac.uk/seren.

The densest infrared dark clouds (IRDCs) may represent the earliest observable stage of high-mass star formation. These clouds are very cold, hence they emit mainly at far-infrared and sub-mm wavelengths. For the first time, Herschel has provided multi-wavelength, spatially resolved observations of cores within IRDCs, which, when combined with radiative transfer modelling, can constrain their properties, such as mass, density profile and dust temperature. We use a 3D, multi-wavelength Monte Carlo radiative transfer code to model in detail the emission from six cores in three typical IRDCs seen in the Hi-GAL survey (G030.50+00.95, G031.03+00.26 and G031.03+00.76), and thereby to determine the properties of these cores and compare them with their low-mass equivalents. We found masses ranging from 90 to 290 solar masses with temperatures from 8 to 11K at the centre of each core and 18 to 28K at the surface. The maximum luminosity of an embedded star within each core was calculated, and we rule out the possibility of significant high mass star formation having yet occurred in three of our cores.

A large fraction of brown dwarfs and low-mass H-burning stars may form by gravitational fragmentation of protostellar discs. We explore the conditions for disc fragmentation and we find that they are satisfied when a disc is large enough (>100 AU) so that its outer regions can cool efficiently, and it has enough mass to be gravitationally unstable, at such radii. We perform radiative hydrodynamic simulations and show that even a disc with mass 0.25 Msun and size 100 AU fragments. The disc mass, radius, and the ratio of disc-to-star mass (Mdisc/Mstar~0.36) are smaller than in previous studies. We find that fragmenting discs decrease in mass and size within a few 10^4 yr of their formation, since a fraction of their mass, especially outside 100 AU is consumed by the new stars and brown dwarfs that form. Fragmenting discs end up with masses ~0.001-0.1 Msun, and sizes ~20-100 AU. On the other hand, discs that are marginally stable live much longer. We produce simulated images of fragmenting discs and find that observing discs that are undergoing fragmentation is possible using current (e.g. IRAM-PdBI) and future (e.g. ALMA) interferometers, but highly improbable due to the short duration of this process. Comparison with observations shows that many observed discs may be remnants of discs that have fragmented at an earlier stage. However, there are only a few candidates that are possibly massive and large enough to currently be gravitationally unstable. The rarity of massive (>0.2 Msun), extended (>100 AU) discs indicates either that such discs are highly transient (i.e. form, increase in mass becoming gravitationally unstable due to infall of material from the surrounding envelope, and quickly fragment), or that their formation is suppressed (e.g. by magnetic fields). We conclude that current observations of early-stage discs cannot exclude the mechanism of disc fragmentation.

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.

We use SPH simulations to investigate the gravitational fragmentation of expanding shells through the linear and non--linear regimes. The results are analysed using spherical harmonic decomposition to capture the initiation of structure during the linear regime; the potential-based method of Smith et al. (2009) to follow the development of clumps in the mildly non-linear regime; and sink particles to capture the properties of the final bound objects during the highly non-linear regime. In the early, mildly non--linear phase of fragmentation, we find that the clump mass function still agrees quite well with the mass function predicted by the analytic model. However, the sink mass function is quite different, in the sense of being skewed towards high-mass objects. This is because, once the growth of a condensation becomes non-linear, it tends to be growing non-competitively from its own essentially separate reservoir; we call this Oligarchic Accretion.

A large fraction of brown dwarfs and low-mass stars may form by gravitational fragmentation of relatively massive (a few 0.1 Msun), extended (a few hundred AU) discs around Sun-like stars. We present an ensemble of radiative hydrodynamic simulations that examine the conditions for disc fragmentation. We demonstrate that this model can explain the low-mass IMF, the brown dwarf desert, and the binary properties of low-mass stars and brown dwarfs. Observing discs that are undergoing fragmentation is possible but very improbable, as the process of disc fragmentation is short lived (discs fragment within a few thousand years).

We present simulations of initially stable isothermal clouds exposed to ionizing radiation from a discrete external source, and identify the conditions that lead to radiatively driven implosion and star formation. We use the Smoothed Particle Hydrodynamics code SEREN (Hubber et al. 2010) and the HEALPix-based photoionization algorithm described in Bisbas et al. (2009). We find that the incident ionizing flux is the critical parameter determining the evolution: high fluxes simply disperse the cloud, whereas low fluxes trigger star formation. We find a clear connection between the intensity of the incident flux and the parameters of star formation.

Most stars are born in clusters and the resulting gravitational interactions between cluster members may significantly affect the evolution of circumstellar discs and therefore the formation of planets and brown dwarfs. Recent findings suggest that tidal perturbations of typical circumstellar discs due to close encounters may inhibit rather than trigger disc fragmentation and so would seem to rule out planet formation by external tidal stimuli. However, the disc models in these calculations were restricted to disc radii of 40 AU and disc masses below 0.1 M_sun. Here we show that even modest encounters can trigger fragmentation around 100 AU in the sorts of massive (~0.5 M_sun), extended (>=100 AU) discs that are observed around young stars. Tidal perturbation alone can do this, no disc-disc collision is required. We also show that very-low-mass binary systems can form through the interaction of objects in the disc. In our computations, otherwise non-fragmenting massive discs, once perturbed, fragment into several objects between about 0.01 and 0.1 M_sun, i.e., over the whole brown dwarf mass range. Typically these orbit on highly eccentric orbits or are even ejected. While probably not suitable for the formation of Jupiter- or Neptune-type planets, our scenario provides a possible formation mechanism for brown dwarfs and very massive planets which, interestingly, leads to a mass distribution consistent with the canonical substellar IMF. As a minor outcome, a possible explanation for the origin of misaligned extrasolar planetary systems is discussed.

We demonstrate the use of the 3D Monte Carlo radiative transfer code PHAETHON to model infrared-dark clouds (IRDCs) that are externally illuminated by the interstellar radiation field (ISRF). These clouds are believed to be the earliest observed phase of high-mass star formation, and may be the high-mass equivalent of lower-mass prestellar cores. We model three different cases as examples of the use of the code, in which we vary the mass, density, radius, morphology and internal velocity field of the IRDC. We show the predicted output of the models at different wavelengths chosen to match the observing wavebands of Herschel and Spitzer. For the wavebands of the long- wavelength SPIRE photometer on Herschel, we also pass the model output through the SPIRE simulator to generate output images that are as close as possible to the ones that would be seen using SPIRE. We then analyse the images as if they were real observations, and compare the results of this analysis with the results of the radiative transfer models. We find that detailed radiative transfer modelling is necessary to accurately determine the physical parameters of IRDCs (e.g. dust temperature, density profile). This method is applied to study G29.55+00.18, an IRDC observed by the Herschel Infrared Galactic Plane survey (Hi-GAL), and in the future it will be used to model a larger sample of IRDCs from the same survey.

We study analytically the development of gravitational instability in an expanding shell having finite thickness. We consider three models for the radial density profile of the shell: (i) an analytic uniform-density model, (ii) a semi-analytic model obtained by numerical solution of the hydrostatic equilibrium equation, and (iii) a 3D hydrodynamic simulation. We show that all three profiles are in close agreement, and this allows us to use the first model to describe fragments in the radial direction of the shell. We then use non-linear equations describing the time-evolution of a uniform oblate spheroid to derive the growth rates of shell fragments having different sizes. This yields a dispersion relation which depends on the shell thickness, and hence on the pressure confining the shell. We compare this dispersion relation with the dispersion relation obtained using the standard thin-shell analysis, and show that, if the confining pressure is low, only large fragments are unstable. On the other hand, if the confining pressure is high, fragments smaller than predicted by the thin-shell analysis become unstable. Finally, we compare the new dispersion relation with the results of 3D hydrodynamic simulations, and show that the two are in good agreement.

The formation process of binary stars and multiple systems is poorly understood. Here, we seek to determine the typical outcome of protostellar collapse and to constrain models of binary formation by core fragmentation during collapse, using high-resolution millimeter continuum imaging of very young (Class 0) protostars observed at the beginning of the main accretion phase. We carried out a pilot high-resolution study of 5 Class 0 objects, using the most extended (A) configuration of the IRAM PdBI at 1.3 mm, which allow us to probe the multiplicity of Class 0 protostars down to separations a ~50 AU and circumstellar mass ratios q ~0.07. We show that our PdBI observations revealed only wide (>1500 AU) protobinary systems and/or outflow-generated features. When combined with previous millimeter interferometric observations of Class 0 protostars, our pilot PdBI study tentatively suggests that the binary fraction in the ~ 75-1000 AU range increases from the Class 0 to the Class I stage. It also seems to argue against purely hydrodynamic models of binary star formation. We briefly discuss possible alternative scenarios to reconcile the low multiplicity rate of Class 0 protostars on small scales with the higher binary fraction observed at later evolutionary stages.

A large fraction of the observed brown dwarfs may form by gravitational fragmentation of unstable discs. This model reproduces the brown dwarf desert, and provides an explanation the existence of planetary-mass objects and for the binary properties of low-mass objects. We have performed an ensemble of radiative hydrodynamic simulations and determined the statistical properties of the low-mass objects produced by gravitational fragmentation of discs. We suggest that there is a population of brown dwarfs loosely bound on wide orbits (100-5000 AU) around Sun-like stars that surveys of brown dwarf companions should target. Our simulations also indicate that planetary-mass companions to Sun-like stars are unlikely to form by disc fragmentation.

Thermodynamics play an important role in determining the way a protostellar disc fragments to form planets, brown dwarfs and low-mass stars. We explore the effect that different treatments of radiative transfer have in simulations of fragmenting discs. Three prescriptions for the radiative transfer are used, (i) the diffusion approximation of Stamatellos et al., (ii) the barotropic equation of state (EOS) of Goodwin et al., and (iii) the barotropic EOS of Bate et al. The barotropic approximations capture the general evolution of the density and temperature at the centre of each proto-fragment but (i) they do not make any adjustments the particular circumstances of a proto-fragment forming in the disc, and (ii) they do not take into account thermal inertia effects that are important for fast-forming proto-fragments in the outer disc region. As a result, the number of fragments formed in the disc and their properties are different, when a barotropic EOS is used. This is important not only for disc studies but also for simulations of collapsing turbulent clouds, as in many cases in such simulations stars form with discs that subsequently fragment. We also examine the difference in the way proto-fragments condense out in the disc at different distances from the central star using the diffusion approximation and following the collapse of each proto-fragment until the formation of the second core (~10^{-3} g/cm3). We find that proto-fragments forming closer to the central star tend to form earlier and evolve faster from the first to the second core than proto-fragments forming in the outer disc region. The former have a large pool of material in the inner disc region that they can accrete from and grow in mass. The latter accrete more slowly and they are hotter because they generally form in a quick abrupt event.

We investigate the collapse and fragmentation of low-mass, trans-sonically turbulent prestellar cores, using SPH simulations. The initial conditions are slightly supercritical Bonnor-Ebert spheres, all with the same density profile, the same mass (M_O=6.1 Msun) and the same radius (R_O=17,000 AU), but having different initial turbulent velocity fields. Four hundred turbulent velocity fields have been generated, all scaled so that the mean Mach number is M=1. Then a subset of these, having a range of net angular momenta, j, has been evolved. The evolution of these turbulent cores is not strongly correlated with j. Instead it is moderated by the formation of filamentary structures due to converging turbulent flows. A high fraction (~ 82%) of the protostars forming from turbulent cores are attended by protostellar accretion discs, but only a very small fraction (~16%) of these discs is sufficiently cool and extended to develop non-linear gravitational instabilities and fragment.

Abridged: We use three-dimensional SPH simulations to investigate the collapse of low-mass prestellar cores and the formation and early evolution of protostellar discs. The initial conditions are slightly supercritical Bonnor-Ebert spheres in rigid rotation. The core mass and initial radius are held fixed at M_O=6.1 M_sun and R_O=17,000 AU, and the only parameter that we vary is the initial angular speed \Omega_O. Protostellar discs forming from cores with \Omega_O<1.35 10d-13 1/s have radii between 100 and 300 AU and are quite centrally concentrated; due to heating by gas infall onto the disc and accretion onto the central object, they are also quite warm, T>100 K, and therefore stable against gravitational fragmentation. In contrast, more rapidly rotating cores form discs which are less concentrated and cooler, and have radii between 400 and 1000 AU; as a consequence they are prone to gravitational fragmentation and the formation of multiple systems. We derive a criterion that predicts whether a rigidly rotating core having given M_O, R_O and \Omega_O will produce a protostellar disc which fragments whilst material is still infalling from the core envelope. We then apply this criterion to core samples for which M_O, R_O and \Omega_O have been estimated observationally. We conclude that the observed cores are stable against fragmentation at this stage, due to their low angular speeds and the heat delivered at the accretion shock where the infalling material hits the disc.

We investigate the gravitational fragmentation of expanding shells in the context of the linear thin--shell analysis. We make use of two very different numerical schemes; the FLASH Adaptive Mesh Refinement code and a version of the Benz Smoothed Particle Hydrodynamics code. We find that the agreement between the two codes is excellent. We use our numerical results to test the thin--shell approximation and we find that the external pressure applied to the shell has a strong effect on the fragmentation process. In cases where shells are not pressure--confined, the shells thicken as they expand and hydrodynamic flows perpendicular to the plane of the shell suppress fragmentation at short wavelengths. If the shells are pressure--confined internally and externally, so that their thickness remains approximately constant during their expansion, the agreement with the analytical solution is better.

We describe a new algorithm for including the dynamical effects of ionizing radiation in SPH simulations, and we present several examples of how the algorithm can be applied to problems in star formation. We use the HEALPix software to tessellate the sky and to solve the equation of ionization equilibrium along a ray towards each of the resulting tesserae. We exploit the hierarchical nature of HEALPix to make the algorithm adaptive, so that fine angular resolution is invoked only where it is needed, and the computational cost is kept low. We present simulations of (i) the spherically symmetric expansion of an HII region inside a uniform-density, non--self-gravitating cloud; (ii) the spherically symmetric expansion of an HII region inside a uniform-density, self-gravitating cloud; (iii) the expansion of an off-centre HII region inside a uniform-density, non--self-gravitating cloud, resulting in rocket acceleration and dispersal of the cloud; and (iv) radiatively driven compression and ablation of a core overrun by an HII region. The new algorithm provides the means to explore and evaluate the role of ionizing radiation in regulating the efficiency and statistics of star formation.

We perform SPH simulations of the collapse and fragmentation of low-mass cores having different initial levels of turbulence (alpha_turb=0.05,0.10,0.25). We use a new treatment of the energy equation which captures the transport of cooling radiation against opacity due to both dust and gas (including the effects of dust sublimation, molecules, and H^- ions). We also perform comparison simulations using a standard barotropic equation of state. We find that -- when compared with the barotropic equation of state -- our more realistic treatment of the energy equation results in more protostellar objects being formed, and a higher proportion of brown dwarfs; the multiplicity frequency is essentially unchanged, but the multiple systems tend to have shorter periods (by a factor ~3), higher eccentricities, and higher mass ratios. The reason for this is that small fragments are able to cool more effectively with the new treatment, as compared with the barotropic equation of state. We find that the process of fragmentation is often bimodal. The first protostar to form is usually, at the end, the most massive, i.e. the primary. However, frequently a disc-like structure subsequently forms round this primary, and then, once it has accumulated sufficient mass, quickly fragments to produce several secondaries. We believe that this delayed fragmentation of a disc-like structure is likely to be an important source of very low-mass hydrogen-burning stars and brown dwarfs.

A new means of incorporating radiative transfer into smoothed particle hydrodynamics (SPH) is introduced, which builds on the success of two previous methods - the polytropic cooling approximation as devised by Stamatellos et al (2007), and flux limited diffusion (e.g. Mayer et al 2007). This hybrid method preserves the strengths of its individual components, while removing the need for atmosphere matching or other boundary conditions to marry optically thick and optically thin regions. The code uses a non-trivial equation of state to calculate temperatures and opacities of SPH particles, which captures the effects of molecular hydrogen dissociation, atomic hydrogen ionisation, He0 and He+ ionisation, ice evaporation, dust sublimation, molecular absorption, bound-free and free-free transitions and electron scattering. The method is tested in several scenarios, including: (1) the evolution of a 0.07 solar mass protoplanetary disc surrounding a 0.5 solar mass star; (2) the collapse of a 1 solar mass protostellar cloud, and (3) the thermal relaxation of temperature fluctuations in a static homogeneous sphere.

We suggest that a high proportion of brown dwarfs are formed by gravitational fragmentation of massive extended discs around Sun-like stars. Such discs should arise frequently, but should be observed infrequently, precisely because they fragment rapidly. By performing an ensemble of radiation-hydrodynamic simulations, we show that such discs fragment within a few thousand years, and produce mainlybrown dwarf (BDs) stars, but also planetary mass (PM) stars and very low-mass hydrogen-burning (HB) stars. Most of the the PM stars and BDs are ejected by mutual interactions. We analyse the statistical properties of these stars, and compare them with observations. After a few hundred thousand years the Sun-like primary is typically left with a close low-mass HB companion, and two much wider companions: a low-mass HB star and a BD star, or a BD-BD binary. There is a BD desert extending out to at least ~100 AU; this is because BDs tend to be formed further out than low-mass HB stars, and then they tend to be scattered even further out, or even into the field. BDs form with discs of a few Mj and radii of a few tens of AU, and they are more likely to retain these discs if they remain bound to the primary star. Binaries form by pairing of the newly-formed stars in the disc, giving a low-mass binary fraction of ~0.16. These binaries include close and wide BD/BD binaries and BD/PM binaries. BDs that remain as companions to Sun-like stars are more likely to be in BD/BD binaries than are BDs ejected into the field. Disc fragmentation is a robust mechanism; even if only a small fraction of Sun-like stars host the required massive extended discs,this mechanism can produce all the PM stars observed, most of the BD stars, and a significant proportion of the very low-mass HB stars.

We suggest that a high proportion of brown dwarfs are formed by gravitational fragmentation of massive, extended discs around Sun-like stars. We argue that such discs should arise frequently, but should be observed infrequently, precisely because they fragment rapidly. By performing an ensemble of radiation-hydrodynamic simulations, we show that such discs typically fragment within a few thousand years to produce mainly brown dwarfs (including planetary-mass brown dwarfs) and low-mass hydrogen-burning stars. Subsequently most of the brown dwarfs are ejected by mutual interactions. We analyse the properties of these objects that form by disc fragmentation, and compare them with observations.

We present a new method of incorporating radiative transfer into Smoothed Particle Hydrodynamics (SPH). There have been many recent attempts at radiative transfer in SPH (Stamatellos et al 2005, 2005, Mayer et al 2007, Whitehouse and Bate 2006), however these are becoming increasingly complex, with some methods requiring the photosphere to be mapped (which is often of non-trivial geometric shape), and extra conditions to be applied there (matching atmospheres as in Cai et al (2008), or specifying cooling at the photosphere as in Mayer et al (2007)). The method of identifying the photosphere is usually a significant addition to the total simulation runtime, and often requires extra free parameters, the changing of which will affect the final results. Our method is not affected by such concerns, as the photosphere is constructed implicitly by the algorithm without the need for extra free parameters. The algorithm used is a synergy of two current formalisms for radiative effects: a) the polytropic cooling formalism proposed by Stamatellos et al (2007), and b) flux-limited diffusion, used by many authors to simulate radiation transport in the optically thick regime (e.g. Mayer et al 2007). We present several tests of this method: (1) The evolution of a 0.07 solar mass protoplanetary disc around a 0.5 solarmass star (Pickett et al 2003, Mejia et al 2005, Boley et al 2006, Cai et al 2008); (2) The collapse of a non-rotating 1 solar mass molecular cloud (Masunaga & Inutsuka 2000, Stamatellos et al 2007); (3) The thermal relaxation of temperature fluctuations in an static homogeneous sphere (Masunaga et al 1998, Spiegel 1957, Stamatellos et al 2007)

Context: It is generally presumed that the outflows from a YSO are directed close to its rotation axis (i.e. along its angular momentum vector and orthogonal to any attendant accretion disc). Many YSOs are formed from dense prestellar cores embedded in filaments, and therefore the relative orientations of outflows and filaments may place useful constraints on the dynamics of core formation. Aims: We explore this possibility, from the viewpoint of what it may tell us about the angular momentum delivered to a core forming in a filament. We stress that we are not here addressing the issue of the relationship of filaments and outflows to the prevailing magnetic field direction, although this is evidently also an interesting issue. Results: Using the limited data available, and neglecting any possible selection effects or correlations between nearby outflows, we infer that the observed outflows have a tendency to be orthogonal to the filaments that contain their driving sources. Specfically, in the cases where the directions of the filaments and outflows are well defined, we infer statistically that 72% of outflows are within $45^{\rm o}$ of being orthogonal to the filament, and only 28% are within $45^{\rm o}$ of being parallel to the filament. Conclusions: This suggests that the prestellar cores which spawned the YSOs driving the observed outflows had angular momenta which were approximately orthogonal to the filaments out of which the cores formed. We briefly discuss the implications of this for two proposed core formation mechanisms.

We present results of hydrodynamic simulations of star formation triggered by cloud-cloud collisions. During the early stages of star formation, low-mass objects form by gravitational instabilities in protostellar discs. A number of these low-mass objects are in the sub-stellar mass range, including a few objects of planetary mass. The disc instabilities that lead to the formation of low-mass objects in our simulations are the product of disc-disc interactions and/or interactions between the discs and their surrounding gas.

Gravitational fragmentation has been proposed as a mechanism for the formation of giant planets in close orbits around solar-type stars. However, it is debatable whether this mechanism can function in the inner regions (R<40 AU) of real discs. We use a newly developed method for treating the energy equation and the equation of state, which accounts for radiative transfer effects in SPH simulations of circumstellar discs. The different chemical and internal states of hydrogen and the properties of dust at different densities and temperatures (ice coated dust grains at low temperatures, ice melting, dust sublimation) are all taken into account by the new method.We present radiative hydrodynamic simulations of the inner regions of massive circumstellar discs and examine two cases: (i) a disc irradiated by a cool background radiation field (T_bgr=10K)and (ii) a disc heated by radiation from its central star (T_bgr~1/R). In neither case does the disc fragment: in the former because it cannot cool fast enough and in the latter because it is not gravitationally unstable. Our results (a) corroborate previous numerical results using different treatments for the hydrodynamics and the radiative transfer, and (b) confirm our own earlier analytic predictions. We conclude that disc fragmentation is unlikely to be able to produce giant planets around solar-type stars at radii <40 AU.

Stars form from dense molecular cores, and the mass function of these cores (the CMF) is often found to be similar to the form of the stellar initial mass function (IMF). This suggests that the form of the IMF is the result of the form of the CMF. However, most stars are thought to form in binary and multiple systems, therefore the relationship between the IMF and the CMF cannot be trivial. We test two star formation scenarios - one in which all stars form as binary or triple systems, and one in which low-mass stars form in a predominantly single mode. We show that from a log-normal CMF, similar to those observed, and expected on theoretical grounds, the model in which all stars form as multiples gives a better fit to the IMF.

We suggest that low-mass hydrogen-burning stars like the Sun should sometimes form with massive extended discs; and we show, by means of radiation hydrodynamic simulations, that the outer parts of such discs (R>100 AU) are likely to fragment on a dynamical timescale (10^3 to $10^4 yr), forming low-mass companions: principally brown dwarfs (BDs), but also very low-mass hydrogen-burning stars and planetary-mass objects. A few of the BDs formed in this way remain attached to the primary star, orbiting at large radii. The majority are released into the field, by interactions amongst themselves; in so doing they acquire only a low velocity dispersion (<2 km/s), and therefore they usually retain small discs, capable of registering an infrared excess and sustaining accretion. Some BDs form close BD/BD binaries, and these binaries can survive ejection into the field. This BD formation mechanism appears to avoid some of the problems associated with the `embryo ejection' scenario, and to answer some of the questions not yet answered by the `turbulent fragmentation' scenario.

We introduce and test a new and highly efficient method for treating the thermal and radiative effects influencing the energy equation in SPH simulations of star formation. The method uses the density, temperature and gravitational potential of each particle to estimate a mean optical depth, which then regulates the particle's heating and cooling. The method captures -- at minimal computational cost -- the effects of (i) the rotational and vibrational degrees of freedom of H2, H2 dissociation, H0 ionisation, (ii) opacity changes due to ice mantle melting, sublimation of dust, molecular lines, H-, bound-free and free-free processes and electron scattering; (iv) external irradiation; and (v) thermal inertia. The new algorithm reproduces the results of previous authors and/or known analytic solutions. The computational cost is comparable to a standard SPH simulation with a simple barotropic equation of state. The method is easy to implement, can be applied to both particle- and grid-based codes, and handles optical depths 0<tau<10^{11}.

This paper describes a James Clerk Maxwell Telescope (JCMT) legacy survey that has been awarded roughly 500 hrs of observing time to be carried out from 2007 to 2009. In this survey we will map with SCUBA-2 (Submillimetre Common User Bolometer Array 2) almost all of the well-known low-mass and intermediate-mass star-forming regions within 0.5 kpc that are accessible from the JCMT. Most of these locations are associated with the Gould Belt. From these observations we will produce a flux-limited snapshot of star formation near the Sun, providing a legacy of images, as well as point-source and extended-source catalogues, over almost 700 square degrees of sky. The resulting images will yield the first catalogue of prestellar and protostellar sources selected by submillimetre continuum emission, and should increase the number of known sources by more than an order of magnitude. We will also obtain CO maps with the array receiver HARP (Heterodyne Array Receiver Programme), in three CO isotopologues, of a large typical sample of prestellar and protostellar sources. We will then map the brightest hundred sources with the SCUBA-2 polarimeter (POL-2), producing the first statistically significant set of polarization maps in the submillimetre. The images and source catalogues will be a powerful reference set for astronomers, providing a detailed legacy archive for future telescopes, including ALMA, Herschel and JWST.

We estimate the dust temperatures of the clumps in the rho Oph main cloud taking into account the 3D geometry of the region, and external heating from the interstellar radiation field and from HD147879,a nearby luminous B2V star, which is believed to dominate the radiation field in the region. We find that the regions where prestellar cores are observed (i.e. at optical visual extinctions >7 mag) are colder than ~10-11K. These dust temperatures are smaller than those which previous studies of the same region have assumed. We use the new dust temperatures to estimate the masses of the prestellar cores in the rho Oph main cloud from mm observations, and we find core masses that are larger than previous estimates by a factor of ~2-3. This affects the core mass function (CMF) of the region; we find that the mass at which the core mass spectrum steepens from a slope alpha~1.5 to a slope alpha~2.5 has moved from ~0.5 Msun to ~1Msun. In contrast with the CMF in other star forming regions (e.g. Orion), there is no indication for a turnover down to the completeness limit (~0.2Msun), but the CMF may flatten at around ~0.4Msun. We generalize our results to the prestellar cores in Taurus and in Orion. In Taurus the ambient radiation field heating the prestellar cores is believed to be weaker than than that in rho Oph. Hence, the dust temperatures of the cores in Taurus are expected to be below ~ 10-11 K.In Orion the radiation field is believed to be 10^3 times stronger than the standard interstellar radiation field. Based on this assumption we estimate that the dust temperatures of the prestellar cores in Orion are around ~20-30K.

We investigate, by means of numerical simulations, the phenomenology of star formation triggered by low-velocity collisions between low-mass molecular clumps. The simulations are performed using an SPH code which satisfies the Jeans condition by invoking On-the-Fly Particle Splitting. Clumps are modelled as stable truncated (non-singular) isothermal, i.e. Bonnor-Ebert, spheres. Collisions are characterised by M_0 (clump mass), b (offset parameter, i.e. ratio of impact parameter to clump radius), and M (Mach Number, i.e. ratio of collision velocity to effective post-shock sound speed). The gas subscribes to a barotropic equation of state, which is intended to capture (i) the scaling of pre-collision internal velocity dispersion with clump mass, (ii) post-shock radiative cooling, and (iii) adiabatic heating in optically thick protostellar fragments. The efficiency of star formation is found to vary between 10% and 30% in the different collisions studied and it appears to increase with decreasing M_0, and/or decreasing b, and/or increasing M. For b<0.5 collisions produce shock compressed layers which fragment into filaments. Protostellar objects then condense out of the filaments and accrete from them. The resulting accretion rates are high, 1 to 5 x 10^{-5} M_sun yr^{-1}, for the first 1 to 3 x 10^4 yrs. The densities in the filaments, n >~ 5 x 10^5 cm^{-3}, are sufficient that they could be mapped in NH_3 or CS line radiation, in nearby star formation regions.

Context: The principal mechanism by which brown dwarfs form, and its relation to the formation of higher-mass (i.e. hydrogen-burning) stars, is poorly understood. Aims: We advocate a new model for the formation of brown dwarfs. Methods: In this model, brown dwarfs are initially binary companions, formed by gravitational fragmentation of the outer parts (R > 100 AU) of protostellar discs around low-mass hydrogen-burning stars. Most of these binaries are then gently disrupted by passing stars to create a largely single population of brown dwarfs and low-mass hydrogen-burning stars. Results: This idea is consistent with the excess of binaries found in low-density pre-main sequence populations, like that in Taurus, where they should survive longer than in denser clusters. Conclusions: If brown dwarfs form in this way, as companions to more massive stars, the difficulty of forming very low-mass prestellar cores is avoided. Since the disrupted binaries will tend to be those involving low-mass components and wide orbits, and since disruption will be due to the gentle tides of passing stars (rather than violent N-body interactions in small-N sub-clusters), the liberated brown dwarfs will have velocity dispersions and spatial distributions very similar to higher-mass stars, and they will be able to retain discs, and thereby to sustain accretion and outflows. Thus the problems associated with the ejection and turbulence mechanisms can be avoided. This model implies that most, possibly all, stars and brown dwarfs form in binary or multiple systems.

Context: There is a need to improve the fidelity of SPH simulations of self-gravitating gas dynamics. Aims: We remind users of SPH that, if smoothing lengths are adjusted so as to keep the number of neighbours, ${\cal N}$, in the range ${\cal N}_{_{\rm NEIB}}\pm\Delta{\cal N}_{_{\rm NEIB}}$, the tolerance, $\Delta{\cal N}_{_{\rm NEIB}}$, should be set to zero, as first noted by Nelson & Papaloizou. We point out that this is a very straightforward and computationally inexpensive constraint to implement. Methods: We demonstrate this by simulating acoustic oscillations of a self-gravitating isentropic monatomic gas-sphere (cf. Lucy), using ${\cal N}_{_{\rm TOT}}\sim6,000$ particles and ${\cal N}_{_{\rm NEIB}}=50$. Results: We show that there is a marked reduction in the rates of numerical dissipation and diffusion as $\Delta{\cal N}_{_{\rm NEIB}}$ is reduced from 10 to zero. Moreover this reduction incurs a very small computational overhead. Conclusions: We propose that this should become a standard test for codes used in simulating star formation. It is a highly relevant test, because pressure waves generated by the switch from approximate isothermality to approximate adiabaticity play a critical role in the fragmentation of collapsing prestellar cores. Since many SPH simulations in the literature use ${\cal N}_{_{\rm NEIB}}=50$ and $\Delta{\cal N}_{_{\rm NEIB}}\geq10$, their results must be viewed with caution.

We conduct Smoothed Particle Hydrodynamics simulations of the `collect and collapse' scenario (Elmegreen and Lada, 1977) for star formation triggered by an expanding HII region. We simulate the evolution of a spherical uniform molecular cloud with an ionising source at its centre. The gas in the cloud is self-gravitating, although the cloud is prevented from globally collapsing. We find that the shell driven by the HII region fragments to form numerous self--gravitating objects. We repeat our calculations at four numerical resolutions to ensure that they are converged. We compare our results to the analytical model of Whitworth et al, 1994 and show that our simulations and the predictions of Whitworth et al are in good agreement in the sense that the shell fragments at the time and radius predicted by Whitworth et al to within 20% and 25% respectively. Most of the fragments produced in our two highest resolution calculations are approximately half the mass of those predicted, but this conclusion is robust against both numerical resolution and the presence of random noise (local fluctuations in density of a factor of $\sim2$) in the initial gas distribution. We conclude that such noise has little impact on the fragmentation process.

Four dense cores, L1582A, L1689A, B133 and B68, classified as prestellar in terms of the absence of detectable NIR emission, are observed at radio wavelengths to investigate whether they nurture very young protostars. No definite young protostars were discovered in any of the four cores observed. A few radio sources were discovered close to the observed cores, but these are most likely extragalactic sources or YSOs unrelated to the cores observed. In L1582A we discovered a weak radio source near the centre of the core with radio characteristics and offset from the peak of the submillimeter emission similar to that of the newly discovered protostar in the core L1014, indicating a possible protostellar nature for this source. This needs to be confirmed with near- and/or mid-infrared observations (e.g. with Spitzer). Hence based on the current observations we are unable to confirm unequivocally that L1582A is starless. In L1689A a possible 4.5-sigma radio source was discovered at the centre of the core, but needs to be confirmed with future observations. In B133 a weak radio source, possibly a protostar, was discovered at the edge of the core on a local peak of the core submm emission, but no source was detected at the centre of the core. Thus, B133 is probably starless, but may have a protostar at its edge. In B68 no radio sources were discovered inside or at the edge of the core, and thus B68 is indeed starless. Four more radio sources with spectral indices characteristic of young protostars were discovered outside the cores but within the extended clouds in which these cores reside. Conclusions: We conclude that the number of cores misclassified as prestellar is probably very small and does not significantly alter the estimated lifetime of the prestellar phase.

Our first aim is to calculate the minimum mass for Primary Fragmentation in a variety of potential star-formation scenarios, i.e. (i) hierarchical fragmentation of a 3-D medium; (ii) one-shot, 2-D fragmentation of a shock-compressed layer; (iii) fragmentation of a circumstellar disc. Our second aim is to evaluate the role of H2 dissociation in facilitating Secondary Fragmentation and thereby producing close, low-mass binaries. Results: (i)For contemporary, local star formation, the minimum mass for Primary Fragmentation is in the range 0.001-0.004Msun, irrespective of the scenario considered. (ii)Circumstellar discs are only able to radiate fast enough to undergo Primary Fragmentation in their cool outer parts (R>100AU). Therefore brown dwarfs (BDs) should have difficulty forming by Primary Fragmentation at R<30AU, explaining the Brown Dwarf Desert.Conversely, Primary Fragmentation at R>100AU could be the source of brown dwarfs in wide orbits, and could explain why massive discs with Rd>100AU are rarely seen.(iii)H2 dissociation can lead to collapse and Secondary Fragmentation, thereby converting primary fragments into close, low-mass binaries, with semi-major axes a~5AU(Msystem/0.1Msun), in good agreement with observation; in this case, the minimum mass for Primary Fragmentation becomes a minimum system mass, rather than a minimum stellar mass.(iv)Any primary fragment can undergo Secondary Fragmentation, producing a close low-mass binary, provided only that the fragment is spinning. Secondary Fragmentation is therefore most likely in fragments formed in the outer parts of discs, and this could explain why a BD in a wide orbit about a Sun-like star has a greater likelihood of having a BD companion than a BD in the field -as seems to be observed.

We present observations taken with SCUBA on the JCMT of the Horsehead Nebula in Orion (B33), at wavelengths of 450 and 850 \mum. We see bright emission from that part of the cloud associated with the photon-dominated region (PDR) at the `top' of the horse's head, which we label B33-SMM1. We characterise the physical parameters of the extended dust responsible for this emission, and find that B33-SMM1 contains a more dense core than was previously suspected. We compare the SCUBA data with data from the Infrared Space Observatory (ISO) and find that the emission at 6.75-\mum is offset towards the west, indicating that the mid-infrared emission is tracing the PDR while the submillimetre emission comes from the molecular cloud core behind the PDR. We calculate the virial balance of this core and find that it is not gravitationally bound but is being confined by the external pressure from the HII region IC434, and that it will either be destroyed by the ionising radiation, or else may undergo triggered star formation. Furthermore we find evidence for a lozenge-shaped clump in the `throat' of the horse, which is not seen in emission at shorter wavelengths. We label this source B33-SMM2 and find that it is brighter at submillimetre wavelengths than B33-SMM1. SMM2 is seen in absorption in the 6.75-\mum ISO data, from which we obtain an independent estimate of the column density in excellent agreement with that calculated from the submillimetre emission. We calculate the stability of this core against collapse and find that it is in approximate gravitational virial equilibrium. This is consistent with it being a pre-existing core in B33, possibly pre-stellar in nature, but that it may also eventually undergo collapse under the effects of the HII region.

We investigate the effect of the turbulent power spectrum (P(k) \propto k^{-n}, with n=3, 4 or 5) on the fragmentation of low-mass cores, by means of SPH simulations. We adopt initial density profiles and low levels of turbulence based on observation, and for each n-value we conduct an ensemble of simulations with different initial seeds for the turbulent velocity field, so as to obtain reasonable statistics. We find that when power is concentrated at larger scales (i.e. for larger n), more protostellar objects form and there is a higher proportion of low-mass stars and brown dwarfs. This is in direct contrast with the recent results of Delgado Donate et al., presumably because they adopted much higher levels of turbulence.

We review five mechanisms for forming brown dwarfs: (i) turbulent fragmentation of molecular clouds, producing very low-mass prestellar cores by shock compression; (ii) collapse and fragmentation of more massive prestellar cores; (iii) disc fragmentation; (iv) premature ejection of protostellar embryos from their natal cores; and (v) photo-erosion of pre-existing cores overrun by HII regions. These mechanisms are not mutually exclusive. Their relative importance probably depends on environment, and should be judged by their ability to reproduce the brown-dwarf IMF, the distribution and kinematics of newly formed brown dwarfs, the binary statistics of brown dwarfs, the ability of brown dwarfs to retain discs, and hence their ability to sustain accretion and outflows. This will require more sophisticated numerical modelling than is presently possible, in particular more realistic initial conditions and more realistic treatments of radiation transport, angular momentum transport and magnetic fields. We discuss the minimum mass for brown dwarfs, and how brown dwarfs should be distinguished from planets.

Lada et al. recently presented a detailed near-infrared extinction map of Globule G2 in the Coalsack molecular cloud complex, showing that this starless core has a well-defined central extinction minimum. We propose a model for G2 in which a rapid increase in external pressure is driving an approximately symmetric compression wave into the core. The rapid increase in external pressure could arise because the core has recently been assimilated by the Coalsack cloud complex, or because the Coalsack has recently been created by two large-scale converging flows. The resulting compression wave has not yet converged on the centre of the core, so there is a central rarefaction. The compression wave has increased the density in the swept-up gas by about a factor of ten, and accelerated it inwards to speeds of order $0.4 {\rm km} {\rm s}^{-1}$. It is shown that even small levels of initial turbulence destroy the ring seen in projection almost completely. In the scenario of strong external compression that we are proposing this implies that the initial turbulent energy in this globule is such that $E_{{\rm turb}} / E_{{\rm grav}} \le 2 %$. Protostar formation should occur in about $40,000 {\rm years}$.

We examine the spatial distribution of brown dwarfs produced by the decay of small-N stellar systems as expected from the embryo ejection scenario. We model a cluster of several hundred stars grouped into 'cores' of a few stars/brown dwarfs. These cores decay, preferentially ejecting their lowest-mass members. Brown dwarfs are found to have a wider spatial distribution than stars, however once the effects of limited survey areas and unresolved binaries are taken into account it can be difficult to distinguish between clusters with many or no ejections. A large difference between the distributions probably indicates that ejections have occurred, however similar distributions sometimes arise even with ejections. Thus the spatial distribution of brown dwarfs is not necessarily a good discriminator between ejection and non-ejection scenarios.

We review four mechanisms for forming brown dwarfs: (i) turbulent fragmentation (producing very low-mass prestellar cores); (ii) gravitational instabilities in discs; (iii) dynamical ejection of stellar embryos from their placental cores; and (iv) photo-erosion of pre-existing cores in HII regions. We argue (a) that these are simply the mechanisms of low-mass star formation, and (b) that they are not mutually exclusive. If, as seems possible, all four mechanisms operate in nature, their relative importance may eventually be constrained by their ability to reproduce the binary statistics of brown dwarfs, but this will require fully 3-D radiative magneto-hydrodynamic simulations.

Jeans showed analytically that, in an infinite uniform-density isothermal gas, plane-wave perturbations collapse to dense sheets if their wavelength, $\lambda$, satisfies $\lambda > \lambda_{_{\rm JEANS}} = (\pi a^2 / G \rho_{_0})^{1/2}$ (where $a$ is the isothermal sound speed and $\rho_{_0}$ is the unperturbed density); in contrast, perturbations with smaller $\lambda$ oscillate about the uniform density state. Here we show that Smoothed Particle Hydrodynamics reproduces these results well, even when the diameters of the SPH particles are twice the wavelength of the perturbation. Our simulations are performed in 3-D with initially settled (i.e. non-crystalline) distributions of particles. Therefore there exists the seed noise for artificial fragmentation, but it does not occur. We conclude that, although there may be -- as with any numerical scheme -- `skeletons in the SPH cupboard', a propensity to fragment artificially is evidently not one of them.

Class 0 objects, which are thought to be the youngest protostars, are identified in terms of NIR or radio emission and/or the presence of molecular outflows. We present combined hydrodynamic and radiative transfer simulations of the collapse of a star-forming molecular core, which suggest two criteria for identifying dense cores with deeply embedded very young protostars that may not be observable in the NIR or radio with current telescopes. We find that cores with protostars are relatively warm (T>15K) with their SEDs peaking at wavelengths <170 micron, and they tend to appear circular.

We study the transition from a prestellar core to a Class 0 protostar, using SPH to simulate the dynamical evolution, and a Monte Carlo radiative transfer code to generate the SED and isophotal maps. For a prestellar core illuminated by the standard interstellar radiation field, the luminosity is low and the SED peaks at ~190 micron. Once a protostar has formed, the luminosity rises (due to a growing contribution from accretion onto the protostar) and the peak of the SED shifts to shorter wavelengths (~80-100 micron). However, by the end of the Class 0 phase, the accretion rate is falling, the luminosity has decreased, and the peak of the SED shifts back towards longer wavelengths (90-150 micron). In our simulations, the density of material around the protostar remains sufficiently high well into the Class 0 phase that the protostar only becomes visible in the NIR if it is displaced from the centre dynamically. Raw submm/mm maps of Class 0 protostars tend to be much more centrally condensed than those of prestellar cores. However, when convolved with a typical telescope beam, the difference in central concentration is less marked, although the Class 0 protostars appear more circular. Our results suggest that, if a core is deemed to be prestellar on the basis of having no associated IRAS source, no cm radio emission, and no outflow, but it has a circular appearance and an SED which peaks at wavelengths below ~170 micron, it may well contain a very young Class 0 protostar.

We combine a Monte Carlo radiative transfer code with an SPH code, so that -- assuming thermal equilibrium -- we can calculate dust-temperature fields, spectral energy distributions, and isophotal maps, for the individual time-frames generated by an SPH simulation. On large scales, the radiative transfer cells (RT cells) are borrowed from the tree structure built by the SPH code, and are chosen so that their size -- and hence the resolution of the calculated temperature field -- is comparable with the resolution of the density field. We refer collectively to these cubic RT cells as the 'global grid'. The code is tested and found to treat externally illuminated dust configurations very well. However, when there are embedded discrete sources, i.e. stars, these produce very steep local temperature gradients which can only be modelled properly if -- in the immediate vicinity of, and centred on, each embedded star -- we supplement the global grid with a 'star grid' of closely spaced concentric RT cells.

The Monte Carlo method is a powerful tool for performing radiative equilibrium calculations, even in complex geometries. The main drawback of the standard Monte Carlo radiative equilibrium methods is that they require iteration, which makes them numerically very demanding. Bjorkman & Wood recently proposed a frequency distribution adjustment scheme, which allows radiative equilibrium Monte Carlo calculations to be performed without iteration, by choosing the frequency of each re-emitted photon such that it corrects for the incorrect spectrum of the previously re-emitted photons. Although the method appears to yield correct results, we argue that its theoretical basis is not completely transparent, and that it is not completely clear whether this technique is an exact rigorous method, or whether it is just a good and convenient approximation. We critically study the general problem of how an already sampled distribution can be adjusted to a new distribution by adding data points sampled from an adjustment distribution. We show that this adjustment is not always possible, and that it depends on the shape of the original and desired distributions, as well as on the relative number of data points that can be added. Applying this theorem to radiative equilibrium Monte Carlo calculations, we provide a firm theoretical basis for the frequency distribution adjustment method of Bjorkman & Wood, and we demonstrate that this method provides the correct frequency distribution through the additional requirement of radiative equilibrium. We discuss the advantages and limitations of this approach, and show that it can easily be combined with the presence of additional heating sources and the concept of photon weighting. However, the method may fail if small dust grains are included... (abridged)

We present a simple model of binary star formation based on the assumption that rotating prestellar cores collapse to form rings and these rings then fragment into protostars. We assume that each ring spawns a small number (N <= 6) of protostars, and that the condensation of the protostars is sufficiently rapid that they can subsequently be treated as point masses. We simulate the dynamical evolution of a ring of N stars and record the properties of the single stars, binaries and higher multiples that form as a result of the dissolution of the ring. The masses of the individual stars in a ring are drawn from a log-normal distribution with dispersion \sigma. This part of the model is perfomed for many different realizations of the ring, to obtain good statistics. The second part of the model is to convolve the results of these dimensionless simulations, first with the distribution of core masses, which yields the distributions of multiplicity,mass ratio and eccentricity, as a function of primary mass; and second with the distribution of core angular momenta, which yields the distributions of semi-major axis and period, again as a function of primary mass. Using the observed distribution of core masses, and the observed distribution of core angular momenta, our model is able to reproduce the observed IMF, the observed high multiplicity frequency of pre-Main Sequence stars, the observed distribution of separations, and -- for long-period systems -- the observed distributions of eccentricity and mass-ratio, provided we invoke N = 4 or 5 and $\sigma = 0.6$. We presume that the reduced multiplicity frequency in the field, compared with young clusters, is the result of dynamical interactions between stars formed in different cores but the same cluster, following ring dissolution.

We review and evaluate four mechanisms for forming brown dwarves: (i) dynamical ejection of a stellar embryo from its placental prestellar core; (ii) opacity-limited fragmentation of a shock-compressed layer; (iii) gravitational instabilities in discs, triggered by impulsive interactions with other discs or naked stars; and (iv) photo-erosion of pre-existing cores. All these mechanisms can produce free-floating brown dwarves, but only (ii) and (iii) are likely to produce brown dwarves in multiple systems, and (i) has difficulty delivering brown dwarves with discs.

We present a new analysis of the minimum mass for star formation, based on opacity-limited fragmentation. Our analysis differs from the standard one, which considers hierarchical fragmentation of a 3-D medium, and yields $M_{_{\rm MIN}} \sim 0.007 {\rm to} 0.010 M_\odot$ for Population I star formation. Instead we analyse the more realistic situation in which there is one-shot fragmentation of a shock-compressed layer, of the sort which arises in turbulent star-forming clouds. In this situation, $M_{_{\rm MIN}}$ can be smaller than $0.003 M_\odot$. Our analysis is more stringent than the standard one in that (a) it requires fragments to have condensation timescales shorter than all competing mass scales, and (b) it takes into acount that a fragment grows by accretion whilst it is condensing out, and has to radiate away the energy dissipated in the associated accretion shock. It also accords with the recent detection, in young star clusters, of free-floating star-like objects having masses as low as $0.003 {\rm M}_\odot$.

We explore the possibility that, in the vicinity of an OB star, a prestellar core which would otherwise have formed an intermediate or low-mass star may form a free-floating brown dwarf or planetary-mass object, because the outer layers of the core are eroded by the ionizing radiation from the OB star before they can accrete onto the protostar at the centre of the core. The masses of objects formed in this way are given approximately by $\sim 0.010 M_\odot (a_{\rm I} / 0.3 {\rm km} {\rm s}^{-1})^6 (\dot{\cal N}_{\rm Lyc} / 10^{50} {\rm s}^{-1})^{-1/3} (n_{\rm 0} / 10^3 {\rm cm}^{-3})^{-1/3} $, where $a_{\rm I}$ is the isothermal sound speed in the neutral gas of the core, $\dot{\cal N}_{\rm Lyc}$ is the rate of emission of Lyman continuum photons from the OB star (or stars), and $n_{\rm 0}$ is the number-density of protons in the HII region surrounding the core. We conclude that the formation of low-mass objects by this mechanism should be quite routine, because the mechanism operates over a wide range of conditions ($10^{50} {\rm s}^{-1} \la \dot{\cal N}_{\rm Lyc} \la 10^{52} {\rm s}^{-1}$, $10 {\rm cm}^{-3} \la n_0 \la 10^5 {\rm cm}^{-3}$, $0.2 {\rm km} {\rm s}^{-1} \la a_{\rm I} \la 0.6 {\rm km} {\rm s}^{-1}$) and is very effective. However, it is also a rather wasteful way of forming low-mass objects, in the sense that it requires a relatively massive initial core to form a single low-mass object. The effectiveness of photo-erosion also implies that that any intermediate-mass protostars which have formed in the vicinity of a group of OB stars must already have been well on the way to formation before the OB stars switched on their ionizing radiation; otherwise these protostars would have been stripped down to extremely low mass.

We examine the initial mass functions (IMFs) of stars produced by different molecular core mass functions. Simulations suggest that more massive cores produce more stars, so we propose a model in which the average number of stars formed in a core is equal to the initial number of Jeans masses in that core. Small-$N$ systems decay through dynamical interactions, ejecting low-mass stars and brown dwarfs which populate the low-mass tail of the IMF. Stars which remain in cores are able to competitively accrete more gas and become more massive. We deduce the forms of the core mass functions required to explain the IMFs of Taurus, Orion, IC 348 and NGC 2547. These core mass functions fall into two categories - one which peaks at a few $M_\odot$ to explain Taurus and NGC 2547, and one that peaks at around $0.2 M_\odot$ to explain Orion and IC 348.

We present 2D Monte Carlo radiative transfer simulations of flattened prestellar cores. We argue the importance of observing prestellar cores near the peak of their emission spectra, and we point out observable characteristic features on isophotal maps of asymmetric cores at FIR wavelengths that are indicative of the cores' density and temperature structure. These features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are needed to observe them. Finally, we discuss the importance of determining the temperatures of prestellar cores with high accuracy.

(Abridged) We explore, by means of a large ensemble of SPH simulations, how the level of turbulence affects the collapse and fragmentation of a star-forming core. All our simulated cores have the same, except that we vary (a) the initial level of turbulence (as measured by the ratio of turbulent to gravitational energy, $\alpha_{\rm turb} \equiv U_{\rm turb}/|\Omega| = 0, 0.01, 0.025, 0.05, 0.10 {\rm and} 0.25$) and (b), for fixed $\alpha_{\rm turb}$, the details of the initial turbulent velocity field (so as to obtain good statistics). A low level of turbulence ($\alpha_{\rm turb} \sim 0.05$) suffices to produce multiple systems. As $\alpha_{\rm turb}$ is increased, the number of objects formed and the companion frequency both increase. The mass function is bimodal, with a flat low-mass segment representing single objects ejected from the core before they can accrete much, and a Gaussian high-mass segment representing objects which because they remain in the core grow by accretion and tend to pair up in multiple systems.

We present 2D Monte Carlo radiative transfer simulations of prestellar cores. We consider two types of asymmetry: disk-like asymmetry, in which the core is denser towards the equatorial plane than towards the poles; and axial asymmetry, in which the core is denser towards the south pole than the north pole. We limit our treatment to cores with mild asymmetries, which are exposed directly to the interstellar radiation field or are embedded inside molecular clouds. The isophotal maps of a core depend strongly on the viewing angle. Maps at wavelengths longer than the peak of the SED (e.g. 850 micron) essentially trace the column-density. Thus, for instance, cores with disk-like asymmetry appear elongated when mapped at 850 micron from close to the equatorial plane. However, at wavelengths near the peak of the SED (e.g. 200 micron), the emissivity is more strongly dependent on the temperature, and therefore, at particular viewing angles, there are characteristic features which reflect a more complicated convolution of the density and temperature fields within the core. These characteristic features are on scales 1/5 to 1/3 of the overall core size, and so high resolution observations are needed to observe them. They are also weaker if the core is embedded in a molecular cloud (because the range of temperature within the core is then smaller), and so high sensitivity is needed to detect them. Herschel, to be launched in 2007, will in principle provide the necessary resolution and sensitivity at 170 to 250 micron.

We review a range of stastistical methods for analyzing the structures of star clusters, and derive a new measure ${\cal Q}$ which both quantifies, and distinguishes between, a (relatively smooth) large-scale radial density gradient and multi-scale (fractal) sub-clustering. Q is derived from the normalised correlation length and the normalised edge length of the minimal spanning tree for each cluster.

In comparison with other well studied star formation regions, Taurus is unusual in several respects. (i) Its stellar initial mass function (IMF) peaks at relatively high mass ($\sim 0.8 M_\odot$), but contains very few stars much more massive than $1 M_\odot$, and is relatively deficient in brown dwarfs. (ii) It has a higher binary fraction, particularly at large separations. (iii) Its core mass function is strongly peaked at a few $M_\odot$, and the cores have extended envelopes and relatively low levels of turbulence. We present here the results of an ensemble of hydrodynamic simulations which suggest that the unusual stellar IMF in Taurus is a direct consequence of the unusual properties of the cores there. By following the collapse and fragmentation of cores having properties typical of Taurus, we find that roughly 50% of the objects formed in a core, predominantly the low-mass ones, are ejected from the core to form a population of low-mass stars and brown dwarfs with a flat mass function. The remaining objects form multiple systems within the core, accreting until their masses approach $1 M_\odot$; this produces a population of intermediate-mass stars whose mass function peaks at $\sim 0.8 M_\odot$. Together these two populations reproduce the IMF in Taurus very well. This demonstrates, for the first time, a direct causal link between the core mass function and the stellar IMF in a star formation region.

We investigate numerically and semi-analytically the collapse of low-mass, rotating prestellar cores. Initially, the cores are in approximate equilibrium with low rotation (the initial ratio of thermal to gravitational energy is $\alpha_0 \simeq 0.5$, and the initial ratio of rotational to gravitational energy is $\beta_0 = 0.02 {\rm to} 0.05$). They are then subjected to a steady increase in external pressure. Fragmentation is promoted -- in the sense that more protostars are formed -- both by more rapid compression, and by higher rotation (larger $\beta_0$). In general, the large-scale collapse is non-homologous, and follows the pattern described in Paper I for non-rotating clouds, viz. a compression wave is driven into the cloud, thereby increasing the density and the inflow velocity. The effects of rotation become important at the centre, where the material with low angular momentum forms a central primary protostar (CPP), whilst the material with higher angular momentum forms an accretion disc around the CPP. More rapid compression drives a stronger compression wave and delivers material more rapidly into the outer parts of the disc.

We simulate the dynamics of fractal star clusters, in order to investigate the evolution of substructure in recently formed clusters. The velocity dispersion is found to be the key parameter determining the survival of substructure. In clusters with a low initial velocity dispersion, the ensuing collapse of the cluster tends to erase substructure, although some substructure may persist beyond the collapse phase. In clusters with virial ratios of 0.5 or higher, initial density substructure survives for several crossing times, in virtually all cases. Even an initially homogeneous cluster can develop substructure, if it is born with coherent velocity dispersion. These results suggest that the simple initial conditions used for many sophisticated $N$-body simulations could be missing a very important and dramatic phase of star cluster evolution.

We present the results of an ensemble of simulations of the collapse and fragmentation of dense star-forming cores. We show that even with very low levels of turbulence the outcome is usually a binary, or higher-order multiple, system. We take as the initial conditions for these simulations a typical low-mass core, based on the average properties of a large sample of observed cores. All the simulated cores start with a mass of $M = 5.4 M_{\odot}$, a flattened central density profile, a ratio of thermal to gravitational energy $\alpha_{\rm therm} = 0.45$ and a ratio of turbulent to gravitational energy $\alpha_{\rm turb} = 0.05 $. Even this low level of turbulence is sufficient to produce multiple star formation in 80% of the cores; the mean number of stars and brown dwarfs formed from a single core is 4.55, and the maximum is 10. At the outset, the cores have no large-scale rotation. The only difference between each individual simulation is the detailed structure of the turbulent velocity field. The multiple systems formed in the simulations have properties consistent with observed multiple systems. Dynamical evolution tends preferentially to eject lower mass stars and brown dwarves whilst hardening the remaining binaries so that the median semi-major axis of binaries formed is $\sim 30$ au. Ejected objects are usually single low-mass stars and brown dwarfs, yielding a strong correlation between mass and multiplicity. Our simulations suggest a natural mechanism for forming binary stars that does not require large-scale rotation, capture, or large amounts of turbulence.

We present millimetre-line imaging of the Galactic HII region Sh 104. We show that it is surrounded by a ring of molecular gas and dust. Four large molecular condensations are regularly spaced around the ring. These condensations are themselves fragmented and contain several massive dense cores. A deeply embedded cluster is observed in the near IR towards the largest condensation. It contains at least one massive star ionizing an ultra-compact HII region. The Sh 104 region is a good illustration of the `collect and collapse' model for star formation triggered by the expansion of an HII region.

We evaluate the intrinsic three dimensional shapes of molecular cores, by analysing their projected shapes. We use the recent catalogue of molecular line observations of Jijina et al. and model the data by the method originally devised for elliptical galaxies. Our analysis broadly supports the conclusion of Jones et al. that molecular cores are better represented by triaxial intrinsic shapes (ellipsoids) than biaxial intrinsic shapes (spheroids). However, we find that the best fit to all of the data is obtained with more extreme axial ratios ($1:0.8:0.4$) than those derived by Jones et al. More surprisingly, we find that starless cores have more extreme axial ratios than protostellar cores -- starless cores appear more `flattened'. This is the opposite of what would be expected from modeling the freefall collapse of triaxial ellipsoids. The collapse of starless cores would be expected to proceed most swiftly along the shortest axis - as has been predicted by every modeller since Zel'dovich - which should produce more flattened cores around protostars, the opposite of what is seen.

We implement a Monte Carlo radiative transfer method, that uses a large number of monochromatic luminosity packets to represent the radiation transported through a system. These packets are injected into the system and interact stochastically with it. We test our code against various benchmark calculations and determine the number of packets required to obtain accurate results under different circumstances. We then use this method to study cores that are directly exposed to the interstellar radiation field (non-embedded cores) and find similar results with previous studies. We also explore a large number of models of cores that are embedded in the centre of a molecular cloud. Our study indicates that the temperature profiles in embedded cores are less steep than those in non-embedded cores. Deeply embedded cores (ambient cloud with visual extinction larger than 15-25) are almost isothermal at around 7-8 K. The temperature inside cores surrounded by an ambient cloud of even moderate thickness (Av~5) is less than 12 K, which is lower than previous studies have assumed. Thus, previous mass calculations of embedded cores (for example in the rho Ophiuchi protocluster), based on mm continuum observations, may underestimate core masses by up to a factor of 2. Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 microns, where the core is quite distinct from the background. We also predict that very sensitive observations (~1-3 MJy/sr) at 170-200 microns can be used to estimate how deeply a core is embedded in its parent molecular cloud. The upcoming HERSCHEL mission (ESA, 2007) will, in principle, be able to detect these features and test our models.

We use our Monte Carlo radiative transfer code to study non-embedded prestellar cores and cores that are embedded at the centre of a molecular cloud. Our study indicates that the temperature inside embedded cores is lower than in isolated non-embedded cores, and generally less than 12 K, even when the cores are surrounded by an ambient cloud of small visual extinction (Av~5). Our study shows that the best wavelength region to observe embedded cores is between 400 and 500 microns, where the core is quite distinct from the background. We also predict that very sensitive observations (~1-3 MJy/sr) at 170-200 microns can be used to estimate how deeply a core is embedded in its parent molecular cloud. Finally, we present preliminary results of asymmetric models of non-embedded cores.

We present numerical simulations of the evolution of low-mass, isothermal, molecular cores which are subjected to an increase in external pressure $P\xt$. If $P\xt$ increases very slowly, the core approaches instability quite quasistatically. However, for larger (but still quite modest) $dP\xt/dt$ a compression wave is driven into the core, thereby triggering collapse from the outside in. If collapse of a core is induced by increasing $P\xt$, this has a number of interesting consequences. (i) The density profile is approximately flat in the centre during the prestellar phase (i.e. before the compression wave converges on the centre creating the central protostar). (ii) During the prestellar phase there are (subsonic) inward velocities in the outer layers of the core, whilst the inner parts are still approximately at rest. (iii) There is an initial short phase of rapid accretion (notionally the Class 0 phase), followed by a longer phase of slower accretion (the Class I phase). All these features accord well with observation, but are at variance with the predictions of the standard theory of star formation based on the inside-out collapse of a singular isothermal sphere. We note that the setting up of a coherent inward velocity field appears to be a generic feature of compression waves; and we speculate that interactions and interference between such velocity fields may play a crucial r\^ole in initiating the fragmentation of cores and the genesis of multiple star systems.

We describe and demonstrate a method for increasing the resolution locally in a Smoothed Particle Hydrodynamic (SPH) simulation, by splitting particles. We show that in simulations of self-gravitating collapse (of the sort which are presumed to occur in star formation) the method is stable, and affords great savings in computer time and memory. When applied to the standard Boss & Bodenheimer test -- which has been shown to depend critically on fulfilment of the Jeans Condition -- the results are comparable both with those obtained using Adaptive Mesh Refinement, and with those obtained using a standard high-resolution SPH simulation, but they are achieved with considerably less computational resource. Further development and testing is required before the method can safely be applied to more general flows.

In order to study the capabilities of SPH in self-gravitating hydrodynamical problems, we have performed a series of three-dimensional numerical simulations of the collapse of a rotating homogeneous cloud, varying the number of particles and the artificial viscosity prescription.

We present hydrodynamical simulations of the formation, structure and evolution of photoionized columns, with parameters based on those observed in the Eagle Nebula. On the basis of these simulations we argue that there is no unequivocal evidence that the dense neutral clumps at heads of the columns were cores in the pre-existing molecular cloud. In our simulations, a variety of initial conditions leads to the formation and maintenance of near-equilibrium columns. Therefore, it is likely that narrow columns will often occur in regions with large-scale inhomogeneities, but that observations of such columns can tell us little about the processes by which they formed. The manner in which the columns in our simulations develop suggests that their evolution may result in extended sequences of radiation-induced star formation.

We propose a new analytic model for the initial conditions of protostellar collapse in relatively isolated regions of star formation. The model is non-magnetic, and is based on a Plummer-like radial density profile as its initial condition. It fits: the observed density profiles of pre-stellar cores and Class 0 protostars; recent observations in pre-stellar cores of roughly constant contraction velocities over a wide range of radii; and the lifetimes and accretion rates derived for Class 0 and Class I protostars. However, the model is very simple, having in effect only 2 free parameters, and so should provide a useful framework for interpreting observations of pre-stellar cores and protostars, and for calculations of radiation transport and time-dependent chemistry. As an example, we model the pre-stellar core L1544.

We analyse and compare the clustering of young stars in Chamaeleon I and Taurus. We compute the mean surface-density of companion stars \bar{N} as a function of angular displacement \theta from each star. We then fit \bar{N}(\theta) with two simultaneous power laws. For Chamaeleon I, the exponents of the power laws are 1.97 and 0.28, with the elbow at ~0.011 degrees. For Taurus, we obtain 2.02 and 0.87, with the elbow at ~0.013 deg. For both star clusters the observational data make quite large systematic excursions from the best fitting curve in the binary regime. These excursions may be attributable to evolutionary effects of the types discussed recently by Nakajima et al. and Bate et al. In the clustering regime the data conform to the best fitting curve very well. We also calculate the box-dimensions for the two star clusters. However, the limited dynamic range makes these estimates simply descriptors of the large-scale clustering, and not admissible evidence for fractality. We also propose two algorithms for objectively generating maps of constant stellar surface-density in young star clusters. Such maps are useful for comparison with molecular-line and dust-continuum maps of star-forming clouds, and with the results of numerical simulations. These maps retain information which is suppressed in the evaluation of \bar{N}(\theta). Analyzing such maps can be used to discriminate between fractal structure and single-level clustering, and to determine the degree of central condensation in small-N clusters. Algorithm I uses a universal smoothing length, and therefore has a restricted dynamic range, but it is implicitly normalized. Algorithm II uses a local smoothing length, which gives it much greater dynamic range, but it has to be normalized explicitly.

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. If protostars are formed individually within a cluster before falling together and interacting, there should be no preferred orientation for such interactions. As star formation within clusters is believed to be coeval, it is probable that during interactions, both protostars possess massive, extended discs. We have used an SPH code to carry out a series of simulations of non-colpanar disc-disc interactions. We find that non-coplanar interactions trigger gravitational instabilities in the discs, which may then fragment to form new companions to the existing stars. (This is different from coplanar interactions, in which most of the new companion stars form after material in the discs has been swept up into a shock layer, and this then fragments.) The original stars may also capture each other, leading to the formation of a small-N cluster. If every star undergoes a randomly oriented disc-disc interaction, then the outcome will be the birth of many new stars. Approximately two-thirds of the stars will end up in multiple systems.

It appears that most stars are born in clusters, and that at birth most stars have circumstellar discs which are comparable in size to the separations between the stars. Interactions between neighbouring stars and discs are therefore likely to play a key role in determining disc lifetimes, stellar masses, and the separations and eccentricities of binary orbits. Such interactions may also cause fragmentation of the discs, thereby triggering the formation of additional stars. We have carried out a series of simulations of disc-star interactions using an SPH code which treats self-gravity, hydrodynamic and viscous forces. We find that interactions between discs and stars provide a mechanism for removing energy from, or adding energy to, the orbits of the stars, and for truncating the discs. However, capture during such encounters is unlikely to be an important binary formation mechanism. A more significant consequence of such encounters is that they can trigger fragmentation of the disc, via tidally and compressionally induced gravitational instabilities, leading to the formation of additional stars. When the disc-spins and stellar orbits are randomly oriented, encounters lead to the formation of new companions to the original star in 20% of encounters. If most encounters are prograde and coplanar, as suggested by simulations of dynamically-triggered star formation, then new companions are formed in approximately 50% of encounters.

It is expected that an average protostar will undergo at least one impulsive interaction with a neighbouring protostar whilst a large fraction of its mass is still in a massive, extended disc. Such interactions must have a significant impact upon the evolution of the protostars and their discs. We have carried out a series of simulations of coplanar encounters between two stars, each possessing a massive circumstellar disc, using an SPH code that models gravitational, hydrodynamic and viscous forces. We find that during a coplanar encounter, disc material is swept up into a shock layer between the two interacting stars, and the layer then fragments to produce new protostellar condensations. The truncated remains of the discs may subsequently fragment; and the outer regions of the discs may be thrown off to form circumbinary disc-like structures around the stars. Thus coplanar disc-disc encounters lead efficiently to the formation of multiple star systems and small-N clusters.

A Lagrangian, particle-based numerical method (tree code gravity plus smoothed particle hydrodynamics) was used to simulate clump-clump collisions occurring within GMCs. The collisions formed shock-compressed layers, out of which condensed approximately co-planar protostellar discs of 7-60 solar masses and 500-1000AU radius. Binary and multiple systems were the usual final state. Lower mass objects were also produced, but commonly underwent disruption or merger. Such objects occasionally survived by being ejected via a three-body slingshot event resulting from an encounter with a binary system. Varying the impact parameter, b, altered the processes by which the protostellar systems formed. At low b a single central disc formed initially, and was then spun-up by an accretion flow, causing it to produce secondaries via rotational instabilities. At mid b the shocked layer w hich formed initially broke up into fragments, and discs were then formed via fragment merger. At large b single objects formed within the compressed leading edge of each clump. These became unbound from each other as b was increased further. The effect of changing numerical factors was examined by : (i) colliding clumps which had been re-oriented before the collision (thus altering the initial particle noise), and (ii) by quadrupling the number of particles in each clump (thus increasing the resolution of the simulation). Both changes were found to affect the small-scale details of a collision, but leave the large scale morphology largely unaltered. It was concluded that clump-clump collisions provide a natural mechanism by which multiple protostellar systems may form.

We suggest that the abrupt switch, from hierarchical clustering on scales larger than 0.04 pc, to binary (and occasionally higher multiple) systems on smaller scales, which Larson has deduced from his analysis of the grouping of pre-Main-Sequence stars in Taurus, arises because pre-protostellar gas becomes thermally coupled to dust at sufficiently high densities. The resulting change from gas-cooling by molecular lines at low densities to gas-cooling by dust at high densities enables the matter to radiate much more efficiently, and hence to undergo dynamical fragmentation. We derive the domain where gas-cooling by dust facilitates dynamical fragmentation. Low-mass (i.e. solar mass) clumps - those supported mainly by thermal pressure - can probably access this domain spontaneously, albeit rather quasistatically, provided they exist in a region where external perturbations are few and far between. More massive clumps probably require an impulsive external perturbation, for instance a supersonic collision with another clump, in order for the gas to reach sufficiently high density to couple thermally to the dust. Impulsive external perturbations should promote fragmentation, by generating highly non-line ar substructures which can then be amplified by gravity during the subsequent collapse.