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

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.9970 across the entire ground truth sample, compared to 0.9871 from the baseline method. Our method's higher purity (0.9966) compared to the baseline (0.9780) 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 17 June, 2024; originally announced June 2024.

Comments: Submitted to MNRAS; 17 pages, 11 figures

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

Dynamics of pairwise motions in the fully non-linear regime in LCDM and Modified Gravity cosmologies

Authors: Mariana Jaber, Wojciech A. Hellwing, Jorge E. García-Farieta, Suhani Gupta, Maciej Bilicki

Abstract: In contrast to our understanding of density field tracers, the modelling of direct statistics pertaining to the cosmic velocity field remains open to significant opportunities for improvement. The lack of accurate modelling for the non-linear domain of pairwise velocities restricts our capacity to fully exploit the information encoded in this observable. We present a robust approach for modelling… ▽ More In contrast to our understanding of density field tracers, the modelling of direct statistics pertaining to the cosmic velocity field remains open to significant opportunities for improvement. The lack of accurate modelling for the non-linear domain of pairwise velocities restricts our capacity to fully exploit the information encoded in this observable. We present a robust approach for modelling the mean infall velocities, $v_{12}(r,a)$, with broad applicability spanning sub-megaparsec scales and cosmologies extending beyond the standard LCDM paradigm. Our approach involves solving the full pair-conservation equation using accurate non-linear power spectrum descriptions. To assess the robustness of our model, we extend it to cosmologies beyond the standard LCDM, in particular, the Hu-Sawicki $f(R)$-gravity and Dvali-Gabadadze-Porrati (DGP) modified gravity models. Remarkably, our predictions for pairwise velocities of dark matter particles at kilo-parsec scales exhibit excellent agreement with N-body simulations throughout the entire dynamical range ($0.1 \lesssim ξ\lesssim 1000$, or $r\geq0.4$Mpc/h). Furthermore, we show that different gravity models leave distinct signatures in the shape and dynamics of the mean pairwise velocities, providing a potent test of cosmological gravity laws. △ Less

Submitted 1 December, 2023; originally announced December 2023.

Comments: 10 pages plus references, 5 figures. Our main findings are in Figs. 1, 2 and 4. Submitted to PRD for publication. Comments welcome

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

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

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

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

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

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

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

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

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

Submitted 11 October, 2023; originally announced October 2023.

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

Search for Eccentric Black Hole Coalescences during the Third Observing Run of LIGO and Virgo

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, A. G. Abac, R. Abbott, H. Abe, F. Acernese, K. Ackley, C. Adamcewicz, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, I. Aguilar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi , et al. (1750 additional authors not shown)

Abstract: Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effect… ▽ More Despite the growing number of confident binary black hole coalescences observed through gravitational waves so far, the astrophysical origin of these binaries remains uncertain. Orbital eccentricity is one of the clearest tracers of binary formation channels. Identifying binary eccentricity, however, remains challenging due to the limited availability of gravitational waveforms that include effects of eccentricity. Here, we present observational results for a waveform-independent search sensitive to eccentric black hole coalescences, covering the third observing run (O3) of the LIGO and Virgo detectors. We identified no new high-significance candidates beyond those that were already identified with searches focusing on quasi-circular binaries. We determine the sensitivity of our search to high-mass (total mass $M>70$ $M_\odot$) binaries covering eccentricities up to 0.3 at 15 Hz orbital frequency, and use this to compare model predictions to search results. Assuming all detections are indeed quasi-circular, for our fiducial population model, we place an upper limit for the merger rate density of high-mass binaries with eccentricities $0 < e \leq 0.3$ at $0.33$ Gpc$^{-3}$ yr$^{-1}$ at 90\% confidence level. △ Less

Submitted 7 August, 2023; originally announced August 2023.

Comments: 24 pages, 5 figures

Report number: LIGO-P2300080

Constraints on dark matter and astrophysics from tomographic $γ$-ray cross-correlations

Authors: Anya Paopiamsap, David Alonso, Deaglan J. Bartlett, Maciej Bilicki

Abstract: We study the cross-correlation between maps of the unresolved $γ$-ray background constructed from the 12-year data release of the Fermi Large-Area Telescope, and the overdensity of galaxies in the redshift range $z\lesssim0.4$ as measured by the 2MASS Photometric Redshift survey and the WISE-SuperCOSMOS photometric survey. A signal is detected at the $8-10σ$ level, which we interpret in terms of b… ▽ More We study the cross-correlation between maps of the unresolved $γ$-ray background constructed from the 12-year data release of the Fermi Large-Area Telescope, and the overdensity of galaxies in the redshift range $z\lesssim0.4$ as measured by the 2MASS Photometric Redshift survey and the WISE-SuperCOSMOS photometric survey. A signal is detected at the $8-10σ$ level, which we interpret in terms of both astrophysical $γ$-ray sources, and WIMP dark matter decay and annihilation. The sensitivity achieved allows us to characterise the energy and redshift dependence of the signal, and we show that the latter is incompatible with a pure dark matter origin. We thus use our measurement to place an upper bound on the WIMP decay rate and the annihilation cross-section, finding constraints that are competitive with those found in other analyses. Our analysis is based on the extraction of clean model-independent observables that can then be used to constrain arbitrary astrophysical and particle physics models. In this sense we produce measurements of the $γ$-ray emissivity as a function of redshift and rest-frame energy $ε$, and of a quantity $F(ε)$ encapsulating all WIMP parameters relevant for dark matter decay or annihilation. We make these measurements, together with a full account of their statistical uncertainties, publicly available. △ Less

Submitted 9 May, 2024; v1 submitted 27 July, 2023; originally announced July 2023.

Comments: 27 pages, 13 figures. Data and code available in https://github.com/anyabua/FermiX

KiDS-1000: Cosmology with improved cosmic shear measurements

Authors: Shun-Sheng Li, Henk Hoekstra, Konrad Kuijken, Marika Asgari, Maciej Bilicki, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Benjamin Joachimi, Lance Miller, Jan Luca van den Busch, Angus H. Wright, Arun Kannawadi, Robert Reischke, HuanYuan Shan

Abstract: We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we inc… ▽ More We present refined cosmological parameter constraints derived from a cosmic shear analysis of the fourth data release of the Kilo-Degree Survey (KiDS-1000). Our main improvements include enhanced galaxy shape measurements made possible by an updated version of the lensfit code and improved shear calibration achieved with a newly developed suite of multi-band image simulations. Additionally, we incorporated recent advancements in cosmological inference from the joint Dark Energy Survey Year 3 and KiDS-1000 cosmic shear analysis. Assuming a spatially flat standard cosmological model, we constrain $S_8\equivσ_8(Ω_{\rm m}/0.3)^{0.5} = 0.776_{-0.027-0.003}^{+0.029+0.002}$, where the second set of uncertainties accounts for the systematic uncertainties within the shear calibration. These systematic uncertainties stem from minor deviations from realism in the image simulations and the sensitivity of the shear measurement algorithm to the morphology of the galaxy sample. Despite these changes, our results align with previous KiDS studies and other weak lensing surveys, and we find a ${\sim}2.3σ$ level of tension with the Planck cosmic microwave background constraints on $S_8$. △ Less

Submitted 2 November, 2023; v1 submitted 19 June, 2023; originally announced June 2023.

Comments: 20 pages, 13 figures, 4 tables, minor revisions to match the final accepted version

Journal ref: A&A 679, A133 (2023)

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

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

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

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

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

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

Wavelength-resolved reverberation mapping of intermediate redshift quasars HE 0413-4031 and HE 0435-4312: Dissecting Mg II, optical Fe II, and UV Fe II emission regions

Authors: Raj Prince, Michal Zajaček, S. Panda, K. Hryniewicz, V. K. Jaiswal, Bożena Czerny, P. Trzcionkowski, M. Bronikowski, M. Rałowski, C. S. Figaredo, M. L. Martinez-Aldama, M. Śniegowska, J. Średzińska, M. Bilicki, M-H Naddaf, A. Pandey, M. Haas, M. J. Sarna, G. Pietrzyński, V. Karas, A. Olejak, R. Przyłuski, R. R. Sefako, A. Genade, H. L. Worters , et al. (2 additional authors not shown)

Abstract: We present the wavelength-resolved reverberation mapping (RM) of combined MgII and UV FeII broad-line emissions for two intermediate redshifts (z$\sim$1), luminous quasars - HE 0413-4031 and HE 0435-4312, monitored by the SALT and 1-m class telescopes between 2012-2022. Through this technique, we aim to disentangle the Mg II and FeII emission regions and to build a radius-luminosity relation for U… ▽ More We present the wavelength-resolved reverberation mapping (RM) of combined MgII and UV FeII broad-line emissions for two intermediate redshifts (z$\sim$1), luminous quasars - HE 0413-4031 and HE 0435-4312, monitored by the SALT and 1-m class telescopes between 2012-2022. Through this technique, we aim to disentangle the Mg II and FeII emission regions and to build a radius-luminosity relation for UV FeII emission, which has so far remained unconstrained. Several methodologies have been applied to constrain the time delays for total MgII and FeII emissions. In addition, this technique is performed to quantify the inflow or outflow of broad-line region gas around the supermassive black hole and to disentangle the emission/emitting regions from lines produced in proximity to each other. The mean total FeII time delay is nearly equal to the mean total Mg II time delay for HE 0435-4312 suggesting a co-spatiality of their emissions. However, in HE 0413-4031, the mean FeII time delay is found to be longer than the mean MgII time delay, suggesting that FeII is produced at longer distances from the black hole. The UV Fe II R-L relation is updated with these two quasars and compared with the optical FeII relation, which suggests that the optical FeII region is located further than the UV FeII by a factor of 1.7-1.9 i.e. $R_{\rm FeII-opt}\sim(1.7-1.9)R_{\rm FeII-UV}$. We detected a weak pattern in the time delay vs. wavelength relation, suggesting that the MgII broad-line originates a bit closer to the SMBH than the UV FeII, however, the difference is not very significant. Comparison of MgII, UV, and optical FeII R-L relations suggests that the difference may be larger for lower-luminosity sources, possibly with the MgII emission originating further from the SMBH. In the future, more RM data will be acquired to put better constraints on these trends, in particular the UV FeII R-L relation. △ Less

Submitted 15 July, 2023; v1 submitted 26 April, 2023; originally announced April 2023.

Comments: 24 pages, 16 figures, 5 tables, Accepted for publication in A&A, in Press

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

ACT-DR5 Sunyaev-Zel'dovich Clusters: weak lensing mass calibration with KiDS

Authors: Naomi Clare Robertson, Cristóbal Sifón, Marika Asgari, Nicholas Battaglia, Maciej Bilicki, J. Richard Bond, Mark J. Devlin, Jo Dunkley, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Matt Hilton, Henk Hoekstra, John P. Hughes, Konrad Kuijken, Thibaut Louis, Maya Mallaby-Kay, Lyman Page, Bruce Partridge, Mario Radovich, Peter Schneider, HuanYuan Shan, David N. Spergel, Tilman Tröster, Edward J. Wollack , et al. (2 additional authors not shown)

