We report the discovery of X-ray emission from CFHQS J142952+544717, the most distant known radio-loud quasar at z=6.18, on Dec. 10--11, 2019 with the eROSITA telescope on board the SRG satellite during its ongoing all-sky survey. The object was identified by cross-matching an intermediate SRG/eROSITA source catalog with the Pan-STARRS1 distant quasar sample at 5.6 < z < 6.7. The measured flux $\sim 8 \times 10^{-14}$ erg cm$^{-2}$ s$^{-1}$ in the 0.3--2 keV energy band corresponds to an X-ray luminosity of $2.6^{+1.7}_{-1.0}\times 10^{46}$ erg s$^{-1}$ in the 2--10 keV rest-frame energy band, which renders CFHQS J142952+544717 the most X-ray luminous quasar ever observed at z > 6. Combining our X-ray measurements with archival and new photometric measurements in other wavebands (radio to optical), we estimate the bolometric luminosity of this quasar at $\sim (2$--$3) \times 10^{47}$ erg s$^{-1}$. Assuming Eddington limited accretion and isotropic emission, we infer a lower limit on the mass of the supermassive black hole of $\sim 2\times 10^9 M_\odot$. The most salient feature of CFHQS J142952+544717 is its X-ray brightness relative to the optical/UV emission. We argue that it may be linked to its radio-loudness (although the object is not a blazar according to its radio properties), specifically to a contribution of inverse Compton scattering of cosmic microwave background photons off relativistic electrons in the jets. If so, CFHQS J142952+544717 might be the tip of the iceberg of high-z quasars with enhanced X-ray emission, and SRG/eROSITA may find many more such objects during its 4 year all-sky survey.

We report on INTEGRAL observations of the soft $\gamma$-ray repeater SGR 1935+2154 performed between 2020 April 28 and May 3. Several short bursts with fluence of $\sim10^{-7}-10^{-6}$ erg cm$^{-2}$ were detected by the IBIS instrument in the 20-200 keV range. The burst with the hardest spectrum, discovered and localized in real time by the INTEGRAL Burst Alert System, was spatially and temporally coincident with a short and very bright radio burst detected by the CHIME and STARE2 radio telescopes at 400-800 MHz and 1.4 GHz, respectively. Its lightcurve shows three narrow peaks separated by $\sim$29 ms time intervals, superimposed on a broad pulse lasting $\sim$0.6 s. The brightest peak had a delay of 6.5$\pm$1.0 ms with respect to the 1.4 GHz radio pulse (that coincides with the second and brightest component seen at lower frequencies). The burst spectrum, an exponentially cut-off power law with photon index $\Gamma=0.7_{-0.2}^{+0.4}$ and peak energy $E_p=65\pm5$ keV, is harder than those of the bursts usually observed from this and other magnetars. By the analysis of an expanding dust scattering ring seen in X-rays with the {\it Neil Gehrels Swift Observatory} XRT instrument, we derived a distance of 4.4$_{-1.3}^{+2.8}$ kpc for SGR 1935+2154, independent of its possible association with the supernova remnant G57.2+0.8. At this distance, the burst 20-200 keV fluence of $(6.1\pm 0.3)\times10^{-7}$ erg cm$^{-2}$ corresponds to an isotropic emitted energy of $\sim1.4\times10^{39}$ erg. This is the first burst with a radio counterpart observed from a soft $\gamma$-ray repeater and it strongly supports models based on magnetars that have been proposed for extragalactic fast radio bursts.

We take advantage of a set of molecular cloud simulations to demonstrate a possibility to uncover statistical properties of the gas density and velocity fields using reflected emission of a short (with duration much less than the cloud's light-crossing time) X-ray flare. Such situation is relevant for the Central Molecular Zone of our Galaxy where several clouds get illuminated by a $\sim110$ yr-old flare from the supermassive black hole Sgr A*. Due to shortness of the flare ($\Delta t\lesssim1.6$ yrs), only a thin slice ($\Delta z\lesssim0.5$ pc) of the molecular gas contributes to the X-ray reflection signal at any given moment, and its surface brightness effectively probes the local gas density. This allows reconstructing the density probability distribution function over a broad range of scales with virtually no influence of attenuation, chemo-dynamical biases and projection effects. Such measurement is key to understanding the structure and star-formation potential of the clouds evolving under extreme conditions in the CMZ. For cloud parameters similar to the currently brightest in X-ray reflection molecular complex Sgr A, the sensitivity level of the best available data is sufficient only for marginal distinction between solenoidal and compressive forcing of turbulence. Future-generation X-ray observatories with large effective area and high spectral resolution will dramatically improve on that by minimising systematic uncertainties due to contaminating signals. Furthermore, measurement of the iron fluorescent line centroid with sub-eV accuracy in combination with the data on molecular line emission will allow direct investigation of the gas velocity field.

The Massive and Distant Clusters of WISE Survey (MaDCoWS) provides a catalog of high-redshift ($0.7\lesssim z\lesssim 1.5$) infrared-selected galaxy clusters. However, the verification of the ionized intracluster medium, indicative of a collapsed and nearly virialized system, is made challenging by the high redshifts of the sample members. The main goal of this work is to test the capabilities of the Atacama Compact Array (ACA; also known as the Morita Array) Band 3 observations, centered at about 97.5 GHz, to provide robust validation of cluster detections via the thermal Sunyaev-Zeldovich (SZ) effect. Using a pilot sample that comprises ten MaDCoWS galaxy clusters, accessible to ACA and representative of the median sample richness, we infer the masses of the selected galaxy clusters and respective detection significance by means of a Bayesian analysis of the interferometric data. Our test of the "Verification with the ACA - Localization and Cluster Analysis" (VACA LoCA) program demonstrates that the ACA can robustly confirm the presence of the virialized intracluster medium in galaxy clusters previously identified in full-sky surveys. In particular, we obtain a significant detection of the SZ effect for seven out of the ten VACA LoCA clusters. We note that this result is independent of the assumed pressure profile. However, the limited angular dynamic range of the ACA in Band 3 alone, short observational integration times, and possible contamination from unresolved sources limit the detailed characterization of the cluster properties and the inference of the cluster masses within scales appropriate for the robust calibration of mass-richness scaling relations.

We show that Compton scattering by electrons of the hot intergalactic gas in galaxy clusters should lead to peculiar distortions of the cosmic background X-ray and soft gamma-ray radiation - an increase in its brightness at E<60-100 keV and a drop at higher energies. The distortions allow the most important cluster parameters to be measured. The spectral shape of the distortions and its dependence on the gas temperature, optical depth, and surface density distribution law have been studied using Monte Carlo computations and confirmed by analytical estimations. In the cluster frame the maximum of the background decrease due to the recoil effect occurs at ~500-600 keV. The photoionization of H- and He-like iron and nickel ions leads to additional distortions in the background spectrum - a strong absorption line with the threshold at ~9 keV (and also to an absorption jump at ~2 keV for cold clusters). The absorption of intrinsic thermal radiation from the cluster gas by these ions also leads to such lines. In nearby (z<1) clusters the line at ~2 keV is noticeably enhanced by absorption in the colder (~10^6 K) plasma of their peripheral (~3 Mpc) regions; moreover, the absorption line at ~1.3 keV splits off from it. The redshift of distant clusters shifts the absorption lines in the background spectrum (at ~2, ~9, and ~500 keV) to lower energies. Thus, in contrast to the microwave background scattering effect, this effect depends on the cluster redshift z, but in a very peculiar way. When observing clusters at z>1, the effect allows one to determine how the X-ray background evolved and how it was "gathered" with z. To detect the effect, the accuracy of measurements should reach ~0.1%. We consider the most promising clusters for observing the effect and discuss the techniques whereby the influence of the thermal gas radiation hindering the detection of background distortions should be minimal.

A neutral beam injection technique is employed in all major TOKAMAK facilities for heating of magnetically confined plasma. The question then arises, whether a similar mechanism might work in astrophysical objects? For instance, a hyper-Eddington galactic binary SS433 possesses baryonic jets, moving at a quarter of the speed of light, and observations revealed signs of gas cooling and recombination on sub-pc scales and equally strong signs of powerful energy deposition on much larger scales $\sim$100 pc. Here we consider a model where neutral atoms transport this energy. A sub-relativistic beam of neutral atoms penetrates the interstellar medium, these atoms gradually get ionised and deposit their energy over a region, whose longitudinal dimension is set by the "ionisation length". The channel, where the energy is deposited, expands sideways and drives a shock in the lateral direction. Once the density in the channel drops, the heating rate by the beam drops accordingly, and the region of the energy release moves along the direction of the beam. We discuss distinct features associated with this scenario and speculate that such configuration might also boost shock acceleration of the "pick-up" protons that arise due to ionisation of neutral atoms both upstream and downstream of the shock.

At the age of about 1 year, the spectra of most Type Ia supernovae are dominated by strong forbidden nebular emission lines of FeII and FeIII. Later observations (at about 2 years) of the nearby SN 2011fe showed an unexpected shift of ionization to FeI and FeII. Spectra of the very nearby SN Ia 2014J at an intermediate phase (1 - 1.5 years) that are presented here show a progressive decline of FeIII emission, while FeI is not yet strong. The decrease in ionization can be explained if the degree of clumping in the ejecta increases significantly at ~1.5 years, at least in the Fe-dominated zone. Models suggest that clumps remain coherent after about one year, behaving like shrapnel. The high density in the clumps, combined with the decreasing heating rate, would then cause recombination. These data may witness the phase of transition from relatively smooth ejecta to the very clumpy morphology that is typical of SN remnants. The origin of the increased clumping may be the development of local magnetic fields.

Reflection of X-ray emission on molecular clouds in the inner $\sim$ 100 pc of our Galaxy reveals that, despite being extremely quiet at the moment, our supermassive black hole Sgr A* should have experienced bright flares of X-ray emission in the recent past. Thanks to the improving characterisation of the reflection signal, we are able to infer parameters of the most recent flare(s) (age, duration and luminosity) and relative line-of-sight disposition of the brightest individual molecular complexes. We show that combining these data with measurements of polarisation in the reflected X-ray continuum will not only justify Sgr A* as the primary source but also allow deriving intrinsic polarisation properties of the flare emission. This will help to identify radiation mechanisms and underlying astrophysical phenomena behind them. For the currently brightest reflecting molecular complex, Sgr A, the required level of sensitivity might be already accessible with upcoming X-ray polarimeters.

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first (and significant) detection of this signal in the tSZ maps, and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales.

The European Space Agency's Planck satellite, which was dedicated to studying the early Universe and its subsequent evolution, was launched on 14 May 2009. It scanned the microwave and submillimetre sky continuously between 12 August 2009 and 23 October 2013, producing deep, high-resolution, all-sky maps in nine frequency bands from 30 to 857GHz. This paper presents the cosmological legacy of Planck, which currently provides our strongest constraints on the parameters of the standard cosmological model and some of the tightest limits available on deviations from that model. The 6-parameter LCDM model continues to provide an excellent fit to the cosmic microwave background data at high and low redshift, describing the cosmological information in over a billion map pixels with just six parameters. With 18 peaks in the temperature and polarization angular power spectra constrained well, Planck measures five of the six parameters to better than 1% (simultaneously), with the best-determined parameter (theta_*) now known to 0.03%. We describe the multi-component sky as seen by Planck, the success of the LCDM model, and the connection to lower-redshift probes of structure formation. We also give a comprehensive summary of the major changes introduced in this 2018 release. The Planck data, alone and in combination with other probes, provide stringent constraints on our models of the early Universe and the large-scale structure within which all astrophysical objects form and evolve. We discuss some lessons learned from the Planck mission, and highlight areas ripe for further experimental advances.

The thermal Sunyaev-Zeldovich (SZ) effect presents a relatively new tool for characterizing galaxy cluster merger shocks, traditionally studied through X-ray observations. Widely regarded as the "textbook example" of a cluster merger bow shock, the western shock front in the Bullet Cluster (1E0657-56) represents the ideal test case for such an SZ study. We aim to reconstruct a parametric model for the shock SZ signal by directly and jointly fitting deep, high-resolution interferometric data from the Atacama Large Millimeter/submillimeter Array (ALMA) and Atacama Compact Array (ACA) in Fourier space. The ALMA+ACA data are primarily sensitive to the electron pressure difference across the shock front. To estimate the shock Mach number $M$, this difference can be combined with the value for the upstream electron pressure derived from an independent Chandra X-ray analysis. In the case of instantaneous electron-ion temperature equilibration, we find $M=2.08^{+0.12}_{-0.12}$, in $\approx 2.4\sigma$ tension with the independent constraint from Chandra, $M_X=2.74\pm0.25$. The assumption of purely adiabatic electron temperature change across the shock leads to $M=2.53^{+0.33}_{-0.25}$, in better agreement with the X-ray estimate $M_X=2.57\pm0.23$ derived for the same heating scenario. We have demonstrated that interferometric observations of the SZ effect provide constraints on the properties of the shock in the Bullet Cluster that are highly complementary to X-ray observations. The combination of X-ray and SZ data yields a powerful probe of the shock properties, capable of measuring $M$ and addressing the question of electron-ion equilibration in cluster shocks. Our analysis is however limited by systematics related to the overall cluster geometry and the complexity of the post-shock gas distribution. To overcome these limitations, a joint analysis of SZ and X-ray data is needed.

We present cosmological parameter results from the final full-mission Planck measurements of the CMB anisotropies. We find good consistency with the standard spatially-flat 6-parameter $\Lambda$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $\Lambda$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $\Omega_c h^2 = 0.120\pm 0.001$, baryon density $\Omega_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $\tau = 0.054\pm 0.007$ (in this abstract we quote $68\,\%$ confidence regions on measured parameters and $95\,\%$ on upper limits). The angular acoustic scale is measured to $0.03\,\%$ precision, with $100\theta_*=1.0411\pm 0.0003$. These results are only weakly dependent on the cosmological model and remain stable, with somewhat increased errors, in many commonly considered extensions. Assuming the base-$\Lambda$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $\Omega_m = 0.315\pm 0.007$; and matter fluctuation amplitude $\sigma_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$\Lambda$CDM model. Combining with BAO we constrain the effective extra relativistic degrees of freedom to be $N_{\rm eff} = 2.99\pm 0.17$, and the neutrino mass is tightly constrained to $\sum m_\nu< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$\Lambda$CDM at over $2\,\sigma$, which pulls some parameters that affect the lensing amplitude away from the base-$\Lambda$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged)

