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)     

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.     

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.     

Hydrostatic equilibrium in a magnetized, warped Galactic disc

Hydrostatic equilibrium of the multiphase interstellar medium in the solar vicinity is reconsidered, with the regular and turbulent magnetic fields treated separately. The regular magnetic field strength required to support the gas is consistent with independent estimates, provided that energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a mid-plane value of B ₀=4 μG for the regular magnetic field near the Sun leads to more attractive models than B ₀=2 μG . The vertical profiles of both the regular and random magnetic fields contain disc and halo components, the parameters of which we have determined. The layer at 1≲| z |≲4 kpc can be overpressured and an outflow at a speed of about 50 km s⁻¹ may occur there, presumably associated with a Galactic fountain flow, if B ₀≃2 μG .
We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in z , in rough agreement with H  i observations in the outer Galaxy. This asymmetry may be a useful diagnostic of the details of the warping mechanism in the Milky Way and other galaxies. We find indications that gas and magnetic field pressures are different above and below the warped midplane in the outer Galaxy, and quantify the difference in terms of turbulent velocity and/or magnetic field strength.     

Milestones in the Observations of Cosmic Magnetic Fields

Magnetic fields are observed everywhere in the universe. In this review, we concentrate on the observational aspects of the magnetic fields of Galactic and extragalactic objects. Readers can follow the milestones in the observations of cosmic magnetic fields obtained from the most important tracers of magnetic fields, namely, the star-light polarization, the Zeeman effect, the rotation measures (RMs, hereafter) of extragalactic radio sources, the pulsar RMs, radio polarization observations, as well as the newly implemented sub-mm and mm polarization capabilities. The magnetic field of the Galaxy was first discovered in 1949 by optical polarization observations. The local magnetic fields within one or two kpc have been well delineated by starlight polarization data. The polarization observations of diffuse Galactic radio background emission in 1962 confirmed unequivocally the existence of a Galactic magnetic field. The bulk of the present information about the magnetic fields in the Galaxy comes from anal     

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.     

The interstellar magnetic field near the Galactic center

