We investigate the possibility of measuring the primordial gravitational wave (GW) signal across 23 decades in frequencies, using the cosmic microwave background (CMB), pulsar timing arrays (PTA), and direct detection with laser and atomic interferometers. For the CMB and PTA experiments we consider the LiteBIRD mission and the Square Kilometer Array (SKA), respectively. For the interferometers we consider space mission proposals including the Laser Interferometer Space Antenna (LISA), the Big Bang Observer (BBO), the Deci-hertz Interferometer Gravitational wave Observatory (DECIGO), the $\mu$Ares experiment, the Decihertz Observatory (DO), and the Atomic Experiment for Dark Matter and Gravity Exploration in Space (AEDGE), as well as the ground-based Einstein Telescope (ET) and Cosmic Explorer (CE) proposals. We implement the mathematics needed to compute sensitivities for both CMB and interferometers, and derive the response functions for the latter from the first principles. We also evaluate the effect of the astrophysical foreground contamination in each experiment. We present binned sensitivity curves and error bars on the energy density parameter, $\Omega_{GW}h^2$, as a function of frequency for two representative classes of models for the stochastic background of primordial GW: the quantum vacuum fluctuation in the metric from single-field slow-roll inflation, and the source-induced tensor perturbation from the spectator axion-SU(2) inflation models. We find excellent prospects for joint measurements of the GW spectrum by CMB and space-borne direct detection mission proposals.

As cosmic structures form, matter density fluctuations collapse gravitationally and baryonic matter is shock-heated and thermalized. We therefore expect a connection between the mean gravitational potential energy density of collapsed halos, $\Omega_{W}^{\rm halo}$, and the mean thermal energy density of baryons, $\Omega_{\rm th}$. These quantities can be obtained using two fundamentally different estimates: we compute $\Omega_{W}^{\rm halo}$ using the theoretical framework of the halo model which is driven by dark matter statistics, and measure $\Omega_{\rm th}$ using the Sunyaev-Zeldovich (SZ) effect which probes the mean thermal pressure of baryons. First, we derive that, at the present time, about 90% of $\Omega_{W}^{\rm halo}$ originates from massive halos with $M>10^{13}\,M_\odot$. Then, using our measurements of the SZ background, we find that $\Omega_{\rm th}$ accounts for about 80% of the kinetic energy of the baryons available for pressure in halos at $z\lesssim 0.5$. This constrains the amount of non-thermal pressure, e.g., due to bulk and turbulent gas motion sourced by mass accretion, to be about $\Omega_{\rm non-th}\simeq 0.4\times 10^{-8}$ at $z=0$.

We develop a strategy to determine the cosmic birefringence and miscalibrated polarization angles simultaneously using the observed $EB$ polarization power spectra of the cosmic microwave background and the Galactic foreground emission. We extend the methodology of Y. Minami et al. (Prog. Theor. Exp. Phys. $\mathbf{2019}$, 083E02, 2019), which was developed for auto frequency power spectra, by including cross frequency spectra. By fitting one global birefringence angle and independent miscalibration angles at different frequency bands, we determine both angles with significantly smaller uncertainties (by more than a factor of two) compared to the auto spectra.

The cosmic thermal history, quantified by the evolution of the mean thermal energy density in the universe, is driven by the growth of structures as baryons get shock-heated in collapsing dark matter halos. This process can be probed by redshift-dependent amplitudes of the thermal Sunyaev-Zeldovich (SZ) effect background. To do so, we cross-correlate eight sky intensity maps in the $\it{Planck}$ and Infrared Astronomical Satellite missions with two million spectroscopic-redshift references in the Sloan Digital Sky Surveys. This delivers snapshot spectra for the far-infrared to microwave background light as a function of redshift up to $z\sim3$. We decompose them into the SZ and thermal dust components. Our SZ measurements directly constrain $\langle bP_{\rm e} \rangle$, the halo bias-weighted mean electron pressure, up to $z\sim 1$. This is the highest redshift achieved to date, with uncorrelated redshift bins thanks to the spectroscopic references. We detect a threefold increase in the density-weighted mean electron temperature $\bar{T}_{\rm{e}}$ from $7\times 10^5~{\rm K}$ at $z=1$ to $2\times 10^6~{\rm K}$ today. Over $z=1$-$0$, we witness the build up of nearly $70\%$ of the present-day mean thermal energy density $\rho_{\rm{th}}$, with the corresponding density parameter $\Omega_{\rm th}$ reaching $1.5 \times10^{-8}$. We find the mass bias parameter of $\it{Planck}$'s universal pressure profile of $B=1.27$ (or $1-b=1/B=0.79$), consistent with the magnitude of non-thermal pressure in gas motion and turbulence from mass assembly. We estimate the redshift-integrated mean Compton parameter $y\sim1.2\times10^{-6}$, which will be tested by future spectral distortion experiments. More than half of which originates from large-scale structure at $z<1$, which we detect directly.

Building on Weinberg's approach to effective field theory for inflation, we construct an effective Lagrangian for a pseudo scalar (axion) inflaton field with shift symmetry. In this Lagrangian we allow the axion field to couple to non-Abelian gauge fields via a Chern-Simons term. We then analyze a class of inflation models driven by kinetic terms. We find that the observational constraints on the amplitudes of curvature perturbations and non-Gaussianity yield a lower bound for the tensor-to-scalar ratio of $r\gtrsim 5\times 10^{-3}$ from the vacuum fluctuation. The sourced gravitational wave from SU(2) gauge fields further increases the tensor-to-scalar ratio and makes the total gravitational wave partially chiral and non-Gaussian, which can be probed by polarization of the cosmic microwave background and direct detection experiments. We discuss constraints on parameter space due to backreaction of spin-2 particles produced by the gauge field.

In this paper, we examine the viability of inflation models with a spectator axion field coupled to both gravitational and SU(2) gauge fields via Chern-Simons couplings. Requiring phenomenological success of the axion-SU(2) sector constrains the coupling strength of the gravitational Chern-Simons term. We find that the impact of this term on the production and propagation of gravitational waves can be as large as fifty percent enhancement for the helicity that is not sourced by the gauge field, if the cut-off scale is as low as $\Lambda$ = 20H. The effect becomes smaller for a larger value of $\Lambda$, while the impact on the helicity sourced by the gauge field is negligible regardless of $\Lambda$.

The key to the phenomenological success of inflation models with axion and SU(2) gauge fields is the isotropic background of the SU(2) field. Previous studies showed that this isotropic background is an attractor solution during inflation starting from anisotropic (Bianchi Type I) spacetime; however, not all possible initial anisotropic parameter space was explored. In this paper, we explore more generic initial conditions without assuming the initial slow-roll dynamics. We find some initial anisotropic parameter space which does not lead to the isotropic background, but to violation of slow-roll conditions, terminating inflation prematurely. The basin of attraction increases when we introduce another scalar field acting as inflaton and make the axion-SU(2) system a spectator sector. Therefore, the spectator axion-SU(2) model is phenomenologically more attractive.

Thermal Sunyaev-Zeldovich (tSZ) power spectrum is a powerful probe of the present-day amplitude of matter density fluctuations, and has been measured up to $\ell\approx 10^3$ from the {\it Planck} data. The largest systematic uncertainty in the interpretation of this data is the so-called "mass bias" parameter $B$, which relates the true halo mass to the mass proxy used by the {\it Planck} team as $M_{\rm 500c}^{Planck}=M_{\rm 500c}^{\rm true}/B$. Since the power spectrum of the cosmic weak lensing shear is also sensitive to the amplitude of matter density fluctuations via $S_8\equiv \sigma_8 \Omega_m^{\alpha}$ with $\alpha\sim 0.5$, we can break the degeneracy between the mass bias and the cosmological parameters by combining the tSZ and cosmic shear power spectra. In this paper, we perform a joint likelihood analysis of the tSZ power spectrum from {\it Planck} and the cosmic shear power spectrum from Subaru Hyper Suprime-Cam. Our analysis does not use the primordial cosmic microwave background (CMB) information. We obtain a new constraint on the mass bias as $B = 1.37 ^{+0.15}_{-0.23}$ or $(1-b) = B^{-1}=0.73^{+0.08}_{-0.13}$ (68\%~C.L.), for $\sigma_8 < 0.9$. This value of $B$ is lower than that needed to reconcile the tSZ data with the primordial CMB and CMB lensing data, i.e., $B = 1.64 \pm 0.19$, but is consistent with the mass bias expected from hydrodynamical simulations, $B = 1.28 \pm 0.20$. Our results thus indicate that the mass bias is consistent with the non-thermal pressure support from mass accretion of galaxy clusters via the cosmic structure formation, and that the cosmologies inferred from the tSZ and the cosmic shear are consistent with each other.

We present measurements of the thermal Sunyaev-Zel'dovich effect (SZE) toward SPT-CL J2334-4243 (the Phoenix galaxy cluster) at z=0.597 by the Atacama Large Millimeter/submillimeter Array (ALMA) in Band 3. The SZE is imaged at 5" resolution (corresponding to the physical scale of 23kpc/h) within 200kpc/h from the central galaxy with the peak signal-to-noise ratio exceeding 11. Combined with the Chandra X-ray image, the ALMA SZE data further allow for non-parametric deprojection of electron temperature, density, and entropy. Our method can minimize contamination by the central AGN and the X-ray absorbing gas within the cluster, both of which largely affect the X-ray spectrum. We find no significant asymmetry or disturbance in the SZE image within the current measurement errors. The detected SZE signal agrees well with the average pressure profile of local cool-core clusters. Unlike any other known clusters, however, gas temperature drops by at least a factor of 5 toward the center. We identify ~6x10^{11} M_sun cool gas with temperature ~3keV in the inner 20kpc/h. Taken together, our results imply that the gas is indeed cooling efficiently and nearly isobarically down to this radius in the Phoenix cluster.

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 examine the impact of fiber assignment on clustering measurements from fiber-fed spectroscopic galaxy surveys. We identify new effects which were absent in previous, relatively shallow galaxy surveys such as Baryon Oscillation Spectroscopic Survey . Specifically, we consider deep surveys covering a wide redshift range from z=0.6 to z=2.4, as in the Subaru Prime Focus Spectrograph survey. Such surveys will have more target galaxies than we can place fibers on. This leads to two effects. First, it eliminates fluctuations with wavelengths longer than the size of the field of view, as the number of observed galaxies per field is nearly fixed to the number of available fibers. We find that we can recover the long-wavelength fluctuation by weighting galaxies in each field by the number of target galaxies. Second, it makes the preferential selection of galaxies in under-dense regions. We mitigate this effect by weighting galaxies using the so-called individual inverse probability. Correcting these two effects, we recover the underlying correlation function at better than 1 percent accuracy on scales greater than 10 Mpc/h.

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints.

Redshifts of an astronomical body measured at multiple epochs (e.g., separated by 10 years) are different due to the cosmic expansion. This so-called Sandage-Loeb test offers a direct measurement of the expansion rate of the Universe. However, acceleration in the motion of Solar System with respect to the cosmic microwave background also changes redshifts measured at multiple epochs. If not accounted for, it yields a biased cosmological inference. To address this, we calculate the acceleration of Solar System with respect to the Local Group of galaxies to quantify the change in the measured redshift due to local motion. Our study is motivated by the recent determination of the mass of Large Magellanic Cloud (LMC), which indicates a significant fraction of the Milky Way mass. We find that the acceleration towards the Galactic Center dominates, which gives a redshift change of 7 cm/s in 10 years, while the accelerations due to LMC and M31 cannot be ignored depending on lines of sight. We create all-sky maps of the expected change in redshift and the corresponding uncertainty, which can be used to correct for this effect.

The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$.

The Cerro Chajnantor Atacama Telescope-prime (CCAT-prime) is a new 6-m, off-axis, low-emissivity, large field-of-view submillimeter telescope scheduled for first light in the last quarter of 2021. In summary, (a) CCAT-prime uniquely combines a large field-of-view (up to 8-deg), low emissivity telescope (< 2%) and excellent atmospheric transmission (5600-m site) to achieve unprecedented survey capability in the submillimeter. (b) Over five years, CCAT-prime first generation science will address the physics of star formation, galaxy evolution, and galaxy cluster formation; probe the re-ionization of the Universe; improve constraints on new particle species; and provide for improved removal of dust foregrounds to aid the search for primordial gravitational waves. (c) The Observatory is being built with non-federal funds (~ \$40M in private and international investments). Public funding is needed for instrumentation (~ \$8M) and operations (\$1-2M/yr). In return, the community will be able to participate in survey planning and gain access to curated data sets. (d) For second generation science, CCAT-prime will be uniquely positioned to contribute high-frequency capabilities to the next generation of CMB surveys in partnership with the CMB-S4 and/or the Simons Observatory projects or revolutionize wide-field, sub-millimetter line intensity mapping surveys.

