No Arabic abstract
The ordered magnetic field observed via polarized 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 magneto-rotational 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 behavior 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 $alphaOmega$~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 galaxys 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.
Context: The interstellar medium (ISM) on all scales is full of structures that can be used as tracers of processes that feed turbulence. Aims: We used HI survey data to derive global properties of the angular power distribution of the local ISM. Methods: HI4PI observations on an nside = 1024 HEALPix grid and Gaussian components representing three phases, the cold, warm, and unstable lukewarm neutral medium (CNM, WNM, and LNM), were used for velocities $|v_{mathrm{LSR}}| leq 25$ kms. For high latitudes $|b| > 20deg$ we generated apodized maps. After beam deconvolution we fitted angular power spectra. Results: Power spectra for observed column densities are exceptionally well defined and straight in log-log presentation with 3D power law indices $gamma geq -3$ for the local gas. For intermediate velocity clouds (IVCs) we derive $gamma = -2.6$ and for high velocity clouds (HVCs) $gamma = -2.0$. Single-phase power distributions for the CNM, LNM, and WNM are highly correlated and shallow with $ gamma sim -2.5$ for multipoles $l leq 100$. Excess power from cold filamentary structures is observed at larger multipoles. The steepest single-channel power spectra for the CNM are found at velocities with large CNM and low WNM phase fractions. Conclusions: The phase space distribution in the local ISM is configured by phase transitions and needs to be described with three distinct different phases, being highly correlated but having distributions with different properties. Phase transitions cause locally hierarchical structures in phase space. The CNM is structured on small scales and is restricted in position-velocity space. The LNM as an interface to the WNM envelops the CNM. It extends to larger scales than the CNM and covers a wider range of velocities. Correlations between the phases are self-similar in velocity.
Supersonic turbulence is a large reservoir of suprathermal energy in the interstellar medium. Its dissipation, because it is intermittent in space and time, can deeply modify the chemistry of the gas. We further explore a hybrid method to compute the chemical and thermal evolution of a magnetized dissipative structure, under the energetic constraints provided by the observed properties of turbulence in the cold neutral medium. For the first time, we model a random line of sight by taking into account the relative duration of the bursts with respect to the thermal and chemical relaxation timescales of the gas. The key parameter is the turbulent rate of strain a due to the ambient turbulence. With the gas density, it controls the size of the dissipative structures, therefore the strength of the burst. For a large range of rates of strain and densities, the models of turbulent dissipation regions (TDR) reproduce the CH+ column densities observed in the diffuse medium and their correlation with highly excited H2. They do so without producing an excess of CH. As a natural consequence, they reproduce the abundance ratios of HCO+/OH and HCO+/H2O, and their dynamic range of about one order of magnitude observed in diffuse gas. Large C2H and CO abundances, also related to those of HCO+, are another outcome of the TDR models that compare well with observed values. The abundances and column densities computed for CN, HCN and HNC are one order of magnitude above PDR model predictions, although still significantly smaller than observed values.
This study is motivated by recent observations on ubiquitous interstellar density filaments and guided by modern theories of compressible magnetohydrodynamic (MHD) turbulence. The interstellar turbulence shapes the observed density structures. As the fundamental dynamics of compressible MHD turbulence, perpendicular turbulent mixing of density fluctuations entails elongated density structures aligned with the local magnetic field, accounting for low-density parallel filaments seen in diffuse atomic and molecular gas. The elongation of low-density parallel filaments depends on the turbulence anisotropy. When taking into account the partial ionization, we find that the minimum width of parallel filaments in the cold neutral medium and molecular clouds is determined by the neutral-ion decoupling scale perpendicular to magnetic field. In highly supersonic MHD turbulence in molecular clouds, both low-density parallel filaments due to anisotropic turbulent mixing and high-density filaments due to shock compression exist.
We study the evolution of dense clumps and provide argument that the existence of the clumps is not limited by the crossing time of the clump. We claim that the lifetimes of the clumps are determined by the turbulent motions on larger scale and predict the correlation of the clump lifetime and its column density. We use numerical simulations and successfully test this relation. In addition, we study the morphological asymmetry and the magnetization of the clumps as a function of their masses.
Understanding the physics of how stars form is a highly-prioritized goal of modern Astrophysics, in part because star formation is linked to both galactic dynamics on large scales and to the formation of planets on small scales. It is well-known that stars form from the gravitational collapse of molecular clouds, which are in turn formed out of the turbulent interstellar medium. Star formation is highly inefficient, with one of the likely culprits being the regulation against gravitational collapse provided by magnetic fields. Measurement of the polarized emission from interstellar dust grains, which are partially aligned with the magnetic field, provides a key tool for understanding the role these fields play in the star formation process. Over the past decade, much progress has been made by the most recent generation of polarimeters operating over a range of wavelengths (from the far-infrared through the millimeter part of the spectrum) and over a range of angular resolutions (from less than an arcsecond through fractions of a degree). Future developments in instrument sensitivity for ground-based, airborne, and space-borne polarimeters operating over range of spatial scales are critical for enabling revolutionary steps forward in our understanding of the magnetized turbulence from which stars are formed.