We study the future prospects of the 21cm forest observations on the axion-like dark matter when the spontaneous breaking of the global Peccei-Quinn (PQ) symmetry occurs after the inflation. The large isocurvature perturbations of order unity sourced from axion-like particles can result in the enhancement of minihalo formation, and the subsequent hierarchical structure formation can affect the minihalo abundance whose masses can exceed ${\cal O}(10^4) M_{\odot}$ relevant for the 21cm forest observations. We show that the 21cm forest observations are capable of probing the axion-like particle mass in the range $10^{-18}\lesssim m_a \lesssim 10^{-12}$ eV for the temperature independent axion mass. For the temperature dependent axion mass, the zero temperature axion mass scale for which the 21cm forest measurements can be affected is extended further to as big as of order $10^{-6}$ eV.

The ultra-light scalar fields can arise ubiquitously, for instance, as a result of the spontaneous breaking of an approximate symmetry such as the axion and more generally the axion-like particles. In addition to the particle physics motivations, these particles can also play a major role in cosmology by contributing to dark matter abundance and affecting the structure formation at sub-Mpc scales. In this paper, we propose to use the 21cm forest observations to probe the nature of ultra-light dark matter. The 21cm forest can probe much smaller scales than the Lyman-$\alpha$ forest, that is, $k\gtrsim 10\mathrm{Mpc}^{-1}$. We explore the range of the ultra-light dark matter mass $m_{u}$ and $f_u$, the fraction of ultra-light dark matter with respect to the total matter, which can be probed by the 21cm forest. We find that 21cm forest can potentially put the dark matter mass lower bound $m_u \gtrsim 10^{-18}$ eV for $f_u=1$, which is 3 orders of magnitude bigger mass scale than those probed by the current Lyman-$\alpha$ forest observations.While the effects of the ultra-light particles on the structure formation become smaller when the dominant component of dark matter is composed of the conventional cold dark matter, we find that the 21cm forest is still powerful enough to probe the sub-component ultra-light dark matter mass up to the order of $10^{-19}$ eV. The Fisher matrix analysis shows that $(m_u,f_u)\sim (10^{-20}\mathrm{eV}, 0.3)$ is the most optimal parameter set which the 21cm forest can probe with the minimal errors for a sub-component ultra-light dark matter scenario.

Recent developments of transition-edge sensors (TESs), based on extensive experience in ground-based experiments, have been making the sensor techniques mature enough for their application on future satellite CMB polarization experiments. LiteBIRD is in the most advanced phase among such future satellites, targeting its launch in Japanese Fiscal Year 2027 (2027FY) with JAXA's H3 rocket. It will accommodate more than 4000 TESs in focal planes of reflective low-frequency and refractive medium-and-high-frequency telescopes in order to detect a signature imprinted on the cosmic microwave background (CMB) by the primordial gravitational waves predicted in cosmic inflation. The total wide frequency coverage between 34GHz and 448GHz enables us to extract such weak spiral polarization patterns through the precise subtraction of our Galaxy's foreground emission by using spectral differences among CMB and foreground signals. Telescopes are cooled down to 5Kelvin for suppressing thermal noise and contain polarization modulators with transmissive half-wave plates at individual apertures for separating sky polarization signals from artificial polarization and for mitigating from instrumental 1/f noise. Passive cooling by using V-grooves supports active cooling with mechanical coolers as well as adiabatic demagnetization refrigerators. Sky observations from the second Sun-Earth Lagrangian point, L2, are planned for three years. An international collaboration between Japan, USA, Canada, and Europe is sharing various roles. In May 2019, the Institute of Space and Astronautical Science (ISAS), JAXA selected LiteBIRD as the strategic large mission No. 2.

We propose a new non-parametric method to constrain the cosmological model through the growth factor of large-scale structure. To constrain the cosmological model from observations such as cosmic microwave background or large-scale structure, the power spectrum or correlation function have been used since they are well described by the linear perturbation theory or little correction is needed for the non-linearity; however, it is always important to find an alternative route to constrain cosmological models in order to avoid systematic effects due to the estimator and make the constraint more robust.Our density fluctuation reconstruction method constrains the cosmological model directly from time evolution of density fluctuation. We propose to apply Kolmogorov-Smirnov test and D'Agostino's K-squared test to the distribution function of density fluctuations, which are observed on a light cone. The density fluctuation recovers the Gaussian distribution only if we reconstruct the density fluctuation with correct linear growth factors. The method is verified by the use of Gaussian random field and N-body simulations. This method will open a new window both to separate the initial three dimensional density fluctuation and the growth factor of large scale structure.

We study the cross-correlation between the 21cm and CMB B mode fluctuations which are induced by the cosmic birefringence when the non-constant scalar field couples to the electromagnetic field strength. Such multi-wavelength signals can potentially probe the reionization history of the Universe and also explore the nature of fundamental physics such as the parity violation and the scalar field dynamics in the early Universe. We illustrate the feasibility to detect such 21cm-B mode cross-correlations through commonly discussed scalar field models, the quintessence-like scalar field which is responsible for the current dark energy and the axion-like scalar field which is responsible for the current dark matter density.

Correlation analysis is obtained among Faraday rotation measure, HI column density, thermal and synchrotron radio brightness using archival all-sky maps of the Galaxy. A method is presented to calculate the magnetic strength and its line-of-sight (LOS) component, volume gas densities, effective LOS depth, effective scale height of the disk) from these data in a hybrid way. Applying the method to archival data, all-sky maps of the local magnetic field strength and its parallel component are obtained, which reveal details of local field orientation.

We show that the cosmic birefringence and miscalibrated polarisation angles can be determined simultaneously by cosmic microwave background (CMB) experiments using the cross-correlation between $E$- and $B$-mode polarisation data. This is possible because polarisation angles of the CMB are rotated by both the cosmic birefringence and miscalibration effects,whereas those of the Galactic foreground emission only by the latter. Our method does not require prior knowledge of the $E$- and $B$-mode power spectra of the foreground emission, but uses only the knowledge of the CMB polarisation spectra. Specifically, we relate the observed $EB$ correlation to the difference between the $\mathit{observed}$ $E$- and $B$-mode spectra in the sky, and use different multipole dependence of the CMB (given by theory) and foreground spectra (with no assumption) to derive the likelihood for the miscalibration angle $\alpha$ and the birefringence angle $\beta$. We show that a future satellite mission similar to LiteBIRD can determine $\beta$ with a precision of ten arcminutes.

We extend the internal template foreground removal method by accounting for spatially varying spectral parameters such as the spectral indices of synchrotron and dust emission and the dust temperature. As the previous algorithm had to assume that the spectral parameters are uniform over the full sky (or some significant fraction of the sky), it resulted in a bias in the tensor-to-scalar ratio parameter $r$ estimated from foreground-cleaned polarization maps of the cosmic microwave background (CMB). The new algorithm, "Delta-map method", accounts for spatially varying spectra to first order in perturbation. The only free parameters are the cosmological parameters such as $r$ and the sky-averaged foreground parameters. We show that a cleaned CMB map is the maximum likelihood solution to first order in perturbation, and derive the posterior distribution of $r$ and the sky-averaged foreground parameters using Bayesian statistics. Applying this to realistic simulated sky maps including polarized CMB, synchrotron and thermal dust emission, we find that the new algorithm removes the bias in $r$ down to undetected level in our noiseless simulation ($r\lesssim 4\times 10^{-4}$). We show that the frequency decorrelation of polarized foregrounds due to averaging of spatially varying spectra can be accounted for to first order in perturbation by using slightly different spectral parameters for the Stokes $Q$ and $U$ maps. Finally, we show that the effect of polarized anomalous microwave emission on foreground removal can be absorbed into the curvature parameter of the synchrotron spectrum.

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.

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.

Observation of the 21cm line signal from neutral hydrogen during the Epoch of Reionization is challenging due to extremely bright Galactic and extragalactic foregrounds and complicated instrumental calibration. A reasonable approach for mitigating these problems is the cross correlation with other observables. In this work, we present the first results of the cross power spectrum (CPS) between radio images observed by the Murchison Widefield Array and the cosmic microwave background (CMB), measured by the Planck experiment. We study the systematics due to the ionospheric activity, the dependence of CPS on group of pointings, and frequency. The resulting CPS is consistent with zero because the error is dominated by the foregrounds in the 21cm observation. Additionally, the variance of the signal indicates the presence of unexpected systematics error at small scales. Furthermore, we reduce the error by one order of magnitude with application of a foreground removal using a polynomial fitting method. Based on the results, we find that the detection of the 21cm-CMB CPS with the MWA Phase I requires more than 99.95% of the foreground signal removed, 2000 hours of deep observation and 50% of the sky fraction coverage.

We explore the possibility that both the suppression of the $\ell = 2$ multipole moment of the power spectrum of cosmic microwave background temperature fluctuations and the possible dip for $\ell = 10-30$ can be explained as well as a possible new dip for $\ell \approx 60$ as the result of the resonant creation of sequential excitations of a fermionic (or bosonic) closed superstring that couples to the inflaton field. We consider a D=26 closed bosonic string with one toroidal compact dimension as an illustration of how string excitations might imprint themselves on the CMB. We analyze the existence of successive momentum states, winding states or oscillations on the string as the source of the three possible dips in the power spectrum. Although the evidence of these dips are of marginal statistical significance, this might constitute the first observational evidence of successive superstring excitations in Nature.

In this paper, we present the effects of dark energy perturbation on the formation and abundance of cosmic voids. We consider dark energy to be a fluid with a negative pressure characterised by a constant equation of state w and speed of sound $c^2_s$. By solving fluid equations for two components, namely, dark matter and dark energy fluids, we quantify the effects of dark energy perturbation on the sizes of top-hat voids. We also explore the effects on the size distribution of voids based on the excursion set theory. We confirm that dark energy perturbation negligibly affects the size evolution of voids;$c^2_s = 0$ varies the size only by 0.1% as compared to the homogeneous dark energy model. We also confirm that dark energy perturbation suppresses the void size when w < -1 and enhances the void size when $w > -1$ (Basse et al. 2011). In contrast to the negligible impact on the size, we find that the size distribution function on scales larger than 10 Mpc/h highly depends on dark energy perturbation; compared to the homogeneous dark energy model, the number of large voids of radius 30Mpc is 25% larger for the model with $w = -0.9$ and $c^2_s = 0$ while they are 20% less abundant for the model with $w = -1.3$ and $c^2_s = 0$.

Inflation is the leading theory of the first instant of the universe. Inflation, which postulates that the universe underwent a period of rapid expansion an instant after its birth, provides convincing explanation for cosmological observations. Recent advancements in detector technology have opened opportunities to explore primordial gravitational waves generated by the inflation through B-mode (divergent-free) polarization pattern embedded in the Cosmic Microwave Background anisotropies. If detected, these signals would provide strong evidence for inflation, point to the correct model for inflation, and open a window to physics at ultra-high energies. LiteBIRD is a satellite mission with a goal of detecting degree-and-larger-angular-scale B-mode polarization. LiteBIRD will observe at the second Lagrange point with a 400 mm diameter telescope and 2,622 detectors. It will survey the entire sky with 15 frequency bands from 40 to 400 GHz to measure and subtract foregrounds. The U.S. LiteBIRD team is proposing to deliver sub-Kelvin instruments that include detectors and readout electronics. A lenslet-coupled sinuous antenna array will cover low-frequency bands (40 GHz to 235 GHz) with four frequency arrangements of trichroic pixels. An orthomode-transducer-coupled corrugated horn array will cover high-frequency bands (280 GHz to 402 GHz) with three types of single frequency detectors. The detectors will be made with Transition Edge Sensor (TES) bolometers cooled to a 100 milli-Kelvin base temperature by an adiabatic demagnetization refrigerator.The TES bolometers will be read out using digital frequency multiplexing with Superconducting QUantum Interference Device (SQUID) amplifiers. Up to 78 bolometers will be multiplexed with a single SQUID amplidier. We report on the sub-Kelvin instrument design and ongoing developments for the LiteBIRD mission.

