Universality class of replica symmetry breaking, scaling behavior, and the low-temperature fixed-point order function of the Sherrington-Kirkpatrick model

Abstract

A scaling theory of replica symmetry breaking (RSB) in the Sherrington-Kirkpatrick (SK) model is presented in the framework of critical phenomena for the scaling regime of large RSB orders kappa , small temperatures T , and small (homogeneous) magnetic fields H . We employ the pseudodynamical picture [R. Oppermann, M. J. Schmidt, and D. Sherrington, Phys. Rev. Lett. 98, 127201 (2007)], where two critical points CP1 and CP2 are associated with the order function's pseudodynamical limits lim_{a-->infinity}q(a)=1 and lim_{a-->0}q(a)=0 at (T=0 , H=0 , 1kappa=0) . CP1 - and CP2 -dominated contributions to the free energy functional F[q(a)] require an unconventional scaling hypothesis. We determine the scaling contributions in accordance with detailed numerical self-consistent solutions for up to 200 orders of RSB. Power laws, scaling functions, and crossover lines are obtained. CP1 -dominated behavior is found for the nonequilibrium susceptibility, which decays like chi_{1}=kappa;{-53}f_{1}(Tkappa;{-53}) , for the entropy, which obeys S(T=0) approximately chi_{1};{2} , and for the subclass of diverging parameters a_{i}=kappa;{53}f_{a_{i}}(Tkappa;{-53}) [describing Parisi box sizes m_{i}(T) identical witha_{i}(T)T ], with f_{1}(zeta) approximately zeta and f_{a_{i}}(zeta) approximately 1zeta for zeta-->infinity , while f(0) is finite. CP2 -dominated behavior, controlled by the magnetic field H while temperature is irrelevant, is retrieved in the plateau height (or width) of the order function q(a) according to q_{pl}(H)=kappa;{-1}f_{pl}(H;{23}kappa;{-1}) with f_{pl}mid R:(zeta)mid R:_{zeta-->infinity} approximately zeta and f_{pl}(0) finite. Divergent characteristic RSB orders kappa_{CP1}(T) approximately T;{-35} and kappa_{CP2}(H) approximately H;{-23} , respectively, describe the crossover from mean field SK- to RSB-critical behavior with rational-valued exponents extracted with high precision from our RSB data. The order function q(a) is obtained as a fixed-point function q(a) of RSB flow, in agreement with integrated fixed-point energy and susceptibility distributions.