[Paper Review] Optimal analysis of the CMB trispectrum

Kendrick M. Smith, Leonardo Senatore+1 more|arXiv (Cornell University)|Feb 2, 2015

Cosmology and Gravitation TheoriesPhysics and Astronomy21 citations

TL;DR

This paper develops an optimal, computationally efficient framework for analyzing the CMB trispectrum—four-point correlation function of primordial density perturbations—using a factorizable representation of three distinct non-Gaussian trispectrum shapes arising from inflationary effective field theory. It introduces a three-parameter model (g^loc_NL, g^{\dot{\sigma}^4}_NL, g^{(\partial\sigma)^4}_NL) and applies it to WMAP data, finding no significant non-Gaussian signal after lensing correction, with constraints at the level of 10^5–10^6.

ABSTRACT

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}^{ m loc}$ of $ζ^3$-type local non-Gaussianity. This three-parameter space is analogous to the parameter space $(f_{NL}^{ m loc}, f_{NL}^{ m equil}, f_{NL}^{ m 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}^{ m loc} = (-3.80 \pm 2.19) imes 10^5$, $g_{NL}^{\dotσ^4} = (-3.20 \pm 3.09) imes 10^6$, and $g_{NL}^{(\partialσ)^4} = (-10.8 \pm 6.33) imes 10^5$ after correcting for the effects of CMB lensing. No evidence of a nonzero inflationary four-point function is seen.

Motivation & Objective

• To develop a general, optimal data analysis framework for the CMB trispectrum, the four-point function of primordial curvature perturbations.
• To identify and parameterize the three leading-order trispectrum shapes from inflationary effective field theory, consistent with symmetries and single-field inflation constraints.
• To enable fast, numerically stable computation of trispectrum estimators and simulations via a factorizable form.
• To apply the framework to WMAP data and derive minimum-variance constraints on inflationary trispectrum parameters, correcting for lensing effects.

Proposed method

• Proposes a three-parameter model for the trispectrum: g^loc_NL (local-type), g^{\dot{\sigma}^4}_NL, and g^{(\partial\sigma)^4}_NL, derived from quartic operators in the inflationary action.
• Introduces a factorizable representation of the trispectrum that enables fast evaluation of CMB estimators and non-Gaussian simulations via efficient matrix operations.
• Uses a Fisher matrix analysis to identify near-degeneracy among the three trispectra, justifying the three-parameter parameterization analogous to f^loc_NL, f^eq_NL, f^orth_NL for the bispectrum.
• Develops an optimal estimator ˆF using a contraction formalism with coefficients (α, β, γ) = (1/16, 9/16, -3/8) chosen to minimize variance and satisfy consistency conditions.
• Implements an end-to-end numerical convergence test using refined tolerance parameters to verify that numerical errors are observationally negligible.
• Applies the framework to WMAP 9-year data using both optimal and Monte Carlo pipelines, correcting for CMB lensing effects via a modified likelihood approach.

Experimental results

Research questions

• RQ1 What are the theoretically motivated, symmetry-consistent trispectrum shapes that arise from inflationary effective field theory at quartic order?
• RQ2 How can the CMB trispectrum be represented in a factorizable form to enable computationally efficient data analysis and simulation?
• RQ3 What are the optimal estimators for the three leading trispectrum parameters, and how can they be constructed to minimize variance?
• RQ4 How do lensing effects distort trispectrum measurements, and how can they be corrected in a data analysis pipeline?
• RQ5 What are the current observational constraints on inflationary trispectrum parameters from WMAP data, and do they show evidence of non-Gaussianity?

Key findings

• The three-parameter trispectrum model (g^loc_NL, g^{\dot{\sigma}^4}_NL, g^{(\partial\sigma)^4}_NL) captures the dominant non-Gaussian signals from inflationary effective field theory at quartic order.
• The Fisher matrix analysis reveals a near-degeneracy among the three trispectra, justifying the use of a three-dimensional parameter space for optimal analysis.
• The optimal estimator ˆF is constructed with coefficients (1/16, 9/16, -3/8) to minimize variance, with a semianalytic justification based on Wick contraction dominance.
• After correcting for CMB lensing, the WMAP data yield the constraints: g^loc_NL = (−3.80 ± 2.19) × 10^5, g^{\dot{\sigma}^4}_NL = (−3.20 ± 3.09) × 10^6, and g^{(\partial\sigma)^4}_NL = (−10.8 ± 6.33) × 10^5.
• No significant non-Gaussian signal is detected in the WMAP data, consistent with the null result expected in single-field inflation, with constraints at the 10^5–10^6 level.