Magnetic excitations

Here we discuss the coherent magnetic excitations known as spin waves. The first section deals with spin waves in a ferromagnet, then we extend the discussion to the much more difficult antiferromagnetic spin waves. We end the chapter by touching the topic of quantum magnets.

*Spin waves in a ferromagnet

Let us assume that the system can be described by a Heisenberg Hamiltonian of the form shown on the Elastic magnetic scattering page.

We assume that the sign of the interaction, \(J_{\delta}\), is dominantly positive; i.e. that the Fourier transform of the spin-spin interaction, \(J({\bf q}')\), is maximal at \(q'=0\). (\(J(\bf q)\) is defined in this equation on the Magnetic neutron scattering page.) In that case, \({\bf Q} = 0\), the ground state of the system is the ferromagnet

\begin{equation} \label{dummy791517836} | GS \rangle = | {\rm FM} \rangle = | \uparrow\uparrow\uparrow\uparrow \rangle , \end{equation}

where all spins point along the positive \(z\)-axis. This quantization axis is in practice defined by the direction of an applied field or some small anisotropy terms, not present in Heisenberg form. We will here neglect anisotropies and just assume that a quantization axis exists. We will describe the elementary excitations in the ferromagnet in two ways: The equation-of-motion method and the stationary operator method, both of which are described below.

Solving the Hamiltonian (\ref{eq:Heisenberg}) quantum mechanically is an easy task for a ferromagnet. Let us first rewrite the spin product term: \begin{equation} {\bf S}_i \cdot {\bf S}_j = s_i^z s_j^ z + \frac{1}{2} (s_i^+ s_j^- + s_i^- s_j^+) , \end{equation} where the spin raising and lowering operators are defined as \begin{align} \label{eq:spinraising} s_{j}^+ &= s_{j}^x + i s_{j}^y , \\ s_{j}^- &= s_{j}^x - i s_{j}^y . \end{align} The raising/lowering operators have the effect of changing \(m\), the eigenvalue of \(s^z\) for one spin, by one unit: \begin{align} s_{j}^+ |S,m\rangle_j &= \sqrt{(S-m)(S+m+1)} |m+1\rangle_j , \\ s_{j}^- |S,m\rangle_j &= \sqrt{(S-m+1)(S+m)} |m-1\rangle_j . \end{align}

We can see that when using the Hamiltonian (\ref{eq:Heisenberg}) on the ferromagnetic state, \(|{\rm FM}\rangle = |\uparrow \uparrow \uparrow \uparrow \rangle\), all spin terms containing \(s^+\) vanish, leaving only the diagonal \(s^z_i s^z_j\) terms. \begin{equation} H | {\rm FM} \rangle = - \sum_{\langle ij \rangle} J_{ij} S^2 | {\rm FM} \rangle . \end{equation} The ferromagnetic state is thus an eigenstate of the Hamiltonian, and is also the ground state. However, for antiferromagnetic couplings, (\ref{eq:Heisenberg}) is a very difficult many-body problem. This can be seen from just considering the two-spin model for spins \(S=1/2\). Here, the classical ground states are \(| {\rm AFM}_1 \rangle = | \uparrow \downarrow \rangle\) and \(| {\rm AFM}_2 \rangle = | \downarrow \uparrow \rangle\). None of these states are eigenstates to (\ref{eq:Heisenberg}), since terms in the Hamiltonian of the type \(s_1^+ s_2^-\) mix the states: \begin{align} H | {\rm AFM}_1 \rangle &= - \frac{J}{4} |{\rm AFM}_1 \rangle + \frac{J}{2} |{\rm AFM}_2 \rangle , \\ H | {\rm AFM}_2 \rangle &= + \frac{J}{2} |{\rm AFM}_1 \rangle - \frac{J}{4} |{\rm AFM}_2 \rangle . \end{align} The resulting ground state is a spin singlet: \begin{equation} |{\rm GS}\rangle = \frac{1}{\sqrt{2}} \left({|\rm AFM}_1 \rangle - |{\rm AFM}_2 \rangle \right) , \end{equation} with a ground state energy of \(-3J/4\), three times lower than (the expectation value of) the energy of the classical solution.

*The equation-of-motion method

Figure 1: Dispersion relation for spin waves in a 1D nearest neighbour ferromagnet in an applied field, given by equation \eqref{eq:omegaq_1Dnn}.

Spin excitation frequencies can in general be determined by the ``equation-of-motion approach \cite{marshall,Blundell}. We assume that the spin components perpendicular to the spin ordering \(z\)-direction, \(s_j^x(t)\) and \(s_j^y(t)\), can be non-zero, although their time average is zero. We shall then show that these perpendicular components will perform periodic motion.

To find the time dependence of the perpendicular components, one must take into account that the spins are angular momenta. Hence, a field (or interaction) trying to align the spin along the \(z\)-axis, instead causes the spins to precess around this axis. We could now write down the classical equations of motion for the spin system, but it turns out to be mathematically simpler to perform a quantum mechanical treatment, which we will therefore do.

We start by considering the spin component along the ordering direction, \(s_j^z(t)\). For the equation of motion of this spin operator, we use the Ehrenfest theorem:

\begin{equation} i \hbar \frac{d}{dt} s_j^z = [s_j^z, \hat{H}] \end{equation} to reach \begin{align} \label{eq:fm_motion_z} i \hbar \frac{d}{dt} s_j^z &= \sum_{j'} \left[s_j^z, J({\bf r}_j-{\bf r}_{j'}) (s_j^x s_{j'}^x + s_j^y s_{j'}^y)\right] \nonumber \\ &= \sum_{j'} J({\bf r}_j-{\bf r}_{j'}) \left(s_{j'}^x s_j^y - s_{j'}^y s_j^x \right) \approx 0. \end{align} We have in the derivation neglected commutator terms that immediately vanish. The last approximation to zero is valid since the terms inside the sum are all small to second order - and in fact can be shown to have a time average of exactly zero.

