In the episodic accretion scenario, a large fraction of the protostellar mass accretes during repeated and large bursts of accretion. Since outbursts on protostars are typically identified at specific wavelengths, interpreting these outbursts requires converting this change in flux to a change in total luminosity. The Class I young stellar object EC 53 in the Serpens Main cloud has undergone repeated increases in brightness at 850 $\mu$m that are likely caused by bursts of accretion. In this study, we perform two- and three-dimensional continuum radiative transfer modeling to quantify the internal luminosity rise in EC 53 that corresponds to the factor of $\sim$1.5 enhancement in flux at 850 $\mu$m. We model the spectral energy distribution and radial intensity profile in both the quiescent and outburst phases. The internal luminosity in the outburst phase is $\sim 3.3$ times brighter than the luminosity in the quiescent phase. The radial intensity profile analysis demonstrates that the detected sub-mm flux variation of EC 53 comes from the heated envelope by the accretion burst. We also find that the role of external heating of the EC 53 envelope by the interstellar radiation field is insignificant.

Context: Around 30 per cent of the observed exoplanets that orbit M dwarf stars are gas giants that are more massive than Jupiter. These planets are prime candidates for formation by disc instability. Aims: We want to determine the conditions for disc fragmentation around M dwarfs and the properties of the planets that are formed by disc instability. Methods: We performed hydrodynamic simulations of M dwarf protostellar discs in order to determine the minimum disc mass required for gravitational fragmentation to occur. Different stellar masses, disc radii, and metallicities were considered. The mass of each protostellar disc was steadily increased until the disc fragmented and a protoplanet was formed. Results: We find that a disc-to-star mass ratio between $\sim 0.3$ and $\sim 0.6$ is required for fragmentation to happen. The minimum mass at which a disc fragments increases with the stellar mass and the disc size. Metallicity does not significantly affect the minimum disc fragmentation mass but high metallicity may suppress fragmentation. Protoplanets form quickly (within a few thousand years) at distances around $\sim50$ AU from the host star, and they are initially very hot; their centres have temperatures similar to the ones expected at the accretion shocks around planets formed by core accretion (up to 12,000K). The final properties of these planets (e.g. mass and orbital radius) are determined through long-term disc-planet or planet-planet interactions. Conclusions: Disc instability is a plausible way to form gas giant planets around M dwarfs provided that discs have at least 30% the mass of their host stars during the initial stages of their formation. Future observations of massive M dwarf discs or planets around very young M dwarfs are required to establish the importance of disc instability for planet formation around low-mass stars.

In the early stages of star formation, a protostar is deeply embedded in an optically thick envelope such that it is not directly observable. Variations in the protostellar accretion rate, however, will cause luminosity changes that are reprocessed by the surrounding envelope and are observable at submillimeter wavelengths. We searched for submillimeter flux variability toward 12 Planck Galactic Cold Clumps detected by the James Clerk Maxwell Telescope (JCMT)-SCUBA-2 Continuum Observations of Pre-protostellar Evolution (SCOPE) survey. These observations were conducted at 850 um using the JCMT/SCUBA-2. Each field was observed three times over about 14 months between 2016 April and 2017 June. We applied a relative flux calibration and achieved a calibration uncertainty of ~ 3.6% on average. We identified 136 clumps across 12 fields and detected four sources with flux variations of ~ 30%. For three of these sources, the variations appear to be primarily due to large-scale contamination, leaving one plausible candidate. The flux change of the candidate may be associated with low- or intermediate-mass star formation assuming a distance of 1.5 kpc, although we cannot completely rule out the possibility that it is a random deviation. Further studies with dedicated monitoring would provide a better understanding of the detailed relationship between submillimeter flux and accretion rate variabilities while enhancing the search for variability in star-forming clumps farther away than the Gould Belt.

