We study shapes and alignments of 45 dark matter (DM) haloes and their brightest cluster galaxies (BCGs) using a sample of 39 massive clusters from Hubble Frontier Field (HFF), Cluster Lensing And Supernova survey with Hubble (CLASH), and Reionization Lensing Cluster Survey (RELICS). We measure shapes of the DM haloes by strong gravitational lensing, whereas BCG shapes are derived from their light profiles in Hubble Space Telescope images. Our measurements from a large sample of massive clusters presented here provide new constraints on dark matter and cluster astrophysics. We find that DM haloes are on average highly elongated with the mean ellipticity of $0.482\pm 0.028$, and position angles of major axes of DM haloes and their BCGs tend to be aligned well with the mean value of alignment angles of $22.2\pm 3.9$ deg. We find that DM haloes in our sample are on average more elongated than their BCGs with the mean difference of their ellipticities of $0.11\pm 0.03$. In contrast, the Horizon-AGN cosmological hydrodynamical simulation predicts on average similar ellipticities between DM haloes and their central galaxies. While such a difference between the observations and the simulation may well be explained by the difference of their halo mass scales, other possibilities include the bias inherent to strong lensing measurements, limited knowledge of baryon physics, or a limitation of cold dark matter.

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 a systematic study of gas density perturbations in cool cores of high-mass galaxy clusters. We select 12 relaxed clusters from the Cluster Lensing And Supernova survey with Hubble (CLASH) sample and analyze their cool core features observed with the Chandra X-ray Observatory. We focus on the X-ray residual image characteristics after subtracting their global profile of the X-ray surface brightness distribution. We find that all the galaxy clusters in our sample have, at least, both one positive and one negative excess regions in the X-ray residual image, indicating the presence of gas density perturbations. We identify and characterize the locally perturbed regions using our detection algorithm, and extract X-ray spectra of the intracluster medium (ICM). The ICM temperature in the positive excess region is lower than that in the negative excess region, whereas the ICM in both regions is in pressure equilibrium in a systematic manner. These results indicate that gas sloshing in cool cores takes place in more than 80% of relaxed clusters (95% CL). We confirm this physical picture by analyzing synthetic X-ray observations of a cool-core cluster from a hydrodynamic simulation, finding that our detection algorithm can accurately extract both the positive and negative excess regions and can reproduce the temperature difference between them. Our findings support the picture that the gas density perturbations are induced by gas sloshing, and a large fraction of cool-core clusters have undergone gas sloshing, indicating that gas sloshing may be capable of suppressing runaway cooling of the ICM.

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.

We present results from simultaneous modeling of high angular resolution GBT/MUSTANG-2 90 GHz Sunyaev-Zel'dovich effect (SZE) measurements and XMM-XXL X-ray images of three rich galaxy clusters selected from the HSC-SSP Survey. The combination of high angular resolution SZE and X-ray imaging enables a spatially resolved multi-component analysis, which is crucial to understand complex distributions of cluster gas properties. The targeted clusters have similar optical richnesses and redshifts, but exhibit different dynamical states in their member galaxy distributions; a single-peaked cluster, a double-peaked cluster, and a cluster belonging to a supercluster. A large-scale residual pattern in both regular Compton-parameter $y$ and X-ray surface brightness distributions is found in the single-peaked cluster, indicating a sloshing mode. The double-peaked cluster shows an X-ray remnant cool core between two SZE peaks associated with galaxy concentrations. The temperatures of the two peaks reach $\sim20-30$ keV in contrast to the cool core component of $\sim2$ keV, indicating a violent merger. The highest SZE signal for the supercluster is elongated along a direction perpendicular to the major axis of the X-ray core, suggesting a minor merger before core passage. The $S_X$ and $y$ distributions are thus perturbed at some level, regardless of the optical properties. We find that an integrated Compton $y$ parameter and a temperature for the major merger are significantly boosted from those expected by the weak-lensing mass and those for the other two clusters show no significant deviations, which is consistent with predictions of numerical simulations.

It is known observationally that the major axes of galaxy clusters and their brightest cluster galaxies are roughly aligned with each other. To understand the origin of the alignment, we identify 40 cluster-sized dark matter (DM) haloes with masses higher than $5\times10^{13}~M_{\odot}$ and their central galaxies (CGs) at $z\approx 0$ in the Horizon-AGN cosmological hydrodynamical simulation. We trace the progenitors at 50 different epochs between $0<z<5$. We then fit their shapes and orientations with a triaxial ellipsoid model. While the orientations of both DM haloes and CGs change significantly due to repeated mergers and mass accretions, their relative orientations are well aligned at each epoch even at high redshifts, $z>1$. The alignment becomes tighter with cosmic time; the major axes of the CGs and their host DM haloes at present are aligned on average within $\sim 30^{\circ}$ in the three dimensional space and $\sim 20^{\circ}$ in the projected plane. The orientations of the major axes of DM haloes on average follow one of the eigen-vectors of the surrounding tidal field that corresponds to the {\it slowest collapsing} (or even stretching) mode, and the alignment with the tidal field also becomes tighter. This implies that the orientations of CGs and DM haloes at the present epoch are largely imprinted in the primordial density field of the Universe, whereas strong dynamical interactions such as mergers are important to explain their mutual alignment at each epoch.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

We present line-of-sight gas sloshing first found in a cool core in a galaxy cluster. The galaxy cluster Abell 907 is identified as a relaxed cluster owing to its global X-ray surface brightness taken by the Chandra X-ray Observatory. The X-ray residual image after removing the global emission of the intracluster medium (ICM), however, shows an arc-like positive excess and a negative excess surrounding the central positive excess in the cluster core, which in turn indicates a disturbance of the ICM. We analyze the X-ray spectra extracted from both regions and find that (1) the ICM temperature and the metal abundance in the positive excess are lower and higher than those in the negative excess, respectively, and (2) the ICM is nearly in pressure equilibrium. We also find a slight redshift difference between the positive and the negative excesses, which corresponds to the velocity shear of $1680^{+1300}_{-920}$ km s$^{-1}$ ($1\sigma$). The X-ray residual image and the ICM properties are consistent with those expected by line-of-sight gas sloshing. Assuming that the gas is moving toward inverse-parallel to each other along the line-of-sight, the shear velocity is expected to be $\sim 800$ km s$^{-1}$. The velocity field of this level is able to provide non-thermal pressure support by $\sim 34\%$ relative to the thermal one. The total kinetic energy inferred from the shear velocity corresponds to $\sim 30\%$ of the bolometric luminosity of the sloshing ICM. Abell 907 is therefore complementary to galaxy clusters in which gas sloshing takes place in the plane of the sky, and is important for understanding gas dynamics driven by sloshing and its influence on the heating to prevent runaway cooling.

High-resolution spectroscopy of the core of the Perseus Cluster of galaxies, using the $Hitomi$ satellite above 2 keV and the $XMM$-$Newton$ Reflection Grating Spectrometer at lower energies, provides reliable constraints on the abundances of O, Ne, Mg, Si, S, Ar, Ca, Cr, Mn, Fe, and Ni. Accounting for all known systematic uncertainties, the Ar/Fe, Ca/Fe, and Ni/Fe ratios are determined with a remarkable precision of less than 10%, while the constraints on Si/Fe, S/Fe, and Cr/Fe are at the 15% level, and Mn/Fe is measured with a 20% uncertainty. The average biases in determining the chemical composition using archival CCD spectra from $XMM$-$Newton$ and $Suzaku$ range typically from 15-40%. A simple model in which the enrichment pattern in the Perseus Cluster core and the proto-solar nebula are identical gives a surprisingly good description of the high-resolution X-ray spectroscopy results, with $\chi^2=10.7$ for 10 d.o.f. However, this pattern is challenging to reproduce with linear combinations of existing supernova nucleosynthesis calculations, particularly given the precise measurements of intermediate $\alpha$-elements enabled by $Hitomi$. We discuss in detail the degeneracies between various supernova progenitor models and explosion mechanisms, and the remaining uncertainties in these theoretical models. We suggest that including neutrino physics in the core-collapse supernova yield calculations may improve the agreement with the observed pattern of $\alpha$-elements in the Perseus Cluster core. Our results provide a complementary benchmark for testing future nucleosynthesis calculations required to understand the origin of chemical elements.

We present results of a search for unidentified line emission and absorption signals in the 2-12keV energy band of spectra extracted from Perseus Cluster core region observations obtained with the 5eV energy resolution Hitomi Soft X-ray Spectrometer. No significant unidentified line emission or absorption is found. Line flux upper limits (1 sigma per resolution element) vary with photon energy and assumed intrinsic width, decreasing from 100 photons cm$^{-2}$ s$^{-1}$ sr $^{-1}$ at 2keV to $<10$ cm$^{-2}$ s$^{-1}$ sr $^{-1}$ over most of the 5-10 keV energy range for a Gaussian line with Doppler broadening of 640 km/s. Limits for narrower and broader lines have a similar energy dependence and are systematically smaller and larger, respectively. These line flux limits are used to constrain the decay rate of hypothetical dark matter candidates. For the sterile neutrino decay rate, new constraints over the the mass range of 4-24 keV with mass resolution better than any previous X-ray analysis are obtained. Additionally, the accuracy of relevant thermal spectral models and atomic data are evaluated. The Perseus cluster spectra may be described by a composite of multi-temperature thermal and AGN power-law continua. Superposed on these, a few line emission signals possibly originating from unmodeled atomic processes (including Si XIV and Fe XXV) are marginally detected and tabulated. Comparisons with previous X-ray upper limits and future prospects for dark matter searches using high-energy resolution spectroscopy are discussed.

We present the results from the Hitomi Soft Gamma-ray Detector (SGD) observation of the Crab nebula. The main part of SGD is a Compton camera, which in addition to being a spectrometer, is capable of measuring polarization of gamma-ray photons. The Crab nebula is one of the brightest X-ray / gamma-ray sources on the sky, and, the only source from which polarized X-ray photons have been detected. SGD observed the Crab nebula during the initial test observation phase of Hitomi. We performed the data analysis of the SGD observation, the SGD background estimation and the SGD Monte Carlo simulations, and, successfully detected polarized gamma-ray emission from the Crab nebula with only about 5 ks exposure time. The obtained polarization fraction of the phase-integrated Crab emission (sum of pulsar and nebula emissions) is (22.1 $\pm$ 10.6)% and, the polarization angle is 110.7$^o$ + 13.2 / $-$13.0$^o$ in the energy range of 60--160 keV (The errors correspond to the 1 sigma deviation). The confidence level of the polarization detection was 99.3%. The polarization angle measured by SGD is about one sigma deviation with the projected spin axis of the pulsar, 124.0$^o$ $\pm$0.1$^o$.

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

