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 develop new numerical schemes for Vlasov--Poisson equations with high-order accuracy. Our methods are based on a spatially monotonicity-preserving (MP) scheme and are modified suitably so that positivity of the distribution function is also preserved. We adopt an efficient semi-Lagrangian time integration scheme that is more accurate and computationally less expensive than the three-stage TVD Runge-Kutta integration. We apply our spatially fifth- and seventh-order schemes to a suite of simulations of collisionless self-gravitating systems and electrostatic plasma simulations, including linear and nonlinear Landau damping in one dimension and Vlasov--Poisson simulations in a six-dimensional phase space. The high-order schemes achieve a significantly improved accuracy in comparison with the third-order positive-flux-conserved scheme adopted in our previous study. With the semi-Lagrangian time integration, the computational cost of our high-order schemes does not significantly increase, but remains roughly the same as that of the third-order scheme. Vlasov--Poisson simulations on $128^3 \times 128^3$ mesh grids have been successfully performed on a massively parallel computer.

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.

X-ray and gamma-ray polarimetry is a promising tool to study the geometry and the magnetic configuration of various celestial objects, such as binary black holes or gamma-ray bursts (GRBs). However, statistically significant polarizations have been detected in few of the brightest objects. Even though future polarimeters using X-ray telescopes are expected to observe weak persistent sources, there are no effective approaches to survey transient and serendipitous sources with a wide field of view (FoV). Here we present an electron-tracking Compton camera (ETCC) as a highly-sensitive gamma-ray imaging polarimeter. The ETCC provides powerful background rejection and a high modulation factor over a FoV of up to 2$\pi$ sr thanks to its excellent imaging based on a well-defined point spread function. Importantly, we demonstrated for the first time the stability of the modulation factor under realistic conditions of off-axis incidence and huge backgrounds using the SPring-8 polarized X-ray beam. The measured modulation factor of the ETCC was 0.65 $\pm$ 0.01 at 150 keV for an off-axis incidence with an oblique angle of 30$^\circ$ and was not degraded compared to the 0.58 $\pm$ 0.02 at 130 keV for on-axis incidence. These measured results are consistent with the simulation results. Consequently, we found that the satellite-ETCC proposed in Tanimori et al. (2015) would provide all-sky surveys of weak persistent sources of 13 mCrab with 10% polarization for a 10$^{7}$ s exposure and over 20 GRBs down to a $6\times10^{-6}$ erg cm$^{-2}$ fluence and 10% polarization during a one-year observation.

Since the discovery of nuclear gamma-rays, its imaging has been limited to pseudo imaging, such as Compton Camera (CC) and coded mask. Pseudo imaging does not keep physical information (intensity, or brightness in Optics) along a ray, and thus is capable of no more than qualitative imaging of bright objects. To attain quantitative imaging, cameras that realize geometrical optics is essential, which would be, for nuclear MeV gammas, only possible via complete reconstruction of the Compton process. Recently we have revealed that "Electron Tracking Compton Camera" (ETCC) provides a well-defined Point Spread Function (PSF). The information of an incoming gamma is kept along a ray with the PSF and that is equivalent to geometrical optics. Here we present an imaging-spectroscopic measurement with the ETCC. Our results highlight the intrinsic difficulty with CCs in performing accurate imaging, and show that the ETCC surmounts this problem. The imaging capability also helps the ETCC suppress the noise level dramatically by ~3 orders of magnitude without a shielding structure. Furthermore, full reconstruction of Compton process with the ETCC provides spectra free of Compton edges. These results mark the first proper imaging of nuclear gammas based on the genuine geometrical optics.

In the past several decades, the standard cosmological model has been established and its parameters have been measured to a high precision, while there are still many of the fundamental questions in cosmology; such as the physics in the very early Universe, the origin of the cosmic acceleration and the nature of the dark matter. The future world's largest radio telescope, Square Kilometre Array (SKA), will be able to open the new frontier of cosmology and will be one of the most powerful tools for cosmology in the next decade. The cosmological surveys conducted by the SKA would have the potential not only to answer these fundamental questions but also deliver the precision cosmology. In this article we briefly review the role of the SKA from the view point of the modern cosmology. The cosmology science led by the SKA-Japan Consortium (SKA-JP) Cosmology Science Working Group is also discussed.

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.

Cosmic reionization is known to be a major phase transition of the gas in the Universe. Since astronomical objects formed in the early Universe, such as the first stars, galaxies and black holes, are expected to have caused cosmic reionization, the formation history and properties of such objects are closely related to the reionization process. In spite of the importance of exploring reionization, our understandings regarding reionization is not sufficient yet. Square Kilometre Array (SKA) is a next-generation large telescope that will be operated in the next decade. Although several programs of next-generation telescopes are currently scheduled, the SKA will be the unique telescope with a potential to directly observe neutral hydrogen up to z~30, and provide us with valuable information on the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The early science with the SKA will start in a few years; it is thus the time for us to elaborate a strategy for CD/EoR Science with the SKA. The purpose of this document is to introduce Japanese scientific interests in the SKA project and to report results of our investigation.

