Next Article in Journal
High Order Split Operators for the Time-Dependent Wavepacket method of Triatomic Reactive Scattering in Hyperspherical Coordinates
Next Article in Special Issue
Spin Glasses in a Field Show a Phase Transition Varying the Distance among Real Replicas (And How to Exploit It to Find the Critical Line in a Field)
Previous Article in Journal
A Complexity-Entropy Based Approach for the Detection of Fish Choruses
Previous Article in Special Issue
Kinetic Models of Discrete Opinion Dynamics on Directed Barabási–Albert Networks
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Hyperscaling Violation in Ising Spin Glasses

1
Laboratoire Charles Coulomb (L2C), UMR 5221 CNRS-Université de Montpellier, 34095 Montpellier, France
2
Department of Mathematics and Mathematical Statistics, Umeå University, SE-901 87 Umeå, Sweden
*
Author to whom correspondence should be addressed.
Entropy 2019, 21(10), 978; https://doi.org/10.3390/e21100978
Submission received: 12 September 2019 / Revised: 4 October 2019 / Accepted: 7 October 2019 / Published: 8 October 2019

Abstract

:
In addition to the standard scaling rules relating critical exponents at second order transitions, hyperscaling rules involve the dimension of the model. It is well known that in canonical Ising models hyperscaling rules are modified above the upper critical dimension. It was shown by M. Schwartz in 1991 that hyperscaling can also break down in Ising systems with quenched random interactions; Random Field Ising models, which are in this class, have been intensively studied. Here, numerical Ising Spin Glass data relating the scaling of the normalized Binder cumulant to that of the reduced correlation length are presented for dimensions 3, 4, 5, and 7. Hyperscaling is clearly violated in dimensions 3 and 4, as well as above the upper critical dimension D = 6 . Estimates are obtained for the “violation of hyperscaling exponent” values in the various models.
PACS:
75.50.Lk; 75.40.Mg; 05.50.+q

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: 2 α = D ν , and 2 Δ = D ν + γ , 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 α + D ν 2 Δ = 2 .
Hyperscaling “collapses” in Ising models in dimensions above the upper critical dimension D = 4 , where the critical exponents take their mean-field values: γ = 1 , ν = 1 / 2 , α = 0 , and Δ = 3 / 2 . 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 2 α = ( D θ ) ν where θ is the “violation of hyperscaling exponent” [5] with θ 1 . 47 in the 3D RFIM [9,10]. Logically, the second hyperscaling relation should simultaneously become 2 Δ = ( D θ ) ν + γ .
Though not conventionally written this way, in the standard Ising models above D = 4 equivalent modified hyperscaling relations 2 α = ( D θ ) ν = 2 and 2 Δ = ( D θ ) ν + γ = 3 can be seen by inspection to be consistent with the mean field exponents plus a violation exponent θ = D 4 .
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 H = i j J i j S i S j . The sum is taken over all nearest-neighbour spins and the J i j are the interaction strengths. For Ising models we let J i j = 1 but for ISGs the J i j are drawn from a symmetric distribution ( J i j = 0 ) normalized to J i j 2 = 1 . Also, while the inverse temperature is β = 1 / T for the Ising case one should keep in mind that for ISGs it is β = ( J i j 2 / T 2 ) 1 / 2 , which is consistent with the Ising case.
The magnetisation m = L d i S i has the role of order parameter for Ising models but for ISGs we use the spin overlap parameter
q = 1 L d i S i A S i B
where A and B indicate two spin systems on the same interaction realizations (samples). The susceptibility is then χ = L d m 2 for Ising models while for ISGs we use χ = L d [ q 2 ] , where [ ] denotes the sample average.
The Binder cumulant [11], a normalized form of the kurtosis, is defined as
g = 1 2 3 m 4 m 2 2
for Ising models and as
g = 1 2 3 [ q 4 ] [ q 2 ] 2
for spin glasses. Finally, the standard finite-L definition for the second-moment correlation length ξ ( β , L ) through the Fourier transformation of the correlation function is used, see Ref. [12] Equation (8),
ξ = χ / F 1 2 sin ( π / L )
and
F k = 1 L d x q x exp ( 2 π i x k / L ) 2
where q x = S x A S x B for ISGs and just S x for Ising models. Here, x k is the kth spin coordinate so we let F = ( F 1 + + F D ) / D 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 t = ( T T c ) / T c as the thermal scaling variable, together with the principal observables χ ( t , L ) , ξ ( t , L ) and g ( t , L ) . This approach is then tailored to the critical region. However, at high temperatures the scaling variable t diverges and ξ ( t , L ) 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 J i j 2 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, L ξ ( β , ) ) and the thermodynamic limit regime (ThL, L ξ ( β , ) ), can be based on the Wegner scaling expression for the bulk Ising susceptibility [16]
χ ( τ ) = C χ τ γ 1 + a χ τ ν ω + b χ τ +
where τ = 1 β / β c 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 τ = 1 , and for all spin- 1 / 2 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 τ = 1 ( β / β c ) 2 , which then can be used in the Wegner expression [14,15,19,20]. In the ThL regime L ξ ( β ) 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 L 7 ξ ( β , L ) 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- 1 / 2 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 Q n ( τ = 1 ) 1 , without the critical limit exponent being modified. For the spin- 1 / 2 [reduced] susceptibility no normalization is required as this condition is automatically fulfilled, with a temperature dependent effective exponent γ ( τ ) = ln [ χ ( τ , L ) ] / ln τ both in Ising models, and in ISGs with the appropriate τ .
The reduced second-moment correlation length is defined as ξ ( τ , L ) / β 1 / 2 in Ising models and as ξ ( τ , L ) / β 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 ν ( τ ) = ln [ ξ ( τ , L ) / β 1 / 2 ] / ln τ in Ising models and ν ( τ ) = ln [ ξ ( τ , L ) / β ] / ln τ 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 γ ( 1 ) = 2 D β c and ν ( 1 ) = D β c for Ising models. For the ISG models they are γ ( 1 ) = 2 D β c 2 and ν ( 1 ) = ( D K / 3 ) β c 2 , where K is the kurtosis of the interaction distribution.

