Source author record

Kendrick M. Smith

Kendrick M. Smith 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

35works
7topics
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

35 published item(s)

preprint2026arXiv

Constraining $f_{NL}$ using the Large-Scale Modulation of Small-Scale Statistics

We implement a novel formalism to constrain primordial non-Gaussianity of the local type from the large-scale modulation of the small-scale power spectrum. Our approach combines information about primordial non-Gaussianity contained in the squeezed bispectrum and the collapsed trispectrum of large-scale structure together in a computationally amenable and consistent way, while avoiding the need to model complicated covariances of higher $N$-point functions. This work generalizes our recent work, which used a neural network estimate of local power, to the more conventional local power spectrum statistics, and explores using both matter field and halo catalogues from the Quijote simulations. We find that higher $N$-point functions of the matter field can provide strong constraints on $f_{NL}$, but higher $N$-point functions of the halo field, at the halo density of Quijote, only marginally improve constraints from the two-point function.

preprint2022arXiv

A High-Time Resolution Search for Compact Objects using Fast Radio Burst Gravitational Lens Interferometry with CHIME/FRB

The gravitational field of compact objects, such as primordial black holes, can create multiple images of background sources. For transients such as fast radio bursts (FRBs), these multiple images can be resolved in the time domain. Under certain circumstances, these images not only have similar burst morphologies but are also phase-coherent at the electric field level. With a novel dechannelization algorithm and a matched filtering technique, we search for repeated copies of the same electric field waveform in observations of FRBs detected by the FRB backend of the Canadian Hydrogen Mapping Intensity Experiment (CHIME). An interference fringe from a coherent gravitational lensing signal will appear in the time-lag domain as a statistically-significant peak in the time-lag autocorrelation function. We calibrate our statistical significance using telescope data containing no FRB signal. Our dataset consists of $\sim$100-ms long recordings of voltage data from 172 FRB events, dechannelized to 1.25-ns time resolution. This coherent search algorithm allows us to search for gravitational lensing signatures from compact objects in the mass range of $10^{-4}-10^{4} ~\mathrm{M_{\odot}}$. After ruling out an anomalous candidate due to diffractive scintillation, we find no significant detections of gravitational lensing in the 172 FRB events that have been analyzed. In a companion work [Leung, Kader+2022], we interpret the constraints on dark matter from this search.

preprint2022arXiv

Constraining Primordial Black Holes using Fast Radio Burst Gravitational-Lens Interferometry with CHIME/FRB

Fast radio bursts (FRBs) represent an exciting frontier in the study of gravitational lensing, due to their brightness, extragalactic nature, and the compact, coherent characteristics of their emission. In a companion work [Kader, Leung+2022], we use a novel interferometric method to search for gravitationally lensed FRBs in the time domain using bursts detected by CHIME/FRB. There, we dechannelize and autocorrelate electric field data at a time resolution of 1.25 ns. This enables a search for FRBs whose emission is coherently deflected by gravitational lensing around a foreground compact object such as a primordial black hole (PBH). Here, we use our non-detection of lensed FRBs to place novel constraints on the PBH abundance outside the Local Group. We use a novel two-screen model to take into account decoherence from scattering screens in our constraints. Our constraints are subject to a single astrophysical model parameter -- the effective distance between an FRB source and the scattering screen, for which we adopt a fiducial distance of 1 parsec. We find that coherent FRB lensing is a sensitive probe of sub-solar mass compact objects. Having observed no lenses in $172$ bursts from $114$ independent sightlines through the cosmic web, we constrain the fraction of dark matter made of compact objects, such as PBHs, to be $f \lesssim 0.8$, if their masses are $\sim 10^{-3} M_{\odot}$.

preprint2022arXiv

Localizing FRBs through VLBI with the Algonquin Radio Observatory 10-m Telescope

