Search

The use of Gaussian Process Regression (GPR) for foregrounds mitigation in data collected by the LOw-Frequency ARray (LOFAR) to measure the high-redshift 21-cm signal power spectrum has been shown to have issues of signal loss when the 21-cm signal covariance is misestimated. To address this problem, we have recently introduced covariance kernels obtained by using a Machine Learning based Variational Auto-Encoder (VAE) algorithm in combination with simulations of the 21-cm signal. In this work, we apply this framework to 141 hours ($\approx 10$ nights) of LOFAR data at $z \approx 9.1$, and report revised upper limits of the 21-cm signal power spectrum. Overall, we agree with past results reporting a 2-$\sigma$ upper limit of $\Delta^2_{21} < (80)^2~\rm mK^2$ at $k = 0.075~h~\rm Mpc^{-1}$. Further, the VAE-based kernel has a smaller correlation with the systematic excess noise, and the overall GPR-based approach is shown to be a good model for the data. Assuming an accurate bias correction for the excess noise, we report a 2-$\sigma$ upper limit of $\Delta^2_{21} < (25)^2~\rm mK^2$ at $k = 0.075~h~\rm Mpc^{-1}$. However, we still caution to take the more conservative approach to jointly report the upper limits of the excess noise and the 21-cm signal components., Comment: 5 pages, 3 figures, 2 tables. Accepted for publication in the Monthly Notices of the Royal Astronomical Society (MNRAS) Letters

The high redshift 21-cm signal promises to be a crucial probe of the state of the intergalactic medium (IGM). Understanding the connection between the observed 21-cm power spectrum and the physical quantities intricately associated with the IGM is crucial to fully understand the evolution of our Universe. In this study, we develop an emulator using artificial neural network (ANN) to predict the 21-cm power spectrum from a given set of IGM properties, namely, the bubble size distribution and the volume averaged ionization fraction. This emulator is implemented within a standard Bayesian framework to constrain the IGM parameters from a given 21-cm power spectrum. We compare the performance of the Bayesian method to an alternate method using ANN to predict the IGM parameters from a given input power spectrum, and find that both methods yield similar levels of accuracy, while the ANN is significantly faster. We also use this ANN method of parameter estimation to predict the IGM parameters from a test set contaminated with noise levels expected from the SKA-LOW instrument after 1000 hours of observation. Finally, we train a separate ANN to predict the source parameters from the IGM parameters directly, at a redshift of $z=9.1$, demonstrating the possibility of a non-analytic inference of the source parameters from the IGM parameters for the first time. We achieve high accuracies, with R2-scores ranging between $0.898-0.978$ for the ANN emulator and between $0.966-0.986$ and $0.817-0.981$ for the predictions of IGM parameters from 21-cm power spectrum and source parameters from IGM parameters, respectively. The predictions of the IGM parameters from the Bayesian method incorporating the ANN emulator leads to tight constraints with error bars around $\pm{0.14}$ on the IGM parameters.

A rather clear problem has remained in black hole physics: localizing black holes. One of the recent theoretical ways proposed to identify black hole mergers' hosts is through multi-messenger gravitational lensing: matching the properties of a lensed galactic host with those of a lensed gravitational wave. This paper reviews the most recent literature and introduces some of the ongoing work on the localization of binary black holes and their host galaxies through lensing of gravitational waves and their electromagnetically-bright hosts., Comment: 17 pages, 4 figures. Proceedings from the Royal Society Meeting for Multi-messenger Lensing (2024). Submitted to Phil. Trans. A

The redshifted 21-cm signal from the epoch of reionization (EoR) directly probes the ionization and thermal states of the intergalactic medium during that period. In particular, the distribution of the ionized regions around the radiating sources during EoR introduces scale-dependent features in the spherically-averaged EoR 21-cm signal power spectrum. The goal is to study these scale-dependent features at different stages of reionization using numerical simulations and build a source model-independent framework to probe the properties of the intergalactic medium using EoR 21-cm signal power spectrum measurements. Under the assumption of high spin temperature, we modelled the redshift evolution of the ratio of EoR 21-cm brightness temperature power spectrum and the corresponding density power spectrum using an ansatz consisting of a set of redshift and scale-independent parameters. This set of eight parameters probes the redshift evolution of the average ionization fraction and the quantities related to the morphology of the ionized regions. We have tested this ansatz on different reionization scenarios generated using different simulation algorithms and found that it is able to recover the redshift evolution of the average neutral fraction within an absolute deviation $\lesssim 0.1$. Our framework allows us to interpret 21-cm signal power spectra in terms of parameters related to the state of the IGM. This source model-independent framework can efficiently constrain reionization scenarios using multi-redshift power spectrum measurements with ongoing and future radio telescopes such as LOFAR, MWA, HERA, and SKA. This will add independent information regarding the EoR IGM properties., Comment: 16 pages, 13 figures, 2 tables, Accepted for publication in Astronomy and Astrophysics

Convolutional neural networks (CNNs) are the state-of-the-art technique for identifying strong gravitational lenses. Although they are highly successful in recovering genuine lens systems with a high true-positive rate, the unbalanced nature of the data set (lens systems are rare), still leads to a high false positive rate. For these techniques to be successful in upcoming surveys (e.g. with Euclid) most emphasis should be set on reducing false positives, rather than on reducing false negatives. In this paper, we introduce densely connected neural networks (DenseNets) as the CNN architecture in a new pipeline-ensemble model containing an ensemble of classification CNNs and regression CNNs to classify and rank-order lenses, respectively. We show that DenseNets achieve comparable true positive rates but considerably lower false positive rates (when compared to residual networks; ResNets). Thus, we recommend DenseNets for future missions involving large data sets, such as Euclid, where low false positive rates play a key role in the automated follow-up and analysis of large numbers of strong gravitational lens candidates when human vetting is no longer feasible

We compare the structural properties and dark matter content of star-forming galaxies taken from the CAMELS cosmological simulations to the observed trends derived from the SPARC sample in the stellar mass range $[10^{9}, 10^{11}]\,\textrm{M}_{\odot}$, to provide constraints on the value of cosmological and astrophysical (SN- and AGN-related) parameters. We consider the size-, internal DM fraction-, internal DM mass- and total-stellar mass relations for all the 1065 simulations from the IllustrisTNG, SIMBA and ASTRID suites of CAMELS, and search for the parameters that minimize the $\chi^{2}$ with respect to the observations. For the IllustrisTNG suite, we find the following constraints for the cosmological parameters: $\Omega_{\textrm{m}} = 0.27_{-0.05}^{+0.01}$, $\sigma_{8} = 0.83_{-0.11}^{+0.08}$ and $S_{8} = 0.78_{-0.09}^{+0.03}$, which are consistent within $1\sigma$ with the results from the nine-year WMAP observations. SN feedback-related astrophysical parameters, which describe the departure of outflow wind energy per unit star formation rate and wind velocity from the reference IllustrisTNG simulations, assume the following values: $A_{\textrm{SN1}} = 0.48_{-0.16}^{+0.25}$ and $A_{\textrm{SN2}} = 1.21_{-0.34}^{+0.03}$, respectively. Therefore, simulations with a lower value of outflow wind energy per unit star formation rate with respect to the reference illustrisTNG simulation better reproduce the observations. Simulations based on SIMBA and ASTRID suites predict central dark matter masses substantially larger than those observed in real galaxies, which can be reconciled with observations only by requiring values of $\Omega_{\textrm{m}}$ inconsistent with cosmological constraints for SIMBA, or simulations characterized by unrealistic galaxy mass distributions for ASTRID., Comment: 24 pages, 10 figures, 9 tables. Accepted by MNRAS for publication; Added a reference to sec. 4.1

