Introduction

The possibility of controlling the magnetic anisotropy of thin ferromagnetic films using a static electric field E is of great interest since it can potentially lead to magnetic random access memory (MRAM) devices which require less energy than spin-torque-transfer random access memory STT-MRAM1,2,3,4,5,6,7. Thin magnetic films with a perpendicular magnetic anisotropy (PMA) are important for applications8,9. That an interfacial internal electric field might be used to engineer such a PMA is also of great interest. Experiment10,11 has indeed shown that such a PMA might, in turn, be modified by an externally applied electric field, however the data is usually interpreted in terms of changes to the electronic contribution to magnetic anisotropy due to the surface doping induced by the applied electric field3,11,12.

The theory of the field-induced changes of the magnetic anisotropy reflecting surface doping is invariably developed in terms of band theory13,14,15,16,17. The results for both the bulk and thin films can be adequately understood in terms of second order perturbation theory18 in which the relevant matrix elements of the spin-orbit interaction are between full and empty states. Large contributions come from regions where different d-bands (almost) cross. That such crossings should be close to the Fermi surface leads to a strong doping dependence in such theories. Nakamura et al.19 pointed out that the strong negative applied field dependence of the PMA for an isolated mono-layer (ML) of Fe(001) arises directly from band splitting rather than from doping. In this case, an E perpendicular to the film breaks reflection symmetry causing a large spin-orbit splitting of d-levels near the Fermi surface. As will be explain below, despite these important theoretical developments, a clear explanation of a number of key experiments is still lacking.

Here we develop a simple analytic theory for the existence and electrical control of the PMA based upon the Rashba spin-orbit interaction20,21,22 and the single band Stoner model of magnetism. We exhibit the somewhat delicate, but very interesting, competition between the Rashba spin-orbit fields and the exchange interaction, reflecting electron correlations. This theory can potentially lead to a very large magnetic anisotropy arising from the internal electric fields Eint which exist at, e.g., ferromagnetic/metal and ferromagnetic/oxide insulator interfaces but modified by the addition of an applied electric field Eext. There is a Rashba splitting of the band structure leading to a quadratic, (Eint + Eext)2, contribution to the magnetic anisotropy, contrasting with a linear in Eext doping effect.

Results

Model

This comprises a band Stoner model with the Rashba interaction added23:

where p is the electron momentum operator, S the order parameter, σ the Pauli matrices and αR = soE the Rashba parameter proportional to ηso, which characterises the spin-orbit coupling. The electric field is taken to be perpendicular to the plane of the system and is perpendicular to and makes an angle θ to the -direction, as in Fig. 1(a).

Figure 1
figure 1

(a) The electric field is perpendicular to the ferromagnet surface while the order parameter direction , is defined by the angle θ relative to . Whatever the direction of k, the Rashba magnetic field BR of direction k × E lies in the plane. (b) The Rashba split bands of a non-magnetic metal. The two Fermi sheets emerge from a “Dirac point” near the bottom of the illustration. For the magnetic case the two Fermi sheets are disconnected. (c) For a perpendicular the electron spins make a constant angle δ to the vertical such that the projection is as in (b). The additional exchange splitting increases as E2. (d) Same but for parallel to the plane. With along the y-direction the majority and minority Fermi seas shift along the x-axis in opposite directions. The tilt of the spin relative to is no longer a constant being zero along the x-axis and a maximum along the y-axis.

Illustrative non-magnetic example

Consider the Rashba effect in non-magnetic two dimensional electron gases or surface states on noble metals, e.g., a surface state of Au. As shown in the methods, the single particle energy

as the spin quantum number σ = ±1. The momentum shift and

identified in the inset of Fig. 1(b), reflects the single particle energy gain relative to zero electric field, i.e., E = 0 and hence αR = 0. For the three dimensional problem and there is no equivalent momentum shift in kz. For the surface state of Au, ER ≈ 3.5 meV24,25 exemplifying the energy scale.

The study of such surface states and differences in chemical potentials, also helps set the scale for E. Between a metal and the vacuum, or dielectric, the E ~ 10 V/nm near the Cu (100) surface26, reflects the electron image potential. For the Ag/Cu(111) system27 an E of similar magnitude occurs at the metal interface. The difference in the (111) chemical potentials28 of Cu(4.96 eV) and Ag(4.74 eV) is reflected by a potential increase from Cu to Ag. The chemical potential of Au(5.31 eV) implies a similar potential change and E between Cu and Au, but of the opposite sign. It is the s-electrons which penetrate the core and determine the spin-orbit parameter ηso. The example Cu/Ag/Au therefore illustrates the higher/lower potential of these s-electrons in the second/third transition series relative to the similar electrons in the 3d elements.

