Pith. sign in

REVIEW 2 major objections 4 minor 43 references

Head-on Collisions of Boson Stars with Bowen-York Type Initial Data

T0 review · 2 major / 4 minor · reviewed 2026-07-13 · grok-4.5

Pith's one-line read A Bowen-York-style initial-data recipe for boson stars recovers known head-on results and shows boson-star binaries radiate more gravitational-wave energy than black-hole binaries, while mixed boson-star–black-hole binaries radiate less.

desk verdict Solid, usable Bowen-York-style initial data for boson-star binaries that cleanly recovers known radiation rankings; incremental but ready for peer review. read the letter →

arxiv 2607.09494 v1 pith:5QBKP4A3 submitted 2026-07-10 gr-qc

classification gr-qc
keywords bosonstarsBowen-Yorkinitialdatanumericalrelativityhead-oncollisionsgravitationalwavesmixedbinariesconstraintequations
verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

This paper introduces a simple way to build initial data for numerical-relativity simulations of boson stars that is deliberately patterned after the classic Bowen-York puncture construction used for black holes. The method solves the momentum constraint with analytic, superposable expressions that encode linear momentum, then solves the Hamiltonian constraint once, so that early post-Newtonian orbital parameters can be imported without solving a full elliptic system for every configuration. After testing a single boosted boson star, the authors evolve equal-mass head-on collisions of two boson stars and of a boson star with a black hole over a range of momenta. The radiated energies and ring-down frequencies match earlier literature: pure boson-star collisions emit a larger fraction of their mass-energy in gravitational waves than the corresponding black-hole binaries, while mixed boson-star–black-hole collisions emit less. The practical payoff is a lightweight initial-data pipeline that can later be extended to spinning, inspiraling mixed binaries while remaining continuous with the post-Newtonian regime.

What carries the argument

Bowen-type analytic solutions of the conformal-transverse-traceless momentum constraint for an extended scalar-field source (Eq. 13 with the boson-star momentum density of Eq. 31), which can be superposed linearly and then fed into a single Hamiltonian-constraint solve for the conformal factor.

What would settle it

A high-resolution head-on or quasi-circular boson-star binary evolved from these initial data whose measured gravitational-wave energy differs systematically from an independent, fully constrained initial-data construction for the same masses and momenta.

Watch

Extended reading notes

Core claim

The authors demonstrate that a Bowen-York-type construction—analytic extrinsic curvature for each compact object superposed with conformally rescaled scalar-field sources, followed by a single Hamiltonian solve—produces constraint-satisfying initial data whose subsequent evolution reproduces the known gravitational-wave hierarchies for head-on boson-star and mixed boson-star–black-hole collisions: boson-star binaries radiate more energy than black-hole binaries of the same mass and momentum, while mixed systems radiate less.

Load-bearing premise

That starting from conformally flat, spherically symmetric scalar profiles plus the analytic Bowen extrinsic curvature still yields physically reliable radiated energies after the inevitable initial oscillations and constraint relaxation.

Share X Bluesky LinkedIn Reddit HN

Editorial analysis

A structured set of objections, weighed in public.

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

Referee Report

2 major / 4 minor

Summary. The manuscript constructs Bowen-York-type initial data for boson stars by deriving an analytic solution of the conformal momentum constraint for a mini-boson-star source (Eqs. 13, 29–32), then solving the Hamiltonian constraint after a compact-support superposition (Eqs. 36–39). Single-star tests recover the expected O(P^{2}) ADM-mass growth and the known normal-mode frequencies of the stationary mini-boson star. Equal-mass head-on BBS and BHBS collisions are evolved for three central amplitudes and three momenta; the extracted (2,0) waveforms, radiated energies, apparent-horizon masses and QNM parameters are reported in Tables II–III. The central claim is that the method is effective because the resulting radiation hierarchy (BBS > BBH > BHBS) and remnant properties reproduce earlier literature results obtained with different initial-data constructions.

