🔭 astrophysics

Non-spherical BUFFALOs: a weak lensing view of the Frontier Field clusters and associated systematics

This study utilizes high-resolution weak lensing data from the BUFFALO survey to quantify systematic biases in mass measurements of the complex Frontier Fields galaxy clusters, revealing that disturbed systems like Abell 2744 are most affected and establishing a crucial framework for future cosmological surveys.

Original authors: A. Niemiec, A. Acebron, B. Beauchesne, M. Jauzac, J. M. Diego, D. Eckert, D. Harvey, A. M. Koekemoer, D. J. Lagattuta, M. Limousin, G. Mahler, N. Patel, S. Tam, J. F. V. Allingham, R. Cen, A. Faisst

Published 2026-02-09

📖 5 min read🧠 Deep dive

CC BY 4.0

Original authors: A. Niemiec, A. Acebron, B. Beauchesne, M. Jauzac, J. M. Diego, D. Eckert, D. Harvey, A. M. Koekemoer, D. J. Lagattuta, M. Limousin, G. Mahler, N. Patel, S. Tam, J. F. V. Allingham, R. Cen, A. Faisst, D. Perera, M. Sereno

Original paper licensed under CC BY 4.0 (http://creativecommons.org/licenses/by/4.0/). ✨ This is an AI-generated explanation of the paper below. It is not written or endorsed by the authors. For technical accuracy, refer to the original paper. Read full disclaimer

The Big Picture: Weighing the Universe's Heavyweights

Imagine the universe is a giant, invisible web. At the intersections of this web sit galaxy clusters—massive cities containing hundreds or thousands of galaxies, held together by gravity. These clusters are the heaviest objects in the universe.

Why do we care about their weight? Because the number of these "cosmic cities" and how heavy they are tells us the recipe of the universe. It helps us understand how much "stuff" (matter) exists and how the universe is expanding. But to get the recipe right, we need to weigh these clusters accurately.

The Problem: The "Distorted Mirror" Effect

The paper focuses on six specific, super-massive clusters (like Abell 2744 and MACS J0717). These aren't calm, round balls of galaxies; they are chaotic, merging messes, often called "non-spherical buffalos" because they are huge, wild, and irregular.

To weigh them, the scientists use a trick called Gravitational Lensing.

The Analogy: Imagine looking at a distant streetlamp through a wavy, distorted piece of glass (like a funhouse mirror). The light from the lamp bends. By measuring how much the light bends, you can calculate how heavy the glass is.
The Reality: Massive galaxy clusters act as that wavy glass. They bend the light from galaxies behind them. By measuring the tiny distortions in the shapes of those background galaxies, scientists can calculate the mass of the cluster in front.

The Data: A High-Definition View

The researchers used the Hubble Space Telescope and a special program called BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations).

The Analogy: Previous maps of these clusters were like looking at a city from a low-flying plane. BUFFALO gave them a view from a high-altitude drone, capturing a much wider area with incredible detail.
The Result: They found about 50 background galaxies per tiny patch of sky (arcminute). This is a very high density, giving them a lot of data points to work with.

The Challenge: The "Messy" Clusters

The paper's main goal was to see if our standard ways of weighing these clusters work when the clusters are messy.

1. The "Center" Problem (Where is the middle?)

The Analogy: If you try to weigh a round pizza, you put it on the scale's center. Easy. But if you try to weigh a pizza that has been torn in half and the two halves are drifting apart, where do you put it on the scale? If you pick the wrong spot, your weight measurement is wrong.
The Finding: For calm clusters, it doesn't matter much where you pick the center. But for the "wild" clusters (like Abell 2744, which is actually two clusters crashing into each other), picking the wrong center (like the brightest galaxy vs. the X-ray peak) can throw off the mass calculation significantly.

2. The "Background Noise" Problem (Who is in the photo?)

