The presence of magnetic fields in many astrophysical objects is due to dynamo action, whereby a part of the kinetic energy is converted into magnetic energy. A turbulent dynamo that produces magnetic field structures on the same scale as the turbulent flow is known as the fluctuation dynamo. We use numerical simulations to explore the nonlinear, statistically steady state of the fluctuation dynamo in driven turbulence. We demonstrate that as the magnetic field growth saturates, its amplification and diffusion are both affected by the back-reaction of the Lorentz force upon the flow. The amplification of the magnetic field is reduced due to stronger alignment between the velocity field, magnetic field, and electric current density. Furthermore, we confirm that the amplification decreases due to a weaker stretching of the magnetic field lines. The enhancement in diffusion relative to the field line stretching is quantified by a decrease in the computed local value of the magnetic Reynolds number. Using the Minkowski functionals, we quantify the shape of the magnetic structures produced by the dynamo as magnetic filaments and ribbons in both kinematic and saturated dynamos and derive the scalings of the typical length, width, and thickness of the magnetic structures with the magnetic Reynolds number. We show that all three of these magnetic length scales increase as the dynamo saturates. The magnetic intermittency, strong in the kinematic dynamo (where the magnetic field strength grows exponentially) persists in the statistically steady state, but intense magnetic filaments and ribbons are more volume-filling.

The origin of large-scale magnetic fields, detected in some low-mass (dwarf and irregular) galaxies via polarised synchrotron emission and Faraday rotation, remained unexplained for a long time. We suggest that mean-field dynamo can be active in galaxies of this class despite their slow rotation because their discs are relatively thick. Earlier assessments of the possibility of the mean-field dynamo action in low-mass galaxies relied on estimates applicable to $\textit{thin}$ discs, such as those in massive spiral galaxies. Using both order-of-magnitude estimates and numerical solutions, we show that the strength of differential rotation required to amplify magnetic field reduces as the aspect ratio of the galactic gas layer increases. Thus, the puzzle of the origin of large-scale magnetic fields appears to be solved. As a result, this class of galaxies provides a new ground for testing our understanding of galactic magnetism.

We explore the effect of magnetic fields on the vertical distribution and multiphase structure of the supernova-driven interstellar medium (ISM) in simulations that admit dynamo action. As the magnetic field is amplified to become dynamically significant, gas becomes cooler and its distribution in the disc becomes more homogeneous. We attribute this to magnetic quenching of vertical velocity, which leads to a decrease in the cooling length of hot gas. A non-monotonic vertical distribution of the large-scale magnetic field strength, with the maximum at |z| $\approx$ 300 pc causes a downward pressure gradient below the maximum which acts against outflow driven by SN explosions, while it provides pressure support above the maximum.

Fluid approximations to cosmic ray (CR) transport are often preferred to kinetic descriptions in studies of the dynamics of the interstellar medium (ISM) of galaxies, because they allow simpler analytical and numerical treatments. Magnetohydrodynamic (MHD) simulations of the ISM usually incorporate CR dynamics as an advection-diffusion equation for CR energy density, with anisotropic, magnetic field-aligned diffusion with the diffusive flux assumed to obey Fick's law. We compare test-particle and fluid simulations of CRs in a random magnetic field. We demonstrate that a non-Fickian prescription of CR diffusion, which corresponds to the telegraph equation for the CR energy density, can be easily calibrated to match the test particle simulations with great accuracy. In particular, we consider a random magnetic field in the fluid simulation that has a lower spatial resolution than that used in the particle simulation to demonstrate that an appropriate choice of the diffusion tensor can account effectively for the unresolved (subgrid) scales of the magnetic field. We show that the characteristic time which appears in the telegraph equation can be physically interpreted as the time required for the particles to reach a diffusive regime and we stress that the Fickian description of the CR fluid is unable to describe complex boundary or initial conditions for the CR energy flux.

The IMAGINE Consortium aims to bring modeling of the magnetic field of the Milky Way to a next level, by using Bayesian inference. IMAGINE includes an open-source modular software pipeline that optimizes parameters in a user-defined Galactic magnetic field model against various selected observational datasets. Bayesian priors can be added as external probabilistic constraints of the model parameters. These conference proceedings describe the science goals of the IMAGINE Consortium, the software pipeline and its inputs, viz observational data sets, Galactic magnetic field models, and Bayesian priors.

A convenient representation of the structure of the large-scale galactic magnetic field is required for the interpretation of polarization data in the sub-mm and radio ranges, in both the Milky Way and external galaxies. We develop a simple and flexible approach to construct parametrised models of the large-scale magnetic field of the Milky Way and other disc galaxies, based on physically justifiable models of magnetic field structure. The resulting models are designed to be optimised against available observational data. Representations for the large-scale magnetic fields in the flared disc and spherical halo of a disc galaxy were obtained in the form of series expansions whose coefficients can be calculated from observable or theoretically known galactic properties. The functional basis for the expansions is derived as eigenfunctions of the mean-field dynamo equation or of the vectorial magnetic diffusion equation. The solutions presented are axially symmetric but the approach can be extended straightforwardly to non-axisymmetric cases. The magnetic fields are solenoidal by construction, can be helical, and are parametrised in terms of observable properties of the host object, such as the rotation curve and the shape of the gaseous disc. The magnetic field in the disc can have a prescribed number of field reversals at any specified radii. Both the disc and halo magnetic fields can separately have either dipolar or quadrupolar symmetry. The model is implemented as a publicly available software package GalMag which allows, in particular, the computation of the synchrotron emission and Faraday rotation produced by the model's magnetic field. The model can be used in interpretations of observations of magnetic fields in the Milky Way and other spiral galaxies, in particular as a prior in Bayesian analyses. (Abridged.)

We study the cosmic evolution of the magnetic fields of a large sample of spiral galaxies in a cosmologically representative volume by employing a semi-analytic galaxy formation model and numerical dynamo solver in tandem. We start by deriving time- and radius-dependent galaxy properties using the galform galaxy formation model, which are then fed into the nonlinear mean-field dynamo equations. These are solved to give the large-scale (mean) field as a function of time and galactocentric radius for a thin disc, assuming axial symmetry. A simple prescription for the evolution of the small-scale (random) magnetic field component is also adopted. We find that, while most massive galaxies are predicted to have large-scale magnetic fields at redshift z=0, a significant fraction of them are expected to contain negligible large-scale field. Our model indicates that, for most of the galaxies containing large-scale magnetic fields today, the mean-field dynamo becomes active at z<3. We compute the radial profiles of pitch angle, and find broad agreement with observational data for nearby galaxies.

We argue that pitch angles of the azimuthally averaged large-scale or mean magnetic fields in nearby spiral galaxies inferred from observations can tentatively be explained with simple galactic dynamo models. Agreement is not perfect, but is reasonable considering the uncertainty in dynamo parameters.

Interpretations of synchrotron observations often assume a tight correlation between magnetic and cosmic ray energy densities. We examine this assumption using both test-particle simulations of cosmic rays and MHD simulations which include cosmic rays as a diffusive fluid. We find no spatial correlation between the cosmic rays and magnetic field energy densities at turbulent scales. Moreover, the cosmic ray number density and magnetic field energy density are statistically independent. Nevertheless, the cosmic ray spatial distribution is highly inhomogeneous, especially at low energies because the particles are trapped between random magnetic mirrors. These results can significantly change the interpretation of synchrotron observations and thus our understanding of the strength and structure of magnetic fields in the Milky Way and nearby spiral galaxies.

To constrain the large-scale magnetic field strengths in cosmologically distant galax- ies, we derive the probability distribution function of Faraday rotation measure (RM) when random lines of sight pass through a sample of disc galaxies, with axisymmetric large-scale magnetic fields. We find that the width of the RM distribution of the galaxy sample is directly related to the mean large-scale field strength of the galaxy population, provided the dispersion within the sample is lower than the mean value. In the absence of additional constraints on parameters describing the magneto-ionic medium of the intervening galaxies, and in the situation where RMs produced in the intervening galaxies have already been statistically isolated from other RM contributions along the lines of sight, our simple model of the magneto-ionic medium in disc galaxies suggests that the mean large-scale magnetic field of the population can be measured to within ~ 50% accuracy.

The role of magnetic fields in the multi-phase interstellar medium (ISM) is explored using magnetohydrodynamic (MHD) simulations that include energy injection by supernova (SN) explosions and allow for dynamo action. Apart from providing additional pressure support of the gas layer, magnetic fields reduce the density contrast between the warm and hot gas phases and quench galactic outflows. A dynamo-generated, self-consistent large-scale magnetic field affects the ISM differently from an artificially imposed, unidirectional magnetic field.