Abstract: We present weak gravitational lensing measurements of a sample of 157 clusters within the Kilo Degree Survey (KiDS), detected with a $>5σ$ thermal Sunyaev-Zel'dovich (SZ) signal by the Atacama Cosmology Telescope (ACT). Using a halo-model approach we constrain the average total cluster mass, $M_{\rm WL}$, accounting for the ACT cluster selection function of the full sample. We find that the SZ clu… ▽ More We present weak gravitational lensing measurements of a sample of 157 clusters within the Kilo Degree Survey (KiDS), detected with a $>5σ$ thermal Sunyaev-Zel'dovich (SZ) signal by the Atacama Cosmology Telescope (ACT). Using a halo-model approach we constrain the average total cluster mass, $M_{\rm WL}$, accounting for the ACT cluster selection function of the full sample. We find that the SZ cluster mass estimate $M_{\rm SZ}$, which was calibrated using X-ray observations, is biased with $M_{\rm SZ}/M_{\rm WL} = (1-b_{\rm SZ}) = 0.65\pm 0.05$. Separating the sample into six mass bins, we find no evidence of a strong mass-dependency for the mass bias, $(1-b_{\rm SZ})$. Adopting this ACT-KiDS SZ mass-calibration would bring the Planck SZ cluster count into agreement with the counts expected from the {\it Planck} cosmic microwave background $Λ$CDM cosmological model, although it should be noted that the cluster sample considered in this work has a lower average mass $M_{\rm SZ, uncor} = 3.64 \times 10^{14} M_{\odot}$ compared to the Planck cluster sample which has an average mass in the range $M_{\rm SZ, uncor} = (5.5-8.5) \times 10^{14} M_{\odot}$, depending on the sub-sample used. △ Less

Submitted 20 April, 2023; originally announced April 2023.

Comments: 12 pages, 7 figures

Search for gravitational-lensing signatures in the full third observing run of the LIGO-Virgo network

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated… ▽ More Gravitational lensing by massive objects along the line of sight to the source causes distortions of gravitational wave-signals; such distortions may reveal information about fundamental physics, cosmology and astrophysics. In this work, we have extended the search for lensing signatures to all binary black hole events from the third observing run of the LIGO--Virgo network. We search for repeated signals from strong lensing by 1) performing targeted searches for subthreshold signals, 2) calculating the degree of overlap amongst the intrinsic parameters and sky location of pairs of signals, 3) comparing the similarities of the spectrograms amongst pairs of signals, and 4) performing dual-signal Bayesian analysis that takes into account selection effects and astrophysical knowledge. We also search for distortions to the gravitational waveform caused by 1) frequency-independent phase shifts in strongly lensed images, and 2) frequency-dependent modulation of the amplitude and phase due to point masses. None of these searches yields significant evidence for lensing. Finally, we use the non-detection of gravitational-wave lensing to constrain the lensing rate based on the latest merger-rate estimates and the fraction of dark matter composed of compact objects. △ Less

Submitted 17 April, 2023; originally announced April 2023.

Comments: 28 pages, 11 figures

Report number: LIGO-P2200031

Impact of modelling galaxy redshift uncertainties on the gravitational-wave dark standard siren measurement of the Hubble constant

Authors: Cezary Turski, Maciej Bilicki, Gergely Dálya, Rachel Gray, Archisman Ghosh

Abstract: Gravitational wave science is a new and rapidly expanding field of observational astronomy. Multimessenger observations of the binary neutron star merger GW170817 have provided some iconic results including the first gravitational-wave standard-siren measurement of the Hubble constant, opening up a new way to probe cosmology. The majority of the compact binary sources observed in gravitational wav… ▽ More Gravitational wave science is a new and rapidly expanding field of observational astronomy. Multimessenger observations of the binary neutron star merger GW170817 have provided some iconic results including the first gravitational-wave standard-siren measurement of the Hubble constant, opening up a new way to probe cosmology. The majority of the compact binary sources observed in gravitational waves are however without bright electromagnetic counterparts. In these cases, one can fall back on the ``dark standard siren'' approach to include information statistically from potential host galaxies. For such a measurement, we need to be cautious about all possible sources of systematic errors. In this paper, we begin to study the possible errors coming from the galaxy catalogue sector, and in particular, look into the effect of galaxy redshift uncertainties for the cases where these are photometry-based. We recalculate the dark standard siren Hubble constant using the latest GWTC-3 events and associated galaxy catalogues, with different galaxy redshift uncertainty models, namely, the standard Gaussian, a modified Lorentzian, and no uncertainty at all. We find that not using redshift uncertainties at all can lead to a potential bias comparable with other potential systematic effects previously considered for the GWTC-3 $H_0$ measurement (however still small compared to the overall statistical error in this measurement). The difference between different uncertainty models leads to small differences in the results for the current data; their impact is much smaller than the current statistical errors and other potential sources of systematic errors which have been considered in previous robustness studies. △ Less

Submitted 10 October, 2023; v1 submitted 23 February, 2023; originally announced February 2023.

Comments: 9 pages, 8 figures, accepted to publication in MNRAS

KiDS-1000: cross-correlation with Planck cosmic microwave background lensing and intrinsic alignment removal with self-calibration

Authors: Ji Yao, Huanyuan Shan, Pengjie Zhang, Xiangkun Liu, Catherine Heymans, Benjamin Joachimi, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Konrad Kuijken, Tilman Tröster, Jan Luca van den Busch, Angus Wright, Ziang Yan

Abstract: Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared… ▽ More Galaxy shear - cosmic microwave background (CMB) lensing convergence cross-correlations contain additional information on cosmology to auto-correlations. While being immune to certain systematic effects, they are affected by the galaxy intrinsic alignments (IA). This may be responsible for the reported low lensing amplitude of the galaxy shear $\times$ CMB convergence cross-correlations, compared to the standard Planck $Λ$CDM (cosmological constant and cold dark matter) cosmology prediction. In this work, we investigate how IA affects the Kilo-Degree Survey (KiDS) galaxy lensing shear - Planck CMB lensing convergence cross-correlation and compare it to previous treatments with or without IA taken into consideration. More specifically, we compare marginalization over IA parameters and the IA self-calibration (SC) method (with additional observables defined only from the source galaxies) and prove that SC can efficiently break the degeneracy between the CMB lensing amplitude $A_{\rm lens}$ and the IA amplitude $A_{\rm IA}$. We further investigate how different systematics affect the resulting $A_{\rm IA}$ and $A_{\rm lens}$, and validate our results with the MICE2 simulation. We find that by including the SC method to constrain IA, the information loss due to the degeneracy between CMB lensing and IA is strongly reduced. The best-fit values are $A_{\rm lens}=0.84^{+0.22}_{-0.22}$ and $A_{\rm IA}=0.60^{+1.03}_{-1.03}$, while different angular scale cuts can affect $A_{\rm lens}$ by $\sim10\%$. We show that appropriate treatment of the boost factor, cosmic magnification, and photometric redshift modeling is important for obtaining the correct IA and cosmological results. △ Less

Submitted 10 April, 2023; v1 submitted 31 January, 2023; originally announced January 2023.

Comments: match version accepted by A&A

Journal ref: A&A 673, A111 (2023)

Improved analytical modeling of the non-linear power spectrum in modified gravity cosmologies

Authors: Suhani Gupta, Wojciech A. Hellwing, Maciej Bilicki

Abstract: Reliable analytical modeling of the non-linear power spectrum (PS) of matter perturbations is among the chief pre-requisites for cosmological analyses from the largest sky surveys. This is especially true for the models that extend the standard general-relativity paradigm by adding the fifth force, where numerical simulations can be prohibitively expensive. Here we present a method for building ac… ▽ More Reliable analytical modeling of the non-linear power spectrum (PS) of matter perturbations is among the chief pre-requisites for cosmological analyses from the largest sky surveys. This is especially true for the models that extend the standard general-relativity paradigm by adding the fifth force, where numerical simulations can be prohibitively expensive. Here we present a method for building accurate PS models for two modified gravity (MG) variants: namely the Hu-Sawicki $f(R)$, and the normal branch of the Dvali-Gabadadze-Porrati (nDGP) braneworld. We start by modifying the standard halo model (HM) with respect to the baseline Lambda-Cold-Dark-Matter ($Λ$CDM) scenario, by using the HM components with specific MG extensions. We find that our $P(k)_{\text{HM}}$ retains 5% accuracy only up to mildly non-linear scales ($k \lesssim 0.3$ $h/\,\mbox{Mpc}$) when compared to PS from numerical simulations. At the same time, our HM prescription much more accurately captures the ratio $Υ(k) = P(k)_{\text{MG}}/P(k)_{Λ\text{CDM}}$ up to non-linear scales. We show that using HM-derived $Υ(k)$ together with a viable non-linear $Λ$CDM $P(k)$ prescription (such as HALOFIT), we render a much better and more accurate PS predictions in MG. The new approach yields considerably improved performance, with modeled $P(k)_{\text{MG}}$ being now accurate to within 5% all the way to non-linear scales of $k \lesssim 2.5-3$ $h/\,\mbox{Mpc}$. The magnitude of deviations from GR as fostered by these MG models is typically $\mathcal{O}(10\%)$ in these regimes. Therefore reaching 5% PS modeling is enough for forecasting constraints on modern-era cosmological observables. △ Less

Submitted 27 January, 2023; originally announced January 2023.

The Faraday Rotation Measure Grid of the LOFAR Two-metre Sky Survey: Data Release 2

Authors: S. P. O'Sullivan, T. W. Shimwell, M. J. Hardcastle, C. Tasse, G. Heald, E. Carretti, M. Brüggen, V. Vacca, C. Sobey, C. L. Van Eck, C. Horellou, R. Beck, M. Bilicki, S. Bourke, A. Botteon, J. H. Croston, A. Drabent, K. Duncan, V. Heesen, S. Ideguchi, M. Kirwan, L. Lawlor, B. Mingo, B. Nikiel-Wroczyński, J. Piotrowska , et al. (2 additional authors not shown)

Abstract: A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polar… ▽ More A Faraday rotation measure (RM) catalogue, or RM Grid, is a valuable resource for the study of cosmic magnetism. Using the second data release (DR2) from the LOFAR Two-metre Sky Survey (LoTSS), we have produced a catalogue of 2461 extragalactic high-precision RM values across 5720 deg$^{2}$ of sky (corresponding to a polarized source areal number density of $\sim$0.43 deg$^{-2}$). The linear polarization and RM properties were derived using RM synthesis from the Stokes $Q$ and $U$ channel images at an angular resolution of 20'' across a frequency range of 120 to 168 MHz with a channel bandwidth of 97.6 kHz. The fraction of total intensity sources ($>1$ mJy beam$^{-1}$) found to be polarized was $\sim$0.2%. The median detection threshold was 0.6 mJy beam$^{-1}$ ($8σ_{QU}$), with a median RM uncertainty of 0.06 rad m$^{-2}$ (although a systematic uncertainty of up to 0.3 rad m$^{-2}$ is possible, after the ionosphere RM correction). The median degree of polarization of the detected sources is 1.8%, with a range of 0.05% to 31%. Comparisons with cm-wavelength RMs indicate minimal amounts of Faraday complexity in the LoTSS detections, making them ideal sources for RM Grid studies. Host galaxy identifications were obtained for 88% of the sources, along with redshifts for 79% (both photometric and spectroscopic), with the median redshift being 0.6. The focus of the current catalogue was on reliability rather than completeness, and we expect future versions of the LoTSS RM Grid to have a higher areal number density. In addition, 25 pulsars were identified, mainly through their high degrees of linear polarization. △ Less

Submitted 18 January, 2023; originally announced January 2023.

Comments: 20 pages, 23 figures, 3 tables. Accepted for publication in MNRAS. Go to https://lofar-mksp.org/data/ to download the catalogue and associated data products

Modeling Strong Lenses from Wide-Field Ground-Based Observations in KiDS and GAMA

