We revisit the relation between magnetic-field strength ($B$) and gas density ($\rho$) for contracting interstellar clouds and fragments (or, cores), which is central in observationally determining the dynamical importance of magnetic fields in cloud evolution and star formation. Recently, it has been claimed that a relation $B \propto \rho^{2/3} $ is statistically preferred over $B \propto \rho^{1/2}$ in molecular clouds, when magnetic field detections and nondetections from Zeeman observations are combined. This finding has unique observational implications on cloud and core geometry: The relation $B \propto \rho^{2/3} $ can only be realized under spherical contraction. However, no indication of spherical geometry can be found for the objects used in the original statistical analysis of the $B-\rho$ relation. We trace the origin of the inconsistency to simplifying assumptions in the statistical model used to arrive at the $B\propto \rho^{2/3}$ conclusion and to an underestimate of observational uncertainties in the determination of cloud and core densities. We show that, when these restrictive assumptions are relaxed, $B \propto \rho^{1/2}$ is the preferred relation for the (self-gravitating) molecular-cloud data, as theoretically predicted four decades ago.

We present a comprehensive study of MHD waves and instabilities in a weakly ionised system, e.g., an interstellar molecular cloud. We determine all the critical wavelengths of perturbations across which the sustainable wave modes can change radically (and so can their decay rates), and various instabilities are present or absent. These critical wavelengths are essential for understanding the effects of MHD waves (or turbulence) on the structure and evolution of molecular clouds. Depending on the angle of propagation relative to the magnetic field and the physical parameters of a cloud, there are wavelength ranges in which no wave can be sustained as such. Yet, for other directions of propagation or different properties of a model cloud, there may exist some wave mode(s) at all wavelengths. For a typical cloud, magnetically-driven ambipolar diffusion leads to removal of any support against gravity that most short-wavelength waves (or turbulence) may have had, and gravitationally-driven ambipolar diffusion sets in and leads to cloud fragmentation into stellar-size masses, as first suggested by Mouschovias three decades ago -- a single-stage fragmentation theory of star formation, distinct from the hierarchical fragmentation picture. Phase velocities, decay times, and eigenvectors (the densities and velocities of neutral particles and the plasma, and components of the magnetic field) are determined as functions of the wavelength of the disturbances and are explained physically. Comparison of the results with nonlinear analytical or numerical calculations is also presented, excellent agreement is found, and confidence in the analytical, linear approach is gained to explore phenomena difficult to study through numerical simulations. Mode splitting and mode merging, which are impossible in single-fluid systems for linear perturbations, occur naturally in multifluid systems.

We had recently re-analyzed in a self-consistent way OH-Zeeman observations in four molecular-cloud envelopes and we had shown that, contrary to claims by Crutcher et al., there is no evidence that the mass-to-flux ratio decreases from the envelopes to the cores of these clouds. The key difference between our data analysis and the earlier one by Crutcher et al. is the relaxation of the overly restrictive assumption made by Crutcher et al, that the magnetic field strength is independent of position in each of the four envelopes. In a more recent paper, Crutcher et al. (1) claim that our analysis is not self-consistent, in that it misses a cosine factor, and (2) present new arguments to support their contention that the magnetic-field strength is indeed independent of position in each of the four envelopes. We show that the claim of the missing cosine factor is false, that the new arguments contain even more serious problems than the Crutcher et al. original data analysis, and we present new observational evidence, independent of the OH-Zeeman data, that suggests significant variations in the magnetic-field strength in the four cloud envelopes.

Recent observations of column densities in molecular clouds find lognormal distributions with power-law high-density tails. These results are often interpreted as indications that supersonic turbulence dominates the dynamics of the observed clouds. We calculate and present the column-density distributions of three clouds, modeled with very different techniques, none of which is dominated by supersonic turbulence. The first star-forming cloud is simulated using smoothed particle hydrodynamics (SPH); in this case gravity, opposed only by thermal-pressure forces, drives the evolution. The second cloud is magnetically subcritical with subsonic turbulence, simulated using nonideal MHD; in this case the evolution is due to gravitationally-driven ambipolar diffusion. The third cloud is isothermal, self-gravitating, and has a smooth density distribution analytically approximated with a uniform inner region and an r^-2 profile at larger radii. We show that in all three cases the column-density distributions are lognormal. Power-law tails develop only at late times (or, in the case of the smooth analytic profile, for strongly centrally concentrated configurations), when gravity dominates all opposing forces. It therefore follows that lognormal column-density distributions are generic features of diverse model clouds, and should not be interpreted as being a consequence of supersonic turbulence.

