Two-fermion spherically symmetric Einstein-Dirac solitons: general spinor ansatz, limiting configuration of the mass-energy spirals, strong-coupling stability, and rigorous global existence

The static, spherically symmetric Einstein-Dirac system for two fermions in a singlet spinor state, introduced by Finster, Smoller, and Yau in 1999, admits an infinite family of particlelike (soliton-like) solutions: normalizable spinor fields in an asymptotically flat spacetime, organized by a rotation number n, whose mass-energy curves form spirals converging to a limiting configuration. The physically motivated coupling to a Maxwell field was taken up by the same authors shortly afterwards [2], and the local, weak-coupling part of the existence question was settled rigorously by Rota Nodari (2009) and, in the almost-Newtonian regime, by Stuart (2010). The surviving open core has four linked components: (i) whether the general spherically symmetric spinor ansatz (both spinors nonzero) admits additional normalizable, asymptotically flat solutions beyond the restricted ansatz actually integrated; (ii) whether the limiting configuration of each spiral is a single point or a nontrivial limit cycle, a dichotomy the source's Conley-index argument leaves undetermined; (iii) a complete stability classification in the strong-coupling spiral regime beyond the qualitative branch results; and (iv) a rigorous analytic proof of the global bifurcation structure, including the infinite accumulation of solutions at the limiting configuration. All four components remain open, with substantial numerical evidence but no proof.

Background

The Einstein-Dirac equations couple Einstein's field equations of general relativity to the Dirac equation, the relativistic wave equation for spin-1/2 fermions, so that the spacetime geometry is sourced by the energy-momentum of the fermionic wave functions. Because the spin of a single fermion singles out a spatial direction, a static, spherically symmetric configuration requires (at least) two fermions of opposite spin combined into a singlet spinor state. Finster, Smoller, and Yau (Phys. Rev. D 59, 104020 (1999)) derived the resulting coupled equations and constructed, by numerical shooting - integrating the equations outward from the origin and tuning the free initial data until the boundary conditions at infinity are met - particlelike solutions: bound states in which the spinor wave functions are normalizable (square-integrable, so the total fermion number is finite) and the spacetime is asymptotically flat (the metric approaches the Minkowski form at large radius, with a well-defined ADM mass, i.e. the total gravitational mass measured at infinity).

After a reduction to two real radial spinor functions $\alpha(r)$, $\beta(r)$ and two metric functions, the solutions are organized by a rotation number $n = 0, 1, 2, \ldots$, the winding number of the plane curve $(\alpha(r), \beta(r))$; $n=0$ is the ground state and $n \geq 1$ are excited states. The equations admit a scaling freedom, which the numerics exploit by temporarily fixing the metric normalization and the fermion mass and rescaling afterwards; the physical solutions are then characterized by the rest mass $m$, the fermion frequency $\omega < |m|$ (so the binding energy $\omega - |m|$ is negative and the state is bound), and the ADM mass. For weak coupling (small $m$) the authors found one linearly stable solution per $n$ (stable under spherically symmetric perturbations). For stronger coupling, plotting the binding energy versus $m$ traces, for each $n$, a spiral curve that tightens toward a limiting configuration $\Gamma_n$; a topological argument based on the Conley index - an algebraic-topological invariant of isolated invariant sets of a flow that is preserved under continuous deformation of parameters - shows that $\Gamma_n$ must be either a single point or a topological circle $S^1$ (a limit cycle), but does not determine which. The same machinery gives only qualitative stability information along the spiral: solutions on the lower branch are stable, those on the upper branch unstable, with changes of the Conley index at the bifurcation points where pairs of solutions appear, while direct numerical linearization of the perturbation equations becomes unreliable precisely in the spiral regime.

