{"id":"8c9a395b-f44e-489d-b2b8-d8c616dc65d1","arxiv_id":"2501.04017","paper_version":1,"verdict":"CONDITIONAL","confidence":"HIGH","novelty_score":6.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":5,"one_line_summary":"In simulated sickle cell blood, stiff marginated sickle red blood cells preferentially enter low-flow branches at bifurcations, while healthy cells keep the usual high-flow preference, and this alters wall shear stress patterns.","lead":"Using computer simulations of blood flowing through branching vessels, this paper shows that stiff, sickled red blood cells that drift to the vessel wall tend to enter the slower-flow branch, opposite to healthy cells. The finding points to a physical mechanism, cell position at the bifurcation, that could explain where vessel-wall damage occurs in sickle cell disease.","discovery_kind":"new_application","skeptic_critique":{"model":"deepseek-v4-flash","headline":"The sickle cell model uses a rest volume ~20% of normal RBCs, far below reported SCD cell volumes; the anti-Zweifach-Fung effect may depend on this exaggerated margination, and no sensitivity test is reported.","rationale":"I read the paper in good faith. It is a mechanistic simulation study that quantitatively reproduces the classical Zweifach-Fung effect for healthy cells and then shows that a marginated stiff-cell population can invert the partitioning because the low-flow branch captures a disproportionate share of the CFL. That reasoning is internally consistent, and the healthy-cell validation against Pries et al. (1989) provides genuine support for the simulation framework. The downstream segregation and WSS statistics are plausible extensions, though they are secondary and would also benefit from error bars. The single most load-bearing assumption is the sickle cell model itself, because the anti-Zweifach-Fung effect is entirely a consequence of strong margination. The reader identified this as the weakest assumption, and I agree. I sharpen it: the stated rest volume of 20% of normal is not just a simplification but a parameter value far outside typical SCD RBC volumes. Since smaller rigid particles marginate more strongly, this choice likely exaggerates the CFL population and thus the headline effect. The paper provides no sensitivity analysis on this parameter, nor on the 10% sickle-cell fraction, so the reader cannot assess whether the effect persists for physiologically realistic cells. This does not invalidate the paper as a mechanistic proof-of-concept, but it fully justifies the original CONDITIONAL verdict. I therefore recommend no verdict change, but the condition should specifically require a volume/stiffness sensitivity study or direct comparison to experimentally measured SCD cell margination profiles.","tokens_in":13433,"tokens_out":7235,"duration_ms":79048,"concrete_test":"Repeat the symmetric-bifurcation partitioning simulations at ηQ = 0.5, 0.7, and 0.9 with the sickle cell rest volume increased from 20% to roughly 70–80% of the normal RBC volume, keeping the same crescent shape and Ca = 0.2, and run at least three independent replicates per condition to obtain error bars. If ηN,sickle remains significantly below ηQ across these flow ratios, the anti-Zweifach-Fung effect is robust to the volume assumption; if it approaches or crosses ηQ, the headline claim is an artifact of the exaggerated 20% volume and should be reframed as a conditional, parameter-dependent prediction.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The headline claim—that sickle RBCs preferentially enter the low-flow branch (anti-Zweifach-Fung effect)—is mechanistically traced to the sickle cells' position in the cell-free layer (CFL) interacting with the shifted separatrix (Fig. 2C and Discussion). That position is set by the sickle cell model: a stiff capsule with a curved prolate spheroidal rest shape and a volume approximately 20% that of the normal RBC model (Formulation). This volume choice is physiologically extreme: normal RBC volume is roughly 90 fL, and even the most dehydrated dense sickle cells are typically reported at 60–80% of normal volume, not 20%. Because smaller, stiffer particles marginate more strongly, this parameter likely inflates the CFL population of sickle cells and therefore the magnitude of the anti-Zweifach-Fung effect. The authors themselves note that real SCD populations are heterogeneous and that their single-class model is 'a substantial simplification,' yet they provide no sensitivity analysis in cell volume, stiffness, or sickle-cell fraction. Since the mechanism is essentially determined by the inlet cross-sectional distribution, a less extreme—and more realistic—sickle cell would marginate less, and ηN for sickle cells could move toward ηQ, weakening or eliminating the headline effect. The healthy-cell Zweifach-Fung validation and the separatrix explanation are internally consistent, but they do not establish that the result survives at physiologically faithful sickle cell parameters.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"This paper presents immersed-boundary simulations of binary red blood cell (RBC) suspensions in model vascular bifurcations, with one population of normal biconcave capsules and one population of stiff, smaller, curved prolate spheroidal capsules intended to represent sickle RBCs. The authors report that normal RBCs reproduce the classical Zweifach-Fung effect, preferentially entering the higher-flow branch, while the model sickle cells instead preferentially enter the lower-flow branch—an 'anti-Zweifach-Fung effect'—which they attribute to the marginated position of the stiff cells in the cell-free layer coupled with the flow-dependent shift of the separatrix. Downstream of the bifurcation, sickle cells accumulate along the outer branch walls, and the simulations show that the SCD suspension produces a higher probability of elevated wall shear stress events than a healthy suspension, especially on the outer side of high-velocity branches. An additional set of simulations in a geometrically asymmetric bifurcation shows that cells preferentially enter the larger-radius branch and that the anti-Zweifach-Fung effect persists.","tokens_in":13706,"tokens_out":3305,"duration_ms":35487,"significance":"If the central finding is robust, it is novel and potentially clinically relevant: it provides a mechanistic, purely hydrodynamic explanation for how a stiff, marginating RBC subpopulation could be routed to low-flow branches in the microcirculation and how this could locally amplify wall shear stress fluctuations, contributing to endothelial dysfunction in sickle cell disease. The paper's strengths include validation of the healthy-cell partitioning behavior against published experimental microcirculation data (Fig. 1B, Ref. 61), a clear mechanistic account of the partitioning via the separatrix and the cell-free layer, and the use of an established simulation method (the authors' prior Science Advances work, Ref. 18). The anti-Zweifach-Fung effect emerges from the simulations rather than being imposed as an input, which is a genuine model output. However, the sickle cell is represented by a single, strongly marginating parameter set, and the quantitative claims are not supported by sensitivity analysis or uncertainty quantification, so the breadth of the physiological conclusion currently outruns the evidence.","major_comments":[{"comment":"The sickle cell model uses a rest volume approximately 20% that of the normal RBC and a membrane shear modulus five times larger. This parameter choice is physiologically extreme: reported sickle RBC volumes are typically 60–80% of normal, not 20%. Because the entire anti-Zweifach-Fung mechanism rests on the near-wall (cell-free layer) location of sickle cells, and because smaller, stiffer particles marginate more strongly, this extreme parameter set likely inflates the reported effect. The authors themselves call the single-class representation 'a substantial simplification,' yet they provide no sensitivity analysis in cell volume, stiffness ratio, or aberrant-cell fraction. Without such a test (e.g., repeating the key partitioning simulations at 50% and 75% of normal volume, or at stiffness ratios of 2 and 3), the reader cannot tell whether the anti-Zweifach-Fung effect is a robust feature of diseased RBCs or an artifact of the chosen idealized parameter point. This is the load-bearing assumption for the paper's headline result, so the omission is a major gap.","section":"Formulation"},{"comment":"No uncertainty quantification or replicate statistics are provided for any of the central quantitative claims. Figures 2(B), 4(D–F), and 5(B) show smooth curves and probability densities from what appear to be single simulations per condition, with no error bars, confidence bands, or indication of run-to-run variability. Given that cell trajectories through a bifurcation are stochastic, the statements that ηN for sickle cells lies below ηQ for ηQ > 0.5 (and above for ηQ < 0.5), and that the SCD suspension raises the probability of high WSS events, cannot be assessed for statistical significance as presented. At a minimum, the authors should report the number of independent simulations or the length of the stationary time window, and provide bootstrap or replicate-based error estimates on ηN and on the WSS probability distributions. This is necessary to support the quantitative, and especially the inverted-branch-preference, claim.","section":"Results, Partitioning and Wall Shear Stress"},{"comment":"The asymmetric-bifurcation result that cells preferentially enter the larger-radius branch even at equal volumetric flow (ηQ = 0.5) is presented without any sensitivity to the degree of geometric asymmetry. Only one pair of daughter radii (14 µm and 11 µm) is reported, and the claim that 'the geometric asymmetry of the bifurcation increases the curvature of the separatrix, further amplifying the uneven distribution' is not backed by a systematic variation of the radius ratio or bifurcation angle. Because the separatrix curvature is the proposed explanation, the lack of a parameter sweep leaves the generality of this secondary conclusion unsupported, even if the symmetric-bifurcation mechanism is accepted.","section":"Results, Geometrically Asymmetric Bifurcation"}],"minor_comments":[{"comment":"The word 'modling' in the first paragraph of the Discussion should be 'modeling'; the same paragraph also contains a missing article ('this is an anti-Zweifach-Fung effect' reads awkwardly in context).","section":"Discussion"},{"comment":"The sentence 'The Chorin projection method is utilized to advance the velocity field u. This method involves solving an advection-diffusion equation...' repeats the earlier statement about the projection method; the redundancy should be removed.","section":"Formulation"},{"comment":"In Fig. 1(B), the experimental data from Ref. 61 are plotted without error bars; reproducing the experimental uncertainty would help the reader judge the agreement quantitatively.","section":"Results, Validation"},{"comment":"The claim that 'the average WSS at point a remains close to zero' should be supported by a numerical value, since the distribution in Fig. 4(D) is broad and symmetric but the mean is not reported.","section":"Results, Wall Shear Stress"},{"comment":"The 'regions of interest' (ROI) in Fig. 3(A) are only defined in the caption; the main text should state the axial extent of the ROI and how the cross-sectional number density is normalized.","section":"Results, Cell Distribution and Segregation"},{"comment":"The verification at Re = 0.05 is mentioned only as a sentence without any results; a supplementary figure or a quantitative statement of the change in ηN would make this verification credible.","section":"Formulation"}],"recommendation":"major_revision","confidential_remarks":"The manuscript stems from a group with a strong track record in this simulation methodology, and the healthy-cell validation against the Zweifach-Fung data is a clear asset. The concern for the editor is that the headline anti-Zweifach-Fung effect is tied to an extreme and explicitly acknowledged simplification of the sickle cell phenotype, and no sensitivity analysis is provided. This is a fixable gap rather than a fundamental error: adding parameter sweeps and uncertainty quantification would substantially strengthen the paper. I would also note that the paper's data availability statement says all data are in the paper, but the authors should specify whether raw simulation trajectories are available for reanalysis, especially given the absence of replicate statistics."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The paper has one central idea: normal RBCs show the classic Zweifach-Fung effect at a bifurcation, while stiff, marginated sickle cells do the opposite and preferentially enter the low-flow branch. That anti-Zweifach effect is mechanistically explained through the cell-free layer and the shifted separatrix, and the explanation is clean and convincing. The authors also show the segregation persists downstream, with sickle cells piling up on outer walls and raising the frequency of high wall shear stress events. That is new context — the malaria work of Li et al. already showed stiff cells going to low-flow branches, but the SCD-specific binary mixture, the ηN versus ηQ mapping, the downstream wall concentration statistics, and the WSS probability distributions are all new content. The healthy-cell validation against the experimental Zweifach-Fung data is solid, and the model is based on an established immersed boundary method.\n\nThe soft spots are in the quantitative extrapolation. The sickle cell model uses a rest volume about 20% of a normal RBC. Real dense sickle cells are typically 60-80% of normal volume, so 20% is extreme, and smaller and stiffer particles marginate more strongly. The headline anti-Zweifach effect depends on that strong margination. The authors state this is a simplification and they do cite the heterogeneity of real SCD populations, but they run no sensitivity analysis on volume, stiffness, or sickled fraction. Without that, we don't know if the magnitude of the effect (or even its existence) holds at physiologically representative sickle cell properties. The qualitative effect is probably robust — the malaria-infected cell simulations also showed stiff cells preferring low-flow branches — but the quantitative curves in Figures 2 and 5 need to be seen as illustrative, not predictive. Also, all partitioning and WSS statistics are presented without error bars or replicate runs, which makes it hard to judge the significance of the differences they report. The code and data are not released, despite a data availability statement that says \"all data needed... are present.\" \n\nThe citation pattern is honest; they explicitly cite Li et al. for the prior stiff-cell-to-low-flow-branch observation and position their contribution as extending it. No signs of overclaiming novelty on that front.\n\nWho should read this? Microcirculation modelers and people who care about the physics of cell segregation in bifurcating vessels. It is a useful computational study, not a paradigm shift. It deserves a serious referee, but the referee should push for sensitivity analysis and error bars. If the authors can show the effect survives at more realistic sickle cell volumes, the paper becomes much stronger.\n\nRecommendation: send to peer review, with a request for sensitivity runs and uncertainty quantification.","headline":"Stiff sickle cells preferring the low-flow branch is a plausible and nicely explained qualitative result, but the extreme 20%-volume cell model and lack of error bars cap how far the quantitative claims should be trusted.","tokens_in":14266,"tokens_out":2307,"would_cite":true,"duration_ms":25729,"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":"Simulations show stiff sickle cells prefer the low-flow branch at vascular bifurcations, reversing the Zweifach-Fung effect for normal red cells.","keywords":["sickle cell disease","red blood cell margination","Zweifach-Fung effect","anti-Zweifach-Fung effect","cell-free layer","wall shear stress","vascular bifurcation","computational simulation"],"falsifier":"Measure, in a microfluidic Y-junction with physiological flow splits and a suspension of stiffened or sickled red cells, the cell partition ratio into the high-flow branch: if the ratio equals or exceeds the fluid partition ratio for the high-flow branch at flow splits of 0.7 and 0.9, the central claim is contradicted. A complementary observation would be imaging the parent-vessel cell distribution and showing that stiff cells do not concentrate in the cell-free layer upstream of the bifurcation.","tokens_in":13194,"feed_emoji":"🩸","tokens_out":4433,"duration_ms":38734,"temperature":0.7,"pith_summary":"This paper uses three-dimensional simulations of a binary red blood cell suspension to argue that stiff sickle-shaped cells reverse a well-known microcirculation rule. Healthy red cells drift to the vessel center and, at a fork, preferentially enter the higher-flow branch—the Zweifach-Fung effect. Sickle cells, being stiffer and smaller, marginate to the cell-free layer near the wall, and when the flow split is unequal they preferentially enter the lower-flow branch instead. The paper further shows this upstream segregation persists downstream, concentrates sickle cells on outer branch walls, and increases the frequency of high wall-shear-stress events, offering a physical route by which blood disorders damage endothelium.","feed_headline":"Sickle cells take the slow branch at blood-vessel forks","feed_subtitle":"Stiff, marginated sickle cells reverse the Zweifach-Fung effect and raise wall shear stress on outer branch walls.","key_machinery":"The load-bearing object is the separatrix, the dividing surface in the parent vessel that marks whether a cell's trajectory ends in the left or right daughter branch, together with the cell-free layer (CFL), the near-wall region depleted of healthy red cells. The argument works because healthy cells are concentrated at the channel center, sickle cells are marginated in the CFL, and shifting the separatrix under unequal flow sends the CFL preferentially into the low-flow branch. The sickle cell itself is modeled as a stiff capsule with a curved prolate spheroidal rest shape, roughly 20% of normal cell volume and five times the membrane shear modulus.","core_discovery":"The central claim is that margination inverts the Zweifach-Fung effect for aberrant red blood cells. In a simulated binary suspension of 90% normal and 10% sickle cells flowing through a symmetric bifurcation, normal cells follow the classical rule (cell partition ratio exceeds fluid partition ratio when a branch receives more flow), while sickle cells show an anti-Zweifach-Fung effect: they preferentially enter the low-flow branch. The mechanism is spatial: upstream of the fork, normal cells occupy the channel center while sickle cells reside in the cell-free layer, and the separatrix—the surface dividing fluid that enters each daughter branch—shifts toward the high-flow side when flow is asymmetric, sweeping the near-wall cell-free layer into the low-flow branch. The segregation persists downstream, with sickle cells accumulating on the outer walls of daughter branches, and this accumulation increases the probability of large wall-shear-stress events, particularly on the outer side of the high-velocity branch. In geometrically asymmetric bifurcations, cells of both types favor the larger-radius branch at equal flow, and sickle cells still exhibit the anti-Zweifach-Fung effect.","pith_inferences":["If real sickle-cell populations contain many cells that only partially stiffen, the magnitude of the anti-Zweifach-Fung effect may be weaker or distributed across a spectrum of partition ratios, since the current model uses a single idealized sickle cell.","The same separatrix-shift logic should apply to other marginated blood elements such as platelets, white cells, or rigid inclusions, so bifurcations may generally sort stiff or small particles into low-flow branches—a possible design principle for microfluidic separators.","A testable extension is to measure the cell partition ratio at varying flow splits in a microfluidic Y-junction with stiffened or sickled cells; the crossover where the cell partition ratio drops below the fluid partition ratio for the high-flow branch would confirm the mechanism.","The wall-shear-stress increase is reported as event frequency at fixed points; mapping these events to endothelial calcium signals or adhesion-molecule expression would connect the physical mechanism to tissue-level pathology."],"forward_implications":["In mixed sickle-cell blood, downstream branches with lower flow will receive a disproportionate share of sickle cells, concentrating stiff cells where flow is already sluggish.","The outer walls of daughter branches, especially the high-velocity branch, see more high wall-shear-stress events in sickle-cell suspensions than in healthy suspensions.","Geometric asymmetry alone biases cell entry toward the larger-radius branch even when flow rates are equal, and this bias compounds the flow-driven partitioning.","The anti-Zweifach-Fung effect provides a purely physical, non-adhesive mechanism connecting margination to endothelial stress, complementing experiments showing that sickle-cell margination upregulates an endothelial dysfunction marker."],"supporting_citations":[{"why":"Defines the Zweifach-Fung effect that the paper extends and then inverts for marginated sickle cells.","marker":"1,2"},{"why":"Supplies the experimental observation that sickle cells marginate and that this margination upregulates endothelial dysfunction markers.","marker":"3,4"},{"why":"Documents the heterogeneous distribution of sickle-cell subpopulations, with the stiffest cells segregating toward channel walls.","marker":"15"},{"why":"Provides the computational model for deformable capsules and the wall-shear-stress framework used throughout the study.","marker":"18"},{"why":"Quantitatively explains the Zweifach-Fung effect through the inlet distribution of cells and the separatrix shift.","marker":"21"},{"why":"Shows that malaria-infected stiffened cells preferentially enter the low-flow branch, an earlier analog of the anti-Zweifach-Fung behavior.","marker":"48"},{"why":"Supplies the experimental data against which the simulation's Zweifach-Fung partitioning is validated.","marker":"61"},{"why":"Reports an inversion of the Zweifach-Fung effect in less deformable red cells at low hematocrit, providing precedent for deviation from the classical rule.","marker":"62"}],"fun_headline_variants":["Sickle cells defy flow rule at vessel forks","Sickle cells take the wrong fork, raising damage risk","Inverted Zweifach-Fung: sickle cells choose slow branch","Marginated sickle cells flip blood-flow partitioning"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"Everything hinges on the idealized sickle cell: stiff, small, and curved enough to marginate strongly into the cell-free layer; if real sickle cells in a patient are softer or heterogeneous enough not to marginate, the anti-Zweifach-Fung effect and the wall-stress consequences would weaken or disappear.","fun_headline_variants_meta":{"raw":{"variants":["Sickle cells defy flow rule at vessel forks","Sickle cells take the wrong fork, raising damage risk","Inverted Zweifach-Fung: sickle cells choose slow branch","Marginated sickle cells flip blood-flow partitioning"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000253,"raw_usage":{"total_tokens":1617,"prompt_tokens":1049,"completion_tokens":568,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":665,"completion_tokens_details":{"reasoning_tokens":504}},"tokens_in":665,"tokens_out":568,"duration_ms":5405,"temperature":1.0,"reasoning_tokens":504,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-11T05:08:20.728957+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"Measure, in a microfluidic Y-junction with physiological flow splits and a suspension of stiffened or sickled red cells, the cell partition ratio into the high-flow branch: if the ratio equals or exceeds the fluid partition ratio for the high-flow branch at flow splits of 0.7 and 0.9, the central claim is contradicted. A complementary observation would be imaging the parent-vessel cell distribution and showing that stiff cells do not concentrate in the cell-free layer upstream of the bifurcation.","supporting_citations":[{"cited_title":"As depicted in Fig","cited_arxiv_id":null,"evidence_quote":"Provides the computational model for deformable capsules and the wall-shear-stress framework used throughout the study."}],"review_version":1}