Density-dependent transport coefficients in two-dimensional cellular aggregates

Subhadip Chakraborti; Vasily Zaburdaev

REVIEW 4 major objections 5 minor 68 references

Density-dependent transport coefficients in two-dimensional cellular aggregates

T0 review · 4 major / 5 minor · reviewed 2026-08-06 · deepseek-v4-flash

Challenge this review Re-run · record.json Download PDF Read on arXiv ↗ Review JSON

Pith's one-line read Two analytical formulas fix the density-dependent bulk diffusivity and conductivity of pili-driven cellular aggregates, and predict that transport slows as colonies form.

desk verdict Solid incremental extension with a real soft spot: the row-reduction to 1D is an assumption that simulation support makes plausible but not rigorous; worth refereeing. read the letter →

arxiv 2507.19919 v1 pith:Z3Y5EEX5 submitted 2025-07-26 cond-mat.stat-mech cond-mat.softphysics.bio-ph

Subhadip Chakraborti , Vasily Zaburdaev This is my paper · ORCID

classification cond-mat.stat-mechcond-mat.softphysics.bio-ph MSC 82C2282C3182C7060K35

keywords two-dimensionalcellularaggregatespili-mediatedaggregationfluctuatinghydrodynamicsbulkdiffusivityconductivityexclusionprocesszero-rangetransportcoefficients

verification ladder T0 review T1 audit T2 compute T3 formal

The pith

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

The reading

Starting from a two-dimensional lattice model of Neisseria gonorrhoeae aggregation, in which hard-core cells hop by substrate-driven diffusion and by short-lived attractive pili pairs, the paper derives closed-form expressions for two macroscopic transport coefficients: bulk diffusivity $D(\rho)$ and conductivity $\chi(\rho)$. These coefficients are the proportionality factors of the hydrodynamic cell current in a density gradient and under an external force, respectively. The formulas show that both coefficients decrease as the cell-cell interaction rate $p$ increases, predicting that transport slows down as colonies form. Monte Carlo simulations of the 2D exclusion model match the analytical curves well for moderate activity, with deviations at the highest activity $p = 0.9$. If the derivation holds, the density dependence of transport in such cellular aggregates is fixed by the microscopic pili parameters rather than left to phenomenological modelling.

What carries the argument

The load-bearing object is the local coarse-grained observable $$G(\rho) = q(1-\rho) + p\rho \sum_{m=0}^\infty m\, v(m)\, P(m|\rho),$$ from which the transport coefficients follow as $D(\rho) = -\frac{1}{4}\partial G/\partial \rho$ and $\chi(\rho) = \frac{1}{4}\rho G(\rho)$. All microscopic detail enters through the jump weight $v(m) = \frac{s_1}{2}(1 + s_2 m/l_0)\exp(-m/l_0)$, the probability that two exponentially distributed pili facing each other across a gap $m$ form a bound pair, with $s_1$ and $s_2$ fitted to simulation. To average $v(m)$ analytically, the exclusion model is mapped to an unbounded mass transfer model (a zero-range process) in which site masses are the gaps between particles and mass-mass correlations vanish beyond nearest neighbours, giving a factorised steady state $P(m) = f(m) z^m / F(z)$. The mapping converts an interacting 2D exclusion problem into a tractable product-measure calculation, and the reverse mapping returns $D$ and $\chi$ to the original cell density $\rho$.

What would settle it

In the 2D agent-based model, measure the bulk diffusivity from the decay rate of a small sinusoidal density perturbation at $p = 0.9$ while increasing system size $L$; the paper already reports a mismatch at that activity, so if the discrepancy does not shrink as $L$ grows, the 1D row-decoupling and factorised steady-state assumptions are the point of failure.

Watch

Extended reading notes

Core claim

The paper's central result is that the fluctuating hydrodynamics of the pili-driven aggregation model is governed by the two transport coefficients $$D(\rho) = \frac{1}{4}\left[q - p \frac{\partial}{\partial \rho}\left(\rho \sum_{m=0}^\infty m\, v(m)\, P(m|\rho)\right)\right]$$ and $$\chi(\rho) = \frac{1}{4}\left[q\rho(1-\rho) + p\$rho^{2}$ \sum_{m=0}^\infty m\, v(m)\, P(m|\rho)\right],$$ where $v(m)$ is the probability that two pili facing each other across a gap $m$ form a bound pair and $P(m|\rho)$ is the steady-state mass distribution of an equivalent unbounded model. The authors derive the observable $G(\rho)$ analytically, avoiding the numerical averages of their earlier 1D study, and use macroscopic fluctuation theory to show that $D$ and $\chi$ enter the hydrodynamic current as $-D\nabla\rho$ and $\chi F$. Both coefficients decrease with increasing $p$, which the authors interpret as transport slowdown during colony formation; the equilibrium limit $p = 0$ recovers $D = 1/4$ and $\chi = \rho(1-\rho)/4$ of the symmetric exclusion process.

