Euclid Quick Data Release (Q1): The Strong Lensing Discovery Engine

Euclid Q1 Strong Lensing Discovery Engine [image to be added on Wed 19th] This repository shares catalogues, images, and modelling results for strong lens candidates found in the Euclid Q1 data release as part of the Strong Lensing Discovery Engine. Lenses were found through a combined machine learning, citizen science, and expert visual inspection campaign. The sample is described and presented in the following five discovery papers: (paper links to be added on Wed 19th). (candidates from paper B are listed separately) If you use Euclid Q1 data in your research, please cite the five core Euclid Q1 papers. Euclid is a space telescope surveying 14k square degrees over the next six years. Q1 is the first major data release, covering 63 square degrees in the Euclid Deep Fields (see the paper). It is representative of the ongoing Euclid Wide Survey. Euclid images include a VIS optical image (400-900nm) and 3 NIR images (Y, J, H). The Euclid VIS PSF FWHM is approx. 0.16 arcsec per pixel. We are continuing to search for strong lenses in the ongoing survey. We expect around 10k candidates in Euclid DR1. Strong Lens Candidate Catalogue The list of strong lens candidates is strong_lens_candidates.csv. This table combines three sets of candidates, identified with the subset column. - Most (250 grade A, 247 grade B) candidates are from the discovery_engine subset. These were identified by our main search and follow the selection reported in our papers. - We also include 56 grade A and 19 grade B candidates identified by Galaxy Zoo Euclid volunteers as part of a campaign to measure visual morphology; these are outside of the Q1 footprint and have a different selection function. - Finally, we include 3 grade A candidates missed by our search but identified serendiptously. These are inside Q1. All subsets have filled values for the following columns subset "discovery_engine" for main search sample (in Q1). "gz_euclid" for lenses found by Galaxy Zoo Euclid volunteers (outside Q1). "q1_serendipitous" for lenses missed by the main search but found serendipitously (in Q1) id_str Primary key for matching with modelling tables below ("discovery_engine" subset only) tile_index Euclid tile index i.e. which MER tile hosts this galaxy right_ascension of the host galaxy, in degrees. May be offset by up to 10''. declination of the host galaxy, in degrees. May be offset by up to 10''. expert_score Numerical score where higher indicates more likely to be a lens grade Suggested binning of expert_score above. A where > 2, B where > 1, C otherwise. expert_total_votes Experts asked per galaxy (typically 10) notes Light notes for interesting lenses e.g. double source plane lens. references Papers who also use Euclid data to note these lenses (e.g. NGC6505) The discovery_engine subset also includes filled values, copied from the MERge Q1 catalogue, for object_id Euclid object id i.e. the MER catalogue identifier for this galaxy segmentation_area Number of pixels included in SourceExtractor++ mask of galaxy (0.1 arcsec/pixel). flux_vis_1fwhm_aper VIS flux within an aperture of radius 1 FWHM (micro-janksy) For the other subsets, each source was not originally drawn from the MERge Q1 catalogue, and so these values are not filled. Lens Modeling We attempted to fit lens models to all candidates identified by the Discovery Engine (discovery_engine in subset column above) and rated above 2.0 (Grade B or better). The rest of this README discusses those results. Modeling was performed by fitting the VIS image as the "primary" wavelength image and then using that mass model to reconstruct the source in the NIR images. The following lens model is used: - A Multi Gaussian Expansion (MGE) lens light model (see He et al. 2024). - An SIE with external shear mass model. - A source galaxy reconstructed using an adaptive Voronoi mesh (see PyAutoLens docs). The Euclid Strong Lens Modeling Pipeline is publicly available here. It is built with the modeling code PyAutoLens. Modelling Parameter Catalogues 336 candidates were successfully fit with a lens model. We hope that nearly all are strong lenses. Catalogues of lens modeling results are provided in five .csv files. Each .csv can be joined with the lens catalogue (q1_discovery_engine_lens_catalog.csv) on id_str. All results for all parameters are given using six different values: - `max_lh`: Maximum likelihood value. - `median_pdf`: Median value of the posterior PDF. - `upper_3_sigma`: Upper 3-sigma value (99.7% percentile). - `lower_3_sigma`: Lower 3-sigma value (0.3% percentile). - `upper_1_sigma`: Upper 1-sigma value (68.3% percentile). - `lower_1_sigma`: Lower 1-sigma value (31.7% percentile). The file `modelling_lens_mass.csv` summarizes the lens mass model properties of all 336 high-quality candidates: einstein_radius_effective The effective Einstein radius of the lens galaxy ($R_{\mathrm{Ein,eff}}$) centre_0 The y-coordinate of the lens galaxy mass centre ($y_{\text{centre}}$) centre_1 The x-coordinate of the lens galaxy mass centre ($x_{\text{centre}}$) ell_comps_0 The first component of the ellipticity of the lens galaxy mass distribution ($\epsilon_{1}$) ell_comps_1 The second component of the ellipticity of the lens galaxy mass distribution ($\epsilon_{2}$) einstein_radius The internal Einstein radius normalization used by PyAutoLens, of the lens galaxy mass distribution ($\theta_{\text{E}}^{\text{mass}}$) shear_gamma_1 The first shear component of the external shear ($\gamma_{1}^{\text{ext}}$) shear_gamma_2 The second shear component of the external shear ($\gamma_{2}^{\text{ext}}$) The file modelling_lens_mge.csv summarizes the MGE lens light model properties. centre_0 The y-coordinate of the MGE lens galaxy light centre ($y_{\text{centre}}$) centre_1 The x-coordinate of the MGE lens galaxy light centre ($x_{\text{centre}}$) ell_comps_0 The first component of the ellipticity of the MGE lens galaxy light distribution ($\epsilon_{1}$): ell_comps_1 The second component of the ellipticity of the MGE lens galaxy light distribution ($\epsilon_{2}$) flux_ratio The flux ratio of the MGE components ($f_{\text{ratio}}$) The file modelling_lens_sersic.csv summarizes the Sersic lens light model properties. centre_0 The y-coordinate of the Sersic lens galaxy light centre ($y_{\text{centre}}$) centre_1 The x-coordinate of the Sersic lens galaxy light centre ($x_{\text{centre}}$) ell_comps_0 The first component of the ellipticity of the Sersic lens galaxy light distribution ($\epsilon_{1}$) ell_comps_1 The second component of the ellipticity of the Sersic lens galaxy light distribution ($\epsilon_{2}$) effective_radius The effective radius of the Sersic lens galaxy light distribution ($R_{\text{eff}}$) sersic_index The Sersic index of the lens galaxy light distribution ($n$) The file modelling_mge_magnitude.csv summarizes the magnitudes of the lens and source galaxies. Includes values for VIS / Y / J / H. max_lensed_source_signal_to_noise_ratio The maximum signal-to-noise ratio of the lensed source ($S/N$) total_lens_flux The total flux of the lens galaxy ($F_{\text{lens}}$) total_lensed_source_flux The total image-plane flux of the lensed source ($F_{\text{lensed}}$) total_source_flux The total source-plane flux of the lensed source ($F_{\text{source}}$) total_magnification The total magnification of the lensed source ($\mu$) lens_magnitude_ab The AB magnitude of the lens galaxy ($M_{\text{lens}}$) lensed_source_magnitude_ab The image-plane AB magnitude of the lensed source ($M_{\text{lensed}}$) source_magnitude_ab The source-plane AB magnitude of the source galaxy ($M_{\text{source}}$) The file modelling_lens_sersic_magnitude.csv contains similar information as above but for the specific lens. Individual Lens Data The four .zip files contain cutouts and modeling results (where successful) for the four sets of candidates described below. Lens The archive lens.zip contains folders for each candidate with a .fits file with the lens's ID ("id_str" in the catalogue), for example `102018665_NEG570040238507752998/102018665_NEG570040238507752998.fits`. Folders also contain various .png and .fits files, including: - `rgb_0.png` and `rgb_1.png`: RGB images of the lens. - `mask_extra_galaxies.fits`: A mask of extra galaxies in the imaging data. - `sie_fit_pix.png`, `source_reconstruction.png`, `sie_fit_mge.png`: Lens modeling results. - `lens_light.fits`, `source_light.fits`, `source_reconstruction.fits`, etc.: Files summarizing model data. Unsuccessful Lens Candidates and Models The archive unsuccess.zip contains images of 144 unsuccessful lens candidates for various reasons, including modeling failures, non-lens status, or ambiguity. Group Scale Lenses The archive group.zip contains images of 41 group-scale lens candidates requiring multiple mass components, which was beyond the scope of the current modeling pipeline. Recenter The archive recenter.zip contains images of lenses which were not centered correctly, breaking lens modeling. We hope to update these in future. Important Caveat: PSF Quality The VIS PSF used for lens modeling is produced by OU-MER, not utilizing the full Euclid PSF model of OU-SHE. While suitable for Einstein radii and source reconstructions, it may be inadequate for more detailed science such as dark matter substructure studies. Appendix A: Modelling Details Geometry Lens modeling results involve several geometric quantities, defined as follows: To begin, a $(y,x)$ Cartesian coordinate system is established, centered on the light or mass profile: $$y_{shift} = y - y_{centre}, \quad x_{shift} = x - x_{centre}$$ The shifted coordinates are then rotated by the position angle $\phi$, defined counter-clockwise from the positive x-axis: $$y_{rot} = y' \cos(\phi) - x' \sin(\phi), \quad x_{rot} = y' \sin(\phi) + x' \cos(\phi)$$ The rotated coordinates are converted to elliptical coordinates using the axis-ratio $q$ of the light or mass profile: $$\xi = \sqrt{(x_{rot})^2 + \frac{(y_{rot})^2}{q^2}}$$ Instead of defining geometry using $\phi$ and $q$, the following "elliptical components" are used: $$\epsilon_{1} = \frac{1 - q}{1 + q} \sin(2\phi), \quad \epsilon_{2} = \frac{1 - q}{1 + q} \cos(2\phi)$$ To convert between these components and the position angle / axis-ratio: $$q = \frac{1 - \text{fac}}{1 + \text{fac}}, \quad \phi = \frac{1}{2} \arctan\left(\frac{\epsilon_{2}}{\epsilon_{1}}\right)$$ where: $$\text{fac} = \sqrt{1 - \epsilon_{1}^2 - \epsilon_{2}^2}$$ Mass Model The SIE lens mass model is defined by: $$\kappa(\xi) = \frac{1}{1 + q^{\text{mass}}} \left( \frac{\theta^{\text{mass}}_{\text{E}}}{\xi} \right)$$ where $\theta^{\text{mass}}_{\text{E}}$ is the Einstein radius of the mass model. Note that this is a scaled Einstein radius to improve modeling efficiency in PyAutoLens. The **effective** Einstein radius, $R_{\text{Ein,eff}}$, which is the definition of Einstein radius more commonly used in the literature, is defined as: $$R_{\text{Ein,eff}} = \sqrt{\frac{A}{\pi}}$$ where $A$ is the area enclosed by the tangential critical curve, ensuring consistency across various mass density profiles. External Shear The external shear field is parameterized by two components $\gamma_{1}^{\text{ext}}$ and $\gamma_{2}^{\text{ext}}$. The magnitude and orientation are given by: $$\gamma^{\text{ext}} = \sqrt{(\gamma_{1}^{\text{ext}})^2 + (\gamma_{2}^{\text{ext}})^2}, \quad \tan(2\phi^{\text{ext}}) = \frac{\gamma_{2}^{\text{ext}}}{\gamma_{1}^{\text{ext}}}$$ Magnitudes For the MGE lens light and Voronoi source reconstruction, refer to the papers and PyAutoLens documentation for detailed descriptions of the profiles. The magnitudes of the lens and source galaxies are calculated in AB magnitudes: $$M_{\text{lens}} = 2.5 \times \log_{10}(F\_{lens}) + Z$$ $$M_{\text{source}} = 2.5 \times \log_{10}(F\_{source}) + Z$$ where $F\_{lens}$ and $F\_{source}$ are the fluxes of the lens and source galaxies, respectively, and $Z$ is the zero-point. The magnification factor $\mu$ is defined as: $$\mu = \frac{F\_{lensed}}{F\_{source}}$$ where $F\_{lensed}$ is the total flux of the lensed source galaxy in the image-plane (e.g. it is magnified by the lens galaxy). The `zero_point` corresponds to the calibration of the data. Sersic Results for a Sersic lens light model are also provided, which is given as: $$I(\xi) = I_{\text{eff}} \exp\left(-b_n \left[ \left( \frac{\xi}{R_{\text{eff}}} \right)^{1/n} - 1 \right] \right)$$ where $I_{\text{eff}}$ is the effective intensity, $R_{\text{eff}}$ is the effective radius, $n$ is the Sersic index, and $b_n$ is a constant that depends on $n$.