In this white paper we introduce the IMAGINE Consortium and its scientific background, goals and structure. Our purpose is to coordinate and facilitate the efforts of a diverse group of researchers in the broad areas of the interstellar medium, Galactic magnetic fields and cosmic rays, and our goal is to develop more comprehensive insights into the structures and roles of interstellar magnetic fields and their interactions with cosmic rays. To achieve a higher level of self-consistency, depth and rigour can only be achieved by the coordinated efforts of experts in diverse areas of astrophysics involved in observational, theoretical and numerical work. We present our view of the present status of this topic, identify its key unsolved problems and suggest a strategy that will underpin our work. The backbone of the consortium is the Interstellar MAGnetic field INference Engine, a publicly available Bayesian platform that employs robust statistical methods to explore the multi-dimensional likelihood space using any number of modular inputs. It provides an interpretation and modelling framework that has the power and flexibility to include a variety of observational, theoretical and numerical lines of evidence into a self-consistent and comprehensive picture of the thermal and non-thermal interstellar media. An important innovation is that a consistent understanding of the phenomena that are directly or indirectly influenced by the Galactic magnetic field, such as the deflection of ultra-high energy cosmic rays or extragalactic backgrounds, is made an integral part of the modelling. The IMAGINE Consortium, which is informal by nature and open to new participants, hereby presents a methodological framework for the modelling and understanding of Galactic magnetic fields that is available to all communities whose research relies on a state-of-the-art solution to this problem. (Abridged.)

We apply Gaussian smoothing to obtain mean density, velocity, magnetic and energy density fields in simulations of the interstellar medium based on three-dimensional magnetohydrodynamic equations in a shearing box $1\times1\times2 \, \rm{kpc}$ in size. Unlike alternative averaging procedures, such as horizontal averaging, Gaussian smoothing retains the three-dimensional structure of the mean fields. Although Gaussian smoothing does not obey the Reynolds rules of averaging, physically meaningful central statistical moments are defined as suggested by Germano (1992). We discuss methods to identify an optimal smoothing scale $\ell$ and the effects of this choice on the results. From spectral analysis of the magnetic, density and velocity fields, we find a suitable smoothing length for all three fields, of $\ell \approx 75 \, \rm{pc}$. We discuss the properties of third-order statistical moments in fluctuations of kinetic energy density in compressible flows and suggest their physical interpretation. The mean magnetic field, amplified by a mean-field dynamo, significantly alters the distribution of kinetic energy in space and between scales, reducing the magnitude of kinetic energy at intermediate scales. This intermediate-scale kinetic energy is a useful diagnostic of the importance of SN-driven outflows.

The predictions of mean-field electrodynamics can now be probed using direct numerical simulations of random flows and magnetic fields. When modelling astrophysical MHD, it is important to verify that such simulations are in agreement with observations. One of the main challenges in this area is to identify robust \it{quantitative} measures to compare structures found in simulations with those inferred from astrophysical observations. A similar challenge is to compare quantitatively results from different simulations. Topological data analysis offers a range of techniques, including the Betti numbers and persistence diagrams, that can be used to facilitate such a comparison. After describing these tools, we first apply them to synthetic random fields and demonstrate that, when the data are standardized in a straightforward manner, some topological measures are insensitive to either large-scale trends or the resolution of the data. Focusing upon one particular astrophysical example, we apply topological data analysis to HI observations of the turbulent interstellar medium (ISM) in the Milky Way and to recent MHD simulations of the random, strongly compressible ISM. We stress that these topological techniques are generic and could be applied to any complex, multi-dimensional random field.

Deriving the Faraday rotation measure (RM) of quasar absorption line systems, which are tracers of high-redshift galaxies intervening background quasars, is a powerful tool for probing magnetic fields in distant galaxies. Statistically comparing the RM distributions of two quasar samples, with and without absorption line systems, allows one to infer magnetic field properties of the intervening galaxy population. Here, we have derived the analytical form of the probability distribution function (PDF) of RM produced by a single galaxy with an axisymmetric large-scale magnetic field. We then further determine the PDF of RM for one random sight line traversing each galaxy in a population with a large-scale magnetic field prescription. We find that the resulting PDF of RM is dominated by a Lorentzian with a width that is directly related to the mean axisymmetric large-scale field strength $\langle B_0 \rangle$ of the galaxy population if the dispersion of $B_0$ within the population is smaller than $\langle B_0 \rangle$. Provided that RMs produced by the intervening galaxies have been successfully isolated from other RM contributions along the line of sight, our simple model suggests that $\langle B_0 \rangle$ in galaxies probed by quasar absorption line systems can be measured within $\approx50$ per cent accuracy without additional constraints on the magneto-ionic medium properties of the galaxies. Finally, we discuss quasar sample selection criteria that are crucial to reliably interpret observations, and argue that within the limitations of the current database of absorption line systems, high-metallicity damped Lyman-$\alpha$ absorbers are best suited to study galactic dynamo action in distant disc galaxies.

The interstellar medium (ISM) is a magnetised system in which transonic or supersonic turbulence is driven by supernova explosions. This leads to the production of intermittent, filamentary structures in the ISM gas density, whilst the associated dynamo action also produces intermittent magnetic fields. The traditional theory of random functions, restricted to second-order statistical moments (or power spectra), does not adequately describe such systems. We apply topological data analysis (TDA), sensitive to all statistical moments and independent of the assumption of Gaussian statistics, to the gas density fluctuations in a magnetohydrodynamic (MHD) simulation of the multi-phase ISM. This simulation admits dynamo action, so produces physically realistic magnetic fields. The topology of the gas distribution, with and without magnetic fields, is quantified in terms of Betti numbers and persistence diagrams. Like the more standard correlation analysis, TDA shows that the ISM gas density is sensitive to the presence of magnetic fields. However, TDA gives us important additional information that cannot be obtained from correlation functions. In particular, the Betti numbers per correlation cell are shown to be physically informative. Magnetic fields make the ISM more homogeneous, reducing the abundance of both isolated gas clouds and cavities, with a stronger effect on the cavities. Remarkably, the modification of the gas distribution by magnetic fields is captured by the Betti numbers even in regions more than 300 pc from the midplane, where the magnetic field is weaker and correlation analysis fails to detect any signatures of magnetic effects.

We apply correlation analysis to random fields in numerical simulations of the supernova-driven interstellar medium (ISM) with the magnetic field produced by dynamo action. We solve the thermo-magnetohydrodynamic (MHD) equations in a shearing, Cartesian box representing a local region of the ISM, subject to thermal and kinetic energy injection by supernova explosions, and parameterized optically-thin radiative cooling. We consider the cold, warm\, and hot phases of the ISM separately; the analysis mostly considers the warm gas, which occupies the bulk of the domain. Various physical variables have different correlation lengths in the warm phase: $40 pc$, $50 pc$, and $60 pc$ for random magnetic field, density, and velocity, respectively, in the midplane. The correlation time of the random velocity is comparable to the eddy turnover time, about $10^7 yr$, although it may be shorter in regions with higher star formation rate. The random magnetic field is anisotropic, with the standard deviations of its components $b_x/b_y/b_z$ having approximate ratios $0.5/0.6/0.6$ in the midplane. The anisotropy is attributed to the global velocity shear from galactic differential rotation, and locally inhomogeneous outflow to the galactic halo. The correlation length of Faraday depth along the $z$-axis, $120 pc$, is greater than for electron density, $60 - 90 pc$, and vertical magnetic field, $60 pc$. Such comparisons may be sensitive to the orientation of the line of sight. Uncertainties of the structure functions of synchrotron intensity rapidly increase with the scale. This feature is hidden in power spectrum analysis, which can undermine the usefulness of power spectra for detailed studies of interstellar turbulence.

Synchrotron radiation from cosmic rays is a key observational probe of the galactic magnetic field. Interpreting synchrotron emission data requires knowledge of the cosmic ray number density, which is often assumed to be in energy equipartition (or otherwise tightly correlated) with the magnetic field energy. However, there is no compelling observational or theoretical reason to expect such tight correlation to hold across all scales. We use test particle simulations, tracing the propagation of charged particles (protons) through a random magnetic field, to study the cosmic ray distribution at scales comparable to the correlation scale of the turbulent flow in the interstellar medium ($\simeq 100\,{\rm pc}$ in spiral galaxies). In these simulations, we find that there is no spatial correlation between the cosmic ray number density and the magnetic field energy density. In fact, their distributions are approximately statistically independent. We find that low-energy cosmic rays can become trapped between magnetic mirrors, whose location depends more on the structure of the field lines than on the field strength.

The propagation of cosmic rays in turbulent magnetic fields is a diffusive process driven by the scattering of the charged particles by random magnetic fluctuations. Such fields are usually highly intermittent, consisting of intense magnetic filaments and ribbons surrounded by weaker, unstructured fluctuations. Studies of cosmic ray propagation have largely overlooked intermittency, instead relying on Gaussian random magnetic fields. Using test particle simulations, we investigate cosmic ray diffusivity in intermittent, dynamo-generated magnetic fields. The results are compared with those obtained from non-intermittent magnetic fields having identical power spectra. The presence of magnetic intermittency significantly enhances cosmic ray diffusion over a wide range of particle energies. We demonstrate that the results can be interpreted in terms of a correlated random walk.

Current astrophysical models of the interstellar medium assume that small scale variation and noise can be modelled as Gaussian random fields or simple transformations thereof, such as lognormal. We use topological methods to investigate this assumption for three regions of the southern sky. We consider Gaussian random fields on two-dimensional lattices and investigate the expected distribution of topological structures quantified through Betti numbers. We demonstrate that there are circumstances where differences in topology can identify differences in distributions when conventional marginal or correlation analyses may not. We propose a non-parametric method for comparing two fields based on the counts of topological features and the geometry of the associated persistence diagrams. When we apply the methods to the astrophysical data, we find strong evidence against a Gaussian random field model for each of the three regions of the interstellar medium that we consider. Further, we show that there are topological differences at a local scale between these different regions.

