(abridged) The nature of early Class 0/I protostellar discs is not clearly understood. Early protostellar discs are needed to drive molecular outflows and jets observed in star forming regions, but there has been some debate to how they form. From a theoretical perspective, the consequences of disc and outflow generation are crucial to understanding the very nature of how stars are assembled. We have performed 3D ideal magnetohydrodynamic (MHD) simulations of collapsing Bonnor-Ebert spheres, employing sink particles with a radius of 3.2 AU alongside an AMR grid and using a cooling function to model radiative cooling of the gas. This has allowed us to explore 2-8x10^4 yr further into the evolution of an early Class 0 disc-outflow system than previous simulations. Our outflow is precessing, kinked, turbulent, contains internal shocks and has a scale of 0.1 pc end-to-end. We form a rotationally dominated disc with a radius of 100 AU embedded inside a transient, unstable, flattened, rotating core extending out to 2000 AU. The larger flattened structure launches a low speed wind (v_r < 1.5 km/s) dominated by B_phi, while the inner disc launches a centrifugally driven jet dominated by B_p with speeds up to 20$ km/s. From the inner disk, the value of dM/dt_out/dM/dt_in ~ 0.1, wheras in the outer core dM/dt_out/dM/dt_in ~ 1.0. The inner disc becomes unstable to a warping instability due to the magnetic structure of the outflow and warps to 30 deg with respect to the z-axis by the end of the simulation. The envelope is cleared out and is less massive than the disc. We measure star formation efficiencies of eta_core=0.63 (and growing), higher than theoretical predictions. This indicates that outflows are not as efficient at expelling envelope mass as some current models estimate. We discuss the relevance of our disc misalignment concerning the formation of mis-aligned hot Jupiters.

Numerical simulations of outflows formed during the collapse of 100 M_sun cloud cores are presented. We derive a generalised criterion from MHD wind theory to analyse the launching mechanism of these outflows. The criterion is successfully applied to the whole outflow structure and cases with sub-Keplerian disc rotation. It allows us to decide whether an outflow is driven centrifugally or by the toroidal magnetic pressure. We show that quantities such as the magnetic field line inclination or the ratio of the toroidal to poloidal magnetic field alone are insufficient to determine the driving mechanism of outflows. By performing 12 runs we are able to study the influence of the initial conditions on the properties of outflows around massive protostars in detail. Our simulations reveal a strong effect of the magnetic field strength on the outflow morphology. In runs with weak fields or high rotational energies, well collimated, fast jets are observed whereas for strong fields poorly collimated, low-velocity outflows are found. We show that the occurrence of a fast jet is coupled to the existence of a Keplerian protostellar disc. Despite the very different morphologies all outflows are launched from the discs by centrifugal acceleration with the toroidal magnetic field increasingly contributing to the gas acceleration further away from the discs. The poor collimation of the outflows in runs with strong fields is a consequence of the weak hoop stresses. This in turn is caused by the slow build-up of a toroidal field due to sub-Keplerian disc rotation. The mass ejection/accretion ratios scatter around a mean of 0.3 in accordance with observational and analytical results. We suggest an evolutionary scenario for the earliest stage of massive star formation in which initially poorly collimated outflows develop which progressively get better collimated due to the generation of fast jets.

We present simulations of collapsing 100 M_\sun mass cores in the context of massive star formation. The effect of variable initial rotational and magnetic energies on the formation of massive stars is studied in detail. We focus on accretion rates and on the question under which conditions massive Keplerian disks can form in the very early evolutionary stage of massive protostars. For this purpose, we perform 12 simulations with different initial conditions extending over a wide range in parameter space. The equations of magnetohydrodynamics (MHD) are solved under the assumption of ideal MHD. We find that the formation of Keplerian disks in the very early stages is suppressed for a mass-to-flux ratio normalised to the critical value \mu below 10, in agreement with a series of low-mass star formation simulations. This is caused by very efficient magnetic braking resulting in a nearly instantaneous removal of angular momentum from the disk. For weak magnetic fields, corresponding to \mu > 10, large-scale, centrifugally supported disks build up with radii exceeding 100 AU. A stability analysis reveals that the disks are supported against gravitationally induced perturbations by the magnetic field and tend to form single stars rather than multiple objects. We find protostellar accretion rates of the order of a few 10^-4 M_\sun yr^-1 which, considering the large range covered by the initial conditions, vary only by a factor of ~ 3 between the different simulations. We attribute this fact to two competing effects of magnetic fields. On the one hand, magnetic braking enhances accretion by removing angular momentum from the disk thus lowering the centrifugal support against gravity. On the other hand, the combined effect of magnetic pressure and magnetic tension counteracts gravity by exerting an outward directed force on the gas in the disk thus reducing the accretion onto the protostars.

