Outline
Comptes Rendus

Effective spin-wave action for ordered Heisenberg antiferromagnets in a magnetic field
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. 60-64.

Abstract

We derive the effective long-wavelength Euclidean action for the antiferromagnetic spin-waves of ordered quantum antiferromagnets subject to a uniform magnetic field. We point out that the magnetic field-dependence of the spin-wave dispersion predicted by the usual O(3)-quantum nonlinear sigma model disagrees with spin-wave theory. We argue that the nonlinear sigma model does not take into account all relevant spin-wave interactions and derive a modified effective action for the long-wavelength spin-waves which contain an additional quartic interaction. At zero temperature the corresponding vertex is relevant in the renormalization group sense below three dimensions.

Metadata
Received:
Accepted:
Published online:
DOI: 10.1016/j.crci.2006.06.008
Keywords: Heisenberg model, Magnetic system, Magnetic field effects

Nils Hasselmann 1; Florian Schütz 1; Ivan Spremo 1; Peter Kopietz 1

1 Institut für Theoretische Physik, Universität Frankfurt, Max-von-Laue-Strasse 1, 60438 Frankfurt/Main, Germany
@article{CRCHIM_2007__10_1-2_60_0,
     author = {Nils Hasselmann and Florian Sch\"utz and Ivan Spremo and Peter Kopietz},
     title = {Effective spin-wave action for ordered {Heisenberg} antiferromagnets in a magnetic field},
     journal = {Comptes Rendus. Chimie},
     pages = {60--64},
     publisher = {Elsevier},
     volume = {10},
     number = {1-2},
     year = {2007},
     doi = {10.1016/j.crci.2006.06.008},
     language = {en},
}
TY  - JOUR
AU  - Nils Hasselmann
AU  - Florian Schütz
AU  - Ivan Spremo
AU  - Peter Kopietz
TI  - Effective spin-wave action for ordered Heisenberg antiferromagnets in a magnetic field
JO  - Comptes Rendus. Chimie
PY  - 2007
SP  - 60
EP  - 64
VL  - 10
IS  - 1-2
PB  - Elsevier
DO  - 10.1016/j.crci.2006.06.008
LA  - en
ID  - CRCHIM_2007__10_1-2_60_0
ER  - 
%0 Journal Article
%A Nils Hasselmann
%A Florian Schütz
%A Ivan Spremo
%A Peter Kopietz
%T Effective spin-wave action for ordered Heisenberg antiferromagnets in a magnetic field
%J Comptes Rendus. Chimie
%D 2007
%P 60-64
%V 10
%N 1-2
%I Elsevier
%R 10.1016/j.crci.2006.06.008
%G en
%F CRCHIM_2007__10_1-2_60_0
Nils Hasselmann; Florian Schütz; Ivan Spremo; Peter Kopietz. Effective spin-wave action for ordered Heisenberg antiferromagnets in a magnetic field. 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. 60-64. doi : 10.1016/j.crci.2006.06.008. https://comptes-rendus.academie-sciences.fr/chimie/articles/10.1016/j.crci.2006.06.008/

Original version of the full text

1 Introduction

Recently, two of us derived a new effective field theory [1] for transverse spin fluctuations in ordered quantum antiferromagnets. Our approach is based on the Holstein–Primakoff transformation [2], which maps the original spin Hamiltonian onto a bosonic many-body Hamiltonian. However, rather than working directly with the canonical Holstein–Primakoff bosons, we express them in terms of Hermitian field operators representing staggered and uniform transverse spin fluctuations [3]. The main advantages of this procedure are as follows: (a) the suppression of the effective interaction between antiferromagnetic magnons at long wavelengths is manifest; (b) the relation between Holstein–Primakoff bosons and the fields of the nonlinear sigma model [4] (NLSM) representing the transverse staggered spin fluctuation is made precise; (c) the zero wave vector modes of the uniform and staggered magnetizations can be easily treated within this approach [3], which is important in the study of finite magnets; (d) the transverse ferromagnetic fluctuations may be eliminated to yield an effective theory involving only transverse antiferromagnetic fluctuations. This is most easily done within a path integral approach; the resulting effective action coincides with the NLSM at the Gaussian level. Although the interaction terms of the effective action differ from those of the NLSM, the renormalization group flows of both theories agree to one-loop order [1].

Here, we extend the Hermitian operator approach to antiferromagnets in a uniform magnetic field. In this case, the differences between this approach and the NLSM become more substantial. Already at the Gaussian level the two theories differ: while the Hermitian operator approach reproduces by construction the correct spin-wave dispersion as known from spin-wave theory, the NLSM does not. More severely, we show below that a quartic interaction term present in the Hermitian operator approach is absent in the NLSM.

We shall constrain the analysis here to a Heisenberg model on a hypercubic lattice with nearest neighbor antiferromagnetic coupling J > 0. It is, however, straightforward to extend the results to general bipartite lattices. This would also cover a recently synthesized metal–organic S = 5/2 quantum antiferromagnet on a distorted honeycomb lattice [5], which is bipartite. In presence of a uniform external field h (measured in units of energy), the Hamiltonian is given by

Hˆ=JijSi·Sjih·Si(1)
Here, Si is spin operator normalized such that Si2=S(S+1), and the first sum is over nearest neighbor pairs ij of a hypercubic lattice with lattice constant a.

Our main interest is a description of the low-energy and long-wavelength properties of the Heisenberg model (1). A standard approach to this problem is the O(3) quantum NLSM whose effective action describes only the relevant staggered spin fluctuations. According to Refs. [6,7], a magnetic field can be taken into account within the NLSM approach by replacing the derivative ∂τ with respect to the imaginary time τ with a covariant derivative (or Lie derivative, see Ref. [8])

ττih×(2)
With this minimal coupling the NLSM in a uniform magnetic field has the form
SNLSM[Ω]=ρ020βτDr[μ=1D(μΩ)2+1c02(τΩih×Ω)2](3)
where the unit vector Ω(τ,r) represents the slowly fluctuating staggered magnetization, ρ0 and c0 are the spin stiffness and spin-wave velocity at temperature T = 1/β = 0, and ∂μ = ∂/∂rμ is the spatial derivative in direction μ = 1,…, D. Numerical values of ρ0 and c0 must be computed microscopically, e.g. by means of a 1/S-expansion [4]. To leading order in 1/S, one has ρ0 ≈ JS2a2−D and c0 ≈ 2D1/2JSa. For later reference, we note that the classical uniform transverse susceptibility is χ0 = ρ0/c02 with χ0 ≈ (4DJaD)−1 to leading order in 1/S. The classical uniform magnetization per volume is M0=χ0h.

2 Spin-wave dispersion

In the absence of a magnetic field, the NLSM predicts a doubly degenerate antiferromagnetic spin-wave mode with long-wavelength dispersion Ek=c0|k|, in agreement with linear spin-wave theory [4]. Rather surprisingly, however, in the presence of a magnetic field the long-wavelength dispersion derived from the NLSM in Eq. (3) disagrees with spin-wave theory. To see this, let us recall the spin-wave dispersion predicted by linear spin-wave theory. In this approach the spin operators are expanded in deviations from the classical ground state spin configuration, using either the Holstein–Primakoff [2] or the Dyson–Maleev [9] transformation. The classical ground state of a quantum antiferromagnet in a uniform external field is canted, where the staggered magnetization is perpendicular to the magnetic field, as shown in Fig. 1. To leading order in 1/S, spin-wave theory [5,10,11] predicts two transverse spin-wave modes, one gapless and one gapped, with long-wavelength dispersions (see Eqs. (17) and (20) below)

Ek+[h2+c+2k2]1/2,Ekc|k|(4)
where the spin-wave velocities are
c+=c013m02,c=c01m02(5)
Here, m0=|m0| is the length of the normalized classical uniform magnetization
m0=aDM0S=h4DJS(6)
which is related to the classical canting angle ϑ0 between the direction of the staggered magnetization and the local magnetic moments by m0 = sin ϑ0. In general, the classical canting angle is renormalized by quantum fluctuations [5,10].

Fig. 1

Spin configuration Si=Smˆi in the classical ground state of a two-sublattice antiferromagnet subject to a uniform magnetic field h = hex, where ri belongs to sublattice A and rj to sublattice B. Here m0=h/(4DJS) and n0=1m02.

Let us compare Eqs. (4) and (5) with the dispersion relation obtained from the NLSM. Denoting by eα the unit vectors in the fixed directions α = x, y, z and assuming h=hex to point along the x-axis, we write

Ω=Πex+Πey+[1Π2Π2]1/2ez(7)
and expand the action SNLSM[Ω] to quadratic order in the fluctuations Π and Π transverse to the direction of the staggered magnetization. One easily finds from Eq. (3) that the spin-wave mode Π polarized parallel to the magnetic field is gapped with energy dispersion
Ek=[h2+c02k2]1/2(8)
while the spin-wave mode Π polarized perpendicular to the magnetic field has a gapless linear dispersion,
Ek=c0|k|(9)

The spin-wave spectrum of the NLSM given in Eqs. (8) and (9) does not reproduce the correct magnetic field-dependence of the spin-wave spectrum in Eqs. (4) and (5) obtained via conventional spin-wave theory. Obviously, the above minimal coupling (2) does not account for order (h/J)2 corrections to the dispersion, even though the qualitative aspects of the dispersion are correctly described.

3 Effective action

In Ref. [1] we have developed a new method to derive the effective action for staggered transverse spin fluctuations. We here sketch the main steps of the derivation adapted to a quantum antiferromagnet subject to a uniform magnetic field [12]. Starting point is the representation of spin operators in terms of canonical boson operators bi and bi using the Holstein–Primakoff transformation [2]. Let us denote by mˆi=Si/S the directions of the spins in the true ground state and introduce a local right-handed orthogonal triad of unit vectors ei(1),ei(2),mˆi. As explained in Ref. [13] there is an U(1) gauge freedom in the choice of the transverse vectors ei(1) and ei(2). Defining the spherical basis vectors eip=ei(1)+ipei(2), p = ±, we express the components of the spin operator Si in terms of canonical boson operators bi and bi as follows [2]:

Si=Simˆi+Si=Simˆi+12p=±Sipeip(10)
with
Si=Sni,ni=bibi(11a)
Si+=(2S)1/2(1ni2S)1/2bi(11b)
and Si=(Si+). The Heisenberg model (1) can then be written as a bosonic many-body Hamiltonian [5]. We expand around the classical ground state configuration (see Fig. 1)
mˆi=ζin0ez+m0ex(12)
and choose
ei(1)=ey,ei(2)=ζin0ex+m0ez(13)
where ζi = 1(ζi = −1) for riA(riB). Here, n0=1m02=cosϑ0. To diagonalize the quadratic part of this Hamiltonian, we first perform a Fourier transformation in the sublattice basis
bi={(2/N)1/2keik·riAk,riA(2/N)1/2keik·riBk,riB(14)
where the momentum sums are over the reduced (antiferromagnetic) Brillouin zone. After introducing the symmetric and antisymmetric combinations Ckσ=21/2[Ak+σBk], σ = ±, we apply the Bogoliubov transformation
(CkσCkσ)=(ukσσvkσσvkσukσ)(ΨˆkσΨˆkσ)(15)
where
ukσ=[1+σm02γk+εkσ2εkσ]1/2(16a)
vkσ=[1+σm02γkεkσ2εkσ]1/2(16b)
with
εkσ=[(1+σm02γk)2(n02γk)2]1/2(17)

Here, γk=D1μcos(k·aμ), where aμ, μ = 1,…,D, are the D primitive lattice vectors of the hypercubic lattice. The quadratic part of the effective boson Hamiltonian assumes then the form of non-interacting harmonic oscillators [5,12],

Hˆ2=k,σEkσ[ΨˆkσΨˆkσ+12](18)
where the operators Ψˆkσ satisfy the canonical bosonic commutation relations,
[Ψˆkσ,Ψˆkσ]=δk,kδσ,σ(19)
and the energy dispersion of the two spin-wave modes is given by
Ekσ=2DJSεkσ(20)

In the limit of long wavelengths Eq. (20) reduces to the energy dispersions given in Eqs. (4) and (5).

To derive the long-wavelength effective action for the staggered transverse spin fluctuations, we need the precise relation between the Bogoliubov quasi-particle operators Ψˆkσ and the field operators Πˆkσ representing the transverse fluctuations of the staggered magnetization. In a Euclidean path integral approach, these operators correspond to continuum fields Πkσ(τ) which are analogs of the fields Πk(τ) and Πk(τ) of the NLSM. The relation between these two different parameterizations of the spin fluctuations is simply [1]

Ψˆk+=i(χ0/2VEk+)1/2[Ek+Πˆk++iχ01Φˆk+](21a)
Ψˆk=(χ0/2VEk)1/2[EkΠˆk+iχ01Φˆk](21b)
where V = NaD is the volume of the system. One easily verifies that the pairs Πˆσ, Φˆσ satisfy commutation relations of canonically conjugate bosonic field operators,
[Πˆkσ,Φˆkσ]=iVδk,kδσ,σ(22)

The field operators Πˆkσ and Φˆkσ have a simple interpretation in terms of transverse staggered and uniform spin components. The Fourier transformed spin operators on the A and B sublattices are given by

SA/B,k=(2/N)1/2riA/Beik·riSi(23)

Using Eqs. (10) and (13) we find that, to leading order in 1/S and m0, the transverse staggered components of the spins are given by1

Sst,k(1)=12(SA,k(1)SB,k(1))(S/aD)N1/2Πˆk(24a)
Sst,k(2)=12(SA,k(2)+SB,k(2))(S/aD)N1/2Πˆk+(24b)

In contrast, the transverse components of the uniform magnetization are related to the operators Φˆkσ,

Sk(1)=12(SA,k(1)+SB,k(1))N1/2Φˆk+(25a)
Sk(2)=12(SA,k(2)SB,k(2))N1/2Φˆk(25b)

In terms of these operators, our quadratic spin-wave Hamiltonian (18) can be written as

Hˆ2=12Vk,σ[χ01ΦˆkσΦˆkσ+χ0Ekσ2ΠˆkσΠˆkσ](26)

The effective Euclidean action Seff[Πσ] for the staggered spin fluctuations can now be obtained by writing the partition function as a phase space path integral over the fields Πkσ(τ) and Φkσ(τ) associated with the above operators and subsequently integrating over the Φσ-fields which represent gapped ferromagnetic fluctuations. At the level of a Gaussian approximation, we obtain in this way Seff ≈ Seff(2) with

Seff(2)[Πσ]=χ02βVK,σ(ωn2+Ekσ2)ΠKσΠKσ(27)

Here, we have combined momenta k and bosonic Matsubara frequencies ωn in a composite label K=(k,iωn). We have further defined

ΠKσ=0βτeiωnτΠkσ(τ)(28)

At long wavelengths the effective action (27) has the same form as the corresponding Gaussian part of the action of the NLSM. However, in contrast to the NLSM (3), our effective action (27) has the correct spin-wave dispersion, even for short wavelengths.

We have calculated the leading corrections to the effective action Seff[Πσ] for staggered spin fluctuations arising from spin-wave interactions in the Holstein–Primakoff approach [12]. Keeping only interaction terms which become relevant below three dimensions and thus omitting terms which are marginal in D = 1, we find in the continuum limit

Seff[Πσ]=χ02βVK,σ(ωn2+cσ2k2+rσ)ΠKσΠKσiχ0h0βτDrΠ+2τΠχ0h280βτDrΠ+2(Π+2+Π2)(29)

Here, the cσ is the Holstein–Primakoff results for the spin-wave velocities, which for large S is given in Eq. (5). To leading orders in 1/S, the values of the gap parameters are r+ = h2 and r = 0. For symmetry reasons, the true spin-wave spectrum of the Π-mode must remain gapless at vanishing wave vector (Goldstone mode) [14], while the spin-wave gap of the Π+-mode is not renormalized [14,15]. In a renormalization group analysis, this requires fine tuning such that at the fixed point these conditions are met. While the cubic term in Eq. (29) corresponds precisely to the cubic (Berry phase) term in the NLSM (3), the quartic interaction between the field components is absent if the magnetic field is included in the NLSM by the minimal coupling (2).

4 Summary and conclusion

In this work we have used the Holstein–Primakoff transformation to derive an effective action for staggered transverse spin fluctuations of quantum Heisenberg antiferromagnets in uniform magnetic fields. Our effective action contains an additional quartic interaction between the field components which is not contained in the NLSM. It is easy to see that at zero temperature the quartic interaction vertex in Eq. (29) is relevant in the renormalization group sense for dimensions D < 3, so that we expect that it generates singularities in perturbation theory [12]. We conclude that the NLSM given in Eq. (3) does not contain all relevant interactions in the ordered phase below three dimensions. The reason for this is that not all approximations that are usually made in the derivation of the NLSM [6] are justified in the presence of a uniform magnetic field. In particular, it is not justified to neglect Umklapp scattering processes between transverse and longitudinal spin fluctuations involving momentum transfers across the boundary of the magnetic Brillouin zone. As an alternative to Eq. (29) it should be possible to analyse the long-wavelength staggered spin fluctuations within a phenomenological Ginzburg–Landau model with a quartic interaction term [16] since this model treats the longitudinal modes as independent degrees of freedom.

Acknowledgements

The authors thank Ian Affleck for illuminating discussions. This work was supported by the DFG via Forschergruppe FOR 412.

1 We use here a basis which is rotated compared to the one used in Ref. [1]. This leads to different appearances of the minus signs in the transformation Eqs. (24) and (25).


References

[1] N. Hasselmann; P. Kopietz Europhys. Lett., 75 (2006), p. 1067

[2] T. Holstein; H. Primakoff Phys. Rev., 58 (1940), p. 1098

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

[4] S. Chakravarty; B.I. Halperin; D. Nelson Phys. Rev. B, 60 (1988), p. 1075

[5] I. Spremo; F. Schütz; P. Kopietz; V. Pashchenko; B. Wolf; M. Lang; J.W. Bats; C. Hu; M.U. Schmidt Phys. Rev. B, 72 (2005), p. 174429

[6] S. Sachdev Quantum Phase Transitions, Cambridge University Press, Cambridge, 1999

[7] D.S. Fisher Phys. Rev. B, 39 (1989), p. 11783

[8] J.L. McCauley, Cambridge University Press, Cambridge, 1997 (See, for example Classical Mechanics p. 212)

[9] F.J. Dyson; S.V. Maleev Sov. Phys. JETP, 102 (1956), p. 1217

[10] M.E. Zhitomirsky; T. Nikuni Phys. Rev. B, 57 (1998), p. 5013

[11] M.E. Zhitomirsky; A.L. Chernychev Phys. Rev. Lett., 82 (1999), p. 4536

[12] N. Hasselmann, P. Kopietz (in preparation).

[13] F. Schütz; M. Kollar; P. Kopietz Phys. Rev. Lett., 91 (2003), p. 017205

[14] D.I. Golosov; A.V. Chubukov Fiz. Tverd. Tela (Leningrad), 30 (1988), p. 893

[15] M. Oshikawa; I. Affleck Phys. Rev. B, 65 (2002), p. 134410

[16] I. Affleck; I. Affleck Phys. Rev. B, 41 (1990), p. 6697


Comments - Policy