We reexamine the unipolar induction mechanism for heating asteroids originally proposed in a classic series of papers by Sonett and collaborators. As originally conceived, induction heating is caused by the "motional electric field" which appears in the frame of an asteroid immersed in a fully-ionized, magnetized solar wind and drives currents through its interior. However we point out that classical induction heating contains a subtle conceptual error, in consequence of which the electric field inside the asteroid was calculated incorrectly. The problem is that the motional electric field used by Sonett et al. is the electric field in the freely streaming plasma far from the asteroid; in fact the motional field vanishes at the asteroid surface for realistic assumptions about the plasma density. In this paper we revisit and improve the induction heating scenario by: (1) correcting the conceptual error by self consistently calculating the electric field in and around the boundary layer at the asteroid-plasma interface; (2) considering weakly-ionized plasmas consistent with current ideas about protoplanetary disks; and (3) considering more realistic scenarios which do not require a fully ionized, powerful T Tauri wind in the disk midplane. We present exemplary solutions for two highly idealized flows which show that the interior electric field can either vanish or be comparable to the fields predicted by classical induction depending on the flow geometry. We term the heating driven by these flows "electrodynamic heating", calculate its upper limits, and compare them to heating produced by short-lived radionuclides.

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

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.

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 present extensive calculations on the efficiency of grain alignment by the Davis-Greenstein mechanism. We model the grains as oblate spheroids with arbitrary axis ratios. Our description of the grain dynamics includes (i) magnetic dissipation and the inverse process driven by thermal fluctuations in the grain magnetization; (ii) gas-grain collisions and thermal evaporation of molecules from the grain surface; (iii) the transformation of rotational energy into heat by the Barnett effect and the inverse process driven by thermal fluctuations; and (iv) rapid Larmor precession of the grain angular momentum about the interstellar magnetic field. For ordinary paramagnetic grains, we calculate the Rayleigh reduction factor, R, for more than 1000 combinations of the 3 dimensionless parameters which characterise the alignment. For superparamagnetic grains, we calculate R from an exact analytic solution for the relevant distribution function. We calibrate the accuracy of a recently-proposed perturbative approximation and show that it yields R values with a mean error of about 17 percent.

The alignment of interstellar dust grains is described by the joint distribution function for certain ``internal'' and ``external'' variables, where the former describe the orientation of a grain's axes with respect to its angular momentum, J, and the latter describe the orientation of J relative to the interstellar magnetic field. I show how the large disparity between the dynamical timescales of the internal and external variables--- which is typically 2--3 orders of magnitude--- can be exploited to greatly simplify calculations of the required distribution. The method is based on an ``adiabatic approximation'' which closely resembles the Born-Oppenheimer approximation in quantum mechanics. The adiabatic approximation prescribes an analytic distribution function for the ``fast'' dynamical variables and a simplified Fokker-Planck equation for the ``slow'' variables which can be solved straightforwardly using various techniques. These solutions are accurate to order epsilon, where epsilon is the ratio of the fast and slow dynamical timescales. As a simple illustration of the method, I derive an analytic solution for the joint distribution established when Barnett relaxation acts in concert with gas damping. The statistics of the analytic solution agree with the results of laborious numerical calculations which do not exploit the adiabatic approximation.

We use polarimetric observations of two stars (HD29647, HD283809) in the general direction of TMC-1 in the Taurus Dark Cloud to investigate grain properties and cloud structure in this region. We show the data to be consistent with a simple two-component model, in which general interstellar polarization in the Taurus Cloud is produced by a widely distributed cloud component with relatively uniform magnetic field orientation; the light from stars close to TMC-1 suffers additional polarization arising in one (or more) subcloud(s) with larger average grain size and different magnetic field directions compared with the general trend. Towards HD29647, in particular, we show that the unusually low degree of visual polarization relative to extinction is due to the presence of distinct cloud components in the line of sight with markedly different magnetic field orientations. Stokes parameter calculations allow us to separate out the polarization characteristics of the individual components. Results are fit with the Serkowski empirical formula to determine the degree and wavelength of maximum polarization. Whereas lambda_max values in the widely distributed material are similar to the average (0.55um) for the diffuse interstellar medium, the subcloud in line of sight to HD~283809, the most heavily reddened star in our study, has lambda_max approx. 0.73um, indicating the presence of grains about 30% larger than this average. Our model also predicts detectable levels of circular polarization toward both HD~29647 and HD~283809.

The recent detection of interstellar polarization in the solid CO feature near 4.67 micron shows that CO-mantled grains can be aligned in cold molecular clouds. These observations conflict with a theory of grain alignment which attributes the polarization in molecular clouds to the effects of cosmic rays: according to this theory, oblate spheroidal grains with H_2O and CO_2-dominated ice mantles are spun up to suprathermal energies by molecular evaporation from cosmic ray impact sites but spin up does not occur for CO-mantled grains. Motivated by this conflict, we reexamine the effects of cosmic rays on the alignment of icy grains. We show that the systematic torques produced by cosmic rays are insufficient to cause suprathermal spin. In principle, the random torques due to cosmic rays can enhance the efficiency of Davis-Greenstein alignment by raising the grain rotational temperature. However, a significant enhancement would require cosmic ray fluxes 6--7 orders of magnitude larger than the flux in a typical cold cloud.

We present an exact formulation of the physics of Barnett relaxation. Our formulation is based on a realistic kinetic model of the relaxation mechanism which includes the alignment of the grain angular momentum in body coordinates by Barnett dissipation, disalignment by thermal fluctuations, and coupling of the angular momentum to the gas via gas damping. We solve the Fokker-Planck equation for the measure of internal alignment using numerical integration of the equivalent Langevin equation for Brownian rotation. The accuracy of our results is calibrated by comparing our numerical solutions with exact analytic results obtained for special cases.We describe an analytic approximation for the measure of alignment which fits our numerical results for cases of practical interest.

We describe an implicit prediction of the accretion disk models constructed by Wardle and Konigl (1990) for the circumnuclear disk (CND) of gas and dust near the Galactic center: supersonic ambipolar diffusion, an essential dynamical ingredient of the Wardle-Konigl disks, will cause the alignment of dust grains due to a process described by Roberge, Hanany, & Messinger (1995). We calculate synthetic maps of the polarized thermal emission which would be caused by ambipolar alignment in the preferred Wardle-Konigl model. Our maps are in reasonable agreement with 100 micron polarimetry of the CND if we assume that the grains have shapes similar to those of grains in nearby molecular clouds and that the CND contains a disordered magnetic field in energy equipartition with its ordered field.