We explore the effects of the multi-phase structure of the interstellar medium (ISM) on galactic magnetic fields. Basing our analysis on compressible magnetohydrodynamic (MHD) simulations of supernova-driven turbulence in the ISM, we investigate the properties of both the mean and fluctuating components of the magnetic field. We find that the mean magnetic field preferentially resides in the warm phase and is generally absent from the hot phase. The fluctuating magnetic field does not show such pronounced sensitivity to the multi-phase structure.

Mean-field galactic dynamo theory is the leading theory to explain the prevalence of regular magnetic fields in spiral galaxies, but its systematic comparison with observations is still incomplete and fragmentary. Here we compare predictions of mean-field dynamo models to observational data on magnetic pitch angle and the strength of the mean magnetic field. We demonstrate that a standard $\alpha^2\Omega$ dynamo model produces pitch angles of the regular magnetic fields of nearby galaxies that are reasonably consistent with available data. The dynamo estimates of the magnetic field strength are generally within a factor of a few of the observational values. Reasonable agreement between theoretical and observed pitch angles generally requires the turbulent correlation time $\tau$ to be in the range 10-20 Myr, in agreement with standard estimates. Moreover, good agreement also requires that the ratio of the ionized gas scale height to root-mean-square turbulent velocity increases with radius. Our results thus widen the possibilities to constrain interstellar medium (ISM) parameters using observations of magnetic fields. This work is a step toward systematic statistical tests of galactic dynamo theory. Such studies are becoming more and more feasible as larger datasets are acquired using current and up-and-coming instruments.

The propagation of charged particles, including cosmic rays, in a partially ordered magnetic field is characterized by a diffusion tensor whose components depend on the particle's Larmor radius $R_L$ and the degree of order in the magnetic field. Most studies of the particle diffusion presuppose a scale separation between the mean and random magnetic fields (e.g., there being a pronounced minimum in the magnetic power spectrum at intermediate scales). Scale separation is often a good approximation in laboratory plasmas, but not in most astrophysical environments such as the interstellar medium (ISM). Modern simulations of the ISM have numerical resolution of order 1 pc, so the Larmor radius of the cosmic rays that dominate in energy density is at least $10^{6}$ times smaller than the resolved scales. Large-scale simulations of cosmic ray propagation in the ISM thus rely on oversimplified forms of the diffusion tensor. We take the first steps towards a more realistic description of cosmic ray diffusion for such simulations, obtaining direct estimates of the diffusion tensor from test particle simulations in random magnetic fields (with the Larmor radius scale being fully resolved), for a range of particle energies corresponding to $10^{-2}\lesssim R_L/l_c \lesssim 10^{3}$, where $l_c$ is the magnetic correlation length. We obtain explicit expressions for the cosmic ray diffusion tensor for $R_L/l_c \ll 1$, that might be used in a sub-grid model of cosmic ray diffusion. The diffusion coefficients obtained are closely connected with existing transport theories that include the random walk of magnetic lines.

We examine the evolution of the Parker instability in galactic disks using 3D numerical simulations. We consider a local Cartesian box section of a galactic disk, where gas, magnetic fields and cosmic rays are all initially in a magnetohydrostatic equilibrium. This is done for different choices of initial cosmic ray density and magnetic field. The growth rates and characteristic scales obtained from the models, as well as their dependences on the density of cosmic rays and magnetic fields, are in broad agreement with previous (linearized, ideal) analytical work. However, this non-ideal instability develops a multi-modal 3D structure, which cannot be quantitatively predicted from the earlier linearized studies. This 3D signature of the instability will be of importance in interpreting observations. As a preliminary step towards such interpretations, we calculate synthetic polarized intensity and Faraday rotation measure maps, and the associated structure functions of the latter, from our simulations; these suggest that the correlation scales inferred from rotation measure maps are a possible probe for the cosmic ray content of a given galaxy. Our calculations highlight the importance of cosmic rays in these measures, making them an essential ingredient of realistic models of the interstellar medium.

Isotropic and anisotropic wavelet transforms are used to decompose the images of the spiral galaxy M83 in various tracers to quantify structures in a range of scales from 0.2 to 10 kpc. We used radio polarization observations at {\lambda}6 cm and 13 cm obtained with the VLA, Effelsberg and ATCA telescopes and APEX sub-mm observations at 870 {\mu}m, which are first published here, together with maps of the emission of warm dust, ionized gas, molecular gas, and atomic gas. The spatial power spectra are similar for the tracers of dust, gas, and total magnetic field, while the spectra of the ordered magnetic field are significantly different. The wavelet cross-correlation between all material tracers and total magnetic field is high, while the structures of the ordered magnetic field are poorly correlated with those of other tracers. -- The magnetic field configuration in M83 contains pronounced magnetic arms. Some of them are displaced from the corresponding material arms, while others overlap with the material arms. The magnetic field vectors at {\lambda}6 cm are aligned with the outer material arms, while significant deviations occur in the inner arms and in the bar region, possibly due to non-axisymmetric gas flows. Outside the bar region, the typical pitch angles of the material and magnetic spiral arms are very close to each other at about 10{\deg}. The typical pitch angle of the magnetic field vectors is about 20{\deg} larger than that of the material spiral arms. One of the main magnetic arms in M83 is displaced from the gaseous arms, while the other main arm overlaps a gaseous arm. We propose that a regular spiral magnetic field generated by a mean-field dynamo is compressed in material arms and partly aligned with them. The interaction of galactic dynamo action with a transient spiral pattern is a promising mechanism for producing such complicated spiral patterns as in M83.

A framework is introduced for coupling the evolution of galactic magnetic fields sustained by the mean-field dynamo with the formation and evolution of galaxies in cold dark matter cosmology. Estimates of the steady-state strength of the large-scale and turbulent magnetic fields from mean-field and fluctuation dynamo models are used together with galaxy properties predicted by semi-analytic models of galaxy formation for a population of spiral galaxies. We find that the field strength is mostly controlled by the evolving gas content of the galaxies. Thus, because of the differences in the implementation of the star formation law, feedback from supernovae and ram-pressure stripping, each of the galaxy formation models considered predicts a distribution of field strengths with unique features. The most prominent of them is the difference in typical magnetic fields strengths obtained for the satellite and central galaxy populations as well as the typical strength of the large-scale magnetic field in galaxies of different mass.

Both observations and modelling of magnetic fields in the diffuse interstellar gas of spiral galaxies are well developed but the theory has been confronted with observations for only a handful of individual galaxies. There is now sufficient data to consider statistical properties of galactic magnetic fields. We have collected data from the literature on the magnetic fields and interstellar media (ISM) of 20 spiral galaxies, and tested for various physically motivated correlations between magnetic field and ISM parameters. Clear correlations emerge between the total magnetic field strength and molecular gas density as well as the star formation rate. The magnetic pitch angle exhibits correlations with the total gas density, the star formation rate and the strength of the axisymmetric component of the mean magnetic field. The total and mean magnetic field strengths exhibit noticeable degree of correlation, suggesting a universal behaviour of the degree of order in galactic magnetic fields. We also compare the predictions of galactic dynamo theory to observed magnetic field parameters and identify directions in which theory and observations might be usefully developed.

We show that both aspect ratios of randomly oriented triaxial ellipsoids (representing isosurfaces of an isotropic 3D random field) can be determined from a single 2D cross-section of their sample using the probability distribution of the filamentarity F of the structures seen in the cross-section (F=0 for a circle and F=1 for a line). The probability distribution of F has a robust form with a sharp maximum and truncation that are sensitive to the ellipsoids' aspect ratios. We show that the aspect ratios of triaxial ellipsoids with randomly distributed dimensions can still be recovered from the probability distribution of F. This method is applicable to many shape recognition and classification problems, here illustrated with neutral hydrogen density in the turbulent interstellar medium of the Milky Way. The gas distribution is shown to be filamentary with the mean axis ratio 1:2:20.

Galactic magnetic arms have been observed between the gaseous arms of some spiral galaxies; their origin remains unclear. We suggest that magnetic spiral arms can be naturally generated in the interarm regions because the galactic fountain flow or wind is likely to be weaker there than in the arms. Galactic outflows lead to two countervailing effects: removal of small-scale magnetic helicity, which helps to avert catastrophic quenching of the dynamo, and advection of the large-scale magnetic field, which suppresses dynamo action. For realistic galactic parameters, the net consequence of outflows being stronger in the gaseous arms is higher saturation large-scale field strengths in the interarm regions as compared to in the arms. By incorporating rather realistic models of spiral structure and evolution into our dynamo models, an interlaced pattern of magnetic and gaseous arms can be produced.

