{"id":"726e4c91-b8cc-4b95-98bc-3b25ed7e8e94","arxiv_id":"1908.10316","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"low","formal_verification":"none","parameter_count":0,"one_line_summary":"A classical relativistic gas in a box has a double-spiral caloric curve that shrinks as the compactness parameter increases, vanishing above ν_max = 0.1764.","lead":"This paper computes the temperature-energy curves (caloric curves) of a hot, star-like gas held in a box and pulled together by gravity according to Einstein's theory. It finds that these curves form two linked spirals, and that there is a maximum compactness above which no equilibrium configuration can exist, so the gas always collapses.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Existence bound νmax=0.1764 assumes the high-Φ0 revived oscillations of Nα(Φ0) never exceed Nmax; this global-supremum claim is unproven and the scanned numerical range is finite.","rationale":"The paper is a careful numerical study, and the independent work of Roupas [167] provides support for the qualitative double-spiral picture and for a maximum compactness near 0.1764. My concern is therefore not that the authors fabricated or mis-plotted results; it is that the strong universal statement 'no equilibrium state exists for ν>νmax' exceeds what the presented numerics can currently certify. The method in Appendix C is a shooting method on a two-parameter family (α,Φ0); the full solution set requires a global maximum of Nα(Φ0). The paper's own footnotes identify the high-Φ0 revived oscillations and call them irrelevant because they are unstable. But the existence claim is not about stability: saddle points of entropy are still solutions of the equilibrium equations and are included on the plotted caloric curves. Unless those revived branches are shown to stay below Nmax for all α, the conclusion 'no equilibrium state is possible' is not established. The same gap underlies the reader's conditional verdict. Because Roupas's independent result reduces the likelihood of a large error, I would keep the conditional verdict rather than reject the paper; the concrete test I propose would settle whether the global-supremum gap is real or merely a missing proof.","tokens_in":50924,"tokens_out":6505,"duration_ms":70481,"concrete_test":"Implement the Appendix C shooting method independently, integrating the TOV equations (7)-(8) with the Maxwell-Juttner equation of state and R=1 at adaptive tolerance ≤1e-12. Scan α around α*=5.012 and more coarsely over α∈[-10^4,10^4]; for each α compute Nα(Φ0) for Φ0+1 spanning at least 10^-12 to 10^12, explicitly resolving the first several revived high-Φ0 oscillations. Test whether any local maximum anywhere exceeds Nmax=0.1764, and perform a step-size/refinement study to see whether the global maximum and α*=5.012 are converged to four digits. If a secondary peak exceeds Nmax, the no-equilibrium bound fails; if the supremum converges to 0.1764 at α*, the concern is resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Appendix C constructs equilibria as intersections of Nα(Φ0) with the line N. The central no-equilibrium claim (abstract and Sec. V.C) requires that, for every α∈(-∞,∞), the maximum of Nα(Φ0) over Φ0∈[-1,∞) is at most Nmax=0.1764 and is attained only at α*=5.012. The paper computes Nα(Φ0) on a finite range and, in footnote 51, dismisses the high-Φ0 revived oscillations as a mathematical curiosity whose solutions are unstable and 'irrelevant'. Instability does not remove existence: these are still extrema of the same TOV plus Tolman-Klein system, and the paper's own caloric curves include unstable branches. If any revived local maximum anywhere in the high-Φ0 region exceeded Nmax, then equilibrium states with N>Nmax would exist at high central densities, contradicting 'no equilibrium state is possible whatever the energy and temperature'. The paper also explicitly states that there is no mathematical proof for the merging of successive branches that defines the spiral topology, and it gives no convergence or error analysis for the numerical peak N(α)=0.1764. Thus the load-bearing point is not the spiral shape as such, but the unproved global supremum over the full Φ0 range, including the branches the paper chooses to ignore.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies box-confined classical self-gravitating gases in general relativity, modeled by the Maxwell-Juttner distribution with the TOV equations and the Tolman-Klein relations. It claims that the caloric curves η(Λ) depend on a single compactness parameter ν = GNm/Rc², and that these curves have a double-spiral structure: a cold spiral (weakly relativistic, generalizing the nonrelativistic curve) and a hot spiral (strongly relativistic, similar to the black-body radiation curve). As ν increases, the spirals approach, merge at ν'_S = 0.128, form a loop above ν_S = 0.1415, reduce to a point at ν_max = 0.1764, and then disappear. The central claim is that for ν > ν_max no equilibrium state exists at any energy or temperature. The paper describes the numerical construction of the caloric curves (Appendix C), the N→0 limits in both normalizations (Sec. VI), and the evolution of critical points (Sec. VII). It also discusses astrophysical consequences, including the distinction between slow thermodynamical and fast dynamical collapse.","tokens_in":51188,"tokens_out":3732,"duration_ms":41533,"significance":"If the result holds, the paper establishes a concrete compactness threshold beyond which a classical self-gravitating gas confined in a box has no isothermal statistical-equilibrium state at all, sharply illustrating how general relativity renders such systems more unstable than Newtonian gravity. The claimed double-spiral topology is a novel and physically suggestive generalization of the known Newtonian spiral and the black-body-radiation spiral. The paper deserves credit for presenting a clearly described numerical method (Appendix C), for openly flagging where results are numerical rather than proven (Appendix C.4, footnote 51, footnote 41), and for making falsifiable quantitative predictions (e.g., ν_max = 0.1764, Λ* = −0.9829, η* = 1.2203). The central quantitative claim is, however, not yet mathematically established; in particular, the global statement that no equilibrium exists for N > Nmax rests on a finite numerical scan and on dismissing high-central-potential branches rather than on a proof of the global supremum of Nα(Φ0).","major_comments":[{"comment":"The central no-equilibrium claim for N>Nmax requires proving that sup_{\\alpha,\\Phi_0} N_\\alpha(\\Phi_0) is attained at Nmax=0.1764 and never exceeded for any \\alpha \\in (-∞,∞) and \\Phi_0 \\in [-1,∞). The numerical construction in Appendix C, however, scans only a finite range of \\Phi_0, and footnote 51 explicitly dismisses the revived high-\\Phi_0 oscillations as a \"mathematical curiosity\" whose solutions are \"unstable\" and hence \"irrelevant\". Instability of those branches does not remove their existence: they are equilibrium solutions of the same TOV plus Tolman-Klein system, and the paper itself includes unstable branches on the caloric curves (e.g., Fig. 3). If any revived local maximum anywhere in the high-\\Phi_0 region exceeded Nmax, equilibrium states with N>Nmax would exist at high central densities, contradicting the statement in the abstract and Sec. V.E that no equilibrium state is possible whatever the energy and temperature. The manuscript explicitly concedes \"we do not have a mathematical proof\" for the merging of branches that defines the spiral topology, so the global supremum claim is a load-bearing assertion that is not presently supported.","section":"Abstract, Sec. V.E, and Appendix C.4"},{"comment":"The quantitative output of the paper — the critical particle number Nmax=0.1764, the critical point (Λ*, η*)=(−0.9829, 1.2203), the scaling laws in Eqs. (35), (38), (41), (44), (45)–(48), and the contrast values — is presented at an \"indicative level\" because \"the numerics is not very accurate close to Nmax\". No convergence test, tolerance, or error estimate is given for the numerical location of the peak of Nα(Φ0), nor for the determination of Λc, Λmin, ηc, ηmin by hand from the caloric curves (Appendix C.4). Since these numbers are the main quantitative results and since the entire existence bound rests on them, a convergence or error analysis is required before the claim can be considered established.","section":"Sec. VII and footnote 41"},{"comment":"The hot-spiral limit values Mmax = 0.24632 and Bmin = 17.809 are cited from the authors' own preprints (refs. [1], [3], [170]) rather than derived in this paper; the same applies to the asymptotic curve in Fig. 19. The comparison with the self-gravitating black-body radiation in Sec. III.B, where Mmax = 0.24632 is quoted from Refs. [162, 163], also relies on those external results. For a paper whose abstract emphasizes the hot spiral and its critical parameters, a self-contained derivation of these limit values (or at least a concise indication of the equations used to obtain them) should be included or made reproducible, otherwise the reader cannot verify a key input to the central scenario.","section":"Sec. VI.B and Sec. III.B"},{"comment":"The assertion that the high-\\Phi_0 revived oscillations are unstable and can therefore be ignored is itself not demonstrated. The paper presents no stability analysis (eigenvalue equation or second-variation calculation) for those branches; the instability claim in footnote 51 appears to be an inference from the fact that the corresponding solutions lie deep in the spirals of the caloric curve, but no explicit check is reported. Since these branches are the ones whose existence could overturn the N>Nmax conclusion, the instability assertion is load-bearing and should be supported by a computation or by an argument that the extrema of the same variational problem are captured by the lower-\\Phi_0 branches.","section":"Appendix C.1 and footnote 51"}],"minor_comments":[{"comment":"The word \"amputed\" appears twice (Sec. V.B and the conclusion); it should be \"amputated\".","section":"Sec. V.B"},{"comment":"The abstract states that the caloric curves \"typically\" have the form of a double spiral, but for NS<N<Nmax the curve is described as a single loop resembling the symbol ∞ (Fig. 13). The phrasing is not contradictory, but the abstract could be clearer that the double spiral is the typical shape only for N<N'S.","section":"Abstract and Sec. IV.A"},{"comment":"The dimensionless relation ν=N with ℏ=c=G=m=g/2=R=1 is stated without comment. Since the figures use N as the control parameter, an explicit sentence reminding the reader that ν=N in these units would improve accessibility.","section":"Eq. (19)"},{"comment":"Several references are listed only as \"preprint\" or \"arXiv\" with no arXiv number or publication status (e.g., refs. [1], [3], [170]). This makes it difficult for the reader to verify the cited values Mmax = 0.24632 and Bmin = 17.809 and the asymptotic caloric curves; the authors should provide arXiv identifiers or complete publication data.","section":"References [1], [3], [170]"},{"comment":"The dashed line in Fig. 16 is mentioned in the caption and text as corresponding to ηc = 2.52, but it is not labeled in the figure itself; adding a small label would aid the reader.","section":"Fig. 16"},{"comment":"The replacement of b0 by Φ0 via Eq. (16) introduces a division by |α|; the paper correctly notes that α=0 gives a singular presentation of Nα(Φ0) as a function of Φ0. This is a useful clarification, but a statement that the physical quantities remain regular at α=0 would be helpful, as the singularity is only in the plotting variable.","section":"Sec. II and Appendix C.2"}],"recommendation":"major_revision","confidential_remarks":"The manuscript relies heavily on the authors' own unpublished preprints (refs. [1], [3], [170]) for key limiting values and asymptotic curves. This is not itself a disqualifying practice, but the editors may wish to ensure that those preprints are publicly available and that the authors provide the missing derivations or a more detailed pointer in the revision. The paper is long and partly review-like, but it is within the scope of the journal; the main concern is the unproved global supremum claim behind ν_max=0.1764, on which the abstract's strongest statement rests. That concern is explicitly acknowledged by the authors in Appendix C.4 and footnote 51, so the manuscript is honest about its limitations; the revision should either supply a proof, a convergence study, or a clearly qualified reformulation of the no-equilibrium claim."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Two things to know. This paper is a confirmation-and-extension of Roupas [167], not a new discovery, and it says so plainly. The double-spiral caloric curve and the maximum compactness nu_max = 0.1764 were already in Roupas. What is actually new: a different normalization (M, B) that makes the hot spiral an asymptotic curve when N->0, explicit scaling laws for the critical points near Nmax, and a branch-by-branch description of how the two spirals merge, loop, and disappear. That material is useful.\n\nThe second thing: the paper's headline claim, 'no equilibrium state exists whatever the energy and temperature for N > Nmax,' is not proven. The stress-test concern is right. Appendix C builds equilibria as intersections of N_alpha(Phi_0) with a horizontal line, and the no-equilibrium claim needs N_alpha(Phi_0) to never exceed Nmax = 0.1764 over the whole Phi_0 range. The paper scans a finite range and then dismisses the revived high-Phi_0 oscillations as irrelevant because those solutions are unstable. But instability does not remove existence: they are still extrema of the same TOV plus Tolman-Klein system. The paper's own caloric curves include unstable branches. The authors explicitly say they have no mathematical proof for the merging of branches that defines the spiral topology, and the numerical peak values come with no convergence or error analysis. These are load-bearing caveats, not minor footnotes, because they sit exactly under the central existence bound.\n\nWhat the paper does well: the method is transparent, the asymptotic cold and hot limits are physically informative, and the comparison with Roupas is honest. The paper is also unusually candid about its limits, which I respect. The weak spots are the dependence on unpublished preprints for key numbers (Mmax = 0.24632 and Bmin = 17.809) and the absence of code or data, which makes the numerics uncheckable from the paper alone. Self-citation is not a sin, but here it does real work: the hot-spiral asymptotics rest on material the reader cannot see.\n\nVerdict: the double spiral and nu_max are very plausibly correct because Roupas independently got them; the issue is proof burden, not plausibility. The paper deserves peer review. I would send it out, with a request that the authors either prove the global-supremum claim or state it as a numerical conjecture, integrate the preprint dependencies, and provide code or data. Readers working on relativistic self-gravitating thermodynamics will find this a useful reference; readers wanting a rigorous existence theorem should look elsewhere.","headline":"Solid numerical extension of Roupas with an honest but unproven global existence bound; worth refereeing, not a new discovery.","tokens_in":51728,"tokens_out":3955,"would_cite":true,"duration_ms":42255,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":["04.40.Dg","05.70.-a","05.70.Fh","95.30.Sf","95.35.+d"],"model":"deepseek-v4-flash","headline":"A self-gravitating gas in a box has no equilibrium state in general relativity once its compactness parameter exceeds 0.1764; below that threshold the temperature-energy curve is a double spiral that shrinks with compactness.","keywords":["caloric curves","self-gravitating gas","general relativity","statistical equilibrium","compactness parameter","double spiral","gravothermal catastrophe","global temperature"],"falsifier":"An independent high-precision integration of the same equilibrium equations that finds any equilibrium state for $\\nu>0.1764$, at any energy and temperature, would refute the claimed maximum compactness; conversely, a relativistic Vlasov-Einstein or N-body simulation starting from a classical gas in a box with $\\nu>0.1764$ that does not collapse would challenge the conclusion.","tokens_in":50680,"feed_emoji":"🕳️","tokens_out":11559,"duration_ms":102287,"temperature":0.7,"pith_summary":"The paper sets out the complete equilibrium diagram of a classical self-gravitating gas confined in a spherical box and treated with general relativity. It argues that the thermodynamics is controlled by a single number, the compactness parameter $\\nu = GNm/(Rc^2)$, and that the caloric curve — inverse global temperature against binding energy — is a double spiral. One spiral extends the known nonrelativistic curve associated with the gravothermal catastrophe; the other is similar to the spiral of self-gravitating radiation but not identical to it. As $\\nu$ increases, the two spirals approach each other, merge at $\\nu'_S=0.128$, form a loop above $\\nu_S=0.1415$, and disappear at $\\nu_{\\max}=0.1764$, so above that compactness no equilibrium exists for any energy or temperature. The consequence is concrete: general relativity makes a self-gravitating gas more unstable than nonrelativistic gravity, with a maximum compactness beyond which collapse is unavoidable.","feed_headline":"General relativity caps how compact a gas can be at 0.1764","feed_subtitle":"Above that compactness, a boxed star cluster has no equilibrium at any energy or temperature, so collapse is inevitable.","key_machinery":"The argument runs on the series of equilibria of the general-relativistic hydrostatic equilibrium equations, closed with the relativistic equilibrium distribution function of a classical gas and parametrised by the uniform ratio $\\alpha$ of chemical potential to temperature and by the central gravitational potential $\\Phi_0$. Equilibrium states for a fixed particle number $N$ are exactly the intersections of the curve $N_\\alpha(\\Phi_0)$ with the horizontal level $N$; varying $\\alpha$ sweeps out branches of the caloric curve, and the turning points of these branches locate the stability changes through the classical turning-point criterion for linear series of equilibria. Two normalisations do the conceptual work: $(\\Lambda,\\eta)$ is adapted to the nonrelativistic limit and exposes the cold spiral, while $(M,B)$ is adapted to the ultrarelativistic limit and exposes the hot spiral as an asymptotic curve when $\\nu\\to0$. The maxima and minima of the curve $N_\\alpha(\\Phi_0)$ directly explain the critical particle numbers $N'_S=0.128$, $N_S=0.1415$ and $N_{\\max}=0.1764$ that organise the changes of topology.","core_discovery":"The central claim is that the caloric curve $\\eta(\\Lambda)$, with $\\eta$ the dimensionless inverse global temperature and $\\Lambda$ the dimensionless binding energy, captures the stability of the gas in both the microcanonical and canonical ensembles, and that for small compactness this curve is a double spiral. The cold spiral is the relativistic generalisation of the classical nonrelativistic spiral; the hot spiral, made visible by an ultrarelativistic normalisation, resembles the spiral of self-gravitating radiation with the same maximum mass-energy $GM_{\\max}/(Rc^2)=0.24632$ in the $\\nu\\to0$ limit. The claimed topology is a one-parameter family: well-separated double spirals for $\\nu<0.128$, touching truncated spirals for $0.128<\\nu<0.1415$, a single loop for $0.1415<\\nu<0.1764$, a single equilibrium point at $\\nu_{\\max}=0.1764$ with $\\Lambda_*=-0.9829$ and $\\eta_*=1.2203$, and no equilibrium at all above. If this is right, the gas has a maximum compactness $\\nu_{\\max}=0.1764$, equivalently a minimum box radius $R_{\\min}=5.67\\,GNm/c^2$, and increasing compactness advances both the low-energy and high-energy collapse thresholds.","pith_inferences":["The square-root scalings $\\Lambda_X-\\Lambda_*\\sim\\pm 6.7\\,(N_{\\max}-N)^{1/2}$ and $\\eta_X-\\eta_*\\sim\\pm 3.9\\,(N_{\\max}-N)^{1/2}$ suggest that the disappearance of equilibria at $\\nu_{\\max}$ is a fold (saddle-node) bifurcation of the equilibrium family; if so, the critical point should be reproducible by a normal-form expansion without resolving every spiral turn.","The 'irrelevant' high-density branches are set aside because they are unstable; if quantum degeneracy or a short-distance repulsion stabilised them, the caloric-curve topology and the maximum compactness would change, so the bound $\\nu_{\\max}=0.1764$ is specific to a purely classical gas.","The same maximum mass-energy $GM_{\\max}/(Rc^2)=0.24632$ appearing for both the classical hot spiral and self-gravitating radiation hints at a universal ultrarelativistic limit; a natural extension is to check whether the hot spiral of a self-gravitating Fermi gas approaches the same value as degeneracy is lowered.","A direct numerical-relativity test is available: a box-confined classical gas started with $\\nu>0.1764$ should collapse at every scanned energy and temperature, and finding even one long-lived equilibrium would contradict the existence bound."],"forward_implications":["For fixed box radius $R$, equilibrium exists only up to the particle number $N_{\\max}=0.1764\\,Rc^2/(Gm)$; above it the gas has no equilibrium at any energy or temperature and is expected to collapse.","For fixed particle number $N$, there is a minimum box radius $R_{\\min}=5.67\\,GNm/c^2$ below which no equilibrium state exists.","Raising the compactness parameter lowers the critical energy and critical temperature of the cold spiral relative to the nonrelativistic values $\\Lambda_c=0.335$ and $\\eta_c=2.52$, so general relativity makes the gas unstable earlier than nonrelativistic gravity does.","The hot spiral endows high positive energies with a new instability: above a maximum energy and above a maximum global temperature the gas collapses, in analogy with the self-gravitating radiation case.","The two collapse routes differ in timescale: the low-energy (cold-spiral) instability is slow and secular, while the high-energy (hot-spiral) instability is fast and dynamical, so the two ends of the caloric curve correspond to different physical collapse mechanisms."],"supporting_citations":[{"why":"Supplies the microcanonical instability threshold (density contrast 709) that the cold spiral must reproduce in the nonrelativistic limit.","marker":"[57]"},{"why":"Defines the nonrelativistic caloric curve and gravothermal catastrophe that the cold spiral generalises.","marker":"[58]"},{"why":"Provides the turning-point criterion used to locate stability changes along the series of equilibria.","marker":"[64]"},{"why":"Establishes the spiralling shape of the nonrelativistic caloric curve and the turning-point analysis the paper extends.","marker":"[76]"},{"why":"Gives the self-gravitating radiation caloric curve whose spiral and maximum mass-energy anchor the hot-spiral comparison.","marker":"[163]"},{"why":"Provides the independent derivation of general-relativistic classical-gas caloric curves that this paper confirms and completes.","marker":"[167]"},{"why":"Supplies the numerical method of intersecting the particle-number curve with the particle-number level, used to build the caloric curves.","marker":"[168]"}],"fun_headline_variants":["Relativity sets a hard limit on star cluster compactness","No equilibrium for gravitating gas beyond ν=0.1764","Double spiral caloric curve collapses as compactness grows","When boxes of stars get too dense, equilibrium vanishes","Gravitational collapse inevitable past compactness limit 0.1764"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is that the numerical construction of the series of equilibria is complete: every equilibrium state corresponds to an intersection of the curve $N_\\alpha(\\Phi_0)$ with the level $N$, and the revived high-density branches at very large central potentials are unstable and can be discarded; the paper states that it has no mathematical proof that successive branches merge exactly at the turning points that define the spiral structure.","fun_headline_variants_meta":{"raw":{"variants":["Relativity sets a hard limit on star cluster compactness","No equilibrium for gravitating gas beyond ν=0.1764","Double spiral caloric curve collapses as compactness grows","When boxes of stars get too dense, equilibrium vanishes","Gravitational collapse inevitable past compactness limit 0.1764"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000313,"raw_usage":{"total_tokens":1870,"prompt_tokens":1125,"completion_tokens":745,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":741,"completion_tokens_details":{"reasoning_tokens":661}},"tokens_in":741,"tokens_out":745,"duration_ms":6878,"temperature":1.0,"reasoning_tokens":661,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-14T10:47:01.591923+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"An independent high-precision integration of the same equilibrium equations that finds any equilibrium state for $\\nu>0.1764$, at any energy and temperature, would refute the claimed maximum compactness; conversely, a relativistic Vlasov-Einstein or N-body simulation starting from a classical gas in a box with $\\nu>0.1764$ that does not collapse would challenge the conclusion.","supporting_citations":[],"review_version":1}