Spatial correlations among proto-planetary disk orientations carry unique information on physics of multiple star formation processes. We select five nearby star-forming regions that comprise a number of proto-planetary disks with spatially-resolved images with ALMA and HST, and search for the mutual alignment of the disk axes. Specifically, we apply the Kuiper test to examine the statistical uniformity of the position angle (PA: the angle of the major axis of the projected disk ellipse measured counter-clockwise from the north) distribution. The disks located in the star-forming regions, except the Lupus clouds, do not show any signature of the alignment, supporting the random orientation. Rotational axes of 16 disks with spectroscopic measurement of PA in the Lupus III cloud, a sub-region of the Lupus field, however, exhibit a weak and possible departure from the random distribution at a $2\sigma$ level, and the inclination angles of the 16 disks are not uniform as well. Furthermore, the mean direction of the disk PAs in the Lupus III cloud is parallel to the direction of its filament structure, and approximately perpendicular to the magnetic field direction. We also confirm the robustness of the estimated PAs in the Lupus clouds by comparing the different observations and estimators based on three different methods including sparse modeling. The absence of the significant alignment of the disk orientation is consistent with the turbulent origin of the disk angular momentum. Further observations are required to confirm/falsify the possible disk alignment in the Lupus III cloud.

Direct-imaging techniques of exoplanets have made significant progress recently, and will eventually enable to monitor photometric and spectroscopic signals of earth-like habitable planets in the future. The presence of clouds, however, would remain as one of the most uncertain components in deciphering such direct-imaged signals of planets. We attempt to examine how the planetary obliquity produce different cloud patterns by performing a series of GCM (General Circulation Model) simulation runs using a set of parameters relevant for our Earth. Then we use the simulated photometric lightcurves to compute their frequency modulation due to the planetary spin-orbit coupling over an entire orbital period, and attempt to see to what extent one can estimate the obliquity of an Earth-twin. We find that it is possible to estimate the obliquity of an Earth-twin within the uncertainty of several degrees with a dedicated 4 m space telescope at 10 pc away from the system if the stellar flux is completely blocked. While our conclusion is based on several idealized assumptions, a frequency modulation of a directly-imaged earth-like planet offers a unique methodology to determine its obliquity.

The Earth viewed from outside the Solar system would be identified merely like a pale blue dot, as coined by Carl Sagan. In order to detect possible signatures of the presence of life on a second earth among several terrestrial planets discovered in a habit-able zone, one has to develop and establish a methodology to characterize the planet as something beyond a mere pale blue dot. We pay particular attention to the periodic change of the color of the dot according to the rotation of the planet. Because of the large-scale inhomogeneous distribution of the planetary surface, the reflected light of the dot comprises different color components corresponding to land, ocean, ice, and cloud that cover the surface of the planet. If we decompose the color of the dot into several principle components, in turn, one can identify the presence of the different surface components. Furthermore, the vegetation on the earth is known to share a remarkable reflection signature; the reflection becomes significantly enhanced at wave-lengths longer than 760nm, which is known as a red-edge of the vegetation. If one can identify the corresponding color signature in a pale blue dot, it can be used as a unique probe of the presence of life. I will describe the feasibility of the methodology for future space missions, and consider the direction towards astrobiology from an astrophysicist's point of view.

A number of ongoing surveys are likely to discover star-black hole binaries in our Galaxy in the near future. A fraction of them may be triple systems comprising an inner binary, instead of a single black hole, which might be progenitors of binary black holes (BBHs) routinely discovered now from the gravitational wave. We extend our previous proposal to locate inner BBHs from the short-term radial-velocity (RV) variation of a tertiary star in coplanar triples, and we consider noncoplanar triples and their long-term RV variations as well. Specifically, we assume coplanar and noncoplanar triples with an inner BBH of the total mass $20~M_\odot$, whose outer and inner orbital periods are 80 days and 10 days, respectively. We perform a series of N-body simulations and compare the results with analytic approximate solutions based on quadrupole perturbation theory. For coplanar triples, the pericenter shift of the outer star can be used to detect the hidden inner BBH. For noncoplanar triples, the total RV semi-amplitude of the outer star is modulated periodically on the order of $100$km/s due to its precession over roughly the Kozai-Lidov oscillation timescale. Such long-term modulations would be detectable within a decade, independent of the short-term RV variations on the order of of $100$ m/s at roughly twice the orbital frequency of the inner binary. Thus the RV monitoring of future star-black hole binary candidates offers a promising method for searching for their inner hidden BBHs in optical bands.

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.

A number of protoplanetary disks observed with ALMA potentially provide direct examples of initial conditions for planetary systems. In particular, the HL Tau disk has been intensively studied, and its rings/gaps are conventionally interpreted to be a result of unseen massive planets embedded in the gaps. Based on this interpretation, we carried out N-body simulations to investigate orbital evolution of planets within the protoplanetary disk and after the disk dispersal. Before the disk dispersal, our N-body simulations include both migration and mass-growth of the planet coupled with evolution of the disk. By varying the disk parameters, we produce a variety of widely-separated planetary systems consisting of three super-Jupiters at the end of disk dispersal. We found the outer planet is more massive than the inner one, and the migration of the innermost planet is inefficient due to the accretion of outer planet(s). We also showed how the final configuration and the final planetary mass depend on disk parameters. The migration is found to be convergent and no planet-pair has a period ratio less than 2. After the disk dispersal, we switch to pure gravitational N-body simulations and integrate the orbits up to 10 Gyr. Most simulated systems remain stable for at least 10 Gyr. We discuss implications of our result in terms of the observed widely-separated planetary systems HR 8799 and PDS 70.

There are several on-going projects to search for stars orbiting around an invisible companion. A fraction of such candidates may be a triple, instead of a binary, consisting of an inner binary black hole (BBH) and an outer orbiting star. In this paper, we propose a methodology to search for a signature of such an inner BBH, possibly a progenitor of gravitational-wave sources discovered by {\it LIGO}, from the precise radial velocity (RV) follow-up of the outer star. We first describe a methodology using an existing approximate RV formula for coplanar circular triples. We apply this method and constrain the parameters of a possible inner binary objects in 2M05215658+4359220, which consists of a red giant and an unseen companion. Next we consider co-planar but non-circular triples. We compute numerically the RV variation of a tertiary star orbiting around an inner BBH, generate mock RV curves, and examine the feasibility of the BBH detection for our fiducial models. We conclude that the short-cadence RV monitoring of a star-BH binary provides an interesting and realistic method to constrain and/or search for possible inner BBHs. Indeed a recent discovery of a star--BH binary system LB-1 may imply that there are a large number of such unknown objects in our Galaxy, which are ideal targets for the methodology proposed here.

We perform a series of three-dimensional smoothed particle hydrodynamics (SPH) simulations to study the evolution of the angle between the protostellar spin and the protoplanetary disk rotation axes (the star-disk angle $\psi_{\rm sd}$) in turbulent molecular cloud cores. While $\psi_{\rm sd}$ at the protostar formation epoch exhibits broad distribution up to $\sim 130^{\circ}$, $\psi_{\rm sd}$ decreases ($\lesssim 20^{\circ}$) in a timescale of $\sim 10^{4}$ yr. This timescale of the star-disk alignment, $t_{\rm alignment}$, corresponds basically to the mass doubling time of the central protostar, in which the protostar forgets its initial spin direction due to the mass accretion from the disk. Values of $\psi_{\rm sd}$ both at $t=10^2$ yr and $t=10^5$ yr after the protostar formation are independent of the ratios of thermal and turbulent energies to gravitational energy of the initial cloud cores: $\alpha=E_{\rm thermal}/|E_{\rm gravity}|$ and $\gamma_{\rm turb}=E_{\rm turbulence}/|E_{\rm gravity}|$. We also find that a warped disk is possibly formed by the turbulent accretion flow from the circumstellar envelope.

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.

This paper presents the second data release of the Hyper Suprime-Cam Subaru Strategic Program, a wide-field optical imaging survey on the 8.2 meter Subaru Telescope. The release includes data from 174 nights of observation through January 2018. The Wide layer data cover about 300 deg^2 in all five broadband filters (grizy) to the nominal survey exposure (10min in gr and 20min in izy). Partially observed areas are also included in the release; about 1100 deg^2 is observed in at least one filter and one exposure. The median seeing in the i-band is 0.6 arcsec, demonstrating the superb image quality of the survey. The Deep (26 deg^2) and UltraDeep (4 deg^2) data are jointly processed and the UltraDeep-COSMOS field reaches an unprecedented depth of i~28 at 5 sigma for point sources. In addition to the broad-bands, narrow-band data are also available in the Deep and UltraDeep fields. This release includes a major update to the processing pipeline, including improved sky subtraction, PSF modeling, object detection, and artifact rejection. The overall data quality has been improved, but this release is not without problems; there is a persistent deblender problem as well as new issues with masks around bright stars. The user is encouraged to review the issue list before utilizing the data for scientific explorations. All the image products as well as catalog products are available for download. The catalogs are also loaded to a database, which provides an easy interface for users to retrieve data for objects of interest. In addition to these main data products, detailed galaxy shape measurements withheld from the Public Data Release 1 (PDR1) are now available to the community. The shape catalog is drawn from the S16A internal release, which has a larger area than PDR1 (160 deg^2). All products are available at the data release site, https://hsc-release.mtk.nao.ac.jp/.

We measure the rotation periods of 19 stars in the {\it Kepler} transiting planetary systems, $P_{\rm rot, astero}$ from asteroseismology and $P_{\rm rot, phot}$ from photometric variation of their lightcurve. Two stars exhibit two clear peaks in the Lomb-Scargle periodogram, neither of which agrees with the seismic rotation period. Other four systems do not show any clear peak, whose stellar rotation period is impossible to estimate reliably from the photometric variation; their stellar equators may be significantly inclined with respect to the planetary orbital plane. For the remaining 13 systems, $P_{\rm rot, astero}$ and $P_{\rm rot, phot}$ agree within 30\%. Interestingly, three out of the 13 systems are in the spin-orbit resonant state in which $P_{\rm orb, b}/P_{\rm rot, astero} \approx 1$ with $P_{\rm orb, b}$ being the orbital period of the inner-most planet of each system. The corresponding chance probability is ($0.2$-$4.7$) \% based on the photometric rotation period data for 464 {\it Kepler} transiting planetary systems. While further analysis of stars with reliable rotation periods is required to examine the statistical significance, the spin-orbit resonance between the star and planets, if confirmed, have important implications for the star-planet tidal interaction, in addition to the origin of the spin-orbit (mis-)alignment of transiting planetary systems.

Kepler-408 is one of the 33 planet-hosting {\it Kepler} stars for which asteroseismology has been used to investigate the orientation of the stellar rotation axis relative to the planetary orbital plane. The transiting "hot Earth," Kepler-408b, has an orbital period of 2.5 days and a radius of $0.86$~$R_\oplus$, making it much smaller than the planets for which spin-orbit alignment has been studied using the Rossiter-McLaughlin effect. Because conflicting asteroseismic results have been reported in the literature, we undertake a thorough re-appraisal of this system and perform numerous checks for consistency and robustness. We find that the conflicting results are due to the different models for the low-frequency noise in the power spectrum. A careful treatment of the background noise resolves these conflicts, and shows that the stellar inclination is $\is=42^{+5}_{-4}$ degrees. Kepler-408b is, by far, the smallest planet known to have a significantly misaligned orbit.

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.

Advances in asteroseismology of solar-like stars, now provide a unique method to estimate the stellar inclination $i_{\star}$. This enables to evaluate the spin-orbit angle of transiting planetary systems, in a complementary fashion to the Rossiter-McLaughlin effect, a well-established method to estimate the projected spin-orbit angle $\lambda$. Although the asteroseismic method has been broadly applied to the Kepler data, its reliability has yet to be assessed intensively. In this work, we evaluate the accuracy of $i_{\star}$ from asteroseismology of solar-like stars using 3000 simulated power spectra. We find that the low signal-to-noise ratio of the power spectra induces a systematic under-estimate (over-estimate) bias for stars with high (low) inclinations. We derive analytical criteria for the reliable asteroseismic estimate, which indicates that reliable measurements are possible in the range of $20^\circ \lesssim i_{\star} \lesssim 80^\circ$ only for stars with high signal-to-noise ratio. We also analyse and measure the stellar inclination of 94 Kepler main-sequence solar-like stars, among which 33 are planetary hosts. According to our reliability criteria, a third of them (9 with planets, 22 without) have accurate stellar inclination. Comparison of our asteroseismic estimate of $v\sin{i_{\star}}$ against spectroscopic measurements indicates that the latter suffers from a large uncertainty possibly due to the modeling of macro-turbulence, especially for stars with projected rotation speed $v\sin{i_{\star}} \lesssim 5$ km/s. This reinforces earlier claims, and the stellar inclination estimated from the combination of measurements from spectroscopy and photometric variation for slowly rotating stars needs to be interpreted with caution.

We perform a systematic search for rings around 168 Kepler planet candidates with sufficient signal-to-noise ratios that are selected from all the short-cadence data. We fit ringed and ringless models to their lightcurves, and compare the fitting results to search for the signatures of planetary rings. First, we identify 29 tentative systems, for which the ringed models exhibit statistically significant improvement over the ringless models. The lightcurves of those systems are individually examined, but we are not able to identify any candidate that indicates evidence for rings. In turn, we find out several mechanisms of false-positives that would produce ring-like signals, and the null detection enables us to place upper limits on the size of rings. Furthermore, assuming the tidal alignment between axes of the planetary rings and orbits, we conclude that the occurrence rate of rings larger than twice the planetary radius is less than 15 percent. Even though the majority of our targets are short-period planets, our null detection provides statistical and quantitative constraints on largely uncertain theoretical models of origin, formation, and evolution of planetary rings.

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.

The Hyper Suprime-Cam Subaru Strategic Program (HSC-SSP) is a three-layered imaging survey aimed at addressing some of the most outstanding questions in astronomy today, including the nature of dark matter and dark energy. The survey has been awarded 300 nights of observing time at the Subaru Telescope and it started in March 2014. This paper presents the first public data release of HSC-SSP. This release includes data taken in the first 1.7 years of observations (61.5 nights) and each of the Wide, Deep, and UltraDeep layers covers about 108, 26, and 4 square degrees down to depths of i~26.4, ~26.5, and ~27.0 mag, respectively (5sigma for point sources). All the layers are observed in five broad bands (grizy), and the Deep and UltraDeep layers are observed in narrow bands as well. We achieve an impressive image quality of 0.6 arcsec in the i-band in the Wide layer. We show that we achieve 1-2 per cent PSF photometry (rms) both internally and externally (against Pan-STARRS1), and ~10 mas and 40 mas internal and external astrometric accuracy, respectively. Both the calibrated images and catalogs are made available to the community through dedicated user interfaces and database servers. In addition to the pipeline products, we also provide value-added products such as photometric redshifts and a collection of public spectroscopic redshifts. Detailed descriptions of all the data can be found online. The data release website is https://hsc-release.mtk.nao.ac.jp/.

