Nikolaos G. Fytas

Strong violation of critical phenomena universality : Wang-Landau study of the two-dimensional Blume-Capel model under bond randomness

We study the pure and random-bond versions of the square lattice ferromagnetic Blume-Capel model, in both the first-order and second-order phase transition regimes of the pure model. Phase transition temperatures, thermal and magnetic critical exponents are determined for lattice sizes in the range L=20-100 via a sophisticated two-stage numerical strategy of entropic sampling in dominant energy subspaces, using mainly the Wang-Landau algorithm. The second-order phase transition, emerging under random bonds from the second-order regime of the pure model, has the same values of critical exponents as the two-dimensional Ising universality class, with the effect of the bond disorder on the specific heat being well described by double-logarithmic corrections, our findings thus supporting the marginal irrelevance of quenched bond randomness. On the other hand, the second-order transition, emerging under bond randomness from the first-order regime of the pure model, has a distinctive universality class with nu=1.30(6) and beta/nu = 0.128(5) . These results amount to a strong violation of universality principle of critical phenomena, since these two second-order transitions, with different sets of critical exponents, are between the same ferromagnetic and paramagnetic phases. Furthermore, the latter of these two sets of results supports an extensive but weak universality, since it has the same magnetic critical exponent (but a different thermal critical exponent) as a wide variety of two-dimensional systems with and without quenched disorder. In the conversion by bond randomness of the first-order transition of the pure system to second order, we detect, by introducing and evaluating connectivity spin densities, a microsegregation that also explains the increase we find in the phase transition temperature under bond randomness.

Multicritical Points and Crossover Mediating the Strong Violation of Universality: Wang-Landau Determinations in the Random-Bond d = 2 Blume-Capel model

The effects of bond randomness on the phase diagram and critical behavior of the square lattice ferromagnetic Blume-Capel model are discussed. The system is studied in both the pure and disordered versions by the same efficient two-stage Wang-Landau method for many values of the crystal field, restricted here in the second-order phase-transition regime of the pure model. For the random-bond version several disorder strengths are considered. We present phase diagram points of both pure and random versions and for a particular disorder strength we locate the emergence of the enhancement of ferromagnetic order observed in an earlier study in the ex-first-order regime. The critical properties of the pure model are contrasted and compared to those of the random model. Accepting, for the weak random version, the assumption of the double-logarithmic scenario for the specific heat we attempt to estimate the range of universality between the pure and random-bond models. The behavior of the strong disorder regime is also discussed and a rather complex and yet not fully understood behavior is observed. It is pointed out that this complexity is related to the ground-state structure of the random-bond version.

Lack of self-averaging of the specific heat in the three-dimensional random-field ising model

We apply the recently developed critical minimum-energy subspace scheme for the investigation of the random-field Ising model. We point out that this method is well suited for the study of this model. The density of states is obtained via the Wang-Landau and broad histogram methods in a unified implementation by employing the N-fold version of the Wang-Landau scheme. The random fields are obtained from a bimodal distribution (hi = +/-2), and the scaling of the specific heat maxima is studied on cubic lattices with sizes ranging from L=4 to L=32. Observing the finite-size scaling behavior of the maxima of the specific heats we examine the question of saturation of the specific heat. The lack of self-averaging of this quantity is fully illustrated, and it is shown that this property may be related to the question mentioned above.

Estimation of critical behavior from the density of states in classical statistical models

We present a simple and efficient approximation scheme which greatly facilitates the extension of Wang-Landau sampling (or similar techniques) in large systems for the estimation of critical behavior. The method, presented in an algorithmic approach, is based on a very simple idea, familiar in statistical mechanics from the notion of thermodynamic equivalence of ensembles and the central limit theorem. It is illustrated that we can predict with high accuracy the critical part of the energy space and by using this restricted part we can extend our simulations to larger systems and improve the accuracy of critical parameters. It is proposed that the extensions of the finite-size critical part of the energy space, determining the specific heat, satisfy a scaling law involving the thermal critical exponent. The method is applied successfully for the estimation of the scaling behavior of specific heat of both square and simple cubic Ising lattices. The proposed scaling law is verified by estimating the thermal critical exponent from the finite-size behavior of the critical part of the energy space. The density of states of the zero-field Ising model on these lattices is obtained via a multirange Wang-Landau sampling.

Entropic sampling via Wang-Landau random walks in dominant energy subspaces.

Dominant energy subspaces of statistical systems are defined with the help of restrictive conditions on various characteristics of the energy distribution, such as the probability density and the fourth order Binders cumulant. Our analysis generalizes the ideas of the critical minimum energy subspace (CRMES) technique, applied previously to study the specific heats finite-size scaling. Here, we illustrate alternatives that are useful for the analysis of further finite-size anomalies and the behavior of the corresponding dominant subspaces is presented for the two-dimensional (2D) Baxter-Wu and the 2D and 3D Ising models. In order to show that a CRMES technique is adequate for the study of magnetic anomalies, we study and test simple methods which provide the means for an accurate determination of the energy-order-parameter (E,M) histograms via Wang-Landau random walks. The 2D Ising model is used as a test case and it is shown that high-level Wang-Landau sampling schemes yield excellent estimates for all magnetic properties. Our estimates compare very well with those of the traditional Metropolis method. The relevant dominant energy subspaces and dominant magnetization subspaces scale as expected with exponents alpha/nu and gamma/nu, respectively. Using the Metropolis method we examine the time evolution of the corresponding dominant magnetization subspaces and we uncover the reasons behind the inadequacy of the Metropolis method to produce a reliable estimation scheme for the tail regime of the order-parameter distribution.

Universality in the three-dimensional random-field ising model

We solve a long-standing puzzle in statistical mechanics of disordered systems. By performing a high-statistics simulation of the D=3 random-field Ising model at zero temperature for different shapes of the random-field distribution, we show that the model is ruled by a single universality class. We compute the complete set of critical exponents for this class, including the correction-to-scaling exponent, and we show, to high numerical accuracy, that scaling is described by two independent exponents. Discrepancies with previous works are explained in terms of strong scaling corrections.

First-order transition features of the 3D bimodal random-field Ising model

Two numerical strategies based on the Wang-Landau and Lee entropic sampling schemes are implemented to investigate the first-order transition features of the 3D bimodal (

Self-assembly of DNA-functionalized colloids

\pm h

Phase diagram of the 3D bimodal random-field Ising model

) random-field Ising model at the strong disorder regime. We consider simple cubic lattices with linear sizes in the range

Scaling and universality in the phase diagram of the 2D Blume-Capel model

L=4-32