Pith. sign in

REVIEW 4 major objections 6 minor 18 references

Charged Dirac stars are gravitationally bound only for q<m; super-critical solutions always remain unbound.

Reviewed by Pith at T0; open to challenge. T0 means a machine referee read the full paper against a public rubric. the ladder, T0–T4 →

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗

T0 review · deepseek-v4-flash

2026-08-01 17:09 UTC pith:DQZSD23A

load-bearing objection A clean 3+1 derivation and numerical survey of charged Dirac stars; the central q<m bound is probably right but phrased more universally than the evidence supports. the 4 major comments →

arxiv 2607.26374 v1 pith:DQZSD23A submitted 2026-07-29 gr-qc

Charged Dirac stars

Maribel Hern\'andez M\'arquez , Miguel Alcubierre This is my paper
classification gr-qc MSC 83C0583C2283C60 PACS 04.20.Ex04.25.Dm95.30.Sf
keywords Einstein–Dirac–MaxwellDirac stars3+1 formalismself-gravitating spinor fieldsbinding energysuper-critical solutionscompactnessboson-star universality
verification ladder T0 review T1 audit T2 compute T3 formal T4 reserved

The pith

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

This paper builds stationary, spherically symmetric solutions of the Einstein–Dirac–Maxwell system — two fermions in a singlet spinor state — using the 3+1 formalism and a numerical shooting method. It claims that gravitationally bound configurations (negative binding energy) exist only when the charge parameter q is smaller than the fermion mass m, and that even when 'super-critical' solutions with q>m exist in general relativity, they are always gravitationally unbound. The mass–frequency curves reproduce the spiral structure familiar from scalar and vector boson stars, indicating that spin does not change the qualitative landscape of charged self-gravitating matter. Some bound solutions reach compactness C~0.1–0.2, comparable to neutron stars.

Core claim

The authors derive the full 3+1 Einstein–Dirac–Maxwell equations, specialize to spherical symmetry and a static metric, and impose a harmonic time dependence on the spinor amplitudes to find stationary solutions. They show that, for each charge q, there is a family of regular solutions parameterized by the central amplitude f0; integrating the ODE system (100)–(105) with shooting on the frequency ω yields exponentially decaying spinors at infinity. The central discovery is that the binding energy EB = M_RN − m Q/q is negative — the configuration is gravitationally bound — only for q<m, and that in every computed solution the total charge-to-mass ratio satisfies Q/M_RN < 1, even when q>m. Thi

What carries the argument

The engine is the six-function ODE system (100)–(105) for the lapse α, radial metric a, Dirac amplitudes f and g, electrostatic potential ϑ=αφ, and radial electric field E, obtained from the 3+1 decomposition of the Einstein–Dirac–Maxwell action with a spherically symmetric, polar-areal metric and a two-spinor singlet ansatz. The harmonic ansatz F=f e^{-iωt}, G=i g e^{-iωt} makes the matter static while allowing the spinors to carry the frequency ω; the requirement m²>ω² for exponential decay turns ω into an eigenvalue found by shooting. The binding energy EB=M_RN−m N (with N=Q/q the effective particle number) is the criterion separating bound (EB<0) from unbound configurations, and the ineq

Load-bearing premise

The universal statement that bound configurations exist only for q<m rests on a finite numerical scan (f0∈[0.01,1], q up to 1.052) plus a cited theorem whose hypotheses are not restated, so a missing branch near q=1, f0→0 — where mass and radius grow without bound — could break the claim.

What would settle it

A targeted search in the region q∈[0.99,1.052], f0∈[0.001,0.05] with an outer boundary at r≥2000, looking for any stationary solution with negative binding energy (equivalently Q/M_RN>q), would falsify the bound claim if found; alternatively, verifying whether the cited no-bound theorem covers the harmonic-time ansatz rather than only strictly static spinors would settle whether the universality is proven.

Watch this falsifier. Get emailed when new claim-graph text bears on it.