Detection of a planetary ring of exoplanets remains as one of the most attractive but challenging goals in the field. We present a methodology of a systematic search for exoplanetary rings via transit photometry of long-period planets. The methodology relies on a precise integration scheme we develop to compute a transit light curve of a ringed planet. We apply the methodology to 89 long-period planet candidates from the Kepler data so as to estimate, and/or set upper limits on, the parameters of possible rings. While a majority of our samples do not have a sufficiently good signal-to-noise ratio for meaningful constraints on ring parameters, we find that six systems with a higher signal-to-noise ratio are inconsistent with the presence of a ring larger than 1.5 times the planetary radius assuming a grazing orbit and a tilted ring. Furthermore, we identify five preliminary candidate systems whose light curves exhibit ring-like features. After removing four false positives due to the contamination from nearby stars, we identify KIC 10403228 as a reasonable candidate for a ringed planet. A systematic parameter fit of its light curve with a ringed planet model indicates two possible solutions corresponding to a Saturn-like planet with a tilted ring. There also remain other two possible scenarios accounting for the data; a circumstellar disk and a hierarchical triple. Due to large uncertain factors, we cannot choose one specific model among the three.

We investigate the formation of close-in planets in near-coplanar eccentric hierarchical triple systems via the secular interaction between an inner planet and an outer perturber (Coplanar High-eccentricity Migration, CHEM). We generalize the previous work on the analytical condition for successful CHEM for point masses interacting only through gravity by taking into account the finite mass effect of the inner planet. We find that efficient CHEM requires that the systems should have m_1<<m_0 and m_1<<m_2. In addition to the gravity for point masses, we examine the importance of the short-range forces, and provide an analytical estimate of the migration time scale. We perform a series of numerical simulations in CHEM for systems consisting of a sun-like central star, giant gas inner planet and planetary outer perturber, including the short-range forces and stellar and planetary dissipative tides. We find that most of such systems end up with a tidal disruption, a small fraction of the systems produce prograde hot Jupiters (HJs), but no retrograde one. In addition, we extend CHEM to super-Earth mass range, and show that the formation of close-in super-Earths in prograde orbits is also possible. Finally, we carry out CHEM simulation for the observed hierarchical triple and counter-orbiting HJ systems. We find that CHEM can explain a part of the former systems, but it is generally very difficult to reproduce counter-orbiting HJ systems.

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.

PFS (Prime Focus Spectrograph), a next generation facility instrument on the 8.2-meter Subaru Telescope, is a very wide-field, massively multiplexed, optical and near-infrared spectrograph. Exploiting the Subaru prime focus, 2394 reconfigurable fibers will be distributed over the 1.3 deg field of view. The spectrograph has been designed with 3 arms of blue, red, and near-infrared cameras to simultaneously observe spectra from 380nm to 1260nm in one exposure at a resolution of ~1.6-2.7A. An international collaboration is developing this instrument under the initiative of Kavli IPMU. The project is now going into the construction phase aiming at undertaking system integration in 2017-2018 and subsequently carrying out engineering operations in 2018-2019. This article gives an overview of the instrument, current project status and future paths forward.

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.

Among a hundred transiting planets with a measured projected spin-orbit angle {\lambda}, several systems are suggested to be counter-orbiting. While they may be due to the projection effect, the mechanism to produce a counter-orbiting planet is not established. A promising scenario for the counter-orbiting planets is the extreme eccentricity evolution in near-coplanar hierarchical triple systems with eccentric inner and outer or- bits. We examine this scenario in detail by performing a series of systematic numerical simulations, and consider the possibility of forming hot Jupiters, especially counter- orbiting one under this mechanism with a distant sub-stellar perturber. We incorporate quadrupole and octupole secular gravitational interaction between the two orbits, and also short-range forces (correction for general relativity, star and inner planetary tide and rotational distortion) simultaneously. We find that most of systems are tidally disrupted and that a small fraction of survived planets turns out to be prograde. The formation of counter-orbiting hot Jupiters in this scenario is possible only in a very restricted parameter region, and thus very unlikely in practice.

Analysis of transit times in exoplanetary systems accurately provides an instantaneous orbital period, $P(t)$, of their member planets. A long-term monitoring of those transiting planetary systems puts limits on the variability of $P(t)$, which are translated into the constraints on the time variation of the gravitational constant $G$. We apply this analysis to $10$ transiting systems observed by the Kepler spacecraft, and find that $\Delta G/G\lesssim 5\times10^{-6}$ for 2009-2013, or $\dot{G}/G \lesssim 10^{-6}\,\mathrm{yr}^{-1}$ if $\dot{G}$ is constant. While the derived limit is weaker than those from other analyses, it is complementary to them and can be improved by analyzing numerous transiting systems that are continuously monitored.

We perform image stacking analysis of Sloan Digital Sky Survey (SDSS) photometric galaxies over the AKARI Far-Infrared Surveyor (FIS) maps at 65{\mu}m, 90{\mu}m, and 140{\mu}m. The resulting image profiles are decomposed into the central galaxy component (single term) and the nearby galaxy component (clustering term), as a function of the r-band magnitude, m_r of the central galaxy. We find that the mean far-infrared (FIR) flux of a galaxy with magnitude m_r is well fitted with f^s_{90{\mu}m}=13*10^{0.306(18-m_r)}[mJy]. The FIR amplitude of the clustering term is consistent with that expected from the angular-correlation function of the SDSS galaxies, but galaxy morphology dependence needs to be taken into account for a more quantitative conclusion. We also fit the spectral energy distribution of stacked galaxies at 65{\mu}m, 90{\mu}m, and 140{\mu}m, and derive a mean dust temperature of ~30K. This is consistent with the typical dust temperature of galaxies that are FIR luminous and individually detected.

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.

We reanalyse the time-variable lightcurves of the transiting planetary system PTFO 8-8695, in which a planet of 3 to 4 Jupiter mass orbits around a rapidly rotating pre-main-sequence star. Both the planetary orbital period of 0.448 days and the stellar spin period less than 0.671 days are unusually short, which makes PTFO 8-8695 an ideal system to check the model of gravity darkening and nodal precession. While the previous analysis of PTFO 8-8695 assumed that the stellar spin and planetary orbital periods are the same, we extend the analysis by discarding the spin-orbit synchronous condition, and find three different classes of solutions roughly corresponding to the nodal precession periods of 199$\pm$16, 475$\pm$21, and 827$\pm$53 days that reproduce the transit lightcurves observed in 2009 and 2010. We compare the predicted lightcurves of the three solutions against the photometry data of a few percent accuracy obtained at Koyama Astronomical Observatory in 2014 and 2015, and find that the solution with the precession period of 199$\pm$16 days is preferred even though preliminary. Future prospect and implications to other transiting systems are briefly discussed.

We propose a novel method to constrain the spatial extent of dust around galaxies through the measurement of dust temperature. Our method combines the dust emission of galaxies from far-infrared (FIR) image stacking analysis and the quasar reddening due to the dust absorption around galaxies. As a specific application of our method, we use the stacked FIR emission profiles of SDSS photometric galaxies over the IRAS 100$\mu$m map, and the recent measurement of the SDSS galaxy-quasar cross-correlation. If we adopt a single-temperature dust model, the resulting temperature is around 18K, which is consistent with a typical dust temperature for a central part of galaxies. If we assume an additional dust component with much lower temperature, the current data imply the temperature of the galactic dust needs to be higher, 20K to 30K. Since the model of the density and temperature distribution of dust adopted in the current paper is very simple, we cannot draw any strong conclusion at this point. Nevertheless our novel method with the elaborated theoretical model and multi-band measurement of dust will offer an interesting constraint on the statistical nature of galactic dust.

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.

The most widely used Galactic extinction map (Schlegel, Finkbeiner, & Davis 1998, SFD) is constructed assuming that the observed FIR fluxes entirely come from the Galactic dust. According to the earlier suggestion by Yahata et al. (2007), we consider how far-infrared (FIR) emission of galaxies affects the SFD map. We first compute the surface number density of SDSS DR7 galaxies as a function of the $r$-band extinction, $A_{r,\rm SFD}$. We confirm that the surface densities of those galaxies positively correlate with $A_{r,\rm SFD}$ for $A_{r,\rm SFD}<0.1$, as first discovered by Yahata et al. (2007) for SDSS DR4 galaxies. Next we construct an analytic model to compute the surface density of galaxies taking account of the contamination of their FIR emission. We adopt a log-normal probability distribution for the ratio of $100\mu {\rm m}$ and $r$-band luminosities of each galaxy, $y \equiv (\nu L)_{100\mu {\rm m}}/(\nu L)_r$. Then we search for the mean and r.m.s values of $y$ that fit the observed anomaly using the analytic model. The required values to reproduce the anomaly are roughly consistent with those measured from the stacking analysis of SDSS galaxies (Kashiwagi, Yahata, & Suto 2013). Due to the limitation of our statistical modeling, we are not yet able to remove the FIR contamination of galaxies from the extinction map. Nevertheless the agreement with the model prediction suggests that the FIR emission of galaxies is mainly responsible for the observed anomaly. While the corresponding systematic error in the Galactic extinction map is 0.1 to 1mmag, it is directly correlated with galaxy clustering, and thus needs to be carefully examined in precision cosmology.

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.

The angle between the stellar spin and the planetary orbit axes (spin-orbit angle) is supposed to carry valuable information on the initial condition of the planet formation and the subsequent migration history. Indeed current observations of the Rossiter- McLaughlin effect have revealed a wide range of spin-orbit misalignments for transiting exoplanets. We examine in detail the tidal evolution of a simple system comprising a Sun-like star and a hot Jupiter adopting the equilibrium tide and the inertial wave dissipation effects simultaneously. We find that the combined tidal model works as a very efficient realignment mechanism; it predicts three distinct states of the spin-orbit angle (i.e., parallel, polar, and anti-parallel orbits) for a while, but the latter two states eventually approach the parallel spin-orbit configuration. The intermediate spin-orbit angles as measured in recent observations are difficult to be achieved. Therefore the current model cannot reproduce the observed broad distribution of the spin-orbit angles, at least in its simple form. This indicates that the observed diversity of the spin-orbit angles may emerge from more complicated interactions with outer planets and/or may be the consequence of the primordial misalignment between the proto-planetary disk and the stellar spin, which requires future detailed studies.

The KOI-94 system is a closely-packed, multi-transiting planetary system discovered by the Kepler space telescope. It is known as the first system that exhibited a rare event called a "planet-planet eclipse (PPE)," in which two planets partially overlap with each other in their double-transit phase. In this paper, we constrain the parameters of the KOI-94 system with an analysis of the transit timing variations (TTVs). Such constraints are independent of the radial velocity (RV) analysis recently performed by Weiss and coworkers, and valuable in examining the reliability of the parameter estimate using TTVs. We numerically fit the observed TTVs of KOI-94c, KOI-94d, and KOI-94e for their masses, eccentricities, and longitudes of periastrons, and obtain the best-fit parameters including $m_{\rm c} = 9.4_{-2.1}^{+2.4} M_{\oplus}$, $m_{\rm d} = 52.1_{-7.1}^{+6.9} M_{\oplus}$, $m_{\rm e} = 13.0_{-2.1}^{+2.5} M_{\oplus}$, and $e \lesssim 0.1$ for all the three planets. While these values are mostly in agreement with the RV result, the mass of KOI-94d estimated from the TTV is significantly smaller than the RV value $m_{\rm d} = 106 \pm 11 M_{\oplus}$. In addition, we find that the TTV of the outermost planet KOI-94e is not well reproduced in the current modeling. We also present analytic modeling of the PPE and derive a simple formula to reconstruct the mutual inclination of the two planets from the observed height, central time, and duration of the brightening caused by the PPE. Based on this model, the implication of the results of TTV analysis for the time of the next PPE is discussed.

The Subaru Prime Focus Spectrograph (PFS) is a massively-multiplexed fiber-fed optical and near-infrared 3-arm spectrograph (N_fiber=2400, 380<lambda<1260nm, 1.3 degree diameter FoV), offering unique opportunities in survey astronomy. Here we summarize the science case feasible for a survey of Subaru 300 nights. We describe plans to constrain the nature of dark energy via a survey of emission line galaxies spanning a comoving volume of 9.3 (Gpc/h)^3 in the redshift range 0.8<z<2.4. In each of 6 redshift bins, the cosmological distances will be measured to 3% precision via BAO, and redshift-space distortions will be used to constrain structure growth to 6% precision. In the GA program, radial velocities and chemical abundances of stars in the Milky Way and M31 will be used to infer the past assembly histories of spiral galaxies and the structure of their dark matter halos. Data will be secured for 10^6 stars in the Galactic thick-disk, halo and tidal streams as faint as V~22, including stars with V < 20 to complement the goals of the Gaia mission. A medium-resolution mode with R = 5000 to be implemented in the red arm will allow the measurement of multiple alpha-element abundances and more precise velocities for Galactic stars, elucidating the detailed chemo-dynamical structure and evolution of each of the main stellar components of the Milky Way Galaxy and of its dwarf spheroidal galaxies. For the extragalactic program, our simulations suggest the wide avelength range will be powerful in probing the galaxy population and its clustering over a wide redshift range. We propose to conduct a color-selected survey of 1<z<2 galaxies and AGN over 16 deg^2 to J~23.4, yielding a fair sample of galaxies with stellar masses above ~10^{10}Ms at z~2. A two-tiered survey of higher redshift LBGs and LAEs will quantify the properties of early systems close to the reionization epoch.

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 variability of major atmospheric absorption features in the disk-integrated spectra of the Earth with future application to Earth-analogs in mind, concentrating on the diurnal timescale. We first analyze observations of the Earth provided by the EPOXI mission, and find 5-20% fractional variation of the absorption depths of H2O and O2 bands, two molecules that have major signatures in the observed range. From a correlation analysis with the cloud map data from the Earth Observing Satellite (EOS), we find that their variation pattern is primarily due to the uneven cloud cover distribution. In order to account for the observed variation quantitatively, we consider a simple opaque cloud model, which assumes that the clouds totally block the spectral influence of the atmosphere below the cloud layer, equivalent to assuming that the incident light is completely scattered at the cloud top level. The model is reasonably successful, and reproduces the EPOXI data from the pixel-level EOS cloud/water vapor data. A difference in the diurnal variability patterns of H2O and O2 bands is ascribed to the differing vertical and horizontal distribution of those molecular species in the atmosphere. On the Earth, the inhomogeneous distribution of atmospheric water vapor is due to the existence of its exchange with liquid and solid phases of H2O on the planet's surface on a timescale short compared to atmospheric mixing times. If such differences in variability patterns were detected in spectra of Earth-analogs, it would provide the information on the inhomogeneous composition of their atmospheres.

