Studies of atomic diffusion in Ni-Pt solid solution by XPCS

Studies of atomic diffusion in Ni-Pt solid solution by X-ray Photon Correlation Spectroscopy

Markus Stana, Michael Leitner, Manuel Ross and Bogdan Sepiol Universität Wien, Fakultät für Physik, Strudlhofgasse 4, 1090 Wien, Austria FRM-II, Technische Universität München, Lichtenbergstraße 1, 85747 Garching, Germany

Atomic scale X-ray Photon Correlation Spectroscopy (aXPCS) was used to study atomic diffusion in the NiPt solid solution with both a single crystal and a polycrystalline sample. Different jump diffusion models are discussed using experimental results and Monte Carlo simulations. The sensitivity of aXPCS experiments to short-range order (in this case governed by a strong Pt-Pt repulsive force) is demonstrated. The activation energy of 2.93(10)eV as well as diffusivities in the range of ms at 830 K agree very well with the results of tracer diffusion studies at much higher temperatures.

66.30.Ny, 78.70.Ck
: J. Phys.: Condens. Matter

1 Introduction

A variety of well established methods can be applied to investigate diffusion in solid state systems [1]. However, studies of diffusion mechanisms on the atomic scale are still challenging. Atomic scale X-ray Photon Correlation Spectroscopy (aXPCS) can provide new insight into this field.

XPCS is the X-ray analogue of the very successful technique of dynamic light scattering or Photon Correlation Spectroscopy [2]. PCS is well established to characterize the equilibrium dynamics of soft condensed matter by determining the intensity autocorrelation function of the scattered laser light as a function of delay time and wave-vector transfer. Similarly, XPCS [3], using coherent synchrotron radiation instead of visible laser light, is being widely used for studies of soft matter dynamics where objects are in the nanometre range (see e.g. [4, 5, 6, 7]). Attention has also been directed towards studies of metal systems as shown by the feasibility test on the (100) Bragg peak of CuAu [8] or of critical dynamics in FeAl [9]. Because of the short X-ray wavelength, XPCS seems to be predestined for studies at atomic length scales. These types of experiments have, however, never been reported in the literature until a successful demonstration of the capability of XPCS to resolve single atomic motion in condensed matter [10]. Particularly, we have clarified the mechanism of atomic dynamics in CuAu [10, 11] and studied dynamics in a metallic glass [12].

The goal of this paper is to describe in detail the atomic motion of dilute platinum atoms in a nickel matrix at relatively low temperatures. The NiPt intermetallic alloy is a perfect candidate for studying diffusion of small quantities of dilute atoms due to its outstanding intensity of diffuse scattering even for low solute concentration levels. Moreover, some profound questions can be addressed regarding interactions between solute atoms and interactions of solute atoms and vacancies. It is well known that vacancies are vehicles of diffusion in alloys [1], but the interactions between vacancy and solute atom can radically influence the motion of the latter. Our objective is to analyse these interactions with the help of aXPCS. These questions were addressed in the past by a number of methods like Mössbauer spectroscopy, nuclear resonant scattering or quasi-elastic neutron scattering [13]. The elementary diffusion steps were investigated for self diffusion [14], for impurity diffusion with consideration to the interaction with vacancies [15], for hydrogen in intermetallic hydrides [16] and in intermetallic B2 [17, 18, 19] and D0 [20] phases. All methods mentioned above have in common that they are limited to only few selected isotopes and have a poor energy resolution which forces measurements at very high temperatures (apart from hydrogen). The time-domain technique of aXPCS is, however, free from these restrictions.

The important aspect of this paper is studying diffusion in single crystal and polycrystalline samples of the same composition. For methods studying spatially resolved atomic motion like aXPCS, single crystal samples provide the possibility to investigate different crystallographic directions and therefore more information. It is, however, sometimes difficult or even impossible to obtain single crystals of satisfactory quality. Showing that studies on polycrystals can, in whole or in part, provide information on a comparable level when using a coherent method, a new field of materials could be investigated with aXPCS. NiPt seems to be a predestined object for such a study.

2 Theory

Atomic structures can be investigated at their elemental scale. Still, monitoring a single atom in a bulk sample on time scales of atomic diffusion it is not possible in real space. It is however possible to measure dynamics in reciprocal space which corresponds to these dynamics in real space.

Usually, reconstruction of real space from a scattering image is very challenging due to the phase problem. Fortunately, one can ignore the details of the exact structure factor, but use the time fluctuations of the intensity pattern to study dynamics of the material with XPCS. This is possible because the intensity autocorrelation function and the autocorrelation function of the electric field are connected by the Siegert relation for Gaussian scattering processes [21]. Therefore it is sufficient to measure the intensity autocorrelation function to gain information about the dynamics in real space


