Study of magnetization reversal and anisotropy of single crystalline ultrathin Fe/MgO (001) film by magneto-optic Kerr effect
Zhang Miao-Ling1, Ye Jun1, Liu Rui1, Mi Shu1, Xie Yong1, Liu Hao-Liang2, Haesendonck Chris Van2, Chen Zi-Yu1, †,
Key Laboratory of Micro-nano Measurement-Manipulation and Physics (Ministry of Education), Beihang University, Beijing 100191, China
Solid State Physics and Magnetism Section, KU Leuven, Celestijnenlaan 200 D, BE-3001 Leuven, Belgium

 

† Corresponding author. E-mail: chenzy@buaa.edu.cn

Project supported by the National Natural Science Foundation of China (Grant Nos. 11274033, 11474015, and 61227902), the Research Fund for the Doctoral Program of Higher Education of China (Grant No. 20131102130005), and the Beijing Key Discipline Foundation of Condensed Matter Physics.

Abstract
Abstract

The magnetization reversal process of Fe/MgO (001) thin film is investigated by combining transverse and longitudinal hysteresis loops. Owing to the competition between domain wall pinning energy and weak uniaxial magnetic anisotropy, the typical magnetization reversal process of Fe ultrathin film can take place via either an “l-jump” process near the easy axis, or a “2-jump” process near the hard axis, depending on the applied field orientation. Besides, the hysteresis loop presents strong asymmetry resulting from the variation of the detected light intensity due to the quadratic magneto-optic effect. Furthermore, we modify the detectable light intensity formula and simulate the hysteresis loops of the Kerr signal. The results show that they are in good agreement with the experimental data.

1. Introduction

The magnetization reversal processes of Fe ultrathin films play an important role in magnetic-storage materials and have attracted sustained interest for both fundamental and technological reasons over the past few decades.[14] The magnetization reversal processes of Fe films occur via coherent rotation and/or domain wall displacement.[57] For Fe (001) single crystalline films, the magnetization reversal processes are primarily dominated by domain wall displacement and can be treated as coherent rotations only at a high external field.[810]

Owing to the high sensitivity down to sub-nm thickness, experimental simplicity and local probing nature, the magneto-optic Kerr effect (MOKE) has been widely employed to investigate the magnetization reversal process.[11] The MOKE exploits the alteration of the polarization state of light reflected from a magnetic surface, which is exceedingly applicable to the thin films.[1214] In the measurements of in-plane magnetization, the Kerr signal depends on the first-order terms of magnetization (ML and MT) as well as the second-order terms of the magnetization (ML MT, , and ).[1518] ML and MT are the longitudinal and transverse magnetization components, respectively. In various conditions, a quadratic magneto–optic effect is unavoidable.[15,19]

A coherent rotation model and a domain wall displacement model have been proposed to account for the magnetization reversal process.[20] Based on the coherent rotation model, the magnetization will change its orientation to minimize the energy of the system to reach an equilibrium state. While in a coercive field, the magnetization will execute reorientation by nucleation and propagation of new domains to lower the total system energy as indicated by the domain wall displacement model. In the coherent rotation stage, the orientation of the magnetization can be obtained by differentiating the equation of system energy, which is extremely dependent on the magnetic anisotropy constant and the applied field. Yan et al.[21] have fitted the hysteresis loops including first-order and second-order magneto–optic terms, which shows that the fitting results are in excellent agreement with the experimental results. However, a systematical and accurate method to determine anisotropy constants is still lacking. The rotational magneto–optic Kerr effect (ROTMOKE) is a powerful method to determine the magnetic anisotropy in magnetic thin films by realizing the coherent rotation of magnetization.[2224]

In the present paper, the magnetization reversal process in Fe/MgO (001) film is investigated by MOKE. Firstly, the normalized angular-dependent longitudinal and transverse magnetization hysteresis loops are measured with the applied field parallel and perpendicular to the incident plane. The result shows that both coherent rotation and domain wall displacement exist during the magnetization reversal process. Besides, the hysteresis loops exhibit strong asymmetry when the field applied is adjacent to the magnetic hard axis. We further employ ROTMOKE to determine the magnetic anisotropy and simulate the hysteresis loops of the Kerr signal by solving the energy equation and modifying the expression of normalized light intensity including first-order and second-order magneto–optic terms. The obtained results validate the contribution of the quadratic magneto–optic effect to the asymmetry of the hysteresis loops.

