The grand-design face-on spiral galaxy M51 is an excellent laboratory for studying magnetic fields in galaxies. We present new observations of M51 using the VLA at the frequency range of S-band (2-4GHz), to shed new light on the transition region between the disk and halo. We present images of the distributions of the total intensity, polarized intensity, degree of polarization, and rotation measure (RM). The RM distribution in S-band shows a fluctuating pattern without any apparent large-scale structure. We discuss a model of the depolarization of synchrotron radiation in a multi-layer magneto-ionic medium and compare the model predictions to the polarization data of M51 between 1-8GHz. Since the model predictions strongly differ within the wavelength range of the S-band, the new data are essential. The parameters of the model are adjusted to fit to the data of polarization fractions in a few selected regions. In three spiral arm regions, the turbulent field in the disk dominates with strengths between 18muG and 24muG, while the regular field strengths are 8-16muG. In one inter-arm region, the regular field strength of 18muG exceeds that of the turbulent field of 11muG. The regular field strengths in the halo are 3-5muG. The observed RMs in the disk-halo transition region are probably dominated by tangled regular fields, as predicted from models of evolving dynamos, and/or vertical fields, as predicted from numerical simulations of Parker instabilities or galactic winds. Both types of magnetic fields have frequent reversals on scales similar to or larger than the beam size (550pc) that contribute to an increase of the RM dispersion and to distortions of any large-scale pattern of the regular field. Our study devises new ways of analyzing and interpreting broadband multi-frequency polarization data that will be applicable to future data from, for example, the Square Kilometre Array.

In this paper we present a detailed analysis of the Faraday depth (FD) spectrum and its clean components obtained through the application of the commonly used technique of Faraday rotation measure synthesis to analyze spectro-polarimetric data. In order to directly compare the Faraday depth spectrum with physical properties of a magneto-ionic medium, we generated synthetic broad-bandwidth spectro-polarimetric observations from magnetohydrodynamic (MHD) simulations of a transonic, isothermal, compressible turbulent medium. We find that correlated magnetic field structures give rise to a combination of spiky, localized peaks at certain FD values, and broad structures in the FD spectrum. Although the majority of these spiky FD structures appear narrow, giving an impression of a Faraday thin medium, we show that they arise from strong synchrotron emissivity at that FD. Strong emissivity at a FD can arise because of both strong spatially-local polarized synchrotron emissivity at a FD or accumulation of weaker emissions along the distance through a medium that have Faraday depths within half the width of the rotation measure spread function. Such a complex Faraday depth spectrum is a natural consequence of MHD turbulence when the lines of sight pass through a few turbulent cells. This therefore complicates the convention of attributing narrow FD peaks to presence of a Faraday rotating medium along the line of sight. Our work shows that it is difficult to extract the FD along a line of sight from the Faraday depth spectrum using standard methods for a turbulent medium in which synchrotron emission and Faraday rotation occur simultaneously.

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.

Much data on the Galactic polarized radio emission has been gathered in the last five decades. All-sky surveys have been made, but only in narrow, widely spaced frequency bands, and the data are inadequate for the characterization of Faraday rotation, the main determinant of the appearance of the polarized radio sky at decimetre wavelengths. We describe a survey of the polarized radio emission from the Southern sky, aiming to characterize the magneto-ionic medium, particularly the strength and configuration of the magnetic field. This work is part of the Global Magneto-Ionic Medium Survey (GMIMS). We have designed and built a feed and receiver covering the band 300 to 900 MHz for the CSIRO Parkes 64-m Telescope. We have surveyed the entire sky between declinations -90 and +20 degrees. We present data covering 300 to 480 MHz with angular resolution 81' to 45'. The survey intensity scale is absolutely calibrated, based on measurements of resistors at known temperatures and on an assumed flux density and spectral index for Taurus A. Data are presented as brightness temperatures. We have applied Rotation Measure Synthesis to the data to obtain a Faraday depth cube of resolution 5.9 radians per metre squared, sensitivity of 60 mK of polarized intensity, and angular resolution 1.35 degrees. The data presented in this paper are available at the Canadian Astronomy Data Centre.

Sharpless 2-27 (Sh2-27) is a nearby ${\rm H\,{\small II}}$ region excited by $\zeta$Oph. We present observations of polarized radio emission from 300 to 480$\,$MHz towards Sh2-27, made with the Parkes 64$\,$m Radio Telescope as part of the Global Magneto-Ionic Medium Survey. These observations have an angular resolution of $1.35^{\circ}$, and the data are uniquely sensitive to magneto-ionic structure on large angular scales. We demonstrate that background polarized emission towards Sh2-27 is totally depolarized in our observations, allowing us to investigate the foreground. We analyse the results of Faraday tomography, mapping the magnetised interstellar medium along the 165$\,$pc path to Sh2-27. The Faraday dispersion function in this direction has peaks at three Faraday depths. We consider both Faraday thick and thin models for this observation, finding that the thin model is preferred. We further model this as Faraday rotation of diffuse synchrotron emission in the Local Bubble and in two foreground neutral clouds. The Local Bubble extends for 80$\,$pc in this direction, and we find a Faraday depth of $-0.8 \pm 0.4\,$rad$\,$m$^{-2}$. This indicates a field directed away from the Sun with a strength of $-2.5\pm1.2\,\mu$G. The near and far neutral clouds are each about 30$\,$pc thick, and we find Faraday depths of $-6.6\pm0.6\,$rad$\,$m$^{-2}$ and $+13.7\pm0.8\,$rad$\,$m$^{-2}$, respectively. We estimate that the line-of-sight magnetic strengths in the near and far cloud are $B_{\parallel, \text{near}} \approx -15\,\mu\text{G}$ and $B_{\parallel, \text{far}} \approx +30\,\mu\text{G}$. Our results demonstrate that Faraday tomography can be used to investigate the magneto-ionic properties of foreground features in front of nearby ${\rm H\,{\small II}}$ regions.

