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.

We report the results of the Australia Telescope Compact Array (ATCA) 15 mm observation of the Phoenix galaxy cluster possessing an extreme star-burst brightest cluster galaxy (BCG) at the cluster center. We spatially resolved radio emission around the BCG, and found diffuse bipolar and bar-shape structures extending from the active galactic nucleus (AGN) of the BCG. They are likely radio jets/lobes, whose sizes are ~10-20 kpc and locations are aligned with X-ray cavities. If we assume that the radio jets/lobes expand with the sound velocity, their ages are estimated to be ~10 Myr. We also found compact radio emissions near the center and suggest that they are another young bipolar jets with ~1 Myr of age. Moreover, we found extended radio emission surrounding the AGN and discussed the possibility that the component is a product of the cooling flow, by considering synchrotron radiation partially absorbed by molecular clumps, free-free emission from the warm ionized gas, and the spinning dust emission from dusty circum-galactic medium.

We present measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward SPT-CL J2334-4243 (the Phoenix galaxy cluster) at z=0.597 by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. The SZE is imaged at 5" resolution (corresponding to the physical scale of 23kpc/h) within 200kpc/h from the central galaxy with the peak signal-to-noise ratio exceeding 11. Combined with the Chandra X-ray image, the ALMA SZE data further allow for non-parametric deprojection of electron temperature, density, and entropy. Our method can minimize contamination by the central AGN and the X-ray absorbing gas within the cluster, both of which largely affect the X-ray spectrum. We find no significant asymmetry or disturbance in the SZE image within the current measurement errors. The detected SZE signal agrees well with the average pressure profile of local cool-core clusters. Unlike any other known clusters, however, gas temperature drops by at least a factor of 5 toward the center. We identify ~6x10^{11} M_sun cool gas with temperature ~3keV in the inner 20kpc/h. Taken together, our results imply that the gas is indeed cooling efficiently and nearly isobarically down to this radius in the Phoenix cluster.

This article summarizes the work presented at the workshop 'The Power of Faraday Tomography: towards 3D mapping of cosmic magnetic fields', held in Miyazaki, Japan, in Spring 2018. We place the various oral and poster presentations given at the workshop in a broader perspective and present some highlight results from every presenter.

We report Atacama Large Millimeter/submillimeter Array (ALMA) observations of the binary system containing the pulsar PSR B1259-63 orbiting around a Be star LS 2883 after the 2017 periastron passage. We detected radio continuum emission from the binary system in the millimeter/submillimeter wavelengths for the first time. At Band 3 (97 GHz), the flux 84 days after the periastron is almost the same as that 71 days after the periastron. Although the binary system showed intense GeV gamma-ray flares during our observations, the Band 3 flux did not indicate any time correlation with them. The Band 3 fluxes are consistent with an extrapolation of the radio spectrum at lower frequencies. Assuming that it is synchrotron emission, we constrain magnetic fields ($\lesssim 0.6$ G) and the high-energy cutoff of the electrons ($\gamma \gtrsim 360$). The flux at Band 7 (343 GHz) 69 days after the periastron shows a significant excess from the extrapolation of the radio spectrum at lower frequencies. The flux may be associated with the circumstellar disk around the Be star. We also present the results of Australian Telescope Compact Array (ATCA) observations at 94 GHz for the 2014 periastron passage, which show that the radio spectrum was relatively soft when the pulsar passed the disk.

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.

Observational study of galactic magnetic fields is hampered by the fact that the observables only probe various projections of the magnetic fields. Comparison with numerical simulations is helpful to understand the real structures, and observational visualization of numerical data is an important task. In this paper, we investigate 8~GHz radio synchrotron emission from spiral galaxies, using the data of global three-dimensional magneto-hydrodynamic simulations. We assume a frequency independent depolarization in our observational visualization. We find that the appearance of the global magnetic field depends on the viewing angle: a face-on view seemingly has hybrid magnetic field types combining axisymmetric modes with higer order modes; at a viewing angle of $\sim 70\degr$, the galaxy seems to contain a ring-like magnetic field structure; while in edge-on view, only field structure parallel to the disk can be seen. The magnetic vector seen at 8~GHz traces the global magnetic field inside the disk. These results indicate that the topology of global magnetic field obtained from the relation between azimuthal angle and Faraday depth strongly depends on the viewing anglue of the galaxy. As one of the examples, we compare our results at a viewing angle of $25\degr$ with the results of IC342. The relation between azimuthal angle and Faraday depth of the numerical result shows a tendency similar to IC342, such as the peak numbers of the Faraday depth.

