Radial symmetry of positive solutions to quasilinear Hardy–Sobolev doubly critical systems

The aim of this paper is to prove radial symmetry results for positive weak solutions with finite energy to the following quasilinear doubly critical system
$$\left\{ \begin{array}{lcc}
-\Delta_p u = \gamma \frac{u^{p-1}}{|x|^p}+u^{p^*-1}+\nu\alpha u^{\alpha-1}v^{\beta} & \text{in} & \mathbb{R}^n \\
-\Delta_p v = \gamma \frac{v^{p-1}}{|x|^p}+v^{p^*-1}+\nu\beta u^{\alpha}v^{\beta-1} & \text{in} & \mathbb{R}^n \end{array}\right. $$
where $1<p<n$, $\gamma \in [0,\Lambda_{n,p})$ with $\Lambda_{n,p} = [(n−p)/p]^p$, $\alpha, \beta > 1$ such that $\alpha+\beta=p^*= np/(n−p)$ and $\nu>0$.