In this paper, we study the nonexistence of positive stable weak solutions to the following problem in the whole space \(\mathbb R^N=\mathbb R^{N_1} \times \mathbb R^{N_2}\) \(\begin{aligned} \operatorname {div}_G (w_1(z)|\nabla _Gu|^{p-2}\nabla _Gu + w_2(z) |\nabla _Gu|^{q-2}\nabla _Gu)= f(z)u^{-r}, \end{aligned}\) where \(\Delta _G=\text {div}_G\circ \nabla _G\) is the Grushin operator, \(q\ge p \ge 2, r>0\) and \(w_1,w_2, f \in L^1_{{{\,\textrm{loc}\,}}}(\mathbb R^N;[0,\infty ))\) satisfy some appropriate conditions at infinity. In particular, our results can be seen as a generalization of some previous results to the double phase problem.