Observational study of galactic magnetic fields is limited by projected observables. Comparison with numerical simulations is helpful to understand the real structures, and observational visualization of numerical data is an important task. Machida et al. (2018) have reported Faraday depth maps obtained from numerical simulations. They showed that the relation between azimuthal angle and Faraday depth depends on the inclination angle. In this paper, we investigate 100MHz to 10GHz radio synchrotron emission from spiral galaxies, using the data of global three-dimensional magneto-hydrodynamic simulations. We model internal and external Faraday depolarization at small scales and assume a frequency independent depolarization. It is found that the internal and external Faraday depolarization becomes comparable inside the disk and the dispersion of Faraday depth becomes about 4rad/m^{2} for face-on view and 40rad/m2 for edge-on view, respectively. The internal depolarization becomes ineffective in the halo. Because of the magnetic turbulence inside the disk, frequency independent depolarization works well and the polarization degree becomes 0.3 at high frequency. When the observed frequency is in the 100 MHz band, polarized intensity vanishes in the disk, while that from the halo can be observed. Because the remaining component of polarized intensity is weak in the halo and the polarization degree is about a few %, it may be difficult to observe that component. These results indicate that the structures of global magnetic fields in spiral galaxies could be elucidated, if broadband polarimetry such as that with the Square Kilometre Array is achieved.

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.

RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5{\sigma}. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL).

The spectral slope of the magnetic energy in supernova remnants (SNRs) can be obtained by analysis of spatial two-point correlation functions of synchrotron intensities. This method has been originally developed for the analysis of magnetic field structure in diffuse interstellar medium and applied when the geometry of the emission region is simple and known. In this paper, applying this correlation analysis to Tycho's SNR, for which the synchrotron emission region is known to be a spherical shell, we find that the magnetic energy spectrum shows the Kolmogorov-like scaling. Our results can be explained by turbulence developed downstream of the shock front via Richtmyer-Meshkov instability or an amplification of upstream magnetic field induced by cosmic rays. They could be discriminated by future observations with a sub arcsecond resolution such as Square Kilometer Array.

The Advanced Telescope for High Energy Astrophysics (Athena) is the X-ray observatory large mission selected by the European Space Agency (ESA), within its Cosmic Vision 2015-2025 programme, to address the "Hot and Energetic Universe" scientific theme (Nandra et al. 2013), and it is provisionally due for launch in the early 2030s. The Square Kilometer Array (SKA) is the next generation radio observatory and consists of two telescopes, one comprised of dishes operating at mid frequencies (SKA1-MID) and located in South Africa, and the other comprised of Log-Periodic antennas operating at low radio frequencies (SKA1-LOW), which will be located in Australia (Braun et al. 2017). The scientific commissioning of the radio telescope is planned to begin in 2021-2022. The SKA-Athena Synergy Team (SAST) has been tasked to single out the potential scientific synergies between Athena and SKA. The astrophysical community was involved in this exercise primarily through a dedicated SKA-Athena Synergy Workshop, which took place on April 24-25, 2017 at SKAO, Jodrell Bank, Manchester. The final result of the synergy exercise, this White Paper, describes in detail a number of scientific opportunities that will be opened up by the combination of Athena and SKA, these include: 1. the Cosmic Dawn; 2. the Evolution of black holes and galaxies; 3. Active galaxy feedback in galaxy clusters; 4. Non-thermal phenomena in galaxy clusters; 5. Detecting the cosmic web; 6. Black-hole accretion physics and astrophysical transients; 7. Galactic astronomy: stars, planets, pulsars and supernovae.

We report the Australia Telescope Compact Array 16 cm observation of CIZA J1358.9-4750. Recent X-ray studies imply that this galaxy cluster is composed of merging, binary clusters. Using the EW367 configuration, we found no significant diffuse radio emission in and around the cluster. An upper limit of the total radio power at 1.4 GHz is $\sim 1.1 \times 10^{22}$ Watt/Hz in $30$ square arcminutes which is a typical size of radio relics. It is known that an empirical relation holds between the total radio power and X-ray luminosity of the host cluster. The upper limit is about one order of magnitude lower than the power expected from the relation. Very young ($\sim 70$ Myr) shocks with low Mach numbers ($\sim 1.3$), which are often seen at an early stage of merger simulations, are suggested by the previous X-ray observation. The shocks may generate cosmic-ray electrons with a steep energy spectrum, which is consistent with non-detection of bright ($>10^{23}$ Watt/Hz) relic in this 16 cm band observation. Based on the assumption of energy equipartition, the upper limit gives the magnetic-field strength below $0.68f(D_{\rm los}/{\rm 1~Mpc})^{-1}(\gamma_{\rm min}/200)^{-1}$ $\mu$G, where $f$ is the cosmic-ray total energy density over the cosmic-ray electron energy density, $D_{\rm los}$ is the depth of the shock wave along the sightline and $\gamma_{\rm min}$ is the lower cutoff Lorentz factor of the cosmic-ray electron energy spectrum.

We report a polarization analysis of the eastern region of W50, observed with the Australia Telescope Compact Array (ATCA) at 1.4 - 3.0 GHz. In order to study the physical structures in the region where the SS433 jet and W50 interact, we obtain an intrinsic magnetic field vector map of that region. We find that the orientation of the intrinsic magnetic field vectors are aligned along the total intensity structures, and that there are characteristic, separate structures related to the jet, the bow shock, and the terminal shock. The Faraday rotation measures (RMs), and the results of Faraday Tomography suggest that a high intensity, filamentary structure in the north-south direction of the eastern-edge region can be separated into at least two parts to the north and south. The results of Faraday Tomography also show that there are multiple components along the line of sight and/or within the beam area. In addition, we also analyze the X-ray ring-like structure observed with XMM-Newton. While the possibility still remains that this X-ray ring is real, it seems that the structure is not ring-like at radio wavelengths. Finally, we suggest that the structure is a part of the helical structure that coils the eastern ear of W50.

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.

