Issue 
A&A
Volume 570, October 2014



Article Number  A94  
Number of page(s)  8  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201424179  
Published online  27 October 2014 
The Komatsu Spergel Wandelt estimator for oscillations in the cosmic microwave background bispectrum
^{1}
Sorbonne Universités, UPMC Univ. Paris 06, UMR 7095, Institut
d’Astrophysique de Paris,
75014
Paris,
France
email:
munchmey@iap.fr
^{2}
CNRS, UMR 7095, Institut d’Astrophysique de Paris,
75014
Paris,
France
^{3}
Lagrange Institute (ILP) 98bis boulevard Arago,
75014
Paris,
France
^{4}
Department of Physics, University of Illinois at
UrbanaChampaign, Urbana, IL
61801,
USA
^{5}
Department of Astronomy, University of Illinois at
UrbanaChampaign, Urbana, IL
61801,
USA
^{6}
Paris Centre for Cosmological Physics and Laboratoire
AstroParticule et Cosmologie, Université Paris 7Denis Diderot,
75205
Paris Cedex 13,
France
^{7}
African Institute for Mathematical Sciences,
7945
Muizenberg, South
Africa
Received:
9
May
2014
Accepted:
1
August
2014
Oscillating shapes of the primordial bispectrum are present in many inflationary models. The Planck experiment has recently published measurements of oscillating shapes, which were, however, limited to the efficient frequency range of the used analysis method. Here, we study the Komatsu Spergel Wandelt (KSW) estimator for oscillations in the cosmic microwave background bispectrum, that examines arbitrary oscillation frequencies for separable oscillating bispectrum shapes. We study the precision with which amplitude, phase, and frequency can be determined with our estimator. An examination of the threepoint function in real space gives further insight into the estimator.
Key words: cosmic background radiation / inflation / early Universe / cosmology: observations
© ESO, 2014
1. Introduction
The cosmic microwave background (CMB) provides the most direct experimental access to the statistics of the primordial curvature perturbations. In standard singlefield, slow roll inflation, these perturbations are Gaussian to a very good approximation (Maldacena 2003). However, more complicated models of inflation often predict detectable amounts of nonGaussianities of various shapes. In particular, the primordial bispectrum B(k_{1},k_{2},k_{3}) arising from the threepoint correlations of the curvature field ⟨ΦΦΦ⟩ is a sensitive probe to discriminate among models of inflation (see e.g. the reviews Liguori et al. 2010; Komatsu 2010; Yadav & Wandelt 2010). High precision measurements of the CMB recently provided by the Planck experiment (Planck Collaboration XXIV 2014) have been consistent with Gaussianity and set stronger limits on primordial nonGaussianities. However, the availability of such limits depend on the specific shape of the bispectra under consideration.
A class of bispectra that has attracted attention in recent years are oscillating shapes. Such bispectrum oscillations can arise in a variety of theoretical models. The authors of Chen et al. (2008) calculated the primordial bispectrum in the presence of features in the inflaton potential of standard single field inflation. They provided two analytical bispectrum shapes that approximate their results. The feature model oscillates linearly with the scale and is induced by sharp features in the inflaton potential. The resonance model includes oscillations with the logarithm of the scale and is induced by periodic features in the inflaton potential. More recently, the authors of Bartolo et al. (2010) used the effective field theory of inflation to examine the influence of sharp features in the inflaton potential, which also provide oscillating bispectrum solutions. Periodical modulations of the inflaton potential appear, for example, in axion monodromy inflation models (Flauger et al. 2010; Flauger & Pajer 2011). Certain bispectrum shapes motivated by NonBunchDavis vacua also include oscillations (e.g. Chen et al. 2007; Meerburg et al. 2009), as do cascade inflation models (Ashoorioon & Krause 2006). A transient reduction in the speed of sound also leads to oscillations (Achucarro et al. 2014b,a). Oscillations in the bispectrum are usually accompanied by oscillations in the primordial power spectrum. Recent searches for power spectrum oscillation were presented in Pahud et al. (2009); Peiris et al. (2013); Planck Collaboration XXII (2014); Easther & Flauger (2014); Meerburg et al. (2014); Meerburg & Spergel (2014), but no statistically significant result has been found. Combining power spectrum and bispectrum measurements can improve the sensitivity.
In the present work, we focus on the feature model shape, since it is simple and approximates some more complicated oscillating shapes. In particular, the feature model has the important property of separability. The Planck paper on nonGaussianities (Planck Collaboration XXIV 2014) already included a targeted search for the feature and resonant shapes with the modal expansion method. In this methodology, the bispectrum under consideration is expanded into a basis of separable shapes (Fergusson et al. 2010), which allows an efficient numerical estimation by the Komatsu Spergel Wandelt (KSW) estimator (Komatsu et al. 2005). However, the separable basis functions used by Planck did not allow to represent high frequency oscillations, limiting the frequency range, which could be searched for oscillations. The situation can be improved by using a set of oscillating basis functions for the modal expansion (Meerburg 2010). In this work, we follow a more direct approach by writing the feature model bispectrum in separable form, which makes it possible to search for oscillations of arbitrary frequency, which are only limited by the resolution of the maps. The gain in computational time with respect to the modal expansion is of the order of the number of modes that would be necessary to approximate the shape. We study the properties of the estimator in detail, including the precision with which phase and frequency of the oscillation can be determined. We also give an illuminating interpretation of the KSW estimator for oscillations in position space.
2. Oscillations in the primordial bispectrum
2.1. Bispectrum shape and experimental constraints
The general translation and rotation invariant primordial bispectrum of the curvature potential Φ can be written as (1)where the bispectrum B_{Φ} is a function of the magnitude of the wave numbers k_{i} and f_{NL} is the amplitude of the bispectrum. The feature model bispectrum that we are primarily interested in is (2)It oscillates linearly with the mean, , of the wave numbers. Here Δ_{Φ} is the primordial power spectrum amplitude and the 1 /k^{6} factor compensates for the phase space factor. The bispectrum is parametrised by the amplitude f_{NL} by the oscillation scale k_{c} and by the phase φ. The oscillation scale k_{c} implies an efficient multipole periodicity of l_{c} ≃ k_{c}[τ − τ_{rec}], where τ − τ_{rec} is the conformal distance to recombination.
Planck has searched for the feature bispectrum shape for sample frequencies in the range 0.01 <k_{c}< 0.1 at four different phases φ = 0,π/ 4,π/ 2,3π/ 4Planck Collaboration XXIV (2014). Here, k_{c} = 0.01 corresponds to an effective multipole periodicity l_{c} = 140. The best fit model has k_{c} = 0.0185 (l_{c} = 260) and phase Φ = 0 with a significance of 3σ. This may correspond to weak hints for oscillation in the full bispectrum reconstruction, which were found for l< 500. However, the statistical significance becomes much lower when one takes the number of statistically uncorrelated feature models into account that were searched. We note that the range of the oscillation frequency was constrained because of the limitations of the analysis method. With the present work, we target the unexplored range k_{c}< 0.01 in particular.
2.2. Separability
For convenience, we rewrite the feature model (2) as a sum of sine and cosine contributions as (3)where and . This can be written in separated form as (4)where we have defined (5)This separability property allows efficient computation and estimation of the bispectrum. One may also include an exponential decay factor of the form , while retaining separability with identical formulas up to trivial replacements.
3. Oscillations in the CMB bispectrum
From the separable expression (4) for the primordial bispectrum, one can calculate the CMB bispectrum with the standard lineofsight integration method Seljak & Zaldarriaga (1996). The reduced bispectrum is then (6)where we have defined the functions and where j_{l} are spherical Bessel functions and Δ_{l} are the CMB transfer functions that we evaluate numerically with CAMB Lewis et al. (2000). We note that these functions depend on k_{c} and have to be evaluated for each oscillation frequency of interest. Both the transfer functions Δ_{l}(k) (in k space) and the Bessel functions j_{l} are highly oscillatory integrals, which must be evaluated with sufficient sampling. Examples of the X(r) and Y(r) functions are given in Fig. 1.
Fig. 1
Function X and Y of the feature model for k_{c} = 0.01 (top) and k_{c} = 0.001 (bottom) as a function of l for r = τ_{0} − τ_{rec}. X and Y have units of Mpc^{1}. 
To calculate the resulting CMB bispectrum, we perform the r integral in Eq. (6). We choose a quadrature of about 2000 points in r with a higher sampling in the range of recombination. The contribution of different values of r is examined in more detail in Sect. 6. To plot the bispectrum, it is convenient to normalise by the constant bispectrum with natural k^{6} scaling, as proposed in Fergusson & Shellard (2009). The constant primordial bispectrum is given by (9)and its large angle SachWolfe CMB solution is Fergusson & Shellard (2009)(10)In Fig. 2, we show a simple 1dimensional visualisation, the equal l bispectrum b_{lll} normalised by for different frequencies and phases of the feature model.
Fig. 2
Equal l reduced bispectrum normalised by the large angle solution of the constant bispectrum for feature models with different parameters. 
4. The KSW estimator for the oscillating bispectrum
4.1. KSW estimator
The optimal KSW estimator for a sum of bispectra b^{i} in the presence of a mask and nonuniform noise is (Komatsu et al. 2005; Babich 2005; Komatsu 2010) (11)with (12)and with the Fisher matrix given by (13)In the present case of an oscillation with a single scale k_{c}, the bispectrum is a sum b = f_{1}b^{1} + f_{2}b^{2} of sine and cosine components, as given by Eq. (6).
To numerically evaluate the KSW estimator terms S_{i} we need the weighted maps (14)The cubic KSW estimator, which is exact in the case of a full sky observation, is given in terms of these maps by (15)Partial sky coverage can be taken into account by incorporating the linear term of the KSW estimator, so that with (16)where the expectation values have to be evaluated by Monte Carlo averaging over Gaussian realisations drawn with the same beam, mask, and noise properties as expected in the data. To make the KSW estimator optimal for a nonuniform sky coverage, it is necessary to perform an inverse covariance weighting with the nondiagonal covariance matrix. This is a computationally challenging problem (Smith et al. 2009; Elsner & Wandelt 2012). It was noted in Planck Collaboration XXIV (2014) that one can also achieve excellent results by assuming a diagonal covariance matrix Ĉ_{l} = C_{l} + N_{l}, where N_{l} assumes homogeneous noise, and by using a diffusive inpainting on the masked areas. In this approximation, the Fisher matrix scales proportionally to the visible fraction of the sky f_{sky}.
For the remainder of this paper, we assume full sky coverage so that (17)where .
From the inverse Fisher matrix (i.e. the covariance), one obtains the correlation between the sine and cosine terms. For example, for l = 1000 and k_{c} = 0.01, the correlation matrix is (18)which shows a weak correlation as expected. From f_{1} and f_{2} of Eq. (6), we can calculate the amplitude and phase of the oscillation. The variance of the quantities f_{NL} and Φ can be calculated by error propagation from f_{1} and f_{2}. For example (19)In the approximation of a diagonal covariance matrix and σ_{f1} = σ_{f2}, this gives σ_{f} = σ_{f1}. In this approximation the sensitivity on the phase depends on f as .
The Fisher matrix allows us to forecast the precision that can be obtained on the bispectrum parameters f_{1},f_{2}. With the approximation that the f_{1},f_{2} covariance matrix is a multiple of the unit matrix, the Fisher forecast on f_{1},f_{2} equals the forecast on f_{NL}. The precision on f_{NL} is then given by . For a noiseless fullsky experiment, the Fisher forecast for the feature model is shown in Fig. 3 for a number of different oscillation frequencies.
Fig. 3
Fisher forecast of σ_{fNL} for the feature model for different k_{c}, assuming a noiseless full sky experiment. 
4.2. Estimating the frequency
The estimator described above explicitly estimates the amplitude f_{NL} and phase φ of the oscillation for a fixed frequency k_{c}. To estimate the oscillation frequency k_{c}, it is necessary to sample the frequency space with the KSW estimator and search for peaks in the significance of the estimated amplitude. We assume the primordial bispectrum to be given by a single oscillation frequency and not a spectrum of contributions. We consider only the sine component of the bispectrum first, meaning that the phase is φ = 0 (see below for the generalisation to the phase). In this case, the estimator for a frequency k_{i} is (20)with covariance (in the usual Gaussian approximation) The Fisher matrix is given by Eq. (13), where the index i now goes over frequency sampling points. An example of the Fisher matrix F_{ij} is shown in Fig. 4 (top). The corresponding correlation matrix is , shown in Fig. 4 (middle).
Fig. 4
Top: fisher matrix F_{ij} for 100 logarithmically spaced frequencies between k_{c} = 0.002 and k_{c} = 0.01. Middle: corresponding correlation matrix . Bottom: correlation matrix onedimensional slice corr(f_{i},f_{j}) for k_{i} = 0.005. 
A onedimensional slice of the correlation matrix is shown in Fig. 4 (bottom) for k_{c} = 0.005. The plot shows strong antipeaks to both sides of the maximum and several small secondary peaks. This is not an artefact of the chosen estimator, but the physical overlap of the CMB bispectra induced by different primordial oscillation frequencies. The frequency sampling must be at least sufficient to resolve the peak structure of this plot. However, the width of the central peak does not directly limit the precision σ_{k} with which the primordial frequency can be determined.
To evaluate the precision with which k_{c} can be estimated, we note that the correlation matrix is identical to the bispectrum correlator, (23)which is the usual measure to discriminate bispectra (see e.g. Fergusson & Shellard 2009). It gives the estimated proportion of f that is recovered when estimating a spectrum B′ when the true underlying spectrum is B. If the underlying bispectum is f_{kc}B_{kc}, the expectation value of the estimated amplitude at a different frequency k is thus ). The variance at each data point is independent of the signal. Thus, the correlation matrix approximates the estimates in the k range where the signal dominates the variance (compare Fig. 4 with the estimation example in Fig. 6).
Knowing the means f(k) and the covariance matrix Cov(k_{1},k_{2}) and using the property that each estimator is Gaussian distributed, we can write the continuum likelihood for the estimated amplitudes , if the true bispectrum is f_{kc}B_{kc}: (24)with mean (25)This likelihood could be explored by Monte Carlo to find maximum likelihood estimates of f_{kc} and k_{c} if a significant peak of f_{NL} is found in the spectrum. The Fisher matrix to forecast optimal precision on the estimated parameters is (26)where θ ∈ { f_{kc},k_{c} }. This can be integrated numerically for any given fiducial parameters.
In the above discussion, we assumed that the bispectrum only has a sine contribution. The generalisation of Eq. (24) to a free phase is (27)with the mean (28)The correlation matrix can be split into terms of f_{1} and f_{2} for efficient evaluation of the likelihood.
As we have seen, the estimation of the frequency requires to run a large number of estimators (around 100 to cover the frequency interval k_{c} = 0.001 to k_{c} = 0.01). For each of these estimators, it is necessary to calculate the linear term in Eq. (16) via a Monte Carlo averaging procedure over Gaussian map realisations with the same mask and noise properties as present in the experiment. Convergence of the linear term is usually achieved with 100 Monte Carlo realisations, although several hundreds can be used to improve accuracy. After the optimisation of the conformal distance integral, that is reviewed in Sect. 6, one can expect, with an angular resolution of l_{max} = 2000, a calculation time of about one day on a single cpu for a single frequency, including the Monte Carlo averaging. On a computation grid, one can thus easily cover the frequency range of interest.
Fig. 5
Parameter estimation histograms for 100 maps with linear model nonGaussianity with k_{c} = 0.01 and l_{max} = 1000. The maps were created with f_{NL} = 500 and φ = 0° (first two rows) and φ = 45° (rows three and four). Top: reconstructed f_{1} (green) and f_{2} (blue) amplitudes for φ = 0°. Second row: corresponding reconstructed amplitudes f_{NL} and phases φ. Third and fourth row: same as above but with phase φ = 45°. 
5. Map making
To verify the implementation and unbiased nature of the estimator, it is useful to be able to generate maps with the bispectrum signature of interest. The authors of Smith & Zaldarriaga (2011) introduced an algorithm that generates maps for arbitrary bispectra in the weak nonGaussian limit. A map is constructed from a linear combination of a Gaussian and a nonGaussian contribution as . The straightforward implementation of the algorithm gives a nonGaussian contribution of the form (29)It is thus easy to create maps of arbitrary phase from the sine and cosine terms.
An example of two sets of 100 maps, created and then estimated with the algorithms presented here, is shown in Fig. 5. The means and variances in these histograms are compatible with their Fisher forecast. We note that the distribution of f_{NL} is not Gaussian but follows a Rayleigh distribution, since f_{NL} represents the magnitude of the vector of the two directional components f_{1} and f_{2}.
An example of a frequency sweep can be found in Fig. 6 for the phase ϕ = 0°. It shows the secondary peaks that are expected from the correlation of different frequencies.
Fig. 6
Left: frequency sweep over a map with k_{c} = 0.005 and f_{NL} = 2000. Right: as left but in units of σ_{fNL}. 
6. Computational speed improvement by optimisation of the rintegral
Due to the large number of frequencies that have to be sampled with the corresponding estimator, the search for oscillations is computationally challenging. The time critical steps are the calculation of the Fisher matrix and the necessity of calculating a large number of r dependent KSW filtered maps. The latter problem becomes even more severe if one has to estimate many Monte Carlo generated maps for the calculation of the linear term. The situation can be improved by an analysis of the Fisher matrix, as shown in Smith & Zaldarriaga (2011).
For a separable shape, the Fisher matrix can also be expressed by a sum over contributions of different r sampling points that arise when numerically evaluating the bispectrum integral in Eq. (6). The total Fisher matrix is then given as a sum , where F_{ij} is the Fisher matrix element between the sampling points i and j and N_{fact} is the number of sampling points. The Fisher matrix elements are then given by (30)We now explicitly consider the sine term of the linear model (the cosine term is analogous), where the contribution of a distance r_{i} is given by (31)To give an impression of the bispectrum contribution of different distances r_{i}, we plot the diagonal elements F_{ii} in Fig. 7. The plots show the contribution of recombination (r = 14 000 Mpc), reionisation (r ≃ 10 500 Mpc), and ISW (r> 5000 Mpc). As expected, the dominant contribution comes from the time around recombination, and one can get a good approximation to the integral by sampling only a window around recombination. This is particularly useful to quickly scan a parameter space, in the present case the oscillation frequency, k_{c}.
In Smith & Zaldarriaga (2011), it was shown that one can go further and optimise the r sampling points to find a quadrature with surprisingly few sampling points that give an almost identical estimator. Their algorithm constructs a new bispectrum B′ from the original bispectrum B by choosing a subsample of points and weighting them so that the Fisher distance between the two is minimised: (32)This means that bispectrum values with a small signaltonoise are allowed to be very different. Using this algorithm, as an example, we obtain an approximate bispectrum B′ consisting of 30 sample points leading to a separability between B and B′ of 0.1σ assuming k_{c} = 0.01 and f_{NL} = 1000.
The most computationally demanding task in the estimation pipeline remains the calculation of the Fisher matrix in Eq. (30). It was also shown in Smith & Zaldarriaga (2011) that one can factorise this equation by inserting the integral representation of the Wigner symbol (33)The integral over z can be computed efficiently by Gauss Legendre integration. However, even with this expression, the calculation of F_{ij} needs many CPU hours depending on the chosen initial quadrature point number. If one does not want to calculate an optimised quadrature at every frequency points, but only wants to know the normalisation F(k_{c}) of the estimators, one can calculate this normalisation on a much wider frequency spacing and interpolate in between. This can be seen from the frequency dependent Fisher matrix in Fig. 4, where the diagonal elements vary slowly.
Fig. 7
F_{ii} as a function of conformal distance r_{i}. Top: k_{c} = 0.01, bottom: k_{c} = 0.005. 
7. A position space interpretation of the KSW estimator for oscillations
In Appendix A, we show that the threepoint function of the feature model in position space peaks for configurations that lie on a circumcircle of radius sin(θ) = 1 / (2k_{0}η), where in the convention of Eq. (2). This suggests searching for bispectrum oscillations in real space by convoluting the CMB map with a ring kernel of varying radius.
For a radially symmetric kernel, the convolution can be done efficiently in harmonic space as s_{lm} = K_{l}r_{lm}, where the kernel is given by the Legendre transformation, (34)and K(z) is a narrow window function in z = cos(θ). The estimate ℰ^{ring}(z) is given by the sum over the pixels of the cube of the convoluted map, (35)Figure 8 shows an example of this estimator for a map that was simulated with the algorithm of the preceding section. The plots show , which means that the estimate from the Gaussian map was subtracted. The red line shows the predicted position of the maximum. We note that this maximum would be harder to locate in real data where the Gaussian contribution cannot simply be subtracted.
Fig. 8
Ring estimate . Top left: k_{c} = 0.01. Top right: k_{c} = 0.005. Bottom: k_{c} = 0.001. The kernel width in all plots was Δθ = 0.005, f_{NL} was 10−20 times the optimal Fisher error. The red line shows the expected position of the maximum. 
It is interesting to compare the “intuitive” ring kernel K(θ) with the KSW kernel that is known to give optimal results. The KSW estimator is of form, (36)The largest contribution to this integral comes from decoupling at r_{rec}. We plot the KSW kernel function at decoupling X_{rec}(θ) in Fig. 9. It shows the maximum at the expected position and has roughly the expected window shape.
Fig. 9
KSW kernel X(θ) at decoupling radius r_{rec}. Top left: k_{c} = 0.01. Top right: k_{c} = 0.005. Bottom: k_{c} = 0.001. The red line is the predicted maximum of the kernel. 
8. Conclusion
In this paper, we have presented and extensively studied the KSW estimator for linear bispectrum oscillations. The main motivation for this approach is that the oscillating bispectrum shapes are difficult to represent with a modal expansion and, thus, have not yet been constrained at high oscillation frequency. We have provided the equations for estimation and map making and validated them with Monte Carlo simulations. Unlike many of the well know bispectum shapes, oscillations have two free parameters in addition to the common amplitude parameter f_{NL}. We have developed the methodology to estimate and constrain the oscillation phase φ and the frequency k_{c}. Our work will therefore allow one to explore a parameter space that was not previously accessible for a theoretically wellmotivated bispectrum shape. Finally, we have found an interesting position space interpretation of the KSW estimator for oscillations, based on an approximate evaluation of the corresponding threepoint function.
Acknowledgments
The authors acknowledge support from NSF Grant NSF AST 0908693 ARRA. B.D.W. acknowledges funding from an ANR Chaire d’Excellence (ANR10CEXC00401), the UPMC Chaire Internationale in Theoretical Cosmology, and NSF grants AST0908 902 and AST0708849. This work made in the ILP LABEX (under reference ANR10LABX63) was supported by French state funds managed by the ANR within the Investissements d’Avenir programme under reference ANR11IDEX000402. M.M. acknowledges funding by Centre National d’Etudes Spatiales (CNES). The authors would like to thank Daan Meerburg for useful discussions and comments.
References
 Achucarro, A., Atal, V., Hu, B., Ortiz, P., & Torrado, J. 2014a, Phys. Rev. D, 90, 023511 [NASA ADS] [CrossRef] [Google Scholar]
 Achucarro, A., Atal, V., Ortiz, P., & Torrado, J. 2014b, Phys. Rev. D, 89, 103006 [NASA ADS] [CrossRef] [Google Scholar]
 Adshead, P., Dvorkin, C., Hu, W., & Lim, E. A. 2012, Phys. Rev. D, 85, 023531 [NASA ADS] [CrossRef] [Google Scholar]
 Ashoorioon, A., & Krause, A. 2006 [arXiv:hepth/0607001] [Google Scholar]
 Babich, D. 2005, Phys. Rev. D, 72, 043003 [NASA ADS] [CrossRef] [Google Scholar]
 Bartolo, N., Fasiello, M., Matarrese, S., & Riotto, A. 2010, JCAP, 1008, 008 [NASA ADS] [CrossRef] [Google Scholar]
 Bashinsky, S., & Bertschinger, E. 2001, Phys. Rev. Lett., 87, 081301 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 Chen, X., Huang, M.X., Kachru, S., & Shiu, G. 2007, JCAP, 0701, 002 [NASA ADS] [CrossRef] [Google Scholar]
 Chen, X., Easther, R., & Lim, E. A. 2008, JCAP, 0804, 010 [NASA ADS] [CrossRef] [Google Scholar]
 Easther, R., & Flauger, R. 2014, JCAP, 1402, 037 [NASA ADS] [CrossRef] [Google Scholar]
 Elsner, F., & Wandelt, B. D. 2012 [arXiv:1211.0585] [Google Scholar]
 Fergusson, J., & Shellard, E. 2009, Phys. Rev. D, 80, 043510 [NASA ADS] [CrossRef] [Google Scholar]
 Fergusson, J., Liguori, M., & Shellard, E. 2010, Phys. Rev. D, 82, 023502 [NASA ADS] [CrossRef] [Google Scholar]
 Flauger, R., & Pajer, E. 2011, JCAP, 1101, 017 [NASA ADS] [CrossRef] [Google Scholar]
 Flauger, R., McAllister, L., Pajer, E., Westphal, A., & Xu, G. 2010, JCAP, 1006, 009 [NASA ADS] [CrossRef] [Google Scholar]
 Jackson, M. G., Wandelt, B., & Bouchet, F. 2014, Phys. Rev. D, 89, 023510 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E. 2010, Class. Quant. Grav., 27, 124010 [NASA ADS] [CrossRef] [Google Scholar]
 Komatsu, E., Spergel, D. N., & Wandelt, B. D. 2005, ApJ, 634, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Lewis, A., Challinor, A., & Lasenby, A. 2000, ApJ, 538, 473 [NASA ADS] [CrossRef] [Google Scholar]
 Liguori, M., Sefusatti, E., Fergusson, J. R., & Shellard, E. 2010, Adv. Astron., 2010, 980523 [NASA ADS] [CrossRef] [Google Scholar]
 Maldacena, J. M. 2003, JHEP, 0305, 013 [NASA ADS] [CrossRef] [Google Scholar]
 Meerburg, P. D. 2010, Phys. Rev. D, 82, 063517 [NASA ADS] [CrossRef] [Google Scholar]
 Meerburg, P. D., & Spergel, D. N. 2014, Phys. Rev. D, 89, 063537 [NASA ADS] [CrossRef] [Google Scholar]
 Meerburg, P. D., Spergel, D. N., & Wandelt, B. D. 2014, Phys. Rev. D, 89, 063536 [NASA ADS] [CrossRef] [Google Scholar]
 Meerburg, P. D., van der Schaar, J. P., & Corasaniti, P. S. 2009, JCAP, 0905, 018 [NASA ADS] [CrossRef] [Google Scholar]
 Pahud, C., Kamionkowski, M., & Liddle, A. R. 2009, Phys. Rev. D, 79, 083503 [NASA ADS] [CrossRef] [Google Scholar]
 Peiris, H., Easther, R., & Flauger, R. 2013, JCAP, 1309, 018 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XXII. 2014, A&A, 571, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIV. 2014, A&A, 571, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Seljak, U., & Zaldarriaga, M. 1996, ApJ, 469, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, K. M., & Zaldarriaga, M. 2011, MNRAS, 417, 2 [NASA ADS] [CrossRef] [Google Scholar]
 Smith, K. M., Senatore, L., & Zaldarriaga, M. 2009, JCAP, 0909, 006 [NASA ADS] [CrossRef] [Google Scholar]
 Yadav, A. P., & Wandelt, B. D. 2010, Adv. Astron., 2010, 565248 [NASA ADS] [CrossRef] [Google Scholar]
