## Abstract

Electron paramagnetic resonance (EPR) is a powerful spectroscopic technique, perfectly suited for determining magnetic anisotropy terms in spin Hamiltonians. Although solid foundations of the EPR theory were established by Kubo and Tomita (KT) more than half a century ago, especially in the last couple of decades, we have witnessed a rapid progress in the field due to the occurrence of enhanced computational capabilities. In this chapter, we overview this progress by summarizing the basic concepts of EPR in exchange-coupled systems. The review builds upon the standard KT theory and the exchange narrowing picture, which is however only suitable at high enough temperatures and for systems with dimensionality exceeding one. We also summarize the predictions of more modern approaches, including exact calculations on finite spin clusters, the Oshikawa-Affleck effective-field theory for 1D systems, and the recently developed EPR-moments approach. Many illuminating examples of the applicability of different approaches are also provided.

### Keywords

- EPR
- ESR
- EMR
- Kubo-Tomita theory
- exact diagonalization
- Oshikawa-Affleck effective-field theory
- EPR moments
- exchange-coupled spin systems
- magnetic anisotropy
- Dzyaloshinskii-Moriya interaction
- anisotropic exchange
- single-ion anisotropy

## 1. Introduction

Since the pioneering demonstration of the electron paramagnetic resonance phenomenon in solids and liquids in 1944 by Zavoisky [1, 2, 3], EPR has become a well-established and broadly spread spectroscopic technique. Although the main principle of detecting microwave absorption by electronic magnetic moments at a fixed frequency and a sweeping applied magnetic fields has not changes from early days, the method has become one of the most sensitive local probes of magnetism.

The term electron paramagnetic resonance in a narrow sense applies to paramagnetic compounds containing transition-metal or rare-earth elements with incomplete inner shells, hence possessing paramagnetic electron moments. In a broader sense, the general term electron magnetic resonance (EMR) stands for magnetic resonance absorption experiments performed on an ensemble of magnetic moments corresponding to localized or itinerant electrons. In addition to paramagnets, EMR thus also covers absorption phenomena in ordinary metals and magnetically ordered systems, as well as absorption by imperfections in insulators and semiconductors, which may trap electrons or holes. In literature, also the term electron spin resonance (ESR) is often encountered. This is usually reserved for cases when the magnetic moment originates primarily from the spin momentum of the electron, like in iron-group metals where the orbital moment of the electron is usually quenched.

The EPR technique provides superior insight into magnetic properties of a particular sample compared to more conventional bulk magnetic techniques, e.g., bulk-magnetization or magnetic-torque measurements. A particular EPR experiment can provide information that help in characterization of local magnetic and electrostatic environments of a magnetic moment, as well as information about development of magnetic correlations and fluctuations [4, 5, 6]. The experiment can also help to determine magnetic coupling with other electronic and nuclear moments, etc. Due to these diverse and detailed information, EPR has earned reputation in various fields of science. Traditionally, it was in the domain of solid-state physics and chemistry, but lately it has become indispensable also in bio-oriented sciences and medical applications. Moreover, it has been recently highlighted for its strength in detecting unconventional magnetic phenomena, such as edge states in topological insulators [7], spinon excitations in spin liquids [8], and spin-nematic states [9].

For a general introduction to EPR the reader is advised to turn to one of many very good EPR monographs and reviews, like the Abragam and Bleaney “Electron Paramagnetic Resonance of Transition Ions” [4], the Pilbrow “Transition Ion Paramagnetic Resonance” [5], or the more recent Weil and Bolton “Electron Paramagnetic Resonance: Elementary Theory and Practical Applications” [6]. The purpose of this chapter is to review a specific problem of EPR in exchange coupled systems. This problem is particularly difficult to treat due to complications induced by the exchange interaction between neighboring moments. These interactions dramatically affect the way the moments respond to the external magnetic fields. In order to model this response properly, the use of modern theoretical concepts and advanced experimental approaches is required. These are review in this chapter.