We investigate the formation and evolution of giant molecular clouds (GMCs) by the collision of convergent warm neutral medium (WNM) streams in the interstellar medium, in the presence of magnetic fields and ambipolar diffusion (AD), focusing on the evolution of the star formation rate (SFR) and efficiency (SFE), as well as of the mass-to-magnetic-flux ratio (M2FR) in the forming clouds. We find that: 1) Clouds formed by supercritical inflow streams proceed directly to collapse, while clouds formed by subcritical streams first contract and then re-expand, oscillating on the scale of tens of Myr. 2) Our suite of simulations with initial magnetic field strength of 2, 3, and 4 $\mu\G$ show that only supercritical or marginal critical streams lead to reasonable star forming rates. 3) The GMC's M2FR is a generally increasing function of time, whose growth rate depends on the details of how mass is added to the GMC from the WNM. 4) The M2FR is a highly fluctuating function of position in the clouds. 5) In our simulations, the SFE approaches stationarity, because mass is added to the GMC at a similar rate at which it converts mass to stars. In such an approximately stationary regime, the SFE provides a proxy of the supercritical mass fraction in the cloud. 6) We observe the occurrence of buoyancy of the low-M2FR regions within the gravitationally-contracting GMCs, so that the latter naturally segregate into a high-density, high-M2FR "core" and a low-density, low-M2FR "envelope", without the intervention of AD. (Abridged)

The pre-stellar cores in which low mass stars form are generally well magnetized. Our simulations show that early protostellar discs are massive and experience strong magnetic torques in the form of magnetic braking and protostellar outflows. Simulations of protostellar disk formation suggest that these torques are strong enough to suppress a rotationally supported structure from forming for near critical values of mass-to-flux. We demonstrate through the use of a 3D adaptive mesh refinement code -- including cooling, sink particles and magnetic fields -- that one produces transient 1000 AU discs while simultaneously generating large outflows which leave the core region, carrying away mass and angular momentum. Early inflow/outflow rates suggest that only a small fraction of the mass is lost in the initial magnetic tower/jet event.

In star formation, magnetic fields act as a cosmic angular momentum extractor that increases mass accretion rates onto protostars and in the process, creates spectacular outflows. However, recently it has been argued that this magnetic brake is so strong that early protostellar disks -- the cradles of planet formation -- cannot form. Our three-dimensional numerical simulations of the early stages of collapse (\lesssim 10^5 yr) of overdense star--forming clouds form early outflows and have magnetically regulated and rotationally dominated disks (inside 10 AU) with high accretion rates, despite the slip of the field through the mostly neutral gas. We find that in three dimensions, magnetic fields suppress gravitationally driven instabilities which would otherwise prevent young, well ordered disks from forming. Our simulations have surprising consequences for the early formation of disks, their density and temperature structure, the mechanism and structure of early outflows, the flash heating of dust grains through ambipolar diffusion, and the origin of planets and binary stars.

Given the importance of simulating hydromagnetic processes that impact star formation, we have earlier developed a 3D adaptive mesh approach that allows us to include hydromagnetic processes during the formation and evolution of cores, discs, and stars in observed regions of star formation. In this paper, we take the next step in this program - namely - to develop a modified version of the 3D adaptive mesh refinement (AMR) code FLASH in which the ambipolar diffusion of the magnetic field in poorly ionized molecular gas is implemented. We approach the problem using a single-fluid approximation to simplify numerical calculations. In this paper, we present a series of test cases including oblique isothermal and non--isothermal C-shocks. We also present a study of the quasi-static collapse of an initial uniform, self-gravitating, magnetized sphere that is initially supported by its magnetic field against collapse (i.e. magnetically subcritical). Applications to the collapse of a pre-stellar Bonnor-Ebert sphere are presented in a companion paper.