We present a Next-to-next-to-leading (fifth or NNLO) order calculation for the covariance matrix of the matter power spectrum, taking into account the effect of survey window functions. Using the grid-based calculation scheme for the standard perturbation theory, GridSPT, we quickly generate multiple realizations of the nonlinear density fields to fifth order in perturbation theory, then estimate the power spectrum and the covariance matrix from the sample. To the end, we have obtained the non-Gaussian covariance originated from the one-loop trispectrum without explicitly computing the trispectrum. By comparing the GridSPT calculations with the N-body results, we show that NNLO GridSPT result reproduces the N-body results on quasi-linear scales where SPT accurately models nonlinear matter power spectrum. Incorporating the survey window function effect to GridSPT is rather straightforward, and the resulting NNLO covariance matrix also matches well with the N-body results.

We present a Fisher information study of the statistical impact of galaxy bias and selection effects on the estimation of key cosmological parameters from galaxy redshift surveys; in particular, the angular diameter distance, Hubble parameter, and linear growth rate at a given redshift, the cold dark matter density, and the tilt and running of the primordial power spectrum. The line-of-sight-dependent selection contributions we include here are known to exist in real galaxy samples. We determine the maximum wavenumber included in the analysis by requiring that the next-order corrections to the galaxy power spectrum or bispectrum, treated here at next-to-leading and leading order, respectively, produce shifts of $\lesssim 0.25\sigma$ on each of the six cosmological parameters. With the galaxy power spectrum alone, selection effects can deteriorate the constraints severely, especially on the linear growth rate. Adding the galaxy bispectrum helps break parameter degeneracies significantly. We find that a joint power spectrum-bispectrum analysis of a Euclid-like survey can still measure the linear growth rate to 10% precision after complete marginalization over selection bias. We also discuss systematic parameter shifts arising from ignoring selection effects and/or other bias parameters, and emphasize that it is necessary to either control selection effects at the percent level or marginalize over them. We obtain similar results for the Roman Space Telescope and HETDEX.

We present a fast implementation of the next-to-leading order (1-loop) redshift-space galaxy power spectrum by using FFTlog-based methods. [V. Desjacques, D. Jeong, and F. Schmidt, JCAP 1812 (12), 035] have shown that the 1-loop galaxy power spectrum in redshift space can be computed with 28 independent loop integrals with 22 bias parameters. Analytical calculation of the angular part of the loop integrals leaves the radial part in the form of a spherical Bessel transformation that is ready to be integrated numerically by using the FFTLog transformation. We find that the original 28 loop integrals can be solved with a total of 85 unique FFTLog transformations, yet leading to a few orders of magnitude speed up over traditional multi-dimensional integration. The code used in this work is publicly available at https://github.com/JosephTomlinson/GeneralBiasPk

Certain anomalies in the CMB bring out a tension between the six-parameter flat $\Lambda$CDM model and the CMB data. We revisit the PLANCK analysis with loop quantum cosmology (LQC) predictions and show that LQC alleviates both the large-scale power anomaly and the tension in the lensing amplitude. These differences arise because, in LQC, the primordial power spectrum is scale dependent for small $k$, with a specific power suppression. We conclude with a prediction of larger optical depth and power suppression in the B-mode polarization power spectrum on large scales.

The discrepancy between the Hubble parameter inferred from local measurements and that from the cosmic microwave background (CMB) has motivated careful scrutiny of the assumptions that enter both analyses. Here we point out that the location of the recombination peak in the CMB B-mode power spectrum is determined by the light horizon at the surface of last scatter and thus provides an alternative early-Universe standard ruler. It can thus be used as a cross-check for the standard ruler inferred from the acoustic peaks in the CMB temperature power spectrum and to test various explanations for the Hubble tension. The measurement can potentially be carried out with a precision of $\lesssim2\%$ with stage-IV B-mode experiments. The measurement can also be used to measure the propagation speed of gravitational waves in the early Universe.

We consider a new technique for neutrino energy and topology reconstruction at DUNE. In particular, we show that when the direction of the incoming neutrino is known, one can use the measured directions of the outgoing leptonic and hadronic particles to reconstruct poorly-measured quantities, such as the hadronic cascade energy. We show that this alternative technique yields an energy resolution which is comparable to current reconstruction methods which sum measured energies. As a proof of concept we apply this new reconstruction method to a search for dark matter annihilation in the Sun. We show that the use of directional information from both the leptonic and hadronic interaction products allows one to effectively reject backgrounds and isolate the signal, giving competitive sensitivities.

The galaxy bispectrum contains a wealth of information about the early universe, gravity, as well as astrophysics such as galaxy bias. In this paper, we study the parity-odd part of the galaxy bispectrum which is hitherto unexplored. In the standard cosmological model, the odd-parity bispectrum is generated by galaxy velocities through redshift-space distortions. While small in the case of General Relativity coupled with smooth dark energy, the signal could be larger in modified gravity scenarios. Thus, apart from being a very useful consistency test of measurements of galaxy clustering, the odd bispectrum offers a novel avenue for searching for new physics.

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.

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.

Nonlinearities in the gravitational evolution, galaxy bias, and redshift-space distortion drive the observed galaxy density fields away from the initial near-Gaussian states. Exploiting such a non-Gaussian galaxy density field requires measuring higher-order correlation functions, or, its Fourier counterpart, polyspectra. Here, we present an efficient parallel algorithm for estimating higher-order polyspectra. Based upon the Scoccimarro estimator, the estimator avoids direct sampling of polygons by using the Fast-Fourier Transform (FFT), and the parallelization overcomes the large memory requirement of the original estimator. In particular, we design the memory layout to minimize the inter-CPU communications, which excels in the code performance.