We compare various models and approximations for non-linear mean-field dynamos in disc galaxies to assess their applicability and accuracy, and thus to suggest a set of simple solutions suitable to model the large-scale galactic magnetic fields in various contexts. The dynamo saturation mechanisms considered are the magnetic helicity balance involving helicity fluxes (the dynamical $\alpha$-quenching) and an algebraic $\alpha$-quenching. The non-linear solutions are then compared with the marginal kinematic and asymptotic solutions. We also discuss the accuracy of the no-$z$ approximation. Although these tools are very different in the degree of approximation and hence complexity, they all lead to remarkably similar solutions for the mean magnetic field. In particular, we show that the algebraic $\alpha$-quenching non-linearity can be obtained from a more physical dynamical $\alpha$-quenching model in the limit of nearly azimuthal magnetic field. This suggests, for instance, that earlier results on galactic disc dynamos based on the simple algebraic non-linearity are likely to be reliable, and that estimates based on simple, even linear models are often a good starting point. We suggest improved no-$z$ and algebraic $\alpha$-quenching models, and also incorporate galactic outflows into a simple analytical dynamo model to show that the outflow can produce leading magnetic spirals near the disc surface. The simple dynamo models developed are applied to estimate the magnetic pitch angle and the arm-interarm contrast in the saturated magnetic field strength for realistic parameter values.

We employ an analytical model that incorporates both wavelength-dependent and wavelength-independent depolarization to describe radio polarimetric observations of polarization at $\lambda \lambda \lambda \, 3.5, 6.2, 20.5$ cm in M51 (NGC 5194). The aim is to constrain both the regular and turbulent magnetic field strengths in the disk and halo, modeled as a two- or three-layer magneto-ionic medium, via differential Faraday rotation and internal Faraday dispersion, along with wavelength-independent depolarization arising from turbulent magnetic fields. A reduced chi-squared analysis is used for the statistical comparison of predicted to observed polarization maps to determine the best-fit magnetic field configuration at each of four radial rings spanning $2.4 - 7.2$ kpc in $1.2$ kpc increments. We find that a two-layer modeling approach provides a better fit to the observations than a three-layer model, where the near and far sides of the halo are taken to be identical, although the resulting best-fit magnetic field strengths are comparable. This implies that all of the signal from the far halo is depolarized at these wavelengths. We find a total magnetic field in the disk of approximately $18~\mu$G and a total magnetic field strength in the halo of $\sim 4-6~\mu$G. Both turbulent and regular magnetic field strengths in the disk exceed those in the halo by a factor of a few. About half of the turbulent magnetic field in the disk is anisotropic, but in the halo all turbulence is only isotropic.

Depolarization of diffuse radio synchrotron emission is classified in terms of wavelength-independent and wavelength-dependent depolarization in the context of regular magnetic fields and of both isotropic and anisotropic turbulent magnetic fields. Previous analytical formulas for depolarization due to differential Faraday rotation are extended to include internal Faraday dispersion concomitantly, for a multilayer synchrotron emitting and Faraday rotating magneto-ionic medium. In particular, depolarization equations for a two- and three-layer system (disk-halo, halo-disk-halo) are explicitly derived. To both serve as a `user's guide' to the theoretical machinery and as an approach for disentangling line-of-sight depolarization contributions in face-on galaxies, the analytical framework is applied to data from a small region in the face-on grand-design spiral galaxy M51. The effectiveness of the multiwavelength observations in constraining the pool of physical depolarization scenarios is illustrated for a two- and three-layer model along with a Faraday screen system for an observationally motivated magnetic field configuration.

We discuss asymptotic solutions of the kinematic $\alpha\omega$-dynamo in a thin disc (slab). Focusing upon the strong dynamo regime, in which the dynamo number $D$ satisfies $|D|\gg1$, we resolve uncertainties in the earlier treatments and conclude that some of the simplifications that have been made in previous studies are questionable. Comparing numerical solutions with asymptotic results obtained for $|D|\gg1$ and $|D|\ll1$ we find that the asymptotic solutions give a reasonably accurate description of the dynamo even far beyond their formal ranges of applicability. Indeed, our results suggest a simple analytical expression for the growth rate of the mean magnetic field that remains accurate in the range $-200< D< -10$ (which is appropriate for dynamos in spiral galaxies and accretion discs). Finally, we analyse the role of various terms in the dynamo equations to clarify the fine details of the dynamo process.

The hydrodynamic state of the interstellar medium (ISM) heated and randomly stirred by supernovae (SNe) is investigated. We use a three-dimensional non-ideal hydrodynamic ISM model in a domain extending 0.5 x 0.5 kpc horizontally and 2 kpc vertically to explore the relative importance of various physical and numerical effects on the multi-phase, turbulent ISM. We include both Type I and II SNe, the latter occurring only in dense regions. First we investigate the role of the thermal instability in the temperature range 300-6100 K, comparing results obtained for two different cooling functions, one susceptible to the instability, the other stable. The presence of thermal instability in the system is mainly visible as the tendency of the gas to avoid the relevant temperature range, as it quickly evolves towards either colder or warmer phases. Nevertheless, the formation of dense structures for both cooling functions appears to be dominated by expanding and colliding supernova remnants, rather than by the thermal instability. As we need to include a finite thermal conduction coefficient to resolve the thermal instability with our uniform grid, we also explore the effects of changing Prandtl number on the system. The purely divergent SN forcing is found to produce significant amounts of vorticity. The relative vorticity is around 0.6-0.7 for the highest Prandtl numbers explored, Pr=40, and is observed to diminish almost by a factor of two for the lowest Prandtl numbers studied, Pr=1. Rotation laws with angular velocity decreasing or increasing outwards are investigated, enabling us to separate the contributions to kinetic helicity due to rotation and shear. When angular velocity decreases outwards, these two contributions partly cancel each other, resulting in a smaller net helicity.

We derive the magnitude of fluctuations in total synchrotron intensity in the Milky Way and M33, from both observations and theory under various assumption about the relation between cosmic rays and interstellar magnetic fields. Given the relative magnitude of the fluctuations in the Galactic magnetic field (the ratio of the rms fluctuations to the mean magnetic field strength) suggested by Faraday rotation and synchrotron polarization, the observations are inconsistent with local energy equipartition between cosmic rays and magnetic fields. Our analysis of relative synchrotron intensity fluctuations indicates that the distribution of cosmic rays is nearly uniform at the scales of the order of and exceeding $100\p$, in contrast to strong fluctuations in the interstellar magnetic field at those scales. A conservative upper limit on the ratio of the the fluctuation magnitude in the cosmic ray number density to its mean value is 0.2--0.4 at scales of order 100\,pc. Our results are consistent with a mild anticorrelation between cosmic-ray and magnetic energy densities at these scales, in both the Milky Way and M33. Energy equipartition between cosmic rays and magnetic fields may still hold, but at scales exceeding 1\,kpc. Therefore, we suggest that equipartition estimates be applied to the observed synchrotron intensity smoothed to a linear scale of kiloparsec order (in spiral galaxies) to obtain the cosmic ray distribution and a large-scale magnetic field. Then the resulting cosmic ray distribution can be used to derive the fluctuating magnetic field strength from the data at the original resolution. The resulting random magnetic field is likely to be significantly stronger than existing estimates.

The exploration of mean-field galactic dynamos affected by a galactic spiral pattern, begun in Chamandy et al. (2013, hereafter Paper I) with numerical simulations, is continued here with an asymptotic solution. The mean-field dynamo model used generalizes the standard theory to include the delayed response of the mean electromotive force to variations of the mean magnetic field and turbulence parameters (the temporal non-locality, or {\tau} effect). The effect of the spiral pattern on the dynamo considered is the enhancement of the {\alpha}-effect in spiral-shaped regions (which may overlap the gaseous spiral arms or be located in the interarm regions). The axisymmetric and enslaved non-axisymmetric modes of the mean magnetic field are studied semi-analytically to clarify and strengthen the numerical results. Good qualitative agreement is obtained between the asymptotic solution and numerical solutions of Paper I for a global, rigidly rotating material spiral (density wave). At all galactocentric distances except for the co-rotation radius, we find magnetic arms displaced in azimuth from the {\alpha}-arms, so that the ridges of magnetic field strength are more tightly wound than the {\alpha}-arms. Moreover, the effect of a finite dynamo relaxation time {\tau} (related to the turbulence correlation time) is to phase-shift the magnetic arms in the direction opposite to the galactic rotation even at the co-rotation radius. This mechanism can be used to explain the phase shifts between magnetic and material arms observed in some spiral galaxies.

We simulate the multi-phase interstellar medium randomly heated and stirred by supernovae, with gravity, differential rotation and other parameters of the solar neighbourhood. Here we describe in detail both numerical and physical aspects of the model, including injection of thermal and kinetic energy by SN explosions, radiative cooling, photoelectric heating and various transport processes. With 3D domain extending 1 kpc^2 horizontally and 2 kpc vertically, the model routinely spans gas number densities 10^-5 - 10^2 cm^-3, temperatures 10-10^8 K, local velocities up to 10^3 km s^-1 (with Mach number up to 25). The thermal structure of the modelled ISM is classified by inspection of the joint probability density of the gas number density and temperature. We confirm that most of the complexity can be captured in terms of just three phases, separated by temperature borderlines at about 10^3 K and 5x10^5 K. The probability distribution of gas density within each phase is approximately lognormal. We clarify the connection between the fractional volume of a phase and its various proxies, and derive an exact relation between the fractional volume and the filling factors defined in terms of the volume and probabilistic averages. These results are discussed in both observational and computational contexts. The correlation scale of the random flows is calculated from the velocity autocorrelation function; it is of order 100 pc and tends to grow with distance from the mid-plane. We use two distinct parameterizations of radiative cooling to show that the multi-phase structure of the gas is robust, as it does not depend significantly on this choice.

The origin and structure of the magnetic fields in the interstellar medium of spiral galaxies is investigated with 3D, non-ideal, compressible MHD simulations, including stratification in the galactic gravity field, differential rotation and radiative cooling. A rectangular domain, 1x1x2 kpc^{3} in size, spans both sides of the galactic mid-plane. Supernova explosions drive transonic turbulence. A seed magnetic field grows exponentially to reach a statistically steady state within 1.6 Gyr. Following Germano (1992) we use volume averaging with a Gaussian kernel to separate magnetic field into a mean field and fluctuations. Such averaging does not satisfy all Reynolds rules, yet allows a formulation of mean-field theory. The mean field thus obtained varies in both space and time. Growth rates differ for the mean-field and fluctuating field and there is clear scale separation between the two elements, whose integral scales are about 0.7 kpc and 0.3 kpc, respectively.

We generalise the theory of mean-field galactic dynamos by allowing for temporal non-locality in the mean electromotive force (emf). This arises in random flows due to a finite response time of the mean emf to changes in the mean magnetic field, and leads to the telegraph equation for the mean field. The resulting dynamo model also includes the nonlinear dynamo effects arising from magnetic helicity balance. Within this framework, coherent large-scale magnetic spiral arms superimposed on the dominant axially symmetric magnetic structure are considered. A non-axisymmetric forcing of the mean-field dynamo by a spiral pattern (either stationary or transient) is invoked, with the aim of explaining the phenomenon of magnetic arms. For a stationary dynamo forcing by a rigidly rotating material spiral, we find corotating non-axisymmetric magnetic modes enslaved to the axisymmetric modes and strongly peaked around the corotation radius. For a forcing by transient material arms wound up by the galactic differential rotation, the magnetic spiral is able to adjust to the winding so that it resembles the material spiral at all times. There are profound effects associated with the temporal non-locality, i.e. finite `dynamo relaxation time'. For the case of a rigidly rotating spiral, a finite relaxation time causes each magnetic arm to mostly lag the corresponding material arm with respect to the rotation. For a transient material spiral that winds up, the finite dynamo relaxation time leads to a large, negative (in the sense of the rotation) phase shift between the magnetic and material arms, similar to that observed in NGC 6946 and other galaxies. We confirm that sufficiently strong random seed fields can lead to global reversals of the regular field along the radius whose long-term survival depends on specific features of a given galaxy.

