Source author record

F. Courbin

F. Courbin appears in the imported research catalog. Authorship, coauthor and topic links are available while profile ownership is still unclaimed.

ResearcherUnclaimed source record

astro-ph.CO astro-ph.GA astro-ph.IM astro-ph.HE astro-ph

Catalog footprint

What is connected

62works

5topics

4close collaborators

Actions

Connect this record

Open graph Browse works

Inspect adjacent papers, topics, institutions and collaborators without losing the researcher page.

Building this map preview

BZPEER is loading the nearby papers, people, topics and institutions for this page.

preprint2026arXiv

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Intrinsic alignments of galaxies are potentially a major contaminant of cosmological analyses of weak gravitational lensing. We construct a semi-analytic model of galaxy ellipticities and alignments in the \Euclid Flagship simulation to predict this contamination in Euclid's weak lensing observations. Galaxy shapes and orientations are determined by the corresponding properties of the host haloes in the underlying $N$-body simulation, as well as the relative positions of galaxies within their halo. Alignment strengths are moderated via stochastic misalignments, separately for central and satellite galaxies and conditional on the galaxy's redshift, luminosity, and rest-frame colour. The resulting model is calibrated against galaxy ellipticity statistics from the COSMOS Survey, selected alignment measurements based on Sloan Digital Sky Survey samples, and galaxy orientations extracted from the Horizon-AGN hydrodynamic simulation at redshift $z=1$. The best-fit model has a total of 12 alignment parameters and generally reproduces the calibration data sets well within the $1σ$ statistical uncertainties of the observations and the \flagship simulation, with notable exceptions for the most luminous sub-samples on small physical scales. The statistical power of the calibration data and the volume of the single \flagship realisation are still too small to provide informative prior ranges for intrinsic alignment amplitudes in relevant galaxy samples. As a first application, we predict that \Euclid end-of-mission tomographic weak gravitational lensing two-point statistics are modified by up to order $10\,\%$ due to intrinsic alignments.

preprint2026arXiv

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

The Euclid satellite will measure spectroscopic redshifts for tens of millions of emission-line galaxies. In the context of Stage-IV surveys, the 3-dimensional clustering of galaxies plays a key role in providing cosmological constraints. In this paper, we conduct a model comparison for the multipole moments of the galaxy 2-point correlation function (2PCF) in redshift space. We test state-of-the-art models, in particular the effective field theory of large-scale structure (EFT), one based on the velocity difference generating function (VDG$_{\infty}$), and different variants of Lagrangian perturbation theory (LPT) models, such as convolutional LPT (CLPT) and its effective-field-theory extension (CLEFT). We analyse the first three even multipoles of the 2PCF in the Flagship 1 simulation, which consists of four snapshots at $z\in\{0.9,1.2,1.5,1.8\}$. We study both template-fitting and full-shape approaches and find that with the template-fitting approach, only the VDG$_{\infty}$ model is able to reach a minimum fitting scale of $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ at $z=0.9$ without biasing the recovered parameters. Indeed, the EFT model becomes inaccurate already at $s_{\rm min}=30\,h^{-1}\,{\rm Mpc}$. Conversely, in the full-shape analysis, the CLEFT and VDG$_{\infty}$ models perform similarly well, but only the CLEFT model can reach $s_{\rm min}=20\,h^{-1}\,{\rm Mpc}$ while the VDG$_{\infty}$ model is unbiased down to $s_{\rm min}=25\,h^{-1}\,{\rm Mpc}$ at the lowest redshift. Overall, in order to achieve the accuracy required by Euclid, non-perturbative modelling such as in the VDG$_{\infty}$ or CLEFT models should be considered. At $z=1.8$, the CLPT model is sufficient to describe the data with high figure of merit. This comparison selects baseline models that perform best in ideal conditions and sets the stage for an optimal analysis of Euclid data in configuration space.

preprint2026arXiv

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

We model intrinsic alignments (IA) in Euclid's Flagship simulation to investigate its impact on Euclid's weak lensing signal. Our IA implementation in the Flagship simulation takes into account photometric properties of galaxies as well as their dark matter host halos. We compare simulations against theory predictions, determining the parameters of two of the most widely used IA models: the Non Linear Alignment (NLA) and the Tidal Alignment and Tidal Torquing (TATT) models. We measure the amplitude of the simulated IA signal as a function of galaxy magnitude and colour in the redshift range $0.1<z<2.1$. We find that both NLA and TATT can accurately describe the IA signal in the simulation down to scales of $6$-$7 \,h^{-1}\,$Mpc. We measure alignment amplitudes for red galaxies comparable to those of the observations, with samples not used in the calibration procedure. For blue galaxies, our constraints are consistent with zero alignments in our first redshift bin $0.1 < z < 0.3$, but we detect a non-negligible signal at higher redshift, which is, however, consistent with the upper limits set by observational constraints. Additionally, several hydrodynamical simulations predict alignment for spiral galaxies, in agreement with our findings. Finally, the evolution of alignment with redshift is realistic and comparable to that determined in the observations. However, we find that the commonly adopted redshift power-law for IA fails to reproduce the simulation alignments above $z=1.1$. A significantly better agreement is obtained when a luminosity dependence is included, capturing the intrinsic luminosity evolution with redshift in magnitude-limited surveys. We conclude that the Flagship IA simulation is a useful tool for translating current IA constraints into predictions for IA contamination of Euclid-like samples.

preprint2026arXiv

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

We present a detailed analysis of the old, extended stellar component of the Local Group dwarf galaxy IC 10 using deep resolved-star photometry in the VIS and NISP bands of the Euclid Early Release Observations. Leveraging Euclid's unique combination of wide field of view and high spatial resolution, we trace red giant branch (RGB) stars out to $\sim$8 kpc from the galaxy centre, reaching azimuthally-averaged surface brightness levels as faint as $μ_{HE}\sim$29 mag arcsec$^{-2}$. Our analysis reveals that IC 10's stellar distribution is significantly more extended than previously thought. After correcting for foreground extinction and subtracting contamination from Milky Way stars and background galaxies, we derive a radial stellar density profile from RGB star counts. The profile shows a marked flattening beyond $\sim$5 kpc, and is best fit by a two-component (Sersic + exponential) model, yielding a total stellar mass in old (age $>$1 Gyr) stars of $M_{\star}=(6.7-8.1)\times10^8 M_{\odot}$. The origin of the outer stellar component is unclear. It may be accreted, even possibly associated with the counter-rotating HI gas in the outer regions of IC 10, or it may represent an ancient in-situ stellar halo. We tentatively detect two symmetric stellar overdensities at the edge of our imagery. These roughly align with the direction of IC 10's orbit around M31, suggesting that they may be signatures of tidal stripping. As part of our analysis, we derive a new distance to IC 10 based on the RGB tip, finding $D=762\pm 20$ kpc and the distance modulus is $(m-M)_0=24.41\pm 0.05$.

preprint2026arXiv

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Weak lensing surveys require accurate correction for the point spread function (PSF) when measuring galaxy shapes. For a diffraction-limited PSF, as arises in space-based missions, this correction depends on each galaxy SED. In the Euclid mission, galaxy SED reconstruction, a tasks of the photometric-redshift processing function (PHZ PF), relies on broad- and medium-band ancillary photometry. The limited wavelength sampling of the Euclid VIS passband and signal-to-noise ratio may affect the reconstruction accuracy and translate into biases in the weak lensing measurements. In this study, we present the methodology, which is employed in the Euclid PHZ PF, for reconstructing galaxy SEDs at 55 wavelengths, sampling the VIS passband every 10 nm, and we assess whether it fulfils the accuracy requirements imposed on the Euclid PSF model. We employ both physics- and data-driven methods, focusing on a new approach of template-based flux correction and Gaussian processes, and we introduce an SED metric whose bias propagates into PSF quadrupole moment errors. Our findings demonstrate that Gaussian processes and template fitting meet the requirements only in specific, but complementary, redshift intervals. We therefore propose a hybrid approach, which leverages both methods. This solution proves to be effective in meeting the Euclid accuracy requirements for most of the redshift range of the survey. Finally, we investigate the impact on the SED reconstruction of a new set of 16 evenly-spaced medium-band filters for the Subaru telescope, providing quasi-spectroscopic coverage of the VIS passband. This study shows promising results, ensuring accurate SED reconstruction and meeting the mission PSF requirements. This work thus provides not only the methodological foundation of galaxy SED reconstruction in the Euclid PHZ PF, but also a roadmap for future improvements using a new medium-band survey.

preprint2026arXiv

Time-delay and lens galaxy redshift in the doubly imaged quasar PS J2305+3714

We present results from a seven-season (2018-2024) monitoring campaign of the gravitationally lensed quasar system PS J2305+3714 using the 1.5-m Maidanak Telescope in the optical R-band. From these data, the amplitude of possible microlensing variability does not exceed 10 mmag on a 7-yr timescale. Additionally, we measure a time delay of $t_{AB}$ = 52.2$\pm$2.5 days, with the brighter image A leading. Long-slit \WHT spectroscopy refines the quasar redshift to $z_s$ = 1.791 and provides the first measurement of the lens redshift, $z_d$ = 0.473. The flux ratios of the quasar images in the MgII $λ$2800 emission line and in the adjacent continuum are nearly identical, indicating minimal microlensing effects in the spectral domain, which is consistent with the very weak microlensing signal in the time domain. Using precise astrometry from recent HST imaging and the MgII flux ratio, we also built two simple mass models for the lens system. The close agreement between the measured delay and those predicted by the mass models, measured redshifts, and a concordance $Λ$CDM cosmology, confirms the robustness of our results and highlights PS J2305+3714 as a promising system for future time-delay cosmography.

preprint2025arXiv

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

