Let $p$ be a prime number and $n$ be positive integer. Let $E_{p,n}:y^2=x^3+p^nx$ be an elliptic curve.
LMFDB reads in the case $(p,n)=(73,3)$ , $\#Sha(E_{p,n})=64$. This is the biggest size of $Sha(E_{p,n})$ I have observed ever. Largest Tate-Shafarevich group with j-invariant $1728$ in LMFDB is this curve.
My question is, is it possible to find more large $Sha(E_{p,n})$?
I searched more than 100 primes with the following code of Magma, but the largest one was above.
A:=EllipticCurve([0,0,0,p^n,0]); MordellWeilShaInformation(A: ShaInfo:=true);
Thank you for your assistance. Even a single example without any background would be also appreciated. I believe the Tate-Shafarevich group for this family might not be well-documented, but if you're aware of any references for this curve's family, I'd greatly appreciate it.