Pulsar search with time-domain observation is very computationally expensive and data volume will be enormous with the next generation telescopes such as the Square Kilometre Array. We apply artificial neural networks (ANNs), a machine learning method, for efficient selection of pulsar candidates from radio continuum surveys, which are much cheaper than time-domain observation. With observed quantities such as radio fluxes, sky position and compactness as inputs, our ANNs output the "score" that indicates the degree of likeliness of an object to be a pulsar. We demonstrate ANNs based on existing survey data by the TIFR GMRT Sky Survey (TGSS) and the NRAO VLA Sky Survey (NVSS) and test their performance. Precision, which is the ratio of the number of pulsars classified correctly as pulsars to that of any objects classified as pulsars, is about 96$\%$. Finally, we apply the trained ANNs to unidentified radio sources and our fiducial ANN with five inputs (the galactic longitude and latitude, the TGSS and NVSS fluxes and compactness) generates 2,436 pulsar candidates from 456,866 unidentified radio sources. These candidates need to be confirmed if they are truly pulsars by time-domain observations. More information such as polarization will narrow the candidates down further.

Faraday tomography is a novel method to probe 3-dimensional structure of magnetic fields of polarized radio sources. In this paper, we investigate Faraday dispersion function (FDF) of disk galaxies extending a simple analytic model of galactic magnetic fields developed in Ideguchi et al. (2017). The model consists of axisymmetric coherent fields and turbulent fields and we consider the effects of inclination, relative amplitude of coherent and turbulent magnetic fields and pitch angle of coherent fields. Our simple model makes it easy to obtain physical interpretation of FDFs and helps understanding observational results. We find FDFs have two peaks when galaxies are observed with non zero inclination. The gap and relative height of two peaks are dependent on the inclination angle and pitch angle but are not affected by relative amplitude of coherent and turbulent magnetic fields so much. These findings give us an important caution that two peaks in observed FDFs do not necessarily imply the presence of two separate radio sources within a beam.

Faraday tomography (or rotation measure synthesis) is a procedure to convert linear polarization spectra into the Faraday dispersion function, which provides us with unique information of magneto-ionic media along the line of sight. Mathematical formulation of Faraday tomography is similar to polarimetric imaging of radio interferometry, where many new methods have been actively developed and shown to outperform the standard CLEAN approaches. In this paper, we propose a sparse reconstruction technique to Faraday tomography. This technique is being developed for interferometric imaging and utilizes computationally less expensive convex regularization functions such as $\ell_1$-norm and total variation (TV) or total squared variation (TSV). The proposed technique solves a convex optimization, and therefore its solution is determined uniquely regardless of the initial condition for given regularization parameters that can be optimized by data themselves. Using a physically-motivated model of turbulent galactic magnetized plasma, we demonstrate that the proposed technique outperforms RM-CLEAN and provides higher-fidelity reconstruction. The proposed technique would be a powerful tool in broadband polarimetry with the Square Kilometre Array (SKA) and its precursors.

Faraday tomography is thought to be a powerful tool to explore cosmic magnetic field. Broadband radio polarimetric data is essential to ensure the quality of Faraday tomography, but such data is not easy to obtain because of radio frequency interferences (RFIs). In this paper, we investigate optimum frequency coverage of Faraday tomography so as to explore Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in filaments of galaxies. We adopt a simple model of the IGMF and estimate confidence intervals of the model parameters using the Fisher information matrix. We find that meaningful constraints for RM due to the IGMF are available with data at multiple narrowbands which are scattered over the ultra-high frequency (UHF, 300 MHz - 3000 MHz). The optimum frequency depends on the Faraday thickness of the Milky Way foreground. These results are obtained for a wide brightness range of the background source including fast radio bursts (FRBs). We discuss the relation between the polarized-intensity spectrum and the optimum frequency.

QU-fitting is a standard model-fitting method to reconstruct distribution of magnetic fields and polarized intensity along a line of sight (LOS) from an observed polarization spectrum. In this paper, we examine the performance of QU-fitting by simulating observations of two polarized sources located along the same LOS, varying the widths of the sources and the gap between them in Faraday depth space, systematically. Markov Chain Monte Carlo (MCMC) approach is used to obtain the best-fit parameters for a fitting model, and Akaike and Bayesian Information Criteria (AIC and BIC, respectively) are adopted to select the best model from four fitting models. We find that the combination of MCMC and AIC/BIC works fairly well in model selection and estimation of model parameters in the cases where two sources have relatively small widths and a larger gap in Faraday depth space. On the other hand, when two sources have large width in Faraday depth space, MCMC chain tends to be trapped in a local maximum so that AIC/BIC cannot select a correct model. We discuss the causes and the tendency of the failure of QU-fitting and suggest a way to improve it.

Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography.