There are many experiments24,25,29,30,31 which put in evidence the Rashba splitting in two dimensional electron gases, surface states of noble metals, bulk layered systems and e.g., of a surface state of ferromagnetic Yb.

Magnetic case - origin of the magnetic anisotropy energy

In the methods it is shown the θ dependent single particle energy, i.e., the equivalent of Eq. (2) for the magnetic case is:

The direction of the now θ dependent momentum shift changes sign as the spin index σ = ±1. These shifts also change sign with for a given σ. This “magnetic Rashba splitting” with is observed for the surface state of Yb31.

Also in the methods the Dzyaloshinskii-Moriya (DM) and pseudo-dipolar (PD) contributions to the magnetic anisotropy are highlighted by contrasting the perpendicular and parallel orientations of order parameter to the plane.

Assuming (J0S)2 > (αRkx)2 and retaining the θ-dependent terms up to the order of E2 in (4), we obtained our principal result:

for the magnetic anisotropy energy, with

where 〈 〉 denotes an average over the Fermi sea (see methods). The Rashba spin-orbit interaction produces a uni-axial anisotropy energy which, as in the Dzyaloshinskii-Moriya theory32,33,34,35, comprises a direct second order in E easy plane pseudo-dipolar interaction and an indirect contribution proportional to E2/J0S reflecting the competition between the first order in E, Rashba-Dzyaloshinskii-Moriya and exchange fields.

Competition between the Dzyaloshinskii-Moriya and pseudo-dipolar contributions

Clearly Eq. (5) implies an E2 dependent PMA results if T > J0S/2, i.e., when the DM is larger than the PD term. Taken literally, the Stoner model Eq. (1), with its quadratic dispersion, predicts the ratio of the DM and PD contributions to the PMA. The result, (see methods), depends upon the spatial dimension. In two dimensions the PD and DM terms cancel although higher order terms (OR4)) lead to a PMA while in three dimensions the DM term is −(4/5) ER cos2 θ and an in-plane magnetisation is favoured. Lastly, a two dimensional system with a highly anisotropic conductivity might be modelled as a series of parallel one dimensional chains. For chains the DM contribution −(4/3) ER cos2 θ which dominates the PD energy ER cos2 θ, appropriate when is in-plane and perpendicular to the chains. Corresponding to the hardest axis, when in-plane but rather parallel to the chains, there is neither a DM or PD contribution to the magnetic anisotropy energy. There is thereby a predicted electric field dependence of the in-plane anisotropy as seen in early experiments3, given the large compressive strain that arises in these experiments.

However, for the real problem of 3d magnets, a quadratic dispersion is not at all realistic and the crystal potential V(r) must be accounted for, see methods. For 3d elements the wavefunction ψ is well localised within the atomic sphere and the averages, e.g., and hence T, are very much increased as compared to the above naïve estimates. In reality, the DM contribution will invariably lead to a PMA.

Discussion

The resulting anisotropy energy can be very large. The work already cited24,25,26,27 on conducting but non-magnetic materials helps set the scales. The value of the scaling prefactor ER in Eq. (5) for the surface state of Au is ~ 3.5 meV or about 35 T in magnetic field units and very much larger than the typical ~ 1 T demagnetising field. If a Au film is polarised by contact with an ultra-thin ferromagnet, the second factor, 2T/J0S, in Eq. (5) for the field inside a Au surface layer can be quite large ~ 5 leading to a PMA and indeed ultra-thin Fe on Au does have such a PMA36,37. Ultra-thin ferromagnetic films in contact with, e.g., Ru, Pd, Pt and Ta, etc., also are found to have a PMA6,12,38,39.

Schematically shown in Fig. 2(a) is the potential seen by electrons in a free standing ultrathin ferromagnetic film. At the surface, the potential reflects an electron's image charge but reaches the vacuum level within a few atomic spacing. As already discussed, this results in a finite large electric field E ~ 10 V/nm at each surface but in opposite senses. Assuming an appreciable spin-orbit coupling in the interface region, this results, in turn, in a Rashba field BR which also changes sign between the two surfaces for a given momentum. Thus, for a perfectly symmetric film, the ferromagnetically polarised electrons see no average field BR. This symmetry can be broken by the application of an external electric field as shown in Fig. 2(b). The electric field is increased at one surface and decreased at the other doubling the net effect. In contrast, for this same symmetric situation, the surface charges are opposite and doping effects must cancel. Experimentally applied fields of 1 V/nm are relatively easy to achieve implying a ~ 10% change in the surface anisotropy. Experiments40 with a 1.5 nm Fe80B20 sandwiched between two MgO layers are perhaps closest to this situation although the thickness 1.5 nm and 2.5 nm of these layers are not equal. Roughly consistent with our estimate [see Fig. 2(c), case(i)], there is40 an ~ 15% symmetric contribution to the magnetic anisotropy for an applied voltage of 2 V.