If this is right

  • The mass–frequency spiral and the q<m bound hold for fermionic matter just as for scalar and vector fields, supporting a spin-independent universality of charged boson-star-like equilibria.
  • Any astrophysical search for compact charged objects should expect gravitationally unbound super-critical branches to be unstable transients, not equilibrium end states.
  • Bound charged Dirac stars with compactness C≈0.1–0.2 are viable neutron-star-like alternatives in mass–radius observations, though their charges would make them electromagnetic emitters.
  • The q=1 branch, with mass and radius both diverging as f0→0, suggests a critical limit worth studying as a possible extremal or horizonless ultracompact object.
  • The bound criterion EB<0 is equivalent to Q/M_RN > q/m; measuring the charge-to-mass ratio of a stationary configuration directly diagnoses whether it is gravitationally bound.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

  • The universality across spins hints that the q<m bound may follow from an energy condition or from the dominant balance of Coulomb and gravitational forces in the Newtonian limit, independent of the spinor structure; a proof from first principles would be a natural next step.
  • A testable extension: repeat the analysis for self-interacting potentials or for a charged scalar with different gauge couplings; if the spiral shape and bound threshold persist, the phenomenon is a gauge-field effect rather than a spin effect.
  • The apparent f0→0 divergence at q=1 may indicate that the limit approaches an extremal Reissner–Nordström solution or a charged naked singularity; a high-resolution study of that limit could connect Dirac stars to black-hole critical phenomena.
  • Because N=Q/q is treated as an effective particle number, a quantum treatment with Pauli exclusion would cut off the two-fermion sector at N=2; classical solutions with N≫2 are therefore not physical multi-fermion stars, and the paper only claims classical-level validity.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit.

Referee Report

4 major / 6 minor

Summary. The paper develops the 3+1 Einstein-Dirac-Maxwell (EDM) system for static, spherically symmetric configurations of two spin-1/2 fields in an opposite-spin singlet state, and reduces it to an ODE system (Eqs. 100–105) for the radial spinor amplitudes f,g, the electric field E, the electric potential ϑ (or ϕ), the lapse α, and the radial metric a. Stationary solutions are found by a shooting method for m=1, q∈[0,1.052], and f0∈[0.01,1]. The authors report families of solutions with a mass–frequency spiral for q<1, loss of the spiral and apparently unbounded mass for q=1 as f0→0, short super-critical branches for 1<q≤1.052, always-positive binding energy for q≥1, and Q/M_RN<1 throughout. They conclude that gravitationally bound charged Dirac stars exist only for q<m, that super-critical solutions are always gravitationally unbound, and that compactness values comparable to neutron stars occur for bound configurations.

Significance. If the results hold, this is a useful contribution: it extends the known charged boson-star and Proca-star phenomenology to spin-1/2 fields, and it provides a 3+1 formulation that can serve as a starting point for dynamical evolutions. The paper is clearly organized and contains explicit derivations of the field equations, boundary conditions, global quantities, and a parameter table of representative solutions. Its main value is the claimed universality across spins s=0,1/2,1. However, the central nonexistence claim is currently backed by a finite numerical scan plus a citation to a theorem whose hypotheses are not stated, and there are no convergence or truncation-error tests for the unusual q=1 branch. The manuscript would be publishable after the numerical evidence is made reproducible or the claims are explicitly qualified.

