The turbulent amplification of cosmic magnetic fields depends upon the material properties of the host plasma. In many hot, dilute astrophysical systems, such as the intracluster medium (ICM) of galaxy clusters, the rarity of particle--particle collisions allows departures from local thermodynamic equilibrium. These departures exert anisotropic viscous stresses on the plasma motions that inhibit their ability to stretch magnetic-field lines. We present a numerical study of the fluctuation dynamo in a weakly collisional plasma using magnetohydrodynamic (MHD) equations endowed with a field-parallel viscous (Braginskii) stress. When the stress is limited to values consistent with a pressure anisotropy regulated by firehose and mirror instabilities, the Braginskii-MHD dynamo largely resembles its MHD counterpart. If instead the parallel viscous stress is left unabated -- a situation relevant to recent kinetic simulations of the fluctuation dynamo and to the early stages of the dynamo in a magnetized ICM -- the dynamo changes its character, amplifying the magnetic field while exhibiting many characteristics of the saturated state of the large-Prandtl-number (${\rm Pm}\gtrsim{1}$) MHD dynamo. We construct an analytic model for the Braginskii-MHD dynamo in this regime, which successfully matches magnetic-energy spectra. A prediction of this model, confirmed by our simulations, is that a Braginskii-MHD plasma without pressure-anisotropy limiters will not support a dynamo if the ratio of perpendicular and parallel viscosities is too small. This ratio reflects the relative allowed rates of field-line stretching and mixing, the latter of which promotes resistive dissipation of the magnetic field. In all cases that do exhibit a dynamo, the generated magnetic field is organized into folds that persist into the saturated state and bias the chaotic flow to acquire a scale-dependent spectral anisotropy.

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. The net effect of the improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is ($3366.6 \pm 2.7$)$\mu$K, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of $\tau = 0.051 \pm 0.006$, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations.

We study cosmological observables on the past light cone of a fixed observer in the context of clustering dark energy. We focus on observables that probe the gravitational field directly, namely the integrated Sachs-Wolfe and non-linear Rees-Sciama effect (ISW-RS), weak gravitational lensing, gravitational redshift and Shapiro time delay. With our purpose-built $N$-body code "$k$-evolution" that tracks the coupled evolution of dark matter particles and the dark energy field, we are able to study the regime of low speed of sound $c_s$ where dark energy perturbations can become quite large. Using ray tracing we produce two-dimensional sky maps for each effect and we compute their angular power spectra. It turns out that the ISW-RS signal is the most promising probe to constrain clustering dark energy properties coded in $w-c_s^2$, as the $\textit{linear}$ clustering of dark energy would change the angular power spectrum by $\sim 30\%$ at low $\ell$ when comparing two different speeds of sound for dark energy. Weak gravitational lensing, Shapiro time-delay and gravitational redshift are less sensitive probes of clustering dark energy, showing variations of a few percent only. The effect of dark energy $\textit{non-linearities}$ in all the power spectra is negligible at low $\ell$, but reaches about $2\%$ and $3\%$, respectively, in the convergence and ISW-RS angular power spectra at multipoles of a few hundred when observed at redshift $\sim 0.85$. Future cosmological surveys achieving percent precision measurements will allow to probe the clustering of dark energy to a high degree of confidence. Clear evidence for clustering dark energy at any scale would rule out a cosmological constant as the leading contender for dark energy, and may help to distinguish between different dark energy models.

We perform null tests of the concordance model, using $H(z)$ measurements that mimic next-generation surveys such as Euclid and the SKA. To this end, we deploy a non-parametric method, so that we make minimal assumptions about the fiducial cosmology as well as the statistical analysis. We produce simulations assuming different cosmological models in order to verify how well we can distinguish between their signatures. We find that SKA- and Euclid-like surveys should be able to discriminate sharply between the concordance and alternative dark energy models that are compatible with the Planck CMB data. We conclude that SKA and Euclid will be able to falsify the concordance model in a statistically significant way, if one of the benchmarks models represents the true Universe, without making assumptions about the underlying cosmology.

The Complete Calibration of the Colour-Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed to empirically calibrate the galaxy colour-redshift relation - P(z|C) to the Euclid depth (i_AB=24.5) and is intimately linked to upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@ VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This paper focuses on high-quality spectroscopic redshifts of high-z galaxies observed with the KMOS spectrograph in the H- and K-bands. A total of 424 highly-reliable z are measured in the 1.3<=z<=2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined z fill 55% of high and 35% of lower priority empty SOM grid cells. We measured Halpha fluxes in a 1."2 radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B-V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z<1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z>2 galaxies.

Hybrid-Vlasov-Maxwell simulations of magnetized plasma turbulence including non-linear electron-inertia effects in a generalized Ohm's law are presented. When fluctuation energy is injected on scales sufficiently close to ion-kinetic scales, the ions efficiently become de-magnetized and electron-scale current sheets largely dominate the distribution of the emerging current structures, in contrast to the usual picture, where a full hierarchy of structure sizes is generally observed. These current sheets are shown to be the sites of electron-only reconnection (e-rec), in which the usual electron exhausts are unaccompanied by ion outflows and which are in qualitative agreement with those recently observed by MMS in the Earth's turbulent magnetosheath, downstream of the bow shock. Some features of the e-rec phenomenology are shown to be consistent with an electron magnetohydrodynamic description. Simulations suggest that this regime of collisionless reconnection may be found in turbulent systems where plasma processes, such as micro-instabilities and/or shocks, overpower the more customary turbulent cascade by directly injecting energy close to the ion-kinetic scales.

Weak lensing, which is the deflection of light by matter along the line of sight, has proven to be an efficient method for constraining models of structure formation and reveal the nature of dark energy. So far, most weak-lensing studies have focused on the shear field that can be measured directly from the ellipticity of background galaxies. However, within the context of forthcoming full-sky weak-lensing surveys such as Euclid, convergence maps (mass maps) offer an important advantage over shear fields in terms of cosmological exploitation. While it carry the same information, the lensing signal is more compressed in the convergence maps than in the shear field. This simplifies otherwise computationally expensive analyses, for instance, non-Gaussianity studies. However, the inversion of the non-local shear field requires accurate control of systematic effects caused by holes in the data field, field borders, shape noise, and the fact that the shear is not a direct observable (reduced shear). We present the two mass-inversion methods that are included in the official Euclid data-processing pipeline: the standard Kaiser & Squires method (KS), and a new mass-inversion method (KS+) that aims to reduce the information loss during the mass inversion. This new method is based on the KS method and includes corrections for mass-mapping systematic effects. The results of the KS+ method are compared to the original implementation of the KS method in its simplest form, using the Euclid Flagship mock galaxy catalogue. In particular, we estimate the quality of the reconstruction by comparing the two-point correlation functions and third- and fourth- order moments obtained from shear and convergence maps, and we analyse each systematic effect independently and simultaneously...

Using analytical theory and hybrid-kinetic numerical simulations, we demonstrate that, in a collisionless plasma, long-wavelength ion-acoustic waves (IAWs) with amplitudes $\delta n/n_0 \gtrsim 2/\beta$ (where $\beta\gg{1}$ is the ratio of thermal to magnetic pressure) generate sufficient pressure anisotropy to destabilize the plasma to firehose and mirror instabilities. These kinetic instabilities grow rapidly to reduce the pressure anisotropy by pitch-angle scattering and trapping particles, respectively, thereby impeding the maintenance of Landau resonances that enable such waves' otherwise potent collisionless damping. The result is wave dynamics that evince a weakly collisional plasma: the ion distribution function is near-Maxwellian, the field-parallel flow of heat resembles its Braginskii form (except in regions where large-amplitude magnetic mirrors strongly suppress particle transport), and the relations between various thermodynamic quantities are more `fluid-like' than kinetic. A nonlinear fluctuation-dissipation relation for self-sustaining IAWs is obtained by solving a plasma-kinetic Langevin problem, which demonstrates suppressed damping, enhanced fluctuation levels, and weakly collisional thermodynamics when IAWs with $\delta n/n_0 \gtrsim 2/\beta$ are stochastically driven. We investigate how our results depend upon the scale separation between the wavelength of the IAW and the Larmor radius of the ions, and discuss briefly their implications for our understanding of turbulence and transport in the solar wind and the intracluster medium of galaxy clusters.

Outside galaxy clusters the competition between the inwards gravitational attraction and the outwards expansion of the Universe leads to a special radius of velocity cancellation, which is called the turn-around radius. Measurements of the turn-around radius hold promises of constraining cosmological parameters, and possibly even properties of gravity. Such a measurement is, however, complicated by the fact that the surroundings of galaxy clusters are not spherical, but instead are a complicated collection of filaments, sheets and voids. In this paper we use the results of numerically simulated universes to quantify realistic error-bars of the measurement of the turn-around radius. We find that for a $\Lambda$CDM cosmology these error-bars are typically of the order of $20\%$. We numerically simulate three different implementations of dark energy models and of a scalar dark sector interaction to address whether the turn-around radius can be used to constrain non-trivial cosmologies, and we find that only rather extreme models can be distinguished from a $\Lambda$CDM universe due to the large error-bars arising from the non-trivial cluster environments.

We introduce $k$-evolution, a relativistic $N$-body code based on $\textit{gevolution}$, which includes clustering dark energy among its cosmological components. To describe dark energy, we use the effective field theory approach. In particular, we focus on $k$-essence with a speed of sound much smaller than unity but we lay down the basis to extend the code to other dark energy and modified gravity models. We develop the formalism including dark energy non-linearities but, as a first step, we implement the equations in the code after dropping non-linear self-coupling in the $k$-essence field. In this simplified setup, we compare $k$-evolution simulations with those of $\texttt{CLASS}$ and $\textit{gevolution}$ 1.2, showing the effect of dark matter and gravitational non-linearities on the power spectrum of dark matter, of dark energy and of the gravitational potential. Moreover, we compare $k$-evolution to Newtonian $N$-body simulations with back-scaled initial conditions and study how dark energy clustering affects massive halos.