Appendix A: Angular correlation function of the feature model
In this appendix, we show that the linear feature model bispectrum peaks in real space for a special class of threepoint function configurations. A similar analysis was presented in Adshead et al. (2012). The corresponding calculation for logarithmic oscillations can be found in Jackson et al. (2014).
In the approximation of instantaneous CMB decoupling at time η, the CMB temperature perturbation is given in terms of the potential φ as Bashinsky & Bertschinger (2001), (A.1)The transfer function T_{rad}(k) is generally a complicated function of scale and cosmological parameters. For simplicity of an analytic answer, which accounts for finite resolution, we take it to be .
The positionspace primordial bispectrum is given by (A.2)where the correlation is of the form,
Since k_{1} + k_{2} + k_{3} = 0, the kspace correlations are categorised by the triangle formed by the k_{i}.
We now examine this integral for the linear model,
The delta function that couples the momenta can be written as
where the factor of η has been included for future convenience. Writing the sine as exponentials and performing the integral over angles, we obtain Defining the dimensionless parameter x_{i} ≡ kη, the momentum integral is This integral can be evaluated exactly in terms of a hypergeometric function, but there is a simplifying limit we can take. The lowmomentum, longdistance approximation allows
The integral then has the simple analytic solution In the ηk_{D} ≫ 1 limit, this maximally peaks when all three products peak near
Squaring then subtracting, we obtain
We can take the three vectors N_{12}, N_{23}, and N_{31} and arrange them in the plane, so that where
We now have the value of w for which the three factors are in resonance for a given configuration (n_{1},n_{2},n_{3}). The value of the factor at the resonance point is largest when is at its minimum, which is the case for ρ = 0. We conclude that the threepoint function peaks for triangle configurations (n_{1},n_{2},n_{3}) that lie on a circle with radius given by .
All Figures
Fig. 1
Function X and Y of the feature model for k_{c} = 0.01 (top) and k_{c} = 0.001 (bottom) as a function of l for r = τ_{0} − τ_{rec}. X and Y have units of Mpc^{1}. 

