Showing 1–50 of 589 results for author: Jarvis, M

NOEMA formIng Cluster survEy (NICE): Characterizing eight massive galaxy groups at $1.5 < z < 4$ in the COSMOS field

Authors: Nikolaj B. Sillassen, Shuowen Jin, Georgios E. Magdis, Emanuele Daddi, Tao Wang, Shiying Lu, Hanwen Sun, Vinod Arumugam, Daizhong Liu, Malte Brinch, Chiara D'Eugenio, Raphael Gobat, Carlos Gómez-Guijarro, Michael Rich, Eva Schinnerer, Veronica Strazzullo, Qinghua Tan, Francesco Valentino, Yijun Wang, Mengyuan Xiao, Luwenjia Zhou, David Blánquez-Sesé, Zheng Cai, Yanmei Chen, Laure Ciesla , et al. (19 additional authors not shown)

Abstract: The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are c… ▽ More The NOEMA formIng Cluster survEy (NICE) is a large program targeting 69 massive galaxy group candidates at $z>2$ in six deep fields. We report spectroscopic confirmation of eight groups at $1.65\leq z\leq3.61$ in COSMOS. Homogeneously selected as significant overdensities of red IRAC sources with red Herschel colors, four groups are confirmed by CO and [CI] with NOEMA 3mm observations, three are confirmed with ALMA, and one is confirmed by H$α$ from Subaru/FMOS. We constructed the integrated FIR SEDs for the eight groups, obtaining total IR SFR $=260-1300~{\rm M_\odot}$~yr$^{-1}$. We adopted six methods to estimate the dark matter masses, including stellar mass to halo mass relations, overdensity with galaxy bias, and NFW profile fitting to radial stellar mass density. We found the radial stellar mass density are consistent with a NFW profile, supporting that they are collapsed structures hosted by a single dark matter halo. The best halo mass estimates are $\log(M_{\rm h}/{\rm M_\odot})=12.8-13.7$ with uncertainty of 0.3 dex. From halo mass estimates, we derive baryonic accretion rate ${\rm BAR}=(1-8)\times10^{3}\,{\rm M_{\odot}/yr}$ for this sample. We find a quasi-linear correlation between the integrated SFR/BAR and the theoretical halo mass limit for cold streams, $M_{\rm stream}/M_{\rm h}$, with ${\rm SFR/BAR}=10^{-0.46\pm0.22}\left({M_{\rm stream}/M_{\rm h}}\right)^{0.71\pm0.16}$ with a scatter of $0.40\,{\rm dex}$. Further, we compare halo masses and stellar masses with simulations, and find all structures are consistent with being progenitors of $M_{\rm h}(z=0)>10^{14}\,{\rm M_{\odot}}$ galaxy clusters, and the most massive central galaxies have stellar masses consistent with brightest cluster galaxies (BCGs) progenitors in the TNG300 simulation. The results strongly suggest these structures are forming massive galaxy clusters via baryonic and dark matter accretion. △ Less

Submitted 3 July, 2024; originally announced July 2024.

Comments: 44 pages (27pp appendix), 32 figures, 18 tables, accepted for publication in A&A

Fast modeling of the shear three-point correlation function

Authors: Sunao Sugiyama, Rafael C. H. Gomes, Mike Jarvis

Abstract: The three-point correlation function (3PCF) of a weak lensing shear field contains information that is complementary to that in the two-point correlation function (2PCF), which can help improve the cosmological parameters and calibrate astrophysical and observational systematics parameters. However, the application of the 3PCF to observed data has been limited due to the computational challenges o… ▽ More The three-point correlation function (3PCF) of a weak lensing shear field contains information that is complementary to that in the two-point correlation function (2PCF), which can help improve the cosmological parameters and calibrate astrophysical and observational systematics parameters. However, the application of the 3PCF to observed data has been limited due to the computational challenges of calculating theoretical predictions of the 3PCF from a bispectrum model. In this paper, we present a new method to compute the shear 3PCF efficiently and accurately. We employ the multipole expansion of the bispectrum to compute the shear 3PCF, and show that the method is substantially more efficient than direct numerical integration. We found that the multipole-based method can compute the shear 3PCF with 5% accuracy in 10 (40) seconds for the single (four) source redshift bin setup. The multipole-based method can be also used to compute the third-order aperture mass statistics quickly and accurately, accounting for the bin-averaging effect on the shear 3PCF. Our method provides a fast and robust tool for probing the underlying cosmological model with third-order statistics of weak lensing shear. △ Less

Submitted 1 July, 2024; originally announced July 2024.

Comments: 20 pages, 9 figures, 2 tables, code available at https://github.com/git-sunao/fastnc

Dark Energy Survey Year 3 Results: Cosmology from galaxy clustering and galaxy-galaxy lensing in harmonic space

Authors: L. Faga, F. Andrade-Oliveira, H. Camacho, R. Rosenfeld, M. Lima, C. Doux, X. Fang, J. Prat, A. Porredon, M. Aguena, A. Alarcon, S. Allam, O. Alves, A. Amon, S. Avila, D. Bacon, K. Bechtol, M. R. Becker, G. M. Bernstein, S. Bocquet, D. Brooks, E. Buckley-Geer, A. Campos, A. Carnero Rosell, M. Carrasco Kind , et al. (78 additional authors not shown)

Abstract: We present the joint tomographic analysis of galaxy-galaxy lensing and galaxy clustering in harmonic space, using galaxy catalogues from the first three years of observations by the Dark Energy Survey (DES Y3). We utilise the redMaGiC and MagLim catalogues as lens galaxies and the METACALIBRATION catalogue as source galaxies. The measurements of angular power spectra are performed using the pseudo… ▽ More We present the joint tomographic analysis of galaxy-galaxy lensing and galaxy clustering in harmonic space, using galaxy catalogues from the first three years of observations by the Dark Energy Survey (DES Y3). We utilise the redMaGiC and MagLim catalogues as lens galaxies and the METACALIBRATION catalogue as source galaxies. The measurements of angular power spectra are performed using the pseudo-$C_\ell$ method, and our theoretical modelling follows the fiducial analyses performed by DES Y3 in configuration space, accounting for galaxy bias, intrinsic alignments, magnification bias, shear magnification bias and photometric redshift uncertainties. We explore different approaches for scale cuts based on non-linear galaxy bias and baryonic effects contamination. Our fiducial covariance matrix is computed analytically, accounting for mask geometry in the Gaussian term, and including non-Gaussian contributions and super-sample covariance terms. To validate our harmonic space pipelines and covariance matrix, we used a suite of 1800 log-normal simulations. We also perform a series of stress tests to gauge the robustness of our harmonic space analysis. In the $Λ$CDM model, the clustering amplitude $S_8 =σ_8(Ω_m/0.3)^{0.5}$ is constrained to $S_8 = 0.704\pm 0.029$ and $S_8 = 0.753\pm 0.024$ ($68\%$ C.L.) for the redMaGiC and MagLim catalogues, respectively. For the $w$CDM, the dark energy equation of state is constrained to $w = -1.28 \pm 0.29$ and $w = -1.26^{+0.34}_{-0.27}$, for redMaGiC and MagLim catalogues, respectively. These results are compatible with the corresponding DES Y3 results in configuration space and pave the way for harmonic space analyses using the DES Y6 data. △ Less

Submitted 18 June, 2024; originally announced June 2024.

Comments: To be submitted to MNRAS

Report number: FERMILAB-PUB-24-0289-PPD

Probabilistic and progressive deblended far-infrared and sub-millimetre point source catalogues I. Methodology and first application in the COSMOS field

Authors: Lingyu Wang, Antonio La Marca, Fangyou Gao, William J. Pearson, Berta Margalef-Bentabol, Matthieu Béthermin, Longji Bing, James Donnellan, Peter D. Hurley, Seb J. Oliver, Catherine L. Hale, Matt J. Jarvis, Lucia Marchetti, Mattia Vaccari, Imogen H. Whittam

Abstract: Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the fi… ▽ More Single-dish far-infrared (far-IR) and sub-millimetre (sub-mm) point source catalogues and their connections with catalogues at other wavelengths are of paramount importance. However, due to the large mismatch in spatial resolution, cross-matching galaxies at different wavelengths is challenging. This work aims to develop the next-generation deblended far-IR and sub-mm catalogues and present the first application in the COSMOS field. Our progressive deblending used the Bayesian probabilistic framework known as XID+. The deblending started from the Spitzer/MIPS 24 micron data, using an initial prior list composed of sources selected from the COSMOS2020 catalogue and radio catalogues from the VLA and the MeerKAT surveys, based on spectral energy distribution modelling which predicts fluxes of the known sources at the deblending wavelength. To speed up flux prediction, we made use of a neural network-based emulator. After deblending the 24 micron data, we proceeded to the Herschel PACS (100 & 160 micron) and SPIRE wavebands (250, 350 & 500 micron). Each time we constructed a tailor-made prior list based on the predicted fluxes of the known sources. Using simulated far-IR and sub-mm sky, we detailed the performance of our deblending pipeline. After validation with simulations, we then deblended the real observations from 24 to 500 micron and compared with blindly extracted catalogues and previous versions of deblended catalogues. As an additional test, we deblended the SCUBA-2 850 micron map and compared our deblended fluxes with ALMA measurements, which demonstrates a higher level of flux accuracy compared to previous results.We publicly release our XID+ deblended point source catalogues. These deblended long-wavelength data are crucial for studies such as deriving the fraction of dust-obscured star formation and better separation of quiescent galaxies from dusty star-forming galaxies. △ Less

Submitted 28 May, 2024; originally announced May 2024.

Comments: 23 pages, 30 figures, accepted for publication in A&A. Catalogues can be downloaded from https://hedam.lam.fr/HELP/dataproducts/dmu26/dmu26_XID+COSMOS2024/

Euclid. I. Overview of the Euclid mission

Authors: Euclid Collaboration, Y. Mellier, Abdurro'uf, J. A. Acevedo Barroso, A. Achúcarro, J. Adamek, R. Adam, G. E. Addison, N. Aghanim, M. Aguena, V. Ajani, Y. Akrami, A. Al-Bahlawan, A. Alavi, I. S. Albuquerque, G. Alestas, G. Alguero, A. Allaoui, S. W. Allen, V. Allevato, A. V. Alonso-Tetilla, B. Altieri, A. Alvarez-Candal, A. Amara, L. Amendola , et al. (1086 additional authors not shown)

Abstract: The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14… ▽ More The current standard model of cosmology successfully describes a variety of measurements, but the nature of its main ingredients, dark matter and dark energy, remains unknown. Euclid is a medium-class mission in the Cosmic Vision 2015-2025 programme of the European Space Agency (ESA) that will provide high-resolution optical imaging, as well as near-infrared imaging and spectroscopy, over about 14,000 deg^2 of extragalactic sky. In addition to accurate weak lensing and clustering measurements that probe structure formation over half of the age of the Universe, its primary probes for cosmology, these exquisite data will enable a wide range of science. This paper provides a high-level overview of the mission, summarising the survey characteristics, the various data-processing steps, and data products. We also highlight the main science objectives and expected performance. △ Less

Submitted 22 May, 2024; originally announced May 2024.

