We describe a numerical scheme for studying time-dependent, multifluid, magnetohydrodynamic shock waves in weakly ionized interstellar clouds and cores. Shocks are modeled as propagating perpendicular to the magnetic field and consist of a neutral molecular fluid plus a fluid of ions and electrons. The scheme is based on operator splitting, wherein time integration of the governing equations is split into separate parts. In one part independent homogeneous Riemann problems for the two fluids are solved using Godunov's method. In the other equations containing the source terms for transfer of mass, momentum, and energy between the fluids are integrated using standard numerical techniques. We show that, for the frequent case where the thermal pressures of the ions and electrons are << magnetic pressure, the Riemann problems for the neutral and ion-electron fluids have a similar mathematical structure which facilitates numerical coding. Implementation of the scheme is discussed and several benchmark tests confirming its accuracy are presented, including (i) MHD wave packets ranging over orders of magnitude in length and time scales; (ii) early evolution of mulitfluid shocks caused by two colliding clouds; and (iii) a multifluid shock with mass transfer between the fluids by cosmic-ray ionization and ion-electron recombination, demonstrating the effect of ion mass loading on magnetic precursors of MHD shocks. An exact solution to a MHD Riemann problem forming the basis for an approximate numerical solver used in the homogeneous part of our scheme is presented, along with derivations of the analytic benchmark solutions and tests showing the convergence of the numerical algorithm.

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

We study the physics of wave propagation in a weakly ionised plasma, as it applies to the formation of multifluid, MHD shock waves. We model the plasma as separate charged and neutral fluids which are coupled by ion-neutral friction. At times much less than the ion-neutral drag time, the fluids are decoupled and so evolve independently. At later times, the evolution is determined by the large inertial mismatch between the charged and neutral particles. The neutral flow continues to evolve independently; the charged flow is driven by and slaved to the neutral flow by friction. We calculate this driven flow analytically by considering the special but realistic case where the charged fluid obeys linearized equations of motion. We carry out an extensive analysis of linear, driven, MHD waves. The physics of driven MHD waves is embodied in certain Green functions which describe wave propagation on short time scales, ambipolar diffusion on long time scales, and transitional behavior at intermediate times. By way of illustration, we give an approximate solution for the formation of a multifluid shock during the collision of two identical interstellar clouds. The collision produces forward- and reverse J shocks in the neutral fluid and a transient in the charged fluid. The latter rapidly evolves into a pair of magnetic precursors on the J shocks, wherein the ions undergo force free motion and the magnetic field grows monotonically with time. The flow appears to be self similar at the time when linear analysis ceases to be valid.

We formulate the problem of the formation and collapse of nonaxisymmetric protostellar cores in weakly ionized, self-gravitating, magnetic molecular clouds. In our formulation, molecular clouds are approximated as isothermal, thin (but with finite thickness) sheets. We present the governing dynamical equations for the multifluid system of neutral gas and ions, including ambipolar diffusion, and also a self-consistent treatment of thermal pressure, gravitational, and magnetic (pressure and tension) forces. The dimensionless free parameters characterizing model clouds are discussed. The response of cloud models to linear perturbations is also examined, with particular emphasis on length and time scales for the growth of gravitational instability in magnetically subcritical and supercritical clouds. We investigate their dependence on a cloud's initial mass-to-magnetic-flux ratio (normalized to the critical value for collapse), the dimensionless initial neutral-ion collision time, and also the relative external pressure exerted on a model cloud. Among our results, we find that nearly-critical model clouds have significantly larger characteristic instability lengthscales than do more distinctly sub- or supercritical models. Another result is that the effect of a greater external pressure is to reduce the critical lengthscale for instability. Numerical simulations showing the evolution of model clouds during the linear regime of evolution are also presented, and compared to the results of the dispersion analysis. They are found to be in agreement with the dispersion results, and confirm the dependence of the characteristic length and time scales on parameters such as the initial mass-to-flux ratio and relative external pressure.

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.

We extend our earlier work on ambipolar diffusion induced formation of protostellar cores in isothermal sheet-like magnetic interstellar clouds, by studying nonaxisymmetric collapse for the physically interesting regime of magnetically critical and supercritical model clouds ($\mui \geq 1$, where $\mui$ is the initial mass-to-magnetic flux ratio in units of the critical value for gravitational collapse). Cores that form in model simulations are effectively triaxial, with shapes that are typically closer to being oblate, rather than prolate. Infall velocities in the critical model ($\mui = 1$) are subsonic; in contrast, a supercritical model ($\mui = 2$) has extended supersonic infall that may be excluded by observations. For the magnetically critical model, ambipolar diffusion forms cores that are supercritical ($\muc > 1$) and embedded within subcritical envelopes ($\muenv < 1$). Cores in our models have density profiles that eventually merge into a near-uniform background, which is suggestive of observed properties of cloud cores.

This is the first in a series of papers on the effects of dust on the formation, propagation, and structure of nonlinear MHD waves and MHD shocks in weakly-ionized plasmas. We model the plasma as a system of 9 interacting fluids, consisting of the neutral gas, ions, electrons, and 6 grain fluids comprised of very small grains or PAHs and classical grains in different charge states. We formulate the governing equations for perpendicular shocks under approximations appropriate for dense molecular clouds, protostellar cores, and protoplanetary disks. We describe a code that obtains numerical solutions using a finite difference method, and estabish its accuracy by comparing numerical and exact solutions for special cases.

