We study the interface tension of the 4-state Potts model in three dimensions using the Wang- Landau algorithm. The interface tension is given by the ratio of the partition function with a twisted boundary condition in one direction and periodic boundary conditions in all other directions over the partition function with periodic boundary conditions in all directions. With the Wang-Landau algorithm we can explicitly calculate both partition functions and obtain the result for all temperatures. We find solid numerical evidence for perfect wetting. Our algorithm is tested by calculating thermodynamic quantities at the phase transition point.
We study the interface tension of the 4-state Potts model in three dimensions using the Wang- Landau algorithm. The interface tension is given by the ratio of the partition function with a twisted boundary condition in one direction and periodic boundary conditions in all other directions over the partition function with periodic boundary conditions in all directions. With the Wang-Landau algorithm we can explicitly calculate both partition functions and obtain the result for all temperatures. We find solid numerical evidence for perfect wetting. Our algorithm is tested by calculating thermodynamic quantities at the phase transition point.
The concept of the partition function zeros is one of the most efficient methods for investigating the phase transitions and the critical phenomena in various physical systems. Estimating the partition function zeros requires information on the density of states Ω( E ) as a function of the energy E . Currently, the Wang-Landau Monte Carlo algorithm is one of the best methods for calculating Ω( E ). The partition function zeros in the complex temperature plane of the Ising model on an L × L square lattice ( L = 10 ~ 80) with a periodic boundary condition have been estimated by using the Wang-Landau Monte Carlo algorithm. The efficiency of the Wang-Landau Monte Carlo algorithm and the accuracies of the partition function zeros have been evaluated for three different, 5%, 10%, and 20%, flatness criteria for the histogram H ( E ).
Phase transitions and critical phenomena are the most universal phenomena in nature.To understand the phase transitions and critical phenomena of a given system as a continuous function of temperature and to obtain the partition function zeros in the complex temperature plane indicating most effectively phase transitions and critical phenomena, we need to calculate the density of states.Currently, Wang-Landau Monte Carlo algorithm is one of the most efficient Monte Carlo methods to calculate the approximate density of states.Using Wang-Landau Monte Carlo algorithm, the density of states for the Ising model on L x L square lattices (L = 4 ~ 32) with periodic boundary conditions is obtained, and the partition function zeros of the Ising model are evaluated in the complex temperature plane.By examining the behavior of the first partition function zero (partition function zero closest to the positive real axis), phase transitions and critical phenomena can be much more accurately analyzed.The approximate first zeros of the Ising ferromagnet, obtained from Wang-Landau algorithm, are quite close to the exact ones, indicating that it is a reliable method for calculating the density of states and the first partition function zeros.
We study interfaces with periodic boundary conditions in the low-temperature phase of the improved Blume-Capel model on the simple cubic lattice. The interface free energy is defined by the difference of the free energy of a system with antiperiodic boundary conditions in one of the directions and that of a system with periodic boundary conditions in all directions. It is obtained by integration of differences of the corresponding internal energies over the inverse temperature. These differences can be computed efficiently by using a variance reduced estimator that is based on the exchange cluster algorithm. The interface tension is obtained from the interface free energy by using predictions based on effective interface models. By using our numerical results for the interface tension σ and the correlation length ξ obtained in previous work, we determine the universal amplitude ratios R_{2nd,+}=σ_{0}f_{2nd,+}^{2}=0.3863(6), R_{2nd,-}=σ_{0}f_{2nd,-}^{2}=0.1028(1), and R_{exp,-}=σ_{0}f_{exp,-}^{2}=0.1077(3). Our results are consistent with those obtained previously for the three-dimensional Ising model, confirming the universality hypothesis.
We consider the exact correlation length calculations for the two-dimensional Potts model at the transition point βt by Klümper, Schadschneider and Zittartz, and by Buffenoir and Wallon. We argue that the correlation length calculated by the latter authors is the correlation length in the disordered phase and then combine their result with duality and the assumption of complete wetting to give an explicit formula for the order-disorder interface tension σod of this model. The result is used to clarify a controversy stemming from different numerical simulations of σod.
In honor of Onsager’s ninetieth birthday, we like to review some exact results obtained so far in the chiral Potts models and to translate these results into language more transparent to physicists, so that experts in Monte Carlo calculations, high and low temperature expansions, and various other methods, can use them. We shall pay special attention to the interfacial tension ǫr between the k state and the k −r state. By examining the ground states, it is seen that the integrable line ends at a superwetting point, on which the relation ǫr = r ǫ1 is satisfied, so that it is energetically neutral to have one interface or more. We present also some partial results on the meaning of the integrable line for low temperatures where it lives in the non-wet regime. We make Baxter’s exact results more explicit for the symmetric case. By performing a Bethe Ansatz calculation with open boundary conditions we confirm a dilogarithm identity for the low-temperature expansion which may be new. We propose a new model for numerical studies. This model has only two variables and exhibits commensurate and incommensurate phase transitions and wetting transitions near zero temperature. It appears to be not integrable, except at one point, and at each temperature there is a point, where it is almost identical with the integrable chiral Potts model. KEY WORDS: Chiral Potts model; chiral clock model; star-triangle equations; Yang-Baxter equations; interfacial tension; wetting; superwetting;
We apply a newly proposed Monte Carlo method, the Wang-Landau algorithm, to the study of the three-dimensional antiferromagnetic q -state Potts models on a simple cubic lattice. We systematically study the phase transition of the models with q = 3, 4, 5 and 6. We obtain the finite-temperature phase transition for q = 3 and 4, whereas the transition temperature is down to zero for q = 5. For q = 6 there exists no order for any temperature. We also study the ground-state properties. The size dependence of the ground-state entropy is investigated. We find that the ground-state entropy is larger than the contribution from the typical configurations of the broken-sublattice-symmetry state for q = 3. The same situations are found for q = 4, 5 and 6.