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.     

Tunable magnetic anisotropy of antiferromagnetic superlattice and resultant exchange bias of ferromagnetic layer on it

Exchange biasing of ferromagnetic layer deposited on the antiferromagnetic superlattice was investigated in (Co₇₀Fe₃₀/Ru)_29.5/Ru/Co₉₀Fe₁₀ multilayers. Uniaxial magnetic anisotropy (K_AF) was induced and tuned in the antiferromagentic superlattice by uniaxial substrate bending method through the inverse effect of magnetostriction. The exchange bias increased and tended to be saturated with increasing the K_AF, while it was not observed at K_AF=0.     

Third and fourth order phase transitions: Exact solution for the Ising model on the Cayley tree

In this work we present the first exact solution of a system of interacting particles with phase transitions of order higher than two. The presented analytical derivation shows that the Ising model on the Cayley tree exhibits a line of third order phase transition points, between temperatures and , and a line of fourth order phase transitions between T_BP and ∞, where k_B is the Boltzmann constant, and J is the nearest-neighbor interaction parameter.     

Phase diagram of uniaxial antiferromagnetic particles: Field perpendicular to the easy axis

We consider a spherical uniaxial antiferromagnetic particle in the presence of an external magnetic field perpendicular to its easy axis. The model is described by a classical Heisenberg Hamiltonian including a single-ion uniaxial anisotropy, where the magnetic moments of the particle are represented by continuous spin vectors. We employ mean-field calculations and Monte Carlo simulations to determine the phase diagram of the system. The phase diagram in the plane field versus temperature is obtained for particles with radii ranging from three up to twelve spacing lattice units. We have seen that a particle with more than nine shells behaves as a true thermodynamic system. We find the explicit dependence of the zero temperature critical field and the Néel temperature on the diameter of the particle. At low temperatures, we have also shown that, for particles with three or more shells, the critical field follows a T² law, which is in agreement with the predictions of the spin-wave theory, when the field is perpendicular to the easy axis.     

Magnetic properties of the mixed ferrimagnetic ternary system with a single-ion anisotropy on the Bethe lattice

The magnetic properties of the ternary system ABC consisting of spins , S=1, and are investigated on the Bethe lattice by using the exact recursion relations. We consider both ferromagnetic and antiferromagnetic exchange interactions. The exact expressions for magnetizations and magnetic susceptibilities are found, and thermal behaviors of magnetizations and susceptibilities are studied. We construct the phase diagrams and find that the system exhibits one, two or even three compensation temperatures depending on the values of the interaction parameters in the Hamiltonian. Moreover, the system undergoes a second-order phase transition for the coordination number q?3 and a second- and first-order phase transitions for q>3; hence the system gives a tricritical point. The system also exhibits the reentrant behaviors.     

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.     

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.     

Effective-field theory of the Ising model with three alternative layers on the honeycomb and square lattices

The Ising model with three alternative layers on the honeycomb and square lattices is studied by using the effective-field theory with correlations. We consider that the nearest-neighbor spins of each layer are coupled ferromagnetically and the adjacent spins of the nearest-neighbor layers are coupled either ferromagnetically or anti-ferromagnetically depending on the sign of the bilinear exchange interactions. We investigate the thermal variations of the magnetizations and present the phase diagrams. The phase diagrams contain the paramagnetic, ferromagnetic and anti-ferromagnetic phases, and the system also exhibits a tricritical behavior.     

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.     

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.     

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.     

Partition function zeros of the antiferromagnetic Ising model on triangular lattice in the complex temperature plane for nonzero magnetic field

The grand partition functions Z(T,B)$Z(T,B)$ of the Ising model on L×L$L\times L$ triangular lattices with fully periodic boundary conditions, as a function of temperature T and magnetic field B  , are evaluated exactly for L<12$L<12$ (using microcanonical transfer matrix) and approximately for L?12$L?12$ (using Wang–Landau Monte Carlo algorithm). From Z(T,B)$Z(T,B)$, the distributions of the partition function zeros of the triangular-lattice Ising model in the complex temperature plane for real B≠0$B\ne 0$ are obtained and discussed for the first time. The critical points aN(x)$a_N(x)$ and the thermal scaling exponents yt(x)$y_t(x)$ of the triangular-lattice Ising antiferromagnet, for various values of x=e−2βB$x=e^{-2\beta B}$, are estimated using the partition function zeros.     

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.     

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.     