By utilizing large-scale graph analytic tools implemented in the modern Big Data platform, Apache Spark, we investigate the topological structure of gravitational clustering in five different universes produced by cosmological $N$-body simulations with varying parameters: (1) a WMAP 5-year compatible $\Lambda$CDM cosmology, (2) two different dark energy equation of state variants, and (3) two different cosmic matter density variants. For the Big Data calculations, we use a custom build of stand-alone Spark/Hadoop cluster at Korea Institute for Advanced Study (KIAS) and Dataproc Compute Engine in Google Cloud Platform (GCP) with the sample size ranging from 7 millions to 200 millions. We find that among the many possible graph-topological measures, three simple ones: (1) the average of number of neighbors (the so-called average vertex degree) $\alpha$, (2) closed-to-connected triple fraction (the so-called transitivity) $\tau_\Delta$, and (3) the cumulative number density $n_{s\ge5}$ of subcomponents with connected component size $s \ge 5$, can effectively discriminate among the five model universes. Since these graph-topological measures are in direct relation with the usual $n$-points correlation functions of the cosmic density field, graph-topological statistics powered by Big Data computational infrastructure opens a new, intuitive, and computationally efficient window into the dark Universe.

This review presents a comprehensive overview of galaxy bias, that is, the statistical relation between the distribution of galaxies and matter. We focus on large scales where cosmic density fields are quasi-linear. On these scales, the clustering of galaxies can be described by a perturbative bias expansion, and the complicated physics of galaxy formation is absorbed by a finite set of coefficients of the expansion, called bias parameters. The review begins with a detailed derivation of this very important result, which forms the basis of the rigorous perturbative description of galaxy clustering, under the assumptions of General Relativity and Gaussian, adiabatic initial conditions. Key components of the bias expansion are all leading local gravitational observables, which include the matter density but also tidal fields and their time derivatives. We hence expand the definition of local bias to encompass all these contributions. This derivation is followed by a presentation of the peak-background split in its general form, which elucidates the physical meaning of the bias parameters, and a detailed description of the connection between bias parameters and galaxy statistics. We then review the excursion-set formalism and peak theory which provide predictions for the values of the bias parameters. In the remainder of the review, we consider the generalizations of galaxy bias required in the presence of various types of cosmological physics that go beyond pressureless matter with adiabatic, Gaussian initial conditions: primordial non-Gaussianity, massive neutrinos, baryon-CDM isocurvature perturbations, dark energy, and modified gravity. Finally, we discuss how the description of galaxy bias in the galaxies' rest frame is related to clustering statistics measured from the observed angular positions and redshifts in actual galaxy catalogs.

We present the complete expression for the next-to-leading (1-loop) order galaxy power spectrum and the leading-order galaxy bispectrum in redshift space in the general bias expansion, or equivalently the effective field theory of biased tracers. We consistently include all line-of-sight dependent selection effects. These are degenerate with many, but not all, of the redshift-space distortion contributions, and have not been consistently derived before. Moreover, we show that, in the framework of effective field theory, a consistent bias expansion in redshift space must include these selection contributions. Physical arguments about the tracer sample considered and its observational selection have to be used to justify neglecting the selection contributions. In summary, the next-to-leading order galaxy power spectrum and leading-order galaxy bispectrum in the general bias expansion are described by 22 parameters, which reduces to 11 parameters if selection effects can be neglected. All contributions to the power spectrum can be written in terms of 28 independent loop integrals.

We explore the possible spectrum of binary mergers of sub-solar mass black holes formed out of dark matter particles interacting via a dark electromagnetism. We estimate the properties of these dark black holes by assuming that their formation process is parallel to Population-III star formation; except that dark molecular cooling can yield smaller opacity limit. We estimate the binary coalescence rates for the Advanced LIGO and Einstein telescope, and find that scenarios compatible with all current constraints could produce dark black holes at rates high enough for detection by Advanced LIGO.

Perturbation theory (PT) calculation of large-scale structure has been used to interpret the observed non-linear statistics of large-scale structure at the quasi-linear regime. In particular, the so-called standard perturbation theory (SPT) provides a basis for the analytical computation of the higher-order quantities of large-scale structure. Here, we present a novel, grid-based algorithm for the SPT calculation, hence named GridSPT, to generate the higher-order density and velocity fields from a given linear power spectrum. Taking advantage of the Fast Fourier Transform, the GridSPT quickly generates the nonlinear density fields at each order, from which we calculate the statistical quantities such as non-linear power spectrum and bispectrum. Comparing the density fields (to fifth order) from GridSPT with those from the full N-body simulations in the configuration space, we find that GridSPT accurately reproduces the N-body result on large scales. The agreement worsens with smaller smoothing radius, particularly for the under-dense regions where we find that 2LPT (second-order Lagrangian perturbation theory) algorithm performs well.

The merger of dark matter halos and the gaseous structures embedded in them, such as proto-galaxies, galaxies, and groups and clusters of galaxies, results in strong shocks that are capable of accelerating cosmic rays (CRs) to $\sim10~\rm PeV$. These shocks will produce high-energy neutrinos and $\gamma$-rays through inelastic $pp$ collisions with ambient gaseous environments. In this work, we study the contributions of these halo mergers to the diffuse neutrino flux measured in IceCube and to the non-blazar portion of the extragalactic $\gamma$-ray background measured by $Fermi$. In order to calculate them, we formulate the redshift dependence of the shock velocity, galactic radius, halo gas content and galactic/intergalactic magnetic fields over the dark matter halo distribution up to a redshift $z=10$. We find that high-redshift mergers contribute a significant amount of the cosmic-ray energy luminosity density, and the resulting neutrino spectra could explain a large part of the observed diffuse neutrino flux above 0.1 PeV up to several PeV. We also show that our model can somewhat alleviate tensions with the extragalactic $\gamma$-ray background. First, since a larger fraction of the CR energy luminosity density comes from high redshifts, the accompanying $\gamma$-rays are more strongly suppressed through $\gamma\gamma$ annihilations with the cosmic microwave background (CMB) and the extragalactic background light (EBL). Second, mildly radiative-cooled shocks may lead to a harder CR spectrum with spectral indices of $1.5\lesssim s\lesssim2.0$. Our study suggests that halo mergers, a fraction of which may also induce starbursts in the merged galaxies, can be promising neutrino emitters without violating the existing $Fermi$ $\gamma$-ray constraints on the non-blazar component of the extragalactic $\gamma$-ray background.

