In this paper, we predict the covariance matrices of both the power spectrum and the bispectrum, including full non-Gaussian contributions, redshift space distortions, linear bias effects and shot-noise corrections, using perturbation theory (PT). To quantify the redshift-space distortion effect, we focus mainly on the monopole and quadrupole components of both the power and bispectra. We, for the first time, compute the 5- and 6-point spectra to predict the cross-covariance between the power and bispectra, and the auto-covariance of the bispectrum in redshift space. We test the validity of our calculations by comparing them with the covariance matrices measured from the MultiDark-Patchy mock catalogues that are designed to reproduce the galaxy clustering measured from the Baryon Oscillation Spectroscopic Survey Data Release 12. We argue that the simple, leading-order perturbation theory works because the shot-noise corrections for the Patchy mocks are more dominant than other higher-order terms we ignore. In the meantime, we confirm some discrepancies in the comparison, especially of the cross-covariance. We discuss potential sources of such discrepancies. We also show that our PT model reproduces well the cumulative signal-to-noise of the power spectrum and the bispectrum as a function of maximum wavenumber, implying that our PT model captures successfully essential contributions to the covariance matrices.

We study an efficient way to enhance the measurability of the galaxy density and/or velocity power spectrum in redshift space. It is based on the angular decomposition with the tripolar spherical harmonic (TripoSH) basis and applicable even to galaxy distributions in wide-angle galaxy surveys. While nontrivial multipole-mode mixings are inevitable in the covariance of the Legendre decomposition coefficient commonly used in the small-angle power spectrum analysis, our analytic computation of the covariance of the TripoSH decomposition coefficient shows that such mixings are absent by virtue of high separability of the TripoSH basis, yielding the minimum variance. Via the simple signal-to-noise ratio assessment, we confirm that the detectability improvement by the TripoSH decomposition approach becomes more significant at higher multipole modes, and, e.g., the hexadecapole of the density power spectrum has two orders of magnitude improvement. The TripoSH decomposition approach is expected to be applied to not only currently available survey data but also forthcoming wide-angle one, and to bring about something new or much more accurate cosmological information.

We determine constraints on spatially-flat tilted dynamical dark energy XCDM and $\phi$CDM inflation models by analyzing Planck 2015 cosmic microwave background (CMB) anisotropy data and baryon acoustic oscillation (BAO) distance measurements. XCDM is a simple and widely used but physically inconsistent parameterization of dynamical dark energy, while the $\phi$CDM model is a physically consistent one in which a scalar field $\phi$ with an inverse power-law potential energy density powers the currently accelerating cosmological expansion. Both these models have one additional parameter compared to standard $\Lambda$CDM and both better fit the TT + lowP + lensing + BAO data than does the standard tilted flat-$\Lambda$CDM model, with $\Delta \chi^2 = -1.26\ (-1.60)$ for the XCDM ($\phi$CDM) model relative to the $\Lambda$CDM model. While this is a 1.1$\sigma$ (1.3$\sigma$) improvement over standard $\Lambda$CDM and so not significant, dynamical dark energy models cannot be ruled out. In addition, both dynamical dark energy models reduce the tension between the Planck 2015 CMB anisotropy and the weak lensing $\sigma_8$ constraints.

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

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 report the discovery of 28 quasars and 7 luminous galaxies at 5.7 $\le$ z $\le$ 7.0. This is the tenth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The total number of spectroscopically identified objects in SHELLQs has now grown to 93 high-z quasars, 31 high-z luminous galaxies, 16 [O III] emitters at z ~ 0.8, and 65 Galactic cool dwarfs (low-mass stars and brown dwarfs). These objects were found over 900 deg2, surveyed by HSC between 2014 March and 2018 January. The full quasar sample includes 18 objects with very strong and narrow Ly alpha emission, whose stacked spectrum is clearly different from that of other quasars or galaxies. While the stacked spectrum shows N V 1240 emission and resembles that of lower-z narrow-line quasars, the small Ly alpha width may suggest a significant contribution from the host galaxies. Thus these objects may be composites of quasars and star-forming galaxies.

To fully extract cosmological information from nonlinear galaxy distribution in redshift space, it is essential to include higher-order statistics beyond the two-point correlation function. In this paper, we propose a new decomposition formalism for computing the anisotropic bispectrum in redshift space and for measuring it from galaxy samples. Our formalism uses tri-polar spherical harmonic decomposition with zero total angular momentum to compress the 3D modes distribution in the redshift-space bispectrum. This approach preserves three fundamental properties of the Universe: statistical homogeneity, isotropy, and parity-symmetry, allowing us to efficiently separate the anisotropic signal induced by redshift-space distortions (RSDs) and the Alcock-Paczy\'{n}ski (AP) effect from the isotropic bispectrum. The relevant expansion coefficients in terms of the anisotropic signal are reduced to one multipole index $L$, and the $L> 0$ modes are induced only by the RSD or AP effects. Our formalism has two advantages: (1) we can make use of Fast Fourier Transforms (FFTs) to measure the bispectrum; (2) it gives a simple expression to correct for the survey geometry, i.e., the survey window function. As a demonstration, we measure the decomposed bispectrum from the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12, and, for the first time, present a $14\sigma$ detection of the anisotropic bispectrum in the $L=2$ mode.

The super-sample tidal effect carries information on long-wavelength fluctuations that we cannot measure directly. It arises from the mode-coupling between short-wavelength and long-wavelength perturbations beyond a finite region of a galaxy survey and violates statistical isotropy of observed galaxy power spectra. In this paper, we propose the use of bipolar spherical harmonic (BipoSH) decomposition formalism to characterize statistically anisotropic power spectra. Using the BipoSH formalism, we perform a comprehensive study of the effect of the super-sample tides on measurements of other cosmological distortions such as the redshift-space distortion (RSD) and Alcock-Paczy\'{n}ski (AP) effects by means of the Fisher information matrix formalism. We find that the BipoSH formalism can break parameter degeneracies among the super-sample tidal, RSD and AP effects, indicating that the super-sample tides have little impact on the measurements of the RSD and AP effects. We also show that the super-sample tides are detectable with an accuracy better than the $\Lambda$CDM prediction without impairing the accuracy of measurements of other anisotropies assuming a SPHEREx-like galaxy survey.

We report the weak lensing measurement of filaments between Sloan Digital Sky Survey (SDSS) III/Baryon Oscillation Spectroscopic Survey (BOSS) CMASS galaxy pairs at $z\sim0.55$, using the Subaru Hyper Suprime-Cam (HSC) first-year galaxy shape catalogue. Despite of the small overlap of $140$ deg$^2$ between these surveys we detect the filament lensing signal at 3.9$\sigma$ significance, which is the highest signal-to-noise lensing measurement of filaments between galaxy-scale halos at this redshift range. We derive a theoretical prediction and covariance using mock catalogues based on full-sky ray-tracing simulations. We find that the intrinsic scatter of filament properties and the fluctuations in large scale structure along the line-of-sight are the primary component of the covariance and the intrinsic shape noise from source galaxies no longer limits our lensing measurement. This fact demonstrates the statistical power of the HSC survey due to its deep observations and high number density of source galaxies. Our result is consistent with the theoretical prediction and supports the "thick" filament model. As the HSC survey area increases, we will be able to study detailed filament properties such as the dark matter distributions and redshift evolution of filaments.

Neutral hydrogen (HI) intensity mapping is capable of measuring redshift evolution of HI density parameter $\Omega_{\rm HI}$, which is an important parameter to understand structure formation in the post-reionization epoch. Future HI observation with Square Kilometre Array (SKA) can significantly improve constraints on the parameter. However, the observation of HI suffers from the contamination from extremely bright foreground emissions, and it is necessary to consider a signal validation method complementary to the measurement of HI auto power spectrum. In this work, we propose to take a cross correlation between 21cm line intensity map and convergence map reconstructed from Planck observation and forecast a constraint on the $\Omega_{\rm HI}$. We find that the SKA-Mid operated as single-dish mode has a sufficient capability to detect the cross correlation and constrain $\Omega_{\rm HI}$ with $10-20\%$ precision at wide range of redshifts $0.5<z<3$, when combined with the Planck constraints on the cosmological parameters.

We report the discovery of a quasar at z = 7.07, which was selected from the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. This quasar, HSC J124353.93+010038.5, has an order of magnitude lower luminosity than do the other known quasars at z > 7. The rest-frame ultraviolet absolute magnitude is M1450 = -24.13 +/- 0.08 mag and the bolometric luminosity is Lbol = (1.4 +/- 0.1) x 10^{46} erg/s. Its spectrum in the optical to near-infrared shows strong emission lines, and shows evidence for a fast gas outflow, as the C IV line is blueshifted and there is indication of broad absorption lines. The Mg II-based black hole mass is Mbh = (3.3 +/- 2.0) x 10^8 Msun, thus indicating a moderate mass accretion rate with an Eddington ratio 0.34 +/- 0.20. It is the first z > 7 quasar with sub-Eddington accretion, besides being the third most distant quasar, known to date. The luminosity and black hole mass are comparable to, or even lower than, those measured for the majority of low-z quasars discovered by the Sloan Digital Sky Survey, and thus this quasar likely represents a z > 7 counterpart to quasars commonly observed in the low-z universe.

We measure cosmic weak lensing shear power spectra with the Subaru Hyper Suprime-Cam (HSC) survey first-year shear catalog covering 137deg$^2$ of the sky. Thanks to the high effective galaxy number density of $\sim$17 arcmin$^{-2}$ even after conservative cuts such as magnitude cut of $i<24.5$ and photometric redshift cut of $0.3\leq z \leq 1.5$, we obtain a high significance measurement of the cosmic shear power spectra in 4 tomographic redshift bins, achieving a total signal-to-noise ratio of 16 in the multipole range $300 \leq \ell \leq 1900$. We carefully account for various uncertainties in our analysis including the intrinsic alignment of galaxies, scatters and biases in photometric redshifts, residual uncertainties in the shear measurement, and modeling of the matter power spectrum. The accuracy of our power spectrum measurement method as well as our analytic model of the covariance matrix are tested against realistic mock shear catalogs. For a flat $\Lambda$ cold dark matter ($\Lambda$CDM) model, we find $S_8\equiv \sigma_8(\Omega_{\rm m}/0.3)^\alpha=0.800^{+0.029}_{-0.028}$ for $\alpha=0.45$ ($S_8=0.780^{+0.030}_{-0.033}$ for $\alpha=0.5$) from our HSC tomographic cosmic shear analysis alone. In comparison with Planck cosmic microwave background constraints, our results prefer slightly lower values of $S_8$, although metrics such as the Bayesian evidence ratio test do not show significant evidence for discordance between these results. We study the effect of possible additional systematic errors that are unaccounted in our fiducial cosmic shear analysis, and find that they can shift the best-fit values of $S_8$ by up to $\sim 0.6\sigma$ in both directions. The full HSC survey data will contain several times more area, and will lead to significantly improved cosmological constraints.

In the present universe, magnetic fields exist with various strengths and on various scales. One possible origin of these cosmic magnetic fields is the primordial magnetic fields (PMFs) generated in the early universe. PMFs are considered to contribute to matter density evolution via Lorentz force and the thermal history of intergalactic medium (IGM) gas due to ambipolar diffusion. Therefore, information about PMFs should be included in the temperature anisotropy of the Cosmic Microwave Background through the thermal Sunyaev-Zel'dovich (tSZ) effect in IGM. In this article, given an initial power spectrum of PMFs, we show the spatial fluctuation of mass density and temperature of the IGM and tSZ angular power spectrum created by the PMFs. Finally, we find that the tSZ angular power spectrum induced by PMFs becomes significant on small scales, even with PMFs below the observational upper limit. Therefore, we conclude that the measurement of tSZ anisotropy on small scales will provide the most stringent constraint on PMFs.

The presence of ubiquitous magnetic fields in the universe is suggested from observations of radiation and cosmic ray from galaxies or the intergalactic medium (IGM). One possible origin of cosmic magnetic fields is the magnetogenesis in the primordial universe. Such magnetic fields are called primordial magnetic fields (PMFs), and are considered to affect the evolution of matter density fluctuations and the thermal history of the IGM gas. Hence the information of PMFs is expected to be imprinted on the anisotropies of the cosmic microwave background (CMB) through the thermal Sunyaev-Zel'dovich (tSZ) effect in the IGM. In this study, given an initial power spectrum of PMFs as $P(k)\propto B_{\rm 1Mpc}^2 k^{n_{B}}$, we calculate dynamical and thermal evolutions of the IGM under the influence of PMFs, and compute the resultant angular power spectrum of the Compton $y$-parameter on the sky. As a result, we find that two physical processes driven by PMFs dominantly determine the power spectrum of the Compton $y$-parameter; (i) the heating due to the ambipolar diffusion effectively works to increase the temperature and the ionization fraction, and (ii) the Lorentz force drastically enhances the density contrast just after the recombination epoch. These facts result in making the tSZ angular power spectrum induced by the PMFs more remarkable at $\ell >10^4$ than that by galaxy clusters even with $B_{\rm 1Mpc}=0.1$ nG and $n_{B}=-1.0$ because the contribution from galaxy clusters decreases with increasing $\ell$. The measurement of the tSZ angular power spectrum on high $\ell$ modes can provide the stringent constraint on PMFs.

We present new measurements of the quasar luminosity function (LF) at $z \sim 6$, over an unprecedentedly wide range of the rest-frame ultraviolet luminosity $M_{1450}$ from $-30$ to $-22$ mag. This is the fifth in a series of publications from the Subaru High-$z$ Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The LF was calculated with a complete sample of 110 quasars at $5.7 \le z \le 6.5$, which includes 48 SHELLQs quasars discovered over 650 deg$^2$, and 63 brighter quasars discovered by the Sloan Digital Sky Survey and the Canada-France-Hawaii Quasar Survey (including one overlapping object). This is the largest sample of $z \sim 6$ quasars with a well-defined selection function constructed to date, and has allowed us to detect significant flattening of the LF at its faint end. A double power-law function fit to the sample yields a faint-end slope $\alpha = -1.23^{+0.44}_{-0.34}$, a bright-end slope $\beta = -2.73^{+0.23}_{-0.31}$, a break magnitude $M_{1450}^* = -24.90^{+0.75}_{-0.90}$, and a characteristic space density $\Phi^* = 10.9^{+10.0}_{-6.8}$ Gpc$^{-3}$ mag$^{-1}$. Integrating this best-fit model over the range $-18 < M_{1450} < -30$ mag, quasars emit ionizing photons at the rate of $\dot{n}_{\rm ion} = 10^{48.8 \pm 0.1}$ s$^{-1}$ Mpc$^{-3}$ at $z = 6.0$. This is less than 10 % of the critical rate necessary to keep the intergalactic medium ionized, which indicates that quasars are not a major contributor to cosmic reionization.

We study Planck 2015 cosmic microwave background (CMB) anisotropy data using the energy density inhomogeneity power spectrum generated by quantum fluctuations during an early epoch of inflation in the non-flat XCDM model. Here dark energy is parameterized using a fluid with a negative equation of state parameter but with the speed of fluid acoustic inhomogeneities set to the speed of light. We use this simple parameterization of dynamical dark energy, that is relatively straightforward to use in a computation, in a first attempt to gain some insight into how dark energy dynamics and non-zero spatial curvature jointly affect the CMB anisotropy data constraints. Unlike earlier analyses of non-flat models, we use a physically consistent power spectrum for energy density inhomogeneities. We find that the Planck 2015 data in conjunction with baryon acoustic oscillation measurements are reasonably well fit by a closed XCDM model in which spatial curvature contributes a percent of the current cosmological energy density budget. In this model, the measured Hubble constant and non-relativistic matter density parameter are in good agreement with values determined using most other data. Depending on parameter values, the closed XCDM model has reduced power, relative to the tilted, spatially-flat $\Lambda$CDM case, and appears to partially alleviate the low multipole CMB temperature anisotropy deficit and can help partially reconcile the CMB anisotropy and weak lensing $\sigma_8$ constraints, at the expense of somewhat worsening the fit to higher multipole CMB temperature anisotropy data. However, the closed XCDM inflation model does not seem to improve the agreement much, if at all, compared to the closed $\Lambda$CDM inflation case, even though it has one more free parameter. Our results are interesting but tentative; a more thorough analysis is needed to properly gauge their significance.

We perform Markov chain Monte Carlo analyses to put constraints on the non-flat $\phi$CDM inflation model using Planck 2015 cosmic microwave background (CMB) anisotropy data and baryon acoustic oscillation distance measurements. The $\phi$CDM model is a consistent dynamical dark energy model in which the currently accelerating cosmological expansion is powered by a scalar field $\phi$ slowly rolling down an inverse power-law potential energy density. We also use a physically consistent power spectrum for energy density inhomogeneities in this non-flat model. We find that, like the closed-$\Lambda$CDM and closed-XCDM models, the closed-$\phi$CDM model provides a better fit to the lower multipole region of the CMB temperature anisotropy data compared to that provided by the tilted flat-$\Lambda$CDM model. Also, like the other closed models, this model reduces the tension between the Planck and the weak lensing $\sigma_8$ constraints. However, the higher multipole region of the CMB temperature anisotropy data are better fit by the tilted flat-$\Lambda$CDM model than by the closed models.

We study Planck 2015 cosmic microwave background (CMB) anisotropy data using the energy density inhomogeneity power spectrum generated by quantum fluctuations during an early epoch of inflation in the non-flat $\Lambda$CDM model. Unlike earlier analyses of non-flat models, which assumed an inconsistent power-law power spectrum of energy density inhomogeneities, we find that the Planck 2015 data alone, and also in conjunction with baryon acoustic oscillation measurements, are reasonably well fit by a closed $\Lambda$CDM model in which spatial curvature contributes a few percent of the current cosmological energy density budget. In this model, the measured Hubble constant and non-relativistic matter density parameter are in good agreement with values determined using most other data. Depending on parameter values, the closed $\Lambda$CDM model has reduced power, relative to the tilted, spatially-flat $\Lambda$CDM case, and can partially alleviate the low multipole CMB temperature anisotropy deficit and can help partially reconcile the CMB anisotropy and weak lensing $\sigma_8$ constraints, at the expense of somewhat worsening the fit to higher multipole CMB temperature anisotropy data. Our results are interesting but tentative; a more thorough analysis is needed to properly gauge their significance.

The 21-cm signal in the vicinity of the first stars is expected to reflect properties of the first stars. In this study we pay special attention to tracing the time evolution of the ionizing photons' escape fraction, which affects the distribution of neutral hydrogen, by performing radiation hydrodynamics (RHD) simulations resolving dense gas in a halo. We find that the radial profile of 21-cm differential brightness temperature is quite sensitive to the stellar and halo masses, which reflects the time evolution of the escape fraction. In the case of a less massive star, ionizing photons hardly escape from its host halo due to the absorption by dense halo gas, thus an deep 21-cm absorption feature at just outside the halo lasts a long time. Whereas photons from a massive star well working to heat the ambient intergalactic medium turn out to cause a spatially extended 21-cm emission signature. Although individual signals are found to be undetectable with the Square Kilometre Array, our analysis using the results from the RHD simulations indicates that the properties of the first stars are imprinted on the 21-cm global signal: its amplitude depends not only on the cosmic star formation rate density, but also on the typical mass of the first stars due to the stellar-mass-dependent heating rate. Thus, we suggest that the initial mass function of the first stars is an essential factor in understanding the global signal.