We derive an observational constraint on a spherical inhomogeneity of the void centered at our position from the angular power spectrum of the cosmic microwave background(CMB) and local measurements of the Hubble parameter. The late time behaviour of the void is assumed to be well described by the so-called $\Lambda$-Lema\^itre-Tolman-Bondi~($\Lambda$LTB) solution. Then, we restrict the models to the asymptotically homogeneous models each of which is approximated by a flat Friedmann-Lema\^itre-Robertson-Walker model. The late time $\Lambda$LTB models are parametrized by four parameters including the value of the cosmological constant and the local Hubble parameter. The other two parameters are used to parametrize the observed distance-redshift relation. Then, the $\Lambda$LTB models are constructed so that they are compatible with the given distance-redshift relation. Including conventional parameters for the CMB analysis, we characterize our models by seven parameters in total. The local Hubble measurements are reflected in the prior distribution of the local Hubble parameter. As a result of a Markov-Chains-Monte-Carlo analysis for the CMB temperature and polarization anisotropies, we found that the inhomogeneous universe models with vanishing cosmological constant are ruled out as is expected. However, a significant under-density around us is still compatible with the angular power spectrum of CMB and the local Hubble parameter.

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 summarize some applications of big bang nucleosythesis (BBN) and the cosmic microwave background (CMB) to constrain the first moments of the creation of matter in the universe. We review the basic elements of BBN and how it constraints physics of the radiation-dominated epoch. In particular, how the existence of higher dimensions impacts the cosmic expansion through the projection of curvature from the higher dimension in the "dark radiation" term. We summarize current constraints from BBN and the CMB on this brane-world dark radiation term. At the same time, the existence of extra dimensions during the earlier inflation impacts the tensor to scalar ratio and the running spectral index as measured in the CMB. We summarize how the constraints on inflation shift when embedded in higher dimensions. Finally, one expects that the universe was born out of a complicated multiverse landscape near the Planck time. In these moments the energy scale of superstrings was obtainable during the early moments of chaotic inflation. We summarize the quest for cosmological evidence of the birth of space-time out of the string theory landscape. We will explore the possibility that a superstring excitations may have made itself known via a coupling to the field of inflation. This may have left an imprint of "dips" in the power spectrum of temperature fluctuations in the cosmic microwave background. The identification of this particle as a superstring is possible because there may be evidence for different oscillator states of the same superstring that appear on different scales on the sky. It will be shown that from this imprint one can deduce the mass, number of oscillations, and coupling constant for the superstring. Although the evidence is marginal, this may constitute the first observation of a superstring in Nature.

Magnetic field is ubiquitous in the Universe and it plays essential roles in various astrophysical phenomena, yet its real origin and evolution are poorly known. This article reviews current understanding of magnetic fields in the interstellar medium, the Milky Way Galaxy, external galaxies, active galactic nuclei, clusters of galaxies, and the cosmic web. Particularly, the review concentrates on the achievements that have been provided by centimeter and meter wavelength radio observations. The article also introduces various methods to analyze linear polarization data, including synchrotron radiation, Faraday rotation, depolarization, and Faraday tomography.

Understanding the detailed process of cosmic reionization is one of the remaining problems in astrophysics and cosmology. Here we construct a model of cosmic reionization that includes contribution from high-$z$ galaxies and active galactic nuclei (AGNs), and calculate the reionization and thermal histories with the model. To keep the model to be general and realistic, we vary the escape fraction of ionizing photons, $f_{\rm esc}$, and the faint-end slope of the AGN luminosity function at high redshifts, $\alpha_{\rm hz}$, within the constraints from the observed cosmic star formation history and the observed bright-end UV luminosity functions at $z \leq 6$. Additionally, we model the spectral energy distribution (SED) of AGNs which depends on the Eddington ratio and the black-hole mass. By comparing the computed reionization histories with the observed \ion{H}{i} fractions and the optical depth for Thomson scattering from Planck, we find that $\alpha_{\rm hz} > -1.5$ and $f_{\rm esc} < 0.15$ are favored when we employ the bright-end luminosity function obtained by Giallongo et al. Our result suggests that an AGN-dominated model with the abundance of faint-AGNs being as large as the estimate by Giallongo et al. is allowed only if the contribution from high-z galaxies is almost negligible, while a galaxy-dominant model is also allowed. We also find that the shape of SED has a significant impact on the thermal history. Therefore it is expected that measurements of the thermal state of the IGM provide useful information on properties of ionizing sources.

Cosmological defects result from cosmological phase transitions in the early Universe and the dynamics reflects their symmetry-breaking mechanisms. These cosmological defects may be probed through weak lensing effects because they interact with ordinary matters only through the gravitational force. In this paper, we investigate global textures by using weak lensing curl and B modes. Non-topological textures are modeled by the non-linear sigma model (NLSM), and induce not only the scalar perturbation but also vector and tensor perturbations in the primordial plasma due to the nonlinearity in the anisotropic stress of scalar fields. We show angular power spectra of curl and B modes from both vector and tensor modes based on the NLSM. Furthermore, we give the analytic estimations for curl and B mode power spectra. The amplitude of weak lensing signals depends on a combined parameter $\epsilon^{2}_{v} = N^{-1}\left( v/m_{\rm pl} \right)^{4}$ where $N$ and $v$ are the number of the scalar fields and the vacuum expectation value, respectively. We discuss the detectability of the curl and B modes with several observation specifications. In the case of the CMB lensing observation without including the instrumental noise, we can reach $\epsilon_{v} \approx 2.7\times 10^{-6}$. This constraint is about 10 times stronger than the current one determined from the Planck. For the cosmic shear observation, we find that the signal-to-noise ratio depends on the mean redshift and the observing number of galaxies as $\propto z^{0.7}_{\rm m}$ and $\propto N^{0.2}_{\rm g}$, respectively. In the study of textures using cosmic shear observations, the mean redshift would be one of the key design parameters.

While the observed nearly scale-invariant initial power spectrum is regarded as one of the favorable evidence of the standard inflationary cosmology, precision observations of the Cosmic Microwave Background (CMB) anisotropies also suggest possible existence of nontrivial features such as those observed around multipoles $\ell\sim120$ by WMAP. Here, we examine the Planck data and investigate the effects of these features on the cosmological parameter estimation performing the Markov-Chain Monte-Carlo (MCMC) analysis. We find that the features exist in the Planck data at the same position as the case of the WMAP data but they do not affect the cosmological parameter estimation significantly.

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 perform a Fisher analysis to estimate the expected constraints on the Epoch of Reionization (EoR) model parameters (i.e., minimum virial temperature, the ionizing efficiency and the mean free path of ionizing photons) taking into account the thermal noise of the exisitting telescopes, MWA and LOFAR. We consider how the inclusion of the 21cm bispectrum improves the constraints compared to using the power spectrum alone. Assuming that we perfectly remove the foreground, we found that the bispectrum, which is calculated by the 21cmFAST code, can constrain the EoR model parameters more tightly than the power spectrum since the bispectrum is more sensitive to the EoR model parameters than the power spectrum. We also found that degeneracy between the EoR model parameters can be broken by combining the bispectrum with the power spectrum.

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

We explore the possible cosmological consequence of the gravitational coupling between the inflaton and axion-like fields. In view of the forthcoming cosmic microwave background (CMB) polarization and lensing data, we study the sensitivity of the CMB data on the cross-correlation between the curvature and axion isocurvature perturbations. Through a concrete example, we illustrate the explicit dependence of the scale dependent cross-correlation power spectrum on the axion parameters.

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.

The power spectrum of the cosmic microwave background from both the {\it Planck} and {\it WMAP} data exhibits a slight dip for multipoles in the range of $l= 10-30$. We show that such a dip could be the result of the resonant creation of massive particles that couple to the inflaton field. For our best-fit models, the epoch of resonant particle creation reenters the horizon at a wave number of $k_* \sim 0.0011 \pm 0.0004 $ ($h$ Mpc$^{-1}$). The amplitude and location of this feature corresponds to the creation of a number of degenerate fermion species of mass $\sim (8-11) /\lambda^{3/2} $ $m_{pl}$ during inflation where $\lambda \sim (1.0 \pm 0.5) N^{-2/5}$ is the coupling constant between the inflaton field and the created fermion species, while $N$ is the number of degenerate species. Although the evidence is of marginal statistical significance, this could constitute new observational hints of unexplored physics beyond the Planck scale

Redshifted 21cm signal is a promising tool to investigate the state of intergalactic medium (IGM) in the Cosmic Dawn (CD) and Epoch of Reionization(EoR). In our previous work (Shimabukuro et al 2015), we studied the variance and skewness to give a clear interpretation of 21cm power spectrum and found that skewness is a good indicator of the epoch when X-ray heating becomes effective. Thus, the non-Gaussian feature of the spatial distribution of the 21cm signal is expected to be useful to investigate the astrophysical effects in the CD and EoR. In this paper, in order to investigate such a non-Gaussian feature in more detail, we focus on the bispectrum of the 21cm signal. It is expected that the 21cm brightness temperature bispectrum is produced by non-gaussianity due to the various astrophysical effects such as Wouthysen-Field (WF) effect, X-ray heating and reionization. We study the various properties of 21cm bispectrum such as scale dependence, shape dependence and redshift evolution. And also we study the contribution from each component of 21cm bispectrum. We find that the contribution from each component has characteristic scale-dependent feature, and it is expected that we could obtain more detailed information on the IGM in the CD and EoR by using the 21cm bispectrum in the future experiments, combined with the power spectrum and skewness.

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

Magnetic fields play essential roles in various astronomical objects. Radio astronomy has revealed that magnetic fields are ubiquitous in our Universe. However, the real origin and evolution of magnetic fields is poorly proven. In order to advance our knowledge of cosmic magnetism in coming decades, the Square Kilometre Array (SKA) should have supreme sensitivity than ever before, which provides numerous observation points in the cosmic space. Furthermore, the SKA should be designed to facilitate wideband polarimetry so as to allow us to examine sightline structures of magnetic fields by means of depolarization and Faraday Tomography. The SKA will be able to drive cosmic magnetism of the interstellar medium, the Milky Way, galaxies, AGN, galaxy clusters, and potentially the cosmic web which may preserve information of the primeval Universe. The Japan SKA Consortium (SKA-JP) Magnetism Science Working Group (SWG) proposes the project "Resolving 4-D Nature of Magnetism with Depolarization and Faraday Tomography", which contains ten scientific use cases.

Great attention is given to the first star formation and the epoch of reionization as main targets of planned large radio interferometries (e.g. Square Kilometre Array). Recently, it is claimed that the supersonic relative velocity between baryons and cold dark matter can suppress the abundance of first stars and impact the cosmological reionization process. Therefore, in order to compare observed results with theoretical predictions it is important to examine the effect of the supersonic relative motion on the small-scale structure formation. In this paper, we investigate this effect on the nonlinear structure formation in the context of the spherical collapse model in order to understand the fundamental physics in a simple configuration. We show the evolution of the dark matter sphere with the relative velocity by both using N-body simulations and numerically calculating the equation of motion for the dark matter mass shell. The effects of the relative motion in the spherical collapse model appear as the delay of the collapse time of dark matter halos and the decrease of the baryon mass fraction within the dark matter sphere. Based on these results, we provide the fitting formula of the critical density contrast for collapses with the relative motion effect and calculate the mass function of dark matter halos in the Press-Schechter formalism. As a result, the relative velocity decreases the abundance of dark matter halos whose mass is smaller than $10^8~M_\odot/h$.

