Showing 1–50 of 400 results for author: Lasenby, A

Balancing ACT: weighing prior dependency and global tensions of DR6 lensing with other datasets

Authors: A. N. Ormondroyd, W. J. Handley, M. P. Hobson, A. N. Lasenby

Abstract: We provide a complementary nested sampling analysis for the Atacama Cosmology Telescope lensing data release 6. This allows the quantification of global consistency statistics between ACT lensing and alternative datasets. In the context of flat $Λ$CDM, we find no inconsistency between ACT, Baryonic Acoustic Oscillations, Planck anisotropies, weak lensing datasets, or NPIPE lensing. As part of our… ▽ More We provide a complementary nested sampling analysis for the Atacama Cosmology Telescope lensing data release 6. This allows the quantification of global consistency statistics between ACT lensing and alternative datasets. In the context of flat $Λ$CDM, we find no inconsistency between ACT, Baryonic Acoustic Oscillations, Planck anisotropies, weak lensing datasets, or NPIPE lensing. As part of our analysis, we also investigate the effect of the prior widths used in the ACT analysis and find that the headline results are quantitatively but not qualitatively affected by the chosen priors. We use both Bayes factors and the suspiciousness statistic to quantify the possibility of tension, and find suspiciousness unsuitable in the case of strong agreement between ACT DR6 and NPIPE. Nested sampling provides a competitive alternative to Metropolis Hastings and we recommend it be used alongside existing analyses. We release the chains and plotting source for the analysis using anesthetic. △ Less

Submitted 3 December, 2023; v1 submitted 12 October, 2023; originally announced October 2023.

Comments: 8 pages, 2 figures, 3 tables. Nested sampling chains and analysis code available at https://zenodo.org/records/10252886. Comments welcome!

QUIJOTE scientific results -- XIII. Intensity and polarization study of supernova remnants in the QUIJOTE-MFI wide survey: CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9

Authors: Carlos Hugo López-Caraballo, Beatriz Ruiz-Granados, Ricardo Genova Santos, Mateo Fernández-Torreiro, Jose Alberto Rubiño-Martin, Mike Peel, Frederick Poidevin, Eduardo Artal, Mark Ashdown, Rita Belen Barreiro, Francisco Javier Casas, Elena de la Hoz, Raul González-González, Federica Guidi, Diego Herranz, Roger Hoyland, Anthony N Lasenby, Enrique Martinez-Gonzalez, Lucio Piccirillo, Rafael Rebolo, Denis Tramonte, Flavien Vansyngel, Patricio Vielva, Robert Watson

Abstract: We use the new QUIJOTE-MFI wide survey (11, 13, 17 and 19 GHz) to produce spectral energy distributions (SEDs), on an angular scale of 1 deg, of the supernova remnants (SNRs) CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9. We provide new measurements of the polarized synchrotron radiation in the microwave range. For each SNR, the intensity and polarization SEDs are obtained and modelled by comb… ▽ More We use the new QUIJOTE-MFI wide survey (11, 13, 17 and 19 GHz) to produce spectral energy distributions (SEDs), on an angular scale of 1 deg, of the supernova remnants (SNRs) CTB 80, Cygnus Loop, HB 21, CTA 1, Tycho and HB 9. We provide new measurements of the polarized synchrotron radiation in the microwave range. For each SNR, the intensity and polarization SEDs are obtained and modelled by combining QUIJOTE-MFI maps with ancillary data. In intensity, we confirm the curved power law spectra of CTB 80 and HB 21 with a break frequency $ν_{\rm b}$ at 2.0$^{+1.2}_{-0.5}$ GHz and 5.0$^{+1.2}_{-1.0}$ GHz respectively; and spectral indices respectively below and above the spectral break of $-0.34\pm0.04$ and $-0.86\pm0.5$ for CTB 80, and $-0.24\pm0.07$ and $-0.60\pm0.05$ for HB 21. In addition, we provide upper limits on the Anomalous Microwave Emission (AME), suggesting that the AME contribution is negligible towards these remnants. From a simultaneous intensity and polarization fit, we recover synchrotron spectral indices as flat as $-0.24$, and the whole sample has a mean and scatter of $-0.44\pm0.12$. The polarization fractions have a mean and scatter of $6.1\pm1.9$\%. When combining our results with the measurements from other QUIJOTE studies of SNRs, we find that radio spectral indices are flatter for mature SNRs, and particularly flatter for CTB 80 ($-0.24^{+0.07}_{-0.06}$) and HB 21 ($-0.34^{+0.04}_{-0.03}$). In addition, the evolution of the spectral indices against the SNRs age is modelled with a power-law function, providing an exponent $-0.07\pm0.03$ and amplitude $-0.49\pm0.02$ (normalised at 10 kyr), which are conservative with respect to previous studies of our Galaxy and the Large Magellanic Cloud. △ Less

Submitted 28 July, 2023; originally announced July 2023.

Comments: 33 pages, 15 figure, 15 tables. Submitted to MNRAS. QUIJOTE data maps available at https://research.iac.es/proyecto/quijote

QUIJOTE Scientific Results -- XVII. Studying the Anomalous Microwave Emission in the Andromeda Galaxy with QUIJOTE-MFI

Authors: M. Fernández-Torreiro, R. T. Génova-Santos, J. A. Rubiño-Martín, C. H. López-Caraballo, M. W. Peel, C. Arce-Tord, R. Rebolo, E. Artal, M. Ashdown, R. B. Barreiro, F. J. Casas, E. de la Hoz, F. Guidi, D. Herranz, R. Hoyland, A. Lasenby, E. Martínez-Gonzalez, L. Piccirillo, F. Poidevin, B. Ruiz-Granados, D. Tramonte, F. Vansyngel, P. Vielva, R. A. Watson

Abstract: The Andromeda Galaxy (M31) is the Local Group galaxy that is most similar to the Milky Way (MW). The similarities between the two galaxies make M31 useful for studying integrated properties common to spiral galaxies. We use the data from the recent QUIJOTE-MFI Wide Survey, together with new raster observations focused on M31, to study its integrated emission. The addition of raster data improves t… ▽ More The Andromeda Galaxy (M31) is the Local Group galaxy that is most similar to the Milky Way (MW). The similarities between the two galaxies make M31 useful for studying integrated properties common to spiral galaxies. We use the data from the recent QUIJOTE-MFI Wide Survey, together with new raster observations focused on M31, to study its integrated emission. The addition of raster data improves the sensitivity of QUIJOTE-MFI maps by almost a factor 3. Our main interest is to confirm if anomalous microwave emission (AME) is present in M31, as previous studies have suggested. To do so, we built the integrated spectral energy distribution of M31 between 0.408 and 3000 GHz. We then performed a component separation analysis taking into account synchrotron, free-free, AME and thermal dust components. AME in M31 is modelled as a log-normal distribution with maximum amplitude, $A_{\rm AME}$, equal to $1.03\pm0.32$ Jy. It peaks at $ν_{\rm AME}=17.2\pm3.2$ GHz with a width of $W_{\rm AME}=0.58\pm0.16$. Both the Akaike and Bayesian Information Criteria find the model without AME to be less than 1 % as probable as the one taking AME into consideration. We find that the AME emissivity per 100 $μ$m intensity in M31 is $ε_{\rm AME}^{\rm 28.4\,GHz}=9.6\pm3.1$ $μ$K/(MJy/sr), similar to that computed for the MW. We also provide the first upper limits for the AME polarization fraction in an extragalactic object. M31 remains the only galaxy where an AME measurement has been made of its integrated spectrum. △ Less

Submitted 13 October, 2023; v1 submitted 15 May, 2023; originally announced May 2023.

Comments: 17 pages, 9 figures. Accepted for publication in MNRAS. QUIJOTE data maps available at https://research.iac.es/proyecto/quijote

QUIJOTE scientific results -- X. Spatial variations of Anomalous Microwave Emission along the Galactic plane

Authors: M. Fernández-Torreiro, J. A. Rubiño-Martín, C. H. López-Caraballo, R. T. Génova-Santos, M. W. Peel, F. Guidi, S. E. Harper, E. Artal, M. Ashdown, R. B. Barreiro, F. J. Casas, E. de la Hoz, D. Herranz, R. Hoyland, A. Lasenby, E. Martínez-Gonzalez, L. Piccirillo, F. Poidevin, R. Rebolo, B. Ruiz-Granados, D. Tramonte, F. Vansyngel, P. Vielva, R. A. Watson

Abstract: Anomalous Microwave Emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquituous in our Galaxy and is observed at a broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17 and 19 GHz to constrain the AME in the Galactic plane ($|b|<10^\circ$) on degree scales. We built the spectral energy… ▽ More Anomalous Microwave Emission (AME) is an important emission component between 10 and 60 GHz that is not yet fully understood. It seems to be ubiquituous in our Galaxy and is observed at a broad range of angular scales. Here we use the new QUIJOTE-MFI wide survey data at 11, 13, 17 and 19 GHz to constrain the AME in the Galactic plane ($|b|<10^\circ$) on degree scales. We built the spectral energy distribution between 0.408 and 3000 GHz for each of the 5309 0.9$^\circ$ pixels in the Galactic plane, and fitted a parametric model by considering five emission components: synchrotron, free-free, AME, thermal dust and CMB anisotropies. We show that not including QUIJOTE-MFI data points leads to the underestimation (up to 50 %) of the AME signal in favour of free-free emission. The parameters describing these components are then intercompared, looking for relations that help to understand AME physical processes. We find median values for the AME width, $W_{\rm AME}$, and for its peak frequency, $ν_{\rm AME}$, respectively of $0.560^{+0.059}_{-0.050}$ and $20.7^{+2.0}_{-1.9}$ GHz, slightly in tension with current theoretical models. We find spatial variations throughout the Galactic plane for $ν_{\rm AME}$, but only with reduced statistical significance. We report correlations of AME parameters with certain ISM properties, such as that between the AME emissivity (which shows variations with the Galactic longitude) and the interstellar radiation field, and that between the AME peak frequency and dust temperature. Finally, we discuss the implications of our results on the possible molecules responsible for AME. △ Less

