arxiv: 2604.13011 · v2 · submitted 2026-04-14 · 🌀 gr-qc

Recognition: unknown

Buchdahl Limit and TOV Equations in Interacting Vacuum Scenarios

Rodrigo Maier

Authors on Pith no claims yet

Pith reviewed 2026-05-10 15:21 UTC · model grok-4.3

classification 🌀 gr-qc

keywords interacting vacuumBuchdahl limitTOV equationsultra-compact starsstellar stabilitygeneral relativityenergy exchangecovariant interaction

0 comments

The pith

Interacting vacuum keeps central pressure finite in stars beyond the Buchdahl compactness limit.

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

The paper modifies the Tolman-Oppenheimer-Volkoff equations to account for energy exchange between matter and an interacting vacuum. Two phenomenological coupling models are considered: one to the matter density gradient and one to spacetime curvature. Numerical solutions demonstrate that the interaction can prevent the central pressure from diverging as the stellar compactness approaches or exceeds the classical Buchdahl bound. This indicates that vacuum interactions offer a way to support stable ultra-compact configurations that would be singular in standard general relativity.

Core claim

Extending the Tolman-Oppenheimer-Volkoff hydrostatic equilibrium equation with a covariant interaction four-vector Q_ν between the fluid and vacuum energy-momentum tensors allows the central pressure to remain finite for compactness parameters that exceed the Buchdahl limit, for appropriate values of the coupling strength in both the density-gradient and curvature-coupled models.

What carries the argument

The interaction term Q_ν, which represents covariant energy exchange and modifies the pressure gradient in the extended TOV equations.

If this is right

Central pressure stays finite where standard GR predicts divergence at the Buchdahl threshold.
Ultra-compact stellar objects become possible without encountering singularities.
Both coupling to matter energy-density gradient and to curvature permit stable solutions in suitable parameter ranges.
The classical geometric bound on compactness can be bypassed through vacuum-fluid interactions.

Where Pith is reading between the lines

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

Such models might allow exotic compact objects that mimic black-hole exteriors without horizons.
Mass-radius relations for neutron stars could shift measurably under these interactions.
Gravitational-wave signals from mergers might carry imprints of the modified internal pressure profile.

Load-bearing premise

The phenomenological interaction terms chosen for the vacuum-fluid coupling yield stable and physically acceptable solutions for some values of the coupling parameter.

What would settle it

Numerical integration of the extended TOV equations with either proposed form of Q_ν showing that central pressure still diverges exactly at the Buchdahl compactness for every coupling value.

Figures

Figures reproduced from arXiv: 2604.13011 by Rodrigo Maier.

**Figure 2.** Figure 2: FIG. 2: Numerical evolution of the internal pressure profile [PITH_FULL_IMAGE:figures/full_fig_p006_2.png] view at source ↗

**Figure 3.** Figure 3: FIG. 3: A diagram illustrating the modified compactness limi [PITH_FULL_IMAGE:figures/full_fig_p007_3.png] view at source ↗

**Figure 4.** Figure 4: FIG. 4: Numerical evolution of the internal pressure profile [PITH_FULL_IMAGE:figures/full_fig_p008_4.png] view at source ↗

read the original abstract

We investigate the stability of ultra-compact stellar configurations in the context of an interacting vacuum component. By extending the Tolman-Oppenheimer-Volkoff equations to include a covariant energy exchange between the fluid and vacuum sectors, we examine how the classical Buchdahl stability limit is modified. We analyze two phenomenological interaction models: a coupling to the matter energy density gradient and a direct coupling to the spacetime curvature. Numerical integration reveals that while standard General Relativity predicts a central pressure divergence as the compactness approaches the Buchdahl threshold, the interaction term $Q_\nu$ relaxes the pressure gradient and maintains a finite, well-behaved central pressure for proper domains of the coupling parameter. These results demonstrate that an interacting vacuum provides a physical mechanism to bypass classical geometric bounds, potentially supporting ultra-compact objects in regimes previously considered singular.

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 numerically shows that two ad-hoc vacuum-matter interaction terms can keep central pressure finite past the Buchdahl limit in the TOV system, but the effect is tied to those specific choices.

read the letter