Load-bearing premise

The calculation assumes that each row and column of the two-dimensional lattice is an independent one-dimensional exclusion process, so that horizontal and vertical currents never couple; if that decoupling fails, the two formulas do not describe the real 2D colony.

Editorial extensions

If this is right

If Eqs. (30) and (31) are correct, the hydrodynamic current in a colony is quantitatively fixed once the density and the microscopic parameters ($p$, $N_p$, $l_0$, $T_0$) are known.
Both $D(\rho)$ and $\chi(\rho)$ fall as the cell-cell interaction rate $p$ increases, so the model predicts that colony formation itself slows further transport, a feedback that could stabilise aggregates.
At $p = 0$ the formulas reduce to the known symmetric exclusion process values $D = 1/4$ and $\chi = \rho(1-\rho)/4$, connecting the active model to its equilibrium baseline.
Conductivity computed from current fluctuations in the unbiased system and from the mean current under a small bias agree with each other and with the analytical formula, supporting the macroscopic fluctuation theory description.
The same fluctuating-hydrodynamics route should apply to other contractile-force aggregates, such as tumour spheroids and organoids.

Reading between the lines

Editorial extensions of the paper, not claims the author makes directly.

Beyond the paper: if vertical-horizontal current coupling becomes measurable at high $p$ in larger simulations, the effective 2D conductivity could fall below the 1D-reduction prediction, since the decoupling assumption is exactly what the paper identifies as strained there.
Beyond the paper: the fitted parameters $s_1$ and $s_2$, which absorb pili lifetime and number, could be estimated directly from pili statistics in experiments; then Eqs. (30)-(31) become parameter-free predictions for colony-scale transport.
Beyond the paper: one could test a density-dependent Einstein-type relation between $D(\rho)$ and $\chi(\rho)$; the paper does not claim it, but its formulas supply both coefficients from the same observable $G(\rho)$.

Share X Bluesky LinkedIn Reddit HN

Signed reviews

No signed human review yet.

Editorial analysis

A structured set of objections, weighed in public.

Desk editor's note, referee report, and a circularity audit.

Desk Editor's Note

read the letter

This is a sensible, incremental step beyond the authors' 1D work, and the main results are the explicit formulas for D(rho) and chi(rho), Eqs. (30)-(31). What the paper does well: it works through the mass-dependent jump weights, proposes a factorized steady state in the UMTM, expresses the transport coefficients in terms of that steady state, and tests D and chi against simulations for p up to 0.9. The agreement is genuinely good for an interacting active lattice gas, and the paper is honest about the high-p mismatch and about the numerical evaluation of P(m).

The soft spot is the 2D-to-1D reduction. The authors argue that x and y currents are independent, so one row (or column) represents the whole system. That is not the same as showing a row is a closed 1D system. In the 2D exclusion model, vertical hops change row occupancy, so the UMTM ring length is time-dependent; the appendix proves x/y separability only in the p=0 SSEP limit. For p>0 the argument rests on the absence of cross-linking between horizontal and vertical pili, which is plausible but not a derivation. So Eqs. (30)-(31) strictly describe an auxiliary 1D chain model. The simulation agreement then gives empirical support for the reduction, but it is support from the same model, and the D and chi comparison is indirect. I would like to see a direct test: measure row occupancy fluctuations and compare against the closed-ring assumption, or probe the 1D row prediction with a perturbation that is not uniform in the transverse direction.

The two fitting constants s1 and s2 are a minor weakness. They are calibrated to microscopic bound-pair counts rather than to the transport coefficients, so not circular, but they temper the claim of a first-principles derivation.

Overall, the central claim holds up at moderate activity and the paper is worth a serious referee. The reduction issue can likely be fixed with additional numerical diagnostics, so my recommendation is to send it to review, with the authors expected to address the open-row problem.

Referee Report

4 major / 5 minor

Summary. This paper develops a fluctuating hydrodynamics description for a two-dimensional lattice model of cellular aggregation mediated by intermittent dipolar (pili) forces. The authors reduce the 2D exclusion model to independent 1D rows/columns, map each row to an unbounded mass transfer model (UMTM), and assume a factorized steady state to derive analytic expressions for the bulk diffusivity D(ρ) and conductivity χ(ρ) [Eqs. (30)-(31)]. They then compare these expressions with Monte Carlo simulations of the full 2D model, obtaining D from the exponential relaxation of a sinusoidal density perturbation and χ from both the variance of the fluctuating current and the mean current under a small external bias. Agreement is good for p up to 0.9, and both coefficients decrease with increasing cell-cell interaction strength p, which the authors interpret as transport slowdown during colony formation.

