This paper was converted on www.awesomepapers.org from LaTeX by an anonymous user.
Want to know more? Visit the Converter page.

\newunicodechar

ʻ

The DES and SPT Collaborations

Joint analysis of DES Year 3 data and CMB lensing from SPT and Planck III:
Combined cosmological constraints

T. M. C. Abbott Cerro Tololo Inter-American Observatory, NSF’s National Optical-Infrared Astronomy Research Laboratory, Casilla 603, La Serena, Chile    M. Aguena Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    A. Alarcon Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA    O. Alves Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    A. Amon Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK    F. Andrade-Oliveira Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA    J. Annis Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    B. Ansarinejad School of Physics, University of Melbourne, Parkville, VIC 3010, Australia    S. Avila Instituto de Fisica Teorica UAM/CSIC, Universidad Autonoma de Madrid, 28049 Madrid, Spain    D. Bacon Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK    E. J. Baxter Institute for Astronomy, University of Hawai’i, 2680 Woodlawn Drive, Honolulu, HI 96822, USA    K. Bechtol Physics Department, 2320 Chamberlin Hall, University of Wisconsin-Madison, 1150 University Avenue Madison, WI 53706-1390    M. R. Becker Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA    B. A. Benson Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    G. M. Bernstein Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    E. Bertin CNRS, UMR 7095, Institut d’Astrophysique de Paris, F-75014, Paris, France Sorbonne Universités, UPMC Univ Paris 06, UMR 7095, Institut d’Astrophysique de Paris, F-75014, Paris, France    J. Blazek Department of Physics, Northeastern University, Boston, MA 02115, USA    L. E. Bleem High-Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue., Argonne, IL, 60439, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    S. Bocquet University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany    D. Brooks Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK    E. Buckley-Geer Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    D. L. Burke Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    H. Camacho Instituto de Física Teórica, Universidade Estadual Paulista, São Paulo, Brazil Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    A. Campos Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA    J. E. Carlstrom Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA High-Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue., Argonne, IL, 60439, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA    A. Carnero Rosell Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Universidad de La Laguna, Dpto. Astrofísica, E-38206 La Laguna, Tenerife, Spain    M. Carrasco Kind Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 W. Green Street, Urbana, IL 61801, USA    J. Carretero Institut de Física d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain    R. Cawthon Physics Department, William Jewell College, Liberty, MO, 64068    C. Chang Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    C. L. Chang High-Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue., Argonne, IL, 60439, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    R. Chen Department of Physics, Duke University Durham, NC 27708, USA    A. Choi California Institute of Technology, 1200 East California Blvd, MC 249-17, Pasadena, CA 91125, USA    R. Chown Department of Physics & Astronomy, The University of Western Ontario, London ON N6A 3K7, Canada Institute for Earth and Space Exploration, The University of Western Ontario, London ON N6A 3K7, Canada    C. Conselice Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK University of Nottingham, School of Physics and Astronomy, Nottingham NG7 2RD, UK    J. Cordero Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK    M. Costanzi Astronomy Unit, Department of Physics, University of Trieste, via Tiepolo 11, I-34131 Trieste, Italy INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, I-34143 Trieste, Italy Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy    T. Crawford Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    A. T. Crites Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA California Institute of Technology, 1200 East California Boulevard., Pasadena, CA, 91125, USA    M. Crocce Institut d’Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain    L. N. da Costa Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    C. Davis Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA    T. M. Davis School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia    T. de Haan High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan Department of Physics, University of California, Berkeley, CA, 94720, USA    J. De Vicente Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain    J. DeRose Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA    S. Desai Department of Physics, IIT Hyderabad, Kandi, Telangana 502285, India    H. T. Diehl Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    M. A. Dobbs Department of Physics and McGill Space Institute, McGill University, 3600 Rue University, Montreal, Quebec H3A 2T8, Canada Canadian Institute for Advanced Research, CIFAR Program in Gravity and the Extreme Universe, Toronto, ON, M5G 1Z8, Canada    S. Dodelson Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA NSF AI Planning Institute for Physics of the Future, Carnegie Mellon University, Pittsburgh, PA 15213, USA    P. Doel Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK    C. Doux Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA Université Grenoble Alpes, CNRS, LPSC-IN2P3, 38000 Grenoble, France    A. Drlica-Wagner Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    K. Eckert Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    T. F. Eifler Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA    F. Elsner Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK    J. Elvin-Poole Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Physics, The Ohio State University, Columbus, OH 43210, USA    S. Everett Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA    W. Everett Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO, 80309, USA    X. Fang Department of Astronomy, University of California, Berkeley, 501 Campbell Hall, Berkeley, CA 94720, USA Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA    I. Ferrero Institute of Theoretical Astrophysics, University of Oslo. P.O. Box 1029 Blindern, NO-0315 Oslo, Norway    A. Ferté Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA    B. Flaugher Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    P. Fosalba Institut d’Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain    O. Friedrich Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK    J. Frieman Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    J. García-Bellido Instituto de Fisica Teorica UAM/CSIC, Universidad Autonoma de Madrid, 28049 Madrid, Spain    M. Gatti Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    E. M. George European Southern Observatory, Karl-Schwarzschild-Straße 2, 85748 Garching, Germany    T. Giannantonio Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK    G. Giannini Institut de Física d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain    D. Gruen University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany    R. A. Gruendl Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 W. Green Street, Urbana, IL 61801, USA    J. Gschwend Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    G. Gutierrez Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    N. W. Halverson Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO, 80309, USA Department of Physics, University of Colorado, Boulder, CO, 80309, USA    I. Harrison Department of Physics, University of Oxford, Denys Wilkinson Building, Keble Road, Oxford OX1 3RH, UK Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK School of Physics and Astronomy, Cardiff University, CF24 3AA, UK    K. Herner Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    S. R. Hinton School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia    G. P. Holder Department of Astronomy, University of Illinois Urbana-Champaign, 1002 West Green Street, Urbana, IL, 61801, USA Department of Physics, University of Illinois Urbana-Champaign, 1110 West Green Street, Urbana, IL, 61801, USA Canadian Institute for Advanced Research, CIFAR Program in Gravity and the Extreme Universe, Toronto, ON, M5G 1Z8, Canada    D. L. Hollowood Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064, USA    W. L. Holzapfel Department of Physics, University of California, Berkeley, CA, 94720, USA    K. Honscheid Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Physics, The Ohio State University, Columbus, OH 43210, USA    J. D. Hrubes University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA    H. Huang Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA Department of Physics, University of Arizona, Tucson, AZ 85721, USA    E. M. Huff Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA    D. Huterer Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA    B. Jain Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    D. J. James Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA    M. Jarvis Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    T. Jeltema Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064, USA    S. Kent Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    L. Knox Department of Physics, University of California, One Shields Avenue, Davis, CA, 95616, USA    A. Kovacs Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain Universidad de La Laguna, Dpto. Astrofísica, E-38206 La Laguna, Tenerife, Spain    E. Krause Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA    K. Kuehn Australian Astronomical Optics, Macquarie University, North Ryde, NSW 2113, Australia Lowell Observatory, 1400 Mars Hill Rd, Flagstaff, AZ 86001, USA    N. Kuropatkin Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    O. Lahav Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK    A. T. Lee Department of Physics, University of California, Berkeley, CA, 94720, USA Physics Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA    P.-F. Leget Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA    P. Lemos Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BN1 9QH, UK    A. R. Liddle Instituto de Astrofísica e Ciências do Espaço, Faculdade de Ciências, Universidade de Lisboa, 1769-016 Lisboa, Portugal    C. Lidman Centre for Gravitational Astrophysics, College of Science, The Australian National University, ACT 2601, Australia The Research School of Astronomy and Astrophysics, Australian National University, ACT 2601, Australia    D. Luong-Van University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA    J. J. McMahon Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA    N. MacCrann Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK    M. March Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    J. L. Marshall George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA    P. Martini Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Astronomy, The Ohio State University, Columbus, OH 43210, USA Radcliffe Institute for Advanced Study, Harvard University, Cambridge, MA 02138    J. McCullough Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA    P. Melchior Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA    F. Menanteau Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA Department of Astronomy, University of Illinois at Urbana-Champaign, 1002 W. Green Street, Urbana, IL 61801, USA    S. S. Meyer Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA    R. Miquel Institució Catalana de Recerca i Estudis Avançats, E-08010 Barcelona, Spain Institut de Física d’Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona) Spain    L. Mocanu Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA    J. J. Mohr Excellence Cluster Universe, Boltzmannstr. 2, 85748 Garching, Germany University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany    R. Morgan Physics Department, 2320 Chamberlin Hall, University of Wisconsin-Madison, 1150 University Avenue Madison, WI 53706-1390    J. Muir Perimeter Institute for Theoretical Physics, 31 Caroline St. North, Waterloo, ON N2L 2Y5, Canada    J. Myles Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    T. Natoli Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada    A. Navarro-Alsina Instituto de Física Gleb Wataghin, Universidade Estadual de Campinas, 13083-859, Campinas, SP, Brazil    R. C. Nichol Department of Physics, University of Surrey, Guildford, UK    Y. Omori Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA    S. Padin California Institute of Technology, 1200 East California Boulevard., Pasadena, CA, 91125, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    S. Pandey Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    Y. Park Kavli Institute for the Physics and Mathematics of the Universe (WPI), UTIAS, The University of Tokyo, Kashiwa, Chiba 277-8583, Japan    F. Paz-Chinchón Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK    M. E. S. Pereira Hamburger Sternwarte, Universität Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany    A. Pieres Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil Observatório Nacional, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil    A. A. Plazas Malagón Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA    A. Porredon Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA Department of Physics, The Ohio State University, Columbus, OH 43210, USA    J. Prat Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    C. Pryke School of Physics and Astronomy, University of Minnesota, 116 Church Street SE Minneapolis, MN, 55455, USA    M. Raveri Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    C. L. Reichardt School of Physics, University of Melbourne, Parkville, VIC 3010, Australia    R. P. Rollins Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK    A. K. Romer Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BN1 9QH, UK    A. Roodman Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    R. Rosenfeld ICTP South American Institute for Fundamental Research
Instituto de Física Teórica, Universidade Estadual Paulista, São Paulo, Brazil
Laboratório Interinstitucional de e-Astronomia - LIneA, Rua Gal. José Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil
   A. J. Ross Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA    J. E. Ruhl Department of Physics, Case Western Reserve University, Cleveland, OH, 44106, USA    E. S. Rykoff Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    C. Sánchez Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA    E. Sanchez Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain    J. Sanchez Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    K. K. Schaffer Liberal Arts Department, School of the Art Institute of Chicago, Chicago, IL, USA 60603 Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA    L. F. Secco Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    I. Sevilla-Noarbe Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas (CIEMAT), Madrid, Spain    E. Sheldon Brookhaven National Laboratory, Bldg 510, Upton, NY 11973, USA    T. Shin Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA    E. Shirokoff Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    M. Smith School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK    Z. Staniszewski Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA Department of Physics, Case Western Reserve University, Cleveland, OH, 44106, USA    A. A. Stark Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA    E. Suchyta Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831    M. E. C. Swanson Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA    G. Tarle Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA    C. To Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA    M. A. Troxel Department of Physics, Duke University Durham, NC 27708, USA    I. Tutusaus Département de Physique Théorique and Center for Astroparticle Physics, Université de Genève, 24 quai Ernest Ansermet, CH-1211 Geneva, Switzerland Institut d’Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain    T. N. Varga Excellence Cluster Origins, Boltzmannstr. 2, 85748 Garching, Germany Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany    J. D. Vieira Department of Astronomy, University of Illinois Urbana-Champaign, 1002 West Green Street, Urbana, IL, 61801, USA Department of Physics, University of Illinois Urbana-Champaign, 1110 West Green Street, Urbana, IL, 61801, USA    N. Weaverdyck Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA    R. H. Wechsler Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    J. Weller Max Planck Institute for Extraterrestrial Physics, Giessenbachstrasse, 85748 Garching, Germany University Observatory, Faculty of Physics, Ludwig-Maximilians-Universität, Scheinerstr. 1, 81679 Munich, Germany    R. Williamson Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91109, USA Department of Astronomy and Astrophysics, University of Chicago, Chicago, IL 60637, USA Kavli Institute for Cosmological Physics, University of Chicago, Chicago, IL 60637, USA    W. L. K. Wu Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA SLAC National Accelerator Laboratory, Menlo Park, CA 94025, USA    B. Yanny Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA    B. Yin Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA    Y. Zhang George P. and Cynthia Woods Mitchell Institute for Fundamental Physics and Astronomy, and Department of Physics and Astronomy, Texas A&M University, College Station, TX 77843, USA    J. Zuntz Institute for Astronomy, University of Edinburgh, Edinburgh EH9 3HJ, UK