Submitted 23 August, 2023; v1 submitted 11 May, 2023; originally announced May 2023.

Comments: 32 pages, 31 figures. Accepted for publication in MNRAS. Derived data products available at https://research.iac.es/proyecto/quijote

Does gravitational confinement sustain flat galactic rotation curves without dark matter?

Authors: W. E. V. Barker, M. P. Hobson, A. N. Lasenby

Abstract: The short answer is $\textit{probably no}$. Specifically, this paper considers a recent body of work which suggests that general relativity requires neither the support of dark matter halos, nor unconventional baryonic profiles, nor any infrared modification, to be consistent after all with the anomalously rapid orbits observed in many galactic discs. In particular, the gravitoelectric flux is all… ▽ More The short answer is $\textit{probably no}$. Specifically, this paper considers a recent body of work which suggests that general relativity requires neither the support of dark matter halos, nor unconventional baryonic profiles, nor any infrared modification, to be consistent after all with the anomalously rapid orbits observed in many galactic discs. In particular, the gravitoelectric flux is alleged to collapse nonlinearly into regions of enhanced force, in an analogue of the colour-confining chromoelectric flux tube model which has yet to be captured by conventional post-Newtonian methods. However, we show that the scalar gravity model underpinning this proposal is wholly inconsistent with the nonlinear Einstein equations, which themselves appear to prohibit the linear confinement-type potentials which could indicate a disordered gravitational phase. Our findings challenge the fidelity of the previous Euclidean lattice analyses. We confirm by direct calculation using a number of perturbation schemes and gauges that the next-to-leading order gravitoelectric correction to the rotation curve of a reasonable baryonic profile would be imperceptible. The `gravitoelectric flux collapse' programme was also supported by using intragalactic lensing near a specific galactic baryon profile as a field strength heuristic. We recalculate this lensing effect, and conclude that it has been overstated by three orders of magnitude. As a by-product, our analysis suggests fresh approaches to (i) the fluid ball conjecture and (ii) gravitational energy localisation, both to be pursued in future work. In summary, whilst it may be interesting to consider the possibility of confinement-type effects in gravity, we may at least conclude here that confinement-type effects $\textit{cannot play any significant part}$ in explaining flat or rising galactic rotation curves without dark matter halos. △ Less

Submitted 20 March, 2023; originally announced March 2023.

Comments: 33 pages, 19 figures

Gravitomagnetism and galaxy rotation curves: a cautionary tale

Authors: A. N. Lasenby, M. P. Hobson, W. E. V. Barker

Abstract: We investigate recent claims that gravitomagnetic effects in linearised general relativity can explain flat and rising rotation curves, such as those observed in galaxies, without the need for dark matter. If one models a galaxy as an axisymmetric, stationary, rotating, non-relativistic and pressureless 'dust' of stars in the gravitoelectromagnetic (GEM) formalism, we show that GEM effects on the… ▽ More We investigate recent claims that gravitomagnetic effects in linearised general relativity can explain flat and rising rotation curves, such as those observed in galaxies, without the need for dark matter. If one models a galaxy as an axisymmetric, stationary, rotating, non-relativistic and pressureless 'dust' of stars in the gravitoelectromagnetic (GEM) formalism, we show that GEM effects on the circular velocity $v$ of a star are $O(10^{-6})$ smaller than the standard Newtonian (gravitoelectric) effects. Moreover, we find that gravitomagnetic effects are $O(10^{-6})$ too small to provide the vertical support necessary to maintain the dynamical equilibrium assumed. These issues are obscured if one constructs a single equation for $v$, as considered previously. We nevertheless solve this equation for a galaxy having a Miyamoto--Nagai density profile. We show that for the values of the mass, $M$, and semi-major and semi-minor axes, $a$ and $b$, typical for a dwarf galaxy, the rotation curve depends only very weakly on $M$. Moreover, for aspect ratios $a/b > 2$, the rotation curves are concave over their entire range, which does not match observations in any galaxy. Most importantly, we show that for the poloidal gravitomagnetic flux $ψ$ to provide the necessary vertical support, it must become singular at the origin. This originates from the unwitting, but forbidden, inclusion of free-space solutions of the Poisson-like equation that determines $ψ$, hence ruling out the methodology as a means of explaining flat galaxy rotation curves. We further show that recent deliberate attempts to leverage such free-space solutions against the rotation curve problem yield no deterministic modification outside the thin disk approximation, and that, in any case, the homogeneous contributions to $ψ$ are ruled out by the boundary value problem posed by any physical axisymmetric galaxy. △ Less

Submitted 11 April, 2023; v1 submitted 10 March, 2023; originally announced March 2023.

Comments: Minor textual changes

QUIJOTE scientific results -- IX. Radio sources in the QUIJOTE-MFI wide survey maps

Authors: D. Herranz, M. López-Caniego, C. H. López-Caraballo, R. T. Génova-Santos, Y. C. Perrott, J. A. Rubiño-Martín, R. Rebolo, E. Artal, M. Ashdown, R. B. Barreiro, F. J. Casas, E. de la Hoz, M. Fernández-Torreiro, F. Guidi, R. J. Hoyland, A. N. Lasenby, E. Martínez-González, M. W. Peel, L. Piccirillo, F. Poidevin, B. Ruiz-Granados, D. Tramonte, F. Vansyngel, P. Vielva, R. A. Watson

Abstract: We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extende… ▽ More We present the catalogue of Q-U-I JOint TEnerife (QUIJOTE) Wide Survey radio sources extracted from the maps of the Multi-Frequency Instrument compiled between 2012 and 2018. The catalogue contains 786 sources observed in intensity and polarization, and is divided into two separate sub-catalogues: one containing 47 bright sources previously studied by the \emph{Planck} collaboration and an extended catalogue of 739 sources either selected from the \emph{Planck} Second Catalogue of Compact Sources or found through a blind search carried out with a Mexican Hat 2 wavelet. A significant fraction of the sources in our catalogue (38.7 per cent) are within the $|b| \leq 20^\circ$ region of the Galactic plane. We determine statistical properties for those sources that are likely to be extragalactic. We find that these statistical properties are compatible with currently available models, with a $\sim$1.8 Jy completeness limit at 11 GHz. We provide the polarimetric properties of (38, 33, 31, 23) sources with P detected above the $99.99\%$ significance level at (11, 13, 17, 19) GHz, respectively. Median polarization fractions are in the $2.8$-$4.7$\% range in the 11-19 GHz frequency interval. We do not distinguish between Galactic and extragalactic sources here. The results presented here are consistent with those reported in the literature for flat- and steep-spectrum radio sources. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 21 pages, 13 figures, 7 tables

QUIJOTE scientific results -- VIII. Diffuse polarized foregrounds from component separation with QUIJOTE-MFI

Authors: E. de la Hoz, R. B. Barreiro, P. Vielva, E. Martínez-González, J. A. Rubiño-Martín, B. Casaponsa, F. Guidi, M. Ashdown, R. T. Génova-Santos, E. Artal, F. J. Casas, R. Fernández-Cobos, M. Fernández-Torreiro, D. Herranz, R. J. Hoyland, A. N. Lasenby, M. López-Caniego, C. H. López-Caraballo, M. W. Peel, L. Piccirillo, F. Poidevin, R. Rebolo, B. Ruiz-Granados, D. Tramonte, F. Vansyngel , et al. (1 additional authors not shown)

Abstract: We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foreground… ▽ More We derive linearly polarized astrophysical component maps in the Northern Sky from the QUIJOTE-MFI data at 11 and 13 GHz in combination with the WMAP K and Ka bands (23 and 33 GHz) and all Planck polarized channels (30-353 GHz), using the parametric component separation method B-SeCRET. The addition of QUIJOTE-MFI data significantly improves the parameter estimation of the low-frequency foregrounds, especially the estimation of the synchrotron spectral index, $β_s$. We present the first detailed $β_s$ map of the Northern Celestial Hemisphere at a smoothing scale of $2^{\circ}$. We find statistically significant spatial variability across the sky. We obtain an average value of $-3.08$ and a dispersion of $0.13$, considering only pixels with reliable goodness-of-fit. The power law model of the synchrotron emission provides a good fit to the data outside the Galactic plane but fails to track the complexity within this region. Moreover, when we assume a synchrotron model with uniform curvature, $c_s$, we find a value of $c_s = -0.0797 \pm 0.0012$. However, there is insufficient statistical significance to determine which model is favoured, either the power law or the power law with uniform curvature. Furthermore, we estimate the thermal dust spectral parameters in polarization. Our CMB, synchrotron, and thermal dust maps are highly correlated with the corresponding products of the PR4 Planck release, although some large-scale differences are observed in the synchrotron emission. Finally, we find that the $β_s$ estimation in the high signal-to-noise synchrotron emission areas is prior-independent while, outside these regions, the prior governs the $β_s$ estimation. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 25 pages, 29 figures, Accepted for publication in MNRAS

QUIJOTE Scientific Results -- VII. Galactic AME sources in the QUIJOTE-MFI Northern Hemisphere Wide-Survey

Authors: F. Poidevin, R. T. Génova-Santos, J. A. Rubiño-Martín, C. H. López-Caraballo, R. A. Watson, E. Artal, M. Ashdown, R. B. Barreiro, F. J. Casas, E. de la Hoz, M. Fernández-Torreiro, F. Guidi, D. Herranz, R. J. Hoyland, A. N. Lasenby, E. Martinez-Gonzalez, M. W. Peel, L. Piccirillo, R. Rebolo, B. Ruiz-Granados, D. Tramonte, F. Vansyngel, P. Vielva

Abstract: The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation… ▽ More The QUIJOTE-MFI Northern Hemisphere Wide-Survey has provided maps of the sky above declinations $-30^\circ$ at 11, 13, 17 and 19$\,$GHz. These data are combined with ancillary data to produce Spectral Energy Distributions in intensity in the frequency range 0.4--3\,000$\,$GHz on a sample of 52 candidate compact sources harbouring anomalous microwave emission (AME). We apply a component separation analysis at 1$^\circ$ scale on the full sample from which we identify 44 sources with high AME significance. We explore correlations between different fitted parameters on this last sample. QUIJOTE-MFI data contribute to notably improve the characterisation of the AME spectrum, and its separation from the other components. In particular, ignoring the 10--20\,GHz data produces on average an underestimation of the AME amplitude, and an overestimation of the free-free component. We find an average AME peak frequency of 23.6 $\pm$ 3.6$\,$GHz, about 4$\,$GHz lower than the value reported in previous studies. The strongest correlation is found between the peak flux density of the thermal dust and of the AME component. A mild correlation is found between the AME emissivity ($A_{\rm AME}/τ_{250}$) and the interstellar radiation field. On the other hand no correlation is found between the AME emissivity and the free-free radiation Emission Measure. Our statistical results suggest that the interstellar radiation field could still be the main driver of the intensity of the AME as regards spinning dust excitation mechanisms. On the other hand, it is not clear whether spinning dust would be most likely associated with cold phases of the interstellar medium rather than with hot phases dominated by free-free radiation. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 36 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

QUIJOTE scientific results -- VI. The Haze as seen by QUIJOTE

Authors: F. Guidi, R. T. Génova-Santos, J. A. Rubiño-Martín, M. W. Peel, M. Fernández-Torreiro, C. H. López-Caraballo, R. Vignaga, E. de la Hoz, P. Vielva, R. A. Watson, M. Ashdown, C. Dickinson, E. Artal, R. B. Barreiro, F. J. Casas, D. Herranz, R. J. Hoyland, A. N. Lasenby, E. Martinez-Gonzalez, L. Piccirillo, F. Poidevin, R. Rebolo, B. Ruiz-Granados, D. Tramonte, F. Vansyngel

Abstract: The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $γ$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provid… ▽ More The Haze is an excess of microwave intensity emission surrounding the Galactic centre. It is spatially correlated with the $γ$-ray Fermi bubbles, and with the S-PASS radio polarization plumes, suggesting a possible common provenance. The models proposed to explain the origin of the Haze, including energetic events at the Galactic centre and dark matter decay in the Galactic halo, do not yet provide a clear physical interpretation. In this paper we present a re-analysis of the Haze including new observations from the Multi-Frequency Instrument (MFI) of the Q-U-I JOint TEnerife (QUIJOTE) experiment, at 11 and 13 GHz. We analyze the Haze in intensity and polarization, characterizing its spectrum. We detect an excess of diffuse intensity signal ascribed to the Haze. The spectrum at frequencies 11$\,\leqν\leq\,$70 GHz is a power-law with spectral index $β^{\rm H}=-2.79\pm0.08$, which is flatter than the Galactic synchrotron in the same region ($β^{\rm S}=-2.98\pm0.04$), but steeper than that obtained from previous works ($β^{\rm H}\sim-2.5$ at 23$\,\leq\,ν\leq\,$70 GHz). We also observe an excess of polarized signal in the QUIJOTE-MFI maps in the Haze area. This is a first hint detection of polarized Haze, or a consequence of curvature of the synchrotron spectrum in that area. Finally, we show that the spectrum of polarized structures associated with Galactic centre activity is steep at low frequencies ($β\sim -3.2$ at 2.3 $\leqν\leq$ 23 GHz), and becomes flatter above 11 GHz. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 31 pages. Accepted for publication in MNRAS

QUIJOTE scientific results -- V. The microwave intensity and polarisation spectra of the Galactic regions W49, W51 and IC443

Authors: D. Tramonte, R. T. Génova-Santos, J. A. Rubiño-Martín, P. Vielva, F. Poidevin, C. H. López-Caraballo, M. W. Peel, M. Ashdown, E. Artal, R. B. Barreiro, F. J. Casas, E. de la Hoz, M. Fernández-Torreiro, F. Guidi, D. Herranz, R. J. Hoyland, A. N. Lasenby, E. Martinez-Gonzalez, L. Piccirillo, R. Rebolo, B. Ruiz-Granados, F. Vansyngel, R. A. Watson

Abstract: We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 $μ$K/beam for total intensity and 13-23 $μ$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ra… ▽ More We present new intensity and polarisation maps obtained with the QUIJOTE experiment towards the Galactic regions W49, W51 and IC443, covering the frequency range from 10 to 20 GHz at $\sim$ 1 deg angular resolution, with a sensitivity in the range 35-79 $μ$K/beam for total intensity and 13-23 $μ$K/beam for polarisation. For each region, we combine QUIJOTE maps with ancillary data at frequencies ranging from 0.4 to 3000 GHz, reconstruct the spectral energy distribution and model it with a combination of known foregrounds. We detect anomalous microwave emission (AME) in total intensity towards W49 at 4.7$σ$ and W51 at 4.0$σ$ with peak frequencies $ν_{AME}$ = (20.0 $\pm$ 1.4) GHz and $ν_{AME}$ = (17.7 $\pm$ 3.6) GHz respectively; this is the first detection of AME towards W51. The contamination from ultra-compact HII regions to the residual AME flux density is estimated at 10% in W49 and 5% in W51, and does not rule out the AME detection. The polarised SEDs reveal a synchrotron contribution with spectral indices $α_s$ = -0.67 $\pm$ 0.10 in W49 and $α_s$ = -0.51 $\pm$ 0.07 in W51, ascribed to the diffuse Galactic emission and to the local supernova remnant respectively. Towards IC443 in total intensity we measure a broken power-law synchrotron spectrum with cut-off frequency $ν_{0,s}$ = (114 $\pm$ 73) GHz, in agreement with previous studies; our analysis, however, rules out any AME contribution which had been previously claimed towards IC443. No evidence of polarised AME emission is detected in this study. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 29 pages, 11 figures. Accepted for publication in MNRAS

QUIJOTE scientific results -- IV. A northern sky survey in intensity and polarization at 10-20GHz with the Multi-Frequency Instrument

Authors: J. A. Rubino-Martin, F. Guidi, R. T. Genova-Santos, S. E. Harper, D. Herranz, R. J. Hoyland, A. N. Lasenby, F. Poidevin, R. Rebolo, B. Ruiz-Granados, F. Vansyngel, P. Vielva, R. A. Watson, E. Artal, M. Ashdown, R. B. Barreiro, J. D. Bilbao-Ahedo, F. J. Casas, B. Casaponsa, R. Cepeda-Arroita, E. de la Hoz, C. Dickinson, R. Fernandez-Cobos, M. Fernandez-Torreiro, R. Gonzalez-Gonzalez , et al. (13 additional authors not shown)

Abstract: We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polariza… ▽ More We present QUIJOTE intensity and polarization maps in four frequency bands centred around 11, 13, 17 and 19GHz, and covering approximately 29000 deg$^2$, including most of the Northern sky region. These maps result from 9000 h of observations taken between May 2013 and June 2018 with the first QUIJOTE instrument (MFI), and have angular resolutions of around $1^\circ$, and sensitivities in polarization within the range 35-40 $μ$K per 1-degree beam, being a factor $\sim 2$-$4$ worse in intensity. We discuss the data processing pipeline employed, and the basic characteristics of the maps in terms of real space statistics and angular power spectra. A number of validation tests have been applied to characterise the accuracy of the calibration and the residual level of systematic effects, finding a conservative overall calibration uncertainty of 5%. We also discuss flux densities for four bright celestial sources (Tau A, Cas A, Cyg A and 3C274) which are often used as calibrators at microwave frequencies. The polarization signal in our maps is dominated by synchrotron emission. The distribution of spectral index values between the 11GHz and WMAP 23GHz map peaks at $β=-3.09$ with a standard deviation of 0.14. The measured BB/EE ratio at scales of $\ell=80$ is $0.26\pm 0.07$ for a Galactic cut $|b|>10^\circ$. We find a positive TE correlation for 11GHz at large angular scales ($\ell \lesssim 50$), while the EB and TB signals are consistent with zero in the multipole range $30 \lesssim \ell \lesssim 150$. The maps discussed in this paper are publicly available. △ Less

Submitted 12 January, 2023; originally announced January 2023.

Comments: 60 pages. Accepted for publication in MNRAS. Maps and derived data products available at https://research.iac.es/proyecto/quijote/

Quantum initial conditions for curved inflating universes

Authors: Mary I. Letey, Zakhar Shumaylov, Fruzsina J. Agocs, Will J. Handley, Michael P. Hobson, Anthony N. Lasenby

Abstract: We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique… ▽ More We discuss the challenges of motivating, constructing, and quantizing a canonically normalized inflationary perturbation in spatially curved universes. We show that this has historically proved challenging due to the interaction of nonadiabaticity with spatial curvature. We construct a novel curvature perturbation that is canonically normalized in the sense of its equation of motion and is unique up to a single scalar parameter. With this construction it becomes possible to set initial conditions invariant under canonical transformations, overcoming known ambiguities in the literature. This corrected quantization has potentially observational consequences via modifications to the primordial power spectrum at large angular scales, as well as theoretical implications for quantization procedures in curved cosmologies filled with a scalar field. △ Less

Submitted 8 July, 2024; v1 submitted 30 November, 2022; originally announced November 2022.

Comments: 15 pages, 3 figures; v3: As published in PRD

Journal ref: Phys. Rev. D 109, 123502 (2024)

Nested sampling for physical scientists

Authors: Greg Ashton, Noam Bernstein, Johannes Buchner, Xi Chen, Gábor Csányi, Andrew Fowlie, Farhan Feroz, Matthew Griffiths, Will Handley, Michael Habeck, Edward Higson, Michael Hobson, Anthony Lasenby, David Parkinson, Livia B. Pártay, Matthew Pitkin, Doris Schneider, Joshua S. Speagle, Leah South, John Veitch, Philipp Wacker, David J. Wales, David Yallup

Abstract: We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the applic… ▽ More We review Skilling's nested sampling (NS) algorithm for Bayesian inference and more broadly multi-dimensional integration. After recapitulating the principles of NS, we survey developments in implementing efficient NS algorithms in practice in high-dimensions, including methods for sampling from the so-called constrained prior. We outline the ways in which NS may be applied and describe the application of NS in three scientific fields in which the algorithm has proved to be useful: cosmology, gravitational-wave astronomy, and materials science. We close by making recommendations for best practice when using NS and by summarizing potential limitations and optimizations of NS. △ Less

Submitted 31 May, 2022; originally announced May 2022.

Comments: 20 pages + supplementary information, 5 figures. preprint version; published version at https://www.nature.com/articles/s43586-022-00121-x

Journal ref: Nature Reviews Methods Primers volume 2, Article number: 39 (2022)

Finite inflation in curved space

Authors: Lukas T. Hergt, Fruzsina J. Agocs, Will J. Handley, Michael P. Hobson, Anthony N. Lasenby

Abstract: We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rend… ▽ More We investigate the effects of non-zero spatial curvature on cosmic inflation in the light of cosmic microwave background (CMB) anisotropy measurements from the Planck 2018 legacy release and from the 2015 observing season of BICEP2 and the Keck Array. Even a small percentage of non-zero curvature today would significantly limit the total number of e-folds of the scale factor during inflation, rendering just-enough inflation scenarios with a kinetically dominated or fast-roll stage prior to slow-roll inflation more likely. Finite inflation leads to oscillations and a cutoff towards large scales in the primordial power spectrum and curvature pushes them into the CMB observable window. Using nested sampling, we carry out Bayesian parameter estimations and model comparisons taking into account constraints from reheating and horizon considerations. We confirm the preference of CMB data for closed universes with Bayesian odds of over $100:1$ and with a posterior on the curvature density parameter of $Ω_{K,0}=-0.051\pm0.017$ for a curvature extension of LCDM and $Ω_{K,0}=-0.031\pm0.014$ for Starobinsky inflation. Model comparisons of various inflation models give similar results as for flat universes with the Starobinsky model outperforming most other models. △ Less

Submitted 15 May, 2022; originally announced May 2022.

Comments: 46 pages, 32 figures, 3 tables, submitted to PRD, from chapter 5 of L. T. Hergt's doctoral thesis at https://www.repository.cam.ac.uk/handle/1810/323129, data available at https://zenodo.org/record/6547872

Perturbations and the Future Conformal Boundary

Authors: A. N. Lasenby, W. J. Handley, D. J. Bartlett, C. S. Negreanu

Abstract: The concordance model of cosmology predicts a universe which finishes in a finite amount of conformal time at a future conformal boundary. We show that for particular cases we study, the background variables and perturbations may be analytically continued beyond this boundary and that the "end of the universe" is not necessarily the end of their physical development. Remarkably, these theoretical… ▽ More The concordance model of cosmology predicts a universe which finishes in a finite amount of conformal time at a future conformal boundary. We show that for particular cases we study, the background variables and perturbations may be analytically continued beyond this boundary and that the "end of the universe" is not necessarily the end of their physical development. Remarkably, these theoretical considerations of the end of the universe might have observable consequences today: perturbation modes consistent with these boundary conditions have a quantised power spectrum which may be relevant to features seen in the large scale cosmic microwave background. Mathematically these cosmological models may either be interpreted as a palindromic universe mirrored in time, a reflecting boundary condition, or a double cover, but are identical with respect to their observational predictions and stand in contrast to the predictions of conformal cyclic cosmologies. △ Less

Submitted 20 April, 2022; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: 29 pages, 32 figures

Journal ref: Phys. Rev. D 105, 083514 - Published 19 April 2022

Improved cosmological fits with quantized primordial power spectra

Authors: D. J. Bartlett, W. J. Handley, A. N. Lasenby

Abstract: We observationally examine cosmological models based on primordial power spectra with quantized wavevectors. Introducing a linearly quantized power spectrum with $k_0=3.225\times10^{-4}\mathrm{Mpc}^{-1}$ and spacing $Δk = 2.257 \times 10^{-4} \mathrm{Mpc}^{-1}$ provides a better fit to the Planck 2018 observations than the concordance baseline, with $Δχ^2 = -8.55$. Extending the results of Lasenby… ▽ More We observationally examine cosmological models based on primordial power spectra with quantized wavevectors. Introducing a linearly quantized power spectrum with $k_0=3.225\times10^{-4}\mathrm{Mpc}^{-1}$ and spacing $Δk = 2.257 \times 10^{-4} \mathrm{Mpc}^{-1}$ provides a better fit to the Planck 2018 observations than the concordance baseline, with $Δχ^2 = -8.55$. Extending the results of Lasenby et al [1], we show that the requirement for perturbations to remain finite beyond the future conformal boundary in a universe containing dark matter and a cosmological constant results in a linearly quantized primordial power spectrum. It is found that the infrared cutoffs for this future conformal boundary quantized cosmology do not provide cosmic microwave background power spectra compatible with observations, but future theories may predict more observationally consistent quantized spectra. △ Less

Submitted 20 April, 2022; v1 submitted 5 April, 2021; originally announced April 2021.

Comments: 14 pages, 8 figures, prepared for submission to PRD

Journal ref: Phys. Rev. D 105, 083515 - Published 19 April 2022

Conformal gravity does not predict flat galaxy rotation curves

Authors: Michael Hobson, Anthony Lasenby

Abstract: We reconsider the widely held view that the Mannheim--Kazanas (MK) vacuum solution for a static, spherically-symmetric system in conformal gravity (CG) predicts flat rotation curves, such as those observed in galaxies, without the need for dark matter. This prediction assumes that test particles have fixed rest mass and follow timelike geodesics in the MK metric in the vacuum region exterior to a… ▽ More We reconsider the widely held view that the Mannheim--Kazanas (MK) vacuum solution for a static, spherically-symmetric system in conformal gravity (CG) predicts flat rotation curves, such as those observed in galaxies, without the need for dark matter. This prediction assumes that test particles have fixed rest mass and follow timelike geodesics in the MK metric in the vacuum region exterior to a spherically-symmetric representation of the galactic mass distribution. Such geodesics are not conformally invariant, however, which leads to an apparent discrepancy with the analogous calculation performed in the conformally-equivalent Schwarzschild-de-Sitter (SdS) metric, where the latter does not predict flat rotation curves. This difference arises since the mass of particles in CG must instead be generated dynamically through interaction with a scalar field. The energy-momentum of this required scalar field means that, in a general conformal frame from the equivalence class of CG solutions outside a static, spherically-symmetric matter distribution, the spacetime is not given by the MK vacuum solution. A unique frame does exist, however, for which the metric retains the MK form, since the scalar field energy-momentum vanishes despite the field being non-zero and radially dependent. Nonetheless, we show that in both this MK frame and the Einstein frame, in which the scalar field is constant, massive particles follow timelike geodesics of the SdS metric, thereby resolving the apparent frame dependence of physical predictions and unambiguously yielding rotation curves with no flat region. We also comment on how our analysis resolves the long-standing uncertainty regarding gravitational lensing in the MK metric. (Abridged) △ Less

Submitted 10 September, 2021; v1 submitted 24 March, 2021; originally announced March 2021.

Comments: 13 pages, no figures

Journal ref: Phys. Rev. D 104, 064014 (2021)

Bayesian evidence for the tensor-to-scalar ratio $r$ and neutrino masses $m_ν$: Effects of uniform vs logarithmic priors

Authors: Lukas T. Hergt, Will J. Handley, Michael P. Hobson, Anthony N. Lasenby

Abstract: We review the effect that the choice of a uniform or logarithmic prior has on the Bayesian evidence and hence on Bayesian model comparisons when data provide only a one-sided bound on a parameter. We investigate two particular examples: the tensor-to-scalar ratio $r$ of primordial perturbations and the mass of individual neutrinos $m_ν$, using the cosmic microwave background temperature and polari… ▽ More We review the effect that the choice of a uniform or logarithmic prior has on the Bayesian evidence and hence on Bayesian model comparisons when data provide only a one-sided bound on a parameter. We investigate two particular examples: the tensor-to-scalar ratio $r$ of primordial perturbations and the mass of individual neutrinos $m_ν$, using the cosmic microwave background temperature and polarisation data from Planck 2018 and the NuFIT 5.0 data from neutrino oscillation experiments. We argue that the Kullback-Leibler divergence, also called the relative entropy, mathematically quantifies the Occam penalty. We further show how the Bayesian evidence stays invariant upon changing the lower prior bound of an upper constrained parameter. While a uniform prior on the tensor-to-scalar ratio disfavours the $r$-extension compared to the base LCDM model with odds of about 1:20, switching to a logarithmic prior renders both models essentially equally likely. LCDM with a single massive neutrino is favoured over an extension with variable neutrino masses with odds of 20:1 in case of a uniform prior on the lightest neutrino mass, which decreases to roughly 2:1 for a logarithmic prior. For both prior options we get only a very slight preference for the normal over the inverted neutrino hierarchy with Bayesian odds of about 3:2 at most. △ Less

Submitted 23 February, 2021; originally announced February 2021.

Comments: 16 pages, 9 figures, 3 tables. Submitted to PRD. Data available at https://doi.org/10.5281/zenodo.4556360

Journal ref: Phys. Rev. D 103, 123511 (2021)

28 -- 40 GHz variability and polarimetry of bright compact sources in the QUIJOTE cosmological fields

Authors: Yvette C. Perrott, Marcos López-Caniego, Ricardo T. Génova-Santos, Jose Alberto Rubiño-Martín, Mark Ashdown, Diego Herranz, Anne Lähteenmäki, Anthony N. Lasenby, Carlos H. López-Caraballo, Frédérick Poidevin, Merja Tornikoski

Abstract: We observed 51 sources in the Q-U-I JOint TEnerife (QUIJOTE) cosmological fields which were brighter than 1 Jy at 30 GHz in the Planck Point Source Catalogue (version 1), with the Very Large Array at 28 -- 40 GHz, in order to characterise their high-radio-frequency variability and polarization properties. We find a roughly log-normal distribution of polarization fractions with a median of 2%, in a… ▽ More We observed 51 sources in the Q-U-I JOint TEnerife (QUIJOTE) cosmological fields which were brighter than 1 Jy at 30 GHz in the Planck Point Source Catalogue (version 1), with the Very Large Array at 28 -- 40 GHz, in order to characterise their high-radio-frequency variability and polarization properties. We find a roughly log-normal distribution of polarization fractions with a median of 2%, in agreement with previous studies, and a median rotation measure (RM) of $\approx$ 1110 rad m$^{-2}$ with one outlier up to $\approx$ 64000 rad m$^{-2}$ which is among the highest RMs measured in quasar cores. We find hints of a correlation between the total intensity flux density and median polarization fraction. We find 59% of sources are variable in total intensity, and 100% in polarization at $3σ$ level, with no apparent correlation between total intensity variability and polarization variability. This indicates that it will be difficult to model these sources without simultaneous polarimetric monitoring observations and they will need to be masked for cosmological analysis. △ Less

Submitted 8 February, 2021; originally announced February 2021.

Comments: 17 pages, 14 figures, accepted to MNRAS

Nonlinear Hamiltonian analysis of new quadratic torsion theories Part I. Cases with curvature-free constraints

Authors: W. E. V. Barker, A. N. Lasenby, M. P. Hobson, W. J. Handley

Abstract: It was recently found that, when linearised in the absence of matter, 58 cases of the general gravitational theory with quadratic curvature and torsion are (i) free from ghosts and tachyons and (ii) power-counting renormalisable. We inspect the nonlinear Hamiltonian structure of the eight cases whose primary constraints do not depend on the curvature tensor. We confirm the particle spectra and uni… ▽ More It was recently found that, when linearised in the absence of matter, 58 cases of the general gravitational theory with quadratic curvature and torsion are (i) free from ghosts and tachyons and (ii) power-counting renormalisable. We inspect the nonlinear Hamiltonian structure of the eight cases whose primary constraints do not depend on the curvature tensor. We confirm the particle spectra and unitarity of all these theories in the linear regime. We uncover qualitative dynamical changes in the nonlinear regimes of all eight cases, suggesting at least a broken gauge symmetry, and possibly the activation of negative kinetic energy spin-parity sectors and acausal behaviour. Two of the cases propagate a pair of massless modes at the linear level, and were interesting as candidate theories of gravity. However, we identify these modes with vector excitations, rather than the tensor polarisations of the graviton. Moreover, we show that these theories do not support a viable cosmological background. △ Less

Submitted 7 January, 2021; originally announced January 2021.

Comments: 20 pages, 0 figures

Journal ref: Phys. Rev. D 104, 084036 (2021)

Planck intermediate results. LV. Reliability and thermal properties of high-frequency sources in the Second Planck Catalogue of Compact Sources

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, C. Burigana, E. Calabrese, P. Carvalho, H. C. Chiang, B. P. Crill, F. Cuttaia, A. de Rosa, G. de Zotti , et al. (95 additional authors not shown)

Abstract: 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 sour… ▽ More 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. △ Less

Submitted 14 September, 2020; originally announced September 2020.

Comments: 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)

Journal ref: A&A 644, A99 (2020)

Dense output for highly oscillatory numerical solutions

Authors: F. J. Agocs, M. P. Hobson, W. J. Handley, A. N. Lasenby

Abstract: We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution… ▽ More We present a method to construct a continuous extension (otherwise known as dense output) for a numerical routine in the special case of the numerical solution being a scalar-valued function exhibiting rapid oscillations. Such cases call for numerical routines that make use of the known global behaviour of the solution, one example being methods using asymptotic expansions to forecast the solution at each step of the independent variable. An example is oscode, numerical routine which uses the Wentzel-Kramers-Brillouin (WKB) approximation when the solution oscillates rapidly and otherwise behaves as a Runge-Kutta (RK) solver. Polynomial interpolation is not suitable for producing the solution at an arbitrary point mid-step, since efficient numerical methods based on the WKB approximation will step through multiple oscillations in a single step. Instead we construct the continuous solution by extending the numerical quadrature used in computing a WKB approximation of the solution with no additional evaluations of the differential equation or terms within, and provide an error estimate on this dense output. Finally, we draw attention to previous work on the continuous extension of Runge-Kutta formulae, and construct an extension to a RK method based on Gauss--Lobatto quadrature nodes, thus describing how to generate dense output from each of the methods underlying oscode. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Comments: 10 pages, 5 figures, 4 tables. Submitted to PRResearch

Planck intermediate results. LVII. Joint Planck LFI and HFI data processing

Authors: Planck Collaboration, Y. Akrami, K. J. Andersen, M. Ashdown, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, C. Burigana, R. C. Butler, E. Calabrese, B. Casaponsa, H. C. Chiang, L. P. L. Colombo, C. Combet, B. P. Crill, F. Cuttaia , et al. (114 additional authors not shown)

Abstract: 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… ▽ More 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$)$μ$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 $τ= 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. △ Less

Submitted 9 July, 2020; originally announced July 2020.

Comments: 97 pages, 93 figures and 16 tables, abstract abridged for arXiv submission, accepted for publication in A&A

Journal ref: A&A 643, A42 (2020)

Mapping Poincaré gauge cosmology to Horndeski theory for emergent dark energy

Authors: W. E. V. Barker, A. N. Lasenby, M. P. Hobson, W. J. Handley

Abstract: The ten-parameter, quadratic Poincaré gauge theory of gravity is a plausible alternative to general relativity. We show that the rich background cosmology of the gauge theory is described by a non-canonical bi-scalar-tensor theory in the Jordan frame: the `metrical analogue'. This provides a unified framework for future investigation by the broader community. For many parameter choices, the non-ca… ▽ More The ten-parameter, quadratic Poincaré gauge theory of gravity is a plausible alternative to general relativity. We show that the rich background cosmology of the gauge theory is described by a non-canonical bi-scalar-tensor theory in the Jordan frame: the `metrical analogue'. This provides a unified framework for future investigation by the broader community. For many parameter choices, the non-canonical term reduces to a Cuscuton field of the form $\smash{\sqrt{|X^{φφ}|}}$. The Einstein-Cartan-Kibble-Sciama theory maps to a pure quadratic Cuscuton, whereas the teleparallel theory maps to the Einstein-Hilbert Lagrangian. We apply the metrical analogue to novel unitary and power-counting-renormalisable cases of Poincaré gauge theory. These theories support the concordance $Λ$CDM background cosmology up to an optional, effective dark radiation component, we explain this behaviour in terms of a stalled Cuscuton. We also obtain two dark energy solutions from one of these cases: accelerated expansion from a negative bare cosmological constant whose magnitude is screened, and emergent dark energy to replace vanishing bare cosmological constant in $Λ$CDM. △ Less