Authors: Shawn Knabel, B. W. Holwerda, J. Nightingale, T. Treu, M. Bilicki, S. Brough, S. Driver, L. Finnerty, L. Haberzettl, S. Hegde, A. M. Hopkins, K. Kuijken, J. Liske, K. A. Pimbblet, R. C. Steele, A. H. Wright

Abstract: Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based i… ▽ More Despite the success of galaxy-scale strong gravitational lens studies with Hubble-quality imaging, the number of well-studied strong lenses remains small. As a result, robust comparisons of the lens models to theoretical predictions are difficult. This motivates our application of automated Bayesian lens modeling methods to observations from public data releases of overlapping large ground-based imaging and spectroscopic surveys: Kilo-Degree Survey (KiDS) and Galaxy and Mass Assembly (GAMA), respectively. We use the open-source lens modeling software PyAutoLens to perform our analysis. We demonstrate the feasibility of strong lens modeling with large-survey data at lower resolution as a complementary avenue to studies that utilize more time-consuming and expensive observations of individual lenses at higher resolution. We discuss advantages and challenges, with special consideration given to determining background source redshifts from single-aperture spectra and to disentangling foreground lens and background source light. High uncertainties in the best-fit parameters for the models due to the limits of optical resolution in ground-based observatories and the small sample size can be improved with future study. We give broadly applicable recommendations for future efforts, and with proper application this approach could yield measurements in the quantities needed for robust statistical inference. △ Less

Submitted 17 January, 2023; v1 submitted 12 January, 2023; originally announced January 2023.

Comments: 25 pages, 19 figures, accepted for publication in MNRAS following peer review

The Hitchhiker's guide to the galaxy catalog approach for gravitational wave cosmology

Authors: Jonathan R. Gair, Archisman Ghosh, Rachel Gray, Daniel E. Holz, Simone Mastrogiovanni, Suvodip Mukherjee, Antonella Palmese, Nicola Tamanini, Tessa Baker, Freija Beirnaert, Maciej Bilicki, Hsin-Yu Chen, Gergely Dálya, Jose Maria Ezquiaga, Will M. Farr, Maya Fishbach, Juan Garcia-Bellido, Tathagata Ghosh, Hsiang-Yu Huang, Christos Karathanasis, Konstantin Leyde, Ignacio Magaña Hernandez, Johannes Noller, Gregoire Pierra, Peter Raffai , et al. (6 additional authors not shown)

Abstract: We outline the ``dark siren'' galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz (… ▽ More We outline the ``dark siren'' galaxy catalog method for cosmological inference using gravitational wave (GW) standard sirens, clarifying some common misconceptions in the implementation of this method. When a confident transient electromagnetic counterpart to a GW event is unavailable, the identification of a unique host galaxy is in general challenging. Instead, as originally proposed by Schutz (1986), one can consult a galaxy catalog and implement a dark siren statistical approach incorporating all potential host galaxies within the localization volume. Trott & Hunterer 2021 recently claimed that this approach results in a biased estimate of the Hubble constant, $H_0$, when implemented on mock data, even if optimistic assumptions are made. We demonstrate explicitly that, as previously shown by multiple independent groups, the dark siren statistical method leads to an unbiased posterior when the method is applied to the data correctly. We highlight common sources of error possible to make in the generation of mock data and implementation of the statistical framework, including the mismodeling of selection effects and inconsistent implementations of the Bayesian framework, which can lead to a spurious bias. △ Less

Submitted 16 December, 2022; originally announced December 2022.

Comments: 19 pages, 11 figures. Submitted to ApJ

Search for subsolar-mass black hole binaries in the second part of Advanced LIGO's and Advanced Virgo's third observing run

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1680 additional authors not shown)

Abstract: We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate t… ▽ More We describe a search for gravitational waves from compact binaries with at least one component with mass 0.2 $M_\odot$ -- $1.0 M_\odot$ and mass ratio $q \geq 0.1$ in Advanced LIGO and Advanced Virgo data collected between 1 November 2019, 15:00 UTC and 27 March 2020, 17:00 UTC. No signals were detected. The most significant candidate has a false alarm rate of 0.2 $\mathrm{yr}^{-1}$. We estimate the sensitivity of our search over the entirety of Advanced LIGO's and Advanced Virgo's third observing run, and present the most stringent limits to date on the merger rate of binary black holes with at least one subsolar-mass component. We use the upper limits to constrain two fiducial scenarios that could produce subsolar-mass black holes: primordial black holes (PBH) and a model of dissipative dark matter. The PBH model uses recent prescriptions for the merger rate of PBH binaries that include a rate suppression factor to effectively account for PBH early binary disruptions. If the PBHs are monochromatically distributed, we can exclude a dark matter fraction in PBHs $f_\mathrm{PBH} \gtrsim 0.6$ (at 90% confidence) in the probed subsolar-mass range. However, if we allow for broad PBH mass distributions we are unable to rule out $f_\mathrm{PBH} = 1$. For the dissipative model, where the dark matter has chemistry that allows a small fraction to cool and collapse into black holes, we find an upper bound $f_{\mathrm{DBH}} < 10^{-5}$ on the fraction of atomic dark matter collapsed into black holes. △ Less

Submitted 26 January, 2024; v1 submitted 2 December, 2022; originally announced December 2022.

Comments: https://dcc.ligo.org/P2200139

KiDS-Legacy calibration: unifying shear and redshift calibration with the SKiLLS multi-band image simulations

Authors: Shun-Sheng Li, Konrad Kuijken, Henk Hoekstra, Lance Miller, Catherine Heymans, Hendrik Hildebrandt, Jan Luca van den Busch, Angus H. Wright, Mijin Yoon, Maciej Bilicki, Matías Bravo, Claudia del P. Lagos

Abstract: We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we int… ▽ More We present SKiLLS, a suite of multi-band image simulations for the weak lensing analysis of the complete Kilo-Degree Survey (KiDS), dubbed KiDS-Legacy analysis. The resulting catalogues enable joint shear and redshift calibration, enhancing the realism and hence accuracy over previous efforts. To create a large volume of simulated galaxies with faithful properties and to a sufficient depth, we integrated cosmological simulations with high-quality imaging observations. We also improved the realism of simulated images by allowing the point spread function (PSF) to differ between CCD images, including stellar density variations and varying noise levels between pointings. Using realistic variable shear fields, we accounted for the impact of blended systems at different redshifts. Although the overall correction is minor, we found a clear redshift-bias correlation in the blending-only variable shear simulations, indicating the non-trivial impact of this higher-order blending effect. We also explored the impact of the PSF modelling errors and found a small yet noticeable effect on the shear bias. Finally, we conducted a series of sensitivity tests, including changing the input galaxy properties. We conclude that our fiducial shape measurement algorithm, lensfit, is robust within the requirements of lensing analyses with KiDS. As for future weak lensing surveys with tighter requirements, we suggest further investments in understanding the impact of blends at different redshifts, improving the PSF modelling algorithm and developing the shape measurement method to be less sensitive to the galaxy properties. △ Less

Submitted 22 January, 2023; v1 submitted 13 October, 2022; originally announced October 2022.

Comments: 28 pages, 31 figures, 2 tables, minor revisions to match the final accepted version

Journal ref: A&A 670, A100 (2023)

KiDS-1000: Combined halo-model cosmology constraints from galaxy abundance, galaxy clustering and galaxy-galaxy lensing

Authors: Andrej Dvornik, Catherine Heymans, Marika Asgari, Constance Mahony, Benjamin Joachimi, Maciej Bilicki, Elisa Chisari, Hendrik Hildebrandt, Henk Hoekstra, Harry Johnston, Konrad Kuijken, Alexander Mead, Hironao Miyatake, Takahiro Nishimichi, Robert Reischke, Sandra Unruh, Angus H. Wright

Abstract: We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matt… ▽ More We present constraints on the flat $Λ$CDM cosmological model through a joint analysis of galaxy abundance, galaxy clustering and galaxy-galaxy lensing observables with the Kilo-Degree Survey. Our theoretical model combines a flexible conditional stellar mass function, to describe the galaxy-halo connection, with a cosmological N-body simulation-calibrated halo model to describe the non-linear matter field. Our magnitude-limited bright galaxy sample combines 9-band optical-to-near-infrared photometry with an extensive and complete spectroscopic training sample to provide accurate redshift and stellar mass estimates. Our faint galaxy sample provides a background of accurately calibrated lensing measurements. We constrain the structure growth parameter $S_8=σ_8\sqrt{Ω_{\mathrm{m}}/0.3}=0.773^{+0.028}_{-0.030}$, and the matter density parameter $Ω_{\mathrm{m}}=0.290^{+0.021}_{-0.017}$. The galaxy-halo connection model adopted in the work is shown to be in agreement with previous studies. Our constraints on cosmological parameters are comparable to, and consistent with, joint $3\times2{\mathrm{pt}}$ clustering-lensing analyses that additionally include a cosmic shear observable. This analysis therefore brings attention to the significant constraining power in the often-excluded non-linear scales for galaxy clustering and galaxy-galaxy lensing observables. By adopting a theoretical model that accounts for non-linear halo bias, halo exclusion, scale-dependent galaxy bias and the impact of baryon feedback, this work demonstrates the potential and a way forward to include non-linear scales in cosmological analyses. Varying the width of the satellite galaxy distribution with an additional parameter yields a strong preference for sub-Poissonian variance, improving the goodness of fit by 0.18 in reduced $χ^{2}$ value compared to a fixed Poisson distribution. △ Less

Submitted 13 June, 2023; v1 submitted 6 October, 2022; originally announced October 2022.

Comments: 25 pages, 16 figures, accepted for publication in A&A

Journal ref: A&A 675, A189 (2023)

Constrained simulations of the local Universe with Modified Gravity

Authors: Krishna Naidoo, Wojciech A. Hellwing, Maciej Bilicki, Noam Libeskind, Simon Pfeifer, Yehuda Hoffman

Abstract: We present a methodology for constructing modified gravity (MG) constrained simulations of the local Universe using positions and peculiar velocities from the CosmicFlows data set. Our analysis focuses on the following MG models: the normal branch of the Dvali-Gabadadze-Porrati (nDGP) model and Hu-Sawicki $f(R)$ model. We develop a model independent methodology for constructing constrained simulat… ▽ More We present a methodology for constructing modified gravity (MG) constrained simulations of the local Universe using positions and peculiar velocities from the CosmicFlows data set. Our analysis focuses on the following MG models: the normal branch of the Dvali-Gabadadze-Porrati (nDGP) model and Hu-Sawicki $f(R)$ model. We develop a model independent methodology for constructing constrained simulations with any given power spectra and numerically calculated linear growth functions. Initial conditions (ICs) for a set of constrained simulations are constructed for the standard cosmological model $Λ$CDM and the MG models. Differences between the model's reconstructed Wiener filtered density and the resultant simulation density are presented showing the importance for the generation of MG constrained ICs to study the subtle effects of MG in the local Universe. These are the first MG constrained simulations ever produced. The current work paves the way to improved approximate methods for models with scale-dependent growth functions, such as $f(R)$, and for high-resolution hydrodynamical MG zoom-in simulations of the local Universe. △ Less

Submitted 24 February, 2023; v1 submitted 28 September, 2022; originally announced September 2022.

Comments: 10 pages (+references), 6 figures, accepted in Phys. Rev. D

Dark matter solution to the $H_{0}$ and $S_{8}$ tensions, and the integrated Sachs-Wolfe void anomaly

Authors: Krishna Naidoo, Mariana Jaber, Wojciech A. Hellwing, Maciej Bilicki