We report discovery of 41 new high-z quasars and luminous galaxies, which were spectroscopically identified at 5.7 < z < 6.9. This is the fourth in a series of papers from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, based on the deep multi-band imaging data collected by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. We selected the photometric candidates by a Bayesian probabilistic algorithm, and then carried out follow-up spectroscopy with the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous papers, we have now spectroscopically identified 137 extremely-red HSC sources over about 650 deg2, which include 64 high-z quasars, 24 high-z luminous galaxies, 6 [O III] emitters at z ~ 0.8, and 43 Galactic cool dwarfs (low-mass stars and brown dwarfs). The new quasars span the luminosity range from M1450 ~ -26 to -22 mag, and continue to populate a few magnitude lower luminosities than have been probed by previous wide-field surveys. In a companion paper, we derive the quasar luminosity function at z ~ 6 over an unprecedentedly wide range of M1450 ~ -28 to -21 mag, exploiting the SHELLQs and other survey outcomes.

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

We present the measurement of the kinematic Sunyaev-Zel'dovich (kSZ) effect in Fourier space, rather than in real space. We measure the density-weighted pairwise kSZ power spectrum, the first use of this promising approach, by cross-correlating a cleaned Cosmic Microwave Background (CMB) temperature map, which jointly uses both Planck Release 2 and Wilkinson Microwave Anisotropy Probe nine-year data, with the two galaxy samples, CMASS and LOWZ, derived fr om the Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 12. With the current data, we constrain the average optical depth $\tau$ multiplied by the ratio of the Hubble parameter at redshift $z$ and the present day, $E=H/H_0$; we find $\tau E = (3.95\pm1.62)\times10^{-5}$ for LOWZ and $\tau E = ( 1.25\pm 1.06)\times10^{-5}$ for CMASS, with the optimal angular radius of an aperture photometry filter to estimate the CMB temperature distortion associ ated with each galaxy. By repeating the pairwise kSZ power analysis for various aperture radii, we measure the optical depth as a function of aperture ra dii. While this analysis results in the kSZ signals with only evidence for a detection, ${\rm S/N}=2.54$ for LOWZ and $1.24$ for CMASS, the combination of future CMB and spectroscopic galaxy surveys should enable precision measurements. We estimate that the combination of CMB-S4 and data from DESI shoul d yield detections of the kSZ signal with ${\rm S/N}=70-100$, depending on the resolution of CMB-S4.

The recent GW170817 measurement favors the simplest dark energy models, such as a single scalar field. Quintessence models can be classified in two classes, freezing and thawing, depending on whether the equation of state decreases towards $-1$ or departs from it. In this paper we put observational constraints on the parameters governing the equations of state of tracking freezing, scaling freezing and thawing models using updated data, from the Planck 2015 release, joint light-curve analysis and baryonic acoustic oscillations. Because of the current tensions on the value of the Hubble parameter $H_0$, unlike previous authors, we let this parameter vary, which modifies significantly the results. Finally, we also derive constraints on neutrino masses in each of these scenarios.

We measure statistically anisotropic signatures imprinted in three-dimensional galaxy clustering using bipolar spherical harmonics (BipoSHs) in both Fourier space and configuration space. We then constrain a well-known quadrupolar anisotropy parameter $g_{2M}$ in the primordial power spectrum, parametrized by $P(\vec{k}) = \bar{P}(k) [ 1 + \sum_{M} g_{2M} Y_{2M}(\hat{k}) ]$, with $M$ determining the direction of the anisotropy. Such an anisotropic signal is easily contaminated by artificial asymmetries due to specific survey geometry. We precisely estimate the contaminated signal and finally subtract it from the data. Using the galaxy samples obtained by the Baryon Oscillation Spectroscopic Survey Data Release 12, we find no evidence for violation of statistical isotropy, $g_{2M}$ for all $M$ to be of zero within the $2\sigma$ level. The $g_{2M}$-type anisotropy can originate from the primordial curvature power spectrum involving a directional-dependent modulation $g_* (\hat{k} \cdot \hat{p})^2$. The bound on $g_{2M}$ is translated into $g_*$ as $-0.09 < g_* < 0.08$ with a $95\%$ confidence level when $\hat{p}$ is marginalized over.

Magnetic fields are ubiquitous in the Universe. They seem to be present at virtually all scales and all epochs. Yet, whether the fields on cosmological scales are of astrophysical or cosmological origin remains an open major problem. Here we focus on an astrophysical mechanism based on the photoionization of the intergalactic medium during the Epoch of Reionization. Building upon previous studies that depicted the physical mechanism around isolated sources of ionization, we present here an analytic model to estimate the level at which this mechanism contributed to the magnetization of the whole Universe, thanks to the distribution of sources, before and alongside early luminous structure formation. This model suggests that the Universe may be globally magnetized to the order of, at least, a few $10^{-20}$~G comoving (i.e. several $10^{-18}$~G during the Epoch of Reionization) by this mechanism, prior to any amplification process.

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

We present spectroscopic identification of 32 new quasars and luminous galaxies discovered at 5.7 < z < 6.8. This is the second in a series of papers presenting the results of the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the deep multi-band imaging data produced by the Hyper Suprime-Cam (HSC) Subaru Strategic Program survey. The photometric candidates were selected by a Bayesian probabilistic algorithm, and then observed with spectrographs on the Gran Telescopio Canarias and the Subaru Telescope. Combined with the sample presented in the previous paper, we have now identified 64 HSC sources over about 430 deg2, which include 33 high-z quasars, 14 high-z luminous galaxies, 2 [O III] emitters at z ~ 0.8, and 15 Galactic brown dwarfs. The new quasars have considerably lower luminosity (M1450 ~ -25 to -22 mag) than most of the previously known high-z quasars. Several of these quasars have luminous (> 10^(43) erg/s) and narrow (< 500 km/s) Ly alpha lines, and also a possible mini broad absorption line system of N V 1240 in the composite spectrum, which clearly separate them from typical quasars. On the other hand, the high-z galaxies have extremely high luminosity (M1450 ~ -24 to -22 mag) compared to other galaxies found at similar redshift. With the discovery of these new classes of objects, we are opening up new parameter spaces in the high-z Universe. Further survey observations and follow-up studies of the identified objects, including the construction of the quasar luminosity function at z ~ 6, are ongoing.

We present cosmological constraints on the scalar-tensor theory of gravity by analyzing the angular power spectrum data of the cosmic microwave background obtained from the Planck 2015 results together with the baryon acoustic oscillations (BAO) data. We find that the inclusion of the BAO data improves the constraints on the time variation of the effective gravitational constant by more than $10\%$, that is, the time variation of the effective gravitational constant between the recombination and the present epochs is constrained as $G_{\rm rec}/G_0-1 <1.9\times 10^{-3}\ (95.45\%\ {\rm C.L.})$ and $G_{\rm rec}/G_0-1 <5.5\times 10^{-3}\ (99.99 \%\ {\rm C.L.})$. We also discuss the dependence of the constraints on the choice of the prior.

We study the properties of dark matter halos that contain star-forming galaxies at $1.43 \le z \le 1.74$ using the FMOS-COSMOS survey. The sample consists of 516 objects with a detection of the H$\alpha$ emission line, that represent the star-forming population at this epoch having a stellar mass range of $10^{9.57}\le M_\ast/M_\odot \lesssim 10^{11.4}$ and a star formation rate range of $15\lesssim \mathrm{SFR}/(M_\odot \mathrm{yr^{-1}}) \lesssim 600$. We measure the projected two-point correlation function while carefully taking into account observational biases, and find a significant clustering amplitude at scales of $0.04$-$10~h^{-1}~\mathrm{cMpc}$, with a correlation length $r_0 = 5.21^{+0.70}_{-0.67}~h^{-1}~\mathrm{cMpc}$ and a bias $b=2.59^{+0.41}_{-0.34}$. We interpret our clustering measurement using a halo occupation distribution model. The sample galaxies appear to reside in halos with mass $M_\mathrm{h} = 4.6^{+1.1}_{-1.6}\times10^{12}~h^{-1}M_\odot$ on average that will likely become present-day halos of mass $M_\mathrm{h} (z=0) \sim2\times10^{13}~h^{-1}M_\odot$, equivalent to the typical halo mass scale of galaxy groups. We then confirm the decline of the stellar-to-halo mass ratio at $M_\mathrm{h}<10^{12}~M_\odot$, finding $M_\ast/M_\mathrm{h} \approx 5\times10^{-3}$ at $M_\mathrm{h}=10^{11.86}~M_\odot$, which is lower by a factor of 2-4 than those measured at higher masses. Finally, we use our results to illustrate the future capabilities of Subaru's Prime-Focus Spectrograph, a next-generation instrument that will provide strong constraints on the galaxy-formation scenario by obtaining precise measurements of galaxy clustering at $z>1$.

