pith. sign in

arxiv: 1907.08193 · v1 · pith:J3TE2EP3new · submitted 2019-07-18 · 🌀 gr-qc · astro-ph.HE

The three dynamical fates of Boson Stars

F. S. Guzman This is my paper

Pith reviewed 2026-05-24 19:46 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HE
keywords boson starsEinstein-Klein-Gordon systemnumerical evolutionstabilityblack hole formationscalar fieldspherical symmetry
0
0 comments X

The pith

Spherically symmetric boson stars evolve into one of three late-time states: stable, collapse to black hole, or explosion.

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

The paper shows through numerical solution of the spherically symmetric Einstein-Klein-Gordon system that boson star initial data lead to three distinct outcomes depending on the parameters chosen. Stable cases remain unchanged at late times, while some unstable cases form black holes and others disperse outward. These behaviors are tracked using a constrained evolution approach that solves the geometry constraints at each step while advancing the scalar field. A sympathetic reader would care because the classification determines whether such objects could persist as long-lived astrophysical candidates or must quickly disappear.

Core claim

The central claim is that spherically symmetric boson stars exhibit three types of late-time behavior: stable configurations, unstable bounded configurations that collapse to form black holes, and unstable unbounded configurations that explode. These results follow from solving the spherically symmetric Einstein-Klein-Gordon system with initial conditions that correspond to equilibrium boson star solutions, using constrained evolution and the method of lines.

What carries the argument

Constrained evolution of the spherically symmetric Einstein-Klein-Gordon system for a complex scalar field, which solves the constraint equations for the geometry on the fly.

If this is right

  • Boson stars with certain initial parameters remain stable and persist indefinitely.
  • Unstable bounded boson stars collapse to form black holes.
  • Unstable unbounded boson stars disperse and explode outward.
  • The specific fate is set by the parameters of the equilibrium boson star solution used as initial data.

Where Pith is reading between the lines

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

  • The three fates narrow the range of boson star parameters that could correspond to long-lived objects in the universe.
  • The explosion case might produce observable scalar radiation or gravitational wave signals in more realistic settings.
  • Analogous three-outcome dynamics could appear in other self-gravitating complex scalar field models beyond the spherical case studied here.

Load-bearing premise

The chosen initial data precisely match equilibrium boson star solutions and the numerical scheme faithfully captures the qualitative dynamics without resolution-dependent artifacts.

What would settle it

A simulation starting from the same boson star initial data but with much higher resolution that produces a qualitatively different late-time outcome, such as persistent oscillation instead of collapse or explosion.

Figures

Figures reproduced from arXiv: 1907.08193 by F. S. Guzman.

