REVIEW 1 major objections 6 minor 62 references
Active sorting to boundaries in active nematic -- passive isotropic fluid mixtures
T0 review · 1 major / 6 minor · reviewed 2026-08-05 · deepseek-v4-flash
Pith's one-line read Active nematic sorting at walls is controlled by anchoring, not wetting
desk verdict Solid simulation study showing anchoring-dependent boundary sorting in active-passive mixtures; the mechanism is convincing within the paranematic model, but the S0=0 assumption leaves a real gap for realistic bulk-nematic systems. 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 central object is the active stress Π_act = -ζϕQ, where ζ is the activity (positive extensile, negative contractile), ϕ is the active-fluid concentration, and Q is the nematic tensor. Its divergence creates body forces on the combined fluid; retaining spatial variations of concentration and nematic order gives a force ∝ ζ∇(ϕS)·(2nn-I) from order and concentration gradients and a force ∝ ζSϕ∇·(nn) from director curvature. The model is a two-fluid formulation in which a center-of-mass fluid and a relative flow are separated, with strong drag between components, plus Beris-Edwards dynamics for the nematic tensor.
What would settle it
Run the same two-fluid simulations with S0 raised to a finite equilibrium bulk order (e.g. S0 = 0.5) in a square box with planar and homeotropic anchoring: the paper's mechanism predicts the boundary enrichment and depletion for extensile activity should weaken or reverse; if the concentration profile is unchanged, the ∇(ϕS) force is not the controlling factor. Experimentally, in a microtubule-kinesin extensile nematic with a fluorescently labeled active component in a chamber with planar surface anchoring, measure the wall concentration: the paper predicts a stationary enrichment above bulk.
Extended reading notes
Core claim
In a two-fluid model of an active nematic mixed with a passive isotropic fluid and confined by walls, the paper claims that the active component sorts to the boundary purely through active stresses. With imposed anchoring, an extensile active fluid enriches the boundary layer when anchoring is planar and depletes it when anchoring is homeotropic; a contractile active fluid does the reverse. The same boundary accumulation occurs without imposed anchoring, because active flows spontaneously create planar active anchoring. The driving force is the divergence of the active stress, which produces forces from gradients of the product of concentration and nematic order, with direction set by the an
Load-bearing premise
The model sets the equilibrium nematic order of the bulk to zero, so the sorting relies on a strong contrast between an anchored, ordered boundary and a disordered bulk; if the bulk were already nematic, the order gradient that drives the effect would shrink and the predicted directions could change.
Editorial extensions
If this is right
- With imposed planar (homeotropic) anchoring, extensile activity enriches (depletes) the boundary layer, and contractile activity does the opposite.
- The same boundary sorting appears when anchoring is not imposed but generated by active flows, so a bare confining wall is enough for extensile material to accumulate.
- The active force normal to the wall, ζ∇(ϕS)(2(m·n)^2 - 1), directly ties anchoring angle and activity sign to the direction of sorting.
- In circular confinement, gradients in the magnitude of nematic order and gradients in director orientation contribute roughly equally to the radial force pushing active material outward.
- If the active species starts as a thermodynamically phase-separated droplet, activity drives it into a spontaneously rotating boundary ring, showing active forces can override equilibrium placement.
Reading between the lines
- If wall anchoring can be patterned, a single active species should be steerable: regions with planar anchoring attract extensile material, homeotropic patches repel it, so confinement geometry alone could route active components without chemical patterning.
- The contractile case may map onto actomyosin cell aggregates, where contractile stresses and homeotropic-like alignment at tissue boundaries would drive inward sorting, providing a mechanical complement to differential adhesion in explaining interior placement.
- Because activity sign is set by motor protein directionality, toggling motor action should reverse the wall concentration in the same chamber, a directly testable prediction of the stress-gradient mechanism.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. This paper presents lattice-Boltzmann simulations of a two-fluid model in which an active nematic component is mixed with a passive isotropic fluid and confined in square, channel, and circular geometries. The central finding is that an extensile active fluid accumulates at a boundary with planar anchoring and is depleted at a homeotropic boundary, with the opposite behavior for contractile activity. The authors attribute sorting to active forces generated by gradients in nematic order and concentration, and they show that the same mechanism works when boundary anchoring arises spontaneously from active flows. A circular-confinement analysis separates contributions from gradients in the magnitude of nematic order and gradients in director orientation, and a final demonstration shows an active droplet spreading into a rotating boundary ring.
Significance. If the proposed mechanism holds, it provides a minimal nonequilibrium route to boundary sorting that does not rely on differential adhesion, and it yields a clear, falsifiable design rule: the sign of activity and the anchoring angle determine whether the active component accumulates or depletes at the boundary. The paper has notable strengths: the result is measured from simulation rather than obtained by fitting; it uses multiple confinement protocols (free-slip and no-slip boxes, channels, circles); it reports time-averaged concentration profiles with standard deviations; it includes parameter sweeps in activity and radius; and the simulation code is publicly available on GitHub. The analytical decomposition in Sec. V uses standard active-stress expressions, and the supplementary movies and figures support the interpretation. The main caveat, discussed below, is that the model fixes the equilibrium nematic order to zero, making the bulk paranematic; the generality of the sorting mechanism outside this limit is not established.
major comments (1)
- [Sec. II, Eq. (16); Sec. IV] The model sets S0=0, so the bulk equilibrium state is isotropic. The sorting mechanism described in Sec. III relies on a gradient of nematic order S between the boundary and the bulk; with S0>0 this gradient is attenuated, and its sign can depend on whether the wall orders more or less strongly than the bulk. The predicted sorting direction, especially for contractile activity, is therefore untested outside the paranematic limit. Indeed, Sec. IV explicitly uses S0=0 to argue that contractile activity produces no active flows when λ>0. Because the abstract and discussion generalize to microtubule-kinesin networks and cell colonies, where bulk nematic order is common, the central claim is not established beyond the paranematic regime. I request either additional simulations with S0>0 or a clear qualification that the results apply to the paranematic limit of the model.
minor comments (6)
- [Sec. II, after Eq. (16)] The text says 'If a>0 and b>0' the two fluids are homogeneously mixed, but the simulation parameters list b=0. Clarify that b=0 (a purely quadratic Landau term) is also used and still corresponds to mixing.
- [Fig. 1(c) and Fig. 3(c)] The horizontal axis is labeled only 'activity'. Please indicate explicitly that positive values are extensile and negative values are contractile, and state the range of ζ plotted.
- [Sec. III, text before Fig. 2] Minor typo: 'as shown in in Fig. 2(a)' contains a duplicated 'in'.
- [Eqs. (19)-(20)] The derivation assumes a constant tilting angle Θ0. The text should note that this is an approximation and that near the boundary Θ varies, as acknowledged later when forder changes sign.
- [Sec. V, Fig. 4] The method for calculating the concentration enhancement near the wall is described in words, but the annulus width (or the criterion for choosing the peak) is not quantified. Please provide the numerical value used.
- [Sec. III, Fig. 2(b)-(c)] The flow arrows in the schematics are illustrative, but the text refers to 'flows due to an extensile active nematic' and 'contractile' without explaining the dipole direction. Adding a sentence defining the dipole force direction would help the reader.
Circularity Check
No significant circularity: the boundary-sorting result is measured in simulations and the analytical explanation derives from the active-stress identity.
full rationale
The paper's central claim—that extensile (contractile) active nematics preferentially accumulate at planar (homeotropic) boundaries and deplete for the opposite anchoring—is established by direct numerical simulation of a two-fluid model, not by fitting parameters to the outcome. The analytical explanation in Sec. III is an interpretation of the observed sorting in terms of the divergence of the active stress, Π_act = −ζϕQ (Eq. 14). The normal force formula quoted from Blow et al. [40] follows immediately from this stress: for a gradient ∇(ϕS) along the boundary normal, the force component is ζ|∇(ϕS)|(2(m·n)^2 − 1)m. This is a parameter-free identity, not an independent assumption that contains the sorting result. The same holds for the tangential active-anchoring force in Sec. IV. The model equations are taken from the authors' prior work ([46,48]), but they are the input dynamical framework, not a conclusion masquerading as support; the boundary sorting is an emergent behavior of these equations, verified by simulation. No fitted quantity is renamed as a prediction: the concentrations ϕ_edge are measured, not tuned. No uniqueness theorem is invoked to forbid alternatives, and no ansatz is smuggled in solely via self-citation. The choice S0=0 (Eq. 16) is a stated modeling assumption that creates the order gradient at the wall, but it is a limitation on the scope of the claims, not a circular step. The paper is self-contained against its own simulations and the analytical arguments are direct consequences of the model equations.
Assumptions & free parameters
free parameters (5)
- Landau mixing coefficients a and b =
a=0.0001, b=0
- Equilibrium nematic order S0 =
0
- Flow-alignment parameter lambda =
0.7
- Drag coefficient gamma =
0.1
- Activity zeta =
0.001 to 0.012 (swept)
assumptions (6)
- domain assumption Relative flow between the two fluids is negligible (|delta u| << |u_c|) because the drag gamma is the fastest relaxation mode.
- domain assumption Active stress on the nematic component is Pi_act = -zeta phi Q (Eq. 14).
- standard math The Beris-Edwards equation (Eq. 10) governs nematic order dynamics.
- domain assumption The boundary anchoring of the director is effectively infinitely strong (perfect anchoring) in Sec. III.
- domain assumption The free energy Eq. 16 with S0=0 and a>0, b=0 gives a homogeneous mixed paranematic equilibrium.
- standard math The expression for the normal active force at a boundary, F_norm = zeta |grad(phi S)| (2(m.n)^2 - 1) m, is valid.
Cite this review
Pith. "Pith review of Active sorting to boundaries in active nematic -- passive isotropic fluid mixtures." pith.science (2026). https://pith.science/paper/2SLKESRX
@misc{pith2026250901523,
author = {Pith},
title = {Pith review of: Active sorting to boundaries in active nematic -- passive isotropic fluid mixtures},
year = {2026},
howpublished = {\url{https://pith.science/paper/2SLKESRX}},
note = {Machine review of arXiv:2509.01523}
}
read the original abstract
We use a two-fluid model to study a confined mixture of an active nematic fluid and a passive isotropic fluid. We find that an extensile active fluid preferentially accumulates at a boundary if the anchoring is planar, whereas its boundary concentration decreases for homeotropic anchoring. These tendencies are reversed if the active fluid is contractile. We argue that the sorting results from gradients in the nematic order, and show that the behaviour can be driven by either imposed boundary anchoring or spontaneous anchoring induced by active flows. Our results can be tested by experiments on microtubule-kinesin motor networks, and may be relevant to sorting to the boundary in cell colonies or cancer spheroids.
Figures
Reference graph
Works this paper leans on
- [1]
-
[2]
Y. L. Huang, C. Shiau, C. Wu, J. E. Segall, and M. Wu, Biophysical Reviews and Letters 15, 131 (2020)
work page 2020
-
[3]
S. F. G. Krens and C.-P. Heisenberg, Current Topics in Developmental Biology 95, 189 (2011)
work page 2011
-
[4]
L. Balasubramaniam, A. Doostmohammadi, T. B. Saw, G. H. N. S. Narayana, R. Mueller, T. Dang, M. Thomas, S. Gupta, S. Sonam, A. S. Yap, Y. Toyama, R.-M. M` ege, J. M. Yeomans, and B. Ladoux, Nature Materials 20, 1156 (2021)
work page 2021
-
[5]
E. R. Oldewurtel, N. Kouzel, L. Dewenter, K. Henseler, and B. Maier, eLife 4, e10811 (2015)
work page 2015
-
[6]
M. S. Steinberg, Journal of Theoretical Biology 55, 431 (1975)
work page 1975
-
[7]
R. A. Foty and M. S. Steinberg, Developmental Biology 278, 255 (2005)
work page 2005
-
[8]
J. D. Tse, R. Moore, Y. Meng, W. Tao, E. R. Smith, and X.-X. Xu, BMC Developmental Biology 21, 2 (2021)
work page 2021
Show all 62 references
-
[9]
Heine, J
P. Heine, J. Lippoldt, G. A. Reddy, P. Katira, and J. A. K¨ as, New Journal of Physics23, 043034 (2021)
2021
-
[10]
Pawlizak, A
S. Pawlizak, A. W. Fritsch, S. Grosser, D. Ahrens, T. Thalheim, S. Riedel, T. R. Kießling, L. Oswald, M. Zink, M. L. Manning, and J. A. K¨ as, New Journal of Physics 17, 083049 (2015)
2015
-
[11]
M. C. Marchetti, J.-F. Joanny, S. Ramaswamy, T. B. Liverpool, J. Prost, M. Rao, and R. A. Simha, Reviews of Modern Physics 85, 1143 (2013)
2013
-
[12]
M. E. Cates and C. Nardini, Reports on Progress in Physics 88, 056601 (2025)
2025
-
[13]
M. E. Cates and J. Tailleur, Annual Review of Condensed Matter Physics 6, 219 (2015)
2015
-
[14]
Zhang, R
J. Zhang, R. Alert, J. Yan, N. S. Wingreen, and S. Granick, Nature Physics 17, 961 (2021)
2021
-
[15]
Gokhale, J
S. Gokhale, J. A. Li, A. Solon, J. Gore, and N. Fakhri, Physical Review E 105, 054605 (2022)
2022
-
[16]
Thutupalli, D
S. Thutupalli, D. Geyer, R. Singh, R. Adhikari, and H. A. Stone, Proceedings of the National Academy of Sciences of the United States of America 115, 5403 (2018)
2018
-
[17]
Furukawa, D
A. Furukawa, D. Marenduzzo, and M. E. Cates, Physical Review E 90, 022303 (2014)
2014
-
[18]
N. J. Lauersdorf, E. Nazockdast, and D. Klotsa, arXiv preprint arXiv:2407.07826 (2024)
2024 arXiv
-
[19]
C. P. Beatrici and L. G. Brunnet, Physical Review E 84, 031927 (2011)
2011
-
[20]
J. M. Belmonte, G. L. Thomas, L. G. Brunnet, R. M. C. de Almeida, and H. Chat´ e, Physical Review Letters100, 248702 (2008)
2008
-
[21]
Stenhammar, A
J. Stenhammar, A. Tiribocchi, R. J. Allen, D. Maren- duzzo, and M. E. Cates, Physical Review Letters 111, 145702 (2013)
2013
-
[22]
Wittkowski, A
R. Wittkowski, A. Tiribocchi, J. Stenhammar, R. J. Allen, D. Marenduzzo, and M. E. Cates, Nature Com- munications 5, 4351 (2014)
2014
-
[23]
C. A. Weber, C. H. Rycroft, and L. Mahadevan, Physical Review Letters 120, 248003 (2018)
2018
-
[24]
Tiribocchi, R
A. Tiribocchi, R. Wittkowski, D. Marenduzzo, and M. E. Cates, Physical Review Letters 115, 188302 (2015)
2015
-
[25]
J. N. Graham, G. M. Zhang, and J. M. Yeomans, Soft Matter 20, 2955 (2024)
2024
-
[26]
Rozman and J
J. Rozman and J. M. Yeomans, Physical Review Letters 133, 248401 (2024)
2024
-
[27]
X. B. Yang, M. L. Manning, and M. C. Marchetti, Soft Matter 10, 6477 (2014)
2014
-
[28]
Neville, J
L. Neville, J. Eggers, and T. B. Liverpool, Soft Matter 20, 8395 (2024)
2024
-
[29]
P´ erez-Gonz´ alez, R
C. P´ erez-Gonz´ alez, R. Alert, C. Blanch-Mercader, M. G´ omez-Gonz´ alez, T. Kolodziej, E. Bazellieres, J. Casademunt, and X. Trepat, Nature Physics 15, 79 (2019)
2019
-
[30]
A. P. Berke, L. Turner, H. C. Berg, and E. Lauga, Phys- ical Review Letters 101, 038102 (2008)
2008
-
[31]
S. E. Spagnolie and E. Lauga, Journal of Fluid Mechanics 700, 105–147 (2012)
2012
-
[32]
Doostmohammadi, J
A. Doostmohammadi, J. Ign´ es-Mullol, J. M. Yeomans, and F. Sagu´ es, Nature Communications9, 3246 (2018)
2018
-
[33]
T. B. Saw, A. Doostmohammadi, V. Nier, L. Kocgozlu, S. Thampi, Y. Toyama, P. Marcq, C. T. Lim, J. M. Yeo- mans, and B. Ladoux, Nature 544, 212 (2017)
2017
-
[34]
Duclos, C
G. Duclos, C. Blanch-Mercader, V. Yashunsky, G. Sal- breux, J.-F. Joanny, J. Prost, and P. Silberzan, Nature Physics 14, 728 (2018)
2018
-
[35]
R. A. Simha and S. Ramaswamy, Physical Review Letters 89, 058101 (2002)
2002
-
[36]
Voituriez, J.-F
R. Voituriez, J.-F. Joanny, and J. Prost, Europhysics Let- ters 70, 404 (2005). 9
2005
-
[37]
T. N. Shendruk, A. Doostmohammadi, K. Thijssen, and J. M. Yeomans, Soft Matter 13, 3853 (2017)
2017
-
[38]
Opathalage, M
A. Opathalage, M. M. Norton, M. P. N. Juniper, B. Langeslay, S. A. Aghvami, S. Fraden, and Z. Dogic, Proceedings of the National Academy of Sciences of the United States of America 116, 4788 (2019)
2019
-
[39]
Caballero, Z
F. Caballero, Z. You, and M. C. Marchetti, Soft Matter 19, 7828–7835 (2023)
2023
-
[40]
M. L. Blow, S. P. Thampi, and J. M. Yeomans, Physical Review Letters 113, 248303 (2014)
2014
-
[41]
H. Soni, W. Luo, R. A. Pelcovits, and T. R. Powers, Soft Matter 15, 6318–6330 (2019)
2019
-
[42]
Gulati, F
P. Gulati, F. Caballero, I. Kolvin, Z. You, and M. C. Marchetti, Soft Matter 20, 7703–7714 (2024)
2024
-
[43]
Adkins, I
R. Adkins, I. Kolvin, Z. H. You, S. Witthaus, M. C. Marchetti, and Z. Dogic, Science 377, 768 (2022)
2022
-
[44]
Alert, C
R. Alert, C. Blanch-Mercader, and J. Casademunt, Phys- ical Review Letters 122, 088104 (2019)
2019
-
[45]
L. Zhao, P. Gulati, F. Caballero, I. Kolvin, R. Adkins, M. C. Marchetti, and Z. Dogic, Proceedings of the Na- tional Academy of Sciences of the United States of Amer- ica 121, 2410345121 (2024)
2024
-
[46]
Bhattacharyya and J
S. Bhattacharyya and J. M. Yeomans, Physical Review Letters 130, 238201 (2023)
2023
-
[47]
J. F. Joanny, F. J¨ ulicher, K. Kruse, and J. Prost, New Journal of Physics 9, 422 (2007)
2007
-
[48]
Bhattacharyya and J
S. Bhattacharyya and J. M. Yeomans, Physical Review E 110, 024607 (2024)
2024
-
[49]
Malevanets and J
A. Malevanets and J. M. Yeomans, Faraday Discussions 112, 237 (1999)
1999
-
[50]
Malevanets and J
A. Malevanets and J. M. Yeomans, Computer Physics Communications 127, 105 (2000)
2000
-
[51]
P. G. De Gennes and J. Prost, The Physics of Liquid Crystals (Oxford University Press, 1993)
1993
-
[52]
A. N. Beris and B. J. Edwards, Thermodynamics of flow- ing systems: with internal microstructure (Oxford Uni- versity Press, 1994)
1994
-
[53]
Pokawanvit, Z
S. Pokawanvit, Z. Chen, Z. You, L. Angheluta, M. C. Marchetti, and M. J. Bowick, Physical Review E 106, 054610 (2022)
2022
-
[54]
Krueger, H
T. Krueger, H. Kusumaatmaja, A. Kuzmin, O. Shardt, G. Silva, and E. Viggen, The Lattice Boltzmann Method: Principles and Practice , Graduate Texts in Physics (Springer, 2016)
2016
-
[55]
R. C. V. Coelho, H. R. J. C. Figueiredo, and M. M. Telo da Gama, Physical Review Research 5, 033165 (2023)
2023
-
[56]
R. C. V. Coelho, N. A. M. Ara´ ujo, and M. M. Telo da Gama, Philosophical Transactions of the Royal Soci- ety A: Mathematical, Physical and Engineering Sciences 379, 20200394 (2021)
2021
-
[57]
Santhosh, M
S. Santhosh, M. R. Nejad, A. Doostmohammadi, J. M. Yeomans, and S. P. Thampi, Journal of Statistical Physics 180, 699 (2020)
2020
-
[58]
F. G. Woodhouse and R. E. Goldstein, Physical Review Letters 109, 168105 (2012)
2012
-
[59]
P. L. Townes and J. Holtfreter, Journal of Experimental Zoology 128, 53 (1955)
1955
-
[60]
Moore, K
R. Moore, K. Q. Cai, D. O. Escudero, and X.-X. Xu, genesis 47, 579 (2009), https://onlinelibrary.wiley.com/doi/pdf/10.1002/dvg.20536
2009 doi
-
[61]
Cochet-Escartin, T
O. Cochet-Escartin, T. T. Locke, W. H. Shi, R. E. Steele, and E.-M. S. Collins, Biophysical Journal 113, 2827 (2017)
2017
-
[62]
Batlle and D
E. Batlle and D. G. Wilkinson, Cold Spring Harbor Per- spectives in Biology 4, a008227 (2012). 10 Supplemental Notes: A1. MOVIE CAPTIONS The eight videos illustrate the results of our simulations, and are described below. Movie 1:Simulations of an extensile active nematic flui...
2012
Reviewed August 5, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.