We present a study of the Milky Way disk and halo magnetic field, determined from observations of Faraday rotation measure (RM) towards 641 polarized extragalactic radio sources in the Galactic longitude range 100-117 degs, within 30 degs of the Galactic plane. For |b| < 15 degs, we observe a symmetric RM distribution about the Galactic plane. This is consistent with a disk field in the Perseus arm of even parity across the Galactic mid-plane. In the range 15<|b|<30 degs, we find median rotation measures of -15+/-4 rad/m^2 and -62+/-5 rad/m^2 in the northern and southern Galactic hemispheres, respectively. If the RM distribution is a signature of the large-scale field parallel to the Galactic plane, this suggests that the halo magnetic field toward the outer Galaxy does not reverse direction across the mid-plane. The variation of RM as a function of Galactic latitude in this longitude range is such that RMs become more negative at larger |b|. This is consistent with an azimuthal magnetic field of strength 2 microGauss (7 microGauss) at a height 0.8-2 kpc above (below) the Galactic plane between the local and the Perseus spiral arm. We propose that the Milky Way could possess spiral-like halo magnetic fields similar to those observed in M51.

We have determined 194 Faraday rotation measures (RMs) of polarized extragalactic radio sources using new, multi-channel polarization observations at frequencies around 1.4 GHz from the Very Large Array (VLA) in the Galactic plane at $17^\circ \leq l \leq 63^\circ$ and $205^\circ \leq l \leq 253^\circ$. This catalog fills in gaps in the RM coverage of the Galactic plane between the Canadian Galactic Plane Survey and Southern Galactic Plane Survey. Using this catalog we have tested the validity of recently-proposed axisymmetric and bisymmetric models of the large-scale (or regular) Galactic magnetic field, and found that of the existing models we tested, an axisymmetric spiral model with reversals occurring in rings (as opposed to along spiral arms) best matched our observations. Building on this, we have performed our own modeling, using RMs from both extragalactic sources and pulsars. By developing independent models for the magnetic field in the outer and inner Galaxy, we conclude that in the inner Galaxy, the magnetic field closely follows the spiral arms, while in the outer Galaxy, the field is consistent with being purely azimuthal.Furthermore, the models contain no reversals in the outer Galaxy, and together seem to suggest the existence of a single reversed region that spirals out from the Galactic center.

(Abridged) We use new multi-wavelength radio observations, made with the VLA and Effelsberg telescopes, to study the magnetic field of the nearby galaxy M51 on scales from $200\pc$ to several $\kpc$. Interferometric and single dish data are combined to obtain new maps at \wwav{3}{6} in total and polarized emission, and earlier \wav{20} data are re-reduced. We compare the spatial distribution of the radio emission with observations of the neutral gas, derive radio spectral index and Faraday depolarization maps, and model the large-scale variation in Faraday rotation in order to deduce the structure of the regular magnetic field. We find that the \wav{20} emission from the disc is severely depolarized and that a dominating fraction of the observed polarized emission at \wav{6} must be due to anisotropic small-scale magnetic fields. Taking this into account, we derive two components for the regular magnetic field in this galaxy: the disc is dominated by a combination of azimuthal modes, $m=0+2$, but in the halo only an $m=1$ mode is required to fit the observations. We disuss how the observed arm-interarm contrast in radio intensities can be reconciled with evidence for strong gas compression in the spiral shocks. The average arm--interam contrast, representative of the radii $r>2\kpc$ where the spiral arms are broader, is not compatible with straightforward compression: lower arm--interarm contrasts than expected may be due to resolution effects and \emph{decompression} of the magnetic field as it leaves the arms. We suggest a simple method to estimate the turbulent scale in the magneto-ionic medium from the dependence of the standard deviation of the observed Faraday rotation measure on resolution. We thus obtain an estimate of $50\pc$ for the size of the turbulent eddies.

We present a study of the vertical magnetic field of the Milky Way towards the Galactic poles, determined from observations of Faraday rotation toward more than 1000 polarized extragalactic radio sources at Galactic latitudes |b| > 77 degs, using the Westerbork Radio Synthesis Telescope and the Australia Telescope Compact Array. We find median rotation measures (RMs) of 0.0 +/- 0.5 rad/m^2 and +6.3 +/- 0.7 rad/m^2 toward the north and south Galactic poles, respectively, demonstrating that there is no coherent vertical magnetic field in the Milky Way at the Sun's position. If this is a global property of the Milky Way's magnetism, then the lack of symmetry across the disk rules out pure dipole or quadrupole geometries for the Galactic magnetic field. The angular fluctuations in RM seen in our data show no preferred scale within the range ~ 0.1 to 25 degs. The observed standard deviation in RM of ~ 9 rad/m^2 then implies an upper limit of ~1microGauss on the strength of the random magnetic field in the warm ionized medium at high Galactic latitudes.

We develop a new model of the fluctuation dynamo in which the magnetic field is confined to thin flux ropes advected by a multi-scale flow which models turbulence. Magnetic dissipation occurs only via reconnections of flux ropes. The model is particularly suitable for rarefied plasma, such as the Solar corona or galactic halos. We investigate the kinetic energy release into heat, mediated by dynamo action, both in our model and by solving the induction equation with the same flow. We find that the flux rope dynamo is more than an order of magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic energy released during reconnections has a power-law form with the slope -3, consistent with the Solar corona heating by nanoflares. We also present a nonlinear extension of the model. This shows that a plausible saturation mechanism of the fluctuation dynamo is the suppression of turbulent magnetic diffusivity, due to suppression of random stretching at the location of the flux ropes. We confirm that the probability distribution function of the magnetic line curvature has a power-law form suggested by (Schekochihin et al., 2002). We argue, however, using our results that this does not imply a persistent folded structure of magnetic field, at least in the nonlinear stage.

We develop a new model of the fluctuation dynamo in which the magnetic field is confined in thin flux ropes advected by a multi-scale flow modeling turbulence. Magnetic dissipation occurs only via reconnection of the flux ropes. We investigate the kinetic energy release into heat, mediated by the dynamo action, both in our model and by solving the induction equation with the same flow. We find that a flux rope dynamo is an order of magnitude more efficient at converting mechanical energy into heat. The probability density of the magnetic energy release in reconnections has a power-law form with the slope -3, consistent with the Solar corona heating by nanoflares.

