Frustrated Spin-Systems









Magnetic interactions are frustrated, if a spin cannot arrange its orientation such that it profits from the interaction with its neighbors. The logo shows the example of Ising spins ( which can only point up or down ) with antiferromagnetic interactions on a triangle. The third spin cannot gain energy, i.e., here spins with antiferro interactions are frustrated.

Frustration can prevent magnetic order; an outstanding example is the Kitaev model on the honeycomb lattice which forms a spin liquid ground state. Here our research focuses on questions as: What is the effect of Heisenberg-type perturbations in this model, or what are the consequences of magnetic fields and defects in the Kitaev-Heisenberg model?

Further research topics include: the Compass-Heisenberg model on the square lattice; moreover we are interested in spin chains formed by copper-oxygen units with frustrated interactions and helical magnetic order.

Another research topic are one-dimensional Wigner lattices and their magnetic interactions. Finally it should be stressed, that there is a close connection to the research on generalized spin-orbital models, which are generically frustrated.

Compass-Heisenberg System / Nanoscale Physics

Frustrated quantum magnetism belongs to the very active research areas in condensed matter theory. In this context we study the compass-Heisenberg model of spins-1/2 on a square lattice, generalizing the quantum compass model via the addition of perturbing Heisenberg interactions between nearest neighbors, and investigate its phase diagram( see Figure 1) and magnetic excitations.

This model has motivations both from the field of strongly correlated systems with orbital degeneracy and from that of solid-state based devices proposed for quantum computing. We find that the high degeneracy of ground states of the compass model is fragile and changes into twofold degenerate ground states for any finite amplitude of Heisenberg coupling. By computing the spin structure factors of finite clusters with Lanczos diagonalization, we evidence a rich variety of phases characterized by Z2 symmetry, that are either ferromagnetic, C-type antiferromagnetic, or of Neel type, and analyze the effects of quantum fluctuations on phase boundaries.

Figure 1: Phase diagram of the compass-Heisenberg model in the plane spanned by the compass interaction Jx/Jz and Heisenberg interaction I/Jz for antiferromagnetic Jz=1. While the quantum compass model (I=0) on a LxL lattice has a 2L degeneracy, this degeneracy is lifted for any finite Heisenberg coupling and a phase with Z2 symmetry is selected.

In the ordered phases the anisotropy of compass interactions leads to a finite excitation gap to spin waves. Furthermore we have shown that for small nanoscale clusters with large anisotropy gap the lowest excitations are column-flip excitations that emerge due to Heisenberg perturbing interactions from the manifold of degenerate ground states of the compass model.

The low energy column-flip or compass-type excitations are of particular interest, as they are robust against decoherence processes and are therefore well designed for storing information in quantum computing. We have also pointed out that the dipolar interactions between nitrogen-vacancy centers forming a rectangular lattice in a diamond matrix may permit a solid-state realization of the anisotropic compass-Heisenberg model.

Kitaev-Heisenberg model on honeycomb lattice

Besides the mathematical beauty of the Kitaev model, recent studies were motivated by its possible relevance for orbitally degenerate systems with strong spin-orbit coupling such as layered iridates Na2IrO3 and Li2IrO3. These applications require to consider the extension, namely the Kitaev-Heisenberg model on honeycomb lattice, which reveals in its phase diagram apart of a liquid phase also several ordered (e.g. stripe or zig-zag) phases.

In particular we have explored the ground state properties of the Kitaev-Heisenberg model in a magnetic field, as well as the evolution of spin correlations in the presence of non-magnetic vacancies. By means of exact diagonalizations, the phase diagram without vacancies was determined as a function of the magnetic field and the ratio between Kitaev and Heisenberg interactions.

Figure 2: Vacancy induced spin density δσjz in the Kitaev-Heisenberg model. The pattern of the z-component of spin around the vacany site (for α=0.9 and an infintesimal magnetic field applied along z) reflects the peculiar spin correlation functions of the model in the Kitaev limit, where only the nearest-neighbor correlations are finite.

In the liquid phase, the magnetization pattern around a single vacancy in a small field is determined (see Figure 2), and its spatial anisotropy is related to that of non-zero further neighbor correlations induced by the field and/or Heisenberg interactions. In the stripe phase, the combination of a vacancy and a small field breaks the six-fold symmetry of the model and stabilizes a particular stripe pattern. Similar symmetry-breaking effects occur even at zero field due to interaction effects between vacancies. This selection mechanism and intrinsic randomness of vacancy positions may lead to spin-glass behavior.

Magnetic frustration in cuprate chain componds