Submitted 7 September, 2020; v1 submitted 5 June, 2020; originally announced June 2020.

Comments: 11 pages, 2 figures. Results previously presented at DAMTP on 14/02/20. Added substantial introduction to Poincaré gauge theory and its perturbative renormalisability. Added various references. Incorporated supplemental material into main text (previously hosted on GitHub). Titular change. Results unaffected. Accepted for publication by Phys. Rev. D

Journal ref: Phys. Rev. D 102, 084002 (2020)

Optical validation and characterisation of Planck PSZ1 sources at the Canary Islands observatories. II. Second year of ITP13 observations

Authors: R. Barrena, A. Ferragamo, J. A. Rubiño-Martín, A. Streblyanska, A. Aguado-Barahona, D. Tramonte, R. T. Génova-Santos, A. Hempel, H. Lietzen, N. Aghanim, H. Böhringer, G. Chon, H. Dahle, M. Douspis, A. N. Lasenby, P. Mazzotta, J. B. Melin, E. Pointecouteau, G. W. Pratt, M. Rossetti

Abstract: We report new galaxy clusters previously unknown included in the PSZ1 catalogue. The results here presented have been achieved during the second year of a 2-year observational programme, the ITP13, developed at the Roque de los Muchachos Observatory (La Palma, Spain). We characterise 75 SZ sources with low SZ significance, SZ S/N$<5.32$ by performing deep optical imaging and spectroscopy in order… ▽ More We report new galaxy clusters previously unknown included in the PSZ1 catalogue. The results here presented have been achieved during the second year of a 2-year observational programme, the ITP13, developed at the Roque de los Muchachos Observatory (La Palma, Spain). We characterise 75 SZ sources with low SZ significance, SZ S/N$<5.32$ by performing deep optical imaging and spectroscopy in order to associate actual galaxy clusters to the SZ Planck source. We adopt robust criteria, based on the 2D-spatial distribution, richness and velocity dispersions to confirm actual optical counterparts up to $z<0.85$. At the end of the ITP13 observational programme, we study 256 SZ sources with $Dec \geq -15^{\circ}$ (212 of them completely unknown), finding optical counterparts for 152 SZ sources. The ITP13 validation programme has allowed us to update the PSZ1 purity, which is now more refined, increasing from 72\% to 83\% in the low SZ S/N regime. Our results are consistent with the predicted purity curve for the full PSZ1 catalogue and with the expected fraction of false detections caused by the non-Gaussian noise of foreground signals. Indeed, we find a strong correlation between the number of unconfirmed sources and the thermal emission of diffuse galactic dust at 857 GHz, thus increasing the fraction of false Planck SZ detections at low galactic latitudes. △ Less

