arxiv: v1 [astro-ph.co] 19 Jul 2018

Transcription

1 Foreground Biases on Primordia Non-Gaussianity Measurements from the CMB Temperature Bispectrum: Impications for Panck and Beyond J. Coin Hi 1, 2 1 Schoo of Natura Sciences, Institute for Advanced Study, Princeton, NJ, USA Center for Computationa Astrophysics, Fatiron Institute, New York, NY, USA arxiv: v1 [astro-ph.co] 19 Ju 2018 The cosmic microwave background (CMB) temperature bispectrum is currenty the most precise too for constraining non-gaussianity (NG) in the primordia curvature perturbations. The Panck temperature data tighty constrain the ampitude of oca-type NG: fnl oc = 2.5 ± 5.7. In this paper, we compute previousy-negected foreground biases in temperature-based fnl oc measurements. We consider signas from the integrated Sachs-Wofe (ISW) effect, gravitationa ensing, the therma (tsz) and kinematic Sunyaev-Ze dovich (ksz) effects, and the cosmic infrared background (CIB). In standard anayses, a significant foreground bias arising from the ISW-ensing bispectrum is subtracted from the fnl oc measurement. However, a number of other terms sourced by the ISW, ensing, tsz, ksz, and CIB fieds are aso present in the temperature bispectrum. We compute the dominant biases on fnl oc arising from these signas, focusing on squeezed bispectrum shapes. Most of the biases are non-backbody in nature, and are thus reduced by mutifrequency component separation methods; however, recent anayses have found that extragaactic foregrounds are present at nonnegigibe eves in the Panck component-separated maps. Moreover, the Panck FFP8 simuations do not incude the correations amongst components that are responsibe for these biases. We compute the biases for individua Panck frequencies, finding that some are comparabe to the statistica error bar on fnl, oc even for the main CMB channes (100, 143, and 217 GHz). For future experiments, they can greaty exceed the statistica error bar (considering temperature data ony). A fu assessment for Panck and other experiments wi require cacuations in tandem with component separation, ideay using simuations. Simiar biases wi aso affict measurements of equiatera and orthogona NG, as we as trispectrum NG. We concude that the search for primordia NG using Panck data may not yet be over. I. INTRODUCTION Primordia non-gaussianity (NG) is a key probe of the physics thought to have generated a structure in our Universe during its eariest moments. The simpest modes of infation (i.e., singe-fied, sow-ro) predict negigiby sma departures from Gaussianity in the primordia curvature perturbations [1, 2], but a rich spectrum of non-gaussian signatures can be produced in more compex infationary scenarios or non-infationary eary- Universe modes (see, e.g., Refs. [3, 4] for reviews). A key quantity of interest is the ampitude of the bispectrum of curvature perturbations in the so-caed squeezed imit (in which one wavenumber is much smaer than the other two, i.e., k 1 k 2, k 3 ), conventionay denoted as fnl oc [5]. Singe-fied, sow-ro infation predicts that fnl oc vanishes exacty (e.g., [6]), moduo sma, higher-order corrections due to the noninearity of gravity [7]. A detection of nonzero fnl oc woud rue out essentiay a singe-fied modes of infation [2, 8]. Athough we focus on oca-type NG in this paper, additiona shapes of NG (equiatera or orthogona) can aso contain a weath of information about the physics of the eary Universe [3, 4]. The most powerfu observabe for constraining primordia NG in current data sets is the bispectrum of temperature fuctuations in the cosmic microwave background jch@ias.edu (CMB) [9 12]. The we-understood, inear physics responsibe for the CMB anisotropy permits an essentiay direct mapping of the primordia curvature perturbations, thus aowing NG tempates to be directy fit to CMB maps. The most stringent current constraint on oca NG is derived from Panck data in this manner, yieding fnl oc = 2.5 ± 5.7 (temperature data ony) or = 0.8±5.0 (temperature and poarization data) [12]. f oc NL Nevertheess, a number of compex probems must be surmounted in order to extract robust NG constraints from CMB data. In this paper, we focus on one such probem: extragaactic foreground contamination in CMB temperature maps. Prior to NG anaysis, maps of the microwave sky at mutipe frequencies must be combined to extract a map of the CMB anisotropy from the mutitude of other sky signas, a process known as component separation (see Refs. [13 15] for an overview of the Panck CMB component separation methods). The goa of such techniques is to minimize the contributions from non-cmb contaminants (incuding both foregrounds and noise) whie preserving the CMB signa. However, such methods are generay imperfect, and some eve of foreground residuas wi propagate into the fina map. These residuas must be weunderstood in order for robust NG constraints to be obtained. We focus in this paper on extragaactic foreground contributions to oca-type NG estimators in CMB temperature maps, as detaied further beow. Most of these contributions are non-backbody in frequency dependence,

