Segregated solutions for a class of systems with Lotka-Volterra interaction

This paper deals with the existence of positive solutions to the system $$-\mathrm{\Delta }{w}_1-\epsilon w_1={\mu }_1{w}_1^p+\beta w_1{w}_2\text{ in }\mathrm{\Omega },-\mathrm{\Delta }{w}_2-\epsilon w_2={\mu }_2{w}_2^p+\beta w_1{w}_2\text{ in }\mathrm{\Omega },w_1=w_2=0\text{ on }\mathrm{\partial }\mathrm{\Omega },$$ where $\mathrm{\Omega }\subseteq {\mathbb{R}}^N$, $N\ge 4$, $p=\frac{N+2}{N-2}$ and $\epsilon$ is positive and sufficiently small. The interaction coefficient $\beta =\beta (\epsilon )\to 0$ as $\epsilon \to 0$. We construct a family of segregated solutions to this system, where each component blows-up at a different critical point of the Robin function as $\varepsilon \to 0. The system lacks a variational formulation due to its specific coupling form, which leads to essentially different behaviors in the subcritical, critical, and supercritical regimes and requires an appropriate functional settings to carry out the construction.