The discrepancy between the amplitudes of matter fluctuations inferred from Sunyaev-Zel'dovich (SZ) cluster number counts, the primary temperature, and the polarization anisotropies of the cosmic microwave background (CMB) measured by the Planck satellite can be reconciled if the local universe is embedded in an under-dense region as shown by Lee, 2014. Here using a simple void model assuming the open Friedmann-Robertson-Walker geometry and a Markov Chain Monte Carlo technique, we investigate how deep the local under-dense region needs to be to resolve this discrepancy. Such local void, if exists, predicts the local Hubble parameter value that is different from the global Hubble constant. We derive the posterior distribution of the local Hubble parameter from a joint fitting of the Planck CMB data and SZ cluster number counts assuming the simple void model. We show that the predicted local Hubble parameter value of $H_{\rm loc}=70.1\pm0.34~{\rm km\,s^{-1}Mpc^{-1}}$ is in better agreement with direct local Hubble parameter measurements, indicating that the local void model may provide a consistent solution to the cluster number counts and Hubble parameter discrepancies.

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.

The power spectrum of the cosmic microwave background from both the Planck and WMAP data exhibits a slight dip in for multipoles in the range of l=10-30. We show that such a dip could be the result of resonant creation of a massive particle that couples to the inflaton field. For our best-fit models, epochs of resonant particle creation reenters the horizon at wave numbers of k* ~ 0.00011 (h/Mpc). The amplitude and location of these features correspond to the creation of a number of degenerate fermion species of mass ~ 15 times the planck mass during inflation with a coupling constant between the inflaton field and the created fermion species of near unity. Although the evidence is marginal, if this interpretation is correct, this could be one of the first observational hints of new physics at the Planck scale.

The vector mode of cosmological perturbation theory imprints characteristic signals on the weak lensing signals such as curl- and B-modes which are never imprinted by the scalar mode. However, the vector mode is neglected in the standard first-order cosmological perturbation theory since it only has a decaying mode. This situation changes if the cosmological perturbation theory is expanded up to second order. The second-order vector and tensor modes are inevitably induced by the product of the first-order scalar modes. We study the effect of the second-order vector mode on the weak lensing curl- and B-modes. We find that the curl-mode induced by the second-order vector mode is comparable to that induced by the primordial gravitational waves with the tensor-to-scalar ratio $r = 0.1$ at $\ell \approx 200$. In this case, the curl-mode induced by the second-order vector mode dominates at $\ell > 200$. Furthermore, the B-mode cosmic shear induced by the second-order vector mode dominates on almost all scales. However, we find that the observational signatures of the second-order vector and tensor modes cannot exceed the expected noise of ongoing and upcoming weak lensing measurements. We conclude that the curl- and B-modes induced by the second-order vector and tensor modes are unlikely to be detected in future experiments.

We summarize three recent efforts to constrain the first few moments of cosmic creation before and during the epoch of inflation. We consider two means to explain a slight dip in the power spectrum of the cosmic microwave background for multipoles in the range of $\ell= 10-30$ from both the {\it Planck} and {\it WMAP} data. We show that such a dip could be the result of resonant creation of a massive particle that couples to the inflaton field. For best-fit models, the epoch of resonant particle creation reenters the horizon at wave numbers of $k_* \sim 0.00011 \pm 0.0004 $ ($h$ Mpc$^{-1}$). The amplitude and location of these features correspond to the creation of a number of degenerate fermion species of mass $\sim 15/\lambda^{3/2} $ $m_{pl}$ during inflation where $\lambda$ is the coupling constant between the inflaton field and the created fermion species. Alternatively, one can explain the existence of such a dip as due to a jump in the inflation generating potential. We show that such a jump can also resolve the excessively large dark flow predicted from the M-theory landscape. Finally, we summarize our efforts to quantify constraints on the cosmic dark flow from a new analysis of the Type Ia supernova distance-redshift relation.

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.

The redshifted 21cm line signal from neutral hydrogens is a promising tool to probe the cosmic dawn and the epoch of reionization (EoR). Ongoing and future low-frequency radio experiments are expected to detect its fluctuations, especially through the power spectrum. In this paper, we give a physical interpretation of the time evolution of the power spectrum of the 21cm brightness temperature fluctuations, which can be decomposed into dark matter density, spin temperature and neutral fraction of hydrogen fluctuations. From the one-point statistics of the fluctuations, such as variance and skewness, we find that the peaks and dips in the time evolution are deeply related to X-ray heating of the intergalactic gas, which controls the spin temperature. We suggest the skewness of the brightness temperature distribution is a key observable to identify the onset of X-ray heating.

Concerted effort is currently ongoing to open up the Epoch of Reionization (EoR) ($z\sim$15-6) for studies with IR and radio telescopes. Whereas IR detections have been made of sources (Lyman-$\alpha$ emitters, quasars and drop-outs) in this redshift regime in relatively small fields of view, no direct detection of neutral hydrogen, via the redshifted 21-cm line, has yet been established. Such a direct detection is expected in the coming years, with ongoing surveys, and could open up the entire universe from $z\sim$6-200 for astrophysical and cosmological studies, opening not only the EoR, but also its preceding Cosmic Dawn ($z\sim$30-15) and possibly even the later phases of the Dark Ages ($z\sim$200-30). All currently ongoing experiments attempt statistical detections of the 21-cm signal during the EoR, with limited signal-to-noise. Direct imaging, except maybe on the largest (degree) scales at lower redshifts, as well as higher redshifts will remain out of reach. The Square Kilometre Array(SKA) will revolutionize the field, allowing direct imaging of neutral hydrogen from scales of arc-minutes to degrees over most of the redshift range $z\sim$6-28 with SKA1-LOW, and possibly even higher redshifts with the SKA2-LOW. In this SKA will be unique, and in parallel provide enormous potential of synergy with other upcoming facilities (e.g. JWST). In this chapter we summarize the physics of 21-cm emission, the different phases the universe is thought to go through, and the observables that the SKA can probe, referring where needed to detailed chapters in this volume (Abridged).

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.

The SKA will build upon early detections of the EoR by precursor instruments, such as MWA, PAPER, and LOFAR, and planned instruments, such as HERA, to make the first high signal-to-noise measurements of fluctuations in the 21 cm brightness temperature from both reionization and the cosmic dawn. This will allow both imaging and statistical maps of the 21cm signal at redshifts z = 6 - 27 and constrain the underlying cosmology and evolution of the density field. This era includes nearly 60% of the (in principle) observable volume of the Universe and many more linear modes than the CMB, presenting an opportunity for SKA to usher in a new level of precision cosmology. This optimistic picture is complicated by the need to understand and remove the effect of astrophysics, so that systematics rather than statistics will limit constraints. This chapter describes the cosmological, as opposed to astrophysical, information available to SKA. Key areas for discussion include: cosmological parameters constraints using 21cm fluctuations as a tracer of the density field; lensing of the 21cm signal, constraints on heating via exotic physics such as decaying or annihilating dark matter; impact of fundamental physics such as non-Gaussianity or warm dark matter on the source population; and constraints on the bulk flows arising from the decoupling of baryons and photons at z = 1000. The chapter explores the path to separating cosmology from astrophysics, for example via velocity space distortions and separation in redshift. We discuss new opportunities for extracting cosmology made possible by the sensitivity of SKA Phase 1 and explores the advances achievable with SKA2.

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.

One working hypothesis on which analyses of cosmological data are based is the zero ensemble mean hypothesis, which is related to the statistical homogeneity of cosmological perturbations. This hypothesis, however, should be tested by observational data in the current era of precision cosmology. Herein, we test the hypothesis by analyzing recent, foreground-reduced cosmic microwave background (CMB) maps, combining the spherical harmonic coefficients of the masked CMB temperature anisotropies in such a way that the combined variables can be treated as statistically independent samples. We find evidence against the zero mean hypothesis in two particular ranges of multipoles, with significance levels of $2.5 \sigma$ and $3.1 \sigma$ in the multipole ranges of $\ell \approx 61$-$86$ and $213$-$256$, respectively, for both the Planck and Wilkinson Microwave Anisotropy Probe maps. The latter signal is consistent with our previous result found by using brute-force Monte-Carlo simulations. However, within the method employed in this paper we conclude that the zero mean hypothesis is consistent with the current CMB data on the basis of Stouffer's weighted Z statistics, which takes multiple testing into account.

There can arise ubiquitous ultra-light scalar fields in the Universe, such as the pseudo-Goldstone bosons from the spontaneous breaking of an approximate symmetry, which can make a partial contribution to the dark matter and affect the large scale structure of the Universe. While the properties of those ultra-light dark matter are heavily model dependent and can vary in a wide range, we develop a model-independent analysis to forecast the constraints on their mass and abundance using futuristic but realistic 21 cm observables as well as CMB fluctuations, including CMB lensing measurements. Avoiding the highly nonlinear regime, the 21 cm emission line spectra are most sensitive to the ultra-light dark matter with mass m ~10^{-26} eV for which the precision attainable on mass and abundance bounds can be of order of a few percent.

Although the cosmological paradigm based on cold dark matter and adiabatic, nearly scale-invariant primordial fluctuations is consistent with a wide variety of existing observations, it has yet to be sufficiently tested on scales smaller than those of massive galaxies, and various alternatives have been proposed that differ significantly in the consequent small-scale power spectrum (SSPS) of large-scale structure. Here we show that a powerful probe of the SSPS at $k\gtrsim 10$ Mpc$^{-1}$ can be provided by the 21 cm forest, that is, systems of narrow absorption lines due to intervening, cold neutral hydrogen in the spectra of high-redshift background radio sources in the cosmic reionization epoch. Such features are expected to be caused predominantly by collapsed gas in starless minihalos, whose mass function can be very sensitive to the SSPS. As specific examples, we consider the effects of neutrino mass, running spectral index (RSI) and warm dark matter (WDM) on the SSPS, and evaluate the expected distribution in optical depth of 21 cm absorbers out to different redshifts. Within the current constraints on quantities such as the sum of neutrino masses $\sum m_\nu$, running of the primordial spectral index $d n_s/d \ln k$ and WDM particle mass $m_{\rm WDM}$, the statistics of the 21 cm forest manifest observationally significant differences that become larger at higher redshifts. In particular, it may be possible to probe the range of $m_{\rm WDM} \gtrsim 10$ keV that may otherwise be inaccessible. Future observations of the 21 cm forest by the Square Kilometer Array may offer a unique and valuable probe of the SSPS, as long as radio sources such as quasars or Population III gamma-ray bursts with sufficient brightness and number exist at redshifts of $z \gtrsim$ 10 - 20, and the astrophysical effects of reionization and heating can be discriminated.