2. Experiments

Fe/MgO (001) film was grown by molecular beam epitaxial (MBE). The base pressure in the sputtering chamber was better than 3 × 10−10 mbar (1 bar = 105 Pa). The substrate was first annealed at 700 °C for one hour and held at 150 °C during deposition. Fe film was deposited on the MgO (001) substrate at a rate of 0.1 Å/s by an electron-beam gun. The Fe beam was incident along the surface normal of the substrate. The nominal thickness of the sample was 20 nm monitored by a calibrated quartz-crystal oscillator.[25]

In order to more comprehensively investigate the magnetization behaviors in-plane, it is necessary to measure the two orthogonal magnetization components ML and MT. The ML and MT are the magnetization components parallel and perpendicular to the applied field, respectively. In this paper, the incident angle of the laser beam on the sample was about 45°. The incident light through the polarizer was S polarized light. The analyzer was set to detect the reflected P-polarized light. In the longitudinal Kerr effect configuration, the hysteresis loops were detected by manipulating the magnetic field parallel or perpendicular to the incident plane.[26] The relative position between the magnetic field and the sample was fixed, which permits a direct comparison between the longitudinal (ML) and transverse (MT) magnetization component loops.

3. Results and discussion

The normalized remanence and coercivity are derived from the hysteresis loops measured by sample rotation with longitudinal Kerr geometry. As respectively shown in Figs. 1(a) and 1(b), the angular dependences of normalized remanence Mr/Ms and coercivity HC on field orientation clearly verify the in-plane four-fold symmetry for Fe/MgO (001) thin film. However, the measured values of both normalized remanence and coercivity appear to be conspicuously different at the equivalent position from the hard axis. The origin of this abnormal phenomenon is that for the hysteresis loops measured around the hard axis their two branches corresponding to the magnetization reversal processes are asymmetric. Figure 1(c) shows the hysteresis loop measured at φH = 31°, where φH is the angle between the applied positive field and the Fe [001] axis. Since the value of the normalized remanence is larger than 1 when we use the value of the Kerr signal at the saturation field as a maximum and a minimum for normalization separately, we ascribe the asymmetry in hysteresis loops to the consequence of the quadratic magneto–optic effect on the shape of hysteresis loop,[18] which is obvious on the Kerr signal at the low field when the value of the transversal magnetization component is large and has no effect on the magnetization reversal process.[17,18] In the following, we will discuss how to clarify the consequence of the quadratic magneto–optic effect from the measured Kerr signal and show the detailed method. However, we first give the calculated hysteresis loop at φH = 31° as shown in Fig. 1(d), in which the normalized value is determined by the linear magneto–optic effect. As clearly shown, the asymmetry in the hysteresis loop, caused by the contribution of the quadratic magneto–optic effect will greatly change the value of the normalized remanence and coercivity. By using the same method, we modify the polar plot for the angular dependences of normalized remanence Mr/Ms and coercivity HC on the field orientation without the contribution of the quadratic magneto–optic effect as given in Figs. 1(e) and 1(f). The result indicates that both the normalized remanence Mr/Ms and coercivity HC decrease as the magnetic field rotates to the hard axis and exhibit minima at the hard axis. Moreover, at the equivalent position from the hard axis the values of the normalized remanence and coercivity are comparable.

Fig. 1. Measured angular dependences of normalized remanence Mr/Ms (a) and coercivity HC (b) on field orientation. The hysteresis loops with (c) and without (d) the contribution of the quadratic magneto–optic effect at φH = 31°, where φH is the angle between the applied positive field and Fe [001] axis; the modified angular dependences of normalized remanence Mr/Ms (e) and coercivity HC (f) on field orientation. The unit 1 Oe = 79.5775 A·m−1.

