An official website of the United States government
Abstract Chemotaxis phenomena govern the directed movement of microorganisms in response to chemical stimuli. In this paper, we investigate two Keller–Segel systems of reaction–advection–diffusion equations modeling chemotaxis on thin networks. The distinction between two systems is driven by the rate of diffusion of the chemo-attractant. The intermediate rate of diffusion is modeled by a coupled pair of parabolic equations, while the rapid rate is described by a parabolic equation coupled with an elliptic one. Assuming the polynomial rate of growth of the chemotaxis sensitivity coefficient, we prove local well-posedness of both systems on compact metric graphs, and, in particular, prove existence of unique classical solutions. This is achieved by constructing sufficiently regular mild solutions via analytic semigroup methods and combinatorial description of the heat kernel on metric graphs. The regularity of mild solutions is shown by applying abstract semigroup results to semi-linear parabolic equations on compact graphs. In addition, for logistic-type Keller–Segel systems we prove global well-posedness and, in some special cases, global uniform boundedness of solutions.
more »
« less
This content will become publicly available on November 5, 2025
Stability and bifurcation for logistic Keller-Segel models on compact graphs
This paper concerns asymptotic stability, instability, and bifurcation of constant steady state solutions of the parabolic-parabolic and parabolic-elliptic chemotaxis models on metric graphs. We determine a threshold value $$\chi^*>0$$ of the chemotaxis sensitivity parameter that separates the regimes of local asymptotic stability and instability, and, in addition, determine the parameter intervals that facilitate global asymptotic convergence of solutions with positive initial data to constant steady states. Moreover, we provide a sequence of bifurcation points for the chemotaxis sensitivity parameter that yields non-constant steady state solutions. In particular, we show that the first bifurcation point coincides with threshold value $$\chi^*$$ for a generic compact metric graph. Finally, we supply numerical computation of bifurcation points for several graphs.
more »
« less
- Award ID(s):
- 2418900
- PAR ID:
- 10621808
- Publisher / Repository:
- Yokohama Publ.
- Date Published:
- Journal Name:
- Pure and applied functional analysis
- Volume:
- 9
- Issue:
- 5
- ISSN:
- 2189-3764
- Format(s):
- Medium: X
- Sponsoring Org:
- National Science Foundation
More Like this
-
-
null (Ed.)Abstract We study steady-state thin films on chemically heterogeneous substrates of finite size, subject to no-flux boundary conditions. Based on the structure of the bifurcation diagram, we classify the 1D steady-state solutions that exist on such substrates into six different branches and develop asymptotic estimates for the steady states on each branch. Using perturbation expansions, we show that leading-order solutions provide good predictions of the steady-state thin films on stepwise-patterned substrates. We show how the analysis in one dimension can be extended to axisymmetric solutions. We also examine the influence of the wettability contrast of the substrate pattern on the linear stability of droplets and the time evolution for dewetting on small domains. Results are also applied to describe 2D droplets on hydrophilic square patches and striped regions used in microfluidic applications.more » « less
-
Abstract We consider steady state solutions of the massive, asymptotically flat, spherically symmetric Einstein–Vlasov system, i.e., relativistic models of galaxies or globular clusters, and steady state solutions of the Einstein–Euler system, i.e., relativistic models of stars. Such steady states are embedded into one-parameter families parameterized by their central redshift$$\kappa >0$$ . We prove their linear instability when$$\kappa $$ is sufficiently large, i.e., when they are strongly relativistic, and prove that the instability is driven by a growing mode. Our work confirms the scenario of dynamic instability proposed in the 1960s by Zel’dovich & Podurets (for the Einstein–Vlasov system) and by Harrison, Thorne, Wakano, & Wheeler (for the Einstein–Euler system). Our results are in sharp contrast to the corresponding non-relativistic, Newtonian setting. We carry out a careful analysis of the linearized dynamics around the above steady states and prove an exponential trichotomy result and the corresponding index theorems for the stable/unstable invariant spaces. Finally, in the case of the Einstein–Euler system we prove a rigorous version of the turning point principle which relates the stability of steady states along the one-parameter family to the winding points of the so-called mass-radius curve.more » « less
-
Abstract When chemotactic bacteria are exposed to a concentration gradient of chemoattractant while flowing along a channel, the bacteria accumulate at the interface between the chemoattractant source and bacterial suspension. Assuming that the interface is no‐slip, we can apply the shear flow approximation near the no‐slip boundary and solve a steady‐state convection‐diffusion model for both chemoattractant and bacterial concentrations. We suggest similarity solutions for the two‐dimensional problem and identify a critical length scaleηcfor bacteria chemotaxis in a given concentration gradient. The analysis identifies three dimensionless groups representing, respectively, chemotactic sensitivity, the chemotaxis receptor constant, and the bacteria diffusion coefficient, which typically show coupled effects in experimental systems. We study the effect of the dimensionless groups separately and provide understanding of the system involving shear flow and chemotaxis.more » « less
-
The construction of bifurcation diagrams is an essential component of understanding nonlinear dynamical systems. The task can be challenging when one knows the equations of the dynamical system and becomes much more difficult if only the underlying data associated with the system are available. In this work, we present a transformer-based method to directly estimate the bifurcation diagram using only noisy data associated with an arbitrary dynamical system. By splitting a bifurcation diagram into segments at bifurcation points, the transformer is trained to simultaneously predict how many segments are present and to minimize the loss with respect to the predicted position, shape, and asymptotic stability of each predicted segment. The trained model is shown, both quantitatively and qualitatively, to reliably estimate the structure of the bifurcation diagram for arbitrarily generated one- and two-dimensional systems experiencing a codimension-one bifurcation with as few as 30 trajectories. We show that the method is robust to noise in both the state variable and the system parameter.more » « less