2 2 and can therefore be reduced (or, in some cases, competey removed [16]) via mutifrequency component separation. However, the extent to which such reduction for extragaactic foregrounds has occurred in the Panck component-separated CMB temperature maps is presenty uncear, and evidence has recenty accumuated that some sma-scae foregrounds may be present at non-negigibe eves [17, 18]. A crucia cross-check coud come from a poarization-ony NG anaysis, where the ony extragaactic foreground is that due to point source emission, which is we-understood and simpe to remove. However, as can be immediatey seen in the resuts quoted above, the Panck fnl oc constraints are strongy dominated by information in the CMB temperature fied. Upcoming experiments, incuding Simons Observatory 1 (SO) and CMB-S4 [19], may be sufficienty sensitive to aow independent measurements from temperature and poarization with comparabe error bars. At present, CMB temperature dominates in NG constraints. Moreover, some contaminants to NG estimators cannot be removed via mutifrequency component separation, as they possess the same backbody frequency dependence as the CMB itsef. Chief amongst these is the CMB temperature bispectrum sourced by the correation between the integrated Sachs-Wofe (ISW) effect [20, 21] and the gravitationa ensing potentia by which CMB photons are defected (see Ref. [22] for a review of CMB ensing). The ISW effect is the change in the temperature of CMB photons due to the decay (or enhancement) of ate-time gravitationa potentias (e.g., as a resut of dark energy). Both the ISW effect and gravitationa ensing do not ater the backbody spectrum of the CMB. In addition, both fieds trace the ate-time gravitationa potentia of the Universe. Finay, since CMB ensing coupes previousy-independent spherica harmonic coefficients of the CMB, the ISW-ensing correation produces a nonzero bispectrum in the CMB temperature fied [23 26]. This ISW-ensing bispectrum has a non-zero projection onto the bispectrum shape of oca-type NG, thereby producing a non-negigibe, irreducibe bias to estimates of fnl oc from CMB temperature maps. For Panck, this bias is substantia and must be subtracted to obtain unbiased constraints: f oc,isw φ NL,P anck = 7.6 (see IV A). In this paper, we point out the existence of additiona foreground biases in CMB temperature-derived fnl oc constraints. To our knowedge, the ISW-reated biases presented here have not been computed esewhere, with the exception of the ISW-ensing bias. Other biases invoving CMB ensing have received some attention, abeit imited [27, 28]. In genera, fnl oc biases are generated in the CMB temperature fied by bispectra invoving the ISW effect, CMB ensing, the cosmic infrared background (CIB), and the therma (tsz) and kinematic Sunyaev- Ze dovich (ksz) effects. 2 The CIB refers to the cumuative infrared emission of dusty, star-forming gaaxies over cosmic time, which has a broad redshift kerne peaking around z 2 (simiar to the CMB ensing kerne). The tsz effect is the inverse-compton scattering of CMB photons off hot, free eectrons, producing a shift in the photon spectrum to higher energies, and thus eaving a non-backbody spectra distortion in the CMB [29, 30]. The ksz effect is the Dopper-boosting of CMB photons scattering off eectrons that have a non-zero ine-of-sight (LOS) veocity in the CMB rest frame [31 33]. Aong with the ISW and CMB ensing signas, these fieds a trace the arge-scae structure of the universe in some way. The resuting correations generate bispectra in CMB temperature maps. Whie these bispectra are interesting in their own right for astrophysica and cosmoogica reasons [34 40], here we focus on their roe in biasing measurements of oca-type primordia NG, akin to the bias due to the ISW-ensing bispectrum described above. The basic concusions of this paper are summarized in Tabe I, which gives the bias on fnl oc sourced by severa types of foreground bispectra at the primary CMB channes of Panck (note that ground-based experiments observe at simiar frequencies). The tabe incudes resuts for Panck noise eves, corresponding to an experiment that is cosmic variance-imited to a maximum mutipoe max = 1590, as we as for future experiments that wi be cosmic variance-imited to max = 3000 (abeit not covering the fu sky). The primary takeaway from these resuts is that even for the main CMB channes of Panck, the biases due to these foregrounds are comparabe to the 1σ error bar on fnl oc from the temperature bispectrum. Of course, component separation wi reduce the non-backbody contributions to some extent (and a contributions wi be summed in a compex manner in the synthesis of the fina CMB map), but a fu modeing of this procedure is beyond the scope of this paper and best impemented via simuations. The other important concusion from Tabe I is that for future experiments aiming to constrain primordia NG from the CMB bispectrum (e.g., SO and CMB-S4), nonbackbody foregrounds must be ceaned very precisey to avoid biases from the tsz and CIB signas. Furthermore, a non-negigibe bias due to the backbody ISW-kSZ-kSZ bispectrum must be subtracted, which has not been computed in the iterature to date. Given these compexities, future NG constraints may be better off reying on CMB poarization, which is much ceaner than temperature on sma scaes. The remainder of this paper is organized as foows. In II, we motivate this work in the context of current measurements and previous cacuations. In III, we f oc NL Bispectra invoving point sources aso generate biases, but are smaer in magnitude, particuary for Panck, and have been considered in previous anayses.

3 3 Biases on f oc NL for Panck T T T : max = 1590, σ(f oc NL) = 5.0/ f sky Frequency ISW-φ ( IV A) tsz-φ ( IV B) CIB-φ ( IV C) ISW-tSZ-tSZ ( V A) ISW-CIB-CIB ( V B) ISW-kSZ-kSZ ( V D) 100 GHz GHz GHz Biases on f oc NL for SO/CMB-S4 T T T : max = 3000, σ(f oc NL) = 2.6/ f sky Frequency ISW-φ ( IV A) tsz-φ ( IV B) CIB-φ ( IV C) ISW-tSZ-tSZ ( V A) ISW-CIB-CIB ( V B) ISW-kSZ-kSZ ( V D) 100 GHz GHz GHz TABLE I: Summary of extragaactic foreground biases on fnl oc measurements from the CMB temperature bispectrum. The top pane gives the biases for an experiment that is cosmic variance-imited to a maximum mutipoe max = 1590 (i.e., Panck), whie the bottom pane contains the resuts for an experiment with max = 3000 (i.e., Simons Observatory or CMB-S4). These correspond to σ(fnl) oc = 5.0 and 2.6, respectivey, assuming a fu-sky measurement (in practice, f sky 0.75 for Panck and f sky 0.4 for SO and CMB-S4). We perform these cacuations using the measured Panck bandpasses centered at 100, 143, and 217 GHz, which are the primary CMB channes of Panck and ground-based experiments, athough higher frequencies are considered ater in the paper as we. For brevity, we do not incude the ISW-tSZ-CIB bispectrum resuts in this tabe, but these biases are computed in V C. provide reevant theoretica background reated to primordia NG and formaism for bispectra. The foowing two sections contain the main resuts of the paper: IV describes biases reated to the correation of CMB ensing with other secondary fieds, whie V describes biases generated by the correation of these fieds with the ISW effect. We discuss the impications of these resuts in VI and concude in VII. We assume a standard fat ΛCDM cosmoogy throughout this paper, adopting the foowing parameters: matter density Ω m = 0.277, Hubbe constant H 0 = 70.2 km/s/mpc, baryon density Ω b = , spectra tit n s = 0.962, and ampitude of density fuctuations σ 8 = Our concusions are weaky sensitive to the assumed vaues of these parameters, but we comment beow on a few notabe exceptions where they are important (reated to the tsz signa). II. MOTIVATION AND CONTEXT Non-Gaussian signas from extragaactic foregrounds in CMB temperature maps have received significant attention in recent years due to the important biases they can generate in reconstructed CMB ensing maps (e.g., [17, 41 45]). These biases arise from both non-backbody (e.g., tsz and CIB) and backbody (e.g., ksz) foregrounds; accordingy, some can be reduced by mutifrequency component separation, whie some cannot be (athough other mitigation methods can be empoyed). Simiary, these non-gaussian signas wi generate biases in CMB temperature bispectrum estimates of primordia NG. To date, the ISW-ensing bispectrum has received significant attention in this context [23 26], but other biases have been ess studied, if at a. For exampe, the Panck 2015 NG anaysis considers ony the ISW-ensing bispectrum, point source bispectrum, and custered CIB auto-bispectrum as contaminants to primordia NG measurements [12]. An eary estimate of the tsz-ensing bispectrum bias on fnl oc was presented in Ref. [27], but it is uncear what observationa frequencies are considered in the anaysis, and theoretica modeing of the tsz signa has significanty evoved in the intervening decade. More recenty, the CIB-ensing bispectrum bias on fnl oc was considered in Ref. [28]; their resuts for specific Panck frequency channes are in quaitative agreement with those presented in IV C beow. However, to our knowedge, no cacuations of ISW-reated biases have been presented in the iterature to date, beyond the ISW-ensing bias. Amongst these contributions is that due to the ISW-kSZ-kSZ bispectrum, which is backbody in frequency dependence and therefore must be subtracted from observationa estimates, ike the ISWensing bias. We provide a first estimate of this bias on in V D. f oc NL However, non-backbody biases are aso of significant concern. Recent anayses indicate that non-negigibe extragaactic foreground contamination has eaked into the Panck component-separated CMB temperature maps, which form the foundation of the Panck NG anaysis. For exampe, Fig. 1 of Ref. [17] demonstrates that the tsz signa of opticay-seected gaaxy custers is present in the Panck SMICA CMB map with an ampitude neary identica to that seen in the Panck 143 GHz map. Simiary, Ref. [18] uses cross-correations with optica gaaxy survey data to detect the presence of tsz residuas in the Panck NILC CMB map at 54σ significance. Their overa estimate is that roughy haf of the 143 GHz tsz