Accretion onto protostars may occur in sharp bursts. Accretion bursts during the embedded phase of young protostars are probably most intense, but can only be inferred indirectly through long-wavelength observations. We perform radiative transfer calculations for young stellar objects (YSOs) formed in hydrodynamic simulations to predict the long wavelength, sub-mm and mm, flux responses to episodic accretion events, taking into account heating from the young protostar and from the interstellar radiation field. We find that the flux increase due to episodic accretion events is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of ~570 increase in the luminosity of the young protostar leads to a flux increase of a factor of 47 at 250 micron but only a factor of 10 at 1.3 mm. Heating from the interstellar radiation field may reduce further the flux increase observed at longer wavelengths. We find that during FU Ori-type outbursts the bolometric temperature and luminosity may incorrectly classify a source as a more evolved YSO, due to a larger fraction of the radiation of the object being emitted at shorter wavelengths

Young stars exhibit variability due to changes in the gas accretion rate onto them, an effect that should be quite significant in the early stages of their formation. As protostars are embedded within their natal cloud, this variability may only be inferred through long wavelength observations. We perform radiative transfer simulations of young stellar objects (YSOs) formed in hydrodynamical simulations, varying the structure and luminosity properties in order to estimate the long-wavelength, sub-mm and mm, variations of their flux. We find that the flux increase due to an outburst event depends on the protostellar structure and is more prominent at sub-mm wavelengths than at mm wavelengths; e.g. a factor of 40 increase in the luminosity of the young protostar leads to a flux increase of a factor of 10 at 250 micron but only a factor of 2.5 at 1.3 mm. We find that the interstellar radiation field dilutes the flux increase but that this effect may be avoided if resolution permits the monitoring of the inner regions of a YSO, where the heating is primarily due to protostellar radiation. We also confirm that the bolometric temperature and luminosity of outbursting protostars may result in an incorrect classification of their evolutionary stage.

We present the observational evidence of a pseudo-disc around the proto-brown dwarf Mayrit 1701117, the driving source of the large-scale HH~1165 jet. Our analysis is based on ALMA $^{12}$CO (2-1) line and 1.37 mm continuum observations at an angular resolution of $\sim$0.4$^{\prime\prime}$. The pseudo-disc is a bright feature in the CO position-velocity diagram (PVD), elongated in a direction perpendicular to the jet axis, with a total (gas+dust) mass of $\sim$0.02 M$_{\odot}$, size of 165-192 AU, and a velocity spread of $\pm$2 km s$^{-1}$. The large velocity gradient is a combination of infalling and rotational motions, indicating a contribution from a pseudo-disc and an unresolved inner Keplerian disc. There is weak emission detected in the H$_{2}$CO (3-2) and N$_{2}$D$^{+}$ (3-2) lines. H$_{2}$CO emission likely probes the inner Keplerian disc where CO is expected to be frozen, while N$_{2}$D$^{+}$ possibly originates from an enhanced clump at the outer edge of the pseudo-disc. We have considered various models (core collapse, disc fragmentation, circum-binary disc) that can fit both the observed CO spectrum and the position-velocity offsets. The observed morphology, velocity structure, and the physical dimensions of the pseudo-disc are consistent with the predictions from the core collapse simulations for brown dwarf formation. From the best model fit, we can constrain the age of the proto-brown dwarf system to be $\sim$30,000-40,000 yr. A comparison of the H$_{2}$ column density derived from the CO line and 1.37 mm continuum emission indicates that only about 2% of the CO is depleted from the gas phase.

Modern astronomy has finally been able to observe protoplanetary disks in reasonable resolution and detail, unveiling the processes happening during planet formation. These observed processes are understood under the framework of disk-planet interaction, a process studied analytically and modeled numerically for over 40 years. Long a theoreticians' game, the wealth of observational data has been allowing for increasingly stringent tests of the theoretical models. Modeling efforts are crucial to support the interpretation of direct imaging analyses, not just for potential detections but also to put meaningful upper limits on mass accretion rates and other physical quantities in current and future large-scale surveys. This white paper addresses the questions of what efforts on the computational side are required in the next decade to advance our theoretical understanding, explain the observational data, and guide new observations. We identified the nature of accretion, ab initio planet formation, early evolution, and circumplanetary disks as major fields of interest in computational planet formation. We recommend that modelers relax the approximations of alpha-viscosity and isothermal equations of state, on the grounds that these models use flawed assumptions, even if they give good visual qualitative agreement with observations. We similarly recommend that population synthesis move away from 1D hydrodynamics. The computational resources to reach these goals should be developed during the next decade, through improvements in algorithms and the hardware for hybrid CPU/GPU clusters. Coupled with high angular resolution and great line sensitivity in ground based interferometers, ELTs and JWST, these advances in computational efforts should allow for large strides in the field in the next decade.

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.