The B-mode polarization spectrum of the Cosmic Microwave Background (CMB) may be the smoking gun of not only the primordial tensor mode but also of the primordial vector mode. If there exist nonzero vector-mode metric perturbations in the early Universe, they are known to be supported by anisotropic stress fluctuations of free-streaming particles such as neutrinos, and to create characteristic signatures on both the CMB temperature, E-mode, and B-mode polarization anisotropies. We place constraints on the properties of the primordial vector mode characterized by the vector-to-scalar ratio $r_{v}$ and the spectral index $n_{v}$ of the vector-shear power spectrum, from the {\it Planck} and BICEP2 B-mode data. We find that, for scale-invariant initial spectra, the $\Lambda$CDM model including the vector mode fits the data better than the model including the tensor mode. The difference in $\chi^{2}$ between the vector and tensor models is $\Delta\chi^{2} = 3.294$, because, on large scales the vector mode generates smaller temperature fluctuations than the tensor mode, which is preferred for the data. In contrast, the tensor mode can fit the data set equally well if we allow a significantly blue-tilted spectrum. We find that the best-fitting tensor mode has a large blue tilt and leads to an indistinct reionization bump on larger angular scales. The slightly red-tilted vector mode supported by the current data set can also create ${\cal O}(10^{-22})$-Gauss magnetic fields at cosmological recombination. Our constraints should motivate research that considers models of the early Universe that involve the vector mode.

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.

Recent observations of magnetic fields in intergalactic void regions and in high redshift galaxies may indicate that large scale magnetic fields have a primordial origin. If primordial magnetic fields were present soon after the recombination epoch, they would have induced density fluctuations on the one hand and dissipated their energy into the primordial gas on the other, and thereby significantly alter the thermal history of the Universe. Here we consider both the effects and calculate the brightness temperature fluctuations of the 21cm line using simple Monte Carlo simulations. We find that the fluctuations of the 21cm line from the energy dissipation appear only on very small scales and those from the density fluctuations always dominate on observationally relevant angular scales.

Quintessence, a scalar field model, has been proposed to account for the acceleration of the Universe at present. We discuss how accurately quintessence models are discriminated by future cosmological surveys, which include experiments of CMB, galaxy clustering, weak lensing, and the type Ia SNe surveys, by making use of the conventional parameterized dark energy models. We can see clear differences between the thawing and the freezing quintessence models at more than $1\sigma$ ($2\sigma$) confidence level as long as the present equation of state for quintessence is away from $-1$ as $w_X \gtrsim -0.95 (-0.90)$. However, it is found to be difficult to probe the effective mass squared for the potential in thawing models, whose signs are different between the quadratic and the cosine-type potentials. This fact may require us to invent a new estimator to distinguish quintessence models beyond the thawing and the freezing ones.

The global positioning system (GPS) is composed of thirty one satellites having atomic clocks with $10^{-15}$ accuracy on board and enables one to calibrate the primary standard for frequency on the ground. Using the fact that oscillators on the ground have been successfully stabilized with high accuracy by receiving radio waves emitted from the GPS satellites, we set a constraint on the strain amplitude of the gravitational wave background $h_{\rm c}$. We find that the GPS has already placed a meaningful constraint, and the constraint on the continuous component of gravitational waves is given as $h_{\rm c}<4.8\times 10^{-12}(1/ f)$ at $10^{-2}\lesssim f \lesssim 10^{0}$ Hz, for stabilized oscillators with ${\Delta \nu}/{\nu}\simeq 10^{-12}$. Thanks to the advantage of the Doppler tracking method, seismic oscillations do not affect the current constraint. Constraints on $h_c$ in the same frequency range from the velocity measurements by the lunar explorers in the Apollo mission are also derived.

Despite the ever growing observational evidence for the existence of the large scale magnetic fields, their origin and the evolution are not fully understood. If the magnetic fields are of primordial origin, they result in the generation of the secondary matter density perturbations and the previous studies show that such density perturbations enhance the number of dark matter halos. We extend the conventional spherical collapse model by including the Lorentz force which has not been implemented in the previous analysis to study the evolution of density perturbations produced by primordial magnetic fields. The critical over-density $\delta_{\rm c}$ characterizing the halo mass function turns out to be a bigger value, $\delta_{\rm c}\simeq 1.78$, than the conventional one $\delta_{\rm c}\simeq 1.69$ for the perturbations evolved only by the gravitational force. The difference in $\delta_{\rm c}$ between our model and the fully matter dominated cosmological model is small at a low redshift and, hence, only the high mass tail of the mass function is affected by the magnetic fields. At a high redshift, on the other hand, the difference in $\delta_{\rm c}$ becomes large enough to suppress the halo abundance over a wide range of mass scales. The halo abundance is reduced for instance by as large a factor as $\sim10^5$ at $z=9$.

Primordial magnetic fields (PMFs), which were generated in the early Universe before recombination, affect the motion of plasma and then the cosmic microwave background and the matter power spectrum. We consider constraints on PMFs with a characteristic correlation length from the observations of the anisotropies of the cosmic microwave background and the matter power spectrum. The spectrum of PMFs is modeled with multi-lognormal distributions, rather than power-law distribution, and we derive constraints on the strength $|\mathbf{B}_k|$ at each wave number $k$ along with the standard cosmological parameters in the flat Universe and the foreground sources. We obtain upper bounds on the field strengths at $k=10^{-1}, 10^{-2}, 10^{-4}$, and $10^{-5}$ Mpc$^{-1}$ as 4.7 nG, 2.1 nG, 5.3 nG, and 10.9 nG ($2\sigma$ C.L.) respectively, while the field strength at $k=10^{-3} $Mpc$^{-1}$ turns out to have a finite value as $|\mathbf{B}_{k = 10^{-3}}| = 6.2 \pm 1.3 $ nG ($1\sigma$ C.L.). This finite value is attributed to the finite values of BB mode data at $\ell > 300$ obtained from the QUAD experiment. If we do not include the BB mode data, we obtain only the upper bound on $B_{k=10^{-3}}$.

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.

Lower bounds are derived on the amplitude B of intergalactic magnetic fields (IGMFs) in the region between Galaxy and the blazar Mrk 421, from constraints on the delayed GeV pair-echo flux that are emitted by secondary e^-e^+ produced in \gamma\gamma interactions between primary TeV gamma-rays and the cosmic infrared background. The distribution of galaxies mapped by the Sloan Digital Sky Survey shows that this region is dominated by a large intergalactic void. We utilize data from long-term, simultaneous GeV-TeV observations by the Fermi Large Area Telescope and the ARGO-YBJ experiment extending over 850 days. For an assumed value of B, we evaluate the daily GeV pair-echo flux expected from the TeV data, select the dates where this exceeds the Fermi 2-\sigma sensitivity, compute the probability that this flux is excluded by the Fermi data for each date, and then combine the probabilities using the inverse normal method. Consequently, we exclude B < 10^{-20.5} G for a field coherence length of 1 kpc at 4- \sigma level, as long as plasma instabilities are unimportant for cooling of the pair beam. This is much more significant than the 2-\sigma bounds we obtained previously from observations of Mrk 501, by virtue of more extensive data from the ARGO-YBJ, as well as improved statistical analysis. Compared with most other studies of IGMF bounds, the evidence we present here for a non-zero IGMF is more robust as it does not rely on unproven assumptions on the primary TeV emission during unobserved periods.

Anisotropic stress perturbations induced by primordial Gaussian vector fields create non-Gaussianity in curvature perturbations. We found that such non-Gaussianity closely resembles the local-type non-Gaussianity parametrized by $f_{\rm NL}$, and generates scale-dependent bias in large-scale structures. We also found a simple relationship between the scale-dependent bias and the power spectrum of the vector fields. When the vector fields are interpreted as primordial magnetic fields, the effective $f_{\rm NL}$ is shown to be always negative. The scale-dependent bias provides a new approach to probing primordial vector fields.

Fast Independent Component Analysis (FastICA) is a component separation algorithm based on the levels of non-Gaussianity. Here we apply the FastICA to the component separation problem of the microwave background including carbon monoxide (CO) line emissions that are found to contaminate the PLANCK High Frequency Instrument (HFI) data. Specifically we prepare 100GHz, 143GHz, and 217GHz mock microwave sky maps including galactic thermal dust, NANTEN CO line, and the Cosmic Microwave Background (CMB) emissions, and then estimate the independent components based on the kurtosis. We find that the FastICA can successfully estimate the CO component as the first independent component in our deflection algorithm as its distribution has the largest degree of non-Gaussianity among the components. By subtracting the CO and the dust components from the original sky maps, we will be able to make an unbiased estimate of the cosmological CMB angular power spectrum.

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.

We study the spherical, top-hat collapse model for a mixed dark matter model including cold dark matter (CDM) and massive neutrinos of mass scales ranging from m_nu= 0.05 to a few 0.1eV, the range of lower- and upper-bounds implied from the neutrino oscillation experiments and the cosmological constraints. To develop this model, we properly take into account relative differences between the density perturbation amplitudes of different components (radiation, baryon, CDM and neutrinos) around the top-hat CDM overdensity region assuming the adiabatic initial conditions. Furthermore, we solve the linearized Boltzmann hierarchy equations to obtain time evolution of the lineariezed neutrino perturbations, yet including the effect of nonlinear gravitational potential due to the nonlinear CDM and baryon overdensities in the late stage. We find that the presence of massive neutrinos slows down the collapse of CDM (plus baryon) overdensity, however, that the neutrinos cannot fully catch up with the the nonlinear CDM perturbation due to its large free-streaming velocity for the ranges of neutrino masses and halo masses we consider. We find that, just like CDM models, the collapse time of CDM overdensity is well monitored by the linear-theory extrapolated overdensity of CDM plus baryon perturbation, smoothed with a given halo mass scale, if taking into account the suppression effect of the massive neutrinos on the linear growth rate. Using these findings, we argue that the presence of massive neutrinos of mass scales 0.05 or 0.1eV may cause a significant decrease in the abundance of massive halos compared to the model without the massive neutrinos; e.g., by 25% or factor 2, respectively, for halos with 10^15Ms and at z=1.

Magnetic fields appear wherever plasma and currents can be found. As such, they thread through all scales in Nature. It is natural, therefore, to suppose that magnetic fields might have been formed within the high temperature environments of the big bang. Such a primordial magnetic field (PMF) would be expected to arise from and/or influence a variety of cosmological phenomena such as inflation, cosmic phase transitions, big bang nucleosynthesis, the cosmic microwave background (CMB) temperature and polarization anisotropies, the cosmic gravity wave background, and the formation of large-scale structure. In this review, we summarize the development of theoretical models for analyzing the observational consequences of a PMF. We also summarize the current state of the art in the search for observational evidence of a PMF. In particular we review the framework needed to calculate the effects of a PMF power spectrum on the CMB and the development of large scale structure. We summarize the current constraints on the PMF amplitude $B_\lambda$ and the power spectral index $n_B$ and discuss prospects for better determining these quantities in the near future.

We investigate bounds on the strength of the primordial magnetic field (PMF) from the cosmic microwave background (CMB) bispectra of the intensity (temperature) modes induced from the auto- and cross-correlated bispectra of the scalar and tensor components of the PMF anisotropic stress. At first, we construct a general formula for the CMB intensity and polarization bispectra from PMFs composed of any type of perturbation. Then we derive an approximate expression which traces the exact shape of the CMB bispectrum in order to reduce the computation time with respect to a large number of the multipole configurations, and also show that the non-Gaussian structure coming from PMFs is classified as the local-type configuration. Computing the signal-to-noise ratio on the basis of the approximate formula with the information of the instrumental noises and resolutions, we find expected upper bounds on the magnetic field strength, when the magnetic spectrum is nearly scale invariant ($n_B = -2.9$), smoothed on $1 {\rm Mpc}$ scale at 95% confidence level from the WMAP and PLANCK experiments as $B_{1 \rm Mpc} < 4.0 - 6.7 {\rm nG}$ and $3.8 - 6.5 {\rm nG}$, respectively, depending on the energy scale of the magnetic field production from $10^{14} {\rm GeV}$ to $10^3 {\rm GeV}$. Our new consequences imply slight overestimations by the previous rough discussions.