SPICA, the cryogenic infrared space telescope recently pre-selected for a `Phase A' concept study as one of the three remaining candidates for ESA's fifth medium class (M5) mission, is foreseen to include a far-infrared polarimetric imager (SPICA-POL, now called B-BOP), which would offer a unique opportunity to resolve major issues in our understanding of the nearby, cold magnetized Universe. This paper presents an overview of the main science drivers for B-BOP, including high dynamic range polarimetric imaging of the cold interstellar medium (ISM) in both our Milky Way and nearby galaxies. Thanks to a cooled telescope, B-BOP will deliver wide-field 100-350 micron images of linearly polarized dust emission in Stokes Q and U with a resolution, signal-to-noise ratio, and both intensity and spatial dynamic ranges comparable to those achieved by Herschel images of the cold ISM in total intensity (Stokes I). The B-BOP 200 micron images will also have a factor ~30 higher resolution than Planck polarization data. This will make B-BOP a unique tool for characterizing the statistical properties of the magnetized interstellar medium and probing the role of magnetic fields in the formation and evolution of the interstellar web of dusty molecular filaments giving birth to most stars in our Galaxy. B-BOP will also be a powerful instrument for studying the magnetism of nearby galaxies and testing galactic dynamo models, constraining the physics of dust grain alignment, informing the problem of the interaction of cosmic rays with molecular clouds, tracing magnetic fields in the inner layers of protoplanetary disks, and monitoring accretion bursts in embedded protostars.

Faraday rotation occurs along every line of sight in the Galaxy; Rotation Measure (RM) synthesis allows a three-dimensional representation of the interstellar magnetic field. This study uses data from the Global Magneto-Ionic Medium Survey, a combination of single-antenna spectro-polarimetric studies, including northern sky data from the DRAO 26-m Telescope (1270-1750 MHz) and southern sky data from the Parkes 64-m Telescope (300-480 MHz). From the synthesized Faraday spectral cubes we compute the zeroth, first, and second moments to find the total polarized emission, mean and RM-width of the polarized emission. From DRAO first moments we find a weak vertical field directed from Galactic North to South, but Parkes data reveal fields directed towards the Sun at high latitudes in both hemispheres: the two surveys clearly sample different volumes. DRAO second moments show feature widths in Faraday spectra increasing with decreasing positive latitudes, implying that longer lines of sight encounter more Faraday rotating medium, but this is not seen at negative latitudes. Parkes data show the opposite: at positive latitudes the second moment decreases with decreasing latitude, but not at negative latitudes. Comparing first moments with RMs of pulsars and extragalactic sources and a study of depolarization together confirm that the DRAO survey samples to larger distances than the Parkes data. Emission regions in the DRAO survey are typically 700 to 1000 pc away, slightly beyond the scale-height of the magneto-ionic medium; emission detected in the Parkes survey is entirely within the magneto-ionic disk, less than 500 pc away.

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.

An excellent laboratory for studying large scale magnetic fields is the grand de- sign face-on spiral galaxy M51. Due to wavelength-dependent Faraday depolarization, linearly polarized synchrotron emission at different radio frequencies gives a picture of the galaxy at different depths: Observations at L-band (1-2 GHz) probe the halo region while at C- and X- band (4-8 GHz) the linearly polarized emission probe the disk region of M51. We present new observations of M51 using the Karl G. Jansky Very Large Array (VLA) at S-band (2-4 GHz), where previously no polarization observations existed, to shed new light on the transition region between the disk and the halo. We discuss a model of the depolarization of synchrotron radiation in a multilayer magneto-ionic medium and compare the model predictions to the multi-frequency polarization data of M51 between 1-8GHz. The new S-band data are essential to distinguish between different models. Our study shows that the initial model parameters, i.e. the total reg- ular and turbulent magnetic field strengths in the disk and halo of M51, need to be adjusted to successfully fit the models to the data.

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 new generation of broad-band radio continuum surveys will provide large data sets with polarization information. New algorithms need to be developed to extract reliable catalogs of linearly polarized sources that can be used to characterize those sources and produce a dense rotation measure (RM) grid to probe magneto-ionized structures along the line of sight via Faraday rotation. The aim of the paper is to develop a computationally efficient and rigorously defined source-finding algorithm for linearly polarized sources. We used a calibrated data set from the LOw Frequency ARray (LOFAR) at 150 MHz centered on the nearby galaxy M51 to search for polarized background sources. With a new imaging software, we re-imaged the field at a resolution of 18''x15'' and cataloged a total of about 3000 continuum sources within 2.5 degrees of the center of M51. We made small Stokes Q and U images centered on each source brighter than 100 mJy in total intensity (201 sources) and used RM synthesis to create corresponding Faraday cubes that were analyzed individually. For each source, the noise distribution function was determined from a subset of measurements at high Faraday depths where no polarization is expected; the peaks in polarized intensity in the Faraday spectrum were identified and the p-value of each source was calculated. Finally, the false discovery rate method was applied to the list of p-values to produce a list of polarized sources and quantify the reliability of the detections. We also analyzed sources fainter than 100 mJy but that were reported as polarized in the literature at at least another radio frequency. Of the 201 sources that were searched for polarization, 6 polarized sources were detected confidently (with a false discovery rate of 5 percent). This corresponds to a number density of 1 polarized source per 3.3 square degrees, or 0.3 source/square degree. [Abridged]

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.

Low-mass galaxies are subject to strong galactic outflows, in which cosmic rays may play an important role, they can be best traced with low-frequency radio continuum observations, which are less affected by spectral ageing. We present a study of the nearby star burst dwarf irregular galaxy IC 10 using observations at 140 MHz with the LOw-Frequency ARray (LOFAR), at 1580 MHz with the Very Large Array (VLA) and at 6200 MHz with the VLA and the 100-m Effelsberg telescope. We find that IC 10 has a low-frequency radio halo, which manifests itself as a second component (thick disc) in the minor axis profiles of the non-thermal radio continuum emission at 140 and 1580 MHz. These profiles are then fitted with 1D cosmic-ray transport models for pure diffusion and advection. We find that a diffusion model fits best, with a diffusion coefficient of $D=(0.4$-$0.8) \times 10^{26}(E/{\rm GeV})^{0.5}~{\rm cm^2\,s^{-1}}$, which is at least an order of magnitude smaller than estimates both from anisotropic diffusion and the diffusion length. In contrast, advection models, which cannot be ruled out due to the mild inclination, while providing poorer fits, result in advection speeds close to the escape velocity of $\approx$$50~\rm km\,s^{-1}$, as expected for a cosmic-ray driven wind. Our favoured model with an accelerating wind provides a self-consistent solution, where the magnetic field is in energy equipartition with both the warm neutral and warm ionized medium with an important contribution from cosmic rays. Consequently, cosmic rays can play a vital role for the launching of galactic winds in the disc--halo interface.

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.

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.

Cosmic ray electrons (CREs) are a crucial part of the ISM and are observed via synchrotron emission. While much modelling has been carried out on the CRE distribution and propagation of the Milky Way, little has been done on normal external star-forming galaxies. Recent spectral data from a new generation of radio telescopes enable us to find more robust estimations of the CRE propagation. We model the synchrotron spectral index of M 51 using the time-dependent diffusion energy-loss equation and to compare the model results with the observed spectral index determined from recent low-frequency observations with LOFAR. This is the first time that this model for CRE propagation has been solved for a realistic distribution of CRE sources, which we derive from the observed star formation rate, in an external galaxy. The radial variation of the synchrotron spectral index and scale-length produced by the model are compared to recent LOFAR and older VLA observational data and also to new observations of M 51 at 325MHz obtained with the GMRT. We find that propagation of CREs by diffusion alone is sufficient to reproduce the observed spectral index distribution in M 51. An isotropic diffusion coefficient with a value of $6.6\pm0.2\,\times 10^{28}\mathrm{cm}^2\,\mathrm{s}^{-1}$ is found to fit best and is similar to what is seen in the Milky Way. We estimate an escape time of $11\,\mathrm{Myr}$ from the central galaxy to $88\,\mathrm{Myr}$ in the extended disk. It is found that an energy dependence of the diffusion coefficient is not important for CRE energies in the range $0.01\,\mathrm{GeV}$--$3\,\mathrm{GeV}$. We are able to reproduce the dependence of the observed synchrotron scale-lengths on frequency, with $l \propto \nu^{-1/4}$ in the outer disk and $l \propto \nu^{-1/8}$ in the inner disk.

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.

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.

We present probability distribution functions (PDFs) of the surface densities of ionized and neutral gas in the nearby spiral galaxies M31 and M51, as well as of dust emission and extinction Av in M31. The PDFs are close to lognormal and those for HI and Av in M31 are nearly identical. However, the PDFs for H2 are wider than the HI PDFs and the M51 PDFs have larger dispersions than those for M31. We use a simple model to determine how the PDFs are changed by variations in the line-of-sight (LOS) pathlength L through the gas, telescope resolution and the volume filling factor of the gas, f_v. In each of these cases the dispersion sigma of the lognormal PDF depends on the variable with a negative power law. We also derive PDFs of mean LOS volume densities of gas components in M31 and M51. Combining these with the volume density PDFs for different components of the ISM in the Milky Way (MW), we find that sigma decreases with increasing length L with an exponent of -0.76 +/- 0.06, which is steeper than expected. We show that the difference is due to variations in f_v. As f_v is similar in M31, M51 and the MW, the density structure in the gas in these galaxies must be similar. Finally, we demonstrate that an increase in f_v with increasing distance to the Galactic plane explains the decrease in sigma with latitude of the PDFs of emission measure and FUV emission observed for the MW.

Magnetic fields play an important role in shaping the structure and evolution of the interstellar medium (ISM) of galaxies, but the details of this relationship remain unclear. With SKA1, the 3D structure of galactic magnetic fields and its connection to star formation will be revealed. A highly sensitive probe of the internal structure of the magnetoionized ISM is the partial depolarization of synchrotron radiation from inside the volume. Different configurations of magnetic field and ionized gas within the resolution element of the telescope lead to frequency-dependent changes in the observed degree of polarization. The results of spectro-polarimetric observations are tied to physical structure in the ISM through comparison with detailed modeling, supplemented with the use of new analysis techniques that are being actively developed and studied within the community such as Rotation Measure Synthesis. The SKA will enable this field to come into its own and begin the study of the detailed structure of the magnetized ISM in a sample of nearby galaxies, thanks to its extraordinary wideband capabilities coupled with the combination of excellent surface brightness sensitivity and angular resolution.

Magnetic fields are an important ingredient of the interstellar medium (ISM). Besides their importance for star formation, they govern the transport of cosmic rays, relevant to the launch and regulation of galactic outflows and winds, which in turn are pivotal in shaping the structure of halo magnetic fields. Mapping the small-scale structure of interstellar magnetic fields in many nearby galaxies is crucial to understand the interaction between gas and magnetic fields, in particular how gas flows are affected. Elucidation of the magnetic role in, e.g., triggering star formation, forming and stabilising spiral arms, driving outflows, gas heating by reconnection and magnetising the intergalactic medium has the potential to revolutionise our physical picture of the ISM and galaxy evolution in general. Radio polarisation observations in the very nearest galaxies at high frequencies (>= 3 GHz) and with high spatial resolution (<= 5") hold the key here. The galaxy survey with SKA1 that we propose will also be a major step to understand the galactic dynamo, which is important for models of galaxy evolution and for astrophysical magnetohydrodynamics in general. Field amplification by turbulent gas motions, which is crucial for efficient dynamo action, has been investigated so far only in simulations, while compelling evidence of turbulent fields from observations is still lacking.

We present a study of the magneto-ionic medium in the Whirlpool galaxy (M51) using new wide-band multi-configuration polarization data at L band (1-2 GHz) obtained at the Karl G. Jansky Very Large Array. By fitting the observed diffuse complex polarization $Q$+$iU$ as a function of wavelength directly to various depolarization models, we find that polarized emission from M51 at 1-2 GHz originates from the top of the synchrotron disk and then experiences Faraday rotation in the near-side thermal halo of the galaxy. Thus, the scale height of the thermal gas must exceed that of the synchrotron emitting gas at L band. The observed Faraday depth distribution at L band is consistent with a halo field that comprises of a plane-parallel bisymmetric component and a vertical component which produces a Faraday rotation of $\sim$ $-$9 rad m$^{-2}$. The derived rotation measure structure functions indicate a characteristic scale of rotation measure fluctuations of less than 560 pc in the disk and approximately 1 kpc in the halo. The outer scale of turbulence of 1 kpc found in the halo of M51 is consistent with superbubbles and the Parker instability being the main energy injection mechanisms in galactic halos.

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.

The grand-design spiral galaxy M51 was observed with the LOFAR High Frequency Antennas (HBA) and imaged in total intensity and polarisation. This observation covered the frequencies between 115 MHz and 175 MHz. We produced an image of total emission of M51 at the mean frequency of 151 MHz with 20 arcsec resolution and 0.3 mJy rms noise, which is the most sensitive image of a galaxy at frequencies below 300 MHz so far. The integrated spectrum of total radio emission is described well by a power law, while flat spectral indices in the central region indicate thermal absorption. We observe that the disk extends out to 16 kpc and see a break in the radial profile near the optical radius of the disk. Our main results, the scale lengths of the inner and outer disks at 151 MHz and 1.4 GHz, arm--interarm contrast, and the break scales of the radio--far-infrared correlations, can be explained consistently by CRE diffusion, leading to a longer propagation length of CRE of lower energy. The distribution of CRE sources drops sharply at about 10 kpc radius, where the star formation rate also decreases sharply. We find evidence that thermal absorption is primarily caused by HII regions. The non-detection of polarisation from M51 at 151 MHz is consistent with the estimates of Faraday depolarisation. Future searches for polarised emission in this frequency range should concentrate on regions with low star formation rates.

Wide-band radio polarization observations offer the possibility to recover information about the magnetic fields in synchrotron sources, such as details of their three-dimensional configuration, that has previously been inaccessible. The key physical process involved is the Faraday rotation of the polarized emission in the source (and elsewhere along the wave's propagation path to the observer). In order to proceed, reliable methods are required for inverting the signals observed in wavelength space into useful data in Faraday space, with robust estimates of their uncertainty. In this paper, we examine how variations of the intrinsic angle of polarized emission $\psi_{0}$ with the Faraday depth $\phi$ within a source affect the observable quantities. Using simple models for the Faraday dispersion $F(\phi)$ and $\psi_{0}(\phi)$, along with the current and planned properties of the main radio interferometers, we demonstrate how degeneracies among the parameters describing the magneto-ionic medium can be minimised by combining observations in different wavebands. We also discuss how depolarization by Faraday dispersion due to a random component of the magnetic field attenuates the variations in the spectral energy distribution of the polarization and shifts its peak towards shorter wavelengths. This additional effect reduces the prospect of recovering the characteristics of the magnetic field helicity in magneto-ionic media dominated by the turbulent component of the magnetic field.

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.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in May 2013 as a large-class mission for investigating within the framework of the ESA Cosmic Vision program a set of important scientific questions that require high resolution, high sensitivity, full-sky observations of the sky emission at wavelengths ranging from millimeter-wave to the far-infrared. PRISM's main objective is to explore the distant universe, probing cosmic history from very early times until now as well as the structures, distribution of matter, and velocity flows throughout our Hubble volume. PRISM will survey the full sky in a large number of frequency bands in both intensity and polarization and will measure the absolute spectrum of sky emission more than three orders of magnitude better than COBE FIRAS. The aim of this Extended White Paper is to provide a more detailed overview of the highlights of the new science that will be made possible by PRISM

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.

PRISM (Polarized Radiation Imaging and Spectroscopy Mission) was proposed to ESA in response to the Call for White Papers for the definition of the L2 and L3 Missions in the ESA Science Programme. PRISM would have two instruments: (1) an imager with a 3.5m mirror (cooled to 4K for high performance in the far-infrared---that is, in the Wien part of the CMB blackbody spectrum), and (2) an Fourier Transform Spectrometer (FTS) somewhat like the COBE FIRAS instrument but over three orders of magnitude more sensitive. Highlights of the new science (beyond the obvious target of B-modes from gravity waves generated during inflation) made possible by these two instruments working in tandem include: (1) the ultimate galaxy cluster survey gathering 10e6 clusters extending to large redshift and measuring their peculiar velocities and temperatures (through the kSZ effect and relativistic corrections to the classic y-distortion spectrum, respectively) (2) a detailed investigation into the nature of the cosmic infrared background (CIB) consisting of at present unresolved dusty high-z galaxies, where most of the star formation in the universe took place, (3) searching for distortions from the perfect CMB blackbody spectrum, which will probe a large number of otherwise inaccessible effects (e.g., energy release through decaying dark matter, the primordial power spectrum on very small scales where measurements today are impossible due to erasure from Silk damping and contamination from non-linear cascading of power from larger length scales). These are but a few of the highlights of the new science that will be made possible with PRISM.

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.

We use previously published high-resolution synchrotron polarization data to perform an angular dispersion analysis with the aim of charactering magnetized turbulence in M51. We first analyze three distinct regions (the center of the galaxy, and the northwest and southwest spiral arms) and can clearly discern the turbulent correlation length scale from the width of the magnetized turbulent correlation function for two regions and detect the imprint of anisotropy in the turbulence for all three. Furthermore, analyzing the galaxy as a whole allows us to determine a two-dimensional Gaussian model for the magnetized turbulence in M51. We measure the turbulent correlation scales parallel and perpendicular to the local mean magnetic field to be, respectively, delta_{para} = 98 +/- 5 pc and delta_{perp} = 54 +/- 3 pc, while the turbulent to ordered magnetic field strength ratio is found to be Bt/B0 = 1.01 +/- 0.04. These results are consistent with those of Fletcher et al. (2011), who performed a Faraday rotation dispersion analysis of the same data, and our detection of anisotropy is consistent with current magnetized turbulence theories.

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.

Measuring radio waves at low frequencies offers a new window to study cosmic magnetism, and LOFAR is the ideal radio telescope to open this window widely. The LOFAR Magnetism Key Science Project (MKSP) draws together expertise from multiple fields of magnetism science and intends to use LOFAR to tackle fundamental questions on cosmic magnetism by exploiting a variety of observational techniques. Surveys will provide diffuse emission from the Milky Way and from nearby galaxies, tracking the propagation of long-lived cosmic-ray electrons through magnetic field structures, to search for radio halos around spiral and dwarf galaxies and for magnetic fields in intergalactic space. Targeted deep-field observations of selected nearby galaxies and suspected intergalactic filaments allow sensitive mapping of weak magnetic fields through Rotation Measure (RM) grids. High-resolution observations of protostellar jets and giant radio galaxies reveal structures on small physical scales and at high redshifts, whilst pulsar RMs map large-scale magnetic structures of the Galactic disk and halo in revolutionary detail. The MKSP is responsible for the development of polarization calibration and processing, thus widening the scientific power of LOFAR.

The probability distribution functions (PDFs) of the average densities of the diffuse ionized gas (DIG) and the diffuse atomic gas are close to lognormal, especially when lines of sight at |b|<5 degree and |b|>5 degree are considered separately. Our results provide strong support for the existence of a lognormal density PDF in the diffuse ISM, consistent with a turbulent origin of density structure in the diffuse gas.

Observations of synchrotron radiation and the Faraday rotation of its polarized component allow us to investigate the magnetic properties of the diffuse interstellar medium in nearby galaxies, on scales down to roughly one hundred parsecs. All disc galaxies seem to have a mean, or regular, magnetic field component that is ordered on length scales comparable to the size of the galaxy as well as a random magnetic field of comparable or greater strength. I present an overview of what is currently known observationally about galactic magnetic fields, focusing on the common features among galaxies that have been studied rather than the distinctive or unusual properties of individual galaxies. Of particular interest are the azimuthal patterns formed by regular magnetic fields and their pitch angles as these quantities can be directly related to the predictions of the mean field dynamo theory, the most promising theoretical explanation for the apparent ubiquitous presence of regular magnetic fields in disc galaxies.

(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 new short-period brown dwarf candidate around the star TYC 1240-00945-1. This candidate was discovered in the first year of the Multi-object APO Radial Velocity Exoplanets Large-area Survey (MARVELS), which is part of the third phase of the Sloan Digital Sky Survey (SDSS-III), and we designate the brown dwarf as MARVELS-1b. MARVELS uses the technique of dispersed fixed-delay interferometery to simultaneously obtain radial velocity measurements for 60 objects per field using a single, custom-built instrument that is fiber fed from the SDSS 2.5-m telescope. From our 20 radial velocity measurements spread over a ~370 d time baseline, we derive a Keplerian orbital fit with semi-amplitude K=2.533+/-0.025 km/s, period P=5.8953+/-0.0004 d, and eccentricity consistent with circular. Independent follow-up radial velocity data confirm the orbit. Adopting a mass of 1.37+/-0.11 M_Sun for the slightly evolved F9 host star, we infer that the companion has a minimum mass of 28.0+/-1.5 M_Jup, a semimajor axis 0.071+/-0.002 AU assuming an edge-on orbit, and is probably tidally synchronized. We find no evidence for coherent instrinsic variability of the host star at the period of the companion at levels greater than a few millimagnitudes. The companion has an a priori transit probability of ~14%. Although we find no evidence for transits, we cannot definitively rule them out for companion radii ~<1 R_Jup.

Rotation measures of pulsars and extragalactic point sources have been known to reveal large-scale antisymmetries in the Faraday rotation sky with respect to the Galactic plane and halo that have been interpreted as signatures of the mean magnetic field in the Galactic halo. We describe Faraday rotation measurements of the diffuse Galactic polarized radio emission over a large region in the northern Galactic hemisphere. Through application of Rotation Measure Synthesis we achieve sensitive Faraday rotation maps with high angular resolution, capable of revealing fine-scale structures of about 1 deg in the Faraday rotation sky. Our analysis suggests that the observed antisymmetry in the Faraday rotation sky at b > 0 deg is dominated by the magnetic field around a local HI bubble at a distance of approx. 100 pc, and not by the magnetic field of the Galactic halo. We derive physical properties of the magnetic field of this shell, which we find to be 20 - 34 uG strong. It is clear that the diffuse polarized radio emission contains important information about the local magneto-ionic medium, which cannot yet be derived from Faraday rotation measures of extragalactic sources or pulsars alone.

The new 1-m f/4 fast-slew Zadko Telescope was installed in June 2008 about 70 km north of Perth, Western Australia. It is the only metre-class optical facility at this southern latitude between the east coast of Australia and South Africa, and can rapidly image optical transients at a longitude not monitored by other similar facilities. We report on first imaging tests of a pilot program of minor planet searches, and Target of Opportunity observations triggered by the Swift satellite. In 12 months, 6 gamma-ray burst afterglows were detected, with estimated magnitudes; two of them, GRB 090205 (z = 4.65) and GRB 090516 (z = 4.11), are among the most distant optical transients imaged by an Australian telescope. Many asteroids were observed in a systematic 3-month search. In September 2009, an automatic telescope control system was installed, which will be used to link the facility to a global robotic telescope network; future targets will include fast optical transients triggered by highenergy satellites, radio transient detections, and LIGO gravitational wave candidate events. We also outline the importance of the facility as a potential tool for education, training, and public outreach.

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.

Using high-resolution data of the linearly polarized intensity and polarization angle at 3.6, 6.2, and 20 cm together with a 3-D model of the regular magnetic field, we study variations of the structure, strength, and energy density of the magnetic field in the Scd galaxy M33. The regular magnetic field consists of a horizontal component (represented by an axisymmetric mode from 1 to 3 kpc radius and a superposition of axisymmetric and bisymmetric modes from 3 to 5 kpc radius) and a vertical component. However, the inferred `vertical field' may be partly due to a galactic warp. We estimate the average total and regular magnetic field strengths as ~ 6.4 and 2.5 $\mu$G, respectively. Generation of interstellar magnetic fields by turbulent gas motion in M33 is indicated as the turbulent and magnetic energy densities are about equal.

