Dynamo saturation in direct simulations of the multi‐phase turbulent interstellar medium

The ordered magnetic field observed via polarised synchrotron emission in nearby disc galaxies can be explained by a mean‐field dynamo operating in the diffuse interstellar medium (ISM). Additionally, vertical‐flux initial conditions are potentially able to influence this dynamo via the occurrence of the magnetorotational instability (MRI). We aim to study the influence of various initial field configurations on the saturated state of the mean‐field dynamo. This is motivated by the observation that different saturation behaviour was previously obtained for different supernova rates. We perform direct numerical simulations (DNS) of three‐dimensional local boxes of the vertically stratified, turbulent interstellar medium, employing shearing‐periodic boundary conditions horizontally. Unlike in our previous work, we also impose a vertical seed magnetic field. We run the simulations until the growth of the magnetic energy becomes negligible. We furthermore perform simulations of equivalent 1D dynamo models, with an algebraic quenching mechanism for the dynamo coefficients. We compare the saturation of the magnetic field in the DNS with the algebraic quenching of a mean‐field dynamo. The final magnetic field strength found in the direct simulation is in excellent agreement with a quenched α) dynamo. For supernova rates representative of the Milky Way, field losses via a Galactic wind are likely responsible for saturation. We conclude that the relative strength of the turbulent and regular magnetic fields in spiral galaxies may depend on the galaxy's star formation rate. We propose that a mean field approach with algebraic quenching may serve as a simple sub‐grid scale model for galaxy evolution simulations including a prescribed feedback from magnetic fields. (© 2015 WILEY‐VCH Verlag GmbH & Co. KGaA, Weinheim)     

The effects of spiral arms on the multi-phase ISM

Statistical parameters of the ISM driven by thermal energy injectionsfrom supernova explosions have been obtained from 3D, nonlinear,magnetohydrodynamic, shearing-box simulations for spiral arm andinterarm regions. The density scale height obtained for the interarm regionsis 50% larger than within the spiral arms because of thehigher gas temperature. The filling factorof the hot gas is also significantly larger between the armsand depends sensitively on magnetic field strength.     

Magnetorotational supernovae. Magnetorotational instability. Jet formation

2D numerical simulations of magnetorotational (MR) supernova mechanism are described. It is shown that magnetic field is amplified due to the differential rotation after core collapse. When magnetic pressure reaches some level, a compression wave starts to move outwards. Moving along steeply decreasing density profile the compression wave transforms quickly into fast MHD shock. The magnetorotational instability (MRI) was found in our simulations. MRI leads to the exponential growth of the components of the magnetic field. The MRI significantly reduces MR supernova explosion time. Configuration of the initial magnetic field qualitatively defines the shape of MR supernova explosion. For the quadrupole-like initial poloidal field the MR supernova explosion develops mainly along equatorial plane, the dipole-like initial field results in MR supernova developing as mildly collimated jet along axis of rotation. The explosion energy of MR supernova found in our simulations is ∼0.5–0.6×10⁵¹ erg.     

Impact of magnetic field on the gas mass fraction of galaxy clusters

Magnetic fields have been observed in galaxy clusters with strengths of the order of  ~ μG. The non-thermal pressure exerted by magnetic fields also contributes to the total pressure in galaxy clusters and can in turn affect the estimates of the gas mass fraction, f_gas. In this paper, we have considered a central magnetic field strength of 5μG, motivated by observations and simulations of galaxy clusters. The profile of the magnetic field has also been taken from the results obtained from simulations and observations. The role of magnetic field has been taken into account in inferring the gas density distribution through the hydrostatic equilibrium condition (HSE) by including the magnetic pressure. We have found that the resultant gas mass fraction is smaller with magnetic field as compared to that without magnetic field. However, this decrease is dependent on the strength and the profile of the magnetic field. We have also determined the total mass using the NFW profile to check for the dependency of f_gas estimates on total mass estimators. From our analysis, we conclude that for the magnetic field strength that galaxy clusters seem to possess, the non-thermal pressure from magnetic fields has an impact of  ≈ 1 % on the gas mass fraction of galaxy clusters. However, with upcoming facilities like Square Kilometre Array (SKA), it can be further expected to improve with more precise observations of the magnetic field strength and profile in galaxy clusters, particularly in the interior region.     

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.     

