Global stabilization of the full attraction-repulsion Keller-Segel system

We are concerned with the following full Attraction-Repulsion Keller-Segel (ARKS) system \begin{equation}\label{ARKS}\tag{$\ast$} \begin{cases} u_t=\Delta u-\nabla\cdot(\chi u\nabla v)+\nabla\cdot(\xi u\nabla w), &x\in \Omega, ~~t>0, v_t=D_1\Delta v+\alpha u-\beta v,& x\in \Omega, ~~t>0, w_t=D_2\Delta w+\gamma u-\delta w, &x\in \Omega, ~~t>0,\ u(x,0)=u_0(x),~v(x,0)= v_0(x), w(x,0)= w_0(x) & x\in \Omega, \end{cases} \end{equation} in a bounded domain $\mathrm{\Omega }\subset {\mathbb{R}}^2$ with smooth boundary subject to homogeneous Neumann boundary conditions. %The parameters $D_1,D_2,\chi ,\xi ,\alpha ,\beta ,\gamma$ and $\delta$ are positive. By constructing an appropriate Lyapunov functions, we establish the boundedness and asymptotical behavior of solutions to the system \eqref{ARKS} with large initial data. Precisely, we show that if the parameters satisfy $\frac{\xi\gamma}{\chi\alpha}\geq \max\Big{\frac{D_1}{D_2},\frac{D_2}{D_1},\frac{\beta}{\delta},\frac{\delta}{\beta}\Big}$forallpositiveparameters$D_1,D_2,\chi,\xi,\alpha,\beta,\gamma$and$\delta$,thesystem\text{eqref}ARKShasauniqueglobalclassicalsolution$(u,v,w)$,whichconvergestotheconstantsteadystate$(\bar{u}_0,\frac{\alpha}{\beta}\bar{u}_0,\frac{\gamma}{\delta}\bar{u}_0)$as$t\to+\infty$,where$\bar{u}0=\frac{1}{|\Omega|}\int\Omega u_0dx$. Furthermore, the decay rate is exponential if $\frac{\xi\gamma}{\chi\alpha}> \max\Big{\frac{\beta}{\delta},\frac{\delta}{\beta}\Big}$. This paper provides the first results on the full ARKS system with unequal chemical diffusion rates (i.e. $D_1\mathrm{\ne }{D}_2$) in multi-dimensions.