We propose an efficient way to test rotational invariance in the cosmological perturbations by use of galaxy correlation functions. In symmetry-breaking cases, the galaxy power spectrum can have extra angular dependence in addition to the usual one due to the redshift-space distortion, $ \hat{k} \cdot \hat{n}$. We confirm that, via the decomposition into not the usual Legendre basis ${\cal L}_\ell(\hat{k} \cdot \hat{n})$ but the bipolar spherical harmonic one $\{Y_{\ell}(\hat{k}) \otimes Y_{\ell'}(\hat{n})\}_{LM}$, the symmetry-breaking signal can be completely distinguished from the usual isotropic one since the former yields nonvanishing $L \geq 1$ modes but the latter is confined to the $L = 0$ one. As a demonstration, we analyze the signatures due to primordial-origin symmetry breakings such as the well-known quadrupolar-type and dipolar-type power asymmetries and find nonzero $L = 2$ and $1$ modes, respectively. Fisher matrix forecasts of their constraints indicate that the $Planck$-level sensitivity could be achieved by the SDSS or BOSS-CMASS data, and an order-of-magnitude improvement is expected in a near future survey as PFS or Euclid by virtue of an increase in accessible Fourier mode. Our methodology is model-independent and hence applicable to the searches for various types of statistically anisotropic fluctuations.

Yes. Future CMB experiments such as Advanced ACTPol and CMB-S4 should achieve measurements with S/N of $> 0.1$ for the typical galaxies in redshift surveys. These measurements will provide complementary measurements of the growth rate of large scale structure $f$ and the expansion rate of the Universe $H$ to galaxy clustering measurements. This paper emphasizes that there is significant information in the anisotropy of the relative pairwise kSZ measurements. We expand the relative pairwise kSZ power spectrum in Legendre polynomials and consider up to its octopole. Assuming that the noise in the filtered maps is uncorrelated between the positions of galaxies in the survey, we derive a simple analytic form for the power spectrum covariance of the relative pairwise kSZ temperature in redshift space. While many previous studies have assumed optimistically that the optical depth of the galaxies $\tau_{\rm T}$ in the survey is known, we marginalize over $\tau_{\rm T}$, to compute constraints on the growth rate $f$ and the expansion rate $H$. For realistic sure parameters, we find that combining kSZ and galaxy redshift survey data reduces the marginalized $1$-$\sigma$ errors on $H$ and $f$ by $\sim50$-$70\%$ compared to the galaxy-only analysis.

We investigate the physical conditions of ionized gas in high-z star-forming galaxies using diagnostic diagrams based on the rest-frame optical emission lines. The sample consists of 701 galaxies with an Ha detection at $1.4\lesssim z\lesssim1.7$, from the FMOS-COSMOS survey, that represent the normal star-forming population over the stellar mass range $10^{9.6} \lesssim M_\ast/M_\odot \lesssim 10^{11.6}$ with those at $M_\ast>10^{11}~M_\odot$ being well sampled. We confirm an offset of the average location of star-forming galaxies in the BPT diagram ([OIII]/Hb vs. [NII]/Ha), primarily towards higher [OIII]/Hb, compared with local galaxies. Based on the [SII] ratio, we measure an electron density ($n_e=220^{+170}_{-130}~\mathrm{cm^{-3}}$), that is higher than that of local galaxies. Based on comparisons to theoretical models, we argue that changes in emission-line ratios, including the offset in the BPT diagram, are caused by a higher ionization parameter both at fixed stellar mass and at fixed metallicity with additional contributions from a higher gas density and possibly a hardening of the ionizing radiation field. Ionization due to AGNs is ruled out as assessed with Chandra. As a consequence, we revisit the mass-metallicity relation using [NII]/Ha and a new calibration including [NII]/[SII] as recently introduced by Dopita et al. Consistent with our previous results, the most massive galaxies ($M_\ast\gtrsim10^{11}~M_\odot$) are fully enriched, while those at lower masses have metallicities lower than local galaxies. Finally, we demonstrate that the stellar masses, metallicities and star formation rates of the FMOS sample are well fit with a physically-motivated model for the chemical evolution of star-forming galaxies.

We study quasiequilibrium solutions of triaxially deformed rotating compact stars -- a generalization of Jacobi ellipsoids under relativistic gravity and compressible equations of state (EOS). For relatively stiff (piecewise) polytropic EOSs, we find supramassive triaxial solutions whose masses exceed the maximum mass of the spherical solution, but are always lower than those of axisymmetric equilibriums. The difference in the maximum masses of triaxial and axisymmetric solutions depends sensitively on the EOS. If the difference turns out to be only about 10%, it will be strong evidence that the EOS of high density matter becomes substantially softer in the core of neutron stars. This finding opens a novel way to probe phase transitions of high density nuclear matter using detections of gravitational waves from new born neutron stars or magnetars under fallback accretion.

We report the discovery of 15 quasars and bright galaxies at 5.7 < z < 6.9. This is the initial result from the Subaru High-z Exploration of Low-Luminosity Quasars (SHELLQs) project, which exploits the exquisite multiband imaging data produced by the Subaru Hyper Suprime-Cam (HSC) Strategic Program survey. The candidate selection is performed by combining several photometric approaches including a Bayesian probabilistic algorithm to reject stars and dwarfs. The spectroscopic identification was carried out with the Gran Telescopio Canarias and the Subaru Telescope for the first 80 deg2 of the survey footprint. The success rate of our photometric selection is quite high, approaching 100 % at the brighter magnitudes (zAB < 23.5 mag). Our selection also recovered all the known high-z quasars on the HSC images. Among the 15 discovered objects, six are likely quasars, while the other six with interstellar absorption lines and in some cases narrow emission lines are likely bright Lyman-break galaxies. The remaining three objects have weak continua and very strong and narrow Ly alpha lines, which may be excited by ultraviolet light from both young stars and quasars. These results indicate that we are starting to see the steep rise of the luminosity function of z > 6 galaxies, compared with that of quasars, at magnitudes fainter than M1450 ~ -22 mag or zAB ~24 mag. Follow-up studies of the discovered objects as well as further survey observations are ongoing.

Cosmological constraints on the scalar-tensor theory of gravity by analyzing the angular power spectrum data of the cosmic microwave background (CMB) obtained from the Planck 2015 results are presented. We consider the harmonic attractor model, in which the scalar field has a harmonic potential with curvature ($\beta$) in the Einstein frame and the theory relaxes toward the Einstein gravity with time. Analyzing the {\it TT}, {\it EE}, {\it TE} and lensing CMB data from Planck by the Markov chain Monte Carlo method, we find that the present-day deviation from the Einstein gravity (${\alpha_0}^2$) is constrained as ${\alpha_0}^2<2.5\times10^{-4-4.5\beta^2}\ (95.45\% {\rm\ C.L.})$ and ${\alpha_0}^2<6.3\times10^{-4-4.5\beta^2}\ (99.99\%\ {\rm C.L.})$ for $0<\beta<0.4$. The time variation of the effective gravitational constant between the recombination and the present epochs is constrained as $G_{\rm rec}/G_0<1.0056\ (95.45\% {\rm\ C.L.})$ and $G_{\rm rec}/G_0<1.0115\ (99.99 \%{\rm\ C.L.})$. We also find that the constraints are little affected by extending to nonflat cosmological models because the diffusion damping effect revealed by Planck breaks the degeneracy of the projection effect.

Observations of the kinetic Sunyaev-Zel'dovich (kSZ) effect measure the density-weighted velocity field, a potentially powerful cosmological probe. This paper presents an analytical method to predict the power spectrum and two-point correlation function of the density-weighted velocity in redshift space, the direct observables in kSZ surveys. We show a simple relation between the density power spectrum and the density-weighted velocity power spectrum that holds for both dark matter and halos. Using this relation, we can then extend familiar perturbation expansion techniques to the kSZ power spectrum. One of the most important features of density-weighted velocity statistics in redshift space is the change in sign of the cross-correlation between the density and density-weighted velocity at mildly small scales due to nonlinear redshift space distortions. Our model can explain this characteristic feature without any free parameters. As a result, our results can precisely predict the non-linear behavior of the density-weighted velocity field in redshift space up to $\sim30\ h^{-1} {\rm Mpc}$ for dark matter particles at the redshifts of $z=0.0$, $0.5$, and $1.0$.

Scattering of cosmic microwave background (CMB) radiation in galaxy clusters induces polarization signals determined by the quadrupole anisotropy in the photon distribution at the location of clusters. This "remote quadrupole" derived from the measurements of the induced polarization in galaxy clusters provides an opportunity of reconstruction of local CMB temperature anisotropies. In this {\em Letter} we develop an algorithm of the reconstruction through the estimation of the underlying primordial gravitational potential, which is the origin of the CMB temperature and polarization fluctuations and CMB induced polarization in galaxy clusters. We found a nice reconstruction for the quadrupole and octopole components of the CMB temperature anisotropies with the assistance of the CMB induced polarization signals. The reconstruction can be an important consistency test on the puzzles of CMB anomaly, especially for the low quadrupole and axis of evil problems reported in WMAP and Planck data.

Cosmic strings are a type of cosmic defect formed by a symmetry-breaking phase transition in the early universe. Individual strings would have gathered to build a network, and their dynamical motion would induce scalar--, vector-- and tensor--type perturbations. In this paper, we focus on the vector mode perturbations arising from the string network based on the one scale model and calculate the time evolution and the power spectrum of the associated magnetic fields. We show that the relative velocity between photon and baryon fluids induced by the string network can generate magnetic fields over a wide range of scales based on standard cosmology. We obtain the magnetic field spectrum before recombination as $a^2B(k,z)\sim4\times10^{-16}G\mu/((1+z)/1000)^{4.25}(k/{\rm Mpc}^{-1})^{3.5}$ Gauss on super-horizon scales, and $a^2B(k,z)\sim2.4\times10^{-17}G\mu/((1+z)/1000)^{3.5}(k/{\rm Mpc}^{-1})^{2.5}$ Gauss on sub-horizon scales in co-moving coordinates. This magnetic field grows up to the end of recombination, and has a final amplitude of approximately $B\sim10^{-17\sim -18} G\mu$ Gauss at the $k\sim1\ {\rm Mpc}^{-1}$ scale today. This field might serve as a seed for cosmological magnetic fields.

We present a spectroscopic survey of galaxies in the COSMOS field using the Fiber Multi-Object Spectrograph (FMOS), a near-infrared instrument on the Subaru Telescope. Our survey is specifically designed to detect the Halpha emission line that falls within the H-band (1.6-1.8 um) spectroscopic window from star-forming galaxies with 1.4 < z < 1.7 and M_stellar>~10^10 Msolar. With the high multiplex capability of FMOS, it is now feasible to construct samples of over one thousand galaxies having spectroscopic redshifts at epochs that were previously challenging. The high-resolution mode (R~2600) effectively separates Halpha and [NII]6585 thus enabling studies of the gas-phase metallicity and photoionization state of the interstellar medium. The primary aim of our program is to establish how star formation depends on stellar mass and environment, both recognized as drivers of galaxy evolution at lower redshifts. In addition to the main galaxy sample, our target selection places priority on those detected in the far-infrared by Herschel/PACS to assess the level of obscured star formation and investigate, in detail, outliers from the star formation rate - stellar mass relation. Galaxies with Halpha detections are followed up with FMOS observations at shorter wavelengths using the J-long (1.11-1.35 um) grating to detect Hbeta and [OIII]5008 that provides an assessment of extinction required to measure star formation rates not hampered by dust, and an indication of embedded Active Galactic Nuclei. With 460 redshifts measured from 1153 spectra, we assess the performance of the instrument with respect to achieving our goals, discuss inherent biases in the sample, and detail the emission-line properties. Our higher-level data products, including catalogs and spectra, are available to the community.

If vector type perturbations are present in the primordial plasma before recombination, the generation of magnetic fields is known to be inevitable through the Harrison mechanism. In the context of the standard cosmological perturbation theory, non-linear couplings of first-order scalar perturbations create second-order vector perturbations, which generate magnetic fields. Here we reinvestigate the generation of magnetic fields at second-order in cosmological perturbations on the basis of our previous study, and extend it by newly taking into account the time evolution of purely second-order vector perturbations with a newly developed second-order Boltzmann code. We confirm that the amplitude of magnetic fields from the product-terms of the first-order scalar modes is consistent with the result in our previous study. However, we find, both numerically and analytically, that the magnetic fields from the purely second-order vector perturbations partially cancel out the magnetic fields from one of the product-terms of the first-order scalar modes, in the tight coupling regime in the radiation dominated era. Therefore, the amplitude of the magnetic fields on small scales, $k \gtrsim 10~h{\rm Mpc}^{-1}$, is smaller than the previous estimates. The amplitude of the generated magnetic fields at cosmological recombination is about $B_{\rm rec} =5.0\times 10^{-24}~{\rm Gauss}$ on $k = 5.0 \times 10^{-1}~h{\rm Mpc}^{-1}$. Finally, we discuss the reason of the discrepancies that exist in estimates of the amplitude of magnetic fields among other authors.

A symmetry-breaking phase transition in the early universe could have led to the formation of cosmic defects. Because these defects dynamically excite not only scalar and tensor type cosmological perturbations but also vector type ones, they may serve as a source of primordial magnetic fields. In this study, we calculate the time evolution and the spectrum of magnetic fields that are generated by a type of cosmic defects, called global textures, using the non-linear sigma (NLSM) model. Based on the standard cosmological perturbation theory, we show, both analytically and numerically, that a vector-mode relative velocity between photon and baryon fluids is induced by textures, which inevitably leads to the generation of magnetic fields over a wide range of scales. We find that the amplitude of the magnetic fields is given by $B\sim{10^{-9}}{((1+z)/10^3)^{-2.5}}({v}/{m_{\rm pl}})^2({k}/{\rm Mpc^{-1}})^{3.5}/{\sqrt{N}}$ Gauss in the radiation dominated era for $k\lesssim 1$ Mpc$^{-1}$, with $v$ being the vacuum expectation value of the O(N) symmetric scalar fields. By extrapolating our numerical result toward smaller scales, we expect that $B\sim {10^{-17}}((1+z)/1000)^{-1/2}({v}/{m_{\rm pl}})^2({k}/{\rm Mpc^{-1}})^{1/2}/{\sqrt{N}}$ Gauss on scales of $k\gtrsim 1$ Mpc$^{-1}$ at redshift $z\gtrsim 1100$. This might be a seed of the magnetic fields observed on large scales today.

Gravitational waves (GWs) are inevitably induced at second-order in cosmological perturbations through non-linear couplings with first order scalar perturbations, whose existence is well established by recent cosmological observations. So far, the evolution and the spectrum of the secondary induced GWs have been derived by taking into account the sources of GWs only from the product of first order scalar perturbations. Here we newly investigate the effects of purely second-order anisotropic stresses of photons and neutrinos on the evolution of GWs, which have been omitted in the literature. We present a full treatment of the Einstein-Boltzmann system to calculate the spectrum of GWs with anisotropic stress based on the formalism of the cosmological perturbation theory. We find that photon anisotropic stress amplifies the amplitude of GWs by about $150 %$ whereas neutrino anisotropic stress suppress that of GWs by about $30 %$ on small scales $k\gtrsim 1.0 h{\rm Mpc}^{-1}$ compared to the case without anisotropic stress. The second order anisotropic stress does not affect GWs with wavenumbers $k\lesssim 1.0 h{\rm Mpc}^{-1}$. The result is in marked contrast with the case at linear order, where the effect of anisotropic stress is damping in amplitude of GWs.

The prospects for direct measurements of inflationary gravitational waves by next generation interferometric detectors inferred from the possible detection of B-mode polarization of the cosmic microwave background are studied. We compute the spectra of the gravitational wave background and the signal-to-noise ratios by two interferometric detectors (DECIGO and BBO) for large-field inflationary models in which the tensor-to-scalar ratio is greater than the order of 0.01. If the reheating temperature $T_{\rm RH}$ of chaotic inflation with the quadratic potential is high ($T_{\rm RH}>7.9\times10^6$ GeV for upgraded DECIGO and $T_{\rm RH}> 1.8\times 10^{6}$ GeV for BBO), it will be possible to reach the sensitivity of the gravitational background in future experiments at $3\sigma$ confidence level. The direct detection is also possible for natural inflation with the potential $V(\phi)=\Lambda^4 [1-\cos(\phi/f)]$, provided that $f>4.2 M_{\rm pl}$ (upgraded DECIGO) and $f>3.6 M_{\rm pl}$ (BBO) with $T_{\rm RH}$ higher than $10^8$ GeV. The quartic potential $V(\phi)=\lambda \phi^4/4$ with a non-minimal coupling $\xi$ between the inflaton field $\phi$ and the Ricci scalar $R$ gives rise to a detectable level of gravitational waves for $|\xi|$ smaller than the order of 0.01, irrespective of the reheating temperature.

Understanding the formation history of brightest cluster galaxies is an important topic in galaxy formation. Utilizing the Planck Sunyaev-Zel'dovich cluster catalog, and applying the Ansatz that the most massive halos at one redshift remain among the most massive ones at a slightly later cosmic epoch, we have constructed cluster samples at redshift z~0.4 and z~0.2 that can be statistically regarded as progenitor-descendant pairs. This allows us to study the stellar mass assembly history of BCGs in these massive clusters at late times, finding the degree of growth between the two epochs is likely at only few percent level, which is far lower compared to the prediction from a state-of-the-art semi-analytic galaxy formation model.

We explore the observability of the neutral hydrogen (HI) and the singly-ionized isotope helium-3 ($^3$HeII) in the intergalactic medium (IGM) from the Epoch of Reionization down to the local Universe. The hyperfine transition of $^3$HeII, which is not as well known as the HI transition, has energy splitting corresponding to 8 cm. It also has a larger spontaneous decay rate than that of neutral hydrogen, whereas its primordial abundance is much smaller. Although both species are mostly ionized in the IGM, the balance between ionization and recombination in moderately high density regions renders them abundant enough to be observed. We estimate the emission signal of both hyperfine transitions from large scale filamentary structures and discuss the prospects for observing them with current and future radio telescopes. We conclude that HI in filaments is possibly observable even with current telescopes after 100 hours of observation. On the other hand, $^3$HeII is only detectable with future telescopes, such as SKA, after the same amount of time.

Decaying dark matter (DDM) is a candidate which can solve the discrepancies between predictions of the concordance $\Lambda$CDM model and observations at small scales such as the number counts of companion galaxies of the Milky Way and the density profile at the center of galaxies. Previous studies are limited to the cases where the decay particles are massless and/or have almost degenerate masses with that of mother particles. Here we expand the DDM models so that one can consider the DDM with arbitrary lifetime and the decay products with arbitrary masses. We calculate the time evolutions of perturbed phase-space distribution functions of decay products for the first time and study effects of DDM on the temperature anisotropy in the cosmic microwave background and the matter power spectrum at present. From a recent observational estimate of $\sigma_{8}$, we derive constraints on the lifetime of DDM and the mass ratio between the decay products and DDM. We also discuss implications of the DDM model for the discrepancy in the measurements of $\sigma_8$ recently claimed by the Planck satellite collaboration.

We investigate the relationships between stellar mass, gas-phase oxygen abundance (metallicity), star formation rate, and dust content of star-forming galaxies at z$\sim$1.6 using Subaru/FMOS spectroscopy in the COSMOS field. The mass-metallicity relation at $z\sim1.6$ is steeper than the relation observed in the local Universe. The steeper MZ relation at $z\sim1.6$ is mainly due to evolution in the stellar mass where the MZ relation begins to turnover and flatten. This turnover mass is 1.2 dex larger at $z\sim1.6$. The most massive galaxies at $z\sim1.6$ ($\sim 10^{11}M_\odot$) are enriched to the level observed in massive galaxies in the local Universe. The mass-metallicity relation we measure at $z\sim1.6$ supports the suggestion of an empirical upper metallicity limit that does not significantly evolve with redshift. We find an anti-correlation between metallicity and star formation rate for galaxies at a fixed stellar mass at $z\sim1.6$ which is similar to trends observed in the local Universe. We do not find a relation between stellar mass, metallicity and star formation rate that is independent of redshift; our data suggest that there is redshift evolution in this relation. We examine the relation between stellar mass, metallicity and dust extinction. We find that at a fixed stellar mass dustier galaxies tend to be more metal rich. From examination of the stellar masses, metallicities, SFRs and dust extinctions we conclude that stellar mass is most closely related to dust extinction.

We explore the Lagrangian perturbation theory (LPT) at 1-loop order with Gaussian initial conditions. We present an expansion method to approximately compute the power spectrum in LPT. Our approximate solution has good convergence in the series expansion and enables us to compute the power spectrum in LPT accurately and quickly. Non-linear corrections in the Lagrangian perturbation theory naturally satisfy the law of conservation of mass because the relation between matter density and the displacement vector of dark matter corresponds to the conservation of mass. By matching the 1-loop solution in LPT to the 2-loop solution in standard perturbation theory, we present an approximate solution of the power spectrum which has higher order corrections than the 2-loop order in standard perturbation theory with the conservation of mass satisfied. With this approximation, we can use LPT to compute a non-linear power spectrum without any free parameters, and this solution agrees with numerical simulations at $k=0.2$ $h{\rm Mpc}^{-1}$ and $z=0.35$ to better than 2%.

We present the first results from a near-IR spectroscopic survey of the COSMOS field, using the Fiber Multi-Object Spectrograph on the Subaru telescope, designed to characterize the star-forming galaxy population at $1.4<z<1.7$. The high-resolution mode is implemented to detect H$\alpha$ in emission between $1.6{\rm -}1.8 \mathrm{\mu m}$ with $f_{\rm H\alpha}\gtrsim4\times10^{-17}$ erg cm$^{-2}$ s$^{-1}$. Here, we specifically focus on 271 sBzK-selected galaxies that yield a H$\alpha$ detection thus providing a redshift and emission line luminosity to establish the relation between star formation rate and stellar mass. With further $J$-band spectroscopy for 89 of these, the level of dust extinction is assessed by measuring the Balmer decrement using co-added spectra. We find that the extinction ($0.6\lesssim A_\mathrm{H\alpha} \lesssim 2.5$) rises with stellar mass and is elevated at high masses compared to low-redshift galaxies. Using this subset of the spectroscopic sample, we further find that the differential extinction between stellar and nebular emission \hbox{$E_\mathrm{star}(B-V)/E_\mathrm{neb}(B-V)$} is 0.7--0.8, dissimilar to that typically seen at low redshift. After correcting for extinction, we derive an H$\alpha$-based main sequence with a slope ($0.81\pm0.04$) and normalization similar to previous studies at these redshifts.

We study the non-linear behavior of the baryon acoustic oscillation in the power spectrum and the correlation function by decomposing the dark matter perturbations into the short- and long-wavelength modes. The evolution of the dark matter fluctuations can be described as a global coordinate transformation caused by the long-wavelength displacement vector acting on short-wavelength matter perturbation undergoing non-linear growth. Using this feature, we investigate the well known cancellation of the high-$k$ solutions in the standard perturbation theory. While the standard perturbation theory naturally satisfies the cancellation of the high-$k$ solutions, some of the recently proposed improved perturbation theories do not guarantee the cancellation. We show that this cancellation clarifies the success of the standard perturbation theory at the 2-loop order in describing the amplitude of the non-linear power spectrum even at high-$k$ regions.We propose an extension of the standard 2-loop level perturbation theory model of the non-linear power spectrum that more accurately models the non-linear evolution of the baryon acoustic oscillation than the standard perturbation theory. The model consists of simple and intuitive parts: the non-linear evolution of the smoothed power spectrum without the baryon acoustic oscillations and the non-linear evolution of the baryon acoustic oscillations due to the large-scale velocity of dark matter and due to the gravitational attraction between dark matter particles. Our extended model predicts the smoothing parameter of the baryon acoustic oscillation peak at $z=0.35$ as $\sim 7.7\ {\rm Mpc}/h$ and describes the small non-linear shift in the peak position due to the galaxy random motions.

We investigate cosmological signatures of uncorrelated isocurvature perturbations whose power spectrum is blue-tilted with spectral index 2<p<4. Such an isocurvature power spectrum can promote early formation of small-scale structure, notably dark matter halos and galaxies, and may thereby resolve the shortage of ionizing photons suggested by observations of galaxies at high redshifts (z~7-8) but that are required to reionize the universe at z~10. We mainly focus on how the formation of dark matter halos can be modified, and explore the connection between the spectral shape of CMB anisotropies and the reionization optical depth as a powerful probe of a highly blue-tilted isocurvature primordial power spectrum. We also study the consequences for 21cm line fluctuations due to neutral hydrogens in minihalos. Combination of measurements of the reionization optical depth and 21cm line fluctuations will provide complementary probes of a highly blue-tilted isocurvature power spectrum.

The 21 cm signal produced by non-evaporating primordial black holes (PBHs) is investigated. X-ray photons emitted by accretion of matter onto a PBH ionize and heat the intergalactic medium (IGM) gas near the PBH. Using a simple analytic model, we show that this X-ray heating can produce an observable differential 21 cm brightness temperature. The region of the observable 21 cm brightness temperature can extend to 1-10 Mpc comoving distance from a PBH depending on the PBH mass. The angular power spectrum of 21 cm fluctuations due to PBHs is also calculated. The peak position of the angular spectrum depends on PBH mass, while the amplitude is independent of PBH mass. Comparing this angular power spectrum with the angular power spectrum caused by primordial density fluctuations, it is found that both of them become comparable if $\Omega_{{\rm PBH}} = 10^{-11} (M/10^{3} M_\odot)^{-0.2}$ at $z=30$ and $10^{-12} (M/10^{3} M_\odot)^{-0.2}$ at $z=20$ for the PBH mass from $10 M_\odot $ to $10^8 M_\odot $. Finally we find that the Square Kilometer Array can detect the signal due to PBHs up to $\Omega_{\rm PBH}=10^{-5} (M/10^{3} M_\odot)^{-0.2}$ at $z=30$ and $10^{-7} (M/10^{3} M_\odot)^{-0.2}$ at $z=20$ for PBHs with mass from $10^2 M_\odot $ to $10^8 M_\odot $.

An optimal method to constrain the non-linearity parameter g_NL of the local-type non-Gaussianity from CMB data is proposed. Our optimal estimator for g_NL is separable and can be efficiently computed in real space. Combining the exact filtering of CMB maps with the full covariance matrix, our method allows us to extract cosmological information from observed data as much as possible and obtain a tighter constraint on g_NL than previous studies. Applying our method to the WMAP 9-year data, we obtain the constraint g_NL = (-3.3 +- 2.2) 10^5, which is a few times tighter than previous ones. We also make a forecast for PLANCK data by using the Fisher matrix analysis.

Recently the lower bounds of the intergalactic magnetic fields $10^{-16} \sim 10^{-20}$ Gauss are set by gamma-ray observations while it is unlikely to generate such large scale magnetic fields through astrophysical processes. It is known that large scale magnetic fields could be generated if there exist cosmological vector mode perturbations in the primordial plasma. The vector mode, however, has only a decaying solution in General Relativity if the plasma consists of perfect fluids. In order to investigate a possible mechanism of magnetogenesis in the primordial plasma, here we consider cosmological perturbations in the Einstein-Aether gravity model, in which the aether field can act as a new source of vector metric perturbations and thus of magnetic fields. We estimate the angular power spectra of temperature and B-mode polarization of the Cosmic Microwave Background (CMB) Anisotropies in this model and put a rough constraint on the aether field parameters from latest observations. We then estimate the power spectrum of associated magnetic fields around the recombination epoch within this limit. It is found that the spectrum has a characteristic peak at $k=0.1 h{\rm Mpc^{-1}}$, and at that scale the amplitude can be as large as $B\sim 10^{-22}$ Gauss where the upper bound comes from CMB temperature anisotropies. The magnetic fields with this amplitude can be seeds of large scale magnetic fields observed today if the sufficient dynamo mechanism takes place. Analytic interpretation for the power spectra is also given.

We develop a novel abundance matching method to construct a mock catalog of luminous red galaxies (LRGs) in SDSS, using catalogs of halos and subhalos in N-body simulations for a LCDM model. Motivated by observations suggesting that LRGs are passively-evolving, massive early-type galaxies with a typical age >5Gyr, we assume that simulated halos at z=2 (z2-halo) are progenitors for LRG-host subhalos observed today, and we label the most tightly bound particles in each progenitor z2-halo as LRG ``stars''. We then identify the subhalos containing these stars to z=0.3 (SDSS redshift) in descending order of the masses of z2-halos until the comoving number density of the matched subhalos becomes comparable to the measured number density of SDSS LRGs, n=10^{-4} (h/Mpc)^3. Once the above prescription is determined, our only free parameter is the number density of halos identified at z=2 and this parameter is fixed to match the observed number density at z = 0.3. By tracing subsequent merging and assembly histories of each progenitor z2-halo, we can directly compute, from the mock catalog, the distributions of central and satellite LRGs and their internal motions in each host halo at z=0.3. While the SDSS LRGs are galaxies selected by the magnitude and color cuts from the SDSS images and are not necessarily a stellar-mass-selected sample, our mock catalog reproduces a host of SDSS measurements: the halo occupation distribution for central and satellite LRGs, the projected auto-correlation function of LRGs, the cross-correlation of LRGs with shapes of background galaxies (LRG-galaxy weak lensing), and the nonlinear redshift-space distortion effect, the Finger-of-God effect, in the angle-averaged redshift-space power spectrum.

We investigate the relation between the regularized multi-propagator method, called "Reg PT", and the standard perturbation theory. Reg PT is one of the most successful models to describe nonlinear evolution of dark matter fluctuations. However, Reg PT is a mathematically unproven interpolation formula between the large-scale solution calculated by the standard perturbation theory and the limiting solution in the small scale calculated by the multi-point propagator method. In this paper, we give an alternative explanation for Reg PT in the context of the standard perturbation theory, showing that Reg PT does not ever have more effective information on nonlinear matter evolution than the standard perturbation theory. In other words, the solutions of the standard perturbation theory reproduce the results of $N$-body simulations better than those of Reg PT, especially at the high-$k$ region. This fact means that the standard perturbation theory at the two-loop level is still one of the best predictions of the nonlinear power spectrum to date. Nevertheless, the standard perturbation theory has not been preferred because of the divergent behavior of the solution at small scales. To solve this problem, we also propose a modified standard perturbation theory which avoids the divergence.

Precise understanding of nonlinear evolution of cosmological perturbations during inflation is necessary for the correct interpretation of measurements of non-Gaussian correlations in the cosmic microwave background and the large-scale structure of the universe. The "{\delta}N formalism" is a popular and powerful technique for computing non-linear evolution of cosmological perturbations on large scales. In particular, it enables us to compute the curvature perturbation, {\zeta}, on large scales without actually solving perturbed field equations. However, people often wonder why this is the case. In order for this approach to be valid, the perturbed Hamiltonian constraint and matter-field equations on large scales must, with a suitable choice of coordinates, take on the same forms as the corresponding unperturbed equations. We find that this is possible when (1) the unperturbed metric is given by a homogeneous and isotropic Friedmann-Lema\^itre-Robertson-Walker metric; and (2) on large scales and with a suitable choice of coordinates, one can ignore the shift vector (g0i) as well as time-dependence of tensor perturbations to gij/a2(t) of the perturbed metric. While the first condition has to be assumed a priori, the second condition can be met when (3) the anisotropic stress becomes negligible on large scales. However, in order to explicitly show that the second condition follows from the third condition, one has to use gravitational field equations, and thus this statement may depend on the details of theory of gravitation. Finally, as the {\delta}N formalism uses only the Hamiltonian constraint and matter-field equations, it does not a priori respect the momentum constraint. We show that the violation of the momentum constraint only yields a decaying mode solution for {\zeta}, and the violation vanishes when the slow-roll conditions are satisfied.

We propose a new approximated expression for non-linear Dark Matter power spectrum much beyond BAO scales. The proposed expression agrees with the result of N-body simulation with the accuracy better than 2 % up to k=1.0 [h/Mpc] and k=0.7 [h/Mpc] at z=3.0 and z=1.0, respectively. Even at z=0.35, the accuracy remains within 10 % up to k=0.8 [h/Mpc]. In doing so, we used an approximation for the kernel functions used in the Standard Perturbation Theory (SPT) which is also used to prove the Reg PT proposed by Bernardeau et al. (2011).

We apply the Wiener Hermite (WH) expansion to the non-linear evolution of Large-Scale Structure, and obtain an approximate expression for the matter power spectrum in full order of the expansion. This method allows us to expand any random function in terms of an orthonormal bases in space of random functions in such a way that the first order of the expansion expresses the Gaussian distribution, and others are the deviation from the Gaussianity. It is proved that the Wiener Hermite expansion is mathematically equivalent with the $\Gamma$-expansion approach in the Renormalized Perturbation Theory (RPT). While an exponential behavior in the high-$k$ limit has been proved for the mass density and velocity fluctuations of Dark Matter in the RPT, we reprove the behavior in the context of the Wiener Hermite expansion using the result of Standard Perturbation Theory (SPT). We propose a new approximate expression for the matter power spectrum which interpolates the low-$k$ expression corresponding to 1-loop level in SPT and the high-$k$ expression obtained by taking a high-$k$ limit of the WH expansion. The validity of our prescription is specifically verified by comparing with the 2-loop solutions of the SPT. The proposed power spectrum agrees with the result of $N$-body simulation with the accuracy better than 1% or 2% in a range of the BAO scales, where the wave number is about $k= 0.2 \sim 0.4$ $h{\rm Mpc^{-1}}$ at $z=0.5-3.0$. This accuracy is comparable to or slightly less than the ones in the Closure Theory whose the fractional difference from N body result is within 1%. A merit of our method is that the computational time is very short because only single and double integrals are involved in our solution.

The angular power spectrum of the cosmic microwave background (CMB) temperature anisotropies is a good probe to look into the primordial density fluctuations at large scales in the universe. Here we re-examine the angular power spectrum of the Wilkinson Microwave Anisotropy Probe data, paying particular attention to the fine structures (oscillations) at $\ell=100 \sim 150$ reported by several authors. Using Monte-Carlo simulations, we confirm that the gap from the simple power law spectrum is a rare event, about 2.5--3$\sigma$, if these fine structures are generated by experimental noise and the cosmic variance. Next, in order to investigate the origin of the structures, we examine frequency and direction dependencies of the fine structures by dividing the observed QUV frequency maps into four sky regions. We find that the structures around $\ell \sim 120$ do not have significant dependences either on frequencies or directions. For the structure around $\ell \sim 140$, however, we find that the characteristic signature found in the all sky power spectrum is attributed to the anomaly only in the South East region.

Shadows of multi-black holes have structures distinct from the mere superposition of the shadow of a single black hole: the eyebrow-like structures outside the main shadows and the deformation of the shadows. We present analytic estimates of these structures using the static multi-black hole solution (Majumdar-Papapetrou solution). We show that the width of the eyebrow is related with the distance between the black holes and that the shadows are deformed into ellipses due to the presence of the second black holes. These results are helpful to understand qualitatively the features of the shadows of colliding black holes. We also present the shadows of colliding/coalescing black holes in the Kastor-Traschen solution.

The consistency relation between non-linear parameters $f_{NL}$ and $\tau_{NL}$ characterizing Non-Gaussianity generated during the inflationary period have been emerged as a useful tool which have a possibility to rule out a large class of inflationary models all at once. In our previous work, we extended the Suyama-Yamaguchi inequality up to 1-loop corrections. In this paper, we derive an inequality extended up to {\it all} loop corrections which has the same expression with the original Suyama-Yamaguchi inequality, $\tau_{NL} \geq (6/5f_{NL})^2$, where the equality is not satisfied in the case of single field models any more.

The stochastic gravitational wave background (GWB) from halo mergers is investigated by a quasi-analytic method. The method we employ consists of two steps. The first step is to construct a merger tree by using the Extended Press-Schechter formalism or the Sheth & Tormen formalism, with Monte-Carlo realizations. This merger tree provides evolution of halo masses. From $N$-body simulation of two-halo mergers, we can estimate the amount of gravitational wave emission induced by the individual merger process. Therefore the second step is to combine this gravitaional wave emission to the merger tree and obtain the amplitude of GWB. We find $\Omega_{GW}\sim 10^{-19}$ for $f\sim 10^{-17}-10^{-16}$ Hz, where $\Omega_{GW}$ is the energy density of the GWB. It turns out that most of the contribution on the GWB comes from halos with masses below $10^{15} M_\odot$ and mergers at low redshift, i.e., $0<z<0.8$.

We study the lensing convergence power spectrum and its covariance for a standard LCDM cosmology. We run 400 cosmological N-body simulations and use the outputs to perform a total of 1000 independent ray-tracing simulations. We compare the simulation results with analytic model predictions. The semi-analytic model based on Smith et al.(2003) fitting formula underestimates the convergence power by ~30% at arc-minute angular scales. For the convergence power spectrum covariance, the halo model reproduces the simulation results remarkably well over a wide range of angular scales and source redshifts. The dominant contribution at small angular scales comes from the sample variance due to the number fluctuations of halos in a finite survey volume. The signal-to-noise ratio for the convergence power spectrum is degraded by the non-Gaussian covariances by up to a factor 5 for a weak lensing survey to z_s ~1. The probability distribution of the convergence power spectrum estimators, among the realizations, is well approximated by a chi-square distribution with broadened variance given by the non-Gaussian covariance, but has a larger positive tail. The skewness and kurtosis have non-negligible values especially for a shallow survey. We argue that a prior knowledge on the full distribution may be needed to obtain an unbiased estimate on the ensemble averaged band power at each angular scale from a finite volume survey.

A primordial vector mode and its associated magnetic field generation are investigated. Firstly, we put an observational constraint on the amount of the primordial vector mode from the seven-year WMAP data. The constraint is found as $r_v \lesssim -\frac{r}{40}+0.012$, where $r_v$ and $r$ are the amounts of vector and tensor perturbation amplitudes with respect to the scalar one, respectively. Secondly, we calculate the spectrum of magnetic fields inevitably created from the primordial vector mode, given the constraint on $r_v$. It is found that the maximum amount of magnetic fields generated from the vector mode is given by $B\lesssim 10^{-22}{\rm G} (\frac{r_v}{0.012})^{1/2} (\frac{k}{0.002})^{(n_v+1)/2}$ with $n_v$ being a spectral index of the vector mode. We find a non-trivial cancellation of the magnetic field generation in the radiation dominated era, which creates a characteristic cut off in the magnetic field spectrum around $k\approx 1.0$ Mpc$^{-1}$.

We compute the shadows of colliding black holes using the Kastor-Traschen cosmological multiblack hole solution that is an exact solution describing the collision of maximally charged black holes with a positive cosmological constant. We find that in addition to the shadow of each black hole, an eyebrowlike structure appears as the black holes come close to each other. These features can be used as probes to find the multiblack hole system at the final stage of its merger process.

Decaying dark matter cosmological models have been proposed to remedy the overproduction problem at small scales in the standard cold dark matter paradigm. We consider a decaying dark matter model in which one CDM mother particle decays into two daughter particles, with arbitrary masses. A complete set of Boltzmann equations of dark matter particles is derived which is necessary to calculate the evolutions of their energy densities and their density perturbations. By comparing the expansion history of the universe in this model and the free-streaming scale of daughter particles with astronomical observational data, we give constraints on the lifetime of the mother particle, $\Gamma^{-1}$, and the mass ratio between the daughter and the mother particles $m_{\rm D}/m_{\rm M}$. From the distance to the last scattering surface of the cosmic microwave background, we obtain $\Gamma^{-1}>$ 30 Gyr in the massless limit of daughter particles and, on the other hand, we obtain $m_{\rm D} >$ 0.97$m_{\rm M}$ in the limit $\Gamma^{-1}\to 0$. The free-streaming constraint tightens the bound on the mass ratio as $(\Gamma^{-1}/10^{-2}{\rm Gyr}) \lesssim ((1-m_{\rm D1}/m_{\rm M})/10^{-2})^{-3/2}$ for $\Gamma^{-1} < H^{-1}(z=3)$.

We report the first observational search for tau neutrinos from gamma ray bursts (GRBs) using one of the Ashra light collectors. The earth-skimming tau-neutrino technique of imaging Cherenkov tau showers was applied as a detection method. We set stringent upper limits on the tau-neutrino fluence in PeV-EeV region for 3780 s (between 2.83 and 1.78 hours before) and another 3780 s (between 21.2 and 22.2 hours after) surrounding GRB081203A triggered by the Swift satellite. This first search for PeV-EeV tau neutrino complements other experiments in energy range and methodology, and suggests the prologue of "multi-particle astronomy" with a precise determination of time and location.

We study the halo mass function in the presence of the kurtosis type of primordial non-Gaussianity. The kurtosis corresponds to the trispectrum as defined in Fourier space. The primordial trispectrum is commonly characterized by two parameters, $\tau_{\rm NL}$ and $g_{\rm NL}$. As applications of the derived non-Gaussian mass function, we consider the effect on the abundance of void structure, the effect on early star formation and on formation of the most massive object at high redshift. We show that by comparing the effects of primordial non-Gaussianity on cluster abundance with that on void abundance, we can distinguish between the skewness and the kurtosis types of primordial non-Gaussianity. As for early star formation, we show that the kurtosis type of primordial non-Gaussianity seems not to affect the reionization history of the Universe on average. However, at high redshifts (up to $z\simeq 20$) such non-Gaussianity does somewhat affect the early stages of reionization.

While detection of the "local form" bispectrum of primordial perturbations would rule out all single-field inflation models, multi-field models would still be allowed. We show that multi-field models described by the $\delta N$ formalism obey an inequality between $f_{\rm NL}$ and one of the local-form {\it trispectrum} amplitudes, $\tau_{\rm NL}$, such that $\tau_{\rm NL}>\frac12(\frac65f_{\rm NL})^2$ with a possible logarithmic scale dependence, provided that 2-loop terms are small. Detection of a violation of this inequality would rule out most of multi-field models, challenging inflation as a mechanism for generating the primoridal perturbations.

A possible origin of the anomalous dip and bump in the primordial power spectrum, which are reconstructed from WMAP data corresponding to the multipole $\ell=100\sim 140$ by using the inversion method, is investigated as a consequence of modification of scalar field dynamics in the inflation era. Utilizing an analytic formula to handle higher order corrections to the slow-roll approximation, we evaluate the relation between a detailed shape of inflaton potential and a fine structure in the primordial power spectrum. We conclude that it is unlikely to generate the observed dip and bump in the power spectrum by adding any features in the inflaton potential. Though we can make a fine enough shape in the power spectrum by controlling the feature of the potential, the amplitude of the dip and bump becomes too small in that case.

The number count of ionized bubbles on a map of 21 cm fluctuations with the primordial non-Gaussianity is investigated. The existence of the primordial non-Gaussianity modifies the reionization process, because the formation of collapsed objects, which could be the source of reionization photons, is affected by the primordial non-Gaussianity. In this paper, the abundance of ionized bubbles is calculated by using a simple analytic model with the local type of the primordial non-Gaussianity, which is parameterized by f_{NL}. In order to take into account the dependence of the number count on the size of ionized bubbles and the resolution of the observation instrument, a threshold parameter B_b which is related to the the surface brightness temperature contrast of an ionized bubble is introduced. We show the potential to put the constraint on f_{NL} from the number count by future observations such as LOFAR and SKA.

We study the potential impact of detecting the inflationary gravitational wave background by the future space-based gravitational wave detectors, such as DECIGO and BBO. The signal-to-noise ratio of each experiment is calculated for chaotic/natural/hybrid inflation models by using the precise predictions of the gravitational wave spectrum based on numerical calculations. We investigate the dependence of each inflation model on the reheating temperature which influences the amplitude and shape of the spectrum, and find that the gravitational waves could be detected for chaotic/natural inflation models with high reheating temperature. From the detection of the gravitational waves, a lower bound on the reheating temperature could be obtained. The implications of this lower bound on the reheating temperature for particle physics are also discussed.

Primordial magnetic fields generated in the very early universe are one of the candidates for the origin of magnetic fields observed in galaxy clusters. After recombination, the Lorentz force acts on the residual ions and electrons to generate density fluctuations of baryons. Accordingly these fluctuations induce the early formation of dark halos which cause the Sunyaev-Zel'dovich (S-Z) effect in cosmic microwave background radiation. This additional S-Z effect due to primordial magnetic fields amplifies the angular power spectrum of cosmic microwave temperature anisotropies on small scales. This amplification depends on the comoving amplitude and the power law index of the primordial magnetic fields spectrum. Comparing with the small scale CMB observations, we obtained the constraints on the primordial magnetic fields, i.e., B < 2.0 nGauss for n=-2.9 or B < 1.0 nGauss for n=-2.6, where B is the comoving amplitude of magnetic fields at h^-1 Mpc and n is the power law index. Future S-Z measurements have the potential to give constraints tighter than those from temperature anisotropies and polarization of cosmic microwave background induced by the magnetic fields at the recombination epoch.

We study the sample variance of the matter power spectrum for the standard Lambda Cold Dark Matter universe. We use a total of 5000 cosmological N-body simulations to study in detail the distribution of best-fit cosmological parameters and the baryon acoustic peak positions. The obtained distribution is compared with the results from the Fisher matrix analysis with and without including non-Gaussian errors. For the Fisher matrix analysis, we compute the derivatives of the matter power spectrum with respect to cosmological parameters using directly full nonlinear simulations. We show that the non-Gaussian errors increase the unmarginalized errors by up to a factor 5 for k_{max}=0.4h/Mpc if there is only one free parameter provided other parameters are well determined by external information. On the other hand, for multi-parameter fitting, the impact of the non-Gaussian errors is significantly mitigated due to severe parameter degeneracies in the power spectrum. The distribution of the acoustic peak positions is well described by a Gaussian distribution, with its width being consistent with the statistical interval predicted from the Fisher matrix. We also examine systematic bias in the best-fit parameter due to the non-Gaussian errors. The bias is found to be smaller than the 1 sigma statistical error for both the cosmological parameters and the acoustic scale positions.

A systematic study of gravitational waves from galaxy mergers, through N-body simulations, was performed. In particular, we investigated the relative importance of galaxy components (disk, bulge and halo) and effects of initial relative velocity, relative angular momentum and mass ratio of the galaxies. We found that the features of light curve of gravitational waves, such as peak width and luminosity, are reliably simulated with particle numbers larger than ~10^4. Dominant contribution to gravitational wave emission came from the halo component, while peak luminosity amounted to 10^31 erg/sec for the collision of two halos with mass 3.8 x10^12Msun/h. We also found that the initial relative velocity in the direction of the initial separation did not significantly affect gravitational wave emission, while the initial relative angular momentum broadened the peak width and suppressed the luminosity. Mass dependence of the peak luminosity was also investigated, and we obtained evidence that the luminosity is proportional to the cubic mass when the scaling relation is satisfied. This behavior was considered by a simple analysis.

We present a clustering analysis of ~60,000 massive (stellar mass Mstar > 10^{11} Msun) galaxies out to z = 1 drawn from 55.2 deg2 of the UKIRT Infrared Deep Sky Survey (UKIDSS) and the Sloan Digital Sky Survey (SDSS) II Supernova Survey. Strong clustering is detected for all the subsamples of massive galaxies characterized by different stellar masses (Mstar = 10^{11.0-11.5} Msun, 10^{11.5-12.0} Msun) or rest-frame colors (blue: U - V < 1.0, red: U - V > 1.0). We find that more mature (more massive or redder) galaxies are more clustered, which implies that more mature galaxies have started stellar-mass assembly earlier within the highly-biased region where the structure formation has also started earlier. By means of halo occupation distribution (HOD) models fitted to the observed angular correlation function, we infer the properties of the underlying host dark halos. We find that the estimated bias factors and host halo masses are systematically larger for galaxies with larger stellar masses, which is consistent with the general agreement that the capability of hosting massive galaxies depends strongly on halo mass. The estimated effective halo masses are ~10^{14} Msun, which gives the stellar-mass to halo-mass ratios of ~0.003. The observed evolution of bias factors indicates rapid evolution of spatial distributions of cold dark matter relative to those traced by the massive galaxies, while the transition of host halo masses might imply that the fractional mass growth rate of halos is less than those of stellar systems. The inferred halo masses and high fractions of central galaxies indicate that the massive galaxies in the current sample are possibly equivalent to central galaxies of galaxy clusters.

We study how direct detection of the inflationary gravitational wave background constrains inflationary parameters and complements CMB polarization measurements. The error ellipsoids calculated using the Fisher information matrix approach with Planck and the direct detection experiment, BBO (Big Bang Observer), show different directions of parameter degeneracy, and the degeneracy is broken when they are combined. For a slow-roll parameterization, we show that BBO could significantly improve the constraints on the tensor-to-scalar ratio compared with Planck alone. We also look at a quadratic and a natural inflation model. In both cases, if the temperature of reheating is also treated as a free parameter, then the addition of BBO can significantly improve the error bars. In the case of natural inflation, we find that the addition of BBO could even partially improve the error bars of a cosmic variance-limited CMB experiment.

We revisit ghost dark matter, the possibility that ghost condensation may serve as an alternative to dark matter. In particular, we investigate the Friedmann-Robertson-Walker (FRW) background evolution and the large-scale structure (LSS) in the $\Lambda$GDM universe, i.e. a late-time universe dominated by a cosmological constant and ghost dark matter. The FRW background of the $\Lambda$GDM universe is indistinguishable from that of the standard $\Lambda$CDM universe if $M\gtrsim 1 {\rm eV}$, where $M$ is the scale of spontaneous Lorentz breaking. From the LSS we find a stronger bound: $M\gtrsim 10 {\rm eV}$. For smaller $M$, ghost dark matter would have non-negligible sound speed after the matter-radiation equality, and thus the matter power spectrum would significantly differ from observation. These bounds are compatible with the phenomenological upper bound $M\lesssim 100 {\rm GeV}$ known in the literature.

We present results of an imaging observation of the central region of a giant radio galaxy B1358+305. The classical, standard scenario of Fanaroff-Riley II radio galaxies suggests that shock produced hot electrons contained in a radio galaxy are a good reservoir of the jet-supplied energy from active nuclei. The aim of our observation is to search for the Sunyaev-Zel'dovich effect induced by these hot electrons. The observation was performed at 21 GHz with the Nobeyama 45-m telescope. Deep imaging observation of a wide region of size 6.7'x6.7' with the beam size theta_HPBW=81.2" enables the most detailed examination of the possible thermal energy of electrons contained in a radio galaxy. The resultant intensity fluctuation is 0.56 mJy/beam (in terms of the Compton y-parameter, y=1.04x 10^{-4}) at a 95 percent confidence level. The intensity fluctuation obtained with imaging analysis sets the most stringent upper limit on the fluctuations in the central region of a giant radio galaxy obtained so far, and our results will be a toehold for future plans of SZE observation in a radio galaxy.

We numerically study the interaction of inflation-produced magnetic fields with gravitational waves, both of which originate from quantum fluctuations during inflation. The resonance between the magnetic field perturbations and the gravitational waves has been suggested as a possible mechanism for magnetic field amplification. However, some analytical studies suggest that the effect of the inflationary gravitational waves is too small to provide significant amplification. Our numerical study shows more clearly how the interaction affects the magnetic fields and confirms the weakness of the influence of the gravitational waves. We present an investigation based on the magnetohydrodynamic approximation and take into account the differences of the Alfven speed.

In this paper, we put constraints on neutrino properties such as mass $m_{\nu}$ and degeneracy parameters $\xi_i$ from WMAP5 data and light element abundances by using a Markov chain Monte Carlo (MCMC) approach. In order to take consistently into account the effects of the degeneracy parameters, we run the Big Bang Nucleosynthesis code for each value of $\xi_i$ and the other cosmological parameters to estimate the Helium abundance, which is then used to calculate CMB anisotropy spectra instead of treating it as a free parameter. We find that the constraint on $m_{\nu}$ is fairly robust and does not vary very much even if the lepton asymmetry is allowed, and is given by $\sum m_\nu < 1.3 \rm eV$ ($95 % \rm C.L.$).

The spectrum of the gravitational wave background originating from quantum fluctuations during inflation is calculated numerically for various inflation models over a wide range of frequencies. We take into account four ingredients : the scalar field dynamics during inflation making no use of the slow-roll approximation, the fermionic decay of the scalar field with a small coupling constant during the reheating process, the change of the effective number of degrees of freedom g_* in the radiation-dominated era, and the anisotropic stress of free-streaming neutrinos. By numerically solving the evolution of gravitational waves during and after inflation up to the present, all of these effects can be examined comprehensively and accurately over a broad spectrum, even at very high frequencies. We find that the spectrum shows (i) a large deviation from the spectrum less accurate obtained by Taylor expanding around the CMB scale using the slow-roll approximation (ii) a characteristic frequency dependence due to the reheating effect, and (iii) damping due to the g_* changes and the neutrino anisotropic stress. We suggest that future analysis of the gravitational wave background should take into consideration the fact that analytical estimates using the Taylor expansion overestimate the amplitude of the spectrum.

We use 5000 cosmological N-body simulations of 1(Gpc/h)^3 box for the concordance LCDM model in order to study the sampling variances of nonlinear matter power spectrum. We show that the non-Gaussian errors can be important even on large length scales relevant for baryon acoustic oscillations (BAO). Our findings are (1) the non-Gaussian errors degrade the cumulative signal-to-noise ratios (S/N) for the power spectrum amplitude by up to a factor of 2 and 4 for redshifts z=1 and 0, respectively. (2) There is little information on the power spectrum amplitudes in the quasi-nonlinear regime, confirming the previous results. (3) The distribution of power spectrum estimators at BAO scales, among the realizations, is well approximated by a Gaussian distribution with variance that is given by the diagonal covariance component. (4) For the redshift-space power spectrum, the degradation in S/N by non-Gaussian errors is mitigated due to nonlinear redshift distortions. (5) For an actual galaxy survey, the additional shot noise contamination compromises the cosmological information inherent in the galaxy power spectrum, but also mitigates the impact of non-Gaussian errors. The S/N is degraded by up to 30% for a WFMOS-type survey. (6) The finite survey volume causes additional non-Gaussian errors via the correlations of long-wavelength fluctuations with the fluctuations we want to measure, further degrading the S/N values by about 30% even at high redshift z=3.

How did the universe evolve? The fine angular scale (l>1000) temperature and polarization anisotropies in the CMB are a Rosetta stone for understanding the evolution of the universe. Through detailed measurements one may address everything from the physics of the birth of the universe to the history of star formation and the process by which galaxies formed. One may in addition track the evolution of the dark energy and discover the net neutrino mass. We are at the dawn of a new era in which hundreds of square degrees of sky can be mapped with arcminute resolution and sensitivities measured in microKelvin. Acquiring these data requires the use of special purpose telescopes such as the Atacama Cosmology Telescope (ACT), located in Chile, and the South Pole Telescope (SPT). These new telescopes are outfitted with a new generation of custom mm-wave kilo-pixel arrays. Additional instruments are in the planning stages.

We calculate the spin temperature and 21 cm brightness of early \HII regions around the first stars. We use outputs from cosmological radiation-hydrodynamics simulations of the formation and evolution of early \HII regions. In the pre-reionization era, \HII regions around massive primordial stars have diameters of a few kpc. The gas within the \HII regions is almost fully ionized, but begins recombining after the central stars die off. The relic \HII regions are then seen as bright {\it emission} sources in hydrogen 21 cm. We make brightness temperature maps of the \HII regions, accounting for radiative coupling with Lyman-$\alpha$ photons in a simplified manner. The spin temperature in the relic \HII region is close to the gas kinetic temperature, generally several hundred to several thousand degrees. We show that the relic \HII region can be as bright as $\delta T_{\rm b} \sim 100$ mK in differential temperature against the cosmic microwave background for an angular resolution of sub-arcseconds. While individual early \HII patches will not be identified by currently planned radio telescopes, the collective fluctuations from early \HII regions might imprint signatures in the 21 cm background.

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.

Quasi-equilibrium models of rapidly rotating triaxially deformed stars are computed in general relativistic gravity, assuming a conformally flat spatial geometry (Isenberg-Wilson-Mathews formulation) and a polytropic equation of state. Highly deformed solutions are calculated on the initial slice covered by spherical coordinate grids, centered at the source, in all angular directions up to a large truncation radius. Constant rest mass sequences are calculated from nearly axisymmetric to maximally deformed triaxial configurations. Selected parameters are to model (proto-) neutron stars; the compactness is M/R = 0.001, 0.1, 0.14, 0.2 for polytropic index n = 0.3 and M/R = 0.001, 0.1, 0.12, 0.14 for n = 0.5. We confirmed that the triaxial solutions exist for these parameters as in the case of Newtonian polytropes. However, it is also found that the triaxial sequences become shorter for higher compactness, and those may disappear at a certain large compactness for the n = 0.5 case. In the scenario of the contraction of proto-neutron stars being subject to strong viscosity and rapid cooling, it is plausible that, once the viscosity driven secular instability sets in during the contraction, the proto-neutron stars are always maximally deformed triaxial configurations, as long as the compactness and the equation of state parameters allow such triaxial sequences. Detection of gravitational waves from such sources may be used as another probe for the nuclear equation of state.

The S-Z effect and Faraday rotation from halos are examined over a wide mass range, an including gas condensation and magnetic field evolution. Contributions to the CMB angular power spectrum are evaluated for galaxy clusters, galaxy groups and galaxies. Smaller mass halos are found to play a more important role than massive halos for the B-mode polarisation associated with the S-Z CMB anisotropies. The B-modes from Faraday rotation dominate the secondary B-modes caused by gravitational lensing at l > 3000. Measurement of B-mode polarisation in combination with the S-Z power spectrum can potentially provide important constraints on intracluster magnetic field and gas evolution at early epochs.

Possible influence of primordial black hole (PBH) evaporations on cosmic microwave background (CMB) is investigated. The spectrum distortions of CMB from the black-body spectrum are described by the chemical potential $\mu$ and the Compton parameter $y$. From COBE/FIRAS limits on $\mu$ and $y$, the power law index $n$ of primordial density fluctuations and the mass fraction of PBHs $\beta$ are constrained by employing the peak theory for the formation process of PBHs. Constraints set here are $n < 1.304$ and $n<1.333$ in the thresholds of peaks $\zeta_{\rm th} =0.7$ and $\zeta_{\rm th} =1.2$, respectively, for the PBH mass range between $2.7\times 10^{11}$g and $1.6 \times 10^{12}$g, and $n < 1.312$ and $n<1.343$ in the thresholds of peaks $\zeta_{\rm th} =0.7$ and $\zeta_{\rm th} =1.2$, respectively, for the PBH mass range between $1.6 \times 10^{12} {\rm g}$ and $3.5\times 10^{13}$g, which correspond to the comoving scales between $3 \times 10^{-18}$Mpc and $ 4\times 10^{-17}$Mpc. The constraint on the PBH fraction, which is the direct probe of the amplitude of density fluctuations in these scales, stays an almost same value as $\beta<10^{-21}$ in these mass ranges. It is also found that, with these constraints, UV photons injected by PBH evaporations are unlikely ionize vast hydrogen atoms after recombination.

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.

Detailed physical processes of magnetic field generation from density fluctuations in the pre-recombination era are studied. Solving Maxwell equations and the generalized Ohm's law, the evolutions of the net charge density, the electric current and the electromagnetic field are solved. Unlike most of previous works, we treat electrons and photons as separate components under the assumption of tight coupling. We find that generation of the magnetic field due to density fluctuations takes place only from the second order of both perturbation theory and the tight coupling approximation.

Non-linear evolution of density fluctuations yields secondary anisotropies in the cosmic microwave background (CMB), which are correlated with the same density fluctuations that can be measured by weak lensing (WL) surveys. We study the CMB-WL correlation using analytical models as well as $N$-body simulations. We show that an analytical model based upon the time derivative of non-linear matter power spectrum agrees with simulations. All-sky cosmic-variance limited CMB and WL surveys allow us to measure the correlation from non-linearity with high significance ($50\sigma$) for $l_{\rm max}=10^4$, whereas the forthcoming missions such as Planck and LSST are expected to yield only marginal detections. The CMB-WL correlation is sensitive to the time derivative of structure growth. We study how this property may be used to constrain the nature of dark energy. While the expected constraints are not very strong, they may provide a cross check of results from other observations.

The angular power spectra of cosmic microwave background are analyzed under the light of the evolutional scenario of the universe based on the renormalizable quantum theory of gravity in four dimensions. The equation of evolution is solved numerically fixing the power law spectrum predicted by the conformal gravity for the initial condition. The equation requires to introduce a dynamical energy scale about 10^{17}GeV, where the inflationary space-time evolution makes a transition to the big-bang of the conventional Friedmann universe. The quality of fit to the three-year data of WMAP implies the possibility to understand the observation by quantum gravity.

A generation of magnetic fields from cosmological density perturbations is investigated. In the primordial plasma before cosmological recombination, all of the materials except dark matter in the universe exist in the form of photons, electrons, and protons (and a small number of light elements). Due to the different scattering nature of photons off electrons and protons, electric currents and electric fields are inevitably induced, and thus magnetic fields are generated. We present a detailed formalism of the generation of cosmological magnetic fields based on the framework of the well-established cosmological perturbation theory following our previous works. We numerically obtain the power spectrum of magnetic fields for a wide range of scales, from $k\sim 10^{-5}$ Mpc$^{-1}$ up to $k\sim 2\times 10^{4}$ Mpc$^{-1}$ and provide its analytic interpretation. Implications of these cosmologically generated magnetic fields are discussed.

A simple analytical model is used to calculate the X-ray heating of the IGM for a range of black hole masses. This process is efficient enough to decouple the spin temperature of the intergalactic medium from the cosmic microwave background (CMB) temperature and produce a differential brightness temperature of the order of $\sim 5-20 \mathrm{mK}$ out to distances as large as a few co-moving Mpc, depending on the redshift, black hole mass and lifetime. We explore the influence of two types of black holes, those with and without ionising UV radiation. The results of the simple analytical model are compared to those of a full spherically symmetric radiative transfer code. Two simple scenarios are proposed for the formation and evolution of black hole mass density in the Universe. The first considers an intermediate mass black hole that form as an end-product of Population III stars, whereas the second considers super-massive black holes that form directly through the collapse of massive halos with low spin parameter. These scenarios are shown not to violate any of the observational constraints, yet produce enough X-ray photons to decouple the spin-temperature from that of the CMB. This is an important issue for future high redshift 21 cm observations.

Primordial magnetic fields possibly generated in the very early universe are one of the candidates for the origin of magnetic fields observed in many galaxies and galaxy clusters. After recombination, the dissipation process of the primordial magnetic fields increases the baryon temperature. The Lorentz force acts on the residual ions and electrons to generate density fluctuations. These effects are imprinted on the cosmic microwave background (CMB) brightness temperature fluctuations produced by the neutral hydrogen 21cm line. We calculate the angular power spectrum of brightness temperature fluctuations for the model with the primordial magnetic fields of a several nano Gauss strength and a power-law spectrum. It is found that the overall amplitude and the shape of the brightness temperature fluctuations depend on the strength and the spectral index of the primordial magnetic fields. Therefore, it is expected that the observations of the CMB brightness temperature fluctuations give us a strong constraint on the primordial magnetic fields.

The early reionization of the intergalactic medium, which is favored from the WMAP temperature-polarization cross-correlations, contests the validity of the standard scenario of structure formation in the cold dark matter cosmogony. It is difficult to achieve early enough star formation without rather extreme assumptions such as very high escape fraction of ionizing photons from proto-galaxies or a top-heavy initial mass function. Here we propose an alternative scenario that is additional fluctuations on small scales induced by primordial magnetic fields trigger the early structure formation. We found that ionizing photons from Population III stars formed in dark haloes can easily reionize the universe by $z \simeq 15$ if the strength of primordial magnetic fields is larger than $0.6 \times 10^{-9}$Gauss.

The origin of the substantial magnetic fields that are found in galaxies and on even larger scales, such as in clusters of galaxies, is yet unclear. If the second-order couplings between photons and electrons are considered, then cosmological density fluctuations, which explain the large scale structure of the universe, can also produce magnetic fields on cosmological scales before the epoch of recombination. By evaluating the power spectrum of these cosmological magnetic fields on a range of scales, we show here that magnetic fields of 10^{-18.1} gauss are generated at a 1 megaparsec scale and can be even stronger at smaller scales (10^{-14.1} gauss at 10 kiloparsecpc). These fields are large enough to seed magnetic fields in galaxies and may therefore have affected primordial star formation in the early universe.

We have modelled the process of reionization of the IGM by photoionization by galaxies, in order to predict the epoch of reionization. We use a sophisticated semi-analytic model to track the formation of galaxies. Our study represents a much more complete and physically consistent modelling of reionization than has been conducted in the past, containing significant improvements in the modelling of the collapse of baryons into dark matter halos, and in the model for the cooling and condensation of gas within halos (including photoheating from a self-consistently computed ionizing background and cooling due to H2). We find that reionization can be achieved by z~10-20 in a LCDM cosmological model with sigma8~0.9. However, a cosmological model with a running spectral index is only able to achieve reionization before z~9 if extreme assumptions are made about the physics of feedback at high redshifts. We also consider the galaxy formation model discussed by Baugh et al., and find that it is able to reionize the Universe by z~12. The previous results assume that all of the ionizing photons produced by stars in galaxies are able to escape and ionize the IGM. If this is not the case, then the redshift of reionization could be substantially reduced. We find that extended periods of partial reionization and double reionizations can occur in models in which the first stars formed via cooling by H2 molecules, are very massive, and in which the escape fraction of ionizing photons ~10-30%. Such models do not fully reionize until z~6-7, but predict an electron scattering optical depth as large as 0.15. Models with lower sigma8=0.7-0.8 as suggested by the recent WMAP three year data have reduced redshifts of reionization, but can be consistent with the lower optical depth also suggested by the WMAP three year data.

In this work we investigate the effects of the primordial magnetic fields on cosmic microwave background anisotropies (CMB). Based on cosmological magnetohydrodynamic (MHD) simulations we calculate the CMB anisotropy spectra and polarization induced by fluid fluctuations (Alfv\'en modes) generated by primordial magnetic fields. The strongest effect on the CMB spectra comes from the transition epoch from a turbulent regime to a viscous regime. The balance between magnetic and kinetic energy until the onset of the viscous regime provides a one to one relation between the comoving coherence length $L$ and the comoving magnetic field strength $B$, such as $L \sim 30 (B/10^{-9}{\rm G})^3 \rm pc$. The resulting CMB temperature and polarization anisotropies are somewhat different from the ones previously obtained by using linear perturbation theory. Our calculation gives a constraint on the magnetic field strength in the intermediate scale of CMB observations. Upper limits are set by WMAP and BOOMERANG results for comoving magnetic field strength of $B < 28 \rm nG$ with a comoving coherence length of $L > 0.7 \rm Mpc$ for the most extreme case, or $B < 30 \rm nG$ and $L> 0.8 \rm Mpc$ for the most conservative case.

We discuss generation of magnetic field from cosmological perturbations. We consider the evolution of three component plasma (electron, proton and photon) evaluating the collision term between elecrons and photons up to the second order. The collision term is shown to induce electric current, which then generate magnetic field. There are three contributions, two of which can be evaluated from the first-order quantities, while the other one is fluid vorticity which is purely second order. We compute numerically the magnitudes of the former contributions and shows that the amplitude of the produced magnetic field is about $\sim 10^{-19} {\rm G}$ at 10kpc comoving scale at present. Compared to astrophysical and inflationary mechanisms for seed-field generation, our study suffers from much less ambiguities concerning unknown physics and/or processes.

On 2004 December 27, a giant flare from the soft gamma repeater 1806$-$20 was observed. The radiation mechanism of the initial peak of the flare would be controversial. In this letter we point out that very high-energy cosmic rays would be produced in the case that the flare was caused by internal shocks, as is usually considered for gamma-ray bursts. The highest energy of cosmic rays can reach $10^{19}$ eV, if the Lorentz factor of the shocks is sufficiently high. Future observations of cosmic rays will inform us about the mechanism of the giant flare.

The age of the universe is obtained in a subset of Cardassian models by using WMAP data. Cardassian expansion is a modification to the Friedmann equation that allows the universe to be flat, matter dominated, and accelerating, without a vacuum component. Since this model changes the evolution of the universe, we should not a priori expect the Cardassian age to be the same as the WMAP Friedmann derived result of 13.7 +/- 0.2 Gyrs. However, in the subset of Cardassian models we consider, we discover that the age of the universe varies from 13.4 - 13.8 Gyr over the range of parameter space we explore, a result close to that of the standard Lambda Cold Dark Matter model. The Hubble constant h, which may also vary in these models, likewise varies little from the Friedmann result.

We classify the future of the universe for general cosmological models including matter and dark energy. If the equation of state of dark energy is less then -1, the age of the universe becomes finite. We compute the rest of the age of the universe for such universe models. The behaviour of the future growth of matter density perturbation is also studied. We find that the collapse of spherical overdensity region is greatly changed if the equation of state of dark energy is less than -1.

Protons can be plausibly accelerated to $\sim 10^{18}$-$10^{19}$ eV in strong accretion shocks surrounding massive clusters of galaxies. Such protons efficiently produce very high energy pairs through interactions with the CMB, which then radiate synchrotron and inverse Compton emission with hard spectra, peaking respectively in hard X-rays and TeV gamma-rays. We discuss the prospects for detecting this emission from nearby clusters with current and future facilities such as HESS and NeXT, which should provide a clear signature of ultra-high energy proton acceleration and an important probe of magnetic fields in the outermost regions of clusters.

We present a Monte-Carlo model of He II reionization by QSOs and its effect on the thermal state of the clumpy intergalactic medium (IGM). The model assumes that patchy reionization develops as a result of the discrete distribution of QSOs. It includes various recipes for the propagation of the ionizing photons, and treats photo-heating self-consistently. The model provides the fraction of He III, the mean temperature in the IGM, and the He II mean optical depth -- all as a function of redshift. It also predicts the evolution of the local temperature versus density relation during reionization. Our findings are as follows: The fraction of He III increases gradually until it becomes close to unity at $z\sim 2.8-3.0$. The He II mean optical depth decreases from $\tau\sim 10$ at $z\geq 3.5$ to $\tau\leq 0.5$ at $z\leq 2.5$. The mean temperature rises gradually between $z\sim 4$ and $z\sim 3$ and declines slowly at lower redshifts. The model predicts a flattening of the temperature-density relation with significant increase in the scatter during reionization at $z\sim 3$. Towards the end of reionization the scatter is reduced and a tight relation is re-established. This scatter should be incorporated in the analysis of the Ly$\alpha$ forest at $z\leq 3$. Comparison with observational results of the optical depth and the mean temperature at moderate redshifts constrains several key physical parameters.

All sufficiently massive clusters of galaxies are expected to be surrounded by strong accretion shocks, where protons can be accelerated to $\sim 10^{18}$-$10^{19}$ eV under plausible conditions. Such protons interact with the cosmic microwave background and efficiently produce very high energy electron-positron pairs, which then radiate synchrotron and inverse Compton emission, peaking respectively at hard X-ray and TeV gamma-ray energies. Characterized by hard spectra (photon indices $\sim 1.5$) and spatial distribution tracing the accretion shock, these can dominate over other nonthermal components depending on the shock magnetic field. HESS and other Cerenkov telescopes may detect the TeV emission from nearby clusters, notwithstanding its extended nature. The hard X-rays may be observable by future imaging facilities such as NeXT, and possibly also by ASTRO-E2/HXD. Such detections will not only provide a clear signature of ultra-high energy proton acceleration, but also an important probe of the accretion shock itself as well as magnetic fields in the outermost regions of clusters.

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts rates owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

We investigate the expected gravitational wave emission from coalescing supermassive black hole (SMBH) binaries resulting from mergers of their host galaxies. When galaxies merge, the SMBHs in the host galaxies sink to the center of the new merged galaxy and form a binary system. We employ a semi-analytic model of galaxy and quasar formation based on the hierarchical clustering scenario to estimate the amplitude of the expected stochastic gravitational wave background owing to inspiraling SMBH binaries and bursts owing to the SMBH binary coalescence events. We find that the characteristic strain amplitude of the background radiation is $h_c(f) \sim 10^{-16} (f/1 \mu {\rm Hz})^{-2/3}$ for $f \lesssim 1 \mu {\rm Hz}$ just below the detection limit from measurements of the pulsar timing provided that SMBHs coalesce simultaneously when host galaxies merge. The main contribution to the total strain amplitude of the background radiation comes from SMBH coalescence events at $0<z<1$. We also find that a future space-based gravitational wave interferometer such as the planned \textit{Laser Interferometer Space Antenna} ({\sl LISA}) might detect intense gravitational wave bursts associated with coalescence of SMBH binaries with total mass $M_{\rm tot} < 10^7 M_{\odot}$ at $z \gtrsim 2$ at a rate $ \sim 1.0 {\rm yr}^{-1}$. Our model predicts that burst signals with a larger amplitude $h_{\rm burst} \sim 10^{-15}$ correspond to coalescence events of massive SMBH binary with total mass $M_{\rm tot} \sim 10^8 M_{\odot}$ at low redshift $ z \lesssim 1$ at a rate $ \sim 0.1 {\rm yr}^{-1}$ whereas those with a smaller amplitude $h_{\rm burst} \sim 10^{-17}$ correspond to coalescence events of less massive SMBH binary with total mass $M_{\rm tot} \sim 10^6 M_{\odot}$ at high redshift $ z \gtrsim 3$.

Positronium (the bound state of electron and positron) has been thought to be formed after proton decay ($>10^{34}$yr) through collisional recombination and then decays by pair annihilation, thereby changing the matter content of the universe. We revisit the issue of the formation of positronium in the long-term future of the universe in light of recent indication that the universe is dominated by dark energy and dark matter. We find that if the equation of state of dark energy $w$ is less than -1/3 (including the cosmological constant $w=-1$), then the formation of positronium would not be possible, while it is possible through bound-bound transitions for $-1/3\siml w\siml-0.2$, or through collisional recombination for $w\simg-0.2$. The radiation from $\epm$ pair annihilation cannot dominate over $\epm$, while that from proton decay will dominate over baryon and $\epm$ for a while but not over dark matter.

We present observational constraints on a scalar-tensor gravity theory by $\chi^2$ test for CMB anisotropy spectrum. We compare the WMAP temperature power spectrum with the harmonic attractor model, in which the scalar field has its harmonic effective potential with curvature $\beta$ in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time. We found that the present value of the scalar coupling, i.e. the present level of deviation from Einstein gravity $(\alpha_0^2)$, is bounded to be smaller than $5\times 10^{-4-7\beta}$ ($2\sigma$), and $10^{-2-7\beta}$ ($4\sigma$) for $0< \beta<0.45$. This constraint is much stronger than the bound from the solar system experiments for large $\beta$ models, i.e., $\beta> 0.2$ and 0.3 in $2\sigma$ and $4\sigma$ limits, respectively. Furthermore, within the framework of this model, the variation of the gravitational constant at the recombination epoch is constrained as $|G(z=z_{rec})-G_0|/G_0 < 0.05(2\sigma)$, and $0.23(4\sigma)$.

The early reionisation of the Universe inferred from the WMAP polarisation results, if confirmed, poses a problem for the hypothesis that scale-invariant adiabatic density fluctuations account for large-scale structure and galaxy formation. One can only generate the required amount of early star formation if extreme assumptions are made about the efficiency and nature of early reionisation. We develop an alternative hypothesis that invokes an additional component of a non-scale-free isocurvature power spectrum together with the scale-free adiabatic power spectrum for inflation-motivated primordial density fluctuations. Such a component is constrained by the Lyman alpha forest observations, can account for the small-scale power required by spectroscopic gravitational lensing, and yields a source of early star formation that can reionise the universe at z~20 yet becomes an inefficient source of ionizing photons by z~10, thereby allowing the conventional adiabatic fluctuation component to reproduce the late thermal history of the intergalactic medium.

We show that UV photons produced by decaying particles can partially reionize the universe and explain the large optical depth observed by Wilkinson Microwave Anisotropy Probe. Together with UV fluxes from early formed stars and quasars, it is possible that the universe is fully ionized at z \lesssim 6 and partially ionized at z \gtrsim 6 as observed by Sloan Digital Sky Survey for large parameter space of the decaying particle. This scenario will be discriminated by future observations, especially by the EE polarization power spectrum of cosmic microwave background radiation.

The first-year Wilkinson Microwave Anisotropy Probe data suggest a high optical depth for Thomson scattering of 0.17 +/- 0.04, implying that the universe was reionized at an early epoch, z ~ 20. Such early reionization is likely to be caused by UV photons from first stars, but it appears that the observed high optical depth can be reconciled within the standard structure formation model only if star-formation in the early universe was extremely efficient. With normal star-formation efficiencies, cosmological models with non-Gaussian density fluctuations may circumvent this conflict as high density peaks collapse at an earlier epoch than in models with Gaussian fluctuations. We study cosmic reionization in non-Gaussian models and explore to what extent, within available constraints, non-Gaussianities affect the reionization history. For mild non-Gaussian fluctuations at redshifts of 30 to 50, the increase in optical depth remains at a level of a few percent and appears unlikely to aid significantly in explaining the measured high optical depth. On the other hand, within available observational constraints, increasing the non-Gaussian nature of density fluctuations can easily reproduce the optical depth and may remain viable in underlying models of non-Gaussianity with a scale-dependence.

We compare three independent, cosmological linear perturbation theory codes to asses the level of agreement between them and to improve upon it by investigating the sources of discrepancy. By eliminating the major sources of numerical instability the final level of agreement between the codes was improved by an order of magnitude. The relative error is now below 0.1% for the dark matter power spectrum. For the cosmic microwave background anisotropies the agreement is below the sampling variance up to l=3000, with close to 0.1% accuracy reached over most of this range of scales. The same level of agreement is also achieved for the polarization spectrum and the temperature-polarization cross-spectrum. Linear perturbation theory codes are thus well prepared for the present and upcoming high precision cosmological observations.

We critically examine how the evolution of the matter density field in cosmological simulations is affected by details of setting up initial conditions. We show that it is non-trivial to realise an initial distribution of matter in N-body/hydrodynamic simulations so that the baryon and dark matter density fluctuations and their velocities evolve consistently as theoretically predicted. We perform a set of cosmological simulations and use them to distinguish and verify an appropriate method for generating initial conditions. We show that a straightforward way of applying the Zel'dovich approximation to each component using distinct transfer functions results in an incorrect growth of density fluctuations and that it is necessary to correct velocities at the initial epoch. The unperturbed uniform particle distribution must be also generated appropriately to avoid tight coupling of the baryonic and dark matter components. We recommend using independent ``glass'' particle distributions, using distinct transfer functions for baryons and dark matter, and taking into account the difference in the velocity fields at the initialisation epoch. The proposed method will be useful for studies of the evolution of the intergalactic medium and the formation of the first cosmological objects using numerical simulations.

(abridged) We use large cosmological simulations to study the origin of primordial star-forming clouds in a Lambda CDM universe, by following the formation of dark matter halos and the cooling of gas within them. To model the physics of chemically pristine gas, we employ a non-equilibrium treatment of the chemistry of 9 species and include cooling by molecular hydrogen. We explore the hierarchical growth of bound structures forming at redshifts z = 25 - 30 with total masses in the range 10^5 - 10^6 Msun. The complex interplay between the gravitational formation of dark halos and the thermodynamic and chemical evolution of the gas clouds compromises analytic estimates of the critical H2 fraction. Dynamical heating from mass accretion and mergers opposes relatively inefficient cooling by molecular hydrogen, delaying the production of star-forming clouds in rapidly growing halos. We also investigate the impact of photo-dissociating ultra-violet (UV) radiation on the formation of primordial gas clouds. We consider two extreme cases by first including a uniform radiation field in the optically thin limit and secondly by accounting for the maximum effect of gas self-shielding in virialized regions. In both the cases we consider, the overall impact can be described by computing an equilibrium H2 abundance for the radiation flux and defining an effective shielding factor. Based on our numerical results, we develop a semi-analytic model of the formation of the first stars, and demonstrate how it can be coupled with large N-body simulations to predict the star formation rate in the early universe.

We reexamine constraints on the spatial size of closed toroidal models with cold dark matter and the cosmological constant from cosmic microwave background. We carry out Bayesian analyses using the Cosmic Background Explorer (COBE) data properly taking into account the statistically anisotropic correlation, i.e., off-diagonal elements in the covariance. We find that the COBE constraint becomes more stringent in comparison with that using only the angular power spectrum, if the likelihood is marginalized over the orientation of the observer. For some limited choices of orientations, the fit to the COBE data is considerably better than that of the infinite counterpart. The best-fit matter normalization is increased because of large-angle suppression in the power and the global anisotropy of the temperature fluctuations. We also study several deformed closed toroidal models in which the fundamental cell is described by a rectangular box. In contrast to the cubic models, the large-angle power can be enhanced in comparison with the infinite counterparts if the cell is sufficiently squashed in a certain direction. It turns out that constraints on some slightly deformed models are less stringent. We comment on how these results affect our understanding of the global topology of our universe.

Intracluster magnetic fields with $ \sim \mu G $ strength induce Faraday rotation on the cosmic microwave background (CMB) polarization. Measurements of this effect can potentially probe the detailed structure of intracluster magnetic fields across clusters, since the CMB polarization is a continuously varying field on the sky, in contrast to the conventional method restricted by the limited number of radio sources behind or inside a cluster. We here construct a method for extracting information on magnetic fields from measurements of the effect, combined with possible observations of the Sunyaev-Zel'dovich effect and X-ray emission for the same cluster which are needed to reconstruct the electron density fields. Employing the high-resolution magneto-hydrodynamic simulations performed by Dolag, Bartelmann & Lesch (1999) as a realistic model of magnetized intracluster gas distribution, we demonstrate how our reconstruction technique can well reproduce the magnetic fields, i.e., the spherically average radial profiles of the field strength and the coherence length.

We use joint likelihood analyses of combinations of fifteen cosmic microwave background (CMB) anisotropy data sets from the DMR, UCSB South Pole 1994, Python I--III, ARGO, MAX 4 and 5, White Dish, OVRO, and SuZIE experiments to constrain cosmogonies. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029) h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the data favor Omega_0 \simeq 0.9--1 (0.4--0.6) flat-Lambda (open) models. The range in deduced Omega_0 values is partially a consequence of the different combinations of smaller-angular-scale CMB anisotropy data sets used in the analyses, but more significantly a consequence of whether the DMR quadrupole moment is accounted for or ignored in the analysis. For both flat-Lambda and open models, after marginalizing over all other parameters, a lower Omega_B h^2 \simeq 0.005--0.009 is favored. This is also marginally at odds with estimates from more recent CMB anisotropy data and some estimates from standard nucleosynthesis theory and observed light element abundances. For both sets of models a younger universe with t_0 \simeq 12--15 Gyr is favored, consistent with other recent non-CMB indicators. We emphasize that since we consider only a small number of data sets, these results are tentative. More importantly, the analyses here do not rule out the currently favored flat-Lambda model with Omega_0 \sim 0.3, nor the larger Omega_B h^2 values favored by some other data.

Evolution of primordial fluctuations in a Brans-Dicke type scalar-tensor gravity theory is comprehensively investigated. The harmonic attractor model, in which the scalar field has its harmonic effective potential in the Einstein conformal frame and the theory relaxes toward Einstein gravity with time, is considered. The evolution of adiabatic initial perturbations in flat SCDM models is examined from the radiation-dominated epoch up to the present. We discuss how the scalar-tensor gravity affects the evolution of metric and matter perturbations, mainly focusing on the observational consequences, i.e., the matter power spectrum and the power spectrum of cosmic microwave background temperature. We find that the early time deviation is characterized only by the large static gravitational constant while the late time behavior is qualitatively different from that in Einstein gravity because the time variation of the gravitational constant and its fluctuation have non-negligible effects. The attracting scalar-tensor gravity affects only small scale modes due to its attracting nature, the degree of which is far beyond the post-Newtonian deviation at the present epoch.

We use Owens Valley Radio Observatory (OVRO) cosmic microwave background (CMB) anisotropy data to constrain cosmological parameters. We account for the OVRO beamwidth and calibration uncertainties, as well as the uncertainty induced by the removal of non-CMB foreground contamination. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029)h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the OVRO data favors an open (spatially-flat-Lambda) model with Omega_0 ~= 0.33 (0.1). At the 2 sigma confidence level model normalizations deduced from the OVRO data are mostly consistent with those deduced from the DMR, UCSB South Pole 1994, Python I-III, ARGO, MAX 4 and 5, White Dish, and SuZIE data sets.

We present predictions for two statistical measures of the hydrogen reionization process at high redshift. The first statistic is the number of neutral segments identified in spectra of high redshift QSOs as a function of their length. The second is the cross-correlation of neutral regions with possible sources of ionizing radiation. These independent probes are sensitive to the topology of the ionized regions. If reionization proceeded from high to low density regions then the cross-correlation will be negative, while if voids were ionized first then we expect a positive correlation and a relatively small number of long neutral segments. We test the sensitivity of these statistics for reionization by stars in high redshift galaxies. The flux of ionizing radiation emitted from stars is estimated by identifying galaxies in an N-body simulation using a semi-analytic galaxy formation model. The spatial distribution of ionized gas is traced in various models for the propagation of the ionization fronts. A model with ionization proceeding from high to low density regions is consistent with the observations of Becker et al. (2001), while models in which ionization begins in the lowest density regions appear to be inconsistent with the present data.

We investigate a cross-correlation between the weak gravitational lensing field of the large-scale structure, $\kappa$, and the secondary temperature fluctuation field, $\Delta$, of cosmic microwave background (CMB) induced by Thomson scattering of CMB photons off the ionized medium in the mildly nonlinear structure. The cross-correlation is expected to observationally unveil the biasing relation between the dark matter and ionized medium distributions in the large-scale structure. We develop a formalism for calculating the cross-correlation function and its angular power spectrum based on the small angle approximation. As a result, we find that leading contribution to the cross-correlation comes from the secondary CMB fluctuation field induced by the quadratic Doppler effects of the bulk velocity field $v$ of ionized medium. The magnitude of the cross-correlation can be estimated as $\kaco{\kappa\Delta}\sim 10^{-10}$ on small angular scales ($l\simgt 1000$) under the currently favored cold CDM model of the structure formation and the simplest scenario of homogeneous reionization after a given redshift of $z_{\rm ion}\simgt 5$. Interestingly, we find that the density-modulated quadratic Doppler effect of $O(\delta_e v^2)$ produces the cross-correlation with the weak lensing field that has {\em linear} dependence on the electron density fluctuation field $\delta_e$ or equivalently on the biasing relation between the dark matter and electron distributions. In other words, the angular power spectrum could be both positive or negative, depending on the positive biasing or antibiasing, respectively. Detection of the cross-correlation thus offers a new opportunity to observationally understand the reionization history of intergalactic medium connected to the dark matter clustering.

The observed magnetic field of microgauss strength in clusters of galaxies should induce the Faraday rotation effect on the linearly polarized cosmic microwave background (CMB) radiation when the CMB radiation propagates through a cluster at low redshift. Employing the Press-Schechter prescription for the cluster abundance under the cold dark matter (CDM) scenario and a plausible isothermal $\beta$-model for the gas distribution, we compute angular power spectra of the CMB polarization fields including the Faraday rotation mixing effect under the simple assumption of uniform magnetic field configuration across a cluster. As a result, we find that a parity-odd B-type polarization pattern is statistically generated on the observed sky, even when the primary polarization only contains the parity-even E-mode component, such as in the case of pure scalar perturbations. The generated B-type polarization has a peak with amplitude of $\sim0.1 \mu{\rm K}(B_0/0.1 \mu{\rm G})(\nu_0/10 {\rm GHz})^{-2}$ at angular scales of $l\approx 1000$ for the currently favored adiabatic $\Lambda$CDM model. This result also implies that, if the magnetic field has a $0.5 \mu{\rm G}$ strength and we observe at lower frequencies such as $\nu_0\simlt 5 {\rm GHz}$, the secondary signal due to the Faraday rotation effect could be comparable to the magnitude of the primary polarization. The frequency dependence of the Faraday rotation can be then used to discriminate the effect from primary and other secondary signals on the CMB polarizations. Our results therefore offer a new empirical opportunity to measure or constrain the intracluster magnetic field in the average sense.

We study the signal in the Cosmic Microwave Background (CMB) polarization anisotropy resulting from patchy reionization. It is well known that the primordial polarization of the CMB is very sensitive to the details of reionization. Combining a semi-analytic model of galaxy formation, in which the optical depth to the reionization epoch is in the range 0.014 to 0.048, with a high resolution N-body simulation we find that reionization generates a peak with amplitude 0.05~0.15 \mu K at large angular scales. The position of this peak reveals the size of the horizon at reionization, whilst its amplitude is a measure of the optical depth to reionization. On small scales, ionized patches prior to full reionization create a second order polarization signal due to the coupling of the free electron density fluctuation with the quadrupole moment of the temperature anisotropy. Careful study reveals that the coupling generates the same power spectra for electric and magnetic modes, whose amplitudes of polarization anisotropies from this process are predicted to be ~10nK. The amplitude depends strongly on the total baryon density $\Omega_b$ and on the spatial correlations of the free electron density in the ionized regions, and weakly on the fraction of ionizing photons able to escape their source galaxy, $f_{esc}$. The first- and second-order signals are therefore sensitive to the details of how the reionization occurred. Detection of these signals will place important constraints on the reionization history of the Universe.

We study evolution of various statistical quantities of smoothed cosmic density and velocity fields using N-body simulations. The parameter $C\equiv <V^2\delta>/(<V^2> <\delta^2>)$ characterizes nonlinear coupling of these two fields and determines behavior of bulk velocity dispersion as a function of local density contrast. It is found that this parameter depends strongly on the smoothing scale even in quasi-linear regimes where the skewness parameter $S_3$ is nearly constant and close to the predicted value by the second-order perturbation theory. We also analyze weakly nonlinear effects caused by an adaptive smoothing known as the gather approach.

We present predictions for the reionization of the intergalactic medium by stars in high-redshift galaxies. We calculate ionizing luminosities of galaxies, including the effects of absorption by interstellar gas and dust on the escape fraction, in order to calculate the filling factor of ionized hydrogen in the IGM. For a $\Lambda{\rm CDM}$ cosmology and a physical calculation of the escape fraction, we find that the hydrogen in the IGM will be reionized at redshift z=6.1 if the IGM has uniform density, but only by z=4.5 if the IGM is clumped. We combine our model with an N-body simulation of the distribution of dark matter in the universe in order to calculate the secondary temperature anisotropies induced in the cosmic microwave background (CMB) by scattering off free electrons. The models predict a spectrum of secondary anisotropies covering a broad range of angular scales, with fractional temperature fluctuations $\sim 10^{-7}-10^{-6}$ on arcminute scales. The amplitude depends somewhat on the geometry of reionization, with models in which the regions of highest gas density are reionized first giving larger CMB fluctuations than the case where galaxies ionize surrounding spherical regions, and models where low density regions reionize first giving the smallest fluctuations. Measurement of these anisotropies can therefore put important constraints on the reionization process, in particular, the redshift evolution of the filling factor, and should be a primary objective of a next generation submillimeter telescope such as the Atacama Large Millimeter Array. (abrdiged)

One of the biggest mystery of the recent observation of cosmic microwave background anisotropies by the boomerang team is insignificance of the second acoustic peak in the angular power spectrum. It is very difficult to explain such a low amplitude without assuming the higher baryon density than predicted by the Big Bang Nucleosynthesis (BBN). Employing the double inflation model in supergravity, we show that the density fluctuations produced by this inflation model can produce a sufficient low second acoustic peak with the standard value of the baryon density from BBN. It is shown that these density fluctuations are also consistent with the observations of cluster abundances and galaxy distributions.

The evolution of baryon density perturbations on very small scales is investigated. In particular, the nonlinear growth induced by the radiation drag force from the shear velocity field on larger scales during the recombination epoch, which is originally proposed by Shaviv in 1998, is studied in detail. It is found that inclusion of the diffusion term which Shaviv neglected in his analysis results in rather mild growth whose growth rate is $\ll 100$ instead of enormous amplification $\sim 10^4$ of Shaviv's original claim since the diffusion suppresses the growth. The growth factor strongly depends on the amplitude of the large scale velocity field. The nonlinear growth mechanism is applied to density perturbations of general adiabatic cold dark matter (CDM) models. In these models, it has been found in the previous works that the baryon density perturbations are not completely erased by diffusion damping if there exists gravitational potential of CDM. With employing the perturbed rate equation which is derived in this paper, the nonlinear evolution of baryon density perturbations is investigated. It is found that: (1) The nonlinear growth is larger for smaller scales. This mechanism only affects the perturbations whose scales are smaller than $\sim 10^2M_\odot$, which are coincident with the stellar scales. (2) The maximum growth factors of baryon density fluctuations for various COBE normalized CDM models are typically less than factor 10 for $3-\sigma$ large scale velocity peaks. (3) The growth factor depends on $\Omega_{\rm b}$.

The cosmological implication of a double inflation model with hybrid + new inflations in supergravity is studied. The hybrid inflation drives an inflaton for new inflation close to the origin through supergravity effects and new inflation naturally occurs. If the total e-fold number of new inflation is smaller than $\sim 60$, both inflations produce cosmologically relevant density fluctuations. Both cluster abundances and galaxy distributions provide strong constraints on the parameters in the double inflation model assuming $\Omega_0=1$ standard cold dark matter scenario. The future satellite experiments to measure the angular power spectrum of the cosmic microwave background will make a precise determination of the model parameters possible.

The effect of a non-trivial topology on the temperature correlations on the cosmic microwave background (CMB) in a small compact hyperbolic universe with volume comparable to the cube of the curvature radius is investigated. Because the bulk of large-angular CMB fluctuations is produced at the late epoch in low $\Omega_0$ models, the effect of a long wavelength cut-off due to the periodic structure does not lead to the significant suppression of large-angular power as in compact flat models. The angular power spectra are consistent with the COBE data for $\Omega_0 \ge 0.1$.

The late-time entropy production by the massive particle decay induces the various cosmological effects in the early epoch and modify the standard scenario. We investigate the thermalization process of the neutrinos after the entropy production by solving the Boltzmann equations numerically. We find that if the large entropy are produced at t $\sim$ 1 sec, the neutrinos are not thermalized very well and do not have the perfect Fermi-Dirac distribution. Then the freeze-out value of the neutron to proton ratio is altered considerably and the produced light elements, especially He4, are drastically changed. Comparing with the observational light element abundances, we find that $T_R$ < 0.7 MeV is excluded at 95 % C.L. We also study the case in which the massive particle has a decay mode into hadrons. Then we find that $T_R$ should be a little higher, i.e. $T_R$ > 2.5 MeV - 4 MeV, for the hadronic branching ratio $B_h = 10^{-2} - 1$. Possible influence of late-time entropy production on the large scale structure formation and temperature anisotropies of cosmic microwave background is studied. It is expected that the future satellite experiments (MAP and PLANCK) to measure anisotropies of cosmic microwave background radiation temperature can detect the vestige of the late-time entropy production as a modification of the effective number of the neutrino species $N_{\nu}^{\rm eff}$.

We estimate the deformation of the cosmic microwave background radiation by the hot region (``cocoon'') around a radio galaxy. A simple model is adopted for cocoon evolution while the jet is on, and a model of evolution is constructed after the jet is off. It is found that at low redshift the phase after the jet is off is longer than the lifetime of the jets. The Compton y-parameter generated by cocoons is calculated with a Press-Schechter number density evolution. The resultant value of y is of the same order as the COBE constraint. The Sunyaev-Zeldovich effect due to cocoons could therefore be a significant foreground source of small angular scale anisotropies in the cosmic microwave background radiation.

We consider the effect of three species of neutrinos with nearly degenerate mass on the cosmic structure formation in a low matter-density universe within a hierarchical clustering scenario with a flat initial perturbation spectrum. The matching condition for fluctuation powers at the COBE scale and at the cluster scale leads to a strong upper limit on neutrino mass. For a flat universe with matter density parameter Omega=0.3, we obtain m_nu<0.6eV for the Hubble constant H_0<80km/s/Mpc. Allowing for the more generous parameter space limited by Omega<0.4, H_0<80km/s/Mpc and age t_0>11.5Gyr, the limit is 0.9eV.

The axion is the most attractive candidate to solve the strong CP problem in QCD. If it exists, the inflationary universe produces axion fluctuations which are mixtures of isocurvature and adiabatic fluctuations in general. We investigate how large isocurvature fluctuations are allowed or favored in order to explain observations of the large scale structure of the present universe. Generic flat universe models with mixed (isocurvature + adiabatic) density fluctuations are studied. It is found that the observations are consistent with the mixed fluctuation model if the ratio \alpha of the power spectrum of isocurvature fluctuations to that of adiabatic fluctuations is less than \sim 0.1. In particular, the mixed fluctuation model with \alpha \sim 0.05, total matter density \Omega_{0} =0.4, and Hubble parameter H_0=70 km/s/Mpc gives a very good fit to the observational data. Since the height of the acoustic peak in the angular power spectrum of the cosmic microwave background (CMB) radiation decreases significantly when the isocurvature fluctuations are present, the mixed fluctuation model can be tested in future satellite experiments. Ratios of the amplitude at the peak location to that at the COBE normalization scale for various models are given. Furthermore, we also obtain the amplitude of isocurvature fluctuations as a function of axion parameters and the Hubble parameter during the inflation. We discuss the axion fluctuations in some realistic inflation models and find that a significant amount of the isocurvature fluctuations are naturally produced.

We investigate cosmological effects concerning the late-time entropy production due to the decay of non-relativistic massive particles. The thermalization process of neutrinos after the entropy production is properly solved by using the Boltzmann equation. If a large entropy production takes place at late time t$\simeq$ 1 sec, it is found that a large fraction of neutrinos cannot be thermalized. This fact loosens the tight constraint on the reheating temperature T_R from the big bang nucleosynthesis and T_R could be as low as 0.5 MeV. The influence on the large scale structure formation and cosmic microwave background anisotropies is also discussed.

We use Python I, II, and III cosmic microwave background anisotropy data to constrain cosmogonies. We account for the Python beamwidth and calibration uncertainties. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with nonrelativistic-mass density parameter Omega_0 in the range 0.1--1, baryonic-mass density parameter Omega_B in the range (0.005--0.029) h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the combined Python data favors an open (spatially-flat-Lambda) model with Omega_0 simeq 0.2 (0.1). At the 2 sigma confidence level model normalizations deduced from the combined Python data are mostly consistent with those drawn from the DMR, UCSB South Pole 1994, ARGO, MAX 4 and 5, White Dish, and SuZIE data sets.

We use combinations of ten small-scale cosmic microwave background anisotropy data sets from the UCSB South Pole 1994, ARGO, MAX 4 and 5, White Dish and SuZIE experiments to constrain cosmogonies. We consider open and spatially-flat-Lambda cold dark matter cosmogonies. Marginalizing over all parameters but $\Omega_0$, the combined data favors an Omega_0 simeq 1 for both open and flat-Lambda models. Excluding the smallest angular scale SuZIE data, an Omega_0 simeq 0.3 (1.0) open (flat-Lambda) model is favored. Considering only multi-frequency data with error bars consistent with sample variance and noise considerations, an Omega_0 simeq 0.1 (1) open (flat-Lambda) model is favored. For both open and flat-$\Lambda$ models, after marginalizing over all the other parameters, a lower Omega_B h^2 (sim 0.005) or younger (t_0 sim 10 Gyr) universe is favored. However, the data do not rule out other values of Omega_0 in the flat-Lambda model and other values of Omega_B h^2 in both models. (abriged)

We study the effect of a void network on the CMB anisotropy in the Einstein-de Sitter background using Thompson &Vishniac's model. We consider comprehensively the Sacks-Wolfe effect, the Rees-Sciama effect and the gravitational lensing effect. Our analysis includes the model of primordial voids existing at recombination, which is realized in some inflationary models associated with a first-order phase transition. If there exist primordial voids whose comoving radius is larger than $\sim10h^{-1}$Mpc at recombination, not only the Sachs-Wolfe effect but also the Rees-Sciama effect is appreciable even for multipoles $l\lsim1000$ of the anisotropy spectrum. The gravitational lensing effect, on the other hand, slightly smoothes the primary anisotropy; quantitatively, our results for the void model are similar to the previous results for a CDM model. All the effects, together, would give some constraints on the configuration or origin of voids with high-resolution data of the CMB anisotropy.

The M-theory, the strong-coupling heterotic string theory, presents various interesting new phenomenologies. The M-theory bulk axion is one of these. The decay constant in this context is estimated as $F_a\simeq 10^{16}$ GeV. Direct searches for the M-theory axion seem impossible because of the large decay constant. However, we point out that large isocurvature fluctuations of the M-theory axion are obtained in a hybrid inflation model, which will most likely be detectable in future satellite experiments on anisotropies of cosmic microwave background radiation.

We study gravitational lensing statistics, matter power spectra and the angular power spectra of the cosmic microwave background (CMB) radiation in x-matter models. We adopt an equation-of-state of x-matter which can express a wide range of matter from the pressureless dust to the cosmological constant. A new ingredient in this model is the sound speed of the x-component in addition to the equation-of-state $w_0=p_{x0}/\rho_{x0}$. Except for the cosmological constant case, the perturbations of x-matter itself are considered. Our primary interest is in the effect of non-zero sound speed on the structure formation and the CMB spectra. It is found that there exist parameter ranges where x-matter models are consistent with all current observations. The x-matter generally leaves imprints in the CMB anisotropy and the matter power spectrum, which should be detectable in future observations.

We use data from the ARGO cosmic microwave background (CMB) anisotropy experiment to constrain cosmogonies. We account for the ARGO beamwidth and calibration uncertainties, and marginalize over the offset removed from the data. Our derived amplitudes of the CMB anisotropy detected by the ARGO experiment are smaller than those derived previously. We consider open and spatially-flat-Lambda cold dark matter cosmogonies, with clustered-mass density parameter Omega_0 in the range 0.1-1, baryonic-mass density parameter Omega_B in the range (0.005-0.029)h^{-2}, and age of the universe t_0 in the range (10--20) Gyr. Marginalizing over all parameters but Omega_0, the ARGO data favors an open (spatially-flat-Lambda) model with Omega_0= 0.23 (0.1). However, these numerical values are model dependent. At the 2 sigma confidence level model normalizations deduced from the ARGO data are consistent with those drawn from the UCSB South Pole 1994, MAX 4+5, White Dish, and SuZIE data sets. The ARGO open model normalizations are also consistent with those deduced from the DMR data. However, for most spatially-flat-Lambda models the DMR normalizations are more than 2 sigma above the ARGO ones.

Equations to investigate fluctuations in cosmic X-ray background radiation due to point-like sources at high-redshift are formulated in a systematic way. The angular power spectrum of X-ray background fluctuations is investigated from large-scales to small-scales in various cosmological models such as open universe models and models with the cosmological constant, assuming a simple evolution model of the sources. The effect of epoch-dependent bias is demonstrated for small-angle fluctuations. The contribution from shot noise fluctuations is also discussed.

Spectra of density perturbations produced during chaotic inflation are calculated, taking both adiabatic and isocurvature modes into account in a class of scalar-tensor theories of gravity in which the dilaton is metrically coupled. Comparing the predicted spectrum of the cosmic microwave background radiation anisotropies with the one observed by the COBE-DMR we calculate constraints on the parameters of these theories, which turn out to be stronger by an order-of-magnitude than those obtained from post-Newtonian experiments.

We have simulated cosmic microwave background (CMB) anisotropy maps for several COBE-DMR-normalized cold dark matter (CDM) cosmogonies, to make predictions for the upcoming MAP experiment. We have studied the sensitivity of the simulated MAP data to cosmology, sky coverage, and instrumental noise. With accurate knowledge of instrumental noise, MAP data will discriminate among the cosmogonies considered, and determine the topology of the initial fluctuations. A correlation function analysis of the simulated MAP data results in a very accurate measurement of the acoustic Hubble radius at decoupling. A low-density open CDM model with Omega_0=0.4 can be distinguished from the Omega_0=1 fiducial CDM model or a Lambda CDM model with > 99% confidence from the location of the acoustic "valley" in the correlation function. A genus analysis of the simulated MAP data indicates that in cosmogonies with Gaussian random-phase initial conditions, a shift of the zero-crossing point of the genus curve near the mean temperature threshold level nu=0 should not exceed |Delta nu| = 0.01 (0.04) when the total effective FWHM smoothing is 0.3 deg (1.0 deg). The asymmetry of the genus curve at the positive and negative threshold levels should not exceed |Delta g/g(nu=1,-1)| = 0.8% (4%) at 0.3 deg (1.0 deg) FWHM smoothing. Deviations of the observed MAP data in excess of these small values will be evidence for non-Gaussian behavior. The amplitude of the genus curve is a measure of the shape of the power spectrum at the smoothing scale. Even with the expected amount of instrumental noise and partial sky coverage (due to the Galaxy), the MAP data should allow discrimination amongst the cosmogonies considered at more than 99% confidence solely from a genus amplitude analysis.

We use data from the White Dish experiment to set limits on cosmic microwave background radiation anisotropies in open and spatially-flat-Lambda cold dark matter cosmogonies. We account for the White Dish calibration uncertainty, and marginalize over the offset and gradient removed from the data. Our 2-sigma upper limits are larger than those derived previously. These upper limits are consistent with those derived from the $COBE$-DMR data for all models tested.

The evolution of small scale cosmological perturbations is carefully re-examined. Through the interaction with photons via electrons, baryon perturbations show interesting behavior in some physical scales. Characteristic features of the evolution of baryon density fluctuations are discussed. In CDM models, it is found a power-law growing phase of the small-scale baryon density fluctuations, which is characterized by the terminal velocity, after the diffusion (Silk) damping and before the decoupling epoch. Then, a transfer function for total matter density fluctuations is studied by taking into account those physical processes. An analytic transfer function is presented, which is applicable for the entire range up to a solar mass scale in the high$-z$ universe, and it is suitable also to the high baryon fraction models.

The fluid equations for the baryon-electron system in an expanding universe are derived from the Boltzmann equation. The effect of the Compton interaction is taken into account properly in order to evaluate the photon-electron collisional term. As an application, the acoustic motions of the baryon-electron system after recombination are investigated. The effective adiabatic index $\gamma$ is computed for sound waves of various wavelengths, assuming the perturbation amplitude is small. The oscillations are found to be dumped when $\gamma$ changes from between 1 (for an isothermal process) to 5/3 (for an adiabatic process).

We account for experimental and observational uncertainties in likelihood analyses of cosmic microwave background (CMB) anisotropy data from the MAX 4 and MAX 5 experiments. These analyses use CMB anisotropy spectra predicted in open and spatially-flat Lambda cold dark matter cosmogonies. Amongst the models considered, the combined MAX data set is most consistent with the CMB anisotropy shape in Omega_0 ~ 0.1-0.2 open models and less so with that in old (t_0 >~ 15 - 16 Gyr, i.e., low h), high baryon density (Omega_B >~ 0.0175/h^2), low density (Omega_0 ~ 0.2 - 0.4), flat-Lambda models. The MAX data alone do not rule out any of the models we consider at the 2-sigma level. Model normalizations deduced from the combined MAX data are consistent with those drawn from the UCSB South Pole 1994 data, except for the flat bandpower model where MAX favours a higher normalization. The combined MAX data normalization for open models with Omega_0 ~ 0.1-0.2 is higher than the upper 2-sigma value of the DMR normalization. The combined MAX data normalization for old (low h), high baryon density, low-density flat-Lambda models is below the lower 2-sigma value of the DMR normalization. Open models with Omega_0 ~ 0.4-0.5 are not far from the shape most favoured by the MAX data, and for these models the MAX and DMR normalizations overlap. The MAX and DMR normalizations also overlap for Omega_0 = 1 and some higher h, lower Omega_B, low-density flat-Lambda models.

Motivated by the possibility of $H_0t_0 > 1$ where $H_0$ and $t_0$ are the Hubble parameter and the age of the universe, respectively, we investigate the cosmology including x-matter. x-matter is expressed by the equation of state $p_x=w_0\rho_{x0} + c_s^2(\rho_x-\rho_{x0})$, where $p_x$, $\rho_x$ and $\rho_{x0}$ are the pressure, the density of x-matter and the density at present, respectively. $w_0$ and $c_s^2$ are functions of $\rho_x$ in general. x-matter has the most general form of the equation of state which is characterized by 1)violation of strong energy condition at present for $w_0 < -1/3$ ; 2) locally stable (i.e. $c_s^2\geq 0$); 3)causality is guaranteed ($c_s\leq 1$). Considering the case that $w_0$ and $c_s^2$ are constants, we find that there is a large parameter space of $(w_0,c_s^2,\Omega_{x0})$ in which the model universe is consistent with the age of the universe and the observations of distant Type I supernovae.

