Non-perturbative effects in spin glasses

We present a numerical study of an Ising spin glass with hierarchical interactions—the hierarchical Edwards-Anderson model with an external magnetic field (HEA). We study the model with Monte Carlo (MC) simulations in the mean-field (MF) and non-mean-field (NMF) regions corresponding to d ≥ 4 and d < 4 for the d-dimensional ferromagnetic Ising model respectively. We compare the MC results with those of a renormalization-group (RG) study where the critical fixed point is treated as a perturbation of the MF one, along the same lines as in the -expansion for the Ising model. The MC and the RG method agree in the MF region, predicting the existence of a transition and compatible values of the critical exponents. Conversely, the two approaches markedly disagree in the NMF case, where the MC data indicates a transition, while the RG analysis predicts that no perturbative critical fixed point exists. Also, the MC estimate of the critical exponent ν in the NMF region is about twice as large as its classical value, even if the analog of the system dimension is within only ~2% from its upper-critical-dimension value. Taken together, these results indicate that the transition in the NMF region is governed by strong non-perturbative effects.

T he renormalization-group (RG) is a powerful method to tackle the complexity of a variety of thermodynamical systems close to their critical point: with the RG, an extensive number of microscopic variables is reduced to a few coarsened degrees of freedom whose behavior can be completely characterized in the thermodynamic limit. The RG theory for homogeneous systems-such as the d-dimensional Ising model on a hypercubic lattice-is based on a perturbative expansion around mean-field (MF) theory. Corrections to the MF fixed point (FP) are described in terms of a power series in E~4{d 1 , and the predictions of this expansion in three dimensions have been shown to be in excellent agreement with both numerical simulations 2 and experimental measurements from a wide variety of systems-binary fluids, superfluid helium, and ferromagnets 3 . For inhomogeneous, disordered systems characterized by a complex energy landscape such as Ising spin glasses 4 , the formulation of a RG theory is more involved: perturbative RG studies based on the E-expansion 5,6 are limited by the intricate diagrammatic structure of the underlying field theory, and the markedly divergent behavior of the E-expansion prevents any physical predictions below the upper critical dimension 7 .
In the last few decades, perturbative RG approaches have been extensively applied to a disordered system of particular relevance in condensed-matter physics: an Ising spin-glass in the presence of an external magnetic field 8 . Indeed, the existence of a spin-glass transition in a field has been raising particular interest 9-12 because evidence in favor or against this transition would provide new insights into the structure of the low-temperature equilibrium states of the system 13 . The spin-glass transition in a field has been extensively studied with RG methods by examining a tentative set of critical FPs [14][15][16] , although the existence of a such a transition below the upper critical dimension is still under debate 16,17 . In this regard, a recent RG study 18 has shown for the first time that the complete set of perturbative FPs in a field can be determined for a spin glass built on a hierarchical lattice 19 : besides the well-known trivial FP in the MF region-corresponding to d $ 4 for a ferromagnetic Ising model in d dimensions-this analysis demonstrated by means of an E-expansion that there is no perturbative FP in the non-mean-field (NMF) region 19 . The latter result provides a natural starting point for testing non-perturbative effects in spin glasses in a field: given that the RG analysis determined the complete set of perturbative FPs, if another FP in the NMF region were found with a non-perturbative approach, then this FP would necessarily lie outside the perturbative regime.
In this paper, we investigate these non-perturbative effects by studying with Monte Carlo (MC) simulations the hierarchical Edwards-Anderson model in a magnetic field (HEA)-a pairwise spin-glass model with binary couplings and magnetic fields where spin couplings are disposed in a hierarchical way. The HEA and the hierarchical model in the RG study above 18 share the same distance dependence of the coupling strength between spins: this common feature between the two models allows us to compare the predictions of MC simulations with those from the RG approach. Specifically, we study the existence of a phase transition by analyzing the temperature dependence of the correlation length and spin-glass susceptibility, and we characterize the value of the critical exponent n 1 . In the MF region, the MC data and the RG approach both predict that there is a transition in a field, and that this transition is characterized by a classical value of the critical exponent. Conversely, in the NMF region the two approaches exhibit a marked discrepancy: The RG method shows that there is no perturbative FP, while the MC data hints that there is a transition. In addition, the MC value of the critical exponent n differs by nearly a factor of two from its classical RG value even if the analog of the system dimension is within only , 2% from its upper-criticaldimension value. Overall, these findings suggest that the phase transition in the NMF region is associated with a FP that cannot be considered as a perturbation of the MF FP, thus hinting at the presence of strong non-perturbative effects.
The rest of the paper is structured as follows: in the Results we introduce the HEA, and in Section A we discuss how a suitable correlation length based on the hierarchical distance can be extracted from the correlation functions. In Section B we present the MC data for the correlation length and overlap probability distribution, and in Section C we discuss the numerical evidence for non-perturbative effects. Finally, in the Conclusions we discuss the physical interpretation of our results and we lay out some topics of future studies.

Results
The HEA is a system of 2 k Ising spins S i 5 61, with Hamiltonian where {J ij } are independent and identically distributed (IID) random variables defined as follows. Given two sites i, j, we denote by d ij the hierarchical distance between i and j, i.e. the number of hierarchical levels that we need to ascend in the hierarchical tree starting from spins i, j to find a root common to i and j, see Fig. 1a. In Eq. (1), we then chose J ij to be different from zero with probability where nonzero J ij s are equal to 61 with equal probability, the parameter s determines how fast spin-spin interactions decrease with distance, and we chose 1/2 , s , 1 because this is the interval where the hierarchical Edwards-Anderson model is expected to have a finite-temperature phase transition 19 . The magnetic fields {h i } in Eq. (1) are also IID random variables: h i is different from zero with probability p h , and nonzero h i s are equal to 61 with equal probability. The choice (2) for the probability of placing a bond implies that the average interaction strength on hierarchical level l . 1 between spins S i and S j is where E½ denotes the average with respect to all random variables. Equation (3) shows that, for a given s, the power-law dependence on l of the average interaction strength is the same as in a HEA with Gaussian couplings studied recently 18 . The observation above suggests that the HEA introduced here and the HEA with Gaussian couplings have the same critical features: this will be the working hypothesis of the rest of the paper. In particular, this assumption implies that the two models share the existence, or the absence, of a spin-glass transition, and that they both have a MF behavior characterized by vanishing order-parameter fluctuations for 1/2 , s # 2/ 3, and a NMF behavior in the region 2/3 , s , 1 where orderparameter fluctuations arise 7,19,20 .
We will now study the existence of a phase transition related to the divergence of the spin-glass susceptibility in the HEA: to do so, we will introduce the correlation length associated with the long-wavelength modes of the spin-glass susceptibility and we will study it as a function of temperature for different system sizes.
A. Correlation length. Since the probability (2) with which non-zero bonds J ij are drawn in the Hamiltonian (1) is a function of the hierarchical distance d ij , the correlation function 21 depends on i and j only through d ij , where AE ae denotes the Boltzmann average with Hamiltonian (1) and inverse temperature b 5 1/T. To obtain the correlation length, we recall that the hierarchical distance d ij is related 22 to the dyadic norm j j 2 by the relation where and the coefficients {a j } are given by the expression in base two of i: Equation (5) implies that the correlation function C ij depends on i, j only through I i ð Þ{I j ð Þ j j 2 : We now relabel the lattice indexes i?I i ð Þ and we perform the Fourier transform of C ij 7 : this leads to the definition of the momentum-dependent spin-glass susceptibility (a) Structure of spin-spin interactions: each dot represents a spin, and arcs represent interactions between spins below them. The green path goes from the bottom to the top of the hierarchical tree until a common arc between the two circled spins is found: since this requires ascending two hierarchical levels, the hierarchical distance between the two spins circled in green is d ij 5 2. (b) Long-wavelength Fourier modes: real part of the Fourier modes exp 2pipI i ð Þ 2 k À Á for the HEA (black solid lines), where the index i 5 0, …, 2 k 2 1 runs over the lattice sites from left to right, and the Fourier modes take values 1 and 21.
We also plot the real part of the Fourier modes exp(2pip i/2 k ) for a onedimensional spin glass where spin-spin interactions decay according to the Euclidean distance | i 2 j | (red dashed lines). Top: first Fourier mode, corresponding to p 5 0 for both the HEA and for the power-law spin glass. Bottom: second Fourier mode, corresponding to p 5 2 k21 for the HEA and to p 5 1 for the power-law spin glass.
www.nature.com/scientificreports SCIENTIFIC REPORTS | 5 : 8697 | DOI: 10.1038/srep08697 where p 5 0, …, 2 k 2 1 is an integer momentum variable, and in the second line we used Eq. (8) and the identity I I i ð Þ ð Þ~i. In spin-glass models with nearest-neighbor interactions, the correlation length is extracted 23 from the Fourier transform of the correlation function by considering the values of p with the smallest norm jpj. For the HEA, the natural definition of distance between sites on the lattice is not the Euclidean norm jpj, but the dyadic norm jpj 2 7 : to extract the long-wavelength behavior, it is thus natural to consider the values of p with the smallest value of jpj 2 . To this end, we observe that the spin-glass susceptibility (9) is the Fourier transform of a function f which depends on i 2 j only through its dyadic norm, thus x SG (p) depends on p only through the dyadic norm jpj 2 of the momentum 22 . The long-wavelength behavior is then encoded into the limit jpj 2 R 0, in which Eq. (9) takes the modified Ornstein-Zernike form 21 To extract j k from Eq. (10), we consider the momentum values with the smallest dyadic norm: these are p 5 0 and p 5 2 k21 , and their dyadic norms read j0j 2 5 0 and j2 k21 j 2 5 2 12k respectively 22 . The correlation length j k is then obtained by considering Eq. (10) for these values of p: We will now discuss a convenient way of computing the correlation length in terms of products of spin overlaps 24 , and thus extract j k from the MC data. We introduce a set of replicasS a È É of the spin configurationS, and we set Q ab p ð Þ: By using Eqs. (9), (12), we have where in what follows the average AE ae of a function ofS 1 ,S 2 , … denotes the Boltzmann average over all replicas 24 , each replica having an independent Boltzmann measure with Hamiltonian (1) and inverse temperature b. We now introduce the overlaps in the left and right half of the lattice Q R ab : and by using Eqs. (12), (14), we have where in the third line we observed that I i ð Þ is even if 0 # i # 2 k21 2 1 and odd if 2 k21 # i # 2 k 2 1, see Eqs. (6), (7). By using Eqs. (11), (13), (15), (16) we obtain the final expression for the correlation length The derivation above provides a numerically convenient way of extracting the correlation length: to compute j k , we do not need the full set of spin-spin correlations {C ij }, but only the products of the left and right overlaps in Eq. (17).
The first two long-wavelength Fourier modes which determine j k , i.e. exp 2pipI i ð Þ 2 k À Á for p 5 0, 2 k21 , are depicted in Fig. 1b: in particular, the second Fourier mode is {1 ð Þ I i ð Þ , and it is equal to 1 and 21 in the left and right half of the lattice respectively, see Eq. (16). This particular form of the Fourier modes is due to the hierarchical structure of interactions in the model: for example, a onedimensional spin glass where interactions depend on the Euclidean distance rather than on the hierarchical distance 6 possesses a different second Fourier mode, given by exp(2pi i/2 k ).
B. Numerical results for the correlation length and overlap probability distribution. We will now study the existence of a phase transition by analyzing the behavior of the correlation length as a function of temperature: namely, we will assume that if there exists a critical point, then for T close to the critical temperature T c and large enough k the following finite-size-scaling relations hold 9,10 where the correlation-length critical exponent exponent n is equal to 1/(2s 2 1) for 1/2 , s # 2/3 7 and g MF , g NMF are two different scaling functions in the MF and NMF region respectively. We computed the correlation length (17) with MC simulations performed with the parallel-tempering algorithm 25 . For zero magnetic field the Hamiltonian (1) is symmetric under a global spin flip, thus the second and third term in both the numerator and denominator of Eq. (17) vanish and j k can be computed by simulating only two replicas per temperature, while for nonzero field four replicas per temperature are needed. We combined parallel tempering with the asynchronous multispin-coding method, where 64 disorder samples are simulated simultaneously by encoding the values of the couplings J ij into a 64-bit integer 26 . This method allowed for equilibrating an extensive number of samples 8 3 10 3 # S # 16 3 10 3 and system sizes 2 6 , …, 2 13 .
Let us first focus on the zero-field case p h 5 0: in the left panel of Fig. 2 we show j k /2 kn/3 vs. T in the MF region, s 5 0.6, while in the right panel we plot j k /2 k vs. T in the NMF region, s 5 0.68. In the left and right insets we plot j k /2 kn/3 and j k /2 k respectively vs. the logarithm of the number of MC sweeps, in order to provide an equilibration test for the data in the main panels. The correlation-length curves cross at a finite critical temperature, showing that there is a spin-glass transition both in the MF and in the NMF case. This occurrence of a spin-glass transition in zero field is in agreement with a previous RG analysis for a HEA with Gaussian couplings and zero magnetic field, where a critical FP has been shown to exist both in the whole MF region 1/2 , s # 2/3 and in the NMF region s~2=3zE by means of an E-expansion 7 .
To study the nonzero-field case, we need to choose a value for the probability p h of placing a local magnetic field: it can be shown that a suitable value of p h is obtained as a tradeoff between two effects. First, for large p h the critical temperature-if any-is too low to achieve equilibration for sufficiently large system sizes. Second, for small p h a phase transition might appear, but this may be a spurious effect of the transition for p h 5 0 shown in Fig. 2 rather than an actual transition in a finite field. In this regard, a natural way of checking whether the apparent transition is an artifact of zero-field effects is to compute the overlap distribution If there is a phase transition, in the thermodynamic limit p(q) must converge to d(q 2 q 0 ) at the critical point 24 , where q 0 :E Q 12 h i ½ , and q 0 is positive because When the system size is finite, the effect of the zero-field transition arises in the shape of a long exponential tail of p(q) for q , 0 9,27 : thus, a natural way of ensuring that we are in the true nonzero-field regime is to choose p h large enough that the support of p(q) is mostly localized for positive q, i.e. s q q 0 =1, where s 2 q :E Q 12 {q 0 ð Þ 2 Â Ã . We have thus chosen the value of p h according to the tradeoff above. In the MF region, s 5 0.6, we take p h 5 0.36: the left panel in Fig. 3 shows that for this value of p h the curves of j k /2 kn/3 vs. T cross at a finite critical temperature T c . In addition, the right panel of Fig. 3 shows that for this value of p h the zero-field effects are weak: According to the top inset, the overlap distribution at T c has the typical left exponential tail due to zero-field effects 27 , nevertheless the main panel illustrates that the support of p(q) is mostly localized at positive q, and the bottom inset shows that s q q 0 =1 for large k. The same qualitative scenario arises in the NMF region, s 5 0.68, see Fig. 4: The left panel indicates that the curves j k /2 k cross at a finite critical temperature, hinting at the existence of a phase transition in a field, and the right panel illustrates that such a transition is not an artifact due to the effects of the zero-field transition.
C. Hallmark of non-perturbative effects. The correlation-curve crossings in the right panel of Fig. 4 suggest that there is a spinglass transition in the NMF region s . 2/3: interestingly, this result is at variance with a recent perturbative RG analysis for a HEA with Gaussian interactions 18 , which showes with an Eexpansion 1 that for s~2=3zE there is no perturbative FP corresponding to a spin-glass transition in a field. The observation above thus raises the possibility that the phase transition in the right panel of Fig. 4 results from non-perturbative effects, a scenario that we will investigate further in what follows.
The RG study mentioned above shows that if there is a spin-glass transition described by a perturbative FP, then at the critical point the rescaled spin-glass susceptibility x SG 2 k(122s) must converge to a finite value for k R ', where x SG ; x SG (0) 18 . The same RG analysis demonstrates that in the MF region there is a critical FP, hence in this region the rescaled susceptibility x SG 2 k(122s) must converge to a finite value at the critical temperature. This scenario is confirmed by the numerical data in the left panel of Fig. 5, which shows that in the MF region x SG 2 k(122s) converges to a finite FP at a finite temperature T x c , which is within ,17% from the critical temperature T c obtained from the correlation length in Fig. 3. The situation is completely different in the NMF region: there, the RG approach shows that there exists no physical FP, thus suggesting that x SG 2 k(122s) does not have a finite FP at any temperature 18 . Still, the right panel of Fig. 5 illustrates that x SG 2 k(122s) converges to a finite FP at a finite temperature T x c . In addition, T x c coincides within error bars with the critical temperature T c obtained from the correlation length in Fig. 4.
Overall, the discrepancies between the predictions of the perturbative RG analysis and the numerical results in the NMF region raise the possibility that the phase transition resulting from the numerics in the NMF region is associated with a non-perturbative FP. In this regard, a general feature characterizing a critical FP is the exponent n governing the exponential departure of RG trajectories as the temperature is moved away from its critical value 1 . To study the existence of a non-perturbative FP, we estimate n by means of the finite-size-scaling equations (18), (19): we chose n in such a way that for large k the curves j k /2 kn/3 vs. 2 k/3 (T2T c ) and j k /2 k vs. 2 k/v (T2T c ) in the MF and NMF region respectively merge into a single curve. In the MF region, s 5 0.6, we have n 5 5.0 6 0.4, see Fig. 6 left, in Left panel: j k /2 kn/3 vs. temperature T in the mean-field region, s 5 0.6 , 2/3, for system sizes 2 k , with k 5 6, …, 13 and n 5 1/(2s 2 1): the curves cross at a critical temperature T c 5 2.38 6 0.24, which is also marked. Inset: j k /2 kn/3 vs. log 2 t for the lowest temperature shown for each system size in the left panel, where t is the number of Monte Carlo sweeps, and for every t we plot j k /2 kn/3 computed from the last t/2 sweeps. Right panel: j k /2 k vs. T in the non-mean-field region, s 5 0.68 . 2/3, with a critical-temperature value T c 5 1.65 6 0.11. Inset: j k /2 k vs. log 2 t for the lowest temperature shown for each system size in the right panel, where j k /2 k is computed as in the left-panel inset.
www.nature.com/scientificreports SCIENTIFIC REPORTS | 5 : 8697 | DOI: 10.1038/srep08697 agreement with the RG classical value n 5 1/(2s 2 1) 5 5 associated with the MF FP 7 . The situation is again very different in the NMF region: if we use the MF FP prediction n 5 1/(2s 2 1) to estimate n for s 5 0.68, we obtain n , 2.78. This values differs by nearly a factor of two from the MC estimate n 5 5.0 6 1.0 in the right panel of Fig. 6, even if s is within only , 2% from its upper-criticaldimension value s 5 2/3. This discrepancy suggests that the FP giving rise to the phase transition in the numerics cannot be regarded as a slight perturbation of a MF FP with classical exponents. Importantly, this last observation may be at the bottom of the overall disagreement between the numerics and the RG analysis: the advance here is that the RG approach may fail to capture the critical behavior in the numerics because it is based on the hypothesis that the FP governing the transition in the NMF region is a perturbation of the MF one.

Conclusions
We investigated the existence of non-perturbative effects in spin glasses with a Monte Carlo (MC) study for an Ising spin glass in a field where pairwise couplings are disposed in a hierarchical waythe hierarchical Edwards-Anderson model in a magnetic field (HEA). The spin-interaction decay in the HEA is controlled by the parameter s, which is reminiscent of the dimension d for a ferromagnetic Ising model 28 : in particular, the mean-field (MF) and nonmean-field (NMF) regions s # 2/3 and s . 2/3 correspond to d $ 4 and d , 4 respectively 19 . The HEA is suitable for assessing the existence of non-perturbative effects for two reasons: First, the binary nature of couplings and magnetic fields makes the HEA fit for non-perturbative numerical approaches, such as MC simulations 26 . Second, the hierarchical structure of the couplings reproduces the 0.16. Left panel: j k /2 k vs. temperature T and system sizes 2 k , k 5 9, …, 13. Smaller sizes are affected by strong finite-size effects, and they are not shown in order to display clearly the curve crossings for larger sizes. The curves cross at a critical temperature T c 5 0.88 6 0.07, which is also marked. Inset: j k /2 k vs. log 2 t for the lowest temperature shown for each system size in the left panel, where t is the number of Monte Carlo sweeps, and for every t we plot j k /2 k computed from the last t/2 sweeps. Right panel: p(q) vs. q for T 5 T c and same system sizes as in the left panel. Top inset: same plot as in the right panel in logarithmic scale on the vertical axis. Bottom inset: ratio between the standard deviation s q and the mean q 0 of p(q) as a function of k. Left panel: j k /2 kn/3 vs. temperature T and system sizes 2 k , with k 5 8, …, 13. Smaller sizes are affected by strong finite-size effects, and they are not shown in order to display clearly the curve crossings for larger sizes. The curves cross at a critical temperature T c 5 1.46 6 0.07, which is also marked. Inset: j k /2 kn/3 vs. log 2 t for the lowest temperature shown for each system size in the left panel, where t is the number of Monte Carlo sweeps, and for every t we plot j k /2 kn/3 computed from the last t/2 sweeps. Right panel: p(q) vs. q for T 5 T c and same system sizes as in the left panel. Top inset: same plot as in the right panel in logarithmic scale on the vertical axis. Bottom inset: ratio between the standard deviation s q and the mean q 0 of p(q) as a function of k.
www.nature.com/scientificreports same average interaction structure as in a HEA with Gaussian couplings recently studied with perturbative renormalization-group (RG) methods 18 . This RG analysis derived for the first time the complete set of perturbative RG fixed points (FPs): in the MF region, it was shown that there is a FP corresponding to a phase transition, while no such FP exists in the NMF region. In the MF region the MC results are in agreement with the RG predictions: the numerics show the existence of a finite-temperature transition where the correlation length and the spin-glass susceptibility become singular at large sizes, providing an estimate of the exponent n in accordance with its classical value 7 . Conversely, the MC and the RG approach markedly disagree in the NMF region, where the MC data for both the correlation length and spin-glass susceptibility support the existence of a finite-temperature transition. Also, the MC estimate of n in the NMF region differs by nearly a factor of two from its classical value, even if s is within only ,2% from the value s 5 2/3 separating the MF and NMF regions. This picture is at variance with other systems where perturbation theory is well-behaved such as the d-dimensional Ising model, where the MC estimate of n in d 5 3, i.e. within 25% from the upper critical dimension, differs by only ,25% from its classical value 2 . Taken together, these findings hint that the spin-glass transition resulting from the numerics in the NMF region cannot be captured with a perturbative framework, raising the possibility that this transition is related to a FP which cannot be described as a perturbation of the MF one.
These findings raise a variety of possible scenarios concerning the structure of the low-temperature phase. On the one hand, the occurrence of a transition in a field is associated with the replicasymmetry-breaking (RSB) scenario, a MF picture characterized by the existence of exponentially many low-lying energy states 24 . Thus, the evidence provided here in favor of such a transition hints that the MF picture is accurate at least in some part of the NMF region. On Figure 5 | Rescaled spin-glass susceptibility as a function of temperature for nonzero magnetic field. Left: x SG 2 k(122s) vs. T in the mean-field region, s 5 0.6 , 2/3, for p h 5 0.36 and system sizes 2 k , with k 5 6, …, 13. Equilibration tests for x SG have been performed along the same lines as in Figs. 2, 3, 4, and they are not shown. The curves cross at a finite temperature T x c~1 :71+0:06, which is also marked. Right: same plot as in the left panel in the non-mean-field region s 5 0.68 . 2/3 for p h 5 0.16 and k 5 6, …, 13. The crossing temperature of two curves with system sizes 2 k and 2 k11 is an increasing function of k, and the estimated infinite-volume critical temperature is T x c~0 :88+0:06. Figure 6 | Estimate of the critical exponent n with a nonzero magnetic field. Left: j k /2 kn/3 vs. 2 k/3 (T 2 T c ) in the mean-field region, s 5 0.6 , 2/3, for system sizes 2 k , k 5 6, …, 13, where n 5 5.0 and T c is determined by Fig. 3. The gray line envelopes the curves for largest system sizes, and it provides an estimate of the finite-size-scaling function g MF . The value of n has been chosen to achieve the best possible overlap between the curves enclosed in the gray envelope, yielding the estimate n 5 5.0 6 0.4. Right: j k /2 k vs. 2 k/n (T 2 T c ) in the non-mean-field region, s 5 0.68 . 2/3, for system sizes 2 k , k 5 6, …, 13, where n 5 5.0 and T c has been determined by Fig. 4. The the finite-size-scaling function g NMF is also shown, and n has been estimated in the same way as in the left panel, providing n 5 5.0 6 1.0.
www.nature.com/scientificreports the other hand, our findings show that the phase transition in the NMF region cannot be captured by an approach where corrections to the MF picture are considered as a perturbation. This fact raises the possibility that the physical features of this transition are markedly different from the MF ones: this possibility is reminiscent of scenarios alternative to the RSB one, such as the droplet picture 13 or other scenarios intermediate between the RSB and the droplet one 29 . In this regard, the structure of the low-lying energy states could be directly analyzed by probing the low-energy excitations above the ground state 29,30 and by characterizing, for example, their stiffness 31 and surface dimension 32 . A interesting possibility would then be to study how the resulting picture for the low-energy landscape relates to the existence of non-perturbative effects. Taken together, the points above provide an interesting direction that is worth exploring in future studies.