We present simulations of a 500 pc$^2$ region, containing gas of mass 4 $\times$ 10$^6$ M$_\odot$, extracted from an entire spiral galaxy simulation, scaled up in resolution, including photoionising feedback from stars of mass > 18 M$_\odot$. Our region is evolved for 10 Myr and shows clustered star formation along the arm generating $\approx$ 5000 cluster sink particles $\approx$ 5% of which contain at least one of the $\approx$ 4000 stars of mass > 18 M$_\odot$. Photoionisation has a noticeable effect on the gas in the region, producing ionised cavities and leading to dense features at the edge of the HII regions. Compared to the no-feedback case, photoionisation produces a larger total mass of clouds and clumps, with around twice as many such objects, which are individually smaller and more broken up. After this we see a rapid decrease in the total mass in clouds and the number of clouds. Unlike studies of isolated clouds, our simulations follow the long range effects of ionisation, with some already-dense gas becoming compressed from multiple sides by neighbouring HII regions. This causes star formation that is both accelerated and partially displaced throughout the spiral arm with up to 30% of our cluster sink particle mass forming at distances > 5 pc from sites of sink formation in the absence of feedback. At later times, the star formation rate decreases to below that of the no-feedback case.

In young multiple stellar systems the gravitational influence of the stars shapes the circumstellar disk, controlling accretion and the material available for planet formation. Our observations of the triple system GW Orionis provide direct evidence for disk tearing, where the gravitational torques due to misaligned disk/orbit planes warp the disk and break the inner disk into precessing rings. We image an eccentric ring that is misaligned with respect to the orbital planes and outer disk, and casts shadows on the strongly warped intermediate disk. Our ring/warp geometry constraints and the fully characterized perturber orbits make the system a potential Rosetta Stone for studying disk hydrodynamics. The ring might offer suitable conditions for planet formation, providing a mechanism for forming wide-separation planets on highly oblique orbits.

We show results from the radiation hydrodynamics (RHD) simulations of tidal disruption of a star on a parabolic orbit by a supermassive black hole (SMBH) based on a three-dimensional smoothed particle hydrodynamics code with radiative transfer. We find that such a tidally disrupted star fragment and form clumps soon after its tidal disruption. The fragmentation results from the endothermic processes of ionization and dissociation that reduce the gas pressure, leading to local gravitational collapse. Radiative cooling is less effective because the stellar debris is still highly optically thick in such an early time. Our simulations reveal that a solar-type star with a stellar density profile of n=3 disrupted by a 10^6 solar mass black hole produces $\sim20$ clumps of masses in the range of 0.1 to 12 Jupiter masses. The mass fallback rate decays with time, with pronounced spikes from early to late time. The spikes provide evidence for the clumps of the returning debris, while the clumps on the unbound debris can be potentially freely-floating planets and brown dwarfs. This ionization and dissociation induced fragmentation on a tidally disrupted star are a promising candidate mechanism to form low-mass stars to planets around an SMBH.

We present results from the first radiation non-ideal magnetohydrodynamics (MHD) simulations of low-mass star cluster formation that resolve the fragmentation process down to the opacity limit. We model 50~M$_\odot$ turbulent clouds initially threaded by a uniform magnetic field with strengths of 3, 5 10 and 20 times the critical mass-to-magnetic flux ratio, and at each strength, we model both an ideal and non-ideal (including Ohmic resistivity, ambipolar diffusion and the Hall effect) MHD cloud. Turbulence and magnetic fields shape the large-scale structure of the cloud, and similar structures form regardless of whether ideal or non-ideal MHD is employed. At high densities ($10^6 \lesssim n_{\rm H} \lesssim 10^{11}$~cm$^{-3}$), all models have a similar magnetic field strength versus density relation, suggesting that the field strength in dense cores is independent of the large-scale environment. Albeit with limited statistics, we find no evidence for the dependence of the initial mass function on the initial magnetic field strength, however, the star formation rate decreases for models with increasing initial field strengths; the exception is the strongest field case where collapse occurs primarily along field lines. Protostellar discs with radii $\gtrsim 20$~au form in all models, suggesting that disc formation is dependent on the gas turbulence rather than on magnetic field strength. We find no evidence for the magnetic braking catastrophe, and find that magnetic fields do not hinder the formation of protostellar discs.

Resolution studies of test problems set baselines and help define minimum resolution requirements, however, resolution studies must also be performed on scientific simulations to determine the effect of resolution on the specific scientific results. We perform a resolution study on the formation of a protostar by modelling the collapse of gas through 14 orders of magnitude in density. This is done using compressible radiative non-ideal magnetohydrodynamics. Our suite consists of an ideal magnetohydrodynamics (MHD) model and two non-ideal MHD models, and we test three resolutions for each model. The resulting structure of the ideal MHD model is approximately independent of resolution, although higher magnetic field strengths are realised in higher resolution models. The non-ideal MHD models are more dependent on resolution, specifically the magnetic field strength and structure. Stronger magnetic fields are realised in higher resolution models, and the evolution of detailed structures such as magnetic walls are only resolved in our highest resolution simulation. In several of the non-ideal MHD models, there is an off-set between the location of the maximum magnetic field strength and the maximum density, which is often obscured or lost at lower resolutions. Thus, understanding the effects of resolution on numerical star formation is imperative for understanding the formation of a star.

The first stable object to develop in the low-mass star formation process has long been predicted to be the first hydrostatic core (FHSC). Despite much effort, it has still yet to be definitively observed in nature. More specific observational signatures are required to enable observers to distinguish the FHSC from young, faint, but more evolved protostars. Here we present synthetic spectral line observations for CO, SO, CS and HCO$^+$ that were calculated from radiation (magneto)hydrodynamical models, chemical modelling and Monte Carlo radiative transfer. HCO$^+$ $(1-0)$ and SO $(8_7 - 7_6)$ spectra of the FHSC show variations for observations at a low inclination which may allow a candidate FHSC to be distinguished from a more evolved object. We find that the FHSC outflow is unlikely to be detectable with ALMA, which would discount the observed sources with slow outflows that are currently identified as candidate FHSCs. We compare the results of simulated ALMA observations with observed candidate FHSCs and recommend Oph A SM1N and N6-mm as the most promising candidates to follow up.

We investigate the formation and fragmentation of discs using a suite of three-dimensional smoothed particle radiative magnetohydrodynamics simulations. Our models are initialised as 1M$_\odot$ rotating Bonnor-Ebert spheres that are threaded with a uniform magnetic field. We examine the effect of including ideal and non-ideal magnetic fields, the orientation and strength of the magnetic field, and the initial rotational rate. We follow the gravitational collapse and early evolution of each system until the final classification of the protostellar disc can be determined. Of our 105 models, 41 fragment, 21 form a spiral structure but do not fragment, and another 12 form smooth discs. Fragmentation is more likely to occur for faster initial rotation rates and weaker magnetic fields. For stronger magnetic field strengths, the inclusion of non-ideal MHD promotes disc formation, and several of these models fragment, whereas their ideal MHD counterparts do not. For the models that fragment, there is no correlation between our parameters and where or when the fragmentation occurs. Bipolar outflows are launched in only 17 models, and these models have strong magnetic fields that are initially parallel to the rotation axis. Counter-rotating envelopes form in four slowly-rotating, strong-field models -- including one ideal MHD model -- indicating they form only in a small fraction of the parameter space investigated.

Substellar multiplicity is a key outcome of the formation process. The biggest challenge for the next decade will be to distinguish between the formation history, environmental conditions, and dynamical evolution leading to the least massive brown dwarfs and the most massive planets at the tail ends of their mass functions. In this white paper, we advocate for a comprehensive characterization of both the statistical distributions of the population of ultracool dwarf multiple systems and the fundamental properties of their individual components as a function of age. A space-based precision astrometry mission in near-infrared wavelengths would provide the necessary measurements to identify and characterize age-calibrated populations of multiple systems.

We report the statistical properties of stars and brown dwarfs obtained from four radiation hydrodynamical simulations of star cluster formation, the metallicities of which span a range from 1/100 to 3 times the solar value. Unlike previous similar investigations of the effects of metallicity on stellar properties, these new calculations treat dust and gas temperatures separately and include a thermochemical model of the diffuse interstellar medium. The more advanced treatment of the interstellar medium gives rise to very different gas and dust temperature distributions in the four calculations, with lower metallicities generally resulting in higher temperatures and a delay in the onset of star formation. Despite this, once star formation begins, all four calculations produce stars at similar rates and many of the statistical properties of their stellar populations are difficult to distinguish from each other and from those of observed stellar systems. We do find, however, that the greater cooling rates at high gas densities due to the lower opacities at low metallicities increase the fragmentation on small spatial scales (disc, filament, and core fragmentation). This produces an anti-correlation between the close binary fraction of low-mass stars and metallicity similar to that which is observed, and an increase in the fraction of protostellar mergers at low metallicities. There are also indications that at lower metallicity close binaries may have lower mass ratios and the abundance of brown dwarfs to stars may increase slightly. However, these latter two effects are quite weak and need to be confirmed with larger samples.

Are the kG-strength magnetic fields observed in young stars a fossil field left over from their formation or are they generated by a dynamo? We use radiation non-ideal magnetohydrodynamics simulations of the gravitational collapse of a rotating, magnetized molecular cloud core over 17 orders of magnitude in density, past the first hydrostatic core to the formation of the second, stellar core, to examine the fossil field hypothesis. Whereas in previous work we found that magnetic fields in excess of 10 kG can be implanted in stars at birth, this assumed ideal magnetohydrodynamics (MHD), i.e. that the gas is coupled to the magnetic field. Here we present non-ideal MHD calculations which include Ohmic resistivity, ambipolar diffusion and the Hall effect. For realistic cosmic ray ionization rates, we find that magnetic field strengths of $\lesssim$ kG are implanted in the stellar core at birth, ruling out a strong fossil field. While these results remain sensitive to resolution, they cautiously provide evidence against a fossil field origin for stellar magnetic fields, suggesting instead that magnetic fields in stars originate in a dynamo process.

We demonstrate the formation of gravitationally unstable discs in magnetized molecular cloud cores with initial mass-to-flux ratios of 5 times the critical value, effectively solving the magnetic braking catastrophe. We model the gravitational collapse through to the formation of the stellar core, using Ohmic resistivity, ambipolar diffusion and the Hall effect and using the canonical cosmic ray ionization rate of $\zeta_\text{cr} = 10^{-17}$ s$^{-1}$. When the magnetic field and rotation axis are initially aligned, a $\lesssim1$~au disc forms after the first core phase, whereas when they are anti-aligned, a gravitationally-unstable 25~au disc forms during the first core phase. The aligned model launches a 3~km~s$^{-1}$ first core outflow, while the anti-aligned model launches only a weak $\lesssim 0.3$~km~s$^{-1}$ first core outflow. Qualitatively, we find that models with $\zeta_\text{cr} = 10^{-17}$ s$^{-1}$ are similar to purely hydrodynamical models if the rotation axis and magnetic field are initially anti-aligned, whereas they are qualitatively similar to ideal magnetohydrodynamical models if initially aligned.

We present a new method for including radiative feedback from sink particles in smoothed particle hydrodynamics simulations of low-mass star formation, and investigate its effects on the formation of small stellar groups. We find that including radiative feedback from sink particles suppresses fragmentation even further than calculations that only include radiative transfer within the gas. This reduces the star-formation rate following the formation of the initial protostars, leading to fewer objects being produced and a lower total stellar mass. The luminosities of sink particles vary due to changes in the accretion rate driven by the dynamics of the cluster gas, leading to different luminosities for protostars of similar mass. Including feedback from sinks also raises the median stellar mass. The median masses of the groups are higher than typically observed values. This may be due to the lack of dynamical interactions and ejections in small groups of protostars compared to those that occur in richer groups. We also find that the temperature distributions in our calculations are in qualitative agreement with recent observations of protostellar heating in Galactic star-forming regions.

We investigate the dependence of stellar properties on the initial mean density of the molecular cloud in which stellar clusters form using radiation hydrodynamical simulations that resolve the opacity limit for fragmentation. We have simulated the formation of three star clusters from the gravitational collapse of molecular clouds whose densities vary by a factor of a hundred. As with previous calculations including radiative feedback, we find that the dependence of the characteristic stellar mass, $M_{\rm c}$, on the initial mean density of the cloud, $\rho$, is weaker than the dependence of the thermal Jeans mass. However, unlike previous calculations, which found no statistically significant variation in the median mass with density, we find a weak dependence approximately of the form $M_{\rm c} \propto \rho^{-1/5}$. The distributions of properties of multiple systems do not vary significantly between the calculations. We compare our results to the result of observational surveys of star-forming regions, and suggest that the similarities between the properties of our lowest density calculation and the nearby Taurus-Auriga region indicate that the apparent excess of solar-type stars observed may be due to the region's low density.

The pace and pattern of star formation leading to rich young stellar clusters is quite uncertain. In this context, we analyze the spatial distribution of ages within 19 young (median t<3 Myr on the Siess et al. (2000) timescale), morphologically simple, isolated, and relatively rich stellar clusters. Our analysis is based on young stellar object samples from the MYStIX and SFiNCs surveys, and a new estimator of pre-main sequence (PMS) stellar ages, AgeJX, derived from X-ray and near-infrared photometric data. Median cluster ages are computed within four annular subregions of the clusters. We confirm and extend the earlier result of Getman et al. (2014): 80% percent of the clusters show age trends where stars in cluster cores are younger than in outer regions. Our cluster stacking analyses establish the existence of an age gradient to high statistical significance in several ways. Time scales vary with the choice of PMS evolutionary model; the inferred median age gradient across the studied clusters ranges from 0.75 Myr/pc to 1.5 Myr/pc. The empirical finding reported in the present study -- late or continuing formation of stars in the cores of star clusters with older stars dispersed in the outer regions -- has a strong foundation with other observational studies and with the astrophysical models like the global hierarchical collapse model of Vazquez-Semadeni et al. (2017).

Photometric detections of dust circumstellar disks around pre-main sequence (PMS) stars, coupled with estimates of stellar ages, provide constraints on the time available for planet formation. Most previous studies on disk longevity, starting with Haisch, Lada & Lada (2001), use star samples from PMS clusters but do not consider datasets with homogeneous photometric sensitivities and/or ages placed on a uniform timescale. Here we conduct the largest study to date of the longevity of inner dust disks using X-ray and 1--8 micrometre infrared photometry from the MYStIX and SFiNCs projects for 69 young clusters in 32 nearby star-forming regions with ages t<=5 Myr. Cluster ages are derived by combining the empirical AgeJX method with PMS evolutionary models, which treat dynamo-generated magnetic fields in different ways. Leveraging X-ray data to identify disk-free objects, we impose similar stellar mass sensitivity limits for disk-bearing and disk-free YSOs while extending the analysis to stellar masses as low as M~0.1 Mo. We find that the disk longevity estimates are strongly affected by the choice of PMS evolutionary model. Assuming a disk fraction of 100% at zero age, the inferred disk half-life changes significantly, from t1/2 ~ 1.3--2 Myr to t1/2 ~ 3.5 Myr when switching from non-magnetic to magnetic PMS models. In addition, we find no statistically significant evidence that disk fraction varies with stellar mass within the first few Myr of life for stars with masses <2 Mo, but our samples may not be complete for more massive stars. The effects of initial disk fraction and star-forming environment are also explored.

