Showing 1–50 of 194 results for author: Bilicki, M

Constraints on ultralight bosons from merging binary and remnant black holes observed during the second and third parts of the fourth LIGO-Virgo-KAGRA observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1786 additional authors not shown)

Abstract: We present constraints on ultralight bosons using binary black hole mergers observed in the second and third parts of the fourth LIGO-Virgo-KAGRA observing run. Directed searches are conducted for long-transient gravitational waves from ultralight vector boson clouds around merger remnants, using a hidden-Markov-model (HMM) tracking scheme. We target the remnant black holes formed in the binary co… ▽ More We present constraints on ultralight bosons using binary black hole mergers observed in the second and third parts of the fourth LIGO-Virgo-KAGRA observing run. Directed searches are conducted for long-transient gravitational waves from ultralight vector boson clouds around merger remnants, using a hidden-Markov-model (HMM) tracking scheme. We target the remnant black holes formed in the binary coalescences that produced GW250114 and GW250207. We find no evidence for such signals from either target. Estimating our search sensitivity at a threshold corresponding to a 1% false alarm probability, we thus disfavor vector boson masses in the range of $[2.80, 3.95]\times 10^{-13}$ eV with greater than 90% confidence. In addition, we derive constraints on ultralight scalar and vector bosons from the inferred high spins of the constituent black holes in three binaries, using events GW240515, GW241113, and GW241225_08. The excluded mass ranges in this approach depend on the assumed black-hole ages. At $10^5$ years, corresponding to typical dynamically formed binaries, we exclude scalar and vector bosons in the ranges $[1.39, 6.94]\times 10^{-13}$ eV and $[0.32, 14.4]\times 10^{-13}$ eV at 90% confidence, respectively. △ Less

Submitted 11 August, 2026; originally announced August 2026.

Comments: 27 pages (12 pages author list, 15 pages main paper), 5 figures

Report number: LIGO-P2600218

Advanced Virgo during the LIGO-Virgo-KAGRA fourth observing run

Authors: Virgo Collaboration, F Acernese, A Agapito, D Agarwal, I-L Ahrend, L Aiello, A Ain, W Ali, A Allocca, W Amar, A Amato, F Amicucci, C Amra, M Andia, T Andri, S Antier, F Arciprete, F Armato, N Arnaud, L Asprea, M Assiduo, S Assis de Souza Melo, P Astone, F Attadio, F Aubin , et al. (524 additional authors not shown)

Abstract: From April 10, 2024 to November 18, 2025 Advanced Virgo participated in the fourth observing run of the network of gravitational-wave detectors, together with Advanced LIGO and KAGRA. For this observing run Advanced Virgo has completed its design optical configuration with the installation of a signal recycling mirror. In this paper we describe the challenges encountered in commissioning this opti… ▽ More From April 10, 2024 to November 18, 2025 Advanced Virgo participated in the fourth observing run of the network of gravitational-wave detectors, together with Advanced LIGO and KAGRA. For this observing run Advanced Virgo has completed its design optical configuration with the installation of a signal recycling mirror. In this paper we describe the challenges encountered in commissioning this optical configuration, alongside the other upgrades performed between the third and fourth observing run. The Virgo detector operated with a 68.9% duty cycle and with an angle-averaged median range to binary neutron star mergers of 53 Mpc. △ Less

Submitted 29 July, 2026; originally announced July 2026.

Comments: 54 pages, 21 figures

Report number: VIR-0259F-26

GWTC-5.0: Tests of General Relativity

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The signals from the LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors allow us to perform sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We present the results of seven tests of GR using the observed binary signals in the fifth GW Transient Catalog (GWTC-5.0), i.e., up to and including the second part of the fourth observing run (O4b).… ▽ More The signals from the LIGO-Virgo-KAGRA network of gravitational-wave (GW) detectors allow us to perform sensitive tests of general relativity (GR) in the dynamical and strong-field regime of gravity. We present the results of seven tests of GR using the observed binary signals in the fifth GW Transient Catalog (GWTC-5.0), i.e., up to and including the second part of the fourth observing run (O4b). We restrict our analysis to the confident signals, henceforth called events, observed by at least two detectors that have estimated false alarm rates $\le 10^{-3} \ \rm{yr}^{-1}$. These include 72 events from O4b and five events from the first part of the fourth observing run that are now analyzed due to their increased significance from updated search results, bringing the total number of events for tests of GR in the cumulative GWTC to 168. After subtracting the best-fit waveforms, we find the residuals are consistent with detector noise for all events considered. We also find no strong evidence for additional polarizations beyond those predicted by GR. We perform tests of GW generation, improving the constraints on deviations from the GR post-Newtonian coefficients by factors of 1.2-2.6. Finally, we find overall consistency of the remnants with GR using both time- and frequency-domain methods. For GW240621_195059, postmerger data are consistent with the dominant quadrupolar ($\ell=|m|=2$) mode of a Kerr black hole and its first overtone, with spurious high-frequency content preventing a spectroscopic constraint of GR. In the frequency-domain ringdown analysis, the GR prediction lies in the tails of the combined results, possibly due to the limited catalog size. However, the combined results indicate improved consistency with GR over GWTC-4.0, owing to the contribution of GW250114 with a network matched-filter signal-to-noise ratio of 76.9. Overall, we find no evidence for physics beyond GR. △ Less

Submitted 21 July, 2026; originally announced July 2026.

Comments: 44 pages and 10 figures

Report number: LIGO-P2500781

KiDS-Legacy: The consistency test of the large-scale structure with Bernardeau-Nishimichi-Taruya transform

Authors: Shiming Gu, Ziang Yan, Ludovic van Waerbeke, Francis Bernardeau, Hendrik Hildebrandt, Angus H. Wright, Maciej Bilicki, Christos Georgiou, Shun-Sheng Li, Laila Linke, Lauro Moscardini, Robert Reischke, Benjamin Stölzner

Abstract: We perform the first $k$-cut cosmic shear analysis of the KiDS-Legacy survey. This method uses the Bernardeau-Nishimichi-Taruya (BNT) transform to construct weak-lensing kernels that are more localised than conventional ones, and remove information from selected physical scales while retaining the constraining power of the targeted range. Removing the scale of $k \geq 0.33~\mathrm{Mpc}^{-1}$ from… ▽ More We perform the first $k$-cut cosmic shear analysis of the KiDS-Legacy survey. This method uses the Bernardeau-Nishimichi-Taruya (BNT) transform to construct weak-lensing kernels that are more localised than conventional ones, and remove information from selected physical scales while retaining the constraining power of the targeted range. Removing the scale of $k \geq 0.33~\mathrm{Mpc}^{-1}$ from the KiDS-Legacy pseudo-$C_\ell$ data vector, and using a covariance matrix whose Gaussian component is computed from the theoretical data vector, we find $S_8 = 0.798 \pm 0.045$. This agrees with both the fiducial KiDS-Legacy bandpower result and our no-$k$-cut pseudo-$C_\ell$ posterior to within $0.1σ$, indicating no significant bias from nonlinear astrophysical feedback at the precision of KiDS-Legacy. We also study the case in which the Gaussian covariance is computed from the observed data vector. In this setup, the same scale cut of $k < 0.33~\mathrm{Mpc}^{-1}$ gives a much lower $S_8=0.717_{-0.046}^{+0.047}$. Further $k$-cut tests reveal a mild scale-dependent trend, with larger physical scales preferring lower $S_8$ values and a maximum low- versus high-$k$ deviation of $1.80σ$. Mock tests show that this behaviour is not produced by the covariance prescription or data vector alone, but may arise from their interplay. These results show that BNT $k$-cuts provide both a mitigation strategy for nonlinear systematics and a diagnostic of weak-lensing inference pipelines. △ Less

Submitted 5 July, 2026; originally announced July 2026.

Comments: 19 pages, 6 plots, comments welcome!

KiDS-Legacy: Joint analysis of second- and third-order cosmic shear

Authors: L. Linke, L. Porth, P. Burger, J. Harnois-Déraps, S. Heydenreich, P. Schneider, M. Asgari, M. Bilicki, C. Georgiou, C. Heymans, H. Hildebrandt, H. Hoekstra, P. Jalan, B. Joachimi, S. Joudaki, K. Kuijken, S. Li, L. Moscardini, M. Radovich, R. Reischke, B. Stölzner, A. H. Wright, Z. Yan, Y. -H. Zhang

Abstract: Weak lensing by large-scale structure is a powerful cosmological probe. While most analyses rely on second-order correlations, these are primarily sensitive to the parameter combination $S_8 = σ_8 (Ω_m/0.3)^{0.5}$, limiting their ability to constrain $Ω_m$ and other cosmological parameters independently. Higher-order statistics capture non-Gaussian features of the density field and can therefore b… ▽ More Weak lensing by large-scale structure is a powerful cosmological probe. While most analyses rely on second-order correlations, these are primarily sensitive to the parameter combination $S_8 = σ_8 (Ω_m/0.3)^{0.5}$, limiting their ability to constrain $Ω_m$ and other cosmological parameters independently. Higher-order statistics capture non-Gaussian features of the density field and can therefore break parameter degeneracies and extract more cosmological information from weak lensing surveys. We present a joint analysis of second- and third-order cosmic shear in the final data release of the Kilo-Degree Survey (KiDS-Legacy). We combine COSEBIs (Complete Orthogonal Sets of E-/B-mode Integrals) at scales between 2' and 300' with third-order aperture mass moments at scales between 4' and 32' to perform a joint analysis of second- and third-order statistics. Compared to previous KiDS analyses, we implement several methodological advances: an intrinsic alignment model with redshift and mass dependence, a baryon correction model validated on multiple hydrodynamical simulations, and corrections for reduced shear and source clustering. Combining COSEBIs with third-order aperture mass statistics in KiDS-Legacy yields $Ω_m = 0.297^{+0.056}_{-0.040}$ and $S_8 = 0.806^{+0.025}_{-0.023}$, significantly tightening the $Ω_m$ constraints and more than doubling the figure of merit in the $Ω_m$--$S_8$ plane compared to the two-point analysis alone. The third-order measurements pass stringent internal consistency tests, are fully compatible with the KiDS-Legacy 2-point constraints, other 2+3-point lensing results and with Planck CMB measurements within $1σ$, providing no evidence for an $S_8$ tension and demonstrating the maturity of 3-point cosmic shear as a key probe for forthcoming surveys. △ Less

Submitted 10 June, 2026; originally announced June 2026.

Comments: 14 pages plus appendix, 15 figures, submitted to A&A

GWTC-5.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1788 additional authors not shown)

Abstract: We employ 236 gravitational-wave (GW) sources in the fifth LIGO--Virgo--KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-5.0) to estimate the Hubble constant $H_0$. We compare the luminosity distance measured from GWs to the redshift inferred i) using features in the mass spectrum, and ii) using statistical host galaxy association. Probing the relationship between source lumino… ▽ More We employ 236 gravitational-wave (GW) sources in the fifth LIGO--Virgo--KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-5.0) to estimate the Hubble constant $H_0$. We compare the luminosity distance measured from GWs to the redshift inferred i) using features in the mass spectrum, and ii) using statistical host galaxy association. Probing the relationship between source luminosity distances and redshifts obtained in this way yields constraints on cosmological parameters. We estimate $H_0 = {71.7}_{-7.5}^{+9.4}\,{\text{km}\,\text{s}^{-1}\,\text{Mpc}^{-1}}$ (median with $68\%$ symmetric credible interval). This combines information from the source-frame mass distribution with the $H_0$ measurement from GW170817 and its electromagnetic counterpart as well as galaxy catalog information from Dark Energy Survey Year 6 (DES-Y6). We improve over the GWTC-4.0 measurement by using more GW sources, some with significantly smaller sky localization volumes, which leads to a reduction by $22.0\%$ of the $H_0$ uncertainty and a reconstructed mass distribution with lower uncertainties. We also constrain deviations from general relativity (GR) which affect GW propagation, specifically that modify the luminosity distance inferred from the GW signal. We find no departures from GR in parameterized tests of GW propagation. △ Less

Submitted 4 August, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: main paper: 20 pages and 7 figures; total with appendices: 51 pages and 13 figures. This article draws heavily from the corresponding GWTC-4 article, arXiv:2509.04348

Report number: LIGO-P2600018

GWTC-5.0: Population Properties of Merging Compact Binaries

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1791 additional authors not shown)

Abstract: We present the population properties of merging compact binaries inferred using 267 mergers from the cumulative Gravitational-Wave Transient Catalog 5.0. As this data set contains no new sources with a neutron star, we primarily focus on the properties of the binary black hole mergers. We infer the merger rate of binary black holes with component masses between $2.5\,\mathrm{M}_\odot $ and… ▽ More We present the population properties of merging compact binaries inferred using 267 mergers from the cumulative Gravitational-Wave Transient Catalog 5.0. As this data set contains no new sources with a neutron star, we primarily focus on the properties of the binary black hole mergers. We infer the merger rate of binary black holes with component masses between $2.5\,\mathrm{M}_\odot $ and $200\,\mathrm{M}_\odot $ to be $27.5\text{--} 49.4 \, \mathrm{Gpc}^{-3}\,\mathrm{yr}^{-1}$ (all intervals at $90\%$ credible levels) at redshift $z = 0.2$. We find evidence for a subpopulation of binary black hole mergers that host a rapidly spinning black hole (dimensionless spins $χ\sim 0.7$), consistent with signatures of hierarchical mergers. We find that these occur at two mass scales, the first at primary masses $\sim 10$--$20\,\mathrm{M}_\odot $ and the second above $\sim 45\,\mathrm{M}_\odot $, and we estimate their total rate at $z=0.2$ to be $0.2\text{--} 3.11 \, {\rm Gpc}^{-3} {\rm yr}^{-1}$. We infer that, above $40\,\mathrm{M}_\odot $, the mass distribution of the less massive (secondary) black hole declines more steeply than that of the more massive (primary) one. This is consistent with a flatter mass-ratio distribution and indicates the prevalence of unequal-mass binaries with large primary masses. We find evidence for two features in the black hole mass spectrum: a peak around $10\,\mathrm{M}_\odot $ and a change of slope at around $35\,\mathrm{M}_\odot $. Black holes of $\sim 35\,\mathrm{M}_\odot $ pair preferentially with companions of similar mass. Additionally, we find that the effective inspiral spin distribution of binary black holes is asymmetric about zero, based on which we infer that at least $9 \%$ of mergers occur in channels with some preference for spin-orbit alignment. We find evidence that... △ Less

