Search

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

Kinetic inductance detectors (KIDs) have been proven as reliable systems for astrophysical observations, especially in the millimetre range. Their compact size enables to optimally fill the focal plane, thus boosting sensitivity. The KISS (KIDs Interferometric Spectral Surveyor) instrument is a millimetre camera that consists of two KID arrays of 316 pixels each coupled to a Martin-Puplett interferometer (MPI). The addition of the MPI grants the KIDs camera the ability to provide spectral information in the 100 and 300 GHz range. In this paper we report the main properties of the KISS instrument and its observations. We also describe the calibration and data analysis procedures used. We present a complete model of the observed data including the sky signal and several identified systematics. We have developed a full photometric and spectroscopic data analysis pipeline that translates our observations into science-ready products. We show examples of the results of this pipeline on selected sources: Moon, Jupiter and Venus. We note the presence of a deficit of response with respect to expectations and laboratory measurements. The detectors noise level is consistent with values obtained during laboratory measurements, pointing to a sub-optimal coupling between the instrument and the telescope as the most probable origin for the problem. This deficit is large enough as to prevent the detection of galaxy clusters, which were KISS main scientific objective. Nevertheless, we have demonstrated the feasibility of this kind of instrument, in the prospect for other KID interferometers (such as the CONCERTO instrument). As this regard, we have developed key instrumental technologies such as optical conception, readout electronics and raw calibration procedures, as well as, adapted data analysis procedures., Comment: 23 pages, 15 figures. Accepted for publication in Publications of the Astronomical Society of the Pacific

Turbulent gas motions are expected to dominate the non-thermal energy budget of the intracluster medium (ICM). The measurement of pressure fluctuations from high angular resolution Sunyaev-Zel'dovich imaging opens a new avenue to study ICM turbulence, complementary to X-ray density fluctuation measures. We develop a methodological framework designed to optimally extract information on the ICM pressure fluctuation power spectrum statistics, and publicly release the associated software named PITSZI. We apply this tool to the NIKA data of the merging cluster MACSJ0717 to measure its pressure fluctuation power spectrum at high significance, and to investigate the implications for its nonthermal content. Depending on the choice of the radial pressure model and the details of the applied methodology, we measure an energy injection scale $L_{inj} \sim 800$ kpc. The power spectrum normalization corresponds to a characteristic amplitude reaching $A(k_{peak}) \sim 0.4$. These results are are obtained assuming that MACSJ0717 can be described as pressure fluctuations on top of a single (smooth) halo, and are dominated by systematics due to the choice of the radial pressure model. Using simulations, we estimate that fitting a radial model to the data can suppress the observed fluctuations by up to 50\%, while a poorly representative radial model can induce spurious fluctuations, which we also quantify. Assuming standard scaling relations between the pressure fluctuations and turbulence, we find that MACSJ0717 presents a turbulent velocity dispersion $\sigma_v \sim 1200$ km/s, a kinetic to kinetic plus thermal pressure fraction $P_{k} / P_{k+th} \sim 20\%$, and we estimate the hydrostatic mass bias to $b_{HSE} \sim 0.3-0.4$. Our results are in excellent agreement with alternative measurements from X-ray surface brightness fluctuations, and in agreement with the fluctuations being adiabatic in nature., Comment: 22 pages, submitted to A&A

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

The cosmic microwave background (CMB), a remnant of the Big Bang, provides unparalleled insights into the primordial universe, its energy content, and the origin of cosmic structures. The success of forthcoming terrestrial and space experiments hinges on meticulously calibrated data. Specifically, the ability to achieve an absolute calibration of the polarization angles with a precision of < 0.1 deg is crucial to identify the signatures of primordial gravitational waves and cosmic birefringence within the CMB polarization. We introduce the COSMOCal project, designed to deploy a polarized source in space for calibrating microwave frequency observations. The project aims to integrate microwave polarization observations from small and large telescopes, ground-based and in space, into a unified scale, enhancing the effectiveness of each observatory and allowing robust combination of data. To demonstrate the feasibility and confirm the observational approach of our project, we developed a prototype instrument that operates in the atmospheric window centered at 260 GHz, specifically tailored for use with the NIKA2 camera at the IRAM 30 m telescope. We present the instrument components and their laboratory characterization. The results of tests performed with the fully assembled prototype using a KIDs-based instrument, similar concept of NIKA2, are also reported. This study paves the way for an observing campaign using the IRAM 30m telescope and contributes to the development of a space-based instrument., Comment: 13 pages, 13 figures. Published