The CHIME/FRB experiment has detected thousands of Fast Radio Bursts (FRBs) due to its sensitivity and wide field of view; however, its low angular resolution prevents it from localizing events to their host galaxies. Very Long Baseline Interferometry (VLBI), triggered by FRB detections from CHIME/FRB will solve the challenge of localization for non-repeating events. Using a refurbished 10-m radio dish at the Algonquin Radio Observatory located in Ontario Canada, we developed a testbed for a VLBI experiment with a theoretical ~<30 masec precision. We provide an overview of the 10-m system and describe its refurbishment, the data acquisition, and a procedure for fringe fitting that simultaneously estimates the geometric delay used for localization and the dispersive delay from the ionosphere. Using single pulses from the Crab pulsar, we validate the system and localization procedure, and analyze the clock stability between sites, which is critical for phase-referencing an FRB event. We find a localization of 50 masec is possible with the performance of the current system. Furthermore, for sources with insufficient signal or restricted wideband to simultaneously measure both geometric and ionospheric delays, we show that the differential ionospheric contribution between the two sites must be measured to a precision of 1e-8 pc/cc to provide a reasonable localization from a detection in the 400--800 MHz band. Finally we show detection of an FRB observed simultaneously in the CHIME and the Algonquin 10-m telescope, the first FRB cross-correlated in this very long baseline. This project serves as a testbed for the forthcoming CHIME/FRB Outriggers project.

preprint2022arXiv

The Hydrogen Intensity and Real-time Analysis eXperiment: 256-Element Array Status and Overview

The Hydrogen Intensity and Real-time Analysis eXperiment (HIRAX) is a radio interferometer array currently in development, with an initial 256-element array to be deployed at the South African Radio Astronomy Observatory (SARAO) Square Kilometer Array (SKA) site in South Africa. Each of the 6m, $f/0.23$ dishes will be instrumented with dual-polarisation feeds operating over a frequency range of 400-800 MHz. Through intensity mapping of the 21 cm emission line of neutral hydrogen, HIRAX will provide a cosmological survey of the distribution of large-scale structure over the redshift range of $0.775 < z < 2.55$ over $\sim$15,000 square degrees of the southern sky. The statistical power of such a survey is sufficient to produce $\sim$7 percent constraints on the dark energy equation of state parameter when combined with measurements from the Planck satellite. Additionally, HIRAX will provide a highly competitive platform for radio transient and HI absorber science while enabling a multitude of cross-correlation studies. In this paper, we describe the science goals of the experiment, overview of the design and status of the sub-components of the telescope system, and describe the expected performance of the initial 256-element array as well as the planned future expansion to the final, 1024-element array.

preprint2020arXiv

Accurate Analytic Model for the Weak Lensing Convergence One-Point Probability Distribution Function and its Auto-Covariance

The one-point probability distribution function (PDF) is a powerful summary statistic for non-Gaussian cosmological fields, such as the weak lensing (WL) convergence reconstructed from galaxy shapes or cosmic microwave background (CMB) maps. Thus far, no analytic model has been developed that successfully describes the high-convergence tail of the WL convergence PDF for small smoothing scales from first principles. Here, we present a halo-model formalism to compute the WL convergence PDF, building upon our previous results for the thermal Sunyaev-Zel'dovich field. Furthermore, we extend our formalism to analytically compute the covariance matrix of the convergence PDF. Comparisons to numerical simulations generally confirm the validity of our formalism in the non-Gaussian, positive tail of the WL convergence PDF, but also reveal the convergence PDF's strong sensitivity to small-scale systematic effects in the simulations (e.g., due to finite resolution). Finally, we present a simple Fisher forecast for a Rubin-Observatory-like survey, based on our new analytic model. Considering the $\{A_s, Ω_m, Σm_ν\}$ parameter space and assuming a Planck CMB prior on $A_s$ only, we forecast a marginalized constraint $σ(Σm_ν) \approx 0.08$ eV from the WL convergence PDF alone, even after marginalizing over parameters describing the halo concentration-mass relation. This error bar on the neutrino mass sum is comparable to the minimum value allowed in the normal hierarchy, illustrating the strong constraining power of the WL convergence PDF. We make our code publicly available at https://github.com/leanderthiele/hmpdf.

preprint2020arXiv

Exploring KSZ velocity reconstruction with $N$-body simulations and the halo model

KSZ velocity reconstruction is a recently proposed method for mapping the largest-scale modes of the universe, by applying a quadratic estimator $\hat{v}_r$ to the small-scale CMB and a galaxy catalog. We implement kSZ velocity reconstruction in an $N$-body simulation pipeline and explore its properties. We find that the reconstruction noise can be larger than the analytic prediction which is usually assumed. We revisit the analytic prediction and find additional noise terms which explain the discrepancy. The new terms are obtained from a six-point halo model calculation, and are analogous to the $N^{(1)}$ and $N^{(3/2)}$ biases in CMB lensing. We implement an MCMC pipeline which estimates $f_{NL}$ from $N$-body kSZ simulations, and show that it recovers unbiased estimates of $f_{NL}$, with statistical errors consistent with a Fisher matrix forecast. Overall, these results confirm that kSZ velocity reconstruction will be a powerful probe of cosmology in the near future, but new terms should be included in the noise power spectrum.