We derive a `Kompaneets equation' for neutrinos, which describes how the distribution function of neutrinos interacting with matter deviates from a Fermi-Dirac distribution with zero chemical potential. To this end, we expand the collision integral in the Boltzmann equation of neutrinos up to the second order in energy transfer between matter and neutrinos. The distortion of the neutrino distribution function changes the rate at which neutrinos heat matter, as the rate is proportional to the mean square energy of neutrinos, $E_\nu^2$. For electron-type neutrinos the enhancement in $E_\nu^2$ over its thermal value is given approximately by $E_\nu^2/E_{\nu,\rm thermal}^2=1+0.086(V/0.1)^2$ where $V$ is the bulk velocity of nucleons, while for the other neutrino species the enhancement is $(1+\delta_v)^3$, where $\delta_v=mV^2/3k_BT$ is the kinetic energy of nucleons divided by the thermal energy. This enhancement has a significant implication for supernova explosions, as it would aid neutrino-driven explosions.

Context. The precise determination of the present-day expansion rate of the Universe, expressed through the Hubble constant $H_0$, is one of the most pressing challenges in modern cosmology. Assuming flat $\Lambda$CDM, $H_0$ inference at high redshift using cosmic-microwave-background data from Planck disagrees at the 4.4$\sigma$ level with measurements based on the local distance ladder made up of parallaxes, Cepheids and Type Ia supernovae (SNe Ia), often referred to as "Hubble tension". Independent, cosmological-model-insensitive ways to infer $H_0$ are of critical importance. Aims. We apply an inverse-distance-ladder approach, combining strong-lensing time-delay-distance measurements with SN Ia data. By themselves, SNe Ia are merely good relative distance indicators, but by anchoring them to strong gravitational lenses one can obtain an $H_0$ measurement that is relatively insensitive to other cosmological parameters. Methods. A cosmological parameter estimate is performed for different cosmological background models, both for strong-lensing data alone and for the combined lensing + SNe Ia data sets. Results. The cosmological-model dependence of strong-lensing $H_0$ measurements is significantly mitigated through the inverse distance ladder. In combination with SN Ia data, the inferred $H_0$ consistently lies around 73-74 km s$^{-1}$ Mpc$^{-1}$, regardless of the assumed cosmological background model. Our results agree nicely with those from the local distance ladder, but there is a >2$\sigma$ tension with Planck results, and a ~1.5$\sigma$ discrepancy with results from an inverse distance ladder including Planck, Baryon Acoustic Oscillations and SNe Ia. Future strong-lensing distance measurements will reduce the uncertainties in $H_0$ from our inverse distance ladder.

The galaxy catalogs generated from low-resolution emission line surveys often contain both foreground and background interlopers due to line misidentification, which can bias the cosmological parameter estimation. In this paper, we present a method for correcting the interloper bias by using the joint-analysis of auto- and cross-power spectra of the main and the interloper samples. In particular, we can measure the interloper fractions from the cross-correlation between the interlopers and survey galaxies, because the true cross-correlation must be negligibly small. The estimated interloper fractions, in turn, remove the interloper bias in the cosmological parameter estimation. For example, in the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) low-redshift ($z<0.5$) [O II] $\lambda3727${\AA} emitters contaminate high-redshift ($1.9<z<3.5$) Lyman-$\alpha$ line emitters. We demonstrate that the joint-analysis method yields a high signal-to-noise ratio measurement of the interloper fractions while only marginally increasing the uncertainties in the cosmological parameters relative to the case without interlopers. We also show the same is true for the high-latitude spectroscopic survey of Wide-Field Infrared Survey Telescope (WFIRST) mission where contamination occurs between the Balmer-$\alpha$ line emitters at lower redshifts ($1.1<z<1.9$) and Oxygen ([O III] $\lambda5007${\AA}) line emitters at higher redshifts ($1.7<z<2.8$).

We perform forecasts for how baryon acoustic oscillation (BAO) scale and redshift-space distortion (RSD) measurements from future spectroscopic emission line galaxy (ELG) surveys such as Euclid are degraded in the presence of spectral line misidentification. Using analytic calculations verified with mock galaxy catalogs from log-normal simulations we find that constraints are degraded in two ways, even when the interloper power spectrum is modeled correctly in the likelihood. Firstly, there is a loss of signal-to-noise ratio for the power spectrum of the target galaxies, which propagates to all cosmological constraints and increases with contamination fraction, $f_c$. Secondly, degeneracies can open up between $f_c$ and cosmological parameters. In our calculations this typically increases BAO scale uncertainties at the 10-20% level when marginalizing over parameters determining the broadband power spectrum shape. External constraints on $f_c$, or parameters determining the shape of the power spectrum, for example from cosmic microwave background (CMB) measurements, can remove this effect. There is a near-perfect degeneracy between $f_c$ and the power spectrum amplitude for low $f_c$ values, where $f_c$ is not well determined from the contaminated sample alone. This has the potential to strongly degrade RSD constraints. The degeneracy can be broken with an external constraint on $f_c$, for example from cross-correlation with a separate galaxy sample containing the misidentified line, or deeper sub-surveys.

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.

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

The presence of a light vector field coupled to a scalar field during inflation makes a distinct prediction: the observed correlation functions of the cosmic microwave background (CMB) become statistically anisotropic. We study the implications of the current bound on statistical anisotropy derived from the Planck 2013 CMB temperature data for such a model. The previous calculations based on the attractor solution indicate that the magnitude of anisotropy in the power spectrum is proportional to $N^2$, where $N$ is the number of $e$-folds of inflation counted from the end of inflation. In this paper, we show that the attractor solution is not compatible with the current bound, and derive new predictions using another branch of anisotropic inflation. In addition, we improve upon the calculation of the mode function of perturbations by including the leading-order slow-roll corrections. We find that the anisotropy is roughly proportional to $[2(\varepsilon_H+4\eta_H)/3-4(c-1)]^{-2}$, where $\varepsilon_H$ and $\eta_H$ are the usual slow-roll parameters and $c$ is the parameter in the model, regardless of the form of potential of an inflaton field. The bound from Planck implies that breaking of rotational invariance during inflation (characterized by the background homogeneous shear divided by the Hubble rate) is limited to be less than ${\cal O}(10^{-9})$. This bound is many orders of magnitude smaller than the amplitude of breaking of time translation invariance, which is observed to be ${\cal O}(10^{-2})$.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

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

The 21-cm signal from the epoch of reionization is non-Gaussian. Current radio telescopes are focused on detecting the 21-cm power spectrum, but in the future the Square Kilometre Array is anticipated to provide a first measurement of the bispectrum. Previous studies have shown that the position-dependent power spectrum is a simple and efficient way to probe the squeezed-limit bispectrum. In this approach, the survey is divided into subvolumes and the correlation between the local power spectrum and the corresponding mean density of the subvolume is computed. This correlation is equivalent to an integral of the bispectrum in the squeezed limit, but is much simpler to implement than the usual bispectrum estimators. It also has a clear physical interpretation: it describes how the small-scale power spectrum of tracers such as galaxies and the 21-cm signal respond to a large-scale environment. Reionization naturally couples large and small scales as ionizing radiation produced by galactic sources can travel up to tens of Megaparsecs through the intergalactic medium during this process. Here we apply the position-dependent power spectrum approach to fluctuations in the 21-cm background from reionization. We show that this statistic has a distinctive evolution in time that can be understood with a simple analytic model. We also show that the statistic can easily distinguish between simple "inside-out" and "outside-in" models of reionization. The position-dependent power spectrum is thus a promising method to validate the reionization signal and to extract higher-order information on this process.

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 discovery of the electromagnetic counterpart to GW170817 severely constrains the tensor mode propagation speed, eliminating a large model space of Horndeski theory. We use the cosmic microwave background data from Planck and the joint analysis of the BICEP2/Keck Array and Planck, galaxy clustering data from the SDSS LRG survey, BOSS baryon acoustic oscillation data, and redshift space distortion measurements to place constraints on the remaining Horndeski parameters. We evolve the Horndeski parameters as power laws with both the amplitude and power law index free. We find a 95% CL upper bound on the present-day coefficient of the Hubble friction term in the cosmological propagation of gravitational waves is 2.38, whereas General Relativity gives 2 at all times. While an enhanced friction suppresses the amplitude of the reionization bump of the primordial B-mode power spectrum at $\ell < 10$, our result limits the suppression to be less than 0.8%. This constraint is primarily due to the scalar integrated Sachs-Wolfe effect in temperature fluctuations at low multipoles.

Non-Abelian gauge fields may exist during inflation. We study the Schwinger effect by an $SU(2)$ gauge field coupled to a charged scalar doublet in a (quasi) de Sitter background and the possible backreaction of the generated charged particles on the homogeneous dynamics. Contrary to the Abelian $U(1)$ case, we find that both the Schwinger pair production and the induced current decrease as the interaction strength increases. The reason for this suppression is the isotropic vacuum expectation value of the $SU(2)$ field which generates a (three times) greater effective mass for the scalar field than the $U(1)$. In the weak interaction limit, the above effect is negligible and both the $SU(2)$ and $U(1)$ cases exhibit a linear increase of the current and a constant conductivity with the interaction strength. We conclude that the Schwinger effect does not pose a threat to the dynamics of inflationary models involving an $SU(2)$ gauge field.

We present a first joint analysis of the power spectra of the thermal Sunyaev-Zeldovich (tSZ) effect measured by the Planck and the number density fluctuations of galaxies in the 2MASS redshift survey (2MRS) catalog, including their cross-correlation. Combining these measurements with the cosmic microwave background (CMB) data and CMB lensing of Planck assuming a flat $\Lambda$CDM model, we constrain the mass bias parameter as $B = 1.54 \pm 0.098\;(1\sigma)$ [$(1-b) = 0.649 \pm 0.041$, where $(1-b) \equiv B^{-1}$], i.e., the Planck cluster mass should be $35\%$ lower than the true mass. The mass bias determined by the 2MRS-tSZ cross-power spectrum alone is consistent with that determined by the tSZ auto-power spectrum alone, suggesting that this large mass bias is not due to obvious systematics in the tSZ data. We find that the 2MRS-tSZ cross-power spectrum is more sensitive to less massive halos than the tSZ auto-power spectrum and it significantly improves a constraint on the mass dependence of the mass bias. The redshift dependence is not strongly constrained since the multipole range in which high redshift clusters mainly contribute to the tSZ auto is dominated by the contaminating sources. We conclude that no strong mass or redshift evolution of the mass bias is needed to explain the data.

RXJ1347.5-1145 (z = 0.451) is one of the most luminous X-ray galaxy clusters, which hosts a prominent cool core and exhibits a signature of a major merger. We present the first direct observational evidence for sub-sonic nature of sloshing motion of the cool core. We find that a residual X-ray image from the Chandra X-ray Observatory after removing the global emission shows a clear dipolar pattern characteristic of gas sloshing, whereas we find no significant residual in the Sunyaev-Zel'dovich effect (SZE) image from the Atacama Large Millimeter/submillimeter Array (ALMA). We estimate the equation of state of perturbations in the gas from the X-ray and SZE residual images. The inferred velocity is 420 +310 -420 km s-1, which is much lower than the adiabatic sound speed of the intracluster medium in the core. We thus conclude that the perturbation is nearly isobaric, and gas sloshing motion is consistent with being in pressure equilibrium. Next, we report evidence for gas stripping of an infalling subcluster, which likely shock-heats gas to high temperature well in excess of 20 keV. Using mass distribution inferred from strong lensing images of the Hubble Space Telescope (HST), we find that the mass peak is located away from the peak position of stripped gas with statistical significance of > 5{\sigma}. Unlike for the gas sloshing, the velocity inferred from the equation of state of the excess hot gas is comparable to the adiabatic sound speed expected for the 20 keV intracluster medium. All of the results support that the southeast substructure is created by a merger. On the other hand, the positional offset between the mass and the gas limits the self-interaction cross section of dark matter to be less than 3.7 h-1 cm2 g-1 (95% CL).

We calculate the bispectrum of scale-invariant tensor modes sourced by spectator SU(2) gauge fields during inflation in a model containing a scalar inflaton, a pseudoscalar axion and SU(2) gauge fields. A large bispectrum is generated in this model at tree-level as the gauge fields contain a tensor degree of freedom, and its production is dominated by self-coupling of the gauge fields. This is a unique feature of non-Abelian gauge theory. The shape of the tensor bispectrum is approximately an equilateral shape for $3\lesssim m_Q\lesssim 4$, where $m_Q$ is an effective dimensionless mass of the SU(2) field normalised by the Hubble expansion rate during inflation. The amplitude of non-Gaussianity of the tensor modes, characterised by the ratio $B_h/P^2_h$, is inversely proportional to the energy density fraction of the gauge field. This ratio can be much greater than unity, whereas the ratio from the vacuum fluctuation of the metric is of order unity. The bispectrum is effective at constraining large $m_Q$ regions of the parameter space, whereas the power spectrum constrains small $m_Q$ regions.

Given observations of B-mode polarization power spectrum of the cosmic microwave background (CMB), we can reconstruct power spectra of primordial tensor modes from the early Universe without assuming their functional form such as a power-law spectrum. Shape of the reconstructed spectra can then be used to probe the origin of tensor modes in a model-independent manner. We use the Fisher matrix to calculate the covariance matrix of tensor power spectra reconstructed in bins. We find that the power spectra are best reconstructed at wavenumbers in the vicinity of $k\approx 6\times 10^{-4}$ and $5\times 10^{-3}~{\rm Mpc}^{-1}$, which correspond to the "reionization bump" at $\ell\lesssim 6$ and "recombination bump" at $\ell\approx 80$ of the CMB B-mode power spectrum, respectively. The error bar between these two wavenumbers is larger because of lack of the signal between the reionization and recombination bumps. The error bars increase sharply towards smaller (larger) wavenumbers because of the cosmic variance (CMB lensing and instrumental noise). To demonstrate utility of the reconstructed power spectra we investigate whether we can distinguish between various sources of tensor modes including those from the vacuum metric fluctuation and SU(2) gauge fields during single-field slow-roll inflation, open inflation and massive gravity inflation. The results depend on the model parameters, but we find that future CMB experiments are sensitive to differences in these models. We make our calculation tool available on-line.

We show that an inflation model in which a spectator axion field is coupled to an SU(2) gauge field produces a large three-point function (bispectrum) of primordial gravitational waves, $B_{h}$, on the scales relevant to the cosmic microwave background experiments. The amplitude of the bispectrum at the equilateral configuration is characterized by $B_{h}/P_h^2=\mathcal{O}(10)\times \Omega_A^{-1}$, where $\Omega_A$ is a fraction of the energy density in the gauge field and $P_h$ is the power spectrum of gravitational waves produced by the gauge field.

We constrain the dark energy equation of state parameter, $w$, using the power spectrum of the thermal Sunyaev-Zeldovich (tSZ) effect. We improve upon previous analyses by taking into account the trispectrum in the covariance matrix and marginalising over the foreground parameters, the correlated noise, the mass bias $B$ in the Planck universal pressure profile, and all the relevant cosmological parameters (i.e., not just $\Omega_{\mathrm{m}}$ and $\sigma_8$). We find that the amplitude of the tSZ power spectrum at $\ell\lesssim 10^3$ depends primarily on $F\equiv \sigma_{8}(\Omega_{{\mathrm{m}}}/B)^{0.40}h^{-0.21}$, where $B$ is related to more commonly used variable $b$ by $B=(1-b)^{-1}$. We measure this parameter with 2.6\% precision, $F=0.460\pm 0.012$ (68% CL). By fixing the bias to $B=1.25$ and adding the local determination of the Hubble constant $H_0$ and the amplitude of the primordial power spectrum constrained by the Planck Cosmic Microwave Background (CMB) data, we find $w=-1.10\pm0.12$, $\sigma_{\mathrm{8}}=0.802\pm0.037$, and $\Omega_{{\mathrm{m}}}=0.265\pm0.022$ (68% CL). Our limit on $w$ is consistent with and is as tight as that from the distance-alone constraint from the CMB and $H_0$. Finally, by combining the tSZ power spectrum and the CMB data we find, in the $\Lambda$ Cold Dark Matter (CDM) model, the mass bias of $B=1.71\pm 0.17$, i.e., $1-b=0.58\pm 0.06$ (68% CL).

The total mass of neutrinos can be constrained in a number of ways using galaxy redshift surveys. Massive neutrinos modify the expansion rate of the Universe, which can be measured using baryon acoustic oscillations (BAOs) or the Alcock-Paczynski (AP) test. Massive neutrinos also change the structure growth rate and the amplitude of the matter power spectrum, which can be measured using redshift-space distortions (RSD). We use the Fisher matrix formalism to disentangle these information sources, to provide projected neutrino mass constraints from each of these probes alone and to determine how sensitive each is to the assumed cosmological model. We isolate the distinctive effect of neutrino free-streaming on the matter power spectrum and structure growth rate as a signal unique to massive neutrinos that can provide the most robust constraints, which are relatively insensitive to extensions to the cosmological model beyond $\Lambda$CDM. We also provide forecasted constraints using all of the information contained in the observed galaxy power spectrum combined, and show that these maximally optimistic constraints are primarily limited by the accuracy to which the optical depth of the cosmic microwave background, $\tau$, is known.

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.

A detection of B-mode polarization of the Cosmic Microwave Background (CMB) anisotropies would confirm the presence of a primordial gravitational wave background (GWB). In the inflation paradigm this would be an unprecedented probe of the energy scale of inflation as it is directly proportional to the power spectrum of the GWB. However, similar tensor perturbations can be produced by the matter fields present during inflation, breaking this simple relationship. It is therefore important to be able to distinguish between different generation mechanisms of the GWB. In this paper, we analyse the detectability of a new axion-SU(2) gauge field model using its chiral, scale-dependent tensor spectrum. We forecast the detectability of the resulting CMB TB and EB cross-correlations by the LiteBIRD satellite, considering the effects of residual foregrounds, gravitational lensing, and for the first time assess the ability of such an experiment to jointly detect primordial TB and EB spectra and self-calibrate its polarimeter. We find that LiteBIRD will be able to detect the chiral signal for $r_*>0.03$ with $r_*$ denoting the tensor-to-scalar ratio at the peak scale, and that the maximum signal-to-noise for $r_*<0.07$ is $\sim 2$. We go on to consider an advanced stage of a LISA-like mission, and find that such experiments would complement CMB observations by providing sensitivity to GWB chirality on scales inaccessible to the CMB. We conclude that in order to use the CMB to distinguish this model from a conventional vacuum fluctuation model two-point statistics provide some power, but to achieve high statistical significance we would require higher order statistics which take advantage of the model's non-Gaussianity. On the other hand, in the case of a spectrum peaked at very small scales, inaccessible to the CMB, a highly significant detection could be made using space-based laser interferometers.

We study cross-correlations of the kinetic Sunyaev-Zel'dovich effect (kSZ) and 21 cm signals during the epoch of reionisation (EoR) to measure the effects of patchy reionisation. Since the kSZ effect is proportional to the line-of-sight velocity, the kSZ-21 cm cross correlation suffers from cancellation at small angular scales. We thus focus on the correlation between the kSZ-squared field (kSZ$^2$) and 21 cm signals. When the global ionisation fraction is low ($x_e\lesssim 0.7$), the kSZ$^2$ fluctuation is dominated by rare ionised bubbles which leads to an anti-correlation with the 21 cm signal. When $0.8\lesssim x_e<1$, the correlation is dominated by small pockets of neutral regions, leading to a positive correlation. However, at very high redshifts when $x_e<0.15$, the spin temperature fluctuations change the sign of the correlation from negative to positive, as weakly ionised regions can have strong 21 cm signals in this case. To extract this correlation, we find that Wiener filtering is effective in removing large signals from the primary CMB anisotropy. The expected signal-to-noise ratios for a $\sim$10-hour integration of upcoming Square Kilometer Array data cross-correlated with maps from the current generation of CMB observatories with 3.4~$\mu$K arcmin noise and 1.7~arcmin beam over 100~deg$^2$ are 51, 60, and 37 for $x_e=0.2$, 0.5, and 0.9, respectively.

We present the measurement and interpretation of the angular power spectrum of nearby galaxies in the 2MASS Redshift Survey catalog with spectroscopic redshifts up to $z\approx 0.1$. We detect the angular power spectrum up to a multipole of $\ell\approx 1000$. We find that the measured power spectrum is dominated by galaxies living inside nearby galaxy clusters and groups. We use the halo occupation distribution (HOD) formalism to model the power spectrum, obtaining a fit with reasonable parameters. These HOD parameters are in agreement with the 2MASS galaxy distribution we measure toward the known nearby galaxy clusters, confirming validity of our analysis.

We present a public code to generate a mock galaxy catalog in redshift space assuming a log-normal probability density function (PDF) of galaxy and matter density fields. We draw galaxies by Poisson-sampling the log-normal field, and calculate the velocity field from the linearised continuity equation of matter fields, assuming zero vorticity. This procedure yields a PDF of the pairwise velocity fields that is qualitatively similar to that of N-body simulations. We check fidelity of the catalog, showing that the measured two-point correlation function and power spectrum in real space agree with the input precisely. We find that a linear bias relation in the power spectrum does not guarantee a linear bias relation in the density contrasts, leading to a cross-correlation coefficient of matter and galaxies deviating from unity on small scales. We also find that linearising the Jacobian of the real-to-redshift space mapping provides a poor model for the two-point statistics in redshift space. That is, non-linear redshift-space distortion is dominated by non-linearity in the Jacobian. The power spectrum in redshift space shows a damping on small scales that is qualitatively similar to that of the well-known Fingers-of-God (FoG) effect due to random velocities, except that the log-normal mock does not include random velocities. This damping is a consequence of non-linearity in the Jacobian, and thus attributing the damping of the power spectrum solely to FoG, as commonly done in the literature, is misleading.

The isotropic gamma-ray background arises from the contribution of unresolved sources, including members of confirmed source classes and proposed gamma-ray emitters such as the radiation induced by dark matter annihilation and decay. Clues about the properties of the contributing sources are imprinted in the anisotropy characteristics of the gamma-ray background. We use 81 months of Pass 7 Reprocessed data from the Fermi Large Area Telescope to perform a measurement of the anisotropy angular power spectrum of the gamma-ray background. We analyze energies between 0.5 and 500 GeV, extending the range considered in the previous measurement based on 22 months of data. We also compute, for the first time, the cross-correlation angular power spectrum between different energy bins. We find that the derived angular spectra are compatible with being Poissonian, i.e. constant in multipole. Moreover, the energy dependence of the anisotropy suggests that the signal is due to two populations of sources, contributing, respectively, below and above 2 GeV. Finally, using data from state-of-the-art numerical simulations to model the dark matter distribution, we constrain the contribution from dark matter annihilation and decay in Galactic and extragalactic structures to the measured anisotropy. These constraints are competitive with those that can be derived from the average intensity of the isotropic gamma-ray background. Data are available at https://www-glast.stanford.edu/pub_data/552.

We show that measurements of the fluctuations in the near-infrared background (NIRB) from the AKARI satellite can be explained by faint galaxy populations at low redshifts. We demonstrate this using reconstructed images from deep galaxy catalogs (HUGS/S-CANDELS) and two independent galaxy population models. In all cases, we find that the NIRB fluctuations measured by AKARI are consistent with faint galaxies and there is no need for a contribution from unknown populations. We find no evidence for a steep Rayleigh-Jeans spectrum for the underlying sources as previously reported. The apparent Rayleigh-Jeans spectrum at large angular scales is likely a consequence of galaxies being removed systematically to deeper levels in the longer wavelength channels.

We calculate the one-point probability density distribution functions (PDF) and the power spectra of the thermal and kinetic Sunyaev-Zeldovich (tSZ and kSZ) effects and the mean Compton Y parameter using the Magneticum Pathfinder simulations, state-of-the-art cosmological hydrodynamical simulations of a large cosmological volume of (896 Mpc/h)^3. These simulations follow in detail the thermal and chemical evolution of the intracluster medium as well as the evolution of super-massive black holes and their associated feedback processes. We construct full-sky maps of tSZ and kSZ from the light-cones out to z=0.17, and one realization of 8.8x8.8 degree wide, deep light-cone out to z=5.2. The local universe at z<0.027 is simulated by a constrained realisation. The tail of the one-point PDF of tSZ from the deep light-cone follows a power-law shape with an index of -3.2. Once convolved with the effective beam of Planck, it agrees with the PDF measured by Planck. The predicted tSZ power spectrum agrees with that of the Planck data at all multipoles up to l~1000, once the calculations are scaled to the Planck 2015 cosmological parameters with \Omega_m=0.308 and \sigma_8=0.8149. Consistent with the results in the literature, however, we continue to find the tSZ power spectrum at l=3000 that is significantly larger than that estimated from the high-resolution ground-based data. The simulation predicts the mean fluctuating Compton Y value of <Y>=1.18x10^{-6} for \Omega_m=0.272 and \sigma_8=0.809. Nearly half (~ 5x10^{-7}) of the signal comes from halos below a virial mass of 10^{13}M_\odot/h. Scaling this to the Planck 2015 parameters, we find <Y>=1.57x10^{-6}. The PDF and the power spectrum of kSZ from our simulation agree broadly with the previous work.

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

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

We present a Bayesian approach to the redshift classification of emission-line galaxies when only a single emission line is detected spectroscopically. We consider the case of surveys for high-redshift Lyman-alpha-emitting galaxies (LAEs), which have traditionally been classified via an inferred rest-frame equivalent width (EW) greater than 20 angstrom. Our Bayesian method relies on known prior probabilities in measured emission-line luminosity functions and equivalent width distributions for the galaxy populations, and returns the probability that an object in question is an LAE given the characteristics observed. This approach will be directly relevant for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX), which seeks to classify ~10^6 emission-line galaxies into LAEs and low-redshift [O II] emitters. For a simulated HETDEX catalog with realistic measurement noise, our Bayesian method recovers 86% of LAEs missed by the traditional EW > 20 angstrom cutoff over 2 < z < 3, outperforming the EW cut in both contamination and incompleteness. This is due to the method's ability to trade off between the two types of binary classification error by adjusting the stringency of the probability requirement for classifying an observed object as an LAE. In our simulations of HETDEX, this method reduces the uncertainty in cosmological distance measurements by 14% with respect to the EW cut, equivalent to recovering 29% more cosmological information. Rather than using binary object labels, this method enables the use of classification probabilities in large-scale structure analyses. It can be applied to narrowband emission-line surveys as well as upcoming large spectroscopic surveys including Euclid and WFIRST.

Strong lensing time-delay systems constrain cosmological parameters via the so-called time-delay distance and the angular diameter distance to the lens. In previous studies, only the former information was used. In this paper, we show that the cosmological constraints improve significantly when the latter information is also included. Specifically, the angular diameter distance plays a crucial role in breaking the degeneracy between the curvature of the Universe and the time-varying equation of state of dark energy. Using a mock sample of 55 bright quadruple lens systems based on expectations for ongoing/future imaging surveys, we find that adding the angular diameter distance information to the time-delay distance information and the cosmic microwave background data of Planck improves the constraint on the constant equation of state by 30%, on the time variation in the equation of state by a factor of two, and on the Hubble constant in the flat $\Lambda$CDM model by a factor of two. Therefore, previous forecasts for the statistical power of time-delay systems were significantly underestimated, i.e., time-delay systems are more powerful than previously appreciated.

Non-thermal pressure in galaxy clusters leads to underestimation of the mass of galaxy clusters based on hydrostatic equilibrium with thermal gas pressure. This occurs even for dynamically relaxed clusters that are used for calibrating the mass-observable scaling relations. We show that the analytical model for non-thermal pressure developed in Shi & Komatsu 2014 can correct for this so-called 'hydrostatic mass bias', if most of the non-thermal pressure comes from bulk and turbulent motions of gas in the intracluster medium. Our correction works for the sample average irrespective of the mass estimation method, or the dynamical state of the clusters. This makes it possible to correct for the bias in the hydrostatic mass estimates from X-ray surface brightness and the Sunyaev-Zel'dovich observations that will be available for clusters in a wide range of redshifts and dynamical states.

Cosmological transverse momentum fields, whose directions are perpendicular to Fourier wave vectors, induce temperature anisotropies in the cosmic microwave background via the kinetic Sunyaev-Zeldovich (kSZ) effect. The transverse momentum power spectrum contains the four-point function of density and velocity fields, $\langle\delta\delta v v\rangle$. In the post-reionization epoch, nonlinear effects dominate in the power spectrum. We use perturbation theory and cosmological $N$-body simulations to calculate this nonlinearity. We derive the next-to-leading order expression for the power spectrum with a particular emphasis on the connected term that has been ignored in the literature. While the contribution from the connected term on small scales ($k>0.1\,h\,\rm{Mpc}^{-1}$) is subdominant relative to the unconnected term, we find that its contribution to the kSZ power spectrum at $\ell = 3000$ at $z<6$ can be as large as ten percent of the unconnected term, which would reduce the allowed contribution from the reionization epoch ($z>6$) by twenty percent. The power spectrum of transverse momentum on large scales is expected to scale as $k^2$ as a consequence of momentum conservation. We show that both the leading and the next-to-leading order terms satisfy this scaling. In particular, we find that both of the unconnected and connected terms are necessary to reproduce $k^2$.

The distance-redshift relation plays a fundamental role in constraining cosmological models. In this paper, we show that measurements of positions and time delays of strongly lensed images of a background galaxy, as well as those of the velocity dispersion and mass profile of a lens galaxy, can be combined to extract the angular diameter distance of the lens galaxy. Physically, as the velocity dispersion and the time delay give a gravitational potential ($GM/r$) and a mass ($GM$) of the lens, respectively, dividing them gives a physical size ($r$) of the lens. Comparing the physical size with the image positions of a lensed galaxy gives the angular diameter distance to the lens. A mismatch between the exact locations at which these measurements are made can be corrected by measuring a local slope of the mass profile. We expand on the original idea put forward by Paraficz and Hjorth, who analyzed singular isothermal lenses, by allowing for an arbitrary slope of a power-law spherical mass density profile, an external convergence, and an anisotropic velocity dispersion. We find that the effect of external convergence cancels out when dividing the time delays and velocity dispersion measurements. We derive a formula for the uncertainty in the angular diameter distance in terms of the uncertainties in the observables. As an application, we use two existing strong lens systems, B1608+656 ($z_{\rm L}=0.6304$) and RXJ1131$-$1231 ($z_{\rm L}=0.295$), to show that the uncertainty in the inferred angular diameter distances is dominated by that in the velocity dispersion, $\sigma^2$, and its anisotropy. We find that the current data on these systems should yield about 16% uncertainty in $D_A$ per object. This improves to 13% when we measure $\sigma^2$ at the so-called sweet-spot radius. Achieving 7% is possible if we can determine $\sigma^2$ with 5% precision.

We report on the first measurement of the three-point function with the position-dependent correlation function from the SDSS-III Baryon Oscillation Spectroscopic Survey (BOSS) Data Release 10 CMASS sample. This new observable measures the correlation between two-point functions of galaxy pairs within different subvolumes, $\hat{\xi}({\rm r},{\rm r}_L)$, where ${\rm r}_L$ is the location of a subvolume, and the corresponding mean overdensities, $\bar{\delta}({\rm r}_L)$. This correlation, which we call the "integrated three-point function", $i\zeta(r)=\langle\hat{\xi}({\rm r},{\rm r}_L)\bar{\delta}({\rm r}_L)\rangle$, measures a three-point function of two short- and one long-wavelength modes, and is generated by nonlinear gravitational evolution and possibly also by the physics of inflation. The $i\zeta(r)$ measured from the BOSS data lies within the scatter of those from the mock galaxy catalogs in redshift space, yielding a ten-percent-level determination of the amplitude of $i\zeta(r)$. The tree-level perturbation theory in redshift space predicts how this amplitude depends on the linear and quadratic nonlinear galaxy bias parameters ($b_1$ and $b_2$), as well as on the amplitude and linear growth rate of matter fluctuations ($\sigma_8$ and $f$). Combining $i\zeta(r)$ with the constraints on $b_1\sigma_8$ and $f\sigma_8$ from the global two-point correlation function and that on $\sigma_8$ from the weak lensing signal of BOSS galaxies, we measure $b_2=0.41\pm0.41$ (68% C.L.) assuming standard perturbation theory at the tree level and the local bias model.

Bayesian model comparison penalizes models with more free parameters that are allowed to vary over a wide range, and thus offers the most robust method to decide whether some given data require new parameters. In this paper, we ask a simple question: do current cosmological data require extensions of the simplest single-field inflation models? Specifically, we calculate the Bayesian evidence of a totally anti-correlated isocurvature perturbation and a running spectral index of the scalar curvature perturbation. These parameters are motivated by recent claims that the observed temperature anisotropy of the cosmic microwave background on large angular scales is too low to be compatible with the simplest inflation models. Both a subdominant, anti-correlated cold dark matter isocurvature component and a negative running index succeed in lowering the large-scale temperature power spectrum. We show that the introduction of isocurvature perturbations is disfavored, whereas that of the running spectral index is only moderately favored, even when the BICEP2 data are included in the analysis without any foreground subtraction.

Damping of magnetic fields via ambipolar diffusion and decay of magnetohydrodynamical (MHD) turbulence in the post decoupling era heats the intergalactic medium (IGM). Delayed recombination of hydrogen atoms in the IGM yields an optical depth to scattering of the cosmic microwave background (CMB). The optical depth generated at $z\gg 10$ does not affect the "reionization bump" of the CMB polarization power spectrum at low multipoles, but affects the temperature and polarization power spectra at high multipoles. Writing the present-day energy density of fields smoothed over the damping scale at the decoupling epoch as $\rho_{B,0}=B_{0}^2/2$, we constrain $B_0$ as a function of the spectral index, $n_B$. Using the Planck 2013 likelihood code that uses the Planck temperature and lensing data together with the WMAP 9-year polarization data, we find the 95% upper bounds of $B_0<0.63$, 0.39, and 0.18~nG for $n_B=-2.9$, $-2.5$, and $-1.5$, respectively. For these spectral indices, the optical depth is dominated by dissipation of the decaying MHD turbulence that occurs shortly after the decoupling epoch. Our limits are stronger than the previous limits ignoring the effects of the fields on ionization history. Inverse Compton scattering of CMB photons off electrons in the heated IGM distorts the thermal spectrum of CMB. Our limits on $B_0$ imply that the $y$-type distortion from dissipation of fields in the post decoupling era should be smaller than $10^{-9}$, $4\times10^{-9}$, and $10^{-9}$, respectively.

We show that in a certain, angle-averaged squeezed limit, the $N$-point function of matter is related to the response of the matter power spectrum to a long-wavelength density perturbation, $P^{-1}d^nP(k|\delta_L)/d\delta_L^n|_{\delta_L=0}$, with $n=N-2$. By performing N-body simulations with a homogeneous overdensity superimposed on a flat Friedmann-Robertson-Lema\^itre-Walker (FRLW) universe using the \emph{separate universe} approach, we obtain measurements of the nonlinear matter power spectrum response up to $n=3$, which is equivalent to measuring the fully nonlinear matter $3-$ to $5-$point function in this squeezed limit. The sub-percent to few percent accuracy of those measurements is unprecedented. We then test the hypothesis that nonlinear $N$-point functions at a given time are a function of the linear power spectrum at that time, which is predicted by standard perturbation theory (SPT) and its variants that are based on the ideal pressureless fluid equations. Specifically, we compare the responses computed from the separate universe simulations and simulations with a rescaled initial (linear) power spectrum amplitude. We find discrepancies of 10\% at $k\simeq 0.2 - 0.5 \,h\,{\rm Mpc}^{-1}$ for $5-$ to $3-$point functions at $z=0$. The discrepancy occurs at higher wavenumbers at $z=2$. Thus, SPT and its variants, carried out to arbitrarily high order, are guaranteed to fail to describe matter $N$-point functions ($N>2$) around that scale.

Turbulent gas motion inside galaxy clusters provides a non-negligible non-thermal pressure support to the intracluster gas. If not corrected, it leads to a systematic bias in the estimation of cluster masses from X-ray and Sunyaev-Zel'dovich (SZ) observations assuming hydrostatic equilibrium, and affects interpretation of measurements of the SZ power spectrum and observations of cluster outskirts from ongoing and upcoming large cluster surveys. Recently, Shi & Komatsu developed an analytical model for predicting the radius, mass, and redshift dependence of the non-thermal pressure contributed by the kinetic random motions of intracluster gas sourced by the cluster mass growth. In this paper, we compare the predictions of this analytical model to a state-of-the-art cosmological hydrodynamics simulation. As different mass growth histories result in different non-thermal pressure, we perform the comparison on 65 simulated galaxy clusters on a cluster-by-cluster basis. We find an excellent agreement between the modelled and simulated non-thermal pressure profiles. Our results open up the possibility of using the analytical model to correct the systematic bias in the mass estimation of galaxy clusters. We also discuss tests of the physical picture underlying the evolution of intracluster non-thermal gas motions, as well as a way to further improve the analytical modeling, which may help achieve a unified understanding of non-thermal phenomena in galaxy clusters.

Temperature anisotropy of the cosmic microwave background offers a test of the fundamental symmetry of spacetime during cosmic inflation. Violation of rotational symmetry yields a distinct signature in the power spectrum of primordial fluctuations as $P({\mathbf k})=P_0(k)[1+g_*(\hat{\mathbf k}\cdot\hat{\mathbf E}_{\rm cl})^2]$, where $\hat{\mathbf E}_{\rm cl}$ is a preferred direction in space and $g_*$ is an amplitude. Using the \textit{Planck} 2013 temperature maps, we find no evidence for violation of rotational symmetry, $g_*=0.002\pm 0.016$ (68% CL), once the known effects of asymmetry of the \textit{Planck} beams and Galactic foreground emission are removed.

The large-scale statistics of observables such as the galaxy density are chiefly determined by their dependence on the local coarse-grained matter density. This dependence can be measured directly and efficiently in N-body simulations by using the fact that a uniform density perturbation with respect to some fiducial background cosmology is equivalent to modifying the background and including curvature, i.e., by simulating a "separate universe". We derive this mapping to fully non-linear order, and provide a step-by-step description of how to perform and analyse the separate universe simulations. This technique can be applied to a wide range of observables. As an example, we calculate the response of the non-linear matter power spectrum to long-wavelength density perturbations, which corresponds to the angle-averaged squeezed limit of the matter bispectrum and higher $n$-point functions. Using only a modest simulation volume, we obtain results with percent-level precision over a wide range of scales.

Cross-correlating gamma-ray maps with locations of galaxies in the low-redshift Universe vastly increases sensitivity to signatures of annihilation of dark matter particles. Low-redshift galaxies are ideal targets, as the largest contribution to anisotropy in the gamma-ray sky from annihilation comes from $z\lesssim 0.1$, where we expect minimal contributions from astrophysical sources such as blazars. Cross-correlating the five-year data of Fermi-LAT with the redshift catalog of the 2MASS survey can detect gamma rays from annihilation if dark matter has the canonical annihilation cross section and its mass is smaller than $\sim$100 GeV.

Non-thermal pressure in the intracluster gas has been found ubiquitously in numerical simulations, and observed indirectly. In this paper we develop an analytical model for intracluster non-thermal pressure in the virial region of relaxed clusters. We write down and solve a first-order differential equation describing the evolution of non-thermal velocity dispersion. This equation is based on insights gained from observations, numerical simulations, and theory of turbulence. The non-thermal energy is sourced, in a self-similar fashion, by the mass growth of clusters via mergers and accretion, and dissipates with a time-scale determined by the turnover time of the largest turbulence eddies. Our model predicts a radial profile of non-thermal pressure for relaxed clusters. The non-thermal fraction increases with radius, redshift, and cluster mass, in agreement with numerical simulations. The radial dependence is due to a rapid increase of the dissipation time-scale with radii, and the mass and redshift dependence comes from the mass growth history. Combing our model for the non-thermal fraction with the Komatsu-Seljak model for the total pressure, we obtain thermal pressure profiles, and compute the hydrostatic mass bias. We find typically 10% bias for the hydrostatic mass enclosed within $r_{500}$.

Fluctuations in the intensity and polarization of the cosmic microwave background (CMB) and the large-scale distribution of matter in the universe each contain clues about the nature of the earliest moments of time. The next generation of CMB and large-scale structure (LSS) experiments are poised to test the leading paradigm for these earliest moments---the theory of cosmic inflation---and to detect the imprints of the inflationary epoch, thereby dramatically increasing our understanding of fundamental physics and the early universe. A future CMB experiment with sufficient angular resolution and frequency coverage that surveys at least 1% of the sky to a depth of 1 uK-arcmin can deliver a constraint on the tensor-to-scalar ratio that will either result in a 5-sigma measurement of the energy scale of inflation or rule out all large-field inflation models, even in the presence of foregrounds and the gravitational lensing B-mode signal. LSS experiments, particularly spectroscopic surveys such as the Dark Energy Spectroscopic Instrument, will complement the CMB effort by improving current constraints on running of the spectral index by up to a factor of four, improving constraints on curvature by a factor of ten, and providing non-Gaussianity constraints that are competitive with the current CMB bounds.

This is a report on the status and prospects of the quantification of neutrino properties through the cosmological neutrino background for the Cosmic Frontier of the Division of Particles and Fields Community Summer Study long-term planning exercise. Experiments planned and underway are prepared to study the cosmological neutrino background in detail via its influence on distance-redshift relations and the growth of structure. The program for the next decade described in this document, including upcoming spectroscopic galaxy surveys eBOSS and DESI and a new Stage-IV CMB polarization experiment CMB-S4, will achieve sigma(sum m_nu) = 16 meV and sigma(N_eff) = 0.020. Such a mass measurement will produce a high significance detection of non-zero sum m_nu, whose lower bound derived from atmospheric and solar neutrino oscillation data is about 58 meV. If neutrinos have a minimal normal mass hierarchy, this measurement will definitively rule out the inverted neutrino mass hierarchy, shedding light on one of the most puzzling aspects of the Standard Model of particle physics --- the origin of mass. This precise a measurement of N_eff will allow for high sensitivity to any light and dark degrees of freedom produced in the big bang and a precision test of the standard cosmological model prediction that N_eff = 3.046.

The Wilkinson Microwave Anisotropy Probe (WMAP) mapped the distribution of temperature and polarization over the entire sky in five microwave frequency bands. These full-sky maps were used to obtain measurements of temperature and polarization anisotropy of the cosmic microwave background with the unprecedented accuracy and precision. The analysis of two-point correlation functions of temperature and polarization data gives determinations of the fundamental cosmological parameters such as the age and composition of the universe, as well as the key parameters describing the physics of inflation, which is further constrained by three-point correlation functions. WMAP observations alone reduced the flat $\Lambda$ cold dark matter ($\Lambda$CDM) cosmological model (six) parameter volume by a factor of >68,000 compared with pre-WMAP measurements. The WMAP observations (sometimes in combination with other astrophysical probes) convincingly show the existence of non-baryonic dark matter, the cosmic neutrino background, flatness of spatial geometry of the universe, a deviation from a scale-invariant spectrum of initial scalar fluctuations, and that the current universe is undergoing an accelerated expansion. The WMAP observations provide the strongest ever support for inflation; namely, the structures we see in the universe originate from quantum fluctuations generated during inflation.

The influence of large-scale density fluctuations on structure formation on small scales is described by the three-point correlation function (bispectrum) in the so-called "squeezed configurations," in which one wavenumber, say $k_3$, is much smaller than the other two, i.e., $k_3\ll k_1\approx k_2$. This bispectrum is generated by non-linear gravitational evolution and possibly also by inflationary physics. In this paper, we use this fact to show that the bispectrum in the squeezed configurations can be measured without employing three-point function estimators. Specifically, we use the "position-dependent power spectrum," i.e., the power spectrum measured in smaller subvolumes of the survey (or simulation box), and correlate it with the mean overdensity of the corresponding subvolume. This correlation directly measures an integral of the bispectrum dominated by the squeezed configurations. Measuring this correlation is only slightly more complex than measuring the power spectrum itself, and sidesteps the considerable complexity of the full bispectrum estimation. We use cosmological $N$-body simulations of collisionless particles with Gaussian initial conditions to show that the measured correlation between the position-dependent power spectrum and the long-wavelength overdensity agrees with the theoretical expectation. The position-dependent power spectrum thus provides a new, efficient, and promising way to measure the squeezed-limit bispectrum from large-scale structure observations such as galaxy redshift surveys.

Soft limits of $N$-point correlation functions, in which one wavenumber is much smaller than the others, play a special role in constraining the physics of inflation. Anisotropic sources such as a vector field during inflation generate distinct angular dependence in all these correlators. In this paper we focus on the four-point correlator (the trispectrum $T$). We adopt a parametrization motivated by models in which the inflaton $\phi$ is coupled to a vector field through a $I^2 \left( \phi \right) F^2$ interaction, namely $T_{\zeta}({\bf k}_1, {\bf k}_2, {\bf k}_3, {\bf k}_4) \equiv \sum_n d_n [ P_n(\hat{\bf k}_1 \cdot \hat{\bf k}_3) + P_n(\hat{\bf k}_1 \cdot \hat{\bf k}_{12}) + P_n(\hat{\bf k}_3 \cdot \hat{\bf k}_{12}) ] P_{\zeta}(k_1) P_{\zeta}(k_3) P_\zeta(k_{12}) + (23~{\rm perm})$, where $P_n$ denotes the Legendre polynomials. This shape is enhanced when the wavenumbers of the diagonals of the quadrilateral are much smaller than the sides, ${\bf k}_i$. The coefficient of the isotropic part, $d_0$, is equal to $\tau_{\rm NL}/6$ discussed in the literature. A $I^2 \left( \phi \right) F^2$ interaction generates $d_2 = 2 d_0$ which is, in turn, related to the quadrupole modulation parameter of the power spectrum, $g_*$, as $d_2 \approx 14 |g_*| N^2$ with $N \approx 60$. We show that $d_0$ and $d_2$ can be equally well-constrained: the expected $68 \%$ CL error bars on these coefficients from a cosmic-variance-limited experiment measuring temperature anisotropy of the cosmic microwave background up to $\ell_{\rm max}=2000$ are $\delta d_2 \approx 4 \delta d_0 = 105$. Therefore, we can reach $|g_*|=10^{-3}$ by measuring the angle-dependent trispectrum. The current upper limit on $\tau_{\rm NL}$ from the ${\it Planck}$ temperature maps yields $|g_*|<0.02$ ($95 \%$ CL).

Survey observations of the three-dimensional locations of galaxies are a powerful approach to measure the distribution of matter in the universe, which can be used to learn about the nature of dark energy, physics of inflation, neutrino masses, etc. A competitive survey, however, requires a large volume (e.g., Vsurvey is roughly 10 Gpc3) to be covered, and thus tends to be expensive. A "sparse sampling" method offers a more affordable solution to this problem: within a survey footprint covering a given survey volume, Vsurvey, we observe only a fraction of the volume. The distribution of observed regions should be chosen such that their separation is smaller than the length scale corresponding to the wavenumber of interest. Then one can recover the power spectrum of galaxies with precision expected for a survey covering a volume of Vsurvey (rather than the volume of the sum of observed regions) with the number density of galaxies given by the total number of observed galaxies divided by Vsurvey (rather than the number density of galaxies within an observed region). We find that regularly-spaced sampling yields an unbiased power spectrum with no window function effect, and deviations from regularly-spaced sampling, which are unavoidable in realistic surveys, introduce calculable window function effects and increase the uncertainties of the recovered power spectrum. While we discuss the sparse sampling method within the context of the forthcoming Hobby-Eberly Telescope Dark Energy Experiment, the method is general and can be applied to other galaxy surveys.

LiteBIRD is a next-generation satellite mission to measure the polarization of the cosmic microwave background (CMB) radiation. On large angular scales the B-mode polarization of the CMB carries the imprint of primordial gravitational waves, and its precise measurement would provide a powerful probe of the epoch of inflation. The goal of LiteBIRD is to achieve a measurement of the characterizing tensor to scalar ratio $r$ to an uncertainty of $\delta r=0.001$. In order to achieve this goal we will employ a kilo-pixel superconducting detector array on a cryogenically cooled sub-Kelvin focal plane with an optical system at a temperature of 4~K. We are currently considering two detector array options; transition edge sensor (TES) bolometers and microwave kinetic inductance detectors (MKID). In this paper we give an overview of LiteBIRD and describe a TES-based polarimeter designed to achieve the target sensitivity of 2~$\mu$K$\cdot$arcmin over the frequency range 50 to 320~GHz.

Observational limits on $y$- and $\mu$-type distortions can constrain properties of magnetic fields in the early universe. For a Gaussian, random, and non-helical field, $\mu$ and $y$ are calculated a function of the present-day strength of the field, $B_0$, smoothed over a certain Gaussian width, $k_c^{-1}$, and spectral index, $n_B$, defined by $P_B(k)\propto k^{n_B}$. For $n_B=-2.9$ and $k_c^{-1}=1 {\rm Mpc}$, the COBE/FIRAS limit on $\mu$ yields $B_0<40$ nG, whereas the projected PIXIE limit on $\mu$ would yield $B_0<0.8$ nG. For non-scale-invariant spectra, constraints can be stronger. For $B_0=1$ nG with $k_c^{-1}=1 {\rm Mpc}$, the COBE/FIRAS limit on $\mu$ excludes a wide range of spectral indices given by $n_B>-2.6$. After decoupling, energy dissipation is due to ambipolar diffusion and decaying MHD turbulence, creating a $y$-type distortion. The distortion is completely dominated by decaying MHD turbulence, and is of order $y\approx 10^{-7}$ for a few nG field smoothed over the damping scale at the decoupling epoch, $k_{d,dec}\approx 290 (B_0/1 {\rm nG})^{-1} {\rm Mpc}^{-1}$. This contribution is as large as those of the known contributions such as reionization and virialized objects at lower redshifts. The projected PIXIE limit on $y$ would exclude $B_0>1.0$ and 0.6 nG for $n_B=-2.9$ and -2.3, respectively, and $B_0>0.6$ nG for $n_B\geq 2$. The current limits on the optical depth to Thomson scattering restrict the predicted $y$-type distortion to be $y\lesssim 10^{-8}$. (Abridged)

We calculate the bispectrum of primordial curvature perturbations, \zeta, generated during "open inflation." Inflation occurs inside a bubble nucleated via quantum tunneling from the background false vacuum state. Our universe lives inside the bubble, which can be described as a Friedman-Lema\^itre-Robertson-Walker (FLRW) universe with negative spatial curvature, undergoing slow-roll inflation. We pay special attention to the issue of an initial state for quantum fluctuations. A "vacuum state" defined by a positive-frequency mode in de Sitter space charted by open coordinates is different from the Euclidean vacuum (which is equivalent to the so-called "Bunch-Davies vacuum"). Quantum tunneling then modifies the initial state away from the original Euclidean vacuum. While most of the previous study on modifications of the initial quantum state introduces, by hand, an initial time at which the quantum state is modified as well as the form of the modification, an effective initial time naturally emerges and the form is fixed by quantum tunneling in open inflation models. Therefore, open inflation enables a self-consistent computation of the effect of a modified initial state on the bispectrum. We find a term which goes as <\zeta_{k_1} \zeta_{k_2} \zeta_{k_3}> \propto 1/k_1^2k_3^4 in the so-called squeezed configurations, k_3\ll k_1\approx k_2, in agreement with the previous study. The bispectrum in the exact folded limit, e.g., k_1=k_2+k_3, is also enhanced and remains finite. However, these terms are exponentially suppressed when the wavelength of \zeta is smaller than the curvature radius of the universe. The leading-order bispectrum is equal to the usual one from single-field slow-roll inflation; the terms specific for open inflation arise only in the sub-leading order when the wavelength of \zeta is smaller than the curvature radius.

In non-attractor single-field inflation models producing a scale-invariant power spectrum, the curvature perturbation on super-horizon scales grows as ${\cal R}\propto a^3$. This is so far the only known class of self-consistent single-field models with a Bunch-Davies initial state that can produce a large squeezed-limit bispectrum violating Maldacena's consistency relation. Given the importance of this result, we calculate the bispectrum with three different methods: using quantum field theory calculations in two different gauges, and classical calculations (the $\delta N$ formalism). All the results agree, giving the local-form bispectrum parameter of $f^{local}_{NL}=5(1+c_s^2)/(4c_s^2)$. This result is valid for arbitrary values of the speed of sound parameter, $c_s$, for a particular non-attractor model we consider in this paper.

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

We present the final nine-year maps and basic results from the WMAP mission. We provide new nine-year full sky temperature maps that were processed to reduce the asymmetry of the effective beams. Temperature and polarization sky maps are examined to separate CMB anisotropy from foreground emission, and both types of signals are analyzed in detail. The WMAP mission has resulted in a highly constrained LCDM cosmological model with precise and accurate parameters in agreement with a host of other cosmological measurements. When WMAP data are combined with finer scale CMB, baryon acoustic oscillation, and Hubble constant measurements, we find that Big Bang nucleosynthesis is well supported and there is no compelling evidence for a non-standard number of neutrino species (3.84+/-0.40). The model fit also implies that the age of the universe is 13.772+/-0.059 Gyr, and the fit Hubble constant is H0 = 69.32+/-0.80 km/s/Mpc. Inflation is also supported: the fluctuations are adiabatic, with Gaussian random phases; the detection of a deviation of the scalar spectral index from unity reported earlier by WMAP now has high statistical significance (n_s = 0.9608+/-0.0080); and the universe is close to flat/Euclidean, Omega_k = -0.0027 (+0.0039/-0.0038). Overall, the WMAP mission has resulted in a reduction of the cosmological parameter volume by a factor of 68,000 for the standard six-parameter LCDM model, based on CMB data alone. For a model including tensors, the allowed seven-parameter volume has been reduced by a factor 117,000. Other cosmological observations are in accord with the CMB predictions, and the combined data reduces the cosmological parameter volume even further. With no significant anomalies and an adequate goodness-of-fit, the inflationary flat LCDM model and its precise and accurate parameters rooted in WMAP data stands as the standard model of cosmology.

We present cosmological parameter constraints based on the final nine-year WMAP data, in conjunction with additional cosmological data sets. The WMAP data alone, and in combination, continue to be remarkably well fit by a six-parameter LCDM model. When WMAP data are combined with measurements of the high-l CMB anisotropy, the BAO scale, and the Hubble constant, the densities, Omegabh2, Omegach2, and Omega_L, are each determined to a precision of ~1.5%. The amplitude of the primordial spectrum is measured to within 3%, and there is now evidence for a tilt in the primordial spectrum at the 5sigma level, confirming the first detection of tilt based on the five-year WMAP data. At the end of the WMAP mission, the nine-year data decrease the allowable volume of the six-dimensional LCDM parameter space by a factor of 68,000 relative to pre-WMAP measurements. We investigate a number of data combinations and show that their LCDM parameter fits are consistent. New limits on deviations from the six-parameter model are presented, for example: the fractional contribution of tensor modes is limited to r<0.13 (95% CL); the spatial curvature parameter is limited to -0.0027 (+0.0039/-0.0038); the summed mass of neutrinos is <0.44 eV (95% CL); and the number of relativistic species is found to be 3.84+/-0.40 when the full data are analyzed. The joint constraint on Neff and the primordial helium abundance agrees with the prediction of standard Big Bang nucleosynthesis. We compare recent PLANCK measurements of the Sunyaev-Zel'dovich effect with our seven-year measurements, and show their mutual agreement. Our analysis of the polarization pattern around temperature extrema is updated. This confirms a fundamental prediction of the standard cosmological model and provides a striking illustration of acoustic oscillations and adiabatic initial conditions in the early universe.

The bispectrum of primordial curvature perturbations in the squeezed configuration, in which one wavenumber, $k_3$, is much smaller than the other two, $k_3\ll k_1\approx k_2$, plays a special role in constraining the physics of inflation. In this paper we study a new phenomenological signature in the squeezed-limit bispectrum: namely, the amplitude of the squeezed-limit bispectrum depends on an angle between ${\bf k}_1$ and ${\bf k}_3$ such that $B_\zeta(k_1, k_2, k_3) \to 2 \sum_L c_L P_L(\hat{\bf k}_1 \cdot \hat{\bf k}_3) P_\zeta(k_1)P_\zeta(k_3)$, where $P_L$ are the Legendre polynomials. While $c_0$ is related to the usual local-form $f_{\rm NL}$ parameter as $c_0=6f_{\rm NL}/5$, the higher-multipole coefficients, $c_1$, $c_2$, etc., have not been constrained by the data. Primordial curvature perturbations sourced by large-scale magnetic fields generate non-vanishing $c_0$, $c_1$, and $c_2$. Inflation models whose action contains a term like $I(\phi)^2 F^2$ generate $c_2=c_0/2$. A recently proposed "solid inflation" model generates $c_2\gg c_0$. A cosmic-variance-limited experiment measuring temperature anisotropy of the cosmic microwave background up to $\ell_{\rm max}=2000$ is able to measure these coefficients down to $\delta c_0=4.4$, $\delta c_1=61$, and $\delta c_2=13$ (68% CL). We also find that $c_0$ and $c_1$, and $c_0$ and $c_2$, are nearly uncorrelated. Measurements of these coefficients will open up a new window into the physics of inflation such as the existence of vector fields during inflation or non-trivial symmetry structure of inflaton fields. Finally, we show that the original form of the Suyama-Yamaguchi inequality does not apply to the case involving higher-spin fields, but a generalized form does.

The Integrated Sachs-Wolfe (ISW) effect produces a secondary temperature anisotropy of CMB. The main contribution comes from z<2, where dark energy leads to a decay of potentials. As the same photons are gravitationally lensed by these decaying potentials, there exists a high degree of correlation between the ISW effect and CMB lensing, leading to a non-zero three-point correlation (bispectrum) of the observed temperature anisotropy. This ISW-lensing bispectrum, whose shape resembles that of the so-called "local-form" primordial bispectrum parametrized by fNL, is known to be the largest contamination of fNL. In order to avoid a spurious detection of primordial non-Gaussianity, we need to remove the ISW-lensing bias. In this work, we investigate three debiasing methods: (I) subtraction of an expected, ensemble average of the ISW-lensing bispectrum; (II) subtraction of a measured ISW-lensing bispectrum; and (III) direct subtraction of an estimated ISW signal from an observed temperature map. One may use an estimation of the ISW map from external non-CMB data or that from the CMB data themselves. As the methods II and III are based on fewer assumptions about the nature of dark energy, they are preferred over the method I. While the methods I and II yield unbiased estimates of fNL with comparable error bars, the method III yields a biased result when the underlying primordial fNL is non-zero and the ISW map is estimated from a lensing potential reconstructed from the observed temperature map. One of the sources of the bias is a lensing reconstruction noise bias which is independent of fNL and can be calculated precisely, but other fNL-dependent terms are difficult to compute reliably. We thus conclude that the method II is the best, model-independent way to remove the ISW-lensing bias of fNL, enabling us to test the physics of inflation with smaller systematic errors.

The detailed origin of the diffuse gamma-ray background is still unknown. However, the contribution of unresolved sources is expected to induce small-scale anisotropies in this emission, which may provide a way to identify and constrain the properties of its contributors. Recent studies have predicted the contributions to the angular power spectrum (APS) from extragalactic and galactic dark matter (DM) annihilation or decay. The Fermi-LAT collaboration reported detection of angular power with a significance larger than $3\sigma$ in the energy range from 1 GeV to 10 GeV on 22 months of data [Ackermann et al. 2012]. For these preliminary results the already published Fermi-LAT APS measurements [Ackermann et al. 2012] are compared to the accurate predictions for DM anisotropies from state-of-the-art cosmological simulations as presented in [Fornasa et al. 2013] to derive constraints on different DM candidates.

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

Annihilation of dark matter particles in cosmological halos (including a halo of the Milky Way) contributes to the diffuse gamma-ray background (DGRB). As this contribution will appear anisotropic in the sky, one can use the angular power spectrum of anisotropies in DGRB to constrain properties of dark matter particles. By comparing the updated analytic model of the angular power spectrum of DGRB from dark matter annihilation with the power spectrum recently measured from the 22-month data of Fermi Large Area Telescope (LAT), we place upper limits on the annihilation cross section of dark matter particles as a function of dark matter masses. We find that the current data exclude <\sigma v> >~ 10^{-25} cm^3 s^{-1} for annihilation into b\bar{b} at the dark matter mass of 10 GeV, which is a factor of three times larger than the canonical cross section. The limits are weaker for larger dark matter masses. The limits can be improved further with more Fermi-LAT data as well as by using the power spectrum at lower multipoles (l <~ 150), which are currently not used due to a potential Galactic foreground contamination.

We calculate the angular power spectrum of the Cosmic Microwave Background (CMB) temperature fluctuations induced by the kinetic Sunyaev-Zel'dovich (kSZ) effect from the epoch of reionization (EOR). We use detailed N-body+radiative transfer simulations to follow inhomogeneous reionization of the intergalactic medium (IGM). For the first time we take into account the "self-regulation" of reionization: star formation in low-mass dwarf galaxies (10^8 M_\sun \lesssim M \lesssim 10^9 M_\sun) or minihalos (10^5 M_\sun \lesssim M \lesssim 10^8 M_\sun) is suppressed if these halos form in the regions that were already ionized or Lyman-Werner dissociated. Some previous work suggested that the amplitude of the kSZ power spectrum from the EOR can be described by a two-parameter family: the epoch of half ionization and the duration of reionization. However, we argue that this picture applies only to simple forms of the reionization history which are roughly symmetric about the half-ionization epoch. In self-regulated reionization, the universe begins to be ionized early, maintains a low level of ionization for an extended period, and then finishes reionization as soon as high-mass atomically-cooling halos dominate. While inclusion of self-regulation affects the amplitude of the kSZ power spectrum only modestly (\sim 10 %), it can change the duration of reionization by a factor of more than two. We conclude that the simple two-parameter family does not capture the effect of a physical, yet complex, reionization history caused by self-regulation. When added to the post-reionization kSZ contribution, our prediction for the total kSZ power spectrum is below the current upper bound from the South Pole Telescope. Therefore, the current upper bound on the kSZ effect from the EOR is consistent with our understanding of the physics of reionization.

Large surveys for Lyman-alpha emitting (LAE) galaxies have been proposed as a new method for measuring clustering of the galaxy population at high redshift with the goal of determining cosmological parameters. However, Lyman-alpha radiative transfer effects may modify the observed clustering of LAE galaxies in a way that mimics gravitational effects, potentially reducing the precision of cosmological constraints. For example, the effect of the linear redshift-space distortion on the power spectrum of LAE galaxies is potentially degenerate with Lyman-alpha radiative transfer effects owing to the dependence of observed flux on intergalactic medium velocity gradients. In this paper, we show that the three-point function (bispectrum) can distinguish between gravitational and non-gravitational effects, and thus breaks these degeneracies, making it possible to recover cosmological parameters from LAE galaxy surveys. Constraints on the angular diameter distance and the Hubble expansion rate can also be improved by combining power spectrum and bispectrum measurements.

Large-scale reionization simulations are described which combine the results of cosmological N-body simulations that model the evolving density and velocity fields and identify the galactic halo sources, with ray-tracing radiative transfer calculations which model the nonequilibrium ionization of the intergalactic medium. These simulations have been used to predict some of the signature effects of reionization on cosmic radiation backgrounds, including the CMB, near-IR, and redshifted 21cm backgrounds. We summarize some of our recent progress in this work, and address the question of whether observations of such signature effects can be used to distinguish the relative contributions of galaxies of different masses to reionization.

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

We place new constraints on the contribution of blazars to the large-scale isotropic gamma-ray background (IGRB) by jointly analyzing the measured source count distribution (logN-logS) of blazars and the measured intensity and anisotropy of the IGRB. We find that these measurements point to a consistent scenario in which unresolved blazars make less than 20% of the IGRB intensity at 1-10 GeV while accounting for the majority of the measured anisotropy in that energy band. These results indicate that the remaining fraction of the IGRB intensity is made by a component with a low level of intrinsic anisotropy. We determine upper limits on the anisotropy from non-blazar sources, adopting the best-fit parameters of the measured source count distribution to calculate the unresolved blazar anisotropy. In addition, we show that the anisotropy measurement excludes some recently proposed models of the unresolved blazar population.

We investigate the phenomenological consequences of a modification of the initial state of a single inflationary field. While single-field inflation with the standard Bunch-Davies initial vacuum state does not generally produce a measurable three-point function (bispectrum) in the squeezed configuration, allowing for a non-standard initial state produces an exception. Here, we calculate the signature of an initial state modification in single-field slow-roll inflation in both the scale-dependent bias of the large-scale structure (LSS) and mu-type distortion in the black-body spectrum of the cosmic microwave background (CMB). We parametrize the initial state modifications and identify certain choices of parameters as natural, though we also note some fine-tuned choices that can yield a larger bispectrum. In both cases, we observe a distinctive k^-3 signature in LSS (as opposed to k^-2 for the local-form). As a non-zero bispectrum in the squeezed configuration correlates a long-wavelength mode with two short-wavelength modes, it induces a correlation between the CMB temperature anisotropy on large scales with the temperature-anisotropy-squared on very small scales; this correlation persists as the small-scale anisotropy-squared is processed into mu-type distortions. While the local-form mu-distortion turns out to be too small to detect in the near future, a modified initial vacuum state enhances the signal by a large factor owing to an extra factor of k_1/k. For example, a proposed absolutely-calibrated experiment, PIXIE, is expected to detect this correlation with a signal-to-noise ratio greater than 10, for an occupation number of about 0.5 in the observable modes. Relatively calibrated experiments such as Planck and LiteBIRD should also be able to measure this effect, provided that the relative calibration between different frequencies meets the required precision. (Abridged)

The bispectrum of the cosmic microwave background (CMB) generated by a correlation between a time-dependent gravitational potential and the weak gravitational lensing effect provides a direct measurement of the influence of dark energy on CMB. This bispectrum is also known to yield the most important contamination of the so-called "local-form" primordial bispectrum, which can be used to rule out all single-field inflation models. In this paper, we reexamine the effect of non-linear matter clustering on this bispectrum. We compare three different approaches: the 3rd-order perturbation theory (3PT), and two empirical fitting formulae available in the literature, finding that detailed modeling of non-linearity appears to be not very important, as most of the signal-to-noise comes from the squeezed triangle, for which the correlation in the linear regime dominates. The expected signal-to-noise ratio for an experiment dominated by the cosmic variance up to $l_{\rm max}=1500$ is about 5, which is much smaller than the previous estimates including non-linearity, but agrees with the estimates based on the linear calculation. We find that the difference between the linear and non-linear predictions is undetectable, and does not alter the contamination of the local-form primordial non-Gaussianity.

Current observations suggest that the universe was reionized sometime before z~6. One way to observe this epoch of the universe is through the Near Infrared Background (NIRB), which contains information about galaxies which may be too faint to be observed individually. We calculate the angular power spectrum (C_l) of the NIRB fluctuations caused by the distribution of these galaxies. Assuming a complete subtraction of any post-reionization component, C_l will be dominated by galaxies responsible for completing reionization (e.g., z~6). The shape of C_l at high l is sensitive to the amount of non-linear bias of dark matter halos hosting galaxies. As the non-linear bias depends on the mass of these halos, we can use the shape of C_l to infer typical masses of dark matter halos responsible for completing reionization. We extend our previous study by using a higher-resolution N-body simulation, which can resolve halos down to 10^8 M_sun. We also include improved radiative transfer, which allows for the suppression of star formation in small-mass halos due to photo-ionization heating. As the non-linear bias enhances the dark-matter-halo power spectrum on small scales, we find that C_l is steeper for the case with a complete suppression of small sources or partial suppression of star formation in small halos (the minimum galaxy mass is M_min=10^9 M_sun in ionized regions and M_min=10^8 M_sun in neutral regions) than the case in which these small halos were unsuppressed. In all cases, we do not see a turn-over toward high l in the shape of l^2 C_l.

The contribution of unresolved sources to the diffuse gamma-ray background could induce anisotropies in this emission on small angular scales. We analyze the angular power spectrum of the diffuse emission measured by the Fermi LAT at Galactic latitudes |b| > 30 deg in four energy bins spanning 1 to 50 GeV. At multipoles \ell \ge 155, corresponding to angular scales \lesssim 2 deg, angular power above the photon noise level is detected at >99.99% CL in the 1-2 GeV, 2-5 GeV, and 5-10 GeV energy bins, and at >99% CL at 10-50 GeV. Within each energy bin the measured angular power takes approximately the same value at all multipoles \ell \ge 155, suggesting that it originates from the contribution of one or more unclustered source populations. The amplitude of the angular power normalized to the mean intensity in each energy bin is consistent with a constant value at all energies, C_P/<I>^2 = 9.05 +/- 0.84 x 10^{-6} sr, while the energy dependence of C_P is consistent with the anisotropy arising from one or more source populations with power-law photon spectra with spectral index \Gamma_s = 2.40 +/- 0.07. We discuss the implications of the measured angular power for gamma-ray source populations that may provide a contribution to the diffuse gamma-ray background.

The small angular scale fluctuations of the (on large scale) isotropic gamma-ray background (IGRB) carry information about the presence of unresolved source classes. A guaranteed contribution to the IGRB is expected from the unresolved gamma-ray AGN while other extragalactic sources, Galactic gamma-ray source populations and dark matter Galactic and extragalactic structures (and sub-structures) are candidate contributors. The IGRB was measured with unprecedented precision by the Large Area Telescope (LAT) on-board of the Fermi gamma-ray observatory, and these data were used for measuring the IGRB angular power spectrum (APS). Detailed Monte Carlo simulations of Fermi-LAT all-sky observations were performed to provide a reference against which to compare the results obtained for the real data set. The Monte Carlo simulations are also a method for performing those detailed studies of the APS contributions of single source populations, which are required in order to identify the actual IGRB contributors. We present preliminary results of an anisotropy search in the IGRB. At angular scales <2deg (e.g. above multipole 155), angular power above the photon noise level is detected, at energies between 1 and 10 GeV in each energy bin, with statistical significance between 7.2 and 4.1 sigmas. The energy not dependence of the fluctuation anisotropy is pointing to the presence of one or more unclustered source populations, while the energy dependence of the intensity anisotropy is consistent with source populations having average photon index 2.40\pm0.07.

We reconsider the pixel-based, "template" polarized foreground removal method within the context of a next-generation, low-noise, low-resolution (0.5 degree FWHM) space-borne experiment measuring the cosmological B-mode polarization signal in the cosmic microwave background (CMB). This method was put forward by the Wilkinson Microwave Anisotropy Probe (WMAP) team and further studied by Efstathiou et al. We need at least 3 frequency channels: one is used for extracting the CMB signal, whereas the other two are used to estimate the spatial distribution of the polarized dust and synchrotron emission. No external template maps are used. We extract the tensor-to-scalar ratio (r) from simulated sky maps consisting of CMB, noise (2 micro K arcmin), and a foreground model, and find that, even for the simplest 3-frequency configuration with 60, 100, and 240 GHz, the residual bias in r is as small as Delta r~0.002. This bias is dominated by the residual synchrotron emission due to spatial variations of the synchrotron spectral index. With an extended mask with fsky=0.5, the bias is reduced further down to <0.001.

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

A long standing problem in weak lensing is about how to construct cosmic shear estimators from galaxy images. Conventional methods average over a single quantity per galaxy to estimate each shear component. We show that any such shear estimators must reduce to a highly nonlinear form when the galaxy image is described by three parameters (pure ellipse), even in the absence of the point spread function (PSF). In the presence of the PSF, we argue that this class of shear estimators do not likely exist. Alternatively, we propose a new way of measuring the cosmic shear: instead of averaging over a single value from each galaxy, we average over two numbers, and then take the ratio to estimate the shear component. In particular, the two numbers correspond to the numerator and denominators which generate the quadrupole moments of the galaxy image in Fourier space, as proposed in Zhang (2008). This yields a statistically unbiased estimate of the shear component. Consequently, measurements of the n-point spatial correlations of the shear fields should also be modified: one needs to take the ratio of two correlation functions to get the desired, unbiased shear correlation.

We study the escape of Ly-alpha photons from Ly-alpha emitting galaxies (LAEs) and the overall galaxy population using a sample of 99 LAEs at 1.9<z<3.8 detected through integral-field spectroscopy of blank fields by the HETDEX Pilot Survey. For 89 LAEs showing counterparts in deep broad-band images we measure the rest-frame UV luminosity and the UV slope, which we use to estimate E(B-V) under the assumption of a constant intrinsic UV slope for LAEs. These two quantities are used to measure the dust-corrected star formation rate (SFR). A comparison between the observed Ly-alpha luminosity and that predicted by the dust-corrected SFR yields the Ly-alpha escape fraction. We also measure the Ly-alpha luminosity function. Integration of the luminosity function provides a measurement of the Ly-alpha luminosity density across our redshift range. We combine our data with that from other surveys at 0.3<z<7.7 to trace the evolution of the Ly-alpha luminosity density. We then compare it to that expected from the star-formation history of the universe in order to characterize the evolution of the Ly-alpha escape fraction of the overall galaxy population [abriged]

(Abridged) A simple six-parameter LCDM model provides a successful fit to WMAP data, both when the data are analyzed alone and in combination with other cosmological data. Even so, it is appropriate to search for any hints of deviations from the now standard model of cosmology, which includes inflation, dark energy, dark matter, baryons, and neutrinos. The cosmological community has subjected the WMAP data to extensive and varied analyses. While there is widespread agreement as to the overall success of the six-parameter LCDM model, various "anomalies" have been reported relative to that model. In this paper we examine potential anomalies and present analyses and assessments of their significance. In most cases we find that claimed anomalies depend on posterior selection of some aspect or subset of the data. Compared with sky simulations based on the best fit model, one can select for low probability features of the WMAP data. Low probability features are expected, but it is not usually straightforward to determine whether any particular low probability feature is the result of the a posteriori selection or of non-standard cosmology. We examine in detail the properties of the power spectrum with respect to the LCDM model. We examine several potential or previously claimed anomalies in the sky maps and power spectra, including cold spots, low quadrupole power, quadropole-octupole alignment, hemispherical or dipole power asymmetry, and quadrupole power asymmetry. We conclude that there is no compelling evidence for deviations from the LCDM model, which is generally an acceptable statistical fit to WMAP and other cosmological data.

(Abridged) We present WMAP seven-year observations of bright sources which are often used as calibrators at microwave frequencies. Ten objects are studied in five frequency bands (23 - 94 GHz): the outer planets (Mars, Jupiter, Saturn, Uranus and Neptune) and five fixed celestial sources (Cas A, Tau A, Cyg A, 3C274 and 3C58). The seven-year analysis of Jupiter provides temperatures which are within 1-sigma of the previously published WMAP five-year values, with slightly tighter constraints on variability with orbital phase, and limits (but no detections) on linear polarization. Scaling factors are provided which, when multiplied by the Wright Mars thermal model predictions at 350 micron, reproduce WMAP seasonally averaged observations of Mars within ~2%. An empirical model is described which fits brightness variations of Saturn due to geometrical effects and can be used to predict the WMAP observations to within 3%. Seven-year mean temperatures for Uranus and Neptune are also tabulated. Uncertainties in Uranus temperatures are 3%-4% in the 41, 61 and 94 GHz bands; the smallest uncertainty for Neptune is ~8% for the 94 GHz band. Intriguingly, the spectrum of Uranus appears to show a dip at ~30 GHz of unidentified origin, although the feature is not of high statistical significance. Flux densities for the five selected fixed celestial sources are derived from the seven-year WMAP sky maps, and are tabulated for Stokes I, Q and U, along with polarization fraction and position angle. Fractional uncertainties for the Stokes I fluxes are typically 1% to 3%. Source variability over the seven-year baseline is also estimated. Significant secular decrease is seen for Cas A and Tau A: our results are consistent with a frequency independent decrease of about 0.53% per year for Cas A and 0.22% per year for Tau A.

The last few years have seen a surge in excitement about measurements of statistics of the primordial fluctuations beyond the power spectrum. New ideas for precision tests of Gaussianity and statistical isotropy in the data are developing simultaneously with proposals for a wide range of new theoretical possibilities. From both the observations and theory, it has become clear that there is a huge discovery potential from upcoming measurements. In this Special Issue of Advances in Astronomy we have collected articles that summarize the theoretical predictions for departures from Gaussianity or statistical isotropy from a variety of potential sources, together with the observational approaches to test these properties using the CMB or large-scale structure. We hope this collection provides an accessible entry point to these topics as they currently stand, indicating what direction future developments may take and demonstrating why these questions are so compelling. The Special Issue is available at http://www.hindawi.com/journals/aa/2010/si.gsiu.html, and individual articles are also available on the arXiv.

[Abridged] We present updated estimates of Galactic foreground emission using seven years of WMAP data. Using the power spectrum of differences between multi-frequency template-cleaned maps, we find no evidence for foreground contamination outside of the updated (KQ85y7) foreground mask. We place a 15 microKelvin upper bound on rms foreground contamination in the cleaned maps used for cosmological analysis. We find no indication in the polarization data of an extra "haze" of hard synchrotron emission from energetic electrons near the Galactic center. We provide an updated map of the cosmic microwave background (CMB) using the internal linear combination (ILC) method, updated foreground masks, and updates to point source catalogs with 62 newly detected sources. Also new are tests of the Markov chain Monte Carlo (MCMC) foreground fitting procedure against systematics in the time-stream data, and tests against the observed beam asymmetry. Within a few degrees of the Galactic plane, WMAP total intensity data show a rapidly steepening spectrum from 20-40 GHz, which may be due to emission from spinning dust grains, steepening synchrotron, or other effects. Comparisons are made to a 1-degree 408 MHz map (Haslam et al.) and the 11-degree ARCADE 2 data (Singal et al.). We find that spinning dust or steepening synchrotron models fit the combination of WMAP and 408 MHz data equally well. ARCADE data appear inconsistent with the steepening synchrotron model, and consistent with the spinning dust model, though some discrepancies remain regarding the relative strength of spinning dust emission. More high-resolution data in the 10-40 GHz range would shed much light on these issues.

In 2004, Creminelli and Zaldarriaga proposed a consistency relation for the primordial curvature perturbation of all single-field inflation models; it related the bispectrum in the squeezed limit to the spectral tilt. We have developed a technique, based in part on the Creminelli and Zaldarriaga argument, that can greatly simplify the calculation of the squeezed-limit bispectrum using the in-in formalism; we were able to arrive at a generic formula that does not rely on a slow-roll approximation. Using our formula, we explicitly tested the consistency relation for power-law inflation and for an exactly scale-invariant model by Starobinsky; for the latter model, Creminelli and Zaldarriaga's argument predicts a vanishing bispectrum whereas our quantum calculation shows a non-zero bispectrum that approaches zero in the long-wavelength limit and for inflation with a large number of e-folds.

(Abridged) The 7-year WMAP data and improved astrophysical data rigorously test the standard cosmological model and its extensions. By combining WMAP with the latest distance measurements from BAO and H0 measurement, we determine the parameters of the simplest LCDM model. The power-law index of the primordial power spectrum is n_s=0.968+-0.012, a measurement that excludes the scale-invariant spectrum by 99.5%CL. The other parameters are also improved from the 5-year results. Notable examples of improved parameters are the total mass of neutrinos, sum(m_nu)<0.58eV, and the effective number of neutrino species, N_eff=4.34+0.86-0.88. We detect the effect of primordial helium on the temperature power spectrum and provide a new test of big bang nucleosynthesis. We detect, and show on the map for the first time, the tangential and radial polarization patterns around hot and cold spots of temperature fluctuations, an important test of physical processes at z=1090 and the dominance of adiabatic scalar fluctuations. With the 7-year TB power spectrum, the limit on a rotation of the polarization plane due to potential parity-violating effects has improved to Delta(alpha)=-1.1+-1.4(stat)+-1.5(syst) degrees. We report significant detections of the SZ effect at the locations of known clusters of galaxies. The measured SZ signal agrees well with the expected signal from the X-ray data. However, it is a factor of 0.5 to 0.7 times the predictions from "universal profile" of Arnaud et al., analytical models, and hydrodynamical simulations. We find, for the first time in the SZ effect, a significant difference between the cooling-flow and non-cooling-flow clusters (or relaxed and non-relaxed clusters), which can explain some of the discrepancy. This lower amplitude is consistent with the lower-than-theoretically-expected SZ power spectrum recently measured by the South Pole Telescope collaboration.

We present a catalog of emission-line galaxies selected solely by their emission-line fluxes using a wide-field integral field spectrograph. This work is partially motivated as a pilot survey for the upcoming Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We describe the observations, reductions, detections, redshift classifications, line fluxes, and counterpart information for 397 emission-line galaxies detected over 169 sq.arcmin with a 3500-5800 Ang. bandpass under 5 Ang. full-width-half-maximum (FWHM) spectral resolution. The survey's best sensitivity for unresolved objects under photometric conditions is between 4-20 E-17 erg/s/sq.cm depending on the wavelength, and Ly-alpha luminosities between 3-6 E42 erg/s are detectable. This survey method complements narrowband and color-selection techniques in the search for high redshift galaxies with its different selection properties and large volume probed. The four survey fields within the COSMOS, GOODS-N, MUNICS, and XMM-LSS areas are rich with existing, complementary data. We find 104 galaxies via their high redshift Ly-alpha emission at 1.9<z<3.8, and the majority of the remainder objects are low redshift [OII]3727 emitters at z<0.56. The classification between low and high redshift objects depends on rest frame equivalent width, as well as other indicators, where available. Based on matches to X-ray catalogs, the active galactic nuclei (AGN) fraction amongst the Ly-alpha emitters (LAEs) is 6%. We also analyze the survey's completeness and contamination properties through simulations. We find five high-z, highly-significant, resolved objects with full-width-half-maximum sizes >44 sq.arcsec which appear to be extended Ly-alpha nebulae. We also find three high-z objects with rest frame Ly-alpha equivalent widths above the level believed to be achievable with normal star formation, EW(rest)>240 Ang.

We study the evolution of linear density fluctuations of free-streaming massive neutrinos at redshift of z<1000, with an explicit justification on the use of a fluid approximation. We solve the collisionless Boltzmann equation in an Einstein de-Sitter (EdS) universe, truncating the Boltzmann hierarchy at lmax=1 and 2, and compare the resulting density contrast of neutrinos, \delta_{\nu}^{fluid}, with that of the exact solutions of the Boltzmann equation that we derive in this paper. Roughly speaking, the fluid approximation is accurate if neutrinos were already non-relativistic when the neutrino density fluctuation of a given wavenumber entered the horizon. We find that the fluid approximation is accurate at few to 25% for massive neutrinos with 0.05<m_{\nu}<0.5eV at the scale of k<0.4~hMpc^{-1} and redshift of z<10. This result quantifies the limitation of the fluid approximation, for the massive neutrinos with m_{\nu}<0.5eV. We also find that the density contrast calculated from fluid equations (i.e., continuity and Euler equations) becomes a better approximation at a lower redshift, and the accuracy can be further improved by including an anisotropic stress term in the Euler equation. The anisotropic stress term effectively increases the pressure term by a factor of 9/5.

We calculate the effects of velocity-dependent dark matter annihilation cross sections on the intensity of the extragalactic gamma-ray background. Our formalism does not assume a locally thermal distribution of dark matter particles in phase space, and is valid for arbitrary velocity-dependent annihilation. As concrete examples, we calculate the effects of p-wave annihilation (with the $v$-weighted cross section of $\sigma v=a+bv^2$) on the mean intensity of extragalactic gamma rays produced in cosmological dark matter halos. This velocity variation makes the shape of the energy spectrum harder, but this change in the shape is too small to see unless $b/a\agt 10^6$. While we find no such models in the parameter space of the Minimal Supersymmetric Standard Model (MSSM), we show that it is possible to find $b/a\agt 10^6$ in the extension MSSM$\otimes U(1)_{B-L}$. However, we find that the most dominant effect of the p-wave annihilation is the suppression of the amplitude of the gamma-ray background. A non-zero $b$ at the dark matter freeze-out epoch requires a smaller value of $a$ in order for the relic density constraint to be satisfied, suppressing the amplitude by a factor as low as $10^{-6}$ for a thermal relic. Non-thermal relics will have weaker amplitude suppression. As another velocity-dependent effect, we calculate the spectrum for s-wave annihilation into fermions enhanced by the attractive Sommerfeld effect. Resonances associated with this effect result in significantly enhanced intensities, with a slightly softer energy spectrum.

A large sub-mm survey with Herschel will enable many exciting science opportunities, especially in an era of wide-field optical and radio surveys and high resolution cosmic microwave background experiments. The Herschel-SPIRE Legacy Survey (HSLS), will lead to imaging data over 4000 sq. degrees at 250, 350, and 500 micron. Major Goals of HSLS are: (a) produce a catalog of 2.5 to 3 million galaxies down to 26, 27 and 33 mJy (50% completeness; 5 sigma confusion noise) at 250, 350 and 500 micron, respectively, in the southern hemisphere (3000 sq. degrees) and in an equatorial strip (1000 sq. degrees), areas which have extensive multi-wavelength coverage and are easily accessible from ALMA. Two thirds of the of the sources are expected to be at z > 1, one third at z > 2 and about a 1000 at z > 5. (b) Remove point source confusion in secondary anisotropy studies with Planck and ground-based CMB data. (c) Find at least 1200 strongly lensed bright sub-mm sources leading to a 2% test of general relativity. (d) Identify 200 proto-cluster regions at z of 2 and perform an unbiased study of the environmental dependence of star formation. (e) Perform an unbiased survey for star formation and dust at high Galactic latitude and make a census of debris disks and dust around AGB stars and white dwarfs.

(Abridged) We present the angular power spectra derived from the 7-year maps and discuss the cosmological conclusions that can be inferred from WMAP data alone. The third acoustic peak in the TT spectrum is now well measured by WMAP. In the context of a flat LambdaCDM model, this improvement allows us to place tighter constraints on the matter density from WMAP data alone, and on the epoch of matter-radiation equality, The temperature-polarization (TE) spectrum is detected in the 7-year data with a significance of 20 sigma, compared to 13 sigma with the 5-year data. The low-l EE spectrum, a measure of the optical depth due to reionization, is detected at 5.5 sigma significance when averaged over l = 2-7. The BB spectrum, an important probe of gravitational waves from inflation, remains consistent with zero. The upper limit on tensor modes from polarization data alone is a factor of 2 lower with the 7-year data than it was using the 5-year data (Komatsu et al. 2010). We test the parameter recovery process for bias and find that the scalar spectral index, ns, is biased high, but only by 0.09 sigma, while the remaining parameters are biased by < 0.15 sigma. The improvement in the third peak measurement leads to tighter lower limits from WMAP on the number of relativistic degrees of freedom (e.g., neutrinos) in the early universe: Neff > 2.7 (95% CL). Also, using WMAP data alone, the primordial helium mass fraction is found to be YHe = 0.28+0.14-0.15, and with data from higher-resolution CMB experiments included, we now establish the existence of pre-stellar helium at > 3 sigma (Komatsu et al. 2010).

To quantify how rare the bullet-cluster-like high-velocity merging systems are in the standard LCDM cosmology, we use a large-volume 27 (Gpc/h)^3 MICE simulation to calculate the distribution of infall velocities of subclusters around massive main clusters. The infall-velocity distribution is given at (1-3)R_{200} of the main cluster (where R_{200} is similar to the virial radius), and thus it gives the distribution of realistic initial velocities of subclusters just before collision. These velocities can be compared with the initial velocities used by the non-cosmological hydrodynamical simulations of 1E0657-56 in the literature. The latest parameter search carried out recently by Mastropietro and Burkert showed that the initial velocity of 3000 km/s at about 2R_{200} is required to explain the observed shock velocity, X-ray brightness ratio of the main and subcluster, and displacement of the X-ray peaks from the mass peaks. We show that such a high infall velocity at 2R_{200} is incompatible with the prediction of a LCDM model: the probability of finding 3000 km/s in (2-3)R_{200} is between 3.3X10^{-11} and 3.6X10^{-9}. It is concluded that the existence of 1E0657-56 is incompatible with the prediction of a LCDM model, unless a lower infall velocity solution for 1E0657-56 with < 1800 km/s at 2R_{200} is found.

Since the first limit on the (local) primordial non-Gaussianity parameter, fNL, was obtained from COBE data in 2002, observations of the CMB have been playing a central role in constraining the amplitudes of various forms of non-Gaussianity in primordial fluctuations. The current 68% limit from the 7-year WMAP data is fNL=32+/-21, and the Planck satellite is expected to reduce the uncertainty by a factor of four in a few years from now. If fNL>>1 is found by Planck with high statistical significance, all single-field models of inflation would be ruled out. Moreover, if the Planck satellite finds fNL=30, then it would be able to test a broad class of multi-field models using the four-point function (trispectrum) test of tauNL>=(6fNL/5)^2. In this article, we review the methods (optimal estimator), results (WMAP 7-year), and challenges (secondary anisotropy, second-order effect, and foreground) of measuring primordial non-Gaussianity from the CMB data, present a science case for the trispectrum, and conclude with future prospects.

(Abridged) New full sky temperature and polarization maps based on seven years of data from WMAP are presented. The new results are consistent with previous results, but have improved due to reduced noise from the additional integration time, improved knowledge of the instrument performance, and improved data analysis procedures. The improvements are described in detail. The seven year data set is well fit by a minimal six-parameter flat Lambda-CDM model. The parameters for this model, using the WMAP data in conjunction with baryon acoustic oscillation data from the Sloan Digital Sky Survey and priors on H_0 from Hubble Space Telescope observations, are: Omega_bh^2 = 0.02260 +-0.00053, Omega_ch^2 = 0.1123 +-0.0035, Omega_Lambda = 0.728 +0.015 -0.016, n_s = 0.963 +-0.012, tau = 0.087 +-0.014 and sigma_8 = 0.809 +-0.024 (68 % CL uncertainties). The temperature power spectrum signal-to-noise ratio per multipole is greater that unity for multipoles < 919, allowing a robust measurement of the third acoustic peak. This measurement results in improved constraints on the matter density, Omega_mh^2 = 0.1334 +0.0056 -0.0055, and the epoch of matter- radiation equality, z_eq = 3196 +134 -133, using WMAP data alone. The new WMAP data, when combined with smaller angular scale microwave background anisotropy data, results in a 3 sigma detection of the abundance of primordial Helium, Y_He = 0.326 +-0.075.The power-law index of the primordial power spectrum is now determined to be n_s = 0.963 +-0.012, excluding the Harrison-Zel'dovich-Peebles spectrum by >3 sigma. These new WMAP measurements provide important tests of Big Bang cosmology.

The Near Infrared Background (NIRB) is one of a few methods that can be used to observe the redshifted light from early stars at a redshift of six and above. Fluctuations of the NIRB can provide information on the first structures, such as halos and their surrounding ionized regions in the IGM. We combine, for the first time, N-body simulations, radiative transfer code, and analytic calculations of luminosity of early structures to predict the angular power spectrum (C_l) of fluctuations in the NIRB. We study the effects of various assumptions about the stellar mass, the initial mass spectrum of stars, metallicity, the star formation efficiency (f_*), the escape fraction of ionizing photons (f_esc), and the star formation timescale (t_SF), on the amplitude as well as the shape of C_l. The power spectrum of NIRB fluctuations is maximized when f_* is the largest (as C_l ~ (f_*)^2) and f_esc is the smallest. A significant uncertainty in the predicted amplitude of C_l exists due to our lack of knowledge of t_SF of these galaxies, which is equivalent to our lack of knowledge of the mass-to-light ratio. We do not see a turnover in the NIRB angular power spectrum of the halo contribution and explain this as the effect of high levels of non-linear bias. This is partly due to our choice of the minimum mass of halos contributing to NIRB, and a smaller minimum mass, which has a smaller non-linear bias, may still exhibit a turn over. Therefore, both the amplitude and shape of the NIRB power spectrum provide important information regarding the nature of sources contributing to the cosmic reionization. The angular power spectrum of the IGM, in most cases, is much smaller than the halo angular power spectrum. In addition, low levels of the observed mean background intensity tend to rule out high values of f_* > 0.2.

We describe a sampling method to estimate the polarized CMB signal from observed maps of the sky. We use a Metropolis-within-Gibbs algorithm to estimate the polarized CMB map, containing Q and U Stokes parameters at each pixel, and its covariance matrix. These can be used as inputs for cosmological analyses. The polarized sky signal is parameterized as the sum of three components: CMB, synchrotron emission, and thermal dust emission. The polarized Galactic components are modeled with spatially varying power law spectral indices for the synchrotron, and a fixed power law for the dust, and their component maps are estimated as by-products. We apply the method to simulated low resolution maps with pixels of side 7.2 degrees, using diagonal and full noise realizations drawn from the WMAP noise matrices. The CMB maps are recovered with goodness of fit consistent with errors. Computing the likelihood of the E-mode power in the maps as a function of optical depth to reionization, tau, for fixed temperature anisotropy power, we recover tau=0.091+-0.019 for a simulation with input tau=0.1, and mean tau=0.098 averaged over 10 simulations. A `null' simulation with no polarized CMB signal has maximum likelihood consistent with tau=0. The method is applied to the five-year WMAP data, using the K, Ka, Q and V channels. We find tau=0.090+-0.019, compared to tau=0.086+-0.016 from the template-cleaned maps used in the primary WMAP analysis. The synchrotron spectral index, beta, averaged over high signal-to-noise pixels with standard deviation sigma(beta)<0.25, but excluding ~6% of the sky masked in the Galactic plane, is -3.03+-0.04. This estimate does not vary significantly with Galactic latitude, although includes an informative prior.

