In this paper we study positive solutions to the following nonlocal system of equations: \begin{equation*} \left\{\begin{aligned} &(-\Delta)^s u = \frac{\alpha}{2_s^*}|u|^{\alpha-2}u|v|^{\beta}+f(x)\;\;\text{in}\;\mathbb{R}^{N}, v \frac{\beta}{2_s^*}|v|^{\beta-2}v|u|^{\alpha}+g(x)\;\;\text{in}\;\mathbb{R}^{N}, & \qquad u, \, >0\, \mbox{ in }\,\mathbb{R}^{N}, \end{aligned} \right. \end{equation*...