This paper discusses the science case for a sensitive spectro-polarimetric survey of the microwave sky. Such a survey would provide a tomographic and dynamic census of the three-dimensional distribution of hot gas, velocity flows, early metals, dust, and mass distribution in the entire Hubble volume, exploit CMB temperature and polarisation anisotropies down to fundamental limits, and track energy injection and absorption into the radiation background across cosmic times by measuring spectral distortions of the CMB blackbody emission. In addition to its exceptional capability for cosmology and fundamental physics, such a survey would provide an unprecedented view of microwave emissions at sub-arcminute to few-arcminute angular resolution in hundreds of frequency channels, a data set that would be of immense legacy value for many branches of astrophysics. We propose that this survey be carried-out with a large space mission featuring a broad-band polarised imager and a moderate resolution spectro-imager at the focus of a 3.5m aperture telescope actively cooled to about 8K, complemented with absolutely-calibrated Fourier Transform Spectrometer modules observing at degree-scale angular resolution in the 10-2000 GHz frequency range. We propose two observing modes: a survey mode to map the entire sky as well as a few selected wide fields, and an observatory mode for deeper observations of regions of specific interest.

This Science White Paper, prepared in response to the ESA Voyage 2050 call for long-term mission planning, aims to describe the various science possibilities that can be realized with an L-class space observatory that is dedicated to the study of the interactions of cosmic microwave background (CMB) photons with the cosmic web. Our aim is specifically to use the CMB as a backlight -- and survey the gas, total mass, and stellar content of the entire observable Universe by means of analyzing the spatial and spectral distortions imprinted on it. These distortions result from two major processes that impact on CMB photons: scattering by electrons (Sunyaev-Zeldovich effect in diverse forms, Rayleigh scattering, resonant scattering) and deflection by gravitational potential (lensing effect). Even though the list of topics collected in this White Paper is not exhaustive, it helps to illustrate the exceptional diversity of major scientific questions that can be addressed by a space mission that will reach an angular resolution of 1.5 arcmin (goal 1 arcmin), have an average sensitivity better than 1 uK-arcmin, and span the microwave frequency range from roughly 50 GHz to 1 THz. The current paper also highlights the synergy of our BACKLIGHT mission concept with several upcoming and proposed ground-based CMB experiments.

Voyage 2050 White Paper highlighting the unique science opportunities using spectral distortions of the cosmic microwave background (CMB). CMB spectral distortions probe many processes throughout the history of the Universe. Precision spectroscopy, possible with existing technology, would provide key tests for processes expected within the cosmological standard model and open an enormous discovery space to new physics. This offers unique scientific opportunities for furthering our understanding of inflation, recombination, reionization and structure formation as well as dark matter and particle physics. A dedicated experimental approach could open this new window to the early Universe in the decades to come, allowing us to turn the long-standing upper distortion limits obtained with COBE/FIRAS some 25 years ago into clear detections of the expected standard distortion signals.

We report on the implications for cosmic inflation of the 2018 Release of the Planck CMB anisotropy measurements. The results are fully consistent with the two previous Planck cosmological releases, but have smaller uncertainties thanks to improvements in the characterization of polarization at low and high multipoles. Planck temperature, polarization, and lensing data determine the spectral index of scalar perturbations to be $n_\mathrm{s}=0.9649\pm 0.0042$ at 68% CL and show no evidence for a scale dependence of $n_\mathrm{s}.$ Spatial flatness is confirmed at a precision of 0.4% at 95% CL with the combination with BAO data. The Planck 95% CL upper limit on the tensor-to-scalar ratio, $r_{0.002}<0.10$, is further tightened by combining with the BICEP2/Keck Array BK15 data to obtain $r_{0.002}<0.056$. In the framework of single-field inflationary models with Einstein gravity, these results imply that: (a) slow-roll models with a concave potential, $V" (\phi) < 0,$ are increasingly favoured by the data; and (b) two different methods for reconstructing the inflaton potential find no evidence for dynamics beyond slow roll. Non-parametric reconstructions of the primordial power spectrum consistently confirm a pure power law. A complementary analysis also finds no evidence for theoretically motivated parameterized features in the Planck power spectrum, a result further strengthened for certain oscillatory models by a new combined analysis that includes Planck bispectrum data. The new Planck polarization data provide a stringent test of the adiabaticity of the initial conditions. The polarization data also provide improved constraints on inflationary models that predict a small statistically anisotropic quadrupolar modulation of the primordial fluctuations. However, the polarization data do not confirm physical models for a scale-dependent dipolar modulation.

This paper describes the 2018 Planck CMB likelihoods, following a hybrid approach similar to the 2015 one, with different approximations at low and high multipoles, and implementing several methodological and analysis refinements. With more realistic simulations, and better correction and modelling of systematics, we can now make full use of the High Frequency Instrument polarization data. The low-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$<800 and $\ell$>800 parameters or the preference for more smoothing of the $C_\ell$ peaks) are shown to be driven by the TT power spectrum and are not significantly modified by the inclusion of polarization. Overall, the legacy Planck CMB likelihoods provide a robust tool for constraining the cosmological model and represent a reference for future CMB observations. (Abridged)

We present measurements of the cosmic microwave background (CMB) lensing potential using the final $\textit{Planck}$ 2018 temperature and polarization data. We increase the significance of the detection of lensing in the polarization maps from $5\,\sigma$ to $9\,\sigma$. Combined with temperature, lensing is detected at $40\,\sigma$. We present an extensive set of tests of the robustness of the lensing-potential power spectrum, and construct a minimum-variance estimator likelihood over lensing multipoles $8 \le L \le 400$. We find good consistency between lensing constraints and the results from the $\textit{Planck}$ CMB power spectra within the $\rm{\Lambda CDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $\sigma_8 \Omega_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,\sigma$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $\sigma_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $\Omega_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $\sigma_8$ to better than $1\,\%$ precision, finding $\sigma_8=0.811\pm 0.006$. We find consistency with the lensing results from the Dark Energy Survey, and give combined lensing-only parameter constraints that are tighter than joint results using galaxy clustering. Using $\textit{Planck}$ cosmic infrared background (CIB) maps we make a combined estimate of the lensing potential over $60\,\%$ of the sky with considerably more small-scale signal. We demonstrate delensing of the $\textit{Planck}$ power spectra, detecting a maximum removal of $40\,\%$ of the lensing-induced power in all spectra. The improvement in the sharpening of the acoustic peaks by including both CIB and the quadratic lensing reconstruction is detected at high significance (abridged).

Analysis of the Planck 2018 data set indicates that the statistical properties of the cosmic microwave background (CMB) temperature anisotropies are in excellent agreement with previous studies using the 2013 and 2015 data releases. In particular, they are consistent with the Gaussian predictions of the $\Lambda$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on large angular scales. The novelty of the current study, however, lies in being a first attempt at a comprehensive analysis of the statistics of the polarization signal over all angular scales, using either maps of the Stokes parameters, $Q$ and $U$, or the $E$-mode signal derived from these using a new methodology (which we describe in an appendix). Although remarkable progress has been made in reducing the systematic effects that contaminated the 2015 polarization maps on large angular scales, it is still the case that residual systematics (and our ability to simulate them) can limit some tests of non-Gaussianity and isotropy. However, a detailed set of null tests applied to the maps indicates that these issues do not dominate the analysis on intermediate and large angular scales (i.e., $\ell \lesssim 400$). In this regime, no unambiguous detections of cosmological non-Gaussianity, or of anomalies corresponding to those seen in temperature, are claimed. Notably, the stacking of CMB polarization signals centred on the positions of temperature hot and cold spots exhibits excellent agreement with the $\Lambda$CDM cosmological model, and also gives a clear indication of how Planck provides state-of-the-art measurements of CMB temperature and polarization on degree scales.

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).

Following the pioneering observations with COBE in the early 1990s, studies of the cosmic microwave background (CMB) have focused on temperature and polarization anisotropies. CMB spectral distortions - tiny departures of the CMB energy spectrum from that of a perfect blackbody - provide a second, independent probe of fundamental physics, with a reach deep into the primordial Universe. The theoretical foundation of spectral distortions has seen major advances in recent years, which highlight the immense potential of this emerging field. Spectral distortions probe a fundamental property of the Universe - its thermal history - thereby providing additional insight into processes within the cosmological standard model (CSM) as well as new physics beyond. Spectral distortions are an important tool for understanding inflation and the nature of dark matter. They shed new light on the physics of recombination and reionization, both prominent stages in the evolution of our Universe, and furnish critical information on baryonic feedback processes, in addition to probing primordial correlation functions at scales inaccessible to other tracers. In principle the range of signals is vast: many orders of magnitude of discovery space could be explored by detailed observations of the CMB energy spectrum. Several CSM signals are predicted and provide clear experimental targets, some of which are already observable with present-day technology. Confirmation of these signals would extend the reach of the CSM by orders of magnitude in physical scale as the Universe evolves from the initial stages to its present form. The absence of these signals would pose a huge theoretical challenge, immediately pointing to new physics.

