InTechOpen uses cookies to offer you the best online experience. By continuing to use our site, you agree to our Privacy Policy.

Physics » Optics and Lasers » "Recent Optical and Photonic Technologies", book edited by Ki Young Kim, ISBN 978-953-7619-71-8, Published: January 1, 2010 under CC BY-NC-SA 3.0 license. © The Author(s).

# Nonreciprocal Phenomena on Reflection of Terahertz Radiation off Antiferromagnets

By T. Dumelow, J. A. P. da Costa, F. Lima and E. L. Albuquerque
DOI: 10.5772/6911

Article top

## Overview

Figure 1. Reflection geometry, showing interchange of incident and reflected beams.

Figure 2. Calculated oblique incidence (θ 1 = 45 ) reflectivity spectrum off MnF2 in an external field of B 0 = +0.1 T (solid curves) and B 0 = −0.1 T (dashed curves). Note that in (a) the two curves are coincident. The symbols B and R in (a) represent the bulk and reststrahl frequency regions respectively.

Figure 3. Experimental (solid curves) and theoretical (dashed curves) oblique incidence reflectivity (θ 1 = 45 ) reflectivity spectra off FeF2 at 4.2 K, in the presence of positive and negative external magnetic fields. After Brown et al. (1994).

Figure 4. Calculated reflected amplitude and phase spectra for oblique incidence (θ=45 ) reflection off MnF2 in an external field of B0 = +0.1 T (solid curves) and B0 = −0.1 T (dashed curves). Note that in (a) the two curves are coincident. The symbols B and R represent the bulk and reststrahl frequency regions respectively.

Figure 5. Angle of refraction θ2 defined by power flow direction.

Figure 6. Calculated θ2 values on oblique incidence reflection (θ1 = 45 ) off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), ignoring damping. Both curves are coincident in the bulk regions. The symbols B and R represent the bulk and reststrahl frequency regions (separated by dashed vertical lines) respectively.

Figure 7. Power flow in the various spectral regions for oblique incidence radiation.

Figure 8. Calculated θ2 values on normal incidence reflection off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), ignoring damping. The symbols B and R represent the bulk and reststrahl frequency regions (separated by dashed vertical lines) respectively. After Lima et al. (2009).

Figure 9. Power flow in the various spectral regions for normal incidence radiation.

Figure 10. Calculated θ2 values on reflection off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), assuming a damping parameter of Γ = 0.0007 cm−1.

Figure 11. Comparison of the normal incidence Goos-Hänchen shift predicted for external reflection off an antiferromagnet with the conventional Goos Hänchen shift for total internal reflection.

Figure 12. Calculated normal incidence reflectivity spectrum and 1ateral shift for reflection off MnF2 in an external field of B0 = +0.1 T. In (a) the solid curve is calculated without damping and the dot-dash curve with Γ = 0.0007 cm−1. In (b) the solid curves are for an external field of B0 = ±0.1 T and the dashed curves for an external field of B0 = −0.1 T. The vertical lines separate the bulk and reststrahl regions. The marked frequencies are A: 9.0769 cm−1; B: 9.0994 cm−1; C: 9.0705 cm−1. After Lima et al. (2009).

Figure 13. Intensity profiles, along the sample surface, in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T. Solid curve: incident beam; dotted curve: reflected beam ignoring damping; dashed curve: reflected beam with Γ = 0.0007 cm−1; vertical solid line: centre of incident beam; vertical dashed line: centre of reflected beam predicted by Equation 44. Note that, in (b), the reflected intensities with and without damping are almost identical, and the curves cannot be separated. After Lima et al. (2009).

Figure 14. Profiles of the amplitudes of the incident (Ei(x,y)), transmitted (Et(x,y)), and reflected (Er(x,y)) fields in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T, ignoring damping. After Lima et al. (2009).

Figure 15. Profiles of the amplitudes of the incident (Ei(x,y)), transmitted (Et(x,y)), and reflected (Er(x,y)) fields in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T, for Γ = 0.0007 cm−1.

Figure 16. Model of lateral displacement of reflected beam in the reststrahl regions in terms of power flow along the antiferromagnet surface.

Figure 17. Overall power intensity and flow for reflection off MnF2 in an external magnetic field of 0.1 T in the upper reststrahl region (frequency A). After Lima et al. (2009).

Figure 18. Overall power intensity and flow for reflection off MnF2 in an external magnetic field of 0.1 T in the bulk regions. Note that the arrows in the bottom half of (c) are simply intended to show the direction of power flow and are not a continuation of the flux lines in the upper part of the figure. After Lima et al. (2009).

Figure 19. Use of a dielectric layer for investigating nonreciprocal phase on reflection off an antiferromagnet.

Figure 20. Calculated oblique incidence reflectivity spectrum off a Si/MnF2 structure of the type shown in Figure 19 in an external field of B0 = +0.1 T (solid curves) and B0 = −0.1 T (dashed curves).

# Nonreciprocal Phenomena on Reflection of Terahertz Radiation off Antiferromagnets

T. Dumelow1, J. A. P. da Costa1, F. Lima2 and E. L. Albuquerque2

## 1. Introduction

There are a number of ways that reciprocity principles in optics may be affected by the presence of a static magnetic field (Potton, 2004). A classic example is Faraday rotation in which a plane polarised electromagnetic beam propagating through a suitable medium is rotated in the presence of a static magnetic field along the direction of propagation. The handedness of this rotation depends on the propagation direction, a nonreciprocal effect usefully applied to the construction of optical isolators (Dötsch et al., 2005). Nonreciprocal effects of this type are closely related to the idea that magnetic fields break time reversal symmetry. Similar nonreciprocal phenomena can occur, in various guises, on reflection off a semi-infinite sample. We discuss such behaviour in the present chapter, in the context of reflection off antiferromagnetic materials. In contrast to nonreciprocal phenomena based on the Faraday effect, our interest is in the Voigt geometry, in which the static magnetic field is perpendicular to the direction of propagation. We consider the well established phenomena of nonreciprocity in the intensity and phase of oblique incidence radiation, but concentrate mainly on recent developments on nonreciprocal power flow and finite beam effects.

We restrict discussion to the simple two dimensional geometry shown in Figure 1. Radiation is reflected, in the xy plane, off a semi-infinite sample, isotropic in this plane, in the presence of a static magnetic field B0 along z (into the page). Note that, in this configuration, we do not have to worry about polarisation effects, since there is no mixing between s-polarised (electromagnetic E field component along z) and p-polarised (electromagnetic H field component along z) radiation.

Now compare Figure 1(a) to Figure 1(b), in which the sign of the incident angle has been reversed. In the absence of the magnetic field (B0 = 0), we can consider Figure 1(b) as the mirror reflection of Figure 1(a) through the yz plane, so we expect no change in the reflection behaviour. In terms of the incident and reflected beam signals, this is a trivial example of the Helmholtz reciprocity principle, which, in the present context, can be interpreted as saying that, in the absence of magnetic fields, an interchange of source and detector should not affect the signal received by the detector (Born & Wolf, 1980). When B0 is nonzero, however, the mirror reflection of Figure 1(a) through the yz plane no longer leads to Figure 1(b), as one might expect. The essential point here is that the static magnetic field B0 is an axial vector, and a mirror symmetry operation through the yz plane would therefore involve reversing the direction of this field, so that it would come out of the page (Scott & Mills, 1977). In fact