We use arcminute-scale data from the Sunyaev-Zel'dovich Infrared Experiment to set limits on anisotropies in the cosmic microwave background radiation in open and spatially-flat-Lambda cold dark matter cosmogonies. There are no 2-sigma detections for the models tested. The upper limits obtained are consistent with the amplitude of anisotropy detected by the COBE/DMR experiment.

We develop methods to account for experimental and observational uncertainties in likelihood analyses of data from cosmic microwave background (CMB) anisotropy experiments and apply them to an analysis of the UCSB South Pole 1994 (SP94) experiment. Observationally motivated open and spatially-flat $\Lambda$, cold dark matter cosmogonies are considered. Among the models we consider, the full SP94 data set is most consistent with $\Omega_0 \sim 0.1-0.2$ open models and less so with old $(t_0 \gap 15 - 16 Gyr)$, high baryon density $(\Omega_B \gap 0.0175 h^{-2})$, low density ($\Omega_0 \sim 0.2 - 0.4), flat-$\Lambda$ models. The SP94 data do not rule out any of the models we consider at the $2\sigma$ level. The SP94 experiment is most sensitive to anisotropies on a somewhat larger, model-dependent, angular scale than the scale at which the window function peaks. For establishing the significance of a detection of CMB anisotropy we derive limits using the highest posterior density (HPD) prescription, since it yields smaller lower limits. Since HPD limits lead to tighter constraints on the CMB amplitude, they also provide for greater discrimination between models. Model normalizations deduced from the SP94 data subsets are mostly consistent with those deduced from the two-year COBE-DMR data, although the Ka-band data prefer a normalization $\sim 1\sigma$ lower than do the Q-band data, the Q and Ka + Q data favour a slightly higher normalization for the $\Omega_0 = 0.1$ open model than does the DMR, and the Ka and Ka + Q data prefer a somewhat lower normalization for the older, higher $\Omega_B$, low-density $\Lambda$ models than does the DMR.

