{"id":"82e324b1-4714-4197-9d7e-3996bb58a138","arxiv_id":"2506.19477","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":5.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":4,"one_line_summary":"A fluid with intrinsic spin forms accretion tori around Kerr black holes whose size, thickness, and density shift with the spin parameter, in opposite directions for co-rotating and counter-rotating disks.","lead":"This paper calculates equilibrium shapes of thick accretion disks made of a fluid whose particles carry intrinsic spin, around rotating (Kerr) black holes, extending earlier work on non-rotating black holes. It shows that the fluid spin can substantially alter the disk's size, thickness, and peak density, with opposite behavior for co-rotating and counter-rotating disks.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Ansatz (26) for the spin density is the load-bearing assumption; without a sensitivity test the reported morphological trends could be closure artifacts.","rationale":"The reader identified the untested spin-density ansatz as the weakest assumption; I agree. It is the single most load-bearing item because every reported morphological effect (sign and magnitude) is produced by solving the field equations with this closure. Without it, the system would be underdetermined, and with a different closure the results could differ. Other issues, such as absent convergence tests or no data/code release, affect confidence in the numerics but not the logical structure of the argument; they are secondary to the closure question. The paper is transparent about the ansatz and about limitations on s0, so the work is not flawed beyond repair. A conditional verdict with a requirement to test sensitivity is appropriate; I would not reject or upgrade to accept without that check.","tokens_in":18465,"tokens_out":23532,"duration_ms":233453,"concrete_test":"Recompute the fiducial co-rotating model (a=0.95, l0=2.36) for s0 = ±0.03 using the alternative closure S = s0 ϵ^{γ−1} k(r,θ) (dropping the (1−Ωl0) factor in Eq. (26)), deriving the new compatibility PDE for k by the same compatibility procedure and solving (27)–(28) with the corresponding modified spin terms. Then compare the sign and magnitude of Δr_out = r_out − r_in, Δr_max, and Δϵ_max relative to the s0=0 torus against the paper's Figures 2 and 7. If any sign flips or the solution ceases to exist for previously allowed s0, the central morphological claims are closure-dependent; if all signs persist, the concern is largely resolved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"All quantitative claims about how s0 changes torus morphology follow from the spin-density ansatz (26), S = s0 ϵ^{γ−1} k(r,θ)(1−Ωl0). This form is not derived from the Weyssenhoff fluid's conservation laws or from microphysics; it is chosen so that the compatibility PDE (29) for k decouples from ϵ and s0. With γ=2, the spin terms in (27)–(28) become independent of ϵ, and k is determined solely by the background spacetime and l0. If a physically different closure, e.g., S = s0 ϵ^{γ−1} k(r,θ) without the (1−Ωl0) factor or S ∝ ϵ^γ, were adopted, the balance near the cusp and the vertical structure would change, potentially altering the sign, magnitude, or existence of the reported effects. The paper presents no sensitivity analysis to this ansatz, and the abstract's 'demonstrate' overstates what is a single-closure construction.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper constructs stationary, axisymmetric, non-self-gravitating equilibrium tori of a neutral Weyssenhoff spin fluid with constant specific orbital angular momentum in the Kerr spacetime, extending an earlier Schwarzschild construction. The model assumes circular flow, the Frenkel supplementary condition, a polytropic equation of state with gamma=2, and a spin-density ansatz. The authors solve a compatibility PDE for the function k(r,theta), integrate the momentum-balance equations numerically, and report how the fluid spin parameter s0 changes the cusp radius, center radius, outer edge, peak energy density, and isodensity surfaces for both co-rotating and counter-rotating disks. They also map the existence region in the (l0,s0) parameter plane and identify bounds beyond which equilibrium solutions do not exist.","tokens_in":18666,"tokens_out":7270,"duration_ms":72685,"significance":"If the results hold, they provide a useful extension of Polish-doughnut models to matter with macroscopic spin, showing that spin-curvature coupling can appreciably alter torus morphology in a rotating Kerr background and that such coupling imposes existence bounds on the spin parameter. The recovery of the standard spinless torus in the s0=0 limit, the use of the previously derived integrability conditions, and the explicit parameter-space exploration are genuine strengths. The central claim is not circular: s0 and l0 are input parameters, not fitted outputs. However, the quantitative trends are tied to the specific spin-density closure and to an arbitrary normalization of k, so the significance of the morphological claims is conditional on those choices.","major_comments":[{"comment":"All reported morphological trends follow from the postulated spin-density ansatz (26), S(r,theta)=s0 epsilon^{gamma-1} k(r,theta)(1-Omega l0), but the paper neither derives this closure from the Weyssenhoff-fluid conservation laws nor tests the sensitivity of the results to alternative closures. With gamma=2 the spin terms in Eqs. (27)-(28) become independent of epsilon and k is fixed solely by the background and l0, so the sign and magnitude of the changes in r_cusp, r_max, r_out and epsilon_max in Figures 4-8 are properties of this particular closure. Please either derive the ansatz from an underlying condition or add sensitivity runs with a different spin-density form (for example, dropping the (1-Omega l0) factor, or using S proportional to epsilon^gamma), and restate the abstract's 'demonstrate' accordingly.","section":"Sec. III A, Eq. (26)"},{"comment":"The choice k(r,pi/2)=1 is an arbitrary normalization, and the text explicitly states that any other normalization can be reabsorbed into s0. As a result, the bounds on s0 shown in Fig. 9 and advertised in the abstract are not invariant under this rescaling; only the product s0 k(r,theta) enters the physical equations. Please state the normalization convention used and express the parameter-space constraints in a normalization-independent form, or clearly characterize them as convention-dependent results.","section":"Sec. III B"},{"comment":"No convergence tests, resolution studies, or error estimates are reported for the method-of-characteristics solution of the PDE (29) or for the RK4 integrations of Eqs. (27)-(28). Without such tests, the reported changes in r_cusp, r_max, r_out and epsilon_max cannot be distinguished from numerical integration error. Please add at least a resolution study and an independent residual check of the PDE and of the equatorial-plane integration.","section":"Sec. III B and Figs. 2-8"},{"comment":"The existence diagram includes a red-dashed region in which the code returns a solution that violates a boundary condition for epsilon, while the white no-solution region is determined by the same code without any stated tolerance or error analysis. As the parameter-space bounds on s0 are one of the paper's central results, the manuscript should specify how the valid/invalid boundary is detected and report the numerical tolerances used.","section":"Sec. IV B, Fig. 9"}],"minor_comments":[{"comment":"The abstract and conclusions state constraints on s0 without noting that the existence analysis in Sec. IV B is explicitly restricted to co-rotating disks; please add this caveat wherever the constraints are advertised.","section":"Abstract and Sec. IV B"},{"comment":"Minus signs are rendered as '□' in several figure labels and axes (for example, 's0 = □0.01' in Fig. 2, the horizontal axis of Fig. 4, and '□0.045' in Figs. 7-8); these should be replaced with proper minus signs.","section":"Figures 2, 4, 7, 8"},{"comment":"The top-left panel of Fig. 9 contains a stray 'htb!' marker that should be removed.","section":"Fig. 9"},{"comment":"Please correct typographical errors, including 'characterstic' in the abstract, 'whith' in Sec. III B, 'compatify' in Sec. III B, 'paramaters' in Sec. III A, and the missing 's0' in 'parameters 0' in the abstract.","section":"Throughout"},{"comment":"The caption's 'a0 = 0.55' and 'a0 = 0.97' should read 'a = 0.55' and 'a = 0.97'.","section":"Fig. 3 caption"}],"recommendation":"major_revision","confidential_remarks":"The manuscript is an incremental extension of the authors' earlier Schwarzschild work [33] to the Kerr background, and the authors overlap with [33]. The text cites [33] prominently, so this is not a disclosure problem, but the editor may wish to ensure that the contribution is framed as an extension rather than as an independent first construction. The two substantive risks are the untested spin-density closure and the normalization dependence of the reported s0 bounds; both are addressable in revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The short version: this is a serious, well-built extension of the authors' Schwarzschild spin-fluid torus work to Kerr. The new content is the Kerr morphology study and the parameter-space existence diagram; the spinless limit recovers standard tori, so the machinery is consistent. Worth a serious referee.\n\nWhat the paper does well: the derivation is traceable. The integrability conditions are stated, the PDE for the spin-density function k is explicit, and the numerical method is described in enough detail to reproduce. The existence diagram with violation regions marked is a genuinely useful addition. The authors are honest about the limits of their model.\n\nThe main soft spot is the spin-density ansatz, Eq. (26). It is not derived from the fluid's microphysics; it is a closure chosen so the compatibility PDE for k decouples from ε. The paper gives no sensitivity analysis to this closure. The stress-test worry is fair: a different closure, such as dropping the (1−Ωl0) factor or changing the ε dependence, could flip the sign, magnitude, or existence of the reported effects. The spinless limit is a necessary sanity check but does not constrain the spin closure, because the spin terms vanish there. So the central claim, that intrinsic spin changes torus morphology, is plausible in direction but not quantitatively pinned down.\n\nMinor issues: no convergence tests or error estimates for the RK4 integrations, and no code or data released, so the numerics are not independently checkable. The abstract also overstates by saying 'demonstrate'; 'indicate' would be more accurate. Some typos ('characterstic') should be fixed.\n\nWho this is for: researchers building initial data for GRMHD simulations, and people working on spin-fluid phenomenology in strong gravity. It is a niche but relevant contribution.\n\nRecommendation: send to peer review. Ask the authors to test sensitivity to the ansatz, report convergence checks, and ideally release the code or data. With those, the quantitative claims would be much more trustworthy.","headline":"A solid Kerr extension of spin-fluid tori whose quantitative claims rest on an untested spin-density ansatz; deserves refereeing with a request for sensitivity analysis.","tokens_in":19198,"tokens_out":3018,"would_cite":false,"duration_ms":31061,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["83C57","83C55","83C10"],"pacs":[],"model":"deepseek-v4-flash","headline":"Macroscopic fluid spin reshapes accretion tori around Kerr black holes.","keywords":["accretion tori","Weyssenhoff spin fluid","Kerr black holes","spin-curvature coupling","equilibrium models","constant specific angular momentum","relativistic hydrodynamics"],"falsifier":"Compute the equilibrium spin-density distribution from a microscopic model of spinning particles with the Frenkel condition in the Kerr background, insert it into the momentum balance equation, and check whether the resulting tori reproduce the predicted shifts of $r_{\\rm max}$, $r_{\\rm cusp}$, and $\\epsilon_{\\rm max}$ with $s_0$; a sign reversal or an order-one deviation in magnitude would falsify the ansatz-based results.","tokens_in":1866,"feed_emoji":"🕳️","tokens_out":2412,"duration_ms":66966,"temperature":0.7,"pith_summary":"The paper constructs stationary, axisymmetric equilibrium models of geometrically thick accretion tori made of a neutral Weyssenhoff spin fluid with constant specific orbital angular momentum around Kerr black holes. It claims that macroscopic fluid spin, encoded in a parameter $s_0$, changes the torus morphology in a systematic way: for co-rotating disks negative $s_0$ enlarges the torus and raises its peak energy density while positive $s_0$ compacts it, and counter-rotating disks show the opposite pattern, with the effects strengthening as the black-hole spin $a$ grows. The authors also map the allowed $(l_0, s_0)$ parameter space and find bounds beyond which equilibrium solutions cease to exist. A sympathetic reader would care because accretion torus models are used as initial data and interpretive tools for observations of supermassive black holes, and intrinsic spin is one physical ingredient usually omitted.","feed_headline":"Fluid spin can enlarge or shrink black-hole accretion tori","feed_subtitle":"Co-rotating fluid spin grows the torus; counter-rotating spin shrinks it, and faster for high-spin Kerr holes.","key_machinery":"The load-bearing object is the spin-density ansatz, Eq. (26): $S(r,\\theta)=s_0\\,\\epsilon^{\\gamma-1} k(r,\\theta)(1-\\Omega l_0)$, with $s_0$ a constant spin parameter and $k$ a shape function fixed by a compatibility PDE, Eq. (29), that is independent of $\\epsilon$ and normalized to $k(r,\\pi/2)=1$. This ansatz makes the mixed-derivative compatibility condition for $\\epsilon$ independent of $\\epsilon$, so $k$ can be solved first by the method of characteristics and the energy density then integrated outward from the cusp. The other essential piece is the integrability framework inherited from the Schwarzschild predecessor, which reduces the momentum balance to an effective-potential equation $W-W_{\\rm in}=\\ln|u_t|-\\ln|u_{t\\rm in}|-\\int \\Omega\\,dl/(1-\\Omega l)$; all of the reported morphological changes trace back to $s_0$ acting through the geometry-fixed characteristic curves of $k$.","core_discovery":"The central discovery is that a neutral Weyssenhoff spin fluid, a continuum fluid whose elements carry intrinsic spin angular momentum, can form closed, non-self-gravitating tori with constant specific angular momentum in the Kerr spacetime, and that the spin-curvature coupling term systematically deforms them. For co-rotating motion, negative fluid spin $s_0$ pushes the density maximum $r_{\\rm max}$ outward, moves the cusp inward, enlarges the radial extent and vertical thickness, and increases the peak energy density; positive $s_0$ does the reverse. Counter-rotating motion inverts the trend. The shifts grow with the Kerr spin parameter and are already visible at $a=0.55$ and pronounced at $a=0.95$. The solutions are built by solving the momentum balance equation with a polytropic equation of state ($\\gamma=2$), a spin-density ansatz that decouples the shape function $k(r,\\theta)$ from the energy density, and the method of characteristics; closed tori exist only in a bounded region of the $(a,l_0,s_0)$ parameter space.","pith_inferences":["I infer that the co- versus counter-rotating sign asymmetry could serve as an observational diagnostic: if real accretion flows carry net macroscopic spin, horizon-scale images of galactic-center black holes should show systematic differences between disks rotating with and against the hole.","A testable extension would be to derive $S(r,\\theta)$ from a kinetic or field-theoretic spin-fluid model instead of postulating the ansatz; if the derived distribution differs, the morphological trends could change sign or magnitude.","I infer from the parameter-space maps that the disappearance of solutions at high $|s_0|$ may signal a real physical limit, spin-curvature forces overpowering pressure gradients, rather than a numerical artifact, meaning spin fluids could be unable to form equilibrium tori around rapidly spinning black holes beyond a threshold."],"forward_implications":["If the claim holds, general-relativistic magnetohydrodynamic simulations around Kerr black holes can be initialized with spin-fluid tori whose size, thickness, cusp position, and density peak depend on $s_0$ as well as on angular momentum and black-hole spin.","For co-rotating disks, negative $s_0$ moves the cusp closer to the horizon and pushes the density maximum outward, which changes where accretion begins and how much material sits at large radius, with direct consequences for the predicted emission region.","For counter-rotating disks the sign flips, so the same fluid spin parameter produces opposite morphological signatures; comparing tori that rotate with and against the black hole could in principle constrain $s_0$.","The existence bounds imply that spin-fluid equilibrium tori exist only in a finite region of $(l_0,s_0)$, so simulations with too large $|s_0|$ have no hydrostatic starting configuration.","Near the marginally stable orbit, spin can make a finite torus where the spinless model gives a point, and near the marginally bound orbit it can make a point where the spinless model gives an infinite disk, so spin changes which constant-$l$ tori are allowed."],"supporting_citations":[{"why":"Provides the integrability conditions, the spin-density ansatz, and the Schwarzschild equilibrium solutions that this work extends to Kerr.","marker":"[33]"},{"why":"Supplies the symmetric energy-momentum tensor and momentum balance equation for the ideal Weyssenhoff spin fluid used as the master equation.","marker":"[29]"},{"why":"Defines the Weyssenhoff spin-fluid model and the Frenkel supplementary condition on which the spin tensor components are based.","marker":"[18]"},{"why":"Documents the spin-curvature shift of circular orbits for spinning particles, used to anticipate and bound the modified $r_{\\rm ms}$ and $r_{\\rm mb}$.","marker":"[34]"},{"why":"Establishes the geometrically thick torus (Polish doughnut) equilibrium framework for non-self-gravitating disks that the spin fluid generalizes.","marker":"[16]"}],"fun_headline_variants":["Spin changes black-hole torus size and shape","Intrinsic fluid spin reshapes accretion tori","Co-rotating spin inflates tori, counter-rotating deflates","Kerr tori stretch or shrink with fluid spin","Fluid spin deforms Kerr accretion tori"],"cache_read_input_tokens":21376,"weakest_assumption_plain":"Everything rests on the assumed spin-density form $S(r,\\theta)=s_0\\,\\epsilon^{\\gamma-1} k(r,\\theta)(1-\\Omega l_0)$, chosen so the equations decouple rather than derived from microphysics; if a real spin fluid distributes its spin differently, the reported enlargement, shrinkage, and even the existence of the tori could change.","fun_headline_variants_meta":{"raw":{"variants":["Spin changes black-hole torus size and shape","Intrinsic fluid spin reshapes accretion tori","Co-rotating spin inflates tori, counter-rotating deflates","Kerr tori stretch or shrink with fluid spin","Fluid spin deforms Kerr accretion tori"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000248,"raw_usage":{"total_tokens":1626,"prompt_tokens":1102,"completion_tokens":524,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":718,"completion_tokens_details":{"reasoning_tokens":445}},"tokens_in":718,"tokens_out":524,"duration_ms":5073,"temperature":1.0,"reasoning_tokens":445,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-15T18:32:29.349828+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the equilibrium spin-density distribution from a microscopic model of spinning particles with the Frenkel condition in the Kerr background, insert it into the momentum balance equation, and check whether the resulting tori reproduce the predicted shifts of $r_{\\rm max}$, $r_{\\rm cusp}$, and $\\epsilon_{\\rm max}$ with $s_0$; a sign reversal or an order-one deviation in magnitude would falsify the ansatz-based results.","supporting_citations":[{"cited_title":"Stationary equilibrium torus supported by Weyssenhoff ideal spin fluid in Schwarzschild spacetime -- I: Case of constant specific angular momentum distribution","cited_arxiv_id":"2307.16292","evidence_quote":"Provides the integrability conditions, the spin-density ansatz, and the Schwarzschild equilibrium solutions that this work extends to Kerr."},{"cited_title":"Obukhov and O","cited_arxiv_id":null,"evidence_quote":"Supplies the symmetric energy-momentum tensor and momentum balance equation for the ideal Weyssenhoff spin fluid used as the master equation."},{"cited_title":"Weyssenhoffand A","cited_arxiv_id":null,"evidence_quote":"Defines the Weyssenhoff spin-fluid model and the Frenkel supplementary condition on which the spin tensor components are based."},{"cited_title":"Abramowicz, M","cited_arxiv_id":null,"evidence_quote":"Establishes the geometrically thick torus (Polish doughnut) equilibrium framework for non-self-gravitating disks that the spin fluid generalizes."}],"review_version":2}