Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of protostars by following the gravitational collapse of 1~M$_\odot$ gas clouds through the first hydrostatic core to stellar densities. The clouds are imposed with both rotational and turbulent velocities, and are threaded with a magnetic field that is parallel/anti-parallel or perpendicular to the rotation axis; we investigate two rotation rates and four Mach numbers. The initial radius and mass of the stellar core are only weakly dependent on the initial parameters. In the models that include ideal MHD, the magnetic field strength implanted in the protostar at birth is much higher than observed, independent of the initial level of turbulence; only non-ideal MHD can reduce this strength to near or below the observed levels. This suggests that not only is ideal MHD an incomplete picture of star formation, but that the magnetic fields in low mass stars are implanted later in life by a dynamo process. Non-ideal MHD suppresses magnetically launched stellar core outflows, but turbulence permits thermally launched outflows to form a few years after stellar core formation.

Non-ideal magnetohydrodynamics (MHD) is the dominant process. We investigate the effect of magnetic fields (ideal and non-ideal) and turbulence (sub- and transsonic) on the formation of circumstellar discs that form nearly simultaneously with the formation of the protostar. This is done by modelling the gravitational collapse of a 1~M$_\odot$ gas cloud that is threaded with a magnetic field and imposed with both rotational and turbulent velocities. We investigate magnetic fields that are parallel/anti-parallel and perpendicular to the rotation axis, two rotation rates and four Mach numbers. Disc formation occurs preferentially in the models that include non-ideal MHD where the magnetic field is anti-parallel or perpendicular to the rotation axis. This is independent of the initial rotation rate and level of turbulence, suggesting that subsonic turbulence plays a minimal role in influencing the formation of discs. Aside from first core outflows which are influenced by the initial level of turbulence, non-ideal MHD processes are more important than turbulent processes during the formation of discs around low-mass stars.

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.

We investigate differences in the molecular abundances between magnetically super- and sub-critical prestellar cores, performing three-dimensional non-ideal magnetohydrodynamical (MHD) simulations with varying densities and magnetic field strengths, and post-processing the results with a time-dependent gas-grain chemical code. Most molecular species show significantly more central depletion in subcritical models, due to the longer duration of collapse. However, the directly observable quantities - the molecule to hydrogen column density ratios - are generally too similar for observational data to discriminate between models. The profiles of N$_2$H$^+$ and HCO$^+$ show qualitative differences between supercritical and subcritical models on scales of $0.01 \, {\rm pc}$, which may allow the two cases to be distinguished. However, this requires knowledge of the hydrogen column density, which is not directly measureable, and predicted line intensity profiles from radiative transfer modelling are similar for these molecules. Other commonly observed species, such as HCN and CH$_3$OH, have line intensity profiles which differ more strongly between models, and so are more promising as tracers of the mechanism of cloud collapse.

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.

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.

Truncated abstract: The formation of a protostellar disc is a natural outcome during the star formation process. As gas in a molecular cloud core collapses under self-gravity, the angular momentum of the gas will slow its collapse on small scales and promote the formation of a protostellar disc. Although the angular momenta of dense star-forming cores remain to be fully characterized observationally, existing data indicates that typical cores have enough angular momenta to form relatively large, rotationally supported discs. However, molecular clouds are observed to be permeated by magnetic fields, which can strongly affect the evolution of angular momentum through magnetic braking. Indeed, in the ideal MHD limit, magnetic braking has been shown to be so efficient as to remove essentially all of the angular momentum of the material close to the forming star such that disc formation is suppressed. This is known as the magnetic braking catastrophe. The catastrophe must be averted in order for the all-important rotationally supported discs to appear, but when and how this happens remains debated. We review the resolutions proposed to date, with emphasis on misalignment, turbulence and especially non-ideal effects. The dissipative non-ideal effects weaken the magnetic field, and the dispersive term redirects it to promote or hinder disc formation. When self-consistently applying non-ideal processes, rotationally supported discs of at least tens of au form, thus preventing the magnetic braking catastrophe. The non-ideal processes are sensitive to the magnetic field strength, cosmic ray ionization rate, and gas and dust grain properties, thus a complete understanding of the host molecular cloud is required. Therefore, the properties of the host molecular cloud -- and especially its magnetic field -- cannot be ignored when numerically modelling the formation and evolution of protostellar discs.

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 Phantom, a fast, parallel, modular and low-memory smoothed particle hydrodynamics and magnetohydrodynamics code developed over the last decade for astrophysical applications in three dimensions. The code has been developed with a focus on stellar, galactic, planetary and high energy astrophysics and has already been used widely for studies of accretion discs and turbulence, from the birth of planets to how black holes accrete. Here we describe and test the core algorithms as well as modules for magnetohydrodynamics, self-gravity, sink particles, H_2 chemistry, dust-gas mixtures, physical viscosity, external forces including numerous galactic potentials as well as implementations of Lense-Thirring precession, Poynting-Robertson drag and stochastic turbulent driving. Phantom is hereby made publicly available.

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

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

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.

