The Square Kilometre Array (SKA) will answer fundamental questions about the origin, evolution, properties, and influence of magnetic fields throughout the Universe. Magnetic fields can illuminate and influence phenomena as diverse as star formation, galactic dynamics, fast radio bursts, active galactic nuclei, large-scale structure, and Dark Matter annihilation. Preparations for the SKA are swiftly continuing worldwide, and the community is making tremendous observational progress in the field of cosmic magnetism using data from a powerful international suite of SKA pathfinder and precursor telescopes. In this contribution, we revisit community plans for magnetism research using the SKA, in the light of these recent rapid developments. We focus in particular on the impact that new radio telescope instrumentation is generating, thus advancing our understanding of key SKA magnetism science areas, as well as the new techniques that are required for processing and interpreting the data. We discuss these recent developments in the context of the ultimate scientific goals for the SKA era.

Measuring the properties of extragalactic magnetic fields through the effect of Faraday rotation provides a means to understand the origin and evolution of cosmic magnetism. Here we use data from the LOFAR Two-Metre Sky Survey (LoTSS) to calculate the Faraday rotation measure (RM) of close pairs of extragalactic radio sources. By considering the RM difference ($\Delta$RM) between physical pairs (e.g. double-lobed radio galaxies) and non-physical pairs (i.e. close projected sources on the sky), we statistically isolate the contribution of extragalactic magnetic fields to $\Delta$RM along the line of sight between non-physical pairs. From our analysis, we find no significant difference between the $\Delta$RM distributions of the physical and non-physical pairs, limiting the excess Faraday rotation contribution to $< 1.9$ rad/m$^2$ ($\sim$$95\%$ confidence). We use this limit with a simple model of an inhomogeneous universe to place an upper limit of 4 nG on the cosmological co-moving magnetic field strength on Mpc scales. We also compare the RM data with a more realistic suite of cosmological MHD simulations, that explore different magnetogenesis scenarios. Both magnetization of the large scale structure by astrophysical processes such as galactic and AGN outflows, and simple primordial scenarios with seed magnetic field strengths $< 0.5$ nG cannot be rejected by the current data; while stronger primordial fields or models with dynamo amplification in filaments are disfavoured.

We present LOFAR observations at 150 MHz of the borderline FRI/FRII giant radio galaxy NGC 6251. This paper presents the most sensitive and highest-resolution images of NGC 6251 at these frequencies to date, revealing for the first time a low-surface-brightness extension to the northern lobe, and a possible backflow associated with the southern lobe. The integrated spectra of components of NGC 6251 are consistent with previous measurements at higher frequencies, similar to results from other LOFAR studies of nearby radio galaxies. We find the outer structures of NGC 6251 to be either at equipartition or slightly electron dominated, similar to those of FRII sources rather than FRIs; but this conclusion remains tentative because of uncertainties associated with the geometry and the extrapolation of X-ray measurements to determine the external pressure distribution on the scale of the outer lobes. We place lower limits on the ages of the extension of the northern lobe and the backflow of the southern lobe of $t \gtrsim 250$ Myr and $t \gtrsim 210$ Myr respectively. We present the first detection of polarisation at 150 MHz in NGC 6251. Taking advantage of the high Faraday resolution of LOFAR, we place an upper limit on the magnetic field in the group of $B < 0.2 (\Lambda_B / 10 {\rm kpc})^{-0.5} \mu$G for a coherence scale of $\Lambda_B < 60 {\rm kpc}$ and $B < 13 \mu$G for $\Lambda_B = 240$ kpc.

Giant radio galaxies (GRGs) are physically large radio sources that extend well beyond their host galaxy environment. Their polarization properties are affected by the poorly constrained magnetic field that permeates the intergalactic medium on Mpc scales. A low frequency ($<$ 200 MHz) polarization study of this class of radio sources is now possible with LOFAR. Here we investigate the polarization properties and Faraday rotation measure (RM) of a catalog of GRGs detected in the LoTSS. This is the first low-frequency polarization study of a large sample of radio galaxies selected on their physical size. We explore the magneto-ionic properties of their under-dense environment and probe intergalactic magnetic fields using the Faraday rotation properties of their radio lobes. We use RM synthesis in the 120-168 MHz band to search for polarized emission and to derive the RM and fractional polarization of each detected source component. We study the depolarization between 1.4 GHz and 144 MHz using images from the NVSS. From a sample of 240 GRGs, we detected 37 sources in polarization, all with a total flux density above 56 mJy. The fractional polarization of the detected GRGs at 1.4 GHz and 144 MHz is consistent with a small amount of Faraday depolarization (a Faraday dispersion $<$ 0.3 rad m$^{-2}$). Our analysis shows that the lobes are expanding into a low-density ($<10^{-5}$ cm$^{-3}$) local environment permeated by weak magnetic fields ($<$0.1 $\mu$G) with fluctuations on scales of 3 to 25 kpc. The presence of foreground galaxy clusters appears to influence the polarization detection rate up to 2R$_{500}$. In general, this work demonstrates the ability of LOFAR to quantify the rarefied environments in which these GRGs exist and highlights them as an excellent statistical sample to use as high precision probes of magnetic fields in the intergalactic medium and the Milky Way.