preprint2019arXiv

Characterizing fast radio bursts through statistical cross-correlations

Understanding the origin of fast radio bursts (FRB's) is a central unsolved problem in astrophysics that is severely hampered by their poorly determined distance scale. Determining the redshift distribution of FRB's appears to require arcsecond angular resolution, in order to associate FRB's with host galaxies. In this paper, we forecast prospects for determining the redshift distribution without host galaxy associations, by cross-correlating FRB's with a galaxy catalog such as the SDSS photometric sample. The forecasts are extremely promising: a survey such as CHIME/FRB that measures catalogs of $\sim 10^3$ FRB's with few-arcminute angular resolution can place strong constraints on the FRB redshift distribution, by measuring the cross-correlation as a function of galaxy redshift $z$ and FRB dispersion measure $D$. In addition, propagation effects from free electron inhomogeneities modulate the observed FRB number density, either by shifting FRB's between dispersion measure (DM) bins or through DM-dependent selection effects. We show that these propagation effects, coupled with the spatial clustering between galaxies and free electrons, can produce FRB-galaxy correlations which are comparable to the intrinsic clustering signal. Such effects can be disentangled based on their angular and $(z, D)$ dependence, providing an opportunity to study not only FRB's but the clustering of free electrons.

preprint2016arXiv

New algorithms for radio pulsar search

The computational cost of searching for new pulsars is a limiting factor for upcoming radio telescopes such as SKA. We introduce four new algorithms: an optimal constant-period search, a coherent tree search which permits optimal searching with O(1) cost per model, a semicoherent search which combines information from coherent subsearches while preserving as much phase information as possible, and a hierarchical search which interpolates between the coherent and semicoherent limits. Taken together, these algorithms improve the computational cost of pulsar search by several orders of magnitude. In this paper, we consider the simple case of a constant-acceleration phase model, but our methods should generalize to more complex search spaces.

preprint2015arXiv

Evidence for the kinematic Sunyaev-Zeľdovich effect with ACTPol and velocity reconstruction from BOSS

We use microwave temperature maps from two seasons of data from the Atacama Cosmology Telescope (ACTPol) at 146 GHz, together with the Constant Mass CMASS galaxy sample from the Baryon Oscillation Spectroscopic Survey to measure the kinematic Sunyaev-Zeľdovich (kSZ) effect over the redshift range z = 0.4 - 0.7. We use galaxy positions and the continuity equation to obtain a reconstruction of the line-of-sight velocity field. We stack the cosmic microwave background temperature at the location of each halo, weighted by the corresponding reconstructed velocity. The resulting best fit kSZ model is preferred over the no-kSZ hypothesis at 3.3sigma and 2.9sigma for two independent velocity reconstruction methods, using 25,537 galaxies over 660 square degrees. The effect of foregrounds that are uncorrelated with the galaxy velocities is expected to be well below our signal, and residual thermal Sunyaev-Zeľdovich contamination is controlled by masking the most massive clusters. Finally, we discuss the systematics involved in converting our measurement of the kSZ amplitude into the mean free electron fraction of the halos in our sample.

preprint2015arXiv

On testing and extending the inflationary consistency relation for tensor modes

If observations confirm BICEP2's claim of a tensor-scalar ratio $r\approx 0.2$ on CMB scales, then the inflationary consistency relation $n_{t}=-r/8$ predicts a small negative value for the tensor spectral index $n_t$. We show that future CMB polarization experiments should be able to confirm this prediction at several sigma. We also show how to properly extend the consistency relation to solar system scales, where the primordial gravitational wave density $Ω_{gw}$ could be measured by proposed experiments such as the Big Bang Observer. This would provide a far more stringent test of the consistency relation and access much more detailed information about the early universe.

preprint2015arXiv

Optimal analysis of the CMB trispectrum

We develop a general framework for data analysis and phenomenology of the CMB four-point function or trispectrum. To lowest order in the derivative expansion, the inflationary action admits three quartic operators consistent with symmetry: $\dotσ^4$, $\dotσ^2 (\partialσ^2)$, and $(\partialσ)^4$. In single field inflation, only the first of these operators can be the leading non-Gaussian signal. A Fisher matrix analysis shows that there is one near-degeneracy among the three CMB trispectra, so we parameterize the trispectrum with two coefficients $g_{NL}^{\dotσ^4}$ and $g_{NL}^{(\partialσ)^4}$, in addition to the coefficient $g_{NL}^{\rm loc}$ of $ζ^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{\rm loc}, f_{NL}^{\rm equil}, f_{NL}^{\rm orth})$ commonly used to parameterize the CMB three-point function. We next turn to data analysis and show how to represent these trispectra in a factorizable form which leads to computationally fast operations such as evaluating a CMB estimator or simulating a non-Gaussian CMB. We discuss practical issues in CMB analysis pipelines, and perform an optimal analysis of WMAP data. Our minimum-variance estimates are $g_{NL}^{\rm loc} = (-3.80 \pm 2.19) \times 10^5$, $g_{NL}^{\dotσ^4} = (-3.20 \pm 3.09) \times 10^6$, and $g_{NL}^{(\partialσ)^4} = (-10.8 \pm 6.33) \times 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.