A convincing detection of primordial non-Gaussianity in the local form of the bispectrum, whose amplitude is given by the fNL parameter, offers a powerful test of inflation. In this paper we calculate the modification of two-point cross-correlation statistics of weak lensing - galaxy-galaxy lensing and galaxy-Cosmic Microwave Background (CMB) cross-correlation - due to fNL. We derive and calculate the covariance matrix of galaxy-galaxy lensing including cosmic variance terms. We focus on large scales (l<100) for which the shape noise of the shear measurement becomes irrelevant and cosmic variance dominates the error budget. For a modest degree of non-Gaussianity, fNL=+/-50, modifications of the galaxy-galaxy lensing signal at the 10% level are seen on scales R~300 Mpc, and grow rapidly toward larger scales as \propto R^2. We also see a clear signature of the baryonic acoustic oscillation feature in the matter power spectrum at ~150 Mpc, which can be measured by next-generation lensing experiments. In addition we can probe the local-form primordial non-Gaussianity in the galaxy-CMB lensing signal by correlating the lensing potential reconstructed from CMB with high-z galaxies. For example, for fNL=+/-50, we find that the galaxy-CMB lensing cross power spectrum is modified by ~10% at l~40, and by a factor of two at l~10, for a population of galaxies at z=2 with a bias of 2. The effect is greater for more highly biased populations at larger z; thus, high-z galaxy surveys cross-correlated with CMB offer a yet another probe of primordial non-Gaussianity.

