{"id":"c455771a-18bb-402b-a6e9-d1776d88d926","arxiv_id":"2411.18703","paper_version":1,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":1,"one_line_summary":"Simulations show that red blood cell aggregation increases wall shear stress fluctuations in microvessels, a possible mechanical cause of endothelial damage in blood disorders.","lead":"Using computer simulations of blood flow in tiny vessels, this study shows that when red blood cells clump together, they create uneven flow near the vessel wall and cause strong, intermittent stress fluctuations. This mechanical stress may help explain why patients with conditions like COVID-19 and sickle cell disease suffer damage to the protective lining of their blood vessels.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Kagg=1 is imposed at the balance point between aggregation and viscous stress without calibration to any pathological condition; the central WSS-fluctuation result rests on this single untested parameter value.","rationale":"The reader identified the calibration of Kagg = 1 as the weakest assumption, and I agree that this is the load-bearing issue. The paper is otherwise carefully constructed: the CFL thickening under aggregation matches experimental observations [34,35], the f proportional to hematocrit scaling is internally consistent, and the event-count estimate independently supports the PSD peak interpretation. However, none of these checks establishes that Kagg = 1 represents aberrant aggregation in disease. Because Kagg = 1 is the exact balance between the aggregation force scale and the imposed wall shear stress, the system sits at a threshold; small changes in De or in local shear rate could suppress the near-wall clusters that produce the hallmark WSS fluctuations. A Kagg sweep or an experimentally calibrated De is therefore necessary before the quantitative claim can be taken as established. The reader's conditional verdict remains appropriate; no reason to reject appears, but the current evidence is not sufficient to fully accept the mechanism as quantitatively linked to pathological aggregation.","tokens_in":16957,"tokens_out":6662,"duration_ms":69203,"concrete_test":"Repeat the straight-tube healthy and SCD simulations at Kagg = 0.25, 0.5, 1, and 2 with otherwise identical setups, and compare the high-WSS tail (e.g., P(tau_w > 2)) and the PSD peak amplitude near f = 0.04. If the effect is absent or much weaker below Kagg = 1, or if the peak shifts non-monotonically, the claimed mechanism is an artifact of the chosen balance point. Ideally, set De from fibrinogen-dependent aggregation energy measurements (e.g., Dasanna et al. 2022 or Deng et al. 2020) and rerun the key cases to check whether pathological conditions map to Kagg near, below, or above unity.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The central claim requires that the simulated aggregation strength actually corresponds to aberrant, pathologically elevated RBC aggregation. The model fixes r0 = 0.49 µm and beta = 3.84 µm^-1 from Zhang et al. [46], but the energy scale De is not calibrated to fibrinogen concentration or to any measured pathological aggregation energy. Instead, De is chosen through the dimensionless condition Kagg = De*beta/(eta*gamma_w) = 1. In Poiseuille flow, the local shear stress at one cell radius from the wall is roughly 0.75 times the wall shear stress, while the maximum attractive traction in the Morse potential is De*beta/2 = 0.5*eta*gamma_w at Kagg = 1. Thus near-wall clusters are operating at or beyond the adhesion threshold, making the predicted near-wall aggregates and the resulting WSS peak sensitive to the exact choice of Kagg. Since every simulation compares only Kagg = 0 with Kagg = 1 and never sweeps intermediate values, the existence, magnitude, and frequency of the f = 0.04 peak are not robustly tied to a pathological aggregation level.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This manuscript presents immersed-boundary simulations of deformable red blood cell capsules in straight and curved cylindrical tubes, with a Morse potential describing cell-cell aggregation. It compares Kagg=0 and Kagg=1 for suspensions of normal RBCs, binary suspensions containing 10% stiff sickle cells, and variations in capillary number and hematocrit. The paper reports that aggregation thickens but destabilizes the cell-free layer, produces large wall shear stress fluctuations, and yields a distinct power spectral density peak near f = 0.04, an effect that is amplified by stiff sickle cells and by vessel curvature. The authors propose that this aggregation-induced wall shear stress environment explains glycocalyx disruption and endothelial damage in disorders with aberrant erythrocyte aggregation.","tokens_in":17168,"tokens_out":5144,"duration_ms":48187,"significance":"If the quantitative claims are robust, the paper provides a plausible mechanical pathway linking pathological RBC aggregation to endothelial injury, complementing recent microfluidic experiments by Druzak et al. and extending the authors' prior work on margination of aberrant RBCs. The model is physically detailed, using the Skalak membrane law, Canham-Helfrich bending energy, and an immersed boundary method, and the paper explicitly tests the effect of hematocrit and capillary number on the aggregation-driven phenomena. The scaling argument f ~ phi*gamma_dot and the observed shift of the PSD peak from f = 0.04 at phi = 0.2 to f = 0.06 at phi = 0.3 is a useful internal consistency check. The main caveats are that the aggregation strength is not calibrated to a pathological condition and that the headline spectral peak and high-stress tail lack statistical uncertainty estimates; these limit the strength of the quantitative conclusions but not the value of the qualitative findings.","major_comments":[{"comment":"The central comparison is between Kagg = 0 and Kagg = 1, but Kagg = 1 is set by the dimensionless balance condition Kagg = De*beta/(eta*gamma_dot) = 1 rather than calibrated to a fibrinogen concentration or to an aggregation energy measured under pathological conditions. Since the near-wall clusters and the f = 0.04 PSD peak arise only at Kagg = 1, the magnitude and even the existence of the predicted wall shear stress fluctuations are not robustly tied to 'aberrant' aggregation. A sweep over Kagg (e.g., 0.5, 1, 2) or a calibration to measured aggregation energies such as those in refs. [46] or [50] is needed to support the quantitative claims.","section":"Section II, Eqs. (1)-(2), and the definition of Kagg"},{"comment":"The high-stress tail in Fig. 3(D) and the PSD peak in Fig. 3(E) are computed from what appears to be a single simulation. The text says the spectrum is 'averaged over many points' but does not state the number of independent samples, provide error bars, or report ensemble averaging over independent initial conditions. The accompanying consistency check (18 events with WSS > 2 over 500 time units) is a rough count, not an estimate of spectral uncertainty. Given that the abstract advertises a distinct peak at f = 0.04, the paper should report confidence intervals or replicate simulations to demonstrate that the peak is not a statistical fluctuation.","section":"Section III.A, Figs. 3(D) and 3(E)"},{"comment":"The text estimates f_w ~ phi*gamma_dot and then interprets the PSD rise below f ~ 0.3 as the effect of individual cells moving past the sampling point. With phi = 0.2 and gamma_dot = 1, this estimate gives f ~ 0.2, yet the claimed aggregation peak is at f = 0.04, a factor of five lower. The paper does not explain this discrepancy; if the peak is attributed to large aggregated regions rather than single cells, the relationship between cluster size or wavelength and the scaling argument should be made explicit. The relative shift from 0.04 to 0.06 with phi = 0.2 to 0.3 is consistent with f ~ phi, but the absolute frequency of the peak remains unexplained.","section":"Section III.A, scaling argument for f_w"}],"minor_comments":[{"comment":"There is a typo in 'substntial' in the sentence about marginated aberrant cells playing a substantial role.","section":"Section I.A"},{"comment":"The two bending moduli in the energy expression are both written as K_B, and the text also says 'KB and KB'; the Gaussian curvature term should have its own modulus, e.g., K_G, to avoid ambiguity.","section":"Section II, membrane energy equation"},{"comment":"The statement that 'the wall shear stress is dimensionless using gamma_dot ~ 100 s^-1' is dimensionally incomplete; presumably tau_w is normalized by eta*gamma_dot, and this should be stated.","section":"Figure 3(D) caption"},{"comment":"There is a typo 'peak frequenct', and the claim that at Ca = 1 the aggregation effect on wall shear stress 'remains limited' is made by visual inspection of PDFs; a quantitative measure such as the variance or tail probability would be more convincing.","section":"Section III.C"},{"comment":"The same Morse potential is applied to sickle-sickle and sickle-normal pairs; the authors acknowledge this simplification, but a brief discussion of how a stronger sickle-specific adhesion would affect the margination and wall stress conclusions would help the reader assess the robustness of the SCD-specific claims.","section":"Section III.B"},{"comment":"The statement that a verification run at Rep = 0.05 produced no changes in conclusions is not documented; please show the comparison or specify where it appears in the Supplementary Materials.","section":"Section II"}],"recommendation":"major_revision","confidential_remarks":"The paper fits the journal's scope as a fluid dynamics study with direct biomedical relevance. The main risk is the uncalibrated aggregation strength and the lack of statistical error analysis for the central spectral claims; both are addressable with modest additional simulations. I do not see a circularity problem: the results are emergent outputs of the simulations, and the scaling comparison is an independent consistency check rather than a fit."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"Colleague,\n\nThis paper uses direct cell-level simulations to argue that pathological RBC aggregation makes the cell-free layer thicker but noisier, and that aggregated clusters brushing the wall create a distinct peak in wall shear stress fluctuations at f≈0.04, amplified by sickle cells and vessel curvature. The core qualitative claim holds up: aggregation is a plausible mechanical pathway from elevated fibrinogen to glycocalyx damage, and it aligns with Druzak et al.'s microfluidic experiments. The genuinely new pieces are the WSS-fluctuation analysis, the enhancement of sickle-cell margination by aggregation, and the curved-vessel results; the existing aggregation literature focuses on rheology and CFL thickness, not the wall stress environment.\n\nWhat they do well: the model is built on a standard Morse-potential aggregation description, the Kagg dimensionless group is a sensible way to compare aggregation to viscous stress, and they run a verification at Rep = 0.05. The scaling argument f ~ ϕγ̇ is neat, and the internal consistency check (18 high-stress events over 500 time units gives f ≈ 0.036, close to the PSD peak at 0.04) shows the peak is a real event signature, not numerical noise.\n\nThe soft spots are real but not fatal. Kagg = 1 is set by a balance argument, not calibrated to any measured fibrinogen or dextran level; the stress-test note is right that near-wall clusters sit near the adhesion threshold, so the amplitude and frequency of the PSD peak could move if the true pathological Kagg is 0.5 or 2. A sweep over Kagg would have made the central claim much more robust, and its absence is the main quantitative weakness. Relatedly, the stress PDFs and PSDs come from single simulations without ensemble averaging or error bars, so the high-stress tail is not statistically characterized. Both issues are fixable and don't undermine the qualitative mechanism.\n\nWho it's for: people working on blood rheology, mechanobiology of the endothelium, and sickle cell or COVID vasculopathy. It deserves a serious referee; I'd send it to peer review and ask for a Kagg sweep plus ensemble statistics in revision rather than desk-rejecting it. I would not cite it in my own work yet—I'd wait for the calibrated version—but it's a worthwhile contribution to the discussion.\n\nBest","headline":"A careful simulation study that gives the aggregation-to-WSS-fluctuation mechanism a plausible physical footing, though the load-bearing parameter Kagg=1 is uncalibrated and the statistics are single-run.","tokens_in":17691,"tokens_out":2581,"would_cite":false,"duration_ms":24311,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["76Z05"],"pacs":[],"model":"deepseek-v4-flash","headline":"Simulations show that red blood cell aggregation alone can create damaging wall shear stress fluctuations in small blood vessels, a mechanical pathway for vascular injury.","keywords":["red blood cell aggregation","cell-free layer","wall shear stress fluctuations","microcirculation","sickle cell disease","margination","immersed boundary simulation","glycocalyx"],"falsifier":"Measure wall shear stress fluctuations in a roughly 32 micrometer microfluidic tube at 20 percent hematocrit and a wall shear rate near 100 per second, with dextran- or fibrinogen-induced aggregation tuned to Kagg approximately 1: the power spectral density should show the predicted peak at f approximately 0.04, about 4 per second dimensional, and an enhanced high-stress tail; if no such peak appears while aggregation is present, the proposed cluster-passage mechanism is not supported.","tokens_in":16744,"feed_emoji":"🩸","tokens_out":5584,"duration_ms":142562,"temperature":0.7,"pith_summary":"This paper uses cell-level simulations of deformable red blood cells in small tubes to establish a mechanical route from aberrant red blood cell aggregation to vascular damage. It shows that when aggregating forces and flow stresses are comparable, red blood cells form clusters and rouleaux that make the cell-free layer near the vessel wall thicker on average but far more variable, with cells sometimes approaching the wall closely. Those near-wall clusters produce large, intermittent spikes in wall shear stress, visible as a distinct spectral peak at a frequency near f=0.04 in the shear-rate-scaled units used here. The paper further shows that stiff sickle cells, which marginate toward the wall, and curved vessel geometry each amplify these shear-stress fluctuations, connecting the simulations to clinical observations of glycocalyx disruption and endothelial inflammation.","feed_headline":"Aggregated red blood cells trigger wall-stress spikes in microvessels","feed_subtitle":"Simulations tie abnormal red blood cell clumping to shear stress spikes on vessel walls, amplified by sickle cells.","key_machinery":"The load-bearing object is a phenomenological Morse potential between surface nodes of neighboring RBCs, with interaction parameters r0=0.49 micrometers and beta=3.84 per micrometer, and a dimensionless aggregation number Kagg=De*beta/(eta*gamma_w) that sets the balance between cell-cell attraction and viscous drag. At Kagg=1 the attraction roughly balances flow-induced separation, and this is the regime studied. The cells themselves are deformable fluid-filled capsules with biconcave resting shapes, while sickle cells are stiffer, smaller, curved-prolate capsules; the suspensions are solved with an immersed boundary method in pressure-driven straight and curved tubes. The analysis then rests on two diagnostics: the local cell-free layer thickness measured from the center of mass of the nearest cell, and wall shear stress sampled at fixed wall points, whose fluctuations are characterized by probability densities and power spectral densities.","core_discovery":"The central claim is that RBC aggregation alone, even without biochemical signaling, can create a damaging hemodynamic environment in the microcirculation. In the simulations, setting the aggregation strength to Kagg=1, where intercellular attraction balances the viscous stress of the roughly 100 per second flow, transforms a fairly uniform suspension into a heterogeneous one: rouleaux and clusters form, the mean cell-free layer thickness increases, and yet the probability of a very thin cell-free layer and of cells nearly touching the wall also rises. Passing aggregates generate wall shear stress excursions whose high-stress tail is orders of magnitude more probable than without aggregation, and the power spectral density develops a distinct peak near f=0.04, corresponding to clusters travelling close to the wall; the frequency of these events scales with hematocrit times shear rate. Adding 10 percent stiff sickle cells deepens margination, so the sickled cells sit closer to the wall and produce an additional higher-frequency peak near f=0.07. Curved vessels accentuate clustering and increase wall shear stress fluctuations, especially on the outer curve.","pith_inferences":["If the f approximately 0.04 peak is confirmed experimentally, wall-shear-stress spectral analysis could become a non-invasive indicator of pathological aggregation and endothelial risk in microfluidic or in vivo settings.","The model implicitly predicts that interventions that reduce aggregation strength should suppress the low-frequency wall stress peak before the mean cell-free layer thickness changes much, because the peak is driven by near-wall cluster passage rather than by average cell-free layer thickness.","The Kagg=1 choice is a single operating point; real aberrant aggregation varies with fibrinogen level and shear rate, so the quantitative threshold for endothelial damage likely depends on both, and mapping Kagg across disease states would test clinical relevance.","Because the same Morse potential is applied to sickle-normal and sickle-sickle pairs, the model isolates mechanical margination; adding specific sickle-cell adhesion could modulate the results and is a natural next step."],"forward_implications":["If aggregation is strong enough to reach the Kagg approximately 1 balance, microvessel walls experience intermittent high shear stress events even when the mean cell-free layer is thicker, so average cell-free layer measurements can mask a more dangerous fluctuation environment.","The reported peak at f approximately 0.04, about 4 per second at a 100 per second wall shear rate, gives a quantitative, testable signature of aggregation-induced wall stress fluctuations in straight vessels.","In sickle cell disease, aggregation and margination act synergistically: the same aggregation that pulls normal cells together pushes stiff sickled cells closer to the wall, raising the chance of wall contact and high stress.","Vessel curvature, as found in capillary networks, promotes RBC clustering and moves the largest shear stress fluctuations to the outer wall of the curve, indicating that geometry should be considered when assessing vascular injury risk.","Higher hematocrit, 30 percent versus 20 percent, shifts the wall shear stress fluctuation peak from f approximately 0.04 to f approximately 0.06, consistent with a linear scaling of passage frequency with volume fraction."],"supporting_citations":[{"why":"Provides the motivating experimental observation that high-fibrinogen blood flow directly damages the glycocalyx, which the simulations aim to explain mechanically.","marker":"[1]"},{"why":"Supplies the Morse potential parameters r0 and beta used to model RBC aggregation in the simulations.","marker":"[46]"},{"why":"Introduces the Morse-potential modeling approach for RBC aggregation that this paper adopts.","marker":"[44]"},{"why":"Prior computational model of marginated aberrant RBCs and wall stress fluctuations that this study extends to include aggregation.","marker":"[29]"},{"why":"Establishes the mechanism of margination through wall-induced migration and hydrodynamic collisions, underlying the cell-free layer and segregation behavior.","marker":"[24]"},{"why":"Demonstrates experimentally that stiff sickle RBCs alone biomechanically induce endothelial inflammation, providing the disease motivation.","marker":"[22]"},{"why":"Experimental observation that pathological aggregation increases cell-free layer thickness variance, a result the simulations reproduce.","marker":"[35]"},{"why":"Experimental evidence that dextran-induced aggregation thickens the arteriolar cell-free layer, supporting the simulated trends.","marker":"[34]"}],"fun_headline_variants":["RBC aggregation triggers wall stress spikes in microvessels","Clumped red cells raise microvessel wall shear stress","Sickle cells enlarge wall stress from RBC clumping","Model links RBC clumping to vessel wall stress bursts","Aberrant RBC aggregation heightens wall shear stress"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The simulations rely on a single adjustable attraction strength between red blood cells, chosen so that aggregation and flow forces are balanced; if the real pathological stickiness is weaker than this choice, the predicted wall-stress fluctuations and their frequency peak would be smaller or absent.","fun_headline_variants_meta":{"raw":{"variants":["RBC aggregation triggers wall stress spikes in microvessels","Clumped red cells raise microvessel wall shear stress","Sickle cells enlarge wall stress from RBC clumping","Model links RBC clumping to vessel wall stress bursts","Aberrant RBC aggregation heightens wall shear stress"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000291,"raw_usage":{"total_tokens":1726,"prompt_tokens":997,"completion_tokens":729,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":613,"completion_tokens_details":{"reasoning_tokens":652}},"tokens_in":613,"tokens_out":729,"duration_ms":18020,"temperature":1.0,"reasoning_tokens":652,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-12T10:57:14.447840+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure wall shear stress fluctuations in a roughly 32 micrometer microfluidic tube at 20 percent hematocrit and a wall shear rate near 100 per second, with dextran- or fibrinogen-induced aggregation tuned to Kagg approximately 1: the power spectral density should show the predicted peak at f approximately 0.04, about 4 per second dimensional, and an enhanced high-stress tail; if no such peak appears while aggregation is present, the proposed cluster-passage mechanism is not supported.","supporting_citations":[{"cited_title":"Druzak, E","cited_arxiv_id":null,"evidence_quote":"Provides the motivating experimental observation that high-fibrinogen blood flow directly damages the glycocalyx, which the simulations aim to explain mechanically."},{"cited_title":"Zhang, P","cited_arxiv_id":null,"evidence_quote":"Supplies the Morse potential parameters r0 and beta used to model RBC aggregation in the simulations."},{"cited_title":"Liu and W","cited_arxiv_id":null,"evidence_quote":"Introduces the Morse-potential modeling approach for RBC aggregation that this paper adopts."},{"cited_title":"Cheng, C","cited_arxiv_id":null,"evidence_quote":"Prior computational model of marginated aberrant RBCs and wall stress fluctuations that this study extends to include aggregation."},{"cited_title":"Mechanism of margination in confined flows of blood and other multicomponent suspensions","cited_arxiv_id":"1208.0943","evidence_quote":"Establishes the mechanism of margination through wall-induced migration and hydrodynamic collisions, underlying the cell-free layer and segregation behavior."},{"cited_title":"Caruso, X","cited_arxiv_id":null,"evidence_quote":"Demonstrates experimentally that stiff sickle RBCs alone biomechanically induce endothelial inflammation, providing the disease motivation."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental observation that pathological aggregation increases cell-free layer thickness variance, a result the simulations reproduce."},{"cited_title":null,"cited_arxiv_id":null,"evidence_quote":"Experimental evidence that dextran-induced aggregation thickens the arteriolar cell-free layer, supporting the simulated trends."}],"review_version":1}