The Global Magneto-Ionic Medium Survey (GMIMS) is a project to map the diffuse polarized emission over the entire sky, Northern and Southern hemispheres, from 300 MHz to 1.8 GHz. With an angular resolution of 30 - 60 arcmin and a frequency resolution of 1 MHz or better, GMIMS will provide the first spectro-polarimetric data set of the large-scale polarized emission over the entire sky, observed with single-dish telescopes. GMIMS will provide an invaluable resource for studies of the magneto-ionic medium of the Galaxy in the local disk, halo, and its transition.

Using the linearly polarized intensity and polarization angle data at 3.6, 6.2 and 20 cm, we determine variations of Faraday rotation and depolarization across the nearby galaxy M33. A 3-D model of the regular magnetic field is fitted to the observed azimuthal distribution of polarization angles. Faraday rotation, measured between 3.6 and 6.2 cm at a linear resolution of 0.7 kpc, shows more variation in the south than in the north of the galaxy. About 10% of the nonthermal emission from M33 at 3.6 cm is polarized. We estimate the average total and regular magnetic field strengths in M33 as ~ 6.4 and 2.5 $\mu$G, respectively. Under the assumption that the disk of M33 is flat, the regular magnetic field consists of horizontal and vertical components: however the inferred vertical field may be partly due to a galactic warp. The horizontal field is represented by an axisymmetric (m=0) mode from 1 to 3 kpc radius and a superposition of axisymmetric and bisymmetric (m=0+1) modes from 3 to 5 kpc radius. An excess of differential Faraday rotation in the southern half together with strong Faraday dispersion in the southern spiral arms seem to be responsible for the north-south asymmetry in the observed wavelength dependent depolarization. The presence of an axisymmetric m=0 mode of the regular magnetic field in each ring suggests that a galactic dynamo is operating in M33. The pitch angles of the spiral regular magnetic field are generally smaller than the pitch angles of the optical spiral arms but are twice as big as simple estimates based on the mean-field dynamo theory and M33's rotation curve. Generation of interstellar magnetic fields from turbulent gas motions in M33 is indicated by the equipartition of turbulent and magnetic energy densities.