Faraday tomography allows astronomers to probe the distribution of magnetic field along the line of sight (LOS), but that can be achieved only after Faraday spectrum is interpreted. However, the interpretation is not straightforward, mainly because Faraday spectrum is complicated due to 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 Faraday spectrum. In this paper, employing "simple toy models" for the magnetic field, we describe numerically as well as analytically the characteristic properties of Faraday spectrum. We show that Faraday spectrum along "multiple loss" can be used to extract the global properties of magnetic field. Specifically, considering face-on spiral galaxies and modeling turbulent magnetic field as a random field with single coherence length, we numerically calculate Faraday spectrum along a number of LOSs and its shape-characterizing parameters, that is, the moments. When multiple LOSs cover a region of $\gtrsim (10\ {\rm coherence\ length)^2}$, the shape of 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 LOS magnetic field and the coherence length. We also consider the turbulent magnetic field modeled with power-law spectrum, and study how the magnetic field is revealed in Faraday spectrum. Our work suggests a way toward obtaining the information of magnetic field from Faraday tomography study.

RM CLEAN is a standard method to reconstruct the distribution of cosmic magnetic fields and polarized sources along the line of sight (LOS) from observed polarization spectrum. This method is similar to the CLEAN algorithm for aperture synthesis radio telescope images but it is rather unclear in what cases RM CLEAN works well. In this paper, we evaluate the performance of RM CLEAN by simulating spectro-polarimetric observations of two compact sources located in the same LOS, varying the relative initial polarization angle and Faraday depth systematically. Especially, we focus on if the two polarized sources can be resolved in the Faraday depth space and how well the source parameters can be estimated. We confirm the previous studies that two sources cannot be resolved when they are closely located in the Faraday depth space for specific values of the relative initial polarization angle. Further, we calculate the chi-square value for the fit between the mock data of polarization spectrum and the one from RM CLEAN. Then we find that the chi-squared value is not always significantly large even when RM CLEAN gives wrong results.

Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides numerous observation points in the cosmic space. Furthermore, the SKA should be designed to facilitate wideband polarimetry so as to allow us to examine sightline structures of magnetic fields by means of depolarization and Faraday Tomography. The SKA will be able to drive cosmic magnetism of the interstellar medium, the Milky Way, galaxies, AGN, galaxy clusters, and potentially the cosmic web which may preserve information of the primeval Universe. The Japan SKA Consortium (SKA-JP) Magnetism Science Working Group (SWG) proposes the project "Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography", which contains ten scientific use cases.

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 the 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. With cosmic magnetism being one of the key science projects of the new generation of radio interferometers we have made predictions for the algorithm's performance on data from the next generation of radio interferometers. Applications to real data are left for future work.

To better understand the origin and properties of cosmological magnetic fields, a detailed knowledge of magnetic fields in the large-scale structure of the Universe (galaxy clusters, filaments) is crucial. We propose a new statistical approach to study magnetic fields on large scales with the rotation measure grid data that will be obtained with the new generation of radio interferometers.

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

The Faraday dispersion function (FDF), which can be derived from an observed polarization spec- trum 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 cosmic- ray 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 pa- rameters. 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).

We investigate the capability of ongoing radio telescopes for probing Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) in the large-scale structure of the universe which is expected to be of order $O(1) {\rm rad/m^2}$. We consider polarization observations of a compact radio source such as quasars behind a diffuse source such as the Galaxy, and calculate Stokes parameters $Q$ and $U$ assuming a simple model of the Faraday dispersion functions with Gaussian shape. Then, we perform the Fisher analysis to estimate the expected errors in the model parameters from QU-fitting of polarization intensity, accounting for sensitivities and frequency bands of Australian Square Kilometer Array Pathfinder, Low Frequency Array, and the Giant Meterwave Radio Telescope. Finally, we examine the condition on the source intensities which are required to detect the IGMF. Our analysis indicates that the QU-fitting is promising for forthcoming wideband polarimetry to explore RM due to the IGMF in filaments of galaxies.