Determination of the volume-specific surface area by using transmission electron tomography for characterization and definition of nanomaterials
© Van Doren et al; licensee BioMed Central Ltd. 2011
Received: 31 January 2011
Accepted: 11 May 2011
Published: 11 May 2011
Transmission electron microscopy (TEM) remains an important technique to investigate the size, shape and surface characteristics of particles at the nanometer scale. Resulting micrographs are two dimensional projections of objects and their interpretation can be difficult. Recently, electron tomography (ET) is increasingly used to reveal the morphology of nanomaterials (NM) in 3D. In this study, we examined the feasibility to visualize and measure silica and gold NM in suspension using conventional bright field electron tomography.
The general morphology of gold and silica NM was visualized in 3D by conventional TEM in bright field mode. In orthoslices of the examined NM the surface features of a NM could be seen and measured without interference of higher or lower lying structures inherent to conventional TEM. Segmentation by isosurface rendering allowed visualizing the 3D information of an electron tomographic reconstruction in greater detail than digital slicing. From the 3D reconstructions, the surface area and the volume of the examined NM could be estimated directly and the volume-specific surface area (VSSA) was calculated. The mean VSSA of all examined NM was significantly larger than the threshold of 60 m2/cm3.
The high correlation between the measured values of area and volume gold nanoparticles with a known spherical morphology and the areas and volumes calculated from the equivalent circle diameter (ECD) of projected nanoparticles (NP) indicates that the values measured from electron tomographic reconstructions are valid for these gold particles.
The characterization and definition of the examined gold and silica NM can benefit from application of conventional bright field electron tomography: the NM can be visualized in 3D, while surface features and the VSSA can be measured.
The number based size distribution of a material and the features of its surface are predominant criteria to classify it as a NM [1, 2]. TEM remains an important technique to measure the size and surface topography of materials at the nanometer level. Because the resulting micrographs are two-dimensional projections of the studied objects, their interpretation can be difficult, particularly when the particles are complex, agglomerated or lack symmetry. In such cases, fine ultrastructural details are blurred due to superposition of projected features. In addition, parameters like the surface area and volume of NM are not accessible by conventional TEM, while the approach to measure the thickness of NM along the projection direction by analyzing focal series in TEM assumes a relatively simple structure . Recently, as data acquisition, alignment and reconstruction software evolves to be more user-friendly; ET is increasingly used to reveal the morphology and to evaluate the three-dimensional characteristics of NP and nanoparticle ensembles [4, 5].
To include also aggregates and agglomerates of primary particles and complex multi-component particles with external dimensions larger than the arbitrarily specified upper size limit of 100 nm, the VSSA is proposed as a complementary qualifier to distinguish a nanostructured material from a non-nanostructured material . The European Commission  proposes to define a material as a NM when it has a specific surface area by volume greater than 60 m2/cm3, excluding materials consisting of particles with a size lower than 1 nm. The VSSA of a material is generally calculated from its bulk density and its mass specific surface area. The latter is usually determined by gas absorption methodology called the BET-method  that allows surface area or porosity measurements as small as 1 nm. From a 3D reconstruction of a NM, its surface area and its volume can, in principle, be estimated directly, such that its VSSA can be calculated, even on a per particle basis. Advanced electron tomography methods were applied advantageously and successfully to characterize NM at a high resolution [4, 5, 7, 8]. Most TEM-facilities do however not dispose of the required expensive equipment and lack the specialized expertise. Conventional electron tomography, where reconstructions are generated from a tilt series recorded in bright-field mode, using a single tilt axes with a tilt range up to ± 70°, becomes however a well-established technique. In this study, we examined the feasibility of three-dimensional visualization of silica and branched gold NM in suspension using conventional bright field (BF) ET. We examined whether such materials can be defined as a NM based on the measurement of their VSSA from its electron tomographic reconstruction. To evaluate the influence of missing wedge artifacts on the reconstruction and on the precision of the estimation of the surface area and volume of such NM, ET analyses of spherical colloidal gold nanoparticles were used as a control.
Suspensions of spherical and branched gold NP were obtained from IMEC (Heverlee, Belgium). Aggregated silica nanomaterials NM-200 and NM-203 are supplied by the European Commission-JRC (Ispra, Italy) as representative reference NM. They are used as well at the OECD Working Party for Manufactured Nanomaterials programme as principal materials and international harmonization standards. The NM were brought on pioloform- and carbon-coated 400 mesh copper grids (Agar Scientific, Essex, England) that were pre-treated with 1% Alcian blue (Fluka, Buchs, Switzerland) to increase hydrophilicity, as described by Mast and Demeestere . Gold NP were used undiluted. NM-200 and NM-203 were suspended in water containing 2% Fetal calf serum (PAA Laboratories GmbH, Pasching, Austria) at a concentration of 0.1 mg/ml and sonicated using a Vibracell™ 75041 sonicator (750 W, 20 kHz, Fisher Bioblock Scientific, Aalst, Belgium) with a 3 mm probe at 40% amplitude (10 W). A total energy of approximately 6200 J was added to the samples.
To obtain a maximal field of view, grids were mounted in a tomography holder (FEI, Eindhoven, The Netherlands) such that the squares were oriented diagonally with respect to the axis of the holder. Only objects in the centre of a grid square were analyzed using a Tecnai Spirit TEM (FEI) with a BioTWIN lens configuration and a LaB6-filament operating at an acceleration voltage of 120 kV.
Series of micrographs (tilt-series) were recorded semi-automatically assisted by the Xplore 3D tomography-module of the microscope control software (FEI) over a tilt range of at least 65°, or highest angle possible, at intervals of 1 degree. Shift and focus changes were corrected at every interval. Electron micrographs were acquired with a 4*4 K Eagle CCD-camera (FEI) at magnifications of 26,500 to 49,000 times and corresponding pixel sizes of 0.49 to 0.22 nm. The tilt series were aligned using the Inspect 3D software, version 2.5 (FEI) by iterative rounds of cross correlation until the alignment shifts were approaching to zero. Because of their higher signal to noise ratio, reconstructions using 10 to 20 iterations of the Simultaneous Iterative Reconstruction Technique (SIRT) algorithm were superior over reconstructions based on weighted back projection (WBP) and on the Algebraic Reconstruction Technique (ART) algorithm (not shown).
For visualization in 3D, the Amira software, version 4.1.2 (Mercury Computer Systems, France) was used. Isosurface rendering was used to compute a triangular approximation of the interfaces between the segmented sections. The segmentation was obtained based on a single threshold. This was chosen such that the obtained surface optimally matches the boundaries of the reconstructed orthogonal digital slices (orthoslices) of the NM in the xy-plane, where resolution is highest. The resulting surface was visualized using pseudo-coloring. To reduce missing wedge artifacts, so-called streaks, the surface was smoothed using a 2 × 2 × 2 averaging of voxels (downsampling). Using the 'Create Surface' function of Amira, a surface was derived from the isosurface, which allowed measurement of the surface area of the reconstructed 3D objects and of their enclosed volume.
Two-dimensional parameters of the reconstructed NP were measured from the TEM micrographs taken at 0° using the AnalySIS Solution of the iTEM software (Olympus, Münster, Germany). Briefly, contrast and brightness of the micrographs were optimized, the involved particles were enclosed in a frame (region of interest) and thresholds were set to separate particles from the background based on their electron density and size. The surface area and volume of individual spherical particles were approximated by the formulas to calculate the surface area (4πr2) and the volume (4/3πr3) of a perfect sphere, where r is replaced by the measured ECD of the projected particle divided by two. The sphericity, describing the 'roundness' of a particle by using central moments, was used to assess the hypothesis that the particle is a sphere in reality.
To measure the strength of correlation between the calculated VSSA and the measured VSSA, the nonparametric Spearman rank order correlation test was calculated using the SigmaPlot software, version 11.0 (Cosinus Computing B.V., Drunen, The Netherlands). To test the hypothesis that the mean VSSA obtained from ET reconstructions equals the threshold of 60 m2/cm3, the one-sample t-test (Sigmaplot) was used.
ET of spherical gold nanoparticles
Mean volume specific surface area of different nanomaterials based on electron tomographic reconstructions
Type of nanomaterial
Volume-specific surface area (m2/cm3)a
316 ± 7
177 ± 29
Precipitated Silica (NM-200)
342 ± 36
Pyrogenic Silica (NM-203)
219 ± 23
The reconstructed gold particles showed no obvious elongation along their z-axis and image analysis of the transmission electron micrographs of the individual particles taken at a tilt angle of 0° resulted in a mean sphericity of 0.86. Hence, it was concluded that these gold particles are almost spherical and that their surface area and volume can be closely approximated by the formulas to calculate the surface area and the volume of a perfect sphere. Figure 1C and 1D show the correlations between the calculated and measured volume and surface area, respectively, for ten ET reconstructions consisting of one to 11 gold NP. Both for the volume and the surface area, the Spearman correlation coefficient was 0.98.
ET of branched gold nanoparticles
ET analyses of silica NM
By electron tomographic reconstruction based on conventional BF TEM, the general morphology of gold and silica NM was visualized in 3D. In orthoslices of the examined NM in the xy-plane, as presented in Figure 2B, the surface can readily be distinguished from the background and from missing wedge artifacts, like streaks. In such orthoslices, the surface features of a NM can be seen and measured without interference of higher or lower lying structures inherent to conventional TEM.
Segmentation by isosurface rendering allows accessing the 3D information of an ET reconstruction in greater detail than digital slicing. Such 3D visualization and measurement of the surface features of NM can contribute to bring the second condition of the definition of a nanomaterial proposed by the European Commission  in practice: structures in one or more dimensions in the size range of 1-100 nm can be shown.
From the 3D reconstructions, the surface area and the volume of the examined NM could be estimated directly and the VSSA was calculated. The mean VSSA of all examined NM was significantly larger than the threshold of 60 m2/cm3 such that these materials can be classified as NP according to the third condition of this definition. As opposed to the BET-method , ET is not limited to powders and/or dry solid materials: it can be applied to a large variety of NM samples, including suspensions of complex particles, provided that the material can be suitably coated on an EM-grid.
To optimally characterize the morphology of a NM by ET reconstruction, it is required that (i) the projection requirement is met ; (ii) missing wedge artifacts are minimal and (iii) isosurface rendering optimally fits the NM surface.
Our results indicate that, in principle, the characterization and definition of NM can benefit from application of conventional BF ET. In the scope of putting this technique in practice for the characterization and definition of gold and silica NM, following approach is suggested to reconcile the limitations of conventional BF ET with the above-described conditions.
(i) The projection requirement states that for an image intensity to be usable for ET reconstruction, it has to be a monotonic function of a projected physical quantity . The examined silica NM were shown to be amorphous and weak scattering such that their mass thickness is the dominant contrast mechanism. The BF images of the tilt series are thus essentially projections on which tomographic reconstructions can be based . In the branched gold particles, and to a small extent in the spherical gold particles, diffraction contributed to image formation and the projection requirement is not fulfilled for the entire tilt series. Certainly, a combination of scanning transmission electron microscopy (STEM) and high angle annular dark field imaging (HAADF) which is insensitive to Bragg diffraction will be preferable over bright field imaging to visualize these NM at high resolution [11–13]. Because diffraction increases the background of the reconstruction and reduces its resolution, BF ET has been suggested to be of only limited value to analyze crystalline nanostructures [11, 12]. However, Ahrenkiel et al.  argumented that conventional BF ET still can provide useful information on the structure of particles with relatively small crystallite size if suitable acquisition conditions are chosen. Hence, the examined material was not embedded to assure positive contrast originating from the specimen at all orientations, while diffraction contrast was minimized while preserving some mass-thickness contrast by using a large objective aperture.
(ii) An important physical limitation of electron tomography arises from the fact that finite specimen thickness and tilt geometry within an electron microscope column prevent the collection of projection images spanning a complete angular range (± 90° tilt series). This results in a "missing wedge" of information in reciprocal space and results in anisotropic resolution in the resulting reconstruction . Only for specific samples that were properly shaped using a focused ion beam and mounted in a special holder, these subsampling effects could be avoided . The missing wedge can be reduced by collecting a so-called dual-axis tilt series in two mutually orthogonal directions [7, 16–18] such that a missing pyramid is obtained. Because above-described approaches require high technicity and are unpractical for routine analyses, the missing wedge artifacts of the reconstructions of the silica and gold NM were only reduced using a small tilt increment and maximizing the tilt range of the electron tomogram. Because gold and silica NM are hardly sensitive to radiation damage, extensive data collection using a high frequency of imaging can be applied. Because our software and hardware limit the amount of data that can be aligned and reconstructed in a timely manner, 4*4 K micrographs were taken with one degree intervals and datasets were reduced to contain only few particles. New developments of soft- and hardware and GPU-based ET implementation  promise faster data processing allowing a combination of smaller intervals and the simultaneous analysis of several hundreds of particles under similar imaging conditions. Moreover, improvement of the quality of reconstructions seems possible by using newly developed reconstruction algorithms like the discrete algebraic reconstruction technique (DART) which suffers less from missing wedge artifacts than SIRT . In our hands, the correction of the tilt axis during alignment appeared very important to reduce reconstruction artifacts like streaks. When the axes were not corrected accurately, the streaks were included in the particle volume resulting in an elongation of the particles which was at least as strong as the missing wedge dependent elongation described in .
(iii) In the examined NM, the value of the threshold was selected such that the obtained surface optimally matches the boundaries of the reconstructed orthogonal digital slices (orthoslices) of the NM in the xy-plane, where resolution is highest. This threshold value was in general close to the minimal value between both peaks of the bimodal curve (Figure 1A) of the histogram representing the number of voxels in function of their grey value. In the future, computational techniques that determine the optimal grey value for thresholding can allow more efficient segmentation, eliminating the subjectivity associated with manual segmentation .
To evaluate the influence of missing wedge artifacts on the reconstruction and, in particular, on the validity of the quantitative results obtained from the reconstructions, the ET based analyses of colloidal gold nanoparticles with a known spherical morphology were evaluated as a control. The high correlation (r = 0.98) between the measured values of area and volume and the areas and volumes calculated from the ECD of projected NP indicates that the values measured from electron tomographic reconstructions are valid for these gold particles. On this basis, it is assumed that the surface and volume measures of the branched gold and silica NM, which lack symmetry, can be relied on also. However, it has to be stressed that the absolute numbers presented in Table 1 should be interpreted with caution since they are based on a very low number of observations. This table illustrates the possibilities of the method in principle but it remains unsure whether the selected particles are representative for the entire examined samples. The latter requires, at least, the ET analysis of larger numbers of randomly selected particles. The evaluation of the VSSA values measured using methods like ET, BET or small-angle X-ray scattering (SAXS), and their corresponding uncertainties and limitations requires a dedicated study. For NM-200, the VSSA estimated by ET (342 ± 36 m2/cm3) is higher than the value obtained by SAXS (270 ± 17 m2/cm3, personal communication Camille Guiot, French Atomic Energy and Alternative Energies Commission, France), but lower than the values obtained by BET (435 m2/cm3, personal communication Rosica Petrova, Institute of Mineralogy and Crystallography - Bulgarian Academy of Sciences, Bulgaria). For NM-203, the VSSA estimated by ET (219 ± 23 m2/cm3) is lower than both the VSSA obtained by SAXS (367 ± 30 m2/cm3, personal communication Camille Guiot) and BET (469 m2/cm3, personal communication Rosica Petrova).
The EU definition expects a resolution of at least 1 nm . It is doubtful whether this resolution can be obtained routinely by BF TEM . For an exact description of the physical characteristics of a NM, atomic resolution within a reconstruction can be important. It is however not an absolute prerequisite to apply the definition if it is taken into account that lack of resolution results in an underestimation of the VSSA because, relatively, nanometer-sized surface features contribute more to an increase of the particle surface than to an increase of its volume.
As a proof of principle, it was shown that application of conventional BF ET allows 3D visualization of the examined gold and silica NM and allows measuring their surface features and VSSA. This approach can hence contribute to bring the second and third condition of the definition of a nanomaterial proposed by the European Commission in practice . Recent technical developments promise for the near future the possibility to analyze large numbers of particles  representative for the sample, a better reconstruction  and less influence of missing wedge artifacts [7, 16–18] such that the characterization of nanomaterials by transmission electron tomography can become more precise and less time-consuming.
This study was funded by the Federal Public Service of Health, Food Chain Safety and Environment (contract RT 09/6223 NANO-TEM) and by the project NANOGENOTOX which has received funding from the European Union, in the framework of the Health Programme. We thank Dr. Christoph Klein (JRC, Ispra, Italy) for providing us the silica NM and Dr. Ir. Bieke Van de Broek and Dr. Ir. Tim Stakenborg (IMEC, Heverlee, Belgium) for providing us the colloidal and branched gold NM. We thank Dr. Ir. Camille Guiot and Dr. Rosica Petrova for providing us with their SAXS and BET results. This publication reflects only the authors' views and the Executive Agency for Health and Consumers is not liable for any use that may be made of the information contained therein.
- SCENIHR: Scientific basis for the definition of the term "Nanomaterial" SCENIHR opinion for public consultation approved by written procedure on 6 July 2010. [http://ec.europa.eu/health/scientific_committees/emerging/docs/scenihr_o_030.pdf].
- European Commission recommendation on the definition of the term "Nanomaterial"(2010). [http://ec.europa.eu/environment/consultations/pdf/recommendation_nano.pdf].
- Alloyeau D, Ricolleau C, Oikawa T, Langlois C, Le Bouar Y, Loiseau A: Comparing electron tomography and HRTEM slicing methods as tools to measure the thickness of nanoparticles. Ultramicroscopy. 2009, 109 (7): 788-796. 10.1016/j.ultramic.2009.02.002.View ArticleGoogle Scholar
- Midgley PA, Ward EP, Hungría AB, Thomas JM: Nanotomography in the chemical, biological and materials sciences. Chem soc Rev. 2007, 36 (9): 1477-1494. 10.1039/b701569k.View ArticleGoogle Scholar
- Friedrich H, de Jongh PE, Verkleij AJ, de Jong KP: Electron tomography for heterogeneous catalysts and related nanostructured materials. Chem Rev. 2009, 109 (5): 1613-1629. 10.1021/cr800434t.View ArticleGoogle Scholar
- Brunauer S, Emmett PH, Teller E: Adsorption of gases in multimolecular layers. J Am Chem Soc. 1938, 60 (2): 309-319. 10.1021/ja01269a023.View ArticleGoogle Scholar
- Lengyel JS, Milne JL, Subramaniam S: Electron tomography in nanoparticle imaging and analysis. Nanomedicine. 2008, 3 (1): 125-131. 10.2217/17435822.214.171.124.View ArticleGoogle Scholar
- Batenburg KJ, Bals S, Sijbers J, Kübel C, Midgley PA, Hernandez JC, Kaiser U, Encina ER, Coronado EA, Van Tendeloo G: 3D imaging of nanomaterials by discrete tomography. Ultramicroscopy. 2009, 109 (6): 730-740. 10.1016/j.ultramic.2009.01.009.View ArticleGoogle Scholar
- Mast J, Demeestere L: Electron tomography of negatively stained complex viruses: application in their diagnosis. Diagnostic Pathol. 2009, 4: 5-10.1186/1746-1596-4-5.View ArticleGoogle Scholar
- Hassellöv M, Readman JW, Ranville JF, Tiede K: Nanoparticle analysis and characterization methodologies in environmental risk assessment of engineered nanoparticles. Ecotoxicology. 2008, 17 (5): 344-361. 10.1007/s10646-008-0225-x.View ArticleGoogle Scholar
- Midgley PA, Weyland M: 3D electron microscopy in the physical sciences: the development of Z-contrast and EFTEM tomography. Ultramicroscopy. 2003, 96 (3-4): 413-431. 10.1016/S0304-3991(03)00105-0.View ArticleGoogle Scholar
- Arslan I, Yates TJ, Browning ND, Midgley PA: Embedded nanostructures revealed in three dimensions. Science. 2005, 309 (5744): 2195-2198. 10.1126/science.1116745.View ArticleGoogle Scholar
- Sueda S, Yoshida K, Tanaka N: Quantification of metallic nanoparticle morphology on TiO2 using HAADF-STEM tomography. Ultramicroscopy 2010.Google Scholar
- Ahrenkiel SP, Yu PR, Murphy JE, Nedeljković JM, Donohoe BS: Nanoparticle shape and configuration analysis by transmission electron tomography. J Microsc. 2008, 230 (Pt 3): 382-387.View ArticleGoogle Scholar
- Ke X, Bals S, Cott D, Hantschel T, Bender H, Van Tendeloo G: Three-dimensional analysis of carbon nanotube networks in interconnects by electron tomography without missing wedge artifacts. Microsc Microanal. 2010, 16 (2): 210-217. 10.1017/S1431927609991371.View ArticleGoogle Scholar
- Penczek P, Marko M, Buttle K, Frank J: Double-tilt electron tomography. Ultramicroscopy. 60 (3): 393-410.
- Mastronarde DN: Dual-axis tomography: an approach with alignment methods that preserve resolution. J Struct Biol. 1997, 120 (3): 343-352. 10.1006/jsbi.1997.3919.View ArticleGoogle Scholar
- Tong J, Midgley P: A novel dual-axis reconstruction algorithm for electron tomography. Journal of Physics: Conference Series. 2006, 26: 33-36.Google Scholar
- Xu W, Xu F, Jones M, Keszthelyi B, Sedat J, Agard D, Mueller K: High-performance iterative electron tomography reconstruction with long-object compensation using graphics processing units (GPUs). J Struct Biol. 2010, 171 (2): 142-153. 10.1016/j.jsb.2010.03.018.View ArticleGoogle Scholar
- Batenburg KJ, Sijbers J: Optimal threshold selection for tomogram segmentation by projection distance minimization. IEEE Trans Med Imaging. 2009, 28 (5): 676-686.View ArticleGoogle Scholar
- Thomas JM, Midgley PA, Yates TJ, Barnard JS, Raja R, Arslan I, Weyland M: The chemical application of high-resolution electron tomography: bright field or dark field?. Angew Chem Int Ed. 2004, 43 (48): 6745-6747. 10.1002/anie.200461453.View ArticleGoogle Scholar
This article is published under license to BioMed Central Ltd. This is an Open Access article distributed under the terms of the Creative Commons Attribution License (http://creativecommons.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.