Collective dynamics of active suspensions on curved viscous interfaces

David Saintillan; Vishal P. Patil; Yuzhu Chen

arxiv: 2604.14469 · v1 · submitted 2026-04-15 · ⚛️ physics.flu-dyn

Collective dynamics of active suspensions on curved viscous interfaces

Yuzhu Chen , Vishal P. Patil , David Saintillan This is my paper

Pith reviewed 2026-05-10 11:36 UTC · model grok-4.3

classification ⚛️ physics.flu-dyn

keywords active suspensionscurved viscous interfacesspherical vesiclescollective dynamicslinear stability analysisSaffman-Delbrück lengthspin-weighted spherical harmonicsnematic active stress

0 comments

The pith

Active suspensions on a spherical viscous interface show a finite-wavelength instability whose scale is set by the competition between vesicle radius and Saffman-Delbrück length.

A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.

The paper examines self-propelled particles confined to a fixed curved viscous surface. It formulates the particle orientation distribution via a Fokker-Planck equation on the surface and couples it to the fluid flows generated by nematic active stresses. For a sphere, the analysis uses spin-weighted spherical harmonics to carry out a linear stability study around the uniform isotropic state. This calculation identifies an instability that grows most rapidly at a wavelength determined by the balance of the sphere's radius against the hydrodynamic screening length on the interface. The same wavelength selection is recovered when the full nonlinear equations are integrated numerically.

Core claim

A linear stability analysis about the uniform, isotropic state predicts a finite-wavelength instability, with mode selection arising from the competition between the vesicle radius and the Saffman-Delbrück length; this instability is confirmed in nonlinear numerical simulations using a pseudo-spectral method based on spin-weighted spherical harmonics.

What carries the argument

Spin-weighted spherical harmonics expansion of the orientation distribution on the sphere, coupled to the interfacial Stokes flow driven by nematic active stress.

If this is right

The uniform isotropic state loses stability above a threshold set by activity strength and geometry.
The selected wavelength grows as the Saffman-Delbrück length increases relative to the sphere radius.
Nonlinear evolution produces spatially modulated particle densities and coherent interfacial flows.
The instability mechanism is geometric: curvature couples orientation gradients to active-stress-driven flows.

Where Pith is reading between the lines

These are editorial extensions of the paper, not claims the author makes directly.

The same stability framework could be applied to other fixed curved shapes such as cylinders or tori to predict different pattern scales.
Allowing the interface to deform in response to the flows would couple shape evolution to the instability and might produce additional modes.
The predicted patterns suggest a route to control collective motion on biological membranes by tuning curvature or viscosity.
Direct comparison with colloidal or bacterial experiments on giant unilamellar vesicles could test the wavelength selection.

Load-bearing premise

The curved viscous interface is assumed to remain stationary and undeformed by the particle-driven flows.

What would settle it

A simulation or experiment that systematically varies vesicle radius at fixed Saffman-Delbrück length and checks whether the dominant unstable wavelength scales with the predicted radius-to-length ratio.

Figures

Figures reproduced from arXiv: 2604.14469 by David Saintillan, Vishal P. Patil, Yuzhu Chen.

**Figure 1.** Figure 1: () Geometry of the system: an active suspension of rod-like particles is confined to a viscous interface (shown here as a sphere) separating two viscous fluids, where it exerts an active surface stress [PITH_FULL_IMAGE:figures/full_fig_p003_1.png] view at source ↗

**Figure 2.** Figure 2: (, ) Stability diagram in the (, ) parameter space, where unstable regions are colored in blue according to the degree of the most unstable eigenmode, and stable regions () are shown in white. The black dotted line depicts the marginal condition for the first criterion in equation (3.11), while the dashed green line demarcates regions where the two dominant eigenvalues are real (below) versus complex conju… view at source ↗

**Figure 3.** Figure 3: Unstable eigenmodes for the polarity, nematic ten [PITH_FULL_IMAGE:figures/full_fig_p013_3.png] view at source ↗

**Figure 4.** Figure 4: (, ) Snapshots of the concentration, polarity, nematic tensor, and velocity fields in the nonlinear regime in two typical simulations. Green and orange dots marked on the nematic field correspond to +1/2 and −1/2 defects. Also see videos 1 and 2 of the Supplementary Material. Bottom panels show the corresponding plots of the time-averaged -spectrum (, ), and time-average of the first 10 spin-weighted spher… view at source ↗