Submitted 1 July, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: Abstract truncated in Arxiv metadata. The paper appendices draw heavily from the corresponding GWTC-4 article, arXiv:2508.18083

Report number: LIGO-P2600045

GWTC-5.0: Observations from the Second Part of the Fourth LIGO-Virgo-KAGRA Observing Run and Updates to the Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1805 additional authors not shown)

Abstract: Version 5.0 of the Gravitational-Wave Transient Catalog (GWTC-5.0) adds new candidates detected by the LIGO Virgo KAGRA network of observatories through the second part of the fourth observing run (O4b: 2024 April 10 15:00:00 to 2025 January 28 17:00:00 UTC) and four days of the preceding engineering run (2024 April 6 to 2024 April 10). We find 161 compact binary coalescence candidates that are id… ▽ More Version 5.0 of the Gravitational-Wave Transient Catalog (GWTC-5.0) adds new candidates detected by the LIGO Virgo KAGRA network of observatories through the second part of the fourth observing run (O4b: 2024 April 10 15:00:00 to 2025 January 28 17:00:00 UTC) and four days of the preceding engineering run (2024 April 6 to 2024 April 10). We find 161 compact binary coalescence candidates that are identified by at least one of our search algorithms with a probability of astrophysical origin $p_\mathrm{astro} \geq 0.5$ and that are not vetoed during event validation. We also provide detailed source property measurements for 104 candidates that have a false-alarm rate < 1yr$^{-1}$. Based on the inferred component masses, all these candidates are consistent with signals from binary black holes. Median inferred component masses in the new candidates range from 5.14$M_\odot$ (GW241109_115924) to 70$M_\odot$ (GW241116_151753). Improvements in detector sensitivity allow us to observe compact binary coalescences with increasing clarity: 5 binary-black-hole signals have network signal-to-noise ratio exceeding 30, with a maximum to date of 76.9 for GW250114_082203. Such loud signals enable more precise studies of properties of their astrophysical sources and tests of general relativity. We also present updated results up to the first part of the fourth observing run, identifying 229 candidates. This brings the total number of transients in the cumulative GWTC having $p_\mathrm{astro} \geq 0.5$ to 390, further expanding the size of the catalog and our view of the gravitational-wave universe. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: main paper: 30 pages, 8 figures; total with appendices: 43 pages, 9 figures. This article draws heavily from the corresponding GWTC-4 article, arXiv:2508.18082

Report number: LIGO-P2600152

GWTC-5.0: Methods for Identifying and Characterizing Gravitational-wave Transients

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The Gravitational-Wave Transient Catalog (GWTC) is a collection of candidate gravitational-wave transient signals identified and characterized by the LIGO-Virgo-KAGRA Collaboration. Producing the contents of the GWTC from detector data requires complex analysis methods. These comprise techniques to model the signal; identify the transients in the data; evaluate the quality of the data and mitigate… ▽ More The Gravitational-Wave Transient Catalog (GWTC) is a collection of candidate gravitational-wave transient signals identified and characterized by the LIGO-Virgo-KAGRA Collaboration. Producing the contents of the GWTC from detector data requires complex analysis methods. These comprise techniques to model the signal; identify the transients in the data; evaluate the quality of the data and mitigate possible instrumental issues; infer the parameters of each transient; compare the data with the waveform models for compact binary coalescences, and handle the large amount of results associated with all these different analyses. In this paper, we describe the methods employed to produce the catalog's fifth release, GWTC-5.0, focusing on the analysis of the second part of the fourth observing run of LIGO, Virgo and KAGRA. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This article supersedes the corresponding GWTC-4 article, arXiv:2508.18081

Report number: LIGO-P2600166

GWTC-5.0: An Introduction to Version 5.0 of the Gravitational-Wave Transient Catalog

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1800 additional authors not shown)

Abstract: The Gravitational-Wave Transient Catalog (GWTC) is a collection of short-duration (transient) gravitational-wave signals identified by the LIGO-Virgo-KAGRA Collaboration in gravitational-wave data produced by the eponymous detectors. The catalog provides information about the identified candidates, such as the arrival time and amplitude of the signal and properties of the signal's source as inferr… ▽ More The Gravitational-Wave Transient Catalog (GWTC) is a collection of short-duration (transient) gravitational-wave signals identified by the LIGO-Virgo-KAGRA Collaboration in gravitational-wave data produced by the eponymous detectors. The catalog provides information about the identified candidates, such as the arrival time and amplitude of the signal and properties of the signal's source as inferred from the observational data. GWTC is the release of this dataset and version 5.0 extends the catalog to include observations made during the second part of the fourth LIGO-Virgo-KAGRA observing run up until 2025 January 28. This paper marks an introduction to a collection of articles related to this version of the catalog, GWTC-5.0. This update significantly increases the number of detected merging binary systems of black holes and neutron stars to over 300, enabling many follow-up studies toward understanding the gravitational-wave universe. The collection of articles accompanying the catalog provides documentation of the methods used to analyze the data, summaries of the catalog of events, observational measurements drawn from the population, and detailed discussions of selected candidates. △ Less

Submitted 23 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This article supersedes the corresponding GWTC-4 article, arXiv:2508.18080

Report number: LIGO-P2500701

Open Data from LIGO, Virgo, and KAGRA through the Second Part of the Fourth Observing Run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, A. Abe, I. Abouelfettouh, F. Acernese, K. Ackley, A. Adam, S. Adhicary, D. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1787 additional authors not shown)

Abstract: LIGO, Virgo, KAGRA, and GEO 600 form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center (GWOSC). This paper describes open data from this network, including the addition of data from the second part of the fourth observing run (O4b) and selected periods from the preceding engineer… ▽ More LIGO, Virgo, KAGRA, and GEO 600 form a network of gravitational-wave observatories. Data and analysis results from this network are made publicly available through the Gravitational Wave Open Science Center (GWOSC). This paper describes open data from this network, including the addition of data from the second part of the fourth observing run (O4b) and selected periods from the preceding engineering run (ER16), which were collected from times spanning April 6th, 2024 to January 28th, 2025. The public data set includes calibrated strain time series for each instrument, data from additional channels used for noise subtraction and detector characterization, and new analysis data products in the online GWOSC release associated with version 5.0 of the Gravitational-Wave Transient Catalog. △ Less

Submitted 17 June, 2026; v1 submitted 26 May, 2026; originally announced May 2026.

Comments: This version has an updated author list, updated references and some text has been modified to improve clarity. This article draws heavily from the corresponding O4a article, arXiv:2508.18079

Report number: LIGO-P2600085

Caught in the Cosmic Web: Environmental Impacts on the Halo Substructure Boosts to Dark Matter Annihilation Signals

Authors: Feven Markos Hunde, Wojciech A. Hellwing, Maciej Bilicki

Abstract: The annihilation of dark matter (DM) particles is expected to produce Standard Model particles, providing a potential indirect signature of DM. The clumpy substructure of DM haloes amplifies the expected annihilation signal, an effect commonly quantified by the subhalo boost factor. Standard semi-analytic models usually treat this boost as a universal function of host-halo mass, neglecting systema… ▽ More The annihilation of dark matter (DM) particles is expected to produce Standard Model particles, providing a potential indirect signature of DM. The clumpy substructure of DM haloes amplifies the expected annihilation signal, an effect commonly quantified by the subhalo boost factor. Standard semi-analytic models usually treat this boost as a universal function of host-halo mass, neglecting systematic variations induced by the large-scale environment. In this work, we extend this framework by incorporating the influence of the cosmic web on subhalo populations. Using simulation-calibrated, environment-dependent ratios for host-halo concentrations, the subhalo mass function, and internal-structure proxies of subhalos based on the $V_{\max}$--$R_{\max}$ relation, we compute environment-conditioned boost predictions for haloes residing in filaments, walls, and voids. Our main result is the boost factor at fixed host-halo mass, expressed relative to the cosmic-mean prediction, $B(M,\mathrm{env})/B_{\mathrm{CM}}(M)$. We find a clear environmental modulation: in the fiducial distance-dependent model, filament haloes show a mass-dependent transition from a $\sim 15\%$ suppression at the low-mass end to a modest enhancement of $\sim 12\%$ for massive hosts, wall haloes remain intermediate, while void haloes stay suppressed by roughly $30$--$33\%$ across the explored host-mass range. These results should be interpreted as deterministic model predictions obtained by propagating environment-dependent ingredient ratios through two standard semi-analytic boost frameworks. We provide an environment-aware prescription for subhalo boosts, together with modular environmental corrections that may also be useful in indirect-detection forecasts, strong-lensing mass modeling, and related halo-population applications. △ Less

Submitted 28 April, 2026; originally announced April 2026.

Comments: 16 pages, 5 figures

Quasar photometric redshifts beyond the spectroscopic coverage: Uncertainty models and redshift distributions in the Kilo-Degree Survey DR5

Authors: Kacper Drabicki, Szymon J. Nakoneczny, Maciej Bilicki

Abstract: Photometric redshifts (photo-$z$) and their distributions underpin cosmology with photometric quasars, tracers in angular clustering and cross-correlations. Progress requires trustworthy uncertainties, especially beyond the spectroscopic training set. We compare how machine-learning frameworks estimate quasar photo-$z$ uncertainties and reconstruct the redshift distribution $n(z)$ under controlled… ▽ More Photometric redshifts (photo-$z$) and their distributions underpin cosmology with photometric quasars, tracers in angular clustering and cross-correlations. Progress requires trustworthy uncertainties, especially beyond the spectroscopic training set. We compare how machine-learning frameworks estimate quasar photo-$z$ uncertainties and reconstruct the redshift distribution $n(z)$ under controlled data-quality shifts. Using Kilo-Degree Survey DR5 photometry and DESI DR1 spectroscopic quasars, we train artificial neural networks (ANNs), Mixture Density Networks (MDNs) and Bayesian Neural Networks (BNNs) with Gaussian-mixture outputs, plus a self-organizing map (SOM) as a direct $n(z)$ estimator. We evaluate them on four subsets, with and without magnitude extrapolation and missing bands, through the negative log-likelihood, probability integral transform, point-estimate accuracy and the bias of binned $n(z)$ moments; degeneracies are sought by clustering the predicted PDFs. At least two mixture components are essential: a single Gaussian is miscalibrated and produces more catastrophic outliers than the ANN. No model performs best everywhere. On well-covered data the five-component MDN, three-component BNN and SOM reconstruct $n(z)$ almost perfectly; under faint extrapolation the BNN gives the best likelihoods, while in the hardest faint-plus-missing regime the single Gaussian becomes the best-calibrated. For point estimates and tomographic binning, uncertainty models outperform the ANN, while the SOM fails out-of-distribution. PDF clustering exposes distinct colour-redshift degeneracies likely to grow for fainter samples. The best model is thus regime- and application-dependent: multi-component MDNs or BNNs are needed for clean binning and are the only intrinsically calibrated choice on a well-covered golden sample, a first step towards a full comparison of quasar photo-$z$ pipelines. △ Less

Submitted 3 August, 2026; v1 submitted 20 March, 2026; originally announced March 2026.

EMU/GAMA: A statistical perspective on active galactic nuclei diagnostics

Authors: J. Prathap, A. M. Hopkins, R. Carvajal, M. Cowley, S. M. Croom, D. Farrah, I. Prandoni, S. S. Shabala, J. Th. van Loon, C. Pappalardo, K. A. Pimbblet, U. T. Ahmed, M. Bilicki, M. J. I. Brown, D. Leahy, A. Mailvaganam, J. R. Marvil, T. Mukherjee, S. F. Rahman, T. Vernstrom, J. Willingham, T. Zafar

Abstract: While it is well known that galaxies are composites of many emission processes, quantifying the various contributions remains challenging. In this work, we use unsupervised machine learning based clustering algorithms to evaluate the agreement between the clustering tools and astrophysical classifications, and hence quantify the fractional contributions of star formation processes and nuclear blac… ▽ More While it is well known that galaxies are composites of many emission processes, quantifying the various contributions remains challenging. In this work, we use unsupervised machine learning based clustering algorithms to evaluate the agreement between the clustering tools and astrophysical classifications, and hence quantify the fractional contributions of star formation processes and nuclear black hole activity to the total galaxy energy budget of radio sources. We perform clustering on the multiwavelength (optical, infrared (IR), and radio) active galactic nuclei (AGN) diagnostic spaces, using the data from the G09 and G23 fields from the Galaxy and Mass Assembly (GAMA) survey, Evolutionary Map of the Universe (EMU) survey, and the Wide-field Infrared Survey Explorer (WISE). We find that the statistical clustering recovers $\approx$ 90 % of the star forming galaxies (SFGs) and $\approx$ 80 % of the AGN. We define a new IR-radio AGN diagnostic scheme that identifies radio AGN from IR SFGs and AGN, corresponding to the KMeans cluster with approximately 90 % reliability. We demonstrate the superior power of radio AGN selection in higher dimensions using a three-dimensional space composed of directly observable parameters ($\rm W_1-W_2$ colour, $\rm W_2$ magnitude, and the 1.4 GHz radio flux density). This novel three dimensional diagnostic shows immense potential in radio AGN selection that is close to 90 % reliable and 90 % complete. We also publish a catalogue of radio sources in the EMU survey with associated probabilities for them to be active in the optical regime, through which we emphasise the philosophy of considering a galaxy to be composed of various fractions rather than a binary classification of SFGs and AGN. △ Less

Submitted 5 March, 2026; originally announced March 2026.

