{"id":"3e98460d-1c78-42fd-afc0-27c086aef233","arxiv_id":"2507.13316","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":2,"one_line_summary":"The 3D-1D Darcy-Poiseuille blood perfusion model is proven to converge to a new 1D Green's function model at rate ε^{1/2}|log ε| as the vessel radius ε→0.","lead":"This paper proves that a simple one-dimensional (1D) 'slender body' model of blood perfusion around a thin vessel is the rigorous ε→0 limit of a more detailed 3D-1D Darcy-Poiseuille model, with an explicit error rate. It also derives a new 1D integrodifferential equation for blood pressure along the vessel and proves the estimates needed to make the convergence argument work.","discovery_kind":"extension","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Proof of Lemma 3.1's near-positivity is invalid as written: (3.17) does not imply diagonal dominance, so the coercivity of B in Theorem 1.2 lacks support.","rationale":"I read the manuscript in good faith: the overall architecture is coherent, the kernel estimates in Section 2 are largely adapted from [51], and the final error estimate in Section 5 closes arithmetically if Theorem 1.2 and Lemma 4.1 are accepted. However, the central existence theory for the 1D integrodifferential equation rests on Lemma 3.1, and the proof of that lemma contains a quantitative error at equation (3.17). The claimed inequality has the correct direction but compares the symmetrized kernel to a quantity much larger than the square root of the diagonal product; therefore diagonal dominance does not follow. This is different from the reader's fragility concern about non-spheroidal tips: even for the admissible spheroidal geometry, the proof as written fails. The statement of Lemma 3.1 may still be true — a Fourier-based positive-definiteness argument might rescue it — but the manuscript does not supply that argument. I also note that Appendix A defers boundedness and coercivity of the 3D-1D bilinear form to [56, Section 4], which is an additional missing foundation; however, the single most load-bearing internal gap is the near-positivity proof, because it underpins existence, uniqueness, and all bounds for pSB. Since the issue is a repairable proof gap rather than an obvious contradiction, rejection is not warranted; the reader's CONDITIONAL verdict stands unchanged.","tokens_in":42822,"tokens_out":33643,"duration_ms":401774,"concrete_test":"Check the algebra in (3.17) against (3.8). On the diagonal, KY(τ, φϵ(τ)) = η/(4π ϵ a*(τ)), hence the product of square roots is η/(4π ϵ)(a*(τ)a*(t))^{-1/2}, not η/(8π)/(ϵ² a*(τ)a*(t)). Then determine whether the intended diagonal-dominance can be replaced by a direct proof that the symmetrized kernel is positive definite on [-Lϵ, Lϵ]² — for example, by showing 1/|x| positive definiteness survives the asymmetric regularization z_i = Y(t_i) + ϵa*(t_i)e_r(t_i, θ) after θ-averaging. If no such replacement is supplied, Lemma 3.1 must be regarded as unproved.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The decisive step in the proof of Lemma 3.1 is the claim that the symmetrized kernel is diagonally dominant. Equation (3.17) asserts that 1/2(KY(τ, φϵ(t)) + KY(t, φϵ(τ))) ≤ η/(8π) · 1/(ϵ² a*(τ)a*(t)) = KY(τ, φϵ(τ))^{1/2} KY(t, φϵ(t))^{1/2}. But from the definition (3.8), KY(τ, φϵ(τ)) = η/(8π²) ∫₀^{2π} dθ/(ϵa*(τ)) = η/(4π ϵ a*(τ)), so the square-root of the diagonal product is η/(4π ϵ)(a*(τ)a*(t))^{-1/2}. The displayed right-hand side exceeds the required diagonal-dominance quantity by a factor 1/(2ϵ√(a*(τ)a*(t))) ≫ 1, so it cannot establish Iϵ,0 ≥ 0. Consequently the asserted positivity (3.18), the lower bound (3.13), the coercivity (3.78) of B, the well-posedness and bounds of Theorem 1.2, and ultimately the rate in Theorem 1.3 lose their stated support. The lemma may be repairable by a different positive-definiteness argument, but as written the proof has a gap at a load-bearing point.","agreement_with_reader":"partial"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a 3D-1D Darcy-Poiseuille system for blood perfusion around a thin vessel and derives from it a 1D slender-body Green's function model in which the interior pressure solves a degenerate integrodifferential equation. The main results are Theorem 1.1 (well-posedness of the 3D-1D system, with proof deferred to an appendix), Theorem 1.2 (well-posedness and weighted a priori estimates for the 1D equation), and Theorem 1.3 (an ε^{1/2}|log ε| error estimate between the two models). Sections 2 and 4 collect integral estimates and compute residuals for the 1D approximation, and Section 5 proves the error estimate assuming Theorem 1.2. Numerical experiments in Section 1.5 illustrate the behavior of the 1D model. The central claim of the paper is the convergence of the 3D-1D solution to the 1D slender-body solution, with the rate controlled by estimates for the 1D integrodifferential equation.","tokens_in":42983,"tokens_out":9050,"duration_ms":94903,"significance":"If the proof chain is completed, this would be a valuable contribution: it gives a rigorous, quantitative sense in which a practical 1D Green's function model is the slender limit of a 3D-1D Darcy-Poiseuille system, and the a priori estimates in Theorem 1.2 are designed to be reusable in Part II for comparison with a 3D-3D Darcy-Stokes system. The paper is careful and detailed in its residual calculations and in the derivation of the 1D equation, and the inclusion of numerical examples helps orient the reader. The main limitations are that the proof of the key near-positivity lemma (Lemma 3.1) contains a serious gap, and that the well-posedness of the 3D-1D system is not proved in this manuscript but is quoted from the companion paper [56].","major_comments":[{"comment":"The proof of Lemma 3.1 fails at the claimed diagonal-dominance step. From (3.8), KY(τ, φϵ(τ)) = η/(4πϵa*(τ)), so the geometric mean of the diagonal entries is η/(4πϵ) (a*(τ)a*(t))^{-1/2}. The expression on the right-hand side of (3.17), η/(8π) 1/(ϵ²a*(τ)a*(t)), is not equal to this geometric mean; it is larger by a factor proportional to 1/(2ϵ√(a*(τ)a*(t))). Moreover, an entrywise bound of the form |K(τ,t)| ≤ √(K(τ,τ)K(t,t)) does not imply positive definiteness of the symmetric kernel, nor does it imply diagonal dominance. Therefore the inference to (3.15) and to Iϵ,0 ≥ 0 is unsupported. Since Lemma 3.1 is used directly in the coercivity estimate (3.78) and again in (3.93)-(3.94), the well-posedness and all weighted bounds of Theorem 1.2, and consequently the rate in Theorem 1.3, lose their stated support.","section":"§3.1.1, Eq. (3.17)"},{"comment":"Theorem 1.1 is announced with a proof in Appendix A, but the appendix does not actually prove the central coercivity of the bilinear form B in (A.5). Instead it states, without proof, that by [56, section 4] the form is bounded and coercive. Since Theorem 1.3 compares the 3D-1D solution (q,p), whose existence is asserted by Theorem 1.1, with the 1D model, the main convergence claim is conditional on an unpublished companion paper unless the coercivity proof is included or the dependency is explicitly disclosed as an assumption. This is a load-bearing missing argument, not a stylistic issue.","section":"Appendix A / Theorem 1.1"},{"comment":"The bound (3.31), stated as |KY(τ,t)| ≥ C/ϵ² for all τ,t ∈ [-1,1], has the wrong direction. Near the diagonal the kernel behaves like C/(ϵa*(τ)), which for a*(τ) near 1 is larger than C/ϵ² only in a different scaling regime; more importantly, the subsequent estimate of Iϵ,3 in (3.32) requires an upper bound |KY(τ,t)| ≤ C/ϵ². As written, (3.31) cannot justify the step that follows. This is a separate defect in the proof of Lemma 3.1, and it reinforces that the lemma's proof needs revision.","section":"§3.1.1, Eq. (3.31)"}],"minor_comments":[{"comment":"The notation [(a⁴v_t)_t]^* is used before it is defined; the defining sentence should be moved earlier to avoid confusion.","section":"§3.2, Eq. (3.72)"},{"comment":"The numerical experiments are descriptive and useful for illustrating the model, but they do not include a convergence study in ϵ or a direct comparison with the 3D-1D solution. A quantitative test of the rate in Theorem 1.3 would strengthen the paper.","section":"§1.5"},{"comment":"The constant d03 = 2 is correct for m=0 and n=3, but the notation d_{mn} is introduced only for the even evaluation lemma; it would be clearer to state the admissible range of m and n explicitly before the display.","section":"§2.1, Lemma 2.4"}],"recommendation":"major_revision","confidential_remarks":"This is a promising manuscript with a clear program, but the proof of Lemma 3.1 is not valid as written and the main theorem depends on it. The dependence on [56] for the coercivity in Theorem 1.1 should be made explicit and, ideally, the proof should be included. I would not recommend rejection, because the gap may be repairable by a different positive-definiteness argument for the kernel or by a different coercivity strategy, but the manuscript cannot be accepted in its present form."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Hi [Colleague],\n\nYou should know about this paper: it proposes a genuinely new coupling between 3D Darcy and 1D Poiseuille models for blood perfusion, derives a closed 1D integrodifferential equation for the interior pressure, and claims a rigorous ε^{1/2}|log ε| convergence rate to a 1D Green's function model. The geometry (spheroidal tip, effective centerline) is handled with real care, and the residual estimates in sections 4-5 are internally consistent as far as I checked.\n\nThe problem is in the proof of Lemma 3.1, where the near-positivity of the kernel is established. The stress-test note is right: equation (3.17) claims an equality that doesn't hold. The bound from (3.16) gives 1/|R| ≤ 1/(ε a*), so the sum of two such terms is 1/(ε a*(τ)) + 1/(ε a*(t)), not 1/(ε² a*(τ)a*(t)). The square root of the diagonal product is η/(4π ε)(a*a*)^{-1/2}, which is much smaller than the claimed bound for small ε. So the inequality (3.17) is too weak to imply positive definiteness of I_{ε,0}. Since I_{ε,0} ≥ 0 is the only source of coercivity in the bilinear form B, the a priori bounds of Theorem 1.2 and the error estimate of Theorem 1.3 lose their stated support. This is a load-bearing gap, not a cosmetic one.\n\nThere are also two secondary concerns. First, Theorem 1.1's well-posedness is deferred to an unpublished companion [56], which strains self-containedness. Second, the numerics show only the 1D solution, never the claimed convergence against the 3D-1D system, so the rate is untested.\n\nThat said, I don't think the underlying idea is wrong. The error is a mis-estimate that might be repaired by a different positive-definiteness argument for the kernel, e.g., by working directly with the θ-integral representation. The rest of the paper is careful and the authors are clearly serious.\n\nMy recommendation: send it to a referee, but flag Lemma 3.1 as a must-fix. If the positivity can be recovered, the paper is a solid contribution. As written, it's not acceptable.\n\nBest,\n[You]","headline":"Promising new 1D perfusion model but a load-bearing gap in the kernel positivity proof means the main theorem lacks support as written.","tokens_in":43777,"tokens_out":5855,"would_cite":false,"duration_ms":55023,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B25","35J05","35R09","76S05","92C35"],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proves that a 1D slender-body equation for blood pressure around a thin vessel is the ε→0 limit of a coupled 3D-1D Darcy–Poiseuille model, with an explicit error of order ε^{1/2}|log ε|.","keywords":["blood perfusion","Darcy–Poiseuille","slender body approximation","integrodifferential equation","Green's function","free-end vessel","3D-1D coupling","asymptotic convergence"],"falsifier":"Compute the quadratic form in Lemma 3.1 for a radius function that vanishes linearly at the tip, for example a(s)=1−s near s=1, which violates (1.4): if the negative lower bound grows worse than C $ε^{{1/2}}$|log ε|^{1/2}, the coercivity step (3.78), and with it Theorem 1.2 and the rate in Theorem 1.3, would fail. Equivalently, solve (1.17) numerically for such a blunt tip and check whether ∥p_{SB}∥_{L∞} grows faster than $ε^{{-1/2}}$.","tokens_in":42351,"feed_emoji":"🩸","tokens_out":6249,"duration_ms":64798,"temperature":0.7,"pith_summary":"The paper tries to prove that a much simpler one-dimensional model captures the behavior of a coupled 3D-1D blood-perfusion model as the vessel radius ε tends to zero. The 1D model expresses the exterior pressure through an explicit half-space Green's function, while the interior pressure solves a new 1D integrodifferential equation whose kernel encodes the vessel geometry. The main work is a detailed analysis of that 1D equation: existence, uniqueness, weighted a priori bounds with explicit ε-dependence, and kernel estimates establishing near-positivity and near-antisymmetry. If the proof is right, one may replace the 3D-1D system by the cheaper 1D equation for curved vessels, with a quantified error, and the bounds feed directly into the companion paper's convergence to the full 3D-3D Darcy–Stokes system.","feed_headline":"1D blood-vessel model proven to match 3D-1D at rate ε^{1/2}|log ε|","feed_subtitle":"A new 1D integrodifferential equation tracks the coupled 3D-1D Darcy–Poiseuille system to within ε^{1/2}|log ε|.","key_machinery":"The load-bearing object is the 1D integrodifferential equation (1.17) for the interior pressure p_{SB}(s), whose kernel K_ε(s,t) is the θ-average of the half-space Neumann Green's function evaluated on the vessel surface. The proof proceeds through three kernel lemmas: near-positivity of the integral operator (Lemma 3.1), near-antisymmetry of its derivatives (Lemma 3.2), and weighted integration by parts (Lemma 3.3). These supply coercivity of the bilinear form B in (3.76) and the weighted a priori bounds (1.20)-(1.21); then the error estimate Theorem 1.3 follows from residual bounds for the normal derivative of q_{SB} (Lemma 4.1) and an energy identity on the difference between the two solutions.","core_discovery":"The central claim is that for a curved vessel whose radius εa(s) satisfies the admissible-radius conditions (1.2)-(1.4), the solution (q_{SB}, p_{SB}) of the 1D slender-body model (1.16)-(1.17) is within C $ε^{{1/2}}$|log ε| |p_0| of the solution (q, p) of the 3D-1D system (1.9a)-(1.9d), measured in the mixed norms of (1.22). This rests on Theorem 1.2, which states that the degenerate 1D integrodifferential equation has a unique solution satisfying ∥p_{SB}∥_{L∞} ≤ C $ε^{{-1/2}}$|p_0| and ∥$a^{{3/2}}$($a^{4}$ p_{SB,s})_{ss}∥_{$L^{2}$} ≤ C |log ε| |p_0|. The paper argues these estimates make precise the sense in which the explicit 1D Green's function model is the asymptotic reduction of the 3D-1D coupling as ε→0.","pith_inferences":["One consequence the paper leaves implicit is that the ε^{1/2}|log ε| rate is largely set by the spheroidal free-end treatment; a blunt or linearly tapered tip would likely destroy the near-positivity of Lemma 3.1 and the coercivity argument.","If the numerically suggested uniform-in-ε L∞ bound for p_{SB} can be proved, it would improve the interior-pressure estimate in Theorem 1.2, though it would probably not change the leading error rate of Theorem 1.3.","Because the 1D model is closed-form and geometry encoded in a kernel, it is a natural surrogate for shape-optimization and inverse problems over vascular arrangements, an application the paper mentions but does not develop."],"forward_implications":["The 1D equation (1.17) can be used in place of the coupled 3D-1D system for curved vessels with spheroidally tipped, non-constant radius, with a guaranteed approximation error O(ε^{1/2}|log ε|) in the norms of Theorem 1.3.","The weighted bounds of Theorem 1.2 serve as the 1D regularity input that the companion Part II uses to compare both reduced models directly with the full 3D-3D Darcy–Stokes system.","The analysis extends to an infinite-slab geometry with prescribed pressures at both ends, as described in Remark 1, at the cost of replacing the single-reflection Green's function by an infinite sum of reflections.","The numerical experiments in Section 1.5 suggest that the ε^{-1/2} bound on ∥p_{SB}∥_{L∞} itself may not be sharp, although the full equation with the integral term does not satisfy a maximum principle near a self-approaching tip."],"supporting_citations":[{"why":"Supplies the slender-body error-analysis framework and the base integral estimates that Section 2 adapts to the present setting.","marker":"[50]"},{"why":"Provides the free-end vessel geometry, the admissible radius conditions (1.2)-(1.4), and most of the integral lemmas reused in the kernel estimates.","marker":"[51]"},{"why":"Companion paper whose Section 4 supplies the coercivity of the 3D-1D bilinear form used in Theorem 1.1, and which consumes the Theorem 1.2 bounds for the Part II convergence result.","marker":"[56]"},{"why":"Gives the half-space Neumann Green's function whose explicit form defines the 1D slender-body operator in (1.15).","marker":"[16]"},{"why":"Provides the positive-definite kernel criterion used to show the main integral I_{\\epsilon,0} in Lemma 3.1 is non-negative.","marker":"[20]"},{"why":"Defines the homogeneous Sobolev space D^{1,2} and the functional framework for the well-posedness statement of Theorem 1.1.","marker":"[26]"}],"fun_headline_variants":["Blood vessel model hierarchy: 3D-1D reduces to 1D at ε^{1/2}|log ε|","1D slender-body model matches 3D-1D Darcy-Poiseuille to ε^{1/2}|log ε|","Proven: 1D Green's function model emerges from 3D-1D Darcy-Poiseuille","Convergence chain: 3D-1D to 1D blood flow with ε^{1/2}|log ε|","New 1D integrodifferential equation captures curved-vessel blood flow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The load-bearing premise is the spheroidal tip condition (1.4), which makes the kernel in Lemma 3.1 nearly positive and closes the coercivity estimate, with the well-posedness of the 3D-1D solution itself deferred to the companion paper.","fun_headline_variants_meta":{"raw":{"variants":["Blood vessel model hierarchy: 3D-1D reduces to 1D at ε^{1/2}|log ε|","1D slender-body model matches 3D-1D Darcy-Poiseuille to ε^{1/2}|log ε|","Proven: 1D Green's function model emerges from 3D-1D Darcy-Poiseuille","Convergence chain: 3D-1D to 1D blood flow with ε^{1/2}|log ε|","New 1D integrodifferential equation captures curved-vessel blood flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.001313,"raw_usage":{"total_tokens":5439,"prompt_tokens":1124,"completion_tokens":4315,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":740,"completion_tokens_details":{"reasoning_tokens":4172}},"tokens_in":740,"tokens_out":4315,"duration_ms":31016,"temperature":1.0,"reasoning_tokens":4172,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:28:44.873414+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Compute the quadratic form in Lemma 3.1 for a radius function that vanishes linearly at the tip, for example a(s)=1−s near s=1, which violates (1.4): if the negative lower bound grows worse than C $ε^{{1/2}}$|log ε|^{1/2}, the coercivity step (3.78), and with it Theorem 1.2 and the rate in Theorem 1.3, would fail. Equivalently, solve (1.17) numerically for such a blunt tip and check whether ∥p_{SB}∥_{L∞} grows faster than $ε^{{-1/2}}$.","supporting_citations":[{"cited_title":"Ohm and S","cited_arxiv_id":null,"evidence_quote":"Companion paper whose Section 4 supplies the coercivity of the 3D-1D bilinear form used in Theorem 1.1, and which consumes the Theorem 1.2 bounds for the Part II convergence result."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Gives the half-space Neumann Green's function whose explicit form defines the 1D slender-body operator in (1.15)."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the positive-definite kernel criterion used to show the main integral I_{\\epsilon,0} in Lemma 3.1 is non-negative."}],"review_version":1}