Observations of bright and variable "reflected" X-ray emission from molecular clouds located within inner hundred parsec of our Galaxy have demonstrated that the central supermassive black hole, Sgr A*, experienced short and powerful flares in the past few hundred years. These flares offer a truly unique opportunity to determine 3D location of the illuminated clouds (with ~10 pc accuracy) and to reveal their internal structure (down to 0.1 pc scales). Short duration of the flare(s), combined with X-rays high penetration power and insensitivity of the reflection signal to thermo- and chemo-dynamical state of the gas, ensures that the provided diagnostics of the density and velocity fields is unbiased and almost free of the projection and opacity effects. Sharp and sensitive snapshots of molecular gas accessible with aid of future X-ray observatories featuring large collecting area and high angular (arcsec-level) and spectral (eV-level) resolution cryogenic bolometers will present invaluable information on properties of the supersonic turbulence inside the illuminated clouds, map their shear velocity field and allow cross-matching between X-ray data and velocity-resolved emission of various molecular species provided by ALMA and other ground-based facilities. This will highlight large and small-scale dynamics of the dense gas and help uncovering specifics of the ISM lifecycle and high-mass star formation under very extreme conditions of galactic centers. While the former is of particular importance for the SMBH feeding and triggering AGN feedback, the latter might be an excellent test case for star formation taking place in high-redshift galaxies.

Sunyaev-Zeldovich (SZ) effects were first proposed in the 1970s as tools to identify the X-ray emitting hot gas inside massive clusters of galaxies and obtain their velocities relative to the cosmic microwave background (CMB). Yet it is only within the last decade that they have begun to significantly impact astronomical research. Thanks to the rapid developments in CMB instrumentation, measurement of the dominant thermal signature of the SZ effects has become a routine tool to find and characterize large samples of galaxy clusters and to seek deeper understanding of several important astrophysical processes via high-resolution imaging studies of many targets. With the notable exception of the Planck satellite and a few combinations of ground-based observatories, much of this "SZ revolution" has happened in the photometric mode, where observations are made at one or two frequencies in the millimeter regime to maximize the cluster detection significance and minimize the foregrounds. Still, there is much more to learn from detailed and systematic analyses of the SZ spectra across multiple wavelengths, specifically in the submillimeter (>300 GHz) domain. The goal of this Science White Paper is to highlight this particular aspect of SZ research, point out what new and potentially groundbreaking insights can be obtained from these studies, and emphasize why the coming decade can be a golden era for SZ spectral measurements.

Achieving a precise understanding of galaxy formation in a cosmological context is one of the great challenges in theoretical astrophysics, due to the vast range of spatial scales involved in the relevant physical processes. Observations in the millimeter bands, particularly those using the cosmic microwave background (CMB) radiation as a "backlight", provide a unique probe of the thermodynamics of these processes, with the capability to directly measure the density, pressure, and temperature of ionized gas. Moreover, these observations have uniquely high sensitivity into the outskirts of the halos of galaxies and clusters, including systems at high redshift. In the next decade, the combination of large spectroscopic and photometric optical galaxy surveys and wide-field, low-noise CMB surveys will transform our understanding of galaxy formation via these probes.

We present 353 GHz full-sky maps of the polarization fraction $p$, angle $\psi$, and dispersion of angles $S$ of Galactic dust thermal emission produced from the 2018 release of Planck data. We confirm that the mean and maximum of $p$ decrease with increasing $N_H$. The uncertainty on the maximum polarization fraction, $p_\mathrm{max}=22.0$% at 80 arcmin resolution, is dominated by the uncertainty on the zero level in total intensity. The observed inverse behaviour between $p$ and $S$ is interpreted with models of the polarized sky that include effects from only the topology of the turbulent Galactic magnetic field. Thus, the statistical properties of $p$, $\psi$, and $S$ mostly reflect the structure of the magnetic field. Nevertheless, we search for potential signatures of varying grain alignment and dust properties. First, we analyse the product map $S \times p$, looking for residual trends. While $p$ decreases by a factor of 3--4 between $N_H=10^{20}$ cm$^{-2}$ and $N_H=2\times 10^{22}$ cm$^{-2}$, $S \times p$ decreases by only about 25%, a systematic trend observed in both the diffuse ISM and molecular clouds. Second, we find no systematic trend of $S \times p$ with the dust temperature, even though in the diffuse ISM lines of sight with high $p$ and low $S$ tend to have colder dust. We also compare Planck data with starlight polarization in the visible at high latitudes. The agreement in polarization angles is remarkable. Two polarization emission-to-extinction ratios that characterize dust optical properties depend only weakly on $N_H$ and converge towards the values previously determined for translucent lines of sight. We determine an upper limit for the polarization fraction in extinction of 13%, compatible with the $p_\mathrm{max}$ observed in emission. These results provide strong constraints for models of Galactic dust in diffuse gas.

The Sunyaev-Zeldovich (SZ) effect was first predicted nearly five decades ago, but has only recently become a mature tool for performing high resolution studies of the warm and hot ionized gas in and between galaxies, groups, and clusters. Galaxy groups and clusters are powerful probes of cosmology, and they also serve as hosts for roughly half of the galaxies in the Universe. In this white paper, we outline the advances in our understanding of thermodynamic and kinematic properties of the warm-hot universe that can come in the next decade through spatially and spectrally resolved measurements of the SZ effects. Many of these advances will be enabled through new/upcoming millimeter/submillimeter (mm/submm) instrumentation on existing facilities, but truly transformative advances will require construction of new facilities with larger fields of view and broad spectral coverage of the mm/submm bands.