We use Minkowski Functionals (MF) to constrain a primordial non-Gaussian contribution to the CMB intensity field as observed in the 150 GHz and 145 GHz BOOMERanG maps from the 1998 and 2003 flights, respectively, performing for the first time a joint analysis of the two datasets. A perturbative expansion of the MF formulae in the limit of a weakly non-Gaussian field yields analytical formulae, derived by Hikage et al. (2006), which can be used to constrain the coupling parameter f_NL without the need for non-Gaussian simulations. We find -1020<f_NL<390 at 95% CL, significantly improving the previous constraints by De Troia et al. (2007) on the BOOMERanG 2003 dataset. These are the best f_NL limits to date for suborbital probes.

We calculate the bispectrum, B_g(k_1,k_2,k_3), Fourier transform of the three-point function of density peaks (e.g., galaxies), using two different methods: the Matarrese-Lucchin-Bonometto formula and the locality of galaxy bias. The bispectrum of peaks is not only sensitive to that of the underlying matter density fluctuations, but also to the four-point function. For a physically-motivated, local form of primordial non-Gaussianity in the curvature perturbation, we show that the galaxy bispectrum contains five physically distinct pieces: (i) non-linear gravitational evolution, (ii) non-linear galaxy bias, (iii) f_nl, (iv) f_nl^2, and (v) \gnl. While (i), (ii), and a part of (iii) have been derived in the literature, (iv) and (v) are derived in this paper for the first time. Our finding suggests that the galaxy bispectrum is more sensitive to f_nl than previously recognized, and is also sensitive to a new term, g_nl. For a more general form of local-type non-Gaussianity, the coefficient \fnl^2 can be interpreted as \tau_nl, which allows us to test multi-field inflation models. The usual terms from Gaussian initial conditions, have the smallest signals in the squeezed configurations, while the others have the largest signals; thus, we can distinguish them easily. We cannot interpret the effects of f_nl on B_g(k_1,k_2,k_3) as a scale-dependent bias, and thus replacing the linear bias in the galaxy bispectrum with the scale-dependent bias known for the power spectrum results in an incorrect prediction. As the importance of primordial non-Gaussianity relative to the non-linear gravity evolution and galaxy bias increases toward higher redshifts, galaxy surveys probing a high-redshift universe are particularly useful for probing the primordial non-Gaussianity.