Abstract: We consider a phenomenological model of dark matter with an equation-of-state $w$ that is negative and changing at late times. We show this couples the $H_{0}$ and $S_{8}$ tensions, alleviating them both simultaneously, reducing the $H_{0}$ tension from $\sim5σ$ to $\sim3σ$ and the $S_{8}$ tension from $\sim3σ$ to $\sim1σ$. Furthermore, the model provides an explanation for the anomalously large i… ▽ More We consider a phenomenological model of dark matter with an equation-of-state $w$ that is negative and changing at late times. We show this couples the $H_{0}$ and $S_{8}$ tensions, alleviating them both simultaneously, reducing the $H_{0}$ tension from $\sim5σ$ to $\sim3σ$ and the $S_{8}$ tension from $\sim3σ$ to $\sim1σ$. Furthermore, the model provides an explanation for the anomalously large integrated Sachs-Wolfe (ISW) effect from cosmic voids, a unique consequence of the changing and negative equation-of-state. Observations of high ISW from cosmic voids may therefore be evidence that dark matter plays a significant role in both the $H_{0}$ and $S_{8}$ tensions. We predict the ISW from cosmic voids to be a factor of up to $\sim2$ greater in this model than what is expected from the standard model $Λ$CDM. These results extend to other degenerate models of dark matter, such as unified or interacting dark matter and dark energy models. △ Less

Submitted 13 April, 2024; v1 submitted 16 September, 2022; originally announced September 2022.

Comments: 9 pages (+references), 7 figures, published in Physical Review D

Model-based cross-correlation search for gravitational waves from the low-mass X-ray binary Scorpius X-1 in LIGO O3 data

Authors: The LIGO Scientific Collaboration, the Virgo Collaboration, the KAGRA Collaboration, R. Abbott, H. Abe, F. Acernese, K. Ackley, S. Adhicary, N. Adhikari, R. X. Adhikari, V. K. Adkins, V. B. Adya, C. Affeldt, D. Agarwal, M. Agathos, O. D. Aguiar, L. Aiello, A. Ain, P. Ajith, T. Akutsu, S. Albanesi, R. A. Alfaidi, C. Alléné, A. Allocca, P. A. Altin , et al. (1670 additional authors not shown)

Abstract: We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to bala… ▽ More We present the results of a model-based search for continuous gravitational waves from the low-mass X-ray binary Scorpius X-1 using LIGO detector data from the third observing run of Advanced LIGO, Advanced Virgo and KAGRA. This is a semicoherent search which uses details of the signal model to coherently combine data separated by less than a specified coherence time, which can be adjusted to balance sensitivity with computing cost. The search covered a range of gravitational-wave frequencies from 25Hz to 1600Hz, as well as ranges in orbital speed, frequency and phase determined from observational constraints. No significant detection candidates were found, and upper limits were set as a function of frequency. The most stringent limits, between 100Hz and 200Hz, correspond to an amplitude h0 of about 1e-25 when marginalized isotropically over the unknown inclination angle of the neutron star's rotation axis, or less than 4e-26 assuming the optimal orientation. The sensitivity of this search is now probing amplitudes predicted by models of torque balance equilibrium. For the usual conservative model assuming accretion at the surface of the neutron star, our isotropically-marginalized upper limits are close to the predicted amplitude from about 70Hz to 100Hz; the limits assuming the neutron star spin is aligned with the most likely orbital angular momentum are below the conservative torque balance predictions from 40Hz to 200Hz. Assuming a broader range of accretion models, our direct limits on gravitational-wave amplitude delve into the relevant parameter space over a wide range of frequencies, to 500Hz or more. △ Less

Submitted 2 January, 2023; v1 submitted 6 September, 2022; originally announced September 2022.

Comments: 19 pages, Open Access Journal PDF

Report number: LIGO-P2100110-v13

Journal ref: The Astrophysical Journal Letters, 941, L30 (2022)

KiDS-1000 Cosmology: Constraints from density split statistics

Authors: Pierre A. Burger, Oliver Friedrich, Joachim Harnois-Déraps, Peter Schneider, Marika Asgari, Maciej Bilicki, Hendrik Hildebrandt, Angus H. Wright, Tiago Castro, Klaus Dolag, Catherine Heymans, Benjamin Joachimi, Konrad Kuijken, Nicolas Martinet, HuanYuan Shan, Tilman Tröster

Abstract: Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based o… ▽ More Context. Weak lensing and clustering statistics beyond two-point functions can capture non-Gaussian information about the matter density field, thereby improving the constraints on cosmological parameters relative to the mainstream methods based on correlation functions and power spectra. Aims. This paper presents a cosmological analysis of the fourth data release of the Kilo Degree Survey based on the density split statistics, which measures the mean shear profiles around regions classified according to foreground densities. The latter is constructed from a bright galaxy sample, which we further split into red and blue samples, allowing us to probe their respective connection to the underlying dark matter density. Methods. We use the state-of-the-art model of the density splitting statistics and validate its robustness against mock data infused with known systematic effects such as intrinsic galaxy alignment and baryonic feedback. Results. After marginalising over the photometric redshift uncertainty and the residual shear calibration bias, we measure for the full KiDS-bright sample a structure growth parameter of $S_8 = σ_8 \sqrt{Ω_\mathrm{m}/0.3} = 0.74^{+0.03}_{-0.02}$ that is competitive to and consistent with two-point cosmic shear results, a matter density of $Ω_\mathrm{m} = 0.28 \pm 0.02$, and a constant galaxy bias of $b = 1.32^{+0.12}_{-0.10}$. △ Less

Submitted 10 January, 2023; v1 submitted 3 August, 2022; originally announced August 2022.

Comments: 18 pages, 23 Figures

Report number: A69

Journal ref: A&A 669, A69 (2023)

Galaxy and Mass Assembly (GAMA): Probing galaxy-group correlations in redshift space with the halo streaming model

Authors: Qianjun Hang, John A. Peacock, Shadab Alam, Yan-Chuan Cai, Katarina Kraljic, Marcel van Daalen, M. Bilicki, B. W. Holwerda, J. Loveday

Abstract: We have studied the galaxy-group cross-correlations in redshift space for the Galaxy And Mass Assembly (GAMA) Survey. We use a set of mock GAMA galaxy and group catalogues to develop and test a novel 'halo streaming' model for redshift-space distortions. This treats 2-halo correlations via the streaming model, plus an empirical 1-halo term derived from the mocks, allowing accurate modelling into t… ▽ More We have studied the galaxy-group cross-correlations in redshift space for the Galaxy And Mass Assembly (GAMA) Survey. We use a set of mock GAMA galaxy and group catalogues to develop and test a novel 'halo streaming' model for redshift-space distortions. This treats 2-halo correlations via the streaming model, plus an empirical 1-halo term derived from the mocks, allowing accurate modelling into the nonlinear regime. In order to probe the robustness of the growth rate inferred from redshift-space distortions, we divide galaxies by colour, and divide groups according to their total stellar mass, calibrated to total mass via gravitational lensing. We fit our model to correlation data, to obtain estimates of the perturbation growth rate, $fσ_8$, validating parameter errors via the dispersion between different mock realizations. In both mocks and real data, we demonstrate that the results are closely consistent between different subsets of the group and galaxy populations, considering the use of correlation data down to some minimum projected radius, $r_{\rm min}$. For the mock data, we can use the halo streaming model to below $r_{\rm min} = 5h^{-1}$ Mpc, finding that all subsets yield growth rates within about 3% of each other, and consistent with the true value. For the actual GAMA data, the results are limited by cosmic variance: $fσ_8=0.29\pm 0.10$ at an effective redshift of 0.20; but there is every reason to expect that this method will yield precise constraints from larger datasets of the same type, such as the DESI bright galaxy survey. △ Less

Submitted 7 September, 2022; v1 submitted 10 June, 2022; originally announced June 2022.

Comments: 17 pages, 12 figures. Matched to the accepted version

Galaxy morphoto-Z with neural Networks (GaZNets). I. Optimized accuracy and outlier fraction from Imaging and Photometry

Authors: Rui Li, Nicola R. Napolitano, Haicheng Feng, Ran Li, Valeria Amaro, Linghua Xie, Crescenzo Tortora, Maciej Bilicki, Massimo Brescia, Stefano Cavuoti, Mario Radovich

Abstract: In the era of large sky surveys, photometric redshifts (photo-z) represent crucial information for galaxy evolution and cosmology studies. In this work, we propose a new Machine Learning (ML) tool called Galaxy morphoto-Z with neural Networks (GaZNet-1), which uses both images and multi-band photometry measurements to predict galaxy redshifts, with accuracy, precision and outlier fraction superior… ▽ More In the era of large sky surveys, photometric redshifts (photo-z) represent crucial information for galaxy evolution and cosmology studies. In this work, we propose a new Machine Learning (ML) tool called Galaxy morphoto-Z with neural Networks (GaZNet-1), which uses both images and multi-band photometry measurements to predict galaxy redshifts, with accuracy, precision and outlier fraction superior to standard methods based on photometry only. As a first application of this tool, we estimate photo-z of a sample of galaxies in the Kilo-Degree Survey (KiDS). GaZNet-1 is trained and tested on $\sim140 000$ galaxies collected from KiDS Data Release 4 (DR4), for which spectroscopic redshifts are available from different surveys. This sample is dominated by bright (MAG$\_$AUTO$<21$) and low redshift ($z < 0.8$) systems, however, we could use $\sim$ 6500 galaxies in the range $0.8 < z < 3$ to effectively extend the training to higher redshift. The inputs are the r-band galaxy images plus the 9-band magnitudes and colours, from the combined catalogs of optical photometry from KiDS and near-infrared photometry from the VISTA Kilo-degree Infrared survey. By combining the images and catalogs, GaZNet-1 can achieve extremely high precision in normalized median absolute deviation (NMAD=0.014 for lower redshift and NMAD=0.041 for higher redshift galaxies) and low fraction of outliers ($0.4$\% for lower and $1.27$\% for higher redshift galaxies). Compared to ML codes using only photometry as input, GaZNet-1 also shows a $\sim 10-35$% improvement in precision at different redshifts and a $\sim$ 45% reduction in the fraction of outliers. We finally discuss that, by correctly separating galaxies from stars and active galactic nuclei, the overall photo-z outlier fraction of galaxies can be cut down to $0.3$\%. △ Less

Submitted 19 July, 2022; v1 submitted 21 May, 2022; originally announced May 2022.

Comments: Accepted for publication in A&A

Journal ref: A&A 666, A85 (2022)

Galaxy and halo angular clustering in LCDM and Modified Gravity cosmologies

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