Gas clumps formed within massive gravitationally unstable circumstellar discs are potential seeds of gas giant planets, brown dwarfs and companion stars. Simulations show that competition between three processes -- migration, gas accretion and tidal disruption -- establishes what grows from a given seed. Here we investigate the robustness of numerical modelling of clump migration and accretion with the codes PHANTOM, GADGET, SPHINX, SEREN, GIZMO-MFM, SPHNG and FARGO. The test problem comprises a clump embedded in a massive disc at an initial separation of 120 AU. There is a general qualitative agreement between the codes, but the quantitative agreement in the planet migration rate ranges from $\sim 10$% to $\sim 50$%, depending on the numerical setup. We find that the artificial viscosity treatment and the sink particle prescription may account for much of the differences between the codes. In order to understand the wider implications of our work, we also attempt to reproduce the planet evolution tracks from our hydrodynamical simulations with prescriptions from three previous population synthesis studies. We find that the disagreement amongst the population synthesis models is far greater than that between our hydrodynamical simulations. The results of our code comparison project are therefore encouraging in that uncertainties in the given problem are probably dominated by the physics not yet included in the codes rather than by how hydrodynamics is modelled in them.

Radiative transfer is an important component of hydrodynamic simulations as it determines the thermal properties of a physical system. It is especially important in cases where heating and cooling regulate significant processes, such as in the collapse of molecular clouds, the development of gravitational instabilities in protostellar discs, disc-planet interactions, and planet migration. We compare two approximate radiative transfer methods which indirectly estimate optical depths within hydrodynamic simulations using two different metrics: (i) the gravitational potential and density of the gas (Stamatellos et al.), and (ii) the pressure scale-height (Lombardi et al.). We find that both methods are accurate for spherical configurations e.g. in collapsing molecular clouds and within clumps that form in protostellar discs. However, the pressure scale-height approach is more accurate in protostellar discs (low and high-mass discs, discs with spiral features, discs with embedded planets). We also investigate the $\beta$-cooling approximation which is commonly used when simulating protostellar discs, and in which the cooling time is proportional to the orbital period of the gas. We demonstrate that the use of a constant $\beta$ cannot capture the wide range of spatial and temporal variations of cooling in protostellar discs, which may affect the development of gravitational instabilities, planet migration, planet mass growth, and the orbital properties of planets.

Gas giant planets may form early-on during the evolution of protostellar discs, while these are relatively massive. We study how Jupiter-mass planet-seeds (termed protoplanets) evolve in massive, but gravitationally stable (Q>1.5), discs using radiative hydrodynamic simulations. We find that the protoplanet initially migrates inwards rapidly, until it opens up a gap in the disc. Thereafter, it either continues to migrate inwards on a much longer timescale or starts migrating outwards. Outward migration occurs when the protoplanet resides within a gap with gravitationally unstable edges, as a high fraction of the accreted gas is high angular momentum gas from outside the protoplanet's orbit. The effect of radiative heating from the protoplanet is critical in determining the direction of the migration and the eccentricity of the protoplanet. Gap opening is facilitated by efficient cooling that may not be captured by the commonly used beta-cooling approximation. The protoplanet initially accretes at a high rate (1e-3Mj/yr), and its accretion luminosity could be a few tenths of the host star's luminosity, making the protoplanet easily observable (albeit only for a short time). Due to the high gas accretion rate, the protoplanet generally grows above the deuterium-burning mass-limit. Protoplanet radiative feedback reduces its mass growth so that its final mass is near the brown dwarf-planet boundary. The fate of a young planet-seed is diverse and could vary from a gas giant planet on a circular orbit at a few AU from the central star to a brown dwarf on an eccentric, wide orbit.

