research papers\(\def\hfill{\hskip 5em}\def\hfil{\hskip 3em}\def\eqno#1{\hfil {#1}}\)

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775

Magnetism in Nickel and Synchrotron Beam Polarization Studied by X-ray Diffraction

aXMaS CRG, ESRF, 38043 Grenoble CEDEX, France, bDepartment of Physics, University of Warwick, Coventry CV4 7AL, UK, cDepartment of Physics, University of Liverpool, Liverpool L69 7ZE, UK, and dDepartment of Physics, University of Keele, Keele, Staffs ST5 5BG, UK
*Correspondence e-mail: csmc@spec.warwick.ac.uk

(Received 9 December 1997; accepted 15 December 1997)

The ratio of the magnetic to the charge form factors of nickel has been determined by white-beam X-­ray diffraction. The measurements were made on the new UK magnetic scattering beamline (XMaS) on a dipole source at the ESRF. The data comprise the three (h,h,0) reflections (4,4,0), (6,6,0) and (8,8,0) and the seven high-order (h,0,0) reflections (6,0,0) to (18,0,0), which doubles the range of wavevectors compared to previous studies. The data have been analysed using Hartree–Fock free-ion wave functions and core electron polarization effects were included. The results support the interpretation of neutron data obtained at lower momentum transfer for the eg and t2g orbital occupancies. The polarization of the dipole source is deduced to vary from 99.88 to 99.83% between 5 and 15 keV, respectively. This high value makes it an extremely suitable source for studies of ferromagnetism.

1. Introduction

As explained below, the degree of polarization of the source is a critical parameter in studies of ferromagnetism. It has been measured in the present study of 3d magnetism in ferromagnetic nickel using a white-beam technique, which extends the measurements of magnetic form factors to higher momentum transfer than any reported in previous neutron or X-ray investigations using X-ray techniques.

The magnetism of the 3d transition metals is thought to be well described by a band model with the exchange interaction causing a difference in the occupancy of the 3d sub-bands. In nickel this leads to a magnetic moment of 0.57 μB per atom at room temperature. There is evidence from polarized neutron measurements three decades ago (Mook, 1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]) to support this picture and also to suggest that there is some negative spin polarization of the 4s band. Those neutron measurements were also able to determine that there was a deviation from spherical symmetry caused by crystal field effects on the charge density of the spin-polarized 3d electrons. The pioneering neutron experiments on iron (Shull & Yamada, 1962[Shull, C. G. & Yamada, Y. (1962). J. Phys. Soc. Jpn Suppl. 7(8), 1-11.]) and nickel (Mook, 1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]) could only be matched by complementary X-ray experiments after the advent of second-generation synchrotron sources which provided high-intensity highly polarized and stable beams of X-rays. The first X-ray form-factor measurements were made at the Synchrotron Radiation Source, Daresbury Laboratory, on iron (Laundy et al., 1991[Laundy, D., Collins, S. P. & Rollason, A. J. (1991). J. Phys. Condens. Matter, 3, 369-375.]; Collins et al., 1992[Collins, S. P., Laundy, D. & Rollason, A. J. (1992). Philos. Mag. B65, 37-46.]) and nickel (Zukowski et al., 1992[Zukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300-310.]).

As discussed below, measurements of this type will yield better results if performed at a low-emittance third-generation source capable of producing highly polarized X-rays, such as the European Synchrotron Radiation Facility (ESRF). The improvement in the polarization of the radiation is a direct consequence of the reduction in the size of the electron beam. The polarization of the X-rays on a second-generation machine is typically 97.0% for a dipole-magnet source, which is much lower than the value reported below for the XMaS beamline at the ESRF. Another advantage of a third-generation source such as the ESRF is the larger critical energy of the X-ray spectrum, which allows higher-order reflections to be studied by the white-beam method. With these advantages in mind, measurements have been made on the XMaS beamline, a new dipole beamline at the ESRF dedicated to magnetic and high-resolution diffraction and operated jointly by the Universities of Liverpool and Warwick for the UK scientific community.

1.1. Principles of the white-beam method

