Search

Context. Weak gravitational lensing offers a powerful method to investigate the projected matter density distribution within galaxy clusters, granting crucial insights into the broader landscape of dark matter on cluster scales. Aims. In this study, we make use of the large photometric galaxy cluster data set derived from the publicly available Third Data Release of the Kilo-Degree Survey, along with the associated shear signal. Our primary objective is to model the peculiar sharp transition in the cluster profile slope, that is what is commonly referred to as the splashback radius. The data set under scrutiny includes 6962 galaxy clusters, selected by AMICO (an optimised detection algorithm of galaxy clusters) on the KiDS-DR3 data, in the redshift range of 0.1 ≤ z ≤ 0.6, all observed at a signal-to-noise ratio greater than 3.5. Methods. Employing a comprehensive Bayesian analysis, we model the stacked excess surface mass density distribution of the clusters. We adopt a model from recent results on numerical simulations that capture the dynamics of both orbiting and infalling materials, separated by the region where the density profile slope undergoes a pronounced deepening. Results. We find that the adopted profile successfully characterizes the cluster masses, consistent with previous works, and models the deepening of the slope of the density profiles measured with weak-lensing data up to the outskirts. Moreover, we measure the splashback radius of galaxy clusters and show that its value is close to the radius within which the enclosed overdensity is 200 times the mean matter density of the Universe, while theoretical models predict a larger value consistent with a low accretion rate. This points to a potential bias of optically selected clusters preferentially characterized by a high density at small scales compared to a pure mass-selected cluster sample. [ABSTRACT FROM AUTHOR]

Although it is well known that the bulk of dark matter (DM) has to be cold, the existence of an additional sub-dominant, hot species remains a valid possibility. In this paper we investigate the potential of the cosmic shear power spectrum to constrain such a mixed (hot plus cold) DM scenario with two additional free parameters, the hot-to-total DM fraction (fhdm) and the thermal mass of the hot component (mhdm). Running a Bayesian inference analysis for both the Kilo-Degree Survey cosmic shear data (KiDS-1000) as well as the cosmic microwave background (CMB) temperature and polarisation data from Planck, we derive new constraints for the mixed DM scenario. We find a 95% confidence limit of fhdm < 0.08 for a very hot species of mhdm ≤ 20 eV. This constraint is weakened to fhdm < 0.25 for mhdm ≤ 80 eV. Scenarios with masses above mhdm ∼ 200 eV remain unconstrained by the data. Next to providing limits, we investigate the potential of mixed DM to address the clustering (or S8) tension between lensing and the CMB. We find a reduction of the 2D (Ωm − S8) tension from 2.9σ to 1.6σ when going from a pure cold DM to a mixed DM scenario. When computing the 1D Gaussian tension on S8 the improvement is milder, from 2.4σ to 2.0σ. [ABSTRACT FROM AUTHOR]

Context. This Letter explores the potential role of primordial black holes (PBHs) to address cosmological tensions as the presence of more massive than expected galaxies at high redshifts, as indicated by recent James Webb Space Telescope observations. Aims. Motivated by inflation models that enhance the power at scales beyond the observable range that produce PBHs with Schechter-like mass functions, we aim to explain the excess of high redshift galaxies via a modification of the Λ cold dark matter power spectrum that consists in adding (i) a blue spectral index nb at kpiv = 10/Mpc and (ii) Poisson and isocurvature contributions from massive PBHs that only make up 0.5% of the dark matter. Methods. We simulated these models using the SWIFT code and find an increased abundance of high redshift galaxies in simulations that include PBHs. We compared these models to estimates from James Webb Space Telescope observations. Results. Unlike the Λ cold dark matter model, the inclusion of PBHs allowed us to reproduce the observations with reasonable values for the star formation efficiency. Furthermore, the power spectra we adopted potentially produce PBHs that can serve as seeds for supermassive black holes with masses 7.57 × 104 M⊙. [ABSTRACT FROM AUTHOR]

Aims. We derived galaxy colour selections from Euclid and ground-based photometry, aiming to accurately define background galaxy samples in cluster weak-lensing analyses. These selections have been implemented in the Euclid data analysis pipelines for galaxy clusters. Methods. Given any set of photometric bands, we developed a method for the calibration of optimal galaxy colour selections that maximises the selection completeness, given a threshold on purity. Such colour selections are expressed as a function of the lens redshift. Results. We calibrated galaxy selections using simulated ground-based griz and EuclidYEJEHE photometry. Both selections produce a purity higher than 97%. The griz selection completeness ranges from 30% to 84% in the lens redshift range zl ∈ [0.2, 0.8]. With the full grizYEJEHE selection, the completeness improves by up to 25 percentage points, and the zl range extends up to zl = 1.5. The calibrated colour selections are stable to changes in the sample limiting magnitudes and redshift, and the selection based on griz bands provides excellent results on real external datasets. Furthermore, the calibrated selections provide stable results using alternative photometric aperture definitions obtained from different ground-based telescopes. The griz selection is also purer at high redshift and more complete at low redshift compared to colour selections found in the literature. We find excellent agreement in terms of purity and completeness between the analysis of an independent, simulated Euclid galaxy catalogue and our calibration sample, except for galaxies at high redshifts, for which we obtain up to 50 percentage points higher completeness. The combination of colour and photo-z selections applied to simulated Euclid data yields up to 95% completeness, while the purity decreases down to 92% at high zl. We show that the calibrated colour selections provide robust results even when observations from a single band are missing from the ground-based data. Finally, we show that colour selections do not disrupt the shear calibration for stage III surveys. The first Euclid data releases will provide further insights into the impact of background selections on the shear calibration. [ABSTRACT FROM AUTHOR]

We present cosmological parameter constraints using maps from the last Planck data release (PR4). In particular, we detail an upgraded version of the cosmic microwave background likelihood, HiLLiPoP, that is based on angular power spectra and relies on a physical modeling of the foreground residuals in the spectral domain. This new version of the likelihood retains a larger sky fraction (up to 75%) and uses an extended multipole range. Using this likelihood, along with low-ℓ measurements from LoLLiPoP, we derived constraints on ΛCDM parameters that are in good agreement with previous Planck 2018 results, but with smaller uncertainties by 10% to 20%. We demonstrate that the foregrounds can be accurately described in the spectral domain, with a negligible impact on ΛCDM parameters. We also derived constraints on single-parameter extensions to ΛCDM, including AL, ΩK, Neff, and ∑mν. Noteworthy results from this updated analysis include a lensing amplitude value of AL = 1.039 ± 0.052, which is more closely aligned with theoretical expectations within the ΛCDM framework. Additionally, our curvature measurement, ΩK = −0.012 ± 0.010, is now fully consistent with a flat universe and our measurement of S8 is closer to the measurements derived from large-scale structure surveys (at the 1.5σ level). We also added constraints from PR4 lensing, making this combination the most tightly constrained data set currently available from Planck. Additionally, we explored the addition of baryon acoustic oscillation data, which tightens the limits on some particular extensions to the standard cosmology. [ABSTRACT FROM AUTHOR]

We investigate the main tensions within the current standard model of cosmology from the perspective of the main statistics of cosmic voids, using the final BOSS DR12 data set. For this purpose, we present the first estimate of the S8 ≡ σ8 Ωm/0.3 and H0 parameters obtained from void number counts and shape distortions. To analyze void counts we relied on an extension of the popular volume-conserving model for the void size function, tailored to the application on data, including geometric and dynamic distortions. We calibrated the two nuisance parameters of this model with the official BOSS Collaboration mock catalogs and propagated their uncertainty through the statistical analysis of the BOSS void number counts. The constraints from void shapes come from the study of the geometric distortions of the stacked void-galaxy cross-correlation function. In this work we focus our analysis on the Ωm − σ8 and Ωm − H0 parameter planes and derive the marginalized constraints S8 = 0.813−0.068+0.093 and H0 = 67.3−9.1+10.0 km s−1 Mpc−1, which are fully compatible with constraints from the literature. These results are expected to notably improve in precision when analyzed jointly with independent probes and will open a new viewing angle on the rising cosmological tensions in the near future. [ABSTRACT FROM AUTHOR]

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

We present ICAROGW 2.0, a pure python code developed to infer the astrophysical and cosmological population properties of noisy, heterogeneous, and incomplete observations. The code has mainly been developed for compact binary coalescence (CBC) population inference with gravitational wave (GW) observations. It contains several models for the masses, spins, and redshift of CBC distributions and it is able to infer population distributions, as well as the cosmological parameters and possible general relativity deviations at cosmological scales. Here, we present the theoretical and computational foundations of ICAROGW 2.0 and describe how the code can be employed for population and cosmological inference using (i) only GWs, (ii) GWs and galaxy surveys, and (iii) GWs with electromagnetic counterparts. We discuss the code performance on GPUs, finding a gain in computation time of about two orders of magnitude when more than 100 GW events are involved in the analysis. We have validated the code by re-analyzing GW population and cosmological studies, finding very good agreement with previous results. [ABSTRACT FROM AUTHOR]

Context. The Euclid mission of the European Space Agency will perform a survey of weak lensing cosmic shear and galaxy clustering in order to constrain cosmological models and fundamental physics. Aims. We expand and adjust the mock Euclid likelihoods of the MontePython software in order to match the exact recipes used in previous Euclid Fisher matrix forecasts for several probes: weak lensing cosmic shear, photometric galaxy clustering, the cross-correlation between the latter observables, and spectroscopic galaxy clustering. We also establish which precision settings are required when running the Einstein–Boltzmann solvers CLASS and CAMB in the context of Euclid. Methods. For the minimal cosmological model, extended to include dynamical dark energy, we perform Fisher matrix forecasts based directly on a numerical evaluation of second derivatives of the likelihood with respect to model parameters. We compare our results with those of previously validated Fisher codes using an independent method based on first derivatives of the Euclid observables. Results. We show that such MontePython forecasts agree very well with previous Fisher forecasts published by the Euclid Collab oration, and also, with new forecasts produced by the CosmicFish code, now interfaced directly with the two Einstein–Boltzmann solvers CAMB and CLASS. Moreover, to establish the validity of the Gaussian approximation, we show that the Fisher matrix marginal error contours coincide with the credible regions obtained when running Monte Carlo Markov chains with MontePython while using the exact same mock likelihoods. Conclusions. The new Euclid forecast pipelines presented here are ready for use with additional cosmological parameters, in order to explore extended cosmological models. [ABSTRACT FROM AUTHOR]

Recent detections of extended Lyman-α halos around Lyα emitters (LAEs) have been reported on a regular basis, but their origin is still under investigation. Simulation studies predict that the outer regions of the extended halos contain a major contribution from the Lyα emission of faint, individually undetected LAEs. To address this matter from an observational angle, we used halo occupation distribution (HOD) modeling to reproduce the clustering of a spectroscopic sample of 1265 LAEs at 3 < z < 5 from the MUSE-Wide survey. We integrated the Lyα luminosity function to estimate the background surface brightness due to discrete faint LAEs. We then extended the HOD statistics inwards towards small separations and computed the factor by which the measured Lyα surface brightness is enhanced by undetected close physical neighbors. We considered various clustering scenarios for the undetected sources and compared the corresponding radial profiles. This enhancement factor from LAE clustering depends strongly on the spectral bandwidth Δv over which the Lyα emission is integrated and this value can amount to ≈20 − 40 for small values of Δv (around 200 − 400 km s−1) as achieved by recent studies utilizing integral-field spectrographic data. The resulting inferred Lyα surface brightness of faint LAEs ranges between (0.4 − 2)×1020 erg s−1 cm−2 arcsec−2, with a very slow radial decline outwards. Our results suggest that the outer regions of observed Lyα halos (R ≳ 50 pkpc) could indeed contain a strong component from external (but physically associated) LAEs, and may even be dominated by them. It is only for a relatively shallow faint-end slope of the Lyα luminosity function that this contribution from clustered LAEs would be rendered insignificant. We also confirm that the observed emission from the inner regions (R ≤ 20 − 30 pkpc) is too bright to be substantially affected by clustering. We compare our findings with predicted profiles from simulations and find good overall agreement. We outline possible future measurements to further constrain the impact of discrete undetected LAEs on observed extended Lyα halos. [ABSTRACT FROM AUTHOR]

We probed the dynamical mass profiles of ten galaxy clusters from the HIghest X-ray FLUx Galaxy Cluster Sample (HIFLUGCS) using galaxy kinematics. We numerically solved the spherical Jeans equation and parameterize the dynamical mass profile and the galaxy velocity anisotropy profile using two general functions to ensure that our results are not biased toward any specific model. The mass-velocity anisotropy degeneracy is ameliorated by using two "virial shape parameters" that depend on the fourth moment of velocity distribution. The resulting velocity anisotropy estimates consistently show a nearly isotropic distribution in the inner regions, with an increasing radial anisotropy toward large radii. We compared our derived dynamical masses with those calculated from X-ray gas data assuming hydrostatic equilibrium, finding that massive and rich relaxed clusters generally present consistent mass measurements, while unrelaxed or low-richness clusters have systematically larger total masses than hydrostatic masses by, on average, 50%. This might help alleviate current tensions in the measurement of σ8, but it also leads to cluster baryon fractions below the cosmic value. Finally, our approach probes accelerations as low as 10−11 m s−2, comparable to the outskirts of individual late-type galaxies. We confirm that galaxy clusters deviate from the radial acceleration relation defined by galaxies. [ABSTRACT FROM AUTHOR]

We have investigated the cosmological constraints that can be expected from measurement of the cross-correlation of galaxies with cosmic voids identified in the Euclid spectroscopic survey, which will include spectroscopic information for tens of millions of galaxies over 15 000 deg2 of the sky in the redshift range 0.9 ≤ z < 1.8. We have done this using simulated measurements obtained from the Flagship mock catalogue, the official Euclid mock that closely matches the expected properties of the spectroscopic dataset. To mitigate anisotropic selection-bias effects, we have used a velocity field reconstruction method to remove large-scale redshift-space distortions from the galaxy field before void-finding. This allowed us to accurately model contributions to the observed anisotropy of the cross-correlation function arising from galaxy velocities around voids as well as from the Alcock–Paczynski effect, and we studied the dependence of constraints on the efficiency of reconstruction. We find that Euclid voids will be able to constrain the ratio of the transverse comoving distance DM and Hubble distance DH to a relative precision of about 0.3%, and the growth rate fσ8 to a precision of between 5% and 8% in each of the four redshift bins covering the full redshift range. In the standard cosmological model, this translates to a statistical uncertainty ΔΩm = ±0.0028 on the matter density parameter from voids, which is better than what can be achieved from either Euclid galaxy clustering and weak lensing individually. We also find that voids alone can measure the dark energy equation of state to a 6% precision. [ABSTRACT FROM AUTHOR]

Context. The Fred Young Submillimeter Telescope (FYST) and the Simons Observatory Large Aperture Telescope (SO LAT) will deliver unprecedented high-resolution measurements of microwave sky emissions. Notably, one of those microwave sky emissions, the thermal Sunyaev–Zeldovich (tSZ) signal, is an essential probe for cluster astrophysics and cosmology. However, an obstacle to its measurement is contamination by the cosmic infrared background (CIB), especially at high frequencies. Aims. Our goal is to assess the detection and purity of tSZ power spectrum measurements from these two telescopes. We demonstrate that FYST's high-frequency coverage helps lower CIB contamination and improves signal detection. Methods. We simulated the various components of the microwave sky at the frequencies, sensitivities, and beam sizes of the upcoming SO LAT and FYST telescopes using full-sky Hierarchical Equal Area isoLatitude Pixelisation (HEALPix) map templates from the Websky simulations and the Python Sky Model (PySM). We used a map-based internal linear combination (ILC) and a constrained ILC (CILC) to extract the tSZ signal and compute residual noises to assess CIB contamination and signal recovery. Results. We find that the CIB's residual noise power spectrum in the ILC-recovered tSZ is lowered by ∼35% on average over the scales ℓ ∈ [500, 5000] when SO LAT and FYST are combined compared to when SO LAT is used alone. We find that when using CILC to deproject CIB, the combined abilities of SO LAT and FYST offer a large ℓ ∈ [1800, 3500] window in which the recovered tSZ power spectrum is not noise dominated. [ABSTRACT FROM AUTHOR]

Context. The observation of primordial B modes in cosmic microwave background (CMB) polarisation data represents the main scientific goal of most of the future CMB experiments. This signal is predicted to be much lower than polarised Galactic emission (foregrounds) in any region of the sky, pointing to the need for effective component separation methods. Aims. Among all the techniques, the blind Needlet Internal Linear Combination (NILC) is of great relevance given our current limited knowledge of the B-mode foregrounds. In this work, we explore the possibility of employing NILC for the analysis of B modes reconstructed from partial-sky data, specifically addressing the complications that such an application yields such as E–B leakage, needlet filtering, and beam convolution. Methods. We consider two complementary simulated datasets of future experiments: the balloon-borne Short Wavelength Instrument for the Polarisation Explorer (SWIPE) of the Large Scale Polarisation Explorer, which targets the observation of both reionisation and recombination peaks of the primordial CMB B-mode angular power spectrum, and the ground-based Small Aperture Telescope of Simons Observatory, which, instead, is designed to observe only the recombination bump at ℓ ∼ 80. We assessed the performance of the following two alternative techniques to correct for the CMB E–B leakage: the recycling technique and the Zhao-Baskaran method. Results. We find that both techniques reduce the E–B leakage residuals at a negligible level given the sensitivity of the considered experiments, except for the recycling method in the SWIPE footprint at ℓ < 20. Thus, we implemented two extensions of the pipeline, the iterative B decomposition and the diffusive inpainting, which enabled us to recover the input CMB B-mode power for ℓ ≥ 5. For the considered experiments, we demonstrate that needlet filtering and beam convolution do not affect the CMB B-mode reconstruction. Finally, with an appropriate masking strategy, we find that NILC foregrounds subtraction allows one to achieve sensitivities on the tensor-to-scalar ratio in agreement with the targets of the considered CMB experiments. [ABSTRACT FROM AUTHOR]

Context. Pressure profiles are sensitive probes of the thermodynamic conditions and the internal structure of galaxy clusters. The intra-cluster gas resides in hydrostatic equilibrium within the dark-matter gravitational potential. However, this equilibrium may be perturbed; for example, as a consequence of thermal energy losses, feedback, and non-thermal pressure supports. Accurate measures of the gas pressure over cosmic time are crucial for constraining cluster evolution as well as the contributions from astrophysical processes. Aims. In this work we present a novel algorithm for deriving the pressure profiles of galaxy clusters from the Sunyaev-Zeldovich (SZ) signal measured on a combination of Planck and South Pole Telescope (SPT) observations. The synergy of the two instruments makes it possible to track the profiles on a wide range of spatial scales. We exploited the sensitivity of the Planck High-Frequency Instrument to the larger scales in order to observe the faint peripheries, and took advantage of the higher spatial resolution of SPT to solve the innermost regions. Methods. We developed a two-step pipeline to take advantage of the specifications of each instrument. We first performed a component separation on the two data sets separately in order to remove the background (CMB) and foreground (Galactic emission) contaminants. We then jointly fitted a parametric pressure profile model on a combination of Planck and SPT data. Results. We validated our technique on a sample of six CHEX-MATE clusters detected by SPT. We compare the results of the SZ analysis with profiles derived from X-ray observations with XMM-Newton. We find excellent agreement between these two independent probes of the gas pressure structure. [ABSTRACT FROM AUTHOR]