Inspired by the discussion in the community on possible hidden systematic errors in late universe cosmological probes and non-trivial physical models developed to reduce the Hubble tension, we investigate the Pantheon and Pantheon+ SNe samples for possible deviations from the original $\Lambda$CDM analysis. To simultaneously account for possible systematics or deviations from $\Lambda$CDM, we adopt Gaussian processes to model additional covariance while making no further assumptions on their origin. We explore both stationary and non-stationarity corrections to the covariance. While small changes in the inferred cosmological parameters $H_0$ and $\Omega_{m}$ can occur, we find no statistically significant evidence for missing covariance. We find an upper limit for the Gaussian processes amplitude $\sigma < 0.031$ mag with $95\%$ confidence, which corresponds to $20\%$ of the average statistical error in the Pantheon+ sample. The strongest effect we find on the inferred cosmological parameter posterior can reduce the statistical significance of the Hubble tension between Pantheon+ and Planck estimates from 5.3$\sigma$ to 4.5$\sigma$. Therefore, we conclude that the SN cosmological parameter inference is robust against the analysis modifications studied in this work., Comment: 16 pages, 6 figures, prepared for submission to JCAP

To self-consistently model galactic properties, reionization of the intergalactic medium, and the associated 21-cm signal, we have developed the algorithm polar by integrating the one-dimensional radiative transfer code grizzly with the semi-analytical galaxy formation code L-Galaxies 2020. Our proof-of-concept results are consistent with observations of the star formation rate history, UV luminosity function and the CMB Thomson scattering optical depth. We then investigate how different galaxy formation models affect UV luminosity functions and 21-cm power spectra, and find that while the former are most sensitive to the parameters describing the merger of halos, the latter have a stronger dependence on the supernovae feedback parameters, and both are affected by the escape fraction model., Comment: MNRAS accepted

Wave-mechanical effects in gravitational lensing have long been predicted, and with the discovery of populations of compact transients such as gravitational wave events and fast radio bursts, may soon be observed. We present an observer's review of the relevant theory underlying wave-mechanical effects in gravitational lensing. Starting from the curved-spacetime scalar wave equation, we derive the Fresnel-Kirchoff diffraction integral, and analyze it in the eikonal and wave optics regimes. We answer the question of what makes interference effects observable in some systems but not in others, and how interference effects allow for complementary information to be extracted from lensing systems as compared to traditional measurements. We end by discussing how diffraction effects affect optical depth forecasts and lensing near caustics, and how compact, low-frequency transients like gravitational waves and fast radio bursts provide promising paths to open up the frontier of interferometric gravitational lensing., Comment: Comments and suggestions for additional references welcome

Interstellar scattering (ISS) of radio pulsar emission can be used as a probe of the ionised interstellar medium (IISM) and causes corruptions in pulsar timing experiments. Two types of ISS phenomena (intensity scintillation and pulse broadening) are caused by electron density fluctuations on small scales (< 0.01 AU). Theory predicts that these are related, and both have been widely employed to study the properties of the IISM. Larger scales ($\sim$1-100\,AU) cause measurable changes in dispersion and these can be correlated with ISS observations to estimate the fluctuation spectrum over a very wide scale range. IISM measurements can often be modeled by a homogeneous power-law spatial spectrum of electron density with the Kolmogorov ($-11/3$) spectral exponent. Here we aim to test the validity of using the Kolmogorov exponent with PSR~J0826+2637. We do so using observations of intensity scintillation, pulse broadening and dispersion variations across a wide fractional bandwidth (20 -- 180\,MHz). We present that the frequency dependence of the intensity scintillation in the high frequency band matches the expectations of a Kolmogorov spectral exponent but the pulse broadening in the low frequency band does not change as rapidly as predicted with this assumption. We show that this behavior is due to an inhomogeneity in the scattering region, specifically that the scattering is dominated by a region of transverse size $\sim$40\,AU. The power spectrum of the electron density, however, maintains the Kolmogorov spectral exponent from spatial scales of 5$\times10^{-6}$\,AU to $\sim$100\,AU., Comment: Accepted for publication in MNRAS, typo fixed

The light-cone (LC) anisotropy arises due to cosmic evolution of the cosmic dawn 21-cm signal along the line-of-sight (LoS) axis of the observation volume. The LC effect makes the signal statistically non-ergodic along the LoS axis. The multi-frequency angular power spectrum (MAPS) provides an unbiased alternative to the popular 3D power spectrum as it does not assume statistical ergodicity along every direction in the signal volume. Unlike the 3D power spectrum which mixes the cosmic evolution of the 21-cm signal along the LoS $k$ modes, MAPS keeps the evolution information disentangled. Here we first study the impact of different underlying physical processes during cosmic dawn on the behaviour of the 21-cm MAPS using simulations of various different scenarios and models. We also make error predictions in 21-cm MAPS measurements considering only the system noise and cosmic variance for mock observations of HERA, NenuFAR and SKA-Low. We find that $100~{\rm h}$ of HERA observations will be able to measure 21-cm MAPS at $\geq 3\sigma$ for $\ell \lesssim 1000$ with $0.1\,{\rm MHz}$ channel-width. The better sensitivity of SKA-Low allows reaching this sensitivity up to $\ell \lesssim 3000$. Note that due to the difference in the frequency coverage of the various experiements, the CD-EoR model considered for NenuFAR is different than those used for the HERA and SKA-Low predictions. Considering NenuFAR with the new model, measurements $\geq 2\sigma$ are possible only for $\ell \lesssim 600$ with $0.2\,{\rm MHz}$ channel-width and for a ten times longer observation time of $t_{\rm obs} = 1000~{\rm h}$. However, for the range $300 \lesssim \ell \lesssim 600$ and $t_{\rm obs}=1000~{\rm h}$ more than $3\sigma$ measurements are still possible for NenuFAR when combining consecutive frequency channels within a $5 ~{\rm MHz}$ band., Comment: 21 pages, 12 figures, 3 tables, accepted for publication in MNRAS

The measurement of the $21$ cm signal from the Cosmic Dawn is a major goal for several existing and upcoming radio interferometers such as NenuFAR and the SKA. During this era before the beginning of the Epoch of Reionization, the signal is more difficult to observe due to brighter foregrounds but reveals additional information on the underlying astrophysical processes encoded in the spatial fluctuations of the spin temperature of hydrogen. To interpret future measurements, controlling the level of accuracy of the Lyman-$\alpha$ flux modelling is mandatory. In this work, we evaluate the impact of various approximations that exist in the main fast modelling approach compared to the results of a costly full radiative transfer simulation. The fast SPINTER code, presented in this work, computes the Lyman-$\alpha$ flux including the effect of wing scatterings for an inhomogeneous emissivity field, but assuming an otherwise homogeneous expanding universe. The LICORICE code computes the full radiative transfer in the Lyman-$\alpha$ line without any substantial approximation. We find that the difference between homogeneous and inhomogeneous gas density and temperature is very small for the computed flux. On the contrary, neglecting the effect of gas velocities produces a significant change in the computed flux. We identify the causes (mainly Doppler shifts due to velocity gradients) and quantify the magnitude of the effect in both an idealised setup and a realistic cosmological situation. We find that the amplitude of the effect, up to a factor of $\sim 2$ on the $21$ cm signal power spectrum on some scales (depending on both other model parameters and the redshift), can be easily discriminated with an SKA-like survey and already be approached, particularly for exotic signals, by the ongoing NenuFAR Cosmic Dawn Key Science Program., Comment: 14 pages, accepted for publication in A&A

In a sky-averaged 21-cm signal experiment, the uncertainty on the extracted signal depends mainly on the covariance between the foreground and 21-cm signal models. In this paper, we construct these models using the modes of variation obtained from the Singular Value Decomposition of a set of simulated foreground and 21-cm signals. We present a strategy to reduce this overlap between the 21-cm and foreground modes by simultaneously fitting the spectra from multiple different antennas, which can be used in combination with the method of utilizing the time dependence of foregrounds while fitting multiple drift scan spectra. To demonstrate this idea, we consider two different foreground models (i) a simple foreground model, where we assume a constant spectral index over the sky, and (ii) a more realistic foreground model, with a spatial variation of the spectral index. For the simple foreground model, with just a single antenna design, we are able to extract the signal with good accuracy if we simultaneously fit the data from multiple time slices. The 21-cm signal extraction is further improved when we simultaneously fit the data from different antennas as well. This improvement becomes more pronounced while using the more realistic mock observations generated from the detailed foreground model. We find that even if we fit multiple time slices, the recovered signal is biased and inaccurate for a single antenna. However, simultaneously fitting the data from different antennas reduces the bias and the uncertainty by a factor of 2-3 on the extracted 21-cm signal., Comment: 12 pages, 13 figures. Accepted for publication in MNRAS. Accompanying code is available https://github.com/anchal-009/SAVED21cm

