A dense dark matter core of the subhalo in the strong lensing system JVAS B1938+666

A non-parametric reconstruction of a dark matter subhalo in the strong-lensing system JVAS B1938+666 reveals that Self-Interacting Dark Matter (SIDM) or Fuzzy/Wave Dark Matter (FDM) profiles provide significantly better fits to the density profile compared to the standard Navarro-Frenk-White (NFW) model, offering insights into potential solutions for the core-cusp discrepancy in dark matter substructure.

Detection and Data Analysis

The study utilized archival near-infrared observations of JVAS B1938+666 obtained with the Keck II telescope’s NIRC2 instrument to analyze a low-mass dark subhalo. The analysis involved a two-step modeling process:

1. Initial main lens and source model fitting using LensCharm, employing a single smooth pseudo-isothermal ellipsoid (PIE) profile for the main galaxy convergence model.

2. Joint optimization including a non-parametric subhalo convergence perturbation, where the total convergence is expressed as: κtot(θ) = κmain(θ) + ∆κsub(θ).

The detection of the subhalo was confirmed by calculating a Bayes information criterion (BIC), yielding a signal detection confidence of approximately 10 σ. The resulting non-parametric dark matter density and mass distribution were listed in Table 3, which includes radius, density ρ at the radius, enclosed mass M(< r) in the inner region of r, and corresponding error σM(<r).

Model Comparison and Fitting

The reconstructed subhalo density profile was fitted to three competing models: (i) the collisionless NFW profile predicted by CDM, (ii) the SIDM profile, and (iii) the solitonic ψDM wave solution. The fitting results across these models are summarized in Table 2. The comparison of the best-fit results shows that:

)&The data moderately favors the SIDM profile model over ψDM with a Bayes factor of 14.44.

)The fitting goodness of the NFW model compared with ψDM is ∆BIC = −186.3, which suggests a significant tension between NFW and data assuming the reconstructed subhalo profile is representative.

Profile Characteristics and Constraints

The analysis of the density profiles yielded specific characteristics for each model:

)&For the ψDM model, we adopted the widely recognized profile from Schive et al. (2014b).

)The dark matter soliton core of ψDM is larger in low-mass halos (i.e., Mh < 109 M⊙) when the dark matter particle mass is fixed.

)The SIDM density profile is described by Equation (A8), where ρ(r) = ρs / (r 4 + r 4c) 1/4 rs-1 + r-2.

Systematic Error Estimation

To assess the robustness of the detection, systematic errors were estimated by simulating two different dark matter subhalos—one NFW and one ψDM—and fitting them to mock strong lensing systems similar to JVAS B1938+666. These simulations demonstrated that:

)The profiles of the two different dark matter halos can be reconstructed.

)The errors listed in Table 3 and plotted in Figure 2 are total errors including systematic errors from the above methods and statistic errors from the final best-fit convergence.

Conclusion on Dark Matter Physics

The current results show that the subhalo of JVAS B1938+666 strong lensing system favours the SIDM or ψDM model with high Bayesian evidences. The best-fit dark matter particle mass for ψDM is found to be mψ = 1.3 + 0.3 − 0.2 × 10−22 eV, which is consistent with results from anomalies in gravitationally lensed images and local dwarf galaxies, but challenged by other observations like the power spectrum of the Lyman-α Forest and the gravitational lensed radio jet. The findings suggest that SIDM or ψDM models might better explain observed flat cores than traditional NFW profiles used in CDM. Furthermore, baryonic feedback remains a viable alternative core formation mechanism, though constraints on stellar mass contribution suggest that for this specific subhalo, the dark matter model is favored. Future high-sensitivity instruments are anticipated to provide stronger constraints on these properties.

How it works

The reconstruction method relies on leveraging purely gravitational effects—as substructures perturb the morphology of lensed arcs—to constrain dark matter properties without relying on baryonic dynamics. This approach allows for a "non-parametric reconstruction of the subhalo mass distribution, in order to be agnostic to specific density profile parametrizations.

As a fastidious and diligent researcher, I have analyzed the provided scientific paper on reconstructing dark matter subhalo profiles in strong lensing systems (JVAS B1938+666) using non-parametric modeling and comparing it against SIDM and Fuzzy Dark Matter (FDM) models.