Low frequency ($\nu\lesssim 150\,$MHz) stellar radio emission is expected to originate in the outer corona at heights comparable to and larger than the stellar radius. Such emission from the Sun has been used to study coronal structure, mass ejections, space-weather conditions around the planets (Schwenn 2006). Searches for low-frequency emission from other stars have only detected a single active flare-star (Lynch et al. 2017) that is not representative of the wider stellar population. Here we report the detection of low-frequency radio emission from a quiescent star, GJ 1151--- a member of the most common stellar type (red dwarf or spectral class M) in the Galaxy. The characteristics of the emission are similar to those of planetary auroral emissions (Zarka 1998) (e.g. Jupiter's decametric emission), suggesting a coronal structure dominated by a global magnetosphere with low plasma density. Our results show that large-scale currents that power radio aurorae operate over a vast range of mass and atmospheric composition, ranging from terrestrial planets to main-sequence stars. The Poynting flux required to produce the observed radio emission cannot be generated by GJ 1151's slow rotation, but can originate in a sub-Alfv\'{e}nic interaction of its magnetospheric plasma with a short-period exoplanet. The emission properties are consistent with theoretical expectations (Zarka 2007; Lanza 2009; Saur et al. 2013; Turnpenney et al. 2018) for interaction with an Earth-size planet in a $\sim 1-5$ day-long orbit.

The Very Large Array Sky Survey (VLASS) is a synoptic, all-sky radio sky survey with a unique combination of high angular resolution ($\approx$2.5"), sensitivity (a 1$\sigma$ goal of 70 $\mu$Jy/beam in the coadded data), full linear Stokes polarimetry, time domain coverage, and wide bandwidth (2-4 GHz). The first observations began in September 2017, and observing for the survey will finish in 2024. VLASS will use approximately 5500 hours of time on the Karl G. Jansky Very Large Array (VLA) to cover the whole sky visible to the VLA (Declination $>-40^{\circ}$), a total of 33,885 deg$^2$. The data will be taken in three epochs to allow the discovery of variable and transient radio sources. The survey is designed to engage radio astronomy experts, multi-wavelength astronomers, and citizen scientists alike. By utilizing an "on the fly" interferometry mode, the observing overheads are much reduced compared to a conventional pointed survey. In this paper, we present the science case and observational strategy for the survey, and also results from early survey observations.

We studied the polarization behavior of the quasar 3C 273 over the 1 mm wavelength band at ALMA with a total bandwidth of 7.5GHz across 223 to 243 GHz at 0.8 arcsec resolution, corresponding to 2.1 kpc at the distance of 3C 273. With these observations we were able to probe the optically thin polarized emission close to the jet base, and constrain the magnetic field structure. We computed the Faraday rotation measure using simple linear fitting and Faraday rotation measure synthesis. In addition, we modeled the broadband behavior of the fractional Stokes Q and U parameters (qu-fitting). The systematic uncertainties in the polarization observations at ALMA were assessed through Monte Carlo simulations. We find the unresolved core of 3C 273 to be 1.8% linearly polarized. We detect a very high rotation measure (RM) of ~5.0 x10^5 rad/m^2 over the 1 mm band when assuming a single polarized component and an external RM screen. This results in a rotation of >40 deg of the intrinsic electric vector position angle, which is significantly higher than typically assumed for millimeter wavelengths. The polarization fraction increases as a function of wavelength, which according to our qu-fitting could be due to multiple polarized components of different Faraday depth within our beam, or internal Faraday rotation. With our limited wavelength coverage, we cannot distinguish between the cases, and additional multifrequency and high angular resolution observations are needed to determine the location and structure of the magnetic field of the Faraday active region. Comparing our RM estimate with values obtained at lower frequencies, the RM increases as a function of observing frequency, following a power law with an index of ~2.0 consistent with a sheath surrounding a conically expanding jet. We also detect ~0.2% circular polarization, although further observations are needed to confirm this result.

We investigate the previously proposed possibility that multi-epoch broadband polarimetry could act as a complement or limited proxy for VLBI observations of blazars, in that the number of polarised emission components in the jet, and some of their properties and those of the foreground environment, might be inferred from the object's time-varying 1D Faraday depth spectrum (FDS) alone. We report on a pilot-scale experiment designed to establish the basic plausibility and utility of this idea. We analyse temporal changes in the complex polarisation spectra of nine spatially unresolved (at arcsecond scales) blazars in two epochs separated by $\sim$5 years, using data taken with the Australia Telescope Compact Array. The data allow for precise modelling, and we demonstrate that all objects in our sample show changes in their polarisation spectrum that cannot be accounted for by uncertainties in calibration or observational effects. By associating polarised emission components across epochs, we infer changes in their number, intrinsic fractional polarisation, intrinsic polarisation angle, rotation measure, and depolarisation characteristics. We attribute these changes to evolution in the structure of the blazar jets, most likely located at distances of up to tens of parsecs from the central active galactic nuclei. Our results suggest that continued work in this area is warranted; in particular, it will be important to determine the frequency ranges and temporal cadence most useful for scientifically exploiting the effects.

The technique of Faraday tomography is a key tool for the study of magnetised plasmas in the new era of broadband radio polarisation observations. In particular, observations at metre-wavelengths provide significantly better Faraday depth accuracies compared to traditional cm-wavelength observations. However, the effect of Faraday depolarisation makes the polarised signal very challenging to detect at metre wavelengths (MHz frequencies). In this work, Faraday tomography is used to characterise the Faraday rotation properties of polarised sources found in data from the LOFAR Two-Metre Sky Survey (LoTSS). Of the 76 extragalactic polarised sources analysed here, we find that all host a radio-loud AGN. The majority of the sources ($\sim$64%) are large FRII radio galaxies with a median projected linear size of 710 kpc and median radio luminosity at 144 MHz of $4\times10^{26}$ W Hz$^{-1}$ (with $\sim$13% of all sources having a linear size $>1$ Mpc). In several cases, both hotspots are detected in polarisation at an angular resolution of $\sim$20". One such case allowed a study of intergalactic magnetic fields on scales of 3.4 Mpc. Other detected source types include an FRI radio galaxy and at least 8 blazars. Most sources display simple Faraday spectra, however, we highlight one blazar that displays a complex Faraday spectrum, with two close peaks in the Faraday dispersion function.

We present observations of linear polarisation in the southern radio lobe of Centaurus A, conducted during commissioning of the Australian Square Kilometre Array Pathfinder (ASKAP) telescope. We used 16 antennas to observe a 30 square degree region in a single 12 hour pointing over a 240 MHz band centred on 913 MHz. Our observations achieve an angular resolution of $26\times33$ arcseconds (480 parsecs), a maximum recoverable angular scale of 30 arcminutes, and a full-band sensitivity of 85 $\muup$Jy beam$^{-1}$. The resulting maps of polarisation and Faraday rotation are amongst the most detailed ever made for radio lobes, with of order 10$^5$ resolution elements covering the source. We describe several as-yet unreported observational features of the lobe, including its detailed peak Faraday depth structure, and intricate networks of depolarised filaments. These results demonstrate the exciting capabilities of ASKAP for widefield radio polarimetry.

Cosmological simulations predict that an intergalactic magnetic field (IGMF) pervades the large scale structure (LSS) of the Universe. Measuring the IGMF is important to determine its origin (i.e. primordial or otherwise). Using data from the LOFAR Two Metre Sky Survey (LoTSS), we present the Faraday rotation measure (RM) and depolarisation properties of the giant radio galaxy J1235+5317, at a redshift of $z = 0.34$ and 3.38 Mpc in size. We find a mean RM difference between the lobes of $2.5\pm0.1$ rad/m$^2$ , in addition to small scale RM variations of ~0.1 rad/m$^2$ . From a catalogue of LSS filaments based on optical spectroscopic observations in the local universe, we find an excess of filaments intersecting the line of sight to only one of the lobes. Associating the entire RM difference to these LSS filaments leads to a gas density-weighted IGMF strength of ~0.3 {\mu}G. However, direct comparison with cosmological simulations of the RM contribution from LSS filaments gives a low probability (~5%) for an RM contribution as large as 2.5 rad/m$^2$ , for the case of IGMF strengths of 10 to 50 nG. It is likely that variations in the RM from the Milky Way (on 11' scales) contribute significantly to the mean RM difference, and a denser RM grid is required to better constrain this contribution. In general, this work demonstrates the potential of the LOFAR telescope to probe the weak signature of the IGMF. Future studies, with thousands of sources with high accuracy RMs from LoTSS, will enable more stringent constraints on the nature of the IGMF.

The LOFAR Two-metre Sky Survey (LoTSS) is an ongoing sensitive, high-resolution 120-168MHz survey of the entire northern sky for which observations are now 20% complete. We present our first full-quality public data release. For this data release 424 square degrees, or 2% of the eventual coverage, in the region of the HETDEX Spring Field (right ascension 10h45m00s to 15h30m00s and declination 45$^\circ$00$'$00$''$ to 57$^\circ$00$'$00$''$) were mapped using a fully automated direction-dependent calibration and imaging pipeline that we developed. A total of 325,694 sources are detected with a signal of at least five times the noise, and the source density is a factor of $\sim 10$ higher than the most sensitive existing very wide-area radio-continuum surveys. The median sensitivity is S$_{\rm 144 MHz} = 71\,\mu$Jy beam$^{-1}$ and the point-source completeness is 90% at an integrated flux density of 0.45mJy. The resolution of the images is 6$''$ and the positional accuracy is within 0.2$''$. This data release consists of a catalogue containing location, flux, and shape estimates together with 58 mosaic images that cover the catalogued area. In this paper we provide an overview of the data release with a focus on the processing of the LOFAR data and the characteristics of the resulting images. In two accompanying papers we provide the radio source associations and deblending and, where possible, the optical identifications of the radio sources together with the photometric redshifts and properties of the host galaxies. These data release papers are published together with a further $\sim$20 articles that highlight the scientific potential of LoTSS.

We constructed a sample of 23,344 radio-loud active galactic nuclei (RLAGN) from the catalogue derived from the LOFAR Two-Metre Sky Survey (LoTSS) survey of the HETDEX Spring field. Although separating AGN from star-forming galaxies remains challenging, the combination of spectroscopic and photometric techniques we used gives us one of the largest available samples of candidate RLAGN. We used the sample, combined with recently developed analytical models, to investigate the lifetime distribution of RLAGN. We show that large or giant powerful RLAGN are probably the old tail of the general RLAGN population, but that the low-luminosity RLAGN candidates in our sample, many of which have sizes $<100$ kpc, either require a very different lifetime distribution or have different jet physics from the more powerful objects. We then used analytical models to develop a method of estimating jet kinetic powers for our candidate objects and constructed a jet kinetic luminosity function based on these estimates. These values can be compared to observational quantities, such as the integrated radiative luminosity of groups and clusters, and to the predictions from models of RLAGN feedback in galaxy formation and evolution. In particular, we show that RLAGN in the local Universe are able to supply all the energy required per comoving unit volume to counterbalance X-ray radiative losses from groups and clusters and thus prevent the hot gas from cooling. Our computation of the kinetic luminosity density of local RLAGN is in good agreement with other recent observational estimates and with models of galaxy formation.

An understanding of the relationship between radio-loud active galaxies and their large-scale environments is essential for realistic modelling of radio-galaxy evolution and environmental impact, for understanding AGN triggering and life cycles, and for calibrating galaxy feedback in cosmological models. We use the LOFAR Two-Metre Sky Survey (LoTSS) Data Release 1 catalogues to investigate this relationship. We cross-matched 8,745 radio-loud AGN with 0.08<z<0.4, selected from LoTSS, with two Sloan Digital Sky Survey cluster catalogues, and find that only 10 percent of LoTSS AGN in this redshift range have an association, so that the majority of low-redshift AGN (including a substantial fraction of the most radio-luminous objects) must inhabit haloes with M < 10^14 M_sun. We find that the probability of a cluster association, and the richness of the associated cluster, is correlated with AGN radio luminosity, and for the cluster population, the number of associated AGN and the radio luminosity of the brightest associated AGN is richness-dependent. We demonstrate that these relations are not driven solely by host-galaxy stellar mass, supporting models in which large-scale environment is influential in driving AGN jet activity. At the lowest radio luminosities we find that the minority of objects with a cluster association are located at larger mean cluster-centre distances than more luminous AGN, which appears to be driven primarily by host-galaxy mass. Finally, we also find that FRI radio galaxies inhabit systematically richer environments than FRIIs. The work presented here demonstrates the potential of LoTSS for AGN environmental studies. In future, the full northern-sky LoTSS catalogue, together with the use of deeper optical/IR imaging data and spectroscopic follow-up with WEAVE-LOFAR, will provide opportunities to extend this type of work to much larger samples and higher redshifts.

We present a polarization and Faraday rotation study of the hybrid morphology radio galaxy NGC 612 (PKS 0131-36), using Australian Telescope Compact Array observations from 1 to 3 GHz. In general, the results are consistent with an external Faraday screen close to the radio source. In the eastern FRII lobe, the RM of the hotspot increases in magnitude towards the leading edge, as well as changing sign (compared to the rest of the lobe). The Faraday depolarization is also ~3 times larger at the hotspot than elsewhere. A plausible explanation for this is significant compression of ambient magnetised gas by the bow shock produced by the advancing hotspot. The western FRI lobe also exhibits some evidence of interaction with local magnetised gas, as a transverse band of high RM coincides with a distinct bend in the lobe. Previous observations of NGC 612 revealed an HI bridge of tidal debris along the direction of the eastern lobe towards the gas-rich companion NGC 619. We find no clear evidence that ionised gas associated with this bridge is either mixing with or lies in the foreground of the radio source. This is consistent with the absence of HI absorption against the hotspot, and indicates that the tidal debris must lie mostly behind the eastern lobe.

Studying the interaction between AGN jets and lobes and their surrounding environment is important in order to understand how they transfer energy to their environment as well as determining the intrinsic physical properties of the sources themselves. This paper presents broadband VLA polarization and Faraday rotation observations of the radio galaxy Coma A (3C 277.3) from 1 to 4 GHz, including archival VLA observations at 4.9 and 15 GHz. Through broadband polarization model-fitting, we find that an external Faraday screen with a turbulent magnetic field provides an appropriate description to most of the data. By combining the polarization and Faraday rotation results with previous H$\alpha$ observations, we identified the H$\alpha$-emitting gas as the Faraday screen responsible for the observed Faraday depolarization. We were able to derive the magnetic field strength in the H$\alpha$-emitting gas, finding typical field strengths of $\sim1$ $\mu$G, which is consistent with studies of the intra-group medium local to other radio galaxies. However, we find a highly depolarized region of the southern lobe coincident with a H$\alpha$ filament that has a field strength comparable to the equipartition field strength in the radio lobe (i.e. $\gtrsim$36 $\mu$G). This implies that the H$\alpha$ filament is internal to the radio emitting plasma. Such clear examples of internal Faraday depolarization are rare, thus providing another key insight into the evolution of radio galaxies and their ability to provide significant feedback on the local gas that would otherwise cool and form stars.

Extended, steep-spectrum radio synchrotron sources are pre-dominantly found in massive galaxy clusters as opposed to groups. LOFAR Two-Metre Sky Survey images have revealed a diffuse, ultra-steep spectrum radio source in the low-mass cluster Abell 1931. The source has a fairly irregular morphology with a largest linear size of about 550 kpc. The source is only seen in LOFAR observations at 143 MHz and GMRT observations at 325 MHz. The spectral index of the total source between 143 MHz and 325 MHz is $\alpha_{143}^{325} = -2.86 \pm 0.36$. The source remains invisible in Very Large Array (1-2 GHz) observations as expected given the spectral index. Chandra X-ray observations of the cluster revealed a bolometric luminosity of $L_X = (1.65 \pm 0.39) \times 10^{43}$ erg s$^{-1}$ and a temperature of $2.92_{-0.87}^{+1.89}$ keV which implies a mass of around $\sim 10^{14} M_{\odot}$. We conclude that the source is a remnant radio galaxy that has shut off around 200 Myr ago. The brightest cluster galaxy, a radio-loud elliptical galaxy, could be the source for this extinct source. Unlike remnant sources studied in the literature, our source has a steep spectrum at low radio frequencies. Studying such remnant radio galaxies at low radio frequencies is important for understanding the scarcity of such sources and their role in feedback processes.

We present observations and analysis of the polarized radio emission from the nearby radio galaxy Fornax A over 1.28--3.1 GHz, using data from the Australia Telescope Compact Array (ATCA). In this, the first of two associated papers, we use modern broadband polarimetric techniques to examine the nature and origin of conspicuous low-polarization (low-$p$) patches in the lobes. We resolve the low-$p$ patches, and find that their low fractional polarization is associated with complicated frequency-dependent interference in the polarized signal generated by Faraday effects along the line of sight. The low-$p$ patches are spatially correlated with interfaces in the magnetic structure of the lobe, across which the line-of-sight-projected magnetic field changes direction. Spatial correlations with the sky-projected magnetic field orientation and structure in total intensity are also identified and discussed. We argue that the low-$p$ patches, along with associated reversals in the line-of-sight magnetic field and other related phenomena, are best explained by the presence of $\mathcal{O}(10^9)$ $M_\odot$ of magnetized thermal plasma in the lobes, structured in shells or filaments, and likely advected from the ISM of NCG 1316 or its surrounding ICM. Our study underscores the power and utility of spatially-resolved, broadband, full-polarization radio observations to reveal new facets of flow behaviors and magneto-ionic structure in radio lobes and their interplay with the surrounding environment.

We present broadband polarimetric observations of a sample of high Faraday rotation measure (RM) AGN using the VLA telescope from 1 to 2 GHz, and 4 to 12 GHz. The sample (14 sources) consists of compact sources (linear resolution smaller than ~ 5 kpc) that are unpolarized at 1.4 GHz in the NVSS. Total intensity data have been modelled using combination of synchrotron components, revealing complex structure in their radio spectra. Depolarization modelling, through the so called qu-fitting, have been performed on the polarized data using an equation that attempts to simplify the process of fitting many different depolarization models that we can divide into two major categories: External Depolarization and Internal Depolarization models. Understanding which of the two mechanisms are the most representative, would help the qualitative understanding of the AGN jet environment, whether it is embedded in a dense external magnetoionic medium or if it is the jet-wind that causes the high RM and strong depolarization. This could help to probe the jet magnetic field geometry (e.g. helical or otherwise). This new high-sensitivity data, shows a complicated behaviour in the total intensity and polarization radio spectrum of individual sources. We observed the presence of several synchrotron components and Faraday components in their total intensity and polarized spectra. For the majority of our targets, (12 sources) the depolarization seems to be caused by a turbulent magnetic field. Thus, our main selection criteria (lack of polarization at 1.4 GHz in the NVSS), results in a sample of sources with very large RMs and depolarization due to turbulent magnetic fields local to the source. We show how the new qu-fitting technique can be used to probe the magnetised radio source environment and to spectrally resolve the polarized components of unresolved radio sources.

We present full-polarisation, broadband observations of the radio galaxy NGC 612 (PKS B0131-637) from 1.3 to 3.1 GHz using the Australia Telescope Compact Array. The relatively large angular scale of the radio galaxy makes it a good candidate with which to investigate the polarisation mechanisms responsible for the observed Faraday depth structure. By fitting complex polarisation models to the polarised spectrum of each pixel, we find that a single polarisation component can adequately describe the observed signal for the majority of the radio galaxy. While we cannot definitively rule out internal Faraday rotation, we argue that the bulk of the Faraday rotation is taking place in a thin skin that girts the polarised emission. Using minimum energy estimates, we find an implied total magnetic field strength of 4.2 microG.

The polarization properties of radio sources at very low frequencies (<200 MHz) have not been widely measured, but the new generation of low-frequency radio telescopes, including the Low Frequency Array (LOFAR: a Square Kilometre Array Low pathfinder), now gives us the opportunity to investigate these properties. In this paper, we report on the preliminary development of a data reduction pipeline to carry out polarization processing and Faraday tomography for data from the LOFAR Two-meter Sky Survey (LOTSS) and present the results of this pipeline from the LOTSS preliminary data release region (10h45m - 15h30m right ascension, 45 - 57 degrees declination, 570 square degrees). We have produced a catalog of 92 polarized radio sources at 150 MHz at 4.3 arcminute resolution and 1 mJy rms sensitivity, which is the largest catalog of polarized sources at such low frequencies. We estimate a lower limit to the polarized source surface density at 150 MHz, with our resolution and sensitivity, of 1 source per 6.2 square degrees. We find that our Faraday depth measurements are in agreement with previous measurements and have significantly smaller errors. Most of our sources show significant depolarization compared to 1.4 GHz, but there is a small population of sources with low depolarization indicating that their polarized emission is highly localized in Faraday depth. We predict that an extension of this work to the full LOTSS data would detect at least 3400 polarized sources using the same methods, and probably considerably more with improved data processing.

We present a low-frequency, broadband polarization study of the FRII radio galaxy PKS J0636-2036 (z = 0.0551), using the Murchison Widefield Array (MWA) from 70 to 230 MHz. The northern and southern hotspots (separated by ~14.5' on the sky) are resolved by the MWA (3'.3 resolution) and both are detected in linear polarization across the full frequency range. A combination of Faraday rotation measure (RM) synthesis and broadband polarization model-fitting are used to constrain the Faraday depolarization properties of the source. For the integrated southern hotspot emission, two RM component models are strongly favoured over a single RM component, and the best-fitting model requires Faraday dispersions of approximately 0.7 and 1.2 rad/m$^2$ (with a mean RM of ~50 rad/m$^2$). High resolution imaging at 5" with the ATCA shows significant sub-structure in the southern hotspot and highlights some of the limitations in the polarization modelling of the MWA data. Based on the observed depolarization, combined with extrapolations of gas density scaling-relations for group environments, we estimate magnetic field strengths in the intergalactic medium between ~0.04 and 0.5 {\mu}G. We also comment on future prospects of detecting more polarized sources at low frequencies.

We present the largest ever sample of radio polarization properties for $z>4$ sources, with 14 sources having significant polarization detections. Using wideband data from the Karl G. Jansky Very Large Array, we obtained the rest-frame total intensity and polarization properties of 37 radio sources, nine of which have spectroscopic redshifts in the range $1 \le z \le 1.4$, with the other 28 having spectroscopic redshifts in the range $3.5 \le z \le 6.21$. Fits are performed for the Stokes $I$ and fractional polarization spectra, and Faraday rotation measures are derived using Rotation measure synthesis and $QU$ fitting. Using archival data of 476 polarized sources, we compare high redshift ($z>3$) source properties to a $15\,$GHz rest-frame luminosity matched sample of low redshift ($z<3$) sources to investigate if the polarization properties of radio sources at high redshifts are intrinsically different than those at low redshift. We find a mean of the rotation measure absolute values, corrected for Galactic rotation, of $50 \pm 22\,$rad m$^{-2}$ for $z>3$ sources and $57 \pm 4\,$rad m$^{-2}$ for $z<3$. Although there is some indication of lower intrinsic rotation measures at high-$z$ possibly due to higher depolarization from the high density environments, using several statistical tests we detect no significant difference between low and high redshift sources. Larger samples are necessary to determine any true physical difference.

We present techniques developed to calibrate and correct Murchison Widefield Array (MWA) low frequency (72-300 MHz) radio observations for polarimetry. The extremely wide field-of-view, excellent instantaneous (u, v)-coverage and sensitivity to degree-scale structure that the MWA provides enable instrumental calibration, removal of instrumental artefacts, and correction for ionospheric Faraday rotation through imaging techniques. With the demonstrated polarimetric capabilities of the MWA, we discuss future directions for polarimetric science at low frequencies to answer outstanding questions relating to polarised source counts, source depolarisation, pulsar science, low-mass stars, exoplanets, the nature of the interstellar and intergalactic media, and the solar environment.

We present the results from a broadband (1 to 3 GHz), spectro-polarimetry study of the integrated emission from 100 extragalactic radio sources with the ATCA, selected to be highly linearly polarized at 1.4 GHz. We use a general purpose, polarization model-fitting procedure that describes the Faraday rotation measure (RM) and intrinsic polarization structure of up to three distinct polarized emission regions or 'RM components' of a source. Overall, 37%/52%/11% of sources are best fit by one/two/three RM components. However, these fractions are dependent on the signal-to-noise ratio (S/N) in polarization (more RM components more likely at higher S/N). In general, our analysis shows that sources with high integrated degrees of polarization at 1.4 GHz have low Faraday depolarization, are typically dominated by a single RM component, have a steep spectral index, and a high intrinsic degree of polarization. After classifying our sample into radiative-mode and jet-mode AGN, we find no significant difference between the Faraday rotation or Faraday depolarization properties of jet-mode and radiative-mode AGN. However, there is a statistically significant difference in the intrinsic degree of polarization between the two types, with the jet-mode sources having more intrinsically ordered magnetic field structures than the radiative-mode sources. We also find a preferred perpendicular orientation of the intrinsic magnetic field structure of jet-mode AGN with respect to the jet direction, while no clear preference is found for the radiative-mode sources.

In previous work, we developed tools for quantifying the tidal distortion of a black hole's event horizon due to an orbiting companion. These tools use techniques which require large mass ratios (companion mass $\mu$ much smaller than black hole mass $M$), but can be used for arbitrary bound orbits, and for any black hole spin. We also showed how to visualize these distorted black holes by embedding their horizons in a global Euclidean 3-space, ${\mathbb{E}}^3$. Such visualizations illustrate interesting and important information about horizon dynamics. Unfortunately, we could not visualize black holes with spin parameter $a_* > \sqrt{3}/2 \approx 0.866$: such holes cannot be globally embedded into ${\mathbb{E}}^3$. In this paper, we overcome this difficulty by showing how to embed the horizons of tidally distorted Kerr black holes in a hyperbolic 3-space, ${\mathbb{H}}^3$. We use black hole perturbation theory to compute the Gaussian curvatures of tidally distorted event horizons, from which we build a two-dimensional metric of their distorted horizons. We develop a numerical method for embedding the tidally distorted horizons in ${\mathbb{H}}^3$. As an application, we give a sequence of embeddings into ${\mathbb{H}}^3$ of a tidally interacting black hole with spin $a_*=0.9999$. A small amplitude, high frequency oscillation seen in previous work shows up particularly clearly in these embeddings.

The second-order extended stability Factorized Runge-Kutta-Chebyshev (FRKC2) class of explicit schemes for the integration of large systems of PDEs with diffusive terms is presented. FRKC2 schemes are straightforward to implement through ordered sequences of forward Euler steps with complex stepsizes, and easily parallelised for large scale problems on distributed architectures. Preserving 7 digits for accuracy at 16 digit precision, the schemes are theoretically capable of maintaining internal stability at acceleration factors in excess of 6000 with respect to standard explicit Runge-Kutta methods. The stability domains have approximately the same extents as those of RKC schemes, and are a third longer than those of RKL2 schemes. Extension of FRKC methods to fourth-order, by both complex splitting and Butcher composition techniques, is discussed. A publicly available implementation of the FRKC2 class of schemes may be obtained from maths.dit.ie/frkc

We present preliminary results of JVLA wideband full polarization observations of a sample of Active Galactic Nuclei (AGN) with very high Rotation Measure (RM) values, a sign of extreme environment. Polarization properties show a complex behaviour such that the polarization angle (PA) and fractional polarization (fp) change dramatically within the wide band. The measured RM is not constant within the wide band. Its complex behaviour reflects the complexity of the medium with the presence of several Faraday components. The depolarization has been studied by modelling the variations of the Stokes parameters Q and U together with the polarization parameters (PA and fp) with wavelength using combinations of the simplest existing depolarization models. With this JVLA study we could spectrally resolve multiple polarized components of unresolved AGN. These preliminary results reveal the complexity of these objects, but improvements to the depolarization modelling are needed to better understand the polarization structure of these sources.

Protogalactic environments are typically identified using quasar absorption lines, and these galactic building blocks can manifest as Damped Lyman-Alpha Absorbers (DLAs) and Lyman Limit Systems (LLSs). We use radio observations of Faraday effects to test whether DLAs and LLSs host a magnetised medium, by combining DLA and LLS detections throughout the literature with 1.4 GHz polarization data from the NRAO VLA Sky Survey (NVSS). We obtain a control, a DLA, and a LLS sample consisting of 114, 19, and 27 lines-of-sight respectively - all of which are polarized at $\ge8\sigma$ to ensure Rician bias is negligible. Using a Bayesian framework, we are unable to detect either coherent or random magnetic fields in DLAs: the regular coherent magnetic fields within the DLAs must be $\le2.8$ $\mu$G, and the lack of depolarization is consistent with the weakly magnetised gas in DLAs being non-turbulent and quiescent. However, we find mild suggestive evidence that LLSs have coherent magnetic fields: after controlling for the redshift-distribution of our data, we find a 71.5% probability that LLSs have a higher RM than a control sample. We also find strong evidence that LLSs host random magnetic fields, with a 95.5% probability that LLS lines-of-sight have lower polarized fractions than a control sample. The regular coherent magnetic fields within the LLSs must be $\le2.4$ $\mu$G, and the magnetised gas must be highly turbulent with a typical scale on the order of $\approx5$-20 pc, which is similar to that of the Milky Way. This is consistent with the standard dynamo pedagogy, whereby magnetic fields in protogalaxies increase in coherence and strength as a function of cosmic time. Our results are consistent with a hierarchical galaxy formation scenario, with the DLAs, LLSs, and strong magnesium II (MgII) systems exploring three different stages of magnetic field evolution in galaxies.

In a previous paper, we developed tools for studying the horizon geometry of a Kerr black hole that is tidally distorted by a binary companion using techniques that require large mass ratios but can be applied to any bound orbit and allow for arbitrary black hole spin. We now apply these tools to generic Kerr black hole orbits. This allows us to investigate horizon dynamics: the tidal field perturbing the horizon's geometry varies over a generic orbit, with significant variations for eccentric orbits. Many of the features of the horizon's behavior found previously carry over to the dynamical case in a natural way. In particular, we find significant offsets between the applied tide and the horizon's response. This leads to bulging in the horizon's geometry which can lag or lead the orbit, depending upon the hole's rotation and the orbit's geometry. An interesting and apparently new feature we find are small-amplitude, high-frequency oscillations in the horizon's response. We have not been able to identify a mechanism for producing these oscillations, but find that they appear most clearly when rapidly rotating black holes are distorted by very strong-field orbits.

There is evidence that magnetized material along the line of sight to distant quasars is detectable in the polarization properties of the background sources. The polarization properties appear to be correlated with the presence of intervening MgII absorption, which is thought to arise in outflowing material from star forming galaxies. In order to investigate this further, we have obtained high spectral resolution polarization measurements, with the VLA and ATCA, of a set of 49 unresolved quasars for which we have high quality optical spectra. These enable us to produce a Faraday Depth spectrum for each source, using Rotation Measure Synthesis. Our new independent radio data confirms that interveners are strongly associated with depolarization. We characterize the complexity of the Faraday Depth spectrum using a number of parameters and show how these are related, or not, to the depolarization and to the presence of MgII absorption along the line of sight. We argue that complexity and structure in the Faraday Depth distribution likely arise from both intervening material and intrinsically to the background source and attempt to separate these. We find that the strong radio depolarization effects associated with intervening material at redshifts out to $z \approx 1$ arise from inhomogeneous Faraday screens producing a dispersion in Rotation Measure across individual sources of around 10~rad/m$^2$. This is likely produced by disordered fields with strengths of at least $3\;\mu$G.

The W50/SS433 system is an unusual Galactic outflow-driven object of debatable origin. We have used the Australia Telescope Compact Array (ATCA) to observe a new 198 pointing mosaic, covering $3^\circ \times 2^\circ$, and present the highest-sensitivity full-Stokes data of W50 to date using wide-field, wide-band imaging over a 2 GHz bandwidth centred at 2.1 GHz. We also present a complementary H$\alpha$ mosaic created using the Isaac Newton Telescope Photometric H$\alpha$ Survey of the Northern Galactic Plane (IPHAS). The magnetic structure of W50 is found to be consistent with the prevailing hypothesis that the nebula is a reanimated shell-like supernova remnant (SNR), that has been re-energised by the jets from SS433. We observe strong depolarization effects that correlate with diffuse H$\alpha$ emission, likely due to spatially-varying Faraday rotation measure (RM) fluctuations of $\ge48$ to 61 rad m$^{-2}$ on scales $\le4.5$ to 6 pc. We also report the discovery of numerous, faint, H$\alpha$ filaments that are unambiguously associated with the central region of W50. These thin filaments are suggestive of a SNR's shock emission, and almost all have a radio counterpart. Furthermore, an RM-gradient is detected across the central region of W50, which we interpret as a loop magnetic field with a symmetry axis offset by $\approx90^{\circ}$ to the east-west jet-alignment axis, and implying that the evolutionary processes of both the jets and the SNR must be coupled. A separate RM-gradient is associated with the termination shock in the Eastern ear, which we interpret as a ring-like field located where the shock of the jet interacts with the circumstellar medium. Future optical observations will be able to use the new H$\alpha$ filaments to probe the kinematics of the shell of W50, potentially allowing for a definitive experiment on W50's formation history.

We present Australia Telescope Compact Array radio data of the dwarf irregular galaxy ESO 324-G024 which is seen in projection against the giant, northern lobe of the radio galaxy Centaurus A (Cen A, NGC 5128). The distorted morphology and kinematics of ESO 324-G024, as observed in the 21 cm spectral line emission of neutral hydrogen, indicate disruptions by external forces. We investigate whether tidal interactions and/or ram pressure stripping are responsible for the formation of the HI tail stretching to the northeast of ESO 324-G024 with the latter being most probable. Furthermore, we closely analyze the sub-structure of Cen A's polarized radio lobes to ascertain whether ESO 324-G024 is located in front, within or behind the northern lobe. Our multi-wavelength, multi-component approach allows us to determine that ESO 324-G024 is most likely behind the northern radio lobe of Cen A. This result helps to constrain the orientation of the lobe, which is likely inclined to our line of sight by approximately 60 degrees if NGC 5128 and ESO 324-G024 are at the same distance.

We use the integrated polarized radio emission at 1.4 GHz ($\Pi_{\rm 1.4\,GHz}$) from a large sample of AGN (796 sources at redshifts $z<0.7$) to study the large-scale magnetic field properties of radio galaxies in relation to the host galaxy accretion state. We find a fundamental difference in $\Pi_{\rm 1.4\,GHz}$ between radiative-mode AGN (i.e. high-excitation radio galaxies, HERGs, and radio-loud QSOs) and jet-mode AGN (i.e. low-excitation radio galaxies, LERGs). While LERGs can achieve a wide range of $\Pi_{\rm 1.4\,GHz}$ (up to $\sim$$30\%$), the HERGs and radio-loud QSOs are limited to $\Pi_{\rm 1.4\,GHz} \lesssim 15\%$. A difference in $\Pi_{\rm 1.4\,GHz}$ is also seen when the sample is divided at 0.5% of the total Eddington-scaled accretion rate, where the weakly accreting sources can attain higher values of $\Pi_{\rm 1.4\,GHz}$. We do not find any clear evidence that this is driven by intrinsic magnetic field differences of the different radio morphological classes. Instead, we attribute the differences in $\Pi_{\rm 1.4\,GHz}$ to the local environments of the radio sources, in terms of both the ambient gas density and the magnetoionic properties of this gas. Thus, not only are different large-scale gaseous environments potentially responsible for the different accretion states of HERGs and LERGs, we argue that the large-scale magnetised environments may also be important for the formation of powerful AGN jets. Upcoming high angular resolution and broadband radio polarization surveys will provide the high precision Faraday rotation measure and depolarization data required to robustly test this claim.

Faraday rotation of polarised background sources is a unique probe of astrophysical magnetic fields in a diverse range of foreground objects. However, to understand the properties of the polarised sources themselves and of depolarising phenomena along the line of sight, we need to complement Faraday rotation data with polarisation observations over very broad bandwidths. Just as it is impossible to properly image a complex source with limited u-v coverage, we can only meaningfully understand the magneto-ionic properties of polarised sources if we have excellent coverage in $\lambda^2$-space. We here propose a set of broadband polarisation surveys with the Square Kilometre Array, which will provide a singular set of scientific insights on the ways in which galaxies and their environments have evolved over cosmic time.

Relativistic jets in active galactic nuclei (AGN) are among the most powerful astrophysical objects discovered to date. Indeed, jetted AGN studies have been considered a prominent science case for SKA, and were included in several different chapters of the previous SKA Science Book (Carilli & Rawlings 2004). Most of the fundamental questions about the physics of relativistic jets still remain unanswered, and await high-sensitivity radio instruments such as SKA to solve them. These questions will be addressed specially through analysis of the massive data sets arising from the deep, all-sky surveys (both total and polarimetric flux) from SKA1. Wide-field very-long-baseline-interferometric survey observations involving SKA1 will serve as a unique tool for distinguishing between extragalactic relativistic jets and star forming galaxies via brightness temperature measurements. Subsequent SKA1 studies of relativistic jets at different resolutions will allow for unprecedented cosmological studies of AGN jets up to the epoch of re-ionization, enabling detailed characterization of the jet composition, magnetic field, particle populations, and plasma properties on all scales. SKA will enable us to study the dependence of jet power and star formation on other properties of the AGN system. SKA1 will enable such studies for large samples of jets, while VLBI observations involving SKA1 will provide the sensitivity for pc-scale imaging, and SKA2 (with its extraordinary sensitivity and dynamic range) will allow us for the first time to resolve and model the weakest radio structures in the most powerful radio-loud AGN.

Tidal coupling between members of a compact binary system can have an interesting and important influence on that binary's dynamical inspiral. Tidal coupling also distorts the binary's members, changing them (at lowest order) from spheres to ellipsoids. At least in the limit of fluid bodies and Newtonian gravity, there are simple connections between the geometry of the distorted ellipsoid and the impact of tides on the orbit's evolution. In this paper, we develop tools for investigating tidal distortions of rapidly rotating black holes using techniques that are good for strong-field, fast-motion binary orbits. We use black hole perturbation theory, so our results assume extreme mass ratios. We develop tools to compute the distortion to a black hole's curvature for any spin parameter, and for tidal fields arising from any bound orbit, in the frequency domain. We also develop tools to visualize the horizon's distortion for black hole spin $a/M \le \sqrt{3}/2$ (leaving the more complicated $a/M > \sqrt{3}/2$ case to a future analysis). We then study how a Kerr black hole's event horizon is distorted by a small body in a circular, equatorial orbit. We find that the connection between the geometry of tidal distortion and the orbit's evolution is not as simple as in the Newtonian limit.

(abridged) We run a Faraday structure determination data challenge to benchmark the currently available algorithms including Faraday synthesis (previously called RM synthesis in the literature), wavelet, compressive sampling and $QU$-fitting. The frequency set is similar to POSSUM/GALFACTS with a 300 MHz bandwidth from 1.1 to 1.4 GHz. We define three figures of merit motivated by the underlying science: a) an average RM weighted by polarized intensity, RMwtd, b) the separation $\Delta\phi$ of two Faraday components and c) the reduced chi-squared. Based on the current test data of signal to noise ratio of about 32, we find that: (1) When only one Faraday thin component is present, most methods perform as expected, with occasional failures where two components are incorrectly found; (2) For two Faraday thin components, QU-fitting routines perform the best, with errors close to the theoretical ones for RMwtd, but with significantly higher errors for $\Delta\phi$. All other methods including standard Faraday synthesis frequently identify only one component when $\Delta\phi$ is below or near the width of the Faraday point spread function; (3) No methods, as currently implemented, work well for Faraday thick components due to the narrow bandwidth; (4) There exist combinations of two Faraday components which produce a large range of acceptable fits and hence large uncertainties in the derived single RMs; in these cases, different RMs lead to the same Q, U behavior, so no method can recover a unique input model.

