Kirchhoff systems involving fractional p-Laplacian and singular nonlinearity

In this work we consider the fractional Kirchhoff equations with singular nonlinearity, $$\displaylines{ M\Big( \int_{\mathbb{R}^{2N}}\frac{|u(x)-u(y)|^p}{|x-y|^{N+sp}}dx dy\Big) (-\Delta)^s_p u = \lambda a(x)|u|^{q-2}u +\frac{1-\alpha}{2-\alpha-\beta} c(x)|u|^{-\alpha}|v|^{1-\beta}, \quad \text{in }\Omega,\cr \int_{\mathbb{R}^{2N}}\frac{|v(x)-v(y)|^p}{|x-y|^{N+sp}}dx v \mu b(x)|v|^{q-2}v +\frac{1-\beta}{2-\alpha-\beta} c(x)|u|^{1-\alpha}|v|^{-\beta}, }\Omega, \cr u=v 0 ,\quad\hbox{in }\mathbb{R}^N\setminus\Omega, }$$ where \(\Omega\) is a bounded domain in \(\mathbb{R}^N\) smooth boundary, \(N> ps\), \(s \in (0,1)\), \(0<\alpha<1\), \(0<\beta<1\), \(2-\alpha-\beta<p\leq p\theta<q<p^*_s\), \(p^*_s=\frac{Np}{N-sp}\) Sobolev exponent, \(\lambda, \mu\) are two parameters, \(a, b, c C(\overline{\Omega})\) non-negative weight functions, \(M(t)=k+lt^{\theta-1}\) \(k>0,l,\theta\geq 1\), and \((-\Delta)^s_p\) p-laplacian operator. We prove existence of multiple solutions by studying nature Nehari manifold respect to parameters \(\lambda\) \(\mu\).