Most protostars have luminosities that are fainter than expected from steady accretion over the protostellar lifetime. The solution to this problem may lie in episodic mass accretion -- prolonged periods of very low accretion punctuated by short bursts of rapid accretion. However, the timescale and amplitude for variability at the protostellar phase is almost entirely unconstrained. In "A JCMT/SCUBA-2 Transient Survey of Protostars in Nearby Star Forming Regions", we are monitoring monthly with SCUBA-2 the sub-mm emission in eight fields within nearby (<500 pc) star forming regions to measure the accretion variability of protostars. The total survey area of ~1.6 sq.deg. includes ~105 peaks with peaks brighter than 0.5 Jy/beam (43 associated with embedded protostars or disks) and 237 peaks of 0.125-0.5 Jy/beam (50 with embedded protostars or disks). Each field has enough bright peaks for flux calibration relative to other peaks in the same field, which improves upon the nominal flux calibration uncertainties of sub-mm observations to reach a precision of ~2-3% rms, and also provides quantified confidence in any measured variability. The timescales and amplitudes of any sub-mm variation will then be converted into variations in accretion rate and subsequently used to infer the physical causes of the variability. This survey is the first dedicated survey for sub-mm variability and complements other transient surveys at optical and near-IR wavelengths, which are not sensitive to accretion variability of deeply embedded protostars.

Young protostellar discs provide the initial conditions for planet formation. The properties of these discs may be different from those of late-phase (T Tauri) discs due to continuing infall from the envelope and protostellar variability resulting from irregular gas accretion. We use a set of hydrodynamic simulations to determine the structure of discs forming in collapsing molecular clouds. We examine how radiative feedback from the host protostar affects the disc properties by examining three regimes: without radiative feedback, with continuous radiative feedback and with episodic feedback, similar to FU Ori-type outbursts. We find that the radial surface density and temperature profiles vary significantly as the disc accretes gas from the infalling envelope. These profiles are sensitive to the presence of spiral structure, induced by gravitational instabilities, and the radiative feedback provided by the protostar, especially in the case when the feedback is episodic. We also investigate whether mass estimates from position-velocity (PV) diagrams are accurate for early-phase discs. We find that the protostellar system mass (i.e. the mass of the protostar and its disc) is underestimated by up to 20%, due to the impact of an enhanced radial pressure gradient on the gas. The mass of early-phase discs is a significant fraction of the mass of the protostar, so position-velocity diagrams cannot accurately provide the mass of the protostar alone. The enhanced radial pressure gradient expected in young discs may lead to an increased rate of dust depletion due to gas drag, and therefore to a reduced dust-to-gas ratio.

Protostellar discs may become massive enough to fragment producing secondary low-mass objects: planets, brown dwarfs and low-mass stars. We study the effect of radiative feedback from such newly-formed secondary objects using radiative hydrodynamic simulations. We compare the results of simulations without any radiative feedback from secondary objects with those where two types of radiative feedback are considered: (i) continuous, and (ii) episodic. We find that: (i) continuous radiative feedback stabilizes the disc and suppresses further fragmentation, reducing the number secondary objects formed; (ii) episodic feedback from secondary objects heats and stabilises the disc when the outburst occurs, but shortly after the outburst stops, the disc becomes unstable and fragments again. However, fewer secondary objects are formed compared to the the case without radiative feedback. We also find that the mass growth of secondary objects is mildly suppressed due to the effect of their radiative feedback. However, their mass growth also depends on where they form in the disc and on their subsequent interactions, such that their final masses are not drastically different from the case without radiative feedback. We find that the masses of secondary objects formed by disc fragmentation are from a few M$_{\rm J}$ to a few 0.1 M$_{\odot}$. Planets formed by fragmentation tend to be ejected from the disc. We conclude that planetary-mass objects on wide orbits (wide-orbit planets) are unlikely to form by disc fragmentation. Nevertheless, disc fragmentation may be a significant source of free-floating planets and brown dwarfs.

