Massive boson stars: Stability and GW emission in head-on mergers

Bo-Xuan Ge

arxiv: 2512.15242 · v4 · pith:SLUTTGCQnew · submitted 2025-12-17 · 🌀 gr-qc · astro-ph.HE

Massive boson stars: Stability and GW emission in head-on mergers

Bo-Xuan Ge This is my paper

Pith reviewed 2026-05-16 22:08 UTC · model grok-4.3

classification 🌀 gr-qc astro-ph.HE

keywords boson starsself-interacting scalar fieldsgravitational waveshead-on mergersnumerical relativitycompact object stabilityblack hole formationtidal deformability

0 comments

The pith

Quartically self-interacting massive boson stars remain stable only up to the first maximum on their mass curve, and equal-mass head-on collisions produce one of three outcomes depending on compactness.

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

The paper builds equilibrium sequences of massive boson stars with quartic self-interaction and evolves them in head-on collisions. The mass curve develops multiple extrema, but stability holds only until the first maximum, after which models are too compact and collapse from small perturbations. Equal-mass mergers fall into three classes: a merged boson star remnant, black holes forming exactly at contact, or each star collapsing to a black hole before contact. Gravitational-wave energy output follows from the tension between greater compactness boosting emission and lower tidal deformability reducing asymmetry, with non-monotonic patterns appearing at strong self-interaction. The simulations supply a catalogue of initial data and waveforms to enable fast surrogate models across parameter space.

Core claim

Along sequences of equilibrium configurations for quartically self-interacting massive boson stars the mass M plotted against central amplitude |φ_c| develops multiple extrema, yet stability is lost only after the first maximum, with subsequent models highly compact and susceptible to collapse from numerical perturbations, sometimes showing brief double-dive behavior near criticality. Head-on collisions between equal-mass stars fall into three classes: formation of a remnant boson star, black hole creation at the moment of contact, or each star collapsing individually to a black hole before contact occurs. The gravitational-wave energy radiated in these events arises from the tension between

What carries the argument

The mass curve M(|φ_c|) along equilibrium sequences that sets stability thresholds, together with numerical dynamical evolution of head-on collisions that classifies the three outcomes and measures gravitational-wave emission.

If this is right

Configurations past the first mass maximum collapse under numerically induced perturbations, restricting the stable region of parameter space.
Gravitational-wave energies emitted in mergers vary with self-interaction strength through the opposing effects of compactness and tidal deformability.
The pre-contact collapse branch displays non-monotonic energy dependence at large self-interaction strengths.
The generated catalogue of initial conditions and waveforms enables construction of surrogate models for rapid gravitational-wave signal prediction.

Where Pith is reading between the lines

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

Detected gravitational-wave events could be compared against these waveforms to place limits on boson star parameters in nature.
The three-outcome classification may extend to spinning or unequal-mass cases, offering additional channels for observable signals.
Surrogate models trained on this catalogue could speed up parameter estimation for exotic compact object mergers in future detector data.

Load-bearing premise

Stability of boson star configurations changes only at the first maximum of the mass curve M as a function of central amplitude, with all more compact solutions being unstable to collapse.

What would settle it

A long-term stable numerical evolution of a configuration lying past the first maximum on the mass curve would falsify the stability threshold.

Figures

Figures reproduced from arXiv: 2512.15242 by Bo-Xuan Ge.