Research of edge-sharing Cu-O chain compounds like LiCuVO4 and NaCu2O2 was stimulated by the discovery of their helical spin structure. The origin of the non-collinear arrangement of spins lies in the frustrated nature of magnetic superexchange interactions J1, J2 and J3 in these compounds. The source of the frustration lies in the relative importance of further neighbor hopping (see Figure 3).

The remarkable magnetoelectric properties are another interesting feature of these compounds. Simultaneously with the helical spin arrangement in LiCuVO4 appears a homogeneous electric polarization field that can be controlled by a magnetic field via the helix! Such a control of electrial polarization by magnetic field or magnetization by an electric field is characteristic for multiferroic compounds.

Figure 3: Effective hopping processes that contribute to superexchange in edgesharing Cu-O chains: In the ground state of the (undoped) chain all Cu ions are in the d9 configuration. By an effective hopping process (which involves intermediate O states) a d9d9 pair is virtually excited into pair formed by a d10 configuration and a d9 ligand-hole pair. The intermediate excitations are singlet configurations and therefore favor antiferromagnetism. Importantly in the edgesharing compounds the 2nd neighbor process, via the Cu-O-O-Cu path, leads to a large AF superexchange integral J2. Due to further processes that involve oxygen J1 is negative. Typically these interactions lead to frustration and a helical spin arrangement.

Excitons in the Mott-Hubbard gap: Edgesharing chains have very small hopping matrix elements as compared to the effective Hubbard U and the longer range Coulomb interactions (V1,V2,V3....). We have shown that in this extreme narrow band case excitons in the Hubbard gap are expected at energies U-V1 and U-V2. This exciton doublet has been seen by ellipsometry experiments in LiCuVO4 and NaCu2O2.

Importantly the intensity of the absorptions is determined by the nearest and next-nearest neighbor spin correlation functions. Thus the optical absorption can be used to study the strength and the temperature dependence of the spin-correlations of the frustrated spin chains.

One-dimensional Wigner lattice

Motivated by the finding that the intrinsically doped edge-sharing Cu-O chain compounds Na3Cu2O4 and Na8Cu5O10, that have been synthesized in Jansen's department, are realizations of one-dimensional (1D) Wigner crystals, we have investigated the electronic properties of these systems. Due to the 90 degree Cu-O-Cu bonds in this structure the kinetic energy is small and the electrons are strongly correlated and their positions controlled by the long-range Coulomb interaction. In particular we studied the optical and the density fluctuation spectra, the single particle spectral functions, and the phase diagram of such systems. Charge excitations in 1D Wigner crystals are described in terms of domain-wall excitations with fractional charge, rather than by particle-hole excitations as in conventional insulators.

A very interesting aspect of the edge-sharing compounds is the relevance of both nearest- and next-nearest neighbor hopping, t1 and t2, which allows to distinguish charge modulations due to Fermi surface instabilities and due to Coulomb interactions, i.e., Wigner lattice formation. In particular we have shown that a soft domain-wall exciton at sufficiently large t2 drives an instability of the Wigner lattice towards a charge-density wave with a modulation period distinct from that of the Wigner lattice.

Figure 4: (a) Schematic representation of Wigner charge order in the edge-sharing CuO chain compound Na5Cu3O6. In this case the unit cell involves 3 Cu ions, where filled circles represent Cu d9 configurations, i.e., sites with spin 1/2, while the open circle corresponds to a non-magnetic Cu site with a hole (Zhang-Rice singlet). The relevant magnetic superexchange interactions J1, J2 and J3 are indicated; (b) and (c) are effective spin models where the sites with holes are eliminated. Such effective spin models, also augmented by relevant interchain couplings, are being used to calculate the magnetic properties.

More recently we discovered that virtual excitations across the Wigner charge gap ΔW contribute to the magnetic superexchange. For example J2 is proportional to t12t2W2 due to constructive interference in the singlet channel. As these terms are odd in t2 they favor effectively the ferromagnetic state for one sign of t2. These exchange interactions are controlled by the Wigner charge gap ΔW and not by the much larger Mott-Hubbard gap ΔMH as the usual AF interactions in the system, thus the new terms become more relevant in systems with decreasing ΔW. Recently, in Jansen's department another doped chain compound Na5Cu3O6 was synthesized whose susceptibility could only be explained by invoking the above mechanism that amplifies the tendency towards ferromagnetism (see Figure 4).

Local experts:

External collaborators:

  • Maria Daghofer (Dresden)

  • R. M. Noack (Marburg)

  • A. M. Oles (Krakau)

  • Jesko Sirker (Kaiserslautern)

Key publications:

Most recent publications of department Metzner

Last modified: 13. August 2013

Go to Editor View