The origin of very low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects at the low-mass end of the initial mass function is not yet fully understood. Gravitational fragmentation of circumstellar discs provides a possible mechanism for the formation of such low-mass objects. The kinematic and binary properties of very low-mass objects formed through disc fragmentation at early times (< 10 Myr) were discussed in Li et al. (2015). In this paper we extend the analysis by following the long-term evolution of disc-fragmented systems, up to an age of 10 Gyr, covering the ages of the stellar and substellar population in the Galactic field. We find that the systems continue to decay, although the rates at which companions escape or collide with each other are substantially lower than during the first 10 Myr, and that dynamical evolution is limited beyond 1 Gyr. By t = 10 Gyr, about one third of the host stars is single, and more than half have only one companion left. Most of the other systems have two companions left that orbit their host star in widely separated orbits. A small fraction of companions have formed binaries that orbit the host star in a hierarchical triple configuration. The majority of such double companion systems have internal orbits that are retrograde with respect to their orbits around their host stars. Our simulations allow a comparison between the predicted outcomes of disc-fragmentation with the observed low-mass hydrogen-burning stars, brown dwarfs, and planetary-mass objects in the Solar neighborhood. Imaging and radial velocity surveys for faint binary companions among nearby stars are necessary for verification or rejection for the formation mechanism proposed in this paper.

Planets form in the discs of gas and dust that surround young stars. It is not known whether gas giant planets on wide orbits form the same way as Jupiter or by fragmentation of gravitationally unstable discs. Here we show that a giant planet, which has formed in the outer regions of a protostellar disc, initially migrates fast towards the central star (migration timescale ~10,000 yr) while accreting gas from the disc. However, in contrast with previous studies, we find that the planet eventually opens up a gap in the disc and the migration is essentially halted. At the same time, accretion-powered radiative feedback from the planet, significantly limits its mass growth, keeping it within the planetary mass regime (i.e. below the deuterium burning limit) at least for the initial stages of disc evolution. Giant planets may therefore be able to survive on wide orbits despite their initial fast inward migration, shaping the environment in which terrestrial planets that may harbour life form.

Theory and simulations suggest that it is possible to form low-mass hydrogen-burning stars, brown dwarfs and planetary-mass objects via disc fragmentation. As disc fragmentation results in the formation of several bodies at comparable distances to the host star, their orbits are generally unstable. Here, we study the dynamical evolution of these objects. We set up the initial conditions based on the outcomes of the SPH simulations of Stamatellos & Whitworth, and for comparison we also study the evolution of systems resulting from lower-mass fragmenting discs. We refer to these two sets of simulations as set 1 and set 2. At 10 Myr, approximately half of the host stars have one companion left, and approximately 22% (set 1) to 9.8% (set 2) of the host stars are single. Systems with multiple secondaries in relatively stable configurations are common (about 30% and 44%, respectively). The majority of the companions are ejected within 1 Myr with velocities mostly below 5 km/s, with some runaway escapers with velocities over 30 km/s. About 6% (set 1) and 2% (set 2) of the companions pair up into very low-mass binary systems. The majority of these pairs escape as very low-mass binaries, while others remain bound to the host star in hierarchical configurations (often with retrograde inner orbits). Physical collisions with the host star (0.43 and 0.18 events per host star for set 1 and set 2) and between companions (0.08 and 0.04 events per host star for set 1 and set 2) are relatively common and their frequency increases with increasing disc mass. Our study predicts observable properties of very low-mass binaries, low-mass hierarchical systems, the brown dwarf desert, and free-floating brown dwarfs and planetary-mass objects in and near young stellar groupings, which can be used to distinguish between different formation scenarios of very low-mass stars, brown dwarfs and planetary-mass objects.

Direct imaging searches have revealed many very low-mass objects, including a small number of planetary mass objects, as wide-orbit companions to young stars. The formation mechanism of these objects remains uncertain. In this paper we present the predictions of the disc fragmentation model regarding the properties of the discs around such low-mass objects. We find that the discs around objects that have formed by fragmentation in discs hosted by Sun-like stars (referred to as 'parent' discs and 'parent' stars) are more massive than expected from the ${M}_{\rm disc}-M_*$ relation (which is derived for stars with masses $M_*>0.2 {\rm M}_{\odot}$). Accordingly, the accretion rates onto these objects are also higher than expected from the $\dot{M}_*-M_*$ relation. Moreover there is no significant correlation between the mass of the brown dwarf or planet with the mass of its disc nor with the accretion rate from the disc onto it. The discs around objects that form by disc fragmentation have larger than expected masses as they accrete gas from the disc of their parent star during the first few kyr after they form. The amount of gas that they accrete and therefore their mass depend on how they move in their parent disc and how they interact with it. Observations of disc masses and accretion rates onto very low-mass objects are consistent with the predictions of the disc fragmentation model. Future observations (e.g. by ALMA) of disc masses and accretion rates onto substellar objects that have even lower masses (young planets and young, low-mass brown dwarfs), where the scaling relations predicted by the disc fragmentation model diverge significantly from the corresponding relations established for higher-mass stars, will test the predictions of this model.

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.