3. Hyperscaling

The second field derivative of the bulk susceptibility χ 4 ( β ) (fourth derivative of log Z , also called the non-linear susceptibility) in a hypercubic lattice is directly related to the ThL Binder cumulant for finite L through
2 g ( β , L ) = χ 4 L D χ 2 = 3 m 2 2 m 4 m 2 2
see Ref. [22] Equation (10.2). It can be noted that in any spin- 1 / 2 Ising system the infinite-temperature limit (i.e., independent spins) for the Binder cumulant is g ( 0 , N ) 1 / N , where N is the number of spins; as N = L D for a hypercubic lattice, at infinite temperature L D g ( τ , L ) 1 . 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 L ξ ( τ ) ) regime, assuming hyperscaling the critical exponent for the second field derivative of the susceptibility χ 4 ( β ) is [2]
γ 4 = γ + 2 Δ = D ν + 2 γ
Thus the bulk χ 4 ( β ) / ( 2 χ ( β ) 2 ) or the ThL normalized Binder cumulant L D g ( β , L ) scales with the critical exponent ( ν D + 2 γ ) 2 γ = ν D , together with appropriate Wegner correction terms.
The standard “dimensionless renormalized coupling constant” can be defined as
G 4 ( β ) = χ 4 ( β ) ξ ( β ) D χ ( β ) 2 = L D g ( β , L ) 2 ξ ( β , L ) D
(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 G 4 ( β c ) at criticality depends strongly on the order in which the limits are taken: the ThL limit ( L , β β c ) or the FSS limit ( β β c , L ) . Also, in order for the infinite-L Ising G 4 to become regular analytic to β = 0 the normalized Ising form ( β c / β ) D / 2 G 4 ( β ) should be used, see Ref. [2] Equation (42). This modification is strictly equivalent to replacing in Equation (9) ξ ( β , L ) by ξ ( β , L ) / β 1 / 2 which is the reduced correlation length introduced above.
From the hyperscaling rule, the critical exponents for L D g ( τ , L ) and ( ξ ( τ , L ) / β 1 / 2 ) D are both D ν . Figure 1 shows a plot of L 3 g ( τ , L ) against ( ξ ( τ , L ) / β 1 / 2 ) 3 for the 3D simple cubic spin- 1 / 2 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 1 and an intercept G 4 β c 3 / 2 / 2 1 . 23 , 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 β c 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 D ν = 5 / 2 but 2 [25], i.e., ( D θ ) ν with a hyperscaling violation exponent θ = 1 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 1 , demonstrating that ( ξ ( β , L ) / β ) 3 ) and L 3 g ( β , L ) 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 ( D θ ) / D = θ / D 1 with θ 0 . 80 , 0 . 70 and 0 . 80 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 θ 0 . 48 and 0 . 52 respectively, rather weaker than in 3D. For the bimodal and Gaussian 5D models, Figure 7 and Figure 8, the limiting slopes are close to 1 and any violation of hyperscaling is too weak to be observed.
Finally, in the bimodal 7D ISG (above the ISG critical dimension D = 6 ), 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 1 so, as for the 5D Ising model, the violation exponent θ 1 . 75 is positive.