Figure 2 shows the longitudinal and transverse Kerr signals depending on the applied field. The hysteresis loops each include two branches, ABCDEF and FGHIJA. On account of the magnetocrystalline anisotropy of Fe/MgO (001), the magnetization reversal process takes place via coherent rotation and/or domain wall placement with the magnetic field applied in the Fe (001) plane. Although both coherent rotation and domain wall displacement can be clearly displayed from the longitudinal or the transverse hysteresis loops, integrated investigation of these two orthogonal components helps to determine the direction of the magnetization vector. As shown in Fig. 2(a), at φH = 44°, when the applied field strength decreases from a positive saturated state to zero, the magnetization changes via coherent rotation from the direction of the applied field to the nearest Fe [100] direction, which corresponds to the A1B1 period. Increasing the negative applied field strength to the first switching field (about −26 Oe, 1 Oe = 79.5775 A·m−1), 90° domains with magnetization orientation close to [0-10] are nucleated and unpinned at the edge of the sample, followed by a speedy 90° domain wall displacement.[21] corresponding to the B1C1 period. At the section C1D1, the magnetization proceeds via coherent rotation towards Fe [0-10]. When the negative magnetic field strength increases up to the second switching field (about −210 Oe), the magnetization rotates closely to [-100] by the second 90° domain wall displacement at the D1E1 section. Further increasing the negative applied field strength, the magnetization rotates towards the direction of the applied field by a coherent rotation before saturation (E1F1 stage). In the F1G1H1I1J1A1, the change of the magnetization is similar to the A1B1C1D1E1F1.[21] The magnetization proceeds with coherent rotation at the stages of F1G1, H1I1, and J1A1, and 90° domain wall displacement takes place at the stages of G1H1 and I1J1.

Fig. 2. Hysteresis loops measured from the longitudinal Kerr signal (hollow circle) and transverse Kerr signal (hollow square). The hysteresis loops are corresponding to φH = 44° (a) and 46° (b), φH being the angle between the applied positive field orientation and Fe [100] axis. The azimuth angle is illustrated in the inset.

The hysteresis loop in Fig. 2(b) shows that the magnetization reversal process of φH = 46° also occurs by coherent rotation and 90° domain wall displacement. However, the magnetization behaviors are different between at φH = 44° and at φH = 46°. As clearly shown in Figs. 2(c) and 2(d), when the field is reduced from the positive maximum value to zero, the magnetization of φH = 44° proceeds via coherent rotation towards [100] while the magnetization of φH = 46° proceeds via coherent rotation towards [0-10]. This result is compatible with the four-fold magnetocrystalline anisotropy of Fe (001).

In order to verify the above description about the magnetization reversal process and clarify the origin of the asymmetry of the hysteresis loops, a quantitative analysis of the magnetization reversal process is performed by simulating the longitudinal Kerr signal.

Generally, four-fold magnetocrystalline anisotropy is accompanied with a weak uniaxial magnetic anisotropy in the Fe (001) plane. We decompose the energy of uniaxial anisotropy to Fe [100] and Fe [110].[6] The total energy density of the system at the applied field H is given by

where Ku1 and Ku2 are the uniaxial magnetic anisotropy constants at Fe [100] and Fe [110] direction, respectively, K1 is the first-order cubic magnetocrystalline anisotropy constant, Ms is the saturation magnetization, H is the applied field strength, and φM is the angle between the direction of magnetization and Fe [100]. According to the Stoner-Wohlfarth model, the magnetization will rotate to a certain direction so that the total energy density of the system reaches a minimum value. Then equation (2) can be obtained by differentiating Eq. (1):

where l(φM) is the normalized torque. As indicated by Eq. (2), the magnitude of l(φM) changes with φHφM, which reflects the change of magnetic anisotropy energy.

Here, ROTMOKE is employed to determine the anisotropy constants.[27] In order to ensure a coherent rotation process, the applied field should be at least larger than the maximum coercive field in the plane. Figures 3(a) and 3(b) show the direction of the magnetization varying with the direction of the applied field and the angular dependence of the torque curve at different strengths of field. At H = 400 Oe, the direction of the magnetization shows four sudden jumps so that the torque curve is broken down into four intervals because of the four occurrences of 90° domain wall displacement process. Further increasing the applied field to 1000 Oe, the domain wall displacement seems to disappear, and the direction of the magnetization changes with the direction of the field as a result of the coherent rotation process. The magnetic anisotropy constants can be obtained by fitting Eq. (2) to the torque curve under an applied field of 1000 Oe. Using the bulk value of 1.76 × 106 A/m as the saturation magnetization of Fe, the obtained results are Ku1 = 1.6 × 103 J/m3, Ku2 = 7.5 × 103 J/m3, and K1 = 5.6 × 104 J/m3.

