Source author record

Donghui Jeong

Donghui Jeong appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

Catalog footprint

What is connected

58works
8topics
4close collaborators

Actions

Connect this record

Log in to claim

Research graph

See the researcher in context

Open full explorer

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

Published work

58 published item(s)

preprint2026arXiv

A New Collision Avoidance Fiber Assignment Algorithm for Robotic Fiber Positioners in Multi-Object Spectroscopy

We present a new fiber assignment algorithm for a robotic fiber positioner system in multi-object spectroscopy. Modern fiber positioner systems typically have overlapping patrol regions, resulting in the number of observable targets being highly dependent on the fiber assignment scheme. To maximize observable targets without fiber collisions, the algorithm proceeds in three steps. First, it assigns the maximum number of targets for a given field of view without considering any collisions between fiber positioners. Then, the fibers in collision are grouped, and the algorithm finds the optimal solution resolving the collision problem within each group. We compare the results from this new algorithm with those from a simple algorithm that assigns targets in descending order of their rank by considering collisions. As a result, we could increase the overall completeness of target assignments by 10% with this new algorithm in comparison with the case using the simple algorithm in a field with 150 fibers. Our new algorithm is designed for the All-sky SPECtroscopic survey of nearby galaxies (A-SPEC) based on the K-SPEC spectrograph system, but can also be applied to similar fiber-based systems with heavily overlapping fiber positioners.

preprint2024arXiv

Absorption Troughs of Lyman Alpha Emitters in HETDEX

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) is designed to detect and measure the redshifts of more than one million Ly$α$ emitting galaxies (LAEs) between $1.88 < z < 3.52$. In addition to its cosmological measurements, these data enable studies of Ly$α$ spectral profiles and the underlying radiative transfer. Using the roughly half a million LAEs in the HETDEX Data Release 3, we stack various subsets to obtain the typical Ly$α$ profile for the $z \sim 2-3$ epoch and to understand their physical properties. We find clear absorption wings around Ly$α$ emission, which extend $\sim 2000$ km $\mathrm{s}^{-1}$ both redward and blueward of the central line. Using far-UV spectra of nearby ($0.002 < z < 0.182$) LAEs in the CLASSY treasury and optical/near-IR spectra of $2.8 < z < 6.7$ LAEs in the MUSE-Wide survey, we observe absorption profiles in both redshift regimes. Dividing the sample by volume density shows that the troughs increase in higher density regions. This trend suggests that the depth of the absorption is dependent on the local density of objects near the LAE, a geometry that is similar to damped Lyman-$α$ systems. Simple simulations of Ly$α$ radiative transfer can produce similar troughs due to absorption of light from background sources by HI gas surrounding the LAEs.

preprint2022arXiv

Molecular Chemistry for Dark Matter

Molecular cooling is essential for studying the formation of sub-structure of dissipative dark-matter halos that may host compact objects such as black holes. Here, we analyze the reaction rates relevant for the formation, dissociation, and transition of hydrogenic molecules while allowing for different values of the physical parameters: the coupling constant, the proton mass, and the electron mass. For all cases, we re-scale the reaction rates for the standard molecular hydrogen, so our results are valid as long as the dark matter is weakly coupled and one of the fermions is much heavier than the other. These results will allow a robust numerical treatment of cosmic structure, in particular for mini-halos for which molecular cooling is important, in a dissipative dark matter scenario.

preprint2022arXiv

Molecular Chemistry for Dark Matter II: Recombination, Molecule Formation, and Halo Mass Function in Atomic Dark Matter