The paper shows that adding two phenomenological interaction terms to the vacuum-fluid system in the TOV equations can keep central pressure finite for stars more compact than the Buchdahl limit. Numerical integration for couplings to the density gradient and to curvature demonstrates the pressure gradient relaxing in the right parameter range. This applies the interacting vacuum concept to stellar hydrostatics in a direct way. The work does the extension of the equations and carries out the integrations to produce explicit examples where the classical divergence is avoided. The calculation is clean on its own terms and the result follows from the extra source term in the conservation law. It illustrates how breaking the separate conservation of the matter tensor changes the geometric bound. The soft spots are the lack of fundamental motivation for the specific forms of Q_ν. Both are introduced as phenomenological choices, so the bypass may not generalize beyond these cases. It would be good to see explicit checks that the solutions obey the weak energy condition and remain stable, along with any convergence tests on the numerical scheme. This paper is aimed at researchers working on compact objects in modified gravity or with extra components. A reader looking for concrete calculations on how interactions affect stellar bounds will get value from the numerical results. I think it deserves peer review. The claim is testable with the provided equations and numerics, and referees can assess whether the physical requirements are met.

Referee Report

2 major / 1 minor

Summary. The manuscript extends the Tolman-Oppenheimer-Volkoff (TOV) equations to incorporate covariant energy exchange between a perfect fluid and an interacting vacuum sector. Two phenomenological interaction terms Q_ν are considered—one proportional to the gradient of the matter energy density and the other to the spacetime curvature scalar. Numerical integration of the modified hydrostatic equilibrium equation demonstrates that, for suitable ranges of the coupling strength, the central pressure remains finite even when the compactness parameter exceeds the classical Buchdahl bound of 8/9, in contrast to the divergence obtained in standard general relativity.

Significance. If the numerical results are robust, the work identifies a concrete mechanism by which an interacting vacuum can relax the hydrostatic pressure gradient and thereby evade the geometric Buchdahl limit. This could be relevant for modeling ultra-compact objects or for exploring the role of vacuum-matter coupling in strong-field regimes. The explicit demonstration that finite central pressure is recovered for nonzero coupling is a clear strength, although the phenomenological status of the interaction terms limits the generality of the conclusion.

major comments (2)

[Sections describing the interaction models and the modified TOV system] The central claim that the interaction term relaxes the pressure gradient and yields finite central pressure rests on the two specific phenomenological forms of Q_ν. Because these forms are introduced by hand without derivation from an action principle or from a well-defined limit in which they reduce to standard GR, it is unclear whether the reported bypass is a generic feature of interacting vacuum or an artifact of the chosen source terms. The classical Buchdahl derivation assumes separate conservation of T_μν; the extra Q_ν term directly violates that assumption, so the result is tied to the arbitrary choice rather than to the presence of an interacting vacuum sector per se.
[Numerical results and integration procedure] The abstract and results section state that numerical integration was performed and that central pressure remains finite for proper coupling domains, yet no boundary conditions at the stellar center, integration scheme, step-size control, convergence tests, or error estimates are supplied. Without these, it is impossible to verify that the reported finite-pressure solutions are not numerical artifacts, which directly undermines the load-bearing claim that the Buchdahl limit is bypassed.

minor comments (1)

[Abstract] The abstract refers to 'proper domains of the coupling parameter' without quoting the numerical interval or the criterion used to define acceptability; this should be stated explicitly.

Simulated Author's Rebuttal

2 responses · 0 unresolved

We are grateful to the referee for the insightful comments, which have helped us improve the clarity and robustness of our manuscript. We provide point-by-point responses to the major comments below.

read point-by-point responses

Referee: The central claim that the interaction term relaxes the pressure gradient and yields finite central pressure rests on the two specific phenomenological forms of Q_ν. Because these forms are introduced by hand without derivation from an action principle or from a well-defined limit in which they reduce to standard GR, it is unclear whether the reported bypass is a generic feature of interacting vacuum or an artifact of the chosen source terms. The classical Buchdahl derivation assumes separate conservation of T_μν; the extra Q_ν term directly violates that assumption, so the result is tied to the arbitrary choice rather than to the presence of an interacting vacuum sector per se.

Authors: We acknowledge that the interaction terms Q_ν are phenomenological, as is common in the literature on interacting dark energy and vacuum models. Our work demonstrates that for these covariant interaction forms, the hydrostatic equilibrium is modified such that the pressure gradient is relaxed, allowing finite central pressure beyond the Buchdahl limit. We agree that this does not prove it is generic to all interacting vacuum scenarios, but rather shows a possible mechanism within this class of models. In the revised manuscript, we will expand the discussion to explicitly state the phenomenological nature and the limitations of our conclusions, including that the classical Buchdahl theorem relies on separate conservation which is violated here by design. revision: partial
Referee: The abstract and results section state that numerical integration was performed and that central pressure remains finite for proper coupling domains, yet no boundary conditions at the stellar center, integration scheme, step-size control, convergence tests, or error estimates are supplied. Without these, it is impossible to verify that the reported finite-pressure solutions are not numerical artifacts, which directly undermines the load-bearing claim that the Buchdahl limit is bypassed.