The standard models of inflation predict statistically homogeneous and isotropic primordial fluc- tuations, which should be tested by observations. In this paper we test the statistical isotropy of the mean of the CMB temperature fluctuations measured by the WMAP in the spherical harmonic space. A classical method to test a mean, like the simple student's-t test, is not appropriate for this purpose because the WMAP data contain anisotropic instrumental noise and suffer from the effect of the mask for the foreground emissions which breaks the statistical independence. Here we perform a band-power analysis with Monte-Carlo simulations in which we take into account the anisotropic noise and the mask. We find evidence of a non-zero mean at 99.93 % confidence level in a particular range of multipoles. The evidence against the zero-mean assumption as a whole is still significant at the 99 % confidence level even if the fact is taken into account that we have tested multiple ranges.

The effect of primordial magnetic fields on X-ray or S-Z galaxy cluster survey is investigated. After recombination, the primordial magnetic fields generate additional density fluctuations. Such density fluctuations enhance the number of galaxy clusters. Taking into account the density fluctuations generated by primordial magnetic fields, we calculate the number of galaxy clusters based on the Press-Schechter formalism. Comparing with the results of Chandra X-ray galaxy cluster survey, we found that the existence of primordial magnetic fields with amplitude larger than 10 nGuass would be inconsistent. Moreover, we show that S-Z cluster surveys also have a sensitivity to constrain primordial magnetic fields. Especially SPT S-Z cluster survey has a potential to constrain the primordial magnetic fields with several nano Gauss.

The primordial non-Gaussianity of local type affects the clustering of dark matter halos, and the planned deep and wide photometric surveys are suitable for examining this class of non-Gaussianity. In our previous paper, we investigated the constraint from the cross correlation between CMB lensing potential and galaxy angular distribution on the primordial non-Gaussianity, without taking into account redshift slicing. To improve our previous analysis, in this paper, we add the galaxy lensing shear into our analysis and take into account redshift slicing to follow the redshift evolution of the clustering. By calculating 81 power spectra and using the Fisher matrix method, we find that the constraint on the primordial non-Gaussianity can be improved from {\Delta}fNL \sim 5.4 to 5.1 by including the galaxy-galaxy lensing shear cross correlations expected from the Hyper Suprime-Cam survey (HSC), in comparison with the constraint without any cross correlations. Moreover, the constraint can go down to {\Delta}fNL \sim 4.8 by including the galaxy-CMB lensing cross correlations from the ACTPol and Planck experiments.

We study the effect of primordial magnetic fields (PMFs) on the anisotropies of the cosmic microwave background (CMB). We assume the spectrum of PMFs is described by log-normal distribution which has a characteristic scale, rather than power-law spectrum. This scale is expected to reflect the generation mechanisms and our analysis is complementary to previous studies with power-law spectrum. We calculate power spectra of energy density and Lorentz force of the log-normal PMFs, and then calculate CMB temperature and polarization angular power spectra from scalar, vector, and tensor modes of perturbations generated from such PMFs. By comparing these spectra with WMAP7, QUaD, CBI, Boomerang, and ACBAR data sets, we find that the current CMB data set places the strongest constraint at $k\simeq 10^{-2.5}$ Mpc$^{-1}$ with the upper limit $B\lesssim 3$ nG.

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

Magnetic fields are everywhere in nature and they play an important role in every astronomical environment which involves the formation of plasma and currents. It is natural therefore to suppose that magnetic fields could be present in the turbulent high temperature environment of the big bang. Such a primordial magnetic field (PMF) would be expected to manifest itself in the cosmic microwave background (CMB) temperature and polarization anisotropies, and also in the formation of large- scale structure. In this review we summarize the theoretical framework which we have developed to calculate the PMF power spectrum to high precision. Using this formulation, we summarize calculations of the effects of a PMF which take accurate quantitative account of the time evolution of the cut off scale. We review the constructed numerical program, which is without approximation, and an improvement over the approach used in a number of previous works for studying the effect of the PMF on the cosmological perturbations. We demonstrate how the PMF is an important cosmological physical process on small scales. We also summarize the current constraints on the PMF amplitude $B_\lambda$ and the power spectral index $n_B$ which have been deduced from the available CMB observational data by using our computational framework.

In the pre-recombination era, cosmological density fluctuations can naturally generate magnetic fields through Thomson scatterings. In previous studies, only the magnetic field generation from the initially-adiabatic fluctuations has been considered. Here we investigate the generation of cosmological magnetic fields sourced by the primordial non-adiabatic fluctuations based on the cosmological perturbation theory, using the tight-coupling approximations between photon and baryon fluids. It is found that the magnetic fields from the non-adiabatic fluctuations can arise at the first-order expansion of the tight coupling approximation. This result is in contrast to the case of adiabatic initial fluctuations, where the magnetic fields can be generated only at the second-order. In a general case where the primordial density perturbations contain small non-adiabatic fluctuations on the top of the dominant adiabatic ones, we show that the leading source of magnetic fields is given by the second-order coupling of the adiabatic and non-adiabatic fluctuations. We calculate the power spectrum of the generated magnetic fields when the non-adiabatic fluctuations have a blue power spectrum, which has been suggested by recent cosmological observations.

We derive lower bounds on intergalactic magnetic fields (IGMFs) from upper limits on the pair echo emission from the blazar Mrk 501, that is, delayed GeV emission from secondary $e^{-}e^{+}$ pairs produced via interactions of primary TeV gamma rays with the cosmic infrared background. Utilizing only simultaneous GeV-TeV light curves observed by VERITAS, MAGIC and {\it Fermi}-LAT during a multiwavelength campaign in 2009 that included a TeV flare, bounds are deduced on the IGMF strength of $B \gtrsim 10^{-20} {\rm G}$ at 90% confidence level for a field coherence length of 1 kpc. Since our analysis is based firmly on the observational data alone and nearly free of assumptions concerning the primary TeV flux in unobserved periods or spectral bands, our evaluation of the pair echo flux is conservative and the evidence for a non-zero IGMF is more robust compared to previous studies.

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 present a detailed calculation of our previous short paper [M. Shiraishi, D. Nitta, S. Yokoyama, K. Ichiki, and K. Takahashi, Phys. Rev. D 82, 121302 (2010).] in which we have investigated a constraint on the magnetic field strength through comic microwave background temperature bispectrum of vector modes induced from primordial magnetic fields. By taking into account full angular dependence of the bispectrum with spin spherical harmonics and Wigner symbols, we explicitly show that the cosmic microwave background bispectrum induced from the statistical-isotropic primordial vector fluctuations can be also described as an angle-averaged form in the rotationally invariant way. We also study the cases with different spectral indices of the power spectrum of the primordial magnetic fields.

If the seed magnetic fields exist in the early Universe, tensor components of their anisotropic stresses are not compensated prior to neutrino decoupling and the tensor metric perturbations generated from them survive passively. Consequently, due to the decay of these metric perturbations after recombination, the so-called integrated Sachs-Wolfe effect, the large-scale fluctuations of CMB radiation are significantly boosted. This kind of CMB anisotropy is called the "tensor passive mode." Because these fluctuations deviate largely from the Gaussian statistics due to the quadratic dependence on the strength of the Gaussian magnetic field, not only the power spectrum but also the higher-order correlations have reasonable signals. With these motives, we compute the CMB bispectrum induced by this mode. When the magnetic spectrum obeys a nearly scale-invariant shape, we obtain an estimation of a typical value of the normalized reduced bispectrum as $\ell_1(\ell_1 + 1)\ell_3(\ell_3+1)|b_{\ell_1\ell_2\ell_3}| \sim (130-6) \times 10^{-16} (B_{1 \rm Mpc} / 4.7 {\rm nG})^6$ depending on the energy scale of the magnetic field production from $10^{14}$GeV to $10^3$GeV. Here, $B_{1 {\rm Mpc}}$ is the strength of the primordial magnetic field smoothed on $1 {\rm Mpc}$. From the above estimation and the current observational constraint on the primordial non-Gaussianity, we get a rough constraint on the magnetic field strength as $B_{1 {\rm Mpc}} < 2.6 - 4.4 {\rm nG}$.

We present an all-sky formalism for the Cosmic Microwave Background (CMB) bispectrum induced by the primordial non-Gaussianities not only in scalar but also in vector and tensor fluctuations. We find that the bispectrum can be formed in an explicitly rationally invariant way by taking into account the angular and polarization dependences of the vector and tensor modes. To demonstrate this and present how to use our formalism, we consider a specific example of the correlation between two scalars and a graviton as the source of non-Gaussianity. As a result, we show that the CMB reduced bispectrum of the intensity anisotropies is evaluated as a function of the multipole and the coupling constant between two scalars and a graviton denoted by $g_{tss}$; $|b_{\ell \ell \ell}| \sim \ell^{-4} \times 8 \times 10^{-18} |g_{tss}|$. By estimating the signal-to-noise ratio, we find that the constraint as $|g_{tss}| < 6$ will be expected from the PLANCK experiment.

We calculate the cosmic microwave background bispectrum of vector modes induced from primordial magnetic fields. We take into account the full angular dependence of the bispectrum and discuss the amplitude and also the shape of the bispectrum in $\ell$ space. In the squeezed limit as $\ell_1 = \ell_2 \gg \ell_3$, we estimate a typical values of the normalized reduced bispectrum as $\ell_1(\ell_1 + 1)\ell_3(\ell_3+1)|b_{\ell_1\ell_2\ell_3}| \sim 2 \times 10^{-19}$, for the strength of the primordial magnetic field smoothed on $1 {\rm Mpc}$ scale $B_{1 \rm Mpc} = 4.7 \rm nG$ assuming a nearly scale-invariant spectrum of magnetic fields. We find that a constraint on the magnetic field strength will be roughly estimated as $B_{1 \rm Mpc} < 10 {\rm nG}$ if PLANCK places a limit on the nonlinearity parameter of local-type configuration as $|f^{\rm local}_{\rm NL}| < 5$.

To estimate cosmological parameters from a given dataset, we need to construct a likelihood function, which sometimes has a complicated functional form. We introduce the copula, a mathematical tool to construct an arbitrary multivariate distribution function from one-dimensional marginal distribution functions with any given dependence structure. It is shown that a likelihood function constructed by the so-called Gaussian copula can reproduce very well the n-dimensional probability distribution of the cosmic shear power spectrum obtained from a large number of ray-tracing simulations. This suggests that the Copula likelihood will be a powerful tool for future weak lensing analyses, instead of the conventional multivariate Gaussian likelihood.

The estimation of cosmological parameters from a given data set requires a construction of a likelihood function which, in general, has a complicated functional form. We adopt a Gaussian copula and constructed a copula likelihood function for the convergence power spectrum from a weak lensing survey. We show that the parameter estimation based on the Gaussian likelihood erroneously introduces a systematic shift in the confidence region, in particular for a parameter of the dark energy equation of state w. Thus, the copula likelihood should be used in future cosmological observations.

Recent studies have shown that the primordial non-Gaussianity affects clustering of dark matter halos through a scale-dependent bias and various constraints on the non-Gaussianity through this scale-dependent bias have been placed. Here we introduce the cross-correlation between the CMB lensing potential and the galaxy angular distribution to effectively extract information about the bias from the galaxy distribution. Then, we estimate the error of non-linear parameter, f_NL, for the on-going CMB experiments and galaxy surveys, such as Planck and Hyper Suprime-Cam (HSC). We found that for the constraint on f_NL with Planck and HSC, the wide field galaxy survey is preferable to the deep one, and the expected error on f_NL can be as small as: \Delta f_NL ~ 20 for b_0 = 2 and \Delta f_NL ~ 10 for b_0 = 4, where b_0 is the linear bias parameter. It is also found that future wide field galaxy survey could achieve \Delta f_NL ~ 5 with CMB prior from Planck if one could observe highly biased objects at higher redshift (z ~ 2).

