Pith. sign in

REVIEW 3 major objections 4 minor 1 cited by

Existence of traveling wave for a coupled incompressible Darcy's free boundary model with undercooling effect and surface tension

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

Pith's one-line read This paper proves that a Darcy free-boundary cell model with membrane undercooling has fixed-area traveling-wave solutions once marker activity exceeds a threshold.

desk verdict A plausible and useful extension with a clean idea, but the bifurcation proof has a false base point and both main theorems outrun the proofs. read the letter →

arxiv 2501.04576 v1 pith:R43SUBBT submitted 2025-01-08 math.AP

classification math.AP MSC 35R3535B3235B3535C0792C17
keywords travelingwavesfreeboundaryproblemDarcyflowcellmotilitybifurcationlinearstabilityundercoolingsurfacetension
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

This paper aims to show that a Darcy free-boundary model of confined cell motility, supplemented by a membrane undercooling force and a marker-driven active force, reproduces spontaneous persistent motion. Below a critical marker strength $\chi_c^*$ the unique resting disk is claimed to be linearly stable, and above it linearly unstable. The main theorem claims that for every $a\in(0,1]$, surface tension, radius, and positive membrane strength, a one-parameter family of fixed-area traveling-wave solutions bifurcates from the disk for marker strengths above the threshold. These traveling waves are the mathematical signature of a cell that polarizes and moves without an external cue, so proving their existence is what makes the model biologically relevant.

What carries the argument

The load-bearing object is the functional $F:\mathbb{R}\times X\times\mathbb{R}\times\mathbb{R}\to Y\times\mathbb{R}\times\mathbb{R}$ defined in (4.3), whose zero set encodes the boundary curvature equation for $x$-symmetric, $2\pi$-periodic perturbations $\rho$ of the disk together with the fixed-area and centering constraints. Its derivative at $(\chi_c^*,0,0,0)$ has kernel $\mathrm{span}\{(0,1,0)\}$, and the Crandall–Rabinowitz transversality condition is verified, so the implicit branch of solutions is produced. The threshold $\chi_c^*=(R_0+\chi_u f_{\mathrm{und}}'(0))/(R_0 a c_0 f_{\mathrm{act}}'(c_0))$ comes from the linearized eigenvalue calculation, and the undercooling term enters both the threshold and the branch-direction formula in Lemma 4.3.

What would settle it

Substitute $\rho=0$, $V=0$, $p_1=0$ into the first component of $F$ in (4.3): it equals $\chi_c f_{\mathrm{act}}(c_0)$, not zero, so the claimed trivial branch does not hold as written. Rerunning the kernel and transversality computations with the resting marker force absorbed into $p_1$ will show whether the one-dimensional kernel and the Crandall–Rabinowitz conditions survive, and a forward simulation of (1.1) just above $\chi_c^*$ would check whether the predicted translating shape actually appears.

Watch

Extended reading notes

Core claim

