Search

We present improvements to the pointing accuracy of the South Pole Telescope (SPT) using machine learning. The ability of the SPT to point accurately at the sky is limited by its structural imperfections, which are impacted by the extreme weather at the South Pole. Pointing accuracy is particularly important during SPT participation in observing campaigns with the Event Horizon Telescope (EHT), which requires stricter accuracy than typical observations with the SPT. We compile a training dataset of historical observations of astronomical sources made with the SPT-3G and EHT receivers on the SPT. We train two XGBoost models to learn a mapping from current weather conditions to two telescope drive control arguments -- one which corrects for errors in azimuth and the other for errors in elevation. Our trained models achieve root mean squared errors on withheld test data of $2.14''$ in cross-elevation and $3.57''$ in elevation, well below our goal of $5''$ along each axis. We deploy our models on the telescope control system and perform further in situ test observations during the EHT observing campaign in 2024 April. Our models result in significantly improved pointing accuracy: for sources within the range of input variables where the models are best trained, average combined pointing error improved 33%, from $15.9''$ to $10.6''$. These improvements, while significant, fall shy of our ultimate goal, but they serve as a proof of concept for the development of future models. Planned upgrades to the EHT receiver on the SPT will necessitate even stricter pointing accuracy which will be achievable with our methods., Comment: 23 pages, 9 figures, to be submitted to Journal of Astronomical Instrumentation (JAI)

Cosmic shear, galaxy clustering, and the abundance of massive halos each probe the large-scale structure of the universe in complementary ways. We present cosmological constraints from the joint analysis of the three probes, building on the latest analyses of the lensing-informed abundance of clusters identified by the South Pole Telescope (SPT) and of the auto- and cross-correlation of galaxy position and weak lensing measurements (3$\times$2pt) in the Dark Energy Survey (DES). We consider the cosmological correlation between the different tracers and we account for the systematic uncertainties that are shared between the large-scale lensing correlation functions and the small-scale lensing-based cluster mass calibration. Marginalized over the remaining $\Lambda$CDM parameters (including the sum of neutrino masses) and 52 astrophysical modeling parameters, we measure $\Omega_\mathrm{m}=0.300\pm0.017$ and $\sigma_8=0.797\pm0.026$. Compared to constraints from Planck primary CMB anisotropies, our constraints are only 15% wider with a probability to exceed of 0.22 ($1.2\sigma$) for the two-parameter difference. We further obtain $S_8\equiv\sigma_8(\Omega_\mathrm{m}/0.3)^{0.5}=0.796\pm0.013$ which is lower than the Planck measurement at the $1.6\sigma$ level. The combined SPT cluster, DES 3$\times$2pt, and Planck datasets mildly prefer a non-zero positive neutrino mass, with a 95% upper limit $\sum m_\nu<0.25~\mathrm{eV}$ on the sum of neutrino masses. Assuming a $w$CDM model, we constrain the dark energy equation of state parameter $w=-1.15^{+0.23}_{-0.17}$ and when combining with Planck primary CMB anisotropies, we recover $w=-1.20^{+0.15}_{-0.09}$, a $1.7\sigma$ difference with a cosmological constant. The precision of our results highlights the benefits of multiwavelength multiprobe cosmology., Comment: Submitted to Phys. Rev. D

From CMB polarization data alone we reconstruct the CMB lensing power spectrum, comparable in overall constraining power to previous temperature-based reconstructions, and an unlensed E-mode power spectrum. The observations, taken in 2019 and 2020 with the South Pole Telescope (SPT) and the SPT-3G camera, cover 1500 deg$^2$ at 95, 150, and 220 GHz with arcminute resolution and roughly 4.9$\mu$K-arcmin coadded noise in polarization. The power spectrum estimates, together with systematic parameter estimates and a joint covariance matrix, follow from a Bayesian analysis using the Marginal Unbiased Score Expansion (MUSE) method. The E-mode spectrum at $\ell>2000$ and lensing spectrum at $L>350$ are the most precise to date. Assuming the $\Lambda$CDM model, and using only these SPT data and priors on $\tau$ and absolute calibration from Planck, we find $H_0=66.81\pm0.81$ km/s/Mpc, comparable in precision to the Planck determination and in 5.4$\sigma$ tension with the most precise $H_0$ inference derived via the distance ladder. We also find $S_8=0.850\pm0.017$, providing further independent evidence of a slight tension with low-redshift structure probes. The $\Lambda$CDM model provides a good simultaneous fit to the combined Planck, ACT, and SPT data, and thus passes a powerful test. Combining these CMB datasets with BAO observations, we find that the effective number of neutrino species, spatial curvature, and primordial helium fraction are consistent with standard model values, and that the 95% confidence upper limit on the neutrino mass sum is 0.075 eV. The SPT data are consistent with the somewhat weak preference for excess lensing power seen in Planck and ACT data relative to predictions of the $\Lambda$CDM model. We also detect at greater than 3$\sigma$ the influence of non-linear evolution in the CMB lensing power spectrum and discuss it in the context of the $S_8$ tension.(abridged), Comment: 28 pages, 21 figures + appendices

The detection of satellite thermal emission at millimeter wavelengths is presented using data from the 3rd-Generation receiver on the South Pole Telescope (SPT-3G). This represents the first reported detection of thermal emission from artificial satellites at millimeter wavelengths. Satellite thermal emission is shown to be detectable at high signal-to-noise on timescales as short as a few tens of milliseconds. An algorithm for downloading orbital information and tracking known satellites given observer constraints and time-ordered observatory pointing is described. Consequences for cosmological surveys and short-duration transient searches are discussed, revealing that the integrated thermal emission from all large satellites does not contribute significantly to the SPT-3G survey intensity map. Measured satellite positions are found to be discrepant from their two-line element (TLE) derived ephemerides up to several arcminutes which may present a difficulty in cross-checking or masking satellites from short-duration transient searches.