The density and temperature structures of dense cores in the L1495 cloud of the Taurus star-forming region are investigated using Herschel SPIRE and PACS images in the 70 $\mu$m, 160 $\mu$m, 250 $\mu$m, 350 $\mu$m and 500 $\mu$m continuum bands. A sample consisting of 20 cores, selected using spectral and spatial criteria, is analysed using a new maximum likelihood technique, COREFIT, which takes full account of the instrumental point spread functions. We obtain central dust temperatures, $T_0$, in the range 6-12 K and find that, in the majority of cases, the radial density falloff at large radial distances is consistent with the $r^{-2}$ variation expected for Bonnor-Ebert spheres. Two of our cores exhibit a significantly steeper falloff, however, and since both appear to be gravitationally unstable, such behaviour may have implications for collapse models. We find a strong negative correlation between $T_0$ and peak column density, as expected if the dust is heated predominantly by the interstellar radiation field. At the temperatures we estimate for the core centres, carbon-bearing molecules freeze out as ice mantles on dust grains, and this behaviour is supported here by the lack of correspondence between our estimated core locations and the previously-published positions of H$^{13}$CO$^+$ peaks. On this basis, our observations suggest a sublimation-zone radius typically $\sim 10^4$ AU. Comparison with previously-published N$_2$H$^+$ data at 8400 AU resolution, however, shows no evidence for N$_2$H$^+$ depletion at that resolution.

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.

Utilizing multi-wavelength dust emission maps acquired with $Herschel$, we reconstruct local volume density and dust temperature profiles for the prestellar cores B68 and L1689B using inverse-Abel transform based technique. We present intrinsic radial dust temperature profiles of starless cores directly from dust continuum emission maps disentangling the effect of temperature variations along the line of sight which was previously limited to the radiative transfer calculations. The reconstructed dust temperature profiles show a significant drop in core center, a flat inner part, and a rising outward trend until the background cloud temperature is reached. The central beam-averaged dust temperatures obtained for B68 and L1689B are 9.3 $\pm$ 0.5 K and 9.8 $\pm$0.5 K, respectively, which are lower than the temperatures of 11.3 K and 11.6 K obtained from direct SED fitting. The best mass estimates derived by integrating the volume density profiles of B68 and L1689B are 1.6 M_sol and 11 M_sol, respectively. Comparing our results for B68 with the near-infrared extinction studies, we find that the dust opacity law adopted by the HGBS project, $\kappa_{\lambda} =0.1(\lambda/300 \mu m)^{-2}$, agrees to within 50% with the dust extinction constraints

I discuss the role that disc fragmentation plays in the formation of gas giant and terrestrial planets, and how this relates to the formation of brown dwarfs and low-mass stars, and ultimately to the process of star formation. Protostellar discs may fragment, if they are massive enough and can cool fast enough, but most of the objects that form by fragmentation are brown dwarfs. It may be possible that planets also form, if the mass growth of a proto-fragment is stopped (e.g. if this fragment is ejected from the disc), or suppressed and even reversed (e.g by tidal stripping). I will discuss if it is possible to distinguish whether a planet has formed by disc fragmentation or core accretion, and mention of a few examples of observed exoplanets that are suggestive of formation by disc fragmentation .

It is estimated that ~60% of all stars (including brown dwarfs) have masses below 0.2Msun. Currently, there is no consensus on how these objects form. I will briefly review the four main theories for the formation of low-mass objects: turbulent fragmentation, ejection of protostellar embryos, disc fragmentation, and photo-erosion of prestellar cores. I will focus on the disc fragmentation theory and discuss how it addresses critical observational constraints, i.e. the low-mass initial mass function, the brown dwarf desert, and the binary statistics of low-mass stars and brown dwarfs. I will examine whether observations may be used to distinguish between different formation mechanisms, and give a few examples of systems that strongly favour a specific formation scenario. Finally, I will argue that it is likely that all mechanisms may play a role in low-mass star and brown dwarf formation.

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.