We present the 2-point function from Fast and Accurate Spherical Bessel Transformation (2-FAST) algorithm for a fast and accurate computation of integrals involving one or two spherical Bessel functions. These types of integrals occur when projecting the galaxy power spectrum $P(k)$ onto the configuration space, $\xi_\ell^\nu(r)$, or spherical harmonic space, $C_\ell(\chi,\chi')$. First, we employ the FFTlog transformation of the power spectrum to divide the calculation into $P(k)$-dependent coefficients and $P(k)$-independent integrations of basis functions multiplied by spherical Bessel functions. We find analytical expressions for the latter integrals in terms of special functions, for which recursion provides a fast and accurate evaluation. The algorithm, therefore, circumvents direct integration of highly oscillating spherical Bessel functions.

We delineate the conditions under which the consistency relation for the non-Gaussian bias and the universality of the halo mass function hold in the context of microscopic Lagrangian descriptions of halos. The former is valid provided that the collapse barrier depends only on the physical fields (instead of fields normalized by their variance for example) and explicitly includes the effect of {\it all} physical fields such as the tidal shear. The latter holds provided that the response of the halo number density to a long-wavelength density fluctuation is equivalent to the response induced by shifting the spherical collapse threshold. Our results apply to any Lagrangian halo bias prescription. Effective "moving" barriers, which are ubiquitous in the literature, do not generally satisfy the consistency relation. Microscopic barriers including the tidal shear lead to two additional, second-order Lagrangian bias parameters which ensure that the consistency relation is satisfied. We provide analytic expressions for them.

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.

We have measured the probability distribution function (PDF) of cosmic matter density field from a suite of N-body simulations. We propose the generalized normal distribution of version 2 (Nv2) as an alternative fitting formula to the well-known log-normal distribution. We find that Nv2 provides significantly better fit than the log-normal distribution for all smoothing radii (2, 5, 10, 25 [Mpc/h]) that we studied. The improvement is substantial in the underdense regions. The development of non- Gaissianities in the cosmic matter density field is captured by continuous evolution of the skewness and shifts parameters of the Nv2 distribution. We present the redshift evolution of these parameters for aforementioned smoothing radii and various background cosmology models. All the PDFs measured from large and high-resolution N-body simulations that we use in this study can be obtained from a Web site at https://astro.kias.re.kr/jhshin.

A tensor-type cosmological perturbation, defined as a transverse and traceless spatial fluctuation, is often interpreted as the gravitational waves. While decoupled from the scalar-type perturbations in linear order, the tensor perturbations can be sourced from the scalar-type in the nonlinear order. The tensor perturbations generated by the quadratic combination of linear scalar-type cosmological perturbation are widely studied in the literature, but all previous studies are based on zero-shear gauge without proper justification. Here, we show that, being second order in perturbation, such an induced tensor perturbation is generically gauge dependent. In particular, the gravitational wave power spectrum depends on the hypersurface (temporal gauge) condition taken for the linear scalar perturbation. We further show that, during the matter-dominated era, the induced tensor modes dominate over the linearly evolved primordial gravitational waves amplitude for $k\gtrsim10^{-2}~[h/{\rm Mpc}]$ even for the gauge that gives lowest induced tensor modes with the optimistic choice of primordial gravitational waves ($r=0.1$). The induced tensor modes, therefore, must be modeled correctly specific to the observational strategy for the measurement of primordial gravitational waves from large-scale structure via, for example, parity-odd mode of weak gravitational lensing, or clustering fossils.

We develop and study the position-dependent bispectrum. It is a generalization of the recently proposed position-dependent power spectrum method of measuring the squeezed-limit bispectrum. The position-dependent bispectrum can similarly be used to measure the squeezed-limit trispectrum in which one of the wavelengths is much longer than the other three. In this work, we will mainly consider the case in which the three smaller wavelengths are nearly the same (the equilateral configuration). We use the Fisher information matrix to forecast constraints on bias parameters and the amplitude of primordial trispectra from the position-dependent bispectrum method. We find that the method can constrain the local-type $g_{\rm NL}$ at a level of $\sigma(g_{\rm NL}^{\rm local}) \approx 3 \times 10^5$ for a large volume SPHEREx-like survey; improvements can be expected by including all the triangular configurations of the bispectra rather than just the equilateral configuration. However, the same measurement would also constrain a much larger family of trispectra than local $g_{\rm NL}$ model. We discuss the implications of the forecasted reach of future surveys in terms of super cosmic variance uncertainties from primordial non-Gaussianities.

We present the leading order non-linear density and velocity power spectra in the complete form; previous studies have omitted the vector- and tensor-type perturbations simultaneously excited by the scalar-type perturbation in nonlinear order. These additional contributions are comparable to the scalar-type purely relativistic perturbations, and thus negligible in the current paradigm of concordance cosmology: i.e., concerning density and velocity perturbations of the pressureless matter in perturbation regime well inside of matter-dominated epoch, we show that pure Einstein's gravity contributions appearing from the third order are entirely negligible (five orders of magnitude smaller than the Newtonian contributions) in all scales. We thus prove that Newtonian perturbation theory is quite reliable in calculating the amplitude of matter fluctuations even in the precision era of cosmology. Therefore, besides the ones imprinted as the initial condition generated in the earlier phase, the other relativistic effect relevant for interpreting observational data must be the projection effect that occurs when mapping galaxies on to the observed coordinate.

We study the parity-odd part (that we shall call Doppler term) of the linear galaxy two-point correlation function that arises from wide-angle, velocity, Doppler lensing and cosmic acceleration effects. As it is important at low redshift and at large angular separations, the Doppler term is usually neglected in the current generation of galaxy surveys. For future wide-angle galaxy surveys such as Euclid, SPHEREx and SKA, however, we show that the Doppler term must be included. The effect of these terms is dominated by the magnification due to relativistic aberration effects and the slope of the galaxy redshift distribution and it generally mimics the effect of the local type primordial non-Gaussianity with the effective nonlinearity parameter $f_{\rm NL}^{\rm eff}$ of a few, we show that this would affect forecasts on measurements of $f_{\rm NL}$ at low-redshift. Our results show that a survey at low redshift with large number density over a wide area of the sky could detect the Doppler term with a signal-to-noise ratio of $\sim 1-20$, depending on survey specifications.

Accurate modeling of nonlinearities in the galaxy bispectrum, the Fourier transform of the galaxy three-point correlation function, is essential to fully exploit it as a cosmological probe. In this paper, we present numerical and theoretical challenges in modeling the nonlinear bispectrum. First, we test the robustness of the matter bispectrum measured from N-body simulations using different initial conditions generators. We run a suite of N-body simulations using the Zel'dovich approximation and second-order Lagrangian perturbation theory (2LPT) at different starting redshifts, and find that transients from initial decaying modes systematically reduce the nonlinearities in the matter bispectrum. To achieve 1% accuracy in the matter bispectrum for $z\le3$ on scales $k<1$ $h$/Mpc, 2LPT initial conditions generator with initial redshift of $z\gtrsim100$ is required. We then compare various analytical formulas and empirical fitting functions for modeling the nonlinear matter bispectrum, and discuss the regimes for which each is valid. We find that the next-to-leading order (one-loop) correction from standard perturbation theory matches with N-body results on quasi-linear scales for $z\ge1$. We find that the fitting formula in Gil-Mar\'{\i}n et al. (2012) accurately predicts the matter bispectrum for $z\le1$ on a wide range of scales, but at higher redshifts, the fitting formula given in Scoccimarro & Couchman (2001) gives the best agreement with measurements from N-body simulations.

Inflation models make specific predictions for a tensor-scalar-scalar three-point correlation, or bispectrum, between one gravitational-wave (tensor) mode and two density-perturbation (scalar) modes. This tensor-scalar-scalar correlation leads to a local power quadrupole, an apparent departure from statistical isotropy in our Universe, as well as characteristic four-point correlations in the current mass distribution in the Universe. So far, the predictions for these observables have been worked out only for single-clock models in which certain consistency conditions between the tensor-scalar-scalar correlation and tensor and scalar power spectra are satisfied. Here we review the requirements on inflation models for these consistency conditions to be satisfied. We then consider several examples of inflation models, such as non-attractor and solid inflation models, in which these conditions are put to the test. In solid inflation the simplest consistency conditions are already violated whilst in the non-attractor model we find that, contrary to the standard scenario, the tensor-scalar-scalar correlator probes directly relevant model-dependent information. We work out the predictions for observables in these models. For non-attractor inflation we find an apparent local quadrupolar departure from statistical isotropy in large-scale structure but that this power quadrupole decreases very rapidly at smaller scales. The consistency of the CMB quadrupole with statistical isotropy then constrains the distance scale that corresponds to the transition from the non-attractor to attractor phase of inflation to be larger than the currently observable horizon. Solid inflation predicts clustering fossils signatures in the current galaxy distribution that may be large enough to be detectable with forthcoming, and possibly even current, galaxy surveys.

The statistics of primordial curvature fluctuations are our window into the period of inflation, where these fluctuations were generated. To date, the cosmic microwave background has been the dominant source of information about these perturbations. Large scale structure is however from where drastic improvements should originate. In this paper, we explain the theoretical motivations for pursuing such measurements and the challenges that lie ahead. In particular, we discuss and identify theoretical targets regarding the measurement of primordial non-Gaussianity. We argue that when quantified in terms of the local (equilateral) template amplitude $f_{\rm NL}^{\rm loc}$ ($f_{\rm NL}^{\rm eq}$), natural target levels of sensitivity are $\Delta f_{\rm NL}^{\rm loc, eq.} \simeq 1$. We highlight that such levels are within reach of future surveys by measuring 2-, 3- and 4-point statistics of the galaxy spatial distribution. This paper summarizes a workshop held at CITA (University of Toronto) on October 23-24, 2014.

We study the non-linear corrections to the matter and velocity power spectra in the synchronous gauge (SG). We consider the perturbations up to third order in a zero-pressure fluid in flat cosmological background, which is relevant for the non-linear growth of cosmic structure. As a result, we point out that the SG is an inappropriate coordinate choice when handling the non-linear growth of the large-scale structure. Although the equations in the SG happen to coincide with those in the comoving gauge (CG) to linear order, they differ from second order. In particular, the second order hydrodynamic equations in the the SG are apparently in the Lagrangian form, whereas those in the CG are in the Eulerian form. Thus, the non-linear power spectra naively presented in the original SG show strange behavior quite different from the result of the Newtonian theory even on sub-horizon scales. The power spectra in the SG show regularized behaviors only after we introduce convective terms in the second order so that the equations in two gauges coincide to the second order.

Future galaxy surveys will measure baryon acoustic oscillations (BAOs) with high significance, and a complete understanding of the anisotropies of BAOs in redshift space will be important to exploit the cosmological information in BAOs. Here we describe the anisotropies that arise in the redshift-space galaxy two-point correlation function (2PCF) and elucidate the origin of features that arise in the dependence of the BAOs on the angle between the orientation of the galaxy pair and the line of sight. We do so with a derivation of the configuration-space 2PCF using streaming model. We find that, contrary to common belief, the locations of BAO peaks in the redshift-space 2PCF are anisotropic even in the linear theory. Anisotropies in BAO depend strongly on the method of extracting the peak, showing maximum 3 % angular variation. We also find that extracting the BAO peak of $r^2\xi(r,\mu)$ significantly reduces the anisotropy to sub-percent level angular variation. When subtracting the tilt due to the broadband behavior of the 2PCF, the BAO bump is enhanced along the line of sight because of local infall velocities toward the BAO bump. Precise measurement of the angular dependence of the redshift-space 2PCF will allow new geometrical tests of dark energy beyond the BAO.

We present a data-driven method to infer the redshift distribution of an arbitrary dataset based on spatial cross-correlation with a reference population and we apply it to various datasets across the electromagnetic spectrum to show its potential and limitations. Our approach advocates the use of clustering measurements on all available scales, in contrast to previous works focusing only on linear scales. We also show how its accuracy can be enhanced by optimally sampling a dataset within its photometric space rather than applying the estimator globally. We show that the ultimate goal of this technique is to characterize the mapping between the space of photometric observables and redshift space as this characterization then allows us to infer the clustering-redshift p.d.f. of a single galaxy. We apply this technique to estimate the redshift distributions of luminous red galaxies and emission line galaxies from the SDSS, infrared sources from WISE and radio sources from FIRST. We show that consistent redshift distributions are found using both quasars and absorber systems as reference populations. This technique brings valuable information on the third dimension of astronomical datasets. It is widely applicable to a large range of extra-galactic surveys.

We review recent studies that rigorously define several key observables of the large-scale structure of the Universe in a general relativistic context. Specifically, we consider i) redshift perturbation of cosmic clock events; ii) distortion of cosmic rulers, including weak lensing shear and magnification; iii) observed number density of tracers of the large-scale structure. We provide covariant and gauge-invariant expressions of these observables. Our expressions are given for a linearly perturbed flat Friedmann-Robertson-Walker metric including scalar, vector, and tensor metric perturbations. While we restrict ourselves to linear order in perturbation theory, the approach can be straightforwardly generalized to higher order.

