Table of Contents
Physics Research International
Volume 2010 (2010), Article ID 284231, 5 pages
Research Article

Bicritical Central Point of Ising Model Phase Diagram

1Laboratoire de Physique de la Matière Condensée, Faculté des Sciences Ben M'sik, Université Hassan II, Casablanca, Morocco
2Laboratoire de Physique de la Matière Condensée, Faculté des Sciences El Jadida, Université Chouaib Doukkali, Casablanca, Morocco

Received 16 September 2010; Accepted 6 November 2010

Academic Editor: Harold Zandvliet

Copyright © 2010 Yahia Boughaleb et al. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.


We deal with a 2D half occupied square lattice with repulsive interactions between first and second neighboring particles. Despite the intensive studies of the present model the central point of the phase diagram for which the ratio of the two interaction strengths is still open. In the present paper we show, using standard Monte Carlo calculations, that the situation corresponds to a phase of mixed ordered structures quantified by an “algebraic” order parameter defined as the sum of densities of the existing ordered clusters. The introduced grandeur also characterizes the transitions towards the known pure ordered phases for the other values of as mentioned by the agreement of our results with those of the literature. The computation of the Cowley short range order parameter against suggests that the central point is bicritical and is a state to cross when passing between the two pure phases.

1. Introduction

The study of frustrated magnetic systems always arouses the interest [17]. It is known that the nature of the ordering at low temperatures regime of these materials depends on the concentration of the magnetic species and on the relative magnitude of competing energy scales [8, 9]. The Ising model offers a valuable framework for dealing with the phenomena involved as has been confirmed by recent discovery of high- superconductivity in pnictides [10, 11] (the model helps in dealing with the magnetism of the square Fe sublattice). However, when the interactions between first and second neighboring particles are repulsive in a square lattice the approximations adopted to solve the Ising model give rise to a controversial phase diagram [1218]. In fact nowadays everyone agrees that the transition line from the disordered to the ordered phase is continuous when the ratio between the second “" and the first “" neighboring interactions . The situation where the degenerate ordered phase holds is still under debate, particularly the extent of interaction rate for which the transition is first order [1, 2, 1417, 19]. In both cases, the characterization of phase transition is done by means of order parameters that quantify the breakdown of lattice occupation symmetry. This cannot be applied to the system with as to this interaction strength corresponds three ground states: the two ordered phases and and a third “2d-honeycomb " structure joining a central particle to one first and to two second neighboring particles.

Numerical calculations using relaxation processes based on the jump of unique particles are far from being sufficient to alter the stability of clusters formed. To overcome the status quo and reach a stable low temperature pure phase, it was suggested to improve the jump algorithm by considering parallel tempering method [20]. Here the finite size scaling analysis on the obtained effective critical temperatures concludes an absence of phase transition in the thermodynamic limit [14, 15]. In the present, we argue that no pure phase holds rather there are clusters with different ordered structures and sizes that interpenetrate each other well. In other words, the low temperature phase lets appear a kind of long range order with mixed structures. By expecting the local environment of a given site, one easily recognizes the ordered structure's the site belongs to and calculates the density of each ordered structure. The sum of densities of the existing structures is a nonconservative stochastic variable we introduce to play the role of an order parameter. Standard Monte Carlo calculations using the Metropolis algorithm are performed and conclude the existence of a nonvanishing critical temperature when freezing the initial disordered lattice. We then generalize the application of the present grandeur to all values of ratio and discuss the agreement with results of the literature. Next we consider the Cowley Short Range Order Parameter “SROP" that we compute versus the strength of the competing interaction energies at fixed temperatures. The calculations lead to the same results except for the central point where the discontinuity in the behavior of SROP reflects that the phase transition is first order. We notice here that the bicriticality does not originate from mixture of different particles as in references [21, 22].

2. Model and Method

