On the critical Schrödinger-Poisson system with $ p $-Laplacian

In this paper we consider the critical quasilinear Schrödinger-Poisson system\begin{document}$ \begin{eqnarray*} \left \{\begin{array}{ll} -\Delta_p u+|u|^{p-2}u+\mu\phi |u|^{p-2}u = \lambda|u|^{q-2}u+|u|^{p^*-2}u,&\mathrm{in} \ \mathbb{R}^3,\\ -\Delta \phi = |u|^p, &\mathrm{in}\ \mathbb{R}^3, \end{array} \right. \end{eqnarray*} $\end{document}where \begin{document}$ \frac{3}{2}<p<3 $\end{document} , \begin{document}$ \Delta_p u = \hbox{div}(|\nabla u|^{p-2}\nabla u) $\end{document} , \begin{document}$ p<q<p^*: = \frac{3p}{3-p} $\end{document} and \begin{document}$ \mu,\lambda>0 $\end{document} . Based upon the variational approach, the ground state solutions and the nontrivial solutions are obtained depending on the parameters \begin{document}$ q $\end{document} , \begin{document}$ \mu $\end{document} and \begin{document}$ \lambda $\end{document} .