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.

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

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 present a joint X-ray, optical and weak-lensing analysis for X-ray luminous galaxy clusters selected from the MCXC (Meta-Catalog of X-Ray Detected Clusters of Galaxies) cluster catalog in the Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) survey field with S16A data, As a pilot study of our planned series papers, we measure hydrostatic equilibrium (H.E.) masses using XMM-Newton data for four clusters in the current coverage area out of a sample of 22 MCXC clusters. We additionally analyze a non-MCXC cluster associated with one MCXC cluster. We show that H.E. masses for the MCXC clusters are correlated with cluster richness from the CAMIRA catalog (Oguri et al. 2017), while that for the non-MCXC cluster deviates from the scaling relation. The mass normalization of the relationship between the cluster richness and H.E. mass is compatible with one inferred by matching CAMIRA cluster abundance with a theoretical halo mass function. The mean gas mass fraction based on H.E. masses for the MCXC clusters is $\langle f_{\rm gas} \rangle = 0.125\pm0.012$ at spherical overdensity $\Delta=500$, which is $\sim80-90$ percent of the cosmic mean baryon fraction, $\Omega_b/\Omega_m$, measured by cosmic microwave background experiments. We find that the mean baryon fraction estimated from X-ray and HSC-SSP optical data is comparable to $\Omega_b/\Omega_m$. A weak-lensing shear catalog of background galaxies, combined with photometric redshifts, is currently available only for three clusters in our sample. Hydrostatic equilibrium masses roughly agree with weak-lensing masses, albeit with large uncertainty. This study demonstrates that further multiwavelength study for a large sample of clusters using X-ray, HSC-SSP optical and weak lensing data will enable us to understand cluster physics and utilize cluster-based cosmology.

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.

The results of Suzaku observations of the outskirts of Abell 3395 including a large-scale structure filament toward Abell 3391 are presented. We measured temperature and abundance distributions from the southern outskirt of Abell 3395 to the north at the virial radius, where a filament structure has been found in the former X-ray and Sunyaev-Zel'dovich effect observations between Abell 3391 and 3395. The overall temperature structure is consistent with the universal profile proposed by Okabe et al.(2014) for relaxed clusters except for the filament region. A hint of the ICM heating is found between the two clusters, which might be due to the interaction of them in the early phase of a cluster merger. Although we obtained relatively low metal abundance of $Z=0.169^{+0.164+0.009+0.018 }_{-0.150-0.004-0.015 }$ solar, where the first, second, and third errors are statistical, cosmic X-ray background systematic, and non X-ray background systematic, respectively, at the virial radius in the filament, our results are still consistent with the former results of other clusters ($Z \sim 0.3$ solar) within errors. Therefore, our results are also consistent with the early enrichment scenario. We estimated Compton $y$ parameters only from X-ray results in the region between Abell 3391 and 3395 assuming a simple geometry. They are smaller than the previous SZ results with Planck satellite. The difference could be attributed to a more elaborate geometry such as a filament inclined to the line-of-sight direction, or underestimation of the X-ray temperature because of the unresolved multi-temperature structures or undetected hot X-ray emission of the shock heated gas.

We present the results of Suzaku and Chandra observations of the galaxy cluster RXC J1053.7+5453 ($z=0.0704$), which contains a radio relic. The radio relic is located at the distance of $\sim 540$ kpc from the X-ray peak toward the west. We measured the temperature of this cluster for the first time. The resultant temperature in the center is $ \sim 1.3$ keV, which is lower than the value expected from the X-ray luminosity - temperature and the velocity dispersion - temperature relation. Though we did not find a significant temperature jump at the outer edge of the relic, our results suggest that the temperature decreases outward across the relic. Assuming the existence of the shock at the relic, its Mach number becomes $M \simeq 1.4 $. A possible spatial variation of Mach number along the relic is suggested. Additionally, a sharp surface brightness edge is found at the distance of $\sim 160$ kpc from the X-ray peak toward the west in the Chandra image. We performed X-ray spectral and surface brightness analyses around the edge with Suzaku and Chandra data, respectively. The obtained surface brightness and temperature profiles suggest that this edge is not a shock but likely a cold front. Alternatively, it cannot be ruled out that thermal pressure is really discontinuous across the edge. In this case, if the pressure across the surface brightness edge is in equilibrium, other forms of pressure sources, such as cosmic-rays, are necessary. We searched for the non-thermal inverse Compton component in the relic region. Assuming the photon index $ \Gamma = 2.0$, the resultant upper limit of the flux is $1.9 \times 10^{-14} {\rm erg \ s^{-1} \ cm^{-2}}$ for $4.50 \times 10^{-3} {\rm \ deg^{2}}$ area in the 0.3-10 keV band, which implies that the lower limit of magnetic field strength becomes $ 0.7 {\rm \ \mu G}$.

The XXL survey currently covers two 25 sq. deg. patches with XMM observations of ~10ks. We summarise the scientific results associated with the first release of the XXL data set, that occurred mid 2016. We review several arguments for increasing the survey depth to 40 ks during the next decade of XMM operations. X-ray (z<2) cluster, (z<4) AGN and cosmic background survey science will then benefit from an extraordinary data reservoir. This, combined with deep multi-$\lambda$ observations, will lead to solid standalone cosmological constraints and provide a wealth of information on the formation and evolution of AGN, clusters and the X-ray background. In particular, it will offer a unique opportunity to pinpoint the z>1 cluster density. It will eventually constitute a reference study and an ideal calibration field for the upcoming eROSITA and Euclid missions.

