Quadratic units and cubic fields

We investigate Eisenstein discriminants, which are squarefree integers $d \equiv 5 \pmod{8}$ such that the fundamental unit $\varepsilon_d$ of the real quadratic field $K=\mathbb{Q}(\sqrt{d})$ satisfies $\varepsilon_d \equiv 1 \pmod{2\mathcal{O}_K}$. These discriminants are related to a classical question of Eisenstein and have connections to the class groups of orders in quadratic fields as well as to real cubic fields. We present numerical computations of Eisenstein discriminants up to $10^{11}$, suggesting that their counting function up to $x$ is approximated by $π_{\mathcal{E}}(x) \approx \frac{1}{3π^2}x - 0.024x^{5/6}$. This supports a conjecture of Stevenhagen while revealing a surprising secondary term, which is similar to (but subtly different from) the secondary term in the counting function of real cubic fields. We include technical details of our computation method, which uses a modified infrastructure approach implemented on GPUs.