#### Figure 1.

Reflection geometry, showing interchange of incident and reflected beams.

there is no symmetry operation that leads us from Figure 1(a) to Figure 1(b), and the two figures are not equivalent. Nonreciprocal behaviour is thus, in principle, possible. Whether or not it occurs in practice, however, depends on the material properties of the sample.

In the present chapter we consider nonreciprocity associated with reflection off a simple uniaxial antiferromagnet. In this case the static field represented by B0 in Figure 1 is an external field, since an antiferromagnet has no intrinsic macroscopic magnetic field. We consider a geometry in which the anisotropy associated with the spin directions, along with the external field B0, is perpendicular to the plane of incidence. This is equivalent to putting the anisotropy along z in Figure 1, thus leaving the antiferromagnet isotropic in the xy plane. The electric component of the electromagnetic field is along z and the magnetic component is in the xy plane (s-polarisation).

In considering nonreciprocity in the intensity and phase of the reflected beam, it is sufficient to simply consider the effect of interchanging the incident and reflected beams (i.e. reversing the sign of θ 1). However, we note that a rotation of Figure 1(b) around the y axis brings us back to Figure 1(a), but with the field direction reversed. It is therefore possible to consider nonreciprocity in terms of a change in optical behaviour when the external field direction is reversed. This turns out to be more convenient when considering nonreciprocal effects inside the antiferromagnet and finite beam effects. It is notable that some of the new phenomena under investigation in this chapter occur at normal incidence, so such a test is simpler to visualise in such cases than a test based on the configurations of Figure 1. Thus our general test for nonreciprocity will be to see what happens when we reverse the sign of B 0.

## 2. Antiferromagnet permeability

The crucial parameter that determines the nonreciprocal optical properties of antiferromagnets is the magnetic permeability in region of the magnon (or spin wave) frequencies (Mills & Burstein, 1974), which typically lie in the terahertz range. We think of an antiferromagnet as two interpenetrating sublattices having opposite spin directions. Waves consisting of spins precessing in opposite directions in the two sublattices are then possible, and magnons of this type can interact with electromagnetic radiation. Their resonant frequencies are linked not only to the anisotropy field B A that tends to align the spins along a preferred axis (the z axis in our coordinate system), but also to the interaction between the spins in the two sublattices. In the long wavelength limit (applicable to terahertz frequencies), and in the absence of any external field, the resonance frequency is given by

 ωr = γ(2BABE+BA2)1/2 (1)

Here B E is the exchange field representing the interaction between the spins of the opposing sublattices and γ is the gyromagnetic ratio. In the presence of an electromagnetic field whose H component lies in the xy plane, the induced magnetisation follows the direction of this field component, since the spins in the two sublattices precess in opposite directions with equal amplitudes. The permeability tensor μ is thus diagonal and of the form

 μ = (μ000μ0001) (2)

The scalar quantity μ is given, at frequency ω, by

 μ =1+2μ0γ2BAMSωr2−ω2+iωΓ (3)

where M S is the sublattice magnetisation and Γ is a damping parameter.

In this study, we are interested in propagation of electromagnetic waves (strictly speaking polariton waves, since the waves include a contribution from the precessing spins in addition to that of the electromagnetic radiation) within the xy plane. We consider the electromagnetic E field component to be directed along z with the corresponding H field component in the xy plane. In this case, for plane waves of the form

 E ( xyt ) = E0exp[i(kxx+k1yy−ωt)] (4)
 H ( xyt ) = H0exp[i(kxx+k1yy−ωt)] (5)

the polaritons follow the familiar dispersion relation

kx2+ky2 =
 μk02 (6)

where k x and k y are wavevector components and is the dielectric constant of the medium. k 0 is the modulus of the free space wavevector, given by

 k0 = ωc (7)

In the presence of an external field B0 along the anisotropy axis, the two sublattices are no longer equivalent. This leads to two effects. Firstly, there are now two resonances instead of one and, secondly, the permeability tensor is no longer diagonal, but gyromagnetic. It thus takes the form (Mills & Burstein, 1974):

 μ = (μ1iμ20−iμ2μ10001) (8)

where

 μ1 =1+μ0γ2BAMS ( Y++Y− ) (9)
 μ2 = μ0γ2BAMS ( Y+−Y− ) (10)
with
 Y± = 1ωr2−(ω±γB0+iΓ)2 (11)

The diagonal elements μ 1 do not depend on the sign of the applied field B0, but μ 2 changes sign when B0 is reversed. This is the basis of the nonreciprocal effects discussed in this chapter. The polariton dispersion relation (Equation 6) is now replaced by

kx2+ky2 =
 μvk02 (12)

where μ v is known as the Voigt permeability, and is given by

 μv = μ12−μ22μ1 (13)

It is straightforward to see that μ v does not depend on the sign of the external field B 0, so the polariton dispersion relation (Equation 12) is similarly unaffected. Thus polariton dispersion corresponding to propagation through an antiferromagnet (as bulk polaritons) is, in the present geometry, totally reciprocal. Nonreciprocal effects only occur in the presence of a surface, as in the case of reflection off an antiferromagnet (Camley, 1987).

## 3. Nonreciprocity in reflection of plane waves

### 3.1. Reflected intensity

As discussed in the introduction, we can regard reflectivity R as nonreciprocal if there is a change in reflected intensity when the incident and reflected beams are interchanged, i.e. R(θ 1) ≠R(−θ 1) where θ 1 is the angle of incidence (see Figure 1), or, equivalently, when the applied field B0 is reversed, i.e. R(B0) ≠R(−B0). The possibility of nonreciprocal reflectivity in the present geometry was first analysed using thermodynamic arguments (Remer et al., 1984, Camley, 1987, Stamps et al., 1991). This analysis shows that reflectivity should be reciprocal in the absence of absorption, but that it need not be in the presence of absorption. Here we demonstrate the same result explicitly in the case of reflection off a uniaxial antiferromagnet, using the arguments outlined by Abraha & Tilley (1996) and Dumelow et al. (1998).

We are interested in reflection from vacuum in s polarisation. The complex reflection coefficient r in this case can be easily worked out from the field continuity conditions at the vacuum/antiferromagnet interface. Written in terms of the E field component of the electromagnetic radiation, the complex reflection coefficient is given by

 r = k1yμv−k2y−ikx ( μ2 / μ1 ) k1yμv+k2y+ikx ( μ2 / μ1 ) (14)

with a corresponding transmission coefficient

 t =1+r = 2k1yμv ( μ2 / μ1 ) k1yμv+k2y+ikx ( μ2 / μ1 ) (15)

k x is the in-plane component of the wavevector, which is continuous in both media and determined by the angle of incidence θ 1:

 kx = k0sinθ1 (16)