Context. Weak lensing mass-mapping is a useful tool for accessing the full distribution of dark matter on the sky, but because of intrinsic galaxy ellipticies, finite fields, and missing data, the recovery of dark matter maps constitutes a challenging, ill-posed inverse problem. Aims. We introduce a novel methodology that enables the effcient sampling of the high-dimensional Bayesian posterior of the weak lensing mass-mapping problem, relying on simulations to define a fully non-Gaussian prior. We aim to demonstrate the accuracy of the method to simulated fields, and then proceed to apply it to the mass reconstruction of the HST/ACS COSMOS field. Methods. The proposed methodology combines elements of Bayesian statistics, analytic theory, and a recent class of deep generative models based on neural score matching. This approach allows us to make full use of analytic cosmological theory to constrain the 2pt statistics of the solution, to understand any differences between this analytic prior and full simulations from cosmological simulations, and to obtain samples from the full Bayesian posterior of the problem for robust uncertainty quantification. Results. We demonstrate the method in the TNG simulations and find that the posterior mean significantly outperfoms previous methods (Kaiser-Squires, Wiener filter, Sparsity priors) both for the root-mean-square error and in terms of the Pearson correlation. We further illustrate the interpretability of the recovered posterior by establishing a close correlation between posterior convergence values and the S/N of the clusters artificially introduced into a field. Finally, we apply the method to the reconstruction of the HST/ACS COSMOS field, which yields the highest-quality convergence map of this field to date. Conclusions. We find the proposed approach to be superior to previous algorithms, scalable, providing uncertainties, and using a fully non-Gaussian prior. [ABSTRACT FROM AUTHOR]

The various Euclid imaging surveys will become a reference for studies of galaxy morphology by delivering imaging over an unprecedented area of 15 000 square degrees with high spatial resolution. In order to understand the capabilities of measuring morphologies from Euclid-detected galaxies and to help implement measurements in the pipeline of the Organisational Unit MER of the Euclid Science Ground Segment, we have conducted the Euclid Morphology Challenge, which we present in two papers. While the companion paper focusses on the analysis of photometry, this paper assesses the accuracy of the parametric galaxy morphology measurements in imaging predicted from within the Euclid Wide Survey. We evaluate the performance of five state-of-the-art surface-brightness-fitting codes, DeepLeGATo, Galapagos-2, Morfometryka, ProFit and SourceXtractor++, on a sample of about 1.5 million simulated galaxies (350 000 above 5σ) resembling reduced observations with the Euclid VIS and NIR instruments. The simulations include analytic Sérsic profiles with one and two components, as well as more realistic galaxies generated with neural networks. We find that, despite some code-specific differences, all methods tend to achieve reliable structural measurements (< 10% scatter on ideal Sérsic simulations) down to an apparent magnitude of about IE = 23 in one component and IE = 21 in two components, which correspond to a signal-to-noise ratio of approximately 1 and 5, respectively. We also show that when tested on non-analytic profiles, the results are typically degraded by a factor of 3, driven by systematics. We conclude that the official Euclid Data Releases will deliver robust structural parameters for at least 400 million galaxies in the Euclid Wide Survey by the end of the mission. We find that a key factor for explaining the different behaviour of the codes at the faint end is the set of adopted priors for the various structural parameters. [ABSTRACT FROM AUTHOR]

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

We investigate the dependence of Lyα emitter (LAE) clustering on Lyα luminosity and connect the clustering properties of ≈L⋆ LAEs with those of much fainter ones, namely, ≈0.04L⋆. We use 1030 LAEs from the MUSE-Wide survey, 679 LAEs from MUSE-Deep, and 367 LAEs from the to-date deepest ever spectroscopic survey, the MUSE Extremely Deep Field. All objects have spectroscopic redshifts of 3 < z < 6 and cover a large dynamic range of Lyα luminosities: 40.15 < log(LLyα/erg s−1) < 43.35. We apply the Adelberger et al. K-estimator as the clustering statistic and fit the measurements with state-of-the-art halo occupation distribution (HOD) models. We find that the large-scale bias factor increases weakly with an increasing line luminosity. For the low-luminosity (log⟨LLyα/[erg s−1]⟩ = 41.22) and intermediate-luminosity (log⟨LLyα/[erg s−1]⟩ = 41.64) LAEs, we compute consistent bias factors blow = 2.43−0.15+0.15 b low = 2. 43 − 0.15 + 0.15 $ b_\mathrm{{low}}=2.43^{+0.15}_{-0.15} $ and binterm. = 2.42−0.09+0.10 b interm. = 2. 42 − 0.09 + 0.10 $ b_\mathrm{{interm.}}=2.42^{+0.10}_{-0.09} $ , whereas for the high-luminosity (log⟨LLyα/[erg s−1]⟩ = 42.34) LAEs we calculated bhigh = 2.65−0.11+0.13 b high = 2. 65 − 0.11 + 0.13 $ b_\mathrm{{high}}=2.65^{+0.13}_{-0.11} $. Consequently, high-luminosity LAEs occupy dark matter halos (DMHs) with typical masses of log(Mh/[h−1 M⊙]) = 11.09−0.09+0.10 log (M h / [ h − 1 M ⊙ ]) = 11. 09 − 0.09 + 0.10 $ \log (M_{\mathrm{h}}/[h^{-1}M_\odot])=11.09^{+0.10}_{-0.09} $ , while low-luminosity LAEs reside in halos of log(Mh/[h−1 M⊙]) = 10.77−0.15+0.13 log (M h / [ h − 1 M ⊙ ]) = 10. 77 − 0.15 + 0.13 $ \log (M_{\mathrm{h}}/[h^{-1}M_\odot])=10.77^{+0.13}_{-0.15} $. The minimum masses to host one central LAE, Mmin, and (on average) one satellite LAE, M1, also vary with Lyα luminosity, growing from log(Mmin/[h−1M⊙]) = 10.3−0.3+0.2 log (M min / [ h − 1 M ⊙ ]) = 10. 3 − 0.3 + 0.2 $ \log (M_\mathrm{{min}}/[h^{-1}M_\odot])=10.3^{+0.2}_{-0.3} $ and log(M1/[h−1 M⊙]) = 11.7−0.2+0.3 log (M 1 / [ h − 1 M ⊙ ]) = 11. 7 − 0.2 + 0.3 $ \log (M_1/[h^{-1}M_\odot])=11.7^{+0.3}_{-0.2} $ to log(Mmin/[h−1 M⊙]) = 10.7−0.3+0.2 log (M min / [ h − 1 M ⊙ ]) = 10. 7 − 0.3 + 0.2 $ \log (M_\mathrm{{min}}/[h^{-1}M_\odot])=10.7^{+0.2}_{-0.3} $ and log(M1/[h−1 M⊙]) = 12.4−0.6+0.4 log (M 1 / [ h − 1 M ⊙ ]) = 12. 4 − 0.6 + 0.4 $ \log (M_1/[h^{-1}M_\odot])=12.4^{+0.4}_{-0.6} $ from low- to high-luminosity samples, respectively. The satellite fractions are ≲10% (≲20%) at 1σ (3σ) confidence level, supporting a scenario in which DMHs typically host one single LAE. We next bisected the three main samples into disjoint subsets to thoroughly explore the dependence of the clustering properties on LLyα. We report a strong (8σ) clustering dependence on Lyα luminosity, not accounting for cosmic variance effects, where the highest luminosity LAE subsample (log(LLyα/erg s−1) ≈ 42.53) clusters more strongly (bhighest = 3.13−0.15+0.08 b highest = 3. 13 − 0.15 + 0.08 $ b_\mathrm{{highest}}=3.13^{+0.08}_{-0.15} $) and resides in more massive DMHs (log(Mh/[h−1M⊙] )= 11.43−0.10+0.04 log (M h / [ h − 1 M ⊙ ]) = 11. 43 − 0.10 + 0.04 $ \log(M_\mathrm{{h}} / [h^{-1}{M}_{\odot}])=11.43^{+0.04}_{-0.10} $) than the lowest luminosity one (log(LLyα/erg s−1) ≈ 40.97), which presents a bias of blowest = 1.79−0.06+0.08 b lowest = 1. 79 − 0.06 + 0.08 $ b_\mathrm{{lowest}}=1.79^{+0.08}_{-0.06} $ and occupies log(Mh/[h−1M⊙]) = 10.00−0.09+0.12 log (M h / [ h − 1 M ⊙ ]) = 10. 00 − 0.09 + 0.12 $ \log(M_\mathrm{{h}} / [h^{-1}{M}_{\odot}])=10.00^{+0.12}_{-0.09} $ halos. We discuss the implications of these results for evolving Lyα luminosity functions, halo mass dependent Lyα escape fractions, and incomplete reionization signatures. [ABSTRACT FROM AUTHOR]