Having a detailed knowledge of the intracluster medium (ICM) to infer the exact cluster physics such as the cluster dynamical state is crucial for cluster-based cosmological studies. This knowledge limits the accuracy and precision of mass estimation, a key parameter for such studies. In this paper, we conduct an in-depth analysis of cluster ACT-CL J0240.0+0116 using a multi-wavelength approach, with a primary focus on high angular resolution Sunyaev-Zeldovich (SZ) thermal component observations obtained under the NIKA2 Sunyaev-Zeldovich Large Programme (LPSZ). We create composite images using NIKA2, X-ray, and optical galaxy number density maps. The results reveal distinct signs of disturbance within the cluster with the distributions of gas and member galaxies that do not overlap. We also find suggestions of an inflow of matter onto the cluster from the southwestern direction. Ultimately, we classify the cluster as disturbed, using morphological indicators derived from its SZ, X-ray, and optical image. The cluster SZ signal is also contaminated by a strong central point source. We adopt different approaches to handling this contaminant and find the estimates of our pressure and hydrostatic mass profiles robust to the point source mitigation model. The cluster hydrostatic mass is estimated at $4.25^{+0.50}_{-0.45\, } \times 10^{14} \,\mathrm{M}_{\odot}$ for the case where the point source was masked. These values are consistent with the mass estimated using only X-ray data and with those from previous SZ studies of the Atacama cosmology telescope (ACT) survey, with improved precision on the mass estimate. Our findings strongly suggest that ACT-CL J0240.0+0116 is a disturbed cluster system, and the detailed observations and derived values serve as a compelling case study for the capabilities of the LPSZ in mapping the cluster ICM with high precision., Comment: 13 pages, 7 figures, submitted to A&A

The millimeter continuum emission from galaxies provides important information about cold dust, its distribution, heating, and role in their InterStellar Medium (ISM). This emission also carries an unknown portion of the free-free and synchrotron radiation. The IRAM 30m Guaranteed Time Large Project, Interpreting Millimeter Emission of Galaxies with IRAM and NIKA2 (IMEGIN) provides a unique opportunity to study the origin of the millimeter emission on angular resolutions of <18" in a sample of nearby galaxies. As a pilot study, we present millimeter observations of two IMEGIN galaxies, NGC 2146 (starburst) and NGC 2976 (peculiar dwarf) at 1.15 mm and 2 mm. Combined with the data taken with Spitzer, Herschel, Plank, WSRT, and the 100m Effelsberg telescopes, we model the infrared-to-radio Spectral Energy Distribution (SED) of these galaxies, both globally and at resolved scales, using a Bayesian approach to 1) dissect different components of the millimeter emission, 2) investigate the physical properties of dust, and 3) explore correlations between millimeter emission, gas, and Star Formation Rate (SFR). We find that cold dust is responsible for most of the 1.15 mm emission in both galaxies and at 2 mm in NGC 2976. The free-free emission emits more importantly in NGC 2146 at 2 mm. The cold dust emissivity index is flatter in the dwarf galaxy ($\beta = 1.3\pm 0.1$) compared to the starburst galaxy ($\beta = 1.7\pm 0.1$). Mapping the dust-to-gas ratio, we find that it changes between 0.004 and 0.01 with a mean of $0.006\pm0.001$ in the dwarf galaxy. In addition, no global balance holds between the formation and dissociation of H$_2$ in this galaxy. We find tight correlations between the millimeter emission and both the SFR and molecular gas mass in both galaxies., Comment: Accepted for publication in Astronomy & Astrophysics

In Sunyaev-Zeldovich (SZ) cluster cosmology, two tools are needed to be able to exploit data from large scale surveys in the millimeter-wave domain. An accurate description of the IntraCluster Medium (ICM) pressure profile is needed along with the scaling relation connecting the SZ brightness to the mass. With its high angular resolution and large field of view, The NIKA2 camera, operating at 150 and 260 GHz, is perfectly suited for precise cluster SZ mapping. The SZ Large Program (LPSZ) of the NIKA2 collaboration is dedicated to the observation of a sample of 38 SZ-selected clusters at intermediate to high redshift and observed both in SZ and X-ray. The current status is that all LPSZ clusters have been observed and the analysis toward the final results is ongoing. We present in detail how NIKA2-LPSZ will obtain a robust estimation of the SZ-Mass scaling relation and how it will be used to obtain cosmological constraints., Comment: 4 pages, 3 figures, contribution to the 2024 Cosmology session of the 58th Rencontres de Moriond

We introduce picasso, a model designed to predict thermodynamic properties of the intracluster medium based on the properties of halos in gravity-only simulations. The predictions result from the combination of an analytical gas model, mapping gas properties to the gravitational potential, and of a machine learning model to predict the model parameters for individual halos based on their scalar properties, such as mass and concentration. Once trained, the model can be applied to make predictions for arbitrary potential distributions, allowing its use with flexible inputs such as N-body particle distributions or radial profiles. We present the model, and train it using pairs of gravity-only and hydrodynamic simulations. We show that when trained to learn the mapping from gravity-only to non-radiative hydrodynamic simulations, picasso can make remarkably accurate and precise predictions of intracluster gas thermodynamics, with percent-level bias and $\sim 20 \%$ scatter for $r / R_{500c} \in [0.1,1]$. Training the model on hydrodynamic simulations including sub-resolution physics modeling yields robust predictions as well, albeit with the introduction of a radius-dependent bias and an increase in scatter. We further show that the model can be trained to make accurate predictions from very minimal halo information, down to mass and concentration, at the cost of modestly reduced precision. picasso is made publicly available as a Python package at https://github.com/fkeruzore/picasso, which includes trained models that can be used to make predictions easily and efficiently, in a fully auto-differentiable and hardware-accelerated framework, Comment: Accepted for publication in the Open Journal of Astrophysics; picasso is publicly available at https://github.com/fkeruzore/picasso