The theory of the scattering of X-rays is given, for example, by Blume (1985[Blume, M. (1985). J. Appl. Phys. 57, 3615-3618.]) and Lovesey (1987[Lovesey, S. W. (1987). J. Phys. C, 34, 5625-5639.]) and detailed in the book by Lovesey & Collins (1996[Lovesey, S. W. & Collins S. P. (1996). X-ray Scattering and Absorption by Magnetic Materials. Oxford University Press.]). The strongest mechanism for coherent scattering of X-rays from bound electrons is Thomson scattering from the electron charge density, ρ(r). The scattering amplitude for a reflection with scattering vector K is proportional to the form factor, f(K), which is given by the formula

[f({\bf K})=\textstyle\int\rho({\bf r})\exp(-i{\bf K.r})\,{\rm d}^3{\bf r}.\eqno(1)]

The magnetic scattering amplitude, which is determined by the interaction between the electronic spin and orbital moments and the electromagnetic field of the X-rays, is weaker than Thomson scattering by a factor EX/mc2, where EX is the energy of the diffracted photons and mc2 is the rest energy of the electron (512 keV). Furthermore, only a small fraction of the electrons have unpaired spins. The magnetic scattered amplitude is proportional to the magnetic form factor f mag(K), the Fourier transform of the charge density, m(r), of the electrons with unpaired moments,

[f^{\rm mag}({\bf K})=\textstyle\int m({\bf r})\exp(-i{\bf K.r})\,{\rm d}^3{\bf r}.\eqno(2)]

The nature of the interaction is such that the magnetic amplitude is phase shifted by π/2 compared with the Thomson amplitude. In ferromagnetic nickel the charge and magnetic lattices are identical and so the two types of reflection are coincident. Thus, by using an incident beam with complex (circular) polarization, there is coupling between the Thomson and the magnetic amplitudes. The measured intensity is then first order in the weak magnetic amplitude and the magnetic signal can be isolated from the charge signal by reversing the moment direction in the sample using an applied magnetic field to change the sign of the cross term.