In a search for the signature of turbulence in the diffuse interstellar medium in gas density distributions, we determined the probability distribution functions (PDFs) of the average volume densities of the diffuse gas. The densities were derived from dispersion measures and HI column densities towards pulsars and stars at known distances. The PDFs of the average densities of the diffuse ionized gas (DIG) and the diffuse atomic gas are close to lognormal, especially when lines of sight at |b|<5 degrees and |b|>=5 degrees are considered separately. The PDF of <n_HI> at high |b| is twice as wide as that at low |b|. The width of the PDF of the DIG is about 30 per cent smaller than that of the warm HI at the same latitudes. The results reported here provide strong support for the existence of a lognormal density PDF in the diffuse ISM, consistent with a turbulent origin of density structure in the diffuse gas.

We have developed a technique of isolating elongated structures in galactic images, such as spiral arms, using anisotropic wavelets and apply this to maps of the CO, infrared and radio continuum emission of the grand-design spiral galaxy M51. Systematic shifts between the ridges of CO, infrared and radio continuum emission that are several $\kpc$ long are identified, as well as large variations in pitch angle along spiral arms, of a few tens of degrees. We find two types of arms of polarized radio emission: one has a ridge close to the ridge of CO, with similar pitch angles for the CO and polarization spirals and the regular magnetic field; the other does not always coincide with the CO arm and its pitch angle differs from the orientation of its regular magnetic field. The offsets between ridges of regular magnetic field, dense gas and warm dust are compatible with the sequence expected from spiral density wave triggered star formation, with a delay of a few tens of millions of years between gas entering the shock and the formation of giant molecular clouds and a similar interval between the formation of the clouds and the emergence of young star clusters. At the position of the CO arms the orientation of the regular magnetic field is the same as the pitch angle of the spiral arm, but away from the gaseous arms the orientation of the regular field varies significantly. Spiral shock compression can explain the generation of one type of arm of strong polarized radio emission but a different mechanism is probably responsible for a second type of polarization arm. (Shortened abstract.)

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.