Comments: Paper submitted as part of the A&A special issue`Euclid on Sky'

Dark Energy Survey Year 3 results: simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps II. Cosmological results

Authors: M. Gatti, G. Campailla, N. Jeffrey, L. Whiteway, A. Porredon, J. Prat, J. Williamson, M. Raveri, B. Jain, V. Ajani, G. Giannini, M. Yamamoto, C. Zhou, J. Blazek, D. Anbajagane, S. Samuroff, T. Kacprzak, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, C. Chang, R. Chen , et al. (77 additional authors not shown)

Abstract: We present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years (Y3) of the Dark Energy Survey (DES). We implement: 1) second and third moments; 2) wavelet phase harmonics; 3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven $νw$CDM cosm… ▽ More We present a simulation-based cosmological analysis using a combination of Gaussian and non-Gaussian statistics of the weak lensing mass (convergence) maps from the first three years (Y3) of the Dark Energy Survey (DES). We implement: 1) second and third moments; 2) wavelet phase harmonics; 3) the scattering transform. Our analysis is fully based on simulations, spans a space of seven $νw$CDM cosmological parameters, and forward models the most relevant sources of systematics inherent in the data: masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration. We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a simulation-based inference approach. Including and combining different non-Gaussian statistics is a powerful tool that strongly improves constraints over Gaussian statistics (in our case, the second moments); in particular, the Figure of Merit $\textrm{FoM}(S_8, Ω_{\textrm{m}})$ is improved by 70 percent ($Λ$CDM) and 90 percent ($w$CDM). When all the summary statistics are combined, we achieve a 2 percent constraint on the amplitude of fluctuations parameter $S_8 \equiv σ_8 (Ω_{\textrm{m}}/0.3)^{0.5}$, obtaining $S_8 = 0.794 \pm 0.017$ ($Λ$CDM) and $S_8 = 0.817 \pm 0.021$ ($w$CDM). The constraints from different statistics are shown to be internally consistent (with a $p$-value>0.1 for all combinations of statistics examined). We compare our results to other weak lensing results from the DES Y3 data, finding good consistency; we also compare with results from external datasets, such as \planck{} constraints from the Cosmic Microwave Background, finding statistical agreement, with discrepancies no greater than $<2.2σ$. △ Less

Submitted 17 May, 2024; originally announced May 2024.

Comments: 24 pages, 13 figures, to be submitted to PRD. Comments welcome!

DEVILS/MIGHTEE/GAMA/DINGO: The Impact of SFR Timescales on the SFR-Radio Luminosity Correlation

Authors: Robin H. W. Cook, Luke J. M. Davies, Jonghwan Rhee, Catherine L. Hale, Sabine Bellstedt, Jessica E. Thorne, Ivan Delvecchio, Jordan D. Collier, Richard Dodson, Simon P. Driver, Benne W. Holwerda, Matt J. Jarvis, Kenda Knowles, Claudia Lagos, Natasha Maddox, Martin Meyer, Aaron S. G. Robotham, Sambit Roychowdhury, Kristof Rozgonyi, Nicholas Seymour, Malgorzata Siudek, Matthew Whiting, Imogen Whittam

Abstract: The tight relationship between infrared luminosity (L$_\mathrm{TIR}$) and 1.4 GHz radio continuum luminosity (L$_\mathrm{1.4GHz}$) has proven useful for understanding star formation free from dust obscuration. Infrared emission in star-forming galaxies typically arises from recently formed, dust-enshrouded stars, whereas radio synchrotron emission is expected from subsequent supernovae. By leverag… ▽ More The tight relationship between infrared luminosity (L$_\mathrm{TIR}$) and 1.4 GHz radio continuum luminosity (L$_\mathrm{1.4GHz}$) has proven useful for understanding star formation free from dust obscuration. Infrared emission in star-forming galaxies typically arises from recently formed, dust-enshrouded stars, whereas radio synchrotron emission is expected from subsequent supernovae. By leveraging the wealth of ancillary far-ultraviolet - far-infrared photometry from the Deep Extragalactic VIsible Legacy Survey (DEVILS) and Galaxy and Mass Assembly (GAMA) surveys, combined with 1.4 GHz observations from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) survey and Deep Investigation of Neutral Gas Origins (DINGO) projects, we investigate the impact of timescale differences between far-ultraviolet - far-infrared and radio-derived star formation rate (SFR) tracers. We examine how the SED-derived star formation histories (SFH) of galaxies can be used to explain discrepancies in these SFR tracers, which are sensitive to different timescales. Galaxies exhibiting an increasing SFH have systematically higher L$_\mathrm{TIR}$ and SED-derived SFRs than predicted from their 1.4 GHz radio luminosity. This indicates that insufficient time has passed for subsequent supernovae-driven radio emission to accumulate. We show that backtracking the SFR(t) of galaxies along their SED-derived SFHs to a time several hundred megayears prior to their observed epoch will both linearise the SFR-L$_\mathrm{1.4GHz}$ relation and reduce the overall scatter. The minimum scatter in the SFR(t)-L$_\mathrm{1.4GHz}$ is reached at 200 - 300 Myr prior, consistent with theoretical predictions for the timescales required to disperse the cosmic ray electrons responsible for the synchrotron emission. △ Less

Submitted 1 May, 2024; originally announced May 2024.

Comments: 23 pages, 13 figures, 2 tables. Accepted for publication in MNRAS

Weak lensing combined with the kinetic Sunyaev Zel'dovich effect: A study of baryonic feedback

Authors: L. Bigwood, A. Amon, A. Schneider, J. Salcido, I. G. McCarthy, C. Preston, D. Sanchez, D. Sijacki, E. Schaan, S. Ferraro, N. Battaglia, A. Chen, S. Dodelson, A. Roodman, A. Pieres, A. Ferte, A. Alarcon, A. Drlica-Wagner, A. Choi, A. Navarro-Alsina, A. Campos, A. J. Ross, A. Carnero Rosell, B. Yin, B. Yanny , et al. (100 additional authors not shown)

Abstract: Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmolo… ▽ More Extracting precise cosmology from weak lensing surveys requires modelling the non-linear matter power spectrum, which is suppressed at small scales due to baryonic feedback processes. However, hydrodynamical galaxy formation simulations make widely varying predictions for the amplitude and extent of this effect. We use measurements of Dark Energy Survey Year 3 weak lensing (WL) and Atacama Cosmology Telescope DR5 kinematic Sunyaev-Zel'dovich (kSZ) to jointly constrain cosmological and astrophysical baryonic feedback parameters using a flexible analytical model, `baryonification'. First, using WL only, we compare the $S_8$ constraints using baryonification to a simulation-calibrated halo model, a simulation-based emulator model and the approach of discarding WL measurements on small angular scales. We find that model flexibility can shift the value of $S_8$ and degrade the uncertainty. The kSZ provides additional constraints on the astrophysical parameters and shifts $S_8$ to $S_8=0.823^{+0.019}_{-0.020}$, a higher value than attained using the WL-only analysis. We measure the suppression of the non-linear matter power spectrum using WL + kSZ and constrain a mean feedback scenario that is more extreme than the predictions from most hydrodynamical simulations. We constrain the baryon fractions and the gas mass fractions and find them to be generally lower than inferred from X-ray observations and simulation predictions. We conclude that the WL + kSZ measurements provide a new and complementary benchmark for building a coherent picture of the impact of gas around galaxies across observations. △ Less

Submitted 9 April, 2024; originally announced April 2024.

Dark Energy Survey Year 3 results: likelihood-free, simulation-based $w$CDM inference with neural compression of weak-lensing map statistics

Authors: N. Jeffrey, L. Whiteway, M. Gatti, J. Williamson, J. Alsing, A. Porredon, J. Prat, C. Doux, B. Jain, C. Chang, T. -Y. Cheng, T. Kacprzak, P. Lemos, A. Alarcon, A. Amon, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Campos, A. Carnero Rosell, R. Chen, A. Choi, J. DeRose, A. Drlica-Wagner, K. Eckert , et al. (66 additional authors not shown)

Abstract: We present simulation-based cosmological $w$CDM inference using Dark Energy Survey Year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled… ▽ More We present simulation-based cosmological $w$CDM inference using Dark Energy Survey Year 3 weak-lensing maps, via neural data compression of weak-lensing map summary statistics: power spectra, peak counts, and direct map-level compression/inference with convolutional neural networks (CNN). Using simulation-based inference, also known as likelihood-free or implicit inference, we use forward-modelled mock data to estimate posterior probability distributions of unknown parameters. This approach allows all statistical assumptions and uncertainties to be propagated through the forward-modelled mock data; these include sky masks, non-Gaussian shape noise, shape measurement bias, source galaxy clustering, photometric redshift uncertainty, intrinsic galaxy alignments, non-Gaussian density fields, neutrinos, and non-linear summary statistics. We include a series of tests to validate our inference results. This paper also describes the Gower Street simulation suite: 791 full-sky PKDGRAV dark matter simulations, with cosmological model parameters sampled with a mixed active-learning strategy, from which we construct over 3000 mock DES lensing data sets. For $w$CDM inference, for which we allow $-1<w<-\frac{1}{3}$, our most constraining result uses power spectra combined with map-level (CNN) inference. Using gravitational lensing data only, this map-level combination gives $Ω_{\rm m} = 0.283^{+0.020}_{-0.027}$, ${S_8 = 0.804^{+0.025}_{-0.017}}$, and $w < -0.80$ (with a 68 per cent credible interval); compared to the power spectrum inference, this is more than a factor of two improvement in dark energy parameter ($Ω_{\rm DE}, w$) precision. △ Less

Submitted 4 March, 2024; originally announced March 2024.

Comments: 19 pages, 15 figures, submitted to Monthly Notices of the Royal Astronomical Society

MIGHTEE-HI: HI galaxy properties in the large scale structure environment at z~0.37 from a stacking experiment

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Alessandro Bianchetti, Mattia Vaccari, Natasha Maddox, Anastasia A. Ponomareva, Bradley S. Frank, Matt J. Jarvis, Barbara Catinella, Luca Cortese, Sambit Roychowdhury, Maarten Baes, Jordan D. Collier, Olivier Ilbert, Ali A. Khostovan, Sushma Kurapati, Hengxing Pan, Isabella Prandoni, Sambatriniaina H. A. Rajohnson, Mara Salvato, Srikrishna Sekhar, Gauri Sharma