Significance. If the derivation holds, this is a valuable contribution: analytic, density-dependent transport coefficients for an interacting active lattice gas are rare, and the paper's closed forms explicitly connect D and χ to the microscopic parameters (p, Np, l0, T0). The numerical verification is unusually complete—two independent measurements of χ and one of D, all consistent with the theory across densities. The main caveat is that the central 2D-to-1D reduction is not rigorously derived for p>0; the appendix proves separability only for the p=0 SSEP limit. Thus the significance is conditional on the validity, or further justification, of that reduction.

major comments (4)

[Sec. III (before III A) and Appendix A] The reduction of the 2D model to independent 1D rows/columns is not established for p>0. Appendix A derives x/y independence of the current only for the p=0 SSEP limit [Eqs. (38)-(41)]. For p>0, the argument based on isotropy and absence of cross-linking between horizontal and vertical pili does not imply that vertical particle exchange decouples from horizontal gap statistics: a vertical hop changes the occupancy of a row and therefore changes the inter-particle gaps m_i that control the interaction weight v(m_i). In the UMTM construction of Sec. III A, the ring length is the instantaneous row particle number N1, which fluctuates because rows are open systems; the 1D master-equation derivation of D and χ assumes a closed ring with fixed particle number. Consequently Eqs. (30)-(31) are derived for an auxiliary 1D chain, and their validity for the coupled 2D system rests on the simulation agreement in Fig. 5, not on a derivation. I recommend either a controlled derivation of the decoupling (e.g., showing that the row dynamics in the 2D model converge to the 1D UMTM in the hydrodynamic limit) or an independent numerical test of the reduction, such as comparing row-resolved transport coefficients with full 2D results, or measuring correlations between horizontal currents and vertical occupancy fluctuations.
[Sec. III B] The assertion that the microscopic jump weight v(m) is independent of ρ and p is stated without supporting data: 'We further check (not shown) that v(m) does not depend on density ρ and cell-cell interaction rate p.' Since v(m) enters every transport coefficient via Eqs. (30)-(31), a hidden p-dependence would mean that part of the p-dependence of D and χ is fitted rather than predicted. Please either show the v(m) data for multiple ρ and p values, or provide a microscopic argument for this independence.
[Appendix B vs. Sec. III C] The rates in Appendix B, Eq. (42), are inconsistent with Eq. (7) in the main text. Eq. (7) assigns a total loss rate (q a_i + p a_i v(m_i))/2 (with each directional rate /4), leading to D = (1/4)∂G/∂ρ and χ = (1/4)G. Appendix B uses q a_i for the q-loss and p a_i v(m_i)/2 for the p-loss, yielding D = (1/2)∂G/∂ρ and χ = (1/2)G in Eq. (50). The main-text factors are the ones that agree with the p=0 SSEP limit (D=1/4) and with simulation; the appendix should be corrected to match the model definition, so the two derivations do not give different answers.
[Sec. III E, Eqs. (20)-(21)] The second-moment equation and its steady-state reduction contain apparent typographical errors: Eq. (20) as printed does not follow from the rates in Eq. (7), and Eq. (21) repeats the term '-2⟨m_i g_i⟩ + ⟨g_i⟩' twice. The final relation Eq. (22) is correct, but the intermediate equations should be fixed for the derivation to be reproducible.

minor comments (5)

[Sec. III A] The notation N1 is used without definition; please define it as the number of particles in the selected row.
[Fig. 4(c) caption] The statement 'correlation vanishes beyond distance r = 0 for any p' is confusing; it should read 'for r > 0'.
[Sec. III E] The statement 'we can now write ⟨BiCj⟩ = ⟨Bi⟩⟨Cj⟩ if i≠j' is justified only by the observed absence of long-range mass correlations; this is an assumption rather than a proven property of the UMTM for p>0, and should be flagged as such in the text.
[Ref. [54]] The footnote [54] 'Please see the appendices for detail' should cite the specific appendix sections (VI A and VI B) rather than a bare reference.
[Sec. II] The definition of MCT as 'each particle, on average, has one opportunity to jump' should be clarified: with the 1/4 hopping weights, the mean number of horizontal jumps per MCT is 1/2 per particle, which is relevant for the relationship between simulation time and the hydrodynamic time τ.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: the only fitted input is the microscopic bound-pair weight v(m), while the claimed predictions D(ρ) and χ(ρ) are derived from it and verified by independent simulations.

full rationale

