{"id":"8163db6a-2144-4683-b944-f3e9430f4aa3","arxiv_id":"2507.13653","paper_version":1,"verdict":"REJECT","confidence":"MODERATE","novelty_score":4.0,"correctness_risk":"high","formal_verification":"none","parameter_count":6,"one_line_summary":"A segmented formula gives capillary flow as a quartic elastic-compression law below 30 mmHg, a collapse-corrected law from 30 to 40 mmHg, and an exponential decay above 40 mmHg.","lead":"This paper proposes a three-phase mathematical formula describing how capillary blood flow drops as external pressure rises, splitting the response into gentle elastic compression, transitional elliptical collapse, and high-pressure closure. The authors intend it as a quantitative tool for dosing compression garments and tourniquets.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Under Table 1's own parameters, the low-pressure branch of Eq. (2) closes the capillary at 7.5 mmHg, so the claimed 0–30 mmHg elastic phase is unphysical and the central three-phase law is internally inconsistent.","rationale":"I read the paper as proposing Eq. (2) as a closed-form dose–response tool for capillary flow versus external pressure, with patient-specific parameters. For that central claim to hold, each branch must be physical across its stated interval and the phase boundaries should not produce unphysical jumps. The load-bearing condition is that the linear-elastic compression branch R(q)=r0(1−q r0/(E h)) is valid for 0≤q<30 mmHg under the paper's own baseline parameters. It is not: the vessel fully closes at 7.5 mmHg, and the quartic branch then increases without bound, so the first phase describes no sensible flow behavior. This is not a minor calibration detail; it undermines the entire three-phase picture. The high-pressure boundary failure compounds it: the left limit at 40− is about 0.00064 Q0 while the right limit is Q0, an upward jump of three orders of magnitude at the supposed collapse threshold. I weighed the paper's honest limitations in Sections 5–6, including unvalidated thresholds, unfitted alpha, and acknowledged low-pressure sensitivity. Those are real but secondary; they do not repair the internal contradiction. Independent support is also weak: no code or formal verification, and the error-bar 'validation' in Fig. 2 uses model-generated data rather than measurements. Thus the central claim is not merely unproven; it is contradicted by the paper's own equations under its own Table 1 parameters. My conclusion matches the Reader's weakest assumption and supports the same REJECT verdict, so no adjustment is needed.","tokens_in":10292,"tokens_out":7546,"duration_ms":88728,"concrete_test":"Evaluate Eq. (2) at the default parameter values using consistent units (E=75 mmHg, r0=10 μm, h=1 μm) at q=7.5, 30, 40− and 40 mmHg. If Q/Q0 equals 0 at 7.5 mmHg, 81 at q=30 mmHg, and at the 40 mmHg boundary the left limit is ≈6.4×10^-4 while the right value is 1, the three-branch model is internally inconsistent. This is a purely arithmetic check from the text and requires no external data.","verdict_should_be":"UNCHANGED","load_bearing_attack":"Eq. (2) is the central quantitative claim: for q<30, Q=Q0(1−q r0/(E h))^4. Using the paper's defaults (E=10 kPa, r0=0.01 mm, h=0.001 mm) and 1 mmHg≈133 Pa, the closure pressure is E h/r0 = (10^4 Pa·10^-6 m)/(10^-5 m)=10^3 Pa≈7.5 mmHg. Hence at 7.5 mmHg the radius reaches zero, and for any larger pressure the negative radius enters a fourth power that grows without bound: at q=30 mmHg the branch gives 81 Q0. This contradicts the model's own definition of 0–30 mmHg as a mild linear-elastic compression phase and the small-strain premise of §2.1. Section 6 flags that some parameter combinations drive the term toward zero 'prematurely,' but these are the paper's baseline values, not an edge case. The second boundary is equally discontinuous: approaching q=40 from below gives Q/Q0=(1−40/50)^4·0.4≈6.4×10^-4, while at q=40 the high-pressure branch Eq. (1) gives e^0=1, a thousand-fold upward jump instead of the 'closure-induced attenuation' the paper describes. These are internal arithmetic inconsistencies in the stated law under the stated parameters, not disagreements with external consensus; they make the piecewise function an invalid dose–response curve as written.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper proposes a three-phase piecewise model for capillary blood flow under external pressure: an elastic-compression regime for q < 30 mmHg, an elliptical-collapse transition for 30 ≤ q < 40 mmHg, and an exponential closure-decay regime for q ≥ 40 mmHg. The central quantitative claim is Eq. (2), which expresses normalized flow Q(q)/Q0 as (1 − q r0/(E h))^4, (1 − q/K_collapse)^4·β(q), and exp(−α(q−40)) in the three regimes. The authors perform sensitivity analyses for the parameters E, K, and α, present simulated error-bar analyses, and compare qualitative trends with published literature to argue that the model captures key determinants of pressure-induced flow reduction and can support dose–response assessment for compression therapies.","tokens_in":10481,"tokens_out":3547,"duration_ms":40209,"significance":"If the model were quantitatively valid, it would provide a closed-form, parameterized tool for predicting capillary perfusion under external pressure, with potential utility in tourniquet and compression-garment dosing and patient-specific tuning. The paper is also explicit about its mechanistic assumptions and provides sensitivity maps that could guide future experiments. However, the central law Eq. (2) is internally inconsistent under the paper's own default parameters: the low-pressure branch becomes unphysical at 7.5 mmHg, predicting growing flow instead of attenuation throughout most of the claimed 0–30 mmHg elastic phase, and the piecewise function has a large upward discontinuity at the 40 mmHg boundary. These are load-bearing arithmetic defects, not merely presentation issues. The validation is also circular because the error bars are generated from the model's own simulated outputs and then used to assert reliability, and the literature comparison is qualitative rather than a test of the specific functional forms. The conceptual three-phase framing may have heuristic value, but the quantitative claim as stated cannot be accepted.","major_comments":[{"comment":"With the paper's own Table 1 defaults (E = 10 kPa, r0 = 0.01 mm, h = 0.001 mm) and 1 mmHg ≈ 133 Pa, the compression term q r0/(E h) equals 1 at q ≈ 7.5 mmHg, so the modeled radius R(q) = r0(1 − q r0/(E h)) reaches zero and then becomes negative inside the claimed 0–30 mmHg elastic phase. For q > 7.5 mmHg the fourth power grows without bound; at q = 30 mmHg the branch would give Q/Q0 = (1 − 30·133·10^−5/(10^4·10^−6))^4 = 81, an increasing flow rather than the monotone attenuation the phase description requires. This contradicts the small-strain (<15%) premise of §2.1 and the statement in §6 that the term only 'approach[es] zero prematurely' for 'certain parameter combinations'; the baseline parameters are exactly one such combination.","section":"§3.2, Eq. (2), first branch"},{"comment":"The piecewise function is discontinuous in the wrong direction at q = 40 mmHg. Approaching from below with the transition-region form and the default values K_collapse = 50 mmHg and β = 0.4 gives Q/Q0 = (1 − 40/50)^4 · 0.4 = 0.00064, while the high-pressure branch at q = 40 gives exp(0) = 1. The flow would thus jump upward by a factor of about 1500 at the moment the vessel supposedly enters the 'closure-induced attenuation' phase, which is inconsistent with the model's own physiological narrative and with the simulated curves in Figures 1–2 that show a continuous decline across this boundary.","section":"§3.4–§3.5, Eq. (2), boundary at q = 40 mmHg"},{"comment":"The error-bar analysis is circular as a validation. The text states that the data are 'simulated data' used to illustrate model prediction variability; the error bars therefore quantify only the sampling variability of the model's own outputs (with n = 20) and cannot support the conclusions that the model has 'good repeatability,' 'high credibility,' or 'reliability' in the respective pressure regimes. Claims of agreement with physiology would require comparison against independent experimental measurements, which are not provided in the paper.","section":"§4, Fig. 2(b,d,f)"},{"comment":"Section 5 offers only qualitative comparisons with selected references, and several of those studies concern negative pressure, pathological states, or pressure ulcer formation rather than the positive-pressure 0–60 mmHg range of the model. Given that the model contains multiple free parameters (K_collapse, β, α, q_yield, and the 30 and 40 mmHg thresholds), a qualitative match to such trends does not test the specific functional forms in Eq. (2), so the claim that the model is 'validated' is overstated.","section":"§5, literature validation"}],"minor_comments":[{"comment":"The endothelial-yield modification Q(q) = Q0(1 − q/q_yield)^4 is introduced in the text but does not appear in the model summary Eq. (2); the relation between this alternative form and the primary low-pressure branch should be made explicit.","section":"§3.2"},{"comment":"β is listed as a 'fitted value' and α as an 'assumption,' but no fitting procedure or data source for these values is described; the reader cannot reproduce the parameter choices.","section":"Table 1"},{"comment":"The x-axis labels in panel (b) appear corrupted (e.g., '258.404932449.5868723372.1805583790.90053'), and in the text 'Figure a' should be 'Figure 2(a)'.","section":"Fig. 2(b)"},{"comment":"The equation for internal pressure P_int = (E h / r0)(r/r0 − 1) is dimensionally inconsistent as written unless the factor (1 − ν^2) is moved to the denominator; please clarify the notation and assumptions for Poisson's ratio.","section":"§2.2"},{"comment":"The reference to Tagliabue (2023) appears incomplete, and several in-text citations (e.g., 'Panula 2022', 'Negosanti et al. 2012') do not follow a consistent format; a careful reference cleanup is needed.","section":"References"}],"recommendation":"reject","confidential_remarks":"The paper's central quantitative law cannot be accepted because of the internal arithmetic inconsistencies under its own stated parameters. The authors do acknowledge a 'premature' approach to zero in the low-pressure regime in Section 6, but they treat it as a parameter-sensitivity issue rather than recognizing that it invalidates the baseline curves and the claimed 0–30 mmHg phase. The validation strategy is also circular. If the authors were able to reformulate the low-pressure branch with a physically grounded compression limit and re-anchor the high-pressure branch to eliminate the discontinuity, the three-phase concept might warrant a fresh look, but that would be a substantial rewrite rather than a minor revision."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Here's the quick take: the paper assembles a tidy three-branch closed form for capillary flow under external pressure, and the prose is honest about several limitations, but the central equation contradicts itself under the paper's own baseline parameters. The low-pressure branch Q=Q0(1 - q r0/(Eh))^4 with Table 1 (E=10 kPa, r0=0.01 mm, h=0.001 mm) closes the vessel at about 7.5 mmHg—well inside the claimed 0–30 mmHg elastic phase—and beyond that the quartic term grows without bound, giving 81 Q0 at 30 mmHg. At the other boundary, the transition branch approaches roughly 6.4e-4 Q0 as q→40, but the high-pressure branch equals Q0 at q=40, a thousand-fold upward jump. These are arithmetic consequences of the stated equations, not subtle physiology.\n\nWhat's genuinely new is the packaging: a single piecewise law with an elastic phase, an elliptical-correction phase, and an exponential closure phase, intended as a dose-response map for compression therapy. Each ingredient is standard—Poiseuille, thin-wall elasticity, elliptical hydraulic resistance, exponential decay—but I don't know of another paper putting them together this compactly. The sensitivity analysis is clearly motivated, and the authors explicitly flag that thresholds are unvalidated, alpha is an assumption, and the low-pressure branch is parameter-sensitive. That honesty counts.\n\nBut the internal inconsistency is load-bearing, not cosmetic. Section 3.2's small-strain justification is violated by the same parameters used to produce Figures 1–2. Section 6 acknowledges that some parameter combinations drive the flow term to zero 'prematurely,' but those are the paper's own baselines. The validation is qualitative trend matching against literature; the error-bar analyses in Figure 2 use simulated data generated by the model itself, so they carry no independent weight. The q_yield modification in Section 6 never makes it back into Eq. (2).\n\nWho is this for? Someone building a simplified clinical compendium of compression-pressure effects might find the three-phase concept useful, but they cannot use Eq. (2) as it stands. With a corrected low-pressure branch and continuous stitching between branches, plus real data, this could become a reasonable engineering tool. As submitted, the central law doesn't describe a sensible pressure-flow curve under its own parameters.\n\nI'd send it to a referee anyway—the flaw is checkable, the topic is clinically relevant, and a reviewer could point the authors to a concrete fix. But I would expect a reject-and-resubmit or major revision, not acceptance. And I wouldn't cite the current version.","headline":"The three-branch formula is internally inconsistent under the paper's own baseline parameters: the elastic branch closes the vessel at 7.5 mmHg and the high-pressure branch jumps upward at 40 mmHg, so the central dose-response law as written is not a valid model.","tokens_in":11196,"tokens_out":2986,"would_cite":false,"duration_ms":32494,"reading_group":"maybe","serious_thinker":"no","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":[],"pacs":[],"model":"deepseek-v4-flash","headline":"This paper proposes that capillary blood flow under external pressure follows a three-branch law: elastic compression below 30 mmHg, elliptical collapse from 30 to 40 mmHg, and exponential decay above 40 mmHg.","keywords":["capillary blood flow","external pressure","vascular collapse","three-phase model","hemodynamics","microcirculation","pressure therapy","dose-response"],"falsifier":"Measure capillary flow under external pressure from 0 to 60 mmHg in a vessel with the Table 1 parameters. The model's first branch predicts $Q=0$ at $q = Eh/r_0 \\approx 7.5$ mmHg; observing measurable flow at, say, 20 mmHg, or a decline that is not quartic, would refute the first branch and with it the three-phase curve.","tokens_in":9889,"feed_emoji":"🩸","tokens_out":11592,"duration_ms":105402,"temperature":0.7,"pith_summary":"This paper proposes a quantitative, piecewise law for how external pressure reduces capillary blood flow. It claims that the pressure–flow curve splits into three regimes: elastic compression below 30 mmHg, elliptical collapse between 30 and 40 mmHg, and exponential decay above 40 mmHg, each with an explicit formula. The payoff is a dose–response tool for pressure-based treatments such as tourniquets and compression garments, with patient-specific parameters for vessel stiffness, collapse threshold, and decay rate. The paper argues that the model captures the key determinants and matches literature trends in each pressure range.","feed_headline":"Capillary flow obeys a three-phase pressure law","feed_subtitle":"Below 30 mmHg perfusion holds; above 40 mmHg flow decays exponentially, setting a safety window for compression therapy.","key_machinery":"The argument rides on a piecewise application of Poiseuille's law, $Q \\propto R^4$, where the effective radius $R(q)$ is replaced by three deformation laws. In the low-pressure branch, $R = r_0(1 - q r_0/(E h))$ from linear elasticity; in the transition branch, the quartic collapse factor $(1 - q/K_{\\text{collapse}})^4$ is multiplied by an elliptical resistance correction $\\beta = (a/b)^2/(1 + (a/b)^2)$; in the high-pressure branch, the lumen closure is modeled as an exponential decay $\\exp(-\\alpha(q-40))$. These three mechanisms carry the entire modeling claim: there is no other fitting or computational machinery in the paper beyond these algebraic expressions and the sensitivity scans around them.","core_discovery":"On the paper's own terms, the central result is Eq. (2): normalized flow $Q(q)/Q_0$ equals $(1 - q r_0/(E h))^4$ for $q < 30$ mmHg, $(1 - q/K_{\text{collapse}})^4 \\beta(q)$ for $30 \\leq q < 40$ mmHg, and $\\exp(-\\alpha(q-40))$ for $q \\geq 40$ mmHg. Here $Q_0$ is the undisturbed flow, $r_0$ the initial radius, $E$ the wall elastic modulus, $h$ the wall thickness, $K_{\\text{collapse}}$ a stiffness parameter, $\\beta$ a function of the elliptical aspect ratio, and $\\alpha$ a decay coefficient. The three branches are claimed to correspond to distinct physical deformation mechanisms: linear elastic compression, shell-buckling into an ellipse, and lumen closure. The paper presents sensitivity analyses and literature comparisons as evidence that this three-phase structure describes capillary flow under external pressure.","pith_inferences":["From the paper's own baseline parameters, the elastic branch reaches zero flow at $q = Eh/r_0 \\approx 7.5$ mmHg, far below the claimed 30 mmHg boundary; if this linear-compression law is taken literally, the first branch must be replaced by a nonlinear elasticity before the three-phase picture can be upheld. This is an inference about internal consistency, not a paper claim.","The threshold structure suggests a direct experiment: if 30 and 40 mmHg are real physiological boundaries, compression garments below 30 mmHg should preserve perfusion while those above 40 mmHg should show a sharp exponential ischemia; laser speckle contrast imaging or OCT could test this without new modeling.","Applying the elastic branch to suction by substituting $q \\rightarrow -q$ predicts flow enhancement at moderate negative pressure; the paper only notes a tentative link to existing observations, so this is a testable extension, not a paper claim.","The transition branch uses $\\beta$ as a fitted constant; a full shell-buckling model that computes the aspect ratio $a/b$ from the applied pressure would remove that fit and make the branch fully predictive."],"forward_implications":["Below 30 mmHg, the model predicts that flow is governed by the quartic of a linearly compressed radius, so higher elastic modulus $E$ or thicker wall $h$ preserves perfusion under the same pressure.","In the 30–40 mmHg band, flow drops steeply and the reduction is controlled jointly by the collapse stiffness $K$ and the elliptical aspect-ratio factor $\\beta$, so vessels that are already elliptical (fibrosis, edema) lose perfusion faster.","Above 40 mmHg, flow decays exponentially with rate $\\alpha$, meaning residual perfusion at high pressure is set by $\\alpha$ and not by wall elasticity.","The piecewise structure implies a safety threshold: compression therapy should stay below roughly 30–40 mmHg to avoid collapse, and brief excursions above 40 mmHg could be used deliberately for hemostasis.","Because $E$, $K$, $\\beta$, and $\\alpha$ are patient-specific, the model yields an individual pressure–flow dose–response curve for tourniquets and compression garments."],"supporting_citations":[{"why":"Supplies the Poiseuille equation $Q = \\pi \\Delta P R^4/(8\\eta L)$ on which the whole model is built.","marker":"Doh et al. 2021"},{"why":"Supplies the elastic shell buckling theory behind the elliptical-collapse geometric relation in the transition branch.","marker":"Han et al. 2004"},{"why":"Provides the elliptical correction factor $\\beta$ for hydraulic resistance in elliptical cross-sections.","marker":"Garg et al. 2024"},{"why":"Provides the low-pressure literature match showing a positive correlation between vessel diameter and blood flow.","marker":"Abdullah 2022"},{"why":"Provides clinical evidence that external pressure induces vascular collapse, anchoring the 30–40 mmHg transition zone.","marker":"Langemo et al. 2008"},{"why":"Provides multiscale modeling evidence that pressure-driven vessel collapse underlies pressure ulcer formation, supporting the high-pressure decay phase.","marker":"Dharmangadan Sree 2019"}],"fun_headline_variants":["Capillary flow: three phases, two critical pressures","Capillary flow obeys a 30–40 mmHg phase switch","Elastic, elliptical, closed: capillary flow's pressure phases","Nonlinear capillary law: hold, buckle, close under pressure","Three-phase pressure law governs capillary blood flow"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The three-phase law rests on the premise that, below 30 mmHg, the capillary radius shrinks linearly as $R(q)=r_0(1 - q r_0/(E h))$; with the paper's own baseline values ($E=10$ kPa, $r_0=0.01$ mm, $h=0.001$ mm), this radius reaches zero at about 7.5 mmHg, well inside the claimed elastic phase.","fun_headline_variants_meta":{"raw":{"variants":["Capillary flow: three phases, two critical pressures","Capillary flow obeys a 30–40 mmHg phase switch","Elastic, elliptical, closed: capillary flow's pressure phases","Nonlinear capillary law: hold, buckle, close under pressure","Three-phase pressure law governs capillary blood flow"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000464,"raw_usage":{"total_tokens":2281,"prompt_tokens":874,"completion_tokens":1407,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":490,"completion_tokens_details":{"reasoning_tokens":1326}},"tokens_in":490,"tokens_out":1407,"duration_ms":14239,"temperature":1.0,"reasoning_tokens":1326,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-06T16:20:23.246272+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure capillary flow under external pressure from 0 to 60 mmHg in a vessel with the Table 1 parameters. The model's first branch predicts $Q=0$ at $q = Eh/r_0 \\approx 7.5$ mmHg; observing measurable flow at, say, 20 mmHg, or a decline that is not quartic, would refute the first branch and with it the three-phase curve.","supporting_citations":[{"cited_title":"Archives of Hand and Microsurgery 26(4), 285–292 (2021) Dharmangadan Sree, V.: Multiscale and multiphysics modeling of pressure driven ischemia and ulcer formation in the skin","cited_arxiv_id":null,"evidence_quote":"Supplies the Poiseuille equation $Q = \\pi \\Delta P R^4/(8\\eta L)$ on which the whole model is built."},{"cited_title":"Composites Part B: Engineering 35(6-8), 591–598 (2004)","cited_arxiv_id":null,"evidence_quote":"Supplies the elastic shell buckling theory behind the elliptical-collapse geometric relation in the transition branch."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Provides the elliptical correction factor $\\beta$ for hydraulic resistance in elliptical cross-sections."},{"cited_title":"International Journal of Nonlinear Analysis and Applications 13(1), 1341–1350 (2022)","cited_arxiv_id":null,"evidence_quote":"Provides the low-pressure literature match showing a positive correlation between vessel diameter and blood flow."},{"cited_title":"Advances in Skin & Wound Care 21(6), 282–292 (2008)","cited_arxiv_id":null,"evidence_quote":"Provides clinical evidence that external pressure induces vascular collapse, anchoring the 30–40 mmHg transition zone."}],"review_version":1}