The outline of the chapter is the following. We will start with a general overview of the Kubo-Tomita EPR theory (Section 2), which will first require the introduction of the spin-Hamiltonian concept. We will pay special attention to the exchange-narrowing limit, which is generally applicable to strongly-exchange-coupled spin systems. Next, a few successful applications of the KT theory will be demonstrated in Section 3. In Section 4, limitations of the KT approach will be summarized. Different approaches that can overcome these limitations and their specific applications will also be given. The concluding Section 5 will summarize this chapter.

## 2. KT theory of EPR in exchange-coupled systems

Dense magnetic insulators, i.e., systems that do not conduct electric current and where magnetic moments are localized at well-defined crystallographic sites (usually occupied by transition metals or rare earths), represent one of the major fields of research in condensed matter, where EPR is particularly powerful [10]. In this chapter, we shall focus on systems that are strongly exchange coupled, i.e., where magnetic moments communicate, and highlight particular information that EPR can provide in such cases.

EPR measures the absorption of microwaves by electrons, i.e., atomic magnetic moments, therefore, it provides a direct insight to the atomic magnetism. This is unlike some other local-probe techniques, such as nuclear magnetic resonance [11, 12] or muon spectroscopy [13] that can only provide indirect information about electron degrees of freedom. However, as we shall see below, this advantage of EPR at the same time turns out to be a drawback, since knowledge of four-spin correlations functions is required to accurately describe the EPR response of exchange-coupled magnetic moments at an arbitrary temperature. On the other hand, for indirect techniques, like NMR, two-spin correction functions suffice. This makes EPR an elaborate technique and prevents a routine analysis of the EPR spectra of exchange-coupled systems.

The beginnings of the EPR theory in exchange-couples magnetic systems go back to the seminal work by Kubo and Tomita (KT) entitled “General Theory of Magnetic Resonance Absorption” [14]. Although it rests on a perturbation approach and is therefore not exact, the KT theory still represents solid foundations in modern times. The EPR theory has seen some progress later on, especially in recent years with the advent of enhanced computational facilities. Within this chapter, we shall make a general overview of the KT theory and its successors that were developed for cases where the KT theory is not valid.

### 2.1. Spin Hamiltonian

We start the body of this review with introducing the concept of the spin Hamiltonian. In this framework the total Hamiltonian of a particular system with all degrees of freedom that are present, i.e., electron orbital, electron spin, nuclear, lattice, etc., is projected onto the spin space of the electrons. In an external magnetic field, the spin Hamiltonian comprises of the following terms [6]:

Here, *g* tensor),

is the exchange Hamiltonian summing terms between electron spins at sites *i* and *j* coupled by an isotropic exchange interaction

includes the dipolar term between electronic spins

where