Comments: Accepted for publication in PASA. 19 pages, 9 figures, and 4 tables

The LOFAR Two-metre Sky Survey: VII. Third Data Release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, A. Drabent, A. Botteon, W. L. Williams, P. N. Best, H. J. A. Röttgering, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, F. De Gasperin, M. Haverkorn, C. Horellou, N. Jackson, G. K. Miley, L. K. Morabito, R. Morganti, S. P. O'Sullivan, D. J. Schwarz, D. J. B. Smith, R. J. van Weeren, H. K. Vedantham , et al. (77 additional authors not shown)

Abstract: We present the third data release of the LOFAR Two-metre Sky Survey (LoTSS-DR3). The survey images cover 88% of the northern sky and were created from 12,950 hrs of data (18.6 PB) accumulated over 10.5 years. The images were produced through direction-independent and direction-dependent calibration pipelines that correct for instrumental effects as well as spatially and temporally varying ionosphe… ▽ More We present the third data release of the LOFAR Two-metre Sky Survey (LoTSS-DR3). The survey images cover 88% of the northern sky and were created from 12,950 hrs of data (18.6 PB) accumulated over 10.5 years. The images were produced through direction-independent and direction-dependent calibration pipelines that correct for instrumental effects as well as spatially and temporally varying ionospheric distortions. In our 120-168 MHz continuum mosaic images with an angular resolution of 6 arcsec (9 arcsec below declination 10$^\circ$) we catalogue 13,667,877 sources, formed from 16,943,656 Gaussian components. The scatter in the astrometric precision approximately follows the expected noise-like behaviour but with an additional systematic component of at least 0.24 arcsec that is likely due to calibration imperfections. The random flux density scale error is 6%, while the systematic offset was previously shown to be within 2%. The median sensitivity of our mosaics is 92$μ$Jy beam$^{-1}$. Completeness simulations, accounting for realistic source models, time- and bandwidth-smearing effects, and astrometric errors, indicate that we detect more than 95% of compact sources with integrated flux densities exceeding 9 times the local root mean square (RMS) noise. However, the recovered source counts in a particular integrated flux density bin do not match the injected counts until flux densities exceed 45 times the local RMS noise. The Euclidean-normalised differential source counts derived from the survey constrain the radio source population over five orders of magnitude and are in good agreement with previous deep and wide-area surveys. All data products are publicly available, including catalogues, individual-field Stokes I, Q, U, and V images, mosaicked Stokes I images, and $uv$ data with associated direction-dependent calibration solutions. △ Less

Submitted 17 February, 2026; originally announced February 2026.

Comments: 15 figures, 2 tables and 22 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

KiDS-Legacy: WIMP dark matter constraints from the cross-correlation of weak lensing and Fermi-LAT gamma rays

Authors: Shiyang Zhang, Hendrik Hildebrandt, Ziang Yan, Tilman Tröster, Athithya Aravinthan, Marika Asgari, Deaglan J. Bartlett, Maciej Bilicki, Dominik Elsässer, Catherine Heymans, Benjamin Joachimi, Lauro Moscardini, Dennis Neumann, Anya Paopiamsap, Robert Reischke, Benjamin Stölzner

Abstract: Dark matter dominates the matter content of the Universe, and its properties can be constrained through large-scale structure probes such as the cross-correlation between the unresolved gamma-ray background (UGRB) and weak gravitational lensing. We analysed 15 years of Fermi-LAT data, constructing UGRB intensity maps in ten energy bins (0.5-1000 GeV), and cross-correlated them with KiDS-Legacy she… ▽ More Dark matter dominates the matter content of the Universe, and its properties can be constrained through large-scale structure probes such as the cross-correlation between the unresolved gamma-ray background (UGRB) and weak gravitational lensing. We analysed 15 years of Fermi-LAT data, constructing UGRB intensity maps in ten energy bins (0.5-1000 GeV), and cross-correlated them with KiDS-Legacy shear in six tomographic bins. The measurements were performed using angular power spectra estimated with the pseudo-$C_\ell$ method. No significant cross-correlation is found. Based on this non-detection, we present 95% upper bounds on the weakly interacting massive particle (WIMP) decay rate $Γ_{\rm dec}$ and velocity-averaged annihilation cross-section $\langleσ_{\rm ann} v\rangle$ as functions of mass. We compare our results with bounds from other cosmological tracers and from local probes, and found them to be complementary, particularly at low masses ($\rm GeV/TeV$). In addition, using a Euclid-like lensing survey cross-correlated with Fermi-LAT, we forecast $\sim$2 times tighter limits, highlighting the potential of forthcoming data to strengthen constraints on dark matter annihilation and decay. △ Less

Submitted 19 January, 2026; v1 submitted 16 January, 2026; originally announced January 2026.

Comments: 18 pages, 13 figures, submitted to A&A

Journal ref: A&A 710, A80 (2026)

KiDS-Legacy: Constraining dark energy, neutrino mass, and curvature

Authors: Robert Reischke, Benjamin Stölzner, Benjamin Joachimi, Angus H. Wright, Marika Asgari, Maciej Bilicki, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Joachim Harnois-Déraps, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Arthur Loureiro, Constance Mahony, Lauro Moscardini, Lucas Porth, Mario Radovich, Tilman Tröster, Maximilian von Wietersheim-Kramsta , et al. (3 additional authors not shown)

Abstract: We constrained minimally extended cosmological models with the cosmic shear analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy) in combination with external probes. Due to the consistency of the KiDS-Legacy analysis with the cosmic microwave background (CMB), we could combine these datasets reliably for the first time. Additionally, we used CMB lensing, galaxy redshift-spa… ▽ More We constrained minimally extended cosmological models with the cosmic shear analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy) in combination with external probes. Due to the consistency of the KiDS-Legacy analysis with the cosmic microwave background (CMB), we could combine these datasets reliably for the first time. Additionally, we used CMB lensing, galaxy redshift-space distortions, and baryon acoustic oscillations. We assessed, in turn, the effects of spatial curvature, varying neutrino masses, and an evolving dark energy component on cosmological constraints from KiDS-Legacy alone and from KiDS-Legacy combined with external probes. We find KiDS-Legacy to be consistent with the fiducial flat $Λ$-cold dark matter ($Λ$CDM) analysis with $c^2 \sum m_ν\leq 1.5\,$eV, $w_0 = -1.0\pm 0.7$, and $w_a = -1.3^{+1.9}_{-2.0}$ while $Ω_K = 0.08^{+0.16}_{-0.17}$ (1$σ$ bounds) with an almost equal goodness of fit. The $w_0w_a$CDM model is not a significant improvement over $Λ$CDM when cosmic shear and CMB lensing are combined, yielding a Bayes factor $B = 0.07$. If all probes are combined, however, $B$ increases to 2.73, corresponding to a $2.6σ$ suspiciousness tension. The constraint on $S_8 = σ_8\sqrt{Ω_\mathrm{m}/0.3}$ is robust to opening up the parameter space for cosmic shear. Adding all external datasets to KiDS-Legacy, we find $S_8 = 0.816 \pm 0.006$ in $Λ$CDM and $S_8 = 0.837 \pm 0.008$ in $w_0 w_a$CDM for all probes combined. △ Less

Submitted 10 May, 2026; v1 submitted 11 December, 2025; originally announced December 2025.

Comments: 14 pages, 8 figures, published in A&A

Journal ref: A&A 709, A82 (2026)

KiDS-Legacy: Constraints on Horndeski gravity from weak lensing combined with galaxy clustering and cosmic microwave background anisotropies

Authors: Benjamin Stölzner, Robert Reischke, Matteo Grasso, Matteo Cataneo, Benjamin Joachimi, Arthur Loureiro, Alessio Spurio Mancini, Angus H. Wright, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Constance Mahony, Lauro Moscardini, Lucas Porth, Mario Radovich, Tilman Tröster, Maximilian von Wietersheim-Kramsta , et al. (3 additional authors not shown)

Abstract: We present constraints on modified gravity from a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy) in combination with DESI measurements of baryon acoustic oscillations, eBOSS observations of redshift space distortions, and cosmic microwave background anisotropies from Planck. We study the Horndeski class of modified gravity models in an effective field theor… ▽ More We present constraints on modified gravity from a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy) in combination with DESI measurements of baryon acoustic oscillations, eBOSS observations of redshift space distortions, and cosmic microwave background anisotropies from Planck. We study the Horndeski class of modified gravity models in an effective field theory framework employing a parameterisation that satisfies stability conditions by construction and, for the first time, present a cosmological analysis in this inherently stable parameter basis. Cosmic shear constrains the Horndeski parameter space significantly, matching or surpassing the CMB contribution. Adopting the de-mixed kinetic term of the scalar field perturbation, $D_{\rm kin}$, and the deviation of the Planck mass from its fiducial value, $ΔM_*^2\equiv M_*^2-1$, as model parameters, we constrain their present values to be $Δ\hat{M}_*^2=0.32^{+0.07}_{-0.21}$ and $\hat{D}_{\rm kin} = 3.74^{+0.69}_{-1.92}$, which deviate from general relativity at $1.5σ$ and $1.9σ$, respectively. We derive constraints on the structure growth parameter $S_8=0.813^{+0.008}_{-0.011}$, which is compatible with the $Λ$CDM constraint at $0.54σ$. We obtain the deviation of the effective Newtonian coupling from the GR value as $Δμ_{\infty,{\rm eff}}=0.066\pm0.023$, corresponding to a $2.9σ$ significance. Although modified gravity provides a slightly better fit to the data, a model comparison shows only a weak preference for modified gravity at the $1.4σ$ level. When adopting a dynamical dark energy model of the background cosmology, the inferred modified gravity parameter constraints are stable with respect to a $Λ$CDM background, while a mild preference at $1.57σ$ for dynamical dark energy remains. △ Less

Submitted 24 March, 2026; v1 submitted 11 December, 2025; originally announced December 2025.

Comments: 21 pages, 13 figures, 5 tables, published in A&A

Journal ref: A&A 707, A323 (2026)

Angular clustering and bias of photometric quasars in the Kilo-Degree Survey Data Release 4

Authors: Anjitha John William, Maciej Bilicki, Wojciech A. Hellwing, Szymon J. Nakoneczny, Priyanka Jalan

Abstract: We investigate the angular clustering and effective bias of photometrically selected quasars in the Kilo-Degree Survey Data Release 4 (KiDS DR4). We update the previous photometric redshifts (photo-$z$s) of the KiDS quasars using Hybrid-z, a deep learning framework combining four-band KiDS images and nine-band KiDS+VIKING magnitudes. Hybrid-z is trained on the latest Dark Energy Spectroscopic Inst… ▽ More We investigate the angular clustering and effective bias of photometrically selected quasars in the Kilo-Degree Survey Data Release 4 (KiDS DR4). We update the previous photometric redshifts (photo-$z$s) of the KiDS quasars using Hybrid-z, a deep learning framework combining four-band KiDS images and nine-band KiDS+VIKING magnitudes. Hybrid-z is trained on the latest Dark Energy Spectroscopic Instrument (DESI) DR1 and Sloan Digital Sky Survey (SDSS) DR17 quasars matching with KiDS, and achieves average bias $\langle δz \rangle < 0.01$ and scatter $\sim 0.04(1 + z)$ on a test sample. The updated catalog of $\sim 157k$ quasars over $777~\mathrm{deg}^2$ is divided into four tomographic bins spanning $0.1 \leq z_{\mathrm{phot}} \leq 2.7$. In each bin, we measure the angular two-point correlation function and compare it with theoretical predictions for dark matter clustering. We estimate the best-fit scale-independent quasar bias, which increases from $b \approx 1.6$ at $z \approx 0.6$ to $b \approx 4.0$ at $z \approx 2.2$, and is well matched by a quadratic relation in redshift. Our clustering analysis indicates that KiDS quasars reside in dark matter halos of mass $\log_{10}(M_{\mathrm{eff}}/h^{-1}M_\odot)$ in the range $\sim 12.7$--$12.9$ and effective peak heights $ν_{\mathrm{eff}}$ rising from $\sim 1.5$ to $2.9$ over our redshift span. We study two systematics that could affect the bias derivation: stellar contamination and the redshift distribution assumed in the theoretical modeling. The former has a negligible effect, whereas the latter significantly impacts the derived $b(z)$, emphasizing the importance of redshift calibration. Our work is the first cosmological application of quasars selected from KiDS and paves the way for future extensions in the final KiDS DR5, the Legacy Survey of Space and Time, or the 4-metre Multi-Object Spectroscopic Telescope. △ Less

Submitted 21 November, 2025; originally announced November 2025.

Comments: 13 pages, 7 figures

Targeting cluster galaxies for the 4MOST CHANCES Low-z sub-survey with photometric redshifts

Authors: Hugo Méndez-Hernández, Ciria Lima-Dias, Antonela Monachesi, Yara L. Jaffé, Christopher P. Haines, Gabriel S. M. Teixeira, Elismar Lösch, Raúl Baier-Soto, Erik V. R. Lima, Amrutha B. M., C. R. Bom, Giuseppe D'Ago, Ricardo Demarco, Alexis Finoguenov, Rodrigo F. Haack, Amanda R. Lopes, C. Mendes de Oliveira, Paola Merluzzi, Franco Piraino-Cerda, Analía V. Smith Castelli, Cristobal Sif'on, Laerte Sodré Jr, Nicolás Tejos, Sergio Torres-Flores, Maria Argudo-Fernández , et al. (32 additional authors not shown)