We report on the results of a $Suzaku$ observation of the plasma in the filament located between the two massive clusters of galaxies Abell 399 and Abell 401. Abell 399 ($z$=0.0724) and Abell 401 ($z$=0.0737) are expected to be in the initial phase of a cluster merger. In the region between the two clusters, we find a clear enhancement in the temperature of the filament plasma from 4 keV (expected value from a typical cluster temperature profile) to $kT\sim$6.5 keV. Our analysis also shows that filament plasma is present out to a radial distance of 15' (1.3 Mpc) from a line connecting the two clusters. The temperature profile is characterized by an almost flat radial shape with $kT\sim$6-7 keV within 10' or $\sim$0.8 Mpc. Across $r$=8'~from the axis, the temperature of the filament plasma shows a drop from 6.3 keV to 5.1 keV, indicating the presence of a shock front. The Mach number based on the temperature drop is estimated to be ${\cal M}\sim$1.3. We also successfully determined the abundance profile up to 15' (1.3 Mpc), showing an almost constant value ($Z$=0.3 solar) at the cluster outskirt. We estimated the Compton $y$-parameter to be $\sim$14.5$\pm1.3\times10^{-6}$, which is in agreement with $Planck$'s results (14-17$\times10^{-6}$ on the filament). The line of sight depth of the filament is $l\sim$1.1 Mpc, indicating that the geometry of filament is likely a pancake shape rather than cylindrical. The total mass of the filamentary structure is $\sim$7.7$\times10^{13}~\rm M_{\odot}$. We discuss a possible interpretation of the drop of X-ray emission at the rim of the filament, which was pushed out by the merging activity and formed by the accretion flow induced by the gravitational force of the filament.

We propose a new method to probe the Warm Hot Intergalactic Medium (WHIM) beyond the virial radius (R_200) of a cluster of galaxies, where X-ray observations are not easily achievable. In this method, we use dispersion measures (DMs) of Fast Radio Bursts (FRBs) that appear behind the cluster and the Sunyaev-Zel'dovich (SZ) effect towards the cluster. The DMs reflect the density of the intracluster medium (ICM) including the WHIM. If we observe a sufficient number of FRBs in the direction of the cluster, we can derive the density profile from the DMs. Similarly, we can derive the pressure profile from the SZ effect. By combining the density and the pressure profiles, the temperature profile can be obtained. Based on mock observations of nearby clusters, we find that the density of the WHIM can be determined even at > 2 R_200 from the cluster center when FRB observations with the Square Kilometre Array (SKA) become available. The temperature can be derived out to r~ 1.5 R_200, and the radius is limited by the current sensitivity of SZ observations.

We present the first image of the thermal Sunyaev-Zel'dovich effect (SZE) obtained by the Atacama Large Millimeter/submillimeter Array (ALMA). Combining 7-m and 12-m arrays in Band 3, we create an SZE map toward a galaxy cluster RXJ1347.5-1145 with 5 arc-second resolution (corresponding to the physical size of 20 kpc/h), the highest angular and physical spatial resolutions achieved to date for imaging the SZE, while retaining extended signals out to 40 arc-seconds. The 1-sigma statistical sensitivity of the image is 0.017 mJy/beam or 0.12 mK_CMB at the 5 arc-second full width at half maximum. The SZE image shows a good agreement with an electron pressure map reconstructed independently from the X-ray data and offers a new probe of the small-scale structure of the intracluster medium. Our results demonstrate that ALMA is a powerful instrument for imaging the SZE in compact galaxy clusters with unprecedented angular resolution and sensitivity. As the first report on the detection of the SZE by ALMA, we present detailed analysis procedures including corrections for the missing flux, to provide guiding methods for analyzing and interpreting future SZE images by ALMA.

We examine the proposal that the dispersion measures (DMs) and Faraday rotation measures (RMs) of extragalactic linearly-polarized fast radio bursts (FRBs) can be used to probe the intergalactic magnetic field (IGMF) in filaments of galaxies. The DM through the cosmic web is dominated by contributions from the warm-hot intergalactic medium (WHIM) in filaments and from the gas in voids. On the other hand, RM is induced mostly by the hot medium in galaxy clusters, and only a fraction of it is produced in the WHIM. We show that if one excludes FRBs whose sightlines pass through galaxy clusters, the line-of-sight strength of the IGMF in filaments, $B_{||}$, is approximately $C(\langle 1+z \rangle/f_{DM})(RM/DM)$, where $C$ is a known constant. Here, the redshift of the FRB is not required to be known; $f_{DM}$ is the fraction of total DM due to the WHIM, while $\langle 1+z \rangle$ is the redshift of interevening gas weighted by the WHIM gas density, both of which can be evaluated for a given cosmology model solely from the DM of an FRB. Using data on structure formation simulations and a model IGMF, we show that $C(\langle 1+z \rangle/f_{DM})(RM/DM)$ closely reproduces the density-weighted line-of-sight strength of the IGMF in filaments of the large-scale structure.

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.

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.

