Search

The current generation of air shower radio arrays has demonstrated that the atmospheric depth of the shower maximum Xmax can be reconstructed with high accuracy. These experiments are now contributing to mass composition studies in the energy range where a transition from galactic to extragalactic cosmic-ray sources is expected. However, we are still far away from an unambiguous interpretation of the data. Here we propose to use radio measurements to derive a new type of constraint on the mass composition, by reconstructing the shower length L. The low-frequency part of the Square Kilometer Array will have an extremely high antenna density of roughly 60.000 antennas within one square kilometer, and is the perfect site for high-resolution studies of air showers. In this contribution, we discuss the impact of being able to reconstruct L, and the unique contribution that SKA can make to cosmic-ray science., Comment: Proceedings 9th International Workshop on Acoustic and Radio EeV Neutrino Detection Activities - ARENA2022, 7-10 June 2022, Santiago de Compostela, Spain (8 pages)

We explore the possibilities of measuring the longitudinal profile of individual air showers beyond $X_{\rm max}$ when using very dense radio arrays such as SKA. The low-frequency part of the Square Kilometre Array, to be built in Australia, features an enormous antenna density of about $50,000$ antennas in the inner core region of radius 500 m, with a frequency band from 50 to 350 MHz. From CoREAS simulations, a SKA-Low antenna model plus noise contributions, and adapted LOFAR analysis scripts, we obtain a resolution in the shower maximum $X_{\rm max}$ and energy that is considerably better than at LOFAR. Already from this setup, we show that at least one additional parameter of the longitudinal profile can be measured. This would improve mass composition analysis by measuring an additional composition-dependent quantity. Moreover, it would offer an opportunity to discriminate between the different predictions of hadronic interaction models, hence contributing to hadronic physics at energy levels beyond man-made accelerators., Comment: 8 pages, 6 figures. Presented at ARENA Conference 2022, Santiago de Compostela, Spain