**Figure 1.** Figure 1: FIG. 1. Illustration of the evolution scheme used in the two-dimensional boson star code. The evolution is performed in [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

**Figure 2.** Figure 2: FIG. 2. Mass curves [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3. The [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4. Evolution of the central amplitude for [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

**Figure 5.** Figure 5: FIG. 5. Time evolution of the central scalar amplitude [PITH_FULL_IMAGE:figures/full_fig_p008_5.png] view at source ↗

**Figure 6.** Figure 6: FIG. 6. Head-on collisions of massive boson stars for four values of the self-coupling parameter [PITH_FULL_IMAGE:figures/full_fig_p010_6.png] view at source ↗

**Figure 7.** Figure 7: FIG. 7. High–compactness regime at [PITH_FULL_IMAGE:figures/full_fig_p012_7.png] view at source ↗

**Figure 8.** Figure 8: FIG. 8. Constraint violation for [PITH_FULL_IMAGE:figures/full_fig_p017_8.png] view at source ↗

**Figure 9.** Figure 9: FIG. 9. Convergence of the gravitational-wave energy for [PITH_FULL_IMAGE:figures/full_fig_p017_9.png] view at source ↗

read the original abstract

We investigate quartically self-interacting massive boson stars by constructing equilibrium sequences and performing dynamical evolutions. The mass curve $M(|\phi_c|)$ along these sequences develops multiple extrema, yet stability changes only at the first maximum; configurations beyond it become highly compact and collapse under numerically induced perturbations, with near-critical models displaying a short-lived double-dive behaviour. Head-on collisions of equal-mass stars yield three distinct outcomes -- boson star remnants, black hole formation at contact, and collapse of each star to a black hole prior to contact. The associated gravitational-wave energies reflect a competition between increasing compactness, which enhances the efficiency of gravitational-wave emission, and decreasing tidal deformability, which suppresses merger asymmetries, and at large self-interaction strengths the collapse-before-contact branch exhibits a pronounced non-monotonic structure. The simulations reported here constitute a substantial catalogue of initial conditions and waveforms, providing a natural basis for constructing surrogate models capable of rapidly predicting gravitational-wave signals across an extended parameter space.

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.

Desk Editor's Note private letter to a colleague

The paper maps three head-on boson-star merger outcomes and a non-monotonic GW energy trend at strong self-interaction, but the stability boundary beyond the first mass maximum rests on untested numerical collapses.

read the letter

The main thing here is that head-on equal-mass collisions of these quartically self-interacting boson stars produce three distinct results: a stable remnant, black-hole formation right at contact, or each star collapsing to a black hole before contact. At large self-interaction the radiated gravitational-wave energy shows a non-monotonic pattern. The authors also supply a catalog of initial data and waveforms that could feed surrogate models later. That is the concrete addition to the literature. The work extends earlier boson-star merger runs by scanning higher self-interaction strengths and catching the pre-contact collapse branch. The interpretation that compactness boosts GW efficiency while lower tidal deformability reduces asymmetries is straightforward and lines up with the trends they report. The numerics themselves look like standard numerical-relativity practice. The soft spot is the stability analysis. The mass curve has multiple extrema, yet the authors tie the stability change to the first maximum only because higher-compactness models collapse when evolved. No convergence tests, isolated-star control runs, or resolution studies are described, so it is not clear whether that collapse is physical or a grid artifact. If the latter, the clean separation into three outcomes and the non-monotonic GW structure both become less reliable. This is the sort of paper that belongs in a numerical-relativity or modified-gravity journal. Readers working on boson stars or alternative compact objects will get usable initial data and waveforms from it. It is not a big conceptual leap, but the results are specific enough that a referee should check the numerical robustness. I would send it to peer review.

Referee Report

2 major / 2 minor

Summary. The paper investigates quartically self-interacting massive boson stars by constructing equilibrium sequences and performing dynamical evolutions. The mass curve M(|φ_c|) develops multiple extrema, yet stability changes only at the first maximum; configurations beyond it become highly compact and collapse under numerically induced perturbations, with near-critical models showing short-lived double-dive behaviour. Head-on collisions of equal-mass stars produce three distinct outcomes—boson star remnants, black hole formation at contact, and collapse of each star to a black hole prior to contact—with associated gravitational-wave energies reflecting a competition between increasing compactness (enhancing emission efficiency) and decreasing tidal deformability (suppressing merger asymmetries). At large self-interaction strengths the collapse-before-contact branch exhibits non-monotonic structure. The simulations provide a catalogue of initial conditions and waveforms for surrogate models.

Significance. If the stability classification and merger outcomes are robust, the work provides a substantial numerical catalogue of boson-star head-on mergers and waveforms that can serve as a basis for surrogate models across an extended parameter space. This is a concrete strength for the field, as it supplies reproducible initial data and signals for exotic compact-object gravitational-wave phenomenology without relying on fitted parameters in the central derivations.

major comments (2)

[Stability analysis and equilibrium sequences] The assertion that stability changes only at the first maximum of M(|φ_c|), with post-maximum configurations collapsing under numerically induced perturbations, underpins the separation into three distinct merger outcomes (including the pre-contact collapse branch). This claim rests on dynamical evolutions showing collapse, yet no resolution studies, isolated-star control runs, or convergence tests are described to establish that the instability is physical rather than a grid artifact. If the collapse is resolution-dependent, the reported GW energy trends and the compactness-versus-tidal-deformability competition become unreliable.
[Dynamical evolutions and gravitational-wave emission] The non-monotonic structure reported for the collapse-before-contact branch at large self-interaction strengths, and the overall GW energy scaling, depend directly on the robustness of the pre-contact versus at-contact black-hole formation distinction. Without quantified convergence for the near-critical double-dive behaviour and the stability threshold, these trends cannot be taken as physical predictions.

minor comments (2)

[Abstract] The abstract introduces 'near-critical models displaying a short-lived double-dive behaviour' without a concise definition or reference to the relevant figure or section; this should be clarified early in the main text.
[Notation and equilibrium construction] Notation for the central scalar amplitude |φ_c| and its relation to the equilibrium sequence parameter should be stated explicitly once in the methods or equilibrium-construction section to avoid ambiguity in later figures.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the careful and constructive review of our manuscript on massive boson stars. The concerns regarding the lack of documented resolution studies for the stability classification and the convergence of the dynamical merger simulations are well taken. We address each major comment below and will revise the manuscript to incorporate additional tests and clarifications.

read point-by-point responses

Referee: [Stability analysis and equilibrium sequences] The assertion that stability changes only at the first maximum of M(|φ_c|), with post-maximum configurations collapsing under numerically induced perturbations, underpins the separation into three distinct merger outcomes (including the pre-contact collapse branch). This claim rests on dynamical evolutions showing collapse, yet no resolution studies, isolated-star control runs, or convergence tests are described to establish that the instability is physical rather than a grid artifact. If the collapse is resolution-dependent, the reported GW energy trends and the compactness-versus-tidal-deformability competition become unreliable.

Authors: We agree that the manuscript does not currently describe resolution studies or convergence tests for the isolated-star evolutions. The collapses observed beyond the first mass maximum are consistent with the turning-point method and occur under small numerical perturbations, but we acknowledge that this requires quantitative support to rule out grid artifacts. In the revised manuscript we will add a dedicated subsection presenting resolution studies for representative models, including near-critical cases exhibiting the short-lived double-dive behaviour. These tests will confirm convergence of the collapse dynamics and thereby strengthen the separation into the three merger outcomes. revision: yes
Referee: [Dynamical evolutions and gravitational-wave emission] The non-monotonic structure reported for the collapse-before-contact branch at large self-interaction strengths, and the overall GW energy scaling, depend directly on the robustness of the pre-contact versus at-contact black-hole formation distinction. Without quantified convergence for the near-critical double-dive behaviour and the stability threshold, these trends cannot be taken as physical predictions.

Authors: We concur that the reported non-monotonic GW energy trends and the distinction between pre-contact and at-contact collapse require demonstrated numerical convergence. The current text does not include quantified convergence tests for the head-on collision runs. We will add convergence analyses for selected models in the collapse-before-contact branch at large self-interaction strengths, showing that the non-monotonic structure persists under increased resolution. This will substantiate the physical competition between compactness and tidal deformability as the driver of the observed GW energy patterns. revision: yes

Circularity Check

0 steps flagged

No circularity: outcomes are direct numerical outputs

full rationale

The paper constructs equilibrium sequences by solving the Einstein-Klein-Gordon system for quartically self-interacting massive boson stars and evolves the initial-value problem in full numerical relativity. Reported outcomes (remnant types, GW energies, pre-contact collapse) follow directly from the time-dependent simulations without any parameter fitted to a subset of data and then renamed as a prediction. The statement that stability changes only at the first mass maximum is an inference from observed collapse under grid perturbations, not a self-definitional or self-citation reduction. No ansatz is smuggled via prior work, no uniqueness theorem is imported from the same authors, and no known empirical pattern is merely relabeled. The derivation chain is therefore self-contained against the numerical evolution code.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 0 invented entities

The central claims rest on the standard Einstein-Klein-Gordon system with a quartic potential. No new entities are postulated; the scalar field and its self-interaction are taken from prior boson-star literature. The only free parameter is the self-interaction coupling strength, which is scanned rather than fitted to data.

free parameters (1)

quartic self-interaction strength
Scanned across a range of values to produce the reported sequences and merger branches; its specific numerical values are not given in the abstract.

axioms (2)

standard math The Einstein-Klein-Gordon equations govern the dynamics of the massive scalar field with quartic self-interaction.
Invoked throughout the equilibrium construction and dynamical evolutions; standard in the boson-star literature.
domain assumption Numerical perturbations are sufficient to trigger collapse of configurations past the first mass maximum.
Used to classify stability; appears in the description of post-maximum behaviour.

pith-pipeline@v0.9.0 · 5462 in / 1653 out tokens · 39849 ms · 2026-05-16T22:08:31.164473+00:00 · methodology

discussion (0)

Lean theorems connected to this paper

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

IndisputableMonolith/Cost/FunctionalEquation.lean washburn_uniqueness_aczel unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

stability changes only at the first maximum; configurations beyond it become highly compact and collapse under numerically induced perturbations
IndisputableMonolith/Foundation/AlexanderDuality.lean alexander_duality_circle_linking unclear

?

unclear
Relation between the paper passage and the cited Recognition theorem.

mass curve M(|φ_c|) develops multiple extrema, yet stability changes only at the first maximum

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.

Forward citations

Cited by 2 Pith papers

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score.

Timing-Window Mechanism for Chain-Like Transients in Collisions of Radially Excited Boson Stars
gr-qc 2026-05 unverdicted novelty 6.0

Chain-like transients in boson star collisions are governed by a timing window set by the binary collision time relative to isolated breathing clocks rather than excitation level alone.
Boson star-black hole binaries: initial data and head-on collisions
gr-qc 2026-04 unverdicted novelty 6.0

A one-body conformal-factor correction stabilizes boson star-black hole initial data, enabling gravitational-wave analysis that shows higher multipoles can discriminate mixed mergers from pure black-hole binaries.