Thus, the spin value along the ordering directions is a constant of motion.

Next, we consider the linear combination \(s_j^+(t) = s_j^x(t) + i s_j^y(t)\), also known as the spin raising operator, defined in (\ref{eq:spinraising}). The Ehrenfest theorem yields:

\begin{align} \label{eq:fm_motion} i \hbar \frac{d}{dt} s_j^+ &= [s_j^+,\hat{H}] \nonumber \\ &= \sum_{j'} \left[(s_j^x + i s_j^y), J({\bf r}_j-{\bf r}_{j'}) (s_j^x s_{j'}^x + s_j^y s_{j'}^y + s_j^z s_{j'}^z - g \mu_{\rm B} B s_j^z )\right] \nonumber \\ &= g \mu_{\rm B} B s_j^+ + \sum_{j'} J({\bf r}_j-{\bf r}_{j'}) \left(s_{j'}^z s_j^+ - s_{j'}^+ s_j^z \right) . \end{align}

As often in solid state physics, we turn to the Fourier transform; this time of the spin operators \begin{equation} S_{\bf q'}^+ = \frac{1}{\sqrt{N}} \sum_j \exp(i {\bf q}'\cdot {\bf r}_j) s_j^+ , \end{equation} where \(N\) is the number of spins. Combined with (\ref{eq:fm_motion}), we arrive at \begin{equation} \label{eq:fm_motion2} \hbar \frac{d}{dt} S_{\bf q'}^+ = g \mu_{\rm B} B S_{\bf q'}^+ + \frac{1}{\sqrt{N}} \sum_{j,j'} \exp(-i {\bf q}'\cdot{\bf r}_j) J({\bf r}_j-{\bf r}_{j'}) \left(s_{j'}^z s_j^+ - s_{j'}^+ s_j^z \right) . \end{equation} Now, in the ground state, the operator \(s_j^z\) has the eigenvalue \(S\). At finite temperatures, we approximate this operator by its thermal mean value, the ordered moment \(\langle S^z \rangle\). Hence, we simplify (\ref{eq:fm_motion2}) \begin{align} \hbar \frac{d}{dt}S_{\bf q'}^+ &= g \mu_{\rm B} B S_{\bf q'}^+ + \frac{\langle S^z \rangle}{\sqrt{N}} \sum_{j,j'} \exp(i {\bf q}'\cdot{\bf r}_j) J({\bf r}_j-{\bf r}_{j'}) \left(s_j^+ - s_{j'}^+ \right) \nonumber \\ &= \big( g \mu_{\rm B} B + \langle S^z \rangle [J({\bf 0})-J({\bf q}')] \big) S_{\bf q'}^+ . \end{align} This differential equation shows that \(S_{\bf q'}^+\) performs a periodic motion with the frequency \begin{equation} \label{eq:fm_omegaq} \hbar\omega_{\bf q'} = g \mu_{\rm B} B + \langle S \rangle [J({\bf 0})-J({\bf q}')] . \end{equation} Now, let us investigate the significance of this result. Assume that \(S_{\bf q'}^+\) is the annihilation operator of one elementary excitation of the system, \begin{equation} S_{\bf q'}^+ |\bf q'\rangle = |{\rm GS}\rangle , \end{equation} where the energy of the groundstate is \(E_0\) and of the excited state \(| {\bf q}'\rangle\) is denoted \(E_0 + \hbar \omega_{\bf q'}\). The time dependent operator, \(S_{\bf q'}^+(t)\), is then found directly by \begin{equation} S_{\bf q'}^+(t) | {\bf q}' \rangle = \exp\left( \frac{i H t}{\hbar} \right) S_{\bf q'}^+ \exp\left( - \frac{i H t}{\hbar} \right) | {\bf q}' \rangle = \exp(-i \omega_{\rm q'} t) S_{\bf q'}^+ | {\bf q'} \rangle . \end{equation} Hence, \(\omega_{\bf q'}\) is the frequency of the spin deviation, and \(\hbar \omega_{\bf q'}\) is the energy of the elementary (spin wave) excitation of the system.

\paragraph{Example.} For nearest neighbour interactions in one dimension, each spin interacts only with its two neighbours. Hence \(J(q') = J(0) \cos(q'd)\), and \(J(0) = 2J\), giving \begin{equation} \label{eq:omegaq_1Dnn} \hbar \omega_{\bf q'} = g \mu_{\rm B} B + 2 J S \big( 1-\cos(q'd) \big) . \end{equation} Hence, for small values of \(|q'|\) (compared to the reciprocal lattice vector, \(2\pi/d\)), the spin wave energy in zero field becomes quadratic in \(q'\). Furthermore, it can be seen that an applied field induces an energy gap of \(\Delta = g \mu_{\rm B} B\). These effects are illustrated in Figure xx--CrossReference--fig:omegaq_FM--xx.

*The stationary operator method

If the ground state can be described (or approximated) in simple terms, it is in general possible to describe the lowest excited states by exact quantum mechanical calculations. For the ferromagnet, this is particularly simple. Let us first define the lowering operator in reciprocal space as

\begin{equation} \label{dummy1791072476} S_{\mathbf q}^- = \dfrac{1}{\sqrt{N}} \displaystyle\sum_j \exp(i {\mathbf q}\cdot {\mathbf r}_j) s_j^- . \end{equation}

It can be shown by direct application of the Hamiltonian form that the state

\begin{equation} \label{dummy1472779595} |{\mathbf q}\rangle = S_{\mathbf q}^- |0\rangle \end{equation}

is an eigenstate of the Hamiltonian,

\begin{equation} \label{dummy1594344458} \hat{H}|{\mathbf q}\rangle = (E_0 + \hbar\omega_{\mathbf q})|{\mathbf q}\rangle. \end{equation}

Reassuringly, for zero temperature \((\langle S^z \rangle = S)\) this gives the same spin wave frequency, \(\omega_{\bf q}\), as the equation-of-motion method in equation \eqref{eq:fm_omegaq}.

*Spin waves in an antiferromagnet

We now proceed to the much more difficult task of calculating the spin waves in a two-sublattice antiferromagnet (AFM), as first described by P. W. Anderson. We will use the equation-of-motion method as described in the section The equation-of-motion method.

The anisotropic Hamiltonian

We assume that the ground state for the AFM is the classical Neél state, where the spins are fully aligned along the positive \(z\)-axis (sublattice m) or along the negative \(z\)-axis (sublattice n), as described in the section Magnetism in materials on the Magnetic neutron scattering page.

We write the magnetic Hamiltonian as an isotropic Heisenberg term, an applied field, \(B\), and an effective anisotropy field, \(B_{\rm A}\), originating from the surrounding crystal. Both fields are parallel to the \(z\)-axis.

\begin{align}\label{dummy1663920410} \hat{H} &= \displaystyle\sum_{\mathbf m,r} J({\mathbf r}) {\mathbf s}_{\mathbf m}\cdot{\mathbf s}_{\mathbf m+r} + \displaystyle\sum_{\mathbf n,r} J({\mathbf r}) {\mathbf s}_{\mathbf n}\cdot{\mathbf s}_{\mathbf n+r} \\ &\quad+ \displaystyle\sum_{\mathbf m,R} J_1({\mathbf R}) {\mathbf s}_{\mathbf m}\cdot{\mathbf s}_{\mathbf m+R} + \displaystyle\sum_{\mathbf n,R} J_1({\mathbf R}) {\mathbf s}_{\mathbf n}\cdot{\mathbf s}_{\mathbf n+R} \nonumber\\ &\quad- g \mu_{\rm B} \left(B+B_{\rm A}\right) \displaystyle\sum_{\mathbf m} s_{\mathbf m}^z - g \mu_{\rm B} \left(B-B_{\rm A}\right) \displaystyle\sum_{\mathbf n} s_{\mathbf n}^z ,\nonumber \end{align}

where \({\bf r}\) connects spins on the same sublattice with coupling strength \(J({\bf r})\), and \({\bf R}\) connects spins on opposite sublattices with coupling \(J_1({\bf R})\). Note that \({\bf B}_{\rm A}\) points in opposite directions in the two sublattices, because the anisotropy is uniaxial, i.e. it aims the spin to be parallel to a certain direction, but does not have preference to one of the two directions.

*The equations of motion

We now simplify the description by transforming the spin operators on the "down" sublattice (n) by

\begin{equation}\label{dummy1097282024} s_{\rm n}^x = t_{\rm n}^x , \quad s_{\rm n}^y = - t_{\rm n}^y , \quad s_{\rm n}^z = - t_{\rm n}^z . \end{equation}

As a start, we recall the commutation relations (\ref{eq:commute_pm}): \begin{align} \label{eq:commute_pm} [ s_j^+ , s_{j'}^- ] &= 2 s_j^z \delta_{j,j'} \approx 2 \langle S \rangle \delta_{j,j'} , \\ [ s_j^z , s_{j'}^\pm ] &= \pm s_j^\pm \delta_{j,j'} , \\ [ s_j^x , s_{j'}^y ] &= i s_j^z \delta_{j,j'} \; {\rm and\, cyclic\, permutations.} \end{align}

where \(\langle S \rangle\) is here the ordered moment of the sublattice. This important approximation is of course valid only when the Neél state is a good approximation to the actual state of the system. We now calculate the equations of motion, keeping only terms linear in the operators:

\begin{align} i \hbar \dfrac{d}{dt}s_{\mathbf m}^+ &= \left[ s_{\mathbf m} , \hat{H} \right] \label{dummy390816099}\\ &\stackrel{\rm (1st\, order)}{=} 2 S \displaystyle\sum_{\mathbf r} J({\mathbf r}) (s_{\mathbf m}^+ + t_{\mathbf m+r}^-) \nonumber\\ &\quad+ 2 S \displaystyle\sum_{\mathbf R} J_1({\mathbf R}) (-s_{\mathbf m}^+ + s_{\mathbf m+R}^+) + g \mu_{\rm B} (B+B_{\rm A}) s_{\mathbf m}^+ . \nonumber\\ i \hbar \dfrac{d}{dt}t_{\mathbf m}^+ &= \left[ t_{\mathbf m} , \hat{H} \right] \label{dummy270224811}\\ &\stackrel{\rm (1st\, order)}{=} - 2 S \displaystyle\sum_{\mathbf r} J({\mathbf r}) (t_{\mathbf m}^+ + s_{\mathbf m+r}^-) \nonumber\\ &\quad- 2 S \displaystyle\sum_{\mathbf R} J_1({\mathbf R}) (-t_{\mathbf m}^+ + t_{\mathbf m+R}^+) + g \mu_{\rm B} (B+B_{\rm A}) t_{\mathbf m}^+ . \nonumber \end{align}

This calculation differs from that of the ferromagnet, in Inelastic magnetic scattering, by the complication that the equations of motion for the \(s^+\) operators involve the \(t^-\) operators and vice versa. Also the \(t^+\) operators and the \(s^-\) operators are mixed (not shown). Below, we show how to deal with this complication.

*Fourier and Boguliobov transformations

We now solve these two coupled equations of motion, describing the spins in the two sublattices. Firstly, Fourier transformation gives

\begin{align} i \hbar \dfrac{d}{dt} S_{\mathbf q}^+ &= \left[ 2 S (J(0)-J_1(0)+J_1({\mathbf q})) + g \mu_{\rm B} (B+B_{\rm A})\right] S_{\mathbf q}^+ \label{dummy567299487}\\ &\quad + 2 S J({\mathbf q}) T_{\mathbf q}^- , \nonumber\\ i \hbar \dfrac{d}{dt} T_{\mathbf q}^- &= \left[ 2 S (-J(0)+J_1(0)-J_1({\mathbf q})) + g \mu_{\rm B} (B+B_{\rm A})\right] T_{\mathbf q}^- \label{dummy635732318}\\ &\quad - 2 S J({\mathbf q}) S_{\mathbf q}^+ .\nonumber \end{align}

Secondly, we perform a Boguliobov transformation, which is a linear transformation of two operators:

\begin{align} S_{\mathbf q}^+ &= u_{\mathbf q} \hat\alpha_{\mathbf q} + v_{\mathbf q} \hat\beta_{\mathbf q}^\dagger , \label{dummy36441344}\\ T_{\mathbf q}^- &= u_{\mathbf q} \hat\beta_{\mathbf q}^\dagger + v_{\mathbf q} \hat\alpha_{\mathbf q} .\label{dummy1987763760} \end{align}

The idea is that by a correct choice of the real coefficients \(u_{\bf q}\) and \(v_{\bf q}\), the new Boson operators, \(\hat\alpha_{\bf q}\) and \(\hat\beta_{\bf q}\) will perform independent oscillations. To determine the coefficients \(u_{\bf q}\) and \(v_{\bf q}\), we seek a number of boundary conditions. Firstly, Boson operators must fulfill the commutation relations:

\begin{equation} \label{dummy1071622307} [\hat\alpha_{\mathbf q}^\dagger,\hat\alpha_{\mathbf q'}] =\delta_{\mathbf q,q'} , \quad [\hat\beta_{\mathbf q}^\dagger,\hat\beta_{\mathbf q'}] =\delta_{\mathbf q,q'} . \end{equation}

Secondly, the "old" spin commutator relations must still be valid:

\begin{equation} \label{dummy2118875371} [ s_{\mathbf m}^+, s_{\mathbf m'}^- ] = 2s_{\mathbf m}^z \delta_{\mathbf m,m'} , \quad [ t_{\mathbf n}^+, t_{\mathbf n-}^- ] = 2t_{\mathbf n}^z \delta_{\mathbf n,n'} . \end{equation}

The resulting calculations are straightforward but tedious, and we will not list them here. We will proceed directly to the result for the energies.

*The spin wave energy

Figure 2: Dispersion relation for spin waves in a 1D nearest neighbour AFM in zero applied field. The anisotropy is given by \(g \mu_{\rm B} B_{\rm A} = 10^{-2} \times 2SJ(0)\), as defined in equation \eqref{eq:omega_zerofield}.

The antiferromagnetic spin waves have two branches, and their energies are

\begin{align} \hbar \omega_0({\mathbf q}) &= + g \mu_{\rm B} B + \Omega({\mathbf q}) , \label{eq:omega_nonzerofield}\\ \hbar \omega_1({\mathbf q}) &= - g \mu_{\rm B} B + \Omega({\mathbf q}) ,\label{dummy196633099} \end{align}

where the zero field energy is

\begin{equation}\label{eq:omega_zerofield} \Omega({\mathbf q}) = \sqrt{\left[2SJ({\mathbf 0})-2S \left(J_1({\mathbf 0})-J_1({\mathbf q}) \right) + g \mu_{\rm B} B_{\rm A} \right]^2 - (2S)^2 J({\mathbf q})^2} . \end{equation}

The coefficients for the Boguliobov transformation are given by

\begin{align} u_{\mathbf q}^2 - v_{\mathbf q}^2 &= 2 S , \label{dummy1668449917}\\ u_{\mathbf q}^2 &= 2S \; \dfrac{\Omega({\mathbf q})+2SJ({\mathbf 0})-2S[J_1({\mathbf 0})-J_1({\mathbf q})]+g \mu_{\rm B}B_{\rm A}} {2\Omega({\mathbf q})} ,\label{dummy210595604} \end{align}

leading to a few helpful expressions we will need later to determine the scattering cross section:

\begin{align} u_{\mathbf q}^2 + v_{\mathbf q}^2 &= 2S \; \dfrac{2SJ({\mathbf 0})-2S[J_1({\mathbf 0})-J_1({\mathbf q})]+g \mu_{\rm B}B_{\rm A}} {\Omega({\mathbf q})} , \label{dummy52712440}\\ u_{\mathbf q} v_{\mathbf q} &= \dfrac{-(2S)^2 J({\mathbf q})}{2 \Omega({\mathbf q})} .\label{dummy228124291} \end{align}

Example

It is worth contemplating equation \eqref{eq:omega_nonzerofield}. Let us consider the simple case of nearest neighbour couplings only, \(J_1({\bf q})=0\), zero applied field, \(B=0\), and a zero anisotropy, \(B_{\rm A}=0\). Then, the spin wave energy becomes

\begin{equation} \label{dummy875378601} \hbar \omega_{\rm simple} = 2 S \sqrt{J(0)^2-J({\mathbf q})^2} . \end{equation}

For small values of \(|{\bf q}|\), this can be approximated by a linear dispersion of the form \(\hbar\omega = c |{\bf q}|\), where \(c\) is called the spin wave velocity. With a small anisotropy, the spin wave energy becomes

\begin{equation} \hbar \omega_{\rm simple, anis} = 2 S \sqrt{(J(0)+g\mu_{\rm B}B_{\rm A}/(2S) )^2-J({\bf q}')^2} , \end{equation}

as illustrated in Figure xx--CrossReference--fig:omegaq_AFM--xx.

*Quantum magnetism

Solving the Hamiltonian quantum mechanically is an easy task for a ferromagnet. Let us first rewrite the spin product term:

\begin{equation} {\bf S}_i \cdot {\bf S}_j = s_i^z s_j^ z + \frac{1}{2} (s_i^+ s_j^- + s_i^- s_j^+) , \end{equation}

where the spin raising and lowering operators are defined as

\begin{align} \label{eq:spinraising} s_{j}^+ &= s_{j}^x + i s_{j}^y , \\ s_{j}^- &= s_{j}^x - i s_{j}^y . \end{align}

The raising/lowering operators have the effect of changing \(m\), the eigenvalue of \(s^z\) for one spin, by one unit:

\begin{align} s_{j}^+ |S,m\rangle_j &= \sqrt{(S-m)(S+m+1)} |m+1\rangle_j , \\ s_{j}^- |S,m\rangle_j &= \sqrt{(S-m+1)(S+m)} |m-1\rangle_j . \end{align}

For completeness, we mention the commutator relations

\begin{align} \label{eq:commute_pm} [ s_j^+ , s_{j'}^- ] &= 2 s_j^z \delta_{j,j'} , \\{} [ s_j^z , s_{j'}^\pm ] &= \pm s_j^\pm \delta_{j,j'} , \\{} [ s_j^x , s_{j'}^y ] &= i s_j^z \delta_{j,j'} \; {\rm and\, cyclic\, permutations.} \end{align}

We can see that when using the Hamiltonian from the Elastic magnetic scattering page on the ferromagnetic state, \(|{\rm FM}\rangle = |\uparrow \uparrow \uparrow \uparrow \rangle\), all spin terms containing \(s^+\) vanish, leaving only the diagonal \(s^z_i s^z_j\) terms.

\begin{equation} H | {\rm FM} \rangle = - \sum_{\langle ij \rangle} J_{ij} S^2 | {\rm FM} \rangle . \end{equation}

The ferromagnetic state is thus an eigenstate of the Hamiltonian, and is also the ground state. However, for antiferromagnetic couplings, the Hamiltonian is a very difficult many-body problem. This can be seen from just considering the two-spin model for spins \(S=1/2\). Here, the classical ground states are \(| {\rm AFM}_1 \rangle = | \uparrow \downarrow \rangle\) and \(| {\rm AFM}_2 \rangle = | \downarrow \uparrow \rangle\). None of these states are eigenstates to the Hamiltonian, since terms in the Hamiltonian of the type \(s_1^+ s_2^-\) mix the states:

\begin{align} H | {\rm AFM}_1 \rangle &= - \frac{J}{4} |{\rm AFM}_1 \rangle + \frac{J}{2} |{\rm AFM}_2 \rangle , \\ H | {\rm AFM}_2 \rangle &= + \frac{J}{2} |{\rm AFM}_1 \rangle - \frac{J}{4} |{\rm AFM}_2 \rangle . \end{align}

The resulting ground state is a spin singlet:

\begin{equation} |{\rm GS}\rangle = \frac{1}{\sqrt{2}} \left({|\rm AFM}_1 \rangle - |{\rm AFM}_2 \rangle \right) , \end{equation}

with a ground state energy of \(-3J/4\), three times lower than (the expectation value of) the energy of the classical solution.

The antiferromagnetic quantum Heisenberg model for macroscopic systems has only been solved in a few simple one-dimensional cases, and only for spin values \(S=1/2\) and \(S=1\). The main conclusion is that there is no magnetic order in the classical sense, and that the ground state is a many-body quantum singlet state with short-range order. Three-dimensional systems are generally well described by a classical state with few modifications, whereas two-dimensional systems may show intermeditate (and highly interesting) behaviour, where the ground state could be described by a mixture of classical long-range order and quantum short-range order.

Magnetic quantum effects are revealed by the strength of the magnetic order, which is often smaller than expected classically (or completely absent). However, the fingerprint of quantum magnetism is found by inelastic magnetic neutron scattering.

Perspectives

Interest in magnetic quantum systems comes from the study of basic many-body quantum mechanics, entanglement of quantum systems with possible applications in quantum computing, in the study of zero-temperature "quantum phase transitions", as models for the string theory description of elementary particles, and as a basis for the study of correlated electron materials like the "heavy-fermion" systems and the high-temperature superconductors.