While various observations measured ellipticities of galaxy clusters and alignments between orientations of the brightest cluster galaxies and their host clusters, there are only a handful of numerical simulations that implement realistic baryon physics to allow direct comparisons with those observations. Here we investigate ellipticities of galaxy clusters and alignments between various components of them and the central galaxies in the state-of-the-art cosmological hydrodynamical simulation Horizon-AGN, which contains dark matter, stellar, and gas components in a large simulation box of $(100 h^{-1}$ Mpc$)^3$ with high spatial resolution ($\sim1$ kpc). We estimate ellipticities of total matter, dark matter, stellar, gas surface mass density distributions, X-ray surface brightness, and the Compton $y$-parameter of the Sunyaev-Zel'dovich effect, as well as alignments between these components and the central galaxies for 120 projected images of galaxy clusters with masses $M_{200}>5\times10^{13}M_{\odot}$. Our results indicate that the distributions of these components are well aligned with the major-axes of the central galaxies, with the root mean square value of differences of their position angles of $\sim 20^\circ$, which vary little from inner to the outer regions. We also estimate alignments of these various components with total matter distributions, and find tighter alignments than those for central galaxies with the root mean square value of $\sim 15^\circ$. We compare our results with previous observations of ellipticities and position angle alignments and find reasonable agreements. The comprehensive analysis presented in this paper provides useful prior information for analyzing stacked lensing signals as well as designing future observations to study ellipticities and alignments of galaxy clusters.

We present results from the Hitomi X-ray observation of a young composite-type supernova remnant (SNR) G21.5$-$0.9, whose emission is dominated by the pulsar wind nebula (PWN) contribution. The X-ray spectra in the 0.8-80 keV range obtained with the Soft X-ray Spectrometer (SXS), Soft X-ray Imager (SXI) and Hard X-ray Imager (HXI) show a significant break in the continuum as previously found with the NuSTAR observation. After taking into account all known emissions from the SNR other than the PWN itself, we find that the Hitomi spectra can be fitted with a broken power law with photon indices of $\Gamma_1=1.74\pm0.02$ and $\Gamma_2=2.14\pm0.01$ below and above the break at $7.1\pm0.3$ keV, which is significantly lower than the NuSTAR result ($\sim9.0$ keV). The spectral break cannot be reproduced by time-dependent particle injection one-zone spectral energy distribution models, which strongly indicates that a more complex emission model is needed, as suggested by recent theoretical models. We also search for narrow emission or absorption lines with the SXS, and perform a timing analysis of PSR J1833$-$1034 with the HXI and SGD. No significant pulsation is found from the pulsar. However, unexpectedly, narrow absorption line features are detected in the SXS data at 4.2345 keV and 9.296 keV with a significance of 3.65 $\sigma$. While the origin of these features is not understood, their mere detection opens up a new field of research and was only possible with the high resolution, sensitivity and ability to measure extended sources provided by an X-ray microcalorimeter.

The present paper investigates the temperature structure of the X-ray emitting plasma in the core of the Perseus cluster using the 1.8--20.0 keV data obtained with the Soft X-ray Spectrometer (SXS) onboard the Hitomi Observatory. A series of four observations were carried out, with a total effective exposure time of 338 ks and covering a central region $\sim7'$ in diameter. The SXS was operated with an energy resolution of $\sim$5 eV (full width at half maximum) at 5.9 keV. Not only fine structures of K-shell lines in He-like ions but also transitions from higher principal quantum numbers are clearly resolved from Si through Fe. This enables us to perform temperature diagnostics using the line ratios of Si, S, Ar, Ca, and Fe, and to provide the first direct measurement of the excitation temperature and ionization temperature in the Perseus cluster. The observed spectrum is roughly reproduced by a single temperature thermal plasma model in collisional ionization equilibrium, but detailed line ratio diagnostics reveal slight deviations from this approximation. In particular, the data exhibit an apparent trend of increasing ionization temperature with increasing atomic mass, as well as small differences between the ionization and excitation temperatures for Fe, the only element for which both temperatures can be measured. The best-fit two-temperature models suggest a combination of 3 and 5 keV gas, which is consistent with the idea that the observed small deviations from a single temperature approximation are due to the effects of projection of the known radial temperature gradient in the cluster core along the line of sight. Comparison with the Chandra/ACIS and the XMM-Newton/RGS results on the other hand suggests that additional lower-temperature components are present in the ICM but not detectable by Hitomi SXS given its 1.8--20 keV energy band.

The Hitomi SXS spectrum of the Perseus cluster, with $\sim$5 eV resolution in the 2-9 keV band, offers an unprecedented benchmark of the atomic modeling and database for hot collisional plasmas. It reveals both successes and challenges of the current atomic codes. The latest versions of AtomDB/APEC (3.0.8), SPEX (3.03.00), and CHIANTI (8.0) all provide reasonable fits to the broad-band spectrum, and are in close agreement on best-fit temperature, emission measure, and abundances of a few elements such as Ni. For the Fe abundance, the APEC and SPEX measurements differ by 16%, which is 17 times higher than the statistical uncertainty. This is mostly attributed to the differences in adopted collisional excitation and dielectronic recombination rates of the strongest emission lines. We further investigate and compare the sensitivity of the derived physical parameters to the astrophysical source modeling and instrumental effects. The Hitomi results show that an accurate atomic code is as important as the astrophysical modeling and instrumental calibration aspects. Substantial updates of atomic databases and targeted laboratory measurements are needed to get the current codes ready for the data from the next Hitomi-level mission.

We present Hitomi observations of N132D, a young, X-ray bright, O-rich core-collapse supernova remnant in the Large Magellanic Cloud (LMC). Despite a very short observation of only 3.7 ks, the Soft X-ray Spectrometer (SXS) easily detects the line complexes of highly ionized S K and Fe K with 16-17 counts in each. The Fe feature is measured for the first time at high spectral resolution. Based on the plausible assumption that the Fe K emission is dominated by He-like ions, we find that the material responsible for this Fe emission is highly redshifted at ~800 km/s compared to the local LMC interstellar medium (ISM), with a 90% credible interval of 50-1500 km/s if a weakly informative prior is placed on possible line broadening. This indicates (1) that the Fe emission arises from the supernova ejecta, and (2) that these ejecta are highly asymmetric, since no blue-shifted component is found. The S K velocity is consistent with the local LMC ISM, and is likely from swept-up ISM material. These results are consistent with spatial mapping that shows the He-like Fe concentrated in the interior of the remnant and the S tracing the outer shell. The results also show that even with a very small number of counts, direct velocity measurements from Doppler-shifted lines detected in extended objects like supernova remnants are now possible. Thanks to the very low SXS background of ~1 event per spectral resolution element per 100 ks, such results are obtainable during short pointed or slew observations with similar instruments. This highlights the power of high-spectral-resolution imaging observations, and demonstrates the new window that has been opened with Hitomi and will be greatly widened with future missions such as the X-ray Astronomy Recovery Mission (XARM) and Athena.

We report a Hitomi observation of IGR J16318-4848, a high-mass X-ray binary system with an extremely strong absorption of N_H~10^{24} cm^{-2}. Previous X-ray studies revealed that its spectrum is dominated by strong fluorescence lines of Fe as well as continuum emission. For physical and geometrical insight into the nature of the reprocessing material, we utilize the high spectroscopic resolving power of the X-ray microcalorimeter (the soft X-ray spectrometer; SXS) and the wide-band sensitivity by the soft and hard X-ray imager (SXI and HXI) aboard Hitomi. Even though photon counts are limited due to unintended off-axis pointing, the SXS spectrum resolves Fe K{\alpha_1} and K{\alpha_2} lines and puts strong constraints on the line centroid and width. The line width corresponds to the velocity of 160^{+300}_{-70} km s^{-1}. This represents the most accurate, and smallest, width measurement of this line made so far from any X-ray binary, much less than the Doppler broadening and shift expected from speeds which are characteristic of similar systems. Combined with the K-shell edge energy measured by the SXI and HXI spectra, the ionization state of Fe is estimated to be in the range of Fe I--IV. Considering the estimated ionization parameter and the distance between the X-ray source and the absorber, the density and thickness of the materials are estimated. The extraordinarily strong absorption and the absence of a Compton shoulder component is confirmed. These characteristics suggest reprocessing materials which are distributed in a narrow solid angle or scattering primarily with warm free electrons or neutral hydrogen.

The origin of the narrow Fe-K{\alpha} fluorescence line at 6.4 keV from active galactic nuclei has long been under debate; some of the possible sites are the outer accretion disk, the broad line region, a molecular torus, or interstellar/intracluster media. In February-March 2016, we performed the first X-ray microcalorimeter spectroscopy with the Soft X-ray Spectrometer (SXS) onboard the Hitomi satellite of the Fanaroff-Riley type I radio galaxy NGC 1275 at the center of the Perseus cluster of galaxies. With the high energy resolution of ~5 eV at 6 keV achieved by Hitomi/SXS, we detected the Fe-K{\alpha} line with ~5.4 {\sigma} significance. The velocity width is constrained to be 500-1600 km s$^{-1}$ (FWHM for Gaussian models) at 90% confidence. The SXS also constrains the continuum level from the NGC 1275 nucleus up to ~20 keV, giving an equivalent width ~20 eV of the 6.4 keV line. Because the velocity width is narrower than that of broad H{\alpha} line of ~2750 km s$^{-1}$, we can exclude a large contribution to the line flux from the accretion disk and the broad line region. Furthermore, we performed pixel map analyses on the Hitomi/SXS data and image analyses on the Chandra archival data, and revealed that the Fe-K{\alpha} line comes from a region within ~1.6 kpc from the NGC 1275 core, where an active galactic nucleus emission dominates, rather than that from intracluster media. Therefore, we suggest that the source of the Fe-K{\alpha} line from NGC 1275 is likely a low-covering fraction molecular torus or a rotating molecular disk which probably extends from a pc to hundreds pc scale in the active galactic nucleus system.