We have performed stacking image analyses of galaxies over the Galactic extinction map constructed by Schlegel, Finkbeiner & Davis (1998). We select ~10^7 galaxies in total from the Sloan Digital Sky Survey (SDSS) DR7 photometric catalog. We detect clear signatures of the enhancement of the extinction in r-band, $\Delta A_r$, around galaxies, indicating that the extinction map is contaminated by their FIR (far infrared) emission. The average amplitude of the contamination per galaxy is well fitted to $\Delta A_r(m_r) = 0.64 \times 10^{0.17(18-m_r)}$ [mmag]. While this value is very small, it is directly associated with galaxies and may have a systematic effect on galaxy statistics. Indeed this correlated contamination leads to a relatively large anomaly of galaxy surface number densities against the SFD extinction A_SFD discovered by Yahata et al. (2007). We model the radial profiles of stacked galaxy images, and find that the FIR signal around each galaxy does not originate from the central galaxy alone, but is dominated by the contributions of nearby galaxies via galaxy angular clustering. The separation of the single galaxy and the clustering terms enables us to infer the statistical relation of the FIR and r-band fluxes of galaxies and also to probe the flux-weighted cross-correlation of galaxies, down to the magnitudes that are difficult to probe directly for individual objects. We repeat the same stacking analysis for SDSS DR6 photometric quasars and discovered the similar signatures but with weaker amplitudes. The implications of the present results for galaxy and quasar statistics and for correction to the Galactic extinction map are briefly discussed.

We report a joint analysis of the Rossiter-McLaughlin (RM) effect with Subaru and the Kepler photometry for Kepler Object of Interest (KOI) 94 system. The system comprises four transiting planet candidates with orbital periods of 22.3 (KOI-94.01), 10.4 (KOI-94.02), 54.3 (KOI-94.03), and 3.7 (KOI-94.04) days from the Kepler photometry. We performed the radial velocity (RV) measurement of the system with the Subaru 8.2 m telescope on August 10, 2012 (UT), covering a complete transit of KOI-94.01 for $\sim 6.7$ hours. The resulting RV variation due to the RM effect spectroscopically confirms that KOI-94.01 is indeed the transiting planet, and implies that its orbital axis is well aligned with the stellar spin axis; the projected spin-orbit angle $\lambda$ is estimated as $-6_{-11}^{+13}$ deg. This is the first measurement of the RM effect for a multiple transiting system. Remarkably, the archived Kepler lightcurve around BJD=2455211.5 (date in UT January 14/15, 2010) indicates a "ouble transit" event of KOI-94.01 and KOI-94.03, in which the two planets transit the stellar disk simultaneously. Moreover, the two planets partially overlap each other, and exhibit a "planet-planet eclipse" around the transit center. This provides a rare opportunity to put tight constraints on the configuration of the two transiting planets by joint analysis with our Subaru RM measurement. Indeed, we find that the projected mutual inclination of KOI-94.01 and KOI-94.03 is estimated to be $\delta = -1.15 \pm 0.55$ deg. Implications for the migration model of multiple planet systems are also discussed.

We present an investigation of spin-orbit angles for planetary system candidates reported by Kepler. By combining the rotational period $P_s$ inferred from the flux variation due to starspots and the projected rotational velocity $V\sin I_s$ and stellar radius obtained by a high resolution spectroscopy, we attempt to estimate the inclination $I_s$ of the stellar spin axis with respect to the line-of-sight. For transiting planetary systems, in which planetary orbits are edge-on seen from us, the stellar inclination $I_s$ can be a useful indicator of a spin-orbit alignment/misalignment. We newly conducted spectroscopic observations with Subaru/HDS for 15 KOI systems, whose lightcurves show periodic flux variations. After detailed analyses of their lightcurves and spectra, it turned out that some of them are binaries, or the flux variations are too coherent to be caused by starspots, probably representing ellipsoidal variations, and consequently we could constrain stellar inclinations $I_s$ for eight systems. Among them, KOI-262 and 280 are in good agreement with $I_s=90^\circ$ suggesting a spin-orbit alignment, while at least one system, KOI-261, shows a possible spin-orbit misalignment. We also obtain a small $I_s$ for KOI-1463, but the transiting companion seems to be a star rather than a planet. The results for KOI-257, 269, 367, and 974 are ambiguous, and can be explained with either misalignments or moderate differential rotation. Since our method can be applied to any system having starspots regardless of the planet size, future observations will allow for the expansion of the parameter space in which the spin-orbit relations are investigated.

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

We report on observations of the Rossiter-McLaughlin (RM) effect for the XO-3 exoplanetary system. The RM effect for the system was previously measured by two different groups, but their results were statistically inconsistent. To obtain a decisive result we observed two full transits of XO-3b with the Subaru 8.2-m telescope. By modeling these data with a new and more accurate analytic formula for the RM effect, we find the projected spin-orbit angle to be \lambda=37.3 deg \pm 3.0 deg, in good agreement with the previous finding by Winn et al. (2009). In addition, an offset of ~22 m/s was observed between the two transit datasets. This offset could be a signal of a third body in the XO-3 system, a possibility that should be checked with future observations. We also attempt to search for a possible signature of the stellar differential rotation in the RM data for the first time, and put weak upper limits on the differential rotation parameters.

We present an improved formula for the anomalous radial velocity of the star during planetary transits due to the Rossiter-McLaughlin (RM) effect. The improvement comes from a more realistic description of the stellar absorption line profiles, taking into account stellar rotation, macroturbulence, thermal broadening, pressure broadening, and instrumental broadening. Although the formula is derived for the case in which radial velocities are measured by cross-correlation, we show through numerical simulations that the formula accurately describes the cases where the radial velocities are measured with the iodine absorption-cell technique. The formula relies on prior knowledge of the parameters describing macroturbulence, instrumental broadening and other broadening mechanisms, but even 30% errors in those parameters do not significantly change the results in typical circumstances. We show that the new analytic formula agrees with previous ones that had been computed on a case-by-case basis via numerical simulations. Finally, as one application of the new formula, we reassess the impact of the differential rotation on the RM velocity anomaly. We show that differential rotation of a rapidly rotating star may have a significant impact on future RM observations.

As a test-bed for future investigations of directly imaged terrestrial exoplanets, we present the recovery of the surface components of the Earth from multi-band diurnal light curves obtained with the EPOXI spacecraft. We find that the presence and longitudinal distribution of ocean, soil and vegetation are reasonably well reproduced by fitting the observed color variations with a simplified model composed of a priori known albedo spectra of ocean, soil, vegetation, snow and clouds. The effect of atmosphere, including clouds, on light scattered from surface components is modeled using a radiative transfer code. The required noise levels for future observations of exoplanets are also determined. Our model-dependent approach allows us to infer the presence of major elements of the planet (in the case of the Earth, clouds and ocean) with observations having S/N $\gtrsim 10$ in most cases and with high confidence if S/N $\gtrsim 20$. In addition, S/N $\gtrsim 100$ enables us to detect the presence of components other than ocean and clouds in a fairly model-independent way. Degradation of our inversion procedure produced by cloud cover is also quantified. While cloud cover significantly dilutes the magnitude of color variations compared to the cloudless case, the pattern of color changes remains. Therefore, the possibility of investigating surface features through light curve fitting remains even for exoplanets with cloud cover similar to the Earth's.

The measurement of the light scattered from extrasolar planets informs atmospheric and formation models. With the discovery of many hot Jupiter planets orbiting nearby stars, this motivates the development of robust methods of characterisation from follow up observations. In this paper we discuss two methods for determining the planetary albedo in transiting systems. First, the most widely used method for measuring the light scattered by hot Jupiters (Collier Cameron et al.) is investigated for application for typical echelle spectra of a transiting planet system, showing that detection requires high signal-to-noise ratio data of bright planets. Secondly a new Fourier analysis method is also presented, which is model-independent and utilises the benefits of the reduced number of unknown parameters in transiting systems. This approach involves solving for the planet and stellar spectra in Fourier space by least-squares. The sensitivities of the methods are determined via Monte Carlo simulations for a range of planet-to-star fluxes. We find the Fourier analysis method to be better suited to the ideal case of typical observations of a well constrained transiting system than the Collier Cameron et al. method. We apply the Fourier analysis method for extracting the light scattered by transiting hot Jupiters from high resolution spectra to echelle spectra of HD 209458 and HD 189733. Unfortunately we are unable to improve on the previous upper limit of the planet-to-star flux for HD 209458b set by space-based observations. A 1{\sigma}upper limit on the planet-to-star flux of HD 189733b is measured in the wavelength range of 558.83-599.56 nm yielding {\epsilon} < 4.5 \times 10-4. Improvement in the measurement of the upper limit of the planet-to-star flux of this system, with ground-based capabilities, requires data with a higher signal-to-noise ratio, and increased stability of the telescope.

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.

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.

Answering well-known fundamental questions is usually regarded as the major goal of science. Discovery of other unknown and fundamental questions is, however, even more important. Recognition that "we didn't know anything" is the basic scientific driver for the next generation. Cosmology indeed enjoys such an exciting epoch. What is the composition of our universe? This is one of the well-known fundamental questions that philosophers, astronomers and physicists have tried to answer for centuries. Around the end of the last century, cosmologists finally recognized that "We didn't know anything". Except for atoms that comprise slightly less than 5% of the universe, our universe is apparently dominated by unknown components; 23% is the known unknown (dark matter), and 72% is the unknown unknown (dark energy). In the course of answering a known fundamental question, we have discovered an unknown, even more fundamental, question: "What is dark matter? What is dark energy?" There are a variety of realistic particle physics models for dark matter, and its experimental detection may be within reach. On the other hand, it is fair to say that there is no widely accepted theoretical framework to describe the nature of dark energy. This is exactly why astronomical observations will play a key role in unveiling its nature. I will review our current understanding of the "dark sky", and then present on-going Japanese project, SuMIRe, to discover even more unexpected questions.

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 obtain analytical expressions for the velocity anomaly due to the Rossiter- McLaughlin effect, for the case when the anomalous radial velocity is obtained by cross-correlation with a stellar template spectrum. In the limit of vanishing width of the stellar absorption lines, our result reduces to the formula derived by Ohta et al. (2005), which is based on the first moment of distorted stellar lines. Our new formula contains a term dependent on the stellar linewidth, which becomes important when rotational line broadening is appreciable. We generate mock transit spectra for four existing exoplanetary systems (HD17156, TrES-2, TrES- 4, and HD209458) following the procedure of Winn et al. (2005), and find that the new formula is in better agreement with the velocity anomaly extracted from the mock data. Thus, our result provides a more reliable analytical description of the velocity anomaly due to the Rossiter-McLaughlin effect, and explains the previously observed dependence of the velocity anomaly on the stellar rotation velocity.

We present an improved measurement of the Rossiter-McLaughlin effect for the exoplanetary system HD 17156, based on radial-velocity data gathered with the Subaru 8.2m telescope throughout the planetary transit of UT 2008 November 7. The data allow for a precise and independent determination of the projected spin-orbit angle of this system: $\lambda = 10.0^{\circ} \pm 5.1^{\circ}$. This result supersedes the previous claim of $\lambda = 62^{\circ} \pm 25^{\circ}$ by Narita et al., which was based on lower-precision data with poor statistics. Thus the stellar spin and planetary orbital axes of the HD 17156 system are likely to be well-aligned, despite the planet's large orbital eccentricity suggesting a history of strong dynamical interactions.

We present photometric and spectroscopic observations of the 2009 Feb. 2 transit of the exoplanet XO-3b. The new data show that the planetary orbital axis and stellar rotation axis are misaligned, as reported earlier by Hebrard and coworkers. We find the angle between the sky projections of the two axes to be 37.3 +/- 3.7 deg, as compared to the previously reported value of 70 +/- 15 deg. The significance of this discrepancy is unclear because there are indications of systematic effects. XO-3b is the first exoplanet known to have a highly inclined orbit relative to the equatorial plane of its parent star, and as such it may fulfill the predictions of some scenarios for the migration of massive planets into close-in orbits. We revisit the statistical analysis of spin-orbit alignment in hot-Jupiter systems. Assuming the stellar obliquities to be drawn from a single Rayleigh distribution, we find the mode of the distribution to be 13^{+5}_{-2} deg. However, it remains the case that a model representing two different migration channels--in which some planets are drawn from a perfectly-aligned distribution and the rest are drawn from an isotropic distribution--is favored over a single Rayleigh distribution.

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

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

We present theoretical predictions for photometric and spectroscopic signatures of rings around transiting extrasolar planets. On the basis of a general formulation for the transiting signature in the stellar light curve and the velocity anomaly due to the Rossiter effect, we compute the expected signals analytically for a face-on ring system, and numerically for more general configurations. We study the detectability of a ring around a transiting planet located at a=3AU for a variety of obliquity and azimuthal angles, and find that it is possible to detect the ring signature both photometrically and spectroscopically unless the ring is almost edge-on (i.e., the obliquity angle of the ring \theta is much less than unity). We also consider the detectability of planetary rings around a close-in planet, HD 209458b ($\theta\approx 90^\circ-i_{\rm orbit} \approx 3^\circ.32$), and Saturn ($\theta \approx 26^\circ.7$) as illustrative examples. While the former is difficult to detect with the current precision (photometric precision of 10^-4 and radial velocity precision of 1 m/s), a marginal detection of the latter is possible photometrically. If the future precision of the radial velocity measurement reaches even below 0.1 m/s, they will be even detectable from the ground-based spectroscopic observations.

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 critically examine how well the evolution of large-scale density perturbations is followed in cosmological $N$-body simulations. We first run a large volume simulation and perform a mode-by-mode analysis in three-dimensional Fourier space. We show that the growth of large-scale fluctuations significantly deviates from linear theory predictions. The deviations are caused by {\it nonlinear} coupling with a small number of modes at largest scales owing to finiteness of the simulation volume. We then develop an analytic model based on second-order perturbation theory to quantify the effect. Our model accurately reproduces the simulation results. For a single realization, the second-order effect appears typically as ``zig-zag'' patterns around the linear-theory prediction, which imprints artificial ``oscillations'' that lie on the real baryon-acoustic oscillations. Although an ensemble average of a number of realizations approaches the linear theory prediction, the dispersions of the realizations remain large even for a large simulation volume of several hundred megaparsecs on a side. For the standard $\Lambda$CDM model, the deviations from linear growth rate are as large as 10 percent for a simulation volume with $L = 500h^{-1}$Mpc and for a bin width in wavenumber of $\Delta k = 0.005h$Mpc$^{-1}$, which are comparable to the intrinsic variance of Gaussian random realizations. We find that the dispersions scales as $\propto L^{-3/2} \Delta k^{-1/2}$ and that the mean dispersion amplitude can be made smaller than a percent only if we use a very large volume of $L > 2h^{-1}$Gpc. The finite box size effect needs to be appropriately taken into account when interpreting results from large-scale structure simulations for future dark energy surveys using baryon acoustic oscillations.