Fig. 3. (a) Variations of magnetization direction with applied field and (b) variations of normalized torque l(φM) with φM at fields of 400 Oe, 600 Oe, and 1000 Oe, respectively.

For fitting the experimental data of the longitudinal Kerr signal, we consider both the linear and the quadratic magneto–optic effect. The formula of the detected normalized light intensity, which is obtained based on the proportional relationship between the output current and the squared modulus of the electric field parallel to the analyzer transmission axis, is as follows:[28]

In Eq. (3), ml = Ml/Ms, mt = Mt/Ms, with Ml and Mt being the magnetization components parallel and perpendicular to the applied field respectively. is tenable for the magnetization in the plane; superscript or subscript l and t represent the longitudinal and transverse effect; ,,, and , are the complex Fresnel coefficients of light reflected from magnetic films;[28,29] c.c is the complex conjugate of the expression and * denotes the complex conjugate of Fresnel coefficients; θp and θa are the orientations of the polarizer and the analyzer with respect to the incident light plane; and are independent of magnetization. Terms and both include only a magnetic term proportional to QMt/Ms, where Q is the Voigt magneto–optical parameter and proportional to the magnetization in the ferromagnet.[28] Terms and are proportional to QMl/Ms. Term corresponds to the sum of terms proportional to QMl/Ms and .

According to the geometrical relationship, Ml = Ms cosθ and Mt = Ms sin θ are obtained, where θ = φHφM is the angle of the positive direction of the applied field with respect to the magnetization vector. Therefore, the expression of normalized detectable light intensity in the longitudinal configuration can be rewritten as

where α, β, and γ are parameters depending on the optic geometry (polarization angle of the incident laser, analyzer angle and incident angle), and the dielectric tensor of the material.[17,28] For a given optic geometry, α, β, and γ are constants.

According to the minimum energy principle, the orientation of the magnetization φM is always changing with applied field H. The calculated anisotropy constants (Ku1, Ku2, K1) are substituted into Eq. (1) with different H values for given φH so as to obtain the value of θ. Then, we substitute θ into Eq. (4) to calculate the value of normalized light intensity and consequently obtain the hysteresis loops with changing magnetic field. The valves of α = 0.9, β = − 0.06822, γ = − 0.2847 show the best fitting of the simulation to the experimental MOKE hysteresis loops for detecting the longitudinal magnetization component.

As illustrated in Fig. 4, the coherent rotation of the magnetization is interrupted by the domain wall displacement at the coercive field. The type of domain wall displacement depends on the direction of the magnetic field. When rotating the field at φH = − 1°, the magnetization reversal process takes place via the “l-jump” process (180° domain wall displacement). At φH = 44° and φH = 46°, the magnetization reversal processes occur via the “2-jump” process (90° domain wall displacement). Besides, the branches of hysteresis loops at φH = 44° and φH = 46° appear to be asymmetric apparently due to the transverse magnetization component contribution to the light intensity that is used for detecting the longitudinal magnetization component. Since the value of α is much larger than those of β and γ, the contribution mainly comes from the quadratic magneto–optic term Ml Mt. When the field is applied closely to the hard axis, the magnetization aligns parallelly to the direction of the applied field at the saturation field. The transverse magnetization component Mt (sin θ) is almost equal to zero. As the field decreases towards the first switching field, the magnetization presents the coherent rotation to the perpendicular direction of the field after a 90° domain wall displacement process. In the coherent rotation process, Ml Mt will acquire a maximum value so that asymmetry originating from the conspicuous quadratic magneto–optic effect can be observed in the hysteresis loop. However, when the field is applied closely to the magnetic easy axis, the magnetization reversal process is governed by the 180° domain wall displacement. In the whole magnetization reversal process, the contribution of Ml Mt is small and symmetric hysteresis loops are observed.