Extending the earlier measurements reported in Hitomi collaboration (2016, Nature, 535, 117), we examine the atmospheric gas motions within the central 100~kpc of the Perseus cluster using observations obtained with the Hitomi satellite. After correcting for the point spread function of the telescope and using optically thin emission lines, we find that the line-of-sight velocity dispersion of the hot gas is remarkably low and mostly uniform. The velocity dispersion reaches maxima of approximately 200~km~s$^{-1}$ toward the central active galactic nucleus (AGN) and toward the AGN inflated north-western `ghost' bubble. Elsewhere within the observed region, the velocity dispersion appears constant around 100~km~s$^{-1}$. We also detect a velocity gradient with a 100~km~s$^{-1}$ amplitude across the cluster core, consistent with large-scale sloshing of the core gas. If the observed gas motions are isotropic, the kinetic pressure support is less than 10\% of the thermal pressure support in the cluster core. The well-resolved optically thin emission lines have Gaussian shapes, indicating that the turbulent driving scale is likely below 100~kpc, which is consistent with the size of the AGN jet inflated bubbles. We also report the first measurement of the ion temperature in the intracluster medium, which we find to be consistent with the electron temperature. In addition, we present a new measurement of the redshift to the brightest cluster galaxy NGC~1275.

Thanks to its high spectral resolution (~5 eV at 6 keV), the Soft X-ray Spectrometer (SXS) on board Hitomi enables us to measure the detailed structure of spatially resolved emission lines from highly ionized ions in galaxy clusters for the first time. In this series of papers, using the SXS we have measured the velocities of gas motions, metallicities and the multi-temperature structure of the gas in the core of the Perseus cluster. Here, we show that when inferring physical properties from line emissivities in systems like Perseus, the resonant scattering (RS) effect should be taken into account. In the Hitomi waveband, RS mostly affects the FeXXV He$\alpha$ line ($w$) - the strongest line in the spectrum. The flux measured by Hitomi in this line is suppressed by a factor ~1.3 in the inner ~30 kpc, compared to predictions for an optically thin plasma; the suppression decreases with the distance from the center. The $w$ line also appears slightly broader than other lines from the same ion. The observed distortions of the $w$ line flux, shape and distance dependence are all consistent with the expected effect of the resonant scattering in the Perseus core. By measuring the ratio of fluxes in optically thick ($w$) and thin (FeXXV forbidden, He$\beta$, Ly$\alpha$) lines, and comparing these ratios with predictions from Monte Carlo radiative transfer simulations, the velocities of gas motions have been obtained. The results are consistent with the direct measurements of gas velocities from line broadening described elsewhere in this series, although the systematic and statistical uncertainties remain significant. Further improvements in the predictions of line emissivities in plasma models, and deeper observations with future X-ray missions will enable RS measurements to provide powerful constraints on the amplitude and anisotropy of clusters gas motions.

To search for giant X-ray pulses correlated with the giant radio pulses (GRPs) from the Crab pulsar, we performed a simultaneous observation of the Crab pulsar with the X-ray satellite Hitomi in the 2 -- 300 keV band and the Kashima NICT radio observatory in the 1.4 -- 1.7 GHz band with a net exposure of about 2 ks on 25 March 2016, just before the loss of the Hitomi mission.The timing performance of the Hitomi instruments was confirmed to meet the timing requirement and about 1,000 and 100 GRPs were simultaneously observed at the main and inter-pulse phases, respectively, and we found no apparent correlation between the giant radio pulses and the X-ray emission in either the main or inter-pulse phases.All variations are within the 2 sigma fluctuations of the X-ray fluxes at the pulse peaks, and the 3 sigma upper limits of variations of main- or inter- pulse GRPs are 22\% or 80\% of the peak flux in a 0.20 phase width, respectively, in the 2 -- 300 keV band.The values become 25\% or 110\% for main or inter-pulse GRPs, respectively, when the phase width is restricted into the 0.03 phase.Among the upper limits from the Hitomi satellite, those in the 4.5-10 keV and the 70-300 keV are obtained for the first time, and those in other bands are consistent with previous reports.Numerically, the upper limits of main- and inter-pulse GRPs in the 0.20 phase width are about (2.4 and 9.3) $\times 10^{-11}$ erg cm$^{-2}$, respectively. No significant variability in pulse profiles implies that the GRPs originated from a local place within the magnetosphere and the number of photon-emitting particles temporally increases.However, the results do not statistically rule out variations correlated with the GRPs, because the possible X-ray enhancement may appear due to a $>0.02$\% brightening of the pulse-peak flux under such conditions.

We present the results of a blind millimeter line emitter search using ALMA Band 6 data with a single frequency tuning toward four gravitational lensing clusters (RXJ1347.5-1145, Abell S0592, MACS J0416.1-2403, and Abell 2744). We construct three-dimensional signal-to-noise ratio (S/N) cubes with 60 MHz and 100 MHz binning, and search for millimeter line emitters. We do not detect any line emitters with a peak S/N > 5, although we do find a line emitter candidate with a peak S/N ~ 4.5. These results provide upper limits to the CO(3-2), CO(4-3), CO(5-4), and [CII] luminosity functions at z ~ 0.3, 0.7, 1.2, and 6, respectively. Because of the magnification effect of gravitational lensing clusters, the new data provide the first constraints on the CO and [CII] luminosity functions at unprecedentedly low luminosity levels, i.e., down to $<~ 10^{-3}-10^{-1} $Mpc$^{-3}$ dex$^{-1}$ at $L'_{CO} ~ 10^8-10^{10}$ K km $s^{-1}$ pc$^2$ and $<~ 10^{-3}-10^{-2}$ Mpc$^{-3}$ dex$^{-1}$ at $L_[CII] ~ 10^8-10^{10} L_\odot$, respectively. Although the constraints to date are not stringent yet, we find that the evolution of the CO and [CII] luminosity functions are broadly consistent with the predictions of semi-analytical models. This study demonstrates that the wide observations with a single frequency tuning toward gravitational lensing clusters are promising for constraining the CO and [CII] luminosity functions.

The Crab nebula originated from a core-collapse supernova (SN) explosion observed in 1054 A.D. When viewed as a supernova remnant (SNR), it has an anomalously low observed ejecta mass and kinetic energy for an Fe-core collapse SN. Intensive searches were made for a massive shell that solves this discrepancy, but none has been detected. An alternative idea is that the SN1054 is an electron-capture (EC) explosion with a lower explosion energy by an order of magnitude than Fe-core collapse SNe. In the X-rays, imaging searches were performed for the plasma emission from the shell in the Crab outskirts to set a stringent upper limit to the X-ray emitting mass. However, the extreme brightness of the source hampers access to its vicinity. We thus employed spectroscopic technique using the X-ray micro-calorimeter onboard the Hitomi satellite. By exploiting its superb energy resolution, we set an upper limit for emission or absorption features from yet undetected thermal plasma in the 2-12 keV range. We also re-evaluated the existing Chandra and XMM-Newton data. By assembling these results, a new upper limit was obtained for the X-ray plasma mass of <~ 1Mo for a wide range of assumed shell radius, size, and plasma temperature both in and out of the collisional equilibrium. To compare with the observation, we further performed hydrodynamic simulations of the Crab SNR for two SN models (Fe-core versus EC) under two SN environments (uniform ISM versus progenitor wind). We found that the observed mass limit can be compatible with both SN models if the SN environment has a low density of <~ 0.03 cm-3 (Fe core) or <~ 0.1 cm-3 (EC) for the uniform density, or a progenitor wind density somewhat less than that provided by a mass loss rate of 10-5 Mo yr-1 at 20 km s-1 for the wind environment.

High-resolution X-ray spectroscopy with Hitomi was expected to resolve the origin of the faint unidentified E=3.5 keV emission line reported in several low-resolution studies of various massive systems, such as galaxies and clusters, including the Perseus cluster. We have analyzed the Hitomi first-light observation of the Perseus cluster. The emission line expected for Perseus based on the XMM-Newton signal from the large cluster sample under the dark matter decay scenario is too faint to be detectable in the Hitomi data. However, the previously reported 3.5 keV flux from Perseus was anomalously high compared to the sample-based prediction. We find no unidentified line at the reported high flux level. Taking into account the XMM measurement uncertainties for this region, the inconsistency with Hitomi is at a 99% significance for a broad dark-matter line and at 99.7% for a narrow line from the gas. We do not find anomalously high fluxes of the nearby faint K line or the Ar satellite line that were proposed as explanations for the earlier 3.5 keV detections. We do find a hint of a broad excess near the energies of high-n transitions of Sxvi (E=3.44 keV rest-frame) -- a possible signature of charge exchange in the molecular nebula and another proposed explanation for the unidentified line. While its energy is consistent with XMM pn detections, it is unlikely to explain the MOS signal. A confirmation of this interesting feature has to wait for a more sensitive observation with a future calorimeter experiment.

We present the properties of intracluster medium (ICM) in the cool core of the massive cluster of galaxies Abell 1835 obtained with the data by $ Chandra$ $X$-$ray$ $Observatory$. We find distinctive spiral patterns with the radius of 70 kpc (or 18 arcsec) as a whole in the residual image of X-ray surface brightness after the 2-dimensional ellipse model of surface brightness is subtracted. The size is smaller by a factor of 2 -- 4 than that of other clusters known to have a similar pattern. The spiral patterns consist of two arms. One of them appears as positive, and the other does as negative excesses in the residual image. Their X-ray spectra show that the ICM temperatures in the positive- and negative-excess regions are $5.09^{+0.12}_{-0.13}$ keV and $6.52^{+0.18}_{-0.15}$ keV, respectively. In contrast, no significant difference is found in the abundance or pressure, the latter of which suggests that the ICM in the two regions of the spiral patterns is in pressure equilibrium or close. The spatially-resolved X-ray spectroscopy of the central region ($r<40$ arcsec) divided into 92 sub-regions indicates that Abell 1835 is a typical cool core cluster. We also find that the spiral patterns extend from the cool core out to the hotter surrounding ICM. The residual image reveals some lumpy sub-structure in the cool core. The line-of-sight component of the disturbance velocity responsible for the sub-structures is estimated to be lower than 600 km/s. Abell 1835 may be now experiencing an off-axis minor merger.

We characterize the non-sphericity of galaxy clusters by the projected axis ratio of spatial distribution of star, dark matter, and X-ray surface brightness (XSB). We select 40 simulated groups and clusters of galaxies with mass larger than 5E13 Msun from the Horizon simulation that fully incorporates the relevant baryon physics, in particular, the AGN feedback. We find that the baryonic physics around the central region of galaxy clusters significantly affects the non-sphericity of dark matter distribution even beyond the central region, approximately up to the half of the virial radius. Therefore it is very difficult to predict the the probability density function (PDF) of the projected axis ratio of XSB from dark-matter only N-body simulations as attempted in previous studies. Indeed we find that the PDF derived from our simulated clusters exhibits much better agreement with that from the observed X-ray clusters. This indicates that our present methodology to estimate the non-sphericity directly from the Horizon simulation is useful and promising. Further improvements in both numerical modeling and observational data will establish the non-sphericity of clusters as a cosmological test complementary to more conventional statistics based on spherically averaged quantities.

We revisit the non-sphericity of cluster-mass scale halos from cosmological N-body simulation on the basis of triaxial modelling. In order to understand the difference between the simulation results and the conventional ellipsoidal collapse model (EC), we first consider the evolution of individual simulated halos. The major difference between EC and the simulation becomes appreciable after the turn-around epoch. Moreover, it is sensitive to the individual evolution history of each halo. Despite such strong dependence on individual halos, the resulting nonsphericity of halos exhibits weak but robust mass dependence in a statistical fashion; massive halos are more spherical up to the turn-around, but gradually become less spherical by z = 0. This is clearly inconsistent with the EC prediction; massive halos are usually more spherical. In addition, at z=0, inner regions of the halos are less spherical than outer regions, i.e., the density distribution inside the halos is highly inhomogeneous and therefore not self-similar. Since most of previous fitting formulae for the PDF of axis ratio of triaxial ellipsoids have been constructed under the self-similarity assumption, they are not accurate. Indeed, we compute the PDF of projected axis ratio a1/a2 directly from the simulation data without the self-similarity assumption, and find that it is very sensitive to the assumption. The latter needs to be carefully taken into account in direct comparison with observations, and therefore we provide an empirical fitting formula for the PDF of a1/a2. Our preliminary analysis suggests that the derived PDF of a1/a2 roughly agrees with the current weak-lensing observations. More importantly, the present results will be useful in future exploration of the non-sphericity of clusters in X-ray and optical observations.

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

Clusters of galaxies are the most massive gravitationally-bound objects in the Universe and are still forming. They are thus important probes of cosmological parameters and a host of astrophysical processes. Knowledge of the dynamics of the pervasive hot gas, which dominates in mass over stars in a cluster, is a crucial missing ingredient. It can enable new insights into mechanical energy injection by the central supermassive black hole and the use of hydrostatic equilibrium for the determination of cluster masses. X-rays from the core of the Perseus cluster are emitted by the 50 million K diffuse hot plasma filling its gravitational potential well. The Active Galactic Nucleus of the central galaxy NGC1275 is pumping jetted energy into the surrounding intracluster medium, creating buoyant bubbles filled with relativistic plasma. These likely induce motions in the intracluster medium and heat the inner gas preventing runaway radiative cooling; a process known as Active Galactic Nucleus Feedback. Here we report on Hitomi X-ray observations of the Perseus cluster core, which reveal a remarkably quiescent atmosphere where the gas has a line-of-sight velocity dispersion of 164+/-10 km/s in a region 30-60 kpc from the central nucleus. A gradient in the line-of-sight velocity of 150+/-70 km/s is found across the 60 kpc image of the cluster core. Turbulent pressure support in the gas is 4% or less of the thermodynamic pressure, with large scale shear at most doubling that estimate. We infer that total cluster masses determined from hydrostatic equilibrium in the central regions need little correction for turbulent pressure.

The top-hat spherical collapse model (TSC) is one of the most fundamental analytical frameworks to describe the non-linear growth of cosmic structure. TSC has motivated, and been widely applied in, various researches even in the current era of precision cosmology. While numerous studies exist to examine its validity against numerical simulations in a statistical fashion, there are few analyses to compare the TSC dynamics in an individual object-wise basis, which is what we attempt in the present paper. We extract 100 halos at z = 0 from a cosmological N-body simulation according to the conventional TSC criterion for the spherical over-density. Then we trace back their spherical counter-parts at earlier epochs. Just prior to the turn-around epoch of the halos, their dynamics is well approximated by TSC, but their turn-around epochs are systematically delayed and the virial radii are larger by ~ 20 percent on average relative to the TSC predictions. We find that this systematic deviation is mainly ascribed to the non-uniformity/inhomogeneity of dark matter density profiles and the non-zero velocity dispersions, both of which are neglected in TSC. In particular, the inside-out-collapse and shell-crossing of dark matter halos play an important role in generating the significant velocity dispersion. The implications of the present result are briefly discussed.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions developed by the Institute of Space and Astronautical Science (ISAS), with a planned launch in 2015. The ASTRO-H mission is equipped with a suite of sensitive instruments with the highest energy resolution ever achieved at E > 3 keV and a wide energy range spanning four decades in energy from soft X-rays to gamma-rays. The simultaneous broad band pass, coupled with the high spectral resolution of Delta E < 7 eV of the micro-calorimeter, will enable a wide variety of important science themes to be pursued. ASTRO-H is expected to provide breakthrough results in scientific areas as diverse as the large-scale structure of the Universe and its evolution, the behavior of matter in the gravitational strong field regime, the physical conditions in sites of cosmic-ray acceleration, and the distribution of dark matter in galaxy clusters at different redshifts.

The next generation X-ray observatory ASTRO-H will open up a new dimension in the study of galaxy clusters by achieving for the first time the spectral resolution required to measure velocities of the intracluster plasma, and extending at the same time the spectral coverage to energies well beyond 10 keV. This white paper provides an overview of the capabilities of ASTRO-H for exploring gas motions in galaxy clusters including their cosmological implications, the physics of AGN feedback, dynamics of cluster mergers as well as associated high-energy processes, chemical enrichment of the intracluster medium, and the nature of missing baryons and unidentified dark matter.

The Sunyaev-Zel'dovich effect provides a useful probe of cosmology and structure formation in the Universe. Recent years have seen rapid progress in both quality and quantity of its measurements. In this review, we overview cosmological and astrophysical implications of recent and near future observations of the effect. They include measuring the evolution of the cosmic microwave background radiation temperature, the distance-redshift relation out to high redshifts, number counts and power spectra of galaxy clusters, distributions and dynamics of intracluster plasma, and large-scale motions of the Universe.

We present the results from Suzaku observations of the hottest Abell galaxy cluster A2163 at $z=0.2$. To study the physics of gas heating in cluster mergers, we investigated hard X-ray emission from the merging cluster A2163, which hosts the brightest synchrotron radio halo. We analyzed hard X-ray spectra accumulated from two-pointed Suzaku observations. Non-thermal hard X-ray emission should result from the inverse Compton (IC) scattering of relativistic electrons by the CMB photons. To measure this emission, the dominant thermal emission in the hard X-ray band must be modeled in detail. To this end, we analyzed the combined broad-band X-ray data of A2163 collected by Suzaku and XMM-Newton, assuming single- and multi-temperature models for thermal emission and the power-law model for non-thermal emission. From the Suzaku data, we detected significant hard X-ray emission from A2163 in the 12-60 keV band at the $28\sigma$ level (or at the $5.5\sigma$ level if a systematic error is considered). The Suzaku HXD spectrum alone is consistent with the single-T thermal model of gas temperature $kT=14$ keV. From the XMM data, we constructed a multi-T model including a very hot ($kT=18$ keV) component in the NE region. Incorporating the multi-T and the power-law models into a two-component model with a radio-band photon index, the 12-60 keV energy flux of non-thermal emission is constrained within $5.3 \pm 0.9 (\pm 3.8)\times 10^{-12}~{\rm erg\, s^{-1} cm^{-2}}$. The 90% upper limit of detected IC emission is marginal ($< 1.2\times 10^{-11}~{\rm erg\, s^{-1} cm^{-2}}$ in the 12-60 keV). The estimated magnetic field in A2163 is $B > 0.098~{\rm \mu G}$. While the present results represent a three-fold increase in the accuracy of the broad band spectral model of A2163, more sensitive hard X-ray observations are needed to decisively test for the presence of hard X-ray emission due to IC emission.

We examine the validity of the hydrostatic equilibrium (HSE) assumption for galaxy clusters using one of the highest-resolution cosmological hydrodynamical simulations. We define and evaluate several effective mass terms corresponding to the Euler equations of the gas dynamics, and quantify the degree of the validity of HSE in terms of the mass estimate. We find that the mass estimated under the HSE assumption (the HSE mass) deviates from the true mass by up to ~ 30 %. This level of departure from HSE is consistent with the previous claims, but our physical interpretation is rather different. We demonstrate that the inertial term in the Euler equations makes a negligible contribution to the total mass, and the overall gravity of the cluster is balanced by the thermal gas pressure gradient and the gas acceleration term. Indeed the deviation from the HSE mass is well explained by the acceleration term at almost all radii. We also clarify the confusion of previous work due to the inappropriate application of the Jeans equations in considering the validity of HSE from the gas dynamics extracted from cosmological hydrodynamical simulations.

We study the formation of low-mass and extremely metal-poor stars in the early universe. Our study is motivated by the recent discovery of a low-mass (M < 0.8 Msun) and extremely metal-poor (Z <= 4.5 x 10^{-5} Zsun) star in the Galactic halo by Caffau et al. We propose a model that early supernova (SN) explosions trigger the formation of low-mass stars via shell fragmentation. We first perform one-dimensional hydrodynamic simulations of the evolution of an early SN remnant. We show that the shocked shell undergoes efficient radiative cooling and then becomes gravitationally unstable to fragment and collapse in about ten million years. We then follow the thermal evolution of the collapsing fragments using a one-zone code. Our one-zone calculation treats chemistry and radiative cooling self-consistently in low-metallicity gas. The collapsing gas cloud evolves roughly isothermally, until it cools rapidly by dust continuum emission at the density 10^{13}-10^{14} /cc. The cloud core then becomes thermally and gravitationally unstable and fragments. We argue that early SNe can trigger the formation of low-mass stars in the extremely metal-poor environment as Caffau et al. discovered recently.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe via a suite of four instruments, covering a very wide energy range, from 0.3 keV to 600 keV. These instruments include a high-resolution, high-throughput spectrometer sensitive over 0.3-2 keV with high spectral resolution of Delta E < 7 eV, enabled by a micro-calorimeter array located in the focal plane of thin-foil X-ray optics; hard X-ray imaging spectrometers covering 5-80 keV, located in the focal plane of multilayer-coated, focusing hard X-ray mirrors; a wide-field imaging spectrometer sensitive over 0.4-12 keV, with an X-ray CCD camera in the focal plane of a soft X-ray telescope; and a non-focusing Compton-camera type soft gamma-ray detector, sensitive in the 40-600 keV band. The simultaneous broad bandpass, coupled with high spectral resolution, will enable the pursuit of a wide variety of important science themes.

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.

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 Stephan's Quintet (SQ, HCG92) was observed with the Far-Infrared Surveyor (FIS) aboard AKARI in four far-infrared (IR) bands at 65, 90, 140, and 160 um. The AKARI four-band images of the SQ show far-IR emission in the intergalactic medium (IGM) of the SQ. In particular, the 160 um band image shows single peak emission in addition to the structure extending in the North-South direction along the shock ridge as seen in the 140 um band, H2 emission and X-ray emission. Whereas most of the far-IR emission in the shocked region comes from the cold dust component, shock-powered [CII]158um emission can significantly contribute to the emission in the 160 um band that shows a single peak at the shocked region. In the shocked region, the observed gas-to-dust mass ratio is in agreement with the Galactic one. The color temperature of the cold dust component (~20 K) is lower than that in surrounding galaxies (~30 K). We discuss a possible origin of the intergalactic dust emission.

The joint JAXA/NASA ASTRO-H mission is the sixth in a series of highly successful X-ray missions initiated by the Institute of Space and Astronautical Science (ISAS). ASTRO-H will investigate the physics of the high-energy universe by performing high-resolution, high-throughput spectroscopy with moderate angular resolution. ASTRO-H covers very wide energy range from 0.3 keV to 600 keV. ASTRO-H allows a combination of wide band X-ray spectroscopy (5-80 keV) provided by multilayer coating, focusing hard X-ray mirrors and hard X-ray imaging detectors, and high energy-resolution soft X-ray spectroscopy (0.3-12 keV) provided by thin-foil X-ray optics and a micro-calorimeter array. The mission will also carry an X-ray CCD camera as a focal plane detector for a soft X-ray telescope (0.4-12 keV) and a non-focusing soft gamma-ray detector (40-600 keV) . The micro-calorimeter system is developed by an international collaboration led by ISAS/JAXA and NASA. The simultaneous broad bandpass, coupled with high spectral resolution of Delta E ~7 eV provided by the micro-calorimeter will enable a wide variety of important science themes to be pursued.

We perform a uniform, systematic analysis of a sample of 38 X-ray galaxy clusters with three different Chandra calibrations. The temperatures change systematically between calibrations. Cluster temperatures change on average by roughly ~6% for the smallest changes and roughly ~13% for the more extreme changes between calibrations. We explore the effects of the changing cluster spectral properties on Sunyaev-Zel'dovich effect (SZE) and X-ray determinations of the Hubble constant. The Hubble parameter changes by +10% and -13% between the current calibration and two previous Chandra calibrations, indicating that changes in the cluster temperature basically explain the entire change in H_0. Although this work focuses on the difference in spectral properties and resultant Hubble parameters between the calibrations, it is intriguing to note that the newer calibrations favor a lower value of the Hubble constant, H_0 ~ 60 km s-1 Mpc-1, typical of results from SZE/X-ray distances. Both galaxy clusters themselves and the details of the instruments must be known precisely to enable reliable precision cosmology with clusters, which will be feasible with combined efforts from ongoing observations and planned missions and observatories covering a wide range of wavelengths.

Spitzer and AKARI observations have found that polycyclic aromatic hydrocarbons (PAHs) are present in nearby elliptical galaxies, but their spatial distributions are still unknown. In order to investigate their distributions, we performed deep spectral mapping observations of the PAH-detected elliptical galaxy NGC4589, a merger remnant with a minor-axis optical dust lane. As a result, we obtain clear evidence that the PAH 11.3 um emission comes predominantly from the dust lane of the galaxy. We also detect molecular hydrogen line emissions from the dust lane. The PAH 17 um emission is distributed differently from the PAH 11.3 um emission, and more similarly to the dust continuum emission. From their distinctive distributions, we suggest that the PAHs responsible for the 11.3 um feature are secondary products through the evolution of the ISM brought in by the merger.

We present a novel scenario for globular cluster (GC) formation, where the ultraviolet (UV) background radiation effectively works so as to produce compact star clusters. Here, we explore the formation of GCs in UV radiation fields. For this purpose, we calculate baryon and dark matter (DM) dynamics in spherical symmetry, incorporating the self-shielding effects by solving the radiative transfer of UV radiation. In addition, we prescribe the star formation in cooled gas components and pursue the dynamics of formed stars. As a result, we find that the evolution of subgalactic objects in UV background radiation are separated into three types, that is, (1) prompt star formation, where less massive clouds ~10^{5-8} M_sun are promptly self-shielded and undergo star formation, (2) delayed star formation, where photoionized massive clouds >10^8 M_sun collapse despite high thermal pressure and are eventually self-shielded to form stars in a delayed fashion, and (3) supersonic infall, where photoionized less massive clouds ~10^{5-8} M_sun contract with supersonic infall velocity and are self-shielded when a compact core forms. In particular, the type (3) is a novel type found in the present simulations, and eventually produces a very compact star cluster. The resultant mass-to-light ratios, half-mass radii, and velocity dispersions for the three types are compared to the observations of GCs, dwarf spheroidals (dSphs), and ultra-compact dwarfs (UCDs). It turns out that the properties of star clusters resulting from supersonic infall match well with those of observed GCs, whereas the other two types are distinct from GCs. Hence, we conclude that supersonic infall in a UV background is a promising mechanism to form GCs.

We have undertaken a search for the infrared emission from the intracluster dust in the Coma cluster of galaxies by the Multiband Imaging Photometer for Spitzer. Our observations yield the deepest mid and far-infrared images of a galaxy cluster ever achieved. In each of the three bands, we have not detected a signature of the central excess component in contrast to the previous report on the detection by Infrared Space Observatory (ISO). We still find that the brightness ratio between 70 and 160 microns shows a marginal sign of the central excess, in qualitative agreement with the ISO result. Our analysis suggests that the excess ratio is more likely due to faint infrared sources lying on fluctuating cirrus foreground. Our observations yield the 2 sigma upper limits on the excess emission within 100 kpc of the cluster center as 5 x 10^-3 MJy/sr, 6 x 10^-2 MJy/sr, and 7 x 10^-2 MJy/sr, at 24, 70, and 160 microns, respectively. These values are in agreement with those found in other galaxy clusters and suggest that dust is deficient near the cluster center by more than 3 orders of magnitude compared to the interstellar medium.

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.

Deep 1.1 mm continuum observations of 1E0657-56 (the "Bullet Cluster") taken with the millimeter-wavelength camera AzTEC on the 10-m Atacama Submillimeter Telescope Experiment (ASTE), have revealed an extremely bright (S$_{\rm{1.1mm}}=15.9$ mJy) unresolved source. This source, MMJ065837-5557.0, lies close to a maximum in the density of underlying mass-distribution, towards the larger of the two interacting clusters as traced by the weak-lensing analysis of Clowe et al. 2006. Using optical--IR colours we argue that MMJ065837-5557.0 lies at a redshift of $z = 2.7 \pm 0.2$. A lensing-derived mass-model for the Bullet Cluster shows a critical-line (caustic) of magnification within a few arcsecs of the AzTEC source, sufficient to amplify the intrinsic millimeter-wavelength flux of the AzTEC galaxy by a factor of $\gg 20$. After subtraction of the foreground cluster emission at 1.1mm due to the Sunyaev-Zel'dovich effect, and correcting for the magnification, the rest-frame FIR luminosity of MMJ065837-5557.0 is $\le 10^{12} \rm L_{\odot}$, characteristic of a luminous infrared galaxy (LIRG). We explore various scenarios to explain the colors, morphologies and positional offsets between the potential optical and IR counterparts, and their relationship with MMJ065837-5557.0. Until higher-resolution and more sensitive (sub)millimeter observations are available, the detection of background galaxies close to the caustics of massive lensing clusters offers the only opportunity to study this intrinsically faint millimeter-galaxy population.

Our previous analysis indicates that small-scale fluctuations in the intracluster medium (ICM) from cosmological hydrodynamic simulations follow the lognormal distribution. In order to test the lognormal nature of the ICM directly against X-ray observations of galaxy clusters, we develop a method of extracting statistical information about the three-dimensional properties of the fluctuations from the two-dimensional X-ray surface brightness. We first create a set of synthetic clusters with lognormal fluctuations. Performing mock observations of these synthetic clusters, we find that the resulting X-ray surface brightness fluctuations also follow the lognormal distribution fairly well. Systematic analysis of the synthetic clusters provides an empirical relation between the density fluctuations and the X-ray surface brightness. We analyze \chandra observations of the galaxy cluster Abell 3667, and find that its X-ray surface brightness fluctuations follow the lognormal distribution. While the lognormal model was originally motivated by cosmological hydrodynamic simulations, this is the first observational confirmation of the lognormal signature in a real cluster. Finally we check the synthetic cluster results against clusters from cosmological hydrodynamic simulations. As a result of the complex structure exhibited by simulated clusters, the empirical relation shows large scatter. Nevertheless we are able to reproduce the true value of the fluctuation amplitude of simulated clusters within a factor of two from their X-ray surface brightness alone. Our current methodology combined with existing observational data is useful in describing and inferring the statistical properties of the three dimensional inhomogeneity in galaxy clusters.

We performed mid-infrared spectroscopic observations of 18 local dusty elliptical galaxies by using the Infrared Spectrograph (IRS) on board Spitzer. We have significantly detected polycyclic aromatic hydrocarbon (PAH) features from 14 out of the 18 galaxies, and thus found that the presence of PAHs is not rare but rather common in dusty elliptical galaxies. Most of these galaxies show an unusually weak 7.7 um emission feature relative to 11.3 um and 17 um emission features. A large fraction of the galaxies also exhibit H2 rotational line and ionic fine-structure line emissions, which have no significant correlation with the PAH emissions. The PAH features are well correlated with the continuum at 35 um, whereas they are not correlated with the continuum at 6 um. We conclude that the PAH emission of the elliptical galaxies is mostly of interstellar origin rather than of stellar origin, and that the unusual PAH interband strength ratios are likely to be due to a large fraction of neutral to ionized PAHs.

The Hubble constant estimated from the combined analysis of the Sunyaev-Zel'dovich effect and X-ray observations of galaxy clusters is systematically lower than those from other methods by 10-15 percent. We examine the origin of the systematic underestimate using an analytic model of the intracluster medium (ICM), and compare the prediction with idealistic triaxial models and with clusters extracted from cosmological hydrodynamical simulations. We identify three important sources for the systematic errors; density and temperature inhomogeneities in the ICM, departures from isothermality, and asphericity. In particular, the combination of the first two leads to the systematic underestimate of the ICM spectroscopic temperature relative to its emission-weighed one. We find that these three systematics well reproduce both the observed bias and the intrinsic dispersions of the Hubble constant estimated from the Sunyaev-Zel'dovich effect.

We carried out a target survey for Lyman break galaxies (LBGs) and Lyman alpha emitters (LAEs) around QSO SDSS J0211-0009 at z=4.87. The deep and wide broadband and narrowband imaging simultaneously revealed the perspective structure of these two high-z populations. The LBGs without Ly-alpha emission form a filamentary structure including the QSO, while the LAEs are distributed around the QSO but avoid it within a distance of ~4.5Mpc. On the other hand, we serendipitously discovered a protocluster with a significant concentration of LBGs and LAEs where no strongly UV ionizing source such as a QSO or radio galaxy is known to exist. In this cluster field, two populations are spatially cross-correlated with each other. The relative spatial distribution of LAEs to LBGs is in stark contrast between the QSO and the cluster fields. We also found a weak trend showing that the number counts based on Ly-alpha and UV continuum fluxes of LAEs in the QSO field are slightly lower than in the cluster field, whereas the number counts of LBGs are almost consistent with each other. The LAEs avoid the nearby region around the QSO where the local UV background radiation could be ~100 times stronger than the average for the epoch. The clustering segregation between LBGs and LAEs seen in the QSO field could be due to either enhanced early galaxy formation in an overdense environment having caused all the LAEs to evolve into LBGs, or local photoionization due to the strong UV radiation from the QSO effectively causing a deficit in low-mass galaxies like LAEs.

(Abridged) We present the results of three-dimensional radiation-hydrodynamics simulations of the formation and evolution of early HII/HeIII regions around the first stars. Cooling (by H2 and HD) and recollapse of the gas in the relic HII region is also followed in a full cosmological context, until second-generation stars are formed. A large HII region with a few kiloparsec diameter is formed, within which a smaller HeIII region is embedded. Radiative feedback effect quenches further star-formation within the halo for a hundred million years. Accretion onto remnant blackholes will be inefficient. Recombination radiation within the HII region is weak, but persists for 50 million years. We also follow the thermal and chemical evolution of the photo-ionized gas in the relic HII region. The gas cools by HD line cooling down to a few tens Kelvin. At high redshifts (z>10), the minimum gas temperature is limited by T_CMB. Because of its low temperature, the characteristic mass of a Jeans-unstable gas clump is ~ 40 Msun, and is significantly smaller than a typical clump mass for early primordial gas clouds. We find no evidence of fragmentation by this epoch. Together with the small cloud mass, this result indicates that massive, rather than very massive, primordial stars may form in the relic HII region. Such stars might be responsible for early metal-enrichment of the interstellar medium from which recently discovered hyper metal-poor stars were born.

The origin of the recently reported systematic bias in the spectroscopic temperature of galaxy clusters is investigated using cosmological hydrodynamical simulations. We find that the local inhomogeneities of the gas temperature and density, after corrected for the global radial profiles, have nearly a universal distribution that resembles the log-normal function. Based on this log-normal approximation for the fluctuations in the intra-cluster medium, we develop an analytical model that explains the bias in the spectroscopic temperature discovered recently. We conclude that the multi-phase nature of the intra-cluster medium not only from the radial profiles but also from the local inhomogeneities plays an essential role in producing the systematic bias.

We report on mid- and far-IR Spitzer observations of 7 nearby dusty elliptical galaxies by using the Multiband Imaging Photometer (MIPS) and Infrared Spectrograph (IRS). Our sample galaxies are known to contain an excessive amount of interstellar dust against sputtering destruction in hot plasma filling the interstellar space of elliptical galaxies. In order to study the origin and the properties of the excess dust in the hot plasma, we selected galaxies with a wide range of X-ray luminosities but similar optical luminos ities for our Spitzer Guest Observers (GO1) program. The 7 galaxies are detected at the MIPS 24 um, 70 um, and 160 um bands; the far- to mid-IR flux ratios of relatively X-ray-bright elliptical galaxies are lower than those of X-ray-faint galaxies. From the IRS spectra, polycyclic aromatic hydrocarbon (PAH) emission features are detected significantly from 5 of the 7 galaxies; t he emission intensities are weaker as the X-ray luminosity of the galaxy is larger. We have found a correlation between the far- to mid-IR flux ratio and the equivalent width of the PAH emission feature. We have obtained apparent spatial correspondence between mid-IR and X-ray distributions in the outer regions for the three X-ray-brightest galaxies in our sample. Possible interpretations for our observational results are discussed.

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 revisit the estimate of the mass fluctuation amplitude, sigma_8, from the observational X-ray cluster abundance. In particular, we examine the effect of the systematic difference between the cluster virial mass estimated from the X-ray spectroscopy, M_{vir, spec}, and the true virial mass of the corresponding halo, M_{vir}. Mazzotta et al. (2004) recently pointed out the possibility that alpha_M = M_{vir, spec}/M_{vir} is systematically lower than unity. We perform the statistical analysis combining the latest X-ray cluster sample and the improved theoretical models and find that sigma_8 \sim 0.76 +/- 0.01 + 0.50 (1-alpha_M) for 0.5 \le alpha_M \le 1, where the quoted errors are statistical only. Thus if alpha_M \sim 0.7, the value of sigma_8 from cluster abundance alone is now in better agreement with other cosmological data including the cosmic microwave background, the galaxy power spectrum and the weak lensing data. The current study also illustrates the importance of possible systematic effects in mapping real clusters to underlying dark halos which changes the interpretation of cluster abundance statistics.

An observational approach is presented to constrain the global structure and evolution of the intracluster medium based on the ROSAT and ASCA distant cluster sample. From statistical analysis of the gas density profile and the connection to the LX-T relation under the beta-model, the scaled gas profile is nearly universal for the outer region and the LX(>0.2r500) is tightly related to the temperature through T^3 rather than T^2. On the other hand, a large density scatter exists in the core region and there is clearly a deviation from the self-similar scaling for clusters with a small core size. A direct link between the core size and the radiative cooling timescale suggest that t_cool is a parameter to control the gas structure and the appearance of small cores in regular clusters may be much connected with the thermal evolution. We derive the luminosity-ambient temperature (T') relation, assuming the universal temperature profile to find the dispersion around the relation significantly decreases: L_1keV is almost constant for a wide range of t_cool. We further examined the LX-Tbeta and LX-T'beta relations and showed a trend that merging clusters segregate from the regular clusters on the planes. A good correlation between t_cool and the X-ray morphology on the L_1keV-t_cool/t_age plane leads us to define three phases according to the different level of cooling, and draw a phenomenological picture: after a cluster collapses and t_cool falls below t_age, the core cools radiatively with quasi-hydrostatic balancing in the gravitational potential, and the central density gradually becomes higher to evolve from an outer-core-dominant cluster to inner-core-dominant cluster.

For 117 clusters of galaxies, we explore the detectability of intracluster dust grains by current and future infrared facilities, taking into account both collisional heating and sputtering of grains by ambient plasma. If the dust grains are injected into the intergalactic space with the amount and size comparable to the Galactic values, the dust-to-gas ratio is typically $10^{-6}$ and the mean dust temperature is $\sim 30$ K near the cluster center. The predicted infrared intensities lie marginally above the detection thresholds for Spitzer Space Telescope, ASTRO-F, Herschel and SPICA missions. For some nearby clusters such as Perseus, A3571, A2319, A3112 and A2204, good detections of intracluster dust signal are expected in the 70 $\mu$m band. Given rather tight constraints on the dust temperature from observed electron density and temperature, the dust mass can be inferred directly from the infrared observations. Further constraints on the size distribution will be obtained once multi-band data are available. They will definitely provide a powerful probe of the dust injection processes and dust-gas interactions in the intergalactic space.

We study the evolution of supernova (SN) remnants of the first stars, taking proper account of the radiative feedback of the progenitor stars on the surroundings. We carry out a series of one-dimensional hydrodynamic simulations with radiative cooling, starting from initial configurations that are drawn from the results of our earlier radiation hydrodynamic simulations of the first HII regions. In low-mass (< 10^6 M_sun) halos, the stellar radiation significantly reduces the ambient gas density prior to the SN explosion. The blastwave quickly propagates over the halo's virial radius, leading to complete evacuation of the gas even with the input energy of 10^50 erg. We find that a large fraction of the remnant's thermal energy is lost in 0.1-10 Myr by line cooling, whereas, for larger explosion energies, the remnant expands even more rapidly with decreasing interior density, and cools predominantly via inverse Compton process. In higher mass halos, the gas density near the explosion site remains high and the SN shock is heavily confined; the thermal energy of the remnant is quickly radiated away by free-free emission, even if the total input energy exceeds the binding energy of halos by two orders of magnitude. We show that the efficiency of halo destruction is determined not only by the explosion energy but also by the gas density profile, and thus controlled by radiative feedback prior to the explosion. Several implications of our results for the formation of first quasars and second-generation stars in the universe are also discussed.

We present an analysis of Chandra X-ray observations of a compact group of galaxies, HCG 80 (z=0.03). The system is a spiral-only group composed of four late-type galaxies, and has a high-velocity dispersion of 309 km/s. With high-sensitivity Chandra observations, we searched for diffuse X-ray emission from the intragroup medium (IGM); however, no significant emission was detected. We place a severe upper limit on the luminosity of the diffuse gas as LX < 6e40 erg/s. On the other hand, significant emission from three of the four members were detected. In particular, we discovered huge halo emission from HCG 80a that extends on a scale of ~30 kpc perpendicular to the galactic disk, whose X-ray temperature and luminosity were measured to be ~0.6 keV and ~4e40 erg/s in the 0.5-2 keV band, respectively. It is most likely to be an outflow powered by intense starburst activity. Based on the results, we discuss possible reasons for the absence of diffuse X-ray emission in the HCG 80 group, suggesting that the system is subject to galaxy interactions, and is possibly at an early stage of IGM evolution.

We study the formation and evolution of HII regions around the first stars formed at z=10-30. We use a one-dimensional Lagrangian hydrodynamics code which self-consistently incorporates radiative transfer and non-equilibrium primordial gas chemistry. The star-forming region is defined as a spherical gas cloud with a Population III star embedded at the center. We explore a large parameter space of host halo mass, gas density profile, and stellar luminosity. The formation of the HII region is characterized by initial slow expansion of a weak D-type ionization front near the center, followed by rapid propagation of an R-type front throughout the outer gas envelope. We find that the transition between the two front types is indeed a critical condition for the complete ionization of halos of cosmological interest. In small mass (< 10^6 M_sun) halos, the transition takes place within a few 10^5 yr, yielding high escape fractions (>80%) of both ionizing and photodissociating photons. In larger mass (> 10^7 M_sun) halos, the ionization front remains to be of D-type over the lifetime of the massive star, the HII region is confined well inside the virial radius, and the escape fractions are essentially zero. We derive an analytic formula that reproduces well the results of our simulations. We discuss immediate implications of the present results for the star formation history and early reionization of the Universe.

We present a method from an X-ray observation of a galaxy cluster to measure the radial profile of the dark matter velocity dispersion, sigma_DM, and to compare the dark matter ``temperature'' defined as mu m_p sigma_DM**2 / k_B with the gas temperature. The method is applied to the XMM-Newton observation of Abell 1795. The ratio between the specific energy of the dark matter and that of the intracluster medium (ICM), which can be defined as beta_DM in analogy with beta_spec, is found to be less than unity everywhere ranging \~0.3-0.8. In other words, the ICM temperature is higher than the dark matter ``temperature'', even in the central region where the radiative cooling time is short. A beta_DM value smaller than unity can most naturally be explained by heating of the ICM. The excess energy of ICM is estimated to be ~1-3 keV per particle.

Quasi-hydrostatic cooling of the intracluster gas is studied. In the quasi-hydrostatic model, work done by gravity on the inflow gas with dP \neq 0, where P is the gas pressure, is taken into account in the thermal balance. The gas flows in from the outer part so as to compensate the pressure loss of the gas undergoing radiative cooling, but the mass flow is so moderate and smooth that the gas is considered to be quasi-hydrostatic. The temperature of the cooling gas decreases toward the cluster center, but, unlike cooling flows with dP = 0, approaches a constant temperature of \sim 1/3 the temperature of the non-cooling ambient gas. This does not mean that gravitational work cancels out radiative cooling, but means that the temperature of the cooling gas appears to approach a constant value toward the cluster center if the gas maintains the quasi-hydrostatic balance. We discuss the mass flow in quasi-hydrostatic cooling, and compare it with the standard isobaric cooling flow model. We also discuss the implication of \dot{M} for the standard cooling flow model.

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.

We consider non-gravitational heating effects on galaxy clusters on the basis of the Monte-Carlo modeling of merging trees of dark matter halos combined with the thermal evolution of gas inside each halo. Under the assumption of hydrostatic equilibrium and the isothermal gas profiles, our model takes account of the metallicity evolution, metallicity-dependent cooling of gas, supernova energy feedback, and heating due to jets of radio galaxies in a consistent manner. The observed properties of galaxy clusters can be explained in models with higher non-gravitational heating efficiency than that in the conventional model. Possibilities include jet heating by the Fanaroff-Riley Type II radio galaxies, and the enhanced star formation efficiency and/or supernova energy feedback, especially at high redshifts.

We study the H2 cooling emission of forming galaxies, and discuss their observability using the future infrared facility SAFIR. Forming galaxies with mass >10^11 Msun emit most of their gravitational energy liberated by contraction in molecular hydrogen line radiation, although a large part of thermal energy at virialization is radiated away by the H Ly alpha emission. For more massive objects, the degree of heating due to dissipation of kinetic energy is so great that the temperature does not drop below 10^4 K and the gravitational energy is emitted mainly by the Ly alpha emission. Therefore, the total H2 luminosity attains the peak value of about 10^42 ergs/s for forming galaxies whose total mass 10^11 Msun. If these sources are situated at redshift z=8, they can be detected by rotational lines of 0-0S(3) at 9.7 micron and 0-0S(1) at 17 micron by SAFIR. An efficient way to find such H2 emitters is to look at the Ly alpha emitters, since the brightest H2 emitters are also luminous in the Ly alpha emission.

We derive constraints on the mass-temperature relation of galaxy clusters from their observed luminosity-temperature relation and X-ray temperature function. Adopting the isothermal gas in hydrostatic equilibrium embedded in the universal density profile of dark matter halos, we compute the X-ray luminosity for clusters as a function of their hosting halo mass. We find that in order to reproduce the two observational statistics, the mass-temperature relation is very tightly constrained as T_{gas} = (1.5~2.0)keV (M_{vir}/10^{14}h_{70}^{-1}M_\odot)^{0.5~0.55}, and a simple self-similar evolution model (T_{gas} \propto M_{vir}^{2/3}) is strongly disfavored. In the cosmological model that we assume (a \Lambda CDM universe with \Omega_0=0.3, \lambda_0=0.7 and h_{70}=1), the derived mass-temperature relation suggests that the mass fluctuation amplitude \sigma_8 is 0.7--0.8.

We examine the reliability of the merger trees generated for the Monte-Carlo modeling of galaxy formation. In particular we focus on the cold gas fraction predicted from the merger trees with different assumptions on the progenitor distribution function, the timestep, and the mass resolution. We show that the cold gas fraction is sensitive to the accuracy of the merger trees at small-mass scales of progenitors at high redshifts. One can reproduce the Press-Schechter prediction to a reasonable degree by adopting a fairly large number of redshift bins, N_{step} ~ 1000 in generating merger trees, which is a factor of ten larger than the canonical value used in previous literature.

We present a newly measured X-ray temperature function of galaxy clusters using a complete flux-limited sample of 61 clusters. The sample is constructed with the total survey area of 8.14 steradians and the flux limit of 1.99e-11 ergs/s/cm^2 in the 0.1-2.4keV band. X-ray temperatures and fluxes of the sample clusters were accurately measured with ASCA and ROSAT data. The derived temperature function covers an unprecedentedly wide temperature range of 1.4-11keV. By fitting these data with theoretically predicted temperature functions given by the Press-Schechter formalism together with a recent formation approximation and the CDM power spectrum, we obtained tight and individual constraints on Omega_m,0 and sigma_8. We also employed the Formation-Epoch model in which the distribution in the formation epoch of clusters as well as the temperature evolution are taken into account, showing significantly different results. Systematics caused by the uncertainty in the mass-temperature relation are studied and found to be as large as the statistical errors.

We explore possibilities of collapse and star formation in Population III objects exposed to the external ultraviolet background (UVB) radiation. Assuming spherical symmetry, we solve self-consistently radiative transfer of photons, non-equilibrium H2 chemistry, and gas hydrodynamics. Although the UVB does suppress the formation of low mass objects, the negative feedback turns out to be weaker than previously suggested. In particular, the cut-off scale of collapse drops significantly below the virial temperature 10^4 K at weak UV intensities, due to both self-shielding of the gas and H2 cooling. Clouds above this cut-off tend to contract highly dynamically, further promoting self-shielding and H2 formation. For plausible radiation intensities and spectra, the collapsing gas can cool efficiently to temperatures well below 10^4 K before rotationally supported and the final H2 fraction reaches 10^{-3}. Our results imply that star formation can take place in low mass objects collapsing in the UVB. The threshold baryon mass for star formation is \sim 10^9 solar mass for clouds collapsing at redshifts z \simlt 3, but drops significantly at higher redshifts. In a conventional cold dark matter universe, the latter coincides roughly with that of the 1 \sigma density fluctuations. Objects near and above this threshold can thus constitute `building blocks' of luminous structures, and we discuss their links to dwarf spheroidal/elliptical galaxies and faint blue objects. These results suggest that the UVB can play a key role in regulating the star formation history of the Universe.