Using the ${\it Planck}$ full-mission data, we present a detection of the temperature (and therefore velocity) dispersion due to the kinetic Sunyaev-Zeldovich (kSZ) effect from clusters of galaxies. To suppress the primary CMB and instrumental noise we derive a matched filter and then convolve it with the ${\it Planck}$ foreground-cleaned `${\tt 2D-ILC\,}$' maps. By using the Meta Catalogue of X-ray detected Clusters of galaxies (MCXC), we determine the normalized ${\it rms}$ dispersion of the temperature fluctuations at the positions of clusters, finding that this shows excess variance compared with the noise expectation. We then build an unbiased statistical estimator of the signal, determining that the normalized mean temperature dispersion of $1526$ clusters is $\langle \left(\Delta T/T \right)^{2} \rangle = (1.64 \pm 0.48) \times 10^{-11}$. However, comparison with analytic calculations and simulations suggest that around $0.7\,\sigma$ of this result is due to cluster lensing rather than the kSZ effect. By correcting this, the temperature dispersion is measured to be $\langle \left(\Delta T/T \right)^{2} \rangle = (1.35 \pm 0.48) \times 10^{-11}$, which gives a detection at the $2.8\,\sigma$ level. We further convert uniform-weight temperature dispersion into a measurement of the line-of-sight velocity dispersion, by using estimates of the optical depth of each cluster (which introduces additional uncertainty into the estimate). We find that the velocity dispersion is $\langle v^{2} \rangle =(123\,000 \pm 71\,000)\,({\rm km}\,{\rm s}^{-1})^{2}$, which is consistent with findings from other large-scale structure studies, and provides direct evidence of statistical homogeneity on scales of $600\,h^{-1}{\rm Mpc}$. Our study shows the promise of using cross-correlations of the kSZ effect with large-scale structure in order to constrain the growth of structure.

This paper presents the High Frequency Instrument (HFI) data processing procedures for the Planck 2018 release. Major improvements in mapmaking have been achieved since the previous 2015 release. They enabled the first significant measurement of the reionization optical depth parameter using HFI data. This paper presents an extensive analysis of systematic effects, including the use of simulations to facilitate their removal and characterize the residuals. The polarized data, which presented a number of known problems in the 2015 Planck release, are very significantly improved. Calibration, based on the CMB dipole, is now extremely accurate and in the frequency range 100 to 353 GHz reduces intensity-to-polarization leakage caused by calibration mismatch. The Solar dipole direction has been determined in the three lowest HFI frequency channels to within one arc minute, and its amplitude has an absolute uncertainty smaller than $0.35\mu$K, an accuracy of order $10^{-4}$. This is a major legacy from the HFI for future CMB experiments. The removal of bandpass leakage has been improved by extracting the bandpass-mismatch coefficients for each detector as part of the mapmaking process; these values in turn improve the intensity maps. This is a major change in the philosophy of "frequency maps", which are now computed from single detector data, all adjusted to the same average bandpass response for the main foregrounds. Simulations reproduce very well the relative gain calibration of detectors, as well as drifts within a frequency induced by the residuals of the main systematic effect. Using these simulations, we measure and correct the small frequency calibration bias induced by this systematic effect at the $10^{-4}$ level. There is no detectable sign of a residual calibration bias between the first and second acoustic peaks in the CMB channels, at the $10^{-3}$ level.

We present new HST-COS G130M spectroscopy for a sightline toward a filament projected 1.9 kpc from the nucleus of M87. The combination of the sensitivity of COS and the proximity of M87 allows us to study the structure of this filament in unparalleled detail. We propose that the filament is composed of many cold clumps, each surrounded by an FUV-emitting boundary layer, with the filament having a radius $r_c \sim 10$ pc and the clumps filling the cylinder with a low volume filling factor. The observed velocity dispersion in emission lines from the filament results from the random motions of these clumps within the filament. We measure fluxes and kinematics for emission lines of Ly$\alpha$, C II $\lambda$1335, and N V $\lambda1238$. We associate these three lines, as well as archival measurements of H$\alpha$, C IV $\lambda$1549, and He II $\lambda$1640, with a multitemperature boundary layer around clumps which are moving with supersonic random motions in the filament. This boundary layer is a significant coolant of the hot gas. We show that the [C II] $\lambda$158$\mu$m flux observed by Herschel-PACS from this region implies the existence of a massive cold ($T \sim 10^3$ K) component in the filament which contains significantly more mass than the FUV-emitting boundary layer. We also detect [Fe XXI] $\lambda$1354 in emission at $4-5\sigma$. This line is emitted from 1 keV ($T \approx 10^7$ K) plasma, and we use it to measure the bulk radial velocity and velocity dispersion of the plasma at this temperature. In contrast to the intermediate-temperature FUV lines, [\ion{Fe}{xxi}] is blueshifted relative to M87 and matches the bulk velocity of a nearby filament to the south. We hypothesize that this line arises from the approaching face of the radio bubble expanding through this sightline, while the filament lies on the receding side of the bubble. (abstract is abridged)

We investigate the potential of large X-ray-selected AGN samples for detecting baryonic acoustic oscillations (BAO). Though AGN selection in X-ray band is very clean and efficient, it does not provide redshift information, and thus needs to be complemented with an optical follow-up. The main focus of this study is (i) to find the requirements needed for the quality of the optical follow-up and (ii) to formulate the optimal strategy of the X-ray survey, in order to detect the BAO. We demonstrate that redshift accuracy of $\sigma_0=10^{-2}$ at $z=1$ and the catastrophic failure rate of $f_{\rm fail}\lesssim 30\%$ are sufficient for a reliable detection of BAO in future X-ray surveys. Spectroscopic quality redshifts ($\sigma_0=10^{-3}$ and $f_{\rm fail}\sim 0$) will boost the confidence level of the BAO detection by a factor of $\sim 2$. For meaningful detection of BAO, X-ray surveys of moderate depth of $F_{\rm lim}\sim {\rm few~} 10^{-15}$ erg/s/cm$^2$ covering sky area from a few hundred to $\sim$ten thousand square degrees are required. The optimal strategy for the BAO detection does not necessarily require full sky coverage. For example, in a 1000-day long survey by an eROSITA type telescope, an optimal strategy would be to survey a sky area of $\sim 9000$ deg$^2$, yielding a $\sim 16\sigma$ BAO detection. A similar detection will be achieved by ATHENA+ or WFXT class telescopes in a survey with a duration of 100 days, covering a similar sky area. XMM-Newton can achieve a marginal BAO detection in a 100-day survey covering $\sim 400$ deg$^2$. These surveys would demand a moderate-to-high cost in terms the optical follow-ups, requiring determination of redshifts of $\sim 10^5$ (XMM-Newton) to $\sim 3\times 10^6$ objects (eROSITA, ATHENA+, and WFXT) in these sky areas.

We present cluster counts and corresponding cosmological constraints from the Planck full mission data set. Our catalogue consists of 439 clusters detected via their Sunyaev-Zeldovich (SZ) signal down to a signal-to-noise ratio of 6, and is more than a factor of 2 larger than the 2013 Planck cluster cosmology sample. The counts are consistent with those from 2013 and yield compatible constraints under the same modelling assumptions. Taking advantage of the larger catalogue, we extend our analysis to the two-dimensional distribution in redshift and signal-to-noise. We use mass estimates from two recent studies of gravitational lensing of background galaxies by Planck clusters to provide priors on the hydrostatic bias parameter, $(1-b)$. In addition, we use lensing of cosmic microwave background (CMB) temperature fluctuations by Planck clusters as an independent constraint on this parameter. These various calibrations imply constraints on the present-day amplitude of matter fluctuations in varying degrees of tension with those from the Planck analysis of primary fluctuations in the CMB; for the lowest estimated values of $(1-b)$ the tension is mild, only a little over one standard deviation, while it remains substantial ($3.7\,\sigma$) for the largest estimated value. We also examine constraints on extensions to the base flat $\Lambda$CDM model by combining the cluster and CMB constraints. The combination appears to favour non-minimal neutrino masses, but this possibility does little to relieve the overall tension because it simultaneously lowers the implied value of the Hubble parameter, thereby exacerbating the discrepancy with most current astrophysical estimates. Improving the precision of cluster mass calibrations from the current 10%-level to 1% would significantly strengthen these combined analyses and provide a stringent test of the base $\Lambda$CDM model.

Thanks to its high orbit and a set of complementary detectors providing continuous coverage of the whole sky, the INTEGRAL satellite has unique capabilities for the identification and study of the electromagnetic radiation associated to gravitational waves signals and, more generally, for multi-messenger astrophysics. Here we briefly review the results obtained during the first two observing runs of the advanced LIGO/Virgo interferometers.

We present the results of optical identifications and spectroscopic redshifts measurements for galaxy clusters from 2-nd Planck catalogue of Sunyaev-Zeldovich sources (PSZ2), located at high redshifts, $z\approx0.7-0.9$. We used the data of optical observations obtained with Russian-Turkish 1.5-m telescope (RTT150), Sayan observatory 1.6-m telescope, Calar Alto 3.5-m telescope and 6-m SAO RAS telescope (Bolshoi Teleskop Alt-azimutalnyi, BTA). Spectroscopic redshift measurements were obtained for seven galaxy clusters, including one cluster, PSZ2 G126.57+51.61, from the cosmological sample of PSZ2 catalogue. In central regions of two clusters, PSZ2 G069.39+68.05 and PSZ2 G087.39-34.58, the strong gravitationally lensed background galaxies are found, one of them at redshift $z=4.262$. The data presented below roughly double the number of known galaxy clusters in the second Planck catalogue of Sunyaev-Zeldovich sources at high redshifts, $z\approx0.8$.

We discuss the evolution of the ratio in number of recombinations due to 2s two photon escape and due to the escape of Lyman-$\alpha$ photons from the resonance during the epoch of cosmological recombination, within the width of the last scattering surface and near its boundaries. We discuss how this ratio evolves in time, and how it defines the profile of the Lyman-$\alpha$ line in the spectrum of CMB. One of the key reasons for explaining its time dependence is the strong overpopulation of the 2p level relative to the 2s level at redshifts $z \la 750$.

We investigate whether the intrinsic spin of neutron stars leaves an observable imprint on the spectral properties of X-ray binaries. To evaluate this we consider a sample of nine NSs for which the spins have been measured that are not accreting pulsars (for which the accretion geometry will be different). For each source, we perform spectroscopy on a majority of RXTE hard state observations. Our sample of sources and observations spans the range of the Eddington ratios Lx/Ledd ~0.005-0.100. We find a clear trend between key Comptonization properties and the NS spin for a given accretion rate. Specifically, at a given Lx/Ledd, for more rapidly rotating NSs we find lower seed photon temperatures and a general increase in Comptonization strength, as parametrised by the Comptonization y parameter and amplification factor A. This is in good agreement with the theoretical scenario whereby less energy is liberated in a boundary layer for more rapidly spinning NSs, resulting in a lower seed photon luminosity and, consequently, less Compton cooling in the corona. This effect in extremis results in the hard states of the most rapidly spinning sources encroaching upon the regime of Comptonization properties occupied by BHs.

We report the e INTernational Gamma-ray Astrophysics Laboratory (INTEGRAL) detection of the short gamma-ray burst GRB 170817A (discovered by Fermi-GBM) with a signal-to-noise ratio of 4.6, and, for the first time, its association with the gravitational waves (GWs) from binary neutron star (BNS) merging event GW170817 detected by the LIGO and Virgo observatories. The significance of association between the gamma-ray burst observed by INTEGRAL and GW170817 is 3.2 $\sigma$, while the association between the Fermi-GBM and INTEGRAL detections is 4.2 $\sigma$. GRB 170817A was detected by the SPI-ACS instrument about 2 s after the end of the gravitational wave event. We measure a fluence of $(1.4 \pm 0.4 \pm 0.6) \times$10$^{-7}$ erg cm$^{-2})$ (75--2000 keV), where, respectively, the statistical error is given at the 1 $\sigma$ confidence level, and the systematic error corresponds to the uncertainty in the spectral model and instrument response. We also report on the pointed follow-up observations carried out by INTEGRAL, starting 19.5 h after the event, and lasting for 5.4 days. We provide a stringent upper limit on any electromagnetic signal in a very broad energy range, from 3 keV to 8 MeV, constraining the soft gamma-ray afterglow flux to $<7.1\times$10$^{-11}$ erg cm$^{-2}$ s$^{-1}$ (80--300 keV). Exploiting the unique capabilities of INTEGRAL, we constrained the gamma-ray line emission from radioactive decays that are expected to be the principal source of the energy behind a kilonova event following a BNS coalescence. Finally, we put a stringent upper limit on any delayed bursting activity, for example from a newly formed magnetar.

We present the implications for cosmic inflation of the Planck measurements of the cosmic microwave background (CMB) anisotropies in both temperature and polarization based on the full Planck survey. The Planck full mission temperature data and a first release of polarization data on large angular scales measure the spectral index of curvature perturbations to be $n_\mathrm{s} = 0.968 \pm 0.006$ and tightly constrain its scale dependence to $d n_s/d \ln k =-0.003 \pm 0.007$ when combined with the Planck lensing likelihood. When the high-$\ell$ polarization data is included, the results are consistent and uncertainties are reduced. The upper bound on the tensor-to-scalar ratio is $r_{0.002} < 0.11$ (95% CL), consistent with the B-mode polarization constraint $r< 0.12$ (95% CL) obtained from a joint BICEP2/Keck Array and Planck analysis. These results imply that $V(\phi) \propto \phi^2$ and natural inflation are now disfavoured compared to models predicting a smaller tensor-to-scalar ratio, such as $R^2$ inflation. Three independent methods reconstructing the primordial power spectrum are investigated. The Planck data are consistent with adiabatic primordial perturbations. We investigate inflationary models producing an anisotropic modulation of the primordial curvature power spectrum as well as generalized models of inflation not governed by a scalar field with a canonical kinetic term. The 2015 results are consistent with the 2013 analysis based on the nominal mission data.

The discovery of a type I X-ray burst from the faint unidentified transient source IGR J17445-2747 in the Galactic bulge by the JEM-X telescope onboard the INTEGRAL observatory is reported. Type I bursts are believed to be associated with thermonuclear explosions of accreted matter on the surface of a neutron star with a weak magnetic field in a low-mass X-ray binary. Thus, this observation allows the nature of this source to be established.

We used data from the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) to set upper-limits on the ${\gamma}$-ray and hard X-ray prompt emission associated with the gravitational wave event GW170104, discovered by the LIGO/Virgo collaboration. The unique omni-directional viewing capability of the instruments on-board INTEGRAL allowed us to examine the full 90% confidence level localization region of the LIGO trigger. Depending on the particular spectral model assumed and the specific position within this region, the upper limits inferred from the INTEGRAL observations range from F${\gamma}$=1.9x10-7 erg cm-2 to F${\gamma}$=10-6 erg cm-2 (75 keV - 2 MeV energy range). This translates into a ratio between the prompt energy released in ${\gamma}$-rays along the direction to the observer and the gravitational wave energy of E${\gamma}$/EGW <2.6x10-5 . Using the INTEGRAL results, we can not confirm the ${\gamma}$-ray proposed counterpart to GW170104 by the AGILE team with the MCAL instrument. The reported flux of the AGILE/MCAL event, E2, is not compatible with the INTEGRAL upper limits within most of the 90% LIGO localization region. There is only a relatively limited portion of the sky where the sensitivity of the INTEGRAL instruments was not optimal and the lowest allowed fluence estimated for E2 would still be compatible with the INTEGRAL results. This region was also observed independently by Fermi/GBM and AstroSAT, from which, as far as we are aware, there are no reports of any significant detection of a prompt high-energy event.

Illumination of dense gas in the Central Molecular Zone (CMZ) by powerful X-ray flares from Sgr A* leads to prominent structures in the reflected emission that can be observed long after the end of the flare. By studying this emission we learn about past activity of the supermassive black hole in our Galactic Center and, at the same time, we obtain unique information on the structure of molecular clouds that is essentially impossible to get by other means. Here we discuss how X-ray data can improve our knowledge of both sides of the problem. Existing data already provide: i) an estimate of the flare age, ii) a model-independent lower limit on the luminosity of Sgr A* during the flare and iii) an estimate of the total emitted energy during Sgr A* flare. On the molecular clouds side, the data clearly show a voids-and-walls structure of the clouds and can provide an almost unbiased probe of the mass/density distribution of the molecular gas with the hydrogen column densities lower than few $10^{23}\;{\rm cm^{-2}}$. For instance, the probability distribution function of the gas density $PDF(\rho)$ can be measured this way. Future high energy resolution X-ray missions will provide the information on the gas velocities, allowing, for example a reconstruction of the velocity field structure functions and cross-matching the X-ray and molecular data based on positions and velocities.

We study the surface brightness fluctuations of the cosmic X-ray background (CXB) using Chandra data of XBOOTES. After masking out resolved sources we compute the power spectrum of fluctuations of the unresolved CXB for angular scales from ~2 arcsec to ~3 deg. The non-trivial large-scale structure (LSS) signal dominates over the shot-noise of unresolved point sources at all scales above ~1 arcmin and is produced mainly by the intracluster medium (ICM) of unresolved clusters and groups of galaxies, as shown in our previous publication. The shot-noise-subtracted power spectrum of CXB fluctuations has a power-law shape with the slope of $\Gamma = 0.96 \pm 0.06$. Its energy spectrum is well described by the redshifted emission spectrum of optically-thin plasma with the best-fit temperature of $T \approx 1.3$ keV and the best-fit redshift of $z \approx 0.40$. They are in good agreement with theoretical expectations based on the X-ray luminosity function and scaling relations of clusters. From these values we estimate the typical mass and luminosity of the objects responsible for CXB fluctuations, $M_{500} \sim 10^{13.6}\,{\rm M}_{\odot}/h$ and $L_{0.5-2.0\,{\rm keV}} \sim 10^{42.5}$ erg/s. On the other hand, the flux-weighted mean temperature and redshift of resolved clusters are $T \approx 2.4$ keV and $z \approx 0.23$, confirming that fluctuations of unresolved CXB are caused by cooler (i.e. less massive) and more distant clusters, as expected. We show that the power spectrum shape is sensitive to the ICM structure all the way to the outskirts, out to $\sim{\rm few}\times R_{500}$. We also look for possible contribution of the warm-hot intergalactic medium (WHIM) to the observed CXB fluctuations. Our results underline the significant diagnostics potential of the CXB fluctuation analysis in studying the ICM structure in clusters.

By using N-body hydrodynamical cosmological simulations in which the chemistry of major metals and molecules is consistently solved for, we study the interaction of metallic fine-structure lines with the CMB. Our analysis shows that the collisional induced emissions in the OI 145 $\mu$m and CII 158 $\mu$m lines during reionization introduce a distortion of the CMB spectrum at low frequencies ($\nu < 300$ GHz) with amplitudes up to $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-8}$-$10^{-7}$, i.e., at the $\sim 0.1$ percent level of FIRAS upper limits. Shorter wavelength fine-structure transitions (OI 63 $\mu$m, FeII 26 $\mu$m, and SiII 35 $\mu$m) typically sample the reionization epoch at higher observing frequencies ($\nu > 400$ GHz). This corresponds to the Wien tail of the CMB spectrum and the distortion level induced by those lines may be as high as $\Delta I_{\nu}/B_{\nu}(T_{\rm CMB})\sim 10^{-4}$. The angular anisotropy produced by these lines should be more relevant at higher frequencies: while practically negligible at $\nu=145 $GHz, signatures from CII 158 $\mu$m and OI 145 $\mu$m should amount to 1%-5% of the anisotropy power measured at $l \sim 5000$ and $\nu=220 $GHz by the ACT and SPT collaborations (after assuming $\Delta \nu_{\rm obs}/\nu_{\rm obs}\simeq 0.005$ for the line observations). Our simulations show that anisotropy maps from different lines (e.g., OI 145 $\mu$m and CII 158 $\mu$m) at the same redshift show a very high degree ($>0.8$) of spatial correlation, allowing for the use of observations at different frequencies to unveil the same snapshot of the reionization epoch. Finally, our simulations demonstrate that line-emission anisotropies extracted in narrow frequency/redshift shells are practically uncorrelated in frequency space, thus enabling standard methods for removal of foregrounds that vary smoothly in frequency, just as in HI 21 cm studies.

X-ray spectra from cores of galaxy clusters can be strongly distorted by resonant scattering of line photons, affecting metal abundance and gas velocity measurements. We introduce simulated spectral models that take into account the resonant scattering effect, radial variations of thermodynamic properties of the hot gas, projection effects and small-scale isotropic gas motions. The key feature of the models is that all these effects are treated self-consistently for the whole spectrum, rather than for individual lines. The model spectra are publicly available and can be used for direct comparison with observed projected spectra. Comparison with the existing XMM-Newton and Chandra data of the Perseus Cluster shows that even though there is no strong evidence for the resonant scattering in Perseus, the low energy resolution of the X-ray CCDs is not sufficient to robustly distinguish spectral distortions due to the resonant scattering, different metal abundance profiles and different levels of gas turbulence. Future Astro-H data will resolve most of the problems we are facing with CCDs. With the help of our models, the resonant scattering analysis can be done self-consistently using the whole spectral information, constraining the level of gas turbulence already with a 100 ks observation with Astro-H.

Future observations of cosmic microwave background (CMB) polarisation have the potential to answer some of the most fundamental questions of modern physics and cosmology. In this paper, we list the requirements for a future CMB polarisation survey addressing these scientific objectives, and discuss the design drivers of the CORE space mission proposed to ESA in answer to the "M5" call for a medium-sized mission. The rationale and options, and the methodologies used to assess the mission's performance, are of interest to other future CMB mission design studies. CORE is designed as a near-ultimate CMB polarisation mission which, for optimal complementarity with ground-based observations, will perform the observations that are known to be essential to CMB polarisation scienceand cannot be obtained by any other means than a dedicated space mission.

A type Ia supernova is thought to be a thermonuclear explosion of either a single carbon-oxygen white dwarf or of a pair of merging white dwarfs. The explosion fuses a large amount of radioactive 56Ni. After the explosion, the decay chain from 56Ni to 56Co to 56Fe generates gamma-ray photons, which are reprocessed in the expanding ejecta and give rise to powerful optical emission. Here we report the detection of 56Co lines at energies of 847 and 1238 keV and a gamma-ray continuum in the 200-400 keV band from the type Ia supernova 2014J in the nearby galaxy M82. The line fluxes suggest that about 0.6 +/- 0.1 solar masses of radioactive 56Ni were synthesized during the explosion. The line broadening gives a characteristic mass-weighted ejecta expansion velocity of 10000 +/- 3000 km/s. The observed gamma-ray properties are in broad agreement with the canonical model of an explosion of a white dwarf just massive enough to be unstable to gravitational collapse, but do not immediately exclude more complicated merger scenarios, which fuse comparable amount of 56Ni.

The International Gamma-Ray Astrophysics Laboratory (INTEGRAL) continues to successfully work in orbit after its launch in 2002. The mission provides the deepest ever survey of hard X-ray sources throughout the Galaxy at energies above 20 keV. We report on a catalogue of new hard X-ray source candidates based on the latest sky maps comprising 14 years of data acquired with the IBIS telescope onboard INTEGRAL in the Galactic Plane (|b|<17.5 deg). The current catalogue includes in total 72 hard X-ray sources detected at S/N>4.7 sigma and not known to previous INTEGRAL surveys. Among them, 31 objects have also been detected in the on-going all-sky survey by the BAT telescope of the Swift observatory. For 26 sources on the list, we suggest possible identifications: 21 active galactic nuclei, two cataclysmic variables, two isolated pulsars or pulsar wind nebulae, and one supernova remnant; 46 sources from the catalogue remain unclassified.

During the first observing run of LIGO, two gravitational wave events and one lower-significance trigger (LVT151012) were reported by the LIGO/Virgo collaboration. At the time of LVT151012, the INTErnational Gamma-Ray Astrophysics Laboratory (INTEGRAL) was pointing at a region of the sky coincident with the high localization probability area of the event and thus permitted us to search for its electromagnetic counterpart (both prompt and afterglow emission). The imaging instruments on-board INTEGRAL (IBIS/ISGRI, IBIS/PICsIT, SPI, and the two JEM-X modules) have been exploited to attempt the detection of any electromagnetic emission associated with LVT151012 over 3 decades in energy (from 3 keV to 8 MeV). The omni-directional instruments on-board the satellite, i.e. the SPI-ACS and IBIS monitored the entire LVT151012 localization region at energies above 75 keV. We did not find any significant transient source that was spatially and/or temporally coincident with LVT151012, obtaining tight upper limits on the associated hard X-ray and $\gamma$-ray radiation. For typical spectral models, the upper limits on the fluence of the emission from any 1 s long-lasting counterpart of LVT151012 ranges from $F_{\gamma}=$3.5$\times$10$^{-8}$ erg cm$^{-2}$ (20 - 200 keV) to $F_{\gamma}$=7.1$\times$10$^{-7}$ erg cm$^{-2}$ (75 - 2000 keV), constraining the ratio of the isotropic equivalent energy released in the electromagnetic emission to the total energy of the gravitational waves: $E_{75-2000~keV}/E_{GW}<$4.4$\times$10$^{-5}$. Finally, we provide an exhaustive summary of the capabilities of all instruments on-board INTEGRAL to hunt for $\gamma$-ray counterparts of gravitational wave events, exploiting both serendipitous and pointed follow-up observations. This will serve as a reference for all future searches.

We use a model of the molecular gas distribution within ~100 pc from the center of the Milky Way (Kruijssen, Dale & Longmore) to simulate time evolution and polarization properties of the reflected X-ray emission, associated with the past outbursts from Sgr A*. While this model is too simple to describe the complexity of the true gas distribution, it illustrates the importance and power of long-term observations of the reflected emission. We show that the variable part of X-ray emission observed by Chandra and XMM from prominent molecular clouds is well described by a pure reflection model, providing strong support of the reflection scenario. While the identification of Sgr A* as a primary source for this reflected emission is already a very appealing hypothesis, a decisive test of this model can be provided by future X-ray polarimetric observations, that will allow placing constraints on the location of the primary source. In addition, X-ray polarimeters (like, e.g., XIPE) have sufficient sensitivity to constrain the line-of-sight positions of molecular complexes, removing major uncertainty in the model.

We present the results of spectroscopic redshift measurements for the galaxy clusters from the first all-sky Planck catalogue of the Sunyaev-Zeldovich sources, that have been mostly identified by means of the optical observations performed previously by our team (Planck Collaboration, 2015a). The data on 13 galaxy clusters at redshifts from z=~0.2 to z=~0.8, including the improved identification and redshift measurement for the cluster PSZ1 G141.73+14.22 at z=0.828, are provided. The measurements were done using the data from Russian-Turkish 1.5-m telescope (RTT-150), 2.2-m Calar Alto Observatory telescope, and 6-m SAO RAS telescope (Bolshoy Teleskop Azimutalnyi, BTA).

This paper presents a study of the ISW effect from the Planck 2015 temperature and polarization data release. The CMB is cross-correlated with different LSS tracers: the NVSS, SDSS and WISE catalogues, and the Planck 2015 lensing map. This cross-correlation yields a detection at $4\,\sigma$, where most of the signal-to-noise is due to the Planck lensing and NVSS. In fact, the ISW effect is detected only from the Planck data (through the ISW-lensing bispectrum) at $\approx 3\,\sigma$, which is similar to the detection level achieved by combining the cross-correlation signal coming from all the catalogues. The ISW signal allow us to detect $\Omega_\Lambda$ at more than $3\,\sigma$. This cross-correlation analysis is performed only with the Planck temperature data, since the polarization scales available in the 2015 release do not permit significant improvement of the CMB-LSS cross-correlation detectability. Nevertheless, polarization data is used to study the anomalously large ISW signal previously reported through the stacking of CMB features at the locations of known superstructures. We find that the current Planck polarization data do not exclude that this signal could be caused by the ISW effect. In addition, the stacking of the Planck lensing map on the locations of superstructures exhibits a positive cross-correlation with these large-scale structures. Finally, we have improved our previous reconstruction of the ISW temperature fluctuations by combining the information encoded in all the previously mentioned LSS tracers. In particular, we construct a map of the ISW secondary anisotropies and the corresponding uncertainties map, obtained from simulations. We also explore the reconstruction of the ISW anisotropies caused by the LSS traced by the 2MPZ survey by directly inverting the density field into the gravitational potential field.

During the scanning observations of the Galactic Center region in August - September 2016 we detected the new outburst of the historical X-ray nova GRS 1739-278, the black hole candidate LMXB system. In this letter we present results of INTEGRAL and Swift-XRT observations taken during the outburst. In hard X-ray band (20-60 keV) the flux from the source raised from $\sim$11 to $\sim$30 mCrab between 3 and 14 of September. For nearly 8 days the source has been observed at this flux level and then faded to $\sim$15 mCrab. The broadband quasi-simultaneous spectrum obtained during the outburst is well described by the absorbed powerlaw with the photon index $\Gamma=1.86\pm0.07$ in broad energy range 0.5-150 keV, with absorption corresponding to ${N_{H}}=2.3\times10^{22}$ cm$^{-2}$ assuming solar abundance. Based on this we can conclude that the source was in the low/hard state. From the lightcurve and spectra we propose that this outburst was `failed', i.e. amount of accreted matter was not sufficient to achieve the high/soft spectral state with dominant soft blackbody component as seen in normal outbursts of black hole candidates.

