Numerical investigation of lower order biases in moment expansions of one parameter families of elliptic curves

For a fixed elliptic curve E without complex multiplication, $a_p≔p+1-\#E\left(F_p\right)$ is $O\left(\sqrt{p}\right)$ and $a_p/2\sqrt{p}$ converges to a semicircular distribution. Michel proved that for a one-parameter family of elliptic curves $y^2=x^3+A\left(T\right)x+B\left(T\right)$ with $A\left(T\right),B\left(T\right)\in Z\left[T\right]$ and non-constant j-invariant, the second moment of $a_p\left(T\right)$ is $p^2+O\left(p^{3/2}\right)$. The size and sign of the lower order terms has applications to the distribution of zeros near the central point of Hasse-Weil L-functions and the Birch and Swinnerton-Dyer conjecture. S. J. Miller conjectured that the highest order term of the lower order terms of the second moment that does not average to zero is on average negative. Previous work on the conjecture has been restricted to a small set of highly nongeneric families. We create a database and a framework to quickly and systematically investigate biases in the second moment of any one-parameter family. When looking at families which have so far been beyond current theory, we find several potential violations of the conjecture for $p\le 250,000$ and discuss new conjectures motivated by the data.