Authors: We thank the referee for pointing out this omission. The manuscript will be revised to include a dedicated subsection on the numerical methods. Specifically, we will detail the central boundary conditions (regularity requiring dp/dr = 0 at r=0, and finite central density and pressure), the integration scheme (fourth-order Runge-Kutta with adaptive step-size control), convergence tests by varying step sizes, and error estimates based on comparison with known GR limits for small coupling. This will ensure the robustness of the finite central pressure solutions. revision: yes

Circularity Check

0 steps flagged

No significant circularity; results follow directly from numerical integration of modified TOV equations with explicit phenomenological source terms.

full rationale

The paper introduces two explicit phenomenological forms for the interaction four-vector Q_ν (coupling to the matter density gradient or to curvature), augments the standard TOV hydrostatic equilibrium equation with the resulting extra source term, and performs numerical integration of the resulting first-order system. The reported finite central pressure at compactness exceeding the classical Buchdahl value is a direct numerical consequence of that added term relaxing the pressure gradient; it is not obtained by fitting any parameter to the target outcome and then relabeling the fit as a prediction. No load-bearing uniqueness theorem, self-citation chain, or ansatz imported from prior work by the same author is invoked to justify the central claim. The derivation therefore remains self-contained: the bypass of the geometric bound is an immediate output of the chosen extended differential equations rather than a reduction to the input assumptions by construction.

Axiom & Free-Parameter Ledger

1 free parameters · 2 axioms · 1 invented entities

The central claim rests on the standard GR derivation of the TOV equations plus two ad-hoc phenomenological interaction terms whose functional forms are chosen rather than derived from a fundamental action.

free parameters (1)

coupling strength parameter
The magnitude of the energy-exchange term Q_ν is varied as a free parameter to identify domains where central pressure remains finite.

axioms (2)

domain assumption The interaction between fluid and vacuum can be modeled covariantly by adding a source term Q_ν to the conservation equations.
Invoked when extending the standard TOV system; no derivation from an action is given in the abstract.
ad hoc to paper The two chosen functional forms of Q_ν (density-gradient coupling and curvature coupling) are representative and physically acceptable.
These specific choices are introduced phenomenologically to explore the effect.

invented entities (1)

interacting vacuum component no independent evidence
purpose: To exchange energy with the stellar fluid and thereby relax the pressure gradient
Introduced as a phenomenological addition; no independent observational evidence is cited in the abstract.

pith-pipeline@v0.9.0 · 5426 in / 1493 out tokens · 57450 ms · 2026-05-10T15:21:03.219842+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

29 extracted references · 24 canonical work pages

[1]

R. C. Tolman, Phys. Rev. 55, 364-373 (1939)

1939
[2]

J. R. Oppenheimer and G. M. Volkoﬀ, Phys. Rev. 55, 374-381 (1939)

1939
[3]

H. A. Buchdahl, Phys. Rev. 116, 1027 (1959)

1959
[4]

B. P. Abbott et al. [LIGO Scientiﬁc and Virgo], Phys. Rev. Lett. 121, no.16, 161101 (2018) [arXiv:1805.11581 [gr-qc]]

work page Pith review arXiv 2018
[5]

M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Ar zoumanian, K. C. Gendreau, S. Guillot, A. K. Harding, W. C. G. Ho and J. M. Lattimer, et al. Astrophys. J. Lett. 887, no.1, L24 (2019) [arXiv:1912.05705 [astro-ph.HE]]. 10

work page arXiv 2019
[6]

M. C. Miller, F. K. Lamb, A. J. Dittmann, S. Bogdanov, Z. Ar zoumanian, K. C. Gendreau, S. Guillot, W. C. G. Ho, J. M. Lattimer and M. Loewenstein, et al. Astrophys. J. Lett. 918, no.2, L28 (2021) [arXiv:2105.06979 [astro-ph.HE]]

work page arXiv 2021
[7]

Weinberg, Rev

S. Weinberg, Rev. Mod. Phys. 61, 1-23 (1989)

1989
[8]

Quintessence, Cosmic Coincidence, and the Cosmological Constant

I. Zlatev, L. M. Wang and P. J. Steinhardt, Phys. Rev. Lett . 82, 896-899 (1999) [arXiv:astro-ph/9807002 [astro-ph]]

work page Pith review arXiv 1999
[9]

Di Valentino, O

E. Di Valentino, O. Mena, S. Pan, L. Visinelli, W. Yang, A. Melchiorri, D. F. Mota, A. G. Riess and J. Silk, Class. Quant. Grav. 38, no.15, 153001 (2021) [arXiv:2103.01183 [astro-ph.CO]]