k 1y and k 2y are the normal components of the wavevector in vacuum and the antiferromagnet respectively, and are given by

 k1y = k02−kx2 (17)
and
 k2y=∈μk02-kx2 (18)

Since k x(θ 1) = −k x(−θ 1) and μ 2(B0) = −μ 2(−B0), the effect of either changing the sign of θ 1 or changing the sign of B 0 is to change the sign of the term ik x(μ 2/μ 1) in both the numerator and the denominator of Equation 14, all other terms in this equation being unaffected.

Let us first consider how this sign change affects the complex reflection coefficient r in the case of zero absorption (Γ = 0). In this case, all the parameters in Equation 14 are real, except for k 2y, which may be either real or imaginary depending on the frequency. When k 2y is real, there is propagation of radiation, as bulk polaritons, into the interior of the sample. When it is imaginary the field within the antiferromagnet is evanescent, decaying away from the interface. We refer to frequencies corresponding to k 2y real as bulk region frequencies and those corresponding to k 2y imaginary as reststrahl region frequencies.

At bulk region frequencies (k 2y real), separation of Equation 14 into real and imaginary parts leads to

 r*(B0)=r(-B0) (19)

The overall reflectivity is given by

 R=rr* (20)

and is therefore reciprocal.

At reststrahl region frequencies (k 2y imaginary), the numerator of Equation 14 is the complex conjugate of the denominator. Thus one can see from Equation 20 that R must be equal to 1, regardless of the sign of B 0, so once again the reflectivity is reciprocal. The result R = 1 is of course what one should expect, since when k 2y is imaginary there is no propagation into the sample, leading to total reflection.

Simulated oblique incidence (θ 1 = 45 ) reflectivity spectra off MnF2 at 4.2 K, in an external magnetic field B0 of magnitude 0.1 T, are shown in Figure 2(a), in which zero damping is assumed. The frequency scale is expressed in terms of wavenumbers ω/2πc, and the MnF2 parameters used in the calculation are (Dumelow & Oliveros, 1997) =5.5, M S =6.0×105 A/m, B A = 0.787 T, B E = 53.0 T and γ = 0.975 cm−1/T, corresponding to ω r = 8.94 cm−1. The curves for B 0 = +0.1 T and B 0 = −0.1 T are coincident at all frequencies, confirming that the reflectivity is reciprocal.

### Figure 2.

Calculated oblique incidence (θ 1 = 45 ) reflectivity spectrum off MnF2 in an external field of B 0 = +0.1 T (solid curves) and B 0 = −0.1 T (dashed curves). Note that in (a) the two curves are coincident. The symbols B and R in (a) represent the bulk and reststrahl frequency regions respectively.

### Figure 3.

Experimental (solid curves) and theoretical (dashed curves) oblique incidence reflectivity (θ 1 = 45 ) reflectivity spectra off FeF2 at 4.2 K, in the presence of positive and negative external magnetic fields. After Brown et al. (1994).

When the damping is nonzero, the above symmetry arguments do not apply, and the reflectivity R is, in general, nonreciprocal. This can be seen from Figure 2(b), in which damping has been included in the calculation, using the experimental value of Γ = 0.0007 cm−1. In this figure, the curves for B 0 = +0.1 T and B 0 = −0.1 T are not coincident. The difference in this case is quite small, since the damping in MnF2 is small. Remer et al. (1986) was able to observe this small difference experimentally using a field scan at fixed frequency, but found a larger nonreciprocity at higher temperature, corresponding to larger damping. 4.2 K frequency scans of the type shown in Figure 2 have been performed on FeF2 (Brown et al., 1994), which has a considerably higher damping parameter than MnF2, using far infrared Fourier transform spectroscopy (Brown et al., 1998). In this case the nonreciprocity in the reflectivity is quite clear, as seen in Figure 3.

### 3.2. Reflected phase

The complex reflection coefficient r given by Equation 14 is commonly expressed in terms of a reflection amplitude ρ r and phase ϕr:

 r = ρrexp(iϕr) (21)

where

 ρr = R1/2 (22)
and
 ϕr = tan−1[Im ( r ) Re ( r ) ] (23)

Although thermodynamic arguments show that, in the absence of damping, the reflected intensity R, and hence the amplitude ρr, should be reciprocal (Remer et al., 1984), such arguments cannot be applied to the reflected phase ϕr. A detailed discussion of nonreciprocity in the reflected phase on reflection off antiferromagnets is given in Dumelow et al. (1998). Here we summarise the main results.

We consider first the case of zero damping (Γ = 0). In the bulk regions, Equation 19 should apply. Thus, since the phase is given by Equation 23, we can see quite straightforwardly that

 ϕr ( B0 ) = −ϕr ( −B0 ) +2πm (24)

where m is an arbitrary integer. We include the term 2πm since it is convenient to plot the phase outside the range −π < ϕr< π.

Equation 24 shows that, in the bulk regions, the reflected phase is nonreciprocal even in the absence of damping. This is also the case in the reststrahl regions, but the phase does not follow a simple symmetry relation of the type given by this equation.

The amplitude and phase for reflection off MnF2 in the absence of damping are shown in Figures 4(a) and 4(c) respectively. The conditions are the same as those used in Figure 2. Note that we have shown the phase as varying within the range π to 3π in order to show that it changes continuously with frequency. The amplitude is reciprocal, in agreement with Figure 2, but nonreciprocity in the reflected phase is quite marked in both the bulk and reststrahl regions, obeying Equation 24 in the bulk regions.

Figures 4(b) and 4(d) show reflection amplitude and phase respectively in the presence of damping. In line with the reflectivity results in the previous subsection, the reflection amplitude now shows slight nonreciprocity. The phase shows the same type of nonreciprocity as seen without damping, although the bulk region symmetry arguments of Equation 24 no longer apply.

#### Figure 4.

Calculated reflected amplitude and phase spectra for oblique incidence (θ=45 ) reflection off MnF2 in an external field of B0 = +0.1 T (solid curves) and B0 = −0.1 T (dashed curves). Note that in (a) the two curves are coincident. The symbols B and R represent the bulk and reststrahl frequency regions respectively.

### 3.3. Power flow

The nonreciprocal phenomena described in the Subsections 3.1 and 3.2 were analysed ten or more years ago, and concern the behaviour of a reflected plane wave. Recently we have started studying nonreciprocal behaviour within the antiferromagnet itself, in particular with respect to the direction of the internal power flow (Lima et al., 2009), represented by the time-averaged Poynting vector (Landau & Lifshitz, 1984),

 〈S〉 =1/2Re(E×H∗) (25)

We consider an angle of refraction in terms of the direction of the time-averaged Poynting vector S2 (which is not necessarily the same as the wavevector direction) in the antiferromagnet, as shown in Figure 5. The angle of refraction θ2, defined in this way, is given by

 tanθ2 = 〈S2x〉〈S2z〉 (26)

