International
Tables for Crystallography Volume C Mathematical, physical and chemical tables Edited by E. Prince © International Union of Crystallography 2006 
International Tables for Crystallography (2006). Vol. C, ch. 2.9, pp. 126146
https://doi.org/10.1107/97809553602060000584 Chapter 2.9. Neutron reflectometryNeutron reflectivity is a powerful tool for the study of surface layers and interfaces. Topics covered in this chapter include: the theory of elastic specular neutron diffraction; polarized neutron reflectivity; surface roughness; experimental methodology; resolution in real space; and applications of neutron reflectometry. 
The neutron reflectivity of a surface is defined as the ratio of the number of neutrons elastically and specularly reflected to the number of incident neutrons. When measured as a function of neutron wavevector transfer, the reflectivity curve contains information regarding the profile of the inplane average of the scatteringlength density (or simply scattering density) normal to that surface. The concentration of a given atomic species at a particular depth can then be inferred. Although similar information can often be extracted from Xray reflectivity data, an additional sensitivity is gained by using neutrons in those cases where elements with nearly the same number of electrons or different isotopes (especially hydrogen and deuterium) need to be distinguished. Furthermore, if the incident neutron beam is polarized and the resultant polarization of the reflected beam analysed, it is possible to determine, in both magnitude and direction, the inplane magnetic moment depth profile. This latter capability is greatly facilitated by the simple correlation of the relative counts of neutron spinflip and nonspinflip scattering to magnetic moment orientation and by the relatively high instrumental polarization and spinflipping efficiencies possible with neutrons. These properties make neutron reflectivity a powerful tool for the study of surface layers and interfaces.
Consider the glancing (smallangle) reflection of a neutron plane wave characterized by a wave vector from a perfectly flat and smooth surface of infinite lateral extent, as depicted schematically in Fig. 2.9.2.1 . Although the density of the material can, in general, vary as a function of depth [along the direction (z) of the surface normal], it is assumed that there are no inplane variations of the density. If the scattering is also elastic, so that the neutron neither gains nor loses energy (i.e.  k_{i}  =  k_{f}  = k = 2π/λ, where the subscripts i and f signify initial and final values, respectively, and λ is the neutron wavelength), then the component of the neutron wave vector parallel to the surface must be conserved. In this case, the magnitude of the wavevector transfer is Q = Q =  k_{f} − k_{i}  = 2ksin() = 2k_{z}, where the angles of incidence and reflection, , are equal, and the scattering is said to be specular. Since at low values of Q the neutrons are strongly scattered from the surface (i.e. the magnitude of the reflectivity approaches 1), the neutron wave function is significantly distorted from its freespace planewave form. The first Born approximation normally applied in the description of highQ crystal diffraction is therefore not valid for the analysis of lowQ reflectivity measurements, and a more accurate, dynamical treatment is required.
Because the inplane component of the neutron wave vector is a constant of the motion in the specular elastic reflection process described above, the appropriate equation of motion is the onedimensional, timeindependent, Schrödinger equation (see, for example, Merzbacher, 1970) where ψ is the neutron wave function [which in free space is proportional to where is the magnitude of the z component of the neutron wave vector in vacuum]. If the infinite planar boundary from which the neutron wave reflects separates vacuum from a medium in which the neutron potential energy is , conservation of the neutron's total energy requires that where is Planck's constant divided by 2π and m is the neutron mass. If 2π/Q has a magnitude much greater than interatomic distances in the medium, then the medium can be treated as if it were a continuum. In this limit, the potential energy can be expressed as (see, for example, Sears, 1989)where , i represents the ith atomic species in the material, is the number density of that species and is the coherent neutron scattering length for the ith atom (which is in general complex if absorption or an effective absorption such as isotopic incoherent scattering exists; magnetic contributions are not accounted for here but will be considered below). The quantity is the effective scattering density. Substituting the expression for given in equation (2.9.2.3) into (2.9.2.2) yields In order to calculate the reflectivity, continuity of the wave function and its first derivative (with respect to z) are imposed. These boundary conditions are a consequence of restrictions on current densities required by particle and momentum conservation. In general, given a sample with layers of varying potentials where the boundaries of the jth layer are at and and the potential, is constant over that layer, it can be shown that (see, for example, Yamada, Ebisawa, Achiwa, Akiyoshi & Okamoto, 1978) where and is the magnitude of the neutron wave vector in the jth layer [equation (2.9.2.4)].
The first equality in (2.9.2.5) relates the wave function at one boundary within the jth layer to the next boundary within the jth layer, whereas the second equality represents the continuity of the wave function and its derivative across the boundary between the jth and (j + 1)th layers. When a neutron plane wave is incident on a multilayer sample, we can take the incident amplitude as unity, set up the coordinate system to have z = 0 at the air/sample interface, and write the wave function in air as the sum of the incident and reflected waves, and the wave function in the substrate as a purely transmitted wave, where is the magnitude of the z component of the neutron wave vector in the substrate. By combining equations (2.9.2.5) through (2.9.2.8), we obtain a working equation for calculating the reflectivity: where total film thickness. The experimentally measured reflection and transmission coefficients R ^{2} and T ^{2} can be computed from (2.9.2.9). The procedure outlined above can be applied in piecewise continuous fashion to arbitrary, smooth potentials, ρ(z), which are approximated to any desired degree of accuracy by an appropriate number of consecutive rectangular slabs, each having its own uniform scattering density, , and thickness, , as depicted in Fig. 2.9.2.2 .
If , then becomes imaginary in the substrate, and total external reflection occurs. In addition, for a single layer deposited on the substrate, the reflectivity will oscillate with a periodicity characteristic of the layer thickness. Fig. 2.9.2.3 compares the ideal Fresnel reflectivity corresponding to an infinite silicon substrate and that of a 1000 Å nickel film deposited on silicon. For a barrier of finite thickness, tunnelling phenomena can also be observed (see, for instance, Merzbacher, 1970; Buttiker, 1983; Nuñez, Majkrzak & Berk, 1993; Steinhauser, Steryl, Scheckenhofer & Malik, 1980).