known as the Dzyaloshinskii-Moriya (DM) interaction (

The exchange anisotropy is a relativistic effect due to the spin-orbit coupling. In transition metals, the Dzyaloshinskii-Moriya interaction is usually the dominant exchange anisotropy term. The reason is that it originates from the first-order perturbation in the spin-orbit coupling, while the symmetric anisotropic exchange results only from the second-order perturbation theory [15, 16]. Consequently, the DM term is proportional to *g-*shift Δ*g* from the free electron value of 2.0023 measures the amount of the orbital momentum in the ground crystal-field state due to mixing of higher crystal-field states. In copper-based magnets, for example, one typically finds

### 2.2. EPR spectrum

In the high-temperature limit, where thermal energy is larger than the Zeeman energy splitting (in the conventional X-band at 9.5 GHz the Zeeman spitting corresponds to the temperature of 0.45 K), the EPR absorption spectrum is determined in the linear-response theory by thermal-averaged (denoted by

where the spin ladder operators are given by *H*′. Rewriting

Eq. (8) reveals an interesting result that the resonant absorption is peaked at the Larmor frequency *ħ* is the reduced Planck constant. Moreover, in the case of no anisotropy, there is no time dependence of the spin correlation functions in the interaction representation (Eq. (8)), therefore the EPR spectrum simply consists of two *δ*-functions. The time dependence of the correlation functions in Eq. (8), which is due to magnetic anisotropy *H*′, is thus solely responsible for finite line widths of the EPR spectra and their shifts from the Larmor frequency. This is an essential results, as it demonstrates that magnetic anisotropy is directly reflected in the shape of the EPR line spectrum, unlike in all other techniques capable of detecting the anisotropy, e.g., inelastic neutron scattering, where the signal is a combined effect of multiple factors. Usually the EPR line width is small compared to the Larmor frequency and one can neglect the contribution peaked at the negative frequency –ω_{0}.

According to the KT theory, the EPR spectrum can be expressed as the Fourier transform of the relaxation function

Thus, spin correlations embedded into the relaxation function determine the EPR spectrum. The calculation of the relaxation function is however nontrivial. Therefore, approximation schemes are required. For Markovian random processes the relaxation function is approximated by [14]

Here, the spin correlation function is defined as

where *A* and *B*. Within the KT theory a Gaussian decay of the spin correlation function is postulated,

where the characteristic spin correlation time is determined by the dominant isotropic exchange *J*,

### 2.3. Exchange narrowing

Let us inspect two limiting cases of the correlation time with respect to the typical EPR time scale given by the Larmor frequency. For slow decay of correlations (*ψ*(*t*) in Eq. (10) can be replaced by its zero-time value, which is proportional to the second moment of the absorption line

This procedure yields a Gaussian relaxation function *φ*(*t*) and, according to Eq. (9), also a Gaussian shaped EPR spectrum, with the width

The fast decay limit (

which leads to an exponential decay of the relaxation function *φ*(*t*). Consequently, the EPR spectrum has a Lorentzian shape. Its line width is *M*_{2}, is opposed by the isotropic exchange interaction causing narrowing of the EPR line. The spin correlation time

which yields the full width at half maximum (FWHM) of the Lorentzian EPR line

The exchange-narrowing limit is typically applicable to real exchange-coupled system, except in cases of small couplings and high Larmor frequencies. We recall again that in the most widespread X-band EPR experiment (

## 3. Applications of the KT theory

Applications of the KT theory to experiments are numerous. Here, we will highlight a few cases from recent literature, where determination of the magnetic anisotropy turned out to be crucial for understanding the magnetic ground state. All examples concern magnetically frustrated spin lattices in 2D, where short-range spin interactions are incompatible with the underlying spin lattices, effectively suppressing long-range spin ordering and leading to unconventional states of matter. In such cases magnetic anisotropy, even if only being a small perturbation to the dominant isotropic exchange interaction, can tip the balance in favor of one or another competing ground state.

### 3.1. Kagome lattice

The first example is the 2D spin lattice in herbertsmithite, ZnCu_{3}(OH)_{6}Cl_{2} [20], a compound that has earned the reputation of being the best experimental realization of a quantum kagome antiferromagnet (QKA) of corner-sharing triangles (Figure 1), where the geometrical frustration is the most severe [21]. Numerous theoretical studies that proposed various different ground states over the last two decades, now seem to have converged on a gapped quantum spin liquid (QSL) – a state that is disordered, yet highly entangled [21]. Experimental signatures of such a state have also been lately advocated, although the bulk of experiments on this and the majority of other known QKA representatives actually speaks in favor of a gapless QSL. This discrepancy may well be related to perturbations beyond the isotropic Heisenberg exchange model on the kagome lattice, such as magnetic anisotropy.

The magnetic anisotropy of herbertsmithite was successfully determined by EPR in Ref. [22]. Based on relatively small *g*-shifts (of the order of 15%, as typical for Cu^{2+} ions [4, 5]), it was argued that the antisymmetric DM interaction (Eq. (5)) dominates the magnetic anisotropy in this compound. The DM vector pattern (Figure 1), which is determined by the symmetry of the kagome lattice, then according to Eq. (16) predicts the following angular dependence of the EPR line width [22]

where *θ* represents the polar angle between the normal to kagome planes and the applied magnetic field direction, while *dp* and *dz* are the in-plane and the out-of-plane components of the DM interaction. This expression is valid only in the infinite-temperature limit, therefore the authors applied it to fit the room-temperature EPR spectrum of herbertsmithite (Figure 1), where the EPR line width was shown to saturate to a constant value [22]. Thus the minute in-plane DM component

In contrast to herbertsmithite, in another QKA representative, vesignieite, BaCu_{3}V_{2}O_{8} (OH)_{2}, a long-range magnetic order was observed [26], which could be due to the fact that this systems is positioned in the ordered part of the above-mentioned phase diagram. In order to verify this conjecture, the same EPR analysis (Figure 2) as the one presented above for herbertsmithite was performed in Ref. [27]. The derived DM components are somewhat different from those in herbertsmithite, as the in-plane component dominates in vesignieite, *J* and should be much more efficient in suppressing quantum fluctuations than the out-of-plane DM component [27]. This could explain why magnetic ordering in vesignieite occurs at surprisingly high temperature for a frustrated system,

### 3.2. Triangular lattice

A regular triangular lattice of edge-sharing triangles is another example of a highly frustrated spin lattice in 2D. Contrary to the kagome lattice, where each spin in surrounded by four nearest neighbors, on the triangular lattice there are six such neighbors, which reduces the amount of frustration. Consequently, the triangular lattice exhibits a magnetically ordered ground state, which is, however, much more complex than on ordinary bi-partite spin lattices.

A slightly more complicated triangular lattice is realized in Fe-langasite, Ba_{3}NbFe_{3}Si_{2}O_{14} [28]. Here the Fe^{3+} (*S* = 5/2) spins reside on vertices of equilateral triangles arranged into a 2D triangular lattice (Figure 3). Quite interestingly, the magnetically ordered ground state below *TN* = 26 K is doubly chiral, as the same 120° spin configuration on each triangle is helically modulated from plane to plane [28].

To identify the anisotropy term that is responsible for such chirality of the magnetic ground state, an EPR study was again conducted [29]. The room-temperature EPR signal was found to exhibit a pronounced angular dependence of the EPR line width and line position. The former could be related either to zero-field-splitting anisotropy (Eq. (6)) or DM exchange anisotropy (Eq. (5)), with the anisotropy patterns as shown in Figure 3. The two models could not be distinguished solely based on the EPR response of the system. However, a combined study of the EPR spectra and antiferromagnetic resonance (AFMR) modes observed below *TN* suggested the DM interaction as the dominant source of anisotropy and thus to be responsible for the observed chiral behavior of Fe-langasite. The out-of-plane DM component

which was fit to the experimentally determined line-width anisotropy at 500 K (Figure 3). Here the constants

A later study combining EPR, AFMR and inelastic neutron scattering refined the anisotropy model in Fe-langasite and showed that actually both the DM anisotropy (

## 4. Pitfalls of the KT theory and alternative approaches

Although the above examples nicely demonstrate the value of the KT theory, this theory should be applied to each particular case with caution, because it is limited in several aspects. Firstly, the KT approach does not take into account a possible hidden symmetry of the DM interaction (Section 4.1) and diffusional decay of spin correlations in low dimensional spin systems (Section 4.2). Secondly, the EPR moments (Eq. (13) and Eq. (15)) implicitly employ four-spin correlation functions, which can be explicitly evaluated only in the infinite-temperature limit, where spin correlations between neighboring sites are negligible. On the other hand, the analysis of the EPR line width at temperatures of the order of the dominant exchange coupling *J* and below requires different approaches, like the Oshikawa-Afflect effective-field-theory (Section 4.3). Lastly, one should keep in mind that the KT theory is perturbative, therefore the cases of large (or even dominant) magnetic anisotropy should be treated with different approaches (Section 4.4).

### 4.1. Reducibility of the DM interaction

It was found theoretically that the DM interaction may in some cases possess a hidden symmetry [31], in the sense that it can be effectively transformed into a term with the symmetry of the anisotropic exchange and with reduced magnitude of *d* are those that sum up to zero within any closed loop on a spin lattice; for example, for the kagome lattice, the in-plane DM component *Dp* is reducible, while the out-of-plane component *Dz* is irreducible. The KT theory remains applicable for the irreducible components of the DM interaction.

### 4.2. Spin diffusion

In low-dimensional magnetic systems it may happen that the Gaussian approximation of the decay of the spin correlation function in Eq. (12) is not justified due to a diffusional contribution to the decay. This dictates slower time dependence of the form [34].

where *D* represents the dimensionality of the spin system. For

When the secular part of the anisotropy Hamiltonian (Eq. (3)), i.e., the part commuting with the Hamiltonian *Γ*. On the other hand, there exists no universal picture for two dimensions. Nevertheless, deviations of the experimentally observed EPR spectra from the Lorentzian shape in 2D were observed in the past and successfully ascribed to the presence of spin diffusion [35]. Quite generally, the spin-diffusion effect in two dimensions is usually much less pronounced than in one dimension.

Finally, we note that the diffusional decay of the electronic spin correlation functions is often not detectable by EPR at all, even in low dimensional systems. Although these systems may be characterized as being low dimensional due to the dominant exchange interaction along a chain or within a plane, also inter-chain/inter-layer exchange couplings can still be large compared to the magnetic anisotropy terms regulating the linewidth of the EPR absorption spectra. In such cases the decay of spin correlations is effectively taking place in three dimensions and the spin-diffusion problem is absent.

### 4.3. Exact calculations on finite clusters

The postulate of the Gaussian decay of the spin correlation function in the KT theory (Eq. (12)) has no theoretical background and is not necessarily valid, as explained in Section 4.1. However, this assumption is not needed at all if the EPR line shape is calculated from the basics, i.e., from Eq. (7). This can be done only on finite clusters of spins. Such a limitation then requires an extrapolation to the thermodynamic limit if these calculations are to be applied to macroscopic samples.

Exact calculations of the EPR line shape on finite clusters were performed by El Shawish et al. [36] for certain 1D and 2D spin lattices. For a spin chain, the results showed a noticeable transformation of the decay of the spin correlations from the Gaussian shape at early times to a much slower decay of diffusional characteristics at longer times. The resulting line broadening and the deviation from the Lorentzian line shape were, however, later shown to be effectively short-cut by inter-chain exchange [37].

The situation is very different in 2D, e.g., for the kagome spin lattice. Namely, the finite-cluster calculations revealed that, at least for the irreducible DM component *dz* the line width indeed scales with

We note that the finite-cluster approach is severely limited, as the extrapolation to the thermodynamic limit, which is usually of interest in experiments, is highly nontrivial and depends on a particular spin lattice [36]. However, since the results are exact, this approach may still be very interesting for small systems, such as molecular magnets. An interesting prediction of a double-peak EPR spectrum was also given (Figure 5). The spectrum should thus strongly differ from the usual Lorentzian line shape, which still awaits experimental confirmation.

### 4.4. Oshikawa-Affleck theory

Exact calculations of the second and fourth moments of the EPR absorption spectra (Eq. (13) and Eq. (15)) are possible within the KT framework for infinite lattices, but only in the limit of infinite temperature. In this case, static spin correlations of the products of spin operators acting on different lattice sites can be neglected. In general, in Eq. (13) and Eq. (15), one is dealing with the computation of four-spin correlation functions since the magnetic anisotropy Hamiltonian is quadratic in spin operators. Therefore, special schemes of disentangling the four-spin correlation functions into products of two-spin correlation functions need to be applied [14, 19]. Further complications emerge at finite temperatures, i.e., at

The problem of how finite temperatures (finite spin correlations) affect the EPR line width is treated within the Oshikawa-Afflect effective-field-theory approach that is applicable to spin chains [38, 39]. The spin diffusion picture, which predicts a non-Lorentzian line shape in 1D, does not apply to the OA theory. In contrast to the KT theory, this approach works well at intermediate and low temperatures,

where the constant *c*_{s} originates from the DM interaction and/or from a staggered *g* factor. This broadening is of the form

The temperature dependence of the DM broadening is inverse to the AE broadening, as the former decreases with increasing temperature while the latter increases. Moreover, while the AE broadening effect is independent of the applied field, the DM broadening increases with the square of the applied field.

If both the AE and the DM term are of similar magnitude in a particular system, one can expect to observe both EPR broadening mechanisms simultaneously. Such is, for instance, the case in the CuSe_{2}O_{5} spin-chain compound [40]. There, simultaneous modeling of the angular, temperature, and frequency-dependent EPR line width with the OA theory (the sum of contributions in Eq. (20) and Eq. (21)) allowed Herak et al. to extract both the AE and the DM anisotropy constants [40]. The simultaneous fits of the AO theory to multiple experimental datasets are presented in Figure 6.

At the end, it should be stressed that the OA approach still relies on the perturbation theory (in magnetic anisotropy). So, cases, where the anisotropy is of the order of the isotropic exchange interactions or larger are untreatable within this theory.

### 4.5. EPR moments

The EPR-moments approach [41] described in this section is more direct, i.e., non-perturbative. Within this approach the line width (and line shift) in the “frequency domain,” where the frequency is varied in a fixed Zeeman field, can be calculated for an arbitrary strength of magnetic anisotropy. Moreover, exact calculations at any temperature are possible for spin chains. In general, the EPR line width is given by the four lowest shifted moments

where

The moments in the frequency domain (Eq. (22)) represent static correlations that can be calculated in the case of 1D spin chains to arbitrary precision for any temperature and applied field [41]. The agreement of this approach with fully numerical calculation for finite chain Hamiltonians is shown in Figure 7.

In a typical EPR experiment, however, the frequency is kept constant and the magnetic field is swept. It turns out that the calculation of the shifted moments in the field domain,

requires the knowledge of infinitely many moments in the frequency domain. Therefore, the experimental line width cannot be calculated exactly at arbitrary temperature and field [40].

We should stress that the fact the EPR spectrum is measured in a field-sweet experiment is actually neglected in almost all theoretical treatments. Furthermore, a complication that arises when applying the EPR-moments approach to an experiment is that “long tails” of the EPR line may considerably contribute to the moments, while these are usually not properly accounted for by the experimentally determined FWHM due to noise [42]. Therefore a cutoff of high-frequency tails is necessary, as it was recently demonstrated for the case of the quasi-one-dimensional magnet Cu(py)_{2}Br_{2} [42].

## 5. Conclusions

This chapter reviews the development of the treatment of the EPR absorption line in strongly exchange-coupled spin systems. The starting point is the Kubo Tomita general theory of magnetic resonance absorption, which demonstrates how the line width can be approximated by two lowest even moments of the EPR line, M_{2} and M_{4}. We note that the knowledge of all the moments,

The KT approach was successfully applied to various spin lattices in the past, including the geometrically frustrated kagome and triangular lattices, which are exemplified here. However, when the theory is applied to a particular system special attentions needs to be made a) to a possible reducibility of the asymmetric Dzyaloshinskii-Moriya exchange anisotropy, b) to the diffusional decay of spin correlations, which may occur in low-dimensional spin systems, c) to finite correlations among spins at different sites, which typically develop below the temperature of the order of the dominant isotropic exchange, and d) to the size of the magnetic anisotropy, which is only treated as a perturbation in the KT theory. All these drawbacks of the KT theory can be overcome, at least in special cases. In this review, special approaches that were developed in this vein have been summarized. These include a) exact calculations of the EPR line on finite clusters, the Oshikawa-Affleck effective-field theory for 1D spin systems, and the recently developed EPR-moments approach. For each approach a representative example has been provided in this review.

## Acknowledgments

The author acknowledges the financial support of the Slovenian Research Agency under the program No. P1-0125.

## Conflict of interest

The author declares no conflict of interests.