REVIEW 4 major objections 5 minor 54 references
Red blood cells aggregates transport for finite concentration
T0 review · 4 major / 5 minor · reviewed 2026-08-10 · deepseek-v4-flash
Pith's one-line read This paper claims that red blood cell transport through a channel is maximized at a moderate, reversible adhesion energy: weak adhesion raises cell flux, stronger adhesion blocks it.
desk verdict A plausible but statistically thin 2D simulation study; the non-monotonic flux result is interesting, but the physiological claim outruns the data. read the letter →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
What carries the argument
The machinery is a two-dimensional vesicle suspension in a straight channel, solved by a lattice Boltzmann method coupled to the membranes through an immersed boundary scheme. Each vesicle membrane combines stretching, area, and bending energies with a Lennard-Jones pair potential for adhesion between cells; the control parameter is $\bar{\epsilon}_{adh} = \epsilon_{adh} R_0^2 / k$, calibrated to fibrinogen levels by AFM measurements. The argument is carried by the competition between two effects of aggregates: small reversible aggregates remove dissipative recirculations between cells and increase flux, while large aggregates span the channel and block it. This competition produces the optimal adhesion energy at which $Q_c/Q_0$ is maximal.
What would settle it
Measure, in a microfluidic straight channel at fixed hematocrit and shear rate, the normalized red-cell flux as fibrinogen concentration is increased: if the flux does not first rise then fall, or if including a cytoskeleton in a 3D simulation removes the peak, the central claim is falsified.
Extended reading notes
Core claim
The central discovery is a non-monotonic dependence: for fixed capillary number and hematocrit, the maximal normalized cell flux $Q_c/Q_0$ first increases with dimensionless adhesion energy $\bar{\epsilon}_{adh}$, reaching a boost of roughly 20–30% over the no-adhesion case, and decreases beyond an optimal adhesion energy. Weak adhesion creates small reversible aggregates that suppress the recirculation zones between cells and lower dissipation, so more cells get carried by the same flow. Stronger adhesion makes aggregates grow laterally until they occupy a substantial part of the channel cross-section and reduce the flux. The effective viscosity shows the mirror image: it drops by up to about 50% as adhesion rises in the moderate range, then increases. The optimal adhesion energy shifts upward with capillary number because stronger flow breaks aggregates and requires stronger adhesion to maintain the benefit.
Load-bearing premise
The load-bearing premise is that a two-dimensional vesicle membrane without a cytoskeleton, with Lennard-Jones adhesion, captures the aggregate shapes and disaggregation behaviour of real red blood cells closely enough that the non-monotonic flux-vs-adhesion trend transfers to three dimensions; the paper itself flags the absent cytoskeleton as a limitation.
Editorial extensions
If this is right
- At fixed hematocrit and capillary number, both too little and too much adhesion reduce RBC transport; the peak flux sits at a nonzero adhesion energy.
- The optimal adhesion energy is larger at higher capillary numbers, because stronger flow dissociates aggregates and needs stronger adhesion to keep the beneficial suppression of recirculations.
- Effective suspension viscosity is non-monotonic in adhesion energy as well, with a minimum that can be about 50% lower than the no-adhesion viscosity.
- When viscosity contrast is high ($\lambda = 10$), aggregates are robust and hardly break under shear, and the cell flux is dramatically reduced compared with the reversible-aggregation regime.
Reading between the lines
- A direct testable extension: in a microfluidic channel at fixed hematocrit and shear rate, raising fibrinogen concentration should first increase and then decrease the measured RBC flux, with the peak shifting to higher fibrinogen at higher shear rates.
- If the mechanism survives in three dimensions, therapies for hyper-aggregation should aim at an optimal aggregation window rather than eliminating aggregation entirely, since moderate adhesion is predicted to be beneficial.
- The blockage term should depend on vessel diameter relative to aggregate size, so the optimal adhesion energy is expected to be lower in narrower vessels where even small aggregates occupy a large fraction of the cross-section; the paper's straight-channel geometry leaves this as a prediction for network studies.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This manuscript reports 2D lattice-Boltzmann simulations of vesicle suspensions as models of red blood cells in a straight channel, with a Lennard-Jones adhesion potential between membranes. The authors vary the dimensionless adhesion energy, capillary number, and volume fraction, and measure the normalized cell flux and effective viscosity. Their central claim is that the maximum cell flux is non-monotonic in adhesion energy: moderate adhesion (claimed to be physiological) increases flux by 20–30% relative to no adhesion, while stronger adhesion reduces it. The proposed mechanism is that small aggregates suppress dissipation in recirculation zones, while large aggregates block the channel. The paper also reports a corresponding minimum in relative viscosity and discusses reversible versus irreversible aggregation at low and high viscosity contrast.
Significance. If the non-monotonic flux result holds, it offers a plausible mechanistic explanation for the presence of reversible RBC aggregation in healthy circulation and a warning that excessive aggregation impairs perfusion. The study is systematic in exploring capillary number and volume fraction, and the no-adhesion limit is checked against the earlier results of Farutin et al. [18]. The simulation setup builds on prior validated work for doublet stability [9] and uses a published conversion from fibrinogen concentration to adhesion energy. The limitations of the 2D vesicle model and the absence of a cytoskeleton are acknowledged in Section IV, which tempers but does not eliminate the physiological relevance claim. The paper's main strength is the counter-intuitive, falsifiable prediction that moderate adhesion can improve transport, which should stimulate further work.
major comments (4)
- [Section III.B, Figs. 5–6] The central non-monotonic flux peak is reported from a single time-averaged simulation per parameter point, with no error bars, no ensemble averages, and no statement of the averaging duration or the number of cells. Since the suspension dynamics are chaotic and collective, the 0.3 to 0.4 increase in Qc/Q0 (quoted as a 30% boost) could be within run-to-run variability. Please provide error bars or confidence intervals from at least several independent initial configurations and describe the time-averaging protocol in Section II.C.
- [Section III.B and abstract/Introduction/Conclusion] The magnitude of the boost is stated as 'about 20%' in the Introduction (for physiological adhesion) and as '30%' in Section III.B and the Conclusion, while the abstract claims the boost occurs 'within the physiological range.' However, the 30% increase is attributed to ε̄adh ≈ 200 at Ca = 25, which exceeds the largest tabulated value ε̄adh = 147.7 in Table I, a value the table color-codes as pathological. Please reconcile these numbers, specify which data point corresponds to the physiological-range claim, and adjust the abstract and conclusion accordingly.
- [Section III.B, Fig. 5] The adhesion-energy scan per capillary number contains only about five hand-picked values, so the existence and the location of the 'optimal' adhesion energy in Fig. 6 are not robustly established. A finer scan with several values near the apparent peak, and across the same capillary numbers, is needed to support the non-monotonic claim and to define the optimal adhesion energy with confidence.
- [Section IV] The physiological conclusion that adhesion 'may be beneficial for perfusion' transfers the 2D vesicle result to real 3D RBCs, but the only support cited is the similarity of the doublet phase diagram in Ref. [9]. The recirculation-suppression and blockage mechanisms are sensitive to aggregate morphology and wall interactions, which differ between 2D and 3D. A concrete test would be to perform 3D simulations for the key peak parameters, or to compare with experimental measurements of flux versus aggregation, to assess whether the non-monotonic effect persists.
minor comments (5)
- [Fig. 9 caption] The caption states 'ϕ = 0.1%', which appears to be a typographical error; the text elsewhere uses volume fractions of 0.1 and 0.4, so this should probably read 'ϕ = 0.1'.
- [Eq. (5)] The equation 'uρ = ...' should be 'ρu = ...' (or 'u = .../ρ'); as written, the left-hand side mixes velocity and density.
- [Eqs. (14) and (15)] The notation is inconsistent: Eq. (14) uses Pc and P0, while Eq. (15) defines [η] with ϕ and the text then defines ϕt. Please standardize the concentration symbol throughout.
- [Section IV] There is a typo: 'cytoskeletton' should be 'cytoskeleton'; also, in Section III.C 'thee center' should be 'the center'.
- [Figs. 2, 3, 7, 8] The color coding is described only in the captions by listing adhesion-energy values; a single consistent legend across figures would help the reader compare the configurations.
Circularity Check
No significant circularity: the central flux curve is a direct simulation observable, not a fitted or self-referential quantity.
full rationale
The paper's central claim, that normalized RBC flux Qc/Q0 is non-monotonic in adhesion energy, comes from counting cells crossing a channel section in lattice-Boltzmann simulations (Figs. 5 and 6). The adhesion energy enters as an input parameter (Eq. 9), and the flux is an output; no parameter of the model is fitted to reproduce the peak, and the non-monotonic curve is not implied by the model definition alone. The comparison with Farutin et al. [18] in the zero-adhesion limit is a reproducible prior-work benchmark, and the use of the authors' earlier doublet study [9] for the adhesion-energy conversion and for the stability of high-viscosity-contrast aggregates is calibration and context rather than a derivation of the new flux result. The paper explicitly labels the 2D-to-3D extrapolation as a limitation rather than a derived consequence. Issues such as the absence of error bars, the sparse adhesion grid, and the 20% versus 30% discrepancy are statistical and internal-consistency concerns, not circularity.
Assumptions & free parameters
free parameters (7)
- Dimensionless adhesion energy scan points =
0, 5.2, 15.1, 27.5, 36.81, 53, 61.1, 73.63, 84.33, 105, 147.27, 158, 220.91, 530.18, 1060.37
- Capillary number values =
0.9, 9.0, 25.0, 90.0
- Hematocrit (volume fraction) values =
roughly 0.1 to 0.46
- Viscosity contrast values =
lambda = 1.0 and 10.0
- Reduced area =
0.65
- Membrane and adhesion micro-parameters =
not reported
- Adhesion energy conversion factor =
epsilon_adh = 1.6862 h epsilon (h not defined here)
assumptions (6)
- domain assumption A 2D vesicle model without a cytoskeleton is a sufficient stand-in for the essential mechanics of 3D red blood cells in this flow regime.
- domain assumption RBC-RBC adhesion is well represented by a Lennard-Jones potential between membrane nodes.
- domain assumption The measured AFM interaction energies at various fibrinogen levels can be mapped onto the 2D simulation's adhesion parameter.
- standard math The LBM plus immersed boundary with no-slip coupling gives an accurate solution of the Navier-Stokes equations for these suspensions.
- domain assumption Time-averaged measurements from a single dynamical run are stationary and representative of the suspension.
- domain assumption Poiseuille flow in a straight 2D channel is an adequate proxy for microcirculatory flow for the purpose of flux comparison.
Cite this review
Pith. "Pith review of Red blood cells aggregates transport for finite concentration." pith.science (2026). https://pith.science/paper/2RC3TCLP
@misc{pith2026250103860,
author = {Pith},
title = {Pith review of: Red blood cells aggregates transport for finite concentration},
year = {2026},
howpublished = {\url{https://pith.science/paper/2RC3TCLP}},
note = {Machine review of arXiv:2501.03860}
}
read the original abstract
Red blood cells (RBCs) are responsible for transporting oxygen and various metabolites to tissues and organs, as well as removing waste. Several cardiovascular diseases can impair these functions. For instance, in diabetes, increased RBC aggregation can lead to blood occlusion, thereby depriving tissues of efficient oxygen delivery. Interestingly, RBC adhesion occurs not only in disease states but also under physiological conditions, with the key difference being that adhesion is reversible in healthy situations. This paper focuses on numerical simulations in 2D, exploring different adhesion energies (both physiological and pathological) alongside varying flow strengths and hematocrit levels. A systematic analysis of RBC flux and viscosity is conducted. A remarkable finding is that moderate adhesion energy (within the physiological range) enhances RBC transport, thereby improving oxygen delivery to tissues. This provides insight into why RBC adhesion is present under normal conditions. Conversely, increasing adhesion energy beyond a certain point causes a collapse in RBC flux, thus reducing oxygen transport. We provide a basic explanation for the non-monotonic effect of adhesion energy on blood flow efficiency. This finding serves as an initial step in understanding the impact (both positive and negative) of RBC adhesion before addressing this issue in complex networks inspired by realistic microvasculatures.
Figures
Figures from the paper (12 more)
Reference graph
Works this paper leans on
-
[18]
A. Farutin, Z. Shen, G. Prado, V. Audemar, H. Ez-Zahraouy, A. Benyoussef, B. Polack, J. Harting, P. M. Vlahovska, T. Podgorski et al. , Physical review fluids , 2018, 3, 103603
work page 2018
- [9]
-
[1]
R. M. Bateman, M. D. Sharpe and C. G. Ellis, Critical care, 2003, 7, 1–15
work page 2003
-
[2]
K. L. M. M. B. L. P. F. C. M. R. M. H. A. K. S. V. B. M. B. G. F. P. F. Di Liberto, G. and P. Bartolucci, American journal of hematology , 2016, 91, 1008–1013
work page 2016
-
[3]
C. P. R. M. V. J. B. V. L. Y. Lemonne, N. and M. Etienne-Julan, American journal of hematology , 2012, 87, E129–E129 20
work page 2012
- [4]
-
[5]
P. Comeglio, S. Fedi, A. A. Liotta, A. P. Cellai, E. Chiarantini, D. Prisco, F. Mecacci, E. Parretti, G. Mello and R. Abbate, Thrombosis research, 1996, 84, 199–202
work page 1996
- [6]
Show all 54 references
-
[7]
Brooks, in Blood Cells, Rheology, and Aging , Springer, 1988, pp
D. Brooks, in Blood Cells, Rheology, and Aging , Springer, 1988, pp. 158–162
1988
-
[8]
Baskurt, B
O. Baskurt, B. Neu and H. J. Meiselman, 2011
2011
-
[10]
Hoore, F
M. Hoore, F. Yaya, T. Podgorski, C. Wagner, G. Gompper and D. A. Fedosov, Soft matter , 2018, 14, 6278–6289
2018
-
[11]
Flormann, O
D. Flormann, O. Aouane, L. Kaestner, C. Ruloff, C. Misbah, T. Podgorski and C. Wagner, Scientific reports, 2017, 7, 1–10
2017
-
[12]
Claver ´ ıa, O
V. Claver ´ ıa, O. Aouane, M. Thi´ ebaud, M. Abkarian, G. Coupier, C. Misbah, T. John and C. Wagner, Soft Matter , 2016, 12, 8235–8245
2016
-
[13]
H. H. Lipowsky, S. Usami and S. Chien, Microvascular research, 1980, 19, 297–319
1980
-
[14]
J. H. Barbee and G. R. Cokelet, Microvascular Research, 1971, 3, 17–21
1971
-
[15]
G. F. Birchard, American Zoologist, 1997, 37, 65–72
1997
-
[16]
M. S. Hedrick, D. A. Duffield and L. H. Cornell, Canadian Journal of Zoology , 1986, 64, 2081–2085
1986
-
[17]
Linderkamp, A
O. Linderkamp, A. A. Stadler and E. P. Zilow, Pediatric research, 1992, 32, 97–102
1992
-
[19]
Z. Gou, H. Zhang, M. Abbasi and C. Misbah, Biophysical Journal, 2021, 120, 4819–4831
2021
-
[20]
Mauer, S
J. Mauer, S. Mendez, L. Lanotte, F. Nicoud, M. Abkarian, G. Gompper and D. A. Fedosov, Physical review letters, 2018, 121, 118103
2018
-
[21]
Lanotte, J
L. Lanotte, J. Mauer, S. Mendez, D. A. Fedosov, J.-M. Fromental, V. Claveria, F. Nicoud, G. Gompper and M. Abkarian, Proceedings of the National Academy of Sciences , 2016, 113, 13289–13294
2016
-
[22]
Fischer and H
T. Fischer and H. Schmid-Sch¨ onbein, inRed Cell Rheology, Springer, 1978, pp. 347–361
1978
-
[23]
Minetti, V
C. Minetti, V. Audemar, T. Podgorski and G. Coupier, Journal of Fluid Mechanics, 2019, 864, 408–448
2019
-
[24]
T. M. Fischer and R. Korzeniewski, Journal of Fluid Mechanics , 2013, 736, 351–365
2013
-
[25]
Dupire, M
J. Dupire, M. Socol and A. Viallat, Proceedings of the National Academy of Sciences , 2012, 109, 20808–20813
2012
-
[26]
Abkarian, M
M. Abkarian, M. Faivre and A. Viallat, Physical review letters , 2007, 98, 188302
2007
-
[27]
Tomaiuolo, M
G. Tomaiuolo, M. Simeone, V. Martinelli, B. Rotoli and S. Guido, Soft Matter , 2009, 5, 3736–3740
2009
-
[28]
Abbasi, A
M. Abbasi, A. Farutin, A. Nait-Ouhra, H. Ez-Zahraouy, A. Benyoussef and C. Misbah, Physical Review Fluids, 2022, 7, 093603
2022
-
[29]
Kaoui, G
B. Kaoui, G. Biros and C. Misbah, Physical review letters , 2009, 103, 188101 21
2009
-
[30]
Kaoui, N
B. Kaoui, N. Tahiri, T. Biben, H. Ez-Zahraouy, A. Benyoussef, G. Biros and C. Misbah, Physical Review E, 2011, 84, 041906
2011
-
[31]
Biben, A
T. Biben, A. Farutin and C. Misbah, Physical Review E , 2011, 83, 031921
2011
-
[32]
Aouane, M
O. Aouane, M. Thi´ ebaud, A. Benyoussef, C. Wagner and C. Misbah, Physical Review E , 2014, 90, 033011
2014
-
[33]
Kaoui, A
B. Kaoui, A. Farutin and C. Misbah, Physical review. E , 2009, 80 6 Pt 1 , 061905
2009
-
[34]
Agarwal and G
D. Agarwal and G. Biros, Physical Review Fluids , 2020, 5, 013603
2020
-
[35]
Trozzo, G
R. Trozzo, G. Boedec, M. Leonetti and M. Jaeger, Journal of Computational Physics , 2015, 289, 62–82
2015
-
[36]
Succi, The lattice Boltzmann equation: for fluid dynamics and beyond , Oxford university press, 2001
S. Succi, The lattice Boltzmann equation: for fluid dynamics and beyond , Oxford university press, 2001
2001
-
[37]
C. K. Aidun and J. R. Clausen, Annual review of fluid mechanics , 2010, 42, 439–472
2010
-
[38]
P. L. Bhatnagar, E. P. Gross and M. Krook, Physical review, 1954, 94, 511
1954
-
[39]
Kr¨ uger, H
T. Kr¨ uger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva and E. M. Viggen,Springer International Publishing, 2017, 10, 4–15
2017
-
[40]
Kaoui and J
B. Kaoui and J. Harting, Rheologica Acta, 2016, 55, 465–475
2016
-
[41]
Guckenberger, A
A. Guckenberger, A. Kihm, T. John, C. Wagner and S. Gekle, Soft Matter , 2018, 14, 2032–2043
2018
-
[42]
Kaoui, A
B. Kaoui, A. Farutin and C. Misbah, Physical Review E , 2009, 80, 061905
2009
-
[43]
Thi´ ebaud, Z
M. Thi´ ebaud, Z. Shen, J. Harting and C. Misbah, Physical review letters , 2014, 112, 238304
2014
-
[44]
Tsubota, S
K.-i. Tsubota, S. Wada and T. Yamaguchi, Computer methods and programs in biomedicine , 2006, 83, 139–14
2006
-
[45]
Tsubota and S
K.-i. Tsubota and S. Wada, Physical Review E , 2010, 81, 011910
2010
-
[46]
Kaoui, J
B. Kaoui, J. Harting and C. Misbah, Physical Review E , 2011, 83, 066319
2011
-
[47]
Zhang, P
J. Zhang, P. C. Johnson and A. S. Popel, Physical biology, 2007, 4, 285
2007
-
[48]
Feng and E
Z.-G. Feng and E. E. Michaelides, Journal of computational physics , 2004, 195, 602–628
2004
-
[49]
Z. Shen, A. Farutin, M. Thi´ ebaud and C. Misbah, Physical Review Fluids , 2017, 2, 103101
2017
-
[50]
F ˚ ahraeus,Physiological reviews, 1929, 9, 241–274
R. F ˚ ahraeus,Physiological reviews, 1929, 9, 241–274
1929
-
[51]
Merrill, E
E. Merrill, E. Gilliland, G. Cokelet, H. Shin, A. Britten and R. Wells Jr, Journal of applied physiology , 1963, 18, 255–260
1963
-
[52]
Chien, The red blood cell , 1975, 2, 1031–1133
S. Chien, The red blood cell , 1975, 2, 1031–1133
1975
-
[53]
P. Menu, W. Bleeker, D. Longrois, A. Caron, B. Faivre-Fiorina, S. Muller, P. Labrude and J.-F. Stoltz, Transfusion, 2000, 40, 1095–1103
2000
-
[54]
J. K. Armstrong, R. B. Wenby, H. J. Meiselman and T. C. Fisher, Biophysical journal , 2004, 87, 4259–4270
2004
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.