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

TL;DR

This paper analyzes a 2D free boundary model of a moving cell with nonlinear diffusion, proving how nonlinearity causes a switch in bifurcation type and deriving formulas to predict this change based on physical parameters.

Contribution

It develops a rigorous framework for bifurcation analysis in a 2D free boundary PDE system with nonlinear diffusion and derives explicit formulas for bifurcation type change.

Findings

Nonlinear diffusion causes a switch between bifurcation types.

Explicit formulas relate bifurcation change to physical parameters.

Results align with numerical and physical literature observations.

Abstract

We introduce a 2D free boundary problem with nonlinear diffusion that models a living cell moving on a substrate. We prove that this nonlinearity results in a qualitative of solution behavior compared to the linear diffusion case (Rybalko et al. TAMS 2023), namely the switch between direct and inverse pitchfork bifurcation. Our objectives are twofold: (i) develop a rigorous framework to prove existence of bifurcation and determining its type (subcritical vs. superctitical) and (ii) the derivation of explicit analytical formulas that control the change of bifurcation type in terms of physical parameters and explain the underlying biophysical mechanisms. While the standard way of applying the Crandall-Rabinowitz theorem via the solution operator seems difficult in our quasilinear PDE system, we apply the theorem directly, by developing a multidimensional, vectorial framework. To determine the bifurcation type, we extract the curvature of the bifurcating curve from the expansion of the solutions around the steady state. The formula for the curvature is obtained via a solvability condition where instead of the Fredholm alternative, we propose a test function trick, suited for free boundary problems. Our rigorous analytical results are in agreement with numerical observations from the physical literature in 1D (Drozdowski et al. Comm. Phys. 2023) and provide the first extension of this phenomenon to a 2D free boundary model.