Context. The Copernican principle, the notion that we are not at a special location in the Universe, is one of the cornerstones of modern cosmology. Its violation would invalidate the Friedmann-Lemaître-Robertson-Walker metric, causing a major change in our understanding of the Universe. Thus, it is of fundamental importance to perform observational tests of this principle. Aims. We determine the precision with which future surveys will be able to test the Copernican principle and their ability to detect any possible violations. Methods. We forecast constraints on the inhomogeneous Lemaître-Tolman-Bondi (LTB) model with a cosmological constant Λ, basically a cosmological constant Λ and cold dark matter (CDM) model but endowed with a spherical inhomogeneity. We consider combinations of currently available data and simulated Euclid data, together with external data products, based on both ΛCDM and ΛLTB fiducial models. These constraints are compared to the expectations from the Copernican principle. Results. When considering the ΛCDM fiducial model, we find that Euclid data, in combination with other current and forthcoming surveys, will improve the constraints on the Copernican principle by about 30%, with ±10% variations depending on the observables and scales considered. On the other hand, when considering a ΛLTB fiducial model, we find that future Euclid data, combined with other current and forthcoming datasets, will be able to detect gigaparsec-scale inhomogeneities of contrast −0.1. Conclusions. Next-generation surveys, such as Euclid, will thoroughly test homogeneity at large scales, tightening the constraints on possible violations of the Copernican principle. [ABSTRACT FROM AUTHOR]

Context. The 21 cm hydrogen line is arguably one of the most powerful probes with which to explore the Universe, from recombination to the present times. To recover it, it is essential to separate the cosmological signal from the much stronger foreground contributions at radio frequencies. The Baryon Acoustic Oscillations from Integrated Neutral Gas Observations (BINGO) radio telescope is designed to measure the 21 cm line and detect baryon acoustic oscillations (BAOs) using the intensity mapping (IM) technique. Aims. This work analyses the performance of the Generalized Needlet Internal Linear Combination (GNILC) method when combined with a power spectrum debiasing procedure. This method was applied to a simulated BINGO mission, building upon previous work from the collaboration. It compares two different synchrotron emission models and different instrumental configurations and takes into account ancillary data in order to optimize both the removal of foreground emission and the recovery of the 21 cm signal across the full BINGO frequency band and to determine an optimal number of frequency (redshift) bands for the signal recovery. Methods. We produced foreground emission maps using the Planck Sky Model (PSM) and generated cosmological HI emission maps using the Full-Sky Log-normal Astro-Fields simulation Kit (FLASK) package. We also created thermal noise maps according to the instrumental setup. We apply the GNILC method to the simulated sky maps to separate the HI plus thermal noise contribution and, through a debiasing procedure, recover an estimate of the noiseless 21 cm power spectrum. Results. We find a near-optimal reconstruction of the HI signal using an 80-bin configuration, which resulted in a power-spectrum reconstruction average error over all frequencies of 3%. Furthermore, our tests show that GNILC is robust against different synchrotron emission models. Finally, adding an extra channel with C-Band All-Sky Survey (CBASS) foregrounds information, we reduced the estimation error of the 21 cm signal. Conclusions. The optimization of our previous work, producing a configuration with an optimal number of channels for binning the data, significantly impacts decisions regarding BINGO hardware configuration before commissioning. We were able to recover the HI signal with good efficiency in the harmonic space, but have yet to investigate the effect of 1/f noise in the data, which will possibly impact the recovery of the HI signal. This issue will be addressed in forthcoming work. [ABSTRACT FROM AUTHOR]

Having two "sibling" Type Ia supernovae (SNe Ia) in the same galaxy offers additional advantages in reducing a variety of systematic errors involved in estimating the Hubble constant, H0. NGC 4414 is a nearby galaxy included in the Hubble Space Telescope Key Project to measure its distance using Cepheid variables. It hosts two sibling SNe Ia: SN 2021J and SN 1974G. This provides the opportunity to improve the precision of the previous estimate of H0, which was based solely on SN 1974G. Here we present new optical BVRI photometry obtained at the Observatorio de Sierra Nevada and complement it with Swift UVOT UBV data, which cover the first 70 days of emission of SN 2021J. A first look at SN 2021J optical spectra obtained with the Gran Telescopio Canarias (GTC) reveals typical SN type Ia features. The main SN luminosity parameters for the two sibling SNe are obtained by using SNooPy, a light curve fitting code based on templates. Using a hierarchical bayesian approach, we build the Hubble diagram with a sample of 96 SNe Ia obtained from the Combined Pantheon Sample in the redshift range z = 0.02 − 0.075, and calibrate the zero point with the two sibling type-Ia SNe in NGC 4414. We report a value of the Hubble constant H0 = 72.19 ± 2.32 (stat.) ±3.42 (syst.) km s−1 Mpc−1. We expect a reduction of the systematic error after a new analysis of the Cepheids period-luminosity relation using the upcoming Gaia DR4 and additional Cepheids from the HST and JWST. [ABSTRACT FROM AUTHOR]

We present a set of maps classifying regions of the sky according to their information gain potential as quantified by Fisher information. These maps can guide the optimal retrieval of relevant physical information with targeted cosmological searches. Specifically, we calculated the response of observed cosmic structures to perturbative changes in the cosmological model and we charted their respective contributions to Fisher information. Our physical forward-modeling machinery transcends the limitations of contemporary analyses based on statistical summaries to yield detailed characterizations of individual 3D structures. We demonstrate this advantage using galaxy counts data and we showcase the potential of our approach by studying the information gain of the Coma cluster. We find that regions in the vicinity of the filaments and cluster core, where mass accretion ensues from gravitational infall, are the most informative with regard to our physical model of structure formation in the Universe. Hence, collecting data in those regions would be most optimal for testing our model predictions. The results presented in this work are the first of their kind to elucidate the inhomogeneous distribution of cosmological information in the Universe. This study paves a new way forward for the performance of efficient targeted searches for the fundamental physics of the Universe, where search strategies are progressively refined with new cosmological data sets within an active learning framework. [ABSTRACT FROM AUTHOR]

Aims. We study a sample of Hβ emission line sources at z ∼ 0.9 to identify the star-forming galaxies sample and characterise them in terms of line luminosity, stellar mass, star formation rate, and morphology. The final aim is to obtain the Hβ luminosity function of the star-forming galaxies at this redshift. Methods. We used the red tunable filter of the instrument Optical System for Imaging low Resolution Integrated Spectroscopy (OSIRIS) at Gran Telescopio de Canarias to obtain the pseudo spectra of emission line sources in the OTELO field. From these pseudo spectra, we identified the objects with Hβ emission. As the resolution of the pseudo spectra allowed us to separate Hβ from [O III], we were able to derive the Hβ flux without contamination from its adjacent line. Using data from the extended OTELO catalogue, we discriminated AGNs and studied the star formation rate, the stellar mass, and the morphology of the star-forming galaxies. Results. We find that our sample is located on the main sequence of star-forming galaxies. The sources are morphologically classified, mostly as disc-like galaxies (76%), and 90% of the sample are low-mass galaxies (M* < 1010 M⊙). The low-mass star-forming galaxies at z ∼ 0.9 that were detected by OTELO present similar properties as low-mass star-forming galaxies in the local universe, suggesting that these kinds of objects do not have a favorite epoch of formation and star formation enhancement from z ∼ 1 to now. Our sample of 40 Hβ star-forming galaxies include the faintest Hβ emitters detected so far. This allows us to constrain the faint end of the luminosity function for the Hβ line alone with a minimum luminosity of log L = 39 erg s−1, which is a hundred times fainter than previous surveys. The dust-corrected OSIRIS Tunable Emission Line Object survey (OTELO) Hβ luminosity function established the faint-end slope as α = −1.36 ± 0.15. We increased the scope of the analysis to the bright end by adding ancillary data from the literature, which was not dust-corrected in this case. The obtained slope for this extended luminosity function is α = −1.43 ± 0.12. [ABSTRACT FROM AUTHOR]

The Javalambre-Physics of the Accelerating Universe Astrophysical Survey (J-PAS) will scan thousands of square degrees of the northern sky with a unique set of 56 filters using the dedicated 2.55 m Javalambre Survey Telescope (JST) at the Javalambre Astrophysical Observatory. Prior to the installation of the main camera (4.2 deg2 field-of-view with 1.2 Gpixels), the JST was equipped with the JPAS-Pathfinder, a one CCD camera with a 0.3 deg2 field-of-view and plate scale of 0.23 arcsec pixel−1. To demonstrate the scientific potential of J-PAS, the JPAS-Pathfinder camera was used to perform miniJPAS, a ∼1 deg2 survey of the AEGIS field (along the Extended Groth Strip). The field was observed with the 56 J-PAS filters, which include 54 narrow band (FWHM ∼ 145 Å) and two broader filters extending to the UV and the near-infrared, complemented by the u, g, r, i SDSS broad band filters. In this miniJPAS survey overview paper, we present the miniJPAS data set (images and catalogs), as we highlight key aspects and applications of these unique spectro-photometric data and describe how to access the public data products. The data parameters reach depths of magAB ≃ 22−23.5 in the 54 narrow band filters and up to 24 in the broader filters (5σ in a 3″ aperture). The miniJPAS primary catalog contains more than 64 000 sources detected in the r band and with matched photometry in all other bands. This catalog is 99% complete at r = 23.6 (r = 22.7) mag for point-like (extended) sources. We show that our photometric redshifts have an accuracy better than 1% for all sources up to r = 22.5, and a precision of ≤0.3% for a subset consisting of about half of the sample. On this basis, we outline several scientific applications of our data, including the study of spatially-resolved stellar populations of nearby galaxies, the analysis of the large scale structure up to z ∼ 0.9, and the detection of large numbers of clusters and groups. Sub-percent redshift precision can also be reached for quasars, allowing for the study of the large-scale structure to be pushed to z > 2. The miniJPAS survey demonstrates the capability of the J-PAS filter system to accurately characterize a broad variety of sources and paves the way for the upcoming arrival of J-PAS, which will multiply this data by three orders of magnitude. [ABSTRACT FROM AUTHOR]