Interfacial instabilities driven by self-gravity: first numerical results

The evolution of the interstellar medium (ISM) is driven by a variety of phenomena, including turbulence, shearing flows, magnetic fields and the thermal properties of the gas. Among the most important forces at work is self-gravity, which ultimately drives protostellar collapse. As part of an ongoing study of instabilities in the ISM, Hunter, Whitaker & Lovelace have discovered another process driven by self-gravity: the instability of an interface of discontinuous density. Theory predicts that this self-gravity driven interfacial instability persists in the static limit and in the absence of a constant background acceleration. Disturbances to a density interface are found to grow on a time-scale of the order of the free-fall time, even when the perturbation wavelength is much less than the Jeans length. Here we present the first numerical simulations of this instability. The theoretical growth rate is confirmed and the non-linear morphology displayed. The self-gravity interfacial instability is shown to be fundamentally different from the Rayleigh–Taylor instability, although both exhibit similar morphologies under the condition of a high density contrast, such as is commonly found in the ISM. Such instabilities are a possible mechanism by which observed features, such as the pillars of gas seen near the boundaries of interstellar clouds, are formed.     

What Drives Star Formation?

Current theoretical models for what drives star formation (especially low-mass star formation) are: (1) magnetic support of self-gravitating clouds with ambipolar diffusion removing support in cores and triggering collapse and (2) compressible turbulence forming self-gravitating clumps that collapse as soon as the turbulent cascade produces insufficient turbulent support. Observations of magnetic fields can distinguish between these two models because of different predictions in three areas: (1) magnetic field morphology, (2) the scaling of field strength with density and non-thermal velocities, and (3) the mass to magnetic flux ratio, M/Φ. We first discuss the techniques and limitations of methods for observing magnetic fields in star formation regions, then describe results for the L1544 prestellar core as an exemplar of the observational results. Application of the three tests leads to the following conclusions. The observational data show that both magnetic fields and turbulence are important in molecular cloud physics. Field lines are generally regular rather than chaotic, implying strong field strengths. But fields are not aligned with the minor axes of oblate spheroidal clouds, suggesting that turbulence is important. Field strengths appear to scale with non-thermal velocity widths, suggesting a significant turbulent support of clouds. Giant Molecular Clouds (GMCs) require mass accumulation over sufficiently large volumes that they would likely have an approximately critical M/Φ. Yet H I clouds are observed to be highly subcritical. If self-gravitating (molecular) clouds form with the subcritical M/Φ of H I clouds, the molecular clouds will be subcritical. However, the observations of molecular cloud cores suggest that they are approximately critical, with no direct evidence for subcritical molecular clouds or cloud envelopes. Hence, the observations remain inconclusive in deciding between the two extreme-case models of what drives star formation. What is needed to further advance our understanding of the role of magnetic fields in the star formation process are additional high sensitivity surveys of magnetic field strengths and other cloud properties in order to further refine the assessment of the importance of magnetic fields in molecular cores and envelopes.     

3D Solar Null Point Reconnection MHD Simulations

