Pith. sign in

REVIEW 3 major objections 4 minor 70 references

A nonlocal-to-local approach to aggregation-diffusion equations

T0 review · 3 major / 4 minor · reviewed 2026-08-15 · deepseek-v4-flash

Pith's one-line read The paper claims that short-range nonlocal adhesion can be replaced by a four-parameter thin-film-like system that still produces the experimentally observed cell-sorting configurations.

desk verdict A clear, honest synthesis of the authors' earlier local adhesion model, but the 'preserves the same phenomenology' claim is asserted rather than tested. read the letter →

arxiv 2505.08443 v1 pith:RJJYK633 submitted 2025-05-13 q-bio.CB math.AP

classification q-bio.CBmath.AP MSC 35Q9292C1535K55
keywords aggregation-diffusioncell-celladhesionthin-filmequationCahn-Hilliarddifferentialhypothesiscellsortingnonlocal-to-locallimitgradientflow
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

The paper claims that the nonlocal equations commonly used for cell-cell adhesion can, in the short-range limit, be replaced by a system of fourth-order thin-film-like equations with only four parameters. These parameters carry biological meaning: they act as relative surface tensions and adhesion strengths between two cell populations. The authors show numerically, in one and two dimensions, that this local system produces the four configurations predicted by the differential adhesion hypothesis: sorting, partial engulfment, engulfment, and mixing. If the claim holds, the simpler local model gives a tractable route to calibrating adhesion models against experiments and to analytical steady-state predictions that nonlocal models rarely offer.

What carries the argument

The load-bearing object is the reduced local system (9), a system of two fourth-order degenerate parabolic equations of thin-film type. It is obtained by writing each nonlocal convolution as a Taylor expansion in the sensing radius and keeping only the zeroth- and second-moment terms; symmetry kills first-order terms and the error is $O(a^4 M_4)$. The machinery does three jobs: it supplies four interpretable parameters (relative self-adhesion $\kappa$, cross-adhesion $\alpha$, and the corresponding population-pressure terms $\mu$ and $\omega$), it preserves the gradient-flow structure through the free energy (11), and it permits explicit stationary bump profiles that match numerics.

What would settle it

Run the nonlocal system (2) and the local system (9) with identical compactly supported kernels mapped through the moment formulas, across the parameter plane of Figure 4; any region where the stationary patterns disagree, such as sorting instead of engulfment near the $\alpha=\sqrt{\kappa}$ boundary, or metastable states that persist only in one model, would falsify the claim that the local limit preserves the phenomenology.

Watch

Extended reading notes

Core claim

On its own terms, the paper establishes that the two-species nonlocal aggregation-diffusion system (2), with convolution kernels $W_{ij}$, formally converges, after Taylor expansion and truncation at second order, to the local system (9): a set of thin-film/Cahn-Hilliard type equations with parameters $\kappa$, $\alpha$, $\mu$, $\omega$ that are ratios of the first and second moments of the interaction potentials plus the repulsion strength. The central finding is that this local limit preserves the phenomenology of differential adhesion: numerical simulations identify parameter regimes matching the four experimentally observed cell-sorting configurations, with the cross-interaction parameters $\alpha$ and $\omega$ controlling the transition from sorting to mixing. The same gradient-flow energy structure of the nonlocal model is retained, giving an energy that decreases along solutions and a basis for existence results reported by the authors in companion work.

Load-bearing premise

The load-bearing premise is that truncating the Taylor expansion of the nonlocal convolution after the second moment, valid for small sensing radius and nicely decaying kernels, does not change the biologically relevant dynamics, so that the local model's patterns are the nonlocal model's patterns.

Editorial extensions

If this is right

  • A nonlocal adhesion model can be simulated and calibrated as a local PDE, avoiding costly convolution evaluations while keeping the qualitative patterns.
  • The four parameters map directly to measurable or inferable quantities, so inference techniques can target parameters rather than whole interaction kernels.
  • The local model admits analytical stationary solutions for individual aggregates, something typically unavailable for nonlocal models with general potentials.
  • Linear stability of the homogeneous state is governed by the sign of $\mu_2$ ($M_0-\epsilon$), matching the nonlocal model's condition for existence of stationary states.
  • Both one- and two-dimensional simulations show the same four experimentally observed configurations, so the simplification is not a one-dimensional artifact.

