{"id":"224d17d1-ccd1-4054-a022-560adc0175ba","arxiv_id":"2607.19815","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"Vlasov gas accretion onto a Schwarzschild black hole is shown explicitly to have coordinate-invariant density, pressures, and accretion rates, with captured particles carrying mean energy m0+kBT and lower specific entropy.","lead":"Physicists derive the accretion of a collisionless gas onto a Schwarzschild black hole in a coordinate-independent form, showing that density, pressure, and accretion rates do not depend on the chosen metric coordinates. The paper also finds that captured particles are on average less energetic than the surrounding gas and carry lower specific entropy.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Entropy-flux definition in Eq. (139) is not the standard Boltzmann current: it uses covariant p_r and omits the +f term. Recomputing with S^r = -k_B∫p^r f(ln f -1)dvol changes the abstract's 3/2 k_B entropy reduction to 1/2 k_B; coordinate-independence claim is unaffected.","rationale":"Reading the paper in good faith, the central claim is coordinate independence of the Vlasov accretion observables. That claim is well supported: n is a scalar by construction, the eigenvalues of T^μ_ν are invariants, and the accretion rates are derived from conserved currents; the explicit h-covariant expansions show the cancellations. I do not see a threat to that core result. The genuinely load-bearing weakness is the thermodynamic claim in the abstract and conclusion: the specific-entropy reduction of 3/2 k_B. The entropy flux defined in Eq. (139) deviates from the standard Boltzmann entropy current in two ways—using p_r instead of p^r and omitting the +f term—while the global entropy quoted in the same section uses the standard convention. This is not just a cosmetic factor: the difference between the two conventions changes the claimed reduction from 3/2 k_B to 1/2 k_B. The reader's weakest assumption identified exactly this issue; my independent check agrees. Since the paper's core coordinate-independence claim survives, the verdict should remain CONDITIONAL: the entropy section must be corrected, but the main theorem does not fall. No ad hominem is involved; this is a technical convention error that is fully checkable by recomputation.","tokens_in":20264,"tokens_out":29239,"duration_ms":262560,"concrete_test":"Recompute the horizon entropy flux with the standard current S^r = -k_B∫_{abs} p^r f(ln f - 1)dvol, using the absorption-branch momentum p^r = -√R0/r^2 and dvol = m0 l/√R0 dE dl dσ. Using \\bar E = ⟨E l_c^2⟩/⟨l_c^2⟩ = m0(1+1/z) from Eq. (138), the absorbed specific entropy becomes k_B(1 - ln A + z + 1) = k_B(-ln A + z + 2). Compare with the global σ∞ = k_B(1 - ln A + z + 3/2) = k_B(-ln A + z + 5/2); the difference is 1/2 k_B, not 3/2 k_B. As a second check, evaluate Eq. (139) in the h=0 and h=1 charts; if it changes, the p_r-based flux is coordinate-dependent and cannot define the physical entropy accretion rate.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing concern is not the coordinate-independence claim, which is sound: n, the mixed-tensor eigenvalues, and the Killing-flux accretion rates are scalars or conserved quantities whose invariance is checked by the explicit expansions. The problem is the entropy result advertised in the abstract and Sec. XII. Eq. (139) defines the entropy flux as 4πr^2 k_B ∫_{abs} f ln f p_r dvol. This differs from the standard relativistic Boltzmann entropy current S^r = -k_B∫ p^r f(ln f - 1) dvol in two ways: it uses the covariant momentum p_r instead of the contravariant p^r, so the flux is not the contraction of S^μ with the radial normal and is not manifestly coordinate-invariant; and it omits the +f term. Equation (140) then substitutes -4π∫ f p_r dvol → dot n and -4π∫ E f p_r dvol → dot E. With this convention the absorbed specific entropy is σ_abs = -k_B ln A + \\bar E/T = k_B(-ln A + z + 1). The global specific entropy quoted in the same section, σ∞ = k_B(-ln A + z + 5/2), uses the standard convention with both p^μ and the +f term. Using the standard current for both the global and absorbed subensembles gives σ_abs = k_B(-ln A + z + 2), which is only 1/2 k_B below the global value, not 3/2 k_B. The advertised 3/2 k_B reduction therefore arises from mixing conventions; fixing the entropy flux invalidates an abstract-level thermodynamic claim, though it leaves the coordinate-independence theorem intact.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies stationary, spherically symmetric accretion of a collisionless Vlasov gas onto a Schwarzschild black hole, using the general stationary metric with off-diagonal g_tr = -h(r). It constructs action-angle variables, solves the Liouville/Vlasov equation, and computes the particle current and stress-energy tensor for Fermi-Dirac, Maxwell-Jüttner, and Bose-Einstein distributions. The central claim is that the particle number density n = sqrt(-g^{μν}J_μJ_ν), the eigenvalues of the mixed stress-energy tensor (energy density ρ, radial pressure p_rad, tangential pressure p_tan), and the accretion rates ˙n and ˙E are independent of the coordinate choice among Schwarzschild, Eddington-Finkelstein, and Painlevé-Gullstrand-type coordinates. Asymptotic expansions at infinity and near the horizon are presented, together with numerical results at finite radii. The paper also claims that in the classical limit captured particles have mean energy m_0 + k_B T and specific entropy lower than the global Maxwell-Jüttner average by 3/2 k_B.","tokens_in":20698,"tokens_out":10841,"duration_ms":105161,"significance":"If the coordinate-invariance result stands, it is a useful formal clarification of the Vlasov accretion formalism: the explicit h_0/H_0 cancellations in Secs. VIII and IX support the claim that the physical observables are well defined, and the low-temperature horizon results reproduce the known Rioseco-Sarbach values. The paper contains no parameter fitting, and the expansions are detailed enough to be checked. However, part of the invariance is by construction—n is a norm of a vector, the pressures are eigenvalues of a tensor, and the accretion rates are conserved Killing currents—so the formal invariance alone is not a new physical principle. The advertised entropy-reduction result is not supported as stated: the entropy flux in Eq. (139) is defined with a nonstandard f ln f current and a covariant momentum component, and using the standard Boltzmann entropy current changes the numerical claim from 3/2 k_B to 1/2 k_B. Since this result appears in the abstract and the concluding section, the paper requires revision before it can be accepted.","major_comments":[{"comment":"The entropy flux is defined as 4πr² k_B ∫_{abs} f ln f p_r dvol. This is not the standard relativistic Boltzmann entropy current S^μ = -k_B ∫ p^μ f(ln f - 1)dvol. Two differences matter: (i) the covariant momentum component p_r is used instead of the contravariant p^r entering S^r, so the flux is not a contraction of the entropy current with the radial normal and is not manifestly coordinate-invariant; (ii) the +f term in ln f - 1 is omitted. When the standard entropy current is used consistently for both the absorbed subensemble and the asymptotic thermal gas, the absorbed specific entropy is k_B(-ln A + z + 2), not k_B(-ln A + z + 1). The global value quoted in the same section, k_B(-ln A + z + 5/2), uses the standard convention. The claimed reduction is therefore 1/2 k_B, not 3/2 k_B. The abstract and Sec. XII's thermodynamic conclusion do not follow from the calculation as written.","section":"Sec. X, Eq. (139)"},{"comment":"The substitution -4π∫ f p_r dvol → ˙n and -4π∫ E f p_r dvol → ˙E is not an identity for the metric (1). The accretion rates are defined through the contravariant component J^r = ∫ p^r f dvol, and p^r = g^{rμ}p_μ = b² p_r + h E when a²b²+h² = 1. Thus the integral over p_r differs from the integral defining ˙n by metric-dependent terms involving ∫ E f dvol. Even aside from the missing +f term, Eq. (139) does not define the physical entropy accretion rate unless this relation is established or the flux is redefined with p^r. The authors should either prove the required cancellation or replace Eq. (139) by the standard S^r = -k_B ∫ p^r f(ln f - 1)dvol.","section":"Sec. X, Eq. (140)"}],"minor_comments":[{"comment":"The expression for n∞ is written with g_{tt}, g_{tr}, g_{rr} as if using covariant components with the metric itself; since the quantity is sqrt(-g^{μν}J_μJ_ν), the inverse metric should appear. Please correct the notation to avoid ambiguity.","section":"Eq. (84)"},{"comment":"The caption lists 'p^H_rad/n^H' twice; one instance is presumably meant to be a different quantity (e.g., p^H_tan/n^H). Please fix.","section":"Fig. 1 caption"},{"comment":"The identity a²b² + h² = 1 is stated for the Schwarzschild metric. Since it is used throughout as a restriction on the metric family, make explicit that it holds for the coordinate representations considered here and is not a property of a generic stationary spherically symmetric spacetime such as Reissner-Nordström.","section":"Sec. II, Eq. (3)"},{"comment":"The phrase 'most general stationary spherically symmetric spacetime' is stronger than what is analyzed: the paper specializes to Schwarzschild with the identity (3) and particular expansions for h(r). Please qualify the scope.","section":"Abstract and Sec. I"},{"comment":"The global specific entropy σ∞ = k_B(-ln A + z + 5/2) is quoted without derivation. Give the standard Maxwell-Jüttner entropy density formula so the comparison with Eq. (141) is transparent and the convention is clear.","section":"Sec. X, after Eq. (141)"},{"comment":"There are minor typographical issues, e.g., 'similer' before Eq. (122), and in Table I the entries for z=10 and z=30 are identical across ε values to many digits; this may be intended but should be commented on.","section":"General"}],"recommendation":"major_revision","confidential_remarks":"The coordinate-independence part of the paper appears sound and useful, but the entropy claim is a headline result stated in the abstract and conclusion, and it is not supported by the calculation because of a nonstandard and inconsistent entropy-flux definition. This is a load-bearing issue that must be fixed, so I recommend major revision rather than rejection. The authors should also be asked to compare their entropy current with standard references and to state clearly which convention they adopt for the absorbed subensemble."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"I read it through. The central claim holds up: the particle number density n, the eigenvalues of the mixed stress-energy tensor, and the conserved accretion rates are coordinate-independent, and the explicit h-coordinate expansions in Secs. VIII–IX are internally consistent. The h0/H0 cancellations work, and the low-temperature accretion rates reproduce Rioseco-Sarbach, which is a good sanity check. Credit where it's due: this is a careful, systematic calculation that puts the coordinate-invariance folklore on explicit footing, even if part of it is definitional (n is a norm, eigenvalues are eigenvalues). The diagonalization argument and the demonstration that -4πr^2 J^r and -4πr^2 T^r_t are the same constants in Schwarzschild, Eddington-Finkelstein, and Painlevé-Gullstrand-type coordinates are useful.\n\nThe real problem is the entropy claim in the abstract and Sec. X. Eq. (139) defines the entropy flux with p_r and f ln f, omitting the +f term that appears in the standard Boltzmann entropy current S^μ = -k_B ∫ p^μ f(ln f - 1) dvol. The global specific entropy quoted just below uses the standard convention. Comparing those two conventions is what produces the advertised 3/2 k_B reduction. If you use one convention consistently, the number changes—the stress-test says 1/2 k_B, and that sounds right. So the abstract-level thermodynamic result is not supported as stated. This does not touch the coordinate-invariance theorem, but it does mean the paper's headline news is partly wrong.\n\nTwo smaller issues. First, 'the most general stationary spherically symmetric spacetime' is overclaimed: the calculations are for Schwarzschild, and the identity a^2b^2+h^2=1 is used everywhere. The claim of ready extension to Reissner-Nordström or modified gravity is not demonstrated. Second, the mean-energy result m0+kBT is fine, but the paper should state explicitly whether it is new or follows immediately from the Rioseco-Sarbach rates—right now it reads as if it might be a trivial ratio of their ˙E and ˙n.\n\nI'd send this to peer review. The coordinate-invariance part is solid, the failure mode is localized, and the authors clearly know the kinetic-theory literature. A competent referee can force the entropy-flux fix and tone down the generalization claims. After that, it's publishable. I wouldn't cite the entropy claim until it's corrected, but I'd cite the h-coordinate treatment if I were working on kinetic accretion.","headline":"The coordinate-invariance theorem is real and the expansions check out, but the abstract's 3/2 k_B entropy reduction is an artifact of mixing entropy-flux conventions; the paper needs revision, not rejection.","tokens_in":734,"tokens_out":879,"would_cite":true,"duration_ms":48007,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"Coordinate choice does not change the density, pressures, or accretion rates of a collisionless gas falling into a Schwarzschild black hole.","keywords":["Schwarzschild black hole","Vlasov gas","accretion","coordinate invariance","particle number density","stress-energy tensor","angular momentum filtering","Maxwell-Jüttner distribution"],"falsifier":"Recompute the entropy accretion rate with the standard Boltzmann entropy current S^μ = −k_B ∫ p^μ f(ln f − 1) dvol; if the specific entropy of accreted particles then differs from the global average by 1/2 k_B instead of 3/2 k_B, the paper's entropy claim is refuted while its coordinate-independence claim is untouched.","tokens_in":20131,"feed_emoji":"🕳️","tokens_out":10788,"duration_ms":92306,"temperature":0.7,"pith_summary":"This paper shows that when a collisionless (Vlasov) gas accretes onto a Schwarzschild black hole in a stationary, spherically symmetric flow, the physically observable quantities—particle number density, energy density, radial and tangential pressures, and the particle and energy accretion rates—are independent of which coordinate system is used to describe the spacetime. Individual components of the particle current and stress-energy tensor do change with the coordinates (for example, in Schwarzschild versus Eddington-Finkelstein versus Painlevé-Gullstrand forms), but certain invariants, built from the metric and these tensors, do not. The paper also finds that in the classical limit the black hole captures particles with mean energy m₀ + k_BT, lower than the m₀ + 3/2 k_BT of the ambient gas, and with a specific entropy lower by 3/2 k_B, because high-energy particles are preferentially scattered away by the centrifugal barrier. If correct, the result makes kinetic accretion theory coordinate-free and sharpens our picture of how black holes grow thermodynamically.","feed_headline":"Coordinate choice cannot change black hole gas accretion","feed_subtitle":"In a collisionless Vlasov gas around a Schwarzschild black hole, density, pressures, and accretion rates stay identical in every coordinate","key_machinery":"The central object is the most general stationary spherically symmetric metric ds² = −b²dt² − 2h dtdr + a²dr² + r²(dθ² + sin²θ dφ²), together with the identity a²b² + h² = 1 that holds for the Schwarzschild spacetime in all coordinate representations considered (orthogonal, Eddington-Finkelstein, and Painlevé-Gullstrand type). This identity makes the volume element √−g = r² sinθ and reduces the conservation equations to ∂_r(r²J^r) = 0 and ∂_r(r²T^r_t) = 0, so the accretion rates take the simple coordinate-invariant forms −4πr²J^r and −4πr²T^r_t. Action-angle variables are constructed to solve the Vlasov equation, and the absorption/scattering split is controlled by the critical angular momen","core_discovery":"For stationary, spherically symmetric Vlasov gas accretion onto a Schwarzschild black hole, the particle number density n = √(−g_{μν}J^μJ^ν), the eigenvalues of the mixed stress-energy tensor T^μ_ν (energy density ρ, radial pressure p_rad, tangential pressure p_tan), and the accretion rates ṅ and Ė are invariant under coordinate changes, even though the components J^μ and T^{μν} are coordinate-dependent. The simple flux formulas −4πr²J^r and −4πr²T^r_t hold in every coordinate system because of the identity a²b² + h² = 1 satisfied by the Schwarzschild metric in the general stationary spherically symmetric form. In the classical (Maxwell–Jüttner) limit, the mean energy of captured particles i","pith_inferences":["Using the standard Boltzmann entropy current (with the −1 term) instead of the truncated flux in Eq. (139) would shift the specific-entropy reduction to 1/2 k_B; the qualitative conclusion that accreted particles carry lower entropy would survive, but the quantitative 3/2 k_B figure would change.","The same angular-momentum filtering should operate in other spherically symmetric spacetimes: whenever an effective potential has a centrifugal barrier, captured particles will be biased toward lower energies, so the mean-energy deficit m₀+k_BT versus m₀+3/2 k_BT may be a generic feature of collisionless capture.","The coordinate-invariant scalars n, ρ, p_rad, p_tan provide natural diagnostic quantities for numerical simulations of kinetic accretion, which often use horizon-penetrating coordinates; comparing these scalars across codes with different coordinate choices would directly test implementation consistency."],"forward_implications":["The accretion theory becomes fully geometric: any coordinate system in the ansatz yields the same density, pressures, and accretion rates, so future work can choose coordinates for convenience (e.g., horizon-penetrating for numerics) without worrying about physical predictions.","The particle and energy accretion rates take the simple, universal forms −4πr²J^r and −4πr²T^r_t at every radius, not just at the horizon, because the metric identity reduces the flux conservation equations to radial constancy.","The angular-momentum filtering means the black hole preferentially absorbs lower-energy, lower-entropy particles; in the classical limit their mean energy is m₀+k_BT instead of m₀+3/2 k_BT, which changes how the black hole's mass and entropy budgets are computed in kinetic models.","At low temperature the three statistics (Fermi-Dirac, Maxwell-Jüttner, Bose-Einstein) converge, so the coordinate-independence and energy-filtering results are robust across quantum statistics."],"fun_headline_variants":["Black hole gas accretion: coordinates don't matter","Vlasov gas accretion rates are coordinate-independent","Schwarzschild accretion: physics same in any frame","Coordinate choice can't change accretion physics","Accretion observables immune to coordinate changes"],"cache_read_input_tokens":2304,"weakest_assumption_plain":"The 3/2 k_B entropy reduction depends on an entropy-flux definition (Eq. 139) that omits the standard '+f' term; with the conventional Boltzmann entropy current used for the global average, the reduction changes to 1/2 k_B, though the coordinate-independence claims are unaffected.","fun_headline_variants_meta":{"raw":{"variants":["Black hole gas accretion: coordinates don't matter","Vlasov gas accretion rates are coordinate-independent","Schwarzschild accretion: physics same in any frame","Coordinate choice can't change accretion physics","Accretion observables immune to coordinate changes"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.0005,"raw_usage":{"total_tokens":2300,"prompt_tokens":777,"completion_tokens":1523,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":521,"completion_tokens_details":{"reasoning_tokens":1467}},"tokens_in":521,"tokens_out":1523,"duration_ms":11797,"temperature":1.0,"reasoning_tokens":1467,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-01T11:38:54.380904+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Recompute the entropy accretion rate with the standard Boltzmann entropy current S^μ = −k_B ∫ p^μ f(ln f − 1) dvol; if the specific entropy of accreted particles then differs from the global average by 1/2 k_B instead of 3/2 k_B, the paper's entropy claim is refuted while its coordinate-independence claim is untouched.","supporting_citations":[],"review_version":1}