Dissipative dark matter predicts rich observable phenomena that can be tested with future large-scale structure surveys. As a specific example, we study atomic dark matter, consisting of a heavy particle and a light particle charged under a dark electromagnetism. In particular, we calculate the cosmological evolution of atomic dark matter focusing on dark recombination and dark-molecule formation. We have obtained the relevant interaction-rate coefficients by re-scaling the rates for normal hydrogen, and evolved the abundances for ionized, atomic, and molecular states using a modified version of Recfast++ (which we have released publicly at https://github.com/jamesgurian/RecfastJulia). We also provide an analytical approximation for the final abundances. We then calculate the effects of the atomic dark matter on the linear power spectrum, which enter through a dark-photon diffusion and dark acoustic oscillations. At the formation time, the atomic dark matter model suppresses halo abundances on scales smaller than the diffusion scale, just like the warm dark matter models suppress the abundance below the free-streaming scale. The subsequent evolution with radiative cooling, however, will alter the halo mass function further.

preprint2022arXiv

Molecular Chemistry for Dark Matter III: DarkKROME

Dark matter that is dissipative may cool sufficiently to form compact objects, including black holes. Determining the abundance and mass spectrum of those objects requires an accurate model of the chemistry relevant for the cooling of the dark matter gas. Here we introduce a chemistry tool for dark matter, DarkKROME, an extension of the KROME software package. DarkKROME is designed to include all atomic and molecular processes relevant for dark matter with two unequal-mass fundamental fermions, interacting via a massless-photon mediated $U(1)$ force. We use DarkKROME to perform one-zone collapse simulations and study the evolution of temperature-density phase diagrams for various dark-sector parameters.

preprint2022arXiv

Surface Brightness Profile of Lyman-$α$ Halos out to 320 kpc in HETDEX

We present the median-stacked Lyman-$α$ surface brightness profile of 968 spectroscopically selected Lyman-$α$ emitting galaxies (LAEs) at redshifts $1.9<z<3.5$ in the early data of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The selected LAEs are high-confidence Lyman-$α$ detections with large signal-to-noise ratios observed with good seeing conditions (point-spread-function full-width-at-half-maximum $<1.4"$), excluding active galactic nuclei (AGN). The Lyman-$α$ luminosities of the LAEs are $10^{42.4}-10^{43}\, \mathrm{erg}\, \mathrm{s}^{-1}$. We detect faint emission in the median-stacked radial profiles at the level of $(3.6\pm 1.3)\times 10^{-20}\,\mathrm{erg}\,\mathrm{s}^{-1}\,\mathrm{cm}^{-2}\,\mathrm{arcsec}^{-2}$ from the surrounding Lyman-$α$ halos out to $r\simeq 160$ kpc (physical). The shape of the median-stacked radial profile is consistent at $r<80\,\mathrm{kpc}$ with that of much fainter LAEs at $3<z<4$ observed with the Multi Unit Spectroscopic Explorer (MUSE), indicating that the median-stacked Lyman-$α$ profiles have similar shapes at redshifts $2<z<4$ and across a factor of $10$ in Lyman-$α$ luminosity. While we agree with the results from the MUSE sample at $r<80\,\mathrm{kpc}$, we extend the profile over a factor of two in radius. At $r>80\,\mathrm{kpc}$, our profile is flatter than the MUSE model. The measured profile agrees at most radii with that of galaxies in the Byrohl et al. (2021) cosmological radiative transfer simulation at $z=3$. This suggests that the surface brightness of a Lyman-$α$ halo at $r\lesssim 100$ kpc is dominated by resonant scattering of Lyman-$α$ photons from star-forming regions in the central galaxy, whereas at $r > 100$ kpc it is dominated by photons from galaxies in surrounding dark matter halos.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) I. Sample selection

We present the first Active Galactic Nuclei (AGN) catalog in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) observed between January 2017 and June 2020. HETDEX is an ongoing spectroscopic survey with no pre-selection based on magnitudes, colors or morphologies, enabling us to select AGN based on their spectral features. Both luminous quasars and low-luminosity Seyferts are found in our catalog. AGN candidates are selected with at least two significant AGN emission lines, such as the LyA and CIV line pair, or with single broad emission lines (FWHM > 1000 km/s). Each source is further confirmed by visual inspections. This catalog contains 5,322 AGN, covering an effective sky coverage of 30.61 deg^2. A total of 3,733 of these AGN have secure redshifts, and we provide redshift estimates for the remaining 1,589 single broad-line AGN with no cross matched spectral redshifts from SDSS DR14Q. The redshift range of the AGN catalog is 0.25 < z < 4.32, with a median of z = 2.1. The bolometric luminosity range is 10^9-10^14 Lsun with a median of 10^12 Lsun. The median r-band magnitude of the AGN is 21.6 mag, with 34% of the AGN have r > 22.5, and 2.6% reaching the detection limit at r ~ 26 mag of the deepest imaging surveys we searched. We also provide a composite spectrum of the AGN sample covering 700 AA - 4400 AA.

preprint2022arXiv

The Active Galactic Nuclei in the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) II. Luminosity Function