We performed an internal-consistency test of the KiDS+VIKING-450 (KV450) cosmic shear analysis with a colour-based split of source galaxies. Utilising the same measurements and calibrations for both sub-samples, we inspected the characteristics of the shear measurements and the performance of the calibration pipelines. On the modelling side, we examined the observational nuisance parameters, specifically those for the redshift calibration and intrinsic alignments, using a Bayesian analysis with dedicated test parameters. We verified that the current nuisance parameters are sufficient for the KV450 data to capture residual systematics, with slight deviations seen in the second and the third redshift tomographic bins. Our test also showcases the degeneracy between the inferred amplitude of intrinsic alignments and the redshift uncertainties in low redshift tomographic bins. The test is rather insensitive to the background cosmology and, therefore, can be implemented before any cosmological inference is made. [ABSTRACT FROM AUTHOR]

The largest temperature anisotropy in the cosmic microwave background (CMB) is the dipole, which has been measured with increasing accuracy for more than three decades, particularly with the Planck satellite. The simplest interpretation of the dipole is that it is due to our motion with respect to the rest frame of the CMB. Since current CMB experiments infer temperature anisotropies from angular intensity variations, the dipole modulates the temperature anisotropies with the same frequency dependence as the thermal Sunyaev-Zeldovich (tSZ) effect. We present the first, and significant, detection of this signal in the tSZ maps and find that it is consistent with direct measurements of the CMB dipole, as expected. The signal contributes power in the tSZ maps, which is modulated in a quadrupolar pattern, and we estimate its contribution to the tSZ bispectrum, noting that it contributes negligible noise to the bispectrum at relevant scales. [ABSTRACT FROM AUTHOR]

We describe an extension of the most recent version of the Planck Catalogue of Compact Sources (PCCS2), produced using a new multi-band Bayesian Extraction and Estimation Package (BeeP). BeeP assumes that the compact sources present in PCCS2 at 857 GHz have a dust-like spectral energy distribution (SED), which leads to emission at both lower and higher frequencies, and adjusts the parameters of the source and its SED to fit the emission observed in Planck's three highest frequency channels at 353, 545, and 857 GHz, as well as the IRIS map at 3000 GHz. In order to reduce confusion regarding diffuse cirrus emission, BeeP's data model includes a description of the background emission surrounding each source, and it adjusts the confidence in the source parameter extraction based on the statistical properties of the spatial distribution of the background emission. BeeP produces the following three new sets of parameters for each source: (a) fits to a modified blackbody (MBB) thermal emission model of the source; (b) SED-independent source flux densities at each frequency considered; and (c) fits to an MBB model of the background in which the source is embedded. BeeP also calculates, for each source, a reliability parameter, which takes into account confusion due to the surrounding cirrus. This parameter can be used to extract sub-samples of high-frequency sources with statistically well-understood properties. We define a high-reliability subset (BeeP/base), containing 26 083 sources (54.1% of the total PCCS2 catalogue), the majority of which have no information on reliability in the PCCS2. We describe the characteristics of this specific high-quality subset of PCCS2 and its validation against other data sets, specifically for: the sub-sample of PCCS2 located in low-cirrus areas; the Planck Catalogue of Galactic Cold Clumps; the Herschel GAMA15-field catalogue; and the temperature- and spectral-index-reconstructed dust maps obtained with Planck's Generalized Needlet Internal Linear Combination method. The results of the BeeP extension of PCCS2, which are made publicly available via the Planck Legacy Archive, will enable the study of the thermal properties of well-defined samples of compact Galactic and extragalactic dusty sources. [ABSTRACT FROM AUTHOR]

As part of an on-going effort to identify, understand and correct for astrophysics biases in the standardization of Type Ia supernovae (SN Ia) for cosmology, we have statistically classified a large sample of nearby SNe Ia into those that are located in predominantly younger or older environments. This classification is based on the specific star formation rate measured within a projected distance of 1 kpc from each SN location (LsSFR). This is an important refinement compared to using the local star formation rate directly, as it provides a normalization for relative numbers of available SN progenitors and is more robust against extinction by dust. We find that the SNe Ia in predominantly younger environments are ΔY = 0.163 ± 0.029 mag (5.7σ) fainter than those in predominantly older environments after conventional light-curve standardization. This is the strongest standardized SN Ia brightness systematic connected to the host-galaxy environment measured to date. The well-established step in standardized brightnesses between SNe Ia in hosts with lower or higher total stellar masses is smaller, at ΔM = 0.119 ± 0.032 mag (4.5σ), for the same set of SNe Ia. When fit simultaneously, the environment-age offset remains very significant, with ΔY = 0.129 ± 0.032 mag (4.0σ), while the global stellar mass step is reduced to ΔM = 0.064  ±  0.029 mag (2.2σ). Thus, approximately 70% of the variance from the stellar mass step is due to an underlying dependence on environment-based progenitor age. Also, we verify that using the local star formation rate alone is not as powerful as LsSFR at sorting SNe Ia into brighter and fainter subsets. Standardization that only uses the SNe Ia in younger environments reduces the total dispersion from 0.142  ±  0.008 mag to 0.120  ±  0.010 mag. We show that as environment-ages evolve with redshift, a strong bias, especially on the measurement of the derivative of the dark energy equation of state, can develop. Fortunately, data that measure and correct for this effect using our local specific star formation rate indicator, are likely to be available for many next-generation SN Ia cosmology experiments. [ABSTRACT FROM AUTHOR]

We review and compare two different cosmic microwave background (CMB) dipole estimators discussed in the literature and assess their performances through Monte Carlo simulations. The first method amounts to simple template regression with partial sky data, while the second method is an optimal Wiener filter (or Gibbs sampling) implementation. The main difference between the two methods is that the latter approach takes into account correlations with higher-order CMB temperature fluctuations that arise from nonorthogonal spherical harmonics on an incomplete sky, which for recent CMB data sets (such as Planck) is the dominant source of uncertainty. For an accepted sky fraction of 81% and an angular CMB power spectrum corresponding to the best-fit Planck 2018 ΛCDM model, we find that the uncertainty on the recovered dipole amplitude is about six times smaller for the Wiener filter approach than for the template approach, corresponding to 0.5 and 3 μK, respectively. Similar relative differences are found for the corresponding directional parameters and other sky fractions. We note that the Wiener filter algorithm is generally applicable to any dipole estimation problem on an incomplete sky, as long as a statistical and computationally tractable model is available for the unmasked higher-order fluctuations. The methodology described in this paper forms the numerical basis for the most recent determination of the CMB solar dipole from Planck, as summarized by Planck Collaboration Int. LVII (2020). [ABSTRACT FROM AUTHOR]

The H0LiCOW collaboration inferred via strong gravitational lensing time delays a Hubble constant value of H0 = 73.3−1.8+1.7 H 0 = 73. 3 − 1.8 + 1.7 $ H_0 = 73.3^{+1.7}_{-1.8} $ km s−1 Mpc−1, describing deflector mass density profiles by either a power-law or stars (constant mass-to-light ratio) plus standard dark matter halos. The mass-sheet transform (MST) that leaves the lensing observables unchanged is considered the dominant source of residual uncertainty in H0. We quantify any potential effect of the MST with a flexible family of mass models, which directly encodes it, and they are hence maximally degenerate with H0. Our calculation is based on a new hierarchical Bayesian approach in which the MST is only constrained by stellar kinematics. The approach is validated on mock lenses, which are generated from hydrodynamic simulations. We first applied the inference to the TDCOSMO sample of seven lenses, six of which are from H0LiCOW, and measured H0 = 74.5−6.1+5.6 H 0 = 74. 5 − 6.1 + 5.6 $ H_0 = 74.5^{+5.6}_{-6.1} $ km s−1 Mpc−1. Secondly, in order to further constrain the deflector mass density profiles, we added imaging and spectroscopy for a set of 33 strong gravitational lenses from the Sloan Lens ACS (SLACS) sample. For nine of the 33 SLAC lenses, we used resolved kinematics to constrain the stellar anisotropy. From the joint hierarchical analysis of the TDCOSMO+SLACS sample, we measured H0 = 67.4−3.2+4.1 H 0 = 67. 4 − 3.2 + 4.1 $ H_0 = 67.4^{+4.1}_{-3.2} $ km s−1 Mpc−1. This measurement assumes that the TDCOSMO and SLACS galaxies are drawn from the same parent population. The blind H0LiCOW, TDCOSMO-only and TDCOSMO+SLACS analyses are in mutual statistical agreement. The TDCOSMO+SLACS analysis prefers marginally shallower mass profiles than H0LiCOW or TDCOSMO-only. Without relying on the form of the mass density profile used by H0LiCOW, we achieve a ∼5% measurement of H0. While our new hierarchical analysis does not statistically invalidate the mass profile assumptions by H0LiCOW – and thus the H0 measurement relying on them – it demonstrates the importance of understanding the mass density profile of elliptical galaxies. The uncertainties on H0 derived in this paper can be reduced by physical or observational priors on the form of the mass profile, or by additional data. [ABSTRACT FROM AUTHOR]

