Issue 
A&A
Volume 594, October 2016
Planck 2015 results



Article Number  A18  
Number of page(s)  21  
Section  Cosmology (including clusters of galaxies)  
DOI  https://doi.org/10.1051/00046361/201525829  
Published online  20 September 2016 
Planck 2015 results
XVIII. Background geometry and topology of the Universe
^{1} APC, AstroParticule et Cosmologie, Université Paris Diderot, CNRS/IN2P3, CEA/lrfu, Observatoire de Paris, Sorbonne Paris Cité, 10 rue Alice Domon et Léonie Duquet, 75205 Paris Cedex 13, France
^{2} Aalto University Metsähovi Radio Observatory and Dept of Radio Science and Engineering, PO Box 13000 00076 Aalto, Finland
^{3} African Institute for Mathematical Sciences, 6−8 Melrose Road, Muizenberg, Cape Town, South Africa
^{4} Agenzia Spaziale Italiana Science Data Center, via del Politecnico snc, 00133 Roma, Italy
^{5} Aix Marseille Université, CNRS, LAM (Laboratoire d’Astrophysique de Marseille) UMR 7326, 13388 Marseille, France
^{6} Astrophysics Group, Cavendish Laboratory, University of Cambridge, J J Thomson Avenue, Cambridge CB3 0HE, UK
^{7} Astrophysics & Cosmology Research Unit, School of Mathematics, Statistics & Computer Science, University of KwaZuluNatal, Westville Campus, Private Bag X54001, 4000 Durban, South Africa
^{8} CGEE, SCS Qd 9, Lote C, Torre C, 4° andar, Ed. Parque Cidade Corporate, CEP 70308200, Brasília, DF, Brazil
^{9} CITA, University of Toronto, 60 St. George St., Toronto, ON M5S 3H8, Canada
^{10} CNRS, IRAP, 9 Av. colonel Roche, BP 44346, 31028 Toulouse Cedex 4, France
^{11} CRANN, Trinity College, 2 Dublin, Ireland
^{12} California Institute of Technology, Pasadena, California, USA
^{13} Centre for Theoretical Cosmology, DAMTP, University of Cambridge, Wilberforce Road, Cambridge CB3 0WA, UK
^{14} Centro de Estudios de Física delCosmos de Aragón (CEFCA), Plaza San Juan, 1, planta 2, 44001 Teruel, Spain
^{15} Computational Cosmology Center, Lawrence Berkeley National Laboratory, Berkeley, CA94720 California, USA
^{16} Consejo Superior de Investigaciones Científicas (CSIC), 28006 Madrid, Spain
^{17} DSM/Irfu/SPP, CEASaclay, 91191 GifsurYvette Cedex, France
^{18} DTU Space, National Space Institute, Technical University of Denmark, Elektrovej 327, 2800 Kgs. Lyngby, Denmark
^{19} Département de Physique Théorique, Université de Genève, 24 Quai E. Ansermet, 1211 Genève 4, Switzerland
^{20} Departamento de Astrofísica, Universidad de La Laguna (ULL), 38206 La Laguna, Tenerife, Spain
^{21} Departamento de Física, Universidad de Oviedo, Avda. Calvo Sotelo s/n, 33003 Oviedo, Spain
^{22} Department of Astronomy and Astrophysics, University of Toronto, 50 Saint George Street, Toronto, ONMSS Ontario, Canada
^{23} Department of Astrophysics/IMAPP, Radboud University Nijmegen, PO Box 9010, 6500 GL Nijmegen, The Netherlands
^{24} Department of Physics & Astronomy, University of British Columbia, 6224 Agricultural Road, Vancouver, British Columbia, Canada
^{25} Department of Physics and Astronomy, Dana and David Dornsife College of Letter, Arts and Sciences, University of Southern California, Los Angeles, CA 90089, USA
^{26} Department of Physics and Astronomy, University College London, London WC1E 6BT, UK
^{27} Department of Physics, Florida State University, Keen Physics Building, 77 Chieftan Way, Tallahassee, Florida, USA
^{28} Department of Physics, Gustaf Hällströmin katu 2a, University of Helsinki, 00100 Helsinki, Finland
^{29} Department of Physics, Princeton University, Princeton, New Jersey NJ 08544, USA
^{30} Department of Physics, University of Alberta, 1132289 Avenue, Edmonton, Alberta, T6G 2G7, Canada
^{31} Department of Physics, University of California, Santa Barbara, California CA 93106, USA
^{32} Department of Physics, University of Illinois at UrbanaChampaign, 1110 West Green Street, Urbana, Illinois, USA
^{33} Dipartimento di Fisica e Astronomia G. Galilei, Università degli Studi di Padova, via Marzolo 8, 35131 Padova, Italy
^{34} Dipartimento di Fisica e Scienze della Terra, Università di Ferrara, via Saragat 1, 44122 Ferrara, Italy
^{35} Dipartimento di Fisica, Università La Sapienza, P.le A. Moro 2, Roma, Italy
^{36} Dipartimento di Fisica, Università degli Studi di Milano, via Celoria, 16 Milano, Italy
^{37} Dipartimento di Fisica, Università degli Studi di Trieste, via A. Valerio 2, Trieste, Italy
^{38} Dipartimento di Matematica, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, Roma, Italy
^{39} Discovery Center, Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
^{40} Discovery Center, Niels Bohr Institute, Copenhagen University, Blegdamsvej 17, Copenhagen, Denmark
^{41} European Space Agency, ESAC, Planck Science Office, Camino bajo del Castillo, s/n, Urbanización Villafranca del Castillo, Villanueva de la Cañada, 28692 Madrid, Spain
^{42} European Space Agency, ESTEC, Keplerlaan 1, 2201 AZ Noordwijk, The Netherlands
^{43} Gran Sasso Science Institute, INFN, viale F. Crispi 7, 67100 L’ Aquila, Italy
^{44} HGSFP and University of Heidelberg, Theoretical Physics Department, Philosophenweg 16, 69120 Heidelberg, Germany
^{45} Helsinki Institute of Physics, Gustaf Hällströmin katu 2, University of Helsinki, Helsinki, Finland
^{46} INAF−Osservatorio Astronomico di Padova, Vicolo dell’Osservatorio 5, Padova, Italy
^{47} INAF−Osservatorio Astronomico di Roma, via di Frascati 33, Monte Porzio Catone, Italy
^{48} INAF−Osservatorio Astronomico di Trieste, via G.B. Tiepolo 11, Trieste, Italy
^{49} INAF/IASF Bologna, via Gobetti 101, Bologna, Italy
^{50} INAF/IASF Milano, via E. Bassini 15, Milano, Italy
^{51} INFN, Sezione di Bologna, viale Berti Pichat 6/2, 40127 Bologna, Italy
^{52} INFN, Sezione di Ferrara, via Saragat 1, 44122 Ferrara, Italy
^{53} INFN, Sezione di Roma 1, Università di Roma Sapienza, Piazzale Aldo Moro 2, 00185, Roma, Italy
^{54} INFN, Sezione di Roma 2, Università di Roma Tor Vergata, via della Ricerca Scientifica 1, Roma, Italy
^{55} INFN/National Institute for Nuclear Physics, via Valerio 2, 34127 Trieste, Italy
^{56} IPAG: Institut de Planétologie et d’Astrophysique de Grenoble, Université Grenoble Alpes, IPAG; CNRS, IPAG, 38000 Grenoble, France
^{57} IUCAA, Post Bag 4, Ganeshkhind, Pune University Campus, 411 007 Pune, India
^{58} Imperial College London, Astrophysics group, Blackett Laboratory, Prince Consort Road, London, SW7 2AZ, UK
^{59} Infrared Processing and Analysis Center, California Institute of Technology, Pasadena, CA 91125, USA
^{60} Institut Néel, CNRS, Université Joseph Fourier Grenoble I, 25 rue des Martyrs, Grenoble, France
^{61} Institut Universitaire de France, 103 Bd SaintMichel, 75005 Paris, France
^{62} Institut d’Astrophysique Spatiale, CNRS, Univ. ParisSud, Université ParisSaclay, Bât. 121, 91405 Orsay Cedex, France
^{63} Institut d’Astrophysique de Paris, CNRS (UMR 7095), 98 bis Boulevard Arago, 75014 Paris, France
^{64} Institut für Theoretische Teilchenphysik und Kosmologie, RWTH Aachen University, 52056 Aachen, Germany
^{65} Institute for Space Sciences, BucharestMagurale, 077125 Bucharest, Romania
^{66} Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK
^{67} Institute of Theoretical Astrophysics, University of Oslo, 1072 Blindern, Oslo, Norway
^{68} Instituto de Astrofísica de Canarias, C/Vía Láctea s/n, La Laguna, 38205 Tenerife, Spain
^{69} Instituto de Física de Cantabria (CSICUniversidad de Cantabria), Avda. de los Castros s/n, 39005 Santander, Spain
^{70} Istituto Nazionale di Fisica Nucleare, Sezione di Padova, via Marzolo 8, 35131 Padova, Italy
^{71} Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Drive, Pasadena, California, USA
^{72} Jodrell Bank Centre for Astrophysics, Alan Turing Building, School of Physics and Astronomy, The University of Manchester, Oxford Road, Manchester, M13 9PL, UK
^{73} Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA
^{74} Kavli Institute for Cosmology Cambridge, Madingley Road, Cambridge, CB3 0HA, UK
^{75} Kazan Federal University, 18 Kremlyovskaya St., 420008 Kazan, Russia
^{76} LAL, Université ParisSud, CNRS/IN2P3, Orsay, France
^{77} LERMA, CNRS, Observatoire de Paris, 61 Avenue de l’Observatoire, Paris, France
^{78} Laboratoire AIM, IRFU/Service d’Astrophysique− CEA/DSM − CNRS − Université Paris Diderot, Bât. 709, CEASaclay, 1191 GifsurYvette Cedex, France
^{79} Laboratoire Traitement et Communication de l’Information, CNRS (UMR 5141) and Télécom ParisTech, 46 rue Barrault, 75634 Paris Cedex 13, France
^{80} Laboratoire de Physique Subatomique et Cosmologie, Université GrenobleAlpes, CNRS/IN2P3, 53 rue des Martyrs, 38026 Grenoble Cedex, France
^{81} Laboratoire de Physique Théorique, Université ParisSud 11 & CNRS, Bâtiment 210, 91405 Orsay, France
^{82} Lawrence Berkeley National Laboratory, Berkeley, CA 94720 California, USA
^{83} Lebedev Physical Institute of the Russian Academy of Sciences, Astro Space Centre, 84/32 Profsoyuznaya st., 117997 Moscow, Russia
^{84} MaxPlanckInstitut für Astrophysik, KarlSchwarzschildStr. 1, 85741 Garching, Germany
^{85} McGill Physics, Ernest Rutherford Physics Building, McGill University, 3600 rue University, Montréal, QC, H3A 2T8, Canada
^{86} Mullard Space Science Laboratory, University College London, Surrey RH5 6NT, UK
^{87} National University of Ireland, Department of Experimental Physics, Maynooth, Co. Kildare, Ireland
^{88} Nicolaus Copernicus Astronomical Center, Bartycka 18, 00716 Warsaw, Poland
^{89} Niels Bohr Institute, Blegdamsvej 17, Copenhagen, Denmark
^{90} Niels Bohr Institute, Copenhagen University, Blegdamsvej 17, Copenhagen, Denmark
^{91} Nordita (Nordic Institute for Theoretical Physics), Roslagstullsbacken 23, 106 91 Stockholm, Sweden
^{92} Optical Science Laboratory, University College London, Gower Street, WC 1E6 BT London, UK
^{93} SISSA, Astrophysics Sector, via Bonomea 265, 34136 Trieste, Italy
^{94} SMARTEST Research Centre, Università degli Studi eCampus, via Isimbardi 10, 22060 Novedrate (CO), Italy
^{95} School of Physics and Astronomy, Cardiff University, Queens Buildings, The Parade, Cardiff, CF24 3AA, UK
^{96} School of Physics and Astronomy, University of Nottingham, Nottingham NG7 2RD, UK
^{97} Sorbonne UniversitéUPMC, UMR7095, Institut d’Astrophysique de Paris, 98 bis Boulevard Arago, 75014 Paris, France
^{98} Space Sciences Laboratory, University of California, Berkeley, CA 94720 California, USA
^{99} Special Astrophysical Observatory, Russian Academy of Sciences, Nizhnij Arkhyz, Zelenchukskiy region, 369167 KarachaiCherkessian Republic, Russia
^{100} Stanford University, Dept of Physics, Varian Physics Bldg, 382 via Pueblo Mall, Stanford, California, USA
^{101} SubDepartment of Astrophysics, University of Oxford, Keble Road, Oxford OX1 3RH, UK
^{102} The Oskar Klein Centre for Cosmoparticle Physics, Department of Physics, Stockholm University, AlbaNova, 106 91 Stockholm, Sweden
^{103} Theory Division, PHTH, CERN, 1211 Geneva 23, Switzerland
^{104} UPMC Univ. Paris 06, UMR7095, 98 bis Boulevard Arago, 75014 Paris, France
^{105} Université de Toulouse, UPSOMP, IRAP, 31028 Toulouse Cedex 4, France
^{106} University of Granada, Departamento de Física Teórica y del Cosmos, Facultad de Ciencias, 18010 Granada, Spain
^{107} University of Granada, Instituto Carlos I de Física Teórica y Computacional, Granada, Spain
^{108} Warsaw University Observatory, Aleje Ujazdowskie 4, 00478 Warszawa, Poland
^{⋆}
Corresponding author: A. H. Jaffe, email: a.jaffe@imperial.ac.uk
Received: 6 February 2015
Accepted: 10 April 2016
Maps of cosmic microwave background (CMB) temperature and polarization from the 2015 release of Planck data provide the highestquality fullsky view of the surface of last scattering available to date. This enables us to detect possible departures from a globally isotropic cosmology. We present the first searches using CMB polarization for correlations induced by a possible nontrivial topology with a fundamental domain that intersects, or nearly intersects, the lastscattering surface (at comoving distance χ_{rec}), both via a direct scan for matched circular patterns at the intersections and by an optimal likelihood calculation for specific topologies. We specialize to flat spaces with cubic toroidal (T3) and slab (T1) topologies, finding that explicit searches for the latter are sensitive to other topologies with antipodal symmetry. These searches yield no detection of a compact topology with a scale below the diameter of the lastscattering surface. The limits on the radius ℛ_{i} of the largest sphere inscribed in the fundamental domain (at loglikelihood ratio Δlnℒ > −5 relative to a simplyconnected flat Planck bestfit model) are: ℛ_{i} > 0.97 χ_{rec} for the T3 cubic torus; and ℛ_{i} > 0.56 χ_{rec} for the T1 slab. The limit for the T3 cubic torus from the matchedcircles search is numerically equivalent, ℛ_{i} > 0.97 χ_{rec} at 99% confidence level from polarization data alone. We also perform a Bayesian search for an anisotropic global Bianchi VII_{h} geometry. In the nonphysical setting, where the Bianchi cosmology is decoupled from the standard cosmology, Planck temperature data favour the inclusion of a Bianchi component with a Bayes factor of at least 2.3 units of logevidence. However, the cosmological parameters that generate this pattern are in strong disagreement with those found from CMB anisotropy data alone. Fitting the induced polarization pattern for this model to the Planck data requires an amplitude of −0.10 ± 0.04 compared to the value of + 1 if the model were to be correct. In the physically motivated setting, where the Bianchi parameters are coupled and fitted simultaneously with the standard cosmological parameters, we find no evidence for a Bianchi VII_{h} cosmology and constrain the vorticity of such models to (ω/H)_{0} < 7.6 × 10^{10} (95% CL).
Key words: cosmic background radiation / cosmology: observations / cosmological parameters / gravitation / methods: data analysis / methods: statistical
© ESO, 2016
1. Introduction
This paper, one of a series associated with the 2015 release of Planck^{1} data, will present limits on departures from the global isotropy of spacetime. We assess anisotropic but homogeneous Bianchi cosmological models and nontrivial global topologies in the light of the latest temperature and polarization data.
In Planck Collaboration XXVI (2014), the limits came from the 2013 Planck cosmological data release: cosmic microwave background (CMB) intensity data collected over approximately one year. This work uses the 2015 Planck data: CMB intensity from the whole mission along with a subset of polarization data. The greater volume of intensity data will allow more restrictive limits on the possibility of topological scales that are slightly larger than the volume enclosed by the lastscattering surface (roughly the Hubble volume), probing the excess anisotropic correlations that would be induced at large angular scales were such a model to obtain. For cubic torus topologies, we can therefore observe explicit repeated patterns (matched circles) when the comoving length of an edge is less than twice the distance to the recombination surface, (using units with c = 1 here and throughout). Polarization, on the other hand, which is largely generated during recombination itself, can provide a more sensitive probe of topological domains smaller than the Hubble volume.
Whereas the analysis of temperature data in multiply connected universes has been treated in some depth in the literature (see Planck Collaboration XXVI 2014, and references therein), the discussion of polarization has been less complete. This paper therefore extends our previous likelihood analysis to polarized data, updates the direct search for matched circles (Cornish et al. 2004) as discussed in Bielewicz et al. (2012), and uses these to present the first limits on global topology from polarized CMB data.
The cosmological properties of Bianchi models (Collins & Hawking 1973; Barrow et al. 1985), were initially discussed in the context of CMB intensity (Barrow 1986; Jaffe et al. 2006c,a; Pontzen 2009). As discussed in Planck Collaboration XXVI (2014), it is by now well known that the observed largescale intensity pattern mimics that of a particular Bianchi VII_{h} model, albeit one with cosmological parameters that are quite different from those needed to reproduce other CMB and cosmological data. More recently the induced polarization patterns have been calculated (Pontzen & Challinor 2007; Pontzen 2009; Pontzen & Challinor 2011). In this paper, we analyse the complete Planck intensity data, and compare the polarization pattern induced by that anisotropic model to Planck polarization data.
We note that the lack of a strong detection of cosmic Bmode polarization already provides some information about the Bianchi models: the induced geometrical focusing does not distinguish between E and B and thus should produce comparable amounts of each (e.g., Pontzen 2009). This does not apply to topological models: the linear evolution of primordial perturbations guarantees that a lack of primordial tensor perturbations results in a lack of Bmode polarization – the transfer function is not altered by topology.
In Sect. 2, we discuss previous limits on anisotropic models from Planck and other experiments. In Sect. 3 we discuss the CMB signals generated in such models, generalized to both temperature and polarization. In Sect. 4 we describe the Planck data and simulations we use in this study, the different methods we apply to those data, and the validation checks performed on those simulations. In Sect. 5 we discuss the results and conclude in Sect. 6 with the outlook for application of these techniques to future data and broader classes of models.
2. Previous results
The first searches for nontrivial topology on cosmic scales looked for repeated patterns or individual objects in the distribution of galaxies (Sokolov & Shvartsman 1974; Fang & Sato 1983; Fagundes & Wichoski 1987; Lehoucq et al. 1996; Roukema 1996; Weatherley et al. 2003; Fujii & Yoshii 2011). Searches for topology using the CMB began with COBE (Bennett et al. 1996) and found no indications of a nontrivial topology on the scale of the lastscattering surface (e.g., Starobinskij 1993; Sokolov 1993; Stevens et al. 1993; de OliveiraCosta & Smoot 1995; Levin et al. 1998; Bond et al. 1998b, 2000; Rocha et al. 2004; but see also Roukema 2000b,a). With the higher resolution and sensitivity of WMAP, there were indications of low power on large scales which could have had a topological origin (Jarosik et al. 2011; Luminet et al. 2003; Caillerie et al. 2007; Aurich 1999; Aurich et al. 2004, 2005, 2006, 2008; Aurich & Lustig 2013; Lew & Roukema 2008; Roukema et al. 2008; Niarchou et al. 2004), but this possibility was not borne out by detailed real and harmonicspace analyses in two dimensions (Cornish et al. 2004; Key et al. 2007; Bielewicz & Riazuelo 2009; Dineen et al. 2005; Kunz et al. 2006; Phillips & Kogut 2006; Niarchou & Jaffe 2007). Most studies, including this work, have emphasized searches for fundamental domains with antipodal correlations; see Vaudrevange et al. (2012) for results from a general search for the patterns induced by nontrivial topology on scales within the volume defined by the lastscattering surface, and, for example, Aurich & Lustig (2014) for a recent discussion of other possible topologies.
For a more complete overview of the field, we direct the reader to Planck Collaboration XXVI (2014). In that work, we applied various techniques to the Planck 2013 intensity data. For topology, we showed that a fundamental topological domain smaller than the Hubble volume is strongly disfavoured. This was done in two ways: first, a direct likelihood calculation of specific topological models; and second, a search for the expected repeated “circles in the sky” (Cornish et al. 2004), calibrated by simplyconnected simulations. Both of these showed that the scale of any possible topology must exceed roughly the distance to the lastscattering surface, χ_{rec}. For the cubic torus, we found that the radius of the largest sphere inscribed in the topological fundamental domain must be ℛ_{i}> 0.92 χ_{rec} (at loglikelihood ratio Δlnℒ > −5 relative to a simplyconnected flat Planck 2013 bestfit model). The matchedcircle limit on topologies predicting backtoback circles was ℛ_{i}> 0.94 χ_{rec} at the 99% confidence level.
Prior to the present work, there have been some extensions of the search for cosmic topology to polarization data. In particular, Bielewicz et al. (2012; see also Riazuelo et al. 2006) extended the direct search for matched circles to polarized data and found that the available WMAP data had insufficient sensitivity to provide useful constraints.
For Bianchi VII_{h} models, in Planck Collaboration XXVI (2014) a full Bayesian analysis of the Planck 2013 temperature data was performed, following the methods of McEwen et al. (2013). It was concluded that a physicallymotivated model was not favoured by the data. If considered as a phenomenological template (for which the parameters common to the standard stochastic CMB and the deterministic Bianchi VII_{h} component are not linked), it was shown that an unphysical Bianchi VII_{h} model is favoured, with a logBayes factor between 1.5 ± 0.1 and 2.8 ± 0.1 – equivalent to an odds ratio of between approximately 1:4 and 1:16 – depending of the component separation technique adopted. Prior to the analysis of Planck Collaboration XXVI (2014), numerous analyses of Bianchi models using COBE (Bennett et al. 1996) and WMAP (Jarosik et al. 2011) data had been performed (Bunn et al. 1996; Kogut et al. 1997; Jaffe et al. 2005, 2006a,c,b; Cayón et al. 2006; Land & Magueijo 2006; McEwen et al. 2006; Bridges et al. 2007; Ghosh et al. 2007; Pontzen & Challinor 2007; Bridges et al. 2008; McEwen et al. 2013), and a similar Bianchi template was found in the WMAP data, first by Jaffe et al. (2005) and then subsequently by others (Bridges et al. 2007; Bridges et al. 2008; McEwen et al. 2013). Pontzen & Challinor (2007) discussed the CMB polarization signal from Bianchi models, and showed some incompatibility with WMAP data due to the large amplitude of both E and Bmode components. For a more detailed review of the analysis of Bianchi models we refer the reader to Planck Collaboration XXVI (2014).
3. CMB signals in anisotropic and multiplyconnected universes
3.1. Topology
There is a long history of studying the possible topological compactification of FriedmannLemaîtreRobertsonWalker (FLRW) cosmologies; we refer readers to overviews such as Levin (2002), LachiezeRey & Luminet (1995), and Riazuelo et al. (2004a,b) for mathematical and physical detail. The effect of a nontrivial topology is equivalent to considering the full (simplyconnected) three dimensional spatial slice of the manifold (the covering space) as being tiled by identical repetitions of a shape which is finite in one or more directions, the fundamental domain. In flat universes, to which we specialize here, there are a finite number of possibilities, each described by one or more continuous parameters describing the size in different directions.
In this paper, we pay special attention to topological models in which the fundamental domain is a rightrectangular prism (the threetorus, also referred to as “T3”), possibly with one or two infinite dimensions (the T2 “chimney” or “rod”, and T1 “slab” models). We limit these models in a number of ways. We explicitly compute the likelihood of the length of the fundamental domain for the cubic torus. Furthermore, we consider the slab model as a proxy for other models in which the matched circles (or excess correlations) are antipodally aligned, similar to the “lens” spaces available in manifolds with constant positive curvature. These models are thus sensitive to tori with varying side lengths, including those with nonrightangle corners. In these cases, the likelihood would have multiple peaks, one for each of the aligned pairs; their sizes correspond to those of the fundamental domains and their relative orientation to the angles. These nonrectangular prisms will be discussed in more detail in Jaffe & Starkman (in prep.).
3.1.1. Computing the covariance matrices
In Planck Collaboration XXVI (2014) we computed the temperaturetemperature (TT) covariance matrices by summing up all modes k_{n} that are present given the boundary conditions imposed by the nontrivial topology. For a cubic torus, we have a threedimensional wave vector k_{n} = (2π/L)n for a triplet of integers n, with unit vector and the harmonicspace covariance matrix (1)where is the temperature radiation transfer function (see, e.g., Bond & Efstathiou 1987; and Seljak & Zaldarriaga 1996).
It is straightforward to extend this method to include polarization, since the cubic topology affects neither the local physics that governs the transfer functions, nor the photon propagation. The only effect is the discretization of the modes. We can therefore simply replace the radiation transfer function for the temperature fluctuations with the one for polarization, and obtain (2)where X,X′ = E,T. We are justified in ignoring the possibility of Bmode polarization as it is sourced only by primordial gravitational radiation even in the presence of nontrivial topology. In this way we obtain three sets of covariance matrices: TT, TE, and EE. In addition, since the publication of Planck Collaboration XXVI (2014) we have optimized the cubic torus calculation by taking into account more of the symmetries. The resulting speedup of about an order of magnitude allowed us to reach a higher resolution of ℓ_{max} = 64.
The fiducial cosmology assumed in the calculation of the covariance matrices is a flat ΛCDM FLRW Universe with Hubble constant H_{0} = 100h km s^{1} Mpc^{1}, where: h = 0.6719; scalar spectral index n_{s} = 0.9635; baryon density Ω_{b}h^{2} = 0.0221; cold dark matter density Ω_{c}h^{2} = 0.1197; and neutrino density Ω_{ν}h^{2} = 0.0006.
3.1.2. Relative information in the matrices
To assess the information content of the covariance matrices, we consider the KullbackLeibler (KL) divergence (see, e.g., Kunz et al. 2006, 2008; and Planck Collaboration XXVI 2014, for further applications of the KL divergence to topology). The KL divergence between two probability distributions p_{1}(x) and p_{2}(x) is given by (3)If the two distributions are Gaussian with covariance matrices C_{1} and C_{2}, this expression simplifies to (4)and is thus an asymmetric measure of the discrepancy between the covariance matrices. The KL divergence can be interpreted as the ensemble average of the loglikelihood ratio Δlnℒ between realizations of the two distributions. Hence, it enables us to probe the ability to tell if, on average, we can distinguish realizations of p_{1} from a fixed p_{2} without having to perform a bruteforce Monte Carlo integration. Thus, the KL divergence is related to ensemble averages of the likelihoodratio plots that we present for simulations (Sect. 4.4.1) and real data (Sect. 5.1) but can be calculated from the covariance matrices alone. Note that with this definition, the KL divergence is minimized for cases with the best match (maximal likelihood).
In Planck Collaboration XXVI (2014) we used the KL divergence to show that the likelihood is robust to differences in the cosmological model and small differences in the topology.
In Fig. 1 we plot the KL divergence relative to an infinite Universe for the slab topology as a function of resolution ℓ_{max} (upper panel) and fundamental domain size (lower panel). Our ability to detect a topology with a fundamental domain smaller than the distance to the lastscattering surface (approximately at the horizon distance , so with sides of length ) grows significantly with the resolution even beyond the cases that we studied. For the noise levels of the 2015 lowP data considered here and defined in Planck Collaboration XIII (2016), polarization maps do not add much information beyond that contained in the temperature maps, although, as also shown in Sect. 4.4.1, the higher sensitivity achievable by the full Planck lowℓ data over all frequencies should enable even stronger constraints on these small fundamental domains.
If, however, the fundamental domain is larger than the horizon (as is the case for ) then the relative information in the covariance matrix saturates quite early and a resolution of ℓ_{max} ≃ 48 is actually sufficient. The main goal is thus to ensure that we have enough discriminatory power right up to the horizon size. In addition, polarization does not add much information in this case, irrespective of the noise level. This is to be expected: polarization is only generated for a short period of time around the surface of last scattering. Once the fundamental domain exceeds the horizon size, the relative information drops rapidly towards zero, and the dependence on ℓ_{max} becomes weak.
In Fig. 2 we plot the KL divergence as a function of the size of the fundamental domain for fixed cube (T3), rod (T2), and slab (T1) topologies, each with fundamental domain size , compared to the slab. Each shows a strong dip at , indicating the ability to detect this topology (although note the presence of a weaker dip around half the correct size, ). The figure also shows that ℓ_{max} = 40 still shows the dip at the correct location, although somewhat more weakly than ℓ_{max} = 80.
Note that the shape of the curves is essentially identical, with the slab likelihood able to detect one or more sets of antipodal matched circles (and their related excess correlations at large angular scales) present in each case. Figure 2 therefore shows that using the covariance matrix for a slab (T1) topology also allows detection of rod (T2) and cubic (T3) topologies: this is advantageous as the slab covariance matrix is considerably easier to calculate than the cube and rod, since it is only discretized in a single direction. Figure 3 shows the KL divergence as a function of the relative rotation of the fundamental domain, showing that, despite the lack of the full set of three pairs of antipodal correlations, we can determine the relative rotation of a single pair. This is exactly how the matchedcircles tests work. Furthermore, as we will demonstrate in Sect. 4.4.1, slab likelihoods are indeed separately sensitive to the different sets of antipodal circles in cubic spaces. We can hence adopt the slab as the most general tool for searching for spaces with antipodal circles.
Fig. 1 KL divergence of slab (T1) topologies relative to an infinite Universe as a function of ℓ_{max} with sizes and (top), and as a function of size L of the fundamental domain for various ℓ_{max} (bottom). A torus with , corresponding to (H_{0}L)^{1}< 0.154, has a fundamental domain that is larger than the distance to the lastscattering surface and leaves only a small trace in the CMB. This is why the KL divergence drops rapidly at this point. Note that the information for continues to rise with ℓ_{max} whereas it levels off for the slightly larger case. In the lower panel we see that there is a slight feature in d_{KL} at about half the horizon distance, which is probably due to harmonic effects. The corresponding figures for cubic (T3) topologies look qualitatively similar except that all d_{KL} values are three times larger. 
Fig. 2 KL divergence of fixed cubic, rod, and slab topologies with fundamental domain side compared to a slab of variable fundamental domain size L. The chimney space T2 dates from the 2013 analysis (Planck Collaboration XXVI 2014) and was computed for the bestfit parameters of that release. In all cases the smallest KL divergence, corresponding to the best fit, appears at , indicating that the slab space can be used to detect other topologies. An additional dip at L ≃ 5.5 / (2H_{0}) may be due to a harmonic effect at half the size of the fundamental domain; it is, however, much smaller than the drop in KL divergence at the size of the fundamental domain. 
Fig. 3 KL divergence of a slab space relative to a cubic topology, as a function of rotation angle of the slab space (blue curve). Both spaces have and ℓ_{max} = 80. The horizontal black dashed line gives the KL divergence of an infinite Universe relative to the cubic topology and illustrates how much better the slab space fits with the correct orientation relative to the cubic torus. 
3.2. Bianchi models
The polarization properties of Bianchi models were first derived in Pontzen & Challinor (2007) and extensively categorized in Pontzen (2009) and Pontzen & Challinor (2011). In these works it was shown that advection in Bianchi universes leads to efficient conversion of Emode polarization to B modes; evidence for a significant Bianchi component found in temperature data would therefore suggest a large Bmode signal (but not necessarily require it; see Pontzen 2009). For examples of the temperature and polarization signatures of Bianchi VII_{h} models we refer the reader to Fig. 1 of Pontzen (2009). Despite the potential for CMB polarization to constrain the Bianchi sector, a full polarization analysis has not yet been carried out. The analysis of Pontzen & Challinor (2007) remains the stateoftheart, where WMAP BB and EB power spectra were used to demonstrate (using a simple χ^{2} analysis) that a Bianchi VII_{h} model derived from temperature data was disfavoured compared to an isotropic model.
The subdominant, deterministic CMB contributions of Bianchi VII_{h} models can be characterized by seven parameters: the matter and dark energy densities, Ω_{m} and Ω_{Λ}, respectively; the present dimensionless vorticity, (ω/H)_{0}; the dimensionless lengthscale parameter, x, which controls the “tightness” of the characteristic Bianchi spirals; and the Euler angles^{2}, (α,β,γ), describing their orientation (i.e., the choice of coordinate system), where H is the Hubble parameter. For further details see Planck Collaboration XXVI (2014), McEwen et al. (2013), Pontzen (2009), Pontzen & Challinor (2007), Jaffe et al. (2006c), Jaffe et al. (2005), and Barrow et al. (1985).
4. Methods
4.1. Data
In this work we use data from the Planck 2015 release. This includes intensity maps from the full mission, along with a subset of polarization data. Specifically, for the likelihood calculations discussed below (Sect. 4.3.1 for application to topology and Sect. 4.3.3 for Bianchi models) which rely on HEALPix maps at N_{side} = 16, we use the data designated “lowT,P”, as defined for the lowℓPlanck likelihood for isotropic models (Planck Collaboration XI 2016; Planck Collaboration XIII 2016): lowP polarization maps based on the LFI 70 GHz channel and lowT temperature maps created by the Commander component separation method, along with the appropriate mask and noise covariance matrix. As in Planck Collaboration XXVI (2014), the intensity noise contribution is negligible on these scales, and diagonal regularizing noise with variance has therefore been added to the intensity portion of the noise covariance matrix. We cut contaminated regions of the sky using the lowℓ mask defined for the Planck isotropic likelihood code (Planck Collaboration XI 2016), retaining 94% of the sky for temperature, and the lowT,P polarization mask, cleaned with the templates created from Planck 30 GHz and 353 GHz data, retaining 47% of the sky for polarization.
The matchedcircle search (Sects. 4.2 and 5.1.1) uses four componentseparated maps (Planck Collaboration IX 2016) which effectively combine both intensity and polarization information from different scales. The maps are smoothed with a Gaussian filter of 30′ and 50′ full width at half maximum (FWHM) for temperature and polarization, respectively, and degraded to N_{side} = 512. Corresponding temperature and polarization common masks for diffuse emission, with a point source cut for the brightest sources, downgraded analogously to the maps, are used. After degradation, and accounting for the needed expansion of the polarization mask due to the conversion of Q and U to E, the temperature map retains 74% of the sky and the polarization map 40%. These Emode maps are calculated using the method of Bielewicz et al. (2012; see also Kim 2011) and correspond to the spherical Laplacian of the scalar E, consequently filtering out power at large angular scales.
4.2. Topology: matched circles
As in Planck Collaboration XXVI (2014), we use the circle comparison statistic of Cornish et al. (1998), optimized for smallscale anisotropies (Cornish et al. 2004), to search for correlated circles in sky maps of the CMB temperature and polarization anisotropy. The circle comparison statistic uses the fact that the intersection of the topological fundamental domain with the surface of last scattering is a circle, potentially viewed from different directions in a multiplyconnected Universe. Contrary to the temperature anisotropy, sourced by multiple terms at the lastscattering surface (i.e., the internal photon density fluctuations combined with the ordinary SachsWolfe and Doppler effects), the CMB polarization anisotropy is sourced only by the quadrupole distribution of radiation scattering from free electrons at the moment of recombination (e.g., Kosowsky 1996). In particular, the recombination signal from polarization is only generated for a short time while there are enough electrons to scatter the photons but few enough for the plasma to be sufficiently transparent. Thus, in a multiconnected Universe the polarization signal does not exhibit the same cancellation of contributions from different terms as in the temperature anisotropy (Bielewicz et al. 2012). Polarization thus can provide a better opportunity for the detection of topological signatures than a temperature anisotropy map. There is a small subtlety here: whereas the intensity is a scalar and thus is unchanged when viewed from different directions, the polarization is a tensor which behaves differently under rotation. The polarization pattern itself depends on the viewing angle; hence, we need to use the coordinateindependent quantities, E and B, which are scalars (or pseudoscalars) and are thus unchanged when viewed from different directions.
The decomposition into E and B of an arbitrary masked CMB polarization map, contaminated by noise, foregrounds, and systematic errors, is itself a computationally demanding task, nonlocal on the sky. Assuming negligible initial B polarization, we use only the E maps produced from componentseparated CMB polarization maps using the same approach as Bielewicz et al. (2012).
Compared with the likelihood method described below, the circles search uses higherresolution maps, and thus is sensitive out to a much higher maximum multipole, ℓ_{max}. It is also potentially less sensitive to largescale systematic errors, as the lowest multipoles are effectively filtered out: the polarization signal is weighted by a factor proportional to ℓ^{2} in the transformation from the Stokes parameters Q and U to an Emode map. From the results of Sect. 3.1.2, this indicates that it uses more of the information available when confronting models with fundamental domains within the lastscattering surface compared to our implementation of the likelihood, limited to ℓ_{max} ≃ 40. As we show in Sect. 4.4.1, this also allows the use of highpass filtered componentseparated maps (as defined in Planck Collaboration IX 2016) without a significant decrease in the ability to detect a multiplyconnected topology.
The matchedcircle statistic is defined by (5)where X_{i,m} and X_{j,m} denote the Fourier coefficients of the temperature or Emode fluctuations around two circles of angular radius α centred at different points on the sky, i and j, respectively, with relative phase φ_{∗}. The mth harmonic of the field anisotropies around the circle is weighted by the factor  m , taking into account the number of degrees of freedom per mode. Such weighting enhances the contribution of smallscale structure relative to largescale fluctuations.
The S^{+} statistic corresponds to pairs of circles with the points ordered in a clockwise direction (phased). For the alternative ordering, when the points are ordered in an anticlockwise direction (antiphased along one of the circles), the Fourier coefficients X_{i,m} are complex conjugated, defining the S^{−} statistic. This allows the detection of both orientable and nonorientable topologies. For orientable topologies the matched circles have antiphased correlations, while for nonorientable topologies they have a mixture of antiphased and phased correlations.
The S^{±} statistics take values over the interval [− 1,1]. Circles that are perfectly matched have S = 1, while uncorrelated circles will have a mean value of S = 0. To find matched circles for each radius α, the maximum value is determined.
Because general searches for matched circles are computationally very intensive, we restrict our analysis to a search for pairs of circles centred around antipodal points, so called backtoback circles. The maps are also downgraded as described in Sect. 4.1. This increases the signaltonoise ratio and greatly speeds up the computations required, but with no significant loss of discriminatory power. Regions most contaminated by Galactic foreground were removed from the analysis using the common temperature or polarization mask. More details on the numerical implementation of the algorithm can be found in Bielewicz & Banday (2011) and Bielewicz et al. (2012).
To draw any conclusions from an analysis based on the statistic , it is very important to correctly estimate the threshold for a statistically significant match of circle pairs. We used 300 Monte Carlo simulations of the PlanckSMICA maps processed in the same way as the data to establish the threshold such that fewer than 1% of simulations would yield a false event. Note that we perform the entire analysis, including the final statistical calibration, separately for temperature and polarization.
4.3. Likelihood
4.3.1. Topology
For the likelihood analysis of the largeangle intensity and polarization data we have generalized the method implemented in Planck Collaboration XXVI (2014) to include polarization. The likelihood, i.e., the probability to find a combined temperature and polarization data map d with associated noise matrix N given a certain topological model T is then given by (6)where now d is a 3N_{pix}component data vector obtained by concatenation of the (I,Q,U) data sets while C and N are 3N_{pix} × 3N_{pix} theoretical signal and noise covariance matrices, arranged in block form as (7)Finally, Θ_{C} is the set of standard cosmological parameters, Θ_{T} is the set of topological parameters (e.g., the size, L, of the fundamental domain), ϕ is the orientation of the topology (e.g., the Euler angles), and A is a single amplitude, scaling the signal covariance matrix (this is equivalent to an overall amplitude in front of the power spectrum in the isotropic case). Working in pixel space allows for the straightforward application of an arbitrary mask, including separate masks for intensity and polarization parts of the data. The masking procedure can also be used to limit the analysis to intensity or polarization only.
Since C + N in pixel space is generally poorly conditioned, we again (following the 2013 procedure) project the data vector and covariance matrices onto a limited set of orthonormal basis vectors, select N_{m} such modes for comparison, and consider the likelihood marginalized over the remainder of the modes, (8)where C and N are restricted to the N_{m} × N_{m} subspace.
The choice of the basis modes and their number N_{m} used for analysis is a compromise between robust invertibility of C + N and the amount of information retained. All the models for which likelihoods are compared must be expanded in the same set of modes. Thus, in Planck Collaboration XXVI (2014) we used the set of eigenmodes of the cutsky covariance matrix of the fiducial bestfit simplyconnected Universe, C_{fid}, as the analysis basis, limiting ourselves to the N_{m} modes with the largest eigenvalues. For comparison with the numbers we use, a fullsky temperature map with maximum multipole ℓ_{max} has (ℓ_{max} + 1)^{2}−4 independent modes (four are removed to account for the unobserved monopole and dipole).
The addition of polarization data, with much lower signaltonoise than the temperature, raises a new question: how is the temperature and polarization data mix reflected in the limited basis set we project onto? The most natural choice is the set of eigenmodes of the signaltonoise matrix C_{fid}N^{1} for the fiducial model, and a restriction of the mode set based on signaltonoise eigenvalues (see, e.g., Bond et al. 1998a). This, however, requires robust invertibility of the noise covariance matrix, which, again, is generally not the case for the smoothed data. Moreover, such a ranking by S/N would inevitably favour the temperature data, and we wish to explore the effect of including polarization data on an equal footing with temperature. We therefore continue to use the eigenmodes of the cutsky fiducial covariance matrix as our basis. By default, we select the first N_{m} = 1085 eigenmodes (corresponding to ℓ_{max} = 32), though we vary the mode count where it is informative to do so.
In Fig. 4 we show I, Q, and U maps of the highesteigenvalue (i.e., highest contribution to the signal covariance) mode for our fiducial simplyconnected model. Note that the scale is different for temperature compared to the two polarization maps: the temperature contribution to the mode is much greater than that of either polarization component. We show modes for the masked sky, although in fact the structure at large scales is similar to the fullsky case, rotated and adjusted somewhat to account for the mask. In Fig. 5 we show the structure of mode 301, with much lower signal amplitude (this particular mode was selected at random to indicate the relative ratios of temperature, polarization, and noise). Temperature remains dominant, although polarization begins to have a greater effect. Note that at this level of signal amplitude, the pattern is aligned with the mask, and shows a strong correlation between temperature and polarization.
Fig. 4 Mode structure plotted as maps for the eigenvector corresponding to the highestsignal eigenvalue of the fiducial simplyconnected model. The top map corresponds to temperature, middle to Q polarization, and bottom to U polarization. Masked pixels are plotted in grey. 
Fig. 5 Mode structure plotted as maps for the eigenvector corresponding to the 301sthighestsignal eigenvalue of the fiducial simplyconnected model. The top map corresponds to temperature, middle to Q polarization, and bottom to U polarization. Masked pixels are plotted in grey. 
4.3.2. Evaluating the topological likelihood
The aim of the topological likelihood analysis is to calculate the likelihood as a function of the parameters pertaining to a particular topology, p(d  Θ_{T},T). To do so, we must marginalize over the other parameters appearing in Eq. (8), namely Θ_{C}, ϕ, and A, as (9)The complexity of the topological covariance matrix calculation precludes a joint examination of the full cosmological and topological parameter spaces. Instead, we adopt the deltafunction prior to fix the cosmological parameters at their fiducial values, (as defined in Sect. 3.1.1), and evaluate the likelihood on a grid of topological parameters using a restricted set of precalculated covariance matrices. We note that, as discussed in Planck Collaboration XXVI (2014), the ability to detect or rule out a multiply connected topology is insensitive to the values of the cosmological parameters adopted for the calculation of the covariance matrices.
In the setting described above, Eq. (9) simplifies to (10)where the likelihood at each gridpoint in topological parameter space, Θ_{T,i}, is equal to the probability of obtaining the data given fixed cosmological and topological parameters and a compactification (i.e., fundamental domain shape and size), marginalized over orientation and amplitude. The calculation therefore reduces to evaluating the Bayesian evidence for a set of gridded topologies. As we focus on cubic torus and slab topologies in this work, we note that the sole topological parameter of interest is the size of the fundamental domain, L.
Even after fixing the cosmological parameters, calculating the Bayesian evidence is a timeconsuming process, and is further complicated by the multimodal likelihood functions typical in nontrivial topologies. We therefore approach the problem on two fronts. We first approximate the likelihood function using a “profile likelihood” approach, as presented in Planck Collaboration XXVI (2014), in which the marginalization in Eq. (10) is replaced with maximization in the fourdimensional space of orientation and amplitude parameters. Specifically, we maximize the likelihood over the three angles defining the orientation of the fundamental domain using a threedimensional Amoeba search (e.g., Press et al. 1992), where at each orientation the likelihood is separately maximized over the amplitude. Due to the complex structure of the likelihood surface in orientation space, we repeat this procedure five times with different starting orientations. This number of repetitions was chosen as a compromise between computational efficiency and assurance of statistical robustness, after testing of various strategies for the number of repetitions and the distribution of starting points, along with explicit extra runs to test outliers. To ensure uniform and nondegenerate coverage, the orientation space is traversed in a Cartesian projection of the northern hemisphere of the threesphere S3 representation of rotations.
The profile likelihood calculation allows rapid evaluation of the likelihood and testing of different models compared with a variety of data and simulations, but it is difficult to interpret in a Bayesian setting. As we show below, however, the numerical results of profiling over this limited set of parameters agree numerically very well with the statistically correct marginalization procedure.
Our second approach explicitly calculates the marginalized likelihood, Eq. (10), allowing full Bayesian inference at the cost of increased computation time. We use the public MultiNest^{3} code (Feroz & Hobson 2008; Feroz et al. 2009, 2013) − optimized for exploring multimodal probability distributions in tens of dimensions − to compute the desired evidence values via nested sampling (Skilling 2004). MultiNest is run in its importance nested sampling mode (Feroz et al. 2013) using 200 live points, with tolerance and efficiency set to their recommended values of 0.5 and 0.3, respectively. The final ingredient needed to calculate the evidence values are priors for the marginalized parameters. We use a log prior for the amplitude, truncated to the range 0.1 ≤ A ≤ 10, and the Euler angles are defined to be uniform in 0 ≤ α< 2π, −1 ≤ cosβ ≤ 1, and 0 ≤ γ< 2π, respectively; MultiNest is able to wrap the priors on α and γ. The combined code will be made public as part of the AniCosmo^{4} package (McEwen et al. 2013).
It is worth noting that this formalism can be extended to compare models with different compactifications (or the simply connected model) using Bayesian model selection: the only additional requirements are priors for the topological parameters. Taking the current slab and cubic torus topologies as examples, by defining a prior on the size of the fundamental domain one can calculate the evidence for each model. Assuming each topology is equally likely a priori, i.e., that p(T_{slab}) = p(T_{cub}), one can then write down the relative probability of the two topologies given the data: (11)Unfortunately, it is difficult to provide a physicallymotivated proper prior distribution for the size of the fundamental domain. Even pleading ignorance and choosing a “naïve” uniform prior would require an arbitrary upper limit to L whose exact value would strongly influence the final conclusion. For this reason, we refrain from extending the formalism to model selection within this manuscript.
4.3.3. Bianchi models
While physically the cosmological densities describing Bianchi models should be identified with their standard ΛCDM counterparts, in previous analyses unphysical models have been considered in which the densities are allowed to differ. The first coherent analysis of Bianchi VII_{h} models was performed by McEwen et al. (2013), where the ΛCDM and Bianchi densities are coupled and all cosmological and Bianchi parameters are fit simultaneously. In the analysis of Planck Collaboration XXVI (2014), in order to compare with all prior studies both coupled and decoupled models were analysed. We consider the same two models here: namely, the physical opencoupledBianchi model where an open cosmology is considered (for consistency with the open Bianchi VII_{h} models), in which the Bianchi densities are coupled to their standard cosmological counterparts; and the phenomenological flatdecoupledBianchi model where a flat cosmology is considered and in which the Bianchi densities are decoupled.
We firstly carry out a full Bayesian analysis for these two Bianchi VII_{h} models, repeating the analysis performed in Planck Collaboration XXVI (2014) with updated Planck temperature data. The methodology is described in detail in McEwen et al. (2013) and summarized in Planck Collaboration XXVI (2014). The complete posterior distribution of all Bianchi and cosmological parameters is sampled and Bayesian evidence values are computed to compare Bianchi VII_{h} models to their concordance counterparts. Bianchi temperature signatures are simulated using the Bianchi2^{5} code (McEwen et al. 2013), while the AniCosmo code is used to perform the analysis, which in turn uses MultiNest to sample the posterior distribution and compute evidence values.
To connect with polarization data, we secondly analyse polarization templates computed using the bestfit parameters from the analysis of temperature data. For the resulting small set of bestfit models, polarization templates are computed using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and have been provided by Pontzen (priv. comm.). These Bianchi VII_{h} simulations are more accurate than those considered for the temperature analyses performed here and in previous works (see, e.g., Planck Collaboration XXVI 2014; McEwen et al. 2013; Bridges et al. 2008; Bridges et al. 2007; Jaffe et al. 2005, 2006a,b,c), since the recombination history is modelled. The overall morphology of the patterns are consistent between the codes; the strongest effect of incorporating the recombination history is its impact on the polarization fraction, although the amplitude of the temperature component can also vary by approximately 5% (which is calibrated in the current analysis, as described below).
Using the simulated Bianchi VII_{h} polarization templates computed following Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.), we perform a maximumlikelihood fit for the amplitude of these templates using Planck polarization data (a full Bayesian evidence calculation of the complete temperature and polarization data set incorporating the more accurate Bianchi models of Pontzen & Challinor 2007; and Pontzen 2009, is left to future work). The likelihood in the Bianchi scenario is identical to that considered in Planck Collaboration XXVI (2014) and McEwen et al. (2013); however, we now consider the Bianchi and cosmological parameters fixed and simply introduce a scaling of the Bianchi template. The resulting likelihood reads: (12)where d denotes the data vector, is the Bianchi template for bestfit Bianchi parameters , C = C(Θ_{C}) is the cosmological covariance matrix for the bestfit cosmological parameters , N is the noise covariance, and λ is the introduced scaling parameter (the effective vorticity of the scaled Bianchi component is simply λ(ω/H)_{0}).
In order to effectively handle noise and partial sky coverage the data are analysed in pixel space. We restrict to polarization data only here since temperature data are used to determine the bestfit Bianchi parameters. The data and template vectors thus contain unmasked Q and U Stokes components only and, correspondingly, the cosmological and noise covariance matrices are given by the polarization (Q and U) subspace of Eq. (7), and again contain unmasked pixels only.
The maximumlikelihood (ML) estimate of the template amplitude is given by and its dispersion by Δλ^{ML} = [t^{†}(C + N)^{1}t] ^{− 1 / 2} (see, e.g., Kogut et al. 1997; Jaffe et al. 2005). If Planck polarization data support the bestfit Bianchi model found from the analysis of temperature data we would expect λ^{ML} ≃ 1. A statistically significant deviation from unity in the fitted amplitude can thus be used to rule out the Bianchi model using polarization data.
As highlighted above, different methods are used to simulate Bianchi temperature and polarization components, where the amplitude of the temperature component may vary by a few percent between methods. We calibrate out this amplitude mismatch by scaling the polarization components by a multiplicative factor fitted so that the temperature components simulated by the two methods match, using a maximumlikelihood template fit again, as described above.
4.4. Simulations and validation
4.4.1. Topology
Matched circles.
Before beginning the search for pairs of matched circles in the Planck data, we validate our algorithm using the same simulations as employed for the Planck Collaboration XXVI (2014) and Bielewicz et al. (2012) papers, i.e., the CMB sky for a Universe with a threetorus topology for which the dimension of the cubic fundamental domain is , well within the lastscattering surface. We computed the a_{ℓm} coefficients up to the multipole of order ℓ = 500 and convolved them with the same smoothing beam profile as used for the PlanckSMICA map. To the map was added noise corresponding to the SMICA map. In particular, we verified that our code is able to find all pairs of matched circles in such a map. The statistic for the Emode map is shown in Fig. 6.
Fig. 6 Example statistic as a function of circle radius α for a simulated CMB Emode map of a Universe with the topology of a cubic 3torus with dimensions . To this map noise was added corresponding to the SMICA map. The thick overlapping curves show the statistic for simulated polarization maps with angular resolution and noise level corresponding to the PlanckSMICA map for three cases: without highpass filtering (solid, red); with filtering (dashed, blue); and with a 5% larger amplitude of noise (dotdashed, green). The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the highpass filtered PlanckSMICA polarization maps, smoothed and masked in the same way as the data, would yield a false event. 
Because for the baseline analysis we use highpass filtered maps, we also show the analysis of the SMICAEmode map highpass filtered so that the lowest order multipoles (ℓ< 20) are removed from the map (the multipoles in the range 20 ≤ ℓ ≤ 40 are apodized between 0 and 1 using a cosine as defined in Planck Collaboration IX 2016). The highpass filtering does not decrease our ability to detect a multiplyconnected topology using the matchedcircle method. This is consistent with the negligible sensitivity of the matchedcircle statistic to the reionization signal, studied by Bielewicz et al. (2012). This is a consequence of the weighting of the polarization data by a factor proportional to ℓ^{2} employed in the transformation from the Stokes parameters Q and U to an Emode map, which effectively filters out the largestscale multipoles from the data. This test shows that the matchedcircle method, contrary to the likelihood method, predominantly exploits the topological signal in the CMB anisotropies at moderate angular scales.
We also checked robustness of detection with respect to noise level in order to account for small discrepancies between the noise level in the Planck FFP8 simulations and the 2015 data (Planck Collaboration XII 2016). We repeated the analysis for the highpass filtered map with added noise with 5% larger amplitude than for the original map. As we can see in Fig. 6, the statistic changes negligibly.
The intersection of the peaks in the matching statistic with the false detection level estimated for the CMB map corresponding to the simplyconnected Universe defines the minimum radius of the correlated circles that can be detected for this map. We estimate the minimum radius by extrapolating the height of the peak with radius 18° seen in Fig. 6 towards smaller radii. This allows for a rough estimation of the radius, with a precision of a few degrees. However, better precision is not required, because for small minimum radius (as obtains here) constraints on the size of the fundamental domain are not very sensitive to differences of the minimum radius of order a few degrees. As we can see in Fig. 6, the minimum radius α_{min} takes a value in the range from 10° to around 15°. To be conservative we use the upper end of this range for the computation of constraints on the size of the fundamental domain, and we thus take α_{min} ≃ 15°.
Likelihood.
To validate and compare the performance of the two likelihood methods, we perform two sets of tests: a null test using a simulation of a simply connected Universe, and a signal test using a simulation of a toroidal Universe with . The two test maps are generated at N_{side} = 16 and are bandlimited using a 640′ Gaussian beam. Diagonal (white) noise is added with pixel variances and , comparable to the expected eventual level of Planck’s 143 GHz channel. For clarity of interpretation, no mask is used in these tests; in this setting, the eigenmodes of the fiducial covariance matrix are linear combinations of the spherical harmonics at fixed wavenumber ℓ. As fully exploring the likelihood is much more timeconsuming than profiling it, we generate a complete set of test results – analyses of the two test maps using cubic torus and slab covariance matrices on a fine grid of fundamental domain scales – using the profilelikelihood code, and aim to verify the main cubic torus results using the marginalized likelihood generated with AniCosmo. Note that to speed up the calculation of the marginalized likelihood we use a slightly smaller bandlimit (ℓ_{max} = 32) than in the profilelikelihood calculation (ℓ_{max} = 40); with our choice of smoothing scale and mode count, and considering the full sky for validation purposes, we obtain the same eigenbasis (and therefore analyse the same projected data) in both cases.
The results for the null test – in the form of the likelihood function for the fundamental domain scale of the assumed topology – are plotted in Fig. 7 for cubic tori and Fig. 8 for slabs. Concentrating initially on the cubic tori, we see that the likelihoods derived from the two codes agree. In both cases, the likelihood is found to be maximal for fundamental domain scales larger than the horizon, and the smallL cubic tori are very strongly disfavoured: . Note that the AniCosmo likelihood curve contains errors on the likelihood at each L considered, but these are orders of magnitude smaller than the changes in likelihood between points (typical MultiNest uncertainties yield errors of order 0.1 in loglikelihood).
In both cases, the profile likelihood exhibits a mild rise around the horizon scale, due to chance alignments along the matched faces of the fundamental domain. It is slightly stronger in the slab case since the probability of such alignments is greater with only a single pair of faces.
Fig. 7 Likelihood function for the fundamental domain scale of a cubic torus derived from simulations of a simplyconnected Universe, calculated through marginalization (black, filled circles) and profiling (pink, empty circles). The horizontal axis gives the inverse of the length of a side of the fundamental domain, relative to the distance to the lastscattering surface. The vertical lines mark the positions where χ_{rec} is equal to various characteristic sizes of the fundamental domain, namely the radius of the largest sphere that can be inscribed in the domain, ℛ_{i} = L/ 2, the smallest sphere in which the domain can be inscribed, , and the intermediate scale . 
Fig. 8 Profile likelihood function for the fundamental domain scale of a slab topology derived from simulations of a simplyconnected Universe. The vertical line marks the position where χ_{rec} is equal to the radius of the largest sphere that can be inscribed in the domain, ℛ_{i} = L/ 2 (for slab spaces, the other two characteristic sizes are infinite). 
The results for the tests on the toroidal simulation are shown in Fig. 9 for toroidal covariance matrices and Fig. 10 for slab covariance matrices. Concentrating first on the results employing toroidal covariance matrices, the correct fundamental domain scale is clearly picked out by both the profile and full likelihood codes, with the simply connected case strongly disfavoured at a likelihood ratio of . Turning to the results derived using slab covariance matrices, we see that – as expected from the KullbackLeibler divergence analysis of Sect. 3.1.2 – the correct fundamental domain scale is also found using the slab profile likelihood. Although, as also expected, the peak is not quite as pronounced when using the wrong covariance matrix, the simply connected Universe is still overwhelmingly disfavoured at a ratio of .
Fig. 9 Likelihood function for the fundamental domain scale of a cubic torus derived from a simulation of a toroidal Universe with or 2χ_{rec}/L = 1.5. The results from the profilelikelihood analysis (pink, clear circles) closely match those from the full marginalized likelihood (black, filled circles). Overlaid are additional profile likelihoods demonstrating the effects of changing the mode count and composition. In order of increasing constraining power, they utilize 837 IQU modes (grey, long dashed), 1085 I modes (blue, dotted), 1085 IQU modes (pink, clear circles), 2170 IQU modes (green, dotdashed), and finally 1085 noiseless IQU modes (purple, dashed). Adding lowℓ (ideally lownoise) polarization greatly increases the constraining power of the data. 
Fig. 10 Profile likelihood function for the fundamental domain scale of a slab derived from a simulation of a toroidal Universe with or 2χ_{rec}/L = 1.5. 
The speed of the profilelikelihood analysis allows for the effects of changing the mode count, composition, and noise level to be investigated. We repeat the toroidal test using intensityonly (I) and full (IQU) covariance matrices, retaining between 837 and 2170 modes at a time. For the smoothing scale employed in our tests, the 838th mode is the first to be dominated by polarization; runs using up to 837 IQU modes are therefore dominated by intensity information. The results of this investigation are contained in Fig. 9. The most striking conclusion is that the impact of adding temperature modes to the analysis is dwarfed by the impact of adding lowℓ polarization information, even though the temperature modes are effectively noiseless. This conclusion is supported by the observation from Fig. 1 that the KL divergence grows most rapidly at low ℓ.
4.4.2. Bianchi
The Bayesian analysis of Bianchi VII_{h} models using temperature data is performed using the AniCosmo code, which has been extensively validated by McEwen et al. (2013), and was used to perform the Bianchi analysis of Planck Collaboration XXVI (2014). The maximumlikelihood template fitting method used to analyse polarization data is straightforward and has been validated on simulations, correctly recovering the amplitude of templates artificially embedded in simulated CMB observations.
Fig. 11 (upper) and (lower) statistics as a function of circle radius α for the PlanckCommander (shortdashed red line), NILC (long dashed orange line), SEVEM (dotdashed green line), and SMICA (three dotdashed blue line) 2015 temperature maps. The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the SMICA CMB temperature map, smoothed and masked in the same way as the data, would yield a false event. The peak at 90° corresponds to a match between two copies of the same circle of radius 90° centred around two antipodal points. 
Fig. 12 (upper) and (lower) statistics as a function of circle radius α for the PlanckCommander (shortdashed red line), NILC (long dashed orange line), SEVEM (dotdashed green line), and SMICA (three dotdashed blue line) Emode maps. The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the SMICA CMB Emode map, smoothed and masked in the same way as the data, would yield a false event. The peak at 90° corresponds to a match between two copies of the same circle of radius 90° centred around two antipodal points. 
5. Results
5.1. Topology
5.1.1. Matched circles
We show the matchedcircle statistic for the CMB temperature and Emode maps in Figs. 11 and 12, respectively. We do not find any statistically significant correlation of circle pairs in any map. Results for the temperature maps are consistent with the Planck Collaboration XXVI (2014) results. As discussed in Sect. 4.4.1, the minimum radius at which the peaks expected for the matching statistic are larger than the false detection level for the polarization map is around α_{min} ≃ 15°. Thus, we can exclude, at the confidence level of 99%, any topology that predicts matching pairs of backtoback circles larger than this radius, assuming that the relative orientation of the fundamental domain and mask allows its detection. This implies that in a flat Universe described otherwise by the fiducial ΛCDM model, a 99% confidencelimit lower bound on the size of the fundamental domain is ℛ_{i} = L/ 2 ≳ χ_{rec}cos(α_{min}) = 0.97 χ_{rec} or . This is slightly stronger than the constraint obtained for the analysis of the 2013 Planck temperature maps, i.e., 0.94 χ_{rec} (Planck Collaboration XXVI 2014). Note that the limits from polarization are at least as strong as those from temperature despite the considerably smaller amount of sky considered in the polarization analysis (40% compared to 74%).
5.1.2. Likelihood
The results of applying the two likelihood codes to the Planck lowℓ data are plotted in Fig. 13 for cubic tori and Fig. 14 for slabs. As with the null test, smallL topologies are strongly ruled out, and the likelihoods are maximized at scales approaching or exceeding the horizon. In the marginalized case, we find that the maximumlikelihood fundamental domain scale of a cubic torus is , with scales of disfavoured at greater than 3.2 σ; for slabs, we find the likelihood to be peaked at (just inside the lastscattering surface), though is allowed at 1.9σ.
Fig. 13 Likelihood for the fundamental domain scale of a cubic torus derived from the Planck lowℓ data set. As in testing, the likelihoods calculated via profiling (pink, clear circles) and marginalization (black, solid circles) agree well. The impact of increasing the polarization content through additional smoothing (purple, dashed) or modes (green, dot dashed) is diminished compared to the test setting due to boosted noise. 
We have investigated the shape of the likelihood as a function of the slab orientation, and find that it is strongly peaked at an orientation such that the induced matched circles lie partially within the large polarization mask (retaining only 47% of the sky). These orientations therefore do not benefit from the extra discriminatory power of the polarization and its correlation with temperature.
As noted when analysing simulations of simplyconnected models (Figs. 7, 8), the profile likelihoods also show a mild rise around the horizon scale for the Planck data: this rise therefore cannot be interpreted as evidence for a multiply connected topology. We found a similar effect with both profile and marginalized likelihoods in Planck Collaboration XXVI (2014) using simulations of simply connected universes and the Planck temperature data, though this peak was considerably more pronounced for the profile likelihood. Further investigation is required to determine whether this extrahorizon rise is still present in the marginalized likelihood with the present data.
Fig. 14 Likelihood for the fundamental domain scale of a slab topology derived from the Planck lowℓ data set via profiling (pink, clear circles) and marginalization (black, solid circles). 
Though the smallL topologies are strongly constrained by the data, we note that the dropoff in likelihood is not as sharp as that observed in our null test (or, indeed, the signal test) presented in Sect. 4.4.1. There are several reasons for this behaviour. The default smoothing scale used in the lowℓPlanck data set (440′ Gaussian for intensity, no smoothing for polarization) means all of the first 1085 modes are in fact temperaturedominated. As our eigenbasis is constructed in decreasingeigenvalue order, the greater the smoothing scale, the more smallscale temperature modes are damped in comparison to largescale polarization modes, and thus the earlier polarization modes appear in the basis. Reducing the smoothing scale from 640′ (as used in testing) to 440′ (as used in the data) means the first polarizationdominated mode no longer appears in the 1085 highesteigenvalue modes we take as our fiducial basis. Exploiting once more the speed of the profile likelihood code, we therefore explore the dependence of the likelihoods on the number of modes and their composition in Fig. 13. We introduce polarizationdominated eigenmodes into the analysis in two ways: by simply doubling the number of modes and by retaining the mode count but applying additional smoothing (to bring the effective Gaussian smoothing FWHM to 640′). Though the constraints on smallL topologies do become tighter when polarizationdominated modes are included, their impact remains weaker than in testing. This is because the noise in the lowℓ data is significantly higher than in our tests – the typical diagonal covariance matrix element is – and, finally, because we must use a sky cut, which can hide excess correlations. This strongly motivates repeating this analysis with the full multifrequency Planck polarization data when they become available.
Natural logBayes factors of Bianchi models relative to equivalent ΛCDM model (positive favours Bianchi model).
Converting these likelihoods into Bayesian constraints on the size of the fundamental domain is not straightforward. Absent a proper prior giving an upper limit on the size of the fundamental domain, there is always infinite parameter volume available for everlarger fundamental domains. Hence, these likelihood plots should be considered the full summary of the 2015 Planck data for these models. Nonetheless, it is often useful to consider a falloff in the likelihood of Δlnℒ < −5 as roughly equivalent to a 3σ–99% confidence level–limit, the location of which we approximate by interpolating between calculated likelihood points. For the cubic T3 torus, we find that the marginal likelihoods of the combined temperature and polarization lowT,P data require that the length of an edge of the fundamental domain satisfies at this significance, or equivalently that the radius of the largest inscribed sphere in the fundamental domain is ℛ_{i}> 0.97 χ_{rec} (recall that is the comoving distance to the lastscattering surface.) The profile likelihood gives the somewhat weaker limit, ℛ_{i}> 0.79 χ_{rec}. For the T1 slab, we have or ℛ_{i}> 0.56 χ_{rec}. Because the temperature data on the relevant scales are still dominated by the cosmological signal, and the polarization noise remains large, these results are only slightly stronger than those presented in Planck Collaboration XXVI (2014).
5.2. Bianchi
Planck temperature data are masked and analysed for evidence of a Bianchi VII_{h} component, using the prior parameter ranges adopted in McEwen et al. (2013) and Planck Collaboration XXVI (2014). Cleaned temperature maps for each of the four component separation techniques are examined, where the mask defined for each technique is applied. The natural logBayes factors for the Bianchi models relative to their standard cosmological counterparts are shown in Table 1. The Bayes factors are broadly consistent across the component separation methods. Most Bayes factors are similar to the values computed in Planck Collaboration XXVI (2014), with the exception of the analysis of the lefthanded flatdecoupledBianchi model with SEVEM data, which has increased, but which is now more consistent with the other component separation methods.
Fig. 15 Marginalized posterior distributions of Bianchi parameters recovered from PlanckSMICA (solid blue curves), SEVEM (dashed green curves), NILC (dotdashed yellow curves), and Commander (dotted red curves) componentseparated data for lefthanded models. Planck data provide evidence in support of a Bianchi component in the phenomenological flatdecoupledBianchi model (panel a)) but not in the physical opencoupledBianchi model (panel b)). Significant differences exist between the posterior distributions shown in panel b) for each component separation method; this model is not favoured by data and parameters are in general poorly constrained. 
Parameters recovered for the lefthanded flatdecoupledBianchi model.
For the phenomenological flatdecoupledBianchi model, evidence in support of a lefthanded Bianchi template is again found (Table 1). The Bayes factors providing evidence for this model range between the values 2.3 ± 0.2 and 3.2 ± 0.2, corresponding to odds ratios of approximately 1:10 and 1:25, respectively (which on the Jeffreys scale are categorized as significant and strong, respectively; Jeffreys 1961). Recovered posterior distributions of the Bianchi parameters of this model for each component separation technique are shown in Fig. 15a. The posterior distributions are consistent across component separation techniques, are similar to those recovered in Planck Collaboration XXVI (2014), and are reasonably well constrained (with the exception of the known – Bianchi parameter degeneracy Jaffe et al. 2006c; Bridges et al. 2007). Recall that the Bianchi densities ( and ) are decoupled from the standard cosmology in the flatdecoupledBianchi model considered here and, as found previously, are inconsistent with standard estimates of the densities. The maximumaposteriori (MAP) and meanposterior Bianchi parameter estimates for this model are given in Table 2, while the corresponding MAP bestfit Bianchi temperature maps are shown in Fig. 16. Note that the maximum of a marginalized onedimensional posterior (e.g., Fig. 15) will not in general coincide with the global MAP estimate for the full set of parameters. The bestfit maps for the lefthanded flatcoupledBianchi model are consistent across component separation techniques and similar to the bestfit maps found in previous Planck (Planck Collaboration XXVI 2014) and WMAP (see, e.g., Jaffe et al. 2005; and McEwen et al. 2013) temperature data.
For the physical opencoupledBianchi model where the Bianchi VII_{h} model is coupled to the standard cosmology, there is again no evidence in support of a Bianchi contribution (Table 1). Recovered posterior distributions of the parameters of this model for each component separation technique are shown in Fig. 15b. The posterior distributions show some similarity across component separation techniques and with the distributions recovered in Planck Collaboration XXVI (2014). However, significant differences exist since the parameters are in general poorly constrained and the model is not favoured by the Bayesian evidence. The MAP and meanposterior Bianchi parameter estimates for this model are given in Table 3, while the corresponding MAP bestfit Bianchi temperature maps are shown in Fig. 17. While the posterior distributions show differences between component separation techniques, the estimated parameters are consistent. Note that, although the mean parameter estimates have not changed markedly from Planck Collaboration XXVI (2014), the MAP estimates have indeed changed. Intriguingly, for each component separation technique, the MAP bestfit maps for the opencoupledBianchi model (Fig. 17), for which there is no evidence in support of a Bianchi contribution, show a similar but not identical morphology to the MAP bestfit maps for the flatdecoupledBianchi model (Fig. 16), for which the Bayesian evidence supports the inclusion of a Bianchi component. This was not the case in the previous Planck analysis (Planck Collaboration XXVI 2014, Fig. 22). A parameter combination is found for the physical opencoupledBianchi model that is broadly consistent with a standard open cosmology and that produces a Bianchi temperature map similar to the one found in the unphysical flatdecoupledBianchi model. This parameter combination lies on the known – Bianchi parameter degeneracy (Jaffe et al. 2006c; Bridges et al. 2007) and lies close to but not directly on the wellknown CMB geometric degeneracy (determined by an independent CMB analysis), since both the cosmological and Bianchi components are fitted simultaneously. It is important to stress that Planck temperature data do not favour the physical opencoupledBianchi model, but neither is it possible to rule out this model using temperature data alone (Planck polarization data are considered subsequently). An overall constraint on the vorticity of Bianchi VII_{h} models, from Planck temperature data alone, of (ω/H)_{0}< 7.6 × 10^{10} (95% confidence level) is obtained from the analysis of the physical opencoupledBianchi model, which is consistent across all component separation techniques, as illustrated in Table 4.
To further constrain Bianchi VII_{h} models using Planck polarization data we simulate Bianchi polarization maps, computed using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.) for the bestfit Bianchi parameters determined from the temperature analysis. E and Bmode Bianchi maps for the bestfit flatdecoupledBianchi model and the opencoupledBianchi model are displayed in Figs. 18 and 19, respectively. As described in Sect. 4.3.3, we estimate the maximumlikelihood amplitude of these Bianchi polarization templates in Planck polarization data, performing the analysis in the pixel space defined by the Q and U Stokes components, where noise and partial sky coverage can be handled effectively. The maximumlikelihood amplitudes estimated for each component separation technique are given in Table 5. The estimated amplitudes are close to zero and consistent across component separation techniques. The difference between the estimated amplitudes and zero is more likely due to small residual foreground contamination than a Bianchi component. Indeed, the amplitude estimates are at least 24 standard deviations from unity, the expected value for the bestfit Bianchi models determined from the temperature analysis. Both the bestfit flatdecoupledBianchi model and the opendecoupledBianchi model are thus strongly disfavoured by the Planck polarization data. This is not surprising since these models produce relatively strong E and Bmode contributions (see Figs. 18 and 19), as highlighted already by Pontzen & Challinor (2007). However, the full freedom of Bianchi models remains to be explored using temperature and polarization data simultaneously, for example through a complete Bayesian analysis, which is left to future work.
Parameters recovered for the lefthanded opencoupledBianchi model.
Upper bounds on vorticity (ω/H)_{0} at 95% confidence level.
Maximumlikelihood amplitude estimates λ^{ML} and 1 σ errors computed using polarization data.
Fig. 16 Bestfit temperature maps for the lefthanded flatdecoupledBianchi model. 
Fig. 17 Bestfit temperature maps for the lefthanded opencoupledBianchi model. 
Fig. 18 Polarization maps for the bestfit lefthanded flatdecoupledBianchi model fitted to temperature data. In each panel E (left) and Bmode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.). 
Fig. 19 Polarization maps for the bestfit lefthanded opencoupledBianchi model fitted to temperature data. In each panel E (left) and Bmode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.). 
6. Discussion
We have used Planck intensity and polarization data to evaluate specific departures from the largescale isotropy of the Universe. Using both frequentist and Bayesian methods applied for the first time to polarization data, we find no evidence for a multiconnected topology with a scale less than roughly the distance to the lastscattering surface. Specifically, a frequentist search for antipodal matched circles on N_{side} = 512 maps finds a lower bound on the size of the fundamental domain of 0.97 χ_{rec} from polarization data alone. Using Bayesian methods applied to lowresolution (N_{side} = 16) maps of both temperature and polarization, we also find a lower limit of 0.97 χ_{rec} for the T3 cubic torus (for the T1 slab, the limit is 0.56 χ_{rec}). These results are both consistent and complementary, giving coincidentally identical limits but with very different statistical foundations and data selections. The addition of polarization data at current levels of accuracy does not significantly improve the limits from intensity alone, but we have found that the polarization sensitivity of the full set of Planck detectors should give quantitative improvements in the limits, decreasing the likelihood of fundamental domains with scales smaller than the distance to the lastscattering surface by many orders of magnitude.
We also find no evidence for a Bianchi VII_{h} model which departs from global isotropy via the presence of both rotation and shear. Although the largescale temperature pattern measured by Planck has some similar features to that induced by focusing in a Bianchi VII_{h} Universe, it requires unphysical parameters. Fixed to those parameters, we have further shown that the polarization pattern induced by such models is strongly disfavoured by the Planck data.
The results outlined here show no evidence for departures from isotropic and simplyconnected models. Improved computational techniques, along with future polarization data from Planck (and beyond), will allow yet stronger checks of even wider classes of models. For the multiconnected case, we can expand to models without antipodal matched circles and with the scale of the fundamental domain closer to (and even slightly beyond) the lastscattering surface. For the likelihood method, this will require computation of the correlation matrix with higherwavenumber modes to capture more of the available information.
For anisotropic models, we can explicitly perform parameter estimation and model selection using polarization data, beyond the simple templatefitting performed here.
Although the evidence thus far corroborates the conventional wisdom that we live in the simplest FLRW Universe, this is likely to be only an approximation vastly beyond the Hubble scale. Detection of a multiplyconnected topology or anisotropic geometry is one of the few ways to probe the global structure of spacetime. We have shown that Planck data give the best handle to date on these possibilities.
Planck (http://www.esa.int/Planck) is a project of the European Space Agency (ESA) with instruments provided by two scientific consortia funded by ESA member states and led by Principal Investigators from France and Italy, telescope reflectors provided through a collaboration between ESA and a scientific consortium led and funded by Denmark, and additional contributions from NASA (USA).
Acknowledgments
The Planck Collaboration acknowledges the support of: ESA; CNES and CNRS/INSUIN2P3INP (France); ASI, CNR, and INAF (Italy); NASA and DoE (USA); STFC and UKSA (UK); CSIC, MINECO, JA, and RES (Spain); Tekes, AoF, and CSC (Finland); DLR and MPG (Germany); CSA (Canada); DTU Space (Denmark); SER/SSO (Switzerland); RCN (Norway); SFI (Ireland); FCT/MCTES (Portugal); ERC and PRACE (EU). A description of the Planck Collaboration and a list of its members, indicating which technical or scientific activities they have been involved in, can be found at http://www.cosmos.esa.int/web/planck/planckcollaboration. The authors thank Andrew Pontzen for computing Bianchi VII_{h} polarization templates for the bestfit models resulting from the analysis of temperature data. We acknowledge the UCL Legion High Performance Computing Facility (Legion@UCL) and associated support services in the completion of this work. Parts of the computations were performed on the Andromeda and Perseus clusters of the University of Geneva, as well as the Carver IBM iDataPlex, the Hopper Cray XE6, and the Edison Cray XC30 at NERSC, and on the GPC supercomputer at the SciNet HPC Consortium. SciNet is funded by: the Canada Foundation for Innovation under the auspices of Compute Canada; the Government of Ontario; Ontario Research Fund − Research Excellence; and the University of Toronto.
References
 Aurich, R. 1999, ApJ, 524, 497 [NASA ADS] [CrossRef] [Google Scholar]
 Aurich, R., & Lustig, S. 2013, MNRAS, 433, 2517 [NASA ADS] [CrossRef] [Google Scholar]
 Aurich, R., & Lustig, S. 2014, Class. Quant. Grav., 31, 165009 [NASA ADS] [CrossRef] [Google Scholar]
 Aurich, R., Lustig, S., Steiner, F., & Then, H. 2004, Class. Quant. Grav., 21, 4901 [NASA ADS] [CrossRef] [Google Scholar]
 Aurich, R., Lustig, S., & Steiner, F. 2005, Class. Quant. Grav., 22, 2061 [Google Scholar]
 Aurich, R., Lustig, S., & Steiner, F. 2006, MNRAS, 369, 240 [NASA ADS] [CrossRef] [Google Scholar]
 Aurich, R., Janzer, H. S., Lustig, S., & Steiner, F. 2008, Class. Quant. Grav., 25, 125006 [NASA ADS] [CrossRef] [Google Scholar]
 Barrow, J. D. 1986, Can. J. Phys., 64, 152 [NASA ADS] [CrossRef] [Google Scholar]
 Barrow, J. D., Juszkiewicz, R., & Sonoda, D. H. 1985, MNRAS, 213, 917 [NASA ADS] [CrossRef] [Google Scholar]
 Bennett, C. L., Banday, A. J., Gorski, K. M., et al. 1996, ApJ, 464, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Bielewicz, P., & Banday, A. J. 2011, MNRAS, 412, 2104 [NASA ADS] [CrossRef] [Google Scholar]
 Bielewicz, P., & Riazuelo, A. 2009, MNRAS, 396, 609 [NASA ADS] [CrossRef] [Google Scholar]
 Bielewicz, P., Banday, A. J., & Górski, K. M. 2012, MNRAS, 421, 1064 [NASA ADS] [CrossRef] [Google Scholar]
 Bond, J. R., & Efstathiou, G. 1987, MNRAS, 226, 655 [NASA ADS] [Google Scholar]
 Bond, J. R., Jaffe, A. H., & Knox, L. E. 1998a, Phys. Rev. D, 57, 2117 [NASA ADS] [CrossRef] [Google Scholar]
 Bond, J. R., Pogosyan, D., & Souradeep, T. 1998b, Class. Quant. Grav., 15, 2671 [NASA ADS] [CrossRef] [Google Scholar]
 Bond, J. R., Pogosyan, D., & Souradeep, T. 2000, Phys. Rev. D, 62, 043006 [NASA ADS] [CrossRef] [Google Scholar]
 Bridges, M., McEwen, J. D., Lasenby, A. N., & Hobson, M. P. 2007, MNRAS, 377, 1473 [NASA ADS] [CrossRef] [MathSciNet] [Google Scholar]
 Bridges, M., McEwen, J. D., Cruz, M., et al. 2008, MNRAS, 390, 1372 [NASA ADS] [Google Scholar]
 Bunn, E. F., Ferreira, P. G., & Silk, J. 1996, Phys. Lett., 77, 2883 [Google Scholar]
 Caillerie, S., LachièzeRey, M., Luminet, J.P., et al. 2007, A&A, 476, 691 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Cayón, L., Banday, A. J., Jaffe, T., et al. 2006, MNRAS, 369, 598 [NASA ADS] [CrossRef] [Google Scholar]
 Collins, C. B., & Hawking, S. W. 1973, MNRAS, 162, 307 [NASA ADS] [CrossRef] [Google Scholar]
 Cornish, N. J., Spergel, D. N., & Starkman, G. D. 1998, Class. Quant. Grav., 15, 2657 [Google Scholar]
 Cornish, N. J., Spergel, D. N., Starkman, G. D., & Komatsu, E. 2004, Phys. Rev. Lett., 92, 201302 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 de OliveiraCosta, A., & Smoot, G. F. 1995, ApJ, 448, 477 [NASA ADS] [CrossRef] [Google Scholar]
 Dineen, P., Rocha, G., & Coles, P. 2005, MNRAS, 358, 1285 [NASA ADS] [CrossRef] [Google Scholar]
 Fagundes, H. V., & Wichoski, U. F. 1987, Nature, 322, L5 [Google Scholar]
 Fang, L.Z., & Sato, H. 1983, Commun. Theor. Phys., 2, 1055 [Google Scholar]
 Feroz, F., & Hobson, M. P. 2008, MNRAS, 384, 449 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., Hobson, M. P., & Bridges, M. 2009, MNRAS, 398, 1601 [NASA ADS] [CrossRef] [Google Scholar]
 Feroz, F., Hobson, M. P., Cameron, E., & Pettitt, A. N. 2013, ArXiv eprints [arXiv:1306.2144] [Google Scholar]
 Fujii, H., & Yoshii, Y. 2011, A&A, 529, A121 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Ghosh, T., Hajian, A., & Souradeep, T. 2007, Phys. Rev. D, 75, 083007 [Google Scholar]
 Jaffe, T. R., Banday, A. J., Eriksen, H. K., Górski, K. M., & Hansen, F. K. 2005, ApJ, 629, L1 [NASA ADS] [CrossRef] [Google Scholar]
 Jaffe, T. R., Banday, A. J., Eriksen, H. K., Górski, K. M., & Hansen, F. K. 2006a, A&A, 460, 393 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Jaffe, T. R., Banday, A. J., Eriksen, H. K., Górski, K. M., & Hansen, F. K. 2006b, ApJ, 643, 616 [NASA ADS] [CrossRef] [Google Scholar]
 Jaffe, T. R., Hervik, S., Banday, A. J., & Górski, K. M. 2006c, ApJ, 644, 701 [NASA ADS] [CrossRef] [Google Scholar]
 Jarosik, N., Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 [NASA ADS] [CrossRef] [Google Scholar]
 Jeffreys, H. 1961, Theory of probability, 3rd edn. (Oxford: Oxford University Press) [Google Scholar]
 Key, J. S., Cornish, N. J., Spergel, D. N., & Starkman, G. D. 2007, Phys. Rev. D, 75, 084034 [NASA ADS] [CrossRef] [Google Scholar]
 Kim, J. 2011, A&A, 531, A32 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Kogut, A., Hinshaw, G., & Banday, A. J. 1997, Phys. Rev. D, 55, 1901 [NASA ADS] [CrossRef] [Google Scholar]
 Kosowsky, A. 1996, Ann. Phys., 246, 49 [NASA ADS] [CrossRef] [Google Scholar]
 Kunz, M., Aghanim, N., Cayon, L., et al. 2006, Phys. Rev. D, D73, 023511 [NASA ADS] [CrossRef] [Google Scholar]
 Kunz, M., Aghanim, N., Riazuelo, A., & Forni, O. 2008, Phys. Rev. D, 77, 23525 [Google Scholar]
 LachiezeRey, M., & Luminet, J. 1995, Phys. Rep., 254, 135 [NASA ADS] [CrossRef] [Google Scholar]
 Land, K., & Magueijo, J. 2006, MNRAS, 367, 1714 [NASA ADS] [CrossRef] [Google Scholar]
 Lehoucq, R., LachiezeRey, M., & Luminet, J. P. 1996, A&A, 313, 339 [NASA ADS] [Google Scholar]
 Levin, J. 2002, Phys. Rep., 365, 251 [NASA ADS] [CrossRef] [Google Scholar]
 Levin, J., Scannapieco, E., & Silk, J. 1998, Nature, 58, 103516 [Google Scholar]
 Lew, B., & Roukema, B. 2008, A&A, 482, 747 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Luminet, J.P., Weeks, J. R., Riazuelo, A., Lehoucq, R., & Uzan, J.P. 2003, Nature, 425, 593 [NASA ADS] [CrossRef] [PubMed] [Google Scholar]
 McEwen, J. D., Hobson, M. P., Lasenby, A. N., & Mortlock, D. J. 2006, MNRAS, 369, 1858 [NASA ADS] [CrossRef] [Google Scholar]
 McEwen, J. D., Josset, T., Feeney, S. M., Peiris, H. V., & Lasenby, A. N. 2013, MNRAS, 436, 3680 [NASA ADS] [CrossRef] [Google Scholar]
 Niarchou, A., & Jaffe, A. H. 2007, Phys. Rev. Lett., 99, 81302 [Google Scholar]
 Niarchou, A., Jaffe, A. H., & Pogosian, L. 2004, Phys. Rev. D, 69, 063515 [Google Scholar]
 Phillips, N. G., & Kogut, A. 2006, ApJ, 645, 820 [NASA ADS] [CrossRef] [Google Scholar]
 Planck Collaboration XXVI. 2014, A&A, 571, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration I. 2016, A&A, 594, A1 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration II. 2016, A&A, 594, A2 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration III. 2016, A&A, 594, A3 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IV. 2016, A&A, 594, A4 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration V. 2016, A&A, 594, A5 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VI. 2016, A&A, 594, A6 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VII. 2016, A&A, 594, A7 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration VIII. 2016, A&A, 594, A8 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration IX. 2016, A&A, 594, A9 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration X. 2016, A&A, 594, A10 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XI. 2016, A&A, 594, A11 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XII. 2016, A&A, 594, A12 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIII. 2016, A&A, 594, A13 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIV. 2016, A&A, 594, A14 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XV. 2016, A&A, 594, A15 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVI. 2016, A&A, 594, A16 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVII. 2016, A&A, 594, A17 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XVIII. 2016, A&A, 594, A18 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XIX. 2016, A&A, 594, A19 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XX. 2016, A&A, 594, A20 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXI. 2016, A&A, 594, A21 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXII. 2016, A&A, 594, A22 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIII. 2016, A&A, 594, A23 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXIV. 2016, A&A, 594, A24 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXV. 2016, A&A, 594, A25 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVI. 2016, A&A, 594, A26 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVII. 2016, A&A, 594, A27 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Planck Collaboration XXVIII. 2016, A&A, 594, A28 [NASA ADS] [CrossRef] [EDP Sciences] [Google Scholar]
 Pontzen, A. 2009, Phys. Rev. D, 79, 103518 [NASA ADS] [CrossRef] [Google Scholar]
 Pontzen, A., & Challinor, A. 2007, MNRAS, 380, 1387 [NASA ADS] [CrossRef] [Google Scholar]
 Pontzen, A., & Challinor, A. 2011, Class. Quant. Grav., 28, 185007 [NASA ADS] [CrossRef] [Google Scholar]
 Press, W. H., Teukolsky, S. A., Vetterling, W. T., & Flannery, B. P. 1992, Numerical Recipes in C: The Art of Scientific Computing, 2nd edn. (New York, NY, USA: Cambridge University Press) [Google Scholar]
 Riazuelo, A., Uzan, J.P., Lehoucq, R., & Weeks, J. 2004a, Phys. Rev. D, 69, 103514 [NASA ADS] [CrossRef] [Google Scholar]
 Riazuelo, A., Weeks, J., Uzan, J.P., Lehoucq, R., & Luminet, J.P. 2004b, Phys. Rev. D, 69, 103518 [NASA ADS] [CrossRef] [Google Scholar]
 Riazuelo, A., Caillerie, S., LachièzeRey, M., Lehoucq, R., & Luminet, J.P. 2006, ArXiv eprints [arXiv:astroph/0601433] [Google Scholar]
 Rocha, G., Cayón, L., Bowen, R., et al. 2004, MNRAS, 351, 769 [NASA ADS] [CrossRef] [Google Scholar]
 Roukema, B. F. 1996, MNRAS, 283, 1147 [NASA ADS] [Google Scholar]
 Roukema, B. F. 2000a, Class. Quant. Grav., 17, 3951 [NASA ADS] [CrossRef] [Google Scholar]
 Roukema, B. F. 2000b, MNRAS, 312, 712 [NASA ADS] [CrossRef] [Google Scholar]
 Roukema, B. F., Buliński, Z., & Gaudin, N. E. 2008, A&A, 492, 657 [Google Scholar]
 Seljak, U., & Zaldarriaga, M. 1996, ApJ, 469, 437 [NASA ADS] [CrossRef] [Google Scholar]
 Skilling, J. 2004, in AIP Conf. Ser. 735, eds. R. Fischer, R. Preuss, & U. V. Toussaint, 395 [Google Scholar]
 Sokolov, I. Y. 1993, Sov. J. Exp. Theor. Phys. Lett., 57, 617 [Google Scholar]
 Sokolov, D. D., & Shvartsman, V. F. 1974, Sov. J. Exp. Theor. Phys., 39, 196 [Google Scholar]
 Starobinskij, A. A. 1993, Sov. J. Exp. Theor. Phys. Lett., 57, 622 [Google Scholar]
 Stevens, D., Scott, D., & Silk, J. 1993, Phys. Rev. Lett., 71, 20 [NASA ADS] [CrossRef] [Google Scholar]
 Vaudrevange, P. M., Starkman, G. D., Cornish, N. J., & Spergel, D. N. 2012, Phys. Rev. D, 86, 083526 [NASA ADS] [CrossRef] [Google Scholar]
 Weatherley, S. J., Warren, S. J., Croom, S. M., et al. 2003, Nature, 342, L9 [Google Scholar]
All Tables
Natural logBayes factors of Bianchi models relative to equivalent ΛCDM model (positive favours Bianchi model).
Maximumlikelihood amplitude estimates λ^{ML} and 1 σ errors computed using polarization data.
All Figures
Fig. 1 KL divergence of slab (T1) topologies relative to an infinite Universe as a function of ℓ_{max} with sizes and (top), and as a function of size L of the fundamental domain for various ℓ_{max} (bottom). A torus with , corresponding to (H_{0}L)^{1}< 0.154, has a fundamental domain that is larger than the distance to the lastscattering surface and leaves only a small trace in the CMB. This is why the KL divergence drops rapidly at this point. Note that the information for continues to rise with ℓ_{max} whereas it levels off for the slightly larger case. In the lower panel we see that there is a slight feature in d_{KL} at about half the horizon distance, which is probably due to harmonic effects. The corresponding figures for cubic (T3) topologies look qualitatively similar except that all d_{KL} values are three times larger. 

In the text 
Fig. 2 KL divergence of fixed cubic, rod, and slab topologies with fundamental domain side compared to a slab of variable fundamental domain size L. The chimney space T2 dates from the 2013 analysis (Planck Collaboration XXVI 2014) and was computed for the bestfit parameters of that release. In all cases the smallest KL divergence, corresponding to the best fit, appears at , indicating that the slab space can be used to detect other topologies. An additional dip at L ≃ 5.5 / (2H_{0}) may be due to a harmonic effect at half the size of the fundamental domain; it is, however, much smaller than the drop in KL divergence at the size of the fundamental domain. 

In the text 
Fig. 3 KL divergence of a slab space relative to a cubic topology, as a function of rotation angle of the slab space (blue curve). Both spaces have and ℓ_{max} = 80. The horizontal black dashed line gives the KL divergence of an infinite Universe relative to the cubic topology and illustrates how much better the slab space fits with the correct orientation relative to the cubic torus. 

In the text 
Fig. 4 Mode structure plotted as maps for the eigenvector corresponding to the highestsignal eigenvalue of the fiducial simplyconnected model. The top map corresponds to temperature, middle to Q polarization, and bottom to U polarization. Masked pixels are plotted in grey. 

In the text 
Fig. 5 Mode structure plotted as maps for the eigenvector corresponding to the 301sthighestsignal eigenvalue of the fiducial simplyconnected model. The top map corresponds to temperature, middle to Q polarization, and bottom to U polarization. Masked pixels are plotted in grey. 

In the text 
Fig. 6 Example statistic as a function of circle radius α for a simulated CMB Emode map of a Universe with the topology of a cubic 3torus with dimensions . To this map noise was added corresponding to the SMICA map. The thick overlapping curves show the statistic for simulated polarization maps with angular resolution and noise level corresponding to the PlanckSMICA map for three cases: without highpass filtering (solid, red); with filtering (dashed, blue); and with a 5% larger amplitude of noise (dotdashed, green). The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the highpass filtered PlanckSMICA polarization maps, smoothed and masked in the same way as the data, would yield a false event. 

In the text 
Fig. 7 Likelihood function for the fundamental domain scale of a cubic torus derived from simulations of a simplyconnected Universe, calculated through marginalization (black, filled circles) and profiling (pink, empty circles). The horizontal axis gives the inverse of the length of a side of the fundamental domain, relative to the distance to the lastscattering surface. The vertical lines mark the positions where χ_{rec} is equal to various characteristic sizes of the fundamental domain, namely the radius of the largest sphere that can be inscribed in the domain, ℛ_{i} = L/ 2, the smallest sphere in which the domain can be inscribed, , and the intermediate scale . 

In the text 
Fig. 8 Profile likelihood function for the fundamental domain scale of a slab topology derived from simulations of a simplyconnected Universe. The vertical line marks the position where χ_{rec} is equal to the radius of the largest sphere that can be inscribed in the domain, ℛ_{i} = L/ 2 (for slab spaces, the other two characteristic sizes are infinite). 

In the text 
Fig. 9 Likelihood function for the fundamental domain scale of a cubic torus derived from a simulation of a toroidal Universe with or 2χ_{rec}/L = 1.5. The results from the profilelikelihood analysis (pink, clear circles) closely match those from the full marginalized likelihood (black, filled circles). Overlaid are additional profile likelihoods demonstrating the effects of changing the mode count and composition. In order of increasing constraining power, they utilize 837 IQU modes (grey, long dashed), 1085 I modes (blue, dotted), 1085 IQU modes (pink, clear circles), 2170 IQU modes (green, dotdashed), and finally 1085 noiseless IQU modes (purple, dashed). Adding lowℓ (ideally lownoise) polarization greatly increases the constraining power of the data. 

In the text 
Fig. 10 Profile likelihood function for the fundamental domain scale of a slab derived from a simulation of a toroidal Universe with or 2χ_{rec}/L = 1.5. 

In the text 
Fig. 11 (upper) and (lower) statistics as a function of circle radius α for the PlanckCommander (shortdashed red line), NILC (long dashed orange line), SEVEM (dotdashed green line), and SMICA (three dotdashed blue line) 2015 temperature maps. The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the SMICA CMB temperature map, smoothed and masked in the same way as the data, would yield a false event. The peak at 90° corresponds to a match between two copies of the same circle of radius 90° centred around two antipodal points. 

In the text 
Fig. 12 (upper) and (lower) statistics as a function of circle radius α for the PlanckCommander (shortdashed red line), NILC (long dashed orange line), SEVEM (dotdashed green line), and SMICA (three dotdashed blue line) Emode maps. The dotted line shows the false detection level established such that fewer than 1% of 300 Monte Carlo simulations of the SMICA CMB Emode map, smoothed and masked in the same way as the data, would yield a false event. The peak at 90° corresponds to a match between two copies of the same circle of radius 90° centred around two antipodal points. 

In the text 
Fig. 13 Likelihood for the fundamental domain scale of a cubic torus derived from the Planck lowℓ data set. As in testing, the likelihoods calculated via profiling (pink, clear circles) and marginalization (black, solid circles) agree well. The impact of increasing the polarization content through additional smoothing (purple, dashed) or modes (green, dot dashed) is diminished compared to the test setting due to boosted noise. 

In the text 
Fig. 14 Likelihood for the fundamental domain scale of a slab topology derived from the Planck lowℓ data set via profiling (pink, clear circles) and marginalization (black, solid circles). 

In the text 
Fig. 15 Marginalized posterior distributions of Bianchi parameters recovered from PlanckSMICA (solid blue curves), SEVEM (dashed green curves), NILC (dotdashed yellow curves), and Commander (dotted red curves) componentseparated data for lefthanded models. Planck data provide evidence in support of a Bianchi component in the phenomenological flatdecoupledBianchi model (panel a)) but not in the physical opencoupledBianchi model (panel b)). Significant differences exist between the posterior distributions shown in panel b) for each component separation method; this model is not favoured by data and parameters are in general poorly constrained. 

In the text 
Fig. 16 Bestfit temperature maps for the lefthanded flatdecoupledBianchi model. 

In the text 
Fig. 17 Bestfit temperature maps for the lefthanded opencoupledBianchi model. 

In the text 
Fig. 18 Polarization maps for the bestfit lefthanded flatdecoupledBianchi model fitted to temperature data. In each panel E (left) and Bmode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.). 

In the text 
Fig. 19 Polarization maps for the bestfit lefthanded opencoupledBianchi model fitted to temperature data. In each panel E (left) and Bmode (right) maps are shown. These polarization maps are simulated using the approach of Pontzen & Challinor (2007) and Pontzen (2009), and provided by Pontzen (priv. comm.). 

In the text 
Current usage metrics show cumulative count of Article Views (fulltext article views including HTML views, PDF and ePub downloads, according to the available data) and Abstracts Views on Vision4Press platform.
Data correspond to usage on the plateform after 2015. The current usage metrics is available 4896 hours after online publication and is updated daily on week days.
Initial download of the metrics may take a while.