We present the LyA emission line luminosity function (LF) of the Active Galactic Nuclei (AGN) in the first release of the Hobby-Eberly Telescope Dark Energy Experiment Survey (HETDEX) AGN catalog (Liu et al. 2022, Paper I). The AGN are selected either by emission-line pairs characteristic of AGN or by single broad emission line, free of any photometric pre-selections (magnitude/color/morphology). The sample consists of 2,346 AGN spanning 1.88<z<3.53, covering an effective area of 30.61 deg^2. Approximately 2.6 of the HETDEX AGN are not detected at $>5σ$ confidence at r~26 in the deepest $r$-band images we have searched. The LyA line luminosity ranges from ~10^42.3 to ~10^45.9 erg s^-1. Our LyA LF shows a turnover luminosity with opposite slopes on the bright end and the faint end: The space density is highest at L_LyA^*=10^43.4 erg s^-1. We explore the evolution of the AGN LF over a broader redshift range (0.8<z<3); constructing the rest-frame ultraviolet (UV) LF with the 1450 AA monochromatic luminosity of the power-law component of the continuum ($\rm M_{1450}$) from M_1450~-18 to ~-27.5. We divide the sample into three redshift bins (z~1.5, 2.1, and 2.6). In all three redshift bins, our UV LFs indicate that the space density of AGN is highest at the turnover luminosity M_1450^* with opposite slopes on the bright end and the faint end. The M_1450 LFs in the three redshift bins can be well-fit with a luminosity-evolution-density-evolution (LEDE) model: the turnover luminosity (M_1450^*) increases and the turnover density (Phi^*) decreases with increasing redshift.

preprint2021arXiv

Covariance of the matter power spectrum including the survey window function effect: N-body simulations vs. fifth-order perturbation theory on grid

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.

preprint2021arXiv

First HETDEX Spectroscopic Determinations of Ly$α$ and UV Luminosity Functions at $z=2-3$: Bridging a Gap Between Faint AGN and Bright Galaxies

We present Ly$α$ and ultraviolet-continuum (UV) luminosity functions (LFs) of galaxies and active galactic nuclei (AGN) at $z=2.0-3.5$ determined by the un-targetted optical spectroscopic survey of the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). We combine deep Subaru imaging with HETDEX spectra resulting in $11.4$ deg$^2$ of fiber-spectra sky coverage, obtaining $18320$ galaxies spectroscopically identified with Ly$α$ emission, $2126$ of which host type 1 AGN showing broad (FWHM$~>1000$ km s$^{-1}$) Ly$α$ emission lines. We derive the Ly$α$ (UV) LF over 2 orders of magnitude covering bright galaxies and AGN in $\log L_\mathrm{Lyα}/\mathrm{[erg~s^{-1}]}=43.3-45.5$ ($-27<M_\mathrm{UV}<-20$) by the $1/V_\mathrm{max}$ estimator. Our results reveal the bright-end hump of the Ly$α$ LF is composed of type 1 AGN. In conjunction with previous spectroscopic results at the faint end, we measure a slope of the best-fit Schechter function to be $α_\mathrm{Sch}=-1.70^{+0.13}_{-0.14}$, which indicates $α_\mathrm{Sch}$ steepens from $z=2-3$ towards high redshift. Our UV LF agrees well with previous AGN UV LFs, and extends to faint-AGN and bright-galaxy regimes. The number fraction of Ly$α$-emitting objects ($X_\mathrm{LAE}$) increases from $M_\mathrm{UV}^*\sim-21$ to bright magnitude due to the contribution of type 1 AGN, while previous studies claim that $X_\mathrm{Lyα}$ decreases from faint magnitude to $M_\mathrm{UV}^*$, suggesting a valley in the $X_\mathrm{Lyα}-$magnitude relation at $M_\mathrm{UV}^*$. Comparing our UV LF of type 1 AGN at $z=2-3$ with those at $z=0$, we find that the number density of faint ($M_\mathrm{UV}>-21$) type 1 AGN increases from $z\sim2$ to $z\sim0$ as opposed to the evolution of bright ($M_\mathrm{UV}<-21$) type 1 AGN, suggesting the AGN downsizing in the rest-frame UV luminosity.