The numerical constructions rest on a restricted spinor ansatz. Appendix A of the paper derives the most general static, spherically symmetric two-fermion ansatz, in which the wave function is governed by the Dirac equations (A.3) and (A.4) for two complex two-spinors $\Phi$ and $\Xi$; the restricted ansatz actually integrated corresponds to the two special cases $\Xi = 0$ and $\Phi = 0$. Whether the general, mixed ansatz admits solutions with both $\Phi$ and $\Xi$ simultaneously normalizable and decaying at infinity - genuinely additional particlelike solutions - was left open. The authors also noted that their existence and multiplicity statements (for instance, that parameter points on $\Gamma_n$ carry infinitely many solutions) rest on numerical integration and topological continuation rather than analytic proof, and they deferred the physically motivated coupling to an electromagnetic (Maxwell) field. The two-fermion system has since become the prototype for an active literature on Dirac stars - solitonic solutions of the Einstein-Dirac equations in which many fermions occupy the same quantum state, the fermionic analogues of boson stars (self-gravitating bound states of a complex scalar field) - whose solution branches inherit the same spiral bifurcation structure.

Problem Statement

Complete the mathematical characterization of the regular, normalizable, asymptotically flat solutions of the static, spherically symmetric two-fermion Einstein-Dirac system of Finster, Smoller, and Yau (Phys. Rev. D 59, 104020 (1999)) - which solutions exist, how they organize into bifurcation curves, and which are stable. Four linked components of the question remain open:

(i) General spinor ansatz. Determine whether the general spherically symmetric spinor ansatz of Appendix A - in which both complex two-spinors $\Phi$ and $\Xi$ are nonzero, i.e. the solution is a superposition of the two restricted ansätze used in the numerical construction - admits solutions of the Einstein-Dirac equations (A.3), (A.4) that are regular at the origin, normalizable, and decaying so that the spacetime is asymptotically flat, and that are not reducible (via the rescaling and phase freedom of the system) to the restricted-ansatz solutions; or prove that no such additional solutions exist.

(ii) Nature of the limiting configuration. For each rotation number $n$, the curve of solutions in the (rest mass, binding energy) plane forms a spiral which, as the shooting parameter $\alpha_1$ (the leading coefficient of the spinor component $\alpha$ at the origin) tends to infinity, converges to a limiting configuration $\Gamma_n$. A Conley-index argument shows that $\Gamma_n$ is either a single point or a nontrivial limit cycle $S^1$, but does not decide which. Determine, for each $n$, which alternative holds.

(iii) Strong-coupling stability. Obtain a complete stability classification of the solutions in the strong-coupling regime where the mass-energy curves form tight spirals: a spectral characterization of the unstable modes of the linearized spherically symmetric perturbation equations along the spiral branches, going beyond the qualitative Conley-index conclusions (lower branch stable, upper branch unstable, index changes at the bifurcation points) and covering the deep spiral regime where the source paper's numerical linearization 'no longer allows simple conclusions.'

(iv) Rigorous global existence. Prove analytically the structural results established numerically and topologically: the infinite accumulation of solutions at the limiting configuration $\Gamma_n$ (infinitely many solutions at parameter points on $\Gamma_n$), and the existence and global shape of the bifurcation curves (spiraling global solution branches), going beyond the perturbative existence of the weak-coupling branch.

These components serve one objective - a complete existence-and-stability picture of the two-fermion Einstein-Dirac soliton spectrum - and are mutually dependent: a stability classification near $\Gamma_n$ presupposes the nature of $\Gamma_n$ settled in (ii), and the accumulation theorem in (iv) requires control of the same global bifurcation curve.

References