We review the current observational knowledge of the interstellar magnetic field within ∼150 pc ofthe Galactic center. We also discuss the various theoretical scenarios that have been put forward to explain the existing observations. Our critical overview leads to two important conclusions: (1) The interstellar magnetic field near the GC is approximately poloidal on average in the diffuse intercloud medium and approximately horizontal in dense interstellar clouds. (2) In the general intercloud medium, the field is relatively weak and probably close to equipartition with cosmic rays (B ∼ (6–20) μ G), but there exist a number of localized filaments where the field is much stronger (some filaments could possibly have B ≳ 1 mG). In dense interstellar clouds, the field is probably rather strong, with typical values ranging between a few 0.1 mG and a few mG (© 2010 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

Can the turbulent galactic dynamo generate large-scale magnetic fields?

Large-scale magnetic fields in galaxies are thought to be generated by a turbulent dynamo. However, the same turbulence also leads to a small-scale dynamo which generates magnetic noise at a more rapid rate. The efficiency of the large-scale dynamo depends on how this noise saturates. We examine this issue, taking into account ambipolar drift, which obtains in a galaxy with significant neutral gas. We argue as follows.
(i) The small-scale dynamo generated field does not fill the volume, but is concentrated into intermittent rope-like structures. The flux ropes are curved on the turbulent eddy scales. Their thickness is set by the diffusive scale determined by the effective ambipolar diffusion.
(ii) For a largely neutral galactic gas, the small-scale dynamo saturates, as a result of inefficient random stretching, when the peak field in a flux rope has grown to a few times the equipartition value.
(iii) The average energy density in the saturated small-scale field is subequipartition, since it does not fill the volume.
(iv) Such fields neither drain significant energy from the turbulence nor convert eddy motion of the turbulence on the outer scale into wave-like motion. The diffusive effects needed for the large-scale dynamo operation are then preserved until the large-scale field itself grows to near equipartition levels.     

A high sampling-density polarization study of the Southern Coalsack

We present a densely sampled map of visual polarimetry of stars in the direction of the Southern Coalsack dark cloud. Our sample consists of new polarimetric observations of 225 stars drawn from the spectrophotometric survey of Seidensticker, and an additional 173 stars, covering the surrounding areas of the cloud, taken from the literature. Because all the target stars have at least spectroscopic parallaxes, we can reliably investigate the spatial origins of the polarization, in three dimensions. We decompose the polarization into three components, due to (i) the wall of the local hot bubble, (ii) the Coalsack cloud and (iii) material in the Carina spiral arm. The polarization due to the Coalsack varies, both in alignment efficiency  ( p / A_V )  and in the dispersion in polarization angle, across the cloud. Using a simplified radiative transfer treatment we show that the measured polarization in background gas is significantly affected by foreground polarization, and specifically that the analysis of the Coalsack polarization must take the effects of the local hot bubble wall into consideration. Correcting for this effect as well as for the internal line-of-sight averaging in the Coalsack, we find, based on a Chandrasekhar–Fermi analysis, a plane-of-the-sky magnetic field for the Coalsack cloud of  〈 B _⊥〉= 93 ± 23 μG  . A systematic error, best described by a multiplicative factor between 0.5 and 1.5, additionally arises from radiative transfer effect uncertainties. We propose that this high value for the magnetic field in the cloud envelope is due to the fact that the Coalsack cloud is embedded in the hot interior of the Upper Centaurus–Lupus superbubble.     

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.     

Whistler wave cascades in solar wind plasma

Non-linear, three-dimensional, time-dependent fluid simulations of whistler wave turbulence are performed to investigate role of whistler waves in solar wind plasma turbulence in which characteristic turbulent fluctuations are characterized typically by the frequency and length-scales that are, respectively, bigger than ion gyrofrequency and smaller than ion gyroradius. The electron inertial length is an intrinsic length-scale in whistler wave turbulence that distinguishably divides the high-frequency solar wind turbulent spectra into scales smaller and bigger than the electron inertial length. Our simulations find that the dispersive whistler modes evolve entirely differently in the two regimes. While the dispersive whistler wave effects are stronger in the large-scale regime, they do not influence the spectral cascades which are describable by a Kolmogorov-like   k ^−7/3  spectrum. By contrast, the small-scale turbulent fluctuations exhibit a Navier–Stokes-like evolution where characteristic turbulent eddies exhibit a typical   k ^−5/3  hydrodynamic turbulent spectrum. By virtue of equipartition between the wave velocity and magnetic fields, we quantify the role of whistler waves in the solar wind plasma fluctuations.     

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.     

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.     

Superbubbles in magnetized interstellar media: blowout or confinement?

The importance of the interstellar magnetic field is studied in relation to the evolution of superbubbles with a three-dimensional (3D) numerical magnetohydrodynamical (MHD) simulation. A superbubble is a large supernova remnant driven by sequential supernova explosions in an OB association. Its evolution is affected by the density stratification in the galactic disc. After the superbubble size reaches 2–3 times the density scaleheight, it expands preferentially in the z -direction, until finally it can punch out a hole in the gas disc (blowout). On the other hand, the magnetic field running parallel to the galactic disc has the effect of preventing it from expanding in the direction perpendicular to the field. The density stratification and the magnetic fields have completely opposite effects on the evolution of the superbubble. We present results of a 3D MHD simulation in which both effects are included. As a result, it is concluded that when the magnetic field has a much larger scaleheight than the density, even for a model in which the bubble would blow out from the disc if the magnetic field were absent, a magnetic field with a strength of 5 μG can confine the bubble in | z |≲300 pc for ≃ 20 Myr (confinement). In a model in which the field strength decreases in the halo as B  ∝ ρ^1/2, the superbubble eventually blows out like a model with B  = 0 even if the magnetic field in the mid-plane is as strong as B  = 5 μG.     

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.     

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.