preprint2021arXiv

Grid-based calculations of redshift-space matter fluctuations from perturbation theory: UV sensitivity and convergence at the field level

Perturbation theory (PT) has been used to interpret the observed nonlinear large-scale structure statistics at the quasi-linear regime. To facilitate the PT-based analysis, we have presented the GridSPT algorithm, a grid-based method to compute the nonlinear density and velocity fields in standard perturbation theory (SPT) from a given linear power spectrum. Here, we further put forward the approach by taking the redshift-space distortions into account. With the new implementation, we have, for the first time, generated the redshift-space density field to the fifth order and computed the next-to-next-to-leading order (2 loop) power spectrum and the next-to-leading order (1 loop) bispectrum of matter clustering in redshift space. By comparing the result with corresponding analytical SPT calculation and $N$-body simulations, we find that the SPT calculation (A) suffers much more from the UV sensitivity due to the higher-derivative operators and (B) deviates from the $N$-body results from the Fourier wavenumber smaller than real space $k_{\rm max}$. Finally, we have shown that while Padé approximation removes spurious features in morphology, it does not improve the modeling of power spectrum and bispectrum.

preprint2021arXiv

The Hobby-Eberly Telescope Dark Energy Experiment (HETDEX) Survey Design, Reductions, and Detections

We describe the survey design, calibration, commissioning, and emission-line detection algorithms for the Hobby-Eberly Telescope Dark Energy Experiment (HETDEX). The goal of HETDEX is to measure the redshifts of over a million Ly$α$ emitting galaxies between 1.88<z<3.52, in a 540 deg^2 area encompassing a co-moving volume of 10.9 Gpc^3. No pre-selection of targets is involved; instead the HETDEX measurements are accomplished via a spectroscopic survey using a suite of wide-field integral field units distributed over the focal plane of the telescope. This survey measures the Hubble expansion parameter and angular diameter distance, with a final expected accuracy of better than 1%. We detail the project's observational strategy, reduction pipeline, source detection, and catalog generation, and present initial results for science verification in the COSMOS, Extended Groth Strip, and GOODS-N fields. We demonstrate that our data reach the required specifications in throughput, astrometric accuracy, flux limit, and object detection, with the end products being a catalog of emission-line sources, their object classifications, and flux-calibrated spectra.

preprint2020arXiv

Alleviating the Tension in the Cosmic Microwave Background using Planck-Scale Physics

Certain anomalies in the CMB bring out a tension between the six-parameter flat $Λ$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.

preprint2020arXiv

Fast Calculation of Nonlinear Redshift-space Galaxy Power Spectrum Including Selection Bias

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

preprint2019arXiv

Constraining Cosmology with Big Data Statistics of Cosmological Graphs

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 $Λ$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) $α$, (2) closed-to-connected triple fraction (the so-called transitivity) $τ_Δ$, 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.

preprint2019arXiv

Gravitational waves, CMB polarization, and the Hubble tension

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.

preprint2019arXiv

The Odd-Parity Galaxy Bispectrum

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.

preprint2016arXiv

Constraining primordial and gravitational mode coupling with the position-dependent bispectrum of the large-scale structure

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

preprint2016arXiv

Cosmological nonlinear density and velocity power spectra including nonlinear vector and tensor modes

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.

preprint2016arXiv

Doppler term in the galaxy two-point correlation function: wide-angle, velocity, Doppler lensing and cosmic acceleration effects

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.

preprint2015arXiv

Inflationary tensor fossils in large-scale structure

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.

preprint2015arXiv

Toward Accurate Modeling of the Nonlinear Matter Bispectrum: Standard Perturbation Theory and Transients from Initial Conditions

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\'ı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.