The unprecedented number of gravitational lenses expected from new-generation facilities such as the ESA Euclid telescope and the Vera Rubin Observatory makes it crucial to rethink our classical approach to lens-modelling. In this paper, we present LeMoN (Lens Modelling with Neural networks): a new machine-learning algorithm able to analyse hundreds of thousands of gravitational lenses in a reasonable amount of time. The algorithm is based on a Bayesian Neural Network: a new generation of neural networks able to associate a reliable confidence interval to each predicted parameter. We train the algorithm to predict the three main parameters of the Singular Isothermal Ellipsoid model (the Einstein radius and the two components of the ellipticity) by employing two simulated datasets built to resemble the imaging capabilities of the Hubble Space Telescope and the forthcoming Euclid satellite. In this work, we assess the accuracy of the algorithm and the reliability of the estimated uncertainties by applying the network to several simulated datasets of 10.000 images each. We obtain accuracies comparable to previous studies present in the current literature and an average modelling time of just 0.5s per lens. Finally, we apply the LeMoN algorithm to a pilot dataset of real lenses observed with HST during the SLACS program, obtaining unbiased estimates of their SIE parameters. The code is publicly available on GitHub (https://github.com/fab-gentile/LeMoN)., Comment: Submitted to MNRAS; The code will be available on GitHub upon acceptance of the paper

To unlock the full spectrum of astrophysical and cosmological applications of gravitational-wave detections, it is essential to localise the associated black-hole mergers to high precision inside their host galaxies. One possible method to achieve this is to compare the properties of multiple detections of gravitationally-lensed binary black-hole merger events with the properties of strong gravitational lens systems located in the joint sky localisation of the gravitational-wave detections. In this work, we simulate the population of binary black-hole mergers lensed by galaxy-scale lenses and detectable by LIGO-Virgo-Kagra in the coming decade and the population of galaxy-scale strong lenses that will be detected by Euclid. We use these simulations to investigate the prospects for localising strongly lensed binary black-hole mergers inside the lensed galaxies of 'Euclid-like' galaxy-scale strong lenses. We find that for 20-50% of strongly lensed gravitational wave events the lens system is detectable with Euclid, if the event falls in its survey footprint. Of these, we expect to correctly identify the strongly-lensed host galaxy as likely (with posterior probability) host galaxy - based on Bayesian evidence ranking of candidate hosts - for 34.6-21.9% of quadruply-lensed gravitational-wave events when given an a-priori 1-5 deg^2 gravitational-wave-only sky localisation. For triply and doubly lensed gravitational-wave events, this becomes 29.8-14.9% and 16.4-6.6% respectively. If successfully identified, however, the localisation can be better than a fraction of the host-galaxy size, i.e. of order milli-arcseconds. A first detection in the coming decade, however, probably requires dedicated deep and high-resolution follow-ups and continued upgrades in the current and planned gravitational-wave detectors., Comment: 23 pages, 14 figures, 3 tables. Accepted for publication in MNRAS

In this paper we will list a few important goals that need to be addressed in the next decade, also taking into account the current discordances between the different cosmological probes, such as the disagreement in the value of the Hubble constant $H_0$, the $\sigma_8$--$S_8$ tension, and other less statistically significant anomalies. While these discordances can still be in part the result of systematic errors, their persistence after several years of accurate analysis strongly hints at cracks in the standard cosmological scenario and the necessity for new physics or generalisations beyond the standard model. In this paper, we focus on the $5.0\,\sigma$ tension between the {\it Planck} CMB estimate of the Hubble constant $H_0$ and the SH0ES collaboration measurements. After showing the $H_0$ evaluations made from different teams using different methods and geometric calibrations, we list a few interesting new physics models that could alleviate this tension and discuss how the next decade's experiments will be crucial. Moreover, we focus on the tension of the {\it Planck} CMB data with weak lensing measurements and redshift surveys, about the value of the matter energy density $\Omega_m$, and the amplitude or rate of the growth of structure ($\sigma_8,f\sigma_8$). We list a few interesting models proposed for alleviating this tension, and we discuss the importance of trying to fit a full array of data with a single model and not just one parameter at a time. Additionally, we present a wide range of other less discussed anomalies at a statistical significance level lower than the $H_0$--$S_8$ tensions which may also constitute hints towards new physics, and we discuss possible generic theoretical approaches that can collectively explain the non-standard nature of these signals.[Abridged], Comment: Contribution to Snowmass 2021. 224 pages, 27 figures. Accepted for publication in JHEAp

Forthcoming large imaging surveys such as Euclid and the Vera Rubin Observatory Legacy Survey of Space and Time are expected to find more than $10^5$ strong gravitational lens systems, including many rare and exotic populations such as compound lenses, but these $10^5$ systems will be interspersed among much larger catalogues of $\sim10^9$ galaxies. This volume of data is too much for visual inspection by volunteers alone to be feasible and gravitational lenses will only appear in a small fraction of these data which could cause a large amount of false positives. Machine learning is the obvious alternative but the algorithms' internal workings are not obviously interpretable, so their selection functions are opaque and it is not clear whether they would select against important rare populations. We design, build, and train several Convolutional Neural Networks (CNNs) to identify strong gravitational lenses using VIS, Y, J, and H bands of simulated data, with F1 scores between 0.83 and 0.91 on 100,000 test set images. We demonstrate for the first time that such CNNs do not select against compound lenses, obtaining recall scores as high as 76\% for compound arcs and 52\% for double rings. We verify this performance using Hubble Space Telescope (HST) and Hyper Suprime-Cam (HSC) data of all known compound lens systems. Finally, we explore for the first time the interpretability of these CNNs using Deep Dream, Guided Grad-CAM, and by exploring the kernels of the convolutional layers, to illuminate why CNNs succeed in compound lens selection., Comment: MNRAS accepted. 16 pages, 17 figures, The code used in this paper is publicly available at github.com/JoshWilde/LensFindery-McLensFinderFace

Reconstructing lens potentials and lensed sources can easily become an underconstrained problem, even when the degrees of freedom are low, due to degeneracies, particularly when potential perturbations superimposed on a smooth lens are included. Regularization has traditionally been used to constrain the solutions where the data failed to do so, e.g. in unlensed parts of the source. In this exploratory work, we go beyond the usual choices of regularization and adopt observationally motivated priors for the source brightness. We also perform a similar comparison when reconstructing lens potential perturbations, which are assumed to be stationary, i.e. permeate the entire field of view. We find that physically motivated priors lead to lower residuals, avoid overfitting, and are decisively preferred within a Bayesian quantitative framework in all the examples considered. For the perturbations, choosing the wrong regularization can have a detrimental effect that even high-quality data cannot correct for, while using a purely smooth lens model can absorb them to a very high degree and lead to biased solutions. Finally, our new implementation of the semi-linear inversion technique provides the first quantitative framework for measuring degeneracies between the source and the potential perturbations., Comment: 25 pages, 21 figures, 5 tables, submitted to MNRAS

We investigate tomography of 21-cm brightness temperature fluctuations during the Dark Ages as a probe for constraining primordial non-Gaussianity. We expand the 21- cm brightness temperature up to cubic order in perturbation theory and improve previous models of the signal by including the effect of the free electron fraction. Using modified standard perturbation theory methods that include baryonic pressure effects we derive an improved secondary bispectrum and for the first time derive the secondary trispectrum of 21-cm brightness temperature fluctuations. We then forecast the amount of information available from the Dark Ages to constrain primordial non-Gaussianity, including the imprints of massive particle exchange during inflation and we determine how much signal is lost due to secondary non-Gaussianity. We find that although secondary non-Gaussianity swamps the primordial signal, primordial non-Gaussianity can still be extracted with signal-to-noise ratios that surpass current and future CMB experiments by several orders of magnitude, depending on the experimental setup. Furthermore, we conclude that for the bi- and trispectra of massive particle exchange marginalizing over other primordial shapes affects signal-to-noise ratios more severely than secondary shapes. Baryonic pressure effects turn out to have a negligible impact on our forecasts, even at scales close to the Jeans scale. The results of this work reinforce the prospects of 21-cm brightness temperature fluctuations from the Dark Ages as the ultimate probe for primordial non-Gaussianity., Comment: 47 pages, 20 figures, code available at https://github.com/tsfloss/21-cm-Dark-Ages-Fisher-Forecast. v2: minor typos corrected. v3: matches published version

NenuFAR, the New Extension in Nancay Upgrading LOFAR, is currently in its early science phase. It is in this context that the Cosmic Filaments and Magnetism Pilot Survey is observing sources with the array as it is still under construction - with 57 (56 core, 1 distant) out of a total planned 102 (96 core, 6 distant) mini-arrays online at the time of observation - to get a first look at the low-frequency sky with NenuFAR. One of its targets is the Coma galaxy cluster: a well-known object, host of the prototype radio halo. It also hosts other features of scientific import, including a radio relic, along with a bridge of emission connecting it with the halo. It is thus a well-studied object. In this paper, we show the first confirmed NenuFAR detection of the radio halo and radio relic of the Coma cluster at 34.4 MHz, with associated intrinsic flux density estimates: we find an integrated flux value of 106.3 +- 3.5 Jy for the radio halo, and 102.0 +- 7.4 Jy for the radio relic. These are upper bound values, as they do not include point-source subtraction. We also give an explanation of the technical difficulties encountered in reducing the data, along with steps taken to resolve them. This will be helpful for other scientific projects which will aim to make use of standalone NenuFAR imaging observations in the future., Comment: 14 pages, 4 figures, special issue of the Galaxies journal

The central dark-matter fraction of galaxies is sensitive to feedback processes during galaxy formation. Strong gravitational lensing has been effective in the precise measurement of the dark-matter fraction inside massive early-type galaxies. Here, we compare the projected dark-matter fraction of early-type galaxies inferred from the SLACS strong-lens survey, with those obtained from the EAGLE, Illustris, and IllustrisTNG hydro-dynamical simulations. Previous comparisons with some simulations revealed a large discrepancy, with considerably higher inferred dark-matter fractions -- by factors 2-3 -- inside half of the effective radius in observed strong-lens galaxies as compared to simulated galaxies. Here, we report good agreement between EAGLE and SLACS for the dark-matter fractions inside both half of the effective radius and the effective radius as a function of the galaxy's stellar mass, effective radius, and total mass-density slope. However, for IllustrisTNG and Illustris, the dark-matter fractions are lower than observed. This work consistently assumes a Chabrier IMF, which suggests that a different IMF (although not excluded) is not necessary to resolve this mismatch. The differences in the stellar feedback model between EAGLE and Illustris and IllustrisTNG, are likely the dominant cause of the difference in their dark-matter fraction, and density slope., Comment: Matches with the published version. Fixed an error in Fig 2. Updated reference. Conclusions unchanged. 7 pages, 4 figures

Strongly lensed quasars can provide measurements of the Hubble constant ($H_{0}$) independent of any other methods. One of the key ingredients is exquisite high-resolution imaging data, such as Hubble Space Telescope (HST) imaging and adaptive-optics (AO) imaging from ground-based telescopes, which provide strong constraints on the mass distribution of the lensing galaxy. In this work, we expand on the previous analysis of three time-delay lenses with AO imaging (RXJ1131-1231, HE0435-1223, and PG1115+080), and perform a joint analysis of J0924+0219 by using AO imaging from the Keck Telescope, obtained as part of the SHARP (Strong lensing at High Angular Resolution Program) AO effort, with HST imaging to constrain the mass distribution of the lensing galaxy. Under the assumption of a flat $\Lambda$CDM model with fixed $\Omega_{\rm m}=0.3$, we show that by marginalizing over two different kinds of mass models (power-law and composite models) and their transformed mass profiles via a mass-sheet transformation, we obtain $\Delta t_{\rm BA}h\hat{\sigma}_{v}^{-2}=6.89\substack{+0.8\\-0.7}$ days, $\Delta t_{\rm CA}h\hat{\sigma}_{v}^{-2}=10.7\substack{+1.6\\-1.2}$ days, and $\Delta t_{\rm DA}h\hat{\sigma}_{v}^{-2}=7.70\substack{+1.0\\-0.9}$ days, where $h=H_{0}/100~\rm km\,s^{-1}\,Mpc^{-1}$ is the dimensionless Hubble constant and $\hat{\sigma}_{v}=\sigma^{\rm ob}_{v}/(280~\rm km\,s^{-1})$ is the scaled dimensionless velocity dispersion. Future measurements of time delays with 10% uncertainty and velocity dispersion with 5% uncertainty would yield a $H_0$ constraint of $\sim15$% precision., Comment: 11 pages, 7 figures

When gravitational waves pass near massive astrophysical objects, they can be gravitationally lensed. The lensing can split them into multiple wave-fronts, magnify them, or imprint beating patterns on the waves. Here we focus on the multiple images produced by strong lensing. In particular, we investigate strong lensing forecasts, the rate of lensing, and the role of lensing statistics in strong lensing searches. Overall, we find a reasonable rate of lensed detections for double, triple, and quadruple images at the LIGO--Virgo--KAGRA design sensitivity. We also report the rates for A+ and LIGO Voyager and briefly comment on potential improvements due to the inclusion of sub-threshold triggers. We find that most galaxy-lensed events originate from redshifts $z \sim 1-4$ and report the expected distribution of lensing parameters for the observed events. Besides forecasts, we investigate the role of lensing forecasts in strong lensing searches, which explore repeated event pairs. One problem associated with the searches is the rising number of event pairs, which leads to a rapidly increasing false alarm probability. We show how knowledge of the expected galaxy lensing time delays in our searches allow us to tackle this problem. Once the time delays are included, the false alarm probability increases linearly (similar to non-lensed searches) instead of quadratically with time, significantly improving the search. For galaxy cluster lenses, the improvement is less significant. The main uncertainty associated with these forecasts are the merger-rate density estimates at high redshift, which may be better resolved in the future.

We present a sample of 16 likely strong gravitational lenses identified in the VST Optical Imaging of the CDFS and ES1 fields (VOICE survey) using Convolutional Neural Networks (CNNs). We train two different CNNs on composite images produced by superimposing simulated gravitational arcs on real Luminous Red Galaxies observed in VOICE. Specifically, the first CNN is trained on single-band images and more easily identifies systems with large Einstein radii, while the second one, trained on composite RGB images, is more accurate in retrieving systems with smaller Einstein radii. We apply both networks to real data from the VOICE survey, taking advantage of the high limiting magnitude (26.1 in the r-band) and low PSF FWHM (0.8" in the r-band) of this deep survey. We analyse $\sim21,200$ images with $mag_r<21.5$, identifying 257 lens candidates. To retrieve a high-confidence sample and to assess the accuracy of our technique, nine of the authors perform a visual inspection. Roughly 75% of the systems are classified as likely lenses by at least one of the authors. Finally, we assemble the LIVE sample (Lenses In VoicE) composed by the 16 systems passing the chosen grading threshold. Three of these candidates show likely lensing features when observed by the Hubble Space Telescope. This work represents a further confirmation of the ability of CNNs to inspect large samples of galaxies searching for gravitational lenses. These algorithms will be crucial to exploit the full scientific potential of forthcoming surveys with the Euclid satellite and the Vera Rubin Observatory, Comment: 15 pages, 12 figures, submitted to MNRAS

High-redshift 21-cm observations will provide crucial insights into the physical processes of the Epoch of Reionization. Next-generation interferometers such as the Square Kilometer Array will have enough sensitivity to directly image the 21-cm fluctuations and trace the evolution of the ionizing fronts. In this work, we develop an inferential approach to recover the sources and IGM properties of the process of reionization using the number and, in particular, the morphological pattern spectra of the ionized regions extracted from realistic mock observations. To do so, we extend the Markov Chain Monte Carlo analysis tool 21CMMC by including these 21-cm tomographic statistics and compare this method to only using the power spectrum. We demonstrate that the evolution of the number-count and morphology of the ionized regions as a function of redshift provides independent information to disentangle multiple reionization scenarios because it probes the average ionizing budget per baryon. Although less precise, we find that constraints inferred using 21-cm tomographic statistics are more robust to the presence of contaminants such as foreground residuals. This work highlights that combining power spectrum and tomographic analyses more accurately recovers the astrophysics of reionization., Comment: 28 pages, 16 figures. Accepted to MNRAS

We derive joint constraints on the warm dark matter (WDM) half-mode scale by combining the analyses of a selection of astrophysical probes: strong gravitational lensing with extended sources, the Lyman-$\alpha$ forest, and the number of luminous satellites in the Milky Way. We derive an upper limit of $\lambda_{\rm hm}=0.089{\rm~Mpc~h^{-1} }$ at the 95 per cent confidence level, which we show to be stable for a broad range of prior choices. Assuming a Planck cosmology and that WDM particles are thermal relics, this corresponds to an upper limit on the half-mode mass of $M_{\rm hm }< 3 \times 10^{7} {\rm~M_{\odot}~h^{-1}}$, and a lower limit on the particle mass of $m_{\rm th }> 6.048 {\rm~keV}$, both at the 95 per cent confidence level. We find that models with $\lambda_{\rm hm}> 0.223 {\rm~Mpc~h^{-1} }$ (corresponding to $m_{\rm th }> 2.552 {\rm~keV}$ and $M_{\rm hm }< 4.8 \times 10^{8} {\rm~M_{\odot}~h^{-1}}$) are ruled out with respect to the maximum likelihood model by a factor $\leq 1/20$. For lepton asymmetries $L_6>10$, we rule out the $7.1 {\rm~keV}$ sterile neutrino dark matter model, which presents a possible explanation to the unidentified $3.55 {\rm~keV}$ line in the Milky Way and clusters of galaxies. The inferred 95 percentiles suggest that we further rule out the ETHOS-4 model of self-interacting DM. Our results highlight the importance of extending the current constraints to lower half-mode scales. We address important sources of systematic errors and provide prospects for how the constraints of these probes can be improved upon in the future., Comment: 16 pages, 3 Figures, minor updates of the results

Using the latest upper limits on the 21-cm power spectrum at $z\approx9.1$ from the Low Frequency Array (LOFAR), we explore regions of parameter space which are inconsistent with the data. We use 21CMMC, a Monte Carlo Markov Chain sampler of 21cmFAST which directly forward models the 3D cosmic 21-cm signal in a fully Bayesian framework. We use the astrophysical parameterisation from 21cmFAST, which includes mass-dependent star formation rates and ionising escape fractions as well as soft-band X-ray luminosities to place limits on the properties of the high-$z$ galaxies. Further, we connect the disfavoured regions of parameter space with existing observational constraints on the Epoch of Reionisation such as ultra-violet (UV) luminosity functions, background UV photoionisation rate, intergalactic medium (IGM) neutral fraction and the electron scattering optical depth. We find that all models exceeding the 21-cm signal limits set by LOFAR at $z\approx9.1$ are excluded at $\gtrsim2\sigma$ by other probes. Finally, we place limits on the IGM spin temperature from LOFAR, disfavouring at 95 per cent confidence spin temperatures below $\sim2.6$ K across an IGM neutral fraction range of $0.15 \lesssim \bar{x}_{H{\scriptscriptstyle I}} \lesssim 0.6$. Note, these limits are only obtained from 141 hrs of data in a single redshift bin. With tighter upper limits, across multiple redshift bins expected in the near future from LOFAR, more viable models will be ruled out. Our approach demonstrates the potential of forward modelling tools such as 21CMMC in combining 21-cm observations with other high-$z$ probes to constrain the astrophysics of galaxies., Comment: 12 pages, 4 pages, 1 table. Submitted to MNRAS, comments welcome

The key challenge in the observation of the redshifted 21-cm signal from cosmic reionization is its separation from the much brighter foreground emission. Such separation relies on the different spectral properties of the two components, although, in real life, the foreground intrinsic spectrum is often corrupted by the instrumental response, inducing systematic effects that can further jeopardize the measurement of the 21-cm signal. In this paper, we use Gaussian Process Regression to model both foreground emission and instrumental systematics in $\sim 2$ hours of data from the Hydrogen Epoch of Reionization Array. We find that a simple co-variance model with three components matches the data well, giving a residual power spectrum with white noise properties. These consist of an "intrinsic" and instrumentally corrupted component with a coherence-scale of 20 MHz and 2.4 MHz respectively (dominating the line of sight power spectrum over scales $k_{\parallel} \le 0.2$ h cMpc$^{-1}$) and a baseline dependent periodic signal with a period of $\sim 1$ MHz (dominating over $k_{\parallel} \sim 0.4 - 0.8$h cMpc$^{-1}$) which should be distinguishable from the 21-cm EoR signal whose typical coherence-scales is $\sim 0.8$ MHz., Comment: 15 pages, 15 figures, 1 table, Accepted to MNRAS

This paper presents the results from one of the first observations of ionospheric scintillation taken using the Low-Frequency Array (LOFAR). The observation was of the strong natural radio source Cas A, taken overnight on 18-19 August 2013, and exhibited moderately strong scattering effects in dynamic spectra of intensity received across an observing bandwidth of 10-80MHz. Delay-Doppler spectra (the 2-D FFT of the dynamic spectrum) from the first hour of observation showed two discrete parabolic arcs, one with a steep curvature and the other shallow, which can be used to provide estimates of the distance to, and velocity of, the scattering plasma. A cross-correlation analysis of data received by the dense array of stations in the LOFAR "core" reveals two different velocities in the scintillation pattern: a primary velocity of ~30m/s with a north-west to south-east direction, associated with the steep parabolic arc and a scattering altitude in the F-region or higher, and a secondary velocity of ~110m/s with a north-east to south-west direction, associated with the shallow arc and a scattering altitude in the D-region. Geomagnetic activity was low in the mid-latitudes at the time, but a weak sub-storm at high latitudes reached its peak at the start of the observation. An analysis of Global Navigation Satellite Systems (GNSS) and ionosonde data from the time reveals a larger-scale travelling ionospheric disturbance (TID), possibly the result of the high-latitude activity, travelling in the north-west to south-east direction, and, simultaneously, a smaller--scale TID travelling in a north-east to south-west direction, which could be associated with atmospheric gravity wave activity. The LOFAR observation shows scattering from both TIDs, at different altitudes and propagating in different directions. To the best of our knowledge this is the first time that such a phenomenon has been reported., Comment: 24 pages, 16 figures. Accepted for open-access publication in the Journal of Space Weather and Space Climate. For associated movie file, see https://www.swsc-journal.org/10.1051/swsc/2020010/olm

Ultra-compact massive galaxies UCMGs, i.e. galaxies with stellar masses $M_{*} > 8 \times 10^{10} M_{\odot}$ and effective radii $R_{e} < 1.5$ kpc, are very rare systems, in particular at low and intermediate redshifts. Their origin as well as their number density across cosmic time are still under scrutiny, especially because of the paucity of spectroscopically confirmed samples. We have started a systematic census of UCMG candidates within the ESO Kilo Degree Survey, together with a large spectroscopic follow-up campaign to build the largest possible sample of confirmed UCMGs. This is the third paper of the series and the second based on the spectroscopic follow-up program. Here, we present photometrical and structural parameters of 33 new candidates at redshifts $0.15 \lesssim z \lesssim 0.5$ and confirm 19 of them as UCMGs, based on their nominal spectroscopically inferred $M_{*}$ and $R_{e}$. This corresponds to a success rate of $\sim 58\%$, nicely consistent with our previous findings. The addition of these 19 newly confirmed objects, allows us to fully assess the systematics on the system selection, and finally reduce the number density uncertainties. Moreover, putting together the results from our current and past observational campaigns and some literature data, we build the largest sample of UCMGs ever collected, comprising 92 spectroscopically confirmed objects at $0.1 \lesssim z \lesssim 0.5$. This number raises to 116, allowing for a $3\sigma$ tolerance on the $M_{*}$ and $R_{e}$ thresholds for the UCMG definition. For all these galaxies we have estimated the velocity dispersion values at the effective radii which have been used to derive a preliminary mass-velocity dispersion correlation.

Persistent tension between low-redshift observations and the Cosmic Microwave Background radiation (CMB), in terms of two fundamental distance scales set by the sound horizon $r_d$ and the Hubble constant $H_0$, suggests new physics beyond the Standard Model or residual systematics. We examine recently updated distance calibrations from Cepheids, gravitational lensing time-delay observations, and the Tip of the Red Giant Branch. Calibrating the Baryon Acoustic Oscillations (BAO) and Type Ia supernovae with combinations of the distance indicators, we obtain a joint and self-consistent measurement of $H_0$ and $r_d$ at low redshift, independent of cosmological models and CMB inference. In an attempt to alleviate the tension between late-time and CMB-based measurements, we consider four extensions of the standard $\Lambda$CDM model. The sound horizon from our different measurements is $r_d=(137\pm3^{stat.}\pm2^{syst.})$~Mpc. Depending on the adopted distance indicators, the $combined$ tension in $H_0$ and $r_d$ ranges between 2.3 and 5.1 $\sigma$. We find that modifications of $\Lambda$CDM that change the physics after recombination fail to solve the problem, for the reason that they only resolve the tension in $H_0$, while the tension in $r_d$ remains unchanged. Pre-recombination extensions (with early dark energy or the effective number of neutrinos $\rm{N}_{\rm{eff}}=3.24 \pm 0.16$) are allowed by the data, unless the calibration from Cepheids is included. Results from time-delay lenses are consistent with those from distance-ladder calibrations and point to a discrepancy between absolute distance scales measured from the CMB (assuming the standard cosmological model) and late-time observations. New proposals to resolve this tension should be examined with respect to reconciling not only the Hubble constant but also the sound horizon derived from the CMB and other cosmological probes., Comment: Accepted for publication in A&A; 14 pages, 8 figures, 4 tables

The local expansion rate of the Universe is parametrized by the Hubble constant, $H_0$, the ratio between recession velocity and distance. Different techniques lead to inconsistent estimates of $H_0$. Observations of Type Ia supernovae (SNe) can be used to measure $H_0$, but this requires an external calibrator to convert relative distances to absolute ones. We use the angular diameter distance to strong gravitational lenses as a suitable calibrator, which is only weakly sensitive to cosmological assumptions. We determine the angular diameter distances to two gravitational lenses, $810^{+160}_{-130}$ and $1230^{+180}_{-150}$~Mpc, at redshifts of $z=0.295$ and $0.6304$. Using these absolute distances to calibrate 740 previously-measured relative distances to SNe, we measure the Hubble constant to be $H_0=82.4^{+8.4}_{-8.3} ~{\rm km\,s^{-1}\,Mpc^{-1}}$., Comment: This paper presents the measurements of angular diameter distances to two time-delay lenses, and the Hubble constant derived only from these two distances and the JLA supernova sample. One of the distance measurements is further used for the cosmological inference in the H0LiCOW XIII paper (arxiv:1907.04869). Published in Science

We present a measurement of the Hubble constant ($H_{0}$) and other cosmological parameters from a joint analysis of six gravitationally lensed quasars with measured time delays. All lenses except the first are analyzed blindly with respect to the cosmological parameters. In a flat $\Lambda$CDM cosmology, we find $H_{0} = 73.3_{-1.8}^{+1.7}$, a 2.4% precision measurement, in agreement with local measurements of $H_{0}$ from type Ia supernovae calibrated by the distance ladder, but in $3.1\sigma$ tension with $Planck$ observations of the cosmic microwave background (CMB). This method is completely independent of both the supernovae and CMB analyses. A combination of time-delay cosmography and the distance ladder results is in $5.3\sigma$ tension with $Planck$ CMB determinations of $H_{0}$ in flat $\Lambda$CDM. We compute Bayes factors to verify that all lenses give statistically consistent results, showing that we are not underestimating our uncertainties and are able to control our systematics. We explore extensions to flat $\Lambda$CDM using constraints from time-delay cosmography alone, as well as combinations with other cosmological probes, including CMB observations from $Planck$, baryon acoustic oscillations, and type Ia supernovae. Time-delay cosmography improves the precision of the other probes, demonstrating the strong complementarity. Allowing for spatial curvature does not resolve the tension with $Planck$. Using the distance constraints from time-delay cosmography to anchor the type Ia supernova distance scale, we reduce the sensitivity of our $H_0$ inference to cosmological model assumptions. For six different cosmological models, our combined inference on $H_{0}$ ranges from $\sim73$-$78~\mathrm{km~s^{-1}~Mpc^{-1}}$, which is consistent with the local distance ladder constraints., Comment: Accepted for publication in MNRAS; 23 pages, 13 figures, 8 tables

We present the measurement of the Hubble Constant, $H_0$, with three strong gravitational lens systems. We describe a blind analysis of both PG1115+080 and HE0435-1223 as well as an extension of our previous analysis of RXJ1131-1231. For each lens, we combine new adaptive optics (AO) imaging from the Keck Telescope, obtained as part of the SHARP AO effort, with Hubble Space Telescope (HST) imaging, velocity dispersion measurements, and a description of the line-of-sight mass distribution to build an accurate and precise lens mass model. This mass model is then combined with the COSMOGRAIL measured time delays in these systems to determine $H_{0}$. We do both an AO-only and an AO+HST analysis of the systems and find that AO and HST results are consistent. After unblinding, the AO-only analysis gives $H_{0}=82.8^{+9.4}_{-8.3}~\rm km\,s^{-1}\,Mpc^{-1}$ for PG1115+080, $H_{0}=70.1^{+5.3}_{-4.5}~\rm km\,s^{-1}\,Mpc^{-1}$ for HE0435-1223, and $H_{0}=77.0^{+4.0}_{-4.6}~\rm km\,s^{-1}\,Mpc^{-1}$ for RXJ1131-1231. The joint AO-only result for the three lenses is $H_{0}=75.6^{+3.2}_{-3.3}~\rm km\,s^{-1}\,Mpc^{-1}$. The joint result of the AO+HST analysis for the three lenses is $H_{0}=76.8^{+2.6}_{-2.6}~\rm km\,s^{-1}\,Mpc^{-1}$. All of the above results assume a flat $\Lambda$ cold dark matter cosmology with a uniform prior on $\Omega_{\textrm{m}}$ in [0.05, 0.5] and $H_{0}$ in [0, 150] $\rm km\,s^{-1}\,Mpc^{-1}$. This work is a collaboration of the SHARP and H0LiCOW teams, and shows that AO data can be used as the high-resolution imaging component in lens-based measurements of $H_0$. The full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper., Comment: 33 pages, 23 figures

Interferometry of the cosmic 21-cm signal is set to revolutionise our understanding of the Epoch of Reionisation (EoR) and the Cosmic Dawn (CD). The culmination of ongoing efforts will be the upcoming Square Kilometre Array (SKA), which will provide tomography of the 21-cm signal from the first billion years of our Universe. Using a galaxy formation model informed by high-$z$ luminosity functions, here we forecast the accuracy with which the first phase of SKA-low (SKA1-low) can constrain the properties of the unseen galaxies driving the astrophysics of the EoR and CD. We consider three observing strategies: (i) deep (1000h on a single field); (ii) medium-deep (100hr on 10 independent fields); and (iii) shallow (10hr on 100 independent fields). Using the 21-cm power spectrum as a summary statistic, and conservatively only using the 21-cm signal above the foreground wedge, we predict that all three observing strategies should recover astrophysical parameters to a fractional precision of $\sim 0.1$ -- 10 per cent. The reionisation history is recovered to an uncertainty of $\Delta z \lesssim 0.1$ (1$\sigma$) for the bulk of its duration. The medium-deep strategy, balancing thermal noise against cosmic variance, results in the tightest constraints, slightly outperforming the deep strategy. The shallow observational strategy performs the worst, with up to a $\sim 10$ -- 60 per cent increase in the recovered uncertainty. We note, however, that non-Gaussian summary statistics, tomography, as well as unbiased foreground removal would likely favour the deep strategy., Comment: 10 pages, 2 figures, submitted to MNRAS, comments welcome

We present the lens mass model of the quadruply-imaged gravitationally lensed quasar WFI2033-4723, and perform a blind cosmographical analysis based on this system. Our analysis combines (1) time-delay measurements from 14 years of data obtained by the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL) collaboration, (2) high-resolution $\textit{Hubble Space Telescope}$ imaging, (3) a measurement of the velocity dispersion of the lens galaxy based on ESO-MUSE data, and (4) multi-band, wide-field imaging and spectroscopy characterizing the lens environment. We account for all known sources of systematics, including the influence of nearby perturbers and complex line-of-sight structure, as well as the parametrization of the light and mass profiles of the lensing galaxy. After unblinding, we determine the effective time-delay distance to be $4784_{-248}^{+399}~\mathrm{Mpc}$, an average precision of $6.6\%$. This translates to a Hubble constant $H_{0} = 71.6_{-4.9}^{+3.8}~\mathrm{km~s^{-1}~Mpc^{-1}}$, assuming a flat $\Lambda$CDM cosmology with a uniform prior on $\Omega_\mathrm{m}$ in the range [0.05, 0.5]. This work is part of the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration, and the full time-delay cosmography results from a total of six strongly lensed systems are presented in a companion paper (H0LiCOW XIII)., Comment: Version accepted by MNRAS. 29 pages including appendix, 17 figures, 6 tables. arXiv admin note: text overlap with arXiv:1607.01403

We present an analysis of seven strongly gravitationally lensed quasars and the corresponding constraints on the properties of dark matter. Our results are derived by modelling the lensed image positions and flux-ratios using a combination of smooth macro models and a population of low-mass haloes within the mass range 10^6 to 10^9 Msun. Our lens models explicitly include higher-order complexity in the form of stellar discs and luminous satellites, as well as low-mass haloes located along the observed lines of sight for the first time. Assuming a Cold Dark Matter (CDM) cosmology, we infer an average total mass fraction in substructure of f_sub = 0.012^{+0.007}_{-0.004} (68 per cent confidence limits), which is in agreement with the predictions from CDM hydrodynamical simulations to within 1 sigma. This result is closer to the predictions than those from previous studies that did not include line-of-sight haloes. Under the assumption of a thermal relic dark matter model, we derive a lower limit on the particle relic mass of m th > 5.58 keV (95 per cent confidence limits), which is consistent with a value of m_th > 5.3 keV from the recent analysis of the Ly-alpha forest. We also identify two main sources of possible systematic errors and conclude that deeper investigations in the complex structure of lens galaxies as well as the size of the background sources should be a priority for this field., Comment: 14 pages, 7 figures, accepted for publication in MNRAS

We use nine different galaxy formation scenarios in ten cosmological simulation boxes from the EAGLE suite of {\Lambda}CDM hydrodynamical simulations to assess the impact of feedback mechanisms in galaxy formation and compare these to observed strong gravitational lenses. To compare observations with simulations, we create strong lenses with $M_\star$ > $10^{11}$ $M_\odot$ with the appropriate resolution and noise level, and model them with an elliptical power-law mass model to constrain their total mass density slope. We also obtain the mass-size relation of the simulated lens-galaxy sample. We find significant variation in the total mass density slope at the Einstein radius and in the projected stellar mass-size relation, mainly due to different implementations of stellar and AGN feedback. We find that for lens selected galaxies, models with either too weak or too strong stellar and/or AGN feedback fail to explain the distribution of observed mass-density slopes, with the counter-intuitive trend that increasing the feedback steepens the mass density slope around the Einstein radius ($\approx$ 3-10 kpc). Models in which stellar feedback becomes inefficient at high gas densities, or weaker AGN feedback with a higher duty cycle, produce strong lenses with total mass density slopes close to isothermal (i.e. -d log({\rho})/d log(r) $\approx$ 2.0) and slope distributions statistically agreeing with observed strong lens galaxies in SLACS and BELLS. Agreement is only slightly worse with the more heterogeneous SL2S lens galaxy sample. Observations of strong-lens selected galaxies thus appear to favor models with relatively weak feedback in massive galaxies., Comment: Accepted in MNRAS

Radio detection of Ultra High Energetic Cosmic Rays and Neutrinos (UHECRv) which hit the Moon has been investigated in recent years. In preparation for near-future lunar science missions, we discuss technical requirements for radio experiments onboard lunar orbiters or on a lunar lander. We also develop an analysis of UHECRv aperture by including UHECv events occurring in the sub-layers of lunar regolith. It is verified that even using a single antenna onboard lunar orbiters or a few meters above the Moon's surface, dozens of lunar UHECRv events are detectable for one-year of observation at energy levels of $10^{18}$ eV to $10^{23}$ eV. Furthermore, it is shown that an antenna 3 meters above the Moon's surface could detect lower energy lunar UHECR events at the level of $10^{15}$ eV to $10^{18}$ eV which might not be detectable from lunar orbiters or ground-based observations., Comment: Accepted for publications in Advances in Space Research(ASR), 27 pages, 17 figures

Time-delay strong lensing provides a unique way to directly measure the Hubble constant ($H_{0}$). The precision of the $H_{0}$ measurement depends on the uncertainties in the time-delay measurements, the mass distribution of the main deflector(s), and the mass distribution along the line of sight. Tie and Kochanek (2018) have proposed a new microlensing effect on time delays based on differential magnification of the coherent accretion disc variability of the lensed quasar. If real, this effect could significantly broaden the uncertainty on the time delay measurements by up to $30\%$ for lens systems such as PG1115+080, which have relatively short time delays and monitoring over several different epochs. In this paper we develop a new technique that uses the time-delay ratios and simulated microlensing maps within a Bayesian framework in order to limit the allowed combinations of microlensing delays and thus to lessen the uncertainties due to the proposed effect. We show that, under the assumption of Tie and Kochanek (2018), the uncertainty on the time-delay distance ($D_{\Delta t}$, which is proportional to 1/$H_{0}$) of short time-delay ($\sim18$ days) lens, PG1115+080, increases from $\sim7\%$ to $\sim10\%$ by simultaneously fitting the three time-delay measurements from the three different datasets across twenty years, while in the case of long time-delay ($\sim90$ days) lens, the microlensing effect on time delays is negligible as the uncertainty on $D_{\Delta t}$ of RXJ1131-1231 only increases from $\sim2.5\%$ to $\sim2.6\%$.

Having an accurate calibration method is crucial for any scientific research done by a radio telescope. The next generation radio telescopes such as the Square Kilometre Array (SKA) will have a large number of receivers which will produce exabytes of data per day. In this paper we propose new direction-dependent and independent calibration algorithms that, while requiring much less storage during calibration, converge very fast. The calibration problem can be formulated as a non-linear least square optimization problem. We show that combining a block-LDU decomposition with Gauss-Newton iterations produces systems of equations with convergent matrices. This allows significant reduction in complexity per iteration and very fast converging algorithms. We also discuss extensions to direction-dependent calibration. The proposed algorithms are evaluated using simulations., Comment: EUSIPCO 2018 conference

In this paper we introduce the SEAGLE (i.e. Simulating EAGLE LEnses) program, that approaches the study of galaxy formation through strong gravitational lensing, using a suite of high-resolution hydrodynamic simulations, Evolution and Assembly of GaLaxies and their Environments (EAGLE) project. We introduce the simulation and analysis pipeline and present the first set of results from our analysis of early-type galaxies. We identify and extract an ensemble of simulated lens galaxies and use the GLAMER ray-tracing lensing code to create mock lenses similar to those observed in the SLACS and SL2S surveys, using a range of source parameters and galaxy orientations, including observational effects such as the Point-Spread-Function (PSF), pixelization and noise levels, representative of single-orbit observations with the Hubble Space Telescope (HST) using the ACS-F814W filter. We subsequently model these mock lenses using the code LENSED, treating them in the same way as observed lenses. We also estimate the mass model parameters directly from the projected surface mass density of the simulated galaxy, using an identical mass model family. We perform a three-way comparison of all the measured quantities with real lenses. We find the average total density slope of EAGLE lenses, $t=2.26\; (0.25\; \rm{rms})$ to be higher than SL2S, $t=2.16$ or SLACS, $t=2.08$. We find a very strong correlation between the external shear ($\gamma$) and the complex ellipticity ($\epsilon$), with $\gamma \sim \epsilon/4$. This correlation indicates a degeneracy in the lens mass modeling. We also see a dispersion between lens modeling and direct fitting results, indicating systematical biases., Comment: Accepted in MNRAS, Highlights: Figure- 6, 12, 13

Large scale imaging surveys will increase the number of galaxy-scale strong lensing candidates by maybe three orders of magnitudes beyond the number known today. Finding these rare objects will require picking them out of at least tens of millions of images and deriving scientific results from them will require quantifying the efficiency and bias of any search method. To achieve these objectives automated methods must be developed. Because gravitational lenses are rare objects reducing false positives will be particularly important. We present a description and results of an open gravitational lens finding challenge. Participants were asked to classify 100,000 candidate objects as to whether they were gravitational lenses or not with the goal of developing better automated methods for finding lenses in large data sets. A variety of methods were used including visual inspection, arc and ring finders, support vector machines (SVM) and convolutional neural networks (CNN). We find that many of the methods will be easily fast enough to analyse the anticipated data flow. In test data, several methods are able to identify upwards of half the lenses after applying some thresholds on the lens characteristics such as lensed image brightness, size or contrast with the lens galaxy without making a single false-positive identification. This is significantly better than direct inspection by humans was able to do. (abridged), Comment: 22 pages, 16 figures, 4 tables, accepted version for A&A

We present the full sample of 118 galaxy-scale strong-lens candidates in the Sloan Lens ACS (SLACS) Survey for the Masses (S4TM) Survey, which are spectroscopically selected from the final data release of the Sloan Digital Sky Survey. Follow-up Hubble Space Telescope (HST) imaging observations confirm that 40 candidates are definite strong lenses with multiple lensed images. The foreground lens galaxies are found to be early-type galaxies (ETGs) at redshifts 0.06 to 0.44, and background sources are emission-line galaxies at redshifts 0.22 to 1.29. As an extension of the SLACS Survey, the S4TM Survey is the first attempt to preferentially search for strong-lens systems with relatively lower lens masses than those in the pre-existing strong-lens samples. By fitting HST data with a singular isothermal ellipsoid model, we find total projected mass within the Einstein radius of the S4TM strong-lens sample ranges from $3 \times10^{10} M_{\odot}$ to $2 \times10^{11} M_{\odot}$. In [Shu15], we have derived the total stellar mass of the S4TM lenses to be $5 \times10^{10} M_{\odot}$ to $1 \times10^{12} M_{\odot}$. Both total enclosed mass and stellar mass of the S4TM lenses are on average almost a factor of 2 smaller than those of the SLACS lenses, which also represent typical mass scales of the current strong-lens samples. The extended mass coverage provided by the S4TM sample can enable a direct test, with the aid of strong lensing, for transitions in scaling relations, kinematic properties, mass structure, and dark-matter content trends of ETGs at intermediate-mass scales as noted in previous studies., Comment: 13 pages, 3 figures, very minor edits to match the ApJ-published version

In the last decade the detection of individual massive dark matter sub-halos has been possible using potential correction formalism in strong gravitational lens imaging. Here we propose a statistical formalism to relate strong gravitational lens surface brightness anomalies to the lens potential fluctuations arising from dark matter distribution in the lens galaxy. We consider these fluctuations as a Gaussian random field in addition to the unperturbed smooth lens model. This is very similar to weak lensing formalism and we show that in this way we can measure the power spectrum of these perturbations to the potential. We test the method by applying it to simulated mock lenses of different geometries and by performing an MCMC analysis of the theoretical power spectra. This method can measure density fluctuations in early type galaxies on scales of 1-10 kpc at typical rms-levels of a percent, using a single lens system observed with the Hubble Space Telescope with typical signal-to-noise ratios obtained in a single orbit.

Direct detection of the Epoch of Reionization (EoR) via the red-shifted 21-cm line will have unprecedented implications on the study of structure formation in the infant Universe. To fulfill this promise, current and future 21-cm experiments need to detect this weak EoR signal in the presence of foregrounds that are several orders of magnitude larger. This requires extreme noise control and improved wide-field high dynamic-range imaging techniques. We propose a new imaging method based on a maximum likelihood framework which solves for the interferometric equation directly on the sphere, or equivalently in the $uvw$-domain. The method uses the one-to-one relation between spherical waves and spherical harmonics (SpH). It consistently handles signals from the entire sky, and does not require a $w$-term correction. The spherical-harmonics coefficients represent the sky-brightness distribution and the visibilities in the $uvw$-domain, and provide a direct estimate of the spatial power spectrum. Using these spectrally-smooth SpH coefficients, bright foregrounds can be removed from the signal, including their side-lobe noise, which is one of the limiting factors in high dynamics range wide-field imaging. Chromatic effects causing the so-called "wedge" are effectively eliminated (i.e. deconvolved) in the cylindrical ($k_{\perp}, k_{\parallel}$) power spectrum, compared to a power spectrum computed directly from the images of the foreground visibilities where the wedge is clearly present. We illustrate our method using simulated LOFAR observations, finding an excellent reconstruction of the input EoR signal with minimal bias., Comment: 13 pages, 8 figures. Replaced to match accepted MNRAS version; few typos corrected & textual clarification added (no changes to results)

This paper is the first in a series of papers describing the impact of antenna instrumental artefacts on the 21-cm cosmology experiments to be carried out by the low frequency instrument (SKA1-LOW) of the Square Kilometre Array telescope (SKA), i.e., the Cosmic Dawn (CD) and the Epoch of Reionization (EoR). The smoothness of the passband response of the current log-periodic antenna being developed for the SKA1-LOW is analyzed using numerical electromagnetic simulations. The amplitude variations over the frequency range are characterized using low-order polynomials defined locally, in order to study the impact of the passband smoothness in the instrument calibration and CD/EoR Science. A solution is offered to correct a fast ripple found at 60~MHz during a test campaign at the SKA site at the Murchison Radio-astronomy Observatory, Western Australia in September 2015 with a minor impact on the telescope's performance and design. A comparison with the Hydrogen Epoch of Reionization Array antenna is also shown demonstrating the potential use of the SKA1-LOW antenna for the Delay Spectrum technique to detect the EoR., Comment: 11 pages, 12 figures, accepted for publication (MNRAS)

We introduce a novel technique, called "granulometry", to characterize and recover the mean size and the size distribution of HII regions from 21-cm tomography. The technique is easy to implement, but places the previously not very well defined concept of morphology on a firm mathematical foundation. The size distribution of the cold spots in 21-cm tomography can be used as a direct tracer of the underlying probability distribution of HII region sizes. We explore the capability of the method using large-scale reionization simulations and mock observational data cubes while considering capabilities of SKA1-low and a future extension to SKA2. We show that the technique allows the recovery of the HII region size distribution with a moderate signal-to-noise ratio from wide-field imaging ($\rm SNR\lesssim3$), for which the statistical uncertainty is sample variance dominated. We address the observational requirements on the angular resolution, the field-of-view, and the thermal noise limit for a successful measurement. To achieve a full scientific return from 21-cm tomography and to exploit a synergy with 21-cm power spectra, we suggest an observing strategy using wide-field imaging (several tens of square degrees) by an interferometric mosaicking/multi-beam observation with additional intermediate baselines (~2-4 km)., Comment: 19 pages, 19 figures, the version accepted in MNRAS

Strong gravitational lenses with measured time delays between the multiple images allow a direct measurement of the time-delay distance to the lens, and thus a measure of cosmological parameters, particularly the Hubble constant, $H_{0}$. We present a blind lens model analysis of the quadruply-imaged quasar lens HE 0435-1223 using deep Hubble Space Telescope imaging, updated time-delay measurements from the COSmological MOnitoring of GRAvItational Lenses (COSMOGRAIL), a measurement of the velocity dispersion of the lens galaxy based on Keck data, and a characterization of the mass distribution along the line of sight. HE 0435-1223 is the third lens analyzed as a part of the $H_{0}$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) project. We account for various sources of systematic uncertainty, including the detailed treatment of nearby perturbers, the parameterization of the galaxy light and mass profile, and the regions used for lens modeling. We constrain the effective time-delay distance to be $D_{\Delta t} = 2612_{-191}^{+208}~\mathrm{Mpc}$, a precision of 7.6%. From HE 0435-1223 alone, we infer a Hubble constant of $H_{0} = 73.1_{-6.0}^{+5.7}~\mathrm{km~s^{-1}~Mpc^{-1}}$ assuming a flat $\Lambda$CDM cosmology. The cosmographic inference based on the three lenses analyzed by H0LiCOW to date is presented in a companion paper (H0LiCOW Paper V)., Comment: Accepted for publication in MNRAS; 22 pages, 9 figures, 6 tables