(Last updated )
Abstract

We present cosmological constraints from the analysis of two-point correlation functions between galaxy positions and galaxy lensing measured in Dark Energy Survey (DES) Year 3 data and measurements of cosmic microwave background (CMB) lensing from the South Pole Telescope (SPT) and Planck. When jointly analyzing the DES-only two-point functions and the DES cross-correlations with SPT+Planck CMB lensing, we find Ωm=0.344±0.030\Omega_{\rm m}=0.344\pm 0.030 and S8σ8(Ωm/0.3)0.5=0.773±0.016S_{8}\equiv\sigma_{8}(\Omega_{\rm m}/0.3)^{0.5}=0.773\pm 0.016, assuming Λ\LambdaCDM. When additionally combining with measurements of the CMB lensing autospectrum, we find Ωm=0.3060.021+0.018\Omega_{\rm m}=0.306^{+0.018}_{-0.021} and S8=0.792±0.012S_{8}=0.792\pm 0.012. The high signal-to-noise of the CMB lensing cross-correlations enables several powerful consistency tests of these results, including comparisons with constraints derived from cross-correlations only, and comparisons designed to test the robustness of the galaxy lensing and clustering measurements from DES. Applying these tests to our measurements, we find no evidence of significant biases in the baseline cosmological constraints from the DES-only analyses or from the joint analyses with CMB lensing cross-correlations. However, the CMB lensing cross-correlations suggest possible problems with the correlation function measurements using alternative lens galaxy samples, in particular the redMaGiC galaxies and high-redshift MagLim galaxies, consistent with the findings of previous studies. We use the CMB lensing cross-correlations to identify directions for further investigating these problems.

preprint: DES-2021-0649preprint: FERMILAB-PUB-22-475-PPD

I Introduction

The late-time large scale structure (LSS) of the Universe is sensitive to a variety of cosmological signals, ranging from the properties of dark energy and dark matter, to the masses of the neutrinos. Galaxy imaging surveys probe this structure using observations of both the positions of galaxies (which trace the LSS) and the shapes of galaxies (which are distorted by the gravitational lensing effects of the LSS). Several galaxy imaging surveys have used two-point correlations between these measurements to place constraints on cosmological models, including the Kilo Degree Survey (KiDS), the Hyper Suprime Cam Subaru Strategic Program (HSC-SSP), and the Dark Energy Survey (DES) (e.g. Heymans et al., 2021; Hamana et al., 2020; DES Collaboration et al., 2018). DES has recently presented cosmological constraints from the joint analysis of the three two-point correlation functions (33×\times22pt{\rm pt}) between measurements of these probes from the first three years (Y3) of DES data DES Collaboration et al. (2022).

Surveys of the cosmic microwave background (CMB) are also able to probe the late-time LSS through the effects of gravitational lensing. Although CMB photons originate from the last scattering surface at redshift z1100z\sim 1100, their paths are perturbed by structure at late times, including the same LSS measured by galaxy surveys. CMB lensing provides a highly complementary probe of structure to galaxy surveys, and cross-correlations between the two have several appealing features. For one, current galaxy imaging surveys (like DES) identify galaxies out to z1z\sim 1, but the galaxy lensing measurements with these surveys do not have significant sensitivity beyond z0.75z\simeq 0.75. Without the high-redshift lensing information, cosmological constraints from galaxy surveys at z0.75z\gtrsim 0.75 are therefore significantly degraded. CMB lensing, however, reaches peak sensitivity at z2z\sim 2. Therefore, by cross-correlating galaxy surveys with CMB lensing measurements, it is possible to obtain high-precision measurements of the evolution of the matter distribution over a broader range of redshifts than by using galaxy surveys alone. Cross-correlations of galaxy surveys with CMB lensing measurements are also expected to be robust to certain types of systematic biases. Because the galaxy survey measurements are so different from the CMB lensing measurements (e.g., they use data measured by different telescopes at different wavelengths, and use different estimators for the lensing signal), biases in the galaxy surveys are unlikely to correlate with biases in CMB lensing, making two-point functions between the two especially robust. Finally, cross-survey correlations often have different parameter dependencies than correlations within a survey, offering the possibility of improved parameter constraints via degeneracy breaking in joint analyses.

The prospect of obtaining tighter and more robust cosmological constraints from the late-time matter distribution via cross-correlations is particularly timely given recent hints of tensions between some cosmological probes. In particular, recent observations of late-time structure from galaxy surveys tend to prefer lower values of S8σ8Ωm/0.3S_{8}\equiv\sigma_{8}\sqrt{\Omega_{\rm m}/0.3} than CMB surveys Battye et al. (2015); MacCrann et al. (2015); Raveri (2016); Raveri and Hu (2019); Aghanim et al. (2020). This tension could result from physics beyond the standard cosmological constant and cold dark matter model (Λ\LambdaCDM), or it could result from systematic biases in the analyses. By cross-correlating galaxy surveys with CMB lensing, we obtain an independent handle on the late-time large scale structure measurements that can be used to investigate the origins of this possible tension Krolewski et al. (2021); Robertson et al. (2021); Chang et al. (2022); White et al. (2022). Recent analyses have also suggested the possibility of systematic biases in galaxy survey measurements Chang et al. (2019). Because cross-correlations between galaxy surveys and CMB lensing are robust to many important sources of systematic error, they provide a powerful way to ensure that late-time measurements of structure are unbiased.

This work presents the joint cosmological analysis of two-point correlations between galaxy positions and galaxy lensing measured in DES data, and CMB lensing measurements from the South Pole Telescope (SPT, Carlstrom et al., 2011) and the Planck satellite (Planck Collaboration, 2011). As part of its 2008-2011 SPT-SZ survey, SPT obtained high-resolution and high-sensitivity maps of the CMB that partially overlap with the full DES footprint (Omori et al., 2017). At somewhat lower sensitivity and resolution, Planck has obtained full-sky maps of the CMB that overlap completely with the DES footprint. Together, these CMB maps enable high signal-to-noise estimation of the CMB lensing signal across the entire DES footprint Omori et al. (2017, 2022), presenting an opportunity for cross-correlation studies.

From the measurements of galaxy positions (used to compute the galaxy overdensity, δg\delta_{\rm g}), galaxy lensing (γ\gamma, or γt\gamma_{\rm t} for the tangential shear), and CMB lensing (κCMB\kappa_{\rm CMB}), it is possible to form six two-point functions: galaxy clustering (δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle), galaxy-galaxy lensing (δgγ\langle\delta_{\rm g}\gamma\rangle), cosmic shear (γγ\langle\gamma\gamma\rangle), galaxy density-CMB lensing cross-correlation (δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle), galaxy shear-CMB lensing cross-correlation (γκCMB\langle\gamma\kappa_{\rm CMB}\rangle), and the CMB lensing auto-correlation (κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle). All six of the above will be considered here (hereafter, we refer to this combination as 66×\times2pt2{\rm pt}). The five two-point functions excluding the CMB lensing auto-correlation (referred to as 55×\times2pt2{\rm pt}) all probe structure below about at z1.25z\lesssim 1.25, and are highly correlated. This combination, which we measure using DES, SPT and Planck data is the primary focus of this work. The CMB lensing autocorrelation measurements used in this analysis are derived from all-sky Planck data, and are minimally correlated with the 55×\times2pt2{\rm pt} measurements owing to their small (fractional) sky overlap and sensitivity to higher redshifts Omori et al. (2022). We therefore treat the CMB lensing autocorrelation as an external probe, and combine it with 55×\times2pt2{\rm pt} at the likelihood level.

As highlighted above, one of the key reasons to consider cross-correlations of galaxy surveys with CMB lensing is to improve robustness to systematic uncertainties. We will therefore also analyze various subsets of the 66×\times2pt2{\rm pt} probes for the purposes of testing robustness and exploring sensitivity to possible systematic errors. Of particular interest for these tests is the unexpected discovery in DES Y3 data of discrepancies in the galaxy bias values preferred by the clustering and lensing measurements. The DES Y3 analysis considered two galaxy samples for the purposes of measuring δg\delta_{g}: MagLim and redMaGiC. The MagLim galaxies at z0.8z\lesssim 0.8 were used for the baseline cosmological results presented in DES Collaboration et al. (2022). Surprisingly, the galaxy bias values inferred for redMaGiC galaxies from their clustering were found to be roughly 10% lower than the bias values inferred from lensing Pandey et al. (2021), with this discrepancy increasing for the highest-redshift galaxies. There is no known physical explanation for this discrepancy, but tests in Pandey et al. (2021) suggest that it may be connected to observational systematics imparting additional clustering power. MagLim galaxies at high redshift (z0.8z\gtrsim 0.8) also showed a discrepancy between clustering and lensing Porredon et al. (2021). Further investigating these discrepancies is one of the main goals of the present analysis.

