Phase transitions in self-dual generalizations of the Baxter–Wu model

We study two types of generalized Baxter–Wu models, by means of transfer-matrix and Monte Carlo techniques. The first generalization allows for different couplings in the up- and down-triangles, and the second generalization is to a q-state spin model with three-spin interactions. Both generalizations lead to self-dual models, so that the probable locations of the phase transitions follow. Our numerical analysis confirms that phase transitions occur at the self-dual points. For both generalizations of the Baxter–Wu model, the phase transitions appear to be discontinuous.     

Phase diagram of the two-dimensional Ising model with random competing interactions

An Ising model with ferromagnetic nearest-neighbor interactions J₁ (J₁>0) and random next-nearest-neighbor interactions [+J₂ with probability p and −J₂ with probability (1−p); J₂>0] is studied within the framework of an effective-field theory based on the differential-operator technique. The order parameters are calculated, considering finite clusters with n=1,2, and 4 spins, using the standard approximation of neglecting correlations. A phase diagram is obtained in the plane temperature versus p, for the particular case J₁=J₂, showing both superantiferromagnetic (low p) and ferromagnetic (higher values of p) orderings at low temperatures.     

Phase diagram of the Ising antiferromagnet with nearest-neighbor and next-nearest-neighbor interactions on a square lattice

The phase diagram of the Ising model in the presence of nearest-neighbor (J₁) and next-nearest-neighbor (J₂) interactions on a square lattice is studied within the framework of the differential operator technique. The Hamiltonian is solved by effective-field theory in finite cluster (we have chosen N=4 spins). We have proposed a functional for the free energy (similar to Landau expansion) to obtain the phase diagram in the (T,α) space (α=J₂/J₁), where the transition line from the superantiferromagnetic (SAF) to the paramagnetic (P) phase is of first-order in the range 1/2<α<0.95 in contrast to previous study of CVM (Cluster Variational Method) that predict first-order transition for α=1.0. Our results for α=1.0 are in accordance with MC (Monte Carlo) simulations, that predict a second-order transition.     

Cluster evaporation transition in finite system

Using Wang-Landau entropic sampling we study the Ising model in the framework of microcanonical ensemble (fixed magnetization). We are working for lattice size up to 1500×1500 in two dimensions and 100×100×100 in three dimensions. As we approach the coexistence curve from inside, varying temperature and keeping the magnetization constant, a first-order phase transition takes place for a temperature near the coexistence curve if the lattice size is large enough. We analyze various features of this transition as well as the scaling behavior of characteristic quantities and we compare our numerical results with existing theories.     

Nonequilibrium phase transition in the kinetic Ising model on a two-layer square lattice under the presence of an oscillating field

The nonequilibrium or dynamic phase transitions are studied, within a mean-field approach, in the kinetic Ising model on a two-layer square lattice consisting of spin- 1/2 ions in the presence of a time varying (sinusoidal) magnetic field has been studied by using Glauber-type stochastic dynamics. The dynamic equations of motion are obtained in terms of the intralayer coupling constants J₁ and J₂ for the first and second layer, respectively, and interlayer coupling constant J₃ between these two layers. The nature (first- or second-order) of the transitions is characterized by investigating the behavior of the thermal variations of the dynamic order parameters. The dynamic phase transitions are obtained and the dynamic phase diagrams are constructed in the plane of the reduced temperature versus the amplitude of the magnetic field and found fourteen fundamental types of phase diagrams. Phase diagrams exhibit one, two or three dynamic tricritical points for various values of J₂/|J₁| and J₃/|J₁|. Besides the paramagnetic (p), ferromagnetic (f) and compensated (c) phases, there were the f+c,f+sf,c+sf,af+p,m+p,f+m and c+af, where the af, sf and m are the antiferromagnetic, surface ferromagnetic and mixed phases respectively. Coexistence phase regions also exist in the system.     

Exact results of potts model on a special kind of Sierpinski carpets

With a way different from renormalization group method and graph theory [1, 2], we have calculated the exact partition function, free energy and spin-spin correlation function of Potts model, other than the Ising model, on a special Sierpinski Carpet (SC). The results indicate no phase transition occurs at any finite temperature for the fractal with finite R(the order of ramification) and thus consist with the conclusion produced by renormalization group method and other physical arguments.     

First order phase transitions in the canonical and the microcanonical ensemble

In the canonical ensemble any singularity of a thermodynamic function at a temperatureT _c is smeared over a temperature range of orderT _T /N. Therefore it is rather difficult to distinguish between a discontinuous and a continuous phase transition on the basis of numerical data obtained for finite systems in the canonical ensemble. It is demonstrated for four model systems that this problem cannot be circumvented by considering higher cumulants of the energy distribution or cumulant ratios. On the other hand, the distinction between first and a second order phase transition is rather direct if based on the microcanonical density of states which is readily obtainable in the dynamical ensemble.     

Monte Carlo analysis of critical properties of the two-dimensional randomly site-diluted Ising model via Wang-Landau algorithm

