The ISW-Lensing Bispectrum & Trispectrum

Due to the integrated Sachs-Wolfe (ISW) effect, cosmic microwave background (CMB) temperature and polarization fluctuations are correlated with the gravitational lensing potential. Famously, this induces a CMB three-point function, whose shape can be used to constrain dark energy and modifications to gravity. An analogous effect occurs at higher-order, producing an ISW-lensing trispectrum whose amplitude is hitherto unconstrained. We present a detailed discussion of this effect, and define minimum-variance estimators for the ISW-lensing three- and four-point functions. These are implemented within the PolySpec code, and bear strong similarities to the quadratic estimators used in lensing analyses. Applying these tools to Planck, we obtain strong detections of the bispectrum amplitude (consistent with previous works), but find only weak constraints on the trispectrum, due to a strong cancellation between the various ISW-induced contributions. We additionally forecast the constraints from future datasets, finding that (a) simple estimators for the ISW-lensing bispectrum will be severely limited by non-Gaussian modifications to the covariance, and (b) the ISW-lensing trispectrum will be very challenging to detect even with high-resolution future experiments. We finally consider the induced bias on primordial non-Gaussianity amplitudes (and lensing itself), which we show to be large for the bispectrum (as expected) but negligible for the trispectrum.