The study is performed in the framework of lattice gas model. We consider a bidimensional square lattice with lateral size for which interaction energy is insured by the following Hamiltonian of static interactions:" is a Boolean occupation number, denotes the ratio between interaction energies of second “” and first “” neighboring particles, and is the chemical potential and number of angle brackets indicates number of the neighbors.

Adsorption/desorption processes on periodically bounded lattices are insured meaning standard Metropolis transition rate applied to the state of single site of the latticedenotes the inverse reduced temperature and corresponds to the energy difference between the final and the initial configurations. We recall that in the Metropolis algorithm a site is allowed to change its state if the global energy of the system is reduced or as a response to a local deformation of the lattice represented by the comparison of the adsorption/desorption length to a number chosen at random. For each given size of the lattice the relaxation is controlled by computing at fixed temperature the energy, the average concentration, and the introduced algebraic order parameter “AOP". The later grandeur characterizes the sum of densities of the ordered clusters on the lattice , where “" designates the structure and / . Here the summation is performed once constraints relative to establishment of a type of order are fulfilled (e.g., a site belongs to the ordered structure whenever its first neighboring sites are all empty) Brackets designate the neighbours and the number of the brackets the number of the neighbours and the symbols and denote perpendicular and parallel direction of the given site.

3. Results and Discussions

Figure 1 reports adsorption isotherm as well as the behavior of AOP for different temperature regimes. For high temperatures adsorption isotherm presents the known Langmuir gas behavior. Absence of ordering is characterized by weak values of AOP. This can be understood by existence of ordered clusters of limited sizes as a consequence of the equal probability of occupying lattice sites.

Figure 1: Variation of adsorption isotherms and density of ordered structures versus chemical potential for different interaction regimes at . Size of the lattice is fixed to .

To the transition to ordering marked by plateau in adsorption isotherm corresponds a clear increase of AOP values. All sites of the lattice contribute to the appearance of the order as can be seen by the collapse of both adsorption isotherm and AOP curves in the interval corresponding to the coverage fixed at . Signature of this kind of long range order constituted by mixture of different types of ordered clusters is also revealed by the behavior of the susceptibility of the system. The later is defined as a measure of the autocorrelation function of the AOP: In Figure 2 we report the susceptibility for different sizes of the lattice as well as the finite size scaling analysis joining the effective critical temperatures (the inset). Maxima of the curve are related according to the formula This leads to a nonvanishing critical temperature at the thermodynamic limit . The obtained critical exponents suggest that the phase transition is nonuniversal and .

Figure 2: Validity of algebraic order parameter hypothesis: behavior of the susceptibility defined as fluctuations of ordered structures density versus reduced temperature for different sizes of the lattice and finite size scaling analysis (the inset: we notice that the error bars have the same size of the symbols used in the graph).

Another proof of the validity of our hypothesis lies in characterizing the phase transition to all values. The phase diagram of Figure 3 shows that obtained critical temperatures agree with the results of studies within Monte Carlo technique [1215, 23]. Once again finite size scaling analysis confirms an Ising-like behavior with universal critical exponents of the system with and that order/disorder phase transition for is nonuniversal and for for example, the value has to be compared with the one in [12, 13]).

Figure 3: Phase Diagram of the half full lattice obtained by means of the algebraic order parameter fluctuations.

Now that transition to ordered phase at has nonnull temperature in the thermodynamic limit, it would be interesting to characterize transition from ordered phase to one as a response of varying the relative magnitude of competing interaction energies. For this reason we compute the Cowley [24] short range order parameter versus intensity of ratio for different temperatures of the system. The situation corresponds to crossing over the phase diagram. One distinguishes plateaus in the low temperature regime (Figure 4) that quantify the appareance of ordered structures in the lattice with a structure when and a phase for .

Figure 4: Behavior of short range order parameter while crossing the phase diagram. The curve shows bicritical character of the remarkable phase diagram central point.

Transition to the disordered phase is signalled by the intersection between the critical curve and the line of constant temperature. Extent of the plateaus increases by cooling while phase with mixed structures constitutes remarkable point. Jump in short range order parameter reflects the transition is first order. Then the central point presents bicritical behaviour and constitutes the only one way between the two ordered pure phases.