A powerful outburst of X-ray radiation from the supermassive black hole Sgr A* at the center of the Milky Way is believed to be responsible for the illumination of molecular clouds in the central ~100 pc of the Galaxy (Sunyaev et al., 1993, Koyama et al., 1996). The reflected/reprocessed radiation comes to us with a delay corresponding to the light propagation time that depends on the 3D position of molecular clouds with respect to Sgr A*. We suggest a novel way of determining the age of the outburst and positions of the clouds by studying characteristic imprints left by the outburst in the spatial and time variations of the reflected emission. We estimated the age of the outburst that illuminates the Sgr A molecular complex to be ~110 yr. This estimate implies that we see the gas located ~10 pc further away from us than Sgr A*. If the Sgr B2 complex is also illuminated by the same outburst, then it is located ~130 pc closer than our Galactic Center. The outburst was short (less than a few years) and the total amount of emitted energy in X-rays is $\displaystyle \sim 10^{48}\rho_3^{-1}$ erg, where $\rho_3$ is the mean hydrogen density of the cloud complex in units of $10^3 {\rm cm^{-3}}$. Energetically, such fluence can be provided by a partial tidal disruption event or even by a capture of a planet. Further progress in more accurate positioning and timing of the outburst should be possible with future X-ray polarimetric observations and long-term systematic observations with Chandra and XMM-Newton. A few hundred-years long X-ray observations would provide a detailed 3D map of the gas density distribution in the central $\sim 100$ pc region.