Reading between the lines

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

  • If the truncation is faithful, the same nonlocal-to-local reduction could be applied to other attraction-repulsion systems, such as crowd dynamics, chemotaxis, or social aggregation, wherever interaction kernels are effectively short-range.
  • A testable extension would be to map experimental adhesion measurements into the four parameters and predict tissue-sorting outcomes quantitatively, rather than only qualitatively.
  • The reported metastable coarsening in the weak cross-adhesion regime suggests the local model could serve as a computationally cheap surrogate for studying long-transient dynamics in nonlocal models, a connection the paper only partially explores.
  • The gradient-flow energy (11) makes this system a convenient testbed for conjectures about Cahn-Hilliard systems with cross-interaction terms, such as wetting-layer behaviour or spinodal decomposition analogies.
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.

Referee Report

3 major / 4 minor

Summary. The paper derives a two-species local aggregation-diffusion system, Eqs. (9), by Taylor-expanding the nonlocal convolution terms in Eqs. (2) to second order, assuming short-range symmetric attractive kernels. It then studies the one-species linear stability and energy dissipation, and uses numerical simulations to argue that the model reproduces the four Steinberg configurations: sorting, partial engulfment, engulfment, and mixing. The paper also surveys recent existence results for the local model and lists open problems.

Significance. If the nonlocal-to-local reduction is faithful at finite sensing radius, the model is a genuinely useful simplification: it has a gradient-flow structure with explicit energy, analytically accessible steady states, four interpretable parameters, cheaper numerical implementation, and a more direct path to calibration than the original nonlocal model. The paper has clear strengths: the formal Taylor derivation and nondimensionalization in Sections 2.1 and 3.1 are internally consistent; the dispersion relation sigma(k) = rho0 |k|^2 (mu2 - |k|^2) recovers the expected Cahn-Hilliard instability condition; the energy-dissipation calculations for F[rho] and F2[rho,eta] are correct under the stated boundary conditions; and the numerical illustrations in Figures 4, 5, and 7, together with the VisualPDE implementations, make the phenomenology accessible. The main unresolved issue is that fidelity to the nonlocal model, which is advertised in the title and abstract, is asserted rather than quantitatively tested.

major comments (3)
  1. [Section 2.1, Eq. (4)] The abstract's claim that the local model 'preserves the same phenomenology' as the nonlocal model is not established. The Taylor expansion truncates W*f at the second moment and drops an O(a^4 M_4) remainder, and the paper explicitly states that its goal is not to compare Eq. (4) with Eq. (1). The cited rigorous results [22,34] prove convergence in suitable limits as a tends to zero, but they do not quantify the finite-a error in metastable dynamics, stationary-pattern selection, or phase boundaries; the paper itself notes in Section 4 that [18] reports significant local/nonlocal differences for non-compact kernels. Since the central claim of the paper rests on this nonlocal-to-local fidelity, the authors should either provide a quantitative truncation estimate in terms of a and the moments M_4, or include a direct numerical comparison of Eqs. (9) with Eqs. (2) on the same parameter sweeps and initial data over the simulated time scales.
  2. [Sections 3.3-3.4, Figures 4, 5, 7] The claim that the model identifies 'parameter regimes' for sorting, partial engulfment, engulfment, and mixing is supported only by one hand-selected (alpha, omega) point per pattern, with kappa = 2 and mu = 4 fixed. The shaded regions in Figure 4 are not specified by precise inequalities, and no sensitivity analysis or phase-diagram computation is presented. Because the system is a degenerate fourth-order gradient flow and may have multiple metastable stationary states, the manuscript should demonstrate that each pattern occurs on an open parameter set and is robust to initial conditions and numerical resolution before claiming regime identification.
  3. [Section 3.3, Eq. (13) and surrounding text] The parameter relations kappa = K11/K22, alpha = K12/K22, mu = ((M0 - eps/K11)/(M0 - eps/K22)) kappa, and omega = ((M0 - eps/K12)/(M0 - eps/K22)) alpha are derived under the assumption phi_ij = K_ij phi, but the numerical simulations do not state which values of M0, epsilon, K11, K22, and K12 produce the chosen (kappa, mu, alpha, omega). Without this mapping, the connection between the 'clear physical interpretation' of the parameters and the actual simulations remains qualitative; the paper should either make this correspondence explicit or clearly separate the mathematical parameter space from the biological calibration.
