Here is a really basic question which I wished I understood better about the primes of the Gaussian field $\mathbb{Z}[i]$. But I was curious about the possibility of generalizing it to other (real quadratic or higher degree) number fields, see below.
Consider the primes $p < X$ which split $\mathbb{Q}(i)$, i.e., the ones $p \equiv 1 \mod{4}$. Then $p = a^2+b^2$ in an essentially unique way: one can switch the signs of $a$ or $b$ or interchange them, giving a total of $8$ points $(a,b)$ on the circle $|z| = \sqrt{p}$ in the Gaussian plane. This corresponds to looking at the two primes $(a+ib)$ and $(a-ib)$ above $p$ in $\mathbb{Z}[i]$, and noting that the unit group is just $\{\pm 1, \pm i\}$. Mark all these $8$ points (so that there is no ambiguity), for all the split primes $p < X$. As $X$ grows, how do these marked points distribute in angular sectors?
This sounds like Sato-Tate involving a curve $y^2 = x^3-x$ with CM by our Gaussian ring, particularly since $a \pm i b$ is a Weil number of weight one. But so also is its reflection $-i(a \pm ib) = b \mp ia$, and the two are not conjugate, hence would correspond to different elliptic isogeny classes / $\mathbb{F}_p$ in Honda's construction. (And I am confused.) Question: What is the relation, and how to solve in the most direct way possible the elementary problem here about sums of two squares?
And the second part of my question, which was also my motivation, concern the possibility of extending this to general number fields $K/\mathbb{Q}$. The prime ideal theorem describes the distribution of the norms of primes of $K$. But these prime ideals sit as lattices in the real vector space $K \otimes_{\mathbb{Q}} \mathbb{R}$, with points representing specific algebraic numbers (elements) from $K$. One could ask about the distribution of other qualities of these lattices than just their covolume (which gives essentially the norm $N (\mathfrak{p})$, as their index in the integer lattice $O_K$). As before, for each prime with $N(\mathfrak{p}) < X$, we could mark the first $\deg{K} = \dim(K \otimes_{\mathbb{Q}} \mathbb{R})$ shortest vectors in the corresponding lattice (if there are repetitions, mark them all), and ask how the marked points distribute as $X \to \infty$. The primes of degree $> 1$ make a negligible contribution, and so for $\mathbb{Q}(i)$ the question reduces to the above elementary problem.
Has anything been written on this type of question (the distribution of lattices of prime ideals, rather than just their covolumes), or is this question uninteresting for some reason that I do not see?