4 4 signa is present in the NILC map. Whie these studies are based on observationa estimates, it woud be more robust to estimate the eakage via end-to-end simuations containing a reevant sky signas, an anaysis which has not yet been performed. In this context, we note that athough the Panck 2015 NG anaysis used the Fu Foca Pane 8 (FFP8) simuations to perform end-to-end tests of their anaysis pipeines, the FFP8 simuations do not contain the signas responsibe for the biases considered in this paper (with the exception of the ISW-ensing bias) [46]. In particuar, the ISW fied is generated ony as part of the primary CMB map (i.e., as a Gaussian random fied), and is not correated with the tsz, ksz, or CIB fieds. Simiary, the CMB ensing fied is not correated with the tsz, ksz, or CIB fieds, nor is the CIB fied correated with the tsz or ksz fieds. The tsz and ksz fieds are partiay correated with one another, but not with any of the other secondary fieds. The ISW-ensing correation is generated by the agorithm with which the primary CMB is gravitationay ensed (because the T φ correation is incuded in CAMB [47] 3 power spectra). Thus, the NG pipeine verification tests run on the FFP8 simuations do not test for any of the biases considered in this paper (except for the ISW-ensing bias). This situation coud be remedied by using sky simuations in which the extragaactic fieds are propery correated with one another (e.g., [48]). Motivated by the existence of non-negigibe extragaactic foreground contamination in the Panck component-separated CMB temperature maps and the absence of neary a reevant foreground bispectra in the FFP8 simuations, we consider the roe that these effects might have on estimates of primordia NG. We note that poarization-ony anayses woud be amost entirey immune to these foreground biases, but the Panck poarization data are not sufficienty sensitive for such a test (σ(fnl oc ) from poarization data aone, in comparison to σ(fnl oc ) 5 6 from temperature data aone [12]). We focus on oca-type NG, but emphasize that simiar cacuations shoud be performed for equiatera- and orthogona-type NG, as we as for NG estimators at the trispectrum eve. The most efficient method woud be to simutaneousy estimate a such foreground biases via simuations. III. PRIMORDIAL NON-GAUSSIANITY AND THE CMB BISPECTRUM In the oca mode of primordia NG, the primordia potentia Φ (where Φ 3 5ζ, and ζ is the adiabatic curvature perturbation) is given by [5, 49, 50] Φ( x) = Φ G ( x) + fnl oc ( Φ 2 G ( x) Φ 2 G ) +, (1) 3 where Φ G is a Gaussian fied and fnl oc is a constant characterizing the owest-order departure from Gaussianity. Muti-fied infationary modes, such as the curvaton mode, or non-infationary eary-universe scenarios, such as the ekpyrotic/cycic mode, can generate ocatype NG [51 56]. More generay, a detection of fnl oc 0 woud fasify singe-fied, sow-ro infation [8]. Current data are consistent with fnl oc = 0 [9, 12, 57, 58]. As discussed above, the tightest error bar (by a factor of a few) comes from the Panck CMB temperature data, thus motivating carefu scrutiny of this particuar observabe. It wi be possibe to further shrink the CMB-derived error bar on fnl oc by a factor of 2 3, but eventuay the cosmic variance (CV) imit wi be reached. Further improvements are then expected to come from arge-scae structure data (e.g., [59]). The non-inear couping in Eq. 1 generates a non-zero bispectrum in the CMB anisotropy. We consider ony the CMB temperature fied in the foowing. The anguar bispectrum B is defined via ( ) a 1m 1 a 2m 2 a 3m 3 = B , (2) m 1 m 2 m 3 where a m are the spherica harmonic ( coefficients ) of the CMB temperature fied and is the m 1 m 2 m 3 Wigner-3j symbo. Eq. 2 assumes ony rotationa invariance; if we additionay assume parity invariance (i.e., B = 0 if is odd, and thus B is invariant under a permutations), then we can define the reduced bispectrum b via B = (21 + 1)( )( ) ( ) 4π b (3) In the oca mode of primordia NG, the rea-space couping in Eq. 1 produces a non-zero Fourier-space bispectrum in the primordia potentia [5]: Φ( k 1 ) Φ( k 2 ) Φ( k 3 ) = 2(2π) 3 δ (3) ( k 1 + k 2 + k 3 )f oc NL P Φ (k 2 )P Φ (k 3 ) + 2 perm., (4) where Φ( k) is the Fourier transform of the primordia potentia, δ (3) is the 3D Dirac deta function, and P Φ (k) is the power spectrum of the primordia potentia. This non-zero bispectrum then yieds a non-zero anguar bispectrum in the CMB anisotropy, B oc The CMB temperature bispectrum can be computed straightforwardy via integras invoving the radiation transfer function, spherica Besse functions, and P Φ (k) (see, e.g., Ref. [5] for expicit formuae). The parameter fnl oc characterizes the ampitude of the bispectrum, as seen in Eq. 4. This bispectrum peaks in the squeezed imit, in which one

