{"id":"07318155-15ce-4cc9-b92f-bffe026696b7","arxiv_id":"2507.16444","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"high","formal_verification":"none","parameter_count":2,"one_line_summary":"A Hamiltonian formulation of relativistic fluid dynamics in curved spacetime in coordinate time is proposed and applied to Schwarzschild radial flows, but the derived equations and stability conclusion are not fully supported.","lead":"This paper proposes a Hamiltonian framework for relativistic fluids in curved spacetime that keeps Newtonian time and familiar density and velocity variables. It applies the framework to radial flows near a Schwarzschild black hole and claims those flows are stable against small perturbations.","discovery_kind":"new_method","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The Eulerian density mapping in Eq. (6) omits metric volume and γ factors, and the continuity equation it generates already disagrees with dust dynamics in flat space; the central 'valid in generic spacetimes' claim is unverified.","rationale":"I read the paper as attempting to construct a Hamiltonian fluid dynamics in coordinate time whose equations reproduce relativistic fluid behavior in curved spacetime. The load-bearing foundation is the Eulerian mapping of Eq. (6), so the correctness of that mapping determines whether the derived continuity and Euler equations are meaningful. The paper never validates the mapping against covariant conservation or a known solution, and the flat-space check exposes an internal discrepancy: the same construction yields different flux factors depending on which displayed equation is used, and none matches the standard dust continuity. This is an internal inconsistency, not merely a disagreement with an alternative formulation. The stability analysis is also problematic because the rates κ and λ are introduced by ansatz rather than determined by boundary conditions, but the mapping issue is more fundamental because it undermines the basic equations before stability is considered. For these reasons I agree with the reader's REJECT verdict; no verdict adjustment is needed.","tokens_in":12556,"tokens_out":16611,"duration_ms":192991,"concrete_test":"Take g_μν=η_μν and recompute ∂_tρ={ρ,H} from the stated bracket (18) and Hamiltonian (20) in flat space. If the resulting continuity equation contains flux ρU^i/√(1+U^2), while the standard dust continuity requires ρU^i (for ρ=nγ), then the Eulerian mapping is internally inconsistent and the central claim fails.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (6) defines ρ(r,t)=ρ0∫d^3x δ(X(x,t)−r) and j^i=ρ0∫d^3x U^i δ(...), and the same mapping is used unchanged in curved spacetime. A density defined this way is not the time component of a conserved current: covariant conservation is ∂_μ(√{−g} n U^μ)=0, and neither √{−g} nor a γ factor appears in the definition. The resulting inconsistency is already visible in flat spacetime. With the flat limit of Eq. (20), H=∫ρ√(1+U^2), and the flat limit of the bracket (18), Hamilton's equation gives ∂_tρ = −∂_i(ρ U^i/√(1+U^2)). This differs from Eq. (24), which has an extra factor √(1+U^2) multiplying the flux, and from Eq. (26), which divides by that factor; all three differ from the standard dust continuity ∂_t(nγ)+∂_i(nγ v^i)=0, or ∂_tρ+∂_i(ρ v^i)=0 if ρ=nγ. The paper never identifies which quantity ρ represents and never checks against ∂_μ(√{−g} n U^μ)=0 or Bondi accretion, so the claim that Eqs. (24)-(27) are the relativistic fluid equations is not supported.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes an Eulerian Hamiltonian formulation of relativistic fluid dynamics in curved spacetime, in which the fluid is described by a density ρ and a coordinate three-velocity U^i evolving in coordinate (Eulerian) time. The authors posit a relativistic Poisson algebra, Eqs. (16)-(18), and a Hamiltonian, Eq. (20), and claim that these yield generalized continuity and Euler equations, Eqs. (21)-(23), valid for generic stationary metrics. They then specialize to the Schwarzschild metric, derive stationary radial background flows, Eqs. (35)-(37), and study first-order perturbations, concluding that the background flow is linearly stable and that vorticity is generated at linear order. The paper presents explicit plots of the background profiles and perturbation eigenfunctions and argues that the coordinate-time description is advantageous for astrophysical applications near horizons.","tokens_in":12929,"tokens_out":5755,"duration_ms":64759,"significance":"If correct, the framework would provide a Hamiltonian route to relativistic fluid equations in strong-field gravity using only coordinate-time variables, potentially useful for accretion modeling and numerical relativity. The paper does contain a concrete construction: an explicit Poisson algebra, an explicit Hamiltonian, an exact stationary radial background solution in Schwarzschild spacetime, and a perturbation analysis with plotted solutions. These are testable and clearly presented. However, the central claim of validity in generic curved spacetimes is not supported: the density mapping is carried over from the non-relativistic theory without metric volume or Lorentz factors, the derived continuity equation disagrees with standard special-relativistic conservation already in flat space, and the derivation of the central equations is not shown. These are load-bearing issues that affect every subsequent result, including the background solution and the stability conclusion. As it stands, the paper cannot be accepted as a reliable derivation of relativistic fluid dynamics.","major_comments":[{"comment":"The Eulerian mapping ρ(r,t)=ρ0∫d^3x δ(X(x,t)-r) and j^i=ρ0∫d^3x U^i δ(...) is used unchanged in curved spacetime, without a √{-g} volume element or a Lorentz factor. This cannot represent the time component of a conserved particle current, since covariant conservation requires ∂_μ(√{-g} n U^μ)=0. The inconsistency is already visible in the flat-space limit: from the bracket (18) and Hamiltonian (20), Hamilton's equation gives ∂_t ρ = -∂_i(ρ U^i/√(1+U^2)), whereas Eq. (24) states a different flux structure. The paper never defines which physical quantity ρ denotes or verifies the equations against standard special-relativistic continuity, so the claim that Eqs. (24)-(27) are the relativistic fluid equations is unsupported.","section":"§2, Eq. (6), and §3, Eqs. (21)-(24)"},{"comment":"The central equations of motion are introduced after the phrase \"after a long and tedious computation\" with no derivation shown. These equations are the foundation for all later results, but the paper provides neither the intermediate steps nor an independent consistency check, such as verifying Jacobi identities for the proposed algebra or demonstrating that the equations reduce to a known relativistic limit beyond the non-relativistic case. Without a verifiable derivation, Eqs. (21)-(23) cannot be considered established.","section":"§3, Eqs. (21)-(23)"},{"comment":"The steady-state continuity equation for the spherically symmetric background, Eq. (33), contains a factor r² inside the radial derivative: ∂_r(r² f ρ U0)=0. However, Eq. (36) drops this factor and asserts ∂_r(f ρ0U0)=0, leading to ρ0=χ/(fU0). The correct stationary solution of the paper's own continuity equation would satisfy r² f ρ0 U0 = constant, not fρ0U0=χ. Consequently, the background density profiles in Figure 3 and all subsequent perturbation results are not solutions of the equations the paper itself writes down.","section":"§3, Eq. (36)"},{"comment":"The perturbations are written as u(r,t)=h(r)e^{-κt} and ω(r,t)=g(r)e^{-λt}, with κ and λ appearing as free constants. The paper then concludes from the exponential time dependence that the perturbations decay and the background is stable. This is circular: a positive κ would give growth, and the analysis never determines the admissible signs or ranges of κ and λ from the dynamical equations. A genuine stability statement requires a dispersion relation or an eigenvalue problem that fixes these constants; the present treatment assumes the decay it claims to establish.","section":"§4, Eqs. (46)-(48) and stability claim (I)"}],"minor_comments":[{"comment":"There are mismatched braces in the displayed brackets, e.g., {j^i(r), ρ(r′} and {U^i(r), ρ(r′}, and the text refers to \"the PB (41)\" although no Eq. (41) exists; the reference should be to Eq. (3) or to the corresponding equation in the text.","section":"§2, Eqs. (16) and (18)"},{"comment":"The phrase \"spaceL time broken-up scenario\" appears to be a typographical error for \"spacetime broken-up scenario.\"","section":"Introduction"},{"comment":"The symbol ω is used both for the generalized vorticity components in Eq. (45) and for the angular perturbation satisfying Eq. (44), which makes the stability subsection unnecessarily confusing.","section":"§4"},{"comment":"The notation √Ψ is used without a clear definition of Ψ; it appears to be defined implicitly by the second line of Eq. (20), but the same symbol is also used in Eq. (19) with a different meaning.","section":"§2, Eq. (20)"},{"comment":"Reference [15] contains stray bracketed items, including \"[12] C.Eckart\" and \"[13]\", which appear to be remnants of an earlier draft; these should be cleaned up.","section":"References"},{"comment":"Figure captions are incomplete, and the figures are not described with enough quantitative detail; for example, Figure 3's caption ends with \"for\" and gives no axis labels or parameter values.","section":"Figures"}],"recommendation":"reject","confidential_remarks":"The paper has a clear and ambitious goal, but the central construction is inconsistent with standard relativistic fluid conservation laws even in flat space, and the derivation of the main equations is not provided. The stability claim is also circular because the decay rates are inserted by hand. These problems are not local presentation issues; they affect the paper's main claims and would require a reformulation of the mapping and a re-derivation of the equations, which is beyond a normal revision. I therefore recommend rejection, despite the paper's useful concrete examples and plots."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Dear [Colleague],\n\nThe short version: the paper sets out to build a coordinate-time Hamiltonian formulation of relativistic fluids in curved spacetime, and the central claim—that it produces the relativistic continuity and Euler equations in generic spacetimes—does not survive contact with flat spacetime. That is the thing to know before reading further.\n\nWhat is genuinely new: the extension of the Jackiw–Nair–Pi–Polychronakos Eulerian bracket construction to curved spacetime, with density and coordinate velocity as the fundamental fields, is a real idea. The Poisson algebra in (16)–(18) and the Hamiltonian in (20) are written down in a way that reduces to the nonrelativistic structure, and the Schwarzschild application produces the expected horizon freezing of coordinate velocity. The generalized vorticity and its appearance at linear order is a nice twist.\n\nThe soft spots are not small. The density mapping (6) uses the flat-space definition ρ(r,t)=ρ0∫d³x δ(X(x,t)−r), carried over unchanged to curved spacetime, with no √−g and no γ factor. That means ρ is not the time component of a conserved current. In the flat limit, Eq. (24) becomes ∂tρ + ∇·(ρ v/√(1+v²)) = 0, which is not the special-relativistic dust continuity equation ∂t(nγ)+∇·(nγv)=0 (or ∂tρ+∇·(ρv)=0 if ρ=nγ). The paper never says which quantity ρ is, and never checks against covariant conservation or Bondi accretion. That is a load-bearing flaw.\n\nThe derivation of (21)–(23) is hidden behind 'a long and tedious computation,' which is not acceptable for equations that are then advertised as generic. The stability analysis assumes perturbations of the form e^{−κt} and e^{−λt} with arbitrary κ and λ, solves for the spatial profile, and then reports decay as a result. That is circular; no mode is ever forced to decay. Finally, the stationary background solutions are obtained from the O(U²)-truncated equations, so calling them 'exact' in the abstract is misleading.\n\nOn balance, this is not a crank paper—the construction is coherent in its own terms—but the physical mapping is wrong at a fundamental level. A would-be reader should treat the Hamiltonian algebra as a formal exercise, not as a derivation of relativistic fluid equations. My recommendation: desk reject in its current form. If the authors can fix the density mapping, verify the flat-space limit, and show a genuine spectral analysis, the paper might become salvageable. As it stands, the central claims are not supported.","headline":"The flat-space limit already breaks the central claim; the Hamiltonian idea is worth a look, but the mapping is wrong and the stability analysis is circular.","tokens_in":13366,"tokens_out":7611,"would_cite":false,"duration_ms":71930,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":false},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C55","76Y05","83C57"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proposes a relativistic Hamiltonian formulation of fluid dynamics in curved spacetime, with dynamics cast in coordinate time, yielding generalized continuity and Euler equations and, in Schwarzschild spacetime, exactly solvable…","keywords":["relativistic fluid dynamics","Hamiltonian formulation","Eulerian description","coordinate time","Schwarzschild spacetime","generalized vorticity","perturbation stability","black hole accretion"],"falsifier":"Compute the standard covariant conservation law for dust in Schwarzschild coordinates, $\\partial_t(\\sqrt{-g}\\,\\rho U^t) + \\partial_r(\\sqrt{-g}\\,\\rho U^r) = 0$ with $\\sqrt{-g} = r^2\\sin\\theta$, for a stationary radial flow and compare the resulting density profile with the paper's equation $\\partial_r(f\\rho_0 U_0)=0$ and the derived velocity profile; any discrepancy would show that the proposed continuity equation does not reproduce relativistic particle conservation.","tokens_in":12326,"feed_emoji":"🕳️","tokens_out":4412,"duration_ms":49274,"temperature":0.7,"pith_summary":"This paper tries to establish that the familiar Hamiltonian structure of non-relativistic fluid mechanics—where density and velocity fields evolve under Poisson brackets and a Hamiltonian—can be lifted to relativistic fluids in curved spacetimes without abandoning coordinate time. The authors construct a relativistic Eulerian Poisson algebra and a fluid Hamiltonian, and derive generalized continuity and Euler equations that hold for generic static metrics. Specializing to Schwarzschild spacetime, they obtain exact stationary radial inflow and outflow solutions whose coordinate velocity freezes at the horizon, and they show that first-order density and velocity perturbations decay exponentially in time. If correct, this offers a first-principles Hamiltonian route to studying fluid behavior near black holes and other strong-field astrophysical settings, with the observable dynamics expressed from the viewpoint of a static observer at fixed position.","feed_headline":"Fluid dynamics in curved spacetime cast in Hamiltonian form","feed_subtitle":"Density and velocity in ordinary time reproduce relativistic continuity and Euler equations, with stable flows around Schwarzschild.","key_machinery":"The key machinery is the Eulerian map from Lagrangian particle variables to fields, $\\rho(\\vec r,t) = \\rho_0\\int d^3x\\, \\delta^3(\\vec X(\\vec x,t)-\\vec r)$ and $j^i = \\rho U^i$, together with the relativistic Poisson bracket algebra derived from $\\{X^\\mu, U^\\nu\\} = (\\Phi/m)\\,g^{\\mu\\nu}$ and $\\{U^\\alpha, U^\\beta\\}$ containing metric derivatives. The Hamiltonian $H = \\int d^3r\\, \\rho\\left(g_{0i}g^{i\\alpha}U^\\alpha/g_{00} + \\sqrt{\\Psi}\\right)$ generates the equations of motion, and the generalized vorticity $\\Omega^{ij} = g^{ik}\\partial_k U^j - g^{jk}\\partial_k U^i$ organizes the coupling terms in the Euler equation, playing no role in the spherically symmetric background but emerging at first order in perturbations to drive angular dynamics.","core_discovery":"The central discovery is a relativistic extension of the Eulerian Hamiltonian formalism: a Poisson bracket algebra for the density field $\\rho(\\vec r,t)$ and velocity field $U^i(\\vec r,t)$, together with a Hamiltonian, such that Hamilton's equations reproduce relativistic continuity and Euler equations in coordinate time. The algebra is obtained by mapping canonical brackets of Lagrangian variables $(X^\\mu, U^\\nu)$ with metric-dependent brackets $\\{X^\\mu, U^\\nu\\} = (\\Phi/m)\\,g^{\\mu\\nu}$, where $\\Phi\\approx\\sqrt{-g_{00}}$, into field brackets. The proposed Hamiltonian is $H = \\int d^3r\\, \\rho\\left[g_{0i}g^{i\\alpha}U^\\alpha/g_{00} + \\sqrt{\\Psi}\\right]$ with $\\Psi$ built from the spatial metric and velocities, which reduces to the non-relativistic kinetic-plus-potential form in the flat limit. For a diagonal static metric the continuity equation becomes $\\dot\\rho + \\partial_i\\left((-g_{00})\\rho U^i/\\sqrt{1+Z}\\right)=0$ with $Z = g_{ij}U^iU^j$, and the Euler equation acquires a generalized vorticity term $\\Omega^{ij} = g^{ik}\\partial_k U^j - g^{jk}\\partial_k U^i$. Specializing to Schwarzschild, stationary radial background flows satisfy $U_0(r) = \\pm\\sqrt{(1-B/r)(-1+e^{B/r})}$ and $f\\rho_0 U_0 = \\text{constant}$, and first-order perturbations about this background decay exponentially in coordinate time, indicating linear stability.","pith_inferences":["The paper does not check whether its continuity equation matches the covariant conservation law $\\nabla_\\mu(\\rho U^\\mu)=0$ when $\\rho$ is interpreted as rest-mass density; a direct comparison with the standard result would either validate the mapping or show that the density field needs a $\\sqrt{-g}$ factor, and this check is a natural next step.","The vanishing of the coordinate velocity at the horizon is an expected coordinate-time artifact (matter takes infinite coordinate time to reach the horizon), so the framework's predictions about the horizon region should be interpreted with proper-time physics in mind.","A testable extension is to verify whether the full Poisson algebra with $\\Phi = \\sqrt{-g_{00}(1+Z)}$ and the complete Hamiltonian $\\sqrt{\\Psi}$ satisfy the Jacobi identities; if they do, the low-velocity truncation used in the paper could be relaxed to obtain higher-order relativistic corrections.","The formalism's explicit dependence on the metric components suggests a direct route toward rotating (Kerr) spacetimes through the $g_{0i}$ terms in the Hamiltonian, as the authors note in their outlook, and such an extension could be tested against known Kerr accretion solutions."],"forward_implications":["If the framework is correct, the derived bracket algebra and Hamiltonian constitute a new, first-principles Hamiltonian description of relativistic ideal fluids in curved spacetimes, valid for any static metric and reducing to the standard non-relativistic equations in the flat limit.","The coordinate-time continuity equation $\\dot\\rho + \\partial_i\\left((-g_{00})\\rho U^i/\\sqrt{1+Z}\\right)=0$ provides a concrete prescription for evolving fluid density and velocity from the viewpoint of a fixed Eulerian observer, which the authors argue is better aligned with astrophysical observations than proper-time covariant formulations.","The exact stationary radial flow solutions in Schwarzschild, with $U_0$ vanishing at the horizon and $f\\rho_0 U_0 = \\text{constant}$, give explicit background profiles that can serve as starting points for studying black-hole accretion in this Hamiltonian language.","The linear stability result—exponential decay of first-order radial and angular perturbations—implies that the spherically symmetric radial background flow is stable under small perturbations, a property that would support the use of these backgrounds in astrophysical modeling.","The generation of vorticity at linear order over an irrotational background indicates that angular motion can be induced purely by perturbations, a feature that could matter for understanding the onset of non-radial structure in accretion flows."],"supporting_citations":[{"why":"Supplies the non-relativistic Eulerian Poisson algebra and the mapping from Lagrangian variables that the paper extends to curved spacetime; this is the template for the relativistic algebra.","marker":"[19]"},{"why":"Provides the covariant canonical Poisson brackets for a relativistic particle in curved spacetime, which the paper uses as the starting point to derive the field algebra.","marker":"[18]"},{"why":"Gives the particle action and Hamiltonian relation $H = -P_0$ with the dispersion relation used to propose the relativistic fluid Hamiltonian.","marker":"[52]"},{"why":"Defines the Bondi accretion boundary condition $U_0 \\to 0$ as $r \\to \\infty$, which the paper adopts to fix the constant in the stationary background flow solution.","marker":"[28]"},{"why":"Establishes the physical context of spherical accretion onto compact objects, which motivates the spherically symmetric background flow studied in the Schwarzschild application.","marker":"[22]"}],"fun_headline_variants":["Relativistic fluid dynamics as Hamiltonian flow in curved space","New Hamiltonian algebra for relativistic fluids in curved spacetime","Curved spacetime fluid equations from Hamiltonian brackets","Relativistic fluid stability via Hamiltonian mechanics in Schwarzschild"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The derivation assumes that the Eulerian density $\\rho(\\vec r,t) = \\rho_0\\int d^3x\\,\\delta^3(\\vec X - \\vec r)$, with no metric volume element or $\\gamma$ factor, together with the particle Hamiltonian $H = -P_0$ evaluated under $\\Phi\\approx\\sqrt{-g_{00}}$, correctly encodes relativistic matter so that the resulting continuity and Euler equations are the true relativistic equations in curved spacetime.","fun_headline_variants_meta":{"raw":{"variants":["Relativistic fluid dynamics as Hamiltonian flow in curved space","New Hamiltonian algebra for relativistic fluids in curved spacetime","Curved spacetime fluid equations from Hamiltonian brackets","Relativistic fluid stability via Hamiltonian mechanics in Schwarzschild"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000374,"raw_usage":{"total_tokens":2110,"prompt_tokens":1169,"completion_tokens":941,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":785,"completion_tokens_details":{"reasoning_tokens":879}},"tokens_in":785,"tokens_out":941,"duration_ms":10967,"temperature":1.0,"reasoning_tokens":879,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T15:10:01.260186+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the standard covariant conservation law for dust in Schwarzschild coordinates, $\\partial_t(\\sqrt{-g}\\,\\rho U^t) + \\partial_r(\\sqrt{-g}\\,\\rho U^r) = 0$ with $\\sqrt{-g} = r^2\\sin\\theta$, for a stationary radial flow and compare the resulting density profile with the paper's equation $\\partial_r(f\\rho_0 U_0)=0$ and the derived velocity profile; any discrepancy would show that the proposed continuity equation does not reproduce relativistic particle conservation.","supporting_citations":[{"cited_title":"Dynamics of Extended Bodies in General Relativity—I. Momentum and An- gular Momentum,","cited_arxiv_id":null,"evidence_quote":"Provides the covariant canonical Poisson brackets for a relativistic particle in curved spacetime, which the paper uses as the starting point to derive the field algebra."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the particle action and Hamiltonian relation $H = -P_0$ with the dispersion relation used to propose the relativistic fluid Hamiltonian."},{"cited_title":"Bondi; ’On spherically symmetrical accretion’, 1952MNRAS.112..195B","cited_arxiv_id":null,"evidence_quote":"Defines the Bondi accretion boundary condition $U_0 \\to 0$ as $r \\to \\infty$, which the paper adopts to fix the constant in the stationary background flow solution."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Establishes the physical context of spherical accretion onto compact objects, which motivates the spherically symmetric background flow studied in the Schwarzschild application."}],"review_version":1}