---
title: 'LensMC: Bayesian Shear Measurement for Euclid'
url: https://www.emergentmind.com/topics/lensmc
type: topic
---

# LensMC: Bayesian Shear Measurement for Euclid

LensMC is a weak-lensing shear measurement method developed for Euclid and Stage-IV surveys. It is based on forward modelling of galaxy images, explicit convolution with the point spread function (PSF), Markov Chain Monte Carlo (MCMC) sampling of galaxy-parameter posteriors, and marginalisation over nuisance parameters. In the Euclid programme, LensMC has a dual role: it is the designated galaxy-shape measurement method for DR1, and it has already been exercised in both pre-launch Euclid-like simulations and on early on-sky data, including the Abell 2390 Early Release Observations and the Euclid Quick Release 1 cluster-lensing catalogue [2405.00669] [2507.07629] [2606.20829].

## 1. Scientific role within Euclid weak lensing

LensMC was designed for the Euclid VIS imager, where the PSF is diffraction-limited, chromatic, spatially varying across the field, and undersampled, with a size comparable to many galaxies. The method targets the accuracy and precision needed for percent-level dark-energy constraints and is formulated to operate at Euclid scale, namely on the order of \(1.5\) billion galaxies. The original Euclid preparation study therefore framed LensMC as a shear-measurement method tailored to explicit PSF treatment, realistic neighbour handling, and statistically grounded uncertainty propagation [2405.00669].

Within Euclid’s operational weak-lensing pipeline, LensMC has moved from simulation validation to survey deployment. In the Abell 2390 Early Release Observations analysis, it was used alongside KSB+ and SourceXtractor++ as one of three independent shape-measurement algorithms, and that study identifies LensMC as Euclid’s designated galaxy-shape measurement method for DR1 [2507.07629]. In Quick Release 1, LensMC was used to generate the first Euclid shear-measurement catalogue produced with the method in anticipation of DR1, covering the three Euclid Deep Field mosaics over approximately \(63\,\mathrm{deg}^2\) [2606.20829].

This positioning is methodologically significant. LensMC is not an internally calibrated shear-response scheme of the metacalibration type, nor a pure moment-based estimator of the KSB family. Instead, it is a Bayesian forward-modelling framework whose calibration strategy relies on dedicated image simulations and, in some data releases, empirical corrections. A plausible implication is that LensMC is intended to be especially useful where PSF convolution, blending, and heterogeneous morphology make pixel-level likelihood modelling preferable to lower-dimensional summaries.

## 2. Forward model, likelihood, and MCMC inference

At its core, LensMC models the observed image as a PSF-convolved galaxy profile plus noise,
\[
I_{\rm obs}(\mathbf{x}) = [P(\mathbf{x}) * M(\mathbf{x};\boldsymbol{\theta})] + n(\mathbf{x}).
\]
The Euclid preparation study describes the galaxy model as a linear mixture of two co-centred circular Sérsic-type components, disc plus bulge, rendered isotropically and then anisotropically distorted. In that baseline description, the disc is exponential, the bulge has fixed \(n_b=1\), the size ratio is fixed to \(r_h/r_e=0.15\), and profiles are truncated at \(r_{\max}/r_e=4.5\). The Abell 2390 analysis also describes LensMC as fitting a two-component bulge+disc profile convolved with the PSF, but there the profile is reported as de Vaucouleurs bulge plus exponential disc [2405.00669] [2507.07629].

The ellipticity enters through the model geometry. In the Euclid preparation paper, the isotropic template is transformed by an anisotropic distortion matrix,
\[
\tens{S}=\frac{\bar{r}_0}{q_\epsilon\,r_0}
\begin{pmatrix}
1-\epsilon_1 & -\epsilon_2 \\
-\epsilon_2 & 1+\epsilon_1
\end{pmatrix},
\qquad q_\epsilon=1-|\epsilon|.
\]
Rendering and convolution are performed in Fourier space. For isotropic \(I(r)\), the 2D Fourier transform reduces to the 1D Hankel transform
\[
\tilde{I}(k)=2\pi\!\int_0^\infty I(r)\,J_0(kr)\,r\,{\rm d}r,
\]
after which shear, centroid shifts, and PSF convolution are applied efficiently before inverse transformation to real space [2405.00669].