5 5 of the wavenumbers is much smaer than the other two (e.g., k 1 k 2, k 3 ) [2, 5]. For weak NG, assuming CV-imited CMB temperature data up to a mutipoe max, we can define the Fisher matrix eement for two bispectra B, B (e.g., [24]): F (B, B ) = 1 6 max B 123 B C 1 2 T T C T T C T T 3, (5) where C T T is the ensed primary CMB power spectrum (see Ref. [26] for usefu comments on impementing the sum in Eq. 5). The error bar on the ampitude of a bispectrum B is then given by the square root of the inverse of F (B, B). For exampe, the error bar on fnl oc (in the Gaussian approximation) is σ(fnl) oc = F 1 (B oc, B oc ). (6) Finay, the bias on the minimum-variance estimator for fnl oc due to a contaminating bispectrum Bcont is given by a ratio of Fisher matrix eements (e.g., [24]): IV. f oc,cont NL = F (Boc, B cont ) F (B oc, B oc ). (7) CMB LENSING-RELATED BIASES We first compute biases to fnl oc associated with the CMB ensing fied. These biases arise from bispectra with a form identica to that of the standard ensing-isw bispectrum [23, 24]: B Xφ = f C Xφ 2 C T T perm., (8) where X {ISW, tsz, CIB}, C Xφ is the cross-power spectrum of X and the CMB ensing potentia, and f is a couping kerne (e.g., [12, 24]): f = 1 2 [ 1( 1 + 1) + 2 ( 2 + 1) + 3 ( 3 + 1)] ( ) (21 + 1)( )( ) (9). 4π Note that B Xφ vanishes if X = ksz due to the equa probabiity of positive or negative LOS veocities (in genera, a bispectra invoving odd numbers of ksz fieds vanish for this reason). In Eq. 8, we use the ensed CMB power spectrum for C T T, as this yieds a more accurate resut for squeezed triange configurations than using the unensed spectrum [25]. The remaining quantity to be computed in Eq. 8 is the cross-power spectrum of X and φ. To owest order, the CMB ensing potentia is a weighted sum of the Newtonian potentia aong the LOS: φ(ˆn) = 2 χ c 2 dχ 0 ( χ χ χ χ ) Ψ(χˆn, χ), (10) where χ(z) is the comoving distance to redshift z, χ is the comoving distance to the surface of ast scattering at z = 1090, and Ψ is the 3D gravitationa potentia. In the foowing subsections, we detai our computation of C Xφ for X {ISW, tsz, CIB}, and cacuate the associated bias on measurements of f oc NL. A. Lensing-ISW Bias The ate-time ISW effect [20, 21] is generated by the decay of gravitationa potentias due to dark energy. In inear theory, this effect produces positive (negative) CMB temperature fuctuations as CMB photons traverse argescae overdensities (underdensities) in the ate-time matter fied. The fractiona temperature shift is given by the LOS integra of the time derivative of the gravitationa potentia: T ISW (ˆn) = 2 T CMB c 2 dt Ψ(ˆn), (11) LOS t where T CMB is the mean CMB temperature today. On sma scaes, noninear growth produces a ate-time ISW effect (the Rees-Sciama effect [21]) with the opposite sign to that sourced by dark energy on arge, inear scaes. However, this effect is much smaer in ampitude than the inear ISW effect due to dark energy. It has been shown that the change to the ensing-isw bispectrum due to the Rees-Sciama effect is essentiay undetectabe in the CMB, and that inear theory is sufficient for precisey computing the associated bias on fnl oc [26, 60]. Thus, we ony consider the inear ISW effect throughout this paper. We comment on instances where this may not suffice, particuary for the ISW-kSZ-kSZ bispectrum in V D. In the Limber approximation [61], the CMB ensing- ISW cross-power spectrum is [62 64] C ISW φ = 2 ( ) χ χ PΨ c 4 dz χ χ 3 z, (12) k=(+1/2)/χ where P Ψ (k, z) is the power spectrum of the 3D gravitationa potentia at wavenumber k and redshift z. Using the Poisson equation, we can express this resut in terms of the inear power spectrum of matter density fuctuations (defined at an arbitrary redshift), P in (k): C ISW φ = (1 + z) d dz 9Ω 2 mh 4 0 2c 4 ( + 1/2) 2 χ ( D(z) a(z) dz χ(z)(χ χ(z)) ) ( ) + 1/2 D(z)P in,(13) χ(z) where a(z) = 1/(1 + z) is the scae factor and D(z) is the inear growth factor, normaized in a manner consistent with the redshift at which P in is defined. We compute the inear matter power spectrum using CAMB. Note that during matter domination, D(z) a(z), and thus it can