The abundance of galaxy clusters with mass and redshift is a well-known cosmological probe. The cluster mass is a key parameter for studies that aim to constrain cosmological parameters using galaxy clusters, making it critical to understand and properly account for the errors in its estimates. Subsequently, it becomes important to correctly calibrate scaling relations between observables like the integrated Compton parameter and the mass of the cluster. The NIKA2 Sunyaev-Zeldovich Large program (LPSZ) enables one to map the intracluster medium profiles in the mm-wavelength band with great details (resolution of $11 \ \mathrm{\&}\ 17^{\prime \prime}$ at $1.2 \ \mathrm{\&}\ 2 $ mm, respectively) and hence, to estimate the cluster hydrostatic mass more precisely than previous SZ observations. However, there are certain systematic effects which can only be accounted for with the use of simulations. For this purpose, we employ THE THREE HUNDRED simulations which have been modelled with a range of physics modules to simulate galaxy clusters. The so-called twin samples are constructed by picking synthetic clusters of galaxies with properties close to the observational targets of the LPSZ. In particular, we use the Compton parameter maps and projected total mass maps of these twin samples along 29 different lines of sight. We investigate the scatter that projection induces on the total masses. Eventually, we consider the statistical values along different lines of sight to construct a kind of 3D scaling law between the integrated Compton parameter, total mass, and overdensity of the galaxy clusters to determine the overdensity that is least impacted by the projection effect., 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 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

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.

The Cosmic Microwave Background (CMB) radiation offers a unique window into the early Universe, facilitating precise examinations of fundamental cosmological theories. However, the quest for detecting B-modes in the CMB, predicted by theoretical models of inflation, faces substantial challenges in terms of calibration and foreground modeling. The COSMOCal (COsmic Survey of Millimeter wavelengths Objects for CMB experiments Calibration) project aims at enhancing the accuracy of the absolute calibration of the polarization angle $\psi$ of current and future CMB experiments. The concept includes the build of a very well known artificial source emitting in the frequency range [20-350] GHz that would act as an absolute calibrator for several polarization facilities on Earth. A feasibility study to place the artificial source in geostationary orbit, in the far field for all the telescopes on Earth, is ongoing. In the meanwhile ongoing hardware work is dedicated to build a prototype to test the technology, the precision and the stability of the polarization recovering in the 1 mm band (220-300 GHz). High-resolution experiments as the NIKA2 camera at the IRAM 30m telescope will be deployed for such use. Once carefully calibrated ($\Delta\psi$ < 0.1 degrees) it will be used to observe astrophysical sources such as the Crab nebula, which is the best candidate in the sky for the absolute calibration of CMB experiments., 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

class_sz is a versatile and robust code in C and Python that can compute theoretical predictions for a wide range of observables relevant to cross-survey science in the Stage IV era. The code is public at https://github.com/CLASS-SZ/class_sz along with a series of tutorial notebooks (https://github.com/CLASS-SZ/notebooks). It will be presented in full detail in paper II. Here we give a brief overview of key features and usage., 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

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

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

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

Cluster number counts will be a key cosmological probe in the next decade thanks to the Euclid satellite mission. For this purpose, cluster detection algorithm performance, which are sensitive to the spatial distribution of the cluster galaxy members and their luminosity function, need to be accurately characterized. Using The Three Hundred hydrodynamical and dark matter only simulations we study a complete sample of massive clusters beyond 7 (5) $\times$ 10$^{14}$ M$_{\odot}$ at redshift 0 (1) on a $(1.48 \ \mathrm{Gpc})^3$ volume. We find that the mass resolution of the current hydrodynamical simulations (1.5 $\times$ 10$^9$ M$_{\odot}$) is not enough to characterize the luminosity function of the sample in the perspective of Euclid data. Nevertheless, these simulations are still useful to characterize the spatial distribution of the cluster substructures assuming a common relative mass threshold for the different flavours and resolutions. By comparing with the dark matter only version of these simulations, we demonstrate that baryonic physics preserves significantly low mass subhalos (galaxies) as have also been observed in previous studies with less statistics. Furthermore, by comparing the hydro simulations with higher resolution dark matter only simulations of the same objects and taking the same limit in subhalo mass we find significantly more cuspy galaxy density profiles towards the center of the clusters, where the low mass substructures would tend to concentrate. We conclude that using dark matter only simulation may lead to some biases on the spatial distribution and density of galaxy cluster members. Based on the preliminary analysis of few high resolution hydro simulations we conclude that a mass resolution of 1.8 $\times$ 10$^8$ h$^{-1}$ M$_{\odot}$ will be needed for The Three Hundred simulations to approach the expected magnitude limits for the Euclid survey.