In an attempt to study whether the formation of brown dwarfs (BDs) takes place as a scaled-down version of low-mass stars, we conducted IRAM30m/MAMBO-II observations at 1.2 mm in a sample of 12 proto-BD candidates selected from Spitzer/IRAC data in the B213-L1495 clouds in Taurus. Subsequent observations with the CSO at 350 micron, VLA at 3.6 and 6 cm, and IRAM30m/EMIR in the 12CO(1-0), 13CO(1-0), and N2H+(1-0) transitions were carried out toward the two most promising Spitzer/IRAC source(s), J042118 and J041757. J042118 is associated with a compact (<10 arcsec or <1400 AU) and faint source at 350 micron, while J041757 is associated with a partially resolved (~16 arcsec or ~2000 AU) and stronger source emitting at centimetre wavelengths with a flat spectral index. The corresponding masses of the dust condensations are ~1 and ~5 Mjup for J042118 and J041757, respectively. In addition, about 40 arcsec to the northeast of J041757 we detect a strong and extended submillimetre source, J041757-NE, which is not associated with NIR/FIR emission down to our detection limits, but is clearly detected in 13CO and N2H+ at ~7 km/s, and for which we estimated a total mass of ~100 Mjup, close to the mass required to be gravitationally bound. In summary, our observational strategy has allowed us to find in B213-L1495 two proto-BD candidates and one pre-substellar core candidate, whose properties seem to be consistent with a scaled-down version of low-mass stars.

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 revised our analysis of the SPIRE observations of 2MASSW J1207334-393254 (2M1207). Recent PACS observations show a bright source located ~25" east of 2M1207. There are issues in terms of the detection/non-detection of the bright source when comparing the Spitzer, WISE, and PACS observations. It is apparently inconsistent, perhaps due to variability or low signal-to-noise of the data. We have looked into the possible misidentification of the target, and have revised the measured SPIRE fluxes and the disc parameters for 2M1207. We have also reviewed which among the various formation mechanisms of this system would still be valid.

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.

We have investigated the brown dwarf (BD) and stellar disc fractions in the Upper Scorpius OB Association (USco) and compared them with several other young regions. We have compiled the most complete sample of of all spectroscopically confirmed BDs in USco, and have made use of the WISE catalog to identify the disc candidates. We report on the discovery of 12 new BD discs in USco, with spectral type (SpT) between M6 and M8.5. The WISE colors for the new discs are similar to the primordial (transition) discs earlier detected in USco. Combining with previous surveys, we find the lowest inner disc fractions (~20-25%) for a wide range in stellar masses (~0.01-4.0 Msun) in the USco association. The low disc fractions for high-mass stars in USco (and the other clusters) are consistent with an evolutionary decline in inner disc frequency with age. However, BD disc fractions are higher than those for the stars in 1-3 Myr clusters, but very low in the ~5 Myr old USco. Also, primordial BD discs are still visible in the ~10 Myr TW Hydrae association, whereas the higher mass stars have all transitioned to the debris stage by this age. The disc frequencies for BDs do not show any dependence on the stellar density or the BD/star number ratio in a cluster. We suggest that the large differences in the observed BD disc fractions between regions may well be due to different BD formation mechanisms and therefore different initial disc fractions/properties. We also present a WISE SED classification scheme, based on the Ks and WISE bands of 3.4-12micron, to distinguish between the Class I/II and the Class III sequence. Our work includes a comparison of the sensitivities of WISE and Spitzer disc surveys. We estimate that WISE can be incomplete for discs at SpT later than M8 in distant clusters such as SOri. WISE should be able to recover the M8-M9 discs in the nearby young 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.

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.

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.

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

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 present the first results from the science demonstration phase for the Hi-GAL survey, the Herschel key-project that will map the inner Galactic Plane of the Milky Way in 5 bands. We outline our data reduction strategy and present some science highlights on the two observed 2{\deg} x 2{\deg} tiles approximately centered at l=30{\deg} and l=59{\deg}. The two regions are extremely rich in intense and highly structured extended emission which shows a widespread organization in filaments. Source SEDs can be built for hundreds of objects in the two fields, and physical parameters can be extracted, for a good fraction of them where the distance could be estimated. The compact sources (which we will call 'cores' in the following) are found for the most part to be associated with the filaments, and the relationship to the local beam-averaged column density of the filament itself shows that a core seems to appear when a threshold around A_V of about 1 is exceeded for the regions in the l=59{\deg} field; a A_V value between 5 and 10 is found for the l=30{\deg} field, likely due to the relatively larger distances of the sources. This outlines an exciting scenario where diffuse clouds first collapse into filaments, which later fragment to cores where the column density has reached a critical level. In spite of core L/M ratios being well in excess of a few for many sources, we find core surface densities between 0.03 and 0.5 g cm-2. Our results are in good agreement with recent MHD numerical simulations of filaments forming from large-scale converging flows.

