10 1063@1 4871098
10 1063@1 4871098
Variable angle magnetometry for exchange-coupled multilayers with in-plane and perpendicular anisotropy
J. Appl. Phys. 110, 123905 (2011); 10.1063/1.3665191
Strong in-plane anisotropy of magneto-optical Kerr effect in corrugated cobalt films deposited on highly ordered
two-dimensional colloidal crystals
Appl. Phys. Lett. 98, 031903 (2011); 10.1063/1.3544582
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
REVIEW OF SCIENTIFIC INSTRUMENTS 85, 053904 (2014)
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
053904-2 Jiménez et al. Rev. Sci. Instrum. 85, 053904 (2014)
2 mm down to 30 μm can be used. The sample is placed on the For the sake of simplicity we use the frame of reference
eucentric goniometer head to ensure a fixed plane of reflection shown in Fig. 1. The x and y directions are lying within the
upon sample rotation. In magneto-optical measurements this film plane, the surface normal is the z direction, and the in-
is important to be able to compare the values of the magneti- cident angle θ i is contained in the xz reflection plane. All
zation components measured at different angles and between calculations presented here are performed using the Jones
different samples. The whole head can be rotated by a step- formalism.28 It is important to note that rss is independent of
ping motor in steps of 0.9◦ . The magnetic field is provided by magnetization. The components rsp and rps are proportional
a home-made air-gap electromagnet. We have developed and while rpp contains an additional constant term. Upon reflec-
calibrated several electromagnets with different gaps and/or tion on the surface of the magnetic sample the electric field
coil turns, depending on the frequency and maximum field vector E of the incident light is reflected and results in E as
required. Magnetic fields up to ±100 mT driven at frequen-
cies up to 1 kHz can be obtained. This enables performing E = rE.
(2)
dynamic studies over 10 orders of magnitudes of applied field Naturally, the final measurement procedure consists of mea-
sweep rates, from 10−4 mT/s to 10+6 mT/s. In the present suring the intensity, i.e., I = E · (E )∗ where the asterisk de-
study, we have used a 10 mm air-gap electromagnet mounted notes the complex conjugate. It is easily seen that the intensity
on a rotatory stage, allowing a controlled way to change from of the reflected s -polarized light due to incoming p-polarized
longitudinal to a transversal MOKE geometry. Polar geome- light (or vise versa) results in a signal that is quadratic in m.
try can also be set. Consequently, such a measurement does not give proper in-
The control unit consists of a computer, an arbitrary func- formation on the magnetic state of the sample. However, the
tion generator driving a bipolar current source, and a four measurement of the reflected p -polarization of incoming p-
channels digital oscilloscope. Via a home-made measurement polarized light gives
software, the computer controls the stepper motor of the go-
niometer as well as the applied magnetic field ramp. The latter Ip = Ep (Ep )∗ = cc∗ + 2Re[cd ∗ ]my + dd ∗ m2y
is driven by the bipolar current source programmed with the
arbitrary function generator. Additionally, the software reads ≈ cc∗ + 2Re[cd ∗ ]my + O m2y , (3)
out the hysteresis data, i.e., the signals from the digital oscillo-
scope triggered by the trigger output of the arbitrary function assuming the typical case where Re[cd ∗ ] dd ∗ . Hence, the
generator. The signals are constant term in rpp allows for easy measurement of my in a
transverse Kerr setup.
1. The applied voltage to the coil: applied field ramp. The presented setup is based on the idea of introducing
2. I : the difference of intensities on the diodes. an additional constant term in the Es component, therefore
3. IDC : the DC component of the sum of intensities. generating a linear term in the corresponding intensity. This
4. IAC : the AC component of the added intensities. is achieved by introducing a λ/2-retarder L2 , rotated by π /8
The software allows adjustment of the parameters to perform off the optical axes, and by reading the two orthogonal com-
specific measurement procedures in order to carry out vecto- ponents of the reflected light. In general, the reflected signal
rial resolved magnetization reversal studies in the whole an- has the form
gular range. These may include major and minor loops, and/or Es Es
first order magnetization curves as well as quasi-static and dy- = L2 r
Ep Ep
namical studies.
i 1 1 rss rsp Es
= √ . (4)
III. THEORETICAL DETAILS OF THE SETUP 2 1 −1 rps rpp Ep
All Kerr measurements are based on the fact that the mag- For simplification let us assume only in-plane magnetization,
netization alters the dielectric tensor such as making an other- i.e., mz = 0, and incoming pure p-polarized light. The electric
wise isotropic material optically anisotropic. Proper consid- field vector after passing the retarder then has the form
erations of the boundary conditions of Maxwell’s equations
enable the calculation of the corresponding reflection matrix Es b1 mx + c + dmy
= Ep , (5)
of a magnetic material.27 Typically the symmetry breaking Ep b1 mx − c − dmy
due to the magnetization results only in small modifications
of the reflection matrix such that the optical anisotropy can and the intensities read
be handled as a first order perturbation. The components of I0 2
Is = |c| + |b1 |2 m2x + |d|2 m2y + 2Re[b1 d ∗ ]mx my
the reflection matrix for a thin magnetic film are given in the 2
Appendix. It is found that in the first order the matrix is of the
form + 2Re[b1 c∗ ]mx + 2Re[dc∗ ]my ,
(6)
I0 2
rss rsp a b1 mx + b2 mz Ip = |c| + |b1 |2 m2x + |d|2 m2y − 2Re[b1 d ∗ ]mx my
r= ≈ . 2
rps rpp −b1 mx + b2 mz c + dmy
(1) − 2Re[b1 c∗ ]mx + 2Re[dc∗ ]my ,
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
053904-4 Jiménez et al. Rev. Sci. Instrum. 85, 053904 (2014)
where I0 is the intensity of the incident beam. Adding and whereas for transversal MOKE geometry, i.e., By ,
subtracting the two intensities gives
IAC
∗ ∗ M = my ∝ , (14)
I = Is − Ip = 2I0 {Re[b1 c ]mx + Re[b1 d ]mx my } IDC
≈ 2I0 Re[b1 c∗ ]mx + O(mx my ), I
p
M⊥ = mx ∝ 2θK ≈ . (15)
IDC
I = Is + Ip = I0 |c|2 + |b1 |2 m2x + |d|2 m2y
Notice that each in-plane magnetization component can be
(7)
+ 2Re[b1 d ∗ ]mx my + 2Re[dc∗ ]my derived independently from Kerr rotations or reflectivity
changes, depending on the MOKE geometry. This will be
≈ I0 {|c|2 + 2Re[dc∗ ]my } used to determine experimentally the scale factor between
both Kerr effects to obtain quantitative information of the in-
+ O m2x , m2y , mx my . plane resolved hysteresis loops. This is further discussed in
Sec IV.
As a result the difference signal is proportional to mx , whereas
the sum is linear in my , but has an additional DC-component,
i.e., IDC = |c|2 = |rpp | . Note that due to this DC component
DC 2 A. Errors due to misaligned optical components
one finally gets
In the previous paragraph the ideal case has been dis-
DC ∗ rsp DC DC ∗
I 2Re[rsp rpp ] 2Re r DC rpp rpp cussed. In experiment, however, one must assume small an-
= =
pp
gular errors in the optical components, i.e., polarizer, λ/2-
IDC r DC 2 r DC 2 retarder, and analyzer. Let us assume that these components
pp pp
have angular errors of α 1 , α 2 , and α 3 . In the following errors
rsp p of the order α i α j (i, j = 1, 2, 3) as well as α i mξ mζ (ξ , ζ = x,
= 2Re ≈ 2θK ∝ mx , (8)
DC
rpp y, z) are neglected. In this approximation the final intensities
including errors have the form
AC
IAC rpp I ∝ I(1) + 2Re[b1 d ∗ ]mx my + 2Re[b2 d ∗ ]mz my
DC
= 2Re DC ∝ my , (9)
I rpp
+ 2α1 (Re[ac∗ ] + Re[ad ∗ ]my )
which are the Kerr rotation and reflectivity change, which are
proportional to the two in-plane components of the magneti- + 2(2α2 − α3 )I(1) , (16)
zation mx and my , respectively.
Similarly for incoming s-polarized light one obtains for and
the difference and the sum terms I ∝ I(1) + |d|2 m2x + |b1 |2 m2y + |b2 |2 m2z
I rps + 2Re[b1 b2∗ ]mx mz
= 2Re DC ≈ 2θKs ∝ mx , (10)
IDC rpp
+ 2α1 (Re[a ∗ b1 − b1 c∗ ]mx + Re[a ∗ b2 + b2 c∗ ]mz )
I I DC
= = Rss , (11)
I0 I0 − α3 I(1) , (17)
the Kerr angle and a constant term, respectively. In this case, where I(1) and I(1)are the error free first order approxima-
the sum, I , does not contain useful information other than tions of the signals. Furthermore, the summation and subtrac-
Rss = |rss |2 , as I ∝ Rss + O(m2y ). Consequently, there is no tion, resulting in I and I , are assumed to be without error.
IAC in the first order. Note that α 2 does not effect I as its effect cancels out in the
Therefore, the combination of p-polarized light and the sum, but mixes I into I . Moreover, α 3 intermixes I and I ;
detection of the two orthogonal components of the reflected however, due to the adding and subtracting of the signals, I
light at the same time allows the determination of both in- is doubly affected. Furthermore, one must keep in mind that
plane magnetization components simultaneously. Notice that the measured intensities, at both diodes, are dominated by the
this is accomplished independently of the MOKE geometry. term Ip, s ≈ I0 /2|c|2 . This term defines the photon noise. As-
In this sense, we define the in-plane magnetization compo- suming a Poisson statistic the resulting photon noise in each
nents parallel M and perpendicular M⊥ with respect to the ex- channel is proportional to |c|. Hence, the linear approxima-
ternal magnetic field direction. In correspondence with Fig. 1, tion becomes better with increasing c and the relative photon
for longitudinal MOKE geometry, i.e., Bx , the in-plane mag- noise decreases with |c|−1 .
netization components can be derived from To minimize the errors of α 1 and α 3 one first omits the
I retarder plate and introduces a non-magnetic reflecting ma-
p
M = mx ∝ 2θK ≈ , (12) terial, i.e., without off-diagonal elements. Assuming that one
IDC
wishes incoming s-polarized light, the according intensity is
IAC at the lowest order approximation of the form
M⊥ = m y ∝ , (13)
IDC Ip ∝ (rpp α1 − rss α3 )2 . (18)
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
053904-5 Jiménez et al. Rev. Sci. Instrum. 85, 053904 (2014)
This equation, however, has infinite solutions. Therefore, the cleaned by repeated cycles of sputtering with 0.5 kV Ar+
optimizing process consists in three steps. First one of the po- ions and annealing up to 800 K. Fe deposition was performed
larizers is repeatedly rotated by 180◦ , such that the beam is at 400 K with a home-made electron-beam evaporator at a
entering from the opposite site. This operation transforms the rate of about 0.4 nm/min. After deposition, the samples were
according angle into its negative and the intensity will remain annealed to about 750 K, in order to reduce possible de-
unchanged only if the angle is 0. In the second step the in- fects formed during the epitaxial growth, and characterized by
tensity in the p-channel must be minimized for the second low energy electron diffraction (LEED). The Fe film grows
polarizer. As the minimum is quadratic, a precise adjustment according to the well-known Fe(001)[110] MgO(001)[100]
can be difficult. This can be avoided by finding two opposite epitaxial relation, as indicated by diffraction LEED patterns
angles with identical intensity, such that the minimum will be (not shown). The film was then capped at room tempera-
at the average of these two angles. This procedure is applied ture with a 3 nm thick Cu film, to prevent oxidation. The
when incoming p-polarized light is required. Finally the re- four-fold crystal symmetry of the Fe(100) film promotes
tarder is reintroduced. The retarder position is then optimized a pure biaxial magnetic anisotropy in the film, as shown
by making the intensities in the Is and Ip channel identical, below.
thereby minimizing I .
IV. EXPERIMENT
B. The analysis procedure
The capabilities of our v-MOKE setup are shown by pre-
senting detailed angular measurements of a 20 nm thick Fe A typical experimental data set for one angular condition
film epitaxially grown on a MgO(100) single crystal sub- is shown in Fig. 3. In this case, longitudinal geometry, incom-
strate in ultra-high-vacuum conditions. The substrate was first ing p-polarized light, and α H = 18◦ were used. α H = 0◦ is
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
053904-6 Jiménez et al. Rev. Sci. Instrum. 85, 053904 (2014)
FIG. 4. M (H) curves acquired at the easy axis in longitudinal (left y-axis)
and transversal (right y-axis) MOKE geometries, as depicted schematically
on top. Notice that the hysteresis curves are obtained from different Kerr
effects, i.e., polarization rotation and reflectivity change, respectively. The
scale factor is 0.4.
FIG. 6. In-plane resolved magnetization hysteresis curves of a 20 nm Fe film grown on MgO(001) acquired at selected angles. Top: Standard-plot M–H
representation. M (H, α H ) and M⊥ (H, α H ) loops are represented by circles and squares, respectively. The descending (forward) and ascending (backward) field
branches are depicted with filled and empty symbols, respectively. The remanences MR, and MR,⊥ and the critical fields HC,1 and HC,2 are indicated. Bottom:
Corresponding polar-plot M⊥ − M representation. Filled and open circle symbols have been used for the forward and the backward branch, respectively. The
circle of radius unity depicted with a solid line indicates the pure rotation of the magnetization vector. The domain wall angles 1 and 2 are indicated. The
letters in the curves denote the sequence of the reversal process as discussed in the text. Solid arrow ended and dashed straight lines indicate the positions of the
e.a. and h.a. directions, respectively.
180◦ or 90◦ away from B, depending on α H . The crossover B. Remanence and magnetic symmetry
between both behaviors is found at αC . For the first case, e.g.,
for α H = 10◦ , after the first irreversible transition the mag- The simplest information that can be extracted from the
netization vector rotates along the circle until the maximum angular evolution is the anisotropy determined by the mag-
field is reached (D). In turn, e.g., for α H = 33◦ , a second irre- netic symmetry of the system. This can be done easily by
versible transition takes place (D), with the return point near plotting the normalized remanence values of the two in-
to the closest e.a. direction (E), which is ca. orthogonal to D. plane magnetization components as a function of α H (see
Thereafter, rotation of magnetization continues until the max- Fig. 7(a)). Both magnetization components display a pro-
imum field is reached (F). Note that the magnetization vector nounced oscillation with periodicity of 90◦ . The parallel com-
is far to be saturated along the field direction even for the ponent follows a |cos 2α H | law dependence, the perpendicu-
largest field used, except for α H = 0◦ . lar component changes the sign when a characteristic e.a. or
Therefore, in a pure fold-fold magnetic symmetry sys- h.a. direction is crossed, and both components are comple-
tem, the reversal proceeds by reversible magnetization rota- mentary, i.e., MS2 = (MR, 2
+ MR,⊥2
). The polar-plots of MR,
tion processes, in order to get the closest e.a. direction, with (Fig. 7(b)) and M⊥,R (Fig. 7(c)) display characteristic “four-
one (two) irreversible transition, related to nucleation of an- leaves clover” and “four-wings windmill” shapes, respec-
tiparallel (orthogonal) magnetic domains and further propaga- tively. This originates from the four-fold magnetic symmetry,
tion of 180◦ (90◦ ) domain walls, when the applied field direc- due to the cubic crystal symmetry of the Fe(100) film. All
tion is close to one of the two easy- (hard-) axis magnetization these features confirm a pure biaxial magnetic anisotropy be-
directions. In addition, the magnetic domains are not oriented havior of the film, where the anisotropy axes are aligned paral-
necessarily parallel to the field direction but rather to the e.a. lel to the in-plane [001] crystallographic directions of the film.
directions. It is worth mentioning that, in the present case for a pure
For a more quantitative analysis, magnetic parameters, biaxial anisotropy system, the information on the magnetic
such as remanence magnetization, coercivity, and domain symmetry does not require the measurement of MR,⊥ (H ).
wall angle, can be readily obtained as a function of α H from However in a magnetic system with competing anisotropies,
the hysteresis loops, as indicated in Fig. 6. the second magnetization component allows for a more
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17
053904-9 Jiménez et al. Rev. Sci. Instrum. 85, 053904 (2014)
ACKNOWLEDGMENTS
We would like to thank the research support services at
UAM, SEGAINVEX, for the development of the detectors
and power supplies, and to gratefully acknowledge the tech-
nical crew of the UAM Condensed Matter Department. P.P.
acknowledges financial support from MINECO through Con-
tract No. JCI-2011-09602. This work has been supported by
FIG. 8. Angular dependence of the critical fields of the four-fold crystal sym-
the MEC through Project No. CSD2007-00010, by MINECO
metry Fe(100) film. HC,1 (solid symbols) and HC,2 (open symbols) are the through Project Nos. FIS2010-18847, MAT2011-25598, and
critical field values of the first and second irreversible transition, respectively. MAT2012-39308, and by the Comunidad de Madrid through
For accuracy reasons, the values of HC,1 are determined from the parallel Project No. S2009/MAT-1726. N.M. acknowledges funding
component, while the values HC,2 are extracted from the perpendicular com-
ponent (see text). Notice that near the easy axes the results coincide in a single
from the French National Research Agency (ANR) under the
value, i.e., single transition by a domain wall of approximately 180◦ , while Project NPOEM No. ANR- 10-NANO-0016.
two transitions are observed around the hard axes (emphasized by the light
gray background).
APPENDIX: REFLECTION COEFFICIENTS
To justify the form of the reflection matrix given in
where the reversal takes place via two consecutive irreversible Eq. (1) we calculated the matrix according to the boundary
transitions. matrix method of Qiu and Bader.4 Within this framework one
Finally, it is worth highlighting the strong influence of multiplies the medium boundary (Ai ) and medium propaga-
the magnetic symmetry on the magnetic properties of the tion matrices (Dj ) to obtain the matrix M defined by
Fe(100) film, and how they are strongly related to each other.
This can be generalized for other systems with more com- G H
M = A−1 i
N −1
m=1 Am Dm Am Af ≡ , (A1)
plex magnetic symmetries.6, 17–22, 25, 26 In this case, the biaxial I J
magnetic anisotropy is promoted by the four-fold crystal sym-
metry of the Fe film and determines hysteretic phenomena, where the Fresnel reflection and transmission coefficients are
remanences, critical fields, domain wall angles, and magneti- defined as
zation reversal processes. tss tsp rss rsp
−1 −1
G = ; IG = . (A2)
tps tpp rps rpp
where ni , nf , and nm indices of refraction of the initial, the fi- 16 J. Camarero, Y. Pennec, J. Vogel, M. Bonfim, S. Pizzini, F. Ernult, F.
nal, and the magnetic medium. The incident angle is θ i , while Fettar, F. Garcia, F. Lanc¸on, L. Billard, B. Dieny, A. Tagliaferri, and N.
B. Brookes, Phys. Rev. Lett. 91, 027201 (2003).
θ m and θ f are the (complex) angles of propagation in the mag- 17 T. Bobek, N. Mikuszeit, J. Camarero, S. Kyrsta, L. Yang, M. A. Niño,
netic layer and the substrate, respectively. The magneto-optic C. Hofer, L. Gridneva, D. Arvanitis, R. Miranda, J. J. de Miguel, C. Te-
effect also depends on the material constant Q—the Voigt ichert, and H. Kurz, Adv. Mater. 19, 4375 (2007).
18 P. Prieto, K. R. Pirota, J. M. Sanz, E. Jiménez, F. Maccherozzi, and F. G.
constant—that accounts for the off diagonal elements in the
Panaccione, Appl. Phys. Lett. 90, 032505 (2007).
dielectric tensor. 19 D. Ecija, E. Jiménez, N. Mikuszeit, N. Sacristán, J. Camarero, J. M.
Sci. Instrum. 71, 1243 (2000). P. Perna, J. M. García-Martín, B. Dieny, and R. Miranda, Appl. Phys. Lett.
5 A. Hubert and R. Schäfer, Magnetic Domains (Springer-Verlag, Berlin). 95, 122508 (2009).
6 J. Camarero, J. Sort, A. Hoffmann, J. M. García-Martín, B. Dieny, 22 P. Perna, C. Rodrigo, E. Jiménez, F. J. Teran, N. Mikuszeit, L. Méchin, J.
R. Miranda, and J. Nogués, Phys. Rev. Lett. 95, 057204 (2005). Camarero, and R. Miranda, J. Appl. Phys. 110, 013919 (2011).
7 J. M. Florczak and E. Dan Dahlberg, J. Appl. Phys. 67, 7520 (1990). 23 J. M. Teixeira, R. Lusche, J. Ventura, R. Fermento, F. Carpinteiro, J. P.
8 C. Daboo, J. A. C. Bland, R. J. Hicken, A. J. R. Ives, M. J. Baird, and M. Araujo, J. B. Sousa, S. Cardoso, and P. P. Freitas, Rev. Sci. Instrum. 82,
J. Walker, Phys. Rev. B 47, 11852 (1993). 043902 (2011).
9 R. P. Cowburn, S. J. Gray, J. Ferr, J. A. C. Bland, and J. Miltat, J. Appl. 24 P. Perna, C. Rodrigo, E. Jiménez, N. Mikuszeit, F. J. Teran, L. Méchin,
Phys. 78, 7210 (1995). J. Camarero, and R. Miranda, J. Appl. Phys. 109, 07B107 (2011).
10 R. M. Osgood, III, B. M. Clemens, and R. L. White, Phys. Rev. B 55, 8990 25 E. Jiménez, J. Camarero, P. Perna, N. Mikuszeit, F. J. Teràn, J. Sort, J.
and J. Wollschläger, J. Phys. D 44, 265003 (2011). (North-Holland Publishing Company, Amsterdam, 1977).
This article is copyrighted as indicated in the article. Reuse of AIP content is subject to the terms at: http://scitationnew.aip.org/termsconditions. Downloaded to IP:
130.63.180.147 On: Mon, 11 Aug 2014 09:08:17