We present two extensive sets of 3500+1000 simulations of dark matter haloes on the past light cone, and two corresponding sets of simulated (`mock') galaxy catalogues that represent the Euclid spectroscopic sample. The simulations were produced with the latest version of the PINOCCHIO code, and provide the largest, public set of simulated skies. Mock galaxy catalogues were obtained by populating haloes with galaxies using an halo occupation distribution (HOD) model extracted from the Flagship galaxy catalogue provided by Euclid Collaboration. The Geppetto set of 3500 simulated skies was obtained by tiling a 1.2 Gpc/h box to cover a light-cone whose sky footprint is a circle of 30 deg radius, for an area of 2763 deg$^2$ and a minimum halo mass of $1.5\times10^{11}$ Msun/h. The relatively small box size makes this set unfit for measuring very large scales. The EuclidLargeBox set consists of 1000 simulations of 3.38 Gpc/h, with the same mass resolution and a footprint that covers half of the sky, excluding the Milky Way zone of avoidance. From this we produced a set of 1000 EuclidLargeMocks on the 30 deg radius footprint, whose comoving volume is fully contained in the simulation box. We validated the two sets of catalogues by analysing number densities, power spectra, and 2-point correlation functions, showing that the Flagship spectroscopic catalogue is consistent with being one of the realisations of the simulated sets, although we noticed small deviations limited to the quadrupole at k>0.2 h/Mpc. We show cosmological parameter inference from these catalogues and demonstrate that using one realisation of EuclidLargeMocks in place of the Flagship mock produces the same posteriors, to within the expected shift given by sample variance. These simulated skies will be used for the galaxy clustering analysis of Euclid's Data Release 1 (DR1).

preprint2023arXiv

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Next generation telescopes, like Euclid, Rubin/LSST, and Roman, will open new windows on the Universe, allowing us to infer physical properties for tens of millions of galaxies. Machine learning methods are increasingly becoming the most efficient tools to handle this enormous amount of data, because they are often faster and more accurate than traditional methods. We investigate how well redshifts, stellar masses, and star-formation rates (SFR) can be measured with deep learning algorithms for observed galaxies within data mimicking the Euclid and Rubin/LSST surveys. We find that Deep Learning Neural Networks and Convolutional Neutral Networks (CNN), which are dependent on the parameter space of the training sample, perform well in measuring the properties of these galaxies and have a better accuracy than methods based on spectral energy distribution fitting. CNNs allow the processing of multi-band magnitudes together with $H_{\scriptscriptstyle\rm E}$-band images. We find that the estimates of stellar masses improve with the use of an image, but those of redshift and SFR do not. Our best results are deriving i) the redshift within a normalised error of less than 0.15 for 99.9$\%$ of the galaxies with S/N>3 in the $H_{\scriptscriptstyle\rm E}$-band; ii) the stellar mass within a factor of two ($\sim0.3 \rm dex$) for 99.5$\%$ of the considered galaxies; iii) the SFR within a factor of two ($\sim0.3 \rm dex$) for $\sim$70$\%$ of the sample. We discuss the implications of our work for application to surveys as well as how measurements of these galaxy parameters can be improved with deep learning.

preprint2023arXiv

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

This work focuses on the pilot run of a simulation campaign aimed at investigating the spectroscopic capabilities of the Euclid Near-Infrared Spectrometer and Photometer (NISP), in terms of continuum and emission line detection in the context of galaxy evolutionary studies. To this purpose we constructed, emulated, and analysed the spectra of 4992 star-forming galaxies at $0.3 \leq z \leq 2.5$ using the NISP pixel-level simulator. We built the spectral library starting from public multi-wavelength galaxy catalogues, with value-added information on spectral energy distribution (SED) fitting results, and from Bruzual and Charlot (2003) stellar population templates. Rest-frame optical and near-IR nebular emission lines were included using empirical and theoretical relations. We inferred the 3.5$σ$ NISP red grism spectroscopic detection limit of the continuum measured in the $H$ band for star-forming galaxies with a median disk half-light radius of \ang{;;0.4} at magnitude $H= 19.5\pm0.2\,$AB$\,$mag for the Euclid Wide Survey and at $H = 20.8\pm0.6\,$AB$\,$mag for the Euclid Deep Survey. We found a very good agreement with the red grism emission line detection limit requirement for the Wide and Deep surveys. We characterised the effect of the galaxy shape on the detection capability of the red grism and highlighted the degradation of the quality of the extracted spectra as the disk size increases. In particular, we found that the extracted emission line signal to noise ratio (SNR) drops by $\sim\,$45$\%$ when the disk size ranges from \ang{;;0.25} to \ang{;;1}. These trends lead to a correlation between the emission line SNR and the stellar mass of the galaxy and we demonstrate the effect in a stacking analysis unveiling emission lines otherwise too faint to detect.

preprint2022arXiv

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

We investigate the importance of lensing magnification for estimates of galaxy clustering and its cross-correlation with shear for the photometric sample of Euclid. Using updated specifications, we study the impact of lensing magnification on the constraints and the shift in the estimation of the best fitting cosmological parameters that we expect if this effect is neglected. We follow the prescriptions of the official Euclid Fisher matrix forecast for the photometric galaxy clustering analysis and the combination of photometric clustering and cosmic shear. The slope of the luminosity function (local count slope), which regulates the amplitude of the lensing magnification, and the galaxy bias have been estimated from the Euclid Flagship simulation.We find that magnification significantly affects both the best-fit estimation of cosmological parameters and the constraints in the galaxy clustering analysis of the photometric sample. In particular, including magnification in the analysis reduces the 1$σ$ errors on $Ω_{\text{m},0}, w_{0}, w_a$ at the level of 20-35%, depending on how well we will be able to independently measure the local count slope. In addition, we find that neglecting magnification in the clustering analysis leads to shifts of up to 1.6$σ$ in the best-fit parameters. In the joint analysis of galaxy clustering, cosmic shear, and galaxy-galaxy lensing, magnification does not improve precision, but it leads to an up to 6$σ$ bias if neglected. Therefore, for all models considered in this work, magnification has to be included in the analysis of galaxy clustering and its cross-correlation with the shear signal ($3\times2$pt analysis) for an accurate parameter estimation.

preprint2022arXiv

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

The Complete Calibration of the Color-Redshift Relation survey (C3R2) is a spectroscopic programme designed to empirically calibrate the galaxy color-redshift relation to the Euclid depth (I_E=24.5), a key ingredient for the success of Stage IV dark energy projects based on weak lensing cosmology. A spectroscopic calibration sample as representative as possible of the galaxies in the Euclid weak lensing sample is being collected, selecting galaxies from a self-organizing map (SOM) representation of the galaxy color space. Here, we present the results of a near-infrared H- and K-bands spectroscopic campaign carried out using the LUCI instruments at the LBT. For a total of 251 galaxies, we present new highly-reliable redshifts in the 1.3<= z <=1.7 and 2<= z<=2.7 ranges. The newly-determined redshifts populate 49 SOM cells which previously contained no spectroscopic measurements and almost double the occupation numbers of an additional 153 SOM cells. A final optical ground-based observational effort is needed to calibrate the missing cells in particular in the redshift range 1.7<= z<=2.7 that lack spectroscopic calibration. In the end, Euclid itself will deliver telluric-free NIR spectra that can complete the calibration.

preprint2022arXiv

Euclid preparation. XVIII. The NISP photometric system

Euclid will be the first space mission to survey most of the extragalactic sky in the 0.95-2.02 $μ$m range, to a 5$σ$ point-source median depth of 24.4 AB mag. This unique photometric data set will find wide use beyond Euclid's core science. In this paper, we present accurate computations of the Euclid Y_E, J_E and H_E passbands used by the Near-Infrared Spectrometer and Photometer (NISP), and the associated photometric system. We pay particular attention to passband variations in the field of view, accounting among others for spatially variable filter transmission, and variations of the angle of incidence on the filter substrate using optical ray tracing. The response curves' cut-on and cut-off wavelengths - and their variation in the field of view - are determined with 0.8 nm accuracy, essential for the photometric redshift accuracy required by Euclid. After computing the photometric zeropoints in the AB mag system, we present linear transformations from and to common ground-based near-infrared photometric systems, for normal stars, red and brown dwarfs, and galaxies separately. A Python tool to compute accurate magnitudes for arbitrary passbands and spectral energy distributions is provided. We discuss various factors from space weathering to material outgassing that may slowly alter Euclid's spectral response. At the absolute flux scale, the Euclid in-flight calibration program connects the NISP photometric system to Hubble Space Telescope spectrophotometric white dwarf standards; at the relative flux scale, the chromatic evolution of the response is tracked at the milli-mag level. In this way, we establish an accurate photometric system that is fully controlled throughout Euclid's lifetime.

preprint2022arXiv

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

The standard cosmological model is based on the fundamental assumptions of a spatially homogeneous and isotropic universe on large scales. An observational detection of a violation of these assumptions at any redshift would immediately indicate the presence of new physics. We quantify the ability of the Euclid mission, together with contemporary surveys, to improve the current sensitivity of null tests of the canonical cosmological constant $Λ$ and the cold dark matter (LCDM) model in the redshift range $0<z<1.8$. We considered both currently available data and simulated Euclid and external data products based on a LCDM fiducial model, an evolving dark energy model assuming the Chevallier-Polarski-Linder (CPL) parameterization or an inhomogeneous Lema\^ıtre-Tolman-Bondi model with a cosmological constant $Λ$ (LLTB), and carried out two separate but complementary analyses: a machine learning reconstruction of the null tests based on genetic algorithms, and a theory-agnostic parametric approach based on Taylor expansion and binning of the data, in order to avoid assumptions about any particular model. We find that in combination with external probes, Euclid can improve current constraints on null tests of the LCDM by approximately a factor of three when using the machine learning approach and by a further factor of two in the case of the parametric approach. However, we also find that in certain cases, the parametric approach may be biased against or missing some features of models far from LCDM. Our analysis highlights the importance of synergies between Euclid and other surveys. These synergies are crucial for providing tighter constraints over an extended redshift range for a plethora of different consistency tests of some of the main assumptions of the current cosmological paradigm.

preprint2021arXiv

Constraining quasar structure using high-frequency microlensing variations and continuum reverberation

Gravitational microlensing is a powerful tool to probe the inner structure of strongly lensed quasars and to constrain parameters of the stellar mass function of lens galaxies. This is done by analysing microlensing light curves between the multiple images of strongly lensed quasars, under the influence of three main variable components: 1- the continuum flux of the source, 2- microlensing by stars in the lens galaxy and 3- reverberation of the continuum by the Broad Line Region (BLR). The latter, ignored by state-of-the-art microlensing techniques, can introduce high-frequency variations which we show carry information on the BLR size. We present a new method which includes all these components simultaneously and fits the power spectrum of the data in the Fourier space, rather than the observed light curve itself. In this new framework, we analyse COSMOGRAIL light curves of the two-image system QJ0158-4325 known to display high-frequency variations. Using exclusively the low frequency part of the power spectrum our constraint on the accretion disk radius agrees with the thin disk model estimate and previous work that fit the microlensing light curves in real space. However, if we also take into account the high-frequency variations, the data favour significantly smaller disk sizes than previous microlensing measurements. In this case, our results are in agreement with the thin disk model prediction only if we assume very low mean masses for the microlens population, i.e. <M> = 0.01 $M_\odot$. Eventually, including the differentially microlensed continuum reverberation by the BLR successfully explains the high frequencies without requiring such low mass microlenses. This allows us to measure, for the first time, the size of the BLR using single-band photometric monitoring, $R_{BLR}$ = $1.6^{+1.5}_{-0.8}\times 10^{17}$cm, in agreement with estimates using the BLR size-luminosity relation.

preprint2021arXiv

Discovery of Strongly Lensed Quasars in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

We report the discovery of five new doubly-imaged lensed quasars from the first 2500 square degrees of the ongoing Canada-France Imaging Survey (CFIS), which is a component of the Ultraviolet Near Infrared Optical Northern Survey (UNIONS), selected from initial catalogues of either Gaia pairs or MILLIQUAS quasars. We take advantage of the deep, 0.6'' median-seeing $r$-band imaging of CFIS to confirm the presence of multiple point sources with similar colour of $u-r$, via convolution of the Laplacian of the point spread function. Requiring similar-colour point sources with flux ratios less than 2.5 mag in $r$-band, reduces the number of candidates from 256314 to 7815. After visual inspection we obtain 30 high-grade candidates, and prioritise spectroscopic follow-up for those showing signs of a lensing galaxy upon subtraction of the point sources. We obtain long-slit spectra for 18 candidates with ALFOSC on the 2.56-m Nordic Optical Telescope (NOT), confirming five new doubly lensed quasars with $1.21<z<3.36$ and angular separations from 0.8'' to 2.5''. One additional system is a probable lensed quasar based on the CFIS imaging and existing SDSS spectrum. We further classify six objects as nearly identical quasars -- still possible lenses but without the detection of a lensing galaxy. Given our recovery rate ($83\%$) of existing optically bright lenses within the CFIS footprint, we expect that a similar strategy, coupled with $u-r$ colour-selection from CFIS alone, will provide an efficient and complete discovery of small-separation lensed quasars of source redshifts below $z=2.7$ within the CFIS $r$-band magnitude limit of 24.1 mag.

preprint2021arXiv

Euclid preparation: I. The Euclid Wide Survey

Euclid is an ESA mission designed to constrain the properties of dark energy and gravity via weak gravitational lensing and galaxy clustering. It will carry out a wide area imaging and spectroscopy survey (EWS) in visible and near-infrared, covering roughly 15,000 square degrees of extragalactic sky on six years. The wide-field telescope and instruments are optimized for pristine PSF and reduced straylight, producing very crisp images. This paper presents the building of the Euclid reference survey: the sequence of pointings of EWS, Deep fields, Auxiliary fields for calibrations, and spacecraft movements followed by Euclid as it operates in a step-and-stare mode from its orbit around the Lagrange point L2. Each EWS pointing has four dithered frames; we simulate the dither pattern at pixel level to analyse the effective coverage. We use up-to-date models for the sky background to define the Euclid region-of-interest (RoI). The building of the reference survey is highly constrained from calibration cadences, spacecraft constraints and background levels; synergies with ground-based coverage are also considered. Via purposely-built software optimized to prioritize best sky areas, produce a compact coverage, and ensure thermal stability, we generate a schedule for the Auxiliary and Deep fields observations and schedule the RoI with EWS transit observations. The resulting reference survey RSD_2021A fulfills all constraints and is a good proxy for the final solution. Its wide survey covers 14,500 square degrees. The limiting AB magnitudes ($5σ$ point-like source) achieved in its footprint are estimated to be 26.2 (visible) and 24.5 (near-infrared); for spectroscopy, the H$_α$ line flux limit is $2\times 10^{-16}$ erg cm$^{-2}$ s$^{-1}$ at 1600 nm; and for diffuse emission the surface brightness limits are 29.8 (visible) and 28.4 (near-infrared) mag arcsec$^{-2}$.

preprint2021arXiv

J1721$+$8842: a gravitationally lensed binary quasar with a proximate damped Lyman-$α$ absorber

High-redshift binary quasars provide key insights into mergers and quasar activity, and are useful tools for probing the spatial kinematics and chemistry of galaxies along the line-of-sight. However, only three sub-10-kpc binaries have been confirmed above $z=1$. Gravitational lensing would provide a way to easily resolve such binaries, study them in higher resolution, and provide more sightlines, though the required alignment with a massive foreground galaxy is rare. Through image deconvolution of StanCam Nordic Optical Telescope (NOT) monitoring data, we reveal two further point sources in the known, $z \approx 2.38$, quadruply lensed quasar (quad), J1721+8842. An ALFOSC/NOT long-slit spectrum shows that the brighter of these two sources is a quasar with $z = 2.369 \pm 0.007$ based on the C III] line, while the C III] redshift of the quad is $z = 2.364 \pm 0.003$. Lens modelling using point source positions rules out a single source model, favouring an isothermal lens mass profile with two quasar sources separated by $\sim6.0$ kpc (0.73$^{\prime \prime}$) in projection. Given the resolving ability from lensing and current lensed quasar statistics, this discovery suggests a large population of undiscovered, unlensed sub-10-kpc binaries. We also analyse spectra of two images of the quad, showing narrow Ly$α$ emission within the trough of a proximate damped Ly$α$ absorber (PDLA). An apparent mismatch between the continuum and narrow line flux ratios provides a new potential tool for simultaneously studying microlensing and the quasar host galaxy. Signs of the PDLA are also seen in the second source, however a deeper spectrum is still required to confirm this. Thanks to the multiple lines-of-sight from lensing and two quasar sources, this system offers simultaneous sub-parsec and kpc-scale probes of a PDLA.

preprint2020arXiv

COSMOGRAIL XIX: Time delays in 18 strongly lensed quasars from 15 years of optical monitoring

We present the results of 15 years of monitoring lensed quasars, which was conducted by the COSMOGRAIL programme at the Leonhard Euler 1.2m Swiss Telescope. The decade-long light curves of 23 lensed systems are presented for the first time. We complement our data set with other monitoring data available in the literature to measure the time delays in 18 systems, among which nine reach a relative precision better than 15% for at least one time delay. To achieve this, we developed an automated version of the curve-shifting toolbox PyCS to ensure robust estimation of the time delay in the presence of microlensing, while accounting for the errors due to the imperfect representation of microlensing. We also re-analysed the previously published time delays of RX J1131$-$1231 and HE 0435$-$1223, by adding six and two new seasons of monitoring, respectively, and confirming the previous time-delay measurements. When the time delay measurement is possible, we corrected the light curves of the lensed images from their time delay and present the difference curves to highlight the microlensing signal contained in the data. To date, this is the largest sample of decade-long lens monitoring data, which is useful to measure $H_0$ and the size of quasar accretion discs with microlensing as well as to study quasar variability.

preprint2020arXiv

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

We investigate the environment and line of sight of the H0LiCOW lens B1608+656 using Subaru Suprime-Cam and the Hubble Space Telescope (HST) to perform a weak lensing analysis. We compare three different methods to reconstruct the mass map of the field, i.e. the standard Kaiser-Squires inversion coupled with inpainting and Gaussian or wavelet filtering, and $\tt{Glimpse}$ a method based on sparse regularization of the shear field. We find no substantial difference between the 2D mass reconstructions, but we find that the ground-based data is less sensitive to small-scale structures than the space-based observations. Marginalising over the results obtained with all the reconstruction techniques applied to the two available HST filters F606W and F814W, we estimate the external convergence, $κ_{\rm ext}$ at the position of B1608+656 is $κ_{\rm ext} = 0.11^{+0.06}_{-0.04}$, where the error bars corresponds respectively to the 16th and 84th quartiles. This result is compatible with previous estimates using the number-counts technique, suggesting that B1608+656 resides in an over-dense line of sight, but with a completely different technique. Using our mass reconstructions, we also compare the convergence at the position of several groups of galaxies in the field of B1608+656 with the mass measurements using various analytical mass profiles, and find that the weak lensing results favor truncated halo models.

preprint2020arXiv

HOLISMOKES -- I. Highly Optimised Lensing Investigations of Supernovae, Microlensing Objects, and Kinematics of Ellipticals and Spirals

We present the HOLISMOKES programme on strong gravitational lensing of supernovae as a probe of supernova (SN) physics and cosmology. We investigate the effects of microlensing on early-phase SN Ia spectra using four different SN explosion models, and find that within 10 rest-frame days after SN explosion, distortions of SN Ia spectra due to microlensing are typically negligible ($<$1% distortion within the 1$σ$ spread, and $\lesssim$10% distortion within the 2$σ$ spread). This shows great prospects of using lensed SNe Ia to obtain intrinsic early-phase SN spectra for deciphering SN Ia progenitors. As a demonstration of the usefulness of lensed SNe Ia for cosmology, we simulate a sample of mock lensed SN Ia systems that are expected to have accurate and precise time-delay measurements in the era of the Rubin Observatory Legacy Survey of Space and Time (LSST). Adopting realistic yet conservative uncertainties on their time-delay distances and lens angular diameter distances (of 6.6% and 5%, respectively), we find that a sample of 20 lensed SNe Ia would allow a constraint on the Hubble constant ($H_0$) with 1.3% uncertainty in the flat $Λ$CDM cosmology. We find a similar constraint on $H_0$ in an open $Λ$CDM cosmology, while the constraint degrades to $3\%$ in a flat $w$CDM cosmology. We anticipate lensed SNe to be an independent and powerful probe of SN physics and cosmology in the upcoming LSST era.

preprint2020arXiv

Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars

Strong-gravitational lens systems with quadruply-imaged quasars (quads) are unique probes to address several fundamental problems in cosmology and astrophysics. Although they are intrinsically very rare, ongoing and planned wide-field deep-sky surveys are set to discover thousands of such systems in the next decade. It is thus paramount to devise a general framework to model strong-lens systems to cope with this large influx without being limited by expert investigator time. We propose such a general modelling framework (implemented with the publicly available software Lenstronomy) and apply it to uniformly model three-band Hubble Space Telescope Wide Field Camera 3 images of 13 quads. This is the largest uniformly modelled sample of quads to date and paves the way for a variety of studies. To illustrate the scientific content of the sample, we investigate the alignment between the mass and light distribution in the deflectors. The position angles of these distributions are well-aligned, except when there is strong external shear. However, we find no correlation between the ellipticity of the light and mass distributions. We also show that the observed flux-ratios between the images depart significantly from the predictions of simple smooth models. The departures are strongest in the bluest band, consistent with microlensing being the dominant cause in addition to millilensing. Future papers will exploit this rich dataset in combination with ground based spectroscopy and time delays to determine quantities such as the Hubble constant, the free streaming length of dark matter, and the normalization of the initial stellar mass function.

preprint2020arXiv

STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354

We present a blind time-delay cosmographic analysis for the lens system DES J0408$-$5354. This system is extraordinary for the presence of two sets of multiple images at different redshifts, which provide the opportunity to obtain more information at the cost of increased modelling complexity with respect to previously analyzed systems. We perform detailed modelling of the mass distribution for this lens system using three band Hubble Space Telescope imaging. We combine the measured time delays, line-of-sight central velocity dispersion of the deflector, and statistically constrained external convergence with our lens models to estimate two cosmological distances. We measure the "effective" time-delay distance corresponding to the redshifts of the deflector and the lensed quasar $D_{Δt}^{\rm eff}=3382^{+146}_{-115}$ Mpc and the angular diameter distance to the deflector $D_{\rm d}=1711^{+376}_{-280}$ Mpc, with covariance between the two distances. From these constraints on the cosmological distances, we infer the Hubble constant $H_0 = 74.2^{+2.7}_{-3.0}$ km s$^{-1}$ Mpc$^{-1}$ assuming a flat $Λ$CDM cosmology and a uniform prior for $Ω_{\rm m}$ as $Ω_{\rm m} \sim \mathcal{U}(0.05, 0.5)$. This measurement gives the most precise constraint on $H_0$ to date from a single lens. Our measurement is consistent with that obtained from the previous sample of six lenses analyzed by the $H_0$ Lenses in COSMOGRAIL's Wellspring (H0LiCOW) collaboration. It is also consistent with measurements of $H_0$ based on the local distance ladder, reinforcing the tension with the inference from early Universe probes, for example, with 2.2$σ$ discrepancy from the cosmic microwave background measurement.

preprint2020arXiv

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

Time-delay cosmography of lensed quasars has achieved 2.4% precision on the measurement of the Hubble constant, $H_0$. As part of an ongoing effort to uncover and control systematic uncertainties, we investigate three potential sources: 1- stellar kinematics, 2- line-of-sight effects, and 3- the deflector mass model. To meet this goal in a quantitative way, we reproduced the H0LiCOW/SHARP/STRIDES (hereafter TDCOSMO) procedures on a set of real and simulated data, and we find the following. First, stellar kinematics cannot be a dominant source of error or bias since we find that a systematic change of 10% of measured velocity dispersion leads to only a 0.7% shift on $H_0$ from the seven lenses analyzed by TDCOSMO. Second, we find no bias to arise from incorrect estimation of the line-of-sight effects. Third, we show that elliptical composite (stars + dark matter halo), power-law, and cored power-law mass profiles have the flexibility to yield a broad range in $H_0$ values. However, the TDCOSMO procedures that model the data with both composite and power-law mass profiles are informative. If the models agree, as we observe in real systems owing to the "bulge-halo" conspiracy, $H_0$ is recovered precisely and accurately by both models. If the two models disagree, as in the case of some pathological models illustrated here, the TDCOSMO procedure either discriminates between them through the goodness of fit, or it accounts for the discrepancy in the final error bars provided by the analysis. This conclusion is consistent with a reanalysis of six of the TDCOSMO (real) lenses: the composite model yields $74.0^{+1.7}_{-1.8}$ $km.s^{-1}.Mpc^{-1}$, while the power-law model yields $H_0=74.2^{+1.6}_{-1.6}$ $km.s^{-1}.Mpc^{-1}$. In conclusion, we find no evidence of bias or errors larger than the current statistical uncertainties reported by TDCOSMO.

preprint2019arXiv

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2017/2018 follow-up campaign: Discovery of 10 lensed quasars and 10 quasar pairs

We report the results of the STRong lensing Insights from the Dark Energy Survey (STRIDES) follow-up campaign of the late 2017/early 2018 season. We obtained spectra of 65 lensed quasar candidates either with EFOSC2 on the NTT or ESI on Keck, which confirm 10 new gravitationally lensed quasars and 10 quasar pairs with similar spectra, but which do not show a lensing galaxy in DES images. Eight lensed quasars are doubly imaged with source redshifts between 0.99 and 2.90, one is triply imaged by a group (DESJ0345-2545, $z=1.68$), and one is quadruply imaged (quad: DESJ0053-2012, $z=3.8$). Singular isothermal ellipsoid models for the doubles, based on high-resolution imaging from SAMI on SOAR or NIRC2 on Keck, give total magnifications between 3.2 and 5.6, and Einstein radii between 0.49 and 1.97 arcseconds. After spectroscopic follow-up, we extract multi-epoch $grizY$ photometry of confirmed lensed quasars and contaminant quasar+star pairs from the first 4 years of DES data using parametric multi-band modelling, and compare variability in each system's components. By measuring the reduced $χ^2$ associated with fitting all epochs to the same magnitude, we find a simple cut on the less variable component that retains all confirmed lensed quasars, while removing 94 per cent of contaminant systems with stellar components. Based on our spectroscopic follow-up, this variability information can improve selection of lensed quasars and quasar pairs from 34-45 per cent to 51-70 per cent, with the majority of remaining contaminants being compact star-forming galaxies. Using mock lensed quasar lightcurves we demonstrate that selection based only on variability will over-represent the quad fraction by 10 per cent over a complete DES magnitude-limited sample (excluding microlensing differences), explained by the magnification bias and hence lower luminosity (more variable) sources in quads.

preprint2019arXiv

Twisted quasar light curves: implications for continuum reverberation mapping of accretion disks

With the advent of high-cadence and multi-band photometric monitoring facilities, continuum reverberation mapping is becoming of increasing importance to measure the physical size of quasar accretion disks. The method is based on the measurement of the time it takes for a signal to propagate from the center to the outer parts of the central engine, assuming the continuum light curve at a given wavelength has a time shift of the order of a few days with respect to light curves obtained at shorter wavelengths. We show that with high-quality light curves, this assumption is not valid anymore and that light curves at different wavelengths are not only shifted in time but also distorted: in the context of the lamp-post model and thin-disk geometry, the multi-band light curves are in fact convolved by a transfer function whose size increase with wavelength. We illustrate the effect with simulated light curves in the LSST ugrizy bands and examine the impact on the delay measurements when using three different methods, namely JAVELIN, CREAM, and PyCS. We find that current accretion disk sizes estimated from JAVELIN and PyCS are underestimated by $\sim30\%$ and that unbiased measurement are only obtained with methods that properly take the skewed transfer functions into account, as the CREAM code does. With the LSST-like light curves, we expect to achieve measurement errors below $5\%$ with typical 2-day photometric cadence.

preprint2016arXiv

ASTErIsM - Application of topometric clustering algorithms in automatic galaxy detection and classification

We present a study on galaxy detection and shape classification using topometric clustering algorithms. We first use the DBSCAN algorithm to extract, from CCD frames, groups of adjacent pixels with significant fluxes and we then apply the DENCLUE algorithm to separate the contributions of overlapping sources. The DENCLUE separation is based on the localization of pattern of local maxima, through an iterative algorithm which associates each pixel to the closest local maximum. Our main classification goal is to take apart elliptical from spiral galaxies. We introduce new sets of features derived from the computation of geometrical invariant moments of the pixel group shape and from the statistics of the spatial distribution of the DENCLUE local maxima patterns. Ellipticals are characterized by a single group of local maxima, related to the galaxy core, while spiral galaxies have additional ones related to segments of spiral arms. We use two different supervised ensemble classification algorithms, Random Forest, and Gradient Boosting. Using a sample of ~ 24000 galaxies taken from the Galaxy Zoo 2 main sample with spectroscopic redshifts, and we test our classification against the Galaxy Zoo 2 catalog. We find that features extracted from our pipeline give on average an accuracy of ~ 93%, when testing on a test set with a size of 20% of our full data set, with features deriving from the angular distribution of density attractor ranking at the top of the discrimination power.

preprint2016arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XV. Assessing the achievability and precision of time-delay measurements

COSMOGRAIL is a long-term photometric monitoring of gravitationally lensed QSOs aimed at implementing Refsdal's time-delay method to measure cosmological parameters, in particular H0. Given long and well sampled light curves of strongly lensed QSOs, time-delay measurements require numerical techniques whose quality must be assessed. To this end, and also in view of future monitoring programs or surveys such as the LSST, a blind signal processing competition named Time Delay Challenge 1 (TDC1) was held in 2014. The aim of the present paper, which is based on the simulated light curves from the TDC1, is double. First, we test the performance of the time-delay measurement techniques currently used in COSMOGRAIL. Second, we analyse the quantity and quality of the harvest of time delays obtained from the TDC1 simulations. To achieve these goals, we first discover time delays through a careful inspection of the light curves via a dedicated visual interface. Our measurement algorithms can then be applied to the data in an automated way. We show that our techniques have no significant biases, and yield adequate uncertainty estimates resulting in reduced chi2 values between 0.5 and 1.0. We provide estimates for the number and precision of time-delay measurements that can be expected from future time-delay monitoring campaigns as a function of the photometric signal-to-noise ratio and of the true time delay. We make our blind measurements on the TDC1 data publicly available

preprint2016arXiv

Firedec: a two-channel finite-resolution image deconvolution algorithm

We present a two-channel deconvolution method that decomposes images into a parametric point-source channel and a pixelized extended-source channel. Based on the central idea of the deconvolution algorithm proposed by Magain, Courbin & Sohy (1998), the method aims at improving the resolution of the data rather than at completely removing the point spread function (PSF). Improvements over the original method include a better regularization of the pixel channel of the image, based on wavelet filtering and multiscale analysis, and a better controlled separation of the point source vs. the extended source. In addition, the method is able to simultaneously deconvolve many individual frames of the same object taken with different instruments under different PSF conditions. For this purpose, we introduce a general geometric transformation between individual images. This transformation allows the combination of the images without having to interpolate them. We illustrate the capability of our algorithm using real and simulated images with complex diffraction-limited PSF.

preprint2016arXiv

Multi-band morpho-Spectral Component Analysis Deblending Tool (MuSCADeT): Deblending colourful objects

We introduce a new algorithm for colour separation and deblending of multi-band astronomical images called MuSCADeT which is based on Morpho-spectral Component Analysis of multi-band images. The MuSCADeT algorithm takes advantage of the sparsity of astronomical objects in morphological dictionaries such as wavelets and their differences in spectral energy distribution (SED) across multi-band observations. This allows us to devise a model independent and automated approach to separate objects with different colours. We show with simulations that we are able to separate highly blended objects and that our algorithm is robust against SED variations of objects across the field of view. To confront our algorithm with real data, we use HST images of the strong lensing galaxy cluster MACS J1149+2223 and we show that MuSCADeT performs better than traditional profile-fitting techniques in deblending the foreground lensing galaxies from background lensed galaxies. Although the main driver for our work is the deblending of strong gravitational lenses, our method is fit to be used for any purpose related to deblending of objects in astronomical images. An example of such an application is the separation of the red and blue stellar populations of a spiral galaxy in the galaxy cluster Abell 2744. We provide a python package along with all simulations and routines used in this paper to contribute to reproducible research efforts. Codes can be found at http://lastro.epfl.ch/page-126973.html.

preprint2016arXiv

Stellar classification from single-band imaging using machine learning

Information on the spectral types of stars is of great interest in view of the exploitation of space-based imaging surveys. In this article, we investigate the classification of stars into spectral types using only the shape of their diffraction pattern in a single broad-band image. We propose a supervised machine learning approach to this endeavour, based on principal component analysis (PCA) for dimensionality reduction, followed by artificial neural networks (ANNs) estimating the spectral type. Our analysis is performed with image simulations mimicking the Hubble Space Telescope (HST) Advanced Camera for Surveys (ACS) in the F606W and F814W bands, as well as the Euclid VIS imager. We first demonstrate this classification in a simple context, assuming perfect knowledge of the point spread function (PSF) model and the possibility of accurately generating mock training data for the machine learning. We then analyse its performance in a fully data-driven situation, in which the training would be performed with a limited subset of bright stars from a survey, and an unknown PSF with spatial variations across the detector. We use simulations of main-sequence stars with flat distributions in spectral type and in signal-to-noise ratio, and classify these stars into 13 spectral subclasses, from O5 to M5. Under these conditions, the algorithm achieves a high success rate both for Euclid and HST images, with typical errors of half a spectral class. Although more detailed simulations would be needed to assess the performance of the algorithm on a specific survey, this shows that stellar classification from single-band images is well possible.

preprint2016arXiv

The PCA Lens-Finder: application to CFHTLS

We present the results of a new search for galaxy-scale strong lensing systems in CFHTLS Wide. Our lens-finding technique involves a preselection of potential lens galaxies, applying simple cuts in size and magnitude. We then perform a Principal Component Analysis of the galaxy images, ensuring a clean removal of the light profile. Lensed features are searched for in the residual images using the clustering topometric algorithm DBSCAN. We find 1098 lens candidates that we inspect visually, leading to a cleaned sample of 109 new lens candidates. Using realistic image simulations we estimate the completeness of our sample and show that it is independent of source surface brightness, Einstein ring size (image separation) or lens redshift. We compare the properties of our sample to previous lens searches in CFHTLS. Including the present search, the total number of lenses found in CFHTLS amounts to 678, which corresponds to ~4 lenses per square degree down to i=24.8. This is equivalent to ~ 60.000 lenses in total in a survey as wide as Euclid, but at the CFHTLS resolution and depth.

preprint2015arXiv

A Consistent Picture Emerges: A Compact X-ray Continuum Emission Region in the Gravitationally Lensed Quasar SDSS J0924+0219

We analyze the optical, UV, and X-ray microlensing variability of the lensed quasar SDSS J0924+0219 using six epochs of Chandra data in two energy bands (spanning 0.4-8.0 keV, or 1-20 keV in the quasar rest frame), 10 epochs of F275W (rest-frame 1089A) Hubble Space Telescope data, and high-cadence R-band (rest-frame 2770A) monitoring spanning eleven years. Our joint analysis provides robust constraints on the extent of the X-ray continuum emission region and the projected area of the accretion disk. The best-fit half-light radius of the soft X-ray continuum emission region is between 5x10^13 and 10^15 cm, and we find an upper limit of 10^15 cm for the hard X-rays. The best-fit soft-band size is about 13 times smaller than the optical size, and roughly 7 GM_BH/c^2 for a 2.8x10^8 M_sol black hole, similar to the results for other systems. We find that the UV emitting region falls in between the optical and X-ray emitting regions at 10^14 cm < r_1/2,UV < 3x10^15 cm. Finally, the optical size is significantly larger, by 1.5*sigma, than the theoretical thin-disk estimate based on the observed, magnification-corrected I-band flux, suggesting a shallower temperature profile than expected for a standard disk.

preprint2015arXiv

Evaluating the effect of stellar multiplicity on the PSF of space-based weak lensing surveys

The next generation of space-based telescopes used for weak lensing surveys will require exquisite point spread function (PSF) determination. Previously negligible effects may become important in the reconstruction of the PSF, in part because of the improved spatial resolution. In this paper, we show that unresolved multiple star systems can affect the ellipticity and size of the PSF and that this effect is not cancelled even when using many stars in the reconstruction process. We estimate the error in the reconstruction of the PSF due to the binaries in the star sample both analytically and with image simulations for different PSFs and stellar populations. The simulations support our analytical finding that the error on the size of the PSF is a function of the multiple stars distribution and of the intrinsic value of the size of the PSF, i.e. if all stars were single. Similarly, the modification of each of the complex ellipticity components (e1,e2) depends on the distribution of multiple stars and on the intrinsic complex ellipticity. Using image simulations, we also show that the predicted error in the PSF shape is a theoretical limit that can be reached only if large number of stars (up to thousands) are used together to build the PSF at any desired spatial position. For a lower number of stars, the PSF reconstruction is worse. Finally, we compute the effect of binarity for different stellar magnitudes and show that bright stars alter the PSF size and ellipticity more than faint stars. This may affect the design of PSF calibration strategies and the choice of the related calibration fields.

preprint2014arXiv

A PCA-based automated finder for galaxy-scale strong lenses

We present an algorithm using Principal Component Analysis (PCA) to subtract galaxies from imaging data, and also two algorithms to find strong, galaxy-scale gravitational lenses in the resulting residual image. The combined method is optimized to find full or partial Einstein rings. Starting from a pre-selection of potential massive galaxies, we first perform a PCA to build a set of basis vectors. The galaxy images are reconstructed using the PCA basis and subtracted from the data. We then filter the residual image with two different methods. The first uses a curvelet (curved wavelets) filter of the residual images to enhance any curved/ring feature. The resulting image is transformed in polar coordinates, centered on the lens galaxy center. In these coordinates, a ring is turned into a line, allowing us to detect very faint rings by taking advantage of the integrated signal-to-noise in the ring (a line in polar coordinates). The second way of analysing the PCA-subtracted images identifies structures in the residual images and assesses whether they are lensed images according to their orientation, multiplicity and elongation. We apply the two methods to a sample of simulated Einstein rings, as they would be observed with the ESA Euclid satellite in the VIS band. The polar coordinates transform allows us to reach a completeness of 90% and a purity of 86%, as soon as the signal-to-noise integrated in the ring is higher than 30, and almost independent of the size of the Einstein ring. Finally, we show with real data that our PCA-based galaxy subtraction scheme performs better than traditional subtraction based on model fitting to the data. Our algorithm can be developed and improved further using machine learning and dictionary learning methods, which would extend the capabilities of the method to more complex and diverse galaxy shapes.

preprint2014arXiv

Cosmology from gravitational lens time delays and Planck data

Under the assumption of a flat Lambda-CDM cosmology, recent data from the Planck satellite point toward a Hubble constant that is in tension with that measured by gravitational lens time delays and by the local distance ladder. Prosaically, this difference could arise from unknown systematic uncertainties in some of the measurements. More interestingly -- if systematics were ruled out -- resolving the tension would require a departure from the flat Lambda-CDM cosmology, introducing for example a modest amount of spatial curvature, or a non-trivial dark energy equation of state. To begin to address these issues, we present here an analysis of the gravitational lens RXJ1131-1231 that is improved in one particular regard: we examine the issue of systematic error introduced by an assumed lens model density profile. We use more flexible gravitational lens models with baryonic and dark matter components, and find that the exquisite Hubble Space Telescope image with thousands of intensity pixels in the Einstein ring and the stellar velocity dispersion of the lens contain sufficient information to constrain these more flexible models. The total uncertainty on the time-delay distance is 6.6% for a single system. We proceed to combine our improved time-delay distance measurements with the WMAP9 and Planck posteriors. In an open Lambda-CDM model, the data for RXJ1131-1231 in combination with Planck favor a flat universe with Omega_k = 0.00+0.01/-0.02 (68% CI). In a flat wCDM model, the combination of RXJ1131-1231 and Planck yields w = -1.52+0.19/-0.20 (68% CI).

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XI. Techniques for time delay measurement in presence of microlensing

Measuring time delays between the multiple images of gravitationally lensed quasars is now recognized as a competitive way to constrain the cosmological parameters, and it is complementary with other cosmological probes. This requires long and well sampled optical light curves of numerous lensed quasars, such as those obtained by the COSMOGRAIL collaboration. High-quality data from our monitoring campaign call for novel numerical techniques to robustly measure the delays, as well as the associated random and systematic uncertainties, even in the presence of microlensing variations. We propose three different point estimators to measure time delays, which are explicitly designed to handle light curves with extrinsic variability. These methods share a common formalism, which enables them to process data from n-image lenses. Since the estimators rely on significantly contrasting ideas, we expect them to be sensitive to different bias sources. For each method and data set, we empirically estimate both the precision and accuracy (bias) of the time delay measurement using simulated light curves with known time delays that closely mimic the observations. Finally, we test the self-consistency of our approach, and we demonstrate that our bias estimation is serviceable. These new methods, including the empirical uncertainty estimator, will represent the standard benchmark for analyzing the COSMOGRAIL light curves.

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XII. Time delays of the doubly lensed quasars SDSS J1206+4332 and HS 2209+1914

Aims. Within the framework of the COSMOGRAIL collaboration we present 7- and 8.5-year-long light curves and time-delay estimates for two gravitationally lensed quasars: SDSS J1206+4332 and HS 2209+1914. Methods. We monitored these doubly lensed quasars in the R-band using four telescopes: the Mercator, Maidanak, Himalayan Chandra, and Euler Telescopes, together spanning a period of 7 to 8.5 observing seasons from mid-2004 to mid-2011. The photometry of the quasar images was obtained through simultaneous deconvolution of these data. The time delays were determined from these resulting light curves using four very different techniques: a dispersion method, a spline fit, a regression difference technique, and a numerical model fit. This minimizes the bias that might be introduced by the use of a single method. Results. The time delay for SDSS J1206+4332 is Delta_t AB = 111.3 +/- 3 days with A leading B, confirming a previously published result within the error bars. For HS 2209+1914 we present a new time delay of Delta_t BA = 20.0 +/- 5 days with B leading A. Conclusions. The combination of data from up to four telescopes have led to well-sampled and nearly 9-season-long light curves, which were necessary to obtain these results, especially for the compact doubly lensed quasar HS 2209+1914.

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIII: Time delays and 9-yr optical monitoring of the lensed quasar RX J1131-1231

We present the results from nine years of optically monitoring the gravitationally lensed z=0.658 quasar RX J1131-1231. The R-band light curves of the four individual images of the quasar were obtained using deconvolution photometry for a total of 707 epochs. Several sharp quasar variability features strongly constrain the time delays between the quasar images. Using three different numerical techniques, we measure these delays for all possible pairs of quasar images while always processing the four light curves simultaneously. For all three methods, the delays between the three close images A, B, and C are compatible with being 0, while we measure the delay of image D to be 91 days, with a fractional uncertainty of 1.5% (1 sigma), including systematic errors. Our analysis of random and systematic errors accounts in a realistic way for the observed quasar variability, fluctuating microlensing magnification over a broad range of temporal scales, noise properties, and seasonal gaps. Finally, we find that our time-delay measurement methods yield compatible results when applied to subsets of the data.

preprint2013arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIV. Time delay of the doubly lensed quasar SDSS J1001+5027

This paper presents optical R-band light curves and the time delay of the doubly imaged gravitationally lensed quasar SDSS J1001+5027 at a redshift of 1.838. We have observed this target for more than six years, between March 2005 and July 2011, using the 1.2-m Mercator Telescope, the 1.5-m telescope of the Maidanak Observatory, and the 2-m Himalayan Chandra Telescope. Our resulting light curves are composed of 443 independent epochs, and show strong intrinsic quasar variability, with an amplitude of the order of 0.2 magnitudes. From this data, we measure the time delay using five different methods, all relying on distinct approaches. One of these techniques is a new development presented in this paper. All our time-delay measurements are perfectly compatible. By combining them, we conclude that image A is leading B by 119.3 +/- 3.3 days (1 sigma, 2.8% uncertainty), including systematic errors. It has been shown recently that such accurate time-delay measurements offer a highly complementary probe of dark energy and spatial curvature, as they independently constrain the Hubble constant. The next mandatory step towards using SDSS J1001+5027 in this context will be the measurement of the velocity dispersion of the lensing galaxy, in combination with deep Hubble Space Telescope imaging.

preprint2013arXiv

Dark energy with gravitational lens time delays

Strong lensing gravitational time delays are a powerful and cost effective probe of dark energy. Recent studies have shown that a single lens can provide a distance measurement with 6-7 % accuracy (including random and systematic uncertainties), provided sufficient data are available to determine the time delay and reconstruct the gravitational potential of the deflector. Gravitational-time delays are a low redshift (z~0-2) probe and thus allow one to break degeneracies in the interpretation of data from higher-redshift probes like the cosmic microwave background in terms of the dark energy equation of state. Current studies are limited by the size of the sample of known lensed quasars, but this situation is about to change. Even in this decade, wide field imaging surveys are likely to discover thousands of lensed quasars, enabling the targeted study of ~100 of these systems and resulting in substantial gains in the dark energy figure of merit. In the next decade, a further order of magnitude improvement will be possible with the 10000 systems expected to be detected and measured with LSST and Euclid. To fully exploit these gains, we identify three priorities. First, support for the development of software required for the analysis of the data. Second, in this decade, small robotic telescopes (1-4m in diameter) dedicated to monitoring of lensed quasars will transform the field by delivering accurate time delays for ~100 systems. Third, in the 2020's, LSST will deliver 1000's of time delays; the bottleneck will instead be the aquisition and analysis of high resolution imaging follow-up. Thus, the top priority for the next decade is to support fast high resolution imaging capabilities, such as those enabled by the James Webb Space Telescope and next generation adaptive optics systems on large ground based telescopes.

preprint2013arXiv

Image Analysis for Cosmology: Results from the GREAT10 Galaxy Challenge

In this paper we present results from the weak lensing shape measurement GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Galaxy Challenge. This marks an order of magnitude step change in the level of scrutiny employed in weak lensing shape measurement analysis. We provide descriptions of each method tested and include 10 evaluation metrics over 24 simulation branches. GREAT10 was the first shape measurement challenge to include variable fields; both the shear field and the Point Spread Function (PSF) vary across the images in a realistic manner. The variable fields enable a variety of metrics that are inaccessible to constant shear simulations including a direct measure of the impact of shape measurement inaccuracies, and the impact of PSF size and ellipticity, on the shear power spectrum. To assess the impact of shape measurement bias for cosmic shear we present a general pseudo-Cl formalism, that propagates spatially varying systematics in cosmic shear through to power spectrum estimates. We also show how one-point estimators of bias can be extracted from variable shear simulations. The GREAT10 Galaxy Challenge received 95 submissions and saw a factor of 3 improvement in the accuracy achieved by shape measurement methods. The best methods achieve sub-percent average biases. We find a strong dependence in accuracy as a function of signal-to-noise, and indications of a weak dependence on galaxy type and size. Some requirements for the most ambitious cosmic shear experiments are met above a signal-to-noise ratio of 20. These results have the caveat that the simulated PSF was a ground-based PSF. Our results are a snapshot of the accuracy of current shape measurement methods and are a benchmark upon which improvement can continue. This provides a foundation for a better understanding of the strengths and limitations of shape measurement methods.

preprint2013arXiv

Probing the inner structure of distant AGNs with gravitational lensing

Microlensing is a powerful technique which can be used to study the continuum and the broad line emitting regions in distant AGNs. After a brief description of the methods and required data, we present recent applications of this technique. We show that microlensing allows one to measure the temperature profile of the accretion disc, estimate the size and study the geometry of the region emitting the broad emission lines.

preprint2013arXiv

Two accurate time-delay distances from strong lensing: Implications for cosmology

Strong gravitational lenses with measured time delays between the multiple images and models of the lens mass distribution allow a one-step determination of the time-delay distance, and thus a measure of cosmological parameters. We present a blind analysis of the gravitational lens RXJ1131-1231 incorporating (1) the newly measured time delays from COSMOGRAIL, (2) archival HST imaging of the lens system, (3) a new velocity-dispersion measurement of the lens galaxy of 323+/-20km/s based on Keck spectroscopy, and (4) a characterization of the line-of-sight structures via observations of the lens' environment and ray tracing through the Millennium Simulation. Our blind analysis is designed to prevent experimenter bias. The joint analysis of the data sets allows a time-delay distance measurement to 6% precision that takes into account all known systematic uncertainties. In combination with the WMAP7 data set in flat wCDM cosmology, our unblinded cosmological constraints for RXJ1131-1231 are: H_0=80.0+5.8/-5.7km/s/Mpc, OmegaDE=0.79+/-0.03 and w=-1.25+0.17/-0.21. We find the results to be statistically consistent with those from the analysis of the gravitational lens B1608+656. The joint constraints from the two lenses and WMAP7 are H_0=75.2+4.4/-4.2km/s/Mpc, OmegaDE=0.76+0.02/-0.03 and w=-1.14+0.17/-0.20 in flat wCDM, and H_0=73.1+2.4/-3.6km/s/Mpc, OmegaL=0.75+0.01/-0.02 and OmegaK=0.003+0.005/-0.006 in open LCDM. Time-delay lenses constrain especially tightly the Hubble constant (5.7% and 4.0% respectively in wCDM and open LCDM) and curvature of the universe. They complement well other cosmological probes, and provide an independent check of unknown systematics. Our measurement of the Hubble constant is completely independent of those based on the local distance ladder method, providing an important consistency check of the standard cosmological model and of general relativity.

preprint2012arXiv

A fast empirical method for galaxy shape measurements in weak lensing surveys

We describe a simple and fast method to correct ellipticity measurements of galaxies from the distortion by the instrumental and atmospheric point spread function (PSF), in view of weak lensing shear measurements. The method performs a classification of galaxies and associated PSFs according to measured shape parameters, and corrects the measured galaxy ellipticites by querying a large lookup table (LUT), built by supervised learning. We have applied this new method to the GREAT10 image analysis challenge, and present in this paper a refined solution that obtains the competitive quality factor of Q = 104, without any shear power spectrum denoising or training. Of particular interest is the efficiency of the method, with a processing time below 3 ms per galaxy on an ordinary CPU.

preprint2012arXiv

Image Analysis for Cosmology: Results from the GREAT10 Star Challenge

We present the results from the first public blind PSF reconstruction challenge, the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) Star Challenge. Reconstruction of a spatially varying PSF, sparsely sampled by stars, at non-star positions is a critical part in the image analysis for weak lensing where inaccuracies in the modelled ellipticity and size-squared can impact the ability to measure the shapes of galaxies. This is of importance because weak lensing is a particularly sensitive probe of dark energy, and can be used to map the mass distribution of large scale structure. Participants in the challenge were presented with 27,500 stars over 1300 images subdivided into 26 sets, where in each set a category change was made in the type or spatial variation of the PSF. Thirty submissions were made by 9 teams. The best methods reconstructed the PSF with an accuracy of ~0.00025 in ellipticity and ~0.00074 in size squared. For a fixed pixel scale narrower PSFs were found to be more difficult to model than larger PSFs, and the PSF reconstruction was severely degraded with the inclusion of an atmospheric turbulence model (although this result is likely to be a strong function of the amplitude of the turbulence power spectrum).

preprint2012arXiv

Interpolating point spread function anisotropy

Planned wide-field weak lensing surveys are expected to reduce the statistical errors on the shear field to unprecedented levels. In contrast, systematic errors like those induced by the convolution with the point spread function (PSF) will not benefit from that scaling effect and will require very accurate modeling and correction. While numerous methods have been devised to carry out the PSF correction itself, modeling of the PSF shape and its spatial variations across the instrument field of view has, so far, attracted much less attention. This step is nevertheless crucial because the PSF is only known at star positions while the correction has to be performed at any position on the sky. A reliable interpolation scheme is therefore mandatory and a popular approach has been to use low-order bivariate polynomials. In the present paper, we evaluate four other classical spatial interpolation methods based on splines (B-splines), inverse distance weighting (IDW), radial basis functions (RBF) and ordinary Kriging (OK). These methods are tested on the Star-challenge part of the GRavitational lEnsing Accuracy Testing 2010 (GREAT10) simulated data and are compared with the classical polynomial fitting (Polyfit). We also test all our interpolation methods independently of the way the PSF is modeled, by interpolating the GREAT10 star fields themselves (i.e., the PSF parameters are known exactly at star positions). We find in that case RBF to be the clear winner, closely followed by the other local methods, IDW and OK. The global methods, Polyfit and B-splines, are largely behind, especially in fields with (ground-based) turbulent PSFs. In fields with non-turbulent PSFs, all interpolators reach a variance on PSF systematics $σ_{sys}^2$ better than the $1\times10^{-7}$ upper bound expected by future space-based surveys, with the local interpolators performing better than the global ones.

preprint2012arXiv

Microlensing of the broad line region in 17 lensed quasars

When an image of a strongly lensed quasar is microlensed, the different components of its spectrum are expected to be differentially magnified owing to the different sizes of the corresponding emitting region. Chromatic changes are expected to be observed in the continuum while the emission lines should be deformed as a function of the size, geometry and kinematics of the regions from which they originate. Microlensing of the emission lines has been reported only in a handful of systems so far. In this paper we search for microlensing deformations of the optical spectra of pairs of images in 17 lensed quasars. This sample is composed of 13 pairs of previously unpublished spectra and four pairs of spectra from literature. Our analysis is based on a spectral decomposition technique which allows us to isolate the microlensed fraction of the flux independently of a detailed modeling of the quasar emission lines. Using this technique, we detect microlensing of the continuum in 85% of the systems. Among them, 80% show microlensing of the broad emission lines. Focusing on the most common lines in our spectra (CIII] and MgII) we detect microlensing of either the blue or the red wing, or of both wings with the same amplitude. This observation implies that the broad line region is not in general spherically symmetric. In addition, the frequent detection of microlensing of the blue and red wings independently but not simultaneously with a different amplitude, does not support existing microlensing simulations of a biconical outflow. Our analysis also provides the intrinsic flux ratio between the lensed images and the magnitude of the microlensing affecting the continuum. These two quantities are particularly relevant for the determination of the fraction of matter in clumpy form in galaxies and for the detection of dark matter substructures via the identification of flux ratio anomalies.

preprint2012arXiv

Spectroscopy of extended Lyα envelopes around z=4.5 quasars

What are the frequency, shape, kinematics, and luminosity of Lyα envelopes surrounding radio-quiet quasars at high redshift, and is the luminosity of these envelopes related to that of the quasar or not? As a first step towards answering these questions, we have searched for Lyα envelopes around six radio-quiet quasars at z~4.5, using deep spectra taken with the FORS2 spectrograph attached to the UT1 of the Very Large Telescope (VLT). Using the multi-slit mode allows us to observe several point spread function stars simultaneously with the quasar, and to remove the point-like emission from the quasar, unveiling the faint underlying Lyα envelope with unprecedented depth. An envelope is detected around four of the six quasars, which suggests that these envelopes are very frequent. Their diameter varies in the range 26<d<64 kpc, their surface brightness in the range 3x10^{-19}<μ<2x10^{-17} erg/s/cm^2/arcsec^2, and their luminosity in the range 10^{42}<L(Lyα)<10^{44} erg/s. Their shape may be strongly asymmetric. The Lyα emission line full width at half maximum (FWHM) is 900<FWHM<2200 km/s and its luminosity correlates with that of the broad line region (BLR) of the quasar, with the notable exception of BR2237-0607, the brightest object in our sample. The same holds for the relation between the envelope Lyα luminosity and the ionizing luminosity of the quasar. While the deep slit spectroscopy presented in this paper is very efficient at detecting very faint Lyα envelopes, narrow-band imaging is now needed to measure accurately their spatial extent, radial luminosity profile, and total luminosity. These observables are crucial to help us discriminate between the three possible radiation processes responsible for the envelope emission: (i) cold accretion, (ii) fluorescence induced by the quasar, and (iii) scattering of the BLR photons by cool gas.

preprint2012arXiv

The Hubble constant and new discoveries in cosmology

We report the outcome of a 3-day workshop on the Hubble constant (H_0) that took place during February 6-8 2012 at the Kavli Institute for Particle Astrophysics and Cosmology, on the campus of Stanford University. The participants met to address the following questions. Are there compelling scientific reasons to obtain more precise and more accurate measurements of H_0 than currently available? If there are, how can we achieve this goal? The answers that emerged from the workshop are (1) better measurements of H_0 provide critical independent constraints on dark energy, spatial curvature of the Universe, neutrino physics, and validity of general relativity, (2) a measurement of H_0 to 1% in both precision and accuracy, supported by rigorous error budgets, is within reach for several methods, and (3) multiple paths to independent determinations of H_0 are needed in order to access and control systematics.

preprint2011arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses IX. Time delays, lens dynamics and baryonic fraction in HE 0435-1223

We present accurate time delays for the quadruply imaged quasar HE 0435-1223. The delays were measured from 575 independent photometric points obtained in the R-band between January 2004 and March 2010. With seven years of data, we clearly show that quasar image A is affected by strong microlensing variations and that the time delays are best expressed relative to quasar image B. We measured Delta_t(BC) = 7.8+/-0.8 days, Delta_t(BD) = -6.5+/-0.7 days and Delta_t_CD = -14.3+/-0.8 days. We spacially deconvolved HST NICMOS2 F160W images to derive accurate astrometry of the quasar images and to infer the light profile of the lensing galaxy. We combined these images with a stellar population fitting of a deep VLT spectrum of the lensing galaxy to estimate the baryonic fraction, $f_b$, in the Einstein radius. We measured f_b = 0.65+0.13-0.10 if the lensing galaxy has a Salpeter IMF and f_b = 0.45+0.04-0.07 if it has a Kroupa IMF. The spectrum also allowed us to estimate the velocity dispersion of the lensing galaxy, sigma_ap = 222+/-34 km/s. We used f_b and sigma_ap to constrain an analytical model of the lensing galaxy composed of an Hernquist plus generalized NFW profile. We solve the Jeans equations numerically for the model and explored the parameter space under the additional requirement that the model must predict the correct astrometry for the quasar images. Given the current error bars on f_b and sigma_ap, we did not constrain H0 yet with high accuracy, i.e., we found a broad range of models with chi^2 < 1. However, narrowing this range is possible, provided a better velocity dispersion measurement becomes available. In addition, increasing the depth of the current HST imaging data of HE 0435-1223 will allow us to combine our constraints with lens reconstruction techniques that make use of the full Einstein ring that is visible in this object.

preprint2011arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses X. Modeling based on high-precision astrometry of a sample of 25 lensed quasars: consequences for ellipticity, shear, and astrometric anomalies

(abridged) Gravitationally lensed quasars can be used as powerful cosmological and astrophysical probes. We can (i) infer the Hubble constant based on the time-delay technique, (ii) unveil substructures along the l.o.s. toward distant galaxies, and (iii) compare the shape and the slope of baryons and dark matter distributions in galaxies. To reach these goals, we need high-accuracy astrometry and morphology measurements of the lens. In this work, we first present new astrometry for 11 lenses with measured time delays. Using MCS deconvolution on NIC2 HST images, we reached an astrometric accuracy of about 1-2.5 mas and an accurate shape measurement of the lens galaxy. Second, we combined these measurements with those of 14 other systems to present new mass models of these lenses. This led to the following results: 1) In 4 double-image quasars, we show that the influence of the lens environment on the time delay can easily be quantified and modeled, hence putting these lenses with high priority for time-delay determination. 2) For quadruple-image quasars, the difficulty often encountered in reproducing the image positions to milli-arcsec accuracy (astrometric anomaly) is overcome by explicitly including the nearest visible galaxy in the model. However, one anomalous system (J1131-1231) does not show any luminous perturber in its vicinity, and three others (WFI2026-4536, WFI2033-4723, and B2045+265) have problematic modeling. These 4 systems are the best candidates for a pertubation by a dark matter substructure. 3) We find a significant correlation between the PA of the light and of the mass distributions in lensing galaxies. In contrast with other studies, we find that the ellipticity of the light and of the mass also correlate well, suggesting that the overall spatial distribution of matter is not very different from the baryon distribution in the inner \sim 5 kpc of lensing galaxies.

preprint2011arXiv

Euclid Definition Study Report

Euclid is a space-based survey mission from the European Space Agency designed to understand the origin of the Universe's accelerating expansion. It will use cosmological probes to investigate the nature of dark energy, dark matter and gravity by tracking their observational signatures on the geometry of the universe and on the cosmic history of structure formation. The mission is optimised for two independent primary cosmological probes: Weak gravitational Lensing (WL) and Baryonic Acoustic Oscillations (BAO). The Euclid payload consists of a 1.2 m Korsch telescope designed to provide a large field of view. It carries two instruments with a common field-of-view of ~0.54 deg2: the visual imager (VIS) and the near infrared instrument (NISP) which contains a slitless spectrometer and a three bands photometer. The Euclid wide survey will cover 15,000 deg2 of the extragalactic sky and is complemented by two 20 deg2 deep fields. For WL, Euclid measures the shapes of 30-40 resolved galaxies per arcmin2 in one broad visible R+I+Z band (550-920 nm). The photometric redshifts for these galaxies reach a precision of dz/(1+z) < 0.05. They are derived from three additional Euclid NIR bands (Y, J, H in the range 0.92-2.0 micron), complemented by ground based photometry in visible bands derived from public data or through engaged collaborations. The BAO are determined from a spectroscopic survey with a redshift accuracy dz/(1+z) =0.001. The slitless spectrometer, with spectral resolution ~250, predominantly detects Ha emission line galaxies. Euclid is a Medium Class mission of the ESA Cosmic Vision 2015-2025 programme, with a foreseen launch date in 2019. This report (also known as the Euclid Red Book) describes the outcome of the Phase A study.

preprint2011arXiv

Evolution of the observed Ly-alpha luminosity function from z = 6.5 to z = 7.7: evidence for the epoch of reionization ?

Aims. Ly-alpha emitters (LAEs) can be detected out to very high redshifts during the epoch of reionization. The evolution of the LAE luminosity function with redshift is a direct probe of the Ly-alpha transmission of the intergalactic medium (IGM), and therefore of the IGM neutral-hydrogen fraction. Measuring the Ly-alpha luminosity function (LF) of LAEs at redshift z = 7.7 therefore allows us to constrain the ionizing state of the Universe at this redshift. Methods. We observed three 7.5'x7.5' fields with the HAWK-I instrument at the VLT with a narrow band filter centred at 1.06 $μ$m and targeting LAEs at redshift z ~ 7.7. The fields were chosen for the availability of multiwavelength data. One field is a galaxy cluster, the Bullet Cluster, which allowed us to use gravitational amplification to probe luminosities that are fainter than in the field. The two other fields are subareas of the GOODS Chandra Deep Field South and CFHTLS-D4 deep field. We selected z=7.7 LAE candidates from a variety of colour criteria, in particular from the absence of detection in the optical bands. Results. We do not find any LAE candidates at z = 7.7 in ~2.4 x 10^4 Mpc^3 down to a narrow band AB magnitude of ~ 26, which allows us to infer robust constraints on the Ly-alpha LAE luminosity function at this redshift. Conclusions. The predicted mean number of objects at z = 6.5, derived from somewhat different LFs of Hu et al. (2010), Ouchi et al. (2010), and Kashikawa et al. (2011) are 2.5, 13.7, and 11.6, respectively. Depending on which of these LFs we refer to, we exclude a scenario with no evolution from z = 6.5 to z = 7.7 at 85% confidence without requiring a strong change in the IGM Ly-alpha transmission, or at 99% confidence with a significant quenching of the IGM Ly-alpha transmission, possibly from a strong increase in the high neutral-hydrogen fraction between these two redshifts.

preprint2011arXiv

On the effect of image denoising on galaxy shape measurements

Weak gravitational lensing is a very sensitive way of measuring cosmological parameters, including dark energy, and of testing current theories of gravitation. In practice, this requires exquisite measurement of the shapes of billions of galaxies over large areas of the sky, as may be obtained with the EUCLID and WFIRST satellites. For a given survey depth, applying image denoising to the data both improves the accuracy of the shape measurements and increases the number density of galaxies with a measurable shape. We perform simple tests of three different denoising techniques, using synthetic data. We propose a new and simple denoising method, based on wavelet decomposition of the data and a Wiener filtering of the resulting wavelet coefficients. When applied to the GREAT08 challenge dataset, this technique allows us to improve the quality factor of the measurement (Q; GREAT08 definition), by up to a factor of two. We demonstrate that the typical pixel size of the EUCLID optical channel will allow us to use image denoising.

preprint2011arXiv

Three QSOs acting as strong gravitational lenses

We report the discovery of three new cases of QSOs acting as strong gravitational lenses on background emission line galaxies: SDSS J0827+5224 (zQSO = 0.293, zs = 0.412), SDSS J0919+2720 (zQSO = 0.209, zs = 0.558), SDSS J1005+4016 (zQSO = 0.230, zs = 0.441). The selection was carried out using a sample of 22,298 SDSS spectra displaying at least four emission lines at a redshift beyond that of the foreground QSO. The lensing nature is confirmed from Keck imaging and spectroscopy, as well as from HST/WFC3 imaging in the F475W and F814W filters. Two of the QSOs have face-on spiral host galaxies and the third is a QSO+galaxy pair. The velocity dispersion of the host galaxies, inferred from simple lens modeling, is between σ_v = 210 and 285 km/s, making these host galaxies comparable in mass with the SLACS sample of early-type strong lenses.

preprint2011arXiv

VLT adaptive optics search for luminous substructures in the lens galaxy towards SDSS J0924+0219

Anomalous flux ratios between quasar images are suspected to be caused by substructures in lens galaxies. We present new deep and high resolution H and Ks imaging of the strongly lensed quasar SDSS J0924+0219 obtained using the ESO VLT with adaptive optics and the Laser Guide Star system. SDSS J0924+0219 is particularly interesting as the observed flux ratio between the quasar images vastly disagree with the predictions from smooth mass models. With our adaptive optics observations we find a luminous object, Object L, located ~0.3" to the North of the lens galaxy, but we show that it can not be responsible for the anomalous flux ratios. Object L as well as a luminous extension of the lens galaxy to the South are seen in the archival HST/ACS image in the F814W filter. This suggests that Object L is part of a bar in the lens galaxy, as also supported by the presence of a significant disk component in the light profile of the lens galaxy. Finally, we do not find evidence for any other luminous substructure that may explain the quasar images flux ratios. However, owe to the persistence of the flux ratio anomaly over time (~ 7 years) a combination of microlensing and milli-lensing is the favorite explanation for the observations.

preprint2011arXiv

Zooming into the broad line region of the gravitationally lensed quasar Q2237+0305 = the Einstein Cross: III. Determination of the size and structure of the CIV and CIII] emitting regions using microlensing

We aim to use microlensing taking place in the lensed quasar Q2237+0305 to study the structure of the broad line region and measure the size of the region emitting the CIV and CIII] lines. Methods: Based on 39 spectrophotometric monitoring data points obtained between Oct. 2004 and Dec. 2007, we derived lightcurves for the CIV and CIII] emission lines. We used three different techniques to analyse the microlensing signal. Different components of the lines (narrow, broad and very broad) are identified and studied. We built a library of simulated microlensing lightcurves that reproduce the signal observed in the continuum and in the lines provided only the source size is changed. A Bayesian analysis scheme is then developed to derive the size of the various components of the BLR. Results: 1. The half-light radius of the region emitting the CIV line is found to be R_CIV ~ 66^{+110}_{-46} lt-days = 0.06$^{+0.09}_{-0.04}$ pc = 1.7$^{+2.8}_{-1.1}$\,10$^{17}$ cm (at 68.3% CI). Similar values are obtained for CIII]. Relative sizes of the carbon-line and V-band continuum emitting-regions are also derived with median values of R(line)/R(cont) in the range 4 to 29, depending of the FWHM of the line component. 2. The size of the CIV emitting region agrees with the Radius-Luminosity relationship derived from reverberation mapping. Using the virial theorem we derive the mass of the black hole in Q2237+0305 to be M_BH ~ 10^{8.3+/-0.3} M_sun. 3. We find that the CIV and CIII] lines are produced in at least 2 spatially distinct regions, the most compact one giving rise to the broadest component of the line. The broad and narrow line profiles are slightly different for CIV and CIII]. 4. Our analysis suggests a different structure for the CIV and FeII+III emitting regions, with the latter produced in the inner part of the BLR or in a less extended emitting region than CIV.

preprint2010arXiv

Direct detection of galaxy stellar halos : NGC 3957 as a test case

We present a direct detection of the stellar halo of the edge-on S0 galaxy NGC 3957, using ultra-deep VLT/VIMOS V and R images. This is achieved with a sky subtraction strategy based on infrared techniques. These observations allow us to reach unprecedented high signal-to-noise ratios up to 15 kpc away from the galaxy center, rendering photon-noise negligible. The 1 sigma detection limits are R = 30.6 mag/arcsec^2 and V = 31.4 mag/arcsec^2. We conduct a thorough analysis of the possible sources of systematic errors that could affect the data: flat-fielding, differences in CCD responses, scaling of the sky background, the extended halo itself, and PSF wings. We conclude that the V-R colour of the NGC 3957 halo, calculated between 5 and 8 kpc above the disc plane where the systematic errors are modest, is consistent with an old and preferentially metal-poor normal stellar population, like that revealed in nearby galaxy halos from studies of their resolved stellar content. We do not find support for the extremely red colours found in earlier studies of diffuse halo emission, which we suggest might have been due to residual systematic errors.

preprint2010arXiv

First case of strong gravitational lensing by a QSO : SDSS J0013+1523 at z = 0.120

We present the first case of strong gravitational lensing by a QSO : SDSS J0013+1523, at z = 0.120. The discovery is the result of a systematic search for emission lines redshifted behind QSOs, among 22298 spectra of the SDSS data release 7. Apart from the z = 0.120 spectral features of the foreground QSO, the spectrum of SDSS J0013+1523 also displays the OII and Hbeta emission lines and the OIII doublet, all at the same redshift, z = 0.640. Using sharp Keck adaptive optics K-band images obtained using laser guide stars, we unveil two objects within a radius of 2 arcsec from the QSO. Deep Keck optical spectroscopy clearly confirms one of these objects at z = 0.640 and shows traces of the OIII, emission line of the second object, also at z = 0.640. Lens modeling suggests that they represent two images of the same z = 0.640 emission-line galaxy. Our Keck spectra also allow us to measure the redshift of an intervening galaxy at z = 0.394, located 3.2 arcsec away from the line of sight to the QSO. If the z = 0.120 QSO host galaxy is modeled as a singular isothermal sphere, its mass within the Einstein radius is M_E(r < 1 kpc) = 2.16e10 M_Sun and its velocity dispersion is sigma_SIS = 169 km/s. This is about 1 sigma away from the velocity dispersion estimated from the width of the QSO Hbeta emission line, sigma_*(M_BH) = 124 +/- 47 km/s. Deep optical HST imaging will be necessary to constrain the total radial mass profile of the QSO host galaxy using the detailed shape of the lensed source. This first case of a QSO acting as a strong lens on a more distant object opens new directions in the study of QSO host galaxies.

preprint2008arXiv

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses VII. Time delays and the Hubble constant from WFI J2033-4723

Gravitationally lensed quasars can be used to map the mass distribution in lensing galaxies and to estimate the Hubble constant H0 by measuring the time delays between the quasar images. Here we report the measurement of two independent time delays in the quadruply imaged quasar WFI J2033-4723 (z = 1.66). Our data consist of R-band images obtained with the Swiss 1.2 m EULER telescope located at La Silla and with the 1.3 m SMARTS telescope located at Cerro Tololo. The light curves have 218 independent epochs spanning 3 full years of monitoring between March 2004 and May 2007, with a mean temporal sampling of one observation every 4th day. We measure the time delays using three different techniques, and we obtain Dt(B-A) = 35.5 +- 1.4 days (3.8%) and Dt(B-C) = 62.6 +4.1/-2.3 days (+6.5%/-3.7%), where A is a composite of the close, merging image pair. After correcting for the time delays, we find R-band flux ratios of F_A/F_B = 2.88 +- 0.04, F_A/F_C = 3.38 +- 0.06, and F_A1/F_A2 = 1.37 +- 0.05 with no evidence for microlensing variability over a time scale of three years. However, these flux ratios do not agree with those measured in the quasar emission lines, suggesting that longer term microlensing is present. Our estimate of H0 agrees with the concordance value: non-parametric modeling of the lensing galaxy predicts H0 = 67 +13/-10 km s-1 Mpc-1, while the Single Isothermal Sphere model yields H0 = 63 +7/-3 km s-1 Mpc-1 (68% confidence level). More complex lens models using a composite de Vaucouleurs plus NFW galaxy mass profile show twisting of the mass isocontours in the lensing galaxy, as do the non-parametric models. As all models also require a significant external shear, this suggests that the lens is a member of the group of galaxies seen in field of view of WFI J2033-4723.

preprint2008arXiv

Slit and integral-field optical spectroscopy of the enigmatic quasar HE0450-2958

Interest in the quasar HE0450-2958 arose following the publication of the non-detection of its expected massive host, leading to various interpretations. This article investigates the gaseous and stellar contents of the system through additional VLT/FORS slit spectra and integral field spectroscopy from VLT/VIMOS. We apply our MCS deconvolution algorithm on slit spectra for the separation of the QSO and diffuse components, and develop a new method to remove the point sources in Integral Field Spectra, allowing extraction of velocity maps, narrow-line images, spatially resolved spectra or ionization diagrams of the surroundings of HE0450-2958. The whole system is embedded in gas, mostly ionized by the QSO radiation field and shocks associated with radio jets. The observed gas and star dynamics are unrelated, revealing a strongly perturbed system. Despite longer spectroscopic observations, the host galaxy remains undetected.

preprint1997arXiv

The geometry of the quadruply imaged quasar PG 1115+080; implications for Ho

Time delay measurements have recently been reported for the lensed quasar PG 1115+080. These measurements can be used to derive Ho, but only if we can constrain the lensing potential. We have applied a recently developed deconvolution technique to analyse sub-arcsecond I band images of PG 1115+080, obtained at the Nordic Optical Telescope (NOT) and the Canada France Hawaii Telescope (CFHT). The high performance of the deconvolution code allows us to derive precise positions and magnitudes for the four lensed images of the quasar, as well as for the lensing galaxy. The new measurement of the galaxy position improves its precision by a factor of 3 and thus strengthens the constraints on the lensing potential. With the new data, a range of models incorporating some of the plausible systematic uncertainties yields Ho = 53 (+10/-7) km/s/mpc.

F. Courbin

What is connected

Connect this record

See the researcher in context

Building this map preview

62 published item(s)

Euclid preparation. Calibrated intrinsic galaxy alignments in the Euclid Flagship simulation

Euclid preparation. Galaxy 2-point correlation function modelling in redshift space

Euclid preparation. Testing analytic models of galaxy intrinsic alignments in the Euclid Flagship simulation

Euclid: Early Release Observations -- The extended stellar component of the IC10 dwarf galaxy

Euclid: Galaxy SED reconstruction in the PHZ processing function: impact on the PSF and the role of medium-band filters

Time-delay and lens galaxy redshift in the doubly imaged quasar PS J2305+3714

Euclid preparation. Simulating thousands of Euclid spectroscopic skies

Euclid preparation: XXIII. Derivation of galaxy physical properties with deep machine learning using mock fluxes and H-band images

Euclid preparation. XXX. Performance assessment of the NISP Red-Grism through spectroscopic simulations for the Wide and Deep surveys

Euclid preparation: XIX. Impact of magnification on photometric galaxy clustering

Euclid preparation: XX. The Complete Calibration of the Color-Redshift Relation survey: LBT observations and data release

Euclid preparation. XVIII. The NISP photometric system

Euclid: Forecast constraints on consistency tests of the $Λ$CDM model

Constraining quasar structure using high-frequency microlensing variations and continuum reverberation

Discovery of Strongly Lensed Quasars in the Ultraviolet Near Infrared Optical Northern Survey (UNIONS)

Euclid preparation: I. The Euclid Wide Survey

J1721$+$8842: a gravitationally lensed binary quasar with a proximate damped Lyman-$α$ absorber

COSMOGRAIL XIX: Time delays in 18 strongly lensed quasars from 15 years of optical monitoring

H0LiCOW XI. A weak lensing measurement of the external convergence in the field of the lensed quasar B1608+656 using HST and Subaru deep imaging

HOLISMOKES -- I. Highly Optimised Lensing Investigations of Supernovae, Microlensing Objects, and Kinematics of Ellipticals and Spirals

Is every strong lens model unhappy in its own way? Uniform modelling of a sample of 13 quadruply+ imaged quasars

STRIDES: a 3.9 per cent measurement of the Hubble constant from the strong lens system DES J0408-5354

TDCOSMO. I. An exploration of systematic uncertainties in the inference of $H_0$ from time-delay cosmography

The STRong lensing Insights into the Dark Energy Survey (STRIDES) 2017/2018 follow-up campaign: Discovery of 10 lensed quasars and 10 quasar pairs

Twisted quasar light curves: implications for continuum reverberation mapping of accretion disks

ASTErIsM - Application of topometric clustering algorithms in automatic galaxy detection and classification

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XV. Assessing the achievability and precision of time-delay measurements

Firedec: a two-channel finite-resolution image deconvolution algorithm

Multi-band morpho-Spectral Component Analysis Deblending Tool (MuSCADeT): Deblending colourful objects

Stellar classification from single-band imaging using machine learning

The PCA Lens-Finder: application to CFHTLS

A Consistent Picture Emerges: A Compact X-ray Continuum Emission Region in the Gravitationally Lensed Quasar SDSS J0924+0219

Evaluating the effect of stellar multiplicity on the PSF of space-based weak lensing surveys

A PCA-based automated finder for galaxy-scale strong lenses

Cosmology from gravitational lens time delays and Planck data

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XI. Techniques for time delay measurement in presence of microlensing

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XII. Time delays of the doubly lensed quasars SDSS J1206+4332 and HS 2209+1914

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIII: Time delays and 9-yr optical monitoring of the lensed quasar RX J1131-1231

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses XIV. Time delay of the doubly lensed quasar SDSS J1001+5027

Dark energy with gravitational lens time delays

Image Analysis for Cosmology: Results from the GREAT10 Galaxy Challenge

Probing the inner structure of distant AGNs with gravitational lensing

Two accurate time-delay distances from strong lensing: Implications for cosmology

A fast empirical method for galaxy shape measurements in weak lensing surveys

Image Analysis for Cosmology: Results from the GREAT10 Star Challenge

Interpolating point spread function anisotropy

Microlensing of the broad line region in 17 lensed quasars

Spectroscopy of extended Lyα envelopes around z=4.5 quasars

The Hubble constant and new discoveries in cosmology

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses IX. Time delays, lens dynamics and baryonic fraction in HE 0435-1223

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses X. Modeling based on high-precision astrometry of a sample of 25 lensed quasars: consequences for ellipticity, shear, and astrometric anomalies

Euclid Definition Study Report

Evolution of the observed Ly-alpha luminosity function from z = 6.5 to z = 7.7: evidence for the epoch of reionization ?

On the effect of image denoising on galaxy shape measurements

Three QSOs acting as strong gravitational lenses

VLT adaptive optics search for luminous substructures in the lens galaxy towards SDSS J0924+0219

Zooming into the broad line region of the gravitationally lensed quasar Q2237+0305 = the Einstein Cross: III. Determination of the size and structure of the CIV and CIII] emitting regions using microlensing

Direct detection of galaxy stellar halos : NGC 3957 as a test case

First case of strong gravitational lensing by a QSO : SDSS J0013+1523 at z = 0.120

COSMOGRAIL: the COSmological MOnitoring of GRAvItational Lenses VII. Time delays and the Hubble constant from WFI J2033-4723

Slit and integral-field optical spectroscopy of the enigmatic quasar HE0450-2958

The geometry of the quadruply imaged quasar PG 1115+080; implications for Ho