We report polarimetry results of a merging cluster of galaxies Abell 2256 with Karl G. Jansky Very Large Array (JVLA). We performed new observations with JVLA at S-band (2051-3947 MHz) and X-band (8051-9947 MHz) in the C array configuration, and detected significant polarized emissions from the radio relic, Source A, and Source B in this cluster. We calculated the total magnetic field strengths toward the radio relic using revised equipartition formula, which is 1.8-5.0 microG. With dispersions of Faraday rotation measure, magnetic-field strengths toward Sources A and B are estimated to be 0.63-1.26 microG and 0.11-0.21 microG, respectively. An extremely high degree of linear polarization, as high as ~ 35 %, about a half of the maximum polarization, was detected toward the radio relic, which indicates highly ordered magnetic lines of force over the beam sizes (~ 52 kpc).The fractional polarization of the radio relic decreases from ~ 35 % to ~ 20 % around 3 GHz as the frequency decreases and is nearly constant between 1.37 and 3 GHz. Both analyses with depolarization models and Faraday tomography suggest multiple depolarization components toward the radio relic and imply the existence of turbulent magnetic fields.

We know that magnetic fields are pervasive across all scales in the Universe and over all of cosmic time and yet our understanding of many of the properties of magnetic fields is still limited. We do not yet know when, where or how the first magnetic fields in the Universe were formed, nor do we fully understand their role in fundamental processes such as galaxy formation or cosmic ray acceleration or how they influence the evolution of astrophysical objects. The greatest challenge to addressing these issues has been a lack of deep, broad bandwidth polarimetric data over large areas of the sky. The Square Kilometre Array will radically improve this situation via an all-sky polarisation survey that delivers both high quality polarisation imaging in combination with observations of 7-14 million extragalactic rotation measures. Here we summarise how this survey will improve our understanding of a range of astrophysical phenomena on scales from individual Galactic objects to the cosmic web.

CIZA J1358.9-4750 is a nearby (z = 0.074) pair of clusters of galaxies located close to the Galactic plane. It consists of two X-Ray extended humps at north-west and south-east separated by 14 arcmin (~ 1.2 Mpc), and an X-Ray bright bridge-like structure in between. With Suzaku, the south-east hump was shown to have a temperature of 5.6 keV and the north-west one 4.6 keV. Neither humps exhibit significant central cool component. The bridge region has a temperature higher than 9 keV at the maximum, and this hot region is distributed almost orthogonal to the bridge axis in agreement with the shock heating seen in numerical simulations at an early phase of a head-on major merger. This resemblance is supported by good positional coincidence between the X-Ray peaks and cD galaxies associated with each cluster. In a short exposure XMM-Newton image, a significant intensity jump was found at a position where the Suzaku-measured temperature exhibits a steep gradient. These properties indicate the presence of a shock discontinuity. The Mach number is estimated to be 1.32 from the temperature difference across the identified shock front, which gives the colliding velocity of approximately 1800 km/s. From optical redshifts of the member galaxies, the two clusters are indicated to be merging nearly on the sky plane. Thus, CIZA J1358.9-4750 is considered as a valuable nearby example of early-phase merger with a clear shock feature.

(abridged) Observations of Faraday rotation for extragalactic sources probe magnetic fields both inside and outside the Milky Way. Building on our earlier estimate of the Galactic contribution, we set out to estimate the extragalactic contributions. We discuss the problems involved; in particular, we point out that taking the difference between the observed values and the Galactic foreground reconstruction is not a good estimate for the extragalactic contributions. We point out a degeneracy between the contributions to the observed values due to extragalactic magnetic fields and observational noise and comment on the dangers of over-interpreting an estimate without taking into account its uncertainty information. To overcome these difficulties, we develop an extended reconstruction algorithm based on the assumption that the observational uncertainties are accurately described for a subset of the data, which can overcome the degeneracy with the extragalactic contributions. We present a probabilistic derivation of the algorithm and demonstrate its performance using a simulation, yielding a high quality reconstruction of the Galactic Faraday rotation foreground, a precise estimate of the typical extragalactic contribution, and a well-defined probabilistic description of the extragalactic contribution for each data point. We then apply this reconstruction technique to a catalog of Faraday rotation observations. We vary our assumptions about the data, showing that the dispersion of extragalactic contributions to observed Faraday depths is most likely lower than 7 rad/m^2, in agreement with earlier results, and that the extragalactic contribution to an individual data point is poorly constrained by the data in most cases.