6 6 2 ( + 1)C Xφ /(2π) [µk] ISW φ (inear theory) CIB φ (143 GHz) CIB φ (217 GHz) CIB φ (353 GHz) tsz φ (100 GHz) [ 1] tsz φ (143 GHz) [ 1] tsz φ (353 GHz) f NL oc ISW-φ bias (higher-precision Limber) ISW-φ bias (standard Limber) CV-imited error bar WMAP Panck max FIG. 1: Cross-power spectra of CMB ensing with other secondary anisotropy fieds, C Xφ, with X {ISW, tsz, CIB}. The ensing-isw cross-power spectrum (soid bue) is computed in inear theory, i.e., the Rees-Sciama effect is not incuded, but this has no measurabe effect on the associated bias on fnl. oc The ensing-tsz (dash-dotted) and ensing-cib cross-power spectra (dashed) are computed via the hao mode as described in IV B and IV C, respectivey. The hao mode cacuations are based on fits to measurements from Panck data in [66] (ensing-tsz) and [67] (ensing-cib). For carity, these cross-power spectra are ony shown for a subset of the Panck frequencies. Note that the tsz signa is negative (positive) at frequencies beow (above) 217 GHz, and vanishes at 217 GHz. be immediatey seen from Eq. 13 that there are no contributions to C ISW φ from this epoch, as expected. Fig. 1 shows the ensing-isw cross-power spectrum (soid bue curve). The signa fas off steepy with due to the ( + 1/2) 2 dependence in Eq. 13, which arises from the reation between the matter density and gravitationa potentia in the Poisson equation. Thus, in Eq. 8, the mutipoe associated with C ISW φ is generay the ong-waveength mode in squeezed triange configurations. Note that the signa is frequency-independent in backbody CMB temperature units. Fig. 2 shows the bias on fnl oc due to the ensing-isw bispectrum, computed via Eq. 7 for a fu-sky, CV-imited experiment up to a maximum mutipoe max. We show the resut cacuated using Eq. 13, as we as a cacuation using a ess-precise impementation of the Limber approximation in which a instances of + 1/2 on the right-hand side of Eq. 13 are repaced with [65]. This choice has a non-negigibe impact on the resuting bias on fnl oc. For Panck, the inferred bias for the higherprecision cacuation is fnl oc = 7.6, which agrees with the vaue used in the Panck 2015 NG anaysis [12]. For WMAP9, the bias is fnl oc = 2.3, sighty smaer than the vaue (2.6) quoted in Ref. [9] (their vaue is consistent with the standard Limber cacuation). Fig. 2 aso shows the 1σ uncertainty on fnl oc for a fu- FIG. 2: Bias on fnl oc from the ensing-isw bispectrum for a fu-sky, CV-imited experiment as a function of the maximum mutipoe incuded in the anaysis, max. This bias is backbody in frequency dependence and cannot be removed via component separation. The soid bue curve shows the bias when the ensing-isw cross-power spectrum is computed with a higher-precision impementation of the Limber approximation (i.e., the repacement + 1/2), whie the dashed orange curve shows the bias for a standard impementation of the Limber approximation. The higher-precision cacuation gives a bias fnl oc = 7.6 for Panck, consistent with the Panck NG anaysis [12]. The ight shaded region shows the 1σ uncertainty on fnl oc as a function of max using ony information in the CMB temperature bispectrum for a fu-sky, CV-imited experiment. The dashed vertica ines indicate the effective max for WMAP9 [9] and Panck 2015 [12]. sky, CV-imited experiment up to max, computed with Eq. 6. The WMAP9 error bar is σ(fnl oc ) = 19.9, whie the Panck 2015 error bar is σ(fnl oc ) = 5.7 (temperature data ony). Taking into account the sky masks used by WMAP9 [9] (f sky = 0.75) and Panck [12] (f sky = 0.76), we infer that max 485 for WMAP9 and max 1590 for Panck, which are potted as dashed vertica ines in Fig. 2. B. Lensing-tSZ Bias The tsz effect is generated by the inverse-compton scattering of CMB photons off hot, free eectrons, which are predominanty ocated in gaaxy groups and custers. Negecting reativistic corrections (e.g., [68]), the tsz signa is characterized by the Compton-y parameter, which is the LOS integra of the eectron pressure [29, 30]: y(ˆn) = σ T m e c 2 dχ a(χ) P e (χˆn, χ), (14) where σ T is the Thomson scattering cross-section, m e c 2 is the eectron rest-mass energy, and P e is the eectron

7 7 pressure. The CMB temperature fuctuation due to the tsz signa at a given frequency ν is then given by T tsz (ˆn) T CMB = g(ν)y(ˆn), (15) where g(ν) is the tsz spectra function: ( x ) g(ν) = x coth 4, (16) 2 with x hν/(k B T CMB ). The cross-correation between the tsz and CMB ensing fieds was first measured using Panck data in Ref. [66]. We adopt a mode consistent with this measurement in the foowing. The ensing-tsz cross-power spectrum, C yφ, can be computed straightforwardy in the hao mode (e.g., [69, 70]), anaogous to the computation of other tsz statistics (e.g., [71, 72]). We summarize the approach here, and refer the reader to Refs. [66, 72] for fu detais of these cacuations. The tota cross-power spectrum is the sum of the one-hao and two-hao terms: where C yφ,1h = and C yφ,2h = C yφ dz d2 V dzdω = C yφ,1h + C yφ,2h, (17) dn(m, z) dm dm ỹ(m, z) φ (M, z), ( ) dz d2 V + 1/2 dzdω P in χ(z), z dm 1 dn(m 1, z) dm 1 b(m 1, z)ỹ (M 1, z) (18) dm 2 dn(m 2, z) dm 2 b(m 2, z) φ (M 2, z). (19) d Here, 2 V dzdω is the comoving voume per steradian, dn/dm is the hao mass function (number of haos per unit mass per unit comoving voume), P in (k, z) D 2 (z)p in (k), b(m, z) is the inear hao bias, and ỹ (M, z) and φ (M, z) are the Fourier transform of the Compton-y and CMB ensing potentia profies, respectivey, of a hao of mass M at redshift z: ỹ (M, z) = σ T 4πr s,y m e c 2 2 s,y dx y x 2 sin(( + 1/2)x y / s,y ) y ( + 1/2)x y / s,y P e (x y r s,y, M, z) ; (20) 2 4πr s,φ φ (M, z) = ( + 1) 2 s,φ dx φ x 2 sin(( + 1/2)x φ / s,φ ) φ ( + 1/2)x φ / s,φ ρ(x φ r s,φ, M, z). (21) Σ crit (z) Here, r s,y is a characteristic scae radius of the eectron pressure profie, s,y = a(z)χ(z)/r s,y = d A (z)/r s,y is the mutipoe moment associated with this scae, and x y r/r s,y is a dimensioness radia variabe for the pressure profie. Anaogousy, r s,φ is a characteristic scae radius of the hao density profie ρ(r, M, z), s,φ = a(z)χ(z)/r s,φ = d A (z)/r s,φ is the mutipoe moment associated with this scae, and x φ r/r s,φ is a dimensioness radia variabe for the density profie. The quantity Σ crit (z) is the critica surface density for CMB ensing: Σ crit (z) = c 2 χ (1 + z) 4πGχ(z) (χ χ(z)). (22) We adopt the eectron pressure profie fitting function from the hydrodynamic simuations of [73, 74], the Navarro-Frenk-White (NFW) density profie [75], the concentration-mass reation of [76], and the fitting functions for the hao mass function and inear hao bias of [77] (updated from [78]). The concentration-mass reation is required in order to convert between mass definitions; we define M to be the viria mass foowing the definition of Ref. [79]. Further detais of this framework can be found in Refs. [66, 72]. The fiducia mass and redshift imits for a hao mode integras in this paper (e.g., Eqs. 18 and 19) are 10 5 M /h < M < M /h and < z < 8, respectivey. The ower redshift imit is imposed to avoid unphysica divergences at z = 0. We verify that a cacuations are converged with these choices. The ensing-tsz cross-power spectrum computed with this mode is shown in the dash-dotted curves in Fig. 1. The frequency dependence is evauated for the Panck channes using the bandpass-integrated tsz spectra function vaues provided in Ref. [40] (note that the signa is negative for ν < 217 GHz, positive for frequencies above this, and effectivey vanishes for the 217 GHz channe). The ensing-tsz cross-power spectrum is comparabe to the ensing-isw cross-power spectrum around 100, and is much arger at higher mutipoes. Fig. 3 shows the bias on fnl oc due to the ensing-tsz bispectrum, computed via Eqs. 17, 8, and 7. We show the bias for the Panck HFI channes from GHz, using the bandpass-integrated tsz spectra function vaues from Ref. [40]. For Panck, the bias is of order 1 2 for these channes, except for the 545 GHz channe, where it is somewhat arger. This concusion appears to agree with the resuts of Ref. [27] (by comparison to their Fig. 2, after removing the ensing- ISW bias). For an experiment with max = 3000, the fnl oc bias is comparabe to the 1σ error bar on fnl oc at the dominant CMB channes (100 and 143 GHz). Unike the ensing-isw bias, the ensing-tsz bias is non-backbody in frequency dependence, and can therefore be mitigated via mutifrequency component separation. In fact, it can be removed exacty using constrained component separation methods, in which the frequency channe weights are required to exacty nu the tsz spectra function in Eq. 16 [16]. However, such