Abstract: The evolution of galaxies is shaped by both internal processes and their external environments. Galaxy clusters and their surroundings provide ideal laboratories to study these effects, particularly mechanisms such as quenching and morphological transformation. The Chilean Cluster galaxy Evolution Survey (CHANCES) Low-z sub-survey is part of the CHileAN Cluster galaxy Evolution Survey, a 4MOST com… ▽ More The evolution of galaxies is shaped by both internal processes and their external environments. Galaxy clusters and their surroundings provide ideal laboratories to study these effects, particularly mechanisms such as quenching and morphological transformation. The Chilean Cluster galaxy Evolution Survey (CHANCES) Low-z sub-survey is part of the CHileAN Cluster galaxy Evolution Survey, a 4MOST community survey designed to uncover the relationship between the formation and evolution of galaxies and hierarchical structure formation as it happens, through deep and wide multi-object spectroscopy. We present the target selection strategy followed to select galaxy cluster candidate members for the CHANCES low-z sub-survey, in and around 50 clusters and two superclusters at z<0.07, out to (5XR200) and down to mr= 20.4. Combining public photometric redshift estimates from the DESI Legacy Imaging Survey and T80S/S-PLUS iDR5, with custom photometric redshifts, we identify likely galaxy cluster candidate members whose photometric redshifts are consistent with being at the known redshift of the cluster and measure the average deviations of their photometric redshifts with respect to the spectroscopic redshift measurements σNMAD. We have successfully compiled our CHANCES-low-redshift catalogues, split into three different sub-surveys: low-z bright (mr<18.5), low-z faint (18.5<=mr<20.4) and low-z faint supplementary, by selecting>= 500,000 galaxy cluster candidate members and including confirmed spectroscopic galaxy cluster members, from which we expect to obtain 4MOST low-resolution (R~6500) spectra for ~320,000 galaxies. The CHANCES Low-z target catalogues form a statistically robust sample for spectroscopic follow-up, allowing studies of galaxy evolution and environmental effects in nearby cluster and supercluster environments. △ Less

Submitted 22 October, 2025; originally announced October 2025.

Journal ref: A&A 706, A34 (2026)

Wide Area VISTA Extragalactic Survey (WAVES): Selection of targets for the Wide survey using decision-tree classification

Authors: G. Kaur, M. Bilicki, S. Bellstedt, E. Tempel, W. A. Hellwing, I. Baldry, B. Bandi, S. Barsanti, S. Driver, N. Guerra-Varas, B. Holwerda, C. Lagos, J. Loveday, A. Robotham

Abstract: The Wide-Area VISTA Extragalactic Survey (WAVES) on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) includes two flux-limited subsurveys with very high (95\%) completeness requirements: Wide over $\sim\!1200$ deg$^2$ and Deep over $\sim\!65$ deg$^2$. Both are $Z$-band selected, respectively as $Z<21.1$ and $Z<21.25$ mag, and additionally redshift-limited, while the true redshifts are not… ▽ More The Wide-Area VISTA Extragalactic Survey (WAVES) on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) includes two flux-limited subsurveys with very high (95\%) completeness requirements: Wide over $\sim\!1200$ deg$^2$ and Deep over $\sim\!65$ deg$^2$. Both are $Z$-band selected, respectively as $Z<21.1$ and $Z<21.25$ mag, and additionally redshift-limited, while the true redshifts are not known a priori but will be only measured by 4MOST. Here, we present a classification-based method to select the targets for WAVES-Wide. Rather than estimating individual redshifts for the input photometric objects, we assign probabilities of them being below $z=0.2$, the redshift limit of the subsurvey. This is done with the supervised machine learning approach of eXtreme Gradient Boosting (XGB), trained on a comprehensive spectroscopic sample overlapping with WAVES fields. Our feature space is composed of nine VST+VISTA magnitudes from $u$ to $K_s$ and all the possible colors, but most relevant for the classification are the $g$-band and the $u-g$, $g-r$ and $J-K_s$ colors. We check the performance of our classifier both for the fiducial WAVES-Wide limits, as well as for a range of neighboring redshift and magnitude thresholds, consistently finding purity and completeness at the level of 94-95\%. We note, however, that this performance deteriorates for sources close to the selection limits, due to deficiencies of the current spectroscopic training sample and the decreasing signal-to-noise of the photometry. We apply the classifier trained on the full spectroscopic sample to 14 million photometric galaxies from the WAVES input catalog, which have all 9 bands measured. Our work demonstrates that a machine-learning classifier could be used to select a flux- and redshift-limited sample from deep photometric data. △ Less

Submitted 13 October, 2025; originally announced October 2025.

Overdispersed radio source counts and excess radio dipole detection

Authors: Lukas Böhme, Dominik J. Schwarz, Prabhakar Tiwari, Morteza Pashapour-Ahmadabadi, Benedict Bahr-Kalus, Maciej Bilicki, Catherine L. Hale, Caroline S. Heneka, Thilo M. Siewert

Abstract: The source count dipole from wide-area radio continuum surveys allows us to test the cosmological standard model. Many radio sources have multiple components, which can cause an overdispersion of the source counts distribution. We account for this effect via a new Bayesian estimator, based on the negative binomial distribution. Combining the two best understood wide-area surveys, NVSS and RACS-low… ▽ More The source count dipole from wide-area radio continuum surveys allows us to test the cosmological standard model. Many radio sources have multiple components, which can cause an overdispersion of the source counts distribution. We account for this effect via a new Bayesian estimator, based on the negative binomial distribution. Combining the two best understood wide-area surveys, NVSS and RACS-low, and the deepest wide-area survey, LoTSS-DR2, we find that the source count dipole exceeds its expected value as the kinematic dipole amplitude from standard cosmology by a factor of $3.67 \pm 0.49$ -- a $5.4σ$ discrepancy. △ Less

Submitted 20 September, 2025; originally announced September 2025.

Comments: 10 pages, 3 figures. Accepted for publication in PRL

Journal ref: Phys. Rev. Lett. 135, 201001 (2025)

Matter power spectrum reconstruction with KiDS-Legacy: Improved internal $Λ$CDM consistency and preference for strong baryonic feedback

Authors: Jeger C. Broxterman, Patrick Simon, Lucas Porth, Konrad Kuijken, Angus H. Wright, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Shun-Sheng Li, Matteo Maturi, Lauro Moscardini, Mario Radovich, Robert Reischke, Maximilian Von Wietersheim-Kramsta

Abstract: Direct measurements of the matter power spectrum, $P_\mathrm{m}(k,z)$, provide a powerful tool to investigate observed tensions between models of structure growth while also testing the internal consistency of cosmological probes. We analyse cosmic shear data from the final data release of the Kilo-Degree Survey (KiDS), presenting a deprojected $P_\mathrm{m}(k,z)$, measured in up to three redshift… ▽ More Direct measurements of the matter power spectrum, $P_\mathrm{m}(k,z)$, provide a powerful tool to investigate observed tensions between models of structure growth while also testing the internal consistency of cosmological probes. We analyse cosmic shear data from the final data release of the Kilo-Degree Survey (KiDS), presenting a deprojected $P_\mathrm{m}(k,z)$, measured in up to three redshift bins. Compared to analyses using previous KiDS releases, we find improved internal consistency in the $z\lesssim0.7$ regime. At large scales, $k\lesssim0.1\,h\,\rm Mpc^{-1}$, our power spectrum reconstruction aligns with $Λ$CDM predictions with a density fluctuation amplitude $σ_8=0.81$. Furthermore, at small scales, $k=3$-$20\,h\,\rm Mpc^{-1}$, the average matter power spectrum is suppressed by $30\%\pm10\%\,{\rm (stat.)}\pm4\%\,{\rm (sys.)}$ with $2.8σ$ significance relative to a dark-matter-only model, consistent with expectations of strong baryonic feedback. △ Less

Submitted 8 October, 2025; v1 submitted 10 September, 2025; originally announced September 2025.

Comments: 8 pages, 5 figures (including the appendices). Accepted for publication in A&A

Journal ref: A&A 703, L3 (2025)

GWTC-4.0: Constraints on the Cosmic Expansion Rate and Modified Gravitational-wave Propagation

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, I. Abouelfettouh, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, D. Adhikari, N. Adhikari, R. X. Adhikari, V. K. Adkins, S. Afroz, A. Agapito, D. Agarwal, M. Agathos, N. Aggarwal, S. Aggarwal, O. D. Aguiar, I. -L. Ahrend, L. Aiello, A. Ain, P. Ajith, T. Akutsu , et al. (1750 additional authors not shown)

Abstract: We analyze data from 142 of the 218 gravitational-wave (GW) sources in the fourth LIGO-Virgo-KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-4.0) to estimate the Hubble constant $H_0$ jointly with the population properties of merging compact binaries. We measure the luminosity distance and redshifted masses of GW sources directly; in contrast, we infer GW source redshifts stat… ▽ More We analyze data from 142 of the 218 gravitational-wave (GW) sources in the fourth LIGO-Virgo-KAGRA Collaboration (LVK) Gravitational-Wave Transient Catalog (GWTC-4.0) to estimate the Hubble constant $H_0$ jointly with the population properties of merging compact binaries. We measure the luminosity distance and redshifted masses of GW sources directly; in contrast, we infer GW source redshifts statistically through i) location of features in the compact object mass spectrum and merger rate evolution, and ii) identifying potential host galaxies in the GW localization volume. Probing the relationship between source luminosity distances and redshifts obtained in this way yields constraints on cosmological parameters. We also constrain parameterized deviations from general relativity which affect GW propagation, specifically those modifying the dependence of a GW signal on the source luminosity distance. Assuming our fiducial model for the source-frame mass distribution and using GW candidates detected up to the end of the fourth observing run (O4a), together with the GLADE+ all-sky galaxy catalog, we estimate $H_0 = 76.6^{+13.0}_{-9.5} (76.6^{+25.2}_{-14.0})$ km s$^{-1}$ Mpc$^{-1}$. This value is reported as a median with 68.3% (90%) symmetric credible interval, and includes combination with the $H_0$ measurement from GW170817 and its electromagnetic counterpart. Using a parametrization of modified GW propagation in terms of the magnitude parameter $Ξ_0$, we estimate $Ξ_0 = 1.2^{+0.8}_{-0.4} (1.2^{+2.4}_{-0.5})$, where $Ξ_0 = 1$ recovers the behavior of general relativity. △ Less

Submitted 7 October, 2025; v1 submitted 4 September, 2025; originally announced September 2025.

Comments: As part of the Astrophysical Journal Letters Focus Issue on the Gravitational Wave Transient Catalog

Report number: LIGO-P2400152

The 4MOST-Cosmology Redshift Survey: Target Selection of Bright Galaxies and Luminous Red Galaxies

Authors: Aurelien Verdier, Antoine Rocher, Behnood Bandi, Johan Richard, Boudewijn Roukema, Jon Loveday, Elmo Tempel, Maciej Bilicki, Jean-Paul Kneib, Mathilde Guitton

Abstract: The Cosmology Redshift Survey of the 4-metre Multi-Object Spectroscopic Telescope (4MOST-CRS) will provide redshift measurements of galaxies and quasars over 5700{\degsq} in the southern hemisphere. As targets for the 4MOST-CRS, we present a selection of an $r<19.25$ magnitude limited sample of Bright Galaxies (BG) and a colour selected sample of Luminous Red Galaxies (LRG) based on DESI Legacy Su… ▽ More The Cosmology Redshift Survey of the 4-metre Multi-Object Spectroscopic Telescope (4MOST-CRS) will provide redshift measurements of galaxies and quasars over 5700{\degsq} in the southern hemisphere. As targets for the 4MOST-CRS, we present a selection of an $r<19.25$ magnitude limited sample of Bright Galaxies (BG) and a colour selected sample of Luminous Red Galaxies (LRG) based on DESI Legacy Survey DR10.1 photometric data, in the redshift ranges $0.1<z<0.5$ and $0.4<z<1$, respectively. These samples are selected using the $g$, $r$, $z$, and $W1$ (from unWISE) photometric bands. For BGs, the star--galaxy separation is performed based on Gaia and Tycho-2 star catalogues. Following 4MOST requirements, the target densities of BGs and LRGs are 250{\mdegsq} and 400{\mdegsq}, respectively. We quantified the stellar contamination to be $<2-4\%$ for both galaxy samples using data from the first DESI data release. Using angular clustering, we show that both samples are robust against imaging systematics, confirming a low stellar contamination. Finally, we provide forecasts of the baryonic acoustic oscillation (BAO) and growth of structure ($σ_8$) measurements from the 4MOST-CRS alone and in combination with DESI. From the 4MOST-CRS, we expect 3\% and 25\% precision on BAO and $σ_8$ measurements, respectively. Combining 4MOST-CRS with DESI improves the constraints from DESI alone by 12$\%$ in the 0.4 < z < 0.8 redshift range, leading to the most stringent constraints on BAO measurements from spectroscopic galaxy clustering. △ Less

Submitted 10 August, 2025; originally announced August 2025.

Journal ref: Mon Not R Astron Soc (2025)

AMICO galaxy clusters in KiDS-1000: cosmological constraints and mass calibration from counts and weak lensing

Authors: G. F. Lesci, F. Marulli, L. Moscardini, M. Maturi, M. Sereno, M. Radovich, M. Romanello, C. Giocoli, A. H. Wright, S. Bardelli, M. Bilicki, G. Castignani, H. Hildebrandt, L. Ingoglia, S. Joudaki, A. Kannawadi, E. Puddu