Deep surveys with the SKA1-MID array offer for the first time the opportunity to systematically explore the polarization properties of the microJy source population. Our knowledge of the polarized sky approaching these levels is still very limited. In total intensity the population will be dominated by star-forming and normal galaxies to intermediate redshifts ($z \sim1-2$), and low-luminosity AGN to high redshift. The polarized emission from these objects is a powerful probe of their intrinsic magnetic fields and of their magnetic environments. For redshift of order 1 and above the broad bandwidth of the mid-bands span the Faraday thick and thin regimes allowing study of the intrinsic polarization properties of these objects as well as depolarization from embedded and foreground plasmas. The deep field polarization images will provide Rotation Measures data with very high solid angle density allowing a sensitive statistical analysis of the angular variation of RM on critical arc-minute scales from a magnetic component of Large Scale Structure of the Universe.

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.

Magnetic fields in the Milky Way are present on a wide variety of sizes and strengths, influencing many processes in the Galactic ecosystem such as star formation, gas dynamics, jets, and evolution of supernova remnants or pulsar wind nebulae. Observation methods are complex and indirect; the most used of these are a grid of rotation measures of unresolved polarized extragalactic sources, and broadband polarimetry of diffuse emission. Current studies of magnetic fields in the Milky Way reveal a global spiral magnetic field with a significant turbulent component; the limited sample of magnetic field measurements in discrete objects such as supernova remnants and HII regions shows a wide variety in field configurations; a few detections of magnetic fields in Young Stellar Object jets have been published; and the magnetic field structure in the Galactic Center is still under debate. The SKA will unravel the 3D structure and configurations of magnetic fields in the Milky Way on sub-parsec to galaxy scales, including field structure in the Galactic Center. The global configuration of the Milky Way disk magnetic field, probed through pulsar RMs, will resolve controversy about reversals in the Galactic plane. Characteristics of interstellar turbulence can be determined from the grid of background RMs. We expect to learn to understand magnetic field structures in protostellar jets, supernova remnants, and other discrete sources, due to the vast increase in sample sizes possible with the SKA. This knowledge of magnetic fields in the Milky Way will not only be crucial in understanding of the evolution and interaction of Galactic structures, but will also help to define and remove Galactic foregrounds for a multitude of extragalactic and cosmological studies.

We investigate the possibility for the SKA to detect and study the magnetic fields in galaxy clusters and in the less dense environments surrounding them using Faraday Rotation Measures. To this end, we produce 3-dimensional magnetic field models for galaxy clusters of different masses and in different stages of their evolution, and derive mock rotation measure observations of background radiogalaxies. According to our results, already in phase I, we will be able to infer the magnetic field properties in galaxy clusters as a function of the cluster mass, down to $10^{13}$ solar-masses. Moreover, using cosmological simulations to model the gas density, we have computed the expected rotation measure through shock-fronts that occur in the intra-cluster medium during cluster mergers. The enhancement in the rotation measure due to the density jump will permit to constraint the magnetic field strength and structure after the shock passage. SKA observations of polarised sources located behind galaxy clusters will answer several questions about the magnetic field strength and structure in galaxy clusters, and its evolution with cosmic time.

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

Rotation measure (RM) grids of extragalactic radio sources have been widely used for studying cosmic magnetism. But their potential for exploring the intergalactic magnetic field (IGMF) in filaments of galaxies is unclear, since other Faraday-rotation media such as the radio source itself, intervening galaxies, and the interstellar medium of our Galaxy are all significant contributors. We study statistical techniques for discriminating the Faraday rotation of filaments from other sources of Faraday rotation in future large-scale surveys of radio polarization. We consider a 30 degree times 30 degree field-of-view toward the south Galactic pole, while varying the number of sources detected in both present and future observations. We select sources located at high redshifts and toward which depolarization and optical absorption systems are not observed, so as to reduce the RM contributions from the sources and intervening galaxies. It is found that a high-pass filter can satisfactorily reduce the RM contribution from the Galaxy, since the angular scale of this component toward high Galactic latitudes would be much larger than that expected for the IGMF. Present observations do not yet provide a sufficient source density to be able to estimate the RM of filaments. However, from the proposed approach with forthcoming surveys, we predict significant residuals of RM that should be ascribable to filaments. The predicted structure of the IGMF down to scales of 0.1 degree should be observable with data from the SKA, if we achieve selections of sources toward which sight lines do not contain intervening galaxies and RM errors are less than a few rad m-2.

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

Unveiling the intergalactic magnetic field (IGMF) in filaments of galaxies is a very important and challenging subject in modern astronomy. In order to probe the IGMF from rotation measures (RMs) of extragalactic radio sources, we need to separate RMs due to other origins such as the source, intervening galaxies, and our Galaxy. In this paper, we discuss observational strategies for the separation by means of Faraday tomography (Faraday RM Synthesis). We consider an observation of a single radio source such as a radio galaxy or a quasar viewed through the Galaxy and the cosmic web. We then compare the observation with another observation of a neighbor source with a small angular separation. Our simulations with simple models of the sources suggest that it would be not easy to detect the RM due to the IGMF of order ~ 1 rad/m/m, an expected value for the IGMF through a single filament. Contrary to it, we find that the RM of at least ~10 rad/m/m could be detected with the SKA or its pathfinders/precursors, if we achieve selections of ideal sources. These results would be improved if we incorporate decomposition techniques such as RMCLEAN and QU-fitting. We discuss feasibility of the strategies for cases with complex Galactic emissions as well as with effects of observational noise and radio frequency interferences.

