Numerical Study of TAP Metastable States in 3-body Ising Spin Glasses

The distribution of solutions of the Thouless-Anderson-Palmer equation is studied by extensive numerical experiments for fully connected 3-body interaction Ising spin glass models in a level of annealed calculation. A recent study predicted that when the equilibrium state of the system is characterized by one-step replica symmetry breaking, the distribution is described by a Becchi-Rouet-Stora-Tyutin (BRST) supersymmetric solution in the relatively low free energy region, whereas the BRST supersymmetry is broken for higher values of free energy (Crisanti et al., Phys. Rev. B 71 (2005) 094202). Our experiments qualitatively reproduce the discriminative behavior of macroscopic variables predicted by the theoretical assessment.