Significance. If the construction extends without major modification to spinning and inspiraling configurations, it supplies a simple, post-Newtonian-compatible route to mixed compact-object binaries that contain boson stars—precisely the class of systems for which constraint-satisfying data remain comparatively scarce. The explicit recovery of the Hawley–Choptuik frequency and the quantitative match of the energy ranking to independent groups constitute genuine validation strengths. The work is therefore a useful methodological contribution even though the physics results themselves are confirmatory rather than novel.

major comments (2)
  1. [Sections V–VI, Tables II–III] Sections V–VI and Tables II–III report radiated energies (and their percentages of ADM mass) and QNM parameters without any resolution study or error bar. The finest grid spacing also changes between pure BBS runs (Δx = 0.125 µ⁻¹) and runs that contain black holes (Δx = 0.03125 µ⁻¹). Because the claimed hierarchy rests on differences of only a few parts in 10⁴ of the ADM energy, at least a two-resolution comparison for one representative BBS and one BHBS case is required to demonstrate that truncation error does not reverse the ordering or shift the remnant masses at the quoted precision.
  2. [Section IV] Section IV shows that conformal flatness plus spherical symmetry excites persistent normal-mode oscillations whose energy content is never quantified relative to the gravitational-wave energy later extracted from the binaries. A short estimate (or a controlled comparison with a non-conformally-flat single-star boost) is needed to confirm that these initial-data artifacts remain sub-dominant for the radiation budgets listed in Tables II–III.
minor comments (4)
  1. [Section III] Notation for the conformally rescaled momentum density oscillates between eSi, ˜Si and Si without a single consistent definition; a short glossary or a uniform choice would improve readability.
  2. [Figures 3 and 6] Figure 3 (bottom panel) and Figure 6 (top panel) would benefit from an explicit statement of the time unit and from a vertical scale that makes the two oscillation frequencies easier to read by eye.
  3. [Abstract and Section VI] The phrase “equivalent black hole binaries” is used repeatedly; a one-sentence clarification that the comparison is performed at equal ADM mass and equal initial linear momentum would remove any ambiguity.
  4. [Section V] Reference [17] is cited for the two-dimensional parameter-space study, yet the present work only samples three discrete points; a brief remark on how the chosen (φ*, P) values sit inside that larger survey would help the reader place the results.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: radiated-energy ranking is an independent numerical outcome, not forced by construction or self-citation.

full rationale

The paper constructs Bowen-York-type initial data for boson stars by solving the standard conformal-transverse-traceless constraints (Eqs. 10-11) with the known analytic Bowen extended-source extrinsic curvature (Eq. 13) and stationary mini-boson-star profiles. Superposition (Eqs. 36-39) is performed only for approximately compact sources; the Hamiltonian constraint is then solved for the conformal factor. Subsequent BSSN evolutions produce radiated energies, remnant masses, and QNM parameters (Tables II-III) that are compared against external literature (Palenzuela et al., Ge et al., Marks et al.) using different initial-data methods. No parameter is fitted to the target radiation hierarchy and then re-presented as a prediction; the single-star test (Figs. 1-4) merely recovers the known O(P^{2}) ADM-mass growth and the known normal-mode frequencies of the mini-BS. The sole self-citation (Clark & Laguna 2016) is methodological background for the neutron-star analogue and is not load-bearing for the BBS/BHBS energy ranking. The derivation chain is therefore self-contained against external benchmarks; circularity score is zero.

Assumptions & free parameters 3 free parameters · 4 assumptions · 0 invented entities

The work rests on standard 3+1 GR, the CTT decomposition with conformal flatness and maximal slicing, the mini-boson-star potential, and the Bowen extended-source solution. No new physical entities are postulated; free parameters are ordinary simulation choices (central amplitude, boost).

free parameters (3)
  • central scalar amplitude φ* = 0.02–0.04
    Chooses the isolated boson-star mass and compactness; scanned over three discrete values (0.02, 0.03, 0.04).
  • linear momentum P/M* = 0.1–0.3
    Sets the initial boost; scanned over 0.1, 0.2, 0.3.
  • initial separation d = 80 µ⁻¹
    Fixed at 80 µ⁻¹ for all runs; controls constraint-violation size and merger time.