Abstract: We present the joint modelling of weak-lensing and count measurements of the galaxy clusters detected with the AMICO code, in the fourth data release of the Kilo Degree Survey (KiDS-1000). The analysed sample comprises about 8000 clusters, covering an effective area of 839 deg$^{2}$ and extending up to a redshift of $z = 0.8$. Stacked cluster weak-lensing and count measurements have been derived i… ▽ More We present the joint modelling of weak-lensing and count measurements of the galaxy clusters detected with the AMICO code, in the fourth data release of the Kilo Degree Survey (KiDS-1000). The analysed sample comprises about 8000 clusters, covering an effective area of 839 deg$^{2}$ and extending up to a redshift of $z = 0.8$. Stacked cluster weak-lensing and count measurements have been derived in bins of redshift and intrinsic richness, $λ^*$. Based on self-organising maps, we reconstructed the true redshift distributions of the background galaxy samples. We accounted for the systematic uncertainties arising from impurities in the background and cluster samples, biases in the cluster $z$ and $λ^*$, projection effects, halo orientation and miscentring, truncation of cluster halo mass distributions, matter correlated with cluster haloes, multiplicative shear bias, baryonic matter, geometric distortions in the lensing profiles, uncertainties in the theoretical halo mass function, and super-sample covariance. We also employed a blinding strategy based on perturbing the cluster sample completeness. The improved statistics and photometry compared to the previous KiDS data release, KiDS-DR3, have led to a halving of the uncertainties on $Ω_{\rm m}$ and $σ_8$, as we obtained $Ω_{\rm m}=0.22\pm0.02$ and $σ_8=0.86\pm0.03$. The constraint on $S_8 \equiv σ_8(Ω_{\rm m}/0.3)^{0.5}$, $S_8=0.74\pm0.03$, is in excellent agreement with recent cluster count and KiDS-1000 cosmic shear analyses, while it shows a $2.8σ$ tension with Planck cosmic microwave background results. The constraints on the $\logλ^*-\log M_{200}$ relation imply a mass precision of 8%, on average. In addition, the result on the intrinsic scatter of the $\logλ^*-\log M_{200}$ relation, $σ_{\rm intr}=0.05\pm0.02$, confirms that $λ^*$ is an excellent mass proxy. △ Less

Submitted 25 September, 2025; v1 submitted 18 July, 2025; originally announced July 2025.

Comments: 26 pages, 19 figures. Submitted to A&A. We addressed the first referee report

Journal ref: A&A 703, A25 (2025)

Anisotropic Counts-in-Cells in Redshift Space: A New Route to Cosmological Constraints from Galaxy Surveys

Authors: Paweł Drozda, Wojciech A. Hellwing, Maciej Bilicki

Abstract: We introduce a novel extension of the volume-averaged correlation function (VACF) framework by replacing the traditional spherical smoothing kernels with anisotropic, ellipsoidal windows. This generalized approach enables the study of shape-dependent clustering statistics and captures directional information encoded in large-scale structure, particularly in redshift space where galaxy distribution… ▽ More We introduce a novel extension of the volume-averaged correlation function (VACF) framework by replacing the traditional spherical smoothing kernels with anisotropic, ellipsoidal windows. This generalized approach enables the study of shape-dependent clustering statistics and captures directional information encoded in large-scale structure, particularly in redshift space where galaxy distribution is distorted by peculiar velocities. We define and compute ellipsoidal VACFs $\barξ_J (r_{\parallel}, r_{\perp})$ and the corresponding reduced cumulants $s_J (r_{\parallel}, r_{\perp})$, allowing for joint sensitivity to both scale and anisotropy across arbitrary statistical order J. Using a suite of COLA N-body simulations spanning a grid of cosmologies with varying $Ω_M$ and $σ_8$, we analyze the behavior of ellipsoidal VACFs and cumulants in both real and redshift space. We find that the shape of the smoothing kernel that maximizes the clustering signal depends strongly on the redshift-space distortion regime: spherical in real space, prolate in the Fingers-of-God-dominated regime, and oblate in the Kaiser squashing-dominated regime. While the standard VACF amplitude is mainly sensitive to $σ_8$, the shape-dependence of redshift-space skewness shows a coherent response to the combined growth parameter $f σ_8$, with a typical sensitivity at the 1-3 $σ$ level between neighboring models. Our results demonstrate that ellipsoidal VACFs offer a computationally efficient and information-rich generalization of counts-in-cells analysis, with promising applications to galaxy survey data, halo catalogs, and cosmological tests of gravity beyond $ΛCDM$. △ Less

Submitted 12 June, 2026; v1 submitted 2 June, 2025; originally announced June 2025.

Comments: 12 pages, 8 figures, 1 table. Code available at: https://github.com/Pawel-96/Avcorr

Journal ref: Phys. Rev. D 112, 063543, 2026

EMU/GAMA: A new approach to characterising radio luminosity functions

Authors: J. Prathap, A. M. Hopkins, J. Afonso, M. Bilicki, M. Cowley, S. M. Croom, Y. Gordon, S. Phillipps, E. M. Sadler, S. S. Shabala, U. T. Ahmed, S. Amarantidis, M. J. I. Brown, R. Carvajal, D. Leahy, J. R. Marvil, T. Mukherjee, J. Willingham, T. Zafar

Abstract: This study characterises the radio luminosity functions (RLFs) for SFGs and AGN using statistical redshift estimation in the absence of comprehensive spectroscopic data. Sensitive radio surveys over large areas detect many sources with faint optical and infrared counterparts, for which redshifts and spectra are unavailable. This challenges our attempt to understand the population of radio sources.… ▽ More This study characterises the radio luminosity functions (RLFs) for SFGs and AGN using statistical redshift estimation in the absence of comprehensive spectroscopic data. Sensitive radio surveys over large areas detect many sources with faint optical and infrared counterparts, for which redshifts and spectra are unavailable. This challenges our attempt to understand the population of radio sources. Statistical tools are often used to model parameters (such as redshift) as an alternative to observational data. Using the data from GAMA G23 and EMU early science observations, we explore simple statistical techniques to estimate the redshifts in order to measure the RLFs of the G23 radio sources as a whole and for SFGs and AGN separately. Redshifts and AGN/SFG classifications are assigned statistically for those radio sources without spectroscopic data. The calculated RLFs are compared with existing studies, and the results suggest that the RLFs match remarkably well for low redshift galaxies with an optical counterpart. We use a more realistic high redshift distribution to model the redshifts of (most likely) high redshift radio sources and find that the LFs from our approach match well with measured LFs. We also look at strategies to compare the RLFs of radio sources without an optical counterpart to existing studies. △ Less

Submitted 16 May, 2025; originally announced May 2025.

Comments: 15 Figures, 2 Tables, accepted for publication in PASA

The Evolutionary Map of the Universe: A new radio atlas for the southern hemisphere sky

Authors: A. M. Hopkins, A. Kapinska, J. Marvil, T. Vernstrom, J. D. Collier, R. P. Norris, Y. A. Gordon, S. W. Duchesne, L. Rudnick, N. Gupta, E. Carretti, C. S. Anderson, S. Dai, G. Gürkan, D. Parkinson, I. Prandoni, S. Riggi, C. S. Saraf, Y. K. Ma, M. D. Filipović, G. Umana, B. Bahr-Kalus, B. S. Koribalski, E. Lenc, A. Ingallinera , et al. (48 additional authors not shown)

Abstract: We present the Evolutionary Map of the Universe (EMU) survey conducted with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU aims to deliver the touchstone radio atlas of the southern hemisphere. We introduce EMU and review its science drivers and key science goals, updated and tailored to the current ASKAP five-year survey plan. The development of the survey strategy and planned sky… ▽ More We present the Evolutionary Map of the Universe (EMU) survey conducted with the Australian Square Kilometre Array Pathfinder (ASKAP). EMU aims to deliver the touchstone radio atlas of the southern hemisphere. We introduce EMU and review its science drivers and key science goals, updated and tailored to the current ASKAP five-year survey plan. The development of the survey strategy and planned sky coverage is presented, along with the operational aspects of the survey and associated data analysis, together with a selection of diagnostics demonstrating the imaging quality and data characteristics. We give a general description of the value-added data pipeline and data products before concluding with a discussion of links to other surveys and projects and an outline of EMU's legacy value. △ Less

Submitted 13 May, 2025; originally announced May 2025.

Comments: 32 pages, 15 figures. Accepted for publication in PASA

EMU: Cross-correlating EMU Pilot Survey 1 with Dark Energy Survey to validate the radio galaxy bias and redshift distribution

Authors: Chandra Shekhar Saraf, David Parkinson, Jacobo Asorey, Catherine L. Hale, Benedict Bahr-Kalus, Maciej Bilicki, Stefano Camera, Andrew M. Hopkins, Konstantinos Tanidis

Abstract: Radio continuum galaxy surveys can provide a relatively fast map of the projected distribution of structure in the Universe, at the cost of lacking information about the radial distribution. We can use these surveys to learn about the growth of structure and the fundamental physics of the Universe, but doing so requires extra information to be provided in the modelling of the redshift distribution… ▽ More Radio continuum galaxy surveys can provide a relatively fast map of the projected distribution of structure in the Universe, at the cost of lacking information about the radial distribution. We can use these surveys to learn about the growth of structure and the fundamental physics of the Universe, but doing so requires extra information to be provided in the modelling of the redshift distribution, $dN/dz$. In this work, we show how the cross-correlation of the two-dimensional radio continuum map with another galaxy map (in this case a photometric optical extragalactic survey), with a known redshift distribution, can be used to determine the redshift distribution through statistical inference. We use data from the Evolutionary Map of the Universe (EMU) Pilot Survey 1 and cross-correlate it with optical data from the Dark Energy Survey to fit the parameters of our $dN/dz$ model. We show that the recovered distribution has a similar shape to the distribution predicted by the current state-of-the-art simulation, and can fit the angular power spectrum data from cross-correlations very well, validating this model. These results will have significance for future cosmological analyses with large-scale radio continuum surveys such as the full EMU, or with the SKAO. △ Less

Submitted 24 June, 2026; v1 submitted 9 May, 2025; originally announced May 2025.

Comments: 18 pages, 8 figures, 5 tables. Prepared for submission to PASA

Skewness as a Probe of Gravity: Real and Redshift Space Counts-In-Cells

Authors: Paweł Drozda, Wojciech A. Hellwing, Maciej Bilicki

Abstract: We study the counts-in-cells reduced skewness $s_3$ for dark matter, halo, and galaxy distributions in both real and redshift space, using the ELEPHANT ($\textit{Extended LEnsing PHysics with ANalytical ray Tracing}$) suite of $N$-body simulations. We compare General Relativity (GR) with two extended (EG) gravity models: $f(R)$ gravity with chameleon screening and the normal-branch Dvali-Gabadadze… ▽ More We study the counts-in-cells reduced skewness $s_3$ for dark matter, halo, and galaxy distributions in both real and redshift space, using the ELEPHANT ($\textit{Extended LEnsing PHysics with ANalytical ray Tracing}$) suite of $N$-body simulations. We compare General Relativity (GR) with two extended (EG) gravity models: $f(R)$ gravity with chameleon screening and the normal-branch Dvali-Gabadadze-Porrati (nDGP) model with Vainshtein screening. We quantify the suppression of $s_3$ by redshift-space distortions (RSD), finding that while small-scale skewness is strongly reduced, the $F5$ model retains a $\sim 4\%$ deviation from GR in galaxy samples, corresponding to a $2σ$ significance. We show that the ratio $s_3^{\mathrm{RSD}}/s_3^{\mathrm{real}}$ is approximately independent of the gravity model across tracers and redshifts. Our results demonstrate that real-space predictions can help reliably infer redshift-space skewness in both GR and extended gravity, providing a new tool for testing gravity with current and forthcoming galaxy redshift surveys. △ Less

Submitted 12 June, 2026; v1 submitted 17 April, 2025; originally announced April 2025.

Comments: 14 pages, 9 figures, 1 table

Journal ref: Phys. Rev. D 112, 063546, 2025

The CosmoVerse White Paper: Addressing observational tensions in cosmology with systematics and fundamental physics

Authors: Eleonora Di Valentino, Jackson Levi Said, Adam Riess, Agnieszka Pollo, Vivian Poulin, Adrià Gómez-Valent, Amanda Weltman, Antonella Palmese, Caroline D. Huang, Carsten van de Bruck, Chandra Shekhar Saraf, Cheng-Yu Kuo, Cora Uhlemann, Daniela Grandón, Dante Paz, Dominique Eckert, Elsa M. Teixeira, Emmanuel N. Saridakis, Eoin Ó Colgáin, Florian Beutler, Florian Niedermann, Francesco Bajardi, Gabriela Barenboim, Giulia Gubitosi, Ilaria Musella , et al. (516 additional authors not shown)

Abstract: The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-t… ▽ More The standard model of cosmology has provided a good phenomenological description of a wide range of observations both at astrophysical and cosmological scales for several decades. This concordance model is constructed by a universal cosmological constant and supported by a matter sector described by the standard model of particle physics and a cold dark matter contribution, as well as very early-time inflationary physics, and underpinned by gravitation through general relativity. There have always been open questions about the soundness of the foundations of the standard model. However, recent years have shown that there may also be questions from the observational sector with the emergence of differences between certain cosmological probes. In this White Paper, we identify the key objectives that need to be addressed over the coming decade together with the core science projects that aim to meet these challenges. These discordances primarily rest on the divergence in the measurement of core cosmological parameters with varying levels of statistical confidence. These possible statistical tensions may be partially accounted for by systematics in various measurements or cosmological probes but there is also a growing indication of potential new physics beyond the standard model. After reviewing the principal probes used in the measurement of cosmological parameters, as well as potential systematics, we discuss the most promising array of potential new physics that may be observable in upcoming surveys. We also discuss the growing set of novel data analysis approaches that go beyond traditional methods to test physical models. [Abridged] △ Less

Submitted 4 August, 2025; v1 submitted 2 April, 2025; originally announced April 2025.

Comments: 416 pages, 81 figures

Journal ref: Phys. Dark Univ. 49 (2025) 101965

KiDS-Legacy: Consistency of cosmic shear measurements and joint cosmological constraints with external probes

Authors: Benjamin Stölzner, Angus H. Wright, Marika Asgari, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Constance Mahony, Robert Reischke, Mijin Yoon, Maciej Bilicki, Pierre Burger, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Joachim Harnois-Déraps, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Giorgio Francesco Lesci, Laila Linke, Arthur Loureiro , et al. (11 additional authors not shown)