Numerical MHD simulations of 3D reconnection events in the solar corona have improved enormously over the last few years, not only in resolution, but also in their complexity, enabling more and more realistic modeling. Various ways to obtain the initial magnetic field, different forms of solar atmospheric models as well as diverse driving speeds and patterns have been employed. This study considers differences between simulations with stratified and non-stratified solar atmospheres, addresses the influence of the driving speed on the plasma flow and energetics, and provides quantitative formulas for mapping electric fields and dissipation levels obtained in numerical simulations to the corresponding solar quantities. The simulations start out from a potential magnetic field containing a null-point, obtained from a Solar and Heliospheric Observatory (SOHO) Michelson Doppler Imager (MDI) magnetogram magnetogram extrapolation approximately 8?hours before a C-class flare was observed. The magnetic field is stressed with a boundary motion pattern similar to?–?although simpler than?–?horizontal motions observed by SOHO during the period preceding the flare. The general behavior is nearly independent of the driving speed, and is also very similar in stratified and non-stratified models, provided only that the boundary motions are slow enough. The boundary motions cause a build-up of current sheets, mainly in the fan-plane of the magnetic null-point, but do not result in a flare-like energy release. The additional free energy required for the flare could have been partly present in non-potential form at the initial state, with subsequent additions from magnetic flux emergence or from components of the boundary motion that were not represented by the idealized driving pattern.     

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.     

Observational Evidence for a Double-Helix Structure in CMEs and Magnetic Clouds

We compare recent observations of a solar eruptive prominence as seen in extreme-UV light on 30 March 2010 by the Solar Dynamics Observatory (SDO) with the multi-tube model for interplanetary magnetic clouds (Osherovich, Fainberg, Stone, Geophys. Res. Lett. 26, 2597, 1999). Our model is based on an exact analytical solution of the plasma equilibrium with magnetic force balanced by a gradient of scalar gas pressure. Topologically, this solution describes two magnetic helices with opposite magnetic polarity embedded in a cylindrical magnetic flux tube that creates magnetic flux inequality between the two helices by enhancing one helix and suppressing the other. The magnetic field in this model is continuous everywhere and has a finite magnetic energy per unit length of the tube. These configurations have been introduced as MHD bounded states (Osherovich, Soln. Dannye 5, 70, 1975). Apparently, the SDO observations depict two non-equal magnetically interacting helices described by this analytical model. We consider magnetic and thermodynamic signatures of multiple magnetic flux ropes inside the same magnetic cloud, using in situ observations. The ratio of magnetic energy density to bulk speed solar wind energy density has been defined as a solar wind quasi-invariant (QI). We analyze the structure of the QI profile to probe the topology of the internal structure of magnetic clouds. From the superposition of 12 magnetically isolated clouds observed by Ulysses, we have found that the corresponding QI is consistent with our double helix model.     

Equatorial propagation of axisymmetric magnetogasdynamic shocks with thermal radiation,I

Similarity solutions, for one-dimensional unsteady of a perfect gas behind a spherical shock wave produced on account of a sudden explosion or driven out by an expanding piston including the effects of radiative cooling and an idealised azimuthal magnetic field, are studied. The shock is assumed to be strong and it is propagating into a transparent medium at rest with varying density. The magnetic field is proportional tor ^?1. The total energy of the wave is time dependent obeying a power law.     

Building the Fast Galactic Dynamo

In this paper we demonstrate the importance of cosmic rays for the dynamics of the interstellar medium. We present the first 3D-MHD numerical simulations of the Parker instability triggered by cosmic rays accelerated in randomly distributed supernova remnants. We show that in the presence of galactic rotation a net radial magnetic field is produced as a result of the cosmic ray injection and Coriolis force. This process provides a possibility of very efficient magnetic field amplification within the general frame of so called fast galactic dynamo proposed by Parker (1992).     

Systematic search for molecular clouds near supernova remnants as sources of very-high-energy γ-ray emission

Supernova remnants accelerate particles up to energies of at least 100 TeV as established by observations in very-high-energy γ-ray astronomy. Molecular clouds in their vicinity provide an increased amount of target material for proton-proton interaction and subsequent neutral pion decay into γ-rays of accelerated hadrons escaping the remnant. Therefore, these molecular clouds are potential γ-ray sources. The γ-ray emission from these clouds provides a unique environment to derive information on the propagation of very-high-energy particles through the interstellar medium as well as on the acceleration of hadrons in supernova remnants. Current Imaging Atmospheric Cherenkov Telescope systems are suitable to explore a large parameter space of the propagation properties depending on the age of the supernova remnant and the distance between the remnant and the nearby molecular cloud.In this paper we present our strategy and results of a systematic search for γ-ray emitting molecular clouds near supernova remnants which are potentially detectable with current experiments in the TeV energy range and explore the prospects of future experiments.     