We present a new numerical scheme to solve the transfer of diffuse radiation on three-dimensional mesh grids which is efficient on processors with highly parallel architecture such as recently popular GPUs and CPUs with multi- and many-core architectures. The scheme is based on the ray-tracing method and the computational cost is proportional to $N_{\rm m}^{5/3}$ where $N_{\rm m}$ is the number of mesh grids, and is devised to compute the radiation transfer along each light-ray completely in parallel with appropriate grouping of the light-rays. We find that the performance of our scheme scales well with the number of adopted CPU cores and GPUs, and also that our scheme is nicely parallelized on a multi-node system by adopting the multiple wave front scheme, and the performance scales well with the amount of the computational resources. As numerical tests to validate our scheme and to give a physical criterion for the angular resolution of our ray-tracing scheme, we perform several numerical simulations of the photo-ionization of neutral hydrogen gas by ionizing radiation sources without the "on-the-spot" approximation, in which the transfer of diffuse radiation by radiative recombination is incorporated in a self-consistent manner.

(Abridged) We have developed a numerical software library for collisionless N-body simulations named "Phantom-GRAPE" which highly accelerates force calculations among particles by use of a new SIMD instruction set extension to the x86 architecture, AVX, an enhanced version of SSE. In our library, not only the Newton's forces, but also central forces with an arbitrary shape f(r), which has a finite cutoff radius r_cut (i.e. f(r)=0 at r>r_cut), can be quickly computed. Using an Intel Core i7--2600 processor, we measure the performance of our library for both the forces. In the case of Newton's forces, we achieve 2 x 10^9 interactions per second with 1 processor core, which is 20 times higher than the performance of an implementation without any explicit use of SIMD instructions, and 2 times than that with the SSE instructions. With 4 processor cores, we obtain the performance of 8 x 10^9 interactions per second. In the case of the arbitrarily shaped forces, we can calculate 1 x 10^9 and 4 x 10^9 interactions per second with 1 and 4 processor cores, respectively. The performance with 1 processor core is 6 times and 2 times higher than those of the implementations without any use of SIMD instructions and with the SSE instructions. These performances depend weakly on the number of particles. It is good contrast with the fact that the performance of force calculations accelerated by GPUs depends strongly on the number of particles. Substantially weak dependence of the performance on the number of particles is suitable to collisionless N-body simulations, since these simulations are usually performed with sophisticated N-body solvers such as Tree- and TreePM-methods combined with an individual timestep scheme. Collisionless N-body simulations accelerated with our library have significant advantage over those accelerated by GPUs, especially on massively parallel environments.

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 present a scheme for numerical simulations of collisionless self-gravitating systems which directly integrates the Vlasov--Poisson equations in six-dimensional phase space. By the results from a suite of large-scale numerical simulations, we demonstrate that the present scheme can simulate collisionless self-gravitating systems properly. The integration scheme is based on the positive flux conservation method recently developed in plasma physics. We test the accuracy of our code by performing several test calculations including the stability of King spheres, the gravitational instability and the Landau damping. We show that the mass and the energy are accurately conserved for all the test cases we study. The results are in good agreement with linear theory predictions and/or analytic solutions. The distribution function keeps the property of positivity and remains non-oscillatory. The largest simulations are run on 64^6 grids. The computation speed scales well with the number of processors, and thus our code performs efficiently on massively parallel supercomputers.

We present two types of numerical prescriptions that accelerate the radiative transfer calculation around point sources within a three-dimensional Cartesian grid by using the oct-tree structure for the distribution of radiation sources. In one prescription, distant radiation sources are grouped as a bright extended source when the group's angular size, $\theta_{\rm s}$, is smaller than a critical value, $\theta_{\rm crit}$, and radiative transfer is solved on supermeshes whose angular sizes are similar to that of the group of sources. The supermesh structure is constructed by coarse-graining the mesh structure. With this method, the computational time scales with $N_{\rm m} \log(N_{\rm m}) \log(N_{\rm s})$ where $N_{\rm m}$ and $N_{\rm s}$ are the number of meshes and that of radiation sources, respectively. While this method is very efficient, it inevitably overestimates the optical depth when a group of sources acts as an extended powerful radiation source and affects distant meshes. In the other prescription, a distant group of sources is treated as a bright point source ignoring the spatial extent of the group and the radiative transfer is solved on the meshes rather than the supermeshes. This prescription is simply a grid-based version of {\scriptsize START} by Hasegawa & Umemura and yields better results in general with slightly more computational cost ($\propto N_{\rm m}^{4/3} \log(N_{\rm s})$) than the supermesh prescription. Our methods can easily be implemented to any grid-based hydrodynamic codes and are well-suited to the adaptive mesh refinement methods.