The SZ effect from clusters of galaxies is a dominant source of secondary CMB anisotropy in the low-redshift universe. We present analytic predictions for the CMB power spectrum from massive halos arising from the SZ effect. Since halos are discrete, the power spectrum consists of a Poisson and a correlation term. The latter is always smaller than the former, which is dominated by nearby bright rich clusters. In practice however, those bright clusters are easy to indentify and can thus be subtracted from the map. After this subtraction, the correlation term dominates degree-scale fluctuations over the Poisson term, as the main contribution to the correlation term comes from distant clusters. We find that the correlation term is detectable by Planck experiment. Since the degree scale spectrum is quite insensitive to the highly uncertain core structures of halos, our predictions are robust on these scales. Measuring the correlation term on degree scales thus cleanly probes the clustering of distant halos. This has not been measured yet, mainly because optical and X-ray surveys are not sufficiently sensitive to include such distant clusters and groups. Our analytic predictions are also compared to adiabatic hydrodynamic simulations. The agreement is remarkably good, down to ten arcminutes scales, indicating that our predictions are robust for the Planck experiment. Below ten arcminute scales, where the details of the core structure dominates the power spectrum, our analytic and simulated predictions might fail. In the near future, interferometer and bolometer array experiments will measure the SZ power spectrum down to arcminutes scales, and yield new insight into the physics of the intrahalo medium.