8 8 C. Lensing-CIB Bias f NL oc WMAP 100 GHz (tsz-φ bias) 143 GHz (tsz-φ bias) 217 GHz (tsz-φ bias) Panck 353 GHz (tsz-φ bias) 545 GHz (tsz-φ bias) CV-imited error bar max FIG. 3: Bias on fnl oc from the ensing-tsz bispectrum as a function of max, as in Fig. 2. The bias is shown for the Panck HFI channes from GHz, with a frequency dependence arising from the tsz spectra function. For Panck, the bias takes vaues of fnl oc = 1.2 (100 GHz), fnl oc = 0.8 (143 GHz), and fnl oc = 1.8 (353 GHz). Due to its nonbackbody nature, this bias coud be removed via component separation, but recent evidence suggests that non-negigibe tsz signa has eaked into the Panck component-separated CMB maps [17, 18]. The ight shaded region and dashed vertica ines are identica to those in Fig. 2. constraints were not appied to the component-separated CMB temperature maps that were used in the Panck 2015 NG anaysis [12, 14]. Recent anayses have presented evidence that these maps have non-negigibe tsz contamination [17, 18]. A precise estimate of the tsz eakage as a function of anguar scae woud be needed to convert the frequency-dependent biases in Fig. 3 into a fina bias for Panck. A simper method woud be to perform the NG anaysis on a component-separated map in which the tsz signa has been nued, athough a statistica penaty in signa-to-noise (S/N) must be paid accordingy. Masking individuay detected gaaxy custers woud reduce the ensing-tsz bias to some extent, athough not by a arge amount, as C yφ is dominated by haos at ower masses and higher redshifts [66, 80] than are present in the Panck tsz cataog [81, 82]. Thus, even if such custers are masked in the Panck NG anaysis, it woud not strongy impact the biases shown in Fig. 3. Finay, note that C yφ has a fairy strong dependence on cosmoogica parameters, particuary σ 8 and Ω m [66], and thus the associated bias on fnl oc wi have a strong dependence as we. If we adopt the Panck 2015 CMB vaues for these parameters (σ 8 = and Ω m = 0.316) [11], the bias woud be 35% arger than shown in Fig. 3, assuming C yφ σ8ω m [66]. The CIB is sourced by the cumuative emission of dusty, star-forming gaaxies over cosmic time. The emission at different observationa frequencies is generated by gaaxies at somewhat different redshift ranges, but in genera the CIB redshift kerne has a broad peak around z 2, corresponding to the peak in the star formation rate density [67, 83, 84]. CIB statistics can be computed in the hao mode, anaogous to the tsz cacuations above, but with a more compicated prescription for the assignment of infrared fux to haos. For the ensing-cib cross-power spectrum, we simpy use the best-fit resuts of the Panck measurement of this quantity at each of the HFI frequencies [67], in ieu of impementing a detaied mode here. 4 This guarantees that our cacuations are consistent with actua measurements of C CIB φ. The cross-power spectra for three of the Panck HFI frequencies are shown as dashed curves in Fig. 1. The strong frequency dependence of the CIB emission is evident; the cross-power spectrum signa at 353 GHz is neary an order of magnitude arger than at 217 GHz. Fig. 4 shows the bias on fnl oc due to the ensing-cib bispectrum, computed using the C CIB φ fits from Ref. [67] in combination with Eqs. 8 and 7. We show the bias for the Panck HFI channes from GHz. The bias is strongy frequency-dependent, a direct resut of the strong frequency dependence shown in Fig. 1. However, unike the ensing-tsz bias, the ensing-cib bias has the same sign at a frequencies (at fixed max ). At 217 GHz, the bias for Panck is fnl oc = 4.3, whie at 353 GHz, it is neary an order of magnitude arger. These resuts are simiar to those presented in Ref. [28] (see their Tabe 3), but up to a factor of 2 3 arger than theirs at some Panck frequencies. In this context, we note that our cacuation of C CIB φ is directy drawn from fits to Panck measurements as described above, rather than a theoretica mode. For an experiment with max = 3000, the ensing-cib bias is comparabe to the 1σ error bar on fnl oc at 100 or 143 GHz, and is much arger than this at higher frequencies. Like the ensing-tsz bias, the ensing-cib bias is nonbackbody in frequency dependence, and can therefore be reduced via mutifrequency component separation. However, unike the ensing-tsz bias, it cannot be fuy eiminated, as the CIB decorreates across frequency channes to some extent (because the redshift kerne of the emission is different at different frequencies) [84, 85]. Given the evidence of tsz eakage into the Panck componentseparated CMB maps, it is pausibe that non-negigibe CIB eakage is aso present, athough the extent of such contamination is presenty uncear. A dedicated 4 We thank Oivier Doré for providing these fits.