We present a high-performance N-body code for self-gravitating collisional systems accelerated with the aid of a new SIMD instruction set extension of the x86 architecture: Advanced Vector eXtensions (AVX), an enhanced version of the Streaming SIMD Extensions (SSE). With one processor core of Intel Core i7-2600 processor (8 MB cache and 3.40 GHz) based on Sandy Bridge micro-architecture, we implemented a fourth-order Hermite scheme with individual timestep scheme (Makino and Aarseth, 1992), and achieved the performance of 20 giga floating point number operations per second (GFLOPS) for double-precision accuracy, which is two times and five times higher than that of the previously developed code implemented with the SSE instructions (Nitadori et al., 2006b), and that of a code implemented without any explicit use of SIMD instructions with the same processor core, respectively. We have parallelized the code by using so-called NINJA scheme (Nitadori et al., 2006a), and achieved 90 GFLOPS for a system containing more than N = 8192 particles with 8 MPI processes on four cores. We expect to achieve about 10 tera FLOPS (TFLOPS) for a self-gravitating collisional system with N 105 on massively parallel systems with at most 800 cores with Sandy Bridge micro-architecture. This performance will be comparable to that of Graphic Processing Unit (GPU) cluster systems, such as the one with about 200 Tesla C1070 GPUs (Spurzem et al., 2010). This paper offers an alternative to collisional N-body simulations with GRAPEs and GPUs.

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.

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.

ORIGIN is a proposal for the M3 mission call of ESA aimed at the study of metal creation from the epoch of cosmic dawn. Using high-spectral resolution in the soft X-ray band, ORIGIN will be able to identify the physical conditions of all abundant elements between C and Ni to red-shifts of z=10, and beyond. The mission will answer questions such as: When were the first metals created? How does the cosmic metal content evolve? Where do most of the metals reside in the Universe? What is the role of metals in structure formation and evolution? To reach out to the early Universe ORIGIN will use Gamma-Ray Bursts (GRBs) to study their local environments in their host galaxies. This requires the capability to slew the satellite in less than a minute to the GRB location. By studying the chemical composition and properties of clusters of galaxies we can extend the range of exploration to lower redshifts (z ~ 0.2). For this task we need a high-resolution spectral imaging instrument with a large field of view. Using the same instrument, we can also study the so far only partially detected baryons in the Warm-Hot Intergalactic Medium (WHIM). The less dense part of the WHIM will be studied using absorption lines at low redshift in the spectra for GRBs.

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.

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 use a series of cosmological N-body simulations and various analytic models to study the evolution of the matter power spectrum in real space in a \Lambda Cold Dark Matter universe. We compare the results of N-body simulations against three analytical model predictions; standard perturbation theory, renormalized perturbation theory, and the closure approximation. We include the effects from finite simulation box size in the comparison. We determine the values of the maximum wavenumbers, k^{lim}_{1%} and k^{lim}_{3%}, below which the analytic models and the simulation results agree to within 1 and 3 percent, respectively. We then provide a simple empirical function which describes the convergence regime determined by comparison between our simulations and the analytical models. We find that if we use the Fourier modes within the convergence regime alone, the characteristic scale of baryon acoustic oscillations can be determined within 1% accuracy from future surveys with a volume of a few h^{-3}Gpc^3 at z\sim1 or z\sim3 in the absence of any systematic distortion of the power spectrum.

We present the results from the analysis of long Suzaku observations of the most X-ray luminous galaxy cluster RX J1347.5-1145 at z=0.451. Aims: We study physical properties of the hot (~20 keV) gas clump in the south-east (SE) region discovered by the Sunyaev-Zel'dovich (SZ) effect observations, to understand the gas physics of a violent cluster merger. We also explore a signature of non-thermal emission using the hard X-ray data. Results: We find that the single-temperature model fails to reproduce the continuum emission and Fe-K lines measured by XIS simultaneously. The two-temperature model with a very hot component improves the fit, although the XIS data can only give a lower bound on its temperature. We detect the hard X-ray emission in the 12-40 keV band at the 7 sigma level; however, the significance becomes marginal when the systematic error in the background estimation is included. With the Suzaku + Chandra joint analysis, we determine the temperature of the SE excess component to be 25.3^{+6.1}_{-4.5} ^{+6.9}_{-9.5} keV (90% statistical and systematic errors), which is in an excellent agreement with the previous SZ + X-ray analysis. This is the first time that the X-ray spectroscopy alone gives a good measurement of the temperature of the hot component in the SE region, which is made possible by Suzaku's unprecedented sensitivity to the wide X-ray band. These results strongly indicate that the cluster has undergone a recent, violent merger. The spectral analysis shows that the SE component is consistent with being thermal. We find the 3 sigma upper limit on the non-thermal flux, F < 8e-12 erg s^{-1} cm^{-2} in the 12-60 keV band. Combining this limit with a recent discovery of the radio mini halo at 1.4 GHz, we find a lower limit on the strength of the intracluster magnetic field, B > 0.007 micro G.

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 discuss the detectability of Warm/Hot Intergalactic medium (WHIM) via the absorption lines toward bright point sources with a future X-ray satellite mission, XEUS. While we consider bright QSOs as specific examples, the methodology can be applied to bright GRB afterglows. We create mock absorption spectra for bright QSOs (more than 20 QSOs over the all sky) using a light-cone output of a cosmological hydrodynamic simulation. We assume that WHIM is under collisional and photo-ionization equilibrium. If WHIM has a constant metallicity of $Z=0.1Z_\odot$, approximately 2 O{\sc vii} absorption line system with $>3\sigma$ will be detected on average along a random line-of-sight toward bright QSOs up to $z=0.3$ for 30 ksec exposure.