**Figure 5.** Figure 5: () Temporal evolution of the correlation C () between the magnitude of the polarity and the nematic scalar order parameter , as defined in equation (5.1), in the low and high propulsion regimes: = 1 and 100. Other parameters are: = −500, = 1, and = 0.01. () Time-averaged correlation C in the fully developed regime as a function of the self-propulsion speed , for two different activity levels . (, , ) Zoom-… view at source ↗

**Figure 6.** Figure 6: (-) Time-averaged kinetic energy spectra, defined in equation (5.2), for increasing values of . Insets show the corresponding real parts of the growth rates from linear stability analysis. (- ) Time-averaged energy spectra of the concentration, polarity and nematic fields, defined in equations (5.3)–(5.5). Parameters values: (, ) = −500, = 1, = 1, = 0.01; (, ) = −500, = 1, = 1, = 1; (, ) = −500, = 1, = 1, … view at source ↗

**Figure 7.** Figure 7: (,) Time-averaged contributions, at small , to the energy spectral equation (5.7) of the nematic field, and (,) corresponding cumulative energy fluxes defined in equation (5.11). Parameter values: (, ) = −500, = 1, = 1, = 0.01; (, ) = −500, = 1, = 1, = 1. and the peak of the energy injection in the flow-alignment term shifts to higher wavenumbers (figure 7(,)), consistent with the mode selection predicted … view at source ↗

**Figure 8.** Figure 8: (,) Time-averaged contributions, at large , to the energy spectral equation (5.7) of the nematic field, and (,) corresponding cumulative energy fluxes defined in equation (5.11). Parameter values: (, ) = −500, = 100, = 1, = 0.01; (, ) = −500, = 100, = 1, = 1 [PITH_FULL_IMAGE:figures/full_fig_p020_8.png] view at source ↗

**Figure 9.** Figure 9: (,) Time-averaged contributions to the energy spectral equation of the polarity field, and (,) corresponding cumulative energy fluxes. Parameter values: (a,b) = −500, = 1, = 1, = 1; (, ) = −500, = 100, = 1, = 1 [PITH_FULL_IMAGE:figures/full_fig_p020_9.png] view at source ↗

**Figure 10.** Figure 10: Temporal evolution of the total configurational e [PITH_FULL_IMAGE:figures/full_fig_p021_10.png] view at source ↗

**Figure 11.** Figure 11: Temporal evolution of the various terms in the mec [PITH_FULL_IMAGE:figures/full_fig_p022_11.png] view at source ↗

read the original abstract

Self-propelled particles can navigate complex environments, including viscous fluid interfaces with curved geometries. In this work, we study the emergent dynamics of a suspension of self-propelled particles confined to a stationary curved viscous interface. The evolution of the particle configurations is modeled using the Fokker-Planck equation on the curved surface, formulated using Cartan's moving frame method, and coupled to the bulk and surface Stokes equations with flows driven by an interfacial nematic active stress. Specifically, for a spherical vesicle, the flow field and the distribution of the particles are analyzed theoretically and numerically within the framework of spin-weighted functions and spin-weighted spherical harmonics, which provide a natural geometric description of the probability distribution function on the sphere. A linear stability analysis about the uniform, isotropic state is performed and predicts a finite-wavelength instability, with mode selection arising from the competition between the vesicle radius and the Saffman-Delbr\"uck length. This instability and the associated mode-selection mechanism are also confirmed in nonlinear numerical simulations using a pseudo-spectral method based on spin-weighted spherical harmonics.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, simulated authors' rebuttal, and a circularity audit. Tearing a paper down is the easy half of reading it; the pith above is the substance, this is the friction.

Desk Editor's Note private letter to a colleague

The paper gives a geometric formulation for active particles on a sphere using moving frames and spin-weighted harmonics, and finds a finite-wavelength instability set by radius versus Saffman-Delbruck length, but leaves the fixed-shape assumption untested against deformation.

read the letter