5. Conclusions

For the canonical simple cubic Ising model, from data presented in the form of the normalized Binder cumulant L 3 g ( τ , L ) against the reduced correlation length to power D, [ ξ / β 1 / 2 ] D , the observed critical-limit scaling is fully consistent with the expected standard hyperscaling log-log slope of 1 (plus mild corrections at high temperatures). The equivalent data plots for ISG models in various dimensions show violations of hyperscaling with violation exponents which evolve regularly with dimension from strongly negative for D = 3 to strongly positive for D = 7 , passing though zero near the upper critical dimension.

Author Contributions

Conceptualization, I.A.C.; Investigation, I.A.C. and P.H.L.; Methodology, I.A.C.; Resources, P.H.L.; Software, P.H.L.; Writing—original draft, I.A.C.; Writing—review & editing, P.H.L.

Funding

This research received no external funding.

Acknowledgments

We would like to thank A. Aharony, P. Butera and R. Kenna for helpful comments, and H. Katzgraber and K. Hukushima for access to their raw numerical data. The computations were performed on resources provided by the Swedish National Infrastructure for Computing (SNIC) at the Chalmers Centre for Computational Science and Engineering (C3SE).

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Simons, B. Phase Transitions and Collective Phenomena; Cambridge University Press: Cambridge, UK, 1997. [Google Scholar]
  2. Butera, P.; Comi, M. Critical universality and hyperscaling revisited for Ising models of general spin using extended high-temperature series. Phys. Rev. B 2002, 65, 144431. [Google Scholar] [CrossRef] [Green Version]
  3. Pelissetto, A.; Vicari, E. Critical phenomena and renormalization-group theory. Phys. Rep. 2002, 368, 549–727. [Google Scholar] [CrossRef] [Green Version]
  4. Schwartz, M. Breakdown of Hyperscaling in Random Systems—An Inequality. Europhys. Lett. 1991, 15, 777. [Google Scholar] [CrossRef]
  5. Aharony, A.; Imry, Y.; Ma, S.K. Lowering of Dimensionality in Phase Transitions with Random Fields. Phys. Rev. Lett. 1976, 37, 1364. [Google Scholar] [CrossRef]
  6. Gofman, M.; Adler, J.; Aharony, A.; Harris, A.B.; Schwartz, M. Evidence for two exponent scaling in the random field Ising model. Phys. Rev. Lett. 1993, 71, 1569. [Google Scholar] [CrossRef] [PubMed]
  7. Vink, R.L.C.; Fischer, T.; Binder, K. Finite-size scaling in Ising-like systems with quenched random fields: Evidence of hyperscaling violation. Phys. Rev. E 2010, 82, 051134. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  8. Fytas, N.G.; Martin-Mayor, V. Universality in the Three-Dimensional Random-Field Ising Model. Phys. Rev. Lett. 2013, 110, 227201. [Google Scholar] [CrossRef] [Green Version]
  9. Middleton, A.A.; Fisher, D.S. Three-dimensional random-field Ising magnet: Interfaces, scaling, and the nature of states. Phys. Rev. B 2002, 65, 134411. [Google Scholar] [CrossRef] [Green Version]
  10. Fernandez, L.A.; Martin-Mayor, V.; Yllanes, D. Critical behavior of the dilute antiferromagnet in a magnetic field. Phys. Rev. B 2011, 84, 100408. [Google Scholar] [CrossRef]
  11. Binder, K. Critical Properties from Monte Carlo Coarse Graining and Renormalization. Phys. Rev. Lett. 1981, 47, 693. [Google Scholar] [CrossRef]
  12. Hasenbusch, M. Finite size scaling study of lattice models in the three-dimensional Ising universality class. Phys. Rev. B 2010, 82, 174433. [Google Scholar] [CrossRef]
  13. Hukushima, K.; Nemoto, K. Exchange Monte Carlo method and application to spin glass simulations. J. Phys. Soc. Jpn. 1996, 65, 1604–1608. [Google Scholar] [CrossRef]
  14. Singh, R.R.P.; Chakravarty, S. Critical Behavior of an Ising Spin-Glass. Phys. Rev. Lett. 1986, 57, 245. [Google Scholar] [CrossRef] [PubMed]
  15. Klein, L.; Adler, J.; Aharony, A.; Harris, A.B.; Meir, Y. Series expansions for the Ising spin glass in general dimension. Phys. Rev. B 1991, 43, 11249. [Google Scholar] [CrossRef] [PubMed]
  16. Wegner, F.J. Corrections to scaling laws. Phys. Rev. B 1972, 5, 4529. [Google Scholar] [CrossRef]
  17. Domb, C.; Sykes, M. Effect of Change of Spin on the Critical Properties of the Ising and Heisenberg Models. Phys. Rev. 1962, 128, 168. [Google Scholar] [CrossRef]
  18. Ferrenberg, A.M.; Xu, J.; Landau, D.P. Pushing the limits of Monte Carlo simulations for the three-dimensional Ising model. Phys. Rev. E 2018, 97, 043301. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Daboul, D.; Chang, I.; Aharony, A. Test of universality in the Ising spin glass using high temperature graph expansion. Eur. Phys. J. B 2004, 41, 231–254. [Google Scholar] [CrossRef]
  20. Campbell, I.A.; Hukushima, K.; Takayama, H. Extended scaling scheme for critically divergent quantities in ferromagnets and spin glasses. Phys. Rev. Lett. 2006, 97, 117202. [Google Scholar] [CrossRef]
  21. Lundow, P.H.; Campbell, I.A. Bimodal and Gaussian Ising spin glasses in dimension two. Phys. Rev. E 2016, 93, 022119. [Google Scholar] [CrossRef] [Green Version]
  22. Privman, V.; Hohenberg, P.C.; Aharony, A. Universal critical-point amplitude relations. In Phase Transitions and Critical Phenomena; Domb, C., Lebowitz, J.L., Eds.; Academic Press: London, UK, 1991; Volume 14, pp. 1–134. [Google Scholar]
  23. Salas, J.; Sokal, A.D. Universal amplitude ratios in the critical two-dimensional Ising model on a torus. J. Stat. Phys. 2000, 98, 551–588. [Google Scholar] [CrossRef]
  24. Campbell, I.A.; Lundow, P.H. Extended scaling analysis of the S=12 Ising ferromagnet on the simple cubic lattice. Phys. Rev. B 2011, 83, 014411. [Google Scholar] [CrossRef]
  25. Lundow, P.H.; Campbell, I.A. Hyperscaling breakdown and Ising spin glasses: The Binder cumulant. Phys. A 2018, 492, 1838–1852. [Google Scholar] [CrossRef] [Green Version]
