{"id":"63d23339-4973-4ae0-9bd2-b898c0c2eb59","arxiv_id":"2508.20852","paper_version":1,"verdict":"REJECT","confidence":"HIGH","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":0,"one_line_summary":"Claims a general curvature-based expression for the streaming operator in curvilinear radiative transport, but the main derivation drops a nonzero shape-operator term, so the general formula is wrong.","lead":"The paper derives a geometric formula for the streaming (transport) term in radiative transfer when space and angle are described in curvilinear coordinates, expressing directional derivatives through the curvature of surfaces and curves. The intended payoff is a more intuitive way to pick coordinate systems for particle transport problems, but the general derivation contains an algebraic error that invalidates it outside special symmetric cases.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Eq (27) is not a valid general streaming formula: the derivation of ∇_Ω μ drops the off-diagonal shape-operator term at Eq (17) and uses the wrong sign for the μ√(1−μ²)(Ω̂·κ_n) term, since Remark 2 defines κ_n = −∇_n n.","rationale":"The reader correctly identifies the Eq (17) cancellation as a real algebraic error; it breaks the equality between the curve-curvature and surface-curvature forms and invalidates the examples that use the curve-curvature form. However, I find an even more direct falsification of the headline Eq (27): the sign of the κ_n term. Substituting Ω∥·∇_n n from Eq (13) into Remark 2's κ_n=−∇_n n gives a minus sign; Eq (16) and Eq (27) use plus. Because κ_n is nonzero in the paper's own translating-paraboloid and ellipsoid examples, this is not a harmless typo in a zero term. Both errors must be fixed before the claimed generality can be assessed; as submitted the central claim is unsupported. The test above is a single finite-difference check that exposes both the sign and cross-term omissions simultaneously.","tokens_in":17805,"tokens_out":20197,"duration_ms":186960,"concrete_test":"Using the paper's translating-paraboloid data with a=1,b=2 at (x,y)=(1,1) and a fixed Ω=(n+t)/|n+t|, numerically evaluate Ω·∇(Ω·n) by finite differences along r0+sΩ. Compare with the RHS of Eq (21) as printed, and with the RHS obtained by (i) keeping the cross term 2(1−μ²)sinω cosω S_tb and (ii) changing the κ_n term to −μ√(1−μ²)(Ω̂·κ_n). The corrected expression matches; Eq (21) as printed does not.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central formula Eq (27) is false as stated. In Section 2.1 the derivation of ∇_Ω μ makes two algebraic errors. First, Eq (17) cancels the sinω cosω term by asserting t·∇_b n + b·∇_t n = 0 via Remark 1. That cancellation is invalid: Remark 1 only gives t·∇_b n = −n·∇_b t and b·∇_t n = −n·∇_t b, which are not negatives of one another. For a surface with normal n and tangent frame (t,b), the Weingarten map is symmetric, so t·∇_b n = b·∇_t n = S_tb; the cross term is 2(1−μ²)sinω cosω S_tb and vanishes only if t,b are principal directions. Integrability of the frame does not imply this, and the paper's own translating-paraboloid example has nonzero S_tb,S_bt (Eqs 56–57). Second, Eq (16) replaces ∇_n n by κ_n even though Remark 2 defines κ_n = −∇_n n; substituting Ω_∥·∇_n n into Eq (13) yields a minus sign in front of μ√(1−μ²)(Ω̂·κ_n). Eq (21)/(27) carry a plus sign. A direct computation of Ω·∇(Ω·n) for the translating paraboloid (e.g., a=1,b=2, x=y=1, Ω along n+t) matches only after flipping that sign. Hence the claimed 'general' streaming expression is not general: it needs the missing S_tb term and the corrected sign, and every example using the curve-curvature form or the κ_n term inherits the error.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The manuscript derives a general expression for the streaming operator Ω·∇Ψ in transport equations when the direction Ω is parameterized by (μ,ω) in a local orthonormal frame (n,t,b). It proposes Eq. (27) as a universal formula involving the surface curvature C(r,ω), the curvature κ_n of the integral curve of n, the curvatures κ_t,κ_b, and a winding term t·∇_n b. The derivation is split into Ω·∇μ and Ω·∇ω, with curve-curvature and surface-curvature forms, and is illustrated on cylindrical, spherical, ellipsoidal, and translating-graph coordinate examples.","tokens_in":18274,"tokens_out":6600,"duration_ms":62037,"significance":"The geometric framework is appealing and, if correct, would be a useful reference for practitioners choosing coordinate systems in kinetic/radiative transfer. The paper also makes a pedagogical contribution by expressing angular derivative terms through curvature of integral curves and shape operators, and the examples are worked in detail. However, the central formula is invalid as stated: the derivation of ∇Ωμ contains a missing off-diagonal shape-operator term and a sign error. These errors change the streaming operator in all but trivial geometries, so the paper's main claim of a general streaming expression is not supported.","major_comments":[{"comment":"The derivation cancels the term sin(ω)cos(ω)(t·∇_b n + b·∇_t n) by invoking Remark 1. This is incorrect: Remark 1 gives t·∇_b n = −n·∇_b t and b·∇_t n = −n·∇_t b, which are not negatives of each other. For a surface with normal n and orthonormal tangents t,b, Weingarten symmetry yields t·∇_b n = b·∇_t n = S_tb, so the cross term is 2(1−μ²)sinω cosω S_tb. It vanishes only if t,b are principal directions. Integrability of the frame does not imply this. The paper's own translating-paraboloid frame has nonzero S_tb and S_bt (Eqs. (56)–(57)), so the cancellation is invalid in general and in the paper's own example. Consequently Eqs. (18), (21), and (27) are missing this term.","section":"§2.1, Eq. (17)"},{"comment":"The sign convention in Remark 2 defines κ_u := −∇_u u. In Eq. (16) the term μ(Ω∥·∇_n n) is replaced by μ√(1−μ²)(Ω̂∥·κ_n). This has the wrong sign: ∇_n n = −κ_n, so Eq. (16) should be −μ√(1−μ²)(Ω̂∥·κ_n). The same sign error propagates into Eqs. (21) and (27). A direct check on the translating paraboloid confirms that the sign must be flipped for the formula to reproduce Ω·∇(Ω·n).","section":"§2.1, Eqs. (15)–(16), (21)"},{"comment":"The translating-paraboloid example is advertised as an illustration of the general framework, but it explicitly violates the assumption needed for the Eq. (17) cancellation: the off-diagonal shape-operator entries S_tb and S_bt are nonzero (Eqs. (56)–(57)). At (x,y)=(1,1) with a=1,b=2, the sum t·∇_b n + b·∇_t n is nonzero, so the sinω cosω term does not vanish. Thus the example cannot be used to validate Eq. (27); it instead demonstrates the missing term.","section":"§3.5, Eqs. (56)–(57)"}],"minor_comments":[{"comment":"Typos: 'Curvelinear coordiantes' should be 'Curvilinear coordinates'.","section":"Keywords"},{"comment":"The sentence 'All that is left is to plug these terms back into the formulas for ∇ΩΨ, which we will refrain from due to the length of the terms' is awkwardly phrased; the authors could either provide the final expression or explicitly say it is omitted for brevity.","section":"§3.4, text after Eq. (48)"},{"comment":"The surface-curvature form of ∇Ωω is introduced without a fully detailed derivation, and the notation (1,0)S_Xb(...) is hard to parse. Since this form is not used in the final examples, the authors should either clarify it or remove it.","section":"§2.2, Eqs. (25)–(26)"},{"comment":"There are several typographical inconsistencies, e.g., 'Froebenius' for 'Frobenius', and some equation numbers are referenced imprecisely (e.g., 'equation (28)' in Section 3 where Eq. (27) is meant).","section":"General"}],"recommendation":"reject","confidential_remarks":"The errors identified in the derivation are fundamental: the missing off-diagonal shape-operator term and the sign error in the κ_n term invalidate the central formula Eq. (27). Because the paper's main contribution is that general streaming expression, a rejection is appropriate. A corrected derivation, if achievable, would require reworking all examples and the concluding claims, which goes beyond minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Punchline: the paper's main result, Eq (27), is not a valid general streaming formula. There are two algebraic errors in Section 2.1. First, the cross term sin(ω)cos(ω)(t·∇_b n + b·∇_t n) in Eq (17) is cancelled by asserting it vanishes via Remark 1. Remark 1 only gives t·∇_b n = −n·∇_b t and b·∇_t n = −n·∇_t b, which are not negatives of one another. For a surface with normal n the Weingarten map is symmetric, so t·∇_b n = b·∇_t n = S_tb, and the sum is 2S_tb—zero only when t and b are principal directions. Integrability does not imply that, and the paper's own translating-paraboloid shape operator (Eqs 56–57) has S_tb ≠ 0. Second, Eq (16) identifies ∇_n n with κ_n, but Remark 2 defines κ_n = −∇_n n, so the µ√(1−µ²)(Ω̂_∥·κ_n) term in (21)/(27) carries the wrong sign. A direct check on the paraboloid confirms both problems. Since these errors enter the central derivation, the ellipse and translating-graph examples inherit them.\n\nCredit where due: the geometric repackaging is genuinely useful. Framing the streaming terms through shape operators, integral-curve curvatures, and the boxed vanishing questions in (28)–(30) gives practitioners a quick way to see which angular derivatives survive. The holonomy argument in A.1 for when ∂ωΨ = 0 is a nice addition. The cylindrical and spherical examples are correct, reproduce the classical streaming operators, and work as sanity checks. The paper honestly positions itself relative to Freimanis (2011). The homothetic-scaling remark and the z-independence of the translating-graph terms are nice practical observations.\n\nSoft spots, in proportion: the two errors above are load-bearing, so the claimed general result is false as stated. The appendix's holonomy classification of surfaces looks garbled—flat surfaces have trivial holonomy, not SO(2)—but that is a side remark. There are typos and the example formulas are hard to check, but the issue is not verification; the derivation itself is wrong.\n\nWho this is for: transport practitioners in nontrivial geometries will find the interpretive framework appealing, but they should not trust Eq (27). The errors are fixable—add the S_tb cross term, correct the κ_n sign, re-derive the non-classical examples—and the framework is worth preserving. I would send this to referees rather than desk-reject; a good referee should catch exactly these issues, and the paper would come back substantially better.","headline":"The central streaming formula (Eq 27) is not correct—the derivation drops a nonzero off-diagonal shape-operator term and mis-signs the κ_n term—but the geometric framework and the classical examples are solid enough that the paper deserves a salvageable referee round.","tokens_in":18749,"tokens_out":17215,"would_cite":false,"duration_ms":134072,"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":"A single curvature formula gives the streaming operator in any curvilinear transport frame","keywords":["radiative transport","streaming operator","curvilinear coordinates","orthonormal frame","curvature","shape operator","integral curves","foliation"],"falsifier":"For the translating-paraboloid frame with a≠b, compute t·∇_b n + b·∇_t n directly at a generic point; the sum is nonzero, so repeating the derivation from Eq. (17) without discarding the sin(ω)cos(ω) term changes the coefficient of ∂Ψ/∂μ. Direct numerical differentiation of Ω·∇μ along the direction Ω̂_∥ would then show the missing cross term.","tokens_in":17679,"feed_emoji":"📐","tokens_out":10866,"duration_ms":102793,"temperature":0.7,"pith_summary":"The paper derives a general expression for the streaming term Ω·∇Ψ in radiative-transport and kinetic equations when the direction Ω is parametrized by a projection µ onto a spatially varying unit vector n and an angle ω in the orthogonal plane spanned by t and b. The expression rewrites the angular derivatives in geometric terms: the drift in µ is controlled by the curvature of a surface normal to n and by the curvature of the integral curves of n, while the drift in ω is controlled by how t and b bend and twist. If correct, this gives practitioners one master formula for writing the streaming operator in any curvilinear coordinate system, replacing the case-by-case geometric derivations common in the literature. It also makes term-by-term vanishing visible: flat surfaces, straight integral curves, geodesics, and non-twisting frames all eliminate terms without further computation.","feed_headline":"Curvature alone fixes streaming in any curvilinear transport frame","feed_subtitle":"General expression ties angular drift to surface and curve curvature, letting practitioners spot zero terms by geometry.","key_machinery":"The machinery is the parametrization Ω(μ,ω)=μn+√(1−μ²)(cos ω t+sin ω b), with (n,t,b) an integrable orthonormal frame allowed to vary in space. Treating μ and ω as spatial fields converts Ω·∇Ψ into a chain rule, and the coefficients are identified with the shape operator of the foliation surface X^n and with curvature vectors κ_n, κ_t, κ_b of the frame's integral curves. Orthonormal-frame identities reduce the μ-coefficient to a squared cosine/sine combination, and the homothetic-scaling remark shows why all curvature terms scale as 1/ρ for self-similar families of surfaces.","core_discovery":"The paper's central claim is Eq. (27), a general formula for the streaming operator when Ω is written in a local orthonormal frame (n,t,b): Ω·∇Ψ equals the spatial derivative Ω·∇_r Ψ, plus a term multiplying ∂Ψ/∂μ built from the normal curvature C(r,ω) of the surface orthogonal to n and the curvature κ_n of the integral curves of n, plus a term multiplying ∂Ψ/∂ω built from the curvatures κ_t, κ_b and the twisting t·∇_n b of the frame around γ_n. Two interchangeable forms are given, one in terms of curve curvatures and one in terms of shape operators, and the paper demonstrates the formula on cylindrical, spherical, elliptical, and translated-graph coordinates, including a homothetic-scaling","pith_inferences":["If the frame's t and b directions are not principal directions of X^n, the sin(ω)cos(ω) cross term kept at Eq. (17) is nonzero; Eq. (27) would then need an extra correction term, so the practical domain of the formula is frames adapted to the surface curvature.","The curve-curvature and shape-operator forms of the angular coefficients can be seen as two gauges for the same geometric drift; switching between them may help build the conservative discretizations the appendix shows are generically unavailable outside planes and spheres.","The translated-graph examples suggest a practical recipe for layered or extruded geometries: since every streaming coefficient depends only on the base graph G0 and not on the translation height z, coefficients can be precomputed once from the cross-section profile.","The appendix's holonomy condition reframes dimensional reduction in ω as a parallel-transport property of the orthogonal foliation, which could guide coordinate selection for kinetic models with toroidal or helical symmetry."],"forward_implications":["In any frame where the surfaces X^n are flat and the integral curves of n are straight, the coefficient of ∂Ψ/∂μ vanishes, simplifying or removing the μ-dependence of the streaming operator.","For frames whose t and b curves are geodesics of X^n and do not twist around γ_n, the entire ∂Ψ/∂ω term drops; this is the geometric reason the sphere example reduces to a single sin(ω)cot(θ)/ρ term.","The formula turns a tedious algebraic chore into a geometric check: each streaming coefficient can be diagnosed as zero from flatness, straightness, geodesy, or non-twisting before any curvature is computed.","For homothetic foliations such as spheres, cylinders, and ellipsoids, all curvature terms scale as 1/ρ, so the streaming coefficients can be computed once at one scale and rescaled.","When X^n are level sets of Ψ with holonomy SO(2), the appendix argument shows ∂Ψ/∂ω vanishes, giving the standard reduced transport equation in μ and space only."],"supporting_citations":[{"why":"Derives the same chain-rule splitting for restricted frame choices; the paper's general-frame result is positioned against it.","marker":"Freimanis (2011)"},{"why":"Classical source for the cylindrical streaming term used as the first example and for the conservation-form discussion.","marker":"Bell and Glasstone (1970)"},{"why":"Example of a helical-symmetry transport derivation that the general formula is meant to subsume.","marker":"Larsen (2007)"},{"why":"Weingarten theorem expressing the shape operator through directional derivatives of the normal, used to identify surface-curvature terms.","marker":"Weingarten (1861)"},{"why":"Textbook treatment of the shape operator and Weingarten map used for the surface-curvature interpretation.","marker":"Kreyszig (2013)"},{"why":"Frobenius theorem cited for the foliation condition V·rot V = 0 that guarantees surfaces orthogonal to a vector field.","marker":"Frobenius (1877)"},{"why":"Modern statement of Frobenius integrability used for the same foliation-existence condition.","marker":"Lee (2000)"},{"why":"Holonomy-group classification cited in the appendix to support the condition under which ∂Ψ/∂ω vanishes.","marker":"Berger (1955)"}],"fun_headline_variants":["Streaming term tied to curvature in curved coordinates","Curvature terms shape angular drift in transport","General streaming formula for curvilinear transport","Curvature reveals when angular terms simplify","Choose coordinates by curvature to simplify streaming"],"cache_read_input_tokens":2688,"weakest_assumption_plain":"The load-bearing step at Eq. (17) assumes t·∇_b n + b·∇_t n = 0, equivalent to t and b being principal directions of the surface X^n; the paper does not state this, the translating-paraboloid example with a≠b violates it, and if it fails Eq. (27) misses the sin(ω)cos(ω) cross term.","fun_headline_variants_meta":{"raw":{"variants":["Streaming term tied to curvature in curved coordinates","Curvature terms shape angular drift in transport","General streaming formula for curvilinear transport","Curvature reveals when angular terms simplify","Choose coordinates by curvature to simplify streaming"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000857,"raw_usage":{"total_tokens":3516,"prompt_tokens":659,"completion_tokens":2857,"prompt_tokens_details":{"cached_tokens":256},"prompt_cache_hit_tokens":256,"prompt_cache_miss_tokens":403,"completion_tokens_details":{"reasoning_tokens":2791}},"tokens_in":403,"tokens_out":2857,"duration_ms":23480,"temperature":1.0,"reasoning_tokens":2791,"cache_read_input_tokens":256,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-05T14:47:34.171613+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"For the translating-paraboloid frame with a≠b, compute t·∇_b n + b·∇_t n directly at a generic point; the sum is nonzero, so repeating the derivation from Eq. (17) without discarding the sin(ω)cos(ω) term changes the coefficient of ∂Ψ/∂μ. Direct numerical differentiation of Ω·∇μ along the direction Ω̂_∥ would then show the missing cross term.","supporting_citations":[{"cited_title":"Sur les groupes d'holonomie homog\\`ene des vari\\'et\\'es \\`a connexion affine et des vari\\'et\\'es riemanniennes","cited_arxiv_id":null,"evidence_quote":"Holonomy-group classification cited in the appendix to support the condition under which ∂Ψ/∂ω vanishes."}],"review_version":1}