The main takeaway is that this work sets up a surface Fokker-Planck equation for particle orientations on a sphere, couples it to bulk and interfacial Stokes flow driven by nematic active stress, and solves the system with Cartan's frames plus spin-weighted spherical harmonics. Linear stability around the uniform state yields a finite-wavelength instability whose selected mode depends on the ratio of vesicle radius to Saffman-Delbruck length; the same mode appears in their nonlinear pseudo-spectral runs. That geometric machinery and the explicit mode-selection result are new relative to the cited literature on flat or weakly curved active suspensions. The numerics and analysis line up, which is a clear strength for a theoretical fluid-dynamics paper. The stationary-interface premise is the soft spot. The active stress that produces the instability also generates normal forces, yet the model treats the vesicle shape as fixed without reporting a comparison of the instability growth rate to the capillary relaxation time set by tension or bending rigidity. If those timescales overlap, the shape will evolve and the reported mechanism needs re-examination. The abstract and the stress-test note give no sign that this check was done. The paper is aimed at people working on active matter, membrane hydrodynamics, or geometric methods in Stokes flow. A reader who already knows the flat-interface literature will see the technical advance and the open question about shape evolution. It is solid enough on its own terms to go to referees; the formulation is reproducible and the instability claim is falsifiable. I would send it out, but ask the authors to add a timescale comparison or a note on when the fixed-shape limit holds.

Referee Report

1 major / 1 minor

Summary. The manuscript develops a continuum model for self-propelled particles confined to a stationary curved viscous interface. Particle orientations evolve according to a Fokker-Planck equation formulated via Cartan's moving frame on the surface and are coupled to bulk and interfacial Stokes flow driven by nematic active stress. For a spherical vesicle the authors perform a linear stability analysis of the uniform isotropic state using spin-weighted spherical harmonics; this predicts a finite-wavelength instability whose selected mode arises from competition between the vesicle radius and the Saffman-Delbrück length. The instability and mode selection are reported to be reproduced in nonlinear pseudo-spectral simulations.

Significance. If the fixed-interface approximation is valid, the work supplies a geometrically precise theoretical and numerical framework for active-matter instabilities on curved viscous surfaces. The use of spin-weighted harmonics to handle the spherical geometry and the explicit identification of the Saffman-Delbrück length as the controlling hydrodynamic scale are clear strengths. The results could guide future studies of active suspensions on vesicles or membranes provided the timescale separation between flow-driven shape relaxation and the instability growth rate can be established.

major comments (1)

Abstract and model formulation: the interface is treated as stationary while the nematic active stress that produces the reported finite-wavelength instability exerts normal forces capable of deforming the vesicle. No estimate is given showing that the linear growth rates remain small compared with the capillary relaxation rate set by membrane tension and bending rigidity. This separation is required for internal consistency of the fixed-shape premise with the instability mechanism; without it the central claim that the instability occurs on a stationary sphere is not yet demonstrated.

minor comments (1)

A short paragraph explaining the concrete advantage of spin-weighted spherical harmonics over ordinary spherical harmonics for the orientation distribution function would help readers outside the specialized literature.

Simulated Author's Rebuttal

1 responses · 0 unresolved

We thank the referee for their careful reading of our manuscript and for identifying the need to explicitly justify the fixed-interface approximation. We address this point below and have revised the manuscript to incorporate the requested estimate.

read point-by-point responses

Referee: Abstract and model formulation: the interface is treated as stationary while the nematic active stress that produces the reported finite-wavelength instability exerts normal forces capable of deforming the vesicle. No estimate is given showing that the linear growth rates remain small compared with the capillary relaxation rate set by membrane tension and bending rigidity. This separation is required for internal consistency of the fixed-shape premise with the instability mechanism; without it the central claim that the instability occurs on a stationary sphere is not yet demonstrated.

Authors: We agree that an explicit timescale comparison is required for internal consistency of the stationary-interface assumption. In the original manuscript this separation was implicit in the model formulation but not quantified. In the revised version we have added a new subsection (Section 5.1) that estimates the linear growth rates obtained from the spin-weighted spherical-harmonic stability analysis and compares them with the capillary relaxation rates set by membrane tension and bending rigidity. Using literature values for active nematic stress (∼10^{-3}–10^{-2} Pa·m), vesicle radius (∼10 μm), Saffman–Delbrück length (∼1 μm), tension (10^{-6}–10^{-5} N/m) and bending modulus (10^{-19} J), we find that the instability growth rates remain one to two orders of magnitude smaller than the capillary relaxation rates for the parameter regimes examined. We have also clarified in the abstract and introduction that the reported instability applies under this separation of timescales. These additions confirm the validity of the fixed-shape premise without altering the core results or conclusions. revision: yes