Inference is Bayesian. With pixel data vector \(\vec{D}\), ellipticity \(\epsilon\), non-linear nuisance parameters \(\theta\), and linear flux parameters \(\phi\), LensMC defines
\[
p(\epsilon,\theta,\phi|\vec{D}) = \frac{p(\vec{D}|\epsilon,\theta,\phi)\,p(\epsilon,\theta,\phi)}{p(\vec{D})}.
\]
Assuming Gaussian pixel noise,
\[
\ln p(\vec{D}|\epsilon,\theta,\phi)= -\frac{1}{2}[\vec{D}-\vec{I}]^\top C^{-1} [\vec{D}-\vec{I}] + \mathrm{const}.
\]
Because the model is linear in the bulge and disc fluxes, those fluxes can be marginalised analytically. LensMC then estimates ellipticity using the posterior mean,
\[
\hat{\epsilon}=\int \epsilon\,p(\epsilon|\vec{D})\,d\epsilon.
\]
This posterior-mean estimator is central to the method’s stated robustness against overfitting and non-Gaussian posteriors [2405.00669].

Sampling is performed with an improved Metropolis–Hastings scheme. The preparation study reports an initial deterministic maximisation, a burn-in phase with a cooling schedule, adaptive proposal-covariance updates every \(100\) samples, \(500\) burn-in samples, and \(N_{MC}=200\) kept samples. Joint likelihoods are used for neighbour groups, typically with \(N=2\), so that overlapping objects are fitted simultaneously rather than circularised by masking or independent treatment [2405.00669].

## 3. Catalogue construction, blending control, and object-level outputs

LensMC produces per-object shape catalogues rather than only stacked shear estimates. In the Abell 2390 analysis, the method returns ellipticity components \(e_1,e_2\), sizes, fluxes or magnitudes, positions, parameter uncertainties, reduced \(\chi^2\), and a shear weight per object, together with quality flags; flagged objects can be assigned zero weight [2507.07629]. In Q1, the merged catalogue includes right ascension and declination, ellipticities, the flux-averaged half-light radius \(r_{\rm hl}\), \(I_E\) magnitude, per-object weight \(w\), quality flags, and selected MER and PHZ columns [2606.20829].

Blending control is an explicit part of the LensMC workflow. The Euclid preparation paper groups neighbours with a friends-of-friends scale \(r_{\rm friend}=1^{\prime\prime}\), uses segmentation maps dilated by one pixel to mask objects outside the group, and jointly fits close pairs to mitigate neighbour bias [2405.00669]. The Abell 2390 implementation likewise groups nearby objects with friends-of-friends at \(1\) arcsec scale, jointly fits blended systems, applies robust sigma-clipping to mitigate residual cosmic rays and image features, uses an external segmentation map and mask, and subtracts local background gradients across \(512\)-pixel postage stamps [2507.07629]. In Q1, deblending and spurious-detection control rely on the MER pipeline flags and a pre-trained spurious classifier, while LensMC is run on all MER detections and science analyses impose catalogue-level quality cuts [2606.20829].

Star–galaxy separation and robustness cuts differ across analyses. In the preparation simulations, star–galaxy separation uses a measured size threshold \(r_{s/g}=0.15^{\prime\prime}\), with true-positive rate approximately \(93\%\) and false-positive rate approximately \(5\%\) at \(I_E<26.5\) [2405.00669]. In the Abell 2390 ERO analysis, objects are retained if the flux-averaged half-light radius exceeds \(0.09^{\prime\prime}\), and faint galaxies with unrealistically large size estimates are removed using the magnitude-dependent bound
\[
r_{\rm eff} < [-0.1875^{\prime\prime}\times(I_E-24)+1.85^{\prime\prime}].
\]
In Q1, star–galaxy separation uses the magnitude-independent threshold \(r_{\rm hl}>0.07^{\prime\prime}\), and analyses may optionally impose \(r_{\rm hl}<2^{\prime\prime}\) to suppress a large-size tail that is mostly spurious or poorly modelled [2507.07629] [2606.20829].

