1. Introduction
It has been tacitly or explicitly assumed that Edwards–Anderson Ising Spin Glasses (ISGs), where the quenched interactions are random, follow the same basic scaling and Universality rules as in the canonical Ising models, whose properties are understood in great detail through Renormalization Group Theory (RGT). Here we will present numerical evidence for hyperscaling violation in ISGs. A textbook definition of hyperscaling is: “Identities obtained from the generalised homogeneity assumption involving the space dimension
D are known as hyperscaling relations.” [
1]. The hyperscaling relations valid in canonical Ising models below the upper critical dimension are:
, and
, where
is the “gap” exponent associated with the critical behavior of the higher field derivatives of the free energy [
2,
3]. The two hyperscaling relations are linked through the Essam–Fisher relation
.
Hyperscaling “collapses” in Ising models in dimensions above the upper critical dimension
, where the critical exponents take their mean-field values:
,
,
, and
. Hyperscaling was predicted by Schwartz to break down also in quenched systems with random interactions [
4]. The breakdown of hyperscaling in the 3D Random Field Ising model (RFIM) has been extensively studied [
5,
6,
7,
8]. The first hyperscaling relation in this model is re-written
where
is the “violation of hyperscaling exponent” [
5] with
in the 3D RFIM [
9,
10]. Logically, the second hyperscaling relation should simultaneously become
.
Though not conventionally written this way, in the standard Ising models above equivalent modified hyperscaling relations and can be seen by inspection to be consistent with the mean field exponents plus a violation exponent .
Ising spin glasses (ISGs) are also systems with quenched randomness in which hyperscaling might be expected to break down, from a generalization of Schwartz’s argument. The exponent in ISGs is always strongly negative and so is very hard to estimate directly; we will be concerned only with the second hyperscaling relation. We are not aware of any tests of hyperscaling in ISGs.
2. Scaling
As usual the Hamiltonian is . The sum is taken over all nearest-neighbour spins and the are the interaction strengths. For Ising models we let but for ISGs the are drawn from a symmetric distribution () normalized to . Also, while the inverse temperature is for the Ising case one should keep in mind that for ISGs it is , which is consistent with the Ising case.
The magnetisation
has the role of order parameter for Ising models but for ISGs we use the spin overlap parameter
where
A and
B indicate two spin systems on the same interaction realizations (samples). The susceptibility is then
for Ising models while for ISGs we use
, where
denotes the sample average.
The Binder cumulant [
11], a normalized form of the kurtosis, is defined as
for Ising models and as
for spin glasses. Finally, the standard finite-
L definition for the second-moment correlation length
through the Fourier transformation of the correlation function is used, see Ref. [
12] Equation (
8),
and
where
for ISGs and just
for Ising models. Here,
is the
kth spin coordinate so we let
be the average over the
D directions. For ISGs we then also take the sample average of
F in Equation (
4). The ISG data were gathered using exchange Monte Carlo [
13] and the Ising data with the Wolff-cluster technique.
The conventional RGT-based approach to analyses of numerical simulation data is to use the reduced temperature
as the thermal scaling variable, together with the principal observables
,
and
. This approach is then tailored to the critical region. However, at high temperatures the scaling variable
t diverges and
tends to zero, so analyses of the entire paramagnetic regime is not possible without introducing diverging correction terms. In addition, in symmetric-interaction ISGs the relevant interaction-strength parameter is
so the ISG thermal scaling variable should depend on the square of the temperature; this basic point was made some thirty years ago [
14,
15] but has been ignored since in most ISG simulation data analyses.
A rational scaling approach which takes in the entire paramagnetic region, thus including both the finite-size scaling regime (FSS,
) and the thermodynamic limit regime (ThL,
), can be based on the Wegner scaling expression for the bulk Ising susceptibility [
16]
where
in Ising models with
the inverse temperature [
17]. (The Wegner expression is often mis-quoted with
t replacing
). The terms inside
are scaling corrections, where
is the leading thermal correction exponent, which is identical for all observables within a universality class.
At infinite temperature
, and for all spin-
models,
tends to 1; hence, for the susceptibility the whole paramagnetic region can generally be covered to good precision when a few mild Wegner correction terms are included. (To obtain infinite precision an infinite number of correction terms are of course needed, just as in standard FSS analyses where perfect precision in principle requires an infinite number of size dependent correction terms [
18]).
For ISG models with symmetric interaction distributions an appropriate thermal scaling variable is
, which then can be used in the Wegner expression [
14,
15,
19,
20]. In the ThL regime
the properties of a finite-size sample, if normalized correctly, are independent of
L and so are the same as those of the infinite-
L model. A standard rule of thumb for the approximate onset of the ThL regime is
and this is easily identified in the simulation data. A huge benefit with this approach is that the ThL numerical data can be readily dovetailed into High Temperature Series Expansion (HTSE) values obtained from exact series terms (in practice limited to a finite number of terms). This is lost when the conventional FSS thermal scaling variable
t is used.
To apply the Wegner formalism to other observables Q than in spin- models, it is convenient to impose the rule that each observable should be normalized in such a way that that the infinite temperature limit is , without the critical limit exponent being modified. For the spin- [reduced] susceptibility no normalization is required as this condition is automatically fulfilled, with a temperature dependent effective exponent both in Ising models, and in ISGs with the appropriate .
The reduced second-moment correlation length is defined as
in Ising models and as
in ISG models [
20]. Note that the critical-limit ThL exponent
is unaltered by this normalization (models with zero-temperature critical points are a special case [
21]). From exact and general HTSE for either model this reduced correlation length tends to 1 at infinite temperature [
2,
19].
The temperature-dependent effective exponent is then defined as in Ising models and in ISG models. A Wegner-like relation again applies with appropriate correction terms. With this definition the temperature-dependent effective usually turns out to remain close to the critical value, except at high temperatures where it may be significantly modified by the corrections.
The temperature-dependent effective exponents and are well-behaved over the entire paramagnetic regime, with exact infinite-temperature hypercubic-lattice limits and for Ising models. For the ISG models they are and , where K is the kurtosis of the interaction distribution.
3. Hyperscaling
The second field derivative of the bulk susceptibility
(fourth derivative of
, also called the non-linear susceptibility) in a hypercubic lattice is directly related to the ThL Binder cumulant for finite
L through
see Ref. [
22] Equation (10.2). It can be noted that in any spin-
Ising system the infinite-temperature limit (i.e., independent spins) for the Binder cumulant is
, where N is the number of spins; as
for a hypercubic lattice, at infinite temperature
. Thus this normalized Binder cumulant also obeys the high-temperature limit rule for normalized observables introduced above.
For Ising models in the ThL (bulk or
) regime, assuming hyperscaling the critical exponent for the second field derivative of the susceptibility
is [
2]
Thus the bulk
or the ThL normalized Binder cumulant
scales with the critical exponent
, together with appropriate Wegner correction terms.
The standard “dimensionless renormalized coupling constant” can be defined as
(other normalizations are also used [
2]). It should be noted that even in the case of the canonical 2D Ising model [
23] the infinite-
L value of
at criticality depends strongly on the order in which the limits are taken: the ThL limit
or the FSS limit
. Also, in order for the infinite-
L Ising
to become regular analytic to
the normalized Ising form
should be used, see Ref. [
2] Equation (42). This modification is strictly equivalent to replacing in Equation (
9)
by
which is the reduced correlation length introduced above.
From the hyperscaling rule, the critical exponents for
and
are both
.
Figure 1 shows a plot of
against
for the 3D simple cubic spin-
Ising model (see Ref. [
24] for details of the simulations). It can be seen that within the statistics the ThL data are consistent with a limiting critical slope
and an intercept
, with corrections coming into play at high temperatures, in full agreement with the hyperscaling rule. It should be noted that in this form of plot neither the critical inverse temperature value
nor the critical exponent value
need to be introduced.
For the 5D Ising model, in an equivalent ThL plot covering the entire paramagnetic temperature regime the breakdown of standard hyperscaling leads to a critical exponent for the normalized Binder cumulant which is not
but 2 [
25], i.e.,
with a hyperscaling violation exponent
as discussed above.
4. Ising Spin Glasses
The normalized Binder cumulant against reduced correlation length to the power
D for ISGs can be displayed in just the same way as for the Ising model data in
Figure 1. We show in
Figure 2 and
Figure 3 data for the standard bimodal and Gaussian interaction ISGs on the 3D simple cubic lattice, and in
Figure 4 for a more exotic model, the Laplacian-interaction ISG on a face-centered cubic lattice.
In each case the ThL data can be seen to have a critical limit constant slope with corrections coming into play at high temperatures. The value of the limit slope is in each case distinctly stronger than the standard hyperscaling value , demonstrating that and have different critical exponents; there is a violation of hyperscaling. With the same formalism as above, the value of the slope in each Figure can be taken to be equal to with , and respectively for the three 3D ISG models.
Data for the bimodal and Gaussian 4D models,
Figure 5 and
Figure 6, also show hyperscaling violations with violation exponents
and
respectively, rather weaker than in 3D. For the bimodal and Gaussian 5D models,
Figure 7 and
Figure 8, the limiting slopes are close to
and any violation of hyperscaling is too weak to be observed.
Finally, in the bimodal 7D ISG (above the ISG critical dimension
), the equivalent plot,
Figure 9, is much more noisy than for the lower dimensions simply because by this dimension the number of spins in each individual sample becomes very large leading to practical limitations, particularly for the Binder cumulant. Nevertheless the slope of the plot can be seen to be lower than
so, as for the 5D Ising model, the violation exponent
is positive.