The influence of random site dilution on the critical properties of the two-dimensional Ising model on a square lattice was explored by Monte Carlo simulations with the Wang-Landau sampling. The lattice linear size was L=20-120 and the concentration of diluted sites q=0.1,0.2,0.3. Its pure version displays a second-order phase transition with a vanishing specific heat critical exponent α, thus, the Harris criterion is inconclusive, in that disorder is a relevant or irrelevant perturbation for the critical behaviour of the pure system. The main effort was focused on the specific heat and magnetic susceptibility. We have also looked at the probability distribution of susceptibility, pseudocritical temperatures and specific heat for assessing self-averaging. The study was carried out in appropriate restricted but dominant energy subspaces. By applying the finite-size scaling analysis, the correlation length exponent ν was found to be greater than one, whereas the ratio of the critical exponents (α/ν) is negative and (γ/ν) retains its pure Ising model value supporting weak universality.     

Kinetics of a mixed spin-1/2 and spin-3/2 Ising ferrimagnetic model

We present a study, within a mean-field approach, of the kinetics of a mixed ferrimagnetic model on a square lattice in which two interpenetrating square sublattices have spins that can take two values, , alternated with spins that can take the four values, . We use the Glauber-type stochastic dynamics to describe the time evolution of the system with a crystal-field interaction in the presence of a time-dependent oscillating external magnetic field. The nature (continuous and discontinuous) of transition is characterized by studying the thermal behaviors of average order parameters in a period. The dynamic phase transition points are obtained and the phase diagrams are presented in the reduced magnetic field amplitude (h) and reduced temperature (T) plane, and in the reduced temperature and interaction parameter planes, namely in the (h, T) and (d, T) planes, d is the reduced crystal-field interaction. The phase diagrams always exhibit a tricritical point in (h, T) plane, but do not exhibit in the (d, T) plane for low values of h. The dynamic multicritical point or dynamic critical end point exist in the (d, T) plane for low values of h. Moreover, phase diagrams contain paramagnetic (p), ferromagnetic (f), ferrimagnetic (i) phases, two coexistence or mixed phase regions, (f+p) and (i+p), that strongly depend on interaction parameters.     

Magnetization and magnetic susceptibility of the Ising ferromagnetic/antiferromagnetic superlattice

A linear cluster mean-field approximation is used to study the magnetic properties of the Ising ferromagnetic/antiferromagnetic superlattice, which is composed of a spin-1/2 ferromagnetic monolayer and a spin-1 antiferromagnetic monolayer with a single-ion anisotropy alternatively. By using the transfer matrix method, we calculate the magnetization and the initial magnetic susceptibility as functions of temperature for different interlayer coupling, single-ion anisotropy. We summarize the changing behaviors of the spin structure in ferromagnetic and antiferromagnetic layers and the characteristics of the corresponding magnetic susceptibilities, give the transition temperature as a function of the interlayer exchange coupling for different single-ion anisotropy, and analyze the features of the magnetization and the magnetic susceptibility.     

Exact results of a dimerized S=1 Ising chain with single-ion anisotropy

The dimerized spin-1 Ising chain with both longitude and transverse single-ion anisotropies D_z and D_x is solved exactly by means of a mapping to the spin- Ising chain with the alternating transverse fields and the Jordan-Wigner transformation. The analytical expressions of the quasi-particles’ spectra Λ_k, the minimal energy gap Δ₀ for exciting a fermion quasi-particle, the minimal energy gap Δ_h for exciting a hole, and the ground-state energy E_g are obtained. The phase diagram of the ground state is also given. The results show that the system exhibits a series of quantum phase transitions depending on the dimerization strength of the crystal fields, while the quantum critical points are determined exactly.     

Dynamic compensation temperature in the kinetic spin-1 Ising model in an oscillating external magnetic field on alternate layers of a hexagonal lattice

The dynamic behavior of a two-sublattice spin-1 Ising model with a crystal-field interaction (D) in the presence of a time-varying magnetic field on a hexagonal lattice is studied by using the Glauber-type stochastic dynamics. The lattice is formed by alternate layers of spins σ=1 and S=1. For this spin arrangement, any spin at one lattice site has two nearest-neighbor spins on the same sublattice, and four on the other sublattice. The intersublattice interaction is antiferromagnetic. We employ the Glauber transition rates to construct the mean-field dynamical equations. Firstly, we study time variations of the average magnetizations in order to find the phases in the system, and the temperature dependence of the average magnetizations in a period, which is also called the dynamic magnetizations, to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (continuous and discontinuous) of transitions. Then, the behavior of the total dynamic magnetization as a function of the temperature is investigated to find the types of the compensation behavior. Dynamic phase diagrams are calculated for both DPT points and dynamic compensation effect. Phase diagrams contain the paramagnetic (p) and antiferromagnetic (af) phases, the p+af and nm+p mixed phases, nm is the non-magnetic phase, and the compensation temperature or the L-type behavior that strongly depend on the interaction parameters. For D<2.835 and H₀>3.8275, H₀ is the magnetic field amplitude, the compensation effect does not appear in the system.     

The Glauber dynamics for a spin-1 metamagnetic Ising system with bilinear and biquadratic interactions

We present a study, within a mean-field approximation, of the dynamics of a spin-1 metamagnetic Ising system with bilinear and biquadratic interactions in the presence of a time-dependent oscillating external magnetic field. First, we employ the Glauber transition rates to construct the set of mean-field dynamic equations. Then, we study the time variation of the average order parameters to find the phases in the system. We also investigate the thermal behavior of dynamic order parameters to characterize the nature (first- or second-order) of the dynamic transitions. The dynamic phase transitions are obtained and the phase diagrams are constructed in two different the planes. The phase diagrams contain a disordered and ordered phases, and four different mixed phases that strongly depend on interaction parameters. Phase diagrams also display one or two dynamic tricritical points, a dynamic double critical end and dynamic quadruple points. A comparison is made with the results of the other metamagnetic Ising systems.     

Mean field study of the mixed Ising model in a random crystal field

The magnetic properties of a mixed Ising ferrimagnetic system, in which the two interacting sublattices have spins σ, (±1/2) and spins S, (±1,0) in the presence of a random crystal field, have been studied with the mean field approach. The obtained results show the existence of some interesting phenomena, such as the appearance of a new ferrimagnetic phase, namely, partly ferrimagnetic phase and consequently the existence of four topologically different types of phase diagrams. Furthermore, compensation behaviour and re-entrant phenomenon are found for appropriate ranges of crystal field. Thermal magnetization behaviour and phase diagrams have been discussed in detail.     

Phase diagram of the three states Potts model with next nearest neighbour interactions on the Bethe lattice

We have found an exact phase diagram of the Potts model with competing nearest neighbor and next nearest neighbor interactions on the Bethe lattice of order two. The diagram consists of five phases: ferromagnetic, paramagnetic, modulated, antiphase and paramodulated, all meeting at the multicritical point . We report on a new phase which we denote as paramodulated, found at low temperatures and characterized by zero average magnetization lying inside the modulated phase. Such a phase, inherent in the Potts model has no analogues in the Ising setting.     

Classical two-site fidelity and quantum phase transitions in the Ising model and the extended Hubbard model

In this Letter, the classical two-site-ground-state fidelity (CTGF) is exploited to identify quantum phase transitions (QPTs) for the transverse field Ising model (TFIM) and the one-dimensional extended Hubbard model (EHM). Our results show that the CTGF exhibits an abrupt change around the regions of criticality and can be used to identify QPTs in spin and fermionic systems. The method is especially convenient when it is connected with the density-matrix renormalization group (DMRG) algorithm.     

Stability conditions for fermionic Ising spin-glass models in the presence of a transverse field

The stability of a spin-glass (SG) phase is analyzed in detail for a fermionic Ising SG (FISG) model in the presence of a magnetic transverse field Γ. The fermionic path integral formalism, replica method and static approach have been used to obtain the thermodynamic potential within one step replica symmetry breaking ansatz. The replica symmetry (RS) results show that the SG phase is always unstable against the replicon. Moreover, the two other eigenvalues λ_± of the Hessian matrix (related to the diagonal elements of the replica matrix) can indicate an additional instability to the SG phase, which enhances when Γ is increased. Therefore, this result suggests that the study of the replicon cannot be enough to guarantee the RS stability in the present quantum FISG model, especially near the quantum critical point. In particular, the FISG model allows changing the occupation number of sites, so one can get a first order transition when the chemical potential exceeds a certain value. In this region, the replicon and the λ_± indicate instability problems for the SG solution close to all ranges of a first order boundary.     

Phase diagrams of a nonequilibrium mixed spin-1/2 and spin-2 Ising ferrimagnetic system under a time-dependent oscillating magnetic field

We present phase diagrams for a nonequilibrium mixed spin-1/2 and spin-2 Ising ferrimagnetic system on a square lattice in the presence of a time dependent oscillating external magnetic field. We employ the Glauber transition rates to construct the mean-field dynamical equations. The time variation of the average magnetizations and the thermal behavior of the dynamic magnetizations are investigated, extensively. The nature (continuous or discontinuous) of the transitions is characterized by studying the thermal behaviors of the dynamic magnetizations. The dynamic phase transition points are obtained and the phase diagrams are presented in two different planes. Phase diagrams contain paramagnetic (p) and ferrimagnetic (i) phases, and one coexistence or mixed phase region, namely the i+p, that strongly depend on interaction parameters. The system exhibits the dynamic tricritical point and the reentrant behaviors.     

A flower-like Ising model. Thermodynamic properties

We consider a flower-like Ising model, in which there are some additional bonds (in the “flower-core”) compared to a pure Ising chain. To understand the behaviour of this system and particularly the competition between ferromagnetic (usual) bonds along the chain and antiferromagnetic (additional) bonds across the chain, we study analytically and iteratively the main thermodynamic quantities. Very interesting is, in the zero-field and zero-temperature limit, the behaviour of the magnetization and the susceptibility, closely related to the ground state configurations and their degeneracies. This degeneracy explains the existence of non-zero entropy at zero temperature, in our results. Also, this model could be useful for the experimental investigations in studying the saturation curves for the enzyme kinetics or the melting curves for DNA-denaturation in some flower-like configurations.