In s polarisation the E field is confined along z, so the Poynting vector is most easily represented in terms of the Ez field, making use of the conversion k × E = ωμ0μH. The resulting time averaged Poynting vector has components

### Figure 5.

Angle of refraction θ2 defined by power flow direction.

 〈S2x〉 = |Ez|22ωμ0Re[kx−ik2y ( μ2 / μ1 ) μv] (27)
 〈S2y〉 = |Ez|22ωμ0Re[k2y+ikx ( μ2 / μ1 ) μv] (28)

The direction of power flow can thus be obtained by substitution into Equation 26.

We now investigate the above expressions in order to search for possible nonreciprocity in the power flow direction in the antiferromagnet, taking power flow to be nonreciprocal if

 θ2 ( −B0 ) ≠θ2 ( B0 ) (29)

In order to consider power flow, we restrict ourselves initially to the case where there is no damping in the system (Γ = 0). The calculated values of θ2 for oblique incident reflection off MnF2 in this case are shown in Figure 6.

### Figure 6.

Calculated θ2 values on oblique incidence reflection (θ1 = 45 ) off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), ignoring damping. Both curves are coincident in the bulk regions. The symbols B and R represent the bulk and reststrahl frequency regions (separated by dashed vertical lines) respectively.

In the case of zero damping, as discussed previously, μ1, μ2, and μv are all wholly real. k2y is real in the bulk regions and imaginary in the reststrahl regions.

First we consider power flow for k2y real (i.e. in the bulk regions). Equations 27 and 28 then give

 〈S2x〉 = |Ez|2kx2ωμ0μv (30)
 〈S2y〉 = |Ez|2k2y2ωμ0μv (31)

Thus, since none of the terms in Equations 30 and 31 depend on the sign of B0, we see that the power flow direction is reciprocal (θ2(−B0) = θ2(B0)). We also see from Equations 26, 30 and 31 that

 tanθ2 = kxk2y (32)

so that S2 is parallel to k2. Since kx is continuous and k2y positive it also follows that refraction must be positive, i.e. θ2 always has the same sign as θ1. Overall, therefore, power flow follows the type of behavior shown in Fig. 7(a). This is confirmed by the calculations of the θ2 shown in Figure 6 for the bulk regions. Note that in these regions radiation may, in principle, flow an infinite distance into the antiferromagnet, and is thus unaffected by the sample surface. Our result that power flow is reciprocal in this case is thus consistent with the idea that radiation should display reciprocal behavior in the interior of a sample.

### Figure 7.

Power flow in the various spectral regions for oblique incidence radiation.

We now consider the case when k2y is imaginary. In this case we get

 〈S2x〉 = |Ez|2[kx−ik2y ( μ2 / μ1 ) ]2ωμ0μv (33)
 〈S2y〉 =0 (34)
There is thus no propagation into the antiferromagnet, but energy can travel along the surface in the form of an evanescent wave within the antiferromagnet. In such cases, which are characterized by θ2 = ±90 , the angle θ2 does not describe refraction in the normal sense of the word, since the associated radiation has an evanescent component. Indeed there is no transfer of radiation across the interface, since power is wholly reflected. The direction of power flow may or may not be reciprocal, depending on the relative magnitudes of the two terms in the numerator of Equation 33. However, if the external field is sufficiently large, the second term dominates within both reststrahl regions. In the present example, for B0 = ±0.1 T, this results in θ2 = ±90 in most of the lower reststrahl region (Figure 7(b)) and θ2 = 90 in most of the upper reststrahl region (Figure 7(c)), as shown in the calculation in Figure 6, and nonreciprocity, as defined by Equation 29, occurs. There is, however, a narrow region of reciprocal power flow (in terms of direction if not of magnitude) near the upper frequency limit of each reststrahl region.Possibly one of the most interesting aspects of the power flow is its behaviour in the case of normally incident radiation (Lima et al., 2009). The calculated θ2 values are shown in Figure 8. In the bulk regions, as expected, power flow is normal to the surface, as shown in Figure 9(a). In the reststrahl regions, however, power flow parallel to the sample surface is induced in the same way as occurs at oblique incidence. The direction of this power flow depends on the sign of the external field, as shown in Figures 9(b) and 9(c). This is distinctly counter- intuitive considering that we are modelling plane waves normally incident on an infinite surface, so that there should be no spatial variation of the field along x. Nevertheless, the

### Figure 8.

Calculated θ2 values on normal incidence reflection off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), ignoring damping. The symbols B and R represent the bulk and reststrahl frequency regions (separated by dashed vertical lines) respectively. After Lima et al. (2009).

### Figure 9.

Power flow in the various spectral regions for normal incidence radiation.

very fact that the field is assumed to extend over an infinite plane means that energy flow along the surface does not violate causality.

We now turn to the case where damping is present. Now k2y should, in general, be complex, with the its imaginary part greater near the reststrahl regions. Since imaginary k2y results in nonreciprocal power flow, we should also expect such nonreciprocity in the case of complex k2y. In Fig. 10 we show the power flow directions both for normal and for oblique incidence, assuming Γ = 0.0007cm−1. There is now no distinct division between reststrahl and bulk regions, and nonreciprocal power flow occurs both inside and outside the nominal reststrahl regions. At normal incidence the power flow directions are now no longer restricted to θ2 = 0 and θ2 = ±90 . Since nonzero θ2 implies nonreciprocity, the associated fields must be some extent be bound to the sample surface in all regions for which θ2 ≠ 0 .

### Figure 10.

Calculated θ2 values on reflection off MnF2 in a field of B0 = +0.1 T (solid curve) and B0 = −0.1 T (dashed curve), assuming a damping parameter of Γ = 0.0007 cm−1.

## 4. Reflection of a finite beam

### 4.1.Introduction

The previous section discusses various phenomena associated with plane wave reflection off an antiferromagnet. However, when the plane wave is replaced by a finite beam, we predict additional effects concerning the profile and position of the reflected beam (Lima et al., 2008, 2009). Such effects are expected with either normal or oblique incidence radiation. However, we concentrate on normal incidence effects firstly for simplicity and secondly because they are more unexpected.

We examine reflection of finite beams in two ways. Firstly, we interpret the reflection using an angular spectrum analysis, in which the incident beam is considered as a Fourier sum of plane waves. We then give a power flow interpretation of the predicted effects.

### 4.2. Angular spectrum analysis

Here we summarise the angular spectrum analysis used in describing reflection of a finite beam normally incident on an antiferromagnet (Lima et al., 2008), using a two dimensional model in which the incident beam, centred at x = 0, is considered as an angular spectrum of plane waves propagating in the xy plane. It can thus be represented in the form

 Ei ( xy ) = ∫−k0k0ψ ( kx ) exp[i(kxx+k1yy)]dkx (35)

where ψ(kx) is a distribution function representing the shape of the beam. At the sample surface, which we place at y = 0, the electric field of the incident beam is

 Ei ( x,0) = ∫−k0k0ψ ( kx ) exp(ikxx)dkx (36)