We present new spectroscopic and photometric observations of the HAT-P-1 planetary system. Spectra obtained during three transits exhibit the Rossiter-McLaughlin effect, allowing us to measure the angle between the sky projections of the stellar spin axis and orbit normal, \lambda = 3.7 +/- 2.1 degrees. The small value of \lambda for this and other systems suggests that the dominant planet migration mechanism preserves spin-orbit alignment. Using two new transit light curves, we refine the transit ephemeris and reduce the uncertainty in the orbital period by an order of magnitude. We find a upper limit on the orbital eccentricity of 0.067, with 99% confidence, by combining our new radial-velocity measurements with those obtained previously.

We monitored the Doppler shift of the G0V star TrES-2 throughout a transit of its giant planet. The anomalous Doppler shift due to stellar rotation (the Rossiter-McLaughlin effect) is discernible in the data, with a signal-to-noise ratio of 2.9, even though the star is a slow rotator. By modeling this effect we find that the planet's trajectory across the face of the star is tilted by -9 +/- 12 degrees relative to the projected stellar equator. With 98% confidence, the orbit is prograde.

We present a new technique for detecting scattered starlight from transiting, close-orbiting extrasolar giant planets (CEGPs) that has the virtues of simplicity, robustness, linearity, and model-independence. Given a series of stellar spectra obtained over various phases of the planetary orbit, the goal is to measure the strength of the component scattered by the planet relative to the component coming directly from the star. We use two complementary strategies, both of which rely on the predictable Doppler shifts of both components and on combining the results from many spectral lines and many exposures. In the first strategy, we identify segments of the stellar spectrum that are free of direct absorption lines and add them after Doppler-shifting into the planetary frame. In the second strategy, we compare the distribution of equivalent-width ratios of the scattered and direct components. Both strategies are calibrated with a ``null test'' in which scrambled Doppler shifts are applied to the spectral segments. As an illustrative test case, we apply our technique to spectra of HD 209458 taken when the planet was near opposition (with orbital phases ranging from 11 to 34$\arcdeg$, where 0$\arcdeg$ is at opposition), finding that the planet-to-star flux ratio is $(1.4 \pm 2.9)\times10^{-4}$ in the wavelength range 554$-$681 nm. This corresponds to a geometric albedo of $0.8 \pm 1.6$, assuming the phase function of a Lambert sphere. Although the result is not statistically significant, the achieved sensitivity and relatively small volume of data upon which it is based are very encouraging for future ground-based spectroscopic studies of scattered light from transiting CEGP systems.

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 report a measurement of the Rossiter--McLaughlin effect in the transiting extrasolar planetary system TrES-1, via simultaneous spectroscopic and photometric observations with the Subaru and MAGNUM telescopes. By modeling the radial velocity anomaly that was observed during a transit, we determine the sky-projected angle between the stellar spin axis and the planetary orbital axis to be $\lambda = 30 \pm 21$ [deg]. This is the third case for which $\lambda$ has been measured in a transiting exoplanetary system, and the first demonstration that such measurements are possible for relatively faint host stars ($V \sim 12$, as compared to $V \sim 8$ for the other systems). We also derive a time of mid-transit, constraints on the eccentricity of the TrES-1b orbit ($e = 0.048 \pm 0.025$), and upper limits on the mass of the Trojan companions ($\lesssim$14 $M_{\oplus}$) at the 3$\sigma$ level.

An acoustic oscillation of the primeval photon-baryon fluid around the decoupling time imprints a characteristic scale in the galaxy distribution today, known as the baryon acoustic oscillation (BAO) scale. Several on-going and/or future galaxy surveys aim at detecting and precisely determining the BAO scale so as to trace the expansion history of the universe. We consider nonlinear and redshift-space distortion effects on the shifts of the BAO scale in $k$-space using perturbation theory. The resulting shifts are indeed sensitive to different choices of the definition of the BAO scale, which needs to be kept in mind in the data analysis. We present a toy model to explain the physical behavior of the shifts. We find that the BAO scale defined as in Percival et al. (2007) indeed shows very small shifts ($\lesssim$ 1%) relative to the prediction in {\it linear theory} in real space. The shifts can be predicted accurately for scales where the perturbation theory is reliable.

We consider the extent to which future imaging surveys of galaxies can distinguish between dark energy and modified gravity models for the origin of the cosmic acceleration. Dynamical dark energy models may have similar expansion rates as models of modified gravity, yet predict different growth of structure histories. We parameterize the cosmic expansion by the two parameters, $w_0$ and $w_a$, and the linear growth rate of density fluctuations by Linder's $\gamma$, independently. Dark energy models generically predict $\gamma \approx 0.55$, while the DGP model $\gamma \approx 0.68$. To determine if future imaging surveys can constrain $\gamma$ within 20 percent (or $\Delta\gamma<0.1$), we perform the Fisher matrix analysis for a weak lensing survey such as the on-going Hyper Suprime-Cam (HSC) project. Under the condition that the total observation time is fixed, we compute the Figure of Merit (FoM) as a function of the exposure time $\texp$. We find that the tomography technique effectively improves the FoM, which has a broad peak around $\texp\simeq {\rm several}\sim 10$ minutes; a shallow and wide survey is preferred to constrain the $\gamma$ parameter. While $\Delta\gamma < 0.1$ cannot be achieved by the HSC weak-lensing survey alone, one can improve the constraints by combining with a follow-up spectroscopic survey like WFMOS and/or future CMB observations.

The short-period exoplanet HD 147506b (also known as HAT-P-2b) has an eccentric orbit, raising the possibility that it migrated through planet-planet scattering or Kozai oscillations accompanied by tidal dissipation. Either of these scenarios could have significantly tilted the orbit relative to the host star's equatorial plane. Here we present spectroscopy of a transit of HD 147506b, and assess the spin-orbit alignment via the Rossiter-McLaughlin effect. We find the sky projections of the stellar spin axis and orbital axis to be aligned within 14 deg. Thus we find no corroborating evidence for scattering or Kozai migration, although these scenarios cannot be ruled out with the present data.

We explore observational constraints on possible deviations from Newtonian gravity by means of large-scale clustering of galaxies. We measure the power spectrum and the bispectrum of Sloan Digital Sky Survey galaxies and compare the result with predictions in an empirical model of modified gravity. Our model assumes an additional Yukawa-like term with two parameters that characterize the amplitude and the length scale of the modified gravity. The model predictions are calculated using two methods; the second-order perturbation theory and direct N-body simulations. These methods allow us to study non-linear evolution of large-scale structure. Using the simulation results, we find that perturbation theory provides reliable estimates for the power spectrum and the bispectrum in the modified Newtonian model. We also construct mock galaxy catalogues from the simulations, and derive constraints on the amplitude and the length scale of deviations from Newtonian gravity. The resulting constraints from power spectrum are consistent with those obtained in our earlier work, indicating the validity of the previous empirical modeling of gravitational nonlinearity in the modified Newtonian model. If linear biasing is adopted, the bispectrum of the SDSS galaxies yields constraints very similar to those from the power spectrum. If we allow for the nonlinear biasing instead, we find that the ratio of the quadratic to linear biasing coefficients, b_2/b_1, should satisfy -0.4 < b_2/b_1<0.3 in the modified Newtonian model.

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 present photometry of HD 189733 during eight transits of its close-in giant planet, and out-of-transit photometry spanning two years. Using the transit photometry, we determine the stellar and planetary radii and the photometric ephemeris. Outside of transits, there are quasiperiodic flux variations with a 13.4 day period that we attribute to stellar rotation. In combination with previous results, we derive upper limits on the orbital eccentricity, and on the true angle between the stellar rotation axis and planetary orbit (as opposed to the angle between the projections of those axes on the sky).

We examine the age-dependence of dark matter halo clustering in an unprecedented accuracy using a set of 7 high-resolution cosmological simulations each with $N=1024^3$ particles. We measure the bias parameters for halos over a large mass range using the cross-power-spectrum method that can effectively suppress the random noise even in the sparse sampling of the most massive halos. This enables us to find, for the first time, that younger halos are more strongly clustered than older ones for halo masses $M>10M_{\ast}$, where $M_\ast$ is the characteristic nonlinear mass scale. For $M<M_{\ast}$, our results confirm the previous finding of Gao et al. that older halos are clustered more strongly than the younger ones. We also study the halo bias as a function of halo concentration, and find that the concentration dependence is weaker than the age dependence for $M<M_{\ast}$, but stronger for $M\ga 50 M_{\ast}$. The accurate and robust measurement of the age dependences of halo bias points to a limitation of the simple excursion set theory which predicts that the formation and structure of a halo of given mass is independent of its environment.

We consider nonlinear biasing models of galaxies with particular attention to a correlation between linear and quadratic biasing coefficients, b_1 and b_2. We first derive perturbative expressions for b_1 and b_2 in halo and peak biasing models. Then we compute power spectra and bispectra of dark matter particles and halos using N-body simulation data and of volume-limited subsamples of Sloan Digital Sky Survey (SDSS) galaxies, and determine their b_1 and b_2. We find that the values of those coefficients at linear regimes (k<0.2h/Mpc) are fairly insensitive to the redshift-space distortion and the survey volume shape. The resulting normalized amplitudes of bispectra, Q, for equilateral triangles, are insensitive to the values of b_1 implying that b_2 indeed correlates with b_1. The present results explain the previous finding of Kayo et al. (2004) for the hierarchical relation of three-point correlation functions of SDSS galaxies. While the relations between b_1 and b_2 are quantitatively different for specific biasing models, their approximately similar correlations indicate a fairly generic outcome of the biasing due to the gravity in primordial Gaussian density fields.

We present spectroscopy of a transit of the exoplanet HD 189733b. By modeling the Rossiter-McLaughlin effect (the anomalous Doppler shift due to the partial eclipse of the rotating stellar surface), we find the angle between the sky projections of the stellar spin axis and orbit normal to be lambda = -1.4 +/- 1.1 deg. This is the third case of a ``hot Jupiter'' for which lambda has been measured. In all three cases lambda is small, ruling out random orientations with 99.96% confidence, and suggesting that the inward migration of hot Jupiters generally preserves spin-orbit alignment.

We discuss how the baryon acoustic oscillation (BAO) signatures in the galaxy power spectrum can distinguish between modified gravity and the cosmological constant as the source of cosmic acceleration. To this end we consider a model characterized by a parameter $n$, which corresponds to the Dvali-Gabadadze-Porrati (DGP) model if $n=2$ and reduces to the standard spatially flat cosmological constant concordance model for $n$ equal to infinity. We find that the different expansion histories of the modified gravity models systematically shifts the peak positions of BAO. A preliminary analysis using the current SDSS LRG sample indicates that the original DGP model is disfavored unless the matter density parameter exceeds 0.3. The constraints will be strongly tightened with future spectroscopic samples of galaxies at high redshifts. We demonstrate that WFMOS, in collaboration with other surveys such as Planck, will powerfully constrain modified gravity alternatives to dark energy as the explanation of cosmic acceleration.

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

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

We compare the most successful and widely used map of Galactic dust extinction, provided by Schlegel, Finkbeiner & Davis (1998; hereafter SFD), to the galaxy number counts in the Sloan Digital Sky Survey (SDSS) photometric/spectroscopic DR4 sample. We divide the SDSS survey area into 69 disjoint subregions according to the dust extinction provided by SFD and compare the surface number density of galaxies in each subregion. As expected, the galaxy surface number density decreases with increasing extinction but only for SFD extinction values above about 0.1 to 0.2 magnitudes (depending on the band). At lower values of the SFD extinction, we find that the sky surface density of galaxies increases with increasing extinction, precisely the opposite of the effect expected from Galactic dust. We suggest that the far infrared (FIR) brightness of the sky in regions of true low dust extinction is significantly ``contaminated'' by the FIR emission from background galaxies. We show that such an explanation is both qualitatively and quantitatively consistent with the available data. Based on this interpretation we conclude that systematic errors in the SFD extinction map due to extragalactic FIR emission are quite small, of order hundredths of a magnitude, but nevertheless statistically detectable. (Abridged)

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 report the results of a prototype photometric search for transiting extrasolar planets using Subaru Suprime-Cam. Out of about 100,000 stars monitored around the Galactic plane (l=90 degree,b=0 degree), we find that 7,700 (27,000) stars show the photometric precision below 1% (3%) for 60 second exposures which is required to detect extrasolar planets by the transit method. Thus Suprime-Cam has the photometric stability and accuracy required for a transiting planet survey. During this observing run, we detected three transiting planetary candidates (i'-band magnitude around 18.5) which exhibit a single full transit-like light curve with a fractional depth of <5%. While future photometric and/or spectroscopic follow-ups remain to be done, the estimated parameters for the three systems are consistent with planetary companions around main-sequence stars. We also found two eclipsing binary candidates and eleven variable stars exhibiting W UMa-like light curves.

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

We present a first analysis of the clustering of SDSS galaxies using the distribution function of the sum of Fourier phases. This statistic was recently proposed by one of authors as a new method to probe phase correlations of cosmological density fields. Since the Fourier phases are statistically independent of Fourier amplitudes, the phase statistic plays a complementary role to the conventional two-point statistics of galaxy clustering. In particular, we focus on the distribution functions of phase sum over three closed wavevectors as a function of the configurations of triangle wavevectors. We find that the observed distribution functions of phase sum are in good agreement with the lowest-order approximation from perturbation theory. For direct comparison with observations, we construct mock catalogs from N-body simulations taking account of several observational effects such as the survey geometry, the redshift distortion, and the discreteness due to the limited number of data. Indeed the observed phase correlations for the galaxies in the range of absolute magnitude -22<Mr<-18 agree well with those for Lambda CDM predictions with sigma8=0.9 evolved from Gaussian initial condition. This agreement implies that the galaxy biasing is approximately linear in redshift space. Instead, assuming that the galaxy biasing is described by a quadratic deterministic function at scale larger than 30Mpc/h, we can constrain the ratio of the quadratic biasing parameter b2 to the linear biasing parameter b1. We find that resulting b2/b1 is mainly dependent on sigma8 and that it is well fitted by b2/b1=0.54(+-0.06)-0.62(+-0.08)sigma8 in CDM models. Indeed, b2/b1 is nearly zero when sigma8=0.9.

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 determine the stellar, planetary, and orbital properties of the transiting planetary system HD 209458, through a joint analysis of high-precision radial velocities, photometry, and timing of the secondary eclipse. Of primary interest is the strong detection of the Rossiter-McLaughlin effect, the alteration of photospheric line profiles that occurs because the planet occults part of the rotating surface of the star. We develop a new technique for modeling this effect, and use it to determine the inclination of the planetary orbit relative to the apparent stellar equator [lambda = (-4.4 +/- 1.4) degrees], and the line-of-sight rotation speed of the star [v*sin(I) = (4.70 +/- 0.16) km/s]. The uncertainty in these quantities has been reduced by an order of magnitude relative to the pioneering measurements by Queloz and collaborators. The small but nonzero misalignment is probably a relic of the planet formation epoch, because the expected timescale for tidal coplanarization is larger than the age of the star. Our determination of v*sin(I) is a rare case in which rotational line broadening has been isolated from other broadening mechanisms.

We have searched for absorption in several common atomic species due to the atmosphere or exosphere of the transiting extrasolar planet HD 209458b, using high precision optical spectra obtained with the Subaru High Dispersion Spectrograph (HDS). Previously we reported an upper limit on H alpha absorption of 0.1% (3 sigma) within a 5.1\AA band. Using the same procedure, we now report upper limits on absorption due to the optical transitions of Na D, Li, H alpha, H beta, H gamma, Fe, and Ca. The 3 sigma upper limit for each transition is approximately 1% within a 0.3\AA band (the core of the line), and a few tenths of a per cent within a 2\AA band (the full line width). The wide-band results are close to the expected limit due to photon-counting (Poisson) statistics, although in the narrow-band case we have encountered unexplained systematic errors at a few times the Poisson level. These results are consistent with all previously reported detections and upper limits, but are significantly more sensitive.

We have identified three multiply imaged galaxies in Hubble Space Telescope images of the redshift z=0.68 cluster responsible for the large-separation quadruply lensed quasar, SDSS J1004+4112. Spectroscopic redshifts have been secured for two of these systems using the Keck I 10m telescope. The most distant lensed galaxy, at z=3.332, forms at least four images, and an Einstein ring encompassing 3.1 times more area than the Einstein ring of the lensed QSO images at z=1.74, due to the greater source distance. For a second multiply imaged galaxy, we identify Ly_alpha emission at a redshift of z=2.74. The cluster mass profile can be constrained from near the center of the brightest cluster galaxy, where we observe both a radial arc and the fifth image of the lensed quasar, to the Einstein radius of the highest redshift galaxy, ~110 kpc. Our preliminary modeling indicates that the mass approximates an elliptical body, with an average projected logarithmic gradient of ~-0.5. The system is potentially useful for a direct measurement of world models in a previously untested redshift range.

We measure the topology of volume-limited galaxy samples selected from a parent sample of 314,050 galaxies in the Sloan Digital Sky Survey(SDSS), which is now complete enough to describe the fully three-dimensional topology and its dependence on galaxy properties. We compare the observed genus statistic G\nu_f) to predictions for a Gaussian random field and to the genus measured for mock surveys constructed from new large-volume simulations of the LambdaCDM cosmology. In this analysis we carefully examine the dependence of the observed genus statistic on the Gaussian smoothing scale R_G from 3.5 to 11 h^-1Mpc and on the luminosity of galaxies over the range -22.50<M_r<-18.5. The void multiplicity A_V is less than unity at all smoothing scales. Because A_V cannot become less than 1 through gravitational evolution, this result provides strong evidence for biased galaxy formation in low density environments. We also find clear evidence of luminosity bias of topology within the volume-limited sub-samples. The shift parameter Delta nu indicates that the genus of brighter galaxies shows a negative shift toward a ``meatball'' (i.e. cluster-dominated) topology, while faint galaxies show a positive shift toward a ``bubble'' (i.e. void-dominated) topology. The transition from negative to positive shift occurs approximately at the characteristic absolute magnitude M_{r*}=-20.4. Even in this analysis of the largest galaxy sample to date, we detect the influence of individual large-scale structures, as the shift parameter Delta nu and cluster multiplicity A_C reflect (at ~3 sigma) the presence of the Sloan Great Wall and a X-shaped structure which runs for several hundred Mpc across the survey volume.

