Pith. sign in

REVIEW 2 major objections 2 minor 38 references

Precise scaling relations for self-interacting bosonic dark matter stars

T0 review · 2 major / 2 minor · reviewed 2026-06-27 · grok-4.3

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

Pith's one-line read Bosonic dark matter stars follow precise scaling relations for their maximum mass, radius, and central density based on the boson mass and self-coupling.

desk verdict The paper fits universal scaling relations for quartic boson stars with small errors, but adds little beyond determining known constants from the scaling symmetry. read the letter →

arxiv 2606.08967 v1 pith:LHPEVTRA submitted 2026-06-08 astro-ph.HE astro-ph.COastro-ph.GA

Zhaorui Zhang , Xudong Wang , Bin Qi This is my paper · ORCID
classification astro-ph.HEastro-ph.COastro-ph.GA
keywords bosonicdarkmatterstarsself-interactingscalarfieldscalingrelationsmaximummassquarticpotentialstellarstructureequationofstate
topics Dark Matter
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

The paper systematically studies the properties of stars made from bosonic dark matter with self-interactions described by a quartic potential. Numerical solutions of the structure equations across a broad range of particle masses and couplings yield scaling relations for the maximum mass, critical radius, and critical central density. These relations are presented with explicit prefactors and hold with relative errors below 4 percent. The authors also supply unified analytical fits for how mass and radius vary with central density on the stable branch, accurate to 0.1 percent. This approach replaces repeated numerical integrations with direct formulas that depend only on the dark matter particle parameters.

What carries the argument

Scaling relations derived from numerical integration of the stellar structure equations using the equation of state for a complex scalar field with quartic self-interaction potential.

What would settle it

A set of numerical stellar models at a boson mass and coupling value within the studied ranges whose maximum mass deviates from the predicted scaling by more than 4 percent would falsify the claimed precision of the relations.

Watch

Extended reading notes

Core claim

Numerical solutions of the stellar structure equations with the equation of state from a complex scalar field with quartic potential produce the scaling relations M_max = 0.1 sqrt(λ)/m_φ² solar masses, R(M_max) = 0.9 sqrt(λ)/m_φ² km, and ε_max = 2.1×10^5 m_φ⁴/λ MeV/fm³ where m_φ is in GeV, with fitting relative error less than 4 percent. The mass-central density and radius-central density relations on the stable branch are described by a single functional form with parameters chosen separately for mass and radius, achieving fitting relative error less than 0.1 percent. A simple quadratic polynomial mass-radius relation is also identified.

Load-bearing premise

The numerical results produce stable prefactors in the scaling relations that can be captured by the chosen functional forms without additional systematic effects across the parameter ranges examined.

Editorial extensions

If this is right

  • Maximum mass and size of bosonic dark matter stars can be calculated directly from the boson mass and coupling without solving the differential equations each time.
  • The stable branch configurations collapse onto universal curves when expressed in terms of the critical values.
  • The mass-radius relation takes a simple quadratic form that follows from the underlying equation of state.
  • These properties hold uniformly for boson masses spanning twelve orders of magnitude and couplings from 0.01π to 100π.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

  • The scaling relations could constrain the allowed range of dark matter particle parameters if any bosonic stars are observed.
  • The unified fitting function might apply to other interaction potentials or be derived analytically in limiting cases.
  • Comparison with the mass-radius relations of neutron stars or other exotic compact objects could distinguish bosonic dark matter stars observationally.
Share X Bluesky LinkedIn Reddit HN

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

2 major / 2 minor

Summary. The manuscript numerically solves the Tolman-Oppenheimer-Volkoff equations using an equation of state derived from a complex scalar field with quartic potential V(φ)=λ/4|φ|^4. It reports scaling relations M_max=0.1 sqrt(λ)/m_φ² M_⊙, R(M_max)=0.9 sqrt(λ)/m_φ² km, and ε_max=2.1×10^5 m_φ⁴/λ MeV/fm³ (with relative fit error <4%) that hold across the scanned ranges 10^{-9} GeV ≤ m_φ ≤ 10^3 GeV and 0.01π ≤ λ ≤ 100π. Global analytical fits of the stable branch are given via a unified function Ỹ=A/[1+(5ε̃)^h]^s for mass-central density and radius-central density (relative error <0.1%), together with a quadratic mass-radius relation.