work page arXiv 2021
[10]

A. V. Astashenok, S. Capozziello and S. D. Odintsov, JCA P 12, 040 (2013) [arXiv:1309.1978 [gr-qc]]

work page arXiv 2013
[11]

Capozziello, M

S. Capozziello, M. De Laurentis, S. D. Odintsov and A. St abile, Phys. Rev. D 83, 064004 (2011) [arXiv:1101.0219 [gr-qc]]

work page arXiv 2011
[12]

M. K. Mak and T. Harko, Proc. Roy. Soc. Lond. A 459, 393-408 (2003) [arXiv:gr-qc/0110103 [gr-qc]]

work page arXiv 2003
[13]

B. S. Guilfoyle, Gen. Rel. Grav. 31, 1645-1673 (1999) [arXiv:gr-qc/9906089 [gr-qc]]

work page arXiv 1999
[14]

B. Wang, E. Abdalla, F. Atrio-Barandela and D. Pavon, Re pt. Prog. Phys. 79, no.9, 096901 (2016) [arXiv:1603.08299 [astro-ph.CO]]

work page Pith review arXiv 2016
[15]

Indications of a late-time interaction in the dark sector

V. Salvatelli, N. Said, M. Bruni, A. Melchiorri and D. Wa nds, Phys. Rev. Lett. 113, no.18, 181301 (2014) [arXiv:1406.7297 [astro-ph.CO]]

work page Pith review arXiv 2014
[16]

Sol` a Peracaula, A

J. Sol` a Peracaula, A. G´ omez-Valent, J. de Cruz Perez a nd C. Moreno-Pulido, EPL 134, no.1, 19001 (2021) [arXiv:2102.12758 [astro-ph.CO]]

work page arXiv 2021
[17]

C. V. R. Board and J. D. Barrow, Phys. Rev. D 96, no.12, 123517 (2017) [erratum: Phys. Rev. D 98, no.12, 129902 (2018)] [arXiv:1709.09501 [gr-qc]]

work page arXiv 2017
[18]

C. G. Boehmer, [arXiv:gr-qc/0308057 [gr-qc]]

work page arXiv
[19]

C. G. Boehmer, Gen. Rel. Grav. 36, 1039-1054 (2004) [arXiv:gr-qc/0312027 [gr-qc]]

work page arXiv 2004
[20]

C. G. Boehmer, N. Dadhich and S. Das, [arXiv:2507.00503 [gr-qc]]

work page arXiv
[21]

Bruni, R

M. Bruni, R. Maier and D. Wands, Phys. Rev. D 105, no.6, 063532 (2022) [arXiv:2111.01765 [gr-qc]]

work page arXiv 2022
[22]

F. C. Alves and R. Maier, Class. Quant. Grav. 41, no.12, 125007 (2024) [arXiv:2405.05936 [gr-qc]]

work page arXiv 2024
[23]

Munyeshyaka, J

A. Munyeshyaka, J. Ntahompagaze, T. Mutabazi, M. R. Mbo nye, A. Ayirwanda, F. Twagirayezu and A. Abebe, Int. J. Geom. Meth. Mod. Phys. 20, no.03, 2350047 (2023) [arXiv:2210.01184 [gr-qc]]

work page arXiv 2023
[24]

Y. Yang, X. Dai and Y. Wang, Phys. Rev. D 111, no.10, 103534 (2025) [arXiv:2505.09879 [astro-ph.CO]]

work page arXiv 2025
[25]

Sola Peracaula, A

J. Sola Peracaula, A. Gomez-Valent, J. de Cruz Perez and C. Moreno-Pulido, Universe 9, no.6, 262 (2023) [arXiv:2304.11157 [astro-ph.CO]]

work page arXiv 2023
[26]

P. H. R. S. Moraes, J. D. V. Arba˜ nil and M. Malheiro, JCAP 06, 005 (2016) [arXiv:1511.06282 [gr-qc]]

work page arXiv 2016
[27]

P. O. Mazur and E. Mottola, Universe 9, no.2, 88 (2023) [arXiv:gr-qc/0109035 [gr-qc]]

work page arXiv 2023
[28]

Beltracchi and P

P. Beltracchi and P. Gondolo, Phys. Rev. D 99, no.4, 044037 (2019) [arXiv:1810.12400 [gr-qc]]

work page arXiv 2019
[29]

Chandrasekhar, Phys

S. Chandrasekhar, Phys. Rev. Lett. 12, 114-116 (1964)

1964