We present the observational results of multi-colour optical monitoring of eight red blazars from 2003 September to 2004 February. The aim of our monitoring is to investigate the spectral variability as well as the flux variations at short and long time scales. The observations were carried out using the 1.0 m robotic telescope of Mt. Lemmon Optical Astronomy Observatory, in Arizona, USA, the 0.6 m telescope of Sobaeksan Optical Astronomy Observatory and the 1.8 m telescope of Bohyunsan Optical Astronomy Observatory, in the Republic of Korea. During the observations, all sources show strong flux variations with amplitudes of larger than 0.5 mag. Variations with amplitudes of over 1 mag are found in four sources. Intraday variations with amplitudes larger than 0.15 mag, and a rapid brightness increase with a rate of ~0.2 mag per day in four days, are detected in S5 0716+71. We investigate the relationship between the colour index and source brightness for each source. We find that two out of three FSRQs tend to be redder when they are brighter, and, conversely, all BL Lac objects tend to be bluer. In particular, we find a significant anti-correlation between the V-I colour index and R magnitude for 3C 454.3. This implies that the spectrum became steeper when the source was brighter, which is opposite to the common trend for blazars. In contrast, significant positive correlations are found in 3C 66A, S5 0716+71, and BL Lac. However, there are only very weak correlations for PKS 0735+17 and OJ 287. We propose that the different relative contributions of the thermal versus non-thermal radiation to the optical emission may be responsible for the different trends of the colour index with brightness in FSRQs and BL Lac objects.

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.

Observations of the southern peculiar galaxy NGC 2442 with the Australia Telescope Compact Array in total and linearly polarized radio continuum at 6cm are presented and compared with previously obtained H$\alpha$ data. The distribution of polarized emission, a signature of regular magnetic fields, reveals some physical phenomena which are unusual among spiral galaxies. We find evidence for tidal interaction and/or ram pressure from the intergalactic medium compressing the magnetic field at the northern and western edges of the galaxy. The radial component of the regular magnetic field in the northern arm is directed away from the centre of the galaxy, a finding which is in contrast to the majority of galaxies studied to date. The oval distortion caused by the interaction generates a sudden jump of the magnetic field pattern upstream of the inner northern spiral arm, similar to galaxies with long bars. An unusual ``island'' of strong regular magnetic field east of the galaxy is probably the brightest part of a magnetic arm frequently seen in spiral galaxies, which appear to be phase-shifted images of the preceding optical arm. The strong magnetic field of the ``island'' indicates a past phase of active star formation when the preceding optical arm was exposed to ram pressure.

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 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.

The radio-loud quasar PMN J0134-0931 was discovered to have an unusual morphology during our search for gravitational lenses. In VLA and MERLIN images, there are 5 compact components with maximum separation 681 millarcseconds. All of these components have the same spectral index from 5 GHz to 43 GHz. In a VLBA image at 1.7 GHz, a curved arc of extended emission joins two of the components in a manner suggestive of gravitational lensing. At least two of the radio components have near-infrared counterparts. We argue that this evidence implies that J0134-0931 is a gravitational lens, although we have not been able to devise a plausible model for the foreground gravitational potential. Like several other radio-loud lenses, the background source has an extraordinarily red optical counterpart.

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.