where is the observed intensity at wave-vector and time and the brackets denote an average over the ensemble. In XPCS and in many other types of experiments, the time average of the investigated quantity is calculated from the measured data.

XPCS is the time Fourier-transformed counterpart of coherent quasi-elastic neutron scattering (QENS). It can be shown that the correlation time corresponds to the energy line broadening [22] via . Therefore concepts developed for QENS can be adopted to XPCS.

For a coherent method like aXPCS two atoms of the same species exchanging their positions do not affect in any way the scattering pattern, as the configuration of the coherently illuminated volume before and after the jump is exactly the same. Compared to non-coherent methods it is, however, necessary to account for the short-range order, i.e. for the influence of the surroundings on the movement of each distinct atom. This phenomenon has first been described qualitatively for coherent neutron scattering in liquids by de Gennes [23]. The physical principle for the so-called de Gennes narrowing is based on the fact that peaks in the structure factor of liquids occur at wave-vectors corresponding to the most probable interatomic separation, i.e. they are responsible for the highly correlated and long living atomic arrangements. The same argument holds for the short-range order intensity () in alloys where the higher diffuse intensity between Bragg peaks corresponds to a longer lifetime of certain configurations of solute atoms.

The simplest case of Lorentzian line-shapes with energy-domain methods [13] corresponds to an exponential decay of the intensity autocorrelation function in time


The intensity autocorrelation function for perfectly coherent radiation would decline from a value of 2 to 1. However, because of an only partially coherent X-ray beam in the experimental setup the amplitude is diminished. The correlation time is given by


This follows from the theory describing the short-range order and coherent dynamics in a linear approximation as derived by Sinha and Ross [25]. It is discussed in more detail in Leitner and Vogl [11]. is the average time between jumps of an atom and therefore corresponds to the inverse of the average jump frequency . The probability for an atom to jump to one of the lattice sites of a certain neighbouring shell is given by , with being the number of the appropriate coordinate shell (1 for nearest-neighbour, 2 for next-nearest-neighbour and so on). The vector points from the origin to the th site in the th shell. It can be easily seen that the inverse correlation time is the sum over the inverse correlation times for each shell corresponding to the particular jump mechanism. Therefore it can also be written as .

In a polycrystalline sample one has to average over all orientations of with respect to the crystal lattice


It is essential for the method at hand to be able to fit the data set with a function depending on . As there is no analytical function which satisfies equation (4), an approximation has to be made. A similar approach was introduced by Chudley and Elliott [26] for liquids. Exchanging integration and exponentiation, an average correlation time can be defined as . For this expression an analytical solution can be found ( was also integrated over ):


Using equation (5) the intensity autocorrelation function can be written as:


Numerical calculations show that equation (6) is a very good approximation for equation (4) away from Bragg peaks.

3 Experimental

Nickel and platinum, both of purity, were melted by induction heating in a “cold boat” in argon atmosphere. For the preparation of the polycrystalline sample the resulting ingot was cut and rolled between stainless steel plates down to a thickness of 17 m. It was then pre-annealed at a temperature of 773 K for 22 h in a vacuum of about  mbar. In order to obtain sufficiently large single crystalline domains the rest of the ingot was heat treated at 1573 K for 140 h resulting in crystallites of about about 5 mm diameter. One domain was oriented, cut with a wire saw and mechanically polished to the final thickness of about 18 m. The composition of the sample was analysed by energy-dispersive X-ray spectroscopy (EDX) showing 3.0(5) at.% platinum content.

Figure 1: Autocorrelation function g in Ni-Pt single crystal orientated with direction parallel to the incoming X-ray beam with azimuthal angle at 855 K (upper plot) and at 815 K (lower plot) fitted with equation (2).