Abstract: Using a suite of $N$-body simulations we study the angular clustering of galaxies, halos, and dark matter in $\mathrm{Λ\text{CDM}}$ and Modified Gravity (MG) scenarios. We consider two general categories of such MG models, one is the $f(R)$ gravity, and the other is the normal branch of the Dvali-Gabadadze-Porrati brane world (nDGP). To measure angular clustering we construct a set of observer-fra… ▽ More Using a suite of $N$-body simulations we study the angular clustering of galaxies, halos, and dark matter in $\mathrm{Λ\text{CDM}}$ and Modified Gravity (MG) scenarios. We consider two general categories of such MG models, one is the $f(R)$ gravity, and the other is the normal branch of the Dvali-Gabadadze-Porrati brane world (nDGP). To measure angular clustering we construct a set of observer-frame lightcones and resulting mock sky catalogs. We focus on the area-averaged angular correlation functions, $W_J$, and the associated reduced cumulants, $S_J\equiv W_J/W_2^{(J-1)}$, and robustly measure them up to the 9th order using counts-in-cells (CIC). We find that $0.15 < z < 0.3$ is the optimal redshift range to maximize the MG signal in our lightcones. Analyzing various scales for the two types of statistics, we identify up to 20\% relative departures in MG measurements from general relativity (GR), with varying signal significance. For the case of halos and galaxies, we find that $3$rd order statistics offer the most sensitive probe of the different structure formation scenarios, with both $W_3$ and the reduced skewness, $S_3$, reaching from $2σ$ to $4σ$ significance at angular scales $θ\sim 0.13 ^\circ$. The MG clustering of the smooth dark matter field is characterized by even stronger deviations ($\stackrel{>}{{}_\sim} 5σ$) from GR, albeit at a bit smaller scales of $θ\sim0.08^\circ$, where baryonic physics is already important. Finally, we stress out that our mock halo and galaxy catalogs are characterized by rather low surface number densities when compared to existing and forthcoming state-of-the-art photometric surveys. This opens up exciting potential for testing GR and MG using angular clustering in future applications, with even higher precision and significance than reported here. △ Less

Submitted 8 August, 2022; v1 submitted 11 April, 2022; originally announced April 2022.

Comments: 15 pages, 8 figures, 1 table

Journal ref: Phys. Rev. D 106, 043513 (2022)

KiDS-1000: cosmic shear with enhanced redshift calibration

Authors: J. L. van den Busch, A. H. Wright, H. Hildebrandt, M. Bilicki, M. Asgari, S. Joudaki, C. Blake, C. Heymans, A. Kannawadi, H. Y. Shan, T. Tröster

Abstract: We present a cosmic shear analysis with an improved redshift calibration for the fourth data release of the Kilo-Degree Survey (KiDS-1000) using self-organising maps (SOMs). Compared to the previous analysis of the KiDS-1000 data, we expand the redshift calibration sample to more than twice its size, now consisting of data of 17 spectroscopic redshift campaigns, and significantly extending the fra… ▽ More We present a cosmic shear analysis with an improved redshift calibration for the fourth data release of the Kilo-Degree Survey (KiDS-1000) using self-organising maps (SOMs). Compared to the previous analysis of the KiDS-1000 data, we expand the redshift calibration sample to more than twice its size, now consisting of data of 17 spectroscopic redshift campaigns, and significantly extending the fraction of KiDS galaxies we are able to calibrate with our SOM redshift methodology. We then enhance the calibration sample with precision photometric redshifts from COSMOS2015 and the Physics of the Accelerated Universe Survey (PAUS), allowing us to fill gaps in the spectroscopic coverage of the KiDS data. Finally we perform a Complete Orthogonal Sets of E/B-Integrals (COSEBIs) cosmic shear analysis of the newly calibrated KiDS sample. We find $S_8 = 0.748_{-0.025}^{+0.021}$, which is in good agreement with previous KiDS studies and increases the tension with measurements of the cosmic microwave background to 3.4σ. We repeat the redshift calibration with different subsets of the full calibration sample and obtain, in all cases, agreement within at most 0.5σ in $S_8$ compared to our fiducial analysis. Including additional photometric redshifts allows us to calibrate an additional 6 % of the source galaxy sample. Even though further systematic testing with simulated data is necessary to quantify the impact of redshift outliers, precision photometric redshifts can be beneficial at high redshifts and to mitigate selection effects commonly found in spectroscopically selected calibration samples. △ Less

Submitted 5 April, 2022; originally announced April 2022.

Comments: 18 pages, 15 figures, 6 tables

Journal ref: A&A 664, A170 (2022)

Cosmic star formation history with tomographic cosmic infrared background-galaxy cross-correlation

Authors: Ziang Yan, Ludovic van Waerbeke, Angus H. Wright, Maciej Bilicki, Shiming Gu, Hendrik Hildebrandt, Abhishek S. Maniyar, Tilman Tröster

Abstract: In this work, we probe the star formation history of the Universe using tomographic cross-correlation between the cosmic infrared background (CIB) and galaxy samples. The galaxy samples are from the Kilo-Degree Survey (KiDS), while the CIB maps are made from \planck\, sky maps. We measure the cross-correlation in harmonic space with a significance of 43$σ$. We model the cross-correlation with a ha… ▽ More In this work, we probe the star formation history of the Universe using tomographic cross-correlation between the cosmic infrared background (CIB) and galaxy samples. The galaxy samples are from the Kilo-Degree Survey (KiDS), while the CIB maps are made from \planck\, sky maps. We measure the cross-correlation in harmonic space with a significance of 43$σ$. We model the cross-correlation with a halo model, which links CIB anisotropies to star formation rates (SFR) and galaxy abundance. We assume that SFR has a lognormal dependence on halo mass, while galaxy abundance follows the halo occupation distribution (HOD) model. The cross-correlations give a best-fit maximum star formation efficiency of $η_{\mathrm{max}}= 0.41^{+0.09}_{-0.14}$ at a halo mass $\log_{10}(M_{\mathrm{peak}}/M_{\odot})= {12.14\pm 0.36}$. The derived star formation rate density (SFRD) is well constrained up to $z\sim 1.5$. The constraining power at high redshift is mainly limited by the KiDS survey depth. A combination with external SFRD measurements from previous studies gives $\log_{10}(M_{\mathrm{peak}}/M_{\odot})=12.42^{+0.35}_{-0.19}$. This tightens the SFRD constraint up to $z=4$, yielding a peak SFRD of $0.09_{-0.004}^{+0.003}\,M_{\odot} \mathrm { year }^{-1} \mathrm{Mpc}^{-3}$ at $z=1.74^{+0.06}_{-0.02}$, corresponding to a lookback time of $10.05^{+0.12}_{-0.03}$ Gyr. Both constraints are consistent, and the derived SFRD agrees with previous studies and simulations. Additionally, we estimate the galaxy bias $b$ of KiDS galaxies from the constrained HOD parameters and yield an increasing bias from $b=1.1_{-0.31}^{+0.17}$ at $z=0$ to $b=1.96_{-0.64}^{+0.18}$ at $z=1.5$. Finally, we provide a forecast for future galaxy surveys and conclude that, due to their considerable depth, future surveys will yield a much tighter constraint on the evolution of the SFRD. △ Less

Submitted 18 November, 2022; v1 submitted 4 April, 2022; originally announced April 2022.

Comments: 22 pages, 14 figures, 3 tables, the abstract is abridged, published on A&A

Journal ref: A&A 665, A52 (2022)

Galaxy And Mass Assembly (GAMA): Self-Organizing Map Application on Nearby Galaxies

Authors: B. W. Holwerda, Dominic Smith, Lori Porter, Chris Henry, Ren Porter-Temple, Kyle Cook, Kevin A. Pimbblet, Andrew M. Hopkins, Maciej Bilicki, Sebastian Turner, Viviana Acquaviva, Lingyu Wang, Angus H. Wright, Lee S. Kelvin, Meiert W. Grootes

Abstract: Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a mor… ▽ More Galaxy populations show bimodality in a variety of properties: stellar mass, colour, specific star-formation rate, size, and Sérsic index. These parameters are our feature space. We use an existing sample of 7556 galaxies from the Galaxy and Mass Assembly (GAMA) survey, represented using five features and the K-means clustering technique, showed that the bimodalities are the manifestation of a more complex population structure, represented by between 2 and 6 clusters. Here we use Self Organizing Maps (SOM), an unsupervised learning technique which can be used to visualize similarity in a higher dimensional space using a 2D representation, to map these five-dimensional clusters in the feature space onto two-dimensional projections. To further analyze these clusters, using the SOM information, we agree with previous results that the sub-populations found in the feature space can be reasonably mapped onto three or five clusters. We explore where the "green valley" galaxies are mapped onto the SOM, indicating multiple interstitial populations within the green valley population. Finally, we use the projection of the SOM to verify whether morphological information provided by GalaxyZoo users, for example, if features are visible, can be mapped onto the SOM-generated map. Voting on whether galaxies are smooth, likely ellipticals, or "featured" can reasonably be separated but smaller morphological features (bar, spiral arms) can not. SOMs promise to be a useful tool to map and identify instructive sub-populations in multidimensional galaxy survey feature space, provided they are large enough. △ Less

Submitted 29 March, 2022; originally announced March 2022.

Comments: 14 pages, 14 figures, accepted by MNRAS

Deep ASKAP EMU Survey of the GAMA23 field: Properties of radio sources

Authors: Gülay Gürkan, I. Prandoni, A. O'Brien, W. Raja, L. Marchetti, M. Vaccari, S. Driver, E. Taylor, T. Franzen, M. J. I. Brown, S. Shabala, H. Andernach, A. M. Hopkins, R. P. Norris, D. Leahy, M. Bilicki, H. Farajollahi, T. Galvin, G. Heald, B. S. Koribalski, T. An, K. Warhurst

Abstract: We present the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the Galaxy and Mass Assembly (GAMA)-23h field. The survey was carried out at 887.5 MHz and covers a 83 square degree field. We imaged the calibrated visibility data, taken as part of the Evolutionary Mapping of Universe (EMU) Early Science Programme, using the latest version of the ASKAPSoft pipeline. The final mos… ▽ More We present the Australian Square Kilometre Array Pathfinder (ASKAP) observations of the Galaxy and Mass Assembly (GAMA)-23h field. The survey was carried out at 887.5 MHz and covers a 83 square degree field. We imaged the calibrated visibility data, taken as part of the Evolutionary Mapping of Universe (EMU) Early Science Programme, using the latest version of the ASKAPSoft pipeline. The final mosaic has an angular resolution of 10 arcsec and a central rms noise of around 38 $μ$Jy beam$^{-1}$. The derived radio source catalogue has 39812 entries above a peak flux density threshold of 5$σ$. We searched for the radio source host galaxy counterparts using the GAMA spectroscopic (with an i-band magnitude limit of 19.2 mag) and multi-wavelength catalogues that are available as part of the collaboration. We identified hosts with GAMA spectroscopic redshifts for 5934 radio sources. We describe the data reduction, imaging, and source identification process, and present the source counts. Thanks to the wide area covered by our survey, we obtain very robust counts down to 0.2 mJy. ASKAP's exceptional survey speed, providing efficient, sensitive and high resolution mapping of large regions of the sky in conjunction with the multi-wavelength data available for the GAMA23 field, allowed us to discover 63 giant radio galaxies. The data presented here demonstrate the excellent capabilities of ASKAP in the pre-SKA era. △ Less

Submitted 28 March, 2022; originally announced March 2022.

Comments: 18 pages, 15 figures, accepted by MNRAS

Galaxy And Mass Assembly (GAMA): Data Release 4 and the z < 0.1 total and z < 0.08 morphological galaxy stellar mass functions

Authors: Simon P. Driver, Sabine Bellstedt, Aaron S. G. Robotham, Ivan K. Baldry, Luke J. Davies, Jochen Liske, Danail Obreschkow, Edward N. Taylor, Angus H. Wright, Mehmet Alpaslan, Steven P. Bamford, Amanda E. Bauer, Joss Bland-Hawthorn, Maciej Bilicki, Matias Bravo, Sarah Brough, Sarah Casura, Michelle E. Cluver, Matthew Colless, Christopher J. Conselice, Scott M. Croom, Jelte de Jong, Franceso D'Eugenio, Roberto De Propris, Burak Dogruel , et al. (45 additional authors not shown)