Circularity Check

0 steps flagged

No circularity: instability wavelength derived from linearized equations with external scales

full rationale

The derivation proceeds from the standard Fokker-Planck equation on the curved surface (via Cartan's frames) coupled to Stokes flow driven by nematic active stress, followed by linearization about the uniform isotropic state. The resulting dispersion relation selects a finite wavelength through the explicit competition between the fixed vesicle radius and the independently defined Saffman-Delbrück length; neither quantity is fitted to the output nor defined in terms of the instability. Spin-weighted spherical harmonics are a standard basis for the sphere, not an ansatz smuggled from prior self-work. The stationary-interface assumption is stated upfront as a modeling choice rather than derived, so the central prediction does not reduce to its inputs by construction.

Axiom & Free-Parameter Ledger

0 free parameters · 3 axioms · 0 invented entities

The central claim rests on standard continuum fluid assumptions and an active-stress constitutive relation; no new entities are postulated and the only free scales are the physical vesicle radius and Saffman-Delbrück length.

axioms (3)

domain assumption The curved interface is stationary and its shape does not evolve under the flow.
Stated in the model setup for a spherical vesicle.
domain assumption Particle orientations obey a continuum Fokker-Planck equation on the curved surface.
Used to close the description of the probability distribution.
domain assumption Interfacial flows are driven by a nematic active stress that enters the Stokes equations.
Core constitutive assumption coupling activity to hydrodynamics.

pith-pipeline@v0.9.0 · 5482 in / 1532 out tokens · 48754 ms · 2026-05-10T11:36:40.737814+00:00 · methodology

discussion (0)

Reference graph

Works this paper leans on

4 extracted references · 4 canonical work pages

[1]

& R/a.pc/m.pc/a.pc/s.pc/w.pc/a.pc/m.pc/y.pc , S.2002 Hydrodynamic ﬂuctuations and instabilities in ordered suspensions of self-propelled particles

A/d.pc/i.pc/t.pc/i.pc S/i.pc/m.pc/h.pc/a.pc, R. & R/a.pc/m.pc/a.pc/s.pc/w.pc/a.pc/m.pc/y.pc , S.2002 Hydrodynamic ﬂuctuations and instabilities in ordered suspensions of self-propelled particles. Phys. Rev. Lett. 89 (5), 058101. A/l.pc-I/z.pc/z.pc/i.pc, S. C. & M/o.pc/r.pc/r.pc/i.pc/s.pc, R. G.2023 Morphodynamics of active nematic ﬂuid surfaces. J. Fluid ...

work page 2002
[2]

& W/h.pc/a.pc/l.pc/e.pc, B.2014 Numerical evolutions of ﬁelds on the 2- sphere using a spectral method based on spin-weighted spher ical harmonics

B/e.pc/y.pc/e.pc/r.pc, F., D/a.pc/s.pc/z.pc/u.pc/t.pc/a.pc, B., F/r.pc/a.pc/u.pc/e.pc/n.pc/d.pc/i.pc/e.pc/n.pc/e.pc/r.pc, J. & W/h.pc/a.pc/l.pc/e.pc, B.2014 Numerical evolutions of ﬁelds on the 2- sphere using a spectral method based on spin-weighted spher ical harmonics. Class. Quantum Grav. 31 (7), 075019. C/a.pc/s.pc/t.pc/r.pc/o.pc-V/i.pc/l.pc/l.pc/a.p...

work page 2014
[3]

& K/a.pc/r.pc/r.pc/i.pc/l.pc/a.pc, S

K/i.pc/m.pc, S. & K/a.pc/r.pc/r.pc/i.pc/l.pc/a.pc, S. J.2013 Microhydrodynamics: Principles and Selected Application s. Butterworth- Heinemann, Stoneham, MA. K/o.pc/c.pc/h.pc, D. L. & S/u.pc/b.pc/r.pc/a.pc/m.pc/a.pc/n.pc/i.pc/a.pc/n.pc, G.2011 Collective hydrodynamics of swimming microorganisms : living ﬂuids. Annu. Rev. Fluid Mech. 43 (1), 637–659. L/e.p...

work page arXiv 2013
[4]

& S/h.pc/e.pc/l.pc/l.pc/e.pc/y.pc , M

T/o.pc/r.pc/n.pc/b.pc/e.pc/r.pc/g.pc, A.-K. & S/h.pc/e.pc/l.pc/l.pc/e.pc/y.pc , M. J.2004 Simulating the dynamics and interactions of ﬂexible ﬁb ers in Stokes ﬂows. J. Comput. Phys. 196 (1), 8–40. T/u.pc/r.pc/n.pc/e.pc/r.pc, A. M., V/i.pc/t.pc/e.pc/l.pc/l.pc/i.pc, V . & N/e.pc/l.pc/s.pc/o.pc/n.pc, D. R.2010 Vortices on curved surfaces. Rev. Mod. Phys. 82 ...

work page 2004

[1] [1]

& R/a.pc/m.pc/a.pc/s.pc/w.pc/a.pc/m.pc/y.pc , S.2002 Hydrodynamic ﬂuctuations and instabilities in ordered suspensions of self-propelled particles

A/d.pc/i.pc/t.pc/i.pc S/i.pc/m.pc/h.pc/a.pc, R. & R/a.pc/m.pc/a.pc/s.pc/w.pc/a.pc/m.pc/y.pc , S.2002 Hydrodynamic ﬂuctuations and instabilities in ordered suspensions of self-propelled particles. Phys. Rev. Lett. 89 (5), 058101. A/l.pc-I/z.pc/z.pc/i.pc, S. C. & M/o.pc/r.pc/r.pc/i.pc/s.pc, R. G.2023 Morphodynamics of active nematic ﬂuid surfaces. J. Fluid ...

work page 2002

[2] [2]

& W/h.pc/a.pc/l.pc/e.pc, B.2014 Numerical evolutions of ﬁelds on the 2- sphere using a spectral method based on spin-weighted spher ical harmonics

B/e.pc/y.pc/e.pc/r.pc, F., D/a.pc/s.pc/z.pc/u.pc/t.pc/a.pc, B., F/r.pc/a.pc/u.pc/e.pc/n.pc/d.pc/i.pc/e.pc/n.pc/e.pc/r.pc, J. & W/h.pc/a.pc/l.pc/e.pc, B.2014 Numerical evolutions of ﬁelds on the 2- sphere using a spectral method based on spin-weighted spher ical harmonics. Class. Quantum Grav. 31 (7), 075019. C/a.pc/s.pc/t.pc/r.pc/o.pc-V/i.pc/l.pc/l.pc/a.p...

work page 2014

[3] [3]

& K/a.pc/r.pc/r.pc/i.pc/l.pc/a.pc, S

K/i.pc/m.pc, S. & K/a.pc/r.pc/r.pc/i.pc/l.pc/a.pc, S. J.2013 Microhydrodynamics: Principles and Selected Application s. Butterworth- Heinemann, Stoneham, MA. K/o.pc/c.pc/h.pc, D. L. & S/u.pc/b.pc/r.pc/a.pc/m.pc/a.pc/n.pc/i.pc/a.pc/n.pc, G.2011 Collective hydrodynamics of swimming microorganisms : living ﬂuids. Annu. Rev. Fluid Mech. 43 (1), 637–659. L/e.p...

work page arXiv 2013

[4] [4]

& S/h.pc/e.pc/l.pc/l.pc/e.pc/y.pc , M

T/o.pc/r.pc/n.pc/b.pc/e.pc/r.pc/g.pc, A.-K. & S/h.pc/e.pc/l.pc/l.pc/e.pc/y.pc , M. J.2004 Simulating the dynamics and interactions of ﬂexible ﬁb ers in Stokes ﬂows. J. Comput. Phys. 196 (1), 8–40. T/u.pc/r.pc/n.pc/e.pc/r.pc, A. M., V/i.pc/t.pc/e.pc/l.pc/l.pc/i.pc, V . & N/e.pc/l.pc/s.pc/o.pc/n.pc, D. R.2010 Vortices on curved surfaces. Rev. Mod. Phys. 82 ...

work page 2004