We present an updated cosmic-ray mass composition analysis in the energy range $10^{16.8}$ to $10^{18.3}$ eV from 334 air showers measured with the LOFAR radio telescope, and selected for minimal bias. In this energy range, the origin of cosmic rays is expected to shift from galactic to extragalactic sources. The analysis is based on an improved method to infer the depth of maximum $X_{\rm max}$ of extensive air showers from radio measurements and air shower simulations. We show results of the average and standard deviation of $X_{\rm max}$ versus primary energy, and analyze the $X_{\rm max}$-dataset at distribution level to estimate the cosmic ray mass composition. Our approach uses an unbinned maximum likelihood analysis, making use of existing parametrizations of $X_{\rm max}$-distributions per element. The analysis has been repeated for three main models of hadronic interactions. Results are consistent with a significant light-mass fraction, at best fit $23$ to $39$ $\%$ protons plus helium, depending on the choice of hadronic interaction model. The fraction of intermediate-mass nuclei dominates. This confirms earlier results from LOFAR, with systematic uncertainties on $X_{\rm max}$ now lowered to 7 to $9$ $\mathrm{g/cm^2}$. We find agreement in mass composition compared to results from Pierre Auger Observatory, within statistical and systematic uncertainties. However, in line with earlier LOFAR results, we find a slightly lower average $X_{\rm max}$. The values are in tension with those found at Pierre Auger Observatory, but agree with results from other cosmic ray observatories based in the Northern hemisphere., Comment: 24 pages, 14 figures. Accepted for publication in Phys. Rev. D

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution, which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1 per cent of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. The results of the BeeP extension of PCCS2, which are made publicly available via the PLA, will enable the study of the thermal properties of well-defined samples of compact Galactic and extra-galactic dusty sources., Comment: 55 pages. Accepted for publication in A&A. The BeeP catalogue will be published in the Planck Legacy Archive (https://pla.esac.esa.int/pla)

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., Comment: 97 pages, 93 figures and 16 tables, abstract abridged for arXiv submission, accepted for publication in A&A

With LOFAR we have been able to image the development of lightning flashes with meter-scale accuracy and unprecedented detail. We discuss the primary steps behind our most recent lightning mapping method. To demonstrate the capabilities of our technique we show and interpret images of the first few milliseconds of two intra-cloud flashes. In all our flashes the negative leaders propagate in the charge layer below the main negative charge. Among several interesting features we show that in about 2~ms after initiation the Primary Initial Leader triggers the formation of a multitude (more than ten) negative leaders in a rather confined area of the atmosphere. From these only one or two continue to propagate after about 30~ms to extend over kilometers horizontally while another may propagate back to the initiation point. We also show that normal negative leaders can transition into an initial-leader like state, potentially in the presence of strong electric fields. In addition, we show some initial breakdown pulses that occurred during the primary initial leader, and even during two "secondary" initial leaders that developed out of stepped leaders., Comment: Submitted to Journal of geophysics research: Atmospheres

We use the Low Frequency ARray (LOFAR) to probe the dynamics of the stepping process of negatively-charged plasma channels (negative leaders) in a lightning discharge. We observe that at each step of a leader, multiple pulses of VHF (30~--~80 MHz) radiation are emitted in short-duration bursts ($<10\ \mu$s). This is evidence for streamer formation during corona flashes that occur with each leader step, which has not been observed before in natural lightning and it could help explain X-ray emission from lightning leaders, as X-rays from laboratory leaders tend to be associated with corona flashes. Surprisingly we find that the stepping length is very similar to what was observed near the ground, however with a stepping time that is considerably larger, which as yet is not understood. These results will help to improve lightning propagation models, and eventually lightning protection models., Comment: 5 pages

The limited knowledge of atmospheric parameters like humidity, pressure, temperature, and the index of refraction has been one of the important systematic uncertainties in reconstructing the depth of the shower maximum from the radio emission of air showers. Current air shower Monte Carlo simulation codes like CORSIKA and the radio plug-in CoREAS use various averaged parameterized atmospheres. However, time-dependent and location-specific atmospheric models are needed for the cosmic ray analysis method used for LOFAR data. There, dedicated simulation sets are used for each detected cosmic ray, to take into account the actual atmospheric conditions at the time of the measurement. Using the Global Data Assimilation System (GDAS), a global atmospheric model, we have implemented time-dependent, realistic atmospheric profiles in CORSIKA and CoREAS. We have produced realistic event-specific atmospheres for all air showers measured with LOFAR, an event set spanning several years and many different weather conditions. A complete re-analysis of our data set shows that for the majority of data, our previous correction factor performed rather well; we found only a small systematic shift of 2 g/cm$^2$ in the reconstructed $X_{\rm max}$. However, under extreme weather conditions, for example, very low air pressure, the shift can be up to 15 g/cm$^2$. We provide a correction formula to determine the shift in $X_{\rm max}$ resulting from a comparison of simulations done using the US-Std atmosphere and the GDAS-based atmosphere., Comment: Accepted for publication in Astroparticle Physics. arXiv admin note: text overlap with arXiv:1911.02859

Cosmic rays are routinely measured at LOFAR, both with a dense array of antennas and with the LOFAR Radboud air shower Array (LORA) which is an array of plastic scintillators. In this paper, we present two results relating to the cosmic-ray energy scale of LOFAR. First, we present the reconstruction of cosmic-ray energy using radio and particle techniques along with a discussion of the event-by-event and absolute scale uncertainties. The resulting energies reconstructed with each method are shown to be in good agreement, and because the radio-based reconstructed energy has smaller uncertainty on an event-to-event basis, LOFAR analyses will use that technique in the future. Second, we present the radiation energy of air showers measured at LOFAR and demonstrate how radiation energy can be used to compare the energy scales of different experiments. The radiation energy scales quadratically with the electromagnetic energy in an air shower, which can in turn be related to the energy of the primary particle. Once the local magnetic field is accounted for, the radiation energy allows for a direct comparison between the LORA particle-based energy scale and that of the Pierre Auger Observatory. They are shown to agree to within (6$\pm$20)% for a radiation energy of 1 MeV, where the uncertainty on the comparison is dominated by the antenna calibrations of each experiment. This study motivates the development of a portable radio array which will be used to cross-calibrate the energy scales of different experiments using radiation energy and the same antennas, thereby significantly reducing the uncertainty on the comparison.

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, which is 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., Comment: 15 pages, 8 figures. Added references, small clarifying and language edits. All results remain the same

This is a collection of papers that have been contributed by the LOFAR Cosmic Ray Key Science Project (CRKSP) to the 36th International Cosmic Ray Conference held in Madison, Wisconsin, on July 24th to August 1st, 2019 (ICRC 2019). All papers contained here have been individually published in PoS(ICRC2019) with paper numbers 205, 362, 352, 416, and 363, in the order they appear in this collection. Minor modifications to the PoS versions have been applied where appropriate., Comment: 42 pages, 17 figures, Collection of papers presented at the ICRC2019

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), Comment: Revised to match version published in Astronomy & Astrophysics

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., Comment: Paper VII of the Planck 2018 release, revised to closely match version published in Astronomy and Astrophysics

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)., Comment: 50 pages, 20 figures