The electric field of the corresponding reflected beam is given, at the surface, by

 Er ( x,0) = ∫−k0k0r ( kx ) ψ ( kx ) exp(ikxx)dkx (37)

Here r(kx) represents the complex reflection coefficient for the relevant plane wave component, and is given by Equation 14. It is convenient to consider this complex coefficient in terms of amplitude ρr(kx) and phase ϕr(kx), in the form of Equation 21, i.e.,

 r ( kx ) = ρr ( kx ) exp[iϕr ( kx ) ] (38)

If the beam is sufficiently wide, there will be a narrow distribution of kx values centred, at normal incidence, around kx = 0. We can thus substitute Equation 38 into Equation 37 and expand ρr(kx) and ϕr(kx) as a Taylor series around kx = 0. If we ignore terms in kx2 and above, this leads to

 Er ( x,0) = r (0) ∫−k0k0ψ ( kx ) exp(ikxX)dkx+r' (0) ∫−k0k0kxψ ( kx ) exp(ikxX)dkx (39)

where

 r (0) = ρr (0) exp[iϕr (0) ] (40)
 r' (0) = dρrdkx | kx =0exp[iϕr (0) ] (41)
 X = x+dϕrdkx | kx =0 (42)

The second term on the right hand side of Equation 39 can, in practice, normally be ignored (Lima et al., 2008). In fact, in the absence of damping, it is identically zero since reciprocity in ρr implies dρr / dkx | kx =0 =0 , and hence r′(0) = 0. Even in the presence of damping, however, the contribution of this term is negligible for a sufficiently narrow kx distribution.

The reflected field is thus given, to a good approximation, by the first term on the right hand side of Equation 39. This term is simply the reflection coefficient r(0) for a normally incident plane wave multiplied by an integral which gives the profile of the field along x. This integral is identical to that of the incident beam (Equation 36) except that x has been replaced by X, given by Equation 42. Thus the shape of the reflected beam is the same as that of the incident beam, but it has been shifted along the surface of the sample by a distance Dr equal to

 Dr = −dϕrdkx | kx =0 (43)

### Figure 11.

Comparison of the normal incidence Goos-Hänchen shift predicted for external reflection off an antiferromagnet with the conventional Goos Hänchen shift for total internal reflection.

as shown in Figure 11(a). In the case of reciprocal reflected phase (i.e. ϕ(−kx) = ϕ(kx)), dϕr / dkx | kx =0 must be zero, so there will be no lateral shift. In the case of reflection off antiferromagnets, however, the reflected phase is, in general, nonreciprocal, as discussed in Section 3.2, and a nonzero displacement of the reflected beam is expected.

Before discussing the lateral shift described by Equation 43 in any detail, we note that the behaviour of the reflected beam is in some ways similar to that of an oblique incidence finite beam which undergoes total internal reflection when passing from an optically denser to a less dense medium. Such a beam also suffers a lateral displacement Dr upon reflection, as shown in Figure 11(b). This displacement was observed experimentally by Goos & Hänchen (1947) and is normally referred to as a Goos-Hänchen shift (Lotsch, 1970). Goos-Hänchen shifts have also been reported in the case of external reflection in specific instances such as reflection off metals (Wild & Giles, 1982, Leung et al., 2007), and the lateral displacement discussed here can be considered as a type of a normal incidence Goos-Hänchen shift. Equation 43 is, indeed, a normal incidence version of the classical expression commonly used to describe Goos-Hänchen shifts (Artmann, 1948), and the angular spectrum analysis used above in deriving the equation is basically the same as that previously used to describe Goos-Hänchen shifts in the case of total internal reflection (Horowitz & Tamir, 1971, McGuirk & Carniglia, 1977).

We now apply Equation 43 to the specific case of reflection off an antiferromagnet. In the absence of damping, the reflection coefficient r can easily be resolved into its real and imaginary parts, allowing explicit evaluation of this equation. This gives

 Dr=2(μ2/μ1)k0(μv-∈) (44)

in both the bulk and the reststrahl regions. Since the sign of μ2 depends on the sign of B0, it is immediately obvious that the direction of the lateral displacement will be reversed if the external field direction is reversed. In the presence of damping, Equation 43 can be evaluated numerically, but the results are found to be almost identical to those with Γ = 0 (Lima et al., 2008).

### Figure 12.

Calculated normal incidence reflectivity spectrum and 1ateral shift for reflection off MnF2 in an external field of B0 = +0.1 T. In (a) the solid curve is calculated without damping and the dot-dash curve with Γ = 0.0007 cm−1. In (b) the solid curves are for an external field of B0 = ±0.1 T and the dashed curves for an external field of B0 = −0.1 T. The vertical lines separate the bulk and reststrahl regions. The marked frequencies are A: 9.0769 cm−1; B: 9.0994 cm−1; C: 9.0705 cm−1. After Lima et al. (2009).

The calculated values of Dr, ignoring damping, are compared with the reflectivity spectrum at upper reststrahl region frequencies in Figure 12. A lateral shift is predicted in both the bulk and the reststrahl regions. Similar shifts are predicted around the lower reststrahl region, but of opposite sign (Lima et al., 2008). There is a divergence in Dr at the reflectivity minimum just below the reststrahl region. At this frequency, the assumptions made in deriving Equation Dr clearly do not apply.

In order ot verify the shifts shown in Figure 12, reflection of a particular incident beam profile may be modelled. For simplicity, we have considered a gaussian beam whose focal plane is the surface of the sample (Lima et al., 2008, 2009). The incident beam can thus be represented by Equation (35) with (Horowitz & Tamir, 1971)

 ψ ( kx ) = g2πexp(−g2kx24) (45)

where 2g represents the beam width at the focal plane. At the sample surface the incident beam profile is thus represented by Equation 36 and the reflected beam profile by Equation 37. The integrals in these two equations can be evaluated numerically, and the profiles of the corresponding E fields thus obtained.

The incident and reflected beam intensities can, in general, be well represented by |E|2. In Figure 13 we show the resulting intensity profiles along x for the three frequencies marked in Figure 12. It is seen that, although the modelled incident beam is very narrow (g = λ, the free space wavelength), all the reflected beams are displaced along the surface in excellent agreement with Equation 44. It is also observed that damping does not noticeably affect this displacement. Explicit simulations have also confirmed the prediction of Equation 44 that beam displacement is, to a very good approximation, independent of beam width (Lima et al., 2008).

It is useful not only to consider the incident and reflected fields at the surface, but also the overall E field distribution in the xy plane. We conveniently consider the electric fields in terms of an incident field Ei(x,y) and a reflected field Er(x,y) in the region x < 0 (vacuum) and a transmitted field Et(x,y) in the region x > 0 (antiferromagnet). Ei(x,y) is given by Equation 35, with Er(x,y) and Et(x,y) given by

### Figure 13.