Fig. 4. Simulated hysteresis loops at φH = − 1° (a), 44° (b), and 46° (c), with the simulating parameters being Ms = 1.76 × 103 Gs (1 Gs = 10−4 T), H = 1000 Oe, Ku1 = 1.6 × 103 J/m3, Ku2 = 7.5 × 103 J/m3, K1 = 5.6 × 104 J/m3, α = 0.9, β = − 0.06822, γ = − 0.2847. α, β, and γ are the coefficients of the magneto–optic terms.
4. Conclusions

The magnetization reversal process of single crystal Fe/MgO (001) thin film is investigated by MOKE. The results show that when the magnetic field is applied along different directions of Fe (001) film, coherent rotation can be interrupted by 90° or 180° domain wall displacement. Based on the longitudinal and transverse hysteresis loops, the magnetization orientation of the Fe (100) plane can be confirmed. Besides, we reveal that the quadratic magneto–optic effect, which comes mainly from the term of ML MT, is not negligible in the Fe/MgO (001) thin film. Therefore, it needs to be identified in analyzing Kerr signal data.

Reference
1Hu J MYang TWang JHuang HZhang JChen L QNan C W 2015 Nano Lett. 15 616
2Tsai J LTzeng H TLin G B 2010 Appl. Phys. Lett. 96 032505
3Wang ZSun KTong WWu MLiu MSun N X 2010 Phys. Rev. 81 064402
4Ye JHe WWu QHu BTang JZhang X QChen Z YCheng Z H 2014 J. Appl. Phys. 115 123910
5Ruiz-Feal IMoore T ALopez-Diaz LBland J A C 2002 Phys. Rev. 65 054409
6Ye JHe WWu QHu BTang JZhang X QChen Z YCheng Z H 2014 Appl. Phys. Lett. 105 102406
7Liu H LHe WDu H FFang Y PWu QZhang X QYang H TCheng Z H 2012 Chin. Phys. 21 077503
8Costa-Krämer J LMenéndez J LCebollada ABriones FGarcía DHernando A 2000 J. Magn. Magn. Mater. 210 341
9Cowburn R PGray S JFerreé JBland J A CMiltat J 1995 J. Appl. Phys. 78 7210
10Zhan Q FVandezande STemst KVan Haesendonck C 2009 New J. Phys. 11 063003
11Jia Y JHe WYe JHu BChen Z YGao Y HZhang X QYang H TCheng Z H 2014 Chin. Phys. 23 017502
12Qiu Z QBader S D 2000 Rev. Sci. Instrum. 71 1243
13Wang Z HCristiani GHabermeier H UBland J A C 2005 Phys. Rev. 72 054407
14Yang Z JScheinfein M R 1993 J. Appl. Phys. 74 6810
15Liang J HCao J ZLi J XWu Y Z 2015 J. Appl. Phys. 117 17E129
16Osgood Iii R MBader S DClemens B MWhite R LMatsuyama H1998J. Magn. Magn. Mater.182297
17Osgood R MClemens B MWhite R L 1997 Phys. Rev. 55 8990
18Postava KJaffres HSchuhl AVan Dau F NGoiran MFert A 1997 J. Magn. Magn. Mater. 172 199
19Buchmeier MSchreiber RBürgler D ESchneider C M 2009 Phys. Rev. 79 064402
20Appino CValsania MBasso V 2000 Physica B: Conden. Matter 275 103
21Yan S SSchreiber RGrünberg PSchäfer R 2000 J. Magn. Magn. Mater. 210 309
22Chen GLi JLiu F ZZhu JHe YWu JQiu Z QWu Y Z 2010 J. Appl. Phys. 108 073905
23Li JJin ESon HTan ACao W NHwang CQiu Z Q 2012 Rev. Sci. Instrum. 83 033906
24Cheng Z HHe WZhang X QSun D LDu H FWu QYe JFang Y PLiu H L 2015 Chin. Phys. 24 077505
25Zhan Q FVandezande STemst KVan Haesendonck C 2009 Phys. Rev. 80 094416
26Daboo CBland J A CHicken R JIves A J RBaird M JWalker M J 1993 Phys. Rev. 47 11852
27Mattheis RQuednau G 1999 J. Magn. Magn. Mater. 205 143
28Florczak J MDahlberg E D 1990 J. Appl. Phys. 67 7520
29Morley N AGibbs M R JAhmad EWill IXu Y B 2006 J. Magn. Magn. Mater. 300 436