Outline
Comptes Rendus

Account / Revue
Frustration-induced exotic properties of magnetic molecules
Comptes Rendus. Chimie, Corrélations de spin et de charge dans les matériaux moléculaires - Propriétés physiques, chimie et aspects « matériaux », Volume 10 (2007) no. 1-2, pp. 15-20.

Abstract

Geometric frustration of interacting spin systems is the driving force of a variety of fascinating phenomena in low-dimensional magnetism. In this contribution, I will review recent results on frustration-induced effects in magnetic molecules, i.e. zero-dimensional magnetic systems, as well as in a recently synthesized frustrated molecule-based spin tube, i.e. a one-dimensional spin system.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2006.07.010
Keywords: Magnetic molecules, Heisenberg model, Antiferromagnetism, Three-leg ladder

Jürgen Schnack 1

1 Universität Osnabrück, Fachbereich Physik, Barbarastr. 7, 49069 Osnabrück, Germany
@article{CRCHIM_2007__10_1-2_15_0,
     author = {J\"urgen Schnack},
     title = {Frustration-induced exotic properties of magnetic molecules},
     journal = {Comptes Rendus. Chimie},
     pages = {15--20},
     publisher = {Elsevier},
     volume = {10},
     number = {1-2},
     year = {2007},
     doi = {10.1016/j.crci.2006.07.010},
     language = {en},
}
TY  - JOUR
AU  - Jürgen Schnack
TI  - Frustration-induced exotic properties of magnetic molecules
JO  - Comptes Rendus. Chimie
PY  - 2007
SP  - 15
EP  - 20
VL  - 10
IS  - 1-2
PB  - Elsevier
DO  - 10.1016/j.crci.2006.07.010
LA  - en
ID  - CRCHIM_2007__10_1-2_15_0
ER  - 
%0 Journal Article
%A Jürgen Schnack
%T Frustration-induced exotic properties of magnetic molecules
%J Comptes Rendus. Chimie
%D 2007
%P 15-20
%V 10
%N 1-2
%I Elsevier
%R 10.1016/j.crci.2006.07.010
%G en
%F CRCHIM_2007__10_1-2_15_0
Jürgen Schnack. Frustration-induced exotic properties of magnetic molecules. Comptes Rendus. Chimie, Corrélations de spin et de charge dans les matériaux moléculaires - Propriétés physiques, chimie et aspects « matériaux », Volume 10 (2007) no. 1-2, pp. 15-20. doi : 10.1016/j.crci.2006.07.010. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2006.07.010/

Version originale du texte intégral

1 Introduction

Antiferromagnets with competing interactions are said to be geometrically frustrated. As a result, they exhibit a variety of fascinating phenomena [1]. In this context, the term frustration describes a situation where, in the ground state of the corresponding classical spin system, not all interactions can be saturated simultaneously. A typical picture for such a situation is a triangle of antiferromagnetically coupled spins, where classically the spins are not in the typical up–down–up configuration, but assume a ground state that is characterized by a relative angle of 120° between neighboring spins. This special type of classical ground state characterizes several frustrated spin systems, among them are giant Keplerate molecules [2–4], the triangular lattice antiferromagnet, and the kagome lattice antiferromagnet.

Among the phenomena caused by frustration are magnetization plateaus and jumps as well as unusual susceptibility minima, as observed for example for the kagome lattice antiferromagnet [5,6]. Some of these effects can also occur in certain strongly frustrated magnetic molecules such as the Keplerate {Mo72Fe30} [7] and in a cuboctahedral molecule [8] which has been already synthesized as {Cu12La8}, but has not yet been magnetically characterized [9]. It has also soon be noticed that the high degeneracy of magnetic levels at the saturation field, where the unusual jumps occur, should also have an enormous impact on the magnetocaloric effect. An adiabatic process should yield a massively enhanced cooling rate in the vicinity of large magnetization jumps. Such a behavior was first anticipated for the classical version of the kagome lattice antiferromagnet [10] and subsequently demonstrated for certain low-dimensional quantum spin systems [11]. In addition it could be shown that this effect is rather stable against perturbations [12]. In the following sections, some of these effects will be reviewed in more detail.

Interlinking the zero-dimensional molecular units to cluster-based networks is a highly interesting new route in molecular magnetism, especially if the degree of magnetic exchange between the cluster entities can be varied over a wide range. Examples for such cluster-based networks are given by chains [13] and square lattices [14,15] of the magnetic Keplerate molecule {Mo72Fe30}. These systems show new combinations of physical properties that stem from both molecular and bulk effects. As an example for a quasi one-dimensional molecule-based frustrated spin system, the recently synthesized three-leg ladder compound [(CuCl2tachH)3Cl]Cl2 will be discussed [16,17].

Throughout the article the spin systems are modeled by an isotropic Heisenberg Hamiltonian augmented with a Zeeman term, i.e.:

H˜=u,vJuvs˜(u)·s˜(v)+gμBBS˜z.(1)
s˜(u) is the individual spin operator at sites u, S˜ is the total spin operator, and S˜z its z-component along the homogeneous magnetic field axis. Juv is the matrix element of the symmetric coupling matrix. In the following, we will consider only antiferromagnetic couplings which are characterized by negative values of Juv.

2 Frustration effects in magnetic molecules

For not too strong frustration, i.e. not too large deviation from bipartiteness [18,19], many spin systems possess properties that are similar to or in accord with their unfrustrated, i.e. bipartite counterparts. This is for instance the case for antiferromagnetically coupled Heisenberg spin rings, which – if even-membered – are covered by such famous rigorous results as the sign rule of Marshall and Peierls [20] and the theorems of Lieb, Schultz, and Mattis [18,19]. A key quantity of interest is the shift quantum number k=0,,N1 associated with the cyclic shift symmetry of the rings. The corresponding crystal momentum is then 2πk/N. For the aforementioned rings with even N (bipartite) one can explain the shift quantum numbers for the relative ground states in subspaces H(M) of total magnetic quantum number M [18–20]. In recent investigations we could numerically verify and prove in some cases that even for frustrated rings with odd N astonishing regularities hold, and thus a generalized sign rule can be established, which holds for chains of arbitrary sizes [21].

An antiferromagnet that can be decomposed into two sublattices has as its lowest excitations the rotation of the Néel vector as well as spin wave excitations [22]. In finite size systems, these excitations are arranged in rotational (parabolic) bands of energy eigenvalues as a function of total spin. Such a behavior is most pronounced for bipartite, i.e. unfrustrated systems [23–25]. It turns out that the idea of rotational bands can be used to approximate the low-lying energy spectrum of the Keplerate molecule {Mo72Fe30} which, due to the huge size of the Hilbert space, practically, cannot be dealt with by other quantum methods [26]. Such an approximation can explain the field dependence of the magnetization at low temperatures [26]. It also predicts a strong resonance in the inelastic neutron scattering cross-section at the separation between the two lowest bands, which was indeed discovered in a recent experiment using inelastic neutron scattering [27].

Nevertheless, the agreement of some properties with those of the corresponding non-frustrated counterparts is only approximate. Frustrated antiferromagnets exhibit their very own characteristic features. One of them is given by independent magnons which occur in certain spin arrays. It has been found that the minimal energies Emin(S) of antiferromagnetic molecules of cuboctahedral and icosidodecahedral structures depend linearly on total spin S above a certain total spin [28]. Such a dependence, which is depicted in Fig. 1 (top) for an icosidodecahedron with intrinsic spin s = 1/2, results in an unusually big jump to saturation, as can be seen in Fig. 1 (bottom). Although first noticed for the Keplerate molecule {Mo72Fe30}, such a behavior is quite common for a whole class of frustrated spin systems, including for instance the kagome or the pyrochlore lattice antiferromagnet [6,29], as well as for other frustrated magnetic molecules [8]. The underlying reason is that, due to the special geometric frustration in such systems – even polygons are surrounded by triangles –, the relative ground states in subspaces H(M) are for big enough M given by product states of independent localized magnons. This is one of the rare examples where eigenstates of an interacting many-body Hamiltonian can be constructed analytically.

Fig. 1

Top: minimal energies of the icosidodecahedron for s = 1/2. The highest four levels fall on a straight line. Bottom: resulting (T = 0) − magnetization curve. The magnetization jump of ΔM = 3 is marked by an arrow.

From a condensed matter point of view, the occurrence of independent and localized magnons can be understood by looking at the band structure in one-magnon space. Fig. 2 shows as an example the energy levels of the sawtooth chain, a structure that can exhibit a truly giant magnetization jump of half of the saturation magnetization [29]. In one-magnon space, the energy levels form two bands: an N/2-times degenerate (flat) ground state band and a cosine shaped excited band. It is clear that the degenerate states of the flat band can be superimposed in order to form localized objects which are independent if well-enough separated. The occurrence of flat bands seems to be rather common [30]: they exist for instance in the magnetic molecules of icosidodecahedral and cuboctahedral structures discussed previously as well as in the kagome and the pyrochlore antiferromagnet [29].

Fig. 2

Top: structure of the sawtooth chain, which hosts independent magnons (V) for J2 = 2J1. Bottom: energy bands in one-magnon space; the flat band consists of N/2 degenerate levels.

As the last point of this section, I want to discuss the observation that hysteresis can occur without anisotropy. Usually the observation of hysteresis effects in magnetic materials is an outcome of their magnetic anisotropy. In a recent article, we could report that magnetic hysteresis (Fig. 3) occurs in a spin system without any anisotropy [31]. Specifically, we investigated an icosahedron where classical spins mounted on the vertices are coupled by antiferromagnetic isotropic nearest-neighbor Heisenberg interaction, giving rise to geometric frustration. At T = 0 this system undergoes a first-order metamagnetic phase transition at a critical field between two distinct families of ground-state configurations. The metastable phases of the system are characterized by temperature and field-dependent survival probability distributions. Our exact classical treatment shows that the abrupt transition at T = 0 originates in the intersection of two energy curves belonging to different families of spin configurations that are ground states below and above the critical field. The minimum of the two energy functions constitutes a non-convex minimal energy function of the spin system and this gives rise to a metamagnetic phase transition [32]. We could also show that the corresponding quantum spin system for sufficiently large spin quantum number s ≥ 4 possesses a non-convex set of lowest energy levels when plotted versus total spin [31], which again leads to a nontrivial magnetization jump. Nevertheless, the investigation of the hysteretic behavior of the quantum systems remains an open problem so far, since the magnetization dynamics should involve the modeling of the relaxation of the quantum spin system in a time-dependent field and coupled to a thermostat. This future investigation will either be performed by solving the complete Schrödinger dynamics of the coupled spin and bath systems or by following effective Master equations (compare Ref. [33]).

Fig. 3

Hysteresis behavior of the classical icosahedron in an applied magnetic field obtained by classical spin dynamics simulations (thick lines) as well as by analytical stability analysis (grey lines). The curves match perfectly [31].

3 Frustration effects of the three-leg ladder compound [(CuCl2tachH)3Cl]Cl2

Recently we reported the magnetic features of a new one-dimensional stack of antiferromagnetically coupled equilateral copper(II) triangles, see top of Fig. 4 [16,17]. High-field magnetization measurements (bottom of Fig. 4) showed that the interaction between the copper triangles is of the same order of magnitude as the intra-triangle exchange, although only coupled via hydrogen bonds. Thus, this molecule-based infinite chain turns out to be an interesting example of a frustrated cylindrical three-leg ladder with competing intra- and inter-triangle interaction.

Fig. 4

Top: schematic structure of the triangular copper chain. The copper ions are placed at the vertices, intra-triangle exchange pathways are drawn by grey lines (J1 = −0.9 K), inter-triangle couplings are given by black lines (J2 = −1.95 K). Bottom: magnetization per copper triangle. Experimental data taken at a cryostat temperature of 0.4 K are given by thick solid lines, the two datasets are practically identical. The theoretical estimates are given for 18 and 24 spins, i.e. 6 or 8 triangles. The saturation field is Bsat = 11 T [17].

A key question concerns the structure of the low-lying levels, and especially the existence of a singlet–triplet and a singlet–singlet gap. Finite-size extrapolations of Lanczos diagonalization results seemed to suggest that the singlet ground state is gapped from the triplet excitation (compare top panel of Fig. 5). For all finite size systems, it turned out that also the first excited singlet is gapped from the ground-state singlet (see bottom panel of Fig. 5). It is amazing that the singlet–singlet gap seems to be nearly independent of the system's size.

Fig. 5

Top: finite-size extrapolation of the singlet–triplet gap Δ1 = Emin(S = 1) − Emin(S = 0) and the pentuplet–single gap Δ2 = Emin(S = 2) − Emin(S = 0) for [(CuCl2tachH)3Cl]Cl2. Both gaps for the infinite chain are estimated to be non-zero. Bottom: finite-size behavior of the singlet–singlet gap for [(CuCl2tachH)3Cl]Cl2. The two leftmost data points for 8 and 10 triangles have an uncertainty of about ± 0.2 K due to the restricted accuracy of the Lanczos procedure for excited states [17]. The dashed line in both figures shows the upper boundary of the gap size according to Ref. [41].

Our findings provoked theoretical investigations of related models [34–40], which yield that this frustrated spin system should not have a gap. The arguments are as follows:

  • (i) It is shown in Ref. [34] that the Hamiltonian of a ladder compound such as [(CuCl2tachH)3Cl]Cl2 can be mapped onto a simpler Hamiltonian of a linear chain of effective spins. For weak inter-triangle couplings J2 each triangle basically remains in the triangular ground state with Striangle = 1/2, and thus the ladder should be an effective s = 1/2 chain, which is augmented by additional chiral degrees of freedom [34]. Such a chain should have a twofold degenerate ground state. In the opposite case where the inter-triangle coupling dominates, the triangular ground state should have Striangle = 3/2, and then the ladder should be an effective s = 3/2 chain, which, since being a half-integer chain, is gapless.
  • (ii) Ref. [35] also provides a numerical argument. The Hamiltonian of [(CuCl2tachH)3Cl]Cl2 is diagonalized for the relevant exchange parameters J1 = −0.9 K, J2 = −1.95 K with the help of Density Matrix Renormalization Group Techniques (DMRG). Using DMRG one can determine the ground state energies in subspaces H(M) of total magnetic quantum number M for much larger chains. Nevertheless, due to technical problems with the boundary conditions which are intrinsic to the method only the gap to the pentuplet Emin(S = 2) − Emin(S = 0) could be determined, which indeed tends to zero.
  • (iii) The most general argument of Ref. [35] states that the three-leg ladder should either have a degenerate ground state or be gapless according to a proof given in Ref. [41]. This proof is applicable for arbitrary ladder compounds with an odd number of sites on the rungs and half-integer intrinsic spin quantum number as long as the system possesses translational symmetry along the ladder. To the best of our knowledge this proof should be correct.

The above arguments suggest that the finite size extrapolation of Ref. [17] might be biased by too short sample systems (up to 10 triangles). In addition the role of the boundary conditions might be crucial for the actual numerical convergence. In order to estimate, to which system size a precise gap estimation must be performed in order to definitely decide about low-lying gaps, the upper bound according to Ref. [41] is displayed as a dashed line in Fig. 5. It seems realistic to reach conclusion about the singlet–singlet gap, whereas the two other gaps Δ1 = Emin(S = 1) − Emin(S = 0) and Δ2 = Emin(S = 2) − Emin(S = 0) approach the theoretical bound for rather large system sizes only.

Summarizing, the existence of low-lying gaps in the new frustrated three-lag ladder appears to be ruled out for the perfect model system. Investigations on the real spin system, both theoretically and experimentally, will continue.

Acknowledgement

It is my pleasure to thank K. Bärwinkel, A. Honecker, P. Kögerler, M. Luban, H. Nojiri, J. Richter, H.-J. Schmidt, and C. Schröder for the fruitful collaboration that produced so many exciting results. I am especially thankful to H.-J. Schmidt for valuable discussions about the three-lag ladder compound. This work was supported by the Ph.D. program of the University of Osnabrück as well as by the Deutsche Forschungsgemeinschaft (grant nos. SCHN 615/5-1, SCHN 615/5-2, and SCHN 615/8-1).


References

[1] J. Greedan J. Mater. Chem., 11 (2001), p. 37

[2] A. Müller; M. Luban; C. Schröder; R. Modler; P. Kögerler; M. Axenovich; J. Schnack; P.C. Canfield; S. Bud'ko; N. Harrison Chem. Phys. Chem., 2 (2001), p. 517

[3] A. Müller; A.M. Todea; J. van Slageren; M. Dressel; H. Bögge; M. Schmidtmann; M. Luban; L. Engelhardt; M. Rusu Angew. Chem., Int. Ed. Engl., 44 (2005), p. 3857

[4] B. Botar; P. Kögerler; C.L. Hill Chem. Commun. (2005), p. 3138

[5] Y. Narumi; K. Katsumata; Z. Honda; J.-C. Domenge; P. Sindzingre; C. Lhuillier; Y. Shimaoka; T.C. Kobayashi; K. Kindo Europhys. Lett., 65 (2004), p. 705

[6] J. Schulenburg; A. Honecker; J. Schnack; J. Richter; H.-J. Schmidt Phys. Rev. Lett., 88 (2002), p. 167207

[7] C. Schröder; H. Nojiri; J. Schnack; P. Hage; M. Luban; P. Kögerler Phys. Rev. Lett., 94 (2005), p. 017205

[8] R. Schmidt; J. Schnack; J. Richter J. Magn. Magn. Mater., 295 (2005), p. 164

[9] A.J. Blake; R.O. Gould; C.M. Grant; P.E.Y. Milne; S. Parsons; R.E.P. Winpenny J. Chem. Soc., Dalton Trans. (1997), p. 485

[10] M.E. Zhitomirsky Phys. Rev., B67 (2003), p. 104421

[11] M.E. Zhitomirsky; A. Honecker J. Stat. Mech. (2004), p. P07012

[12] O. Derzhko; J. Richter Phys. Rev., B70 (2004), p. 104415

[13] A. Müller; S.K. Das; M.O. Talismanova; H. Bögge; P. Kögerler; M. Schmidtmann; S.S. Talismanov; M. Luban; E. Krickemeyer Angew Chem., Int. Ed. Engl., 41 (2000), p. 579

[14] A. Müller; E. Krickemeyer; S.K. Das; P. Kögerler; S. Sarkar; H. Bögge; M. Schmidtmann; S. Sarkar Angew. Chem., Int. Ed. Engl., 39 (2000), p. 1612

[15] A. Müller; S.K. Das; E. Krickemeyer; P. Kögerler; H. Bögge; M. Schmidtmann Solid-State Sci., 2 (2000), p. 847

[16] G. Seeber; P. Kögerler; B.M. Kariuki; L. Cronin Chem. Commun. (2004), p. 1580

[17] J. Schnack; H. Nojiri; P. Kögerler; G.J.T. Cooper; L. Cronin Phys. Rev., B70 (2004), p. 174420

[18] E.H. Lieb; T. Schultz; D.C. Mattis Ann. Phys. (N.Y.), 16 (1961), p. 407

[19] E.H. Lieb; D.C. Mattis J. Math. Phys., 3 (1962), p. 749

[20] W. Marshall Proc. R. Soc. A (Lond.), 232 (1955), p. 48

[21] K. Bärwinkel; P. Hage; H.-J. Schmidt; J. Schnack Phys. Rev., B68 (2003), p. 054422

[22] P.W. Anderson Phys. Rev., 86 (1952), p. 694

[23] J. Schnack; M. Luban Phys. Rev., B63 (2001), p. 014418

[24] O. Waldmann Phys. Rev., B65 (2002), p. 024424

[25] O. Waldmann; T. Guidi; S. Carretta; C. Mondelli; A.L. Dearden Phys. Rev. Lett., 91 (2003), p. 237202

[26] J. Schnack; M. Luban; R. Modler Europhys. Lett., 56 (2001), p. 863

[27] V.O. Garlea; S.E. Nagler; J.L. Zarestky; C. Stassis; D. Vaknin; P. Kögerler; D.F. McMorrow; C. Niedermayer; D.A. Tennant; B. Lake; Y. Qiu; M. Exler; J. Schnack; M. Luban Phys. Rev., B73 (2006), p. 024414

[28] J. Schnack; H.-J. Schmidt; J. Richter; J. Schulenburg Eur. Phys. J., B24 (2001), p. 475

[29] J. Richter; J. Schulenburg; A. Honecker; J. Schnack; H.-J. Schmidt J. Phys.: Condens. Matter, 16 (2004), p. S779

[30] S. Miyahara; K. Kubo; H. Ono; Y. Shimomura; N. Furukawa J. Phys. Soc. Jpn, 74 (2005), p. 1918

[31] C. Schröder; H.-J. Schmidt; J. Schnack; M. Luban Phys. Rev. Lett., 94 (2005), p. 207203

[32] C. Lhuillier; G. Misguich High Magnetic Fields, Springer, Berlin, 2002 (p. 161)

[33] M.J. Henrich; M. Michel; M. Hartmann; G. Mahler; J. Gemmer Phys. Rev., E72 (2005), p. 026104

[34] A. Lüscher; R.M. Noack; G. Misguich; V.N. Kotov; F. Mila Phys. Rev., B70 (2004), p. 060405 (R)

[35] J.-B. Fouet; A. Läuchli; S. Pilgram; R.M. Noack; F. Mila Phys. Rev., B73 (2006), p. 014409 (2005)

[36] K. Okunishi; S. Yoshikawa; T. Sakai; S. Miyashita Prog. Theor. Phys. Suppl., 159 (2005), p. 297

[37] T. Sakai; M. Matsumoto; K. Okunishi; K. Okamoto; M. Sato Physica E, 29 (2005), p. 633

[38] M. Matsumoto; T. Sakai; M. Sato; H. Takayama; S. Todo Physica E, 29 (2005), p. 660

[39] M. Sato Phys. Rev., B72 (2005), p. 104438

[40] M. Sato J. Phys. Chem. Solids, 66 (2005), p. 1454

[41] A.G. Rojo Phys. Rev., B53 (1996), p. 9172


Comments - Policy