Submitted 16 April, 2020; originally announced April 2020.

Comments: Accepted for publications in A&A. 15 pages including 14 figures and 1 long table

Journal ref: A&A 638, A146 (2020)

Addressing $H_0$ tension with emergent dark radiation in unitary gravity

Authors: W. E. V. Barker, A. N. Lasenby, M. P. Hobson, W. J. Handley

Abstract: We propose a one-parameter extension to $Λ$CDM, expected to strongly affect cosmological tensions. An effective dark radiation component in the early universe redshifts away as hot dark matter, then quintessence, tracking the dominant equation-of-state parameter and leaving a falsifiable torsion field in the current epoch. This picture results from a new Poincaré gauge theory (PGT), one of the mos… ▽ More We propose a one-parameter extension to $Λ$CDM, expected to strongly affect cosmological tensions. An effective dark radiation component in the early universe redshifts away as hot dark matter, then quintessence, tracking the dominant equation-of-state parameter and leaving a falsifiable torsion field in the current epoch. This picture results from a new Poincaré gauge theory (PGT), one of the most promising among the latest batch of 58 PGTs found to be both power-counting renormalisable and free from ghosts and tachyons. We systematically categorise the cosmologies of 33 of these PGTs, as special cases of the most general parity-preserving, Ostrogradsky-stable PGT with a purely Yang-Mills action. The theory we consider contains two propagating massless gravitons, which may be $J^P=2^+$ (long-range gravitation and gravitational waves). A conspiracy among the coupling constants eliminates the spatial curvature $k\in\{\pm 1,0\}$ from the field equations. We show that this `$k$-screening' is not restricted to conformal gravity theories. The flat Friedmann equations are then emergent, with potentially tension-resolving freedom at the early scale-invariant epoch that reliably gives way to an attractor-like state of modern $Λ$CDM evolution. We compare with related theories and promising special cases, such as $k$-screened theories with negative-definite effective $k$, and more traditional theories with effective $Λ$ and a $J^P=0^-$ massive graviton (dark matter candidate). As a bonus, we analyse similarly constrained actions in the new extended Weyl gauge theory (eWGT). We show that in cosmology, PGT and eWGT span exactly the same classical phenomenology up to a linear map between their coupling constants, hinting at a deeper relationship between the two. △ Less

Submitted 13 July, 2020; v1 submitted 5 March, 2020; originally announced March 2020.

Comments: 32 pages, 5 figures, 2 tables. Titular change. Updated ref. [85] and ref. [96]. Version to appear in Phys. Rev. D

Journal ref: Phys. Rev. D 102, 024048 (2020)

Quantum initial conditions for inflation and canonical invariance

Authors: F. J. Agocs, L. T. Hergt, W. J. Handley, A. N. Lasenby, M. P. Hobson

Abstract: We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since… ▽ More We investigate the transformation of initial conditions for primordial curvature perturbations under two types of transformations of the associated action: simultaneous redefinition of time and the field to be quantised, and the addition of surface terms. The latter encompasses all canonical transformations, whilst the time- and field-redefinition is a distinct, non-canonical transformation since the initial and destination systems use different times. Actions related to each other via such transformations yield identical equations of motion and preserve the commutator structure. They further preserve the time-evolution of expectation values of quantum operators unless the vacuum state also changes under the transformation. These properties suggest that it is of interest to investigate vacuum prescriptions that also remain unchanged under canonical transformations. We find that initial conditions derived via minimising the vacuum expectation value of the Hamiltonian and those obtained using the Danielsson vacuum prescription are not invariant under these transformations, whereas those obtained by minimising the local energy density are. We derive the range of physically distinct initial conditions obtainable by Hamiltonian diagonalisation, and illustrate their effect on the scalar primordial power spectrum and the Cosmic Microwave Background under the just enough inflation model. We also generalise the analogy between the dynamics of a quantum scalar field on a curved, time-dependent spacetime and the gauge-invariant curvature perturbation. We argue that the invariance of the vacuum prescription obtained by minimising the renormalised stress--energy tensor should make it the preferred procedure for setting initial conditions for primordial perturbations. All other procedures reviewed in this work yield ambiguous initial conditions, which is problematic both in theory and practice. △ Less

Submitted 27 May, 2020; v1 submitted 14 February, 2020; originally announced February 2020.

Comments: Revisions after review added. 17 pages, 8 figures. Accepted to PRD

Journal ref: Phys. Rev. D 102, 023507 (2020)

Blueshifted absorption lines from X-ray reflection in IRAS 13224-3809

Authors: A. C. Fabian, C. S. Reynolds, J. Jiang, C. Pinto, L. C. Gallo, M. L. Parker, A. N. Lasenby, W. N. Alston, D. J. K Buisson, E. M. Cackett, B. De Marco, J. Garcia, E. Kara, P. Kosec, M. J. Middleton, J. M. Miller, G. Miniutti, D. J. Walton, D. R. Wilkins, A. J. Young

Abstract: We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorpti… ▽ More We explore a disc origin for the highly-blueshifted, variable absorption lines seen in the X-ray spectrum of the Narrow Line Seyfert 1 galaxy IRAS13224-3809. The blueshift corresponds to a velocity of about 0.25c. Such features in other Active Galactic Nuclei are often interpreted as UltraFast Outflows (UFOs). The velocity is of course present in the orbital motions of the inner disk. The absorption lines in IRAS13224-3809 are best seen when the flux is low and the reflection component of the disk is strong relative to the power-law continuum. The spectra are consistent with a model in which the reflection component passes through a thin, highly-ionized absorbing layer at the surface of the inner disc, the blue-shifted side of which dominates the flux due to relativistic aberration (the disc inclination is about 70 deg). No fast outflow need occur beyond the disc. △ Less

Submitted 15 February, 2020; originally announced February 2020.

Comments: 6 pages, 5 figures, accepted for publication in MNRAS

Detection of Spectral Variations of Anomalous Microwave Emission with QUIJOTE and C-BASS

Authors: R. Cepeda-Arroita, S. Harper, C. Dickinson, J. A. Rubiño-Martín, R. T. Génova-Santos, Angela C. Taylor, T. J. Pearson, M. Ashdown, A. Barr, R. B. Barreiro, B. Casaponsa, F. J. Casas, H. C. Chiang, R. Fernandez-Cobos, R. D. P. Grumitt, F. Guidi, H. M. Heilgendorff, D. Herranz, L. R. P. Jew, J. L. Jonas, Michael E. Jones, A. Lasenby, J. Leech, J. P. Leahy, E. Martínez-González , et al. (10 additional authors not shown)

Abstract: Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $λ$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and… ▽ More Anomalous Microwave Emission (AME) is a significant component of Galactic diffuse emission in the frequency range $10$-$60\,$GHz and a new window into the properties of sub-nanometre-sized grains in the interstellar medium. We investigate the morphology of AME in the $\approx10^{\circ}$ diameter $λ$ Orionis ring by combining intensity data from the QUIJOTE experiment at $11$, $13$, $17$ and $19\,$GHz and the C-Band All Sky Survey (C-BASS) at $4.76\,$GHz, together with 19 ancillary datasets between $1.42$ and $3000\,$GHz. Maps of physical parameters at $1^{\circ}$ resolution are produced through Markov Chain Monte Carlo (MCMC) fits of spectral energy distributions (SEDs), approximating the AME component with a log-normal distribution. AME is detected in excess of $20\,σ$ at degree-scales around the entirety of the ring along photodissociation regions (PDRs), with three primary bright regions containing dark clouds. A radial decrease is observed in the AME peak frequency from $\approx35\,$GHz near the free-free region to $\approx21\,$GHz in the outer regions of the ring, which is the first detection of AME spectral variations across a single region. A strong correlation between AME peak frequency, emission measure and dust temperature is an indication for the dependence of the AME peak frequency on the local radiation field. The AME amplitude normalised by the optical depth is also strongly correlated with the radiation field, giving an overall picture consistent with spinning dust where the local radiation field plays a key role. △ Less

Submitted 25 February, 2021; v1 submitted 20 January, 2020; originally announced January 2020.

Comments: 19 pages, 7 figures, accepted for publication by MNRAS

Astrometric Effects of Gravitational Wave Backgrounds with non-Luminal Propagation Speeds

Authors: Deyan P. Mihaylov, Christopher J. Moore, Jonathan Gair, Anthony Lasenby, Gerard Gilmore

Abstract: A passing gravitational wave causes a deflection in the apparent astrometric positions of distant stars. The effect of the speed of the gravitational wave on this astrometric shift is discussed. A stochastic background of gravitational waves would result in a pattern of astrometric deflections which are correlated on large angular scales. These correlations are quantified and investigated for back… ▽ More A passing gravitational wave causes a deflection in the apparent astrometric positions of distant stars. The effect of the speed of the gravitational wave on this astrometric shift is discussed. A stochastic background of gravitational waves would result in a pattern of astrometric deflections which are correlated on large angular scales. These correlations are quantified and investigated for backgrounds of gravitational waves with sub- and super-luminal group velocities. The statistical properties of the correlations are depicted in two equivalent and related ways: as correlation curves and as angular power spectra. Sub-(super-)luminal gravitational wave backgrounds have the effect of enhancing (suppressing) the power in low-order angular modes. Analytical representations of the redshift-redshift and redshift-astrometry correlations are also derived. The potential for using this effect for constraining the speed of gravity is discussed. △ Less

Submitted 23 November, 2019; originally announced November 2019.

Comments: 14 pages, 16 figures, accepted for publication in Physical Review D

Journal ref: Phys. Rev. D 101, 024038 (2020)

Astrophysical Black Holes

Authors: Andrew C. Fabian, Anthony N. Lasenby

Abstract: Black holes are a common feature of the Universe. They are observed as stellar mass black holes spread throughout galaxies and as supermassive objects in their centres. Observations of stars orbiting close to the centre of our Galaxy provide detailed clear evidence for the presence of a 4 million Solar mass black hole. Gas accreting onto distant supermassive black holes produces the most luminous… ▽ More Black holes are a common feature of the Universe. They are observed as stellar mass black holes spread throughout galaxies and as supermassive objects in their centres. Observations of stars orbiting close to the centre of our Galaxy provide detailed clear evidence for the presence of a 4 million Solar mass black hole. Gas accreting onto distant supermassive black holes produces the most luminous persistent sources of radiation observed, outshining galaxies as quasars. The energy generated by such displays may even profoundly affect the fate of a galaxy. We briefly review the history of black holes and relativistic astrophysics before exploring the observational evidence for black holes and reviewing current observations including black hole mass and spin. In parallel we outline the general relativistic derivation of the physical properties of black holes relevant to observation. Finally we speculate on future observations and touch on black hole thermodynamics and the extraction of energy from rotating black holes. △ Less

Submitted 11 November, 2019; originally announced November 2019.

Comments: Published in 2015 as Chapter 1 of 'General Relativity: The most beautiful of theories: Applications and trends after 100 years', Edited by Carlo Rovelli (pub. De Gruyter)

Bayesian inflationary reconstructions from Planck 2018 data

Authors: Will J. Handley, Anthony N. Lasenby, Hiranya V. Peiris, Michael P. Hobson

Abstract: We present three non-parametric Bayesian primordial reconstructions using Planck 2018 polarization data: linear spline primordial power spectrum reconstructions, cubic spline inflationary potential reconstructions and sharp-featured primordial power spectrum reconstructions. All three methods conditionally show hints of an oscillatory feature in the primordial power spectrum in the multipole range… ▽ More We present three non-parametric Bayesian primordial reconstructions using Planck 2018 polarization data: linear spline primordial power spectrum reconstructions, cubic spline inflationary potential reconstructions and sharp-featured primordial power spectrum reconstructions. All three methods conditionally show hints of an oscillatory feature in the primordial power spectrum in the multipole range $\ell\sim20$ to $\ell\sim50$, which is to some extent preserved upon marginalization. We find no evidence for deviations from a pure power law across a broad observable window ($50\lesssim\ell\lesssim2000$), but find that parameterizations are preferred which are able to account for lack of resolution at large angular scales due to cosmic variance, and at small angular scales due to Planck instrument noise. Furthermore, the late-time cosmological parameters are unperturbed by these extensions to the primordial power spectrum. This work is intended to provide a background and give more details of the Bayesian primordial reconstruction work found in the Planck 2018 papers. △ Less

Submitted 11 November, 2019; v1 submitted 2 August, 2019; originally announced August 2019.

Comments: 25 Pages, 15 Figures

Journal ref: Phys. Rev. D 100, 103511 (2019)

Planck 2018 results. V. CMB power spectra and likelihoods

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (143 additional authors not shown)

Abstract: 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… ▽ More 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 $τ$ 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 $τ$ 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 $Λ$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$σ$ 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) △ Less

Submitted 15 September, 2020; v1 submitted 30 July, 2019; originally announced July 2019.

Comments: Revised to match version published in Astronomy & Astrophysics

Journal ref: A&A 641, A5 (2020)

Planck 2018 results. VII. Isotropy and Statistics of the CMB

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso, B. Casaponsa, H. C. Chiang , et al. (125 additional authors not shown)

Abstract: 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 $Λ$CDM cosmological model, yet also confirm the presence of several so-called "anomalies" on la… ▽ More 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 $Λ$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 $Λ$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. △ Less

Submitted 14 September, 2020; v1 submitted 6 June, 2019; originally announced June 2019.

Comments: Paper VII of the Planck 2018 release, revised to closely match version published in Astronomy and Astrophysics

Journal ref: A&A 641, A7 (2020)

An efficient method for solving highly oscillatory ordinary differential equations with applications to physical systems

Authors: F. J. Agocs, W. J. Handley, A. N. Lasenby, M. P. Hobson

Abstract: We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequenc… ▽ More We present a novel numerical routine (oscode) with a C++ and Python interface for the efficient solution of one-dimensional, second-order, ordinary differential equations with rapidly oscillating solutions. The method is based on a Runge-Kutta-like stepping procedure that makes use of the Wentzel-Kramers-Brillouin (WKB) approximation to skip regions of integration where the characteristic frequency varies slowly. In regions where this is not the case, the method is able to switch to a made-to-measure Runge-Kutta integrator that minimises the total number of function evaluations. We demonstrate the effectiveness of the method with example solutions of the Airy equation and an equation exhibiting a burst of oscillations, discussing the error properties of the method in detail. We then show the method applied to physical systems. First, the one-dimensional, time-independent Schrödinger equation is solved as part of a shooting method to search for the energy eigenvalues for a potential with quartic anharmonicity. Then, the method is used to solve the Mukhanov-Sasaki equation describing the evolution of cosmological perturbations, and the primordial power spectrum of the perturbations is computed in different cosmological scenarios. We compare the performance of our solver in calculating a primordial power spectrum of scalar perturbations to that of BINGO, an efficient code specifically designed for such applications. △ Less

Submitted 13 December, 2019; v1 submitted 30 May, 2019; originally announced June 2019.

Comments: Physical Review Research, accepted. 23 pages, 15 figures. The associated code is available online at https://github.com/fruzsinaagocs/oscode

Journal ref: Phys. Rev. Research 2, 013030 (2020)

Planck 2018 results. IX. Constraints on primordial non-Gaussianity

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso, B. Casaponsa, A. Challinor , et al. (135 additional authors not shown)

Abstract: 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 polariz… ▽ More 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). △ Less

Submitted 14 May, 2019; originally announced May 2019.

Comments: 50 pages, 20 figures

Sunyaev--Zel'dovich profile fitting with joint AMI-Planck analysis

Authors: Yvette C. Perrott, Kamran Javid, Pedro Carvalho, Patrick J. Elwood, Michael P. Hobson, Anthony N. Lasenby, Malak Olamaie, Richard D. E. Saunders

Abstract: We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the an… ▽ More We develop a Bayesian method of analysing Sunyaev-Zel'dovich measurements of galaxy clusters obtained from the Arcminute Microkelvin Imager (AMI) radio interferometer system and from the Planck satellite, using a joint likelihood function for the data from both instruments. Our method is applicable to any combination of Planck data with interferometric data from one or more arrays. We apply the analysis to simulated clusters and find that when the cluster pressure profile is known a-priori, the joint dataset provides precise and accurate constraints on the cluster parameters, removing the need for external information to reduce the parameter degeneracy. When the pressure profile deviates from that assumed for the fit, the constraints become biased. Allowing the pressure profile shape parameters to vary in the analysis allows an unbiased recovery of the integrated cluster signal and produces constraints on some shape parameters, depending on the angular size of the cluster. When applied to real data from Planck-detected cluster PSZ2 G063.80+11.42, our method resolves the discrepancy between the AMI and Planck $Y$-estimates and usefully constrains the gas pressure profile shape parameters at intermediate and large radii. △ Less

Submitted 9 February, 2020; v1 submitted 28 January, 2019; originally announced January 2019.

Comments: 13 pages, 10 figures. Submitted to MNRAS. Post-referee version

Journal ref: MNRAS 486 (2019) 2116-2128

Logolinear series expansions with applications to primordial cosmology

Authors: Will Handley, Anthony Lasenby, Mike Hobson

Abstract: We develop a method for computing series expansions for solutions to ordinary differential equations when the asymptotic form contains both linear and logarithmic terms. Such situations are common in primordial cosmology when considering series expansions out of a singularity in the equations arising from a pre-inflationary phase of the universe. We develop mathematical techniques for generating t… ▽ More We develop a method for computing series expansions for solutions to ordinary differential equations when the asymptotic form contains both linear and logarithmic terms. Such situations are common in primordial cosmology when considering series expansions out of a singularity in the equations arising from a pre-inflationary phase of the universe. We develop mathematical techniques for generating these series expansions, and apply them to polynomial and Starobinsky inflationary potentials with kinetic initial conditions. Code for analytic and numerical computation of logolinear series is provided on GitHub. △ Less

Submitted 14 June, 2019; v1 submitted 21 January, 2019; originally announced January 2019.

Comments: 12 pages, 4 figures

Journal ref: Phys. Rev. D 99, 123512 (2019)

QUIJOTE Scientific results. III. Microwave spectrum of intensity and polarization in the Taurus molecular cloud Complex and L1527

Authors: F. Poidevin, J. A. Rubiño-Martín, C. Dickinson, R. Génova-Santos, S. Harper, R. Rebolo, B. Casaponsa, A. Peláez-Santos, R. Vignaga, F. Guidi, B. Ruiz-Granados, D. Tramonte, F. Vansyngel, M. Ashdown, D. Herranz, R. Hoyland, A. Lasenby, E. Martínez-González, L. Piccirillo, R. A. Watson

Abstract: We present new intensity and polarization observations of the Taurus molecular cloud (TMC) region in the frequency range 10-20 GHz with the Multi-Frequency Instrument (MFI) mounted on the first telescope of the QUIJOTE experiment. From the combination of the QUIJOTE data with the WMAP 9-yr data release, the Planck second data release, the DIRBE maps and ancillary data, we detect an anomalous micro… ▽ More We present new intensity and polarization observations of the Taurus molecular cloud (TMC) region in the frequency range 10-20 GHz with the Multi-Frequency Instrument (MFI) mounted on the first telescope of the QUIJOTE experiment. From the combination of the QUIJOTE data with the WMAP 9-yr data release, the Planck second data release, the DIRBE maps and ancillary data, we detect an anomalous microwave emission (AME) component with flux density $S_{\rm AME, peak} = 43.0 \pm 7.9\,$Jy in the Taurus Molecular Cloud (TMC) and $S_{\rm AME, peak} = 10.7 \pm 2.7\,$Jy in the dark cloud nebula L1527, which is part of the TMC. In the TMC the diffuse AME emission peaks around a frequency of 19 GHz, compared with an emission peak about a frequency of 25 GHz in L1527. In the TMC, the best constraint on the level of AME polarisation is obtained at the Planck channel of 28.4 GHz, with an upper limit $π_{\rm AME}<$4.2$\,\%$ (95$\,\%$ C. L.), which reduces to $π_{\rm AME} <$3.8$\,\%$ (95$\,\%$ C.L.) if the intensity of all the free-free, synchrotron and thermal dust components are negligible at this frequency. The same analysis in L1527 leads to $π_{\rm AME}<$5.3$\%$ (95$\,\%$C.L.), or $π_{\rm AME}<$4.5$\,\%$ (95$\%$C.L.) under the same assumption. We find that in the TMC and L1527 on average about $80\%$ of the HII gas should be mixed with thermal dust. Our analysis shows how the QUIJOTE-MFI 10-20 GHz data provides key information to properly separate the synchrotron, free-free and AME components. △ Less

Submitted 14 January, 2019; v1 submitted 19 December, 2018; originally announced December 2018.

Comments: Accepted for publication in the MNRAS, 27 pages, 19 Figures, 12 Tables

Constraining the kinetically dominated Universe

Authors: L. T. Hergt, W. J. Handley, M. P. Hobson, A. N. Lasenby

Abstract: We present cosmological constraints from Planck 2015 data for a universe that is kinetically dominated at very early times. We perform a Markov chain Monte Carlo analysis to estimate parameters and use nested sampling to determine the evidence for a model comparison of the single-field quadratic and Starobinsky inflationary models with the standard $Λ$CDM cosmology. In particular we investigate ho… ▽ More We present cosmological constraints from Planck 2015 data for a universe that is kinetically dominated at very early times. We perform a Markov chain Monte Carlo analysis to estimate parameters and use nested sampling to determine the evidence for a model comparison of the single-field quadratic and Starobinsky inflationary models with the standard $Λ$CDM cosmology. In particular we investigate how different amounts of inflation before and after horizon exit affect the primordial power spectrum and subsequently the power spectrum of the cosmic microwave background. We find that the model using kinetically dominated initial conditions for inflation performs similarly well in terms of Bayesian evidence as a model directly starting out in the slow-roll phase, despite having an additional parameter. The data show a slight preference for a cutoff at large scales in the primordial and temperature power spectra. △ Less

Submitted 20 September, 2018; originally announced September 2018.

Comments: Submitted to Physical Review D, 15 pages, 12 figures, 1 table

Journal ref: Phys. Rev. D 100, 023501 (2019)

A case for kinetically dominated initial conditions for inflation

Authors: L. T. Hergt, W. J. Handley, M. P. Hobson, A. N. Lasenby

Abstract: We make a case for setting initial conditions for inflation at the Planck epoch in the kinetically dominated regime. For inflationary potentials with a plateau or a hill, i.e. potentials that are bounded from above within a certain region of interest, we cannot claim complete ignorance of the energy distribution between kinetic and potential energy, and equipartition of energy at the Planck epoch… ▽ More We make a case for setting initial conditions for inflation at the Planck epoch in the kinetically dominated regime. For inflationary potentials with a plateau or a hill, i.e. potentials that are bounded from above within a certain region of interest, we cannot claim complete ignorance of the energy distribution between kinetic and potential energy, and equipartition of energy at the Planck epoch becomes questionable. We analyse different classes of potentials in phase-space and quantify the fraction of the Planck surface that is kinetically dominated. For the small amplitudes of the potentials as suggested by current data, the Planck surface lies in the region of kinetic dominance for almost all values of interest of the inflaton field. △ Less

Submitted 19 September, 2018; originally announced September 2018.

Comments: Submitted to PRL, 6 pages, 4 figures

Journal ref: Phys. Rev. D 100, 023502 (2019)

Bayesian sparse reconstruction: a brute-force approach to astronomical imaging and machine learning

Authors: Edward Higson, Will Handley, Michael Hobson, Anthony Lasenby

Abstract: We present a principled Bayesian framework for signal reconstruction, in which the signal is modelled by basis functions whose number (and form, if required) is determined by the data themselves. This approach is based on a Bayesian interpretation of conventional sparse reconstruction and regularisation techniques, in which sparsity is imposed through priors via Bayesian model selection. We demons… ▽ More We present a principled Bayesian framework for signal reconstruction, in which the signal is modelled by basis functions whose number (and form, if required) is determined by the data themselves. This approach is based on a Bayesian interpretation of conventional sparse reconstruction and regularisation techniques, in which sparsity is imposed through priors via Bayesian model selection. We demonstrate our method for noisy 1- and 2-dimensional signals, including astronomical images. Furthermore, by using a product-space approach, the number and type of basis functions can be treated as integer parameters and their posterior distributions sampled directly. We show that order-of-magnitude increases in computational efficiency are possible from this technique compared to calculating the Bayesian evidences separately, and that further computational gains are possible using it in combination with dynamic nested sampling. Our approach can also be readily applied to neural networks, where it allows the network architecture to be determined by the data in a principled Bayesian manner by treating the number of nodes and hidden layers as parameters. △ Less

Submitted 25 November, 2018; v1 submitted 12 September, 2018; originally announced September 2018.

Comments: 18 pages + appendix, 19 figures, minor updates to text and layout

Journal ref: Mon. Notices Royal Astron. Soc. 483, 4 (2019) p4828-4846

Planck 2018 results. XII. Galactic astrophysics using polarized dust emission

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. I. R. Alves, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, A. Bracco, M. Bucher, C. Burigana, E. Calabrese , et al. (138 additional authors not shown)

Abstract: We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, 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… ▽ More We present 353 GHz full-sky maps of the polarization fraction $p$, angle $ψ$, 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$, $ψ$, 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. △ Less

Submitted 7 March, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication in Astronomy & Astrophysics

Journal ref: A&A 641, A12 (2020)

Planck 2018 results. X. Constraints on inflation

Authors: Planck Collaboration, Y. Akrami, F. Arroja, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (151 additional authors not shown)

Abstract: 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… ▽ More 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" (φ) < 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. △ Less

Submitted 2 August, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: References added and minor improvements. BICEP2/Keck Array BK15 is used in the place of BICEP2/Keck Array BK14

Journal ref: A&A 641, A10 (2020)

Planck 2018 results. VIII. Gravitational lensing

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron , et al. (133 additional authors not shown)

Abstract: 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\,σ$ to $9\,σ$. Combined with temperature, lensing is detected at $40\,σ$. We present an extensive set of tests of the robustness of the lensing-potential… ▽ More 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\,σ$ to $9\,σ$. Combined with temperature, lensing is detected at $40\,σ$. 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{ΛCDM}$ model. Combined with baryon density and other weak priors, the lensing analysis alone constrains $σ_8 Ω_{\rm m}^{0.25}=0.589\pm 0.020$ ($1\,σ$ errors). Also combining with baryon acoustic oscillation (BAO) data, we find tight individual parameter constraints, $σ_8=0.811\pm0.019$, $H_0=67.9_{-1.3}^{+1.2}\,\text{km}\,\text{s}^{-1}\,\rm{Mpc}^{-1}$, and $Ω_{\rm m}=0.303^{+0.016}_{-0.018}$. Combining with $\textit{Planck}$ CMB power spectrum data, we measure $σ_8$ to better than $1\,\%$ precision, finding $σ_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). △ Less

Submitted 29 July, 2019; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 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

Journal ref: A&A 641, A8 (2020)

Planck 2018 results. VI. Cosmological parameters

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, R. Battye, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. J. Bock, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese , et al. (157 additional authors not shown)

Abstract: 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 $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis g… ▽ More 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 $Λ$CDM cosmology having a power-law spectrum of adiabatic scalar perturbations (denoted "base $Λ$CDM" in this paper), from polarization, temperature, and lensing, separately and in combination. A combined analysis gives dark matter density $Ω_c h^2 = 0.120\pm 0.001$, baryon density $Ω_b h^2 = 0.0224\pm 0.0001$, scalar spectral index $n_s = 0.965\pm 0.004$, and optical depth $τ= 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θ_*=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-$Λ$CDM cosmology, the inferred late-Universe parameters are: Hubble constant $H_0 = (67.4\pm 0.5)$km/s/Mpc; matter density parameter $Ω_m = 0.315\pm 0.007$; and matter fluctuation amplitude $σ_8 = 0.811\pm 0.006$. We find no compelling evidence for extensions to the base-$Λ$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_ν< 0.12$eV. The CMB spectra continue to prefer higher lensing amplitudes than predicted in base -$Λ$CDM at over $2\,σ$, which pulls some parameters that affect the lensing amplitude away from the base-$Λ$CDM model; however, this is not supported by the lensing reconstruction or (in models that also change the background geometry) BAO data. (Abridged) △ Less

Submitted 9 August, 2021; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 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

Journal ref: A&A 641, A6 (2020)

Planck 2018 results. IV. Diffuse component separation

Authors: Planck Collaboration, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron, B. Casaponsa, A. Challinor, L. P. L. Colombo , et al. (128 additional authors not shown)

Abstract: 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 (Comm… ▽ More 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. △ Less

Submitted 26 September, 2020; v1 submitted 17 July, 2018; originally announced July 2018.

Comments: 74 pages, A&A, 641, A4

Planck 2018 results. III. High Frequency Instrument data processing and frequency maps

Authors: Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, E. Calabrese, J. -F. Cardoso, J. Carron, A. Challinor , et al. (130 additional authors not shown)

Abstract: 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… ▽ More 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μ$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. △ Less

Submitted 17 July, 2018; originally announced July 2018.

Comments: Accepted for publication on A&A (AA/2018/32909)

Planck 2018 results. II. Low Frequency Instrument data processing

Authors: Planck Collaboration, Y. Akrami, F. Argüeso, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, S. Basak, K. Benabed, J. -P. Bernard, M. Bersanelli, P. Bielewicz, L. Bonavera, J. R. Bond, J. Borrill, F. R. Bouchet, F. Boulanger, M. Bucher, C. Burigana, R. C. Butler, E. Calabrese, J. -F. Cardoso , et al. (126 additional authors not shown)

Abstract: 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.… ▽ More 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. △ Less

Submitted 11 September, 2018; v1 submitted 17 July, 2018; originally announced July 2018.

Journal ref: A&A 641, A2 (2020)