Abstract: In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19… ▽ More In Galaxy And Mass Assembly Data Release 4 (GAMA DR4), we make available our full spectroscopic redshift sample. This includes 248682 galaxy spectra, and, in combination with earlier surveys, results in 330542 redshifts across five sky regions covering ~250deg^2. The redshift density, is the highest available over such a sustained area, has exceptionally high completeness (95 per cent to r_KIDS=19.65mag), and is well suited for the study of galaxy mergers, galaxy groups, and the low redshift (z<0.25) galaxy population. DR4 includes 32 value-added tables or Data Management Units (DMUs) that provide a number of measured and derived data products including GALEX, ESO KiDS, ESO VIKING, WISE and Herschel Space Observatory imaging. Within this release, we provide visual morphologies for 15330 galaxies to z<0.08, photometric redshift estimates for all 18million objects to r_KIDS~25mag, and stellar velocity dispersions for 111830 galaxies. We conclude by deriving the total galaxy stellar mass function (GSMF) and its sub-division by morphological class (elliptical, compact-bulge and disc, diffuse-bulge and disc, and disc only). This extends our previous measurement of the total GSMF down to 10^6.75 M_sol h^-2_70 and we find a total stellar mass density of rho_*=(2.97+/-0.04)x10^8 M_sol h_70 Mpc^-3 or Omega_*=(2.17+/-0.03)x10^-3 h^-1_70. We conclude that at z<0.1, the Universe has converted 4.9+/-0.1 per cent of the baryonic mass implied by Big Bang Nucleosynthesis into stars that are gravitationally bound within the galaxy population. △ Less

Submitted 16 March, 2022; originally announced March 2022.

Comments: Accepted for publication in MNRAS. GAMA Data Release 4 is available at: http://www.gama-survey.org/dr4/

The LOFAR Two-metre Sky Survey -- V. Second data release

Authors: T. W. Shimwell, M. J. Hardcastle, C. Tasse, P. N. Best, H. J. A. Röttgering, W. L. Williams, A. Botteon, A. Drabent, A. Mechev, A. Shulevski, R. J. van Weeren, L. Bester, M. Brüggen, G. Brunetti, J. R. Callingham, K. T. Chyży, J. E. Conway, T. J. Dijkema, K. Duncan, F. de Gasperin, C. L. Hale, M. Haverkorn, B. Hugo, N. Jackson, M. Mevius , et al. (81 additional authors not shown)

Abstract: In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for th… ▽ More In this data release from the LOFAR Two-metre Sky Survey (LoTSS) we present 120-168MHz images covering 27% of the northern sky. Our coverage is split into two regions centred at approximately 12h45m +44$^\circ$30' and 1h00m +28$^\circ$00' and spanning 4178 and 1457 square degrees respectively. The images were derived from 3,451hrs (7.6PB) of LOFAR High Band Antenna data which were corrected for the direction-independent instrumental properties as well as direction-dependent ionospheric distortions during extensive, but fully automated, data processing. A catalogue of 4,396,228 radio sources is derived from our total intensity (Stokes I) maps, where the majority of these have never been detected at radio wavelengths before. At 6" resolution, our full bandwidth Stokes I continuum maps with a central frequency of 144MHz have: a median rms sensitivity of 83$μ$Jy/beam; a flux density scale accuracy of approximately 10%; an astrometric accuracy of 0.2"; and we estimate the point-source completeness to be 90% at a peak brightness of 0.8mJy/beam. By creating three 16MHz bandwidth images across the band we are able to measure the in-band spectral index of many sources, albeit with an error on the derived spectral index of +/-0.2 which is a consequence of our flux-density scale accuracy and small fractional bandwidth. Our circular polarisation (Stokes V) 20" resolution 120-168MHz continuum images have a median rms sensitivity of 95$μ$Jy/beam, and we estimate a Stokes I to Stokes V leakage of 0.056%. Our linear polarisation (Stokes Q and Stokes U) image cubes consist of 480 x 97.6 kHz wide planes and have a median rms sensitivity per plane of 10.8mJy/beam at 4' and 2.2mJy/beam at 20"; we estimate the Stokes I to Stokes Q/U leakage to be approximately 0.2%. Here we characterise and publicly release our Stokes I, Q, U and V images in addition to the calibrated uv-data. △ Less

Submitted 23 February, 2022; originally announced February 2022.

Comments: 23 figures, 1 table and 29 pages. The catalogues, images and uv-data associated with this data release are publicly available via https://lofar-surveys.org/

Consistent lensing and clustering in a low-$S_8$ Universe with BOSS, DES Year 3, HSC Year 1 and KiDS-1000

Authors: A. Amon, N. C. Robertson, H. Miyatake, C. Heymans, M. White, J. DeRose, S. Yuan, R. H. Wechsler, T. N. Varga, S. Bocquet, A. Dvornik, S. More, A. J. Ross, H. Hoekstra, A. Alarcon, M. Asgari, J. Blazek, A. Campos, R. Chen, A. Choi, M. Crocce, H. T. Diehl, C. Doux, K. Eckert, J. Elvin-Poole , et al. (83 additional authors not shown)

Abstract: We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data a… ▽ More We evaluate the consistency between lensing and clustering probes of large-scale structure based on measurements of projected galaxy clustering from BOSS combined with overlapping galaxy-galaxy lensing from three surveys: DES Y3, HSC Y1, and KiDS-1000. An intra-lensing-survey study finds good agreement between these lensing data. We model the observations using the Dark Emulator and fit the data at two fixed cosmologies: Planck, with $S_8=0.83$, and a Lensing cosmology with $S_8=0.76$. For a joint analysis limited to scales with $R>5.25h^{-1}$Mpc, we find that both cosmologies provide an acceptable fit to the data. Full utilisation of the small-scale clustering and lensing measurements is hindered by uncertainty in the impact of baryon feedback and assembly bias, which we account for with a reasoned theoretical error budget. We incorporate a systematic scaling parameter for each redshift bin, $A$, that decouples the lensing and clustering to capture any inconsistency. When a wide range of scales ($0.15<R<60h^{-1}$Mpc) are incorporated, we find different results for the consistency of clustering and lensing between the two cosmologies. Limiting the analysis to the bins for which the impact of the selection of the lens sample is expected to be minimal, for the low-$S_8$ Lensing cosmology, the measurements are consistent with $A$=1; $A=0.91\pm0.04$ using DES+KiDS and $A=0.97\pm0.06$ using HSC. For the Planck cosmology case, we find a discrepancy: $A=0.79\pm0.03$ using DES+KiDS and $A=0.84\pm0.05$ using HSC. We demonstrate that a kSZ-based estimate for baryonic effects alleviates some of the discrepancy in the Planck cosmology. This analysis demonstrates the statistical power of these small-scale measurements, but also indicates that caution is still warranted given current uncertainties in modelling baryonic effects, assembly bias, and selection effects in the foreground sample. △ Less

Submitted 13 October, 2022; v1 submitted 15 February, 2022; originally announced February 2022.

Comments: 28 pages, 11 figures

Wavelength-resolved Reverberation Mapping of quasar CTSC30.10: Dissecting MgII and FeII emission regions

Authors: Raj Prince, Michal Zajaček, Bożena Czerny, P. Trzcionkowski, M. Bronikowski, C. S. Figaredo, S. Panda, M. L. Martinez-Aldama, K. Hryniewicz, V. K. Jaiswal, M. Śniegowska, M. Naddaf, M. Bilicki, M. Haas, M. J. Sarna, V. Karas, A. Olejak, R. Przyłuski, M. Rałowski, A. Udalski, R. R. Sefako, A. Genade, H. L. Worters

Abstract: We present the results of the reverberation monitoring aimed at MgII broad line and FeII pseudocontinuum for the luminous quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope covering the years 2012-2021. We aimed at disentangling the MgII and UV FeII variability and the first measurement of UV FeII time delay for a distant quasar. We used several methods for time-delay measur… ▽ More We present the results of the reverberation monitoring aimed at MgII broad line and FeII pseudocontinuum for the luminous quasar CTS C30.10 (z = 0.90052) with the Southern African Large Telescope covering the years 2012-2021. We aimed at disentangling the MgII and UV FeII variability and the first measurement of UV FeII time delay for a distant quasar. We used several methods for time-delay measurements and determined both FeII and MgII time delays as well as performed a wavelength-resolved time delay study for a combination of MgII and FeII in the 2700 - 2900 Årestframe wavelength range. We obtain the time delay for MgII of $275.5^{+12.4}_{-19.5}$ days in the rest frame, while for FeII we have two possible solutions of $270.0^{+13.8}_{-25.3}$ days and $180.3^{+26.6}_{-30.0}$ in the rest frame. Combining this result with the old measurement of FeII UV time delay for NGC 5548 we discuss for the first time the radius-luminosity relation for UV FeII with the slope consistent with $0.5$ within uncertainties. Since FeII time delay has a shorter time-delay component but lines are narrower than MgII, we propose that the line delay measurement is biased towards the BLR part facing the observer, with the bulk of the Fe II emission may arise from the more distant BLR region, one that is shielded from the observer. △ Less

Submitted 26 January, 2022; originally announced January 2022.

Comments: 22 pages, 19 Figures, 6 Tables, Submitted to Astronomy and Astrophysics, Comments are welcome

Journal ref: A&A 667, A42 (2022)

Galaxy and Mass Assembly (GAMA): Mid-infrared properties as tracers of galaxy environment

Authors: U. Sureshkumar, A. Durkalec, A. Pollo, M. Bilicki, M. E. Cluver, S. Bellstedt, D. J. Farrow, J. Loveday, E. N. Taylor, J. Bland-Hawthorn

Abstract: We investigate how different mid-infrared (mid-IR) properties of galaxies trace the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 $μ$m and redshift. Then, we look into the environmental dependence of mid-IR luminosities and the galaxy properties derived from these luminosities. We also explore ho… ▽ More We investigate how different mid-infrared (mid-IR) properties of galaxies trace the environment in which the galaxies are located. For this purpose, we first study the dependence of galaxy clustering on the absolute magnitude at 3.4 $μ$m and redshift. Then, we look into the environmental dependence of mid-IR luminosities and the galaxy properties derived from these luminosities. We also explore how various infrared galaxy luminosity selections influence the galaxy clustering measurements. We use a set of W1 (3.4 $μ$m) absolute magnitude ($M_\text{W1}$) selected samples from the Galaxy and Mass Assembly (GAMA) survey matched with mid-IR properties from the Wide-field Infrared Survey Explorer (WISE) in the redshift range $0.07 \leq z < 0.43$. We compute the galaxy two-point correlation function (2pCF) and compare the clustering lengths between subsamples binned in $M_\text{W1}$ and in redshift. We also measure the marked correlation functions (MCFs) using the luminosities in the WISE W1 to W4 (3.4 to 22 $μ$m) bands as marks. Additionally, we compare the measurements of MCFs with different estimates of stellar mass and star formation rate used as marks. Finally, we check how different selections applied to the sample affect the clustering measurements. We show strong clustering dependence on the W1 absolute magnitude: galaxies brighter in the W1 band are more strongly clustered than their fainter counterparts. We also observe a lack of significant redshift dependence of clustering in the redshift range $0.07 \leq z < 0.43$. We show that although W1 and W2 bands are direct indicators of stellar mass, a galaxy sample selected based on W1 or W2 bands does not perfectly show the clustering behaviour of a stellar mass selected sample. Similar is the case with the proxy relation between W3 and W4 bands and star formation rate. △ Less

Submitted 4 January, 2023; v1 submitted 25 January, 2022; originally announced January 2022.