Local analysis of frustration based on Kagomé lattices

Magnetic frustration in the framework of the Edwards-Anderson model is studied for the ground level (T=0) of the Kagomé lattice (KL). A sample consists of a realization of a random distribution of ferromagnetic (F) and antiferromagnetic (AF) interactions of the same strength along lines connecting nearest neighbors. Our goal is to compare two methods to calculate the following parameters: ground state energy, frustration and average frustration segment as functions of x, the concentration of F bonds. In doing so we make use of topological concepts such as plaquettes and frustration segments. The probability of a plaquette being unfrustrated (or flat) is ℘_p(x), while the probability of a plaquette being frustrated (or curved) is ℘_c(x). The analysis is done locally on a representative portion of the lattice which is called cell; cells of two different sizes are used in the present work. One method (which is simpler) is based on the probability function of any plaquette configuration ?(℘_p,℘_c). The other method is more exact but also more complex and increasingly difficult to use for large cells; it is based on the probability of bond configurations ψ(x). These methods are compared between themselves noting that the simpler method can be enough for most of the range for x. In addition, numerical simulations for many random samples at different concentrations x for a size given by 75 spins with periodic boundary conditions were done. This provides reference lines to compare with the properties under study. The local frustration analysis to obtain both ?(℘_p,℘_c) and ψ(x) is done over two cells of different size. Robustness of the criteria used in the local frustration analysis is also investigated.     

Thermodynamics of the quantum and classical Ising models with skew magnetic field

The thermodynamics of the unitary (normalized spin) quantum and classical Ising models with skew magnetic field, for |J|β?0.9, is derived for the ferromagnetic and antiferromagnetic cases. The high-temperature expansion (β-expansion) of the Helmholtz free energy is calculated up to order β⁷ for the quantum version (spin S≥1/2) and up to order β¹⁹ for the classical version. In contrast to the S=1/2 case, the thermodynamics of the transverse Ising and that of the XY model for S>1/2 are not equivalent. Moreover, the critical line of the T=0 classical antiferromagnetic Ising model with skew magnetic field is absent from this classical model, at least in the temperature range of |J|β?0.9.     

Frequency dependence of the complex susceptibility for a spin-1 Ising model

The complex susceptibility or the dynamic susceptibility (χ(ω)=χ′(ω)−iχ″(ω)) for a spin-1 Ising system with bilinear and biquadratic interactions is obtained on the basis of Onsager theory of irreversible processes. If the logarithm of the susceptibilities is plotted as a function of the logarithm of frequency, then the real part (χ′) displays a sequence of plateau regions and the imaginary part (χ″) has a sequence of maxima in the ordered or ferromagnetic phase. On the other hand, only one plateau region in χ′ and one maximum in χ″ is observed in the disordered or paramagnetic phase. Argand or Cole-Cole plots (χ″−χ′) for a selection of temperatures are also shown, and a sequence of semicircles is illustrated in the ordered phase and only one semicircle for the disordered phase in these plots.     

Phase diagrams of a nonequilibrium mixed spin-3/2 and spin-2 Ising system in an oscillating magnetic field

The phase diagrams of the nonequilibrium mixed spin-3/2 and spin-2 Ising ferrimagnetic system on square lattice under a time-dependent external magnetic field are presented by using the Glauber-type stochastic dynamics. The model system consists of two interpenetrating sublattices of spins σ=3/2 and S=2, and we take only nearest-neighbor interactions between pairs of spins. The system is in contact with a heat bath at absolute temperature T_abs and the exchange of energy with the heat bath occurs via one-spin flip of the Glauber dynamics. First, we investigate the time variations of average order parameters to find the phases in the system and then the thermal behavior of the dynamic order parameters to obtain the dynamic phase transition (DPT) points as well as to characterize the nature (first- or second-order) phase transitions. The dynamic phase diagrams are presented in two different planes. Phase diagrams contain paramagnetic (p), ferrimagnetic (i₁, i₂, i₃) phases, and three coexistence or mixed phase regions, namely i₁+p, i₂+p and i₃+p mixed phases that strongly depend on interaction parameters.     

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.