Strong singly-ionized magnesium (MgII) absorption lines in quasar spectra typically serve as a proxy for intervening galaxies along the line of sight. Previous studies have found a correlation between the number of these MgII absorbers and the Faraday rotation measure (RM) at $\approx5$ GHz. We cross-match a sample of 35,752 optically-identified non-intrinsic MgII absorption systems with 25,649 polarized background radio sources for which we have measurements of both the spectral index and RM at 1.4 GHz. We use the spectral index to split the resulting sample of 599 sources into flat-spectrum and steep-spectrum subsamples. We find that our flat-spectrum sample shows significant ($\sim3.5\sigma$) evidence for a correlation between MgII absorption and RM at 1.4 GHz, while our steep-spectrum sample shows no such correlation. We argue that such an effect cannot be explained by either luminosity or other observational effects, by evolution in another confounding variable, by wavelength-dependent polarization structure in an active galactic nucleus, by the Galactic foreground, by cosmological expansion, or by partial coverage models. We conclude that our data are most consistent with intervenors directly contributing to the Faraday rotation along the line of sight, and that the intervening systems must therefore have coherent magnetic fields of substantial strength ($\bar{B}=1.8\pm0.4$ $\mu$G). Nevertheless, the weak nature of the correlation will require future high-resolution and broadband radio observations in order to place it on a much firmer statistical footing.

A VLA Sky Survey of the extragalactic sky at S band (2-4 GHz) with polarization information can uniquely probe the magneto-ionic medium in a wide range of astrophysical environments over cosmic time. For a shallow all-sky survey, we expect to detect over 4 million sources in total intensity $>$ 0.45 mJy beam$^{-1}$ and over 2.2$\times$10$^5$ sources in polarized intensity. With these new observations, we expect to discover new classes of polarized radio sources in very turbulent astrophysical environments and those with extreme values of Faraday depth. Moreover, by determining reliable Faraday depths and by modeling depolarization effects, we can derive properties of the magneto-ionic medium associated with AGNs, absorption line systems and galaxies, addressing the following unresolved questions: (1) What is the covering fraction, the degree of turbulence and the origin of absorption line systems? (2) What is the thermal content in AGNs and radio galaxies? (3) How do AGNs and galaxies evolve over cosmic time? (4) What causes the increase in percentage polarization with decreasing flux densities at the low flux density end of the polarized source count? (5) What is the growth rate of large-scale magnetic fields in galaxies?

We present full-Stokes radio polarization observations of the quasar PKS B2126-158 (z=3.268) from 1 to 10 GHz using the Australia Telescope Compact Array. The source has large fractional circular polarization, m_c \equiv |V|/I, detected at high significance across the entire band (from 15 to 90\sigma per 128 MHz sub-band). This allows us to construct the most robust circular polarization (CP) spectrum of an AGN jet to date. We find m_c \propto \nu^{+0.60\pm0.03} from 1.5 to 6.5 GHz, with a peak of m_c ~ 1% before the spectrum turns over somewhere between 6.5 and 8 GHz, above which m_c \propto \nu^{-3.0\pm0.4}. The fractional linear polarization (p) varies from <~0.2% to ~1% across our frequency range and is strongly anti-correlated with the fractional CP, with a best-fit power law giving m_c \propto p^{-0.24\pm0.03}. This is the first clear relation between the observed linear and circular polarisations of an AGN jet, revealing the action of Faraday conversion of linear polarization (LP) to CP within the jet. More detailed modelling in conjunction with high-spatial resolution observations are required to determine the true driving force behind the conversion (i.e. magnetic twist or internal Faraday rotation). In particular determining whether the observed Faraday rotation is internal or entirely external to the jet is key to this goal. The simplest interpretation of our observations favours some internal Faraday rotation, implying that Faraday rotation-driven conversion of LP to CP is the dominant CP generation mechanism. In this case, a small amount of vector-ordered magnetic field along the jet axis is required, along with internal Faraday rotation from the low energy end of the relativistic electron energy spectrum in an electron-proton dominated jet.