preprint2014arXiv

New natural shapes of non-Gaussianity from high-derivative interactions and their optimal limits from WMAP 9-year data

Given the fantastic experimental effort, it is important to thoroughly explore the signature space of inflationary models. The fact that higher derivative operators do not renormalize lower derivative ones allows us to find a large class of technically natural single-clock inflationary models where, in the context of the Effective Field Theory of Inflation, the leading interactions have many derivatives. We systematically explore the 3-point function induced by these models and their overlap with the standard equilateral and orthogonal templates. We find that in order to satisfactorily cover the signature space of these models, two new additional templates need to be included. We then perform the optimal analysis of the WMAP 9-year data for the resulting four templates, finding that the overall significance of a non-zero signal is between 2--2.5$σ$, depending on the choice of parameter space, partially driven by the preference for nonzero $f_{NL}^{\rm orth}$ in WMAP9.

preprint2014arXiv

On quantifying and resolving the BICEP2/Planck tension over gravitational waves

The recent BICEP2 measurement of primordial gravity waves (r = 0.2^{+0.07}_{-0.05}) appears to be in tension with the upper limit from WMAP (r<0.13 at 95% CL) and Planck (r<0.11 at 95% CL). We carefully quantify the level of tension and show that it is very significant (around 0.1% unlikely) when the observed deficit of large-scale temperature power is taken into account. We show that measurements of TE and EE power spectra in the near future will discriminate between the hypotheses that this tension is either a statistical fluke, or a sign of new physics. We also discuss extensions of the standard cosmological model that relieve the tension, and some novel ways to constrain them.

preprint2014arXiv

The Atacama Cosmology Telescope: A Measurement of the Thermal Sunyaev-Zel'dovich One-Point PDF

We present a measurement of the one-point probability distribution function (PDF) of the thermal Sunyaev-Zel'dovich (tSZ) decrement in the pixel temperature histogram of filtered 148 GHz sky maps from the Atacama Cosmology Telescope (ACT). The PDF includes the signal from all galaxy clusters in the map, including objects below the signal-to-noise threshold for individual detection, making it a particularly sensitive probe of the amplitude of matter density perturbations, $σ_8$. We use a combination of analytic halo model calculations and numerical simulations to compute the theoretical tSZ PDF and its covariance matrix, accounting for all noise sources and including relativistic corrections. From the measured ACT 148 GHz PDF alone, we find $σ_8 = 0.793 \pm 0.018$, with additional systematic errors of $\pm 0.017$ due to uncertainty in intracluster medium gas physics and $\pm 0.006$ due to uncertainty in infrared point source contamination. Using effectively the same data set, the statistical error here is a factor of two lower than that found in ACT's previous $σ_8$ determination based solely on the skewness of the tSZ signal. In future temperature maps with higher sensitivity, the tSZ PDF will break the degeneracy between intracluster medium gas physics and cosmological parameters.

preprint2014arXiv

Using Large Scale Structure to test Multifield Inflation