The analysis presented here makes several significant improvements relative to previous cross-correlation analyses between DES and SPT/Planck measurements of CMB lensing Baxter et al. (2016); Kirk et al. (2016); Giannantonio et al. (2016); Baxter et al. (2018); Omori et al. (2019a, b); Abbott et al. (2019). First, the DES data have significantly expanded in going from Y1 observations to Y3, covering roughly a factor of three larger area. Second, the CMB lensing maps from SPT/Planck have been remade with several improvements (described in more detail in Omori et al. (2022)). Foremost among these is that we have used the CMB lensing estimator from Madhavacheril and Hill (2018) to reduce contamination in the lensing maps from the thermal Sunyaev-Zel’dovich (tSZ) effect. This contamination was the dominant source of systematic uncertainty for the analysis of (Abbott et al., 2019), and required us to remove a significant fraction of the small-scale measurements from our analysis to ensure that our results were unbiased. As a result, the total signal-to-noise of the CMB lensing cross-correlations was significantly reduced. Using a CMB lensing map that is immune to contamination from the tSZ effect allows us to extract signal from a wider range of angular scales and hence improve our signal-to-noise ratio. Finally, we have also implemented several improvements to the modeling of the correlation functions, which are described in more detail in Omori et al. (2022).

The analysis presented here is the last in a series of three papers: In (Omori et al., 2022, hereafter Paper I) we described the construction of the combined, tSZ-cleaned SPT+Planck CMB lensing map and the methodology of the cosmological analysis. In (Chang et al., 2022, hereafter Paper II), we presented the measurements of the cross-correlation probes δgκCMB+γtκCMB\langle\delta_{g}\kappa_{\rm CMB}\rangle+\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle, a series of diagnostic tests of the measurements, and cosmological constraints from this cross-correlation combination. In this paper (Paper III), we present the joint cosmological constraints from all the 66×\times2pt2{\rm pt} probes, and tests of consistency between various combinations of two-point functions.

The plan of the paper is as follows. In §II we describe the data sets from DES, SPT and Planck that we use in this analysis, and in §III we provide an abridged summary of our model for the correlation function measurements. In §IV, we present cosmological constraints from the joint analysis of cross-correlations between DES and CMB lensing measurements from SPT and Planck, and discuss several tests of the robustness of these constraints enabled by the cross-correlation measurements. We conclude in §V.

II Data from DES, SPT and Planck

DES (Flaugher, 2005) is a photometric survey in five broadband filters (grizYgrizY), with a footprint of nearly 5000deg25000\;{\rm deg}^{2} of the southern sky, imaging hundreds of millions of galaxies. It employs the 570-megapixel Dark Energy Camera (DECam, Flaugher et al., 2015) on the Cerro Tololo Inter-American Observatory (CTIO) 4m Blanco telescope in Chile. We use data from the first three years (Y3) of DES observations. The foundation of the various DES Y3 data products is the Y3 Gold catalog described in Sevilla-Noarbe et al. (2021), which achieves a depth of S/N\sim10 for extended objects up to i\sim23.0 over an unmasked area of 4143deg24143\;{\rm deg}^{2}. In this work, we consider two types of galaxy samples: lens galaxies that are used as biased tracers of the underlying density field, and source galaxies which are used to measure the shape-distorting effects of gravitational lensing. We use the same galaxy samples as in the DES 33×\times22pt{\rm pt} analysis (DES Collaboration et al., 2022). That is, the lens galaxies are taken from the four-redshift bin MagLim sample described in (Porredon et al., 2021), and the source galaxy shapes are taken from the four-redshift bin MetaCalibration sample described in (Gatti et al., 2021). We will additionally consider lens galaxies from the redMaGiC sample described in (Pandey et al., 2021). In particular, we will investigate the potential systematic biases that led to that sample not being used as the baseline cosmology sample in DES Collaboration et al. (2022). The redshift distributions for the MagLim, redMaGiC, Metacalibration samples are shown in Fig. 1.

As mentioned above, we use two CMB lensing maps in this work: one covering the SPT-SZ footprint that uses data from SPT-SZ and Planck (with an overlapping area of 1800\sim 1800 deg2), and a second that covers the northern part of the DES survey that uses only Planck data (with an overlapping area of 2200\sim 2200 deg2). Together, these two CMB lensing maps cover the full DES Y3 survey region. Since the noise levels and beam sizes of SPT-SZ and Planck are different, the resulting CMB lensing maps must be treated separately in our analysis. In Paper II we tested the consistency between the cosmological constraints from these two patches, finding good agreement.

Refer to caption
Figure 1: Redshift distributions of the galaxy samples considered in this work. The MagLim (top panel) and redMaGiC (second from top) lens galaxy samples are used to measure the galaxy overdensity, while the Metacalibration (third from top) source galaxy samples are used to measure weak lensing. Our main cosmological results use only the first four bins of the MagLim sample (solid lines). We perform tests with alternate samples (dashed lines) for exploratory and diagnostic purposes. In the bottom panel we show the lensing kernels (Eq. 3) corresponding to the source galaxies (blue). The grey band in every panel represents the CMB lensing kernel (Eq. 4).

III Modeling and Measurements

The theoretical framework we use in this analysis is laid out in Paper I and Krause et al. (2021). The full 66×\times2pt2{\rm pt} data vector consists of six two-point functions. Since there is little correlation between 55×\times2pt2{\rm pt} and the CMB lensing autocorrelation measurements from Planck, we combine the corresponding constraints at the likelihood level; this approximation is validated in Paper I.

We fit the 66×\times2pt2{\rm pt} data to two different cosmological models: a spatially flat, cosmological constant and cold dark matter model, and a cosmological model where the equation of state parameter of dark energy, ww, is additionally allowed to vary. Following the DES convention, we will refer to these models as Λ\LambdaCDM and wwCDM; note, though, that we allow the sum of the neutrino masses to vary in both of these analyses.

The modeling of the 55×\times2pt2{\rm pt} correlations begins with the auto and cross-power spectrum of the three fields (δg\delta_{\rm g}, γ\gamma, κCMB\kappa_{\rm CMB}). For the correlation functions other than galaxy clustering, we use the Limber approximation Limber (1953):

CXiYj()=𝑑χqXi(χ)qYj(χ)χ2PNL(+1/2χ,z(χ)),C^{X^{i}Y^{j}}(\ell)=\int d\chi\frac{q^{i}_{X}(\chi)q^{j}_{Y}(\chi)}{\chi^{2}}P_{\rm NL}\left(\frac{\ell+1/2}{\chi},z(\chi)\right), (1)

where X,Y{δg,γ,κCMB}X,Y\in\{\delta_{g},\gamma,\kappa_{\rm CMB}\}, i,ji,j labels the redshift bin, PNL(k,z)P_{\rm NL}(k,z) is the non-linear matter power spectrum, which we compute using CAMB and Halofit (Lewis et al., 2000; Takahashi et al., 2012), χ\chi is the comoving distance from the observer, and z(χ)z(\chi) is the redshift corresponding to χ\chi. The weighting functions, q(χ)q(\chi), describe how the different probes respond to large-scale structure at different distances, and are given by

qδgi(χ)=bi(k,z(χ))nδgi(z(χ))dzdχq^{i}_{\delta_{g}}(\chi)=b^{i}(k,z(\chi))n_{\delta_{\rm g}}^{i}(z(\chi))\frac{dz}{d\chi} (2)
qγi(χ)=3H02Ωm2c2χa(χ)χχh𝑑χnγi(z(χ))dzdχχχχ,q^{i}_{\gamma}(\chi)=\frac{3H_{0}^{2}\Omega_{\rm m}}{2c^{2}}\frac{\chi}{a(\chi)}\int_{\chi}^{\chi_{h}}d\chi^{\prime}n_{\gamma}^{i}(z(\chi^{\prime}))\frac{dz}{d\chi^{\prime}}\frac{\chi^{\prime}-\chi}{\chi^{\prime}}, (3)
qκCMB(χ)=3H02Ωm2c2χa(χ)χχχ,q_{\kappa_{\rm CMB}}(\chi)=\frac{3H_{0}^{2}\Omega_{\rm m}}{2c^{2}}\frac{\chi}{a(\chi)}\frac{\chi^{*}-\chi}{\chi^{*}}, (4)

where H0H_{0} and Ωm\Omega_{\rm m} are the Hubble constant and matter density parameters, respectively, a(χ)a(\chi) is the scale factor corresponding to comoving distance χ\chi, b(k,z)b(k,z) is galaxy bias as a function of scale (kk) and redshift, nδg/γi(z)n^{i}_{{\delta_{\rm g}}/\gamma}(z) are the normalized redshift distributions of the lens/source galaxies in bin ii. χ\chi^{*} denotes the comoving distance to the CMB last scattering surface. The model for galaxy clustering is computed without using the Limber approximation, as described in Krause et al. (2021). The angular-space correlation functions are then computed from the auto- and cross-spectra as described in Krause et al. (2021); Omori et al. (2022).

In addition to the basic modeling, described above, we also consider several other physical and observational effects. We list these below but refer the readers to Paper I and Krause et al. (2021) for details.

  • Galaxy bias: Our baseline model assumes linear galaxy bias, but we also explore the potential improvement from using a nonlinear galaxy bias model and including smaller angular scales in our analysis, as described in (Pandey et al., 2021) and Paper I.

  • Intrinsic alignments (IA): We use the Tidal Alignment and Tidal Torquing (TATT, Blazek et al., 2019) model to describe the effect of galaxy intrinsic alignments. We consider an alternate IA model in Appendix C.

  • Lens magnification: Gravitational lensing by foreground mass changes the observed projected number density of lens galaxies as a result of geometric dilution and modulation of galaxy flux and size. We model this effect based on measurements in simulations as described in (Elvin-Poole et al., 2022; Krause et al., 2021).

  • Redshift uncertainties: There are uncertainties associated with the estimation of the redshift distributions of the lens and the source sample, which we model as described in (Myles et al., 2021; Gatti et al., 2022; Cawthon et al., 2020). In (Cordero et al., 2022), an alternate approach to marginalizing over uncertainties in the redshift distributions was also considered, which we explore in Appendix C.

  • Shear calibration uncertainties: We include a prescription for uncertainties in shear calibration as described in DES Collaboration et al. (2022). We estimate uncertainties in the shear measurements using realistic image simulations as described in (MacCrann et al., 2022).

  • CMB map filtering: In order to suppress very small-scale noise in the CMB lensing cross-correlations, we apply filtering to the CMB lensing maps. This filtering is included in the model as described in Paper I.

  • Point mass marginalization: The correlation functions at small scales are impacted by baryonic effects that are challenging to model, such as galaxy formation. This is particularly problematic for δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle: changes in e.g. the masses of the lens galaxies at very small scales can impact the large-scale δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle because tangential shear is a non-local quantity. To reduce sensitivity of our analysis to small-scale effects in δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle, we therefore adopt the point mass marginalization approach of (MacCrann et al., 2020), which involves modifying the covariance matrix of δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle.

