 1. Introduction
 2. Model for the magnetic microstructure of bulk ferromagnets
 3. Magnetic SANS theory of bulk ferromagnets – unpolarized neutrons
 4. Correlation function of the spinmisalignment SANS cross section
 5. Models for SH and SM
 6. Results and discussion
 7. Summary and conclusion
 References
 1. Introduction
 2. Model for the magnetic microstructure of bulk ferromagnets
 3. Magnetic SANS theory of bulk ferromagnets – unpolarized neutrons
 4. Correlation function of the spinmisalignment SANS cross section
 5. Models for SH and SM
 6. Results and discussion
 7. Summary and conclusion
 References
research papers
Smallangle neutron scattering correlation functions of bulk magnetic materials
^{a}Physics and Materials Science Research Unit, University of Luxembourg, 162A Avenue de la Faïencerie, L1511 Luxembourg, Luxembourg
^{*}Correspondence email: andreas.michels@uni.lu
On the basis of the continuum theory of micromagnetics, the correlation function of the spinmisalignment smallangle neutron scattering e.g. elemental polycrystalline soft and hard magnetic nanocomposites, nanoporous or magnetic steels) is computed. For such materials, the spin disorder which is related to spatial variations in the saturation magnetization and magnetic anisotropy field results in strong spinmisalignment scattering dΣ_{M}/dΩ along the forward direction. When the applied magnetic field is perpendicular to the incoming neutron beam, the characteristics of dΣ_{M}/dΩ (e.g. the angular anisotropy on a twodimensional detector or the asymptotic powerlaw exponent) are determined by the ratio of magnetic anisotropy field strength H_{p} to the jump ΔM in the saturation magnetization at internal interfaces. Here, the corresponding one and twodimensional realspace correlations are analyzed as a function of applied magnetic field, the ratio H_{p}/ΔM, the singleparticle form factor and the particle Finally, the theoretical results for the correlation function are compared with experimental data on nanocrystalline cobalt and nickel.
of bulk (1. Introduction
Smallangle neutron scattering (SANS) is a very popular method for investigating nanoscale structural and magnetic inhomogeneities in the bulk of materials. In most situations, SANS data are analyzed in ; Hansen, 2000; Fritz & Glatter, 2006; Hansen, 2012), which has recently been extended to allow for the analysis of twodimensional smallangle scattering patterns of oriented samples (FritzPopovski, 2013; FritzPopovski, 2015). For dilute, monodisperse and uniform particle–matrix systems, several analytical expressions for the density–density autocorrelation function or, likewise, for the distance distribution function have been derived (see e.g. Svergun & Koch, 2003); this is a well established procedure in smallangle Xray scattering and in nuclear SANS, e.g. in the analysis of polymers (Mortensen & Pedersen, 1993) or in the study of the formation of magnetic nanocrystals in glass ceramics (Lembke et al., 1999).
by fitting a particular model to the experimental SANS An alternative realspace approach to analyzing SANS data is the computation of the (auto)correlation function of the system, for instance by means of the indirect Fourier transformation technique (Glatter, 1977In the context of realspace analysis of scattering data, it is also worth mentioning the recent progress made in the computation of the magnetic pair distribution function (Frandsen et al., 2014), which is obtained via Fourier transformation of the magnetic neutron scattering This approach permits the analysis of long and shortrange magnetic correlations of a wide range of magnetic structures such as spindensity waves, spinice compounds or molecular magnets.
We have recently provided a theory of magnetic SANS of polycrystalline bulk ), which was successfully employed in order to analyze the magnetic microstructure of ironbased twophase nanocomposites (Honecker et al., 2013). In addition to nanocomposites, the theory is also applicable to the study of elemental nanoporous magnets or ferromagnetic steels; it provides information on the exchangestiffness constant, as well as on the strength and spatial structure of the magnetic anisotropy and magnetostatic field.
(Honecker & Michels, 2013Magnetic SANS of statistically isotropic bulk i.e. the magnetic SANS depends not only on the magnitude but also on the orientation of the momentumtransfer vector. The results for the Fourier coefficients of the magnetization (Honecker & Michels, 2013) demonstrate the unmistakable impact of the magnetodipolar interaction on magnetic SANS. Magnetostatics is essential for understanding the complex magnetic fielddependent angular anisotropies which may be observed on a twodimensional positionsensitive detector; these anisotropies go beyond the well known `' anisotropy of magnetic SANS. Furthermore, the classical particle–matrix concept of smallangle scattering is not adapted to the complex magnetic textures that may form inside the bulk of magnetic media [see discussion in the introduction of Michels (2014)]; for such materials, the continuum theory of micromagnetics (Brown, 1963) provides the proper theoretical framework for computing the magnetic SANS It is the purpose of this paper to provide a discussion of the predictions of our micromagnetic SANS theory in real space by calculating the correlation function of the spinmisalignment SANS cross section.
is, in contrast to nuclear SANS on such structures, highly anisotropic,The paper is organized as follows: §2 introduces the model for the magnetic microstructure of bulk which underlies our magnetic SANS theory; in §3, we summarize the main expressions for the unpolarized magnetic SANS in §4, we define the correlation function of the spinmisalignment SANS and we compare its definition with the corresponding result from nuclear SANS theory; §5 details the models for the anisotropy field and longitudinal magnetization Fourier coefficient; in §6, we discuss the results for the correlation functions and correlation lengths, and we provide a comparison with experimental data; §7 summarizes the main findings of this study.
2. Model for the magnetic microstructure of bulk ferromagnets
We consider polycrystalline statistically isotropic bulk et al., 2004; Löffler et al., 2005; Michels et al., 2008, 2009; Döbrich et al., 2012), soft magnetic twophase nanocomposites from the FINEMET (VITROPERM) or NANOPERM family of alloys (Ohnuma et al., 2000; Heinemann et al., 2000; Michels et al., 2006), NdFeBbased permanent magnets (Bick et al., 2013; Périgo et al., 2015), and magnetic steels (Coppola et al., 1998; Bischof et al., 2007; Michaud et al., 2007; Alinger et al., 2009; Bergner et al., 2013). Fig. 1(a) shows a sketch of the nuclear (grain) microstructure of such a material, and Fig. 1(b) displays qualitatively the magnetic (spin) distribution at a nearly saturating applied magnetic field.
Examples of such materials are inertgas condensed singlephase elemental (WeissmüllerOn the basis of the continuum theory of micromagnetics (Brown, 1963), we have provided (Honecker & Michels, 2013) a firstorder theory for the magnetic spinmisalignment SANS of weakly inhomogeneous bulk which accounts for spatial variations in the magnetic anisotropy and saturation magnetization. The theory, valid close to magnetic saturation, is based on the solution of the well known balanceoftorques equation,
which expresses the fact that at static equilibrium the torque on the magnetization vector field due to an effective magnetic field vanishes everywhere inside the material. The effective field
is composed of a uniform applied magnetic field , the magnetostatic field , the magnetic anisotropy field and the exchange field . The general solution of equation (1) for the transverse magnetization Fourier coefficients (in the highfield limit) is given in Appendix A. Metlov & Michels (2015) extended the firstorder theory to second order in the amplitudes of the inhomogeneities (including fluctuations in the exchange interaction), and the corresponding magnetic SANS was computed up to the third order. For the sake of a selfcontained presentation, we summarize in §3 the main results for the magnetic SANS of bulk magnetic materials.
3. Magnetic SANS theory of bulk – unpolarized neutrons
Since the spinmisalignment scattering of bulk ), halfpolarized (SANSPOL) experiments on bulk magnetic materials do not provide significantly more information regarding the spinmisalignment SANS than can already be learned from the analysis of unpolarized data; this is because the SANSPOL `spinup' and `spindown' cross sections differ essentially only by a nuclearmagnetic interference term , which is usually small and weakly field dependent as compared to the spinmisalignment SANS. In order to demonstrate the main effects, we concentrate in the following on the unpolarized magnetic SANS of bulk in the two scattering geometries that have the applied magnetic field either perpendicular or parallel to the incident neutronbeam direction (see Fig. 2). The corresponding equations for polarized SANS and, in particular, the spinflip (POLARIS) equations are given by Honecker et al. (2010) and Michels (2014).
is independent of the polarization of the incident neutron beam, it is sufficient to restrict the considerations for the correlation function to the unpolarized As discussed by Michels (20143.1. k_{0} ⊥ H_{0}
For the scattering geometry where the applied magnetic field is perpendicular to the wavevector of the incoming neutron beam (compare Fig. 2a), the elastic unpolarized SANS at scattering vector can be written as (Michels, 2014)
, where ψ is half the scattering angle and λ is the wavelength of the incident radiation, V is the scattering volume, b_{H} = 2.91 ×10^{8} A^{1} m^{1} relates the atomic to the and denote, respectively, the Fourier coefficients of the nuclear scattering length density and of the magnetization , and θ represents the angle between and ; the asterisks (^{*}) mark the complex conjugate quantity, and the atomic magnetic form factor (in the expression for b_{H}) is approximated to unity (forward scattering).
As shown by Honecker & Michels (2013), near magnetic saturation and for a weakly inhomogeneous bulk ferromagnet, can be evaluated by means of micromagnetic theory. In particular,
where
represents the nuclear and magnetic residual SANS
which is measured at complete magnetic saturation (infinite field), andis the spinmisalignment SANS S_{H} R_{H} due to perturbing magnetic anisotropy fields and a part S_{M} R_{M} related to magnetostatic fields. The micromagnetic SANS theory considers a uniform exchange interaction and a random distribution of magnetic easy axes, but takes explicitly into account variations in the magnitude of the magnetization [via the function S_{M}, see equation (8) below].
The magnetic scattering due to transverse spin components, with related Fourier amplitudes and , is contained in , which decomposes into a contributionThe anisotropyfield scattering function (in units of cm^{1} sr^{1})
depends on the Fourier coefficient of the magnetic anisotropy field, whereas the scattering function of the longitudinal magnetization (in units of cm^{1} sr^{1})
provides information on the magnitude of the magnetization jump at internal (e.g. particle–matrix) interfaces. The corresponding (dimensionless) micromagnetic response functions can be expressed as
and
where
is a (dimensionless) function. The effective magnetic field
depends on the internal magnetic field
on and on the exchange length
(M_{s}: saturation magnetization; A: exchangestiffness parameter; : demagnetizing factor; 10^{7} T m A^{1}). The θ dependence of R_{H} and R_{M} is essentially a consequence of the magnetodipolar interaction. Depending on the values of q and H_{i}, and on the ratio , a variety of angular anisotropies may be seen on a twodimensional detector (see e.g. Fig. 11 in §6.2 below) (Michels et al., 2014; Michels, 2014).
By assuming that the functions , and h depend only on the magnitude of the scattering vector, one can perform an azimuthal average of equation (4), i.e. . The resulting expressions for the response functions then read (see Fig. 3)
and
so that the azimuthally averaged total nuclear and magnetic unpolarized SANS
of a bulk ferromagnet can be written aswhere
and
3.2. k_{0} ⊥ H_{0}
For the scattering geometry where the external magnetic field is parallel to the incidentbeam direction (compare Fig. 2b), the total unpolarized SANS can be written as (Michels, 2014)
where . Using linearized micromagnetic theory, the azimuthally averaged version of equation (20) can be expressed as
where the residual SANS
explicitly readsand the spinmisalignment SANS equals
with
S_{H}(q) is given by equation (7), and we note that in this geometry does not depend on fluctuations and equals the expression for the singlephase material case (Weissmüller et al., 1999), in other words, inhomogeneities in the saturation magnetization are (for ) only contained in and not in .
4. Correlation function of the spinmisalignment SANS cross section
Before addressing the magnetic correlation functions, we will briefly recall the corresponding well known results from nuclear SANS theory (Guinier & Fournet, 1955; Porod, 1982; Feigin & Svergun, 1987). The nuclear SANS
can be expressed in terms of the autocorrelation function of the (excess) nuclear scattering length density (in units of m^{2}) as
where
and
The function denotes the socalled excess scattering length density, where is the (constant) average scattering length density, which only gives a contribution to at . The backtransform of equation (26) is
which for isotropic systems reduces to
In analogy to the above formalism, one may define the autocorrelation function of the spin misalignment as (Michels et al., 2003; Weissmüller et al., 2004; Michels, 2010)
where denotes the deviation of the local magnetization vector field from the mean magnetization . Alternatively, can be expressed as
where is the Fourier transform of . In the highfield limit, is nearly parallel to the applied magnetic field with , so that and
Note that in our theory of magnetic SANS (Honecker & Michels, 2013) the magnetization components are all considered to be real valued.
Comparison of equations (27) and (31) reveals an important difference between nuclear and magnetic scattering [besides the fact that is a scalar and a vector quantity]: while the nuclear SANS is directly proportional to the Fourier transform of , the function [being the Fourier transfrom of ] does not represent the experimentally , which, according to equations (3) and (20), is a weighted sum of the Cartesian Fourier components of the magnetization.
Therefore, we define the correlation function of the spinmisalignment SANS
as the Fourier transform of , for which we have a theory, according toThe normalized version of equation (34),
forms the basis for the calculations of the present work. We emphasize that the that is defined in this way is not an autocorrelation function, as are and . Likewise, the well known result that the evaluation of and at the origin yields, respectively, the meansquared density fluctuation (Porod invariant) and the meansquared magnetization fluctuation does not pertain to ; the integral of over
does not provide an obvious invariant of the spinmisalignment SANS.We remind the reader that at a particular applied magnetic field H_{i} can be (approximately) obtained by subtracting the total nuclear and magnetic scattering at a saturating field from the measurement of the total at the particular H_{i}.
The spinmisalignment SANS q and the direction θ of the scattering vector on the detector (see e.g. Fig. 11 in §6.2 below). The θ dependence of is a consequence of the magnetodipolar interaction – via the Fourier coefficients (Erokhin et al., 2012; Honecker & Michels, 2013; Michels et al., 2014) – and of the trigonometric functions which are explicitly contained in the [equation (3)] and are due to the dipolar nature of the neutron–magnetic interaction. The final expression for the (azimuthally) θaveraged [equation (19)] contains the averages over these Since from a practical point of view it is easier to work with onedimensional data, i.e. with , equation (35) may be simplified to
for the perpendicular scattering geometry depends on both the magnitudewhere j_{0}(x) = sin(x)/x denotes the zerothorder spherical Bessel function. Note that spherical Bessel functions are denoted with a lowercase `j', whereas Bessel functions are represented with an uppercase `J'. Equation (36), which from now on is called the `onedimensional' correlation function of the spinmisalignment SANS has the same mathematical structure as the corresponding equation (30) for nuclear SANS.
Since for statistically isotropic bulk θ) (Michels et al., 2014), equation (36) also applies to .
in the parallel scattering geometry is isotropic (independent of the angleIn a SANS experiment, only the components of the momentumtransfer vector perpendicular to the incidentbeam direction () are effectively probed, which from a mathematical point of view means that the measured ). In §6.2 below, we will also study (for ) the case of anisotropic twodimensional correlations by considering the following expression for c(y, z) (Šaroun, 2000):
already represents an average over the incidentbeam direction. For , this implies that , whereas for (compare Fig. 2Because , the c(y, z) that is computed according to equation (37) represents a projection (average) of the threedimensional correlation function c(x, y, z) along the direction of the incident neutron beam (FritzPopovski, 2013, 2015).
Equation (37) can be transformed into polar coordinates, which results in
where the angle φ specifies the orientation of in the yz plane. By introducing the nthorder Bessel function (Watson, 1966),
where n is an integer and the last equation is valid for any angle δ, we can obtain an average of over all angles φ in the detector plane:
Since the integration with respect to the angle θ can be taken analytically [compare arguments leading to equations (15) and (16)], it follows that
Equation (41) is called the averaged `twodimensional' correlation function of the spinmisalignment SANS Note that this expression differs from equation (36) which is obtained after threedimensional integration of the θaveraged .
In Appendix A, we provide a comparison between the autocorrelation function of the spin misalignment, c_{SM}(r), and the correlation functions of the spinmisalignment SANS cross sections, equations (36) and (41).
5. Models for S_{H} and S_{M}
In order to solve equation (36) [or equation (41)], we have to specify certain models for the anisotropyfield scattering function [equation (7)] and for the scattering function of the longitudinal magnetization [equation (8)] in the expression for . As outlined in §2, we consider a statistically isotropic nearly saturated bulk ferromagnet which exhibits (weak) spatial fluctuations of the saturation magnetization and the magnetic anisotropy field. For such a system, the functions h^{2} and depend only on the magnitude q of the momentumtransfer vector . Furthermore, we assume a monodisperse scattering system and that both functions h^{2}(q) and can be written as the product of the same singleparticle form factor P(q) and S(q) (Pedersen, 1997), i.e.
and
where V_{p} is the particle volume. Later on in the calculations, we will use (for illustration purposes) the Percus–Yevick hardsphere for S(q) (Kinning & Thomas, 1984) and (unless stated otherwise) the sphere form factor for P(q),
where j_{1}(x) denotes the spherical Bessel function of first order. Any other particle form factor or may be straightforwardly implemented (see below). We also note that the characteristic structure sizes of h^{2} and need not be identical; these are related, respectively, to the spatial extent of regions with uniform magnetic anisotropy field and saturation magnetization.
Under these assumptions (same size and shape), h^{2}(q) and differ only by constant prefactors, i.e. the magnitude H_{p} of the mean magnetic anisotropy field and the jump of the magnitude of the magnetization at internal interfaces. In fact, it is the ratio of which determines the angular anisotropy and the asymptotic powerlaw dependence of as well as the characteristic decay length of spinmisalignment fluctuations (Honecker & Michels, 2013).
In agreement with the assumption of a sharp interface in the nuclear (grain) microstructure (compare Fig. 1) both h^{2}(q) and vary asymptotically as q^{4}. Together with the micromagnetic response functions which, respectively, vary as and [compare equations (15) and (16), and see Fig. 3], this results in with n ranging between 6 and 8 (Honecker & Michels, 2013). We emphasize that other models for the anisotropyfield microstructure may result in different powerlaw exponents of ; in particular, the h^{2}(q) that are related to the longrange stress fields of dislocations are expected to give rise to asymptotic power laws that are different from the Porod exponent (Seeger, 1959; Heuser, 1994; Thomson et al., 1999; Maxelon et al., 2001; Long & Levine, 2005). This is, however, the subject of further investigations.
By inserting equations (42) and (43) into the θaveraged spinmisalignment SANS cross sections [equations (19) and (23)], we can express the onedimensional correlation functions of the spinmisalignment SANS [equation (36)] as
for and
for . Note that c(r) for the parallel geometry is [in contrast to c(r) for the perpendicular case] independent of both H_{p} and ; the dependence of c(r) on the applied magnetic field H_{i} and on the magnetic interactions (A, M_{s}) is contained in the function p(q, H_{i}) [compare equation (11)]. We also reemphasize that we have assumed that both Fourier coefficients h^{2}(q) and can be written as the product of the same form factor P(q) and S(q); this assumption might be relaxed, e.g. when studying diffusion zones or core–shelltype nanoparticle structures with reduced surface magnetization (Heinemann et al., 2000). The averaged twodimensional correlation function [equation (41)] is obtained by making the corresponding replacements in equations (45) and (46).
6. Results and discussion
The following materials parameters were used in the calculations: saturation magnetization , exchangestiffness constant A = 2.5 ×10^{11} J m^{1} and R = 5 nm for the particle radius in the sphere form factor P(q) [equation (44)].
6.1. Onedimensional correlation functions
All results in this section are obtained by numerical integration of equations (45) and (46), which are based on the onedimensional correlation function equation (36). In the first set of calculations, we concentrate on the dependence of the correlation functions on the applied magnetic field H_{i}, scattering geometry ( and ), ratio , singleparticle form factor P(q) and S(q).
Fig. 4 displays the results for c(r) at several values of H_{i} and for both scattering geometries, assuming a dilute scattering system [S(q) = 1] and . The dotted horizontal lines indicate the value of the correlation length l_{C} of the spin misalignment, which can be taken as a measure of the size of inhomogeneously magnetized regions around defects. l_{C} is defined as the exp(1) decay length, i.e. c(r = l_{C}) = exp(1). Note, however, that this definition does not imply that the correlations decay exponentially. In fact, it is readily verified that the spinmisalignment correlations that are investigated in this study do not decay exponentially. We would also like to mention that an alternative route to extracting a spinmisalignment length may be realized by the computation of moments of the correlation function; for instance, for exponentially decaying c(r) the above definition and are equivalent.
Increasing H_{i} results in both scattering geometries in the suppression of transverse spinmisalignment fluctuations and in a concomitant reduction of the c(r) and reduced l_{C} values. At small fields, l_{C} may take on values of the order of 100 nm, which decrease to values of the order of the assumed particle size, here R = 5 nm, for fields larger than a few tesla [see also dotted horizontal line in Fig. 5(b)]. For the chosen limiting case of , the difference between the c(r) and the l_{C}(H_{i}) in the two scattering geometries is only minor (see Fig. 5). However, noting that c(r) in the parallel geometry is independent of and with reference to Figs. 6 and 7, we note that this difference increases with decreasing value of .
Within the framework of our micromagnetic SANS theory of bulk et al., 2013; Metlov & Michels, 2015), the magnetic microstructure in real space, , corresponds to a complicated convolution product between the magnetic anisotropyfield microstructure and micromagnetic functions. As a result, smoothly varying magnetization profiles are at the origin of the related spinmisalignment scattering. In agreement with the absence of a sharp interface in the magnetic microstructure (compare Fig. 1b), we note that the correlation functions of bulk enter the origin r = 0 with zero slope (Bick et al., 2013), so that
(Honeckerfor (where k is a constant). This observation may be compared to the well known result for scattering, where (for isolated uniform particles) the first derivative of c(r) evaluated at r = 0 is related to the particle surface. In particular, for small r, the correlation function can be expanded as (Porod, 1982)
where the `differential' parameters a, b, c are related to the size and shape of the particle; for example, for a uniform sphere one finds a =  3/(4 R), b = 0 and c = 1/(16 R^{3}).
The effect of the ratio on the correlation functions and on the l_{C} values is shown in Figs. 6 and 7 [for and S(q) = 1]. Perturbations in the spin microstructure that are dominated by fluctuations of the magnetic anisotropy field () decay on a larger length scale than magnetostatically dominated () perturbations.
For soft magnets (with low crystalline anisotropy), the following relation for l_{C}(H_{i}) has previously been suggested (Michels, 2014):
Equation (49) provides an excellent description of the fielddependent correlations [solid lines in Figs. 5(b) and 7 with R = 5 nm, A = 2.5 ×10^{11} J m^{−1} and T]. At large fields, when the spinmisalignment SANS is small and the exchange length l_{H} takes on values of a few nanometres, l_{C} reflects, irrespective of , the size of the (in this case spherical) defect.
For the perpendicular scattering geometry, Fig. 8 displays (for T) the correlation function for different singleparticle form factors P(q), ignoring interparticle interactions [S(q) = 1]. In addition to the sphere form factor [equation (44)], we have used in the expressions for h^{2}(q) and the cylinder form factor (radius: R; length: L) (Pedersen, 1997),
and the form factor of an ellipsoid of revolution (semiaxes: R, R, ),
J_{1}(x) denotes the firstorder Bessel function, j_{1}(x) is the firstorder spherical Bessel function and . Note that equation (51) reduces to the sphere form factor for . Besides the cylinder and ellipsoid of revolution form factor we have also used other form factors (data not shown); the above form factors were chosen because they allow one to investigate different limiting cases (from thin circular discs to elongated spheroids and elongated thin rods). Examples for bulk magnetic materials with elongated cylindrically or elliptically shaped precipitates are alnico magnets (Zhou et al., 2014), which are nanostructured alloys composed of Fe, Al, Ni and Co.
It is seen in Fig. 8 that for a given form factor the shape of the correlation function and the value of the correlation length depend (as expected) on the particle dimensions. Isotropically distributed cylinders (dashed lines) with a radius equal to the radius of the ellipsoid of revolution and a length result in nearly the same (slightly larger) correlation functions as the ellipsoid of revolution. l_{C} at large fields appears to be related to the smallest dimension of the particle, although the precise dependency of on the particle dimensions is not clear to us. Note also that for the case of very thin discs and oblate spheroids () the correlation function still approaches the origin with zero slope (which becomes visible only for small r).
Finally, Fig. 9 illustrates the effect of interparticle interactions on the correlation function (Fig. 9a) and correlation length (Fig. 9b). In order to model the effect of dense packing, we have used the Percus–Yevick hardsphere for S(q) (Kinning & Thomas, 1984) in equations (45) and (46) and, as before, the sphere form factor for P(q). Note also that the hardsphere interaction radius R_{HS} in S(q) was set equal to the sphere radius R.
It is clearly seen that with increasing particle η the range of the correlations decreases. However, the characteristic features of the only become visible at relatively large values of η (above about 20%), while at the lower end of η values both c(r) and l_{C}(H_{i}) are smoothly decaying functions. Furthermore, we note that with increasing η, i.e. with increasing interparticle interactions, we progressively introduce, in addition to the original (diffuse) spinmisalignment length l_{C}, a second structural correlation length into the system (compare e.g. the hump in l_{C} at around 50 mT for ).
The field dependence of this feature is depicted in Fig. 10, where we show c(r) for several H_{i} and for ; here, we see that slight changes in H_{i} result in relatively large jumps in l_{C} [l_{C}(0.08 T) = 7.7 nm and l_{C}(0.05 T) = 13.8 nm]. This is an artifact which is clearly related to the strong structural correlations, and the determined correlation length now represents a fielddependent (unknown) average over the structural and the magnetic spinmisalignment correlation lengths. We note that by using other definitions for l_{C}, for instance in terms of some integral weight over c(r), the position of the artifact on the H_{i} axis may be different but the effect of S(q) will still become visible.
6.2. Twodimensional correlation functions
Since the spinmisalignment SANS θ dependence of is a consequence of the trigonometric functions in the (which are due to the dipolar neutron–magnetic interaction) and of the θ dependence of the magnetization Fourier coefficients (which is due to the internal magnetostatic interaction) (Michels, 2014). Figs. 11(a)–11(d) show [equation (6)] at selected applied magnetic fields [and for S(q) = 1]. The change in the angular anisotropy that becomes visible in Figs. 11(a)–11(d), from a spiketype anisotropy at low fields (a) to a cloverleafshaped anisotropy at large fields (d), is related to the field dependence of the Fourier coefficients and demonstrates that different terms in the response functions [equations (9) and (10)] dominate in different field regimes. For instance, the spike anisotropy (Fig. 11a) was recently observed in an isotropic sintered Nd–Fe–B magnet (Périgo et al., 2014); it is related to magnetostatic terms in the denominator of the response functions.
is highly anisotropic for , the corresponding correlation function may also be anisotropic. We reemphasize that the angularThe corresponding twodimensional correlation functions, computed according to equation (37), are displayed in Figs. 11(e)–11(h), where we plot the c(y, z) at the same fields as the . While the spinmisalignment SANS at small fields (Figs. 11a and 11b) is enhanced parallel to the appliedfield direction, the correlation function exhibits maxima in the direction perpendicular to the field; the range of the correlations extends to several hundreds of nanometres (Figs. 11e and 11f). Increasing the field results in the suppression of the correlations. At the largest field possesses a nearly fourfold anisotropy with maxima along the detector diagonals and minima along the horizontal and vertical axes (Fig. 11d), which translate into the corresponding extrema in c(y, z) (Fig. 11h).
In Fig. 12(a), we depict the correlation function along different directions: while the correlation length at 1.2 T varies only relatively little with direction (from 8.8 to 10.9 nm), the functional dependencies of the c(y, z) are significantly different, with the correlation function along the horizontal z direction becoming negative at nm; the curves in Fig. 12(a) were obtained by solving equation (38) for , , . In nuclear SANS, negative values of the distance distribution function p(r) are attributed to distances that connect regions with opposite sign of the scattering length density more frequently than regions with the same sign (Glatter & Kratky, 1982). However, for magnetic SANS, such an easily accessible interpretation of the correlation function c(r) of the spinmisalignment SANS in terms of a specific magnetization distribution is not straightforward; this is mainly related to the (above mentioned) fact that c(r) does not directly represent the correlations in the magnetic microstructure (as does C_{SM}), but also includes the magnetodipolar interaction of the neutrons with the sample (via the trigonometric functions and the cross term in the cross section). The anisotropy of the correlations is further depicted in Fig. 12(b), where we show a contour plot for several values of H_{i}. This graph reveals a relatively weak anisotropy of l_{C}. At small fields, the correlations along the vertical (y) direction decay on a larger length scale than along the horizontal (z) direction; with increasing field, the anisotropy becomes less pronounced.
Fig. 13 compares (for ) the results for the onedimensional [equation (36)] and the averaged twodimensional [equation (41)] correlation functions of the spinmisalignment SANS. We recall that the former is obtained by threedimensional integration of the azimuthally averaged , and the latter by twodimensional integration of (compare §4). At small fields, the results for c(r) and l_{C}(H_{i}) differ considerably, whereas for T both equations yield almost the same correlation lengths.
The question may arise as to which one of these correlation functions should be used in order to analyze experimental data. From an experimental point of view, the averaged twodimensional equation (41) reflects the dataanalysis procedure, namely that the measured is a function of only two independent components of the scattering vector; in fact, in the smallangle approximation only probes correlations in the directions perpendicular to the incident beam. Reconstruction (from experimental ) of the onedimensional c(r) (which is an orientation average of the threedimensional correlation function) is always an extrapolation.
6.3. Comparison with experimental data
In order to test our magnetic SANS theory, we depict in Fig. 14 a comparison between experiment and theory; in particular, we have fitted equations (36) and (41) [using, in each case, equation (19) for ] to experimental data for the correlation function of the spinmisalignment SANS of nanocrystalline Co and Ni (Michels et al., 2003). These C(r) data have previously been analyzed by Michels & Bick (2013) using a simple approach based on the autocorrelation function of the spin misalignment, neglecting terms due to spatial fluctuations of the saturation magnetization. Such contributions are included in the present theory via the term S_{M} R_{M} in equation (19). The nanocrystalline Co and Ni samples constitute fully dense polycrystalline metals with average crystallite sizes of D = 10 nm (Co) and D = 49 nm (Ni) (Weissmüller et al., 2001). The experimental SANS data of both samples were recorded between nm^{−1} and nm^{−1}. The correlation functions were then obtained by direct Fourier transformation according to equation (36), so that this expression should actually also be used for the data analysis. Nevertheless, we have also employed the twodimensional equation (41) for fitting the experimental C(r) data, which is motivated by the fact that for larger applied fields the difference between the two correlation functions is only minor (compare Fig. 13). In the following discussion, one should therefore keep in mind that for the analysis of this particular C(r) data set equation (36) represents the proper theoretical model.
In the fitting procedure, the integrals in equations (36) and (41) were approximated by discrete sums, where the upper integration limit of `' was taken as –10 nm^{−1} and the typical q resolution was set to –0.02 nm^{−1}. The resulting expressions were fitted by means of a nonlinear (Levenberg–Marquardt) fitting routine to the experimental C(r) data. We have treated the exchangestiffness constant A, the ratio and R as global fit parameters. Since we work with unnormalized C(r) data, we have introduced fielddependent local scaling constants K_{1}, K_{2}, K_{3} and K_{4} (one for each data set); M_{s} = 1434 kA m^{−1} for Co and M_{s} = 522 kA m^{−1} for Ni were held constant. Since the experimental SANS data [e.g. Fig. 1 of Michels et al. (2003)] do not give a visible indication of a strong impact of dense packing, we have for simplicity decided to set S = 1. The results for the global fit parameters are summarized in Table 1. The data analysis was restricted to r values below about 50 nm and to fields larger than 50 mT, where the magnetization of both samples approaches saturation (Weissmüller et al., 2001).

As is seen in Fig. 14 (solid and dashed lines), both equations provide a reasonable global description of the fielddependent correlations. The obtained values for the anisotropyfield radii R of both materials are in the range 8–13 nm, slightly smaller than the ones estimated previously (Michels & Bick, 2013). The parameter R characterizes the length scale over which the magnetic anisotropy field is uniform; for singlecrystal grains, R is sensibly related to the average crystallite size (compare Fig. 1a). Therefore, the finding –10 nm for Co suggests that the magnetic anisotropy field is approximately homogeneous on a length scale of the order of the average grain size of 10 nm, whereas for Ni nonuniformities in exist on a scale smaller than the average crystallite size of 49 nm, presumably related to twin faults or to the defect cores of grain boundaries (Michels et al., 2003). While the obtained values for the exchangestiffness constant of Ni (using both equations) are larger by a factor of about two than the ones reported in the literature (Kronmüller & Fähnle, 2003), the A value for Co using equation (41) agrees excellently with literature data and with the result of our previous SANS data analysis (in Fourier space) (Michels et al., 2003). Values for the ratio of have not been determined previously for these materials, but our results suggest [except for the case of Ni using equation (41)] that perturbations in the spin microstructure due to spatially fluctuating magnetic anisotropy fields dominate over magnetostatic fluctuations. This might be expected, since in singlephase variations in M_{s} are relatively small, compared to e.g. nanocomposites (Michels et al., 2006). Overall, the good agreement between experiment and theory suggests that equation (36) may be used for the analysis of realspace correlations of bulk magnetic materials; equation (41) may also be employed for the analysis of experimental data, provided that the original has been Fourier transformed according to equation (41).
7. Summary and conclusion
On the basis of a recent micromagnetic theory for the magnetic SANS c(r) of the spinmisalignment SANS depends on the applied magnetic field and, for , on the ratio of magnetic anisotropy field strength H_{p} to magnetization jump at internal interfaces. Additional in c(r) relate to the particle (anisotropyfield) form factor or to the inclusion of interparticle correlations via a The result for c(r) (for ) [equation (45)] demonstrates a strong impact of on the shape and range of the correlations: magnetostatically dominated correlations () decay on a rather short length scale, whereas anisotropyfielddominated correlations () are characterized by a longrange decay, which is reasonably described by equation (49). The difference between the correlation functions in the two scattering geometries ( and ) increases with decreasing ratio of . The correlation functions do not decay exponentially and approach the origin with zero slope; as far as equation (36) is concerned, this is consistent with the absence of a sharp interface in the magnetic microstructure. Experimental data for the correlation function of the spinmisalignment SANS of nanocrystalline Co and Ni have been successfully analyzed using the here presented theoretical expressions. It would also be of interest to employ the present approach for studying longrange magnetic correlations, as accessible on a USANS instrument (Jericha et al., 2013), or the magnetic microstructure of stateoftheart nanocrystalline NdFeBbased permanent magnets (Bick et al., 2013; Yano et al., 2014; Périgo et al., 2015; Saito et al., 2015).
of inhomogeneous bulk we have studied the corresponding magnetic fielddependent spinmisalignment correlations in real space. The correlation functionAPPENDIX A
Autocorrelation function of the spin misalignment
In the highfield limit and for a general orientation of the wavevector , the solution, in Fourier space, of the linearized balanceoftorques equation (1) can be written as (Honecker & Michels, 2013)
For q_{x} = 0 or q_{z} = 0, one obtains the expressions for and , respectively, which enter the equations for the perpendicular or the parallel SANS cross sections [equations (3) and (20)]. Averaging the expression for over the orientation (angle β) of the magnetic anisotropyfield Fourier coefficient,
results in
where and p = p(q, H_{i}) is given by equation (11). Inserting this function into the normalized version of equation (33) allows us to obtain the autocorrelation function of the spin misalignment, c_{SM}(r), by threedimensional integration. Fig. 13 displays c_{SM}(r) and l_{C}(H_{i}) and compares the results with the correlation functions of the spinmisalignment SANS equations (36) and (41).
Acknowledgements
We thank the National Research Fund of Luxembourg for financial support (CORE/INTER project No. INTER/DFG/12/07 and ATTRACT project No. FNR/A09/01). Critical reading of the manuscript by Élio Périgo, Dirk Honecker, Sergey Erokhin and Konstantin Metlov is gratefully acknowledged.
References
Alinger, M. J., Odette, G. R. & Hoelzer, D. T. (2009). Acta Mater. 57, 392–406. Web of Science CrossRef CAS Google Scholar
Bergner, F., Pareige, C., Kuksenko, V., Malerba, L., Pareige, P., Ulbricht, A. & Wagner, A. (2013). J. Nucl. Mater. 442, 463–469. Web of Science CrossRef CAS Google Scholar
Bick, J.P., Honecker, D., Döbrich, F., Suzuki, K., Gilbert, E. P., Frielinghaus, H., Kohlbrecher, J., Gavilano, J., Forgan, E. M., Schweins, R., Lindner, P., Birringer, R. & Michels, A. (2013). Appl. Phys. Lett. 102, 022415. Web of Science CrossRef Google Scholar
Bischof, M., Staron, P., Michels, A., Granitzer, P., Rumpf, K., Leitner, H., Scheu, C. & Clemens, H. (2007). Acta Mater. 55, 2637–2646. Web of Science CrossRef CAS Google Scholar
Brown, W. F. Jr (1963). Micromagnetics. New York: Interscience Publishers. Google Scholar
Coppola, R., Kampmann, R., Magnani, M. & Staron, P. (1998). Acta Mater. 46, 5447–5456. Web of Science CrossRef CAS Google Scholar
Döbrich, F., Kohlbrecher, J., Sharp, M., Eckerlebe, H., Birringer, R. & Michels, A. (2012). Phys. Rev. B, 85, 094411. Google Scholar
Erokhin, S., Berkov, D., Gorn, N. & Michels, A. (2012). Phys. Rev. B, 85, 024410. Web of Science CrossRef Google Scholar
Feigin, L. A. & Svergun, D. I. (1987). Structure Analysis by SmallAngle Xray and Neutron Scattering. New York: Plenum Press. Google Scholar
Frandsen, B., Yang, X. & Billinge, S. J. L. (2014). Acta Cryst. A70, 3–11. Web of Science CrossRef IUCr Journals Google Scholar
Fritz, G. & Glatter, O. (2006). J. Phys. Condens. Matter, 18, S2403. Web of Science CrossRef Google Scholar
FritzPopovski, G. (2013). J. Appl. Cryst. 46, 1447–1454. Web of Science CrossRef CAS IUCr Journals Google Scholar
FritzPopovski, G. (2015). J. Appl. Cryst. 48, 44–51. Web of Science CrossRef CAS IUCr Journals Google Scholar
Glatter, O. (1977). J. Appl. Cryst. 10, 415–421. CrossRef IUCr Journals Web of Science Google Scholar
Glatter, O. & Kratky, O. (1982). Small Angle Xray Scattering. London: Academic Press. Google Scholar
Guinier, A. & Fournet, G. (1955). SmallAngle Scattering of Xrays. New York: Wiley. Google Scholar
Hansen, S. (2000). J. Appl. Cryst. 33, 1415–1421. Web of Science CrossRef CAS IUCr Journals Google Scholar
Hansen, S. (2012). Estimation of Distribution Functions from SmallAngle Scattering Data, https://www.bayesapp.org/. Google Scholar
Heinemann, A., Hermann, H., Wiedenmann, A., Mattern, N. & Wetzig, K. (2000). J. Appl. Cryst. 33, 1386–1392. Web of Science CrossRef CAS IUCr Journals Google Scholar
Heuser, B. J. (1994). J. Appl. Cryst. 27, 1020–1029. CrossRef CAS Web of Science IUCr Journals Google Scholar
Honecker, D., Dewhurst, C. D., Suzuki, K., Erokhin, S. & Michels, A. (2013). Phys. Rev. B, 88, 094428. Web of Science CrossRef Google Scholar
Honecker, D., Ferdinand, A., Döbrich, F., Dewhurst, C. D., Wiedenmann, A., GómezPolo, C., Suzuki, K. & Michels, A. (2010). Eur. Phys. J. B, 76, 209–213. Web of Science CrossRef CAS Google Scholar
Honecker, D. & Michels, A. (2013). Phys. Rev. B, 87, 224426. Web of Science CrossRef Google Scholar
Jericha, E., Badurek, G. & Gösselsberger, C. (2013). Phys. Proc. 42, 58–65. CrossRef CAS Google Scholar
Kinning, D. J. & Thomas, E. L. (1984). Macromolecules, 17, 1712–1718. CrossRef CAS Web of Science Google Scholar
Kronmüller, H. & Fähnle, M. (2003). Micromagnetism and the Microstructure of Ferromagnetic Solids. Cambridge University Press. Google Scholar
Lembke, U., Hoell, A., Kranold, R., Müller, R., Schüppel, W., Goerigk, G., Gilles, R. & Wiedenmann, A. (1999). J. Appl. Phys. 85, 2279–2286. Web of Science CrossRef CAS Google Scholar
Löffler, J. F., Braun, H. B., Wagner, W., Kostorz, G. & Wiedenmann, A. (2005). Phys. Rev. B, 71, 134410. Google Scholar
Long, G. G. & Levine, L. E. (2005). Acta Cryst. A61, 557–567. Web of Science CrossRef CAS IUCr Journals Google Scholar
Maxelon, M., Pundt, A., PyckhoutHintzen, W., Barker, J. & Kirchheim, R. (2001). Acta Mater. 49, 2625–2634. Web of Science CrossRef CAS Google Scholar
Metlov, K. L. & Michels, A. (2015). Phys. Rev. B, 91, 054404. Web of Science CrossRef Google Scholar
Michaud, P., Delagnes, D., Lamesle, P., Mathon, M. H. & Levaillant, C. (2007). Acta Mater. 55, 4877–4889. Web of Science CrossRef CAS Google Scholar
Michels, A. (2010). Phys. Rev. B, 82, 024433. Web of Science CrossRef Google Scholar
Michels, A. (2014). J. Phys. Condens. Matter, 26, 383201. Web of Science CrossRef PubMed Google Scholar
Michels, A. & Bick, J.P. (2013). J. Appl. Cryst. 46, 788–790. Web of Science CrossRef CAS IUCr Journals Google Scholar
Michels, A., Döbrich, F., Elmas, M., Ferdinand, A., Markmann, J., Sharp, M., Eckerlebe, H., Kohlbrecher, J. & Birringer, R. (2008). EPL, 81, 66003. Google Scholar
Michels, A., Elmas, M., Döbrich, F., Ames, M., Markmann, J., Sharp, M., Eckerlebe, H., Kohlbrecher, J. & Birringer, R. (2009). EPL, 85, 47003. Google Scholar
Michels, A., Erokhin, S., Berkov, D. & Gorn, N. (2014). J. Magn. Magn. Mater. 350, 55–68. Web of Science CrossRef CAS Google Scholar
Michels, A., Vecchini, C., Moze, O., Suzuki, K., Pranzas, P. K., Kohlbrecher, J. & Weissmüller, J. (2006). Phys. Rev. B, 74, 134407. Web of Science CrossRef Google Scholar
Michels, A., Viswanath, R. N., Barker, J. G., Birringer, R. & Weissmüller, J. (2003). Phys. Rev. Lett. 91, 267204. Web of Science CrossRef PubMed Google Scholar
Mortensen, K. & Pedersen, J. S. (1993). Macromolecules, 26, 805–812. CrossRef CAS Web of Science Google Scholar
Ohnuma, M., Hono, K., Linderoth, S., Pedersen, J. S., Yoshizawa, Y. & Onodera, H. (2000). Acta Mater. 48, 4783–4790. Web of Science CrossRef CAS Google Scholar
Pedersen, J. S. (1997). Adv. Colloid Interface Sci. 70, 171–210. CrossRef CAS Web of Science Google Scholar
Périgo, E. A., Gilbert, E. P., Metlov, K. L. & Michels, A. (2014). New. J. Phys. 16, 123031. Google Scholar
Périgo, E. A., Gilbert, E. P. & Michels, A. (2015). Acta Mater. 87, 142–149. Google Scholar
Porod, G. (1982). Small Angle Xray Scattering, edited by O. Glatter & O. Kratky, pp. 17–51. London: Academic Press. Google Scholar
Saito, K., Ueno, T., Yano, M., Harada, M., Shoji, T., Sakuma, N., Manabe, A., Kato, A., Keiderling, U. & Ono, K. (2015). J. Appl. Phys. 117, 17B302. Web of Science CrossRef Google Scholar
Šaroun, J. (2000). J. Appl. Cryst. 33, 824–828. Web of Science CrossRef IUCr Journals Google Scholar
Seeger, A. K. (1959). J. Appl. Phys. 30, 629–637. CrossRef CAS Web of Science Google Scholar
Svergun, D. I. & Koch, M. H. J. (2003). Rep. Prog. Phys. 66, 1735–1782. Web of Science CrossRef CAS Google Scholar
Thomson, R., Levine, L. E. & Long, G. G. (1999). Acta Cryst. A55, 433–447. Web of Science CrossRef CAS IUCr Journals Google Scholar
Watson, G. N. (1966). A Treatise on the Theory of Bessel Functions, 2nd ed. Cambridge University Press. Google Scholar
Weissmüller, J., McMichael, R. D., Michels, A. & Shull, R. D. (1999). J. Res. Natl Inst. Stand. Technol. 104, 261–275. Google Scholar
Weissmüller, J., Michels, A., Barker, J. G., Wiedenmann, A., Erb, U. & Shull, R. D. (2001). Phys. Rev. B, 63, 214414. Google Scholar
Weissmüller, J., Michels, A., Michels, D., Wiedenmann, A., Krill, C. E. III, Sauer, H. M. & Birringer, R. (2004). Phys. Rev. B, 69, 054402. Google Scholar
Yano, M., Ono, K., Harada, M., Manabe, A., Shoji, T., Kato, A. & Kohlbrecher, J. (2014). J. Appl. Phys. 115, 17A730. Web of Science CrossRef Google Scholar
Zhou, L., Miller, M., Lu, P., Ke, L., Skomski, R., Dillon, H., Xing, Q., Palasyuk, A., McCartney, M., Smith, D., Constantinides, S., McCallum, R., Anderson, I., Antropov, V. & Kramer, M. (2014). Acta Mater. 74, 224–233. Web of Science CrossRef CAS Google Scholar
This is an openaccess article distributed under the terms of the Creative Commons Attribution (CCBY) Licence, which permits unrestricted use, distribution, and reproduction in any medium, provided the original authors and source are cited.