We present an analytic expression for the axis ratio distribution of triaxial dark matter halos driven from physical principles. Adopting the picture of triaxial collapse based on the Zel'dovich approximation, we derive analytically both the minor-to-major and the conditional intermediate-to-major axis ratio distributions, and examine how they depend on the halo mass, redshift, and cosmology. Our analytic model is tested against the simulation data given by Jing & Suto in 2002, and found to reproduce the conditional intermediate-to-major axis ratio distribution successfully and the minor-to-major axis ratio distribution approximately. However, the trends of our analytic axis-ratio distributions with mass and redshift are opposite to what is found in N-body simulations. This failure of our analytic model puts a limitation on analytic approaches based on the Lagrangian theory to the halo ellipticity. Given the overall agreement with the simulation results, our model provides a new theoretical step toward using the axis-ratio distribution of dark halos as a cosmological probe. We also discuss several possibilities to improve the model.

We present a first morphometric investigation of a preliminary sample from the SDSS of 154287 galaxies with apparent magnitude 14.5<m_r<17.5 and redshift 0.001<z<0.4. We measure the Minkowski Functionals, which are a complete set of morphological descriptors. To account for the complicated wedge--like geometry of the present survey data, we construct isodensity contour surfaces from the galaxy positions in redshift space and employ two complementary methods of computing the Minkowski Functionals. We find that the observed Minkowski Functionals for SDSS galaxies are consistent with the prediction of a Lambda--dominated spatially--flat Cold Dark Matter model with random--Gaussian initial conditions, within the cosmic variance estimated from the corresponding mock catalogue. We expect that future releases of the SDSS survey will allow us to distinguish morphological differences in the galaxy distribution with regard to different morphological type and luminosity ranges.

We explore a signature of phase correlations in Fourier modes of dark matter density fields induced by nonlinear gravitational clustering. We compute the distribution function of the phase sum of the Fourier modes, theta_k1+theta_k2+theta_k3, for triangle wavevectors satisfying k1+k2+k3=0, and compare with the analytic prediction in perturbation theory recently derived by one of us. Using a series of cosmological N-body simulations, we extensively examine the time evolution and the dependence on the configuration of triangles and the sampling volume. Overall we find that the numerical results are remarkably consistent with the analytic formula from the perturbation theory. Interestingly the validity of the perturbation theory at a scale k is determined by P(k)/V_samp, the ratio of the power spectrum P(k) and the sampling volume V_samp, not by k^3P(k) as in the case of the conventional cosmological perturbation theory. Consequently this statistics of phase correlations is sensitive to the size of the sampling volume itself. This feature does not show up in more conventional cosmological statistics including the one-point density distribution function and the two-point correlation functions except as a sample-to-sample variation. Similarly if the sampling volume size V_samp is fixed, the stronger phase correlation emerges first at the wavevector where P(k) becomes largest, i.e., in linear regimes according to the standard cosmological perturbation theory, while the distribution of the phase sum stays fairly uniform in nonlinear regimes. The above feature can be naturally understood from the corresponding density structures in real space as we discuss in detail.

We report the first result of the clustering analysis of Sloan Digital Sky Survey (SDSS) quasars. We compute the two-point correlation function (2PCF) of SDSS quasars in redshift space at $8h^{-1}{\rm Mpc} < s < 500h^{-1}{\rm Mpc}$, with particular attention to its baryonic signature. Our sample consists of 19986 quasars extracted from the SDSS Data Release 4 (DR4). The redshift range of the sample is $0.72 \le z \le 2.24$ (the mean redshift is $\bar z = 1.46$) and the reddening-corrected $i$-band apparent magnitude range is $15.0 \le m_{i,{\rm rc}} \le 19.1$. Due to the relatively low number density of the quasar sample, the bump in the power spectrum due to the baryon density, $\Omega_{\rm b}$, is not clearly visible. The effect of the baryon density is, however, to distort the overall shape of the 2PCF.The degree of distortion makes it an interesting alternate measure of the baryonic signature. Assuming a scale-independent linear bias and the spatially flat universe, i.e., $\Omega_{\rm b} + \Omega_{\rm d} + \Omega_\Lambda =1$, where $\Omega_{\rm d}$ and $\Omega_\Lambda$ denote the density parameters of dark matter and the cosmological constant, we combine the observed quasar 2PCF and the predicted matter 2PCF to put constraints on $\Omega_{\rm b}$ and $\Omega_\Lambda$. Our result is fitted as $0.80- 2.8\Omega_{\rm b} < \Omega_\Lambda < 0.90 - 1.4\Omega_{\rm b}$ at the 2$\sigma$ confidence level, which is consistent with results from other cosmological observations such as WMAP. (abridged)

In spite of the growing observational evidence for dark matter and dark energy in the universe, their physical nature is largely unknown. In fact, several authors have proposed modifications of Newton's law of gravity at cosmological scales to account for the apparent acceleration of the cosmic expansion. Inspired by such suggestions, we attempt to constrain possible deviations from Newtonian gravity by means of the clustering of SDSS (Sloan Digital Sky Survey) galaxies. To be specific, we assume a simple modification of Newton's law with an additional Yukawa-type term characterized by the amplitude \alpha and the length scale \lambda. Adopting spatially-flat universes dominated by cold dark matter and/or dark energy, we solve a linear perturbation equation for the growth of density fluctuations. In particular, we find an exact analytic solution for the Einstein -- de Sitter case. Following the Peacock-Dodds prescription, we compute the nonlinear power spectra of mass fluctuations, perform a statistical comparison with the SDSS galaxy data, and derive constraints in the \alpha-\lambda plane; for instance, we obtain the constraints of -0.5<\alpha <0.6 and -0.8<\alpha<0.9 (99.7% confidence level) for \lambda=5h^{-1}Mpc and 10h^{-1}Mpc, respectively. We also discuss several future possibilities for improving our analysis.

A transiting extrasolar planet sequentially blocks off the light coming from the different parts of the disk of the host star in a time dependent manner. Due to the spin of the star, this produces an asymmetric distortion in the line profiles of the stellar spectrum, leading to an apparent anomaly of the radial velocity curves, known as the Rossiter - McLaughlin effect. Here, we derive approximate but accurate analytic formulae for the anomaly of radial velocity curves taking account of the stellar limb darkening. The formulae are particularly useful in extracting information of the projected angle between the planetary orbit axis and the stellar spin axis, \lambda, and the projected stellar spin velocity, V sin I_s. We create mock samples for the radial curves for the transiting extrasolar system HD209458, and demonstrate that constraints on the spin parameters (V sin I_s, \lambda) may be significantly improved by combining our analytic template formulae and the precision velocity curves from high-resolution spectroscopic observations with 8-10 m class telescopes. Thus future observational exploration of transiting systems using the Rossiter - McLaughlin effect is one of the most important probes to better understanding of the origin of extrasolar planetary systems, especially the origin of their angular momentum.

We present a simple model that accurately describes various statistical properties of peculiar velocities of dark matter halos. We pay particular attention to the following two effects; first, the evolution of the halo peculiar velocity depends on the local matter density, instead of the global density. Second, dark matter halos are biased tracers of the underlying mass distribution, thus halos tend to be located preferentially at high density regions. For the former, we develop an empirical model calibrated with N-body simulations, while for the latter, we use a conventional halo bias models based on the extended Press-Schechter model combined with an empirical log-normal probability distribution function of the mass density distribution. We find that compared with linear theory, the present model significantly improves the accuracy of predictions of statistical properties of the halo peculiar velocity field including the velocity dispersion, the probability distribution function, and the pairwise velocity dispersion at large separations. Thus our model predictions may be useful in analyzing future observations of the peculiar velocities of galaxy clusters.

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.