The mass density field as extracted from peculiar velocities in our cosmological neighborhood is mapped back in time to the CMB in two ways. First, the density power spectrum ($P_k$) is translated into a temperature angular power spectrum of sub-degree resolution ($C_l$) and compared to observations. Second, the local density field is translated into a temperature map in a patch on the last-scattering surface of a distant observer. A likelihood analysis of the Mark III peculiar velocity data have constrained the range of parameters for $P_k$ within the family of COBE-normalized CDM models (Zaroubi et al 1996), favoring a slight tilt in the initial spectrum, $n<1$. The corresponding range of $C_l$'s is plotted against current observations, indicating that the CMB data can tighten the constraints further: only models with ``small'' tilt ($n\sim 0.9$) and ``high'' baryonic content ($\Omega_b \sim 0.1$) could survive the two data sets simultaneously. The local mass density field that has been recovered from the velocities via a Wiener method is convolved with a Boltzmann calculation to recover $10'$ resolution temperature maps as viewed from different directions. The extent of the CMB patch and the amplitude of fluctuations depend on the choice of cosmological parameters, e.g., the local $100\hmpc$ sphere corresponds to $90'$ to $30'$ at the CMB for $\Omega$ between 1 and 0 respectively. The phases of the temperature map are correlated with those of the density field, contrary to the contribution of the Sachs-Wolfe effect alone. This correlation suggests the possibility of an inverse reconstruction of the underlying density field from CMB data with interesting theoretical implications.