Intensity profiles, along the sample surface, in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T. Solid curve: incident beam; dotted curve: reflected beam ignoring damping; dashed curve: reflected beam with Γ = 0.0007 cm−1; vertical solid line: centre of incident beam; vertical dashed line: centre of reflected beam predicted by Equation 44. Note that, in (b), the reflected intensities with and without damping are almost identical, and the curves cannot be separated. After Lima et al. (2009).

 Er ( xy ) = ∫−k0k0r ( kx ) ψ ( kx ) exp[i(kxx−k1yy)]dkx (46)
and
 Et ( xy ) = ∫−k0k0t ( kx ) ψ ( kx ) exp[i(kxx+k2yy)]dkx (47)

respectively. The resulting profiles at the three frequencies indicated in Figure 12 are shown in Figures 14 and 15. Figure 14 shows the profile in the absence of damping and Figure 15 shows the profile when damping is included. The left hand panels show the incident and transmitted fields, whilst the right hand panels show the reflected field.

In all cases the reflected fields are displaced along x in accordance with Figures 12 and 13. In addition, we see that the transmitted field is also displaced. This result can be anticipated from the fact that the phase ϕt of the complex transmission coefficient t is nonreciprocal. An analysis equivalent to that used in determining Dr then gives a lateral displacement Dt of the transmitted field profile given by

 Dt = −dϕtdkx | kx =0 (48)

### Figure 14.

Profiles of the amplitudes of the incident (Ei(x,y)), transmitted (Et(x,y)), and reflected (Er(x,y)) fields in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T, ignoring damping. After Lima et al. (2009).

### Figure 15.

Profiles of the amplitudes of the incident (Ei(x,y)), transmitted (Et(x,y)), and reflected (Er(x,y)) fields in the case of reflection of a normally incident gaussian beam off MnF2 in an external magnetic field of 0.1 T, for Γ = 0.0007 cm−1.

It is noticeable that, in the case of the bulk frequency B, near the top of the reststrahl band, there is a much larger displacement in the transmitted field than in that of the reflected beam. In the absence of damping, the field decays away from the interface in the reststrahl region (Figure 14(a)), but propagates into the sample, perpendicular to the surface, in the bulk regions (Figures 14(b) and 14(c)) . When damping is included, the reststrahl region behaviour is largely unchanged, but there may now be significant decay in the bulk regions. At frequency B, just above the top of the reststrahl region, this decay is fairly small, but at frequency C, near the bottom of the reststrahl region, it is similar to that in the reststrahl region itself.

### 4.3. Power flow analysis

#### 4.3.1. Reststrahl region

In the above analysis we have considered the displacement of the reflected beam as an interference effect. It is also useful, however, to consider this effect in terms of power flow. In the reststrahl regions, in the absence of absorption, energy conservation principles can be used in a straightforward way to analyse the lateral shift (Lima et al., 2009). The analysis is similar to that used by Renard (1964) in the case of the conventional Goos-Hänchen shift for total internal reflection.

We make use of the result of subsection 3.3 that, in the reststrahl regions and in the absence of damping, a normally incident plane wave reflected off an antiferromagnet induces power flow parallel to the surface within the antiferromagnet, as shown in Figures 9(b) and 9(c). We consider this to be the behaviour in the centre portion of a wide finite incident beam, such as that represented in Figure 16. The central portion of this incident beam lies between x2 and x3, and it gradually decays away to zero between x2 and x1 and between x3 and x4. The internal energy flux associated with plane wave reflection in the central portion is represented by P2. Energy conservation therefore requires that there is a net flux P1 entering the antiferromagnet near one edge of the beam and a net flux leaving near the other edge. This is equivalent to a lateral shift Dr of the reflected beam with respect to the incident beam. Thus P1 enters in the region between x1 and x2 + Dr and P3 leaves in the region between x3 and x4 + Dr. We consider energy flow within a slice, of thickness Δz, in the xy plane. Within this slice we have P1 = P2 = P3.

P1 is the difference between the incident and reflected flux between X2 and x2 + Dr, and can be written as

 P1 = ∫x1x2〈Si ( x ) 〉Δzdx+∫x2x2+Dr〈Si ( x ) 〉Δzdx+∫x1+Drx2+Dr〈Sr ( x ) 〉Δzdx (49)

where Si(x) and Sr(x) are the intensities, represented in terms of time averaged Poynting vectors along y, of the incident and reflected waves respectively. Since we expect the reflected beam to have the same shape as the incident beam, we can write

 〈Si ( x ) 〉 = −〈Sr ( x+Dr ) 〉 (50)

so that the first and last terms on the right hand side of Equation 49 cancel. Thus only the second term contributes to P1, and the integral is only performed between x2 and x2 + Dr, in the central portion of the beam. Within this interval Si(x) has a constant value, which we denote as Smax. Equation 49 thus becomes

### Figure 16.

Model of lateral displacement of reflected beam in the reststrahl regions in terms of power flow along the antiferromagnet surface.

 P1 = 〈Smax〉ΔzDr (51)

Using Equation 25 and a standard application of Maxwell’s equations, Smax can be expressed in terms of the incident field Emax within the central region of the beam as

 〈Smax〉 = k0|Emax|22ωμ0 (52)

The flux P2 within the antiferromagnet is given by

 P2 = ∫0∞〈S2x ( y ) 〉Δzdy (53)

S2x(y) can be expressed in terms of the electric field Ez(y) within the antiferromagnet using Equation 27. Ez(y) itself can be related to the field Emax of the incident beam by the equation

 Ez ( y ) = tEmaxexp(ik2yy) (54)

so P2 may also be obtained as a function of Emax. On putting P1 = P2 in the above equations, and solving for Dr we get exactly the same result as obtained using the angular spectrum analysis (Equation 44). Thus simple energy conservation principles may be used to predict the Goos-Hänchen shift in this case.

We may examine explicitly the power flow behaviour in the xy plane using the type of calculation used in obtaining Figures 14 and 15. Figure 17(a) shows the overall power intensity and flux lines (Lai et al., 2000) for reflection off MnF2 at frequency A, corresponding to the upper reststrahl region. The behaviour is similar to that predicted from Figure 16 (although in the present case the shift is negative). Thus, for x > −0.01 cm, the incident intensity is greater than the reflected intensity, so that the overall power flow is to the right. At x −0.01 cm, the incident and reflected beams cancel, whereas for x < −0.01 cm, the reflected beam dominates and the overall power flow is to the left. Flux continuity is thus preserved in the manner illustrated in Figure 16.

### Figure 17.

Overall power intensity and flow for reflection off MnF2 in an external magnetic field of 0.1 T in the upper reststrahl region (frequency A). After Lima et al. (2009).

Figure 17(b) shows power flow when damping is included. In this case the flux continuity argument no longer applies (hence we have not attempted to show continuous flux lines), but the overall behaviour is not significantly changed.

It is interesting to note that the power flow explanation of the normal incidence Goos Hänchen shift gives a very clear example of the breaking of time reversal symmetry by a static magnetic field. In the example shown in Figure 16, there is power flow along the positive x direction within the antiferromagnet, and the reflected beam is thus displaced in this direction. If, however, time reversal were applied to the reflected beam, so that it became a new incident beam, without reversing the direction of B0, the flow of energy within the antiferromagnet would not retrace its path along the negative x direction, but would once again flow along positive x, and there would be a further displacement of the reflected beam in this direction.