The Analogy: Imagine trying to measure the distortion of a mirror by looking at reflections. If you accidentally include reflections of people standing in front of the mirror (who aren't being distorted), your math will be wrong. You need to make sure you are only looking at the people behind the mirror.
The Finding: The team had to be very careful to filter out galaxies that were too close to the cluster (foreground) or part of the cluster itself. If they didn't, the signal gets "diluted," and the cluster looks lighter than it really is. They developed a "boost factor" to fix this mathematically.

3. The "Shape" Problem (Are they round?)

The Analogy: Most weighing formulas assume the object is a perfect sphere. But these clusters are more like lumpy potatoes or flattened pancakes.
The Finding: The team tested if assuming a sphere was okay. They found that for the total mass, the shape doesn't matter as much as we thought (the error is small). However, the assumption that the cluster is a single, smooth object breaks down for the wildest clusters.

The "Double-Cluster" Dilemma

The most interesting discovery involves Abell 2744.

The Analogy: Imagine two cars crashing and sticking together. Is it one giant, heavy wreck, or two smaller cars?
The Finding: When the scientists modeled Abell 2744 as one single object, they got a certain weight. When they modeled it as two separate objects crashing, they got a total weight that was surprisingly similar. This raises a philosophical question for astronomers: How do we define a "cluster"? Is a merging pair one cluster or two? This matters because if we count them wrong in our universe surveys, our cosmological calculations could be off.

The Conclusion: What Did They Learn?

The paper concludes that while our standard methods work well for calm, round clusters, they get tricky with the "wild" ones.

The biggest errors come from the most disturbed, merging clusters.
The solution isn't just better telescopes, but better math that accounts for messy centers and merging shapes.

This work is a "stress test" for the tools astronomers will use in the future with massive new surveys (like the Euclid mission and the Rubin Observatory). Before we can trust the results of those future surveys, we need to know that our "scales" work even when the objects we are weighing are chaotic, non-spherical, and crashing into each other.

In short: The scientists used high-definition space photos to weigh six chaotic galaxy clusters. They found that while we can weigh them, we have to be very careful about where we center our measurement and how we define the object, especially when the clusters are in the middle of a cosmic crash.

Problem
Galaxy clusters serve as critical tracers of large-scale structure, with their mass function evolution constraining key cosmological parameters such as matter density and the amplitude of density fluctuations ( $\sigma_8$ ). However, achieving precision cosmology requires accurate total mass estimates, which are challenging for clusters with complex, disturbed morphologies. Standard weak-lensing analyses often rely on simplifying assumptions, such as spherical symmetry and unimodal mass distributions (e.g., fitting Navarro-Frenk-White or NFW profiles to azimuthally averaged shear). While effective for relaxed systems, these assumptions may introduce significant model biases in highly disturbed, merging clusters. As upcoming wide-field surveys (e.g., Euclid, LSST) prepare to detect tens to hundreds of thousands of clusters, understanding the systematic uncertainties inherent in these standard methods for complex systems is essential.

Methodology
The authors analyze a sample of six extremely massive and morphologically complex galaxy clusters from the Hubble Frontier Fields (HFF), extended by the BUFFALO (Beyond Ultra-deep Frontier Fields And Legacy Observations) survey. The study utilizes high-resolution, multi-band HST imaging (ACS and WFC3) covering several Mpc around each cluster, combined with ancillary X-ray data and spectroscopy.

The methodology proceeds through the following steps:

Catalogue Construction: Weak-lensing catalogues were generated using the pyRRG code on ACS/F814W coadded images. This involved source extraction, star/galaxy separation via random forest classifiers, and precise shape measurements corrected for the Point Spread Function (PSF) using Tiny Tim models.
Source Selection: Background galaxies were selected using calibrated color-color and color-magnitude cuts to remove foreground and cluster-member contaminants. For galaxies without spectroscopic redshifts, photometric redshifts were utilized where available, though the study found color cuts to be more robust for signal purity.
Signal Measurement: The excess surface mass density ( $\Delta\Sigma$ ) and reduced excess surface mass density ( $\Delta\Sigma_r$ ) profiles were measured in radial bins. The analysis accounted for missing redshift information by drawing random redshifts from an empirically derived $N(z)$ distribution, performing 100 realizations with bootstrapping to estimate statistical uncertainties.
Mass Modeling: Cluster masses ( $M_{200}$ $M_{200}$ ) were estimated by fitting NFW profiles to the shear signals. The study tested various modeling assumptions, including:
- Different cluster centers (X-ray peaks vs. Brightest Cluster Galaxies).
- Mass-concentration ( $M-c$ ) relations (DF23 vs. Duffy et al. 2008) to constrain the concentration parameter.
- Combined Strong + Weak lensing (SL+WL) models using detailed core mass maps.
- Alternative density profiles (Singular Isothermal Sphere, SIS).
- Corrections for foreground contamination (boost factors) and the validity of the weak-lensing approximation.

Key Contributions

High-Density Catalogues: The paper produces and makes available weak-lensing catalogues with exceptionally high source densities ( $\sim50$ sources/arcmin $^2$ ), surpassing the projected densities of future stage-IV surveys.
Systematic Analysis: It provides a comprehensive, data-driven assessment of systematic biases affecting mass estimates in complex clusters. The study isolates specific sources of error, including centering offsets, redshift distribution uncertainties, foreground contamination, and the impact of multi-modal mass distributions.
Model Validation: By applying standard parametric fitting methods to a sample of known complex systems, the work quantifies the deviation between simple models and the true mass distributions of merging clusters.

Results

Centring Sensitivity: The choice of cluster center significantly impacts mass estimates, particularly for highly disturbed, multi-modal systems like Abell 2744 (A2744) and MACS J0717.5+3745 (M0717). For these systems, offsets between the X-ray peak and the true mass center (or BCGs) introduce measurable biases. However, the study demonstrates that these biases can be mitigated by fitting off-centered halo models.
Redshift Distribution: The mass estimates are robust to variations in the assumed background redshift distribution $N(z)$ , provided the distribution is reasonably accurate. Using a simple mean redshift or a log-normal fit yields consistent results within $1\sigma$ , whereas poor fits (e.g., original Smail & Dickinson parameters) can introduce bias.
Contamination: Foreground and cluster-member contamination dilutes the lensing signal, leading to underestimated masses. The study confirms that applying a boost factor correction effectively recovers the fiducial mass estimates.
Model Bias: For the most disturbed clusters, single-halo NFW models may not fully capture the mass distribution. For A2744, a two-halo model was tested, revealing that the total mass of the merging system is consistent with a single-halo fit centered on the primary substructure, raising questions about the definition of a "cluster" in post-merger states.
Approximation Validity: The weak-lensing approximation (ignoring convergence in the reduced shear) leads to slight mass overestimates for the most massive clusters, though the effect is generally small compared to other systematics.
Comparison with SL: Strong-lensing-only masses are consistently higher than weak-lensing or combined SL+WL masses, likely due to the difficulty of extrapolating inner-core constraints to larger radii.

Significance
The paper establishes a framework for assessing the impact of systematics on weak-lensing mass estimates, specifically for unrelaxed clusters. It identifies that the largest sources of systematic bias arise in the most disturbed, merging systems. These findings are critical for forthcoming cosmological analyses based on wide-field surveys like Euclid and the Rubin Observatory Legacy Survey of Space and Time (LSST), where detailed modeling of individual systems will be impractical, and statistical stacking will be the primary tool. The work suggests that while standard methods perform well for relaxed clusters, the treatment of complex, merging systems requires careful consideration of centering and morphology to avoid propagating biases into cosmological constraints. The authors conclude that the definition of a cluster in such systems (single vs. multiple halos) remains an open question with potential implications for cosmological parameter derivation.