The experimental geometry is chosen to maximize the ratio of the magnetically modulated term to the Thomson intensity (see, for example, Zukowski et al., 1992[Zukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300-310.]). The sample is aligned to diffract through 90° in a horizontal plane to take advantage of the polarization dependence of the Thomson scattering and the linear polarization of the synchrotron radiation in the orbital plane of the storage ring, i.e. to diminish the charge scattering. The incident beam is not precisely in the orbital plane but at a small angle above or below it in order that the beam contains both linearly and circularly polarized components. The moments are oriented by an applied magnetic field parallel to the scattered beam. If I+ and I are the measured intensities for the field applied first parallel and then antiparallel to the scattered beam, the X-ray `flipping ratio', R, is defined as

[R=(I_+-I_-)/(I_++I_-),\eqno(3)]

and in this geometry it is related to the charge f(K) and magnetic form factors f mag(K) (Laundy et al., 1991[Laundy, D., Collins, S. P. & Rollason, A. J. (1991). J. Phys. Condens. Matter, 3, 369-375.]) by the equation

[R=(E_X/mc^2)[P_C/(1-P_L)]\;f^{\rm mag}({\bf K})/f({\bf K}),\eqno(4)]

where PC and PL are the degrees of circular and linear polarization, respectively, of the incident X-ray beam. The ratio R depends on the polarization of the incident beam through the factor fP = PC/(1 − PL). The polarization of the synchrotron radiation varies with the take-off angle, φ, above or below the orbital plane. By varying φ the factor fP can be maximized. Fig. 1[link] shows PC, PL and fP, plotted as a function of φ, as deduced from the analysis described below. The optimum take-off angle for the XMaS dipole magnet beamline, D28, at the ESRF is seen to be φ ≃ 3 µrad, which corresponds to a height at the sample position, 25 m from the source, of only approximately 100 µm above or below the beam centre.

[Figure 1]
Figure 1
(a) The linear polarization, PL, calculated from measurements made on the XMaS beamline; (b) circular polarization, PC; (c) the factor fP = PC/(1 − PL), which enhances the ratio of magnetic to charge scattering in the white-beam diffraction method.

2. Experimental details

The experimental method follows that of previous experiments (Laundy et al., 1991[Laundy, D., Collins, S. P. & Rollason, A. J. (1991). J. Phys. Condens. Matter, 3, 369-375.]; Collins et al., 1992[Collins, S. P., Laundy, D. & Rollason, A. J. (1992). Philos. Mag. B65, 37-46.]; Zukowski et al., 1992[Zukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300-310.]). The fact that the Bragg angle is fixed (2θB = 90° in the horizontal diffraction plane) means that white beam must be used. An intrinsic germanium semiconductor detector with sufficient resolution (∼200 eV at 6 keV, ∼400 eV at 60 keV) to separate the Bragg peaks from each other and from fluorescence lines was used. The sample was a single crystal of nickel with dimensions 3 × 3 × 4 mm with four (110) and two (100) polished rectangular faces. It was mounted on a goniometer which, in turn, sat on a rotary table. The sample was aligned by taking Laue photographs of the required reflection in back reflection and then rotated horizontally through 45° to set the scattering angle to 90°. The output from the detector was amplified and fed into an analogue-to-digital converter and then to a multichannel analyser. An electromagnet with a hole bored through one of the poles provided a 0.4 T field to magnetize the sample in the direction of the diffracted beam. The magnetic field was reversed under computer control following an asynchronous cycle with a basic period of 1 s. During each period the data from the detector were accumulated and stored separately for each of the reflection orders being studied.

The experiment was set up in the optics hutch of the XMaS beamline. The X-ray energy range available for the experiment was determined by the critical energy of the dipole source, which is 9 keV. In order to take advantage of the high beam polarization it was necessary to locate the beam centre accurately and then to translate the apparatus a known distance above or below the beam centre. It was also necessary to know the linear and circular polarization of the beam at the chosen position. To this end the entire apparatus was mounted on a table with a motorized vertical translation. The beam centre was determined by scanning the apparatus vertically through a range of approximately −250 µm to +250 µm about the beam centre and measuring the diffracted intensities, the minimum corresponding to the orbital plane. The variation of intensity with the azimuthal angle allowed the electron beam source to be characterized so that the X-ray polarization could be calculated. The multichannel analyser recorded the reflection set (h,h,0) or (h,0,0) depending on the orientation of the nickel crystal. Corrections were made for the deadtime of the detector.

The low-order nickel reflections were not measured in this experiment for two reasons: (i) they had been studied by Zukowski et al. (1992[Zukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300-310.]) and (ii) this investigation was made as part of the commissioning of the optics hutch at a time when a long air path was unavoidable. This meant that these low-energy reflections were in any case severely attenuated in air. Other details of the measurement technique have been reported by Zukowski et al. (1992[Zukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300-310.]).

3. Results

3.1. Beam polarization

For a single electron, the synchrotron radiation flux and polarization may be calculated using standard formulae (Schwinger, 1949[Schwinger, J. (1949). Phys. Rev. 75, 1912-1925.]). In a real electron beam containing many electrons there will be a range of electron trajectories such that at any point in the orbit the beam has a finite size and a range of velocities. It has been shown by Laundy (1990[Laundy, D. (1990). Nucl. Instrum. Methods, A290, 248-253.]) that the modification to the formulae describing the emitted synchrotron radiation can be accomplished by the introduction of a single parameter, σγ, which characterizes the distribution of electron trajectories at a given point in the orbit.

In this experiment σγ was determined by measuring the Bragg intensity, I′(φ, σγ), as a function of the angle above, and below, the orbital plane for four of the strongest diffraction peaks,

[I^\prime(\varphi,\sigma_ \gamma)=A[1-P_L(\varphi, \sigma_\gamma)]\,I(\varphi, \sigma_\gamma),\eqno(5)]

where A is a constant for a given reflection and φ is the azimuthal angle. The functions PL(φ, σγ) and I(φ, σγ) involve Bessel functions (for details, see Laundy, 1990[Laundy, D. (1990). Nucl. Instrum. Methods, A290, 248-253.]). The X-ray linear polarization, PL, peaks close to unity on the beam orbit, i.e. when φ = 0, and hence the Bragg intensity is at a minimum. This locates accurately the beam centre. The sharpness of the minima depends strongly on σγ which can therefore be determined to a high degree of accuracy by fitting I′(φ, σγ) to the measured Bragg intensity with A and σγ as variables. We found σγ = 33.4 ± 0.6 nrad. From this result, PL and PC, together with their angle and energy dependence, can be determined from the formulae developed by Laundy (1990[Laundy, D. (1990). Nucl. Instrum. Methods, A290, 248-253.]).

The linear polarization is shown in Fig. 1[link](a) for an angular range of 6 µrad above and below the orbital plane and a photon energy of 10 keV which is close to the critical energy of the dipole magnet source. Its calculated value is 99.85 (1)% at the beam centre. The corresponding value at 5 keV is 99.88 (1)% and at 15 keV it is 99.83 (1)%. The high precision of this result arises because it is essentially 1 − PL(φ, σγ) that is measured and consequently the random error in PL is no more than ±0.01%. The value is dependent on the electron beam parameters which determine σγ. At the ESRF these are specified to remain within ±10% of their nominal value. The corresponding change in the linear polarization, PL, on orbit is only ±0.01%. The extremely high degree of linear polarization determined here is crucial in magnetic form factor measurements. Fig. 1[link] also shows the angular variation of PC and the factor fP = PC/(1 − PL) which enhances the ratio of magnetic to charge form factors in (4)[link].

3.2. Magnetic form factor

The original analysis of neutron measurements of the magnetic form factor of nickel (Mook, 1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]) used a model in which the magnetization density is derived from a free-ion Ni++ density plus a delocalized reverse-polarized 4s spin density. The Ni 3d orbitals are split by the cubic crystalline field into triply degenerate, t2g, and doubly degenerate, eg, orbitals. The occupancy ratio, r, of eg to t2g orbitals determines the degree of asphericity of the 3d charge cloud. When r = 2/5, the d orbitals are spherically symmetric. The total magnetic form factor, fmag(K), can be expressed in terms of the 3d spin form factor, fS(K), the 3d orbital form factor, fL(K), and the magnetic form factor produced by polarization of the core electrons, fcore(K). In this treatment we use form factors which are normalized at |K| = 0 to the total magnetic moment as opposed to the unitary neutron form factors which are normalized to unity at K = 0,

[f^{\rm mag}({\bf K})=1+\alpha[1-\delta({\bf K})]\;f^S({\bf K})+f^L({\bf K})+f^{\rm core}({\bf K}),\eqno(6)]

where the delta function, δ(K), insures that the uniform 4s density only contributes at K = 0; the fraction of the total moment derived from the reverse polarized spin density is represented by the parameter α.

Weiss & Freeman (1958[Weiss, R. J. & Freeman, A. J. (1958). J. Phys. Chem. Solids, 10, 147-161.]) showed that the form factor can be written in terms of integrals of the form

[\langle\;j_n\rangle_{ij}=\textstyle\int\limits^\infty_0 U_i(r)U_j(r)j_n(Kr)\,{\rm d}r,\eqno(7)]

where K = 4πsinθ/λ, Ui(r) is the radial part of the wave function for the ith orbital and jn(Kr) are the nth-order spherical Bessel functions. The ratio of the spin to orbital moments is determined by the spectroscopic splitting factor, g, which in nickel is 2.2 (Scott, 1962[Scott, G. G. (1962). J. Phys. Soc. Jpn Suppl. 7(B), 372-375.]). The 3d spin and orbital form factor can be written as

[f^S({\bf K})=\left[\langle \;j_0\rangle+A_{hkl}(5r/2-1) \langle \;j_4\rangle\right]2\mu/g, \eqno(8)]

[f^L(K)=\left(\langle \;j_0\rangle+\langle \;j_2\rangle\right) [(g-2)/g]\mu,\eqno(9)]

where the factor Ahkl is a function of the Miller indices hkl,

[A_{hkl}={h^4+k^4+l^4- 3(h^2k^2+h^2l^2+k^2l^2) \over (h^2+k^2+l^2)^2}.\eqno(10)]

In the white-beam X-ray method it is possible to measure reflections to higher order than with neutrons for two reasons. First, the X-ray cross section for magnetic scattering is proportional to the photon energy, which in the energy-dispersive method increases with the momentum transfer |K|. Second, the experiment measures the ratio of the magnetic to charge form factors. Although both decrease with K, their ratio decreases less rapidly than the absolute value of either. The highest-order X-ray reflections studied here are at |K| > 2.5 Å−1 whereas the reflections measured using neutrons all fall below 1.2 Å−1. The Ni charge form factor was calculated in the high-K regime by taking Ni++ radial wave functions published by Watson & Freeman (1960a[Watson, R. E. & Freeman, A. J. (1960a). Phys. Rev. 120, 1125-1134.],b[Watson, R. E. & Freeman, A. J. (1960b). Phys. Rev. 120, 1134-1141.]) and transforming them to give f(K), fmag(K) and fcore(K). To calculate f(K) and fcore(K), the symmetric 〈joii term was summed over the appropriate orbitals, while, for the 3d magnetic contribution, equations (6)[link], (8)[link] and (9)[link] were used. We took g = 2.2, γ = 0.19 and a = 0.19 from Mook's analysis (Mook, 1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]) of the neutron data. The experimental data and calculated values are summarized in Table 1[link], and are compared in Fig. 2[link]. As can be seen, the agreement between the calculated and measured magnetic form factors is excellent over the entire wavevector range. The graph includes three low-order form factors from the neutron study: they are in good agreement with these X-ray data.

Table 1
X-ray magnetic form factor of nickel

The calculated values are based on the Hartree–Fock wave functions of Watson & Freeman (1960a[Watson, R. E. & Freeman, A. J. (1960a). Phys. Rev. 120, 1125-1134.],b[Watson, R. E. & Freeman, A. J. (1960b). Phys. Rev. 120, 1134-1141.]).

    Experiment Calculation 
Reflection (hkl) sin θ/λ−1) fmag(K)/f(K) (%) fmag(K) f(K) fcore(K) fmag(K)
440 0.803 +0.31 (6) +0.025 (5) 8.211 −0.0027 +0.0279
600 0.851 −0.19 (3) −0.018 (3) 7.809 −0.0030 −0.0173
800 1.135 −0.46 (1) −0.035 (1) 6.313 −0.0040 −0.0335
660 1.204 −0.07 (2) −0.004 (1) 6.057 −0.0041 −0.0038
1000 1.419 −0.41 (1) −0.027 (1) 5.332 −0.0041 −0.0284
880 1.605 −0.10 (2) −0.005 (1) 4.735 −0.0038 −0.0051
1200 1.703 −0.36 (2) −0.019 (1) 4.432 −0.0036 −0.0206
1400 1.985 −0.26 (8) −0.012 (3) 3.609 −0.0029 −0.0140
1600 2.270 −0.15 (7) −0.005 (3) 2.922 −0.0021 −0.0093
1800 2.554 −0.19 (7) −0.005 (2) 2.390 −0.0015 −0.0060
[Figure 2]
Figure 2
Measured and calculated magnetic form factors for nickel at large momentum-transfer K. The filled black circles with error bars are the experimental values and lines represent the values predicted from the radial wave functions of Watson & Freeman (1960a[Watson, R. E. & Freeman, A. J. (1960a). Phys. Rev. 120, 1125-1134.],b[Watson, R. E. & Freeman, A. J. (1960b). Phys. Rev. 120, 1134-1141.]). Three neutron results are available for comparison, at low momentum transfer values; they are shown by open squares and are taken from Mook (1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]).

4. Conclusions

The degree of polarization of the ESRF dipole source on the XMaS beamline varies between 99.88 (1) and 99.83 (1)% between 5 and 15 keV, respectively. These values demonstrate that it is an excellent source for studies of ferromagnetism.

The experimental data fit the calculated magnetic form factor within experimental error and are therefore consistent with the interpretation of magnetization density in nickel given by Mook (1966[Mook, H. A. (1966). Phys. Rev. 148, 495-501.]). While there are still some questions over the validity of the form-factor fitting method (see, for example, Van Laar et al., 1982[Van Laar, B., Maniawski, F. & Kaprzyk, S. (1982). J. Phys. C, 7, 113-118.]), it is interesting that for these high-order reflections in nickel, which have not previously been measured, the X-ray data support the asphericity of the charge density with 19% of the 3d electrons in eg and 81% in t2g orbitals. They are also consistent with a negatively polarized component of the spin density amounting to 19% of the 3d moment, although these high-order reflections are less sensitive to this delocalized contribution than, for example, magnetic Compton ­scattering studies. The core polarization is significant for these high-order reflections although the 3d magnetic form factor is still dominant in this momentum transfer range.

Footnotes

Now at Daresbury Laboratory, CLRC, Warrington WA4 4AD, UK.

Acknowledgements

The authors are grateful to the EPSRC for funding the XMaS construction-phase project; these measurements formed part of the commissioning process. We thank the CLRC Daresbury Laboratory for the loan of equipment. The study also formed part of a programme supported by the European Union under Contact ERCBCHRX CT 930135. We are grateful to M. Hart, J. B. Forsyth and P. J Brown for invaluable guidance and advice.

References

First citationBlume, M. (1985). J. Appl. Phys. 57, 3615–3618. CrossRef CAS Web of Science
First citationCollins, S. P., Laundy, D. & Rollason, A. J. (1992). Philos. Mag. B65, 37–46. CrossRef Web of Science
First citationLaundy, D. (1990). Nucl. Instrum. Methods, A290, 248–253.  CrossRef CAS Web of Science
First citationLaundy, D., Collins, S. P. & Rollason, A. J. (1991). J. Phys. Condens. Matter, 3, 369–375. CrossRef CAS Web of Science
First citationLovesey, S. W. (1987). J. Phys. C, 34, 5625–5639.  CrossRef Web of Science
First citationLovesey, S. W. & Collins S. P. (1996). X-ray Scattering and Absorption by Magnetic Materials. Oxford University Press.
First citationMook, H. A. (1966). Phys. Rev. 148, 495–501. CrossRef CAS Web of Science
First citationSchwinger, J. (1949). Phys. Rev. 75, 1912–1925.  CrossRef Web of Science
First citationScott, G. G. (1962). J. Phys. Soc. Jpn Suppl. 7(B), 372–375.
First citationShull, C. G. & Yamada, Y. (1962). J. Phys. Soc. Jpn Suppl. 7(8), 1–11.
First citationVan Laar, B., Maniawski, F. & Kaprzyk, S. (1982). J. Phys. C, 7, 113–118.
First citationWatson, R. E. & Freeman, A. J. (1960a). Phys. Rev. 120, 1125–1134. CrossRef CAS Web of Science
First citationWatson, R. E. & Freeman, A. J. (1960b). Phys. Rev. 120, 1134–1141. CrossRef CAS Web of Science
First citationWeiss, R. J. & Freeman, A. J. (1958). J. Phys. Chem. Solids, 10, 147–161.  CrossRef Web of Science
First citationZukowski, E., Cooper, M. J., Armstrong, R., Ito, M., Collins, S. P., Laundy, D. & Andrejczuk, A. J. (1992). J. X-ray Sci. Technol. 3, 300–310.  CAS PubMed

© International Union of Crystallography. Prior permission is not required to reproduce short quotations, tables and figures from this article, provided the original authors and source are cited. For more information, click here.

Journal logoJOURNAL OF
SYNCHROTRON
RADIATION
ISSN: 1600-5775
Follow J. Synchrotron Rad.
Sign up for e-alerts
Follow J. Synchrotron Rad. on Twitter
Follow us on facebook
Sign up for RSS feeds