A reproducible number-based sizing method for pigment-grade titanium dioxide

Summary A strong demand for reliable characterization methods of particulate materials is triggered by the prospect of forthcoming national and international regulations concerning the classification of nanomaterials. Scientific efforts towards standardized number-based sizing methods have so far been concentrated on model systems, such as spherical gold or silica nanoparticles. However, for industrial particulate materials, which are typically targets of regulatory efforts, characterisation is in most cases complicated by irregular particle shapes, broad size distributions and a strong tendency to agglomeration. Reliable sizing methods that overcome these obstacles, and are practical for industrial use, are still lacking. By using the example of titanium dioxide, this paper shows that both necessities are well met by the sophisticated counting algorithm presented here, which is based on the imaging of polished sections of embedded particles and subsequent automated image analysis. The data presented demonstrate that the typical difficulties of sizing processes are overcome by the proposed method of sample preparation and image analysis. In other words, a robust, reproducible and statistically reliable method is presented, which leads to a number-based size distribution of pigment-grade titanium dioxide, for example, and therefore allows reliable classification of this material according to forthcoming regulations.


Introduction
Titanium dioxide is among the ten most abundant materials on the Earth [1]. In the form of a fine powder, it is used as white pigment in many application systems such as paints, plastics, paper and building materials. It is also used in cosmetics, foods and pharmaceuticals. Its superior properties as white pigment are based on its high refractive index, leading to maximum whiteness and opacity, if its particle size distribution is opti-mized for best scattering efficiency according to Mie's theory [2,3]. The most common industrial processing routes are the sulfate process and the chloride process. In the sulfate process, for example, ilmenite ore is dissolved in sulfuric acid, iron and titanium are separated by controlled precipitation, and colouring transition metals are removed in a bleaching process prior to calcination. In the chloride process, rutile in the form of sand or The images (a) through (e) demonstrate that the projected area of the pigments varies greatly, depending on viewing direction. Moreover, they illustrate the difficulty with particle detection in projection images in the event of overlapping particles. slag, for example, is treated with gaseous chlorine to form titanium tetrachloride. The titanium tetrachloride is purified by distillation and then transformed into solid TiO 2 in a combustion process. After both processes, the resulting powders are de-agglomerated by standard milling procedures and then subjected to a finishing process, which is generally followed by a final de-agglomeration step. In the ideal case, the size (distribution) of the final products is optimised for the optical pigment properties, which are described by using Mie's theory. Mie's theory states that the optimum particle size for refracting light is just about half the wavelength it is meant to interact with. Speaking of visible light, a mean wavelength of 500 nm is reasonable to assume. Therefore, a mean pigment size of 250 nm is close to what the titanium dioxide makers are producing.
The calcination temperature, as the final step of the powder synthesis of the sulfate process, gives the choice of producing either the low-temperature crystallographic phase, anatase, or the high-temperature crystallographic phase, rutile. The chloride process leads to formation of rutile phase particles due to the high temperatures of the combustion process. Both phases have a number of advantages and disadvantages, which lead to their typical applications. Preferred fields of application for rutile pigments are coatings, paints, plastics and building materials, whereas anatase pigments are mainly used in cosmetics, pharmaceuticals or food. One of the most important properties of titanium dioxide is its UV absorption, which protects human skin against sunburn and skin cancer. Optically transparent TiO 2 is the most important ingredient of any commercially available sun cream. The most common way of achieving a transparent, highly UV-absorbent sun cream is to use intentionally manufactured nano-sized titanium dioxide.
An urgent demand for reliable methods for the characterization of particulate materials is triggered by the prospect of forthcoming national and international regulations concerning the classification of nanomaterials [3][4][5][6]. Scientific efforts to establish standardized methods for determining number-based size distributions have so far been focused on model systems, such as spherical gold particles, spherical latex particles or spherical silica particles [7][8][9][10]. The typical obstacles in the characterization of industrial materials, such as irregular particle shapes [11], a broad size distribution [12] and strong aggregation and agglomeration effects [13,14], have not been addressed successfully. But as industrial materials are the goal of all regulatory efforts, robust, reproducible and reliable methods for the determination of the number-based size distributions for a wide range of industrial materials are urgently needed [3].
This study focuses on two commercially available pigments, KRONOS K2360 and KRONOS K1171. The former is a pigment with a rutile structure for use in coatings and paints, while the latter is a food-grade pigment with an anatase structure. We decided to develop a procedure that requires a minimum of subjective, user-based decisions. The major task to tackle on the way is to prepare a sample suitable for automated particle detection. Since automated detection routines are primarily based on grey-value thresholding, overlapping particles are a serious problem, which is illustrated in Figure 1. As an example, Figure 1 shows five different projections of one agglomerate found in K2360. The number of obviously visible primary particles varies for the different projections, as well as the size of the projected area, which is generally evaluated. These TEM tomogram-based data illustrate the fundamental problem associated with number-based particle size measurements: A true size distribution of the pigment will only be obtained if the full, three-dimensional shape of the particles is evaluated, e.g., as the equivalent sphere diameter or a minimum 3D Feret diameter. All known evaluation methods based on the evaluation of projections or slices contain systematic deviations [15][16][17][18]. As for every particle sizing method, besides the physical principles applied, the decision as to how to measure the particle size and its distribution is somewhat arbitrary and generally restricted by practicality aspects.
As the first selective decision, the present study uses polished sections in which overlapping particles are ruled out. As a consequence, the measurement used here does not correspond to a projected area of the particles, but to a section through the particles. The correlation between a measured size distribution and the "real" size distribution is the subject of so-called stereologic correction: Based on the assumption of a common, known particle shape, stereologic corrections may be used to estimate the "real" size distribution, either from the measured projection sizes or from the measured section sizes [15,19,20]. It is a fact that there is nothing like a common or known particle shape for titanium dioxide pigments. The presence of arbitrary shapes in industrially produced titanium dioxide makes stereologic correction impractical, both for the evaluation of projections and for the evaluation of sections. Consequently, the choice of the measuring technique is somewhat arbitrary. The technique presented in this study was chosen because of its high degree of automation and reproducibility. The effect of lacking stereologic correction, if significant, leads to a result which is slightly smaller than the "real" particle size. Therefore, the uncorrected section-size distributions presented in this publication give lower estimates of the real particle size distributions of the analysed pigments. For the intended purpose, i.e., classification of the material according to the recommendation of the European Commission [4], the method gives a conservative estimate of the particle size distribution.