minor comments (4)
  1. [Section 2.1] The sentence 'the energy is non-decreasing in time' should read 'non-increasing', since the calculation gives dF/dt <= 0.
  2. [Section 2.3] The finite-volume scheme is only referenced; the paper should state the specific discretization or point to equation numbers in [39], and add at least one grid-refinement check for a two-species steady state to support the numerical claims.
  3. [Section 3.3, Figure 4] The 'weak' and 'strong' cross-adhesion regions should be defined by explicit inequalities rather than only by a sketch, so that the claimed regimes are reproducible by readers.
  4. [Figure 6 caption] The energy label in the caption appears garbled in the rendered manuscript; please verify the mathematical notation in the final version.

Circularity Check

0 steps flagged · score 0.0 of 10

No significant circularity: Eqs. (9) are obtained by a direct Taylor expansion of the nonlocal convolution, and the numerical reproduction of Steinberg patterns is a model-consistency test rather than a fitted prediction.

full rationale

The formal derivation of Eq. (9) from Eqs. (2) is a direct Taylor expansion of the convolution, keeping the zeroth and second moments and dropping the O(a^4) remainder; radial symmetry makes the odd moments vanish, and no target quantity is inserted into the expansion to force the result. The later numerical studies do not fit parameters to data and then relabel the fit as a prediction: they choose representative parameter values in the weak and strong cross-adhesion regimes and demonstrate that the local model can reproduce the four Steinberg configurations. This is a capability or consistency demonstration, not a statistically forced prediction. The abstract's claim that the local model 'preserves the same phenomenology' is asserted rather than quantitatively verified, and the paper explicitly disclaims a local-versus-nonlocal comparison ('the goal of our paper is not to compare (4) with the nonlocal model (1)'), but an unverified assertion is a validation gap, not circularity. The cited works [22], [23], and [39] include self-citations, but they are used for rigorous convergence results, existence theory, and companion technical details; they are not invoked as the source of the central derivation, and the rigorous limits are stated with assumptions that are independent of the target pattern-reproduction claim. No equation was found that reduces to its own input by construction, and no fitted parameter was renamed as a prediction.

Assumptions & free parameters 5 free parameters · 5 assumptions · 0 invented entities

The model's parameters are not inferred from data; alpha and omega are hand-picked to reach each of the four target patterns, so the numerical 'confirmation' is a demonstration of regimes rather than a predictive test. The formal derivation relies on moment assumptions and on cited rigorous limits [22,34] and existence results [23], some authored by the present group. No invented entities are posited.

free parameters (5)
  • mu2 (one-species attraction parameter) = mu2=1 in one-species simulations
    Controls linear instability of the homogeneous state in Eq. (5); fixed in simulations and not inferred from data.
  • kappa (relative self-adhesion of rho to eta) = kappa=2 in all two-species simulations
    Chosen so that eta is the less cohesive population; not measured from experiments.
  • mu (relative population pressure of rho to eta) = mu=4 in all two-species simulations
    Chosen to satisfy mu>kappa>1, matching the assumed adhesion ordering; not inferred from data.
  • alpha (cross-diffusion / surface-tension coupling) = sorting alpha=0; partial engulfment alpha=0.8 (1D) or 0.5 (2D); engulfment alpha=1.3; mixing alpha=1.4
    Hand-picked to produce each of the four Steinberg patterns in Figures 4-7; the paper does not calibrate alpha to independent data.
  • omega (cross-attraction / repulsion strength) = sorting omega=-1; partial engulfment omega=0.2 (1D) or -0.02 (2D); engulfment omega=2; mixing omega=6 (1D) or 8 (2D)
    Chosen by hand to move the system from sorting to mixing; this is a parameter-regime demonstration, not a calibrated prediction.