We report the detection of an X-ray absorption feature near the galaxy M86 in the Virgo cluster. The absorber has a column density of 2-3 E20 /cm^2, and its position coincides with the peak of an intracluster HI cloud which was removed from the galaxy NGC 4388 presumably by ram pressure. These results indicate that the HI cloud is located in front of M86 along the line-of-sight, and suggest that the stripping was primarily created by an interaction between NGC 4388 and the hot plasmas of the Virgo cluster, not the M86 halo. By calculating an X-ray temperature map, we further detected an X-ray counterpart of the HI cloud up to about 3' south of M86. It has a temperature of 0.89 keV and a mass of ~4.5 E8 solar mass, exceeding the estimated HI gas mass. The high hot-to-cold gas ratio in the cloud indicates a significant evaporation of the HI gas, probably by thermal conduction from the hotter cluster plasma with a sub-Spitzer rate.

With multiband photometric data in public archives, we detected four intracluster star-forming regions in the Virgo cluster. Two of them were at a projected distance of 35 kpc away from NGC4388, and the other two were 66 kpc away. Our new spectroscopic observation revealed that their recession velocities were comparable to the ram-pressure-stripped tail of NGC4388 and confirmed their association. The stellar mass of the star-forming regions ranged from 10^4 - 10^4.5 M_sun except for that of the faintest one which would be <10^3 M_sun. The metallicity was comparable to the solar abundance, and the age of the stars was ~ 10^6.8 years. Their young stellar age meant that the star formation should have started after the gas was stripped from NGC4388. This implied in situ condensation of the stripped gas. We also found that two star-forming regions lay near the leading edge of a filamentary dark cloud. The extinction of the filament was smaller than that derived from the Balmer decrement of the star-forming regions. It implied that the dust in the filament would be locally dense around the star-forming regions.

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.

We study the Faraday rotation measure (RM) due to the Galactic magnetic field (GMF) toward high Galactic latitudes. The RM arises from the global, regular component as well as from the turbulent, random component of the GMF. We model the former based on observations and the latter using the data of magnetohydrodynamic turbulence simulations. For a large number of different GMF models, we produce mock RM maps around the Galactic poles and calculate various statistical quantities with the RM maps. We find that the observed medians of RMs toward the north and south Galactic poles, $\sim 0.0\pm 0.5 {\rm rad m^{-2}}$ and $\sim +6.3\pm 0.5 {\rm rad m^{-2}}$, are difficult to explain with any of our many alternate GMF models. The standard deviation of observed RMs, $\sim 9 {\rm rad m^{-2}}$, is clearly larger than that of simulated RMs. The second-order structure function of observed RMs is substantially larger than that of simulated RMs, especially at small angular scales. We discuss other possible contributions to RM toward high Galactic latitudes. Besides observational errors and the intrinsic RM of background radio sources against which RM is observed, we suggest that the RM due to the intergalactic magnetic field may account for a substantial fraction of the observed RM. Finally we note that reproducing the observed medians may require additional components or/and structures of the GMF that are not present in our models.

We carried out global three-dimensional magneto-hydrodynamic simulations of dynamo activities in galactic gaseous disks without assuming equatorial symmetry. Numerical results indicate the growth of azimuthal magnetic fields non-symmetric to the equatorial plane. As magneto-rotational instability (MRI) grows, the mean strength of magnetic fields is amplified until the magnetic pressure becomes as large as 10% of the gas pressure. When the local plasma $\beta$ ($ = p_{\rm gas}/p_{\rm mag}$) becomes less than 5 near the disk surface, magnetic flux escapes from the disk by Parker instability within one rotation period of the disk. The buoyant escape of coherent magnetic fields drives dynamo activities by generating disk magnetic fields with opposite polarity to satisfy the magnetic flux conservation. The flotation of the azimuthal magnetic flux from the disk and the subsequent amplification of disk magnetic field by MRI drive quasi-periodic reversal of azimuthal magnetic fields in timescale of 10 rotation period. Since the rotation speed decreases with radius, the interval between the reversal of azimuthal magnetic fields increases with radius. The rotation measure computed from the numerical results shows symmetry corresponding to a dipole field.

We present imaging simulations of the Sunyaev-Zel'dovich effect of galaxy clusters for the Atacama Large Millimeter/submillimeter Array (ALMA) including the Atacama Compact Array (ACA). In its most compact configuration at 90GHz, ALMA will resolve the intracluster medium with an effective angular resolution of 5 arcsec. It will provide a unique probe of shock fronts and relativistic electrons produced during cluster mergers at high redshifts, that are hard to spatially resolve by current and near-future X-ray detectors. Quality of image reconstruction is poor with the 12m array alone but improved significantly by adding ACA; expected sensitivity of the 12m array based on the thermal noise is not valid for the Sunyaev-Zel'dovich effect mapping unless accompanied by an ACA observation of at least equal duration. The observations above 100 GHz will become excessively time-consuming owing to the narrower beam size and the higher system temperature. On the other hand, significant improvement of the observing efficiency is expected once Band 1 is implemented in the future.