Abstract: We present the first measurement of HI mass of star-forming galaxies in different large scale structure environments from a blind survey at $z\sim 0.37$. In particular, we carry out a spectral line stacking analysis considering $2875$ spectra of colour-selected star-forming galaxies undetected in HI at $0.23 < z < 0.49$ in the COSMOS field, extracted from the MIGHTEE-HI Early Science datacubes, ac… ▽ More We present the first measurement of HI mass of star-forming galaxies in different large scale structure environments from a blind survey at $z\sim 0.37$. In particular, we carry out a spectral line stacking analysis considering $2875$ spectra of colour-selected star-forming galaxies undetected in HI at $0.23 < z < 0.49$ in the COSMOS field, extracted from the MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies belonging to different subsamples depending on three different definitions of large scale structure environment: local galaxy overdensity, position inside the host dark matter halo (central, satellite, or isolated), and cosmic web type (field, filament, or knot). We first stack the full star-forming galaxy sample and find a robust HI detection yielding an average galaxy HI mass of $M_{\rm HI}=(8.12\pm 0.75)\times 10^9\, {\rm M}_\odot$ at $\sim 11.8σ$. Next, we investigate the different subsamples finding a negligible difference in $M_{\rm HI}$ as a function of the galaxy overdensity. We report an HI excess compared to the full sample in satellite galaxies ($M_{\rm HI}=(11.31\pm1.22)\times 10^9$, at $\sim 10.2 σ$) and in filaments ($M_{\rm HI}=(11.62\pm 0.90)\times 10^9$. Conversely, we report non-detections for the central and knot galaxies subsamples, which appear to be HI-deficient. We find the same qualitative results also when stacking in units of HI fraction ($f_{\rm HI}$). We conclude that the HI amount in star-forming galaxies at the studied redshifts correlates with the large scale structure environment. △ Less

Submitted 1 March, 2024; originally announced March 2024.

Comments: Accepted for publication in MNRAS. 15 figures, 3 tables

The sizes of bright Lyman-break galaxies at $z\simeq3-5$ with JWST PRIMER

Authors: R. G. Varadaraj, R. A. A. Bowler, M. J. Jarvis, N. J. Adams, N. Choustikov, A. M. Koekemoer, A. C. Carnall, D. J. McLeod, J. S. Dunlop, C. T. Donnan, N. A. Grogin

Abstract: We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distribution… ▽ More We use data from the JWST Public Release IMaging for Extragalactic Research (PRIMER) survey to measure the size scaling relations of 1668 rest-frame UV-bright Lyman-break galaxies (LBGs) at $z=3-5$ with stellar masses $\mathrm{log}_{10}(M_{\star}/M_{\odot}) > 9$. The sample was selected from seeing-dominated ground-based data, presenting an unbiased sampling of the morphology and size distributions of luminous sources. We fit Sérsic profiles to eight NIRCam bands and also measure a non-parametric half-light radius. We find that the size distributions with both measurements are well-fit by a log-normal distribution at all redshifts, consistent with disk formation models where size is governed by host dark-matter halo angular momentum. We find a size-redshift evolution of $R_{e} = 3.51(1+z)^{-0.60\pm0.22}$ kpc, in agreement with JWST studies. When considering the typical (modal) size over $z=3-5$, we find little evolution with bright LBGs remaining compact at $R_{e}\simeq0.7-0.9$ kpc. Simultaneously, we find evidence for a build-up of large ($R_{e} > 2$ kpc) galaxies by $z=3$. We find some evidence for a negatively sloped size-mass relation at $z=5$ when Sérsic profiles are used to fit the data in F200W. The intrinsic scatter in our size-mass relations increases at higher redshifts. Additionally, measurements probing the rest-UV (F200W) show larger intrinsic scatter than those probing the rest-optical (F356W). Finally, we leverage rest-UV and rest-optical photometry to show that disky galaxies are well established by $z=5$, but are beginning to undergo dissipative processes, such as mergers, by $z=3$. The agreement of our size-mass and size-luminosity relations with simulations provides tentative evidence for centrally concentrated star formation at high-redshift. △ Less

Submitted 29 January, 2024; originally announced January 2024.

Comments: 15 pages, 9 figures, 3 tables. Submitted to MNRAS

The VLBA CANDELS GOODS-North Survey. I -- Survey Design, Processing, Data Products, and Source Counts

Authors: Roger P. Deane, Jack F. Radcliffe, Ann Njeri, Alexander Akoto-Danso, Gianni Bernardi, Oleg M. Smirnov, Rob Beswick, Michael A. Garrett, Matt J. Jarvis, Imogen H. Whittam, Stephen Bourke, Zsolt Paragi

Abstract: The past decade has seen significant advances in wide-field cm-wave very long baseline interferometry (VLBI), which is timely given the wide-area, synoptic survey-driven strategy of major facilities across the electromagnetic spectrum. While wide-field VLBI poses significant post-processing challenges that can severely curtail its potential scientific yield, many developments in the km-scale conne… ▽ More The past decade has seen significant advances in wide-field cm-wave very long baseline interferometry (VLBI), which is timely given the wide-area, synoptic survey-driven strategy of major facilities across the electromagnetic spectrum. While wide-field VLBI poses significant post-processing challenges that can severely curtail its potential scientific yield, many developments in the km-scale connected-element interferometer sphere are directly applicable to addressing these. Here we present the design, processing, data products, and source counts from a deep (11 $μ$Jy beam$^{-1}$), quasi-uniform sensitivity, contiguous wide-field (160 arcmin$^2$) 1.6 GHz VLBI survey of the CANDELS GOODS-North field. This is one of the best-studied extragalactic fields at milli-arcsecond resolution and, therefore, is well-suited as a comparative study for our Tera-pixel VLBI image. The derived VLBI source counts show consistency with those measured in the COSMOS field, which broadly traces the AGN population detected in arcsecond-scale radio surveys. However, there is a distinctive flattening in the $ S_{\rm 1.4GHz}\sim$100-500 $μ$Jy flux density range, which suggests a transition in the population of compact faint radio sources, qualitatively consistent with the excess source counts at 15 GHz that is argued to be an unmodelled population of radio cores. This survey approach will assist in deriving robust VLBI source counts and broadening the discovery space for future wide-field VLBI surveys, including VLBI with the Square Kilometre Array, which will include new large field-of-view antennas on the African continent at $\gtrsim$1000~km baselines. In addition, it may be useful in the design of both monitoring and/or rapidly triggered VLBI transient programmes. △ Less

Submitted 22 January, 2024; originally announced January 2024.

Comments: 15 pages, 8 figures. Accepted for publication in Monthly Notices of the Royal Astronomical Society

SPT Clusters with DES and HST Weak Lensing. II. Cosmological Constraints from the Abundance of Massive Halos

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, T. Schrabback, T. M. C. Abbott, P. A. R. Ade, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, A. J. Anderson, J. Annis, B. Ansarinejad, J. E. Austermann, S. Avila, D. Bacon, M. Bayliss, J. A. Beall, K. Bechtol, M. R. Becker, A. N. Bender , et al. (171 additional authors not shown)

Abstract: 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… ▽ More 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<z<1.7$. The weak-lensing measurements enable robust mass measurements of sample clusters and allow us to empirically constrain the SZ observable--mass relation. For a flat $Λ$CDM cosmology, and marginalizing over the sum of massive neutrinos, we measure $Ω_\mathrm{m}=0.286\pm0.032$, $σ_8=0.817\pm0.026$, and the parameter combination $σ_8\,(Ω_\mathrm{m}/0.3)^{0.25}=0.805\pm0.016$. Our measurement of $S_8\equivσ_8\,\sqrt{Ω_\mathrm{m}/0.3}=0.795\pm0.029$ and the constraint from Planck CMB anisotropies (2018 TT,TE,EE+lowE) differ by $1.1σ$. In combination with that Planck dataset, we place a 95% upper limit on the sum of neutrino masses $\sum m_ν<0.18$ eV. When additionally allowing the dark energy equation of state parameter $w$ to vary, we obtain $w=-1.45\pm0.31$ from our cluster-based analysis. In combination with Planck data, we measure $w=-1.34^{+0.22}_{-0.15}$, or a $2.2σ$ difference with a cosmological constant. We use the cluster abundance to measure $σ_8$ in five redshift bins between 0.25 and 1.8, and we find the results to be consistent with structure growth as predicted by the $Λ$CDM model fit to Planck primary CMB data. △ Less

Submitted 21 June, 2024; v1 submitted 4 January, 2024; originally announced January 2024.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

MIGHTEE Polarization Early Science Fields: The Deep Polarized Sky

Authors: A. R. Taylor, S. Sekhar, L. Heino, A. M. M. Scaife, J. Stil, M. Bowles, M. Jarvis, I. Heywood, J. D. Collier

Abstract: The MeerKAT International GigaHertz Tiered Extragalactic Exploration (MIGHTEE) is one of the MeerKAT large survey projects, designed to pathfind SKA key science. MIGHTEE is undertaking deep radio imaging of four well observed fields (COSMOS, XMM-LSS, ELAIS S1 and CDFS) totaling 20 square degrees to $μ$Jy sensitivities. Broadband imaging observations between 880--1690 MHz yield total intensity cont… ▽ More The MeerKAT International GigaHertz Tiered Extragalactic Exploration (MIGHTEE) is one of the MeerKAT large survey projects, designed to pathfind SKA key science. MIGHTEE is undertaking deep radio imaging of four well observed fields (COSMOS, XMM-LSS, ELAIS S1 and CDFS) totaling 20 square degrees to $μ$Jy sensitivities. Broadband imaging observations between 880--1690 MHz yield total intensity continuum, spectro-polarimetry, and atomic hydrogen spectral imaging. Early science data from MIGHTEE are being released from initial observations of COSMOS and XMM-LSS. This paper describes the spectro-polarimetric observations, the polarization data processing of the MIGHTEE early science fields, and presents polarization data images and catalogues. The catalogues include radio spectral index, redshift information and Faraday rotation measure synthesis results for 13,271 total intensity radio sources down to a polarized intensity detection limit of $\sim$20 $μ$Jy\,bm$^{-1}$. Polarized signals were detected from 324 sources. For the polarized detections we include a catalogue of Faraday Depth from both Faraday Synthesis and $Q$, $U$ fitting, as well as total intensity and polarization spectral indices. The distribution of redshift of the total radio sources and detected polarized sources are the same, with median redshifts of 0.86 and 0.82 respectively. Depolarization of the emission at longer-wavelengths is seen to increase with decreasing total-intensity spectral index, implying that depolarisation is intrinsic to the radio sources. No evidence is seen for a redshift dependence of the variance of Faraday Depth. △ Less

Submitted 20 December, 2023; originally announced December 2023.

Comments: 12 pages, 16 figures

The discovery of a z=0.7092 OH megamaser with the MIGHTEE survey

Authors: Matt J. Jarvis, Ian Heywood, Sophie M. Jewell, Roger P. Deane, H. -R. Klöckner, Anastasia A. Ponomareva, Natasha Maddox, Andrew J. Baker, Alessandro Bianchetti, Kelley M. Hess, Hayley Roberts, Giulia Rodighiero, Ilaria Ruffa, Francesco Sinigaglia, R. G. Varadaraj, I. H. Whittam, Elizabeth A. K. Adams, Maarten Baes, Eric J. Murphy, Hengxing Pan, Mattia Vaccari

Abstract: We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both… ▽ More We present the discovery of the most distant OH megamaser to be observed in the main lines, using data from the MeerKAT International Giga-Hertz Tiered Extragalactic Exploration (MIGHTEE) survey. At a newly measured redshift of $z = 0.7092$, the system has strong emission in both the 1665MHz ($L \approx 2500$ L$_{\odot}$) and 1667 MHz ($L \approx 4.5\times10^4$ L$_{\odot}$) transitions, with both narrow and broad components. We interpret the broad line as a high-velocity-dispersion component of the 1667 MHz transition, with velocity $v \sim 330$km s$^{-1}$ with respect to the systemic velocity. The host galaxy has a stellar mass of $M_{\star} = 2.95 \times 10^{10}$ M$_{\odot}$ and a star-formation rate of SFR = 371 M$_{\odot}$yr$^{-1}$, placing it $\sim 1.5$dex above the main sequence for star-forming galaxies at this redshift, and can be classified as an ultra-luminous infrared galaxy. Alongside the optical imaging data, which exhibits evidence for a tidal tail, this suggests that the OH megamaser arises from a system that is currently undergoing a merger, which is stimulating star formation and providing the necessary conditions for pumping the OH molecule to saturation. The OHM is likely to be lensed, with a magnification factor of $\sim 2.5$, and perhaps more if the maser emitting region is compact and suitably offset relative to the centroid of its host galaxy's optical light. This discovery demonstrates that spectral line mapping with the new generation of radio interferometers may provide important information on the cosmic merger history of galaxies. △ Less

Submitted 7 December, 2023; originally announced December 2023.

Comments: 11 pages, accepted for publication in MNRAS

Extragalactic Magnetism with SOFIA (SALSA Legacy Program). VII. A Tomographic View of Far-infrared and Radio Polarimetric Observations through MHD Simulations of Galaxies

Authors: Sergio Martin-Alvarez, Enrique Lopez-Rodriguez, Tara Dacunha, Susan E. Clark, Alejandro S. Borlaff, Rainer Beck, Francisco Rodríguez Montero, S. Lyla Jung, Julien Devriendt, Adrianne Slyz, Julia Roman-Duval, Evangelia Ntormousi, Mehrnoosh Tahani, Kandaswamy Subramanian, Daniel A. Dale, Pamela M. Marcum, Konstantinos Tassis, Ignacio del Moral-Castro, Le Ngoc Tram, Matt J. Jarvis

Abstract: The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy t… ▽ More The structure of magnetic fields in galaxies remains poorly constrained, despite the importance of magnetism in the evolution of galaxies. Radio synchrotron and far-infrared (FIR) polarization and polarimetric observations are the best methods to measure galactic scale properties of magnetic fields in galaxies beyond the Milky Way. We use synthetic polarimetric observations of a simulated galaxy to identify and quantify the regions, scales, and interstellar medium (ISM) phases probed at FIR and radio wavelengths. Our studied suite of magnetohydrodynamical cosmological zoom-in simulations features high-resolutions (10 pc full-cell size) and multiple magnetization models. Our synthetic observations have a striking resemblance to those of observed galaxies. We find that the total and polarized radio emission extends to approximately double the altitude above the galactic disk (half-intensity disk thickness of $h_\text{I radio} \sim h_\text{PI radio} = 0.23 \pm 0.03$ kpc) relative to the total FIR and polarized emission that are concentrated in the disk midplane ($h_\text{I FIR} \sim h_\text{PI FIR} = 0.11 \pm 0.01$ kpc). Radio emission traces magnetic fields at scales of $\gtrsim 300$ pc, whereas FIR emission probes magnetic fields at the smallest scales of our simulations. These scales are comparable to our spatial resolution and well below the spatial resolution ($<300$ pc) of existing FIR polarimetric measurements. Finally, we confirm that synchrotron emission traces a combination of the warm neutral and cold neutral gas phases, whereas FIR emission follows the densest gas in the cold neutral phase in the simulation. These results are independent of the ISM magnetic field strength. The complementarity we measure between radio and FIR wavelengths motivates future multiwavelength polarimetric observations to advance our knowledge of extragalactic magnetism. △ Less

Submitted 24 March, 2024; v1 submitted 10 November, 2023; originally announced November 2023.

Comments: Accepted for publication in the ApJ. Replaced to match final version. 35 pages, 17 figures, 4 tables

Dark Energy Survey Year 3 results: simulation-based cosmological inference with wavelet harmonics, scattering transforms, and moments of weak lensing mass maps I: validation on simulations

Authors: M. Gatti, N. Jeffrey, L. Whiteway, J. Williamson, B. Jain, V. Ajani, D. Anbajagane, G. Giannini, C. Zhou, A. Porredon, J. Prat, M. Yamamoto, J. Blazek, T. Kacprzak, S. Samuroff, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, C. Chang, R. Chen, A. Choi, C. Davis , et al. (76 additional authors not shown)

Abstract: Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulation-based analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogues of the Dark Energy Survey. We implement: 1) second and third moments; 2) wavelet phase… ▽ More Beyond-two-point statistics contain additional information on cosmological as well as astrophysical and observational (systematics) parameters. In this methodology paper we provide an end-to-end simulation-based analysis of a set of Gaussian and non-Gaussian weak lensing statistics using detailed mock catalogues of the Dark Energy Survey. We implement: 1) second and third moments; 2) wavelet phase harmonics (WPH); 3) the scattering transform (ST). Our analysis is fully based on simulations, it spans a space of seven $νw$CDM cosmological parameters, and it forward models the most relevant sources of systematics of the data (masks, noise variations, clustering of the sources, intrinsic alignments, and shear and redshift calibration). We implement a neural network compression of the summary statistics, and we estimate the parameter posteriors using a likelihood-free-inference approach. We validate the pipeline extensively, and we find that WPH exhibits the strongest performance when combined with second moments, followed by ST. and then by third moments. The combination of all the different statistics further enhances constraints with respect to second moments, up to 25 per cent, 15 per cent, and 90 per cent for $S_8$, $Ω_{\rm m}$, and the Figure-Of-Merit ${\rm FoM_{S_8,Ω_{\rm m}}}$, respectively. We further find that non-Gaussian statistics improve constraints on $w$ and on the amplitude of intrinsic alignment with respect to second moments constraints. The methodological advances presented here are suitable for application to Stage IV surveys from Euclid, Rubin-LSST, and Roman with additional validation on mock catalogues for each survey. In a companion paper we present an application to DES Year 3 data. △ Less

Submitted 4 November, 2023; v1 submitted 26 October, 2023; originally announced October 2023.

Comments: 25 pages, 18 figures. Comments welcome!

MIGHTEE: multi-wavelength counterparts in the COSMOS field

Authors: I. H. Whittam, M. Prescott, C. L. Hale, M . J. Jarvis, I. Heywood, Fangxia An, M. Glowacki, N. Maddox, L. Marchetti, L. K. Morabito, N. J. Adams, R. A. A. Bowler, P. W. Hatfield, R. G. Varadaraj, J. Collier, B. Frank, A. R. Taylor, M. G. Santos, M. Vaccari, J. Afonso, Y. Ao, J. Delhaize, K. Knowles, S. Kolwa, S. M. Randriamampandry , et al. (4 additional authors not shown)

Abstract: In this paper we combine the Early Science radio continuum data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey, with optical and near-infrared data and release the cross-matched catalogues. The radio data used in this work covers $0.86$ deg$^2$ of the COSMOS field, reaches a thermal noise of $1.7$ $μ$Jy/beam and contains $6102$ radio components. We visually in… ▽ More In this paper we combine the Early Science radio continuum data from the MeerKAT International GHz Tiered Extragalactic Exploration (MIGHTEE) Survey, with optical and near-infrared data and release the cross-matched catalogues. The radio data used in this work covers $0.86$ deg$^2$ of the COSMOS field, reaches a thermal noise of $1.7$ $μ$Jy/beam and contains $6102$ radio components. We visually inspect and cross-match the radio sample with optical and near-infrared data from the Hyper Suprime-Cam (HSC) and UltraVISTA surveys. This allows the properties of active galactic nuclei and star-forming populations of galaxies to be probed out to $z \approx 5$. Additionally, we use the likelihood ratio method to automatically cross-match the radio and optical catalogues and compare this to the visually cross-matched catalogue. We find that 94 per cent of our radio source catalogue can be matched with this method, with a reliability of $95$ per cent. We proceed to show that visual classification will still remain an essential process for the cross-matching of complex and extended radio sources. In the near future, the MIGHTEE survey will be expanded in area to cover a total of $\sim$20~deg$^2$; thus the combination of automated and visual identification will be critical. We compare redshift distribution of SFG and AGN to the SKADS and T-RECS simulations and find more AGN than predicted at $z \sim 1$. △ Less

Submitted 26 October, 2023; originally announced October 2023.

Comments: 15 pages, 15 figures. Accepted for publication in MNRAS