Comments: 18 pages, 12 figures, matches the version published in A&A

Journal ref: A&A 669, A27 (2023)

COWS: A filament finder for Hessian cosmic web identifiers

Authors: Simon Pfeifer, Noam I. Libeskind, Yehuda Hoffman, Wojciech A. Hellwing, Maciej Bilicki, Krishna Naidoo

Abstract: The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells… ▽ More The large scale galaxy and matter distribution is often described by means of the cosmic web made up of voids, sheets, filaments and knots. Many different recipes exist for identifying this cosmic web. Here we focus on a sub-class of cosmic web identifiers, based on the analysis of the Hessian matrix, and proposed a method, called COsmic Web Skeleton (COWS), of separating a set of filaments cells into an ensemble of individual discreet filaments. Specifically, a thinning algorithm is applied to velocity shear tensor based cosmic web (V-web) to identify the spine of the filaments. This results in a set of filaments with well defined end-point and length. It is confirmed that these sit at local density ridges and align with the appropriate direction defined by the underlying velocity field. The radial density profile of these curved cylindrical filaments, as well as the distribution of their lengths is also examined. The robustness of all results is checked against resolution and the V-web threshold. The code for the COWS method has been made publicly available. △ Less

Submitted 12 January, 2022; originally announced January 2022.

Comments: 10 pages, 10 figures, submitted to MNRAS

Universality of the halo mass function in modified gravity cosmologies

Authors: Suhani Gupta, Wojciech A. Hellwing, Maciej Bilicki, Jorge Enrique García-Farieta

Abstract: We study the halo mass function (HMF) in modified gravity (MG) models using a set of large $N$-body simulations -- the ELEPHANT suite. We consider two popular beyond-general relativity scenarios: the Hu-Sawicki chameleon $f(R)$ model and the normal branch of the Dvali-Gabadadze-Porrati braneworld (nDGP). We show that in MG, analytic formulation based on the Press-Schechter framework offers a gross… ▽ More We study the halo mass function (HMF) in modified gravity (MG) models using a set of large $N$-body simulations -- the ELEPHANT suite. We consider two popular beyond-general relativity scenarios: the Hu-Sawicki chameleon $f(R)$ model and the normal branch of the Dvali-Gabadadze-Porrati braneworld (nDGP). We show that in MG, analytic formulation based on the Press-Schechter framework offers a grossly inaccurate description of the HMF. We find, however, that once the HMF is expressed in terms of the dimensionless multiplicity function, it approximately assumes a redshift-independent universal character for all the models. Exploiting this property, we propose universal fits for the MG HMF in terms of their fractional departures from the $\mathrm{Λ\text{CDM}}$ case. We find two enclosed formulas, one for $f(R)$ and another for nDGP, that provide a reliable description of the HMF over the mass range covered by the simulations. These are accurate to a few percent with respect to the $N$-body data. We test the extrapolation potential of our fits against separate simulations with a different cosmological background and mass resolution, and find very good accuracy, within $\sim 10\%$. A particularly interesting finding from our analysis is a Gaussian-like shape of the HMF deviation that seems to appear universally across the whole $f(R)$ family, peaking at a mass variance scale characteristic for each $f(R)$ variant. We attribute this behavior to the specific physics of the environmentally-dependent chameleon screening models. △ Less

Submitted 18 March, 2022; v1 submitted 7 December, 2021; originally announced December 2021.

Comments: 16 pages and 8 figures

Journal ref: Physical Review D, Volume 105, Issue 4, article id.043538 (2022)

The dark matter halo masses of elliptical galaxies as a function of observationally robust quantities

Authors: Alessandro Sonnenfeld, Crescenzo Tortora, Henk Hoekstra, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Koen Kuijken, Nicola R. Napolitano, Nivya Roy, Edwin Valentijn, Angus H. Wright

Abstract: Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing ob… ▽ More Context. The assembly history of the stellar component of a massive elliptical galaxy is closely related to that of its dark matter halo. Measuring how the properties of galaxies correlate with their halo mass can help understand their evolution. Aims. We investigate how the dark matter halo mass of elliptical galaxies varies as a function of their properties, using weak gravitational lensing observations. To minimise the chances of biases, we focus on galaxy properties that can be determined robustly: the surface brightness profile and the colour. Methods. We selected 2409 central massive elliptical galaxies from the SDSS spectroscopic sample. We first measured their surface brightness profile and colours by fitting Sersic models to photometric data from the Kilo-Degree Survey (KiDS). We fitted their halo mass distribution as a function of redshift, rest-frame $r-$band luminosity, half-light radius and rest-frame $u-g$ colour, using KiDS weak lensing data and a Bayesian hierarchical approach. For the sake of robustness to assumptions on the large-radii behaviour of the surface brightness, we repeated the analysis replacing total luminosity and half-light radius with the luminosity within a 10~kpc aperture, $L_{r,10}$, and the light-weighted surface brightness slope, $Γ_{10}$. Results. We did not detect any correlation between halo mass and either half-light radius or colour, at fixed redshift and luminosity. Conclusions. Our results indicate that the average star formation efficiency of massive elliptical galaxies has little dependence on their final size or colour. This suggests that the origin of the diversity in the size and colour distribution of these objects lies with properties other than halo mass. △ Less

Submitted 15 April, 2022; v1 submitted 22 October, 2021; originally announced October 2021.

Comments: In press on Astronomy & Astrophysics. A short summary video is available at https://youtu.be/LlvwQNjLT7k

Journal ref: A&A 662, A55 (2022)

KiDS & Euclid: Cosmological implications of a pseudo angular power spectrum analysis of KiDS-1000 cosmic shear tomography

Authors: A. Loureiro, L. Whittaker, A. Spurio Mancini, B. Joachimi, A. Cuceu, M. Asgari, B. Stölzner, T. Tröster, A. H. Wright, M. Bilicki, A. Dvornik, B. Giblin, C. Heymans, H. Hildebrandt, H. Shan, A. Amara, N. Auricchio, C. Bodendorf, D. Bonino, E. Branchini, M. Brescia, V. Capobianco, C. Carbone, J. Carretero, M. Castellano , et al. (89 additional authors not shown)

Abstract: We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using… ▽ More We present a tomographic weak lensing analysis of the Kilo Degree Survey Data Release 4 (KiDS-1000), using a new pseudo angular power spectrum estimator (pseudo-$C_{\ell}$) under development for the ESA Euclid mission. Over 21 million galaxies with shape information are divided into five tomographic redshift bins, ranging from 0.1 to 1.2 in photometric redshift. We measured pseudo-$C_{\ell}$ using eight bands in the multipole range $76<\ell<1500$ for auto- and cross-power spectra between the tomographic bins. A series of tests were carried out to check for systematic contamination from a variety of observational sources including stellar number density, variations in survey depth, and point spread function properties. While some marginal correlations with these systematic tracers were observed, there is no evidence of bias in the cosmological inference. B-mode power spectra are consistent with zero signal, with no significant residual contamination from E/B-mode leakage. We performed a Bayesian analysis of the pseudo-$C_{\ell}$ estimates by forward modelling the effects of the mask. Assuming a spatially flat $Λ$CDM cosmology, we constrained the structure growth parameter $S_8 = σ_8(Ω_{\rm m}/0.3)^{1/2} = 0.754_{-0.029}^{+0.027}$. When combining cosmic shear from KiDS-1000 with baryon acoustic oscillation and redshift space distortion data from recent Sloan Digital Sky Survey (SDSS) measurements of luminous red galaxies, as well as the Lyman-$α$ forest and its cross-correlation with quasars, we tightened these constraints to $S_8 = 0.771^{+0.006}_{-0.032}$. These results are in very good agreement with previous KiDS-1000 and SDSS analyses and confirm a $\sim 3σ$ tension with early-Universe constraints from cosmic microwave background experiments. △ Less

Submitted 4 July, 2022; v1 submitted 13 October, 2021; originally announced October 2021.

Comments: 24 pages, 17 figures, Accepted by A&A

Journal ref: A&A 665, A56 (2022)

GLADE+: An Extended Galaxy Catalogue for Multimessenger Searches with Advanced Gravitational-wave Detectors

Authors: G. Dálya, R. Díaz, F. R. Bouchet, Z. Frei, J. Jasche, G. Lavaux, R. Macas, S. Mukherjee, M. Pálfi, R. S. de Souza, B. D. Wandelt, M. Bilicki, P. Raffai

Abstract: We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-fram… ▽ More We present GLADE+, an extended version of the GLADE galaxy catalogue introduced in our previous paper for multimessenger searches with advanced gravitational-wave detectors. GLADE+ combines data from six separate but not independent astronomical catalogues: the GWGC, 2MPZ, 2MASS XSC, HyperLEDA, and WISExSCOSPZ galaxy catalogues, and the SDSS-DR16Q quasar catalogue. To allow corrections of CMB-frame redshifts for peculiar motions, we calculated peculiar velocities along with their standard deviations of all galaxies having $B$-band magnitude data within redshift $z=0.05$ using the "Bayesian Origin Reconstruction from Galaxies" formalism. GLADE+ is complete up to luminosity distance $d_L=47^{+4}_{-2}$ Mpc in terms of the total expected $B$-band luminosity of galaxies, and contains all of the brightest galaxies giving 90\% of the total $B$-band and $K$-band luminosity up to $d_L\simeq 130$ Mpc. We include estimations of stellar masses and individual binary neutron star merger rates for galaxies with $W1$ magnitudes. These parameters can help in ranking galaxies in a given gravitational wave localization volume in terms of their likelihood of being hosts, thereby possibly reducing the number of pointings and total integration time needed to find the electromagnetic counterpart. △ Less

Submitted 2 June, 2022; v1 submitted 12 October, 2021; originally announced October 2021.

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

Galaxy And Mass Assembly (GAMA): $\mathbf{z \sim 0}$ Galaxy Luminosity Function down to $\mathbf{L \sim 10^{6}~L_\odot}$ via Clustering Based Redshift Inference

Authors: Geray S. Karademir, Edward N. Taylor, Chris Blake, Ivan K. Baldry, Sabine Bellstedt, Maciej Bilicki, Michael J. I. Brown, Michelle E. Cluver, Simon P. Driver, Hendrik Hildebrandt, Benne W. Holwerda, Andrew M. Hopkins, Jonathan Loveday, Steven Phillipps, Angus H. Wright

Abstract: In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of… ▽ More In this study we present a new experimental design using clustering-based redshift inference to measure the evolving galaxy luminosity function (GLF) spanning 5.5 decades from $L \sim 10^{11.5}$ to $ 10^6 ~ \mathrm{L}_\odot$. We use data from the Galaxy And Mass Assembly (GAMA) survey and the Kilo-Degree Survey (KiDS). We derive redshift distributions in bins of apparent magnitude to the limits of the GAMA-KiDS photometric catalogue: $m_r \lesssim 23$; more than a decade in luminosity beyond the limits of the GAMA spectroscopic redshift sample via clustering-based redshift inference. This technique uses spatial cross-correlation statistics for a reference set with known redshifts (in our case, the main GAMA sample) to derive the redshift distribution for the target ensemble. For the calibration of the redshift distribution we use a simple parametrisation with an adaptive normalisation factor over the interval $0.005 < z < 0.48$ to derive the clustering redshift results. We find that the GLF has a relatively constant power-law slope $α\approx -1.2$ for $-17 \lesssim M_r \lesssim -13$, and then appears to steepen sharply for $-13 \lesssim M_r \lesssim -10$. This upturn appears to be where Globular Clusters (GCs) take over to dominate the source counts as a function of luminosity. Thus we have mapped the GLF across the full range of the $z \sim 0$ field galaxy population from the most luminous galaxies down to the GC scale. △ Less

