skip to main content
US FlagAn official website of the United States government
dot gov icon
Official websites use .gov
A .gov website belongs to an official government organization in the United States.
https lock icon
Secure .gov websites use HTTPS
A lock ( lock ) or https:// means you've safely connected to the .gov website. Share sensitive information only on official, secure websites.


Title: Multiscale Convergence of the Inverse Problem for Chemotaxis in the Bayesian Setting
Chemotaxis describes the movement of an organism, such as single or multi-cellular organisms and bacteria, in response to a chemical stimulus. Two widely used models to describe the phenomenon are the celebrated Keller–Segel equation and a chemotaxis kinetic equation. These two equations describe the organism’s movement at the macro- and mesoscopic level, respectively, and are asymptotically equivalent in the parabolic regime. The way in which the organism responds to a chemical stimulus is embedded in the diffusion/advection coefficients of the Keller–Segel equation or the turning kernel of the chemotaxis kinetic equation. Experiments are conducted to measure the time dynamics of the organisms’ population level movement when reacting to certain stimulation. From this, one infers the chemotaxis response, which constitutes an inverse problem. In this paper, we discuss the relation between both the macro- and mesoscopic inverse problems, each of which is associated with two different forward models. The discussion is presented in the Bayesian framework, where the posterior distribution of the turning kernel of the organism population is sought. We prove the asymptotic equivalence of the two posterior distributions.  more » « less
Award ID(s):
1750488
PAR ID:
10332545
Author(s) / Creator(s):
; ; ;
Date Published:
Journal Name:
Computation
Volume:
9
Issue:
11
ISSN:
2079-3197
Page Range / eLocation ID:
119
Format(s):
Medium: X
Sponsoring Org:
National Science Foundation
More Like this
  1. 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
  2. null (Ed.)
    We consider a hyperbolic–parabolic system arising from a chemotaxis model in tumor angiogenesis, which is described by a Keller–Segel equation with singular sensitivity. It is known to allow viscous shocks (so-called traveling waves). We introduce a relative entropy of the system, which can capture how close a solution at a given time is to a given shock wave in almost [Formula: see text]-sense. When the shock strength is small enough, we show the functional is non-increasing in time for any large initial perturbation. The contraction property holds independently of the strength of the diffusion. 
    more » « less
  3. D. M. Kilgour et al. (Ed.)
    We consider a Keller-Segel type chemotaxis model with logarithmic sensitivity and density-dependent production/consumption rate. It is a 2 × 2 reaction-diffusion system describing the interaction of cells and a chemical signal. We study Cauchy problem for the original system and its transformed system, which is one of hyperbolic-parabolic conservation laws. In both cases of diffusive and non-diffusive chemical,we obtain optimal L^2 time decay rates for the solution. Our results improve those in Li et al. (Nonlinearity 28:2181-2210, 2015 [5]), Martinez et al. (Indiana Univ Math J 67:1383-1424, 2018 [7]). 
    more » « less
  4. We consider a [Formula: see text] system of hyperbolic balance laws that is the converted form under inverse Hopf–Cole transformation of a Keller–Segel type chemotaxis model. We study Cauchy problem when Cauchy data connect two different end-states as [Formula: see text]. The background wave is a diffusive contact wave of the reduced system. We establish global existence of solution and study the time asymptotic behavior. In the special case where the cellular population initially approaches its stable equilibrium value as [Formula: see text], we obtain nonlinear stability of the diffusive contact wave under smallness assumption. In the general case where the population initially does not approach to its stable equilibrium value at least at one of the far fields, we use a correction function in the time asymptotic ansatz, and show that the population approaches logistically to its stable equilibrium value. Our result shows two significant differences when comparing to Euler equations with damping. The first one is the existence of a secondary wave in the time asymptotic ansatz. This implies that our solutions converge to the diffusive contact wave slower than those of Euler equations with damping. The second one is that the correction function logistically grows rather than exponentially decays. 
    more » « less
  5. We consider a model of congestion dynamics with chemotaxis, where the density of cells follows the chemical signal it generates, while observing an incompressibility constraint (incompressible parabolic-elliptic Patlak-Keller-Segel model). We show that when the chemical diffuses slowly and attracts the cells strongly, then the dynamics of the congested cells is well approximated by a surface-tension driven free boundary problem. More precisely, we rigorously establish the convergence of the solution to the characteristic function of a set whose evolution is determined by the classical Hele-Shaw free boundary problem with surface tension. The problem is set in a bounded domain, which leads to an interesting analysis on the limiting boundary conditions. Namely, we prove that the assumption of Robin boundary conditions for the chemical potential leads to a contact angle condition for the free interface (in particular Neumann boundary conditions lead to an orthogonal contact angle condition, while Dirichlet boundary conditions lead to a tangential contact angle condition). 
    more » « less