[abridged] We report on Suzaku observations of selected regions within the Southern giant lobe of the radio galaxy Centaurus A. We focus on distinct X-ray features likely associated with fine radio structure of the lobe. We find that the spectral properties of the detected X-ray features are equally consistent with thermal emission from hot gas, or with a power-law radiation continuum. However, the plasma parameters implied by these different models favor a synchrotron origin for the analyzed X-ray spots, indicating that a very efficient acceleration of electrons is taking place within the giant structure of Centaurus A, albeit only in isolated and compact regions. We also present a detailed analysis of the diffuse X-ray emission, resulting in a tentative detection of a soft excess component best fitted by a thermal model with a temperature of 0.5 keV. The exact origin of the observed excess remains uncertain, although energetic considerations point to thermal gas filling the bulk of the volume of the lobe and mixed with the non-thermal plasma. The corresponding pressure of the thermal gas in such a case appears to be in almost exact equipartition with the pressure provided by the radio-emitting electrons and the magnetic field. Although tentative, our findings potentially imply that the structure of the extended lobes in active galaxies is likely to be highly inhomogeneous, with magnetic reconnection processes continuously converting magnetic energy to internal energy of the plasma particles, leading to spatial and temporal variations in the plasma parameters around the equilibrium condition.

We present a Faraday rotation measure (RM) study of the diffuse, polarized, radio emission from the giant lobes of the nearest radio galaxy, Centaurus A. After removal of the smooth Galactic foreground RM component, using an ensemble of background source RMs located outside the giant lobes, we are left with a residual RM signal associated with the giant lobes. We find the most likely origin of this residual RM is from thermal material mixed throughout the relativistic lobe plasma. The alternative possibility of a thin-skin/boundary layer of magnetoionic material swept up by the expansion of the lobes is highly unlikely since it requires, at least, an order of magnitude enhancement of the swept up gas over the expected intragroup density on these scales. Strong depolarisation observed from 2.3 to 0.96 GHz also supports the presence of a significant amount of thermal gas within the lobes; although depolarisation solely due to RM fluctuations in a foreground Faraday screen on scales smaller than the beam cannot be ruled out. Considering the internal Faraday rotation scenario, we find a thermal gas number density of ~10^{-4} cm^{-3} implying a total gas mass of ~10^{10} M_solar within the lobes. The thermal pressure associated with this gas (with temperature kT ~ 0.5 keV, obtained from recent X-ray results) is approximately equal to the non-thermal pressure, indicating that over the volume of the lobes, there is approximate equipartition between the thermal gas, radio-emitting electrons and magnetic field (and potentially any relativistic protons present).