The LOFAR radio telescope is able to measure the radio emission from cosmic ray induced air showers with hundreds of individual antennas. This allows for precision testing of the emission mechanisms for the radio signal as well as determination of the depth of shower maximum $X_{\max}$, the shower observable most sensitive to the mass of the primary cosmic ray, to better than 20 g/cm$^2$. With a densely instrumented circular area of roughly 320 m$^2$, LOFAR is targeting for cosmic ray astrophysics in the energy range $10^{16}$ - $10^{18}$ eV. In this contribution we give an overview of the status, recent results, and future plans of cosmic ray detection with the LOFAR radio telescope., Comment: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

Cosmic particles hitting Earth's moon produce radio emission via the Askaryan effect. If the resulting radio ns-pulse can be detected by radio telescopes, this technique potentially increases the available collective area for ZeV scale particles by several orders of magnitude compared to current experiments. The LOw Frequency ARray (LOFAR) is the largest radio telescope operating in the optimum frequency regime for this technique. In this contribution, we report on the status of the implementation of the lunar detection mode at LOFAR., Comment: Proceedings of the 26th Extended European Cosmic Ray Symposium (ECRS), Barnaul/Belokurikha, 2018

The LOw-Frequency ARray (LOFAR) is used to make precise measurements of radio emission from extensive air showers, yielding information about the primary cosmic ray. Interpreting the measured data requires an absolute and frequency-dependent calibration of the LOFAR system response. This is particularly important for spectral analyses, because the shape of the detected signal holds information about the shower development. We revisit the calibration of the LOFAR antennas in the range of 30 - 80 MHz. Using the Galactic emission and a detailed model of the LOFAR signal chain, we find an improved calibration that provides an absolute energy scale and allows for the study of frequency-dependent features in measured signals. With the new calibration, systematic uncertainties of 13% are reached, and comparisons of the spectral shape of calibrated data with simulations show promising agreement., Comment: 23 pages, 10 figures

The advent of the Fermi-GST with its unprecedented capability to monitor the entire 4 pi sky within less than 2-3 hours, introduced new standard in time domain gamma-ray astronomy. To explore this new avenue of extragalactic physics the F-GAMMA programme undertook the task of conducting nearly monthly, broadband radio monitoring of selected blazars from January 2007 to January 2015. In this work we release all the light curves at 2.64, 4.85, 8.35, 10.45, 14.6, 23.05, 32, and 43 GHz and present first order derivative data products after all necessary post-measurement corrections and quality checks; that is flux density moments and spectral indices. The release includes 155 sources. The effective cadence after the quality flagging is around one radio SED every 1.3 months. The coherence of each radio SED is around 40 minutes. The released dataset includes more than $4\times10^4$ measurements. The median fractional error at the lowest frequencies (2.64-10.45 GHz) is below 2%. At the highest frequencies (14.6-43 GHz) with limiting factor of the atmospheric conditions, the errors range from 3% to 9%, respectively., Comment: Accepted for publication in Section: Catalogs and data of Astronomy & Astrophysics