We present the results of deep 140 ks Suzaku X-ray observations of the north-east (NE) radio relic of the merging galaxy cluster Abell2255. The temperature structure of Abell2255 is measured out to 0.9 times the virial radius (1.9 Mpc) in the NE direction for the first time. The Suzaku temperature map of the central region suggests a complex temperature distribution, which agrees with previous work. Additionally, on a larger-scale, we confirm that the temperature drops from 6 keV around the cluster center to 3 keV at the outskirts, with two discontinuities at {\it r}$\sim$5\arcmin~(450 kpc) and $\sim$12\arcmin~(1100 kpc) from the cluster center. Their locations coincide with surface brightness discontinuities marginally detected in the XMM-Newton image, which indicates the presence of shock structures. From the temperature drop, we estimate the Mach numbers to be ${\cal M}_{\rm inner}\sim$1.2 and, ${\cal M}_{\rm outer}\sim$1.4. The first structure is most likely related to the large cluster core region ($\sim$350--430 kpc), and its Mach number is consistent with the XMM-Newton observation (${\cal M}\sim$1.24: Sakelliou & Ponman 2006). Our detection of the second temperature jump, based on the Suzaku key project observation, shows the presence of a shock structure across the NE radio relic. This indicates a connection between the shock structure and the relativistic electrons that generate radio emission. Across the NE radio relic, however, we find a significantly lower temperature ratio ($T_1/T_2\sim1.44\pm0.16$ corresponds to~${\cal M}_{\rm X-ray}\sim1.4$) than the value expected from radio wavelengths, based on the standard diffusive shock acceleration mechanism ($T_1/T_2>$ 3.2 or ${\cal M}_{\rm Radio}>$ 2.8).

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.

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.

Abell 548W, one of the galaxy clusters located in the Abell 548 region, has about an order of magnitude lower X-ray luminosity compared to ordinal clusters in view of the well known intracluster medium (ICM) temperature vs X-ray luminosity (kT-L_X) relation. The cluster hosts a pair of diffuse radio sources to the north west and north, both about 10' apart from the cluster center. They are candidate radio relics, frequently associated with merging clusters. A Suzaku deep observation with exposure of 84.4 ks was performed to search signatures for merging in this cluster. The XIS detectors successfully detected the ICM emission out to 16' from the cluster center. The temperature is ~3.6 keV around its center, and ~2 keV at the outermost regions. The hot region (~6 keV) aside the relic candidates shifted to the cluster center reported by XMM-Newton was not seen in the Suzaku data, although its temperature of 3.6 keV itself is higher than the average temperature of 2.5 keV around the radio sources. In addition, a signature of a cool (kT ~0.9 keV) component was found around the north west source. A marginal temperature jump at its outer-edge was also found, consistent with the canonical idea of shock acceleration origin of the radio relics. The cluster has among the highest central entropy of ~400 keV cm^2 and is one of the so-called low surface brightness clusters. Taking into account the fact that its shape itself is relatively circular and smooth and also its temperature structure is nearly flat, possible scenarios for merging is discussed.

Radio relics are synchrotron emission found on the periphery of galaxy clusters. From the position and the morphology, it is often believed that the relics are generated by cosmic-ray (CR) electrons accelerated at shocks through a diffusive shock acceleration (DSA) mechanism. However, some radio relics have harder spectra than the prediction of the standard DSA model. One example is observed in the cluster 1RXS J0603.3+4214, which is often called the ``Toothbrush Cluster''. Interestingly, the position of the relic is shifted from that of a possible shock. In this study, we show that these discrepancies in the spectrum and the position can be solved if turbulent (re)acceleration is very effective behind the shock. This means that for some relics turbulent reacceleration may be the main mechanism to produce high-energy electrons, contrary to the common belief that it is the DSA. Moreover, we show that for efficient reacceleration, the effective mean free path of the electrons has to be much smaller than their Coulomb mean free path. We also study the merging cluster 1E 0657-56 or the ``Bullet Cluster'', in which a radio relic has not been found at the position of the prominent shock ahead of the bullet. We indicate that a possible relic at the shock is obscured by the observed large radio halo that is generated by strong turbulence behind the shock. We propose a simple explanation if the morphological differences of radio emission among the Toothbrush, the Bullet, and the Sausage (CIZA J2242.8+5301) Clusters.

We perform a weak-lensing study of the nearby cool-core galaxy clusters, Hydra A ($z=0.0538$) and A478 ($z=0.0881$), of which brightest cluster galaxies (BCGs) host powerful activities of active galactic nuclei (AGNs). For each cluster, the observed tangential shear profile is well described either by a single Navarro--Frenk--White model or a two-component model including the BCG as an unresolved point mass. For A478, we determine the BCG and its host-halo masses from a joint fit to weak-lensing and stellar photometry measurements. We find that the choice of initial mass functions (IMFs) can introduce a factor of two uncertainty in the BCG mass, whereas the BCG host halo mass is well constrained by data. We perform a joint analysis of weak-lensing and stellar kinematics data available for the Hydra A cluster, which allows us to constrain the central mass profile without assuming specific IMFs. We find that the central mass profile ($r<300$kpc) determined from the joint analysis is in excellent agreement with those from independent measurements, including dynamical masses estimated from the cold gas disk component, X-ray hydrostatic total mass estimates, and the central stellar mass estimated based on the Salpeter IMF. The observed dark-matter fraction around the BCG for Hydra A is found to be smaller than those predicted by adiabatic contraction models, suggesting the importance of other physical processes, such as the the AGN feedback and/or dissipationless mergers.