We report data for $I$ band Surface Brightness Fluctuation (SBF) magnitudes, V-I colors, and distance moduli for 300 galaxies. The Survey contains E, S0 and early-type spiral galaxies in the proportions of 49:42:9, and is essentially complete for E galaxies to Hubble velocities of 2000 km/s, with a substantial sampling of E galaxies out to 4000 km/s. The median error in distance modulus is 0.22 mag. We also present two new results from the Survey. (1) We compare the mean peculiar flow velocity (bulk flow) implied by our distances with predictions of typical cold dark matter transfer functions as a function of scale, and find very good agreement with cold, dark matter cosmologies if the transfer function scale parameter $\Gamma$, and the power spectrum normalization $\sigma_8$ are related by $\sigma_8 \Gamma^{-0.5} \approx 2\pm0.5$. Derived directly from velocities, this result is independent of the distribution of galaxies or models for biasing. The modest bulk flow contradicts reports of large-scale, large-amplitude flows in the $\sim200$ Mpc diameter volume surrounding our Survey volume. (2) We present a distance-independent measure of absolute galaxy luminosity, \Nbar, and show how it correlates with galaxy properties such as color and velocity dispersion, demonstrating its utility for measuring galaxy distances through large and unknown extinction.

The structure of the radio source RC B1146+052 of the ``Cold'' catalogue is investigated by data of the MIT-GB-VLA survey at 4850 MHz. This source belongs to the steep spectrum radio sources subsample of the RC catalogue. Its spectral index is $\alpha$ = -1.04. The optical image of this source obtained with 6m telescope is analysed. The radio source center is situated in a group of 8 galaxies of about 24$^m$ in the R-filter. The possible explanations of the complex structure of radio components are considered.