We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor), consisting of a 4-stage X-ray telescope and an array of TES microcalorimeters, cooled with mechanical coolers, with a total weight of about 400 kg. The mission will perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines, with the energy coverage up to 1.5 keV. The wide field of view of about 50' diameter, superior energy resolution close to 2 eV FWHM, and very low background will together enable us a wide range of science for diffuse X-ray sources. We briefly describe the design of the satellite, performance of the subsystems and the expected results.

We compare the highly clustered populations of very high redshift galaxies with proto-clusters identified numerically in a standard $\Lambda$CDM universe ($\Omega_0=0.3, \lambda_0=0.7$) simulation. We evolve 256^3 dark matter particles in a comoving box of side 150h^{-1}Mpc. By the present day there are 63 cluster sized objects of mass in excess of 10^{14}h^{-1}Mo in this box. We trace these clusters back to higher redshift finding that their progenitors at z=4--5 are extended regions of typically 20--40 Mpc (comoving) in size, with dark halos of mass in excess of 10^{12}h^{-1}Mo and are overdense by typically 1.3--13 times the cosmological mean density. Comparison with the observation of Lyman alpha emitting (LAEs) galaxies at z=4.86 and at z=4.1 indicates that the observed excess clustering is consistent with that expected for a proto-cluster region if LAEs typically correspond to massive dark halos of more than 10^{12}h^{-1}Mo. We give a brief discussion on the relation between high redshift concentration of massive dark halos and present day rich clusters of galaxies.

Time evolution of the ionization state of metals in the cosmic baryons is investigated in a cosmological context without the assumption of ionization equilibrium. We find that a significant fraction of ionized oxygen ions (OVII and OVIII) in the warm-hot intergalactic medium (WHIM) is not in the ionization equilibrium state at a redshift of z~0. We also investigate the effect on the detectability and observables of WHIM as a consequence of such deviation from ionization equilibrium. It is found that the detectability of WHIM is not altered very much both through its emission and absorption signatures, but line ratios between OVII and OVIII are significantly different from those in the ionization equilibrium state.

We present Particle-Particle-Particle-Mesh (PPPM) and Tree Particle-Mesh (TreePM) implementations on GRAPE-5 and GRAPE-6A systems, special-purpose hardware accelerators for gravitational many-body simulations. In our PPPM and TreePM implementations on GRAPE, the computational time is significantly reduced compared with the conventional implementations without GRAPE, especially under the strong particle clustering, and almost constant irrespective of the degree of particle clustering. We carry out the survey of two simulation parameters, the PM grid spacing and the opening parameter for the most optimal combination of force accuracy and computational speed. We also describe the parallelization of these implementations on a PC-GRAPE cluster, in which each node has one GRAPE board, and present the optimal configuration of simulation parameters for good parallel scalability.

We present an analysis of mock spectral observation of warm-hot intergalactic medium (WHIM) using a constrained simulation of the local universe. The simulated map of oxygen emission lines from local WHIM reproduces well the observed structures traced by galaxies in the real local universe. We further attempt to perform mock observations of outer parts of simulated Coma cluster and A3627 adopting the expected performance of DIOS (Diffuse Intergalactic Oxygen Surveyor), which is proposed as a dedicated soft X-ray mission to search for cosmic missing baryons. We find that WHIMs surrounding nearby clusters are detectable with a typical exposure time of a day, and thus constitute realistic and promising targets for DIOS. We also find that an X-ray emitting clump in front of Coma cluster, recently reported in the XMM-Newton observation, has a counterpart in the simulated local universe, and its observed spectrum can be well reproduced in the simulated local universe if the gas temperature is set to the observationally estimated value.

We present our proposal for a small X-ray mission DIOS (Diffuse Intergalactic Oxygen Surveyor) to perform survey observations of warm-hot intergalactic medium using OVII and OVIII emission lines. This will be proposed to a small satellite program planned by ISAS/JAXA in Japan for a launch in 2008. The instrument consists of an array of TES microcalorimeters with an energy resolution 2 eV, cooled by mechanical coolers. The X-ray telescope will employ 4-stage reflection mirrors with a focal length as short as 70 cm and an angular resolution 2'. In addition to DIOS, we briefly describe the NeXT (New X-ray Telescope) mission, which is a larger Japanese X-ray observatory to be launched in 2010 and plans to explore non-thermal processes in the universe.

Foreseeing the era of high spatial resolution measurements of the Sunyaev-Zel'dovich effect (SZE) in clusters of galaxies, we present a prototype analysis of this sort combined with Chandra X-ray data. It is applied specifically to RX J1347-1145 at z=0.451, the most X-ray-luminous galaxy cluster known, for which the highest resolution SZE and X-ray images are currently available. We demonstrate that the combined analysis yields a unique probe of complex structures in the intracluster medium, offering determinations of their temperature, density, and line-of-sight extent. For a subclump in RX J1347-1145, previously discovered in our SZE map, the temperature inferred after removing the foreground and background components is well in excess of 20 keV, indicating that the cluster has recently undergone a violent merger. Excluding the region around this subclump, the SZE signals in submillimeter to centimeter bands (350, 150, and 21 GHz) are all consistent with those expected from Chandra X-ray observations. We further present a temperature deprojection technique based on the SZE and X-ray images, without any knowledge of spatially resolved X-ray spectroscopy. The methodology presented here will be applicable to a statistical sample of clusters available in the future SZE surveys.

Approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. In order to carry out a direct and homogeneous survey of elusive cosmic missing baryons, we propose a dedicated soft-X-ray mission, {\it DIOS} (Diffuse Intergalactic Oxygen Surveyor). The unprecedented energy resolution ($\sim 2$eV) of the XSA (X-ray Spectrometer Array) on-board {\it DIOS} enables us to identify WHIM with gas temperature $T=10^{6-7}$K and overdensity $\delta=10-100$ located at $z<0.3$ through emission lines of O{\sc vii} and O{\sc viii}. {\it DIOS}, hopefully launched in several years, promises to open a new window of detection and characterization of cosmic missing baryons, and to provide yet another important and complementary tool to trace the large-scale structure of the universe.

We study the angular momentum profiles both for dark matter and for gas within virialized halos, using a statistical sample of halos drawn from cosmological hydrodynamics simulations. Three simulations have been analyzed, one is the ``non-radiative'' simulation, and the other two have radiative cooling. We find that the gas component on average has a larger spin and contains a smaller fraction of mass with negative angular momentum than its dark matter counterpart in the non-radiative model. As to the cooling models, the gas component shares approximately the same spin parameter as its dark matter counterpart, but the hot gas has a higher spin and is more aligned in angular momentum than dark matter, while the opposite holds for the cold gas. After the mass of negative angular momentum is excluded, the angular momentum profile of the hot gas component approximately follows the universal function originally proposed by Bullock et al. for dark matter, though the shape parameter $\mu$ is much larger for hot gas and is comfortably in the range required by observations of disk galaxies. Since disk formation is related to the distribution of hot gas that will cool, our study may explain the fact that the disk component of observed galaxies contains a smaller fraction of low angular momentum material than dark matter in halos.

Most of cosmic baryons predicted by the big-bang nucleosynthesis has evaded the direct detection. Recent numerical simulations indicate that approximately 30 to 50 percent of the total baryons in the present universe is supposed to take a form of warm/hot intergalactic medium (WHIM) whose X-ray continuum emission is very weak. To identify those missing baryons, we consider in detail the detectability of WHIM directly through emission lines of OVII (561, 568, 574, 665eV) and OVIII (653eV). For this purpose, we create mock spectra of the emission lines of WHIM using a light-cone output of the cosmological hydrodynamic simulations. Since the predicted fluxes are generally below the current detection limit, the unambiguous detection requires a dedicated X-ray satellite mission that we also discuss in detail. We find that our proposed mission is sensitive to the WHIM with gas temperature $T=10^{6-7}$K and overdensity $\delta=10-100$ up to a redshift of 0.3 without being significantly contaminated by the cosmic X-ray background and the Galactic emissions. Thus such a mission provides a unique and important tool to identify a large fraction of otherwise elusive baryons in the universe.

We estimated the luminosity function (LF) of IRAS galaxies in the PSCz catalogue. The faint end of the PSCz LF is slightly steeper than that of the LF derived by Saunders et al. (1990; S90). Using an analytical form for the LF used by S90, we obtain the following parameters: \alpha = 1.23 \pm 0.04, L_*=(8.85 \pm 1.75) \times 10^8 h^{-2} L_\odot, \sigma =0.724 \pm 0.010, and \phi_* = (2.34 \pm 0.30) \times 10^{-2} h^3 Mpc^{-3}. We also examined the evolution in the sample by a simple assumption \phi_*(z) \propto (1+z)^P, and found P=3.40 \pm 0.70. It does not affect the three parameters, \alpha, L_*, and \sigma, but \phi_*(z=0) is overestimated up to \sim 15% if we ignore evolution. We estimated the temperature dependence of the LF. The LFs of warm and cool galaxies are quite different: the LF of warm galaxies has a very steep faint end with \alpha =1.37. We also discuss a lump found at the brightest end of the LF.

We present velocity statistics of galaxies and their biases inferred from the statistics of the underlying dark matter using a cosmological hydrodynamic simulation of galaxy formation in low-density and spatially flat cold dark matter cosmogony. We find that the pairwise velocity dispersion (PVD) of all galaxies is significantly lower than that of the dark matter particles, and that the PVD of the young galaxies is lower than that of the old types, and even of all galaxies together, especially at small separations. These results are in reasonable agreement with the recent measurements of PVDs in the Las Campanas redshift survey, the PSCz catalogue and the SDSS data. We also find that the low PVD of young galaxies is due to the effects of dynamical friction as well as the different spatial distribution. We also consider the mean infall velocity and the POTENT density reconstruction that are often used to measure the cosmological parameters, and investigate the effects of spatial bias and dynamical friction. In our simulation, the mean infall velocity of young galaxies is significantly lower than that of all the galaxies or of the old galaxies, and the dynamical bias becomes important on scales less than 3Mpc/h. The mass density field reconstructed from the velocity field of young galaxies using the POTENT-style method suffers in accuracy both from the spatial bias and the dynamical friction on the smoothing scale of R_s=8Mpc/h. On the other hand, in the case of R_s=12Mpc/h, which is typically adopted in the actual POTENT analysis, the density reconstruction based on various tracers of galaxies is reasonably accurate.