Results and Discussion
In order to establish the proposed method, especially for the sizing of pigment-grade titanium dioxide, the reproducibility of the method was primarily tested. Each pigment was prepared several times according to the procedure described in the experimental section. The first preparation run included four samples, labelled Ra01 through Ra04 for the rutile pigment and Aa01 trough Aa04 for the anatase pigment. A second preparation run was done several days later. It included the rutile samples Rb05 through Rb10 and the anatase samples Ab05 through Ab08. The first measurement of each sample is labelled M1; the second measurement was performed by a different operator and is labelled M2. The measured results in terms of minimum Feret diameter and equivalent circle diameter (ECD) are summarised in Tables 1 to 4. The first column gives the number of detected primary particles, followed by mean and standard deviation. Meaningful quantiles (d 10 , d 16 , d 25 , d 50 and d 84 ) are given in the next columns. The data of the different samples are treated without any assumption of a model for the distribution function. Consequently, no standard error based on the standard deviation and the number of particles counted is given. Instead, the standard error is calculated as the standard deviation of the repeated measurements, as done for the mean and all evaluated quantiles. The standard error is used to calculate the 95% confidence interval.
The high reproducibility of the measurement, which leads to a standard error below 2 nm for all evaluated measurements, justifies use of the non-parametric rank sum test of Siegel and Tukey [21] to test whether the measured distribution can be attributed to the same population. The comparison of all measured samples proves that this is the case, both for KRONOS K2360 and for KRONOS K1171. No significant differences between the standard errors for the ECD and minimum Feret diameter measurements are obtained for the rutile pigment. For the anatase pigment, lower standard    errors are found for the minimum Feret diameter. This is attributed to a less regular shape as a consequence of the production process.
The data prove that the method presented for primary particle sizing of pigmentary titanium dioxide is highly robust and reproducible. This is the consequence of a rigorous procedure, which primarily targets reliability: a representative, macroscopic amount of sample to start with, preparation that eliminates the possibility of overlapping particles, defined standard operating procedures for sample preparation, measurement and evaluation, and a minimum influence of operator-based uncertainties. The result is an obtained standard error of less than 2 nm and a relative coefficient of variation below 1.6% for all measured quantities.
However, the method has several limitations, which also need to be addressed and discussed here. First of all, the procedure was developed for pigment-grade titanium dioxide prepared as a cross-section. An adjustment of the method may possibly be needed for its application to other particulate materials with other size distributions. The presented results are, accurately speaking, the distribution of ECD and minimum Feret diameter of the pigment sections. Stereologic correction has not been performed, and is not feasible for this application [15,20]. Errors introduced by stereologic effects may gain importance if the size of the measured objects is large compared to the penetration depth of the electron beam. The size distributions measured from polished sections can be expected to differ from those of projected areas. However, the high reproducibility of the proposed method will certainly allow comparison with any other sizing method, as long as a reference or standard for this purpose is present.