assumptions (4)
  • domain assumption Conformal flatness (γ̃ij = ηij) and maximal slicing (K = 0)
    Standard CTT simplifications that decouple the momentum constraint; known to induce spurious oscillations (explicitly noted in §IV).
  • domain assumption Mini-boson-star potential V = ½ µ² |Φ|²
    Limits the study to the free massive scalar; solitonic potentials are mentioned but not used.
  • domain assumption Bowen extended-source solution for Ãij remains valid when the source is the scalar momentum density
    Derived in §III.B by matching σ to Re(Π* ∇Φ); assumes the conformal rescalings preserve the form of the momentum constraint.
  • ad hoc to paper Superposition of two isolated solutions plus analytic Ãij yields usable initial data after solving only the Hamiltonian constraint
    Matter sources are non-linear; the paper asserts that approximate compact support makes the error tolerable (§V).

how reviews work

0 comments
Cite this review

Pith. "Pith review of Head-on Collisions of Boson Stars with Bowen-York Type Initial Data." pith.science (2026). https://pith.science/paper/5QBKP4A3

@misc{pith2026260709494,
  author       = {Pith},
  title        = {Pith review of: Head-on Collisions of Boson Stars with Bowen-York Type Initial Data},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/5QBKP4A3}},
  note         = {Machine review of arXiv:2607.09494}
}
read the original abstract

We present a numerical relativity study of head-on collisions involving boson stars using initial data inspired by the Bowen-York initial data used to model black hole binaries with punctures. The initial data method preserves the simplicity of the Bowen-York approach, thus allowing incorporating information from the early, post-Newtonian inspiral phase in binary coalescences. We test the method on a single boson star with linear momentum. We present results from head-on collisions of boson stars as well as encounters of boson stars with black holes. In general, the results are consistent with previous studies, demonstrating the effectiveness of the initial data method. In particular, we show that boson star head-on collisions emit more energy in gravitational waves than the equivalent black hole binaries. On the other hand, head-on collisions of a boson star with a black hole radiate less than their black hole binary counterparts.

Figures

Figures reproduced from arXiv: 2607.09494 by the authors.