We calculate the bispectrum of the Cosmic Microwave Background (CMB) temperature anisotropies induced by the second-order fluctuations in the Boltzmann equation. In this paper, which is one of a series of papers on the numerical calculation of the bispectrum from the second-order fluctuations, we consider the terms that are products of the first-order perturbations, and leave intrinsically second-order terms and perturbations in the recombination history to the subsequent papers. We show that the bispectrum has the maximum signal in the squeezed triangles, similar to the local-type primordial bispectrum, as both types generate non-linearities via products of the first-order terms in position space. However, detailed calculations show that their shapes are sufficiently different: the cross-correlation coefficient reaches 0.5 at the maximum multipole of l_{max}~ 200, and then weakens to 0.3 at l_{max}~ 2000. The differences in shape arise from (i) the way the acoustic oscillations affect the bispectrum, and (ii) the second-order effects not being scale-invariant. This implies that the contamination of the primordial bispectrum due to the second-order effects (from the products of the first-order terms) is small. The expected signal-to-noise ratio of the products of the first-order terms is ~ 0.4 at l_{max}~ 2000 for a full-sky, cosmic variance limited experiment. We therefore conclude that the products of the first-order terms may be safely ignored in the analysis of the future CMB experiments. The expected contamination of the local-form f_{NL} is f^{local}_{NL}~ 0.9 at l_{max}~ 200, and f^{local}_{NL}~ 0.5 at l_{max}~ 2000.

