Helical Fields and Filamentary Molecular Clouds

We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first derive a new virial equation appropriate for magnetized filamentary clouds, which includes the effects of non-thermal motions and the turbulent pressure of the surrounding ISM. When compared with the data, we find that many filamentary clouds have a mass per unit length that is significantly reduced by the effects of external pressure, and that toroidal fields play a significant role in squeezing such clouds. We also develop exact numerical MHD models of filamentary molecular clouds with more general helical field configurations than have previously been considered. We also examine the effects of the equation of state by comparing ``isothermal'' filaments, with constant total (thermal plus turbulent) velocity dispersion, with equilibria constructed using a logatropic equation of state. We perform a Monte Carlo exploration of our parameter space to determine which choices of parameters result in models that agree with the available observational constraints. We find that both equations of state result in equilibria that agree with the observational results. Moreover, we find that models with helical fields have more realistic density profiles than either unmagnetized models or those with purely poloidal fields; we find that most isothermal models have density distributions that fall off as r^{-1.8} to r^{-2}, while logatropes have density profiles that range from r^{-1} to r^{-1.8}. We find that purely poloidal fields produce filaments with steep density gradients that not allowed by the observations.Comment: 21 pages, 8 eps figures, submitted to MNRAS. Significant streamlining of tex

Helical Fields and Filamentary Molecular Clouds II - Axisymmetric Stability and Fragmentation

In Paper I (Fiege & Pudritz, 1999), we constructed models of filamentary molecular clouds that are truncated by a realistic external pressure and contain a rather general helical magnetic field. We address the stability of our models to gravitational fragmentation and axisymmetric MHD-driven instabilities. By calculating the dominant modes of axisymmetric instability, we determine the dominant length scales and growth rates for fragmentation. We find that the role of pressure truncation is to decrease the growth rate of gravitational instabilities by decreasing the self-gravitating mass per unit length. Purely poloidal and toroidal fields also help to stabilize filamentary clouds against fragmentation. The overall effect of helical fields is to stabilize gravity-driven modes, so that the growth rates are significantly reduced below what is expected for unmagnetized clouds. However, MHD ``sausage'' instabilities are triggered in models whose toroidal flux to mass ratio exceeds the poloidal flux to mass ratio by more than a factor of $\sim 2$. We find that observed filaments appear to lie in a physical regime where the growth rates of both gravitational fragmentation and axisymmetric MHD-driven modes are at a minimum.Comment: 16 pages with 18 eps figures. Submitted to MNRA

Protostellar Evolution during Time Dependent, Anisotropic Collapse

The formation and collapse of a protostar involves the simultaneous infall and outflow of material in the presence of magnetic fields, self-gravity, and rotation. We use self-similar techniques to self-consistently model the anisotropic collapse and outflow by a set of angle-separated self-similar equations. The outflow is quite strong in our model, with the velocity increasing in proportion to radius, and material formally escaping to infinity in the finite time required for the central singularity to develop. Analytically tractable collapse models have been limited mainly to spherically symmetric collapse, with neither magnetic field nor rotation. Other analyses usually employ extensive numerical simulations, or either perturbative or quasistatic techniques. Our model is unique as an exact solution to the non-stationary equations of self-gravitating MHD, which features co-existing regions of infall and outflow. The velocity and magnetic topology of our model is quadrupolar, although dipolar solutions may also exist. We provide a qualitative model for the origin and subsequent evolution of such a state. However, a central singularity forms at late times, and we expect the late time behaviour to be dominated by the singularity rather than to depend on the details of its initial state. Our solution may, therefore, have the character of an attractor among a much more general class of self-similarity.Comment: 11 pages, 5 figures, To appear in MNRAS, Memorial paper for M. Aburiha

Anomalous velocity distributions in active Brownian suspensions

Large scale simulations and analytical theory have been combined to obtain the non-equilibrium velocity distribution, $f(v)$, of randomly accelerated particles in suspension. The simulations are based on an event-driven algorithm, generalised to include friction. They reveal strongly anomalous but largely universal distributions which are independent of volume fraction and collision processes, which suggests a one-particle model should capture all the essential features. We have formulated this one-particle model and solved it analytically in the limit of strong damping, where we find that $f(v)$ decays as $1/v$ for multiple decades, eventually crossing over to a Gaussian decay for the largest velocities. Many particle simulations and numerical solution of the one-particle model agree for all values of the damping.Comment: 6 pages, 5 figure