Multiple patterns formation for an aggregation/diffusion predator-prey system

We investigate existence of stationary solutions to an aggregation/diffusion system of PDEs, modelling a two species predator-prey interaction. In the model this interaction is described by non-local potentials that are mutually proportional by a negative constant \begin{document}$ -\alpha $\end{document} , with \begin{document}$ \alpha>0 $\end{document} . Each species is also subject to non-local self-attraction forces together with quadratic diffusion effects. The competition between the aforementioned mechanisms produce a rich asymptotic behavior, namely the formation of steady states that are composed of multiple bumps, i.e. sums of Barenblatt-type profiles. The existence of such stationary states, under some conditions on the positions of the bumps and the proportionality constant \begin{document}$ \alpha $\end{document} , is showed for small diffusion, by using the functional version of the Implicit Function Theorem. We complement our results with some numerical simulations, that suggest a large variety in the possible strategies the two species use in order to interact each other.