The experiment was performed at beamline ID10A of the ESRF using a standard coherent setup with 8 keV photons at a monochromaticity of . Rollerblade slits were used to select an m beam section and speckle patterns were recorded with a direct-illumination CCD camera with pixels at a size of m. Prior experience has shown that the CCD camera placed close to the sample (about 0.6 m) is advantageous to maximize the measured signal. It was mounted in a way that allowed it to be moved along a scattering angle and an azimuthal angle . Note that in a scattering experiment the length of the scattering vector is connected to the scattering angle via , with being the wave vector of the incoming beam. As one can move freely on the Ewald sphere in reciprocal space by changing the detector position in real space, the intensity autocorrelation function was measured for several . In the polycrystalline sample only the scattering angle 2 (i.e. ) can influence the correlation time. Therefore measurements at arbitrary azimuthal angles are equivalent. For the single crystal sample an additional second scan for fixed and different azimuthal angles was performed. Measurements were performed in transmission geometry in a resistively heated vacuum furnace ( mbar). The temperature was stabilized by a PID temperature controller within a range of 0.1 K. The measurement time at a fixed scattering vector was between 20 and 40 minutes at a frame rate of one frame every 10.9 seconds (exposure plus readout time). Due to the limited intensity of the synchrotron source and low scattering intensity at large scattering angles in the diffuse regime it is beneficial to make the data collected by the CCD camera subject to the so-called droplet algorithm [27] to detect single photon events. Our implementation of the droplet algorithm is described in [24]. These corrected images were evaluated according to equation (1), finding the correlations in the fluctuating speckle intensities (see figure 1) where the autocorrelation function was averaged over all pixels to obtain sufficient statistical accuracy.

All measurements were carried out in the diffuse regime, away from Bragg peaks. Any kind of strain broadening / asymmetry resulting from dislocations, planar faults etc. are constrained to -regions very close to the Bragg reflections [28]. Therefore such effects do not influence the intensities obtained in an aXPCS experiment. The number of vacancies relative to the number of atoms in Ni at 830 K is in the order of [29]. Therefore the total scattering power of the vacancies and their induced displacements is negligible. An estimation using the Debye model at this temperature gives a contribution in the order of 10 originating from thermal diffuse scattering. This contribution can therefore also be neglected. As a result the diffuse intensity is only due to the stochastic occupations of lattice sites by Ni or Pt, and the measured correlation decay therefore corresponds to exchanges of atoms.

Small absolute values of reciprocal vectors correspond to very long correlation times. Instabilities of the experimental setup result in additional intensity variations. This can lead to fault when overlapping with the slowly decaying autocorrelation function caused by the sample itself but is not important if the time-scale of both processes is clearly different. Correlation times measured at very small 2 angles were therefore not used for the evaluation of diffusion models.

It is vitally important for aXPCS measurements to be performed in thermal equilibrium. To ensure that order relaxation is negligible, multiple measurements were made repeatedly at different detector positions for each sample. Usually the sample relaxes faster when starting from a higher temperature and shows relaxation times of about one hour. After that time the single crystal stayed stable during the whole experiment. The polycrystalline sample, however seemingly relaxed again after about 4 h at 830 K, leading to still longer correlation times.

Note that the correlation time as measured in an aXPCS experiment quantifies the timescale over which the atomic arrangement changes by a sizeable amount. Even though effects such as enhanced grain boundary diffusion can affect the mass transport over macroscopic distances significantly, the contribution of the grain boundaries to the scattered signal and therefore to the correlation decay is negligible. As a consequence, aXPCS is inherently sensitive only to bulk diffusion.

To investigate whether there were crystallite growth effects in the polycrystalline sample during the experiment, the sample was studied with electron microscopy. Unfortunately the results were inconclusive. Therefore we assumed that such effects can be neglected as long as no relaxation in the correlation time is detectable.

4 Monte Carlo simulations

To simulate a simple diffusion model, Monte Carlo simulation of diffusion is a valuable approach. Two kinds of atoms were arranged in a proportion of over an fcc-type lattice with periodic boundary conditions. One vacancy was introduced by randomly emptying a lattice site. The simulation consists in iteratively choosing one of the 12 lattice sites of the vacancy’s first coordination shell randomly, then, according to the Metropolis algorithm, calculating the probability , where denote the system energies calculated with pair potentials for nearest-neighbours before () and after () the jump, respectively, where is the Boltzmann constant and the absolute temperature. A decision to perform a jump or not was taken by comparing with a random number. A unit time (MC step) was defined as a number of trials equal to the number of lattice sites. After equilibrating the system, the actual diffusion simulation was started. After every half Monte Carlo step a discrete Fourier transform of the lattice was carried out and the amplitudes were used to calculate the scattering intensity and the amplitude autocorrelation function , t) for certain values. Finally, using Siegerts relation, the intensity autocorrelation function was calculated and compared with the experimental correlation time.

The same method was used to calculate short-range order intensities from a reciprocal space map of 64 voxels with an atom configuration of 432 atoms averaging the intensities over 10000 Fourier transforms.

5 Results

In order to distinguish between different diffusion models, measurements of intensity autocorrelation times in polycrystalline and/or in single crystal samples alone are insufficient. As shown in section 2 short-range order intensity plays a crucial role in the interpretation of the measured correlation times.

