We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT3G the thirdgeneration camera on the South Pole Telescope. The surveys cover separate regions of 1500 2650 and 6000 degSUP2SUP to different depths in total observing 25 of the sky. These regions will be measured to white noise levels of roughly 2.5 9 and inlineformula mmlmath overflowscrollmmlmn12mmlmnmmlmspace width0.33emmmlmspacemmlmimmlmimmlmi mathvariantnormalKmmlmimmlmath inlineformulaarmin respectively in cosmic microwave background CMB temperature units at 150 GHz by the end of 2024. The survey also includes measurements at 95 and 220 GHz which have noise levels a factor of 1.2 and 3.5 times higher than 150 GHz respectively with each band having a polarization noise level inlineformula mmlmath overflowscrollmmlmommlmommlmsqrtmmlmrowmmlmn2mmlmnmmlmrowmmlmsqrtmmlmath inlineformula times higher than the temperature noise. We use a novel approach to obtain the covariance matrices for jointly and optimally estimated gravitational lensing potential band powers and unlensed CMB temperature and polarization band powers. We demonstrate the ability to test the CDM model via the consistency of cosmological parameters constrained independently from SPT3G and Planck data and consider the improvement in constraints on CDM extension parameters from a joint analysis of SPT3G and Planck data. The CDM cosmological parameters are typically constrained with uncertainties up to 2 times smaller with SPT3G data compared to Planck with the two data sets measuring significantly different angular scales and polarization levels providing additional tests of the standard cosmological model.
['Cosmic microwave background radiation', 'Cosmological parameters', '322', '339', 'Astrophysics - Cosmology and Nongalactic Astrophysics']
Testing the CDM Cosmological Model with Forthcoming Measurements of the Cosmic Microwave Background with SPT3G
{'Testing the ΛCDM Cosmological Model with Forthcoming Measurements of the Cosmic Microwave Background with SPT-3G': "K. Prabhu, 1 S. Raghunathan, 2 M. Millea, 3 G. P. Lynch, 1 P. A. R. Ade, 4 E. Anderes, 5 A. J. Anderson, 6, 7, 8 B. Ansarinejad, 9 M. Archipley, 10, 2 L. Balkenhol, 11 K. Benabed, 11 A. N. Bender, 12, 7, 8 B. A. Benson, 6, 7, 8 F. Bianchini, 13, 14, 15 L. E. Bleem, 12, 7 F. R. Bouchet, 11 L. Bryant, 16 E. Camphuis, 11 J. E. Carlstrom, 7, 16, 17, 12, 8 T. W. Cecil, 12 C. L. Chang, 12, 7, 8 P. Chaubal, 9 P. M. Chichura, 17, 7 A. Chokshi, 18 T.-L. Chou, 17, 7 A. Coerver, 3 T. M. Crawford, 7, 8 A. Cukierman, 13, 15, 14 C. Daley, 10 T. de Haan, 19 K. R. Dibert, 8, 7 M. A. Dobbs, 20, 21 A. Doussot, 11 D. Dutcher, 22 W. Everett, 23 C. Feng, 24 K. R. Ferguson, 25 K. Fichman, 17, 7 A. Foster, 22 S. Galli, 11 A. E. Gambrel, 7 R. W. Gardner, 16 F. Ge, 1 N. Goeckner-Wald, 14, 13 R. Gualtieri, 26 F. Guidi, 11 S. Guns, 3 N. W. Halverson, 27, 28 E. Hivon, 11 G. P. Holder, 24 W. L. Holzapfel, 3 J. C. Hood, 7 A. Hryciuk, 17, 7 N. Huang, 3 F. K'eruzor'e, 12 L. Knox, 1 M. Korman, 29 K. Kornoelje, 8, 7 C.-L. Kuo, 13, 14, 15 A. T. Lee, 3, 30 K. Levy, 9 A. E. Lowitz, 7 C. Lu, 24 A. Maniyar, 13, 14, 15 F. Menanteau, 10, 2 J. Montgomery, 20 Y. Nakato, 14 T. Natoli, 7 G. I. Noble, 31, 32 V. Novosad, 33 Y. Omori, 8, 7 S. Padin, 7, 34 Z. Pan, 12, 7, 17 P. Paschos, 16 K. A. Phadke, 10, 2 A. W. Pollak, 18 W. Quan, 17, 7 M. Rahimi, 9 A. Rahlin, 6, 7 C. L. Reichardt, 9 M. Rouble, 20 J. E. Ruhl, 29 E. Schiappucci, 9 G. Smecher, 35 J. A. Sobrin, 6, 7 A. A. Stark, 36 J. Stephen, 16 A. Suzuki, 30 C. Tandoi, 10 K. L. Thompson, 13, 14, 15 B. Thorne, 1 C. Trendafilova, 2 C. Tucker, 4 C. Umilta, 24 A. Vitrier, 11 J. D. Vieira, 10, 24, 2 Y. Wan, 10, 2 G. Wang, 12 N. Whitehorn, 37 W. L. K. Wu, 13, 15 V. Yefremenko, 12 M. R. Young, 6, 7 and J. A. Zebrowski 7, 8, 6  \n- 1 Department of Physics & Astronomy, University of California, One Shields Avenue, Davis, CA 95616, USA\n- 2 Center for AstroPhysical Surveys, National Center for Supercomputing Applications, Urbana, IL, 61801, USA 3 Department of Physics, University of California, Berkeley, CA, 94720, USA\n- 4 School of Physics and Astronomy, Cardiff University, Cardiff CF24 3YB, United Kingdom  \n5  \nDepartment of Statistics, University of California, One Shields Avenue, Davis, CA 95616, USA  \n6  \nFermi National Accelerator Laboratory, MS209, P.O. Box 500, Batavia, IL, 60510, USA  \n- 7 Kavli Institute for Cosmological Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA\n- 8 Department of Astronomy and Astrophysics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA 9 School of Physics, University of Melbourne, Parkville, VIC 3010, Australia\n- 10 Department of Astronomy, University of Illinois Urbana-Champaign, 1002 West Green Street, Urbana, IL, 61801, USA  \n11  \nSorbonne Universit'e, CNRS, UMR 7095, Institut d'Astrophysique de Paris, 98 bis bd Arago, 75014 Paris, France  \n- 12 High-Energy Physics Division, Argonne National Laboratory, 9700 South Cass Avenue., Lemont, IL, 60439, USA\n- 13 Kavli Institute for Particle Astrophysics and Cosmology, Stanford University, 452 Lomita Mall, Stanford, CA, 94305, USA\n- 14 Department of Physics, Stanford University, 382 Via Pueblo Mall, Stanford, CA, 94305, USA\n- 15 SLAC National Accelerator Laboratory, 2575 Sand Hill Road, Menlo Park, CA, 94025, USA\n- 16 Enrico Fermi Institute, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA\n- 17 Department of Physics, University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA\n- 18 University of Chicago, 5640 South Ellis Avenue, Chicago, IL, 60637, USA\n- 19 High Energy Accelerator Research Organization (KEK), Tsukuba, Ibaraki 305-0801, Japan\n- 20 Department of Physics and McGill Space Institute, McGill University, 3600 Rue University, Montreal, Quebec H3A 2T8, Canada  \n21  \nCanadian Institute for Advanced Research, CIFAR Program in Gravity and the Extreme Universe, Toronto, ON, M5G 1Z8, Canada  \n- 22 Joseph Henry Laboratories of Physics, Jadwin Hall, Princeton University, Princeton, NJ 08544, USA\n- 23 Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO, 80309, USA  \nDepartment of Physics, University of Illinois Urbana-Champaign, 1110 West Green Street, Urbana, IL, 61801, USA  \n25  \nDepartment of Physics and Astronomy, University of California, Los Angeles, CA, 90095, USA  \n- 26 Department of Physics and Astronomy, Northwestern University, 633 Clark St, Evanston, IL, 60208, USA\n- 27 CASA, Department of Astrophysical and Planetary Sciences, University of Colorado, Boulder, CO, 80309, USA 28 Department of Physics, University of Colorado, Boulder, CO, 80309, USA\n- 29 Department of Physics, Case Western Reserve University, Cleveland, OH, 44106, USA  \n30 Physics Division, Lawrence Berkeley National Laboratory, Berkeley, CA, 94720, USA  \nCorresponding author: K. Prabhu  \nkarthikprabhu22@gmail.com  \nCorresponding author: S. Raghunathan  \nsrinirag@illinois.edu  \n24  \n- 31 Dunlap Institute for Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada 32 David A. Dunlap Department of Astronomy & Astrophysics, University of Toronto, 50 St. George Street, Toronto, ON, M5S 3H4, Canada\n- 33 Materials Sciences Division, Argonne National Laboratory, 9700 South Cass Avenue, Lemont, IL, 60439, USA\n- 34 California Institute of Technology, 1200 East California Boulevard., Pasadena, CA, 91125, USA 35 Three-Speed Logic, Inc., Victoria, B.C., V8S 3Z5, Canada\n- 36 Harvard-Smithsonian Center for Astrophysics, 60 Garden Street, Cambridge, MA, 02138, USA\n- 37 Department of Physics and Astronomy, Michigan State University, East Lansing, MI 48824, USA", 'ABSTRACT': 'We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 deg 2 to different depths, in total observing 25% of the sky. These regions will be measured to white noise levels of roughly 2 . 5, 9, and 12 µ K -arcmin, respectively, in CMB temperature units at 150 GHz by the end of 2024. The survey also includes measurements at 95 and 220 GHz, which have noise levels a factor of ∼ 1.2 and 3.5 times higher than 150 GHz, respectively, with each band having a polarization noise level ∼ √ 2 times higher than the temperature noise. We use a novel approach to obtain the covariance matrices for jointly and optimally estimated gravitational lensing potential bandpowers and unlensed CMB temperature and polarization bandpowers. We demonstrate the ability to test the ΛCDM model via the consistency of cosmological parameters constrained independently from SPT-3G and Planck data, and consider the improvement in constraints on ΛCDM extension parameters from a joint analysis of SPT-3G and Planck data. The ΛCDM cosmological parameters are typically constrained with uncertainties up to ∼ 2 times smaller with SPT-3G data, compared to Planck , with the two data sets measuring significantly different angular scales and polarization levels, providing additional tests of the standard cosmological model.', '1. INTRODUCTION': "Observations of the cosmic microwave background (CMB) temperature and polarization anisotropies have played a crucial role in the field of cosmology, particularly in the establishment of a six-parameter standard cosmological model, ΛCDM. This model which provides an excellent fit to CMB data from Planck (Planck Collaboration et al. 2020), the Atacama Cosmology Telescope (ACT) (Aiola et al. 2020; Choi et al. 2020; Madhavacheril et al. 2024) and the South Pole Telescope (SPT) (Simard et al. 2018; Dutcher et al. 2021; Balkenhol et al. 2021; Pan et al. 2023), as well as a great variety of other astrophysical datasets (Alam et al. 2021; DES Collaboration 2024; Abbott et al. 2022; Kuijken et al. 2019), has been enormously successful in providing a consistent framework for understanding the universe.  \nThese successes are perhaps surprising, since critical ingredients to the model (dark matter, dark energy, and the generator of primordial fluctuations) appear disconnected from the physics we know of through laboratory experiments. The desire for clues that could deepen our understanding of these ingredients motivates us to continue our testing of the model. One particularly compelling line of testing emerges from the fact that the standard cosmological model, conditioned on data from the Planck satellite, makes extremely precise predic-  \ntions for the CMB temperature, polarization, and lensing power spectra on angular scales not well measured by Planck .  \nAdditional motivation for testing ΛCDM (beyond our ignorance about dark matter, dark energy, and the generator of primordial fluctuations) comes from cosmological tensions and anomalies as these might be evidence for beyond-ΛCDM physics. A prominent issue is the Hubble tension - a notable discrepancy between classical distance ladder determinations (Murakami et al. 2023; Riess et al. 2022) and ΛCDM-dependent determinations (Planck Collaboration et al. 2020) of the current expansion rate of the universe. An additional, though statistically weaker tension, is in the values of σ 8 , the root mean square of the density field smoothed over 8 h -1 Mpc radius spheres, derived from weak lensing observations in optical galaxy surveys and those inferred from CMB data 1 (Asgari et al. 2021; Joudaki et al. 2020; Abbott et al. 2022). There are also peculiar patterns in the CMB temperature anisotropies on large angular scales, which, by some calculations, are extremely unlikely in a ΛCDM universe (Copi et al. 2010; Givans & Kamionkowski 2023).  \nAltogether these discrepancies offer some evidence that ΛCDM is not the final word in cosmology, encouraging us to search for beyond-ΛCDM signals in our surveys. Indeed, there are alternative models that reduce the σ 8 or H 0 tensions and make different predictions for the power spectra than is the case for ΛCDM. In this context, it becomes pertinent to investigate the extent to which upcoming observational data can distinguish between the standard ΛCDM model and some of these alternatives.  \nSPT-3G, the third-generation camera on the SPT, has the potential to test the ΛCDM model by precisely mapping the primary and lensing anisotropies of the CMB. The full SPT-3G data set will comprise three surveys, which, altogether, cover 10000 deg 2 , or ∼ 25% of the sky. The SPT-3G Main survey (1500 deg 2 ), observed over five austral winters, is approaching the target depths of the next generation CMB-S4 Deep and Wide survey of 70% of the sky planned to begin next decade (Abazajian et al. 2022). Next-deepest is the SPT-3G Summer survey, a completed collection of three fields observed during four austral summers totaling 2650 deg 2 . The shallowest is the 6000 deg 2 SPT-3G Wide survey, which, after one observing season, will still be deeper than surveys from any other existing high-resolution CMB experiment of comparable survey size.  \nBesides mapping the primary CMB and lensing anisotropies, SPT-3G will provide a powerful set of data for studying a broad range of topics in cosmology and astrophysics. These include the delensing of lensinginduced B-modes (Ade et al. 2021); production of masslimited catalogs of galaxy clusters out to high redshifts using the thermal Sunyaev-Zel'dovich (tSZ) signature (Bleem et al. 2015; Huang et al. 2020; Bleem et al. 2020, 2024), and constraining cosmology with these catalogs (Bocquet et al. 2024); robust measurements of the kinematic Sunyaev-Zel'dovich (kSZ) signal for velocity reconstruction and constraining the epoch of reionization (Schiappucci et al. 2023; Raghunathan & Omori 2023; Raghunathan et al. 2024); production of mm-wave point source catalogues (Everett et al. 2020); and the detection and monitoring of mm-wave galactic, extragalactic, and solar system transient, variable, and moving objects (Whitehorn et al. 2016; Guns et al. 2021; Chichura et al. 2022; Hood et al. 2023; Tandoi et al. 2024).  \nHere we demonstrate the capacity of the SPT-3G surveys to test ΛCDM in several ways, focusing exclusively on the constraining power of CMB and lensing power spectrum measurements. We forecast how well we will measure the temperature and E -mode polarization power spectra, their cross power spectrum, and the power spectrum of the lensing potential.  \nWe also propagate the power spectrum uncertainties forward to constraints on cosmological parameters. From SPT-3G data alone we can obtain constraints on the parameters of the ΛCDM model that are comparable to, and in some cases, better than those from Planck . When compared to Planck , a greater fraction of the weight of these constraints comes from polarization (see Galli et al. 2014, for example), from smaller angular scales, and from lensing. Therefore, comparison of these parameter estimates with those from Planck will be an excellent test of the ΛCDM model.  \nThe outcome of such a test could provide evidence for physics beyond ΛCDM. Looking ahead to this scenario, we also forecast constraints on single- and doubleparameter extensions of ΛCDM from the combination of Planck and SPT-3G data and report how these constraints improve upon those from Planck alone. The test could equally well reveal consistency with ΛCDM, an outcome that should not be dismissed given the model's track record. If this is the case, joint Planck +SPT-3G ΛCDM parameter constraints will be of interest, and we provide forecasts for the combined constraints on ΛCDM as well.  \nIn addition to our comparison with Planck , we compare with forecasts we make for the high-resolution survey of the Simons Observatory (SO, Simons Observatory Collaboration 2019). SO is expected to survey roughly 40% of the sky (much more sky is available to survey from mid-latitude sites such as the SO site in the Atacama Plateau in Chile compared to the South Pole) to noise levels similar to or slightly better than the SPT-3G Summer fields, leading to constraints on ΛCDM parameters that are comparable to what we expect from the combined SPT-3G surveys.  \nFor an ongoing analysis of the first two full years of data from the SPT-3G Main survey, one of our analysis approaches is a Bayesian one, which is in principle unbiased and optimal (information lossless) despite the presence of gravitational lensing, and naturally results in joint constraints on the lensing power spectrum and unlensed CMB power spectra. To speed up estimation of the unlensed and lensing power spectra, from this likelihood, we employ the Marginal Unbiased Score Expansion (MUSE) technique (Millea & Seljak 2021). We refer to this whole inference pipeline from maps to the bandpowers and their covariance matrix as MUSE. Since MUSE can be used to obtain a covariance matrix for all the relevant spectra, it is quite convenient for us to propagate these uncertainties to a parameter error covariance matrix, and compare the results with more traditional forecasting techniques.  \nThis paper is structured as follows: In Sec. 2, we describe the specifications of the SPT-3G instrument, the survey regions and our assumptions about the noise and foregrounds. In Sec. 3 we describe the methodology used to obtain these forecasts. We report the forecasts of the parameter constraints for ΛCDM with SPT-3G alone and for extensions with SPT-3G combined with Planck in Sec. 4. In Sec. 5, we discuss prospects for SPT-3G to detect departures from ΛCDM predictions and test alternate models. Finally, we summarize our findings in Sec. 6.  \nIn this work, we set the fiducial values of the baseline ΛCDM cosmology and its extensions to Planck 2018 measurements (TT, TE, EE + lowE + lensing in Table 2 of Planck Collaboration et al. 2020) and compute the theoretical CMB power spectra using CAMB software (Lewis et al. 2000).", '2. EXPERIMENT SPECIFICATIONS': 'In this section, we briefly describe the specifications of the SPT-3G experiment along with the details of the three survey regions and a description of the noise characteristics.', '2.1. Survey Specifications': "SPT is a 10-meter diameter telescope specifically engineered for low-noise and high-resolution measurements of the millimeter-wave sky. The telescope is currently equipped with SPT-3G, the third CMB camera installed on SPT. SPT-3G is a significant upgrade compared to previous cameras, incorporating a polarization sensitive tri-chroic focal plane with almost 16,000 detectors, operating in three frequency bands, 95, 150, and 220 GHz. For more details about the instrument, we refer the readers to Benson et al. (2014); Bender et al. (2018); Anderson et al. (2018); Sobrin et al. (2018, 2022).  \nOne of the primary science goals of SPT-3G is to create a template of the gravitational-lensing-induced B modes in the patch of sky observed by the BICEPKeck (BK) family of experiments. This template can be subtracted from the BK data, reducing the primary source of variance in the BK B-mode analysis and potentially improving the sensitivity to primordial gravitational waves by a factor of three. To that end, since 2019, SPT-3G has been used primarily (eight months per year) to observe the 1500 deg 2 SPT-3G Main survey, which overlaps completely with the BICEP3 and BICEP Array sky coverage (see Fig. 1). During the austral summer season (roughly at the beginning of December), when the Sun is above the horizon and gets  \nsufficiently close to the SPT-3G Main survey region, it is picked up by the telescope far sidelobes . Therefore, until sunset, we switch to observing three fields that are safe from sun contamination. These comprise the 2650 deg 2 SPT-3G Summer survey (Fig. 1). We refer to these three fields as Summer-A, Summer-B, and Summer-C.  \nEven with the added sky coverage of the Summer survey, the constraints on ΛCDM parameters from SPT3G data are limited by sample variance. For this reason, we have paused our standard observing strategy for one year to undertake a new survey of all the sky viewable from the South Pole at instrumentally feasible observing elevations and with minimal Galactic foreground contamination. We have empirically determined that, during the winter, we are able to observe with SPT-3G between elevations of roughly 20 · and 80 · (limited by detector linearity at the low end and refrigerator performance on the high end). The intersection of these elevation limits (which, at the South Pole, correspond directly to declination limits) and the Planck GAL080 (Planck Collaboration et al. 2016b) Galactic mask defines our new SPT-3G Wide survey, which totals 6000 deg 2 and is shown in Fig. 1. Observations of the SPT3G Wide survey commenced in December of 2023 and will continue through the austral winter of 2024. We will refer to the combination of the SPT-3G Main and the SPT-3G Summer surveys as the Ext-4k survey. Similarly, we will refer to the combination of SPT-3G Main, SPT-3G Summer, and SPT-3G Wide surveys as the Ext10k survey. In Table 1, we summarize the areas and depths for the SPT-3G Main, Summer, and Wide surveys.  \nThe survey footprints for SPT-3G and other surveys are shown in Fig. 1. As is evident from the figure, all the three SPT-3G surveys have an excellent overlap with SO, which will enable a variety of cross-checks to be performed between the two surveys. We list the noise levels and the field centers for the SPT-3G surveys in Table 1. In this work, we produce forecasts both for current noise levels (2019-2023) and for the ones expected from future SPT-3G observations (2024-2026). Besides instrumental noise, our data is also affected by atmospheric noise, which increases towards large angular scales. We model the temperature noise power spectrum as  \nN ℓ = ∆ 2 T [ 1 + ( ℓ ℓ knee ,T ) α T ] , (1)  \nwhere ∆ T represents the white noise level at a given band in µ K -arcmin, ℓ knee ,T and α T are used to parameterize the atmospheric noise. For the SPT-3G Main survey, we adopt the values of ℓ knee and α in Table 2 from previous analyses, but we find that the results of  \nFigure 1. Footprints of the three SPT-3G surveys: SPT-3G Main (green), SPT-3G Summer (yellow), and SPT-3G Wide (red). Also shown are the footprints of other surveys: BICEP3 (pink dash-dotted), DES (black dash-dotted), and SO/LSST (blue dash-dotted). The background is the map of the galactic dust produced by the Planck satellite (Planck Collaboration et al. 2016a).  \n<!-- image -->  \nTable 1. Sky fraction and temperature white noise levels in the three bands for different SPT surveys considered in this work. The polarization white noise levels are expected to be ∼ √ 2 times higher.  \nour forecasting are insensitive to the exact values of ℓ knee and α . For the SPT-3G Summer survey, during which the atmospheric opacity and precipitable water vapor levels are higher, we fit real noise data to the model in Eq. 1, and the parameters in the Table 2 are the result of those fits. We also adopt the SPT-3G Summer survey parameters for the SPT-3G Wide survey (which is likely pessimistic for the higher-elevation Wide fields). We extrapolate the white noise levels for deeper integration  \ntimes. We also fit the noise model parameters for the polarization maps in a similar way. We note that the white noise levels in the polarization maps are roughly √ 2 × higher than in the temperature maps. In Table 2, we list the ℓ knee and α values adopted for different bands for both temperature and polarization. Finally, we deconvolve the experimental beam window function B ℓ from the noise spectra as N ' ℓ ≡ N ℓ B 2 ℓ where the beam window  \nFigure 2. Expected TT and EE noise power spectra for the three SPT-3G surveys: Beam-deconvolved noise curves in the three bands for SPT-3G Main (5 years), Summer (4 years) and Wide (1 year) surveys are shown as colored curves. The Minimum Variance (MV) combination of the LC residuals is the thick black solid curve in all the panels. The gray curves represent the CMB power spectra. Top panels are for temperature and the bottom panels are for polarization. The LC residuals are dominated by foregrounds for temperature for SPT-3G Main while for polarization they are roughly equal to the inverse variance weighted noise spectra from the three bands.  \n<!-- image -->  \nTable 2. Atmospheric noise specifications for the three SPT bands.  \nfunction at each frequency is estimated as detailed in Dutcher et al. (2021). We show the model noise spectra in different bands (90 GHz in green, 150 GHz in yellow, and 220 GHz in red) for all the SPT-3G surveys in Fig. 2.  \nThese spectra assume integration times of 5 years for the SPT-3G Main survey, 4 years for the SPT-3G Summersurvey, and 1 year for the SPT-3G Wide survey.", '3. FORECASTING METHODOLOGY': 'In this section, we outline the two methods used for calculating the bandpower covariance matrix from which a Fisher matrix for cosmological parameters (Tegmark et al. 1997) can be derived. The first method relies on an analytic approximation (See § 3.2) for the covariance matrix. This approach has been widely adopted in various forecasting studies, and is quite simple and flexible. Our second method employs MUSE (see § 3.3), which is one of the analysis pipelines developed for SPT-3G data. MUSE is more accurate as it readily incorporates non-ideal factors such as covariance between multipoles, covariance between different types of power spectra, and impacts of delensing. Although previous works (for example, Hotinli et al. 2022; Trendafilova 2023), suggest that these factors may not significantly alter our forecasts, adopting both these methods improves the robustness and reliability of our forecasts. We also describe the method we use to combine Planck data with SPT-3G.  \nFigure 3. S/N per multipole on the primary and lensing CMB power spectra for the SPT-3G Main (green) and the Ext-10k (orange) SPT-3G surveys, with values for Planck (black) plotted for comparison. The constraining power of the full Ext-10k survey will be larger than Planck at ℓ ≳ 1800 in TT , ℓ ≳ 800 in TE , ℓ ≳ 450 in EE and effectively for all the modes in ϕϕ . The non-shaded area corresponds to the range of multipoles used in the forecasts.  \n<!-- image -->  \nWe begin with a brief description of how we optimally combine data from multiple frequency bands.', '3.1. Linear Combination of Frequency Bands': 'We combine data from multiple frequency bands listed in Table 1 using a scale-dependent linear combination (LC) technique (Cardoso et al. 2008; Planck Collaboration et al. 2014). In harmonic space denoted by subscript ℓ , this corresponds to  \nS ℓm = N bands ∑ i =1 w i ℓ M i ℓm , (2)  \nwhere S is the desired sky signal (in our case the CMB), w ℓ is the weight, M ℓm is the spherical harmonic transform of the input map, and the subscript i runs over all the frequency bands. The multipole-dependent weights w ℓ are tuned in order to minimize the overall variance  \nfrom noise and foregrounds. They are derived as  \nw ℓ = C -1 ℓ A s A s † C -1 ℓ A s . (3)  \nHere C ℓ is a N bands × N bands matrix containing the covariance of noise and foregrounds across different frequencies at a given multipole ℓ ; A s = [1 1 .. 1] is the frequency response of the CMB in different bands and is a N bands × 1 vector. The noise modeling is described in § 2.1.  \nThe foreground signal receives contributions from both the galactic and extragalactic signals. In our baseline approach, we ignore the galactic foreground signals by introducing sky cuts but we perform tests to assess their impact using pySM simulations in § 4.1.2. The extragalactic foreground can be decomposed into signals from radio galaxies, dusty star forming galaxies (DSFG), tSZ and kSZ signals. We model these using measurements from Reichardt et al. (2021). To this end, we assume a conservative point source masking threshold, similar to Reichardt et al. (2021), and remove sources with flux above S 150 = 6 mJy. For the radio galaxies, we assume a spectral index of α radio = -0 . 76 (Reichardt et al. 2021) to model their frequency dependence. On the other hand, the DSFG contribution is split into clustering and Poisson components, which are both modeled using a modified blackbody model with T CIB = 25K, and β CIB -Clus = 2 . 2 and β CIB -Poi = 1 . 5 for the clustering and Poisson terms respectively. For the polarization power spectrum EE , we scale the above TT power spectrum of sources as C EE ℓ = P 2 C TT ℓ assuming polarization fraction of P = 3% ( P = 2%) for radio galaxies (DSFG) based on Datta et al. (2019); Gupta et al. (2019). We also remove the contribution from clusters detected at S/N ≥ 5 which roughly corresponds to M 500 c ≳ 10 14 M ⊙ (see Fig. 5 of Raghunathan 2022). For kSZ, we assume a flat D ℓ = 3 µ K 2 for all the bands. The tSZ and kSZ signals are assumed to be unpolarized, and hence we ignore the contribution from clusters on C EE ℓ and C TE ℓ .  \nWe also note here that the instrumental bandpass of the individual frequency bands are not delta functions and are roughly 25, 30, and 50 GHz for 95, 150, and 220 GHz bands respectively (Sobrin et al. 2022). The frequency responses of the detectors are calculated using a Fourier transform spectrometer. With these bandpass measurements, we calculate the effective frequencies ν eff for each of the foreground signals described above to compute their respective power spectra. For simplicity, we assume the values of ν eff to be the same as the ones reported in Reichardt et al. (2021) although the bandpasses of the current SPT-3G receiver are slightly different (see Fig. 7 of Sobrin et al. 2022). We checked  \nthe impact of the difference in the bandpasses, and it has negligible impact on the foreground modeling and the constraints reported in this work.  \nIn Fig. 2, we show the LC residuals for all three SPT3G surveys as the thick black solid curve. The top and bottom panels in the figure are for CMB temperature ( TT ) and polarization ( EE ) power spectra, denoted as gray curves. In the absence of foregrounds, the LC residuals simply correspond to the inverse variance combination of noise from different frequency bands, which is almost true for EE (bottom panels).  \nIn Fig. 3, we present signal-to-noise ( S/N ) per multipole for TT/EE/TE/ϕϕ spectra for SPT-3G Main (green) and Ext-10k (orange) surveys along with Planck (black). The power spectrum errors were calculated using the analytical formula (see Eq.5) (Knox 1995; Jungman et al. 1996; Zaldarriaga et al. 1997). Planck has a higher S/N on large scales because of the larger sky coverage. However, the S/N s of SPT-3G Main and Ext-10k surveys are much higher than Planck on small scales. In particular, the S/N s for the Ext-10k survey will sur-  \nΣ ℓ = 2 (2 ℓ +1) f sky       ( ˜ C TT ℓ ) 2 ( C TE ℓ ) 2 ( C TE ℓ ) 2 ( ˜ C EE ℓ ) 2 ˜ C TT ℓ C TE ℓ ˜ C EE ℓ C TE ℓ ( C Tϕ ℓ ) 2 ( C Eϕ ℓ ) 2  \n˜ C TT ℓ C TE ℓ ( C Tϕ ℓ ) 2 ˜ C EE ℓ C TE ℓ ( C Eϕ ℓ ) 2 1 2 [ ( C TE ℓ ) 2 + ˜ C TT ℓ ˜ C EE ℓ ] C Tϕ ℓ C Eϕ ℓ C Tϕ ℓ C Eϕ ℓ ( ˜ C ϕϕ ℓ ) 2       , (5)  \nwhere f sky is the fraction of the sky for the survey under consideration and  \n˜ C X ℓ = C X ℓ + N X ℓ , (6)  \nX = { TT,EE,ϕϕ } . N TT ℓ and N EE ℓ are the residual noise and foreground power after combining data from all the bands using the LC method. N ϕϕ ℓ is the lensing reconstruction noise derived using Quadratic Lensing Estimator (Hu & Okamoto 2002). In Table 3 we list the fiducial values of the cosmological parameters at which the Fisher matrix is evaluated and associated priors.  \nWe calculate the Fisher matrix for each of the SPT3G surveys given in Table 1 and combine them to obtain the final SPT-3G Fisher matrices as  \nF Summer = F Summer -A + F Summer -B + F Summer -C F Ext -4k = F Main + F Summer F Ext -10k = F Main + F Summer + F Wide . (7)  \nThe parameter covariance matrices are then obtained by inverting the Fisher matrix. We fix the ℓ ranges for all the SPT-3G surveys to the following for the Fisher matrix calculation: 300 ≤ ℓ ≤ 3500 for TT/TE ; 300 ≤ ℓ ≤ 4000 for EE ; and 30 ≤ L ≤ 3500 for ϕϕ . For lensing reconstruction, we use 300 ≤ ℓ ≤ 3500 for T and 300 ≤ ℓ ≤ 4000 for P. We ignore ℓ > 3500 in TT/TE because of the difficulties in modeling the extragalactic foreground signals in the temperature data. We have checked the correlations between the SPT-3G Main and the contiguous SPT-3G Summer surveys and we have found that they are negligible. This is due to the fact that we do not include the large-scale modes ℓ ≤ 300. Since the extragalactic foregrounds are largely unpolarized (Datta et al. 2019; Gupta et al. 2019), we can set a higher ℓ max for EE . We find that extending ℓ EE max > 4000 results in marginal improvements, indicating those modes are noise dominated, as can be inferred from Fig. 3.  \n3.3. MUSE  \npass Planck beyond ℓ = 1800, 800, and 450 for TT , TE , and EE respectively. For ϕϕ , the S/N are better than Planck on all scales. As a result, we can expect significant improvement on cosmological parameters that control structure formation and that are sensitive to the damping tail from SPT-3G compared to Planck as we show later in this work.', '3.2. Analytic method': "Given a set of power spectra of the lensed CMB and the lensing potential, X,Y ∈ [ TT,EE,TE,ϕϕ ], the Fisher matrix is given by,  \nF ij = ∑ X,Y ∑ ℓ ∂C X ℓ ∂θ i · ( Σ -1 ℓ ) XY · ∂C Y ℓ ∂θ j , (4)  \nwhere θ s correspond to the cosmological parameters being constrained. The partial derivatives of the spectra with respect to the parameters θ , denoted as ∂C ℓ /∂θ j , are obtained using the finite difference method. Σ ℓ is the covariance of the power spectra, given by Eq.(5)  \nTable 3. Fiducial values of the parameters and priors used in this work. All the applied priors are Gaussian with the widths given below in the table.  \nExtensions:  \nMUSE (Millea & Seljak 2021) is a general algorithm for approximate marginalization over arbitrary latent parameters which yields Gaussianized, asymptotically unbiased and near-optimal constraints on parameters of interest. In our application of MUSE, the parameters of interest are the unlensed TT , TE , EE , and ϕϕ bandpowers, which are weighted averages of the respective power spectra across defined multipole bins. These bandpowers depend on the latent parameters, which are the pixel values of the lensing map and the unlensed CMB field maps.  \nAs previously demonstrated by Millea et al. (2021) using SPTpol data, the standard Quadratic Lensing Estimator is sub-optimal for lensing reconstruction at noise levels below ∆ T ∼ 5 µ K -arcmin (Hirata & Seljak 2003a,b) and hence, using MUSE has an advantage in terms of the lensing reconstruction. Another advantage of MUSE is that it performs optimal delensing on the lensed CMB maps and produces estimates of the unlensed power spectra, lensing potential power spectra and their joint covariance. This covariance naturally includes lensing-induced correlations between different bandpowers ( ℓ -ℓ ' ); and also between delensed CMB  \nTT/EE/TE and lensing ϕϕ spectra (Peloton et al. 2017; Trendafilova 2023).  \nDespite the above advantages, we do not expect a significant improvement in the constraints from using MUSE over QE methods. This is because the lensing bandpower errors are dominated by sample variance over the angular scales for which MUSE reduces ϕ map noise, and the lensing-induced correlations are not important for the current noise levels. Moreover, the optimality of MUSE is non-negligibly better only for the SPT-3G Main survey, and not for the others given their higher noise levels (see Table 1).  \nThe optimality of the MUSE estimates is useful for other purposes like delensing the lensing-induced Bmodes and cross correlations, at least for the SPT-3G Main survey with its low noise levels (see Table 1). We demonstrate this in Fig. 4 which shows the lensing reconstruction noise N (0 ,ϕϕ ) L for different surveys. For the SPT-3G Main survey shown in green, the lensing reconstruction noise goes down by roughly × 2 . 5 ( × 1 . 9) for polarization-only (including temperature information) lensing estimators when we replace QE (solid) with MUSE (dash-dotted). As mentioned above, the difference between QE and MUSE is negligible for other surveys and hence, not shown.  \nIn MUSE, we express the lensing and unlensed CMB power spectra as scaled versions of the fiducial power spectra, which, for our purposes, are based on the Planck 2018 spectra (Planck Collaboration et al. 2020). This is given by  \nC ℓ = A ℓ C fid ℓ , (8)  \nwhere C fid ℓ represents the fiducial power spectrum, and A ℓ is the scaling factor applied to this fiducial spectrum. Here, C ℓ is used as a general notation that encompasses all types of power spectra, including unlensed TT/TE/EE and ϕϕ .  \nIn our application of MUSE, we do not estimate the power at each individual multipole but instead bandpower parameters, A b , that each scale a fiducial power spectrum across a band of width ∆ ℓ = 100. As a result, MUSE produces a covariance matrix of these scaling coefficients Σ A b ≡ ⟨ A b A ∗ b ⟩ . The A b can be understood as the binned version of A ℓ .  \nAs described below, this covariance is subsequently utilized to construct the cosmological parameter-Fisher matrix, from which we derive the constraints on parameters. We provide more details about MUSE in the appendices. Specifically, we describe the derivation of the MUSE covariance matrix in Appendix A, the simulations used for estimating the covariance matrix in Appendix A.1, and the lensing reconstruction noise in Appendix A.2.  \n0  \nFigure 4. Left panel: Lensing reconstruction noise N (0 ,ϕϕ ) L power spectra for the SPT-3G surveys along with Planck , SO and ACT. For the SPT-3G Main survey, we show two curves: solid green for QE and dash-dotted green for MUSE. For SPT-3G Main survey noise levels, lensing noise from MUSE is roughly × 2 better than QE at large scales. For other surveys, the difference between MUSE and QE is negligible and hence, not shown. Right panel: Bandpower errors for the lensing power spectrum (gray) from different surveys. While the lensing reconstruction noise N (0 ,ϕϕ ) L is lower on large scales for MUSE, the bandpower errors are not significantly different as they are dominated by the sample variance.  \n<!-- image -->  \nThe standard Fisher formula described in the previous section Eq.(4) can be adapted to incorporate the MUSE covariance matrix Σ A b . Let the covariance matrix of the power spectra be represented by Σ C ℓ and the covariance matrix of the amplitudes be represented by Σ A ℓ , then the Fisher matrix can be written as,  \nF ij = ∑ XY ∑ ℓℓ ' ∂C X ℓ ∂θ i ( Σ -1 C ℓ ) ℓℓ ' XY ∂C Y ℓ ' ∂θ j (9) = 1 C fid ℓ ∂C ℓ ∂θ i · ( Σ -1 A ℓ ) · 1 C fid ℓ ∂C ℓ ∂θ j (10) ≈ ∂A b ∂θ i · ( Σ -1 A b ) · ∂A b ∂θ j . (11)  \nFor ease of notation, we drop the summation symbols and transition to using matrix products in the second line. In the last line we assume that the Fisher matrix is approximately equivalent to its binned version. We obtain the derivatives of A b with respect to the cosmological parameters via a minimum-variance binning of ∂C ℓ /∂θ given by  \n∂A b ∂θ = 1 ∑ ℓ ∈ b 1 σ 2 A ℓ ∑ ℓ ∈ b 1 σ 2 A ℓ ∂C ℓ /∂θ C fid ℓ , (12)  \nwhere we use an analytic ansatz for the weighting function σ 2 A ℓ . For X ∈ [ TT,EE,ϕϕ ] this is given by  \nσ 2 A X ℓ = σ 2 C X ℓ ( C X, fid ℓ ) 2 ≡ 2 2 ℓ +1 ( C X ℓ + N X ℓ ) 2 ( C X, fid ℓ ) 2 (13)  \nand for TE it is  \nσ 2 A TE ℓ = 1 2 ℓ +1 ( C TE ℓ ) 2 +( C TT ℓ + N TT ℓ )( C EE ℓ + N EE ℓ ) ( C TE, fid ℓ ) 2 . (14)  \nIn Fig. 5, we show the ratio of parameter constraints obtained from MUSE and the analytic covariance for SPT-3G Main (green) and the SPT-3G Wide (orange) surveys. We find that these methods are consistent with each other within a 10% margin, which we deem sufficient to satisfy our validation criteria.  \nThere are a number of reasons to expect these errors to differ to some degree, which we already brought up when describing the advantages of MUSE. Another reason is that the analytic results are not binned, but we found that for both the Wide and Main surveys, binning to ∆ ℓ = 100 only leads to at most 2% degradation in constraints. The largest difference is the 8% tighter constraint on Ω b h 2 from the analytic forecast for the SPT-3G Wide survey. This is because the analytic estimate does not include the ℓ -byℓ lensing-induced corre-  \nFigure 5. Ratio of the ΛCDM parameter constraints derived from MUSE to those derived from analytical covariance matrix for the SPT-3G Main survey (dark green circle) and the SPT-3G Wide survey (orange triangle). We note that they agree with each other within 10%.  \n<!-- image -->  \nong different CMB multipoles, and that can lead to overly-optimistic constraints. By artificially removing off-diagonal elements from the MUSE covariance matrix 2 , we were able to largely reproduce this reduced analytic error for both Ω b h 2 and n s .", '3.4. Combining with Planck data': "Comparison of the SPT-only ΛCDM constraints with those from Planck , as we will see, will enable a powerful test of the ΛCDM model. If these tests are consistent with ΛCDM then we will be interested in using the combined SPT-3G+ Planck data to estimate ΛCDM parameters, as Planck provides complementary information on large angular scales. Consistent or not, we will also be interested in using the combined data to constrain ΛCDM extensions. Hence, we report the expected constraints from the combined datasets under the assumption of ΛCDM as well as under the assumption of various extensions.  \nThe standard approach to obtain this combination is to add the Fisher matrices from SPT-3G and Planck similar to Eq.(7) as  \nF SPT+ Planck = F ̂ SPT + F ̂ Planck , (15)  \nwhere F ̂ SPT is the modified SPT-3G Fisher matrix obtained by taking the inverse variance combination of SPT-3G and Planck noise in the field under consideration; and F ̂ Planck is the modified Planck Fisher matrix after scaling to remove the overlapping sky region be-  \nFigure 6. Ratio of power spectrum errors for Planck to those of SPT-3G for the same sky fraction. ℓ min , the multipole value at which we switch from Planck to SPT-3G in combined constraints (shown as the vertical dash-dotted lines), is chosen to be where this ratio roughly crosses 2. The blue, green, and orange colors correspond to TT,TE,EE respectively. For ϕϕ this ratio is greater than 2 for all the multipoles used in the forecast, and therefore it is not shown. The darker curves correspond to the SPT-3G Main survey whereas the lighter curves represent the Ext-10k survey. The choice of ℓ min does not change significantly for different SPT3G surveys.  \n<!-- image -->  \neen SPT-3G and Planck . For example, if f sky = 1 for F Planck , then for F ̂ Planck , we use f sky = 1 -f SPT sky .  \nHowever, obtaining Eq.(15) is non-trivial as it requires us to know the sky fraction used by Planck , which varies for different bands and power spectra (Planck Collaboration et al. 2020). Instead, we follow a different approach which involves two steps: First, we derive the Planck Fisher matrix directly from the publicly available Planck PR3 chains (Planck Collaboration et al. 2020); second, we combine the above Planck Fisher matrix with SPT3G Fisher matrix, obtained using the analytic method. We modify the SPT-3G Fisher matrix by removing the large-scale modes from it, to avoid double counting information in the overlapping region.  \nWe use BASE PLIKHM TTTEEE lowL lowE lensing 3 chains for the first step unless specified otherwise. These chains, which contain samples from the probability distribution of cosmological parameters, are used to construct a parameter covariance matrix. By inverting this covariance matrix, we obtain the Planck Fisher matrix. For the second step, we restrict the information we use from SPT power spectra to multipoles above some ℓ min or L min threshold, set by the approximate multipole at which the S/N for Planck is × 2 lower than that for SPT  \nin the same patch. We find that the choice of ℓ min does not change significantly for different SPT-3G surveys. As can be seen from Fig. 6, which shows the ratio of Planck TT , TE , and EE power spectrum uncertainties to those projected for SPT-3G, this threshold is roughly 2000, 1000, and 750 for TT , TE , and EE respectively, and we adopt these values for the joint forecasting. For lensing power spectra, we set L min = 30 (the same as our fiducial setting) because Planck 's lensing S/N is much lower than SPT on all scales, as can be seen from Fig. 3. However, we find our forecasts to be insensitive to the choice of L min . For example, with L min = 400 for SPT, our constraints only weaken by ≲ 5%.  \nWe validate this approach of using ℓ min cuts rather than the f sky scaling in Eq.(15) by creating a toy model for a Planck -like Fisher matrix with f sky = 0 . 5. For this toy model, we use the following multipole ranges: ℓ ∈ [2 , 2000] , [30 , 2000] , [30 , 2000] respectively for TT , TE , and EE and L ∈ [30 , 400] for ϕϕ . Next, we compare the Fisher forecasts from Eq.(15) and the ℓ min -cut approach discussed above. We find the results from the ℓ min cut approach to be ≲ 8% worse than the traditional method. Given that the two methods agree well and given that our ℓ min -cut-based forecasts are slightly on the conservative side, we choose to use this approach for the rest of this work.", '4. RESULTS': 'In this section, we begin by presenting our forecasts on cosmological parameters using all three SPT-3G surveys, assuming the 6-parameter ΛCDM model. Subsequently, we project the extent to which combined SPT-3G and Planck data can constrain both singleand double-parameter extensions to the ΛCDM model. We also qualitatively investigate how well SPT-3G data can differentiate between ΛCDM and alternative models that have been proposed in the literature to address current cosmological tensions (for example, Cyr-Racine et al. 2022; Schoneberg et al. 2022; Meiers et al. 2023; Hughes et al. 2024; Khalife et al. 2024; Schoneberg et al. 2022). The models we consider are fully consistent with current CMB and lensing power spectra measurements but make different predictions, from those of ΛCDM, for signals that can be measured well by SPT-3G.', '4.1. Forecasts of constraints on ΛCDM parameters': "We start first with forecasts of constraints on the six parameters of the ΛCDM model, namely baryon density Ω b h 2 , cold dark matter density Ω c h 2 , Hubble constant H 0 , optical depth to reionization τ re , amplitude (defined at the pivot scale k 0 = 0 . 05 Mpc -1 ) and spectral index of scalar fluctuations, A s and n s . For these constraints  \nFigure 7. Figure of merit (FoM) for ΛCDM compared to Planck : We present the FoM for SPT-3G Main (green), Ext4k (yellow), and Ext-10k (red). The transparent yellow circle represents the improvement in Ext-4k, if we were to continue the SPT-3G Main survey during 2024 as opposed to conducting the SPT-3G Wide survey. The evident increase in FoM with the new strategy emphasizes its potential efficacy over continuing the SPT-3G Main survey for cosmological parameter constraints. For reference, we also show the FoM expected for the two SO configurations (Baseline in blue and Goal in orange) under the assumption that the observations will begin in early 2025.  \n<!-- image -->  \nFigure 8. Constraints on the standard ΛCDM parameters expected from combinations of the three SPT-3G surveys (SPT-3G Main in green, Ext-4k in yellow, and Ext-10k in red) compared to Planck (grey shaded region). We expect significant improvements on all the parameters except n s and τ re which depend on the largest scales that ground-based surveys cannot easily measure.  \n<!-- image -->  \nwe fix the neutrino mass to the minimal value of 0 . 06 eV. We also include a prior on σ ( τ re ) = 0 . 007 from Planck (Planck Collaboration et al. 2020). We use the MUSEapproach as described in § 3.3 to derive these constraints, and remind the reader that there is agreement between MUSE and the traditional analytic approach to better than 10%, as shown in Fig. 5.  \nAs we will see below, SPT-3G can produce constraints on ΛCDM parameters that are tighter than those we have from Planck . Since these constraints are obtained from largely different signals than those from Planck , they provide a robust consistency test for the ΛCDM model. As mentioned earlier in § 3.4, we also present ΛCDM forecasts by combining SPT-3G and Planck .  \nIn Fig. 7, we present the projected Figure of Merit (FoM), which is inversely proportional to the volume of allowed 6-dimensional ΛCDM parameter space. It is calculated as  \nFoM = √ 1 / det( F -1 ) (16)  \nfor different SPT-3G surveys (SPT-3G Main in green, Ext-4k in yellow, and Ext-10k in orange-red) along with SO (blue). The FoM is inversely proportional to the 6dimensional volume of the 68% confidence region. 4 We plot it divided by the FoM for Planck .  \nAs is evident from the figure, SPT-3G surveys significantly increase the FoM relative to that from Planck . From just the first two years (2019-2020) of SPT-3G Main survey data, the FoM increases by × 2 . 5 over Planck , and that jumps to × 13 with the addition of the first two years of SPT-3G Summer data. Both of these are further improved by roughly × 1 . 7 with the inclusion of SPT-3G Main and SPT-3G Summer survey data through 2023.  \nThe figure also demonstrates the benefit of expanding sky coverage rather than continuing to go deeper on the SPT-3G Main survey patch. One can see this from the minimal improvement to Ext-4k ( × 1 . 1) that we anticipate if we continue with the SPT-3G Main survey patch observations for the 2024 winter, indicated by the semi-transparent yellow circle. On the other hand, switching to the SPT-3G Wide observations in the 2024 season improves the constraining power by × 5 compared to Ext-4k and by × 130 compared to Planck .  \nWere we to observe the SPT-3G Wide patch for one more year we would not notice a significant increase in the FoM. Instead of doing this we plan to switch back to the SPT-3G Main survey during the 2025 and 2026  \nwinter seasons. Further noise reduction in the SPT-3G Main survey will not increase the FoM much, as can be seen in Fig. 7, but will significantly improve the measurements of the secondary CMB anisotropies including lensing, which will further improve the subtraction of lensed B-modes and hence, the search for B-modes from primordial gravitational waves (Ade et al. 2021).  \nIn Table 4 we list our forecasts of ΛCDM parameter errors for different SPT-3G datasets. We also present them in Fig. 8, where the bars show the ratios of the values of σ ( θ ) from Planck to those from SPT.  \nLet us begin by noting two things relevant for interpretation of our A s and τ re results. First, the calibration of SPT-3G maps relies on data from Planck (Dutcher et al. 2021), and, consequently, SPT-3G constraints on the normalization parameter A s are not entirely independent of Planck . Second, for SPT-only constraints in this work, we include a Planck -like prior on the optical depth, σ ( τ re ) = 0 . 007, since τ re is primarily constrained by the large-scale ℓ ≲ 10 reionization bump in EE/TE , not accessible by SPT. All the τ re constraints are dominated by this prior. These two points are not unrelated since the combination A s e -2 τ re is very tightly constrained by CMB data at ℓ > 20, so improvements in τ re determination translate into improvements in A s determination, and vice versa. The small improvements that we do see in τ re are due to the dependence of the amplitude of gravitational lensing on A s , an amplitude unaffected by τ re .  \nFocusing on the other four parameters, we find that the SPT-only constraints are generally comparable to, or better than, those obtained from Planck alone, for each of our three SPT-3G survey combinations. For the case of Ext-10k, the improvements in H 0 and Ω b h 2 are nearly a factor of 2, and for Ω c h 2 they are just more than a factor of 1.5. For n s , the relative errors from SPT-3G Main survey and Ext-4k survey are larger compared to those from Planck due to Planck 's significantly broader multipole range. Nonetheless, the Ext-10k survey still manages to achieve slightly tighter constraints on n s than what we have from Planck .", '4.1.1. Comparison with SO': "In Fig. 7, we also show the expected constraining power from the Simons Observatory for two configurations: SO-Baseline (blue dash-dotted) and SO-Goal (orange dash-dotted), based on the noise levels in Simons Observatory Collaboration (2019). We use the original  \nTable 4. Forecasted 1 σ errors on ΛCDM parameters for various survey regions using SPT data. The numbers in the parentheses correspond to the constraints obtained when SPT-3G is combined with Planck . Note that when SPT-3G is combined with Planck , we remove information below ℓ min = 2000 , 1000 , 750 , 30 for TT,TE,EE and ϕϕ power spectra respectively from SPT-3G, to avoid double counting information.  \nSO noise model code 5 to get the noise curves for the individual frequency bands and combine information from all the bands using the LC technique described in § 3.1. We use the same multipole ranges described in § 3.2 for forecasting constraints from SO. As a check, we compare the constraints that we get from our forecasts for SO against what is reported in Simons Observatory Collaboration (2019) but note that we do not include Planck unlike SO's original forecasts.  \nComparing our forecasts for SO to Table 3 of Simons Observatory Collaboration (2019), our forecasts are slightly better for H 0 and Ω c h 2 but they agree within 10%. For n s and Ω b h 2 , we obtain 15-35% worse constraints, which is presumably because we do not add the large-scale information from Planck . On the other hand, we obtain roughly 35% (20%) better constraints on A s ( τ re ). When we increase the prior uncertainty on optical depth to σ ( τ re ) = 0 . 009 to match SO's forecasting code (see caption of Table 3 in Simons Observatory Collaboration 2019), our results on A s and τ re match the numbers quoted in Simons Observatory Collaboration (2019). Switching to SO-Goal configuration reduces the errors on all parameters by 5% except for Ω c h 2 where the constraining power improves by 15%.  \nWe also note from Fig. 7, the constraining power of the SPT-3G Ext-10k survey for ΛCDM parameters will be roughly similar to the first two years of SO-Baseline. Such similar levels of constraining power between SPT and SO along with significant sky overlap (see Fig. 1) will allow for comparisons of the two datasets at both the power spectrum and map level (such as the ACTPlanck and SPTPlanck comparisons shown in Louis et al. 2014 and Hou et al. 2018) which will provide useful opportu-  \nities for identifying and/or limiting systematic errors in individual datasets.", '4.1.2. Real-world effects': 'In this section, we examine how our forecasts change as we vary our baseline modeling assumptions, and introduce additional non-ideal aspects into our model of the data on which the forecasts are based. We modify our noise curves to account for filtering effects and to include a higher level of complexity in the noise due to the correlated atmospheric signal on large scales in the temperature data. We also investigate the degradation in constraints that could come from galactic foregrounds. We also check the impact of including uncertainties in the calibration of our temperature and polarization maps. As detailed below, we find the degradation of parameter errors relative to the baseline case, from all of these effects, is ≲ 10%. Finally, we build further confidence in our forecasts by comparing them to constraints from an ongoing MUSE analysis of the two-year SPT-3G Main data.  \nAtmospheric noise -The first thing we examine is the impact of the filtering applied to the detector time-ordered data (TOD) before it is processed into maps. To avoid excess noise in the scan direction in the resulting maps, primarily as a result of the large-scale atmospheric noise, we apply filters to remove the long-timescale information from the detector TOD. While this improves the quality of the maps, it results in a reduction of data volume that increases the final power spectrum errors by a factor equivalent to raising map noise levels by roughly 10%. We test the impact of this by increasing the white noise levels in Table 1 by 10%. We found that this results in changes to the parameter constraints at the ≤ 1% level. Secondly, since SPT-3G uses multichroic detectors that can simultaneously observe in multiple frequencies, the large-scale atmospheric noise is expected to be corre-  \nted between bands. From the 2019-20 dataset, we estimate the correlation at the largest scales to be ∼ 90% between the 150 and 220 GHz bands, and around 70% between the 90 and 150 GHz bands as well as between the 90 and 220 GHz bands. This degree of correlation diminishes as we move to small scales. To take this into account, we set the cross-frequency noise to be correlated across bands and rerun the LC step to get the residuals. This correlated noise, in fact, reduces the large-scale LC residuals slightly. When we propagate these to parameter constraints, we find changes to the parameter constraints at the ≤ 1% level, indicating that our forecasts are robust to these details about the noise spectra.  \nGalactic foregrounds and mask -Given that the SPT-3G Wide survey extends close to the plane of the galaxy (see Fig. 1), it is crucial to understand the impact of galactic foregrounds on the parameter constraints. To assess this, we make use of the PySM simulations (Thorne et al. 2017) of galactic foregrounds and the galactic foreground masks ( GAL070, GAL080, GAL090 ) produced by the Planck collaboration (Planck Collaboration et al. 2016b). The names of the above masks correspond to the sky retained after masking the galactic plane. For instance, GAL090 masks 10% of the region with a high level of foregrounds close to the galactic plane while GAL070 removes 30% of the sky and is the most conservative amongst the three.  \nWe use these masks and calculate the level of galactic dust and synchrotron signals in PySM simulations in the following regions: Case (A) GAL080-GAL070 ; Case (B) GAL090-GAL070 ; and Case (C) GAL090-GAL080 . The three cases correspond to the differences between the respective masks; for instance, GAL080-GAL070 represents the sky area included in mask GAL080 but excluded in GAL070 . As expected, we find the galactic foreground in Case (C) to be the highest followed by Cases (B) and (A).  \nWe perform a conservative test by assuming that the galactic foreground levels across the entire SPT-3G Wide survey match one of these three cases. With this assumption, we derive the LC residuals and the subsequent impact on the parameter constraints. Even for Case (C), where the galactic foregrounds are the highest, we find the constraints degrade by only 8%. This is not surprising given that the galactic foregrounds are mostly important on the largest scales and the constraining power of SPT-3G is predominately driven by the highℓ region.  \nIn the preceding sections, we mentioned that the observing strategy to minimize galactic foreground contamination in the SPT-3G Wide involves observing  \nwithin the confines of the Planck GAL080 mask, which corresponds to roughly 6000 deg 2 of sky coverage. In the scenario where we might need to adopt a more conservative Planck GAL070 mask, we found that this would lead to a modest degradation in parameter constraints, at 4-5%. While we do not address potential biases due to mismodeling of the galactic foregrounds, we remain confident that any such biases would be minimal, given that the levels of galactic foregrounds are significantly lower than the CMB power spectra levels at the multipoles relevant to this analysis.  \nTemperature and polarization calibration -Next, we assess the impact of calibration uncertainty on the temperature and polarization maps by including two absolute calibration parameters, T cal and P cal (as defined in Dutcher et al. 2021), in our forecasts and marginalizing over them. In the case where we do not apply a prior, we find a significant × 2 hit on some of the parameters, which is not surprising since parameters like A s are fully degenerate with the calibration parameters. By comparing SPT-3G Main survey maps from 20192020 to Planck maps in that same region, we are able to constrain T cal and P cal to roughly 0.2% and 1%, respectively (where the errors are a combination of statistical uncertainty in the comparison and uncertainties in the Planck calibration). If we adopt these as priors, the resulting degradation in the parameter constraints are ≤ 5% for all parameters. If we loosen the T cal prior to 1%, the degradation increases to ∼ 20%. If we assume either parameter is estimated perfectly and allow for uncertainty in only one calibration parameter, we find the resulting degradation to be negligible, even without any prior, consistent with the results from Galli et al. (2021).  \nFinally, we replace the covariance matrix for the 2year SPT-3G Main survey from this work with the one from an ongoing MUSE analysis performed on the real data from the same survey. The latter takes into account certain effects that are present in the real data but not taken into account in our forecasts. These include anisotropic noise, anisotropic filtering, and masking. With this covariance swap, we find that the degradation in the parameter constraints is ≤ 10%.', '4.2. Single-parameter extensions to ΛCDM': 'Now we turn our attention towards extensions to the ΛCDM model. For extensions, we always combine SPT3G with Planck , and use the analytic method. We consider 7 different cases by varying the following parameters one at a time along with the 6 ΛCDM parameters: dark-energy equation of state ( w 0 ), Helium mass fraction ( Y P ), running of the spectral index ( n run ≡ d n s / dln k ), sum of neutrino masses ( ∑ m ν ), effective  \nFigure 9. Ratio of 1 σ constraints on single-parameter extensions from Planck to those forecasted from the combination of Ext-10k and Planck . The gray shaded region corresponds to the current best constraints from Planck alone. The actual errors are expressed in the text boxes.  \n<!-- image -->  \n+  \nFigure 10. Degradation in the ΛCDM constraints due to single-parameter extensions compared to the original constraints reported in Fig. 8 for the Ext-10k+ Planck dataset. The horizontal axis shows the parameter that is varied in our extensions. Except for H 0 , which is affected by a geometric parameter degeneracy, the degradation in constraints for other ΛCDM parameters is quite mild.  \n<!-- image -->  \n∑  \nnumber of neutrino species ( N eff ) and the spatial curvature parameter (Ω k ). The forecasted constraints for Ext-10k are reported in Fig. 9. As before, we show the improvement in constraints compared to Planck using the bars. The actual parameter constraints are given in the text boxes. The constraints from other SPT-3G surveys are presented in Table 5 of Appendix B.  \nWe find a significant ( × 2 -× 3) reduction of the errors in many of these parameters, namely Ω k , w 0 , Y P , and  \nN eff . For other parameters, namely ∑ m ν , and n run , the improvement is × 1 . 4 -× 1 . 7 compared to Planck . Some of these improvements are driven by the inclusion of lensing (Ω k , w 0 ) while others come from higher-redshift information in the damping tail of the CMB, whose extraction benefits from wide sky coverage ( n run , Y P , and N eff ).  \nWe also investigate the degradation in ΛCDM constraints due to degeneracies with the extension parameters. We present these error degradations in Fig. 10 for Ext-10k+ Planck ; the full list can be found in Table 6 of Appendix B. As can be inferred from the figure, the constraints on the Hubble constant H 0 degrade significantly when we free θ ∈ [ ∑ m ν , N eff , Ω k , w 0 ]. Note that these degradations largely disappear with the addition of external datasets that constrain the redshift-distance relationship down to z ≲ 0 . 1, namely supernova brightness measurements. This is the reason why these extensions are not viable solutions to the H 0 tension (Knox & Millea 2020). The degradation is less severe ( ∼ × 1 . 6) for n s , Ω b h 2 when we vary θ ∈ [ N eff , Y P ]. For the rest, the constraints are robust and the degradation is ≲ × 1 . 3.  \nIn expanding the parameter space from ΛCDM to ΛCDM+ w 0 , we notice an unexpected tightening in the constraints for parameters Ω b h 2 , Ω c h 2 , and n s . We suspect that this arises from our use of a Gaussian approximation for a posterior that is significantly non-Gaussian. The non-Gaussianity is visually evident in the one and two-dimensional marginal posteriors derived from the appropriate Planck chains. With the addition of an H 0 prior, the range of acceptable w values tightens significantly, reducing the degree of non-Gaussianity, and this anomalous behavior of our forecasts no longer appears.', '4.3. Double-parameter extensions to ΛCDM': 'We also explore two-parameter extensions similar to those considered in Planck Collaboration et al. (2020). These are (a) N eff and Y P , (b) N eff and n run , and (c) N eff and ∑ m ν . In Fig. 11, we present the marginalized twodimensional contours for the above extensions from both Planck (black) and Ext-10k+ Planck (orange). The FoM for the two-parameter extensions, calculated similar to Fig. 7, is given in the legend. In most cases, we find an improvement of greater than a factor of 2 relative to Planck alone.  \nIn case (a), we vary the effective number of neutrino species N eff and the fraction of baryonic mass in helium Y P simultaneously. Such constraints are useful for constraining scenarios in which the value of N eff at big bang nucleosynthesis differs from that at recombination (Millea et al. 2015; Cyburt et al. 2016; Particle Data Group et al. 2020). Since both these parameters af-  \n∑ m ν Figure 11. Marginalized two-dimensional parameter contours for the three pairs of double-parameter extensions to ΛCDM considered in this work. The red contours represent the joint constraints from the Ext-10k+ Planck , and the black contours show the constraints from Planck alone. The inverse of the allowed volume of the two-dimensional parameter space (FoM) are given in the legend. Compared to the constraints from Planck alone, we find reductions in individual parameter errors that approach or exceed a factor of two except for n run , where we obtain × 1 . 3 reduction.  \n<!-- image -->  \nfect the damping tail of the CMB, they are partially degenerate and, as expected, the constraints on both parameters degrade by more than × 2 compared to the single-parameter extensions in Fig. 9.  \nIn case (b), we vary N eff and n run which are also degenerate to some level leading to × 1 . 5 weaker constraints compared to Fig. 9.  \nIn case (c), we vary N eff and the sum of neutrino masses ∑ m ν since one can expect some level of correlation between the two (Planck Collaboration et al. 2020). However, this correlation depends on the datasets being considered. In our case, ∑ m ν is primarily constrained by lensing while N eff constraints are dominated by the small-scale TT , EE , and TE power spectra. Subsequently, the degradation in the constraints on these two parameters is negligible compared to what we obtain for the single-parameter extensions in Fig. 9.', '5. POWER SPECTRUM PREDICTIONS FROM ΛCDM AND ALTERNATIVES': "In this section we discuss, in a qualitative manner, prospects for using SPT-3G data to detect departures from ΛCDM predictions for CMB power spectra conditioned on Planck data. Many alternative models that resolve cosmological tensions, such as the H 0 and σ 8 tensions, predict power spectra that differ from ΛCDM predictions at small angular scales, and we use a selection of such models to guide our discussion. In Fig. 12, we present the ΛCDM predictions for angular power spectra given Planck data, the predictions of alternative models  \nthat fit the Planck data well while either having a high H 0 or a low σ 8 , and the forecasted Ext-10k error bars.  \nThe selected alternative models are (a) a phenomenological model of non-perturbatively modified recombination (ModRec) (Lynch et al. 2024); (b) an atomic dark matter (ADM) model arising from a hidden dark sector charged under a U(1) gauge interaction (Hughes et al. 2024); and (c) a model featuring an additional dark fluid with an equation of state parameter that is a free function of redshift, hereafter 'high dimensional generalized dark matter' (HDGDM) (Meiers et al. 2023).  \nWe present one set of spectra from each of these models in Fig. 12 as a relative shift ∆ C X ℓ = C X ℓ -C X ℓ, ref compared to the reference ΛCDM mean Planck C X ℓ, ref , normalized as described in the figure caption. Here X corresponds to one of the primary lensed CMB TT/EE/TE or the lensing ϕϕ power spectra. Each of these spectra are drawn from posterior distributions conditioned on Planck and BAO data, although the HDGDM spectra also includes uncalibrated supernovae data in the fit. It should be noted that the specific BAO datasets used in these constraints differ: please see the appropriate references for details.  \nThe left-hand panels of Fig. 12 show these spectra in relation to existing CMB measurements from Planck , SPTpol (Bianchini et al. 2020), SPT-3G 2018 (Dutcher et al. 2021) and ACT DR6 (Qu et al. 2024). For angular scales that have been probed by multiple experiments, only the most constraining dataset has been plotted. This set of panels highlights that although these models make predictions that differ from ΛCDM, they are  \nFigure 12. Left panels: Current data from CMB experiments and some alternative models with either higher H 0 or lower σ 8 (see text). Except for lensing data, only the most constraining observations in a given multipole range have been plotted for clarity. Each spectrum has the best-fit ΛCDM spectrum given Planck subtracted as a reference and is normalized in an ℓ -dependent manner to reduce dynamic range, as described in the y -axis label. Note that for the ADM model, neither the model prediction nor the reference model includes non-linear corrections to the lensing power spectrum, as existing non-linear correction procedures do not account for the modified dark sector physics. All spectra are taken from chains constrained using Planck and BAO data, but the specific BAO datasets used differ between models. The HDGDM constraints also used uncalibrated supernovae data as well as the SH0ES H 0 measurement. Right panels: The same models as in the left panels, forecasted bandpower errors for Ext-10k, and the 68% confidence region given Planck and the assumption of ΛCDM. The bandpower errors have been binned with ∆ ℓ = 100. Note that the ΛCDM predictions conditioned on Planck are very tight, even in regions where Planck data on their own are completely unconstraining.  \n<!-- image -->  \ncurrently not ruled out by the data. To the right, we show the 68% confidence region for angular power spectrum predictions conditioned on Planck data and assuming ΛCDM, along with forecasted bandpower errors for Ext-10k in black.  \nThe Planck data, along with the assumption of ΛCDM, make tight predictions for the power spectrum residuals even at angular scales not probed by Planck . This is evident from the righthand panels of Figure 12, where the gray confidence interval extends to no more than about 2% for the TT , TE , and EE power spectra. This range represents the expected region that the power spectra observed by SPT-3G will fall in if measurements continue to be consistent with ΛCDM. As is evident in these panels, the alternative models presented here exhibit deviations which are larger than this expected range, and larger than the expected precision of SPT-3G. We now discuss each of these models in turn.  \nThe potential role of changes to recombination in alleviating the H 0 tension has been explored in Chiang & Slosar (2018); Hart & Chluba (2020a); Sekiguchi & Takahashi (2021); Lee et al. (2023); Cyr-Racine et al. (2022). Physical models which have a modified recombination epoch include models with varying fundamental constants (Ade et al. 2015; Hart & Chluba 2018; Hart & Chluba 2020b; Sekiguchi & Takahashi 2021), primordial magnetic fields (Jedamzik & Pogosian 2020; Galli et al. 2022; Rashkovetskyi et al. 2021), or energy injection from decaying dark matter (Galli et al. 2013; Slatyer &Wu2017; Finkbeiner et al. 2012). The ModRec model presented here is an extension of ΛCDM that allows for a purely phenomenological modification of the ionization fraction as a function of redshift through the epoch of hydrogen recombination, in order to study the role of these changes independently of the details of particular physical models. We note that the different modified ionization fractions predicted by the above physical models can be well approximated within the ModRec model using seven parameters characterizing deviations from the standard recombination scenario.  \nThis additional freedom allows for higher values of H 0 to provide good fits to the Planck data. The ModRec model accommodates a high value of H 0 by altering the recombination process, shifting last scattering to higher redshifts and increasing the width of the visibility function. This reduces the comoving size of the sound horizon at recombination, so that the shorter distance to last-scattering associated with a higher H 0 does not shift the well-determined angular size of the sound horizon. It also impacts the damping scale, an effect ameliorated by the change to the visibility function width.  \nIn Fig. 12 we show predictions within the ModRec model for a set of parameters which fit Planck and eBOSS BAO data at least as well as a ΛCDM best-fit, while having H 0 = 70 . 48 km/s/Mpc. To maintain the fit with Planck data, residuals with respect to the fiducial cosmology get shifted to higher multipoles, which, while unprobed by Planck , are expected to be measured with high SNR by SPT-3G. We see that SPT-3G has clear potential to rule out such a model, particularly through measurements of the TE spectrum. We note that without including BAO data, much higher values of H 0 are achievable within the ModRec model while staying consistent with Planck . In this case, if SPT-3G measurements remain consistent with ΛCDM, solutions to the H 0 tension which modify recombination will be seriously challenged by CMB data alone.  \nAnother approach to resolving the H 0 tension is to leave the recombination history unmodified and to instead introduce more freedom in the dark sector. The HDGDMmodel is another phenomenological model that uses the generalized dark matter (GDM) model of Hu (1998) and allows the equation of state parameter for the GDM to be an arbitrary function of redshift. Meiers et al. (2023) used such a model to approximate, and perturb around, a more physical model, the 'step-like' model of Aloni et al. (2022). These authors introduce a massive particle which becomes non-relativistic before decaying, in equilibrium, into a massless final state. The equation of state parameter w ( z ) for such a species begins at 1/3, drops as the particle becomes non-relativistic, and returns to 1/3 in the final massless state. This new component affects the expansion history in a non-trivial way, and alters the time-dependent gravitational driving of the plasma. They found that a transition at z ∼ 20000 could accommodate a higher H 0 . In Figure 12, we present the best-fit spectrum of the 'steplike restricted' model of Meiers et al. (2023) constrained using the ' D +' data combination in that work, which includes Planck BAO, uncalibrated supernovae data, as well a Gaussian likelihood centered at the SH0ES measurement H 0 = 73 . 04 ± 1 . 04 km/s/Mpc (Riess et al. 2022). This model has H 0 = 71 . 88 km/s/Mpc, illustrating how the HDGDM model can accommodate higher values of H 0 . However, the TT damping tail predictions from this model differ from the ΛCDM prediction using the Planck mean parameters at a more than 2% for ℓ > 3000, which should be distinguishable with SPT-3G.  \nIn addition to the H 0 tension, in recent years there has been increased attention to the potential discrepancy between different measurements of the amplitude of matter fluctuations on 8 h -1 Mpc scales, σ 8 . Inferences from Planck , that depend on assuming the ΛCDM  \nmodel, give a value of σ 8 = 0 . 811 ± 0 . 006, whereas the KiDS-1000 survey, measuring cosmic shear, gives a value of σ 8 = 0 . 76 +0 . 025 -0 . 020 (Planck Collaboration et al. 2020; Heymans et al. 2021). Although, DES sees less of a σ 8 tension compared to KiDS (Dark Energy Survey and Kilo-Degree Survey Collaboration 2023). In the ADM model considered here, this tension is addressed through the introduction of two dark fermions interacting with a dark photon, where the dark fermions make up some fraction f ADM of the dark matter. This results in dark sector dynamics analogous to the dynamics of the baryon-photon plasma. The pressure support from the dark photon at early times, prior to the formation of dark atoms, suppresses the growth of structure. In general, this leads to effects in the CMB power spectra observable by Planck but if the dark sector temperature is cool enough, the epoch of significant dark photon pressure support is early enough, and the Planck observables are relatively unaffected. More direct inferences of σ 8 are impacted by smaller scales than those to which Planck is sensitive; so, at the right dark photon temperature with a suppressed σ 8 , Planck observables are unaffected. In Fig. 12, we present the best-fit model identified in Hughes et al. (2024), coming from a region of parameter space consistent with Planck data, with σ 8 = 0 . 7869, a low dark photon temperature, and high f ADM . For this fit, both CMB data from Planck were used, as were BAO measurements (see Hughes et al. 2024 for details).  \nFrom this survey of alternatives to ΛCDM we conclude that measurements of temperature, polarization, and lensing power spectra from SPT-3G are capable of clarifying existing cosmological tensions.", '6. CONCLUSIONS': 'We presented forecasts of constraints on cosmological parameters that can be derived from the estimates of lensing, temperature and E-mode polarization power spectra to come from completed and ongoing SPT-3G sky surveys. We made forecasts for three combinations of the three SPT-3G surveys of varying sizes and depths, with particular emphasis on the combination of all three of these surveys, called Ext-10k, covering ∼ 10000 deg 2 (25% of the sky). We showed that SPT-3G will enable some powerful tests of the ΛCDM model.  \nWe showed that the constraining power of the full SPT-3G dataset will surpass that of Planck for the power spectra at modes ℓ ≳ 1800 in TT , ℓ ≳ 800 in TE , at ℓ ≳ 450 for EE and at L ≳ 30 in ϕϕ . We calculated the predictions of the ΛCDM model, conditioned on Planck data, and found that they are highly precise  \nin the power spectra on these scales, setting us up for a powerful test of ΛCDM.  \nWealso demonstrated that SPT-3G will facilitate testing of the ΛCDM model at the level of cosmological parameters. By propagating the uncertainties on power spectra to uncertainties on cosmological parameters, our analysis suggests that data from SPT-3G will be able to achieve similar, or in some cases, tighter constraints than those from Planck . Since these constraints incorporate higher weights from polarization, from small angular scales in temperature, and from lensing when compared to those from Planck , a consistency check between the two will be a powerful test of ΛCDM. The constraints on the current expansion rate ( H 0 ) and baryon density (Ω b h 2 ) improve by roughly a factor of 2 when compared with the constraints from Planck , while for the dark matter density (Ω c h 2 ), the improvement is a factor of 1.5. We find that the full SPT-3G dataset is expected to achieve an improvement in the ΛCDM FoM of 130; i.e., SPT-3G will reduce the allowed six-dimensional parameter volume by a factor of 130 compared to Planck .  \nThis reduction in volume is comparable to what is expected from the first two years of the SO-Baseline dataset. This similarity in constraining power, plus the substantial amount of overlap in sky coverage, provides an opportunity for significant map-level consistency tests. These can be used to detect, or limit, sources of systematic error, increasing the robustness of the parameter determinations from both SO and SPT3G.  \nWe also explored how well we can test ΛCDM by forecasting constraints on single- and double- parameter extensions with the combination of Planck +SPT-3G. We found that the combined data will result in considerable improvements over Planck alone. For example, we found that the constraints on single-parameter extensions improve by more than × 2 for θ ∈ [Ω k , N eff , w 0 , Y P ] and roughly × 1 . 5 for θ ∈ [ n run , ∑ m ν ] compared to Planck -only constraints. We also showed that including these extensions does not significantly degrade the constraints on the ΛCDM model parameters, except for H 0 , which is highly degenerate with some of the extensions like w 0 .  \nIn addition to showing the highly-precise predictions of the ΛCDM model on angular scales to be measured by SPT-3G, we showed that there are also alternative models that are consistent with current data, with significantly different predictions on these same angular scales. We considered three different models that could address existing cosmological tensions. While all of these models are consistent with Planck data, they exhibit deviations from the ΛCDM predictions-deviations that are larger than the forecasted SPT-3G error bars. This set of alter-  \ntive models illuminates the potential discovery space of SPT-3G observations.  \nThese SPT-3G tests of ΛCDM will begin soon and then improve in precision over the next five years. Later this year, we expect to release the first cosmological constraints from the first two years of the SPT-3G Main survey observations, and soon afterwards, from the first two years of Ext-4k survey. We forecast that this early dataset will offer constraining power on the ΛCDM model parameters that is comparable to that achieved by Planck , improving the ΛCDM FoM by a factor of 2 (SPT-3G Main) and 13 (Ext-4k) relative to Planck . A significant portion of the SPT-3G Wide survey has already been completed and we could expect Ext-10k results as early as sometime in 2025.', 'ACKNOWLEDGMENTS': "We thank the anonymous referee for their comments which has helped in shaping this manuscript better.  \nSR acknowledges support by the Illinois Survey Science Fellowship from the Center for AstroPhysical Surveys at the National Center for Supercomputing Applications. Work at Argonne National Lab is supported by UChicago Argonne LLC, Operator of Argonne National Laboratory (Argonne). Argonne, a U.S. Department of Energy Office of Science Laboratory, is operated under  \ncontract no. DE-AC02-06CH11357. The IAP group acknowledges funding from the European Research Council (ERC) under the European Union's Horizon 2020 research and innovation programme (grant agreement No 101001897).  \nThe South Pole Telescope program is supported by the National Science Foundation (NSF) through award OPP-1852617. Partial support is also provided by the Kavli Institute of Cosmological Physics at the University of Chicago.  \nThis work made use of the following computing resources: Illinois Campus Cluster, a computing resource that is operated by the Illinois Campus Cluster Program (ICCP) in conjunction with the National Center for Supercomputing Applications (NCSA) and which is supported by funds from the University of Illinois at Urbana-Champaign; the computational and storage services associated with the Hoffman2 Shared Cluster provided by UCLA Institute for Digital Research and Education's Research Technology Group; the computing resources provided on Crossover, a high-performance computing cluster operated by the Laboratory Computing Resource Center at Argonne National Laboratory; and the National Energy Research Scientific Computing Center, which is supported by the Office of Science of the U.S. Department of Energy under Contract No. DE-AC02-05CH11231.", 'A. DERIVATION OF THE MUSE COVARIANCE MATRIX': "In the context of parameter inferences from CMB data, the observed maps ( x ), do not depend just on cosmological parameters ( θ ) but also on unobserved latent variables ( z ) such as the unlensed CMB field and the lensing potential field. In such a scenario, the likelihood ( L ( x | θ )) involves marginalizing the joint likelihood L ( x, z | θ ) over the latent space.  \nL ( x | θ ) = ∫ d N z L ( x, z | θ ) = ∫ d N z L ( x | z, θ ) L ( z | θ ) , (A1)  \nwhere L ( z | θ ) represents the likelihood of the unlensed CMB field or the lensing potential z , given a set of cosmological parameters θ . For parameter inferences, we are often interested in the gradient of this marginal likelihood, a quantity called the marginal score s i .  \ns i ( θ, x ) ≡ d dθ i log L ( x | θ ) . (A2)  \nHowever, performing the integral over the maps of unlensed CMB and lensing potential which contain an order of million pixels is a computationally challenging task and analytic solutions only exist for the simplest cases. MUSE gets around this problem by approximating the marginal score. The key ingredient of the MUSE score is a quantity referred to as s MAP , which is the gradient of the joint likelihood evaluated at the maximum a posteriori (MAP) estimate of the unlensed CMB field and the lensing potential, ˆ z θ,x ≡ argmax z log L ( x, z | θ )  \ns MAP i ( θ, x ) ≡ d dθ i log L ( x, ˆ z θ,x | θ ) (A3)  \nand the MUSE score is defined as  \ns MUSE i ( θ, x ) ≡ s MAP i ( θ, x ) -〈 s MAP i ( θ, x ' ) 〉 x ' ∼P ( x ' | θ ) , (A4)  \nwhere the second term is calculated by averaging over simulations x ' drawn from a probability distribution generated at the specified value of θ . Once we have the MUSE score, we can define an estimator for the bandpowers θ as the root of the equation s MUSE i ( θ MUSE , x ) = 0  \nand the bandpower covariance can be constructed by two related matrices, J and H defined as:  \nJ ij = 〈 s MAP i ( θ ∗ , x ) s MAP j ( θ ∗ , x ) 〉 x ∼P ( x | θ ∗ ) -〈 s MAP i ( θ ∗ , x ) 〉 〈 s MAP j ( θ ∗ , x ) 〉 x ∼P ( x | θ ∗ ) H ij = d dθ j [ 〈 s MAP i ( θ ∗ , x ) 〉 x ∼P ( x | θ ) ] ∣ ∣ ∣ θ = θ ∗ . (A5)  \n∣ J is the covariance of the MAP gradients at the true value of bandpowers θ ∗ . We use simulations of the data in order to calculate this quantity. H can be thought of as the response matrix of the MAP gradients with respect to small changes in parameters that generate the simulated data. The bandpower covariance is constructed as follows:  \nΣ MUSE ij ≡ 〈 ∆ ˆ θ MUSE i ∆ ˆ θ MUSE j 〉 = ( H -1 JH -1 † ) ij . (A6)  \nA.1. Simulations  \nWe create simulations of the data, which consists of the lensed CMB along with noise as described in Table 1. Utilizing the best-fit Planck '15 spectrum, we produce a suite of 1000 flat-sky realizations of unlensed T, Q, and U CMB maps on a 512x512 pixel grid, spanning an area of 300 deg 2 , which are subsequently lensed by the lenseflow algorithm (Millea et al. 2019). We do not simulate the entire observation patch to limit the computational cost, thereby allowing us to run a large number of simulations. We do not expect this to have a significant impact as lensing is a local operation and the induced deflections are coherent only across a few degrees. Subsequently, we scaled this covariance matrix using an appropriate factor that corresponds to the ratio between our survey area and the simulated patch size. For the purposes of this forecast, we use a simplified version of the data model, which will be used in the actual analysis, whereby we ignore the effect of masking and filtering as is generally done with traditional fisher forecasts. We don't anticipate these to have a significant impact on the parameter constraints.  \nThe data model used in this work can be written as:  \nd = F · B · L ( ϕ ) f + n f ∼ C f ( C TT / TE / EE ℓ ) ϕ ∼ C ϕ ( C ϕϕ ℓ ) n ∼ C n ( N ℓ ) ,  \nwhere, d is the data; f is a Gaussian realization of the unlensed CMB map; ϕ is a Gaussian realization of the gravitational lensing potential; L ( ϕ ) is a linear lensing operator acting on f ; B is the instrumental beam function; F is a mid-pass Fourier filter that cuts off ℓ s below 300 and above 3500(4000) for TT,TE ( EE ), and N ℓ is the white+1 /f noise power spectrum as defined in Table 1.", 'A.2. Effective lensing reconstruction map noise': 'The effective lensing reconstruction noise is defined such that for an unbiased estimate of the lensing potential harmonic coefficients, ˆ ϕ LM , the per-mode variance of the residual to the true ϕ map is  \nN ϕϕ L = 〈 ( ˆ ϕ LM -ϕ true LM ) 2 〉 . (A7)  \nFor the quadratic estimate, for example, this is given by the sum of the N 0 , N 1 , N 3 / 2 , etc.. noise bias terms (Hu & Okamoto 2002; Madhavacheril et al. 2020). Note that if one has a biased estimate of the lensing potential, ¯ ϕ = Aϕ + n with some unknown normalization bias, A , and effective noise, n , the noise power can be written as  \nN ϕϕ L = C ϕϕ L ( 1 ρ L -1 ) , (A8)  \nwhere,  \nis the average correlation coefficient between the true and estimated lensing potential. This allows computing the effective noise without ever explicitly debiasing the lensing estimate to remove the bias, A .  \nTable 5. Forecasted 1 σ errors on single-parameter extensions to the ΛCDM for various SPT3G surveys in combination with Planck .  \nρ L = ⟨ ϕ true LM ¯ ϕ LM ⟩ √ ⟨ ϕ true LM ϕ true LM ⟩⟨ ¯ ϕ LM ¯ ϕ LM ⟩ (A9)  \nThe per-mode noise enters the error bars on the estimated lensing power spectrum as,  \n〈 (∆ C ϕϕ L ) 2 〉 = 2( C ϕϕ L + N ϕϕ L ) 2 ( f sky (2 L +1)) . (A10)  \nOur MUSE forecasts do not compute an unbiased lensing map estimate, ˆ ϕ LM , nor do they explicitly compute N ϕϕ L or use Eq.(A10). Instead, MUSE directly computes the total posterior bandpower covariance, including signal and noise contributions, using Eq.(A6). It is useful however to explore the effective noise levels which are reached. This can be conveniently performed by computing the averages in ρ L via Monte Carlo methods, where the biased lensing map estimate, ¯ ϕ , is the MAP estimate of ϕ , which is a byproduct of the MUSE inference procedure.', 'B. ADDITIONAL DETAILS ABOUT PARAMETER CONSTRAINTS': 'In Table 5 we present the constraints on the single-parameter extensions to ΛCDM from all the SPT-3G datasets. A comparison of the Ext-10k dataset with Planck is presented in Fig. 9. In Table 6, we present the degradation in cosmological parameters when we add one more parameter to the ΛCDM model. This is the same as the information presented in Fig. 10.  \nIn Fig. 13, we demonstrate the parameter constraints from Planck in black, and SPT-3G Main in green and  \nExt-10k in orange using the MUSE method. For SPT3G, we show the constraints that we get from delensed TT/EE/TE CMB-only as the dash-dotted and the ones with the inclusion of lensing ϕϕ as the solid curves. It is evident from the figure that the inclusion of lensing, breaks parameter degeneracies, which helps in improving the constraints on many cosmological parameters.', 'REFERENCES': "Abazajian, K., Addison, G. E., Adshead, P., et al. 2022, ApJ, 926, 54, doi: 10.3847/1538-4357/ac1596 Abbott, T. M. C., Aguena, M., Alarcon, A., et al. 2022, PhRvD, 105, 023520, doi: 10.1103/PhysRevD.105.023520 Ade, P. A. R., et al. 2015, Astron. Astrophys., 580, A22, doi: 10.1051/0004-6361/201424496  \nAde, P. A. R., Ahmed, Z., Amiri, M., et al. 2021, PhRvL, 127, 151301, doi: 10.1103/PhysRevLett.127.151301  \nAiola, S., Calabrese, E., Maurin, L., et al. 2020, Journal of Cosmology and Astroparticle Physics, 2020, 047-047,  \ndoi: 10.1088/1475-7516/2020/12/047  \n2  \nFigure 13. Marginalized 1-dimensional and 2-dimensional 68% contours for Hubble constant( H 0 ), scalar spectral index ( n s ), Cold dark matter density (Ω c h 2 ), and optical depth ( τ re ), using the MUSE method, for SPT-3G Main (dark green), Ext-10k (orange) and Planck (black). The dash-dotted curves show the constraints from delensed TT/TE/EE spectra whereas the solid curves also include lensing power spectrum information. The inclusion of the lensing spectrum significantly helps to break the degeneracies especially for the SPT-3G Main survey. The 1 σ errors are listed in the legend.  \n<!-- image -->  \nAlam, S., Aubert, M., Avila, S., et al. 2021, PhRvD, 103,  \n083533, doi: 10.1103/PhysRevD.103.083533  \nAloni, D., Berlin, A., Joseph, M., Schmaltz, M., & Weiner,  \nN. 2022, PhRvD, 105, 123516,  \ndoi: 10.1103/PhysRevD.105.123516  \nAnderson, A. J., Ade, P. A. R., Ahmed, Z., et al. 2018,  \nJournal of Low Temperature Physics,  \ndoi: 10.1007/s10909-018-2007-z  \nAsgari, M., Lin, C.-A., Joachimi, B., et al. 2021, A&A, 645, A104, doi: 10.1051/0004-6361/202039070  \nBalkenhol, L., Dutcher, D., Ade, P. A. R., et al. 2021, PhRvD, 104, 083509, doi: 10.1103/PhysRevD.104.083509 Bender, A. N., Ade, P. A. R., Ahmed, Z., et al. 2018, in Proc. SPIE, Vol. 10708, Proc. SPIE, 1070803, doi: 10.1117/12.2312426 Benson, B. A., Ade, P. A. R., Ahmed, Z., et al. 2014, in Proc. SPIE, Vol. 9153, Millimeter, Submillimeter, and Far-Infrared Detectors and Instrumentation for Astronomy VII, 91531P, doi: 10.1117/12.2057305 Bianchini, F., et al. 2020, Astrophys. J., 888, 119, doi: 10.3847/1538-4357/ab6082  \nH  \nTable 6. Degradation in constraints on standard cosmological parameters when the ΛCDM model is extended by single additional parameters. Note that H 0 is poorly constrained when w 0 is allowed to vary.  \nBleem, L. E., Stalder, B., de Haan, T., et al. 2015, ApJS, 216, 27, doi: 10.1088/0067-0049/216/2/27  \nBleem, L. E., Bocquet, S., Stalder, B., et al. 2020, ApJS, 247, 25, doi: 10.3847/1538-4365/ab6993  \nBleem, L. E., Klein, M., Abbot, T. M. C., et al. 2024, The Open Journal of Astrophysics, 7, 13,  \ndoi: 10.21105/astro.2311.07512  \nBocquet, S., Grandis, S., Bleem, L. E., et al. 2024, SPT Clusters with DES and HST Weak Lensing. II. Cosmological Constraints from the Abundance of Massive Halos. https://arxiv.org/abs/2401.02075  \n- Cardoso, J.-F., Le Jeune, M., Delabrouille, J., Betoule, M., & Patanchon, G. 2008, IEEE J. Sel. Top. Signal Process., 2, 735, doi: 10.1109/JSTSP.2008.2005346\n- Chiang, C.-T., & Slosar, A. 2018, arXiv e-prints, arXiv:1811.03624, doi: 10.48550/arXiv.1811.03624\n- Chichura, P. M., Foster, A., Patel, C., et al. 2022, The Astrophysical Journal, 936, 173,  \ndoi: 10.3847/1538-4357/ac89ec  \nChoi, S. K., Hasselfield, M., Ho, S.-P. P., et al. 2020, JCAP, 2020, 045, doi: 10.1088/1475-7516/2020/12/045  \n- Copi, C. J., Huterer, D., Schwarz, D. J., & Starkman, G. D. 2010, Advances in Astronomy, 2010, 1-17,  \ndoi: 10.1155/2010/847541  \nCyburt, R. H., Fields, B. D., Olive, K. A., & Yeh, T.-H.  \n2016, Reviews of Modern Physics, 88, 015004,  \ndoi: 10.1103/RevModPhys.88.015004  \nCyr-Racine, F.-Y., Ge, F., & Knox, L. 2022, Phys. Rev.  \nLett., 128, 201301, doi: 10.1103/PhysRevLett.128.201301  \nDark Energy Survey and Kilo-Degree Survey Collaboration. 2023, The Open Journal of Astrophysics, 6, 36,  \ndoi: 10.21105/astro.2305.17173  \nDatta, R., Aiola, S., Choi, S. K., et al. 2019, MNRAS, 486, 5239, doi: 10.1093/mnras/sty2934  \nDES Collaboration. 2024, The Dark Energy Survey: Cosmology Results With 1500 New High-redshift Type Ia Supernovae Using The Full 5-year Dataset. https://arxiv.org/abs/2401.02929  \nDutcher, D., Balkenhol, L., Ade, P. A. R., et al. 2021, PhRvD, 104, 022003, doi: 10.1103/PhysRevD.104.022003 Dutcher, D., et al. 2021, Phys. Rev. D, 104, 022003, doi: 10.1103/PhysRevD.104.022003  \n- Everett, W. B., Zhang, L., Crawford, T. M., et al. 2020, ApJ, 900, 55, doi: 10.3847/1538-4357/ab9df7  \nFinkbeiner, D. P., Galli, S., Lin, T., & Slatyer, T. R. 2012, Phys. Rev. D, 85, 043522,  \ndoi: 10.1103/PhysRevD.85.043522  \n- Galli, S., Pogosian, L., Jedamzik, K., & Balkenhol, L. 2022,\n- Phys. Rev. D, 105, 023513, doi: 10.1103/PhysRevD.105.023513\n- Galli, S., Slatyer, T. R., Valdes, M., & Iocco, F. 2013, Phys.\n- Rev. D, 88, 063502, doi: 10.1103/PhysRevD.88.063502\n- Galli, S., Wu, W. L. K., Benabed, K., et al. 2021, PhRvD, 104, 023518, doi: 10.1103/PhysRevD.104.023518\n- Galli, S., Benabed, K., Bouchet, F., et al. 2014, PhRvD, 90, 063504, doi: 10.1103/PhysRevD.90.063504\n- Givans, J. J., & Kamionkowski, M. 2023, arXiv e-prints, arXiv:2311.06196, doi: 10.48550/arXiv.2311.06196  \nGuns, S., Foster, A., Daley, C., et al. 2021, The Astrophysical Journal, 916, 98,  \ndoi: 10.3847/1538-4357/ac06a3  \nGupta, N., Reichardt, C. L., Ade, P. A. R., et al. 2019, MNRAS, 490, 5712, doi: 10.1093/mnras/stz2905  \nHart, L., & Chluba, J. 2018, Mon. Not. Roy. Astron. Soc., 474, 1850, doi: 10.1093/mnras/stx2783  \nHart, L., & Chluba, J. 2020a, MNRAS, 493, 3255,  \ndoi: 10.1093/mnras/staa412  \n-. 2020b, MNRAS, 495, 4210,  \ndoi: 10.1093/mnras/staa1426  \n- Schoneberg, N., Franco Abell'an, G., P'erez S'anchez, A., et al. 2022, Phys.Rept., 984, 1,\n- doi: 10.1016/j.physrep.2022.07.001\n- Sekiguchi, T., & Takahashi, T. 2021, PhRvD, 103, 083507,\n- doi: 10.1103/PhysRevD.103.083507\n- Sekiguchi, T., & Takahashi, T. 2021, Phys. Rev. D, 103, 083507, doi: 10.1103/PhysRevD.103.083507\n- Simard, G., Omori, Y., Aylor, K., et al. 2018, ApJ, 860, 137, doi: 10.3847/1538-4357/aac264\n- Simons Observatory Collaboration. 2019, JCAP, 2019, 056, doi: 10.1088/1475-7516/2019/02/056\n- Slatyer, T. R., & Wu, C.-L. 2017, Phys. Rev. D, 95, 023010, doi: 10.1103/PhysRevD.95.023010\n- Sobrin, J. A., Ade, P. A. R., Ahmed, Z., et al. 2018, in Proc. SPIE, Vol. 10708, Proc. SPIE, 107081H, doi: 10.1117/12.2314366  \n- Sobrin, J. A., Anderson, A. J., Bender, A. N., et al. 2022, ApJS, 258, 42, doi: 10.3847/1538-4365/ac374f\n- Tandoi, C., Guns, S., Foster, A., et al. 2024, Flaring Stars in a Non-targeted mm-wave Survey with SPT-3G. https://arxiv.org/abs/2401.13525\n- Tegmark, M., Taylor, A. N., & Heavens, A. F. 1997, ApJ, 480, 22\n- Thorne, B., Dunkley, J., Alonso, D., & Næss, S. 2017, MNRAS, 469, 2821, doi: 10.1093/mnras/stx949  \nTrendafilova, C. 2023, JCAP, 2023, 071,  \ndoi: 10.1088/1475-7516/2023/10/071  \n- Tristram, M., Banday, A. J., Douspis, M., et al. 2024, Astronomy &; Astrophysics, 682, A37, doi: 10.1051/0004-6361/202348015\n- Whitehorn, N., Natoli, T., Ade, P. A. R., et al. 2016, ApJ, 830, 143, doi: 10.3847/0004-637X/830/2/143\n- Zaldarriaga, M., Spergel, D. N., & Seljak, U. 1997, ApJ, 488, 1, doi: 10.1086/304692"}