We analyse the spectra of black hole (BH) and neutron star (NS) X-ray binaries (XBs) in the hard state using archival RXTE observations. We find that there is a clear dichotomy in the strength of Comptonisation between NS and BH sources, as measured by both the Compton y-parameter and amplification factor A, with distinct groups of BH and NS XBs separated at y~0.9 and A~3. The electron temperature kTe can occupy a broad range in BH systems, from kTe~30-200 keV, whereas for NSs kTe is peaked at ~15-25 keV, but can extend to higher values. The difference between BHs and NSs in y implies that kTe is higher at a given optical depth for BH XBs. Our results also imply that for NS systems the accreting material loses ~1/2-2/3 of its energy through Comptonisation in the corona. The remaining energy is released on the surface of the neutron star, making it a powerful source of soft radiation, which alters the properties of the Comptonising corona. Finally, we find evidence at the 2.4 sigma confidence level that Comptonisation parameters may be correlated with the neutron star spin, whereas no correlation with the BH spin is found. Our results highlight a further observational distinction between BH and NS XBs that is a consequence of NSs possessing a physical surface.

The secondary cosmic microwave background (CMB) $B$-modes stem from the post-decoupling distortion of the polarization $E$-modes due to the gravitational lensing effect of large-scale structures. These lensing-induced $B$-modes constitute both a valuable probe of the dark matter distribution and an important contaminant for the extraction of the primary CMB $B$-modes from inflation. Planck provides accurate nearly all-sky measurements of both the polarization $E$-modes and the integrated mass distribution via the reconstruction of the CMB lensing potential. By combining these two data products, we have produced an all-sky template map of the lensing-induced $B$-modes using a real-space algorithm that minimizes the impact of sky masks. The cross-correlation of this template with an observed (primordial and secondary) $B$-mode map can be used to measure the lensing $B$-mode power spectrum at multipoles up to $2000$. In particular, when cross-correlating with the $B$-mode contribution directly derived from the Planck polarization maps, we obtain lensing-induced $B$-mode power spectrum measurement at a significance level of $12\,\sigma$, which agrees with the theoretical expectation derived from the Planck best-fit $\Lambda$CDM model. This unique nearly all-sky secondary $B$-mode template, which includes the lensing-induced information from intermediate to small ($10\lesssim \ell\lesssim 1000$) angular scales, is delivered as part of the Planck 2015 public data release. It will be particularly useful for experiments searching for primordial $B$-modes, such as BICEP2/Keck Array or LiteBIRD, since it will enable an estimate to be made of the lensing-induced contribution to the measured total CMB $B$-modes.

We present the most significant measurement of the cosmic microwave background (CMB) lensing potential to date (at a level of 40 sigma), using temperature and polarization data from the Planck 2015 full-mission release. Using a polarization-only estimator we detect lensing at a significance of 5 sigma. We cross-check the accuracy of our measurement using the wide frequency coverage and complementarity of the temperature and polarization measurements. Public products based on this measurement include an estimate of the lensing potential over approximately 70% of the sky, an estimate of the lensing potential power spectrum in bandpowers for the multipole range 40<L<400 and an associated likelihood for cosmological parameter constraints. We find good agreement between our measurement of the lensing potential power spectrum and that found in the best-fitting LCDM model based on the Planck temperature and polarization power spectra. Using the lensing likelihood alone we obtain a percent-level measurement of the parameter combination $\sigma_8 \Omega_m^{0.25} = 0.591\pm 0.021$. We combine our determination of the lensing potential with the E-mode polarization also measured by Planck to generate an estimate of the lensing B-mode. We show that this lensing B-mode estimate is correlated with the B-modes observed directly by Planck at the expected level and with a statistical significance of 10 sigma, confirming Planck's sensitivity to this known sky signal. We also correlate our lensing potential estimate with the large-scale temperature anisotropies, detecting a cross-correlation at the 3 sigma level, as expected due to dark energy in the concordance LCDM model.

Fluctuations of the surface brightness of cosmic X-ray background (CXB) carry unique information about faint and low luminosity source populations, which is inaccessible for conventional large-scale structure (LSS) studies based on resolved sources. We used Chandra data of the XBOOTES field ($\sim9\,\mathrm{deg^2}$) to conduct the most accurate measurement to date of the power spectrum of fluctuations of the unresolved CXB on the angular scales of $\sim3\,$arcsec $-$ $\sim17\,$arcmin. We find that at sub-arcmin angular scales, the power spectrum is consistent with the AGN shot noise, without much need for any significant contribution from their one-halo term. This is consistent with the theoretical expectation that low-luminosity AGN reside alone in their dark matter halos. However, at larger angular scales we detect a significant LSS signal above the AGN shot noise. Its power spectrum, obtained after subtracting the AGN shot noise, follows a power law with the slope of $-0.8\pm0.1$ and its amplitude is much larger than what can be plausibly explained by the two-halo term of AGN. We demonstrate that the detected LSS signal is produced by unresolved clusters and groups of galaxies. For the flux limit of the XBOOTES survey, their flux-weighted mean redshift equals $\left<z\right>\sim0.3$, and the mean temperature of their intracluster medium (ICM), $\left<T\right>\approx 1.4$ keV, corresponds to the mass of $M_{500} \sim 10^{13.5}\,\mathrm{M}_\odot$. The power spectrum of CXB fluctuations carries information about the redshift distribution of these objects and the spatial structure of their ICM on the linear scales of up to $\sim$Mpc, i.e. of the order of the virial radius.

We investigate constraints on cosmic reionization extracted from the Planck cosmic microwave background (CMB) data. We combine the Planck CMB anisotropy data in temperature with the low-multipole polarization data to fit LCDM models with various parameterizations of the reionization history. We obtain a Thomson optical depth tau=0.058 +/- 0.012 for the commonly adopted instantaneous reionization model. This confirms, with only data from CMB anisotropies, the low value suggested by combining Planck 2015 results with other data sets and also reduces the uncertainties. We reconstruct the history of the ionization fraction using either a symmetric or an asymmetric model for the transition between the neutral and ionized phases. To determine better constraints on the duration of the reionization process, we also make use of measurements of the amplitude of the kinetic Sunyaev-Zeldovich (kSZ) effect using additional information from the high resolution Atacama Cosmology Telescope and South Pole Telescope experiments. The average redshift at which reionization occurs is found to lie between z=7.8 and 8.8, depending on the model of reionization adopted. Using kSZ constraints and a redshift-symmetric reionization model, we find an upper limit to the width of the reionization period of Dz < 2.8. In all cases, we find that the Universe is ionized at less than the 10% level at redshifts above z~10. This suggests that an early onset of reionization is strongly disfavoured by the Planck data. We show that this result also reduces the tension between CMB-based analyses and constraints from other astrophysical sources.

We calculate the one-point probability density distribution functions (PDF) and the power spectra of the thermal and kinetic Sunyaev-Zeldovich (tSZ and kSZ) effects and the mean Compton Y parameter using the Magneticum Pathfinder simulations, state-of-the-art cosmological hydrodynamical simulations of a large cosmological volume of (896 Mpc/h)^3. These simulations follow in detail the thermal and chemical evolution of the intracluster medium as well as the evolution of super-massive black holes and their associated feedback processes. We construct full-sky maps of tSZ and kSZ from the light-cones out to z=0.17, and one realization of 8.8x8.8 degree wide, deep light-cone out to z=5.2. The local universe at z<0.027 is simulated by a constrained realisation. The tail of the one-point PDF of tSZ from the deep light-cone follows a power-law shape with an index of -3.2. Once convolved with the effective beam of Planck, it agrees with the PDF measured by Planck. The predicted tSZ power spectrum agrees with that of the Planck data at all multipoles up to l~1000, once the calculations are scaled to the Planck 2015 cosmological parameters with \Omega_m=0.308 and \sigma_8=0.8149. Consistent with the results in the literature, however, we continue to find the tSZ power spectrum at l=3000 that is significantly larger than that estimated from the high-resolution ground-based data. The simulation predicts the mean fluctuating Compton Y value of <Y>=1.18x10^{-6} for \Omega_m=0.272 and \sigma_8=0.809. Nearly half (~ 5x10^{-7}) of the signal comes from halos below a virial mass of 10^{13}M_\odot/h. Scaling this to the Planck 2015 parameters, we find <Y>=1.57x10^{-6}. The PDF and the power spectrum of kSZ from our simulation agree broadly with the previous work.

The Planck full mission cosmic microwave background(CMB) temperature and E-mode polarization maps are analysed to obtain constraints on primordial non-Gaussianity(NG). Using three classes of optimal bispectrum estimators - separable template-fitting (KSW), binned, and modal - we obtain consistent values for the local, equilateral, and orthogonal bispectrum amplitudes, quoting as our final result from temperature alone fNL^local=2.5+\-5.7, fNL^equil=-16+\-70 and fNL^ortho=-34+\-33(68%CL). Combining temperature and polarization data we obtain fNL^local=0.8+\-5.0, fNL^equil=-4+\-43 and fNL^ortho=-26+\-21 (68%CL). The results are based on cross-validation of these estimators on simulations, are stable across component separation techniques, pass an extensive suite of tests, and are consistent with Minkowski functionals based measurements. The effect of time-domain de-glitching systematics on the bispectrum is negligible. In spite of these test outcomes we conservatively label the results including polarization data as preliminary, owing to a known mismatch of the noise model in simulations and the data. Beyond fNL estimates, we present model-independent reconstructions of the CMB bispectrum and derive constraints on early universe scenarios that generate NG, including general single-field and axion inflation, initial state modifications, parity-violating tensor bispectra, and directionally dependent vector models. We also present a wide survey of scale-dependent oscillatory bispectra, and we look for isocurvature NG. Our constraint on the local primordial trispectrum amplitude is gNL^local=(-9.0+\-7.7)x10^4 (68%CL), and we perform an analysis of additional trispectrum shapes. The global picture is one of consistency with the premises of the LambdaCDM cosmology, namely that the structure we observe today was sourced by adiabatic, passive, Gaussian, and primordial seed perturbations.[abridged]