5.1 Short-range order

The particular importance of local atomic arrangements in alloys for the interpretation of aXPCS measurements follows from equation (3). Detailed data about SRO like in the case of CuAu [30] is rarely available. SRO must be actually measured for exactly the same alloy composition and temperature as the sample under study. It is very difficult to deduce short-range order data from experimentally measured intensity data. A detailed description of the diffuse scattering theory can be found in several articles and books (see e.g. [31, 32] and literature cited therein). The factor which actually makes the analysis of diffuse scattering so difficult and time consuming is the separation of SRO scattering (which is of interest to us) and atomic displacement scattering. Simulations showed that in a system with a very high scattering contrast and low solute atom concentration like NiPt, the magnitude of those two terms can be of the same order. Also terms of higher order in the expansion of diffuse intensity, e.g. caused by static displacements (Huang scattering) or thermal diffuse scattering, can play a role.

Phase diagrams of the Ni-Pt system [33] show three ordered equilibrium phases similar to Cu-Au: NiPt (L1 structure), NiPt, and NiPt (both L1 structures). There was a weak tendency found for a CuAu-type ordering for NiPt with composition [34] which is also indicated by the diffuse intensity for as shown in the inset in figure 2. There is, however, no data for short-range order configuration on a system with available in the literature. Note that our main goal remains an atomistic diffusion mechanism and the details of the diffuse scattering intensity are not vitally important. Therefore we choose two simple models: (m1) a model with no interaction between the atoms and (m2) a model with strong Pt-Pt repulsing force. The latter was chosen because a L1 ordering type tendency resulting from strong Pt-Pt nearest-neighbour repulsion was assumed to be more likely than from Pt-Pt next-nearest-neighbour attraction or further ranging interactions 111A more thorough investigation of the interaction parameters in the system NiPt based on ab initio calculation will be subject to a following paper.. The two model systems m1 and m2 were simulated using a Metropolis Monte Carlo algorithm (see section 4). From the atomic configuration we calculated the Warren-Cowley parameters [35]. In case of 3% solute atoms with strong repulsive force the calculated parameters were , and . These parameters can be used to calculate the short-range order intensity in Laue units according to


where again is the family of vectors for a certain shell as discussed in section 2. As shown in figure 2 the so calculated short-range order intensity is in good agreement with the intensities calculated directly from the Fourier transform of the real-space configuration.

Figure 2: For a single crystal in orientation with azimuthal angle the short-range order intensity was calculated without interaction between the atoms (red triangles) and for strong Pt-Pt repulsion (blue disks). For the latter the short-range order was also calculated from Warren-Cowley parameters gained from the simulation (grey line). Inset shows measured diffuse intensity as a function of the wave-vector transfer.

5.2 Diffusion mechanism

The main objective of this paper is to discuss the mechanism which governs diffusion in a Ni-Pt intermetallic alloy. An essential advantage of a microscopic (atomistic) method in diffusion investigation is the possibility to identify such a diffusion mechanism. In order to verify which diffusion mechanism operates in a specific intermetallic phase, it is necessary to develop a mathematical model for the diffusion process. This mathematical idea is generated from the knowledge of the system or it is simply a right guess. The mathematical model can be (at least in certain cases) solved analytically and yields the wave-vector dependent correlation time in a single crystal sample (see equation (3)), or the correlation time dependent on the magnitude of the wave-vector in a polycrystalline sample (see equation (5)). A model is determined by the jump probability parameters . The only free parameter fitted to the experimental data is the inverse average residence time which corresponds to the jump frequency.

Needless to say it is not particularly difficult to predict the atomistic mechanism of diffusion in Ni-Pt. It is simply an exchange of a nickel or platinum atom with a vacancy. As the scattering length of platinum considerably exceeds that of nickel and vacancies are very rare in the alloy at the temperatures of measurement, the measurable speckle fluctuations result predominantly from the exchanges of platinum atoms with nickel atoms. This very simple picture enables studies of further subtleties of the diffusion mechanism, such as the influence of interaction between solute atoms or between solute atoms and vacancies.

A simple and useful tool when discussing a diffusion model involving a vacancy is the so-called encounter approximation [36]. It has been successfully applied in the interpretation of diffusion in nuclear magnetic resonance and in Mössbauer spectroscopy studies[17, 18]. An encounter is here defined as the sum of all exchanges of one and the same atom with one and the same vacancy. As the number of vacancies is much smaller than the number of atoms, the time between one and the same vacancy jumping twice is much smaller than the time between consecutive jumps of one and the same atom. Details during an encounter are thus “invisible” for aXPCS and only a net result of one or more exchanges with a vacancy, which can carry an atom to nearest or to further neighbour shells, is accounted. One can assume that encounters are mutually independent and calculate e.g. diffusivity as a simple random walk of the vacancy.

In an fcc lattice the first neighbour shell is at a distance of , the second of , the third of and fourth of , where is the lattice constant. The nearest-neighbour shell can be reached by only one exchange between a solute atom and a vacancy. The second, third and fourth shells can, however all be reached after two exchanges with a vacancy. Still farther shells () need three or more exchanges and can therefore be neglected. Also encounters that do not lead to an effective jump are invisible for an aXPCS experiment and are therefore neglected from here on ( = 0). We used this to re-normalize the jump probabilities for self diffusion given in [37] and gained jump probabilities of: = 92.6 %, = 2.4 %, = 3.9 % and = 0.8 % (the neglected probability for further jumps sums up to = 0.3 %). This encounter model (m1), however, describes only systems with negligible interactions between matrix atoms, solvent atoms and vacancies. Therefore we extend the encounter model according to equation (3), using the SRO intensity for strong Pt-Pt repulsion as discussed in section 5.1. We use this as a second model for the description of Pt diffusion and call it modified encounter model (m2).

In order to distinguish between the two models, correlation times were measured at a constant temperature of 830 K and at different detector positions corresponding to different and for the single crystal and the polycrystalline sample respectively.

(a) Radial scan of single crystal sample
(b) Azimuthal scan of single crystal sample
Figure 3: Reciprocal correlation times in single crystal sample orientated in direction at 830 K. The direction lies at = . The radial scan (a) was taken at = . The azimuthal scan given in (b) was taken at 2 = . Both datasets were fitted simultaneously with the classical encounter model with no short-range order (m1) and with the modified encounter model with correction for the short-range order calculated for strong Pt-Pt repulsion (m2). For comparison, a model with only next-nearest-neighbour jumps (m3) corrected for SRO is also shown. Monte Carlo simulations without interaction potential (s1), with Pt-Pt nearest-neighbour repulsion (s2) and, again for comparison, with Pt-Pt attraction (s3) were carried out.

MC simulations sites for different interaction energies between platinum atoms were carried out. One for no interaction between the atoms (s1), one with strong repulsive potential between Pt atoms (s2) and one with an attractive potential between Pt atoms (s3) within the first neighbour shell ().

The experimentally measured data as well as the data from MC simulations and the models for a diffusion mechanism in the single crystal are shown in figure 3. Note that the only adjustment parameter is the jump frequency, which is simply a multiplicative constant. The figure shows that the MC simulation with no interaction potential (s1) complies with the encounter model for a system with no SRO (m1). The model with jump probabilities according to the encounter model and short-range order correction according to equation (3) (m2) gives the same values as the simulated data for strong Pt-Pt repulsion (s2) and describes the experimental data sufficiently well. For this model the fit produces a jump frequency of s. Unfortunately the quality of the measured data was insufficient to draw sound conclusions about solute-vacancy interaction. Obviously a model with only next-nearest-neighbour jumps ( = 100 %, = 0) (m3) as well as a model with attractive force between Pt atoms (s3) is wrong.

The polycrystalline sample was fitted by equation (6) and with the same jump probabilities as in the case of the single crystal. The short-range order correction was calculated by numerically integrating equation (7) over all orientations of . As mentioned in section 3, this sample was only quasi-stable for a limited time. It can therefore not be assumed that a sufficient number of small crystallites was covered by the beam during the measurement. For averaging over a small number of crystallites with unknown orientation the Chudley-Elliott model [26] is, however, still the best guess. The calculated jump frequency using the modified encounter model as shown in figure 4 was s. The difference to the single crystal value arises most likely from a difference in sample temperature as a different furnace with a thermocouple at a different position with respect to the sample was used for this measurement. Again a model with only next-nearest neighbour jumps (m3) can obviously not appropriately describe the data measured.

Figure 4: Reciprocal correlation times in polycrystalline sample at 830 K fitted with the classical encounter model (m1), modified encounter model (m2) and for comparison a model with only next-nearest-neighbour jumps (m3).

5.3 Activation energy and diffusion constant

The activation energy is the sum of a vacancy formation and migration energy. The relationship between diffusivity, temperature and activation energy can be described fairly accurately by the Arrhenius law [1]. A number of measurements were performed at different sample temperatures. Again measurements for different temperatures where repeated in order to guarantee thermal equilibrium of the sample.

