Turbulence-Induced Fluctuating Interfaces in Heterogeneously-Active Suspensions
Pith reviewed 2026-05-23 02:31 UTC · model grok-4.3
The pith
Heterogeneous activity in bacterial suspensions generates fluctuating hydrodynamic interfaces between turbulent and jammed regions.
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
Core claim
In a hydrodynamical model of dense bacterial suspensions with quenched heterogeneous activity patterns, the evolution of the bacterial velocity field produces hydrodynamic interfaces that separate spatially localized turbulence from jammed frictional surroundings. These interfaces are characterized by intermittent and multiscale fluctuations, and the activity heterogeneity affects mixing as measured by the residence times of Lagrangian tracers.
What carries the argument
Hydrodynamic interfaces emerging from velocity field evolution under quenched activity patterning, which separate localized turbulence from jammed frictional surroundings.
Load-bearing premise
The chosen hydrodynamical model for dense bacterial suspensions with fixed-in-time activity patterns is sufficient to capture the relevant physics of real heterogeneous suspensions.
What would settle it
An experiment or simulation with heterogeneous activity in bacterial suspensions that fails to produce fluctuating interfaces separating turbulent zones from jammed regions would falsify the claim.
Figures
read the original abstract
We investigate the effects of heterogeneous (spatially varying) activity in a hydrodynamical model for dense bacterial suspensions, confining ourselves to experimentally realizable, simple, quenched, activity patterns. We show that the evolution of the bacterial velocity field under such activity patterning leads to the emergence of hydrodynamic interfaces separating spatially localized turbulence from jammed frictional surroundings. We characterise the intermittent and multiscale fluctuations of this interface and also investigate how heterogeneity influences mixing via the residence times of Lagrangian tracers. This work reveals how naturally occurring heterogeneities could decisively steer active flows into more complex configurations than those typically studied, opening up parallels to droplet dynamics, front propagation and turbulent mixing layers.
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
Summary. The manuscript investigates the effects of spatially heterogeneous but time-quenched activity patterns within a hydrodynamic model of dense bacterial suspensions. It reports that the evolution of the bacterial velocity field produces fluctuating hydrodynamic interfaces separating localized turbulent regions from jammed frictional surroundings. The work characterizes the intermittent, multiscale statistics of these interfaces and examines their consequences for mixing through the residence-time statistics of Lagrangian tracers. The central claim is that such experimentally realizable heterogeneities can steer active flows into configurations more complex than those in homogeneous systems, with suggested parallels to droplet dynamics and turbulent mixing layers.
Significance. If the reported interfaces and fluctuation statistics are robust, the result would be significant for understanding how natural spatial heterogeneities control the structure and transport properties of active turbulence. The restriction to simple, quenched patterns that are experimentally realizable is a constructive choice that keeps the study focused and potentially testable. The Lagrangian mixing analysis adds a practical dimension. However, the significance is tempered by the absence of direct validation against either experiments or self-consistent (time-evolving) activity fields.
major comments (3)
- [§3] §3 (Model and activity patterning): The emergence of the hydrodynamic interfaces is demonstrated exclusively for quenched (fixed-in-time) activity fields. No auxiliary simulations with slowly varying or self-consistent activity are presented to test whether the interface formation and its fluctuation spectrum survive when the activity pattern is allowed to evolve, which is required to establish that the reported interfaces are not an artifact of the quenching approximation.
- [§5] §5 (Lagrangian mixing): The residence-time distributions used to quantify mixing are shown for a single realization of the quenched pattern. Without ensemble averaging over multiple independent activity patterns or reported error bars on the residence-time statistics, it is not possible to assess the statistical significance of the claimed influence of heterogeneity on mixing.
- [§2] §2 (Hydrodynamic model): The continuum equations are stated without a direct comparison to experimental velocity statistics or interface widths in heterogeneous bacterial suspensions. Because the central claim rests on the model faithfully capturing steric, frictional, and polarity dynamics at the relevant scales, a quantitative benchmark against existing experimental data for even one heterogeneous pattern is needed to support the extrapolation to real suspensions.
minor comments (3)
- [Figure 1] Figure 1: the color bar for the activity field is not labeled with units or the precise functional form used to generate the pattern.
- [§4] Notation: the symbol for the interface position is introduced without an explicit definition in the text preceding the first use in §4.
- [Discussion] References: several recent experimental papers on heterogeneous active suspensions are cited only in passing; a short dedicated paragraph comparing the present quenched-pattern results to those experiments would improve context.
Simulated Author's Rebuttal
We thank the referee for their constructive comments. Our study is deliberately restricted to quenched activity patterns, as stated in the abstract and introduction, to examine experimentally realizable heterogeneities. We respond point-by-point to the major comments below.
read point-by-point responses
-
Referee: [§3] §3 (Model and activity patterning): The emergence of the hydrodynamic interfaces is demonstrated exclusively for quenched (fixed-in-time) activity fields. No auxiliary simulations with slowly varying or self-consistent activity are presented to test whether the interface formation and its fluctuation spectrum survive when the activity pattern is allowed to evolve, which is required to establish that the reported interfaces are not an artifact of the quenching approximation.
Authors: The manuscript explicitly confines the analysis to quenched patterns because these are simple and experimentally realizable, allowing isolation of spatial heterogeneity effects on the velocity field. The interfaces emerge from the hydrodynamic evolution under fixed activity; we make no claim about persistence under dynamic activity. Adding self-consistent simulations would require a substantially different model and lies outside the stated scope of the work. revision: no
-
Referee: [§5] §5 (Lagrangian mixing): The residence-time distributions used to quantify mixing are shown for a single realization of the quenched pattern. Without ensemble averaging over multiple independent activity patterns or reported error bars on the residence-time statistics, it is not possible to assess the statistical significance of the claimed influence of heterogeneity on mixing.
Authors: The residence-time results illustrate the qualitative effect of heterogeneity for a representative quenched pattern. We agree that ensemble statistics over multiple patterns would strengthen the quantitative claims. We will attempt to add results from additional independent realizations together with error bars if the computational data can be generated within the revision timeframe. revision: partial
-
Referee: [§2] §2 (Hydrodynamic model): The continuum equations are stated without a direct comparison to experimental velocity statistics or interface widths in heterogeneous bacterial suspensions. Because the central claim rests on the model faithfully capturing steric, frictional, and polarity dynamics at the relevant scales, a quantitative benchmark against existing experimental data for even one heterogeneous pattern is needed to support the extrapolation to real suspensions.
Authors: The continuum model is the standard hydrodynamic description used in prior literature on dense bacterial suspensions, with parameters chosen to lie within experimentally reported ranges. Direct quantitative benchmarks against heterogeneous experimental data are not included because such targeted data for quenched patterns remain limited. The work is positioned as a theoretical exploration of realizable heterogeneities rather than a direct experimental validation study. revision: no
Circularity Check
No circularity; results are direct outputs of numerical integration
full rationale
The paper reports outcomes from direct numerical simulations of a standard hydrodynamical model (with quenched activity patterns) rather than any derivation that reduces the target quantities (interface fluctuations, residence times) to fitted parameters or self-referential definitions. No equations are presented that equate a prediction to its own input by construction, and no self-citation chain is invoked as the sole justification for the central claim. The simulation itself constitutes an independent computation whose outputs are not forced by the reported statistics.
Axiom & Free-Parameter Ledger
Lean theorems connected to this paper
-
IndisputableMonolith/Cost/FunctionalEquation.leanwashburn_uniqueness_aczel unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
the mean bacterial velocity field u(x, t) evolves according to the Toner-Tu Swift-Hohenberg (TTSH) model ... ∂tu + λu · ∇u = −∇p − Γ0∇2u − Γ2∇4u − (α + β|u|2)u
-
IndisputableMonolith/Foundation/RealityFromDistinction.leanreality_from_one_distinction unclear?
unclearRelation between the paper passage and the cited Recognition theorem.
quenched activity patterns ... α ≡ α(x) ... circular or rectangularly striped patches of activity αL < 0 surrounded by ... αS > 0
What do these tags mean?
- matches
- The paper's claim is directly supported by a theorem in the formal canon.
- supports
- The theorem supports part of the paper's argument, but the paper may add assumptions or extra steps.
- extends
- The paper goes beyond the formal theorem; the theorem is a base layer rather than the whole result.
- uses
- The paper appears to rely on the theorem as machinery.
- contradicts
- The paper's claim conflicts with a theorem or certificate in the canon.
- unclear
- Pith found a possible connection, but the passage is too broad, indirect, or ambiguous to say the theorem truly supports the claim.
Reference graph
Works this paper leans on
-
[1]
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)
work page 2013
- [2]
-
[3]
I. S. Aranson, Reports on Progress in Physics 85, 076601 (2022)
work page 2022
- [4]
- [5]
- [6]
- [7]
- [8]
-
[9]
T. Sanchez, D. T. Chen, S. J. DeCamp, M. Heymann, and Z. Dogic, Nature 491, 431 (2012)
work page 2012
-
[10]
H. H. Wensink, J. Dunkel, S. Heidenreich, K. Drescher, R. E. Goldstein, H. L¨ owen, and J. M. Yeomans, Pro- ceedings of the national academy of sciences 109, 14308 (2012)
work page 2012
- [11]
-
[12]
A. Doostmohammadi, M. F. Adamer, S. P. Thampi, and J. M. Yeomans, Nature communications 7, 10557 (2016)
work page 2016
-
[13]
A. Doostmohammadi, T. N. Shendruk, K. Thijssen, and J. M. Yeomans, Nature communications 8, 15326 (2017)
work page 2017
-
[14]
C. Chen, S. Liu, X.-q. Shi, H. Chat´ e, and Y. Wu, Nature 542, 210 (2017)
work page 2017
-
[15]
B. Mart´ ınez-Prat, J. Ign´ es-Mullol, J. Casademunt, and F. Sagu´ es, Nature physics15, 362 (2019)
work page 2019
-
[16]
S. Liu, S. Shankar, M. C. Marchetti, and Y. Wu, Nature 590, 80 (2021)
work page 2021
- [17]
- [18]
- [19]
- [20]
-
[21]
S. Mukherjee, R. K. Singh, M. James, and S. S. Ray, Nature Physics 19, 891 (2023)
work page 2023
-
[22]
V. Bratanov, F. Jenko, and E. Frey, Proceedings of the National Academy of Sciences 112, 15048 (2015)
work page 2015
- [23]
- [24]
- [25]
- [26]
-
[27]
S. Mukherjee, R. K. Singh, M. James, and S. S. Ray, Physical Review Letters 127, 118001 (2021)
work page 2021
- [28]
-
[29]
R. K. Singh, S. Mukherjee, and S. S. Ray, Physical Re- view Fluids 7, 033101 (2022)
work page 2022
-
[30]
H. Wioland, E. Lushi, and R. E. Goldstein, New Journal of Physics 18, 075002 (2016)
work page 2016
-
[31]
K.-T. Wu, J. B. Hishamunda, D. T. Chen, S. J. De- Camp, Y.-W. Chang, A. Fern´ andez-Nieves, S. Fraden, and Z. Dogic, Science 355, eaal1979 (2017)
work page 2017
-
[32]
S. Chandragiri, A. Doostmohammadi, J. M. Yeomans, and S. P. Thampi, Physical Review Letters 125, 148002 (2020)
work page 2020
-
[33]
D. Nishiguchi, S. Shiratani, K. A. Takeuchi, and I. S. Aranson, arXiv preprint arXiv:2407.05269 (2024)
-
[34]
S. P. Thampi, R. Golestanian, and J. M. Yeomans, Phys- ical Review E 90, 062307 (2014)
work page 2014
-
[35]
W. D. Hoff, M. A. van der Horst, C. B. Nudel, and K. J. Hellingwerf, Chemotaxis: Methods and Protocols pp. 25– 49 (2009)
work page 2009
- [36]
-
[37]
J. Yang, P. E. Arratia, A. E. Patteson, and A. Gopinath, Journal of the Royal Society Interface 16, 20180960 (2019)
work page 2019
-
[38]
R. Thar and M. K¨ uhl, Proceedings of the National Academy of Sciences 100, 5748 (2003)
work page 2003
- [39]
-
[40]
A. P. Petroff, X.-L. Wu, and A. Libchaber, Physical re- 10 view letters 114, 158102 (2015)
work page 2015
- [41]
-
[42]
I. Engelhardt, D. Patko, Y. Liu, M. Mimault, G. de Las Heras Martinez, T. George, M. MacDonald, M. Ptash- nyk, T. Sukhodub, N. Stanley-Wall, et al., The ISME Journal 16, 2337 (2022)
work page 2022
-
[43]
D. Nishiguchi, I. S. Aranson, A. Snezhko, and A. Sokolov, Nature communications 9, 4486 (2018)
work page 2018
-
[44]
H. Reinken, D. Nishiguchi, S. Heidenreich, A. Sokolov, M. B¨ ar, S. H. Klapp, and I. S. Aranson, Communications Physics 3, 76 (2020)
work page 2020
-
[45]
H. Reinken, S. Heidenreich, M. B¨ ar, and S. H. Klapp, Physical Review Letters 128, 048004 (2022)
work page 2022
- [46]
-
[47]
J. Palacci, S. Sacanna, A. P. Steinberg, D. J. Pine, and P. M. Chaikin, Science 339, 936 (2013)
work page 2013
-
[48]
M. Schuppler, F. C. Keber, M. Kr¨ oger, and A. R. Bausch, Nature communications 7, 13120 (2016)
work page 2016
-
[49]
T. D. Ross, H. J. Lee, Z. Qu, R. A. Banks, R. Phillips, and M. Thomson, Nature 572, 224 (2019)
work page 2019
- [50]
- [51]
-
[52]
S. Shankar, A. Souslov, M. J. Bowick, M. C. Marchetti, and V. Vitelli, Nature Reviews Physics 4, 380 (2022)
work page 2022
-
[53]
J. Arlt, V. A. Martinez, A. Dawson, T. Pilizota, and W. C. Poon, Nature communications 9, 768 (2018)
work page 2018
-
[54]
G. Frangipane, D. Dell’Arciprete, S. Petracchini, C. Maggi, F. Saglimbeni, S. Bianchi, G. Vizsnyiczai, M. L. Bernardini, and R. Di Leonardo, Elife 7, e36608 (2018)
work page 2018
-
[55]
J. Arlt, V. A. Martinez, A. Dawson, T. Pilizota, and W. C. Poon, Nature communications 10, 2321 (2019)
work page 2019
-
[56]
See animation showing the development of the vor- ticity field for an activity quench over a circu- lar geometry https://www.youtube.com/watch?v=VJJg_ SMyfw8, along with a composite comparison of the emerging interface https://www.youtube.com/shorts/ a1vjBuiQSLM. Made with Processing [91–93]
-
[57]
See https://youtube.com/shorts/AIcsM2GP3d4 for an animation showing the formation of a large vortex, simi- lar to what is obtained for highly active turbulence under circular confinement
-
[58]
L. Puggioni, G. Boffetta, and S. Musacchio, Physical Re- view E 106, 055103 (2022)
work page 2022
-
[59]
See https://www.youtube.com/shorts/l6AVRfRtIyo for an animation showing the evolution of the interface sep- arating the highly active flow from a region low activity
- [60]
-
[61]
N. Cao, S. Chen, and G. D. Doolen, Physics of Fluids 11, 2235 (1999)
work page 1999
-
[62]
G. Boffetta and R. E. Ecke, Annual review of fluid me- chanics 44, 427 (2012)
work page 2012
-
[63]
Frisch, Turbulence: the legacy of AN Kolmogorov (Cambridge university press, 1995)
U. Frisch, Turbulence: the legacy of AN Kolmogorov (Cambridge university press, 1995)
work page 1995
-
[64]
C. Meneveau and K. Sreenivasan, Physical review letters 59, 1424 (1987)
work page 1987
-
[65]
C. Meneveau and K. Sreenivasan, Journal of Fluid Me- chanics 224, 429 (1991)
work page 1991
-
[66]
S. Mukherjee, S. D. Murugan, R. Mukherjee, and S. S. Ray, Phys. Rev. Lett. 132, 184002 (2024)
work page 2024
-
[67]
P. Grassberger and I. Procaccia, Physical review letters 50, 346 (1983)
work page 1983
-
[68]
See https://www.youtube.com/shorts/AzAhfmAAwyc for an animation showing the evolution of the vorticity field for a striped activity pattern
-
[69]
N. Pal, P. Perlekar, A. Gupta, and R. Pandit, Physical Review E 93, 063115 (2016)
work page 2016
-
[70]
N. B. Padhan and R. Pandit, Physical Review Research 5, L032013 (2023)
work page 2023
-
[71]
P. Perlekar, S. S. Ray, D. Mitra, and R. Pandit, Phys. Rev. Lett. 106, 054501 (2011)
work page 2011
-
[72]
A. Bhatnagar, A. Gupta, D. Mitra, R. Pandit, and P. Perlekar, Phys. Rev. E 94, 053119 (2016)
work page 2016
- [73]
-
[74]
L. Puggioni, G. Boffetta, and S. Musacchio, Physical Re- view E 107, 055107 (2023)
work page 2023
-
[75]
Alexakis, Journal of Fluid Mechanics 977, R1 (2023)
A. Alexakis, Journal of Fluid Mechanics 977, R1 (2023)
work page 2023
-
[76]
T. Heus and H. J. Jonker, Journal of the Atmospheric Sciences 65, 1003 (2008)
work page 2008
-
[77]
J. Westerweel, C. Fukushima, J. M. Pedersen, and J. C. Hunt, Journal of Fluid Mechanics 631, 199 (2009)
work page 2009
-
[78]
K. Chauhan, J. Philip, C. M. De Silva, N. Hutchins, and I. Marusic, Journal of Fluid Mechanics 742, 119 (2014)
work page 2014
-
[79]
T. Watanabe, Y. Sakai, K. Nagata, Y. Ito, and T. Hayase, Physics of Fluids 27 (2015)
work page 2015
-
[80]
G. Borrell and J. Jim´ enez, Journal of Fluid Mechanics 801, 554 (2016)
work page 2016
discussion (0)
Sign in with ORCID, Apple, or X to comment. Anyone can read and Pith papers without signing in.