We investigate how cluster morphology is affected by the cosmological constant in low-density universes. Using high-resolution cosmological N-body/SPH simulations of flat (\Omega_0 = 0.3, \lambda_0 = 0.7, \Lambda CDM) and open (\Omega_0 = 0.3, \lambda_0 = 0, OCDM) cold dark matter universes, we calculate statistical indicators to quantify the irregularity of the cluster morphologies. We study axial ratios, center shifts, cluster clumpiness, and multipole moment power ratios as indicators for the simulated clusters at z=0 and 0.5. Some of these indicators are calculated for both the X-ray surface brightness and projected mass distributions. In \Lambda CDM all these indicators tend to be larger than those in OCDM at z=0. This result is consistent with the analytical prediction of Richstone, Loeb, & Turner, that is, clusters in \Lambda CDM are formed later than in OCDM, and have more substructure at z=0. We make a Kolmogorov-Smirnov test on each indicator for these two models. We then find that the results for the multipole moment power ratios and the center shifts for the X-ray surface brightness are under the significance level (5%). We results also show that these two cosmological models can be distinguished more clearly at z=0 than z = 0.5 by these indicators.

Before the end of 2002 will be launched the GALEX satellite (a NASA/SMEX project) which will observe all the sky in Ultraviolet (UV) through filters at 1500 and 2300 A down to m(AB) 21. In 2004 will be launched the ASTRO-F satellite which will perform an all sky survey at Far-Infrared (FIR) wavelengths. The cross-correlation of both surveys will lead to very large samples of galaxies for which FIR and UV fluxes will be available. Using the FIR to UV flux ratio as a quantitative tracer of the dust extinction we will be able to measure the extinction in the nearby universe (z<0.2) and to perform a statistically significant analysis of the extinction as a function of galactic properties. Of particular interest is the construction of pure FIR and UV selected samples for which the extinction will be measured as templates for the observation of high redshift galaxies.

Cosmological evolution of the QSO luminosity functions (LFs) at NIR/optical/X-ray bands for 1.3 < z < 3.5 is investigated based on the realistic QSO spectra. The accretion-disk theory predicts that although QSO luminosities only depend on mass-accretion rate, \Mdot, QSO spectra have a dependence on black-hole mass, M_{BH}, as well. The smaller M_{BH} is and/or the larger \Mdot is, the harder becomes the QSO NIR/optical/UV spectrum. We model disk spectra which can reproduce these features and calculated LFs for redshift z ~ 3 with the assumption of new-born QSOs being shining at the Eddington luminosity. The main results are: (i) the observed LFs at optical and X-rays can be simultaneously reproduced. (ii) LFs at optical and X-ray bands are not sensitive to M_{BH}, while LFs at NIR bands are; about one order of magnitude difference is expected in volume number densities at L_{I, J} ~ 10^{46} erg s^{-1} between the case that all QSOs would have the same spectral shape as that of M_{BH} = 10^{9} M_sun and the case with M_{BH} = 10^{11} M_sun. (iii) The resultant LFs at NIR are dominated by 10^{7} M_sun black-holes at L_{I, J} ~ 10^{44} erg s^{-1}, and by 10^{11} M_sun black-holes at L_{I, J} \~ 10^{46} erg s^{-1}. Future infrared observations from space(e.g.NGST) will probe cosmological evolution of black hole masses. For redshift z < 3, on the other hand, the observed optical/X-ray LFs can be fitted, if the initial QSO luminosity L_0 is below the Eddington luminosity. Interestingly, the best fitting values of l = L_0/L_{Edd} are different in B- and X-ray bands; l_B ~ 2.5 l_X. The reason for this discrepancy is briefly discussed.

We perform an extensive analysis of nonlinear and stochastic biasing of galaxies and dark halos in spatially flat low-density CDM universe using cosmological hydrodynamic simulations. We compare their biasing properties with the predictions of an analytic halo biasing model. Dark halos in our simulations exhibit reasonable agreement with the predictions only on scales larger than 10h^{-1}Mpc, and on smaller scales the volume exclusion effect of halos due to their finite size becomes substantial. Interestingly the biasing properties of galaxies are better described by extrapolating the halo biasing model predictions. We also find the clear dependence of galaxy biasing on their formation epoch; the distribution of old populations of galaxies tightly correlates with the underlying mass density field, while that of young populations is slightly more stochastic and anti-biased relative to dark matter. The amplitude of two-point correlation function of old populations is about 3 times larger than that of the young populations. Furthermore, the old population of galaxies reside within massive dark halos while the young galaxies are preferentially formed in smaller dark halos. Assuming that the observed early and late-type galaxies correspond to the simulated old and young populations of galaxies, respectively, all of these segregations of galaxies are consistent with observational ones for the early and late-type of galaxies such as the morphology--density relation of galaxies.