Primordial non-Gaussianity of local type is known to produce a scale-dependent contribution to the galaxy bias. Several classes of multi-field inflationary models predict non-Gaussian bias which is stochastic, in the sense that dark matter and halos don't trace each other perfectly on large scales. In this work, we forecast the ability of next-generation Large Scale Structure surveys to constrain common types of primordial non-Gaussianity like $f_{NL}$, $g_{NL}$ and $τ_{NL}$ using halo bias, including stochastic contributions. We provide fitting functions for statistical errors on these parameters which can be used for rapid forecasting or survey optimization. A next-generation survey with volume $V = 25 h^{-3}$Gpc$^3$, median redshift $z = 0.7$ and mean bias $b_g = 2.5$, can achieve $σ(f_{NL}) = 6$, $σ(g_{NL}) = 10^5$ and $σ(τ_{NL}) = 10^3$ if no mass information is available. If halo masses are available, we show that optimally weighting the halo field in order to reduce sample variance can achieve $σ(f_{NL}) = 1.5$, $σ(g_{NL}) = 10^4$ and $σ(τ_{NL}) = 100$ if halos with mass down to $M_{min} = 10^{11}$ $h^{-1} M_\odot $ are resolved, outperforming Planck by a factor of 4 on $f_{NL}$ and nearly an order of magnitude on $g_{NL}$ and $τ_{NL}$. Finally, we study the effect of photometric redshift errors and discuss degeneracies between different non-Gaussian parameters, as well as the impact of marginalizing Gaussian bias and shot noise.

preprint2013arXiv

Non-Gaussian structure of the lensed CMB power spectra covariance matrix

Gravitational lensing of the Cosmic Microwave Background (CMB) encodes cosmological information in the observed anisotropies of temperature and polarization. Accurate extraction of this additional information requires precise modeling of the covariance matrix of the power spectra of observed CMB fields. We introduce a new analytical model to describe the non-Gaussian structure of this covariance matrix and display the importance of second-order terms that were previously neglected. When compared with direct numerical simulations our model captures parameter errors to better than a few percent for cases where the non-Gaussianity causes an order unity degradation in errors. We also provide a detailed comparison between the information content of lensed CMB power spectra and ideal reconstruction of the lensing potential. We illustrate the impact of the non-Gaussian terms in the power spectrum covariance by providing Fisher errors on the sum of the masses of the neutrinos, the dark energy equation of state, and the curvature of the Universe.

preprint2012arXiv

On the Equivalence of Barrier Crossing, Peak-Background Split, and Local Biasing

Several, apparently distinct, formalisms exist in the literature for predicting the clustering of dark matter halos. It has been noticed on a case-by-case basis that the predictions of these different methods agree in specific examples, but there is no general proof that they are equivalent. In this paper, we give a simple proof of the mathematical equivalence of barrier crossing, peak-background split, and local biasing.

preprint2012arXiv

Stochastic Bias from Non-Gaussian Initial Conditions

In this article, we show that a stochastic form of scale-dependent halo bias arises in multi-source inflationary models, where multiple fields determine the initial curvature perturbation. We derive this effect for general non-Gaussian initial conditions and study various examples, such as curvaton models and quasi-single field inflation. We present a general formula for both the stochastic and the non-stochastic parts of the halo bias, in terms of the N-point cumulants of the curvature perturbation at the end of inflation. At lowest order, the stochasticity arises if the collapsed limit of the four-point function is boosted relative to the square of the three-point function in the squeezed limit. We derive all our results in two ways, using the barrier crossing formalism and the peak-background split method. In a companion paper, we prove that these two approaches are mathematically equivalent.

preprint2011arXiv

A universal bound on N-point correlations from inflation

Models of inflation in which non-Gaussianity is generated outside the horizon, such as curvaton models, generate distinctive higher-order correlation functions in the CMB and other cosmological observables. Testing for violation of the Suyama-Yamaguchi inequality tauNL >= (6/5 fNL)^2, where fNL and tauNL denote the amplitude of the three-point and four-point functions in certain limits, has been proposed as a way to distinguish qualitative classes of models. This inequality has been proved for a wide range of models, but only weaker versions have been proved in general. In this paper, we give a proof that the Suyama-Yamaguchi inequality is always satisfied. We discuss scenarios in which the inequality may appear to be violated in an experiment such as Planck, and how this apparent violation should be interpreted. We analyze a specific example, the "ungaussiton" model, in which leading-order scaling relations suggest that the Suyama-Yamaguchi inequality is eventually violated, and show that the inequality always holds.

preprint2011arXiv

Algorithms for bispectra: forecasting, optimal analysis, and simulation