Significance. The quartic potential admits an exact scaling symmetry that renders the dimensionless maximum mass, radius, and central density universal constants independent of m_φ and λ. Accurate numerical determination of these constants therefore supplies ready-to-use formulas that eliminate the need to re-integrate the structure equations for each particle-physics parameter choice. The work thereby supplies a practical tool for rapid estimates of bosonic dark-matter star properties in the self-interacting regime.

major comments (2)
  1. [Abstract, §3] Abstract and §3 (numerical procedure): the quoted prefactors 0.1, 0.9 and 2.1×10^5 are obtained from numerical integration, yet the manuscript supplies no information on the integration scheme, radial grid resolution, convergence tests, or cross-checks against independent codes. Because these prefactors are the central quantitative results, the absence of such documentation prevents independent verification of the stated <4% fitting error.
  2. [Abstract] Abstract, unified-function paragraph: the specific exponents h=-2 (mass) and h=1 (radius) together with amplitudes A=1 and A=1.634 are presented as empirical fits. No derivation or physical motivation is given for these functional choices, nor is it shown that alternative forms (e.g., polytropic or Lane-Emden inspired) yield comparable or worse residuals across the full range of ε̃.
minor comments (2)
  1. [Abstract] Notation: the definition ε̃ ≡ ε_0/ε_max is introduced only in the unified-function paragraph; an explicit statement earlier in the text would improve readability.
  2. [Abstract] The mass-radius quadratic polynomial is mentioned but neither its coefficients nor its fitting domain are stated; these should be supplied explicitly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We thank the referee for the constructive report and positive recommendation. We address each major comment below. The revisions will focus on adding the requested numerical documentation and functional-form justification without altering the core results.

read point-by-point responses

  1. Referee: [Abstract, §3] Abstract and §3 (numerical procedure): the quoted prefactors 0.1, 0.9 and 2.1×10^5 are obtained from numerical integration, yet the manuscript supplies no information on the integration scheme, radial grid resolution, convergence tests, or cross-checks against independent codes. Because these prefactors are the central quantitative results, the absence of such documentation prevents independent verification of the stated <4% fitting error.

    Authors: We agree that additional documentation of the numerical methods is required for reproducibility. In the revised manuscript we will expand §3 with a new subsection detailing the integration scheme (fourth-order Runge-Kutta with adaptive step-size control), the radial grid (typically 5000–10000 points with adaptive refinement near the surface), explicit convergence tests (results stable to <0.1% when resolution is doubled), and cross-checks against the λ=0 analytic limit and published boson-star codes. These additions will directly support the quoted fit accuracy. revision: yes

  2. Referee: [Abstract] Abstract, unified-function paragraph: the specific exponents h=-2 (mass) and h=1 (radius) together with amplitudes A=1 and A=1.634 are presented as empirical fits. No derivation or physical motivation is given for these functional choices, nor is it shown that alternative forms (e.g., polytropic or Lane-Emden inspired) yield comparable or worse residuals across the full range of ε̃.

    Authors: The exponents and amplitudes were selected after systematic trials because they simultaneously reproduce the low-density asymptotic scaling (M ∝ ε̃ for small ε̃) and the high-density behavior while keeping the functional form compact. In the revision we will add a paragraph explaining this rationale and include a short comparison demonstrating that a polytropic-inspired power-law alternative produces residuals >5% over parts of the range, whereas the chosen form stays below 0.1%. revision: yes

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity identified

full rationale

The paper derives the EOS from the complex scalar with quartic potential, solves the stellar structure equations numerically over the stated ranges of m_φ and λ, and reports the resulting maxima and curves. The scaling symmetry (ξ = m_φ r, σ = √λ φ / m_φ) renders all dimensionless combinations independent of the parameters, so the quoted prefactors (0.1, 0.9, 2.1×10^5) are simply the numerically determined values of the universal constants M_max m_φ²/√λ etc.; the <4% and <0.1% errors quantify the quality of the subsequent analytic fits to those computed curves. No step equates a claimed prediction to its own input by construction, no self-citation is load-bearing, and the central results remain independent numerical outputs rather than tautological re-expressions of fitted parameters.