We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium in the merging galaxy cluster, 1E0657-56 (the Bullet cluster), using a series of N-body and hydrodynamic simulations. We find that the electron temperature at the shock layer associated with the X-ray sub peak (bullet) is quite different depending on the thermal relaxation model between electrons and ions; ~25 keV for the Coulomb thermal relaxation model and ~45 keV for the instantaneous thermal relaxation model in the simulations which reproduce the observed X-ray morphology. Furthermore, both of Fe xxv and Fe xxvi are overabundant compared with the ionization equilibrium state around the shock layer, and thus, the intensity ratio between Fe xxv and Fe xxvi K alpha lines are significantly altered from that in the ionization equilibrium state. We also carry out the simulations with various sets of merger parameters, and discuss a possible range of the non-equilibrium effects in this system. Our results could be tested with future X-ray observations such as Astro-H with better sensitivity in high energy band.

We investigate the Faraday rotation measure (RM) due to the intergalactic magnetic field (IGMF) through the cosmic web up to cosmological distances, using a model IGMF based on turbulence dynamo in the large-scale structure of the universe. By stacking the IGMF and gas density data up to redshift $z=5$ and taking account of the redshift distribution of polarized background radio sources against which the RM is measured, we simulate the sky map of the RM. The contribution from galaxy clusters is subtracted from the map, based on several different criteria of X-ray brightness and temperature. Our findings are as follows. The distribution of RM for radio sources of different redshifts shows that the root-mean-square (rms) value increases with redshift and saturates for $z \ga 1$. The saturated value is RM$_{\rm rms} \approx$ several ${\rm rad m^{-2}}$. The probability distribution function of $|{\rm RM}|$ follows the lognormal distribution. The power spectrum has a broad plateau over the angular scale of $\sim 1 - 0.1^\circ$ with a peak around $\sim 0.15^\circ$. The second-order structure function has a flat profile in the angular separation of $\ga 0.2^\circ$. Our results could provide useful insights for surveys to explore the IGMF with the Square Kilometer Array (SKA) and upcoming SKA pathfinders.

High-frequency, high-resolution imaging of the Sunyaev-Zel'dovich (SZ) effect is an important technique to study the complex structures of the atmospheres of merging galaxy clusters. Such observations are sensitive to the details of the electron spectrum. We show that the morphology of the SZ intensity maps in simulated galaxy clusters observed at 345 GHz, 600 GHz, and 857 GHz are significantly different because of SZ relativistic corrections. These differences can be revealed by high-resolution imaging instruments. We calculate relativistically corrected SZ intensity maps of a simulated, massive, merging galaxy cluster and of the massive, merging clusters 1E0657-558 (the Bullet Cluster) and Abell 2219. The morphologies of the SZ intensity maps are remarkably different between 345 GHz and 857 GHz for each merging cluster. We show that high-resolution imaging observations of the SZ intensity maps at these frequencies, obtainable with the LABOCA and HERSCHEL-SPIRE instruments, allow to fully exploit the astrophysical relevance of the predicted SZ morphological effect.

Measurements of the Sunyaev-Zel'dovich (hereafter SZ) effect distortion of the cosmic microwave background provide methods to derive the gas pressure and temperature of galaxy clusters. Here we study the ability of SZ effect observations to derive the electron distribution function (DF) in massive galaxy clusters. Our calculations of the SZ effect include relativistic corrections considered within the framework of the Wright formalism and use a decomposition technique of electron DFs into Fourier series. Using multi-frequency measurements of the SZ effect, we find the solution of a linear system of equations that is used to derive the Fourier coefficients; we further analyze different frequency samples to decrease uncertainties in Fourier coefficient estimations. We propose a method to derive DFs of electrons using SZ multi-frequency observations of massive galaxy clusters. We found that the best frequency sample to derive an electron DF includes high frequencies $\nu$=375, 600, 700, 857 GHz. We show that it is possible to distinguish a Juttner DF from a Maxwell-Bolzman DF as well as from a Juttner DF with the second electron population by means of SZ observations for the best frequency sample if the precision of SZ intensity measurements is less than 0.1%. We demonstrate by means of 3D hydrodynamic numerical simulations of a hot merging galaxy cluster that the morphologies of SZ intensity maps are different for frequencies $\nu$=375, 600, 700, 857 GHz. We stress that measurements of SZ intensities at these frequencies are a promising tool for studying electron distribution functions in galaxy clusters.

The Sunyaev-Zel'dovich (hereafter SZ) effect is a promising tool to derive the gas temperature of galaxy clusters. The approximation of a spherically symmetric gas distribution is usually used to determine the temperature structure of galaxy clusters, but this approximation cannot properly describe merging galaxy clusters. The methods used so far, which do not assume the spherically symmetric distribution, permit us to derive 2D temperature maps of merging galaxy clusters. In this paper, we propose a method to derive the standard temperature deviation and temperature variance along the line-of-sight, which permits us to analyze the 3D temperature structure of galaxy clusters by means of the thermal SZ effect. We also propose a method to reveal merger shock waves in galaxy clusters by analyzing the presence of temperature inhomogeneities along the line-of-sight.

