Hyper Suprime-Cam (HSC) is a wide-field imaging camera on the prime focus of the 8.2m Subaru telescope on the summit of Maunakea in Hawaii. A team of scientists from Japan, Taiwan and Princeton University is using HSC to carry out a 300-night multi-band imaging survey of the high-latitude sky. The survey includes three layers: the Wide layer will cover 1400 deg$^2$ in five broad bands ($grizy$), with a $5\,\sigma$ point-source depth of $r \approx 26$. The Deep layer covers a total of 26~deg$^2$ in four fields, going roughly a magnitude fainter, while the UltraDeep layer goes almost a magnitude fainter still in two pointings of HSC (a total of 3.5 deg$^2$). Here we describe the instrument, the science goals of the survey, and the survey strategy and data processing. This paper serves as an introduction to a special issue of the Publications of the Astronomical Society of Japan, which includes a large number of technical and scientific papers describing results from the early phases of this survey.

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

We perform a detailed comparison of the phase-space density traced by the particle distribution in Gadget simulations to the result obtained with a spherical Vlasov solver using the splitting algorithm. The systems considered are apodized H\'enon spheres with two values of the virial ratio, R ~ 0.1 and 0.5. After checking that spherical symmetry is well preserved by the N-body simulations, visual and quantitative comparisons are performed. In particular we introduce new statistics, correlators and entropic estimators, based on the likelihood of whether N-body simulations actually trace randomly the Vlasov phase-space density. When taking into account the limits of both the N-body and the Vlasov codes, namely collective effects due to the particle shot noise in the first case and diffusion and possible nonlinear instabilities due to finite resolution of the phase-space grid in the second case, we find a spectacular agreement between both methods, even in regions of phase-space where nontrivial physical instabilities develop. However, in the colder case, R=0.1, it was not possible to prove actual numerical convergence of the N-body results after a number of dynamical times, even with N=10$^8$ particles.

We develop a detailed methodology of determining three-dimensionally the angle between the stellar spin and the planetary orbit axis vectors, $\psi$, for transiting planetary systems. The determination of $\psi$ requires the independent estimates of the inclination angles of the stellar spin axis and of the planetary orbital axis with respect to the line-of-sight, $i_\star$ and $i_{\rm orb}$, and the projection of the spin--orbit angle onto the plane of the sky, $\lambda$. These are mainly derived from asteroseismology, transit lightcurve and the Rossiter-McLaughlin effect, respectively. The detailed joint analysis of those three datasets enables an accurate and precise determination of the numerous parameters characterizing the planetary system, in addition to $\psi$. We demonstrate the power of the joint analysis for the two specific systems, HAT-P-7 and Kepler-25. HAT-P-7b is the first exoplanet suspected to be a retrograde (or polar) planet because of the significant misalignment $\lambda \approx 180^\circ$. Our joint analysis indicates $i_\star \approx {30D}$ and $\psi \approx 120^\circ$, suggesting that the planetary orbit is closer to polar rather than retrograde. Kepler-25 is one of the few multi-transiting planetary systems with measured $\lambda$, and hosts two short-period transiting planets and one outer non-transiting planet. The projected spin--orbit angle of the larger transiting planet, Kepler-25c, has been measured to be $\lambda \approx 0^\circ$, implying that the system is well-aligned. With the help of the tight constraint from asteroseismology, however, we obtain $i_\star={65.4}^{+{10.6}}_{-{6.4}}$ and $\psi={26.9}^{+{7.0}}_{-{9.2}}$, and thus find that the system is actually mildly misaligned.

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.

Characterizing the surfaces of rocky exoplanets via the scattered light will be an essential challenge to investigate the existence of life on habitable exoplanets. We present a simple reconstruction method for fractional areas of different surface types from photometric variations, or colors, of a second Earth. We create mock light curves for Earth without clouds using empirical data. Then these light curves are fitted to the isotropic scattering model consisting of 4 surface types: ocean, soil, snow and vegetation. In an idealized situation where the photometric errors are only photon shot noise, we are able to reproduce the fractional areas of those components fairly well. We may be even able to detect a signature of vegetation from the distinct feature of photosynthesis on the Earth, known as the red edge. In our reconstruction method, Rayleigh scattering due to the atmosphere has an important effect, and for terrestrial exoplanets with atmosphere similar to our Earth, it is possible to estimate the presence of oceans and an atmosphere simultaneously.

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