We present the results of Suzaku observations of the galaxy cluster 1RXS J0603.3+4214 with "toothbrush" radio relic. Although a shock with Mach number $M \simeq 4$ is expected at the outer edge of the relic from the radio observation, our temperature measurements of the intracluster medium indicate a weaker temperature difference than what is expected. The Mach number estimated from the temperature difference at the outer edge of the relic is $M \simeq 1.5$, which is significantly lower than the value estimated from the radio data even considering both statistical and systematic errors. This suggests that a diffusive shock acceleration theory in the linear test particle regime, which is commonly used to link the radio spectral index to the Mach number, is invalid for this relic. We also measured the temperature difference across the western part of the relic, where a shock with $M \simeq 1.6$ is suggested from the X-ray surface brightness analysis of the XMM-Newton data, and obtained consistent results in an independent way. We searched for the non-thermal inverse Compton component in the relic region and the resultant upper limit on the flux is $2.4 \times 10^{-13}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3-10 keV band. The lower limit of the magnetic field strength becomes 1.6 $\mu$G, which means that magnetic energy density could be more than a few $\% $ of the thermal energy.

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.

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.

We conduct a joint X-ray and weak-lensing study of four relaxed galaxy clusters (Hydra A, A478, A1689 and A1835) observed by both Suzaku and Subaru out to virial radii, with an aim to understand recently-discovered unexpected feature of the ICM in cluster outskirts. We show that the average hydrostatic-to-lensing total mass ratio for the four clusters decreases from \sim 70% to \sim 40% as the overdensity contrast decreases from 500 to the virial value.The average gas mass fraction from lensing total mass estimates increases with cluster radius and agrees with the cosmic mean baryon fraction within the virial radius, whereas the X-ray-based gas fraction considerably exceeds the cosmic values due to underestimation of the hydrostatic mass. We also develop a new advanced method for determining normalized cluster radial profiles for multiple X-ray observables by simultaneously taking into account both their radial dependence and multivariate scaling relations with weak-lensing masses. Although the four clusters span a range of halo mass, concentration, X-ray luminosity and redshift, we find that the gas entropy, pressure, temperature and density profiles are all remarkably self-similar when scaled with the lensing M_200 mass and r_200 radius.The entropy monotonically increases out to \sim 0.5r_200 following the accretion shock heating model K(r)\propto r^1.1, and flattens at \simgt 0.5r_200.The universality of the scaled entropy profiles indicates that the thermalization mechanism over the entire cluster region (>0.1r_200) is controlled by gravitation in a common to all clusters, although the heating efficiency in the outskirts needs to be modified from the standard law.The bivariate scaling functions of the gas density and temperature reveal that the flattening of the outskirts entropy profile is caused by the steepening of the temperature, rather than the flattening of the gas density.

The central cD galaxy of the Hydra A cluster has one of the most powerful active galactic nuclei (AGNs) in the nearby Universe (z <0.2). We report on the discovery of a dust lane in the cD galaxy using Subaru telescope. The i'-band image shows the existence of a dark band of the size of 3.6"x 0.7" (4 kpc x 0.8 kpc), which appears to be quite similar to the dust lane observed in Centaurus A. The morphology indicates that the cold disk that seen as the dust lane is almost edge-on and rotates around the AGN. Since the minor axis of the dust lane is nearly parallel to the radio jets emerging from the AGN, the disk is probably feeding its gas into the central black hole. From the absorption, we estimate the hydrogen column density of the lane is N_H=2.0 x 10^21 cm^-2, and the mass of the disk is ~8 x 10^7 M_sun. The column density is consistent with constraints obtained from Chandra X-ray observations. The age of the disk is >~ 4 x 10^7 yr. The position angle of the disk and the galaxy's photometric axis are misaligned, which may imply that the cold gas in the disk is brought via galaxy mergers. Our observations may indicate that the supply of cold gas by galaxy mergers is required for the most intensive feedback from AGNs.

We present the results from Suzaku observations of the Perseus galaxy cluster, which is relatively close, the brightest in the X-ray sky and a relaxed object with a cool core. A number of exposures of central regions and offset pointing with the X-ray Imaging Spectrometer cover a region within radii of 20'-30'. The central data are used to evaluate the instrumental energy-scale calibration with accuracy confirmed to within around 300 km/s, by the spatial and temporal variation of the instruments. These deep and well-calibrated data are used to measure X-ray redshifts of the intracluster medium. A hint of gas bulk motion, with radial velocity of about -(150-300) km/s, relative to the main system was found at 2-4 arcmin (45-90kpc) west of the cluster center, where an X-ray excess and a cold front were found previously. No other velocity structure was discovered. Over spatial scales of 50-100kpc and within 200kpc radii of the center, the gas-radial-velocity variation is below 300 km/s, while over scales of 400 kpc within 600 kpc radii, the variation is below 600 km/s. These X-ray redshift distributions are compared spatially with those of optical member galaxies for the first time in galaxy clusters. Based on X-ray line widths gas turbulent velocities within these regions are also constrained within 1000-3000 km/s. These results of gas dynamics in the core and larger scales in association with cluster merger activities are discussed and future potential of high-energy resolution spectroscopy with ASTRO-H is considered.