We present hydrodynamic simulations of the common envelope binary interaction between a giant star and a compact companion carried out with the adaptive mesh refinement code ENZO and the smooth particle hydrodynamics code PHANTOM. These simulations mimic the parameters of one of the simulations by Passy et al., but assess the impact of a larger, more realistic initial orbital separation on the simulation outcome. We conclude that for both codes the post-common envelope separation is somewhat larger and the amount of unbound mass slightly greater when the initial separation is wide enough that the giant does not yet overflow or just overflows its Roche lobe. PHANTOM has been adapted to the common envelope problem here for the first time and a full comparison with ENZO is presented, including an investigation of convergence as well as energy and angular momentum conservation. We also set our simulations in the context of past simulations. This comparison reveals that it is the expansion of the giant before rapid in-spiral and not spinning up of the star that causes a larger final separation. We also suggest that the large range in unbound mass for different simulations is difficult to explain and may have something to do with simulations that are not fully converged.

In this paper, we introduce Nicil: Non-Ideal magnetohydrodynamics Coefficients and Ionisation Library. Nicil is a stand-alone Fortran90 module that calculates the ionisation values and the coefficients of the non-ideal magnetohydrodynamics terms of Ohmic resistivity, the Hall effect, and ambipolar diffusion. The module is fully parameterised such that the user can decide which processes to include and decide upon the values of the free parameters, making this a versatile and customisable code. The module includes both cosmic ray and thermal ionisation; the former includes two ion species and three species of dust grains (positively charged, negatively charged and neutral), and the latter includes five elements which can be doubly ionised. We demonstrate tests of the module, and then describe how to implement it into an existing numerical code.

We perform adaptive mesh refinement (AMR) and smoothed particle hydrodynamics (SPH) cosmological zoom simulations of a region around a forming galaxy cluster, comparing the ability of the methods to handle successively more complex baryonic physics. In the simplest, non-radiative case, the two methods are in good agreement with each other, but the SPH simulations generate central cores with slightly lower entropies and virial shocks at slightly larger radii, consistent with what has been seen in previous studies. The inclusion of radiative cooling, star formation, and stellar feedback leads to much larger differences between the two methods. Most dramatically, at z=5, rapid cooling in the AMR case moves the accretion shock well within the virial radius, while this shock remains near the virial radius in the SPH case, due to excess heating, coupled with poorer capturing of the shock width. On the other hand, the addition of feedback from active galactic nuclei (AGN) to the simulations results in much better agreement between the methods. In this case both simulations display halo gas entropies of 100 keV cm^2, similar decrements in the star-formation rate, and a drop in the halo baryon content of roughly 30%. This is consistent with AGN growth being self-regulated, regardless of the numerical method. However, the simulations with AGN feedback continue to differ in aspects that are not self-regulated, such that in SPH a larger volume of gas is impacted by feedback, and the cluster still has a lower entropy central core.

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 partially ionised plasmas, the magnetic field can become decoupled from the neutral gas and diffuse through it in a process known as ambipolar diffusion. Although ambipolar diffusion has been implemented in several grid codes, we here provide an implementation in smoothed particle magnetohydrodynamics (SPMHD). We use the strong coupling approximation in which the ion density is negligible, allowing a single fluid approach. The equations are derived to conserve energy, and to provide a positive definite contribution to the entropy. We test the implementation in both a simple 1D SPMHD code and the fully 3D code PHANTOM. The wave damping test yields agreement within 0.03-2 per cent of the analytical result, depending on the value of the collisional coupling constant. The oblique C-shocks test yields results that typically agree within 4 per cent of the semi-analytical result. Our algorithm is therefore suitable for exploring the effect ambipolar diffusion has on physical processes, such as the formation of stars from molecular clouds.