Figure 5: Chemical diffusion coefficients for polycrystalline (□) and single crystal sample (▽) compared to results from tracer measurements [38] (○). The green line (dashed) represents the fit for the single crystal sample, the blue line (solid) represents polycrystalline data and the red line (dash-dotted) represents the results of the tracer experiment.

Using the Einstein relation [39] and the sum of partial diffusion coefficients into farther shells one can write:


where is the radius of neighbour shell . For fixed values of one can use equation (3) to calculate the residence time from the measured correlation times for each temperature. The lattice constant for NiPt was calculated using Vegard’s law [40] and corrected for the thermal expansion at 830 K ( Å). With this value and using the modified encounter model the chemical diffusion coefficients can be calculated for each temperature as shown in figure 5.

One should notice that a diffusion coefficient measured by the aXPCS technique is, owing to the coherency of the method, the chemical diffusion coefficient . As already emphasized, exchanges between atoms of the same species do not influence the speckle pattern and are not recorded. However, in the case of our particular Ni-Pt solid-solution sample, chemical diffusion and self-diffusion coefficients of platinum are equal in a first approximation. This is simply a result of a low concentration of platinum atoms which makes direct exchanges between them highly improbable.

Using the Arrhenius law the activation energy as well as the pre-exponential factor can be calculated. Calculated activation energies for the single crystal and the polycrystalline sample, respectively were very similar with  eV and  eV. These values are in excellent agreement with literature values for interdiffusion [38], where the activation energy measured at much higher temperatures (between 900 and 1300C) is (291.2 3.7) kJ/mol = (3.02 0.04) eV. The good agreement between our measurements and the data given in [38] allows us to conclude that there was no grain boundary diffusion contribution in the data given in that reference (apart from the point measured at the lowest temperature). The calculated pre-exponential factors are  = 1.5 (+4.6/-1.1) ms for the single crystal and  = 3.7 (+9.2/-3.4) ms for the polycrystalline data. Again the value of  = 9.2 (+3.4/-2.5) ms found in [38] is in agreement with our results.

6 Conclusions

We could show that aXPCS allows the observation of low concentrations of solute atoms if the electronic form-factor contrast between solvent and solute atoms is high enough. On the other hand it is the question of the coherent synchrotron brilliance to be able to measure still more dilute samples. Correlated diffusion according to the encounter model with correction for short-range order effects is a simple atomic diffusion process and the measured factors are reliable and in very good agreement with the results of tracer methods. Diffusion rates measured by aXPCS are, however, in the order of 10ms and therefore on the lower border of what can be reached by the best tracer measurements. For measuring diffusivities this low, aXPCS is the only method capable of determining atomic jump vectors.

The Pt-Pt repulsion indicated by the Monte Carlo simulation is in agreement with the NiPt phase [33] for low platinum concentration. Although for an atomic concentration as small as 3% one can not speak of a NiPt phase, the characteristic atomic interactions which lead to its formation are present. It was, however, due to insufficient accuracy not yet possible to draw conclusions about solute-vacancy interaction from the comparison with simulations. With brighter synchrotron sources and better statistical significance coming along with them, it should be possible to investigate these features in the near future. Using Monte Carlo simulations on the bases of atomic interaction potentials and/or cluster expansion, aXPCS shows promise in being a valid tool to confirm and test ab initio calculated potentials.

X-ray photon correlation spectroscopy allows one to draw conclusions about the dominating diffusion mechanism. To distinguish subtle differences e.g. whether the encounter model or a model with slightly different jump probabilities takes place requires much better statistics than we achieved in our experiment.

The experimental data from the polycrystalline sample showed generally good agreement in activation energy and diffusion rates with the tracer method. It seems also possible to reveal the dominating diffusion mechanism (in this case nearest-neighbour jumps) with these data. In order to prepare polycrystalline samples allowing for exact measurements of diffusion rates, efforts have to be made to prepare samples with sufficiently small crystallites and to eliminate crystallite growth effects. Especially aXPCS experiments on powder samples seem to be particularly promising in this regard.

Overall we conclude that aXPCS is a very reliable method for investigating atomic diffusion for all types of alloys and provides a reasonable reference for ab initio simulations.