Observations show that magnetic fields in the interstellar medium (ISM) often do not respond to increases in gas density as would be naively expected for a frozen-in field. This may suggest that the magnetic field in the diffuse gas becomes detached from dense clouds as they form. We have investigated this possibility using theoretical estimates, a simple magneto-hydrodynamic model of a flow without mass conservation and numerical simulations of a thermally unstable flow. Our results show that significant magnetic flux can be shed from dense clouds as they form in the diffuse ISM, leaving behind a magnetically dominated diffuse gas.

Fluctuation dynamo action has been obtained for a flow previously used to model fluid turbulence, where the sensitivity of the magnetic field parameters to the kinetic energy spectrum can be explored. We apply quantitative morphology diagnostics, based on the Minkowski functionals, to magnetic fields produced by the kinematic fluctuation dynamo to show that magnetic structures are predominantly filamentary rather than sheet-like. Our results suggest that the thickness, width and length of the structures scale differently with magnetic Reynolds number.

We present a simple semi-analytical model of nonlinear, mean-field galactic dynamos and use it to study the effects of various magnetic helicity fluxes. The dynamo equations are reduced using the `no-$z$' approximation to a nonlinear system of ordinary differential equations in time; we demonstrate that the model reproduces accurately earlier results, including those where nonlinear behaviour is driven by a magnetic helicity flux. We discuss the implications and interplay of two types of magnetic helicity flux, one produced by advection (e.g., due to the galactic fountain or wind) and the other, arising from anisotropy of turbulence as suggested by Vishniac & Cho(2001). We argue that the latter is significant if the galactic differential rotation is strong enough: in our model, for $\Rw\la-10$ in terms of the corresponding turbulent magnetic Reynolds number. We confirm that the intensity of gas outflow from the galactic disc optimal for the dynamo action is close to that expected for normal spiral galaxies. The steady-state strength of the large-scale magnetic field supported by the helicity advection is still weaker than that corresponding to equipartition with the turbulent energy. However, the Vishniac-Cho helicity flux can boost magnetic field further to achieve energy equipartition with turbulence. For stronger outflows that may occur in starburst galaxies, the Vishniac-Cho flux can be essential for the dynamo action. However, this mechanism requires a large-scale magnetic field of at least $\simeq1\mkG$ to be launched, so that it has to be preceded by a conventional dynamo assisted by the advection of magnetic helicity by the fountain or wind.

We present a model of the global magnetic field in the barred galaxy NGC 1365 based jointly on the large-scale velocity field of interstellar gas fitted to HI and CO observations of this galaxy and on mean-field dynamo theory. The aim of the paper is to present a detailed quantitative comparison of a galactic dynamo model with independent radio observations. We consider several gas dynamical and nonlinear dynamo models that include plausible variations of parameters that are poorly known. Models of cosmic ray distribution in the galaxy are introduced to produce synthetic radio polarization maps allowing direct comparison with those observed at 3.5cm and 6.2cm. We show that the dynamo model is robust in that the most important magnetic features are controlled by the relatively well established properties of the density distribution and gas velocity field. The optimal agreement between the synthetic polarization maps and observations is obtained when a uniform cosmic ray distribution is adopted. We find some indirect evidence for enhanced turbulence in the regions of strong velocity shear and within 1-2kpc of the galactic centre. We confirm that magnetic stresses can drive an inflow of gas into the inner 1kpc of the galaxy at a rate of a few solar masses per year. The dynamo models are successful to some extent in modelling the large scale regular magnetic field in this galaxy. Our results demonstrate that dynamo models and synthetic polarization maps can provide information about both the gas dynamical models and conditions in the interstellar medium. We demonstrate that the dynamical effects of magnetic fields cannot be everywhere ignored in galaxy modelling.

The macroscopic behaviour of cosmic rays in turbulent magnetic fields is discussed. An implementation of anisotropic diffusion of cosmic rays with respect to the magnetic field in a non-conservative, high-order, finite-difference magnetohydrodynamic code is discussed. It is shown that the standard implementation fails near singular X-points of the magnetic field, which are common if the field is random. A modification to the diffusion model for cosmic rays is described and the resulting telegraph equation (implemented by solving a dynamic equation for the diffusive flux of cosmic rays) is used; it is argued that this modification may better describe the physics of cosmic ray diffusion. The present model reproduces several processes important for the propagation and local confinement of cosmic rays, including spreading perpendicular to the local large-scale magnetic field, controlled by the random-to-total magnetic field ratio, and the balance between cosmic ray pressure and magnetic tension. Cosmic ray diffusion is discussed in the context of a random magnetic field produced by turbulent dynamo action. It is argued that energy equipartition between cosmic rays and other constituents of the interstellar medium do not necessarily imply that cosmic rays play a significant role in the balance of forces.

Narrow depolarized canals are common in maps of the polarized synchrotron emission of the Milky Way. Two physical effects that can produce these canals have been identified: the presence of Faraday rotation measure ($\RM$) gradients in a foreground screen and the cumulative cancellation of polarization known as differential Faraday rotation. We show that the behaviour of the Stokes parameters $Q$ and $U$ in the vicinity of a canal can be used to identify its origin. In the case of canals produced by a Faraday screen we demonstrate that, if the polarization angle changes by $90\degr$ across the canal, as is observed in all fields to-date, the gradients in $\RM$ must be discontinuous. Shocks are an obvious source of such discontinuities and we derive a relation of the expected mean separation of canals to the abundance and Mach number of supernova driven shocks, and compare this with recent observations by \citet{Haverkorn03}. We also predict the existence of less common canals with polarization angle changes other than $90\degr$. Differential Faraday rotation can produce canals in a uniform magneto-ionic medium, but as the emitting layer becomes less uniform the canals will disappear. We show that for moderate differences in emissivity in a two-layer medium, of up to 1/2, and for Faraday depth fluctuations of standard deviation $\lesssim 1 \mathrm{rad}$, canals produced by differential rotation will still be visible.

Recent radio polarization observations have revealed a plethora of unexpected features in the polarized Galactic radio background that arise from propagation effects in the random (turbulent) interstellar medium. The canals are especially striking among them, a random network of very dark, narrow regions clearly visible in many directions against a bright polarized Galactic synchrotron background. There are no obvious physical structures in the ISM that may have caused the canals, and so they have been called Faraday ghosts. They evidently carry information about interstellar turbulence but only now is it becoming clear how this information can be extracted. Two theories for the origin of the canals have been proposed; both attribute the canals to Faraday rotation, but one invokes strong gradients in Faraday rotation in the sky plane (specifically, in a foreground Faraday screen) and the other only relies on line-of-sight effects (differential Faraday rotation). In this review we discuss the physical nature of the canals and how they can be used to explore statistical properties of interstellar turbulence. This opens studies of magnetized interstellar turbulence to new methods of analysis, such as contour statistics and related techniques of computational geometry and topology. In particular, we can hope to measure such elusive quantities as the Taylor microscale and the effective magnetic Reynolds number of interstellar MHD turbulence.

Random motions can occur in the intergalactic gas of galaxy clusters at all stages of their evolution. Depending on the poorly known value of the Reynolds number, these motions can or cannot become turbulent, but in any case they can generate random magnetic fields via dynamo action. We argue that magnetic fields inferred observationally for the intracluster medium require dynamo action, and then estimate parameters of random flows and magnetic fields at various stages of the cluster evolution. Polarization in cluster radio halos predicted by the model would be detectable with the SKA.

Nonlinear behaviour of galactic dynamos is studied, allowing for magnetic helicity removal by the galactic fountain flow. A suitable advection speed is estimated, and a one-dimensional mean-field dynamo model with dynamic alpha-effect is explored. It is shown that the galactic fountain flow is efficient in removing magnetic helicity from galactic discs. This alleviates the constraint on the galactic mean-field dynamo resulting from magnetic helicity conservation and thereby allows the mean magnetic field to saturate at a strength comparable to equipartition with the turbulent kinetic energy.

We discuss, using simple analytical models and MHD simulations, the origin and parameters of turbulence and magnetic fields in galaxy clusters. Three physically distinct regimes can be identified in the evolution of cluster turbulence and magnetic fields. Firstly, the fluctuation dynamo will produce microgauss-strong, random magnetic fields during cluster formation and major mergers. Turbulent velocity of about 300 km/s can be maintained at scales 100-200 kpc. The magnetic field is intermittent, has a smaller scale of 20-30 kpc and average strength of 2 microgauss. Secondly, when major mergers end, turbulent speed and magnetic field undergo a power-law decay, decreasing in strength but increasing in scale by a factor of about two. Thirdly, smaller-mass subclusters and cluster galaxies produce turbulent wakes, with turbulent speeds and magnetic field strengths similar to those quoted above. The velocity scales are about 200 kpc and 10 kpc respectively, and the magnetic field scale is about 6 times smaller. Although these wakes may fill only a small fraction of the cluster volume, their area covering factor can be close to unity. So one can potentially reconcile observations that indicate the coexistence of turbulence with ordered filamentary gas structures, as in the Perseus cluster. Random Faraday rotation measure is estimated to be typically 100-200 rad/m^2, in agreement with observations. We predict detectable synchrotron polarization from cluster radio halos at wavelengths 3-6 cm, if observed at sufficiently high resolution (abridged).

We present 3.5cm and 6.2cm radio continuum maps in total and polarized intensity of the barred galaxies NGC 1097 and NGC 1365. Both galaxies exhibit radio ridges roughly overlapping with the massive dust lanes in the bar region. The contrast in total intensity across the radio ridges is compatible with compression and shear of an isotropic random magnetic field. The contrast in polarized intensity is significantly smaller than that expected from compression and shearing of the regular magnetic field; this could be the result of decoupling of the regular field from the dense molecular clouds. The regular field in the ridge is probably strong enough to reduce significantly shear in the diffuse gas (to which it is coupled) and hence to reduce magnetic field amplification by shearing. This contributes to the misalignment of the observed field orientation with respect to the velocity vectors of the dense gas. Our observations, for the first time, indicate that magnetic forces can control the flow of the diffuse interstellar gas at kiloparsec scales. The total radio intensity reaches its maximum in the circumnuclear starburst regions, where the equipartition field strength is about 60\mu G, amongst the strongest fields detected in spiral galaxies so far. The regular field in the inner region has a spiral shape with large pitch angle, indicating the action of a dynamo. Magnetic stress leads to mass inflow towards the centre, sufficient to feed the active nucleus in NGC 1097. We detected diffuse X-ray emission, possibly forming a halo of hot gas around NGC 1097.

A review of galactic magnetic fields is given with emphasis on simple analytical solutions of relevant equations and detailed comparison with observations. It is shown that mean-field dynamo theory, in its simplest form, is in a remarkably good agreement with observations of galactic magnetic fields. It is argued that the disc-halo connection can play crucial role in the nonlinear evolution of regular magnetic fields in spiral galaxies by driving flux of magnetic helicity from the disc, thus supporting healthy dynamo action.

The thin-disc global asymptotics are discussed for axisymmetric mean-field dynamos with vacuum boundary conditions allowing for non-local terms arising from a finite radial component of the mean magnetic field at the disc surface. This leads to an integro-differential operator in the equation for the radial distribution of the mean magnetic field strength, $Q(r)$ in the disc plane at a distance $r$ from its centre; an asymptotic form of its solution at large distances from the dynamo active region is obtained. Numerical solutions of the integro-differential equation confirm that the non-local effects act similarly to an enhanced magnetic diffusion. This leads to a wider radial distribution of the eigensolution and faster propagation of magnetic fronts, compared to solutions with the radial surface field neglected. Another result of non-local effects is a slowly decaying algebraic tail of the eigenfunctions outside the dynamo active region, $Q(r)\sim r^{-4}$, which is shown to persist in nonlinear solutions where $\alpha$-quenching is included. The non-local nature of the solutions can affect the radial profile of the regular magnetic field in spiral galaxies and accretion discs at large distances from the centre.

Numerical simulations of the multi-phase interstellar medium have been carried out, using a 3D, nonlinear, magnetohydrodynamic, shearing-box model, with random motions driven by supernova explosions. These calculations incorporate the effects of magnetic fields and rotation in 3D; these play important dynamical roles in the galaxy, but are neglected in many other simulations. The supernovae driving the motions are not arbitrarily imposed, but occur where gas accumulates into cold, dense clouds; their implementation uses a physically motivated model for the evolution of such clouds. The process is self-regulating, and produces mean supernova rates as part of the solution. Simulations with differing mean density show a power law relation between the supernova rate and density, with exponent 1.7; this value is within the range suggested from observations (taking star formation rate as a proxy for supernova rate). The global structure of the supernova driven medium is strongly affected by the presence of magnetic fields; e.g. for one solution the filling factor of hot gas is found to vary from 0.19 (with no field) to 0.12 (with initial mid-plane field B_0=6 muG).

The configuration of the regular magnetic field in M31 is deduced from radio polarization observations at the wavelengths 6, 11 and 20 cm. By fitting the observed azimuthal distribution of polarization angles, we find that the regular magnetic field, averaged over scales 1--3 kpc, is almost perfectly axisymmetric in the radial range 8 to 14 kpc, and follows a spiral pattern with pitch angles of $p\simeq -19\degr$ to $p\simeq -8\degr$. In the ring between 6 and 8 kpc a perturbation of the dominant axisymmetric mode may be present, having the azimuthal wave number m=2. A systematic analysis of the observed depolarization allows us to identify the main mechanism for wavelength dependent depolarization -- Faraday rotation measure gradients arising in a magneto-ionic screen above the synchrotron disk. Modelling of the depolarization leads to constraints on the relative scale heights of the thermal and synchrotron emitting layers in M31; the thermal layer is found to be up to three times thicker than the synchrotron disk. The regular magnetic field must be coherent over a vertical scale at least similar to the scale height of the thermal layer, estimated to be $h\therm\simeq 1$ kpc. Faraday effects offer a powerful method to detect thick magneto-ionic disks or halos around spiral galaxies.

We reconsider thin-disc global asymptotics for kinematic, axisymmetric mean-field dynamos with vacuum boundary conditions. Non-local terms arising from a small but finite radial field component at the disc surface are consistently taken into account for quadrupole modes. As in earlier approaches, the solution splits into a local part describing the field distribution along the vertical direction and a radial part describing the radial (global) variation of the eigenfunction. However, the radial part of the eigenfunction is now governed by an integro-differential equation whose kernel has a weak (logarithmic) singularity. The integral term arises from non-local interactions of magnetic fields at different radii through vacuum outside the disc. The non-local interaction can have a stronger effect on the solution than the local radial diffusion in a thin disc, however the effect of the integral term is still qualitatively similar to magnetic diffusion.

Faraday rotation of the polarization plane in magnetized thermal plasma provides one of the most efficient methods to deduce regular magnetic fields from radio astronomical observations. Since the Faraday rotation measure RM is proportional to an integral, along the line of sight, of magnetic field weighted with thermal electron density, RM is believed to yield the regular magnetic field averaged over large volume. Here we show that this is not the case in a turbulent medium where fluctuations in magnetic field and electron density are not statistically independent, and so contribute to RM. For example, in the case of pressure equilibrium, magnetic field can be anticorrelated with plasma density to produce a negative contribution. As a result, the strength of the regular magnetic field obtained from RM can be underestimated if the fluctuations in electron density and magnetic field are neglected. The anticorrelation also reduces the standard deviation of RM. We further discuss the effect of the positive correlations where the standard treatment of RM leads to an overestimated magnetic field. Because of the anisotropy of the turbulent magnetic field, the regular magnetic fields strength, obtained from synchrotron emission using standard formulae, can be overestimated. A positive correlation between cosmic-ray number density and magnetic field leads to an overestimate of the strengths of the regular and total fields. These effects can explain the difference between the strengths of the regular Galactic magnetic field as indicated by RM and synchrotron emissivity data and reconcile the magnetic field strength in the Solar vicinity with typical strength of regular magnetic fields in external galaxies.

An axisymmetric model of a cool, dynamo-active accretion disc is applied to protostellar discs. Thermally and magnetically driven outflows develop that are not collimated within 0.1 AU. In the presence of a central magnetic field from the protostar, accretion onto the protostar is highly episodic, which is in agreement with earlier work.

Narrow, elongated regions of very low polarized intensity -- so-called canals -- have recently been observed by several authors at decimeter wavelengths in various directions in the Milky Way, but their origin remains enigmatic. We show that the canals arise from depolarization by differential Faraday rotation in the interstellar medium and that they represent level lines of Faraday rotation measure RM, a random function of position in the sky. Statistical properties of the separation of canals depend on the autocorrelation function of RM, and so provide a useful tool for studies of interstellar turbulence.

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

We present an axisymmetric numerical model of a dynamo active accretion disc. If the dynamo-generated magnetic field in the disc is sufficiently strong (close to equipartition with thermal energy), a fast magneto-centrifugally driven outflow develops within a conical shell near the rotation axis, together with a slower pressure driven outflow from the outer parts of the disc as well as around the axis. Our results show that a dynamo active accretion disc can contribute to driving an outflow even without any external magnetic field. The fast outflow in the conical shell is confined by the azimuthal field which is produced by the dynamo in the disc and advected to the disc corona. This part of the outflow has high angular momentum and is cooler and less dense than its surroundings. The conical shell's half-opening angle is typically about 30 degrees near the disc and decreases slightly with height. The slow outflow is hotter and denser.

The total and polarized radio continuum emission of 20 barred galaxies was observed with the Very Large Array (VLA) at 3, 6, 18 and 22 cm and with the Australia Telescope Compact Array (ATCA) at 6 cm and 13 cm. Maps at 30 arcsec angular resolution are presented here. Polarized emission (and therefore a large-scale regular magnetic field) was detected in 17 galaxies. Most galaxies of our sample are similar to non-barred galaxies with respect to the radio/far-infrared flux correlation and equipartition strength of the total magnetic field. Galaxies with highly elongated bars are not always radio-bright. We discuss the correlation of radio properties with the aspect ratio of the bar and other measures of the bar strength. We introduce a new measure of the bar strength, \Lambda, related to the quadrupole moment of the bar's gravitational potential. The radio surface brightness I of the barred galaxies in our sample is correlated with \Lambda, I \propto \Lambda^0.4+/-0.1, and thus is highest in galaxies with a long bar where the velocity field is distorted by the bar over a large fraction of the disc. In these galaxies, the pattern of the regular field is significantly different from that in non-barred galaxies. In particular, field enhancements occur upstream of the dust lanes where the field lines are oriented at large angles to the bar's major axis. Polarized radio emission seems to be a good indicator of large-scale non-axisymmetric motions.

The self-excitation of magnetic field by a spiral Couette flow between two coaxial cylinders is considered. We solve numerically the fully nonlinear, three-dimensional MHD equations for magnetic Prandtl numbers Pm (ratio of kinematic viscosity to magnetic diffusivity) between 0.14 and 10 and kinematic and magnetic Reynolds numbers up to about 2000. In the initial stage of exponential field growth (kinematic dynamo regime), we find that the dynamo switches from one distinct regime to another as the radial width \delta r_B$ of the magnetic field distribution becomes smaller than the separation of the field maximum from the flow boundary. The saturation of magnetic field growth is due to a reduction in the velocity shear resulting mainly from the azimuthally and axially averaged part of the Lorentz force, which agrees with an asymptotic result for the limit of Pm<<1. In the parameter regime considered, the magnetic energy decreases with kinematic Reynolds number as Re^{-0.84}, which is approximately as predicted by the nonlinear asymptotic theory (\sim Re^{-1}). However, when the velocity field is maintained by a volume force (rather than by viscous stress) the dependence of magnetic energy on the kinematic Reynolds number is much weaker.

We suggest a two-dimensional wavelet devised to deduce the large-scale structure of a physical field (e.g., the Galactic magnetic field) from its integrals along straight paths from irregularly spaced data points to a fixed interior point (the observer). The method can be applied to the analysis of pulsar rotation and dispersion measures in terms of the large-scale Galactic magnetic field and electron density. The method does not use any 'a priori' assumptions about the physical field and can be considered as an algorithm of wavelet differentiation. We argue that a certain combination of the wavelet transformation with model fitting would be most efficient in the interpretation of the available pulsar RM data.

We study the generation and maintenance of large-scale magnetic fields in barred galaxies. We take a velocity field (with strong noncircular components) from a published gas dynamical simulation of Athanassoula (1992), and use this as input to a galactic dynamo calculation. Our work is largely motivated by recent high quality VLA radio observations of the barred galaxy NGC 1097, and we compare our results in detail with the regular magnetic fields deduced from these observations. We are able to reproduce most of the conspicuous large-scale features of the observed regular field, including the field structure in the central regions, by using a simple mean-field dynamo model in which the intensity of interstellar turbulence (more precisely, the turbulent diffusivity) is enhanced by a factor of 2-6 in the dust lanes and near the circumnuclear ring. We argue that magnetic fields can be dynamically important, and therefore should be included in models of gas flow in barred galaxies.

We examine the behaviour of an axisymmetric galactic dynamo model with a radial accretion flow in the disc. We also introduce a vertical magnetic flux through the galactic midplane, to simulate the presence of a large scale magnetic field trapped by the galaxy when forming. The trapped vertical flux is conserved and advected towards the disc centre by the radial flow. We confirm that accretion flows of magnitude several km/s through a significant part of the galactic disc can markedly inhibit dynamo action. Moreover, advection of the vertical flux in general results in mixed parity galactic fields. However, the effect is nonlinear and non-additive -- global magnetic field energies are usually significantly smaller that the sum of purely dynamo generated and purely advected field energies. For large inflow speeds, a form of `semi-dynamo' action may occur. We apply our results to the accumulation and redistribution, by a radial inflow, of a vertical magnetic flux captured by the Galactic disc. Taking representative values, it appears difficult to obtain mean vertical fields near the centre of the Milky Way that are much in excess of 10 microgauss, largely because the galactic dynamo and turbulent magnetic diffusion modify the external magnetic field before it can reach the disc centre.