Abstract: We present a cosmic shear consistency analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy). By adopting three tiers of consistency metrics, we compare cosmological constraints between subsets of the KiDS-Legacy dataset split by redshift, angular scale, galaxy colour and spatial region. We also review a range of two-point cosmic shear statistics. With the data passing all ou… ▽ More We present a cosmic shear consistency analysis of the final data release from the Kilo-Degree Survey (KiDS-Legacy). By adopting three tiers of consistency metrics, we compare cosmological constraints between subsets of the KiDS-Legacy dataset split by redshift, angular scale, galaxy colour and spatial region. We also review a range of two-point cosmic shear statistics. With the data passing all our consistency metric tests, we demonstrate that KiDS-Legacy is the most internally consistent KiDS catalogue to date. In a joint cosmological analysis of KiDS-Legacy and DES Y3 cosmic shear, combined with data from the Pantheon+ Type Ia supernovae compilation and baryon acoustic oscillations from DESI Y1, we find constraints consistent with Planck measurements of the cosmic microwave background with $S_8\equiv σ_8\sqrt{Ω_{\rm m}/0.3} = 0.814^{+0.011}_{-0.012}$ and $σ_8 = 0.802^{+0.022}_{-0.018}$. △ Less

Submitted 20 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 25 pages, 16 figures, 7 tables, published in A&A

Journal ref: A&A 702, A169 (2025)

KiDS-Legacy: Cosmological constraints from cosmic shear with the complete Kilo-Degree Survey

Authors: Angus H. Wright, Benjamin Stölzner, Marika Asgari, Maciej Bilicki, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Henk Hoekstra, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Robert Reischke, Maximilian von Wietersheim-Kramsta, Mijin Yoon, Pierre Burger, Nora Elisa Chisari, Jelte de Jong, Andrej Dvornik, Christos Georgiou, Joachim Harnois-Déraps, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Giorgio Francesco Lesci, Laila Linke , et al. (13 additional authors not shown)

Abstract: We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-in… ▽ More We present cosmic shear constraints from the completed Kilo-Degree Survey (KiDS), where the cosmological parameter $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3} = 0.815^{+0.016}_{-0.021}$, is found to be in agreement ($0.73σ$) with results from the Planck Legacy cosmic microwave background experiment. The final KiDS footprint spans $1347$ square degrees of deep nine-band imaging across the optical and near-infrared, along with an extra $23$ square degrees of KiDS-like calibration observations of deep spectroscopic surveys. Improvements in our redshift distribution estimation methodology, combined with our enhanced calibration data and multi-band image simulations, allow us to extend our lensed sample out to a photometric redshift of $z_{\rm B}\leq2.0$. Compared to previous KiDS analyses, the increased survey area and redshift depth results in a $\sim32\%$ improvement in constraining power in terms of $Σ_8\equivσ_8\left(Ω_{\rm m}/0.3\right)^α= 0.821^{+0.014}_{-0.016}$, where $α= 0.58$ has been optimised to match the revised degeneracy direction of $σ_8$ and $Ω_{\rm m}$. We adopt a new physically motivated intrinsic alignment model that depends jointly on the galaxy sample's halo mass and spectral type distributions, and that is informed by previous direct alignment measurements. We also marginalise over our uncertainty on the impact of baryon feedback on the non-linear matter power spectrum. Comparing to previous KiDS analyses, we conclude that the increase seen in $S_8$ primarily results from our improved redshift distribution estimation and calibration, as well as new survey area and improved image reduction. Our companion paper Stölzner et al. (submitted) presents a full suite of internal and external consistency tests, finding the KiDS-Legacy data set to be the most internally robust sample produced by KiDS to date. △ Less

Submitted 21 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 38 pages, 24 figures, 7 tables, accepted for publication in A&A

Journal ref: A&A 703, A158 (2025)

KiDS-Legacy: Redshift distributions and their calibration

Authors: Angus H. Wright, Hendrik Hildebrandt, Jan Luca van den Busch, Maciej Bilicki, Catherine Heymans, Benjamin Joachimi, Constance Mahony, Robert Reischke, Benjamin Stölzner, Anna Wittje, Marika Asgari, Nora Elisa Chisari, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Henk Hoekstra, Priyanka Jalan, Anjitha John William, Shahab Joudaki, Konrad Kuijken, Giorgio Francesco Lesci, Shun-Sheng Li, Laila Linke, Arthur Loureiro, Matteo Maturi , et al. (8 additional authors not shown)

Abstract: We present the redshift calibration methodology and bias estimates for the cosmic shear analysis of the fifth and final data release (DR5) of the Kilo-Degree Survey (KiDS). KiDS-DR5 includes a greatly expanded compilation of calibrating spectra, drawn from $27$ square degrees of dedicated optical and near-IR imaging taken over deep spectroscopic fields. The redshift distribution calibration levera… ▽ More We present the redshift calibration methodology and bias estimates for the cosmic shear analysis of the fifth and final data release (DR5) of the Kilo-Degree Survey (KiDS). KiDS-DR5 includes a greatly expanded compilation of calibrating spectra, drawn from $27$ square degrees of dedicated optical and near-IR imaging taken over deep spectroscopic fields. The redshift distribution calibration leverages a range of new methods and updated simulations to produce the most precise $N(z)$ bias estimates used by KiDS to date. Improvements to our colour-based redshift distribution measurement method (SOM) mean that we are able to use many more sources per tomographic bin for our cosmological analyses, and better estimate the representation of our source sample given the available spec-$z$. We validate our colour-based redshift distribution estimates with spectroscopic cross-correlations (CC). We find that improvements to our cross-correlation redshift distribution measurement methods mean that redshift distribution biases estimated between the SOM and CC methods are fully consistent on simulations, and the data calibration is consistent to better than $2σ$ in all tomographic bins. △ Less

Submitted 21 October, 2025; v1 submitted 25 March, 2025; originally announced March 2025.

Comments: 26 pages, 18 figures, 8 tables, accepted for publication in A&A

Journal ref: A&A 703, A144 (2025)

The fifth data release of the Kilo Degree Survey: Multi-epoch optical/NIR imaging covering wide and legacy-calibration fields

Authors: Angus H. Wright, Konrad Kuijken, Hendrik Hildebrandt, Mario Radovich, Maciej Bilicki, Andrej Dvornik, Fedor Getman, Catherine Heymans, Henk Hoekstra, Shun-Sheng Li, Lance Miller, Nicola R. Napolitano, Qianli Xia, Marika Asgari, Massimo Brescia, Hugo Buddelmeijer, Pierre Burger, Gianluca Castignani, Stefano Cavuoti, Jelte de Jong, Alastair Edge, Benjamin Giblin, Carlo Giocoli, Joachim Harnois-Déraps, Priyanka Jalan , et al. (29 additional authors not shown)

Abstract: We present the final data release of the Kilo-Degree Survey (KiDS-DR5), a public European Southern Observatory (ESO) wide-field imaging survey optimised for weak gravitational lensing studies. We combined matched-depth multi-wavelength observations from the VLT Survey Telescope and the VISTA Kilo-degree INfrared Galaxy (VIKING) survey to create a nine-band optical-to-near-infrared survey spanning… ▽ More We present the final data release of the Kilo-Degree Survey (KiDS-DR5), a public European Southern Observatory (ESO) wide-field imaging survey optimised for weak gravitational lensing studies. We combined matched-depth multi-wavelength observations from the VLT Survey Telescope and the VISTA Kilo-degree INfrared Galaxy (VIKING) survey to create a nine-band optical-to-near-infrared survey spanning $1347$ deg$^2$. The median $r$-band $5σ$ limiting magnitude is 24.8 with median seeing $0.7^{\prime\prime}$. The main survey footprint includes $4$ deg$^2$ of overlap with existing deep spectroscopic surveys. We complemented these data in DR5 with a targeted campaign to secure an additional $23$ deg$^2$ of KiDS- and VIKING-like imaging over a range of additional deep spectroscopic survey fields. From these fields, we extracted a catalogue of $126\,085$ sources with both spectroscopic and photometric redshift information, which enables the robust calibration of photometric redshifts across the full survey footprint. In comparison to previous releases, DR5 represents a $34\%$ areal extension and includes an $i$-band re-observation of the full footprint, thereby increasing the effective $i$-band depth by $0.4$ magnitudes and enabling multi-epoch science. Our processed nine-band imaging, single- and multi-band catalogues with masks, and homogenised photometry and photometric redshifts can be accessed through the ESO Archive Science Portal. △ Less

Submitted 25 March, 2025; originally announced March 2025.

Comments: 60 pages, 43 figures, 20 tables, published in A&A

Journal ref: A&A, 686, A170 (2024)

Target Selection for the Redshift-Limited WAVES-Wide with Machine Learning

Authors: Gursharanjit Kaur, Maciej Bilicki, Wojciech Hellwing, the WAVES team

Abstract: The forthcoming Wide Area Vista Extragalactic Survey (WAVES) on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) has a key science goal of probing the halo mass function to lower limits than possible with previous surveys. For that purpose, in its Wide component, galaxies targetted by WAVES will be flux-limited to $Z<21.1$ mag and will cover the redshift range of $z<0.2$, at a spectroscopi… ▽ More The forthcoming Wide Area Vista Extragalactic Survey (WAVES) on the 4-metre Multi-Object Spectroscopic Telescope (4MOST) has a key science goal of probing the halo mass function to lower limits than possible with previous surveys. For that purpose, in its Wide component, galaxies targetted by WAVES will be flux-limited to $Z<21.1$ mag and will cover the redshift range of $z<0.2$, at a spectroscopic success rate of $\sim95\%$. Meeting this completeness requirement, when the redshift is unknown a priori, is a challenge. We solve this problem with supervised machine learning to predict the probability of a galaxy falling within the WAVES-Wide redshift limit, rather than estimate each object's redshift. This is done by training an XGBoost tree-based classifier to decide if a galaxy should be a target or not. Our photometric data come from 9-band VST+VISTA observations, including KiDS+VIKING surveys. The redshift labels for calibration are derived from an extensive spectroscopic sample overlapping with KiDS and ancillary fields. Our current results indicate that with our approach, we should be able to achieve the completeness of $\sim95\%$, which is the WAVES success criterion. △ Less

Submitted 28 February, 2025; originally announced February 2025.

Comments: 5 pages, 1 figure, ML4ASTRO2 conference proceedings

Intrinsic galaxy alignments in the KiDS-1000 bright sample: dependence on colour, luminosity, morphology, and galaxy scale

Authors: Christos Georgiou, Nora Elisa Chisari, Maciej Bilicki, Francesco La Barbera, Nicola R. Napolitano, Nivya Roy, Crescenzo Tortora

Abstract: The intrinsic alignment (IA) of galaxies is a major astrophysical contaminant to weak gravitational lensing measurements, and the study of its dependence on galaxy properties helps provide meaningful physical priors that aid cosmological analyses. This work studied for the first time the dependence of IA on galaxy structural parameters. We measured the IA of bright galaxies, selected on apparent r… ▽ More The intrinsic alignment (IA) of galaxies is a major astrophysical contaminant to weak gravitational lensing measurements, and the study of its dependence on galaxy properties helps provide meaningful physical priors that aid cosmological analyses. This work studied for the first time the dependence of IA on galaxy structural parameters. We measured the IA of bright galaxies, selected on apparent r-band magnitude r<20, in the Kilo-Degree Survey (KiDS). Machine-learning-based photometric redshift estimates are available for this galaxy sample that helped us obtain a clean measurement of its IA signal. We supplemented this sample with a catalogue of structural parameters from Sersic profile fits to the surface-brightness profiles of the galaxies. We split the sample on galaxy intrinsic colour, luminosity, and Sersic index, and we fitted the non-linear linear alignment model to galaxy position-shape projected correlation function measurements on large scales. We observe a power-law luminosity dependence of the large-scale IA amplitude, $A_{IA}$, for both the red and high-Sersic-index ($n_s>2.5$) samples, and find no significant difference between the two. We measure an $\sim1.5σ$ lower $A_{IA}$ for red galaxies that also have a Sersic index of $n_s<4$ compared to the expected amplitude predicted using the sample's luminosity. We also probe the IA of red galaxies as a function of galaxy scale by varying the radial weight employed in the shape measurement. On large scales (above 6 Mpc/$h$), we do not detect a significant difference in the alignment. On smaller scales, we observe that IA increase with galaxy scale, with outer galaxy regions showing stronger alignments than inner regions. Finally, for intrinsically blue galaxies, we find $A_{IA}=-0.67\pm1.00$, which is consistent with previous works, and we find IA to be consistent with zero for the low-Sersic-index ($n_s<2.5$) sample. △ Less

Submitted 2 July, 2025; v1 submitted 13 February, 2025; originally announced February 2025.

Comments: 14 pages, 9 figures

Journal ref: A&A 699, A252 (2025)

Hybrid-z: Enhancing Kilo-Degree Survey bright galaxy sample photometric redshifts with deep learning

Authors: Anjitha John William, Priyanka Jalan, Maciej Bilicki, Wojciech A. Hellwing, Hareesh Thuruthipilly, Szymon J. Nakoneczny

Abstract: We employ deep learning (DL) to improve photometric redshifts (photo-$z$s) in the Kilo-Degree Survey Data Release 4 Bright galaxy sample (KiDS-Bright DR4). This dataset, used as a foreground for KiDS lensing and clustering studies, is flux-limited to $r<20$ mag with mean $z=0.23$ and covers 1000 deg$^2$. Its photo-$z$s were previously derived with artificial neural networks from the ANNz2 package,… ▽ More We employ deep learning (DL) to improve photometric redshifts (photo-$z$s) in the Kilo-Degree Survey Data Release 4 Bright galaxy sample (KiDS-Bright DR4). This dataset, used as a foreground for KiDS lensing and clustering studies, is flux-limited to $r<20$ mag with mean $z=0.23$ and covers 1000 deg$^2$. Its photo-$z$s were previously derived with artificial neural networks from the ANNz2 package, trained on the Galaxy And Mass Assembly (GAMA) spectroscopy. Here we considerably improve over these previous redshift estimations by building a DL model, Hybrid-z, which combines four-band KiDS images with nine-band magnitudes from KiDS+VIKING. The Hybrid-z framework provides photo-$z$s for KiDS-Bright, with negligible mean residuals of O($10^{-4}$) and scatter at the level of $0.014(1+z)$ -- reduction by 20% over the previous nine-band derivations with ANNz2. We check our photo-$z$ model performance on test data drawn from GAMA, as well as from other KiDS-overlapping wide-angle spectroscopic surveys, namely SDSS, 2dFLenS, and 2dFGRS. We find stable behavior and consistent improvement over ANNz2 throughout. We finally apply Hybrid-z trained on GAMA to the entire KiDS-Bright DR4 sample of 1.2 million galaxies. For these final predictions, we design a method of smoothing the input redshift distribution of the training set, to avoid propagation of features present in GAMA, related to its small sky area and large-scale structure imprint in its fields. Our work paves the way towards the best-possible photo-$z$s achievable with machine learning for any galaxy type both for the final KiDS-Bright DR5 data and for future deeper imaging, such as from the Legacy Survey of Space and Time. △ Less