We measure the two-point angular correlation functions of the data using the fast tree-based algorithm TreeCorr Jarvis et al. (2004) as described in (Rodríguez-Monroy et al., 2022; Prat et al., 2022; Amon et al., 2022; Secco et al., 2022; Chang et al., 2022). The shear measurements define a spin-2 field on the sky, and there are several ways of decomposing this field for the purposes of measuring two-point functions. For measuring γγ\langle\gamma\gamma\rangle, we use the ξ+\xi_{+} and ξ\xi_{-} decomposition, while for measuring δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle, we consider the correlation only with tangential shear, γt\gamma_{\rm t} (Bartelmann and Schneider, 2001). The covariance matrix associated with these measurements is constructed by combining an analytical halo model covariance, analytical lognormal covariance, and empirical noise estimation from simulations (Friedrich et al., 2021; Omori et al., 2022).

For the final parameter inference, we assume a Gaussian likelihood.111See e.g. (Lin et al., 2020) for tests of the validity of this assumption in the context of cosmic shear, which would also apply here. The priors imposed on the model parameters are shown in Table 2 in Appendix A. The modeling and likelihood framework is built within the CosmoSIS package Zuntz et al. (2015). We generate parameter samples using the nested sampler PolyChord Handley et al. (2015).

Due to uncertainties in the modeling of the correlation functions on small scales (e.g., nonlinear galaxy bias and baryonic effects on the matter power spectrum), in our likelihood analysis we remove the small-scale measurements that could potentially bias our cosmological constraints. The procedure of determining these “scale cuts” is described in (Krause et al., 2021) and Paper I. Note that the choice of angular scales used in the analysis varies somewhat depending on whether we assume a linear or nonlinear galaxy bias model. We focus on the results with linear bias, but consider the results from the nonlinear bias analysis in Appendix B.

In each of the cosmological analyses performed in this work, we include a separate likelihood constructed using a set of ratios of galaxy-galaxy lensing measurements on small scales (Sánchez et al., 2022). These lensing ratios are found to primarily constrain parameters describing the intrinsic alignment model and redshift biases, and are effectively independent of the 55×\times2pt2{\rm pt} data vector.

We utilize two different statistical metrics to assess the consistency of the DES and CMB-lensing cross-correlation measurements, both internally (i.e. between the different two point functions that we measure) and with other cosmological probes. To assess internal consistency, we primarily rely on the posterior predictive distribution (PPD) methods described in Doux et al. (2021). For these assessments, we will quote pp-values, with p<0.01p<0.01 taken as significant evidence of inconsistency. To assess external consistency, we rely on the parameter difference methods developed in Raveri and Doux (2021). For this metric, we will quote differences between parameter constraints in terms of effective σ\sigma values, corresponding to the probability values obtained from the non-Gaussian parameter difference metric. When computing the goodness of fit of our measurements to a particular model, we again rely on the PPD methodology, as discussed in Doux et al. (2021). In this case, the associated pp-values can be thought of as a generalization of the classical pp-value computed from the χ2\chi^{2} statistic that correctly marginalizes over parameter uncertainty.

IV Cosmological constraints

Probe σ8\sigma_{8} Ωm\Omega_{\rm m} S8S_{8} GoF pp-value Comments
33×\times22pt{\rm pt} 0.7330.049+0.0390.733^{+0.039}_{-0.049} 0.3390.031+0.0320.339^{+0.032}_{-0.031} 0.776±0.0170.776\pm 0.017 0.023 DES Collaboration et al. (2022)
δgκCMB+γtκCMB\langle\delta_{g}\kappa_{\rm CMB}\rangle+\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle 0.78±0.070.78\pm 0.07 0.270.05+0.030.27^{+0.03}_{-0.05} 0.74±0.030.74\pm 0.03 0.50 CMB lensing cross-correlations, Paper II
δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle+δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle 0.768±0.0710.768\pm 0.071 0.3030.059+0.0360.303^{+0.036}_{-0.059} 0.765±0.0250.765\pm 0.025 0.063 All cross-correlations, §IV.2.2
55×\times2pt2{\rm pt} 0.7240.043+0.0380.724^{+0.038}_{-0.043} 0.344±0.0300.344\pm 0.030 0.773±0.0160.773\pm 0.016 0.062 §IV.1
66×\times2pt2{\rm pt} 0.785±0.0290.785\pm 0.029 0.306±0.0180.306\pm 0.018 0.792±0.0120.792\pm 0.012 §IV.1
Table 1: Λ\LambdaCDM constraints on Ωm\Omega_{\rm m}, σ8\sigma_{8} and S8σ8(Ωm/0.3)0.5S_{8}\equiv\sigma_{8}(\Omega_{\rm m}/0.3)^{0.5} using different subsets of the 66×\times2pt2{\rm pt} two-point functions. The pp-values correspond to the goodness of fit, as calculated using the PPD methodology. All results here use 4-bin MagLim lens sample and linear galaxy bias. For the 66×\times2pt2{\rm pt} combination, we do not quote a goodness of fit because the CMB lensing autospectrum is treated as an external probe. Rather, we use the parameter difference metric to assess tension between 55×\times2pt2{\rm pt} and κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle (see §III).

IV.1 Baseline cosmological constraints

IV.1.1 Λ\LambdaCDM

We first present constraints on Λ\LambdaCDM from the joint analysis of two-point functions involving DES galaxy position and lensing measurements, and measurements of CMB lensing from SPT and Planck. Following DES Collaboration et al. (2022), all of the results in this subsection use the four redshift bin MagLim lens galaxy sample.

Fig. 2 shows how the constraints from the CMB lensing cross-correlations δgκCMB+γtκCMB\langle\delta_{g}\kappa_{\rm CMB}\rangle+\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle compare to those from 33×\times22pt{\rm pt}. The resulting 68% credible intervals on σ8\sigma_{8}, S8S_{8}, and Ωm\Omega_{\rm m} computed from the marginalized 33×\times22pt{\rm pt} and δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle posteriors are summarized in Table 1. In the same table, we list the goodness of fit pp-values for 33×\times22pt{\rm pt} and δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle + γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle, computed using the PPD formalism. As noted in DES Collaboration et al. (2022), the goodness of fit for 33×\times22pt{\rm pt} alone is not particularly high, but is still above our threshold of p=0.01p=0.01. The goodness of fit for δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle + γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle is acceptable, as described in Chang et al. (2022). While the cross-correlations prefer somewhat lower Ωm\Omega_{\rm m} and higher σ8\sigma_{8}, they are statistically consistent with 33×\times22pt{\rm pt}. Using the PPD formalism, we find p=0.347p=0.347 when comparing the two data subsets, indicating acceptable consistency. We are therefore justified in combining the constraints to form 55×\times2pt2{\rm pt}, shown with the teal contours in the figure.

Given the weaker constraining power of δgκCMB+γtκCMB\langle\delta_{g}\kappa_{\rm CMB}\rangle+\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle relative to 33×\times22pt{\rm pt}, the 55×\times2pt2{\rm pt} constraints are not much tighter than the 33×\times22pt{\rm pt} constraints: we find an improvement of roughly 10% in the precision of the marginalized constraints on Ωm\Omega_{m} and S8S_{8} (see Table 1). The goodness of fit for the full 55×\times2pt2{\rm pt} data vector is p=0.062p=0.062, indicating an acceptable fit.

In Fig. 3 we compare the constraints from 55×\times2pt2{\rm pt} with those from the CMB lensing autospectrum κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle. Owing to the high redshift of the CMB source plane, the CMB lensing-only contour has a different degeneracy direction than 55×\times2pt2{\rm pt}, resulting in a weaker constraint when projecting to the Ωm\Omega_{\rm m} direction, but a comparable constraint in the σ8\sigma_{8} direction. While the CMB lensing autospectrum prefers somewhat higher σ8\sigma_{8} than 55×\times2pt2{\rm pt}, the constraints are generally consistent. Because the CMB lensing autospectrum measurements are treated as an independent probe, we quantify the tension between these measurements and 55×\times2pt2{\rm pt} using the parameter shift metric, finding a difference of 0.8σ0.8\sigma, indicating no evidence of significant tension. We therefore combine the two to generate constraints from all six two-point functions, 66×\times2pt2{\rm pt}, shown with the orange contour in the figure. Due to degeneracy breaking, the joint analysis leads to notably tighter constraints on both Ωm\Omega_{\rm m} and σ8\sigma_{8}. The 1D posterior constraints on these parameters from 66×\times2pt2{\rm pt} are summarized in Table 1. Fig. 3 also shows constraints from Planck measurements of CMB temperature and polarization fluctuations Aghanim et al. (2020). We will assess consistency between our measurements and the Planck measurements in §IV.3.

Refer to caption
Figure 2: Λ\LambdaCDM constraints from the DES Y3 33×\times22pt{\rm pt} measurements (red), cross-correlations between DES Y3 galaxies and shears with SPT+Planck CMB lensing (grey), and from the joint analysis of all five two-point functions (teal). The constraints from 33×\times22pt{\rm pt} are in acceptable agreement with the CMB lensing cross-correlations, justifying the joint analysis of 55×\times2pt2{\rm pt}.
Refer to caption
Figure 3: Λ\LambdaCDM constraints from our 55×\times2pt2{\rm pt} analysis (teal) are compared to those from the Planck CMB lensing autospectrum measurements (grey). The two are in acceptable agreement, justifying the joint analysis of 66×\times2pt2{\rm pt} (orange), which yields significantly tighter constraints due to degeneracy breaking. Also shown are parameter constraints from Planck measurements of primary CMB fluctuations (TT+TE+EE+lowE, dark red).
Refer to caption
Figure 4: Constraints on wwCDM from different combinations of two-point functions. The 55×\times2pt2{\rm pt} constraints (teal) on this model are essentially identical to those of 33×\times22pt{\rm pt} (red). Adding the CMB lensing autospectrum information in the joint 66×\times2pt2{\rm pt} analysis (orange) significantly improves the parameter constraints on wwCDM.