SPT Clusters with DES and HST Weak Lensing. I. Cluster Lensing and Bayesian Population Modeling of Multi-Wavelength Cluster Datasets

Authors: S. Bocquet, S. Grandis, L. E. Bleem, M. Klein, J. J. Mohr, M. Aguena, A. Alarcon, S. Allam, S. W. Allen, O. Alves, A. Amon, B. Ansarinejad, D. Bacon, M. Bayliss, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, M. Brodwin, D. Brooks, A. Campos, R. E. A. Canning, J. E. Carlstrom, A. Carnero Rosell, M. Carrasco Kind , et al. (108 additional authors not shown)

Abstract: 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 calibrati… ▽ More 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. △ Less

Submitted 21 June, 2024; v1 submitted 18 October, 2023; originally announced October 2023.

Comments: Accepted for publication in Phys. Rev. D. arXiv v2 corresponds to published article

Cosmology from LOFAR Two-metre Sky Survey Data Release 2: Cross-correlation with the cosmic microwave background

Authors: S. J. Nakoneczny, D. Alonso, M. Bilicki, D. J. Schwarz, C. L. Hale, A. Pollo, C. Heneka, P. Tiwari, J. Zheng, M. Brüggen, M. J. Jarvis, T. W. Shimwell

Abstract: We combine the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) catalogue with gravitational lensing maps from the Cosmic Microwave Background (CMB) to place constraints on the bias evolution of LoTSS radio galaxies, and on the amplitude of matter perturbations. We construct a flux-limited catalogue, and analyse its harmonic-space cross-correlation with CMB lensin… ▽ More We combine the LOw-Frequency ARray (LOFAR) Two-metre Sky Survey (LoTSS) second data release (DR2) catalogue with gravitational lensing maps from the Cosmic Microwave Background (CMB) to place constraints on the bias evolution of LoTSS radio galaxies, and on the amplitude of matter perturbations. We construct a flux-limited catalogue, and analyse its harmonic-space cross-correlation with CMB lensing maps from Planck, $C_\ell^{gκ}$, as well as its auto-correlation, $C_\ell^{gg}$. We explore the models describing the redshift evolution of the large-scale radio galaxy bias, discriminating between them through the combination of both $C_\ell^{gκ}$ and $C_\ell^{gg}$. Fixing the bias evolution, we then use these data to place constraints on the amplitude of large scale density fluctuations. We report the significance of the $C_\ell^{gκ}$ signal at a level of $26.6σ$. We determine that a linear bias evolution of the form $b_g(z) = b_{g,D} / D(z)$, where $D(z)$ is the growth rate, is able to provide a good description of the data, and measure $b_{g,D} = 1.41 \pm 0.06$ for a sample flux-limited at $1.5\,{\rm mJy}$, for scales $\ell < 250$ for $C_\ell^{gg}$, and $\ell < 500$ for $C_\ell^{gκ}$. At the sample's median redshift, we obtain $b(z = 0.82) = 2.34 \pm 0.10$. Using $σ_8$ as a free parameter, while keeping other cosmological parameters fixed to the Planck values, we find fluctuations of $σ_8 = 0.75^{+0.05}_{-0.04}$. The result is in agreement with weak lensing surveys, and at $1σ$ difference with Planck CMB constraints. We also attempt to detect the late-time integrated Sachs-Wolfe effect with LOFAR, but with the current sky coverage, the cross-correlation with CMB temperature maps is consistent with zero. Our results are an important step towards constraining cosmology with radio continuum surveys from LOFAR and other future large radio surveys. △ Less

Submitted 15 May, 2024; v1 submitted 11 October, 2023; originally announced October 2023.

Comments: The code is available at https://github.com/snakoneczny/cosmo-pipe

Journal ref: A&A, 681, A105 (2024)

Cosmology from LOFAR Two-metre Sky Survey Data Release 2: Angular Clustering of Radio Sources

Authors: C. L. Hale, D. J. Schwarz, P. N. Best, S. J. Nakoneczny, D. Alonso, D. Bacon, L. Böhme, N. Bhardwaj, M. Bilicki, S. Camera, C. S. Heneka, M. Pashapour-Ahmadabadi, P. Tiwari, J. Zheng, K. J. Duncan, M. J. Jarvis, R. Kondapally, M. Magliocchetti, H. J. A. Rottgering, T. W. Shimwell