We present a complete set of formulae for calculating the bispectra of CMB temperature and polarization anisotropies generated from non-Gaussianity in the vector and tensor mode perturbations. In the all-sky analysis it is found that the bispectrum formulae for the tensor and vector-mode non-Gaussianity formally take complicated forms compared to the scalar mode one because the photon transfer functions in the tensor and vector modes depend on the azimuthal angle between the direction of the wave number vector of the photon's perturbation and that of the line of sight. We demonstrate that flat-sky approximations remove this difficulty because this kind of azimuthal angle dependence apparently vanishes in the flat-sky limit. Through the flat-sky analysis, we also find that the vector or tensor bispectrum of $B$-mode polarization vanishes in the squeezed limit, unless the cosmological parity is violated at the nonlinear level.

We discuss the expected properties of pair echo emission from gamma-ray bursts (GRBs) at high redshifts ($z \gtrsim 5$), their detectability, and the consequent implications for probing intergalactic magnetic fields (IGMFs) at early epochs. Pair echos comprise inverse Compton emission by secondary electron-positron pairs produced via interactions between primary gamma-rays from the GRB and low-energy photons of the diffuse intergalactic radiation, arriving with a time delay that depends on the nature of the intervening IGMFs. At sufficiently high $z$, the IGMFs are unlikely to have been significantly contaminated by astrophysical outflows, and the relevant intergalactic radiation may be dominated by the well-understood cosmic microwave background (CMB). Pair echoes from luminous GRBs at $z \sim 5-10$ may be detectable by future facilities such as the Cherenkov Telescope Array or the Advanced Gamma-ray Imaging System, as long as the GRB primary emission extends to multi-TeV energies, the comoving IGMFs at these redshifts are $B \sim 10^{-16}-10^{-15}~{\rm Gauss}$, and the non-CMB component of the diffuse intergalactic radiation is relatively low. Observations of pair echos from high-$z$ GRBs can provide a unique, in-situ probe of weak IGMFs during the epochs of early structure formation and cosmic reionization.

We have made an analysis of limits on the neutrino mass based upon the formation of large-scale structure in the presence of a primordial magnetic field. We find that a new upper bound on the neutrino mass is possible based upon fits to the cosmic microwave background and matter power spectrum when the existing independent constraints on the matter density fluctuation parameter $\sigma_8$ and the primordial magnetic field are taken into account.

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 the newest statistical and numerical analysis of the matter and cosmic microwave background power spectrum with effects of the primordial magnetic field (PMF) included. New limits to the PMF strength and power spectral index are obtained based upon the accumulated data for both the matter and CMB power spectra on small angular scales. We find that a maximum develops in the probability distribution for a magnitude of the PMF of $|B_\lambda| = 0.85 \pm 1.25(\pm 1\sigma)$ nG on a comoving scale of at 1 Mpc, corresponding to upper limits of $\mathbf{< 2.10 nG} (68% \mathrm{CL})$ and $\mathbf{< 2.98 nG} (95% \mathrm{CL})$. While for the power spectral index we find $n_\mathrm{B}= -2.37^{+0.88}_{-0.73}(\pm 1\sigma)$, corresponding to upper limits of $ \mathbf{< -1.19} (68% \mathrm{CL})$ and $ \mathbf{< -0.25} (95% \mathrm{CL})$. This result provides new constraints on models for magnetic field generation and the physics of the early universe. We conclude that future observational programs for the CMB and matter power spectrum will likely provide not only upper limits but also lower limits to the PMF parameters.

As a test of the standard inflationary cosmology, which generically predicts nearly scale-invariant spectrum of primordial curvature fluctuations, we perform Markov-Chain Monte-Carlo analysis to search for possible modulations in the power spectrum and determine its shape together with the cosmological parameters using cosmic microwave background radiation data. By incorporating various three-parameter features on the simple power-law spectrum, we find an oscillatory modulation localized around the comoving wavenumber $k\simeq 0.009\mathrm{Mpc}^{-1}$ at 99.995% confidence level which improves the log-likelihood as much as $-\Delta 2\ln {\cal L}\equiv \Delta \chi^2_{\rm eff}=-22$. This feature can be detected even if we use only the cross correlation between the temperature and the E-mode polarization anisotropies.

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

Possible existence of the primordial magnetic fields has affected the structure formation of the universe. In this paper it is shown that the initial conditions for density perturbations with magnetic fields derived in previous works are inconsistent with Einstein equations. We find that this inconsistency arises due to the unwanted cancellation of contributions from the magnetic fields and primordial radiations. A complete set of equations and consistent initial conditions in the long wavelength limit are given with an explicit derivation in the covariant approach with CDM frame, by newly taking into account a non-relativistic matter contribution in the radiation dominated era. By solving these equations numerically, we derive the angular spectrum of cosmic microwave background anisotropies and the matter power spectrum with magnetic fields. We find that the amplitude of the angular power spectrum of CMB anisotropies can alter at most a order of magnitude at $l \lesssim 4000$ compared with the previous results in the literature.

Recent observations of the cosmic microwave background (CMB) have extended the measured power spectrum to higher multipoles $l\gtrsim$1000, and there appears to be possible evidence for excess power on small angular scales. The primordial magnetic field (PMF) can strongly affect the CMB power spectrum and the formation of large scale structure. In this paper, we calculate the CMB temperature anisotropies generated by including a power-law magnetic field at the photon last-scattering surface (PLSS). We then deduce an upper limit on the PMF based on our theoretical analysis of the power excess on small angular scales. We have taken into account several important effects such as the modified matter sound speed in the presence of a magnetic field. An upper limit to the field strength of $|B_\lambda|\lesssim$ 4.7 nG at the present scale of 1 Mpc is deduced. This is obtained by comparing the calculated theoretical result including the Sunyaev-Zeldovich (SZ) effect with recent observed data on the small-scale CMB anisotropies from the $Wilkinson Microwave Anisotropy Probe$ (WMAP), the Cosmic Background Imager (CBI), and the Arcminute Cosmology Bolometer Array Receiver (ACBAR). We discuss several possible mechanisms for the generation and evolution of the PMF.

The weak lensing (WL) distortions of distant galaxy images are sensitive to neutrino masses by probing the suppression effect on clustering strengths of total matter in large-scale structure. We use the latest measurement of WL correlations, the CFHTLS data, to explore constraints on neutrino masses. We find that, while the WL data alone cannot place a stringent limit on neutrino masses due to parameter degeneracies, the constraint can be significantly improved when combined with other cosmological probes, the WMAP 5-year (WMAP5) data and the distance measurements of type-Ia supernovae (SNe) and baryon acoustic oscillations (BAO). The upper bounds on the sum of neutrino masses are m_tot = 1.1, 0.76 and 0.54 eV (95% CL) for WL+WMAP5, WMAP5+SNe+BAO, and WL+WMAP5+SNe+BAO, respectively, assuming a flat LCDM model with finite-mass neutrinos. In deriving these constraints, our analysis includes the non-Gaussian covariances of the WL correlation functions to properly take into account significant correlations between different angles.

A primordial magnetic field (PMF) can affect the evolution of density field fluctuations in the early universe.In this paper we constrain the PMF amplitude $B_\lambda$ and power spectral index $n_\mathrm{B}$ by comparing calculated density field fluctuations with observational data, i.e. the number density fluctuation of galaxies.We show that the observational constraints on cosmological density fluctuations, as parameterized by $\sigma_8$, lead to strong constraints on the amplitude and spectral index of the PMF.

More than a dozen blazars are known to be emitters of multi-TeV gamma rays, often with strong and rapid flaring activity. By interacting with photons of the cosmic microwave and infrared backgrounds, these gamma rays inevitably produce electron-positron pairs, which in turn radiate secondary inverse Compton gamma rays in the GeV-TeV range with a characteristic time delay that depends on the properties of the intergalactic magnetic field (IGMF). For sufficiently weak IGMF, such "pair echo" emission may be detectable by the Gamma-ray Large Area Space Telescope (GLAST), providing valuable information on the IGMF. We perform detailed calculations of the time-dependent spectra of pair echos from flaring TeV blazars such as Mrk 501 and PKS 2155-304, taking proper account of the echo geometry and other crucial effects. In some cases, the presence of a weak but non-zero IGMF may enhance the detectability of echos. We discuss the quantitative constraints that can be imposed on the IGMF from GLAST observations, including the case of non-detections.

High-energy emission from gamma-ray bursts (GRBs) can give rise to pair echos, i.e. delayed inverse Compton emission from secondary $e^{\pm}$ pairs produced in $\gamma-\gamma$ interactions with intergalactic background radiation. We investigate the detectability of such emission with modern-day gamma-ray telescopes. The spectra and light curves are calculated for a wide range of parameters, applying the formalism recently developed by Ichiki et al. The flux depends strongly on the unknown magnitude and coherence length of intergalactic magnetic fields, and we delineate the range of field strength and redshift that allow detectable echos. Relevant uncertainties such as the high-energy cutoff of the primary gamma-ray spectrum and the intensity of the cosmic infrared background are addressed. GLAST and MAGIC may be able to detect pair echo emission from GRBs with redshift $\lesssim 1$ if the primary spectra extend to $\sim 10 ~ {\rm TeV}$.

We investigate whether interaction between massive neutrinos and quintessence scalar field is the origin of the late time accelerated expansion of the universe. We present explicit formulas of the cosmological linear perturbation theory in the neutrinos probes of dark-energy model, and calculate cosmic microwave background anisotropies and matter power spectra. In these models, the evolution of the mass of neutrinos is determined by the quintessence scalar field, which is responsible for a varying effective equation of states: $\omega_{eff}(z)$ goes down -1. We consider several types of scalar field potential and put constraints on the coupling parameter between neutrinos and dark energy. By combining data from cosmic microwave background (CMB) experiments including the WMAP 3-year results, large scale structure with 2dFGRS data sets, we constrain the hypothesis of massive neutrinos in the mass-varying neutrino scenario. Assuming the flatness of the universe, the constraint we can derive from the current observation is $\sum m_{\nu} < 0.45$ eV at 1$\sigma$ (0.87 eV at 2$\sigma$) confidence level for the sum over three species of neutrinos. The dynamics of scalar field and the impact of scalar field perturbations on cosmic microwave background anisotropies are discussed. We also discuss on the instability issue of the our model and confirm that neutrinos are stable against the density fluctuation.

If a primordial magnetic field (PMF) is present during photon decoupling and afterward, a finite neutrino mass can affect all modes of the CMB. In this work, we expand on earlier studies of the scalar mode effects by constructing the vector and tensor mode equations in the presence of massive neutrinos and a PMF. We compute the power spectrum of the various modes in an illustrative example and find that the neutrino mass can significantly affect the vector and tensor modes when a PMF exists, while the effects are negligible for no PMF. The most prominent result of the present analysis is the behavior of the EE component of the tensor mode at low multipoles. For massive neutrinos the EE mode can become comparable to the observed primary anisotropy. Therefore, if and when the EE mode power spectrum is measured at low multipoles the possibility exists to place a strong constraint on the sum of the neutrino masses.