While repeating fast radio bursts (FRBs) remain scarce in number, they provide a unique opportunity for follow-up observations that enhance our knowledge of their sources and potentially of the FRB population as a whole. Attaining more burst spectra could lead to a better understanding of the origin of these bright, millisecond-duration radio pulses. We therefore performed $\sim$20 hr of simultaneous observations on FRB 121102 with the Effelsberg 100-m radio telescope and the Low Frequency Array (LOFAR) to constrain the spectral behaviour of bursts from FRB 121102 at 1.4 GHz and 150 MHz. This campaign resulted in the detection of nine new bursts at 1.4 GHz but no simultaneous detections with LOFAR. Assuming that the ratio of the fluence at two frequencies scales as a power law, we placed a lower limit of $\alpha$ > -1.2 $\pm$ 0.4 on the spectral index for the fluence of the instantaneous broad band emission of FRB 121102. For the derivation of this limit, a realistic fluence detection threshold for LOFAR was determined empirically assuming a burst would be scattered as predicted by the NE2001 model. A significant variation was observed in the burst repeat rate R at L-band. During observations in September 2016, nine bursts were detected, giving R = 1.1 $\pm$ 0.4 hr$^{-1}$, while in November no bursts were detected, yielding R < 0.3 hr$^{-1}$ (95% confidence limit). This variation is consistent with earlier seen episodic emission of FRB 121102. In a blind and targeted search, no bursts were found with LOFAR at 150 MHz, resulting in a repeat rate limit of R < 0.16 hr$^{-1}$ (95% confidence limit). Burst repeat rate ratios of FRB 121102 at 3, 2, 1.4, and 0.15 GHz are consistent within the uncertainties with a flattening of its spectrum below 1 GHz., Comment: 8 pages, 2 tables, 1 figure

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 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., Comment: 74 pages, A&A, 641, A4

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., Comment: Accepted for publication on A&A (AA/2018/32909)

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., Comment: Accepted for publication in Astronomy & Astrophysics

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., Comment: 61 pages, 40 figures, matches version accepted by A&A

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.

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)., Comment: Abstract abridged for arxiv submission. Lensing data products available at https://wiki.cosmos.esa.int/planck-legacy-archive/index.php/Lensing. Matches version accepted by A&A, with minor updates from v1

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., Comment: References added and minor improvements. BICEP2/Keck Array BK15 is used in the place of BICEP2/Keck Array BK14

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), Comment: 73 pages; Updated with published reionization result corrigendum on p59. Parameter tables and chains available at https://wiki.cosmos.esa.int/planck-legacy-archive/index.php/Cosmological_Parameters

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., Comment: 24 pages, 15 figures. Accepted for publication in Astronomy & Astrophysics

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., Comment: Final version to appear in A&A

Contributions of the LOFAR Cosmic Ray Key Science Project to the 35th International Cosmic Ray Conference (ICRC 2017), Comment: compendium of 8 proceedings for the ICRC 2017

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., Comment: 20 pages, 12 figures and 8 tables, A&A in press

We observe a large fraction of circular polarization in radio emission from extensive air showers recorded during thunderstorms, much higher than in the emission from air showers measured during fair-weather circumstances. We show that the circular polarization of the air showers measured during thunderstorms can be explained by the change in the direction of the transverse current as a function of altitude induced by atmospheric electric fields. Thus by using the full set of Stokes parameters for these events, we obtain a good characterization of the electric fields in thunderclouds. We also measure a large horizontal component of the electric fields in the two events that we have analysed.

For the interpretation of measurements of radio emission from extensive air showers, an important systematic uncertainty arises from natural variations of the atmospheric refractive index $n$. At a given altitude, the refractivity $N=10^6\, (n-1)$ can have relative variations on the order of $10 \%$ depending on temperature, humidity, and air pressure. Typical corrections to be applied to $N$ are about $4\%$. Using CoREAS simulations of radio emission from air showers, we have evaluated the effect of varying $N$ on measurements of the depth of shower maximum $X_{\rm max}$. For an observation band of 30 to 80 MHz, a difference of $4 \%$ in refractivity gives rise to a systematic error in the inferred $X_{\rm max}$ between 3.5 and 11 $\mathrm{g/cm^2}$, for proton showers with zenith angles ranging from 15 to 50 degrees. At higher frequencies, from 120 to 250 MHz, the offset ranges from 10 to 22 $\mathrm{g/cm^2}$. These offsets were found to be proportional to the geometric distance to $X_{\rm max}$. We have compared the results to a simple model based on the Cherenkov angle. For the 120 to 250 MHz band, the model is in qualitative agreement with the simulations. In typical circumstances, we find a slight decrease in $X_{\rm max}$ compared to the default refractivity treatment in CoREAS. While this is within commonly treated systematic uncertainties, accounting for it explicitly improves the accuracy of $X_{\rm max}$ measurements., Comment: 13 pages, 5 figures. Accepted for publication in Astroparticle Physics

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., Comment: 20 pages, 14 figures, abstract abridged for arXiv submission

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. Even lower fluxes of neutrinos with energies above $10^{22}$ eV are predicted in certain Grand-Unifying-Theories (GUTs) and e.g.\ models for super-heavy dark matter (SHDM). The significant increase in detector volume required to detect these particles can be achieved by searching for the nano-second radio pulses that are emitted when a particle interacts in Earth's moon with current and future radio telescopes. In this contribution we present the design of an online analysis and trigger pipeline for the detection of nano-second pulses with the LOFAR radio telescope. The most important steps of the processing pipeline are digital focusing of the antennas towards the Moon, correction of the signal for ionospheric dispersion, and synthesis of the time-domain signal from the polyphased-filtered signal in frequency domain. The implementation of the pipeline on a GPU/CPU cluster will be discussed together with the computing performance of the prototype., Comment: Proceedings of the 22nd International Conference on Computing in High Energy and Nuclear Physics (CHEP2016), USA