We investigate what we can learn about galaxy formation and evolution from the data which will be obtained by the forthcoming large submillimeter/millimeter facilities, mainly by the Atacama Submillimeter Telescope Experiment (ASTE) and the Atacama Large Millimeter Array/Large Millimeter and Submillimeter Array (ALMA/LMSA). We first calculate the source counts from 350 um to 3 mm using the empirical infrared galaxy number count model of Takeuchi et al. (2001). Based on the number counts, we evaluate the source confusion and determine the confusion limit at various wavebands as a function of the characteristic beam size. At submillimeter wavelengths, source confusion with the 10-15-m class facilities becomes severe at 0.1 to 1 mJy level, and astrometry and flux measurements will be difficult. However, we show that very a large-area survey of submillimeter sources brighter than 10-50 mJy can provide a unique constraint on infrared galaxy evolution at z = 1-2. We also find that the 5\sigma-confusion limit of the LMSA is fainter than 1 \mu Jy. The source counts at such faint flux levels give important information of the epoch of galaxy formation. We then show that multiband photometry from the infrared (by ASTRO-F) to the millimeter can be utilized as a redshift estimator. In addition, we compare the observed 1.4, 5, and 8-GHz source counts with our model counts to examine the contribution of star forming galaxies to the faint radio galaxies.

Recently reported infrared (IR) galaxy number counts and cosmic infrared background (CIRB) all suggest that galaxies have experienced a strong evolution sometime in their lifetime. We statistically estimate the galaxy evolution history from these data. We find that an order of magnitude increase of the far-infrared (FIR) luminosity at redshift z = 0.5 - 1.0 is necessary to reproduce the very high CIRB intensity at 140 um reported by Hauser et al. (1998). z \sim 0.75 and decreases to, even at most, a factor of 10 toward z \sim 5, though many variants are allowed within these constraints. This evolution history also satisfies the constraints from the galaxy number counts obtained by IRAS, ISO and, roughly, SCUBA. The rapid evolution of the comoving IR luminosity density required from the CIRB well reproduces the very steep slope of galaxy number counts obtained by ISO. We also estimate the cosmic star formation history (SFH) from the obtained FIR luminosity density, considering the effect of the metal enrichment in galaxies. The derived SFH increases steeply with redshift in 0 < z < 0.75, and becomes flat or even declines at z > 0.75. This is consistent with the SFH estimated from the reported ultraviolet luminosity density. In addition, we present the performance of the Japanese ASTRO-F FIR galaxy survey. We show the expected number counts in the survey. We also evaluate how large a sky area is necessary to derive a secure information of galaxy evolution up to z \sim 1 from the survey, and find that at least 50 - 300 deg^2 is required.

We report on mapping observations toward the region of the most luminous X-ray cluster RXJ1347-1145 (z=0.45) through the Sunyaev-Zel'dovich effect at 21GHz and 150GHz with the Nobeyama 45-m telescope. While a low angular resolution image at 21GHz (beam-size of 76'') shows a consistent feature with the ROSAT/HRI X-ray image, a higher angular resolution image (13'') at 150GHz reveals complex morphological structures of the cluster region, which cannot be simply described by the spherical isothermal beta-model. If such inhomogeneous morphological features prove to be generic for high redshift clusters, distance measurements to the clusters based on their Sunyaev-Zel'dovich data with low angular resolution imaging should be interpreted with caution.

We studied the statistical methods for the estimation of the luminosity function (LF) of galaxies. We focused on four nonparametric estimators: $1/V_{\rm max}$ estimator, maximum-likelihood estimator of Efstathiou et al. (1988), Cho{\l}oniewski's estimator, and improved Lynden-Bell's estimator. The performance of the $1/V_{\rm max}$ estimator has been recently questioned, especially for the faint-end estimation of the LF. We improved these estimators for the studies of the distant Universe, and examined their performances for various classes of functional forms by Monte Carlo simulations. We also applied these estimation methods to the mock 2dF redshift survey catalog prepared by Cole et al. (1998). We found that $1/V_{\rm max}$ estimator yields a completely unbiased result if there is no inhomogeneity, but is not robust against clusters or voids. This is consistent with the well-known results, and we did not confirm the bias trend of $1/V_{\rm max}$ estimator claimed by Willmer (1997) in the case of homogeneous sample. We also found that the other three maximum-likelihood type estimators are quite robust and give consistent results with each other. In practice we recommend Cho{\l}oniewski's estimator for two reasons: 1. it simultaneously provides the shape and normalization of the LF; 2. it is the fastest among these four estimators, because of the algorithmic simplicity. Then, we analyzed the photometric redshift data of the Hubble Deep Field prepared by Fern\'{a}ndez-Soto et al. (1999) using the above four methods. We also derived luminosity density $\rho_{\rm L}$ at $B$- and $I$-band. Our $B$-band estimation is roughly consistent with that of Sawicki, Lin, & Yee (1997), but a few times lower at $2.0 < z < 3.0$. The evolution of $\rho_{\rm L}(I)$ is found to be less prominent.