We constrain the thermal evolution of the universe with a decaying cosmological term by using the method of the analysis for the Wilkinson Microwave Anisotropy Probe (WMAP) observation data. The cosmological term is assumed to be a function of the scale factor that increases toward the early universe, and the radiation energy density is lower compared to that in the model with the standard cosmological "constant" (LCDM). The decrease in the radiation density affects the thermal history of the universe; e.g. the photon decoupling occurs at higher-z compared to the case of the standard LCDM model. As a consequence, a decaying cosmological term affects the cosmic microwave background anisotropy. Thanks to the Markov-chain Monte Carlo method, we compare the angular power spectrum in the decaying LCDM model with the CMB data, and we get severe constraints on parameters of the model.

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.

We investigate the delayed, secondary GeV-TeV emission of gamma-ray bursts and its potential to probe the nature of intergalactic magnetic fields. Geometrical effects are properly taken into account for the time delay between primary high energy photons and secondary inverse Compton photons from electron-positron pairs, which are produced in $\gamma$-$\gamma$ interactions with background radiation fields and deflected by intervening magnetic fields. The time-dependent spectra of the delayed emission are evaluated for a wide range of magnetic field strengths and redshifts. The typical flux and delay time of secondary photons from bursts at $z \sim 1$ are respectively $\sim 10^{-8}$ GeV cm$^{-2}$ s$^{-1}$ and $\sim 10^4$ s if the field strengths are $\sim 10^{-18}$ G, as might be the case in intergalactic void regions. We find crucial differences between the cases of coherent and tangled magnetic fields, as well as dependences on the field coherence length.

We investigate whether interaction between massive neutrinos and quintessence scalar field is the origin of the late time accelerated expansion of the universe. We present cosmological perturbation theory in neutrinos probe interacting dark-energy models, and calculate cosmic microwave background anisotropies and matter power spectrum. In these models, the evolution of the mass of neutrinos is determined by the quintessence scalar field, which is responsible for the cosmic acceleration today. We consider several types of scalar field potentials and put constraints on the coupling parameter between neutrinos and dark energy. Assuming the flatness of the universe, the constraint we can derive from the current observation is $\sum m_{\nu} < 0.87 eV$ at the 95 % confidence level for the sum over three species of neutrinos. We also discuss on the stability issue of the our model and on the impact of the scattering term in Boltzmann equation from the mass-varying neutrinos.

We present cosmological perturbation theory in neutrinos probe interacting dark-energy models, and calculate cosmic microwave background anisotropies and matter power spectrum. In these models, the evolution of the mass of neutrinos is determined by the quintessence scalar field, which is responsible for the cosmic acceleration today. We consider several types of scalar field potentials and put constraints on the coupling parameter between neutrinos and dark energy. Assuming the flatness of the universe, the constraint we can derive from the current observation is $\sum m_{\nu} < 0.87 eV$ at the 95 % confidence level for the sum over three species of neutrinos. We also discuss on the stability issue of the our model and on the impact of the scattering term in Boltzmann equation from the mass-varying neutrinos.

Weak gravitational lensing by foreground density perturbations generates a gradient mode in the shear of background images. In contrast, cosmological tensor perturbations induce a non-zero curl mode associated with image rotations. In this note, we study the lensing signatures of both primordial gravitational waves from inflation and second-order gravitational waves generated from the observed spectrum of primordial density fluctuations. We derive the curl mode for galaxy lensing surveys at redshifts of 1 to 3 and for lensing of the cosmic microwave background (CMB) at a redshift of 1100. We find that the curl mode angular power spectrum associated with secondary tensor modes for galaxy lensing surveys dominates over the corresponding signal generated by primary gravitational waves from inflation. However, both tensor contributions to the shear curl mode spectrum are below the projected noise levels of upcoming galaxy and CMB lensing surveys and therefore are unlikely to be detectable.

The existence of a primordial magnetic field (PMF) would affect both the temperature and polarization anisotropies of the cosmic microwave background (CMB). It also provides a plausible explanation for the possible disparity between observations and theoretical fits to the CMB power spectrum. Here we report on calculations of not only the numerical CMB power spectrum from the PMF, but also the correlations between the CMB power spectrum from the PMF and the primary curvature perturbations. We then deduce a precise estimate of the PMF effect on all modes of perturbations. We find that the PMF affects not only the CMB TT and TE modes on small angular scales, but also on large angular scales. The introduction of a PMF leads to a better fit to the CMB power spectrum for the higher multipoles, and the fit at lowest multipoles can be used to constrain the correlation of the PMF with the density fluctuations for large negative values of the spectral index. Our prediction for the BB mode for a PMF average field strength $|B_\lambda| =4.0$ nG is consistent with the upper limit on the BB mode deduced from the latest CMB observations. We find that the BB mode is dominated by the vector mode of the PMF for higher multipoles. We also show that by fitting the complete power spectrum one can break the degeneracy between the PMF amplitude and its power spectral index.

We discuss the polarization signature of primordial gravitational waves imprinted in cosmic microwave background (CMB) anisotropies. The high-energy physics motivated by superstring theory or M-theory generically yield parity violating terms, which may produce a circularly polarized gravitational wave background (GWB) during inflation. In contrast to the standard prediction of inflation with un-polarized GWB, circularly polarized GWB generates non-vanishing TB and EB-mode power spectra of CMB anisotropies. We evaluate the TB and EB-mode power spectra taking into account the secondary effects and investigate the dependence of cosmological parameters. We then discuss current constraints on the circularly polarized GWB from large angular scales (l < 16) of the three year WMAP data. Prospects for future CMB experiments are also investigated based on a Monte Carlo analysis of parameter estimation, showing that the circular polarization degree, varepsilon, which is the asymmetry of the tensor power spectra between right- and left-handed modes normalized by the total amplitude, can be measured down to |varepsilon| 0.35(r/0.05)^{-0.6}.

Effects of a generalized dark energy fluid is investigated on cosmic density fluctuations such as cosmic microwave background. As a general dark energy fluid, we take into consideration the possibility of the anisotropic stress for dark energy, which has not been discussed much in the literature. We comprehensively study its effects on the evolution of density fluctuations along with that of non-adiabatic pressure fluctuation of dark energy, then give constraints on such a generalized dark energy from current observations. We show that, though we cannot find any stringent limits on the anisotropic stress or the non-adiabatic pressure fluctuation themselves, the constrains on the equation of state of dark energy can be affected in some cases by the nature of dark energy fluctuation characterized by these properties. This may have important implications to the strategy to study the nature of dark energy.

We derive the complete spectrum of gravitational waves induced by primordial scalar perturbations ranging over all observable wavelengths. This scalar-induced contribution can be computed directly from the observed scalar perturbations and general relativity and is, in this sense, independent of the cosmological model for generating the perturbations. The spectrum is scale-invariant on small scales, but has an interesting scale-dependence on large and intermediate scales, where scalar-induced gravitational waves do not redshift and are hence enhanced relative to the background density of the Universe. This contribution to the tensor spectrum is significantly different in form from the direct model-dependent primordial tensor spectrum and, although small in magnitude, it dominates the primordial signal for some cosmological models. We confirm our analytical results by direct numerical integration of the equations of motion.

We investigate the evolution of a universe with a decaying cosmological term (vacuum energy) that is assumed to be a function of the scale factor. In this model, while the cosmological term increases to the early universe, the radiation energy density is lower than the model with the cosmological "constant". We find that the effects of the decaying cosmological term on the expansion rate at the redshift z<2 is negligible. However, the decrease in the radiation density affects on the thermal history of the universe; e.g. the photon decoupling occurs at higher $z$ compared to the case of the standard \Lambda CDM model. As a consequence, a decaying cosmological term affects on the cosmic microwave background anisotropy. We show the angular power spectrum in D\Lambda CDM model and compare with the Wilkinson Microwave Anisotropy Probe (WMAP) data.

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.

Magnetic field plays a very important role in many astronomical phenomena at various scales of the universe. It is no exception in the early universe. Since the energy density, pressure, and tension of the primordial magnetic field affect gravitational collapses of plasma, the formation of seeds for large scale structures should be influenced by them. Here we numerically investigate the effects of stochastic primordial magnetic field on the seeds of large scale structures in the universe in detail. We found that the amplitude ratio between the density spectra with and without PMF ($|P(k)/P_0(k)|$ at $k>0.2$ Mpc$^{-1}$) lies between 75% and 130% at present for the range of PMF strengths 0.5 nG $< B_\lambda < 1.0$ nG, depending on the spectral index of PMF and the correlation between the matter density and the PMF distributions.

Effects of neutrino free streaming is evaluated on the primordial spectrum of gravitational radiation taking both neutrino chemical potential and masses into account. The former or the lepton asymmetry induces two competitive effects, namely, to increase anisotropic pressure, which damps the gravitational wave more, and to delay the matter-radiation equality time, which reduces the damping. The latter effect is more prominent and a large lepton asymmetry would reduce the damping. We may thereby be able to measure the magnitude of lepton asymmetry from the primordial gravitational wave spectrum.

The primordial magnetic field (PMF) can strongly affect the cosmic microwave background (CMB) power spectrum and the formation of large scale structure. In this presentation, we calculate the CMB temperature anisotropies generated by including a power-law magnetic field at the photon last scattering surface (PLSS). We then deduce an upper limit on the primordial magnetic field based upon our theoretical analysis of the power excess on small angular scales. We have taken into account several important effects such as the modified matter sound speed in the presence of a magnetic field. An upper limit to the field strength of $|B_\lambda|\lesssim$ 4.7 nG at the present scale of 1 Mpc is deduced. This is obtained by comparing the calculated theoretical result including the Sunyaev-Zeldovich (SZ) effect with recent observed data on the small scale CMB anisotropies from the Wilkinson Microwave Anisotropy Probe (WMAP), the Cosmic Background Imager (CBI) and the Arcminute Cosmology Bolometer Array Receiver (ACBAR). We discuss several possible mechanisms for the generation and evolution of the PMF before, during and after the BBN.

We explore observational constraints on a cosmological brane-world scenario in which the bulk is not empty. Rather, exchange of mass-energy between the bulk and the bane is allowed. The evolution of matter fields to an observer on the brane is then modified due to new terms in the energy momentum tensor describing this exchange. We investigate the constraints from various cosmological observations on the flow of matter from the bulk into the brane. Interestingly, we show that it is possible to have a $\Lambda = 0$ cosmology to an observer in the brane which satisfies standard cosmological constraints including the CMB temperature fluctuations, Type Ia supernovae at high redshift, and the matter power spectrum. This model even accounts for the observed suppression of the CMB power spectrum at low multipoles. In this cosmology, the observed cosmic acceleration is attributable to the flow of matter from the bulk to the brane. A peculiar aspect of this cosmology is that the present dark-matter content of the universe may be significantly larger than that of a $\Lambda$CDM cosmology. Its influence, however, is offset by the dark-radiation term. Possible additional observational tests of this new cosmological paradigm are suggested.

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 studied the generation of magnetic fields by the Biermann mechanism in the pair-instability supernovae explosions of the first stars. The Biermann mechanism produces magnetic fields in the shocked region between the bubble and interstellar medium (ISM), even if magnetic fields are absent initially. We have performed a series of two-dimensional magnetohydrodynamic simulations with the Biermann term and estimate the amplitude and total energy of the produced magnetic fields. We find that magnetic fields with amplitude $10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of magnetic fields generated depend on specific values of initial conditions. This corresponds to magnetic fields of $10^{28}-10^{31}$ ergs per each supernova remnant, which is strong enough to be the seed magnetic field for a galactic and/or interstellar dynamo.

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.

In this letter, 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 estimate the magnitudes of the former contributions and shows that the amplitude of the produced magnetic field is about $\sim 10^{-19} {\rm G}$ at 10Mpc comoving scale at the recombination. Compared to astrophysical and inflationary mechanisms for seed-field generation, our study suffers from much less ambiguities concerning unknown physics and/or processes.