Figure 1. (Color on line) Dimension 3 simple cubic Ising model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β 1 / 2 ) 3 . Sample sizes: L = 32 , 24, 16, 8, 6 (top to bottom). Blue straight line: slope 1 . 00 . In this and all following figures each line for fixed L begins to bend over towards horizontal when it leaves the ThL regime L > > ξ ( β ) .
Figure 1. (Color on line) Dimension 3 simple cubic Ising model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β 1 / 2 ) 3 . Sample sizes: L = 32 , 24, 16, 8, 6 (top to bottom). Blue straight line: slope 1 . 00 . In this and all following figures each line for fixed L begins to bend over towards horizontal when it leaves the ThL regime L > > ξ ( β ) .
Entropy 21 00978 g001
Figure 2. (Color on line) Dimension 3 simple cubic bimodal ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 32 , 24, 16, 12, 8, 6 (top to bottom). Green straight line: slope 1 . 27 .
Figure 2. (Color on line) Dimension 3 simple cubic bimodal ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 32 , 24, 16, 12, 8, 6 (top to bottom). Green straight line: slope 1 . 27 .
Entropy 21 00978 g002
Figure 3. (Color on line) Dimension 3 simple cubic Gaussian ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 32 , 24, 16, 12, 8, 6 (top to bottom). Green straight line: slope 1 . 23 .
Figure 3. (Color on line) Dimension 3 simple cubic Gaussian ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 32 , 24, 16, 12, 8, 6 (top to bottom). Green straight line: slope 1 . 23 .
Entropy 21 00978 g003
Figure 4. (Color on line) Dimension 3 face-centered cubic Laplacian ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 28 , 24, 20, 18, 16, 14, 12, 10, 8, 6 (top to bottom). Green straight line: slope 1 . 27 .
Figure 4. (Color on line) Dimension 3 face-centered cubic Laplacian ISG model. Normalized Binder cumulant L 3 g ( β , L ) against reduced correlation length to the power 3, 1 / ( ξ ( β , L ) / β ) 3 . Sample sizes: L = 28 , 24, 20, 18, 16, 14, 12, 10, 8, 6 (top to bottom). Green straight line: slope 1 . 27 .
Entropy 21 00978 g004
Figure 5. (Color on line) Dimension 4 hypercubic bimodal ISG model. Normalized Binder cumulant L 4 g ( β , L ) against reduced correlation length to the power 4, 1 / ( ξ ( β , L ) / β ) 4 . Sample sizes: L = 14 , 12, 10, 8, 6, 4 (top to bottom). Green straight line: slope 1 . 12 .
Figure 5. (Color on line) Dimension 4 hypercubic bimodal ISG model. Normalized Binder cumulant L 4 g ( β , L ) against reduced correlation length to the power 4, 1 / ( ξ ( β , L ) / β ) 4 . Sample sizes: L = 14 , 12, 10, 8, 6, 4 (top to bottom). Green straight line: slope 1 . 12 .
Entropy 21 00978 g005
Figure 6. (Color on line) Dimension 4 hypercubic Gaussian ISG model. Normalized Binder cumulant L 4 g ( β , L ) against reduced correlation length to the power 4, 1 / ( ξ ( β , L ) / β ) 4 . Sample sizes: L = 14 , 12, 10, 8, 6, 4 (top to bottom). Green straight line: slope 1 . 13 .
Figure 6. (Color on line) Dimension 4 hypercubic Gaussian ISG model. Normalized Binder cumulant L 4 g ( β , L ) against reduced correlation length to the power 4, 1 / ( ξ ( β , L ) / β ) 4 . Sample sizes: L = 14 , 12, 10, 8, 6, 4 (top to bottom). Green straight line: slope 1 . 13 .
Entropy 21 00978 g006
Figure 7. (Color on line) Dimension 5 hypercubic bimodal ISG model. Normalized Binder cumulant L 5 g ( β , L ) against reduced correlation length to the power 5, 1 / ( ξ ( β , L ) / β ) 5 . Sample sizes: L = 10 , 8, 6, 4 (top to bottom). Green straight line: slope 1 . 00 .
Figure 7. (Color on line) Dimension 5 hypercubic bimodal ISG model. Normalized Binder cumulant L 5 g ( β , L ) against reduced correlation length to the power 5, 1 / ( ξ ( β , L ) / β ) 5 . Sample sizes: L = 10 , 8, 6, 4 (top to bottom). Green straight line: slope 1 . 00 .
Entropy 21 00978 g007
Figure 8. (Color on line) Dimension 5 hypercubic Gaussian ISG model. Normalized Binder cumulant L 5 g ( β , L ) against reduced correlation length to the power 5, 1 / ( ξ ( β , L ) / β ) 5 . Sample sizes: L = 10 , 8, 6, 4 (top to bottom). Green straight line: slope 1 . 00 .
Figure 8. (Color on line) Dimension 5 hypercubic Gaussian ISG model. Normalized Binder cumulant L 5 g ( β , L ) against reduced correlation length to the power 5, 1 / ( ξ ( β , L ) / β ) 5 . Sample sizes: L = 10 , 8, 6, 4 (top to bottom). Green straight line: slope 1 . 00 .
Entropy 21 00978 g008
Figure 9. (Color on line) Dimension 7 hypercubic bimodal ISG model. Normalized Binder cumulant L 7 g ( β , L ) against reduced correlation length to the power 7, 1 / ( ξ ( β , L ) / β ) 7 . Sample sizes: L = 7 , 6, 5, 4, 3 (top to bottom). Green straight line: slope 0 . 75 .
Figure 9. (Color on line) Dimension 7 hypercubic bimodal ISG model. Normalized Binder cumulant L 7 g ( β , L ) against reduced correlation length to the power 7, 1 / ( ξ ( β , L ) / β ) 7 . Sample sizes: L = 7 , 6, 5, 4, 3 (top to bottom). Green straight line: slope 0 . 75 .
Entropy 21 00978 g009

Share and Cite

MDPI and ACS Style

Campbell, I.A.; Lundow, P.H. Hyperscaling Violation in Ising Spin Glasses. Entropy 2019, 21, 978. https://doi.org/10.3390/e21100978

AMA Style

Campbell IA, Lundow PH. Hyperscaling Violation in Ising Spin Glasses. Entropy. 2019; 21(10):978. https://doi.org/10.3390/e21100978

Chicago/Turabian Style

Campbell, Ian A., and Per H. Lundow. 2019. "Hyperscaling Violation in Ising Spin Glasses" Entropy 21, no. 10: 978. https://doi.org/10.3390/e21100978

APA Style

Campbell, I. A., & Lundow, P. H. (2019). Hyperscaling Violation in Ising Spin Glasses. Entropy, 21(10), 978. https://doi.org/10.3390/e21100978

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop