We review recent results of non-ideal magnetohydrodynamic models for the fragmentation of molecular clouds and the collapse of cloud cores to form protostar-disk systems. Thin-disk models can elucidate many aspects of the physical problem and allow the calculation of large dynamic range of time and length scales.

Context: Ideal MHD simulations have revealed catastrophic magnetic braking (MB) in the protostellar phase, which prevents the formation of a centrifugal disk around a nascent protostar. Aims: We determine if non-ideal MHD, including the effects of ambipolar diffusion and Ohmic dissipation determined from a detailed chemical network model, allows for disk formation at the earliest stages of star formation (SF). Methods: We employ the axisymmetric thin-disk approximation in order to resolve a dynamic range of 9 orders of magnitude in length and 16 in density, while also calculating partial ionization using up to 19 species in a detailed chemical equilibrium model. MB is applied using a steady-state approximation, and a barotropic relation is used to capture the thermal evolution. Results: We resolve the formation of the first and second cores, with expansion waves at the periphery of each, a magnetic diffusion shock, and prestellar infall profiles at larger radii. Power-law profiles in each region can be understood analytically. After the formation of the second core, centrifugal support rises rapidly and a low-mass disk of radius ~10 R_Sun is formed, when the second core has mass ~0.001 M_Sun. The mass-to-flux ratio is ~10,000 times the critical value in the central region. Conclusions: A small centrifugal disk can form in the earliest stage of SF, due to a shut-off of MB caused by magnetic field dissipation in the first core region. There is enough angular momentum loss to allow the second collapse to occur directly, and a low-mass stellar core to form with a surrounding disk. The disk mass and size will depend upon how the angular momentum transport mechanisms within the disk can keep up with mass infall onto the disk. We estimate that the disk will remain <~10 AU, undetectable even by ALMA, in the early Class 0 phase.

A non-ideal MHD collapse calculation employing the axisymmetric thin-disk approximation is used to resolve cloud core collapse down to the scales of the second (stellar) core. Rotation and a magnetic braking torque are included in the model, and the partial ionization resulting in ambipolar diffusion and Ohmic dissipation is calculated from a detailed chemical network. We find that a centrifugal disk can indeed form in the earliest stage of star formation, due to a shut-off of magnetic braking caused by magnetic field diffusion in the first core region. Thus, there is no catastrophic magnetic braking in a model with realistic non-ideal MHD.

Magnetic diffusion plays a vital role in star formation. We trace its influence from interstellar cloud scales down to star-disk scales. On both scales, we find that magnetic diffusion can be significantly enhanced by the buildup of strong gradients in magnetic field structure. Large scale nonlinear flows can create compressed cloud layers within which ambipolar diffusion occurs rapidly. However, in the flux-freezing limit that may be applicable to photoionized molecular cloud envelopes, supersonic motions can persist for long times if driven by an externally generated magnetic field that corresponds to a subcritical mass-to-flux ratio. In the case of protostellar accretion, rapid magnetic diffusion (through Ohmic dissipation with additional support from ambipolar diffusion) near the protostar causes dramatic magnetic flux loss. By doing so, it also allows the formation of a centrifugal disk, thereby avoiding the magnetic braking catastrophe.

We perform axisymmetric resistive MHD calculations that demonstrate that centrifugal disks can indeed form around Class 0 objects despite magnetic braking. We follow the evolution of a prestellar core all the way to near-stellar densities and stellar radii. Under flux-freezing, the core is braked and disk formation is inhibited, while Ohmic dissipation renders magnetic braking ineffective within the first core. In agreement with observations that do not show evidence for large disks around Class 0 objects, the resultant disk forms in close proximity to the second core and has a radius of only $\approx 10~R_{\odot}$ early on. Disk formation does not require enhanced resistivity. We speculate that the disks can grow to the sizes observed around Class II stars over time under the influence of both Ohmic dissipation and ambipolar diffusion, as well as internal angular momentum redistribution.

We calculate and analyze the longevity of magnetohydrodynamic (MHD) wave modes that occur in the plane of a magnetic thin sheet. Initial turbulent conditions applied to a magnetically subcritical cloud are shown to lead to relatively rapid energy decay if ambipolar diffusion is introduced at a level corresponding to partial ionization primarily by cosmic rays. However, in the flux-freezing limit, as may be applicable to photoionized molecular cloud envelopes, the turbulence persists at "nonlinear" levels in comparison with the isothermal sound speed $\cs$, with one-dimensional rms material motions in the range of $\approx 2\,\cs -5\,\cs$ for cloud sizes in the range of $\approx 2\,\pc - 16\,\pc$. These fluctuations persist indefinitely, maintaining a significant portion of the initial turbulent kinetic energy. We find the analytic explanation for these persistent fluctuations. They are magnetic-tension-driven modes associated with the interaction of the sheet with the external magnetic field. The phase speed of such modes is quite large, allowing residual motions to persist without dissipation in the flux-freezing limit, even as they are nonlinear with respect to the sound speed. We speculate that long-lived large-scale MHD modes such as these may provide the key to understanding observed supersonic motions in molecular clouds.

We present a new analytical three-parameter formula to fit observed column density profiles of prestellar cores. It represents a line-of-sight integral through a spherically symmetric or disc-like isothermal cloud. The underlying model resembles the Bonnor-Ebert model in that it features a flat central region leading into a power-law decline \propto r^{-2} in density, and a well-defined outer radius. However, we do not assume that the cloud is in equilibrium, and can instead make qualitative statements about its dynamical state (expansion, equilibrium, collapse) using the size of the flat region as a proxy. Instead of having temperature as a fitting parameter, our model includes it as input, and thus avoids possible inconsistencies. It is significantly easier to fit to observational data than the Bonnor-Ebert sphere. We apply this model to L1689B and B68. We show that L1689B cannot be in equilibrium but instead appears to be collapsing, while our model verifies that B68 is not far from being a hydrostatic object.

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.