major comments (4)
  1. [Abstract, Sec. IX] The universal claim that gravitationally bound configurations exist only for q<m is not established by the finite scan (q≤1.052, f0∈[0.01,1]) alone, and the cited theorem [18] is never stated. The solutions have harmonic time dependence e^{-iωt} (Eq. 127), and it is not clear that [18], which the authors do not summarize, applies to this stationary ansatz rather than to strictly static spinors. Please restate the theorem, verify that its hypotheses cover the present ansatz, or explicitly qualify the claim to the scanned parameter range.
  2. [Sec. VII (q=1 branch)] The q=1, f0→0 branch is central to the claim that the mass grows without bound and that no bound states exist at q=1. The outer boundary is moved from r=300 to r=1000 as f0 decreases, yet no convergence test with respect to Δr or r_max is reported. Since M_RN is extracted from a(r) at the outer boundary through Eq. (124), an insufficiently distant outer boundary could artificially produce an apparent divergence. Please provide systematic convergence data or an analytic estimate for this limit.
  3. [Sec. VI, Eqs. (120) and (124)] The paper states that the integrated mass (120) and the Reissner–Nordström mass (124) coincide asymptotically, and that this was 'verified', but no quantitative comparison is provided. The binding-energy sign in Eq. (125) is sensitive to M_RN, so the claimed E_B<0 intervals depend on the accuracy of the extraction. Please report the difference between the two mass definitions at the outer boundaries used for each family, or justify why Eq. (124) is accurate at those radii.
  4. [Sec. VII and Sec. IX] The statement that super-critical solutions are 'always gravitationally unbound' is operationalized by E_B>0, but the paper does not show that E_B>0 implies dynamical instability for these stationary solutions. If the claim is only about the binding-energy criterion, it should be stated as such; if it is a stability claim, it needs dynamical perturbation analysis or a clear reference to a theorem covering this case.
minor comments (6)
  1. [Sec. V, Eq. (116)] The notation ∂_r E = 2qf_0^2/3 should be written as ∂_r E(0)=2qf_0^2/3 to avoid ambiguity.
  2. [Table II] The header contains a duplicated Q: '(f0,ω,Q,M_RN,Q,N,R99,C)' should be '(f0,ω,M_RN,Q,N,R99,C)'.
  3. [Sec. VII B] R99 is defined as 99% of total charge for q>0, but for q=0 the paper switches to 99% of total mass. This mixed definition should be stated explicitly in the main text and figure/table captions, since compactness comparisons across q=0 and q>0 rely on it.
  4. [Sec. V] The gauge transformation (119) is described only briefly. Please state more explicitly how ϑ_∞ and the sign in ω→ω+q∂_tθ enter the final physical frequency and potential, since this rescaling and gauge step affect all reported values of ω and ϕ.
  5. [Global presentation] The paper would benefit from a data-availability or reproducibility section: no code, data files, or grid-convergence tables are provided. This is not required for all journals, but would strongly increase confidence in the numerical results.
  6. [Throughout] There are several typographical and formatting errors (e.g., 'adjuct', 'con gurations', inconsistent spacing in equations and references). A careful proofread is needed before final submission.

Circularity Check

0 steps flagged

No significant circularity: central predictions are computed outputs; self-citations are methodological rather than load-bearing.

full rationale

The derivation chain is self-contained at the level of the claimed predictions. Starting from the EDM action (Eqs. 1-13), the paper derives the 3+1 equations, specializes to spherical symmetry, imposes the stationary harmonic ansatz F=f(r)e^{-i\omega t}, G=i g(r)e^{-i\omega t} (Eq. 81), and solves the resulting ODE system (100)-(105) by shooting on \omega so that the spinor fields decay at infinity. The global quantities M_RN, Q and E_B are then computed from the solutions (Eqs. 122, 124, 125). The binding-energy sign is an output: E_B = M_RN[1-(m/q)(Q/M_RN)] (Eq. 126), and the statement that q\ge m solutions are unbound follows from the numerically computed Q<M_RN (Fig. 7), not from imposing E_B>0. The mass-frequency spiral and compactness are also outputs of the shooting families. There is no fitted parameter that is relabeled as a prediction. The self-citations [15], [5], [16] are used for the 3+1 Dirac formalism, gauge/rescaling conventions, and comparison with boson/Proca stars; the charged Dirac extension is derived and solved in this paper, so the central result does not reduce to those citations. Citation [18] is an external consistency check, but its hypotheses are not restated; this is a rigor concern, not circularity. The main caveats are numerical: the scan covers only f0\in[0.01,1] and q\le 1.052 with \Delta r=0.01, and for q=1 the paper itself notes that 'the total mass seems to increase without bound as we decrease the value of f0' (Sec. VII C), with outer boundaries placed at finite radii. These limit the strength of the universal 'only q<m' claim but do not make it circular.

Axiom & Free-Parameter Ledger

3 free parameters · 7 axioms · 0 invented entities