We review theoretical models of the early stages of star formation, in which gravitational collapse is strongly regulated by magnetic fields and the associated process of ambipolar diffusion. We discuss results of numerical simulations and analytical studies of core formation and collapse, which can be directly tested against observation. We also focus on recent data which are relevant to this theory of star formation, such as: observations of extended infall in protostellar cores, estimates of evolutionary timescales ~ 1 Myr for cores, measured mass-to-flux ratios of cores, and the relative alignment of polarization vectors with apparent cloud elongation. It is shown that in all of these areas, the data lie within the observationally allowable range of parameter space. Other areas of interest (protostellar accretion rates and the presence of core edges) and issues that remain unresolved or under study (the role of non-thermal motions and cluster formation) are also discussed. Moreover, we highlight some differences between our model predictions and those of highly turbulent star formation models, and discuss how these differences can be distinguished observationally.

We revisit the problem of the formation of dense protostellar cores due to ambipolar diffusion within magnetically supported molecular clouds, and derive an analytical expression for the core formation timescale. The resulting expression is similar to the canonical expression = t_{ff}^2/t_{ni} ~ 10 t_{ni} (where t_{ff} is the free-fall time and t_{ni} is the neutral-ion collision time), except that it is multiplied by a numerical factor C(\mu_{c0}), where \mu_{c0} is the initial central mass-to-flux ratio normalized to the critical value for gravitational collapse. C(\mu_{c0}) is typically ~ 1 in highly subcritical clouds (\mu_{c0} << 1), although certain conditions allow C(\mu_{c0}) >> 1. For clouds that are not highly subcritical, C(\mu_{c0}) can be much less than unity, with C(\mu_{c0}) --> 0 for \mu_{c0} --> 1, significantly reducing the time required to form a supercritical core. This, along with recent observations of clouds with mass-to-flux ratios close to the critical value, may reconcile the results of ambipolar diffusion models with statistical analyses of cores and YSO's which suggest an evolutionary timescale \~ 1 Myr for objects of mean density ~ 10^4 cm^{-3}. We compare our analytical relation to the results of numerical simulations, and also discuss the effects of dust grains on the core formation timescale.

Recent high-resolution studies of the L1544 protostellar core by Tafalla et al. and Williams et al. reveal the structure and kinematics of the gas. The observations of this prestellar core provide a natural test for theoretical models of core formation and evolution. Based on their results, the above authors claim a discrepancy with the implied infall motions from ambipolar diffusion models. In this paper, we reexamine the earlier ambipolar diffusion models, and conclude that the L1544 core can be understood to be a magnetically supercritical core undergoing magnetically diluted collapse. We also present a new model specifically designed to simulate the formation and evolution of the L1544 core. This model, which uses reasonable input parameters, yields mass and radial density distributions, as well as neutral and ion infall speed profiles, that are in very good agreement with physical values deduced by observations. The lifetime of the core is also in good agreement with prestellar core lifetimes estimated from statistics of an ensemble of cores. The observational input can act to constrain other currently unobserved quantities such as the degree of ionization, and the background magnetic field strength and orientation near the L1544 core.

We obtain self-similar solutions that describe the gravitational collapse of nonrotating, isothermal, magnetic molecular cloud cores. We use simplifying assumptions but explicitly include the induction equation, and the semianalytic solutions we derive are the first to account for the effects of ambipolar diffusion following the formation of a central point mass. Our results demonstrate that, after the protostar first forms, ambipolar diffusion causes the magnetic flux to decouple in a growing region around the center. The decoupled field lines remain approximately stationary and drive a hydromagnetic C-shock that moves outward at a fraction of the speed of sound (typically a few tenths of a kilometer per second), reaching a distance of a few thousand AU at the end of the main accretion phase for a solar-mass star. We also show that, in the absence of field diffusivity, a contracting core will not give rise to a shock if, as is likely to be the case, the inflow speed near the origin is nonzero at the time of point-mass formation. Although the evolution of realistic molecular cloud cores will not be exactly self-similar, our results reproduce the main qualitative features found in detailed core-collapse simulations (Ciolek & Konigl 1998)

We present a numerical simulation of the dynamical collapse of a nonrotating magnetic molecular cloud core and follow the core's evolution through the formation of a central point mass and its subsequent growth to a 1 solar-mass protostar. The epoch of point-mass formation (PMF) is investigated by a self- consistent extension of previously presented models of core formation and contraction in axisymmetric, self-gravitating, isothermal, magnetically supported interstellar molecular clouds. Prior to PMF, the core is dynamically contracting and is not well approximated by a quasistatic equilibrium model. Ambipolar diffusion, which plays a key role in the early evolution of the core, is unimportant during the dynamical pre-PMF collapse phase. However, the appearance of a central mass, through its effect on the gravitational field in the inner core regions, leads to a "revitalization" of ambipolar diffusion in the weakly ionized gas surrounding the central protostar. This process is so efficient that it leads to a decoupling of the field from the matter and results in an outward propagating hydromagnetic C-type shock. The existence of an ambipolar diffusion-mediated shock was predicted by Li & McKee (1996), and we find that the basic shock structure given by their analytical model is well reproduced by our more accurate numerical results. Our calculation also demonstrates that ambipolar diffusion, rather than Ohmic diffusivity operating in the innermost core region, is the main field decoupling mechanism responsible for driving the shock after PMF.