Figure 1
Figure 1. Figure 1: FIG. 1: Sequences of equilibrium configurations for two val [PITH_FULL_IMAGE:figures/full_fig_p002_1.png] view at source ↗
Figure 2
Figure 2. Figure 2: FIG. 2: (top) Convergence of [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. Figure 3: FIG. 3: (Top) Maximum of the metric function [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figure 4
Figure 4. Figure 4: FIG. 4: (Top) Snapshots of the metric functions for configura [PITH_FULL_IMAGE:figures/full_fig_p006_4.png] view at source ↗
read the original abstract

In this manuscript the three types of late-time behavior of spherically symmetric Boson Stars are presented, these are: stable configurations, unstable bounded that collapse to form black holes and unstable unbounded that explode. The results are found by solving the spherically symmetric Einstein-Klein-Gordon system of equations for a complex scalar field with initial conditions corresponding to a Boson Star. The solution is based on a constrained evolution that uses the method of lines for the scalar field and solves the constraint equations for the geometry on the fly.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Referee Report

2 major / 1 minor

Summary. The manuscript numerically evolves the spherically symmetric Einstein-Klein-Gordon system for complex scalar fields with initial data taken from equilibrium boson-star solutions. It reports three distinct late-time fates: stable equilibria that persist, unstable bounded configurations that collapse to black holes, and unstable unbounded configurations that disperse or explode. The evolution employs a constrained scheme in which the scalar field is advanced via the method of lines while the Hamiltonian and momentum constraints are solved at each step for the metric variables.

Significance. If the numerical results prove robust under refinement and constraint monitoring, the classification supplies a concrete dynamical taxonomy for boson stars that complements existing linear stability analyses and may inform studies of scalar-field compact objects or dark-matter candidates.

major comments (2)
  1. [Numerical Methods] Numerical Methods section: the constrained evolution is described without any reported convergence tests (spatial or temporal resolution series), constraint-violation residual plots, or early-time comparisons of the evolved fields against the known stationary boson-star solutions. These diagnostics are required to establish that the reported collapse and explosion outcomes are not artifacts of discretization error, artificial dissipation, or secular constraint drift.
  2. [Results] Results section: the three fates are asserted on the basis of the chosen initial data, yet no quantitative measure is supplied of how closely the initial data satisfy the equilibrium ODEs to machine precision, nor are error bars or resolution dependence shown for the critical thresholds separating the three regimes.
minor comments (1)
  1. The abstract and introduction could clarify the precise meaning of 'unbounded configurations that explode' (e.g., whether this refers to dispersal with the scalar field amplitude decaying at large radii or to a different diagnostic).

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for their thorough review and valuable suggestions. We will revise the manuscript to include the requested numerical diagnostics and quantitative measures as detailed in our point-by-point responses below.

read point-by-point responses

  1. Referee: [Numerical Methods] Numerical Methods section: the constrained evolution is described without any reported convergence tests (spatial or temporal resolution series), constraint-violation residual plots, or early-time comparisons of the evolved fields against the known stationary boson-star solutions. These diagnostics are required to establish that the reported collapse and explosion outcomes are not artifacts of discretization error, artificial dissipation, or secular constraint drift.

    Authors: We agree with the referee that these numerical validation tests are crucial. Although the manuscript focuses on the physical classification, we acknowledge the omission of these diagnostics in the current version. In the revised manuscript, we will add convergence tests with multiple spatial resolutions, plots of constraint violations, and comparisons of the evolved fields at early times to the initial data. This will confirm that the reported fates are not numerical artifacts. revision: yes

  2. Referee: [Results] Results section: the three fates are asserted on the basis of the chosen initial data, yet no quantitative measure is supplied of how closely the initial data satisfy the equilibrium ODEs to machine precision, nor are error bars or resolution dependence shown for the critical thresholds separating the three regimes.

    Authors: We will include in the revision a quantitative assessment of the initial data accuracy by reporting the residual of the equilibrium equations. For the critical thresholds, we will present results from simulations at different resolutions to demonstrate convergence and provide estimates of the uncertainty in the separating values. revision: yes

Circularity Check

0 steps flagged

No circularity: direct numerical evolution of EKG system yields emergent fates

full rationale

The manuscript reports three late-time behaviors observed from constrained numerical integration of the spherically symmetric Einstein-Klein-Gordon system, initialized with equilibrium boson-star solutions. These outcomes (stable, collapse to BH, dispersal) are dynamical results of the evolution equations rather than quantities fitted to data or defined in terms of themselves. No load-bearing self-citation, ansatz smuggling, or renaming of known results appears in the provided text; the central claim rests on the numerical scheme applied to external analytic equilibria, which is independent of the reported fates. This is the expected non-finding for a pure simulation study.

Axiom & Free-Parameter Ledger

0 free parameters · 2 axioms · 0 invented entities

The paper rests on the standard Einstein-Klein-Gordon system and the assumption that supplied initial data are exact boson-star equilibria; no free parameters, new entities, or ad-hoc axioms are introduced in the abstract.

axioms (2)
  • standard math The Einstein-Klein-Gordon system governs the dynamics of a complex scalar field minimally coupled to gravity.
    Invoked as the governing equations for the evolution.
  • domain assumption Initial data correspond to equilibrium boson-star solutions.
    Stated directly in the abstract as the starting point for the simulations.

pith-pipeline@v0.9.0 · 5599 in / 1181 out tokens · 23425 ms · 2026-05-24T19:46:14.418443+00:00 · methodology

discussion (0)

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

Lean theorems connected to this paper

Citations machine-checked in the Pith Canon. Every link opens the source theorem in the public Lean library.

What do these tags mean?
matches
The paper's claim is directly supported by a theorem in the formal canon.
supports
The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
extends
The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
uses
The paper appears to rely on the theorem as machinery.
contradicts
The paper's claim conflicts with a theorem or certificate in the canon.
unclear
Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.

Reference graph

Works this paper leans on

22 extracted references · 22 canonical work pages · 6 internal anchors

  1. [1]

    F. S. Guzm´ an and L. A. Ure˜ na-L´ opez, ApJ 645, 814 (2006). ArXiV: astro-ph/0603613

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  2. [2]

    Y-F. Tuan, R. Narayan and M. J. Rees, ApJ 606, 1112 (2004)

    work page 2004

  3. [3]

    D. F. Torres, Nucl. Phys. B 26 , 377 (2002)

    work page 2002

  4. [4]

    F. S. Guzm´ an, Phys. Rev. D 73 , 021501 (2006). ArXiV: gr-qc/0512081

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  5. [5]

    Scalar Field Dark Matter: non-spherical collapse and late time behavior

    A. Bernal and F. S. Guzm´ an, Phys. Rev. D 74 , 063504 (2006). ArXiV: astro-ph/0608523

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  6. [6]

    Scalar Field Dark Matter: head-on interaction between two structures

    A. Bernal and F. S. Guzm´ an, Phys. Rev. D 74 , 103002 (2006). ArXiV: astro-ph/0610682

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  7. [7]

    E. W. Mielke and J. A. V´ elez-P´ erez, Phys. Rev. D 75 043504 (2007). Ibid, Phys.Lett. B 671 174 (2009)

    work page 2007

  8. [8]

    Jetzer, Phys

    Ph. Jetzer, Phys. Rep. 220 163, 1992

    work page 1992

  9. [9]

    F. E. Schunck, E. W. Mielke, Class. Quantum Grav. 20 R301-R356

    work page

  10. [10]

    Evolution of 3D Boson Stars with Waveform Extraction

    J. Balakrishna, R. Bondarescu, G. Daues, F. S. Guzm´ an, E. Seidel, Class. Quantum Grav. 23, 2631-2652 (2006). ArXiV: gr-qc/0602078

    work page internal anchor Pith review Pith/arXiv arXiv 2006

  11. [11]

    Head-on collisions of boson stars

    C. Palenzuela, I. Olabarrieta, L. Lehner, S. Liebling, Phys. Rev. D 75 , 064005, 2007. ArXiV: gr-qc/0612067

    work page internal anchor Pith review Pith/arXiv arXiv 2007

  12. [12]

    Seidel and W-M

    E. Seidel and W-M. Suen, Phys. Rev. D 42 , 384 (1990)

    work page 1990

  13. [13]

    Balakrishna, E

    J. Balakrishna, E. Seidel and W-M. Suen, Phys. Rev. D 58, 104004 (1998)

    work page 1998

  14. [14]

    F. S. Guzm´ an, Phys. Rev. D 70 , 044033 (2004)

    work page 2004

  15. [15]

    Gleiser, Phys

    M. Gleiser, Phys. Rev. D 38 , 2376 (1988)

    work page 1988

  16. [16]

    F. V. Kusmartsev, et al., Phys. Rev. D 43, 3895 - 3901 (1991)

    work page 1991

  17. [17]

    F. E. Schunck and D. F. Torres, Int. J. Mod. Phys. D 9 (2000) 601

    work page 2000

  18. [18]

    Ruffini and S

    R. Ruffini and S. Bonazolla, Phys. Rev. 187, 1767 (1969)

    work page 1969

  19. [19]

    S. H. Hawley and M. W. Choptuik, Phys. Rev. D 62 , 104024 (2000)

    work page 2000

  20. [20]

    Alcubierre, Introduction to 3+1 Numerical Relativity , Oxford University Press, 2008

    M. Alcubierre, Introduction to 3+1 Numerical Relativity , Oxford University Press, 2008

    work page 2008

  21. [21]

    F. S. Guzm´ an, Phys. Rev. D 73 021501(R) (2006)

    work page 2006

  22. [22]

    F. S. Guzm´ an, J. M. Rueda Becerril, Phys. Rev. D 80 , 084023 (2009)

    work page 2009