1. F. Finster, J. Smoller, S.-T. Yau, Particlelike solutions of the Einstein-Dirac equations, Phys. Rev. D 59, 104020 (1999), DOI 10.1103/PhysRevD.59.104020, arXiv:gr-qc/9801079, https://arxiv.org/abs/gr-qc/9801079
2. F. Finster, J. Smoller, S.-T. Yau, Particle-Like Solutions of the Einstein-Dirac-Maxwell Equations, Phys. Lett. A 259, 431-436 (1999), DOI 10.1016/S0375-9601(99)00457-0, arXiv:gr-qc/9802012, https://arxiv.org/abs/gr-qc/9802012
3. S. Rota Nodari, Perturbation Method for Particle-like Solutions of the Einstein-Dirac Equations, Ann. Henri Poincaré 10, 1377-1393 (2009), DOI 10.1007/s00023-009-0015-x, arXiv:0906.5043, https://arxiv.org/abs/0906.5043
4. S. Rota Nodari, Perturbation method for particle-like solutions of the Einstein-Dirac-Maxwell equations, C. R. Math. Acad. Sci. Paris 348, 791-794 (2010), DOI 10.1016/j.crma.2010.06.003, arXiv:0912.4246, https://arxiv.org/abs/0912.4246
5. D. Stuart, Existence and Newtonian limit of nonlinear bound states in the Einstein-Dirac system, J. Math. Phys. 51, 032501 (2010), DOI 10.1063/1.3294085, arXiv:0906.4661, https://arxiv.org/abs/0906.4661
6. E. Daka, N. N. Phan, B. Kain, Perturbing the ground state of Dirac stars, Phys. Rev. D 100, 084042 (2019), DOI 10.1103/PhysRevD.100.084042, arXiv:1910.09415, https://arxiv.org/abs/1910.09415
7. D. Bakucz Canário, S. Lloyd, K. Horne, C. A. Hooley, Infinite-redshift localized states of Dirac fermions under Einsteinian gravity, Phys. Rev. D 102, 084049 (2020), DOI 10.1103/PhysRevD.102.084049, arXiv:2008.10455, https://arxiv.org/abs/2008.10455
8. P. E. D. Leith, C. A. Hooley, K. Horne, D. G. Dritschel, Nonlinear effects in the excited states of many-fermion Einstein-Dirac solitons, Phys. Rev. D 104, 046024 (2021), DOI 10.1103/PhysRevD.104.046024, arXiv:2105.12672, https://arxiv.org/abs/2105.12672
9. P. E. D. Leith, C. A. Hooley, K. Horne, D. G. Dritschel, Fermion self-trapping in the optical geometry of Einstein-Dirac solitons, Phys. Rev. D 101, 106012 (2020), DOI 10.1103/PhysRevD.101.106012, arXiv:2002.02747, https://arxiv.org/abs/2002.02747
10. P. E. D. Leith, C. A. Hooley, K. Horne, D. G. Dritschel, 'Stealth' singularities from self-gravitating fermions, Phys. Rev. D 113, L121507 (2026), DOI 10.1103/81mh-mdt5, arXiv:2508.10604, https://arxiv.org/abs/2508.10604
11. J. L. Blázquez-Salcedo, C. Knoll, Constructing spherically symmetric Einstein-Dirac systems with multiple spinors: Ansatz, wormholes and other analytical solutions, Eur. Phys. J. C 80, 174 (2020), DOI 10.1140/epjc/s10052-020-7706-3, https://doi.org/10.1140/epjc/s10052-020-7706-3
12. B. Kain, Einstein-Dirac system in semiclassical gravity, Phys. Rev. D 107, 124001 (2023), DOI 10.1103/PhysRevD.107.124001, arXiv:2304.10627, https://arxiv.org/abs/2304.10627
13. M. Hernández Márquez, M. Alcubierre, Charged Dirac stars, arXiv:2607.26374 (2026), https://arxiv.org/abs/2607.26374
14. F. Finster, J. Smoller, S.-T. Yau, Non-Existence of Black Hole Solutions for a Spherically Symmetric, Static Einstein-Dirac-Maxwell System, Commun. Math. Phys. 205, 249-262 (1999), DOI 10.1007/s002200050675, https://doi.org/10.1007/s002200050675
15. C. Liang, S.-X. Sun, J.-R. Ren, Y.-Q. Wang, Multi-state Dirac stars, arXiv:2306.11437 (2023), https://arxiv.org/abs/2306.11437
16. S. Vignolo, L. Fabbri, A note on spinor fields in spherical symmetry, Int. J. Geom. Methods Mod. Phys. 23, 2650054 (2026), DOI 10.1142/S0219887826500544, arXiv:2604.13582, https://arxiv.org/abs/2604.13582