CONCERTO is the first experiment to perform a [CII] line intensity mapping survey to target $z>5.2$. Measuring the [CII] power spectrum allows us to study the role of dusty star-forming galaxies in the star formation history during the Reionization and post-Reionization. The main obstacle to this measurement is the contamination by bright foregrounds. We evaluate our ability to retrieve the [CII] signal in mock observations using the Simulated Infrared Dusty Extragalactic Sky. We compared two methods for dealing with the dust continuum emission from galaxies: the standard PCA and the arPLS method. For line interlopers, the strategy relies on masking low-redshift galaxies using external catalogues. As we do not have observations of CO or classical CO proxies ,we relied on the COSMOS stellar mass catalogue. To measure the power spectrum of masked data, we adapted the P of K EstimatoR and discuss its use on LIM data. The arPLS method achieves a reduction of the continuum background to a sub-dominant level of the [CII] at z=7 by a factor of>70. When using PCA, this factor is only 0.7. The masking lowers the power amplitude of line contamination down to $2 \times 10^2 Jy^2/sr$ This residual level is dominated by faint undetected sources. For our [CII] model, this results in a detection at z = 5.2 with a power ratio [CII]/(residual interlopers) = $62 \pm 32$ for a 22 % area survey loss. However, at z = 7, [C II ] / (residual interlopers)$=2.0 \pm 1.4$. Thanks to the large area covered by SIDES-Uchuu, we show that the power amplitude of line residuals varies by 12-15% for z=5.2-7. We present an end-to-end simulation of the extragalactic foreground removal that we ran to detect the [CII] at high redshift via its power spectrum. We show that dust continuum emission are not a limiting foreground for [CII] LIM. Residual CO and [CI] limits our ability to measure the [CII] power spectrum at z>7., Comment: 15 pages, 12 figures, to be published in Astronomy & Astrophysics

The determination of the mass of galaxy clusters from observations is subject to systematic uncertainties. Beyond the errors due to instrumental and observational systematic effects, in this work we investigate the bias introduced by modelling assumptions. In particular, we consider the reconstruction of the mass of galaxy clusters from convergence maps employing spherical mass density models. We made use of The Three Hundred simulations, selecting clusters in the same redshift and mass range as the NIKA2 Sunyaev-Zel'dovich Large Programme sample: $3 \leq M_{500}/ 10^{14} \mathrm{M}_{\odot} \leq 10$ and $0.5 \leq z \leq 0.9$. We studied different modelling and intrinsic uncertainties that should be accounted for when using the single cluster mass estimates for scaling relations. We confirm that the orientation of clusters and the radial ranges considered for the fit have an important impact on the mass bias. The effect of the projection adds uncertainties to the order of $10\%$ to $16\%$ to the mass estimates. We also find that the scatter from cluster to cluster in the mass bias when using spherical mass models is less than $9\%$ of the true mass of the clusters.

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

In the submm regime, spectral line scans and line intensity mapping (LIM) are new promising probes for the cold gas content and star formation rate of galaxies across cosmic time. However, both of these two measurements suffer from field-to-field variance. We study the effect of field-to-field variance on the predicted CO and [CII] power spectra from future LIM experiments such as CONCERTO, as well as on the line luminosity functions (LFs) and the cosmic molecular gas mass density that are currently derived from spectral line scans. We combined a 117 $\rm deg^2$ dark matter lightcone from the Uchuu cosmological simulation with the simulated infrared dusty extragalactic sky (SIDES) approach. We find that in order to constrain the CO LF with an uncertainty below 20%, we need survey sizes of at least 0.1 $\rm deg^2$. Furthermore, accounting for the field-to-field variance using only the Poisson variance can underestimate the total variance by up to 80%. The lower the luminosity is and the larger the survey size is, the higher the level of underestimate. At $z$<3, the impact of field-to-field variance on the cosmic molecular gas density can be as high as 40% for the 4.6 arcmin$^2$ field, but drops below 10% for areas larger than 0.2 deg$^2$. However, at $z>3$ the variance decreases more slowly with survey size and for example drops below 10% for 1 deg$^2$ fields. Finally, we find that the CO and [CII] LIM power spectra can vary by up to 50% in $\rm 1 deg^2$ fields. This limits the accuracy of the constraints provided by the first 1 deg$^2$ surveys. The level of the shot noise power is always dominated by the sources that are just below the detection thresholds. We provide an analytical formula to estimate the field-to-field variance of current or future LIM experiments. The code and the full SIDES-Uchuu products (catalogs, cubes, and maps) are publicly available.

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.

The CONCERTO spectral-imaging instrument was installed at the Atacama Pathfinder EXperiment (APEX) 12-meter telescope in April 2021. It has been designed to look at radiation emitted by ionised carbon atoms, [CII], and use the "intensity Mapping" technique to set the first constraints on the power spectrum of dusty star-forming galaxies. The instrument features two arrays of 2152 pixels constituted of Lumped Element Kinectic Inductance Detectors (LEKID) operated at cryogenic temperatures, cold optics and a fast Fourier Transform Spectrometer (FTS). To readout and operate the instrument, a newly designed electronic system hosted in five microTCA crates and composed of twelve readout boards and two control boards was designed and commissioned. The architecture and the performances are presented in this paper.

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning., Comment: accepted for publication in MNRAS, 21 pages, 22 figures, 6 tables

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