The improvements suggested for AI systems are centered around leveraging the methodology, data processing pipeline, model comparison framework, and uncertainty quantification presented in this research.

Here are the specific improvements for AI systems:

)

1. AI Systems could be improved by integrating a fully automated, multi-stage gravitational lensing analysis pipeline based on the LensCharm framework (Rüstig et al., 2024). This system would automate the two-step modeling process:

2. Automated Main Lens/Source Fitting (Step 1) using parametric models like Sersic profiles and Pseudo-Isothermal Ellipsoids (PIE), followed by automated subtraction of the lens light.

3. Automated Non-Parametric Subhalo Reconstruction (Step 2) using a Correlated Field Model that employs Gaussian Processes controlled by Wiener processes, allowing for the dynamic reconstruction of subhalo convergence perturbations without pre-assuming a specific density profile parametrization (NFW, SIDM, or FDM).

4. AI could be trained to perform Bayesian model comparison automatically by calculating the Bayes Factor (BF) based on the difference in Bayesian Information Criterion (BIC) between competing models (e.g., NFW vs. SIDM vs. ψDM), as demonstrated by the calculation of a BF = 14.44 in this study.

5. The system should incorporate robust systematic error estimation capabilities, specifically by simulating mock systems with known subhalo profiles (NFW and ψDM) under observational noise conditions to quantify the relative convergence errors and density profile uncertainties, as detailed in Section B.2 and Table 3 of the paper.

6. AI can be specialized in parameter inference using Markov Chain Monte Carlo (MCMC) methods, specifically implementing code like emcee (Foreman-Mackey et al., 2013) to derive posterior distributions for complex dark matter parameters (e.g., particle mass, core radii, scale densities), as shown in Figure 8 and Table 1.

7. The system should be capable of performing model-independent density profile fitting by testing the reconstructed subhalo against analytic profiles from different theoretical frameworks (NFW, SIDM, ψDM) to determine which model provides the best fit based on maximized log-likelihood (e.g., comparing ln(Lmax) values).

8. AI can be enhanced to perform substructure discrimination by analyzing the resulting density and slope ratios (e.g., stellar mass-to-halo mass ratio constraints shown in Figure 3), which helps distinguish between different physical mechanisms like baryonic feedback versus alternative dark matter physics (SIDM/FDM).

The improved AI system can perform the following specific tasks:

1. Perform high-precision, model-independent reconstruction of dark matter density profiles within strong lensing systems, yielding a reconstructed mass distribution (like Table 3) that is agnostic to assumed density profile parametrizations.

2. Quantify and differentiate between competing dark matter models (NFW, SIDM, ψDM) by calculating statistical evidence metrics like the Bayes Factor and BIC to determine which physical model best fits the observed data.

3. Generate comprehensive uncertainty quantification for the reconstructed subhalo mass distribution by incorporating both statistical errors from fitting and systematic errors derived from simulated mock systems (e.g., Figure 7).

4. Infer fundamental dark matter particle properties (like ψDM mass, mψ) by performing parameter estimation via MCMC methods on the posterior distributions of the fitted models.

5. Discriminate between physical mechanisms driving halo cores—specifically determining whether a flat core is better explained by alternative dark matter models (SIDM/ψDM) or baryonic feedback processes based on derived constraints (e.g., M∗/Mhalo limits).

The nature of dark matter remains unknown, motivating the study of fuzzy/wave dark matter (FDM/ ψ DM) and self-interacting dark matter (SIDM) as alternative frameworks to address small-scale discrepancies in halo profiles inferred from observations. This study presents a non-parametric reconstruction of the mass distribution of the previously-found, dark subhalo in the strong-lensing system JVAS B1938+666. Compared with the standard Navarro-Frenk-White (NFW) profile, both SIDM and ψ DM (m ψ=1.32+0.22-0.31 times 10-22, eV) provide significantly better fits to the resulting density profile. Moreover, the SIDM model is favored over ψ DM with a Bayes factor of 14.44. The reconstructed density profile features a characteristic kiloparsec-scale core (r c about 0.5, kpc) with central density ρ c about 2.5 times 10 7, M, kpc-3, exhibiting remarkable consistency with the core-halo mass scaling relations observed in Local Group dwarf spheroidals. These findings offer insights that may help address the core-cusp discrepancy in Λ CDM substructure predictions.