This research was funded by the Austrian Science Fund (FWF) contract P-22402. We thank F. Zontone and the whole team of beamline ID10A at the ESRF in Grenoble, S. Puchegger from Faculty Center for Nano Structure Research University of Vienna for the EDX analysis, D. Geist and C. Rentenberger from Physics of Nanostructured Materials for TEM analysis, J. Akbarzadeh for her help with the SAXS and M. Rohrer for support with metallography. Furthermore our thanks go to F. Gröstlinger for his assistance during beamtimes and to Martin Leitner for preliminary ab initio calculations of the interaction constants. Also we thank H. Sassik from the Technical University of Vienna for providing the cold boat.



  • [1] H. Mehrer. Diffusion in solids: fundamentals, methods, materials, diffusion-controlled processes. Springer, New York, 2007.
  • [2] W Brown. Dynamic light scattering: the method and some applications, volume 49. Oxford University Press, USA, 1993.
  • [3] Mark Sutton. A review of X-ray intensity fluctuation spectroscopy. C. R. Physique, 9(5-6):657–667, 2008.
  • [4] G Grübel, A Madsen, and A Robert. X-Ray Photon Correlation Spectroscopy (XPCS). In R Borsali and R. Pecora, editors, Soft-Matter Characterization, pages 935–995. Springer, 2008.
  • [5] X Lu, S G J Mochrie, S Narayanan, A R Sandy, and M Sprung. Temperature-dependent structural arrest of silica colloids in a water–lutidine binary mixture. Soft Matter, 6(24):6160–6177, 2010.
  • [6] L Cipelletti and E C Weeks. Glassy dynamics and dynamical heterogeneity in colloids. In L. Berthier, G. Biroli, J.-P. Bouchaud, L. Cipelletti, and W. van Saarlos, editors, Dynamical Heterogeneities in Glasses, Colloids and Granular Media, pages 110–151. Oxford University Press, 2011.
  • [7] Y Chushkin, C Caronna, and A Madsen. A novel event correlation scheme for X-ray photon correlation spectroscopy. J. Appl. Crystallogr., 45(Part 4):807–813, 2012.
  • [8] M Sutton, S G J Mochrie, T Greytak, S E Nagler, L E Berman, G A Held, and G B Stephenson. Observation of Speckle by Diffraction with Coherent X-Rays. Nature, 352(6336):608–610, 1991.
  • [9] S Brauer, G B Stephenson, M Sutton, R Bruning, E Dufresne, S G J Mochrie, G Grübel, J Als-Nielsen, and D L Abernathy. X-Ray Intensity Fluctuation Spectroscopy Observations of Critical Dynamics in FeAl. Phys. Rev. Lett., 74(11):2010–2013, 1995.
  • [10] M Leitner, B Sepiol, L-M Stadler, B Pfau, and G Vogl. Atomic diffusion studied with coherent X-rays. Nature Materials, 8(9):717–720, 2009.
  • [11] M Leitner and G Vogl. Quasi-elastic scattering under short-range order: the linear regime and beyond. J. Phys.-Condens. Mat., 23:254206, 2011.
  • [12] Michael Leitner, Bogdan Sepiol, Lorenz-Mathias Stadler, and Bastian Pfau. Time-resolved study of the crystallization dynamics in a metallic glass. Phys. Rev. B, 86:064202, Aug 2012.
  • [13] G Vogl and B Sepiol. The Elementary Diffusion Step in Metals Studied by the Interefrence of Gamma-Rays, X-Rays and Neutrons. In P. Heitjans and J. Kärger, editors, Diffusion in Condensed Matter, pages 65–91. Springer, 2005.
  • [14] G Vogl, W Petry, Th. Flottmann, and A Heiming. Direct determination of the self-diffusion mechanism in bcc beta -titanium. Phys. Rev. B, 39(8):5025–5034, 1989.
  • [15] K H Steinmetz, G Vogl, W Petry, and K Schroeder. Diffusion of iron in copper studied by Mössbauer spectroscopy on single crystals. Phys. Rev. B, 34(1):107–116, 1986.
  • [16] A V Skripov, T J Udovic, J C Cook, R Hempelmann, A A Rempel, and A I Gusev. Quasielastic neutron scattering study of hydrogen motion in NbC0.71H0.28. J. Phys.-Condens. Mat., 21(17), 2009.
  • [17] R Feldwisch, B Sepiol, and G Vogl. Elementary diffusion jump of iron atoms in intermetallic phases studied by Mössbauer spectroscopy - I. Fe-Al close to equiatomic stoichiometry. Acta metall. mater., 42(9):3175 – 3181, 1994.
  • [18] R Feldwisch, B Sepiol, and G Vogl. Elementary diffusion jump of iron atoms in intermetallic phases studied by Mössbauer spectroscopy - II. from order to disorder. Acta metall. mater., 43(5):2033–2039, 1995.
  • [19] M Kaisermayr, J Combet, H Ipser, H Schicketanz, B Sepiol, and G Vogl. Nickel diffusion in B2-NiGa studied with quasielastic neutron scattering. Phys. Rev. B, 61(18):12038–12044, MAY 1 2000.
  • [20] B Sepiol, A Meyer, G Vogl, H Franz, and R Rüffer. Diffusion in a crystal lattice with nuclear resonant scattering of synchrotron radiation. Phys. Rev. B, 57(17):10433–10439, 1998.
  • [21] P A Lemieux and D J Durian. Investigating non-Gaussian scattering processes by using nth-order intensity correlation functions. J. Opt. Soc. Am., 16(7):1651–1664, 1999.
  • [22] R Hempelmann. Quasielastic Neutron Scattering and Solid State Diffusion. Oxford University Press, Oxford, 2000.
  • [23] P G De Gennes. Liquid dynamics and inelastic scattering of neutrons. Physica, 25(7-12):825–839, 1959.
  • [24] M Leitner. Studying Atomic Dynamics with Coherent X-rays. Springer, 2012.
  • [25] S K Sinha and D K Ross. Self-consistent density response function method for dynamics of light interstitials in crystals. Physica B+C, 149(1–3):51 – 56, 1988.
  • [26] C T Chudley and R J Elliott. Neutron scattering from a liquid on a jump diffusion model. Proc. Phys. Soc., 77:353, 1961.
  • [27] F Livet, F Bley, J Mainville, R Caudron, S G J Mochrie, E Geissler, G Dolino, D Abernathy, G Grübel, and M Sutton. Using direct illumination CCDs as high-resolution area detectors for X-ray scattering. Nucl. Instrum. Meth. A, 451(3):596–609, 2000.
  • [28] M.A. Krivoglaz. X-ray and neutron diffraction in nonideal crystals. Springer, 1996.
  • [29] A.M. Monti M. De Koning, S. Ramos de Debiaggi. Vacancy-Formation Thermodynamics in Aluminium and Nickel: a Computational Study. Defect Diffus Forum, 224–225:59–74, 2003.
  • [30] B Schönfeld, M J Portmann, S Y Yu, and G Kostorz. The type of order in Cu-10 at.% Au–evidence from the diffuse scattering of X-rays. Acta mater., 47(5):1413–1416, 1999.
  • [31] L H Schwartz and J B Cohen. Diffraction from materials. Springer, New York, 1987.
  • [32] B Schönfeld. Local atomic arrangements in binary alloys. Prog. Mater. Sci., 44(5):435 – 543, 1999.
  • [33] M C Cadeville, C E Dahmani, and F Kern. Magnetism and spatial order in Ni-Pt and Co-Pt alloys. J. Magn. Magn. Mater., 54:1055–1056, 1986.
  • [34] R E Parra and J W Cable. Neutron study of magnetic-moment distribution in Ni-Pt alloys. Phys. Rev. B, 21:5494–5504, Jun 1980.
  • [35] J M Cowley. Short- and Long-Range Order Parameters in Disordered Solid Solutions. Phys. Rev., 120:1648–1657, Dec2 1960.
  • [36] D Wolf. Theory of Mössbauer line broadening due to correlated diffusion in crystals. Appl. Phys. Lett., 30(12):617–619, 1977.
  • [37] C A Sholl. Diffusion correlation factors and atomic displacements for the vacancy mechanism. J. Phys. C, 14:2723–2729, 1981.
  • [38] M S A Karunaratne and R C Reed. Interdiffusion of the platinum-group metals in nickel at elevated temperatures. Acta mater., 51(10):2905–2919, 2003.
  • [39] A Einstein. Über die von der molekularkinetischen Theorie der Wärme geforderte Bewegung von in ruhenden Flüssigkeiten suspendierten Teilchen. Ann. Phys., 17:549–560, 1905.
  • [40] L. Vegard. Die Konstitution der Mischkristalle und die Raumfüllung der Atome. Z. Phys, 5(1):17–26, 1921.
Comments 0
Request Comment
You are adding the first comment!
How to quickly get a good reply:
  • Give credit where it’s due by listing out the positive aspects of a paper before getting into which changes should be made.
  • Be specific in your critique, and provide supporting evidence with appropriate references to substantiate general statements.
  • Your comment should inspire ideas to flow and help the author improves the paper.

The better we are at sharing our knowledge with each other, the faster we move forward.
The feedback must be of minimum 40 characters and the title a minimum of 5 characters
Add comment
Loading ...
This is a comment super asjknd jkasnjk adsnkj
The feedback must be of minumum 40 characters
The feedback must be of minumum 40 characters

You are asking your first question!
How to quickly get a good answer:
  • Keep your question short and to the point
  • Check for grammar or spelling errors.
  • Phrase it like a question
Test description