In this work we consider the following fractional Kirchhoff equations with singular nonlinearity: M(∫R2N|u(x)-u(y)|2|x-y|N+2sdxdy)(-Δ)su=λa(x)|u|q-2u+1-α2-α-βc(x)|u|-α|v|1-β,inΩ,M(∫R2N|v(x)-v(y)|2|x-y|N+2sdxdy)(-Δ)sv=μb(x)|v|q-2v+1-β2-α-βc(x)|u|1-α|v|-β,inΩ,u=v=0,inRN\Ω,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\begin{aligned} \left\{ \begin{array}{ll} M\Big ( \int _{\mathbb {R}^{2N}}\frac{|u(x)-u(y)|^{2}}{|x-y|^{N+2s}}dx dy\Big )(-\Delta )^s 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 ,\\ M\Big ( \int _{\mathbb {R}^{2N}}\frac{|v(x)-v(y)|^{2}}{|x-y|^{N+2s}}dx dy\Big ) (-\Delta )^s v= \mu b(x)|v|^{q-2}v +\frac{1-\beta }{2-\alpha -\beta } c(x)|u|^{1-\alpha }|v|^{-\beta }, \quad \text {in }\Omega ,\\ u=v = 0,\;\; \text{ in } \,\mathbb {R}^N\setminus \Omega , \end{array} \right. \end{aligned}$$\end{document}where Ω\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\Omega $$\end{document} is a bounded domain in Rn\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mathbb {R}^n$$\end{document} with smooth boundary ∂Ω\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\partial \Omega $$\end{document}, N>2s\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$N> 2s$$\end{document}, s∈(0,1)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$s \in (0,1)$$\end{document}, 0<α<1,0<β<1,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0<\alpha<1,\;0<\beta <1,$$\end{document}0<α+β<2θ<q<2s∗,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$0<\alpha +\beta<2\theta<q<2^*_s,$$\end{document}2s∗=2NN-2s\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$2^*_s=\frac{2N}{N-2s}$$\end{document} is the fractional Sobolev exponent, λ,μ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda , \mu $$\end{document} are two parameters, a,b,c∈C(Ω¯)\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$a,\, b, \,c \in C({\overline{\Omega }})$$\end{document} are non-negative weight functions, M is a continuous function, given by M(t)=k+ltθ-1\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$M(t)=k+lt^{\theta -1}$$\end{document}k>0,l,θ≥1,\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$k>0,\,l,\,\theta \ge 1,$$\end{document} and (-Δ)s\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$(-\Delta )^s$$\end{document} is the fractional Laplacien operator. We use the Nehari manifold approach and some variational techniques in order to show the existence and multiplicity of positive solutions of the above problem with respect to the parameter λ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\lambda $$\end{document} and μ\documentclass[12pt]{minimal}
\usepackage{amsmath}
\usepackage{wasysym}
\usepackage{amsfonts}
\usepackage{amssymb}
\usepackage{amsbsy}
\usepackage{mathrsfs}
\usepackage{upgreek}
\setlength{\oddsidemargin}{-69pt}
\begin{document}$$\mu $$\end{document}.