IV.1.2 wwCDM

We now consider constraints on wwCDM, the cosmological model with a constant equation-of-state parameter of dark energy, ww. The constraints from 33×\times22pt{\rm pt}, 55×\times2pt2{\rm pt} and 66×\times2pt2{\rm pt} are shown in Fig. 4. We find that there is little improvement in constraining power on wwCDM when adding the CMB lensing cross-correlations to 33×\times22pt{\rm pt}. Adding the κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle correlation, however, significantly impacts the constraints, presumably because this correlation function adds additional information about structure at z1z\gtrsim 1. The 66×\times2pt2{\rm pt} analysis yields w=0.750.14+0.20w=-0.75^{+0.20}_{-0.14}, S8=0.801±0.013S_{8}=0.801\pm 0.013, and Ωm=0.3540.035+0.041\Omega_{\rm m}=0.354^{+0.041}_{-0.035}. Therefore, the constraints on the dark-energy equation of state parameter are largely consistent with the cosmological-constant scenario of w=1w=-1, and the constraints on Ωm\Omega_{\rm m} and S8S_{8} are consistent with those obtained assuming Λ\LambdaCDM.

IV.2 Robustness tests

In addition to improving cosmological constraints relative to the DES-only 33×\times22pt{\rm pt} analysis, a significant motivation for cross-correlating DES with CMB lensing is to test the robustness of the DES-only constraints. The cross-correlations probe the same large-scale structure as the DES 33×\times22pt{\rm pt} analysis, but with sensitivity to different potential sources of systematic bias, making them powerful cross-checks on the DES results. In this section, we subject the 66×\times2pt2{\rm pt} data vector to several tests of internal consistency.

Refer to caption
Figure 5: Comparison of constraints on Λ\LambdaCDM from 33×\times22pt{\rm pt} (red) with the constraints from the other probes of 66×\times2pt2{\rm pt}, i.e. δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle+κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle(grey). The joint analysis of both (66×\times2pt2{\rm pt}) is shown in orange. The two subsets of the full 66×\times2pt2{\rm pt} analysis are in reasonable agreement. The 66×\times2pt2{\rm pt} analysis prefers higher S8S_{8} than either of the two subsets.

IV.2.1 33×\times22pt{\rm pt} vs. 3×2pt¯\overline{3\!\times\!2{\rm pt}}

We first assess the internal consistency of the 66×\times2pt2{\rm pt} combination of probes by comparing constraints from 33×\times22pt{\rm pt} to the other three two-point functions making up 66×\times2pt2{\rm pt}, which we call 3×2pt¯\overline{3\!\times\!2{\rm pt}} (i.e. δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle+κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle). This comparison is shown in Fig. 5. We find that the constraining power from 3×2pt¯\overline{3\!\times\!2{\rm pt}} is very similar to that of 33×\times22pt{\rm pt}. Because 3×2pt¯\overline{3\!\times\!2{\rm pt}} does not constrain galaxy bias or intrinsic alignment parameters very well, applying the PPD methodology to test consistency between 33×\times22pt{\rm pt} and 3×2pt¯\overline{3\!\times\!2{\rm pt}} is not well motivated. However, we note that we have already tested the consistency of 33×\times22pt{\rm pt} with δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle + γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle (i.e. part of 3×2pt¯\overline{3\!\times\!2{\rm pt}}), finding acceptable agreement (p=0.347p=0.347).

Fig. 5 makes it clear why 66×\times2pt2{\rm pt} prefers a somewhat higher value of S8S_{8} than 33×\times22pt{\rm pt}. It is not the case that 3×2pt¯\overline{3\!\times\!2{\rm pt}} prefers a higher value of S8S_{8} than 33×\times22pt{\rm pt}; indeed, the opposite is true. Rather, the slightly high value of S8S_{8} found for 66×\times2pt2{\rm pt} is caused by the fact that 33×\times22pt{\rm pt} and 3×2pt¯\overline{3\!\times\!2{\rm pt}} have somewhat different degeneracy directions, and intersect at a high value of S8S_{8} for both probes.

IV.2.2 Cross-correlations

Cross-correlations between different observables are generally expected to be more robust to systematic biases than auto-correlations of those observables. Additive systematics that impact a single observable are expected to drop out of cross-correlations with another observable that has uncorrelated systematics. In Fig. 6 we compare the cosmological constraints obtained from only cross-correlations to those from the full 55×\times2pt2{\rm pt}. It is clear that removing the information from the auto-correlations — particularly cosmic shear — degrades the constraints somewhat. However, we find that the value of S8S_{8} inferred only from cross-correlations is consistent with that inferred from the full 55×\times2pt2{\rm pt} analysis. This suggests that additive biases are unlikely to be having a major impact on the DES 33×\times22pt{\rm pt} cosmology results. Using the PPD formalism to evaluate the goodness of fit of the cross-correlations conditioned on the posterior from 55×\times2pt2{\rm pt}, we find p=0.054p=0.054, indicating an acceptable level of consistency between the 55×\times2pt2{\rm pt} constraints and the cross-correlations measurements.

Refer to caption
Figure 6: Comparison of constraints on Λ\LambdaCDM resulting from 55×\times2pt2{\rm pt} (teal) to those that result from only cross-correlations between δg\delta_{g}, γ\gamma and κCMB\kappa_{\rm CMB} (grey). Cross-correlations are expected to be robust to additive systematics that impact only a single field. While some constraining power is lost by removing the auto-correlations, the resulting constraints on S8S_{8} are consistent with those of the baseline analysis, providing a powerful robustness test.

IV.2.3 Lensing only

Refer to caption
Figure 7: Comparison of baseline 55×\times2pt2{\rm pt} constraints on Λ\LambdaCDM (teal) to constraints from various combinations of probes that only involve gravitational lensing. The lensing-only constraints are consistent with our baseline result, suggesting that any systematics which might be impacting the galaxy overdensity measurements are not dramatically biasing our cosmological constraints.

The relationship between galaxy overdensity and the underlying matter field — galaxy bias — presents a significant challenge for analyses of the galaxy distribution. The baseline 33×\times22pt{\rm pt} results presented in DES Collaboration et al. (2022) and the baseline cross-correlation results presented here assume a linear galaxy bias relation when modeling the galaxy field. This model is known to break down at small scales, as investigated for the DES galaxy samples in (Pandey et al., 2021). More complex bias models, such as the perturbation theory-motivated model developed in Pandey et al. (2020), are also expected to have a limited range of validity. There is therefore value in performing analyses that use only lensing information.

Another motivation to consider lensing-only analyses is that the DES galaxy overdensity measurements made with the redMaGiC and high-redshift MagLim galaxies show evidence of systematic biases (i.e. the samples shown with dashed lines in Fig. 1). Measurements of galaxy-galaxy lensing with the redMaGiC galaxies were shown to be inconsistent with clustering measurements using those galaxies Pandey et al. (2021). This inconsistency suggests a potential problem with the redMaGiC overdensity measurements, although it is not clear whether such issues could be impacting the galaxy-galaxy lensing measurements, clustering measurements, or both. Similarly, galaxy-galaxy lensing and clustering measurements with the high-redshift MagLim galaxies were also found to be mutually inconsistent, contributing to a very poor goodness of fit to any of the cosmological models considered. For these reasons, the high-redshift MagLim galaxies were removed from the cosmological analysis in DES Collaboration et al. (2022). These issues, which we investigate further in §IV.2.6, further motivate a cosmological analysis that does not rely on galaxy overdensity measurements.

In Fig. 7, we present cosmological constraints from gravitational lensing only, namely the two-point functions of galaxy lensing and CMB lensing, and their cross-correlation. The lensing-only analysis obtains cosmological constraints that are of comparable precision to those from the full 55×\times2pt2{\rm pt} analysis. We find that the lensing-only analysis yields a constraint on S8S_{8} that is in excellent agreement with the baseline analysis.

IV.2.4 No galaxy lensing

Refer to caption
Figure 8: Comparison of baseline 55×\times2pt2{\rm pt} constraints on Λ\LambdaCDM (teal) to constraints from those combinations of probes that do not rely on galaxy lensing (grey and purple). For reference, we also show the lensing-only constraints — excluding the κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle, which is sensitive to higher redshifts — with the orange curve.

We also consider the constraints that result from those probes that do not involve galaxy lensing. The galaxy lensing measurements could in principle be biased by systematic errors in photometric redshifts of the source galaxies, shear calibration, or an incorrect intrinsic alignment model. Such issues could bias constraints involving galaxy lensing, but would not impact the galaxy overdensity or CMB lensing measurements. Fig. 8 shows the constraints that result only from probes that do not include galaxy lensing (i.e. δg\delta_{g} and κCMB\kappa_{\rm CMB}). Again, we find that the results are consistent with those of 55×\times2pt2{\rm pt}. Fig. 8 also shows the γγ\langle\gamma\gamma\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle constraints for comparison (i.e. lensing only, but excluding κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle, which receives contributions from higher redshifts than the other two-point functions). We find that the constraints involving lens galaxy overdensities are consistent with the lensing-only constraints.

Refer to caption
Refer to caption
Figure 9: Left: Constraints on S8S_{8} and the shear calibration parameters, mim_{i}, from 33×\times22pt{\rm pt} and 55×\times2pt2{\rm pt} using different priors on mim_{i}. With the nominal tight priors on these parameters, 33×\times22pt{\rm pt} (red dashed) and 55×\times2pt2{\rm pt} (teal dashed) yield comparable cosmological constraints. However, when the priors on mim_{i} are substantially weakened, the constraints from 55×\times2pt2{\rm pt} (teal solid) become significantly tighter than those from 33×\times22pt{\rm pt} (red solid). Similarly, the 55×\times2pt2{\rm pt} analysis obtains tighter constraints on the mim_{i} parameters themselves. Right: Same as left panel, but showing constraints on S8S_{8} and Ωm\Omega_{\rm m}. Using the broad mim_{i} priors significantly weakens the cosmological constraints from 33×\times22pt{\rm pt}, but has less of an impact on 55×\times2pt2{\rm pt}.

IV.2.5 Shear calibration

A potentially significant source of systematic uncertainty impacting cosmological constraints from cosmic shear is biases in shear estimation (Hirata and Seljak, 2003). Typically, estimators of lensing shear are calibrated via application to simulated lensed galaxy images. For the DES Year 3 cosmological analysis, calibration of shear biases is described in (MacCrann et al., 2022). While this approach can be used to place tight constraints on shear biases, it has the disadvantage of relying on simulated data. A mismatch between the simulated galaxies used to calibrate the shear estimators and real galaxies could potentially introduce systematic bias.