Studying the nature and origin of the intergalactic magnetic field (IGMF) is an outstanding problem of cosmology. Measuring Faraday rotation would be a promising method to explore the IGMF in the large-scale structure (LSS) of the universe. We investigated the Faraday rotation measure (RM) due to the IGMF in filaments of galaxies using simulations for cosmological structure formation. We employed a model IGMF based on turbulence dynamo in the LSS of the universe; it has an average strength of $< B > \sim 10$ nG and a coherence length of several $\times\ 100\ h^{-1}$ kpc in filaments. With the coherence length smaller than path length, the inducement of RM would be a random walk process, and we found that the resultant RM is dominantly contributed by the density peak along line of sight. The rms of RM through filaments at the present universe was predicted to be $\sim 1\ {\rm rad\ m^{-2}}$. In addition, we predicted that the probability distribution function of $|{\rm RM}|$ through filaments follows the log-normal distribution, and the power spectrum of RM in the local universe peaks at a scale of $\sim 1\ h^{-1}$ Mpc. Our prediction of RM could be tested with future instruments.

We investigate a non-equilibrium ionization state and an electron-ion two-temperature structure of the intracluster medium (ICM) in merging galaxy clusters using a series of N-body and hydrodynamic simulations. Mergers with various sets of mass ratios and impact parameters are systematically investigated, and it is found that, in most cases, ICM significantly departs from the ionization equilibrium state at the shock layers with a Mach number of ~1.5-2.0 in the outskirts of the clusters, and the shock layers with a Mach number of ~2-4 in front of the ICM cores. Accordingly, the intensity ratio between Fe xxv and Fe xxvi K alpha line emissions is significantly altered from that in the ionization equilibrium state. If the effect of the two-temperature structure of ICM is incorporated, the electron temperature is ~10-20 % and ~30-50 % lower than the mean temperature of ICM at the shock layers in the outskirts and in front of the ICM cores, respectively, and the deviation from the ionization equilibrium state becomes larger. We also address the dependence of the intensity ratio on the viewing angle with respect to the merging plane.

We investigate a non-equilibrium ionization state and two-temperature structure of the intracluster medium in the linked region of Abell 399/401, using a series of N-body + SPH simulations, and find that there exist significant shock layers at the edge of the linked region, and that the ionization state of iron departs from the ionization equilibrium state at the shock layers and around the center of the linked region. As for the two-temperature structure, an obvious difference of temperature between electrons and ions is found in the edge of the linked regions. K alpha line emissions of Fe xxiv and Fe xxv are not severely affected by the deviation from the ionization equilibrium state around the center of the linked region, suggesting that the detection of relatively high metallicity in this area cannot be ascribed to the non-equilibrium ionization state of the intracluster medium. On the other hand, the K alpha emissions are significantly deviated from the equilibrium values at the shock layers, and the intensity ratio of K alpha lines between Fe xxiv-xxv and Fe xxvi is found to be significantly altered from that in the ionization equilibrium state.

We investigate the core structure of radiatively cooling intracluster gas, using a hydrodynamics code. We calculate evolution of model clusters of the initial core radii 160--300 kpc until the initial central cooling time, and analyze the resultant clusters using the double beta-model as done by observational studies. It is found that the core-size distribution thus obtained shows two peaks \sim 30--100 kpc and \sim 100--300 kpc and marginally can reproduce the observed distribution which exhibits two distinct peaks around \sim 50 kpc and \sim 200 kpc. This result may suggest radiative-cooling origin for small cores, while cooling is yet insignificant in the clusters of large cores. It should be noted that the small core peak is reproduced by clusters that are still keeping quasi-hydrostatic balance before the initial central cooling time has elapsed.

We investigate core structures of X-ray emitting intracluster gas based on the so-called beta-model, which is an isothermal hydrostatic model often used in observational studies. We reconsider the beta-model and find that the virial temperature T_vir of a cluster may be represented better by beta T_X than T_X, where beta is the parameter obtained from the X-ray surface brightness and T_X is the emission-weighted mean temperature of the gas. We investigate 121 clusters observed by ROSAT and ASCA and find that the luminosity-temperature relation L_X - beta T_X is less steep than L_X - T_X. We classify the clusters into two core-size groups in order to investigate their properties in detail. While in the larger core group the core radius is marginally proportional to the virial radius, no significant relation is found for the smaller core group. This may suggest that the smaller cores reflect the presence of cD galaxies, effect of radiative cooling or asymmetry in the surface brightness. We examine such possibilities, and find that the clusters of smaller cores have shorter cooling time than the Hubble time, while no significant correlation is found with cD or asymmetry. We carry out hydrodynamical calculations to simulate the beta-model, intending to see the behavior of the isothermal gas under the gravitational potential including the dark matter and galaxies with or without a central cD galaxy. Calculations show r_c \propto r_vir and T_vir \simeq beta T consistently with our consideration to the beta-model. Also is found from calculations that the presence of a large cD galaxy may form a gas core \sim 40 kpc, which seems too small to account for the range of the core sizes, 40--80 kpc, of the smaller core group.