A new and powerful probe of the origin and evolution of structures in the Universe has emerged and been actively developed over the last decade. In the coming decade, non-Gaussianity, i.e., the study of non-Gaussian contributions to the correlations of cosmological fluctuations, will become an important probe of both the early and the late Universe. Specifically, it will play a leading role in furthering our understanding of two fundamental aspects of cosmology and astrophysics: (i) the physics of the very early universe that created the primordial seeds for large-scale structures, and (ii) the subsequent growth of structures via gravitational instability and gas physics at later times. To date, observations of fluctuations in the Cosmic Microwave Background (CMB) and the Large-Scale Structure of the Universe (LSS) have focused largely on the Gaussian contribution as measured by the two-point correlations (or the power spectrum) of density fluctuations. However, an even greater amount of information is contained in non-Gaussianity and a large discovery space therefore still remains to be explored. Many observational probes can be used to measure non-Gaussianity, including CMB, LSS, gravitational lensing, Lyman-alpha forest, 21-cm fluctuations, and the abundance of rare objects such as clusters of galaxies and high-redshift galaxies. Not only does the study of non-Gaussianity maximize the science return from a plethora of present and future cosmological experiments and observations, but it also carries great potential for important discoveries in the coming decade.

We calculate the non-linear matter power spectrum using the 3rd-order perturbation theory without ignoring the pressure gradient term. We consider a semi-realistic system consisting of two matter components with and without pressure, and both are expanded into the 3rd order in perturbations in a self-consistent manner, for the first time. While the pressured component may be identified with baryons or neutrinos, in this paper we mainly explore the physics of the non-linear pressure effect using a toy model in which the Jeans length does not depend on time, i.e., the sound speed decreases as 1/a^{1/2}, where a is the scale factor. The linear analysis shows that the power spectrum below the so-called filtering scale is suppressed relative to the power spectrum of the cold dark matter. Our non-linear calculation shows that the actual filtering scale for a given sound speed is smaller than the linear filtering scale by a factor depending on the redshift and the Jeans length. A ~40% change is common, and our results suggest that, when applied to baryons, the temperature of the Inter-galactic Medium inferred from the filtering scale observed in the flux power spectrum of Lyman-alpha forests would be underestimated by a factor of two, if one used the linear filtering scale to interpret the data. The filtering mass, which is proportional to the filtering scale cubed, can also be significantly smaller than the linear theory prediction especially at low redshift, where the actual filtering mass can be smaller than the linear prediction by a factor of three. Finally, when applied to neutrinos, we find that neutrino perturbations deviate significantly from linear perturbations even below the free-streaming scales, and thus neutrinos cannot be treated as linear perturbations.

We summarize the utility of precise cosmic microwave background (CMB) polarization measurements as probes of the physics of inflation. We focus on the prospects for using CMB measurements to differentiate various inflationary mechanisms. In particular, a detection of primordial B-mode polarization would demonstrate that inflation occurred at a very high energy scale, and that the inflaton traversed a super-Planckian distance in field space. We explain how such a detection or constraint would illuminate aspects of physics at the Planck scale. Moreover, CMB measurements can constrain the scale-dependence and non-Gaussianity of the primordial fluctuations and limit the possibility of a significant isocurvature contribution. Each such limit provides crucial information on the underlying inflationary dynamics. Finally, we quantify these considerations by presenting forecasts for the sensitivities of a future satellite experiment to the inflationary parameters.

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

(Brief Summary) What is the total radiative content of the Universe since the epoch of recombination? The extragalactic background light (EBL) spectrum captures the redshifted energy released from the first stellar objects, protogalaxies, and galaxies throughout cosmic history. Yet, we have not determined the brightness of the extragalactic sky from UV/optical to far-infrared wavelengths with sufficient accuracy to establish the radiative content of the Universe to better than an order of magnitude. Among many science topics, an accurate measurement of the EBL spectrum from optical to far-IR wavelengths, will address: What is the total energy released by stellar nucleosynthesis over cosmic history? Was significant energy released by non-stellar processes? Is there a diffuse component to the EBL anywhere from optical to sub-millimeter? When did first stars appear and how luminous was the reionization epoch? Absolute optical to mid-IR EBL spectrum to an astrophysically interesting accuracy can be established by wide field imagingat a distance of 5 AU or above the ecliptic plane where the zodiacal foreground is reduced by more than two orders of magnitude.

The hidden scalar field, which couples to the visible sector only through Planck-suppressed interactions, is a candidate for dark matter owing to its long lifetime. Decay of such a scalar field offers observational tests of this scenario. We show that decay of the hidden scalar field can explain the observed excess of high-energy positrons/electrons observed by ATIC/PPB-BETS, for a suitable choice of the mass and the vacuum expectation value of the field. We also show that the same choice of the parameters gives the observed dark matter abundance. Such a remarkable coincidence suggests that the Planck-suppressed interactions may be responsible for the observed excess in the cosmic-ray positrons/electrons.

We summarize existing constraints on the epoch of reionization and discuss the observational probes that are sensitive to the process. We focus on the role large scale polarization can play. Polarization probes the integrated optical depth across the entire epoch of reionization. Future missions such as Planck and CMBPol will greatly enhance our knowledge of the reionization history, allowing us to measure the time evolution of the ionization fraction. As large scale polarization probes high redshift activity, it can best constrain models where the Universe was fully or partially ionized at early times. In fact, large scale polarization could be our only probe of the highest redshifts.

We present a new estimate of foreground emission in the WMAP data, using a Markov chain Monte Carlo (MCMC) method. The new technique delivers maps of each foreground component for a variety of foreground models, error estimates of the uncertainty of each foreground component, and provides an overall goodness-of-fit measurement. The resulting foreground maps are in broad agreement with those from previous techniques used both within the collaboration and by other authors. We find that for WMAP data, a simple model with power-law synchrotron, free-free, and thermal dust components fits 90% of the sky with a reduced chi-squared of 1.14. However, the model does not work well inside the Galactic plane. The addition of either synchrotron steepening or a modified spinning dust model improves the fit. This component may account for up to 14% of the total flux at Ka-band (33 GHz). We find no evidence for foreground contamination of the CMB temperature map in the 85% of the sky used for cosmological analysis.

This paper focuses on cosmological constraints derived from analysis of WMAP data alone. A simple LCDM cosmological model fits the five-year WMAP temperature and polarization data. The basic parameters of the model are consistent with the three-year data and now better constrained: Omega_b h^2 = 0.02273+-0.00062, Omega_c h^2 = 0.1099+-0.0062, Omega_L = 0.742+-0.030, n_s = 0.963+0.014- 0.015, tau = 0.087+-0.017, sigma_8 = 0.796+-0.036. With five years of polarization data, we have measured the optical depth to reionization, tau>0, at 5 sigma significance. The redshift of an instantaneous reionization is constrained to be z_reion = 11.0+-1.4 with 68% confidence. This excludes a sudden reionization of the universe at z=6 at more than 3.5 sigma significance, suggesting that reionization was an extended process. Using two methods for polarized foreground cleaning we get consistent estimates for the optical depth, indicating an error due to foreground treatment of tau~0.01. This cosmological model also fits small-scale CMB data, and a range of astronomical data measuring the expansion rate and clustering of matter in the universe. We find evidence for the first time in the CMB power spectrum for a non-zero cosmic neutrino background, or a background of relativistic species, with the standard three light neutrino species preferred over the best-fit LCDM model with N_eff=0 at >99.5% confidence, and N_eff > 2.3 (95% CL) when varied. The five-year WMAP data improve the upper limit on the tensor-to-scalar ratio to r < 0.43 (95% CL), for power-law models. With longer integration we find no evidence for a running spectral index, with dn_s/dlnk = -0.037+-0.028.

We present new full-sky temperature and polarization maps in five frequency bands from 23 to 94 GHz, based on data from the first five years of the WMAP sky survey. The five-year maps incorporate several improvements in data processing made possible by the additional years of data and by a more complete analysis of the instrument calibration and in-flight beam response. We present several new tests for systematic errors in the polarization data and conclude that Ka band data (33 GHz) is suitable for use in cosmological analysis, after foreground cleaning. This significantly reduces the overall polarization uncertainty. With the 5 year WMAP data, we detect no convincing deviations from the minimal 6-parameter LCDM model: a flat universe dominated by a cosmological constant, with adiabatic and nearly scale-invariant Gaussian fluctuations. Using WMAP data combined with measurements of Type Ia supernovae and Baryon Acoustic Oscillations, we find (68% CL uncertainties): Omega_bh^2 = 0.02267 \pm 0.00059, Omega_ch^2 = 0.1131 \pm 0.0034, Omega_Lambda = 0.726 \pm 0.015, n_s = 0.960 \pm 0.013, tau = 0.084 \pm 0.016, and Delta_R^2 = (2.445 \pm 0.096) x 10^-9. From these we derive: sigma_8 = 0.812 \pm 0.026, H_0 = 70.5 \pm 1.3 km/s/Mpc, z_{reion} = 10.9 \pm 1.4, and t_0 = 13.72 \pm 0.12 Gyr. The new limit on the tensor-to-scalar ratio is r < 0.22 (95% CL). We obtain tight, simultaneous limits on the (constant) dark energy equation of state and spatial curvature: -0.14 < 1+w < 0.12 and -0.0179 < Omega_k < 0.0081 (both 95% CL). The number of relativistic degrees of freedom (e.g. neutrinos) is found to be N_{eff} = 4.4 \pm 1.5, consistent with the standard value of 3.04. Models with N_{eff} = 0 are disfavored at >99.5% confidence.

We present the temperature and polarization angular power spectra of the cosmic microwave background (CMB) derived from the first 5 years of WMAP data. The 5-year temperature (TT) spectrum is cosmic variance limited up to multipole l=530, and individual l-modes have S/N>1 for l<920. The best fitting six-parameter LambdaCDM model has a reduced chi^2 for l=33-1000 of chi^2/nu=1.06, with a probability to exceed of 9.3%. There is now significantly improved data near the third peak which leads to improved cosmological constraints. The temperature-polarization correlation (TE) is seen with high significance. After accounting for foreground emission, the low-l reionization feature in the EE power spectrum is preferred by \Delta\chi^2=19.6 for optical depth tau=0.089 by the EE data alone, and is now largely cosmic variance limited for l=2-6. There is no evidence for cosmic signal in the BB, TB, or EB spectra after accounting for foreground emission. We find that, when averaged over l=2-6, l(l+1)C^{BB}_l/2\pi < 0.15 uK^2 (95% CL).

We present the list of point sources found in the WMAP 5-year maps. The technique used in the first-year and three-year analysis now finds 390 point sources, and the five-year source catalog is complete for regions of the sky away from the galactic plane to a 2 Jy limit, with SNR > 4.7 in all bands in the least covered parts of the sky. The noise at high frequencies is still mainly radiometer noise, but at low frequencies the CMB anisotropy is the largest uncertainty. A separate search of CMB-free V-W maps finds 99 sources of which all but one can be identified with known radio sources. The sources seen by WMAP are not strongly polarized. Many of the WMAP sources show significant variability from year to year, with more than a 2:1 range between the minimum and maximum fluxes.

(Abridged) The WMAP 5-year data strongly limit deviations from the minimal LCDM model. We constrain the physics of inflation via Gaussianity, adiabaticity, the power spectrum shape, gravitational waves, and spatial curvature. We also constrain the properties of dark energy, parity-violation, and neutrinos. We detect no convincing deviations from the minimal model. The parameters of the LCDM model, derived from WMAP combined with the distance measurements from the Type Ia supernovae (SN) and the Baryon Acoustic Oscillations (BAO), are: Omega_b=0.0456+-0.0015, Omega_c=0.228+-0.013, Omega_Lambda=0.726+-0.015, H_0=70.5+-1.3 km/s/Mpc, n_s=0.960+-0.013, tau=0.084+-0.016, and sigma_8=0.812+-0.026. With WMAP+BAO+SN, we find the tensor-to-scalar ratio r<0.22 (95% CL), and n_s>1 is disfavored regardless of r. We obtain tight, simultaneous limits on the (constant) equation of state of dark energy and curvature. We provide a set of "WMAP distance priors," to test a variety of dark energy models. We test a time-dependent w with a present value constrained as -0.33<1+w_0<0.21 (95% CL). Temperature and matter fluctuations obey the adiabatic relation to within 8.9% and 2.1% for the axion and curvaton-type dark matter, respectively. The TE and EB spectra constrain cosmic parity-violation. We find the limit on the total mass of neutrinos, sum(m_nu)<0.67 eV (95% CL), which is free from the uncertainty in the normalization of the large-scale structure data. The effective number of neutrino species is constrained as N_{eff} = 4.4+-1.5 (68%), consistent with the standard value of 3.04. Finally, limits on primordial non-Gaussianity are -9<f_{NL}^{local}<111 and -151<f_{NL}^{equil}<253 (95% CL) for the local and equilateral models, respectively.

We present a method for extracting the angular diameter distances, $D_A$, and the expansion rates, $H$, of the universe from the {\it two-dimensional} Baryon Acoustic Oscillations (BAO) in the galaxy power spectrum. Our method builds upon the existing algorithm called the "fit-and-extract" (FITEX) method, which allows one to extract only $D_A^2/H$ from a spherically averaged one-dimensional power spectrum. We develop the FITEX-2d method, an extension of the FITEX method, to include the two-dimensional information, which allows us to extract $D_A$ and $H$ simultaneously. We test the FITEX-2d method using the Millennium Simulation as well as simplified Monte Carlo simulations with a bigger volume. The BAOs, however, contain only a limited amount of information. We show that the full modeling, including the overall shape of the power spectrum, yields much better determinations of $D_A$ and $H$, hence the dark energy equation of state parameters such as $w_0$ and $w_a$, than the BAO-only analysis by more than a factor of two, provided that non-linear effects are under control.

We calculate the non-linear galaxy power spectrum in real space, including non-linear distortion of the Baryon Acoustic Oscillations, using the standard 3rd-order perturbation theory (PT). The calculation is based upon the assumption that the number density of galaxies is a local function of the underlying, non-linear density field. The galaxy bias is allowed to be both non-linear and stochastic. We show that the PT calculation agrees with the galaxy power spectrum estimated from the Millennium Simulation, in the weakly non-linear regime (defined by the matter power spectrum) at high redshifts, $1\le z\le6$. We also show that, once 3 free parameters characterizing galaxy bias are marginalized over, the PT power spectrum fit to the Millennium Simulation data yields unbiased estimates of the distance scale, $D$, to within the statistical error. This distance scale corresponds to the angular diameter distance, $D_A(z)$, and the expansion rate, $H(z)$, in real galaxy surveys. Our results presented in this paper are still restricted to real space. The future work should include the effects of non-linear redshift space distortion. Nevertheless, our results indicate that non-linear galaxy bias in the weakly non-linear regime at high redshifts is reasonably under control.

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

We develop a new interpolation scheme, based on harmonic inpainting, for reconstructing the cosmic microwave background temperature data within the Galaxy mask from the data outside the mask. We find that, for scale-invariant isotropic random Gaussian fluctuations, the developed algorithm reduces the errors in the reconstructed map for the odd-parity modes significantly for azimuthally symmetric masks with constant galactic latitudes. For a more realistic Galaxy mask, we find a modest improvement in the even-parity modes as well.

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) will outfit the 10 m HET with a new wide field and an array of 150 integral-field spectrographs to survey a 420 sq. deg. area in the north Galactic cap. Each fiber-coupled unit spectrograph will cover 350-550 nm, simultaneously. This instrument, called VIRUS, will produce ~34,000 spectra per exposure, and will open up the emission-line universe to large surveys for the first time. The survey will detect 0.8 million Lyman-alpha emitting (LAE) galaxies with 1.9<z<3.5 and more than a million [OII] emitting galaxies with z<0.5. The 3-D map of LAE galaxies in 9 cubic Gpc volume will be used to measure the expansion history at this early epoch using baryonic acoustic oscillations and the shape of the power spectrum. The aim of HETDEX is to provide a direct detection of dark energy at z~3. The measurement will constrain the evolution of dark energy and will also provide 0.1%-level accuracy on the curvature of the Universe, ten times better than current. The prototype of the VIRUS unit spectrograph (VIRUS-P) is a powerful instrument in its own right. Used on the McDonald 2.7 m, it covers the largest area of any integral field spectrograph, and reaches wavelengths down to 340 nm. VIRUS-P is being used for a pilot survey to better measure the properties of LAE galaxies in support of HETDEX. We report initial results from this survey.

We show that reheating of the universe occurs spontaneously in a broad class of inflation models with f(phi)R gravity (phi is inflaton). The model does not require explicit couplings between phi and bosonic or fermionic matter fields. The couplings arise spontaneously when phi settles in the vacuum expectation value (vev) and oscillates, with coupling constants given by derivatives of f(phi) at the vev and the mass of resulting bosonic or fermionic fields. This mechanism allows inflaton quanta to decay into any fields which are not conformally invariant in f(phi)R gravity theories.

We study implications of the large-N species solution to the hierarchy problem, proposed by G. Dvali, for reheating of the universe after inflation. Dvali's proposal contains additional N~10^{32} Z_2-conserved quantum fields beyond the Standard Model particles with mass ~1 TeV, which weaken gravity by a factor of 1/N, and thus explain the hierarchy between the Plank scale and the electroweak scale. We show that, in this scenario, the decay rates of inflaton fields through gravitational decay channels are enhanced by a factor of N, and thus they decay into N species of the quantum fields very efficiently, in the limit that quantum gravity effects are unimportant for the gravitational decay rate. In order not to over-reheat the universe, inflaton mass, vacuum expectation value of inflaton, or non-minimal gravitational coupling should be tightly fine-tuned. Our conclusion holds even when the gravitational decay is prohibited by some symmetry of the theory; the universe may still be over-reheated via annihilation of inflatons, if the number density of inflaton quanta is greater than the critical value.

In our recent paper (Yadav et al. 2007) we described a fast cubic (bispectrum) estimator of the amplitude of primordial non-Gaussianity of local type, f_{NL}, from a combined analysis of the Cosmic Microwave Background (CMB) temperature and E-polarization observations. In this paper we generalize the estimator to deal with a partial sky coverage as well as inhomogeneous noise. Our generalized estimator is still computationally efficient, scaling as O(N^3/2) compared to the O(N^5/2) scaling of the brute force bispectrum calculation for sky maps with N pixels. Upcoming CMB experiments are expected to yield high-sensitivity temperature and E-polarization data. Our generalized estimator will allow us to optimally utilize the combined CMB temperature and E-polarization information from these realistic experiments, and to constrain primordial non-Gaussianity.

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

Using a simple method to interpret the luminosity function of Ly-alpha emitters, we explore properties of Ly-alpha emitters from 5.7 < z < 8.8 with various assumptions about metallicity and stellar mass spectra. We constrain a mass-to-'observed' light ratio, M_h/L_band. For narrow-band surveys, L_band is simply related to the intrinsic Ly-alpha luminosity with a survival fraction of Ly-alpha photons, alpha_esc. The mass-to-'bolometric light', M_h/L_bol, can also be deduced, once the metallicity and stellar mass spectrum are given. The inferred M_h/L_bol is more sensitive to metallicity than to the mass spectrum. We find the following constraints on a mass-to-light ratio of Ly-alpha emitters from 5.7 < z < 7: (M_h/L_bol)(alpha_{esc}epsilon^{1/gamma})^{-1}=21-38, 14-26, and 9-17 for Z=0, 1/50, and 1 Z_sun, respectively, where epsilon is the 'duty cycle' of Ly-alpha emitters, and gamma ~ 2 is a local slope of the cumulative luminosity function. Only weak lower limits are obtained for z=8.8. Therefore, Ly-alpha emitters are consistent with either starburst galaxies M_h/L_bol ~ 0.1-1 with a smaller Ly-alpha survival fraction, alpha_{esc}epsilon^{1/gamma} ~0.01-0.05, or normal populations (M_h/L_bol ~ 10) if a good fraction of Ly-alpha photons survived, alpha_{esc}epsilon^{1/gamma} ~ 0.5-1. We find no evidence for the end of reionization in the luminosity functions of Ly-alpha emitters discovered in the current Ly-alpha surveys, including recent discovery of one Ly-alpha emitter at z=7. The data are consistent with no evolution of intrinsic properties of Ly-alpha emitters or neutral fraction in the intergalactic medium up to z=7. No detection of sources at z=8.8 does not yield a significant constraint yet. We also show that the lack of detection at z=8.8 does not rule out the high-z galaxies being the origin of the excess NIRB.

The forthcoming Planck experiment will provide high sensitivity polarization measurements that will allow us to further tighten the f_NL bounds from the temperature data. Monte Carlo simulations of non-Gaussian CMB maps have been used as a fundamental tool to characterize non-Gaussian signatures in the data, as they allow us to calibrate any statistical estimators and understand the effect of systematics, foregrounds and other contaminants. We describe an algorithm to generate high-angular resolution simulations of non-Gaussian CMB maps in temperature and polarization. We consider non-Gaussianities of the local type, for which the level of non-Gaussianity is defined by the dimensionless parameter, f_NL. We then apply the temperature and polarization fast cubic statistics recently developed by Yadav et al. to a set of non-Gaussian temperature and polarization simulations. We compare our results to theoretical expectations based on a Fisher matrix analysis, test the unbiasedness of the estimator, and study the dependence of the error bars on f_NL. All our results are in very good agreement with theoretical predictions, thus confirming the reliability of both the simulation algorithm and the fast cubic temperature and polarization estimator.

The greatest challenge in the interpretation of galaxy clustering data from any surveys is galaxy bias. Using a simple Fisher matrix analysis, we show that the bispectrum provides an excellent determination of linear and non-linear bias parameters of intermediate and high-z galaxies, when all measurable triangle configurations down to mildly non-linear scales, where perturbation theory is still valid, are included. The bispectrum is also a powerful probe of primordial non-Gaussianity. The planned galaxy surveys at z>2 should yield constraints on non-Gaussian parameters, f_{NL}^{loc.} and f_{NL}^{eq.}, that are comparable to, or even better than, those from CMB experiments. We study how these constraints improve with volume, redshift range, as well as the number density of galaxies. Finally, we show that a halo occupation distribution may be used to improve these constraints further by lifting degeneracies between gravity, bias, and primordial non-Gaussianity.

We present a full-sky model of polarized Galactic microwave emission based on three years of observations by the Wilkinson Microwave Anisotropy Probe (WMAP) at frequencies from 23 to 94 GHz. The model compares maps of the Stokes Q and U components from each of the 5 WMAP frequency bands in order to separate synchrotron from dust emission, taking into account the spatial and frequency dependence of the synchrotron and dust components. This simple two-component model of the interstellar medium accounts for at least 97% of the polarized emission in the WMAP maps of the microwave sky. Synchrotron emission dominates the polarized foregrounds at frequencies below 50 GHz, and is comparable to the dust contribution at 65 GHz. The spectral index of the synchrotron component, derived solely from polarization data, is -3.2 averaged over the full sky, with a modestly flatter index on the Galactic plane. The synchrotron emission has mean polarization fraction 2--4% in the Galactic plane and rising to over 20% at high latitude, with prominent features such as the North Galactic Spur more polarized than the diffuse component. Thermal dust emission has polarization fraction 1% near the Galactic center, rising to 6% at the anti-center. Diffuse emission from high-latitude dust is also polarized with mean fractional polarization 0.036 +/- 0.011.