Centaurus B is a nearby radio galaxy positioned in the Southern hemisphere close to the Galactic plane. Here we present a detailed analysis of about 43 months of accumulated Fermi-LAT data of the gamma-ray counterpart of the source initially reported in the 2nd Fermi-LAT catalog, and of newly acquired Suzaku X-ray data. We confirm its detection at GeV photon energies, and analyze the extension and variability of the gamma-ray source in the LAT dataset, in which it appears as a steady gamma-ray emitter. The X-ray core of Centaurus B is detected as a bright source of a continuum radiation. We do not detect however any diffuse X-ray emission from the known radio lobes, with the provided upper limit only marginally consistent with the previously claimed ASCA flux. Two scenarios that connect the X-ray and gamma-ray properties are considered. In the first one, we assume that the diffuse non-thermal X-ray emission component is not significantly below the derived Suzaku upper limit. In this case, modeling the inverse-Compton emission shows that the observed gamma-ray flux of the source may in principle be produced within the lobes. This association would imply that efficient in-situ acceleration of the radiating electrons is occurring and that the lobes are dominated by the pressure from the relativistic particles. In the second scenario, with the diffuse X-ray emission well below the Suzaku upper limits, the lobes in the system are instead dominated by the magnetic pressure. In this case, the observed gamma-ray flux is not likely to be produced within the lobes, but instead within the nuclear parts of the jet. By means of synchrotron self-Compton modeling we show that this possibility could be consistent with the broad-band data collected for the unresolved core of Centaurus B, including the newly derived Suzaku spectrum.

We aim to summarize the current state of knowledge regarding Galactic Faraday rotation in an all-sky map of the Galactic Faraday depth. For this we have assembled the most extensive catalog of Faraday rotation data of compact extragalactic polarized radio sources to date. In the map making procedure we use a recently developed algorithm that reconstructs the map and the power spectrum of a statistically isotropic and homogeneous field while taking into account uncertainties in the noise statistics. This procedure is able to identify some rotation angles that are offset by an integer multiple of pi. The resulting map can be seen as an improved version of earlier such maps and is made publicly available, along with a map of its uncertainty. For the angular power spectrum we find a power law behavior with a power law index of -2.14 for a Faraday sky where an overall variance profile as a function of Galactic latitude has been removed, in agreement with earlier work. We show that this is in accordance with a 3D Fourier power spectrum P(k) proportional to k^-2.14 of the underlying field n_e times B_r under simplifying geometrical and statistical assumptions.

We present a detailed study of the Faraday depth structure of four bright (> 1 Jy), strongly polarized, unresolved, radio-loud quasars. The Australia Telescope Compact Array (ATCA) was used to observe these sources with 2 GHz of instantaneous bandwidth from 1.1 to 3.1 GHz. This allowed us to spectrally resolve the polarization structure of spatially unresolved radio sources, and by fitting various Faraday rotation models to the data, we conclusively demonstrate that two of the sources cannot be described by a simple rotation measure (RM) component modified by depolarization from a foreground Faraday screen. Our results have important implications for using background extragalactic radio sources as probes of the Galactic and intergalactic magneto-ionic media as we show how RM estimations from narrow-bandwidth observations can give erroneous results in the presence of multiple interfering Faraday components. We postulate that the additional RM components arise from polarized structure in the compact inner regions of the radio source itself and not from polarized emission from Galactic or intergalactic foreground regions. We further suggest that this may contribute significantly to any RM time-variability seen in RM studies on these angular scales. Follow-up, high-sensitivity VLBI observations of these sources will directly test our predictions.

Global VLBI (EVN+VLBA) polarization observations at 5 and 8.4 GHz of ten high redshift (z > 3) quasars are presented. The core and jet brightness temperatures are found through modelling the self-calibrated uv-data with Gaussian components, which provide reliable estimates of the flux density and size of individual components. The observed high core brightness temperatures (median $T_{\rm b,\,core}=4\times10^{11}$ K) are consistent with Doppler boosted emission from a relativistic jet orientated close to the line-of-sight. This can also explain the dramatic jet bends observed for some of our sources since small intrinsic bends can be significantly amplified due to projection effects in a highly beamed relativistic jet. We also model-fit the polarized emission and, by taking the minimum angle separation between the model-fitted polarization angles at 5 and 8.4 GHz, we calculate the minimum inferred Faraday rotation measure (RM$_{\rm min}$) for each component. We also calculate the minimum intrinsic RM in the rest frame of the AGN (RM$_{\rm min}^{\rm intr}$ = RM$_{\rm min} (1+z)^2$), first subtracting the integrated (presumed foreground) RM in those cases where we felt we could do this reliably. The resulting mean core $|$RM$_{\rm min}^{\rm intr}|$ is 5580 rad m$^{-2}$, with a standard deviation of 3390 rad m$^{-2}$, for four high-z quasars for which we believe we could reliably remove the foreground RM. We find relatively steep core and jet spectral index values, with a median core spectral index of -0.3 and a median jet spectral index of -1.0. Comparing our results with RM observations of more nearby Active Galactic Nuclei at similar emitted frequencies does not provide any significant evidence for dependence of the quasar nuclear environment with redshift.

Using the effective-one-body (EOB) formalism and a time-domain Teukolsky code, we generate inspiral, merger, and ringdown waveforms in the small mass-ratio limit. We use EOB inspiral and plunge trajectories to build the Teukolsky equation source term, and compute full coalescence waveforms for a range of black hole spins. By comparing EOB waveforms that were recently developed for comparable mass binary black holes to these Teukolsky waveforms, we improve the EOB model for the (2,2), (2,1), (3,3), and (4,4) modes. Our results can be used to quickly and accurately extract useful information about merger waves for binaries with spin, and should be useful for improving analytic models of such binaries.

A 45 deg2 radio continuum imaging campaign of the nearest radio galaxy, Centaurus A, is reported. Using the Australia Telescope Compact Array and the Parkes 64m radio telescope at 1.4 GHz, the spatial resolution of the resultant image is ~600 pc, resolving the ~500 kpc giant radio lobes with approximately five times better physical resolution compared to any previous image, and making this the most detailed radio continuum image of any radio galaxy to date. In this paper, we present these new data and discuss briefly some of the most interesting morphological features that we have discovered in the images... abstract truncated.

It is generally believed that turbulence has a significant impact on the dynamics and evolution of molecular clouds and the star formation which occurs within them. Non-ideal magnetohydrodynamic effects are known to influence the nature of this turbulence. We present the results of a suite of 512-cubed resolution simulations of the decay of initially super-Alfvenic and supersonic fully multifluid MHD turbulence. We find that ambipolar diffusion increases the rate of decay of the turbulence while the Hall effect has virtually no impact. The decay of the kinetic energy can be fitted as a power-law in time and the exponent is found to be -1.34 for fully multifluid MHD turbulence. The power spectra of density, velocity and magnetic field are all steepened significantly by the inclusion of non-ideal terms. The dominant reason for this steepening is ambipolar diffusion with the Hall effect again playing a minimal role except at short length scales where it creates extra structure in the magnetic field. Interestingly we find that, at least at these resolutions, the majority of the physics of multifluid turbulence can be captured by simply introducing fixed (in time and space) resistive terms into the induction equation without the need for a full multifluid MHD treatment. The velocity dispersion is also examined and, in common with previously published results, it is found not to be power-law in nature.

We use observations of six "blazars" with the Very Long Baseline Array (VLBA), at eight frequencies (4.6, 5.1, 7.9, 8.9, 12.9, 15.4, 22.2, 43.1 GHz), to investigate the frequency-dependent position of their VLBI cores ("core-shift") and their overall jet spectral distribution. By cross-correlating the optically thin jet emission, we are able to accurately align the multi-frequency images of three of the jets (1418+546, 2007+777, 2200+420), whose core-shifts and spectra we find consistent with the equipartition regime of the Blandford & Konigl conical jet model, where the position of the radio core from the base of the jet follows $r_{core}\propto\nu^{-1}$. For the jet of 0954+658, we align the higher frequency images using our lower frequency measurements assuming equipartition in the radio core from 4.6-43 GHz. The jet emission of the other two sources in our sample (1156+295, 1749+096) is too sparse for our alignment technique to work. Using our measured core-shifts, we calculate equipartition magnetic field strengths of the order of 10's to 100's of mG in the radio cores of these four AGN from 4.6-43 GHz. Extrapolating our results back to the accretion disk and black hole jet-launching distances, we find magnetic field strengths consistent with those expected from theoretical models of magnetically powered jets.

It is well known that non-ideal magnetohydrodynamic effects are important in the dynamics of molecular clouds: both ambipolar diffusion and possibly the Hall effect have been identified as significant. We present the results of a suite of simulations with a resolution of 512-cubed of turbulent decay in molecular clouds incorporating a simplified form of both ambipolar diffusion and the Hall effect simultaneously. The initial velocity field in the turbulence is varied from being super-Alfv\'enic and hypersonic, through to trans-Alfv\'enic but still supersonic. We find that ambipolar diffusion increases the rate of decay of the turbulence increasing the decay from $t^{-1.25}$ to $t^{-1.4}$. The Hall effect has virtually no impact in this regard. The power spectra of density, velocity and the magnetic field are all affected by the non-ideal terms, being steepened significantly when compared with ideal MHD turbulence with exponents. The density power spectra components change from about 1.4 to about 2.1 for the ideal and non-ideal simulations respectively, and power spectra of the other variables all show similar modifications when non-ideal effects are considered. Again, the dominant source of these changes is ambipolar diffusion rather than the Hall effect. There is also a decoupling between the velocity field and the magnetic field at short length scales. The Hall effect leads to enhanced magnetic reconnection, and hence less power, at short length scales. The dependence of the velocity dispersion on the characteristic length scale is studied and found not to be power-law in nature.

We investigate the acceleration of particles via the second-order Fermi process in the lobes of giant radio galaxies. Such sites are candidates for the accelerators of ultra-high energy cosmic rays. We focus on the nearby FR I radio galaxy Centaurus A. This is motivated by the coincidence of its position with the arrival direction of several of the highest energy Auger events. The conditions necessary for consistency with the acceleration timescales predicted by quasi-linear theory are reviewed. Test particle calculations are performed in fields which guarantee electric fields with no component parallel to the local magnetic field. The results of quasilinear theory are, to order of magnitude, found to be accurate at low turbulence levels for non-relativistic Alfven waves and at both low and high turbulence levels in the mildly relativistic case. We conclude that for pure stochastic acceleration to be plausible as the generator of ultra-high energy cosmic rays in Centaurus A, the baryon number density would need to be several orders of magnitude below currently held upper-limits.