We report the discovery of a binary quasar, MGC2214+3550 A,B, whose components have similar optical spectra at a redshift z=0.88. The quasars are separated on the sky by 3.0", and have a magnitude difference of Delta(m_I)=0.5 mag. The VLA radio map at 3.6 cm shows a single 47mJy radio source with a core-jet morphology that is coincident with the brighter optical quasar A. Gravitational lensing is ruled out by the lack of radio emission from quasar B, and the lack of any visible galaxies to act as the lens. We conclude that MGC 2214+3550 A and B are physically associated. With a projected separation of 12.7h^{-1} kpc (Omega_0=1), MGC 2214+3550 A,B is one of the smallest z>0.5 binary quasars.

We present radio and optical observations of MG0248+0641, which contains a kiloparsec-scale ring-like structure in one of its radio lobes. The radio observations show a typical core-double morphology: a central core between two lobes, each of which has a hotspot. The western radio lobe appears as a nearly continuous ring, with linear polarization electric field vectors which are oriented in a radial direction from the ring center. We consider several different interpretations for the nature of this ring, including gravitational lensing of a normal jet by a foreground galaxy. Even though simple lensing models can describe the ring morphology reasonably well, the high linear polarization seen around the ring cannot be easily explained. The chance interposition of a galactic supernova remnant, nova, planetary nebula, or H II region, has been ruled out. The highly polarized ring of MG0248+0641 is much like the prominent ring seen in 3C219, and the multiple ones in 3C310 and Hercules A, suggesting that similar physical processes are producing shell structures in these radio galaxies. The ring in MG0248+0641 may be caused by the formation of ``bubbles'', as a result of instabilities in the energy flow down the western radio jet. It may also be possible that the required instabilities are triggered by the infall of gas, via tidal interaction of the central source with a nearby galaxy. This scenario may be indicated by our marginal detection of an optical source close to the western hotspot.

We report the discovery of a new gravitationally lensed radio source. Radio maps of MG0751+2716 show four lensed images, which, at higher resolution, are resolved into long arcs of emission. A group of galaxies is present in optical images, including the principal lensing galaxy, with a much brighter galaxy just a few arcseconds away. We have measured the redshift of this brighter galaxy. No optical counterpart to the background source has been detected. Lens models that can readily reproduce the lensed image positions all require a substantial shear component. However, neither the very elongated lens nor the bright nearby galaxy are correctly positioned to explain the shear. Lens models which associate the mass with the light of galaxies in the group can produce an acceptable fit, but only with an extreme mass-to-light ratio in one of the minor group members.