On the paper's own terms, the central discovery is a bifurcation result: at $\chi_c=\chi_c^*$, the linearization of the traveling-wave equations around the resting disk has a one-dimensional kernel spanned by the translation mode, and the Crandall–Rabinowitz transversality condition holds, so a branch of even, fixed-area perturbation shapes exists. Along the branch the velocity satisfies $V(s)=s+o(s)$, the shape deformation starts at second order, and in the moving frame the pressure and marker fields are explicit: $P=p_1-Vx$ and $c=c_1e^{-aVx}$. The same analysis identifies the stabilizing role of undercooling: the term $\chi_u f_{\mathrm{und}}(V_n)$ raises the threshold by $\chi_u f_{\mathrm{und}}'(0)/(R_0 a c_0 f_{\mathrm{act}}'(c_0))$ relative to the model without membrane friction, so the membrane delays the onset of motion rather than preventing it.

Load-bearing premise

The argument presupposes that the resting disk, with zero velocity and zero added pressure, is an exact solution of the bifurcation functional for every marker strength, but substituting it into (4.3) leaves the nonzero marker force $\chi_c f_{\mathrm{act}}(c_0)$ unaccounted for unless that force is shifted into the pressure.

Editorial extensions

If this is right

  • For $\chi_c>\chi_c^*$ the resting disk is linearly unstable, giving a concrete symmetry-breaking mechanism for the onset of cell polarization.
  • A fixed-area traveling-wave family exists for all $a\in(0,1]$, $\gamma>0$, $R_0>0$ and $\chi_u>0$, so the membrane undercooling term does not eliminate the motility produced by the active marker force; it only shifts the threshold.
  • In the moving frame $P$ and $c$ are explicit, $P=p_1-Vx$ and $c=c_1e^{-aVx}$, so the full existence problem reduces to solving the scalar curvature equation (1.7) on the boundary.
  • The threshold formula separates the roles of parameters: larger membrane friction raises the required marker activity, while larger marker adsorption $a$ lowers the threshold.

Reading between the lines

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

  • The sign of $\chi_c''(0)$ in Lemma 4.3 controls which side of $\chi_c^*$ the branch lies on; because it depends on the second and third derivatives of $f_{\mathrm{act}}$ and the third derivative of $f_{\mathrm{und}}$, evaluating it for natural saturating choices such as $f_{\mathrm{act}}(u)=\tanh(u)$ would test the asserted parametrization over $(\chi_c^*,\infty)$.
  • A natural next step, not pursued in the paper, is to use $M$, $R_0$ or $a$ as the bifurcation parameter in the same Crandall–Rabinowitz setup; this would yield analogous thresholds and connect the onset of motion to marker mass or cell size.
  • In the formal limit $a\to0$ the threshold $\chi_c^*$ diverges, suggesting that a nonzero coupling between marker advection and boundary motion is necessary for the instability; the restriction $a\in(0,1]$ in the theorem may be more than technical.
  • Stability of the bifurcating traveling-wave branch is not addressed; applying the same spectral machinery along the branch would determine whether the predicted persistent motion is observable in a time-dependent simulation.
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 studies a two-dimensional incompressible Darcy free-boundary model for cell motility in which the boundary condition involves an undercooling membrane term χ_u f_und(V_n) and an active polarity-marker term χ_c f_act(c). After deriving the traveling-wave formulation and identifying a unique radial resting state, the authors perform a linear stability analysis (Theorem 1.5) and then use a Crandall-Rabinowitz bifurcation argument (Section 4) to claim a one-parameter family of non-trivial traveling waves for every χ_c > χ_c^* (Theorem 1.6).

Significance. If correct, the paper would extend the existence theory for Darcy-type cell motility models beyond the χ_u = 0 case in Alazard et al. (2022), giving a rigorous example of spontaneous persistent motion with membrane undercooling. The stability threshold χ_c^* and the stabilizing role of χ_u are concrete and physically meaningful. The paper also contains useful independent computations, including the traveling-wave characterization in Proposition 1.2 and the spectral decomposition in Section 3. However, the current proof of the main existence theorem is invalid as written: the claimed trivial branch is not a solution of the bifurcation functional, and the global-in-χ statement is not supported by the local bifurcation argument. The linear stability theorem also relies on an explicitly admitted unproved extension.

major comments (3)
  1. [§4, Eq. (4.3), Lemma 4.1(1)] Lemma 4.1(1) is false as stated. Substituting ρ = 0, V = 0, p1 = 0 into F in (4.3) gives F_1 = γ/R0 + χ_c fact(c0) + χ_u fund(0) + 0 - 0 - γ/R0 = χ_c fact(c0), since κ(0) = 1/R0, c1(0,0) = M/(πR0^2) = c0, and fund(0) = 0. For M > 0 this does not vanish because fact is increasing with fact(0) = 0. Thus F(χ_c,0,0,0) ≠ 0, so the point (χ_c^*,0,0,0) is not on a trivial solution branch and the Crandall-Rabinowitz theorem cannot be applied at that point. The actual resting state in Proposition 1.4 has p1 = γ/R0 + χ_c fact(c0), not p1 = 0. After the necessary shift p1 ↦ p1 + χ_c fact(c0), the parametrization (4.9) and the expansion in Lemma 4.2 change: the term χ'_c(0) fact(c0) that the authors use to force χ'_c(0) = 0 cancels against the shifted constant, so the pitchfork-direction conclusion is not established. The proof of Theorem 1.6 therefore has no valid base point as written.
  2. [§4, Theorem 1.6] The theorem claims a one-parameter family for all χ in (χ_c^*, +∞), but the proof applies the Crandall-Rabinowitz theorem, which yields a local curve near (χ_c^*,0,0,0) for a parameter s in some interval (-ε,ε). This gives χ_c(s) only in a neighborhood of χ_c^*. No continuation argument, global bifurcation theorem, or a priori bound is provided to extend the branch to the entire supercritical half-line. The statement of Theorem 1.6 is therefore not supported by the proof; at most a local bifurcation statement could follow if the base point issue in Lemma 4.1 were repaired.
  3. [§3, Remark 3.5 and Theorem 1.5] Theorem 1.5 asserts linear stability for all χ_c < χ_c^*, but Lemma 3.4 proves non-positive real parts of the eigenvalues only under the stronger restriction χ_c ≤ 1/(a c0 f'_act(c0)). Remark 3.5 explicitly states that the extension to all χ_c < χ_c^* is not proved and is only suggested by the later bifurcation study. Since χ_c^* = (R0 + χ_u f'_und(0))/(R0 a c0 f'_act(c0)) is strictly larger than 1/(a c0 f'_act(c0)) when χ_u > 0, the stability half of Theorem 1.5 rests on an admitted unproved assertion. The authors should either prove the extension or adjust the theorem to state only the interval for which the proof is valid.
minor comments (4)
  1. [§2, proof of Proposition 1.2] In the displayed formula for c1, the integrand contains a spurious factor c1: it should read c1 = M / ∫_Ω e^{-aV x'} dx' dy'. The current expression is dimensionally inconsistent and appears to be a typographical error.
  2. [§3, proof of Lemma 3.4] In the computation after equation (3.5), the notation f'_act(˜c) appears once with a tilde; this should be f'_act(c0) to match the rest of the argument.
  3. [§4, Lemma 4.3] The word 'Differenciating' should be 'Differentiating'. There are also several misspellings such as 'soution' in the proof of Lemma 4.1 item 3.
  4. [References] Alazard et al. (2022) is cited as an arXiv preprint throughout. If a peer-reviewed version has appeared, the reference should be updated to the published venue.

Circularity Check

0 steps flagged · score 1.0 of 10

No significant circularity: the bifurcation proof is self-contained, and the shared-author citation supplies independent computational details rather than the target result.

full rationale

The paper's central claim (Theorem 1.6) is a bifurcation existence result derived from the stated free-boundary problem (1.1) via the characterization in Proposition 1.2 and the Crandall-Rabinowitz theorem. No fitted parameter is renamed as a prediction: all quantities (chi_c^*, the bifurcation curve, the threshold) are computed from model parameters. The only self-citation with overlapping authorship is Alazard et al. (2022), used for the shape of the argument and for appendix C integral identities in Lemma 4.3; those identities are concrete, verifiable computations and do not assume the existence of traveling waves for chi_u > 0, so under the review rule they are independent support and do not raise the circularity score. The reader-identified issue that F(chi_c,0,0,0) = chi_c fact(c0) rather than 0 at p1 = 0 (Lemma 4.1(1)) is an internal base-point/correctness concern, not a circularity: it does not reduce Theorem 1.6 to its own inputs by definition or by a self-citation chain. There is no evidence of self-definitional reasoning, fitted parameters called predictions, or imported uniqueness theorems.

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

The central claims rest on the model assumptions, on the bifurcation setup, and on an unproved stability extension. No numbers are fitted to data; the parameters and functions are model inputs. The unproved assertions are recorded in the axioms and red flags.

assumptions (3)
  • domain assumption Model assumptions (1.2) and (1.3) on the functions fact and fund, constant total marker mass M, and fixed domain area.
    These assumptions define the class of models for which Theorems 1.5 and 1.6 are stated; they are inputs from the biological modeling context rather than derived results.
  • ad hoc to paper The bifurcation functional F admits a trivial branch F(chi_c,0,0,0)=0 for every chi_c, and the local bifurcation curve extends to cover all chi>chi_c^*.
    Lemma 4.1 asserts the trivial branch, but direct substitution into (4.3) contradicts it. The global range in Theorem 1.6 is also beyond the local Crandall-Rabinowitz conclusion.
  • ad hoc to paper Linear stability of the resting state holds for the full interval chi_c<chi_c^*, not only for the proved subinterval chi_c <= 1/(a c0 f'_act(c0)).
    Remark 3.5 explicitly states that the full result is not proved, yet Proposition 3.8 item 1 and Theorem 1.5 use it.

how reviews work

0 comments
Cite this review

Pith. "Pith review of Existence of traveling wave for a coupled incompressible Darcy's free boundary model with undercooling effect and surface tension." pith.science (2026). https://pith.science/paper/R43SUBBT

@misc{pith2026250104576,
  author       = {Pith},
  title        = {Pith review of: Existence of traveling wave for a coupled incompressible Darcy's free boundary model with undercooling effect and surface tension},
  year         = {2026},
  howpublished = {\url{https://pith.science/paper/R43SUBBT}},
  note         = {Machine review of arXiv:2501.04576}
}
read the original abstract

In this paper, we present a cell motility model that takes into account the cell membrane effect. The model introduced is an incompressible Darcy free boundary problem. This model involves a nonlinear term in the boundary condition to model the action of the membrane. This term can be seen as a undercooling effect of the membrane on the cell. It also implies a destabilizing nonlinear term in the boundary condition, depending on polarity markers and modeling the active character of the cytoskeleton. First, we study the linear stability of the steady state and prove that above a threshold, the disk is linearly unstable. This analysis highlights the stabilizing effect of undercooling. Then, using a bifurcation argument, we prove the existence of traveling waves that describe a persistent motion in cell migration and justify the relevance of the model.

Discussion (0). Continue with ORCID to comment.

Forward citations

Cited by 1 Pith paper

Reviewed papers in the Pith corpus that reference this work. Sorted by Pith novelty score. Full citation record

  1. Change of bifurcation type in 2D free boundary model of a moving cell with nonlinear diffusion

    math.AP 2025-06 conditional novelty 7.0 of 10

    A 2D free boundary cell motility model with nonlinear diffusion yields an explicit curvature formula K2 that is claimed to determine whether the pitchfork bifurcation to traveling waves is direct or inverse.

Reference graph

Works this paper leans on

28 extracted references · 27 canonical work pages · cited by 1 Pith paper

  1. [1]

    Abercrombie, M. (1980). The Croonian Lecture, 1978 - The cra wling movement of metazoan cells. Proceedings of the Royal Society of London. Series B. Biological Sciences, 207(1167):129–147

  2. [2]

    Alazard, T., Magliocca, M., and Meunier, N. (2022). Traveli ng wave solution for a coupled incompressible Darcy’s free boundary problem with surface tension. arXiv:2205.04365 [math]

  3. [3]

    M., Theriot, J

    Barnhart, E., Lee, K.-C., Allen, G. M., Theriot, J. A., and Mog ilner, A. (2015). Balance between cell-substrate adhesion and myosin contraction determine s the frequency of motility initiation in fish keratocytes. Proceedings of the National Academy of Sciences, 112(16):5045–5050

  4. [4]

    Berlyand, L., Fuhrmann, J., and Rybalko, V. (2018). Bifurcati on of traveling waves in a keller- segel type free boundary model of cell motility. Communications in Mathematical Sciences, 16(3):735–762

  5. [5]

    Berlyand, L., Potomkin, M., and Rybalko, V. (2016). Phase-fie ld model of cell motility: Traveling waves and sharp interface limit. Comptes Rendus. Mathématique, 354(10):986–992

  6. [6]

    Berlyand, L., Potomkin, M., and Rybalko, V. (2017). Sharp int erface limit in a phase field model of cell motility. Networks and Heterogeneous Media, 12(4):551–590

  7. [7]

    and Rybalko, V

    Berlyand, L. and Rybalko, V. (2019). Stability of steady stat es and bifurcation to traveling waves in a free boundary model of cell motility

  8. [8]

    Crandall, M. G. and Rabinowitz, P. H. (1971). Bifurcation fro m simple eigenvalues. Journal of Functional Analysis, 8(2):321–340

Show all 28 references
  1. [9]

    Cucchi, A., Mellet, A., and Meunier, N. (2020). A cahn–hilli ard model for cell motility. SIAM Journal on Mathematical Analysis, 52(4):3843–3880. 22

  2. [10]

    Cucchi, A., Mellet, A., and Meunier, N. (2022). Self polariz ation and traveling wave in a model for cell crawling migration. Discrete and Continuous Dynamical Systems, 42(5):2381–2407

  3. [11]

    Danuser, G., Allard, J., and Mogilner, A. (2013). Mathemati cal Modeling of Eukaryotic Cell Migration: Insights Beyond Experiments. Annual Review of Cell and Developmental Biology, 29(1):501–528

  4. [12]

    Etchegaray, C., Meunier, N., and Voituriez, R. (2017). Anal ysis of a Nonlocal and Nonlinear Fokker–Planck Model for Cell Crawling Migration. SIAM Journal on Applied Mathematics, 77(6):2040–2065

  5. [13]

    Friedman, A. (2004). A hierarchy of cancer models and their m athematical challenges. Discrete and Continuous Dynamical Systems - B, 4(1):147–159

  6. [14]

    and Hu, B

    Friedman, A. and Hu, B. (2006). Asymptotic stability for a fre e boundary problem arising in a tumor model. Journal of Differential Equations, 227(2):598–639

  7. [15]

    and Reitich, F

    Friedman, A. and Reitich, F. (2001). Symmetry-breaking bif urcation of analytic solutions to free boundary problems: An application to a model of tumor growth . Trans. Amer. Math. Soc., 353:1587–1634. Günther, M. and Prokert, G. (2008). On travelling-wave solu tions for a movi...

  8. [16]

    M., Barnhart, E

    Keren, K., Pincus, Z., Allen, G. M., Barnhart, E. L., Marriott , G., Mogilner, A., and Theriot, J. A. (2008). Mechanism of shape determination in motile cel ls. Nature, 453(7194):475–480

  9. [17]

    Lavi, I., Meunier, N., and Pantz, O. (2023). Implicit like ti me discretization for the one-phase Hele-Shaw problem with surface tension. arXiv:2305.06180 [cs, math]

  10. [18]

    Mizuhara, M. S. and Zhang, P. (2019). Uniqueness and traveli ng waves in a cell motility model. Discrete and Continuous Dynamical Systems - B, 24(6):2811–2835

  11. [19]

    Muller, N., Piel, M., Calvez, V., Voituriez, R., Gonçalves- Sá, J., Guo, C.-L., Jiang, X., Murray, A., and Meunier, N. (2016). A predictive model for yeast cell polarization in pheromone gradients. PLOS Computational Biology, 12(4):1–21

  12. [20]

    Putelat, T., Recho, P., and Truskinovsky, L. (2018). Mechan ical stress as a regulator of cell motility. Phys. Rev. E, 97:012410

  13. [21]

    Recho, P., Putelat, T., and Truskinovsky, L. (2013). Contra ction-driven cell motility. Phys. Rev. Lett., 111:108102. 23

  14. [22]

    Recho, P., Putelat, T., and Truskinovsky, L. (2015). Mechan ics of motility initiation and motility arrest in crawling cells. Journal of the Mechanics and Physics of Solids, 84:469–505

  15. [23]

    and Truskinovsky, L

    Recho, P. and Truskinovsky, L. (2013). Asymmetry between pu shing and pulling for crawling cells. Phys. Rev. E, 87:022720

  16. [24]

    Reig, G., Pulgar, E., and Concha, M. L. (2014). Cell migratio n: from tissue culture to embryos. Development, 141(10):1999–2013

  17. [25]

    and Berlyand, L

    Rybalko, V. and Berlyand, L. (2022). Emergence of traveling w aves and their stability in a free boundary model of cell motility. Trans. Amer. Math. Soc., 376 (2023):1799–1844

  18. [26]

    Shao, D., Levine, H., and Rappel, W.-J. (2012). Coupling act in flow, adhesion, and morphology in a computational cell motility model. Proceedings of the National Academy of Sciences, 109(18):6851–6856

  19. [27]

    Thiam, H.-R. (2014). Cell migration under confinement: how c an a cell squeeze through narrow gaps?

  20. [28]

    and Aranson, I

    Ziebert, F. and Aranson, I. S. (2016). Computational approa ches to substrate-based cell motility. npj Computational Materials, 2(1):16019. 24

Pith tools

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