We propose a factorizability ansatz for angular bispectra which permits fast algorithms for forecasting, analysis, and simulation, yet is general enough to encompass many interesting CMB bispectra. We describe a suite of general algorithms which apply to any bispectrum which can be represented in factorizable form. First, we present algorithms for Fisher matrix forecasts and the related problem of optimizing the factorizable representation, giving a Fisher forecast for Planck as an example. We show that the CMB can give independent constraints on the amplitude of primordial bispectra of both local and equilateral shape as well as those created by secondary anisotropies. We also show that the ISW-lensing bispectrum should be detected by Planck and could bias estimates of the local type of non-Gaussianity if not properly accounted for. Second, we implement a bispectrum estimator which is fully optimal in the presence of sky cuts and inhomogeneous noise, extends the generality of fast estimators which have been limited to a few specific forms of the bispectrum, and improves the running time of existing implementations by several orders of magnitude. Third, we give an algorithm for simulating random, weakly non-Gaussian maps with prescribed power spectrum and factorizable bispectrum.

preprint2011arXiv

Beyond the power spectrum: primordial and secondary non-Gaussianity in the microwave background

Cosmic microwave background observations are most commonly analyzed by estimating the power spectrum. In the limit where the CMB statistics are perfectly Gaussian, this extracts all the information, but the CMB also contains detectable non-Gaussian contributions from secondary, and possibly primordial, sources. We review possible sources of CMB non-Gaussianity and describe statistical techniques which are optimized for measuring them, complementing the power spectrum analysis. The machinery of $N$-point correlation functions provides a unifying framework for optimal estimation of primordial non-Gaussian signals or gravitational lensing. We review recent results from applying these estimators to data from the WMAP satellite mission.

preprint2011arXiv

Halo clustering and g_{NL}-type primordial non-Gaussianity

A wide range of multifield inflationary models generate non-Gaussian initial conditions in which the initial adiabatic fluctuation is of the form (zeta_G + g_{NL} zeta_G^3). We study halo clustering in these models using two different analytic methods: the peak-background split framework, and brute force calculation in a barrier crossing model, obtaining agreement between the two. We find a simple, theoretically motivated expression for halo bias which agrees with N-body simulations and can be used to constrain g_{NL} from observations. We discuss practical caveats to constraining g_{NL} using only observable properties of a tracer population, and argue that constraints obtained from populations whose observed bias is <~ 2.5 are generally not robust to uncertainties in modeling the halo occupation distribution of the population.

preprint2011arXiv

Local stochastic non-Gaussianity and N-body simulations

Large-scale clustering of highly biased tracers of large-scale structure has emerged as one of the best observational probes of primordial non-Gaussianity of the local type (i.e. f_{NL}^{local}). This type of non-Gaussianity can be generated in multifield models of inflation such as the curvaton model. Recently, Tseliakhovich, Hirata, and Slosar showed that the clustering statistics depend qualitatively on the ratio of inflaton to curvaton power ξafter reheating, a free parameter of the model. If ξis significantly different from zero, so that the inflaton makes a non-negligible contribution to the primordial adiabatic curvature, then the peak-background split ansatz predicts that the halo bias will be stochastic on large scales. In this paper, we test this prediction in N-body simulations. We find that large-scale stochasticity is generated, in qualitative agreement with the prediction, but that the level of stochasticity is overpredicted by ~30%. Other predictions, such as ξindependence of the halo bias, are confirmed by the simulations. Surprisingly, even in the Gaussian case we do not find that halo model predictions for stochasticity agree consistently with simulations, suggesting that semi-analytic modeling of stochasticity is generally more difficult than modeling halo bias.

preprint2011arXiv

No evidence for the cold spot in the NVSS radio survey

We revisit recent claims that there is a "cold spot" in both number counts and brightness of radio sources in the NVSS survey, with location coincident with the previously detected cold spot in WMAP. Such matching cold spots would be difficult if not impossible to explain in the standard LCDM cosmological model. Contrary to the claim, we find no significant evidence for the radio cold spot, after including systematic effects in NVSS, and carefully accounting for the effect of a posteriori choices when assessing statistical significance.

preprint2011arXiv

Non-Gaussianities in Single Field Inflation and their Optimal Limits from the WMAP 5-year Data