We test the operation of two methods for selective application of Artificial Viscosity (AV) in SPH simulations of Keplerian Accretion Disks, using a ring spreading test to quantify effective viscosity, and a correlation coefficient technique to measure the formation of unwanted prograde alignments of particles. Neither the Balsara Switch nor Time Dependent Viscosity work effectively, as they leave AV active in areas of smooth shearing flow, and do not eliminate the accumulation of alignments of particles in the prograde direction. The effect of both switches is periodic, the periodicity dependent on radius and unaffected by the density of particles. We demonstrate that a very simple algorithm activates AV only when truly convergent flow is detected and reduces the unwanted formation of prograde alignments. The new switch works by testing whether all the neighbours of a particle are in Keplerian orbit around the same point, rather than calculating the divergence of the velocity field, which is very strongly affected by Poisson noise in the positions of the SPH particles.

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 present observations of L1155 and L1148 in the Cepheus molecular cloud, taken using the FIS instrument on the Akari satellite. We compare these data to submillimetre data taken using the SCUBA camera on the JCMT, and far-infrared data taken with the ISOPHOT camera on board the ISO satellite. All of the data show a relation between the position of the peak of emission and the wavelength for the core of L1155. We interpret this as a temperature gradient. We fit modified blackbody curves to the spectral energy distributions at two positions in the core and see that the central core in L1155 (L1155C) is approximately 2 degrees warmer at one edge than it is in the centre. We consider a number of possible heating sources and conclude that the A6V star BD+67 1263 is the most likely candidate. This star is at a distance of 0.7 pc from the front of L1155C in the plane of the sky. We carry out radiative transfer modelling of the L1155C core including the effects from the nearby star. We find that we can generate a good fit to the observed data at all wavelengths, and demonstrate that the different morphologies of the core at different wavelengths can be explained by the observed 2 degree temperature gradient. The L1148 core exhibits a similar morphology to that of L1155C, and the data are also consistent with a temperature gradient across the core. In this case, the most likely heating source is the star BD197053. Our findings illustrate very clearly that the apparent observed morphology of a pre-stellar core can be highly dependent on the wavelength of the observation, and that temperature gradients must be taken into account before converting images into column density distributions. This is important to note when interpreting Akari and Spitzer data and will also be significant for Herschel data.

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)

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.

We present continuum data from the Submillimetre Common-User Bolometer Array (SCUBA) on the James Clerk Maxwell Telescope (JCMT), and the Mid-Infrared Photometer for Spitzer (MIPS) on the Spitzer Space Telescope, at submillimetre and infrared wavelengths respectively. We study the Taurus molecular cloud 1 (TMC1), and in particular the region of the Taurus Molecular Ring (TMR). In the continuum data we see no real evidence for a ring, but rather we see one side of it only, appearing as a filament. We name the filament `the bull's tail'. The filament is seen in emission at 850, 450 and 160um, and in absorption at 70um. We compare the data with archive data from the Infra-Red Astronomical Satellite (IRAS) at 12, 25, 60, 100um, in which the filament is also seen in absorption. We find that the emission from the filament consists of two components: a narrow, cold (~8K), central core; and a broader, slightly warmer (~12K), shoulder of emission. We use a radiative transfer code to model the filament's appearance, either in emission or absorption, simultaneously at each of the different wavelengths. Our best fit model uses a Plummer-like density profile and a homogeneous interstellar dust grain population. Unlike previous work on a similar, but different filament in Taurus, we require no grain coagulation to explain our data.

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.

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.

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

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