Broken Power Law Mass Model in Lensing

• Broken power law (BPL) mass model is a piecewise analytic density profile with distinct inner and outer slopes defined by a break radius.
• It is used in extragalactic astrophysics and gravitational lensing to model galaxy and dark matter halo profiles with high accuracy.
• The model enables precise predictions of lensing deflection, magnification, and kinematics, reducing biases in cosmological parameter estimates.

The broken power law (BPL) mass model is a class of analytic, piecewise-defined mass distributions widely employed in extragalactic astrophysics and gravitational lensing to represent systems whose density profile exhibits a clear transition between distinct power-law slopes at a finite radius. The flexible parameterization captures deviations from single power-law (SPL) models—including central cores, excesses, and sharp changes in slope—while maintaining analytic tractability for lensing and dynamical calculations. BPL models are central to studies of galaxy mass reconstructions, time-delay cosmography, and the assembly history of supermassive black holes, and have been extensively validated against numerical simulations and observational data.

1. Mathematical Formulation

The core of the BPL model is a density or surface density profile comprised of two power-law branches, joined at a break radius $r_b$ (or projected break radius $\theta_{\rm B}, r_c$), with independent slopes inside and outside:

• 3D Spherical Density (generic):

$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$

where $\rho_c = \rho(r_c)$ is the normalization, $\alpha_c$ and $\alpha$ are the inner and outer slopes respectively, and $r_c$ is the break or core radius (Du et al., 2019, Du et al., 2023, Rui et al., 14 Jan 2026).

• Surface Mass Density (projected, elliptical BPL):

$$\kappa(R) = \begin{cases} \frac{3-\alpha}{2}\left(\frac{b}{R}\right)^{\alpha-1} - \frac{3-\alpha}{\mathcal B(\alpha)} \left(\frac{b}{r_c}\right)^{\alpha-1} \tilde{z}\left[{}_2F_1\left(\tfrac{\alpha_c}{2},1;\tfrac32;\tilde{z}^2\right) - {}_2F_1\left(\tfrac{\alpha}{2},1;\tfrac32;\tilde{z}^2\right)\right], & R < r_c \ \frac{3-\alpha}{2}\left(\frac{b}{R}\right)^{\alpha-1}, & R \geq r_c \end{cases}$$

with $\tilde{z} = \sqrt{1 - R^2/r_c^2}$ and normalization $b$. Elliptical coordinates are introduced via $\theta_{\rm B}, r_c$0 (axis ratio $\theta_{\rm B}, r_c$1) (Rui et al., 14 Jan 2026, Du et al., 2023).

• Generalized Smooth Transitions:

Some variants employ a smooth break using a "sharpness" parameter $\theta_{\rm B}, r_c$2:

$\theta_{\rm B}, r_c$3

where $\theta_{\rm B}, r_c$4, $\theta_{\rm B}, r_c$5 are the inner and outer logarithmic slopes, and $\theta_{\rm B}, r_c$6 controls the sharpness of the transition ($\theta_{\rm B}, r_c$7 yields a perfectly sharp BPL) (Baes et al., 2021).