Over half of GeV gamma-ray sources observed by the EGRET experiment have not yet been identified as known astronomical objects. There is an isotropic component of such unidentified sources, whose number is about 60 in the whole sky. Here we calculate the expected number of dynamically forming clusters of galaxies emitting gamma-rays by high energy electrons accelerated in the shock wave when they form, in the framework of the standard theory of structure formation. We find that a few tens of such forming clusters should be detectable by EGRET and hence a considerable fraction of the isotropic unidentified sources can be accounted for, if about 5% of the shock energy is going into electron acceleration. We argue that these clusters are very difficult to detect in x-ray or optical surveys compared with the conventional clusters, because of their extended angular size of about 1 degree. Hence they define a new population of ``gamma-ray clusters''. If this hypothesis is true, the next generation gamma-ray telescopes such as GLAST will detect more than a few thousands of gamma-ray clusters. It would provide a new tracer of dynamically evolving structures in the universe, in contrast to the x-ray clusters as a tracer of hydrodynamically stabilized systems. We also derive the strength of magnetic field required for the extragalactic gamma-ray background by structure formation to extend up to 100 GeV as observed, that is about 10^{-5} of the shock-heated baryon energy density.

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.

The collapse of marginally Jeans unstable primordial gas clouds in the presence of UV radiation field is discussed. Assuming that the dynamical collapse proceeds approximately in an isothermal self-similar fashion, we investigate the thermal evolution of collapsing central core until H$_2$ cooling dominates photoheating and the temperature drops to below $10^4$K. Consequently, the mass of the cooled core is evaluated as $M_{cool}=3.6\times 10^6 M_\odot(I_{21}/1)^{-0.32}$. This scale depends only on the incident UV intensity, and provides a lower limit to the mass of collapsed objects in the UV radiation field.

