P.G.Walsh proved in this paper that the diophantine equation $\,(2^{x_1}-1)(3^{x_2}-1)=y^2\,$ has no solution in positive integers $\,x_1$, $\,x_2\,$ and $\,y$.
If we generalize the previous equation to the following (in $\,n+1$ variables: $\,x_1,\,x_2,\,...,\,x_n,\,y$) $$\prod_{k=1}^n(p_k^{x_k}-1)=y^2$$ ($p_k$ being the prime of index $\,k$) positive integer solutions appear. For example: $$n=3\;\;\rightarrow\;\;(2,2,2)\;\;\;\;y=24$$ $$n=4\;\;\rightarrow\;\;(2,1,1,1)\;\;\;\;y=12$$ $$n=5\;\;\rightarrow\;\;(2,4,1,1,1)\;\;\;\;y=240$$ $$n=6\;\;\rightarrow\;\;(1,4,1,1,1,1)\;\;\;\;y=480$$ $$n=7\;\;\rightarrow\;\;(1,4,1,1,1,1,1)\;\;\;\;y=1920$$ Does the diophantine equation $$\prod_{k=1}^n(p_k^{x_k}-1)=y^2$$ have always (at least) a solution for $\,n\gt2\,$?
Many thanks.