Faraday rotation measures (RMs) and more general Faraday structures are key parameters for studying cosmic magnetism and also are sensitive probes of faint ionized thermal gas. There is a need to define what derived quantities are required for various scientific studies, and then to address the challenges in determining Faraday structures. A wide variety of algorithms have been proposed to reconstruct these structures. In preparation for the Polarization Sky Survey of the Universe's Magnetism (POSSUM) to be conducted with the Australian Square Kilometre Array Pathfinder (ASKAP) and the ongoing Galactic Arecibo L-band Feeds Array Continuum Transit Survey (GALFACTS), 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 input models include sources with one Faraday thin component, two Faraday thin components and one Faraday thick component. The frequency set is similar to POS-SUM/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, RM wtd , b) the separation ∆φ of two Faraday components and c) the reduced chi-squared χ 2 r . 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 RM wtd , but with significantly higher errors for ∆φ. All other methods including standard Faraday synthesis frequently identify only one component when ∆φ 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. Further exploration of all these issues is required before upcoming surveys will be able to provide reliable results on Faraday structures.

The Faraday dispersion function (FDF), which can be derived from an observed polarization spectrum by Faraday rotation measure synthesis, is a profile of polarized emissions as a function of Faraday depth. We study intrinsic FDFs along sight lines through face-on, Milky-Way-like galaxies by means of a sophisticated galactic model incorporating 3D MHD turbulence, and investigate how much the FDF contains information intrinsically. Since the FDF reflects distributions of thermal and cosmicray electrons as well as magnetic fields, it has been expected that the FDF could be a new probe to examine internal structures of galaxies. We, however, find that an intrinsic FDF along a sight line through a galaxy is very complicated, depending significantly on actual configurations of turbulence. We perform 800 realizations of turbulence, and find no universal shape of the FDF even if we fix the global parameters of the model. We calculate the probability distribution functions of the standard deviation, skewness, and kurtosis of FDFs and compare them for models with different global parameters. Our models predict that the presence of vertical magnetic fields and large scale-height of cosmic-ray electrons tend to make the standard deviation relatively large. Contrastingly, differences in skewness and kurtosis are relatively less significant.

Faraday tomography is a powerful method to diagnose polarizations and Faraday rotations along the line of sight. The quality of Faraday tomography is, however, limited by several conditions. Recently, it is reported that Faraday tomography indicates false signals in some specific situations. In this paper, we systematically investigate the condition of the appearance of false signals in Faraday tomography. We study this by pseudo-observing two sources within a beam, and change in the intrinsic polarization angles, rotation measures, intensities, and frequency coverage. We find that false signals arise when rotation measure between the sources is less than 1.5 times the full width at half maximum of the rotation measure spread function. False signals also depend on the intensity ratio between the sources and are reduced for large ratio. On the other hand, the appearance of false signals does not depend on frequency coverage, meaning that the uncertainty should be correctly understood and taken into consideration even with future wide-band observations such as Square Kilometer Array (SKA).

Determining magnetic field properties in different environments of the cosmic large-scale structure as well as their evolution over redshift is a fundamental step toward uncovering the origin of cosmic magnetic fields. Radio observations permit the study of extragalactic magnetic fields via measurements of the Faraday depth of extragalactic radio sources. Our aim is to investigate how much different extragalactic environments contribute to the Faraday depth variance of these sources. We develop a Bayesian algorithm to distinguish statistically Faraday depth variance contributions intrinsic to the source from those due to the medium between the source and the observer. In our algorithm the Galactic foreground and measurement noise are taken into account as the uncertainty correlations of the Galactic model. Additionally, our algorithm allows for the investigation of possible redshift evolution of the extragalactic contribution. This work presents the derivation of the algorithm and tests performed on mock observations. Because cosmic magnetism is one of the key science projects of the new generation of radio interferometers, we have predicted the performance of our algorithm on mock data collected with these instruments. According to our tests, high-quality catalogs of a few thousands of sources should already enable us to investigate magnetic fields in the cosmic structure.

Faraday tomography allows astronomers to probe the distribution of the magnetic field along the line of sight (LOS), but that can be achieved only after the Faraday spectrum is interpreted. However, the interpretation is not straightforward, mainly because the Faraday spectrum is complicated due to a turbulent magnetic field; it ruins the one-to-one relation between the Faraday depth and the physical depth, and appears as many small-scale features in the Faraday spectrum. In this paper, by employing "simple toy models" for the magnetic field, we describe numerically as well as analytically the characteristic properties of the Faraday spectrum. We show that the Faraday spectrum along "multiple LOSs" can be used to extract the global properties of the magnetic field. Specifically, considering face-on spiral galaxies and modeling turbulent magnetic field as a random field with a single coherence length, we numerically calculate the Faraday spectrum along a number of LOSs and its shape-characterizing parameters, that is, the moments. When multiple LOSs cover a region of (10 coherence length) 2 , the shape of the Faraday spectrum becomes smooth and the shape-characterizing parameters are well specified. With the Faraday spectrum constructed as a sum of Gaussian functions with different means and variances, we analytically show that the parameters are expressed in terms of the regular and turbulent components of the LOS magnetic field and the coherence length. We also consider the turbulent magnetic field modeled with a power-law spectrum, and study how the magnetic field is revealed in the Faraday spectrum. Our work suggests a way to obtain information on the magnetic field from a Faraday tomography study.