#### 4.3.2. Bulk regions

The above power conservation principles used in analysing the normal incidence Goos- Hänchen shift in the reststrahl regions are based on the principle that all the incident energy will be reflected back from the antiferromagnet surface. This is clearly not the case in the bulk regions. Furthermore, in the absence of damping, the power flow resulting from normally incident plane waves is transmitted normal to the surface, as seen in Figure 9(a). The reststrahl region analysis requires power flow parallel to the interface, so one might expect that there would be no lateral shift of the reflected beam in the bulk regions. The angular spectrum analysis, however, predicts that such a shift should occur in both the reststrahl and the bulk regions, so it is important to understand, in terms of energy flow, how such a shift is possible in the bulk regions.

In order to analyse the power flow, we once again take the example of reflection of a Gaussian beam and examine the power flow profile in the xy plane. The resulting overall power intensity and flux at the two bulk frequencies B and C are shown in Figure 18. Before analysing the lateral shift in any detail, it is worth noting that the profile of the power intensity within the antiferromagnet is different from that of the electric field amplitudes

### Figure 18.

Overall power intensity and flow for reflection off MnF2 in an external magnetic field of 0.1 T in the bulk regions. Note that the arrows in the bottom half of (c) are simply intended to show the direction of power flow and are not a continuation of the flux lines in the upper part of the figure. After Lima et al. (2009).

shown in Figures 14 and 15. This is because the power intensity depends on the magnetic field component of the electromagnetic field as well as its electric field component. The magnetic field profile is in fact very different from that of the electric field (Lima et al., 2009), leading to a different power intensity profile.

At frequency B, ignoring damping (Figure 18(a)), overall power flow is to the right both in the vacuum region y < 0 and within the antiferromagnet y > 0, and remains perpendicular to the surface. There is flux continuity across the interface, as expected, and the overall profile is slightly displaced in the positive x direction with respect to the incident beam, which is centred at x = 0. We can understand this if we recall that the reflected beam is slightly displaced in the negative x direction and is less intense than the incident beam (see Figure 13(b)). The overall flux for y < 0, within the vacuum region, is the incident beam minus the reflected beam, so it is shifted in the positive x direction with respect to the pure incident beam. Flux continuity requires that this displacement is transferred to the transmitted beam, i.e. in the bulk regions the transmitted beam is displaced along x in the direction opposite to the displacement of the reflected beam. The behaviour in the presence of damping (Figure 18(b)) is not appreciably different.

At frequency C in the absence of damping (Figure 18(c)) the situation is similar to that at frequency B except that there is now a resultant power to the left for x < −0.08 cm−1. This is related to the fact that the reflected beam is more intense than the incident beam at these values of x (see Figure 13(c)). This also means that there must be power flow to the left within the antiferromagnet. This presents an apparent problem since there is no obvious source for the associated energy. Nevertheless, detailed calculations of power flow several centimeters into the antiferromagnet show that a portion of the incident energy does indeed return to the left in this region (Lima et al., 2009). This is easier to see when damping is present (Figure 18(d)), in which case some of the incident energy returns to the left without penetrating a long distance into the antiferromagnet. In fact, in this case, the overall power flow is somewhat similar to that observed in the reststrahl regions. This is another illustration of the concept that, in the presence of damping, the bulk and reststrahl regions can be thought of as merging into one another.

In the above examples, frequency B (Figures 18(a) and 18(b)) can be regarded as representative of most of the bulk frequencies. Frequency C (Figures 18(c) and 18(d)), in contrast, displays rather peculiar behaviour particular to frequencies between the reststrahl region and the reflection minimum (see Figure 12).

## 5. Conclusions and future prospects

In this chapter we have examined various nonreciprocal effects associated with reflection of terahertz radiation off antiferromagnets. Of these effects, only nonreciprocity in the reflectivity has, to our knowledge, been investigated experimentally at the time of writing (Remer et al., 1986, Brown et al., 1994).

A simple, if slightly indirect, way of observing the nonreciprocal reflected phase has been suggested (Dumelow & Camley, 1996, Dumelow et al., 1998), and uses the configuration shown in Figure 19. Here, a dielectric layer is deposited onto the surface of an antiferromagnet and the overall reflectivity off the overall structure measured. In this setup, there is reciprocal reflection from the vacuum/dielectric interface, but the phase of the radiation reflected from the dielectric/antiferromagnet interface is nonreciprocal. Interference between these partial waves is thus nonreciprocal, leading to a nonreciprocal

#### Figure 19.

Use of a dielectric layer for investigating nonreciprocal phase on reflection off an antiferromagnet.

#### Figure 20.

Calculated oblique incidence reflectivity spectrum off a Si/MnF2 structure of the type shown in Figure 19 in an external field of B0 = +0.1 T (solid curves) and B0 = −0.1 T (dashed curves).

overall reflectivity which depends on the dielectric layer thickness, as shown in Figure 20. This is true even when the reflectivity off the pure antiferromagnet is close to reciprocal, as is the case for MnF2.

The discussion of nonreciprocity in the power flow is concerned with power flow behaviour within the interior of an antiferromagnet. Obviously it is not straightforward to measure this experimentally. It appears more reasonable to investigate the effect of this nonreciprocal power flow on the radiation interacting with a finite sized sample of a given shape. The analysis presented in this chapter does not extend to this type of system, since the antiferromagnet is considered to be infinite along x. However, other techniques such as the finite difference time domain (FDTD) method should help clarify the expected behaviour.

The lateral shift predicted in the case of reflection of a finite beam off an antiferromagnet should in principle be measurable given a suitable coherent source such as a far infrared laser (Rosenbluh et al., 1976), backward wave oscillator (Dobroiu et al., 2004), or YIG oscillator with frequency multiplied output (Kurtz et al., 2005). In order to observe the normal incidence shift, a beam splitting arrangement appears necessary. It is also, however, important to consider the effect at oblique incidence, both theoretically and experimentally. In this case the effect should be observable directly without the use of a beamsplitter.

In this chapter we have only discussed phenomena in the Voigt configuration with the external field aligned along the anisotropy axis, deliberately avoiding the more complicated configurations in which the external field makes an angle with the anisotropy field (Almeida & Mills, 1988), or in which these axes are not perpendicular to the plane of incidence. However, we should point out that theoretical works on the reflected amplitude and phase do exist for more complex geometries (Stamps et al., 1991, Dumelow et al., 1998), and, in the case of reflectivity, there is some experimental work (Abraha et al., 1994, Brown et al., 1995).

Finally, we stress that, although we have concentrated on reflection off antiferromagnets in this chapter, the basic priciples involved stem from the form of the permeability tensor given in Equation 8. However, there are other types of material, such as ferromagnets or ferrimagnets, that also have a gyromagnetic permeability of this form. We therefore expect similar phenomena for these materials, although some of the symmetry arguments have to be looked at in a slightly different way since, in general, such materials have their own internal macroscopic magnetic field. One can also have a dielectric tensor of this form, such as that associated with magnetoplasma excitations. In this case, p-polarisation radiation should give results similar to those presented here for s-polarisation reflection off antiferromagnets (Remer et al., 1984).