We report here on a novel analysis of the complete set of four Stokes parameters that uniquely determine the linear and/or circular polarization of the radio signal for an extensive air shower. The observed dependency of the circular polarization on azimuth angle and distance to the shower axis is a clear signature of the interfering contributions from two different radiation mechanisms, a main contribution due to a geomagnetically-induced transverse current and a secondary component due to the build-up of excess charge at the shower front. The data, as measured at LOFAR, agree very well with a calculation from first principles. This opens the possibility to use circular polarization as an investigative tool in the analysis of air shower structure, such as for the determination of atmospheric electric fields., Comment: Accepted for publication in Phys. Rev. D

The low flux of the ultra-high energy cosmic rays (UHECR) at the highest energies provides a challenge to answer the long standing question about their origin and nature. A significant increase in the number of detected UHECR is expected to be achieved by employing Earth's moon as detector, and search for short radio pulses that are emitted when a particle interacts in the lunar rock. Observation of these short pulses with current and future radio telescopes also allows to search for the even lower fluxes of neutrinos with energies above $10^{22}$ eV, that are predicted in certain Grand-Unifying-Theories (GUTs), and e.g. models for super-heavy dark matter (SHDM). In this contribution we present the initial design for such a search with the LOFAR radio telescope., Comment: To be published in the Proceedings of the ARENA2016 conference, Groningen, The Netherlands

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], Comment: 22 pages, 17 figures, abstract abridged for Arxiv submission

The characterization of the Galactic foregrounds has been shown to be the main obstacle in the challenging quest to detect primordial B-modes in the polarized microwave sky. We make use of the Planck-HFI 2015 data release at high frequencies to place new constraints on the properties of the polarized thermal dust emission at high Galactic latitudes. Here, we specifically study the spatial variability of the dust polarized spectral energy distribution, and its potential impact on the determination of the tensor-to-scalar ratio. We use the correlation ratio of the $C_\ell^{BB}$ angular power spectra between the 217- and 353-GHz channels as a tracer of these potential variations, computed on different high Galactic latitude regions, ranging from 80% to 20% of the sky. The new insight from Planck data is a departure of the correlation ratio from unity that cannot be attributed to a spurious decorrelation due to the cosmic microwave background, instrumental noise, or instrumental systematics. The effect is marginally detected on each region, but the statistical combination of all the regions gives more than 99% confidence for this variation in polarized dust properties. In addition, we show that the decorrelation increases when there is a decrease in the mean column density of the region of the sky being considered, and we propose a simple power-law empirical model for this dependence, which matches what is seen in the Planck data. We explore the effect that this measured decorrelation has on simulations of the BICEP2-Keck Array/Planck analysis and show that the 2015 constraints from those data still allow a decorrelation between the dust at 150 and 353GHz of the order of the one we measure. Finally we show that either spatial variation of the dust SED or of the dust polarization angle could produce decorrelations between 217- and 353-GHz data similar to those we observe in the data.