Conclusion
A highly reproducible, statistically tested method for the sizing of pigment-grade rutile and anatase is established. The standard error of the method is shown to be below 2 nm for all measured quantities, the relative coefficient of variation being below 1.6%. The presence of a systematic bias due to the lack of stereologic correction cannot be verified for the time being, since a certified standard for titanium dioxide is still lacking. The reproducibility of the method is based on the use of a representative, macroscopic amount of sample material, a high degree of automation, the elimination of detection errors due to overlapping particles and a transparent filtering procedure for detected particles. As a consequence, the method presented is well suited to classifying pigment-grade titanium dioxide according to the recommendation of the EU Commission of 18 October 2011 on the definition of particulate nanomaterial, or for other forthcoming regulations in the future.

Experimental
The experimental procedures described below are in full agreement with the practical guide for particle size characterization published by NIST [22].

Sample preparation
A well-prepared sample is the basis for reliable evaluation of the size distribution of a pigment. The limited quantity of samples used in electron microscopy emphasises this point.
Starting with approximately 10 g pigment, which is assumed to be a representative sample, the following steps are taken: 1. The "cone and quartering" method is repeatedly used to divide the sample in half until a quantity of 2 g is left. 2. The 2 g of pigment are thoroughly mixed with 4 g hotmounting resin and 22.5 g ZrO 2 milling beads. Mixing is done for 4 min at a moderate frequency of 20 Hz, using a Retsch Mixer Mill MM400. The milling beads are removed after the mixing process. 3. The mixture is then hot-mounted by using a Struers Citopress under standard conditions of 180 °C and 3 bar pressure.

4.
A five-step polishing process, with active oxide polishing as the final step, leads to a smooth, scratch-free section through the embedded pigments. 5. As the final step before SEM investigation, an 8 nm Au/ Pd conductive layer is deposited on top of the sample surface.
The sample preparation procedure summarised above leads to a sample with a high pigment density and allows imaging of several thousand pigments within a set of approx. 50 images. The "cone and quartering" method ensures that the embedded sample is representative. A macroscopic amount is taken for the mixing and embedding procedure. Since the whole representative amount is embedded, effects of flushing are ruled out and the low viscosity of the resin ensures that sedimentation is not a problem. The section finally analysed can therefore be reasonably assumed to be representative. The high area density of the pigment in all images ensures that constant conditions for the automated post-processing and detection procedure are assured.

