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

Universality in the two-dimensional dilute Baxter-Wu model

We study the question of universality in the two-dimensional spin-1 Baxter-Wu model in the presence of a crystal field Δ. We employ extensive numerical simulations of two types, providing us with complementary results: Wang-Landau sampling at fixed values of Δ and a parallelized variant of the multicanonical approach performed at constant temperature T. A detailed finite-size scaling analysis in the regime of second-order phase transitions in the (Δ,T) phase diagram indicates that the transition belongs to the universality class of the four-state Potts model. Previous controversies with respect to the nature of the transition are discussed and attributed to the presence of strong finite-size effects, especially as one approaches the pentacritical point of the model.

First-principles calculation of the configurational energy density of states for a solid-state ion conductor with a variant of the Wang and Landau algorithm

In this work, a variant of the Wang and Landau algorithm for calculation of the configurational energy density of states is proposed. The algorithm was developed for the purpose of using first-principles simulations, such as density functional theory, to calculate the partition function of disordered sublattices in crystal materials. The expensive calculations of first-principles methods make a parallel algorithm necessary for a practical computation of the configurational energy density of states within a supercell approximation of a solid-state material. The algorithm developed in this work is tested with the two-dimensional (2d) Ising model to bench mark the algorithm and to help provide insight for implementation to a materials science application. Tests with the 2d Ising model revealed that the algorithm has good performance compared to the original Wang and Landau algorithm and the 1/t algorithm, in particular the short iteration performance. A proof of convergence is presented within an adiabatic assumption, and the analysis is able to correctly predict the time dependence of the modification factor to the density of states. The algorithm was then applied to the lithium and lanthanum sublattice of the solid-state lithium ion conductor Li_{0.5}La_{0.5}TiO_{3}. This was done to help understand the disordered nature of the lithium and lanthanum. The results find, overall, that the algorithm performs very well for the 2d Ising model and that the results for Li_{0.5}La_{0.5}TiO_{3} are consistent with experiment while providing additional insight into the lithium and lanthanum ordering in the material. The primary result is that the lithium and lanthanum become more mixed between layers along the c axis for increasing temperature. In part, the simulation of the disordered Li_{0.5}La_{0.5}TiO_{3} system serves as a benchmark for what size systems are currently and in the near future practical to calculate with density functional theory methods.

Statistical Physics Approach to Liquid Crystals: Dynamics of Mobile Potts Model Leading to Smectic Phase, Phase Transition by Wang-Landau Method

We study the nature of the smectic-isotropic phase transition using a mobile 6-state Potts model. Each Potts state represents a molecular orientation. We show that with the choice of an appropriate microscopic Hamiltonian describing the interaction between individual molecules modeled by a mobile 6-state Potts spins, we observe the smectic phase dynamically formed when we cool the molecules from the isotropic phase to low temperatures (T). In order to elucidate the order of the transition and the low-T properties, we use the high-performance Wang-Landau flat energy-histogram technique. We show that the smectic phase goes to the liquid (isotropic) phase by melting/evaporating layer by layer starting from the film surface with increasing T. At a higher T, the whole remaining layers become orientationally disordered. The melting of each layer is characterized by a peak of the specific heat. Such a succession of partial transitions cannot be seen by the Metropolis algorithm. The successive layer meltings/evaporations at low T are found to have a first-order character by examining the energy histogram. These results are in agreement with experiments performed on some smectic liquid crystals.

Monte Carlo study of the triangular Blume-Capel model under bond randomness

The effects of bond randomness on the universality aspects of a two-dimensional (d = 2) Blume-Capel model embedded in the triangular lattice are discussed. The system is studied numerically in both its first- and second-order phase-transition regimes by a comprehensive finite-size scaling analysis for a particularly suitable value of the disorder strength. We find that our data for the second-order phase transition, emerging under random bonds from the second-order regime of the pure model, are compatible with the universality class of the two-dimensional (2D) random Ising model. Furthermore, we find evidence that, the second-order transition emerging under bond randomness from the first-order regime of the pure model, belongs again to the same universality class. Although the first finding reinforces the scenario of strong universality in the 2D Ising model with quenched disorder, the second is in difference from the critical behavior, emerging under randomness, in the cases of the ex-first-order transitions of the Potts model. Finally, our results verify previous renormalization-group calculations on the Blume-Capel model with disorder in the crystal-field coupling.

Complete high-precision entropic sampling

Monte Carlo simulations using entropic sampling to estimate the number of configurations of a given energy are a valuable alternative to traditional methods. We introduce tomographic entropic sampling, a scheme which uses multiple studies, starting from different regions of configuration space, to yield precise estimates of the number of configurations over the full range of energies, without dividing the latter into subsets or windows. Applied to the Ising model on the square lattice, the method yields the critical temperature to an accuracy of about 0.01%, and critical exponents to 1% or better. Predictions for system sizes L=10-160, for the temperature of the specific heat maximum, and of the specific heat at the critical temperature, are in very close agreement with exact results. For the Ising model on the simple cubic lattice the critical temperature is given to within 0.003% of the best available estimate; the exponent ratios β/ν and γ/ν are given to within about 0.04% and 1%, respectively, of the literature values. In both two and three dimensions, results for the antiferromagnetic critical point are fully consistent with those of the ferromagnetic transition. Application to the lattice gas with nearest-neighbor exclusion on the square lattice again yields the critical chemical potential and exponent ratios β/ν and γ/ν to good precision.

Flat energy-histogram simulation of the phase transition in an Ising fully frustrated lattice