Continuum spectra covering centimetre to submillimetre wavelengths are presented for a northern sample of 104 extragalactic radio sources, mainly active galactic nuclei, based on four-epoch Planck data. The nine Planck frequencies, from 30 to 857 GHz, are complemented by a set of simultaneous ground-based radio observations between 1.1 and 37 GHz. The single-survey Planck data confirm that the flattest high-frequency radio spectral indices are close to zero, indicating that the original accelerated electron energy spectrum is much harder than commonly thought, with power-law index around 1.5 instead of the canonical 2.5. The radio spectra peak at high frequencies and exhibit a variety of shapes. For a small set of low-z sources, we find a spectral upturn at high frequencies, indicating the presence of intrinsic cold dust. Variability can generally be approximated by achromatic variations, while sources with clear signatures of evolving shocks appear to be limited to the strongest outbursts., Comment: 39 pages, 113 figures. Accepted by A&A

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., Comment: 26 pages, 25 figures (reduced in quality for arXiv), 1 table. Updated to match version accepted by A&A

Parity violating extensions of the standard electromagnetic theory cause in vacuo rotation of the plane of polarization of propagating photons. This effect, also known as cosmic birefringence, impacts the cosmic microwave background (CMB) anisotropy angular power spectra, producing non-vanishing $T$--$B$ and $E$--$B$ correlations that are otherwise null when parity is a symmetry. Here we present new constraints on an isotropic rotation, parametrized by the angle $\alpha$, derived from Planck 2015 CMB polarization data. To increase the robustness of our analyses, we employ two complementary approaches, in harmonic space and in map space, the latter based on a peak stacking technique. The two approaches provide estimates for $\alpha$ that are in agreement within statistical uncertainties and very stable against several consistency tests. Considering the $T$--$B$ and $E$--$B$ information jointly, we find $\alpha = 0.31^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the harmonic analysis and $\alpha = 0.35^{\circ} \pm 0.05^{\circ} \, ({\rm stat.})\, \pm 0.28^{\circ} \, ({\rm syst.})$ from the stacking approach. These constraints are compatible with no parity violation and are dominated by the systematic uncertainty in the orientation of Planck's polarization-sensitive bolometers., Comment: 15 pages, 8 figures. Accepted for publication in Astronomy and Astrophysics

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., Comment: 19 pages, 18 figures. accepted in A&A

We have carried out a detailed study to understand the observed energy spectrum and composition of cosmic rays with energies up to ~10^18 eV. Our study shows that a single Galactic component with subsequent energy cut-offs in the individual spectra of different elements, optimised to explain the observed spectra below ~10^14 eV and the knee in the all-particle spectrum, cannot explain the observed all-particle spectrum above ~2x10^16 eV. We discuss two approaches for a second component of Galactic cosmic rays -- re-acceleration at a Galactic wind termination shock, and supernova explosions of Wolf-Rayet stars, and show that the latter scenario can explain almost all observed features in the all-particle spectrum and the composition up to ~10^18 eV, when combined with a canonical extra-galactic spectrum expected from strong radio galaxies or a source population with similar cosmological evolution. In this two-component Galactic model, the knee at ~ 3x10^15 eV and the second knee at ~10^17 eV in the all-particle spectrum are due to the cut-offs in the first and second components, respectively. We also discuss several variations of the extra-galactic component, from a minimal contribution to scenarios with a significant component below the ankle (at ~4x10^18 eV), and find that extra-galactic contributions in excess of regular source evolution are neither indicated nor in conflict with the existing data. Our main result is that the second Galactic component predicts a composition of Galactic cosmic rays at and above the second knee that largely consists of helium or a mixture of helium and CNO nuclei, with a weak or essentially vanishing iron fraction, in contrast to most common assumptions. This prediction is in agreement with new measurements from LOFAR and the Pierre Auger Observatory which indicate a strong light component and a rather low iron fraction between ~10^17 and 10^18 eV., Comment: Added Table 4; Published in A&A, 595 (2016) A33 (Highlight paper)

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