As pointed out in Vallinotto (2012); Baxter et al. (2016); Schaan et al. (2017), joint analyses of cross-correlations between galaxy surveys and CMB lensing measurements offer the potential of constraining shear calibration biases using only the data. To explore this idea, we repeat our analysis of the 33×\times22pt{\rm pt} and 55×\times2pt2{\rm pt} data vectors using very wide, flat priors on the shear calibration parameters: mi(0.5,0.5)m_{i}\in(-0.5,0.5).

The results of this analysis are shown in Fig. 9. Removing the tight priors on the mim_{i} significantly weakens the cosmological constraints from 33×\times22pt{\rm pt}, especially the constraint on S8S_{8}. This is because both mm and S8S_{8} impact the amplitude of the lensing correlation functions, leading to strong degeneracy between the two. The shear calibration parameters mim_{i} are also very poorly constrained without the tight priors. However, when the CMB lensing cross-correlations are analyzed jointly with 33×\times22pt{\rm pt} (i.e. forming 55×\times2pt2{\rm pt}), the analysis becomes significantly more robust to shear calibration. Removing the priors on mim_{i} weakens the cosmological constraints, but not nearly as much as for 33×\times22pt{\rm pt}: Removing the mm priors degrades the constraints on S8S_{8} by a factor of 4.7 for 33×\times22pt{\rm pt}, but only by a factor of 2.3 for 55×\times2pt2{\rm pt} (see right panel of Fig. 9). The resulting cosmological constraints are consistent with those in the baseline analysis, providing evidence that the DES Y3 33×\times22pt{\rm pt} and 55×\times2pt2{\rm pt} constraints are robust to shear calibration biases. We also find that the 55×\times2pt2{\rm pt} data vector achieves constraints on mm at roughly the 5–10% level depending on the redshift bin, roughly a factor of two improvement over the Y1 analysis presented in (Abbott et al., 2019).

Refer to caption
Figure 10: Posteriors on the linear bias parameters for the MagLim galaxies resulting from different combinations of probes. The parameter bib_{i} represents the linear bias for the iith redshift bin. For the two highest redshift bins (excluded in the baseline cosmology analysis), galaxy clustering (δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle) and galaxy-galaxy lensing (δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle) prefer somewhat different values of the bias, with δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle more in line with the values preferred by clustering.
Refer to caption
Figure 11: Same as Fig. 10, but for redMaGiC galaxies. The bias values preferred by δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle are in good agreement with those preferred by δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle, but show a preference for lower bias values than δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle across all redshifts.

IV.2.6 Investigating the XlensX_{\rm lens} systematic

As noted previously, analyses of δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle and δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle measured with the DES Y3 MagLim Porredon et al. (2021) and redMaGiC Pandey et al. (2021) galaxy samples uncovered discrepancies between the values of galaxy bias preferred by these two correlation functions. The δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle measurements with MagLim galaxies in the two highest redshift bins (i.e. those shown with the dashed lines in the top panel of Fig. 1) prefer higher bias values than δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle by roughly 40 to 60%. Measurements of δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle with the redMaGiC galaxies, on the other hand, show roughly 10% higher bias values than the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle measurements for the first four redshift bins, with this discrepancy increasing to roughly 40% for the highest redshift bin. In principle, some difference between the bias values inferred from δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle and δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle could result from stochastic biasing (e.g. Baumann et al., 2013). However, the amplitude of the difference seen for the redMaGiC galaxies and the high-redshift MagLim galaxies (roughly 10 to 40% percent) is significantly larger than expected from stochasticity (a few percent) Desjacques et al. (2018); Pandey et al. (2020). In Pandey et al. (2021), a new parameter, XlensX_{\rm lens}, was introduced to explore this effect:

Xlensi=bδgγti/bδgδgi,X_{\rm lens}^{i}=b^{i}_{\langle\delta_{g}\gamma_{\rm t}\rangle}/b^{i}_{\langle\delta_{g}\delta_{g}\rangle}, (5)

where bδgγtib^{i}_{\langle\delta_{g}\gamma_{\rm t}\rangle} (bδgδgib^{i}_{\langle\delta_{g}\delta_{g}\rangle}) is the bias parameter for δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle (δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle) in lens galaxy redshift bin ii. The finding that some of the clustering measurements prefer a higher value of galaxy bias than the galaxy-galaxy lensing measurements amounts to a preference for Xlensi<1X_{\rm lens}^{i}<1 when we expect Xlensi=1X_{\rm lens}^{i}=1.

The galaxy-CMB lensing cross-correlations also constrain galaxy bias, providing another handle on the anomalous values of the XlensX_{\rm lens} parameter seen with the redMaGiC and high-redshift MagLim galaxies. We show constraints on the galaxy bias parameters of the MagLim and redMaGiC galaxies from three combinations of probes in Fig. 10 and Fig. 11, respectively. Each of the plotted constraints uses the combination of γγ\langle\gamma\gamma\rangle and γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle — which are effectively independent of the lens galaxies — to constrain the cosmology. The remaining probe is then chosen to be δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle, δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle, or δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle, and this probe is used to constrain the galaxy bias.222This analysis is similar to that presented in Chang et al. (2022), but differs in that we have allowed cosmological parameters to vary, and have included the γγ\langle\gamma\gamma\rangle and γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle measurements (in effect letting the data constrain the cosmological model).

For the two highest redshift bins of MagLim galaxies, we see from Fig. 10 that the δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle measurements prefer higher values of galaxy bias than the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle measurements, consistent with the preference for Xlensi<1X_{\rm lens}^{i}<1 described above. Interestingly, it appears that the δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle measurements prefer galaxy bias values more in line with the δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle measurements. This suggest that the preference for Xlensi<1X_{\rm lens}^{i}<1 is likely driven by δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle. This is perhaps not surprising, given the large residuals of the model fits to δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle seen in Porredon et al. (2021). However, note that there is no obvious reason for a possible failure of the baseline model to fit δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle. As a cross-correlation, the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle measurements are expected to be quite robust to many observational systematics. Moreover, any systematic impacting δg\delta_{\rm g} would likely show up even more strongly in δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle, and any systematic impacting γ\gamma would likely show up more strongly in γγ\langle\gamma\gamma\rangle. Another possibility is a failure in modeling some physical effect. One such effect is lens magnification, which is known to have a significant impact on the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle correlations at high redshifts Porredon et al. (2021).

Fig. 11 shows the analogous bias constraints for redMaGiC galaxies. In this case, we see that δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle and δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle measurements both prefer consistently lower values of galaxy bias than δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle, with this difference particularly pronounced in the last redshift bin. This suggests that a possible cause of the redMaGiC preference for Xilens<1X^{\rm lens}_{i}<1 is in the δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle measurements. In the case of δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle, it is possible that some observational systematic is modulating the redMaGiC galaxy overdensity field, resulting in a higher than expected clustering amplitude and thus a preference for higher galaxy bias. Such a systematic in the δg\delta_{\rm g} measurements would be expected to have a less noticeable impact on δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle. At the same time, it should be emphasized that the analysis of Rodríguez-Monroy et al. (2022) extensively tested the redMaGiC sample for possible contamination by various observational systematics. While some correlation of known systematics with galaxy density is detected, this correlation is corrected using galaxy re-weighting. It therefore appears to be difficult to explain the anomalous XlensX^{\rm lens} values with any known observational systematic.

Refer to caption
Figure 12: Cosmological constraints on wwCDM from the 33×\times22pt{\rm pt} data vector measured with the MagLim (red) and redMaGiC (grey) lens galaxy samples. The constraints from redMaGiC prefer surprisingly less negative ww, as discussed in DES Collaboration et al. (2022). However, when the redMaGiC clustering measurements (δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle) are replaced by δgκCMB+γtκCMB\langle\delta_{g}\kappa_{\rm CMB}\rangle+\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle to form a combination of four two-point functions (orange), the constraints agree better with those of MagLim.

The interpretation of the redMaGiC preference for Xlens<1X^{\rm lens}<1 in terms of a systematic impacting redMaGiC δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle measurements is supported by tests with a modified redMaGiC galaxy sample presented in Pandey et al. (2021). The nominal redMaGiC galaxy sample is selected by requiring that galaxies match a red sequence template, as measured by χ2\chi^{2}. In Pandey et al. (2021), an alternative, “broad χ2\chi^{2}” sample of galaxies was selected by relaxing the χ2\chi^{2} threshold for selection. One would expect that if an observational systematic is modulating the photometry of galaxies, it should have a smaller impact on the “broad χ2\chi^{2}” sample than on the nominal sample. Indeed, it was found that for this alternate sample, the preference for Xlensi<1X_{\rm lens}^{i}<1 seen for the first four redshift bins disappears. While it might seem surprising that the preference for Xlensi<1X_{\rm lens}^{i}<1 is possibly driven by two different factors for MagLim and redMaGiC galaxies, this interpretation seems consistent with the observed redshift trends. It may be that observational systematics in redMaGiC galaxy selection are impacting the bias values inferred from δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle at low redshift, while problems in modeling δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle are impacting the bias values inferred for MagLim from δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle at high redshift. The redMaGiC galaxies may be less affected by this latter systematic, as they do not extend to the high redshifts probed by the last two redshift bins of the MagLim sample. We note, though, that even for redMaGiC, the CMB lensing cross-correlations prefer higher galaxy bias than δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle in the highest redshift bin; this could be suggesting that the same problem impacting the high-redshift MagLim galaxies is impacting the high-redshift redMaGiC galaxies. This interpretation would be consistent with mismodeling of δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle at high redshift.

The impact of the apparent systematic in the redMaGiC sample is also noticeable when the cosmological model is changed from Λ\LambdaCDM to wwCDM. While the redMaGiC 33×\times22pt{\rm pt} constraints on Λ\LambdaCDM are quite robust to allowing the XlensX_{\rm lens} parameter to vary, the constraints on wwCDM shift significantly when this additional freedom is introduced. This is perhaps not surprising given that the systematic biases with redMaGiC appear to be redshift-dependent, and might therefore be somewhat degenerate with the effects of ww.

Refer to caption
Figure 13: Comparison of the cosmological constraints resulting from different combinations of two-point functions involving DES measurements of galaxy positions and lensing, and SPT+Planck measurements of CMB lensing. We also show (bottom row) constraints from Planck-only measurements of the primary CMB fluctuations.

