{"id":"7b9d3f49-7a5e-4a3c-8be6-636488767b50","arxiv_id":"2506.03138","paper_version":2,"verdict":"CONDITIONAL","confidence":"MODERATE","novelty_score":7.0,"correctness_risk":"medium","formal_verification":"none","parameter_count":0,"one_line_summary":"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.","lead":"This math paper derives a formula for the curvature of the traveling-wave bifurcation curve in a 2D free boundary model of a crawling cell with nonlinear myosin diffusion. The formula is meant to predict when the onset of motion switches from a direct to an inverse pitchfork, a qualitative change previously seen numerically only in 1D.","discovery_kind":"first_principles","skeptic_critique":{"model":"deepseek-v4-flash","headline":"Corollary 1's advertised direct/inverse switch rests on an undefined '1D version' of (19) and uncomputed coefficients A_i; the eA* = 0.5990 value and sign pattern are unsupported, and non-degeneracy (16) is unchecked for the stated parameters.","rationale":"I read the paper as aiming to prove both a 2D bifurcation-type formula and a concrete prediction of the direct/inverse switch for the van der Waals myosin model. The 2D development is substantive: the direct Crandall-Rabinowitz application, the explicit first-order eigenfunction, the test-function method for K2, and the structurally derived D-dependence are nontrivial and could well be correct. However, the theorem's hypothesis (16) and the second-order systems are not carried through to a numerical evaluation, and the corollary switches to 1D without derivation. Since the abstract advertises agreement with 1D numerics and the first extension to 2D, the unsupported corollary is the most load-bearing part: the public claim is the switch at eA* = 0.5990, not merely the existence of a formula. The reader's weakest assumption identifies the same gap. This is a gap of missing derivation and computation, not an internal inconsistency in the 2D argument, so the honest verdict remains conditional: accept after the 1D reduction and coefficient evaluations are supplied, or explicitly narrow the claims. I do not see a reason to change the reader's conditional verdict.","tokens_in":27614,"tokens_out":7285,"duration_ms":87298,"concrete_test":"Ask the authors to supply the 1D analogue of (19) by repeating the expansion of Section 2 on an interval (or by proving a suitable reduction from the disk), then: (i) evaluate the left-hand side of (16) at P = 0.1, Z = 1.25, m_infinity = 10 and confirm it is nonzero; (ii) numerically solve the second-order systems in Lemma 4, compute A1..A4 from (38)-(41) for these parameters over eA in [0, 0.7], and (iii) re-solve K2(eA) = 0 and compare with eA* = 0.5990 and with the claimed sign pattern. If the recomputed eA* differs by more than 1%, or if any A_i changes sign, Corollary 1 should be withdrawn or reclassified as an unverified numerical observation.","verdict_should_be":"UNCHANGED","load_bearing_attack":"The load-bearing point is not the 2D formula (19) itself but its only physical application. Corollary 1 states a '1D version of (19)' with sign pattern A1, A3, A4 > 0, A2 < 0 and a critical value eA* = 0.5990, yet no 1D analogue of Theorem 2 or of Lemmas 3-5 is stated or proved. All expansion work in Appendices B-C is 2D: the unknowns are Fourier modes on the disk, the test function U(r) is the radial Bessel J1 solution of (26)/(169), and the boundary conditions involve curvature and Hanzawa coordinates. A 1D interval model has none of these structures, so the restriction of (19) to 1D is not a routine 'set cos(theta) terms aside.' Moreover, the advertised coefficients are not evaluated: (42) defines A_i as ratios of integrals involving f1..f4 in (161)-(164), and those f_i contain m20A, m20B, m22A, m22B, sigma20A, sigma22A, sigma22B from the second-order system in Lemma 4, which the paper does not solve. No computation producing eA* = 0.5990 or the asserted signs is shown. Finally, Theorem 2 and Corollary 1 both assume non-degeneracy (16); the paper never checks it for P = 0.1, Z = 1.25. If (16) fails there, or if the 1D coefficients change sign, the headline switch does not follow. This is exactly the unsupported link between the 2D theorem and the advertised physics.","agreement_with_reader":"agree"},"referee_report":{"model":"deepseek-v4-flash","summary":"The paper studies a two-dimensional free boundary model of cell motility with nonlinear myosin diffusion, extending the linear-diffusion model of Rybalko et al. It proves, via a direct application of the Crandall-Rabinowitz theorem to the quasilinear PDE system, the existence of traveling wave solutions bifurcating from the stationary radial solution, and derives an explicit formula for the curvature K2 of the bifurcating branch. This curvature determines whether the pitchfork bifurcation is direct or inverse. The authors then apply a claimed one-dimensional version of this formula to the van der Waals diffusion coefficient and predict a switch between direct and inverse bifurcation at eA ≈ 0.5990, matching numerical observations of Drozdowski et al. The paper's advertised contribution is the rigorous derivation of the 2D bifurcation-type formula and its physical consequence.","tokens_in":27931,"tokens_out":3956,"duration_ms":47132,"significance":"If the results are correct, the paper would provide a novel rigorous framework for the Crandall-Rabinowitz theorem in quasilinear free-boundary problems, bypassing the solution-operator construction, and would give the first mathematically derived formula for the change of pitchfork bifurcation type in a 2D moving-cell model. The proposed test-function trick for extracting the curvature without a full Fredholm alternative is an interesting methodological idea, and the explicit dependence of K2 on the nonlinear diffusion coefficient D(m) would be a useful quantitative tool. However, the advertised physical application rests on a one-dimensional reduction that is neither stated as a theorem nor proved, and the coefficients in the general formula are not evaluated. These gaps substantially reduce the current significance of the paper's main claims.","major_comments":[{"comment":"Corollary 1 asserts that for P=0.1, Z=1.25, m∞=10, the curvature K2(eA) is given by a '1D version of (19)' with A1, A3, A4 > 0 and A2 < 0, and that the bifurcation switches from direct to inverse at eA* = 0.5990. No one-dimensional analogue of Theorem 2, Lemma 3, Lemma 4, or Lemma 5 is stated or proved. The 2D derivation in Appendix B uses Fourier modes on the disk, the Bessel test function U(r) of (169), and boundary conditions involving curvature and the Hanzawa transform; it is not a routine restriction. Moreover, the coefficients A_i in (42) are not evaluated anywhere, the numerical value eA* = 0.5990 is not accompanied by any computation, and the non-degeneracy condition (16) is not checked for the stated parameters. Since the change-of-bifurcation-type claim is the paper's headline result, this unsupported link is load-bearing.","section":"§1.3, Corollary 1"},{"comment":"The proof that F_x(x0,K0) has one-dimensional kernel and codimension-one range is only sketched. After reducing to the system (58), the text asserts that 'from the ellipticity of the problem' one can solve for m in terms of σ and that the resulting problem (61) has a solution 'provided compatibility conditions that are linear in V', and that one can always choose V so that a solution exists except for a special case deferred to Step 3. This does not establish the claimed codim(Range(F_x))=1, which is needed for condition (iii) of the Crandall-Rabinowitz theorem. The argument also does not address the Fredholm index of F_x in the chosen spaces X and Y in (95)-(96). This is a central technical step and requires a complete proof.","section":"§3, Step 2 (Fredholm property and codimension-one range)"},{"comment":"The transversality argument rests on the explicit solution (70) and then asserts the derivative formula (71). The transition from (70) to (71) is not shown, and the formula involves the constant B whose definition is not made explicit in the displayed equations. Since the contradiction with the boundary condition (69d) is exactly what forces the non-degeneracy condition (16), this omitted computation is not cosmetic. The authors should provide the detailed differentiation and account for all terms, including the A(V) term, before invoking (16).","section":"§3, Eqs. (70)-(71)"},{"comment":"The claim that formula (19) is explicit is undercut by the fact that the functions f_i(r) in (161)-(164) depend on m20A, m20B, m22A, m22B, σ20A, σ22A, σ22B, which are introduced in Lemma 4 through the systems (149)-(150) but are never solved or even shown to exist uniquely. Consequently the 'explicitly given' coefficients A_i in (42) remain formal integrals over unknown functions. The proof of Lemma 4 in Appendix C stops after writing the systems (149)-(150); no closed-form or qualitative information about these second-order coefficients is provided. The sign pattern asserted in Corollary 1 therefore cannot be derived from the material in the paper.","section":"§2, Lemma 5 and Eq. (42)"}],"minor_comments":[{"comment":"There are several typos: 'BIFURCA TION' in the title, 'superctitical' in the abstract, and 'qualitative of solution behavior' should be 'qualitative change of solution behavior'.","section":"Title and Abstract"},{"comment":"For the van der Waals diffusion coefficient D(m) = m∞^2/(m∞−m)^2 − eA m, the paper does not specify the admissible range of m relative to m∞, nor does it state the value of m0 for the parameters used in Corollary 1. Since D becomes singular at m = m∞, this is needed to verify positivity and the differentiability assumptions of Theorem 2.","section":"§1.3, Eq. (20)"},{"comment":"Figure 3 plots K2(eA) but contains no data points or reproducible numerical values; it cannot serve as a substitute for the missing computation of the coefficients and the critical value eA*.","section":"Fig. 3"},{"comment":"The integration by parts leading to (171) produces boundary terms at r = 0; the paper does not justify their vanishing, e.g., by stating the regularity of σ31 near the origin. This should be made explicit.","section":"§2, Lemma 5"}],"recommendation":"major_revision","confidential_remarks":"The core 2D framework is potentially valuable, and the direct Crandall-Rabinowitz approach is interesting. However, the advertised physical result in Corollary 1 is not supported by the text: the 1D reduction is absent, the coefficients are not evaluated, and the non-degeneracy condition is unchecked. This is fixable in principle, but it requires substantial additional work, not merely a revision of wording. I recommend major revision, with the expectation that the authors either provide a complete 1D theorem and computation of eA*, or substantially weaken the claims about the van der Waals application."},"author_rebuttal":null,"desk_editor":{"model":"deepseek-v4-flash","letter":"The thing to know: the 2D machinery is genuinely new and worth engaging with, but the 1D switch that drives the abstract and the biophysical punchline is not, as written, a proved consequence. Applying Crandall-Rabinowitz directly to the PDE system, avoiding the solution operator that is so natural in the linear case, plus the test-function trick for the third-order solvability condition, is a real improvement over the earlier linear diffusion papers. The formula K2 = A1 D''/D^2 + ... with D-independent coefficients is exactly the kind of explicit output that makes this line useful. The paper is also honest about its structure: lemmas are stated, expansions are relegated to appendices, and the 1D reduction is not hidden—it just isn't there.\n\nSoft spots, in proportion. First, the Fredholm property and codimension-one range condition in Step 2 of Theorem 1 are asserted with a sketch rather than established; the line about 'from the ellipticity of the problem' does a lot of work, and a referee will want a real argument. Second, the jump from the explicit solution (70) to the derivative (71), and the claim that non-degeneracy (16) gives the transversality contradiction, is compressed to the point of being uncheckable. That is fixable but needs to be written out.\n\nThe serious problem is Corollary 1. It invokes a '1D version of (19)' but no 1D analogue of Theorem 2 or of Lemmas 3-5 is stated or proved. The 2D expansions are built on Fourier modes on the disk, Bessel test functions, Hanzawa coordinates, and curvature; none of these reduces by setting cos(theta) terms aside. The coefficients A_i are defined as ratios of integrals involving f1..f4, and those f_i depend on second-order quantities m20A, m20B, etc. that are never solved. No computation produces eA* = 0.5990 or the sign pattern A1, A3, A4 > 0, A2 < 0. Non-degeneracy (16) is never checked for P = 0.1, Z = 1.25. So the advertised switch—the thing that makes the paper matter to physicists—is unsupported. The 2D theorem may well be correct; the corollary needs either a full derivation or an explicit caveat that it is a conjecture/computational observation.\n\nThe paper deserves a serious referee and likely major revision. I would not cite the corollary in my own work until the 1D derivation or narrowed claims appear, but the 2D method itself is worth a look.","headline":"Novel direct Crandall-Rabinowitz framework for a quasilinear free boundary problem, but the advertised 1D pitchfork-switch corollary is currently an unproved claim.","tokens_in":28504,"tokens_out":2935,"would_cite":false,"duration_ms":35228,"reading_group":"maybe","serious_thinker":"yes","would_accept_peer_review":true},"rs_alignment":null,"lean_confirmation":null,"pith_extraction":{"msc":["35B32","70K50","92C17"],"pacs":[],"model":"deepseek-v4-flash","headline":"The paper proves that nonlinear diffusion can flip a 2D free-boundary cell-motility model from a direct to an inverse pitchfork bifurcation, and derives an explicit curvature formula that locates the switch.","keywords":["free boundary problem","nonlinear diffusion","pitchfork bifurcation","cell motility","traveling waves","Crandall-Rabinowitz theorem","van der Waals diffusion","active matter"],"falsifier":"A concrete check would be to numerically compute the traveling-wave branch of the full 2D system for the van der Waals diffusion coefficient at $P=0.1$, $Z=1.25$ across $e_A$ from 0.55 to 0.65 and verify whether the sign of the numerically extracted branch curvature at the bifurcation point crosses zero near $e_A \\approx 0.5990$; if it does not, the quantitative prediction of Corollary 1 fails.","tokens_in":27338,"feed_emoji":"🧫","tokens_out":6575,"duration_ms":73948,"temperature":0.7,"pith_summary":"This paper shows that adding nonlinear diffusion to a two-dimensional free-boundary model of a crawling cell changes the way motion begins. In the linear-diffusion case only a direct (supercritical) pitchfork bifurcation appears, but the authors prove that with nonlinear diffusion the bifurcation can become inverse (subcritical), allowing bistability between a resting and a moving cell. To establish this, they apply the Crandall-Rabinowitz theorem directly to the quasilinear PDE system and extract the curvature of the bifurcating branch via a test-function trick. A closed formula expresses this curvature in terms of the diffusion coefficient and its derivatives at the steady state, and a 1D application to a van der Waals diffusion coefficient predicts a switch at a critical parameter value, matching numerical observations.","feed_headline":"Curvature formula reveals switch in cell-motion bifurcation","feed_subtitle":"A new explicit K2 formula tells whether a crawling cell's onset is direct or inverse pitchfork.","key_machinery":"The protagonist is the curvature $K_2$ of the bifurcating curve of traveling waves, defined through the expansion $K(V) = K_0 + K_2 V^2 + \\dots$. The machinery is a third-order asymptotic expansion of the PDE system in the small speed $V$, combined with a specially constructed test function: a solution of the formal adjoint Helmholtz-type equation $Z\\Delta u + (\\alpha/R_0)^2 u = 0$ with boundary value $\\cos\\theta$ on the disk. Testing the overdetermined third-order system with this function turns the solvability condition into an explicit equation for $K_2$, avoiding the Fredholm alternative that is difficult to apply in the presence of the free boundary.","core_discovery":"The central discovery is that the type of pitchfork bifurcation in this 2D free-boundary cell-motility model is controlled by an explicit formula for the curvature of the traveling-wave branch. Expanding the Peclet number as $K = K_0 + K_2 V^2 + \\dots$, the paper proves that $K_2 = A_1 \\frac{D''(m_0)}{D(m_0)^2} + A_2 \\frac{D'(m_0)^2}{D(m_0)^3} + A_3 \\frac{D'(m_0)}{D(m_0)^2} + A_4 \\frac{1}{D(m_0)}$, where the coefficients $A_i$ depend only on the physical parameters $P, Z, \\gamma$, not on the diffusion function $D$. The sign of $K_2$ determines the bifurcation type: positive $K_2$ gives a direct (supercritical) pitchfork, negative $K_2$ gives an inverse (subcritical) pitchfork. The paper also proves existence of the traveling-wave branch through a direct application of the Crandall-Rabinowitz theorem under a stated non-degeneracy condition, and it applies a 1D analogue of the curvature formula to the van der Waals diffusion coefficient to predict a switch at $e_A^* \\approx 0.5990$.","pith_inferences":["Editorial inference: the coefficient signs stated for the 1D van der Waals case suggest a general heuristic that decreasing diffusion at the resting density ($D'(m_0) < 0$) is what drives the switch to an inverse pitchfork; testing this against other diffusion profiles would clarify how universal the mechanism is.","Editorial inference: the same curvature formula could be used to detect whether alternative bifurcation branches, such as higher-mode or asymmetric traveling waves, become preferred when $K_2$ changes sign, since the sign of the curvature controls the local geometry of the branch.","Editorial inference: a direct numerical computation of the 2D traveling-wave branch for the van der Waals coefficient would provide a check of whether the 1D critical value $e_A^* \\approx 0.5990$ persists in 2D, which the paper does not explicitly compute."],"forward_implications":["Nonlinear myosin diffusion can produce a subcritical pitchfork bifurcation, meaning the onset of cell motion can be abrupt and can coexist with a stable resting state, a behavior impossible in the linear-diffusion model.","The explicit formula for $K_2$ lets one predict the direct-versus-inverse character of the bifurcation from the shape of the diffusion coefficient $D(m)$ and the physical parameters $P, Z, \\gamma$, without solving the full nonlinear problem.","For the van der Waals diffusion coefficient with the stated parameters, the model predicts a sharp transition at $e_A^* \\approx 0.5990$, below which motion begins supercritically and above which it begins subcritically.","The test-function method provides a route for determining bifurcation type in free-boundary problems where the standard solution-operator approach and the Fredholm alternative are not available."],"supporting_citations":[{"why":"Supplies the Crandall-Rabinowitz theorem whose hypotheses the paper verifies directly on the PDE system.","marker":"[16]"},{"why":"Establishes the linear-diffusion case with only direct pitchfork bifurcation, providing the contrast that nonlinear diffusion overturns.","marker":"[33]"},{"why":"Provides the 1D van der Waals model and numerical observations that Corollary 1 is designed to match.","marker":"[20]"},{"why":"Gives the non-dimensionalization and the resting-state parameters used to set up the free-boundary model.","marker":"[31]"},{"why":"Introduces the nonlocal boundary condition and the linear stability framework that the traveling-wave analysis builds on.","marker":"[34]"}],"fun_headline_variants":["Curvature sign reveals cell-motion bifurcation type","Nonlinear diffusion flips 2D cell bifurcation","Explicit K2 formula fixes pitchfork type in cell model","Test-function trick yields bifurcation in 2D motility","Supercritical or subcritical: K2 formula decides"],"cache_read_input_tokens":3200,"weakest_assumption_plain":"The 1D application of the curvature formula rests on an unstated 1D analogue of the 2D formula (19) with a specific sign pattern for the coefficients $A_i$, and this reduction is not derived in the paper.","fun_headline_variants_meta":{"raw":{"variants":["Curvature sign reveals cell-motion bifurcation type","Nonlinear diffusion flips 2D cell bifurcation","Explicit K2 formula fixes pitchfork type in cell model","Test-function trick yields bifurcation in 2D motility","Supercritical or subcritical: K2 formula decides"]},"model":"deepseek-v4-flash","effort":"low","cost_usd":0.000211,"raw_usage":{"total_tokens":1489,"prompt_tokens":1093,"completion_tokens":396,"prompt_tokens_details":{"cached_tokens":384},"prompt_cache_hit_tokens":384,"prompt_cache_miss_tokens":709,"completion_tokens_details":{"reasoning_tokens":315}},"tokens_in":709,"tokens_out":396,"duration_ms":5111,"temperature":1.0,"reasoning_tokens":315,"cache_read_input_tokens":384,"cache_creation_input_tokens":0},"cache_creation_input_tokens":0},"created_at":"2026-08-07T11:08:31.575819+00:00","model_set":{"reader":"deepseek-v4-flash"},"falsifier":"A concrete check would be to numerically compute the traveling-wave branch of the full 2D system for the van der Waals diffusion coefficient at $P=0.1$, $Z=1.25$ across $e_A$ from 0.55 to 0.65 and verify whether the sign of the numerically extracted branch curvature at the bifurcation point crosses zero near $e_A \\approx 0.5990$; if it does not, the quantitative prediction of Corollary 1 fails.","supporting_citations":[{"cited_title":"Bifurcation from simple eigen- values","cited_arxiv_id":null,"evidence_quote":"Supplies the Crandall-Rabinowitz theorem whose hypotheses the paper verifies directly on the PDE system."},{"cited_title":"Emergence of traveling waves and their stability in a free boundary model of cell motility","cited_arxiv_id":null,"evidence_quote":"Establishes the linear-diffusion case with only direct pitchfork bifurcation, providing the contrast that nonlinear diffusion overturns."},{"cited_title":"Optogenetic control of migration of contractile cells predicted by an active gel model","cited_arxiv_id":null,"evidence_quote":"Provides the 1D van der Waals model and numerical observations that Corollary 1 is designed to match."},{"cited_title":"Mechanics of motility initiation and motility arrest in crawling cells","cited_arxiv_id":null,"evidence_quote":"Gives the non-dimensionalization and the resting-state parameters used to set up the free-boundary model."},{"cited_title":"Asymptotic sta- bility of contraction-driven cell motion","cited_arxiv_id":null,"evidence_quote":"Introduces the nonlocal boundary condition and the linear stability framework that the traveling-wave analysis builds on."}],"review_version":1}