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 →
The pith
A machine-rendered reading of the paper's core claim, the machinery that carries it, and where it could break.
The reading
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.
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
- 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.
Signed reviews
Editorial analysis
A structured set of objections, weighed in public.
Referee Report
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)
- [§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.
- [§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, 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)
- [§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.
- [§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.
- [§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.
- [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
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
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.
- 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^*.
- 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)).
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.
Forward citations
Cited by 1 Pith paper
-
Change of bifurcation type in 2D free boundary model of a moving cell with nonlinear diffusion
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
-
[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
work page 1980
-
[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]
arXiv 2022
-
[3]
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
work page 2015
-
[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
work page 2018
-
[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
work page 2016
-
[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
work page 2017
-
[7]
Berlyand, L. and Rybalko, V. (2019). Stability of steady stat es and bifurcation to traveling waves in a free boundary model of cell motility
work page 2019
-
[8]
Crandall, M. G. and Rabinowitz, P. H. (1971). Bifurcation fro m simple eigenvalues. Journal of Functional Analysis, 8(2):321–340
work page 1971
Show all 28 references
-
[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
2020
-
[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
2022
-
[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
2013
-
[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
2017
-
[13]
Friedman, A. (2004). A hierarchy of cancer models and their m athematical challenges. Discrete and Continuous Dynamical Systems - B, 4(1):147–159
2004
-
[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
2006
-
[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...
2001
-
[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
2008
-
[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]
2023 arXiv
-
[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
2019
-
[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
2016
-
[20]
Putelat, T., Recho, P., and Truskinovsky, L. (2018). Mechan ical stress as a regulator of cell motility. Phys. Rev. E, 97:012410
2018
-
[21]
Recho, P., Putelat, T., and Truskinovsky, L. (2013). Contra ction-driven cell motility. Phys. Rev. Lett., 111:108102. 23
2013
-
[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
2015
-
[23]
and Truskinovsky, L
Recho, P. and Truskinovsky, L. (2013). Asymmetry between pu shing and pulling for crawling cells. Phys. Rev. E, 87:022720
2013
-
[24]
Reig, G., Pulgar, E., and Concha, M. L. (2014). Cell migratio n: from tissue culture to embryos. Development, 141(10):1999–2013
2014
-
[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
2022
-
[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
2012
-
[27]
Thiam, H.-R. (2014). Cell migration under confinement: how c an a cell squeeze through narrow gaps?
2014
-
[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
2016
Reviewed August 10, 2026 · model on record in the stance chip above.
Discussion (0). Continue with ORCID to comment.