Using the recently developed effective field theory of inflation, we argue that the size and the shape of the non-Gaussianities generated by single-field inflation are generically well described by two parameters: f_NL^equil, which characterizes the size of the signal that is peaked on equilateral configurations, and f_NL^orthog, which instead characterizes the size of the signal which is peaked both on equilateral configurations and flat-triangle configurations (with opposite signs). The shape of non-Gaussianities associated with f_NL^orthog is orthogonal to the one associated to f_NL^equil, and former analysis have been mostly blind to it. We perform the optimal analysis of the WMAP 5-year data for both of these parameters. We find no evidence of non-Gaussianity, and we have the following constraints: -125 < f_NL^equil < 435, -369 < f_NL^orthog < 71 at 95% CL. We show that both of these constraints can be translated into limits on parameters of the Lagrangian of single-field inflation. For one of them, the speed of sound of the inflaton fluctuations, we find that it is either bounded to be c_s > 0.011 at 95% CL. or alternatively to be so small that the higher-derivative kinetic term dominate at horizon crossing. We are able to put similar constraints on the other operators of the inflaton Lagrangian.

preprint2011arXiv

Supersonic baryon-CDM velocities and CMB B-mode polarization

It has recently been shown that supersonic relative velocities between dark matter and baryonic matter can have a significant effect on formation of the first structures in the universe. If this effect is still non-negligible during the epoch of hydrogen reionization, it generates large-scale anisotropy in the free electron density, which gives rise to a CMB B-mode. We compute the B-mode power spectrum and find a characteristic shape with acoustic peaks at l ~ 200, 400, ... The amplitude of this signal is a free parameter which is related to the dependence of the ionization fraction on the relative baryon-CDM velocity during the epoch of reionization. However, we find that the B-mode signal is undetectably small for currently favored reionization models in which hydrogen is reionized promptly at z ~ 10, although constraints on this signal by future experiments may help constrain models in which partial reionization occurs at higher redshift, e.g. by accretion onto primordial black holes.

preprint2011arXiv

Testing Inflation with Dark Matter Halos

Cosmic inflation provides a mechanism for generating the early density perturbations that seeded the large-scale structures we see today. Primordial non-Gaussianity is among the most promising of few observational tests of physics at this epoch. At present non-Gaussianity is best constrained by the cosmic microwave background, but in the near term large-scale structure data may be competitive so long as the effects of primordial non-Gaussianity can be modeled through the non-linear process of structure formation. We discuss recent work modeling effects of a few types of primordial non-Gaussianity on the large-scale halo clustering and the halo mass function. More specifically, we compare analytic and N-body results for two variants of the curvaton model of inflation: (i) a "tauNL" scenario in which the curvaton and inflaton contribute equally to the primordial curvature perturbation and (ii) a "gNL" model where the usual quadratic fNL term in the potential cancels, but a large cubic term remains.

preprint2011arXiv

The Non-Gaussian Halo Mass Function with fNL, gNL and tauNL

Primordial non-Gaussianity has emerged as one of the most promising probes of the inflationary epoch. While the cosmic microwave background and large-scale halo bias currently provide the most stringent constraints on the non-Gaussian parameter fNL, the abundance of dark matter halos is a complementary probe which may allow tests of Gaussianity which are independent of the precise form of non-Gaussian initial conditions. We study the halo mass function in N-body simulations with a range of non-Gaussian initial conditions. In addition to the usual fNL model, we consider gNL Phi^3-type non-Gaussianity and models where the 4-point amplitude tauNL is an independent parameter. We introduce a new analytic form for the halo mass function in the presence of primordial non-Gaussianity, the "log-Edgeworth" mass function, and find good agreement with the N-body simulations. The log-Edgeworth mass function introduces no free parameters and can be constructed from first principles for any model of primordial non-Gaussianity.

preprint2010arXiv

CMB lensing and primordial non-Gaussianity

We study the effects of gravitational lensing on the estimation of non-Gaussianity from the bispectrum of the cosmic microwave background (CMB) temperature anisotropies. We find that the effect of lensing on the bispectrum may qualitatively be described as a smoothing of the acoustic features analogous to the temperature power spectrum. In contrast to previous results, for a Planck-like experiment which is cosmic-variance limited to L=2000, we find that lensing causes no significant degradation of our ability to constrain the non-Gaussianity amplitude fNL for both local and equilateral configurations, provided that the biases due to the cross correlation between the lensing potential and the integrated-Sachs-Wolfe (ISW) contribution to the CMB temperature are adequately understood. With numerical simulations, we also verify that low-order Taylor approximations to the lensed bispectrum and ISW-lensing biases are accurate.

preprint2010arXiv

Delensing CMB Polarization with External Datasets