This paper presents the Planck 2015 likelihoods, statistical descriptions of the 2-point correlations of CMB data, using the hybrid approach employed previously: pixel-based at $\ell<30$ and a Gaussian approximation to the distribution of spectra at higher $\ell$. The main improvements are the use of more and better processed data and of Planck polarization data, and more detailed foreground and instrumental models, allowing further checks and enhanced immunity to systematics. Progress in foreground modelling enables a larger sky fraction. Improvements in processing and instrumental models further reduce uncertainties. For temperature, we perform an analysis of end-to-end instrumental simulations fed into the data processing pipeline; this does not reveal biases from residual instrumental systematics. The $\Lambda$CDM cosmological model continues to offer a very good fit to Planck data. The slope of primordial scalar fluctuations, $n_s$, is confirmed smaller than unity at more than 5{\sigma} from Planck alone. We further validate robustness against specific extensions to the baseline cosmology. E.g., the effective number of neutrino species remains compatible with the canonical value of 3.046. This first detailed analysis of Planck polarization concentrates on E modes. At low $\ell$ we use temperature at all frequencies and a subset of polarization. The frequency range improves CMB-foreground separation. Within the baseline model this requires a reionization optical depth $\tau=0.078\pm0.019$, significantly lower than without high-frequency data for explicit dust monitoring. At high $\ell$ we detect residual errors in E, typically O($\mu$K$^2$); we recommend temperature alone as the high-$\ell$ baseline. Nevertheless, Planck high-$\ell$ polarization allows a separate determination of $\Lambda$CDM parameters consistent with those from temperature alone.

We present results based on full-mission Planck observations of temperature and polarization anisotropies of the CMB. These data are consistent with the six-parameter inflationary LCDM cosmology. From the Planck temperature and lensing data, for this cosmology we find a Hubble constant, H0= (67.8 +/- 0.9) km/s/Mpc, a matter density parameter Omega_m = 0.308 +/- 0.012 and a scalar spectral index with n_s = 0.968 +/- 0.006. (We quote 68% errors on measured parameters and 95% limits on other parameters.) Combined with Planck temperature and lensing data, Planck LFI polarization measurements lead to a reionization optical depth of tau = 0.066 +/- 0.016. Combining Planck with other astrophysical data we find N_ eff = 3.15 +/- 0.23 for the effective number of relativistic degrees of freedom and the sum of neutrino masses is constrained to < 0.23 eV. Spatial curvature is found to be |Omega_K| < 0.005. For LCDM we find a limit on the tensor-to-scalar ratio of r <0.11 consistent with the B-mode constraints from an analysis of BICEP2, Keck Array, and Planck (BKP) data. Adding the BKP data leads to a tighter constraint of r < 0.09. We find no evidence for isocurvature perturbations or cosmic defects. The equation of state of dark energy is constrained to w = -1.006 +/- 0.045. Standard big bang nucleosynthesis predictions for the Planck LCDM cosmology are in excellent agreement with observations. We investigate annihilating dark matter and deviations from standard recombination, finding no evidence for new physics. The Planck results for base LCDM are in agreement with BAO data and with the JLA SNe sample. However the amplitude of the fluctuations is found to be higher than inferred from rich cluster counts and weak gravitational lensing. Apart from these tensions, the base LCDM cosmology provides an excellent description of the Planck CMB observations and many other astrophysical data sets.

(abridged) We discuss the Galactic foreground emission between 20 and 100GHz based on observations by Planck/WMAP. The Commander component-separation tool has been used to separate the various astrophysical processes in total intensity. Comparison with RRL templates verifies the recovery of the free-free emission along the Galactic plane. Comparison of the high-latitude Halpha emission with our free-free map shows residuals that correlate with dust optical depth, consistent with a fraction (~30%) of Halpha having been scattered by high-latitude dust. We highlight a number of diffuse spinning dust morphological features at high latitude. There is substantial spatial variation in the spinning dust spectrum, with the emission peak ranging from below 20GHz to more than 50GHz. There is a strong tendency for the spinning dust component near many prominent HII regions to have a higher peak frequency, suggesting that this increase in peak frequency is associated with dust in the photodissociation regions around the nebulae. The emissivity of spinning dust in these diffuse regions is of the same order as previous detections in the literature. Over the entire sky, the commander solution finds more anomalous microwave emission than the WMAP component maps, at the expense of synchrotron and free-free emission. This can be explained by the difficulty in separating multiple broadband components with a limited number of frequency maps. Future surveys (5-20GHz), will greatly improve the separation by constraining the synchrotron spectrum. We combine Planck/WMAP data to make the highest S/N ratio maps yet of the intensity of the all-sky polarized synchrotron emission at frequencies above a few GHz. Most of the high-latitude polarized emission is associated with distinct large-scale loops and spurs, and we re-discuss their structure...

This paper describes the identification, modelling, and removal of previously unexplained systematic effects in the polarization data of the Planck High Frequency Instrument (HFI) on large angular scales, including new mapmaking and calibration procedures, new and more complete end-to-end simulations, and a set of robust internal consistency checks on the resulting maps. These maps, at 100, 143, 217, and 353 GHz, are early versions of those that will be released in final form later in 2016. The improvements allow us to determine the cosmic reionization optical depth $\tau$ using, for the first time, the low-multipole $EE$ data from HFI, reducing significantly the central value and uncertainty, and hence the upper limit. Two different likelihood procedures are used to constrain $\tau$ from two estimators of the CMB $E$- and $B$-mode angular power spectra at 100 and 143 GHz, after debiasing the spectra from a small remaining systematic contamination. These all give fully consistent results. A further consistency test is performed using cross-correlations derived from the Low Frequency Instrument maps of the Planck 2015 data release and the new HFI data. For this purpose, end-to-end analyses of systematic effects from the two instruments are used to demonstrate the near independence of their dominant systematic error residuals. The tightest result comes from the HFI-based $\tau$ posterior distribution using the maximum likelihood power spectrum estimator from $EE$ data only, giving a value $0.055\pm 0.009$. In a companion paper these results are discussed in the context of the best-fit Planck $\Lambda$CDM cosmological model and recent models of reionization.

Recent models for the large-scale Galactic magnetic fields in the literature have been largely constrained by synchrotron emission and Faraday rotation measures. We use three different but representative models to compare their predicted polarized synchrotron and dust emission with that measured by the Planck satellite. We first update these models to match the Planck synchrotron products using a common model for the cosmic-ray leptons. We discuss the impact on this analysis of the ongoing problems of component separation in the Planck microwave bands and of the uncertain cosmic-ray spectrum. In particular, the inferred degree of ordering in the magnetic fields is sensitive to these systematic uncertainties, and we further show the importance of considering the expected variations in the observables in addition to their mean morphology. We then compare the resulting simulated emission to the observed dust polarization and find that the dust predictions do not match the morphology in the Planck data but underpredict the dust polarization away from the plane. We modify one of the models to roughly match both observables at high latitudes by increasing the field ordering in the thin disc near the observer. Though this specific analysis is dependent on the component separation issues, we present the improved model as a proof of concept for how these studies can be advanced in future using complementary information from ongoing and planned observational projects.

We study the implications of Planck data for models of dark energy (DE) and modified gravity (MG), beyond the cosmological constant scenario. We start with cases where the DE only directly affects the background evolution, considering Taylor expansions of the equation of state, principal component analysis and parameterizations related to the potential of a minimally coupled DE scalar field. When estimating the density of DE at early times, we significantly improve present constraints. We then move to general parameterizations of the DE or MG perturbations that encompass both effective field theories and the phenomenology of gravitational potentials in MG models. Lastly, we test a range of specific models, such as k-essence, f(R) theories and coupled DE. In addition to the latest Planck data, for our main analyses we use baryonic acoustic oscillations, type-Ia supernovae and local measurements of the Hubble constant. We further show the impact of measurements of the cosmological perturbations, such as redshift-space distortions and weak gravitational lensing. These additional probes are important tools for testing MG models and for breaking degeneracies that are still present in the combination of Planck and background data sets. All results that include only background parameterizations are in agreement with LCDM. When testing models that also change perturbations (even when the background is fixed to LCDM), some tensions appear in a few scenarios: the maximum one found is \sim 2 sigma for Planck TT+lowP when parameterizing observables related to the gravitational potentials with a chosen time dependence; the tension increases to at most 3 sigma when external data sets are included. It however disappears when including CMB lensing.

Non-thermal pressure from turbulence and bulk flows is a fundamental ingredient in hot gaseous halos, and in the intracluster medium it will be measured through emission line kinematics with calorimeters on future X-ray spacecraft. In this paper we present a complementary method for measuring these effects, using forbidden FUV emission lines of highly ionized Iron which trace $10^7$ K gas. The brightest of these is [Fe XXI] $\lambda$1354.1. We search for these lines in archival HST-COS spectra from the well-known elliptical galaxies M87 and NGC4696, which harbor large reservoirs of $10^7$ K gas. We report a 2.2$\sigma$ feature which we attribute to [Fe XXI] from a filament in M87, and positive residuals in the nuclei of M87 and NGC4696, for which the 90\% upper limits on the line flux are close to the predicted fluxes based on X-ray observations. In a newer reduction of the data from the Hubble Spectroscopic Legacy Archive, these limits become tighter and the [Fe XXI] feature reaches a formal significance of 5.3$\sigma$, neglecting uncertainty in fitting the continuum. Using our constraints, we perform emission measure analysis, constraining the characteristic path length and column density of the $\sim10^7$ K gas. We also examine several sightlines towards filaments or cooling flows in other galaxy clusters, for which the fraction of gas at $10^7$ K is unknown, and place upper limits on its emission measure in each case. A medium-resolution HST-COS observation of the M87 filament for $\sim$10 orbits would confirm our detection of [Fe XXI] and measure its width.

Although infrared (IR) overall dust emission from clusters of galaxies has been statistically detected using data from the Infrared Astronomical Satellite (IRAS), it has not been possible to sample the spectral energy distribution (SED) of this emission over its peak, and thus to break the degeneracy between dust temperature and mass. By complementing the IRAS spectral coverage with Planck satellite data from 100 to 857 GHz, we provide new constraints on the IR spectrum of thermal dust emission in clusters of galaxies. We achieve this by using a stacking approach for a sample of several hundred objects from the Planck cluster sample; this procedure averages out fluctuations from the IR sky, allowing us to reach a significant detection of the faint cluster contribution. We also use the large frequency range probed by Planck, together with component-separation techniques, to remove the contamination from both cosmic microwave background anisotropies and the thermal Sunyaev-Zeldovich effect (tSZ) signal, which dominate below 353 GHz. By excluding dominant spurious signals or systematic effects, averaged detections are reported at frequencies between 353 and 5000 GHz. We confirm the presence of dust in clusters of galaxies at low and intermediate redshifts, yielding an SED with a shape similar to that of the Milky Way. Planck's beam does not allow us to investigate the detailed spatial distribution of this emission (e.g., whether it comes from intergalactic dust or simply the dust content of the cluster galaxies), but the radial distribution of the emission appears to follow that of the stacked SZ signal, and thus the extent of the clusters. The recovered SED allows us to constrain the dust mass responsible for the signal, as well as its temperature. We additionally explore the evolution of the IR emission as a function of cluster mass and redshift.

We report the results of the deepest Galactic Plane (|b| < 17.5 deg) survey in the 67.9 and 78.4 keV nuclear de-excitation lines of titanium-44 (44Ti) performed using the data acquired with the IBIS/ISGRI instrument onboard the INTEGRAL satellite during 12 years of operation. The peak sensitivity of our survey reached an unprecedented level of $4.8\times10^{-6}$ erg s$^{-1}$ cm$^{-2}$ (3 sigma) that improves the sensitivity of the survey done by CGRO/COMPTEL by a factor of ~5. As a result, constraining upper limits for all sources from the catalog of Galactic supernova remnants (SNRs; Green 2014) are derived. These upper limits can be used to estimate the exposure needed to detect 44Ti emission from any known SNR using existing and prospective X- and gamma-ray telescopes. Among the youngest Galactic SNRs, only Cas A shows significant 44Ti emission flux in good agreement with the NuSTAR measurements. We did not detect any other sources of titanium emission in the Galactic Plane at significance level higher than 5 sigma confirming previous claims of the rarity of such 44Ti-producing SNRs.