Cut-sky orthogonal mode analyses of the COBE-DMR 53 and 90 GHz sky maps are used to determine the normalization of a variety of open cosmogonical models based on the cold dark matter scenario. To constrain the allowed cosmological-parameter range for these open cosmogonies, the predictions of the DMR-normalized models are compared to various observational measures of cosmography and large-scale structure, viz.: the age of the universe; small-scale dynamical estimates of the clustered-mass density parameter (Omega_0); constraints on the Hubble parameter (h), the X-ray cluster baryonic-mass fraction (Omega_B/Omega_0), and the matter power spectrum shape parameter; estimates of the mass perturbation amplitude; and constraints on the large-scale peculiar velocity field. The open-bubble inflation model is consistent with current determinations of the 95% confidence level range of these observational constraints.

The $ee\gamma\gamma +\rlap/E_T$ event observed by the CDF at Fermilab is naturally explained by dynamically supersymmetry breaking models and suggests the presence of the light gravitino which can be a warm dark matter. We consider large scale structure of the universe in the worm dark matter model and find that the warm dark matter plays almost the same role in the formation of the large scale structure as a cold dark matter if its mass is about $0.5$keV. We also study the Ly~$\alpha$ absorption systems which are presumed to be galaxies at high redshifts and show that the baryon density in the damped Ly~$\alpha$ absorption systems predicted by the warm dark matter model is quite consistent with the present observation.