Submitted 16 July, 2025; v1 submitted 3 January, 2025; originally announced January 2025.

Comments: 14 pages, 7 figures

Journal ref: A&A 698, A276 (2025)

Cross-correlating the EMU Pilot Survey 1 with CMB lensing: Constraints on cosmology and galaxy bias with harmonic-space power spectra

Authors: K. Tanidis, J. Asorey, C. S. Saraf, C. L. Hale, B. Bahr-Kalus, D. Parkinson, S. Camera, R. P. Norris, A. M. Hopkins, M. Bilicki, N. Gupta

Abstract: We measured the harmonic-space power spectrum of galaxy clustering auto-correlation from the Evolutionary Map of the Universe Pilot Survey 1 data (EMU PS1) and its cross-correlation with the lensing convergence map of cosmic microwave background (CMB) from Planck Public Release 4 at the linear scale range from $\ell=2$ to 500. We applied two flux density cuts at $0.18$ and $0.4$mJy on the radio ga… ▽ More We measured the harmonic-space power spectrum of galaxy clustering auto-correlation from the Evolutionary Map of the Universe Pilot Survey 1 data (EMU PS1) and its cross-correlation with the lensing convergence map of cosmic microwave background (CMB) from Planck Public Release 4 at the linear scale range from $\ell=2$ to 500. We applied two flux density cuts at $0.18$ and $0.4$mJy on the radio galaxies observed at 944MHz and considered two source detection algorithms. We found the auto-correlation measurements from the two algorithms at the 0.18mJy cut to deviate for $\ell\gtrsim250$ due to the different criteria assumed on the source detection and decided to ignore data above this scale. We report a cross-correlation detection of EMU PS1 with CMB lensing at $\sim$5.5$σ$, irrespective of flux density cut. In our theoretical modelling we considered two redshift distribution simulation models that yield consistent results, a linear and a non-linear matter power spectrum, and two linear galaxy bias models. That is a constant redshift-independent galaxy bias $b(z)=b_g$ and a constant amplitude galaxy bias $b(z)=b_g/D(z)$. By fixing a cosmology model and considering a non-linear matter power spectrum, we measured a constant galaxy bias at $0.18$mJy ($0.4$mJy) with $b_g=2.32^{+0.41}_{-0.33}$ ($2.18^{+0.17}_{-0.25}$) and a constant amplitude bias with $b_g=1.72^{+0.31}_{-0.21}$ ($1.78^{+0.22}_{-0.15}$). When $σ_8$ is a free parameter for the same models at $0.18$mJy ($0.4$mJy) with the constant model we found $σ_8=0.68^{+0.16}_{-0.14}$ ($0.82\pm0.10$), while with the constant amplitude model we measured $σ_8=0.61^{+0.18}_{-0.20}$ ($0.78^{+0.11}_{-0.09}$), respectively. Our results agree at $1σ$ with the measurements from Planck CMB and the weak lensing surveys and also show the potential of cosmology studies with future radio continuum survey data. △ Less

Submitted 20 November, 2025; v1 submitted 8 November, 2024; originally announced November 2024.

Comments: 23 pages, 9 figures, 4 tables. Version matching publication at journal

Journal ref: Publications of the Astronomical Society of Australia, Volume 42, 2025, e062

KiDS-Legacy: Angular galaxy clustering from deep surveys with complex selection effects

Authors: Ziang Yan, Angus H. Wright, Nora Elisa Chisari, Christos Georgiou, Shahab Joudaki, Arthur Loureiro, Robert Reischke, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Catherine Heymans, Hendrik Hildebrandt, Priyanka Jalan, Benjamin Joachimi, Giorgio Francesco Lesci, Shun-Sheng Li, Laila Linke, Constance Mahony, Lauro Moscardini, Nicola R. Napolitano, Benjamin Stoelzner, Maximilian Von Wietersheim-Kramsta, Mijin Yoon

Abstract: Photometric galaxy surveys, despite their limited resolution along the line of sight, encode rich information about the large-scale structure (LSS) of the Universe thanks to the high number density and extensive depth of the data. However, the complicated selection effects in wide and deep surveys can potentially cause significant bias in the angular two-point correlation function (2PCF) measured… ▽ More Photometric galaxy surveys, despite their limited resolution along the line of sight, encode rich information about the large-scale structure (LSS) of the Universe thanks to the high number density and extensive depth of the data. However, the complicated selection effects in wide and deep surveys can potentially cause significant bias in the angular two-point correlation function (2PCF) measured from those surveys. In this paper, we measure the 2PCF from the newly published KiDS-Legacy sample. Given an $r$-band $5σ$ magnitude limit of $24.8$ and survey footprint of $1347$ deg$^2$, it achieves an excellent combination of sky coverage and depth for such a measurement. We find that complex selection effects, primarily induced by varying seeing, introduce over-estimation of the 2PCF by approximately an order of magnitude. To correct for such effects, we apply a machine learning-based method to recover an organised random (OR) that presents the same selection pattern as the galaxy sample. The basic idea is to find the selection-induced clustering of galaxies using a combination of self-organising maps (SOMs) and hierarchical clustering (HC). This unsupervised machine learning method is able to recover complicated selection effects without specifying their functional forms. We validate this SOM+HC method on mock deep galaxy samples with realistic systematics and selections derived from the KiDS-Legacy catalogue. Using mock data, we demonstrate that the OR delivers unbiased 2PCF cosmological parameter constraints, removing the $27σ$ offset in the galaxy bias parameter that is recovered when adopting uniform randoms. Blinded measurements on the real KiDS-Legacy data show that the corrected 2PCF is robust to the SOM+HC configuration near the optimal set-up suggested by the mock tests. △ Less

Submitted 3 February, 2025; v1 submitted 30 October, 2024; originally announced October 2024.

Comments: 29 pages, 27 figures, 4 tables; Accepted for publication on Astronomy & Astrophysics; The code used for this work is published on https://github.com/yanzastro/tiaogeng

Journal ref: A&A 694, A259 (2025)

KiDS-Legacy: Covariance validation and the unified OneCovariance framework for projected large-scale structure observables

Authors: Robert Reischke, Sandra Unruh, Marika Asgari, Andrej Dvornik, Hendrik Hildebrandt, Benjamin Joachimi, Lucas Porth, Maximilian von Wietersheim-Kramsta, Jan Luca van den Busch, Benjamin Stölzner, Angus H. Wright, Ziang Yan, Maciej Bilicki, Pierre Burger, Nora Elisa Chisari, Joachim Harnois-Deraps, Christos Georgiou, Catherine Heymans, Priyanka Jalan, Shahab Joudaki, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Constance Mahony, Davide Sciotti , et al. (2 additional authors not shown)

Abstract: We introduce OneCovariance, an open-source software designed to accurately compute covariance matrices for an arbitrary set of two-point summary statistics across a variety of large-scale structure tracers. Utilising the halo model, we estimated the statistical properties of matter and biased tracer fields, incorporating all Gaussian, non-Gaussian, and super-sample covariance terms. The flexible c… ▽ More We introduce OneCovariance, an open-source software designed to accurately compute covariance matrices for an arbitrary set of two-point summary statistics across a variety of large-scale structure tracers. Utilising the halo model, we estimated the statistical properties of matter and biased tracer fields, incorporating all Gaussian, non-Gaussian, and super-sample covariance terms. The flexible configuration permits user-specific parameters, such as the complexity of survey geometry, the halo occupation distribution employed to define each galaxy sample, or the form of the real-space and/or Fourier space statistics to be analysed. We illustrate the capabilities of OneCovariance within the context of a cosmic shear analysis of the final data release of the Kilo-Degree Survey (KiDS-Legacy). Upon comparing our estimated covariance with measurements from mock data and calculations from independent software, we ascertain that OneCovariance achieves accuracy at the per cent level. When assessing the impact of ignoring complex survey geometry in the cosmic shear covariance computation, we discover misestimations at approximately the $10\%$ level for cosmic variance terms. Nonetheless, these discrepancies do not significantly affect the KiDS-Legacy recovery of cosmological parameters. We derive the cross-covariance between real-space correlation functions, bandpowers, and COSEBIs, facilitating future consistency tests among these three cosmic shear statistics. Additionally, we calculate the covariance matrix of photometric-spectroscopic galaxy clustering measurements, validating the jackknife covariance estimates for calibrating KiDS-Legacy redshift distributions. The OneCovariance can be found on GitHub, together with comprehensive documentation and examples. △ Less

Submitted 16 September, 2025; v1 submitted 9 October, 2024; originally announced October 2024.

Comments: 37 pages, 11 figures, published in A&A, code available at https://github.com/rreischke/OneCovariance

Journal ref: A&A 699, A124 (2025)

6x2pt: Forecasting gains from joint weak lensing and galaxy clustering analyses with spectroscopic-photometric galaxy cross-correlations

Authors: Harry Johnston, Nora Elisa Chisari, Shahab Joudaki, Robert Reischke, Benjamin Stölzner, Arthur Loureiro, Constance Mahony, Sandra Unruh, Angus H. Wright, Marika Asgari, Maciej Bilicki, Pierre Burger, Andrej Dvornik, Christos Georgiou, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Konrad Kuijken, Shun-Sheng Li, Laila Linke, Lucas Porth, HuanYuan Shan, Tilman Tröster, Jan Luca van den Busch , et al. (3 additional authors not shown)

Abstract: We explore the enhanced self-calibration of photometric galaxy redshift distributions, $n(z)$, through the combination of up to six two-point functions. Our $\rm 3\times2pt$ configuration is comprised of photometric shear, spectroscopic galaxy clustering, and spectroscopic-photometric galaxy-galaxy lensing (GGL). We further include spectroscopic-photometric cross-clustering; photometric GGL; and p… ▽ More We explore the enhanced self-calibration of photometric galaxy redshift distributions, $n(z)$, through the combination of up to six two-point functions. Our $\rm 3\times2pt$ configuration is comprised of photometric shear, spectroscopic galaxy clustering, and spectroscopic-photometric galaxy-galaxy lensing (GGL). We further include spectroscopic-photometric cross-clustering; photometric GGL; and photometric auto-clustering, using the photometric shear sample as density tracer. We perform simulated likelihood forecasts of the cosmological and nuisance parameter constraints for Stage-III- and Stage-IV-like surveys. For the Stage-III-like case, we employ realistic but perturbed redshift distributions, and distinguish between "coherent" shifting in one direction, versus more internal scattering and full-shape errors. For perfectly known $n(z)$, a $\rm 6\times2pt$ analysis gains $\sim40\%$ in Figure of Merit (FoM) in the $S_8\equivσ_8\sqrt{Ω_{\rm m}/0.3}$ and $Ω_{\rm m}$ plane relative to the $\rm 3\times2pt$ analysis. If untreated, coherent and incoherent redshift errors lead to inaccurate inferences of $S_8$ and $Ω_{\rm m}$, respectively. Employing bin-wise scalar shifts $δ{z}_i$ in the tomographic mean redshifts reduces cosmological parameter biases, with a $\rm 6x2pt$ analysis constraining the shift parameters with $2-4$ times the precision of a photometric $\rm 3^{ph}\times2pt$ analysis. For the Stage-IV-like survey, a $\rm 6\times2pt$ analysis doubles the FoM($σ_8{-}Ω_{\rm m}$) compared to any $\rm 3\times2pt$ or $\rm 3^{ph}\times2pt$ analysis, and is only $8\%$ less constraining than if the $n(z)$ were perfectly known. A Gaussian mixture model for the $n(z)$ reduces mean-redshift errors and preserves the $n(z)$ shape. It also yields the most accurate and precise cosmological constraints for any $N\rm\times2pt$ configuration given $n(z)$ biases. △ Less

Submitted 25 September, 2024; originally announced September 2024.

Comments: 38 pages, 20 figures, to be submitted to A&A

Journal ref: A&A 699, A127 (2025)

KiDS-1000: Weak lensing and intrinsic alignment around luminous red galaxies

Authors: Maria Cristina Fortuna, Andrej Dvornik, Henk Hoekstra, Nora Elisa Chisari, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Koen Kuijken, Angus H. Wright, Ji Yao

Abstract: We study the properties of luminous red galaxies (LRGs) selected from the fourth data release of the Kilo Degree Survey (KiDS-1000) via galaxy-galaxy lensing of the background galaxies from KiDS-1000. We used a halo model formalism to interpret our measurements and obtain estimates of the halo masses as well as the satellite fractions of the LRGs, resulting in halo masses of… ▽ More We study the properties of luminous red galaxies (LRGs) selected from the fourth data release of the Kilo Degree Survey (KiDS-1000) via galaxy-galaxy lensing of the background galaxies from KiDS-1000. We used a halo model formalism to interpret our measurements and obtain estimates of the halo masses as well as the satellite fractions of the LRGs, resulting in halo masses of $2.7 \times 10^{12} h^{-1} {\rm M}_{\odot}<M_{\rm h}< 2.6 \times 10^{13} h^{-1} {\rm M}_{\odot}$. We studied the strength of intrinsic alignments (IA) using the position-shape correlations as a function of LRG luminosity, where we used a double power law to describe the relation between luminosity and halo mass to allow for a comparison with previous works. Here, we directly linked the observed IA of the (central) galaxy to the mass of the hosting halo, which is expected to be a fundamental quantity in establishing the alignment. We find that the dependence of the IA amplitude on halo mass is described well by a single power law, with an amplitude of $A = 5.74\pm{0.32}$ and slope of $β_M = 0.44\pm{0.04}$, in the range of $1.9 \times 10^{12}h^{-1} {\rm M}_{\odot}<M_{\rm h}<3.7 \times 10^{14} h^{-1} {\rm M}_{\odot}$. We also find that both red and blue galaxies from the source sample associated with the LRGs are randomly oriented, with respect to the LRGs, although our detection significance is limited by the uncertainty in our photometric redshifts. △ Less