We examine the correlation between the star formation rate (SFR) and black hole accretion rate (BHAR) across a suite of different AGN feedback models, using the time evolution of a merger simulation. By considering three different stages of evolution, and a distinction between the nuclear and outer regions of star formation, we consider 63 different cases. Despite many of the feedback models fitting the M-\sigma\ relationship well, there are often distinct differences in the SFR-BHAR correlations, with close to linear trends only being present after the merger. Some of the models also show evolution in the SFR-BHAR parameter space that is at times directly across the long-term averaged SFR-BHAR correlation. This suggests that the observational SFR-BHAR correlation found for ensembles of galaxies is an approximate statistical trend, as suggested by Hickox et al. Decomposing the SFR into nuclear and outer components also highlights notable differences between models and there is only modest agreement with observational studies examining this in Seyfert galaxies. For the fraction of the black hole mass growth from the merger event relative to the final black hole mass, we find as much as a factor of three variation among models. This also translates into a similar variation in the post-starburst black hole mass growth. Overall, we find that while qualitative features are often similar amongst models, precise quantitative analysis shows there can be quite distinct differences.

We examine simulations of isolated galaxies to analyse the effects of localised feedback on the formation and evolution of molecular clouds. Feedback contributes to turbulence and the destruction of clouds, leading to a population of clouds that is younger, less massive, and with more retrograde rotation. We investigate the evolution of clouds as they interact with each other and the diffuse ISM, and determine that the role of cloud interactions differs strongly with the presence of feedback: in models without feedback, scattering events dramatically increase the retrograde fraction, but in models with feedback, mergers between clouds may slightly increase the prograde fraction. We also produce an estimate of the viscous time-scale due to cloud-cloud collisions, which increases with increasing strength of feedback (~20 Gyr vs ~10 Gyr), but is still much smaller than previous estimates (~1000 Gyr); although collisions become more frequent with feedback, less energy is lost in each collision than in the models without feedback.

Modelling AGN feedback in numerical simulations is both technically and theoretically challenging, with numerous approaches having been published in the literature. We present a study of five distinct approaches to modelling AGN feedback within gravitohydrodynamic simulations of major mergers of Milky Way-sized galaxies. To constrain differences to only be between AGN feedback models, all simulations start from the same initial conditions and use the same star formation algorithm. Most AGN feedback algorithms have five key aspects: black hole accretion rate, energy feedback rate and method, particle accretion algorithm, black hole advection algorithm, and black hole merger algorithm. All models follow different accretion histories, with accretion rates that differ by up to three orders of magnitude at any given time. We consider models with either thermal or kinetic feedback, with the associated energy deposited locally around the black hole. Each feedback algorithm modifies the gas properties near the black hole to different extents. The particle accretion algorithms usually (but not always) maintain good agreement between the mass accreted by \dot{M} dt and the mass of gas particles removed from the simulation. The black hole advection algorithms dampen inappropriate dragging of the black holes by two-body interactions. Advecting the black hole a limited distance based upon local mass distributions has many desirably properties. The black holes merge when given criteria are met, and we find a range of merger times for different criteria. Using the M_{BH}-\sigma relation as a diagnostic of the remnants yields three models that lie within the one-sigma scatter of the observed relation and two that fall below it. The wide variation in accretion behaviours of the models reinforces the fact that there remains much to be learnt about the evolution of galactic nuclei. (abridged)