In this paper, we quantify the non-linear effects from $k$-essence dark energy through an effective parameter $\mu$ that encodes the additional contribution of a dark energy fluid or a modification of gravity to the Poisson equation. This is a first step toward quantifying non-linear effects of dark energy/modified gravity models in a more general approach. We compare our $N$-body simulation results from $k$-evolution with predictions from the linear Boltzmann code $\texttt{CLASS}$, and we show that for the $k$-essence model one can safely neglect the difference between the two potentials, $ \Phi -\Psi$, and short wave corrections appearing as higher order terms in the Poisson equation, which allows us to use single parameter $\mu$ for characterizing this model. We also show that for a large $k$-essence speed of sound the $\texttt{CLASS}$ results are sufficiently accurate, while for a low speed of sound non-linearities in matter and in the $k$-essence field are non-negligible. We propose a $\tanh$-based parameterisation for $\mu$, motivated by the results for two cases with low ($c_s^2=10^{-7}$) and high ($c_s^2=10^{-4}$) speed of sound, to include the non-linear effects based on the simulation results. This parametric form of $\mu$ can be used to improve Fisher forecasts or Newtonian $N$-body simulations for $k$-essence models.

The data from the Euclid mission will enable the measurement of the photometric redshifts, angular positions, and weak lensing shapes for over a billion galaxies. This large dataset will allow for cosmological analyses using the angular clustering of galaxies and cosmic shear. The cross-correlation (XC) between these probes can tighten constraints and it is therefore important to quantify their impact for Euclid. In this study we carefully quantify the impact of XC not only on the final parameter constraints for different cosmological models, but also on the nuisance parameters. In particular, we aim at understanding the amount of additional information that XC can provide for parameters encoding systematic effects, such as galaxy bias or intrinsic alignments (IA). We follow the formalism presented in Euclid Collaboration: Blanchard et al. (2019) and make use of the codes validated therein. We show that XC improves the dark energy Figure of Merit (FoM) by a factor $\sim 5$, whilst it also reduces the uncertainties on galaxy bias by $\sim 17\%$ and the uncertainties on IA by a factor $\sim 4$. We observe that the role of XC on the final parameter constraints is qualitatively the same irrespective of the galaxy bias model used. We also show that XC can help in distinguishing between different IA models, and that if IA terms are neglected then this can lead to significant biases on the cosmological parameters. We find that the XC terms are necessary to extract the full information content from the data in future analyses. They help in better constraining the cosmological model, and lead to a better understanding of the systematic effects that contaminate these probes. Furthermore, we find that XC helps in constraining the mean of the photometric-redshift distributions, but it requires a more precise knowledge of this mean in order not to degrade the final FoM. [Abridged]

Context. Future weak lensing surveys, such as the Euclid mission, will attempt to measure the shapes of billions of galaxies in order to derive cosmological information. These surveys will attain very low levels of statistical error, and systematic errors must be extremely well controlled. In particular, the point spread function (PSF) must be estimated using stars in the field, and recovered with high accuracy. Aims. The aims of this paper are twofold. Firstly, we took steps toward a nonparametric method to address the issue of recovering the PSF field, namely that of finding the correct PSF at the position of any galaxy in the field, applicable to Euclid. Our approach relies solely on the data, as opposed to parametric methods that make use of our knowledge of the instrument. Secondly, we studied the impact of imperfect PSF models on the shape measurement of galaxies themselves, and whether common assumptions about this impact hold true in an Euclid scenario. Methods. We extended the recently proposed resolved components analysis approach, which performs super-resolution on a field of under-sampled observations of a spatially varying, image-valued function. We added a spatial interpolation component to the method, making it a true 2-dimensional PSF model. We compared our approach to PSFEx, then quantified the impact of PSF recovery errors on galaxy shape measurements through image simulations. Results. Our approach yields an improvement over PSFEx in terms of the PSF model and on observed galaxy shape errors, though it is at present far from reaching the required Euclid accuracy. We also find that the usual formalism used for the propagation of PSF model errors to weak lensing quantities no longer holds in the case of an Euclid-like PSF. In particular, different shape measurement approaches can react differently to the same PSF modeling errors.

Stage IV weak lensing experiments will offer more than an order of magnitude leap in precision. We must therefore ensure that our analyses remain accurate in this new era. Accordingly, previously ignored systematic effects must be addressed. In this work, we evaluate the impact of the reduced shear approximation and magnification bias, on the information obtained from the angular power spectrum. To first-order, the statistics of reduced shear, a combination of shear and convergence, are taken to be equal to those of shear. However, this approximation can induce a bias in the cosmological parameters that can no longer be neglected. A separate bias arises from the statistics of shear being altered by the preferential selection of galaxies and the dilution of their surface densities, in high-magnification regions. The corrections for these systematic effects take similar forms, allowing them to be treated together. We calculated the impact of neglecting these effects on the cosmological parameters that would be determined from Euclid, using cosmic shear tomography. To do so, we employed the Fisher matrix formalism, and included the impact of the super-sample covariance. We also demonstrate how the reduced shear correction can be calculated using a lognormal field forward modelling approach. These effects cause significant biases in Omega_m, sigma_8, n_s, Omega_DE, w_0, and w_a of -0.53 sigma, 0.43 sigma, -0.34 sigma, 1.36 sigma, -0.68 sigma, and 1.21 sigma, respectively. We then show that these lensing biases interact with another systematic: the intrinsic alignment of galaxies. Accordingly, we develop the formalism for an intrinsic alignment-enhanced lensing bias correction. Applying this to Euclid, we find that the additional terms introduced by this correction are sub-dominant.