Submitted 4 February, 2025; v1 submitted 23 September, 2024; originally announced September 2024.

Comments: 14 pages, 6 figures, accepted for publication by A&A

Journal ref: A&A 694, A322 (2025)

Enhancing Photometric Redshift Catalogs Through Color-Space Analysis: Application to KiDS-Bright Galaxies

Authors: Priyanka Jalan, Maciej Bilicki, Wojciech A. Hellwing, Angus H. Wright, Andrej Dvornik, Catherine Heymans, Hendrik Hildebrandt, Shahab Joudaki, Konrad Kuijken, Constance Mahony, Szymon Jan Nakoneczny, Mario Radovich, Jan Luca van den Busch, Mijin Yoon

Abstract: We present a method to refine photometric redshift galaxy catalogs by comparing their color-space matching with overlapping spectroscopic calibration data. We focus on cases where photometric redshifts (photo-$z$) are estimated empirically. Identifying galaxies that are poorly represented in spectroscopic data is crucial, as their photo-$z$ may be unreliable due to extrapolation beyond the trainin… ▽ More We present a method to refine photometric redshift galaxy catalogs by comparing their color-space matching with overlapping spectroscopic calibration data. We focus on cases where photometric redshifts (photo-$z$) are estimated empirically. Identifying galaxies that are poorly represented in spectroscopic data is crucial, as their photo-$z$ may be unreliable due to extrapolation beyond the training sample. Our approach uses a self-organizing map (SOM) to project a multi-dimensional parameter space of magnitudes and colors onto a 2-D manifold, allowing us to analyze the resulting patterns as a function of various galaxy properties. Using SOM, we compare the Kilo-Degree Survey bright galaxy sample (KiDS-Bright), limited to $r<20$ mag, with various spectroscopic samples, including the Galaxy And Mass Assembly (GAMA). Our analysis reveals that GAMA under-represents KiDS-Bright at its faintest ($r\gtrsim19.5$) and highest-redshift ($z\gtrsim0.4$) ranges, however no strong trends in color or stellar mass. By incorporating additional spectroscopic data from the SDSS, 2dF, and early DESI, we identify SOM cells where photo-$z$ are estimated suboptimally. We derive a set of SOM-based criteria to refine the photometric sample and improve photo-$z$ statistics. For the KiDS-Bright sample, this improvement is modest: exclusion of the least represented 20% of the sample reduces photo-$z$ scatter by less than 10%. We conclude that GAMA, used for KiDS-Bright photo-$z$ training, is sufficiently representative for reliable redshift estimation across most of the color space. Future spectroscopic data from surveys such as DESI should be better suited for exploiting the full improvement potential of our method. △ Less

Submitted 19 December, 2024; v1 submitted 23 September, 2024; originally announced September 2024.

Comments: Published in A&A

Caught in the cosmic web: environmental effects on subhalo abundance and internal density profiles

Authors: Feven Markos Hunde, Oliver Newton, Wojciech A. Hellwing, Maciej Bilicki, Krishna Naidoo

Abstract: Using the high-resolution $N$-body cosmological simulation COLOR, we explore the cosmic web (CW) environmental effects on subhalo populations and their internal properties. We use CaCTus, which incorporates an implementation of the state-of-the-art segmentation method NEXUS+, to delineate the simulation volume into nodes, filaments, walls, and voids. We group host haloes by virial mass and segment… ▽ More Using the high-resolution $N$-body cosmological simulation COLOR, we explore the cosmic web (CW) environmental effects on subhalo populations and their internal properties. We use CaCTus, which incorporates an implementation of the state-of-the-art segmentation method NEXUS+, to delineate the simulation volume into nodes, filaments, walls, and voids. We group host haloes by virial mass and segment each mass bin into consecutive CW elements. This reveals that subhalo populations in hosts within specific environments differ on average from the cosmic mean. The subhalo mass function is affected strongly, where hosts in filaments typically contain more subhaloes ($5$ to $20\%$), while hosts in voids are subhalo-poor, with $25\%$ fewer subhaloes. We find that the abundance of the most massive subhaloes, with reduced masses of $μ\equiv M_\mathrm{sub}/M_{200}$ is most sensitive to the CW environment. A corresponding picture emerges when looking at subhalo mass fractions, $f_\mathrm{sub}$, where the filament hosts are significantly more `granular' (having higher $f_\mathrm{sub}$) than the cosmic mean, while the void hosts have much smoother density distributions (with $f_\mathrm{sub}$ lower by $2$ - $20\%$ than the mean). Finally, when we look at the subhalo internal kinematic $V_{\rm max}$-$R_{\rm max}$ relations, we find that subhaloes located in the void and wall hosts exhibit density profiles with lower concentrations than the mean, while the filament hosts demonstrate much more concentrated mass profiles. Across all our samples, the effect of the CW environment generally strengthens with decreasing host halo virial mass. Our results show that host location in the large-scale CW introduces significant systematic effects on internal subhalo properties and population statistics. △ Less

Submitted 26 June, 2025; v1 submitted 13 September, 2024; originally announced September 2024.

Comments: 14 pages, 5 figures, and 4 tables; accepted for publication in Astronomy & Astrophysics after major revision

Journal ref: Astronomy & Astrophysics, Volume 700, Article A65 (2025)

Wide Area VISTA Extra-galactic Survey (WAVES): Unsupervised star-galaxy separation on the WAVES-Wide photometric input catalogue using UMAP and ${\rm{\scriptsize HDBSCAN}}$

Authors: Todd L. Cook, Behnood Bandi, Sam Philipsborn, Jon Loveday, Sabine Bellstedt, Simon P. Driver, Aaron S. G. Robotham, Maciej Bilicki, Gursharanjit Kaur, Elmo Tempel, Ivan Baldry, Daniel Gruen, Marcella Longhetti, Angela Iovino, Benne W. Holwerda, Ricardo Demarco

Abstract: Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an… ▽ More Star-galaxy separation is a crucial step in creating target catalogues for extragalactic spectroscopic surveys. A classifier biased towards inclusivity risks including spurious stars, wasting fibre hours, while a more conservative classifier might overlook galaxies, compromising completeness and hence survey objectives. To avoid bias introduced by a training set in supervised methods, we employ an unsupervised machine learning approach. Using photometry from the Wide Area VISTA Extragalactic Survey (WAVES)-Wide catalogue comprising 9-band $u-K_s$ data, we create a feature space with colours, fluxes, and apparent size information extracted by ${\rm P{\scriptsize RO} F{\scriptsize OUND}}$. We apply the non-linear dimensionality reduction method UMAP (Uniform Manifold Approximation and Projection) combined with the classifier ${\rm{\scriptsize HDBSCAN}}$ to classify stars and galaxies. Our method is verified against a baseline colour and morphological method using a truth catalogue from Gaia, SDSS, GAMA, and DESI. We correctly identify 99.72% of galaxies within the AB magnitude limit of $Z = 21.2$, with an F1 score of $0.9971 \pm 0.0018$ across the entire ground truth sample, compared to $0.9879 \pm 0.0088$ from the baseline method. Our method's higher purity ($0.9967 \pm 0.0021$) compared to the baseline ($0.9795 \pm 0.0172$) increases efficiency, identifying 11% fewer galaxy or ambiguous sources, saving approximately 70,000 fibre hours on the 4MOST instrument. We achieve reliable classification statistics for challenging sources including quasars, compact galaxies, and low surface brightness galaxies, retrieving 95.1%, 84.6%, and 99.5% of them respectively. Angular clustering analysis validates our classifications, showing consistency with expected galaxy clustering, regardless of the baseline classification. △ Less

Submitted 18 November, 2024; v1 submitted 17 June, 2024; originally announced June 2024.

Comments: Accepted by MNRAS; 20 pages, 14 figures

Journal ref: MNRAS 535(3), 2024, 2129 2148

EMU/GAMA: A Technique for Detecting Active Galactic Nuclei in Low Mass Systems

Authors: Jahang Prathap, Andrew M. Hopkins, Aaron S. G. Robotham, Sabine Bellstedt, José Afonso, Ummee T. Ahmed, Maciej Bilicki, Malcolm N. Bremer, Sarah Brough, Michael J. I. Brown, Yjan Gordon, Benne W. Holwerda, Denis Leahy, Ángel R. López-Sánchez, Joshua R. Marvil, Tamal Mukherjee, Isabella Prandoni, Stanislav S. Shabala, Tessa Vernstrom, Tayyaba Zafar

Abstract: We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe… ▽ More We propose a new method for identifying active galactic nuclei (AGN) in low mass ($\rm M_*\leq10^{10}M_\odot$) galaxies. This method relies on spectral energy distribution (SED) fitting to identify galaxies whose radio flux density has an excess over that expected from star formation alone. Combining data in the Galaxy and Mass Assembly (GAMA) G23 region from GAMA, Evolutionary Map of the Universe (EMU) early science observations, and Wide-field Infrared Survey Explorer (WISE), we compare this technique with a selection of different AGN diagnostics to explore the similarities and differences in AGN classification. We find that diagnostics based on optical and near-infrared criteria (the standard BPT diagram, the WISE colour criterion, and the mass-excitation, or MEx diagram) tend to favour detection of AGN in high mass, high luminosity systems, while the ``ProSpect'' SED fitting tool can identify AGN efficiently in low mass systems. We investigate an explanation for this result in the context of proportionally lower mass black holes in lower mass galaxies compared to higher mass galaxies and differing proportions of emission from AGN and star formation dominating the light at optical and infrared wavelengths as a function of galaxy stellar mass. We conclude that SED-derived AGN classification is an efficient approach to identify low mass hosts with low radio luminosity AGN. △ Less

Submitted 18 February, 2024; originally announced February 2024.

Comments: 16 pages, 12 figures, 5 tables. Accepted for publication in PASA

EMU/GAMA: Radio detected galaxies are more obscured than optically selected galaxies

Authors: U. T. Ahmed, A. M. Hopkins, J. Ware, Y. A. Gordon, M. Bilicki, M. J. I. Brown, M. Cluver, G. Gürkan, Á. R. López-Sánchez, D. A. Leahy, L. Marchetti, S. Phillipps, I. Prandoni, N. Seymour, E. N. Taylor, E. Vardoulaki

Abstract: We demonstrate the importance of radio selection in probing heavily obscured galaxy populations. We combine Evolutionary Map of the Universe (EMU) Early Science data in the Galaxy and Mass Assembly (GAMA) G23 field with the GAMA data, providing optical photometry and spectral line measurements, together with Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry, providing IR luminosi… ▽ More We demonstrate the importance of radio selection in probing heavily obscured galaxy populations. We combine Evolutionary Map of the Universe (EMU) Early Science data in the Galaxy and Mass Assembly (GAMA) G23 field with the GAMA data, providing optical photometry and spectral line measurements, together with Wide-field Infrared Survey Explorer (WISE) infrared (IR) photometry, providing IR luminosities and colours. We investigate the degree of obscuration in star forming galaxies, based on the Balmer decrement (BD), and explore how this trend varies, over a redshift range of 0<z<0.345. We demonstrate that the radio detected population has on average higher levels of obscuration than the parent optical sample, arising through missing the lowest BD and lowest mass galaxies, which are also the lower star formation rate (SFR) and metallicity systems. We discuss possible explanations for this result, including speculation around whether it might arise from steeper stellar initial mass functions in low mass, low SFR galaxies. △ Less

Submitted 19 December, 2023; originally announced December 2023.

Comments: Accepted for publication in PASA, 17 pages, 14 figures, 3 tables

Deep learning based photometric redshifts for the Kilo-Degree Survey Bright Galaxy Sample

Authors: Anjitha John William, Priyanka Jalan, Maciej Bilicki, Wojciech A. Hellwing

Abstract: In cosmological analyses, precise redshift determination remains pivotal for understanding cosmic evolution. However, with only a fraction of galaxies having spectroscopic redshifts (spec-$z$s), the challenge lies in estimating redshifts for a larger number. To address this, photometry-based redshift (photo-$z$) estimation, employing machine learning algorithms, is a viable solution. Identifying t… ▽ More In cosmological analyses, precise redshift determination remains pivotal for understanding cosmic evolution. However, with only a fraction of galaxies having spectroscopic redshifts (spec-$z$s), the challenge lies in estimating redshifts for a larger number. To address this, photometry-based redshift (photo-$z$) estimation, employing machine learning algorithms, is a viable solution. Identifying the limitations of previous methods, this study focuses on implementing deep learning (DL) techniques within the Kilo-Degree Survey (KiDS) Bright Galaxy Sample for more accurate photo-$z$ estimations. Comparing our new DL-based model against prior `shallow' neural networks, we showcase improvements in redshift accuracy. Our model gives mean photo-$z$ bias $\langle Δz\rangle= 10^{-3}$ and scatter $\mathrm{SMAD}(Δz)=0.016$, where $Δz = (z_\mathrm{phot}-z_\mathrm{spec})/(1+z_\mathrm{spec})$. This research highlights the promising role of DL in revolutionizing photo-$z$ estimation. △ Less

Submitted 14 December, 2023; v1 submitted 13 December, 2023; originally announced December 2023.

Comments: Proceedings of the XLI Polish Astronomical Society Meeting, 11-15 September 2023, Toruń, Poland

Journal ref: Proceedings of the Polish Astronomical Society, vol. 13, 243-248 (2024), https://www.pta.edu.pl/proc/v13p243