Submitted 3 November, 2021; v1 submitted 8 September, 2021; originally announced September 2021.

Comments: 21 pages, 13 figures, accepted by MNRAS

Joint constraints on cosmology and the impact of baryon feedback: combining KiDS-1000 lensing with the thermal Sunyaev-Zeldovich effect from Planck and ACT

Authors: Tilman Tröster, Alexander J. Mead, Catherine Heymans, Ziang Yan, David Alonso, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Hendrik Hildebrandt, Benjamin Joachimi, Arun Kannawadi, Konrad Kuijken, Peter Schneider, HuanYuan Shan, Ludovic van Waerbeke, Angus H. Wright

Abstract: We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive const… ▽ More We conduct a pseudo-$C_\ell$ analysis of the tomographic cross-correlation between 1000 deg$^2$ of weak lensing data from the Kilo-Degree Survey (KiDS-1000) and the thermal Sunyaev-Zeldovich (tSZ) effect measured by Planck and the Atacama Cosmology Telescope (ACT). Using HMx, a halo-model-based approach that consistently models the gas, star, and dark matter components, we are able to derive constraints on both cosmology and baryon feedback for the first time from these data, marginalising over redshift uncertainties, intrinsic alignment of galaxies, and contamination by the cosmic infrared background (CIB). We find our results to be insensitive to the CIB, while intrinsic alignment provides a small but significant contribution to the lensing--tSZ cross-correlation. The cosmological constraints are consistent with those of other low-redshift probes and prefer strong baryon feedback. The inferred amplitude of the lensing--tSZ cross-correlation signal, which scales as $σ_8(Ω_\mathrm{m}/0.3)^{0.2}$, is low by $\sim 2\,σ$ compared to the primary cosmic microwave background constraints by Planck. The lensing--tSZ measurements are then combined with pseudo-$C_\ell$ measurements of KiDS-1000 cosmic shear into a novel joint analysis, accounting for the full cross-covariance between the probes, providing tight cosmological constraints by breaking parameter degeneracies inherent to both probes. The joint analysis gives an improvement of 40% on the constraint of $S_8=σ_8\sqrt{Ω_\mathrm{m}/0.3}$ over cosmic shear alone, while providing constraints on baryon feedback consistent with hydrodynamical simulations, demonstrating the potential of such joint analyses with baryonic tracers such as the tSZ effect. We discuss remaining modelling challenges that need to be addressed if these baryonic probes are to be included in future precision-cosmology analyses. △ Less

Submitted 8 February, 2022; v1 submitted 9 September, 2021; originally announced September 2021.

Comments: 21 pages, 16 figures. Accepted in A&A. Comments welcome. Code, data vectors, and covariances are available at https://github.com/tilmantroester/KiDS-1000xtSZ

Journal ref: A&A 660, A27 (2022)

KiDS-1000: Constraints on the intrinsic alignment of luminous red galaxies

Authors: Maria Cristina Fortuna, Henk Hoekstra, Harry Johnston, Mohammadjavad Vakili, Arun Kannawadi, Christos Georgiou, Benjamin Joachimi, Angus H. Wright, Marika Asgari, Maciej Bilicki, Catherine Heymans, Hendrik Hildebrandt, Konrad Kuijken, Maximilian Von Wietersheim-Kramsta

Abstract: We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of… ▽ More We constrain the luminosity and redshift dependence of the intrinsic alignment (IA) of a nearly volume-limited sample of luminous red galaxies selected from the fourth public data release of the Kilo-Degree Survey (KiDS-1000). To measure the shapes of the galaxies, we used two complementary algorithms, finding consistent IA measurements for the overlapping galaxy sample. The global significance of IA detection across our two independent luminous red galaxy samples, with our favoured method of shape estimation, is $\sim10.7σ$. We find no significant dependence with redshift of the IA signal in the range $0.2<z<0.8$, nor a dependence with luminosity below $L_r\lesssim 2.9 \times 10^{10} h^{-2} L_{r,\odot}$. Above this luminosity, however, we find that the IA signal increases as a power law, although our results are also compatible with linear growth within the current uncertainties. This behaviour motivates the use of a broken power law model when accounting for the luminosity dependence of IA contamination in cosmic shear studies. △ Less

Submitted 6 September, 2021; originally announced September 2021.

Comments: 21 pages, 14 figures. Accepted for publication in A&A

The Weak Lensing Radial Acceleration Relation: Constraining Modified Gravity and Cold Dark Matter theories with KiDS-1000

Authors: Margot M. Brouwer, Kyle A. Oman, Edwin A. Valentijn, Maciej Bilicki, Catherine Heymans, Henk Hoekstra, Nicola R. Napolitano, Nivya Roy, Crescenzo Tortora, Angus H. Wright, Marika Asgari, Jan Luca van den Busch, Andrej Dvornik, Thomas Erben, Benjamin Giblin, Alister W. Graham, Hendrik Hildebrandt, Andrew M. Hopkins, Arun Kannawadi, Konrad Kuijken, Jochen Liske, HuanYuan Shan, Tilman Tröster, Erik Verlinde, Manus Visser

Abstract: We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-accelera… ▽ More We present measurements of the radial gravitational acceleration around isolated galaxies, comparing the expected gravitational acceleration given the baryonic matter with the observed gravitational acceleration, using weak lensing measurements from the fourth data release of the Kilo-Degree Survey. These measurements extend the radial acceleration relation (RAR) by 2 decades into the low-acceleration regime beyond the outskirts of the observable galaxy. We compare our RAR measurements to the predictions of two modified gravity (MG) theories: MOND and Verlinde's emergent gravity. We find that the measured RAR agrees well with the MG predictions. In addition, we find a difference of at least $6σ$ between the RARs of early- and late-type galaxies (split by Sérsic index and $u-r$ colour) with the same stellar mass. Current MG theories involve a gravity modification that is independent of other galaxy properties, which would be unable to explain this behaviour. The difference might be explained if only the early-type galaxies have significant ($M_{gas} \approx M_*$) circumgalactic gaseous haloes. The observed behaviour is also expected in $Λ$CDM models where the galaxy-to-halo mass relation depends on the galaxy formation history. We find that MICE, a $Λ$CDM simulation with hybrid halo occupation distribution modelling and abundance matching, reproduces the observed RAR but significantly differs from BAHAMAS, a hydrodynamical cosmological galaxy formation simulation. Our results are sensitive to the amount of circumgalactic gas; current observational constraints indicate that the resulting corrections are likely moderate. Measurements of the lensing RAR with future cosmological surveys will be able to further distinguish between MG and $Λ$CDM models if systematic uncertainties in the baryonic mass distribution around galaxies are reduced. △ Less

Submitted 22 June, 2021; originally announced June 2021.

Comments: 31 pages, 15 figures, published in A&A

Journal ref: Astronomy & Astrophysics Volume 650 (2021) ArticleID A113

Deep Extragalactic VIsible Legacy Survey (DEVILS): Consistent multi-wavelength photometry for the DEVILS regions (COSMOS, XMMLSS & ECDFS)

Authors: L. J. M. Davies, J. E. Thorne, A. S. G. Robotham, S. Bellstedt, S. P. Driver, N. J. Adams, M. Bilicki, R. A. A. Bowler, M. Bravo, L. Cortese, C. Foster, M. W. Grootes, B. Häußler, A. Hashemizadeh, B. W. Holwerda, P. Hurley, M. J. Jarvis, C. Lidman, N. Maddox, M. Meyer, M. Paolillo, S. Phillipps, M. Radovich, M. Siudek, M. Vaccari , et al. (1 additional authors not shown)

Abstract: The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy… ▽ More The Deep Extragalactic VIsible Legacy Survey (DEVILS) is an ongoing high-completeness, deep spectroscopic survey of $\sim$60,000 galaxies to Y$<$21.2 mag, over $\sim$6 deg2 in three well-studied deep extragalactic fields: D10 (COSMOS), D02 (XMM-LSS) and D03 (ECDFS). Numerous DEVILS projects all require consistent, uniformly-derived and state-of-the-art photometric data with which to measure galaxy properties. Existing photometric catalogues in these regions either use varied photometric measurement techniques for different facilities/wavelengths leading to inconsistencies, older imaging data and/or rely on source detection and photometry techniques with known problems. Here we use the ProFound image analysis package and state-of-the-art imaging datasets (including Subaru-HSC, VST-VOICE, VISTA-VIDEO and UltraVISTA-DR4) to derive matched-source photometry in 22 bands from the FUV to 500μm. This photometry is found to be consistent, or better, in colour-analysis to previous approaches using fixed-size apertures (which are specifically tuned to derive colours), but produces superior total source photometry, essential for the derivation of stellar masses, star-formation rates, star-formation histories, etc. Our photometric catalogue is described in detail and, after internal DEVILS team projects, will be publicly released for use by the broader scientific community. △ Less

Submitted 11 June, 2021; originally announced June 2021.

Comments: 33 pages, 22 figures, Accepted to MNRAS

Geometry vs growth: Internal consistency of the flat ΛCDM model with KiDS-1000

Authors: J. Ruiz-Zapatero, Benjamin Stölzner, Benjamin Joachimi, Marika Asgari, Maciej Bilicki, Andrej Dvornik, Benjamin Giblin, Catherine Heymans, Hendrik Hildebrandt, Arun Kannawadi, Konrad Kuijken, TilmanTröster, Jan Luca van den Busch, Angus H. Wright

Abstract: We carry out a multi-probe self-consistency test of the flat $Λ$CDM model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of Lambda Cold… ▽ More We carry out a multi-probe self-consistency test of the flat $Λ$CDM model with the aim of exploring potential causes of the reported tensions between high- and low-redshift cosmological observations. We divide the model into two theory regimes determined by the smooth background (geometry) and the evolution of matter density fluctuations (growth), each governed by an independent set of Lambda Cold Dark Matter ($Λ$CDM) cosmological parameters. This extended model is constrained by a combination of weak gravitational lensing measurements from the Kilo-Degree Survey, galaxy clustering signatures extracted from Sloan Digital Sky Survey campaigns and the Six-Degree Field Galaxy Survey, and the angular baryon acoustic scale and the primordial scalar fluctuation power spectrum measured in $\textit{Planck}$ cosmic microwave background (CMB) data. We find strong consistency between the geometry and growth parameters, and with the posterior of standard $Λ$CDM analysis. Tension in the amplitude of matter density fluctuations as measured by the parameter $S_8$ persists at around 3$σ$, with a $1.5\,\%$ constraint of $S_8 = 0.776_{-0.008}^{+0.016}$ for the combined probes. We also observe a less significant preference (at least $2σ$) for higher values of the Hubble constant, $H_0 = 70.5^{+0.7}_{-1.5}\,{\rm km\, s^{-1} Mpc^{-1}}$, as well as for lower values of the total matter density parameter $Ω_{\rm{m}} = 0.289^{+0.007}_{-0.005}$ compared to the full $\textit{Planck}$ analysis. Including the subset of the CMB information in the probe combination enhances these differences rather than alleviate them, which we link to the discrepancy between low and high multipoles in $\textit{Planck}$ data. △ Less

Submitted 8 September, 2021; v1 submitted 20 May, 2021; originally announced May 2021.

Comments: 25 pages, 14 figures, updated after being accepted for publication to Astronomy and Astrophysics

Journal ref: A&A 655, A11 (2021)