The Environment of YSO Jets

It is commonly accepted that stars form in molecular clouds by the gravitational collapse of dense gas. However, it is precisely not the infalling but the outflowing material that is primarily observed. Outflow motions prevail around both low and high mass young stellar objects. We present here results from a family of self-similar models that could possibly help to understand this paradox. The models take into account the heating of the central protostar for the deflection and acceleration of the gas. The models make room for all the ingredients observed around the central objects, i.e. molecular outflows, fast jets, accretion disks and infalling envelopes. We suggest that radiative heating and magnetic field may ultimately be the main energy sources driving outflows for both low and high mass stars. The models show that the ambient medium surrounding the jet is unhomogeneous in density, velocity, magnetic field. Consequently, we suggest that jets and outflows have a prehistory that is inprinted in their environment, and that this should have direct consequences on the setting of jet numerical simulations.     

Effects of cosmic rays and density inhomogeneity of the circumstellar medium on the formation of a supernova shell

We consider the self-similar problem of a supernova explosion in a radially inhomogeneous medium by taking into account the generation of accelerated relativistic particles. The initial density of the medium is assumed to decrease with distance from the explosion center as a power law, ρ ₀ = A/r ^θ. We use a two-fluid approach in which the total pressure in the medium is the sum of the circumstellar gas pressure and the relativistic particle pressure. The relativistic particle pressure at the shock front is specified as an external parameter. This approach is applicable in the case where the diffusion coefficient of accelerated particles is small and the thickness of the shock front is much smaller than its radius. We have numerically solved a system of ordinary differential equations for the dimensionless quantities that describe the velocity and density behind the shock front as well as the nonrelativistic gas and relativistic particle pressures for various parameters of the inhomogeneity of the medium and various compression ratios of the medium at the shock front. We have established that the shock acceleration of cosmic rays affects most strongly the formation of a supernova shell (making it thinner) in a homogeneous circumstellar medium. A decrease in the circumstellar matter density with distance from the explosion center causes the effect of shock-accelerated relativistic particles on the supernova shell formation to weaken considerably. Inhomogeneity of the medium makes the shell thicker and less dense, while an increase in the compression ratio of the medium at the shock front causes the shell to become thinner and denser. As the relativistic particle density increases, the effect of circumstellar matter inhomogeneity on the shell formation becomes weaker.     

The Role of Magnetic Fields in Transient Seismic Emission Driven by Atmospheric Heating in Flares

Transient seismic emission in flares remains largely mysterious. Its discoverers proposed that seismic transients are driven by impulsive heating of the flaring chromosphere. Simulations of such heating show strong shocks, but these are damped by heavy radiative losses as they proceed downward. Because compression of the gas the shock enters both heats it and increases its density, the radiative losses increase radically with the strength of the shock, leaving doubt that sufficient energy can penetrate into the solar interior to explain helioseismic signatures. We note that simulations to date have no account for strong, inclined magnetic fields characteristic of transient-seismic-source environments. A strong horizontal magnetic field, for example, greatly increases the compressional modulus of the chromospheric medium, greatly reducing compression of the gas, hence radiative losses. Inclined magnetic fields, then, must be fundamental to the role of impulsive heating in transient seismic emission.     

Dissipation of magnetic flux and magnetic energy in partially ionized gases

Macroscopic equations of motion are used to derive several forms of the generalized Ohm's law for partially ionized ternary gases in magnetic fields, and a conductivity σ is defined that is independent of the magnetic field. A flux theorem is derived using a velocityu _H that can be defined to be the velocity of magnetic field lines;u _H is only slightly different from the velocity of the electron component of the gas. It is shown that σ is the conductivity relevant to the decay of magnetic flux through any surface moving everywhere with velocityu _H . The rate of increase of the thermal energy density of the gas arising through collisions between particles of different species can be resolved into Joule heating at the ratej ²/σ, wherej is the current density, and heating associated with ambipolar drift. The latter, contrary to what has been claimed by some authors, is not necessarily fully compensated by a decrease in the energy of the electromagnetic field. In many applications such compensation does occur, but it may not in interstellar clouds where large amounts of gravitational energy can be made available by collapse, and then both heating and an increase in electromagnetic field energy may occur.     

