🤖 AI Summary
Conventional modeling of dark matter halo density profiles relies heavily on N-body simulations and baryonic physics assumptions, introducing theoretical biases. Method: We propose an entirely data-driven approach based on Exhaustive Symbolic Regression (ESR) to infer halo density profiles directly from weak gravitational lensing observations, objectively balancing fitting accuracy and model simplicity without presupposing functional forms. Contribution/Results: Using simulated galaxy cluster surface mass density data, we systematically assess the impact of measurement uncertainty and sample size: with 5% observational error, as few as 20 clusters suffice to robustly recover the Navarro–Frenk–White (NFW) profile at high confidence; under typical current-survey uncertainties, although simpler functions are preferred, the NFW profile remains statistically competitive. This is the first method to construct halo structural models solely from observational data—free of prior functional assumptions—thereby establishing a new paradigm for empirically testing fundamental dark matter properties.
📝 Abstract
Dark matter haloes are typically characterised by radial density profiles with fixed forms motivated by simulations (e.g. NFW). However, simulation predictions depend on uncertain dark matter physics and baryonic modelling. Here, we present a method to constrain halo density profiles directly from observations using Exhaustive Symbolic Regression (ESR), a technique that searches the space of analytic expressions for the function that best balances accuracy and simplicity for a given dataset. We test the approach on mock weak lensing excess surface density (ESD) data of synthetic clusters with NFW profiles. Motivated by real data, we assign each ESD data point a constant fractional uncertainty and vary this uncertainty and the number of clusters to probe how data precision and sample size affect model selection. For fractional errors around 5%, ESR recovers the NFW profile even from samples as small as 20 clusters. At higher uncertainties representative of current surveys, simpler functions are favoured over NFW, though it remains competitive. This preference arises because weak lensing errors are smallest in the outskirts, causing the fits to be dominated by the outer profile. ESR therefore provides a robust, simulation-independent framework both for testing mass models and determining which features of a halo's density profile are genuinely constrained by the data.