We present cosmological constraints from the abundance of galaxy clusters selected via the thermal SunyaevZeldovich SZ effect in South Pole Telescope SPT data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey DES and the Hubble Space Telescope HST. The cluster sample is constructed from the combined SPTSZ SPTpol ECS and SPTpol 500d surveys and comprises 1005 confirmed clusters in the redshift range 0.251.78 over a total sky area of inlineformulammlmath displayinlinemmlmrowmmlmn5200mmlmnmmlmtext mmlmtextmmlmtext mmlmtextmmlmrowmmlmsupmmlmrowmmlmidegmmlmimmlmrowmmlmrowmmlmn2mmlmnmmlmrowmmlmsupmmlmrowmmlmrowmmlmathinlineformula. We use DES Year 3 weaklensing data for 688 clusters with redshifts inlineformulammlmath displayinlinemmlmizmmlmimmlmoltmmlmommlmn0.95mmlmnmmlmathinlineformula and HST weaklensing data for 39 clusters with inlineformulammlmath displayinlinemmlmn0.6mmlmnmmlmoltmmlmommlmizmmlmimmlmoltmmlmommlmn1.7mmlmnmmlmathinlineformula. The weaklensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observablemass relation without having to make strong assumptions about e.g. the hydrodynamical state of the clusters. For a flat inlineformulammlmath displayinlinemmlmi mathvariantnormalmmlmimmlmiCDMmmlmimmlmathinlineformula cosmology and marginalizing over the sum of massive neutrinos we measure inlineformulammlmath displayinlinemmlmsubmmlmi mathvariantnormalmmlmimmlmi mathvariantnormalmmmlmimmlmsubmmlmommlmommlmn0.286mmlmnmmlmommlmommlmn0.032mmlmnmmlmathinlineformula inlineformulammlmath displayinlinemmlmsubmmlmimmlmimmlmn8mmlmnmmlmsubmmlmommlmommlmn0.817mmlmnmmlmommlmommlmn0.026mmlmnmmlmathinlineformula and the parameter combination inlineformulammlmath displayinlinemmlmrowmmlmsubmmlmrowmmlmimmlmimmlmrowmmlmrowmmlmn8mmlmnmmlmrowmmlmsubmmlmo stretchyfalsemmlmommlmsubmmlmrowmmlmi mathvariantnormalmmlmimmlmrowmmlmrowmmlmi mathvariantnormalmmmlmimmlmrowmmlmsubmmlmommlmommlmn0.3mmlmnmmlmsupmmlmrowmmlmo stretchyfalsemmlmommlmrowmmlmrowmmlmn0.25mmlmnmmlmrowmmlmsupmmlmommlmommlmn0.805mmlmnmmlmommlmommlmn0.016mmlmnmmlmrowmmlmathinlineformula. Our measurement of inlineformulammlmath displayinlinemmlmsubmmlmiSmmlmimmlmn8mmlmnmmlmsubmmlmommlmommlmsubmmlmimmlmimmlmn8mmlmnmmlmsubmmlmsqrtmmlmrowmmlmsubmmlmi mathvariantnormalmmlmimmlmi mathvariantnormalmmmlmimmlmsubmmlmommlmommlmn0.3mmlmnmmlmrowmmlmsqrtmmlmommlmommlmn0.795mmlmnmmlmommlmommlmn0.029mmlmnmmlmathinlineformula and the constraint from Planck CMB anisotropies 2018 TT TE inlineformulammlmath displayinlinemmlmrowmmlmiEEmmlmimmlmommlmommlmilowEmmlmimmlmrowmmlmathinlineformula differ by inlineformulammlmath displayinlinemmlmn1.1mmlmnmmlmimmlmimmlmathinlineformula. In combination with that Planck dataset we place a 95 upper limit on the sum of neutrino masses inlineformulammlmath displayinlinemmlmommlmommlmsubmmlmimmmlmimmlmimmlmimmlmsubmmlmoltmmlmommlmn0.18mmlmnmmlmtext mmlmtextmmlmtext mmlmtextmmlmieVmmlmimmlmathinlineformula. When additionally allowing the dark energy equation of state parameter inlineformulammlmath displayinlinemmlmiwmmlmimmlmathinlineformula to vary we obtain inlineformulammlmath displayinlinemmlmiwmmlmimmlmommlmommlmommlmommlmn1.45mmlmnmmlmommlmommlmn0.31mmlmnmmlmathinlineformula from our clusterbased analysis. In combination with Planck data we measure inlineformulammlmath displayinlinemmlmiwmmlmimmlmommlmommlmommlmommlmn1.3mmlmnmmlmsubsupmmlmn4mmlmnmmlmrowmmlmommlmommlmn0.15mmlmnmmlmrowmmlmrowmmlmommlmommlmn0.22mmlmnmmlmrowmmlmsubsupmmlmathinlineformula or a inlineformulammlmath displayinlinemmlmn2.2mmlmnmmlmimmlmimmlmathinlineformula difference with a cosmological constant. We use the cluster abundance to measure inlineformulammlmath displayinlinemmlmsubmmlmimmlmimmlmn8mmlmnmmlmsubmmlmathinlineformula in five redshift bins between 0.25 and 1.8 and we find the results to be consistent with structure growth as predicted by the inlineformulammlmath displayinlinemmlmi mathvariantnormalmmlmimmlmiCDMmmlmimmlmathinlineformula model fit to Planck primary CMB data.
['Cosmology', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
SPT clusters with DES and HST weak lensing. II. Cosmological constraints from the abundance of massive halos
{'SPT Clusters with DES and HST Weak Lensing. II. Cosmological Constraints from the Abundance of Massive Halos': "S. Bocquet, 1, ∗ S. Grandis, 2, 1 L. E. Bleem, 3, 4 M. Klein, 1 J. J. Mohr, 1, 5 T. Schrabback, 2, 6 T. M. C. Abbott, 7 P. A. R. Ade, 8 M. Aguena, 9 A. Alarcon, 3 S. Allam, 10 S. W. Allen, 11, 12, 13 O. Alves, 14 A. Amon, 15, 16 A. J. Anderson, 10 J. Annis, 10 B. Ansarinejad, 17 J. E. Austermann, 18, 19 S. Avila, 20 D. Bacon, 21 M. Bayliss, 22 J. A. Beall, 18 K. Bechtol, 23 M. R. Becker, 3 A. N. Bender, 3, 4, 24 B. A. Benson, 24, 4, 10 G. M. Bernstein, 25 S. Bhargava, 26 F. Bianchini, 11, 12, 13 M. Brodwin, 27 D. Brooks, 28 L. Bryant, 29 A. Campos, 30 R. E. A. Canning, 31 J. E. Carlstrom, 24, 4, 32, 3, 29 A. Carnero Rosell, 33, 9, 34 M. Carrasco Kind, 35, 36 J. Carretero, 20 F. J. Castander, 37, 38 R. Cawthon, 39 C. L. Chang, 4, 3, 24 C. Chang, 24, 4 P. Chaubal, 17 R. Chen, 40 H. C. Chiang, 41, 42 A. Choi, 43 T-L. Chou, 4, 32 R. Citron, 44 C. Corbett Moran, 45 J. Cordero, 46 M. Costanzi, 47, 48, 49 T. M. Crawford, 4, 24 A. T. Crites, 50 L. N. da Costa, 9 M. E. S. Pereira, 51 C. Davis, 11 T. M. Davis, 52 J. DeRose, 53 S. Desai, 54 T. de Haan, 55, 56 H. T. Diehl, 10 M. A. Dobbs, 41, 57 S. Dodelson, 30, 58 C. Doux, 25, 59 A. Drlica-Wagner, 24, 10, 4 K. Eckert, 25 J. Elvin-Poole, 60 S. Everett, 61 W. Everett, 62 I. Ferrero, 63 A. Fert'e, 13 A. M. Flores, 12, 11 J. Frieman, 10, 4 J. Gallicchio, 4, 64 J. Garc'ıa-Bellido, 65 M. Gatti, 25 E. M. George, 66 G. Giannini, 20, 4 M. D. Gladders, 24, 4 D. Gruen, 1 R. A. Gruendl, 35, 36 N. Gupta, 67 G. Gutierrez, 10 N. W. Halverson, 62, 19 I. Harrison, 68 W. G. Hartley, 69 K. Herner, 10 S. R. Hinton, 52 G. P. Holder, 36, 70, 57 D. L. Hollowood, 71 W. L. Holzapfel, 72 K. Honscheid, 73, 74 J. D. Hrubes, 44 N. Huang, 72 J. Hubmayr, 18 E. M. Huff, 61 D. Huterer, 14 K. D. Irwin, 13, 12 D. J. James, 75 M. Jarvis, 25 G. Khullar, 4, 24 K. Kim, 22 L. Knox, 76 R. Kraft, 75 E. Krause, 77 K. Kuehn, 78, 79 N. Kuropatkin, 10 F. K'eruzor'e, 3 O. Lahav, 28 A. T. Lee, 72, 80 P.-F. Leget, 81 D. Li, 18, 13 H. Lin, 10 A. Lowitz, 24 N. MacCrann, 82 G. Mahler, 83, 84 A. Mantz, 11, 12 J. L. Marshall, 85 J. McCullough, 81 M. McDonald, 86 J. J. McMahon, 4, 32, 24 J. Mena-Fern'andez, 87 F. Menanteau, 35, 36 S. S. Meyer, 4, 32, 24, 29 R. Miquel, 88, 20 J. Montgomery, 41 J. Myles, 89 T. Natoli, 24, 4 A. Navarro-Alsina, 90 J. P. Nibarger, 18 G. I. Noble, 91 V. Novosad, 92 R. L. C. Ogando, 93 Y. Omori, 4 S. Padin, 94 S. Pandey, 25 P. Paschos, 29 S. Patil, 17 A. Pieres, 9, 93 A. A. Plazas Malag'on, 81, 13 A. Porredon, 95 J. Prat, 24, 4 C. Pryke, 96 M. Raveri, 97 C. L. Reichardt, 17 J. Roberson, 22 R. P. Rollins, 46 C. Romero, 75 A. Roodman, 81, 13 J. E. Ruhl, 98 E. S. Rykoff, 81, 13 B. R. Saliwanchik, 99 L. Salvati, 100, 101, 102 C. S'anchez, 25 E. Sanchez, 103 D. Sanchez Cid, 103 A. Saro, 104, 102, 101, 105, 106 K. K. Schaffer, 4, 29, 107 L. F. Secco, 4 I. Sevilla-Noarbe, 103 K. Sharon, 108 E. Sheldon, 109 T. Shin, 110 C. Sievers, 44 G. Smecher, 41, 111 M. Smith, 112 T. Somboonpanyakul, 113 M. Sommer, 6 B. Stalder, 75 A. A. Stark, 75 J. Stephen, 29 V. Strazzullo, 101, 102 E. Suchyta, 114 G. Tarle, 14 C. To, 73 M. A. Troxel, 40 C. Tucker, 8 I. Tutusaus, 115 T. N. Varga, 116, 5, 117 T. Veach, 118 J. D. Vieira, 36, 70 A. Vikhlinin, 75 A. von der Linden, 110 G. Wang, 3 N. Weaverdyck, 14, 53 J. Weller, 5, 117 N. Whitehorn, 119 W. L. K. Wu, 13 B. Yanny, 10 V. Yefremenko, 3 B. Yin, 30 M. Young, 91 J. A. Zebrowski, 4, 24, 10 Y. Zhang, 7 H. Zohren, 6 and J. Zuntz 120  \n(the SPT and DES Collaborations)  \n1  \nUniversity Observatory, Faculty of Physics, Ludwig-Maximilians-Universitat, Scheinerstr. 1, 81679 Munich, Germany 2 Universitat Innsbruck, Institut fur Astro- und Teilchenphysik, Technikerstr. 25/8, 6020 Innsbruck, Austria 3 High-Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA 4 Kavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 5 Max Planck Institute for Extraterrestrial Physics, Gießenbachstr. 1, 85748 Garching, Germany 6 Argelander-Institut fur Astronomie, Auf dem Hugel 71, 53121 Bonn, Germany 7 Cerro Tololo Inter-American Observatory, NSF's National Optical-Infrared Astronomy Research Laboratory, Casilla 603, La Serena, Chile 8 School of Physics and Astronomy, Cardiff University, Cardiff CF24 3YB, United Kingdom 9 Laborat'orio Interinstitucional de e-Astronomia - LIneA, Rua Gal. Jos'e Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil 10 Fermi National Accelerator Laboratory, P. O. Box 500, Batavia, IL 60510, USA 11 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA 94305, USA 12 Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA 94305, USA 13 SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA 94025, USA 14 Department of Physics, University of Michigan, Ann Arbor, MI 48109, USA 15 Institute of Astronomy, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK 16 Kavli Institute for Cosmology, University of Cambridge, Madingley Road, Cambridge CB3 0HA, UK 17 School of Physics, University of Melbourne, Parkville, VIC 3010, Australia 18 NIST Quantum Devices Group, 325 Broadway Mailcode 817.03, Boulder, CO 80305, USA 19 Department of Physics, University of Colorado, Boulder, CO 80309, USA  \n56  \n20 Institut de F'ısica d'Altes Energies (IFAE), The Barcelona Institute of Science and Technology, Campus UAB, 08193 Bellaterra (Barcelona), Spain 21 Institute of Cosmology and Gravitation, University of Portsmouth, Portsmouth, PO1 3FX, UK 22 Department of Physics, University of Cincinnati, Cincinnati, OH 45221, USA 23 Physics Department, 2320 Chamberlin Hall, University of Wisconsin-Madison, 1150 University Avenue Madison, WI 53706-1390, USA 24 Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 25 Department of Physics and Astronomy, University of Pennsylvania, Philadelphia, PA 19104, USA 26 Department of Physics and Astronomy, Pevensey Building, University of Sussex, Brighton, BN1 9QH, UK 27 Department of Physics and Astronomy, University of Missouri, 5110 Rockhill Road, Kansas City, MO 64110, USA 28 Department of Physics & Astronomy, University College London, Gower Street, London, WC1E 6BT, UK 29 Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 30 Department of Physics, Carnegie Mellon University, Pittsburgh, Pennsylvania 15312, USA 31 Institute of Cosmology & Gravitation, University of Portsmouth, Dennis Sciama Building, Portsmouth, PO1 3FX, UK 32 Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 33 Instituto de Astrofisica de Canarias, E-38205 La Laguna, Tenerife, Spain 34 Universidad de La Laguna, Dpto. Astrof'ısica, E-38206 La Laguna, Tenerife, Spain 35 Center for Astrophysical Surveys, National Center for Supercomputing Applications, 1205 West Clark St., Urbana, IL 61801, USA 36 Department of Astronomy, University of Illinois Urbana-Champaign, 1002 West Green Street, Urbana, IL 61801, USA 37 Institut d'Estudis Espacials de Catalunya (IEEC), 08034 Barcelona, Spain 38 Institute of Space Sciences (ICE, CSIC), Campus UAB, Carrer de Can Magrans, s/n, 08193 Barcelona, Spain 39 Physics Department, William Jewell College, Liberty, MO, 64068 40 Department of Physics, Duke University Durham, NC 27708, USA 41 Department of Physics and McGill Space Institute, McGill University, 3600 Rue University, Montreal, Quebec H3A 2T8, Canada 42 School of Mathematics, Statistics & Computer Science, University of KwaZulu-Natal, Durban, South Africa 43 NASA Goddard Space Flight Center, 8800 Greenbelt Rd, Greenbelt, MD 20771, USA 44 University of Chicago, 5640 South Ellis Avenue, Chicago, IL 60637, USA 45 Jet Propulsion Laboratory, California Institute of Technology, Pasadena, CA 91011, USA 46 Jodrell Bank Center for Astrophysics, School of Physics and Astronomy, University of Manchester, Oxford Road, Manchester, M13 9PL, UK 47 Astronomy Unit, Department of Physics, University of Trieste, via Tiepolo 11, I-34131 Trieste, Italy 48 INAF-Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, I-34143 Trieste, Italy 49 Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy 50 Cornell University, Ithaca, NY 14853, USA 51 Hamburger Sternwarte, Universitat Hamburg, Gojenbergsweg 112, 21029 Hamburg, Germany 52 School of Mathematics and Physics, University of Queensland, Brisbane, QLD 4072, Australia 53 Lawrence Berkeley National Laboratory, 1 Cyclotron Road, Berkeley, CA 94720, USA 54 Department of Physics, IIT Hyderabad, Kandi, Telangana 502285, India 55 Institute of Particle and Nuclear Studies (IPNS), High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan International Center for Quantum-field Measurement Systems for Studies of the Universe and Particles (QUP), High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan 57 Canadian Institute for Advanced Research, CIFAR Program in Gravity and the Extreme Universe, Toronto, ON, M5G 1Z8, Canada 58 NSF AI Planning Institute for Physics of the Future, Carnegie Mellon University, Pittsburgh, PA 15213, USA 59 Universit'e Grenoble Alpes, CNRS, LPSC-IN2P3, 38000 Grenoble, France 60 Department of Physics and Astronomy, University of Waterloo, 200 University Ave W, Waterloo, ON N2L 3G1, Canada 61 Jet Propulsion Laboratory, California Institute of Technology, 4800 Oak Grove Dr., Pasadena, CA 91109, USA 62 Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO 80309, USA 63 Institute of Theoretical Astrophysics, University of Oslo. P.O. Box 1029 Blindern, NO-0315 Oslo, Norway 64  \nHarvey Mudd College, 301 Platt Boulevard, Claremont, CA 91711, USA  \n106  \n65 Instituto de Fisica Teorica UAM/CSIC, Universidad Autonoma de Madrid, 28049 Madrid, Spain 66 European Southern Observatory, Karl-Schwarzschild-Str.,DE-85748 Garching b. Munchen, Germany 67 CSIRO Space & Astronomy, PO Box 1130, Bentley WA 6102, Australia 68 School of Physics and Astronomy, Cardiff University, CF24 3AA, UK 69 Department of Astronomy, University of Geneva, ch. d' ' Ecogia 16, CH-1290 Versoix, Switzerland 70 Department of Physics, University of Illinois Urbana-Champaign, 1110 West Green Street, Urbana, IL 61801, USA 71 Santa Cruz Institute for Particle Physics, Santa Cruz, CA 95064, USA 72 Department of Physics, University of California, Berkeley, CA 94720, USA 73 Center for Cosmology and Astro-Particle Physics, The Ohio State University, Columbus, OH 43210, USA 74 Department of Physics, The Ohio State University, Columbus, OH 43210, USA 75 Center for Astrophysics | Harvard & Smithsonian, 60 Garden Street, Cambridge, MA 02138, USA 76 Department of Physics, University of California, One Shields Avenue, Davis, CA 95616, USA 77 Department of Astronomy/Steward Observatory, University of Arizona, 933 North Cherry Avenue, Tucson, AZ 85721-0065, USA 78 Australian Astronomical Optics, Macquarie University, North Ryde, NSW 2113, Australia 79 Lowell Observatory, 1400 Mars Hill Rd, Flagstaff, AZ 86001, USA 80 Physics Division, Lawrence Berkeley National Laboratory, Berkeley, CA 94720, USA 81 Kavli Institute for Particle Astrophysics & Cosmology, P. O. Box 2450, Stanford University, Stanford, CA 94305, USA 82 Department of Applied Mathematics and Theoretical Physics, University of Cambridge, Cambridge CB3 0WA, UK 83 Centre for Extragalactic Astronomy, Durham University, South Road, Durham DH1 3LE, UK 84 Institute for Computational Cosmology, Durham University, South Road, Durham DH1 3LE, UK 85 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 86 Kavli Institute for Astrophysics and Space Research, Massachusetts Institute of Technology, 77 Massachusetts Avenue, Cambridge, MA 02139, USA 87 LPSC Grenoble - 53, Avenue des Martyrs 38026 Grenoble, France 88 Instituci'o Catalana de Recerca i Estudis Avan¸cats, E-08010 Barcelona, Spain 89 Department of Astrophysical Sciences, Princeton University, Peyton Hall, Princeton, NJ 08544, USA 90 Instituto de F'ısica Gleb Wataghin, Universidade Estadual de Campinas, 13083-859, Campinas, SP, Brazil 91 David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada 92 Materials Sciences Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL 60439, USA 93 Observat'orio Nacional, Rua Gal. Jos'e Cristino 77, Rio de Janeiro, RJ - 20921-400, Brazil 94 California Institute of Technology, 1200 East California Boulevard., Pasadena, CA 91125, USA 95 Ruhr University Bochum, Faculty of Physics and Astronomy, Astronomical Institute, German Centre for Cosmological Lensing, 44780 Bochum, Germany 96 School of Physics and Astronomy, University of Minnesota, 116 Church Street SE Minneapolis, MN 55455, USA 97 Department of Physics, University of Genova and INFN, Via Dodecaneso 33, 16146, Genova, Italy 98 Department of Physics, Case Western Reserve University, Cleveland, OH 44106, USA 99 Brookhaven National Laboratory, Upton, NY 11973, USA 100 Universit'e Paris-Saclay, CNRS, Institut d'Astrophysique Spatiale, 91405, Orsay, France 101 INAF - Osservatorio Astronomico di Trieste, via G. B. Tiepolo 11, 34143 Trieste, Italy 102 IFPU - Institute for Fundamental Physics of the Universe, Via Beirut 2, 34014 Trieste, Italy 103 Centro de Investigaciones Energ'eticas, Medioambientales y Tecnol'ogicas (CIEMAT), Madrid, Spain 104 Astronomy Unit, Department of Physics, University of Trieste, via Tiepolo 11, 34131 Trieste, Italy 105 INFN - National Institute for Nuclear Physics, Via Valerio 2, I-34127 Trieste, Italy ICSC - Italian Research Center on High Performance Computing, Big Data and Quantum Computing, Italy 107 Liberal Arts Department, School of the Art Institute of Chicago, 112 South Michigan Avenue, Chicago, IL 60603, USA 108 Department of Astronomy, University of Michigan, 1085 S. University Ave, Ann Arbor, MI 48109, USA 109 Brookhaven National Laboratory, Bldg 510, Upton, NY 11973, USA 110 Department of Physics and Astronomy, Stony Brook University, Stony Brook, NY 11794, USA 111 Three-Speed Logic, Inc., Victoria, B.C., V8S 3Z5, Canada 112 School of Physics and Astronomy, University of Southampton, Southampton, SO17 1BJ, UK 113 Department of Physics, Faculty of Science, Chulalongkorn University, 254 Phayathai Road, Pathumwan, Bangkok 10330, Thailand  \n114 Computer Science and Mathematics Division, Oak Ridge National Laboratory, Oak Ridge, TN 37831, USA  \n115 Institut de Recherche en Astrophysique et Plan'etologie (IRAP), Universit'e de Toulouse, CNRS, UPS, CNES, 14 Av. Edouard Belin, 31400 Toulouse, France 116 Excellence Cluster Origins, Boltzmannstr. 2, 85748 Garching, Germany 117 Universitats-Sternwarte, Fakultat fur Physik, Ludwig-Maximilians Universitat Munchen, Scheinerstr. 1, 81679 Munchen, Germany 118 Space Science and Engineering Division, Southwest Research Institute, San Antonio, TX 78238, USA 119 Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA 120 Institute for Astronomy, University of Edinburgh, Edinburgh EH9 3HJ, UK (Phys. Rev. D accepted 24 May 2024)  \nWe present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range 0 . 25 -1 . 78 over a total sky area of 5,200 deg 2 . We use DES Year 3 weak-lensing data for 688 clusters with redshifts z < 0 . 95 and HST weak-lensing data for 39 clusters with 0 . 6 < z < 1 . 7. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable-mass relation without having to make strong assumptions about, e.g., the hydrodynamical state of the clusters. For a flat ΛCDM cosmology, and marginalizing over the sum of massive neutrinos, we measure Ω m = 0 . 286 ± 0 . 032, σ 8 = 0 . 817 ± 0 . 026, and the parameter combination σ 8 (Ω m / 0 . 3) 0 . 25 = 0 . 805 ± 0 . 016. Our measurement of S 8 ≡ σ 8 √ Ω m / 0 . 3 = 0 . 795 ± 0 . 029 and the constraint from Planck CMB anisotropies (2018 TT,TE,EE+lowE) differ by 1 . 1 σ . In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses ∑ m ν < 0 . 18 eV. When additionally allowing the dark energy equation of state parameter w to vary, we obtain w = -1 . 45 ± 0 . 31 from our cluster-based analysis. In combination with Planck data, we measure w = -1 . 34 +0 . 22 -0 . 15 , or a 2 . 2 σ difference with a cosmological constant. We use the cluster abundance to measure σ 8 in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the ΛCDM model fit to Planck primary CMB data.", 'I. INTRODUCTION': "Understanding the cause for the accelerated expansion of our universe is currently one of the biggest challenges in physics. To address this robustly, a multitude of different probes and techniques have been proposed and are being used to observationally constrain the cosmological parameters. Interestingly, the comparison of results obtained from these probes have recently uncovered two anomalies: First, the present-day value of the Hubble parameter determined from primary anisotropies of the cosmic microwave background (CMB) at high redshift is in at least 4 σ tension with measurements in the local universe using the distance ladder e.g., [1-4]. Second, measurements of the S 8 -parameter 1 by low-redshift probes have tended to be lower than predicted by the ΛCDM model with parameters set by high-redshift CMB data e.g., [5-13].  \nOver recent years, new large-scale structure probes such as joint analyses of galaxy clustering and weaklensing 2-point correlation functions (so-called 3 × 2 pt analyses) [8, 9, 14], the CMB lensing power spectrum [1517], or the combination of all of these [18] have emerged. Measurements of the abundance of massive dark mat-  \ner halos (and the galaxy clusters they host) probe the matter density field through a different mechanism and are subject to different systematics, thus offering an ideal and powerful complement [19, 20]. However, for a study of halo abundance to be successful, the relationship between the underlying, un-observable halo mass and the observable properties of the clusters needs to be known (see e.g., reviews [21, 22]).  \nIn this work, we present cosmological constraints derived from a sample of 1,005 galaxy clusters detected using the South Pole Telescope (SPT) [23] in combination with galaxy and weak-lensing data from the Dark Energy Survey (DES) [24-26], photometric data from the Wide-field Infrared Survey Explorer (WISE) [27], and targeted weak-lensing measurements from the Hubble Space Telescope (HST). We select the galaxy clusters via the thermal Sunyaev-Zel'dovich (hereafter SZ) [28] effect in data from the first two SPT cameras, SPT-SZ and SPTpol (SPT-3G being the third, currently operating camera). The SZ effect is sourced by the electron pressure in the hot gas component of the cluster (the intra-cluster medium, or ICM) and arises when CMB photons scatter with these high-energy electrons, leading to a distinct spectral signature. The SZ effect signal is thus a clean tracer of the hot ICM and enables the detection of galaxy clusters out to the highest redshifts at which massive halos exist, thereby creating an essentially mass-limited sample whose limiting mass does not evolve strongly with redshift. Optical and near-infrared data from DES, WISE, and other targeted observation pro-  \ngrams are then used to confirm clusters (and to exclude false detections) and to assign redshifts to all confirmed systems. With these follow-up data taken from widefield surveys, we can robustly calibrate the probability of chance association of random optical structures with an SZ noise fluctuations through measurements along random lines of sight [29, 30].  \nTo relate the strength of the SZ effect signal to the underlying halo mass, we cannot rely on first principles because modeling the dynamical state of the ICM is highly complex. In particular, the ICM cannot be assumed to be in hydrostatic equilibrium. Instead, we use measurements of the weak gravitational shear, which depends on the total halo mass irrespective of its dynamical state, and which can thus be mapped to halo mass with excellent control over systematic uncertainties (for a review on cluster lensing, see [31]). We now leverage the strategic overlap of the SPT survey and DES footprints andfor the first time-use DES lensing data for mass calibration in a cosmological analysis of SPT-selected clusters. Within our modeling framework, we empirically constrain the SZ-mass relation in a robust, weak-lensingbased way, thereby addressing and overcoming the main challenge in cluster cosmology: the accurate calibration of cluster masses.  \nUsing 5,200 deg 2 of SPT cluster surveys, we select 1,005 clusters at redshift z > 0 . 25. In the 3,567 deg 2 overlap region between the DES and SPT surveys, we obtain measurements of weak-lensing shear from DES Year 3 data (hereafter DES Y3) for 688 clusters. Since DES lensing runs out of constraining power for highredshift ( z ≳ 0 . 9) clusters, we supplement the DES lensing dataset with targeted lensing observations of 39 clusters in the redshift range 0 . 6 < z < 1 . 7 using the HST. This joint cluster and lensing dataset represents a significant improvement over the previous analysis of SPT clusters [7], which was based on 344 clusters in the SPT-SZ survey [32]. In that analysis, we used cluster lensing data from targeted programs at Magellan (19 clusters, [33]) and HST (13 clusters, [34]). The reduced statistical and systematic uncertainties of the dataset we consider here enables significantly improved cosmological constraints.  \nIn this paper, we present cosmological constraints from the joint analysis of the SPT cluster abundance with a weak-lensing mass calibration using DES and HST data. The data products, the cluster lensing measurements using DES Y3 data, the modeling choices, the analysis pipeline, and its validation are presented in detail in a companion paper ([35], hereafter Paper I). We briefly review the dataset in Sec. II. We summarize the analysis method in Sec. III. We discuss our blinding strategy and the robustness tests that were performed prior to unblinding of the results in Sec. IV. We present our cosmological constraints in Sec. V, and close with a summary in Sec. VI.  \nThroughout this analysis, we assume spatial flatness. The parameter combination S 8 ≡ σ 8 (Ω m / 0 . 3) 0 . 5 is defined to be tightly constrained by cosmic shear. We  \nempirically determine that our cluster dataset optimally constrains a parameter combination with an exponent of 0.225. However, for comparability, we adopt S opt 8 ≡ σ 8 (Ω m / 0 . 3) 0 . 25 , which is the optimal combination for CMB lensing [17]. 2 Halo masses M 200c refer to the mass enclosed within a sphere of radius r 200c , within which the mean density is 200 times larger than the critical density ρ c ( z ) at the cluster redshift z . We express the (multivariate) normal distribution with mean µ and (co)variance K as N ( µ , K ).", 'II. DATA': 'We summarize the data products used in this analysis. For a detailed description, we refer the reader to Paper I.', 'A. SPT Cluster Catalog': 'We use the cluster catalogs from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys [30, 32, 36, 37]. Combined, they cover over 5,200 deg 2 of the southern sky. Note that the deep SPTpol 500d survey lies within the footprint of the SPT-SZ survey; in the overlap area, we use only data from SPTpol 500d. The cluster candidate list is selected on the basis of SPT detection significance ξ and redshift z .  \nOver the part of the SPT survey that is not covered by DES (27% of the total SPT-SZ and SPTpol cluster survey area), we confirm clusters and assign redshifts using targeted optical observations (using among others, the PISCO imager [38]) as described in [30, 32]. We select all candidates with  \nξ > 5 , z > 0 . 25 , (1)  \nwhich, after follow up, results in a highly pure cluster sample ( ≳ 95%) [30, 32].  \nOver the 3,567 deg 2 of the SPT survey that is covered by DES (73% of the total survey area), we obtain measurements of optical richness λ and of redshift z using the multi-component matched filter cluster confirmation tool (MCMF) [29, 36]. At high cluster redshifts z > 1 . 1, which are beyond the reach of DES, we use data from WISE to compute richness and redshift. To characterize the probability of chance associations, we also compute richness and redshift for random lines of sight in the DES footprint. With this information, we can define a minimum richness λ min above which we consider an SZ cluster candidate as confirmed. In practice, λ min depends on redshift, and we define λ min ( z ) such that the  \n<!-- image -->  \n<!-- image -->  \nFIG. 1. Projection of the 1,005 clusters in the SPT sample along the dimensions of redshift, SPT detection significance, and optical richness. Orange histograms show the data. Thick lines show the mean recovered model in an analysis of the cluster abundance and cluster lensing, thin solid lines show the 1 σ uncertainty on the mean, and dotted lines show the 1 σ Poisson uncertainty (statistical shot noise). There is qualitative agreement between the model and the data.  \n<!-- image -->  \nfinal sample purity is > 98% (see Fig. 2 in Paper I). 3 At fixed SPT detection significance, the different SPT surveys have different levels of purity because of the varying depths. To keep the overall sample purity approximately constant, over the joint SPT-DES area, we select clusters according to  \nξ > 4 . 25 / 4 . 5 / 5 (500d / SZ / ECS) , λ > λ min ( z ) , z > 0 . 25 . (2)  \nThe final sample comprises 1,005 confirmed clusters with redshift measurements. 4 All clusters in the DES region also have richness measurements and optically determined center positions (from DES or WISE). Figure 1 shows the distribution of clusters as a function of redshift, SPT detection significance, and optical richness.', 'B. DES Y3 Weak-Lensing Data': 'The DES covers 5,000 deg 2 of the southern sky in the g , r , i , z , and Y bands. The DES Y3 cosmology dataset covers 4,143 deg 2 after masking, of which 3,567 deg 2 overlap with the SPT cluster surveys. The weak-lensing shape catalog [39] is created by applying the Metacalibration pipeline [40, 41] to data from the r , i , and z bands. Detailed information about the photometric dataset [42], the modeling of the point-spread function [43], and image and survey simulations [44, 45] can be found in dedicated  \nDES Y3 publications. Following the 3 × 2 pt analysis [9], we select lensing source galaxies in four tomographic bins. The redshift distributions of these source bins are calibrated using self-organizing maps [46].  \nFor every SPT cluster in the DES footprint, we extract a weak-lensing shear profile within the radial range 0 . 5 < r/ ( h -1 Mpc) < 3 . 2 (1 + z cluster ) -1 , centered on the optically determined cluster center. 5 This radial cut avoids the problematic central region of the cluster (affected by feedback from active galactic nuclei, miscentering, blending, cluster member contamination, non-linear shear) and ensures that only the 1-halo term regime is considered [47]. We restrict the use of DES weak-lensing data to clusters below redshift z = 0 . 95, which corresponds to the median redshift of the highest-redshift source bin. We extract 688 cluster shear profiles from a total of 555,912 source galaxies. For illustrative purposes, we show stacked shear profiles in bins of cluster redshift and SPT detection significance in Fig. 2.  \nFor the purpose of robustness checks (see Sec. IV), we also measure shear profiles in a more conservative radial range 0 . 8 < r/ ( h -1 Mpc) < 3 . 2 (1 + z cluster ) -1 . As another alternative, we measure shear around the cluster centers as determined in the SZ analysis. We estimate the level of cluster member contamination using two different redshift estimators ( DNF [48] as the default, and BPZ [49] as an alternative).  \nWe quantify and discuss the relevant sources of systematic and statistical uncertainties related to the weaklensing measurements (cluster member contamination, miscentering of the shear profile, shear and photoz cali-  \nFIG. 2. Stacked DES Y3 lensing shear profiles. Upper panel: Tangential shear, along with the mean model prediction from an analysis of the cluster abundance and cluster lensing. The likelihood analysis is performed on a cluster-by-cluster basis, the stacks are only shown for validation purposes. The dark shading shows the radial ranges used in the analysis, light shaded data points outside that range are only shown for reference. For the 26 data points included in the analysis, we measure a signal-to-noise ratio of 32.2 and we obtain a reduced χ 2 red = 1 . 48 between the data and the model. Lower panel: Within the radial ranges used in the analysis, the cross shear is consistent with null with χ 2 red = 0 . 90.  \n<!-- image -->  \nbration, halo mass modeling, impact of large-scale structure) in detail in Paper I, Sec. V.', 'C. High-Redshift HST Weak-Lensing Data': 'The DES lensing data does not have substantial constraining power for lenses with z ≳ 0 . 9, and we use HST data to enable weak-lensing mass calibration at higher redshifts. Using the HST-39 dataset presented and analyzed in [34, 50-53], we have a sample of 39 clusters in the redshift range 0 . 6 -1 . 7 with space-based weak-lensing measurements.', 'III. ANALYSIS METHOD': 'The analysis method is introduced and described in detail in Paper I. Here, we briefly review the key features.', 'A. Summary of Analysis Strategy': "We set up a Bayesian (hierarchical) population model to describe our dataset (see also [54-56]). We explicitly model the selection criteria that we adopted to construct the cluster catalog. Specifically, these are cuts in the SPT detection significance ξ , cuts in the optical richness λ (over the footprint that is common with DES), and a cut in redshift [see Eqs. (1) and (2)]. Clusters with weak-lensing data additionally have measurements of tangential shear profiles.  \nWe follow analytical and simulation-based models for observable-mass relations e.g., [57, 58], and describe the scaling relation between the unbiased SPT detection significance ζ and mass as  \n⟨ ln ζ ⟩ =ln A SZ + B SZ ln ( M 200c 3 × 10 14 h -1 M ⊙ ) + C SZ ln ( E ( z ) E (0 . 6) ) (3)  \nwith lognormal intrinsic scatter of width σ ln ζ , and relate the detection significance ξ to ζ as  \nP ( ξ | ζ ) = N ( √ ζ 2 +3 , 1 ) (4)  \nto account for the maximization bias when measuring ξ in noisy data [59]. We account for the varying depth within the SPT survey by rescaling A SZ and C SZ for each individual field [30, 32, 37]. For the SPTpol ECS survey, the overall normalization γ ECS of A SZ is difficult to calibrate (owing to detector linearity changes under high atmospheric loading at the elevation of these fields, see Sec. 2.2 [30]), and we thus consider γ ECS as an additional fit parameter.  \nWe choose the same ansatz for the relation between the intrinsic richness ˜ λ and mass  \n⟨ ln ˜ λ ⟩ =ln A λ + B λ ln ( M 200c 3 × 10 14 h -1 M ⊙ ) + C λ ln ( 1 + z 1 . 6 ) , (5)  \nwith lognormal intrinsic scatter of width σ ln ˜ λ , and relate the observed richness λ to ˜ λ as  \nP (ln λ | ln ˜ λ ) = N (ln ˜ λ, 1 / ˜ λ ) . (6)  \nWe consider two different types of measured richness; those based on DES data, which we use for clusters with redshift z < 1 . 1, and those based on WISE, which we use for clusters with redshifts beyond 1.1. Rather than attempting to match the two types of richness within an overlap range in redshift, we fit a separate observablemass relation for each of the two types of richness, with a clean transition from DES-based to WISE-based richness at redshift z = 1 . 1.  \nUsing the relations defined in Eqs. (3)-(6), we can model the ξ -λ -z sample by convolving the halo mass  \nfunction with the observable-mass relations (accounting for the covariance in intrinsic scatter ρ SZ , ˜ λ σ ln ζ σ ln ˜ λ ) and applying the sample selection cuts. There is no external information on the parameters A,B,C of the scaling relations or on the intrinsic scatter σ ln ζ , σ ln ˜ λ , and we use weak-lensing data of sample clusters to calibrate these observable-mass relation parameters empirically. Because we can infer the halo mass from lensing data with excellent control over systematic uncertainties, this is a robust analysis strategy.  \nWe introduce the modeling framework with which we analyze the DES Y3 lensing data in Paper I. Here, we briefly review the method. We adopt a simple model for the projected halo mass distribution Σ based on a modified Navarro-Frenk-White profile (NFW) [60]. 6 We introduce a latent variable that we call the 'weak-lensing mass' M WL , which we define such that the reduced tangential shear is  \ng t ( r, M WL ) = ∆Σ( r, M WL ) Σ -1 crit 1 -Σ( r, M WL ) Σ -1 crit ( 1 -f cl ( r ) ) . (7)  \nThe density contrast ∆Σ( r ) ≡ ⟨ Σ( < r ) ⟩ -Σ( r ) is computed from Σ, and the lensing efficiency Σ -1 crit is computed from the distribution of source redshifts and the lens redshift. Note that Σ -1 crit explicitly depends on cosmology. The cluster member contamination f cl ( r ) accounts for contaminants in the lensing source sample that are not sheared and thus bias the measurement low [therefore, (1 -f cl ) -1 is often referred to as the boost factor].  \nThe cluster lensing model we just introduced is not perfect (e.g., halos do not individually match the NFW profile) and M WL inferred from the shear profile g t is thus a biased and noisy estimator of the true halo mass M halo [61, 62]. To account for this, we establish a mean relationship,  \n〈 ln ( M WL M 0 )〉 = b WL ( z ) + b WL ,M ln ( M 200c M 0 ) , (8)  \nand describe the width of the lognormal scatter around the mean as  \nln σ ln M WL = 1 2 [ s WL ( z ) + s WL ,M ln ( M 200c M 0 )] , (9)  \nwith a pivot mass M 0 = 2 × 10 14 h -1 M ⊙ .  \nWe create synthetic cluster shear maps by applying the source redshift distribution, cluster miscentering, and cluster member contamination to halo mass maps from numerical simulations, and we then use these maps to calibrate the free parameters of the M WL -M halo relation  \n[47]. 7 Note that we use the mass maps from full-physics hydrodynamical simulations and establish the relationship between M WL and the 'gravity-only mass' M halo as measured in paired gravity-only simulations with identical initial conditions. This has the benefit that our model includes the effect of baryons, but we are still able to use the arguably more robust halo mass function predictions from gravity-only simulations [such as 63] (including emulators e.g., [64-66]). In practice, we apply this method to the Magneticum [67-70] and the Illustris TNG simulations [71-76]. The recovered constraints on the M WL -M halo relation differ somewhat, which we interpret as an uncertainty in the modeling of baryonic effects. We inflate the uncertainty on all parameters of the M WL -M halo relation by this additional uncertainty given in Table 2 [47] that, e.g., for the amplitude b WL , amounts to 2%. As we will show below, this level of systematic uncertainty is smaller than the current level of statistical uncertainty in the lensing dataset, and our analysis therefore does not strongly depend on the assumption that the Magneticum or Illustris TNG simulation correctly reproduces the real universe.  \nFor the HST-39 dataset, we use a similar modeling approach, and also adopt M WL as a latent variable. Note, however, that several analysis choices differ. For example, the selection of HST lensing sources via color cuts leads to a much smaller level of cluster member contamination than in the DES Y3 analysis. Therefore, we do not explicitly model this contamination, and we do not explicitly model the effect of miscentering, either. We do, of course, take all sources of uncertainty into account in the M WL -M halo . The full HST cluster lensing model is presented in detail in the original analyses [34, 52, 53].", 'B. Likelihood, Pipeline, Priors, and Sampling': "We infer parameters p from our data assuming a cluster population model. Bayes' theorem states  \nP ( p | data , model) ∝ P (data | p , model) P ( p | model) , (10)  \nwhere P ( p | model) is the prior probability distribution of the parameters. The factor P (data | p , model) is the likelihood L of the data given the (fixed) model and the parameters. A detailed discussion of our likelihood function and its numerical implementation can be found in Paper I, Sec. VII. Briefly, we adopt a log-likelihood  \nln L =lnPoisson ( sample { ξ i , λ i , z i } N clusters i =1 ∣ ∣ ∣ ∣ d 3 N ( p ) d ξ d λ d z ) + ∑ i ln P ( λ i , g t ,i | λ i > λ min ( z ) , ξ i , z i , p ) ,  \n(11)  \nwhere the first term is the unbinned Poisson likelihood of the cluster sample { ξ, λ, z } and its distribution d 3 N ( p ) d ξ d λ d z in ξ -λ -z space, and the second term is the 'mass calibration likelihood': the conditional probability to observe the tangential shear profile g t and richness λ for a given cluster in the sample with measured ξ and z , and whose richness is above the minimum richness λ min .  \nWe implement the likelihood function as a Python module in the CosmoSIS framework [77]. 8 We validate the implementation of the analysis pipeline by analyzing several statistically independent mock catalogs that we create from the model. Because creating those mocks is less challenging than analyzing them, this is a meaningful test of the analysis pipeline (rather than a test of the mock creation). For further details, we refer the reader to Paper I, Sec. VIII.  \nThe parameters that we vary in the baseline analysis are defined in Table I, along with the informative priors where applied. The model for the DES cluster lensing mass bias and scatter are functions of cluster redshift, and the parameters σ ln b WL , 1 / 2 and s WL modulate the uncertainty of the model (therefore, the priors are centered on 0 with unit width) [see Paper I, Sec. V B]. We follow the same approach for the HST lensing data, where each cluster has its own bias and scatter, and the shared uncertainties are modulated via b HST and s HST . The priors on the scatter ( > 0 . 05) improve the numerical stability of the code. The priors on the correlation coefficients U ( -0 . 5 , 0 . 5) are chosen such that the resulting correlation matrix is nonsingular for all parameter combinations. The parameters h , Ω b h 2 , and n s are not constrained by the cluster dataset, and are not significantly correlated with the parameters of interest. We apply Gaussian priors on Ω b h 2 and n s from Planck [2], and a prior on the Hubble parameter h ∼ N (0 . 70 , 0 . 05 2 ) that is wide enough to encompass all currently favored constraints. When analyzing the cluster dataset jointly with Planck , we do not apply these three Gaussian priors, and we additionally consider the optical depth to reionization τ as a free parameter.  \nWe explore the parameter space using either the multinest or the nautilus sampler [78, 79]. The potential amount by which these samplers may or may not underestimate the parameter uncertainties is small and does not change our findings in any qualitatively relevant way [80].", 'IV. BLINDING, ROBUSTNESS TESTS, AND UNBLINDING': 'In this section, we describe our blinding strategy and the robustness tests that were performed prior to unblinding of the results. We check for internal consistency  \nTABLE I. Fit parameters in the baseline analysis of the abundance of SPT clusters with DES Y3 and HST-39 lensing. A missing entry (-) for the prior indicates that a uniform prior is applied that is wide enough to not be informative. Priors with a lower limit ( > ) are uniform above that bound with an un-informative upper bound. The parameters of the lensing models are prior-dominated and inform the empirical calibration of the other observable-mass relations. The prior on Ω ν h 2 corresponds to a prior on the sum of neutrino masses ∑ m ν ∼ U (0 , 0 . 6) eV.  \nof the dataset and verify the robustness to several analysis choices in the processing of the DES Y3 lensing data.', 'A. Blinding Strategy': 'To avoid compromising our analysis with any expectations we may have regarding what the results should or should not look like, we proceed as follows. First, the SPT cluster catalogs and the DES Y3 lensing data were finalized before we started this analysis. Second, we did not compare any weak-lensing mass estimates or preliminary SZ-mass relations with literature values. Third, we blinded the astrophysical and cosmological constraints obtained in this work until the robustness of the results against specific analysis choices was demonstrated. To do so, we created unknown, random, but constant blinding offset parameters with which we shifted the parameters when reading in any MCMC chain. It was thus possible to compare different blinded analyses with each other to check for relative shifts in the recovered parameters (or changes in the parameter uncertainties) without knowing the actual parameter values.  \nThe HST-39 cluster lensing dataset was finalized, and its analysis published before we started this analysis. Therefore, the parameter constraints favored by this dataset were known, and we preferred to not use these data in the blinded phase of the analysis.  \nThe blinded runs were set up assuming a ΛCDM model with wide flat priors on the fit parameters. Failure to choose sufficiently wide priors might lead to the blinded results showing hard cuts, thereby effectively unblinding the respective parameter constraint. Note that arbitrarily wide priors are not always possible, because, e.g., the scatter parameters or Ω m need to be strictly positive. All but the first test are summarized in Fig. 3, which shows the constraints on key parameters. The results in the first row (blue) are obtained from the full dataset (without HST-39), marginalized over all parameters in Table I. The second row (black) shows the constraints from the same analysis, except that we fix a selection of model parameters to their fiducial values: the parameters of the M WL -M halo relations, the correlation coefficients on intrinsic scatter, Ω b h 2 , Ω ν h 2 , and n s . The parameter shifts are small, which also suggests that our constraints are not limited by the systematic uncertainty in the DES Y3 lensing model. To save computational resources, and to make the robustness tests somewhat more stringent, all cross-check analyses discussed in the next section are run in this restricted parameter space, with a selection of parameters fixed as just discussed.  \nBefore running the blinded robustness tests, we had agreed that we would unblind the results if no test reveals a statistically significant difference with the baseline analysis.', 'B. Tests Performed prior to Unblinding': 'We now summarize the tests that were conducted with the blinding setup discussed above.', '1. Combination of Mass Calibration and Abundance': "We derive constraints from just the cluster abundance [first term in the likelihood Eq. (11)] and just mass calibration (second term in the equation). We do not expect interesting constraints on the cosmological parameters from these individual analyses, and Fig. 4 shows broad uncertainties in the cosmological parameters and elongated parameter degeneracies (gray and red contours). 9 The two posterior distributions overlap, and we combine the two analyses to obtain constraints from the cluster abundance with weak-lensing-informed mass calibration (blue contours in Fig. 4). The constraints on select parameters are shown labeled as 'full sample' in Fig. 3. We note that the situation is similar for the parameters of the richness-mass relations, for which the individual constraints overlap, and tighten when we combine the abundance with the lensing mass calibration.", '2. Goodness of Fit': 'We validate that the model is an adequate description of the data by extracting posterior predictive distributions from the blinded MCMC chains (of the cluster abundance and lensing analysis) and comparing them with the data. In Fig. 1, blue lines show the mean predicted model for the cluster abundance. Similarly, in Fig. 2, we show the mean predicted stacked shear profiles and the error on the mean. We compute the statistical difference between each stack and the mean model, and obtain values for χ 2 ∈ { 11 . 3 , 8 . 1 , 7 . 3 , 10 . 6 } for { 7 , 7 , 6 , 6 } data points in the { ( z < 0 . 5 , ξ ≥ 5 . 6) , ( z < 0 . 5 , ξ < 5 . 6) , ( z ≥ 0 . 5 , ξ ≥ 5 . 6) , ( z ≥ 0 . 5 , ξ < 5 . 6) } bins. For the full dataset, we measure χ 2 = 37 . 4 for 26 data points, and conclude that the model adequately describes the data.', '3. Uncertainty in the Halo Mass Function': "We marginalize over a 3.5% uncertainty in the amplitude and a 2% uncertainty in the slope of the mass  \nFIG. 3. Impact of different analysis choices on the constraints on key cosmological and SZ-mass relation parameters (mean and 68% credible interval), assuming a ΛCDM model. During the blinded phase of the analysis, the parameter values were artificially offset by an unknown amount. The vertical gray bands show the 68% credible interval of the 'full sample' analysis (second row) for reference. All results shown here are based on mass calibration using only the DES Y3 lensing data, the HST-39 lensing data are not used. Only the results in the first row shows constraints that are marginalized over all nuisance parameters; for simplicity, for all other tests, we fix the parameters of the M WL -M halo relations, the intrinsic scatter correlation coefficients, and Ω ν h 2 , Ω b h 2 , and n s to their nominal values. DNF and BPZ are two redshift estimators that we use to determine the amount of cluster member contamination in the lensing signal. MCMF and SPT refer to cluster centers as measured in optical data by MCMF or as determined in the SZ analysis of SPT data.  \n<!-- image -->  \n8  \nfunction (following [81]). The lack of parameter shifts in Fig. 3 suggests that we can neglect uncertainties in the halo mass function, which we will do from here on out. 10", '4. Restricting the Analysis to Low Redshift z < 1': 'At redshifts beyond z ≳ 1, it remains unclear how strong the correlated astrophysical emission, and thus the contamination of the cluster SZ signal is, although observational studies favor an effect that is only at the few-percent level e.g., [82-85]. Also recall that we use DES data for cluster confirmation up to redshift z = 1 . 1, beyond which we use data from WISE. When limiting the analysis to z < 1, we thus test the robustness of using WISE data for cluster confirmation and redshift assignment jointly with the robustness against astrophysical contamination of the SZ signal. We caution that disagreement between the fiducial sample and the lowredshift subset could also indicate interesting new physics  \n(e.g., different growth of structure beyond z ≳ 1) rather than pointing to the effects we might worry about.  \nThe constraints from the 0 . 25 < z cluster < 1 dataset (869 clusters compared to 1,005 in the full sample) show negligible shifts compared to the analysis of the full sample (see Fig. 3). Since the remaining robustness checks are not expected to be sensitive to the redshift range of the sample, we perform the following tests on the 0 . 25 < z cluster < 1 sample.', '5. The ξ -ζ Relation': "The ξ -ζ relation [Eq. (4)] cannot hold down to arbitrarily low significance as the formalism breaks down for ξ < √ 3. As in previous work, we define a limit ζ min , below which the maximization bias in the ξ -ζ relation does not apply. Using simulations of the SZ sky, we estimate ζ min = 1. Here, we confirm that raising the threshold to ζ min = 2 (which was applied in, e.g., [86]) only has a small impact on the recovered parameters (' ζ min : 1 → 2' in Fig. 3).  \nFIG. 4. Constraints (68% and 95% credible regions) on select scaling relation and cosmology parameters, as obtained from just the DES Y3 weak-lensing-based mass calibration (gray) or just the cluster abundance (red). The abundance and mass calibration results do not show signs of significant disagreement, and we combine them to obtain cosmological constraints from the joint analysis of cluster abundance and lensing (blue).  \n<!-- image -->", '6. Selection Cut in SPT Detection Significance ξ': "Enabled by the DES data and the MCMF analysis, in this study, we select clusters down to lower detection significances than previously done [see Eq. (2)]. We compare the constraints obtained using this selection scheme  \n(869 clusters with z < 1) to the constraints obtained using the more conservative ξ > 5 selection (654 clusters with z < 1). As shown in Fig. 3 as ' ξ > 5' in purple, the constraint on S opt 8 is robust to this change, but the SZ mass slope B SZ and the redshift evolution C SZ both shift to higher values, and, due to parameter degeneracies, Ω m also shifts to larger values.  \nWe quantify these shifts more carefully. Since S opt 8 is effectively unchanged, and because of the parameter degeneracies, it is sufficient to focus on the constraints on Ω m . First, we consider the analysis of the 'low-redshift' 0 . 25 < z cluster < 1 dataset (green in Fig. 3), and examine the impact of additionally applying ξ > 5. Naively, one would compare the difference in the mean recovered values for Ω m with the expected variance, computed as the sum of the variances from each analysis. Doing so, we observe a shift in Ω m that is at the 0 . 5 σ level. However, this simple estimate neglects the fact that the ξ > 5 sample is a subset of the larger catalog. In this case, 'the covariance of the parameter differences between the partial and full analyses is simply the difference of the covariances' [87]. With this estimator, we report that Ω m shifts at the 0 . 9 σ level.  \nTo empirically validate that this level of parameter shift is expected, we analyze pairs of mock catalogs (restricted to z < 1) for which one catalog is additionally restricted to ξ > 5. 11 After running two such mock analyses, we do indeed observe parameter shifts that are very similar to the one obtained from the real data.  \nFinally, we compute the consistency of the Ω m measurement between the full cluster sample ( z > 0 . 25) and the ξ > 5 , 0 . 25 < z < 1 subsample (purple vs. black in Fig. 3). Naively adding the variances gives a shift at the 0 . 9 σ level, and correctly subtracting the variances yields 1 . 2 σ . We caution the reader that in this comparison, we consider the combined effects of restricting ourselves to low redshifts (which shifts Ω m toward higher values) and to high significances (which further shifts Ω m high). To summarize, we find that, while the observed shifts are interesting and might reveal new properties of SZ-selected clusters in upcoming analyses of larger samples, there is no statistically significant evidence of internal tension in the dataset.", '7. Cluster Member Contamination': "In the DES Y3 lensing analysis, we account for cluster member contamination of the lensing source catalog as determined by applying the P ( z ) decomposition method to DNF source redshifts (see Paper I, Sec. IV D). However, the level of contamination derived using BPZ redshifts is significantly different. We test the robustness of our cosmological results to adopting the cluster member contamination as derived using BPZ instead of the fiducial one derived using DNF. The resulting parameter shifts are ≲ 1 σ (see 'DES WL boost: DNF → BPZ' in Fig. 3). While this level of shifts is acceptable given the overall level of uncertainty, this test highlights the relative importance of the cluster member contamination  \nin our analysis. Future lensing surveys are expected to enable a cleaner source selection, e.g., [88], which would reduce that relative importance.", '8. Centers Adopted for the Lensing Analysis': "Our default analysis uses lensing data extracted around cluster positions as determined using the DES imaging data. We constrain the optical miscentering distribution and account for it in the lensing analysis (see Paper I, Secs. IV C and V B). As an alternative, we use the cluster centers as determined in the SZ analysis, and account for the corresponding miscentering distribution in the lensing model. The recovered constraints on A SZ and S opt 8 shift by ≲ 1 σ , while the other parameters remain essentially unchanged ('DES WL center: MCMF → SPT' in Fig. 3). We conclude that our miscentering models and corrections are adequate.", '9. Radial Range in the Lensing Analysis': "Our fiducial setup ( r > 500 h -1 kpc) is a compromise between avoiding the cluster central regions (where miscentering, cluster member contamination, baryonic feedback, and possible non-linear shear responses are particularly challenging to model), and including scales with high weak-lensing signal-to-noise ratios [47]. We validate this choice by comparing the fiducial results to an analysis with an even more conservative radial cut, r > 800 h -1 kpc. The outer limit remains unchanged: r < 3 . 2 (1 + z cluster ) -1 h -1 Mpc. As expected, the recovered parameter constraints slightly broaden compared to the fiducial setup, but the central values remain almost identical (last row in Fig. 3, 'DES WL r min : 500 → 800 h -1 Mpc). This confirms that our model can adequately describe cluster shear profiles down to r = 500 h -1 kpc.", 'C. Unblinding': 'The validation tests discussed in the previous section were designed to probe what we expected to be the most probable points of failure of our model. We performed these tests while being blinded to the actual values of the mean recovered parameter constraints. Since no alternative analysis led to significant parameter shifts, we proceeded and removed the so-far unknown blinding parameter offsets. Throughout this work, we present the results as obtained after unblinding, with no further modification of the data or the analysis.  \nFIG. 5. Constraints on Ω m and σ 8 from the abundance of SPT-selected clusters (68% and 95% credible regions). In comparison with the precursor analysis of 343 SPT-SZ clusters (gray contours; [7]), the larger cluster sample and the improved weak-lensing mass calibration enable tighter constraints (red contours). We also show the constraints obtained from the SPT cluster abundance with DES Y3 lensing, without using the HST-39 high-redshift cluster lensing dataset (blue contours). Green contours show the constraints obtained when restricting the cluster sample to ξ > 5, as was done in precursor analyses (this result also corresponds to item 6 in Sec. IV).  \n<!-- image -->', 'V. RESULTS': 'We present and discuss constraints on the ΛCDM and w CDM models, and on the scaling relation parameters. Our baseline setup includes the sum of neutrino masses ∑ m ν as a free parameter; however, we can only constrain it in combination with measurements of primary CMB anisotropies from Planck . We also present constraints on an empirical extension of the ΛCDM model, in which we allow the rate of structure growth to deviate from the standard prediction. We implement this test as a direct measurements of σ 8 at different redshifts between 0.25 and 1.8.', 'A. ΛCDM with Massive Neutrinos': "In Fig. 5, we show the constraints on Ω m and σ 8 as derived from the SPT cluster abundance. We note that supplementing the DES lensing data with high-redshift cluster lensing data from HST leads to a slight shift of the constraints along the degeneracy axis. The results obtained from the fiducial, combined analysis setup-  \nFIG. 6. Constraints on Ω m and σ 8 from various analyses of the cluster abundance (68% and 95% credible regions). Along with our work, we also show results obtained from X-ray selected clusters (ROSAT [89], eFEDS [90]) and optically selected clusters (DES Y1 [81, 91], KiDS [92]). Our analysis places the tightest constraints on σ 8 and S opt 8 .  \n<!-- image -->  \nSPT(SZ+pol) cluster abundance with (DES Y3 + HST39) lensing-are summarized in Table II. As demonstrated in Fig. 5, our new analysis significantly improves over the precursor study, which used 343 clusters in the SPT-SZ survey and lensing information from 32 targeted observations [7]. Our new results are shifted toward higher values of σ 8 , which is caused by using DES Y3 lensing data for mass calibration. Remember that our analysis also uses a larger cluster sample than the predecessor analysis, but the SPT cluster samples are consistent with one another, ruling out the possibility that the shift in σ 8 can be explained with changes in the cluster sample [30, 37]. 12 Finally, we note that neglecting the systematic uncertainties in the lensing model only has a small impact on the recovered constraints (first vs. second row in Fig. 3). This implies that our constraints are still limited by the precision of the lensing measurements, rather than the uncertainties in the modeling. The statistical uncertainty will be improved in future analyses by increasing the number of lenses (larger cluster sample), and/or by using lensing datasets with higher source densities and higher lensing efficiencies.  \nTABLE II. Parameter constraints for the ΛCDM and w CDM models (mean and 68% credible intervals, or 95% limit). The SPTclusters+lensing dataset corresponds to SPT(SZ+pol) clusters + (DES Y3 + HST-39) WL. The parameters n s and Ω b h 2 are not meaningfully constrained by the cluster dataset and the results are thus omitted here. Similarly, the cluster dataset cannot constrain h and ∑ m ν on its own, and we only quote the results recovered in combination with Planck 2018 TT,TE,EE+lowE [2].  \nWe compare our results with other analyses of the cluster abundance and lensing mass calibration (see Fig. 6). The results from the DES Y1 cluster abundance and lensing are shifted toward lower values of Ω m than expected from most cosmological probes [81], and are wider than the constraints we report here. The same cluster sample, cross-correlated with clusters, galaxies, and shear on large scales ('DES Y1 clusters + 4 × 2 pt'), yields constraints that largely overlap with ours [91]. The analyses of optically selected clusters in the KiDS survey with KiDS lensing [92], and the analysis of clusters selected in the eROSITA Science Verification data (eFEDS) with HSC lensing [90, 93] both exhibit constraints that are less tightly constrained along the direction of the Ω m -σ 8 degeneracy than our analysis, but they overlap with our results. The result obtained from ROSAT clusters with Subaru/HSC lensing 'Weighing the Giants' (WtG) is almost centered on ours, but is somewhat less constraining along the width of the degeneracy [89]. We note that their constraint also includes a constraint on Ω m through measurements of gas mass fractions f gas , whereas our constraints are solely derived from the cluster abundance.  \nIn Fig. 7, we compare our constraints in the Ω m -σ 8 plane with those obtained from other cosmological probes. Current measurements on CMB lensing from ACTand SPT cannot break the Ω m -σ 8 parameter degeneracy, but their degeneracy goes through our constraints [16, 94, 95]. Our measurement S opt 8 = 0 . 805 ± 0 . 016 and S opt 8 = 0 . 818 ± 0 . 022 from ACT DR6 CMB lensing alone [17] agree to within 0 . 5 σ . As an example of joint analyses of weak lensing and galaxy clustering, we show the DES Y3 3 × 2 pt results [9]. They constrain S opt 8 = 0 . 754 ± 0 . 031, which is 1 . 5 σ lower than our measurement. In the Ω m -σ 8 plane, the difference between their results and ours has a probability to exceed of 0.25 (1 . 1 σ ). 13 Finally, Planck 2018 TT,TE,EE+lowE data (without lensing) constrain S opt 8 = 0 . 820 ± 0 . 016 [2],  \nFIG. 7. Constraints on Ω m and σ 8 (68% and 95% credible regions) in a ΛCDM universe with massive neutrinos from the abundance of SPT clusters, a selection of large-scale structure lensing probes [9, 17], Planck primary CMB measurements (TT,TE,EE+lowE) [2], and BAO [96-99]. Our measurement of the cluster abundance places competitive constraints especially on S opt 8 .  \n<!-- image -->  \nwhich differs from our measurement by 0 . 7 σ , and with an uncertainty that is equal to ours. In the Ω m -σ 8 plane, the difference is essentially identical (probability to exceed of 0 . 47 , 0 . 7 σ ). For the S 8 parameter, the difference is 1 . 1 σ , which does not strengthen the claim of an ' S 8 tension' between large-scale structure measurements at late times and the primary CMB. We also show constraints from baryonic acoustic oscillations (BAO), which constrain a combination of Ω m and H 0 , but are unable to constrain  \nFIG. 8. Constraints (68% and 95% credible regions) on Ω m , σ 8 , and on the sum of neutrino masses ∑ m ν from the abundance of SPT clusters and its combination with Planck primary CMB anisotropies (TT,TE,EE+lowE) [2]. The joint analysis places an upper limit ∑ m ν < 0 . 18 eV (95% limit).  \n<!-- image -->", 'σ 8 . 14': 'We proceed and combine our cluster abundance measurements with primary CMB anisotropies ( Planck 2018 TT,TE,EE+lowE). As shown in Fig. 8, the Planck analysis exhibits a well-known degeneracy between the sum of neutrino masses ∑ m ν , Ω m , and σ 8 . Conversely, the cluster abundance-like other probes of the largescale structure-constrains σ 8 essentially independently of ∑ m ν . The combined analysis therefore leads to improved constraints by breaking that degeneracy. Our measurement of σ 8 fully encompasses the Planck measurement, and the cluster measurements only help to exclude the high-∑ m ν tail. We report an improved constraint on the upper limit,  \n∑ m ν < 0 . 18 eV (95% limit) , (12)  \ncompared to ∑ m ν < 0 . 26 eV from Planck data alone. The constraints on the other cosmological parameters only improve mildly (see Table II).', 'B. Constraints on the Observable-Mass Relations': 'In our analysis, we simultaneously calibrate the observable-mass relations while fitting for cosmology. In Table III, we present the recovered constraints on the parameters of the relations defined in Eqs. (3) and (5). In the Appendix, we show marginalized parameter constraints for all scaling relation and cosmological parameters that are not prior-dominated.', '1. The SZ-Mass Relation': 'The SZ parameter constraints are similar to, but tighter than the ones presented in the precursor SPT cluster analysis based on targeted lensing observations [7], with the exception of the amplitude A SZ . 15 Our new amplitude is somewhat lower than the old result, which also explains the shift toward larger σ 8 seen in Fig. 5 (a lower A SZ implies that the observed clusters are more massive). As discussed in Sec. V A, this shift is caused by the updated weak-lensing mass calibration using DES Y3 data. Finally, our constraint on γ ECS , the relative depth of the SPTpol ECS survey, is centered on the value obtained from an SZ-richness scaling relation analysis (see Sec. 6.1.3 in [30]), but we recover a tighter constraint, most likely due to simultaneously analysing the cluster abundance and multi-observable scaling relations.', '2. The Richness-Mass Relation': 'Throughout this work, we use richness measurements as computed by MCMF (see Sec. II). At low cluster redshift ( z cluster < 1 . 1), where our richness measurements are based on DES data, the recovered richness-mass relation implies that a halo of mass M 200c = 3 × 10 14 h -1 M ⊙ has a typical richness of 40. We find a power-law dependence with mass that is steeper than unity at > 3 σ , and a dependence with redshift that is consistent with no evolution at < 1 σ . The intrinsic scatter in ln ˜ λ is small (0 . 18 ± 0 . 03), but note that the relation between observed richness λ and mass contains an additional lognormal scatter of width ˜ λ -1 / 2 [see Eq. (6)]. 16  \nThe parameters of the mean richness-mass relation at high redshift z > 1 . 1, where we use WISE data, do not quite match those of the low-redshift richness-mass relation. This is not entirely surprising, given how challenging it is to measure richness at these high redshifts and with the modest angular resolution of WISE. Interestingly though, the intrinsic scatter is smaller than for  \nTABLE III. Constraints on the parameters of the observablemass relations (mean and 68% credible interval or 95% limit). We do not show the prior-dominated parameters of the M WL -M halo relation. Note that the constraints on the correlation coefficients are heavily affected by the prior U ( -0 . 5 , 0 . 5).  \nDES-based richness at 1 . 5 σ . The mass slope is consistent with unity within 1 σ , whereas the preferred redshift evolution is negative at the ∼ 3 σ level. The significant difference between the richness measurements from the two observatories motivates the separate empirical calibration of each relation, as performed in this work.  \nWe discuss our low-redshift richness-mass relation in the light of existing DES results. These results are based on richness measurements using redMaPPer [100], whereas our richness measurements are based on MCMF, with median ( λ MCMF /λ redMaPPer ) = 1 . 09 and 24% scatter [101]. The redMaPPer richness-mass relation was determined by cross-matching with SPT-confirmed clusters, and calibrating the SZ-mass relation assuming a fixed fiducial cosmology [30, 102], and through an analysis of the DES Y1 cluster abundance with mass calibration enabled by cross-matching with SPT clusters, and using lensing data for SPT clusters [103]. Qualitatively, our richness-mass relation behaves as in these analyses, with residual differences being explained by the different richness measurements and by the fact that in this work, we calibrate the observable-mass relations using different lensing data (that is, DES Y3). The DES Y1 weaklensing based calibration of the redMaPPer richnessmass relation [104] was shown to yield a significantly  \nshallower mass slope ( B λ ≈ 1 / 1 . 356 = 0 . 74) than the SPT-based analyses (e.g., B λ = 1 . 020 ± 0 . 080 in [30]). Within the limits of comparability, it appears that our constraint on the mass slope is in significant tension with the DES Y1 lensing result [104].', '3. The Correlation Coefficients of Intrinsic Scatter': 'The constraints on the correlation coefficients (bottom rows in Table III) are heavily affected by the uniform priors ρ ∼ U ( -0 . 5 , 0 . 5) that we apply to make the covariance matrix nonsingular. Note that as part of the blinding tests, we investigated the impact of fixing a set of fit parameters. This included setting the correlation coefficients to 0. The recovered scaling relation and cosmology constraints did not exhibit any significant shifts compared to the full analysis (see Fig. 3). We therefore conclude that the exact range of the priors that we apply to the correlation coefficients does not have a significant impact on the recovered constraints on other parameters.  \nThe posterior distributions of the two correlation coefficients that involve the SZ observable ( ρ SZ , ˜ λ and ρ SZ , WL ) both peak at the lower bound ρ = -0 . 5 and we present the 95% upper limits in Table III. A negative correlation coefficient ρ SZ , ˜ λ < 0 is also found in the analysis of the LoCuSS cluster sample [105], whereas simulation works tend to favor positive correlation coefficients between the SZ signal and other observables, e.g., [58, 106, 107]. Negative correlation coefficients would imply that cluster triaxiality-which causes positive correlation-is not driving the joint scatter. Our constraint on the correlation coefficient between lensing and richness ( ρ WL , ˜ λ ) peaks within the prior range, but the quoted constraints are also strongly affected by the prior.  \nIn our hierarchical modeling of the cluster population, the correlation coefficients are degenerate with other parameters of the observable-mass relations. Within our approach of empirical mass calibration, we do not expect strong constraints on the correlated scatter from current datasets (see also Sec. 4.1 in [108]).', '4. Observable-Mass Relations from Joint Analysis with Planck primary CMB': 'In the last column in Table III, we show the constraints obtained from the joint analysis of the cluster dataset and Planck primary CMB anisotropies. Some parameters ( B SZ , C SZ , C λ, DES ) shift within 1 σ while most remain essentially unchanged. This is explained by the fact that our cluster-based cosmology constraints overlap with the Planck results (see Fig. 7). Therefore, their combination does not pull the parameters away from their mean in the clusters-only analysis. The uncertainties on most scaling relation parameters only tighten mildly when adding Planck CMB data. However, there is a noticeable improvement on the constraints on the redshift evolutions  \nC SZ and C λ, DES . These parameters are degenerate with Ω m (see Fig. 11), and the uncertainty on that parameter significantly tightens when adding Planck CMB data, as shown in Table II and Fig. 8.', 'C. w CDM with Massive Neutrinos': 'We consider an extension of the ΛCDM model by opening up the dark energy equation of state parameter w as an additional degree of freedom. As before, we also marginalize over the sum of neutrino masses. The constraints are summarized in Table II and shown in Fig. 9. Considering only the cluster dataset, we note that the constraints on Ω m and σ 8 are very similar to those obtained assuming the ΛCDM model. We obtain w = -1 . 45 ± 0 . 31, a constraint that is more negative than that corresponding to a cosmological constant, with a probability to exceed of 0.10 (1 . 7 σ ) (note that the posterior distribution is affected by the hard cut w > -2). 17 We point out that the uncertainty is comparable to the DES Y3 3 × 2 pt result ( w = -0 . 98 +0 . 32 -0 . 20 ), underscoring the constraining power of cluster abundance measurements. Figure 9 also shows wide degeneracies between Ω m , σ 8 , and w in the analysis of Planck primary CMB data. In combination with the cluster dataset, some of these degeneracies can be broken, and the constraints are tightened. 18 The joint constraints are qualitatively similar to those obtained from the cluster dataset alone, except for w , for which the uncertainty reduces significantly. The mean recovered value for w is less negative for the joint analysis, but the reduced uncertainty leads to a larger statistical difference with a cosmological constant (probability to exceed of 0.03, 2 . 2 σ ) than for the clusters-only analysis.  \nOverall, the recovered constraints on w are still quite weak. As expected, w exhibits a rather narrow degeneracy with the parameter of the SZ redshift evolution C SZ (see Fig. 11), which is only poorly constrained with the current weak-lensing data. In w CDM, we observe a shift to larger values ( C SZ = 1 . 00 ± 0 . 32), and away from its preferred value in ΛCDM ( C SZ = 0 . 50 ± 0 . 27). More powerful lensing data, especially at intermediate and high cluster redshifts, or a tight informative prior on the evolution of the ICM with redshift (as applied in the WtG analysis [89]) would enable tighter constraints on C SZ , and thus on w . Alternatively, a measurement of the gas fraction f gas (as also applied in the WtG analysis) would constrain Ω m independently, and through the Ω m -w parameter degeneracy, reduce the uncertainty on w .  \nw  \nFIG. 9. Constraints on Ω m , σ 8 , and the dark energy equation of state parameter w (68% and 95% credible regions) in a universe with massive neutrinos. We show results obtained from the analysis of the SPT cluster abundance, Planck primary CMB anisotropies, and their combination. The cluster dataset prefers a value that is more negative than a cosmological constant w = -1 at the 1 . 7 σ level. Note that the constraints on Ω m and σ 8 are only mildly affected by opening up w as an additional degree of freedom.  \n<!-- image -->', 'D. Modified Growth of Structure: σ 8 ( z )': 'We consider a second, and final, extension of the ΛCDM model, in which we allow the rate of structure growth to depart from the standard prediction. We do so by directly measuring σ 8 at different redshifts. Because the SPT cluster sample spans a wide range of redshifts, this exercise allows us to map the growth of cosmic structures since the universe was about 3 billion years old.  \nWe define five bins in redshift such that each bin contains an equal number of clusters. The six bin edges are z ∈ { 0 . 25 , 0 . 42 , 0 . 57 , 0 . 71 , 0 . 93 , 1 . 81 } . Then, we modify the computation of the matter power spectrum such that σ 8 at each bin center is normalized to one of five new fit parameters σ 8 ( z i ) (see Table IV). The observablemass relations and the lensing data are treated the same as in the rest of this work, and we also marginalize over the sum of neutrino masses. When analyzing the cluster dataset in this way, the recovered constraint on Ω m degrades significantly compared to the ΛCDM or w CDM analyses. 19 Therefore, we consider an additional  \nTABLE IV. Constraints on Ω m and σ 8 ( z ) for a ΛCDM cosmology with massive neutrinos (mean and 68% credible interval). Using the abundance of SPT clusters with DES Y3 + HST weak-lensing mass calibration, and a prior on the sound horizon at recombination θ ∗ from Planck , we constrain σ 8 at discrete redshifts. We also quote σ 8 ( z ) for ΛCDM and parameters from Planck 2018 TT,TE,EE+lowE.  \nconstraint on the matter density through the measurement of the sound horizon at recombination 100 θ ∗ = 1 . 04109 ± 0 . 00030 by Planck [2]. 20  \nWe summarize our results in Fig. 10. The gray band shows the 68% credible interval of σ 8 as a function of redshift, assuming the ΛCDM model with massive neutrinos, and parameters as measured by Planck primary CMB anisotropies. Red data points show direct measurements of σ 8 at five redshifts, as probed by the cluster abundance in combination with the θ ∗ prior. Note that the θ ∗ prior represents a subset of the information contained in the Planck CMB data, and we are thus comparing growth measurements for similar background cosmologies (thereby avoiding regimes where differences in σ 8 ( z ) could be caused by differences in Ω m instead of differences in structure growth). Blue data points are obtained from a cross-correlation analysis of DESI LRGs and Planck CMB lensing, with θ ∗ fixed [109]. 21 Brown data points show the constraint from DES Y3 galaxy clustering and lensing (3 × 2 pt) and shear ratios [110].  \nOur measurements of σ 8 ( z ) presented in Table IV result in a difference with the ΛCDM prediction with Planck parameters of χ 2 = 2 . 8 for five parameters, which indicates good agreement. However, the recovered parameter constraints are degenerate [the correlation coefficients between the individual σ 8 ( z i ) are all > 0 . 76], and taking these parameter correlations into account increases the χ 2 value. We anticipate that a more precise determination of the observable-mass relations, enabled by future lensing dataset, will reduce the level of correlation among the measured σ 8 ( z i ) and will therefore enable  \nFIG. 10. Evolution of σ 8 with time. Red diamonds show constraints from the SPT cluster abundance (mean and 68% credible interval). Blue dots are obtained from the crosscorrelation of DESI LRGs and Planck CMB lensing [109]. Both analyses include a Planck -based prior on the sound horizon at recombination θ ∗ , with the value of θ ∗ fixed to the mean Planck value for DESI × Planck , and marginalized over the Planck posterior on θ ∗ for SPT clusters. Brown squares show the constraints from DES Y3 3 × 2 pt + shear ratios [110]. The gray band shows the 68% credible interval of the ΛCDM prediction, assuming parameters determined by Planck . Our measurement of the growth of structure is consistent with that prediction.  \n<!-- image -->  \nmore robust tests of the evolution of growth of structure.', 'VI. SUMMARY': "In this work, we present cosmological constraints derived from the abundance of galaxy clusters selected in the SPT-SZ and SPTpol surveys [30, 32, 36, 37], and a simultaneous weak-lensing based mass calibration using DES Y3 and HST data. The cluster sample (1,005 confirmed clusters above z > 0 . 25), the analysis of DES weak-lensing data for 688 sample clusters, the Bayesian analysis framework, and the validation of the analysis pipeline are presented in companion Paper I. The HST lensing measurements for 39 SPT clusters and their analysis were presented in earlier works [34, 50-53].  \nIn a first phase of this work, we performed a blinded analysis in which we artificially offset the parameter constraints from their true output values by unknown blinding offset parameters (see Sec. IV). This allowed us to perform a series of robustness tests with which we investigated the relative impact of alternative analysis setups on the results without knowing the actual parameter values. These tests included various cuts in the cluster sample (excluding high-redshift and/or low-mass objects), details of the modeling of the SPT observable, and three alternative analyses of the DES Y3 weak-lensing data.  \nSince none of these alternative analyses resulted in constraints that were significantly offset from the fiducial ones, we unblinded the results. We did not further modify the data or the analysis, and we present the results as obtained after unblinding.  \nAssuming a flat ΛCDM cosmology with massive neutrinos, we find that our constraints are consistent with, but tighter than, results obtained from previous weaklensing calibrated cluster abundance measurements (see Fig. 6 and Table II). Our constraints are competitive with other leading probes of the large-scale structure such as cosmic shear, joint 3 × 2 pt analyses, and the CMB lensing power spectrum. Our dataset best constrains the parameter combination σ 8 (Ω m / 0 . 3) 0 . 25 = 0 . 805 ± 0 . 016, with an uncertainty that is equal to the one recovered from Planck 2018 TT,TE,EE+lowE data. The difference between our constraint on S 8 ≡ σ 8 √ Ω m / 0 . 3 and the one from Planck is 1 . 1 σ , which does not strengthen the claim of an ' S 8 tension'. The difference in the Ω m -σ 8 plane has a probability to exceed of 0.47 (0 . 7 σ ). The combination with the Planck dataset further tightens parameter constraints, and we place a 95% upper limit on the sum of neutrino masses ∑ m ν < 0 . 18 eV.  \nAdditionally allowing the dark energy equation of state parameter w to vary does not change the recovered constraints on Ω m and σ 8 in any significant way. We recover w = -1 . 45 ± 0 . 31. The analysis of Planck primary CMB data exhibits a degeneracy between w and Ω m , σ 8 , and ∑ m ν , which is broken in a joint analysis with our cluster dataset. The resulting constraint on w is still low, and we report w = -1 . 34 +0 . 22 -0 . 15 , which differs from a cosmological constant by 2 . 2 σ .  \nFinally, we consider a phenomenological extension of the ΛCDM model, in which we directly constrain the growth of structure. Using the cluster dataset and a prior on the sound horizon at recombination θ ∗ , we constrain σ 8 at five discrete redshifts between 0.34 and 1.35. Our measurements of σ 8 ( z ) are consistent with the ΛCDM prediction with Planck 2018 TT,TE,EE+lowE parameters (see Fig. 10 and Table IV), but we note that our individual measurements are highly correlated.  \nAdistinctive feature of this analysis is that, for the first time, we leverage the exquisite quality of weak-lensing data as provided by wide-field galaxy lensing surveys to perform a robust and precise weak-lensing informed measurement of the abundance of SZ-selected clusters. In future work, we will use our cluster and lensing dataset to calibrate other cluster observables (X-rays, dynamics of cluster member galaxies) and to investigate further extensions of the ΛCDM model, such as dark matterdark radiation interactions and modified gravity models. We also note that our analysis is the first cluster cosmology work that uses DES Y3 lensing data. This lensing dataset is also being used for mass calibration of redMaPPer-selected clusters in the DES Y3 data and for mass calibration of further ICM-selected cluster samples. Therefore, the DES lensing data play an important role in enabling cosmological constraints from the abundance  \nof galaxy clusters.  \nOur measurements are still limited by size of the cluster sample, and by shape noise in the lensing measurements. Therefore, increasing the cluster sample (thereby also increasing the number of lenses), increasing the size (and depth) of the lensing source sample, and the additional analysis of cluster lensing of the CMB e.g., [111113], will lead to improved cosmological constraints [114]. However, our analysis also highlights areas that deserve further attention (e.g., low-significance clusters, cluster member contamination, miscentering). While the relative importance of these effects is currently at the ≲ 1 σ level (see Fig. 3), more work is needed to ensure that these do not become limiting factors in future analyses. Opportunities for near-term improvements are numerous, with large ICM-selected cluster samples from the Atacama Cosmology Telescope [115], eROSITA [116], the SPT-3G survey [117], and the upcoming Simons Observatory [118]. These cluster samples will be analyzed jointly with a weak-lensing mass calibration using DES Year 6 data, soon to be followed by galaxy and lensing data from Euclid [119] and the Legacy Survey of Space and Time (LSST) [120] at the Vera C. Rubin Observatory, whose combination is expected to be particularly powerful for extending sensitive cluster mass calibrations to higher redshifts [121]. Another promising avenue is to additionally consider the cluster 2-pt correlation function (and cluster-galaxy 2 pt), which further boosts the constraining power [91, 122]. Finally, the CMB-S4 survey will improve upon the aforementioned mm-wave surveys, and will extend analyses of SZ-selected cluster samples to well into the 2030s [123].", 'ACKNOWLEDGMENTS': "This research was supported by the Excellence Cluster ORIGINS, which is funded by the Deutsche Forschungsgemeinschaft (DFG, German Research Foundation) under Germany's Excellence Strategy -EXC-2094390783311, the MPG Faculty Fellowship program and the Ludwig-Maximilians-Universitat Munchen. Parts of the MCMC computations were carried out on the computing facilities of the Computational Center for Particle and Astrophysics (C2PAP). The Bonn and Innsbruck authors acknowledge support from the German Federal Ministry for Economic Affairs and Energy (BMWi) provided through DLR under projects 50OR2002 and 50OR2302, from the German Research Foundation (DFG) under grant 415537506, and the Austrian Research Promotion Agency (FFG) and the Federal Ministry of the Republic of Austria for Climate Action, Environment, Mobility, Innovation and Technology (BMK) via grants 899537 and 900565. The Melbourne authors acknowledge support from the Australian Research Council's Discovery Projects scheme (No. DP200101068).  \nThe South Pole Telescope program is supported by the National Science Foundation (NSF) through the Grant  \nNo. OPP-1852617. Partial support is also provided by the Kavli Institute of Cosmological Physics at the University of Chicago. PISCO observations were supported by US NSF grant AST-0126090. Work at Argonne National Laboratory was supported by the U.S. Department of Energy, Office of High Energy Physics, under Contract No. DE-AC02-06CH11357.  \nFunding 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¸c˜ao Carlos Chagas Filho de Amparo 'a Pesquisa do Estado do Rio de Janeiro, Conselho Nacional de Desenvolvimento Cient'ıfico e Tecnol'ogico and the Minist'erio da Ciˆencia, Tecnologia e Inova¸c˜ao, the Deutsche Forschungsgemeinschaft and the Collaborating Institutions in the Dark Energy Survey.  \nThe Collaborating Institutions are Argonne National Laboratory, the University of California at Santa Cruz, the University of Cambridge, Centro de Investigaciones Energ'eticas, Medioambientales y Tecnol'ogicas-Madrid, the University of Chicago, University College London, the DES-Brazil Consortium, the University of Edinburgh, the Eidgenossische Technische Hochschule (ETH) Zurich, Fermi National Accelerator Laboratory, the University of Illinois at Urbana-Champaign, the Institut de Ci'encies de l'Espai (IEEC/CSIC), the Institut de F'ısica d'Altes Energies, Lawrence Berkeley National Laboratory, the Ludwig-Maximilians-Universitat Munchen and the associated Excellence Cluster Origins, the University of Michigan, NSF'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.  \nBased in part on observations at Cerro Tololo InterAmerican 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.  \nThe 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, PGC2018094773, PGC2018-102021, SEV-2016-0588, SEV-20160597, and MDM-2015-0509, some of which include ERDF funds from the European Union. IFAE is partially funded  \nby 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/20072013) including ERC grant agreements 240672, 291329, and 306478. We acknowledge support from the Brazilian Instituto Nacional de Ciˆencia e Tecnologia (INCT) do e-Universo (CNPq grant 465376/2014-2).  \nThis work is based on observations made with the NASA/ESA Hubble Space Telescope , using imaging data from the SPT follow-up GO programs 12246 (PI: C. Stubbs), 12477 (PI: F. W. High), 13412 (PI: T. Schrabback), 14252 (PI: V. Strazzullo), 14352 (PI: J. Hlavacek-Larrondo), and 14677 (PI: T. Schrabback). STScI is operated by the Association of Universities for Research in Astronomy, Inc. under NASA contract NAS 5-26555. It is also based on observations made with ESO Telescopes at the La Silla Paranal Observatory under programs 086.A-0741 (PI: Bazin), 088.A0796 (PI: Bazin), 088.A-0889 (PI: Mohr), 089.A-0824 (PI: Mohr), 0100.A-0204 (PI: Schrabback), 0100.A-0217 (PI: Hern'andez-Mart'ın), 0101.A-0694 (PI: Zohren), and 0102.A-0189 (PI: Zohren). It is also based on observations obtained at the Gemini Observatory, which is operated by the Association of Universities for Research in Astronomy, Inc., under a cooperative agreement with the NSF on behalf of the Gemini partnership: the National Science Foundation (United States), National Research Council (Canada), CONICYT (Chile), Ministerio de Ciencia, Tecnolog'ıa e Innovaci'on Productiva (Argentina), Minist'erio da Ciˆencia, Tecnologia e Inova¸c˜ao (Brazil), and Korea Astronomy and Space Science Institute (Republic of Korea), under programs 2014B-0338 and 2016B-0176 (PI: B. Benson).  \nThis 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. This research has made use of the SAO/NASA Astrophysics Data System and of adstex. 22 We thank Ludwig, the SPT support cat for this analysis. Inquiries about SPT support cats shall be directed to T.C.", 'Appendix A: Full Parameter Posterior Distributions': "In Fig. 11, we show the marginalized one- and two-dimensional parameter constraints obtained in the ΛCDM and w CDM analyses. We show all parameters that are not strongly prior-dominated. We highlight the joint constraint on Ω m and σ 8 , which is quite robust to the extension with w as a free parameter. We also highlight the degeneracy between C SZ and w , which implies that a tighter constraint on the redshift evolution  \nstraints on the dark energy equation of state.  \n- [1] A. G. Riess, S. Casertano, W. Yuan, L. Macri, J. Anderson, J. W. MacKenty, J. B. Bowers, K. I. Clubb, A. V. Filippenko, D. O. Jones, et al., ApJ 855 , 136 (2018), arXiv:1801.01120.\n- [2] Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, et al., A&A 641 , A6 (2020), arXiv:1807.06209.\n- [3] A. G. Riess, Nature Reviews Physics 2 , 10 (2020), arXiv:2001.03624.\n- [4] A. G. Riess, W. Yuan, L. M. Macri, D. Scolnic, D. Brout, S. Casertano, D. O. Jones, Y. Murakami, G. S. Anand, L. Breuval, et al., ApJ 934 , L7 (2022), arXiv:2112.04510.\n- [5] DES Collaboration, T. M. C. Abbott, F. B. Abdalla, A. Alarcon, J. Aleksi'c, S. Allam, S. Allen, A. Amara, J. Annis, J. Asorey, et al., Phys. Rev. D 98 , 043526 (2018), arXiv:1708.01530.\n- [6] E. van Uitert, B. Joachimi, S. Joudaki, A. Amon, C. Heymans, F. Kohlinger, M. Asgari, C. Blake, A. Choi, T. Erben, et al., MNRAS 476 , 4662 (2018), arXiv:1706.05004.\n- [7] S. Bocquet, J. P. Dietrich, T. Schrabback, L. E. Bleem, M. Klein, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, et al., ApJ 878 , 55 (2019), arXiv:1812.01679.\n- [8] C. Heymans, T. Troster, M. Asgari, C. Blake, H. Hildebrandt, B. Joachimi, K. Kuijken, C.-A. Lin, A. G. S'anchez, J. L. van den Busch, et al., A&A 646 , A140 (2021), arXiv:2007.15632.\n- [9] DES Collaboration, T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, F. AndradeOliveira, J. Annis, S. Avila, et al., Phys. Rev. D 105 , 023520 (2022), arXiv:2105.13549.\n- [10] Y. Kobayashi, T. Nishimichi, M. Takada, and H. Miyatake, Phys. Rev. D 105 , 083517 (2022), arXiv:2110.06969.\n- [11] R. Dalal, X. Li, A. Nicola, J. Zuntz, M. A. Strauss, S. Sugiyama, T. Zhang, M. M. Rau, R. Mandelbaum, M. Takada, et al., Phys. Rev. D 108 , 123519 (2023), arXiv:2304.00701.\n- [12] X. Li, T. Zhang, S. Sugiyama, R. Dalal, R. Terasawa, M. M. Rau, R. Mandelbaum, M. Takada, S. More, M. A. Strauss, et al., Phys. Rev. D 108 , 123518 (2023), arXiv:2304.00702.\n- [13] Dark Energy Survey and Kilo-Degree Survey Collaboration, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, M. Asgari, S. Avila, D. Bacon, et al., The Open Journal of Astrophysics 6 , 36 (2023), arXiv:2305.17173.\n- [14] S. More, S. Sugiyama, H. Miyatake, M. M. Rau, M. Shirasaki, X. Li, A. J. Nishizawa, K. Osato, T. Zhang, M. Takada, et al., Phys. Rev. D 108 , 123520 (2023), arXiv:2304.00703.\n- [15] Planck Collaboration, N. Aghanim, Y. Akrami, M. Ashdown, J. Aumont, C. Baccigalupi, M. Ballardini, A. J. Banday, R. B. Barreiro, N. Bartolo, et al., A&A 641 , A8 (2020), arXiv:1807.06210.\n- [16] F. Bianchini, W. L. K. Wu, P. A. R. Ade, A. J. Anderson, J. E. Austermann, J. S. Avva, J. A. Beall, A. N. Bender, B. A. Benson, L. E. Bleem, et al., ApJ 888 , 119 (2020), arXiv:1910.07157.\n- [17] F. J. Qu, B. D. Sherwin, M. S. Madhavacheril, D. Han, K. T. Crowley, I. Abril-Cabezas, P. A. R. Ade, S. Aiola, T. Alford, M. Amiri, et al., ApJ 962 , 112 (2024), arXiv:2304.05202.\n- [18] C. Chang, Y. Omori, E. J. Baxter, C. Doux, A. Choi, S. Pandey, A. Alarcon, O. Alves, A. Amon, F. AndradeOliveira, et al., Phys. Rev. D 107 , 023530 (2023), arXiv:2203.12440.\n- [19] Z. Haiman, J. J. Mohr, and G. P. Holder, ApJ 553 , 545 (2001), arXiv:astro-ph/0002336.\n- [20] S. Dodelson, K. Heitmann, C. Hirata, K. Honscheid, A. Roodman, U. Seljak, A. Slosar, and M. Trodden (2016), arXiv:1604.07626.\n- [21] S. W. Allen, A. E. Evrard, and A. B. Mantz, ARA&A 49 , 409 (2011), arXiv:1103.4829.\n- [22] G. W. Pratt, M. Arnaud, A. Biviano, D. Eckert, S. Ettori, D. Nagai, N. Okabe, and T. H. Reiprich, Space Sci. Rev. 215 , 25 (2019), arXiv:1902.10837.\n- [23] J. E. Carlstrom, P. A. R. Ade, K. A. Aird, B. A. Benson, L. E. Bleem, S. Busetti, C. L. Chang, E. Chauvin, H. M. Cho, T. M. Crawford, et al., PASP 123 , 568 (2011), arXiv:0907.4445.\n- [24] B. Flaugher, H. T. Diehl, K. Honscheid, T. M. C. Abbott, O. Alvarez, R. Angstadt, J. T. Annis, M. Antonik, O. Ballester, L. Beaufore, et al., AJ 150 , 150 (2015), arXiv:1504.02900.\n- [25] Dark Energy Survey Collaboration, T. Abbott, F. B. Abdalla, J. Aleksi'c, S. Allam, A. Amara, D. Bacon, E. Balbinot, M. Banerji, K. Bechtol, et al., MNRAS 460 , 1270 (2016), arXiv:1601.00329.\n- [26] T. M. C. Abbott and et al., ApJS 239 , 18 (2018), arXiv:1801.03181.\n- [27] E. L. Wright, P. R. M. Eisenhardt, A. K. Mainzer, M. E. Ressler, R. M. Cutri, T. Jarrett, J. D. Kirkpatrick, D. Padgett, R. S. McMillan, M. Skrutskie, et al., AJ 140 , 1868 (2010), arXiv:1008.0031.\n- [28] R. A. Sunyaev and Y. B. Zeldovich, Comments on Astrophysics and Space Physics 4 , 173 (1972).\n- [29] M. Klein, J. J. Mohr, S. Desai, H. Israel, S. Allam, A. Benoit-L'evy, D. Brooks, E. Buckley-Geer, A. Carnero Rosell, M. Carrasco Kind, et al., MNRAS 474 , 3324 (2018), arXiv:1706.06577.\n- [30] L. E. Bleem, S. Bocquet, B. Stalder, M. D. Gladders, P. A. R. Ade, S. W. Allen, A. J. Anderson, J. Annis, M. L. N. Ashby, J. E. Austermann, et al., ApJS 247 , 25 (2020), arXiv:1910.04121.\n- [31] K. Umetsu, A&A Rev. 28 , 7 (2020), arXiv:2007.00506.\n- [32] L. E. Bleem, B. Stalder, T. de Haan, K. A. Aird, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, B. A. Benson, et al., ApJS 216 , 27 (2015), arXiv:1409.0850.\n- [33] J. P. Dietrich, S. Bocquet, T. Schrabback, D. Applegate, H. Hoekstra, S. Grandis, J. J. Mohr, S. W. Allen, M. B. Bayliss, B. A. Benson, et al., MNRAS 483 , 2871 (2019),  \n2  \nFIG. 11. Parameter constraints on all parameters that are not strongly prior-dominated in the ΛCDM and w CDM analyses (68% and 95% credible regions). Between the two analyses, only the constraints on the redshift evolution parameters of the observable-mass relations C SZ and C λ shift by an appreciable amount. Conversely, the joint constraints on Ω m and σ 8 are very similar (highlighted toward the bottom-right of the plot). The constraint on w shows some significant degeneracies with other parameters, and in particular with the redshift evolution parameters C SZ (highlighted at the bottom-left of the plot) of the SZ-mass relation, and C λ of the richness-mass relation.  \n<!-- image -->", 'arXiv:1711.05344.': "[34] T. Schrabback, D. Applegate, J. P. Dietrich, H. Hoekstra, S. Bocquet, A. H. Gonzalez, A. von der Linden, M. McDonald, C. B. Morrison, S. F. Raihan, et al., MNRAS 474 , 2635 (2018), arXiv:1611.03866.  \n[35] S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, et al., Phys. Rev. D (2024), arXiv:2310.12213. [36] M. Klein, J. J. Mohr, S. Bocquet, M. Aguena, S. W. Allen, O. Alves, B. Ansarinejad, M. L. N.  \n- Ashby, D. Bacon, M. Bayliss, et al., MNRAS (2024), arXiv:2309.09908.\n- [37] L. E. Bleem, M. Klein, T. M. C. Abbot, P. A. R. Ade, M. Aguena, O. Alves, A. J. Anderson, F. AndradeOliveira, B. Ansarinejad, M. Archipley, et al., The Open Journal of Astrophysics 7 , 13 (2024), arXiv:2311.07512.\n- [38] B. Stalder, A. A. Stark, S. M. Amato, J. Geary, S. A. Shectman, C. W. Stubbs, and A. Szentgyorgyi, in Ground-based and Airborne Instrumentation for Astronomy V , edited by S. K. Ramsay, I. S. McLean, and H. Takami (2014), vol. 9147 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series , p. 91473Y.\n- [39] M. Gatti, E. Sheldon, A. Amon, M. Becker, M. Troxel, A. Choi, C. Doux, N. MacCrann, A. NavarroAlsina, I. Harrison, et al., MNRAS 504 , 4312 (2021), arXiv:2011.03408.\n- [40] E. Huff and R. Mandelbaum (2017), arXiv:1702.02600.\n- [41] E. S. Sheldon and E. M. Huff, ApJ 841 , 24 (2017), arXiv:1702.02601.\n- [42] I. Sevilla-Noarbe, K. Bechtol, M. Carrasco Kind, A. Carnero Rosell, M. R. Becker, A. Drlica-Wagner, R. A. Gruendl, E. S. Rykoff, E. Sheldon, B. Yanny, et al., ApJS 254 , 24 (2021), arXiv:2011.03407.\n- [43] M. Jarvis, G. M. Bernstein, A. Amon, C. Davis, P. F. L'eget, K. Bechtol, I. Harrison, M. Gatti, A. Roodman, C. Chang, et al., MNRAS 501 , 1282 (2021), arXiv:2011.03409.\n- [44] N. MacCrann, M. R. Becker, J. McCullough, A. Amon, D. Gruen, M. Jarvis, A. Choi, M. A. Troxel, E. Sheldon, B. Yanny, et al., MNRAS 509 , 3371 (2022), arXiv:2012.08567.\n- [45] S. Everett, B. Yanny, N. Kuropatkin, E. M. Huff, Y. Zhang, J. Myles, A. Masegian, J. Elvin-Poole, S. Allam, G. M. Bernstein, et al., ApJS 258 , 15 (2022), arXiv:2012.12825.\n- [46] J. Myles, A. Alarcon, A. Amon, C. S'anchez, S. Everett, J. DeRose, J. McCullough, D. Gruen, G. M. Bernstein, M. A. Troxel, et al., MNRAS 505 , 4249 (2021), arXiv:2012.08566.\n- [47] S. Grandis, S. Bocquet, J. J. Mohr, M. Klein, and K. Dolag, MNRAS 507 , 5671 (2021), arXiv:2103.16212.\n- [48] J. De Vicente, E. S'anchez, and I. Sevilla-Noarbe, MNRAS 459 , 3078 (2016), arXiv:1511.07623.\n- [49] N. Ben'ıtez, ApJ 536 , 571 (2000), arXiv:astroph/9811189.\n- [50] S. F. Raihan, T. Schrabback, H. Hildebrandt, D. Applegate, and G. Mahler, MNRAS 497 , 1404 (2020), arXiv:2007.01211.\n- [51] B. Hern'andez-Mart'ın, T. Schrabback, H. Hoekstra, N. Martinet, J. Hlavacek-Larrondo, L. E. Bleem, M. D. Gladders, B. Stalder, A. A. Stark, and M. Bayliss, A&A 640 , A117 (2020), arXiv:2007.00386.\n- [52] T. Schrabback, S. Bocquet, M. Sommer, H. Zohren, J. L. van den Busch, B. Hern'andez-Mart'ın, H. Hoekstra, S. F. Raihan, M. Schirmer, D. Applegate, et al., MNRAS 505 , 3923 (2021), arXiv:2009.07591.\n- [53] H. Zohren, T. Schrabback, S. Bocquet, M. Sommer, F. Raihan, B. Hern'andez-Mart'ın, O. Marggraf, B. Ansarinejad, M. B. Bayliss, L. E. Bleem, et al., A&A 668 , A18 (2022), arXiv:2208.10232.\n- [54] A. Mantz, S. W. Allen, D. Rapetti, and H. Ebeling, MNRAS 406 , 1759 (2010), arXiv:0909.3098.  \n- [55] B. A. Benson, T. de Haan, J. P. Dudley, C. L. Reichardt, K. A. Aird, K. Andersson, R. Armstrong, M. L. N. Ashby, M. Bautz, M. Bayliss, et al., ApJ 763 , 147 (2013), arXiv:1112.5435.\n- [56] S. Bocquet, A. Saro, J. J. Mohr, K. A. Aird, M. L. N. Ashby, M. Bautz, M. Bayliss, G. Bazin, B. A. Benson, L. E. Bleem, et al., ApJ 799 , 214 (2015), arXiv:1407.2942.\n- [57] N. Kaiser, MNRAS 222 , 323 (1986).\n- [58] R. E. Angulo, V. Springel, S. D. M. White, A. Jenkins, C. M. Baugh, and C. S. Frenk, MNRAS 426 , 2046 (2012), arXiv:1203.3216.\n- [59] K. Vanderlinde, T. M. Crawford, T. de Haan, J. P. Dudley, L. Shaw, P. A. R. Ade, K. A. Aird, B. A. Benson, L. E. Bleem, M. Brodwin, et al., ApJ 722 , 1180 (2010), arXiv:1003.0003.\n- [60] J. F. Navarro, C. S. Frenk, and S. D. M. White, ApJ 490 , 493 (1997), arXiv:astro-ph/9611107.\n- [61] M. R. Becker and A. V. Kravtsov, ApJ 740 , 25 (2011), arXiv:1011.1681.\n- [62] M. Oguri and T. Hamana, MNRAS 414 , 1851 (2011), arXiv:1101.0650.\n- [63] J. Tinker, A. V. Kravtsov, A. Klypin, K. Abazajian, M. Warren, G. Yepes, S. Gottlober, and D. E. Holz, ApJ 688 , 709 (2008), arXiv:0803.2706.\n- [64] T. McClintock, E. Rozo, M. R. Becker, J. DeRose, Y.Y. Mao, S. McLaughlin, J. L. Tinker, R. H. Wechsler, and Z. Zhai, ApJ 872 , 53 (2019), arXiv:1804.05866.\n- [65] T. Nishimichi, M. Takada, R. Takahashi, K. Osato, M. Shirasaki, T. Oogi, H. Miyatake, M. Oguri, R. Murata, Y. Kobayashi, et al., ApJ 884 , 29 (2019), arXiv:1811.09504.\n- [66] S. Bocquet, K. Heitmann, S. Habib, E. Lawrence, T. Uram, N. Frontiere, A. Pope, and H. Finkel, ApJ 901 , 5 (2020), arXiv:2003.12116.\n- [67] M. Hirschmann, K. Dolag, A. Saro, L. Bachmann, S. Borgani, and A. Burkert, MNRAS 442 , 2304 (2014), arXiv:1308.0333.\n- [68] A. F. Teklu, R.-S. Remus, K. Dolag, A. M. Beck, A. Burkert, A. S. Schmidt, F. Schulze, and L. K. Steinborn, ApJ 812 , 29 (2015), arXiv:1503.03501.\n- [69] A. M. Beck, G. Murante, A. Arth, R. S. Remus, A. F. Teklu, J. M. F. Donnert, S. Planelles, M. C. Beck, P. Forster, M. Imgrund, et al., MNRAS 455 , 2110 (2016), arXiv:1502.07358.\n- [70] K. Dolag, E. Mevius, and R.-S. Remus, Galaxies 5 , 35 (2017), arXiv:1708.00027.\n- [71] A. Pillepich, D. Nelson, L. Hernquist, V. Springel, R. Pakmor, P. Torrey, R. Weinberger, S. Genel, J. P. Naiman, F. Marinacci, et al., MNRAS 475 , 648 (2018), arXiv:1707.03406.\n- [72] F. Marinacci, M. Vogelsberger, R. Pakmor, P. Torrey, V. Springel, L. Hernquist, D. Nelson, R. Weinberger, A. Pillepich, J. Naiman, et al., MNRAS 480 , 5113 (2018), arXiv:1707.03396.\n- [73] V. Springel, R. Pakmor, A. Pillepich, R. Weinberger, D. Nelson, L. Hernquist, M. Vogelsberger, S. Genel, P. Torrey, F. Marinacci, et al., MNRAS 475 , 676 (2018), arXiv:1707.03397.\n- [74] D. Nelson, A. Pillepich, V. Springel, R. Weinberger, L. Hernquist, R. Pakmor, S. Genel, P. Torrey, M. Vogelsberger, G. Kauffmann, et al., MNRAS 475 , 624 (2018), arXiv:1707.03395.\n- [75] J. P. Naiman, A. Pillepich, V. Springel, E. RamirezRuiz, P. Torrey, M. Vogelsberger, R. Pakmor, D. Nelson, F. Marinacci, L. Hernquist, et al., MNRAS 477 , 1206 (2018), arXiv:1707.03401.\n- [76] D. Nelson, V. Springel, A. Pillepich, V. RodriguezGomez, P. Torrey, S. Genel, M. Vogelsberger, R. Pakmor, F. Marinacci, R. Weinberger, et al., Computational Astrophysics and Cosmology 6 , 2 (2019), arXiv:1812.05609.\n- [77] J. Zuntz, M. Paterno, E. Jennings, D. Rudd, A. Manzotti, S. Dodelson, S. Bridle, S. Sehrish, and J. Kowalkowski, Astronomy and Computing 12 , 45 (2015), arXiv:1409.3409.\n- [78] F. Feroz, M. P. Hobson, and M. Bridges, MNRAS 398 , 1601 (2009), arXiv:0809.3437.\n- [79] J. U. Lange, MNRAS 525 , 3181 (2023), arXiv:2306.16923.\n- [80] P. Lemos, N. Weaverdyck, R. P. Rollins, J. Muir, A. Fert'e, A. R. Liddle, A. Campos, D. Huterer, M. Raveri, J. Zuntz, et al., MNRAS 521 , 1184 (2023), arXiv:2202.08233.\n- [81] DES Collaboration, T. M. C. Abbott, M. Aguena, A. Alarcon, S. Allam, S. Allen, J. Annis, S. Avila, D. Bacon, K. Bechtol, et al., Phys. Rev. D 102 , 023509 (2020), arXiv:2002.11124.\n- [82] N. Gupta, A. Saro, J. J. Mohr, B. A. Benson, S. Bocquet, R. Capasso, J. E. Carlstrom, I. Chiu, T. M. Crawford, T. de Haan, et al., MNRAS 467 , 3737 (2017), arXiv:1605.05329.\n- [83] J. B. Melin, J. G. Bartlett, Z. Y. Cai, G. De Zotti, J. Delabrouille, M. Roman, and A. Bonaldi, A&A 617 , A75 (2018), arXiv:1808.06807.\n- [84] ' I. Zubeldia, J. Chluba, and R. Battye, MNRAS 522 , 5123 (2023), arXiv:2212.07410.\n- [85] A. Hashiguchi, Y. Toba, N. Ota, M. Oguri, N. Okabe, Y. Ueda, M. Imanishi, S. Yamada, T. Goto, S. Koyama, et al., PASJ 75 , 1246 (2023), arXiv:2309.01926.\n- [86] T. de Haan, B. A. Benson, L. E. Bleem, S. W. Allen, D. E. Applegate, M. L. N. Ashby, M. Bautz, M. Bayliss, S. Bocquet, M. Brodwin, et al., ApJ 832 , 95 (2016), arXiv:1603.06522.\n- [87] S. Gratton and A. Challinor, MNRAS 499 , 3410 (2020), arXiv:1911.07754.\n- [88] Euclid Collaboration, G. F. Lesci, M. Sereno, M. Radovich, G. Castignani, L. Bisigello, F. Marulli, L. Moscardini, L. Baumont, G. Covone, et al., A&A 684 , A139 (2024), arXiv:2311.16239.\n- [89] A. B. Mantz, A. von der Linden, S. W. Allen, D. E. Applegate, P. L. Kelly, R. G. Morris, D. A. Rapetti, R. W. Schmidt, S. Adhikari, M. T. Allen, et al., MNRAS 446 , 2205 (2015), arXiv:1407.4516.\n- [90] I. N. Chiu, M. Klein, J. Mohr, and S. Bocquet, MNRAS 522 , 1601 (2023), arXiv:2207.12429.\n- [91] C. To, E. Krause, E. Rozo, H. Wu, D. Gruen, R. H. Wechsler, T. F. Eifler, E. S. Rykoff, M. Costanzi, M. R. Becker, et al., Phys. Rev. Lett. 126 , 141301 (2021), arXiv:2010.01138.\n- [92] G. F. Lesci, F. Marulli, L. Moscardini, M. Sereno, A. Veropalumbo, M. Maturi, C. Giocoli, M. Radovich, F. Bellagamba, M. Roncarelli, et al., A&A 659 , A88 (2022), arXiv:2012.12273.\n- [93] A. Liu, E. Bulbul, V. Ghirardini, T. Liu, M. Klein, N. Clerc, Y. Ozsoy, M. E. Ramos-Ceja, F. Pacaud,  \n- J. Comparat, et al., A&A 661 , A2 (2022), arXiv:2106.14518.\n- [94] M. S. Madhavacheril, F. J. Qu, B. D. Sherwin, N. MacCrann, Y. Li, I. Abril-Cabezas, P. A. R. Ade, S. Aiola, T. Alford, M. Amiri, et al., ApJ 962 , 113 (2024), arXiv:2304.05203.\n- [95] Z. Pan, F. Bianchini, W. L. K. Wu, P. A. R. Ade, Z. Ahmed, E. Anderes, A. J. Anderson, B. Ansarinejad, M. Archipley, K. Aylor, et al., Phys. Rev. D 108 , 122005 (2023), arXiv:2308.11608.\n- [96] F. Beutler, C. Blake, M. Colless, D. H. Jones, L. Staveley-Smith, L. Campbell, Q. Parker, W. Saunders, and F. Watson, MNRAS 416 , 3017 (2011), arXiv:1106.3366.\n- [97] A. J. Ross, L. Samushia, C. Howlett, W. J. Percival, A. Burden, and M. Manera, MNRAS 449 , 835 (2015), arXiv:1409.3242.\n- [98] S. Alam, M. Ata, S. Bailey, F. Beutler, D. Bizyaev, J. A. Blazek, A. S. Bolton, J. R. Brownstein, A. Burden, C.-H. Chuang, et al., MNRAS 470 , 2617 (2017), arXiv:1607.03155.\n- [99] S. Alam, M. Aubert, S. Avila, C. Balland, J. E. Bautista, M. A. Bershady, D. Bizyaev, M. R. Blanton, A. S. Bolton, J. Bovy, et al., Phys. Rev. D 103 , 083533 (2021), arXiv:2007.08991.\n- [100] E. S. Rykoff, E. Rozo, M. T. Busha, C. E. Cunha, A. Finoguenov, A. Evrard, J. Hao, B. P. Koester, A. Leauthaud, B. Nord, et al., ApJ 785 , 104 (2014), arXiv:1303.3562.\n- [101] M. Klein, S. Grandis, J. J. Mohr, M. Paulus, T. M. C. Abbott, J. Annis, S. Avila, E. Bertin, D. Brooks, E. Buckley-Geer, et al., MNRAS 488 , 739 (2019), arXiv:1812.09956.\n- [102] A. Saro, S. Bocquet, E. Rozo, B. A. Benson, J. Mohr, E. S. Rykoff, M. Soares-Santos, L. Bleem, S. Dodelson, P. Melchior, et al., MNRAS 454 , 2305 (2015), arXiv:1506.07814.\n- [103] M. Costanzi, A. Saro, S. Bocquet, T. M. C. Abbott, M. Aguena, S. Allam, A. Amara, J. Annis, S. Avila, D. Bacon, et al., Phys. Rev. D 103 , 043522 (2021), arXiv:2010.13800.\n- [104] T. McClintock, T. N. Varga, D. Gruen, E. Rozo, E. S. Rykoff, T. Shin, P. Melchior, J. DeRose, S. Seitz, J. P. Dietrich, et al., MNRAS 482 , 1352 (2019), arXiv:1805.00039.\n- [105] A. Farahi, S. L. Mulroy, A. E. Evrard, G. P. Smith, A. Finoguenov, H. Bourdin, J. E. Carlstrom, C. P. Haines, D. P. Marrone, R. Martino, et al., Nature Communications 10 , 2504 (2019), arXiv:1907.02502.\n- [106] R. Stanek, E. Rasia, A. E. Evrard, F. Pearce, and L. Gazzola, ApJ 715 , 1508 (2010), arXiv:0910.1599.\n- [107] M. Shirasaki, D. Nagai, and E. T. Lau, MNRAS 460 , 3913 (2016), arXiv:1603.08609.\n- [108] S. Grandis, J. J. Mohr, M. Costanzi, A. Saro, S. Bocquet, M. Klein, M. Aguena, S. Allam, J. Annis, B. Ansarinejad, et al., MNRAS 504 , 1253 (2021), arXiv:2101.04984.\n- [109] M. White, R. Zhou, J. DeRose, S. Ferraro, S.-F. Chen, N. Kokron, S. Bailey, D. Brooks, J. Garc'ıa-Bellido, J. Guy, et al., J. Cosmology Astropart. Phys. 2022 , 007 (2022), arXiv:2111.09898.\n- [110] T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, J. Annis, S. Avila, D. Bacon, E. Baxter, et al., Phys. Rev. D 107 , 083504  \n(2023), arXiv:2207.05766.  \n- [111] E. J. Baxter, R. Keisler, S. Dodelson, K. A. Aird, S. W. Allen, M. L. N. Ashby, M. Bautz, M. Bayliss, B. A. Benson, L. E. Bleem, et al., ApJ 806 , 247 (2015), arXiv:1412.7521.\n- [112] ' I. Zubeldia and A. Challinor, MNRAS 489 , 401 (2019), arXiv:1904.07887.\n- [113] S. Raghunathan, S. Patil, E. Baxter, B. A. Benson, L. E. Bleem, T. L. Chou, T. M. Crawford, G. P. Holder, T. McClintock, C. L. Reichardt, et al., ApJ 872 , 170 (2019), arXiv:1810.10998.\n- [114] P. S. Chaubal, C. L. Reichardt, N. Gupta, B. Ansarinejad, K. Aylor, L. Balkenhol, E. J. Baxter, F. Bianchini, B. A. Benson, L. E. Bleem, et al., ApJ 931 , 139 (2022), arXiv:2111.07491.\n- [115] M. Hilton, M. Hasselfield, C. Sif'on, N. Battaglia, S. Aiola, V. Bharadwaj, J. R. Bond, S. K. Choi, D. Crichton, R. Datta, et al., ApJS 235 , 20 (2018), arXiv:1709.05600.\n- [116] P. Predehl, R. Andritschke, V. Arefiev, V. Babyshkin, O. Batanov, W. Becker, H. Bohringer, A. Bogomolov, T. Boller, K. Borm, et al., A&A 647 , A1 (2021), arXiv:2010.03477.\n- [117] B. A. Benson, P. A. R. Ade, Z. Ahmed, S. W. Allen, K. Arnold, J. E. Austermann, A. N. Bender, L. E. Bleem, J. E. Carlstrom, C. L. Chang, et al., in Mil-  \n- limeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII , edited by W. S. Holland and J. Zmuidzinas (2014), vol. 9153 of Society of Photo-Optical Instrumentation Engineers (SPIE) Conference Series , p. 91531P, arXiv:1407.2973.\n- [118] The Simons Observatory Collaboration, P. Ade, J. Aguirre, Z. Ahmed, S. Aiola, A. Ali, D. Alonso, M. A. Alvarez, K. Arnold, P. Ashton, et al., J. Cosmology Astropart. Phys. 2019 , 056 (2019), arXiv:1808.07445.\n- [119] R. Laureijs, J. Amiaux, S. Arduini, J. L. Augu'eres, J. Brinchmann, R. Cole, M. Cropper, C. Dabin, L. Duvet, A. Ealet, et al. (2011), arXiv:1110.3193.\n- [120] LSST Science Collaboration, P. A. Abell, J. Allison, S. F. Anderson, J. R. Andrew, J. R. P. Angel, L. Armus, D. Arnett, S. J. Asztalos, T. S. Axelrod, et al. (2009), arXiv:0912.0201.\n- [121] J. Rhodes, R. C. Nichol, ' E. Aubourg, R. Bean, D. Boutigny, M. N. Bremer, P. Capak, V. Cardone, B. Carry, C. J. Conselice, et al., ApJS 233 , 21 (2017), arXiv:1710.08489.\n- [122] A. Fumagalli, M. Costanzi, A. Saro, T. Castro, and S. Borgani, A&A 682 , A148 (2024), arXiv:2310.09146.\n- [123] K. N. Abazajian, P. Adshead, Z. Ahmed, S. W. Allen, D. Alonso, K. S. Arnold, C. Baccigalupi, J. G. Bartlett, N. Battaglia, B. A. Benson, et al. (2016), arXiv:1610.02743."}