Cosmic microwave background anisotropies provide a vast amount of cosmological information. Their full physical content and detailed structure can be understood in a simple and intuitive fashion through a systematic investigation of the individual mechanisms for anisotropy formation.

Through analytic techniques verified by numerical calculations, we establish general relations between the matter and cosmic microwave background (CMB) power spectra and their dependence on cosmological parameters on small scales. Fluctuations in the CMB, baryons, cold dark matter (CDM), and neutrinos receive a boost at horizon crossing. Baryon drag on the photons causes alternating acoustic peak heights in the CMB and is uncovered in its bare form under the photon diffusion scale. Decoupling of the photons at last scattering and of the baryons at the end of the Compton drag epoch, freezes the diffusion-damped acoustic oscillations into the CMB and matter power spectra at different scales. We determine the dependence of the respective acoustic amplitudes and damping lengths on fundamental cosmological parameters. The baryonic oscillations, enhanced by the velocity overshoot effect, compete with CDM fluctuations in the present matter power spectrum. We present new exact analytic solutions for the cold dark matter fluctuations in the presence of a growth- inhibiting radiation {\it and} baryon background. Combined with the acoustic contributions and baryonic infall into CDM potential wells, this provides a highly accurate analytic form of the small-scale transfer function in the general case.

A second--order perturbation theory approach is developed to calculate temperature anisotropies in the cosmic microwave background. Results are given for open universes and fluctuations corresponding to CDM models with either Harrison-Zeldovich (HZ) or Lyth-Stewart-Ratra-Peebles (LSRP) primordial energy--density fluctuation power spectrum. Our perturbation theory approach provides a distinctive multipole contribution as compared to the primary one, the amplitude of the effect being very dependent on normalization. For low--$\Omega$ models, the contribution of the secondary multipoles to the radiation power spectrum is negligible both for standard recombination and reionized scenarios, with the 2--year COBE--DMR normalization. For a flat universe this contribution is $\approx 0.1-10\%$ depending on the reionization history of the universe and on the normalization of the power spectrum.