Figure 2
figure 2

(a) There is an electric field E in the surface region of a ferromagnet, however for a given wave vector k, the Rashba field BR, proportional to k × E, has an opposite sign at the two surfaces and the average field is zero. (b) With a finite external field this symmetry is broken and there is a net Rashba field acting upon the electrons. (c) The gate voltage dependence of the anisotropy energy. The internal electric field causes the shift of the parabola in the lateral axis as indicated by V0 for case i). For cases ii) and iii) the internal field shift is far beyond the external field range and nearly linear E-dependence arises. (d) The symmetry is also broken for a insulator-ferromagnet-metal sandwich. Also despite the electric field being smaller at the right surface, for a suitable metal, the spin-orbit coupling is larger and hence the metal interface can still dominate the net Rashba field. (e) Here the work function is larger for the metal than for the ferromagnet and the field for that surface is reversed. Now the Rashba fields at the two surfaces add. (f), (g) Applying positive gate voltages decreases the Rashba field at the insulating surface which, for this case, causes a net increase/decrease in the average Rashba field.

Clearly the intrinsic Rashba field BR is modified when the materials adjacent to a 3d ferromagnet (F) are different. In a number of experiments an insulator I, often MgO, lies to one side and a normal metal (N), e.g., Au, Pt, Pd, Ta, or Ru, “cap” completes a tri-layer system. The potential, Figs. 2(d) and (e), will increase in passing from Fe to MgO but in passing from the F to N-layer the potential will either increase, Fig. 2(d), or decrease, Fig. 2(e). As discussed above, the relevant potential would be expected to increase, Fig. 2(d), in passing from a 3d element, e.g., Co, Fe, Ni, or Cu, to a 4d transition metal such as Ru, Pd, Ag but decrease, Fig. 2(e), for the 5d elements, e.g., Au, Pt, Ta. The latter case is particularly favourable since the intrinsic Rashba fields have the same sense and add. In addition, the 5d elements have a larger spin-orbit coupling, resulting in a larger αR and hence are more likely to produce a sizeable PMA. If the electric field decreases at the FI interface, the average Rashba field increases in the first case [Fig. 2(f)] when the effects of the surfaces tend to cancel and, as illustrated in Fig. 2(g), decreases in the second case when the inverse is true [see Fig. 2c, case (ii) and (iii)]. Experiment12 indeed shows an opposite field dependence for such systems with Pd(4d) and Pt(5d) N-layers. That the sign of the electric field contribution to the PMA reflects the N-layer whereas the field is applied to the opposite surface between F and I supports the current Rashba model. This is in stark contrast with the popular surface doping model3,12, for which the effects of surface doping are limited by the (possibly magnetically modified) Fermi-Thomas screening length. In reality12, the screening length is estimated to be much less than 1 nm and much too short for there to be an appreciable doping effect of the Pt or Pd layers that are typically distant by a few nanometers.

Simulating a large applied electric field E, the required asymmetry might be controlled in NFN tri-layers by varying in a systematic manner, at the mono-layer level, the thickness of one of the normal metal layers and by using metals with different spin-orbit couplings. In reality the effect of the substrate transmitted to and through, the bottom normal metal will imply an asymmetry even for largish N-layer thickness. Indeed the PMA surface term for Au/Fe(110)/Au(111) structures does show an non-monotonic dependence on the top Au layer thickness36. Experiments37 for Fe layers on vicinal Ag(001) and Au(001) surfaces and which undergo a symmetry breaking (5 × 20) surface reconstruction manifest an in-plane surface term reflecting this broken symmetry and which is larger for Au, with its stronger spin-orbit coupling, than for Ag.

It is predicted that the surface coercivity field Hc is proportional to (Eint + Eext)2 where Eint is the internal electric field corresponding to the zero-bias Rashba contribution to the anisotropy. Such a non-linear field dependence is observed, e.g., for the in-plane contribution for a (Ge,Mn)As/ZrO2 surface3. In other experiments10 with CoFeB/MgO/CoFeB structures there is qualitative difference between the E dependence of the anisotropy field Hc of the, “top” and “bottom”, CoFeB layers of this three layer structure, even when they have similar thicknesses. The bottom layer has a larger Hc and is roughly linear while Hc becomes highly non-linear as Hc → 0 as would be expected as Eext → −Eint.