Since our analysis above suggests that the problems with redMaGiC may be isolated to the clustering measurements, in Fig. 12 we present constraints on wwCDM from the 55×\times2pt2{\rm pt} combination of probes without the clustering measurements. Interestingly, we see that there is a significant shift in the constraints on ww relative to the 33×\times22pt{\rm pt} analysis. The constraints without the clustering measurements are in good agreement with the MagLim constraints. This lends additional support to the idea that the redMaGiC clustering measurements may be systematically biased.

To summarize the above discussion, our analysis with CMB lensing cross-correlations suggests that there may be two different sources for the XlensX_{\rm lens} systematic seen with redMaGiC and MagLim galaxies. For redMaGiC galaxies, our analysis suggests a possible bias in the clustering measurements across all redshift bins. Such a bias could conceivably be caused by some observational systematic impacting the redMaGiC selection, which would be consistent with tests performed in Pandey et al. (2021). At the same time, high-redshift MagLim galaxies (and possibly high-redshift redMaGiC galaxies as well) show evidence of a potentially different systematic error that favors a problem with the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle fits. Such an issue could conceivably be caused by a problem with the δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle modeling, such as an incorrect prescription for magnification effects, which become more pronounced at high redshifts.

IV.3 Consistency with Planck primary CMB measurements

As seen in Fig. 3, we find that the cosmological constraints on Λ\LambdaCDM from 55×\times2pt2{\rm pt} and 66×\times2pt2{\rm pt} are not in significant tension with the constraints from the primary CMB measurements of Planck. In particular, we compare our constraints to those from the combination of Planck TTTT, TETE, EEEE and low-\ell EE-mode polarization measurements (Planck TT+TE+EE+lowE) Aghanim et al. (2020). Note that we do not include Planck measurements of the CMB lensing power spectrum in this combination. Using the tension metric of Raveri and Doux (2021), we find that the 33×\times22pt{\rm pt}, 55×\times2pt2{\rm pt}, and 66×\times2pt2{\rm pt} constraints are in agreement with Planck at the level of 1.5σ1.5\sigma, 1.4σ1.4\sigma, and 1.4σ1.4\sigma, respectively. The fact that 33×\times22pt{\rm pt} and 55×\times2pt2{\rm pt} are roughly equally consistent with Planck is not surprising, given that the 55×\times2pt2{\rm pt} constraints are quite close to those of 33×\times22pt{\rm pt}. Interestingly, while the 66×\times2pt2{\rm pt} constraints are significantly tighter than 55×\times2pt2{\rm pt}, the level of consistency with Planck remains roughly the same. This results from the preference by κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle for somewhat higher values of σ8\sigma_{8}, as seen in Fig. 3. Fig. 13 directly compares the S8S_{8} and Ωm\Omega_{\rm m} constraints from these and other two-point function combinations, assuming Λ\LambdaCDM. We note that for consistency with our analysis, we vary the sum of the neutrino masses and impose the priors shown in Table 2 when generating the Planck primary CMB constraints shown in this figure.

V Summary

We have presented cosmological constraints from an analysis of two-point correlation functions between measurements of galaxy positions and galaxy lensing from DES Y3 data, and CMB lensing measurements from SPT and Planck. Our main cosmological constraints are summarized in Table 1.

The high signal-to-noise of the CMB lensing cross-correlation measurements using DES Y3, SPT-SZ and Planck data enables powerful robustness tests of our cosmological constraints. The results of several of these tests are shown in Fig. 13. We summarize the main findings of these tests below:

  • The goodness of fit of Λ\LambdaCDM to the 55×\times2pt2{\rm pt} data vector is acceptable (p=0.062p=0.062), and the corresponding parameter constraints are consistent with those from κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle measurements by Planck.

  • Using only cross-correlations between DES and CMB lensing, we obtain constraints on S8S_{8} that are comparable in precision and consistent with the baseline 55×\times2pt2{\rm pt} results. This result suggests that additive systematics are not significantly impacting the 55×\times2pt2{\rm pt} cosmological constraints.

  • Using only gravitational lensing (i.e. no information from galaxy overdensities) yields constraints in agreement with the baseline results. This result suggests that potential systematics impacting the DES galaxy samples, as well as modeling of galaxy bias, are not significantly biasing the 55×\times2pt2{\rm pt} cosmological constraints.

  • The cosmological constraints from two-point functions of MagLim galaxy overdensity measurements and CMB lensing are generally consistent with the baseline 55×\times2pt2{\rm pt} analysis. This result suggests that shear systematics and modeling of galaxy lensing are not significantly biasing the 55×\times2pt2{\rm pt} cosmological constraints. We do, however, observe a low-significance increase in S8S_{8} when considering only those two-point functions that do not involve galaxy lensing. This shift is driven by the intersection of the δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle + δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle and κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle constraints, and is not present when considering δgδg\langle\delta_{\rm g}\delta_{\rm g}\rangle + δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle alone.

  • Without priors on shear calibration, the cosmological constraints on S8S_{8} from 55×\times2pt2{\rm pt} are in good agreement with the baseline 55×\times2pt2{\rm pt} results. The data calibrate the shear bias parameters at the 5–10% level, and yield constraints consistent with our nominal priors. These results suggest that shear calibration biases are not significantly impacting the 55×\times2pt2{\rm pt} cosmological constraints.

  • The constraints on Ωm\Omega_{\rm m} from the different analysis variations are generally consistent. Although the analysis of δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle+κCMBκCMB\langle\kappa_{\rm CMB}\kappa_{\rm CMB}\rangle prefers a somewhat lower value of Ωm\Omega_{\rm m}, this combination of probes is statistically consistent with 33×\times22pt{\rm pt}.

The cosmological constraints from the 33×\times22pt{\rm pt}, 55×\times2pt2{\rm pt}, and 66×\times2pt2{\rm pt} analyses therefore appear remarkably robust to possible systematic biases.

Assessing the consistency between our constraint on Λ\LambdaCDM and those of Planck, we find that the 55×\times2pt2{\rm pt} and 66×\times2pt2{\rm pt} constraints are statistically consistent with Planck at the 1.4σ1.4\sigma level, as assessed using the full, multi-dimensional posteriors from these measurements. As seen in Fig. 13, however, essentially all combinations of two point functions that we consider prefer lower S8S_{8} values than Planck. Note, though, that there is significant covariance between some of these measurements.

We have also investigated possible issues with the analysis of alternate lens galaxy samples, namely the high-redshift MagLim galaxies and the redMaGiC galaxies. Evidence for biases when analyzing correlation functions measured with these samples was found previously in Porredon et al. (2021), Pandey et al. (2021), Elvin-Poole et al. (2022), and DES Collaboration et al. (2022). The CMB lensing cross-correlations considered here provide a powerful way to probe the sources of these biases. In the context of Λ\LambdaCDM, our analysis of CMB lensing cross-correlations suggests a possible problem in the modeling of δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle at high redshift for the MagLim galaxies, and possibly the redMaGiC galaxies as well. At the same time, the δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle measurements with redMaGiC suggest a possible observational systematic that impacts redMaGiC galaxy clustering across all redshifts. This interpretation is supported by tests with an alternate redMaGiC galaxy sample in Pandey et al. (2021). In the context of wwCDM, the 33×\times22pt{\rm pt} measurements with redMaGiC have previously shown to yield constraints inconsistent with the MagLim analysis, and a preference for surprisingly less negative ww. We show that analysis of γγ\langle\gamma\gamma\rangle+δgγt\langle\delta_{\rm g}\gamma_{\rm t}\rangle+δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle+γtκCMB\langle\gamma_{\rm t}\kappa_{\rm CMB}\rangle (i.e. two-point functions between DES and CMB lensing, excluding galaxy clustering) measured with redMaGiC yields cosmological constraints that are in better agreement with MagLim, and do not show a strong preference for w>1w>-1. Finally, we note that while the analyses presented here suggest possible interpretations of the XlensX_{\rm lens} bias, more work with current and future DES data is needed to clarify the true source of this systematic uncertainty.

As the data volume and quality from cosmological surveys continue to improve, we expect similar cross-correlation analyses between galaxy surveys and CMB lensing measurements to play an important role in constraining late-time large scale structure. Excitingly, we expect constraints from such measurements to improve dramatically in the very near future with Year 6 data from DES and new CMB lensing maps from SPT-3G Sobrin et al. (2022) and AdvACT Henderson et al. (2016). These measurements should help to provide a clearer picture of any possible S8S_{8} tension. Looking farther forward, cross-correlations between surveys such as the Vera Rubin Observatory Legacy Survey of Space and Time Ivezić et al. (2019); LSST Science Collaboration et al. (2009), the Nancy Grace Roman Space Telescope Doré et al. (2019) , the ESA Euclid mission Laureijs et al. (2011), Simons Observatory Ade et al. (2019), and CMB-S4 Abazajian et al. (2016) will enable significantly more powerful cross-correlation studies that will deliver some of the most precise and accurate cosmological constraints, and that will allow us to continue stress-testing the concordance Λ\LambdaCDM model.

