Abstract
Kaolin M1 and M2 studied by X-ray diffraction focus on the mullite phase, which is the main phase present in both products. The Williamson–Hall and Warren–Averbach methods for determining the crystallite size and microstrains of integral breadth β are calculated by the FullProf program. The integral breadth (β) is a mixture resulting from the microstrains and size effect, so this should be taken into account during the calculation. The Williamson–Hall chart determines whether the sample is affected by grain size or microstrain. It appears very clearly that the principal phase of the various sintered kaolins, mullite, is free from internal microstrains. It is the case of the mixtures fritted at low temperature (1200 °C) during 1 h and also the case of the mixtures of the type chamotte cooks with 1350 °C during very long times (several weeks). This result is very significant as it gives an element of explanation to a very significant quality of mullite: its mechanical resistance during uses at high temperature remains.
Introduction
When recording a diffraction pattern using powders, it is common to observe that the shape or width of the lines do not correspond to the instrumental resolution function that would be measured, for example, using a reference sample. This additional intrinsic enlargement is due to the microstructure of the sample and comes mainly from two effects size and microstrain.
To find the true profile of the sample f (x), a deconvolution operation is used. Deconvolution is an operation that corrects the observed line patterns of instrumental broadening provides the function f (2θ).
In addition to the analytical methods that use rigorous deconvolution, simple methods are often used to separate the widths of the convoluted functions. The intrinsic error of these methods has been treated in particular by several methods, as related by authors.1–3 To be reliable, the quality of the recording, with sufficiently long times, is required. This makes it difficult to study low symmetry compounds where the overlap of the lines is greater. Moreover, using these methods, the deconvolution is always carried out based on the hypotheses of the profiles of the components involved. The microstructural results can differ significantly from one hypothesis to another.1–15
Methods of microstructural analysis
The intrinsic diffraction line pattern is mainly caused by the size and internal stresses of the sample. The methods used for the evaluation of these two phenomena are as below.
The Scherrer method is used when the integral breadth (β) results only from the crystallite size effect (when the matter does not contain microstrains). The grain size is given by the well-known Scherrer relation (equation (1)) 16
For widths specified in 2
For spherical particles, it is generally given as (equation (2))
For widths in units where reciprocal θ is the Bragg angle, Df is the average grain size, λ is the wavelength, and β* is the integral width characterizing enlargement, expressed as the reciprocal lattice.
The Warren–Averbach method
If the sample does not contain a microstrain, then
∆s = λ∕2 (sinθ2 sinθ1) and θ1 and θ2 are the values of the limit angles that define the profile.

Fourier coefficients
In Figure 2, the point of intersection of the initial slope of the curve with the x-axis (L) provides the size Df.
The Williamson–Hall method
This method applies to Voightian or pseudo-Voightian profile types. It uses the commonly called Willamson–Hall plot. 20 This graph, as a straight line, gives valuable indications on the presence or absence of microstresses in the sample. The slope of the line shows the presence of the latter (Figure 2).

Williamson–Hall plot.
If the profiles are purely Gaussian or purely Voightian, this method provides a very simple and quick way to evaluate the effects of stress and size. For the other pseudo-Voightian profiles, approximations are necessary for the separation of these two effects.21,22
Results and discussion
Since the data is obtained from point-by-point survey, they do not correspond to that obtained by the Stokes method, which works on constant steps. For this reason, we have written a program that performs this adaptation of the data. The profile parameters obtained for the different kaolins are recorded in Tables 1 and 2.
Parameters of the X-ray diffraction profiles for kaolin DD at 1200 °C.
Parameters of the X-ray diffraction profiles for KT2 kaolin at 1350 °C.
In order to adapt the standard, a least squares method has been established in the laboratory. This makes it possible to approach the widths of the diffraction peaks according to the Caglioti formula.
Moreover, from the hypothesis of a pseudo-Voightian shape of quartz peaks, we automated obtaining the standard corresponding to each peak of the studied sample. This makes it much easier to use the standard in analyses of different samples. It is therefore interesting to use this second method for microstructural analysis.
The results of the application of the Stokes method were evaluated each time by comparing the calculated profiles with the experimental profiles. Figures 3 and 4 show the true profile and the experimental profile for a peak chosen from the diffractogram of each kaolin DD studied. As can be seen in these figures, the calculated and experimental profiles are sufficiently superposed. From the results obtained, the microstructural parameters are deduced.

Profile of the peak (001) of kaolin DD at 1200 °C.

Profile of the peak (210) of kaolin DD at 1200 °C.
The application of the Stokes method provides the profile of the true function. Figures 5 and 6 show two profiles corresponding to the peaks, respectively, at 2θ = 16.36° and 2θ = 31°. In these latter figures and taking into account the difference introduced by the intrinsic effect of the sample, it is easy to see that the calculated and experimental profiles are sufficiently superimposed.

Profile of the KT2 kaolin peak (100).

Profile of the KT2 kaolin peak (110).
Williamson–Hall plot
The Williamson–Hall diagram of each sample of kaolin DD is given in Figure 7. Note that the curve is a straight line with an almost zero slope. The negative value of this slope is probably due to the accuracy of the data obtained from the diffraction and fitting data. It can be deduced that the sample has no microstrains.

Williamson–Hall plot of kaolin DD at 1200 °C.
Figure 8 shows the Williamson–Hall diagram of kaolin KT2. It was observed that the curve is a straight line with a negligible and negative slope.

Williamson–Hall plot of KT2 kaolin.
Evaluation of kaolin’s DD crystallite size
Since the sample does not have any microstrains, the Warren–Averbach method was applied to a non-distorted case. Figures 9 and 10 show the variation of the Fourier coefficients as a function of L (harmonic number) for two peaks of kaolin obtained at 1200 °C (peaks (100) and (001)). The slope originally obtained by extrapolation provides an average value for the crystallite size in the direction perpendicular to the diffractive plane. So (equation (4))
where L harmonic number is
with
2θ1 and 2θ2 are the values of the limit angles that define the profile.

Variation of the Fourier coefficients of the peak (001) as a function of L (DD 1200 °C).

Variation of the Fourier coefficients of the peak (210) as a function of L (DD 1200 °C).
The values obtained for the size of the crystallites of the different kaolin DD studied are reported in Table 3. As can be seen, the crystallites are of quasi-nanometric dimensions. This result is consistent with observations on other DD kaolin. The values found can be somewhat confirmed by the Williamson–Hall plot. Indeed the average size can be evaluated as (equation (5))
where β* is the integral width in reciprocal coordinates. For kaolin at 1200 °C, the average value found from the Williamson–Hall plot is about 190 Å (compared with 114 Å found by the Warren–Averbach method.
The crystallite sizes of kaolin DD obtained at 1200 °C.
Figures 11 and 12 show the variation of the Fourier coefficients as a function of L (harmonic number).

Variation of the Fourier coefficients of the peak (100) as a function of L (KT2).

Variation of the Fourier coefficients of the peak (001) as a function of L (KT2).
The average value of the size of the crystallites in the direction perpendicular to the diffracting plane is provided by the slope at the origin.
The values obtained for the average crystallite size are reported in Table 4.
Average size of KT2 crystallites from different diffraction peaks.
As can be seen in Table 4, the dimensions of the crystallites are of nanometric order.
To confirm, to some extent, the values found, the average crystallite size from the Williamson–Hall diagram was evaluated. The value found is 230 Å, which is close with the average size determined using the Warren–Averbach method (245 Å).
Size distribution in the main phase of kaolin DD and KT2
A study of the size distribution within the main phase of kaolin DD was performed. This distribution was calculated using a parabolic approach. The values of the dominant size for each peak are given in Table 5.
Dominant size of mullite in kaolin DD (1200 °C).
The values of the dominant size evaluated from the relation of the second derivative of the different peaks of the principal phase of kaolin KT2 are given in Table 6.
Dominant size of the main phase of kaolin KT2.
According to Table 6, the average arithmetic size of the crystallites is 164 Å. From the Wialliamson–Hall diagram, we can estimate this average size (this time in the least squares sense) as 229 Å. The relative difference between the two values is less than 40%, which, in our opinion, is acceptable because the size of a monocrystallite can vary greatly depending on the direction of the diffracting planes and, moreover, there are several crystallites of varying size diffracting in the same direction.
Conclusion
By studying the contours of the lines, it seems very clear that the main phase of many sintered kaolin and mullite is devoid of internal microstrains. This is the case where the mixtures are whipped not only at low temperature (1200 °C) for 1 h, but also at a temperature of 1350 °C for very long periods (several weeks).
The Williamson–Hall plot was constructed and used to determine whether a sample is affected by deformation and to determine the size of the crystals. Since the specimen had no small folds, we applied the Scherrer relation. The Warren–Averbach method was also used in an undistorted condition.
This result is very important because it gives an element of explanation for a very important quality of mullite: its mechanical resistance during repeated use at high temperatures.
Footnotes
Acknowledgements
The authors thank all members of the material science department, University of Biskra, for their help and facilities.
Declaration of conflicting interests
The author(s) declared no potential conflicts of interest with respect to the research, authorship, and/or publication of this article.
Funding
The author(s) received no financial support for the research, authorship, and/or publication of this article.