We present results of four-pointing Suzaku X-ray observations (total ~200 ks) of the intracluster medium (ICM) in the Abell 1835 galaxy cluster (kT ~ 8 keV, z = 0.253) out to the virial radius (r_vir ~ 2.9 Mpc) and beyond. Faint X-ray emission from the ICM out to r_vir is detected. The temperature gradually decreases with radius from ~8 keV in the inner region to ~2 keV at r_vir. The entropy profile is shown to flatten beyond r_500, in disagreement with the r_1.1 dependence predicted from the accretion shock heating model. The thermal pressure profile in the range 0.3r_500 < r < r_vir agrees well with that obtained from the stacked Sunyaev-Zel'dovich effect observations with the Planck satellite. The hydrostatic mass profile in the cluster outskirts (r_500 < r < r_vir) falls well short of the weak lensing one derived from Subaru/Suprime-Cam observations, showing an unphysical decrease with radius. The gas mass fraction at r_vir defined with the lensing total mass agrees with the cosmic baryon fraction from the WMAP 7-year data. All these results indicate, rather than the gas-clumping effect, that the bulk of the ICM in the cluster outskirts is far from hydrostatic equilibrium and infalling matter retained some of its kinetic energy. Finally, combining with our recent Suzaku and lensing analysis of Abell 1689, a cluster of similar mass, temperature, and redshift, we show that the cluster temperature distribution in the outskirts is significantly correlated with the galaxy density field in the surrounding large-scale environment at (1-2)r_vir.

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 report Suzaku observations of the northern half of the Hydra A cluster out to ~1.4 Mpc, reaching the virial radius. This is the first Suzaku observations of a medium-size (kT ~3 keV) cluster out to the virial radius. Two observations were conducted, north-west and north-east offsets, which continue in a filament direction and a void direction of the large-scale structure of the Universe, respectively. The X-ray emission and distribution of galaxies elongate in the filament direction. The temperature profiles in the two directions are mostly consistent with each other within the error bars and drop to 1.5 keV at 1.5 r_500. As observed by Suzaku in hot clusters, the entropy profile becomes flatter beyond r_500, in disagreement with the r^1.1 relationship that is expected from accretion shock heating models. When scaled with the average intracluster medium (ICM) temperature, the entropy profiles of clusters observed with Suzaku are universal and do not depend on system mass. The hydrostatic mass values in the void and filament directions are in good agreement, and the Navarro, Frenk, and White universal mass profile represents the hydrostatic mass distribution up to ~ 2 r_500. Beyond r_500, the ratio of gas mass to hydrostatic mass exceeds the result of the Wilkinson microwave anisotropy probe, and at r_100, these ratios in the filament and void directions reach 0.4 and 0.3, respectively. We discuss possible deviations from hydrostatic equilibrium at cluster outskirts. We derived radial profiles of the gasmass- to-light ratio and iron-mass-to-light ratio out to the virial radius. Within r_500, the iron-mass-to-light ratio of the Hydra A cluster was compared with those in other clusters observed with Suzaku.

We report on a Suzaku measurement of the shock feature associated with the western radio relic in the merging cluster A3376. The temperature profile is characterized by an almost flat radial shape with kT ~ 4 keV within 0.5 r200 and a rise by about 1 keV inside the radio relic. Across the relic region (0.6-0.8 r200), the temperature shows a remarkable drop from about 4.7 keV to 1.3 keV. This is a clear evidence that the radio relic really corresponds to a shock front possibly caused by a past major merger. The observed sharp changes of the temperature and electron density indicate the Mach number M~3. The radial entropy profile is flatter than the prediction (r^1.1) of numerical simulations within 0.5 r200}, and becomes steeper around the relic region. These observed features and time-scale estimation consistently imply that the ICM around the radio relic has experienced a merger shock and is in the middle of the process of dynamical and thermal relaxation.

We present the results of Suzaku X-ray follow-up observation of weak-lensing-detected halos in the field around galaxy cluster ZwCl0823.2+0425. We clearly detected X-ray emission associated with most of these halos and determined their detailed physical parameters such as X-ray luminosity, temperature, and metal abundance, for the first time. We find that the X-ray luminosity - temperature relation for these halos agrees with former typical results. With mass determined from the weak gravitational lensing data, the mass-temperature relation for them is also investigated and found to be consistent with the prediction from a simple self-similar model and results of the previous studies with both lensing and X-ray data. We would like to emphasize that the self-similar scaling relation of mass and temperature is shown here for the first time using a weak-lensing selected sample, whereas previous studies of the mass scaling relation used X-ray-selected samples of clusters. Therefore, our study demonstrates importance of X-ray follow-up observations of shear-selected clusters, and shows that a joint X-ray and lensing analysis will be crucial for clusters discovered by the forthcoming weak-lensing surveys, such as the one planned with Subaru/Hyper-Suprime-Cam.

We investigate the impact of mergers on the mass estimation of galaxy clusters using $N$-body + hydrodynamical simulation data. We estimate virial mass from these data and compare it with real mass. When the smaller subcluster's mass is larger than a quarter of that of the larger one, virial mass can be larger than twice of the real mass. The results strongly depend on the observational directions, because of anisotropic velocity distribution of the member galaxies. We also make the X-ray surface brightness and spectroscopic-like temperature maps from the simulation data. The mass profile is estimated from these data on the assumption of hydrostatic equilibrium. In general, mass estimation with X-ray data gives us better results than virial mass estimation. The dependence upon observational directions is weaker than in case of virial mass estimation. When the system is observed along the collision axis, the projected mass tends to be underestimated. This fact should be noted especially when the virial and/or X-ray mass are compared with gravitational lensing results.