To elucidate the effects of the UV background radiation on the collapse of pregalactic clouds, we implement a radiation-hydrodynamical calculation, combining one-dimensional spherical hydrodynamics with an accurate treatment of the radiative transfer of ionizing photons. Both absorption and scattering of UV photons are explicitly taken into account. It turns out that a gas cloud contracting within the dark matter potential does not settle into hydrostatic equilibrium, but undergoes run-away collapse even under the presence of the external UV field. The cloud center is shown to become self-shielded against ionizing photons by radiative transfer effects before shrinking to the rotation barrier. Based on our simulation results, we further discuss the possibility of H2 cooling and subsequent star formation in a run-away collapsing core. The present results are closely relevant to the survival of subgalactic Population III objects as well as to metal injection into intergalactic space.

We have developed a model to describe two-point correlation functions of clusters of galaxies in X-ray flux-limited surveys. Our model properly takes account of nonlinear gravitational evolution of mass fluctuations, redshift-space distortion due to linear peculiar velocity field and to finger-of-god, cluster abundance and bias evolution on the basis of the Press -- Schechter theory, the light-cone effect, and the selection function due to the X-ray flux, temperature and luminosity limits. Applying this model in representative cosmological models, we have presented quantitative predictions for X-ray selected samples feasible from the future surveys with the X-ray satellites including Astro-E, Chandra, and XMM. The comparison of these predictions and the observed cluster clustering will place important cosmological constraints which are complementary to the cluster abundance and the cosmic microwave background.