We show in this paper the results of the phase transition of the so-called fully frustrated simple cubic lattice with the Ising spin model. We use here the Monte Carlo method with the flat energy-histogram Wang-Landau technique which is very powerful for detecting weak first-order phase transition. We show that the phase transition is clearly of first order, providing an answer to a question raised 25 years ago.

First-order transition in the XY model on a fully frustrated simple cubic lattice

We study the nature of the phase transition in the fully frustrated simple cubic lattice with the XY spin model. This system is the Villain's model generalized in three dimensions. The ground state is very particular with a 12-fold degeneracy. Previous studies have shown unusual critical properties. With the powerful Wang-Landau flat-histogram Monte Carlo method, we carry out in this work intensive simulations with very large lattice sizes. We show that the phase transition is clearly of first order, putting an end to the uncertainty which has lasted for more than 20 years.

Critical behavior of the pure and random-bond two-dimensional triangular Ising ferromagnet

We investigate the effects of quenched bond randomness on the critical properties of the two-dimensional ferromagnetic Ising model embedded in a triangular lattice. The system is studied in both the pure and disordered versions by the same efficient two-stage Wang-Landau method. In the first part of our study, we present the finite-size scaling behavior of the pure model, for which we calculate the critical amplitude of the specific heat's logarithmic expansion. For the disordered system, the numerical data and the relevant detailed finite-size scaling analysis along the lines of the two well-known scenarios-logarithmic corrections versus weak universality--strongly support the field-theoretically predicted scenario of logarithmic corrections. A particular interest is paid to the sample-to-sample fluctuations of the random model and their scaling behavior that are used as a successful alternative approach to criticality.

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.

Crossover from first- to second-order transition in frustrated Ising antiferromagnetic films

In the bulk state, the Ising face-centered-cubic (fcc) antiferromagnet is fully frustrated and is known to have a very strong first-order transition. In this paper, we study the nature of this phase transition in the case of a thin film as a function of the film thickness. Using Monte Carlo simulations, we show that the transition remains first order down to a thickness of four fcc cells (eight atomic layers). It becomes clearly second order at a thickness of two fcc cells, i.e., four atomic layers. It is also interesting to note that the presence of the surface reduces the ground-state degeneracy found in the bulk. For the two-cell thickness, the surface magnetization is larger than the interior one. It undergoes a second-order phase transition at a temperature TC while interior spins become disordered at a lower temperature TD. This loss of order is characterized by a peak of the interior spins susceptibility and a peak of the specific heat which do not depend on the lattice size suggesting that either it is not a real transition or it is a Kosterlitz-Thouless nature. The surface transition, on the other hand, is shown to be of second order with critical exponents deviated from those of pure two-dimensional Ising universality class. We also show results obtained from the Green's function method. A discussion is given.

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.

Phase transition in Heisenberg stacked triangular antiferromagnets: end of a controversy

By using the Wang-Landau flat-histogram Monte Carlo (MC) method for very large lattice sizes never simulated before, we show that the phase transition in the frustrated Heisenberg stacked triangular antiferromagnet is of first order, contrary to results of earlier MC simulations using old-fashioned methods. Our result lends support to the conclusion of a nonperturbative renormalization group performed on an effective Hamiltonian. It puts an end to a 20-year -long controversial issue.

Critical endpoint behavior in an asymmetric Ising model: application of Wang-Landau sampling to calculate the density of states

Using the Wang-Landau sampling method with a two-dimensional random walk we determine the density of states for an asymmetric Ising model with two- and three-body interactions on a triangular lattice, in the presence of an external field. With an accurate density of states we were able to map out the phase diagram accurately and perform quantitative finite-size analyses at, and away from, the critical endpoint. We observe a clear divergence of the curvature of the spectator phase boundary and of the magnetization coexistence diameter derivative at the critical endpoint, and the exponents for both divergences agree well with previous theoretical predictions.

Universal features and tail analysis of the order-parameter distribution of the two-dimensional Ising model: an entropic sampling Monte Carlo study

We present a numerical study of the order-parameter probability density function (PDF) of the square Ising model for lattices with linear sizes L = 80 - 140. A recent efficient entropic sampling scheme, combining the Wang-Landau and broad histogram methods and based on the high levels of the Wang-Landau process in dominant energy subspaces is employed. We find that for large lattices there exists a stable window of the scaled order-parameter in which the full ansatz including the pre-exponential factor for the tail regime of the universal PDF is well obeyed. This window is used to estimate the equation of state exponent and to observe the behavior of the universal constants implicit in the functional form of the universal PDF. The probability densities are used to estimate the universal Privman-Fisher coefficient and to investigate whether one could obtain reliable estimates of the universal constants controlling the asymptotic behavior of the tail regime.

Performances of Wang-Landau algorithms for continuous systems

The relative performances of different implementations of the Wang-Landau method are assessed on two classes of systems with continuous degrees of freedom, namely, two polypeptides and two atomic Lennard-Jones clusters. Parallel tempering Monte Carlo simulations serve as a reference, and we pay particular attention to the variations of the multiplicative factor f during the course of the simulation. For the systems studied, the Wang-Landau method is found to be of comparable accuracy as parallel tempering, but has significant difficulties in reproducing low-temperature transitions exhibited by the Lennard-Jones clusters at low temperature. Using a complementary order parameter and calculating a two-dimensional joint density of states significantly improves the situation, especially for the notoriously difficult LJ(38) system. However, while parallel tempering easily converges for LJ(31), we have not been able to get data of comparable accuracy with Wang-Landau multicanonical sampling.

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.