Binder cumulant vs temperature for different lattice sizes . The values of obey the Fibonacci sequence. We estimated the critical temperature by averaging the numerical values of the temperatures where the curves intersect each other. We have, for this model, a phase transition from a paramagnetic phase to a spin-glass phase by decreasing the temperature.
The EA order parameter as a function of temperature for different lattice sizes . The values of obey the Fibonacci sequence. The curves suggest a second order phase transition.
Critical behavior of at as a function of lattice size obtained from Eq.(6) of our original article. Alongside the points we show the error bars on the same scale. The curve slope gives the exponent ratio . The exponent ratio differs from the pure model and this change of the universality class is induced by the quasiperiodic ordering.
The susceptibility as a function of temperature for different lattice sizes . The values of obey the Fibonacci sequence. The susceptibility diverges at in the large lattice size limit suggesting a second order phase transition.
Critical behavior of at as a function of lattice size obtained from Eq.(7) of our original article. Alongside the points we show the error bars on the same scale. The curve slope gives the exponent ratio , differing from the pure Ising 2D case.
Critical behavior of susceptibility maximum temperatures as a function of lattice size obtained from Eq.(10) of our original article. The curve slope gives the exponent , differing from the pure Ising 2D case.
Specific heat as a function of temperature for different lattice sizes . The values of obey the Fibonacci sequence. When increasing the lattice size, we observe a crescent maximum, suggesting a logarithm divergence or a negative exponent divergence at the critical temperature .
Data collapse of EA order parameter and susceptibility . The thermodynamic parameters as functions of lattice sizes collapse for , and next to the critical temperature according to the scale forms given in Eqs. (6) and (7) of our original article differing from the pure Ising 2D case.