assumptions (5)
  • domain assumption The mean-field limit from the particle system to the nonlocal aggregation-diffusion equation (1) holds.
    Invoked in Section 1 as the starting point of the model; cites Oelschlager [59] for the rigorous scaling limit.
  • domain assumption The interaction potential W is radially symmetric, attractive, smooth, with W=-a^{-d} phi(x/a) for small a, and phi is non-increasing with suitably decaying moments.
    Required for the Taylor expansion and for vanishing odd-order terms in Section 2.1.
  • ad hoc to paper The remainder O(a^4 M4) is negligible, so truncating after the second moment preserves the dynamics.
    The paper notes this is only valid when higher moments are small and a is small, but it does not quantify the truncation error for the parameter regimes used in the simulations.
  • domain assumption Cross-interaction potentials are symmetric (W12=W21) and all potentials share a common shape phi_ij=K_ij phi.
    Used in Sections 3.1 and 3.3 to reduce the model to four parameters; follows common practice in [5,26] but restricts generality.
  • domain assumption Rigorous nonlocal-to-local convergence theorems [22,34] and existence results [23] are accepted as background.
    The paper relies on these cited results without proof; two of them are by the present authors.

how reviews work

0 comments
Cite this review

Pith. "Pith review of A nonlocal-to-local approach to aggregation-diffusion equations." pith.science (2026). https://pith.science/paper/RJJYK633

@misc{pith2026250508443,
  author       = {Pith},
  title        = {Pith review of: A nonlocal-to-local approach to aggregation-diffusion equations},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/RJJYK633}},
  note         = {Machine review of arXiv:2505.08443}
}
read the original abstract

Over the past decades, nonlocal models have been widely used to describe aggregation phenomena in biology, physics, engineering, and the social sciences. These are often derived as mean-field limits of attraction-repulsion agent-based models, and consist of systems of nonlocal partial differential equations. Using differential adhesion between cells as a biological case study, we introduce a novel local model of aggregation-diffusion phenomena. This system of local aggregation-diffusion equations is fourth-order, resembling thin-film or Cahn-Hilliard type equations. In this framework, cell sorting phenomena are explained through relative surface tensions between distinct cell types. The local model emerges as a limiting case of short-range interactions, providing a significant simplification of earlier nonlocal models, while preserving the same phenomenology. This simplification makes the model easier to implement numerically and more amenable to calibration to quantitative data. Additionally, we discuss recent analytical results based on the gradient-flow structure of the model, along with open problems and future research directions.

Figures

Figures reproduced from arXiv: 2505.08443 by the authors.

Figure 1
Figure 1. Aggregation is possible in the local model as long as [PITH_FULL_IMAGE:figures/full_fig_p008_1.png] view at source ↗
Figure 2
Figure 2. (left) Aggregation in the two-dimensional local model. Initial data is [PITH_FULL_IMAGE:figures/full_fig_p010_2.png] view at source ↗
Figure 3
Figure 3. Possible configurations for Steinberg experiments in terms of the cross-adhesion and the self [PITH_FULL_IMAGE:figures/full_fig_p011_3.png] view at source ↗
Figures from the paper (4 more)
Figure 4
Figure 4. Figure 4: Understanding the impact of changing model parameters. Imposing that [PITH_FULL_IMAGE:figures/full_fig_p015_4.png]
Figure 5
Figure 5. Figure 5: Solutions of the local model using model parameters related to the Steinberg experiments. Each [PITH_FULL_IMAGE:figures/full_fig_p016_5.png]
Figure 6
Figure 6. Figure 6: (left) Energy decay given by Eq. (11) for the numerical solutions in [PITH_FULL_IMAGE:figures/full_fig_p017_6.png]
Figure 7
Figure 7. Figure 7: Numerical solutions of the local model using model parameters related to Steinberg experiments. [PITH_FULL_IMAGE:figures/full_fig_p018_7.png]

