--- title: Attraction-Repulsion Chemotaxis Model url: https://www.emergentmind.com/topics/attraction-repulsion-chemotaxis-model type: topic --- # Attraction-Repulsion Chemotaxis Model Attraction–repulsion chemotaxis models are PDE systems for a cell density \(u(x,t)\) coupled to at least two chemical fields, typically an attractant \(v(x,t)\) and a repellent \(w(x,t)\), in which directed motion toward higher \(v\) and away from higher \(w\) competes with diffusion, reaction kinetics, and, in many formulations, logistic damping or nonlinear production and consumption. A canonical form is \[ u_t=\Delta u-\chi \nabla\!\cdot(u\nabla v)+\xi \nabla\!\cdot(u\nabla w), \] with \(\chi,\xi>0\), where the term \(-\chi \nabla\!\cdot(u\nabla v)\) induces aggregation and \(+\xi \nabla\!\cdot(u\nabla w)\) counteracts it. Across the literature, the same attraction–repulsion architecture appears in fully parabolic, parabolic–elliptic, quasilinear, nonlocal, logistic, fractional, and fluid-coupled variants, and the central analytical question is whether the competition among attraction, repulsion, diffusion, and damping yields global boundedness, stabilization, or finite-time blow-up [2103.02241]. ## 1. Canonical formulations and modeling conventions The basic attraction–repulsion chemotaxis system is posed on a domain \(\Omega\subset\mathbb{R}^n\), usually bounded and smooth, with homogeneous Neumann boundary conditions \[ \partial_\nu u=\partial_\nu v=\partial_\nu w=0, \] so that diffusive and chemotactic fluxes are confined within \(\Omega\). In the standard local formulation, \(u\) denotes cell density, \(v\) a chemoattractant, and \(w\) a chemorepellent. The simplest linear-diffusion version takes the form \[ u_t=\Delta u-\chi \nabla\!\cdot(u\nabla v)+\xi \nabla\!\cdot(u\nabla w), \] supplemented by either parabolic signal equations such as \[ v_t=\Delta v-v+u,\qquad w_t=\Delta w-w+u, \] or elliptic signal equations such as \[ 0=\Delta v+\alpha u-\beta v,\qquad 0=\Delta w+\gamma u-\delta w. \] These two regimes are commonly described as fully parabolic and parabolic–elliptic, respectively [2103.02241]. A broad class of nonlinear generalizations replaces linear diffusion and linear sensitivities by density-dependent coefficients. One representative quasilinear system is \[ u_t=\nabla\cdot \big((u+1)^{m-1}\nabla u -\chi u(u+1)^{p-2}\nabla v +\xi u(u+1)^{q-2}\nabla w\big), \] with elliptic equations for \(v\) and \(w\). Here \((u+1)^{m-1}\) models density-dependent diffusion, while \(u(u+1)^{p-2}\) and \(u(u+1)^{q-2}\) encode density-dependent sensitivities to attraction and repulsion [2203.10240]. Another fully nonlinear local model includes a source \[ \lambda u^\rho-\mu u^k-c|\nabla u|^\gamma, \] together with nonlinear productions \(f_1(u)\) and \(f_2(u)\) for the chemical fields [2405.03586]. The chemical kinetics vary substantially across applications. In some models both signals are produced by the cells; in others both are consumed; and in mixed production–consumption systems the chemoattractant is produced while the chemorepellent is consumed, or conversely [2211.03608]. In the “double saturation” setting, \[ v_t=\Delta v-uv,\qquad w_t=\Delta w-uw, \] both chemicals are consumed, and the saturation terminology refers to the fact that the concentrations remain uniformly bounded by the maximum principle, not to sensitivity saturation [2106.07038]. Nonlocal variants replace local productions by mean-subtracted source terms, so that \(v\) and \(w\) represent deviations from their spatial averages and satisfy zero-mean constraints [2208.05679]. ## 2. Competition between attraction and repulsion The defining structural feature of the model is the competition between the attractive flux and the repulsive flux. In the most common sign convention, attraction is encoded by \[ -\chi \nabla\!\cdot(u\nabla v), \] which pulls cells up the attractant gradient, whereas repulsion is encoded by \[ +\xi \nabla\!\cdot(u\nabla w), \] which pushes cells down the repellent gradient. This competition is often quantified not by \(\chi\) and \(\xi\) alone, but by effective combinations involving production and decay parameters of the chemicals [2103.17044]. In linear production models with \[ \tau v_t=\Delta v+\alpha u-\beta v,\qquad \tau w_t=\Delta w+\gamma u-\delta w, \] the quantity \[ \chi\alpha-\xi\gamma \] is the natural measure of attraction versus repulsion. The regime \(\chi\alpha-\xi\gamma>0\) is described as attraction-dominated, while \(\chi\alpha-\xi\gamma<0\) is repulsion-dominated [2103.17044]. In the parabolic–parabolic and parabolic–elliptic systems studied with measure-valued initial data, the transformed coupling \[ \zeta:=\xi\gamma-\chi\alpha \] is used instead, so that “strong repulsion” corresponds to \(\zeta>0\) [2210.12208]. A recurrent misconception is that repulsion alone precludes singularity formation. The literature does not support that general statement. Finite-time blow-up persists in attraction-dominated fully parabolic systems in three dimensions, even though repulsion is present [2103.17044]. Conversely, sufficiently strong repulsion, or repulsion combined with damping, can ensure global boundedness and stabilization in many regimes [2203.10240]. This suggests that the analytically relevant issue is not the mere presence of a repellent, but the quantitative balance among attraction, repulsion, diffusion, production, consumption, and logistic or gradient-dependent damping. The same balance appears in one-dimensional half-line models with logistic growth. There the effective attraction and repulsion strengths are \[ \widehat{\chi}_1=\frac{\chi_1\mu_1}{\lambda_1},\qquad \widehat{\chi}_2=\frac{\chi_2\mu_2}{\lambda_2}, \] and global boundedness is obtained under conditions of the form \[ b_{\inf}>\widehat{\chi}_1-\widehat{\chi}_2+M, \] where \(M\) is a resolvent-based constant measuring the mismatch of the two elliptic signal operators [1812.09860]. ## 3. Reduction to effective Keller–Segel structures One of the most important analytical devices in attraction–repulsion chemotaxis is the reduction of the two-signal system to a single effective signal. In the symmetric fully parabolic model \[ u_t=\Delta u-\chi\nabla\!\cdot(u\nabla v)+\xi\nabla\!\cdot(u\nabla w),\qquad v_t=\Delta v-v+u,\qquad w_t=\Delta w-w+u, \] the linear transformation \[ z:=\chi v-\xi w \] yields \[ u_t=\Delta u-\nabla\!\cdot(u\nabla z),\qquad z_t=\Delta z-z+(\chi-\xi)u. \] This is exactly a parabolic–parabolic Keller–Segel-like system with effective coupling \(\chi-\xi\) [2103.02241]. The structural significance of this reduction is twofold. First, it converts the competing chemotactic drifts into the single advection structure \(-\nabla\!\cdot(u\nabla z)\). Second, it imports the classical Keller–Segel energy machinery. The canonical functional becomes \[ \mathcal{F}(u,z):=\frac{1}{2}\!\int_\Omega |\nabla z|^2+\frac{1}{2}\!\int_\Omega z^2-(\chi-\xi)\!\int_\Omega uz+(\chi-\xi)\!\int_\Omega u\ln u, \] and satisfies the dissipation inequality \[ \frac{d}{dt}\,\mathcal{F}(u(\cdot,t),z(\cdot,t))\le -\mathcal{D}(u(\cdot,t),z(\cdot,t)), \] with \[ \mathcal{D}(u,z):=\int_\Omega z_t^2+(\chi-\xi)\!\int_\Omega u\,|\nabla \ln u-\nabla z|^2. \] Under \(\chi>\xi\), this reproduces the higher-dimensional Keller–Segel blow-up mechanism in the attraction–repulsion setting [2103.02241]. A related transformation is used in the three-dimensional attraction-dominated fully parabolic system with unequal damping and production parameters. Writing \[ w := (\chi\alpha-\xi\gamma)u,\qquad z:=\chi v_1-\xi v_2, \] the system becomes \[ w_t=\Delta w-\nabla\!\cdot(w\nabla z),\qquad z_t=\Delta z-az+bv+w,\qquad v_t=\Delta v-cv+dw. \] When \(b=0\), this is the classical parabolic–parabolic Keller–Segel structure; when \(b\neq 0\), the exact Lyapunov structure is broken, but an energy method remains possible [2103.17044]. The reduction idea also appears in repulsion-dominant quasilinear parabolic–elliptic–elliptic systems, but there it no longer closes because the nonlinear sensitivity factors cannot be written as a single \(u\nabla z\) drift. The stabilization analysis therefore uses a tailored entropy functional rather than a direct Keller–Segel reduction [2203.10240]. ## 4. Global boundedness mechanisms The boundedness theory is organized around several distinct mechanisms. One is the weakening of signal growth through consumption. In the double-saturation model \[ v_t=\Delta v-uv,\qquad w_t=\Delta w-uw, \] the maximum principle yields \[ 0\le v(x,t)\le \|v_0\|_{L^\infty(\Omega)},\qquad 0\le w(x,t)\le \|w_0\|_{L^\infty(\Omega)}, \] and a weighted functional \[ \mathcal{E}_k(t)=\int_\Omega u^k e^{\beta^2 v^2+\gamma^2 w^2} \] provides uniform-in-time \(L^k\) bounds on \(u\) when \[ 0<\chi<\frac{1}{5n\|v_0\|_{L^\infty(\Omega)}},\qquad 0<\xi<\frac{1}{5n\|w_0\|_{L^\infty(\Omega)}}. \] These \(L^k\) bounds with \(k>n/2\) imply \(L^\infty\)-boundedness by heat-semigroup smoothing and bootstrap [2106.07038]. Another mechanism is sublinear or weaker-than-competing production. In the parabolic–elliptic–elliptic system \[ 0=\Delta v+\alpha u^p-\beta v,\qquad 0=\Delta w+\gamma u-\delta w, \] with \(00\), any \(0\tfrac12 n(\alpha-1)\quad\text{and}\quad \ell\le \min\{\alpha-1,\rho\}, \] or \[ \beta>\tfrac12\big(n\ell+2(\ell-\alpha+1)\big)\quad\text{and}\quad \alpha-1<\min\{\ell,\rho\}. \] For fully parabolic signals, the corresponding thresholds involve \(\max\{\rho,\ell\}\) and reflect the added difficulty of time-dependent signal dynamics [2602.14227]. Further boundedness mechanisms include density-dependent diffusion, repulsion-dominant sensitivity growth, and gradient-dependent damping. In the local and nonlocal fully nonlinear systems with source \[ \lambda u^\rho-\mu u^k-c|\nabla u|^\gamma, \] global boundedness follows under \[ 1\le \rho\xi\). More precisely, for any \(m>0\) and \(A>0\), there exist \(T=T(m,A)>0\) and \(K=K(m,A)>0\) such that if the initial data belong to the class \(\mathcal{C}(m,A)\) with negative energy \[ \mathcal{G}(u_0,v_0,w_0)\le -K, \] then the maximal existence time satisfies \(T_{\max}\le T\) and \[ \lim_{t\nearrow T_{\max}}\|u(\cdot,t)\|_{L^\infty(\Omega)}=\infty. \] This closes the previously open \(n\ge 4\) gap in the symmetric fully parabolic attraction–repulsion setting [2103.02241]. In the three-dimensional fully parabolic attraction-dominated system \[ u_t=\Delta u-\chi\nabla\!\cdot(u\nabla v_1)+\xi\nabla\!\cdot(u\nabla v_2),\quad \partial_t v_1=\Delta v_1-\beta v_1+\alpha u,\quad \partial_t v_2=\Delta v_2-\delta v_2+\gamma u, \] posed on a ball with Neumann boundary conditions, finite-time blow-up occurs whenever \[ \chi\alpha-\xi\gamma>0. \] The result is stated by approximation: for every radial positive triple \((u_0,v_{1,0},v_{2,0})\) and every \(p\in(1,3/2)\), there exist arbitrarily close radial initial data in \(L^p\times L^2\times L^2\) with the same cell mass whose corresponding classical solutions blow up in finite time [2103.17044]. The blow-up analysis in fully parabolic systems also yields explicit lower bounds on blow-up time. For \[ u_t=\Delta u-\chi \nabla \cdot (u \nabla v)+\xi \nabla \cdot (u \nabla w)+g(u),\quad v_t=\Delta v-\alpha v+\beta u,\quad w_t=\Delta w-\gamma w+\delta u, \] with \[ g(u)=\mu_1 u-\mu_2u^k, \] the energy \[ E_{p,q}(t)=\frac1p\int_\Omega u^p+\frac1q\int_\Omega |\nabla v|^q+\frac1q\int_\Omega |\nabla w|^q \] satisfies a superlinear ODI, leading to the explicit estimate \[ T_{\max} \geq \int_{E_{p,q}(0)}^{\infty} \left ( m \sum_{i=0}^3 s^{\eta_i}+ \sum_{i=0}^3 m_{i} s^{k(\eta_i)}+\mu_1 s+ c \right )^{-1} \, ds, \] for admissible \(p,q,s_1,s_2\). This result does not establish blow-up, but quantifies the minimal blow-up time when blow-up occurs [2203.06251]. ## 6. Extensions, regularization, and broader variants The attraction–repulsion architecture extends well beyond the basic local Keller–Segel system. One notable direction concerns regularization from singular initial data. In the two-dimensional parabolic–parabolic system and the two- and three-dimensional parabolic–elliptic system, \[ u_0\in\mathcal{M}^+(\Omega) \] may be a positive Radon measure, including sums of Dirac masses. Under sufficiently strong repulsion, classical solutions become instantaneously smooth. The precise conditions are \[ \zeta=\xi\gamma-\chi\alpha>0 \] for the parabolic–parabolic system in \(n=2\), \[ \zeta>-C_\Omega \] for the parabolic–elliptic system in \(n=2\), and \[ \zeta\ge 0,\qquad u_0\in L^K(\Omega)\ \text{for some }K\in(1,2) \] for the parabolic–elliptic system in \(n=3\). The solution satisfies \[ \int_\Omega u(x,t)\phi(x)\,dx\to \int_\Omega \phi(x)\,du_0(x) \] for every \(\phi\in C^0(\Omega)\) as \(t\downarrow 0\) [2210.12208]. Another direction adds fractional diffusion and whole-space propagation. In \[ u_t=-(-\Delta)^\alpha u-\chi_1\nabla\!\cdot(u\nabla v)+\chi_2\nabla\!\cdot(u\nabla w)+au-bu^\gamma, \] with elliptic signals \[ 0=\Delta v-\lambda_1 v+\mu_1u^k,\qquad 0=\Delta w-\lambda_2 w+\mu_2u^k, \] global boundedness follows under conditions involving \[ M := \min \left\{ \frac{(\chi_2 \mu_2 \lambda_2 - \chi_1 \mu_1 \lambda_1)_+ + \chi_2 \mu_2 (\lambda_1 - \lambda_2)_+}{\lambda_1},\ \frac{(\chi_2 \mu_2 \lambda_2 - \chi_1 \mu_1 \lambda_1)_+ + \chi_1 \mu_1 (\lambda_1 - \lambda_2)_+}{\lambda_2} \right\}, \] and, under balanced intensities \[ \chi_2\mu_2=\chi_1\mu_1,\qquad \lambda_1=\lambda_2, \] the exact spreading speed is \[ c^*=\frac{a}{N+2\alpha}. \] This reflects the accelerated propagation induced by fractional diffusion [2603.26148]. Attraction–repulsion models also couple to fluid motion. In the chemotaxis–Navier–Stokes system \[ n_t+u\cdot\nabla n=\Delta n-\chi\nabla\cdot(n\nabla c)+\xi\nabla\cdot(n\nabla v)+\varsigma n-\mu n^2, \] with transport equations for the chemical signals and incompressible fluid velocity, global mild solutions exist on bounded domains in two and three dimensions for small initial data in \(L^p\)-spaces [1709.04032]. At a different scale, a field-theoretic model for run-and-tumble particles uses a density field \(\rho\), polarization field \(P\), and chemical field \(c\), with chemotaxis entering through a biased tumbling rate \[ \alpha(\hat{\mathbf{u}})=\tilde{\alpha}_0-\chi(c,\tilde v)\nabla c\cdot \hat{\mathbf{u}}, \] thereby recovering Keller–Segel-type drift in an appropriate closure while retaining particle-resolution, elasticity, and soft-core repulsion [2011.04580]. Across these variants, a coherent picture emerges. Repulsion can regularize measure-valued data, enlarge boundedness regimes, or determine spreading behavior, but it does not universally eliminate blow-up. A plausible implication is that attraction–repulsion chemotaxis should be understood not as a single model with one threshold, but as a family of competing-drift systems whose long-time behavior depends on which structural mechanism—effective attraction, repulsion, sublinear production, signal consumption, nonlinear diffusion, logistic damping, nonlocal feedback, or gradient dissipation—dominates in the relevant scaling regime [2210.12208]. Source: https://www.emergentmind.com/topics/attraction-repulsion-chemotaxis-model