This white paper summarizes major scientific challenges and opportunities in understanding magnetic reconnection and related explosive phenomena as a fundamental plasma process.

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 Euclid mission will observe well over a billion galaxies out to $z\sim6$ and beyond. This will offer an unrivalled opportunity to investigate several key questions for understanding galaxy formation and evolution. The first step for many of these studies will be the selection of a sample of quiescent and star-forming galaxies, as is often done in the literature by using well known colour techniques such as the `UVJ' diagram. However, given the limited number of filters available for the Euclid telescope, the recovery of such rest-frame colours will be challenging. We therefore investigate the use of observed Euclid colours, on their own and together with ground-based u-band observations, for selecting quiescent and star-forming galaxies. The most efficient colour combination, among the ones tested in this work, consists of the (u-VIS) and (VIS-J) colours. We find that this combination allows users to select a sample of quiescent galaxies complete to above $\sim70\%$ and with less than 15$\%$ contamination at redshifts in the range $0.75<z<1$. For galaxies at high-z or without the u-band complementary observations, the (VIS-Y) and (J-H) colours represent a valid alternative, with $>65\%$ completeness level and contamination below 20$\%$ at $1<z<2$ for finding quiescent galaxies. In comparison, the sample of quiescent galaxies selected with the traditional UVJ technique is only $\sim20\%$ complete at $z<3$, when recovering the rest-frame colours using mock Euclid observations. This shows that our new methodology is the most suitable one when only Euclid bands, along with u-band imaging, are available.

We revisit the dark degeneracy that arises from the Einstein equations relating geometry to the total cosmic substratum but not resolving its individual components separately. We establish the explicit conditions for the dark degeneracy in the fluid description of the dark sector. At the background level, this degeneracy can be formally understood in terms of a unified dark sector Equation of State (EoS) that depends both on the dynamical nature of the dark energy (DE) as well as on its interaction with the pressureless dark matter. For linear perturbations, the degeneracy arises for specified DE pressure perturbations (or sound speed, equivalently) and DE anisotropic stress. Specializing to the degeneracy between non-interacting dynamical DE and interacting vacuum DE models, we perform a parameter estimation analysis for a range of dynamical DE parametrizations, where for illustration we explicitly break the degeneracy at the linear level by adopting a luminal sound speed for both scenarios. We conduct this analysis using cosmological background data alone and in combination with Planck measurements of the cosmic microwave background radiation. We find that although the overall phenomenology between the dynamical DE and interacting approaches is similar, there are some intriguing differences. In particular, there is an ambiguity in the strength of constraints on $\Omega_{m0}$ and $\sigma_8$, which are considerably weakened for interacting vacuum DE, indicating that the dark degeneracy can change the significance of tensions in cosmological parameters inferred from different data sets.

The Universe is not perfectly homogeneous, the large scale structure forms overdense regions and voids. In this paper, we consider the possibility that we occupy a special position in our Universe, close to the center of a local underdense region that we model as an LTB void embedded in a homogeneous and isotropic Universe. The CMB sky measured by an off-center observer in this void is not statistically isotropic. In addition to the non-stochastic CMB anisotropies due to the geometry of the model, we also observe a lensing-like distortion of the CMB anisotropies. In this article, we propose a framework to forecast the precision with which we can measure the amplitude of the lensing-like deformation of the CMB temperature anisotropies. For illustrative purposes, we apply this method to a couple of large-scale void models differing over the matter density profile and we show that the CMB temperature data from the Planck satellite is potentially capable to detect the effect for the large voids chosen here. A companion paper will be dedicated to a systematic exploration of different realistic void models.

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.

We provide predictions of the yield of $7<z<9$ quasars from the Euclid wide survey, updating the calculation presented in the Euclid Red Book in several ways. We account for revisions to the Euclid near-infrared filter wavelengths; we adopt steeper rates of decline of the quasar luminosity function (QLF; $\Phi$) with redshift, $\Phi\propto10^{k(z-6)}$, $k=-0.72$, and a further steeper rate of decline, $k=-0.92$; we use better models of the contaminating populations (MLT dwarfs and compact early-type galaxies); and we use an improved Bayesian selection method, compared to the colour cuts used for the Red Book calculation, allowing the identification of fainter quasars, down to $J_{AB}\sim23$. Quasars at $z>8$ may be selected from Euclid $OYJH$ photometry alone, but selection over the redshift interval $7<z<8$ is greatly improved by the addition of $z$-band data from, e.g., Pan-STARRS and LSST. We calculate predicted quasar yields for the assumed values of the rate of decline of the QLF beyond $z=6$. For the case that the decline of the QLF accelerates beyond $z=6$, with $k=-0.92$, Euclid should nevertheless find over 100 quasars with $7.0<z<7.5$, and $\sim25$ quasars beyond the current record of $z=7.5$, including $\sim8$ beyond $z=8.0$. The first Euclid quasars at $z>7.5$ should be found in the DR1 data release, expected in 2024. It will be possible to determine the bright-end slope of the QLF, $7<z<8$, $M_{1450}<-25$, using 8m class telescopes to confirm candidates, but follow-up with JWST or E-ELT will be required to measure the faint-end slope. Contamination of the candidate lists is predicted to be modest even at $J_{AB}\sim23$. The precision with which $k$ can be determined over $7<z<8$ depends on the value of $k$, but assuming $k=-0.72$ it can be measured to a 1 sigma uncertainty of 0.07.

Our aim is to quantify the impact of systematic effects on the inference of cosmological parameters from cosmic shear. We present an end-to-end approach that introduces sources of bias in a modelled weak lensing survey on a galaxy-by-galaxy level. Residual biases are propagated through a pipeline from galaxy properties (one end) through to cosmic shear power spectra and cosmological parameter estimates (the other end), to quantify how imperfect knowledge of the pipeline changes the maximum likelihood values of dark energy parameters. We quantify the impact of an imperfect correction for charge transfer inefficiency (CTI) and modelling uncertainties of the point spread function (PSF) for Euclid, and find that the biases introduced can be corrected to acceptable levels.

The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. Estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on different methodologies and numerical implementations, developed for different observational probes and for their combination. In this paper we present validated forecasts, that combine both theoretical and observational expertise for different cosmological probes. This is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. We describe in detail the methodology adopted for Fisher matrix forecasts, applied to galaxy clustering, weak lensing and their combination. We estimate the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations that can be used by the reader to validate their own codes if required. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setup, i.e. flat or non-flat spatial cosmologies, or different cuts at nonlinear scales. The validated numerical implementations can now be reliably used for any setup. We present results for an optimistic and a pessimistic choice of such settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy Figure of Merit by at least a factor of three.

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)

Anomaly detection is challenging, especially for large datasets in high dimensions. Here we explore a general anomaly detection framework based on dimensionality reduction and unsupervised clustering. We release DRAMA, a general python package that implements the general framework with a wide range of built-in options. We test DRAMA on a wide variety of simulated and real datasets, in up to 3000 dimensions, and find it robust and highly competitive with commonly-used anomaly detection algorithms, especially in high dimensions. The flexibility of the DRAMA framework allows for significant optimization once some examples of anomalies are available, making it ideal for online anomaly detection, active learning and highly unbalanced datasets.

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.

We propose a new tomographic estimator for the gravitational lensing potential, based on a combination of intensity mapping (IM) and galaxy number counts. The estimator can be written schematically as IM$\times$galaxy $-$ galaxy$\times$IM; this combination allows to greatly reduce the contamination by density-density correlations, thus isolating the lensing signal. As a pure cross-correlation estimator, it is additionally less susceptible to systematic effects. We show that the new estimator strongly suppresses cosmic variance and consequently improves the signal-to-noise ratio (SNR) for the detection of lensing, especially on linear scales and intermediate redshifts. For cosmic variance dominated surveys the SNR of our estimator is a factor 30 larger than the SNR obtained from the correlation of galaxy number counts only. Shot noise and interferometer noise reduce the SNR. For the specific example of DES cross-correlated with HIRAX, the SNR is around 4, whereas for Euclid cross-correlated with HIRAX it reaches 52. This corresponds to an improvement of a factor 4-5 compared to the SNR from DES alone. For Euclid cross-correlated with HIRAX the improvement with respect to Euclid alone strongly depends on the redshift. We find that the improvement is particularly important for redshifts below 1.6, where it reaches 5. This makes our estimator especially valuable to test dark energy and modified gravity, that are expected to leave an impact at low and intermediate redshifts.

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)

Future galaxy clustering surveys will probe small scales where non-linearities become important. Since the number of modes accessible on intermediate to small scales is very high, having a precise model at these scales is important especially in the context of discriminating alternative cosmological models from the standard one. In the mildly non-linear regime, such models typically differ from each other, and galaxy clustering data will become very precise on these scales in the near future. As the observable quantity is the angular power spectrum in redshift space, it is important to study the effects of non-linear density and redshift space distortion (RSD) in the angular power spectrum. We compute non-linear contributions to the angular power spectrum using a flat-sky approximation that we introduce in this work, and compare the results of different perturbative approaches with $N$-body simulations. We find that the TNS perturbative approach is significantly closer to the $N$-body result than Eulerian or Lagrangian 1-loop approximations, effective field theory of large scale structure or a halofit-inspired model. However, none of these prescriptions is accurate enough to model the angular power spectrum well into the non-linear regime. In addition, for narrow redshift bins, $\Delta z \lesssim 0.01$, the angular power spectrum acquires non-linear contributions on all scales, right down to $\ell=2$, and is hence not a reliable tool at this time. To overcome this problem, we need to model non-linear RSD terms, for example as TNS does, but for a matter power spectrum that remains reasonably accurate well into the deeply non-linear regime, such as halofit.

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

This is a science white paper submitted to the Astro-2020 and Plasma-2020 Decadal Surveys. The paper describes the present status and emerging opportunities in Extreme Plasma Astrophysics -- a study of astrophysically-relevant plasma processes taking place under extreme conditions that necessitate taking into account relativistic, radiation, and QED effects.

An important part of cosmological model fitting relies on correlating distance indicators of objects (for example type Ia supernovae) with their redshift, often illustrated on a Hubble diagram. Comparing the observed correlation with a homogeneous model is one of the key pieces of evidence for dark energy. The presence of cosmic structures introduces a bias and scatter, mainly due to gravitational lensing and peculiar velocities, but also due to smaller non-linear relativistic contributions which are more difficult to account for. For the first time we perform ray tracing onto halos in a relativistic N-body simulation. Our simulation is the largest that takes into account all leading-order corrections from general relativity in the evolution of structure, and we present a novel methodology for working out the non-linear projection of that structure onto the observer's past light cone. We show that the mean of the bias in the Hubble diagram is indeed as small as expected from perturbation theory. However, the distribution of sources is significantly skewed with a very long tail of highly magnified objects and we illustrate that the bias of cosmological parameters strongly depends on the function of distance which we consider.

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 present three-dimensional, hybrid-kinetic numerical simulations of driven Alfv\'{e}n-wave turbulence of relevance to the collisionless near-Earth solar wind. Special attention is paid to the spectral transition that occurs near the ion-Larmor scale and to the origins of preferential perpendicular ion heating and of non-thermal wings in the parallel distribution function. Several novel diagnostics are used to show that the ion heating rate increases as the kinetic-Alfv\'{e}n-wave fluctuations, which comprise the majority of the sub-ion-Larmor turbulent cascade, attain near-ion-cyclotron frequencies. We find that ${\approx}75$-$80\%$ of the cascade energy goes into heating the ions, broadly consistent with the near-Earth solar wind. This heating is accompanied by clear velocity-space signatures in the particle energization rates and the distribution functions, including a flattened core in the perpendicular-velocity distribution and non-Maxwellian wings in the parallel-velocity distribution. The latter are attributed to transit-time damping and the pitch-angle scattering of perpendicularly heated particles into the parallel direction. Accompanying these features is a steepening of the spectral index of sub-ion-Larmor magnetic-field fluctuations beyond the canonical $-2.8$, as field energy is transferred to thermal energy. These predictions may be tested by measurements in the near-Earth solar wind.

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

We present results from the first simulations of networks of Type I Abelian Higgs cosmic strings to include both matter and radiation eras and Cosmic Microwave Background (CMB) constraints. In Type I strings, the string tension is a slowly decreasing function of the ratio of the scalar and gauge mass-squared, $\beta$. We find that the mean string separation shows no dependence on $\beta$, and that the energy-momentum tensor correlators decrease approximately in proportion to the square of the string tension, with additional O(1) correction factors which asymptote to constants below $\beta \lesssim 0.01$. Strings in models with low self-couplings can therefore satisfy current CMB bounds at higher symmetry-breaking scales. This is particularly relevant for models where the gauge symmetry is broken in a supersymmetric flat direction, for which the effective self-coupling can be extremely small. If our results can be extrapolated to $\beta \simeq 10^{-15}$, even strings formed at $10^{16}$ GeV (approximately the grand unification scale in supersymmetric extensions of the Standard Model) can be compatible with CMB constraints.

We present a systematic shearing-box investigation of MRI-driven turbulence in a weakly collisional plasma by including the effects of an anisotropic pressure stress, i.e. anisotropic (Braginskii) viscosity. We constrain the pressure anisotropy ($\Delta p$) to lie within the stability bounds that would be otherwise imposed by kinetic microinstabilities. We explore a broad region of parameter space by considering different Reynolds numbers and magnetic-field configurations, including net vertical flux, net toroidal-vertical flux and zero net flux. Remarkably, we find that the level of turbulence and angular-momentum transport are not greatly affected by large anisotropic viscosities: the Maxwell and Reynolds stresses do not differ much from the MHD result. Angular-momentum transport in Braginskii MHD still depends strongly on isotropic dissipation, e.g., the isotropic magnetic Prandtl number, even when the anisotropic viscosity is orders of magnitude larger than the isotropic diffusivities. Braginskii viscosity nevertheless changes the flow structure, rearranging the turbulence to largely counter the parallel rate of strain from the background shear. We also show that the volume-averaged pressure anisotropy and anisotropic viscous transport decrease with increasing isotropic Reynolds number (${\rm Re}$); e.g., in simulations with net vertical field, the ratio of anisotropic to Maxwell stress ($\alpha_{\rm A} / \alpha_{\rm M}$) decreases from $\sim 0.5$ to $\sim 0.1$ as we move from ${\rm Re} \sim 10^3$ to ${\rm Re} \sim 10^4$. Anisotropic transport may thus become negligible at high ${\rm Re}$. Anisotropic viscosity nevertheless becomes the dominant source of heating at large ${\rm Re}$, accounting for $\gtrsim 50 \%$ of the plasma heating. We conclude by briefly discussing the implications of our results for RIAFs onto black holes.

Galaxy clusters are the largest and most massive bound objects resulting from cosmic hierarchical structure formation. Baryons account for somewhat more than 10% of that mass, with roughly 90% of the baryonic matter distributed throughout the clusters as hot ($T>1$ keV), high-$\beta$, very weakly collisional plasma; the so-called "intracluster medium" (ICM). Cluster mergers, close gravitational encounters and accretion, along with violent feedback from galaxies and relativistic jets from active galactic nuclei, drive winds, gravity waves, turbulence and shocks within the ICM. Those dynamics, in turn, generate cluster-scale magnetic fields and accelerate and mediate the transport of high-energy charged particles. Kinetic-scale, collective plasma processes define the basic character and fundamental signatures of these ICM phenomena, which are observed primarily by X-ray and radio astronomers.

Galaxy clusters are massive dark matter-dominated systems filled with X-ray emitting, optically thin plasma. Their large size and relative simplicity (at least as astrophysical objects go) make them a unique laboratory to measure some of the interesting plasma properties that are inaccessible by other means but fundamentally important for understanding and modeling many astrophysical phenomena -- from solar flares to black hole accretion to galaxy formation and the emergence of the cosmological Large Scale Structure. While every cluster astrophysicist is eagerly anticipating the direct gas velocity measurements from the forthcoming microcalorimeters onboard XRISM, Athena and future missions such as Lynx, a number of those plasma properties can best be probed by high-resolution X-ray imaging of galaxy clusters. Chandra has obtained some trailblazing results, but only grazed the surface of such studies. In this white paper, we discuss why we need arcsecond-resolution, high collecting area, low relative background X-ray imagers (with modest spectral resolution), such as the proposed AXIS and the imaging detector of Lynx.

This white paper submitted for 2020 Decadal Assessment of Plasma Science concerns the importance of multi-spacecraft missions to address fundamental questions concerning plasma turbulence. Plasma turbulence is ubiquitous in the universe, and it is responsible for the transport of mass, momentum, and energy in such diverse systems as the solar corona and wind, accretion discs, planet formation, and laboratory fusion devices. Turbulence is an inherently multi-scale and multi-process phenomenon, coupling the largest scales of a system to sub-electron scales via a cascade of energy, while simultaneously generating reconnecting current layers, shocks, and a myriad of instabilities and waves. The solar wind is humankind's best resource for studying the naturally occurring turbulent plasmas that permeate the universe. Since launching our first major scientific spacecraft mission, Explorer 1, in 1958, we have made significant progress characterizing solar wind turbulence. Yet, due to the severe limitations imposed by single point measurements, we are unable to characterize sufficiently the spatial and temporal properties of the solar wind, leaving many fundamental questions about plasma turbulence unanswered. Therefore, the time has now come wherein making significant additional progress to determine the dynamical nature of solar wind turbulence requires multi-spacecraft missions spanning a wide range of scales simultaneously. A dedicated multi-spacecraft mission concurrently covering a wide range of scales in the solar wind would not only allow us to directly determine the spatial and temporal structure of plasma turbulence, but it would also mitigate the limitations that current multi-spacecraft missions face, such as non-ideal orbits for observing solar wind turbulence. Some of the fundamentally important questions that can only be addressed by in situ multipoint measurements are discussed.

This white paper, submitted for the Plasma 2020 Decadal Survey, concerns the physics of weakly collisional, high-beta plasmas -- plasmas in which the thermal pressure dominates over the magnetic pressure and in which the inter-particle collision time is comparable to the characteristic timescales of bulk motions. This state of matter, although widespread in the Universe, remains poorly understood: we lack a predictive theory for how it responds to perturbations, how it transports momentum and energy, and how it generates and amplifies magnetic fields. Such topics are foundational to the scientific study of plasmas, and are of intrinsic interest to those who regard plasma physics as a fundamental physics discipline. But these topics are also of extrinsic interest: addressing them directly informs upon our understanding of a wide variety of space and astrophysical systems, including accretion flows around supermassive black holes, the intracluster medium (ICM) between galaxies in clusters, and regions of the near-Earth solar wind. Specific recommendations to advance this field of study are discussed.

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.

We propose that pressure anisotropy causes weakly collisional turbulent plasmas to self-organize so as to resist changes in magnetic-field strength. We term this effect "magneto-immutability" by analogy with incompressibility (resistance to changes in pressure). The effect is important when the pressure anisotropy becomes comparable to the magnetic pressure, suggesting that in collisionless, weakly magnetized (high-$\beta$) plasmas its dynamical relevance is similar to that of incompressibility. Simulations of magnetized turbulence using the weakly collisional Braginskii model show that magneto-immutable turbulence is surprisingly similar, in most statistical measures, to critically balanced MHD turbulence. However, in order to minimize magnetic-field variation, the flow direction becomes more constrained than in MHD, and the turbulence is more strongly dominated by magnetic energy (a nonzero "residual energy"). These effects represent key differences between pressure-anisotropic and fluid turbulence, and should be observable in the $\beta\gtrsim1$ turbulent solar wind.

In a magnetized, collisionless plasma, the magnetic moment of the constituent particles is an adiabatic invariant. An increase in the magnetic-field strength in such a plasma thus leads to an increase in the thermal pressure perpendicular to the field lines. Above a $\beta$-dependent threshold (where $\beta$ is the ratio of thermal to magnetic pressure), this pressure anisotropy drives the mirror instability, producing strong distortions in the field lines on ion-Larmor scales. The impact of this instability on magnetic reconnection is investigated using a simple analytical model for the formation of a current sheet (CS) and the associated production of pressure anisotropy. The difficulty in maintaining an isotropic, Maxwellian particle distribution during the formation and subsequent thinning of a CS in a collisionless plasma, coupled with the low threshold for the mirror instability in a high-$\beta$ plasma, imply that the geometry of reconnecting magnetic fields can differ radically from the standard Harris-sheet profile often used in simulations of collisionless reconnection. As a result, depending on the rate of CS formation and the initial CS thickness, tearing modes whose growth rates and wavenumbers are boosted by this difference may disrupt the mirror-infested CS before standard tearing modes can develop. A quantitative theory is developed to illustrate this process, which may find application in the tearing-mediated disruption of kinetic magnetorotational "channel" modes.

The recent discovery of a $\gamma$-ray counterpart to a gravitational wave event has put extremely stringent constraints on the speed of gravitational waves at the present epoch. In turn, these constraints place strong theoretical pressure on potential modifications of gravity, essentially allowing only a conformally-coupled scalar to be active in the present Universe. In this paper, we show that direct detection of gravitational waves from optically identified sources can also measure or constrain the strength of the conformal coupling in scalar--tensor models through the time variation of the Planck mass. As a first rough estimate, we find that the LISA satellite can measure the dimensionless time variation of the Planck mass (the so-called parameter $\alpha_M$) at redshift around 1.5 with an error of about 0.03 to 0.13, depending on the assumptions concerning future observations. Stronger constraints can be achieved once reliable distance indicators at $z>2$ are developed, or with GW detectors that extend the capabilities of LISA, like the proposed Big Bang Observer. We emphasize that, just like the constraints on the gravitational speed, the bound on $\alpha_M$ is independent of the cosmological model.

The goal of this short report is to summarise some key results based on our previous works on model independent tests of gravity at large scales in the Universe, their connection with the properties of gravitational waves, and the implications of the recent measurement of the speed of tensors for the phenomenology of general families of gravity models for dark energy.

A persistent theme in the study of dark energy is the question of whether it really exists or not. It is often claimed hat we are mis-calculating the cosmological model by neglecting the effects associated with averaging over large-scale structures. In the Newtonian approximation this is clear: there is no effect. Within the full relativistic picture this remains an important open question however, owing to the complex mathematics involved. We study this issue using particle numerical simulations which account for all relevant relativistic effects without any problems from shell crossing. In this context we show for the first time that the backreaction from structure can differ by many orders of magnitude depending upon the slicing of spacetime one chooses to average over. In the worst case, where smoothing is carried out in synchronous spatial surfaces, the corrections can reach ten percent and more. However, when smoothing on the constant time hypersurface of the Newtonian gauge, backreaction contributions remain 3-5 orders of magnitude smaller.

The study of polarized dust emission has become entwined with the analysis of the cosmic microwave background (CMB) polarization. We use new Planck maps to characterize Galactic dust emission as a foreground to the CMB polarization. We present Planck EE, BB, and TE power spectra of dust polarization at 353 GHz for six nested sky regions covering from 24 to 71 % of the sky. We present power-law fits to the angular power spectra, yielding evidence for statistically significant variations of the exponents over sky regions and a difference between the values for the EE and BB spectra. The TE correlation and E/B power asymmetry extend to low multipoles that were not included in earlier Planck polarization papers. We also report evidence for a positive TB dust signal. Combining data from Planck and WMAP, we determine the amplitudes and spectral energy distributions (SEDs) of polarized foregrounds, including the correlation between dust and synchrotron polarized emission, for the six sky regions as a function of multipole. This quantifies the challenge of the component separation procedure required for detecting the reionization and recombination peaks of primordial CMB B modes. The SED of polarized dust emission is fit well by a single-temperature modified blackbody emission law from 353 GHz to below 70 GHz. For a dust temperature of 19.6 K, the mean spectral index for dust polarization is $\beta_{\rm d}^{P} = 1.53\pm0.02 $. By fitting multi-frequency cross-spectra, we examine the correlation of the dust polarization maps across frequency. We find no evidence for decorrelation. If the Planck limit for the largest sky region applies to the smaller sky regions observed by sub-orbital experiments, then decorrelation might not be a problem for CMB experiments aiming at a primordial B-mode detection limit on the tensor-to-scalar ratio $r\simeq0.01$ at the recombination peak.

This paper presents the Planck Multi-frequency Catalogue of Non-thermal (i.e. synchrotron-dominated) Sources (PCNT) observed between 30 and 857 GHz by the ESA Planck mission. This catalogue was constructed by selecting objects detected in the full mission all-sky temperature maps at 30 and 143 GHz, with a signal-to-noise ratio (S/N)>3 in at least one of the two channels after filtering with a particular Mexican hat wavelet. As a result, 29400 source candidates were selected. Then, a multi-frequency analysis was performed using the Matrix Filters methodology at the position of these objects, and flux densities and errors were calculated for all of them in the nine Planck channels. The present catalogue is the first unbiased, full-sky catalogue of synchrotron-dominated sources published at millimetre and submillimetre wavelengths and constitutes a powerful database for statistical studies of non-thermal extragalactic sources, whose emission is dominated by the central active galactic nucleus. Together with the full multi-frequency catalogue, we also define the Bright Planck Multi-frequency Catalogue of Non-thermal Sources PCNTb, where only those objects with a S/N>4 at both 30 and 143 GHz were selected. In this catalogue 1146 compact sources are detected outside the adopted Planck GAL070 mask; thus, these sources constitute a highly reliable sample of extragalactic radio sources. We also flag the high-significance subsample PCNThs, a subset of 151 sources that are detected with S/N>4 in all nine Planck channels, 75 of which are found outside the Planck mask adopted here. The remaining 76 sources inside the Galactic mask are very likely Galactic objects.

We present a final description of the data-processing pipeline for the Planck, Low Frequency Instrument (LFI), implemented for the 2018 data release. Several improvements have been made with respect to the previous release, especially in the calibration process and in the correction of instrumental features such as the effects of nonlinearity in the response of the analogue-to-digital converters. We provide a brief pedagogical introduction to the complete pipeline, as well as a detailed description of the important changes implemented. Self-consistency of the pipeline is demonstrated using dedicated simulations and null tests. We present the final version of the LFI full sky maps at 30, 44, and 70 GHz, both in temperature and polarization, together with a refined estimate of the Solar dipole and a final assessment of the main LFI instrumental parameters.

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.

Results from a numerical study of fluctuation dynamo in a collisionless, weakly magnetized plasma are presented. The key difference between this dynamo and its magnetohydrodynamic (MHD) counterpart is the adiabatic production of magnetic-field-aligned pressure anisotropy by the amplification of a weak seed field. This in turn drives kinetic instabilities on the ion-Larmor scale---namely, firehose and mirror---which sever the adiabatic link between the thermal and magnetic pressures, thereby allowing the dynamo to proceed. After an initial phase of rapid growth driven by these instabilities, the magnetic energy grows exponentially and exhibits a $k^{3/2}$ spectrum that peaks near the resistive scale, similar to the large-magnetic-Prandtl-number ($\mathrm{Pm}\gg{1}$) MHD dynamo. The magnetic field self-organizes into a folded-sheet topology, with direction reversals at the resistive scale and field lines curved at the parallel scale of the flow. The effective $\mathrm{Pm}$ is determined by whether the ion-Larmor scale is above or below the field-reversing scale: in the former case, particles undergo Bohm-like diffusion; in the latter case, particles scatter primarily off firehose fluctuations residing at the ends of the magnetic folds, and the viscosity becomes anisotropic. The magnetic field ultimately saturates at dynamical strengths, with its spectral peak migrating towards larger scales. This feature, along with an anti-correlation of magnetic-field strength and field-line curvature and a gradual thinning of magnetic sheets into ribbons, resemble the saturated state of the large-$\mathrm{Pm}$ dynamo, the primary differences manifesting in firehose/mirror-unstable regions. These results have implications for magnetic-field growth in the weakly collisional intracluster medium of galaxy clusters.

We present full-sky maps of the cosmic microwave background (CMB) and polarized synchrotron and thermal dust emission, derived from the third set of Planck frequency maps. These products have significantly lower contamination from instrumental systematic effects than previous versions. The methodologies used to derive these maps follow closely those described in earlier papers, adopting four methods (Commander, NILC, SEVEM, and SMICA) to extract the CMB component, as well as three methods (Commander, GNILC, and SMICA) to extract astrophysical components. Our revised CMB temperature maps agree with corresponding products in the Planck 2015 delivery, whereas the polarization maps exhibit significantly lower large-scale power, reflecting the improved data processing described in companion papers; however, the noise properties of the resulting data products are complicated, and the best available end-to-end simulations exhibit relative biases with respect to the data at the few percent level. Using these maps, we are for the first time able to fit the spectral index of thermal dust independently over 3 degree regions. We derive a conservative estimate of the mean spectral index of polarized thermal dust emission of beta_d = 1.55 +/- 0.05, where the uncertainty marginalizes both over all known systematic uncertainties and different estimation techniques. For polarized synchrotron emission, we find a mean spectral index of beta_s = -3.1 +/- 0.1, consistent with previously reported measurements. We note that the current data processing does not allow for construction of unbiased single-bolometer maps, and this limits our ability to extract CO emission and correlated components. The foreground results for intensity derived in this paper therefore do not supersede corresponding Planck 2015 products. For polarization the new results supersede the corresponding 2015 products in all respects.

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.

In this paper, we study lensing of 21cm intensity mapping (IM). Like in the cosmic microwave background (CMB), there is no first order lensing in intensity mapping. The first effects in the power spectrum are therefore of second and third order. Despite this, lensing of the CMB power spectrum is an important effect that needs to be taken into account, which motivates the study of the impact of lensing on the IM power spectrum. We derive a general formula up to third order in perturbation theory including all the terms with two derivatives of the gravitational potential, i.e. the dominant terms on sub-Hubble scales. We then show that in intensity mapping there is a new lensing term which is not present in the CMB. We obtain that the signal-to-noise of 21 cm lensing for futuristic surveys like SKA2 is about 10. We find that surveys probing only large scales, lmax < 700, can safely neglect the lensing of the intensity mapping power spectrum, but that otherwise this effect should be included.

We study the impact on the galaxy correlation function of the presence of a vector component in the tracers' peculiar velocities, in the case in which statistical isotropy is violated. We present a general framework - based on the bipolar spherical harmonics expansion - to study this effect in a model independent way, without any hypothesis on the origin or the properties of these vector modes. We construct six new observables, that can be directly measured in galaxy catalogs in addition to the standard monopole, quadrupole and hexadecapole, and we show that they completely describe any deviations from isotropy. We then perform a Fisher analysis in order to quantify the constraining power of future galaxy surveys. As an example, we show that the SKA2 would be able to detect anisotropic rotational velocities with amplitudes as low as 1% of that of the vorticity generated during shell-crossing in standard dark matter scenarios.

We give an analytical interpretation of how subsample-based internal covariance estimators lead to biased estimates of the covariance, due to underestimating the super-sample covariance (SSC). This includes the jackknife and bootstrap methods as estimators for the full survey area, and subsampling as an estimator of the covariance of subsamples. The limitations of the jackknife covariance have been previously presented in the literature because it is effectively a rescaling of the covariance of the subsample area. However we point out that subsampling is also biased, but for a different reason: the subsamples are not independent, and the corresponding lack of power results in SSC underprediction. We develop the formalism in the case of cluster counts that allows the bias of each covariance estimator to be exactly predicted. We find significant effects for a small-scale area or when a low number of subsamples is used, with auto-redshift biases ranging from 0.4% to 15% for subsampling and from 5% to 75% for jackknife covariance estimates. The cross-redshift covariance is even more affected; biases range from 8% to 25% for subsampling and from 50% to 90% for jackknife. Owing to the redshift evolution of the probe, the covariances cannot be debiased by a simple rescaling factor, and an exact debiasing has the same requirements as the full SSC prediction. These results thus disfavour the use of internal covariance estimators on data itself or a single simulation, leaving analytical prediction and simulations suites as possible SSC predictors.

We present a theoretical framework for describing electromagnetic kinetic turbulence in a multi-species, magnetized, pressure-anisotropic plasma. Turbulent fluctuations are assumed to be small compared to the mean field, to be spatially anisotropic with respect to it, and to have frequencies small compared to the ion cyclotron frequency. At scales above the ion Larmor radius, the theory reduces to the pressure-anisotropic generalization of kinetic reduced magnetohydrodynamics (KRMHD) formulated by Kunz et al. (2015). At scales at and below the ion Larmor radius, three main objectives are achieved. First, we analyse the linear response of the pressure-anisotropic gyrokinetic system, and show it to be a generalisation of previously explored limits. The effects of pressure anisotropy on the stability and collisionless damping of Alfvenic and compressive fluctuations are highlighted, with attention paid to the spectral location and width of the frequency jump that occurs as Alfven waves transition into kinetic Alfven waves. Secondly, we derive and discuss a general free-energy conservation law, which captures both the KRMHD free-energy conservation at long wavelengths and dual cascades of kinetic Alfven waves and ion entropy at sub-ion-Larmor scales. We show that non-Maxwellian features in the distribution function change the amount of phase mixing and the efficiency of magnetic stresses, and thus influence the partitioning of free energy amongst the cascade channels. Thirdly, a simple model is used to show that pressure anisotropy can cause large variations in the ion-to-electron heating ratio due to the dissipation of Alfvenic turbulence. Our theory provides a foundation for determining how pressure anisotropy affects the turbulent fluctuation spectra, the differential heating of particle species, and the ratio of parallel and perpendicular phase mixing in space and astrophysical plasmas.

To explain energy dissipation via turbulence in collisionless, magnetized plasmas, the existence of a dual real- and velocity-space cascade of ion-entropy fluctuations below the ion gyroradius has been proposed. Such a dual cascade, predicted by the gyrokinetic theory, has previously been observed in gyrokinetic simulations of two-dimensional, electrostatic turbulence. For the first time we show evidence for a dual phase-space cascade of ion-entropy fluctuations in a three-dimensional simulation of hybrid-kinetic, electromagnetic turbulence. Some of the scalings observed in the energy spectra are consistent with a generalized theory for the cascade that accounts for the spectral anisotropy of critically balanced, intermittent, sub-ion-Larmor-scale fluctuations. The observed velocity-space cascade is also anisotropic with respect to the magnetic-field direction, with linear phase mixing along magnetic-field lines proceeding mainly at spatial scales above the ion gyroradius and nonlinear phase mixing across magnetic-field lines proceeding at perpendicular scales below the ion gyroradius. Such phase-space anisotropy could be sought in heliospheric and magnetospheric data of solar-wind turbulence and has far-reaching implications for the dissipation of turbulence in weakly collisional astrophysical plasmas.

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.

We explore the use of random forest and gradient boosting, two powerful tree-based machine learning algorithms, for the detection of cosmic strings in maps of the cosmic microwave background (CMB), through their unique Gott-Kaiser-Stebbins effect on the temperature anisotropies.The information in the maps is compressed into feature vectors before being passed to the learning units. The feature vectors contain various statistical measures of processed CMB maps that boost the cosmic string detectability. Our proposed classifiers, after training, give results improved over or similar to the claimed detectability levels of the existing methods for string tension, $G\mu$. They can make $3\sigma$ detection of strings with $G\mu \gtrsim 2.1\times 10^{-10}$ for noise-free, $0.9'$-resolution CMB observations. The minimum detectable tension increases to $G\mu \gtrsim 3.0\times 10^{-8}$ for a more realistic, CMB S4-like (II) strategy, still a significant improvement over the previous results.

Using two-dimensional hybrid-kinetic simulations, we explore the nonlinear "interruption" of standing and traveling shear-Alfv\'en waves in collisionless plasmas. Interruption involves a self-generated pressure anisotropy removing the restoring force of a linearly polarized Alfv\'enic perturbation, and occurs for wave amplitudes $\delta B_{\perp}/B_{0}\gtrsim \beta^{\,-1/2}$ (where $\beta$ is the ratio of thermal to magnetic pressure). We use highly elongated domains to obtain maximal scale separation between the wave and the ion gyroscale. For standing waves above the amplitude limit, we find that the large-scale magnetic field of the wave decays rapidly. The dynamics are strongly affected by the excitation of oblique firehose modes, which transition into long-lived parallel fluctuations at the ion gyroscale and cause significant particle scattering. Traveling waves are damped more slowly, but are also influenced by small-scale parallel fluctuations created by the decay of firehose modes. Our results demonstrate that collisionless plasmas cannot support linearly polarized Alfv\'en waves above $\delta B_{\perp}/B_{0}\sim \beta^{\,-1/2}$. They also provide a vivid illustration of two key aspects of low-collisionality plasma dynamics: (i) the importance of velocity-space instabilities in regulating plasma dynamics at high $\beta$, and (ii) how nonlinear collisionless processes can transfer mechanical energy directly from the largest scales into thermal energy and microscale fluctuations, without the need for a scale-by-scale turbulent cascade.

We compute a general expression for the contribution of vector perturbations to the redshift-space distortion of galaxy surveys. We show that they contribute to the same multipoles of the correlation function as scalar perturbations and should thus in principle be taken into account in data analysis. We derive constraints for next-generation surveys on the amplitude of two sources of vector perturbations, namely non-linear clustering and topological defects. While topological defects leave a very small imprint on redshift-space distortions, we show that the multipoles of the correlation function are sensitive to vorticity induced by non-linear clustering. Therefore future redshift surveys such as DESI or the SKA should be capable of measuring such vector modes, especially with the hexadecapole which appears to be the most sensitive to the presence of vorticity.

In collisionless and weakly collisional plasmas, such as hot accretion flows onto compact objects, the magnetorotational instability (MRI) can differ significantly from the standard (collisional) MRI. In particular, pressure anisotropy with respect to the local magnetic-field direction can both change the linear MRI dispersion relation and cause nonlinear modifications to the mode structure and growth rate, even when the field and flow perturbations are small. This work studies these pressure-anisotropy-induced nonlinearities in the weakly nonlinear, high-ion-beta regime, before the MRI saturates into strong turbulence. Our goal is to better understand how the saturation of the MRI in a low collisionality plasma might differ from that in the collisional regime. We focus on two key effects: (i) the direct impact of self-induced pressure-anisotropy nonlinearities on the evolution of an MRI mode, and (ii) the influence of pressure anisotropy on the "parasitic instabilities" that are suspected to cause the mode to break up into turbulence. Our main conclusions are: (i) The mirror instability regulates the pressure anisotropy in such a way that the linear MRI in a collisionless plasma is an approximate nonlinear solution once the mode amplitude becomes larger than the background field (just as in MHD). This implies that differences between the collisionless and collisional MRI become unimportant at large amplitudes. (ii) The break up of large amplitude MRI modes into turbulence via parasitic instabilities is similar in collisionless and collisional plasmas. Together, these conclusions suggest that the route to magnetorotational turbulence in a collisionless plasma may well be similar to that in a collisional plasma, as suggested by recent kinetic simulations. As a supplement to these findings, we offer guidance for the design of future kinetic simulations of magnetorotational turbulence.

The coincident detection of gravitational waves (GW) and a gamma-ray burst from a merger of neutron stars has placed an extremely stringent bound on the speed of GW. We showed previously that the presence of gravitational slip ($\eta$) in cosmology is intimately tied to modifications of GW propagation. This new constraint implies that the only remaining viable source of gravitational slip is a conformal coupling to gravity in scalar-tensor theories, while viable vector-tensor theories cannot now generate gravitational slip at all. We discuss structure formation in the remaining viable models, demonstrating that (i) the dark-matter growth rate must now be at least as fast as in GR, with the possible exception of the beyond Horndeski model. (ii) If there is any scale-dependence at all in the slip parameter, it is such that it takes the GR value at large scales. We show a consistency relation which must be violated if gravity is modified.

Some of the dark matter in the Universe is made up of massive neutrinos. Their impact on the formation of large scale structure can be used to determine their absolute mass scale from cosmology, but to this end accurate numerical simulations have to be developed. Due to their relativistic nature, neutrinos pose additional challenges when one tries to include them in N-body simulations that are traditionally based on Newtonian physics. Here we present the first numerical study of massive neutrinos that uses a fully relativistic approach. Our N-body code, gevolution, is based on a weak-field formulation of general relativity that naturally provides a self-consistent framework for relativistic particle species. This allows us to model neutrinos from first principles, without invoking any ad-hoc recipes. Our simulation suite comprises some of the largest neutrino simulations performed to date. We study the effect of massive neutrinos on the nonlinear power spectra and the halo mass function, focusing on the interesting mass range between 0.06 eV and 0.3 eV and including a case for an inverted mass hierarchy.

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.

We discuss the effects on the CMB, CIB, and thermal SZ effect due to the peculiar motion of an observer with respect to the CMB rest frame, which induces boosting effects. We investigate the scientific perspectives opened by future CMB space missions, focussing on the CORE proposal. The improvements in sensitivity offered by a mission like CORE, together with its high resolution over a wide frequency range, will provide a more accurate estimate of the CMB dipole. The extension of boosting effects to polarization and cross-correlations will enable a more robust determination of purely velocity-driven effects that are not degenerate with the intrinsic CMB dipole, allowing us to achieve a S/N ratio of 13; this improves on the Planck detection and essentially equals that of an ideal cosmic-variance-limited experiment up to a multipole l of 2000. Precise inter-frequency calibration will offer the opportunity to constrain or even detect CMB spectral distortions, particularly from the cosmological reionization, because of the frequency dependence of the dipole spectrum, without resorting to precise absolute calibration. The expected improvement with respect to COBE-FIRAS in the recovery of distortion parameters (in principle, a factor of several hundred for an ideal experiment with the CORE configuration) ranges from a factor of several up to about 50, depending on the quality of foreground removal and relative calibration. Even for 1% accuracy in both foreground removal and relative calibration at an angular scale of 1 deg, we find that dipole analyses for a mission like CORE will be able to improve the recovery of the CIB spectrum amplitude by a factor of 17 in comparison with current results based on FIRAS. In addition to the scientific potential of a mission like CORE for these analyses, synergies with other planned and ongoing projects are also discussed.

We investigate cosmic string networks in the Abelian Higgs model using data from a campaign of large-scale numerical simulations on lattices of up to $4096^3$ grid points. We observe scaling or self-similarity of the networks over a wide range of scales, and estimate the asymptotic values of the mean string separation in horizon length units $\dot{\xi}$ and of the mean square string velocity $\bar v^2$ in the continuum and large time limits. The scaling occurs because the strings lose energy into classical radiation of the scalar and gauge fields of the Abelian Higgs model. We quantify the energy loss with a dimensionless radiative efficiency parameter, and show that it does not vary significantly with lattice spacing or string separation. This implies that the radiative energy loss underlying the scaling behaviour is not a lattice artefact, and justifies the extrapolation of measured network properties to large times for computations of cosmological perturbations. We also show that the core growth method, which increases the defect core width with time to extend the dynamic range of simulations, does not introduce significant systematic error. We compare $\dot{\xi}$ and $\bar v^2$ to values measured in simulations using the Nambu-Goto approximation, finding that the latter underestimate the mean string separation by about 25%, and overestimate $\bar v^2$ by about 10%. The scaling of the string separation implies that string loops decay by the emission of massive radiation within a Hubble time in field theory simulations, in contrast to the Nambu-Goto scenario which neglects this energy loss mechanism. String loops surviving for only one Hubble time emit much less gravitational radiation than in the Nambu-Goto scenario, and are consequently subject to much weaker gravitational wave constraints on their tension.

We present an analysis of the main systematic effects that could impact the measurement of CMB polarization with the proposed CORE space mission. We employ timeline-to-map simulations to verify that the CORE instrumental set-up and scanning strategy allow us to measure sky polarization to a level of accuracy adequate to the mission science goals. We also show how the CORE observations can be processed to mitigate the level of contamination by potentially worrying systematics, including intensity-to-polarization leakage due to bandpass mismatch, asymmetric main beams, pointing errors and correlated noise. We use analysis techniques that are well validated on data from current missions such as Planck to demonstrate how the residual contamination of the measurements by these effects can be brought to a level low enough not to hamper the scientific capability of the mission, nor significantly increase the overall error budget. We also present a prototype of the CORE photometric calibration pipeline, based on that used for Planck, and discuss its robustness to systematics, showing how CORE can achieve its calibration requirements. While a fine-grained assessment of the impact of systematics requires a level of knowledge of the system that can only be achieved in a future study phase, the analysis presented here strongly suggests that the main areas of concern for the CORE mission can be addressed using existing knowledge, techniques and algorithms.

Lensing of the CMB is now a well-developed probe of large-scale clustering over a broad range of redshifts. By exploiting the non-Gaussian imprints of lensing in the polarization of the CMB, the CORE mission can produce a clean map of the lensing deflections over nearly the full-sky. The number of high-S/N modes in this map will exceed current CMB lensing maps by a factor of 40, and the measurement will be sample-variance limited on all scales where linear theory is valid. Here, we summarise this mission product and discuss the science that it will enable. For example, the summed mass of neutrinos will be determined to an accuracy of 17 meV combining CORE lensing and CMB two-point information with contemporaneous BAO measurements, three times smaller than the minimum total mass allowed by neutrino oscillations. In the search for B-mode polarization from primordial gravitational waves with CORE, lens-induced B-modes will dominate over instrument noise, limiting constraints on the gravitational wave power spectrum amplitude. With lensing reconstructed by CORE, one can "delens" the observed polarization internally, reducing the lensing B-mode power by 60%. This improves to 70% by combining lensing and CIB measurements from CORE, reducing the error on the gravitational wave amplitude by 2.5 compared to no delensing (in the null hypothesis). Lensing measurements from CORE will allow calibration of the halo masses of the 40000 galaxy clusters that it will find, with constraints dominated by the clean polarization-based estimators. CORE can accurately remove Galactic emission from CMB maps with its 19 frequency channels. We present initial findings that show that residual Galactic foreground contamination will not be a significant source of bias for lensing power spectrum measurements with CORE. [abridged]

We demonstrate that, for the baseline design of the CORE satellite mission, the polarized foregrounds can be controlled at the level required to allow the detection of the primordial cosmic microwave background (CMB) $B$-mode polarization with the desired accuracy at both reionization and recombination scales, for tensor-to-scalar ratio values of ${r\gtrsim 5\times 10^{-3}}$. We consider detailed sky simulations based on state-of-the-art CMB observations that consist of CMB polarization with $\tau=0.055$ and tensor-to-scalar values ranging from $r=10^{-2}$ to $10^{-3}$, Galactic synchrotron, and thermal dust polarization with variable spectral indices over the sky, polarized anomalous microwave emission, polarized infrared and radio sources, and gravitational lensing effects. Using both parametric and blind approaches, we perform full component separation and likelihood analysis of the simulations, allowing us to quantify both uncertainties and biases on the reconstructed primordial $B$-modes. Under the assumption of perfect control of lensing effects, CORE would measure an unbiased estimate of $r=\left(5 \pm 0.4\right)\times 10^{-3}$ after foreground cleaning. In the presence of both gravitational lensing effects and astrophysical foregrounds, the significance of the detection is lowered, with CORE achieving a $4\sigma$-measurement of $r=5\times 10^{-3}$ after foreground cleaning and $60$% delensing. For lower tensor-to-scalar ratios ($r=10^{-3}$) the overall uncertainty on $r$ is dominated by foreground residuals, not by the 40% residual of lensing cosmic variance. Moreover, the residual contribution of unprocessed polarized point-sources can be the dominant foreground contamination to primordial B-modes at this $r$ level, even on relatively large angular scales, $\ell \sim 50$. Finally, we report two sources of potential bias for the detection of the primordial $B$-modes.[abridged]

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.

We describe a space-borne, multi-band, multi-beam polarimeter aiming at a precise and accurate measurement of the polarization of the Cosmic Microwave Background. The instrument is optimized to be compatible with the strict budget requirements of a medium-size space mission within the Cosmic Vision Programme of the European Space Agency. The instrument has no moving parts, and uses arrays of diffraction-limited Kinetic Inductance Detectors to cover the frequency range from 60 GHz to 600 GHz in 19 wide bands, in the focal plane of a 1.2 m aperture telescope cooled at 40 K, allowing for an accurate extraction of the CMB signal from polarized foreground emission. The projected CMB polarization survey sensitivity of this instrument, after foregrounds removal, is 1.7 {\mu}K$\cdot$arcmin. The design is robust enough to allow, if needed, a downscoped version of the instrument covering the 100 GHz to 600 GHz range with a 0.8 m aperture telescope cooled at 85 K, with a projected CMB polarization survey sensitivity of 3.2 {\mu}K$\cdot$arcmin.

We discuss the potential of a next generation space-borne Cosmic Microwave Background (CMB) experiment for studies of extragalactic sources. Our analysis has particular bearing on the definition of the future space project, CORE, that has been submitted in response to ESA's call for a Medium-size mission opportunity as the successor of the Planck satellite. Even though the effective telescope size will be somewhat smaller than that of Planck, CORE will have a considerably better angular resolution at its highest frequencies, since, in contrast with Planck, it will be diffraction limited at all frequencies. The improved resolution implies a considerable decrease of the source confusion, i.e. substantially fainter detection limits. In particular, CORE will detect thousands of strongly lensed high-z galaxies distributed over the full sky. The extreme brightness of these galaxies will make it possible to study them, via follow-up observations, in extraordinary detail. Also, the CORE resolution matches the typical sizes of high-z galaxy proto-clusters much better than the Planck resolution, resulting in a much higher detection efficiency; these objects will be caught in an evolutionary phase beyond the reach of surveys in other wavebands. Furthermore, CORE will provide unique information on the evolution of the star formation in virialized groups and clusters of galaxies up to the highest possible redshifts. Finally, thanks to its very high sensitivity, CORE will detect the polarized emission of thousands of radio sources and, for the first time, of dusty galaxies, at mm and sub-mm wavelengths, respectively.

The six parameters of the standard $\Lambda$CDM model have best-fit values derived from the Planck temperature power spectrum that are shifted somewhat from the best-fit values derived from WMAP data. These shifts are driven by features in the Planck temperature power spectrum at angular scales that had never before been measured to cosmic-variance level precision. We investigate these shifts to determine whether they are within the range of expectation and to understand their origin in the data. Taking our parameter set to be the optical depth of the reionized intergalactic medium $\tau$, the baryon density $\omega_{\rm b}$, the matter density $\omega_{\rm m}$, the angular size of the sound horizon $\theta_*$, the spectral index of the primordial power spectrum, $n_{\rm s}$, and $A_{\rm s}e^{-2\tau}$ (where $A_{\rm s}$ is the amplitude of the primordial power spectrum), we examine the change in best-fit values between a WMAP-like large angular-scale data set (with multipole moment $\ell<800$ in the Planck temperature power spectrum) and an all angular-scale data set ($\ell<2500$ Planck temperature power spectrum), each with a prior on $\tau$ of $0.07\pm0.02$. We find that the shifts, in units of the 1$\sigma$ expected dispersion for each parameter, are $\{\Delta \tau, \Delta A_{\rm s} e^{-2\tau}, \Delta n_{\rm s}, \Delta \omega_{\rm m}, \Delta \omega_{\rm b}, \Delta \theta_*\} = \{-1.7, -2.2, 1.2, -2.0, 1.1, 0.9\}$, with a $\chi^2$ value of 8.0. We find that this $\chi^2$ value is exceeded in 15% of our simulated data sets, and that a parameter deviates by more than 2.2$\sigma$ in 9% of simulated data sets, meaning that the shifts are not unusually large. Comparing $\ell<800$ instead to $\ell>800$, or splitting at a different multipole, yields similar results. We examine the $\ell<800$ model residuals in the $\ell>800$ power spectrum data and find that the features there... [abridged]

We forecast the scientific capabilities to improve our understanding of cosmic inflation of CORE, a proposed CMB space satellite submitted in response to the ESA fifth call for a medium-size mission opportunity. The CORE satellite will map the CMB anisotropies in temperature and polarization in 19 frequency channels spanning the range 60-600 GHz. CORE will have an aggregate noise sensitivity of $1.7 \mu$K$\cdot \,$arcmin and an angular resolution of 5' at 200 GHz. We explore the impact of telescope size and noise sensitivity on the inflation science return by making forecasts for several instrumental configurations. This study assumes that the lower and higher frequency channels suffice to remove foreground contaminations and complements other related studies of component separation and systematic effects, which will be reported in other papers of the series "Exploring Cosmic Origins with CORE." We forecast the capability to determine key inflationary parameters, to lower the detection limit for the tensor-to-scalar ratio down to the $10^{-3}$ level, to chart the landscape of single field slow-roll inflationary models, to constrain the epoch of reheating, thus connecting inflation to the standard radiation-matter dominated Big Bang era, to reconstruct the primordial power spectrum, to constrain the contribution from isocurvature perturbations to the $10^{-3}$ level, to improve constraints on the cosmic string tension to a level below the presumptive GUT scale, and to improve the current measurements of primordial non-Gaussianities down to the $f_{NL}^{\rm local} < 1$ level. For all the models explored, CORE alone will improve significantly on the present constraints on the physics of inflation. Its capabilities will be further enhanced by combining with complementary future cosmological observations.

We forecast the main cosmological parameter constraints achievable with the CORE space mission which is dedicated to mapping the polarisation of the Cosmic Microwave Background (CMB). CORE was recently submitted in response to ESA's fifth call for medium-sized mission proposals (M5). Here we report the results from our pre-submission study of the impact of various instrumental options, in particular the telescope size and sensitivity level, and review the great, transformative potential of the mission as proposed. Specifically, we assess the impact on a broad range of fundamental parameters of our Universe as a function of the expected CMB characteristics, with other papers in the series focusing on controlling astrophysical and instrumental residual systematics. In this paper, we assume that only a few central CORE frequency channels are usable for our purpose, all others being devoted to the cleaning of astrophysical contaminants. On the theoretical side, we assume LCDM as our general framework and quantify the improvement provided by CORE over the current constraints from the Planck 2015 release. We also study the joint sensitivity of CORE and of future Baryon Acoustic Oscillation and Large Scale Structure experiments like DESI and Euclid. Specific constraints on the physics of inflation are presented in another paper of the series. In addition to the six parameters of the base LCDM, which describe the matter content of a spatially flat universe with adiabatic and scalar primordial fluctuations from inflation, we derive the precision achievable on parameters like those describing curvature, neutrino physics, extra light relics, primordial helium abundance, dark matter annihilation, recombination physics, variation of fundamental constants, dark energy, modified gravity, reionization and cosmic birefringence. (ABRIDGED)

In many generalized models of gravity, perfect fluids in cosmology give rise to gravitational slip. Simultaneously, in very broad classes of such models, the propagation of gravitational waves is altered. We investigate the extent to which there is a one-to-one relationship between these two properties in three classes of models with one extra degree of freedom: scalar (Horndeski and beyond), vector (Einstein-Aether) and tensor (bimetric). We prove that in bimetric gravity and Einstein-Aether, it is impossible to dynamically hide the gravitational slip on all scales whenever the propagation of gravitational waves is modified. Horndeski models are much more flexible, but it is nonetheless only possible to hide gravitational slip dynamically when the action for perturbations is tuned to evolve in time toward a divergent kinetic term. These results provide an explicit, theoretical argument for the interpretation of future observations if they disfavoured the presence of gravitational slip.

We examine the cosmological constraints that can be achieved with a galaxy cluster survey with the future CORE space mission. Using realistic simulations of the millimeter sky, produced with the latest version of the Planck Sky Model, we characterize the CORE cluster catalogues as a function of the main mission performance parameters. We pay particular attention to telescope size, key to improved angular resolution, and discuss the comparison and the complementarity of CORE with ambitious future ground-based CMB experiments that could be deployed in the next decade. A possible CORE mission concept with a 150 cm diameter primary mirror can detect of the order of 50,000 clusters through the thermal Sunyaev-Zeldovich effect (SZE). The total yield increases (decreases) by 25% when increasing (decreasing) the mirror diameter by 30 cm. The 150 cm telescope configuration will detect the most massive clusters ($>10^{14}\, M_\odot$) at redshift $z>1.5$ over the whole sky, although the exact number above this redshift is tied to the uncertain evolution of the cluster SZE flux-mass relation; assuming self-similar evolution, CORE will detect $\sim 500$ clusters at redshift $z>1.5$. This changes to 800 (200) when increasing (decreasing) the mirror size by 30 cm. CORE will be able to measure individual cluster halo masses through lensing of the cosmic microwave background anisotropies with a 1-$\sigma$ sensitivity of $4\times10^{14} M_\odot$, for a 120 cm aperture telescope, and $10^{14} M_\odot$ for a 180 cm one. [abridged]

We re-analyse recent Cepheid data to estimate the Hubble parameter $H_0$ by using Bayesian hyper-parameters (HPs). We consider the two data sets from Riess et al 2011 and 2016 (labelled R11 and R16, with R11 containing less than half the data of R16) and include the available anchor distances (megamaser system NGC4258, detached eclipsing binary distances to LMC and M31, and MW Cepheids with parallaxes), use a weak metallicity prior and no period cut for Cepheids. We find that part of the R11 data is down-weighted by the HPs but that R16 is mostly consistent with expectations for a Gaussian distribution, meaning that there is no need to down-weight the R16 data set. For R16, we find a value of $H_0 = 73.75 \pm 2.11 \, \mathrm{km} \, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ if we use HPs for all data points (including Cepheid stars, supernovae type Ia, and the available anchor distances), which is about 2.6 $\sigma$ larger than the Planck 2015 value of $H_0 = 67.81 \pm 0.92 \,\mathrm{km}\, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$ and about 3.1 $\sigma$ larger than the updated Planck 2016 value $66.93 \pm 0.62 \,\mathrm{km}\, \mathrm{s}^{-1} \, \mathrm{Mpc}^{-1}$. We test the effect of different assumptions, and find that the choice of anchor distances affects the final value significantly. If we exclude the Milky Way from the anchors, then the value of $H_0$ decreases. We find however no evident reason to exclude the MW data. The HP method used here avoids subjective rejection criteria for outliers and offers a way to test datasets for unknown systematics.

Modified Gravity theories generally affect the Poisson equation and the gravitational slip (effective anisotropic stress) in an observable way, that can be parameterized by two generic functions ($\eta$ and $\mu$) of time and space. We bin the time dependence of these functions in redshift and present forecasts on each bin for future surveys like Euclid. We consider both Galaxy Clustering and Weak Lensing surveys, showing the impact of the non-linear regime, treated with two different semi-analytical approximations. In addition to these future observables, we use a prior covariance matrix derived from the Planck observations of the Cosmic Microwave Background. Our results show that $\eta$ and $\mu$ in different redshift bins are significantly correlated, but including non-linear scales reduces or even eliminates the correlation, breaking the degeneracy between Modified Gravity parameters and the overall amplitude of the matter power spectrum. We further decorrelate parameters with a Zero-phase Component Analysis and identify which combinations of the Modified Gravity parameter amplitudes, in different redshift bins, are best constrained by future surveys. We also extend the analysis to two particular parameterizations of the time evolution of $\mu$ and $\eta$ and consider, in addition to Euclid, also SKA1, SKA2, DESI: we find in this case that future surveys will be able to constrain the current values of $\eta$ and $\mu$ at the $2-5\%$ level when using only linear scales (wavevector k < 0.15 h/Mpc), depending on the specific time parameterization; sensitivity improves to about $1\%$ when non-linearities are included.

We present cosmic microwave background (CMB) power spectra from recent numerical simulations of cosmic strings in the Abelian Higgs model and compare them to CMB power spectra measured by Planck. We obtain revised constraints on the cosmic string tension parameter $G\mu$. For example, in the $\Lambda$CDM model with the addition of strings and no primordial tensor perturbations, we find $G\mu < 2.0 \times 10^{-7}$ at 95% confidence, about 20% lower than the value obtained from previous simulations, which had 1/64 of the spatial volume. We investigate the source of the difference, showing that the main cause is an improved treatment of the string evolution across the radiation-matter transition. The increased computational volume also makes possible to simulate fully the physical equations of motion, in which the string cores shrink in comoving coordinates. This, and the larger dynamic range, changes the amplitude of the power spectra by only about 10%, demonstrating that field theory simulations of cosmic strings have now reached the required dynamic range for CMB calculations.

We report on the energy-momentum correlators obtained with recent numerical simulations of the Abelian Higgs model, essential for the computation of cosmic microwave background and matter perturbations of cosmic strings. Due to significant improvements both in raw computing power and in our parallel simulation framework, the dynamical range of the simulations has increased four-fold both in space and time, and for the first time we are able to simulate strings with a constant physical width in both the radiation and matter eras. The new simulations improve the accuracy of the measurements of the correlation functions at the horizon scale and confirm the shape around the peak. The normalization is slightly higher in the high wave-number tails, due to a small increase in the string density. We study for the first time the behaviour of the correlators across cosmological transitions, and discover that the correlation functions evolve adiabatically, ie the network adapts quickly to changes in the expansion rate. We propose a new method for constructing source functions for Einstein-Boltzmann integrators, comparing it with two other methods previously used. The new method is more consistent, easier to implement, and significantly more accurate.

Measurements of flux density are described for five planets, Mars, Jupiter, Saturn, Uranus, and Neptune, across the six Planck High Frequency Instrument frequency bands (100-857 GHz) and these are then compared with models and existing data. In our analysis, we have also included estimates of the brightness of Jupiter and Saturn at the three frequencies of the Planck Low Frequency Instrument (30, 44, and 70 GHz). The results provide constraints on the intrinsic brightness and the brightness time-variability of these planets. The majority of the planet flux density estimates are limited by systematic errors, but still yield better than 1% measurements in many cases. Applying data from Planck HFI, the Wilkinson Microwave Anisotropy Probe (WMAP), and the Atacama Cosmology Telescope (ACT) to a model that incorporates contributions from Saturn's rings to the planet's total flux density suggests a best fit value for the spectral index of Saturn's ring system of $\beta _\mathrm{ring} = 2.30\pm0.03$ over the 30-1000 GHz frequency range. The average ratio between the Planck-HFI measurements and the adopted model predictions for all five planets (excluding Jupiter observations for 353 GHz) is 0.997, 0.997, 1.018, and 1.032 for 100, 143, 217, and 353 GHz, respectively. Model predictions for planet thermodynamic temperatures are therefore consistent with the absolute calibration of Planck-HFI detectors at about the three-percent-level. We compare our measurements with published results from recent cosmic microwave background experiments. In particular, we observe that the flux densities measured by Planck HFI and WMAP agree to within 2%. These results allow experiments operating in the mm-wavelength range to cross-calibrate against Planck and improve models of radiative transport used in planetary science.

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.

We present results from the first 3D kinetic numerical simulation of magnetorotational turbulence and dynamo, using the local shearing-box model of a collisionless accretion disc. The kinetic magnetorotational instability grows from a subthermal magnetic field having zero net flux over the computational domain to generate self-sustained turbulence and outward angular-momentum transport. Significant Maxwell and Reynolds stresses are accompanied by comparable viscous stresses produced by field-aligned ion pressure anisotropy, which is regulated primarily by the mirror and ion-cyclotron instabilities through particle trapping and pitch-angle scattering. The latter endow the plasma with an effective viscosity that is biased with respect to the magnetic-field direction and spatio-temporally variable. Energy spectra suggest an Alfv\'en-wave cascade at large scales and a kinetic-Alfv\'en-wave cascade at small scales, with strong small-scale density fluctuations and weak non-axisymmetric density waves. Ions undergo non-thermal particle acceleration, their distribution accurately described by a kappa distribution. These results have implications for the properties of low-collisionality accretion flows, such as that near the black hole at the Galactic center.

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.

We outline a new method to compute the Bayes Factor for model selection which bypasses the Bayesian Evidence. Our method combines multiple models into a single, nested, Supermodel using one or more hyperparameters. Since the models are now nested the Bayes Factors between the models can be efficiently computed using the Savage-Dickey Density Ratio (SDDR). In this way model selection becomes a problem of parameter estimation. We consider two ways of constructing the supermodel in detail: one based on combined models, and a second based on combined likelihoods. We report on these two approaches for a Gaussian linear model for which the Bayesian evidence can be calculated analytically and a toy nonlinear problem. Unlike the combined model approach, where a standard Monte Carlo Markov Chain (MCMC) struggles, the combined-likelihood approach fares much better in providing a reliable estimate of the log-Bayes Factor. This scheme potentially opens the way to computationally efficient ways to compute Bayes Factors in high dimensions that exploit the good scaling properties of MCMC, as compared to methods such as nested sampling that fail for high dimensions.

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.

The stability of a collisionless, magnetised plasma to local convective disturbances is examined, with a focus on kinetic and finite-Larmor-radius effects. Specific application is made to the outskirts of galaxy clusters, which contain hot and tenuous plasma whose temperature increases in the direction of gravity. At long wavelengths (the "drift-kinetic" limit), we obtain the kinetic version of the magnetothermal instability (MTI) and its Alfv\'enic counterpart (Alfv\'enic MTI), which were previously discovered and analysed using a magnetofluid (i.e. Braginskii) description. At sub-ion-Larmor scales, we discover an overstability driven by the electron temperature gradient of kinetic-Alfv\'en drift waves -- the electron MTI (eMTI) -- whose growth rate is even larger than the standard MTI. At intermediate scales, we find that ion finite-Larmor-radius effects tend to stabilise the plasma. We discuss the physical interpretation of these instabilities in detail, and compare them both with previous work on magnetised convection in a collisional plasma and with temperature-gradient-driven drift-wave instabilities well-known to the magnetic-confinement-fusion community. The implications of having both fluid and kinetic scales simultaneously driven unstable by the same temperature gradient are briefly discussed.

We demonstrate the importance of including the lensing contribution in galaxy clustering analyses with large galaxy redshift surveys. It is well known that radial cross-correlations between different redshift bins of galaxy surveys are dominated by lensing. But we show here that also neglecting lensing in the autocorrelations within one bin severely biases cosmological parameter estimation with redshift surveys. It leads to significant shifts for several cosmological parameters, most notably the scalar spectral index and the neutrino mass scale. Especially the latter parameter is one of the main targets of future galaxy surveys.

Using the Planck 2015 data release (PR2) temperature maps, we separate Galactic thermal dust emission from cosmic infrared background (CIB) anisotropies. For this purpose, we implement a specifically tailored component-separation method, the so-called generalized needlet internal linear combination (GNILC) method, which uses spatial information (the angular power spectra) to disentangle the Galactic dust emission and CIB anisotropies. We produce significantly improved all-sky maps of Planck thermal dust emission, with reduced CIB contamination, at 353, 545, and 857 GHz. By reducing the CIB contamination of the thermal dust maps, we provide more accurate estimates of the local dust temperature and dust spectral index over the sky with reduced dispersion, especially at high Galactic latitudes above $b = \pm 20{\deg}$. We find that the dust temperature is $T = (19.4 \pm 1.3)$ K and the dust spectral index is $\beta = 1.6 \pm 0.1$ averaged over the whole sky, while $T = (19.4 \pm 1.5)$ K and $\beta = 1.6 \pm 0.2$ on 21 % of the sky at high latitudes. Moreover, subtracting the new CIB-removed thermal dust maps from the CMB-removed Planck maps gives access to the CIB anisotropies over 60 % of the sky at Galactic latitudes $|b| > 20{\deg}$. Because they are a significant improvement over previous Planck products, the GNILC maps are recommended for thermal dust science. The new CIB maps can be regarded as indirect tracers of the dark matter and they are recommended for exploring cross-correlations with lensing and large-scale structure optical surveys. The reconstructed GNILC thermal dust and CIB maps are delivered as Planck products.