We compute and investigate four types of imprint of a stochastic background of primordial magnetic fields (PMFs) on the cosmic microwave background (CMB) anisotropies: the impact of PMFs on the CMB spectra; the effect on CMB polarization induced by Faraday rotation; the impact of PMFs on the ionization history; magnetically-induced non-Gaussianities; and the magnetically-induced breaking of statistical isotropy. Overall, Planck data constrain the amplitude of PMFs to less than a few nanogauss. In particular, individual limits coming from the analysis of the CMB angular power spectra, using the Planck likelihood, are $B_{1\,\mathrm{Mpc}}< 4.4$ nG (where $B_{1\,\mathrm{Mpc}}$ is the comoving field amplitude at a scale of 1 Mpc) at 95% confidence level, assuming zero helicity, and $B_{1\,\mathrm{Mpc}}< 5.6$ nG for a maximally helical field.For nearly scale-invariant PMFs we obtain $B_{1\,\mathrm{Mpc}}<2.0$ nG and $B_{1\,\mathrm{Mpc}}<0.9$ nG if the impact of PMFs on the ionization history of the Universe is included. From the analysis of magnetically-induced non-Gaussianity we obtain three different values, corresponding to three applied methods, all below 5 nG. The constraint from the magnetically-induced passive-tensor bispectrum is $B_{1\,\mathrm{Mpc}}< 2.8$ nG. A search for preferred directions in the magnetically-induced passive bispectrum yields $B_{1\,\mathrm{Mpc}}< 4.5$ nG, whereas the the compensated-scalar bispectrum gives $B_{1\,\mathrm{Mpc}}< 3$ nG. The analysis of the Faraday rotation of CMB polarization by PMFs uses the Planck power spectra in $EE$ and $BB$ at 70 GHz and gives $B_{1\,\mathrm{Mpc}}< 1380$ nG. In our final analysis, we consider the harmonic-space correlations produced by Alfv\'en waves, finding no significant evidence for the presence of these waves. Together, these results comprise a comprehensive set of constraints on possible PMFs with Planck data.

We present all-sky modelling of the high resolution Planck, IRAS, and WISE infrared (IR) observations using the physical dust model presented by Draine and Li in 2007 (DL). We study the performance and results of this model, and discuss implications for future dust modelling. The present work extends the DL dust modelling carried out on nearby galaxies using Herschel and Spitzer data to Galactic dust emission. We employ the DL dust model to generate maps of the dust mass surface density, the optical extinction Av, and the starlight intensity parametrized by Umin. The DL model reproduces the observed spectral energy distribution (SED) satisfactorily over most of the sky, with small deviations in the inner Galactic disk and in low ecliptic latitude areas. We compare the DL optical extinction Av for the diffuse interstellar medium with optical estimates for 2 10^5 quasi-stellar objects (QSOs) observed in the Sloan digital sky survey. The DL Av estimates are larger than those determined towards QSOs by a factor of about 2, which depends on Umin. The DL fitting parameter Umin, effectively determined by the wavelength where the SED peaks, appears to trace variations in the far-IR opacity of the dust grains per unit Av, and not only in the starlight intensity. To circumvent the model deficiency, we propose an empirical renormalization of the DL Av estimate, dependent of Umin, which compensates for the systematic differences found with QSO observations. This renormalization also brings into agreement the DL Av estimates with those derived for molecular clouds from the near-IR colours of stars in the 2 micron all sky survey. The DL model and the QSOs data are used to compress the spectral information in the Planck and IRAS observations for the diffuse ISM to a family of 20 SEDs normalized per Av, parameterized by Umin, which may be used to test and empirically calibrate dust models.

The Second Planck Catalogue of Compact Sources is a catalogue of sources detected in single-frequency maps from the full duration of the Planck mission and supersedes previous versions of the Planck compact source catalogues. It consists of compact sources, both Galactic and extragalactic, detected over the entire sky. Compact sources detected in the lower frequency channels are assigned to the PCCS2, while at higher frequencies they are assigned to one of two sub-catalogues, the PCCS2 or PCCS2E, depending on their location on the sky. The first of these catalogues covers most of the sky and allows the user to produce subsamples at higher reliabilities than the target 80% integral reliability of the catalogue. The PCCS2E contains sources detected in sky regions where the diffuse emission makes it difficult to quantify the reliability of the detections. Both the PCCS2 and PCCS2E include polarization measurements, in the form of polarized flux densities, or upper limits, and orientation angles for all seven polarization-sensitive Planck channels. The improved data-processing of the full-mission maps and their reduced noise levels allow us to increase the number of objects in the catalogue, improving its completeness for the target 80 % reliability as compared with the previous versions, the PCCS and ERCSC catalogues.

We present the results of approximately three years of observations of Planck Sunyaev-Zeldovich (SZ) sources with telescopes at the Canary Islands observatories as part of the general optical follow-up programme undertaken by the Planck collaboration. In total, 78 SZ sources are discussed. Deep-imaging observations were obtained for most of these sources; spectroscopic observations in either in long-slit or multi-object modes were obtained for many. We effectively used 37.5 clear nights. We found optical counterparts for 73 of the 78 candidates. This sample includes 53 spectroscopic redshift determinations, 20 of them obtained with a multi-object spectroscopic mode. The sample contains new redshifts for 27 Planck clusters that were not included in the first Planck SZ source catalogue (PSZ1).

The optical light curve of Type Ia supernovae (SNIa) is powered by thermalized gamma-rays produced by the decay of 56Ni and 56Co, the main radioactive isotopes synthesized by the thermonuclear explosion of a C/O white dwarf. Gamma-rays escaping the ejecta can be used as a diagnostic tool for studying the characteristics of the explosion. In particular, it is expected that the analysis of the early gamma emission, near the maximum of the optical light curve, could provide information about the distribution of the radioactive elements in the debris. In this paper, the gamma data obtained from SN2014J in M82 by the instruments on board of INTEGRAL are analyzed taking special care of the impact that the detailed spectral response has on the measurements of the intensity of the lines. The 158 keV emission of 56Ni has been detected in SN2014J at ~5 sigma at low energy with both ISGRI and SPI around the maximum of the optical light curve. After correcting the spectral response of the detector, the fluxes in the lines suggest that, in addition to the bulk of radioactive elements buried in the central layers of the debris, there is a plume of 56Ni, with a significance of ~3 sigma, moving at high velocity and receding from the observer. The mass of the plume is in the range of ~ 0.03-0.08 solar masses. No SNIa explosion model had predicted the mass and geometrical distribution of 56Ni suggested here. According to its optical properties, SN2014J looks as a normal SNIa. So it is extremely important to discern if it is also representative in the gamma-ray band.

We present results of deep surveys of three extragalactic fields, M81 (exposure of 9.7 Ms), LMC (6.8 Ms) and 3C 273/Coma (9.3 Ms), in the hard X-ray (17-60 keV) energy band with the IBIS telescope onboard the INTEGRAL observatory, based on 12 years of observations (2003-2015). The combined survey reaches a $4\sigma$ peak sensitivity of 0.18 mCrab (2.6$\times$10$^{-12}$ erg s$^{-1}$ cm$^{-2}$) and sensitivity better than 0.25 and 0.87 mCrab over 10% and 90% of its full area of 4900 deg$^{2}$, respectively. We have detected in total 147 sources at $S/N>4\sigma$, including 37 sources observed in hard X-rays for the first time. The survey is dominated by extragalactic sources, mostly by active galactic nuclei (AGN). The sample of identified sources contains 98 AGN (including 64 Seyfert galaxies, 7 LINERs, 3 XBONGs, 16 blazars and 8 AGN of unclear optical class), two galaxy clusters (Coma and Abell 3266), 17 objects located in the Large and Small Magellanic Clouds (13 high- and 2 low-mass X-ray binaries and 2 X-ray pulsars), three Galactic cataclysmic variables, one ultraluminous X-ray source (ULX, M82 X-1) and one blended source (SWIFT J1105.7+5854). The nature of 25 sources remains unknown, so that the surveys identification is currently complete at 83%. We have constructed AGN number-flux relations ($\log{N}$-$\log{S}$) and calculated AGN number densities in the local Universe for the entire survey and for each of the three extragalactic fields.

We test the statistical isotropy and Gaussianity of the cosmic microwave background (CMB) anisotropies using observations made by the Planck satellite. Our results are based mainly on the full Planck mission for temperature, but also include some polarization measurements. In particular, we consider the CMB anisotropy maps derived from the multi-frequency Planck data by several component-separation methods. For the temperature anisotropies, we find excellent agreement between results based on these sky maps over both a very large fraction of the sky and a broad range of angular scales, establishing that potential foreground residuals do not affect our studies. Tests of skewness, kurtosis, multi-normality, N-point functions, and Minkowski functionals indicate consistency with Gaussianity, while a power deficit at large angular scales is manifested in several ways, for example low map variance. The results of a peak statistics analysis are consistent with the expectations of a Gaussian random field. The "Cold Spot" is detected with several methods, including map kurtosis, peak statistics, and mean temperature profile. We thoroughly probe the large-scale dipolar power asymmetry, detecting it with several independent tests, and address the subject of a posteriori correction. Tests of directionality suggest the presence of angular clustering from large to small scales, but at a significance that is dependent on the details of the approach. We perform the first examination of polarization data, finding the morphology of stacked peaks to be consistent with the expectations of statistically isotropic simulations. Where they overlap, these results are consistent with the Planck 2013 analysis based on the nominal mission data and provide our most thorough view of the statistics of the CMB fluctuations to date.

Cores of relaxed galaxy clusters are often disturbed by AGN. Their Chandra observations revealed a wealth of structures induced by shocks, subsonic gas motions, bubbles of relativistic plasma, etc. In this paper, we determine the nature and energy content of gas fluctuations in the Perseus core by probing statistical properties of emissivity fluctuations imprinted in the soft- and hard-band X-ray images. About 80 per cent of the total variance of perturbations on ~ 8-70 kpc scales in the inner region have an isobaric nature, i.e., are consistent with slow displacements of the gas in pressure equilibrium with ambient medium. Observed variance translates to the ratio of non-thermal to thermal energy of ~13 per cent. In the region dominated by weak "ripples", about half of the total variance is also associated with isobaric perturbations on scales ~ a few tens of kpc. If these isobaric perturbations are induced by buoyantly rising bubbles, then these results suggest that most of the AGN-injected energy should first go into bubbles rather than into shocks. Using simulations of a shock propagating through the Perseus atmosphere, we found that models reproducing the observed features of a central shock have more than 50 per cent of the AGN-injected energy associated with the bubble enthalpy and only about 20 per cent is carried away with the shock. Such energy partition is consistent with the AGN-feedback model, mediated by bubbles of relativistic plasma, and supports the importance of turbulence in the balance between gas heating and radiative cooling.

The Virgo cluster is the largest Sunyaev-Zeldovich (SZ) source in the sky, both in terms of angular size and total integrated flux. Planck's wide angular scale and frequency coverage, together with its high sensitivity, allow a detailed study of this large object through the SZ effect. Virgo is well resolved by Planck, showing an elongated structure, which correlates well with the morphology observed from X-rays, but extends beyond the observed X-ray signal. We find a good agreement between the SZ signal (or Compton paranmeter, y_c) observed by Planck and the expected signal inferred from X-ray observations and simple analytical models. Due to its proximity to us, the gas beyond the virial radius can be studied with unprecedented sensitivity by integrating the SZ signal over tens of square degrees. We study the signal in the outskirts of Virgo and compare it with analytical models and a constrained simulation of the environment of Virgo. Planck data suggest that significant amounts of low-density plasma surround Virgo out to twice the virial radius. We find the SZ signal in the outskirts of Virgo to be consistent with a simple model that extrapolates the inferred pressure at lower radii while assuming that the temperature stays in the keV range beyond the virial radius. The observed signal is also consistent with simulations and points to a shallow pressure profile in the outskirts of the cluster. This reservoir of gas at large radii can be linked with the hottest phase of the elusive warm/hot intergalactic medium. Taking the lack of symmetry of Virgo into account, we find that a prolate model is favoured by the combination of SZ and X-ray data, in agreement with predictions.