No new particles, forces, or conserved quantities are invented. The central assumptions are the standard EDM model, the spherical two-spinor ansatz, the harmonic-time form of the Dirac field, and the numerical convergence of the shooting method. The effective particle number N=Q/q is an interpretation of the classical Noether charge, not a new entity.

free parameters (3)
  • q (charge parameter) = scanned values 0, 0.1, 0.3, 0.6, 0.9, 0.95, 0.99, 1, 1.01-1.052
    Electric coupling constant of the Dirac field; varied by hand to generate solution families. Physical input, not fitted to external data.
  • m (Dirac field mass) = 1 (rescalable)
    The authors set m=1 and note the system can be rescaled to arbitrary m; a physical parameter, not fitted.
  • f0 (central spinor amplitude) = scanned over [0.01,1]
    Parametrizes the shooting family; chosen by hand at the origin and used as continuation parameter.
axioms (7)
  • domain assumption Minimal-coupling Einstein-Dirac-Maxwell action and field equations (Eq. 1), with G=c=hbar=1 and Lorentzian signature.
    The physical model the paper solves; not derived within the paper.
  • domain assumption Static, spherically symmetric spacetime with areal radial coordinate, zero shift, and Fermi-Walker propagated triad (Q_IJ=0).
    Used to reduce the 3+1 system to ODEs in Secs. III-IV.
  • domain assumption Two-spinor singlet ansatz with opposite half-integer azimuthal modes restores spherical symmetry of the stress-energy tensor.
    Required for strictly spherically symmetric spinor configurations; standard in Dirac-star literature.
  • domain assumption Harmonic time dependence F=f e^{-i omega t}, G=i g e^{-i omega t} with real omega, and decaying boundary condition requiring m^2>omega^2.
    Defines stationary rather than static spinor fields; Sec. IV-V.
  • domain assumption Regularity expansions at the origin: f even, g odd, a=1+O(r^2), alpha=alpha_0+O(r^2), E=O(r).
    Boundary conditions (115)-(116) used to start the shooting integration.
  • domain assumption Cited theorem [18] that no gravitationally bound q>=1 configurations exist.
    Used in Sec. IX to convert a finite numerical scan into a universal claim; hypotheses not restated.
  • ad hoc to paper Numerical shooting on grids with Delta r=0.01 and outer boundaries up to r=1000 gives the true exponentially decaying ground-state solutions over the scanned ranges.
    No convergence or truncation-error study is presented; this is the weakest support for the universal claims.

pith-pipeline@v1.3.0-daily-deepseek · 20588 in / 19915 out tokens · 227890 ms · 2026-08-01T17:09:57.648496+00:00 · methodology

0 comments
read the original abstract

In this work we solve the coupled Einstein-Dirac-Maxwell (EDM) system for static spherically symmetric configurations of two fermions in a singlet spinor state within the $3+1$ formalism of general relativity. We find different families of stationary self-gravitating solutions for the Dirac field through a numerical shooting method for different values of the electric charge parameter q. Furthermore, we investigate the effect of the charge q on the binding energy, mass, radius, and compactness of the solutions. We show that gravitationally bound configurations exist only for $q<m$, with m the mass of the Dirac field, and that the mass frequency relation exhibits the characteristic spiral structure previously found for bosonic fields of spin $s=0$ and $s=1$. We are able to show that some of these gravitationally bound configurations have a compactness comparable to that of neutron stars. With these results, we conclude that at least at the classical level, self-gravitating fields with different spins $s=0,1/2,1$ share some common characteristics when they are coupled to gravity. As has been previously shown for the case of bosonic stars, we also find some super-critical solutions with q slightly larger than m. Such super-critical solutions correspond to configurations such that in the Newtonian regime the Coulomb repulsion overcomes the gravitational attraction, and as such they would not be expected to exist. Nevertheless, even if they do exist in the general relativistic case for a limited range of values of $q>m$, we find that they are always gravitationally unbound.

Figures

Figures reproduced from arXiv: 2607.26374 by Maribel Hern\'andez M\'arquez, Miguel Alcubierre.