We present the detection and characterization of fluctuations in linearly polarized emission from the atmosphere above the South Pole. These measurements make use of Austral winter survey data from the SPT-3G receiver on the South Pole Telescope in three frequency bands centered at 95, 150, and 220 GHz. We use the cross-correlation between detectors to produce an unbiased estimate of the power in Stokes I, Q, and U parameters on large angular scales. Our results are consistent with the polarized signal being produced by the combination of Rayleigh scattering of thermal radiation from the ground and thermal emission from a population of horizontally aligned ice crystals with an anisotropic distribution described by Kolmogorov turbulence. The signal is most significant at large angular scales, high observing frequency, and low elevation angle. Polarized atmospheric emission has the potential to significantly impact observations on the large angular scales being targeted by searches for inflationary B-mode CMB polarization. We present the distribution of measured angular power spectrum amplitudes in Stokes Q and I for 4 years of winter observations, which can be used to simulate the impact of atmospheric polarization and intensity fluctuations at the South Pole on a specified experiment and observation strategy. For the SPT-3G data, downweighting the small fraction of significantly contaminated observations is an effective mitigation strategy. In addition, we present a strategy for further improving sensitivity on large angular scales where maps made in the 220 GHz band are used to measure and subtract the polarized atmosphere signal from the 150 GHz band maps. In observations with the SPT-3G instrument at the South Pole, the polarized atmospheric signal is a well-understood and sub-dominant contribution to the measured noise after implementing the mitigation strategies described here., Comment: 32 pages, 28 figures

We measure the stacked lensing signal in the direction of galaxy clusters in the Dark Energy Survey Year 3 (DES Y3) redMaPPer sample, using cosmic microwave background (CMB) temperature data from SPT-3G, the third-generation CMB camera on the South Pole Telescope (SPT). We estimate the lensing signal using temperature maps constructed from the initial 2 years of data from the SPT-3G 'Main' survey, covering 1500 deg$^2$ of the Southern sky. We then use this signal as a proxy for the mean cluster mass of the DES sample. In this work, we employ three versions of the redMaPPer catalogue: a Flux-Limited sample containing 8865 clusters, a Volume-Limited sample with 5391 clusters, and a Volume&Redshift-Limited sample with 4450 clusters. For the three samples, we find the mean cluster masses to be ${M}_{200{\rm{m}}}=1.66\pm0.13$ [stat.]$\pm0.03$ [sys.], $1.97\pm0.18$ [stat.]$\pm0.05$ [sys.], and $2.11\pm0.20$ [stat.]$\pm0.05$ [sys.]$\times{10}^{14}\ {\rm{M}}_{\odot }$, respectively. This is a factor of $\sim2$ improvement relative to the precision of measurements with previous generations of SPT surveys and the most constraining cluster mass measurements using CMB cluster lensing to date. Overall, we find no significant tensions between our results and masses given by redMaPPer mass-richness scaling relations of previous works, which were calibrated using CMB cluster lensing, optical weak lensing, and velocity dispersion measurements from various combinations of DES, SDSS and Planck data. We then divide our sample into 3 redshift and 3 richness bins, finding no significant tensions with optical weak-lensing calibrated masses in these bins. We forecast a $5.7\%$ constraint on the mean cluster mass of the DES Y3 sample with the complete SPT-3G surveys when using both temperature and polarization data and including an additional $\sim1400$ deg$^2$ of observations from the 'Extended' SPT-3G survey., Comment: 23 pages, 9 figures, accepted for publication in JCAP. Minor changes and corrections have been made relative to v1

