Pith. Sign in

REVIEW

Self-Bound Anisotropic Stars in $\mathcal{R}+\alpha\mathcal{R}^2$ Gravity with a Hayward Exterior: Junction Conditions and NICER Constraints

PaperPDF from arXiv

data files

ReviewThe objections and verdict, as JSON Paper breakdownThe authors' claims, conditions and measurements, as JSON Signed receiptProof this review has not been altered
X Bluesky LinkedIn Reddit Hacker News
Furqan Nazeer, Tooba Feroze

General Relativity and Quantum Cosmology · Submitted 29 Oct 2026 · 85 references

Pith's verdict: Solid algebra and honest limitations, but the headline bounds rest on a Hayward exterior the theory does not solve, and the NICER 'reproduction' is a consistency check.

T0 · AI review

How this review was made: An AI referee (deepseek/deepseek-v4-flash) read the paper on October 30, 2026 and answered a fixed list of questions about it. T0 is the starting level every paper gets. The levels describe how deeply Pith has checked the paper, never the authors.

  1. T0 AI review An AI read the paper and answered a fixed list of questions.
  2. T1 Audit Automatic checks that the references and files the paper relies on exist and match.
  3. T2 Compute The paper's code or calculations re-run with the same results.
  4. T3 Formal statement The claims written out in a formal math language a computer can read.
  5. T4 Proof Those claims proved, and the proof checked by computer.

How Pith checks papers →

· 2 major objections · 3 minor · Ask about this paper · The compiled record · Read the full review →
Challenge this review Comment (0) Read on arXiv ↗
Cite ▾
arXiv ID2610.39916 (version 1)
Pith IDR4XWGSN6, Pith's permanent record of this paper and its review
CitationFurqan Nazeer, Tooba Feroze. "Self-Bound Anisotropic Stars in $\mathcal{R}+\alpha\mathcal{R}^2$ Gravity with a Hayward Exterior: Junction Conditions and NICER Constraints." arXiv:2610.39916 (2026). https://pith.science/paper/R4XWGSN6
Cite this reviewBibTeX for Pith's review of this paper
Claim paper

Why Pith says this

Highlight any text to comment on it.

Referee report

How Pith checks papers → · Full audit of this paper →

The paper

The paper claims

In R + αR² gravity, matching a star to a Hayward exterior with the theory's own junction conditions fixes every interior constant in closed form and leaves a surface stress that bounds the two parameters.

The reading

Most analytic star models in modified gravity borrow the matching rules of general relativity and demand that the radial pressure vanish at the stellar surface. This paper argues that in quadratic $F(\mathcal{R}) = \mathcal{R} + \alpha\mathcal{R}^2$ gravity that is the wrong condition: the theory's junction conditions instead require continuity of the trace of the extrinsic curvature, $[K^\mu{}_\mu] = 0$, and imposing this fixes every interior constant in closed form from the mass, the radius, and the Hayward deformation length alone. The price is a small negative radial stress at the surface and a thin shell where the Ricci scalar jumps; demanding that the stress stay below one percent of the central pressure bounds the two free parameters, $\alpha \lesssim 0.2$ km$^2$ and $\ell \lesssim 0.3M$. With those values, a single surface energy density $\rho_s = (1.37 \pm 0.16)\rho_0$ reproduces four NICER pulsars within 1.2$\sigma$ and predicts a 1.4-solar-mass radius of 11.5 km and a maximum mass near $2.9\,M_\odot$; the authors note that per-pulsar fits cannot by themselves confirm the model, because the constants can always adjust to a measured $(M,R)$, so the joint fit is the real test. If correct, the paper gives a parameter-free route from a chosen interior ansatz and exterior metric to testable mass–radius predictions, and it exposes a mismatch in earlier models that imposed $p_r(R) = 0$.

What carries the argument

