Helical fields and filamentary molecular clouds – I

We study the equilibrium of pressure truncated, filamentary molecular clouds that are threaded by rather general helical magnetic fields. We first apply the virial theorem to filamentary molecular clouds, including 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 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.
Our theoretical models involve three parameters: two to describe the mass loading of the toroidal and poloidal fields, and a third that describes the radial concentration of the filament. We thoroughly explore 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 ⁻², while logatropes have density profiles that range from r ⁻¹ to r ^−1.8. We find that purely poloidal fields produce filaments with steep radial density gradients that are not allowed by the observations.     

Collapse and fragmentation of rotating magnetized clouds – I. Magnetic flux–spin relation

We discuss the evolution of the magnetic flux density and angular velocity in a molecular cloud core, on the basis of three-dimensional numerical simulations, in which a rotating magnetized cloud fragments and collapses to form a very dense optically thick core of  >5 × 10¹⁰ cm⁻³  . As the density increases towards the formation of the optically thick core, the magnetic flux density and angular velocity converge towards a single relationship between the two quantities. If the core is magnetically dominated its magnetic flux density approaches  1.5( n /5 × 10¹⁰ cm⁻³)^1/2 mG  , while if the core is rotationally dominated the angular velocity approaches  2.57 × 10⁻³ ( n /5 × 10¹⁰ cm⁻³)^1/2 yr⁻¹  , where n is the density of the gas. We also find that the ratio of the angular velocity to the magnetic flux density remains nearly constant until the density exceeds  5 × 10¹⁰ cm⁻³  . Fragmentation of the very dense core and emergence of outflows from fragments will be shown in the subsequent paper.     

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 using 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 is 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 magnetohydrodynamics (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 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.     

Formation of interstellar clouds: Parker instability with phase transitions

We numerically follow the nonlinear evolution of the Parker instability in the presence of phase transitions from a warm to a cold H  i interstellar medium in two spatial dimensions. The nonlinear evolution of the system favours modes that allow the magnetic field lines to cross the galactic plane. Cold H  i clouds form with typical masses  ≃10⁵ M_⊙  , mean densities  ≃20 cm⁻³  , mean magnetic-field strengths  ≃4.3 μG  (rms field strengths  ≃6.4 μG  ), mass-to-flux ratios  ≃0.1–0.3  relative to critical, temperatures  ≃50 K  , (two-dimensional) turbulent velocity dispersions  ≃1.6 km s⁻¹  and separations  ≃500 pc  , in agreement with observations. The maximum density and magnetic-field strength are  ≃10³ cm⁻³  and  ≃20 μG  , respectively. Approximately 60 per cent of all H  i 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. 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.     

Helicoidal magneto-electron waves in interstellar molecular clouds

We investigate the evolution of the magnetic flux density in a magnetically supported molecular cloud driven by Hall and Ohmic components of the electric field generated by the flows of thermal electrons. Particular attention is given to the wave transport of the magnetic field in a cloud whose gas dynamics is dominated by electron flows; the mobility of neutrals and ions is regarded as heavily suppressed. It is shown that electromagnetic waves penetrating such a cloud can be converted into helicons – weakly damped, circularly polarized waves in which the densities of the magnetic flux and the electron current undergo coherent oscillations. These waves are interesting in their own right, because for electron magnetohydrodynamics the low-frequency helicoidal waves have the same physical significance as the transverse Alfvén waves do for a single-component magnetohydrodynamics. The latter, as is known, are considered to be responsible for the widths of molecular lines detected in dark, magnetically supported clouds. From our numerical estimates for the group velocity and the rate of dissipation of helicons it follows that a possible contribution of these waves to the broadening of molecular lines is consistent with the conditions typical of dark molecular clouds.     

A three-dimensional numerical method for modelling weakly ionized plasmas

Astrophysical fluids under the influence of magnetic fields are often subjected to single- or two-fluid approximations. In the case of weakly ionized plasmas, however, this can be inappropriate due to distinct responses from the multiple constituent species to both collisional and non-collisional forces. As a result, in dense molecular clouds and protostellar accretion discs, for instance, the conductivity of the plasma may be highly anisotropic leading to phenomena such as Hall and ambipolar diffusion strongly influencing the dynamics.
Diffusive processes are known to restrict the stability of conventional numerical schemes which are not implicit in nature. Furthermore, recent work establishes that a large Hall term can impose an additional severe stability limit on standard explicit schemes. Following a previous paper, which presented the one-dimensional case, we describe a fully three-dimensional method which relaxes the normal restrictions on explicit schemes for multifluid processes. This is achieved by applying the little-known Super TimeStepping technique to the symmetric (ambipolar) component of the evolution operator for the magnetic field in the local plasma rest frame, and the new Hall Diffusion Scheme to the skew-symmetric (Hall) component.     

Three-dimensional simulations of molecular cloud fragmentation regulated by magnetic fields and ambipolar diffusion

We employ the first fully three-dimensional simulation to study the role of magnetic fields and ion–neutral friction in regulating gravitationally driven fragmentation of molecular clouds. The cores in an initially subcritical cloud develop gradually over an ambipolar diffusion time while the cores in an initially supercritical cloud develop in a dynamical time. The infalling speeds on to cores are subsonic in the case of an initially subcritical cloud, while an extended (≳0.1 pc) region of supersonic infall exists in the case of an initially supercritical cloud. These results are consistent with previous two-dimensional simulations. We also found that a snapshot of the relation between density (ρ) and the strength of the magnetic field ( B ) at different spatial points of the cloud coincides with the evolutionary track of an individual core. When the density becomes large, both the relations tend to   B ∝ρ^0.5  .     

On the influence of cooling and heating processes on the Parker instability – II. Numerical simulations

We investigate the Parker instability (PI) influenced by thermal processes in a non-adiabatic, gravitationally stratified interstellar medium and discuss a model including the photoionization heating together with the supplemental heating mechanisms postulated by Reynolds, Haffner and Tufte. A cooling rate due to radiative losses is described by an approximation to the realistic cooling function of Dalgarno and McCray for ionized interstellar gas. An unperturbed initial state of the system simultaneously represents both a magnetohydrostatic and thermal equilibrium, and is thermally stable. We perform a set of 3D numerical magnetohydrodynamic simulations using the zeusmp code. We find that PI developing in the presence of non-adiabatic effects promotes a transition of gas in magnetic valleys to a thermally unstable regime. We find that the region of initially enhanced density due to PI starts to condense more as the result of thermal instability action. The density in this region rises above the classical isothermal limit of two times the equilibrium value at the mid-plane. The maximum density in an evolved system reaches 10–40 times the equilibrium value at the mid-plane, and the structures so formed attain oval shapes. These results lead to the conclusion that PI, operating in the presence of realistic cooling and heating processes, can trigger the formation of dense clouds, which may give rise to giant molecular complexes.     

Magnetic fields and the dynamics of spiral galaxies

We investigate the dynamics of magnetic fields in spiral galaxies by performing 3D magnetohydrodynamics simulations of galactic discs subject to a spiral potential using cold gas, warm gas and a two-phase mixture of both. Recent hydrodynamic simulations have demonstrated the formation of interarm spurs as well as spiral arm molecular clouds, provided the interstellar medium model includes a cold H  i phase. We find that the main effect of adding a magnetic field to these calculations is to inhibit the formation of structure in the disc. However, provided a cold phase is included, spurs and spiral arm clumps are still present if β≳ 0.1 in the cold gas. A caveat to the two-phase calculations though is that by assuming a uniform initial distribution, β≳ 10 in the warm gas, emphasizing that models with more consistent initial conditions and thermodynamics are required. Our simulations with only warm gas do not show such structure, irrespective of the magnetic field strength.
Furthermore, we find that the introduction of a cold H  i phase naturally produces the observed degree of disorder in the magnetic field, which is again absent from simulations using only warm gas. Whilst the global magnetic field follows the large-scale gas flow, the magnetic field also contains a substantial random component that is produced by the velocity dispersion induced in the cold gas during the passage through a spiral shock. Without any cold gas, the magnetic field in the warm phase remains relatively well ordered apart from becoming compressed in the spiral shocks. Our results provide a natural explanation for the observed high proportions of disordered magnetic field in spiral galaxies and we thus predict that the relative strengths of the random and ordered components of the magnetic field observed in spiral galaxies will depend on the dynamics of spiral shocks.     

Resolved shocks in clumpy media

We study the structure of shocks in clumpy media, using a multifluid formalism. As expected, shocks broaden as they weaken: for sufficiently weak shocks, no viscous subshock appears in the structure. This has significant implications for the survival of dense clouds in regions overrun by shocks in a wide range of astrophysical circumstances, from planetary nebulae to the nuclei of starburst galaxies.     

Driven waves in a two-fluid plasma

We study the physics of wave propagation in a weakly ionized plasma, as it applies to the formation of multifluid, magnetohydrodynamics (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 behaviour 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.     

On the meaning and inapplicability of the Zeldovich relations of magnetohydrodynamics

Considering a plasma with an initially weak large scale field subject to nonhelical turbulent stirring, Zeldovich (1957), for two‐dimensions, followed by others for three dimensions, have presented formulae of the form 〈b²〉 = f(R_M) . Such “Zeldovich relations” have sometimes been interpreted to provide steady‐state relations between the energy associated with the fluctuating magnetic field and that associated with a large scale or mean field multiplied by a function f that depends on spatial dimension and a magnetic Reynolds number R_M. Here we dissect the origin of these relations and pinpoint pitfalls that show why they are inapplicable to realistic, dynamical MHD turbulence and that they disagree with many numerical simulations. For 2D, we show that when the total magnetic field is determined by a vector potential, the standard Zeldovich relation applies only transiently, characterizing a maximum possible value that the field energy can reach before necessarily decaying. In 3D, we show that the standard Zeldovich relations are derived by balancing subdominant terms. In contrast, balancing the dominant terms shows that the fluctuating field can grow to a value independent of R_M and the initially imposed , as seen in numerical simulations. We also emphasize that these Zeldovich relations of nonhelical turbulence imply nothing about the amount mean field growth in a helical dynamo. In short, by re‐analyzing the origin of the Zeldovich relations, we highlight that they are inapplicable to realistic steady‐states of large R_M MHD turbulence. (© 2005 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Magnetic flux expulsion in powerful superbubble explosions and the α–Ω dynamo

The possibility of magnetic flux expulsion from the Galaxy in superbubble (SB) explosions, important for the α –Ω dynamo, is considered. Special emphasis is put on investigation of the downsliding of the matter from the top of the shell formed by the SB explosion, which is able to influence the kinematics of the shell. It is shown that either Galactic gravity or the development of the Rayleigh–Taylor instabilities in the shell limit the SB expansion, thus making magnetic flux expulsion impossible. The effect of cosmic rays in the shell on the sliding is considered, and it is shown that it is negligible compared with Galactic gravity. Thus the question of the possible mechanism of flux expulsion in the α –Ω dynamo remains open.     

Dynamical Molecular Regions in Interstellar Clouds

Work on three problems concerning the chemistry and dynamics of star forming fregions is addressed. Most of the work was done in collaboration with David A. Williams. A chemical model of core D ofTMC-1 is described; it is consistent with that dense core having evolved to its present physical state in the last 10⁵years. Slow-mode wave excitation by the nonlinear evolution of fast-mode waves is examined as a possible mechanism for the production of structures, including densecores, in molecular sources. Sufficient atomic hydrogen may be present at the births of some dense cores that the radiative association of O with H is important in the initiation of the gas phase chemistry. This revised version was published online in July 2006 with corrections to the Cover Date.     

Mean-field dynamo in partially ionized plasmas – I

There are several astrophysical situations where one needs to study the dynamics of magnetic flux in partially ionized turbulent plasmas. In a partially ionized plasma, the magnetic induction is subjected to the ambipolar diffusion and the Hall effect in addition to the usual resistive dissipation. In this paper, we initiate the study of the kinematic dynamo in a partially ionized turbulent plasma. The Hall effect arises from the treatment of the electrons and the ions as two separate fluids and the ambipolar diffusion due to the inclusion of neutrals as the third fluid. It is shown that these non-ideal effects modify the so-called α effect and the turbulent diffusion coefficient β in a rather substantial way. The Hall effect may enhance or quench the dynamo action altogether. The ambipolar diffusion brings in an α which depends on the mean magnetic field. The new correlations embodying the coupling of the charged fluids and the neutral fluid appear in a decisive manner. The turbulence is necessarily magnetohydrodynamic with new spatial and time-scales. The nature of the new correlations is demonstrated by taking the Alfvénic turbulence as an example.     

On non-axisymmetric magnetic equilibria in stars

In previous work, stable approximately axisymmetric equilibrium configurations for magnetic stars were found by numerical simulation. Here, I investigate the conditions under which more complex, non-axisymmetric configurations can form. I present numerical simulations of the formation of stable equilibria from turbulent initial conditions and demonstrate the existence of non-axisymmetric equilibria consisting of twisted flux tubes lying horizontally below the surface of the star, meandering around the star in random patterns. Whether such a non-axisymmetric equilibrium or a simple axisymmetric equilibrium forms depends on the radial profile of the strength of the initial magnetic field. The results could explain observations of non-dipolar fields on stars such as the B0.2 main-sequence star τ Sco or the pulsar 1E 1207.4-5209. The secular evolution of these equilibria due to Ohmic and buoyancy processes is also examined.     

Non-linear magnetic diffusivity in mean-field electrodynamics

We consider non-linear transport and drift processes caused by an inhomogeneous magnetic field in a turbulent fluid. The coefficients of magnetic diffusivity and drift velocity are calculated by making use of the second-order correlation approximation. Transport processes in the presence of a sufficiently strong magnetic field become anisotropic with larger diffusion rate and turbulent electrical resistivity across the field than along the field. Non-linear effects also lead to a drift of the magnetic field away from the regions with a higher magnetic energy.