We 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 present results from Chandra observations of SDSS J1004+4112, a strongly lensed quasar system with a maximum image separation of 15". All four bright images of the quasar, as well as resolved X-ray emission originating from the lensing cluster, are clearly detected. The emission from the lensing cluster extends out to approximately 1.5 arcmin. We measure the bolometric X-ray luminosity and temperature of the lensing cluster to be 4.7e44 erg s^-1 and 6.4 keV, consistent with the luminosity-temperature relation for distant clusters. The mass estimated from the X-ray observation shows excellent agreement with the mass derived from gravitational lensing. The X-ray flux ratios of the quasar images differ markedly from the optical flux ratios, and the combined X-ray spectrum of the images possesses an unusually strong Fe Kalpha emission line, both of which are indicative of microlensing.

We present measurements of the normalised redshift-space three-point correlation function (Q_z) of galaxies from the Sloan Digital Sky Survey (SDSS) main galaxy sample. We have applied our "npt" algorithm to both a volume-limited (36738 galaxies) and magnitude-limited sample (134741 galaxies) of SDSS galaxies, and find consistent results between the two samples, thus confirming the weak luminosity dependence of Q_z recently seen by other authors. We compare our results to other Q_z measurements in the literature and find it to be consistent within the full jack-knife error estimates. However, we find these errors are significantly increased by the presence of the ``Sloan Great Wall'' (at z ~ 0.08) within these two SDSS datasets, which changes the 3-point correlation function (3PCF) by 70% on large scales (s>=10h^-1 Mpc). If we exclude this supercluster, our observed Q_z is in better agreement with that obtained from the 2dFGRS by other authors, thus demonstrating the sensitivity of these higher-order correlation functions to large-scale structures in the Universe. This analysis highlights that the SDSS datasets used here are not ``fair samples'' of the Universe for the estimation of higher-order clustering statistics and larger volumes are required. We study the shape-dependence of Q_z(s,q,theta) as one expects this measurement to depend on scale if the large scale structure in the Universe has grown via gravitational instability from Gaussian initial conditions. On small scales (s <= 6h^-1 Mpc), we see some evidence for shape-dependence in Q_z, but at present our measurements are consistent with a constant within the errors (Q_z ~ 0.75 +/- 0.05). On scales >10h^-1 Mpc, we see considerable shape-dependence in Q_z.

We study the luminosity and color dependence of the galaxy 2-point correlation function in the Sloan Digital Sky Survey, starting from a sample of 200,000 galaxies over 2500 deg^2. We concentrate on the projected correlation function w(r_p), which is directly related to the real space \xi(r). The amplitude of w(r_p) grows continuously with luminosity, rising more steeply above the characteristic luminosity L_*. Redder galaxies exhibit a higher amplitude and steeper correlation function at all luminosities. The correlation amplitude of blue galaxies increases continuously with luminosity, but the luminosity dependence for red galaxies is less regular, with bright red galaxies more strongly clustered at large scales and faint red galaxies more strongly clustered at small scales. We interpret these results using halo occupation distribution (HOD) models assuming concordance cosmological parameters. For most samples, an HOD model with two adjustable parameters fits the w(r_p) data better than a power-law, explaining inflections at r_p ~ 1-3 Mpc/h as the transition between the 1-halo and 2-halo regimes of \xi(r). The implied minimum mass for a halo hosting a central galaxy above a luminosity threshold L grows as M_min ~ L at low luminosities and more steeply above L_*. The mass at which an average halo has one satellite galaxy brighter than L is M_1 ~ 23 M_min(L). These results imply a conditional luminosity function (at fixed halo mass) in which central galaxies lie far above a Schechter function extrapolation of the satellite population. HOD models nicely explain the joint luminosity-color dependence of w(r_p) in terms of the color fractions of central and satellite populations as a function of halo mass. The inferred HOD properties are in good qualitative agreement with theoretical predictions.

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

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

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