(Abridged) We present results of Suzaku observations of the intracluster medium (ICM) in Abell 1689, combined with complementary analysis of the SDSS data and weak and strong lensing analysis of Subaru/Suprime-Cam and HST/ACS observations. Faint X-ray emission from the ICM around the virial radius is detected at 4.0 sigma significance. We find anisotropic gas temperature and entropy distributions in cluster outskirts correlated with large-scale structure of galaxies. The high temperature and entropy region in the northeastern (NE) outskirts is connected to an overdense filamentary structure. The outskirt regions in contact with low density void environments have low gas temperatures and entropies, deviating from hydrostatic equilibrium. These results suggest that thermalization of the ICM occurs faster along the filamentary structures than the void regions. A joint X-ray and lensing analysis shows that the hydrostatic mass is $\sim60-90%$ of spherical lensing one but comparable to a triaxial halo mass within errors in $0.6r_{2500} \simlt r \simlt 0.8r_{500}$, and that it is significantly biased as low as $\simlt60%$ within $0.4r_{2500}$, irrespective of mass models. The thermal gas pressure within $r_{500}$ is, at most, $\sim50$--60% of the total pressure to balance fully the gravity of the spherical lensing mass, and $\sim30$--40% around the virial radius. Although these constitute lower limits when one considers the possible halo triaxiality, these small relative contributions of thermal pressure would require additional sources of pressure, such as bulk and/or turbulent motions.

We present the results of Suzaku observation of the radio halo cluster Abell 2319. The metal abundance in the central cool region is found to be higher than the surrounding region, which was not resolved in the former studies. We confirm that the line-of-sight velocities of the intracluster medium in the observed region are consistent with those of the member galaxies of entire A2319 and A2319A subgroup for the first time, though any velocity difference within the region is not detected. On the other hand, we do not find any signs of gas motion relevant to A2319B subgroup. Hard X-ray emission from the cluster is clearly detected, but its spectrum is likely thermal. Assuming a simple single temperature model for the thermal component, we find that the upper limit of the non-thermal inverse Compton component becomes $2.6 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ in the 10-40 keV band, which means that the lower limit of the magnetic field is 0.19 $\mu$G with the radio spectral index 0.92. Although the results slightly depend on the detailed spectral modeling, it is robust that the upper limit of the power-law component flux and lower limit of the magnetic field strength become $\sim 3 \times 10^{-11}$ erg s$^{-1}$ cm$^{-2}$ and $\sim 0.2 \mu$G, respectively. Considering the lack of a significant amount of very hot ($\sim 20$ keV) gas and the strong bulk flow motion, it is more likely that the relativistic non-thermal electrons responsible for the radio halo are accelerated through the intracluster turbulence rather than the shocks.

Wide-band Suzaku data on the merging cluster Abell 3667 were examined for hard X-ray emission in excess to the known thermal component. Suzaku detected X-ray signals in the wide energy band from 0.5 to 40 keV. The hard X-ray (> 10 keV) flux observed by the HXD around the cluster center cannot be explained by a simple extension of the thermal emission with average temperature of ~7 keV. The emission is most likely an emission from a very hot (kT > 13.2 keV) thermal component around the cluster center, produced via a strong heating process in the merger. In the north-west radio relic, no signature of non-thermal emission was observed. Using the HXD, the overall upper-limit flux within a 34'x34' field-of-view around the relic is derived to be 5.3e-12 erg s-1 cm-2 in the 10-40 keV band, after subtracting the ICM contribution estimated using the XIS or the XMM-Newton spectra. Directly on the relic region, the upper limit is further tightened by the XIS data to be less than 7.3e-13 erg s-1 cm-2, when converted into the 10--40 keV band. The latter value suggest that the average magnetic field within the relic is higher than 1.6 uG. The non-thermal pressure due to magnetic fields and relativistic electrons may be as large as ~20% of the thermal pressure in the region.

We present results from N-body + magnetohydrodynamical simulations of merging clusters of galaxies. We find that cluster mergers cause various characteristic magnetic field structures because of the strong bulk flows in the intracluster medium. The moving substructures result in cool regions surrounded by the magnetic field. These will be recognized as magnetized cold fronts in the observational point of view. A relatively ordered magnetic field structure is generated just behind the moving substructure. Eddy-like field configurations are also formed by Kelvin-Helmholtz instabilities. These features are similarly seen even in off-center mergers though the detailed structures change slightly. The above-mentioned characteristic magnetic field structures are partly recognized in Faraday rotation measure maps. The higher absolute values of the rotation measure are expected when observed along the collision axis, because of the elongated density distribution and relatively ordered field structure along the axis. The rotation measure maps on the cosmic microwave background radiation, which covers clusters entirely, could be useful probes of not only the magnetic field structures but also the internal dynamics of the intracluster medium.