Goodness-of-fit statistics are used to quantitatively establish the compatibility of CMB anisotropy predictions in a wide range of DMR-normalized, open and spatially-flat $\Lambda$, CDM cosmogonies with the set of all presently available small-scale CMB anisotropy detection data. Conclusions regarding model viability depend sensitively on the prescription used to account for the 1$\sigma$ uncertainty in the assumed value of the DMR normalization, except for low-density, $\Omega_0 \sim 0.3$ -- 0.4, open models which are compatible with the data for all prescriptions used. While large baryon-density ($\Omega_B \gap 0.0175 h^{-2}$), old ($t_0 \gap 15$ -- 16 Gyr), low-density ($\Omega_0 \sim 0.2$ -- 0.4), flat-$\Lambda$ models might be incompatible, no model is incompatible with the data for all prescriptions. In fact, some open models seem to fit the data better than should be expected, and this might be an indication that some error bars are mildly overconservative.

We compute the cosmic microwave background (CMB) anisotropy in a low-density, flat, cosmological constant, cold dark matter model which is normalized to the two-year COBE DMR sky map. Although conclusions regarding model viability must remain tentative until systematic effects are better understood, there are mild indications that these models have more intermediate scale power than is indicated by presently available CMB anisotropy observational data, with old ($t_0 \geq 15-16$Gyr), high baryon density ($\Omega_B \geq 0.0175 h^{-2}$), low density ($\Omega_0 \sim 0.2-0.4$) models doing the worst.

We compute the CMB anisotropy in an open inflation CDM model which is normalized to the two-year DMR sky map. Presently available CMB data is consistent with $\Omega_0 \sim 0.3 - 0.4$, but does not strongly disfavour either $\Omega_0 \sim 0.1$ or $\Omega_0 \sim 1$. Accumulating CMB data should strengthen the observational discriminative power, however, a better understanding of systematic effects will be needed before it will be possible to draw robust conclusions about model viability.

In the chaotic inflation models, quantum fluctuations for axion fields lead to the overproduction of domain walls and too large isocurvature fluctuations which is inconsistent with the observations of cosmic microwave background anisotropies. These problems are solved by assuming a very flat potential for the Peccei-Quinn scalar. As the simplest possibility, we consider a model where the Peccei-Quinn scalar is an inflaton itself and show that the isocurvature fluctuations can be comparable with the adiabatic ones. We investigate cosmological implications in the case that both adiabatic and isocurvature fluctuations exist and find that the amplitude of the matter spectrum becomes smaller than that for the pure adiabatic case. This leads to relatively high bias parameter ($b \simeq 2$) which is favoured by the current observations.

We have made a likelihood statistical analysis of the angular correlations in the {\it COBE}-DMR two-year sky maps by Monte Carlo simulation of the temperature fluctuations. We assume an open universe and consider as primordial power spectrum the Harrison-Zeldovich one, $P(k)=Ak$. We find that the flatness of the universe is not implied by the data. The quadrupole normalization amplitude, $Q_{rms-PS}$, is related to the density parameter, $\Omega$, by $Q_{rms-PS} = 10.67 + 55.81 \Omega - 128.59 \Omega^2 + 81.26 \Omega^3\ \mu$K. We have determined the p.d.f. of $\Omega$ due to cosmic plus sampling (i.e. $20^\circ$ galactic cut) variance which generically shows a bimodal shape. The uncertainty as given by the r.m.s. is $\approx 0.35$, therefore to better constrain $\Omega$ experiments sensitive to higher multipoles ($l>20$) should be considered.

As the data on CMB anisotropies improve and potential cosmological applications are realized, it will be increasingly important for theoretical calculations to be as accurate as possible. All modern calculations for inflationary-inspired fluctuations involve the numerical solution of coupled Boltzmann equations. There are many assumptions and choices to be made when carrying out such calculations. We go through each in turn, pointing out the best selections to make, and the level of inaccuracy expected through incorrect choice: (1) neglecting the effects of neutrinos or polarization has a 10% effect; (2) varying radiation temperature and He fraction can have smaller, but noticeable effects; (3) numerical issues, such as k-range and smoothing are discussed; (4) short-cut methods, e.g. free-streaming and tilt approximations, are generally inadequate at the few % level; (5) at the 1% level somewhat baroque effects are important, such as He recombination and even minimal reionization; (6) at smaller angular scales there are secondary effects which will ultimately have to be considered. Extracting information from the subsidiary acoustic peaks and the damping region will be an extremely challenging problem. However, given the prospect of measuring just such information on the sky, it will be important to meet this challenge. In principle we can obtain details about reionization history, neutrino contribution, helium abundance, non-power-law initial conditions etc.

Cosmic microwave background anisotropies provide a vast amount of information on both structure formation in the universe and the background dynamics and geometry. The full physical content and detailed structure of anisotropies can be understood in a simple and intuitive fashion through a systematic investigation of the individual mechanisms for anisotropy formation, based on elementary gravitational and fluid dynamics.

Fourier analysis on a cut sky of the two-year COBE DMR 53 and 90 GHz sky maps (in galactic coordinates) is used to determine the normalisation of an open inflation model based on the cold dark matter scenario. The normalised model is compared to measures of large-scale structure in the universe. Although the DMR data alone does not provide sufficient discriminative power to prefer a particular value of the mass density parameter, the open model appears to be reasonably consistent with observations when Omega_0 ~ 0.3 - 0.4 and merits further study.

We consider Bayesian constraints on standard isocurvature baryon models from the slope and normalization of the anisotropy power spectrum detected by the COBE DMR experiment in their two year maps. In conjunction with either the amplitude of matter fluctuations $\sigma_8$ {\it or} its slope, all open models are ruled out at greater than 95\% confidence, whereas cosmological constant $\Lambda$ dominated models are constrained to be highly ionized. By including the COBE FIRAS 95\% confidence upper limit on spectral distortion, we further reduce the available parameter space for $\Lambda$ models by excluding these highly ionized models. These constraints define a single class of barely viable standard models which makes definite and testable predictions for degree scale anisotropies and large scale structure.

Cosmic microwave background (CMB) anisotropies and density fluctuations are calculated for flat cold dark matter (CDM) models with a wide range of parameters, i.e., $\Omega_0, h$ and $\Omega_B$ for both standard recombination and various epochs of reionization. Tables of the power spectrum of CMB anisotropies in the form of $C_\ell$'s as a function of $\ell$ are presented. Although the Harrison-Zeldovich initial spectrum is assumed in these tables, we present simple approximations for obtaining the $C_\ell$'s corresponding to a tilted spectrum from those with a Harrison-Zeldovich spectrum. The $\sigma_8$ values are obtained for the matter density spectrum, with $\sigma(10^\circ)$, fixed $Q_{rms-PS}$ and COBE DMR 2 year normalizations. Simple modifications of the fitting formula of the density transfer function which are applicable for models with high baryon density are given. By using both numerical results and these fitting formulae, we calculate the relation between $\sigma_8$ and $Q_{rms-PS}$, and find good agreement. Velocity fields are also calculated.

Working toward a model independent understanding of cosmic microwave background (CMB) anisotropies and their significance, we undertake a comprehensive and self-contained study of scalar perturbation theory. Initial conditions, evolution, thermal history, matter content, background dynamics, and geometry all play a role in determining the anisotropy. By employing {\it analytic} techniques to illuminate the numerical results, we are able to separate and identify each contribution. We thus bring out the nature of the {\it total} Sachs-Wolfe effect, acoustic oscillations, diffusion damping, Doppler shifts, and reionization, as well as their particular manifestation in a critical, curvature, or cosmological constant dominated universe. By studying the full angular {\it and} spatial content of the resultant anisotropies, we isolate the signature of these effects from the dependence on initial conditions. Whereas structure in the Sachs-Wolfe anisotropy depends strongly on the underlying power spectra, the acoustic oscillations provide features which are nearly model independent. This may allow for future determination of the matter content of the universe as well as the adiabatic and/or isocurvature nature of the initial fluctuations.

We compare the two-year COBE DMR sky maps with the predictions of cosmological-constant cold dark matter models. Using a Bayesian analysis, we find that the most likely value of the cosmological constant in such a model is Lambda = 0. The data set an upper limit on Lambda of 0.78 (0.85) at 90% confidence, and 0.86 (0.92) at 95% confidence with (without) the quadrupole anisotropy.

We introduce a simple yet powerful {\it analytic} method which obtains the structure of cosmic microwave background anisotropies to better than 5-10\% in temperature fluctuations on {\it all} scales. It is applicable to {\it any} model in which the potential fluctuations at recombination are both linear and known. Moreover, it recovers and explains the presence of the ``Doppler peaks'' at degree scales as {\it driven} acoustic oscillations of the photon-baryon fluid. We treat in detail such subtleties as the time dependence of the gravitational driving force, anisotropic stress from the neutrino quadrupole, and damping during the recombination process, again all from an analytic standpoint. We apply this formalism to the standard cold dark matter model to gain physical insight into the anisotropies, including the dependence of the peak locations and heights on cosmological parameters such as $\Omega_b$ and $h$, as well as model parameters such as the ionization history. Damping due to the finite thickness of the last scattering surface and photon diffusion are further more shown to be identical. In addition to being a powerful probe into the nature of anisotropies, this treatment can be used in place of the standard Boltzmann code where 5-10\% accuracy in temperature fluctuations is satisfactory and/or speed is essential. Equally importantly, it can be used as a portable standard by which numerical codes can be tested and compared.

We compute the cosmic background radiation anisotropy, produced by energy-density fluctuations generated during an early epoch of inflation, in an open cosmological model based on the cold dark matter scenario. At $\Omega_0 \sim$ 0.3 -- 0.4, the COBE normalized open model appears to be consistent with most observations.

We investigate the dependence of large angular scale cosmic microwave background anisotropies on various initial conditions, including both adiabatic and isocurvature perturbations and the initial power law index $n$, in a variety of low-$\Omega$ cosmological models. Cosmological constant- dominated flat models and inflationary open adiabatic models, with $n_{\rm eff}<1$, and open isocurvature models, with $n_{\rm eff} \simeq 2$, are significantly constrained by the present observations.

We present a thorough and detailed investigation of baryon isocurvature models including realistic thermal histories, \ie, late or partial reionization of the universe and compact object formation after standard recombination. Constraints on these models are imposed from spectral distortion in the cosmic microwave background (CMB), number fluctuations in galaxy counts $\sigma_8$, and recent CMB anisotropy observations. Since the status of degree scale anisotropies is as yet unclear, the lack of spectral distortions is the most serious constraint on the spectral index of initial fluctuations: $-1.2 \simlt n \simlt -0.5$ in partially ionized models and $-1.2 \simlt n \simlt -0.9$ in compact object dominated universes in order to avoid significant bias in galaxy formation. Full ionization of the universe before $z = 300$ is also forbidden by this constraint. Intermediate scale CMB fluctuations which are significantly larger than standard CDM is a {\it robust} prediction of these models. Many models are consequently in conflict with the recent degree scale detections by Python, ARGO, and MSAM, \eg, compact object dominated models with $\Omega_0 \simlt 0.2$. Moreover, {\it all} models will be ruled out if the low CMB fluctuations on degree scales as detected by the SP91 13 point scan are confirmed. On the other hand, most models fare well compared with the high Tenerife and MAX detections.

( to appear in: Proceedings of the Nishonomiya Yukawa Memorial Symposium Edited by M. Sasaki) There are several models for generating fluctuations in an open universe that are compatible with the microwave background fluctuations detected by COBE {\it and} observations of large scale structure. Topological defects, such as strings and textures, appear to be more successful in a low-$\Omega$ model than in an $\Omega = 1$ model. These models predict that the cosmic microwave background (CMB) fluctuations are non-Gaussian with distinctive signatures that may lead to their confirmation. If there exists a mechanism for generating scale invariant mass fluctuations, $(\delta M/M)^2 \propto M^{-4/3}$, then these models are also compatible with COBE. Both models predict mass fluctuations on the 8/h Mpc scale, $\sigma_8$, $\sim 0.5 - 1$ for $\Omega h \sim 0.2 - 0.3$. A variety of observations of large-scale structure, clusters and galaxies, as well as the age problem, suggest that the low-$\Omega$ models are attractive.

We perform detailed calculations of cosmic microwave background (CMB) anisotropies in a CDM-dominated open universe with primordial adiabatic density perturbations for a variety of reionization histories. We show that to a great extent, the CMB anisotropies depend only on the geometry of the Universe, which in a matter dominated universe is determined by $\Omega$, and the optical depth to the surface of last scattering. In particular, the location of the primary Doppler peak depends primarily on $\Omega$ and is fairly insensitive to the other unknown parameters, such as $\Omega_b$, $h$, $\Lambda$, and the shape of the power spectrum. Therefore, measurements of CMB anisotropies on small scales may be used to determine $\Omega$.

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.

The effects of reionization, occurring after standard recombination in cold dark matter-dominated models, on CMB anisotropies are investigated. Late-time reionization reduces the CMB anisotropies, in particular, on degree scales. It is found that constraints on cold dark matter-dominated models from the highest frequency channel of the 9-point South Pole data are significantly relaxed for models which are consistent with Big Bang nucleosynthesis if reionization is assumed to have occurred by redshift $\sim 20.$

The integrated Sachs-Wolfe effect (ISW) can be an important factor in the generation of Cosmic Microwave Background anisotropies on all scales, especially in a reionized curvature or lambda dominated universe. We present an analytic treatment of the ISW effect, which is analogous to thick last scattering surface techniques for the Doppler effect, that compares quite well with the full numerical calculations. The power spectrum of temperature fluctuations due to the small scale ISW effect has wave number dependence k^{-5} times that of the matter power spectrum.