Coherent structures in the distribution of the Faraday rotation measure of extragalactic radio sources are isolated using wavelet transformation technique. Signatures of the magnetic fields in the local (Orion) spiral arm, the Sagittarius and Carina arms, the synchrotron Loop I and, possibly, the Perseus arm have been revealed using the RM catalogues of Simard-Normandin et al. (1981, 551 source)and Broten et al. (1988, 663 sources). Unlike earlier analyses of the RM sky, our approach has allowed us to assess the stability of the results with respect to modifications of the data sample. We estimate the mean magnetic field strength in the Orion arm as 1.4\pm0.3 mkGauss and find clear indication of a magnetic field reversal at a distance of 0.6-1 kpc towards the Galactic centre. Our analysis has revealed for the first time the extension of the reversal in the Carina arm. The direction of magnetic field in the Perseus arm is reversed with respect to that in the Orion arm. The average pitch angle of magnetic field in the nearby spiral arms is 15^\circ, and the mean field strength in the Sagittarius-Carina and Perseus arms is 1.7\pm0.3 mkGauss and 1.4\pm1.2 mkGauss, respectively. The line-of-sight magnetic field in Loop I is estimated as 0.9\pm0.3 mkGauss. We find firm evidence of a dominant even symmetry of the local mean field with respect to the galactic equator.

