Harry S.
Geddes
a,
Henry D.
Hutchinson
a,
Alex R.
Ha
a,
Nicholas P.
Funnell
b and
Andrew L.
Goodwin
*a
aInorganic Chemistry Laboratory, Department of Chemistry, University of Oxford, South Parks Road, Oxford, OX1 3QR, UK. E-mail: andrew.goodwin@chem.ox.ac.uk
bISIS Facility, Rutherford Appleton Laboratory, Harwell Science and Innovation Campus, Didcot, OX11 0QX, UK
First published on 23rd July 2021
Using a non-negative matrix factorisation (NMF) approach, we show how the pair distribution function (PDF) of complex mixtures can be deconvolved into the contributions from the individual phase components and also the interface between phases. Our focus is on the model system Fe∥Fe3O4. We establish proof-of-concept using idealised PDF data generated from established theory-driven models of the Fe∥Fe3O4 interface. Using X-ray total scattering measurements for corroded Fe samples, and employing our newly-developed NMF analysis, we extract the experimental interface PDF (‘iPDF’) for this same system. We find excellent agreement between theory and experiment. The implications of our results in the broader context of interface characterisation for complex functional materials are discussed.
By their very nature, interfaces are notoriously more difficult to characterise than the bulk phases they connect.9 NMR,10 XAS,11 electron microscopy,12 and vibrational spectroscopy13,14 measurements all offer to varying degrees some experimental sensitivity to interface structure. Yet the process of translating these measurements into a picture of atomic-scale structure remains a significant general challenge.15 It is no surprise then that the field relies heavily on computational methods—including both ab initio13 and empirical16,17 approaches—to develop a collective understanding of interfaces and surface science. This is a mature science, but one that would benefit nonetheless from access to a more diverse range of experimental methods against which its results might be compared.
Here we explore the possibility that pair distribution function (PDF) measurements offer sensitivity to interface structure in a way that is strongly complementary to existing experimental and computational approaches. The PDF is determined experimentally by the inverse Fourier transform of the X-ray or neutron total scattering function—essentially a form of powder diffraction pattern that contains both Bragg and diffuse components.18 It represents a scattering-weighted histogram of interatomic separations. A key advantage over many other experimental techniques is that the function is sensitive at once to both short- and long-range structural correlations.19–22 By convention, the PDFs of complex mixtures are usually interpreted as linear sums of the individual component PDFs;23 such is the basis of the widely-used ‘differential PDF’ approach, for example.24–29 Yet this interpretation remains an approximation because it assumes the interface contribution is weak and incoherent—and so can be neglected.
When might this approximation break down? Consider the simple illustrative case of the Fe∥Fe3O4 interface [Fig. 1(a)] that lies at the heart of corrosion science and is industrially relevant to e.g. the design of pressurised heavy water reactors and industrial descaling processes.30,31 Whereas the conventional (Bragg) diffraction pattern of this two-phase mixture appears at face value to be a straightforward sum of the contributions from crystalline Fe and Fe3O4, its PDF cannot be accounted for entirely in terms of the Fe and Fe3O4 PDFs [Fig. 1(b)]. This is because both intra- and inter-phase correlations contribute, and the latter can be significant if a sufficient fraction of the mixture lies within an appropriate distance of the interface. For example, in a spherical particle of 30 nm diameter—typical for corrosion products32—nearly 35% of atoms lie within 2 nm of the surface; this fraction increases further with surface roughness. Hence composites containing nanocrystalline or amorphous phases—e.g. ferrihydrite33—are an obvious case where the interface contribution might not safely be discounted.
In this study, we show that this interface contribution to the PDF (we call it the ‘iPDF’) can actually be extracted using non-negative matrix factorisation (NMF) analysis34—a method only recently applied to PDF data as a means of deconvolving the function into its individual component factors.35–37 We begin by establishing proof-of-concept, using idealised PDF data generated from suitably-constructed models of the Fe∥Fe3O4 interface. We then report experimental X-ray total scattering measurements for controlled corrosion of polycrystalline Fe powder. We demonstrate, as anticipated, that the PDFs derived from these data contain varying contributions from three components—Fe, Fe3O4, and the Fe∥Fe3O4 iPDF. By extracting the interface PDF component and comparing against that calculated from our synthetic data, we are able to validate the interface structure model proposed on the basis of ab initio DFT calculations in ref. 38. Our paper concludes with a discussion of the relevance of our results to the broader study of interfaces in other families of functional materials.
We used the DFT-relaxed structures of ref. 38 to build a series of atomistic models of the Fe∥Fe3O4 interface with varying Fe/Fe3O4 fractions and layer thicknesses. One step in this process is to align the in-plane lattice constants for the two phases in order to guarantee structural registry. Full details of the various models are given as ESI;† the configuration shown in Fig. 1(a) corresponds to one of these. The X-ray PDFs for all our models, together with those for pure Fe and pure Fe3O4, were calculated using PDFGui,43 transformed to the always-positive G′(r) normalisation of ref. 44 (we omit the prime notation hereafter for ease) and are shown in Fig. 2(a). Least squares fits to these data using linear combinations of the Fe and Fe3O4 PDFs give weighting fractions that deviate from those calculable from the known compositions [Fig. 2(b)], and also a series of difference functions with common features [Fig. 2(c)]. Our hypothesis, which we proceed to test, is that this data set can only be accounted for in terms of varying linear combinations of three PDFs: those of the two bulk-phase PDFs, and also the iPDF.
We use the Metropolis Monte Carlo NMF implementation developed in ref. 35 to deconvolve our synthetic PDF data set into its constituent components. NMF is a nonlinear fitting process that has as its variables the component PDFs and also the component weightings for each input PDF. These variables can be constrained, if so wished, in a variety of sensible ways. For example, it is possible to fix one component PDF to be that of bulk Fe, and/or another to be that of bulk Fe3O4. Likewise the component fractions can be constrained such that the bulk phase PDFs are fitted using a single NMF component. We used various combinations of constraints (including none; see ESI† for further details), and rigorously derived the results shown in Fig. 2(d and e). NMF indeed extracts three components with significant weightings, including the newly-characterised iPDF. The weighting of this component is largest for intermediate compositions; note that asymmetry of the weighting fractions in Fig. 2(e) reflects the difference in X-ray scattering powers of Fe and Fe3O4 phases. In particular, the NMF weights describe the relative contributions to the X-ray PDF at each composition, rather than e.g. the volume fraction of each component.
The Fe∥Fe3O4 iPDF contains features of both the Fe and Fe3O4 PDFs, but is not a linear combination of the two. This is why the iPDF fractions determined in Fig. 2(e) are much larger than the magnitude of the difference function shown in Fig. 1(b) might suggest: when fitting a composite PDF using only bulk phases, the interface contribution is approximated by a sum of the bulk-phase PDFs. By contrast, NMF allows extraction of this iPDF component in its own right. An interesting distinction between the bulk-phase and interface PDFs is that the former are independent of rmax, but the latter (and its weighting) is affected by the real-space range used during the fitting process. This is sensible: the value of rmax essentially defines the interface depth to be characterised and the iPDF will be subtly different for different depths. Some further discussion on this point is given in the ESI.†
Following the very same approach applied to our synthetic data set, we used Metropolis NMF to deconvolve our experimental X-ray PDF data into their constituent parts. Again, we required three NMF components to obtain satisfactory fits to data; the corresponding component fractions are given in Fig. 3(b). Note the very large interface fraction, which implies substantial interface roughness and/or nanoparticulate Fe3O4 corrosion products. Both points are borne out by scanning electron microscopy (SEM) images measured from our samples [Fig. 3(d)].
The experimental iPDF extracted from our X-ray total scattering measurements is shown in Fig. 3(c). To the best of our knowledge, this is the first experimental characterisation of the interatomic correlations at the Fe∥Fe3O4 interface. What is remarkable (in our view) is the extremely close agreement to the iPDF determined from the synthetic data set described above. Hence our data make clear that experiment is entirely consistent with theory in this case—and that the DFT model developed in ref. 38 represents a realistic picture of the interface structure in this industrially important system. In principle, measurements of the type made here might be used to examine the evolution of corrosion product growth. Our data are too coarse to draw any strong conclusions in this regard, but we note the suggestion that the interface fraction grows more rapidly than the bulk Fe3O4 phase fraction early on in the corrosion process. Crucially, these local phase fraction evolutions are related to, but distinct from, the phase fractions determined using conventional analysis of the Bragg scattering intensities.
If, as we show here, the iPDF can be calculated and measured, then in principle it can be used to refine atomistic models of interface structure. Real-space refinement packages such as PDFGui43 or TOPAS47 might straightforwardly be extended to carry out such a refinement. Doing so would bypass the need for NMF altogether: instead the PDF determined for a composite mixture might be interpreted directly in terms of crystallographic models for each component, and a description of the interface orientation. It is in this latter respect that such a refinement would differ from a conventional two-phase analysis of PDF data. There may also be cases where analysis of the interface contribution is more robust in reciprocal—rather than real—space, and we include as ESI† some discussion of the ways in which NMF might be used in this regard.
Our analysis has been understandably simplistic in many respects, and in future one might wish to explore further the role of interface shape, phase size, and even surface effects. Likewise the Fe∥Fe3O4 interface is particularly ‘clean’ in the sense that the atomic structure on either side resembles closely that of the bulk phases. It will be an interesting challenge to test measurement and interpretation of the iPDF for composites where the interface structure is very distinct over large distances; in such cases the r-range over which the iPDF is determined may be particularly important. Nevertheless, we hope to have demonstrated that the interface contribution to PDF data, conventionally discounted, might in fact provide enormously valuable insight into the complex structures of many functional composite materials.
Footnote |
† Electronic supplementary information (ESI) available: Computational modelling of Fe||Fe3O4, PDF calculations, experimental methods, NMF details. See DOI: 10.1039/d1nr01922h |
This journal is © The Royal Society of Chemistry 2021 |