Quasilinear Elliptic Cooperative and Competitive Systems

We study the existence and multiplicity of weak solutions for the following quasilinear elliptic system: $$\{\begin{array}{ll}-\mathrm{d}\mathrm{i}\mathrm{v}(A_1(x,u_1)\mathrm{\nabla }{u}_1)+{\displaystyle \frac{1}{2}{D}_{u_1}{A}_1(x,u_1)\mathrm{\nabla }{u}_1\cdot \mathrm{\nabla }{u}_1={\lambda }_1{u}_1+g_{\beta ,1}(u)}& \text{in }\mathrm{\Omega },\\ -\mathrm{d}\mathrm{i}\mathrm{v}(A_2(x,u_2)\mathrm{\nabla }{u}_2)+{\displaystyle \frac{1}{2}{D}_{u_2}{A}_2(x,u_2)\mathrm{\nabla }{u}_2\cdot \mathrm{\nabla }{u}_2={\lambda }_2{u}_2+g_{\beta ,2}(u)}& \text{in }\mathrm{\Omega },\\ u_1=u_2=0& \text{on }\mathrm{\partial }\mathrm{\Omega },\end{array}$$ where ${\lambda }_1,{\lambda }_2<{\mu }_1$, the first Dirichlet eigenvalue of the Laplacian, and $\mathrm{\Omega }$ is a bounded domain. The nonlinearity derives from a potential $G_{\beta }$ with subcritical growth. Due to the lack of differentiability of the associated energy functional, we employ nonsmooth critical point theory and variational methods based on the concept of weak slope. We prove the existence of least energy solutions in both the cooperative ($\beta >0$) and competitive ($\beta <0$) regimes.