We forecast constraints on cosmological parameters enabled by three surveys conducted with SPT-3G, the third-generation camera on the South Pole Telescope. The surveys cover separate regions of 1500, 2650, and 6000 ${\rm 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 $\mu{\rm 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 ~$\sqrt{\text{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 $\Lambda{\rm CDM}$ model via the consistency of cosmological parameters constrained independently from SPT-3G and Planck data, and consider the improvement in constraints on $\Lambda{\rm CDM}$ extension parameters from a joint analysis of SPT-3G and Planck data. The $\Lambda{\rm 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., Comment: 26 pages; 13 figures; Accepted for publication in ApJ; Minor edits have been made

We report results from an analysis aimed at detecting the trispectrum of the kinematic Sunyaev-Zel{'}dovich (kSZ) effect by combining data from the South Pole Telescope (SPT) and {\it Herschel}-SPIRE experiments over a 100 ${\rm deg}^{2}$ field. The SPT observations combine data from the previous and current surveys, namely SPTpol and SPT-3G, to achieve depths of 4.5, 3, and 16 $\mu {\rm K-arcmin}$ in bands centered at 95, 150, and 220 GHz. For SPIRE, we include data from the 600 and 857 GHz bands. We reconstruct the velocity-induced large-scale correlation of the small-scale kSZ signal with a quadratic estimator that uses two cosmic microwave background (CMB) temperature maps, constructed by optimally combining data from all the frequency bands. We reject the null hypothesis of a zero trispectrum at $10.3\sigma$ level. However, the measured trispectrum contains contributions from both the kSZ and other undesired components, such as CMB lensing and astrophysical foregrounds, with kSZ being sub-dominant. We use the \textsc{Agora} simulations to estimate the expected signal from CMB lensing and astrophysical foregrounds. After accounting for the contributions from CMB lensing and foreground signals, we do not detect an excess kSZ-only trispectrum and use this non-detection to set constraints on reionization. By applying a prior based on observations of the Gunn-Peterson trough, we obtain an upper limit on the duration of reionization of $\Delta z_{\rm re, 50} < 4.5$ (95\% C.L). We find these constraints are fairly robust to foregrounds assumptions. This trispectrum measurement is independent of, but consistent with, {\it Planck}'s optical depth measurement. This result is the first constraint on the epoch of reionization using the non-Gaussian nature of the kSZ signal., Comment: 15 pages, 5 figures (3 in main text and 2 in Appendix); Accepted for publication in PRL; Some texts have been moved to Appendix; Minor change in Fig. 2 to include nomalization; Data products and plotting scripts can be downloaded from https://github.com/sriniraghunathan/kSZ_4pt_SPT_SPIRE

We develop a new framework to constrain the source redshift. The method jointly accounts for the detection/non-detection of spectral lines and the prior information from the photometric redshift and total infrared luminosity from spectral energy distribution analysis. The method uses the estimated total infrared luminosity to predict the line fluxes at given redshifts and generates model spectra. The redshift-dependent spectral models are then compared with the observed spectra to find the redshift. Results. We apply the aforementioned joint redshift analysis method to four high-z dusty star-forming galaxy candidates selected from the NIKA2 observations of the HLSJ091828.6+514223 (HLS) field, and further observed by NOEMA with blind spectral scans. These sources only have SPIRE/Herschel photometry as ancillary data. They were selected because of very faint or no SPIRE counterparts, as to bias the sample towards the highest redshift candidates. The method finds the spectroscopic redshift of 4 in the 5 NOEMA-counterpart detected sources, with z>3. Based on these measurements, we derive the CO/[CI] lines and millimeter continuum fluxes from the NOEMA data and study their ISM and star-formation properties. We find cold dust temperatures in some of the HLS sources compared to the general population of sub-millimeter galaxies, which might be related to the bias introduced by the SPIRE-dropout selection. Our sources, but one, have short gas depletion time of a few hundred Myrs, which is typical among high-z sub-millimeter galaxies. The only exception shows a longer gas depletion time, up to a few Gyrs, comparable to that of main-sequence galaxies at the same redshift. Furthermore, we identify a possible over-density of dusty star-forming galaxies at z=5.2, traced by two sources in our sample, as well as the lensed galaxy HLSJ091828.6+514223. (abridged), Comment: A&A in press

We present a flare star catalog from four years of non-targeted millimeter-wave survey data from the South Pole Telescope (SPT). The data were taken with the SPT-3G camera and cover a 1500-square-degree region of the sky from $20^{h}40^{m}0^{s}$ to $3^{h}20^{m}0^{s}$ in right ascension and $-42^{\circ}$ to $-70^{\circ}$ in declination. This region was observed on a nearly daily cadence from 2019-2022 and chosen to avoid the plane of the galaxy. A short-duration transient search of this survey yields 111 flaring events from 66 stars, increasing the number of both flaring events and detected flare stars by an order of magnitude from the previous SPT-3G data release. We provide cross-matching to Gaia DR3, as well as matches to X-ray point sources found in the second ROSAT all-sky survey. We have detected flaring stars across the main sequence, from early-type A stars to M dwarfs, as well as a large population of evolved stars. These stars are mostly nearby, spanning 10 to 1000 parsecs in distance. Most of the flare spectral indices are constant or gently rising as a function of frequency at 95/150/220 GHz. The timescale of these events can range from minutes to hours, and the peak $\nu L_{\nu}$ luminosities range from $10^{27}$ to $10^{31}$ erg s$^{-1}$ in the SPT-3G frequency bands.

We present cosmological constraints from the abundance of galaxy clusters selected via the thermal Sunyaev-Zel'dovich (SZ) effect in South Pole Telescope (SPT) data with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). The cluster sample is constructed from the combined SPT-SZ, SPTpol ECS, and SPTpol 500d surveys, and comprises 1,005 confirmed clusters in the redshift range $0.25-1.78$ over a total sky area of 5,200 deg$^2$. We use DES Year 3 weak-lensing data for 688 clusters with redshifts $z<0.95$ and HST weak-lensing data for 39 clusters with $0.6

An accurate reconstruction of galaxy cluster masses is key to use this population of objects as a cosmological probe. In this work we present a study on the hydrostatic-to-lensing mass scaling relation for a sample of 53 clusters whose masses were reconstructed homogeneously in a redshift range between $z= 0.05$ and $1.07$. The $M_{500}$ mass for each cluster was indeed inferred from the mass profiles extracted from the X-ray and lensing data, without using a priori observable-mass scaling relations. We assessed the systematic dispersion of the masses estimated with our reference analyses with respect to other published mass estimates. Accounting for this systematic scatter does not change our main results, but enables the propagation of the uncertainties related to the mass reconstruction method or used dataset. Our analysis gives a hydrostatic-to-lensing mass bias of $(1-b) =0.739^{+0.075}_{-0.070}$ and no evidence of evolution with redshift. These results are robust against possible subsample differences.

We present a catalog of 689 galaxy cluster candidates detected at significance $\xi>4$ via their thermal Sunyaev-Zel'dovich (SZ) effect signature in 95 and 150 GHz data from the 500-square-degree SPTpol survey. We use optical and infrared data from the Dark Energy Camera and the Wide-field Infrared Survey Explorer (WISE) and \spitzer \ satellites, to confirm 544 of these candidates as clusters with $\sim94\%$ purity. The sample has an approximately redshift-independent mass threshold at redshift $z>0.25$ and spans $1.5 \times 10^{14} < M_{500c} < 9.1 \times 10^{14}$ $M_\odot/h_{70}$ \ and $0.031$. We use external radio data from the Sydney University Molonglo Sky Survey (SUMSS) to estimate contamination to the SZ signal from synchrotron sources. The contamination reduces the recovered $\xi$ by a median value of 0.032, or $\sim0.8\%$ of the $\xi=4$ threshold value, and $\sim7\%$ of candidates have a predicted contamination greater than $\Delta \xi = 1$. With the exception of a small number of systems $(<1\%)$, an analysis of clusters detected in single-frequency 95 and 150 GHz data shows no significant contamination of the SZ signal by emission from dusty or synchrotron sources. This cluster sample will be a key component in upcoming astrophysical and cosmological analyses of clusters. The SPTpol millimeter-wave maps and associated data products used to produce this sample are available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html, and the NASA LAMBDA website. An interactive sky server with the SPTpol maps and Dark Energy Survey data release 2 images is also available at NCSA https://skyviewer.ncsa.illinois.edu., Comment: Matches version accepted by OJA. 19 pages + references, 14 figures, cluster candidate table provided in Appendix. Data products available at https://pole.uchicago.edu/public/data/sptpol_500d_clusters/index.html and an interactive sky server at https://skyviewer.ncsa.illinois.edu

To understand the evolution of dust properties in molecular clouds in the course of the star formation process, we constrain the changes in the dust emissivity index from star-forming filaments to prestellar and protostellar cores to T Tauri stars. Using the NIKA2 continuum camera on the IRAM 30~m telescope, we observed the Taurus B211/B213 filament at 1.2\,mm and 2\,mm with unprecedented sensitivity and used the resulting maps to derive the dust emissivity index $\beta$. Our sample of 105 objects detected in the $\beta$ map of the B211/B213 filament indicates that, overall, $\beta$ decreases from filament and prestellar cores ($\beta \sim 2\pm0.5$) to protostellar cores ($\beta \sim 1.2 \pm 0.2$) to T-Tauri protoplanetary disk ($\beta < 1$). The averaged dust emissivity index $\beta$ across the B211/B213 filament exhibits a flat ($\beta \sim 2\pm0.3$) profile. This may imply that dust grain sizes are rather homogeneous in the filament, start to grow significantly in size only after the onset of the gravitational contraction/collapse of prestellar cores to protostars, reaching big sizes in T Tauri protoplanetary disks. This evolution from the parent filament to T-Tauri disks happens on a timescale of about 1-2~Myr., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

We present a Bayesian population modeling method to analyze the abundance of galaxy clusters identified by the South Pole Telescope (SPT) with a simultaneous mass calibration using weak gravitational lensing data from the Dark Energy Survey (DES) and the Hubble Space Telescope (HST). We discuss and validate the modeling choices with a particular focus on a robust, weak-lensing-based mass calibration using DES data. For the DES Year 3 data, we report a systematic uncertainty in weak-lensing mass calibration that increases from 1% at $z=0.25$ to 10% at $z=0.95$, to which we add 2% in quadrature to account for uncertainties in the impact of baryonic effects. We implement an analysis pipeline that joins the cluster abundance likelihood with a multi-observable likelihood for the Sunyaev-Zel'dovich effect, optical richness, and weak-lensing measurements for each individual cluster. We validate that our analysis pipeline can recover unbiased cosmological constraints by analyzing mocks that closely resemble the cluster sample extracted from the SPT-SZ, SPTpol ECS, and SPTpol 500d surveys and the DES Year 3 and HST-39 weak-lensing datasets. This work represents a crucial prerequisite for the subsequent cosmological analysis of the real dataset., Comment: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

High-resolution mapping of the hot gas in galaxy clusters is a key tool for cluster-based cosmological analyses. Taking advantage of the NIKA2 millimeter camera operated at the IRAM 30-m telescope, the NIKA2 SZ Large Program seeks to get a high-resolution follow-up of 38 galaxy clusters covering a wide mass range at intermediate to high redshift. The measured SZ fluxes will be essential to calibrate the SZ scaling relation and the galaxy clusters mean pressure profile, needed for the cosmological exploitation of SZ surveys. We present in this study a method to infer a mean pressure profile from cluster observations. We have designed a pipeline encompassing the map-making and the thermodynamical properties estimates from maps. We then combine all the individual fits, propagating the uncertainties on integrated quantities, such as $R_{500}$ or $P_{500}$, and the intrinsic scatter coming from the deviation to the standard self-similar model. We validate the proposed method on realistic LPSZ-like cluster simulations., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

We present the first preliminary results of the project \textit{ICED}, focusing on the face-on galaxy NGC4254. We use the millimetre maps observed with NIKA2 at IRAM-30m, as part of the IMEGIN Guaranteed Time Large Program, and of a wide collection of ancillary data (multi-wavelength photometry and gas phase spectral lines) that are publicly available. We derive the global and local properties of interstellar dust grains through infrared-to-radio spectral energy distribution fitting, using the hierarchical Bayesian code HerBIE, which includes the grain properties of the state-of-the-art dust model, THEMIS. Our method allows us to get the following dust parameters: dust mass, average interstellar radiation field, and fraction of small grains. Also, it is effective in retrieving the intrinsic correlations between dust parameters and interstellar medium properties. We find an evident anti-correlation between the interstellar radiation field and the fraction of small grains in the centre of NGC4254, meaning that, at strong radiation field intensities, very small amorphous carbon grains are efficiently destroyed by the ultra-violet photons coming from newly formed stars, through photo-desorption and sublimation. We observe a flattening of the anti-correlation at larger radial distances, which may be driven by the steep metallicity gradient measured in NGC4254., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

Three galaxy clusters selected from the XXL X-ray survey at high redshift and low mass ($z\sim1$ and $M_{500} \sim 1-2 \times 10^{14}$ M$_{\odot}$) were observed with NIKA2 to image their Sunyaev-Zel'dovich effect (SZ) signal. They all present an SZ morphology, together with the comparison with X-ray and optical data, that indicates dynamical activity related to merging events. Despite their disturbed intracluster medium, their high redshifts, and their low masses, the three clusters follow remarkably well the pressure profile and the SZ flux-mass relation expected from standard evolution. This suggests that the physics that drives cluster formation is already in place at $z \sim 1$ down to $M_{500} \sim 10^{14}$ M$_{\odot}$., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

The thermodynamical properties of the intracluster medium (ICM) are driven by scale-free gravitational collapse, but they also reflect the rich astrophysical processes at play in galaxy clusters. At low masses ($\sim 10^{14}$ M$_{\odot}$) and high redshift ($z \gtrsim 1$), these properties remain poorly constrained observationally, due to the difficulty in obtaining resolved and sensitive data. This paper aims at investigating the inner structure of the ICM as seen through the Sunyaev-Zel'dovich (SZ) effect in this regime of mass and redshift. Focus is set on the thermal pressure profile and the scaling relation between SZ flux and mass, namely the $Y_{\rm SZ} - M$ scaling relation. The three galaxy clusters XLSSC~072 ($z=1.002$), XLSSC~100 ($z=0.915$), and XLSSC~102 ($z=0.969$), with $M_{500} \sim 2 \times 10^{14}$ M$_{\odot}$, were selected from the XXL X-ray survey and observed with the NIKA2 millimeter camera to image their SZ signal. XMM-Newton X-ray data were used in complement to the NIKA2 data to derive masses based on the $Y_X - M$ relation and the hydrostatic equilibrium. The SZ images of the three clusters, along with the X-ray and optical data, indicate dynamical activity related to merging events. The pressure profile is consistent with that expected for morphologically disturbed systems, with a relatively flat core and a shallow outer slope. Despite significant disturbances in the ICM, the three high-redshift low-mass clusters follow remarkably well the $Y_{\rm SZ}-M$ relation expected from standard evolution. These results indicate that the dominant physics that drives cluster evolution is already in place by $z \sim 1$, at least for systems with masses above $M_{500} \sim 10^{14}$ M$_{\odot}$., Comment: 24 pages, published in A&A. Note that the title number had to be changed. This version matches the one from A&A

The NIKA2 camera operating at the IRAM 30 m telescope excels in high-angular resolution mapping of the thermal Sunyaev-Zeldovich effect towards galaxy clusters at intermediate and high-redshift. As part of the NIKA2 guaranteed time, the SZ Large Program (LPSZ) aims at tSZ-mapping a representative sample of SZ-selected galaxy clusters in the catalogues of the Planck satellite and of the Atacama Cosmology Telescope, and also observed in X-ray with XMM Newton or Chandra. Having completed observations in January 2023, we present tSZ maps of 38 clusters spanning the targeted mass ($3 < M_{500}/10^{14} M_{\odot} < 10$) and redshift ($0.5 < z < 0.9$) ranges. The first in depth studies of individual clusters highlight the potential of combining tSZ and X-ray observations at similar angular resolution for accurate mass measurements. These were milestones for the development of a standard data analysis pipeline to go from NIKA2 raw data to the thermodynamic properties of galaxy clusters for the upcoming LPSZ data release. Final products will include unprecedented measurements of the mean pressure profile and mass observable scaling relation using a distinctive SZ-selected sample, which will be key for ultimately improving the accuracy of cluster based cosmology., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

In the framework of the IMEGIN Large Program, we used the NIKA2 camera on the IRAM 30-m telescope to observe the edge-on galaxy NGC 891 at 1.15 mm and 2 mm and at a FWHM of 11.1" and 17.6", respectively. Multiwavelength data enriched with the new NIKA2 observations fitted by the HerBIE SED code (coupled with the THEMIS dust model) were used to constrain the physical properties of the ISM. Emission originating from the diffuse dust disk is detected at all wavelengths from mid-IR to mm, while mid-IR observations reveal warm dust emission from compact HII regions. Indications of mm excess emission have also been found in the outer parts of the galactic disk. Furthermore, our SED fitting analysis constrained the mass fraction of the small (< 15 Angstrom) dust grains. We found that small grains constitute 9.5% of the total dust mass in the galactic plane, but this fraction increases up to ~ 20% at large distances (|z| > 3 kpc) from the galactic plane., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

This study presents the first millimeter continuum mapping observations of two nearby galaxies, the starburst spiral galaxy NGC2146 and the dwarf galaxy NGC2976, at 1.15 mm and 2 mm using the NIKA2 camera on the IRAM 30m telescope, as part of the Guaranteed Time Large Project IMEGIN. These observations provide robust resolved information about the physical properties of dust in nearby galaxies by constraining their FIR-radio SED in the millimeter domain. After subtracting the contribution from the CO line emission, the SEDs are modeled spatially using a Bayesian approach. Maps of dust mass surface density, temperature, emissivity index, and thermal radio component of the galaxies are presented, allowing for a study of the relations between the dust properties and star formation activity (using observations at 24$\mu$m as a tracer). We report that dust temperature is correlated with star formation rate in both galaxies. The effect of star formation activity on dust temperature is stronger in NGC2976, an indication of the thinner interstellar medium of dwarf galaxies. Moreover, an anti-correlation trend is reported between the dust emissivity index and temperature in both galaxies., Comment: To appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

In cluster cosmology, cluster masses are the main parameter of interest. They are needed to constrain cosmological parameters through the cluster number count. As the mass is not an observable, a scaling relation is needed to link cluster masses to the integrated Compton parameters Y, i.e. the Sunyaev-Zeldovich observable (SZ). Planck cosmological results obtained with cluster number counts are based on a scaling relation measured with clusters at low redshift ($z$<0.5) observed in SZ and X-ray. In the SZ Large Program (LPSZ) of the NIKA2 collaboration, the scaling relation will be obtained with a sample of 38 clusters at intermediate to high redshift ($0.5

Dusty starless cores play an important role in regulating the initial phases of the formation of stars and planets. In their interiors, dust grains coagulate and ice mantles form, thereby changing the millimeter emissivities and hence the ability to cool. We mapped four regions with more than a dozen cores in the nearby Galactic filaments of Taurus and Perseus using the NIKA2 camera at the IRAM 30-meter telescope. Combining the 1mm to 2mm flux ratio maps with dust temperature maps from Herschel allowed to create maps of the dust emissivity index $\beta_{1,2}$ at resolutions of 2430 and 5600 a.u. in Taurus and Perseus, respectively. Here, we study the variation with total column densities and environment. $\beta_{1,2}$ values at the core centers ($A_V=12-19$mag) vary significantly between $\sim1.1$ and $2.3$. Several cores show a strong rise of $\beta_{1,2}$ from the outskirts at $\sim4$mag to the peaks of optical extinctions, consistent with the predictions of grain models and the gradual build-up of ice mantles on coagulated grains in the dense interiors of starless cores., Comment: to appear in Proc. of the mm Universe 2023 conference, Grenoble (France), June 2023, published by F. Mayet et al. (Eds), EPJ Web of conferences, EDP Sciences

We present an extension to a Sunyaev-Zel'dovich Effect (SZE) selected cluster catalog based on observations from the South Pole Telescope (SPT); this catalog extends to lower signal-to-noise than the previous SPT-SZ catalog and therefore includes lower mass clusters. Optically derived redshifts, centers, richnesses and morphological parameters together with catalog contamination and completeness statistics are extracted using the multi-component matched filter algorithm (MCMF) applied to the S/N>4 SPT-SZ candidate list and the Dark Energy Survey (DES) photometric galaxy catalog. The main catalog contains 811 sources above S/N=4, has 91% purity and is 95% complete with respect to the original SZE selection. It contains 50% more total clusters and twice as many clusters above z=0.8 in comparison to the original SPT-SZ sample. The MCMF algorithm allows us to define subsamples of the desired purity with traceable impact on catalog completeness. As an example, we provide two subsamples with S/N>4.25 and S/N>4.5 for which the sample contamination and cleaning-induced incompleteness are both as low as the expected Poisson noise for samples of their size. The subsample with S/N>4.5 has 98% purity and 96% completeness, and will be included in a combined SPT cluster and DES weak-lensing cosmological analysis. We measure the number of false detections in the SPT-SZ candidate list as function of S/N, finding that it follows that expected from assuming Gaussian noise, but with a lower amplitude compared to previous estimates from simulations., Comment: 16 pages, 17 figures, submitted to MNRAS

As the millimeter wavelength range remains a largely unexplored spectral region for galaxies, the IMEGIN large program aims to map the millimeter continuum emission of 22 nearby galaxies at 1.15 and 2 mm. Using the high-resolution maps produced by the NIKA2 camera, we explore the existence of very cold dust and take possible contamination by free-free and synchrotron emission into account. We study the IR-to-radio emission coming from different regions along the galactic plane and at large vertical distances. New observations of NGC 891, using the NIKA2 camera on the IRAM 30m telescope, along with a suite of observations at other wavelengths were used to perform a multiwavelength study of the spectral energy distribution in the interstellar medium in this galaxy. This analysis was performed globally and locally, using the advanced hierarchical Bayesian fitting code, HerBIE, coupled with the THEMIS dust model. Our dust modeling is able to reproduce the near-IR to millimeter emission of NGC 891, with the exception of an excess at a level of 25% obtained by the NIKA2 observations in the outermost parts of the disk. The radio continuum and thermal dust emission are distributed differently in the disk and galaxy halo. Different dusty environments are also revealed by a multiwavelength investigation of the emission features. Our detailed decomposition at millimeter and centimeter wavelengths shows that emission at 1 mm is purely originated by dust. Radio components become progressively important with increasing wavelengths. Finally, we find that emission arising from small dust grains accounts for ~ 9.5% of the total dust mass, reaching up to 20% at large galactic latitudes. Shock waves in the outflows that shatter the dust grains might explain this higher fraction of small grains in the halo., Comment: 15 pages, 10 figures, Accepted for publication in A&A

We present a measurement of gravitational lensing over 1500 deg$^2$ of the Southern sky using SPT-3G temperature data at 95 and 150 GHz taken in 2018. The lensing amplitude relative to a fiducial Planck 2018 $\Lambda$CDM cosmology is found to be $1.020\pm0.060$, excluding instrumental and astrophysical systematic uncertainties. We conduct extensive systematic and null tests to check the robustness of the lensing measurements, and report a minimum-variance combined lensing power spectrum over angular multipoles of $50

Synthetic datasets generated from large-volume gravity-only simulations are an important tool in the calibration of cosmological analyses. Their creation often requires accurate inference of baryonic observables from the dark matter field. We explore the effectiveness of a baryon pasting algorithm in providing precise estimations of three-dimensional gas thermodynamic properties based on gravity-only simulations. We use the Borg Cube, a pair of simulations originating from identical initial conditions, with one run evolved as a gravity-only simulation, and the other incorporating non-radiative hydrodynamics. Matching halos in both simulations enables comparisons of gas properties on an individual halo basis. This comparative analysis allows us to fit for the model parameters that yield the closest agreement between the gas properties in both runs. To capture the redshift evolution of these parameters, we perform the analysis at five distinct redshift steps, spanning from $z=0$ to $2$. We find that the investigated algorithm, utilizing information solely from the gravity-only simulation, achieves few-percent accuracy in reproducing the median intracluster gas pressure and density, albeit with a scatter of approximately 20%, for cluster-scale objects up to $z=2$. We measure the scaling relation between integrated Compton parameter and cluster mass ($Y_{500c} | M_{500c}$), and find that the imprecision of baryon pasting adds less than 5% to the intrinsic scatter measured in the hydrodynamic simulation. We provide best-fitting values and their redshift evolution, and discuss future investigations that will be undertaken to extend this work., Comment: 14 pages, 8 figures, 3 tables; for submission to the Open Journal of Astrophysics

Aims. Deep millimeter surveys are necessary to probe the dust-obscured galaxies at high redshift. We conducted a large observing program at 1.2 and 2 mm with the NIKA2 camera installed on the IRAM 30-meter telescope. This NIKA2 Cosmological Legacy Survey (N2CLS) covers two emblematic fields: GOODS-N and COSMOS. We introduce the N2CLS survey and present new 1.2 and 2 mm number count measurements based on the tiered N2CLS observations from October 2017 to May 2021. Methods. We develop an end-to-end simulation that combines an input sky model with the instrument noise and data reduction pipeline artifacts. This simulation is used to compute the sample purity, flux boosting, pipeline transfer function, completeness, and effective area of the survey. We used the 117 deg$^2$ SIDES simulations as the sky model, which include the galaxy clustering. Our formalism allows us to correct the source number counts to obtain galaxy number counts, the difference between the two being due to resolution effects caused by the blending of several galaxies inside the large beam of single-dish instruments. Results. The N2CLS-May2021 survey reaches an average 1-$\sigma$ noise level of 0.17 and 0.048 mJy on GOODS-N over 159 arcmin$^2$, and 0.46 and 0.14 mJy on COSMOS over 1010 arcmin$^2$, at 1.2 and 2 mm, respectively. For a purity threshold of 80%, we detect 120 and 67 sources in GOODS-N and 195 and 76 sources in COSMOS, at 1.2 and 2 mm, respectively. Our measurement connects the bright single-dish to the deep interferometric number counts. After correcting for resolution effects, our results reconcile the single-dish and interferometric number counts and are further accurately compared with model predictions., Comment: Accepted by A&A. 23 pages, 12 figures

We present panco2, an open-source Python library designed to extract galaxy cluster pressure profiles from maps of the thermal Sunyaev-Zeldovich effect. The extraction is based on forward modeling of the total observed signal, allowing to take into account usual features of millimeter observations, such as beam smearing, data processing filtering, and point source contamination. panco2 offers a large flexibility in the inputs that can be handled and in the analysis options, enabling refined analyses and studies of systematic effects. We detail the functionalities of the code, the algorithm used to infer pressure profile measurements, and the typical data products. We present examples of running sequences, and the validation on simulated inputs. The code is available on GitHub at https://github.com/fkeruzore/panco2, and comes with an extensive technical documentation to complement this paper at https://panco2.readthedocs.io., Comment: 21 pages, 17 figures, 2 tables; Accepted for publication in the Open Journal of Astrophysics

Distinctive large-scale structures have been identified in the spatial distribution of optical galaxies up to redshift z ~ 1. In the more distant universe, the relationship between the dust-obscured population of star-forming galaxies observed at millimetre wavelengths and the network of cosmic filaments of dark matter apparent in all cosmological hydrodynamical simulations is still under study. Using the NIKA2 dual-band millimetre camera, we mapped a field of ~ 90 arcminutes^2 in the direction of the star GJ526 simultaneously in its 1.15-mm and 2.0-mm continuum wavebands to investigate the nature of the quasi-alignment of five sources found ten years earlier with the MAMBO camera at 1.2 mm. We find that these sources are not clumps of a circumstellar debris disc around this star as initially hypothesized. Rather, they must be dust-obscured star-forming galaxies, or sub-millimetre galaxies (SMGs), in the distant background. The new NIKA2 map at 1.15 mm reveals a total of seven SMGs distributed in projection on the sky along a filament-like structure crossing the whole observed field. Furthermore, we show that the NIKA2 and supplemental Herschel photometric data are compatible with a model of the spectral energy distributions (SEDs) of these sources when a common redshift of 2.5 and typical values of the dust parameters for SMGs are adopted. Hence, we speculate that these SMGs might be located in a filament of the distant `cosmic web'. The length of this candidate cosmic filament crossing the whole map is at least 4 cMpc (comoving), and the separations between sources are between 0.25 cMpc and 1.25 cMpc at this redshift, in line with expectations from cosmological simulations. Nonetheless, further observations to determine the precise spectroscopic redshifts of these sources are required to definitively support this hypothesis of SMGs embedded in a cosmic filament of dark matter., Comment: Article accepted in Astronomy and Astrophysics ; 9 pages, 6 figures

The precise estimation of the mass of galaxy clusters is a major issue for cosmology. Large galaxy cluster surveys rely on scaling laws that relate cluster observables to their masses. From the high resolution observations of ~ 45 galaxy clusters with NIKA2 and XMM-Newton instruments, the NIKA2 SZ Large Program should provide an accurate scaling relation between the thermal Sunyaev-Zel'dovich effect and the hydrostatic mass. In this paper, we present an exhaustive analysis of the hydrostatic mass of the well known galaxy cluster CL J1226.9+3332, the highest-redshift cluster in the NIKA2 SZ Large Program at z = 0.89. We combine the NIKA2 observations with thermal Sunyaev-Zel'dovich data from NIKA, Bolocam and MUSTANG instruments and XMM-Newton X-ray observations and test the impact of the systematic effects on the mass reconstruction. We conclude that slight differences in the shape of the mass profile can be crucial when defining the integrated mass at R500, which demonstrates the importance of the modeling in the mass determination. We prove the robustness of our hydrostatic mass estimates by showing the agreement with all the results found in the literature. Another key information for cosmology is the bias of the masses estimated assuming hydrostatic equilibrium hypothesis. Based on the lensing convergence maps from the Cluster Lensing And Supernova survey with Hubble (CLASH) data, we obtain the lensing mass estimate for CL J1226.9+3332. From this we are able to measure the hydrostatic-to-lensing mass bias for this cluster, that spans from 1 - bHSE/lens ~ 0.7 to 1, presenting the impact of data-sets and mass reconstruction models on the bias.

PSZ2 G091.83+26.11 is a galaxy cluster with M500 = 7.43 x 10^14 Msun at z = 0.822 1. This object exhibits a complex morphology with a clear bimodality observed in X-rays. However, it was detected and analysed in the Planck sample as a single, spherical cluster following a universal profile 2. This model can lead to miscalculations of thermodynamical quantities, like the pressure profile. As future multiwavelength cluster experiments will detect more and more objects at high redshifts, it is crucial to quantify this systematic effect. In this work, we use high-resolution observations of the NIKA2 camera3,4,5,6 to integrate the morphological characteristics of the cluster in our modelling. This is achieved by fitting a two-halo model to the SZ image and then by reconstruction of the resulting projected pressure profile. We then compare these results with the spherical assumption., Comment: 4 pages, 3 figures Contribution to the 2022 Cosmology session of the 56th Rencontres de Moriond. arXiv admin note: substantial text overlap with arXiv:2111.05082