I have come across the following surface: let $X$ be the double covering of $\mathbb{P}_\mathbb{Z}^2$ defined by the equation $$y^2=x_0^6+x_1^6+x_2^6$$ where $y$ is a variable of degree 3.
There is an obvious action of $\mu_6 \times \mu_6$.
Less obvious, I have made some short calculations that show a great link between this surface over $\mathbb{F}_p$ and $E^2$ over $\mathbb{F}_p$, for any prime $p$, where $E$ is the well-known CM elliptic curve $$y^2=x^3+1$$
Without getting too much into specifics, the surface's Brauer group seems to give rise to a modular representation coming from a weight $3$ cusp form, that seems to be a modular form also attached to $E^2$ (in the same way), maybe up to a character (of degree 6). This is just some very simple mod $p$ calculations and strong multiplicity one. But the question isn't about the validity of these calculations, so I don't want you to dwell on this.
The described link makes me guess that the surface is geometrically (and not only arithmetically) related to the squared elliptic curve. Are they birationally equivalent? Maybe after dividing $E^2$ by some subgroup of $\mu_6 \times \mu_6$? How does one begin to check this in sage?