Zero-temperature TAP equations for the Ghatak-Sherrington model

The zero-temperature TAP equations for the spin-1 Ghatak-Sherrington model are investigated. The spin-glass energy density (ground state) is determined as a function of the anisotropy crystal field $D$ for a large number of spins. This allows us to locate a first-order transition between the spin-glass and paramagnetic phases within a good accuracy. The total number of solutions is also determined as a function of $D$.