No Arabic abstract
We argue that when the number of spins $N$ in the SK model is finite, the Parisi scheme can be terminated after $K$ replica-symmetry breaking steps, where $K(N) propto N^{1/6}$. We have checked this idea by Monte Carlo simulations: we expect the typical number of peaks and features $R$ in the (non-bond averaged) Parisi overlap function $P_J(q)$ to be of order $2K(N)$, and our counting (for samples of size $N$ up to 4096 spins) gives results which are consistent with our arguments. We can estimate the leading finite size correction for any thermodynamic quantity by finding its $K$ dependence in the Parisi scheme and then replacing $K$ by K(N). Our predictions of how the Edwards-Anderson order parameter and the internal energy of the system approach their thermodynamic limit compare well with the results of our Monte Carlo simulations. The $N$-dependence of the sample-to-sample fluctuations of thermodynamic quantities can also be obtained; the total internal energy should have sample-to-sample fluctuations of order $N^{1/6}$, which is again consistent with the results of our numerical simulations.
We analyze many-body localization (MBL) to delocalization transition in Sherrington-Kirkpatrick (SK) model of Ising spin glass (SG) in the presence of a transverse field $Gamma$. Based on energy resolved analysis, which is of relevance for a closed quantum system, we show that the quantum SK model has many-body mobility edges separating MBL phase which is non-ergodic and non-thermal from the delocalized phase which is ergodic and thermal. The range of the delocalized regime increases with increase in the strength of $Gamma$ and eventually for $Gamma$ larger than $Gamma_{CP}$ the entire many-body spectrum is delocalized. We show that the Renyi entropy is almost independent of the system size in the MBL phase, hinting towards an area law in this infinite range model while the delocalized phase shows volume law scaling of Renyi entropy. We further obtain spin glass transition curve in energy density $epsilon$-$Gamma$ plane from the collapse of eigenstate spin susceptibility. We demonstrate that in most of the parameter regime SG transition occurs close to the MBL transition indicating that the SG phase is non-ergodic and non-thermal while the paramagnetic phase is delocalized and thermal.
Some recent results concerning the Sherrington-Kirkpatrick model are reported. For $T$ near the critical temperature $T_c$, the replica free energy of the Sherrington-Kirkpatrick model is taken as the starting point of an expansion in powers of $delta Q_{ab} = (Q_{ab} - Q_{ab}^{rm RS})$ about the Replica Symmetric solution $Q_{ab}^{rm RS}$. The expansion is kept up to 4-th order in $delta{bm Q}$ where a Parisi solution $Q_{ab} = Q(x)$ emerges, but only if one remains close enough to $T_c$. For $T$ near zero we show how to separate contributions from $xll Tll 1$ where the Hessian maintains the standard structure of Parisi Replica Symmetry Breaking with bands of eigenvalues bounded below by zero modes. For $Tll x leq 1$ the bands collapse and only two eigenvalues, a null one and a positive one, are found. In this region the solution stands in what can be called a {sl droplet-like} regime.
We develop a simple method to study the high temperature, or high external field, behavior of the Sherrington-Kirkpatrick mean field spin glass model. The basic idea is to couple two different replicas with a quadratic term, trying to push out the two replica overlap from its replica symmetric value. In the case of zero external field, our results reproduce the well known validity of the annealed approximation, up to the known critical value for the temperature. In the case of nontrivial external field, we prove the validity of the Sherrington-Kirkpatrick replica symmetric solution up to a line, which falls short of the Almeida-Thouless line, associated to the onset of the spontaneous replica symmetry breaking, in the Parisi Ansatz. The main difference with the method, recently developed by Michel Talagrand, is that we employ a quadratic coupling, and not a linear one. The resulting flow equations, with respect to the parameters of the model, turn out to be much simpler, and more tractable. By applying the cavity method, we show also how to determine free energy and overlap fluctuations, in the region where replica symmetry has been shown to hold.
An expansion for the free energy functional of the Sherrington-Kirkpatrick (SK) model, around the Replica Symmetric SK solution $Q^{({rm RS})}_{ab} = delta_{ab} + q(1-delta_{ab})$ is investigated. In particular, when the expansion is truncated to fourth order in. $Q_{ab} - Q^{({rm RS})}_{ab}$. The Full Replica Symmetry Broken (FRSB) solution is explicitly found but it turns out to exist only in the range of temperature $0.549...leq Tleq T_c=1$, not including T=0. On the other hand an expansion around the paramagnetic solution $Q^{({rm PM})}_{ab} = delta_{ab}$ up to fourth order yields a FRSB solution that exists in a limited temperature range $0.915...leq T leq T_c=1$.
To test the stability of the Parisi solution near T=0, we study the spectrum of the Hessian of the Sherrington-Kirkpatrick model near T=0, whose eigenvalues are the masses of the bare propagators in the expansion around the mean-field solution. In the limit $Tll 1$ two regions can be identified. In the first region, for $x$ close to 0, where $x$ is the Parisi replica symmetry breaking scheme parameter, the spectrum of the Hessian is not trivial and maintains the structure of the full replica symmetry breaking state found at higher temperatures. In the second region $Tll x leq 1$ as $Tto 0$, the components of the Hessian become insensitive to changes of the overlaps and the bands typical of the full replica symmetry breaking state collapse. In this region only two eigenvalues are found: a null one and a positive one, ensuring stability for $Tll 1$. In the limit $Tto 0$ the width of the first region shrinks to zero and only the positive and null eigenvalues survive. As byproduct we enlighten the close analogy between the static Parisi replica symmetry breaking scheme and the multiple time-scales approach of dynamics, and compute the static susceptibility showing that it equals the static limit of the dynamic susceptibility computed via the modified fluctuation dissipation theorem.