We present angular power spectra of the cosmic microwave background radiation anisotropy due to fluctuations of the Sunyaev-Zel'dovich (SZ) effect through clusters of galaxies. A contribution from the correlation among clusters is especially focused on, which has been neglected in the previous analyses. Employing the evolving linear bias factor based on the Press-Schechter formalism, we find that the clustering contribution amounts to 20-30% of the Poissonian one at degree angular scales. If we exclude clusters in the local universe, it even exceeds the Poissonian noise, and makes dominant contribution to the angular power spectrum. As a concrete example, we demonstrate the subtraction of the ROSAT X-ray flux-limited cluster samples. It indicates that we should include the clustering effect in the analysis of the SZ fluctuations. We further find that the degree scale spectra essentially depend upon the normalization of the density fluctuations, i.e., \sigma_8, and the gas mass fraction of the cluster, rather than the density parameter of the universe and details of cluster evolution models. Our results show that the SZ fluctuations at the degree scale will provide a possible measure of \sigma_8, while the arc-minute spectra a probe of the cluster evolution. In addition, the clustering spectrum will give us valuable information on the bias at high redshift, if we can detect it by removing X-ray luminous clusters.

The effects of the UV background radiation on the formation of sub-galactic clouds are studied by means of one-dimensional hydrodynamical simulations. The radiative transfer of the ionizing photons due to the absorption by HI, HeI and HeII, neglecting the emission, is explicitly taken into account. We find that the complete suppression of collapse occurs for the clouds with circular velocities typically in the range V_c \sim 15-40 km/s and the 50% reduction in the cooled gas mass with V_c \sim 20-55 km/s. These values depend most sensitively on the collapse epoch of the cloud, the shape of the UV spectrum, and the evolution of the UV intensity. Compared to the optically thin case, previously investigated by Thoul & Weinberg (1996), the absorption of the external UV photon by the intervening medium systematically lowers the above threshold values by \Delta V_c \sim 5 km/s. Whether the gas can contract or keeps expanding is roughly determined by the balance between the gravitational force and the thermal pressure gradient when it is maximally exposed to the external UV flux. Based on our simulation results, we discuss a number of implications on galaxy formation, cosmic star formation history, and the observations of quasar absorption lines. In Appendix, we derive analytical formulae for the photoionization coefficients and heating rates, which incorporate the frequency/direction-dependent transfer of external photons.

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.

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.

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. We pay particular attention to a set of theoretical models which well reproduce the {\it ROSAT} 0.5-2 keV band Log $N$ -- Log $S$, and explore possibilities to break the degeneracy among the viable cosmological models.