One of the primary scientific targets of current and future CMB polarization experiments is the search for a stochastic background of gravity waves in the early universe. As instrumental sensitivity improves, the limiting factor will eventually be B-mode power generated by gravitational lensing, which can be removed through use of so-called delensing algorithms. We forecast prospects for delensing using lensing maps which are obtained externally to CMB polarization: either from large-scale structure observations, or from high-resolution maps of CMB temperature. We conclude that the forecasts in either case are not encouraging, and that significantly delensing large-scale CMB polarization requires high-resolution polarization maps with sufficient sensitivity to measure the lensing B-mode. We also present a simple formalism for including delensing in CMB forecasts which is computationally fast and agrees well with Monte Carlos.

preprint2008arXiv

CMBPol Mission Concept Study: Gravitational Lensing

Gravitational lensing of the cosmic microwave background by large-scale structure in the late universe is both a source of cosmological information and a potential contaminant of primordial gravity waves. Because lensing imprints growth of structure in the late universe on the CMB, measurements of CMB lensing will constrain parameters to which the CMB would not otherwise be sensitive, such as neutrino mass. If the instrumental noise is sufficiently small (<~ 5 uK-arcmin), the gravitational lensing contribution to the large-scale B-mode will be the limiting source of contamination when constraining a stochastic background of gravity waves in the early universe, one of the most exciting prospects for future CMB polarization experiments. High-sensitivity measurements of small-scale B-modes can reduce this contamination through a lens reconstruction technique that separates the lensing and primordial contributions to the B-mode on large scales. A fundamental design decision for a future CMB polarization experiment such as CMBpol is whether to have coarse angular resolution so that only the large-scale B-mode (and the large-scale E-mode from reionization) is measured, or high resolution to additionally measure CMB lensing. The purpose of this white paper is to evaluate the science case for CMB lensing in polarization: constraints on cosmological parameters, increased sensitivity to the gravity wave B-mode via lens reconstruction, expected level of contamination from non-CMB foregrounds, and required control of beam systematics.

preprint2006arXiv

Pseudo-$C_\ell$ estimators which do not mix E and B modes

Pseudo-$C_\ell$ quadratic estimators for CMB temperature and polarization power spectra have been used in the analysis pipelines of many CMB experiments, such as WMAP and Boomerang. In the polarization case, these estimators mix E and B modes, in the sense that the estimated B-mode power is nonzero for a noiseless CMB realization which contains only E-modes. Recently, Challinor & Chon showed that for moderately sized surveys ($f_{sky} \sim 0.01$), this mixing limits the gravity wave B-mode signal which can be detected using pseudo-$C_\ell$ estimators to $T/S \sim 0.05$. We modify the pseudo-$C_\ell$ construction, defining ``pure'' pseudo-$C_\ell$ estimators, which do not mix E and B modes in this sense. We study these estimators in detail for a survey geometry similar to that which has been proposed for the QUIET experiment, for a variety of noise levels, and both homogeneous and inhomogeneous noise. For noise levels $\simle 20$ $μ$K-arcmin, our modification significantly improves the B-mode power spectrum errors obtained using pseudo-$C_\ell$ estimators. In the homogeneous case, we compute optimal power spectrum errors using a Fisher matrix approach, and show that our pure pseudo-$C_\ell$ estimators are roughly 80% of optimal, across a wide range of noise levels. There is no limit, imposed by the estimators alone, to the value of $T/S$ which can be detected.

preprint2006arXiv

Pure pseudo-C_l estimators for CMB B-modes

Fast heuristically weighted, or pseudo-C_l, estimators are a frequently used method for estimating power spectra in CMB surveys with large numbers of pixels. Recently, Challinor & Chon showed that the E-B mixing in these estimators can become a dominant contaminant at low noise levels, ultimately limiting the gravity wave signal which can be detected on a finite patch of sky. We define a modified version of the estimators which eliminates E-B mixing and is near-optimal at all noise levels.

preprint2004arXiv

Weak Lensing of the CMB: Sampling Errors on B-Modes

The B modes generated by the lensing of CMB polarization are a primary target for the upcoming generation of experiments and can potentially constrain quantities such as the neutrino mass and dark energy equation of state. The net sample variance on the small scale B modes out to l=2000 exceeds Gaussian expectations by a factor of 10 reflecting the variance of the larger scale lenses that generate them. It manifests itself as highly correlated band powers with correlation coefficients approaching 70% for wide bands of Delta l/l \~0.25. It will double the total variance for experiments that achieve a sensitivity of approximately 4 uK-arcmin and a beam of several arcminutes or better. This non-Gaussianity must be taken into account in the analysis of experiments that go beyond first detection.