Assumptions & free parameters 6 free parameters · 2 assumptions · 0 invented entities

The scaling relations rest on fitted numerical coefficients and the choice of quartic potential; no independent evidence is supplied for the model beyond the internal numerical consistency.

free parameters (6)
  • M_max prefactor = 0.1
    Numerical coefficient 0.1 fitted to maximum-mass data
  • R(M_max) prefactor = 0.9
    Numerical coefficient 0.9 fitted to critical-radius data
  • ε_max prefactor = 2.1e5
    Numerical coefficient 2.1e5 fitted to critical-density data
  • Unified R amplitude A = 1.634
    Fitted amplitude 1.634 in the radius unified function
  • Unified R exponent s = 0.28
    Fitted exponent 0.28 in the radius unified function
  • Unified M exponent s = 0.42
    Fitted exponent 0.42 in the mass unified function
assumptions (2)
  • domain assumption Bosonic dark matter is described by a complex scalar field with quartic self-interaction potential V(φ) = λ/4 |φ|^4
    This potential is used to derive the equation of state for the stars.
  • standard math Stellar structure obeys the Tolman-Oppenheimer-Volkoff equation of general relativity
    Standard assumption invoked for compact-object modeling.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Precise scaling relations for self-interacting bosonic dark matter stars." pith.science (2026). https://pith.science/paper/LHPEVTRA

@misc{pith2026260608967,
  author       = {Pith},
  title        = {Pith review of: Precise scaling relations for self-interacting bosonic dark matter stars},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/LHPEVTRA}},
  note         = {Machine review of arXiv:2606.08967}
}
abstract

The structural properties of bosonic dark matter stars are systematically investigated, presenting precise scaling relations for the mass, radius, central density, and the properties of dark matter particles. The dark matter equation of state is derived from a complex scalar field theory with a quartic self-interaction potential $V(\phi) = \frac{\lambda}{4} |\phi|^4$, considering boson masses $m_{\phi}$ ranging from $10^{-9}$ to $10^{3}$ GeV and self-coupling constants $\lambda$ ranging from $0.01\pi$ to $100\pi$. The scaling relation for the maximum mass of bosonic dark matter stars, the corresponding critical radius and critical central density are obtained as \[ M_{\text{max}} = 0.1 \frac{\sqrt{\lambda}}{m_\phi^2} M_\odot, \qquad R(M_{\text{max}}) = 0.9 \frac{\sqrt{\lambda}}{m_\phi^2} \ \text{km}, \qquad \varepsilon_{\text{max}} = 2.1 \times 10^5 \frac{m_\phi^4}{\lambda} \ \mathrm{MeV/fm^3}, \] where $m_\phi$ is in GeV, the relations for $R(M_{\text{max}})$ and $\varepsilon_{\text{max}}$ are first put forward. The fitting relative error is less than $4\%$. Based on these scaling relations, we further provide global analytical fits for the stable branch. The relationships between mass and central density as well as radius and central density can be described by a unified function of the form: \[ \tilde{Y} = \frac{A}{\left[1 + \left(5\tilde{\varepsilon}\right)^h\right]^s}, \] where for $Y=M$, $\tilde{M} \equiv M/M_{\text{max}}$, $A=1$, $h=-2$, $s=0.42$; for $Y=R$, $\tilde{R} \equiv R/R(M_{\text{max}})$, $A=1.634$, $h=1$, $s=0.28$; and $\tilde{\varepsilon} \equiv \varepsilon_0/\varepsilon_{\text{max}}$. The fitting relative error is less than $0.1\%$. Furthermore, we find a simple quadratic polynomial mass-radius relation for bosonic dark matter stars.

Figures

Figures reproduced from arXiv: 2606.08967 by the authors.