The origin of the cosmic gamma-ray background (CGB) is a longstanding mystery in high-energy astrophysics. Possible candidates include ordinary astrophysical objects such as unresolved blazars, as well as more exotic processes such as dark matter annihilation. While it would be difficult to distinguish them from the mean intensity data alone, one can use anisotropy data instead. We investigate the CGB anisotropy both from unresolved blazars and dark matter annihilation (including contributions from dark matter substructures), and find that the angular power spectra from these sources are very different. We then focus on detectability of dark matter annihilation signals using the anisotropy data, by treating the unresolved blazar component as a known background. We find that the dark matter signature should be detectable in the angular power spectrum of the CGB from two-year all-sky observations with the Gamma Ray Large Area Space Telescope (GLAST), as long as the dark matter annihilation contributes to a reasonable fraction, e.g., >~ 0.3, of the CGB at around 10 GeV. We conclude that the anisotropy measurement of the CGB with GLAST should be a powerful tool for revealing the CGB origin, and potentially for the first detection of dark matter annihilation.

We present new full-sky temperature maps in five frequency bands from 23 to 94 GHz, based on the first three years of the WMAP sky survey. The new maps, which are consistent with the first-year maps and more sensitive, incorporate improvements in data processing made possible by the additional years of data and by a more complete analysis of the polarization signal. These include refinements in the gain calibration and beam response models. We employ two forms of multi-frequency analysis to separate astrophysical foreground signals from the CMB, each of which improves on our first-year analyses. First, we form an improved 'Internal Linear Combination' map, based solely on WMAP data, by adding a bias correction step and by quantifying residual uncertainties in the resulting map. Second, we fit and subtract new spatial templates that trace Galactic emission; in particular, we now use low-frequency WMAP data to trace synchrotron emission. The WMAP point source catalog is updated to include 115 new sources. We derive the angular power spectrum of the temperature anisotropy using a hybrid approach that combines a maximum likelihood estimate at low l (large angular scales) with a quadratic cross-power estimate for l>30. Our best estimate of the CMB power spectrum is derived by averaging cross-power spectra from 153 statistically independent channel pairs. The combined spectrum is cosmic variance limited to l=400, and the signal-to-noise ratio per l-mode exceeds unity up to l=850. The first two acoustic peaks are seen at l=220.8 +- 0.7 and l=530.9 +- 3.8, respectively, while the first two troughs are seen at l=412.4 +- 1.9 and l=675.1 +- 11.1, respectively. The rise to the third peak is unambiguous; when the WMAP data are combined with higher resolution CMB measurements, the existence of a third acoustic peak is well established.

A simple cosmological model with only six parameters (matter density, Omega_m h^2, baryon density, Omega_b h^2, Hubble Constant, H_0, amplitude of fluctuations, sigma_8, optical depth, tau, and a slope for the scalar perturbation spectrum, n_s) fits not only the three year WMAP temperature and polarization data, but also small scale CMB data, light element abundances, large-scale structure observations, and the supernova luminosity/distance relationship. Using WMAP data only, the best fit values for cosmological parameters for the power-law flat LCDM model are (Omega_m h^2, Omega_b h^2, h, n_s, tau, sigma_8) = 0.1277+0.0080-0.0079, 0.02229+-0.00073, 0.732+0.031-0.032, 0.958+-0.016, 0.089+-0.030, 0.761+0.049-0.048). The three year data dramatically shrink the allowed volume in this six dimensional parameter space. Assuming that the primordial fluctuations are adiabatic with a power law spectrum, the WMAP data_alone_ require dark matter, and favor a spectral index that is significantly less than the Harrison-Zel'dovich-Peebles scale-invariant spectrum (n_s=1, r=0). Models that suppress large-scale power through a running spectral index or a large-scale cut-off in the power spectrum are a better fit to the WMAP and small scale CMB data than the power-law LCDM model; however, the improvement in the fit to the WMAP data is only Delta chi^2 = 3 for 1 extra degree of freedom. The combination of WMAP and other astronomical data yields significant constraints on the geometry of the universe, the equation of state of the dark energy, the gravitational wave energy density, and neutrino properties. Consistent with the predictions of simple inflationary theories, we detect no significant deviations from Gaussianity in the CMB maps.

The WMAP satellite has completed 3 years of observations of the cosmic microwave background radiation. The 3-year data products include several sets of full sky maps of the Stokes I, Q and U parameters in 5 frequency bands, spanning 23 to 94 GHz, and supporting items, such as beam window functions and noise covariance matrices. The processing used to produce the current sky maps and supporting products represents a significant advancement over the first year analysis, and is described herein. Improvements to the pointing reconstruction, radiometer gain modeling, window function determination and radiometer spectral noise parametrization are presented. A detailed description of the updated data processing that produces maximum likelihood sky map estimates is presented, along with the methods used to produce reduced resolution maps and corresponding noise covariance matrices. Finally two methods used to evaluate the noise of the full resolution sky maps are presented along with several representative year-to-year null tests, demonstrating that sky maps produced from data from different observational epochs are consistent.

The Wilkinson Microwave Anisotropy Probe WMAP has mapped the entire sky in five frequency bands between 23 and 94 GHz with polarization sensitive radiometers. We present three-year full-sky maps of the polarization and analyze them for foreground emission and cosmological implications. These observations open up a new window for understanding the universe. WMAP observes significant levels of polarized foreground emission due to both Galactic synchrotron radiation and thermal dust emission. The least contaminated channel is at 61 GHz. Informed by a model of the Galactic foreground emission, we subtract the foreground emission from the maps. In the foreground corrected maps, for l=2-6, we detect l(l+1) C_l^{EE} / (2 pi) = 0.086 +-0.029 microkelvin^2. This is interpreted as the result of rescattering of the CMB by free electrons released during reionization and corresponds to an optical depth of tau = 0.10 +- 0.03. We see no evidence for B-modes, limiting them to l(l+1) C_l^{BB} / (2 pi) = -0.04 +- 0.03 microkelvin^2. We find that the limit from the polarization signals alone is r<2.2 (95% CL) corresponding to a limit on the cosmic density of gravitational waves of Omega_{GW}h^2 < 5 times 10^{-12}. From the full WMAP analysis, we find r<0.55 (95% CL) corresponding to a limit of Omega_{GW}h^2 < 10^{-12} (95% CL).

Measurements of primordial non-Gaussianity ($f_{NL}$) open a new window onto the physics of inflation. We describe a fast cubic (bispectrum) estimator of $f_{NL}$, using a combined analysis of temperature and polarization observations. The speed of our estimator allows us to use a sufficient number of Monte Carlo simulations to characterize its statistical properties in the presence of real world issues such as instrumental effects, partial sky coverage, and foreground contamination. We find that our estimator is optimal, where optimality is defined by saturation of the Cramer Rao bound, if noise is homogeneous. Our estimator is also computationally efficient, scaling as $O(N^{3/2})$ compared to the $O(N^{5/2})$ scaling of the brute force bispectrum calculation for sky maps with $N$ pixels. For Planck this translates into a speed-up by factors of millions, reducing the required computing time from thousands of years to just hours and thus making $f_{NL}$ estimation feasible for future surveys. Our estimator in its current form is optimal if noise is homogeneous. In future work our fast polarized bispectrum estimator should be extended to deal with inhomogeneous noise in an analogous way to how the existing fast temperature estimator was generalized.

Blazars, a beamed population of active galactic nuclei, radiate high-energy gamma-rays, and thus are a good target for the Gamma Ray Large Area Space Telescope (GLAST). As the blazars trace the large-scale structure of the universe, one may observe spatial clustering of blazars. We calculate the angular power spectrum of blazars that would be detected by GLAST. We show that we have the best chance of detecting their clustering at large angular scales, \theta >~ 10 deg, where shot noise is less important, and the dominant contribution to the correlation comes from relatively low redshift, z <~ 0.1. The GLAST can detect the correlation signal, if the blazars detected by GLAST trace the distribution of low-z quasars observed by optical galaxy surveys, which have the bias of unity. If the bias of blazars is greater than 1.5, GLAST will detect the correlation signal unambiguously. We also find that GLAST may detect spatial clustering of clusters of galaxies in gamma-rays. The shape of the angular power spectrum is different for blazars and clusters of galaxies; thus, we can separate these two contributions on the basis of the shape of the power spectrum.

We derive analytical formulae for the Minkowski Functions of the cosmic microwave background (CMB) and large-scale structure (LSS) from primordial non-Gaussianity. These formulae enable us to estimate a non-linear coupling parameter, f_NL, directly from the CMB and LSS data without relying on numerical simulations of non-Gaussian primordial fluctuations. One can use these formulae to estimate statistical errors on f_NL from Gaussian realizations, which are much faster to generate than non-Gaussian ones, fully taking into account the cosmic/sampling variance, beam smearing, survey mask, etc. We show that the CMB data from the Wilkinson Microwave Anisotropy Probe should be sensitive to |f_NL|\simeq 40 at the 68% confidence level. The Planck data should be sensitive to |f_NL|\simeq 20. As for the LSS data, the late-time non-Gaussianity arising from gravitational instability and galaxy biasing makes it more challenging to detect primordial non-Gaussianity at low redshifts. The late-time effects obscure the primordial signals at small spatial scales. High-redshift galaxy surveys at z>2 covering \sim 10Gpc^3 volume would be required for the LSS data to detect |f_NL|\simeq 100. Minkowski Functionals are nicely complementary to the bispectrum because the Minkowski Functionals are defined in real space and the bispectrum is defined in Fourier space. This property makes the Minksowski Functionals a useful tool in the presence of real-world issues such as anisotropic noise, foreground and survey masks. Our formalism can be extended to scale-dependent f_NL easily.

We compare the non-linear matter power spectrum in real space calculated analytically from 3rd-order perturbation theory with N-body simulations at 1<z<6. We find that the perturbation theory prediction agrees with the simulations to better than 1% accuracy in the weakly non-linear regime where the dimensionless power spectrum, Delta^2(k)=k^3P(k)/2pi^2, which approximately gives variance of matter density field at a given k, is less than 0.4. While the baryonic acoustic oscillation features are preserved in the weakly non-linear regime at z>1, the shape of oscillations is distorted from the linear theory prediction. Nevertheless, our results suggest that one can correct the distortion caused by non-linearity almost exactly. We also find that perturbation theory, which does not contain any free parameters, provides a significantly better fit to the simulations than the conventional approaches based on empirical fitting functions to simulations. The future work would include perturbation theory calculations of non-linearity in redshift space distortion and halo biasing in the weakly non-linear regime.

(abridged) High-z galaxy redshift surveys open up exciting possibilities for precision determinations of neutrino masses and inflationary models. The high-z surveys are more useful for cosmology than low-z ones owing to much weaker non-linearities in matter clustering, redshift-space distortion and galaxy bias. We can then utilize the two-dimensional information of the linear power spectrum in angular and redshift space to measure the scale-dependent suppression of matter clustering due to neutrino free-streaming as well as the shape of the primordial power spectrum. To illustrate capabilities of high-z surveys for constraining neutrino masses and the primordial power spectrum, we compare three future redshift surveys covering 300 square degrees at 0.5<z<2, 2<z<4, and 3.5<z<6.5. We find that, combined with the cosmic microwave background data expected from the Planck satellite, these surveys allow precision determination of the total neutrino mass with the projected errors of sigma(m_nu)=0.059, 0.043, and 0.025 eV, respectively, thus yielding a positive detection of the neutrino mass rather than an upper limit, as sigma(m_nu) is smaller than the lower limits to the neutrino masses implied from the neutrino oscillation experiments. The accuracies of constraining the tilt and running index of the primordial power spectrum, sigma(n_s)=(3.8, 3.7, 3.0)x10^-3, and sigma(alpha_s)=(5.9, 5.7, 2.4)x10^-3, respectively, are smaller than the current uncertainties by more than an order of magnitude, which will allow us to discriminate between candidate inflationary models. In particular, the error on alpha_s from the highest redshift survey is not very far away from the prediction of a class of simple inflationary models driven by a massive scalar field with self-coupling, alpha_s=-(0.8-1.2)x10^-3.

We calculate the cosmic microwave background (CMB) angular trispectrum, spherical harmonic transform of the four-point correlation function, from primordial non-Gaussianity in primordial curvature perturbations characterized by a constant non-linear coupling parameter, $f_{\rm NL}$. We fully take into account the effect of the radiation transfer function, and thus provide the most accurate estimate of the signal-to-noise ratio of the angular trispectrum of CMB temperature anisotropy. We find that the predicted signal-to-noise ratio of the trispectrum summed up to a given $l$ is approximately a power-law, $(S/N)(<l)\sim 2.2\times 10^{-9}f^2_{\rm NL}l^2$, up to the maximum multipole that we have reached in our numerical calculation, $l=1200$, assuming that the error is dominated by cosmic variance. Our results indicate that the signal-to-noise ratio of the temperature trispectrum exceeds that of the bispectrum at the critical multipole, $l_c \sim 1500~(50/|f_{\rm NL}|)$. Therefore, the trispectrum of the Planck data is more sensitive to primordial non-Gaussianity than the bispectrum for $|f_{\rm NL}|\gtrsim 50$. We also report the predicted constraints on the amplitude of trispectrum, which may be useful for other non-Gaussian models such as curvaton models.

The redshifted ultraviolet light from early stars at z ~ 10 contributes to the cosmic near infrared background. We present detailed calculations of its spectrum with various assumptions about metallicity and mass spectrum of early stars. We show that if the near infrared background has a stellar origin, metal-free stars are not the only explanation of the excess near infrared background; stars with metals (e.g. Z=1/50 Z_sun) can produce the same amount of background intensity as the metal-free stars. We quantitatively show that the predicted average intensity at 1-2 microns is essentially determined by the efficiency of nuclear burning in stars, which is not very sensitive to metallicity. We predict \nu I_\nu / \dot{\rho}_* ~ 4-8 nW m^-2 sr^-1, where \dot{\rho_*} is the mean star formation rate at z=7-15 (in units of M_sun yr^-1 Mpc^-3) for stars more massive than 5 M_sun. On the other hand, since we have very little knowledge about the form of mass spectrum of early stars, uncertainty in the average intensity due to the mass spectrum could be large. An accurate determination of the near infrared background allows us to probe formation history of early stars, which is difficult to constrain by other means. While the star formation rate at z=7-15 inferred from the current data is significantly higher than the local rate at z<5, it does not rule out the stellar origin of the cosmic near infrared background. In addition, we show that a reasonable initial mass function, coupled with this star formation rate, does not over-produce metals in the universe in most cases, and may produce as little as less than 1 % of the metals observed in the universe today.

We show that the epoch(s) of reionization when the ionization fraction of the universe is about half can be determined by correlating Cosmic Microwave Background (CMB) temperature maps with 21-cm line maps at degree scales ($l\sim 100$). During reionization peculiar motion of free electrons induces the Doppler anisotropy of the CMB, while density fluctuations of neutral hydrogen induce the 21-cm line anisotropy. In our simplified model of inhomogeneous reionization, a positive correlation arises as the universe reionizes whereas a negative correlation arises as it recombines; thus, the sign of the correlation provides information on the reionization history which cannot be obtained by presently. The signal comes mainly from large scales (k~0.01 Mpc^-1) where linear perturbation theory is valid and complexity due to patchy reionization is averaged out. Since the Doppler signal comes from ionized regions and the 21-cm comes from neutral ones, the correlation has a well-defined peak(s) in redshift when the ionization fraction of the universe is about half. Furthermore, the cross-correlation is much less sensitive to systematic errors, especially foreground emission, than the auto-correlation of 21-cm lines: this is analogous to the temperature-polarization correlation of the CMB being more immune to systematic errors than the polarization-polarization. Therefore, we argue that the Doppler-21cm correlation provides a robust measurement of the 21-cm anisotropy, which can also be used as a diagnostic tool for detected signals in the 21-cm data -- detection of the cross-correlation provides the strongest confirmation that the signal is of cosmological origin. We show that the Square Kilometer Array can easily measure the predicted correlation signal for 1~year of survey observation.

We show that the energy density spectrum of the primordial gravitational waves has characteristic features due to the successive changes in the relativistic degrees of freedom during the radiation era. These changes make the evolution of radiation energy density deviate from the conventional adiabatic evolution, \rho_r~ a^{-4}, and thus cause the expansion rate of the universe to change suddenly at each transition which, in turn, modifies the spectrum of primordial gravitational waves. We take into account all the particles in the Standard Model of elementary particles. In addition, free-streaming of neutrinos damps the amplitude of gravitational waves, leaving characteristic features in the energy density spectrum. Our calculations are solely based on the standard model of cosmology and particle physics, and therefore these features must exist. Our calculations significantly improve the previous ones which ignored these effects and predicted a smooth, featureless spectrum.

High-energy photons from dark matter annihilation contribute to the cosmic gamma-ray background (CGB). Since dark matter particles are weakly interacting, annihilation can happen only in high density regions such as dark matter halos. The precise shape of the energy spectrum of CGB depends on the nature of dark matter particles, as well as the cosmological evolution of dark matter halos. In order to discriminate between the signals from dark matter annihilation and other astrophysical sources, however, the information from the energy spectrum may not be sufficient. We show that dark matter annihilation also produces a characteristic anisotropy of the CGB, which provides a powerful tool for testing the origin. We develop the formalism based on a halo model approach to calculate the three-dimensional power spectrum of dark matter clumping, which determines the power spectrum of annihilation signals. We show that the sensitivity of future gamma-ray detectors such as GLAST should allow us to measure the angular power spectrum of CGB anisotropy, if dark matter particles are supersymmetric neutralinos and they account for most of the observed mean intensity of CGB in GeV region. On the other hand, if dark matter has a relatively small mass, and accounts for most of the CGB in MeV region, then the future Advanced Compton Telescope should be able to measure the anisotropy in MeV region. As the intensity of photons from annihilation is proportional to the density squared, we show that the predicted shape of the angular power spectrum of gamma rays from dark matter annihilation is different from that due to other astrophysical sources such as blazars. Therefore, the angular power spectrum of the CGB provides a smoking-gun signature of dark matter annihilation.

The origin of the cosmic gamma-ray background at 1-20 MeV remains a mystery. We show that gamma-ray emission accompanying annihilation of 20 MeV dark matter particles explains most of the observed signal. Our model satisfies all of the current observational constraints, and naturally provides the origin of "missing" gamma-ray background at 1-20 MeV and 511 keV line emission from the Galactic center. We conclude that gamma-ray observations support the existence of 20 MeV dark matter particles. Improved measurements of the gamma-ray background in this energy band undoubtedly test our proposal.

Recent second-order perturbation computations have provided an accurate prediction for the primordial gravitational potential $\Phi(x)$ in scenarios in which cosmological perturbations are generated either during or after inflation. This enables us to make realistic predictions for a non-Gaussian part of $\Phi(x)$, which is generically written in momentum space as a double convolution of its Gaussian part with a suitable kernel, f_NL(k1,k2). This kernel defines the amplitude and angular structure of the non-Gaussian signals and originates from the evolution of second-order perturbations after the generation of the curvature perturbation. We derive a generic formula for the CMB angular bispectrum with arbitrary f_NL(k1,k2) and examine the detectability of the primordial non-Gaussian signals from various scenarios such as single-field inflation, inhomogeneous reheating, and curvaton scenarios. Our results show that in the standard slow-roll inflation scenario the signal actually comes from the momentum-dependent part of f_NL(k1,k2), and thus it is important to include the momentum dependence in the data analysis. In the other scenarios the primordial non-Gaussianity is comparable to or larger than these post-inflationary effects. We find that WMAP cannot detect non-Gaussian signals generated by these models. Numerical calculations for l>500 are still computationally expensive, and we are not yet able to extend our calculations to Planck's angular resolution; however, there is an encouraging trend which shows that Planck may be able to detect these non-Gaussian signals.

We derive a fast way for measuring primordial non-Gaussianity in a nearly full-sky map of the cosmic microwave background. We find a cubic combination of sky maps combining bispectrum configurations to capture a quadratic term in primordial fluctuations. Our method takes only N^1.5 operations rather than N^2.5 of the bispectrum analysis (1000 times faster for l=512), retaining the same sensitivity. A key component is a map of underlying primordial fluctuations, which can be more sensitive to the primordial non-Gaussianity than a temperature map. We also derive a fast and accurate statistic for measuring non-Gaussian signals from foreground point sources. The statistic is 10^6 times faster than the full bispectrum analysis, and can be used to estimate contamination from the sources. Our algorithm has been successfully applied to the Wilkinson Microwave Anisotropy Probe sky maps by Komatsu et al. (2003).

The observed cosmic gamma-ray background at ~ MeV has often been attributed to Type Ia supernovae (SNIa). Since SNIa is close to a standard candle, one can calculate the gamma-ray intensity of SNIa integrated over redshifts fairly accurately, once the evolution of the SNIa rate is known. The latest SNIa rate measured at z < ~1.6 (Dahlen et al. 2004) indicates that the previous calculations of the gamma-ray background consistently overestimated the SNIa rate. With the new rate, we find that the SNIa contribution is an order of magnitude smaller than observed, and thus new population(s) of sources should be invoked.

Motivated by a recent detection of 511 keV photons from the center of our Galaxy, we calculate the spectrum of the soft gamma-ray background of the redshifted 511 keV photons from cosmological halos. Annihilation of dark matter particles into electron-positron pairs makes a substantial contribution to the gamma-ray background. Mass of such dark matter particles must be <~ 100 MeV so that resulting electron-positron pairs are on-relativistic. On the other hand, we show that in order for the annihilation not to exceed the observed background, the dark matter mass needs to be >~ 20 MeV. We include the contribution from the active galactic nuclei and supernovae. The halo substructures may increase the lower bound to >~ 60 MeV.