In a previous paper we formulated the problem of the formation and evolution of fragments (or cores) in magnetically-supported, self-gravitating molecular clouds in axisymmetric geometry, accounting for the effects of ambipolar diffusion and Ohmic dissipation, grain chemistry and dynamics, and radiative transfer. Here we present results of star formation simulations that accurately track the evolution of a protostellar fragment over eleven orders of magnitude in density (from 300 cm^-3 to \approx 10^14 cm^-3), i.e., from the early ambipolar-diffusion--initiated fragmentation phase, through the magnetically supercritical, dynamical-contraction phase and the subsequent magnetic decoupling stage, to the formation of a protostellar core in near hydrostatic equilibrium. As found by Fiedler & Mouschovias (1993), gravitationally-driven ambipolar diffusion leads to the formation and subsequent dynamic contraction of a magnetically supercritical core. Moreover, we find that ambipolar diffusion, not Ohmic dissipation, is responsible for decoupling all the species except the electrons from the magnetic field, by a density \approx 3 x 10^12 cm^-3. Magnetic decoupling precedes the formation of a central stellar object and ultimately gives rise to a concentration of magnetic flux (a `magnetic wall') outside the hydrostatic core --- as also found by Tassis & Mouschovias (2005a,b) through a different approach. At approximately the same density at which Ohmic dissipation becomes more important than ambipolar diffusion (\gtrsim 7 x 10^12 cm^-3), the grains carry most of the electric charge as well as the electric current. The prestellar core remains disclike down to radii ~ 10 AU, inside which thermal pressure becomes important. The magnetic flux problem of star formation is resolved for at least strongly magnetic newborn stars by this stage of the evolution, i.e., by a central density \approx 10^14 cm^-3. The hydrostatic core has radius \approx 2 AU, density \approx 10^14 cm^-3, temperature \approx 300 K, magnetic field strength \approx 0.2 G, magnetic flux \approx 5 x 10^18 Wb, luminosity ~ 10^-3 L_\odot, and mass ~ 10^-2 M_\odot.

The ambipolar-diffusion theory of star formation predicts the formation of fragments in molecular clouds with mass-to-flux ratios greater than that of the parent-cloud envelope. By contrast, scenarios of turbulence-induced fragmentation do not yield such a robust prediction. Based on this property, Crutcher et al. (2009) proposed an observational test that could potentially discriminate between fragmentation theories. However, the analysis applied to the data severely restricts the discriminative power of the test: the authors conclude that they can only constrain what they refer to as the "idealized" ambipolar-diffusion theory that assumes initially straight-parallel magnetic field lines in the parent cloud. We present an original, self-consistent analysis of the same data taking into account the nonuniformity of the magnetic field in the cloud envelopes, which is suggested by the data themselves, and we discuss important geometrical effects that must be accounted for in using this test. We show quantitatively that the quality of current data does not allow for a strong conclusion about any fragmentation theory. Given the discriminative potential of the test, we urge for more and better-quality data.

We show that the ambipolar-diffusion--initiated fragmentation of molecular clouds leads simply and naturally to an initial core mass function (CMF) which is very similar to the initial stellar mass function (IMF) and in excellent agreement with existing observations. This agreement is robust provided that the three (input) free parameters remain within their range of values suggested by observations. Other, observationally testable, predictions are made.

We follow numerically the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold HI interstellar medium in two spatial dimensions. The nonlinear evolution of the system favors modes that allow the magnetic field lines to cross the galactic plane. Cold HI clouds form with typical masses ~= 10^5 M_sun, mean densities ~= 20 cm^-3, mean magnetic field strengths ~= 4.3 muG (rms field strengths ~= 6.4 muG), mass-to-flux ratios ~= 0.1 - 0.3 relative to critical, temperatures ~= 50 K, (two-dimensional) turbulent velocity dispersions ~= 1.6 km s^-1, and separations ~= 500 pc, in agreement with observations. The maximum density and magnetic field strength are ~= 10^3 cm^-3 and ~= 20 muG, respectively. Approximately 60% of all HI mass is in the warm neutral medium. The cold neutral medium is arranged into sheet-like structures both perpendicular and parallel to the galactic plane, but it is also found almost everywhere in the galactic plane, with the density being highest in valleys of the magnetic field lines. `Cloudlets' also form whose physical properties are in quantitative agreement with those observed for such objects by Heiles (1967). The nonlinear phase of the evolution takes ~< 30 Myr, so that, if the instability is triggered by a nonlinear perturbation such as a spiral density shock wave, interstellar clouds can form within a time suggested by observations.

We formulate the problem of the formation and subsequent evolution of fragments (or cores) in magnetically-supported, self-gravitating molecular clouds in two spatial dimensions. The six-fluid (neutrals, electrons, molecular and atomic ions, positively-charged, negatively-charged, and neutral grains) physical system is governed by the radiative, nonideal magnetohydrodynamic (RMHD) equations. The magnetic flux is not assumed to be frozen in any of the charged species. Its evolution is determined by a newly-derived generalized Ohm's law, which accounts for the contributions of both elastic and inelastic collisions to ambipolar diffusion and Ohmic dissipation. The species abundances are calculated using an extensive chemical-equilibrium network. Both MRN and uniform grain size distributions are considered. The thermal evolution of the protostellar core and its effect on the dynamics are followed by employing the grey flux-limited diffusion approximation. Realistic temperature-dependent grain opacities are used that account for a variety of grain compositions. We have augmented the publicly-available Zeus-MP code to take into consideration all these effects and have modified several of its algorithms to improve convergence, accuracy and efficiency. Results of magnetic star formation simulations that accurately track the evolution of a protostellar fragment from a density ~10^3 cm^-3 to a density ~10^15 cm^-3, while rigorously accounting for both nonideal MHD processes and radiative transfer, are presented in a separate paper.

Until recently, many of the dozens of quantitative predictions of the ambipolar-diffusion theory of gravitational fragmentation (or core formation) of molecular clouds have been confirmed by observations and, just as importantly, no prediction has been contradicted by any observation. A recent paper, however, claims that measurements of the variation of the mass-to-flux ratio from envelopes to cores in four clouds {\it de}creases, in direct contrast to a prediction of the theory but in agreement with turbulent fragmentation (in the absence of gravity) and, therefore, the ambipolar-diffusion theory is invalid (Crutcher et al 2008). The paper treats magnetic-field nondetections as if they were detections. We show that the analysis of the data is fundamentally flawed and, moreover, the comparison with the theoretical prediction ignores major geometrical effects, suggested by the data themselves if taken at face value. The magnetic fluxes of the envelopes are also miscalculated. We carry out a proper error analysis and treatment of the nondetections and we show that the claimed measurement of the variation of the mass-to-flux ratio from envelopes to cores is not valid, no contradiction with the ambipolar-diffusion theory can be concluded, and no theory can be tested on the basis of these data.

Using interferometric data from BIMA observations, combined with detailed modeling in Fourier space of the physical structures predicted by models, we constrain the circumstellar envelope parameters for four Class 0 young stellar objects, as well as their embedded circumstellar disks. The envelopes of these objects are still undergoing collapse, and theoretical collapse models can be compared to the observations. Since it has been suggested in a previous study that both the Larson-Penston and Shu similarity solutions underestimate the age of the system, we adopt Tassis & Mouschovias' model of the collapse process, which includes all relevant magnetic fields effects. The results of the model fitting show a good consistency between theory and data; furthermore, no age problem exists since the Tassis & Mouschovias' model is age independent for the first 255 kyr. Although the majority of the continuum dust emission arises from the circumstellar envelopes, these objects have well known outflows, which suggest the presence of circumstellar disks. At the highest resolution, most of the large-scale envelope emission is resolved out by interferometry, but the small-scale residual emission remains, making it difficult to observe only the compact disk component. By modeling the emission of the envelope and subtracting it from the total emission, we constrain the disk masses in our four systems to be comparable to or smaller than the typical disk masses for T Tauri systems.

We follow the ambipolar-diffusion--driven formation and evolution of a fragment in a magnetically supported molecular cloud, until a hydrostatic protostellar core forms at its center. This problem was formulated in Paper I. We determine the density, velocity and magnetic field as functions of space and time, and the contribution of ambipolar diffusion and Ohmic dissipation to the resolution of the magnetic flux problem of star formation. The issue of whether the magnetic field ever decouples from the (neutral) matter is also addressed. We also find that the electrons do not decouple from the field lines before thermal ionization becomes important and recouples the magnetic field to the neutral matter. Ohmic dissipation becomes more effective than ambipolar diffusion as a flux reduction mechanism only at the highest densities (a few times 10^12 particles per cubic cm). In the high-density central parts of the core, the magnetic field acquires an almost spatially uniform structure, with a value that, at the end of the calculation (number density ~ 5 times 10^14 particles per cubic cm), is found to be in excellent agreement with meteoritic measurements of magnetic fields in the protosolar nebula. Outside the hydrostatic protostellar core, a concentration of magnetic flux (a "magnetic wall") forms, which gives rise to a magnetic shock. This magnetic shock is the precursor of the repeated shocks found by Tassis & Mouschovias (2005) which cause spasmodic accretion onto the hydrostatic core at later times.

We formulate the problem of the formation of magnetically supercritical cores in magnetically subcritical parent molecular clouds, and the subsequent collapse of the cores to high densities, past the detachment of ions from magnetic field lines and into the opaque regime. We employ the six-fluid MHD equations, accounting for the effects of grains (negative, positive and neutral) including their inelastic collisions with other species. We do not assume that the magnetic flux is frozen in any of the charged species. We derive a generalized Ohm's law that explicitly distinguishes between flux advection (and the associated process of ambipolar diffusion) and Ohmic dissipation, in order to assess the contribution of each mechanism to the increase of the mass-to-flux ratio of the central parts of a collapsing core and possibly to the resolution of the magnetic flux problem of star formation. We show how our formulation is related to and can be transformed into the traditional, directional formulation of the generalized Ohm's law, and we derive formulae for the perpendicular, parallel and Hall conductivities entering the latter, which include, for the first time, the effect of inelastic collisions between grains. In addition, we present a general (valid in any geometry) solution for the velocities of charged species as functions of the velocity of the neutrals and of the effective flux velocity (which can in turn be calculated from the dynamics of the system and Faraday's law). The last two sets of formulae can be adapted for use in any general non-ideal MHD code to study phenomena beyond star formation in magnetic clouds. The results, including a detailed parameter study, are presented in two accompanying papers.

In two previous papers we formulated and solved, for a fiducial set of free parameters, the problem of the formation and evolution of a magnetically supercritical core inside a magnetically subcritical parent cloud. In this paper we present a parameter study to assess the sensitivity of the results (1) to the density at which the equation of state becomes adiabatic; (2) to the initial mass-to-flux ratio of the parent cloud; and (3) to ionization by radioactive decay of different nuclei (40K and 26Al) at high densities (number density > 10^12 particles per cubic cm). We find that (1) the results depend only slightly on the density at which the onset of adiabaticity occurs; (2) memory of the initial mass-to-flux ratio is completely lost at late times, which emphasizes the relevance of this work, idependently of the adopted theory of core formation; and (3) the precise source of radioactive ionization alters the degree of attachment of the electrons to the field lines (at high densities), and the relative importance of ambipolar diffusion and Ohmic dissipation in reducing the magnetic flux of the protostar. The value of the magnetic field at the end of the runs is insensitive to the values of the free parameters and in excellent agreement with meteoritic measurements of the protosolar nebula magnetic field. The magnetic flux problem of star formation is resolved for at least strongly magnetic newborn stars. A complete detachment of the magnetic field from the matter is unlikely. The formation of a "magnetic wall" (with an associated magnetic shock) is independent of the assumed equation of state, although the process is enhanced and accelerated by the formation of a central hydrostatic core.

We revisit the problem of the star formation timescale and the ages of molecular clouds. The apparent overabundance of star-forming molecular clouds over clouds without active star formation has been thought to indicate that molecular clouds are "short-lived" and that star formation is "rapid". We show that this statistical argument lacks self-consistency and, even within the rapid star-formation scenario, implies cloud lifetimes of approximately 10 Myr. We discuss additional observational evidence from external galaxies that indicate lifetimes of molecular clouds and a timescale of star formation of approximately 10 Myr . These long cloud lifetimes in conjunction with the rapid (approximately 1 Myr) decay of supersonic turbulence present severe difficulties for the scenario of turbulence-controlled star formation. By contrast, we show that all 31 existing observations of objects for which the linewidth, the size, and the magnetic field strength have been reliably measured are in excellent quantitative agreement with the predictions of the ambipolar-diffusion theory. Within the ambipolar-diffusion-controlled star formation theory the linewidths may be attributed to large-scale non-radial cloud oscillations (essentially standing large-amplitude, long-wavelength Alfven waves), and the predicted relation between the linewidth, the size, and the magnetic field is a natural consequence of magnetic support of self-gravitating clouds.

We formulate the problem of the late accretion phase of the evolution of an isothermal magnetic disk surrounding a forming star. The evolution is described by the six-fluid MHD equations, accounting for the presence of neutrals, atomic and molecular ions, electrons, and neutral, positively, and negatively charged grains. Only the electron fluid is assumed to be attached to the magnetic field, in order to investigate the effect of the detachment of the ions from the magnetic field lines that begins at densities as low as 10^8 cm^-3. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 10^11 cm^-3. In this way, the structure and evolution of the isothermal disk surrounding the forming star can be studied at late times without having to implement the numerically costly radiative transfer required by the physics of the opaque core. The mass and magnetic flux accumulating in the forming star arecalculated, as are their effects on the structure & evolution of the surrounding disk. The numerical method of solution first uses an adaptive grid and later, after a central region a few AU in radius becomes opaque, switches to a stationary but nonuniform grid with a central sink cell. It also involves an implicit time integrator, an advective difference scheme that possesses the transportive property, a second-order difference approximation of forces inside a cell, an integral approximation of the gravitational and magnetic fields, and tensor artificial viscosity that permits an accurate investigation of the formation and evolution of shocks in the neutral fluid.

The problem of the late accretion phase of the evolution of an axisymmetric, isothermal magnetic disk surrounding a forming star has been formulated in a companion paper. The "central sink approximation" is used to circumvent the problem of describing the evolution inside the opaque central region for densities greater than 10^11 cm^-3 and radii smaller than a few AUs. Only the electrons are assumed to be attached to the magnetic field lines, and the effects of both negatively and positively charged grains are accounted for. After a mass of 0.1 solar mass accumulates in the central cell (forming star), a series of magnetically driven outflows and associated outward propagating shocks form in a quasi-periodic fashion. As a result, mass accretion onto the protostar occurs in magnetically controlled bursts. We refer to this process as spasmodic accretion. The shocks propagate outward with supermagnetosonic speeds. The period of dissipation and revival of the outflow decreases in time, as the mass accumulated in the central sink increases. We evaluate the contribution of ambipolar diffusion to the resolution of the magnetic flux problem of star formation during the accretion phase, and we find it to be very significant although not sufficient to resolve the entire problem yet. Ohmic dissipation is completely negligible in the disk during this phase of the evolution. The protostellar disk is found to be stable against interchange-like instabilities, despite the fact that the mass-to-flux ratio has temporary local maxima.

We re-examine critically the estimates of the duration of different phases of star formation and the lifetimes of molecular clouds, based on the ages of T-Tauri stars, age spreads of stars in clusters, and statistics of pre-stellar cores. We show that all available observational data are consistent with lifetimes of molecular clouds comparable to 10 Myr, as well as with the predictions of the theory of self-initiated, ambipolar-diffusion--controlled star formation. We conclude that there exists no observational support for either "young" molecular clouds or "rapid" star formation.

This is the second in a series of papers on the effects of dust on multifluid, MHD shock waves in weakly ionized molecular gas. We investigate the influence of dust on the critical shock speed, v_crit, above which C shocks cease to exist. Chernoff showed that v_crit cannot exceed the grain magnetosound speed, v_gms, if dust grains are dynamically well coupled to the magnetic field. We present numerical simulations of steady shocks where the grains may be well- or poorly coupled to the field. We use a time-dependent, multifluid MHD code that models the plasma as a system of interacting fluids: neutral particles, ions, electrons, and various ``dust fluids'' comprised of grains with different sizes and charges. Our simulations include grain inertia and grain charge fluctuations but to highlight the essential physics we assume adiabatic flow, single-size grains, and neglect the effects of chemistry. We show that the existence of a phase speed v_phi does not necessarily mean that C shocks will form for all shock speeds v_s less than v_phi. When the grains are weakly coupled to the field, steady, adiabatic shocks resemble shocks with no dust: the transition to J type flow occurs at v_crit = 2.76 v_nA, where v_nA is the neutral Alfven speed, and steady shocks with v_s > 2.76 v_nA are J shocks with magnetic precursors in the ion-electron fluid. When the grains are strongly coupled to the field, v_crit = min(2.76 v_nA, v_gms). Shocks with v_crit < v_s < v_gms have magnetic precursors in the ion-electron-dust fluid. Shocks with v_s > v_gms have no magnetic precursor in any fluid. We present time-dependent calculations to study the formation of steady multifluid shocks. The dynamics differ qualitatively depending on whether or not the grains and field are well coupled.