The recent indications for a tensor-to-scalar ratio $ r \simeq 0.2$ from BICEP2 measurements of the cosmic microwave background (CMB) B-mode polarization present some tension with upper limits $r \lesssim 0.1$ from measurements of CMB temperature fluctuations. Here we point out that tensor perturbations can induce B modes in the CMB polarization without inducing any temperature fluctuations nor E-mode polarization whatsoever, but only, at the expense of violating the Copernican principle. We present this mathematical possibility as a new ingredient for the model-builder's toolkit in case the tension between B modes and temperature fluctuations cannot be resolved with more conventional ideas.

For the simplest inflation models, the recent detection of a large primordial B-mode polarization signal by the BICEP2 experiment indicates a slight tension with the upper limit on the tensor-to-scalar ratio, $r$, from the Planck satellite. Here, we discuss spatially varying $r$ as a possible explanation for this discrepancy. This idea seems attractive since it may also explain part of the hemispherical temperature power asymmetry seen by WMAP and Planck at large angular scales. If these two aspects are indeed connected, the model suggests that in the Northern hemisphere $r$ should be much smaller, a hypothesis that could be confirmed with future B-mode experiments, providing a test for the stationarity of primordial tensor contributions across the sky. The BICEP2 measurement furthermore rules out that a simple dipolar modulation of $r$ alone can be responsible for the full hemispherical power asymmetry.

