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.

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.

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.

SS433, located at the center of the supernova remnant W50, is a close proximity binary system consisting of a compact star and a normal star. Jets of material are directed outwards from the vicinity of the compact star symmetrically to the east and west. Non-thermal hard X-ray emission is detected from lobes lying on both sides. Shock accelerated electrons are expected to generate sub-TeV gamma rays through the inverse-Compton process in the lobes. Observations of the western X-ray lobe region of SS433/W50 system have been performed to detect sub-TeV gamma-rays using the 10m CANGAROO-II telescope in August and September, 2001, and July and September, 2002. The total observation times are 85.2 hours for ON source, and 80.8 hours for OFF source data. No significant excess of sub-TeV gamma rays has been found at 3 regions of the western X-ray lobe of SS433/W50 system. We have derived 99% confidence level upper limits to the fluxes of gamma rays and have set constraints on the strengths of the magnetic fields assuming the synchrotron/inverse-Compton model for the wide energy range of photon spectrum from radio to TeV. The derived lower limits are 4.3 microgauss for the center of the brightest X-ray emission region and 6.3 microgauss for the far end from SS433 in the western X-ray lobe. In addition, we suggest that the spot-like X-ray emission may provide a major contribution to the hardest X-ray spectrum in the lobe.