Abstract: Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies wit… ▽ More Covering $\sim$5600 deg$^2$ to rms sensitivities of $\sim$70$-$100 $μ$Jy beam$^{-1}$, the LOFAR Two-metre Sky Survey Data Release 2 (LoTSS-DR2) provides the largest low-frequency ($\sim$150 MHz) radio catalogue to date, making it an excellent tool for large-area radio cosmology studies. In this work, we use LoTSS-DR2 sources to investigate the angular two-point correlation function of galaxies within the survey. We discuss systematics in the data and an improved methodology for generating random catalogues, compared to that used for LoTSS-DR1, before presenting the angular clustering for $\sim$900,000 sources $\geq$$1.5$ mJy and a peak signal-to-noise $\geq$$7.5$ across $\sim$$80\%$ of the observed area. Using the clustering we infer the bias assuming two evolutionary models. When fitting {angular scales of $0.5 \leqθ<5\,°$, using a linear bias model, we find LoTSS-DR2 sources are biased tracers of the underlying matter, with a bias of $b_{C}= 2.14^{+0.22}_{-0.20}$ (assuming constant bias) and $b_{E}(z=0)= 1.79^{+0.15}_{-0.14}$ (for an evolving model, inversely proportional to the growth factor), corresponding to $b_E= 2.81^{+0.24}_{-0.22}$ at the median redshift of our sample, assuming the LoTSS Deep Fields redshift distribution is representative of our data. This reduces to $b_{C}= 2.02^{+0.17}_{-0.16}$ and $b_{E}(z=0)= 1.67^{+0.12}_{-0.12}$ when allowing preferential redshift distributions from the Deep Fields to model our data. Whilst the clustering amplitude is slightly lower than LoTSS-DR1 ($\geq$2 mJy), our study benefits from larger samples and improved redshift estimates. △ Less

Submitted 11 October, 2023; originally announced October 2023.

Comments: Accepted for publication in MNRAS. 29 pages, 24 figures

Cosmological shocks around galaxy clusters: A coherent investigation with DES, SPT & ACT

Authors: D. Anbajagane, C. Chang, E. J. Baxter, S. Charney, M. Lokken, M. Aguena, S. Allam, O. Alves, A. Amon, R. An, F. Andrade-Oliveira, D. Bacon, N. Battaglia, K. Bechtol, M. R. Becker, B. A. Benson, G. M. Bernstein, L. Bleem, S. Bocquet, J. R. Bond, D. Brooks, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi , et al. (89 additional authors not shown)

Abstract: We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around… ▽ More We search for signatures of cosmological shocks in gas pressure profiles of galaxy clusters using the cluster catalogs from three surveys: the Dark Energy Survey (DES) Year 3, the South Pole Telescope (SPT) SZ survey, and the Atacama Cosmology Telescope (ACT) data releases 4, 5, and 6, and using thermal Sunyaev-Zeldovich (SZ) maps from SPT and ACT. The combined cluster sample contains around $10^5$ clusters with mass and redshift ranges $10^{13.7} < M_{\rm 200m}/M_\odot < 10^{15.5}$ and $0.1 < z < 2$, and the total sky coverage of the maps is $\approx 15,000 \,\,{\rm deg}^2$. We find a clear pressure deficit at $R/R_{\rm 200m}\approx 1.1$ in SZ profiles around both ACT and SPT clusters, estimated at $6σ$ significance, which is qualitatively consistent with a shock-induced thermal non-equilibrium between electrons and ions. The feature is not as clearly determined in profiles around DES clusters. We verify that measurements using SPT or ACT maps are consistent across all scales, including in the deficit feature. The SZ profiles of optically selected and SZ-selected clusters are also consistent for higher mass clusters. Those of less massive, optically selected clusters are suppressed on small scales by factors of 2-5 compared to predictions, and we discuss possible interpretations of this behavior. An oriented stacking of clusters -- where the orientation is inferred from the SZ image, the brightest cluster galaxy, or the surrounding large-scale structure measured using galaxy catalogs -- shows the normalization of the one-halo and two-halo terms vary with orientation. Finally, the location of the pressure deficit feature is statistically consistent with existing estimates of the splashback radius. △ Less

Submitted 12 December, 2023; v1 submitted 29 September, 2023; originally announced October 2023.

Comments: [v2]: Version accepted to MNRAS

A Differentiable Perturbation-based Weak Lensing Shear Estimator

Authors: Xiangchong Li, Rachel Mandelbaum, Mike Jarvis, Yin Li, Andy Park, Tianqing Zhang

Abstract: Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-percent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (FPFS), using JAX, which is publicly available as part of a new suite of analytic shear algorithms called AnaCal. This code can… ▽ More Upcoming imaging surveys will use weak gravitational lensing to study the large-scale structure of the Universe, demanding sub-percent accuracy for precise cosmic shear measurements. We present a new differentiable implementation of our perturbation-based shear estimator (FPFS), using JAX, which is publicly available as part of a new suite of analytic shear algorithms called AnaCal. This code can analytically calibrate the shear response of any nonlinear observable constructed with the FPFS shapelets and detection modes utilizing auto-differentiation (AD), generalizing the formalism to include a family of shear estimators with corrections for detection and selection biases. Using the AD capability of JAX, it calculates the full Hessian matrix of the non-linear observables, which improves the previously presented second-order noise bias correction in the shear estimation. As an illustration of the power of the new AnaCal framework, we optimize the effective galaxy number density in the space of the generalized shear estimators using an LSST-like galaxy image simulation for the ten-year LSST. For the generic shear estimator, the magnitude of the multiplicative bias $|m|$ is below $3\times 10^{-3}$ (99.7% confidence interval), and the effective galaxy number density is improved by 5%. We also discuss some planned future additions to the AnaCal software suite to extend its applicability beyond the FPFS measurements. △ Less

Submitted 21 December, 2023; v1 submitted 12 September, 2023; originally announced September 2023.

Comments: 9 pages, 7 figures, Accepted for publication in MNRAS

Cosmology from Cross-Correlation of ACT-DR4 CMB Lensing and DES-Y3 Cosmic Shear

Authors: S. Shaikh, I. Harrison, A. van Engelen, G. A. Marques, T. M. C. Abbott, M. Aguena, O. Alves, A. Amon, R. An, D. Bacon, N. Battaglia, M. R. Becker, G. M. Bernstein, E. Bertin, J. Blazek, J. R. Bond, D. Brooks, D. L. Burke, E. Calabrese, A. Carnero Rosell, J. Carretero, R. Cawthon, C. Chang, R. Chen, A. Choi , et al. (83 additional authors not shown)

Abstract: Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy… ▽ More Cross-correlation between weak lensing of the Cosmic Microwave Background (CMB) and weak lensing of galaxies offers a way to place robust constraints on cosmological and astrophysical parameters with reduced sensitivity to certain systematic effects affecting individual surveys. We measure the angular cross-power spectrum between the Atacama Cosmology Telescope (ACT) DR4 CMB lensing and the galaxy weak lensing measured by the Dark Energy Survey (DES) Y3 data. Our baseline analysis uses the CMB convergence map derived from ACT-DR4 and $\textit{Planck}$ data, where most of the contamination due to the thermal Sunyaev Zel'dovich effect is removed, thus avoiding important systematics in the cross-correlation. In our modelling, we consider the nuisance parameters of the photometric uncertainty, multiplicative shear bias and intrinsic alignment of galaxies. The resulting cross-power spectrum has a signal-to-noise ratio $= 7.1$ and passes a set of null tests. We use it to infer the amplitude of the fluctuations in the matter distribution ($S_8 \equiv σ_8 (Ω_{\rm m}/0.3)^{0.5} = 0.782\pm 0.059$) with informative but well-motivated priors on the nuisance parameters. We also investigate the validity of these priors by significantly relaxing them and checking the consistency of the resulting posteriors, finding them consistent, albeit only with relatively weak constraints. This cross-correlation measurement will improve significantly with the new ACT-DR6 lensing map and form a key component of the joint 6x2pt analysis between DES and ACT. △ Less

Submitted 8 September, 2023; originally announced September 2023.

Comments: 26 pages, 30 figures (including appendices). Data associated with this article is available at https://github.com/itrharrison/actdr4kappa-x-desy3gamma-data

Report number: FERMILAB-PUB-23-432-PPD

Probing magnetic fields in the circumgalactic medium using polarization data from MIGHTEE

Authors: K. Böckmann, M. Brüggen, V. Heesen, A. Basu, S. P. O'Sullivan, I. Heywood, M. Jarvis, A. Scaife, J. Stil, R. Taylor, N. J. Adams, R. A. A. Bowler, M. N. Tudorache

Abstract: The detection and study of magnetic fields surrounding galaxies is important to understand galaxy evolution since magnetic fields are tracers for dynamical processes in the circumgalactic medium (CGM) and can have a significant impact on the evolution of the CGM. The Faraday rotation measure (RM) of the polarized light of background radio sources passing through the magnetized CGM of intervening g… ▽ More The detection and study of magnetic fields surrounding galaxies is important to understand galaxy evolution since magnetic fields are tracers for dynamical processes in the circumgalactic medium (CGM) and can have a significant impact on the evolution of the CGM. The Faraday rotation measure (RM) of the polarized light of background radio sources passing through the magnetized CGM of intervening galaxies can be used as a tracer for the strength and extent of magnetic fields around galaxies. We use rotation measures observed by the MIGHTEE-POL (MeerKAT International GHz Tiered Extragalactic Exploration POLarisation) survey by MeerKAT in the XMM-LSS and COSMOS fields to investigate the RM around foreground star-forming galaxies. We use spectroscopic catalogs of star-forming and blue cloud galaxies to measure the RM of MIGHTEE-POL sources as a function of the impact parameter from the intervening galaxy. We then repeat this procedure using a deeper galaxy catalog with photometric redshifts. For the spectroscopic star-forming sample we find a redshift-corrected |RM| excess of 5.6 +/- 2.3 rad m-2 which corresponds to a 2.5 sigma significance around galaxies with a median redshift of z = 0.46 for impact parameters below 130 kpc only selecting the intervenor with the smallest impact parameter. Making use of a photometric galaxy catalog and taking into account all intervenors with Mg < -13.6 mag, the signal disappears. We find no indication for a correlation between redshift and RM, nor do we find a connection between the total number of intervenors to the total |RM| . We have presented tentative evidence that the CGM of star-forming galaxies is permeated by coherent magnetic fields within the virial radius. We conclude that mostly bright, star-forming galaxies with impact parameters less than 130 kpc significantly contribute to the RM of the background radio source. △ Less

Submitted 22 August, 2023; originally announced August 2023.

Comments: 11 pages, 8 figures, accepted for publication in A&A

Journal ref: A&A 678, A56 (2023)

Beyond the 3rd moment: A practical study of using lensing convergence CDFs for cosmology with DES Y3

Authors: D. Anbajagane, C. Chang, A. Banerjee, T. Abel, M. Gatti, V. Ajani, A. Alarcon, A. Amon, E. J. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Chen, A. Choi, C. Davis, J. DeRose, H. T. Diehl, S. Dodelson, C. Doux, A. Drlica-Wagner, K. Eckert, J. Elvin-Poole , et al. (73 additional authors not shown)

Abstract: Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at… ▽ More Widefield surveys of the sky probe many clustered scalar fields -- such as galaxy counts, lensing potential, gas pressure, etc. -- that are sensitive to different cosmological and astrophysical processes. Our ability to constrain such processes from these fields depends crucially on the statistics chosen to summarize the field. In this work, we explore the cumulative distribution function (CDF) at multiple scales as a summary of the galaxy lensing convergence field. Using a suite of N-body lightcone simulations, we show the CDFs' constraining power is modestly better than that of the 2nd and 3rd moments of the field, as they approximately capture the information from all moments of the field in a concise data vector. We then study the practical aspects of applying the CDFs to observational data, using the first three years of the Dark Energy Survey (DES Y3) data as an example, and compute the impact of different systematics on the CDFs. The contributions from the point spread function are 2-3 orders of magnitude below the cosmological signal, while those from reduced shear approximation contribute $\lesssim 1\%$ to the signal. Source clustering effects and baryon imprints contribute $1-10\%$. Enforcing scale cuts to limit systematics-driven biases in parameter constraints degrades these constraints a noticeable amount, and this degradation is similar for the CDFs and the moments. We also detect correlations between the observed convergence field and the shape noise field at $13σ$. We find that the non-Gaussian correlations in the noise field must be modeled accurately to use the CDFs, or other statistics sensitive to all moments, as a rigorous cosmology tool. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 21 pages, 12 figures

Detection of the significant impact of source clustering on higher-order statistics with DES Year 3 weak gravitational lensing data

Authors: M. Gatti, N. Jeffrey, L. Whiteway, V. Ajani, T. Kacprzak, D. Zürcher, C. Chang, B. Jain, J. Blazek, E. Krause, A. Alarcon, A. Amon, K. Bechtol, M. Becker, G. Bernstein, A. Campos, R. Chen, A. Choi, C. Davis, J. Derose, H. T. Diehl, S. Dodelson, C. Doux, K. Eckert, J. Elvin-Poole , et al. (76 additional authors not shown)

Abstract: We demonstrate and measure the impact of source galaxy clustering on higher-order summary statistics of weak gravitational lensing data. By comparing simulated data with galaxies that either trace or do not trace the underlying density field, we show this effect can exceed measurement uncertainties for common higher-order statistics for certain analysis choices. Source clustering effects are large… ▽ More We demonstrate and measure the impact of source galaxy clustering on higher-order summary statistics of weak gravitational lensing data. By comparing simulated data with galaxies that either trace or do not trace the underlying density field, we show this effect can exceed measurement uncertainties for common higher-order statistics for certain analysis choices. Source clustering effects are larger at small scales and for statistics applied to combinations of low and high redshift samples, and diminish at high redshift. We evaluate the impact on different weak lensing observables, finding that third moments and wavelet phase harmonics are more affected than peak count statistics. Using Dark Energy Survey Year 3 data we construct null tests for the source-clustering-free case, finding a $p$-value of $p=4\times10^{-3}$ (2.6 $σ$) using third-order map moments and $p=3\times10^{-11}$ (6.5 $σ$) using wavelet phase harmonics. The impact of source clustering on cosmological inference can be either be included in the model or minimized through \textit{ad-hoc} procedures (e.g. scale cuts). We verify that the procedures adopted in existing DES Y3 cosmological analyses (using map moments and peaks) were sufficient to render this effect negligible. Failing to account for source clustering can significantly impact cosmological inference from higher-order gravitational lensing statistics, e.g. higher-order N-point functions, wavelet-moment observables (including phase harmonics and scattering transforms), and deep learning or field level summary statistics of weak lensing maps. We provide recipes both to minimise the impact of source clustering and to incorporate source clustering effects into forward-modelled mock data. △ Less

Submitted 27 July, 2023; v1 submitted 25 July, 2023; originally announced July 2023.

Comments: 5 pages, 2 figures, submitted to MNRAS Letters

DES Y3 + KiDS-1000: Consistent cosmology combining cosmic shear surveys

Authors: Dark Energy Survey, Kilo-Degree Survey Collaboration, :, T. M. C. Abbott, M. Aguena, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, M. Asgari, S. Avila, D. Bacon, K. Bechtol, M. R. Becker, G. M. Bernstein, E. Bertin, M. Bilicki, J. Blazek, S. Bocquet, D. Brooks, P. Burger, D. L. Burke, H. Camacho, A. Campos, A. Carnero Rosell , et al. (138 additional authors not shown)

Abstract: We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of… ▽ More We present a joint cosmic shear analysis of the Dark Energy Survey (DES Y3) and the Kilo-Degree Survey (KiDS-1000) in a collaborative effort between the two survey teams. We find consistent cosmological parameter constraints between DES Y3 and KiDS-1000 which, when combined in a joint-survey analysis, constrain the parameter $S_8 = σ_8 \sqrt{Ω_{\rm m}/0.3}$ with a mean value of $0.790^{+0.018}_{-0.014}$. The mean marginal is lower than the maximum a posteriori estimate, $S_8=0.801$, owing to skewness in the marginal distribution and projection effects in the multi-dimensional parameter space. Our results are consistent with $S_8$ constraints from observations of the cosmic microwave background by Planck, with agreement at the $1.7σ$ level. We use a Hybrid analysis pipeline, defined from a mock survey study quantifying the impact of the different analysis choices originally adopted by each survey team. We review intrinsic alignment models, baryon feedback mitigation strategies, priors, samplers and models of the non-linear matter power spectrum. △ Less

Submitted 19 October, 2023; v1 submitted 26 May, 2023; originally announced May 2023.

Comments: 40 pages, 21 figures, 15 tables, accepted Open Journal of Astrophysics. Download the chains from https://des.ncsa.illinois.edu/releases/y3a2/Y3key-joint-des-kids or create your own chains with CosmoSIS using https://github.com/joezuntz/cosmosis-standard-library/blob/main/examples/des-y3_and_kids-1000.ini Watch the core team discuss this analysis at https://cosmologytalks.com/2023/05/26/des-kids

Report number: FERMILAB-PUB-23-267-PPD

Spectral age distribution for radio-loud active galaxies in the XMM-LSS field

Authors: Siddhant Pinjarkar, Martin J. Hardcastle, Jeremy J. Harwood, Dharam V. Lal, Peter W. Hatfield, Matt J. Jarvis, Zara Randriamanakoto, Imogen H. Whittam

Abstract: Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require hi… ▽ More Jets of energetic particles, as seen in FR type-I and FR type-II sources, ejected from the center of Radio-Loud AGN affect the sources surrounding intracluster medium/intergalactic medium. Placing constraints on the age of such sources is important in order to measure the jet powers and determine the effects on feedback. To evaluate the age of these sources using spectral age models, we require high-resolution multi-wavelength data. The new sensitive and high-resolution MIGHTEE survey of the XMM-LSS field along with data from the Low Frequency Array (LOFAR) and the Giant Metrewave Radio Telescope (GMRT) provide data taken at different frequencies with similar resolution, which enables us to determine the spectral age distribution for radio loud AGN in the survey field. In this study we present a sample of 28 radio galaxies with their best fitting spectral age distribution analyzed using the Jaffe-Perola (JP) model on a pixel-by-pixel basis. Fits are generally good and objects in our sample show maximum ages within the range of 2.8 Myr to 115 Myr with a median of 8.71 Myr. High-resolution maps over a range of frequencies are required to observe detailed age distributions for small sources and high-sensitivity maps will be needed in order to observe fainter extended emission. We do not observe any correlation between the total physical size of the sources and their age and we speculate both dynamical models and the approach to spectral age analysis may need some modification to account for our observations. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 20 pages, 9 figures

The Spitzer Extragalactic Representative Volume Survey and DeepDrill extension: clustering of near-infrared galaxies

Authors: Eelco van Kampen, Mark Lacy, Duncan Farrah, Claudia del P. Lagos, Matt Jarvis, Claudia Maraston, Kristina Nyland, Seb Oliver, Jason Surace, Jessica Thorne

Abstract: We have measured the angular auto-correlation function of near-infrared galaxies in SERVS+DeepDrill, the Spitzer Extragalactic Representative Volume Survey and its follow-up survey of the Deep Drilling Fields, in three large fields totalling over 20 sq. deg on the sky, observed in two bands centred on 3.6 and 4.5 micron. We performed this analysis on the full sample as well as on sources selected… ▽ More We have measured the angular auto-correlation function of near-infrared galaxies in SERVS+DeepDrill, the Spitzer Extragalactic Representative Volume Survey and its follow-up survey of the Deep Drilling Fields, in three large fields totalling over 20 sq. deg on the sky, observed in two bands centred on 3.6 and 4.5 micron. We performed this analysis on the full sample as well as on sources selected by [3.6]-[4.5] colour in order to probe clustering for different redshift regimes. We estimated the spatial correlation strength as well, using the redshift distribution from S-COSMOS with the same source selection. The strongest clustering was found for our bluest subsample, with z~0.7, which has the narrowest redshift distribution of all our subsamples. We compare these estimates to previous results from the literature, but also to estimates derived from mock samples, selected in the same way as the observational data, using deep light-cones generated from the SHARK semi-analytical model of galaxy formation. For all simulated (sub)samples we find a slightly steeper slope than for the corresponding observed ones, but the spatial clustering length is comparable in most cases. △ Less

Submitted 12 May, 2023; originally announced May 2023.

Comments: 20 pages, 15 figures, 6 tables, accepted for publication in MNRAS

The LOFAR Two-metre Sky Survey Deep Fields Data Release 1: V. Survey description, source classifications and host galaxy properties

Authors: P. N. Best, R. Kondapally, W. L. Williams, R. K. Cochrane, K. J. Duncan, C. L. Hale, P. Haskell, K. Malek, I. McCheyne, D. J. B. Smith, L. Wang, A. Botteon, M. Bonato, M. Bondi, G. Calistro Rivera, F. Gao, G. Gurkan, M. J. Hardcastle, M. J. Jarvis, B. Mingo, H. Miraghaei, L. K. Morabito, D. Nisbet, I. Prandoni, H. J. A. Rottgering , et al. (4 additional authors not shown)

Abstract: Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried… ▽ More Source classifications, stellar masses and star formation rates are presented for 80,000 radio sources from the first data release of the Low Frequency Array Two-metre Sky Survey (LoTSS) Deep Fields, which represents the widest deep radio survey ever undertaken. Using deep multi-wavelength data spanning from the ultraviolet to the far-infrared, spectral energy distribution (SED) fitting is carried out for all of the LoTSS-Deep host galaxies using four different SED codes, two of which include modelling of the contributions from an active galactic nucleus (AGN). Comparing the results of the four codes, galaxies that host a radiative AGN are identified, and an optimised consensus estimate of the stellar mass and star-formation rate for each galaxy is derived. Those galaxies with an excess of radio emission over that expected from star formation are then identified, and the LoTSS-Deep sources are divided into four classes: star-forming galaxies, radio-quiet AGN, and radio-loud high-excitation and low-excitation AGN. Ninety-five per cent of the sources can be reliably classified, of which more than two-thirds are star-forming galaxies, ranging from normal galaxies in the nearby Universe to highly-starbursting systems at z>4. Star-forming galaxies become the dominant population below 150-MHz flux densities of about 1 mJy, accounting for 90 per cent of sources at a 150-MHz flux density of 100 microJy. Radio-quiet AGN comprise around 10 per cent of the overall population. Results are compared against the predictions of the SKADS and T-RECS radio sky simulations, and improvements to the simulations are suggested. △ Less

Submitted 9 May, 2023; originally announced May 2023.

Comments: Accepted for publication in MNRAS. Catalogues available at www.lofar-surveys.org/deepfields.html. 27 pages

MIGHTEE-HI: The first MeerKAT HI mass function from an untargeted interferometric survey

Authors: Anastasia A. Ponomareva, Matt J. Jarvis, Hengxing Pan, Natasha Maddox, Michael G. Jones, Bradley S. Frank, Sambatriniaina H. A. Rajohnson, Wanga Mulaudzi, Martin Meyer, Elizabeth A. K. Adams, Maarten Baes, Kelley M. Hess, Sushma Kurapati, Isabella Prandoni, Francesco Sinigaglia, Kristine Spekkens, Madalina Tudorache, Ian Heywood, Jordan D. Collier, Srikrishna Sekhar

Abstract: We present the first measurement of the HI mass function (HIMF) using data from MeerKAT, based on 276 direct detections from the MIGHTEE Survey Early Science data covering a period of approximately a billion years ($0 \leq z \leq 0.084 $). This is the first HIMF measured using interferometric data over non-group or cluster field, i.e. a deep blank field. We constrain the parameters of the Schechte… ▽ More We present the first measurement of the HI mass function (HIMF) using data from MeerKAT, based on 276 direct detections from the MIGHTEE Survey Early Science data covering a period of approximately a billion years ($0 \leq z \leq 0.084 $). This is the first HIMF measured using interferometric data over non-group or cluster field, i.e. a deep blank field. We constrain the parameters of the Schechter function which describes the HIMF with two different methods: $1/\rm V_{\rm max}$ and Modified Maximum Likelihood (MML). We find a low-mass slope $α=-1.29^{+0.37}_{-0.26}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.07^{+0.24}_{-0.24}$ and normalisation $\log_{10}(φ_{*}/\rm Mpc^{-3})=-2.34^{+0.32}_{-0.36}$ ($H_0 = 67.4$ kms$^{-1}$ Mpc$^{-1}$) for $1/\rm V_{\rm max}$ and $α=-1.44^{+0.13}_{-0.10}$, `knee' mass $\log_{10}(M_{*}/{\rm M_{\odot}}) = 10.22^{+0.10}_{-0.13}$ and normalisation $\log_{10}(φ_{*}/\rm Mpc^{-3})=-2.52^{+0.19}_{-0.14}$ for MML. When using $1/\rm V_{\rm max}$ we find both the low-mass slope and `knee' mass to be consistent within $1σ$ with previous studies based on single-dish surveys. The cosmological mass density of HI is found to be slightly larger than previously reported: $Ω_{\rm HI}=5.46^{+0.94}_{-0.99} \times 10^{-4}h^{-1}_{67.4}$ from $1/\rm V_{\rm max}$ and $Ω_{\rm HI}=6.31^{+0.31}_{-0.31} \times 10^{-4}h^{-1}_{67.4}$ from MML but consistent within the uncertainties. We find no evidence for evolution of the HIMF over the last billion years. △ Less

Submitted 25 April, 2023; originally announced April 2023.

Comments: 13 pages, 9 figures, accepted for publication in MNRAS

The bright end of the galaxy luminosity function at $z \simeq 7$ from the VISTA VIDEO survey

Authors: R. G. Varadaraj, R. A. A. Bowler, M. J. Jarvis, N. J. Adams, B. Häußler

Abstract: We have conducted a search for $z\simeq7$ Lyman break galaxies over 8.2 square degrees of near-infrared imaging from the VISTA Deep Extragalactic Observations (VIDEO) survey in the XMM-Newton - Large Scale Structure (XMM-LSS) and the Extended Chandra Deep Field South (ECDF-S) fields. Candidate galaxies were selected from a full photometric redshift analysis down to a $Y+J$ depth of 25.3 ($5σ$), ut… ▽ More We have conducted a search for $z\simeq7$ Lyman break galaxies over 8.2 square degrees of near-infrared imaging from the VISTA Deep Extragalactic Observations (VIDEO) survey in the XMM-Newton - Large Scale Structure (XMM-LSS) and the Extended Chandra Deep Field South (ECDF-S) fields. Candidate galaxies were selected from a full photometric redshift analysis down to a $Y+J$ depth of 25.3 ($5σ$), utilizing deep auxiliary optical and Spitzer/IRAC data to remove brown dwarf and red interloper galaxy contaminants. Our final sample consists of 28 candidate galaxies at $6.5\le z \le7.5$ with $-23.5 \le M_{\mathrm{UV}} \le -21.6$. We derive stellar masses of $9.1 \le \mathrm{log}_{10}(M/M_{\odot}) \le 10.9$ for the sample, suggesting that these candidates represent some of the most massive galaxies known at this epoch. We measure the rest-frame UV luminosity function (LF) at $z\simeq7$, confirming previous findings of a gradual decline in number density at the bright-end ($M_{\mathrm{UV}} < -22$) that is well described by a double-power law (DPL). We show that quasar contamination in this magnitude range is expected to be minimal, in contrast to conclusions from recent pure-parallel Hubble studies. Our results are up to a factor of ten lower than previous determinations from optical-only ground-based studies at $M_{\rm UV} \lesssim - 23$. We find that the inclusion of $YJHK_{s}$ photometry is vital for removing brown-dwarf contaminants, and $z \simeq 7$ samples based on red-optical data alone could be highly contaminated ($\gtrsim 50$ per cent). In comparison with other robust $z > 5$ samples, our results further support little evolution in the very bright-end of the rest-frame UV LF from $z = 5-10$, potentially signalling a lack of mass quenching and/or dust obscuration in the most massive galaxies in the first Gyr. △ Less

Submitted 10 July, 2023; v1 submitted 5 April, 2023; originally announced April 2023.

Comments: 15 pages, 6 figures, 5 tables (plus additional figures/tables in Appendix). Accepted for publication in MNRAS

The Dark Energy Survey Supernova Program: Corrections on photometry due to wavelength-dependent atmospheric effects

Authors: J. Lee, M. Acevedo, M. Sako, M. Vincenzi, D. Brout, B. Sanchez, R. Chen, T. M. Davis, M. Jarvis, D. Scolnic, H. Qu, L. Galbany, R. Kessler, J. Lasker, M. Sullivan, P. Wiseman, M. Aguena, S. Allam, O. Alves, F. Andrade-Oliveira, E. Bertin, S. Bocquet, D. Brooks, D. L. Burke, A. Carnero Rosell , et al. (42 additional authors not shown)

Abstract: Wavelength-dependent atmospheric effects impact photometric supernova flux measurements for ground-based observations. We present corrections on supernova flux measurements from the Dark Energy Survey Supernova Program's 5YR sample (DES-SN5YR) for differential chromatic refraction (DCR) and wavelength-dependent seeing, and we show their impact on the cosmological parameters $w$ and $Ω_m$. We use… ▽ More Wavelength-dependent atmospheric effects impact photometric supernova flux measurements for ground-based observations. We present corrections on supernova flux measurements from the Dark Energy Survey Supernova Program's 5YR sample (DES-SN5YR) for differential chromatic refraction (DCR) and wavelength-dependent seeing, and we show their impact on the cosmological parameters $w$ and $Ω_m$. We use $g-i$ colors of Type Ia supernovae (SNe Ia) to quantify astrometric offsets caused by DCR and simulate point spread functions (PSFs) using the GalSIM package to predict the shapes of the PSFs with DCR and wavelength-dependent seeing. We calculate the magnitude corrections and apply them to the magnitudes computed by the DES-SN5YR photometric pipeline. We find that for the DES-SN5YR analysis, not accounting for the astrometric offsets and changes in the PSF shape cause an average bias of $+0.2$ mmag and $-0.3$ mmag respectively, with standard deviations of $0.7$ mmag and $2.7$ mmag across all DES observing bands (\textit{griz}) throughout all redshifts. When the DCR and seeing effects are not accounted for, we find that $w$ and $Ω_m$ are lower by less than $0.004\pm0.02$ and $0.001\pm0.01$ respectively, with $0.02$ and $0.01$ being the $1σ$ statistical uncertainties. Although we find that these biases do not limit the constraints of the DES-SN5YR sample, future surveys with much higher statistics, lower systematics, and especially those that observe in the $u$ band will require these corrections as wavelength-dependent atmospheric effects are larger at shorter wavelengths. We also discuss limitations of our method and how they can be better accounted for in future surveys. △ Less

Submitted 4 April, 2023; originally announced April 2023.

Comments: 15 pages, 13 figures, accepted by the Astronomical Journal

Report number: DES-2022-0740, FERMILAB-PUB-23-040-PPD

MIGHTEE-\HI: Possible interactions with the galaxy NGC~895

Authors: Brenda Namumba, Javier Román, Jesus Falcon Barroso, Johan H. Knapen, Ianjamasimanana Roger, Elizabeth Naluminsa, Gyula I. G. Jozsa, Marie Korsaga, Natasha Maddox, Brad Frank, Sinenhlanhla Sikhosana, Samuel Legodi, Claude Carignan, Anastasia A. Ponomareva, Tom Jarrett, Danielle Lucero, Oleg M. Smirnov, Thijs van der Hulst, D. J. Pisano, kasia Malek, Lucia Marchetti, Mattia Vaccari, Matt Jarvis, Maarten Baes, Martin Meyer , et al. (7 additional authors not shown)

Abstract: The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction fea… ▽ More The transformation and evolution of a galaxy is strongly influenced by interactions with its environment. Neutral hydrogen (\HI) is an excellent way to trace these interactions. Here, we present \HI\ observations of the spiral galaxy NGC~895, which was previously thought to be isolated. High-sensitivity \HI\ observations from the MeerKAT large survey project MIGHTEE reveal possible interaction features, such as extended spiral arms, and the two newly discovered \HI\ companions, that drive us to change the narrative that it is an isolated galaxy. We combine these observations with deep optical images from the Hyper Suprime Camera to show an absence of tidal debris between NGC 895 and its companions. We do find an excess of light in the outer parts of the companion galaxy MGTH$\_$J022138.1-052631 which could be an indication of external perturbation and thus possible sign of interactions. Our analysis shows that NGC~895 is an actively star-forming galaxy with a SFR of $\mathrm{1.75 \pm 0.09 [M_{\odot}/yr]}$, a value typical for high stellar mass galaxies on the star forming main sequence. It is reasonable to state that different mechanisms may have contributed to the observed features in NGC~895 and this emphasizes the need to revisit the target with more detailed observations. Our work shows the high potential and synergy of using state-of-the-art data in both \HI\ and optical to reveal a more complete picture of galaxy environments. △ Less

Submitted 16 March, 2023; originally announced March 2023.

Comments: 14 pages, 10 figures. Accepted for publication in MNRAS

Simulating image coaddition with the Nancy Grace Roman Space Telescope: II. Analysis of the simulated images and implications for weak lensing

Authors: Masaya Yamamoto, Katherine Laliotis, Emily Macbeth, Tianqing Zhang, Christopher M. Hirata, M. A. Troxel, Kaili Cao, Ami Choi, Jahmour Givans, Katrin Heitmann, Mustapha Ishak, Mike Jarvis, Eve Kovacs, Heyang Long, Rachel Mandelbaum, Andy Park, Anna Porredon, Christopher W. Walter, W. Michael Wood-Vasey

Abstract: One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the ou… ▽ More One challenge for applying current weak lensing analysis tools to the Nancy Grace Roman Space Telescope is that individual images will be undersampled. Our companion paper presented an initial application of Imcom - an algorithm that builds an optimal mapping from input to output pixels to reconstruct a fully sampled combined image - on the Roman image simulations. In this paper, we measure the output noise power spectra, identify the sources of the major features in the power spectra, and show that simple analytic models that ignore sampling effects underestimate the power spectra of the coadded noise images. We compute the moments of both idealized injected stars and fully simulated stars in the coadded images, and their 1- and 2-point statistics. We show that the idealized injected stars have root-mean-square ellipticity errors (1 - 6) x 10-4 per component depending on the band; the correlation functions are >= 2 orders of magnitude below requirements, indicating that the image combination step itself is using a small fraction of the overall Roman 2nd moment error budget, although the 4th moments are larger and warrant further investigation. The stars in the simulated sky images, which include blending and chromaticity effects, have correlation functions near the requirement level (and below the requirement level in a wide-band image constructed by stacking all 4 filters). We evaluate the noise-induced biases in the ellipticities of injected stars, and explain the resulting trends with an analytical model. We conclude by enumerating the next steps in developing an image coaddition pipeline for Roman. △ Less

Submitted 12 January, 2024; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 25 pages, 20 figures. Submitted to MNRAS

Simulating image coaddition with the Nancy Grace Roman Space Telescope: I. Simulation methodology and general results

Authors: Christopher M. Hirata, Masaya Yamamoto, Katherine Laliotis, Emily Macbeth, M. A. Troxel, Tianqing Zhang, Kaili Cao, Ami Choi, Jahmour Givans, Katrin Heitmann, Mustapha Ishak, Mike Jarvis, Eve Kovacs, Heyang Long, Rachel Mandelbaum, Andy Park, Anna Porredon, Christopher W. Walter, W. Michael Wood-Vasey

Abstract: The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that… ▽ More The upcoming Nancy Grace Roman Space Telescope will carry out a wide-area survey in the near infrared. A key science objective is the measurement of cosmic structure via weak gravitational lensing. Roman data will be undersampled, which introduces new challenges in the measurement of source galaxy shapes; a potential solution is to use linear algebra-based coaddition techniques such as Imcom that combine multiple undersampled images to produce a single oversampled output mosaic with a desired "target" point spread function (PSF). We present here an initial application of Imcom to 0.64 square degrees of simulated Roman data, based on the Roman branch of the Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2) simulation. We show that Imcom runs successfully on simulated data that includes features such as plate scale distortions, chip gaps, detector defects, and cosmic ray masks. We simultaneously propagate grids of injected sources and simulated noise fields as well as the full simulation. We quantify the residual deviations of the PSF from the target (the "leakage"), as well as noise properties of the output images; we discuss how the overall tiling pattern as well as Moiré patterns appear in the final leakage and noise maps. We include appendices on interpolation algorithms and the interaction of undersampling with image processing operations that may be of broader applicability. The companion paper ("Paper II") explores the implications for weak lensing analyses. △ Less