We constructed a model of infrared and sub-mm (hereafter IR) galaxy number count and estimated history of the IR luminosity density. We treat the evolutionary change of galaxy luminosities as a stepwise nonparametric form, in order to explore the most suitable evolutionary history which reproduces the present observational results. We found the evolutionary patterns which satisfy both constraints required from Cosmic Infrared Background (CIRB) and IR galaxy number counts. One order of magnitude increase of luminosity at redshift $z=0.75 - 1.0$ was found in IR $60 \mu$m luminosity density evolution. We also found that a large number of galaxies ($\sim 10^7$ in the whole sky) will be detected in all-sky survey at far-infrared by Infrared Imaging Surveyor (IRIS); Japanese infrared satellite project Astro-F.

We present results from a series of cosmological SPH (smoothed particle hydrodynamics) simulations coupled with the P3M (Particle-Particle-Particle-Mesh) solver for the gravitational force. The simulations are designed to predict the statistical properties of X-ray clusters of galaxies as well as to study the formation of galaxies. We have seven simulation runs with different assumptions on the thermal state of the intracluster gas. Following the recent work by Pearce et al., we modify our SPH algorithm so as to phenomenologically incorporate the galaxy formation by decoupling the cooled gas particles from the hot gas particles. All the simulations employ 128^3 particles both for dark matter and for gas components, and thus constitute the largest systematic catalogues of simulated clusters in the SPH method performed so far. These enable us to compare the analytical predictions on statistical properties of X-ray clusters against our direct simulation results in an unbiased manner. We find that the luminosities of the simulated clusters are quite sensitive to the thermal history and also to the numerical resolution of the simulations, and thus are not reliable. On the other hand, the mass-temperature relation for the simulated clusters is fairly insensitive to the assumptions of the thermal state of the intracluster gas, robust against the numerical resolution, and in fact agrees well with the analytic prediction. Therefore the prediction for the X-ray temperature function of clusters on the basis of the Press-Schechter mass function and the virial equilibrium is fairly reliable.

Cosmological implications of clusters of galaxies are discussed with particular attention to their importance in probing the cosmological parameters. More specifically we compute the number counts of clusters of galaxies, Log $N$ -- Log $S$ relation, in X-ray and submm bands on the basis of the Press--Schechter theory. As an important step toward breaking the degeneracy among the viable cosmological models, we observed the most luminous X-ray cluster RXJ1347-1145 in three bands (21 and 43 GHz in the Nobeyama Radio Observatory, Japan, and 350 GHz in the J. C. Maxwell telescope at Mauna Kea, Hawaii). We report on the preliminary results which are in good agreement with the profile of the Sunayev -- Zel'dovich effect predicted on the basis of the previous X-ray observation of the cluster.

We report on the detection of the Sunyaev -- Zel'dovich (SZ) signals toward the most luminous X-ray cluster RXJ1347-1145 at Nobeyama Radio Observatory (21 and 43 GHz) and at James Clerk Maxwell Telescope (350 GHz). In particular the latter is the first successful detection of the SZ temperature increment in the submillimeter band which resolved the profile of a cluster of galaxies. Both the observed spectral dependence and the radial profile of the SZ signals are fully consistent with those expected from the X-ray observation of the cluster. The combined analysis of 21GHz and 350GHz data reproduces the temperature and core-radius of the cluster determined with the ROSAT and ASCA satellites when we adopt the slope of the density profile from the X-ray observations. Therefore our present data provide the strongest and most convincing case for the detection of the submillimeter SZ signal from the cluster, as well as in the Rayleigh -- Jeans regime. We also discuss briefly the cosmological implications of the present results.

We describe a non-parametric method to reconstruct gas density and temperature profiles of galaxy clusters from observations of X-ray surface brightness, emission-weighted temperature and the y-parameter through the Sunyaev-Zel'dovich effect. This method is based on the inversion of the projected profiles under the spherically symmetric approximation, which is independent of the equation of state of cluster gas and does not assume that the gas is in hydrostatic equilibrium. In particular we examine the reliability of the reconstruction method assuming a few theoretical models for cluster gas and assigning the statistical errors expected for future X-ray and radio observations. We also discuss briefly the effect of non-sphericity on the basis of simulated clusters.

We have performed a series of simulations of clusters of galaxies on the basis of the smoothed particle hydrodynamics technique in a spatially-flat cold dark matter universe with $\Omega=0.3$, $\lambda=0.7$, and H_0=70 km/s/Mpc as one of the most successful representative cosmological scenarios. In particular, we focus on the Sunyaev--Zel'dovich effect in submm and mm bands, and estimate the reliability of the estimates of the global Hubble constant $H_0$ and the peculiar velocity of clusters $v_r$. Our simulations indicate that fractional uncertainties of the estimates of $H_0$ amount to $\sim 20%$ mainly due to the departure from the isothermal and spherical gas density distribution. We find a systematic underestimate bias of $H_0$ by $\sim 20%$ for clusters $z\approx 1$, but not at $z\approx 0$. The gas temperature drop in the central regions of our simulated clusters leads to the underestimate bias of $v_r$ by $\sim 5%$ at $z\approx 0$ and by $\sim 15%$ at $z\approx 1$ in addition to the statistical errors of the comparable amount due to the non-spherical gas profile.