Reported source densities illustrate the method’s survey role. In Euclid-like simulations LensMC measured objects with density approximately \(90\,\mathrm{arcmin}^{-2}\) for \(I_E<26.5\) across \(4500\,\mathrm{deg}^2\) [2405.00669]. In Abell 2390, after masking and selections, the catalogue reached \(110\,\mathrm{arcmin}^{-2}\) in total, \(82\,\mathrm{arcmin}^{-2}\) in the weak-lensing magnitude range \(22<I_E<26.5\), and \(22\,\mathrm{arcmin}^{-2}\) for \(I_E<24.5\) [2507.07629]. In Q1, the catalogue is reported as complete to \(I_E \approx 25\), with \(26\,\mathrm{arcmin}^{-2}\) for \(I_E<24.5\) and \(75\,\mathrm{arcmin}^{-2}\) for \(I_E<27\) after quality cuts [2606.20829].

## 4. Bias model, calibration strategy, and systematics

LensMC adopts the standard linear shear-bias model
\[
\hat{g}_i=(1+m_i)\,g_i+c_i+n_i,
\]
with multiplicative bias \(m_i\), additive bias \(c_i\), and noise \(n_i\). PSF leakage is parameterised by
\[
\alpha_i=\frac{d c_i}{d \epsilon_{{\rm PSF},i}}.
\]
The Euclid preparation study explicitly separates measurement biases from detection and selection biases, rather than treating the full pipeline response as a single number [2405.00669].

In Euclid-like Flagship and VIS emulation, LensMC’s measurement-only biases were reported as
\[
m_1=(-3.6\pm0.2)\times10^{-3},\qquad
m_2=(-4.3\pm0.2)\times10^{-3},
\]
\[
c_1=(-1.78\pm0.03)\times10^{-4},\qquad
c_2=(0.09\pm0.03)\times10^{-4}.
\]
The same study found a large detection bias with multiplicative component \(1.2\times10^{-2}\) and additive component \(-3\times10^{-4}\), together with measurement PSF leakage
\[
\alpha_1=(-9\pm3)\times10^{-4},\qquad
\alpha_2=(2\pm3)\times10^{-4}.
\]
A dominant contribution to the measurement multiplicative bias comes from undetected faint galaxies, with \(m_{\rm faint}\approx-5\times10^{-3}\). Morphology mismatch can add approximately \(-8\times10^{-3}\) of multiplicative bias when full bulge variability and a bulge-only subpopulation are allowed, but the study characterises this model bias as straightforward to calibrate because of weak sensitivity to the assumed distributions [2405.00669].

For the Abell 2390 cluster-regime application, calibration was performed with dedicated Flagship-based image simulations extending to \(|g| \le 0.2\). In that analysis, LensMC used a refined linear multiplicative correction appropriate for native PSF sampling, with multiplicative bias \(m \approx 0.0498\) and additive bias \(c \approx 1.24\times10^{-3}\). A conservative multiplicative-bias uncertainty of \(3.0\%\) was adopted, with \(1.2\%\) added for PSF SED-dependence not modelled, giving a total shear-calibration uncertainty of \(3.2\%\); the paper states that this dominates the weak-lensing systematic budget for that single-target study [2507.07629].

The Q1 catalogue uses an empirical additive correction rather than a full per-object multiplicative calibration. Field-of-view averaged additive biases are reported as \(c_1 \approx -2.35\times10^{-3}\) and \(c_2 \approx +1.79\times10^{-3}\), with no significant spatial dependence detected over the Q1 footprint. The adopted correction is a map \(c=c(\mathrm{size},I_E)\), built by binning in \(r_{\rm hl}\) and \(I_E\); after correction, two-point statistics show suppressed excess power at \(\theta \gtrsim 10\) arcmin, and B-modes are an order of magnitude below E-modes in the aperture-mass test. By contrast, Q1 does not derive a per-object \(m(\mathrm{size},I_E)\) calibration, and for cluster science the impact is described as sub-dominant [2606.20829].