Submitted 12 January, 2024; v1 submitted 15 March, 2023; originally announced March 2023.

Comments: 28 pages, 19 figures, matches version accepted by MNRAS

Metadetection Weak Lensing for the Vera C. Rubin Observatory

Authors: Erin S. Sheldon, Matthew R. Becker, Michael Jarvis, Robert Armstrong, The LSST Dark Energy Science Collaboration

Abstract: Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Scien… ▽ More Forthcoming astronomical imaging surveys will use weak gravitational lensing shear as a primary probe to study dark energy, with accuracy requirements at the 0.1% level. We present an implementation of the Metadetection shear measurement algorithm for use with the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). This new code works with the data products produced by the LSST Science Pipelines, and uses the pipeline algorithms when possible. We tested the code using a new set of simulations designed to mimic LSST imaging data. The simulated images contained semi-realistic galaxies, stars with representative distributions of magnitudes and galactic spatial density, cosmic rays, bad CCD columns and spatially variable point spread functions. Bright stars were saturated and simulated ``bleed trails'' were drawn. Problem areas were interpolated, and the images were coadded into small cells, excluding images not fully covering the cell to guarantee a continuous point spread function. In all our tests the measured shear was accurate within the LSST requirements. △ Less

Submitted 11 May, 2023; v1 submitted 6 March, 2023; originally announced March 2023.