9 9 f NL oc WMAP 100 GHz (CIB-φ bias) 143 GHz (CIB-φ bias) 217 GHz (CIB-φ bias) Panck 353 GHz (CIB-φ bias) / GHz (CIB-φ bias) / 80 CV-imited error bar max FIG. 4: Bias on fnl oc from the ensing-cib bispectrum as a function of max, as in Fig. 2. The bias is shown for the Panck HFI channes from GHz, with a strong frequency dependence due to the increase in dust emission intensity at high frequencies. Note that the 353 GHz and 545 GHz resuts are divided by factors of 8 and 80, respectivey, to reduce the dynamic range of the pot. For Panck, the bias takes vaues of fnl oc = 1.4 (143 GHz), fnl oc = 4.3 (217 GHz), and fnl oc = 38 (353 GHz). This bias is non-backbody in frequency dependence, and is thus reduced by mutifrequency component separation techniques (however, it cannot be fuy eiminated due to CIB decorreation). The extent of this reduction for the Panck NG anaysis is currenty uncear. The ight shaded region and dashed vertica ines are identica to those in Fig. 2. simuation-based anaysis is necessary to quantify the tota bias on fnl oc resuting from the ensing-cib crosscorreation. V. ISW-RELATED BIASES strong contributions in the squeezed imit, as the ISW signa peaks on arge scaes (i.e., comprising the ongwaveength mode of the triange) whie the other fieds peak on sma scaes (i.e., comprising the two shortwaveength modes of the triange). A physica interpretation of these bispectra is that the ong-waveength ISW fied is moduating the ampitude of the short-waveength power spectra of the other fieds (e.g., the sma-scae tsz power spectrum). A ong-waveength overdensity (underdensity) producing a positive (negative) ISW fuctuation wi aso contain more (fewer) massive haos, thereby corresponding to a higher (ower) ampitude of the smascae tsz/cib/ksz power spectra. We compute these bispectra in the hao mode, working throughout in the Limber approximation for bispectra [86, 87]. In genera, the fu bispectrum wi contain three-hao, two-hao, and one-hao terms, but here we focus ony on the contributions that are expected to dominate in squeezed configurations reevant to fnl oc. The primary such contribution arises from a two- hao term, in which one mutipoe corresponds to a ong-waveength ISW fuctuation (hence, we use the term hao oosey here) and the other two mutipoes correspond to short-waveength tsz/cib/ksz fuctuations. Note that throughout we consider ony the inear-theory ISW effect, i.e., the Rees-Sciama effect is negected. Thus, there is effectivey no one-hao contribution to these bispectra (noninear growth woud generate a Rees-Sciama one-hao term, but this is much smaer than the inear-theory ISW signa). We aso negect the three-hao term, which is sourced by the tree-eve bispectrum, as it is not expected to contribute strongy to squeezed configurations. One exception to this may be the three-hao contribution to the ISW-kSZ-kSZ bispectrum (for which we aso negect other potentiay important contributions see V D). We defer a fu cacuation to future work, and focus ony on the two-hao contributions in the foowing. The ISW effect is a tracer of the ate-time gravitationa potentia, and therefore it is correated with other such tracers that source secondary CMB anisotropies. For our purposes here, we are primariy concerned with bispectra invoving one ISW eg and two egs drawn from the tsz, CIB, or ksz fieds. These bispectra have A. ISW-tSZ-tSZ Bias For CMB temperature maps at three frequencies ν 1, ν 2, and ν 3, the two- hao contribution to the ISW-tSZtSZ reduced bispectrum is: b (Tν 1 Tν 2 Tν 3 ),2h 1 2 3,(ISW tsz tsz) = g(ν 2 )g(ν 3 ) + P in ( χ(z) ( P in χ(z) + 2 perm., dz I 1 (z) ) ) [ ( ) 2 P in χ(z) dm dn dm b(m, z)ỹ 2 (M, z) dm dn dm b(m, z)ỹ 3 (M, z) dm dn dm b(m, z)ỹ 2 (M, z)ỹ 3 (M, z) dm dn dm b(m, z)ỹ 3 (M, z) M ρ m dm dn dm b(m, z)ỹ 2 (M, z) M ρ m ] (23)

10 10 where I (z) = 3Ω mh 2 0 c 2 ( + 1 2) 2 χ 2 (z)d(z) d dz ( ) D(z), (24) a(z) and the additiona permutations correspond to cases in which the ISW mutipoe is either 2 or 3. Here, we have made the approximation that the contribution from the interna structure of the hao density profie can be negected in the ISW factors, i.e., k1 ISW 0, so that the Fourier transform of the density profie simpy yieds a factor of M. This approximation is accurate due to the rapid decine of the ISW signa as increases. Stated differenty, the ISW signa is effectivey sourced ony by inear modes of the density fied. For computationa efficiency, we set the ISW signa to zero above ISW,cut = 200 in a of the foowing cacuations. We verify that our resuts are converged with this choice, i.e., higher vaues of ISW,cut do not change the derived bias on fnl oc for any of these bispectra. Finay, note that the first term in Eq. 23, in which both tsz mutipoes beong to the same hao, dominates by a factor of over the atter two terms, in which the tsz mutipoes are in two distinct haos, except for configurations in which a three mutipoes are very sma (i.e., on very arge scaes). This is precisey anaogous to the dominance of the one-hao term over the two-hao term in the tsz power spectrum for a 10 [71, 72, 88]. Nevertheess, we incude a terms in the foowing cacuations. We use the same modes to compute Eq. 23 as used for the tsz cacuations described in IV B. 5 Fig. 5 shows a sice through the ISW-tSZ-tSZ bispectrum (considering ony the contributions in Eq. 23), as we as the oca bispectrum (with fnl oc = 1) and the ensing-isw bispectrum. The atter two bispectra dispay acoustic osciations arising from the radiation transfer functions, whereas the ISW-tSZ-tSZ bispectrum is smooth (the primary temperature power spectrum does not appear in Eq. 23, in contrast to Eq. 8). For this particuar sice, the ISW-tSZ-tSZ bispectrum becomes comparabe in ampitude to the oca bispectrum at 1500, and is much arger at higher mutipoes. This pot aso iustrates the origin of the osciatory behavior seen for the ensing-isw bias on fnl oc in Fig. 2 (and to some extent in Figs. 3 and 4): the ensing-isw and oca bispectra have osciations that are not exacty in phase, and the ensing-isw bispectrum furthermore osciates between positive and negative vaues. These effects ead to osciations in the inner product in the numerator of Eq. 7. In contrast, the smooth shape of the ISW-tSZ-tSZ bispectrum suggests that the associated bias on fnl oc wi be a smoothy increasing function of max, which indeed is the case (see Fig. 6). Finay, whie Fig. 5 shows the ISW-tSZ-tSZ bispectrum at an 5 As a cross-check, we aso compute the ISW-tSZ cross-power spectrum and obtain resuts in genera agreement with those of Ref. [89]. ( 1) 4 B10,,+10 [µk 3 ] ISW-tSZ-tSZ (148 GHz) oca (fnl = 1) 1 ISW-φ FIG. 5: A sice through the ISW-tSZ-tSZ bispectrum at 148 GHz computed via Eq. 23 (bue), the oca primordia bispectrum with fnl oc = 1 (orange), and the ensing-isw bispectrum (green), for squeezed configurations with 1 = 10, 2, and 3 = In contrast to the atter two bispectra, the ISW-tSZ-tSZ bispectrum is a smoothy increasing function of. expected. Note that it peaks in the squeezed imit, as exampe frequency of 148 GHz, this bispectrum is aways positive when evauated at a singe frequency, due to the quadratic tsz spectra function dependence (when evauated for a set of different frequencies, it coud be negative or positive, but never crosses zero). Fig. 6 shows the bias on fnl oc due to the ISW-tSZ-tSZ bispectrum, computed via Eqs. 23 and 7. We show the bias for the Panck HFI channes from GHz, using the bandpass-integrated tsz spectra function vaues from Ref. [40]. We do not pot any cross-frequency biases (i.e., invoving different vaues of ν 2 and ν 3 in Eq. 23), athough these are present and can be of positive or negative sign. In contrast, the singe-frequency biases are aways negative, as shown in Fig. 6. For Panck, the bias is fnl oc oc = 4.5 (100 GHz), fnl = 2.1 (143 GHz), and fnl oc = 11 (353 GHz). If a non-negigibe fraction of tsz signa has eaked into the component-separated CMB maps used in the Panck NG anaysis, the ISW-tSZ-tSZ bias coud thus yied a shift of order 1σ in the inferred vaue of fnl oc. For an experiment with max = 3000, the bias is many times arger than the 1σ error bar on fnl oc at a of the frequencies considered (except 217 GHz, where the tsz nu occurs), incuding the dominant CMB channes (100 and 143 GHz). As discussed in IV B for the ensing-tsz bias, the ISWtSZ-tSZ bias is non-backbody in nature, and in fact can be exacty removed via constrained component separation techniques (at a cost in S/N). It can aso be modeed and subtracted, athough modeing it through the process of component separation ikey requires fu simuations. The bias wi aso be reduced to some extent