Discussion (0). Continue with ORCID to comment.

Reference graph

Works this paper leans on

70 extracted references · 67 canonical work pages

  1. [39]

    Falc´ o, R

    C. Falc´ o, R. E. Baker, and J. A. Carrillo. A local continuum model of cell-cell adhesion SIAM Journal on Applied Mathematics , 84(3):S17–S42, 2024

  2. [23]

    J. A. Carrillo, A. Esposito, C. Falc´ o and A. Fern´ andez-Jim´ enez. Competing effects in fourth-order aggregation–diffusion equations. Proceedings of the London Mathematical Society, 129:e12623, 2024

  3. [18]

    Buttensch¨ on, S

    A. Buttensch¨ on, S. Sinclair, and L. Edelstein-Keshet. How cells stay together: a mechanism for mainte- nance of a robust cluster explored by local and non-local continuum models. Bulletin of Mathematical Biology, 86(11):129, 2024

  4. [1]

    Alert and X

    R. Alert and X. Trepat. Physical models of collective cell migration. Annual Review of Condensed Matter Physics, 11(1):77–101, 2020

  5. [2]

    S. Alt, P. Ganguly, and G. Salbreux. Vertex models: from cell mechanics to tissue morphogenesis. Philosophical Transactions of the Royal Society B: Biological Sciences , 372:20150520, 05 2017

  6. [3]

    J. D. Amack and M. L. Manning. Knowing the boundaries: extending the differential adhesion hypothesis in embryonic cell sorting. Science, 338(6104):212–215, 2012

  7. [4]

    Anguige and C

    K. Anguige and C. Schmeiser. A one-dimensional model of cell diffusion and aggregation, incorporating volume filling and cell-to-cell adhesion. Journal of Mathematical Biology , 58:395–427, 07 2008

  8. [5]

    N. J. Armstrong, K. J. Painter, and J. A. Sherratt. A continuum approach to modelling cell–cell adhesion. Journal of Theoretical Biology, 243(1):98–113, 2006