The SFiNCs (Star Formation in Nearby Clouds) project is an X-ray/infrared study of the young stellar populations in 22 star forming regions with distances <=1 kpc designed to extend our earlier MYStIX survey of more distant clusters. Our central goal is to give empirical constraints on cluster formation mechanisms. Using parametric mixture models applied homogeneously to the catalog of SFiNCs young stars, we identify 52 SFiNCs clusters and 19 unclustered stellar structures. The procedure gives cluster properties including location, population, morphology, association to molecular clouds, absorption, age (AgeJX), and infrared spectral energy distribution (SED) slope. Absorption, SED slope, and AgeJX are age indicators. SFiNCs clusters are examined individually, and collectively with MYStIX clusters, to give the following results. (1) SFiNCs is dominated by smaller, younger, and more heavily obscured clusters than MYStIX. (2) SFiNCs cloud-associated clusters have the high ellipticities aligned with their host molecular filaments indicating morphology inherited from their parental clouds. (3) The effect of cluster expansion is evident from the radius-age, radius-absorption, and radius-SED correlations. Core radii increase dramatically from ~0.08 to ~0.9 pc over the age range 1--3.5 Myr. Inferred gas removal timescales are longer than 1 Myr. (4) Rich, spatially distributed stellar populations are present in SFiNCs clouds representing early generations of star formation. An Appendix compares the performance of the mixture models and nonparametric Minimum Spanning Tree to identify clusters. This work is a foundation for future SFiNCs/MYStIX studies including disk longevity, age gradients, and dynamical modeling.

We present the results of eighteen magnetohydrodynamical calculations of the collapse of a molecular cloud core to form a protostar. Some calculations include radiative transfer in the flux-limited diffusion approximation while others employ a barotropic equation of state. We cover a wide parameter space, with mass-to-flux ratios ranging from $\mu = 5$ to $20$; initial turbulent amplitudes ranging from a laminar calculation (i.e. where the Mach number, $\mathscr{M} = 0$) to transonic $\mathscr{M} = 1$; and initial rotation rates from $\beta_\mathrm{rot} = 0.005$ to $0.02$. We first show that using a radiative transfer scheme produces warmer pseudo-discs than the barotropic equation of state, making them more stable. We then `shake' the core by increasing the initial turbulent velocity field, and find that at all three mass-to-flux ratios transonic cores are weakly bound and do not produce pseudo-discs; $\mathscr{M} = 0.3$ cores produce very disrupted discs; and $\mathscr{M} = 0.1$ cores produce discs broadly comparable to a laminar core. In our previous paper (arXiv:1701.08741), we showed that a pseudo-disc coupled with sufficient magnetic field is necessary to form a bipolar outflow. Here we show that only weakly turbulent cores exhibit collimated jets. We finally take the $\mathscr{M} = 1.0$, $\mu = 5$ core and `stir' it by increasing the initial angular momentum, finding that once the degree of rotational energy exceeds the turbulent energy in the core the disc returns, with a corresponding (though slower), outflow. These conclusions place constraints on the initial mixtures of rotation and turbulence in molecular cloud cores which are conducive to the formation of bipolar outflows early in the star formation process.

What cosmic ray ionisation rate is required such that a non-ideal magnetohydrodynamics (MHD) simulation of a collapsing molecular cloud will follow the same evolutionary path as an ideal MHD simulation or as a purely hydrodynamics simulation? To investigate this question, we perform three-dimensional smoothed particle non-ideal magnetohydrodynamics simulations of the gravitational collapse of rotating, one solar mass, magnetised molecular cloud cores, that include Ohmic resistivity, ambipolar diffusion, and the Hall effect. We assume a uniform grain size of $a_\text{g} = 0.1\mu$m, and our free parameter is the cosmic ray ionisation rate, $\zeta_\text{cr}$. We evolve our models, where possible, until they have produced a first hydrostatic core. Models with $\zeta_\text{cr}\gtrsim10^{-13}$ s$^{-1}$ are indistinguishable from ideal MHD models and the evolution of the model with $\zeta_\text{cr}=10^{-14}$ s$^{-1}$ matches the evolution of the ideal MHD model within one per cent when considering maximum density, magnetic energy, and maximum magnetic field strength as a function of time; these results are independent of $a_\text{g}$. Models with very low ionisation rates ($\zeta_\text{cr}\lesssim10^{-24}$ s$^{-1}$) are required to approach hydrodynamical collapse, and even lower ionisation rates may be required for larger $a_\text{g}$. Thus, it is possible to reproduce ideal MHD and purely hydrodynamical collapses using non-ideal MHD given an appropriate cosmic ray ionisation rate. However, realistic cosmic ray ionisation rates approach neither limit, thus non-ideal MHD cannot be neglected in star formation simulations.

We present results from the first population synthesis study of protostellar discs. We analyse the evolution and properties of a large sample of protostellar discs formed in a radiation hydrodynamical simulation of star cluster formation. Due to the chaotic nature of the star formation process, we find an enormous diversity of young protostellar discs, including misaligned discs, and discs whose orientations vary with time. Star-disc interactions truncate discs and produce multiple systems. Discs may be destroyed in dynamical encounters and/or through ram-pressure stripping, but reform by later gas accretion. We quantify the distributions of disc mass and radii for protostellar ages up to $\approx 10^5$ yrs. For low-mass protostars, disc masses tend to increase with both age and protostellar mass. Disc radii range from of order ten to a few hundred au, grow in size on timescales $\le 10^4$ yr, and are smaller around lower-mass protostars. The radial surface density profiles of isolated protostellar discs are flatter than the minimum mass solar nebula model, typically scaling as $\Sigma \propto r^{-1}$. Disc to protostar mass ratios rarely exceed two, with a typical range of $M_{\rm d}/M_* = 0.1-1$ to ages $\le 10^4$ yrs and decreasing thereafter. We quantify the relative orientation angles of circumstellar discs and the orbit of bound pairs of protostars, finding a preference for alignment that strengths with decreasing separation. We also investigate how the orientations of the outer parts of discs differ from the protostellar and inner disc spins for isolated protostars and pairs.

We present results from radiation non-ideal magnetohydrodynamics (MHD) calculations that follow the collapse of rotating, magnetised, molecular cloud cores to stellar densities. These are the first such calculations to include all three non-ideal effects: ambipolar diffusion, Ohmic resistivity and the Hall effect. We employ an ionisation model in which cosmic ray ionisation dominates at low temperatures and thermal ionisation takes over at high temperatures. We explore the effects of varying the cosmic ray ionisation rate from $\zeta_\text{cr}= 10^{-10}$ to $10^{-16}$ s$^{-1}$. Models with ionisation rates $\gtrsim 10^{-12}$ s$^{-1}$ produce results that are indistinguishable from ideal MHD. Decreasing the cosmic ray ionisation rate extends the lifetime of the first hydrostatic core up to a factor of two, but the lifetimes are still substantially shorter than those obtained without magnetic fields. Outflows from the first hydrostatic core phase are launched in all models, but the outflows become broader and slower as the ionisation rate is reduced. The outflow morphology following stellar core formation is complex and strongly dependent on the cosmic ray ionisation rate. Calculations with high ionisation rates quickly produce a fast (~14km s$^{-1}$) bipolar outflow that is distinct from the first core outflow, but with the lowest ionisation rate a slower (~3-4 km s$^{-1}$) conical outflow develops gradually and seamlessly merges into the first core outflow.