In the text 
Fig. 2
Equal l reduced bispectrum normalised by the large angle solution of the constant bispectrum for feature models with different parameters. 

In the text 
Fig. 3
Fisher forecast of σ_{fNL} for the feature model for different k_{c}, assuming a noiseless full sky experiment. 

In the text 
Fig. 4
Top: fisher matrix F_{ij} for 100 logarithmically spaced frequencies between k_{c} = 0.002 and k_{c} = 0.01. Middle: corresponding correlation matrix . Bottom: correlation matrix onedimensional slice corr(f_{i},f_{j}) for k_{i} = 0.005. 

In the text 
Fig. 5
Parameter estimation histograms for 100 maps with linear model nonGaussianity with k_{c} = 0.01 and l_{max} = 1000. The maps were created with f_{NL} = 500 and φ = 0° (first two rows) and φ = 45° (rows three and four). Top: reconstructed f_{1} (green) and f_{2} (blue) amplitudes for φ = 0°. Second row: corresponding reconstructed amplitudes f_{NL} and phases φ. Third and fourth row: same as above but with phase φ = 45°. 

In the text 
Fig. 6
Left: frequency sweep over a map with k_{c} = 0.005 and f_{NL} = 2000. Right: as left but in units of σ_{fNL}. 

In the text 
Fig. 7
F_{ii} as a function of conformal distance r_{i}. Top: k_{c} = 0.01, bottom: k_{c} = 0.005. 

In the text 
Fig. 8
Ring estimate . Top left: k_{c} = 0.01. Top right: k_{c} = 0.005. Bottom: k_{c} = 0.001. The kernel width in all plots was Δθ = 0.005, f_{NL} was 10−20 times the optimal Fisher error. The red line shows the expected position of the maximum. 

In the text 
Fig. 9
KSW kernel X(θ) at decoupling radius r_{rec}. Top left: k_{c} = 0.01. Top right: k_{c} = 0.005. Bottom: k_{c} = 0.001. The red line is the predicted maximum of the kernel. 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.