Figure 1
Figure 1. FIG. 1. ADM mass [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Relative differences with respect to the stationary [PITH_FULL_IMAGE:figures/full_fig_p004_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. For a single BS with linear momentum, top panels depict, as a function of [PITH_FULL_IMAGE:figures/full_fig_p005_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Fourier transform of ( [PITH_FULL_IMAGE:figures/full_fig_p005_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Mode [PITH_FULL_IMAGE:figures/full_fig_p006_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Top panel depicts the maximum of the scalar field [PITH_FULL_IMAGE:figures/full_fig_p007_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Mode [PITH_FULL_IMAGE:figures/full_fig_p008_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

43 extracted references · 2 canonical work pages

  1. [1]

    Rev.97(2) 511–536 URLhttps: //link.aps.org/doi/10.1103/PhysRev.97.511

    Wheeler J A 1955Phys. Rev.97(2) 511–536 URLhttps: //link.aps.org/doi/10.1103/PhysRev.97.511

  2. [2]

    Rev.172(5) 1331–1342 URL https://link.aps.org/doi/10.1103/PhysRev.172

    Kaup D J 1968Phys. Rev.172(5) 1331–1342 URL https://link.aps.org/doi/10.1103/PhysRev.172. 9 1331

  3. [3]

    Seidel E and Suen W M 1990Phys. Rev. D 42(2) 384–403 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.42.384

  4. [4]

    Kain B 2021Phys. Rev. D103(12) 123003 URLhttps: //link.aps.org/doi/10.1103/PhysRevD.103.123003

  5. [5]

    Hawley S H and Choptuik M W 2000Phys. Rev. D62(10) 104024 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.62.104024

  6. [6]

    thesis University of British Columbia URLhttps://open.library.ubc.ca/ collections/ubctheses/831/items/1.0091568

    Rousseau B 2003Axisymmetric boson stars in the con- formally flat approximationPh.D. thesis University of British Columbia URLhttps://open.library.ubc.ca/ collections/ubctheses/831/items/1.0091568

  7. [7]

    thesis University of British Columbia URLhttps://open.library.ubc.ca/collections/ ubctheses/831/items/1.0085582

    Lai C W K 2004A numerical study of boson starsPh.D. thesis University of British Columbia URLhttps://open.library.ubc.ca/collections/ ubctheses/831/items/1.0085582

  8. [9]

    semanticscholar.org/CorpusID:195316543

    Alcubierre M, Barranco J, Bernal A, Degollado J C, Diez- Tejedor A, Megevand M, N’unez D and Sarbach O 2019 Classical and Quantum Gravity36URLhttps://api. semanticscholar.org/CorpusID:195316543

Show all 43 references
  1. [10]

    Palenzuela C, Olabarrieta I, Lehner L and Liebling S L 2007Phys. Rev. D75064005 (Preprintgr-qc/0612067)

  2. [11]

    Choptuik M W and Pretorius F 2010Phys. Rev. Lett. 104(11) 111101 URLhttps://link.aps.org/doi/10. 1103/PhysRevLett.104.111101

  3. [12]

    Mundim B C 2010 A numerical study of boson star bina- ries (Preprint1003.0239) URLhttps://arxiv.org/abs/ 1003.0239

  4. [13]

    Helfer T, Sperhake U, Croft R, Radia M, Ge B X and Lim E A 2022Classical and Quantum Gravity39074001

  5. [14]

    Evstafyeva T, Sperhake U, Helfer T, Croft R, Radia M, Ge B X and Lim E A 2023Class. Quant. Grav.40085009 (Preprint2212.08023)

  6. [15]

    Jaramillo V, Sanchis-Gual N, Barranco J, Bernal A, De- gollado J C, Herdeiro C, Megevand M and N´ u˜ nez D 2022 Phys. Rev. D105(10) 104057 URLhttps://link.aps. org/doi/10.1103/PhysRevD.105.104057

  7. [16]

    Atteneder F, R¨ uter H R, Cors D, Rosca-Mead R, Hilditch D and Br¨ ugmann B 2024Phys. Rev. D 109(4) 044058 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.109.044058

  8. [17]

    Ge B X, Lim E A, Sperhake U, Evstafyeva T, Cors D, de Jong E, Croft R and Helfer T 2025Phys. Rev. D 112(12) 124080 URLhttps://link.aps.org/doi/10. 1103/2dhs-phl4

  9. [18]

    Brito M, Herdeiro C, Radu E, Sanchis-Gual N and Zilh˜ ao M 2026Physical Review D113ISSN 2470-0029 URL http://dx.doi.org/10.1103/8bwz-wqhp

  10. [19]

    Palloni G, Sanchis-Gual N, Font J A, Santos-P´ erez S, Cordero-Carri´ on I, Cerd´ a-Dur´ an P and Lazarte C 2026 Constraint-satisfying binary boson star initial data via xcfc (Preprint2605.20888) URLhttps://arxiv.org/ abs/2605.20888

  11. [20]

    Palenzuela C, Lehner L and Liebling S L 2008Phys. Rev. D77044036 (Preprint0706.2435)

  12. [21]

    Palenzuela C, Pani P, Bezares M, Cardoso V, Lehner L and Liebling S 2017Phys. Rev. D96(10) 104058 URLhttps://link.aps.org/doi/10.1103/PhysRevD. 96.104058

  13. [22]

    Bezares M, Palenzuela C and Bona C 2017Phys. Rev. D95(12) 124005 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.95.124005

  14. [23]

    Bezares M, Boˇ skovi´ c M, Liebling S, Palenzuela C, Pani P and Barausse E 2022Phys. Rev. D 105(6) 064067 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.105.064067

  15. [24]

    Siemonsen N and East W E 2023Phys. Rev. D 107(12) 124018 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.107.124018

  16. [25]

    1088/1361-6382/acace4

    Croft R, Helfer T, Ge B X, Radia M, Evstafyeva T, Lim E A, Sperhake U and Clough K 2023Classical and Quantum Gravity40065001 URLhttps://doi.org/10. 1088/1361-6382/acace4

  17. [26]

    Siemonsen N and East W E 2023Phys. Rev. D 108(12) 124015 URLhttps://link.aps.org/doi/10. 1103/PhysRevD.108.124015

  18. [27]

    Evstafyeva T, Seifert A, Sperhake U, Moore C J and Jain T 2026 Lessons from binary dynamics of inspiralling equal-mass boson-star mergers (Preprint2604.25582) URLhttps://arxiv.org/abs/2604.25582

  19. [28]

    1088/1361-6382/aaf43e

    Dietrich T, Ossokine S and Clough K 2018Classical and Quantum Gravity36025002 URLhttps://doi.org/10. 1088/1361-6382/aaf43e

  20. [29]

    Cardoso V, Ikeda T, Zhong Z and Zilh˜ ao M 2022Phys. Rev. D106(4) 044030 URLhttps://link.aps.org/ doi/10.1103/PhysRevD.106.044030

  21. [30]

    Zhong Z, Cardoso V, Ikeda T and Zilh˜ ao M 2023Phys. Rev. D108(8) 084051 URLhttps://link.aps.org/ doi/10.1103/PhysRevD.108.084051

  22. [31]

    org/abs/2604.06312

    Marks G A, Staelens S J and Sperhake U 2026 Black hole- boson star binaries: Gravitational wave signals and tidal disruption (Preprint2604.06312) URLhttps://arxiv. org/abs/2604.06312

  23. [32]

    Ning Z 2026 Boson star-black hole binaries: initial data and head-on collisions (Preprint2604.15240) URL https://arxiv.org/abs/2604.15240

  24. [33]

    Goodale T, Allen G, Lanfermann G, Masso J, Radke T, Seidel E and Shalf J 2003 The Cactus Framework and Toolkit: Design and ApplicationsVECPAVector and Parallel Processing R’2002, 5th International Conference (Berlin: Springer)

  25. [34]

    Husa S, Hinder I and Lechner C 2006Comput. Phys. Commun.174983–1004 (Preprintgr-qc/0404023)

  26. [35]

    Haas R, Shcherbakov R V, Bode T and Laguna P 2012 The Astrophysical Journal749117 ISSN 0004-637X

  27. [36]

    Evans C, Laguna P and Eracleous M 2015The Astro- physical Journal Letters805L19 ISSN 2041-8205

  28. [37]

    Clark M and Laguna P 2016Physical Review D94 064058

  29. [38]

    Jani K, Healy J, Clark J A, London L, Laguna P and Shoemaker D 2016Classical and Quantum Gravity33 204001 ISSN 0264-9381

  30. [39]

    Rizzo M, Haas R, Brandt S R, Etienne Z, Ferguson D, Sanches L T, Tsao B J, Werneck L, Boyer D, Bozzola G, Cheng C H, Cupp S, Diener P, Jacques T P, Ji L, Macpherson H, Markin I, Schnetter E, Tichy W, Tootle S, Xu Y, Zilh˜ ao M, Zlochower Y, Alcubierre M, Alic D, Allen G, Ansor...

  31. [40]

    Baumgarte T W and Shapiro S L 2010Numerical Rel- ativity: Solving Einstein’s Equations on the Computer (Cambridge University Press)

  32. [41]

    Smarr L L (ed) 1979Proceedings, Sources of Gravita- tional Radiation: Seattle, WA, USA, July 24 - August 4, 1978(Cambridge: Cambridge Univ. Press)

  33. [42]

    Bowen J M and York Jr J W 1980Phys. Rev. D21 2047–2056

  34. [43]

    1007/BF00762132

    Bowen J M 1979General Relativity and Gravitation 11227–231 ISSN 1572-9532 URLhttps://doi.org/10. 1007/BF00762132

  35. [44]

    Clark M and Laguna P 2016Phys. Rev. D94064058 (Preprint1606.04881)

Pith tools

Reviewed July 13, 2026 · model on record in the stance chip above.