research papers
On the use of grazing-incidence small-angle X-ray scattering (GISAXS) in the morphological study of ion-implanted materials
aINFM, OGG, c/o ESRF, GILDA CRG, BP 220, F-38043 Grenoble, France, bUniversità di Venezia, Dipartimento di Chimica Fisica, Dorsoduro 2137, I-30123 Venezia, Italy, cUniversità di Padova, Dipartimento di Fisica, Via Marzolo 8, I-35131 Padova, Italy, dINFM, OGG, c/o ILL, BP 156, F-38042 Grenoble, France, eICTPN, CNR, Via Ugo La Malfa 153, I-90146 Palermo, Italy, and fUniversità di Palermo, Dipartimento di Chimica Inorganica, Viale delle Scienze, I-90128 Palermo, Italy
*Correspondence e-mail: dacapito@esrf.fr
Grazing-incidence small-angle X-ray scattering has become a widely used technique for the morphological analysis of surface systems. Here it is show how this technique can be applied to a buried system, like metallic clusters in glass obtained by ion implantation. The optimization of the data-collection geometry is described as well as the details of the quantitative data analysis. An experimental example on Cu + Au-implanted glasses shows the potentiality of the technique.
Keywords: grazing-incidence small-angle X-ray scattering; ion implantation; metallic nanoparticles; glass–nanocluster composites; small-angle neutron scattering.
1. Introduction
Composite materials consisting of nanometer-scale metallic particles embedded in glass have recently received great attention for their peculiar physical properties (Gonella & Mazzoldi, 2000). In particular, the reduced particle size creates a number of unusual effects, like quantum confinement of the electronic wavefunctions (Alivisatos, 1996; Halperlin, 1986) or dielectric confinement (Haglund, 1995), that strongly modify the behavior of such systems with respect to the corresponding bulk forms. Applications for these materials have been proposed in optoelectronics for fast switching (Gonella & Mazzoldi, 2000) or in magnetic recording (Haglund, 1995), to cite only a few. A key point in the construction of a nanocluster-based device is the control of the morphology in terms of cluster shape, dimension, size spread and cluster arrangement in the matrix. Laboratory techniques, like (TEM), can provide only some of the cited parameters and in an extremely spatially reduced portion of the sample. Grazing-incidence small-angle X-ray scattering (GISAXS) is a recently developed technique (Levine et al., 1989; Naudon & Thiaudiere, 1997) for the analysis of surface systems. Detailed descriptions of systems of particles deposited on a surface have been carried out (Robach et al., 1999; Babonneau et al., 1999, 2001; Naudon et al., 2000; Stemmer et al., 2003; Renaud et al., 2003) and a complete treatment of the GISAXS theory and analysis methods for this kind of material has been recently presented (Lazzari, 2002).
GISAXS has also found a particularly useful application in the study of ion-implanted materials, with particular attention to metal implants in glasses (d'Acapito et al., 1998; Cattaruzza et al., 2000; Maurizio, Longo et al., 2003), semiconductors in glasses and crystalline silicon (Dubcek et al., 2003; Desnica et al., 2003), and rare gases in silicon (Dubcek et al., 2004). It must be noted that, differently from surface-deposited particles, in the case of implanted materials the layer under analysis is typically a few hundreds of Å below the surface. This leads to the necessity of a non-trivial choice of the experimental conditions, in order to maximize the signal from the implanted species without being overwhelmed by that from the substrate. In this contribution we will focus on the use of GISAXS on ion-implanted materials, describing (i) the procedure to optimize the data collection, and (ii) the quantitative analysis, that takes into account the particular physical processes driving the formation of clusters in implanted substrates.
2. Cluster formation by ion implantation
An efficient way of obtaining nanoparticles in matrices is the ion-implantation technique (Cattaruzza, 2000). In this way, problems related to low atomic solubility or diffusion coefficients are easily overcome, practically permitting the doping of any matrix with any atomic species. In an ion-implantation process, ions with energies in the keV range enter the matrix and start to lose energy. In the initial part of the path, the ions transfer their to the matrix mainly via electronic processes. Then, when the ion has reached a sufficiently low speed, nuclear interactions dominate. Ions stop in the substrate following a distribution function, that can be parametrized by two quantities: the projected range Rp and the straggling Rs. Both depend on the implantation energy E and on the of the implanted species Z, in a way that higher E leads to higher Rp and Rs values. An opposite effect on these parameters is played by Z. By using simulation programs based on Monte Carlo procedures it is possible to determine these values ab initio (Biersack & Haggmark, 1980). A typical distribution of implants in a matrix is shown in Fig. 1 and was obtained using TRIM (Biersack & Haggmark, 1980).
Another fundamental parameter in implantation processes is the fluence that gives the number of ions per unit surface deposited in the material. As Rp and Rs are independent of (at least for low fluences), it is also possible to calculate the volume concentration of the implant given a suitable set of initial conditions. The formation of clusters in the matrix can be seen, in a simplified way, as a diffusion process followed by a chemical interaction. The implanted atomic species, after deposition, start to migrate in the matrix at a rate determined by their diffusion constant Dc. When two diffusing atoms come in contact and the metal Gibbs free energy is lower than that of the related oxides (if implanting an oxide matrix), they start to form a nucleation center for the future metallic cluster. This process, which already exists during the implantation, can be continued with a subsequent annealing of the material. The first step, of kinematical nature, can be described by stating that the average distance l between two atoms (that depend on their volume density ) must be smaller than the diffusion length l,
l is determined by the Dc and by the treatment time t, as stated in (1).
The second condition means that the compound we wish to form in the matrix (metallic cluster, namely) is thermodynamically more stable than other possible (oxides) compounds, i.e.
This picture represents an idealized model for the matrix–implant interaction and accounts for the simplest processes. The following more complex phenomena are not accounted for:
(i) The matrix damage owing to the implant can locally alter the distribution of chemical species leading to a modified behavior of the implant from what is foreseen for the unperturbed matrix. An example of this effect is described by d'Acapito et al. (2000).
(ii) The surface atoms are usually a considerable fraction of the total atoms in the cluster, so phenomena like formation of core–shell clusters or de-alloying of mixed component particles are likely to be formed during particular treatments (Mattei et al. 2003).
3. The GISAXS technique applied to implanted materials
The method is based on the observation that X-rays have a refraction index n of less than 1 (Parratt, 1954),
where and are the so-called n is smaller than 1 can be seen, in a semiclassical way, as originating in the lag in the response of the bound electrons to the electric field that oscillates at a frequency higher than their typical resonance value (Feynman et al., 1970). The common laws for refraction and reflection still hold, but the material plays the role of the low n zone and the vacuum that of the high n zone. In particular, the total external reflection of a beam impinging with an angle lower than a critical value is possible. In this condition the components of the wavevector k of the refracted beam parallel to the surface (kx, ky) are real, whereas the component perpendicular to the surface (kz) is imaginary. Thus, the beam propagates parallel to the surface and penetrates in the material with an exponentially damped wave described by a damping parameter z1/e. By increasing the incidence angle, kz starts to develop a non-zero real part and propagation sets up also inside the material (Parratt, 1954).
corrections. The fact that the real part ofAs already pointed out, when dealing with implanted materials the interesting part is only a thin layer embedded in a much thicker (about four orders of magnitude) substrate. The problem is then to confine the probe beam into a zone slightly below the surface, where the atomic density of the implanted species is maximum, minimizing the contribution of the substrate. Indeed, when working in normal incidence the substrate signal overwhelms the signal from the implanted layer, whereas, when working in total ) in the case of X-ray diffraction. For this reason, when dealing with ion-implanted materials, the incidence angle has to be optimized to correctly carry out this kind of experiment. This is achieved by analyzing a typical reflectivity curve, shown in Fig. 2, relative to a heavy ion implant in a light matrix. In the low-angle region the reflectivity is maximum as this is the total reflection regime for the glass/air interface. Above 0.2° the signal starts to decrease, indicating that the beam is penetrating the glass. However, it is still reflected at the glass/implanted-region interface, as evidenced by the presence of Kiessig fringes (Kiessig, 1931) originated by the interference between the beam partially reflected at the air/glass interface and that reflected at the glass/implanted-region interface. The observation of the fringes in this kind of sample would seem surprising, since the glass/implanted-region interface is not sharp but has a roughly Gaussian shape, as shown in Fig. 1. However, it must be noted that, owing to the low scattering angle, the value of the exchanged wavevector q is also low, of the order of ∼0.03 Å-1 in the present case, corresponding to a spatial resolution of ∼210 Å. This value is well comparable with the straggling of the implant profile, so the probe beam sees the system as an almost sharp interface. Clearly, owing to the diffuseness of the interface, the intensity of the Kiessig fringes is much lower than that observed in abrupt interfaces. The presence of the maxima indicates that, in the related angular zone, the beam is confined to the implanted region and can be used for structural studies with enhanced spatial sensitivity to the region of interest. Confirmation of this can be found in Fig. 3, which represents the intensity (in the substrate) of the refracted electric field in a two-dimensional incidence-angle-versus-depth space. Here, a high field intensity zone can be seen at a depth corresponding to the implanted region when the incidence angle is at a value near to the first fringe. Thus, by choosing the incidence angle on this point, a maximum sensitivity to the implanted region is obtained permitting experiments of X-ray diffraction or scattering. d'Acapito & Zontone (1999) report the effect of the choice of the incidence angle on the diffraction spectrum from an Ag-implanted sample, showing the maximization of the signal when working on the first fringe.
the limited penetration prevents the probe beam from joining the implanted region. Examples of these phenomena have already been shown by d'Acapito & Zontone (19994. Data collection and analysis
A typical experimental setup for the collection of GISAXS data consists of a precise sample mounting stage that permits the alignment of the sample in total reflection condition and a detector for specular reflectivity. The sample has to be aligned in a suitable way, as shown in the previous section, in order to maximize the sensitivity to the implanted region. After the collection of a reflectivity curve, and the choice of the correct incidence angle, the small-angle scattering spectrum can be collected. Both linear or two-dimensional detectors are suitable for this purpose even if, using the latter kind, a more efficient use of beam time is obtained, since the entire scattering pattern is recorded in a single shot. Typical integration times on an undulator beamline from a third-generation source (here, the European Synchrotron Radiation Facility, ESRF) are of the order of 10 min. Typical GISAXS spectra from ion-implanted glasses are shown in Figs. 4 and 5.
This kind of data permits the extraction of detailed information on the morphology of the particles, as already pointed out in the case of surface-deposited particles (Robach et al., 1999) or oriented ellipsoids in glass (Polizzi et al., 1998). When dealing with ion-implanted systems, different spatial distributions of the particles in the substrate can be observed. For example, the ring visible in the scattering pattern of Fig. 4 illustrates the presence of a well defined interparticle spatial correlation length; the radial symmetry indicates a uniform distribution of particles in space. When the particles are arranged in a layer, only in-plane correlations are visible and the scattering pattern takes the form shown in Fig. 5. In most cases, however, ion implantation in glassy systems leads to the formation of spherical particles distributed in a three-dimensional way. A consequence of this is that the obtained scattering patterns have a radial symmetry, and this can be exploited in the data analysis. The first step consists of reducing the two-dimensional data to a conventional one-dimensional scattering curve. For this reason, data have to be integrated at constant exchanged wavevector q over the various spatial orientations. The origin of the q scale can be determined by finding the center of the particle correlation ring using routines available on image-analysis programs (Hammersley et al., 1996), or by trials by angularly integrating two sectors opposite to the symmetry axis to obtain two curves with coinciding maxima. At this point we have to consider how refraction processes enter into the data-collection procedure. First of all, the direction of the probe beam (refracted in the substrate) is slightly different from the direction of the direct beam. Moreover, the outgoing (scattered) beams undergo a refraction at the glass–air interface, mainly at low scattering angles. The result of these two facts is that the q scale, in the so-calculated monodimensional scattering curve, needs to be suitably corrected in order to obtain reliable quantitative results. The problem of GISAXS data corrections related to refraction processes has been presented in the literature and procedures have been published by Kutsch et al. (1997) and Martorana et al. (2001).
A typical example of a monodimensional scattering curve is shown in Fig. 6. The ensemble of particles, created in a substrate by ion implantation, has peculiar features that must be kept in mind for a correct quantitative analysis:
(i) The particles have a spherical shape and grow by a diffusion process. They reasonably form an ensemble of hard non-penetrating spheres.
(ii) The average particle–particle correlation distance is of the same order of magnitude as the particle size. Both the particle S and form factor contribute in the same q region.
(iii) The particle distribution size is quite wide, its width being a considerable fraction of the mean value. From TEM images a correlation between particle size and depth in the substrate is evident, its details depending on the preparation parameters (Fig. 7).
A suitable model for data analysis has been developed, taking into account the previous points. The model proposed by our group, see the analysis of Cu + Ni-implanted samples by Cattaruzza et al. (2000), describes the scattered intensity I(q) as follows,
where is the particle form factor, is the particle ensemble N(R) is the particle size distribution. For each size R, we calculate a scattering curve I(R) as the product between and S and then we integrate over the size distribution. This corresponds to the so-called `local monodisperse approximation' (LMA) (Pedersen, 1993, 1994), which provides a better agreement with experimental data than the `decoupled approximation' (DA) (Cattaruzza et al., 2000). The observation recalled in (iii) justifies the use of this model. Looking in more detail, we use the sphere form factor (Guinier & Fournet, 1955)
andwhereas the particle ensemble is described by the well known Percus–Yevick ; Ashcroft & Lekner, 1966; Pedersen, 1993; Maurizio, Longo et al., 2003). Here, is the hard sphere (HS) packing fraction (particles per unit volume) and RHS is their radius. It should be noted that RHS is different from the particle radius R, and, in particular, is greater. This is because, as the particles grow with a diffusion process, i.e. the metal ions leave the glass matrix to migrate in the particle, a glass corona strongly depleted in metal ions tends to form around bigger particles. So, in an ideal fully aggregated system, all the particles are separated by spaces of pure substrate. The ratio between the particle radius and the hard-sphere radii (particle + depletion corona) can be calculated by supposing that all the ions originally contained in the hard-sphere volume (determined by the local density achieved at the end of the implant process) collapse in the central cluster. This cluster exhibits the typical metal density and from this we derive
(Percus & Yevick, 1958This ratio can be determined a priori by the density estimated by simulations of the implant and varied only slightly during the fitting process. Eventually the scattering curves calculated in this way are averaged over the particle size distribution, which in our case is described by a Weibull function (Pedersen, 1993) as shown here,
where Rm is the average value of the particle size and b is a parameter that controls the distribution width and shape (Pedersen, 1993). All these parameters have different effects on the calculated scattering curve. Different values of Rm affect the slope of the ln[I(q)] versus q2 curve at low exchanged momentum data above the correlation peak (bigger values of Rm lead to steeper curves), whereas the value of affects the position of the correlation peak (for higher packing fractions the peak shifts at higher q values). The b parameter affects the width of the particle size distribution function: for higher values a narrower distribution is obtained around Rm, and shape oscillations start to show up in the high-q part of the scattering curve. Correlations among the parameters are easily checked and accounted for in the error estimation, by using suitable codes [like MINUIT (1994)] for the curve fitting. In the data-fitting process, a further scale factor S and a constant background B are added to the modelling function. In principle it would be possible to use absolute intensities as is usual in measurements collected in transmission mode. However, a number of poorly controlled parameters contribute to prevent a reliable determination of absolute intensities when working at the critical angle:
(i) The scattered intensity depends on the intensity of the probe beam. This, in turn, strongly depends on the incidence angle at values near the critical value and on the details of the electronic density of the material. Both quantities are known only at an approximate level in most cases.
(ii) The surface roughness can contribute to modify the intensity of the refracted beam as well as to add a background to the scattering pattern.
For this reason we preferred to perform the data analysis without considering absolute intensities. It is worth noticing, at this point, that no scattering et al., 2003).
of the particle system and glass are cited in the whole treatment, making it suitable to both neutron and X-ray scattering. The experiment in the former case is strongly dependent on the achievement of a sufficient contrast between the particles and the matrix, a condition that can be actually satisfied in some cases as we will show in the next section. The model discussed can be used to process GISAXS (and SAXS) data to obtain some structural parameters from the investigated samples, such as the cluster size distribution and the intercluster mean distance. In principle this technique presents some advantages with respect to common TEM analysis, because it is not destructive, it averages over a macroscopic number of particles and permits a direct determination of the interparticle correlation distance. Moreover, time-resolved experiments (morphological evolution under thermal annealing or similar) can be carried out considering the relatively short acquisition time (Renaud5. Experimental examples
Here we show an example of how GISAXS can be used in the investigation of ion-implanted glasses. Silica substrates were implanted with Cu and Au ions at energies of 90 keV and 190 keV, respectively, and at different fluences . These implantation parameters determine a good overlap of the two implantation profiles in the sample. Three samples are described: the first (Au3Cu3) implanted with the same amount of Cu and Au ( = 3 × 1016 ions cm-2, the second (Au3Cu6) rich in Cu (6 × 1016 Cu+ cm-2 + 3 × 1016 Au+ cm-2) and the third (Au3Cu6H) identical to the second and annealed, after implantation, in a reducing atmosphere [Ar(0.96)H2(0.04) at 1173 K for 1 h].
The fitting to the data was carried out in the framework of the LMA model, as described in the previous section, by using as physical-free parameters the hard-sphere packing fraction , the radius of the metallic sphere Rm, the ratio between the hard-sphere radius RHS and Rm, kc, and the Weibull size distribution parameter b. Added to these, two non-physical free parameters, namely a background and a scaling factor, were also used. An initial guess for kc was derived by knowing that the ratio between the densities of the metals in the glass and bulk form was roughly 6:1 (kc = 1.8). The code uses the MINUIT CERN routine (MINUIT, 1994) to minimize the sum of the squares of the differences between the experimental data and the model; linear or logarithmic scale fits can be chosen. Errors were derived from a function analysis at a 68% confidence level. The results of the analysis are shown in Table 1. Here, , kc, Rm and b are the actual fit parameters, the others (the hard-sphere radius RHS, the metallic-sphere packing fraction , the metallic-sphere interparticle distance D) are derived from these values using the following equations,
We recall that the values for the FCC, BCC, diamond and gas at normal conditions are 0.74, 0.68, 0.34 and 0.0007, respectively. The sequential Au + Cu ion implantation in silica (samples Au3Cu3, Au3Cu6H) determines the formation of Au–Cu alloy nanoclusters (Gonella et al., 1999). The increase of the fluence of the second implanted ion leads to the formation of larger and more separated alloy clusters (see Table 1). There are different facts that explain this result. First, the high-fluence implantation determines an increase in the Cu concentration in the implanted region, thus promoting the dopant aggregation; and second, the irradiation process itself enhances the diffusion of the Cu atoms into the alloy clusters. In the Au3Cu3 sample the alloy clusters are poor in Cu because a significant fraction of Cu atoms are dispersed and oxidized in the matrix (Gonella et al., 1999; Maurizio, Mattei et al., 2003), whereas in the case of higher Cu implantation (sample Au3Cu6) the alloy cluster average composition is rich in Cu (Maurizio, Mattei et al., 2003), indicating that the fraction of Cu atoms that participate in the alloy is higher than in the first case. The effect of the annealing in H2 (4%) atmosphere is to promote the cluster growth; moreover, the clusters after 1 h annealing (sample Au3Cu6) are less packed than in the as-implanted sample (sample Au3Cu6H, see Table 1). This behavior is what is expected from a diffusive growth where larger units grow at the expense of the smaller ones. In this case the annealing induces both a chemical reduction of the Cu atoms dispersed in the matrix and a diffusion and aggregation of the dopant. As a final point, it is worth noting that differences are generally evidenced between the mean cluster size obtained by TEM and GISAXS analyses. These are due not only to the fact that the two techniques investigate different population of clusters, but also to the fact that in the expression (4) of the scattered intensity the number of particles N(R) with size R is weighted with the sixth power of the cluster radius, so that the contribution of the larger particles to the scattering pattern is dominant. For a detailed comparison of GISAXS and TEM results for ion-implanted silica, the reader is referred to Maurizio, Longo et al. (2003).
|
6. Application to neutron scattering
In the formulation of the model presented above, no explicit reference to scattering factors was made owing to the fact that, with this technique, the absolute scattering intensities are difficult to measure. Thus the model uses an arbitrary scale factor to reproduce the data. This also means that data coming from a neutron-scattering experiment could in principle be treated in the same way. A test experiment was carried out in transmission mode, on a 6 × 1016 Ni cm-2 implanted silica glass to check the feasibility of the technique. The Ni implant was chosen because of the relatively large scattering length of Ni (18.5 barns atom−1) compared with those of Si (2.617 barns atom−1) and O (4.232 barns atom−1). Even if the restrictions on the scattering lengths represent a serious limitation to the applicability of this technique, the possibility of performing spin-dependent scattering opens new scenarios in the study of magnetic clusters (de Julian et al., 2001; Maurizio, Longo et al., 2003). The experiment was carried out at the D22 beamline at the Institut Laue Langevin (ILL). The neutron wavelength was = 8 Å and the detector was placed 3 m from the sample. The scattering data were quantitatively analyzed with the LMA model and the fit quality is shown in Fig. 8. We obtained a metallic-sphere radius Rm = 11 Å, an interparticle distance D = 75 Å and a packing fraction = 0.024.
7. Conclusion
In this paper we have presented an introduction to the use of the grazing-incidence small-angle X-ray scattering technique in the morphological study of metal–glass composites obtained by ion implantation. The buried nature of these samples requires a particular optimization of the data-acquisition system in order to collect meaningful data. By suitably choosing the incident angle, integrating the two-dimensional detector data and applying the refraction corrections, scattering functions can be obtained from which quantitative morphological parameters can be extracted. The model for data analysis is based on a Percus–Yevick distribution of polydispersed hard spheres. The size distribution is of the Weibull type and the local monodispersed approximation is used. In the test case of Cu + Au-implanted glasses, parameters like the metallic particle radius, the interparticle distance and packing fraction have been determined. We have shown that the model is also suitable for the analysis of neutron-scattering data, which will open the possibility of performing studies on magnetic structures on magnetic clusters.
Acknowledgements
The authors are indebted to the staff of the ID10 and ID01 beamlines at the ESRF (D. Thiaudiere, D. Smilgies, O. Konovalov, F. Zontone, P. Boesecke and T. H. Metzger) and to R. May of the D22 beamline at ILL for skillful support and fruitful discussions.
References
d'Acapito, F., Mobilio, S., Battaglin, G., Cattaruzza, E., Gonella, F., Caccavale, F., Mazzoldi, P. & Regnard, J. R. (2000). J. Appl. Phys. 87, 1819–1824. Web of Science CrossRef CAS
d'Acapito, F., Thiaudiere, D., Zontone, F. & Regnard, J. R. (1998). Mater. Sci. Forum, 278/281, 891–896.
d'Acapito, F. & Zontone, F. (1999). J. Appl. Cryst. 32, 234–240. Web of Science CrossRef CAS IUCr Journals
Alivisatos, P. (1996). Science, 271, 933–937. CrossRef CAS Web of Science
Ashcroft, N. W. & Lekner, J. (1966). Phys. Rev. 145, 83–90. CrossRef CAS Web of Science
Babonneau, D., Naudon, A., Thiaudiere, D. & Lequien, S. (1999). J. Appl. Cryst. 32, 226–233. Web of Science CrossRef CAS IUCr Journals
Babonneau, D., Videnovic, I. R., Garnier, M. & Oelhafen, P. (2001). Phys. Rev. B, 63, 195401–195405. Web of Science CrossRef
Biersack, J. P. & Haggmark, L. G. (1980). Nucl. Instrum. Methods, 174, 257–269. CrossRef CAS Web of Science
Cattaruzza, E. (2000). Nucl. Instrum. Methods B, 169, 141–155. Web of Science CrossRef CAS
Cattaruzza, E., d'Acapito, F., Gonella, F., Longo, A., Martorana, A., Mattei, G., Maurizio, C. & Thiaudiere, D. (2000). J. Appl. Cryst. 33, 740–743. Web of Science CrossRef CAS IUCr Journals
Desnica, U. V., Dubcek, P., Desnica-Frankovic, I. D., Buljan, M., Salamon, K., Milat, O.., Bernstoff, S. & White, C. W. (2003). Nucl. Instrum. Methods B, 200, 191–195. Web of Science CrossRef CAS
Dubcek, P., Desnica, U. V., Desnica-Frankovic, I. D., Bernstoff, S. & Meldrum, A. (2003). Nucl. Instrum. Methods B, 200, 138–141. CrossRef CAS
Dubcek, P., Pivac, B., Milat, O., Bernstoff, S. & Zulim, I. (2004). Nucl. Instrum. Methods B, 215, 122–128. Web of Science CrossRef CAS
Feynman, R. P., Leighton, R. B. & Sands, M. (1970). The Feynman Lectures on Physics. New York: Addison-Wesley.
Gonella, F., Mattei, G., Mazzoldi, P., Sada, C., Battaglin, G. & Cattaruzza, E. (1999). Appl. Phys. Lett. 75, 55–57. Web of Science CrossRef CAS
Gonella, F. & Mazzoldi, P. (2000). Handbook of Nanostructured Materials and Nanotechnology, Vol. 4. San Diego: Academic Press.
Guinier, A. & Fournet, G. (1955). Small Angle Scattering of X-rays. New York: Wiley.
Haglund, R. (1995). Handbook of Optical Properties, Vol. II. Boca Raton: CRC Press.
Halperlin, W. (1986). Rev. Mod. Phys. 58, 533–606. CrossRef Web of Science
Hammersley, A. P., Svensson, S., Hanfland, M., Fitch, A. & Hausermann, D. (1996). High Press. Res. 14, 235–248. CrossRef Web of Science
Julian, C. de, Sangregorio, C., Mattei, G., Battaglin, G., Cattaruzza, E., Gonella, F., Russo, S. L., D'Orazio, F., Lucari, F., De, G., Gatteschi, D. & Mazzoldi, P. (2001). J. Magn. Magn. Mater. 226/230, 1912.
Kiessig, H. (1931). Ann. Phys. 10, 769–788. CrossRef CAS
Kutsch, B., Lyon, O., Schmitt, M., Mennig, M. & Schmidt, H. (1997). J. Appl. Cryst. 30, 948–956. Web of Science CrossRef CAS IUCr Journals
Lazzari, R. (2002). J. Appl. Cryst. 35, 406–421. Web of Science CrossRef CAS IUCr Journals
Levine, J. R., Cohen, J. B., Chung, Y. W. & Georgopoulos, P. (1989). J. Appl. Cryst. 22, 528–532. CrossRef CAS Web of Science IUCr Journals
Martorana, A., Longo, A., d'Acapito, F., Maurizio, C., Cattaruzza, E. & Gonella, F. (2001). J. Appl. Cryst. 34, 152–156. Web of Science CrossRef CAS IUCr Journals
Mattei, G., Marchi, G. D., Maurizio, C., Mazzoldi, P., Sada, C., Bello, V. & Battaglin, G. (2003). Phys. Rev. Lett. 90, 85502–85505. Web of Science CrossRef CAS
Maurizio, C., Longo, A., Martorana, A., Cattaruzza, E., d'Acapito, F., Gonella, F., de Julian, C., Mattei, G., Mazzoldi, P., Padovani, S. & Boesecke, P. (2003). J. Appl. Cryst. 36, 732–735. Web of Science CrossRef CAS IUCr Journals
Maurizio, C., Mattei, G., Mazzoldi, P., Padovani, S., Cattaruzza, E., Gonella, F., d'Acapito, F. & Zontone, F. (2003). Nucl. Instrum. Methods B, 200, 178–184. Web of Science CrossRef CAS
MINUIT (1994). CERN Program Library Long Writeup D506. CERN, Geneva, Switzerland.
Naudon, A., Babonneau, D., Thiaudiere, D. & Luien, S. (2000). Physica B, 283, 69–74. Web of Science CrossRef CAS
Naudon, A. & Thiaudiere, D. (1997). J. Appl. Cryst. 30, 822–827. Web of Science CrossRef CAS IUCr Journals
Parratt, L. (1954). Phys. Rev. 95, 359–369. CrossRef Web of Science
Pedersen, J. (1993). Phys. Rev. B, 47, 657–665. CrossRef Web of Science
Pedersen, J. (1994). J. Appl. Cryst. 27, 595–608. CrossRef Web of Science IUCr Journals
Percus, J. & Yevick, G. (1958). Phys. Rev. 110, 1–13. CrossRef CAS Web of Science
Polizzi, S., Riello, P., Fagherazzi, G. & Borrelli, N. (1998). J. Non-Cryst. Solids, 232/234, 147–154. Web of Science CrossRef
Renaud, G., Lazzari, R., Revenant, C., Barbier, A., Noblet, M., Ulrich, O., Leroy, F., Jupille, J., Borensztein, Y., Henry, C. R., Deville, J. P., Scheurer, F., Mane-Mane, J. & Fruchart, O. (2003). Science, 300, 1416–1419. Web of Science CrossRef PubMed CAS
Robach, O., Renaud, G. & Barbier, A. (1999). Phys. Rev. B, 60, 5858–5861. Web of Science CrossRef CAS
Stemmer, S., Li, Y., Foran, B., Lysaght, P. S., Streiffer, S. K., Fuoss, P. & Seifert, S. (2003). Appl. Phys. Lett. 83, 3141–3143. Web of Science CrossRef CAS
Windt, D. L. (1998). Comput. Phys. 12, 360–370. Web of Science CrossRef CAS
© 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.