We study the position-dependent power spectrum and the integrated bispectrum statistic for 2D cosmological fields on the sphere (integrated angular bispectrum). First, we derive a useful, $m$-independent, formula for the full-sky integrated angular bispectrum, based on the construction of azimuthally symmetric patches. We then implement a pipeline for integrated angular bispectrum estimation, including a mean-field correction to account for spurious isotropy-breaking effects in realistic conditions (e.g., inhomogenous noise, sky masking). Finally, we show examples of applications of this estimator to CMB analysis, both using simulations and actual Planck data. Such examples include $f_\mathrm{NL}$ estimation, analyses of non-Gaussianity from secondary anisotropies (ISW-lensing and ISW-tSZ-tSZ bispectra) and studies of non-Gaussian signatures from foreground contamination.

We analyse the Planck full-mission cosmic microwave background (CMB) temperature and E-mode polarization maps to obtain constraints on primordial non-Gaussianity (NG). We compare estimates obtained from separable template-fitting, binned, and modal bispectrum estimators, finding consistent values for the local, equilateral, and orthogonal bispectrum amplitudes. Our combined temperature and polarization analysis produces the following results: f_NL^local = -0.9 +\- 5.1; f_NL^equil = -26 +\- 47; and f_NL^ortho = - 38 +\- 24 (68%CL, statistical). These results include the low-multipole (4 <= l < 40) polarization data, not included in our previous analysis, pass an extensive battery of tests, and are stable with respect to our 2015 measurements. Polarization bispectra display a significant improvement in robustness; they can now be used independently to set NG constraints. We consider a large number of additional cases, e.g. scale-dependent feature and resonance bispectra, isocurvature primordial NG, and parity-breaking models, where we also place tight constraints but do not detect any signal. The non-primordial lensing bispectrum is detected with an improved significance compared to 2015, excluding the null hypothesis at 3.5 sigma. We present model-independent reconstructions and analyses of the CMB bispectrum. Our final constraint on the local trispectrum shape is g_NLl^local = (-5.8 +\-6.5) x 10^4 (68%CL, statistical), while constraints for other trispectra are also determined. We constrain the parameter space of different early-Universe scenarios, including general single-field models of inflation, multi-field and axion field parity-breaking models. Our results provide a high-precision test for structure-formation scenarios, in complete agreement with the basic picture of the LambdaCDM cosmology regarding the statistics of the initial conditions (abridged).

We use the binned bispectrum estimator to determine the bispectra of the dust, free-free, synchrotron, and AME galactic foregrounds using maps produced by the Commander component separation method from Planck 2015 data. We find that all of these peak in the squeezed configuration, allowing for potential confusion with in particular the local primordial shape. Applying an additional functionality implemented in the binned bispectrum estimator code, we then use these galactic bispectra as templates in an $f_\mathrm{NL}$ analysis of other maps. After testing and validating the method and code with simulations, we show that we detect the dust in the raw 143 GHz map with the expected amplitude (the other galactic foregrounds are too weak at 143 GHz to be detected) and that no galactic residuals are detected in the cleaned CMB map. We also investigate the effect of the mask on the templates and the effect of the choice of binning on a joint dust-primordial $f_\mathrm{NL}$ analysis.

We use the long-wavelength formalism to investigate the level of bispectral non-Gaussianity produced in two-field inflation models with standard kinetic terms. Even though the Planck satellite has so far not detected any primordial non-Gaussianity, it has tightened the constraints significantly, and it is important to better understand what regions of inflation model space have been ruled out, as well as prepare for the next generation of experiments that might reach the important milestone of Delta f_NL(local) = 1. We derive an alternative formulation of the previously derived integral expression for f_NL, which makes it easier to physically interpret the result and see which types of potentials can produce large non-Gaussianity. We apply this to the case of a sum potential and show that it is very difficult to satisfy simultaneously the conditions for a large f_NL and the observational constraints on the spectral index n_s. In the case of the sum of two monomial potentials and a constant we explicitly show in which small region of parameter space this is possible, and we show how to construct such a model. Finally, the new general expression for f_NL also allows us to prove that for the sum potential the explicit expressions derived within the slow-roll approximation remain valid even when the slow-roll approximation is broken during the turn of the field trajectory (as long as only the epsilon slow-roll parameter remains small).

The spectrum of the diffuse isotropic component of cosmic X-rays over the 13-180 keV range was determined by the UCSD/MIT Hard X-Ray and Gamma-Ray instrument (HEAO A4) on the High Energy Astronomical Observatory-1 (HEAO-1). The instrument consists of a complex of actively shielded and collimated scintillation counters, including the Low Energy Detector set from which the data reported here were obtained. These data join smoothly with the spectrum at lower energies reported by the GSFC HEAO A2 instrument and with that measured to 400 keV by the HEAO A4 Medium Energy Detectors. The HEAO data set also joins the recent results from COMPTEL on the Compton Gamma-Ray Observatory in the 1-10 MeV range, which failed to confirm the existence of an "MeV bump" in this range. Although the spectrum over the entire range 3 keV < E < 100 GeV can be fit by a simple empirical analytic expression, the origin is likely due to a number of distinct source components. The prevailing idea for the origin is that the hard X-ray spectrum is due to X-rays from various AGN components, particularly Seyfert galaxies extending to cosmological distances, and that the low energy gamma-rays may be due to emission from type 1a Supernovae, also integrated to cosmological distances. The higher energy gamma-ray spectrum defined by EGRET, also on the CGRO, may be due to unresolved gamma-ray emitting blazars. Models of production by these source components, extrapolated to the present epoch, must reproduce the observationally derived spectrum.

We present the results of OSSE observations of the soft gamma ray continuum emission from the Galactic plane at longitude 95 degrees. Emission is detected between 50 and 600 keV where the spectrum is fit well by a power law with photon index -2.6+-0.3 and flux (4.0+-0.5) 10^{-2} photons/s/cm^2/rad/MeV at 100 keV. This spectral shape in this range is similar to that found for the continuum emission from the inner Galaxy but the amplitude is lower by a factor of four. This emission is either due to unresolved and previously unknown point sources or it is of diffuse origin, or a combination of the two. Simultaneous observations with OSSE and smaller field of view instruments operating in the soft gamma ray energy band, such as XTE or SAX, would help resolve this issue. If it is primarily diffuse emission due to nonthermal electron bremsstrahlung, as is the >1 MeV Galactic ridge continuum, then the power in low energy cosmic ray electrons exceeds that of the nuclear component of the cosmic rays by an order of magnitude. This would have profound implications for the origin of cosmic rays and the energetics of the interstellar medium. Alternatively, if the emission is diffuse and thermal, then there must be a component of the interstellar medium at temperatures near 10^9 K.