The radial drift problem constitutes one of the most fundamental problems in planet formation theory, as it predicts particles to drift into the star before they are able to grow to planetesimal size. Dust-trapping vortices have been proposed as a possible solution to this problem, as they might be able to trap particles over millions of years, allowing them to grow beyond the radial drift barrier. Here, we present ALMA 0.04"-resolution imaging of the pre-transitional disk of V1247 Orionis that reveals an asymmetric ring as well as a sharply-confined crescent structure, resembling morphologies seen in theoretical models of vortex formation. The asymmetric ring (at 0.17"=54 au separation from the star) and the crescent (at 0.38"=120 au) seem smoothly connected through a one-armed spiral arm structure that has been found previously in scattered light. We propose a physical scenario with a planet orbiting at $\sim0.3$"$\approx$100 au, where the one-armed spiral arm detected in polarised light traces the accretion stream feeding the protoplanet. The dynamical influence of the planet clears the gap between the ring and the crescent and triggers two vortices that trap mm-sized particles, namely the crescent and the bright asymmetry seen in the ring. We conducted dedicated hydrodynamics simulations of a disk with an embedded planet, which results in similar spiral-arm morphologies as seen in our scattered light images. At the position of the spiral wake and the crescent we also observe $^{12}$CO (3-2) and H$^{12}$CO$^{+}$ (4-3) excess line emission, likely tracing the increased scale-height in these disk regions.

The first hydrostatic core (FHSC) is the first stable object to form in simulations of star formation. This stage has yet to be observed definitively, although several candidate FHSCs have been reported. We have produced synthetic spectral energy distributions (SEDs) from 3D hydrodynamical simulations of pre-stellar cores undergoing gravitational collapse for a variety of initial conditions. Variations in the initial rotation rate, radius and mass lead to differences in the location of the SED peak and far-infrared flux. Secondly, we attempt to fit the SEDs of five FHSC candidates from the literature and five newly identified FHSC candidates located in the Serpens South molecular cloud with simulated SEDs. The most promising FHSC candidates are fitted by a limited number of model SEDs with consistent properties, which suggests the SED can be useful for placing constraints on the age and rotation rate of the source. The sources we consider most likely to be in FHSC phase are B1-bN, CB17-MMS, Aqu-MM1 and Serpens South candidate K242. We were unable to fit SerpS-MM22, Per-Bolo 58 and Chamaeleon-MMS1 with reasonable parameters, which indicates that they are likely to be more evolved.

In numerical simulations, artificial terms are applied to the evolution equations for stability. To prove their validity, these terms are thoroughly tested in test problems where the results are well known. However, they are seldom tested in production-quality simulations at high resolution where they interact with a plethora of physical and numerical algorithms. We test three artificial resistivities in both the Orszag-Tang vortex and in a star formation simulation. From the Orszag-Tang vortex, the Price et. al. (2017) artificial resistivity is the least dissipative thus captures the density and magnetic features; in the star formation algorithm, each artificial resistivity algorithm interacts differently with the sink particle to produce various results, including gas bubbles, dense discs, and migrating sink particles. The star formation simulations suggest that it is important to rely upon physical resistivity rather than artificial resistivity for convergence.

We report results from twelve simulations of the collapse of a molecular cloud core to form one or more protostars, comprising three field strengths (mass-to-flux ratios, {\mu}, of 5, 10, and 20) and four field geometries (with values of the angle between the field and rotation axes, {\theta}, of 0{\deg}, 20{\deg}, 45{\deg}, and 90{\deg}), using a smoothed particle magnetohydrodynamics method. We find that the values of both parameters have a strong effect on the resultant protostellar system and outflows. This ranges from the formation of binary systems when {\mu} = 20 to strikingly differing outflow structures for differing values of {\theta}, in particular highly suppressed outflows when {\theta} = 90{\deg}. Misaligned magnetic fields can also produce warped pseudo-discs where the outer regions align perpendicular to the magnetic field but the innermost region re-orientates to be perpendicular to the rotation axis. We follow the collapse to sizes comparable to those of first cores and find that none of the outflow speeds exceed 8 km s$^{-1}$. These results may place constraints on both observed protostellar outflows, and also on which molecular cloud cores may eventually form either single stars and binaries: a sufficiently weak magnetic field may allow for disc fragmentation, whilst conversely the greater angular momentum transport of a strong field may inhibit disc fragmentation.

High-mass multiples might form via fragmentation of self-gravitational disks or alternative scenarios such as disk-assisted capture. However, only few observational constraints exist on the architecture and disk structure of high-mass protobinaries and their accretion properties. Here we report the discovery of a close ($57.9\pm0.2$mas=170au) high-mass protobinary, IRAS17216-3801, where our VLTI/GRAVITY+AMBER near-infrared interferometry allows us to image the circumstellar disks around the individual components with 3 milliarcsecond resolution. We estimate the component masses to $\sim20$ and $\sim18 M_{\odot}$ and find that the radial intensity profiles can be reproduced with an irradiated disk model, where the inner regions are excavated of dust, likely tracing the dust sublimation region in these disks. The circumstellar disks are strongly misaligned with respect to the binary separation vector, which indicates that the tidal forces did not have time to realign the disks, pointing towards a young dynamical age of the system. We constrain the distribution of the Br$\gamma$ and CO-emitting gas using VLTI/GRAVITY spectro-interferometry and VLT/CRIRES spectro-astrometry and find that the secondary is accreting at a higher rate than the primary. VLT/NACO imaging shows $L'$-band emission on 3-4 times larger scales than the binary separation, matching the expected dynamical truncation radius for the circumbinary disk. The IRAS17216-3801 system is $\sim3\times$ more massive and $\sim5\times$ more compact than other high-mass multiplies imaged at infrared wavelengths and the first high-mass protobinary system where circumstellar and circumbinary dust disks could be spatially resolved. This opens exciting new opportunities for studying star-disk interactions and the role of multiplicity in high-mass star formation.

The Star Formation in Nearby Clouds (SFiNCs) project is aimed at providing detailed study of the young stellar populations and star cluster formation in nearby 22 star forming regions (SFRs) for comparison with our earlier MYStIX survey of richer, more distant clusters. As a foundation for the SFiNCs science studies, here, homogeneous data analyses of the Chandra X-ray and Spitzer mid-infrared archival SFiNCs data are described, and the resulting catalogs of over 15300 X-ray and over 1630000 mid-infrared point sources are presented. On the basis of their X-ray/infrared properties and spatial distributions, nearly 8500 point sources have been identified as probable young stellar members of the SFiNCs regions. Compared to the existing X-ray/mid-infrared publications, the SFiNCs member list increases the census of YSO members by 6-200% for individual SFRs and by 40% for the merged sample of all 22 SFiNCs SFRs.

We investigate the effect of non-ideal magnetohydrodynamics (MHD) on the formation of binary stars using a suite of three-dimensional smoothed particle magnetohydrodynamics simulations of the gravitational collapse of one solar mass, rotating, perturbed molecular cloud cores. Alongside the role of Ohmic resistivity, ambipolar diffusion and the Hall effect, we also examine the effects of magnetic field strength, orientation and amplitude of the density perturbation. When modelling sub-critical cores, ideal MHD models do not collapse whereas non-ideal MHD models collapse to form single protostars. In super-critical ideal MHD models, increasing the magnetic field strength or decreasing the initial density perturbation amplitude decreases the initial binary separation. Strong magnetic fields initially perpendicular to the rotation axis suppress the formation of binaries and yield discs with magnetic fields ~10 times stronger than if the magnetic field was initially aligned with the rotation axis. When non-ideal MHD is included, the resulting discs are larger and more massive, and the binary forms on a wider orbit. Small differences in the super-critical cores caused by non-ideal MHD effects are amplified by the binary interaction near periastron. Overall, the non-ideal effects have only a small impact on binary formation and early evolution, with the initial conditions playing the dominant role.

The dynamics of dust and gas can be quite different from each other when the dust is poorly coupled to the gas. In protoplanetary discs, it is well known that this decoupling of the dust and gas can lead to diverse spatial structures and dust-to-gas ratios. In this paper, we study the dynamics of dust and gas during the earlier phase of protostellar collapse, before a protoplanetary disc is formed. We find that for dust grains with sizes < 10 micron, the dust is well coupled during the collapse of a rotating, pre-stellar core and there is little variation of the dust-to-gas ratio during the collapse. However, if larger grains are present, they may have trajectories that are very different from the gas during the collapse, leading to mid-plane settling and/or oscillations of the dust grains through the mid-plane. This may produce variations in the dust-to-gas ratio and very different distributions of large and small dust grains at the very earliest stages of star formation, if large grains are present in pre-stellar cores.

We test the hypothesis that the initial mass function (IMF) is determined by the density probability distribution function (PDF) produced by supersonic turbulence. We compare 14 simulations of star cluster formation in 50 solar mass molecular cloud cores where the initial turbulence contains either purely solenoidal or purely compressive modes, in each case resolving fragmentation to the opacity limit to determine the resultant IMF. We find statistically indistinguishable IMFs between the two sets of calculations, despite a factor of two difference in the star formation rate and in the standard deviation of $\log(\rho)$. This suggests that the density PDF, while determining the star formation rate, is not the primary driver of the IMF.

The Planet Formation Imager (PFI) project aims to provide a strong scientific vision for ground-based optical astronomy beyond the upcoming generation of Extremely Large Telescopes. We make the case that a breakthrough in angular resolution imaging capabilities is required in order to unravel the processes involved in planet formation. PFI will be optimised to provide a complete census of the protoplanet population at all stellocentric radii and over the age range from 0.1 to about 100 Myr. Within this age period, planetary systems undergo dramatic changes and the final architecture of planetary systems is determined. Our goal is to study the planetary birth on the natural spatial scale where the material is assembled, which is the "Hill Sphere" of the forming planet, and to characterise the protoplanetary cores by measuring their masses and physical properties. Our science working group has investigated the observational characteristics of these young protoplanets as well as the migration mechanisms that might alter the system architecture. We simulated the imprints that the planets leave in the disk and study how PFI could revolutionise areas ranging from exoplanet to extragalactic science. In this contribution we outline the key science drivers of PFI and discuss the requirements that will guide the technology choices, the site selection, and potential science/technology tradeoffs.

We study the evolution of galactic magnetic fields using 3D smoothed particle magnetohydrodynamics (SPMHD) simulations of galaxies with an imposed spiral potential. We consider the appearance of reversals of the field, and amplification of the field. We find magnetic field reversals occur when the velocity jump across the spiral shock is above $\approx$20km s$^{-1}$, occurring where the velocity change is highest, typically at the inner Lindblad resonance (ILR) in our models. Reversals also occur at corotation, where the direction of the velocity field reverses in the co-rotating frame of a spiral arm. They occur earlier with a stronger amplitude spiral potential, and later or not at all with weaker or no spiral arms. The presence of a reversal at a radii of around 4--6 kpc in our fiducial model is consistent with a reversal identified in the Milky Way, though we caution that alternative Galaxy models could give a similar reversal. We find that relatively high resolution, a few million particles in SPMHD, is required to produce consistent behaviour of the magnetic field. Amplification of the magnetic field occurs in the models, and while some may be genuinely attributable to differential rotation or spiral arms, some may be a numerical artefact. We check our results using Athena, finding reversals but less amplification of the field, suggesting that some of the amplification of the field with SPMHD is numerical.

We present an updated constrained hyperbolic/parabolic divergence cleaning algorithm for smoothed particle magnetohydrodynamics (SPMHD) that remains conservative with wave cleaning speeds which vary in space and time. This is accomplished by evolving the quantity $\psi / c_h$ instead of $\psi$. Doing so allows each particle to carry an individual wave cleaning speed, $c_h$, that can evolve in time without needing an explicit prescription for how it should evolve, preventing circumstances which we demonstrate could lead to runaway energy growth related to variable wave cleaning speeds. This modification requires only a minor adjustment to the cleaning equations and is trivial to adopt in existing codes. Finally, we demonstrate that our constrained hyperbolic/parabolic divergence cleaning algorithm, run for a large number of iterations, can reduce the divergence of the field to an arbitrarily small value, achieving $\nabla \cdot B=0$ to machine precision.

Obtaining a stable magnetohydrodynamical (MHD) formalism in SPH - i.e. smoothed particle magnetohydrodynamics (SPMHD) - has proven remarkably difficult. To implement MHD requires two steps: a modification to the momentum equation and an induction equation, and both present challenges. We first provide an overview of how SPMHD is implemented, and then discuss how this implementation fails and the limitation of various corrective methods - with particular reference to the effects of particle disorder. Although there are many problems for which, with careful choice of corrective measures, good results can be obtained, we then show that, at the very limits of the state of the art, the ability to perform stable MHD calculations in SPH is curtailed by numerical issues.

We investigate whether or not the low ionisation fractions in molecular cloud cores can solve the `magnetic braking catastrophe', where magnetic fields prevent the formation of circumstellar discs around young stars. We perform three-dimensional smoothed particle non-ideal magnetohydrodynamics (MHD) simulations of the gravitational collapse of one solar mass molecular cloud cores, incorporating the effects of ambipolar diffusion, Ohmic resistivity and the Hall effect alongside a self-consistent calculation of the ionisation chemistry assuming 0.1 micron grains. When including only ambipolar diffusion or Ohmic resistivity, discs do not form in the presence of strong magnetic fields, similar to the cases using ideal MHD. With the Hall effect included, disc formation depends on the direction of the magnetic field with respect to the rotation vector of the gas cloud. When the vectors are aligned, strong magnetic braking occurs and no disc is formed. When the vectors are anti-aligned, a disc with radius of 13AU can form even in strong magnetic when all three non-ideal terms are present, and a disc of 38 AU can form when only the Hall effect is present; in both cases, a counter-rotating envelope forms around the first hydrostatic core. For weaker, anti-aligned fields, the Hall effect produces massive discs comparable to those produced in the absence of magnetic fields, suggesting that planet formation via gravitational instability may depend on the sign of the magnetic field in the precursor molecular cloud core.

In an earlier letter, we reported that dust settling in protoplanetary discs may lead to a dynamical dust-gas instability that produces global toroidal vortices. In this letter, we investigate the evolution of a dusty protoplanetary disc with two different dust species (1 mm and 50 cm dust grains), under the presence of the instability. We show how toroidal vortices, triggered by the interaction of mm grains with the gas, stop the radial migration of metre-sized dust, potentially offering a natural and efficient solution to the dust migration problem.

We present an improved version of the Loren-Aguilar & Bate (2014) method to integrate the two-fluid dust/gas equations that correctly captures the limiting velocity of small grains in the presence of net differences (excluding the drag force) between the accelerations of the dust and the gas. A series of accelerated DUSTYBOX tests and a simulation of dust-settling in a protoplanetary disc are performed comparing the performance of the new and old methods. The modified method can accurately capture the correct limiting velocity while preserving all the conservation properties of the original method.

We identify a new hydrodynamical instability in protoplanetary discs that may arise due to variations in the dust-to-gas ratio and may lead to concentration of dust grains within a disc. The instability can arise due to dust settling, which produces a vertical compositional entropy gradient. The entropy gradient drives a baroclinic instability that is capable of creating toroidal gas vortices that gather dust into rings. Such dust rings are potentially observable via continuum emission of the dust or scattered light. Indeed, this instability may offer an explanation for the rings recently observed in the discs around the young stars HL Tau and TW Hya that does not rely on clearing by protoplanets. The instability may also have wider ramifications, potentially aiding dust agglomeration, altering the radial migration of larger planetesimals, and modifying angular momentum transport within a disc.

We analyze the physical properties of stellar clusters that are detected in massive star-forming regions in the MYStIX project--a comparative, multiwavelength study of young stellar clusters within 3.6 kpc that contain at least one O-type star. Tabulated properties of subclusters in these regions include physical sizes and shapes, intrinsic numbers of stars, absorptions by the molecular clouds, and median subcluster ages. Physical signs of dynamical evolution are present in the relations of these properties, including statistically significant correlations between subcluster size, central density, and age, which are likely the result of cluster expansion after gas removal. We argue that many of the subclusters identified in Paper I are gravitationally bound because their radii are significantly less than what would be expected from freely expanding clumps of stars with a typical initial stellar velocity dispersion of ~3 km/s for star-forming regions. We explore a model for cluster formation in which structurally simpler clusters are built up hierarchically through the mergers of subclusters--subcluster mergers are indicated by an inverse relation between the numbers of stars in a subcluster and their central densities (also seen as a density vs. radius relation that is less steep than would be expected from pure expansion). We discuss implications of these effects for the dynamical relaxation of young stellar clusters.

The equations of smoothed particle magnetohydrodynamics (SPMHD), even with the various corrections to instabilities so far proposed, have been observed to be unstable when a very steep density gradient is necessarily combined with a variable smoothing length formalism. Here we consider in more detail the modifications made to the SPMHD equations in LBP2015 that resolve this instability by replacing the smoothing length in the induction and anisotropic force equations with an average smoothing length term. We then explore the choice of average used and compare the effects on a test `cylinder-in-a-box' problem and the collapse of a magnetised molecular cloud core. We find that, aside from some benign numerical effects at low resolutions for the quadratic mean, the formalism is robust as to the choice of average but that in complicated models it is essential to apply the average to both equations; in particular, all four averages considered exhibit similar conservation properties. This improved formalism allows for arbitrarily small sink particles and field geometries to be explored, vastly expanding the range of astronomical problems that can be modeled using SPMHD.

We have developed a modified form of the equations of smoothed particle magnetohydrodynamics which are stable in the presence of very steep density gradients. Using this formalism, we have performed simulations of the collapse of magnetised molecular cloud cores to form protostars and drive outflows. Our stable formalism allows for smaller sink particles (< 5 AU) than used previously and the investigation of the effect of varying the angle, {\theta}, between the initial field axis and the rotation axis. The nature of the outflows depends strongly on this angle: jet-like outflows are not produced at all when {\theta} > 30{\deg}, and a collimated outflow is not sustained when {\theta} > 10{\deg}. No substantial outflows of any kind are produced when {\theta} > 60{\deg}. This may place constraints on the geometry of the magnetic field in molecular clouds where bipolar outflows are seen.

We present a method for modelling star-forming clouds that combines two different models of the thermal evolution of the interstellar medium (ISM). In the combined model, where the densities are low enough that at least some part of the spectrum is optically thin, a model of the thermodynamics of the diffuse ISM is more significant in setting the temperatures. Where the densities are high enough to be optically thick across the spectrum, a model of flux limited diffusion is more appropriate. Previous methods either model the low-density interstellar medium and ignore the thermal behaviour at high densities (e.g. inside collapsing molecular cloud cores), or model the thermal behaviour near protostars but assume a fixed background temperature (e.g. approximately 10 K) on large-scales. Our new method treats both regimes. It also captures the different thermal evolution of the gas, dust, and radiation separately. We compare our results with those from the literature, and investigate the dependence of the thermal behaviour of the gas on the various model parameters. This new method should allow us to model the ISM across a wide range of densities and, thus, develop a more complete and consistent understanding of the role of thermodynamics in the star formation process.

The arm structure of the Milky Way remains somewhat of an unknown, with observational studies hindered by our location within the Galactic disc. In the work presented here we use smoothed particle hydrodynamics (SPH) and radiative transfer to create synthetic longitude-velocity observations. Our aim is to reverse-engineer a top down map of the Galaxy by comparing synthetic longitude-velocity maps to those observed. We set up a system of N-body particles to represent the disc and bulge, allowing for dynamic creation of spiral features. Interstellar gas, and the molecular content, is evolved alongside the stellar system. A 3D-radiative transfer code is then used to compare the models to observational data. The resulting models display arm features that are a good reproduction of many of the observed emission structures of the Milky Way. These arms however are dynamic and transient, allowing for a wide range of morphologies not possible with standard density wave theory. The best fitting models are a much better match than previous work using fixed potentials. They favour a 4-armed model with a pitch angle of approximately 20 degrees, though with a pattern speed that decreases with increasing Galactic radius. Inner bars are lacking however, which appear required to fully reproduce the central molecular zone.

Among the most fascinating and hotly-debated areas in contemporary astrophysics are the means by which planetary systems are assembled from the large rotating disks of gas and dust which attend a stellar birth. Although important work has already been, and is still being done both in theory and observation, a full understanding of the physics of planet formation can only be achieved by opening observational windows able to directly witness the process in action. The key requirement is then to probe planet-forming systems at the natural spatial scales over which material is being assembled. By definition, this is the so-called Hill Sphere which delineates the region of influence of a gravitating body within its surrounding environment. The Planet Formation Imager project (PFI) has crystallized around this challenging goal: to deliver resolved images of Hill-Sphere-sized structures within candidate planet-hosting disks in the nearest star-forming regions. In this contribution we outline the primary science case of PFI. For this purpose, we briefly review our knowledge about the planet-formation process and discuss recent observational results that have been obtained on the class of transition disks. Spectro-photometric and multi-wavelength interferometric studies of these systems revealed the presence of extended gaps and complex density inhomogeneities that might be triggered by orbiting planets. We present detailed 3-D radiation-hydrodynamic simulations of disks with single and multiple embedded planets, from which we compute synthetic images at near-infrared, mid-infrared, far-infrared, and sub-millimeter wavelengths, enabling a direct comparison of the signatures that are detectable with PFI and complementary facilities such as ALMA. From these simulations, we derive some preliminary specifications that will guide the array design and technology roadmap of the facility.

A method to avoid the explicit time integration of small dust grains in the two fluid gas/dust smoothed particle hydrodynamics (SPH) approach is proposed. By assuming a very simple exponential decay model for the relative velocity between the gas and dust components, all the effective characteristics of the drag force can be reproduced. A series of tests has been performed to compare the accuracy of the method with analytical and explicit integration results. We find that the method performs well on a wide range of tests, and can provide large speed ups over explicit integration when the dust stopping time is small. We have also found that the method is much less dissipative than conventional explicit or implicit two-fluid SPH approaches when modelling dusty shocks.

We report the statistical properties of stars and brown dwarfs obtained from four radiation hydrodynamical simulations of star cluster formation that resolve masses down to the opacity limit for fragmentation. The calculations are identical except for their dust and gas opacities. Assuming dust opacity is proportional to metallicity, the calculations span a range of metallicities from 1/100 to 3 times solar, although we emphasise that changing the metallicity has other thermodynamic effects that the calculations do not capture (e.g. on the thermal coupling between gas and dust). All four calculations produce stellar populations whose statistical properties are difficult to distinguish from observed stellar systems, and we find no significant dependence of stellar properties on opacity. The mass functions and properties of multiple stellar systems are consistent with each other. However, we find protostellar mergers are more common with lower opacities. Combining the results from the three calculations with the highest opacities, we obtain a stellar population consisting of more than 500 stars and brown dwarfs. Many of the statistical properties of this population are in good agreement with those observed in our Galaxy, implying that gravity, hydrodynamics, and radiative feedback may be the primary ingredients for determining the statistical properties of low-mass stars. However, we do find indications that the calculations may be slightly too dissipative. Although further calculations will be required to understand all of the effects of metallicity on stellar properties, we conclude that stellar properties are surprisingly resilient to variations of the dust and gas opacities.

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.

The clusters of young stars in massive star-forming regions show a wide range of sizes, morphologies, and numbers of stars. Their highly subclustered structures are revealed by the MYStIX project's sample of 31,754 young stars in nearby sites of star formation (regions at distances <3.6 kpc that contain at least one O-type star.) In 17 of the regions surveyed by MYStIX, we identify subclusters of young stars using finite mixture models -- collections of isothermal ellipsoids that model individual subclusters. Maximum likelihood estimation is used to estimate the model parameters, and the Akaike Information Criterion is used to determine the number of subclusters. This procedure often successfully finds famous subclusters, such as the BN/KL complex behind the Orion Nebula Cluster and the KW-object complex in M17. A catalog of 142 subclusters is presented, with 1 to 20 subclusters per region. The subcluster core radius distribution for this sample is peaked at 0.17 pc with a standard deviation of 0.43 dex, and subcluster core radius is negatively correlated with gas/dust absorption of the stars -- a possible age effect. Based on the morphological arrangements of subclusters, we identify four classes of spatial structure: long chains of subclusters, clumpy structures, isolated clusters with a core-halo structure, and isolated clusters well fit by a single isothermal ellipsoid.

Stars do not generally form in isolation. Instead, they form in clusters, and in these clustered environments newborn stars can have profound effects on one another and on their parent gas clouds. Feedback from clustered stars is almost certainly responsible for a number of otherwise puzzling facts about star formation: that it is an inefficient process that proceeds slowly when averaged over galactic scales; that most stars disperse from their birth sites and dissolve into the galactic field over timescales $\ll 1$ Gyr; and that newborn stars follow an initial mass function (IMF) with a distinct peak in the range $0.1 - 1$ $M_\odot$, rather than an IMF dominated by brown dwarfs. In this review we summarize current observational constraints and theoretical models for the complex interplay between clustered star formation and feedback.

We presents results from Smoothed Particle Magnetohydrodynamics simulations of collapsing molecular cloud cores, and dynamo amplification of the magnetic field in the presence of Mach 10 magnetised turbulence. Our star formation simulations have produced, for the first time ever, highly collimated magnetised protostellar jets from the first hydrostatic core phase. Up to 40% of the initial core mass may be ejected through this outflow. The primary difficulty in performing these simulations is maintaining the divergence free constraint of the magnetic field, and to address this issue, we have developed a new divergence cleaning method which has allowed us to stably follow the evolution of these protostellar jets for long periods. The simulations performed of supersonic MHD turbulence are able to exponentially amplify magnetic energy by up to 10 orders of magnitude via turbulent dynamo. To reduce numerical dissipation, a new shock detection algorithm is utilised which is able to track magnetic shocks throughout a large range of magnetic field strengths.

We have performed smoothed particle radiation magnetohydrodynamics (SPRMHD) simulations of the collapse of rotating, magnetised molecular cloud cores to form protostars. The calculations follow the formation and evolution of the first hydrostatic core, the collapse to form a stellar core, the launching of outflows from both the first hydrostatic core and stellar cores, and the breakout of the stellar outflow from the remnant of the first core. We investigate the roles of magnetic fields and thermal feedback on the outflow launching process, finding that both magnetic and thermal forces contribute to the launching of the stellar outflow. We also follow the stellar cores until they grow to masses of up to 20 Jupiter-masses, and determine their properties. We find that at this early stage, before fusion begins, the stellar cores have radii of $\approx 3$ R$_\odot$ with radial entropy profiles that increase outward (i.e. are convectively stable) and minimum entropies per baryon of $s/k_{\rm B} \approx 14$ in their interiors. The structure of the stellar cores is found to be insensitive to variations in the initial magnetic field strength. With reasonably strong initial magnetic fields, accretion on to the stellar cores occurs through inspiralling magnetised pseudo-discs with negligible radiative losses, as opposed to first cores which effectively radiate away the energy liberated in the accretion shocks at their surfaces. We find that magnetic field strengths of >10 kG can be implanted in stellar cores at birth.

MYStIX (Massive Young Star-Forming Complex Study in Infrared and Xray) seeks to characterize 20 OB-dominated young clusters and their environs at distances d < 4 kpc using imaging detectors on the Chandra X-ray Observatory, Spitzer Space Telescope, and the United Kingdom InfraRed Telescope. The observational goals are to construct catalogs of star-forming complex stellar members with well-defined criteria, and maps of nebular gas (particularly of hot X-ray emitting plasma) and dust. A catalog of MYStIX Probable Complex Members (MPCMs) with several hundred OB stars and > 30, 000 low mass premain sequence is assembled. This sample and related data products will be used to seek new empirical constraints on theoretical models of cluster formation and dynamics, mass segregation, OB star formation, star formation triggering on the periphery of HII regions, the survivability of protoplanetary disks in HII regions. This paper give an introduction and overview of the project, covering the data analysis methodology and application to two star forming regions, NGC 2264 and the Trifid Nebula.

We carry out simulations of gravitationally unstable discs using a Smoothed Particle Hydrodynamics (SPH) code and a grid-based hydrodynamics code, FARGO, to understand the previous non-convergent results reported by Meru & Bate (2011a). We obtain evidence that convergence with increasing resolution occurs with both SPH and FARGO and in both cases we find that the critical cooling timescale is larger than previously thought. We show that SPH has a first-order convergence rate while FARGO converges with a second-order rate. We show that the convergence of the critical cooling timescale for fragmentation depends largely on the numerical viscosity employed in both SPH and FARGO. With SPH, particle velocity dispersion may also play a role. We show that reducing the dissipation from the numerical viscosity leads to larger values of the critical cooling time at a given resolution. For SPH, we find that the effect of the dissipation due to the numerical viscosity is somewhat larger than had previously been appreciated. In particular, we show that using a quadratic term in the SPH artificial viscosity (beta_{SPH}) that is too low appears to lead to excess dissipation in gravitationally unstable discs, which may affect any results that sensitively depend on the thermodynamics, such as disc fragmentation. We show that the two codes converge to values of the critical cooling timescale, beta_{crit} > 20 (for a ratio of specific heats of gamma=5/3), and perhaps even as large as beta_{crit} \approx 30. These are approximately 3-5 times larger than has been found by most previous studies. This is equivalent to a maximum gravitational stress that a disc can withstand without fragmenting of alpha_{GI,crit} \approx 0.013-0.02, which is much smaller than the values typically used in the literature. It is therefore easier for self-gravitating discs to fragment than has been concluded from most past studies.

We have conducted three-dimensional self-gravitating radiation hydrodynamical models of gas accretion onto high mass cores (15-33 Earth masses) over hundreds of orbits. Of these models, one case accretes more than a third of a Jupiter mass of gas, before eventually undergoing a hydrodynamic collapse. This collapse causes the density near the core to increase by more than an order of magnitude, and the outer envelope to evolve into a circumplanetary disc. A small reduction in the mass within the Hill radius (R_H) accompanies this collapse as a shock propagates outwards. This collapse leads to a new hydrostatic equilibrium for the protoplanetary envelope, at which point 97 per cent of the mass contained within the Hill radius is within the inner 0.03 R_H which had previously contained less than 40 per cent. Following this collapse the protoplanet resumes accretion at its prior rate. The net flow of mass towards this dense protoplanet is predominantly from high latitudes, whilst at the outer edge of the circumplanetary disc there is net outflow of gas along the midplane. We also find a turnover of gas deep within the bound envelope that may be caused by the establishment of convection cells.

We have performed three-dimensional two-fluid (gas-dust) hydrodynamical models of circumstellar discs with embedded protoplanets (3 - 333 M\oplu) and small solid bodies (radii 10cm to 10m). We find that high mass planets (\gtrsim Saturn mass) open sufficiently deep gaps in the gas disc such that the density maximum at the outer edge of the gap can very efficiently trap metre-sized solid bodies. This allows the accumulation of solids at the outer edge of the gap as solids from large radii spiral inwards to the trapping region. This process of accumulation occurs fastest for those bodies that spiral inwards most rapidly, typically metre-sized boulders, whilst smaller and larger objects will not migrate sufficiently rapidly in the discs lifetime to benefit from the process. Around a Jupiter mass planet we find that bound clumps of solid material, as large as several Earth masses, may form, potentially collapsing under self-gravity to form planets or planetesimals. These results are in agreement with Lyra et al. (2009), supporting their finding that the formation of a second generation of planetesimals or of terrestrial mass planets may be triggered by the presence of a high mass planet.

We have performed Smoothed Particle Magnetohydrodynamics (SPMHD) simulations demonstrating the production of collimated jets during collapse of 1 solar mass molecular cloud cores to form the `first hydrostatic core' in low mass star formation. Recently a number of candidate first core objects have been observed, including L1448 IRS2E, L1451-mm and Per Bolo 58, although it is not yet clear that these are first hydrostatic cores. Recent observations of Per Bolo 58 in particular appear to show collimated, bipolar outflows which are inconsistent with previous theoretical expectations. We show that low mass first cores can indeed produce tightly collimated jets (opening angles <~ 10 degrees) with speeds of ~2-7 km/s, consistent with some of the observed candidates. We have also demonstrated, for the first time, that such phenomena can be successfully captured in SPMHD simulations.

We report the statistical properties of stars, brown dwarfs and multiple systems obtained from the largest radiation hydrodynamical simulation of star cluster formation to date that resolves masses down to the opacity limit for fragmentation (a few Jupiter masses). The initial conditions are identical to those of previous barotropic calculations published by Bate, but this time the calculation is performed using a realistic equation of state and radiation hydrodynamics. The calculation uses sink particles to model 183 stars and brown dwarfs, including 28 binaries and 12 higher-order multiple systems, the properties of which are compared the results from observational surveys. We find that the radiation hydrodynamical/sink particle simulation reproduces many observed stellar properties very well. In particular, whereas using a barotropic equation of state produces more brown dwarfs than stars, the inclusion of radiative feedback results in a stellar mass function and a ratio of brown dwarfs to stars in good agreement with observations of Galactic star-forming regions. In addition, many of the other statistical properties of the stars and brown dwarfs are in reasonable agreement with observations, including multiplicity as a function of primary mass, the frequency of very-low-mass binaries, and general trends for the mass ratio and separation distributions of binaries. We also examine the velocity dispersion of the stars, the distributions of disc truncation radii due to dynamical interactions, and coplanarity of orbits and sink particle spins in multiple systems. Overall, the calculation produces a cluster of stars whose statistical properties are difficult to distinguish from observed systems, implying that gravity, hydrodynamics, and radiative feedback are the primary ingredients for determining the origin of the statistical properties of low-mass stars.

We report results from radiation hydrodynamical simulations of the collapse of molecular cloud cores to form protostars. The calculations follow the formation and evolution of the first hydrostatic core/disc, the collapse to form a stellar core, and effect of stellar core formation on the surrounding disc and envelope. Past barotropic calculations have shown that rapidly-rotating first cores evolve into `pre-stellar discs' with radii up to ~100 AU that may last thousands of years before a stellar core forms. We investigate how the inclusion of a realistic equation of state and radiative transfer alters this behaviour, finding that the qualitative behaviour is similar, but that the pre-stellar discs may last 1.5-3 times longer in the more realistic calculations. The masses, radii, and lifetimes of the discs increase for initial molecular cloud cores with faster rotation rates. In the most extreme case we model, a pre-stellar disc with a mass of 0.22 Msun and a radius of ~100 AU can form in a solar-mass cloud and last several thousand years before a stellar core is formed. Such large, massive objects may be imaged using ALMA. Fragmentation of these massive discs may also provide an effective route to binary and multiple star formation, before radiative feedback from accretion onto the stellar core can inhibit fragmentation. Once collapse to form a stellar core occurs within the pre-stellar disc, the radiation hydrodynamical simulations produce qualitatively different behaviour from the barotropic calculations due to the accretion energy released. This drives a shock wave through the circumstellar disc and launches a bipolar outflow even in the absence of magnetic fields.

Dynamical interactions between planets and debris disks may sculpt the disk structure and impact planetary orbits, but only a few systems with both imaged planets and spatially resolved debris disks are known. With the Caltech Submm Observatory (CSO), we have observed the HR 8799 debris disk at 350{\mu}m. The 350{\mu}m map is the first spatially resolved measurement of the debris disk encircling the HR 8799 planetary system at this wavelength. Both the flux and size of the emission are consistent with a Kuiper belt of dust extending from ~100-300 AU. Although the resolution of the current map is limited, the map shows an indication of offset asymmetric emission, and several scenarios for this possibility are explored with radiative transfer calculations of a star-disk system and N-body numerical simulations of planet-disk interactions with parameters representative of the HR 8799 system.

We perform three-dimensional self-gravitating radiative transfer simulations of protoplanet migration in circumstellar discs to explore the impact upon migration of the radial temperature profiles in these discs. We model protoplanets with masses ranging between 10-100 M\bigoplus, in discs with surface density profiles of {\Sigma} \varpropto r^-1/2, and temperature profiles of the form T \varpropto r^-{\beta}, where {\beta} ranges 0-2. We find that steep ({\beta} > 1) temperature profiles lead to outward migration of low mass protoplanets in interstellar grain opacity discs, but in more optically thin discs the migration is always inwards. The trend in migration rates with changing {\beta} obtained from our models shows good agreement with those obtained using recent analytic descriptions which include consideration of the co-orbital torques and their saturation. We find that switching between two models of the protoplanet, one in which accretion acts by evacuating gas and one in which gas piles up on a surface to form an atmosphere, leads to a small shift in the migration rates. If comparing these models in discs with conditions which lead to a marginally inward migration, the small shift can lead to outward migration. However, the direction and speed of migration is dominated by disc conditions rather than by the specific prescription used to model the flow near the protoplanet.

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics (SPH) simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results certainly suggest that beta_{crit} is larger than previously thought. However, the absence of convergence also questions whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

We carry out a resolution study on the fragmentation boundary of self-gravitating discs. We perform three-dimensional Smoothed Particle Hydrodynamics simulations of discs to determine whether the critical value of the cooling timescale in units of the orbital timescale, beta_{crit}, converges with increasing resolution. Using particle numbers ranging from 31,250 to 16 million (the highest resolution simulations to date) we do not find convergence. Instead, fragmentation occurs for longer cooling timescales as the resolution is increased. These results suggest that at the very least, the critical value of the cooling timescale is longer than previously thought. However, the absence of convergence also raises the question of whether or not a critical value exists. In light of these results, we caution against using cooling timescale or gravitational stress arguments to deduce whether gravitational instability may or may not have been the formation mechanism for observed planetary systems.

We investigate the formation of both clustered and distributed populations of young stars in a single molecular cloud. We present a numerical simulation of a 10,000 solar mass elongated, turbulent, molecular cloud and the formation of over 2500 stars. The stars form both in stellar clusters and in a distributed mode which is determined by the local gravitational binding of the cloud. A density gradient along the major axis of the cloud produces bound regions that form stellar clusters and unbound regions that form a more distributed population. The initial mass function also depends on the local gravitational binding of the cloud with bound regions forming full IMFs whereas in the unbound, distributed regions the stellar masses cluster around the local Jeans mass and lack both the high-mass and the low-mass stars. The overall efficiency of star formation is ~ 15 % in the cloud when the calculation is terminated, but varies from less than 1 % in the the regions of distributed star formation to ~ 40 % in regions containing large stellar clusters. Considering that large scale surveys are likely to catch clouds at all evolutionary stages, estimates of the (time-averaged) star formation efficiency for the giant molecular cloud reported here is only ~ 4 %. This would lead to the erroneous conclusion of 'slow' star formation when in fact it is occurring on a dynamical timescale.

We investigate the fragmentation criterion in massive self-gravitating discs. We present new analysis of the fragmentation conditions which we test by carrying out global three-dimensional numerical simulations. Whilst previous work has placed emphasis on the cooling timescale in units of the orbital timescale, \beta , we find that at a given radius the surface mass density (i.e. disc mass and profile) and star mass also play a crucial role in determining whether a disc fragments or not as well as where in the disc fragments form. We find that for shallow surface mass density profiles (p<2, where \Sigma \propto R^{-p}), fragments form in the outer regions of the disc. However for steep surface mass density profiles (p \gtrsim 2), fragments form in the inner regions of a disc. In addition, we also find that the critical value of the cooling timescale in units of the orbital timescale found in previous simulations is only applicable to certain disc surface mass density profiles and for particular disc radii and is not a general rule for all discs. We find an empirical fragmentation criteria between the cooling timescale in units of the orbital timescale, \beta , the surface mass density, the star mass and the radius.

We calculate radial migration rates of protoplanets in laminar minimum mass solar nebula discs using three-dimensional self-gravitating radiation hydrodynamical (RHD) models. The protoplanets are free to migrate, whereupon their migration rates are measured. For low mass protoplanets (10-50 M_\oplus) we find increases in the migration timescales of up to an order of magnitude between locally-isothermal and RHD models. In the high-mass regime the migration rates are changed very little. These results are arrived at by calculating migration rates in locally-isothermal models, before sequentially introducing self-gravity, and radiative transfer, allowing us to isolate the effects of the additional physics. We find that using a locally-isothermal equation of state, without self-gravity, we reproduce the migration rates obtained by previous analytic and numerical models. We explore the impact of different protoplanet models, and changes to their assumed radii, upon migration. The introduction of self-gravity gives a slight reduction of the migration rates, whilst the inertial mass problem, which has been proposed for high mass protoplanets with circumplanetary discs, is reproduced. Upon introducing radiative transfer to models of low mass protoplanets (\approx 10 M_\oplus), modelled as small radius accreting point masses, we find outward migration with a rate of approximately twice the analytic inward rate. However, when modelling such a protoplanet in a more realistic manner, with a surface which enables the formation of a deep envelope, this outward migration is not seen.

We present the first non-LTE, co-moving frame molecular line calculations of a star-forming cluster simulated using smoothed particle hydrodynamics (SPH), from which we derive high-resolution synthetic observations. We have resampled a particle representation onto an adaptive mesh and self-consistently solved the equations of statistical equilibrium in the co-moving frame, using TORUS, a three-dimensional adaptive mesh refined (AMR) radiative transfer (RT) code. We verified the applicability of the code to the conditions of the SPH simulation by testing its output against other codes. We find that the level populations obtained for optically thick and thin scenarios closely match the ensemble average of the other codes. We have used the code to obtain non-LTE level populations of multiple molecular species throughout the cluster and have created three-dimensional velocity-resolved spatial maps of the emergent intensity. Line profiles of cores traced by N2H+ (1-0) are compared to probes of low density gas, 13CO (1-0) and C18O (1-0), surrounding the cores along the line of sight. The relative differences of the line-centre velocities are shown to be small compared to the velocity dispersion, matching recent observations. We conclude that one cannot reject competitive accretion as a viable theory of star formation based on observed velocity profiles.

We review the statistical properties of stars and brown dwarfs obtained from the first hydrodynamical simulation of star cluster formation to produce more than a thousand stars and brown dwarfs while simultaneously resolving the lowest mass brown dwarfs (those with masses set by the opacity limit for fragmentation), binaries with separations down to 1 AU, and discs with radii greater than 10 AU. In particular, we present the eccentricity distribution of the calculation's very-low-mass and brown dwarf binaries which has not been previously published.

We carry out global three-dimensional radiation hydrodynamical simulations of self-gravitating accretion discs to determine if, and under what conditions, a disc may fragment to form giant planets. We explore the parameter space (in terms of the disc opacity, temperature and size) and include the effect of stellar irradiation. We find that the disc opacity plays a vital role in determining whether a disc fragments. Specifically, opacities that are smaller than interstellar Rosseland mean values promote fragmentation (even at small radii, R < 25AU) since low opacities allow a disc to cool quickly. This may occur if a disc has a low metallicity or if grain growth has occurred. With specific reference to the HR 8799 planetary system, given its star is metal-poor, our results suggest that the formation of its imaged planetary system could potentially have occurred by gravitational instability. We also find that the presence of stellar irradiation generally acts to inhibit fragmentation (since the discs can only cool to the temperature defined by stellar irradiation). However, fragmentation may occur if the irradiation is sufficiently weak that it allows the disc to attain a low Toomre stability parameter.

We present results from the first three-dimensional radiation hydrodynamical calculations to follow the collapse of a molecular cloud core beyond the formation of the stellar core. We find the energy released by the formation of the stellar core, within the optically-thick first hydrostatic core, is comparable to the binding energy of the disc-like first core. This heats the inner regions of the disc, drives a shock wave through the disc, dramatically decreases the accretion rate on to the stellar core, and launches a temporary bipolar outflow perpendicular to the rotation axis that travels in excess of 50 AU into the infalling envelope. This outburst may assist the young protostar in launching a conventional magnetic jet. Furthermore, if these events are cyclic, they may provide a mechanism for intense bursts of accretion separated by long periods of relatively quiescent accretion which can potentially solve both the protostellar luminosity problem and the apparent age spread of stars in young clusters. Such outbursts may also provide a formation mechanism for the chondrules found in meteorites, with the outflow transporting them to large distances in the circumstellar disc.

Star formation is a complex process involving the interplay of many physical effects, including gravity, turbulent gas dynamics, magnetic fields and radiation. Our understanding of the process has improved substantially in recent years, primarily as a result of our increased ability to incorporate the relevant physics in numerical calculations of the star formation process. In this contribution we present an overview of our recent studies of star cluster formation in turbulent, magnetised clouds using self-gravitating radiation-magnetohydrodynamics calculations (Price and Bate 2008, 2009). Our incorporation of magnetic fields and radiative transfer into the Smoothed Particle Hydrodynamics method are discussed. We highlight how magnetic fields and radiative heating of the gas around newborn stars can solve several of the key puzzles in star formation, including an explanation for why star formation is such a slow and inefficient process. However, the presence of magnetic fields at observed strengths in collapsing protostellar cores also leads to problems on smaller scales, including a difficulty in forming protostellar discs and binary stars (Price and Bate 2007, Hennebelle and Teyssier 2008), which suggests that our understanding of the role of magnetic fields in star formation is not yet complete.

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.

We investigate the evolution, following gas dispersal, of a star cluster produced from a hydrodynamical calculation. We find that when the gas, initially comprising 60% of the mass, is removed, the system settles into a bound cluster containing 30-40% of the stellar mass surrounding by an expanding halo of ejected stars. The bound cluster expands from an initial radius of <0.05 pc to 1-2 pc over 4-10 Myr, depending on how quickly the gas is removed, implying that stellar clusters may begin with far higher stellar densities than usually assumed. With rapid gas dispersal the most massive stars are found to be mass segregated for the first ~1 Myr of evolution, but classical mass segregation only develops for cases with long gas removal timescales. Eventually, many of the most massive stars are expelled from the bound cluster. Despite the high initial stellar density and the extensive dynamical evolution of the system, we find that the stellar multiplicity is almost constant during the 10 Myr of evolution. This is because the primordial multiple systems are formed in a clustered environment and, thus, by their nature are already resistant to further evolution. The majority of multiple system evolution is confined to the decay of high-order systems and the formation of a significant population of very wide (10^4-10^5 AU) multiple systems in the expanding halo. This formation mechanism for wide binaries potentially solves the problem of how most stars apparently form in clusters and yet a substantial population of wide binaries exist in the field. Many of these wide binaries and the binaries produced by the decay of high-order multiple systems have unequal mass components, potentially solving the problem that hydrodynamical simulations of star formation are found to under-produce unequal-mass solar-type binaries.

We present hydrodynamical models of the grand design spiral M51 (NGC 5194), and its interaction with its companion NGC 5195. Despite the simplicity of our models, our simulations capture the present day spiral structure of M51 remarkably well, and even reproduce details such as a kink along one spiral arm, and spiral arm bifurcations. We investigate the offset between the stellar and gaseous spiral arms, and find at most times (including the present day) there is no offset between the stars and gas to within our error bars. We also compare our simulations with recent observational analysis of M51. We compute the pattern speed versus radius, and like the observations, find no single global pattern speed. We also show that the spiral arms cannot be fitted well by logarithmic spirals. We interpret these findings as evidence that M51 does not exhibit a quasi-steady density wave, as would be predicted by density wave theory. The internal structure of M51 derives from the complicated and dynamical interaction with its companion, resulting in spiral arms showing considerable structure in the form of short-lived kinks and bifurcations. Rather than trying to model such galaxies in terms of global spiral modes with fixed pattern speeds, it is more realistic to start from a picture in which the spiral arms, while not being simple material arms, are the result of tidally induced kinematic density `waves' or density patterns, which wind up slowly over time.

We investigate the evolution of the relative angle between the stellar rotation axis and the circumstellar disc axis of a star that forms in a stellar cluster from the collapse of a turbulent molecular cloud. This is an inherently chaotic environment with variable accretion, both in terms of rate and the angular momentum of the material, and dynamical interactions between stars. We find that the final stellar rotation axis and disc spin axis can be strongly misaligned, but this occurs primarily when the disc is truncated by a dynamical encounter so that the final disc rotation axis depends simply on what fell in last. This may lead to planetary systems with orbits that are misaligned with the stellar rotation axis, but only if the final disc contains enough mass to form planets. We also investigate the time variability of the inner disc spin axis, which is likely to determine the direction of a protostellar jet. We find that the jet direction varies more strongly for lighter discs, such as those that have been truncated by dynamical interactions or have suffered a period of rapid accretion. Finally, we note that variability of the angular momentum of the material accreting by a star implies that the internal velocity field of such stars may be more complicated than that of aligned differential rotation.

The M51 galaxy is one of best examples of grand design structure in galaxies. The most likely origin of the spiral arms in M51 is the ongoing interaction with NGC 5195. Here we report recent calculations which model the orbit of M51 and NGC 5195. These calculations, which for the first time focus on the gas dynamics, are able to reproduce even detailed features in M51, as well as demonstrate the nature of the spiral structure of M51.

We investigate the dependence of stellar properties on the initial kinematic structure of the gas in star-forming molecular clouds. We compare the results from two large-scale hydrodynamical simulations of star cluster formation that resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell and Bromm. The initial conditions of the two calculations are identical, but in the new simulation the power spectrum of the velocity field imposed on the cloud initially and allowed to decay is biased in favour of large-scale motions. Whereas the calculation of Bate et al. began with a power spectrum P(k) ~ k^{-4} to match the Larson scaling relations for the turbulent motions observed in molecular clouds, the new calculation begins with a power spectrum P(k) ~ k^{-6}. Despite this change to the initial motions in the cloud and the resulting density structure of the molecular cloud, the resulting stellar properties resulting from the two calculations are indistinguishable. This demonstrates that the results of such hydrodynamical calculations of star cluster formation are relatively insensitive to the initial conditions. It is also consistent with the fact that the statistical properties of stars and brown dwarfs (e.g. the stellar initial mass function) are observed to be relatively invariant within our Galaxy and do not appear to depend on environment.

We investigate the effects of magnetic fields and radiative protostellar feedback on the star formation process using self-gravitating radiation magnetohydrodynamical calculations. We present results from a series of calculations of the collapse of 50 solar mass molecular clouds with various magnetic field strengths and with and without radiative transfer. We find that both magnetic fields and radiation have a dramatic impact on star formation, though the two effects are in many ways complementary. Magnetic fields primarily provide support on large scales to low density gas, whereas radiation is found to strongly suppress small-scale fragmentation by increasing the temperature in the high-density material near the protostars. With strong magnetic fields and radiative feedback the net result is an inefficient star formation process with a star formation rate of ~< 10% per free-fall time that approaches the observed rate, although we have only been able to follow the calculations for ~1/3 of a free-fall time beyond the onset of star formation.

We present accretion rates obtained from three-dimensional self-gravitating radiation hydrodynamical models of giant planet growth. We investigate the dependence of accretion rates upon grain opacity and core/protoplanet mass. The accretion rates found for low mass cores are inline with the results of previous one-dimensional models that include radiative transfer.

We present preliminary results on the ability of self-gravitating discs to cool in response to their internal heating. These discs are modelled using a Smoothed Particle Hydrodynamics (SPH) code with radiative transfer (Whitehouse, Bate & Monaghan 2005) and we investigate the ability of these discs to maintain a state of thermal equilibrium with their boundaries as an indication of their likelihood to fragment.

We investigate the effect of radiative feedback on the star formation process using radiation hydrodynamical simulations. We repeat the previous hydrodynamical star cluster formation simulations of Bate et al., and Bate & Bonnell, but we use a realistic gas equation of state and radiative transfer in the flux-limited diffusion approximation rather than the original barotropic equation of state. Whereas star formation in the barotropic simulations continued unabated until the simulation was stopped, we find that radiative feedback essentially terminates the production of new objects after roughly one local dynamical time. Radiative feedback also dramatically decreases the propensity of massive circumstellar discs to fragment and inhibits fragmentation of other dense gas (e.g. filaments) close to existing protostellar objects. These two effects decrease the numbers of protostars formed by a factor of approximately 4. Whereas the original simulations produced more brown dwarfs than stars, the radiative feedback results in a ratio of stars to brown dwarfs of approximately 5:1, in much better agreement with observations. Most importantly, we find that although the characteristic stellar mass in the original calculations scaled linearly with the initial mean Jeans mass, when radiative feedback is included the characteristic stellar mass is indistinguishable for the two calculations. We thus propose that the reason the observed initial mass function appears to be universal is due to self-regulation of the star formation process by radiative feedback. We present an analytic argument showing how a characteristic mass may be derived that is relatively independent of initial conditions such as the cloud's density. (Abridged)

We investigate the properties of circumplanetary discs formed in three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by protoplanets. We determine disc sizes, scaleheights, and density and temperature profiles for different protoplanet masses, in solar nebulae of differing grain opacities. We find that the analytical prediction of circumplanetary disc radii in an evacuated gap (R_Hill/3) from Quillen & Trilling (1998) yields a good estimate for discs formed by high mass protoplanets. The radial density profiles of the circumplanetary discs may be described by power-laws between r^-2 and r^-3/2. We find no evidence for the ring-like density enhancements that have been found in some previous models of circumplanetary discs. Temperature profiles follow a ~r^-7/10 power-law regardless of protoplanet mass or nebula grain opacity. The discs invariably have large scaleheights (H/r > 0.2), making them thick in comparison with their encompassing circumstellar discs, and they show no flaring.

We present results from three-dimensional, self-gravitating radiation hydrodynamical models of gas accretion by planetary cores. In some cases, the accretion flow is resolved down to the surface of the solid core -- the first time such simulations have been performed. We investigate the dependence of the gas accretion rate upon the planetary core mass, and the surface density and opacity of the encompassing protoplanetary disc. Accretion of planetesimals is neglected. We find that high-mass protoplanets are surrounded by thick circumplanetary discs during their gas accretion phase but, contrary to locally-isothermal calculations, discs do not form around accreting protoplanets with masses ~< 50M_Earth when radiation hydrodynamical simulations are performed, even if the grain opacity is reduced from interstellar values by a factor of 100. We find that the opacity of the gas plays a large role in determining the accretion rates for low-mass planetary cores. For example, reducing the opacities from interstellar values by a factor of 100 leads to roughly an order of magnitude increase in the accretion rates for 10-20M_Earth protoplanets. The dependence on opacity becomes less important in determining the accretion rate for more massive cores where gravity dominates the effects of thermal support and the protoplanet is essentially accreting at the runaway rate. Finally, for low-mass planetary cores (~< 20M_Earth), we obtain accretion rates that are in agreement with previous one-dimensional quasi-static models. This indicates that three-dimensional hydrodynamical effects may not significantly alter the gas accretion timescales that have been obtained from quasi-static models.

We report the statistical properties of stars, brown dwarfs and multiple systems obtained from the largest hydrodynamical simulation of star cluster formation to date that resolves masses down to the opacity limit for fragmentation (a few Jupiter masses). The simulation is essentially identical to that of Bate, Bonnell & Bromm except that the initial molecular cloud is larger and more massive. It produces more than 1250 stars and brown dwarfs, providing unprecedented statistical information that can be compared with observational surveys. We find that hydrodynamical/sink particle simulations can reproduce many of the observed stellar properties very well. Binarity as a function of primary mass, the frequency of very-low-mass (VLM) binaries, general trends for the separation and mass ratio distributions of binaries, and the relative orbital orientations of triples systems are all in reasonable agreement with observations. We also examine the radial variations of binarity, velocity dispersion, and mass function in the resulting stellar cluster and the distributions of disc truncation radii due to dynamical interactions. For VLM binaries, we find that their frequency when using small accretion radii and gravitational softening is similar to that expected from observational surveys (approximately 20 percent). We also find that VLM binaries evolve from wide, unequal-mass systems towards close equal-mass systems as they form. The two main deficiencies of the calculations are that they over produce brown dwarfs relative to stars and that there are too few unequal mass binaries with K and G-dwarf primaries. [Abridged]

We investigate the formation of brown dwarfs and very low-mass stars through the gravitational fragmentation of infalling gas into stellar clusters. The gravitational potential of a forming stellar cluster provides the focus that attracts gas from the surrounding molecular cloud. Structures present in the gas grow, forming filaments flowing into the cluster centre. These filaments attain high gas densities due to the combination of the cluster potential and local self-gravity. The resultant Jeans masses are low, allowing the formation of very low-mass fragments. The tidal shear and high velocity dispersion present in the cluster preclude any subsequent accretion thus resulting in the formation of brown dwarfs or very low-mass stars. Ejections are not required as the brown dwarfs enter the cluster with high relative velocities, suggesting that their disc and binary properties should be similar to that of low-mass stars. This mechanism requires the presence of a strong gravitational potential due to the stellar cluster implying that brown dwarf formation should be more frequent in stellar clusters than in distributed populations of young stars. Brown dwarfs formed in isolation would require another formation mechanism such as due to turbulent fragmentation.

Magnetic fields are important at every scale in the star formation process: from the dynamics of the ISM in galaxies, to the collapse of turbulent molecular clouds to form stars and in the fragmentation of individual star forming cores. The recent development of a robust algorithm for MHD in the Smoothed Particle Hydrodynamics method has enabled us to perform simulations of star formation including magnetic fields at each of these scales. This paper focusses on three questions in particular: What is the effect of magnetic fields on fragmentation in star forming cores? How do magnetic fields affect the collapse of turbulent molecular clouds to form stars? and: What effect do magnetic fields have on the dynamics of the interstellar medium?

We examine the effect of magnetic fields on star cluster formation by performing simulations following the self-gravitating collapse of a turbulent molecular cloud to form stars in ideal MHD. The collapse of the cloud is computed for global mass-to-flux ratios of infinity, 20, 10, 5 and 3, that is using both weak and strong magnetic fields. Whilst even at very low strengths the magnetic field is able to significantly influence the star formation process, for magnetic fields with plasma beta < 1 the results are substantially different to the hydrodynamic case. In these cases we find large-scale magnetically-supported voids imprinted in the cloud structure; anisotropic turbulent motions and column density structure aligned with the magnetic field lines, both of which have recently been observed in the Taurus molecular cloud. We also find strongly suppressed accretion in the magnetised runs, leading to up to a 75% reduction in the amount of mass converted into stars over the course of the calculations and a more quiescent mode of star formation. There is also some indication that the relative formation efficiency of brown dwarfs is lower in the strongly magnetised runs due to the reduction in the importance of protostellar ejections.

We present first results of our simulations of magnetic fields in the formation of single and binary stars using a recently developed method for incorporating Magnetohydrodynamics (MHD) into the Smoothed Particle Hydrodynamics (SPH) method. An overview of the method is presented before discussing the effect of magnetic fields on the formation of circumstellar discs around young stars. We find that the presence of magnetic fields during the disc formation process can lead to significantly smaller and less massive discs which are much less prone to gravitational instability. Similarly in the case of binary star formation we find that magnetic fields, overall, suppress fragmentation. However these effects are found to be largely driven by magnetic pressure. The relative importance of magnetic tension is dependent on the orientation of the field with respect to the rotation axis, but can, with the right orientation, lead to a dilution of the magnetic pressure-driven suppression of fragmentation.

We have performed magnetohydrodynamic (MHD) simulations of the collapse and fragmentation of molecular cloud cores using a new algorithm for MHD within the smoothed particle hydrodynamics (SPH) method, that enforces the zero magnetic divergence constraint. We find that the support provided by magnetic fields over thermal pressure alone has several important effects on fragmentation and the formation of binary and multiple systems, and on the properties of massive circumstellar discs. The extra support suppresses the tendency of molecular cloud cores to fragment due to either initial density perturbations or disc fragmentation. Furthermore, unlike most previous studies, we find that magnetic pressure plays the dominant role in inhibiting fragmentation rather than magnetic tension or magnetic braking. In particular, we find that if the magnetic field is aligned with the rotation axis of the molecular cloud core, the effects of the magnetic field on fragmentation and disc structure are almost entirely due to magnetic pressure, while if the rotation axis is initially perpendicular to the magnetic field, magnetic tension plays a greater role and can actually aid fragmentation. Despite these effects, and contrary to several past studies, we find that strongly-perturbed molecular cloud cores are able to fragment to form wide binary systems even in the presence of quite strong magnetic fields. For massive circumstellar discs, we find that slowing of the collapse caused by the magnetic support decreases the mass infall rate on to the disc and, thus, weakens gravitational instabilities in young massive circumstellar discs. This not only reduces the likelihood that they will fragment, but also decreases the importance of spiral density waves in providing angular momentum transport and in promoting planet formation.

Hydrodynamical simulations of star formation indicate that the motions of protostars through their natal molecular clouds may be crucial in determining the properties of stars through competitive accretion and dynamical interactions. Walsh, Myers & Burton recently investigated whether such motions might be observable in the earliest stages of star formation by measuring the relative shifts of line-centre velocities of low- and high-density tracers of low-mass star-forming cores. They found very small (~0.1 km/s) relative motions. In this paper, we analyse the hydrodynamical simulation of Bate, Bonnell & Bromm and find that it also gives small relative velocities between high-density cores and low-density envelopes, despite the fact that competitive accretion and dynamical interactions occur between protostars in the simulation. Thus, the simulation is consistent with the observations in this respect. However, we also find some differences between the simulation and the observations. Overall, we find that the high-density gas has a higher velocity dispersion than that observed by Walsh et al. We explore this by examining the dependence of the gas velocity dispersion on density and its evolution with time during the simulation. We find that early in the simulation the gas velocity dispersion decreases monotonically with increasing density, while later in the simulation, when the dense cores have formed multiple objects, the velocity dispersion of the high-density gas increases. Thus, the simulation is in best agreement with the observations early on, before many objects have formed in each dense core.

We investigate the interaction between a giant planet and a viscous circumstellar disk by means of high-resolution, two-dimensional hydrodynamical simulations. We consider planet masses that range from 1 to 3 Jupiter masses (Mjup) and initial orbital eccentricities that range from 0 to 0.4. We find that a planet can cause eccentricity growth in a disk region adjacent to the planet's orbit, even if the planet's orbit is circular. Disk-planet interactions lead to growth in a planet's orbital eccentricity. The orbital eccentricities of a 2 Mjup and a 3 Mjup planet increase from 0 to 0.11 within about 3000 orbits. Over a similar time period, the orbital eccentricity of a 1 Mjup planet grows from 0 to 0.02. For a case of a 1 Mjup planet with an initial eccentricity of 0.01, the orbital eccentricity grows to 0.09 over 4000 orbits. Radial migration is directed inwards, but slows considerably as a planet's orbit becomes eccentric. If a planet's orbital eccentricity becomes sufficiently large, e > ~0.2, migration can reverse and so be directed outwards. The accretion rate towards a planet depends on both the disk and the planet orbital eccentricity and is pulsed over the orbital period. Planet mass growth rates increase with planet orbital eccentricity. For e~0.2 the mass growth rate of a planet increases by approximately 30% above the value for e=0. For e > ~0.1, most of the accretion within the planet's Roche lobe occurs when the planet is near the apocenter. Similar accretion modulation occurs for flow at the inner disk boundary which represents accretion toward the star.

Competitive accretion, a process to explain the origin of the IMF, occurs when stars in a common gravitational potential accrete from a distributed gaseous component. We show that concerns recently raised on the efficiency of competitive accretion are incorrect as they use globally averaged properties which are inappropriate for the detailed physics of a forming stellar cluster. A full treatment requires a realistic treatment of the cluster potential, the distribution of turbulent velocities and gas densities. Accreting gas does not travel at the global virial velocity of the system due to the velocity-sizescale relation inherent in turbulent gas and due to the lower velocity dispersion of small-N clusters in which much of the accretion occurs. Stars located in the gas-rich centres of such systems initially accrete from low relative velocity gas attaining larger masses before needing to accrete the higher velocity gas. Stars not in the centres of such potentials, or that enter the cluster later when the velocity dispersion is higher, do not accrete significantly and thus retain their low-masses. In competitive accretion, most stars do not continue to accrete significantly such that their masses are set from the fragmentation process. It is the few stars which continue to accrete that become higher-mass stars. Competitive accretion is therefore likely to be responsible for the formation of higher-mass stars and can explain the mass distribution, mass segregation and binary frequency of these stars. Global kinematics of competitive accretion models include large-scale mass infall, with mean inflow velocities of order 0.5 km/s at scales of 0.5 pc, but infall signatures are likely to be confused by the large tangential velocities and the velocity dispersion present.

We use numerical simulations of the fragmentation of a 1000 solar mass molecular cloud and the formation of a stellar cluster to study how the initial conditions for star formation affect the resulting initial mass function (IMF). In particular, we are interested in the relation between the thermal Jeans mass in a cloud and the knee of the initial mass function, i.e. the mass separating the region with a flat IMF slope from that typified by a steeper, Salpeter-like, slope. In three isothermal simulations with M_J=1 solar mass, M_J=2 solar masses and M_J=5 solar masses, the number of stars formed, at comparable dynamical times, scales roughly with the number of initial Jeans masses in the cloud. The mean stellar mass also increases (though less than linearly) with the initial Jeans mass in the cloud. It is found that the IMF in each case displays a prominent knee, located roughly at the mass scale of the initial Jeans mass. Thus clouds with higher initial Jeans masses produce IMFs which are shallow to higher masses. This implies that a universal IMF requires a physical mechanism that sets the Jeans mass to be near 1 solar mass. Simulations including a barotropic equation of state as suggested by Larson, with cooling at low densities followed by gentle heating at higher densities, are able to produce realistic IMFs with the knee located at approximately 1 solar mass, even with an initial M_J=5 solar masses. We therefore suggest that the observed universality of the IMF in the local Universe does not require any fine tuning of the initial conditions in star forming clouds but is instead imprinted by details of the cooling physics of the collapsing gas.

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.

We present the results of a series of calculations studying the collapse of molecular cloud cores performed using a three-dimensional smoothed particle hydr odynamics code with radiative transfer in the flux-limited diffusion approximation. The opacities and specific heat capacities are identical for each calculation. However, we find that the temperature evolution during the simulations varies significantly when starting from different initial conditions. Even spherically-symmetric clouds with different initial densities show markedly different development. We conclude that simple barotropic equations of state like those used in some previous calculations provide at best a crude approximation to the thermal behaviour of the gas. Radiative transfer is necessary to obtain accurate temperatures.

We describe a new, faster implicit algorithm for solving the radiation hydrodynamics equations in the flux-limited diffusion approximation for smoothed particle hydrodynamics. This improves on the method elucidated in Whitehouse & Bate by using a Gauss-Seidel iterative method rather than iterating over the exchange of energy between pairs of particles. The new algorithm is typically many thousands of times faster than the old one, which will enable more complex problems to be solved. The new algorithm is tested using the same tests performed by Turner & Stone for ZEUS-2D, and repeated by Whitehouse & Bate.

We present a model for the formation of high-mass close binary systems in the context of forming massive stars through gas accretion in the centres of stellar clusters. A low-mass wide binary evolves under mass accretion towards a high-mass close binary, attaining system masses of order 30-50 solar masses at separations of order 1 AU. The resulting high frequency of binary systems with two massive components is in agreement with observations. These systems are typically highly eccentric and may evolve to have periastron separations less than their stellar radii. Mergers of these binary systems are therefore likely and can lead to the formation of the most massive stars, circumventing the problem of radiation pressure stopping the accretion. The stellar density required to induce binary mergers is one million stars per cubic parsec, or approximately one percent that required for direct stellar collisions.

We investigate the formation of star clusters in an unbound GMC, where the supporting kinetic energy is twice as large as the cloud's self-gravity. This cloud manages to form a series of star clusters and disperse, all within roughly 2 crossing times (10 Myr), supporting recent claims that star formation is a rapid process. Simple assumptions about the nature of the star formation occurring in the clusters allows us to place an estimate for the star formation efficiency at about 5 to 10 %, consistent with observations. We also propose that unbound clouds can act as a mechanism for forming OB associations. The clusters that form in the cloud behave as OB subgroups. These clusters are naturally expanding from one another due to unbound nature of the flows that create them. The properties of the cloud we present here are are consistent with those of classic OB associations.

We present the first ever hydrodynamic calculations of star cluster formation that incorporate the effect of feedback from ionising radiation. In our simulations, the ionising source forms in the cluster core at the intersection of several dense filaments of inflowing gas. We show that these filaments collimate ionised outflows and suggest such an environmental origin for at least some observed outflows in regions of massive star formation. Our simulations show both positive feedback (i.e. promotion of star formation in neutral gas compressed by expanding HII regions) and negative feedback (i.e. suppression of the accretion flow in to the central regions). We show that the volume filling factor of ionised gas is very different in our simulations than would result from the case where the central source interacted with an azimuthally smoothed gas density distribution. As expected, gas density is the key parameter in determining whether clusters are unbound by photoionising radiation. Nevertheless, we find - on account of the acceleration of a small fraction of the gas to high velocities in the outflows - that the deposition in the gas of an energy that exceeds the binding energy of the cluster is not a sufficient criterion for unbinding the bulk of the cluster mass.

We investigate the migration rates of high-mass protoplanets embedded in accretion discs via two and three-dimensional hydrodynamical simulations. The simulations follow the planet's radial motion and employ a nested-grid code that allows for high resolution close to the planet. We concentrate on the possible role of the coorbital torques in affecting migration rates. We analyse two cases: (a) a Jupiter-mass planet in a low-mass disc and (b) a Saturn-mass planet in a high-mass disc. Planet migration in case (b) is much more susceptible to coorbital torques than in case (a). We find that the coorbital torques in both cases do not depend sensitively on whether the planet is allowed to migrate through the disc or is held on a fixed orbit. We also examine the dependence of the planet's migration rate on the numerical resolution near the planet. For case (a), numerical convergence is relatively easy to obtain, even when including torques arising from deep within the planet's Hill sphere. The migration rate in this case is numerically on order of the Type II migration rate and much smaller than the Type I rate, if the disc has 0.01 solar-masses inside 26 AU. Torques from within the Hill sphere provide a substantial opposing contribution to the migration rate. In case (b), the gas mass within the Hill sphere is larger than the planet's mass and convergence is more difficult to obtain. If the torques within the Hill sphere are ignored, convergence is more easily achieved but the migration rate is artificially large. At our highest resolution, the migration rate for case (b) is much less than the Type I rate, but somewhat larger than the Type II rate.

We investigate the dependence of stellar properties on the mean thermal Jeans mass in molecular clouds. We compare the results from the two largest hydrodynamical simulations of star formation to resolve the fragmentation process down to the opacity limit, the first of which was reported by Bate, Bonnell & Bromm. The initial conditions of the two calculations are identical except for the radii of the clouds, which are chosen so that the mean densities and mean thermal Jeans masses of the clouds differ by factors of nine and three, respectively. We find that the denser cloud, with the lower mean thermal Jeans mass, produces a higher proportion of brown dwarfs and has a lower characteristic (median) mass of the stars and brown dwarfs. This dependence of the initial mass function (IMF) on the density of the cloud may explain the observation that the Taurus star-forming region appears to be deficient in brown dwarfs when compared with the Orion Trapezium cluster. The new calculation also produces wide binaries (separations >20 AU), one of which is a wide binary brown dwarf system. Based on the hydrodynamical calculations, we develop a simple accretion/ejection model for the origin of the IMF. In the model, all stars and brown dwarfs begin with the same mass (set by the opacity limit for fragmentation) and grow in mass until their accretion is terminated stochastically by their ejection from the cloud through dynamically interactions. The model predicts that the main variation of the IMF in different star-forming environments should be in the location of the peak (due to variations in the mean thermal Jeans mass of the cloud) and in the substellar regime. However, the slope of the IMF at high-masses may depend on the dispersion in the accretion rates of protostars.

We describe the implementation and testing of a smoothed particle hydrodynamics (SPH) code that solves the equations of radiation hydrodynamics in the flux-limited diffusion (FLD) approximation. The SPH equations of radiation hydrodynamics for an explicit integration scheme are derived and tested. We also discuss the implementation of an implicit numerical scheme for solving the radiation equations that allows the system to be evolved using timesteps much longer than the explicit radiation timestep. The code is tested on a variety of one-dimensional radiation hydrodynamics problems including radiation propagating in an optically thin medium, optically thick and thin shocks, subcritical and supercritical radiating shocks, and a radiation dominated shock. Many of the tests were also performed by Turner and Stone (2001) to test their implementation of a FLD module for the ZEUS-2D code. The SPH code performs at least as well as the ZEUS-2D code in these tests.

We present three-dimensional Monte Carlo radiative transfer models of a very young (<10^5 years old) low mass (50 Msun) stellar cluster containing 23 stars and 27 brown dwarfs. The models use the density and the stellar mass distributions from the large-scale smoothed particle hydrodynamics (SPH) simulation of the formation of a low-mass stellar cluster by Bate, Bonnell and Bromm. We compute the spectral energy distributions of the whole cluster and the individual objects, simulated far-infrared Spitzer Space Telescope (SST) images and construct colour-colour diagrams (near-infrared HKL and SST mid-infrared bands). The presence of accretion discs around the light sources influences the morphology of the dust temperature structure on a large scale (up to a several 10^4 au). A considerable fraction of the interstellar dust is underheated compared to a model without the accretion discs because the radiation from the light sources is blocked/shadowed by the discs. The spectral energy distribution (SED) of the model cluster with accretion discs shows excess emission at lambda = 3-30 micron and lambda > 500 micron compared to that without accretion discs. Our model with accretion discs around each object shows a similar distribution of spectral index (2.2-20 micron) values (i.e. Class 0--III sources) as seen in the Rho Ophiuchus cloud. We confirm that the best diagnostics for identifying objects with accretion discs are mid-infrared (lambda = 3-10 micron) colours (e.g. SST IRAC bands) rather than HKL colours.

Abridged/ We present numerical results on the properties of young binary and multiple stellar systems. Our analysis is based on a series of SPH + Nbody simulations of the fragmentation of small molecular clouds, that fully resolve the opacity limit for fragmentation. We have produced a statistically significant number of stable multiple systems, with components separations in the range 1-1000 AU. At the end of the hydrodynamical evolution (0.5 Myr) we find that ~60% of stars and brown dwarfs are members of multiples systems, with about a third of these being low mass, weakly bound outliers in wide eccentric orbits. Our results imply that in the stellar regime most stars are in multiples (~80%) and that this fraction is an increasing function of primary mass. After Nbody integration to 10.5 Myr, the percentage of bound objects has dropped to ~40%, as most very low mass stars and brown dwarfs have been released to the field. Brown dwarfs are never found to be very close companions to stars (brown dwarf desert at very small separations), but one case exists of a brown dwarf companion at intermediate separations (10 AU). Our simulations can accommodate the existence of brown dwarf companions at large separations, but only if the primaries of these systems are themselves multiples. We have compared the outcome of our simulations with the properties of real stellar systems as deduced from the IR CM diagram of the Praesepe cluster and from spectroscopic and high-resolution imaging surveys of young clusters and the field.

We investigate the physical processes which lead to the formation of massive stars. Using a numerical simulation of the formation of a stellar cluster from a turbulent molecular cloud, we evaluate the relevant contributions of fragmentation and competitive accretion in determining the masses of the more massive stars. We find no correlation between the final mass of a massive star, and the mass of the clump from which it forms. Instead, we find that the bulk of the mass of massive stars comes from subsequent competitive accretion in a clustered environment. In fact, the majority of this mass infalls onto a pre-existing stellar cluster. Furthermore, the mass of the most massive star in a system increases as the system grows in numbers of stars and in total mass. This arises as the infalling gas is accompanied by newly formed stars, resulting in a larger cluster around a more massive star. High-mass stars gain mass as they gain companions, implying a direct causal relationship between the cluster formation process, and the formation of higher-mass stars therein.

Self-gravitating protostellar discs are unstable to fragmentation if the gas can cool on a time scale that is short compared to the orbital period. We use a combination of hydrodynamic simulations and N-body orbit integrations to study the long term evolution of a fragmenting disc with an initial mass ratio to the star of M_disc/M_star = 0.1. For a disc which is initially unstable across a range of radii, a combination of collapse and subsequent accretion yields substellar objects with a spectrum of masses extending (for a Solar mass star) up to ~0.01 M_sun. Subsequent gravitational evolution ejects most of the lower mass objects within a few million years, leaving a small number of very massive planets or brown dwarfs in eccentric orbits at moderately small radii. Based on these results, systems such as HD 168443 -- in which the companions are close to or beyond the deuterium burning limit -- appear to be the best candidates to have formed via gravitational instability. If massive substellar companions originate from disc fragmentation, while lower-mass planetary companions originate from core accretion, the metallicity distribution of stars which host massive substellar companions at radii of ~1 au should differ from that of stars with lower mass planetary companions.

Abridged: We have undertaken a series of hydrodynamical simulations of multiple star formation in small turbulent molecular clouds. Our goal is to determine the sensitivity of the properties of the resulting stars and brown dwarfs to variations in the initial conditions imposed. In this paper we report on the results obtained by applying two different initial turbulent velocity fields. The slope of the turbulent power-law spectrum alpha is set to -3 in half of the calculations and to -5 in the other half. We find that, whereas the stellar mass function seems to only be weakly dependent on the value of alpha, the sub-stellar mass function turns out to be more sensitive to the initial slope of the velocity field. We argue that, since the role of turbulence is to create substructure from which gravitational instabilities may grow, variations in other initial conditions that also determine the fragmentation process are likely to affect the shape of the sub-stellar mass function as well. The absence of many planetary mass `free-floaters' in our simulations, especially in the mass range 1-10 MJ, suggests that, if these objects are abundant, they are likely to form by similar mechanisms to those thought to operate in quiescent accretion discs, instead of via instabilities in gravitationally unstable discs. Finally, we find that multiple and single stars share comparable kinematical properties, both populations being able to attain velocities in the range 1-10 km/s. From these values we draw the conclusion that only low-mass star-forming regions such as Taurus-Auriga or Ophiuchus, where the escape speed is low, might have suffered some depletion of its single and binary stellar population.

It has already been shown, using a local model, that accretion discs with cooling times t_cool <= 3 Omega^-1 fragment into gravitationally bound objects, while those with cooling times t_cool > 3 Omega^-1 evolve into a quasi-steady state. We present results of three-dimensional simulations that test if the local result still holds globally. We find that the fragmentation boundary is close to that determined using the local model, but that fragmentation may occur for longer cooling times when the disc is more massive or when the mass is distributed in such a way as to make a particular region of the disc more susceptible to the growth of the gravitational instability. These results have significant implications for the formation of of gaseous planets in protoplanetary discs and also for the redistribution of angular momentum which could be driven by the presence of relatively massive, bound objects within the disc.

We present an analytical method for constraining the substellar initial mass function in globular clusters, based on the observed frequency of transit events. Globular clusters typically have very high stellar densities where close encounters are relatively common, and thus tidal capture can occur to form close binary systems. Encounters between main sequence stars and lower-mass objects can result in tidal capture if the mass ratio is > 0.01. If brown dwarfs exist in significant numbers, they too will be found in close binaries, and some fraction of their number should be revealed as they transit their stellar companions. We calculate the rate of tidal capture of brown dwarfs in both segregated and unsegregated clusters, and find that the tidal capture is more likely to occur over an initial relaxation time before equipartition occurs. The lack of any such transits in recent HST monitoring of 47 Tuc implies an upper limit on the frequency of brown dwarfs (< 15 % relative to stars) which is significantly below that measured in the galactic field and young clusters.

Recent surveys of star forming regions have shown that most stars, and probably all massive stars, are born in dense stellar clusters. The mechanism by which a molecular cloud fragments to form several hundred to thousands of individual stars has remained elusive. Here, we use a numerical simulation to follow the fragmentation of a turbulent molecular cloud and the subsequent formation and early evolution of a stellar cluster containing more than 400 stars. We show that the stellar cluster forms through the hierarchical fragmentation of a turbulent molecular cloud. This leads to the formation of many small subclusters which interact and merge to form the final stellar cluster. The hierarchical nature of the cluster formation has serious implications in terms of the properties of the new-born stars. The higher number-density of stars in subclusters, compared to a more uniform distribution arising from a monolithic formation, results in closer and more frequent dynamical interactions. Such close interactions can truncate circumstellar discs, harden existing binaries, and potentially liberate a population of planets. We estimate that at least one-third of all stars, and most massive stars, suffer such disruptive interactions.

We present results from high-resolution hydrodynamical simulations which explore the effects of small scale clustering in star-forming regions. A large ensemble of small-N clusters with 5 stellar seeds have been modelled and the resulting properties of stars and brown dwarfs statistically derived and compared with observational data. Close dynamical interactions between the protostars and competitive accretion driven by the cloud collapse are shown to produce a distribution of final masses which is bimodal, with most of the mass residing in the binary components. When convolved with a suitable core mass function, the final distribution of masses resembles the observed IMF, both in the stellar and sub-stellar regime. Binaries and single stars are found to constitute two kinematically distinct populations, with about half of the singles attaining velocities > 2 km/s, which might deprive low mass star-forming regions of their lightest members in a few crossing times. The eccentricity distribution of binaries and multiples is found to follow a distribution similar to that of observed long period (uncircularized) binaries. We predict that many close binary stars should have wide brown dwarfs companions. Brown dwarfs and, in general, very low mass stars, would be rare as pure binary companions. The binary fraction should be a decreasing function of primary mass, with low-mass or sub-stellar primaries being scarce. Where such binaries exist, they are either expected to be close enough (semi-major axis ~ 10 AU) to survive strong interactions with more massive binaries or else born in very small molecular cloud cores.

We analyse the non-linear, three-dimensional response of a gaseous, viscous protoplanetary disc to the presence of a planet of mass ranging from one Earth mass (1 M$_e$) to one Jupiter mass (1 M$_J$) by using the ZEUS hydrodynamics code. We determine the gas flow pattern, and the accretion and migration rates of the planet. The planet is assumed to be in a fixed circular orbit about the central star. It is also assumed to be able to accrete gas without expansion on the scale of its Roche radius. Only planets with masses $M \gsim 0.1$ M$_J$ produce significant perturbations in the disc's surface density. The flow within the Roche lobe of the planet is fully three-dimensional. Gas streams generally enter the Roche lobe close to the disc midplane, but produce much weaker shocks than the streams in two-dimensional models. The streams supply material to a circumplanetary disc that rotates in the same sense as the planet's orbit. Much of the mass supply to the circumplanetary disc comes from non-coplanar flow. The accretion rate peaks with a planet mass of approximately 0.1 M$_J$ and is highly efficient, occurring at the local viscous rate. The migration timescales for planets of mass less than 0.1 M$_J$, based on torques from disc material outside the planets' Roche lobes, are in excellent agreement with the linear theory of Type I (non-gap) migration for three-dimensional discs. The transition from Type I to Type II (gap) migration is smooth, with changes in migration times of about a factor of 2. Starting with a core which can undergo runaway growth, a planet can gain up to a few M$_J$ with little migration. Planets with final masses of order 10 M$_J$ would undergo large migration, which makes formation and survival difficult.

We present a full hydrodynamical investigation of a system of self-contained, interacting gas `cloudlets'. Calculations are performed using SPH, and cloudlets are allowed to collide, dissipate kinetic energy, merge and undergo local gravitational collapse. The numerical feasibility of such a technique is explored, and the resolution requirements examined. We examine the effect of the collision time and the velocity field on the timescale for evolution of the system, and show that almost all star formation will be confined to dense regions much smaller than the entire cloud. We also investigate the possibility, discussed by various authors, that such a system could lead to a power-law IMF by a process of repeated cloudlet coagulation, finding that in fact the inter-cloudlet collisions occur at too high a Mach number for merging to play an important part.

We present results from the first hydrodynamical star formation calculation to demonstrate that close binary stellar systems (separations $\lsim 10$ AU) need not be formed directly by fragmentation. Instead, a high frequency of close binaries can be produced through a combination of dynamical interactions in unstable multiple systems and the orbital decay of initially wider binaries. Orbital decay may occur due to gas accretion and/or the interaction of a binary with its circumbinary disc. These three mechanisms avoid the problems associated with the fragmentation of optically-thick gas to form close systems directly. They also result in a preference for close binaries to have roughly equal-mass components because dynamical exchange interactions and the accretion of gas with high specific angular momentum drive mass ratios towards unity. Furthermore, due to the importance of dynamical interactions, we find that stars with greater masses ought to have a higher frequency of close companions, and that many close binaries ought to have wide companions. These properties are in good agreement with the results of observational surveys.

We present results from the largest numerical simulation of star formation to resolve the fragmentation process down to the opacity limit. The simulation follows the collapse and fragmentation of a large-scale turbulent molecular cloud to form a stellar cluster and, simultaneously, the formation of circumstellar discs and binary stars. This large range of scales enables us to predict a wide variety of stellar properties for comparison with observations. The calculation clearly demonstrates that star formation is a highly-dynamic and chaotic process. Star-disc encounters form binaries and truncate discs. Stellar encounters disrupt bound multiple systems. The cloud produces roughly equal numbers of stars and brown dwarfs, with masses down to the opacity limit for fragmentation (~5 Jupiter masses). The initial mass function is consistent with a Salpeter slope (Gamma=-1.35) above 0.5 Msun, a roughly flat distribution (Gamma=0) in the range 0.006-0.5 Msun, and a sharp cutoff below ~0.005 Msun. This is consistent with recent observational surveys. The brown dwarfs form by the dynamical ejection of low-mass fragments from dynamically unstable multiple systems before the fragments have been able to accrete to stellar masses. Close binary systems (with separations <10 AU) are not formed by fragmentation in situ. Rather, they are produced by hardening of initially wider multiple systems through a combination of dynamical encounters, gas accretion, and/or the interaction with circumbinary and circumtriple discs. Finally, we find that the majority of circumstellar discs have radii less than 20 AU due to truncation in dynamical encounters. This is consistent with observations of the Orion Trapezium Cluster and implies that most stars and brown dwarfs do not form large planetary systems.

Using a local model Gammie (2001) has shown that accretion discs with cooling times t_cool <= 3/Omega fragment into gravitationally bound objects, while those with cooling times t_cool > 3/Omega evolve into a quasi-steady state. We use three-dimensional smoothed particle hydrodynamic simulations of protoplanetary accretion discs to test if the local results hold globally. We find that for disc masses appropriate for T Tauri discs, the fragmentation boundary still occurs at a cooling time close to t_cool = 3/Omega. For more massive discs, which are likely to be present at an earlier stage of the star formation process, fragmentation occurs for longer cooling times, but still within a factor of two of that predicted using a local model. These results have implications not only for planet formation in protoplanetary discs and star formation in AGN discs, but also for the redistribution of angular momentum which could be driven by the presence of relatively massive objects within the accretion disc.

We present three-dimensional numerical simulations of the interaction of a circular-orbit planet with a protoplanetary disk. We calculate the flow pattern, the accretion rate, and torques on the planet. We consider planet masses ranging from 1 Earth mass to 1 Jupiter mass.

We use high resolution numerical simulations to study whether gravitational instabilities within circumstellar discs can produce astrometrically detectable motion of the central star. For discs with masses of M_disc = 0.1 M_star, which are permanantly stable against fragmentation, we find that the magnitude of the astrometric signal depends upon the efficiency of disc cooling. Short cooling times produce prominent filamentary spiral structures in the disc, and lead to stellar motions that are potentially observable with future high precision astrometric experiments. For a disc that is marginally unstable within radii of \~10 au, we estimate astrometric displacements of 10-100 microarcsec on decade timescales for a star at a distance of 100 pc. The predicted displacement is suppressed by a factor of several in more stable discs in which the cooling time exceeds the local dynamical time by an order of magnitude. We find that the largest contribution comes from material in the outer regions of the disc and hence, in the most pessimistic scenario, the stellar motions caused by the disc could confuse astrometric searches for low mass planets orbiting at large radii. They are, however, unlikely to present any complications in searches for embedded planets orbiting at small radii, relative to the disc size, or Jupiter mass planets or greater orbiting at large radii.

We study the viscous effects on the interaction between the coplanar Be-star disc and the neutron star in Be/X-ray binaries, using a three-dimensional, smoothed particle hydrodynamics code. For simplicity, we assume the Be disc to be isothermal at the temperature of half the stellar effective temperature. In order to mimic the gas ejection process from the Be star, we inject particles with the Keplerian rotation velocity at a radius just outside the star. Both Be star and neutron star are treated as point masses. We find that the Be-star disc is effectively truncated if the Shakura-Sunyaev viscosity parameter alpha_SS >> 1, which confirms the previous semi-analytical result. In the truncated disc, the material decreted from the Be star accumulates, so that the disc becomes denser more rapidly than if around an isolated Be star. The resonant truncation of the Be disc results in a significant reduction of the amount of gas captured by the neutron star and a strong dependence of the mass capture rate on the orbital phase. We also find that an eccentric mode is excited in the Be disc through direct driving due to a one-armed bar potential of the binary. The strength of the mode becomes greater in the case of a smaller viscosity. In a high-resolution simulation with alpha_SS=0.1, the eccentric mode is found to precess in a prograde sense. The mass capture rate by the neutron star modulates as the mode precesses.

We present results from the first hydrodynamical star formation calculation to demonstrate that brown dwarfs are a natural and frequent product of the collapse and fragmentation of a turbulent molecular cloud. The brown dwarfs form via the fragmentation of dense molecular gas in unstable multiple systems and are ejected from the dense gas before they have been able to accrete to stellar masses. Thus, they can be viewed as `failed stars'. Approximately three quarters of the brown dwarfs form in gravitationally-unstable circumstellar discs while the remainder form in collapsing filaments of molecular gas. These formation mechanisms are very efficient, producing roughly the same number of brown dwarfs as stars, in agreement with recent observations. However, because close dynamical interactions are involved in their formation, we find a very low frequency of binary brown dwarf systems ($\lsim 5$%) and that those binary brown dwarf systems that do exist must be close $\lsim 10$ AU. Similarly, we find that young brown dwarfs with large circumstellar discs (radii $\gsim 10$ AU) are rare ($\approx 5$%).

We analyse the non-linear propagation and dissipation of axisymmetric waves in accretion discs using the ZEUS-2D hydrodynamics code. The waves are numerically resolved in the vertical and radial directions. Both vertically isothermal and thermally stratified accretion discs are considered. The waves are generated by means of resonant forcing and several forms of forcing are considered. Compressional motions are taken to be locally adiabatic ($\gamma = 5/3$). Prior to non-linear dissipation, the numerical results are in excellent agreement with the linear theory of wave channelling in predicting the types of modes that are excited, the energy flux by carried by each mode, and the vertical wave energy distribution as a function of radius. In all cases, waves are excited that propagate on both sides of the resonance (inwards and outwards). For vertically isothermal discs, non-linear dissipation occurs primarily through shocks that result from the classical steepening of acoustic waves. For discs that are substantially thermally stratified, wave channelling is the primary mechanism for shock generation. Wave channelling boosts the Mach number of the wave by vertically confining the wave to a small cool region at the base of the disc atmosphere. In general, outwardly propagating waves with Mach numbers near resonance ${\cal M}_{\rm r} \ga 0.01$ undergo shocks within a distance of order the resonance radius.

We performed high-resolution simulations of two stellar collisions relevant for stars in globular clusters. We considered one head-on collision and one off-axis collision between two 0.6 M_sun main sequence stars. We show that a resolution of about 100 000 particles is sufficient for most studies of the structure and evolution of blue stragglers. We demonstrate conclusively that collision products between main-sequence stars in globular clusters do not have surface convection zones larger than 0.004 M_sun after the collision, nor do they develop convection zones during the `pre-main-sequence' thermal relaxation phase of their post-collision evolution. Therefore, any mechanism which requires a surface convection zone (i.e. chemical mixing or angular momentum loss via a magnetic wind) cannot operate in these stars. We show that no disk of material surrounding the collision product is produced in off-axis collisions. The lack of both a convection zone and a disk proves a continuing problem for the angular momentum evolution of blue stragglers in globular clusters.

We present a simple physical mechanism that can account for the observed stellar mass spectrum for masses $\ms \simgreat 0.5 \solm$. The model depends solely on the competitive accretion that occurs in stellar clusters where each star's accretion rate depends on the local gas density and the square of the accretion radius. In a stellar cluster, there are two different regimes depending on whether the gas or the stars dominate the gravitational potential. When the cluster is dominated by cold gas, the accretion radius is given by a tidal-lobe radius. This occurs as the cluster collapses towards a $\rho\propto R^{-2}$ distribution. Accretion in this regime results in a mass spectrum with an asymptotic limit of $\gamma=-3/2$ (where Salpeter is $\gamma=-2.35$). Once the stars dominate the potential and are virialised, which occurs first in the cluster core, the accretion radius is the Bondi-Hoyle radius. The resultant mass spectrum has an asymptotic limit of $\gamma=-2$ with slightly steeper slopes ($\gamma\approx-2.5$) if the stars are already mass-segr egated. Simulations of accretion onto clusters containing 1000 stars show that as expected, the low-mass stars accumulate the majority of their masses during the gas dominated phase whereas the high-mass stars accumulate the majority of their massed during the stellar dominated phase. This results in a mass spectrum with a relatively shallow $\gamma\approx 3/2$ power-law for low-mass stars and a steeper, power-law for high-mass stars $ -2.5\simless\gamma\le -2$. This competitive accretion model also results in a mass segregated cluster.

We investigate the physics of gas accretion in young stellar clusters. Accretion in clusters is a dynamic phenomenon as both the stars and the gas respond to the same gravitational potential. Accretion rates are highly non-uniform with stars nearer the centre of the cluster, where gas densities are higher, accreting more than others. This competitive accretion naturally results in both initial mass segregation and a spectrum of stellar masses. Accretion in gas-dominated clusters is well modelled using a tidal-lobe radius instead of the commonly used Bondi-Hoyle accretion radius. This works as both the stellar and gas velocities are under the influence of the same gravitational potential and are thus comparable. The low relative velocity that results means that the tidal radius is smaller than the Bondi-Hoyle radius in these systems. In contrast, when the stars dominate the potential and are virialised, the Bondi-Hoyle radius is smaller than the tidal radius and thus Bondi-Hoyle accretion is a better fit to the accretion rates.

We consider the dynamics of a protostellar disc in a binary system where the disc is misaligned with the orbital plane of the binary, with the aim of determining the observational consequences for such systems. The disc wobbles with a period approximately equal to half the binary's orbital period and precesses on a longer timescale. We determine the characteristic timescale for realignment of the disc with the orbital plane due to dissipation. If the dissipation is determined by a simple isotropic viscosity then we find, in line with previous studies, that the alignment timescale is of order the viscous evolution timescale. However, for typical protostellar disc parameters, if the disc tilt exceeds the opening angle of the disc, then tidally induced shearing within the disc is transonic. In general, hydrodynamic instabilities associated with the internally driven shear result in extra dissipation which is expected to drastically reduce the alignment timescale. For large disc tilts the alignment timescale is then comparable to the precession timescale, while for smaller tilt angles $\delta$, the alignment timescale varies as $(\sin \delta)^{-1}$. We discuss the consequences of the wobbling, precession and rapid realignment for observations of protostellar jets and the implications for binary star formation mechanisms.

We review recent advances in the study of binary star formation that have been made using the smoothed particle hydrodynamics (SPH) technique. We discuss the importance of ensuring that the Jeans mass is always resolved during a self-gravitating SPH calculation. Even obeying this criterion, we show that it is now possible to perform calculations which follow the collapse of a molecular cloud core all the way to stellar densities. These calculations resolve 7 orders of magnitude in spatial extent and 17 orders of magnitude in density contrast. Finally, we discuss the progress that has been made in understanding the mass-ratio distributions of binary stellar systems by studying the evolution of protobinary systems as they accrete gas from their surrounding envelopes.

We investigate the formation of binary stellar systems. We consider a model where a `seed' protobinary system forms, via fragmentation, within a collapsing molecular cloud core and evolves to its final mass by accreting material from an infalling gaseous envelope. This accretion alters the mass ratio and orbit of the binary, and is largely responsible for forming the circumstellar and/or circumbinary discs. Given this model for binary formation, we predict the properties of binary systems and how they depend on the initial conditions within the molecular cloud core. We predict that there should be a continuous trend such that closer binaries are more likely to have equal mass components and are more likely to have circumbinary discs than wider systems. Comparing our results to observations, we find that the observed mass-ratio distributions of binaries and the frequency of circumbinary discs as a function of separation are most easily reproduced if the progenitor molecular cloud cores have radial density profiles between uniform and 1/r (e.g. Gaussian) with near uniform-rotation. This is in good agreement with the observed properties of pre-stellar cores. Conversely, we find that the observed properties of binaries cannot be reproduced if the cloud cores are in solid-body rotation and have initial density profiles which are strongly centrally condensed. Finally, in agreement with the radial-velocity searches for extra-solar planets, we find that it is very difficult to form a brown dwarf companion to a solar-type star with a separation < 10 AU, but that the frequency of brown dwarf companions should increase with larger separations or lower mass primaries.

We present results from the first three-dimensional calculations ever to follow the collapse of a molecular cloud core (~ 10^{-18} g cm^{-3}) to stellar densities (> 0.01 g cm^{-3}). The calculations resolve structures over 7 orders of magnitude in spatial extent (~ 5000 AU - 0.1 R_\odot), and over 17 orders of magnitude in density contrast. With these calculations, we consider whether fragmentation to form a close binary stellar system can occur during the second collapse phase. We find that, if the quasistatic core that forms before the second collapse phase is dynamically unstable to the growth of non-axisymmetric perturbations, the angular momentum extracted from the central regions of the core, via gravitational torques, is sufficient to prevent fragmentation and the formation of a close binary during the subsequent second collapse.

The isothermal gravitational collapse and fragmentation of a region within a molecular cloud and the subsequent formation of a protostellar cluster is investigated numerically. The clump mass spectrum which forms during the fragmentation phase can be well approximated by a power law distribution dN/dM ~ M^(-1.5). In contrast, the mass spectrum of protostellar cores that form in the centers of Jeans-unstable clumps and evolve through accretion and N-body interactions is described by a log-normal distribution with a width that is in excellent agreement with observations of multiple stellar systems.

We study the interpretation of the mean surface density of stellar companions as a function of separation (or the two point correlation function of stars) in star-forming regions. First, we consider the form of the functions for various simple stellar distributions (binaries, global density profiles, clusters, and fractals) and the effects of survey boundaries. Following this, we study the dependencies of the separation at which a transition from the binary to the large-scale clustering regime occurs. Larson (1995) found that the mean surface density of companions follows different power-law functions of separation in the two regimes. He identified the transition separation with the typical Jeans length in the molecular cloud. However, we show that this is valid only for special cases. In general, the transition separation depends on the volume density of stars, the depth of the star-forming region, the volume-filling nature of the stellar distribution, and on the parameters of the binaries. Furthermore, the transition separation evolves with time. We then apply these results to the Taurus-Auriga, Ophiuchus, and Orion Trapezium star-forming regions. We find that while the transition separation in the Taurus-Auriga star-forming region may indicate a typical Jeans length, this is not true of the Orion Trapezium Cluster. We also show that there is no evidence for sub-clustering or fractal structure in the stars of the Orion Trapezium Cluster. This is consistent with the fact that, if such structure were present when the stars formed, it would have been erased by the current age of the cluster due to the stellar velocity dispersion. [Abstract abridged]

We present a model for the formation of massive ($M > 10 M_\odot$) stars through accretion-induced collisions in the cores of embedded dense stellar clusters. This model circumvents the problem of accreting onto a star whose luminosity is sufficient to reverse the infall of gas. Instead, the central core of the cluster accretes from the surrounding gas, thereby decreasing its radius until collisions between individual components become significant. These components are, in general, intermediate-mass stars that have formed through accretion onto low-mass protostars. Once a sufficiently massive star has formed to expel the remaining gas, the cluster expands in accordance with this loss of mass, halting further collisions. This process implies a critical stellar density for the formation of massive stars, and a high rate of binaries formed by tidal capture.

Hydrodynamical calculations in three space dimensions of the collapse of an isothermal, rotating 1 M\sol protostellar cloud are presented. The initial density stratification is a power law with density $\rho \propto r^{-p}$, with $p=1$. The case of the singular isothermal sphere ($p=2$) is not considered; however $p=1$ has been shown observationally to be a good representation of the density distribution in molecular cloud cores just before the beginning of collapse. The collapse is studied with two independent numerical methods, an SPH code with 200,000 particles, and a finite-difference code with nested grids which give high spatial resolution in the inner regions. Although previous numerical studies have indicated that such a power-law distribution would not result in fragmentation into a binary system, both codes show, in contrast, that multiple fragmentation does occur in the central regions of the protostar. Thus the process of binary formation by fragmentation is shown to be consistent with the fact that a large fraction of young stars are observed to be in binary or multiple systems.

A method for following fragmentation simulations further in time using smoothed particle hydrodynamics (SPH) is presented. In a normal SPH simulation of the collapse and fragmentation of a molecular cloud, high-density regions of gas that form protostars are represented by many particles with small separations. These high-density regions require small time steps, limiting the time for which the simulation can be followed. Thus, the end result of the fragmentation can never be definitively ascertained, and comparisons between cloud fragmentation calculations and the observed characteristics of stellar systems cannot be made. In this paper, each high-density region is replaced by a single, non-gaseous particle, with appropriate boundary conditions, which contains all the mass in the region and accretes any infalling mass. This enables the evolution of the cloud and the resulting protostars to be followed for many orbits or until most of the original cloud mass has been accreted.

A viable solution to the origin of close binary systems, unaccounted for in recent theories, is presented. Fragmentation, occurring at the end of the secondary collapse phase (during which molecular hydrogen is dissociating), can form binary systems with separations less than 1 au. Two fragmentation modes are found to occur after the collapse is halted. The first consists of the fragmentation of a protostellar disc due to rotational instabilities in a protostellar core, involving both an $m=1$ and an $m=2$ mode. This fragmentation mechanism is found to be insensitive to the initial density distribution: it can occur in both centrally condensed and uniform initial conditions. The second fragmentation mode involves the formation of a rapidly rotating core at the end of the collapse phase which is unstable to the axisymmetric perturbations. This core bounces into a ring which quickly fragments into several components. The binary systems thus formed contain less than 1 per cent of a solar mass and therefore will need to accrete most of their final mass if they are to form a binary star system. Their orbital properties will thus be determined by the properties of the accreted matter.