We present the analysis of a Suzaku observation of the Ophiuchus galaxy cluster. We confirmed that the cluster has a cool core. While the temperature of the intracluster medium (ICM) decreases toward the center, the metal abundance increases. Except for the core (r<~50 kpc), the cluster is hot (~9-10 keV) and is almost isothermal for r<~1 Mpc; the latter contradicts a previous study. We do not detect the variation of the redshift of the ICM in the cluster; the upper limit of the velocity difference is 3000 km s^-1. The iron line ratios in X-ray spectra indicate that the ICM has reached the ionization equilibrium state. From these results, we conclude that the Ophiuchus cluster is not a major merger cluster but one of the hottest clusters with a cool core. We obtain the upper limit of non-thermal emission from the cluster, which is consistent with both the recent claimed detection with INTEGRAL and the recent upper limits with the Swift/BAT. If the cluster has bright non-thermal emission as suggested by the INTEGRAL measurement, it is probably not due to a recent major cluster merger.

We present an analysis of a Suzaku observation of the link region between the galaxy clusters A399 and A401. We obtained the metallicity of the intracluster medium (ICM) up to the cluster virial radii for the first time. We determine the metallicity where the virial radii of the two clusters cross each other (~2 Mpc away from their centers) and found that it is comparable to that in their inner regions (~0.2 Zsun). It is unlikely that the uniformity of metallicity up to the virial radii is due to mixing caused by a cluster collision. Since the ram-pressure is too small to strip the interstellar medium of galaxies around the virial radius of a cluster, the fairly high metallicity that we found there indicates that the metals in the ICM are not transported from member galaxies by ram-pressure stripping. Instead, the uniformity suggests that the proto-cluster region was extensively polluted with metals by extremely powerful outflows (superwinds) from galaxies before the clusters formed. We also searched for the oxygen emission from the warm--hot intergalactic medium in that region and obtained a strict upper limit of the hydrogen density (nH<4.1x10^-5 cm^-3).

We investigate the X-ray and mass distribution in the merging galaxy cluster 1E 0657-56. We study head-on collisions of two virialized clusters with an NFW density profile in the $\Lambda$CDM universe using an N-body + hydrodynamical code. A clear off-set of an X-ray peak from a mass peak, which is like what is reported in Clow et al. (2004), is first reproduced in the N-body + hydrodynamical simulations. We estimate the ram pressure-stripping conditions of the substructure in mergers of two NFW dark halos using a simple analytical model. We find that the ram pressure dominates the gravity of the substructure when the smaller cluster's mass is less than approximately one tenth of the larger cluster's mass. The characteristic X-ray and mass structures found in 1E 0657-56 suggest that neither the ram pressure nor the gravitational bound force overwhelms the other and that the mass ratio between the progenitors is near the critical value mentioned above.

We perform three dimensional hydrodynamical simulations of a moving substructure in a cluster of galaxies. We investigate dynamical evolution of the intracluster medium (ICM) in and around the substructure moving radially in the larger cluster's gravitational potential, and its observational consequences. After the substructure passes the larger cluster's center, a bow shock and clear contact discontinuity form in front of it. The contact discontinuity looks like a sharp cold front in the X-ray image synthesized from the simulation results. This agrees with a structure found in 1E 0657-56. The flow structure remains laminar before the first turnaround because the ram-pressure stripping is dominant over the development of Kelvin-Helmholtz instabilities on the boundary between the substructure and the ambient ICM. When a subcluster oscillates radially around the larger cluster's center, both Kelvin-Helmholtz and Rayleigh-Taylor instabilities develop well and the flow structure becomes highly turbulent. Around the turnaround, the subcluster's cold gas is pushed out of its potential well. Therefore, the cold gas component appears to be in front of the subcluster. A relatively blunt cold front appears in the simulated X-ray image because the contact discontinuity is perturbed by Rayleigh-Taylor instabilities. This can explain the ICM structure found in A168.

Clusters of galaxies form through major merger and/or absorption of smaller groups. In fact, some characteristic structures such as cold fronts, which are likely relevant to moving substructures, are found by {\it Chandra}. It is expected that moving substructures generate turbulence in the intracluster medium (ICM). Such turbulence probably plays a crucial role in mixture and transport of gas energy and heavy elements, and particle acceleration. The {\it Astro-E2} satellite, which is planned to be launched in 2005, will detect broadened lines due to turbulent motion. In order to explore the above-mentioned issues, it is important to investigate the generation processes and structure of ICM turbulence. We investigate the ICM dynamical evolution in and around a moving substructure with three-dimensional hydrodynamical simulations. Eddy-like structures develop near the boundary between the substructure and the ambient ICM through Kelvin-Helmholtz instabilities. Because of these structures, characteristic patterns appear in the line-of-sight velocity distribution of the ICM.

We present the results of a {\it Chandra} observation of the central region of Abell 3112. This cluster has a powerful radio source in the center and was believed to have a strong cooling flow. The X-ray image shows that the intracluster medium (ICM) is distributed smoothly on large scales, but has significant deviations from a simple concentric elliptical isophotal model near the center. Regions of excess emission appear to surround two lobe-like radio-emitting regions. This structure probably indicates that hot X-ray gas and radio lobes are interacting. From an analysis of the X-ray spectra in annuli, we found clear evidence for a temperature decrease and abundance increase toward the center. The X-ray spectrum of the central region is consistent with a single-temperature thermal plasma model. The contribution of X-ray emission from a multiphase cooling flow component with gas cooling to very low temperatures locally is limited to less than 10% of the total emission. However, the whole cluster spectrum indicates that the ICM is cooling significantly as a whole, but in only a limited temperature range ($\geq 2$ keV). Inside the cooling radius, the conduction timescales based on the Spitzer conductivity are shorter than the cooling timescales. We detect an X-ray point source in the cluster center which is coincident with the optical nucleus of the central cD galaxy and the core of the associated radio source. The X-ray spectrum of the central point source can be fit by a 1.3 keV thermal plasma and a power-law component whose photon index is 1.9. The thermal component is probably plasma associated with the cD galaxy. We attribute the power-law component to the central AGN.

