Critical blow-up curve in a two-species chemotaxis system with two chemicals involving flux-limitation

We investigate the following two-species chemotaxis system with two chemicals involving flux-limitation \begin{align}\tag{$\star$} \begin{cases} u_t = Δu - \nabla \cdot \left(u(1+|\nabla v|^2)^{-\frac{p}{2}}\nabla v\right), & x \in Ω, \ t > 0, \\ 0 = Δv - μ_w + w, \quad μ_{w}=f_Ω w, & x \in Ω, \ t > 0, \\ w_t = Δw - \nabla \cdot \left(w (1+|\nabla z|^2)^{-\frac{q}{2}} \nabla z\right), & x \in Ω, \ t > 0, \\ 0 = Δz - μ_u + u, \quad μ_{u}=f_Ω u, & x \in Ω, \ t > 0, \\ \frac{\partial u}{\partial ν} = \frac{\partial v}{\partial ν} = \frac{\partial w}{\partial ν} = \frac{\partial z}{\partial ν} = 0, & x \in \partial Ω, \ t > 0, \\ u(x, 0) = u_0(x), \quad w(x, 0) = w_0(x), & x \in Ω, \end{cases} \end{align} where $p,q \in \mathbb{R}$ and $Ω\subset \mathbb{R}^n$ is a smooth bounded domain. In this paper, we identify a critical blow-up curve for system ($\star$) with $n\geq 3$. If $p<\frac{n-2}{n-1}$ and $q<\frac{n-2}{n-1}$, and $Ω=B_R(0) \subset \mathbb{R}^n$ with $n\geq 3$, there exist radially symmetric initial data such that the corresponding solution blows up in finite time; if either $p>\frac{n-2}{n-1}$ or $q>\frac{n-2}{n-1}$ with $n\geq 2$, then solutions exist globally and remain bounded.