Galaxy redshift surveys have achieved significant progress over the last couple of decades. Those surveys tell us in the most straightforward way what our local universe looks like. While the galaxy distribution traces the bright side of the universe, detailed quantitative analyses of the data have even revealed the dark side of the universe dominated by non-baryonic dark matter as well as more mysterious dark energy (or Einstein's cosmological constant). We describe several methodologies of using galaxy redshift surveys as cosmological probes, and then summarize the recent results from the existing surveys. Finally we present our views on the future of redshift surveys in the era of Precision Cosmology.

We present measurements of the redshift--space three-point correlation function of galaxies in the Sloan Digital Sky Survey (SDSS). For the first time, we analyze the dependence of this statistic on galaxy morphology, color and luminosity. In order to control systematics due to selection effects, we used $r$--band, volume-limited samples of galaxies, constructed from the magnitude-limited SDSS data ($14.5<r<17.5$), and further divided the samples into two morphological types (early and late) or two color populations (red and blue). The three-point correlation function of SDSS galaxies follow the hierarchical relation well and the reduced three-point amplitudes in redshift--space are almost scale-independent ($Q_z=0.5\sim1.0$). In addition, their dependence on the morphology, color and luminosity is not statistically significant. Given the robust morphological, color and luminosity dependences of the two-point correlation function, this implies that galaxy biasing is complex on weakly non-linear to non-linear scales. We show that simple deterministic linear relation with the underlying mass could not explain our measurements on these scales.

There is evidence that the transiting planet HD 209458b has a large exosphere of neutral hydrogen, based on a 15% decrement in Lyman-alpha flux that was observed by Vidal-Madjar et al. during transits. Here we report upper limits on H-alpha absorption by the exosphere. The results are based on optical spectra of the parent star obtained with the Subaru High Dispersion Spectrograph. Comparison of the spectra taken inside and outside of transit reveals no exospheric H-alpha signal greater than 0.1% within a 5.1A band (chosen to have the same Delta_lambda/lambda as the 15% Ly-alpha absorption). The corresponding limit on the column density of n=2 neutral hydrogen is N_2 <~ 10^9 cm^{-2}. This limit constrains proposed models involving a hot (~10^4 K) and hydrodynamically escaping exosphere.

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

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

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

Statistics of lensed arcs in clusters of galaxies serve as a powerful probe of both the non-sphericity and the inner slope of dark matter halos. We develop a semi-analytic method to compute the number of arcs in triaxial dark matter halos. This combines the lensing cross section from the Monte Carlo ray-tracing simulations, and the probability distribution function (PDF) of the axis ratios evaluated from cosmological N-body simulations. This approach enables one to incorporate both asymmetries in the projected mass density and elongations along the line-of-sight analytically, for the first time in cosmological lensed arc statistics. As expected, triaxial dark matter halos significantly increase the number of arcs relative to spherical models; the difference amounts to more than one order of magnitude while the value of enhancement depends on the specific properties of density profiles. Then we compare our theoretical predictions with the observed number of arcs from 38 X-ray selected clusters. In contrast to the previous claims, our triaxial dark matter halos with inner density profile \rho \propto r^{-1.5} in a Lambda-dominated cold dark matter (CDM) universe reproduces well the observation. Since both the central mass concentration and triaxial axis ratios (minor to major axis ratio ~0.5) required to account for the observed data are consistent with cosmological N-body simulations, our result may be interpreted to lend strong support for the collisionless CDM paradigm at the mass scale of clusters.

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.

Simulations of large-scale structure in the universe have played a vital role in observational cosmology since 1980's in particular. Their important role will definitely continue to be true in the 21st century. Rather the requirements for simulations in the precision cosmology era will become more progressively demanding; they are supposed to fill the missing link in an accurate and reliable manner between the ``initial'' condition at z=1000 revealed by WMAP and the galaxy/quasar distribution at z=0 - 6 surveyed by 2dF and SDSS. In this review, I will summarize what we have learned so far from the previous cosmological simulations, and discuss several remaining problems for the new millennium.

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.

While the use of galaxy clusters as {\it tools} to probe cosmology is established, their conventional description still relies on the spherical and/or isothermal models that were proposed more than 20 years ago. We present, instead, a deprojection method to extract their intrinsic properties from X-ray and Sunyaev--Zel'dovich effect observations in order to improve our understanding of cluster physics. First we develop a theoretical model for the intra-cluster gas in hydrostatic equilibrium in a triaxial dark matter halo with a constant axis ratio. In this theoretical model, the gas density profiles are expressed in terms of the intrinsic properties of the dark matter halos. Then, we incorporate the projection effect into the gas profiles, and show that the gas surface brightness profiles are expressed in terms of the eccentricities and the orientation angles of the dark halos. For the practical purpose of our theoretical model, we provide several empirical fitting formulae for the gas density and temperature profiles, and also for the surface brightness profiles relevant to X-ray and Sunyaev--Zel'dovich effect observations. Finally, we construct a numerical algorithm to determine the halo eccentricities and orientation angles using our model, and demonstrate that it is possible in principle to reconstruct the 3D structures of the dark halos from the X-ray and/or Sunyaev-Zel'dovich effect cluster data alone without requiring priors such as weak lensing informations and without relying on such restrictive assumptions as the halo axial symmetry about the line-of-sight.

We explore empirical relations between three different populations of high-redshift galaxies and their hosting dark halos employing the halo model approach. Specifically we consider Lyman-break galaxies at z\sim4 and at z\sim5, and Lyman-Alpha emitters at z\simeq 4.86. We adopt a halo occupation function prescription to parameterize the properties of their hosting halos and the efficiency of halo-dependent star formation. We find that the two LBG samples are well described by the halo model with an appropriate HOF. We obtain constraints on properties of their hosting halos. A typical mass of hosting halos for LBGs is 5\times10^{11}h^{-1}M_\odot and the expected number of LBGs per halo is \sim0.5, therefore there is an approximate one-to-one correspondence between halos and LBGs. We also find a sign of the minimum mass of LBG hosting halos decreasing with time. We discuss implications of these findings on the star formation history of LBGs. On the other hand, for LAEs, our simple HOF prescription fails to reproduce simultaneously the observed angular correlation function and the number density. This might imply either that the distribution of LAEs within hosting halos differs from that of dark matter, or that the strong large-scale correlation is due to the existence of an unusual, large overdense region, and so the survey region is not a representative of the z\sim5 universe, the definite answer should wait for a much wider survey of LAEs at high redshifts.

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

We 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 present a numerical analysis of genus statistics for dark matter halo catalogs from the Hubble volume simulation. The huge box-size of the Hubble volume simulation enables us to carry out a reliable statistical analysis of the biasing properties of halos at a Gaussian smoothing scale of R_G>30Mpc/h with a cluster-mass scale of between 7*10^{13}Msolar/h and 6*10^{15}Msolar/h. A detailed comparison of the genus for dark matter halos with that for the mass distribution shows that the non-Gaussianity induced by the halo biasing is comparable to that by nonlinear gravitational evolution, and both the shape and the amplitude of the genus are almost insensitive to the halo mass at R_G>30Mpc/h. In order to characterize the biasing effect on the genus, we apply a perturbative formula developed by Matsubara (1994). We find that the perturbative formula well describes the simulated halo genus at R_G>50Mpc/h. The result indicates that the biasing effect on the halo genus is well approximated by nonlinear deterministic biasing up to the second-order term in the mass density fluctuation. The two parameters describing the linear and quadratic terms in the nonlinear biasing accurately specify the genus for galaxy clusters.

We present the first physical model for the non-spherical intra-cluster gas distribution in hydrostatic equilibrium under the gravity of triaxial dark matter halos. Adopting the concentric triaxial density profiles of the dark halos with constant axis ratios proposed by Jing & Suto (2002), we derive an analytical expression for the triaxial halo potential on the basis of the perturbation theory, and find the hydrostatic solutions for the gas density and temperature profiles both in isothermal and polytropic equations of state. The resulting iso-potential surfaces are well approximated by triaxial ellipsoids with the eccentricities dependent on the radial distance. We also find a formula for the eccentricity ratio between the intra-cluster gas and the underlying dark halo. Our results allow one to determine the shapes of the underlying dark halos from the observed intra-cluster gas through the X-ray and/or the Sunyaev-Zel'dovich effects clusters.

Strong gravitational lensing of distant supernovae (SNe), particularly Type Ia's, has some exploitable properties not available when other sorts of cosmologically distant sources are lensed. One such property is that the ``standard candle'' nature of the SN at peak brightness allows a direct determination of the lensing magnification factor for each well observed image. Another is that the duration of an SN event is of the same order as the differential time delays between the various lens images. The comparable time scales of the photometric event and the time delay invite a variety of applications, including high precision measurements of the delay and the targeting of especially interesting phases of the explosion for intensive observation when they appear in trailing images. As an initial exploration of these possibilities we present calculations of SN lensing statistics in a ``concordance cosmology''. We emphasize magnification and time delay effects. Plausible SN surveys, such as the proposed SNAP space mission, would discover several to some tens of strongly lensed SNe Ia per year, and at least a few of these will be at redshifts well beyond those that would be accessible via unlensed events. The total number of such anomalously high redshift SNe Ia will be a useful test of high redshift star formation models. SN surveys of finite duration will, of course, miss the appearance of some images; we quantify this selection bias. Finally, we investigate how well the appearance of trailing images can be predicted based on various amounts of available information on the lensing event. Knowledge of the magnification factor for the leading image makes it possible to predict the appearance of a trailing image relatively accurately if the lens redshift is also known. (Abridged)

We present a detailed non-spherical modeling of dark matter halos on the basis of a combined analysis of the high-resolution halo simulations (12 halos with $N\sim 10^6$ particles within their virial radius) and the large cosmological simulations (5 realizations with $N=512^3$ particles in a $100h^{-1}$Mpc boxsize). The density profiles of those simulated halos are well approximated by a sequence of the concentric triaxial distribution with their axis directions being fairly aligned. We characterize the triaxial model quantitatively by generalizing the universal density profile which has previously been discussed only in the framework of the spherical model. We obtain a series of practically useful fitting formulae in applying the triaxial model; the mass and redshift dependence of the axis ratio, the mean of the concentration parameter, and the probability distribution functions of the the axis ratio and the concentration parameter. These accurate fitting formulae form a complete description of the triaxial density profiles of halos in Cold Dark Matter models. Our current description of the dark halos will be particularly useful in predicting a variety of nonsphericity effects, to a reasonably reliable degree, including the weak and strong lens statistics, the orbital evolution of galactic satellites and triaxiality of galactic halos, and the non-linear clustering of dark matter. In addition, this provides a useful framework for the non-spherical modeling of the intra-cluster gas, which is crucial in discussing the gas and temperature profiles of X-ray clusters and the Hubble constant estimated via the Sunyaev -- Zel'dovich effect.

We present the first analysis of three-dimensional genus statistics for the SDSS EDR galaxy sample. Due to the complicated survey volume and the selection function, analytic predictions of the genus statistics for this sample are not feasible, therefore we construct extensive mock catalogs from N-body simulations in order to compare the observed data with model predictions. This comparison allows us to evaluate the effects of a variety of observational systematics on the estimated genus for the SDSS sample, including the shape of the survey volume, the redshift distortion effect, and the radial selection function due to the magnitude limit. The observed genus for the SDSS EDR galaxy sample is consistent with that predicted by simulations of a $\Lambda$-dominated spatially-flat cold dark matter model. Standard ($\Omega_0=1$) cold dark matter model predictions do not match the observations. We discuss how future SDSS galaxy samples will yield improved estimates of the genus.

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.

Density profiles of cosmological virialized systems, or dark halos, have recently attracted much attention. I first present a brief historical review of numerical simulations to quantify the halo density profiles. Then I describe the latest results on the universal density profile and their observational confrontation. Finally I discuss a clustering model of those halos with particular emphasis on the cosmological light-cone effect.

We derive the pairwise peculiar velocity distribution function of dark matter particles applying the dark matter halo approach. Unlike the previous work, we do not assume a Gaussian velocity distribution function of dark matter in a single halo, but compute it self-consistently with the assumed density profile for dark matter halo. The resulting distribution function is well approximated by an exponential distribution which is consistent with the previous observational, numerical and theoretical results. We also compute the pairwise peculiar velocity dispersion for different density profiles, and provide a practical fitting formula. We apply an empirical biasing scheme into our model and present prediction for pairwise peculiar velocity dispersion of {\it galaxies}, and reproduce the previous results of simulations using our semi-analytical method.

We critically examine the reliability of the dark matter clustering in high-resolution cosmological $N$-body simulations on scales below the mean separation length of particles. The particle discreteness effect imposes the two fundamental limitations on those scales; the lack of the initial fluctuation power and the finite mass resolution. We address this problem applying the dark halo approach and are able to discuss separately how those two limitations affect the dark matter clustering in $N$-body simulations at early epochs. We find that limitations of the dark matter clustering are primarily determined by the mass of particles. By a detailed comparison with three major cosmological simulations, we also find that in order to reproduce a proper amplitude of the dark matter clustering on small scales, halos with a characteristic nonlinear mass, $M_{\rm NL}(z)$ defined by $\sigma_R(M_{\rm NL}; z)=1$, must be resolved in the simulation. This leads to a critical redshift $z_{\rm crit}$ determined by $M_{\rm NL}(z_{\rm crit})=n_{\rm halo} m_{\rm part}$ where $n_{\rm halo}$ is the number of particles necessary to resolve the typical nonlinear mass halo ($\sim 10$). We conclude that, at least as far as the two-point correlation functions are concerned, the dark matter clustering in high-resolution $N$-body simulations on scales below the mean particle separation is reliable down to the gravitational force resolution length only for $z<z_{\rm crit}$, while it is strongly affected by the finite mass resolution for $z>z_{\rm crit}$.

Understanding of the structure evolution in the Universe has been greatly advanced with a rapid progress in high-resolution cosmological simulations. In this contribution, we report a new set of cosmological simulations with $512^3$ ($\sim 1.4\times 10^8$) particles in a simulation box of $100\mpc$. With this unprecedented resolution, we have successfully overcomed the over-merging problem while retaining a cosmological volume size that is crucial for the statistical analysis. The simulations are being applied to studying important astrophysical problems that are sensitive to the simulation resolution.

We report on the measurement of the two-point correlation function and the pairwise peculiar velocity of galaxies in the IRAS PSCz survey. The real space two-point correlation function can be fitted to a power law $\xi(r) = (r_0/r)^{\gamma}$ with $\gamma=1.69$ and $r_0=3.70 \mpc$. The pairwise peculiar velocity dispersion $\sigma_{12}(r_p)$ is close to $400 \kms$ at $r_p=3\mpc$ and decreases to about $150 \kms$ at $r_p \approx 0.2 \mpc$. These values are significantly lower than those obtained from the Las Campanas Redshift Survey, but agree very well with the results of blue galaxies reported by the SDSS team later on. We have constructed mock samples from N-body simulations with a cluster-weighted bias and from the theoretically constructed GIF catalog. We find that the two-point correlation function of the mock galaxies can be brought into agreemnt with the observed result, but the model does not reduce the velocity dispersions of galaxies to the level measured in the PSCz data. Thus we conclude that the peculiar velocity dispersions of the PSCz galaxies require a biasing model which substantially reduces the peculiar velocity dispersion on small scales relative to their spatial clustering. The results imply that either the cosmogony model needs to be revised or the velocity bias is important for the velocity dispersion of the IRAS galaxies.

The distribution of differential time delays \Delta t between images produced by strong gravitational lensing contains information on the mass distributions in the lensing objects as well as on cosmological parameters such as H_0. We derive an explicit expression for the conditional probability distribution function of time delays P(\Delta t|\theta), given an image separation between multiple images \theta, and related statistics. We consider lensing halos described by the singular isothermal sphere (SIS) approximation and by its generalization as proposed by Navarro, Frenk, & White (NFW) which has a density profile \rho \propto r^{-\alpha} in the innermost region. The time delay distribution is very sensitive to these profiles; steeper inner slopes tend to produce larger time delays. For example, if H_0=70 km/s/Mpc, a \Lambda-dominated cosmology and a source redshift z_S=1.27 are assumed, lenses with \theta=5'' produce a time delay of \Delta t[yr]=1.5^{+1.7}_{-0.9}, 0.39^{+0.37}_{-0.22}, 0.15^{+0.11}_{-0.09}, and 0.071^{+0.054}_{-0.038} (50% confidence interval), for SIS, generalized NFW with \alpha=1.5, \alpha=1.0, and \alpha=0.5, respectively. At a fixed image separation, the time delay is determined by the difference in the lensing potential between the position of the two images, which typically occur at different impact parameters. Although the values of \Delta t are proportional to the inverse of H_0, P(\Delta t|\theta) is rather insensitive to all other cosmological model parameters, source redshifts, magnification biases and so on. A knowledge of P(\Delta t|\theta) will also be useful in designing the observing program of future large scale synoptic variability surveys and for evaluating possible selection biases operating against large splitting lens systems.

The well-known self-similar solution for the two-point correlation function of the density field is valid only in an Einstein-de Sitter universe. We attempt to extend the solution for non-Einstein-de Sitter universes. For this purpose we introduce an idea of quasi-self-similar evolution; this approach is based on the assumption that the evolution of the two-point correlation is a succession of stages of evolution, each of which spans a short enough period to be considered approximately self-similar. In addition we assume that clustering is stable on scales where a physically motivated `virialization condition' is satisfied. These assumptions lead to a definite prediction for the behavior of the two-point correlation function in the strongly nonlinear regime. We show that the prediction agrees well with N-body simulations in non-Einstein-de Sitter cases, and discuss some remaining problems.

We develop a phenomenological model to predict the clustering of dark matter halos on the light-cone by combining several existing theoretical models. Assuming that the velocity field of halos on large scales is approximated by linear theory, we propose an empirical prescription of a scale-, mass-, and time-dependence of halo biasing. We test our model against the Hubble Volume $N$-body simulation and examine its validity and limitations. We find a good agreement in two-point correlation functions of dark matter halos between the phenomenological model predictions and measurements from the simulation for $R>5h^{-1}$Mpc both in the real and redshift spaces. Although calibrated on the mass scale of groups and clusters and for redshifts up to $z\sim2$, the model is quite general and can be applied to a wider range of astrophysical objects, such as galaxies and quasars, if the relation between dark halos and visible objects is specified.

A phenomenological model for the clustering of dark matter halos on the light-cone is presented. In particular, an empirical prescription for the scale-, mass-, and time-dependence of halo biasing is described in detail. A comparison of the model predictions against the light-cone output from the Hubble Volume $N$-body simulation indicates that the present model is fairly accurate for scale above $\sim 5h^{-1}$Mpc. Then I argue that the practical limitation in applying this model comes from the fact that we have not yet fully understood what are clusters of galaxies, especially at high redshifts. This point of view may turn out to be too pessimistic after all, but should be kept in mind in attempting {\it precision cosmology} with clusters of galaxies.

We report on the measurement of the two-point correlation function, and the pairwise peculiar velocity of galaxies in the IRAS PSCz survey. We compute these statistics first in redshift space, and then obtain the projected functions which have simple relations to the real-space correlation functions on the basis of the method developed earlier in analyzing the Las Campanas Redshift Survey (LCRS) by Jing, Mo, & B\"orner (1998). We find that the real space two-point correlation function can be fitted to a power law $\xi(r) = (r_0/r)^{\gamma}$ with $\gamma=1.69$ and $r_0=3.70 \mpc$. The pairwise peculiar velocity dispersion $\sigma_{12}(r_p)$ is close to $400 \kms$ at $r_p=3\mpc$ and decreases to about $150 \kms$ at $r_p \approx 0.2 \mpc$. These values are significantly lower than those obtained from the LCRS. In order to understand the implications of those measurements on the galaxy biasing, we construct mock samples for a low density spatially-flat cold dark matter model ($\Omega_0 = 0.3$, $\lambda_0=0.7$, $\Gamma=0.2$, $\sigma_8=1$) using a set of high-resolution N-body simulations. Applying a stronger cluster-underweight biasing ($\propto M^{-0.25}$) than for the LCRS ($\propto M^{-0.08}$), we are able to reproduce these observational data, except for the strong decrease of the pairwise peculiar velocity at small separations. This is qualitatively ascribed to the different morphological mixture of galaxies in the two catalogues. Disk-dominated galaxy samples drawn from the theoretically constructed GIF catalog yield results rather similar to our mock samples with the simple cluster-underweight biasing.

We quantitatively study the probability distribution function (PDF) of cosmological nonlinear density fluctuations from N-body simulations with Gaussian initial condition. In particular, we examine the validity and limitations of one-point and two-point log-normal PDF models against those directly estimated from the simulations. We find that the one-point log-normal PDF describes very accurately the cosmological density distribution even in the nonlinear regime (the rms variance \sigma_{nl} \simlt 4 and the over-density \delta \simlt 100). Furthermore the two-point log-normal PDFs are also in good agreement with the simulation data from linear to fairly nonlinear regime, while slightly deviate from them for \delta \simlt -0.5. Thus the log-normal PDF can be used as a useful empirical model for the cosmological density fluctuations. While this conclusion is fairly insensitive to the shape of the underlying power spectrum of density fluctuations P(k), models with substantial power on large scales, i.e., n\equiv d\ln P(k)/d \ln k \simlt -1, are better described by the log-normal PDF. On the other hand, we note that the one-to-one mapping of the initial and the evolved density fields consistent with the log-normal model does not approximate the broad distribution of their mutual correlation even on average. Thus the origin of the phenomenological log-normal PDF approximation still remains to be understood.

We report on the detailed analysis of the mean-pairwise peculiar velocity profile in high-resolution cosmological N-body simulations ($N=8.8\times 10^6$ particles in a sphere of $50 \sim 200$Mpc radius). In particular we examine the validity and limitations of the stable condition which states that the mean {\it physical} separation of particle pairs is constant on small scales. We find a significant time-variation (irregular oscillatory behavior) of the mean-pairwise peculiar velocity in nonlinear regimes. We argue that this behavior is not due to any numerical artifact, but a natural consequence of the continuous merging processes in the hierarchical clustering universe. While such a time-variation is significant in a relatively local patch of the universe, the global average over a huge spatial volume $\simgt$(200Mpc)$^3$ does not reveal any {\it systematic} departure from the stable condition. Thus we conclude that the mean-pairwise peculiar velocity is rather unstable statistics but still satisfies the stable condition when averaged over the cosmological volume.

We study the arc statistics of gravitational lensing generated by dark matter halos in order to probe their density profile. We characterize the halo profile by two parameters, the inner slope of the central cusp $\alpha$, and the median amplitude of the concentration parameter, $c_{\rm norm}$, for a halo of mass $10^{14}h^{-1}M_\odot$ at $z=0$, and compute the numbers of tangential and radial arcs produced by gravitational lensing of galaxy clusters. We find that the number of arcs divided by the number of halos is a good statistic which is sensitive to both $c_{\rm norm}$ and $\alpha$ with very weak dependence on the cosmological parameters. If the arc samples with well-defined selection criteria for the clusters become available, one can strongly constrain both $c_{\rm norm}$ and $\alpha$. While our tentative comparison with the existing observational data indicates that the inner density profile of dark halos is indeed as steep as predicted by recent simulations ($\alpha \sim 1.5$), the homogeneous samples of tangential and radial arcs are required for more quantitative discussions.

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

We derive an analytic formula of genus for dark matter halos assuming that the primordial mass density field obeys the random-Gaussian statistics. In particular, we for the first time take account of the nonlinear nature of the halo biasing, in addition to the nonlinear gravitational evolution of the underlying mass distribution. We examine in detail the extent to which the predicted genus for halos depends on the redshift, smoothing scale of the density field, and the mass of those halos in representative cold dark matter models. In addition to the full model predictions, we derive an explicit perturbation formula for the halo genus which can be applied to a wider class of biasing model for dark matter halos. With the empirical relation between the halo mass and the luminosity/temperature of galaxy clusters, our prediction can be compared with the upcoming cluster redshift catalogues in optical and X-ray bands. Furthermore, we discuss the detectability of the biasing and cosmological model dependence from the future data.

We quantify the degree of nonlinearity and stochasticity of the clustering of biased objects, using cosmological N-body simulations. Adopting the peaks and the halos as representative biasing models, we focus on the two-point correlation of the biased objects, dark matter and their cross-correlation. Especially, we take account of the effect of redshift-space distortion and attempt to clarify the scale-dependence and the time-dependence by analyzing the biasing factor and the cross-correlation factor. On small scales, stochasticity and nonlinearity become appreciable and strongly object-dependent, especially in redshift space due to the pair-wise velocity dispersion of the biased objects. Nevertheless, an approximation of deterministic linear biasing $\delta_{\rm obj}\simeq b_{\rm obj} \delta_{\rm mass}$ works reasonably well even in the quasi-linear regime $r > 10 h^{-1}$ Mpc, and linear redshift-space distortion explains the clustering amplitudes in redshift space in this regime.

We examine the light-cone effect on the two-point correlation functions using numerical simulations for the first time. Specifically, we generate several sets of dark matter particle distributions on the light-cone up to z=0.4 and z=2 over the field-of-view of \pi degree^2 from cosmological N-body simulations. Then we apply the selection function to the dark matter distribution according to the galaxy and QSO luminosity functions. Finally we compute the two-point correlation functions on the light-cone both in real and in redshift spaces using the pair-count estimator and compare with the theoretical predictions. We find that the previous theoretical modeling for nonlinear gravitational evolution, linear and nonlinear redshift-distortion, and the light-cone effect including the selection function is in good agreement with our numerical results, and thus is an accurate and reliable description of the clustering in the universe on the light-cone.

I critically review a methodology of using clusters of galaxies as cosmological probes. The understanding of the abundances and spatial correlations of dark matter halos has been significantly advanced especially for a last few years. Nevertheless such dark matter halos are not necessarily identical to clusters of galaxies in the universe. This is a quite obvious fact, but it seems that the resulting systematic errors, on the cosmological parameter estimates for instance, are often neglected in the literature.

We report on a wide-area (48' x 9') imaging survey of faint galaxies in R and I bands toward the 1338+27 field where an unusual concentration of five QSOs at z\sim1.1, embedded in a larger-scale clustering of 23 QSOs, is known to exist. Using a quite homogeneous galaxy catalog with a detection completeness limit of I~23.5, we detect a significant clustering signature of faint red galaxies with I > 21 and R-I > 1.2 over a scale extending to ~20h^{-1}_{50} Mpc. Close examination of the color-magnitude diagram, the luminosity function, and the angular correlation function indeed suggests that those galaxies are located at z~1.1 and trace the underlying large-scale structure at that epoch, together with the group of 5 QSOs. Since the whole extent of the cluster of 23 QSOs (~70h^{-1}_{50} Mpc) is roughly similar to the local ``Great Wall'', the area may contain a high-redshift counterpart of superclusters in the local universe.

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

We propose a physical model for nonlinear stochastic biasing of one-point statistics resulting from the formation epoch distribution of dark halos. In contrast to previous works on the basis of extensive numerical simulations, our model provides for the first time an analytic expression for the joint probability function. Specifically we derive the joint probability function of halo and mass density contrasts from the extended Press-Schechter theory. Since this function is derived in the framework of the standard gravitational instability theory assuming the random-Gaussianity of the primordial density field alone, we expect that the basic features of the nonlinear and stochastic biasing predicted from our model are fairly generic. As representative examples, we compute the various biasing parameters in cold dark matter models as a function of a redshift and a smoothing length. Our major findings are (1) the biasing of the variance evolves strongly as redshift while its scale-dependence is generally weak and a simple linear biasing model provides a reasonable approximation roughly at $R\simgt 2(1+z)\himpc$, and (2) the stochasticity exhibits moderate scale-dependence especially on $R\simlt 20\himpc$, but is almost independent of $z$. Comparison with the previous numerical simulations shows good agreement with the above behavior, indicating that the nonlinear and stochastic nature of the halo biasing is essentially understood by taking account of the distribution of the halo mass and the formation epoch.

We summarize a series of our recent work concerning the cosmological redshift-space distortion and light-cone effects. After briefly describing the theoretical formalism, we show how those effects are sensitive to the cosmological parameters. Then we apply this formalism to predict the two-point correlation functions and power spectra for the X-ray clusters, galaxies and quasars in future surveys, and discuss their cosmological implications.

We report a series of high-resolution N-body simulations designed to examine the internal physical properties of dark matter halos. A total of fifteen halos, each represented by $\sim 1$ million particles within the virial radius, have been simulated covering the mass range of $2\times10^{12}\sim 5\times10^{14} h^{-1}{\rm M_\odot}$. As the first application of these simulations, we have examined the density profiles of the halos. We found a clear systematic correlation between the halo mass and the slope of the density profile at one percent of the virial radius, in addition to the variations of the slope among halos of the similar mass. More specifically, the slope is $\sim - 1.5$, -1.3, and -1.1 for galaxy, group, and cluster mass halos, respectively. Thus we conclude that the dark matter density profiles, especially in the inner region, are not universal.

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

We present a theoretical formalism to predict the two-point clustering statistics (the power spectrum and the two-point correlation function), simultaneously taking account of the linear velocity distortion, the nonlinear velocity distortion (finger-of-god), the cosmological redshift-space distortion and the light-cone effect. To demonstrate the importance of these effects in exploring the clustering of objects at high redshifts, we show several model predictions for magnitude-limited surveys of galaxies and quasars. This methodology provides a quantitative tool to confront theoretical models against the upcoming precision data on clustering in the universe.

We perform a series of high-resolution N-body simulations designed to examine the density profiles of dark matter halos. From 12 simulated halos ranging the mass of $2\times10^{12}\sim 5\times10^{14} h^{-1}{\rm M_\odot}$ (represented by $\sim 1$ million particles within the virial radius), we find a clear systematic correlation between the halo mass and the slope of the density profile at 1% of the virial radius, in addition to the variations of the slope among halos of the similar mass. More specifically, the slope is $\sim -1.5$, -1.3, and -1.1 for galaxy, group, and cluster mass halos, respectively. While we confirm the earlier simulation results that the inner slope is steeper than the {\it universal} profile originally proposed by Navarro, Frenk & White, this mass dependence is inconsistent with the several analytical arguments attempting to link the inner slope with the primordial index of the fluctuation spectrum. Thus we conclude that the dark matter density profiles, especially in the inner region, are not universal.

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 examine the cosmological redshift-space distortion effect on the power spectrum of the objects at high-redshifts, which is an unavoidable observational contamination in general relativistic cosmology. In particular, we consider the nonlinear effects of density and velocity evolution using high-resolution N-body simulations in cold dark matter models. We find that the theoretical modeling on the basis of the fitting formulae of nonlinear density and velocity fields accurately describes the numerical results, especially in quasi-nonlinear regimes. These corrections for nonlinear effects are essential in order to use the the cosmological redshift-space distortion as a cosmological test. We perform a feasibility test to derive constraints from the future catalogues of high-redshift quasars using our theoretical modeling and results of N-body simulations. Applying the present methodology to the future data from on-going surveys of high-redshift galaxies and quasars will provide a useful tool to constrain a geometry of the universe.

We examine observational consequences of the cosmological light-cone effect on the power spectrum of the distribution of galaxies and quasars from upcoming redshift surveys. First we derive an expression for the power spectrum of cosmological objects in real space on a light cone, $P^{\rm LC}_{\rm R,lin}(k)$, which is exact in linear theory of density perturbations. Next we incorporate corrections for the nonlinear density evolution and redshift-space distortion in the formula in a phenomenological manner which is consistent with recent numerical simulations. On the basis of this formula, we predict the power spectrum of galaxies and quasars on the light cone for future redshift surveys taking account of the selection function properly. We demonstrate that this formula provides a reliable and useful method to compute the power spectrum on the light cone given an evolution model of bias.

The three-dimensional distribution of astronomical objects observed in redshift space significantly differs from the true distribution since the distance to each object cannot be determined by its redshift $z$ only; for $z \ll 1$ the peculiar velocity field contaminates the true recession velocity of the Hubble flow, while the true distance for objects at $z < 1$ sensitively depends on the (unknown and thus assumed) cosmological parameters. This hampers the effort to understand the true distribution of large-scale structure of the universe. In addition, all cosmological observations are carried out on a light-cone, the null hypersurface of an observer at $z=0$. This implies that their intrinsic properties and clustering statistics should change even within the survey volume. Therefore a proper comparison taking account of the light-cone effect is important to extract any cosmological information from redshift catalogues, especially for $z < 1$. We present recent theoretical development on the two effects -- the cosmological light-cone effect and the cosmological redshift-space distortion -- which should play key roles in observational cosmology in the 21st century.

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.

While all the cosmological observations are carried out on a light-cone, the null hypersurface of an observer at z=0, the clustering statistics has been properly defined only on the constant-time hypersurface. We develop a theoretical formulation for a two-point correlation function on the light-cone, and derive a practical approximate expression relevant to the discussion of clustering of high-redshift objects at large separations. As a specific example, we present predictions of the two-point correlation function for the Durham/AAT, SDSS and 2dF quasar catalogues. We also briefly discuss the effects of adopted luminosity function, cosmological parameters and bias models on the correlation function on the light-cone.

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

We describe a theoretical framework to compute the cluster gas distribution in hydrostatic equilibrium embedded in a class of spherical dark matter halo potentials. Unlike the conventional isothermal $\beta$-model, the present method provides a physical basis to directly probe the shape of dark matter halo from the observed X-ray surface brightness and temperature profiles of clusters of galaxies. Specifically, we examine the extent to which the resulting gas density and X-ray surface brightness profiles are sensitive to the inner slope of the dark matter halo density and other more realistic effects including the self-gravity of the gas and the polytropic equation of state. We also discuss a practical strategy to apply the present methodology to the actual cluster profiles from future X-ray observations.

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

We perform a detailed analysis of the statistical significance of a concentration of Lyman break galaxies at $z \sim 3$ recently discovered by Steidel et al. (1997), using a series of N-body simulations with $N=256^3$ particles in a $(100\himpc)^3$ comoving box. While the observed number density of Lyman break galaxies at $z\sim3$ implies that they correspond to systems with dark matter halos of $\simlt 10^{12}M_\odot$, the resulting clustering of such objects on average is not strong enough to be reconciled with the concentration if it is fairly common; we predict one similar concentration approximately per ($6\sim 10$) fields in three representative cold dark matter models. Considering the current observational uncertainty of the frequency of such clustering at $z\sim3$, it would be premature to rule out the models, but the future spectroscopic surveys in a dozen fields could definitely challenge all the existing cosmological models a posteriori fitted to the $z=0$ universe.

We explore some of the effects that discreteness and two-body scattering may have on N-body simulations with ``realistic'' cosmological initial conditions. We use an identical subset of particles from the initial conditions for a $128^3$ Particle-Mesh (PM) calculation as the initial conditions for a variety P$^3$M and Tree code runs. We investigate the effect of mass resolution (the mean interparticle separation) since most ``high resolution'' codes only have high resolution in gravitational force. The phase-insensitive two--point statistics, such as the power spectrum (autocorrelation) are somewhat affected by these variations, but phase-sensitive statistics show greater differences. Results converge at the mean interparticle separation scale of the lowest mass-resolution code. As more particles are added, but the force resolution is held constant, the P$^3$M and the Tree runs agree more and more strongly with each other and with the PM run which had the same initial conditions. This shows high particle density is necessary for correct time evolution, since many different results cannot all be correct. However, they do not so converge to a PM run which continued the fluctuations to small scales. Our results show that ignoring them is a major source of error on comoving scales of the missing wavelengths. This can be resolved by putting in a high particle density. Since the codes never agree well on scales below the mean comoving interparticle separation, we find little justification for quantitative predictions on this scale. Some measures vary by 50%, but others can be off by a factor of three or more. Our results suggest possible problems with the density of galaxy halos, formation of early generation objects such as QSO absorber clouds, etc.

The X-ray cluster gas density distribution in hydrostatic equilibrium is computed from the universal density profile of the dark matter halo recently proposed by Navarro, Frenk and White (1996, 1997). If one assumes the isothermality, the resulting distribution is well approximated by the conventional $\beta$-model. We predict the core radius $r_{c}$, the $\beta$-parameter, and the X-ray luminosity of clusters as a function of the temperature $T_{X}$ of clusters in some representative cosmological models, and compare them with observations and results of numerical simulations. The predicted size of $r_{c}$ is a factor of ($3 \sim 10$) smaller than the average of observed values. If both the universal density profile and the hydrostatic equilibrium are reasonable approximation to the truth, then this suggests either that the previous X-ray observations systematically overestimate the core radius of gas densities in clusters of galaxies, or that some important physical mechanisms, which significantly increase the core radius, is still missing.

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 have evaluated a systematic effect on counts-in-cells analysis of deep, wide-field galaxy catalogues induced by the evolution of clustering within the survey volume. A multiplicative correction factor is explicitly presented, which can be applied after the higher order correlation functions have been extracted in the usual way, without taking into account the evolution. The general theory of this effect combined with the ansatz describing the non-linear evolution of clustering in simulations enables us to estimate the magnitude of the correction factor in different cosmologies. In a series of numerical calculations assuming an array of cold dark matter models, it is found that, as long as galaxies are unbiased tracers of underlying density field, the effect is relatively small ($ \simeq 10%$) for the shallow surveys ($z <0.2$), while it becomes significant (order of unity) in deep surveys ($z \sim1$). Depending on the scales of interest, the required correction is comparable to or smaller than the expected errors of on-going wide-field galaxy surveys such as the SDSS and 2dF. Therefore at present, the effect has to be taken into account for high precision measurements at very small scales only, while in future deep surveys it amounts to a significant correction.

The observed two-point correlation functions of galaxies in redshift space become anisotropic due to the geometry of the universe as well as due to the presence of the peculiar velocity field. On the basis of linear perturbation theory, we expand the induced anisotropies of the correlation functions with respect to the redshift $z$, and obtain analytic formulae to infer the deceleration parameter $q_0$, the density parameter $\Omega_0$ and the derivative of the bias parameter $d\ln b/dz$ at $z=0$ in terms of the observable statistical quantities. The present method does not require any assumption of the shape and amplitude of the underlying fluctuation spectrum, and thus can be applied to future redshift surveys of galaxies including the Sloan Digital Sky Survey. We also evaluate quantitatively the systematic error in estimating the value of $\beta_0 \equiv \Omega_0^{0.6}/b$ from a galaxy redshift survey on the basis of a conventional estimator for $\beta_0$ which neglects both the geometrical distortion effect and the time evolution of the parameter $\beta(z)$. If the magnitude limit of the survey is as faint as 18.5 (in B-band) as in the case of the Sloan Digital Sky Survey, the systematic error ranges between -20% and 10% depending on the cosmological parameters. Although such systematic errors are smaller than the statistical errors in the current surveys, they will dominate the expected statistical error for future surveys.

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.

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.

Constraints on cosmological models from strong gravitational lensing statistics are investigated. We pay particular attention to the role of the velocity function in the calculation of the lensing probability. The velocity function derived from the observed galaxy luminosity function, which is used in most previous work, is unable to predict the large separation lensing events. In this paper, we also use the Press-Schechter theory to construct a velocity function theoretically. Model predictions are compared with the observed velocity function and the HST snapshot survey. Comparison with the latter observation shows that the predictions based on the theoretical velocity function are consistent with the observed large separation events in COBE normalized low-density models, especially with a non-vanishing cosmological constant. Adopting the COBE normalization, however, we could find no model which simultaneously satisfies both the observed velocity function and the HST snapshot survey. We systematically investigate various uncertainties in the gravitational lensing statistics including finite core radius, the distance formula, magnification bias, and dust obscuration. The results are very sensitive to these effects as well as theoretical models for the velocity function, implying that current limits on the cosmological parameters should be interpreted with caution. Predictions for future surveys are also presented.

We performed a systematic and comprehensive estimate of the Hubble constant $H_0$ via the Sunyaev-Zel'dovich temperature decrement and the X-ray measurements for clusters up to the redshift $z=0.541$, with particular attention to its dependence on the density parameter $\Omega_0$, and the cosmological constant $\lambda_0$. The resulting values of $H_0$ are largely consistent with that derived from the Cepheid distance for nearby galaxies. We discuss how the angular diameter distance $d_A$ vs. $z$ diagram from the accurate SZ measurements of high-$z$ clusters can distinguish the values of $\Omega_0$ and $\lambda_0$.

We discuss the effect of the finite size of galaxies on estimating small-scale relative pairwise peculiar velocity dispersions from the cosmic virial theorem (CVT). Specifically we evaluate the effect by incorporating the finite core radius $r_c$ in the two-point correlation function of mass, i.e. $\xi_\rho(r) \propto (r+r_c)^{-\gamma}$ and the effective gravitational force softening $r_s$ on small scales. We analytically obtain the lowest-order correction term for $\gamma <2$ which is in quantitative agreement with the full numerical evaluation. With a nonzero $r_s$ and/or $r_c$ the cosmic virial theorem is no longer limited to the case of $\gamma<2$. We present accurate fitting formulae for the CVT predicted pairwise velocity dispersion for the case of $\gamma>2$. Compared with the idealistic point-mass approximation ($r_s=r_c=0$), the finite size effect can significantly reduce the small-scale velocity dispersions of galaxies at scales much larger than $r_s$ and $r_c$. Even without considering the finite size of galaxies, nonzero values for $r_c$ are generally expected, for instance, for cold dark matter (CDM) models with a scale-invariant primordial spectrum. For these CDM models, a reasonable force softening $r_s\le 100 \hikpc$ would have rather tiny effect. We present the CVT predictions for the small-scale pairwise velocity dispersion in the CDM models normalized by the COBE observation. The implication of our results for confrontation of observations of galaxy pair-wise velocity dispersions and theoretical predictions of the CVT is also discussed.

I discuss two cosmological tests to determine the cosmological density parameter \Omega_0 the cosmological constant \lambda_0, which make use of the anisotropy of the two-point correlation functions due to the peculiar velocity field and the cosmological redshift-space distortion.

We discuss the two potentially important effects which should be taken into account in the analysis of the cosmological redshift-space distortion especially at high redshifts; the effect of inhomogeneities in the light propagation and the evolution of bias. Although such inhomogeneities affect significantly the degree of the cosmological redshift distortion, one can determine the density parameter $\Omega_0$ from the distortion pattern at $z=0$, and the cosmological constant $\lambda_0$ from that at $z = (1\sim 2)$. In particular we show that in low-density ($\Omega_0\ll 1$) universes one can determine the value of $\lambda_0$ and even the evolution of the bias parameter almost insensitively to the degree of the clumpiness of the matter distribution in the universe.

We propose cosmological redshift-space distortion of correlation functions of galaxies and quasars as a probe of both the density parameter $\Omega_0$ and the cosmological constant $\lambda_0$. In particular, we show that redshift-space distortion of quasar correlation functions at $z\sim2$ can in principle set a constraint on the value of $\lambda_0$. This is in contrast to the popular analysis of galaxy correlation functions in redshift space which basically determines $\Omega_0^{0.6}/b$, where $b$ is the bias parameter, but is insensitive to $\lambda_0$. For specific applications, we present redshift-space distortion of correlation functions both in cold dark matter models and in power-law correlation function models, and discuss the extent to which one can discriminate between the different $\lambda_0$ models.

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

Since recent X-ray observations have revealed that most clusters of galaxies are surrounded by an X-ray emitting gaseous halo, it is reasonable to expect that the Local Group of galaxies has its own X-ray halo. We show that such a halo, with temperature $\sim 1\keV$ and column density $\sim O(10^{21}) {\rm cm}^{-2}$, is a possible source for the excess low-energy component in the X-ray background. The halo should also generate temperature anisotropies in the microwave background via the Sunyaev-Zel'dovich effect. Assuming an isothermal spherical halo with the above temperature and density, the amplitude of the induced quadrupole turns out to be comparable to the COBE data without violating the upper limit on the $y$-parameter.

Since recent X-ray observations have revealed that most clusters of galaxies are surrounded by an X-ray emitting gaseous halo, it is reasonable to expect that the Local Group of galaxies has its own X-ray halo. We show that such a halo, with temperature $\sim 1$keV and column density $\sim O(10^{21}) {\rm cm}^{-2}$, is a possible source for the excess low-energy component in the X-ray background. The halo should also generate temperature anisotropies in the microwave background via the Sunyaev-Zel'dovich effect. Assuming an isothermal spherical halo with the above temperature and density, the amplitude of the induced quadrupole turns out to be comparable to the COBE data without violating the upper limit on the $y$-parameter. The induced dipole is negligible compared to the peculiar velocity of the Local Group, and multipoles higher than quadrupole are generally much smaller than the observed ones. However non-sphericity and/or clumpiness of the halo will produce a stronger effect. Therefore the gaseous halo of the Local Group, if it exists, will affect the estimate of the primordial spectral index $n$ and the amplitude of the density fluctuations deduced from the {\sl COBE} data.

The genus statistics is studied using large $N$-body simulations for several cosmological models. We consider the effects of nonlinear gravitational evolution, smoothing the particle data in fully nonlinear regime, and the redshift-space distortion on the genus curve. Detailed comparison between the theoretical prediction in weakly nonlinear theory and the appropriate simulation results shows that the analytic formula describes the behavior of genus in weakly nonlinear regime fairly accurately. We also find that the redshift-space distortion on genus statistics is small in linear and weakly nonlinear regimes. We conclude that if weakly nonlinear theory and direct numerical simulations are combined, the normalized genus curve $G(\nu)/G(0)$ is a powerful tool to directly check the random-Gaussian paradigm of the origin of the large-scale structure in the universe.

We present a method to compute the probability distribution function of the (true) Hubble constant and the age of the universe, given the estimate of the Hubble constant in our nearby galaxy samples. Our method takes into account both the observational errors and the cosmic variance, and enables to quantitatively compute the constraints on the cosmological models. Based on the present local observation $H_0=80\pm17$ km/s/Mpc, the probability of $H_0<50$ km/s/Mpc is about 6\% for the Einstein -- de Sitter universe ($\Omega_0=1$) and 4\% for an open ($\Omega_0=0.2$) universe. These probabilities are reduced to 0.8\% and 0.03\%, respectively, if the accuracy of the observational uncertainty is improved within 10\%. Similar probabilistic constraints on $t_0$ are also discussed.

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.

We compute the angular two-point correlation functions of the gamma-ray bursts at cosmological distances. Since the gamma-ray burst emission mechanism is not yet established, we simply assume that the gamma-ray burst sources are associated with high-redshift galaxies in some way. Then on the basis of several simple models for the evolution of galaxy spatial correlations, we calculate the amplitude of angular two-point correlation functions on scales appropriate for the Compton Gamma Ray Observatory data. we find that in most cases the predicted correlations are difficult to detect with the current data rate and the angular resolution, but models in which the bursts preferentially occur at relatively low redshift ($z \simlt 0.5$) predict correlation amplitudes on $\theta \sim 5^\circ$ which will be marginally detectable with the Gamma Ray Observatory data in several years. If future observations detect a signal of angular correlations, it will imply useful information on the correlation of galaxies at high redshifts provided that the gamma-ray bursts are cosmological.

We have performed the most comprehensive predictions of the temperature fluctuations \dtt in the primeval isocurvature baryon models to see whether or not the models are consistent with the recent data on the cosmic microwave background anisotropies. More specifically, we computed the \dtt corresponding to the experimental set-up by the South-Pole and the Owens Valley experiments as well as the COBE satellite. The amplitudes of the predicted \dtt are normalized by means of the COBE 10$^\circ$ data. The resulting constraints on the models are presented on $n - \Omega_b$ plane in the case of $\lambda_0=1-\Omega_b$ (flat models) and $\lambda_0=0$ (open models), where $n$ is the primordial spectral index of entropy fluctuations and $\Omega_b$ is the present baryon density parameter. Our results imply that the PIB models cannot be reconciled with the current observations for any reasonable set of cosmological parameters.