The load-bearing mechanism is the pair of junction conditions for quadratic $F(\mathcal{R})$ gravity, $[h_{ab}] = 0$ and $[K^\mu{}_\mu] = 0$, which replace the general-relativity condition $p_r(R) = 0$. Because the angular component of the extrinsic curvature is automatically continuous once the induced metrics agree, the trace condition reduces to $[A'/A] = 0$, and this fixes the integration constant $b$ in closed form from the Hayward profile $C(r) = 1 - 2Mr^2/(r^3 + 2\ell^2 M)$; the other two constants then follow from metric continuity. The Hayward metric is a controlled one-parameter deformation of Schwarzschild whose correction to surface observables enters only at order $(\ell/M)^2(M/R)^3$, which is why the results survive in the Schwarzschild limit. The quantity that carries the physics is the closed-form surface stress of Eq. (23), $p_r(R)$: its smallness sets the bounds on $\alpha$ and $\ell$, while the shell and double-layer stress tensor of Eq. (24) describes the surface layer that the Ricci-scalar jump produces.

Extended reading notes

Core claim

The central claim is that the matching condition itself, not any equation of state, determines the model. For a Durgapal–Bannerji interior of embedding class one (the Karmarkar condition) matched to a Hayward exterior through the mandatory conditions $[h_{ab}] = 0$ and $[K^\mu{}_\mu] = 0$, the constants $a$, $b$, $c$ are fixed in closed form by $M$, $R$ and the Hayward length $\ell$ alone, with $b = f_2(R) C'(R)/(6R\sqrt{C(R)})$; none of them depends on $\alpha$. Because $b$ is fixed by the exterior, the radial pressure no longer vanishes at the surface: $p_r(R)$ is small and negative, vanishing instead 17–67 m below the surface, and the jump of the Ricci scalar creates a thin shell with negative surface energy density for $\alpha > 0$, a gravitational double layer, and a shell mass of 0.1–0.2% of the stellar mass. Requiring $|p_r(R)| \le 0.01\,p_c$ gives $\alpha \lesssim 0.2$ km$^2$ and $\ell \lesssim 0.3M$, roughly ten times tighter than the bulk admissibility tests allow. The observational consequence is that one shared surface density, $\rho_s = (1.37 \pm 0.16)\rho_0$, fits PSR J0740+6620, J0437−4715, J0030+0451 and J0614−3329 within 1.2$\sigma$, yielding $R_{1.4} = 11.5$ km, a maximum mass of $2.90^{+0.20}_{-0.16}\,M_\odot$, and a radius that grows with mass between 1 and 2 $M_\odot$.

Load-bearing premise

The load-bearing premise is that the Hayward metric can serve as the star's exterior even though it is not a solution of the $R + \alpha R^2$ field equations the interior is built in — it comes from general relativity with a different matter source — so the closed-form constants, the surface stress, and both parameter bounds are computed against an exterior the theory itself does not justify; the authors state this limitation and defer a self-consistent scalaron exterior to future work.

Editorial extensions

If this is right

  • Any stellar model in quadratic $F(\mathcal{R})$ gravity that fixes its last constant by $p_r(R) = 0$ at a non-Schwarzschild exterior violates a mandatory junction condition; the authors recommend re-checking such models against $[K^\mu{}_\mu] = 0$.
  • The surface stress, not the bulk admissibility tests, is the tightest constraint: it limits $\alpha$ to about 0.2 km$^2$ and the Hayward length to $\ell \lesssim 0.3M$, roughly ten times stricter than the values at which the interior matter profiles first fail.
  • Self-bound matter with one surface density, $\rho_s = (1.37 \pm 0.16)\rho_0$, describes four NICER pulsars within 1.2$\sigma$ and predicts $R_{1.4} = 11.5$ km and a maximum mass of $2.90^{+0.20}_{-0.16}\,M_\odot$ (3.31 at the causal limit).
  • The radius increases with mass between 1 and 2 $M_\odot$, so a prior-independent radius for PSR J1231−1411 near 12.6 km, or a firm neutron-star maximum-mass limit near $2.3\,M_\odot$, would each count against the model.
  • Because the Hayward length enters surface observables only at order $(\ell/M)^2(M/R)^3$, the NICER constraints remain valid for a Schwarzschild exterior, and $\ell$ is constrained by the junction conditions rather than by the radii.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • Read through the bag-model relation $\rho_s = 4B$ that the paper quotes, the fitted surface density corresponds to an effective bag constant of about 50–60 MeV/fm$^3$, comparable to standard strange-quark-star values; the paper does not draw this connection, but it reinforces the self-bound interpretation.
  • The bounds $\alpha \lesssim 0.2$ km$^2$ and $\ell \lesssim 0.3M$ inherit the paper's chosen 1% tolerance on the surface stress — relaxing it to 5% loosens the $\alpha$ bound to about 1.1 km$^2$ — so the quoted limits are modeling conventions rather than sharp physical thresholds until an independent criterion fixes the acceptable shell strength.
  • The same junction-condition logic extends to other regular exteriors: the paper's estimate for the Bardeen metric gives a correction of order $(g/R)^2$, suggesting that any one-parameter regular deformation of Schwarzschild should be admitted only after its own surface-stress bound is imposed.
  • If a self-consistent scalaron exterior were constructed, the Ricci scalar would no longer jump, so the thin shell and double layer would disappear; because the shell's negative energy density is what violates the weak energy condition, such an exterior would likely shift the parameter bounds as well, a route the authors flag as the natural next step.
Original abstract

Analytic compact-star models in modified gravity usually borrow the junction conditions of general relativity and fix their constants with radii derived for a particular stellar model. Here we avoid both. We study $\mathcal{F}(\mathcal{R})=\mathcal{R}+\alpha\mathcal{R}^{2}$ gravity with a Durgapal-Bannerji interior of embedding class one (Karmarkar condition) and a Hayward exterior. Without Birkhoff's theorem, the exterior is an assumption, and we measure its effect with the Hayward metric, which deforms Schwarzschild by a length $\ell$. The junction conditions of quadratic $\mathcal{F}(\mathcal{R})$ gravity fix all constants in closed form in terms of $M$, $R$ and $\ell$. The condition $p_r(R)=0$ is replaced by continuity of the trace of the extrinsic curvature. The Ricci scalar jumps at the surface, so there is a thin shell with negative surface energy density for $\alpha>0$, a gravitational double layer, and a small radial stress. For this stress to be below 1% of the central pressure we need $\alpha\lesssim0.2$ km$^{2}$ and $\ell\lesssim0.3M$. The interior behaves as a self-bound star. A single surface energy density, $\rho_s=(1.37\pm0.16)\rho_0$, reproduces PSR J0740+6620, J0437-4715, J0030+0451 and J0614-3329, and for $\alpha=0.1$ km$^{2}$ the model predicts $R_{1.4}=11.5$ km and a maximum mass of $2.9 M_\odot$. Because $\ell$ enters surface observables only at order $(\ell/M)^{2}(M/R)^{3}$, the results carry over to a Schwarzschild exterior. Calibrating instead with strange-star radii of eclipsing X-ray pulsars gives $R_{1.4}\approx9$ km, which current data rule out.

The compiled record

The paper's own sentences, arranged, kept beside the review and separate from it. Each label is drawn from a closed vocabulary; the number beside it is the probability the classifier gave that label. Numbers in the tables are the authors', checked by code. No sentence in this record was written by a model.

Compiled by classifier-1.13.0 via TypeSafe API, question bank abstract-p1-20260923, 2026-10-30: 1 request, 8,697 input tokens.

What the abstract itself states, as probabilities: claims first 0.06 · states limitation 0.91 · human participants 0.01 · controlled comparison 0.14 · randomized assignment 0.01 · states quantitative result 0.98 · compares with prior results 0.84 · mentions code or data availability 0.02.

The abstract, sentence by sentence

arXiv listing abstract v1. Results and the main claim are underlined; .

Analytic compact-star models in modified gravity usually borrow the junction conditions of general relativity and fix their constants with radii derived for a particular stellar model. Here we avoid both. We study $\mathcal{F}(\mathcal{R})=\mathcal{R}+\alpha\mathcal{R}^{2}$ gravity with a Durgapal-Bannerji interior of embedding class one (Karmarkar condition) and a Hayward exterior. Without Birkhoff's theorem, the exterior is an assumption, and we measure its effect with the Hayward metric, which deforms Schwarzschild by a length $\ell$. The junction conditions of quadratic $\mathcal{F}(\mathcal{R})$ gravity fix all constants in closed form in terms of $M$, $R$ and $\ell$. The condition $p_r(R)=0$ is replaced by continuity of the trace of the extrinsic curvature. The Ricci scalar jumps at the surface, so there is a thin shell with negative surface energy density for $\alpha>0$, a gravitational double layer, and a small radial stress. For this stress to be below 1% of the central pressure we need $\alpha\lesssim0.2$ km$^{2}$ and $\ell\lesssim0.3M$. The interior behaves as a self-bound star. A single surface energy density, $\rho_s=(1.37\pm0.16)\rho_0$, reproduces PSR J0740+6620, J0437-4715, J0030+0451 and J0614-3329, and for $\alpha=0.1$ km$^{2}$ the model predicts $R_{1.4}=11.5$ km and a maximum mass of $2.9 M_\odot$. Because $\ell$ enters surface observables only at order $(\ell/M)^{2}(M/R)^{3}$, the results carry over to a Schwarzschild exterior. Calibrating instead with strange-star radii of eclipsing X-ray pulsars gives $R_{1.4}\approx9$ km, which current data rule out.

background contribution method result interpretation main claim

What the classifier reports for each sentence. The number is the probability it gave the label.
#sentencerolethe authors' own workstates the main conclusionstrength
s0 Analytic compact-star models in modified gravity usually borrow the junction conditions of general relativity and fix their constants with radii derived for a particular stellar model. Here we avoid both. background 0.53 0.40 0.06 asserted 0.99
s1 We study $\mathcal{F}(\mathcal{R})=\mathcal{R}+\alpha\mathcal{R}^{2}$ gravity with a Durgapal-Bannerji interior of embedding class one (Karmarkar condition) and a Hayward exterior. contribution 0.68 0.79 0.04 asserted 0.98
s2 Without Birkhoff's theorem, the exterior is an assumption, and we measure its effect with the Hayward metric, which deforms Schwarzschild by a length $\ell$. method 0.68 0.71 0.04 asserted 1.00
s3 The junction conditions of quadratic $\mathcal{F}(\mathcal{R})$ gravity fix all constants in closed form in terms of $M$, $R$ and $\ell$. result 0.44 0.91 0.19 asserted 1.00
s4 The condition $p_r(R)=0$ is replaced by continuity of the trace of the extrinsic curvature. method 0.72 0.76 0.08 asserted 1.00
s5 The Ricci scalar jumps at the surface, so there is a thin shell with negative surface energy density for $\alpha>0$, a gravitational double layer, and a small radial stress. result 0.99 0.89 0.12 asserted 1.00
s6 For this stress to be below 1% of the central pressure we need $\alpha\lesssim0.2$ km$^{2}$ and $\ell\lesssim0.3M$. result 0.80 0.89 0.12 asserted 0.95
s7 The interior behaves as a self-bound star. result 0.67 0.84 0.19 asserted 1.00
s8 A single surface energy density, $\rho_s=(1.37\pm0.16)\rho_0$, reproduces PSR J0740+6620, J0437-4715, J0030+0451 and J0614-3329, and for $\alpha=0.1$ km$^{2}$ the model predicts $R_{1.4}=11.5$ km and a maximum mass of $2.9 M_\odot$. result 1.00 0.89 0.37 asserted 0.99
s9 Because $\ell$ enters surface observables only at order $(\ell/M)^{2}(M/R)^{3}$, the results carry over to a Schwarzschild exterior. interpretation 0.80 0.82 0.13 asserted 1.00
s10 Calibrating instead with strange-star radii of eclipsing X-ray pulsars gives $R_{1.4}\approx9$ km, which current data rule out. result 0.81 0.76 0.08 asserted 0.99

What kind of paper the classifier reports this is: theoretical result 1.00 · evidence from proof or derivation 0.98.

What was examined, and what was not

sources read
  • arXiv listing abstract v1
not examined
  • full text
  • figures and tables
  • references and citing sentences
  • relations to other papers
levels run
  • P abstract
not run
  • S sections
  • T body sentences
  • C citations
  • R cross-paper relations

Every quoted sentence on this record is the authors'. Every label is from Pith's closed vocabulary, with the probability shown beside it. Every number in the measurements table is the authors', checked by code. Record as JSON

Figures
Figure 1. Metric potentials (a) A2 (r) and (b) B2 (r) as functions of r/R for the pulsars of
Figure 2. (a) Energy density, (b) radial pressure and (c) tangential pressure as functions of r/R. The radial pressure vanishes at r0, 17–67 m below the surface (
Figure 3. (a) Anisotropy ∆ = pt − pr and (b) normalized trace anomaly ∆tr = 1/3 − (pr + 2pt)/(3ρ) as functions of r/R; the dashed line marks the conformal limit ∆tr = 0. than one everywhere inside the star. The two ratios, wr = pr/ρ and wt = pt/ρ, are plotted in
Figure 4. Equation-of-state ratios (a) wr = pr/ρ and (b) wt = pt/ρ as functions of r/R. Both stay below 0.16, so the dominant energy condition is satisfied with a wide margin. The surface redshift is fixed by the exterior metric, zs = C(R) −1/2 − 1 = B(R) − 1. (26) For the configurations of
Figure 5. (a) Radial and (b) tangential sound speeds and (c) the cracking indicator v 2 r −v 2 t as functions of r/R. Dashed lines mark the conformal value 1/3 in (a) and (b); in (c) the indicator never leaves the interval [0, 1]. Another important quantity, the adiabatic index Γi , is associated with the stability and stiffness of matter distribution. Following Heintzmann and Hillebrandt [66], it is defined as Γi = ρ + pi pi dpi dρ . (28) For relativistic stars the Newtonian bound Γ > 4/3 is not sufficient, so we use Eq. (25) instead. As
Figure 6. Adiabatic indices (a) Γr and (b) Γt as functions of r/R. The gray band covers the values of Γcrit = 4/3 + (19/21)M/R for the five pulsars. • SEC: ρ + pi ≥ 0, ρ + pr + 2pt ≥ 0, Since ρ > 0 and pt > 0 everywhere, pr ≥ 0 below r0, and |pr| ≤ 0.0085 pc in the thin layer r0 < r ≤ R, the null, weak and strong energy conditions are satisfied throughout the star. The dominant energy condition follows from wr, wt ≤ 0.16 (
Figure 7. Force densities Fg, Fh and Fa of Eq. (29) and their sum for PSR J0740+6620 (α = 0.1 km2 , β = 0.1). reaches 0.34. Two-solar-mass stars require such a stiffening of the equation of state [70]. For this reason we interpret these objects as self-bound compact stars rather than neutron stars with a crust. It has recently been argued that the small radius of PSR J0614−3329 favors precisely this class of stars [71]. In Sec. 5.4 the surface density ρs becomes the parameter shared by all pulsars. 5 Comparison with NICER and Radio Observations NICER measures the masses and radii of…
Figure 8. Allowed values of α and β = ℓ/M. For each pulsar, the radial stress at the surface satisfies |pr(R)| ≤ 0.01 pc in the region between its curve and the origin. The star shows the values α = β = 0.1 used in this work. The bulk tests of Sec. 4 alone would allow α up to 2–3 km2 , which is outside the plotted range. such stars the surface density is set by the properties of the matter itself (ρs = 4B in the MIT bag model with bag constant B), so it is the same for every star. This suggests using the surface energy density ρs ≡ ρ(r = R) (30) as a parameter common to all the pulsars. Given…
Figure 9. (a) Mass–radius relation of the model for a fixed surface energy density ρs = (1.37 ± 0.16) ρ0 (green line and band; α = 0.1 km2 , β = 0.1). Solid: configurations that pass all tests of Sec. 4. Dashed: causal configurations that fail the no-cracking test (high mass) or the surface-stress test of Sec. 5.3 (low mass); the dot marks the causal limit. Orange: calibration with the masses of eclipsing X-ray pulsars [80] and their strange-star radii [39] (open squares; ρs = 2.68 ρ0). Black: NICER measurements (68%; the J1231−1411 radius depends on its prior). Gray bands: radio masses of…

Figures reproduced from arXiv:2610.39916 by the authors.

References

Discussion

Disagree with the review, or know something it missed? Leave a comment, highlight any passage above to comment on it, or challenge this review formally.