Analytic formulae for lensing deflection, magnification, and velocity dispersions are expressible in terms of hypergeometric and Beta functions, preserving computational efficiency and flexibility for forward modeling (Du et al., 2019, O'Riordan et al., 2020, Rui et al., 14 Jan 2026).

2. Physical Motivation and Applications

BPL mass models arise wherever the underlying density structure exhibits a genuine transition between two distinct power-law regimes, with applications including:

• Galaxy and Dark Matter Halos:

Lensing galaxies exhibit central mass deficits or surpluses relative to the extrapolation of their outer profiles, motivating the BPL model as an extension of the singular power-law and SIS/SIE families (Du et al., 2019, Du et al., 2023).

• Strong Lensing and Cosmography:

The 2D BPL (2DBPL) model provides key flexibility required for unbiased inference of the Hubble constant $\theta_{\rm B}, r_c$8 from lensed time delays, with additional control over radial degrees of freedom critical for mitigating mass-profile systematics (O'Riordan et al., 2020, Rui et al., 14 Jan 2026).

• Supermassive Black Hole Seed Distributions:

The BPL form naturally arises in self-similar models of black hole mass functions shaped by exponentially growing formation rates and rapid accretion, yielding an intermediate-mass power-law followed by a high-mass cutoff (Basu et al., 2019).

• Kinematic Constraints:

The analytic aperture-weighted velocity dispersion predictions derived via the Jeans equation enable joint fits to lensing and kinematic observables, crucial for breaking mass–anisotropy degeneracies (Du et al., 2019, Du et al., 2023).

3. Key Properties, Parameters, and Special Cases

BPL models are typically parameterized by the following quantities:

Parameter	Description	Typical Range
$\theta_{\rm B}, r_c$9	Inner 3D (or projected) slope	$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$0
$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$1	Outer 3D (or projected) slope	$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$2
$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$3	Break/core radius	$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$4–$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$5 kpc typical
$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$6	Scale normalization (often set by Einstein radius)	$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$7 arcsec (lens scales)
$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$8	Axis ratio (for ellipticity)	$$\rho(r) = \begin{cases} \rho_c\left(\frac{r}{r_c}\right)^{-\alpha_c}, & r \leq r_c \ \rho_c\left(\frac{r}{r_c}\right)^{-\alpha}, & r > r_c \end{cases}$$9
$\rho_c = \rho(r_c)$0	Central mass excess/deficit	derived from $\rho_c = \rho(r_c)$1, $\rho_c = \rho(r_c)$2

Special limiting cases are:

• Singular Isothermal Ellipsoid (SIE): $\rho_c = \rho(r_c)$3, $\rho_c = \rho(r_c)$4;
• Single Power Law (SPL): $\rho_c = \rho(r_c)$5, arbitrary $\rho_c = \rho(r_c)$6 (irrelevant);
• Broken Power Law (BPL): $\rho_c = \rho(r_c)$7, finite $\rho_c = \rho(r_c)$8.

The normalization $\rho_c = \rho(r_c)$9 is typically set by requiring the mean surface density inside the Einstein radius equals the critical density, or by empirical fitting (Du et al., 2023, O'Riordan et al., 2020).

4. Lensing Implementation and Analytical Properties

The BPL framework enables closed-form or semi-analytic expressions for all core lensing quantities:

• Deflection Angles and Potentials:

For elliptically symmetric lenses, all deflection and potential integrals reduce to expressions involving generalized hypergeometric functions; the model admits fast and robust evaluation (O'Riordan et al., 2020, Rui et al., 14 Jan 2026).

• Shear and Magnification:

Shear and magnification are computed via Wirtinger derivatives of the deflection, with explicit expressions for both branches. Magnification critical curves and caustics are readily mapped (Du et al., 2019).

• Numerical Implementation:

BPL models have been implemented in major lensing analysis codes such as Lenstronomy, allowing for joint modeling of image positions, extended arc morphologies, and time delays. Numerical validation against analytic limits confirms accuracy to machine precision (Rui et al., 14 Jan 2026).

• Parameter Degeneracies:

The inner slope $\alpha_c$0 is partially degenerate with ellipticity $\alpha_c$1; the outer slope $\alpha_c$2 trades off with the break radius $\alpha_c$3 in fitting the outer profile. Empirically motivated or light-profile-based priors are used to stabilize fitting (O'Riordan et al., 2020, Du et al., 2023).

5. Performance in Simulations and Observational Applications

• Galaxy-scale Strong Lensing:

Mock lensing studies utilizing Illustris-1 galaxies demonstrate that BPL models with empirically motivated "rigid priors" achieve $\alpha_c$4 median bias in surface mass density within the Einstein radius, outperforming SIE or SPL models (which yield $\alpha_c$5 bias). The Einstein radius itself can be robustly determined from imaging alone, with $\alpha_c$6 median bias (Du et al., 2023).

• Time-delay Cosmography:

BPL modeling consistently reduces systematic biases in time-delay $\alpha_c$7 inference, yielding unbiased results when the true mass profile possesses an inner break. In analyses of quad systems such as WGD 2038-4008, the difference in $\alpha_c$8 derived from BPL and EPL models—that is, $\alpha_c$9 versus $\alpha$0 km s$\alpha$1 Mpc$\alpha$2—is comparable to the Hubble tension, and is dominated by uncertainties in the inner ($\alpha$3) mass profile (Rui et al., 14 Jan 2026). This underscores the relevance of radial modeling freedom and PSF systematics at small radii.

• Joint Lensing-Kinematics Fits:

The BPL approach, combined with analytic velocity dispersion predictions and mild ellipticity-dependent corrections, achieves $\alpha$4 scatter in line-of-sight velocity dispersion recovery, matching or exceeding dedicated kinematic modeling approaches (Du et al., 2019, Du et al., 2023).

6. Theoretical Constraints and Limitations

• Physicality of the 3D Distribution:

The sharpest (idealized) BPL transition ($\alpha$5) cannot represent a fully physically consistent phase-space distribution under isotropic or radially anisotropic orbits: the isotropic distribution function $\alpha$6 is negative just above the binding energy at the break (Baes et al., 2021). Tangential velocity anisotropy could in principle restore positivity, but is generally dynamically disfavored in cosmological halos.

• Regularization and Usage:

In practice, smoothing the radial break with a finite sharpness parameter ($\alpha$7) suffices for a physically reasonable model compatible with realistic dark matter halos. For lensing applications, the projected mass and deflection are always analytic, but caution must be applied in interpreting the implied 3D dynamics (Baes et al., 2021).

• Parameter Choice and Model Selection:

Model selection based solely on image $\alpha$8 can favor over-flexible parameterizations or fail to penalize physically incorrect profiles if lens–light confusion is present. Empirical priors and external constraints (e.g., stellar kinematics, light-profile scaling) are necessary for reliable inference (Du et al., 2023).

7. Broader Context and Variants

• Alternative Broken Power-Law Contexts:

The BPL formalism emerges in the context of black-hole mass functions for DCBH scenarios, where the exponent $\alpha$9 relates directly to the ratio of exponential birth to accretion rates, and the location of the high-mass break to the total accretion phase duration (parametrized by $r_c$0) (Basu et al., 2019).

• Comparison to Double Power Law (DPL) and Zhao Models:

The BPL is a limiting case of the DPL/Zhao family, widely used for dark matter halo modeling. Not all combinations of DPL parameters are physically consistent: sharp transitions (pure BPL) are the most susceptible to dynamical inconsistency (Baes et al., 2021). Softening the transition restores physicality and is standard practice in both lensing and dynamical modeling.

• Implementation in Analysis Pipelines:

BPL models are now implemented in forward-modeling pipelines for imaging, time-delay, and kinematic analysis, with analytic and numerical routines for all quantities of interest (Rui et al., 14 Jan 2026).

---

References: (Basu et al., 2019, Du et al., 2019, O'Riordan et al., 2020, Baes et al., 2021, Du et al., 2023, Rui et al., 14 Jan 2026).