We consider nonthermal emission from clusters of galaxies produced by particle acceleration by resonant scattering of Alfv\'{e}n waves driven by fluid turbulence through the Lighthill mechanism in the intracluster medium. We assume that the turbulence is driven by cluster mergers. We find that the resonant Alfv\'{e}n waves can accelerate electrons up to \gamma ~10^5 through resonant scattering. We also find that the turbulent resonant acceleration can give enough energy to electrons to produce the observed diffuse radio relic emission from clusters if the clusters have a pool of electrons with \gamma ~10^3. This mechanism can also explain the observed hard X-ray emission from clusters if the magnetic field in a cluster is small enough (~< \mu G) or the fluid turbulence spectrum is flatter than the Kolmogorov law. The fluid turbulence could be observed with Astro-E2 in the regions where diffuse radio emission is observed. Although non-gravitational heating before cluster formation (preheating) steepens a relation between radio luminosity and X-ray temperature, our predicted relation is still flatter than the observed one.

We developed a magneto-turbulent model for the cosmic ray (CR) electrons seen in the radio halo clusters of galaxies. Steady state momentum distribution functions of the CR electrons are calculated for given spectra of the turbulent Alfven waves. The radio spectrum produced by the obtained CR electron distribution is compared to the observed radio spectrum of the Coma radio halo. We find that the observed radio spectrum of the Coma cluster is well reproduced when the spectral index of the turbulent Alfven waves is 2.8. The obtained energy spectrum of the turbulent Alfven waves is steeper than that expected from the turbulence theory, suggesting back reaction of the particle acceleration. The fitting procedure constraints the amplitude of the turbulent Alfven waves as well as the spectral index. Then we estimate the dissipation of the turbulent Alfven waves, which is found to be less than the cooling rate by X-ray radiation. We suggest that the turbulence which is sufficient for particle acceleration is developed in the clusters containing the radio halo. It is most likely that cluster mergers create the turbulence and seed relativistic electrons.

We have calculated evolution of a non-thermal electron population from super-thermal but weakly relativistic to highly relativistic energy range in clusters of galaxies. We investigate evolution of hard X-ray radiation due to both bremsstrahlung and inverse Compton scattering of the cosmic microwave background photons. The bremsstrahlung component is more significant than the inverse Compton scattering one when the momentum spectra of electron sources are steeper than $\sim P_e^{-3.0}$ and vice versa in the case of Coma, where $P_e$ is an electron momentum. The resultant hard X-ray spectra are flatter when the bremsstrahlung component is dominant. When the spectral indices of the source term are in the intermediate range ($ -2.5 \sim -3.5$), too much extreme ultraviolet emission is produced. Inverse Compton dominant models can reproduce Coma cluster results with reasonable injection rates, which are possible in cluster mergers and/or ambient gas accretion.

We have investigated evolution of non-thermal emission from relativistic electrons accelerated at around the shock fronts during merger of clusters of galaxies. We estimate synchrotron radio emission and inverse Compton scattering of cosmic microwave background photons from extreme ultraviolet (EUV) to hard X-ray range. The hard X-ray emission is most luminous in the later stage of merger. Both hard X-ray and radio emissions are luminous only while signatures of merging events are clearly seen in thermal intracluster medium (ICM). On the other hand, EUV radiation is still luminous after the system has relaxed. Propagation of shock waves and bulk-flow motion of ICM play crucial roles to extend radio halos. In the contracting phase, radio halos are located at the hot region of ICM, or between two substructures. In the expanding phase, on the other hand, radio halos are located between two ICM hot regions and shows rather diffuse distribution.

We investigate the dynamical evolution of clusters of galaxies and their observational consequences during off-center mergers, explicitly considering the relaxation process between ions and electrons in intracluster medium by N-body and hydrodynamical simulations. In the contracting phase a bow shock is formed between the two subclusters. The observed temperature between two peaks in this phase depends on the viewing angle even if the geometry of the system seems to be very simple like head-on collisions. Around the most contracting epoch, when we observe merging clusters nearly along the collision axis, they look like spherical relaxed clusters with large temperature gradients. In the expanding phase, spiral bow shocks occur. As in head-on mergers, the electron temperature is significantly lower than the plasma mean one especially in the post-shock regions in the expanding phase. When the systems have relatively large angular momentum, double-peak structures in the X-ray images can survive even after the most contracting epoch. Morphological features in both X-ray images and electron temperature distribution characteristic to off-center mergers are seriously affected by the viewing angle. When the clusters are observed nearly along the collision axis, the distribution of galaxies' line-of-sight (LOS) velocities is a good indicator of mergers. In the contracting phase, an negative kurtosis and a large skewness are expected for nearly equal mass collisions and rather different mass ones, respectively. To obtain statistically significant results, about 1000 galaxies' LOS velocities are required. For nearby clusters ($z<0.05$), large redshift surveys such as 2dF will enable us to study merger dynamics.