## Acknowledgements

This work was partially financed by the Brazilian Research Agency CNPq (projects Universal 482238/2007-0, CT-ENERG 554889/2006-4, and CNPq-Rede NanoBioestruturas 555183/2005-0).

## References

1 - K. Abraha, D. E. Brown, T. Dumelow, T. J. Parker, D. R. Tilley, 1994 Oblique incidence far-infrared reflectivity study of the uniaxial antiferromagnet fef2, Phys. Rev.B 50 10 6808 6816 .
2 - K. Abraha, D. R. Tilley, 1996 Theory of far infrared properties of magnetic surfaces, films and superlattices, Surf. Sci. Rep. 24(5-6): 129-222.
3 - N. S. Almeida, D. L. Mills, 1988 Dynamical response of antiferromagnets in an oblique magnetic field: Application to surface magnons, Phys. Rev. B 37 7 3400 3408 .
4 - K. Artmann, 1948 Berechnung der seitenversetzung des totalrelektierten strahles, Ann. Physik 437 87 102 .
5 - M. . Born, E. Wolf, 1980 Principles of Optics, sixth edn, Pergamon Press, Oxford.
6 - D. E. Brown, T. Dumelow, T. J. Parker, K. Abraha, D. R. Tilley, 1994 Nonreciprocal reflection by magnons in FeF2: a high resolution study, Phys. Rev. B 49 17 12266 12269 .
7 - D. E. Brown, T. Dumelow, T. J. Parker, K. Abraha, D. R. Tilley, 1995 Far infrared reflectivity off FeF2, J. Magn. Magn. Mater. 140-144: 181-182.
8 - D. E. Brown, T. Dumelow, T. J. Parker, M. R. F. Jensen, D. R. Tilley, 1998 A high resolution fourier transform spectrometer for far infrared magneto-optic spectroscopy of magnetic materials, Infrared Physics and Technology 40 219 230 .
9 - R. E. Camley, 1987 Nonreciprocal surface modes, Surf. Sci. Rep. 7(3-4): 103-188.
10 - A. Dobroiu, M. Yamashita, Y. N. Ohshima, Y. Morita, C. Otani, K. Kawase, 2004 Terahertz imaging system based on a backward-wave oscillator, Appl. Opt. 43 30 5637 5646 .
11 - H. Dötsch, N. Bahlmann, O. Zhuromskyy, M. Hammer, L. Wilkens, R. Gerhardt, P. Hertel, A. F. Popkov, 2005 Applications of magneto-optical waveguides in integrated optics: review, J. Opt. Soc. Am. B 22 1 240 253 .
12 - T. Dumelow, R. E. Camley, 1996 Nonreciprocal reflection of infrared radiation from structures with antiferromagnets and dielectrics, Phys. Rev. B 54 17 12232 12237 .
13 - T. Dumelow, R. E. Camley, K. Abraha, D. R. Tilley, 1998 Nonreciprocal phase behavior in reflection of electromagnetic waves from magnetic materials, Phys. Rev. B 58 2 897 908 .
14 - T. Dumelow, M. C. Oliveros, 1997 Continuum model of confined magnon polaritons in superlattices of antiferromagnets, Phys. Rev. B 55 2 994 1005 .
15 - F. Goos, H. Hänchen, 1947 Ein neuer und fundamentaler versuch zur totalreflexion, Ann. Physik 436 6 333 346 .
16 - B. R. Horowitz, T. Tamir, 1971 Lateral displacement of a light beamat a dielectric interface, J. Opt. Soc. Am. 61 5 586 594 .
17 - D. Kurtz, T. Crowe, J. Hesier, D. Porterfield, V. Inc, V. Charlottesville, 2005 Frequency domain terahertz spectroscopy, Infrared and Millimeter Waves and 13th International Conference on Terahertz Electronics, 2005., 1 76 77 .
18 - H. M. Lai, C. W. Kwok, Y. W. Loo, B. Y. Xu, 2000 Energy-flux pattern an the Goos-Hänchen effect, Phys. Rev. E 62 5 7330 7339 .
19 - L. D. Landau, E. M. Lifshitz, 1984 Electrodynamics of Continuous Media, second edn, Pergamon Press, Oxford.
20 - P. T. Leung, Z. W. Chen, H. Chiang, P. , 2007 Large negative Goos-Hänchen shift at metal surfaces, Opt. Commun. 276 2 206 208 .
21 - F. Lima, T. Dumelow, E. L. Albuquerque, Costa. J. A. P. da, 2009 Power flow associated with the Goos-Hänchen shift of a normally incident electromagnetic beam reflected off an antiferromagnet, Phys. Rev. B 79(15): 155124.
22 - F. Lima, T. Dumelow, Costa. J. A. P. da, E. L. Albuquerque, 2008 Lateral shift on normal incidence reflection off an antiferromagnet, Europhys. Lett. 83(1): 17003.
23 - H. K. V. Lotsch, 1970 Beam displacement at total reflection: the Goos-Hänchen effect, Optik 32: 116,189,299,553.
24 - M. Mc Guirk, C. K. Carniglia, 1977 An angular spectrum approach to the Goos-Hänchen shift, J. Opt. Soc. Am. 67 1 103 107 .
25 - D. L. Mills, E. Burstein, 1974 Polaritons: the electromagnetic modes of media, Rep. Prog. Phys. 37 7 817 926 .
26 - R. Potton, 2004 Reciprocity in optics, Rep. Progr. Phys. 67 5 717 753 .
27 - L. Remer, B. Lüthi, H. Sauer, R. Geick, R. E. Camley, 1986 Nonreciprocal optical reflection of the uniaxial antiferromagnet MnF2, Phys. Rev. Lett. 56 25 2752 2754 .
28 - L. Remer, E. Mohler, W. Grill, B. Lüthi, 1984 Nonreciprocity in the optical reflection of magnetoplasmas, Phys. Rev. B 30 6 3277 3282 .
29 - R. H. Renard, 1964 Total reflection: a new evaluation of the Goos-Hänchen shift, J. Opt. Soc. Am. 54 10 1190 1197 .
30 - M. Rosenbluh, R. J. Temkin, K. J. Button, 1976 Submillimeter laser wavelength tables, Appl. Opt. 15 11 2635 2644 .
31 - R. Q. . Scott, D. L. Mills, 1977 Propagation of surface magnetostatic waves on ferromagnetic crystal structures, Phys. Rev. B 15 7 3545 3557 .
32 - R. L. Stamps, B. L. Johnson, R. E. Camley, 1991 Nonreciprocal reflection from semi-infinite antiferromagnets, Phys. Rev. B 43 4 3626 3636 .
33 - W. J. Wild, C. L. Giles, 1982 Goos-Hänchen shifts from absorbing media, Phys. Rev. A 25 4 2099 2101 .