We study the relationship between dark-matter haloes and matter in the MIP $N$-body simulation ensemble, which allows precision measurements of this relationship, even deeply into voids. What enables this is a lack of discreteness, stochasticity, and exclusion, achieved by averaging over hundreds of possible sets of initial small-scale modes, while holding fixed large-scale modes that give the cosmic web. We find (i) that dark-matter-halo formation is greatly suppressed in voids; there is an exponential downturn at low densities in the otherwise power-law matter-to-halo density bias function. Thus, the rarity of haloes in voids is akin to the rarity of the largest clusters, and their abundance is quite sensitive to cosmological parameters. The exponential downturn appears both in an excursion-set model, and in a model in which fluctuations evolve in voids as in an open universe with an effective $\Omega_m$ proportional to a large-scale density. We also find that (ii) haloes typically populate the average halo-density field in a super-Poisson way, i.e. with a variance exceeding the mean; and (iii) the rank-order-Gaussianized halo and dark-matter fields are impressively similar in Fourier space. We compare both their power spectra and cross-correlation, supporting the conclusion that one is roughly a strictly-increasing mapping of the other. The MIP ensemble especially reveals how halo abundance varies with `environmental' quantities beyond the local matter density; (iv) we find a visual suggestion that at fixed matter density, filaments are more populated by haloes than clusters.

We study the dissipation of small-scale adiabatic perturbations at early times when the Universe is hotter than T ~ 0.5 keV. When the wavelength falls below the damping scale 1/kD, the acoustic modes diffuse and thermalize, causing entropy production. Before neutrino decoupling, kD is primarily set by the neutrino shear viscosity, and we study the effect of acoustic damping on the relic neutrino number, primordial nucleosynthesis, dark-matter freeze-out, and baryogenesis. This sets a new limit on the amplitude of primordial fluctuations of DeltaR^2 < 0.007 at 10^4/Mpc< k < 10^5/Mpc and a model dependent limit of DeltaR^2 < 0.3 at k < 10^{20-25}/Mpc.

We show that the relativistic effects are negligibly small in the non-linear density and velocity bispectra. Although the non-linearities of Einstein equation introduce additional non-linear terms to the Newtonian fluid equations, the corrections to the bispectrum only show up on super-horizon scales. We show this with the next-to-leading order non-linear bispectrum for a pressureless fluid in a flat Friedmann-Robertson-Walker background, by calculating the density and velocity fields up to fourth order. We work in the comoving gauge, where the dynamics is identical to the Newtonian up to second order. We also discuss the leading order matter bispectrum in various gauges, and show yet another relativistic effect near horizon scales that the matter bispectrum strongly depends on the gauge choice.

Departures of the Cosmic Microwave Background (CMB) frequency spectrum from a blackbody - commonly referred to as spectral distortions - encode information about the thermal history of the early Universe (redshift z < few x 10^6). While the signal is usually characterized as mu- and y-type distortion, a smaller residual (non-y/non-mu) distortion can also be created at intermediate redshifts 10^4 < z < 3 x 10^5. Here, we construct a new set of observables, mu_k, that describes the principal components of this residual distortion. The principal components are orthogonal to temperature shift, y- and mu-type distortion, and ranked by their detectability, thereby delivering a compression of all valuable information offered by the CMB spectrum. This method provides an efficient way of analyzing the spectral distortion for given experimental settings, and can be applied to a wide range of energy-release scenarios. As an illustration, we discuss the analysis of the spectral distortion signatures caused by dissipation of small-scale acoustic waves and decaying/annihilating particles for a PIXIE-type experiments. We provide forecasts for the expected measurement uncertainties of model parameters and detections limits in each case. We furthermore show that a PIXIE-type experiments can in principle distinguish dissipative energy release from particle decays for a nearly scale-invariant primordial power spectrum with small running. Future CMB spectroscopy thus offers a unique probe of physical processes in the primordial Universe.

The B-(curl-)mode of the correlation of galaxy ellipticities (shear) can be used to detect a stochastic gravitational wave background, such as that predicted by inflation. In this paper, we derive the tensor mode contributions to shear from both gravitational lensing and intrinsic alignments, using the gauge-invariant, full-sky results of arXiv:1204.3625. We find that the intrinsic alignment contribution, calculated using the linear alignment model, is larger than the lensing contribution by an order of magnitude or more, if the alignment strength for tensor modes is of the same order as for scalar modes. This contribution also extends to higher multipoles. These results make the prospects for probing tensor modes using galaxy surveys less pessimistic than previously thought, though still very challenging.

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.

Our motion relative to the cosmic-microwave-background (CMB) rest frame deflects light rays giving rise to shifts as large as L -> L(1+beta), where beta=0.00123 is our velocity (in units of the speed of light) on measurements CMB fluctuations. Here we present a novel harmonic-space approach to this CMB aberration that improves upon prior work by allowing us to (i) go to higher orders in beta, thus extending the validity of the analysis to measurements at L > 1/beta ~ 800; and (ii) treat the effects of window functions and pixelization in a more accurate and computationally efficient manner. We calculate precisely the magnitude of the systematic bias in the power spectrum inferred from the partial sky, and show that aberration shifts the multipole moment by Delta L/L ~ beta<cos(theta)>, with <cos(theta)> averaged over the survey footprint. Such a shift, if ignored, would bias the measurement of the sound-horizon size theta_* at the 0.01%-level, which is comparable to the measurement uncertainties of Planck. The bias can then propagate into cosmological parameters such as the angular-diameter distance, Hubble parameter and dark-energy equation of state. We study the effect of aberration for current Planck, South Pole Telescope (SPT) and Atacama Cosmology Telescope (ACT) data and show that the bias cannot be neglected. We suggest that the small tension between Planck, ACT, and SPT may be due partially to aberration. An Appendix shows how the near constancy of the full-sky power spectrum under aberration follows from unitarity of the aberration kernel.

Inflationary models predict a correlation between primordial density perturbations (scalar metric perturbations) and gravitational waves (tensor metric perturbations) in the form of a scalar-scalar-tensor three-point correlation, or bispectrum in Fourier space. The squeezed limit of this bispectrum implies a quadrupolar asymmetry in the observed local power spectrum for matter and galaxies. Here we show (like others before) that an infrared divergence in the amplitude of this power asymmetry predicted in single-field slow-roll models is canceled by projection effects when considering the observed power spectrum. We then further evaluate the nonzero, but finite, residual quadrupolar power asymmetry that remains after the divergences are canceled. While the quadrupolar power asymmetry is small, it is conceptually important. Our calculation moreover clarifies how the predictions for this power asymmetry may change in models with different scalar-scalar-tensor bispectra, and shows that convincing detection of the quadrupolar power asymmetry would rule out the single-field slow-roll models of inflation.

Thanks to the rapid progress in precision cosmology in the last few years, we now have access to physical observables that may constrain the theory of inflation through the non-Gaussianity (NG) signatures in the cosmic microwave background radiation and the distribution of large-scale structure. Numerical modeling of the NG signals from different inflation models is essential to correctly interpret current and near future data from large-scale structure surveys. In this study, we use high-resolution cosmological hydrodynamical simulations to investigate the effects of primordial NG on the formation and evolution of galaxies from the cosmic dawn to the present day. Focusing on the local type primordial NG, we find that it may affect the formation history of stars and black holes in galaxies, and their distribution. Compared to the Gaussian case, large non-Gaussian potential with $f_{NL} \gtrsim 10^3$ leads to earlier collapse of the first structures, more massive galaxies especially at high redshifts, stronger clustering of galaxies, and higher halo bias. However, for smaller NG with $f_{NL} \lesssim 10^2$, the effect is significantly weaker. Observations of the distribution and properties of high-redshift, rare objects such as the first galaxies and quasars may provide further constraints on the primordial NG.

We present a derivation of two-point correlations of general tracers in the peak-background split (PBS) framework by way of a rigorous definition of the PBS argument. Our expressions only depend on connected matter correlators and "renormalized" bias parameters with clear physical interpretation, and are independent of any coarse-graining scale. This result should be contrasted with the naive expression derived from a local bias expansion of the tracer number density with respect to the matter density perturbation \delta_L coarse-grained on a scale R_L. In the latter case, the predicted tracer correlation function receives contributions of order <\delta_L^n> at each perturbative order n, whereas, in our formalism, these are absorbed in the PBS bias parameters at all orders. Further, this approach naturally predicts both a scale-dependent bias ~ k^2 such as found for peaks of the density field, and the scale-dependent bias induced by primordial non-Gaussianity in the initial conditions. The only assumption made about the tracers is that their abundance at a given position depends solely on the matter distribution within a finite region around that position.

If during inflation the inflaton couples to a "fossil" field, some new scalar, vector, or tensor field, it typically induces a scalar-scalar-fossil bispectrum. Even if the fossil field leaves no direct physical trace after inflation, it gives rise to correlations between different Fourier modes of the curvature or, equivalently, a nonzero curvature trispectrum, but without a curvature bispectrum. Here we quantify the effects of a fossil field on the cosmic microwave background (CMB) temperature fluctuations in terms of bipolar spherical harmonics (BiPoSHs). The effects of vector and tensor fossils can be distinguished geometrically from those of scalars through the parity of the BiPoSHs they induce. However, the two-dimensional nature of the CMB sky does not allow vectors to be distinguished geometrically from tensors. We estimate the detectability of a signal in terms of the scalar-scalar-fossil coupling for scalar, vector, and tensor fossils, assuming a local-type coupling. We comment on a divergence that arises in the quadrupolar BiPoSH from the scalar-scalar-tensor correlation in single-field slow-roll inflation.

We derive general covariant expressions for the six independent observable modes of distortion of ideal standard rulers in a perturbed Friedmann-Robertson-Walker spacetime. Our expressions are gauge-invariant and valid on the full sky. These six modes are most naturally classified in terms of their rotational properties on the sphere, yielding two scalars, two vector (spin-1), and two tensor (spin-2) components. One scalar corresponds to the magnification, while the spin-2 components correspond to the shear. The vector components allow for a polar/axial decomposition analogous to the E/B-decomposition for the shear. Scalar modes do not contribute to the axial (B-)vector, opening a new avenue to probing tensor modes. Our results apply, but are not limited to, the distortion of correlation functions (of the CMB, 21cm emission, or galaxies) as well as to weak lensing shear and magnification, all of which can be seen as methods relying on "standard rulers".

In a perturbed Universe, comoving tracers on a two-dimensional surface of constant observed redshift are at different proper time since the Big Bang. For tracers whose age is known independently, one can measure these perturbations of the proper time. Examples of such sources include cosmic events which only happen during a short period of cosmic history, as well as evolving standard candles and standard rulers. In this paper we derive a general gauge-invariant linear expression for this perturbation in terms of space-time perturbations. As an example, we show that the observed temperature perturbations of the cosmic microwave background (CMB) on large scales are exactly given by these proper time perturbations. Together with the six ruler perturbations derived in Schmidt and Jeong (2012), this completes the set of independent observables which can be measured with standard rulers and candles.

Physical models for the hemispherical power asymmetry in the cosmic microwave background (CMB) reported by the Planck Collaboration must satisfy CMB constraints to the homogeneity of the Universe and quasar constraints to power asymmetries. We survey a variety of models for the power asymmetry and show that consistent models include a modulated scale-dependent isocurvature contribution to the matter power spectrum or a modulation of the reionization optical depth, gravitational-wave amplitude, or scalar spectral index. We propose further tests to distinguish between the different scenarios.

Total-angular-momentum (TAM) waves provide a set of basis functions for scalar, vector, and tensor fields that can be used in place of plane waves and that reflect the rotational symmetry of the spherical sky. Here we discuss three-point correlation functions, or bispectra in harmonic space, for scalar, vector, and tensor fields in terms of TAM waves. The Wigner-Eckart theorem dictates that the expectation value, assuming statistical isotropy, of the product of three TAM waves is the product of a Clebsch-Gordan coefficient (or Wigner-3j symbol) times a function only of the total-angular-momentum quantum numbers. Here we show how this works, and we provide explicit expressions relating the bispectra for TAM waves in terms of the more commonly used Fourier-space bispectra. This formalism will be useful to simplify calculations of projections of three-dimensional bispectra onto the spherical sky.

Most calculations in cosmological perturbation theorydecompose those perturbations into plane waves (Fourier modes). However, for some calculations, particularly those involving observations performed on a spherical sky, a decomposition into waves of fixed total angular momentum (TAM) may be more appropriate. Here we introduce TAM waves, solutions of fixed total angular momentum to the Helmholtz equation, for three-dimensional scalar, vector, and tensor fields. The vector TAM waves of given total angular momentum can be decomposed further into a set of three basis functions of fixed orbital angular momentum (OAM), a set of fixed helicity, or a basis consisting of a longitudinal (L) and two transverse (E and B) TAM waves. The symmetric traceless rank-2 tensor TAM waves can be similarly decomposed into a basis of fixed OAM or fixed helicity, or a basis that consists of a longitudinal (L), two vector (VE and VB, of opposite parity), and two tensor (TE and TB, of opposite parity) waves. We show how all of the vector and tensor TAM waves can be obtained by applying derivative operators to scalar TAM waves. This operator approach then allows one to decompose a vector field into three covariant scalar fields for the L, E, and B components and symmetric-traceless-tensor fields into five covariant scalar fields for the L, VE, VB, TE, and TB components. We provide projections of the vector and tensor TAM waves onto vector and tensor spherical harmonics. We provide calculational detail to facilitate the assimilation of this formalism into cosmological calculations. As an example, we calculate the power spectra of the deflection angle for gravitational lensing by density perturbations and by gravitational waves. We comment on an alternative approach to CMB fluctuations based on TAM waves. Our work may have applications elsewhere in field theory and in general relativity.

Characterizing the level of primordial non-Gaussianity (PNG) in the initial conditions for structure formation is one of the most promising ways to test inflation and differentiate among different scenarios. The scale-dependent imprint of PNG on the large-scale clustering of galaxies and quasars has already been used to place significant constraints on the level of PNG in our observed Universe. Such measurements depend upon an accurate and robust theory of how PNG affects the bias of galactic halos relative to the underlying matter density field. We improve upon previous work by employing a more general analytical method - the path-integral extension of the excursion set formalism - which is able to account for the non-Markovianity caused by PNG in the random-walk model used to identify halos in the initial density field. This non-Markovianity encodes information about environmental effects on halo formation which have so far not been taken into account in analytical bias calculations. We compute both scale-dependent and -independent corrections to the halo bias, along the way presenting an expression for the conditional collapsed fraction for the first time, and a new expression for the conditional halo mass function. To leading order in our perturbative calculation, we recover the halo bias results of Desjacques et. al. (2011), including the new scale-dependent correction reported there. However, we show that the non-Markovian dynamics from PNG can lead to marked differences in halo bias when next-to-leading order terms are included. We quantify these differences here. [abridged]

Observed angular positions and redshifts of large-scale structure tracers such as galaxies are affected by gravitational waves through volume distortion and magnification effects. Thus, a gravitational wave background can in principle be probed through clustering statistics of large-scale structure. We calculate the observed angular clustering of galaxies in the presence of a gravitational wave background at linear order including all relativistic effects. For a scale-invariant spectrum of gravitational waves, the effects are most significant at the smallest multipoles (2 <= l <= 5), but typically suppressed by six or more orders of magnitude with respect to scalar contributions for currently allowed amplitudes of the inflationary gravitational wave background. We also discuss the most relevant second-order terms, corresponding to the distortion of tracer correlation functions by gravitational waves. These provide a natural application of the approach recently developed in arXiv:1204.3625.

Several recent studies have shown how to properly calculate the observed clustering of galaxies in a relativistic context, and uncovered corrections to the Newtonian calculation that become significant on scales near the horizon. Here, we retrace these calculations and show that, on scales approaching the horizon, the observed galaxy power spectrum depends strongly on which gauge is assumed to relate the intrinsic fluctuations in galaxy density to matter perturbations through a linear bias relation. Starting from simple physical assumptions, we derive a gauge-invariant expression relating galaxy density perturbations to matter density perturbations on large scales, and show that it reduces to a linear bias relation in synchronous-comoving gauge, corroborating an assumption made in several recent papers. We evaluate the resulting observed galaxy power spectrum, and show that it leads to corrections similar to an effective non-Gaussian bias corresponding to a local (effective) fNL < 0.5. This number can serve as a guideline as to which surveys need to take into account relativistic effects. We also discuss the scale-dependent bias induced by primordial non-Gaussianity in the relativistic context, which again is simplest in synchronous-comoving gauge.

Many inflationary theories introduce new scalar, vector, or tensor degrees of freedom that may then affect the generation of primordial density perturbations. Here we show how to search a galaxy (or 21-cm) survey for the imprint of primordial scalar, vector, and tensor fields. These new fields induce local departures to an otherwise statistically isotropic two-point correlation function, or equivalently, nontrivial four-point correlation functions (or trispectra, in Fourier space), that can be decomposed into scalar, vector, and tensor components. We write down the optimal estimators for these various components and show how the sensitivity to these modes depends on the galaxy-survey parameters. New probes of parity-violating early-Universe physics are also presented.

The large-scale clustering of galaxies can serve as a probe of primordial non-Gaussianity in the Universe competitive with the anisotropies of the CMB. Here, we present improved theoretical predictions which include an important, previously overlooked correction to the bias. We demonstrate that the new predictions are able to reproduce the results of N-body simulations, explaining the significant departures seen from previous theoretical results. These refined predictions open the way to accurate constraints on primordial physics with large-scale structure surveys.

Recent results of N-body simulations have shown that current theoretical models are not able to correctly predict the amplitude of the scale-dependent halo bias induced by primordial non-Gaussianity, for models going beyond the simplest, local quadratic case. Motivated by these discrepancies, we carefully examine three theoretical approaches based on (1) the statistics of thresholded regions, (2) a peak-background split method based on separation of scales, and (3) a peak-background split method using the conditional mass function. We first demonstrate that the statistics of thresholded regions, which is shown to be equivalent at leading order to a local bias expansion, cannot explain the mass-dependent deviation between theory and N-body simulations. In the two formulations of the peak-background split on the other hand, we identify an important, but previously overlooked, correction to the non-Gaussian bias that strongly depends on halo mass. This new term is in general significant for any primordial non-Gaussianity going beyond the simplest local fNL model. In a separate paper, we compare these new theoretical predictions with N-body simulations, showing good agreement for all simulated types of non-Gaussianity.

Estimators for weak lensing observables such as shear and convergence generally have non-linear corrections, which, in principle, make weak lensing power spectra sensitive to primordial non-Gaussianity. In this paper, we quantitatively evaluate these contributions for weak lensing auto- and cross-correlation power spectra, and show that they are strongly suppressed by projection effects. This is a consequence of the central limit theorem, which suppresses departures from Gaussianity when the projection reaches over several correlation lengths of the density field, L_P~55 [Mpc/h]. Furthermore, the typical scales that contribute to projected bispectra are generally smaller than those that contribute to projected power spectra. Both of these effects are not specific to lensing, and thus affect the statistics of non-linear tracers (e.g., peaks) of any projected density field. Thus, the clustering of biased tracers of the three-dimensional density field is generically more sensitive to non-Gaussianity than observables constructed from projected density fields.

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]

Non-linear nature of Einstein equation introduces genuine relativistic higher order corrections to the usual Newtonian fluid equations describing the evolution of cosmological perturbations. We study the effect of such novel non-linearities on the next-to-leading order matter and velocity power spectra for the case of pressureless, irrotational fluid in a flat Friedmann background. We find that pure general relativistic corrections are negligibly small over all scales. Our result guarantees that, in the current paradigm of standard cosmology, one can safely use Newtonian cosmology even in non-linear regimes.

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.

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 probe the pure Einstein's gravity contributions to the second-order density power spectrum. In the small-scale, we discover that the Einstein's gravity contribution is negligibly small. This guarantees that Newton's gravity is sufficient to handle the baryon acoustic oscillation scale. In the large scale, however, we discover that the Einstein's gravity contribution to the second-order power spectrum dominates the linear-order power spectrum. Thus, pure Einstein gravity contribution appearing in the third-order perturbation leads to an infrared divergence in the power spectrum.

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 study the evolution of the curvature perturbation on the super-horizon scales starting from the inflationary epoch until there remains only a single dynamical degree of freedom, presureless matter, in the universe. We consider the cosmic inflation driven by a multiple number of the inflaton fields, which decay into both radiation and pressureless matter components. We present a complete set of the exact background and perturbation equations which describe the evolution of the universe throughout its history. By applying these equations to the simple but reasonable model of multi-field chaotic inflation, we explicitly show that the total curvature perturbation is continuously varying because of the non-adiabatic components of the curvature perturbation generated by the multiple inflaton fields throughout the whole evolution of the universe. We also provide an useful analytic estimation of the total as well as matter and radiation curvature perturbations, assuming that matter is completely decoupled from radiation from the beginning. The resulting isocurvature perturbation between matter and radiation is at most sub-percent level when the masses of the inflaton fields are distributed between $10^{-6}m_{pl}$ and $10^{-5}m_{pl}$. We find that this result is robust unless we use non-trivial decay rates, and that thus, in general, it is hard to obtain large matter-radiation isocurvature perturbation. Also, by using the $\delta{N}$ formalism, we point out that the inflationary calculation, especially when involving multiple inflaton fields, is likely to lose the potentially important post-inflationary evolution which can modify the resulting curvature perturbation.

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