Measurement, pigment detection and size analysis
A working distance of 7 mm and an acceleration voltage of 5 kV with an Everhart-Thornley detector are chosen as the standard imaging conditions for SEM imaging. The microscope used is a Leo 1530VP SEM. A magnification of 20.000× (5.7 nm/pix) is chosen for the rutile sample, which displays a larger primary particle size, and a magnification of 30.000× (3.8 nm/pix) for the anatase sample, which has a significantly smaller primary particle size. In both cases, this corresponds to approx. 20 pixels for the equivalent circle diameter (ECD) of the 10% quantile d 10 [23]. A silicon wafer with etched structures is used for calibration of the instrument. The size of the structures is certified by the PTB under Serial No. IMS-HR 08 3641-01 490 and the calibration mark 44049 11 PTB. An automated protocol (macro) for image post-processing and particle detection was developed. The procedure was realised with the commercial software package "analysis" from Olympus-SIS. All filters used are standard filters implemented in the software. The detailed settings are defined in the standard operating procedure and are listed below. This allows maximum transparency of the procedure and maximum comparability and reproducibility of the results, e.g., when dealing with regulatory authorities or their contracting measuring partners.
The full procedure includes unification of the grey values, masking of the unified image, particle detection and filtering of the results, as shown in the images of Figure 2. The first step ensures comparability, not only within the set of images of one sample, but also between different samples or even between samples measured on different instruments with varying detector settings or noise levels. The comparability is sufficient, as long as the penetration depth of the electrons, which is determined by the acceleration voltage, remains constant. In the second step, the images are binarised by using automated greyvalue threshold determination [22]. A morphological particle separation filter, the watershed transform, is applied to the binary image. The resulting binary image with the separated particles is used as a mask for the unified image. The final detection of the primary particles is done on the masked unified image, which allows for the detection of the particles, including their grey values. The grey values of the particles are important for the subsequent filtering of the results obtained. It is used in order to remove particles that are located below the polished surface, but give rise to an increased intensity compared to the background. In general, these falsely detected particles show the lowest detected intensities and, therefore, also the lowest standard deviations. A combination of both values has proven adequate for grey-value filtering. Additional filtering is based on the pigment shape. As for the grey values, combined filtering of two shape characteristics has proven advantages, namely the so-called convexity and the shape factor: 1. Unifying the grey values as the basis for final pigment detection (Figure 2b) • automated contrast and brightness adjustment • median filter for noise suppression • rank filter for noise suppression with a size of 5 pixels and a rank of 50% • multiplicative shading correction for a size of 255 × 255 pixels and a cut-off level of 10% 2. Preparing and applying a mask image, followed by particle detection (Figure 2c,d) • automated grey-value thresholding and binarisation • repeated removal of pixels connecting only to three or less neighbouring pixels, using the connectivity filter • applying the watershed transform in order to separate aggregated or adjacent particles (Figure 2c   • removing all particles with a "shape-factor" below 0.86 as non-primary particles • removing all particles with a convexity of less than 0.90 as non-primary particles Figure 3 visualises the particle detection characteristics that are intrinsic to the analysis of a polished section. The section through the particles can in principle be located close to the centre of mass, in which case the section gives the maximum area along the viewing direction for a convex particle (plane 3 in Figure 3a). But it can also be directed above the particle (plane 1 and 2 in Figure 3a) or below the centre of mass with respect to the viewing direction (plane 4 in Figure 3a), which leads to a smaller detected area. The cases of planes 3 and 4 in Figure 3 are essentially indistinguishable -they provide sharp particle edges. The cases of planes 1 and 2 give rise to a rim of lower intensities at the particle edges. Figure 3b shows how the different cases look in the SEM image, the numbers given corresponding to the section planes in Figure 3a. Only the particles with a coloured boundary are the ones finally detected after applying all filters. It can be seen that grey-value filtering effectively filters those particles that are located fully below section plane 1. Particles sectioned above the centre of mass, labelled with a 2, are only partially filtered out. The boundaries are detected reasonably well for the counted particles. The detection of particles sectioned close to the centre of mass or below the centre of mass is precise, both cases are labelled with a 3 in Figure 3b, since they are indistinguishable in the SEM image.
Uncounted particles labelled with a 3 are filtered out as non-primary particles due to their shape.