We present a new shear calibration method based on machine learning. The method estimates the individual shear responses of the objects from the combination of several measured properties on the images using supervised learning. The supervised learning uses the true individual shear responses obtained from copies of the image simulations with different shear values. On simulated GREAT3 data, we obtain a residual bias after the calibration compatible with 0 and beyond Euclid requirements for a signal-to-noise ratio > 20 within ∼15 CPU hours of training using only ∼105 objects. This efficient machine-learning approach can use a smaller data set because the method avoids the contribution from shape noise. The low dimensionality of the input data also leads to simple neural network architectures. We compare it to the recently described method Metacalibration, which shows similar performances. The different methods and systematics suggest that the two methods are very good complementary methods. Our method can therefore be applied without much effort to any survey such as Euclid or the Vera C. Rubin Observatory, with fewer than a million images to simulate to learn the calibration function. [ABSTRACT FROM AUTHOR]

We report the first statistical detection of X-ray emission from cosmic web filaments in ROSAT data. We selected 15 165 filaments at 0.2 <  z <  0.6 ranging from 30 Mpc to 100 Mpc in length, identified in the Sloan Digital Sky Survey survey. We stacked the X-ray count-rate maps from ROSAT around the filaments, excluding resolved galaxy groups and clusters above the mass of ∼3 × 1013 M⊙ as well as the detected X-ray point sources from the ROSAT, Chandra, and XMM-Newton observations. The stacked signal results in the detection of the X-ray emission from the cosmic filaments at a significance of 4.2σ in the energy band of 0.56−1.21 keV. The signal is interpreted, assuming the Astrophysical Plasma Emission Code model, as an emission from the hot gas in the filament-core regions with an average gas temperature of 0.9−0.6+1.0 0. 9 − 0.6 + 1.0 $ 0.9^{+1.0}_{-0.6} $ keV and a gas overdensity of δ ∼ 30 at the center of the filaments. Furthermore, we show that stacking the SRG/eROSITA data for ∼2000 filaments only would lead to a ≳5σ detection of their X-ray signal, even with an average gas temperature as low as ∼0.3 keV. [ABSTRACT FROM AUTHOR]

Aims. We present the detailed characterisation of a sample of 56 sources serendipitously detected in ALMA band 7 as part of the ALMA Large Program to INvestigate CII at Early Times (ALPINE). These sources, detected in COSMOS and ECDFS, have been used to derive the total infrared luminosity function (LF) and to estimate the cosmic star formation rate density (SFRD) up to z ≃ 6. Methods. We looked for counterparts of the ALMA sources in all the available multi-wavelength (from HST to VLA) and photometric redshift catalogues. We also made use of deeper UltraVISTA and Spitzer source lists and maps to identify optically dark sources with no matches in the public catalogues. We used the sources with estimated redshifts to derive the 250 μm rest-frame and total infrared (8–1000 μm) LFs from z ≃ 0.5 to 6. Results. Our ALMA blind survey (860 μm flux density range: ∼0.3–12.5 mJy) allows us to further push the study of the nature and evolution of dusty galaxies at high-z, identifying luminous and massive sources to redshifts and faint luminosities never probed before by any far-infrared surveys. The ALPINE data are the first ones to sample the faint end of the infrared LF, showing little evolution from z ≃ 2.5 to z ≃ 6, and a "flat" slope up to the highest redshifts (i.e. 4.5 <  z <  6). The SFRD obtained by integrating the luminosity function remains almost constant between z ≃ 2 and z ≃ 6, and significantly higher than the optical or ultra-violet derivations, showing a significant contribution of dusty galaxies and obscured star formation at high-z. About 14% of all the ALPINE serendipitous continuum sources are found to be optically and near-infrared (near-IR) dark (to a depth Ks ∼ 24.9 mag). Six show a counterpart only in the mid-IR and no HST or near-IR identification, while two are detected as [C II] emitters at z ≃ 5. The six HST+near-IR dark galaxies with mid-IR counterparts are found to contribute about 17% of the total SFRD at z ≃ 5 and to dominate the high-mass end of the stellar mass function at z >  3. [ABSTRACT FROM AUTHOR]

We present the NPIPE processing pipeline, which produces calibrated frequency maps in temperature and polarization from data from the Planck Low Frequency Instrument (LFI) and High Frequency Instrument (HFI) using high-performance computers. NPIPE represents a natural evolution of previous Planck analysis efforts, and combines some of the most powerful features of the separate LFI and HFI analysis pipelines. For example, following the LFI 2018 processing procedure, NPIPE uses foreground polarization priors during the calibration stage in order to break scanning-induced degeneracies. Similarly, NPIPE employs the HFI 2018 time-domain processing methodology to correct for bandpass mismatch at all frequencies. In addition, NPIPE introduces several improvements, including, but not limited to: inclusion of the 8% of data collected during repointing manoeuvres; smoothing of the LFI reference load data streams; in-flight estimation of detector polarization parameters; and construction of maximally independent detector-set split maps. For component-separation purposes, important improvements include: maps that retain the CMB Solar dipole, allowing for high-precision relative calibration in higher-level analyses; well-defined single-detector maps, allowing for robust CO extraction; and HFI temperature maps between 217 and 857 GHz that are binned into 0′.9 pixels (Nside = 4096), ensuring that the full angular information in the data is represented in the maps even at the highest Planck resolutions. The net effect of these improvements is lower levels of noise and systematics in both frequency and component maps at essentially all angular scales, as well as notably improved internal consistency between the various frequency channels. Based on the NPIPE maps, we present the first estimate of the Solar dipole determined through component separation across all nine Planck frequencies. The amplitude is (3366.6 ± 2.7) μK, consistent with, albeit slightly higher than, earlier estimates. From the large-scale polarization data, we derive an updated estimate of the optical depth of reionization of τ = 0.051 ± 0.006, which appears robust with respect to data and sky cuts. There are 600 complete signal, noise and systematics simulations of the full-frequency and detector-set maps. As a Planck first, these simulations include full time-domain processing of the beam-convolved CMB anisotropies. The release of NPIPE maps and simulations is accompanied with a complete suite of raw and processed time-ordered data and the software, scripts, auxiliary data, and parameter files needed to improve further on the analysis and to run matching simulations. [ABSTRACT FROM AUTHOR]

Context. The LOFAR Two-metre Sky Survey (LoTSS) will eventually map the complete Northern sky and provide an excellent opportunity to study the distribution and evolution of the large-scale structure of the Universe. Aims. We test the quality of LoTSS observations through a statistical comparison of the LoTSS first data release (DR1) catalogues to expectations from the established cosmological model of a statistically isotropic and homogeneous Universe. Methods. We study the point-source completeness and define several quality cuts, in order to determine the count-in-cell statistics and differential source count statistics, and measure the angular two-point correlation function. We use the photometric redshift estimates, which are available for about half of the LoTSS-DR1 radio sources, to compare the clustering throughout the history of the Universe. Results. For the masked LoTSS-DR1 value-added source catalogue, we find a point-source completeness of 99% above flux densities of 0.8 mJy. The counts-in-cell statistic reveals that the distribution of radio sources cannot be described by a spatial Poisson process. Instead, a good fit is provided by a compound Poisson distribution. The differential source counts are in good agreement with previous findings in deep fields at low radio frequencies and with simulated catalogues from the SKA Design Study and the Tiered Radio Extragalactic Continuum Simulation. Restricting the value added source catalogue to low-noise regions and applying a flux density threshold of 2 mJy provides our most reliable estimate of the angular two-point correlation. Based on the distribution of photometric redshifts and the Planck 2018 best-fit cosmological model, the theoretically predicted angular two-point correlation between 0.1 deg and 6 deg agrees reasonably well with the measured clustering for the sub-sample of radio sources with redshift information. Conclusions. The deviation from a Poissonian distribution might be a consequence of the multi-component nature of a large number of resolved radio sources and/or of uncertainties on the flux density calibration. The angular two-point correlation function is < 10−2 at angular scales > 1 deg and up to the largest scales probed. At a 2 mJy flux density threshold and at a pivot angle of 1 deg, we find a clustering amplitude of A = (5.1 ± 0.6) × 10−3 with a slope parameter of γ = 0.74 ± 0.16. For smaller flux density thresholds, systematic issues are identified, which are most likely related to the flux density calibration of the individual pointings. We conclude that we find agreement with the expectation of large-scale statistical isotropy of the radio sky at the per cent level. The angular two-point correlation agrees well with the expectation of the cosmological standard model. [ABSTRACT FROM AUTHOR]

We consider the observational aspects of the value of dark energy density from quantum vacuum fluctuations based initially on the Gurzadyan–Xue model. We reduce the Djorgovski–Gurzadyan integral equation to a differential equation for the co-moving horizon and then, by means of the obtained explicit form for the luminosity distance, we construct the Hubble diagram for two classes of observational samples. For supernova and gamma-ray burst data we show that this approach provides viable predictions for distances up to z ≃ 9, quantitatively at least as good as those provided by the Λ cold dark matter model. The Hubble parameter dependence H(z) of the two models also reveals mutual crossing at z = 0.4018, the interpretation of which is less evident. [ABSTRACT FROM AUTHOR]