Comments: 10 pages, 5 figures; minor text updates and additional references

The Intrinsic Alignment of Red Galaxies in DES Y1 redMaPPer Galaxy Clusters

Authors: C. Zhou, A. Tong, M. A. Troxel, J. Blazek, C. Lin, D. Bacon, L. Bleem, A. Carnero Rosell, C. Chang, M. Costanzi, J. DeRose, J. P. Dietrich, A. Drlica-Wagner, D. Gruen, R. A. Gruendl, B. Hoyle, M. Jarvis, N. MacCrann, B. Mawdsley, T. McClintock, P. Melchior, J. Prat, A. Pujol, E. Rozo, E. S. Rykoff , et al. (57 additional authors not shown)

Abstract: Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We meas… ▽ More Clusters of galaxies are sensitive to the most nonlinear peaks in the cosmic density field. The weak gravitational lensing of background galaxies by clusters can allow us to infer their masses. However, galaxies associated with the local environment of the cluster can also be intrinsically aligned due to the local tidal gradient, contaminating any cosmology derived from the lensing signal. We measure this intrinsic alignment in Dark Energy Survey (DES) Year 1 redMaPPer clusters. We find evidence of a non-zero mean radial alignment of galaxies within clusters between redshift 0.1-0.7. We find a significant systematic in the measured ellipticities of cluster satellite galaxies that we attribute to the central galaxy flux and other intracluster light. We attempt to correct this signal, and fit a simple model for intrinsic alignment amplitude ($A_{\textrm{IA}}$) to the measurement, finding $A_{\textrm{IA}}=0.15\pm 0.04$, when excluding data near the edge of the cluster. We find a significantly stronger alignment of the central galaxy with the cluster dark matter halo at low redshift and with higher richness and central galaxy absolute magnitude (proxies for cluster mass). This is an important demonstration of the ability of large photometric data sets like DES to provide direct constraints on the intrinsic alignment of galaxies within clusters. These measurements can inform improvements to small-scale modeling and simulation of the intrinsic alignment of galaxies to help improve the separation of the intrinsic alignment signal in weak lensing studies. △ Less

Submitted 5 September, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 14 pages, 13 figures. Accepted to MNRAS

The Dark Energy Survey Year 3 and eBOSS: constraining galaxy intrinsic alignments across luminosity and colour space

Authors: S. Samuroff, R. Mandelbaum, J. Blazek, A. Campos, N. MacCrann, G. Zacharegkas, A. Amon, J. Prat, S. Singh, J. Elvin-Poole, A. J. Ross, A. Alarcon, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, C. Chang, R. Chen, A. Choi, M. Crocce, C. Davis, J. DeRose , et al. (82 additional authors not shown)

Abstract: We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS a… ▽ More We present direct constraints on galaxy intrinsic alignments using the Dark Energy Survey Year 3 (DES Y3), the Extended Baryon Oscillation Spectroscopic Survey (eBOSS) and its precursor, the Baryon Oscillation Spectroscopic Survey (BOSS). Our measurements incorporate photometric red sequence (redMaGiC) galaxies from DES with median redshift $z\sim0.2-1.0$, luminous red galaxies (LRGs) from eBOSS at $z\sim0.8$, and also a SDSS-III BOSS CMASS sample at $z\sim0.5$. We measure two point intrinsic alignment correlations, which we fit using a model that includes lensing, magnification and photometric redshift error. Fitting on scales $6<r_{\rm p} < 70$ Mpc$/h$, we make a detection of intrinsic alignments in each sample, at $5σ-22σ$ (assuming a simple one parameter model for IAs). Using these red samples, we measure the IA-luminosity relation. Our results are statistically consistent with previous results, but offer a significant improvement in constraining power, particularly at low luminosity. With this improved precision, we see detectable dependence on colour between broadly defined red samples. It is likely that a more sophisticated approach than a binary red/blue split, which jointly considers colour and luminosity dependence in the IA signal, will be needed in future. We also compare the various signal components at the best fitting point in parameter space for each sample, and find that magnification and lensing contribute $\sim2-18\%$ of the total signal. As precision continues to improve, it will certainly be necessary to account for these effects in future direct IA measurements. Finally, we make equivalent measurements on a sample of Emission Line Galaxies (ELGs) from eBOSS at $z\sim 0.8$. We report a null detection, constraining the IA amplitude (assuming the nonlinear alignment model) to be $A_1=0.07^{+0.32}_{-0.42}$ ($|A_1|<0.78$ at $95\%$ CL). △ Less

Submitted 21 December, 2022; originally announced December 2022.

Comments: Submitted to MNRAS; 22 pages, 11 figures

Report number: FERMILAB-PUB-22-914-PPD

The catalog-to-cosmology framework for weak lensing and galaxy clustering for LSST

Authors: J. Prat, J. Zuntz, Y. Omori, C. Chang, T. Tröster, E. Pedersen, C. García-García, E. Phillips-Longley, J. Sanchez, D. Alonso, X. Fang, E. Gawiser, K. Heitmann, M. Ishak, M. Jarvis, E. Kovacs, P. Larsen, Y. -Y. Mao, L. Medina Varela, M. Paterno, S. D. Vitenti, Z. Zhang, The LSST Dark Energy Science Collaboration

Abstract: We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Sci… ▽ More We present TXPipe, a modular, automated and reproducible pipeline for ingesting catalog data and performing all the calculations required to obtain quality-assured two-point measurements of lensing and clustering, and their covariances, with the metadata necessary for parameter estimation. The pipeline is developed within the Rubin Observatory Legacy Survey of Space and Time (LSST) Dark Energy Science Collaboration (DESC), and designed for cosmology analyses using LSST data. In this paper, we present the pipeline for the so-called 3x2pt analysis -- a combination of three two-point functions that measure the auto- and cross-correlation between galaxy density and shapes. We perform the analysis both in real and harmonic space using TXPipe and other LSST-DESC tools. We validate the pipeline using Gaussian simulations and show that it accurately measures data vectors and recovers the input cosmology to the accuracy level required for the first year of LSST data under this simplified scenario. We also apply the pipeline to a realistic mock galaxy sample extracted from the CosmoDC2 simulation suite (Korytov et al. 2019). TXPipe establishes a baseline framework that can be built upon as the LSST survey proceeds. Furthermore, the pipeline is designed to be easily extended to science probes beyond the 3x2pt analysis. △ Less

Submitted 21 April, 2023; v1 submitted 19 December, 2022; originally announced December 2022.

Comments: 20+11 pages, 10+11 figures. Version accepted in The Open Journal of Astrophysics

The wide-field, multiplexed, spectroscopic facility WEAVE: Survey design, overview, and simulated implementation

Authors: Shoko Jin, Scott C. Trager, Gavin B. Dalton, J. Alfonso L. Aguerri, J. E. Drew, Jesús Falcón-Barroso, Boris T. Gänsicke, Vanessa Hill, Angela Iovino, Matthew M. Pieri, Bianca M. Poggianti, D. J. B. Smith, Antonella Vallenari, Don Carlos Abrams, David S. Aguado, Teresa Antoja, Alfonso Aragón-Salamanca, Yago Ascasibar, Carine Babusiaux, Marc Balcells, R. Barrena, Giuseppina Battaglia, Vasily Belokurov, Thomas Bensby, Piercarlo Bonifacio , et al. (190 additional authors not shown)

Abstract: WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrogr… ▽ More WEAVE, the new wide-field, massively multiplexed spectroscopic survey facility for the William Herschel Telescope, will see first light in late 2022. WEAVE comprises a new 2-degree field-of-view prime-focus corrector system, a nearly 1000-multiplex fibre positioner, 20 individually deployable 'mini' integral field units (IFUs), and a single large IFU. These fibre systems feed a dual-beam spectrograph covering the wavelength range 366$-$959\,nm at $R\sim5000$, or two shorter ranges at $R\sim20\,000$. After summarising the design and implementation of WEAVE and its data systems, we present the organisation, science drivers and design of a five- to seven-year programme of eight individual surveys to: (i) study our Galaxy's origins by completing Gaia's phase-space information, providing metallicities to its limiting magnitude for $\sim$3 million stars and detailed abundances for $\sim1.5$ million brighter field and open-cluster stars; (ii) survey $\sim0.4$ million Galactic-plane OBA stars, young stellar objects and nearby gas to understand the evolution of young stars and their environments; (iii) perform an extensive spectral survey of white dwarfs; (iv) survey $\sim400$ neutral-hydrogen-selected galaxies with the IFUs; (v) study properties and kinematics of stellar populations and ionised gas in $z<0.5$ cluster galaxies; (vi) survey stellar populations and kinematics in $\sim25\,000$ field galaxies at $0.3\lesssim z \lesssim 0.7$; (vii) study the cosmic evolution of accretion and star formation using $>1$ million spectra of LOFAR-selected radio sources; (viii) trace structures using intergalactic/circumgalactic gas at $z>2$. Finally, we describe the WEAVE Operational Rehearsals using the WEAVE Simulator. △ Less

Submitted 31 October, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 41 pages, 27 figures, accepted for publication by MNRAS; updated version including information on individual grants in a revised Acknowledgements section, corrections to the affiliation list, and an updated references list

Non-local contribution from small scales in galaxy-galaxy lensing: Comparison of mitigation schemes

Authors: J. Prat, G. Zacharegkas, Y. Park, N. MacCrann, E. R. Switzer, S. Pandey, C. Chang, J. Blazek, R. Miquel, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, K. Bechtol, M. R. Becker, G. M. Bernstein, R. Chen, A. Choi, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, R. Cawthon, J. Cordero, M. Crocce , et al. (90 additional authors not shown)

Abstract: Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3$\times$2pt, had to discard a lot of signal-to-noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy-galaxy lensing, or the position-shear correlation between lens and source galaxies, is one of the three two-point correlation functi… ▽ More Recent cosmological analyses with large-scale structure and weak lensing measurements, usually referred to as 3$\times$2pt, had to discard a lot of signal-to-noise from small scales due to our inability to accurately model non-linearities and baryonic effects. Galaxy-galaxy lensing, or the position-shear correlation between lens and source galaxies, is one of the three two-point correlation functions that are included in such analyses, usually estimated with the mean tangential shear. However, tangential shear measurements at a given angular scale $θ$ or physical scale $R$ carry information from all scales below that, forcing the scale cuts applied in real data to be significantly larger than the scale at which theoretical uncertainties become problematic. Recently there have been a few independent efforts that aim to mitigate the non-locality of the galaxy-galaxy lensing signal. Here we perform a comparison of the different methods, including the Y-transformation, the Point-Mass marginalization methodology and the Annular Differential Surface Density statistic. We do the comparison at the cosmological constraints level in a combined galaxy clustering and galaxy-galaxy lensing analysis. We find that all the estimators yield equivalent cosmological results assuming a simulated Rubin Observatory Legacy Survey of Space and Time (LSST) Year 1 like setup and also when applied to DES Y3 data. With the LSST Y1 setup, we find that the mitigation schemes yield $\sim$1.3 times more constraining $S_8$ results than applying larger scale cuts without using any mitigation scheme. △ Less

Submitted 4 April, 2023; v1 submitted 7 December, 2022; originally announced December 2022.

Comments: 11+4 pages, 4+4 figures. Matches the accepted version in MNRAS

Report number: FERMILAB-PUB-22-909-PPD

Journal ref: MNRAS, Volume 522, Issue 1, June 2023, Pages 412-425

MIGHTEE: Deep 1.4 GHz Source Counts and the Sky Temperature Contribution of Star Forming Galaxies and Active Galactic Nuclei

Authors: C. L. Hale, I. H. Whittam, M. J. Jarvis, P. N. Best, N. L. Thomas, I. Heywood, M. Prescott, N. Adams, J. Afonso, Fangxia An, R. A. A. Bowler, J. D. Collier, R. H. W. Cook, R. Davé, B. S. Frank, M. Glowacki, P. W. Hatfield, S. Kolwa C. C. Lovell, N. Maddox, L. Marchetti, L. K. Morabito, E. Murphy, I. Prandoni, Z. Randriamanakoto, A. R. Taylor