Show all 70 references
  1. [6]

    Bailo, J

    R. Bailo, J. A. Carrillo, and P. Degond. Pedestrian models based on rational behaviour. Crowd Dynamics, Volume 1: Theory, Models, and Safety Problems, 259–292

  2. [7]

    Bailo, J

    R. Bailo, J. A. Carrillo, D. G´ omez-Castro. Aggregation-diffusion equations for collective behaviour in the sciences. arxiv preprint arXiv:2405.16679, 2024

  3. [8]

    Bailo, J

    R. Bailo, J. A. Carrillo, H. Murakawa, and M. Schmidtchen. Convergence of a fully discrete and energy- dissipating finite-volume scheme for aggregation-diffusion equations. Mathematical Models and Methods in Applied Sciences, 30(13):2487–2522, 2020

  4. [9]

    Bailo, J

    R. Bailo, J. A. Carrillo, S. Kalliadasis, and S. P. Perez. Unconditional bound-preserving and energy- dissipating finite-volume schemes for the Cahn-Hilliard equation. Communications in Computational Physics, 34(3):713–748, 2023

  5. [10]

    R. E. Baker and M. J. Simpson. Models of collective cell motion for cell populations with different aspect ratio: diffusion, proliferation and travelling waves. Physica A: Statistical Mechanics and its Applications, 391(14):3729–3750, 2012

  6. [11]

    J. W. Barrett, J. F. Blowey, and H. Garcke. On fully practical finite element approximations of degenerate Cahn-Hilliard systems. ESAIM: Mathematical Modelling and Numerical Analysis , 35(4):713–748, 2001

  7. [12]

    Bernoff and C

    A. Bernoff and C. Topaz. Biological aggregation driven by social and environmental factors: a nonlocal nodel and its degenerate Cahn-Hilliard approximation. SIAM Journal on Applied Dynamical Systems , 15(3):1528–1562, 2016

  8. [13]

    A. L. Bertozzi and M. C. Pugh. The lubrication approximation for thin viscous films: regularity and long-time behavior of weak solutions. Communications on Pure and Applied Mathematics, 49(2):85–123, 1996. 21

  9. [14]

    A. L. Bertozzi and M. C. Pugh. Long-wave instabilities and saturation in thin film equations. Commu- nications on Pure and Applied Mathematics , 51(6):625–661, 1998

  10. [15]

    Brodland

    G. Brodland. The Differential Interfacial Tension Hypothesis (DITH): A Comprehensive theory for the self-rearrangement of embryonic cells and tissues. Journal of Biomechanical Engineering , 124:188–97, 2002

  11. [16]

    Burger, M

    M. Burger, M. Di Francesco, and M. Franek. Stationary states of quadratic diffusion equations with long-range attraction. Communications in Mathematical Sciences , 11:709–738, 2013

  12. [17]

    Buttensch¨ on, and L

    A. Buttensch¨ on, and L. Edelstein-Keshet. Bridging from single to collective cell migration: A review of models and links to experiments. PLOS Computational Biology , 16(12):e1008411, 2020

  13. [19]

    Calvez and J

    V. Calvez and J. A. Carrillo. Volume effects in the Keller–Segel model: energy estimates preventing blow-up. Journal de Math´ ematiques Pures et Appliqu´ ees, 86(2):155–175, 2006

  14. [20]

    J. A. Carrillo, A. Chertock, and Y. Huang. A finite-volume method for nonlinear nonlocal equations with a gradient flow structure. Communications in Computational Physics , 7(1):233–258, 2015

  15. [21]

    J. A. Carrillo, A. Colombi, and M. Scianna. Adhesion and volume constraints via nonlocal interactions determine cell organisation and migration profiles. Journal of Theoretical Biology, 445:75–91, 2018

  16. [22]

    J. A. Carrillo, C. Elbar, and J. Skrzeczkowski. Degenerate Cahn-Hilliard systems: from nonlocal to local. arxiv preprint arXiv:2303.11929, to appear in Communications in Contemporary Mathematics

  17. [24]

    J. A. Carrillo, Y. Huang, and M. Schmidtchen. Zoology of a nonlocal cross-diffusion model for two species. SIAM Journal on Applied Mathematics , 78(2):1078–1104, 2018

  18. [25]

    J. A. Carrillo, R. J. McCann, and C. Villani. Kinetic equilibration rates for granular media and related equations: entropy dissipation and mass transportation estimates. Revista Matematica Iberoamericana, 19(3):971–1018, 2003

  19. [26]

    J. A. Carrillo, H. Murakawa, M. Sato, H. Togashi, and O. Trush. A population dynamics model of cell- cell adhesion incorporating population pressure and density saturation. Journal of Theoretical Biology, 474:14–24, 2019 22

  20. [27]

    G. L. Celora, M. G. Hennessy, A. M¨ unch, B. Wagner, and S. L. Waters. The dynamics of a collapsing polyelectrolyte gel. SIAM Journal on Applied Mathematics , 83(3):1146–1171, 2023

  21. [28]

    L. Chen, K. Painter, C. Surulescu, and A. Zhigun. Mathematical models for cell migration: a non-local perspective. Philosophical Transactions of the Royal Society B , 375(1807):20190379, 2020

  22. [29]

    M. G. Delgadino. Convergence of a One-Dimensional Cahn–Hilliard Equation with Degenerate Mobility. SIAM Journal on Mathematical Analysis , 50(4):4457–4482, 2018

  23. [30]

    Di Francesco and S

    M. Di Francesco and S. Fagioli. Measure solutions for non-local interaction PDEs with two species. Nonlinearity, 26(10):2777, 2013

  24. [31]

    Duguay, R

    D. Duguay, R. A. Foty, and M. S. Steinberg. Cadherin-mediated cell adhesion and tissue segregation: qualitative and quantitative determinants. Developmental Biology, 253(2):309–323, 2003

  25. [32]

    Dyson and R

    L. Dyson and R. E. Baker. The importance of volume exclusion in modelling cellular migration. Journal of Mathematical Biology, 71:691–711, 2015

  26. [33]

    Dyson, P

    L. Dyson, P. K. Maini, and R. E. Baker. Macroscopic limits of individual-based models for motile cell populations with volume exclusion. Physical Review E, 86:031903, 2012

  27. [34]

    Elbar and J

    C. Elbar and J. Skrzeczkowski. Degenerate Cahn-Hilliard equation: From nonlocal to local. Journal of Differential Equations, 364:576–611, 2023

  28. [35]

    Ellefsen and N

    E. Ellefsen and N. Rodriguez. On equilibrium solutions to nonlocal mechanistic models in ecology. Journal of Applied Analysis & Computation , 11(6):2664–2686, 2021

  29. [36]

    C. M. Elliott. The Cahn–Hilliard model for the kinetics of phase separation. In Mathematical models for phase change problems , pages 35–73. Springer, 1989

  30. [37]

    C. M. Elliott and H. Garcke. On the Cahn–Hilliard equation with degenerate mobility. SIAM Journal on Mathematical Analysis, 27(2):404–423, 1996

  31. [38]

    C. Falc´ o. From random walks on networks to nonlinear diffusion. Physical Review E, 106:054103, 2022

  32. [40]

    Falc´ o, R

    C. Falc´ o, R. E. Baker, and J. A. Carrillo. A local continuum model of cell-cell adhesion: figshare me- dia, 2022. https://figshare.com/projects/A_local_continuum_model_of_cell-cell_adhesion/ 142427

  33. [41]

    Falc´ o, D

    C. Falc´ o, D. J. Cohen, J. A. Carrillo, and R. E. Baker. Quantifying tissue growth, shape and collision via continuum models and Bayesian inference. Journal of the Royal Society Interface , 20:20230184, 2023 23

  34. [42]

    R. A. Foty and M. S. Steinberg. Cadherin-mediated cell-cell adhesion and tissue segregation in relation to malignancy. International Journal of Developmental Biology , 48(5-6):397–409, 2004

  35. [43]

    R. A. Foty and M. S. Steinberg. The differential adhesion hypothesis: a direct evaluation. Developmental Biology, 278(1):255–263, 2005

  36. [44]

    Di Francesco, A

    M. Di Francesco, A. Esposito, and S. Fagioli. Nonlinear degenerate cross-diffusion systems with nonlocal interaction. Nonlinear Analysis, Theory, Methods and Applications , 169:94–117, 2018

  37. [45]

    R. E. Goldstein, A. I. Pesci, and M. J. Shelley. Topology transitions and singularities in viscous flows. Physical Review Letters, 70:3043–3046, 1993

  38. [46]

    G´ omez-Castro, Beginner’s guide to aggregation-diffusion equations.SeMA Journal, 81:531–587, 2024

    D. G´ omez-Castro, Beginner’s guide to aggregation-diffusion equations.SeMA Journal, 81:531–587, 2024

  39. [47]

    R. E. Goldstein, A. I. Pesci, and M. J. Shelley. Instabilities and singularities in Hele-Shaw flow. Physics of Fluids, 10(11):2701–2723, 1998

  40. [48]

    M. E. Gurtin and R. C. MacCamy. On the diffusion of biological populations. Mathematical Biosciences, 33(1-2):35–49, 1977

  41. [49]

    Hashimoto, A

    A. Hashimoto, A. Nagao, and S. Okuda. Topological graph description of multicellular dynamics based on vertex model. Journal of Theoretical Biology, 437:187–201, 2018

  42. [50]

    Hirashima, E

    T. Hirashima, E. G. Rens, and R. M. Merks. Cellular Potts modeling of complex multicellular behaviors in tissue morphogenesis. Development, Growth & Differentiation , 59(5):329–339, 2017

  43. [51]

    Krieg, Y

    M. Krieg, Y. Arboleda-Estudillo, P.-H. Puech, J. K¨ afer, F. Graner, D. M¨ uller, and C.-P. Heisenberg. Tensile forces govern germ-layer organization in zebrafish. Nature Cell Biology , 10(4):429–436, 2008

  44. [52]

    Laugesen and M

    R. Laugesen and M. Pugh. Linear stability of steady states for thin film and Cahn-Hilliard type equations. Archive for Rational Mechanics and Analysis , 154:3–51, 2000

  45. [53]

    Laugesen and M

    R. Laugesen and M. Pugh. Energy levels of steady states for Thin-Film-type equations. Journal of Differential Equations, 182(2):377–415, 2002

  46. [54]

    Laugesen and M

    R. Laugesen and M. Pugh. Properties of steady states for thin film equations. European Journal of Applied Mathematics, 11:293–351, 2000

  47. [55]

    M. L. Manning, R. A. Foty, M. S. Steinberg, and E.-M. Schoetz. Coaction of intercellular adhesion and cortical tension specifies tissue surface tension. Proceedings of the National Academy of Sciences , 107(28):12517–12522, 2010

  48. [56]

    Matthes, R

    D. Matthes, R. J. McCann, and G. Savar´ e. A family of nonlinear fourth order equations of gradient flow type. Communications in Partial Differential Equations , 34(11):1352–1397, 2009. 24

  49. [57]

    Murakawa and H

    H. Murakawa and H. Togashi. Continuous models for cell–cell adhesion. Journal of Theoretical Biology, 374:1–12, 2015

  50. [58]

    T. G. Myers. Thin films with high surface tension. SIAM Review, 40(3):441–462, 1998

  51. [59]

    Oelschl¨ ager

    K. Oelschl¨ ager. Large systems of interacting particles and the porous medium equation. Journal of Differential Equations, 88:294–346, 1990

  52. [60]

    J. R. Potts and K. J. Painter. Distinguishing between long-transient and asymptotic states in a biological aggregation model. Bulletin of Mathematical Biology , 86(3):28, 2024

  53. [61]

    Santambrogio

    F. Santambrogio. Optimal Transport for Applied Mathematicians. Birk¨ auser, NY, 55(58-63):94, 2015

  54. [62]

    Slepˇ cev

    D. Slepˇ cev. Linear stability of selfsimilar solutions of unstable thin-film equations. Interfaces and Free Boundaries, 11(3):375–398, 2009

  55. [63]

    M. S. Steinberg. Mechanism of tissue reconstruction by dissociated cells, II: Time-course of events. Science, 137(3532):762–763, 1962

  56. [64]

    M. S. Steinberg. On the mechanism of tissue reconstruction by dissociated cells, I. Population kinetics, differential adhesiveness, and the absence of directed migration. Proceedings of the National Academy of Sciences, 48(9):1577, 1962

  57. [65]

    M. S. Steinberg. On the mechanism of tissue reconstruction by dissociated cells, III. Free energy relations and the reorganization of fused, heteronomic tissue fragments. Proceedings of the National Academy of Sciences, 48(10):1769, 1962

  58. [66]

    M. S. Steinberg. Reconstruction of tissues by dissociated cells: some morphogenetic tissue movements and the sorting out of embryonic cells may have a common explanation. Science, 141(3579):401–408, 1963

  59. [67]

    T. Y.-C. Tsai, R. M. Garner, and S. G. Megason. Adhesion-based self-organization in tissue patterning. Annual Review of Cell and Developmental Biology , 38(1):349–374, 2022

  60. [68]

    Volkening, M

    A. Volkening, M. R. Abbott, N. Chandra, B. Dubois, F. Lim, D. Sexton, and B. Sandstede. Modeling stripe formation on growing zebrafish tailfins. Bulletin of Mathematical Biology , 82:1–33, 2020

  61. [69]

    B. J. Walker, A. K. Townsend, A. K. Chudasama, and A. L. Krause. VisualPDE: Rapid Interactive Simulations of Partial Differential Equations. Bulletin of Mathematical Biology , 85(113), 2023

  62. [70]

    Youssef, A

    J. Youssef, A. K. Nurse, L. Freund, and J. R. Morgan. Quantification of the forces driving self-assembly of three-dimensional microtissues. Proceedings of the National Academy of Sciences, 108(17):6993–6998, 2011. 25

Pith tools

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