Measuring cosmic shear in wide-field imaging surveys requires accurate knowledge of the redshift distribution of all sources. The clustering-redshift technique exploits the angular cross-correlation of a target galaxy sample with unknown redshifts and a reference sample with known redshifts. It represents an attractive alternative to colour-based methods of redshift calibration. Here we test the performance of such clustering redshift measurements using mock catalogues that resemble the Kilo-Degree Survey (KiDS). These mocks are created from the MICE simulation and closely mimic the properties of the KiDS source sample and the overlapping spectroscopic reference samples. We quantify the performance of the clustering redshifts by comparing the cross-correlation results with the true redshift distributions in each of the five KiDS photometric redshift bins. Such a comparison to an informative model is necessary due to the incompleteness of the reference samples at high redshifts. Clustering mean redshifts are unbiased at |Δz|< 0.006 under these conditions. The redshift evolution of the galaxy bias of the reference and target samples represents one of the most important systematic errors when estimating clustering redshifts. It can be reliably mitigated at this level of precision using auto-correlation measurements and self-consistency relations, and will not become a dominant source of systematic error until the arrival of Stage-IV cosmic shear surveys. Using redshift distributions from a direct colour-based estimate instead of the true redshift distributions as a model for comparison with the clustering redshifts increases the biases in the mean to up to |Δz|∼0.04. This indicates that the interpretation of clustering redshifts in real-world applications will require more sophisticated (parameterised) models of the redshift distribution in the future. If such better models are available, the clustering-redshift technique promises to be a highly complementary alternative to other methods of redshift calibration. [ABSTRACT FROM AUTHOR]

Context. Characterising the large-scale structure in the Universe from present times to the high redshift epoch of reionisation is essential to constraining the cosmology, the history of star formation, and reionisation, to measuring the gas content of the Universe, and to obtaining a better understanding of the physical processes that drive galaxy formation and evolution. Using the integrated emission from unresolved galaxies or gas clouds, line intensity mapping (LIM) provides a new observational window to measure the larger properties of structures. This very promising technique motivates the community to plan for LIM experiments. Aims. We describe the development of a large field-of-view instrument, named CONCERTO (for CarbON CII line in post-rEionisation and ReionisaTiOn epoch), operating in the range 130–310 GHz from the APEX 12-m telescope (5100 m above sea level). CONCERTO is a low-resolution spectrometer based on the lumped element kinetic inductance detectors (LEKID) technology. Spectra are obtained using a fast Fourier transform spectrometer (FTS), coupled to a dilution cryostat with a base temperature of 0.1 K. Two two kilo-pixel arrays of LEKID are mounted inside the cryostat that also contains the cold optics and the front-end electronics. Methods. We present, in detail, the technological choices leading to the instrumental concept, together with the design and fabrication of the instrument and preliminary laboratory tests on the detectors. We also give our best estimates for CONCERTO sensitivity and give predictions for two of the main scientific goals of CONCERTO, that is, a [CII]-intensity mapping survey and observations of galaxy clusters. Results. We provide a detailed description of the instrument design. Based on realistic comparisons with existing instruments developed by our group (NIKA, NIKA2, and KISS), and on the laboratory characterisation of our detectors, we provide an estimate for CONCERTO sensitivity on the sky. Finally, we describe, in detail, two of the main scientific goals offered by CONCERTO at APEX. [ABSTRACT FROM AUTHOR]

The Complete Calibration of the Colour–Redshift Relation survey (C3R2) is a spectroscopic effort involving ESO and Keck facilities designed specifically to empirically calibrate the galaxy colour–redshift relation – P(z|C) to the Euclid depth (iAB = 24.5) and is intimately linked to the success of upcoming Stage IV dark energy missions based on weak lensing cosmology. The aim is to build a spectroscopic calibration sample that is as representative as possible of the galaxies of the Euclid weak lensing sample. In order to minimise the number of spectroscopic observations necessary to fill the gaps in current knowledge of the P(z|C), self-organising map (SOM) representations of the galaxy colour space have been constructed. Here we present the first results of an ESO@VLT Large Programme approved in the context of C3R2, which makes use of the two VLT optical and near-infrared multi-object spectrographs, FORS2 and KMOS. This data release paper focuses on high-quality spectroscopic redshifts of high-redshift galaxies observed with the KMOS spectrograph in the near-infrared H- and K-bands. A total of 424 highly-reliable redshifts are measured in the 1.3 ≤ z ≤ 2.5 range, with total success rates of 60.7% in the H-band and 32.8% in the K-band. The newly determined redshifts fill 55% of high (mainly regions with no spectroscopic measurements) and 35% of lower (regions with low-resolution/low-quality spectroscopic measurements) priority empty SOM grid cells. We measured Hα fluxes in a 1.″2 1. ″ 2 $ 1{{\overset{\prime\prime}{.}}}2 $ radius aperture from the spectra of the spectroscopically confirmed galaxies and converted them into star formation rates. In addition, we performed an SED fitting analysis on the same sample in order to derive stellar masses, E(B − V), total magnitudes, and SFRs. We combine the results obtained from the spectra with those derived via SED fitting, and we show that the spectroscopic failures come from either weakly star-forming galaxies (at z <  1.7, i.e. in the H-band) or low S/N spectra (in the K-band) of z >  2 galaxies. [ABSTRACT FROM AUTHOR]

Aims. The Euclid space telescope will measure the shapes and redshifts of galaxies to reconstruct the expansion history of the Universe and the growth of cosmic structures. The estimation of the expected performance of the experiment, in terms of predicted constraints on cosmological parameters, has so far relied on various individual methodologies and numerical implementations, which were developed for different observational probes and for the combination thereof. In this paper we present validated forecasts, which combine both theoretical and observational ingredients for different cosmological probes. This work is presented to provide the community with reliable numerical codes and methods for Euclid cosmological forecasts. Methods. We describe in detail the methods adopted for Fisher matrix forecasts, which were applied to galaxy clustering, weak lensing, and the combination thereof. We estimated the required accuracy for Euclid forecasts and outline a methodology for their development. We then compare and improve different numerical implementations, reaching uncertainties on the errors of cosmological parameters that are less than the required precision in all cases. Furthermore, we provide details on the validated implementations, some of which are made publicly available, in different programming languages, together with a reference training-set of input and output matrices for a set of specific models. These can be used by the reader to validate their own implementations if required. Results. We present new cosmological forecasts for Euclid. We find that results depend on the specific cosmological model and remaining freedom in each setting, for example flat or non-flat spatial cosmologies, or different cuts at non-linear scales. The numerical implementations are now reliable for these settings. We present the results for an optimistic and a pessimistic choice for these types of settings. We demonstrate that the impact of cross-correlations is particularly relevant for models beyond a cosmological constant and may allow us to increase the dark energy figure of merit by at least a factor of three. [ABSTRACT FROM AUTHOR]

Context. Lyman α blobs (LABs) are large-scale radio-quiet Lyman α (Lyα) nebula at high-z that occur predominantly in overdense proto-cluster regions. In particular, there is the prototypical SSA22a-LAB1 at z = 3.1, which has become an observational reference for LABs across the electromagnetic spectrum. Aims. We want to understand the powering mechanisms that drive the LAB so that we may gain empirical insights into the galaxy-formation processes within a rare dense environment at high-z. Thus, we need to infer the distribution, the dynamics, and the ionisation state of LAB 1's Lyα emitting gas. Methods. LAB 1 was observed for 17.2 h with the VLT/MUSE integral-field spectrograph. We produced optimally extracted narrow band images, in Lyαλ1216, He IIλ1640, and we tried to detect C IVλ1549 emission. By utilising a moment-based analysis, we mapped the kinematics and the line profile characteristics of the blob. We also linked the inferences from the line profile analysis to previous results from imaging polarimetry. Results. We map Lyα emission from the blob down to surface-brightness limits of ≈6 × 10−19 erg s−1 cm−2 arcsec−2. At this depth, we reveal a bridge between LAB 1 and its northern neighbour LAB 8, as well as a shell-like filament towards the south of LAB 1. The complexity and morphology of the Lyα profile vary strongly throughout the blob. Despite the complexity, we find a coherent large-scale east-west velocity gradient of ∼1000 km s−1 that is aligned perpendicular to the major axis of the blob. Moreover, we observe a negative correlation of Lyα polarisation fraction with Lyα line width and a positive correlation with absolute line-of-sight velocity. Finally, we reveal He II emission in three distinct regions within the blob, however, we can only provide upper limits for C IV. Conclusions. Various gas excitation mechanisms are at play in LAB 1: ionising radiation and feedback effects dominate near the embedded galaxies, while Lyα scattering contributes at larger distances. However, He II/Lyα ratios combined with upper limits on C IV/Lyα are not able to discriminate between active galactic nucleus ionisation and feedback- driven shocks. The alignment of the angular momentum vector parallel to the morphological principal axis appears to be at odds with the predicted norm for high-mass halos, but this most likely reflects that LAB 1 resides at a node of multiple intersecting filaments of the cosmic web. LAB 1 can thus be thought of as a progenitor of a present-day massive elliptical within a galaxy cluster. [ABSTRACT FROM AUTHOR]

Context. Type Ia supernovae (SNe Ia) are widely used to measure the expansion of the Universe. Improving distance measurements of SNe Ia is one technique to better constrain the acceleration of expansion and determine its physical nature. Aims. This document develops a new SNe Ia spectral energy distribution (SED) model, called the SUpernova Generator And Reconstructor (SUGAR), which improves the spectral description of SNe Ia, and consequently could improve the distance measurements. Methods. This model was constructed from SNe Ia spectral properties and spectrophotometric data from the Nearby Supernova Factory collaboration. In a first step, a principal component analysis-like method was used on spectral features measured at maximum light, which allowed us to extract the intrinsic properties of SNe Ia. Next, the intrinsic properties were used to extract the average extinction curve. Third, an interpolation using Gaussian processes facilitated using data taken at different epochs during the lifetime of an SN Ia and then projecting the data on a fixed time grid. Finally, the three steps were combined to build the SED model as a function of time and wavelength. This is the SUGAR model. Results. The main advancement in SUGAR is the addition of two additional parameters to characterize SNe Ia variability. The first is tied to the properties of SNe Ia ejecta velocity and the second correlates with their calcium lines. The addition of these parameters, as well as the high quality of the Nearby Supernova Factory data, makes SUGAR an accurate and efficient model for describing the spectra of normal SNe Ia as they brighten and fade. Conclusions. The performance of this model makes it an excellent SED model for experiments like the Zwicky Transient Facility, the Large Synoptic Survey Telescope, or the Wide Field Infrared Survey Telescope. [ABSTRACT FROM AUTHOR]