A recent observation with the Hipparcos satellite and some numerical simulations imply that the interaction between an oblate halo and a disc is inappropriate for the persistence of galactic warps. Then, we have compared the time evolution of galactic warps in a prolate halo with that in an oblate halo. The haloes were approximated as fixed potentials, while the discs were represented by N-body particles. We have found that the warping in the oblate halo continues to wind up, and finally disappears. On the other hand, for the prolate halo model, the precession rate of the outer disc increases when the precession of the outer disc recedes from that of the inner disc, and vice versa. Consequently, the warping in the prolate halo persisted to the end of the simulation by retaining the alignment of the line of nodes of the warped disc. Therefore, our results suggest that prolate haloes could sustain galactic warps. The physical mechanism of the persistence of warp is discussed on the basis of the torque between a halo and a disc and that between the inner and outer regions of the disc.

We have investigated the effect of pressure from intracluster medium (ICM) on disk galaxies in merging clusters. The ram-pressure on the galaxies rapidly increases when two clusters collide. This leads to stripping of the interstellar medium (ISM) and decrease of star formation rate (SFR) of the galaxies. On the other hand, the increase of SFR caused by compression of ISM is less significant. Thus, cluster merger does not trigger, but weakens star formation activity of the galaxies. In the central region of the colliding clusters, blue galaxies with high velocity should exist, although most of galaxies become red. Following the decrease of blue galaxy fraction in the clusters, the fraction of post-starburst galaxies increases. After merger, many galaxies in the cluster restart star formation activity and the segregation of blue and red galaxies becomes prominent.

We investigate the evolution of intracluster medium during a cluster merger, explicitly considering the relaxation process between the ions and electrons by N-body and hydrodynamical simulations. When two subclusters collide each other, a bow shock is formed between the centers of two substructures and propagate in both directions along the collision axis. The shock primarily heats the ions because the kinetic energy of an ion entering the shock is larger than that of an electron by the ratio of masses. In the post-shock region the energy is transported from the ions to electrons via Coulomb coupling. However, since the energy exchange timescale depends both on the gas density and temperature, distribution of electron temperature becomes more complex than that of the plasma mean temperature, especially in the expanding phase. After the collision of two subclusters, gas outflow occurs not only along the collision axis but also in its perpendicular direction. The gas which is originally located in the central part of the subclusters moves both in the parallel and perpendicular directions. Since the equilibrium timescale of the gas along these directions is relatively short, temperature difference between ions and electrons is larger in the directions tilted by the angles of $\pm 45^\circ$ with respect to the collision axis. The electron temperature could be significantly lower that the plasma mean temperature by $\sim 50 %$ at most. The significance of our results in the interpretation of X-ray observations is briefly discussed.

We investigate evolution of the intracluster medium (ICM), considering the relaxation process between the ions and electrons. According to the standard scenario of structure formation, ICM is heated by the shock in the accretion flow to the gravitational potential well of the dark halo. The shock primarily heats the ions because the kinetic energy of an ion entering the shock is larger than that of an electron by the ratio of masses. Then the electrons and ions exchange the energy through coulomb collisions and reach the equilibrium. From simple order estimation we find that the region where the electron temperature is considerably lower than the ion temperature spreads out on a Mpc scale. We then calculate the ion and electron temperature profiles by combining the adiabatic model of two-temperature plasma by Fox & Loeb (1997) with spherically symmetric N-body and hydrodynamic simulations based on three different cosmological models. It is found that the electron temperature is about a half of the mean temperature at radii $\sim$ 1 Mpc. This could lead to an about 50 % underestimation in the total mass contained within $\sim$ 1 Mpc when the electron temperature profiles are used. The polytropic indices of the electron temperature profiles are $\simeq 1.5$ whereas those of mean temperature $\simeq 1.3$ for $r \geq 1$ Mpc. This result is consistent both with the X-ray observations on electron temperature profiles and with some theoretical and numerical predictions about mean temperature profiles.

Evolutions of X-ray clusters of galaxies are studied by N-body (shell model) + mesh code (TVD) simulations on the assumption of spherical symmetry. We consider a density perturbation of $10^{15} M_{\odot}$ composed of dark matter and gas in cold dark matter dominated universe with the cosmological density parameter, $\Omega_0 = 1$ or $0.2$. A shock front appears during its initial collapse, moving outward as ambient gas accretes towards cluster center. The shock front separates the inner X-ray emitting, hot region, where gas is almost in hydrostatic equilibrium but with small radial infall ($\sim 100$km s$^{-1}$) being left, from the outer cool region. Gas inside the shock is strongly compressed and heated by shock so that X-ray luminosity rapidly rises in the early stage. In the late stage, however, the X-ray luminosity rises only gradually due partly to the expansion of the inner high temperature region and partly to the increase of X-ray emissivity of gas as the results of continuous adiabatic compression inside the shock. The shock front structure, which was not clearly resolved in the previous SPH simulations, is clearly captured by the present simulations. In addition, we find a sound wave propagating outward, thereby producing spatial modulations with amplitudes of $\sim$ 10 % in the radial temperature and density profiles and time variations in the strength of the shock. Such modulations, if observed, could be used as a probe to investigate the structure of clusters.

Violent relaxation process of spherical stellar systems is examined by numerical simulations of shell model. The collapse of uniform density sphere both with and without external force is investigated. It is found that time variation of mean gravitational potential field induced by external force makes the system closer to Lynden-Bell distribution. Our result imply that systems which are gravitationary affected by other systems may approach the Lynden-Bell distribution.