Abstract: We present deep 1.4 GHz source counts from $\sim$5 deg$^2$ of the continuum Early Science data release of the MeerKAT International Gigahertz Tiered Extragalactic Exploration (MIGHTEE) survey down to $S_{1.4\textrm{GHz}}\sim$15 $μ$Jy. Using observations over two extragalactic fields (COSMOS and XMM-LSS), we provide a comprehensive investigation into correcting the incompleteness of the raw source… ▽ More We present deep 1.4 GHz source counts from $\sim$5 deg$^2$ of the continuum Early Science data release of the MeerKAT International Gigahertz Tiered Extragalactic Exploration (MIGHTEE) survey down to $S_{1.4\textrm{GHz}}\sim$15 $μ$Jy. Using observations over two extragalactic fields (COSMOS and XMM-LSS), we provide a comprehensive investigation into correcting the incompleteness of the raw source counts within the survey to understand the true underlying source count population. We use a variety of simulations that account for: errors in source detection and characterisation, clustering, and variations in the assumed source model used to simulate sources within the field and characterise source count incompleteness. We present these deep source count distributions and use them to investigate the contribution of extragalactic sources to the sky background temperature at 1.4 GHz using a relatively large sky area. We then use the wealth of ancillary data covering{a subset of the COSMOS field to investigate the specific contributions from both active galactic nuclei (AGN) and star forming galaxies (SFGs) to the source counts and sky background temperature. We find, similar to previous deep studies, that we are unable to reconcile the sky temperature observed by the ARCADE 2 experiment. We show that AGN provide the majority contribution to the sky temperature contribution from radio sources, but the relative contribution of SFGs rises sharply below 1 mJy, reaching an approximate 15-25% contribution to the total sky background temperature ($T_b\sim$100 mK) at $\sim$15 $μ$Jy. △ Less

Submitted 10 November, 2022; originally announced November 2022.

Comments: 24 pages, 12 figures; Accepted for publication in MNRAS

MIGHTEE-HI: The HI mass-stellar mass relation over the last billion years

Authors: Hengxing Pan, Matt J. Jarvis, Mario G. Santos, Natasha Maddox, Bradley S. Frank, Anastasia A. Ponomareva, Isabella Prandoni, Sushma Kurapati, Maarten Baes, Pavel E. Mancera Piña, Giulia Rodighiero, Martin J. Meyer, Romeel Davé, Gauri Sharma, Sambatriniaina H. A. Rajohnson, Nathan J. Adams, Rebecca A. A. Bowler, Francesco Sinigaglia, Thijs van der Hulst, Peter W. Hatfield, Srikrishna Sekhar, Jordan D. Collier

Abstract: We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and fin… ▽ More We study the $M_{\rm HI}-M_{\star}$ relation over the last billion years using the MIGHTEE-HI sample. We first model the upper envelope of the $M_{\rm HI}-M_{\star}$ relation with a Bayesian technique applied to a total number of 249 HI-selected galaxies, without binning the datasets, while taking account of the intrinsic scatter. We fit the envelope with both linear and non-linear models, and find that the non-linear model is preferred over the linear one with a measured transition stellar mass of $\log_{10}(M_\star$/$M_{\odot})$ = $9.15\pm0.87$, beyond which the slope flattens. This finding supports the view that the lack of HI gas is ultimately responsible for the decreasing star formation rate observed in the massive main-sequence galaxies. For spirals alone, which are biased towards the massive galaxies in our sample, the slope beyond the transition mass is shallower than for the full sample, indicative of distinct gas processes ongoing for the spirals/high-mass galaxies from other types with lower stellar masses. We then create mock catalogues for the MIGHTEE-HI detections and non-detections with two main galaxy populations of late- and early-type galaxies to measure the underlying $M_{\rm HI}-M_{\star}$ relation. We find that the turnover in this relation persists whether considering the two galaxy populations as a whole or separately. We note that an underlying linear relation could mimic this turnover in the observed scaling relation, but a model with a turnover is strongly preferred. Measurements on the logarithmic average of HI masses against the stellar mass are provided as a benchmark for future studies. △ Less

Submitted 27 July, 2023; v1 submitted 10 October, 2022; originally announced October 2022.

Comments: Revised manuscript accepted for publication in MNRAS

Dark Energy Survey Year 3 results: Magnification modeling and impact on cosmological constraints from galaxy clustering and galaxy-galaxy lensing

Authors: J. Elvin-Poole, N. MacCrann, S. Everett, J. Prat, E. S. Rykoff, J. De Vicente, B. Yanny, K. Herner, A. Ferté, E. Di Valentino, A. Choi, D. L. Burke, I. Sevilla-Noarbe, A. Alarcon, O. Alves, A. Amon, F. Andrade-Oliveira, E. Baxter, K. Bechtol, M. R. Becker, G. M. Bernstein, J. Blazek, H. Camacho, A. Campos, A. Carnero Rosell , et al. (71 additional authors not shown)

Abstract: We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects usin… ▽ More We study the effect of magnification in the Dark Energy Survey Year 3 analysis of galaxy clustering and galaxy-galaxy lensing, using two different lens samples: a sample of Luminous red galaxies, redMaGiC, and a sample with a redshift-dependent magnitude limit, MagLim. We account for the effect of magnification on both the flux and size selection of galaxies, accounting for systematic effects using the Balrog image simulations. We estimate the impact of magnification on the galaxy clustering and galaxy-galaxy lensing cosmology analysis, finding it to be a significant systematic for the MagLim sample. We show cosmological constraints from the galaxy clustering auto-correlation and galaxy-galaxy lensing signal with different magnifications priors, finding broad consistency in cosmological parameters in $Λ$CDM and $w$CDM. However, when magnification bias amplitude is allowed to be free, we find the two-point correlations functions prefer a different amplitude to the fiducial input derived from the image simulations. We validate the magnification analysis by comparing the cross-clustering between lens bins with the prediction from the baseline analysis, which uses only the auto-correlation of the lens bins, indicating systematics other than magnification may be the cause of the discrepancy. We show adding the cross-clustering between lens redshift bins to the fit significantly improves the constraints on lens magnification parameters and allows uninformative priors to be used on magnification coefficients, without any loss of constraining power or prior volume concerns. △ Less

Submitted 26 May, 2023; v1 submitted 20 September, 2022; originally announced September 2022.

Comments: Version accepted for publication in MNRAS. 21 pages, 13 figures, See this https://www.darkenergysurvey.org/des-year-3-cosmology-results-papers/ URL for the full DES Y3 cosmology release

PSFs of coadded images

Authors: Rachel Mandelbaum, Mike Jarvis, Robert H. Lupton, James Bosch, Arun Kannawadi, Michael D. Murphy, Tianqing Zhang, the LSST Dark Energy Science Collaboration

Abstract: We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows: 1. We provide an illustration of how linear coaddition… ▽ More We provide a detailed exploration of the connection between choice of coaddition schemes and the point-spread function (PSF) of the resulting coadded images. In particular, we investigate what properties of the coaddition algorithm lead to the final coadded image having a well-defined PSF. The key elements of this discussion are as follows: 1. We provide an illustration of how linear coaddition schemes can produce a coadd that lacks a well-defined PSF even for relatively simple scenarios and choices of weight functions. 2. We provide a more formal demonstration of the fact that a linear coadd only has a well-defined PSF in the case that either (a) each input image has the same PSF or (b) the coadd is produced with weights that are independent of the signal. 3. We discuss some reasons that two plausible nonlinear coaddition algorithms (median and clipped-mean) fail to produce a consistent PSF profile for stars. 4. We demonstrate that all nonlinear coaddition procedures fail to produce a well-defined PSF for extended objects. In the end, we conclude that, for any purpose where a well-defined PSF is desired, one should use a linear coaddition scheme with weights that do not correlate with the signal and are approximately uniform across typical objects of interest. △ Less

Submitted 2 February, 2023; v1 submitted 19 September, 2022; originally announced September 2022.

Comments: 13 pages, 4 figures; pedagogical article; v2 accepted for publication in the Open Journal of Astrophysics

Journal ref: The Open Journal of Astrophysics, Vol. 6, 2023

A Joint Roman Space Telescope and Rubin Observatory Synthetic Wide-Field Imaging Survey

Authors: M. A. Troxel, C. Lin, A. Park, C. Hirata, R. Mandelbaum, M. Jarvis, A. Choi, J. Givans, M. Higgins, B. Sanchez, M. Yamamoto, H. Awan, J. Chiang, O. Dore, C. W. Walter, T. Zhang, J. Cohen-Tanugi, E. Gawiser, A. Hearin, K. Heitmann, M. Ishak, E. Kovacs, Y. -Y. Mao, M. Wood-Vasey, the LSST Dark Energy Science Collaboration

Abstract: We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging prod… ▽ More We present and validate 20 deg$^2$ of overlapping synthetic imaging surveys representing the full depth of the Nancy Grace Roman Space Telescope High-Latitude Imaging Survey (HLIS) and five years of observations of the Vera C. Rubin Observatory Legacy Survey of Space and Time (LSST). The two synthetic surveys are summarized, with reference to the existing 300 deg$^2$ of LSST simulated imaging produced as part of Dark Energy Science Collaboration (DESC) Data Challenge 2 (DC2). Both synthetic surveys observe the same simulated DESC DC2 universe. For the synthetic Roman survey, we simulate for the first time fully chromatic images along with the detailed physics of the Sensor Chip Assemblies derived from lab measurements using the flight detectors. The simulated imaging and resulting pixel-level measurements of photometric properties of objects span a wavelength range of $\sim$0.3 to 2.0 $μ$m. We also describe updates to the Roman simulation pipeline, changes in how astrophysical objects are simulated relative to the original DC2 simulations, and the resulting simulated Roman data products. We use these simulations to explore the relative fraction of unrecognized blends in LSST images, finding that 20-30% of objects identified in LSST images with $i$-band magnitudes brighter than 25 can be identified as multiple objects in Roman images. These simulations provide a unique testing ground for the development and validation of joint pixel-level analysis techniques of ground- and space-based imaging data sets in the second half of the 2020s -- in particular the case of joint Roman--LSST analyses. △ Less

Submitted 14 September, 2022; originally announced September 2022.

Dark Energy Survey Year 3 Results: Redshift Calibration of the MagLim Lens Sample from the combination of SOMPZ and clustering and its impact on Cosmology

Authors: G. Giannini, A. Alarcon, M. Gatti, A. Porredon, M. Crocce, G. M. Bernstein, R. Cawthon, C. Sánchez, C. Doux, J. Elvin-Poole, M. Raveri, J. Myles, A. Amon, S. Allam, O. Alves, F. Andrade-Oliveira, E. Baxter, K. Bechtol, M. R. Becker, J. Blazek, H. Camacho, A. Campos, A. Carnero Rosell, M. Carrasco Kind, A. Choi , et al. (89 additional authors not shown)

Abstract: We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift ca… ▽ More We present an alternative calibration of the MagLim lens sample redshift distributions from the Dark Energy Survey (DES) first three years of data (Y3). The new calibration is based on a combination of a Self-Organising Maps based scheme and clustering redshifts to estimate redshift distributions and inherent uncertainties, which is expected to be more accurate than the original DES Y3 redshift calibration of the lens sample. We describe in detail the methodology, we validate it on simulations and discuss the main effects dominating our error budget. The new calibration is in fair agreement with the fiducial DES Y3 redshift distributions calibration, with only mild differences ($<3σ$) in the means and widths of the distributions. We study the impact of this new calibration on cosmological constraints, analysing DES Y3 galaxy clustering and galaxy-galaxy lensing measurements, assuming a $Λ$CDM cosmology. We obtain $Ω_{\rm m} = 0.30\pm 0.04$, $σ_8 = 0.81\pm 0.07 $ and $S_8 = 0.81\pm 0.04$, which implies a $\sim 0.4σ$ shift in the $Ω_{\rm}-S_8$ plane compared to the fiducial DES Y3 results, highlighting the importance of the redshift calibration of the lens sample in multi-probe cosmological analyses. △ Less

Submitted 18 October, 2023; v1 submitted 13 September, 2022; originally announced September 2022.

Journal ref: MNRAS, stad2945 (2023)

MIGHTEE-HI: Evolution of HI scaling relations of star-forming galaxies at $z<0.5$

Authors: Francesco Sinigaglia, Giulia Rodighiero, Ed Elson, Mattia Vaccari, Natasha Maddox, Bradley S. Frank, Matt J. Jarvis, Tom Oosterloo, Romeel Davé, Mara Salvato, Maarten Baes, Sabine Bellstedt, Laura Bisigello, Jordan D. Collier, Robin H. W. Cook, Luke J. M. Davies, Jacinta Delhaize, Simon P. Driver, Caroline Foster, Sushma Kurapati, Claudia del P. Lagos, Christopher Lidman, Pavel E. Mancera Piña, Martin J. Meyer, K. Moses Mogotsi , et al. (11 additional authors not shown)

Abstract: We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with… ▽ More We present the first measurements of HI galaxy scaling relations from a blind survey at $z>0.15$. We perform spectral stacking of 9023 spectra of star-forming galaxies undetected in HI at $0.23<z<0.49$, extracted from MIGHTEE-HI Early Science datacubes, acquired with the MeerKAT radio telescope. We stack galaxies in bins of galaxy properties ($M_*$, SFR, and sSFR, with ${\rm sSFR}\equiv M_*/{\rm SFR}$), obtaining $\gtrsim 5σ$ detections in most cases, the strongest HI-stacking detections to date in this redshift range. With these detections, we are able to measure scaling relations in the probed redshift interval, finding evidence for a moderate evolution from the median redshift of our sample $z_{\rm med}\sim 0.37$ to $z\sim 0$. In particular, low-$M_*$ galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 9$) experience a strong HI depletion ($\sim 0.5$ dex in $\log_{10}(M_{\rm HI}/{\rm M}_\odot)$), while massive galaxies ($\log_{10}(M_*/{\rm M_\odot})\sim 11$) keep their HI mass nearly unchanged. When looking at the star formation activity, highly star-forming galaxies evolve significantly in $M_{\rm HI}$ ($f_{\rm HI}$, where $f_{\rm HI}\equiv M_{\rm}/M_*$) at fixed SFR (sSFR), while at the lowest probed SFR (sSFR) the scaling relations show no evolution. These findings suggest a scenario in which low-$M_*$ galaxies have experienced a strong HI depletion during the last $\sim4$ Gyr, while massive galaxies have undergone a significant HI replenishment through some accretion mechanism, possibly minor mergers. Interestingly, our results are in good agreement with the predictions of the SIMBA simulation. We conclude that this work sets novel important observational constraints on galaxy scaling relations. △ Less

Submitted 1 August, 2022; originally announced August 2022.

Comments: 12 pages, 5 figures. Accepted for publication in ApJ Letters on 30 July 2022