Showing papers by "Anthony Lasenby published in 2019"

Planck 2018 results. V. CMB power spectra and likelihoods

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-multipole 100x143 GHz EE cross-spectrum constrains the reionization optical-depth parameter $\tau$ to better than 15% (in combination with with the other low- and high-$\ell$ likelihoods). We also update the 2015 baseline low-$\ell$ joint TEB likelihood based on the Low Frequency Instrument data, which provides a weaker $\tau$ constraint. At high multipoles, a better model of the temperature-to-polarization leakage and corrections for the effective calibrations of the polarization channels (polarization efficiency or PE) allow us to fully use the polarization spectra, improving the constraints on the $\Lambda$CDM parameters by 20 to 30% compared to TT-only constraints. Tests on the modelling of the polarization demonstrate good consistency, with some residual modelling uncertainties, the accuracy of the PE modelling being the main limitation. Using our various tests, simulations, and comparison between different high-$\ell$ implementations, we estimate the consistency of the results to be better than the 0.5$\sigma$ level. Minor curiosities already present before (differences between $\ell$ 800 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)

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

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

NESTCHECK: diagnostic tests for nested sampling calculations

Abstract: Nested sampling is an increasingly popular technique for Bayesian computation, in particular for multimodal, degenerate problems of moderate to high dimensionality Without appropriate settings, however, nested sampling software may fail to explore such posteriors correctly; for example producing correlated samples or missing important modes This paper introduces new diagnostic tests to assess the reliability both of parameter estimation and evidence calculations using nested sampling software, and demonstrates them empirically We present two new diagnostic plots for nested sampling, and give practical advice for nested sampling software users in astronomy and beyond Our diagnostic tests and diagrams are implemented in nestcheck: a publicly available Python package for analysing nested sampling calculations, which is compatible with output from MultiNest, PolyChord and dyPolyChord

Ghost and tachyon free Poincaré gauge theories: A systematic approach

Abstract: A systematic method is presented for determining the conditions on the parameters in the action of a parity-preserving gauge theory of gravity for it to contain no ghost or tachyon particles. The technique naturally accommodates critical cases in which the parameter values lead to additional gauge invariances. The method is implemented as a computer program, and is used here to investigate the particle content of parity-conserving Poincare gauge theory, which we compare with previous results in the literature. We find 450 critical cases that are free of ghosts and tachyons, and we further identify 10 of these that are also power-counting renormalizable, of which four have only massless tordion propagating particles and the remaining six have only a massive tordion propagating mode.

Constraining the kinetically dominated universe

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 $\mathrm{\ensuremath{\Lambda}}\mathrm{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.

Case for kinetically dominated initial conditions for inflation

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 becomes questionable. We analyze different classes of potentials in phase space and quantify the fraction of the Planck surface that is kinetically dominated. Considering bounded potentials with very small amplitudes as favored by current data and restricting ourselves to the domains of phase space that are of interest to cosmic inflation, we find that initial conditions of the inflaton field should be set in the kinetically dominated regime regardless of any choice of prior.

Bayesian inflationary reconstructions from Planck 2018 data

Abstract: We present three nonparametric 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 $\ensuremath{\ell}\ensuremath{\sim}20$ to $\ensuremath{\ell}\ensuremath{\sim}50$, 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\ensuremath{\lesssim}\ensuremath{\ell}\ensuremath{\lesssim}2000$), but find that parametrizations 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.

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

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

An alternative approach to modelling a cosmic void and its effect on the cosmic microwave background

Abstract: We apply our tetrad-based approach for constructing spherically-symmetric solutions in general relativity to modelling a void, and compare it with the standard Lemaitre-Tolman-Bondi (LTB) formalism. In particular, we construct models for the void observed in the direction of Draco in the WISE-2MASS galaxy survey, and a corresponding cosmic microwave background (CMB) temperature decrement in the Planck data in the same direction. We find that the present-day density and velocity profiles of the void are not well constrained by the existing data, so that void models produced from the two approaches can differ substantially while remaining broadly consistent with the observations. We highlight the importance of considering the velocity as well as the density profile in constraining voids.

QUIJOTE scientific results - III. Microwave spectrum of intensity and polarization in the Taurus Molecular Cloud complex and L1527

Abstract: This work has been partially funded by the Spanish Ministry of Economy and Competitiveness (MINECO) under the projects AYA2007-68058-C03-01, AYA2010-21766-C03-02, AYA2012-39475-C02-01, AYA2014-60438-P: ESP2015- 70646.C2-1-R, AYA2015-64508-P and the Consolider-Ingenio project CSD2010-00064 (EPI: Exploring the Physics of Inflation).

Sunyaev–Zel’dovich profile fitting with joint AMI-Planck analysis

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

Logolinear series expansions with applications to primordial cosmology

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 preinflationary phase of the universe. We develop mathematical techniques for generating these series expansions, and we 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.

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

Abstract: Context. The second legacy catalog of Planck Sunyaev–Zeldovich (SZ) sources, hereafter PSZ2, provides the largest galaxy cluster sample selected by means of the SZ signature of the clusters in a full sky survey. In order to fully characterize this PSZ2 sample for cosmological studies, all the members should be validated and the physical properties of the clusters, including mass and redshift, should be derived. However, at the time of its publication, roughly 21% of the 1653 PSZ2 members had no known counterpart at other wavelengths.Aims. Here, we present the second and last year of observations of our optical follow-up program 128-MULTIPLE-16/15B (hereafter LP15), which has been developed with the aim of validating all the unidentified PSZ2 sources in the northern sky with declinations higher than −15° that have no correspondence in the first Planck catalog PSZ1. The description of the program and the first year of observations have been presented previously.Methods. The LP15 program was awarded 44 observing nights that were spread over two years with the Isaac Newton Telescope (INT), the Telescopio Nazionale Galileo (TNG), and the Gran Telescopio Canarias (GTC), all at Roque de los Muchachos Observatory (La Palma). Following the same method as described previously, we performed deep optical imaging for more than 200 sources with the INT and spectroscopy for almost 100 sources with the TNG and GTC at the end of the LP15 program. We adopted robust confirmation criteria based on velocity dispersion and richness estimates for the final classification of the new galaxy clusters as the optical counterparts of the PSZ2 detections.Results. Here, we present the observations of the second year of LP15, as well as the final results of the program. The full LP15 sample comprises 190 previously unidentified PSZ2 sources. Of these, 106 objects were studied before, while the remaining sample (except for 6 candidates) has been completed in the second year and is discussed here. In addition to the LP15 sample, we here study 42 additional PSZ2 objects that were originally validated as real clusters because they matched a WISE or PSZ1 counterpart, but they had no measured spectroscopic redshift. In total, we confirm the optical counterparts for 81 PSZ2 sources after the full LP15 program, 55 of them with new spectroscopic information. Forty of these 81 clusters are presented in this paper. After the LP15 observational program the purity of the PSZ2 catalog has increased from 76.7% originally to 86.2%. In addition, we study the possible reasons for false detection, and we report a clear correlation between the number of unconfirmed sources and galactic thermal dust emission.

Kinetically dominated curved universes: Logolinear series expansions

Abstract: We develop a method for computing series expansions out of a singularity 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 expanding out of a singularity in a pre-inflationary phase of the universe. We develop mathematical techniques for generating these expansions, and apply them to polynomial and Starobinsky inflationary potentials. This paper is the first in a programme of work on kinetically dominated curved universes, for which such power series are essential. Code for analytic and numerical computation of logolinear series is provided on GitHub.