Acknowledgements.
The South Pole Telescope program is supported by the National Science Foundation (NSF) through the grant OPP-1852617. Partial support is also provided by the Kavli Institute of Cosmological Physics at the University of Chicago. Argonne National Laboratory’s work was supported by the U.S. Department of Energy, Office of Science, Office of High Energy Physics, under contract DE-AC02- 06CH11357. Work at Fermi National Accelerator Laboratory, a DOE-OS, HEP User Facility managed by the Fermi Research Alliance, LLC, was supported under Contract No. DE-AC02- 07CH11359. The Melbourne authors acknowledge support from the Australian Research Council’s Discovery Projects scheme (DP210102386). The McGill authors acknowledge funding from the Natural Sciences and Engineering Research Council of Canada, Canadian Institute for Advanced research, and the Fonds de recherche du Quúbec Nature et technologies. The CU Boulder group acknowledges support from NSF AST-0956135. The Munich group acknowledges the support by the ORIGINS Cluster (funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany’s Excellence Strategy – EXC-2094 – 390783311), the MaxPlanck-Gesellschaft Faculty Fellowship Program, and the Ludwig-Maximilians-Universität München. JV acknowledges support from the Sloan Foundation. Funding for the DES Projects has been provided by the U.S. Department of Energy, the U.S. National Science Foundation, the Ministry of Science and Education of Spain, the Science and Technology Facilities Council of the United Kingdom, the Higher Education Funding Council for England, the National Center for Supercomputing Applications at the University of Illinois at Urbana-Champaign, the Kavli Institute of Cosmological Physics at the University of Chicago, the Center for Cosmology and Astro-Particle Physics at the Ohio State University, the Mitchell Institute for Fundamental Physics and Astronomy at Texas A&M University, Financiadora de Estudos e Projetos, Fundação Carlos Chagas Filho de Amparo à Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Científico e Tecnológico and the Ministério da Ciência, Tecnologia e Inovação, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey. The Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energéticas, Medioambientales y Tecnológicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenössische Technische Hochschule (ETH) Zürich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ciències de l’Espai (IEEC/CSIC), the Institut de Física d’Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians Universität München and the associated Excellence Cluster Universe, the University of Michigan, NFS’s NOIRLab, the University of Nottingham, The Ohio State University, the University of Pennsylvania, the University of Portsmouth, SLAC National Accelerator Laboratory, Stanford University, the University of Sussex, Texas A&M University, and the OzDES Membership Consortium. Based in part on observations at Cerro Tololo Inter-American Observatory at NSF’s NOIRLab (NOIRLab Prop. ID 2012B-0001; PI: J. Frieman), which is managed by the Association of Universities for Research in Astronomy (AURA) under a cooperative agreement with the National Science Foundation. The DES data management system is supported by the National Science Foundation under Grant Numbers AST-1138766 and AST-1536171. The DES participants from Spanish institutions are partially supported by MICINN under grants ESP2017-89838, PGC2018-094773, PGC2018-102021, SEV-2016-0588, SEV-2016-0597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded by the CERCA program of the Generalitat de Catalunya. Research leading to these results has received funding from the European Research Council under the European Union’s Seventh Framework Program (FP7/2007-2013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciência e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2). This manuscript has been authored by Fermi Research Alliance, LLC under Contract No. DE-AC02-07CH11359 with the U.S. Department of Energy, Office of Science, Office of High Energy Physics. We gratefully acknowledge the computing resources provided on Crossover ( and/or Bebop and/or Swing and/or Blues), a high-performance computing cluster operated by the Laboratory Computing Resource Center at Argonne National Laboratory.

Appendix A Parameter priors

In Table 2 we list the priors used in our analysis.

Parameter Prior
Ωm\Omega_{\rm m} 𝒰[0.1,0.9]\mathcal{U}[0.1,0.9]
As×109A_{\rm{s}}\times 10^{9} 𝒰[0.5,5.0]\mathcal{U}[0.5,5.0]
nsn_{\rm{s}} 𝒰[0.87,1.07]\mathcal{U}[0.87,1.07]
Ωb\Omega_{\rm{b}} 𝒰[0.03,0.07]\mathcal{U}[0.03,0.07]
hh 𝒰[0.55,0.91]\mathcal{U}[0.55,0.91]
Ωνh2×104\Omega_{\nu}h^{2}\times 10^{4} 𝒰[6.0,64.4]\mathcal{U}[6.0,64.4]
a1a_{1} 𝒰[5.0,5.0]\mathcal{U}[-5.0,5.0]
a2a_{2} 𝒰[5.0,5.0]\mathcal{U}[-5.0,5.0]
η1\eta_{1} 𝒰[5.0,5.0]\mathcal{U}[-5.0,5.0]
η2\eta_{2} 𝒰[5.0,5.0]\mathcal{U}[-5.0,5.0]
btab_{\rm{ta}} 𝒰[0.0,2.0]\mathcal{U}[0.0,2.0]
MagLim
b16b^{1\cdots 6} 𝒰[0.8,3.0]\mathcal{U}[0.8,3.0]
b116b_{1}^{1\cdots 6} 𝒰[0.67,3.0]\mathcal{U}[0.67,3.0]
b216b_{2}^{1\cdots 6} 𝒰[4.2,4.2]\mathcal{U}[-4.2,4.2]
Cl16C_{\rm l}^{1\cdots 6} δ(0.42)\delta(0.42), δ(0.3)\delta(0.3), δ(1.76)\delta(1.76), δ(1.94)\delta(1.94), δ(1.56)\delta(1.56), δ(2.96)\delta(2.96)
Δz16×102\Delta_{z}^{1...6}\times 10^{2} 𝒩[0.9,0.7]\mathcal{N}[-0.9,0.7], 𝒩[3.5,1.1]\mathcal{N}[-3.5,1.1], 𝒩[0.5,0.6]\mathcal{N}[-0.5,0.6], 𝒩[0.7,0.6]\mathcal{N}[-0.7,0.6], 𝒩[0.2,0.7]\mathcal{N}[0.2,0.7] , 𝒩[0.2,0.8]\mathcal{N}[0.2,0.8]
σz16\sigma_{z}^{1...6} 𝒩[0.98,0.062]\mathcal{N}[0.98,0.062], 𝒩[1.31,0.093]\mathcal{N}[1.31,0.093], 𝒩[0.87,0.054]\mathcal{N}[0.87,0.054], 𝒩[0.92,0.05]\mathcal{N}[0.92,0.05], 𝒩[1.08,0.067]\mathcal{N}[1.08,0.067], 𝒩[0.845,0.073]\mathcal{N}[0.845,0.073]
redMaGiC
b15b^{1\cdots 5} 𝒰[0.8,3.0]\mathcal{U}[0.8,3.0]
b115b_{1}^{1\cdots 5} 𝒰[0.67,2.52]\mathcal{U}[0.67,2.52]
b215b_{2}^{1\cdots 5} 𝒰[3.5,3.5]\mathcal{U}[-3.5,3.5]
Cl15C_{\rm l}^{1\cdots 5} δ(0.62)\delta(0.62), δ(3.04)\delta(-3.04), δ(1.32)\delta(-1.32), δ(2.5)\delta(2.5), δ(1.94)\delta(1.94)
Δz15×102\Delta_{z}^{1...5}\times 10^{2} 𝒩[0.6,0.4]\mathcal{N}[0.6,0.4], 𝒩[0.1,0.3]\mathcal{N}[0.1,0.3], 𝒩[0.4,0.3]\mathcal{N}[0.4,0.3], 𝒩[0.2,0.5]\mathcal{N}[-0.2,0.5], 𝒩[0.7,1.0]\mathcal{N}[-0.7,1.0]
σz14\sigma_{z}^{1...4} δ(1)\delta(1), δ(1)\delta(1), δ(1)\delta(1), δ(1)\delta(1), 𝒩[1.23,0.054]\mathcal{N}[1.23,0.054]
MetaCalibration
m14×103m^{1...4}\times 10^{3} 𝒩[6.0,9.1]\mathcal{N}[-6.0,9.1], 𝒩[20.0,7.8]\mathcal{N}[-20.0,7.8], 𝒩[24.0,7.6]\mathcal{N}[-24.0,7.6], 𝒩[37.0,7.6]\mathcal{N}[-37.0,7.6]
Δz14×102\Delta_{z}^{1...4}\times 10^{-2} 𝒩[0.0,1.8]\mathcal{N}[0.0,1.8], 𝒩[0.0,1.5]\mathcal{N}[0.0,1.5], 𝒩[0.0,1.1]\mathcal{N}[0.0,1.1], 𝒩[0.0,1.7]\mathcal{N}[0.0,1.7]
Table 2: Prior values for cosmological and nuisance parameters included in our model. For the priors, 𝒰[a,b]\mathcal{U}[a,b] indicates a uniform prior between aa and bb, while 𝒩[a,b]\mathcal{N}[a,b] indicates a Gaussian prior with mean aa and standard deviation bb. δ(a)\delta(a) is a Dirac Delta function at value aa, which effectively means that the parameter is fixed at aa. Note that the fiducial lens sample is the first 4 bins of the MagLim sample. The two high-redshift MagLim bins and the redMaGiC sample are shown in grey to indicate they are not part of the fiducial analysis.

Appendix B Adding small-scale information with nonlinear galaxy bias

Our baseline analysis adopts a linear galaxy bias model to describe the relationship between the galaxy overdensity and the underlying matter field. At small scales, this description of galaxy biasing is known to break down. The breakdown in linear galaxy bias drives our choice of angular scales used to analyzing the δgκCMB\langle\delta_{\rm g}\kappa_{\rm CMB}\rangle correlation, as described in Paper I. By adopting a higher-order bias model, it is possible to include smaller angular scales in the cosmological analysis and potentially improve parameter constraints. At the same time, a more complex bias model necessitates more free parameters, which degrades the parameter constraints to some extent. We now consider the parameter constraints from 55×\times2pt2{\rm pt} using the nonlinear galaxy bias model described in Pandey et al. (2021).

The constraints from this analysis are presented in Fig. 14. We find that adopting a nonlinear description of galaxy bias improves the precision of the constraint on S8S_{8} by roughly 10%, and the precision of the constraints on both Ωm\Omega_{m} and S8S_{8} by roughly 10%.

Refer to caption
Figure 14: Parameter constraints obtained when using a nonlinear galaxy bias model to analyze the 55×\times2pt2{\rm pt} data vector (grey) compared to our baseline 55×\times2pt2{\rm pt} analysis (teal), which adopts a linear bias model. The nonlinear bias analysis can be used to fit smaller scales of measured correlation functions resulting in improved constraints.
Refer to caption
Figure 15: Parameter constraints obtained when using alternative prescriptions for modeling photometric redshift biases and intrinsic alignments. The teal curves show our baseline results, while the grey dashed curves show results assuming the hyperrank method for calibrating the source galaxy redshift distributions, and the grey solid curves show results assuming the NLA intrinsic alignment model (rather than the baseline TATT model). In both cases, there are minimal shifts relative to our baseline results.

Appendix C Alternative redshift calibration and IA model

Our baseline analysis assumes that uncertainties in the source galaxy redshift distributions are characterized by shift and stretch parameters, as described in (Myles et al., 2021; Gatti et al., 2022). An alternative approach to characterizing the uncertainties in the redshift distributions is hyperrank, described in Cordero et al. (2022). Rather than attempt to parameterize biases in the redshift distributions, hyperrank provides a way to sample over realizations of the full posteriors on these distributions. Repeating our analysis of the 55×\times2pt2{\rm pt} data using this alternative redshift uncertainty prescription yields the constraints shown in Fig. 15. Although there is a small shift in S8S_{8}, it is well within our uncertainties.

The intrinsic alignment (IA) model that we adopt in our baseline analysis is TATT (TATT, Blazek et al., 2019). In Fig. 15, we show the results of instead adopting the nonlinear alignment model (NLA, Bridle and King, 2007). The NLA model is more restrictive than TATT in the sense that the latter becomes equivalent to the former in the limit that a2=η2=bta=0a_{2}=\eta_{2}=b_{\rm ta}=0. We find that switching to NLA results in minimal changes to the parameter constraints from 55×\times2pt2{\rm pt}.

References