4. Conclusion

We can conclude from the results obtained that within the definition of the AOP it is possible to envisage—for the frustrated square lattice with repulsive interactions—a transition to an ordered phase with a mixture of structures constituted by ground state configurations of a same energy. This phase transition presents a bicritical behavior depending on whether the establishment phase is carried out by cooling at or by varying the ratio of exchange constants at fixed temperature. AOP grandeur can also be applied to characterize ordinary ordered/disordered phase transition. The critical universal behavior for is found. For the calculations lead to nonuniversal critical behavior and support then the absence of any first-order phase transition.


The authors are grateful to Pr. J. Prost for fruitful discussions and suggestions.


  1. A. Malakis, S. S. Martinos, I. A. Hadjiagapiou, N. G. Fytas, and P. Kalozoumis, “Entropic sampling via Wang-Landau random walks in dominant energy subspaces,” Physical Review E, vol. 72, no. 6, Article ID 066120, 2005. View at Publisher · View at Google Scholar
  2. A. Malakis, A. N. Berker, I. A. Hadjiagapiou, and N. G. Fytas, “Strong violation of critical phenomena universality: Wang-Landau study of the two-dimensional Blume-Capel model under bond randomness,” Physical Review E, vol. 79, no. 1, Article ID 011125, 2009. View at Publisher · View at Google Scholar
  3. P. Carretta, M. Filibian, R. Nath, C. Geibel, and P. J. C. King, “Fluctuations and correlations in a frustrated S = 1/2 square lattice with competing ferromagnetic and antiferromagnetic interactions studied by muon-spin relaxation,” Physical Review B, vol. 79, no. 22, Article ID 224432, 2009. View at Publisher · View at Google Scholar · View at Scopus
  4. N. G. Fytas, A. Malakis, and I. Georgiou, “Wang-Landau study of the random bond square Ising model with nearest-and next-nearest-neighbor interactions,” Journal of Statistical Mechanics: Theory and Experiment, vol. 2008, no. 7, Article ID L07001, 2008. View at Publisher · View at Google Scholar · View at Scopus
  5. O. R. Salmon, J. R. de Sousa, and F. D. Nobre, “Phase diagram of the two-dimensional Ising model with random competing interactions,” Physics Letters A, vol. 373, no. 30, pp. 2525–2528, 2009. View at Publisher · View at Google Scholar · View at Scopus
  6. L. M. Volkova and S. A. Polyshchuk, “The study of nearest- and next-nearest-neighbour magnetic interactions in seven tetragonal compounds of V(IV) containing linear chains and square lattices,” Journal of Physics Condensed Matter, vol. 18, no. 49, pp. 11177–11193, 2006. View at Publisher · View at Google Scholar · View at Scopus
  7. A. A. Tsirlin and H. Rosner, “Extension of the spin- 1/2 frustrated square lattice model: the case of layered vanadium phosphates,” Physical Review B, vol. 79, no. 21, Article ID 214417, 2009. View at Publisher · View at Google Scholar · View at Scopus
  8. H. V. Löhneysen, A. Rosch, M. Vojta, and P. Wölfle, “Fermi-liquid instabilities at magnetic quantum phase transitions,” Reviews of Modern Physics, vol. 79, no. 3, pp. 1015–1075, 2007. View at Publisher · View at Google Scholar · View at Scopus
  9. S. T. Bramwell, M. J. P. Gingras, and P. C. W. Holdsmorth, “Spin Ice,” in Frustrated Spin Systems, H. T. Diep, Ed., p. 367, World Scientific, Singapore, 2005. View at Google Scholar
  10. Y. Kamihara, T. Watanabe, M. Hirano, and H. Hosono, “Iron-based layered superconductor La[O1xFx]FeAs(x=0.050.12) with Tc=26 K,” Journal of the American Chemical Society, vol. 130, no. 11, pp. 3296–3297, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  11. X. H. Chen, T. Wu, G. Wu, R. H. Liu, H. Chen, and D. F. Fang, “Superconductivity at 43 K in SmFeAsO1xFx,” Nature, vol. 453, no. 7196, pp. 761–762, 2008. View at Publisher · View at Google Scholar · View at PubMed · View at Scopus
  12. D. P. Landau, “Phase transitions in the Ising square lattice with next-nearest-neighbor interactions,” Physical Review B, vol. 21, no. 3, pp. 1285–1297, 1980. View at Publisher · View at Google Scholar
  13. K. Binder and D. P. Landau, “Phase diagrams and critical behavior in Ising square lattices with nearest- and next-nearest-neighbor interactions,” Physical Review B, vol. 21, no. 5, pp. 1941–1962, 1980. View at Publisher · View at Google Scholar
  14. A. Kalz, A. Honecker, S. Fuchs, and T. Pruschke, “Phase diagram of the Ising square lattice with competing interactions,” European Physical Journal B, vol. 65, no. 4, pp. 533–537, 2008. View at Publisher · View at Google Scholar
  15. A. Kalz, A. Honecker, S. Fuchs, and T. Pruschke, “Monte Carlo studies of the ising square lattice with competing interactions,” in Proceedings of the Conference on Highly Frustrated Magnetism, Journal of Physics: Conference Series 145, (2009)012051, September 2008.
  16. A. Malakis, P. Kalozoumis, and N. Tyraskis, “Monte Carlo studies of the square Ising model with next-nearest-neighbor interactions,” European Physical Journal B, vol. 50, no. 1-2, pp. 63–67, 2006. View at Publisher · View at Google Scholar
  17. J. L. Monroe and S.-Y. Kim, “Phase diagram and critical exponent ν for the nearest-neighbor and next-nearest-neighbor interaction Ising model,” Physical Review E, vol. 76, no. 2, Article ID 021123, 2007. View at Publisher · View at Google Scholar
  18. R. A. dos Anjos, J. Roberto Viana, and J. Ricardo de Sousa, “Phase diagram of the Ising antiferromagnet with nearest-neighbor and next-nearest-neighbor interactions on a square lattice,” Physics Letters A, vol. 372, no. 8, pp. 1180–1184, 2008. View at Publisher · View at Google Scholar
  19. J. L. Moran-Lopez, F. Aguilera-Granja, and J. M. Sanchez, “Phase transitions in Ising square antiferromagnets with first- and second-neighbour interactions,” Journal of Physics: Condensed Matter, vol. 6, no. 45, pp. 9759–9772, 1994. View at Publisher · View at Google Scholar
  20. D. J. Earl and M. W. Deem, “Parallel tempering: theory, applications, and new perspectives,” Physical Chemistry Chemical Physics, vol. 7, no. 23, pp. 3910–3916, 2005. View at Publisher · View at Google Scholar
  21. J. Strečka, “Weak universality, bicritical points and reentrant transitions in the critical behaviour of a mixed spin-1/2 and spin-3/2 Ising model on the Union Jack (centered square) lattice,” Physica Status Solidi (B), vol. 243, no. 3, pp. 708–715, 2006. View at Publisher · View at Google Scholar
  22. J. Strečka, L. Čanová, and J. Dely, “Weak universal critical behaviour of the mixed spin-(1/2, S) Ising model on the Union Jack (centered square) lattice: integer versus half-odd-integer spin-S case,” Physica Status Solidi (B), vol. 243, no. 8, pp. 1946–1955, 2006. View at Publisher · View at Google Scholar
  23. M. Noureddine, R. Nassif, K. Saadouni, and Y. Boughaleb, Inernational journal of Modern Physics B, vol. 24, no. 23, Article ID 05530, 2010.
  24. J. M. Cowley, “An approximate theory of order in alloys,” Physical Review, vol. 77, no. 5, pp. 669–675, 1950. View at Publisher · View at Google Scholar