A recurrent methodological point is that LensMC’s calibration philosophy differs from internally calibrated methods. Relative to metacalibration and BFD-style approaches, it relies on realistic external simulations, and relative to KSB-style moment methods it absorbs PSF convolution and undersampling directly into the forward likelihood. This suggests that the method is best understood as a calibrated model-fitting estimator rather than a self-calibrating one [2405.00669] [2606.20829].

## 5. Tomographic weak-lensing analysis of Abell 2390

The Abell 2390 Early Release Observations paper provides the first end-to-end demonstration of LensMC in Euclid tomographic cluster weak lensing. The analysis combines Euclid VIS imaging in the broad \(I_E\) band over \(540\)–\(920\) nm and NISP \(Y,J,H\) imaging with Subaru/Suprime-Cam \(B,V,R_c,i,I_c,z^\prime\) data and CFHT Megacam \(u\)-band imaging, although the \(u\)-band is dropped from the photometric-redshift estimation because of unstable calibration under cirrus. The Euclid VIS data were obtained in three dithered ROS sequences for a total VIS exposure of \(7932\) s; VIS stacks have \(0.1^{\prime\prime}\) pixels and \(5\sigma\) depth \(I_E \approx 27.0\), while NISP stacks have \(0.3^{\prime\prime}\) pixels with \(5\sigma\) depths \(Y \approx 25.18\), \(J \approx 25.22\), and \(H \approx 25.12\). The overlap area provides complete azimuthal coverage out to approximately \(3.07\) Mpc, and Galactic cirrus was removed from all stacks with DeNeb [2507.07629].

Photometric redshifts are estimated with Phosphoros using the NISP and Suprime-Cam bands, with the VIS \(I_E\) band excluded because of chromatic detrending in ERO and the CFHT \(u\) band excluded because of unstable offsets under reddening. Systematic photometric offsets are calibrated against limited spectroscopic redshifts at or near the cluster, with an added per-band fractional flux uncertainty of approximately \(5\%\). The reported spectroscopic performance is \(\mathrm{NMAD}\approx0.03\) and outlier rate \(\eta\approx3\%\). Tomographic selection uses magnitude bins \(22<I_E<24.5\) and \(24.5<I_E<26.5\), and redshift bins \((0.2,0.3]\), \((0.3,0.6]\), \((0.6,0.9]\), \((0.9,1.5]\), \((1.5,2.8]\), and \((2.8,6]\); the cluster mass analysis uses the four central bins, omitting the lowest and highest [2507.07629].

The source redshift distributions are calibrated with a self-organising map approach using COSMOS2020. One SOM is trained on the Abell 2390 photometric space, COSMOS calibrators are assigned to the same cells, and per-cell mini-\(n(z)\) distributions are weighted by the sum of LensMC shear weights in the cell. For the bins used in the mass analysis, the resulting mean lensing efficiencies span \(\langle\beta\rangle \approx 0.45\)–\(0.80\) and \(\langle\beta^2\rangle \approx 0.24\)–\(0.65\) [2507.07629].

Cluster-member contamination, obscuration, and magnification are treated explicitly. Source-density profiles show strong central excess at low photometric redshift and central depletion at high photometric redshift. Obscuration is mapped through large-scale image injections at approximately \(29.4\,\mathrm{arcmin}^{-2}\), and magnification is modelled by applying a reference NFW magnification to a Flagship population. The contamination boost factor is fitted to the corrected source-density profiles using
\[
\frac{n(r)}{n_{\rm outer}} = 1 + A \exp(1-r/r_S),
\]
with joint scale radius \(r_S = 350 \pm 39\) kpc and bin-dependent amplitudes \(A\). The uncertainty in the boost correction contributes approximately \(0.7\%\) to the mass systematic [2507.07629].