11 11 f NL oc WMAP Panck 100 GHz (ISW-tSZ-tSZ bias) 143 GHz (ISW-tSZ-tSZ bias) 217 GHz (ISW-tSZ-tSZ bias) 353 GHz (ISW-tSZ-tSZ bias) 545 GHz (ISW-tSZ-tSZ bias) CV-imited error bar max FIG. 6: Bias on fnl oc from the ISW-tSZ-tSZ bispectrum as a function of max, as in Fig. 2. The bias is shown for the Panck HFI channes from GHz, with a frequency dependence arising from the tsz spectra function. For Panck, the bias takes vaues of fnl oc = 4.5 (100 GHz), fnl oc = 2.1 (143 GHz), and fnl oc = 11 (353 GHz). Due to its nonbackbody nature, this bias coud be removed via component separation, but recent evidence suggests that non-negigibe tsz signa has eaked into the Panck component-separated CMB maps [17, 18]. The ight shaded region and dashed vertica ines are identica to those in Fig. 2. by masking known gaaxy custers, but since it is effectivey sourced by the sma-scae tsz power spectrum, the bias is mosty generated by custers that are beow the mass threshod for individua detection in Panck (see, e.g., Refs. [72, 74] for breakdowns of the hao mass and redshift contributions to the tsz power spectrum). Amongst the biases considered in this paper, the ISWtSZ-tSZ bias is the most sensitive to cosmoogica parameters, as it inherits the strong dependence of the tsz power spectrum on σ 8 and Ω m (e.g., [71, 72]). If we adopt the Panck 2015 CMB vaues for these parameters [11], the bias woud be 70% arger than shown in Fig. 6, assuming that the bispectrum in Eq. 23 foows the tsz power spectrum parameter dependence: C yy σ8ω 8 3 m [72, 90]. In addition, it is sensitive to the modeing of astrophysica processes in the intracuster medium. Given these parameter and modeing dependences, it seems best to simpy nu the tsz spectra function via constrained component separation when constructing CMB maps for NG anaysis. B. ISW-CIB-CIB Bias The two-hao contribution to the ISW-CIB-CIB reduced bispectrum is identica to Eq. 23, with the repacement g(ν)ỹ (M, z) ĨCIB,ν(1+z) (M, z), i.e., the tsz profie of each hao is repaced by its redshifted infrared emission profie. 6 In ieu of impementing a detaied mode for the infrared emission of each hao, here we adopt a simper, approximate approach reying on the very high correation coefficient between the CIB fied (at GHz) and the CMB ensing fied [67, 91, 92]. If the correation coefficient between these fieds were unity at a mutipoes, a CIB map woud simpy be a rescaed version of a CMB ensing map, with an -dependent rescaing factor f. Empiricay, the correation coefficient is 80% at a mutipoes up to 2000 for most Panck HFI frequencies [67]. For the purpose of approximatey estimating biases to fnl oc, we consider this sufficienty cose to unity to simpy approximate ĨCIB,ν(1+z) (M, z) f (ν) φ (M, z), where φ (M, z) is given by Eq. 21. We determine the rescaing factor f (ν) for a given frequency channe using the ensing-cib cross-power spectrum resuts of Ref. [67], i.e., the same fits to C CIB φ used in IV C, in combination with a theoretica cacuation of the CMB ensing auto-power spectrum: f (ν) = CCIBν φ C φφ. (25) Thus, we approximate the ISW-CIB-CIB bispectrum by computing Eq. 23 with g(ν)ỹ (M, z) f (ν) φ (M, z), for both ν 2 and ν 3. Ceary this approach wi not yied percent-eve accuracy, but it suffices to assess the order of magnitude of the bias on fnl oc. Fig. 7 shows the bias on fnl oc due to the ISW-CIB-CIB bispectrum, computed via Eqs. 23, 25, and 7. Due to the approximate nature of this cacuation, we ony show resuts for three of the Panck HFI channes (143, 217, and 353 GHz). We do not pot any cross-frequency biases (i.e., invoving different vaues of ν 2 and ν 3 in Eq. 23), athough these aso exist. In contrast to the ISW-tSZ-tSZ bias, the ISW-CIB-CIB bias is aways negative, since the CIB signa is positive at a frequencies. For Panck, the (approximate) ISW-CIB-CIB bias is fnl oc = 0.4 (217 GHz) and fnl oc = 33 (353 GHz). The steep frequency dependence of this bias is expected due to the strong frequency dependence of the CIB. The resut is thus very sensitive to an assessment of the degree to which CIB emission has eaked into the component-separated Panck CMB maps. For an experiment with max = 3000, the ISW-CIB-CIB bias is many times arger than the 1σ error bar on fnl oc, even at 217 GHz. C. ISW-tSZ-CIB Bias Using the modes described in the previous two subsections, we can readiy compute the ISW-tSZ-CIB bispectrum via Eq. 23. We simpy repace ony one of the tsz 6 Additiona shot noise terms may aso be incuded, depending on the detais of the underying CIB hao mode.