This is a review of models of inflation and of their predictions for the primordial non-Gaussianity in the density perturbations which are thought to be at the origin of structures in the Universe. Non-Gaussianity emerges as a key observable to discriminate among competing scenarios for the generation of cosmological perturbations and is one of the primary targets of present and future Cosmic Microwave Background satellite missions. We give a detailed presentation of the state-of-the-art of the subject of non-Gaussianity, both from the theoretical and the observational point of view, and provide all the tools necessary to compute at second order in perturbation theory the level of non-Gaussianity in any model of cosmological perturbations. We discuss the new wave of models of inflation, which are firmly rooted in modern particle physics theory and predict a significant amount of non-Gaussianity. The review is addressed to both astrophysicists and particle physicists and contains useful tables which summarize the theoretical and observational results regarding non-Gaussianity.

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

The first year data from the Wilkinson Microwave Anisotropy Probe are used to place stringent constraints on the topology of the Universe. We search for pairs of circles on the sky with similar temperature patterns along each circle. We restrict the search to back-to-back circle pairs, and to nearly back-to-back circle pairs, as this covers the majority of the topologies that one might hope to detect in a nearly flat universe. We do not find any matched circles with radius greater than 25 degrees. For a wide class of models, the non-detection rules out the possibility that we live in a universe with topology scale smaller than 24 Gpc.

WMAP precision data enables accurate testing of cosmological models. We find that the emerging standard model of cosmology, a flat Lambda-dominated universe seeded by nearly scale-invariant adiabatic Gaussian fluctuations, fits the WMAP data. With parameters fixed only by WMAP data, we can fit finer scale CMB measurements and measurements of large scle structure (galaxy surveys and the Lyman alpha forest). This simple model is also consistent with a host of other astronomical measurements. We then fit the model parameters to a combination of WMAP data with other finer scale CMB experiments (ACBAR and CBI), 2dFGRS measurements and Lyman alpha forest data to find the model's best fit cosmological parameters: h=0.71+0.04-0.03, Omega_b h^2=0.0224+-0.0009, Omega_m h^2=0.135+0.008-0.009, tau=0.17+-0.06, n_s(0.05/Mpc)=0.93+-0.03, and sigma_8=0.84+-0.04. WMAP's best determination of tau=0.17+-0.04 arises directly from the TE data and not from this model fit, but they are consistent. These parameters imply that the age of the universe is 13.7+-0.2 Gyr. The data favors but does not require a slowly varying spectral index. By combining WMAP data with other astronomical data sets, we constrain the geometry of the universe, Omega_tot = 1.02 +- 0.02, the equation of state of the dark energy w < -0.78 (95% confidence limit assuming w >= -1), and the energy density in stable neutrinos, Omega_nu h^2 < 0.0076 (95% confidence limit). For 3 degenerate neutrino species, this limit implies that their mass is less than 0.23 eV (95% confidence limit). The WMAP detection of early reionization rules out warm dark matter.

We present limits to the amplitude of non-Gaussian primordial fluctuations in the WMAP 1-year cosmic microwave background sky maps. A non-linear coupling parameter, f_NL, characterizes the amplitude of a quadratic term in the primordial potential. We use two statistics: one is a cubic statistic which measures phase correlations of temperature fluctuations after combining all configurations of the angular bispectrum. The other uses the Minkowski functionals to measure the morphology of the sky maps. Both methods find the WMAP data consistent with Gaussian primordial fluctuations and establish limits, -58<f_NL<134, at 95% confidence. There is no significant frequency or scale dependence of f_NL. The WMAP limit is 30 times better than COBE, and validates that the power spectrum can fully characterize statistical properties of CMB anisotropy in the WMAP data to high degree of accuracy. Our results also validate the use of a Gaussian theory for predicting the abundance of clusters in the local universe. We detect a point-source contribution to the bispectrum at 41 GHz, b_src = (9.5+-4.4) X 1e-5 uK^3 sr^2, which gives a power spectrum from point sources of c_src = (15+-6) X 1e-3 uK^2 sr in thermodynamic temperature units. This value agrees well with independent estimates of source number counts and the power spectrum at 41 GHz, indicating that b_src directly measures residual source contributions.

We present full sky microwave maps in five bands (23 to 94 GHz) from the WMAP first year sky survey. Calibration errors are <0.5% and the low systematic error level is well specified. The 2<l<900 anisotropy power spectrum is cosmic variance limited for l<354 with a signal-to-noise ratio >1 per mode to l=658. The temperature-polarization cross-power spectrum reveals both acoustic features and a large angle correlation from reionization. The optical depth of reionization is 0.17 +/- 0.04, which implies a reionization epoch of 180+220-80 Myr (95% CL) after the Big Bang at a redshift of 20+10-9 (95% CL) for a range of ionization scenarios. This early reionization is incompatible with the presence of a significant warm dark matter density. The age of the best-fit universe is 13.7 +/- 0.2 Gyr old. Decoupling was 379+8-7 kyr after the Big Bang at a redshift of 1089 +/- 1. The thickness of the decoupling surface was dz=195 +/- 2. The matter density is Omega_m h^2 = 0.135 +0.008 -0.009, the baryon density is Omega_b h^2 = 0.0224 +/- 0.0009, and the total mass-energy of the universe is Omega_tot = 1.02 +/- 0.02. The spectral index of scalar fluctuations is fit as n_s = 0.93 +/- 0.03 at wavenumber k_0 = 0.05 Mpc^-1, with a running index slope of dn_s/d ln k = -0.031 +0.016 -0.018 in the best-fit model. This flat universe model is composed of 4.4% baryons, 22% dark matter and 73% dark energy. The dark energy equation of state is limited to w<-0.78 (95% CL). Inflation theory is supported with n_s~1, Omega_tot~1, Gaussian random phases of the CMB anisotropy, and superhorizon fluctuations. An admixture of isocurvature modes does not improve the fit. The tensor-to-scalar ratio is r(k_0=0.002 Mpc^-1)<0.90 (95% CL).

We confront predictions of inflationary scenarios with the WMAP data, in combination with complementary small-scale CMB measurements and large-scale structure data. The WMAP detection of a large-angle anti-correlation in the temperature--polarization cross-power spectrum is the signature of adiabatic superhorizon fluctuations at the time of decoupling. The WMAP data are described by pure adiabatic fluctuations: we place an upper limit on a correlated CDM isocurvature component. Using WMAP constraints on the shape of the scalar power spectrum and the amplitude of gravity waves, we explore the parameter space of inflationary models that is consistent with the data. We place limits on inflationary models; for example, a minimally-coupled lambda phi^4 is disfavored at more than 3-sigma using WMAP data in combination with smaller scale CMB and large scale structure survey data. The limits on the primordial parameters using WMAP data alone are: n_s(k_0=0.002 Mpc^{-1})=1.20_{-0.11}^{+0.12}, dn/dlnk=-0.077^{+0.050}_{- 0.052}, A(k_0=0.002 Mpc}^{-1})=0.71^{+0.10}_{-0.11} (68% CL), and r(k_0=0.002 Mpc^{-1})<1.28 (95% CL).

We present the angular power spectrum derived from the first-year Wilkinson Microwave Anisotropy Probe (WMAP) sky maps. We study a variety of power spectrum estimation methods and data combinations and demonstrate that the results are robust. The data are modestly contaminated by diffuse Galactic foreground emission, but we show that a simple Galactic template model is sufficient to remove the signal. Point sources produce a modest contamination in the low frequency data. After masking ~700 known bright sources from the maps, we estimate residual sources contribute ~3500 uK^2 at 41 GHz, and ~130 uK^2 at 94 GHz, to the power spectrum l*(l+1)*C_l/(2*pi) at l=1000. Systematic errors are negligible compared to the (modest) level of foreground emission. Our best estimate of the power spectrum is derived from 28 cross-power spectra of statistically independent channels. The final spectrum is essentially independent of the noise properties of an individual radiometer. The resulting spectrum provides a definitive measurement of the CMB power spectrum, with uncertainties limited by cosmic variance, up to l~350. The spectrum clearly exhibits a first acoustic peak at l=220 and a second acoustic peak at l~540 and it provides strong support for adiabatic initial conditions. Kogut et al. (2003) analyze the C_l^TE power spectrum, and present evidence for a relatively high optical depth, and an early period of cosmic reionization. Among other things, this implies that the temperature power spectrum has been suppressed by \~30% on degree angular scales, due to secondary scattering.

The angular power spectrum of the SZ effect, C_l, is a powerful probe of cosmology. It is easier to detect than individual clusters in the field, is insensitive to observational selection effects and does not require a calibration between cluster mass and flux, reducing the systematic errors which dominate the cluster-counting constraints. It receives a dominant contribution from cluster region between 20-40% of the virial radius and is thus insensitive to the poorly known gas physics in the cluster centre, such as cooling or (pre)heating. In this paper we derive a refined analytic prediction for C_l using the universal gas-density and temperature profile and the dark-matter halo mass function. The predicted C_l has no free parameters and fits all of the published hydrodynamic simulation results to better than a factor of two around l=3000. We find that C_l scales as (sigma_8)^7 times (Omega_b h)^2 and is almost independent of all of the other cosmological parameters. This differs from the local cluster abundance studies, which give a relation between sigma_8 and Omega_m. We also compute the covariance matrix of C_l using the halo model and find a good agreement relative to the simulations. We estimate how well we can determine sigma_8 with sampling-variance-limited observations and find that for a several-square-degree survey with 1-2 arcminute resolution one should be able to determine sigma_8 to within a few percent, with the remaining uncertainty dominated by theoretical modelling. If the recent excess of the CMB power on small scales reported by the CBI experiment is due to the SZ effect, then we find sigma_8(Omega_b h/0.029)^0.3 = 1.04 +- 0.12 at the 95% confidence level (statistical) and with a residual 10% systematic (theoretical) uncertainty.

We present theoretical and observational studies of non-Gaussian fluctuations in CMB, by using the angular bispectrum and trispectrum. We predict the primary angular bispectrum from inflation, and forecast how well we can measure the primordial non-Gaussian signal. In addition to that, secondary anisotropy sources in the low-redshift universe also produce non-Gaussianity, so do foreground emissions from extragalactic or interstellar microwave sources. We study how well we can measure these non-Gaussian signals, including the primordial signal. We find that slow-roll inflation produces too small bispectrum to be detected by any experiments; thus, any detection strongly constrains this class of models. We also find that the secondary bispectrum from coupling between the SZ effect and the weak lensing effect, and the foreground bispectrum from extragalactic point sources, give detectable non-Gaussian signals on small angular scales. We test Gaussianity of the COBE DMR sky maps, by measuring all the modes of the angular bispectrum down to the DMR beam size. We find no significant signal of the bispectrum. We also find that the previously reported detection of the bispectrum is consistent with a statistical fluctuation. By fitting the theoretical prediction to the data for the primary bispectrum, we put a constraint on non-linearity in inflation. We conclude that the angular bispectrum finds no significant non-Gaussian signals in the DMR data. We present the first measurement of the angular trispectrum on the DMR sky maps, further testing Gaussianity of the DMR data. We find no significant non-Gaussian signals in the trispectrum. Therefore, the angular bispectrum and trispectrum show that the DMR sky map is comfortably consistent with Gaussianity.

We have performed a combined analysis of X-ray and Sunyaev-Zel'dovich data in the direction of the distant galaxy cluster, RX J2228+2037. Fitting a $\beta$-model to the high-resolution HRI data gives $r_c = 103 \pm 12 h_{70}^{-1}$ kpc and $\beta=0.54 \pm 0.03$. The dependency of the Sunyaev-Zel'dovich effect with respect to the gas temperature allows us, through the additional use of the 21 GHz data of the cluster, to determine $k_B T_e=10.4 \pm 1.8 h_{70}^{1/2}$ keV. Extrapolating the gas density profile out to the virial radius ($R_v=r_{178}=2.9$ Mpc), we derived a gas mass of $M_{g}(r<R_v)=(4.0\pm 0.2)\times 10^{14} h_{70}^{-5/2} \rm{M}_\odot$. Within the hypothesis of hydrostatic equilibrium, the corresponding extrapolated total mass for this source is: $M_{tot}(r<R_v)=(1.8 \pm 0.4)\times 10^{15} h^{-1} \rm{M}_\odot$, which corresponds to a gas fraction of $f_{gas}=0.22\pm 0.06 h_{70}^{-3/2}$. Our results on the temperature and on the cluster mass classify RX J2228+2037 among the distant, hot and very massive galaxy clusters. Our work highlights the power of the association of galaxy cluster mapping observations in X-ray and the SZ effect to derive the cluster's physical properties, even without X-ray spectroscopy.

We measure the angular bispectrum of the cosmic microwave background (CMB) radiation anisotropy from the COBE Differential Microwave Radiometer (DMR) four-year sky maps. The angular bispectrum is the harmonic transform of the three-point correlation function, analogous to the angular power spectrum, the harmonic transform of the two-point correlation function. First, we study statistical properties of the bispectrum and the normalized bispectrum. We find the latter more useful for statistical analysis; the distribution of the normalized bispectrum is very much Gaussian, while the bare bispectrum distribution is highly non-Gaussian. Then, we measure 466 modes of the normalized bispectrum, all independent combinations of three-point configurations up to a maximum multipole of 20, the mode corresponding to the DMR beam size. By measuring 10 times as many modes as the sum of previous work, we test Gaussianity of the DMR maps. We compare the data with the simulated Gaussian realizations, finding no significant detection of the normalized bispectrum on the mode-by-mode basis. We also find that the previously reported detection of the normalized bispectrum is consistent with a statistical fluctuation. By fitting a theoretical prediction to the data for the primary CMB bispectrum, which is motivated by slow-roll inflation, we put a weak constraint on a parameter characterizing non-linearity in inflation. Simultaneously fitting the foreground bispectra estimated from interstellar dust and synchrotron template maps shows that neither dust nor synchrotron emissions significantly contribute to the bispectrum at high Galactic latitude. We conclude that the DMR map is consistent with Gaussianity.

We present an analytic approach to predict gas density and temperature profiles in dark matter haloes. We assume that the gas density profile traces the dark matter density profile in outer parts of the haloes, as suggested by many hydrodynamic simulations. Under this assumption, the hydrostatic equilibrium uniquely fixes the two free parameters in the model, the mass-temperature normalization and the polytropic index, that determine the temperature profile. This enables us to predict a universal gas profile from any universal dark matter profile. Our results show that gas tracing dark matter in the outer parts of the haloes is inconsistent with gas being isothermal; on the contrary, it requires temperature to decrease slowly with radius, in agreement with observations. We compare our predictions for X-ray surface brightness profiles of the haloes and the mass-temperature relation with observations. We find that they are generally in a good agreement. We compare the universal profile with the beta profile and find that, although the beta profile gives a reasonable fit to our predicted profiles, the deviation from it can explain many of the observed trends, once we take into account the observational selection effects. Our model predicts that the mass-temperature relation does not follow the self-similar relation because of the mass-dependent halo concentration. We also predict surface brightness profiles of the Sunyaev-Zel'dovich (SZ) effect. We find that fitted to the beta profile the core radii and beta inferred from the SZ effect are systematically larger than those from the X-ray measurement.

Why is non-Gaussianity interesting? One of generic predictions from inflationary scenarios is that primordial fluctuations are exactly Gaussian in linear order; however, the non-linearity in the inflation will produce weak non-Gaussianity. Thus, measuring the non-Gaussianity in the cosmic microwave background radiation anisotropy is a probe of the non-linear physics in the very early universe. Since the angular three point function is zero for the Gaussian field, it is sensitive to the non-Gaussianity. We predict its harmonic transform counterpart, the angular {\it Bispectrum}, down to arcminutes angular scales, including the full effect of the radiation transfer function. We find that even the Planck experiment cannot detect the primary bispectrum from the inflation, as long as the single field slow-roll inflation is right. Non-linearities in the low redshift universe also produce the non-Gaussianity. We find that secondary bispectra are detectable by both MAP and Planck experiments. The secondary bispectra probe the non-linear physics of the low-redshift universe. Although this could be a contaminant to the primary signal, MAP and Planck experiments are found to be able to separate the primary from secondary effects well. We present a tentative comparison of the primary bispectrum to the published COBE 4 year bispectrum. The data put a weak constraint on the parameter, and the constraint would become much tighter when we use all modes available in the COBE data, and certainly forthcoming satellite experiments. As a conclusion, the bispectrum is a key measure to confirm or destroy the simple inflationary scenario in non-linear order that seems quite successful in linear order.

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

If the primordial fluctuations are non-Gaussian, then this non-Gaussianity will be apparent in the cosmic microwave background (CMB) sky. With their sensitive all-sky observation, MAP and Planck satellites should be able to detect weak non-Gaussianity in the CMB sky. On large angular scale, there is a simple relationship between the CMB temperature and the primordial curvature perturbation. On smaller scales; however, the radiation transfer function becomes more complex. In this paper, we present the angular bispectrum of the primary CMB anisotropy that uses the full transfer function. We find that the bispectrum has a series of acoustic peaks that change a sign, and a period of acoustic oscillations is twice as long as that of the angular power spectrum. Using a single non-linear coupling parameter to characterize the amplitude of the bispectrum, we estimate the expected signal-to-noise ratio for COBE, MAP, and Planck experiments. We find that the detection of the primary bispectrum by any kind of experiments should be problematic for the simple slow-roll inflationary scenarios. We compare the sensitivity of the primary bispectrum to the primary skewness and conclude that when we can compute the predicted form of the bispectrum, it becomes a ``matched filter'' for detecting the non-Gaussianity in the data, and much more powerful tool than the skewness. We also show that MAP and Planck can separate the primary bispectrum from various secondary bispectra on the basis of the shape difference. The primary CMB bispectrum is a test of the inflationary scenario, and also a probe of the non-linear physics in the very early universe.

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

We investigate the weak gravitational lensing effect due to the large-scale structure of the universe on two-point correlations of local maxima ({\em hotspots}) in the 2D sky map of the cosmic microwave background (CMB) anisotropy. According to the Gaussian random statistics as most inflationary scenarios predict, the hotspots are discretely distributed with some {\em characteristic} angular separations on the last scattering surface owing to oscillations of the CMB angular power spectrum. The weak lensing then causes pairs of hotspots which are separated with the characteristic scale to be observed with various separations. We found that the lensing fairly smoothes the oscillatory features of the two-point correlation function of hotspots. This indicates that the hotspots correlations can be a new statistical tool for measuring shape and normalization of the power spectrum of matter fluctuations from the lensing signatures.

(Abridged) The hot gas in the IGM produces anisotropies in the Cosmic Microwave Background (CMB) through the thermal Sunyaev-Zel'dovich (SZ) effect. The SZ effect is a powerful probe of large-scale structure in the universe and must be carefully subtracted from measurements of the primary CMB anisotropies. We use moving-mesh hydrodynamical simulations to study the 3-dimensional statistics of the gas, and compute the mean comptonization parameter and the angular power spectrum of the SZ fluctuations, for different cosmologies. We compare these results with predictions using the Press-Schechter formalism. We find that the two methods agree approximately, but differ in details. We discuss this discrepancy, and show that resolution limits the reliability of our results to the 200<l<2000 range. For cluster- normalized CDM models, the SZ power spectrum is comparable to the primordial power spectrum around l=2000. We show that groups and filaments (kT<5 keV) contribute about 50% of the SZ power spectrum at l=500. About half of the SZ power spectrum on these scales is produced at redshifts z<0.1, and can thus be detected and removed using existing catalogs of galaxies and X-ray clusters. We discuss the implications of these results for the future MAP and Planck Surveyor missions.

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

Cosmological implications of clusters of galaxies are discussed with particular attention to their importance in probing the cosmological parameters. More specifically we compute the number counts of clusters of galaxies, Log $N$ -- Log $S$ relation, in X-ray and submm bands on the basis of the Press--Schechter theory. As an important step toward breaking the degeneracy among the viable cosmological models, we observed the most luminous X-ray cluster RXJ1347-1145 in three bands (21 and 43 GHz in the Nobeyama Radio Observatory, Japan, and 350 GHz in the J. C. Maxwell telescope at Mauna Kea, Hawaii). We report on the preliminary results which are in good agreement with the profile of the Sunayev -- Zel'dovich effect predicted on the basis of the previous X-ray observation of the cluster.

We report on the detection of the Sunyaev -- Zel'dovich (SZ) signals toward the most luminous X-ray cluster RXJ1347-1145 at Nobeyama Radio Observatory (21 and 43 GHz) and at James Clerk Maxwell Telescope (350 GHz). In particular the latter is the first successful detection of the SZ temperature increment in the submillimeter band which resolved the profile of a cluster of galaxies. Both the observed spectral dependence and the radial profile of the SZ signals are fully consistent with those expected from the X-ray observation of the cluster. The combined analysis of 21GHz and 350GHz data reproduces the temperature and core-radius of the cluster determined with the ROSAT and ASCA satellites when we adopt the slope of the density profile from the X-ray observations. Therefore our present data provide the strongest and most convincing case for the detection of the submillimeter SZ signal from the cluster, as well as in the Rayleigh -- Jeans regime. We also discuss briefly the cosmological implications of the present results.

We present complete constraints imposed from observations of the cosmic microwave background radiation (CMBR) on the chaotic inflationary scenario with a nonminimally coupled inflaton field proposed by Fakir and Unruh (FU). Our constraints are complete in the sense that we investigate both the scalar density perturbation and the tensor gravitational wave in the Jordan frame, as well as in the Einstein frame. This makes the constraints extremely strong without any ambiguities due to the choice of frames. We find that the FU scenario generates tiny tensor contributions to the CMBR relative to chaotic models in minimal coupling theory, in spite of its spectral index of scalar perturbation being slightly tilted. This means that the FU scenario will be excluded if any tensor contributions to CMBR are detected by the forthcoming satellite missions. Conversely, if no tensor nature is detected despite the tilted spectrum, a minimal chaotic scenario will be hard to explain and the FU scenario will be supported.

We investigate a possibility to restrict the chaotic inflationary scenario with "large" nonminimally coupled inflaton("phi") field considered by Fakir and Unruh by means of the observed CMBR anisotropy. We calculate the contributions of the long wavelength gravitational waves generated in the inflationary period to CMBR anisotropy quadrupole moment. We obtain the constraint lambda/(xi^2) = 1.8 X 10^(-8)/psi_i, where "lambda" is the self coupling, "xi" is the nonminimal coupling constant and "psi_i" means the initial value of psi = 8piG(phi^2). Combining this with previously obtained constraint that 2.1 times squre root of lambda/(xi^2) = dT/T = 1.1 X 10^(-5), we conclude that the initial value has to be psi_i = 1.6 X 10^3. If the self-coupling has a reasonably values of order 10^(-2), then xi = 10^4 and phi_i = 10^(-1) times Planck mass.