A recent approach to simulating localized feedback from active galactic nuclei by Power et al. (2011) uses an accretion disc particle to represent both the black hole and its accretion disc. We have extrapolated and adapted this approach to simulations of Milky Way-sized galaxy mergers containing black holes and explored the impact of the various parameters in this model as well as its resolution dependence. The two key parameters in the model are an effective accretion radius, which determines the radius within which gas particles are added to the accretion disc, and a viscous time-scale which determines how long it takes for material in the accretion disc to accrete on to the black hole itself. We find that there is a limited range of permitted accretion radii and viscous time-scales, with unphysical results produced outside this range. For permitted model parameters, the nuclear regions of simulations with the same resolution follow similar evolutionary paths, producing final black hole masses that are consistent within a factor of two. When comparing the resolution dependence of the model, there is a trend towards higher resolution producing slightly lower mass black holes, but values for the two resolutions studied again agree within a factor of two. We also compare these results to two other AGN feedback algorithms found in the literature. While the evolution of the systems vary, most notably the intermediate total black hole mass, the final black hole masses differ by less than a factor of five amongst all of our models, and the remnants exhibit similar structural parameters. The implication of this accretion model is that, unlike most accretion algorithms, a decoupling of the accretion rate on to the black hole and the local gas properties is permitted and obtained; this allows for black hole growth even after feedback has prevented additional accretion events on to the disc.

We present a parameter study of simulations of fragmentation regulated by gravity, magnetic fields, ambipolar diffusion, and nonlinear flows. The thin-sheet approximation is employed with periodic lateral boundary conditions, and the nonlinear flow field ("turbulence") is allowed to freely decay. In agreement with previous results in the literature, our results show that the onset of runaway collapse (formation of the first star) in subcritical clouds is significantly accelerated by nonlinear flows in which a large-scale wave mode dominates the power spectrum. In addition, we find that a power spectrum with equal energy on all scales also accelerates collapse, but by a lesser amount. For a highly super-Alfvenic initial velocity field with most power on the largest scales, the collapse occurs promptly during the initial compression wave. However, for trans-Alfvenic perturbations, a subcritical magnetic field causes a rebound from the initial compression, and the system undergoes several oscillations before runaway collapse occurs. Models that undergo prompt collapse have highly supersonic infall motions at the core boundaries. Cores in magnetically subcritical models with trans-Alfvenic initial perturbations also pick up significant systematic speeds by inheriting motions associated with magnetically-driven oscillations. Core mass distributions are much broader than in models with small-amplitude initial perturbations, although the disturbed structure of cores that form due to nonlinear flows does not guarantee subsequent monolithic collapse. Our simulations also demonstrate that significant power can (if present initially) be maintained with negligible dissipation in large-scale compressive modes of a magnetic thin sheet, in the limit of perfect flux freezing.

We present results from an extensive set of simulations of gravitational fragmentation in the presence of magnetic fields and ambipolar diffusion. The average fragmentation spacing in the nonlinear phase of evolution is in excellent agreement with the prediction of linear perturbation theory. The time scale for nonlinear growth and runaway of the first core is $\approx 10$ times the calculated growth time $\taugm$ of the eigenmode with minimum growth time, when starting from a uniform background state with small-amplitude white-noise perturbations. Subcritical and transcritical models typically evolve on a significantly longer time scale than the supercritical models. Infall motions in the nonlinear fully-developed contracting cores are subsonic on the core scale in subcritical and transcritical clouds, but are somewhat supersonic in supercritical clouds. Core mass distributions are sharply peaked with a steep decline to large masses, consistent with the existence of a preferred mass scale for each unique set of dimensionless free parameters. However, a sum total of results for various initial mass-to-flux ratios yields a broad distribution reminiscent of observed core mass distributions. Based on our results, we conclude that fragmentation spacings, magnitude of infall motions, core shapes, and, especially, the curvature of magnetic field morphology, may serve as indirect observational means of determining a cloud's ambient mass-to-flux ratio.