We present results from simulations of the extragalactic polarized sky at 1.4 GHz. As the basis for our polarization models, we use a semi-empirical simulation of the extragalactic total intensity (Stokes I) continuum sky developed at the University of Oxford (http://scubed.physics.ox.ac.uk) under the European SKA Design Study (SKADS) initiative, and polarization distributions derived from analysis of polarization observations. By considering a luminosity dependence for the polarization of AGN, we are able to fit the 1.4 GHz polarized source counts derived from the NVSS and the DRAO ELAIS N1 deep field survey down to approximately 1 mJy. This trend is confirmed by analysis of the polarization of a complete sample of bright AGN. We are unable to fit the additional flattening of the polarized source counts from the deepest observations of the ELAIS N1 survey, which go down to ~0.5 mJy. Below 1 mJy in Stokes I at 1.4 GHz, starforming galaxies become an increasingly important fraction of all radio sources. We use a spiral galaxy integrated polarization model to make realistic predictions of the number of polarized sources at microJy levels in polarized flux density and hence, realistic predictions of what the next generation radio telescopes such as ASKAP, other SKA pathfinders and the SKA itself will see.

The parsec-scale Faraday rotation measure (RM) distribution of six "blazars" is investigated using multi-frequency (4.6--43 GHz) polarization observations taken on 2006 July 2 with the VLBA. Analysis of the RM provides the direction of the line-of-sight (LoS) magnetic field component, as well as the intrinsic 2-D polarization distribution on the plane of the sky. Our results show that the magnitude of the core RM increases systematically with frequency, and is well described by a power-law, where |RM_{core}| \propto \nu^a. Our measured values of $a$ vary from 0.9 to 3.8, providing information on the assumed power-law fall-off in the electron density with distance from the central engine for each source. RM gradients were detected across the jets of three sources, supporting the presence of helical magnetic fields in a sheath or boundary layer surrounding their jets. We find a bi-modal distribution of the intrinsic jet polarization orientation; either aligned or orthogonal to the jet direction. A helical magnetic field geometry can neatly explain both the bi-model distribution of the jet polarization orientation and the ordered polarization structure detected on these scales. In half the sources, we find that the core RM changes sign with distance from the central engine. We provide an explanation for this by considering a boundary layer of Faraday rotating material threaded by a helical magnetic field, where bends in the relativistic jet or accelerating/decelerating flows give rise to changes in the dominant LoS components of the magnetic field, which in turn gives rise to different signs of the RM. (abridged)

Accurate measurement of the frequency-dependent shift of the self-absorbed radio core is required for multi-frequency analysis of VLBI data since absolute positional information is lost as a result of phase self-calibration. We use the cross-correlation technique of Croke & Gabuzda (2008) on the optically thin jet emission to align our VLBA images. Our results are consistent with those obtained from the phase-referencing method, as well as alignment by model-fitted optically thin jet components. Physical parameters of the compact jet regions, such as the magnetic field strength (B_core) and the distance of the radio core to the jet origin (r_core), can be calculated from these measurements. For the source Mrk 501, we find a magnetic field strength of 0.15\pm0.04 G in the 8.4-GHz core at a distance of 0.8\pm0.2 pc from the base of the jet. By extrapolating our 4.6 to 15.4 GHz results for BL Lac (2200+420), we estimate magnetic field strengths of the order of 1 G in the millimetre VLBI core. Using our core-shift measurement between 1.6 and 4.8 GHz for 1803+784, we find B_core(4.8 GHz) = 0.11\pm0.02 G and r_core(4.8 GHz) = 20\pm5 pc. The phase-referencing observations of this source at 8.4 and 43 GHz by Jim\'enez-Monferrer et al. (2008) imply B_core(43 GHz) = 1.0\pm0.4 G and r_core(43 GHz) = 2.0\pm0.9 pc.

We present the results of parsec-scale circular polarization (CP) measurements based on VLBA data for a number of radio-bright, core-dominated AGN obtained simultaneously at 15, 22 and 43 GHz. The degrees of CP $m_c$ for the VLBI core region at 15 GHz are similar to values reported earlier at this wavelength, with typical values of a few tenths of a percent.We find that $m_c$ as often rises as falls with increasing frequency between 15 and 22 GHz, while $m_c$ at 43 GHz is in all cases higher than at 22 and 15 GHz. This behaviour seems contrary to expectations, since the degree of CP from both synchrotron radiation and Faraday conversion of linear to circular polarization should {\em decrease} towards higher frequencies if the source is homogeneous. The increase in $m_c$ at 43 GHz may be due to the presence of regions of both positive and negative CP with different frequency dependences on small scales within the core region; alternatively, it may be associated with the intrinsic inhomogeneity of a Blandford-K\"onigl-like jet. In several objects, the detected CO appears to be near, but not coincident with the core, although further observations are needed to confirm this. We find several cases of changes in sign with frequency, most often between 22 and 43 GHz. We find tentative evidence for transverse structure in the CP of 1055+018 and 1334$-$127 that is consistent with their being generated by either the synchrotron mechanism or Faraday conversion in a helical magnetic field. Our results confirm the earlier finding that the sign of the CP at a given observing frequency is generally consistent across epochs separated by several years or more, suggesting stability of the magnetic field orientation in the innermost jets.

We investigate the so-called non-resonant cosmic-ray streaming instability, first discussed by Bell (2004). The extent to which thermal damping and ion-neutral collisions reduce the growth of this instability is calculated. Limits on the growth of the non-resonant mode in SN1006 and RX J1713.7-3946 are presented.

The process of diffusive shock acceleration relies on the efficacy with which hydromagnetic waves can scatter charged particles in the precursor of a shock. The growth of self-generated waves is driven by both resonant and non-resonant processes. We perform high-resolution magnetohydrodynamic simulations of the non-resonant cosmic-ray driven instability, in which the unstable waves are excited beyond the linear regime. In a snapshot of the resultant field, particle transport simulations are carried out. The use of a static snapshot of the field is reasonable given that the Larmor period for particles is typically very short relative to the instability growth time. The diffusion rate is found to be close to, or below, the Bohm limit for a range of energies. This provides the first explicit demonstration that self-excited turbulence reduces the diffusion coefficient and has important implications for cosmic ray transport and acceleration in supernova remnants.

We are using "broadband" (4.6 to 43 GHz) multi-frequency VLBA polarization observations of compact AGN to investigate the 3-D structure of their jet magnetic (B) fields. Observing at several frequencies, separated by short and long intervals, enables reliable determination of the distribution of Faraday Rotation, and thereby the intrinsic B field structure. Transverse Rotation Measure (RM) gradients were detected in the jets of 0954+658 and 1418+546, providing evidence for the presence of a helical B field surrounding the jet. The RM in the core regions of 2200+420 (BL Lac), 0954+658 and 1418+546 display different signs in different frequency-intervals (on different spatial scales); we suggest an explanation for this in terms of modest bends in a helical B field surrounding their jets. In future, polarization observations with a combination of VSOP-2 at 8, 22 and 43 GHz and ground arrays at frequencies with corresponding resolution will help map out the distributions of Faraday rotation, spectral index and the 3-D B field structure both across the jet and closer to the central engine, providing strong constraints for any jet B field models.

Multi-frequency (4.6, 5, 5.5, 8, 8.8, 13, 15, 22 & 43 GHz) polarization observations of 6 "blazars" were obtained on the American Very Long Baseline Array (VLBA) over a 24-hr period on 2 July 2006. Observing at several frequencies, separated by short and long intervals, enabled reliable determination of the distribution of Faraday Rotation on a range of scales. In all cases the magnitude of the RM increases in the higher frequency observations, implying that the electron density and/or magnetic field strength is increasing as we get closer to the central engine. After correcting for Faraday rotation, the polarization orientation in the jet is either parallel or perpendicular to the jet direction. A transverse Rotation Measure (RM) gradient was detected in the jet of 0954+658, providing evidence for the presence of a helical magnetic field surrounding the jet. For three of the sources (0954+658, 1418+546, 2200+420), the sign of the RM in the core region changes in different frequency-intervals, indicating that the line-of-sight component of the magnetic field is changing with distance from the base of the jet. We suggest an explanation for this in terms of bends in a relativistic jet surrounded by a helical magnetic field; where there is no clear evidence for pc-scale bends, the same effect can be explained by an accelerating/decelerating jet.

Homan & Lister (2006) have recently published circular-polarization (CP) detections for 34 objects in the MOJAVE sample - a set of bright, compact AGN being monitored by the Very Long Baseline Array at 15 GHz. We report the detection of 15-GHz parsec-scale CP in two more AGN (3C345 and 2231+114), and confirm the MOJAVE detection of CP in 1633+382. It is generally believed that the most likely mechanism for the generation of this CP is Faraday conversion of linear polarization to CP. A helical jet magnetic-field (B-field) geometry can facilitate this process - linearly polarized emission from the far side of the jet is converted to CP as it passes through the magnetised plasma at the front side of the jet on its way toward the observer. In this case, the sign of the generated CP is essentially determined by the pitch angle and helicity of the helical B field. We have determined the pitch-angle regimes and helicities of the helical jet B fields in 8 AGN for which parsec-scale CP has been detected, and used them to predict the expected CP signs for these AGN if the CP is generated via conversion in these helical fields. We have obtained the intriguing result that our predictions agree with the observed signs in all eight cases, provided that the longitudinal B-field components in the jets correspond to South magnetic poles. This clearly non-random pattern demonstrates that the observed CP in AGN is directly associated with the presence of helical jet B fields. These results suggest that helical B fields are ubiquitous in AGN jets.

We investigate potential models that could explain why multiple proto-stellar systems predominantly show single jets. During their formation, stars most frequently produce energetic outflows and jets. However, binary jets have only been observed in a very small number of systems. We model numerically 3D binary jets for various outflow parameters. We also model the propagation of jets from a specific source, namely L1551 IRS 5, known to have two jets, using recent observations as constraints for simulations with a new MHD code. We examine their morphology and dynamics, and produce synthetic emission maps. We find that the two jets interfere up to the stage where one of them is almost destroyed or engulfed into the second one. We are able to reproduce some of the observational features of L1551 such as the bending of the secondary jet. While the effects of orbital motion are negligible over the jets dynamical timeline, their interaction has significant impact on their morphology. If the jets are not strictly parallel, as in most observed cases, we show that the magnetic field can help the collimation and refocusing of both of the two jets.

Preliminary VLBA polarisation results on 6 ``blazars'' from 6.5-cm to 7-mm are presented here. Observing at several different wavelengths, separated by short and long intervals, enabled reliable information about the magnetic (B) field structure to be obtained and for the effect of Faraday Rotation to be determined and corrected. For all sources the magnitude of the core Rotation Measure (RM) derived from the shorter wavelength data was greater than that derived from the longer wavelength data, consistent with a higher electron density and/or B-field strength closer to the central engine. A transverse RM gradient was detected in the jet of 0954+658, providing evidence for the presence of a helical B-field surrounding the jet. The RM in the core region of 2200+420 (BL Lac) displays sign changes in different wavelength intervals (on different spatial scales); we suggest an explanation for this in terms of modest bends in a helical B-field surrounding the jet.

Astrophysical fluids under the influence of magnetic fields are often subjected to single-fluid or two-fluid approximations. In the case of weakly ionized plasmas however, this can be inappropriate due to distinct responses from the multiple constituent species to both collisional and non-collisional forces. As a result, in dense molecular clouds and proto-stellar accretion discs for instance, the conductivity of the plasma may be highly anisotropic leading to phenomena such as Hall and ambipolar diffusion strongly influencing the dynamics. Diffusive processes are known to restrict the stability of conventional numerical schemes which are not implicit in nature. Furthermore, recent work establishes that a large Hall term can impose an additional severe stability limit on standard explicit schemes. Following a previous paper which presented the one-dimensional case, we describe a fully three-dimensional method which relaxes the normal restrictions on explicit schemes for multifluid processes. This is achieved by applying the little known Super TimeStepping technique to the symmetric (ambipolar) component of the evolution operator for the magnetic field in the local plasma rest-frame, and the new Hall Diffusion Scheme to the skew-symmetric (Hall) component.

When modeling astrophysical fluid flows, it is often appropriate to discard the canonical magnetohydrodynamic approximation thereby freeing the magnetic field to diffuse with respect to the bulk velocity field. As a consequence, however, the induction equation can become problematic to solve via standard explicit techniques. In particular, the Hall diffusion term admits fast-moving whistler waves which can impose a vanishing timestep limit. Within an explicit differencing framework, a multifluid scheme for weakly ionised plasmas is presented which relies upon a new approach to integrating the induction equation efficiently. The first component of this approach is a relatively unknown method of accelerating the integration of parabolic systems by enforcing stability over large compound timesteps rather than over each of the constituent substeps. This method, Super Time Stepping, proves to be very effective in applying a part of the Hall term up to a known critical value. The excess of the Hall term above this critical value is then included via a new scheme for pure Hall diffusion.