The lensing analysis uses the reduced shear relation
\[
g=\frac{\gamma}{1-\kappa},
\]
the tangential ellipticity estimator
\[
e_t = -(e_1 \cos 2\phi + e_2 \sin 2\phi), \qquad \langle e_t\rangle \approx g_t,
\]
and the critical surface density
\[
\Sigma_{\rm crit} = \frac{c^2}{4\pi G}\frac{D_s}{D_d D_{ds}}.
\]
Mass modelling adopts a spherical NFW profile,
\[
\rho(r)=\frac{\rho_s}{(r/r_s)(1+r/r_s)^2},
\]
with fixed \(c_{200c}=4.0\). The tangential reduced shear is jointly fitted in eight bin combinations, corresponding to four photometric-redshift bins times two magnitude bins, over \(0.5 \le r \le 3.3\) Mpc. The centre is fixed to the strong-lensing/MCMC mass centre at \(\mathrm{RA}(J2000)=328.40494^\circ\), \(\mathrm{Dec}(J2000)=17.69460^\circ\), approximately \(6^{\prime\prime}\) south-east of the brightest cluster galaxy. A Wiener-filtered convergence reconstruction shows a high-significance E-mode detection aligned south-east to north-west and peaking at the BCG, with negligible B-mode residual; the paper emphasises, however, that quantitative constraints come from the \(g_t\) profile fits [2507.07629].

The inferred masses are as follows.

| Method | \(M_{200c}\) | \(M_{500c}\) |
|---|---:|---:|
| LensMC | \((16.0 \pm 1.2 \pm 1.7)\times10^{14}\,M_\odot\) | \((11.1 \pm 0.8 \pm 1.2)\times10^{14}\,M_\odot\) |
| KSB+ | \((14.9 \pm 1.2 \pm 1.6)\times10^{14}\,M_\odot\) | \((10.3 \pm 0.9 \pm 1.1)\times10^{14}\,M_\odot\) |
| SourceXtractor++ | \((14.5 \pm 1.1 \pm 1.6)\times10^{14}\,M_\odot\) | \((10.0 \pm 0.7 \pm 1.1)\times10^{14}\,M_\odot\) |

The LensMC result is consistent with KSB+ at approximately \(1.2\sigma\) once shear-calibration systematics are included, and all three methods agree within the quoted uncertainties. Individual tomographic-bin mass fits scatter around the joint best fit without significant trends versus photometric redshift or magnitude; one bin, the bright \((1.5,2.8]\) sample, lies approximately \(2\sigma\) high, which the paper describes as not unexpected among eight bins. The cross-component \(g_\times\) is consistent with zero for all methods. For Abell 2390 at \(z_d=0.228\), the LensMC mass is also in good agreement with earlier wide-field weak-lensing results from WtG, LoCuSS, and CCCP/MENeaCS, while the Euclid analysis reaches tighter mass uncertainties, approximately \(14\%\) total statistical, primarily because of higher source density and tomographic selection [2507.07629].

## 6. Quick Release 1 catalogue, validation, and near-term outlook

The Q1 LensMC catalogue extends the method from a single-cluster demonstration to a survey-scale cluster-lensing resource. It covers the Euclid Deep Field North, South, and Fornax mosaics, amounting to approximately \(63\,\mathrm{deg}^2\) of VIS stacked images, with accompanying photometric redshifts. In Q1, the PSF is reconstructed from the released VIS stacked-image PSF grids, with grid pitch approximately \(12\) arcmin over approximately \(32\)-arcmin stacks. The adopted strategy uses PSF cutouts, flagging around the core, \(3\times\) oversampling with splines to avoid undersampling bias, and bilinear interpolation, falling back to nearest-neighbour interpolation when needed [2606.20829].

