Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices

Phase transitions of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on several decorated planar lattices consisting of interconnected diamonds are investigated within the framework of the generalized decoration-iteration transformation. The main attention ispaid to the systematic...

Ausführliche Beschreibung

Gespeichert in:
Bibliographische Detailangaben
Veröffentlicht in:Condensed Matter Physics
Datum:2011
Hauptverfasser: Gálisová, L., Strečka, J.
Format: Artikel
Sprache:English
Veröffentlicht: Інститут фізики конденсованих систем НАН України 2011
Online Zugang:https://nasplib.isofts.kiev.ua/handle/123456789/119900
Tags: Tag hinzufügen
Keine Tags, Fügen Sie den ersten Tag hinzu!
Назва журналу:Digital Library of Periodicals of National Academy of Sciences of Ukraine
Zitieren:Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices / Lucia Gálisová, Jozef Strečka // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13002: 1-11. — Бібліогр.: 39 назв. — англ.

Institution

Digital Library of Periodicals of National Academy of Sciences of Ukraine
id nasplib_isofts_kiev_ua-123456789-119900
record_format dspace
spelling Gálisová, L.
Strečka, J.
2017-06-10T10:43:44Z
2017-06-10T10:43:44Z
2011
Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices / Lucia Gálisová, Jozef Strečka // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13002: 1-11. — Бібліогр.: 39 назв. — англ.
1607-324X
PACS: 05.50.+q, 75.10.Hk, 75.10.Jm, 68.35.Rh
DOI:10.5488/CMP.14.13002
arXiv:1106.4687
https://nasplib.isofts.kiev.ua/handle/123456789/119900
Phase transitions of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on several decorated planar lattices consisting of interconnected diamonds are investigated within the framework of the generalized decoration-iteration transformation. The main attention ispaid to the systematic study of the finite-temperature phase diagrams in dependence on the lattice topology. The critical behaviour of the hybrid quantum-classical Ising-Heisenberg model is compared with the relevant behaviour of its semi-classical Ising analogue. It is shown that both models on diamond-like decorated planar lattices exhibit a striking critical behaviour including reentrant phase transitions. The higher the lattice coordination number is, the more pronounced reentrance may be detected.
Дослiджуються фазовi переходи у змiшанiй спiн 1/2 i спiн 1 моделi Iзiнга-Гайзенберга на декiлькох декорованих плоских ґратках, що складаються зi сполучених ромб i в рамках узагальненого декорацiйно-iтерацiйного перетворення. Основна увага придiляється систематичному вивченню фазових дiаграм при скiнчених температурах в залежностi вiд топологiї гратки. Критична поведiнка гiбридної квантово-класичної моделi Iзiнга-Гайзенберга порiвнюється з поведiнкою її напiвкласичного iзінгiвського аналога. Показано, що обидвi моделi на ромбiчноподiбних декорованих плоских ґратках демонструють цiкаву критичну поведiнку, включаючи зворотнi фазовi переходи. Чим вище координацiйне число ґратки, тим помiтнiшою може бути зворотнiсть.
This work was partially supported by the internal grant of P.J. Saf ́arik University under the contract No. VVGS 1/10–11 and by European Union European regional development fond (ERDFEU) grant under the contract No. ITMS26220120005 (activity 3.2.).
en
Інститут фізики конденсованих систем НАН України
Condensed Matter Physics
Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
Фазовi переходи у геометрично фрустрованi й змiшанiй спiн -1/2 i спiн -1 моделi Iзiнга - Гайзенберга на ромбiчноподiбних декорованих плоских ґратках
Article
published earlier
institution Digital Library of Periodicals of National Academy of Sciences of Ukraine
collection DSpace DC
title Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
spellingShingle Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
Gálisová, L.
Strečka, J.
title_short Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
title_full Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
title_fullStr Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
title_full_unstemmed Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices
title_sort phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 ising-heisenberg model on diamond-like decorated planar lattices
author Gálisová, L.
Strečka, J.
author_facet Gálisová, L.
Strečka, J.
publishDate 2011
language English
container_title Condensed Matter Physics
publisher Інститут фізики конденсованих систем НАН України
format Article
title_alt Фазовi переходи у геометрично фрустрованi й змiшанiй спiн -1/2 i спiн -1 моделi Iзiнга - Гайзенберга на ромбiчноподiбних декорованих плоских ґратках
description Phase transitions of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on several decorated planar lattices consisting of interconnected diamonds are investigated within the framework of the generalized decoration-iteration transformation. The main attention ispaid to the systematic study of the finite-temperature phase diagrams in dependence on the lattice topology. The critical behaviour of the hybrid quantum-classical Ising-Heisenberg model is compared with the relevant behaviour of its semi-classical Ising analogue. It is shown that both models on diamond-like decorated planar lattices exhibit a striking critical behaviour including reentrant phase transitions. The higher the lattice coordination number is, the more pronounced reentrance may be detected. Дослiджуються фазовi переходи у змiшанiй спiн 1/2 i спiн 1 моделi Iзiнга-Гайзенберга на декiлькох декорованих плоских ґратках, що складаються зi сполучених ромб i в рамках узагальненого декорацiйно-iтерацiйного перетворення. Основна увага придiляється систематичному вивченню фазових дiаграм при скiнчених температурах в залежностi вiд топологiї гратки. Критична поведiнка гiбридної квантово-класичної моделi Iзiнга-Гайзенберга порiвнюється з поведiнкою її напiвкласичного iзінгiвського аналога. Показано, що обидвi моделi на ромбiчноподiбних декорованих плоских ґратках демонструють цiкаву критичну поведiнку, включаючи зворотнi фазовi переходи. Чим вище координацiйне число ґратки, тим помiтнiшою може бути зворотнiсть.
issn 1607-324X
url https://nasplib.isofts.kiev.ua/handle/123456789/119900
citation_txt Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices / Lucia Gálisová, Jozef Strečka // Condensed Matter Physics. — 2011. — Т. 14, № 1. — С. 13002: 1-11. — Бібліогр.: 39 назв. — англ.
work_keys_str_mv AT galisoval phasetransitionsofgeometricallyfrustratedmixedspin12andspin1isingheisenbergmodelondiamondlikedecoratedplanarlattices
AT streckaj phasetransitionsofgeometricallyfrustratedmixedspin12andspin1isingheisenbergmodelondiamondlikedecoratedplanarlattices
AT galisoval fazoviperehodiugeometričnofrustrovaniizmišaniispin12ispin1modeliizingagaizenberganarombičnopodibnihdekorovanihploskihgratkah
AT streckaj fazoviperehodiugeometričnofrustrovaniizmišaniispin12ispin1modeliizingagaizenberganarombičnopodibnihdekorovanihploskihgratkah
first_indexed 2025-11-24T16:25:17Z
last_indexed 2025-11-24T16:25:17Z
_version_ 1850482388872724480
fulltext Condensed Matter Physics 2011, Vol. 14, No 1, 13002: 1–11 DOI:10.5488/CMP.14.13002 http://www.icmp.lviv.ua/journal Phase transitions of geometrically frustrated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices Lucia Gálisová1, Jozef Strečka2 1 Department of Applied Mathematics, Faculty of Mechanical Engineering, Technical University, Letná 9, 04200 Košice, Slovak Republic 2 Department of Theoretical Physics and Astrophysics, Faculty of Science, P. J. Šafárik University, Park Angelinum 9, 04001 Košice, Slovak Republic Received July 26, 2010, in final form October 21, 2010 Phase transitions of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on several decorated planar lat- tices consisting of interconnected diamonds are investigated within the framework of the generalized de- coration-iteration transformation. The main attention is paid to the systematic study of the finite-temperature phase diagrams in dependence on the lattice topology. The critical behaviour of the hybrid quantum-classical Ising-Heisenberg model is compared with the relevant behaviour of its semi-classical Ising analogue. It is shown that both models on diamond-like decorated planar lattices exhibit a striking critical behaviour including reentrant phase transitions. The higher the lattice coordination number is, the more pronounced reentrance may be detected. Key words: Ising-Heisenberg model, decoration-iteration transformation, geometric frustration, reentrant phase transitions PACS: 05.50.+q, 75.10.Hk, 75.10.Jm, 68.35.Rh 1. Introduction Low-dimensional quantum spin models provide an excellent playground for theoretical studies of cooperative and quantum phenomena. It should be pointed out that quantum spin systems, which are prone to a mutual interplay of geometric frustration, thermal and quantum fluctuations, reveal the most remarkable properties [1–3]. Under this specific constraint, one often encounters a strikingly rich physics manifested in diverse unusual ordered and disordered ground states [1], the order-from-disorder effect [1–3], the chirality [4, 5], the enhanced magnetocaloric effect [6–9] or the quantized magnetization plateaus observable at low enough temperatures [2, 10]. It is also well-known dictum that the theoretical investigation of quantum spin models is fre- quently accompanied by rather cumbersome and sophisticated mathematics, which precludes an exact treatment of the most (even simple-minded) spin systems. Hence, it follows that one usually has to rely on an application of some approximative method(s) when treating the vast majority of low-dimensional quantum spin models. However, the approximative methods might have a pro- found deficiency in that they need not be reliable enough in order to decide whether the observed (often subtle) quantum phenomena are of a real physical significance or they merely arose as an artefact of the applied approximation(s). From this point of view, it is highly desirable to search for artificial but exactly tractable quantum spin models, which might display non-trivial quantum effects without a danger of over-interpretation inherent in any approximation. The hybrid Ising-Heisenberg models on diamond-like decorated lattices [11–19] belong to the simplest exactly solved quantum spin models, which were envisaged for describing lattice-statistical systems composed of the semi-classical Ising and the quantum Heisenberg spins. These simplified quantum-classical models can be rigorously treated within the framework of generalized decoration- iteration mapping transformation [20–22], because the nodal Ising spins represent a barrier for c© L. Gálisová, J. Strečka 13002-1 http://dx.doi.org/10.5488/CMP.14.13002 http://www.icmp.lviv.ua/journal L. Gálisová, J. Strečka quantum fluctuations that are consequently restricted to elementary diamond-shaped units only. It is also worth mentioning that the Ising-Heisenberg models on the diamond-like decorated lattices have turned out to be a very useful testing ground for elucidating several typical quantum features. Indeed, these interesting but still exactly tractable spin systems may exhibit diverse quantum ordered and disordered ground states [11–15, 18, 19], the multi-step magnetization process with quantized intermediate magnetization plateaus [14, 17, 18], the enhanced magnetocaloric effect [14, 18], as well as the non-trivial criticality [13, 16, 19]. In our recent works [11, 16, 19], we have examined in detail the ferromagnetic mixed spin- 1/2 and spin-1 Ising-Heisenberg model on several diamond-like decorated planar lattices. It has been demonstrated that this hybrid classical-quantum model exhibits a rather rich ground-state and finite-temperature phase diagrams on account of the competition between the easy-axis Ising and the easy-plane Heisenberg interaction. The main purpose of the present work is to investigate critical properties of the antiferromagnetic mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices, which might even display more complex and non-trivial criticality owing to an interplay of the geometric frustration, thermal and quantum fluctuations as evidenced by our preliminary ground-state analysis [15]. The organization of this paper is as follows. In the next section, we shall define the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated lattices and we shall also briefly review the most important steps of an exact calculation based on the decoration-iteration procedure. The section 3 deals with a discussion of the most interesting numerical results for finite- temperature phase diagrams and temperature variations of the total magnetization both for the classical-quantum Ising-Heisenberg model and its semi-classical Ising analogue. In this section, our main attention will be focused on a possibility of observing reentrant phenomenon. Finally, some concluding remarks are drawn in section 4. 2. Formulation Let us consider two-dimensional lattices composed of inter-connected diamonds as is illustrated in figure 1 for honeycomb, square and triangular lattices. In this figure, the empty circles denote lattice positions of the Ising spins σ = 1/2 that interact with other spins through the interaction JI and the full ones represent lattice positions of the decorating Heisenberg spins S = 1 that interact among themselves via the anisotropic XXZ coupling JH(∆). Note that the parameter ∆ allows to control the interaction JH between the easy-axis (∆ < 1) and easy-plane (∆ > 1) type, as well as, to obtain the Ising model as a special limiting case when assuming ∆ = 0. Providing that the parameter D stands for the axial zero-field splitting (AZFS) parameter [23, 24], which acts on the decorating Heisenberg spins S = 1 only, the total Hamiltonian of the investigated model can be (c)(b)(a) Figure 1. The mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb (figure 1a), square (figure 1b) and triangular (figure 1c) lattices. White circles rep- resent the nodal Ising spins σ = 1/2 and black ones denote the decorating Heisenberg spins S = 1. 13002-2 Phase transitions of geometrically frustrated Ising-Heisenberg model written as Ĥ = −JH Nq/2∑ (i,j) [∆(Ŝx i Ŝ x j + Ŝy i Ŝ y j ) + Ŝz i Ŝ z j ]− JI 2Nq∑ (i,n) Ŝz i σ̂ z n −D Nq∑ i=1 (Ŝz i ) 2. (2.1) Here, the first summation is carried out over all interactions between pairs of the nearest-neigh- bouring Heisenberg spins, the second summation takes into account the interaction between the nearest-neighbouring Ising and Heisenberg spins and the last summation runs over all lattice sites occupied by the decorating Heisenberg spins. The spin variables σ̂z n and Ŝγ i (γ = x, y, z) denote spatial components of the spin-1/2 and spin-1 operators located at nth and ith lattice position, respectively. Finally, N represents the total number of the Ising spins, and q is the coordination number of the lattice, which implies that the total number of all spins (lattice sites) is Ntot = N(q + 1). It is worth mentioning that the quantum-classical model defined in this way is exactly solvable within the framework of a generalized decoration-iteration mapping transformation [20–22] (for more computational details see our recent work [19], where we presented a detailed procedure for an exact solution of this model). As a result of this mapping, one obtains a simple relation between the partition function Z of the investigated mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated lattices and the partition function ZIsing of the simple spin-1/2 Ising model on corresponding undecorated lattices with the nearest-neighbour coupling R: Z(T, JI, JH,∆, D) = ANq/2ZIsing(T,R). (2.2) Note that both mapping parameters A and R are unambiguously determined by a self-consistency condition of the applied decoration-iteration mapping transformation and their explicit forms are given by equations (5)–(7) of reference [19]. At this stage, it is worthwhile to remark that the mapping relation (2.2) is universal and valid regardless of the lattice topology and spatial dimensionality of the model system. In addition, it also permits a rather comprehensive analysis of its critical behaviour as well as some basic thermodynamic quantities. Indeed, it directly follows from equation (2.2) that the investigated mixed-spin Ising-Heisenberg model becomes critical if and only if the spin-1/2 Ising model with the effective coupling R on the corresponding undecorated lattice becomes critical as well. Bearing this in mind, the exact critical temperature of the spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb, square and triangular lattices can be straightforwardly obtained from the conditions βcR = 2 ln(2 + √ 3) , βcR = 2 ln(1 + √ 2) , βcR = ln 3, (2.3) that represent exact relationships between the effective temperature-dependent coupling βR and relevant critical temperatures of the simple spin-1/2 Ising model on a corresponding undecorated honeycomb [25], square [26] and triangular [27] lattices, respectively. In the above, βc = 1/(kBTc) and Tc is the critical temperature of the studied quantum mixed-spin Ising-Heisenberg model. Similarly, the sub-lattice and total magnetization can also be derived from the exact mapping equivalence (2.2) between the partition functions Z and ZIsing. More specifically, by combining equation (2.2) with the exact mapping theorems developed by Barry et al. [28–30] and the gener- alized Callen-Suzuki spin identity [31–33], the sub-lattice magnetization mz i , m z h reduced per one Ising and Heisenberg spin, respectively, can be directly computed from the precise relations: mz i ≡ 〈σ̂z k1〉 = 〈σ̂z k1〉0 ≡ m0, (2.4) mz h ≡ 〈Ŝz k1〉 = 4m0 F (βJI), (2.5) where F (x) = [exp(2βD + βJH) sinh(2x) + exp(βD) sinh(x) cosh(βJH∆)] /W1 , the parameter W1 is defined by equation (6) of reference [19] and the symbols 〈. . .〉 and 〈. . .〉0 represent standard canonical averages performed over the ensemble defined by the mixed-spin Ising-Heisenberg model on the diamond-like decorated lattice and the spin-1/2 Ising model on the corresponding lattice, respectively, andm0 labels the single-site magnetization of the corresponding Ising model. In view of 13002-3 L. Gálisová, J. Strečka this notation, the total magnetization reduced per spin of the investigated lattice, which represents the order parameter of the mixed-spin Ising-Heisenberg model on diamond-like decorated lattices, can be expressed as follows: m = (mz i + qmz h)/(q + 1). 3. Results and discussion In this part, let us proceed to a discussion of the most interesting results obtained for the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb, square and triangular lattices. Before doing this, however, it is worth mentioning that all the results obtained in section 2 are universal since they hold true regardless of whether ferromagnetic or antiferromagnetic interactions JI and JH are assumed, and are independent of the lattice topology or spatial dimensionality of the investigated system. Moreover, the effective nearest-neighbour interaction R of the spin-1/2 Ising model on the corre- sponding undecorated lattice is invariant under the transformation JI → −JI [see equations (5)–(7) in reference [19]]. This observation leads to the conclusion that the critical temperature as well as other thermodynamic quantities of the model under investigation remain unchanged under the change of the nature of the Ising interaction JI. Indeed, a change of the ferromagnetic interaction JI > 0 to the antiferromagnetic one JI < 0 causes just a trivial change in the local alignment of the nodal Ising spins with respect to their nearest Heisenberg neighbours. By contrast, there are some fundamental differences between magnetic behaviour of models with distinct nature of the Heisen- berg interaction JH(∆) (see our preliminary reports [11, 15, 16]). Taking this fact into account, we have restricted our recent work [19] to the particular case with the ferromagnetic Heisenberg and Ising interactions (JH > 0, JI > 0) only. Besides a rather complex ground state composed of two unusual quantum phases, a striking critical behaviour including reentrant phase transitions with two or three consecutive critical points has been discussed in this work. With this background, our attention in this paper will be focused on the analysis of finite-temperature behaviour of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices, in which both the interaction constants JI and JH are supposed to be antiferromagnetic (JI < 0, JH < 0). 3.1. Summary of preliminary results Before proceeding further, let us briefly summarize our previous numerical results published in reference [15]. This work was devoted to the comparison of the ground-state properties of the quantum antiferromagnetic spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices and its semi-classical Ising analogue. The obtained results revealed that the inves- tigated model system has a rich ground-state phase diagram, which consists of a semi-classically ordered ferrimagnetic phase FRI1 with the perfect antiparallel alignment between the nearest- neighbouring Ising and Heisenberg spins and two ferrimagnetic phases FRI2 and FRI3, which differ one from another just in a quantum entanglement of the Heisenberg spin pairs described by the antisymmetric wave function (|1, 0〉 − |0, 1〉)/ √ 2 that emerges just in FRI2. The latter ferrimagnetic phase FRI3 represents the classical ferrimagnetic phase, where the decorating spin pairs may reside in one of the two possible spin states: either |1, 0〉 or |0, 1〉. According to this ambiguous order of the decorating spins, one may regard the phase FRI3 as a partly degenerate state with the non-zero residual entropy S0/NtotkB = [q/(2q + 2)] ln 2, which is proportional to the coordination number q of the lattice (total number of decorating spin pairs). Finally, several geometrically frustrated phases can also be found in the ground state depending on whether the hybrid quantum-classical Ising-Heisenberg model or its semi-classical Ising version is considered. In the former case, the ground state is formed by a frustrated phase FRU, where all nodal Ising spins are frustrated due to a quantum superposition of three spin states |0, 0〉, |1,−1〉 and | − 1, 1〉 of the decorating Heisenberg spins, whose relative probabilities depend on a mutual ratio between interaction parameters. On the other hand, two different frustrated phases FRU1 and FRU2 can be detected for the semi-classical Ising analogue of the model. In the former phase FRU1 , the geomet- ric frustration of the nodal Ising spins is caused by “non-magnetic” nature of the Heisenberg spin 13002-4 Phase transitions of geometrically frustrated Ising-Heisenberg model dimers |0, 0〉, while the geometric frustration in the latter phase FRU2 comes from antiferromag- netic spin states (either |1,−1〉 or | − 1, 1〉) of the Heisenberg spin pairs. As it has been shown in reference [15], both these phases can be regarded as special limiting cases of the unique frustrated phase FRU. 3.2. Finite-temperature behaviour of the semi-classical Ising model Now, let us proceed to the discussion of the finite-temperature behaviour of the antiferromag- netic spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices. To enable a direct comparison with the ground-state analysis published in reference [15], we start first with the discussion of finite-temperature phase diagrams, which are displayed in figure 2 in the re- duced units t = kBT/|JI|, d = D/|JI| and α = |JH|/|JI|, describing the dimensionless temperature, the relative strength of the AZFS parameter and the strength of the Heisenberg interaction normal- ized with respect to the Ising interaction, respectively. This figure shows the critical temperature of the semi-classical Ising version of the investigated mixed-spin model on diamond-like decorated honeycomb (figure 2a), square (figure 2b) and triangular (figure 2c) lattices as a function of the AZFS parameter d for several values of the interaction ratio α. Note that solid lines depicted in this figure are the unique solutions of the critical conditions (2.3) and, as a consequence, they rep- resent the lines of the second-order phase transitions separating the spontaneously ordered phases (FRI1 or FRI3) from the disordered paramagnetic one. As one can clearly see from figure 2, the overall critical behaviour of the system very sensitively depends on the strength of the interaction parameters α and d, as well as, the topology (coordination number) of the lattice; for α < 1, the critical temperature tc either monotonously decreases upon decrease of the AZFS parameter until it tends towards zero temperature at the boundary value d = −1 (see the curves labeled as α = 0.1 and 0.5), or it exhibits an interesting non-monotonous dependence to be closely related to the FRI1 → FRI3 phase transition when the interaction ratio α is sufficiently close to the value α = 1 (see e.g. the case JH/JI = 0.9 and figure 2(a) in reference [15] for clarity). On the other hand, if one considers 1 6 α 6 2, the non-zero critical temperature may be observed just around the interval of the AZFS parameters d ∈ (−1, 1− α), where the partly degenerate ferrimagnetic phase FRI3 constitutes the ground state. Furthermore, several interesting regions with reentrant phase transitions can also be observed in finite-temperature phase diagrams shown in figure 2. More specifically, the diamond-like decorated honeycomb lattice exhibits the reentrant behaviour with either two or three consecutive critical points before the critical line tends to zero temperature when the interaction ratio α acquires very small values, i.e., when the antiferromagnetic interaction JH between decorating spins is much less than the competitive antiferromagnetic interaction JI between the nodal and decorating spins (see the detail of figure 2a). This non-trivial behaviour quickly vanishes as the competitive interaction JH strengthens (see the cases labeled as α = 0.01, 0.05 and 0.1), which suggests that it appears just as a result of a mutual competition between the interaction JI favouring the antiferromagnetic arrangement of nodal Ising spins with respect to their nearest spin neighbours and the AZFS parameter D < 0 that tends to lower the spin states of decorating spins. By contrast, other two diamond-like decorated square and triangular lattices exhibit the reentrant behaviour in the left neighbourhood of the boundary value d = −1, where the frustrated phase FRU1 constitutes the ground state, but for any α ∈ (0, 2) (see figures 2b and c). In these parts of finite-temperature phase diagrams, both the lattices start from the disordered ground state FRU1 before entering the partly ordered phase FRI3 at lower critical temperature tc1, whose spontaneous order disappears due to strong thermal fluctuations at the upper critical temperature tc2. Note that the observed non-trivial behaviour fully corresponds to the condition that the coexistence of the spin order and disorder in the ground state of frustrated Ising systems leads to the occurrence of the reentrant phenomenon at non-zero temperatures [34]. However, this condition is only necessary but not sufficient for the occurrence of the reentrant phenomenon in Ising systems, which can be directly seen from the regions around the boundary points d = 1− α (α > 1) and d = α − 1 (α < 1) in figures 2b, c. Obviously, just the square and triangular lattices show the reentrance in these parts of phase diagrams (see details of figures 2b and c). As the de- tailed ground-state analysis has revealed, the ground-state degeneracy of the ferrimagnetic phase 13002-5 L. Gálisová, J. Strečka Figure 2. The critical temperature of the semi-classical Ising limit of the mixed spin-1/2 and spin- 1 Ising-Heisenberg model on diamond-like decorated honeycomb (figure 2a), square (figure 2b) and triangular (figure 2c) lattices as a function of the AZFS parameter d for several values of the interaction ratio α. The right column illustrates in detail the regions, where the system exhibits reentrant phase transitions with two or three consecutive critical points. 13002-6 Phase transitions of geometrically frustrated Ising-Heisenberg model FRI3 is not high enough in order to cause the reentrant phenomenon in the case of honeycomb lattice (remember that the residual entropy of the honeycomb lattice is S0/NtotkB = 0.2599, while the residual entropy of the square and triangular lattices is S0/NtotkB = 0.2773 and 0.2971, respec- tively). Actually, the square lattice exhibits the reentrant behaviour only if competitive exchange 0.00 0.05 0.10 0.15 0.20 0.25 0.0 0.1 0.2 0.3 = 1.01; d = -0.0130 = 1.02; d = -0.0275 |m | t (a) (b) (c) (d) = 0 q = 4 0.00 0.05 0.10 0.15 0.20 0.25 0.0 0.2 0.4 0.6 = 0.99; d = 0.03250 = 0.98; d = 0.05988 |m | t = 0 q = 4 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.1 0.2 0.3 0.4 0.5 = 1.1; d = 0.00 = 1.5; d = -0.47 = 1.8; d = -0.78 = 0 q = 6 |m | t 0.0 0.1 0.2 0.3 0.4 0.5 0.6 0.0 0.2 0.4 0.6 0.8 1.0 = 0.98; d = 0.50 = 0.97; d = 0.62 |m | t = 0 q = 6 Figure 3. Temperature dependencies of the total magnetization of the decorated square (fig- ures 3a, b) and triangular (figures 3c, d) lattices for such combination of the interaction param- eters α and d, where the system exhibits the reentrant behaviour. interactions between the nearest-neighbouring spins are equal or approximately equal to each other (i.e. if α ≃ 1). On the other hand, the triangular lattice shows the reentrance in this region for the interaction ratio α . 1 as well as for all α ∈ 〈1, 2) (see figure 2c and its detail). It is noteworthy that the region with reentrant phase transitions gradually vanishes as the difference between the considered ratio α and the value 1 increases. These results can also be convincingly evidenced by thermal dependencies of the absolute value of the total magnetization |m| displayed in figure 3 for diamond-like decorated square (figure 3a, b) and triangular (figure 3c, d) lattices. 3.3. Finite-temperature behaviour of the isotropic quantum version of the model To provide a deeper insight into the finite-temperature behaviour of the quantum version of the investigated model system, let us turn our attention to the phase diagrams displayed in fig- ure 4. This figure shows the critical temperature of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb, square and triangular lattices as a function of the interaction ratio α and the AZFS parameter d for the fixed the exchange anisotropy ∆ = 1. Apart from the expected monotonous decrease of tc with the increasing ratio α and/or decreasing d, one may also find here more interesting non-monotonic dependencies of tc as well as several striking regions with reentrant phase transitions. More specifically, the honeycomb lattice shows the reen- trant behaviour with either two or three consecutive critical points just for very small values of the interaction ratio α and for the particular values of the AZFS paremeter d & −1. It is quite obvious that the origin of this non-trivial phenomenon lies in a mutual competition between the antiferromagnetic Ising interaction JI and the easy-plane AZFS parameter D < 0. As could be expected, the parameter region with reentrant phase transitions enlarges with the increasing coor- 13002-7 L. Gálisová, J. Strečka Figure 4. The critical temperature of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb (figures 4a, b), square (figures 4c, d) and triangular (figures 4e, f) lattices the fixed the exchange anisotropy ∆ = 1 as a function of the interaction ratio α for several values of the AZFS parameter d and as a function of the AZFS parameter d for several values of α. 13002-8 Phase transitions of geometrically frustrated Ising-Heisenberg model dination number of the lattice. Indeed, this becomes quite clear from figures 4c and e, the spin-1/2 and spin-1 Ising-Heisenberg model on the decorated square and triangular lattices exhibits two or three critical points for D > 0 and JH(∆), whereas the effect of the former interaction parameter is supported by the Ising interaction JI. Finally, the remarkable parts of finite-temperature phase di- agrams represent particular dependencies of tc forming closed loops (see figure 2 in reference [35]). In these regions of phase diagrams, the system starts with the disordered ground state before en- tering the spontaneously ordered ferrimagnetic phase FRI2 at lower critical temperature tc1, which subsequently disappears due to strong thermal fluctuations at upper critical temperature tc2 (see the curves d = −1.35 in figure 4e and α = 1.15 in figure 4f). The results displayed in figure 4 can 0.0 0.2 0.4 0.6 0.8 1.0 0.0 0.2 0.4 0.6 0.8 1.0 |m | t (a) (b) = 1 = 0.5 q = 6 d = 3.0 0.50.0 -0.2 -0.5 -1.0 -1.2-1.3 -1.35 -0.25(171/2+ 1) 0.0 0.3 0.6 0.9 1.2 0.0 0.2 0.4 0.6 0.8 1.0 |m | t = 1 d = 1.0 q = 6 = 0.1 0.5 0.7 0.78 0.8 0.81 0.83 Figure 5. Temperature dependencies of the total magnetization of the mixed-spin Ising- Heisenberg model on diamond-like decorated triangular lattice with the fixed the exchange anisotropy ∆ = 1 for the interaction ratio α = 0.5 and several values of the AZFS parameter d (figure 5a) and for the AZFS parameter d = 1 and several values of the ratio α. The depen- dencies depicted as broken curves correspond to the two-phase coexistence at t = 0 (see the discussion). be also convincingly evidenced by thermal dependencies of the total magnetization. Since the most diverse finite-temperature behaviour has been observed for the diamond-like decorated triangular lattice, in figure 5 we have depicted some typical temperature variations of the absolute value of total magnetization |m| for this particular model system. As one can clearly see, the magnetization exhibits reentrant transitions with two consecutive critical points for interaction ratio α = 0.5 and the AZFS parameters d . −0.25( √ 17+ 1) as depicted in figure 5a or for the particular case of the AZFS parameter d = 1 and the interaction ratio α & 0.8 as depicted in figure 5b. Note that these values of the interaction parameters α and d are from a close vicinity of the ground-state phase transition between the ferrimagnetically ordered phase (FRI1 or FRI2) and the disordered FRU phase. On the other hand, if one considers the parameter space d > −0.25( √ 17 + 1) for the inter- action ratio α = 0.5 or α < 0.8 for the AZFS parameter d = 1, then the displayed magnetization curves |m(t)| exhibit a single critical point only (see, e.g., the curves d = −0.5, 0.5 in figure 5a and α = 0.1, 0.5, 0.7 in figure 5b). 4. Concluding remarks In the present paper, the critical behaviour of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated honeycomb, square and triangular lattices has been investigated within the framework of the generalized decoration-iteration mapping transformation. Using this rigorous procedure, the exact solution for the investigated mixed-spin model has been obtained by establishing a precise mapping equivalence with the spin-1/2 Ising model on the corresponding undecorated planar lattice with the known exact solution. 13002-9 L. Gálisová, J. Strečka The main emphasis of this work was laid upon the systematic study of the finite-temperature phase diagrams that basically depend on the lattice topology and upon the comparison of the critical behaviour of the hybrid quantum-classical Ising-Heisenberg model and its semi-classical Ising variant. The most interesting result consists in the exact evidence of the existence of reentrant phase transitions with either two or three consecutive critical points. The presence of this non-trivial phenomenon has also been evidenced by temperature dependencies of the total magnetization. It has been proved that the higher the coordination number of the lattice is, the more pronounced reentrance is and the more diverse critical behaviour may be found. Finally, it is worthwhile to remark that even though our theoretical investigation of the mixed spin-1/2 and spin-1 Ising-Heisenberg model on diamond-like decorated planar lattices has been mainly aimed at providing a deeper insight into cooperative and quantum features of this ex- actly solvable model, we hope that our results might stimulate research on possible experimen- tal realizations of this interesting quantum spin model and confirm our theoretical predictions. From this perspective, the most promising approach in experimental realizations of our model system might represent a targeted design of cyano-based polymeric coordination compounds or their isostructural analogues. For instance, the series of bimetallic polymeric coordination com- pounds {Cu(L)}3[Fe(CN)6](ClO4)2 · nH2O, L = N-(3 -aminopropyl)-1, 3 -propanediamine [35, 36] or N-(2 -aminoethyl)-1, 3 -propanediamine [37], which possess the diamond-like decorated honey- comb network structure, may represent a useful starting point for this rational synthesis. In this series, the divalent CuII and FeII metal atoms reside decorating and nodal sites of the diamond-like decorated honeycomb lattice, respectively (see figure 1a). Unfortunately, the divalent FeII atoms are due to a strong ligand field of the cyano group in the diamagnetic low-spin state with S = 0. Similarly, the bimetallic polymeric coordination compounds [W{(CN)4Fe(H2O)2}2] ·nH2O [38] and [W{(CN)4Co(H2O)2}2] · nH2O [39] with the diamond-like decorated square network structure (see figure 1b) have also been reported quite recently, but the tetravalent WIV metal atoms residing nodal sites of the diamond-like decorated square lattice are diamagnetic due to a strong ligand field of the cyano group in those systems. Acknowledgements This work was partially supported by the internal grant of P.J. Šafárik University under the contract No. VVGS 1/10–11 and by European Union European regional development fond (ERDF EU) grant under the contract No. ITMS26220120005 (activity 3.2.). References 1. Lhuillier C., Misguich G., in: High Magnetic Fields: Applications in Condensed Matter Physics and Spectroscopy, Edited by Berthier C., Lévy L.P., Martinez G., Springer-Verlag, Berlin, 2002; doi:10.1007/3-540-45649-X 6 2. Richter J., Schulenburg J., Honecker A., Quantum Magnetism, Lecture Notes in Physics, Vol. 645, Edited by Schollwöck U., Richter J., Farnell D.J.J., Bishop R.F., Springer-Verlag, Berlin, 2004; doi:10.1007/BFb0119592. 3. Misguich G., Lhuillier C., Frustrated Spin Systems, Edited by Diep H.T., World Scientific, Singapore, 2004. 4. Villian J., J. Phys. C: Solid State Phys., 1977, 10, 1717; doi:10.1088/0022-3719/10/10/014. 5. Villain J. et al., J. Physique, 1980, 41, 1263; doi:10.1051/jphys:0198000410110126300. 6. Zhitomirsky M.E., Phys. Rev. B, 2003, 67, 104421; doi:10.1103/PhysRevB.67.104421. 7. Honecker A., Richter J., Condens. Matter Phys., 2005, 8, 813. 8. Derzhko O., Richter J., Eur. Phys. J. B, 2006, 52, 23; doi:10.1140/epjb/e2006-00273-y. 9. Schnack J., Schmidt R., Richter J., Phys. Rev. B, 2007, 76, 054413; doi:10.1103/PhysRevB.76.054413. 10. Honecker A., Schulenburg J., Richter J., J. Phys.: Condens. Matter, 2004, 16, S749; doi:10.1088/0953-8984/16/11/025. 11. Strečka J., Jaščur M., Phys. Status Solidi B, 2002, 233, R12; doi:10.1002/1521-3951(200210)233:3<R12::AID-PSSB999912>3.0.CO;2-2 12. Čanová L., Strečka J., Jaščur M., Czech. J. Phys., 2004, 54, D579. 13002-10 http://dx.doi.org/10.1007/3-540-45649-X_6 http://dx.doi.org/10.1007/BFb0119592 http://dx.doi.org/10.1088/0022-3719/10/10/014 http://dx.doi.org/10.1051/jphys:0198000410110126300 http://dx.doi.org/10.1103/PhysRevB.67.104421 http://dx.doi.org/10.1140/epjb/e2006-00273-y http://dx.doi.org/10.1103/PhysRevB.76.054413 http://dx.doi.org/10.1088/0953-8984/16/11/025 http://dx.doi.org/10.1002/1521-3951(200210)233:3<R12::AID-PSSB999912>3.0.CO;2-2 Phase transitions of geometrically frustrated Ising-Heisenberg model 13. Strečka J., Jaščur M., Acta Phys. Slovaca, 2006, 56, 65. 14. Čanová L., Strečka J., Jaščur M., J. Phys.: Condens. Matter, 2006, 18, 4967; doi:10.1088/0953-8984/18/20/020. 15. Čanová L., Strečka J., Dely J., Jaščur M., Acta Phys. Pol., 2008, 113, 449. 16. Jaščur M., Strečka J., Čanová L., Acta Phys. Pol., 2008, 113, 453. 17. Strečka J., Čanová L., Lučivjanský T., Jaščur M., J. Phys.: Conf. Ser., 2009, 145, 012058; doi:10.1088/1742-6596/145/1/012058. 18. Čanová L., Strečka J., Lučivjanský T., Condens. Matter Phys., 2009, 12, 353. 19. Čanová L., Strečka J., Phys. Status Solidi B, 2010, 247, 433; doi:10.1002/pssb.200945444. 20. Fisher M.E., Phys. Rev., 1959, 113, 969; doi:10.1103/PhysRev.113.969. 21. Syozi I., Phase Transition and Critical Phenomena, Vol. 1, Edited by Domb C., Green M.S., Academic Press, New–York, 1972. 22. Rojas O., Valverde J.S., de Sousa S.M., Physica A, 2009, 388, 1419; doi:10.1016/j.physa.2008.12.063. 23. Rudowicz C., Physica B, 2008, 403, 1882; doi:10.1016/j.physb.2007.10.219. 24. Rudowicz C., Physica B, 2008, 403, 2312; doi:10.1016/j.physb.2007.12.011. 25. Houtappel R.M.F., Physica (Amsterdam), 1950, 16, 425; doi:10.1016/0031-8914(50)90130-3. 26. Onsager L., Phys. Rev., 1944, 65, 117; doi:10.1103/PhysRev.65.117. 27. Baxter R.J., Choy T.C., Proc. R. Soc. London, Ser. A, 1989, 423, 279; doi:10.1098/rspa.1989.0055. 28. Barry J.H., Khatun M., Tanaka T., Phys. Rev. B, 1988, 37, 5193; doi:10.1103/PhysRevB.37.5193. 29. Khatun M., Barry J.H., Tanaka T., Phys. Rev. B, 1990, 42, 4398; doi:10.1103/PhysRevB.42.4398. 30. Barry J.H., Tanaka T., Khatun M., Múnera C.H., Phys. Rev. B, 1981, 44, 2595; doi:10.1103/PhysRevB.44.2595. 31. Callen H.B., Phys. Lett., 1963, 4, 161; doi:10.1016/0031-9163(63)90344-5. 32. Suzuki M., Phys. Lett., 1965, 19, 267; doi:10.1016/0031-9163(65)90978-9. 33. Balcerzak T., J. Magn. Magn. Mater., 2002, 246, 213; doi:10.1016/S0304-8853(02)00056-2. 34. Diep H.T., Giacomini H., Frustrated Spin Systems, Edited by Diep H.T., World Scientific, Singapure, 2004. 35. Hong Ch.S., You Y.S., Inorg. Chim. Acta, 2004, 357, 3271; doi:10.1016/j.ica.2004.04.004. 36. Zhang H.-X., et al., J. Organomet. Chem., 2000, 598, 63; doi:10.1016/S0022-328X(99)00679-8. 37. Trávniček Z., Smékal Z., Escuer A., Marek J., New J. Chem., 2001, 25, 655; doi:10.1039/b006741p. 38. Pilkington M. et al., J. Solid St. Chem., 2001, 159, 262; doi:10.1006/jssc.2001.9155. 39. Herrera J.M. et al., Inorg. Chem., 2003, 42, 7052; doi:10.1021/ic034188+. Фазовi переходи у геометрично фрустрованiй змiшанiй спiн-1/2 i спiн-1 моделi Iзiнга-Гайзенберга на ромбiчноподiбних декорованих плоских ґратках Л. Ґалiсова1, Й. Стречка2 1 Механiко-iнженерний факультет, Технiчний унiверситет, Кошiце, Словацька республiка 2 Природничий факультет, Унiверситет iм. П.Й. Шафарика, Кошiце, Словацька республiка Дослiджуються фазовi переходи у змiшанiй спiн-1/2 i спiн-1 моделi Iзiнга-Гайзенберга на декiлькох декорованих плоских ґратках, що складаються зi сполучених ромбiв в рамках узагальненого декорацiйно-iтерацiйного перетворення. Основна увага придiляється систематичному вивченню фазових дiаграм при скiнчених температурах в залежностi вiд топологiї гратки. Критична поведiнка гiбридної квантово-класичної моделi Iзiнга-Гайзенберга порiвнюється з поведiнкою її напiвкласичного iзiнгiвського аналога. Показано, що обидвi моделi на ромбiчноподiбних декорованих плоских ґратках демонструють цiкаву критичну поведiнку, включаючи зворотнi фазовi переходи. Чим вище координацiйне число ґратки, тим помiтнiшою може бути зворотнiсть. Ключовi слова: модель Iзiгнга-Гайзенберга, декорацiйно-iтерацiйне перетворення, геометрична фрустрацiя, зворотнi фазовi переходи 13002-11 http://dx.doi.org/10.1088/0953-8984/18/20/020 http://dx.doi.org/10.1088/1742-6596/145/1/012058 http://dx.doi.org/10.1002/pssb.200945444 http://dx.doi.org/10.1103/PhysRev.113.969 http://dx.doi.org/10.1016/j.physa.2008.12.063 http://dx.doi.org/10.1016/j.physb.2007.10.219 http://dx.doi.org/10.1016/j.physb.2007.12.011 http://dx.doi.org/10.1016/0031-8914(50)90130-3 http://dx.doi.org/10.1103/PhysRev.65.117 http://dx.doi.org/10.1098/rspa.1989.0055 http://dx.doi.org/10.1103/PhysRevB.37.5193 http://dx.doi.org/10.1103/PhysRevB.42.4398 http://dx.doi.org/10.1103/PhysRevB.44.2595 http://dx.doi.org/10.1016/0031-9163(63)90344-5 http://dx.doi.org/10.1016/0031-9163(65)90978-9 http://dx.doi.org/10.1016/S0304-8853(02)00056-2 http://dx.doi.org/10.1016/j.ica.2004.04.004 http://dx.doi.org/10.1016/S0022-328X(99)00679-8 http://dx.doi.org/10.1039/b006741p http://dx.doi.org/10.1006/jssc.2001.9155 http://dx.doi.org/10.1021/ic034188+ Introduction Formulation Results and discussion Summary of preliminary results Finite-temperature behaviour of the semi-classical Ising model Finite-temperature behaviour of the isotropic quantum version of the model Concluding remarks