3D global simulations of galactic magnetic fields and gas flows

The evolution of three-dimensional (3D), dynamo excited galactic magnetic fields under the influence of a time-dependent gas flow in spiral arms is already well investigated. Our principal goal is to check how the dynamo-driven turbulent magnetic fields affect the gas flows. Numerical solutions of the full set of 3D MHD equations for dynamos in spiral galaxies are presented. Further we try to investigate the nonlinear evolution of magnetic instabilities in a global galactic model. The model includes differential rotation, eddy diffusivity and tensorial alpha-effect. In a first step the flow is driven by a prescribed gravitational potential. The vertical density stratification and the radial-azimutal spiral pattern are taken closely to observational data. We use a modified variant of the highly parallelized time-stepping ZeusMP code for the simulations of global galactic magnetic fields and gas flows. This revised version was published online in August 2006 with corrections to the Cover Date.     

Storm-intensity criteria for several classes of the driving interplanetary structures

In this paper we discuss the main-phase evolution of intense magnetic storms, associated with the passage of different interplanetary magnetic structures. It is shown that their evolution, driven by intense magnetic fields in the sheath region just behind interplanetary shocks, evolves faster (implying physically different magnetospheric configurations) than that associated with intense magnetic fields in the ejecta itself and in corotating streams. The estimated ring-current injection rate for the main phase of intense magnetic storms caused by sheath fields is ～10 times greater than the estimated injection rate for N–S magnetic clouds. Based on these results, we propose storm-intensity criteria for several classes of the driving interplanetary structures. The time necessary to reach a Dst/SYM index threshold level is an important parameter for a space weather forecast.     

Hydrodynamic and turbulent motions in the galactic disk

Statistics in absorption 21-cm data show two main types of clouds at low galactic latitudes: dense small clouds, many of them with molecular cores, with dispersions σ≈1.5 km s⁻¹ and large clouds forming the fine features of the spiral arms (the shingle like features) with a dispersion range α≈3–4 km s⁻¹. Sizes and dispersions of both types of clouds are compatible with the Kolmogorov law of turbulence: σ∞d ^1/3. The large clouds forming the shingle-like features can be considered as the largest clouds of a Kolmogorov spectrum (the initial vortices), or as the hydrodynamic features with minimum sizes in the Galaxy. In order to define hydrodynamic motions in the same sense as given by Ogrodnikov (1965) we use here the tensorial form of the Helmholtz theorem to obtain an approximation for the hydrodynamic motions depending on distances and seen from the local standard of rest:V _r ∞r. The intermediate range of sizes between turbulent motions and hydrodynamic motions is 100<d<300 pc which is also the range of sizes of the large clouds forming the fine features of the spiral arms. A classification on of motions in the Galaxy is postulated: (a) a basic rotation motion given by an smooth unperturbed curveΘ _b (R) associated to the old disk population. (b) Systematic motions of the spiral arms. (c) Systematic motions in the fine structure of the arms. For scale sizes smaller than these fine features one has turbulent motions according to the Kolmogorov law. The densities and sizes of the turbulent clouds behave asn _H ∝d ⁻² in a range of sizes 7 pc<d<300 pc. The obtained gas densities of the clouds are confirmed with the dust densities from photometric studies. The conditions for gravitational binding of the clouds are analyzed. Factors as the geometry and the magnetic field within the clouds increases the critic densities for gravitational binding. When we consider these factors we find that the wide component clouds have densities below such a critical value. The narrow component clouds have densities similar or above the critical value; but the real fraction of collapsing clouds remains unknown as far as the factor of geometry and the inner magnetic field of each cloud are not determinated.