Figure 1
Figure 1. FIG. 1. Equation of state for self-interacting bosonic dark matter with a quartic potential. The figure shows the pressure [PITH_FULL_IMAGE:figures/full_fig_p004_1.png] view at source ↗
Figure 2
Figure 2. FIG. 2. Fixed [PITH_FULL_IMAGE:figures/full_fig_p005_2.png] view at source ↗
Figure 3
Figure 3. FIG. 3. Fixed self-coupling constant [PITH_FULL_IMAGE:figures/full_fig_p006_3.png] view at source ↗
Figures from the paper (7 more)
Figure 4
Figure 4. Figure 4: FIG. 4. Fixed mass [PITH_FULL_IMAGE:figures/full_fig_p007_4.png]
Figure 5
Figure 5. Figure 5: FIG. 5. Fixed mass [PITH_FULL_IMAGE:figures/full_fig_p008_5.png]
Figure 6
Figure 6. Figure 6: FIG. 6. Fixed self-coupling constant [PITH_FULL_IMAGE:figures/full_fig_p008_6.png]
Figure 7
Figure 7. Figure 7: FIG. 7. Fixed self-coupling constant [PITH_FULL_IMAGE:figures/full_fig_p009_7.png]
Figure 8
Figure 8. Figure 8: FIG. 8. Dimensionless mass [PITH_FULL_IMAGE:figures/full_fig_p011_8.png]
Figure 9
Figure 9. Figure 9: FIG. 9. Dimensionless mass-radius relation on the stable branch. Scattered points are numerically calculated results from the [PITH_FULL_IMAGE:figures/full_fig_p012_9.png]
Figure 10
Figure 10. Figure 10: FIG. 10. Allowed parameter space for Sgr A* as a bosonic dark matter star. The blue solid line is the lower bound from [PITH_FULL_IMAGE:figures/full_fig_p014_10.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

38 extracted references · 1 canonical work pages

  1. [1]

    Scaling relation for characteristic density:ε max ∝m 4 ϕ/λ, together withY max ∝ √ λ/m2 ϕ, reflects the self- similarity of the system

  2. [2]

    Difference in shape parameters: for the mass fith=−2<0 while for the radius fith= 1>0, indicating that the contraction behavior of radius with central density and the growth behavior of mass have different functional characteristics, consistent with the trends observed in Figs. 2 and 4

  3. [3]

    The radius varies by a factor of about 1.6, quantitatively characterizing the evolution of the boson star radius with increasing central density

    Variation of radius within the stable branch: for the radius fit, within the stable branch (ε 0 ≤ε max), when ε0 ≪ε max, ˜R≈A= 1.634, i.e.,R≈1.634R(M max); at the critical point, ˜R= 1, i.e.,R=R(M max). The radius varies by a factor of about 1.6, quantitatively characterizing the evolution of the boson star radius with increasing central density. B. Unive...

  4. [4]

    Bertone, D

    G. Bertone, D. Hooper, and J. Silk,Phys. Rep.405, 279 (2005)

    2005

  5. [5]

    Arcadi, M

    G. Arcadi, M. Dutra, P. Ghosh, M. Lindner, Y. Mambrini, M. Pierre, S. Profumo, and F. S. Queiroz,Eur. Phys. J. C78, 203 (2018)

    2018

  6. [6]

    Schumann,J

    M. Schumann,J. Phys. G46, 103003 (2019). 15

    2019

  7. [7]

    Clowe, M

    D. Clowe, M. Bradaˇ c, A. H. Gonzalez, M. Markevitch, S. W. Randall, C. Jones, and D. Zaritsky,Astrophys. J. Lett.648, L109 (2006)

    2006

  8. [8]

    Aghanimet al.,Astron

    Planck Collaboration, N. Aghanimet al.,Astron. Astrophys.641, A6 (2020)

    2020

Show all 38 references
  1. [9]

    D. S. Akeribet al.(LUX Collaboration),Phys. Rev. Lett.118, 021303 (2017)

    2017

  2. [10]

    Aprileet al.(XENON Collaboration),Phys

    E. Aprileet al.(XENON Collaboration),Phys. Rev. Lett.121, 111302 (2018)

    2018

  3. [11]

    D. J. E. Marsh,Phys. Rep.643, 1 (2016)

    2016

  4. [12]

    L. Hui, J. P. Ostriker, S. Tremaine, and E. Witten,Phys. Rev. D95, 043541 (2017)

    2017

  5. [13]

    A. M. Green and B. J. Kavanagh,J. Phys. G48, 043001 (2021)

    2021

  6. [14]

    D. J. Kaup,Phys. Rev.172, 1331 (1968)

    1968

  7. [15]

    Ruffini and S

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

    1969

  8. [16]

    J. A. Wheeler,Phys. Rev.97, 511 (1955)

    1955

  9. [17]

    Chandrasekhar,Astrophys

    S. Chandrasekhar,Astrophys. J.74, 81 (1931)

    1931

  10. [18]

    F. E. Schunck and E. W. Mielke,Class. Quantum Grav.20, R301 (2003)

    2003

  11. [19]

    Jetzer,Phys

    P. Jetzer,Phys. Rep.220, 163 (1992)

    1992

  12. [20]

    Colpi, S

    M. Colpi, S. L. Shapiro, and I. Wasserman,Phys. Rev. Lett.57, 2485 (1986)

    1986

  13. [21]

    S. L. Liebling and C. Palenzuela,Living Rev. Rel.15, 6 (2012)

    2012

  14. [22]

    Olivares, Z

    H. Olivares, Z. Younsi, C. M. Fromm, M. De Laurentis, O. Porth, Y. Mizuno, H. Falcke, M. Kramer, and L. Rezzolla, Mon. Not. R. Astron. Soc.497, 521 (2020)

    2020

  15. [23]

    F. H. Vincent, M. Wielgus, M. A. Abramowicz, E. Gourgoulhon, J.-P. Lasota, T. Paumard, and G. Perrin,Astron. Astrophys.646, A37 (2021)

    2021

  16. [24]

    Jaramillo and S.-Y

    V. Jaramillo and S.-Y. Zhou,Complex structures of boson stars and anisotropic distribution of satellite galaxies, arXiv:2407.12084 (2024)

    2024

  17. [25]

    P. H. Chavanis,Phys. Rev. D107, 103503 (2023)

    2023

  18. [26]

    P. H. Chavanis and T. Harko,Phys. Rev. D86, 064011 (2012)

    2012

  19. [27]

    R. F. Tooper,Astrophys. J.140, 434 (1964)

    1964

  20. [28]

    J. R. Oppenheimer and G. M. Volkoff,Phys. Rev.55, 374 (1939)

    1939

  21. [29]

    B. K. Harrison, K. S. Thorne, M. Wakano, and J. A. Wheeler,Gravitation Theory and Gravitational Collapse(University of Chicago Press, 1965)

    1965

  22. [30]

    J. M. Bardeen, K. S. Thorne, and D. W. Meltzer,Astrophys. J.145, 505 (1966)

    1966

  23. [31]

    J. M. Lattimer and M. Prakash,Astrophys. J.550, 426 (2001)

    2001

  24. [32]

    J. M. Lattimer and M. Prakash,Science304, 536 (2004)

    2004

  25. [33]

    J. M. Lattimer and M. Prakash,Phys. Rep.442, 109 (2007)

    2007

  26. [34]

    Yagi and N

    K. Yagi and N. Yunes,Phys. Rep.681, 1 (2017)

    2017

  27. [35]

    Abuteret al.,Astron

    GRAVITY Collaboration, R. Abuteret al.,Astron. Astrophys.657, L12 (2022)

    2022

  28. [36]

    Akiyamaet al.,Astrophys

    Event Horizon Telescope Collaboration, K. Akiyamaet al.,Astrophys. J. Lett.930, L12 (2022)

    2022

  29. [37]

    Akiyamaet al.,Astrophys

    Event Horizon Telescope Collaboration, K. Akiyamaet al.,Astrophys. J. Lett.930, L13 (2022)

    2022

  30. [38]

    Agnihotri, J

    P. Agnihotri, J. Schaffner-Bielich, and I. N. Mishustin,Phys. Rev. D79, 084033 (2009)

    2009

Pith tools

Reviewed June 27, 2026 · model on record in the stance chip above.