We study the statistical properties of interstellar dust polarization at high Galactic latitude, using the Stokes parameter Planck maps at 353 GHz. Our aim is to advance the understanding of the magnetized interstellar medium (ISM), and to provide a model of the polarized dust foreground for cosmic microwave background component-separation procedures. Focusing on the southern Galactic cap, we examine the statistical distributions of the polarization fraction ($p$) and angle ($\psi$) to characterize the ordered and turbulent components of the Galactic magnetic field (GMF) in the solar neighbourhood. We relate patterns at large angular scales in polarization to the orientation of the mean (ordered) GMF towards Galactic coordinates $(l_0,b_0)=(70^\circ \pm 5^\circ,24^\circ \pm 5^\circ)$. The histogram of $p$ shows a wide dispersion up to 25 %. The histogram of $\psi$ has a standard deviation of $12^\circ$ about the regular pattern expected from the ordered GMF. We use these histograms to build a phenomenological model of the turbulent component of the GMF, assuming a uniform effective polarization fraction ($p_0$) of dust emission. To model the Stokes parameters, we approximate the integration along the line of sight (LOS) as a sum over a set of $N$ independent polarization layers, in each of which the turbulent component of the GMF is obtained from Gaussian realizations of a power-law power spectrum. We are able to reproduce the observed $p$ and $\psi$ distributions using: a $p_0$ value of (26 $\pm$ 3)%; a ratio of 0.9 $\pm$ 0.1 between the strengths of the turbulent and mean components of the GMF; and a small value of $N$. We relate the polarization layers to the density structure and to the correlation length of the GMF along the LOS.

We describe a method for spectral cleaning and timing calibration of short voltage time series data from individual radio interferometer receivers. It makes use of the phase differences in Fast Fourier Transform (FFT) spectra across antenna pairs. For strong, localized terrestrial sources these are stable over time, while being approximately uniform-random for a sum over many sources or for noise. Using only milliseconds-long datasets, the method finds the strongest interfering transmitters, a first-order solution for relative timing calibrations, and faulty data channels. No knowledge of gain response or quiescent noise levels of the receivers is required. With relatively small data volumes, this approach is suitable for use in an online system monitoring setup for interferometric arrays. We have applied the method to our cosmic-ray data collection, a collection of measurements of short pulses from extensive air showers, recorded by the LOFAR radio telescope. Per air shower, we have collected 2 ms of raw time series data for each receiver. The spectral cleaning has a calculated optimal sensitivity corresponding to a power signal-to-noise ratio of 0.08 (or -11 dB) in a spectral window of 25 kHz, for 2 ms of data in 48 antennas. This is well sufficient for our application. Timing calibration across individual antenna pairs has been performed at 0.4 ns precision; for calibration of signal clocks across stations of 48 antennas the precision is 0.1 ns. Monitoring differences in timing calibration per antenna pair over the course of the period 2011 to 2015 shows a precision of 0.08 ns, which is useful for monitoring and correcting drifts in signal path synchronizations. A cross-check method for timing calibration is presented, using a pulse transmitter carried by a drone flying over the array. Timing precision is similar, 0.3 ns., Comment: Accepted for publication in Astronomy & Astrophysics; 10 pages, 7 figures

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

Cosmic rays are the highest energy particles found in nature. Measurements of the mass composition of cosmic rays between 10^{17} eV and 10^{18} eV are essential to understand whether this energy range is dominated by Galactic or extragalactic sources. It has also been proposed that the astrophysical neutrino signal comes from accelerators capable of producing cosmic rays of these energies. Cosmic rays initiate cascades of secondary particles (air showers) in the atmosphere and their masses are inferred from measurements of the atmospheric depth of the shower maximum, Xmax, or the composition of shower particles reaching the ground. Current measurements suffer from either low precision, or a low duty cycle and a high energy threshold. Radio detection of cosmic rays is a rapidly developing technique, suitable for determination of Xmax with a duty cycle of in principle nearly 100%. The radiation is generated by the separation of relativistic charged particles in the geomagnetic field and a negative charge excess in the shower front. Here we report radio measurements of Xmax with a mean precision of 16 g/cm^2 between 10^{17}-10^{17.5} eV. Because of the high resolution in $Xmax we can determine the mass spectrum and find a mixed composition, containing a light mass fraction of ~80%. Unless the extragalactic component becomes significant already below 10^{17.5} eV, our measurements indicate an additional Galactic component dominating at this energy range., Comment: 35 pages, 11 figures, updated version: Pierre Auger Observatory data ICRC 2015 added to Fig 2

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