Neutron reflectivities calculated for an infinite Si substrate (dashed line) and 1000 Å Ni film on an Si substrate (solid line). 
With the matrix method described above, the reflectivity of any model scatteringdensity profile can be calculated with quantitative accuracy over many orders of magnitude. Unfortunately, the inverse computation of an unknown scattering density profile corresponding to a given reflectivity curve can be exceedingly difficult, in part due to the the lack of phase information on R(Q), which forces one to use highly nonlinear relations between R(Q)^{2} and ρ(z). Often, parameterized model scatteringdensity profiles are fit to the experimental data (Felcher & Russell, 1991). Recently, several authors have described modelindependent methods for obtaining ρ(z) from measured reflectivity curves (Zhou & Chen, 1993; Pedersen & Hamley, 1994; Berk & Majkrzak, 1995).
Reflectivity measurements with polarized neutrons can reveal the inplane magnetizationvector depth profile in magnetic thin films and multilayers. The interaction between the neutron and atomic magnetic moments is dependent upon their relative orientations. Two important yet simple selection rules apply in the case where the neutron polarization axis (defined by an applied magnetic field at the sample position) is perpendicular to Q. Any component of the inplane magnetization parallel to this quantization axis gives rise to nonspinflip (NSF) neutron scattering, which interferes with scattering due to the nuclear potential: any perpendicular magnetic component creates spinflip (SF) scattering, which is purely magnetic. Consequently, the atomic magnetic moment's direction can be inferred by measuring the two NSF (++, −−) and two SF (+−, −+) reflectivities (where +− refers to a reflection measurement in which the incident neutron magnetic moment is parallel to the applied field and only neutrons with their magnetic moment antiparallel to the applied field are measured, etc.), in addition to its absolute magnitude (which is proportional to a magnetic scattering density ρ_{M} = Np, where p is a magnetic scattering length). The matrix formalism described earlier to obtain the reflectivity in the nonpolarizedbeam case can be extended to treat polarized beams as well. The resulting transfer matrix is, however, in the latter instance a 4 × 4 matrix relating two spindependent reflectivities and transmissions for each of two possible incidentneutron spin states. The matrix elements are given in Felcher, Hilleke, Crawford, Haumann, Kleb & Ostrowski (1987), and more detailed discussions of the method of polarized neutron reflectometry can be found in Majkrzak (1991) and Majkrzak, Ankner, Berk & Gibbs (1994).
Up to this point, we have only considered reflection from smooth, flat surfaces. In reality, however, all surfaces have microscopic or mesoscopic imperfections such as steps, facets and rough hills and valleys. In this case, the potential must be represented by a threedimensional function instead of the simple onedimensional example discussed above. In addition, the roughness may not be confined to the outer surface or substrate, but the imperfections may propagate through several layers. This roughness at the interfaces modifies the specularly reflected beam and adds a diffuse component to the scattered beam (i.e. neutrons scattered at angles other than the incident angle). Theories based on the distortedwave, Born approximation have been developed to describe this type of scattering (Nevot & Croce, 1980; Sinha, Sirota, Garroff & Stanley, 1988; Pynn, 1992; Sears, 1993; Holy, Kubena, Ohlidal, Lischka & Plotz, 1993; de Boer, 1994) for a microscopically rough surface. These theories give results consistent with the earlier work (Nevot & Croce, 1980) for the modification to the specular scattering due to a single rough surface. The reader is referred to Sinha, Sirota, Garroff & Stanley (1988) and Pynn (1992) for a more complete discussion of diffuse scattering.
In order to fit the specular scattering from a rough surface, two simple methods have been employed. First, using the matrix method discussed above, the rough interface can be modelled as a smoothly varying scattering density approximated as a series of steps. This has the advantage that complex interfaces that are combinations of rough surfaces and intermixed layers can be approximated. The other method is to extend the results of Nevot to each successive interface while iteratively calculating the reflection amplitude. This method works well for simple interfaces of Gaussian roughness and is faster, in general, than the matrix method, since fewer calculations are needed for each interface. However, this latter technique suffers from frequently yielding unphysical answers (i.e. surface widths greater than adjacent layer thicknesses). Both of these methods are inadequate, in that there is no way to separate the effects of graded interfaces from rough surfaces. This can only be done with a simultaneous examination of both the diffuse and specular scattering.
Neutron reflectivity measurements can be carried out in two principal ways: either (1) with a monochromatic incident beam of narrow angular divergence in the plane of reflection (defined by and where λ is constant, and Q is varied by changing the glancing angle of incidence, θ, relative to the sample surface; or (2) using a pulsed polychromatic incident beam, also of narrow angular divergence at fixed , and obtaining data over a range of Q values simultaneously by performing timeofflight analysis on the reflected neutrons. For either method, the instrumental resolution is simply given as where Δ is the angular divergence of the reflected beam, and Δλ is the wavelength spread. In the case of a steadystate source, the wavelength resolution is determined by the monochromator, whereas the timing and moderator characteristics determine the wavelength resolution on a timeofflight instrument. Although the second term in equation (2.9.5.1) is standard in scattering, it has a unique characteristic, in that the angular divergence of the reflected beam determines the resolution. This is the case because the sample is a δfunction scatterer, so that the angle of the incident beam can be determined precisely by knowing the reflected angle (Hamilton, Hayter & Smith, 1994). For a more complete description of both types of neutron reflectometry instrumentation, see Russell (1990).
From Fig. 2.9.2.3, the period δQ of the reflectivity oscillation (in the region where the Born approximation becomes valid, sufficiently far away from the critical angle) is inversely proportional to the thickness t of the film. That is, 2π/(δQ) = t. Consequently, in order to be able to resolve reflectivity oscillations for a film of thickness t, the instrumental Q resolution ΔQ [from equation (2.9.5.1)] must be approximately 2π/t or smaller. With sufficiently good instrumental resolution, even the thickness of a film with nonabrupt interfaces can be accurately determined, as demonstrated by the hypothetical case depicted in Fig. 2.9.6.1 (where the instrumental resolution is taken to be perfect): an overall filmthickness difference of 2 Å (between 42 and 40 Å films) is clearly resolved at a Q of about 0.2 Å^{−1}. In practice, differences even less than this can be distinguished. Note, however, that to `see' more detailed features in the scatteringdensity profile (such as the oscillation on top of the plateau shown for the longdash profile in the inset of Fig. 2.9.6.1), other than the overall film thickness, it can be necessary to make reflectivity measurements at values of Q corresponding to 2π/(characteristic dimension of the feature).
One of the simplest, yet powerful, examples of the use of neutron reflectivity is in the study of selfdiffusion. Most techniques to measure diffusion coefficients rely on chemical and mechanical methods to measure density profiles after a sample has been annealed. Then a model for the diffusion is assumed, and the coefficients are calculated. Using standard techniques, researchers are unable to detect the movement of an atom through a sample of like atoms. However, using single bilayers of amorphous ^{10}B and ^{11}B, it was shown (Smith, Hamilton, Fitzsimmons, Baker, Hubbard, Nastasi, Hirvonen & Zocco, 1992) through neutronreflectivity measurements that the diffusion of boron in boron could be measured by studying the density profile (see Figs. 2.9.7.1 and 2.9.7.2 ) of one isotope in the other as a function of annealing time. Also, because of the sensitivity of the technique to the interfacial density profile, it was found that standard Fickian diffusion models could not explain the measured density profiles.
In order to understand interlayer coupling mechanisms, it is necessary to know what the magnetic superstructure is for a given nonmagnetic spacer layer thickness and/or applied field strength. Fig. 2.9.7.3 shows the spindependent reflectivities for a Co/Cu (111) multilayer along with the nuclear (Nb) and magnetic (Np) scatteringdensity profiles deduced from the data of Schreyer, Zeidler, Morawe, Metoki, Zabel, Ankner & Majkrzak (1993). In this specific case, the inplane ferromagnetic Co layers are themselves coupled ferromagnetically across the nonmagnetic Cu, all at a constant angle.
There are a substantial number of applications of neutron reflectometry in the study of hydrogenous films and multilayers, including diblock copolymer, surfactant, Langmuir–Blodgett, selfassembled monolayer, and lipid bilayer films. Reviews of the extensive research that has already been done have been written by Russell (1990) and Penfold & Thomas (1990). Only one specific example will be given here.
Fig. 2.9.7.4 shows neutron reflectivity data and the corresponding density profile for a Langmuir–Blodgett film composed of alternating bilayers of deuterated and hydrogenated stearic acid [after Wiesler, Feigin, Majkrzak, Ankner, Berzina & Troitsky (1995)]. Also shown in Fig. 2.9.7.4 is the scatteringdensity profile for the same sample as seen by Xrays. It is obvious that the Xrays are more sensitive to the highZ barium in the head group, whereas the neutrons are especially good at distinguishing the degree of mixing between adjacent hydrogenated and deuterated hydrocarbon tails. This is a good example of the complementary nature of neutron and Xray reflectivities.
References
Berk, N. F. & Majkrzak, C. F. (1995). Using parametric Bsplines to fit specular reflectivities. Phys. Rev. B, 51, 11296–11309.Boer, D. K. G. de (1994). Influence of the roughness profile on the specular reflectivity of Xrays and neutrons. Phys. Rev. B, 49, 5817–5820.
Buttiker, M. (1983). Larmor precession and the traversal time for tunneling. Phys. Rev. B,27, 6178–6188.
Felcher, G. P., Hilleke, R. O., Crawford, R. K., Haumann, J., Kleb, R. & Ostrowski, G. (1987). Polarized neutron reflectometer: a new instrument to measure magnetic depth profiles. Rev. Sci. Instrum. 58, 609–619.
Felcher, G. P. & Russell, T. P. (1991). Editors. Physica (Utrecht), B173, 1–210.
Hamilton, W. A., Hayter, J. B. & Smith, G. S. (1994). Neutron reflectometry as optical imaging. J. Neutron Res. 2, 1–19.
Holy, V., Kubena, J., Ohlidal, I., Lischka, K. & Plotz, W. (1993). Xray reflection from rough layered systems. Phys. Rev. B, 47, 15896–15903.
Majkrzak, C. F. (1991). Polarized neutron reflectometry. Physica (Utrecht), B173, 75–88.
Majkrzak, C. F., Ankner, J. F., Berk, N. F. & Gibbs, D. (1994). Magnetic multilayers, edited by L. H. Bennett & R. E. Watson, pp. 299–354. Singapore: World Scientific.
Merzbacher, E. (1970). Quantum mechanics, 2nd ed. New York: John Wiley.
Nevot, L. & Croce, P. (1980). Caracterisation des surfaces par reflexion rasante de rayons X. Application a l'etude du polissage de quelques verres silicates. Rev. Phys. Appl. 15, 761–779.
Nuñez, V., Majkrzak, C. F. & Berk, N. F. (1993). Dynamical scattering of polarized neutrons by thin magnetic films. MRS Symp. Proc. 313, 431–436.
Pedersen, J. S. & Hamley, I. W. (1994). Analysis of neutron and Xray reflectivity data by constrained leastsquares methods. Physica (Utrecht), B198, 16–23.
Penfold, J. & Thomas, R. K. (1990). The application of specular reflection of neutrons to the study of surfaces and interfaces. J. Phys. Condens. Matter, 2, 1369–1412.
Pynn, R. (1992). Neutron scattering by rough surfaces at grazing incidence. Phys. Rev. B, 45, 602–612.
Russell, T. P. (1990). Xray and neutron reflectivity for the investigation of polymers. Mater. Sci. Rep. 5, 171–271.
Schreyer, A., Zeidler, T., Morawe, C., Metoki, N., Zabel, H., Ankner, J. F. & Majkrzak, C. F. (1993). Spin polarized neutron reflectivity study of a Co/Cu superlattice. J. Appl. Phys. 73, 7616–7621.
Sears, V. F. (1989). Neutron optics. Oxford University Press.
Sears, V. F. (1993). Generalized distortedwave Born approximation for neutron reflection. Phys. Rev. B, 48, 17477–17485.
Sinha, S. K., Sirota, E. B., Garroff, S. & Stanley, H. B. (1988). Xray and neutron scattering from rough surfaces. Phys. Rev. B, 38, 2297–2311.
Smith, G. S., Hamilton, W., Fitzsimmons, M., Baker, S. M., Hubbard, K. M., Nastasi, M., Hirvonen, J.P. & Zocco, T. G. (1992). Neutron reflectivity study of thermallyinduced boron diffusion in amorphous elemental boron. SPIE Proc. Ser. 1738, 246–253.
Steinhauser, K. A., Steryl, A., Scheckenhofer, H. & Malik, S. S. (1980). Observation of quasibound states of the neutron in matter. Phys. Rev. Lett. 44, 1306–1309.
Wiesler, D. G., Feigin, L. A., Majkrzak, C. F., Ankner, J. F., Berzina, T. S. & Troitsky, V. I. (1995). Neutron and Xray reflectivity study of Ba salts of alternating bilayers of deuterated and hydrogenated stearic acid. Thin Solid Films, 266, 69–77.
Yamada, S., Ebisawa, T., Achiwa, N., Akiyoshi, T. & Okamoto, S. (1978). Neutronoptical properties of a multilayer system. Annu. Rep. Res. React. Inst. Kyoto Univ. 11, 8–27.
Zhou, X.L. & Chen, S.H. (1993). Modelindependent method for reconstruction of scatteringlengthdensity profiles using neutron or Xray reflectivity data. Phys. Rev. E, 47, 3174–3190.