We compute the number counts of clusters of galaxies, the logN-logS relation, in several X-ray and submm bands on the basis of the Press-Schechter theory. We pay particular attention to a set of theoretical models which well reproduce the ROSAT 0.5-2 keV band logN-logS, and explore possibilities to further constrain the models from future observations with ASCA and/or at submm bands. The latter is closely related to the European PLANCK mission and the Japanese Large Millimeter and Submillimeter Array (LMSA) project. We exhibit that one can break the degeneracy in an acceptable parameter region on the $\Omega_0 - \sigma_8$ plane by combining the ROSAT logN-logS and the submm number counts. Models which reproduce the ROSAT band logN-logS will have $N(>S) \sim (150-300) (S/10^{-12} erg cm^{-2} s^{-1})^{-1.3}$ str$^{-1}$ at $S > 10^{-12} erg cm^{-2} s^{-1}$ in the ASCA 2-10 keV band, and $N(>S_\nu) \sim (10^2-10^4) (S_\nu/100 mJy)^{-1.5} str^{-1}$ at $S_\nu > 100 mJy$ in the submm (0.85mm) band. The amplitude of the logN-logS is very sensitive to the model parameters in the submm band. We also compute the redshift evolution of the cluster number counts and compare with that of the X-ray brightest Abell-type clusters. The results, although still preliminary, point to low density ($\Omega_0\sim 0.3$) universes. The contribution of clusters to the X-ray and submm background radiations is shown to be insignificant in any model compatible with the ROSAT logN-logS.

We find that the observed log N - log S relation of X-ray clusters can be reproduced remarkably well with a certain range of values for the fluctuation amplitude $\sigma_8$ and the cosmological density parameter $\Omega_0$ in cold dark matter (CDM) universes. The $1\sigma$ confidence limits on $\sigma_8$ in the CDM models with $n=1$ and $h = 0.7$ are expressed as $(0.54 \pm 0.02) \Omega_0^{-0.35-0.82\Omega_0+0.55\Omega_0^2}$ ($\lambda_0=1-\Omega_0$) and $(0.54 \pm 0.02) \Omega_0^{-0.28-0.91\Omega_0+0.68\Omega_0^2}$ ($\lambda_0=0$), where $n$ is the primordial spectral index, and $h$ and $\lambda_0$ are the dimensionless Hubble and cosmological constants. The errors quoted above indicate the statistical ones from the observed log N - log S only, and the systematic uncertainty from our theoretical modelling of X-ray flux in the best-fit value of $\sigma_8$ is about 15%. In the case of $n=1$, we find that the CDM models with $(\Omega_0,\lambda_0,h,\sigma_8) \simeq (0.3,0.7,0.7,1)$ and $(0.45, 0, 0.7, 0.8)$ simultaneously account for the cluster log N - log$S$, X-ray temperature functions, and the normalization from the COBE 4 year data. The derived values assume the observations are without systematic errors, and we discuss in details other theoretical uncertainties which may change the limits on $\Omega_0$ and $\sigma_8$ from the log N - log S relation. We have shown the power of this new approach which will become a strong tool as the observations attain more precision.

Temperature and luminosity functions of X-ray clusters are computed semi-analytically, combining a simple model for the cluster gas properties with the distribution functions of halo formation epochs proposed by Lacey & Cole (1993) and Kitayama & Suto (1996). In contrast to several previous approaches which apply the Press--Schechter mass function in a straightforward manner, our method can explicitly take into account the temperature and luminosity evolution of clusters. In order to make quantitative predictions in a specific cosmological context, we adopt cold dark matter (CDM) universes. Assuming the baryon density parameter $\Omega_{\rm B}=0.0125 h^{-2}$ ($h$ is the Hubble constant in units of 100 km$\cdot$sec$^{-1}\cdot$Mpc$^{-1}$) and the {\it COBE} normalization of matter fluctuations, temperature and luminosity functions of X-ray clusters depend sensitively on the density parameter $\Omega_0$. Allowing for several uncertainties in observational data as well as in our simplified assumptions, we conclude that $\Omega_0 \sim 0.2-0.5$ and $h\sim 0.7$ CDM models with/without the cosmological constant reproduce simultaneously the observed temperature and luminosity functions of X-ray clusters at redshift $z\sim0$.

Analytical expressions for the rates of formation and destruction of gravitationally bound systems are derived assuming that they are originated from primordial random-Gaussian density fluctuations. The resulting formulae reproduce the time derivative of the Press-Schechter mass function in a certain limit. Combining a theoretical model for the evolution of structures with the formation rate, we can make various cosmological predictions which are to be compared with observations. As an example to elucidate such applicability, we evaluate the contribution of clusters of galaxies to the cosmic X-ray background radiation. With the {\it COBE} normalization, we find that the significant fraction of the observed soft X-ray background is accounted for by clusters of galaxies in a cold dark matter universe with $\Omega_0 \sim 0.2$, $\lambda_0 =1- \Omega_0$ and $h \sim 0.8$.