The derivation chain is not circular in the sense defined by the review criteria. The only empirically fitted quantities are the two constants s1=0.4 and s2=0.8 entering the microscopic bound-pair weight v(m) in Eq. (6), which the paper transparently calibrates against the measured scaled mean number of bound pili pairs (Fig. 4a). This is a microscopic input, not the predicted output: the transport coefficients D(ρ) and χ(ρ) in Eqs. (30)-(31) are then obtained from a master-equation/MFT calculation (Secs. III C-D), and the comparison in Fig. 5 uses independent measurements—sinusoidal-profile relaxation for D and biased/unbiased current fluctuations for χ—that were not used to set s1 and s2. The self-citations to Ref. [36] occur for the standard EM-to-UMTM mapping and as background; the mapping itself is also cited to the external zero-range-process literature [40], so the self-citation is not load-bearing. The unproven step is the 2D-to-1D row/column reduction (Sec. III), which the paper states as an assumption supported by isotropy and absence of pilus cross-linking, and whose breakdown at high p is explicitly acknowledged in Sec. IV. That is a correctness/validity risk, not a circular reduction: no equation is defined in terms of its own prediction, and no fitted transport coefficient is relabeled as a prediction. The paper is self-contained against its own simulations as external benchmarks, so the honest finding is no significant circularity.

Assumptions & free parameters 2 free parameters · 6 assumptions · 0 invented entities

The derivation rests on the UMTM mapping, factorization, local steady state, row and column decoupling, and the fitted interaction kernel v(m). No genuinely new physical entities are introduced; the UMTM is a mathematical mapping of vacancies to masses, not a new force or particle.

free parameters (2)

s1 = 0.4 for l0=2, T0=10, Np=20
Fitted to the scaled mean number of bound pili pairs v(m) from 2D EM simulations in Eq. (6) and Fig. 4(a). It absorbs the influence of Np, T0, and the pili binding rate, but is treated as a constant for all densities and activities.
s2 = 0.8 for l0=2 (0.85 for l0=4, 0.9 for l0=8)
Fitted alongside s1 to reproduce the distance dependence of v(m) in Eq. (6) and Fig. 4(a). The paper sets s2=0.8 for all subsequent calculations.

assumptions (6)

domain assumption Factorization of the steady-state measure in the UMTM, P({mi}) = Z^{-1} prod_i f(mi), with f(m)=1/prod_{k=1}^m u(k).
Assumed in Sec. III.E, Eq. (15), based on vanishing correlations in Fig. 4(c). This is standard for zero-range processes under certain conditions, but its validity for p>0 in this non-equilibrium model is verified numerically, not proven.
domain assumption Local steady state assumption: local observables can be evaluated using the global steady-state mass distribution at the local density.
Invoked in Sec. III.E to compute G(rho) from single-site averages. It is a common hydrodynamic closure but breaks down in strongly clustered states, as the paper acknowledges for high p.
ad hoc to paper Row and column decoupling in 2D: the hydrodynamics of the full 2D system equals the 1D UMTM result for a representative row or column.
Stated in Sec. III without a derivation: 'we can derive hydrodynamics for any row or column of the 2D system which effectively will represent the hydrodynamics of the entire system.' This is the weakest load-bearing premise; it is supported only by comparison with 2D simulations.
domain assumption Pili lengths are exponentially distributed with mean l0 and pili lifetimes are exponentially distributed with mean T0.
Used in Sec. II and in Eqs. (2)-(5) to compute the pairing probability that defines v(m). This is a biologically motivated simplification, not derived from data in this paper.
ad hoc to paper The interaction kernel v(m) does not depend on global density rho or cell-cell interaction rate p.
Stated in Sec. III.B with 'We further check (not shown)'. This allows s1 and s2 to be fixed once for all rho and p, but the supporting data are not displayed.
standard math A small external bias F modifies jump rates through a local detailed balance form, expanded to first order in F.
Used in Eq. (12) and in Appendix B to derive the biased hydrodynamic equation and identify chi(rho). This is the standard MFT weak-bias assumption for gradient-type processes.

how reviews work

0 comments

Cite this review

Pith. "Pith review of Density-dependent transport coefficients in two-dimensional cellular aggregates." pith.science (2026). https://pith.science/paper/Z3Y5EEX5

@misc{pith2026250719919,
  author       = {Pith},
  title        = {Pith review of: Density-dependent transport coefficients in two-dimensional cellular aggregates},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/Z3Y5EEX5}},
  note         = {Machine review of arXiv:2507.19919}
}

read the original abstract

The large-scale collective behavior of biological systems can be characterized by macroscopic transport, which arises from the non-equilibrium microscopic interactions among individual constituents. A prominent example is the formation of dynamic aggregates by motile eukaryotic cells or bacteria mediated by active contractile forces. In this work, we develop the two-dimensional fluctuating hydrodynamics theory based on the microscopic dynamics of a model system of aggregation by \textit{Neisseria gonorrhoeae} bacteria. The derivation of two macroscopic transport coefficients of bulk diffusivity and conductivity which determine hydrodynamic current of cells is the central result of this work. By showing how transport coefficients depend on cell density and microscopic parameters of the system we predict transport slowdown during the colony formation process. This study provides valuable analytical tools for quantifying hydrodynamic transport in experimental systems involving cellular aggregation occurring due to intermittent contractile dipole forces.