We present spatially resolved maps of six individually-detected Lyman ff haloes (LAHs) as well as a first statistical analysis of the Lyman ff (Lyff) spectral signature in the circum-galactic medium of high-redshift star-forming galaxies (-17:5 > MUV > -21:5) using the Multi-Unit Spectroscopic Explorer. Our resolved spectroscopic analysis of the LAHs reveals significant intrahalo variations of the Lyff line profile. Using a three-dimensional two-component model for the Lyff emission, we measured the full width at half maximum (FWHM), the peak velocity shift, and the asymmetry of the Lyff line in the core and in the halo of 19 galaxies. We find that the Lyff line shape is statistically different in the halo compared to the core (in terms of width, peak wavelength, and asymmetry) for 40% of our galaxies. Similarly to object-by-object based studies and a recent resolved study using lensing, we find a correlation between the peak velocity shift and the width of the Lyff line both at the interstellar and circum-galactic scales. This trend has been predicted by radiative transfer simulations of galactic winds as a result of resonant scattering in outflows. While there is a lack of correlation between the spectral properties and the spatial scale lengths of our LAHs, we find a correlation between the width of the line in the LAH and the halo flux fraction. Interestingly, UV bright galaxies (MUV < -20) show broader, more redshifted, and less asymmetric Lyff lines in their haloes. The most significant correlation found is for the FWHM of the line and the UV continuum slope of the galaxy, suggesting that the redder galaxies have broader Lyff lines. The generally broad and red line shapes found in the halo component suggest that the Lyff haloes are powered either by scattering processes through an outflowing medium, fluorescent emission from outflowing cold clumps of gas, or a mix of both. Considering the large diversity of the Lyff line profiles observed in our sample and the lack of strong correlation, the interpretation of our results is still broadly open and underlines the need for realistic spatially resolved models of the LAHs. [ABSTRACT FROM AUTHOR]

Mapping the behaviour of dark energy is a pressing task for observational cosmology. Phenomenological classification divides dynamical dark energy models into freezing and thawing, depending on whether the dark energy equation of state is approaching or moving away from w = p/p = -1. Moreover, in realistic dynamical dark energy models the dynamical degree of freedom is expected to couple to the electromagnetic sector, leading to variations of the fine-structure constant α. We discuss the feasibility of distinguishing between the freezing and thawing classes of models with current and forthcoming observational facilities and using a parametrisation of the dark energy equation of state, which can have either behaviour, introduced by Mukhanov as fiducial paradigm. We illustrate how freezing and thawing models lead to different redshift dependencies of α, and use a combination of current astrophysical observations and local experiments to constrain this class of models, improving the constraints on the key coupling parameter by more than a factor of two, despite considering a more extended parameter space than the one used in previous studies. We also briefly discuss the improvements expected from future facilities and comment on the practical limitations of this class of parametrisations. In particular, we show that sufficiently sensitive data can distinguish between freezing and thawing models, at least if one assumes that the relevant parameter space does not include phantom dark energy models. [ABSTRACT FROM AUTHOR]

Observations reveal a "bulk flow" in the local Universe which is faster and extends to much larger scales than are expected around a typical observer in the standard ΛCDM cosmology. This is expected to result in a scale-dependent dipolar modulation of the acceleration of the expansion rate inferred from observations of objects within the bulk flow. From a maximum-likelihood analysis of the Joint Light-curve Analysis catalogue of Type Ia supernovae, we find that the deceleration parameter, in addition to a small monopole, indeed has a much bigger dipole component aligned with the cosmic microwave background dipole, which falls exponentially with redshift z: q0 = qm + qd.n̂ exp(-z/S). The best fit to data yields qd = −8.03 and S = 0.0262 (⇒d ∼ 100 Mpc), rejecting isotropy (qd = 0) with 3.9σ statistical significance, while qm = −0.157 and consistent with no acceleration (qm = 0) at 1.4σ. Thus the cosmic acceleration deduced from supernovae may be an artefact of our being non-Copernican observers, rather than evidence for a dominant component of "dark energy" in the Universe. [ABSTRACT FROM AUTHOR]

Context. HI intensity mapping (IM) will be used to do precision cosmology, using many existing and upcoming radio observatories. It will measure the integrated HI 21 cm emission signal from "voxels" of the sky at different redshifts. The signal will be contaminated due to absorption, the largest component of which will be the flux absorbed by the HI emitting sources themselves from the potentially bright flux incident on them from background radio continuum sources. Aims. We, for the first time, provide a quantitative estimate of the magnitude of the absorbed flux compared to the emitted HI flux. The ratio of the two fluxes was calculated for various voxels placed at redshifts between 0.1 and 2.5. Methods. We used a cosmological sky simulation of the atomic HI emission line, and summed over the emitted and absorbed fluxes for all sources within voxels at different redshifts. In order to determine the absorbed flux, for each HI source the flux incident from background radio continuum sources was estimated by determining the numbers, sizes, and redshift distribution of radio continuum sources that lie behind it, based on existing observations and simulations. The amount of this incident flux that is absorbed by each HI source was calculated using a relation between integrated optical depth with HI column density determined using observations of damped Lyman-α systems (DLAs) and sub-DLAs. Results. We find that for the same co-moving volume of sky, the HI emission decreases quickly with increasing redshift, while the absorption varies much less with redshift and follows the redshift distribution of faint sources that dominate the number counts of radio continuum sources. This results in the fraction of absorption compared to emission to be negligible in the nearby Universe (up to a redshift of ∼0.5), increases to about 10% at a redshift of one, and continues to increase to about 30% up to a redshift of 2.5. These numbers can vary significantly due to the uncertainty on the exact form of the following relations: firstly, the number counts of radio continuum sources at sub-mJy flux densities; secondly, the relation between integrated optical depth and HI column density of HI sources; and thirdly, the redshift distribution of radio continuum sources up to the highest redshifts. Conclusions. Absorption of the flux incident from background radio continuum sources might become an important contaminant to HI IM signals beyond redshifts of 0.5. The impact of absorption needs to be quantified more accurately using inputs from upcoming deep surveys of radio continuum sources, H I absorption, and HI emission with the Square Kilometre Array and its precursors. [ABSTRACT FROM AUTHOR]

We present an improved determination of the total mass distribution of three massive clusters from the Cluster Lensing and Supernova Survey with Hubble and Hubble Frontier Fields, MACS J1206.2−0847 (z = 0.44), MACS J0416.1−2403 (z = 0.40), Abell S1063 (z = 0.35). We specifically reconstructed the sub-halo mass component with robust stellar kinematics information of cluster galaxies, in combination with precise strong lensing models based on large samples of spectroscopically identified multiple images. We used integral-field spectroscopy in the cluster cores, from the Multi Unit Spectroscopic Explorer on the Very Large Telescope, to measure the stellar velocity dispersion, σ, of 40−60 member galaxies per cluster, covering four to five magnitudes to mF160W ≃ 21.5. We verified the robustness and quantified the accuracy of the velocity dispersion measurements with extensive spectral simulations. With these data, we determined the normalization and slope of the galaxy L–σ Faber–Jackson relation in each cluster and used these parameters as a prior for the scaling relations of the sub-halo population in the mass distribution modeling. When compared to our previous lens models, the inclusion of member galaxies' kinematics provides a similar precision in reproducing the positions of the multiple images. However, the inherent degeneracy between the central effective velocity dispersion, σ0, and truncation radius, rcut, of sub-halos is strongly reduced, thus significantly alleviating possible systematics in the measurements of sub-halo masses. The three independent determinations of the σ0 − rcut scaling relation in each cluster are found to be fully consistent, enabling a statistical determination of sub-halo sizes as a function of σ0, or halo masses. Finally, we derived the galaxy central velocity dispersion functions of the three clusters projected within 16% of their virial radius, finding that they are well in agreement with each other. We argue that such a methodology, when applied to high-quality kinematics and strong lensing data, allows the sub-halo mass functions to be determined and compared with those obtained from cosmological simulations. [ABSTRACT FROM AUTHOR]

The two-point correlation function of the galaxy distribution is a key cosmological observable that allows us to constrain the dynamical and geometrical state of our Universe. To measure the correlation function we need to know both the galaxy positions and the expected galaxy density field. The expected field is commonly specified using a Monte-Carlo sampling of the volume covered by the survey and, to minimize additional sampling errors, this random catalog has to be much larger than the data catalog. Correlation function estimators compare data–data pair counts to data–random and random–random pair counts, where random–random pairs usually dominate the computational cost. Future redshift surveys will deliver spectroscopic catalogs of tens of millions of galaxies. Given the large number of random objects required to guarantee sub-percent accuracy, it is of paramount importance to improve the efficiency of the algorithm without degrading its precision. We show both analytically and numerically that splitting the random catalog into a number of subcatalogs of the same size as the data catalog when calculating random–random pairs and excluding pairs across different subcatalogs provides the optimal error at fixed computational cost. For a random catalog fifty times larger than the data catalog, this reduces the computation time by a factor of more than ten without affecting estimator variance or bias. [ABSTRACT FROM AUTHOR]