We study the recently discovered gravitational lens SDSS J1004+4112, the first quasar lensed by a cluster of galaxies. It consists of four images with a maximum separation of 14.62''. The system has been confirmed as a lensed quasar at z=1.734 on the basis of deep imaging and spectroscopic follow-up observations. We present color-magnitude relations for galaxies near the lens plus spectroscopy of three central cluster members, which unambiguously confirm that a cluster at z=0.68 is responsible for the large image separation. We find a wide range of lens models consistent with the data, but they suggest four general conclusions: (1) the brightest cluster galaxy and the center of the cluster potential well appear to be offset by several kpc; (2) the cluster mass distribution must be elongated in the North--South direction, which is consistent with the observed distribution of cluster galaxies; (3) the inference of a large tidal shear (~0.2) suggests significant substructure in the cluster; and (4) enormous uncertainty in the predicted time delays between the images means that measuring the delays would greatly improve constraints on the models. We also compute the probability of such large separation lensing in the SDSS quasar sample, on the basis of the CDM model. The lack of large separation lenses in previous surveys and the discovery of one in SDSS together imply a mass fluctuation normalization \sigma_8=1.0^{+0.4}_{-0.2} (95% CL), if cluster dark matter halos have an inner slope -1.5. Shallower profiles would require higher values of \sigma_8. Although the statistical conclusion might be somewhat dependent on the degree of the complexity of the lens potential, the discovery is consistent with the predictions of the abundance of cluster-scale halos in the CDM scenario. (Abridged)

Gravitational lensing is a powerful tool for the study of the distribution of dark matter in the Universe. The cold-dark-matter model of the formation of large-scale structures predicts the existence of quasars gravitationally lensed by concentrations of dark matter so massive that the quasar images would be split by over 7 arcsec. Numerous searches for large-separation lensed quasars have, however, been unsuccessful. All of the roughly 70 lensed quasars known, including the first lensed quasar discovered, have smaller separations that can be explained in terms of galaxy-scale concentrations of baryonic matter. Although gravitationally lensed galaxies with large separations are known, quasars are more useful cosmological probes because of the simplicity of the resulting lens systems. Here we report the discovery of a lensed quasar, SDSS J1004+4112, which has a maximum separation between the components of 14.62 arcsec. Such a large separation means that the lensing object must be dominated by dark matter. Our results are fully consistent with theoretical expectations based on the cold-dark-matter model.

Two-body scattering and other discreteness effects are unimportant in cosmological gravitational clustering in most scenarios, since the dark matter has a small particle mass. The collective field should determine evolution: Two-body scattering in simulations violates the Poisson-Vlasov equations. We test this in PM, P$^3$M, Tree, and NGPM codes, noting that a collisionless code will preserve the one-dimensional character of plane wave collapse. We find collisionality vanishing as the softening parameter approaches the mean interparticle separation. Solutions for the problem are suggested, involving greater computer power, PM-based nested grid codes, and a more conservative approach to resolution claims.

While the recent discovery of the Cepheid variables in the Virgo cluster galaxies puts additional support for the Hubble constant $H_0 \sim 80$km/sec/Mpc, a relatively lower value $H_0 \sim 50$km/sec/Mpc is suggested by other distance indicators based on the Sunyaev-Zel'dovich effect and the gravitational lens which probe the universe at higher redshifts $z=(0.1\sim 1)$. In order to reconcile the possible discrepancy between the estimates of the Hubble constants from nearby galaxy samples and high-redshift clusters, we consider a model of locally open universe embedded in the spatially flat universe. We find analytic expressions for the lower limit on the global Hubble constant $\hg$, and the upper limit on the age of the universe with a given value for the Hubble constant $\hl$ in the local universe. We conclude that it is quite unlikely that the above difference in the estimates of the Hubble constant is explained within the framework of the gravitational instability picture.