Figures

Figures reproduced from arXiv: 2507.19919 by the authors.

**Figure 2.** FIG. 2. Clustering transition. We consider [PITH_FULL_IMAGE:figures/full_fig_p003_2.png] view at source ↗

**Figure 3.** FIG. 3. Characterization of clustering transition. Panel (a): [PITH_FULL_IMAGE:figures/full_fig_p004_3.png] view at source ↗

Figures from the paper (2 more)

**Figure 4.** Figure 4: FIG. 4. Mass-dependent jump weights and distributions in [PITH_FULL_IMAGE:figures/full_fig_p006_4.png]

**Figure 5.** Figure 5: FIG. 5. Determination and demonstration of transport coeffi [PITH_FULL_IMAGE:figures/full_fig_p009_5.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

68 extracted references · 62 canonical work pages

[1]

Ben-Jacob, I

E. Ben-Jacob, I. Cohen, and D. L. Gutnick, Coopera- tive organization of bacterial colonies: From genotype to morphotype, Annu. Rev. Microbiol. 52, 779 (1998)

work page 1998
[2]

Taktikos, Y

J. Taktikos, Y. T. Lin, H. Stark, N. Biais, and V. Zabur- daev, Pili-induced clustering of n. gonorrhoeae bacteria, PLoS ONE 10, e0137661 (2015)

work page 2015
[3]

Mueller-Klieser, Multicellular spheroids: A review on cellular aggregates in cancer research, J

W. Mueller-Klieser, Multicellular spheroids: A review on cellular aggregates in cancer research, J. Cancer. Res. Clin. Oncol. 113, 101 (1987)

work page 1987
[4]

q − p ∂ ∂ρ ( ρ ∞X m=0 mv(m)P (m | ρ) )# , (30) χ(ρ) = 1 4

= 0 .85 and s2(l0 = 8) = 0 .9. We further check (not shown) that v(m) does not depend on density ρ and cell- cell interaction rate p. Throughout this work we consider T0 = 10 MCT, Np = 20 and l0 = 2, therefore, we set s1 = 0.4 and s2 = 0.8 for the rest of this paper. The expression for jump weight v(m) given by Eq. (6) is the key starting point allowing t...

work page
[5]

Friedl and K

P. Friedl and K. Wolf, Tumour-cell invasion and migra- tion: diversity and escape mechanisms, Nat. Rev. Cancer 3, 362 (2003)

work page 2003
[6]

T. M. Koch, S. Munster, N. Bonakdar, J. P. Butler, and B. Fabry, 3d traction forces in cancer cell invasion, PLoS ONE 7, e33476 (2012)

work page 2012
[7]

Mark et al., Collective forces of tumor spheroids in three-dimensional biopolymer networks, eLife 9, e51912 (2020)

C. Mark et al., Collective forces of tumor spheroids in three-dimensional biopolymer networks, eLife 9, e51912 (2020)

work page 2020
[8]

Karow et al., Direct pericyte-to-neuron reprogram- ming via unfolding of a neural stem cell-like program, Nat

M. Karow et al., Direct pericyte-to-neuron reprogram- ming via unfolding of a neural stem cell-like program, Nat. Neurosci. 21, 932 (2018)

work page 2018

Show all 68 references

[9]

S. Falk, D. Han, and M. Karow, Cellular identity through the lens of direct lineage reprogramming, Curr. Opin. Genet. Dev. 70, 97 (2021)

2021
[10]

Lancaster et al., Cerebral organoids model human brain development and microcephaly, Nature 501, 373 (2013)

M. Lancaster et al., Cerebral organoids model human brain development and microcephaly, Nature 501, 373 (2013)

2013
[11]

S. P. Pasca, The rise of three-dimensional human brain cultures, Nature 553, 437 (2018)

2018
[12]

P¨ onisch, C

W. P¨ onisch, C. A. Weber, G. Juckeland, N. Biais, and V. Zaburdaev, Multiscale modeling of bacterial colonies: how pili mediate the dynamics of single cells and cellular aggregates, New J. Phys. 19, 015003 (2017)

2017
[13]

Fern´ andez et al

P. Fern´ andez et al. , Surface-tension-induced bud- ding drives alveologenesis in human mammary gland organoids, Nat. Phys. 17, 1130 (2021)

2021
[14]

M. M. Onesto, J. il Kim, and S. P. Pasca, Assembloid models of cell-cell interaction to study tissue and disease biology, Cell Stem Cell 31, 1563 (2024)

2024
[15]

Wang et al., Molecular and cellular dynamics of the de- veloping human neocortex, Nature 10.1038/s41586-024- 08351-7 (2025)

L. Wang et al., Molecular and cellular dynamics of the de- veloping human neocortex, Nature 10.1038/s41586-024- 08351-7 (2025)

2025 doi
[16]

Sanchez et al., Spontaneous motion in hierarchically assembled active matter, Nature 491, 431 (2012)

T. Sanchez et al., Spontaneous motion in hierarchically assembled active matter, Nature 491, 431 (2012)

2012
[17]

Alert and X

R. Alert and X. Trepat, Physical models of collective cell migration, Annu. Rev. Condens. Matter Phys. 11, 77 (2020)

2020
[18]

Onsager, Reciprocal relations in irreversible processes

L. Onsager, Reciprocal relations in irreversible processes. ii., Phys. Rev. 38, 2265 (1931)

1931
[19]

R. K. Pathria and P. D. Beale, Statistical Mechanics, 3rd ed. (Elsevier, Amsterdam, 2011)

2011
[20]

Bertini, A

L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and C. Landim, Macroscopic fluctuation theory, Rev. Mod. Phys. 87, 593 (2015)

2015
[21]

C. W. Gardiner, Handbook of Stochastic Methods: For Physics, Chemistry and the Natural Sciences(Springer Berlin Heidelberg, 1990)

1990
[22]

F. S. Gnesotto1, F. Mura1, J. Gladrow, and C. P. Broed- ersz, Broken detailed balance and non-equilibrium dy- namics in living systems: a review, Rep. Prog. Phys. 81, 066601 (2018)

2018
[23]

Thapa et al., Nonequilibrium probability currents in optically-driven colloidal suspensions, SciPost Phys

S. Thapa et al., Nonequilibrium probability currents in optically-driven colloidal suspensions, SciPost Phys. 17, 096 (2024)

2024
[24]

Ascione et al., Gradient-induced instability in tumour spheroids unveils the impact of microenvironmental nu- trient changes, Sci

F. Ascione et al., Gradient-induced instability in tumour spheroids unveils the impact of microenvironmental nu- trient changes, Sci. Rep. 14, 20837 (2024)

2024
[25]

Xu et al., Self-organization of swimmers drives long- range fluid transport in bacterial colonies, Nat

H. Xu et al., Self-organization of swimmers drives long- range fluid transport in bacterial colonies, Nat. Commun. 10, 1792 (2019)

2019
[26]

P¨ onisch, K

W. P¨ onisch, K. B. Eckenrode, K. Alzurqa, H. Nasrollahi, C. Weber, V. Zaburdaev, and N. Biais, Pili mediated intercellular forces shape heterogeneous bacterial micro- colonies prior to multicellular differentiation, Sci. Rep. 8, 16567 (2018)

2018
[27]

Tortorella, C

I. Tortorella, C. Argentati, C. Emiliani, S. Martino, and F. Morena, The role of physical cues in the develop- ment of stem cell-derived organoids, Eur. Biophys. J. 51, 105–117 (2022)

2022
[28]

W. I. M. B. Cronenberg T, Hennes M, Antibiotics modu- late attractive interactions in bacterial colonies affecting survivability under combined treatment, PLoS Pathog. 17, e1009251 (2021)

2021
[29]

Craig, M

L. Craig, M. E. Pique, and J. A. Tainer, Type iv pilus structure and bacterial pathogenicity, Nat. Rev. Micro- 12 biol. 2, 363 (2004)

2004
[30]

J. S. Mattick, Type iv pili and twitching motility, Annu. Rev. Microbiol. 56, 289 (2002)

2002
[31]

C. Holz, D. Opitz, L. Greune, R. Kurre, M. Koomey, M. A. Schmidt, and B. Maier, Multiple pilus motors co- operate for persistent bacterial movement in two dimen- sions, Phys. Rev. Lett. 104, 178104 (2010)

2010
[32]

A. Merz, M. So, and M. Sheetz, Pilus retraction powers bacterial twitching motility, Nature 407, 98 (2000)

2000
[33]

J. M. Skerker and H. C. Berg, Direct observation of ex- tension and retraction of type iv pili, Proc. Natl. Acad. Sci. U.S.A. 98, 6901 (2001)

2001
[34]

Craig, K

L. Craig, K. T. Forest, and B. Maier, Type iv pili: dynamics, biophysics and functional consequences, Nat. Rev. Microbiol. 17, 429 (2019)

2019
[35]

Bertini, A

L. Bertini, A. De Sole, D. Gabrielli, G. Jona-Lasinio, and C. Landim, Fluctuations in stationary nonequilib- rium states of irreversible processes, Phys. Rev. Lett. 87, 040601 (2001)

2001
[36]

Derrida, Non-equilibrium steady states: fluctuations and large deviations of the density and of the current, J

B. Derrida, Non-equilibrium steady states: fluctuations and large deviations of the density and of the current, J. Stat. Mech. 2007, P07023 (2007)

2007
[37]

Chakraborti and V

S. Chakraborti and V. Zaburdaev, Transport in cellu- lar aggregates described by fluctuating hydrodynamics, Phys. Rev. Res. 6, 043064 (2024)

2024
[38]

M. E. Cates and J. Tailleur, Motility-induced phase sepa- ration, Annu. Rev. Condens. Matter Phys. 6, 219 (2015)

2015
[39]

C. K. Ellison, G. B. Whitfield, and Y. V. Brun, Type IV Pili: dynamic bacterial nanomachines, FEMS Microbiol- ogy Reviews 46, fuab053 (2021)

2021
[40]

Henrichsen, Twitching motility, Annual Review of Mi- crobiology 37, 81 (1983), pMID: 6139059

J. Henrichsen, Twitching motility, Annual Review of Mi- crobiology 37, 81 (1983), pMID: 6139059

1983
[41]

M. R. Evans and T. Hanney, Nonequilibrium statistical mechanics of the zero-range process and related models, J. Phys. A: Math. Gen. 38, R195 (2005)

2005
[42]

C. A. Weber, Y. T. Lin, N. Biais, and V. Zaburdaev, For- mation and dissolution of bacterial colonies, Phys. Rev. E 92, 032704 (2015)

2015
[43]

H.-S. Kuan, W. P¨ onisch, F. J¨ ulicher, and V. Zaburdaev, Continuum theory of active phase separation in cellular aggregates, Phys. Rev. Lett. 126, 018102 (2021)

2021
[44]

M. E. Cates and J. Tailleur, When are active brown- ian particles and run-and-tumble particles equivalent? consequences for motility-induced phase separation, EPL 101, 20010 (2013)

2013
[45]

Caprini, U

L. Caprini, U. Marini Bettolo Marconi, and A. Puglisi, Spontaneous velocity alignment in motility-induced phase separation, Phys. Rev. Lett. 124, 078001 (2020)

2020
[46]

C. B. Caporusso, P. Digregorio, D. Levis, L. F. Cuglian- dolo, and G. Gonnella, Motility-induced microphase and macrophase separation in a two-dimensional active brow- nian particle system, Phys. Rev. Lett. 125, 178004 (2020)

2020
[47]

Dandekar, S

R. Dandekar, S. Chakraborti, and R. Rajesh, Hard core run and tumble particles on a one-dimensional lattice, Phys. Rev. E 102, 062111 (2020)

2020
[48]

G. S. Redner, A. Baskaran, and M. F. Hagan, Reentrant phase behavior in active colloids with attraction, Phys. Rev. E 88, 012305 (2013)

2013
[49]

Caprini and H

L. Caprini and H. L¨ owen, Flocking without alignment interactions in attractive active brownian particles, Phys. Rev. Lett. 130, 148202 (2023)

2023
[50]

Alston, A

H. Alston, A. O. Parry, R. Voituriez, and T. Bertrand, Intermittent attractive interactions lead to microphase separation in nonmotile active matter, Phys. Rev. E106, 034603 (2022)

2022
[51]

S. N. Majumdar, S. Krishnamurthy, and M. Barma, Nonequilibrium phase transitions in models of aggrega- tion, adsorption, and dissociation, Phys. Rev. Lett. 81, 3691 (1998)

1998
[52]

Chakraborti, T

S. Chakraborti, T. Chakraborty, A. Das, R. Dandekar, and P. Pradhan, Transport and fluctuations in mass ag- gregation processes: Mobility-driven clustering, Phys. Rev. E 103, 042133 (2021)

2021
[53]

Bodineau and B

T. Bodineau and B. Derrida, Current fluctuations in nonequilibrium diffusive systems: An additivity princi- ple, Phys. Rev. Lett. 92, 180601 (2004)

2004
[54]

Arita, P

C. Arita, P. L. Krapivsky, and K. Mallick, Generalized exclusion processes: Transport coefficients, Phys. Rev. E 90, 052108 (2014)

2014
[55]

Please see the appendices for detail,
[56]

Alexander and P

S. Alexander and P. Pincus, Diffusion of labeled particles on one-dimensional chains, Phys. Rev. B18, 2011 (1978)

1978
[57]

P. L. Krapivsky, K. Mallick, and T. Sadhu, Tagged parti- cle in single-file diffusion, J. Stat. Phys. 160, 885 (2015)

2015
[58]

P. L. Krapivsky, K. Mallick, and T. Sadhu, Large devia- tions in single-file diffusion, Phys. Rev. Lett. 113, 078101 (2014)

2014
[59]

Derrida and A

B. Derrida and A. Gerschenfeld, Current fluctuations in one dimensional diffusive systems with a step initial den- sity profile, J. Stat. Phys. 137, 978 (2009)

2009
[60]

Sadhu and B

T. Sadhu and B. Derrida, Correlations of the density and of the current in non-equilibrium diffusive systems, J. Stat. Mech. 2016, 113202 (2016)

2016
[61]

P. L. Krapivsky and B. Meerson, Fluctuations of current in nonstationary diffusive lattice gases, Phys. Rev. E 86, 031106 (2012)

2012
[62]

A. Das, A. Kundu, and P. Pradhan, Einstein relation and hydrodynamics of nonequilibrium mass transport pro- cesses, Phys. Rev. E 95, 062128 (2017)

2017
[63]

Chakraborty, S

T. Chakraborty, S. Chakraborti, A. Das, and P. Pradhan, Hydrodynamics, superfluidity, and giant number fluctu- ations in a model of self-propelled particles, Phys. Rev. E 101, 052611 (2020)

2020
[64]

Agranov, S

T. Agranov, S. Ro, Y. Kafri, and V. Lecomte, Macro- scopic fluctuation theory and current fluctuations in ac- tive lattice gases, SciPost Phys. 14, 045 (2023)

2023
[65]

S. Jose, R. Dandekar, and K. Ramola, Current fluctua- tions in an interacting active lattice gas, J. Stat. Mech. 2023, 083208 (2023)

2023
[66]

Chakraborty and P

T. Chakraborty and P. Pradhan, Time-dependent prop- erties of run-and-tumble particles. ii. current fluctua- tions, Phys. Rev. E 109, 044135 (2024)

2024
[67]

Godr` eche, Dynamics of condensation in zero-range processes, J

C. Godr` eche, Dynamics of condensation in zero-range processes, J. Phys. A: Math. Gen. 36, 6313 (2003)

2003
[68]

Rizkallah, A

P. Rizkallah, A. Grabsch, P. Illien, and O. B´ enichou, Du- ality relations in single-file diffusion, J. Stat. Mech.2023, 013202 (2023). 13 VI. APPENDIX A. Derivation of hydrodynamics of 2D SSEP in EM Our two-dimensional model with hard core occupancy can be described with occ...

2023

Pith tools

Reviewed August 6, 2026 · model on record in the stance chip above.

[1] [1]

Ben-Jacob, I

E. Ben-Jacob, I. Cohen, and D. L. Gutnick, Coopera- tive organization of bacterial colonies: From genotype to morphotype, Annu. Rev. Microbiol. 52, 779 (1998)

work page 1998

[2] [2]

Taktikos, Y

J. Taktikos, Y. T. Lin, H. Stark, N. Biais, and V. Zabur- daev, Pili-induced clustering of n. gonorrhoeae bacteria, PLoS ONE 10, e0137661 (2015)

work page 2015

[3] [3]

Mueller-Klieser, Multicellular spheroids: A review on cellular aggregates in cancer research, J

W. Mueller-Klieser, Multicellular spheroids: A review on cellular aggregates in cancer research, J. Cancer. Res. Clin. Oncol. 113, 101 (1987)

work page 1987

[4] [4]

q − p ∂ ∂ρ ( ρ ∞X m=0 mv(m)P (m | ρ) )# , (30) χ(ρ) = 1 4

= 0 .85 and s2(l0 = 8) = 0 .9. We further check (not shown) that v(m) does not depend on density ρ and cell- cell interaction rate p. Throughout this work we consider T0 = 10 MCT, Np = 20 and l0 = 2, therefore, we set s1 = 0.4 and s2 = 0.8 for the rest of this paper. The expression for jump weight v(m) given by Eq. (6) is the key starting point allowing t...

work page

[5] [5]

Friedl and K

P. Friedl and K. Wolf, Tumour-cell invasion and migra- tion: diversity and escape mechanisms, Nat. Rev. Cancer 3, 362 (2003)

work page 2003

[6] [6]

T. M. Koch, S. Munster, N. Bonakdar, J. P. Butler, and B. Fabry, 3d traction forces in cancer cell invasion, PLoS ONE 7, e33476 (2012)

work page 2012

[7] [7]

Mark et al., Collective forces of tumor spheroids in three-dimensional biopolymer networks, eLife 9, e51912 (2020)

C. Mark et al., Collective forces of tumor spheroids in three-dimensional biopolymer networks, eLife 9, e51912 (2020)

work page 2020

[8] [8]

Karow et al., Direct pericyte-to-neuron reprogram- ming via unfolding of a neural stem cell-like program, Nat

M. Karow et al., Direct pericyte-to-neuron reprogram- ming via unfolding of a neural stem cell-like program, Nat. Neurosci. 21, 932 (2018)

work page 2018