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...
Gespeichert in:
| Veröffentlicht in: | Condensed Matter Physics |
|---|---|
| Datum: | 2011 |
| Hauptverfasser: | , |
| 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
|