The most direct experimental test of the model is the observation of the band splittings for a model Rashba system with a variable contact with an itinerant ferromagnet. This can result in giant magnetic anisotropy (GMA) energies. For example an ER ~ 100 meV (or ~ 1000 T) is reported in angle-resolved photoemission spectroscopy (ARPES) measurements41 on bulk BiTeI. For a thin film of this, or similar material, in contact with an itinerant ferromagnet such as Fe, a suitable exchange splitting J0S, tuned to the order of ER, might be induced and a GMA will result. ARPES performed as a function of the direction of the magnetisation m might determine both ER and the momentum dependence of the exchange splitting leading to estimates of both the PM and DM contributions and which might be directly compared with magnetisation and magnetic resonance measurements. The electrical control of such a GMA has evident important application for non-volatile memory applications. There are clearly many more complicated embodiments of such a device.

In conclusion, it is suggested that the Rashba magnetic field due to the internal electric field in the surface region of an ultra-thin ferromagnet can make an important contribution to the perpendicular magnetic anisotropy. Such surface fields might be modified by application of an applied electric field. Since the internal fields at two surfaces tend to cancel, an asymmetry between the surfaces is important. Such an asymmetry is caused by different metal and insulator caping layers. These ideas are consistent with a large number of experiments.

Methods

The non-magnetic case

This corresponds to Eq. (1) with J0 = 0. It is solved by taking the axis of quantisation to be perpendicular to the in-plane k as in Fig. 1(b). The eigenstates are eik·r|s〉 and , where the Rashba magnetic field in energy units is defined as , with μB the Bohr magneton and g the g-factor, leaving the spin state |s〉 to be determined. There are two concentric Fermi surfaces. The energy splitting , where Δ is the value for , with kF↑,↓ the Fermi wave number for the spin up/down (σ = ±1) band. For the surface state of Au, Δ ≈ 110 meV while EF ≈ 420 meV giving the ER ≈ 3.5 meV cited in the text. The magnetic case. The full Eq. (1) is solved by defining axes such that lies in the yz-plane and . The total field, which defines the axis of quantisation, . It is assumed that, for a 3d ferromagnet J0S ~ 0.5–1.0 eV and BBR < J0S, i.e., the Rashba is smaller than the exchange splitting. To second order in BBR, where and where differs in direction from by a small angle δ where tan δ ≈ αR(kx2 cos2θ + ky2)1/2/J0S. The linear in kx term, αRkx sin θ, causes a shift in Fermi sea to give the the single particle energy Eq. (4).

With perpendicular to the plane, i.e., , the exchange and Rashba fields are orthogonal and hence the net energy for a single electron Eq. (4) is

The axis of quantisation is tilted by away from the z-axis as shown in Fig. 1(c). The σ = ±1 electrons gain/lose an energy that is even in E. This arises from the competition of the Rashba field, perpendicular to , with the exchange field. Such a competition generates a second order in E contribution to the magnetic anisotropy and is identified with the Dzyaloshinskii-Moriya (DM) mechanism32,33,34,35.

Now take parallel to the y-axis, i.e., . The y-component of BR is parallel to the exchange field and is combined with the kinetic energy. The Fermi sea is shifted along the x-axis and lowered by ER as shown in Fig. 1(d). This energy gain corresponds to a pseudo-dipolar (PD) contribution to anisotropy energy35 which favours an in-plane magnetisation. On the other hand, the x-component of BR, which is perpendicular to , gives rise to a correction to the effective exchange field. The direction of the moment tilts away from the y-axis in the direction perpendicular to the wave vector by as shown in Fig. 1(d). The single particle energy, Eq. (4), is now,

where the shift k0 is the same as in Eq. (2) but only along the x-axis.

The effective exchange field in Eq. (8) is smaller than that in Eq. (7) due to the absence of a kx2 term. This indicates that the overall DM contribution favours a perpendicular while the PD term favours an in-plane . This exchange field changes sign with σ = ±1.

Evaluation of the Dzyaloshinskii-Moriya and pseudo-dipolar contributions

Needed for T in Eq. (5) are the Fermi sea averages and , determined analytically, for quadratic dispersion, in one, two and three dimensions respectively. For an isotropic system these averages are related to J0S via

which determines the ratio 2T/J0S in the principal result, Eq. (5), given in the text.

Role of the crystal potential

The effects of the crystal potential are exhibited by considering a wave function which is a linear combination of plane waves, where K are the reciprocal lattice vectors and the aK are determined by V (r). While not convenient for 3d electrons, at least in principle, such an expansion in the true, rather than crystal, momentum states is always possible. The PD contribution, ER cos2 θ is independent of the momentum k + K. However , where 〈 〉BZ is the average over the first Brillouin zone. For 3d electrons, the average and hence T, are dominated by the aK for largish K. It follows T is significantly increased with the consequences discussed in the text.