For stacked cluster lensing, Q1 defines the tangential and cross ellipticity components in the field-of-view aligned frame as
\[
e_t = - e_1\cos(2\phi) - e_2\sin(2\phi),\qquad
e_\times = e_1\sin(2\phi) - e_2\cos(2\phi),
\]
and the annular estimators
\[
\gamma_t(R)=\frac{\sum_i w_i e_{t,i}}{\sum_i w_i},\qquad
\gamma_\times(R)=\frac{\sum_i w_i e_{\times,i}}{\sum_i w_i}.
\]
For excess surface density in comoving units, the analysis uses
\[
\Sigma_{\rm crit}^{(\mathrm{com})}=\frac{c^2}{4\pi G}\frac{D_s}{D_d D_{ds}}(1+z_d)^2,
\]
\[
w_{\Delta\Sigma,i}=\Sigma_{{\rm crit},i}^{-2} w_i,
\]
and
\[
\Delta\Sigma_{g,t}(R)=\frac{\sum_i w_i \Sigma_{{\rm crit},i}^{-1} e_{t,i}}{\sum_i w_i \Sigma_{{\rm crit},i}^{-2}}.
\]
The recommended baseline cuts are \(\mathrm{DET\_QUALITY\_FLAG}\in\{0,1,2,512\}\), \(\mathrm{SPURIOUS\_FLAG}=0\), \(r_{\rm hl}>0.07^{\prime\prime}\), optional \(r_{\rm hl}<2^{\prime\prime}\), and \(\mathrm{PHZ\_FLAGS}\in\{0,12\}\), together with a background condition such as \(z_s>z_d+0.1(1+z_d)\) [2606.20829].

Validation in Q1 proceeds through both internal and external tests. Two-point statistics on \(I_E<24.5\) sources recover the expected trends, although excess \(\xi_-\) power remains at large angular scales when using uncalibrated photo-\(z\) values and nominal cosmology; the additive correction suppresses but does not fully eliminate that excess. In stacked cluster analyses, \(495\) MaDCoWS2 candidates in the Q1 footprint show good consistency between raw and corrected catalogues out to \(R_c \approx 20\) Mpc, with \(\gamma_\times\) consistent with zero. Cross-validation against DES redMaPPer clusters yields excellent agreement of LensMC and DES tangential shear over \(0.2\) to several Mpc; the differential biases are consistent with zero, \(\delta m = (-5 \pm 8)\times10^{-2}\) and \(\delta c = (2 \pm 1)\times10^{-4}\), and improve further when the farthest bins are removed. Random cluster stacks are null in most bins, but low-redshift large-radius bins show residuals consistent with edge effects, incomplete azimuthal coverage, or PSF residuals, so the Q1 paper recommends downweighting or omitting those bins in precision analyses [2606.20829].

The scientific reach of the Q1 catalogue is already substantial. The paper states that, thanks to Euclid image resolution, depth, and overall control of systematic errors, stacked lensing profiles of clusters with masses \(10^{14}M_\odot\) can be constrained out to \(z\approx2\), spanning nearly \(10\) Gyr of evolution history. At the same time, the catalogue is explicitly transitional: Q1 PSF modelling is stack-based and grid-interpolated; no per-object multiplicative calibration is provided; and residual systematics remain at low redshift and large radius. The stated DR1 direction is a wavefront-based PSF forward model per exposure, full multiplicative-bias calibration from image simulations, refined magnitude–size and star–galaxy selections, and improved photo-\(z\) calibration and \(\Sigma_{\rm crit}\) uncertainty propagation, with survey area projected to reach approximately \(1900\,\mathrm{deg}^2\) by the end of 2026 [2606.20829].

Taken together, the simulation paper, the Abell 2390 ERO study, and the Q1 catalogue delineate LensMC as a Euclid weak-lensing framework with three defining properties: explicit pixel-level forward modelling, MCMC posterior sampling with neighbour-aware fitting, and a calibration strategy anchored in realistic simulations plus release-specific empirical corrections. The available evidence shows that this combination yields high-density shape catalogues, internally consistent tomographic cluster masses, and robust stacked cluster-lensing measurements, while also making clear that shear calibration, PSF modelling, and low-level selection effects remain the dominant technical levers for DR1 and subsequent releases [2405.00669] [2507.07629] [2606.20829].

Source: https://www.emergentmind.com/topics/lensmc