We study generation of magnetic fields by the Biermann mechanism in the pair-instability supernovae explosions of first stars. The Biermann mechanism produces magnetic fields in the shocked region between the bubble and interstellar medium (ISM), even if magnetic fields are absent initially. We perform a series of two-dimensional magnetohydrodynamic simulations with the Biermann term and estimate the amplitude and total energy of the produced magnetic fields. We find that magnetic fields with amplitude $10^{-14}-10^{-17}$ G are generated inside the bubble, though the amount of magnetic fields generated depend on specific values of initial conditions. This corresponds to magnetic fields of $10^{28}-10^{31}$ erg per each supernova remnant, which is strong enough to be the seed magnetic field for galactic and/or interstellar dynamo.

The cosmological magnetic field is one of the important physical quantities which affect strongly the cosmic microwave background (CMB) power spectrum. Recent CMB observations have been extended to higher multipoles $l\gtrsim$1000, and they resultantly exhibit an excess power than the standard model prediction in cosmological theory which best fits the Wilkinson Microwave Anisotropy Probe (WMAP) data at lower multipoles $l\lesssim$900. We calculate the CMB temperature anisotropies generated by the power-law magnetic field at the last scattering surface (LSS) in order to remove the tension between theory and observation at higher multipoles and also place an upper limit on primordial magnetic field. In our present calculation we take account of the effect of ionization ratio exactly without approximation. This effect is very crucial to precisely estimate the effect of the magnetic field on CMB power spectrum. We consider both effects of the scalar and vector modes of magnetic field on the CMB anisotropies, where current data are known to be insensitive to the tensor mode which we ignore in the present study. In order to constrain the primordial magnetic field, we evaluate likelihood function of the WMAP data in a wide range of parameters of the magnetic field strength $|\mathbf{B}|_\lambda$ and the power-law spectral index $n$, along with six cosmological parameters in flat Universe models, using the technique of the Markov Chain Monte Carlo(MCMC) method. We find that the upper limit at $2\sigma$ C.L. turns out to be $|\mathbf{B}_\lambda|\lesssim 3.9 $nG at 1 Mpc for any $n_B$ values, which is obtained by comparing the calculated result including the Sunyaev-Zeldovich(SZ) effect with recent WMAP data of the CMB anisotropies.

We present cosmological constraints on deviations of Newton's constant at large scales, analyzing latest cosmic microwave background (CMB) anisotropies and primordial abundances of light elements synthesized by big bang nucleosynthesis (BBN). BBN limits the possible deviation at typical scales of BBN epoch, say at 10^8 \sim 10^12m, to lie between -5% and +1% of the experimental value, and CMB restricts the deviation at larger scales 10^2 \sim 10^9pc to be between -26% and +66% at the 2\sigma confidence level. The cosmological constraints are compared with the astronomical one from the evolution of isochrone of globular clusters.

We construct an analytic model for the mass distributions of superclusters and clusters in each supercluster. Our model is a modification of the Press-Schechter theory, and defines superclusters as the regions that have some overdensity smaller than those for usual virialized objects. We compare the mass functions with a catalog of superclusters in the Sloan Digital Sky Survey Early Data Release and found that they are in reasonable agreement with each other.

We analyze the limits on resonant particle production during inflation based upon the power spectrum of fluctuations in matter and the cosmic microwave background. We show that such a model is consistent with features observed in the matter power spectrum deduced from galaxy surveys and damped Lyman-alpha systems at high redshift. It also provides an alternative explanation for the excess power observed in the power spectrum of the cosmic microwave background fluctuations in the range of 1000 < l < 3500. For our best-fit models, epochs of resonant particle creation reenter the horizon at wave numbers ~ 0.4 and/or 0.2 (h/Mpc). The amplitude and location of these features correspond to the creation of fermion species of mass ~ 1-2 Mpl during inflation with a coupling constant between the inflaton field and the created fermion species of near unity. Although the evidence is marginal, if this interpretation is correct, this could be one of the first observational hints of new physics at the Planck scale.

The causal structure of the flat brane universe of RSII type is re-investigated to clarify the boundary conditions for stochastic gravitational waves. In terms of the Gaussian normal coordinate of the brane, a singularity of the equation for gravitational waves appears in the bulk. We show that this singularity corresponds to the ``seam singularity'' which is a singular subspace on the brane universe. Based upon the causal structure, we discuss the boundary conditions for gravitational waves in the bulk. Introducing a null coordinate, we propose a numerical procedure to solve gravitational waves with appropriate boundary conditions and show some examples of our numerical results. This procedure could be also applied in scalar type perturbations. The problem in the choice of the initial condition for gravitational waves is briefly discussed.

We re-formulate cosmological perturbations in the decaying cold dark matter model, and calculate cosmological microwave background (CMB) anisotropies. By comparing our theoretical predictions with recent observational data from the Wilkinson Microwave Anisotropy Probe(WMAP), we derive a new bound on the abundance and lifetime of decaying dark matter particles. We show that the data of WMAP alone do not prefer the decaying cold dark matter model: the lifetime is constrained to \Gamma^{-1}> 123 Gyr at 68% C.L. (52 Gyr at 95.4% C.L.) when cold dark matter consists only of such decaying particles. We also consider a more general case in which cold dark matter consists of both stable and decaying particles, and show that the constraint generalizes down to \Omega_{DDM}h^2 < -0.5 (Gamma^{-1}/1 Gyr)^{-1}+0.12 for \Gamma^{-1} > 5 Gyr at 95.4% C.L.. These bounds are robust and widely applicable, because they are derived only from gravitational effects and therefore do not rely on the details of decay channels or decay products.

We investigate the cosmological evolution of gravitational waves in Friedman-Robertson-Walker brane world models embedded in a five dimensional anti de-Sitter spacetime. To predict the spectrum of stochastic gravitational background at present generated in the inflationary phase, we numerically calculate the evolution of gravitational waves according to the method developed in our previous paper [K. Ichiki and K. Nakamura, arXiv:hep-th/0310282]. The resulting logarithmic energy spectrum is significantly altered from that of the standard four dimensional inflationary models in the frequency range f > f_{AdS} ~ 10^{-4} Hz, which corresponds to the curvature scale we set of anti de-Sitter spacetime, $l ~ 1$mm.

Recent studies of structures of galaxies and clusters imply that dark matter might be unstable and decay with lifetime $\Gamma^{-1}$ about the age of universe. We study the effects of the decay of cold dark matter on the Sunyaev-Zel'dovich (SZ) power spectrum. We analytically calculate the SZ power spectrum taking finite lifetime of cold dark matter into account. We find the finite lifetime of dark matter decreases the power at large scale ($l < 4000$) and increases at small scale ($l > 4000$). This is in marked contrast with the dependence of other cosmological parameters such as the amplitude of mass fluctuations $\sigma_{8}$ and the cosmological constant $\Omega_{\lambda}$ (under the assumption of a flat universe) which mainly change the normalization of the angular power spectrum. This difference allows one to determine the lifetime and other cosmological parameters rather separately. We also investigate sensitivity of a future SZ survey to the cosmological parameters including the life time, assuming a fiducial model $\Gamma^{-1} = 10 h^{-1} {\rm Gyr}$, $\sigma_{8} = 1.0$, and $\Omega_{\lambda} = 0.7$. We show that future SZ surveys such as ACT, AMIBA, and BOLOCAM can determine the lifetime within factor of two even if $\sigma_{8}$ and $\Omega_{\lambda}$ are marginalized.

We explore cosmological implications of dark matter as massive particles trapped on a brane embedded in a Randall-Sundrum noncompact higher dimension $AdS_5$ space. It is an unavoidable consequence of this cosmology that massive particles are metastable and can disappear into the bulk dimension. Here, we show that a massive dark matter particle (e.g. the lightest supersymmetric particle) is likely to have the shortest lifetime for disappearing into the bulk. We examine cosmological constraints on this new paradigm and show that disappearing dark matter is consistent (at the 95% confidence level) with all cosmological constraints, i.e. present observations of Type Ia supernovae at the highest redshift, trends in the mass-to-light ratios of galaxy clusters with redshift, the fraction of X-ray emitting gas in rich clusters, and the spectrum of power fluctuations in the cosmic microwave background. A best $2 \sigma$ concordance region is identified corresponding to a mean lifetime for dark matter disappearance of $15 \le \Gamma^{-1} \le 80$ Gyr. The implication of these results for brane-world physics is discussed.

It has recently been suggested that the formation of horizon-size primordial black hole (PBH) from pre-existing density fluctuations is effective during the cosmic QCD phase transition. In this Letter we discuss the dependence of PBH formation on effective relativistic degrees of freedom, $g_{\rm eff}$ during the cosmic QCD phase transition. Our finding is important in the light of recent cosmological arguments of several new classes of cosmic energy that appear from universal neutrino degeneracy, quintessential inflation, and dark radiation in brane world cosmology. Extra-energy component from the standard value in these new cosmological theories is represented as an effective radiation in terms of $g_{\rm eff}$. We conclude that the PBH formation during QCD phase transition becomes more efficient if negative extra-component of the cosmic energy is allowed because of the increase of the duration of the QCD phase transition, which leads to smaller mass scale of PBHs. This suggests larger probability of finding more PBHs if the dark radiation exists as allowed in the brane world cosmology.

We analyze the observational constraints on brane-world cosmology whereby the universe is described as a three-brane embedded in a five-dimensional anti-de Sitter space. In this brane-universe cosmology, the Friedmann equation is modified by the appearance of extra terms which derive from existence of the extra dimensions. In the present work we concentrate on the ``dark radiation'' term which diminishes with cosmic scale factor as $a^{-4}$. We show that, although the observational constraints from primordial abundances allow only a small contribution when this term is positive, a much wider range of negative values is allowed. Furthermore, such a negative contribution can reconcile the tension between the observed primordial $\he4$ and D abundances. We also discuss the possible constraints on this term from the power spectrum of CMB anisotropies in the limit of negligible cosmological perturbation on the brane world. We show that BBN limits the possible contribution from dark radiation just before the nucleosynthesis epoch to lie between -65% and $+5%$ of the background photon energy density. Combining this with the CMB constraint reduces this range to between -24% and $+3.5%$ at the $2\sigma$ confidence level.

A generic feature of cosmic quintessence models is the possibility of significant energy density in the scalar field during the radiation dominated epoch. This possibility is even greater if the quintessence field begins in a kinetic-dominated regime, for example as might be generated at the end of "quintessential inflation." As such, these models can be constrained by primordial nucleosynthesis and the epoch of photon decoupling. Here, we analyze both kinetic dominated and power-law quintessence fields (with and without a supergravity correction). We quantify the allowed initial conditions and effective-potential parameters. We also deduce constraints on the epoch of matter creation at the end of quintessential inflation.

Recently, attempts have been made to understand the apparent near coincidence of the present dark energy and matter energy in terms of a dynamical attractor-like solution for the evolution of a scalar field. In these models the field couples with the dominant constituent and only acts like a cosmological constant after the onset of the matter-diminated epoch. A generic feature of such solutions, however, is the possibility of significant energy density in the scalar field during the radiation-dominated epoch. This possibility is even greater if the quintessence field begins in a kinetic-dominated regime generated at the end of quintessential inflation. As such, these models can affect and be constrained by primordial nucleosynthesis and the epoch of of photon decoupling. Here, we analyze one popular form of the quintessence field (with and without a supergravity correction) and quantify constraints on the allowed initial conditions and parameters for the effective potential. We also deduce constraints on the epoch of matter creation at the end of quintessential inflation.