preprint2014arXiv

A halo bias function measured deeply into voids without stochasticity

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

preprint2014arXiv

Clustering-based redshift estimation: method and application to data

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.

preprint2014arXiv

Large-Scale Structure Observables in General Relativity

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.

preprint2014arXiv

Linking the BICEP2 result and the hemispherical power asymmetry through spatial variation of $r$

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.

preprint2014arXiv

Non-linear matter bispectrum in general relativity

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.

preprint2014arXiv

Non-linear power spectra in the synchronous gauge

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.

preprint2014arXiv

Silk damping at a redshift of a billion: a new limit on small-scale adiabatic perturbations

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.

preprint2014arXiv

Tensor-induced B modes with no temperature fluctuations

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.

preprint2014arXiv

Testing Inflation with Large Scale Structure: Connecting Hopes with Reality

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

preprint2014arXiv

The redshift-space galaxy two-point correlation function and baryon acoustic oscillations

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ξ(r,μ)$ 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.

preprint2013arXiv

Anisotropic imprint of long-wavelength tensor perturbations on cosmic structure

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.

preprint2013arXiv

Cosmic Clocks

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.

preprint2013arXiv

Cosmic Rulers

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

preprint2013arXiv

Galaxy redshift surveys with sparse sampling

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.

preprint2013arXiv

Large-Scale Structure with Gravitational Waves II: Shear

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.

preprint2013arXiv

Peak-Background Split, Renormalization, and Galaxy Clustering

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 δ_L coarse-grained on a scale R_L. In the latter case, the predicted tracer correlation function receives contributions of order <δ_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.

preprint2013arXiv

Seeking Inflation Fossils in the Cosmic Microwave Background

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.

preprint2013arXiv

Teasing bits of information out of the CMB energy spectrum

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.

preprint2013arXiv

The effect of aberration on partial-sky measurements of the cosmic microwave background temperature power spectrum

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.

preprint2013arXiv

The Effects of Local Primordial Non-Gaussianity on the Formation and Evolution of Galaxies

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.

preprint2013arXiv

The Pesky Power Asymmetry

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.

preprint2013arXiv

Wigner-Eckart theorem in cosmology: Bispectra for total-angular-momentum waves

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.

preprint2012arXiv

Clustering Fossils from the Early Universe

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.

preprint2012arXiv

Halo statistics in non-Gaussian cosmologies: the collapsed fraction, conditional mass function, and halo bias from the path-integral excursion set method

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]

preprint2012arXiv

Large-scale clustering of galaxies in general relativity

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.

preprint2012arXiv

Large-Scale Structure with Gravitational Waves I: Galaxy Clustering

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.

preprint2012arXiv

Total Angular Momentum Waves for Scalar, Vector, and Tensor Fields

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.

preprint2011arXiv

Accurate Predictions for the Scale-Dependent Galaxy Bias from Primordial Non-Gaussianity

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.

preprint2011arXiv

Non-Gaussian Halo Bias Re-examined: Mass-dependent Amplitude from the Peak-Background Split and Thresholding

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.

preprint2011arXiv

Primordial Non-Gaussianity and the Statistics of Weak Lensing and other Projected Density Fields

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.

preprint2011arXiv

The HETDEX Pilot Survey. II. The Evolution of the Ly-alpha Escape Fraction from the UV Slope and Luminosity Function of 1.9<z<3.8 LAEs

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]

preprint2010arXiv

General relativistic effects on non-linear power spectra

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.

preprint2010arXiv

HETDEX pilot survey for emission-line galaxies - I. Survey design, performance, and catalog

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.

preprint2009arXiv

Galaxy-CMB and galaxy-galaxy lensing on large scales: sensitivity to primordial non-Gaussianity

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.

preprint2008arXiv

Extracting Angular Diameter Distance and Expansion Rate of the Universe from Two-dimensional Galaxy Power Spectrum at High Redshifts: Baryon Acoustic Oscillation Fitting versus Full Modeling

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.

preprint2006arXiv

Perturbation Theory Reloaded: Analytical Calculation of Non-linearity in Baryonic Oscillations in the Real Space Matter Power Spectrum

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.