Figure 1
Figure 1. Figure 1: Frequency ω as a function of f0 for dierent values of q. For subcritical solutions with q ≤ 1 the lowest value of the frequency is ω = 07327 for q = 0 and f0 = 059. We do nd super-critical solutions for 1 < q ≤ 1052, but only for a limited range of values of f0. outer boundary is xed at r = 300 for f0  [005, 1], at r = 400 for f0 = (003, 004), at r = 500 for f0 = 002, and at r = 1000 for f0 = 0… view at source ↗
Figure 2
Figure 2. Figure 2: Eective radius R99 as a function of f0 for dierent values of q. For q ≤ 1, the radius becomes larger as f0 → 0, but decreases with increasing f0 until reaching a minimum value, and then increases again as f0 becomes even larger. Also, as the charge q increases the values of R99 increase for the same value of f0. of f0, which is to be expected since the Coulomb repulsion works against gravity. In particul… view at source ↗
Figure 3
Figure 3. Figure 3: Left panel: Total mass MRN with respect to ω for dierent values of q. For q ≤ 095 the solutions describe a clear spiral, while for q = 099 the spiral shrinks and moves to larger values of ω. For q = 1 the spiral structure is lost, and the mass appears to grow without bound as ω → 1 and f0 → 0. The curves for q = 0 and q = 01 are nearly indistinguishable. Right panel: We show only the super-critical cas… view at source ↗
Figure 4
Figure 4. Figure 4: Total mass MRN with respect to f0 for dierent values of q. For the cases with q < 1 the maximum value of the mass increases and moves to smaller values of f0. On the other hand, for the special case with q = 1 the mass seems to increase without bound as f0 approaches 0. a limit of C ∼ 0302. We can compare these values with the compactness of a black hole C = 05, and with that of a neutron star C ∼ 03 [… view at source ↗
Figure 5
Figure 5. Figure 5: Left panel: Compactness C with respect to f0 for dierent values of q. For q < 1 the compactness has a clear maximum that becomes larger and moves to the left as q increases. On the other hand, for q = 1 the maximum compactness is reached as f0 → 0, with a value of C ∼ 0302. Right panel: We show only the cases with q ≥ 1 for clarity. increases. In Table I we show the values of (f0, ω, MRN , Q, N, R99, C) … view at source ↗
Figure 6
Figure 6. Figure 6: Upper left panel: Binding energy EB as a function of f0 for the dierent values of q. For q < 1 there exist intervals of f0 for which EB < 0, and the range of these intervals decreases as the charge increases. However, for q ≥ 1 we always have EB > 0. Upper right panel: Variation of EB with respect to ω. Lower left panel: Variation MRN with respect to binding energy EB. When the binding energy is a minimum… view at source ↗
Figure 7
Figure 7. Figure 7: Left panel: QMRN as a function of f0 for q ≤ 1, we can see that we always have QMRN < 1. Right panel: QMRN as a function of f0 for q ≥ 1. Although for these cases we have q ≥ m = 1, we always nd Q ≤ MRN . According to equation (72), the spinor solutions are given by: ψ± = ei(±φ/2ωt) √4π   f(r)y± (θ) ±if(r) y(θ) ig(r) y±(θ) ±g(r) y(θ)    (127) The radial proles of the spinor amplitudes f(r) a… view at source ↗
Figure 8
Figure 8. Figure 8: Spinor amplitudes f(r) and g(r) with for the dierent values of q and f0. We found that for 0 ≤ q ≤ 1052, and f0  [001, 1], the Dirac eld admits stationary solutions of the form given by equation (127), with the corresponding values of ω shown in [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗
Figure 9
Figure 9. Figure 9: Left panel: Lapse function (r). Right panel: Radial metric function a(r). s = 0, 12, 1 share a wide range of common characteristics when they coupled to gravity. ACKNOWLEDGMENTS The authors wish to thank Carlos Joaquin, Yahir Mio and Elishah Candanosa for many useful discussions and comments. This work was partially supported by DGAPA-UNAM project IN100523. MH-M also acknowledges nancial support from a … view at source ↗
Figure 10
Figure 10. Figure 10: Left panel: Electric potential ϕ(r). Right panel: Electric eld E(r) [PITH_FULL_IMAGE:figures/full_fig_p022_10.png] view at source ↗
Figure 11
Figure 11. Figure 11: Total energy density ρ(r) for four values of q = (060, 090, 100, 101), and four values of the central amplitude of the Dirac eld f0 = (006, 015, 025, 060). For f0 = 006 the density is very small but is dierent from zero. In all cases the energy density is concentrated near the origin [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗

discussion (0)

Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.

Reference graph

Works this paper leans on

18 extracted references · 4 linked inside Pith

  1. [1]

    Pugliese, H

    D. Pugliese, H. Quevedo, J. A. Rueda H, and R. Runi, Physical Review D—Particles, Fields, Gravitation, and Cosmology 88, 024053 (2013)

    2013

  2. [2]

    Jetzer, Physics Reports220, 163 (1992)

    P. Jetzer, Physics Reports220, 163 (1992)

    1992

  3. [3]

    Brito, V

    R. Brito, V. Cardoso, C. A. Herdeiro, and E. Radu, Physics Letters B752, 291 (2016)

    2016

  4. [4]

    I. S. Landea and F. Garc ´ıa, Physical Review D94, 104006 (2016)

    2016

  5. [5]

    Mio and M

    Y. Mio and M. Alcubierre, General Relativity and Gravitation57, 1 (2025)

    2025

  6. [6]

    Finster, J

    F. Finster, J. Smoller, and S.-T. Yau, Physical Review D59, 104020 (1999)

    1999

  7. [7]

    Sun, S.-Y

    S.-X. Sun, S.-Y. Cui, L.-X. Huang, T.-F. Fang, and Y.-Q. Wang, European Physical Journal C84, 699 (2024)

    2024

  8. [8]

    Liang, J.-R

    C. Liang, J.-R. Ren, S.-X. Sun, and Y.-Q. Wang, European Physical Journal C84, 14 (2024)

    2024

  9. [9]

    C. A. R. Herdeiro and E. Radu, Phys. Lett. B777, 8 (2018) , arXiv:1708.05674 [gr-qc]

    Pith/arXiv arXiv 2018

  10. [10]

    Dzhunushaliev and V

    V. Dzhunushaliev and V. Folomeev, Physical Review D99, 084030 (2019) , arXiv:1811.07500

    Pith/arXiv arXiv 2019

  11. [11]

    Finster, J

    F. Finster, J. Smoller, and S.-T. Yau, Communications in Mathematical Physics205, 249 (1999)

    1999

  12. [12]

    C. A. R. Herdeiro, I. Perapechka, E. Radu, and Y. Shnir, Phys. Lett. B824, 136811 (2022) , arXiv:2111.14475 [gr-qc]

    Pith/arXiv arXiv 2022

  13. [13]

    Sun, L.-X

    S.-X. Sun, L.-X. Huang, Z.-H. Zhao, and Y.-Q. Wang, arXiv preprint arXiv:2411.12969 (2024). 22 Figure 10:Left panel: Electric potentialϕ(r). Right panel: ElectriceldE(r). Figure 11:Total energy densityρ(r) for four values ofq= (060,090,100,101), and four values of the central amplitude of the Diraceldf 0 = (006,015,025,060). Forf 0 = 006 the de...

    Pith/arXiv arXiv 2024

  14. [14]

    Alcubierre,Introduction to 3+ 1 numerical relativity, Vol

    M. Alcubierre,Introduction to 3+ 1 numerical relativity, Vol. 140 (OUP Oxford, 2008)

    2008

  15. [15]

    Alcubierre, General Relativity and Gravitation57, 122 (2025)

    M. Alcubierre, General Relativity and Gravitation57, 122 (2025)

    2025

  16. [16]

    J. D. L´ opez and M. Alcubierre, General Relativity and Gravitation55, 72 (2023)

    2023

  17. [17]

    Rezzolla and C

    L. Rezzolla and C. Ecker, SciPost Physics20, 014 (2026)

    2026

  18. [18]

    Finster, J

    F. Finster, J. Smoller, and S.-T. Yau, Physics Letters A259, 431 (1999)

    1999