We consider the effect of vertical outflows on the mean-field dynamo in a thin disk. These outflows could be due to winds or magnetic buoyancy. We analyse both two-dimensional finite-difference numerical solutions of the axisymmetric dynamo equations and a free-decay mode expansion using the thin-disk approximation. Contrary to expectations, a vertical velocity can enhance dynamo action, provided it is not too strong. In the nonlinear regime this can lead to super-exponential growth of the magnetic field.

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 energy equipartition is maintained between turbulence and random magnetic fields. Our results indicate that a midplane value of $B_{0}=4\mkG$ for the regular magnetic field near the Sun leads to more attractive models than $B_{0}=2\mkG$. The vertical profiles of both the regular and random magnetic fields contain disc and halo components whose parameters we have determined. The layer at $1\la|z|\la4\kpc$ can be overpressured and an outflow at a speed of about $50\kms$ may occur there, presumably associated with a Galactic fountain flow, if $B_{0}\simeq 2\mkG$. We show that hydrostatic equilibrium in a warped disc must produce asymmetric density distributions in $z$, in rough agreement with \HI 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.

Observational evidence for dynamo action in spiral galaxies is reviewed, and the capabilities of various theories in explaining the basic features of galactic magnetic fields are discussed. Mean-field dynamo models appear to be unique in providing a coherent explanation of a wide variety of magnetic features in spiral galaxies. Reliable modelling of global magnetic structures, such as the magnetic ring in M31, requires detailed knowledge of the rotation curve, the magnitude and radial profile of turbulent and noncircular systematic velocities, the scale height of the warm ionized layer, the total gas density, the turbulent scale and their variations with galactocentric radius. More detailed models involving the effects of the spiral arms on magnetic field require the knowledge of the arm-interarm contrasts in the above quantities.

Recent numerical simulations of the interstellar medium driven by energy input from supernovae and stellar winds indicate that HI clouds can be formed by compression in shock waves and colliding turbulent streams without any help from thermal and gravitational instabilities. The filling factor of the hot phase in these models does not exceed 20-40% at the midplane. The hot gas is involved in a systematic vertical outflow at |z|<1-3 kpc, similar to that expected for galactic fountains, whereas the warm component may remain in hydrostatic equilibrium. The turbulent velocity is larger in the warmer phases, being 3, 10 and 40 km/s in the cool, warm and hot phases, respectively, according to one of the simulations. The models exhibit global variability in the total kinetic and thermal energy and star formation rate at periods of (0.4-4)x10^8 yr. Current models are still unable to reproduce dynamo action in the interstellar gas; we briefly discuss implications of the dynamo theory for turbulent interstellar magnetic fields.

We study the effects of vertical outflows on mean-field dynamos in disks. These outflows could be due to thermal winds or magnetic buoyancy. We analyse numerical solutions of the nonlinear mean-field dynamo equations using a two-dimensional finite-difference model. Contrary to expectations, a modest vertical velocity can enhance dynamo action. This can lead to super-exponential growth of the magnetic field and to higher magnetic energies at saturation in the nonlinear regime.

We discuss the regular magnetic field in the galaxy M51 on the basis of a new interpretation of polarization angles observed at the wavelengths 2.8, 6.2, 18.0 and 20.5 cm. We found a magneto-ionic halo in M51 with a radial extent of about 10 kpc. The regular magnetic fields in the disk and the halo have different structures. The regular magnetic field in the halo is axisymmetric and horizontal. Its field lines are spirals pointing inwards and generally opposite to those in the disk. The azimuthal structure of the magnetic field in the disk is neither axisymme- tric nor bisymmetric but can be represented by a superposition of these two basic harmonics with about equal weights. Magnetic lines of the regular field in the disk are spirals generally directed outwards. They are well aligned with the optical spiral arms as traced by the dust lanes. The regular magnetic field strength averaged in 3 kpc wide rings is about 5-10 muG in the disk and 3muG in the radial range 3--6 kpc in the halo. These values are in good agreement with independent estimates from the synchrotron emission. The general features of the magnetic patterns revealed in the disk and the halo seem to be in agreement with predictions of dynamo theory. Our results are also based on an analysis of the magneto-ionic medium in M51. For a range of radii in the galactic disk, we estimated the scale heights of synchrotron and thermal disks, electron densities and filling factors as well as the amount of depolarization.