3.1. Classification of Au Nanoplates and Multiply Twinned Au Nanoparticles
The Au NPs possess exceptional optical properties, which are caused by the localized surface plasmon resonance (LSPR) effect [
45,
46]. The LSPR effect is a result of the collective resonant oscillations of conduction electrons, which are excited by the electromagnetic field of the incident light [
7,
47]. An important consequence of the LSPR effect is the appearance of optical extinction at a specific frequency that depends mainly on the size [
48] and shape [
49] of the Au NPs, and that is affected by the kind and density of the crystal structure defects [
12]. The Au NPs can grow as small faceted single crystals having truncated octahedral (TOh) shape [
21], as multiply twinned particles (MTPs) [
50] with either icosahedral (Ih) or decahedral (Dh) shape, or as nanoplates (NPLs) [
17]. The Au NPs contain typically planar defects that are mainly stacking faults (SFs) and twins (TWs). The shape of Au NPs depends usually on the conditions of the Au NP synthesis [
21]. However, as the transition between the individual types of the Au NPs is smooth, Au NPs with different properties can occur within the same sample [
25]. Therefore, an experimental technique for statistically reliable classification of the Au NPs is required, in particular if the structural characteristics of the NPs are to be correlated with their properties. This section illustrates, how the NP shape can be correlated with the kind of the crystal structure defects and how the combination of HRTEM and HAADF-STEM improves the statistical quality of the Au NP classification.
HRTEM images of individual Au NPs (
Figure 1) reveal a coexistence of NPLs and MTPs. While NPLs appear as truncated triangles (upper panels in
Figure 1a and
Figure 1b) with almost perfect crystal structure, MTPs look like disks with a high density of planar defects (upper panels in
Figure 1c and
Figure 1d). According to FFT/HRTEM (lower panels in
Figure 1a and
Figure 1b), almost all Au NPLs possess a
orientation along the direction of the primary electron beam. The truncated triangles are terminated by the lattice planes
. The diffraction spots highlighted by dotted circles that resemble the reflections
are produced by the stacking faults that are located on the lattice planes
, which are perpendicular to the direction of the primary electron beam. Such SFs are not visible directly in the HRTEM images, but they produce truncation rods
, which intersect the Ewald sphere at the reciprocal space vector having the size of
where
is the interplanar spacing of the lattice planes
,
Å the lattice parameter of Au and
or
the angle between the faulted lattice planes
and the crystallographically equivalent lattice planes
. As
for both
angles, the
positions of the truncation rods
,
,
,
,
and
are almost equal to one third of the
positions of the reciprocal lattice points
,
,
,
,
and
,
i.e.,
In the HRTEM images of multiply twinned particles (
Figure 1c and
Figure 1d), the stacking faults are visible directly, because the normal directions to the faulted lattice planes are not perpendicular to the diffraction vector. The projected form of the particles appears almost circular for both MTPs classes,
i.e., for Dh shaped (
Figure 1c) and Ih shaped (
Figure 1d) MTPs. Still, these particle shapes can be distinguished by means of FFT/HRTEM. The FFT/HRTEM of the Dh MTPs shows a pattern with a `five-fold symmetry’ consisting of ten 111 diffraction spots (marked by solid circles in
Figure 1c) and ten 200 diffraction spots (marked by solid boxes in
Figure 1c) that have approximately equidistant azimuthal positions. The FFT/HRTEM patterns of the Ih MTPs contain additional diffraction spots (marked by dashed circles in
Figure 1d), which are caused by the Moiré pattern effect [
51].
As illustrated above, several kinds of Au NPs,
e.g., `flat triangular’ NPLs with different degree of truncation and `spherical’ MTPs having the Dh or Ih shapes, can be differentiated using the HRTEM imaging. Using FFT/HRTEM, the NP edges can be assigned to the crystallographic directions and the NP shape can be correlated with the kind and orientation of the planar defects [
18,
40]. However, these techniques operating on the atomic scale are not suitable to reveal a statistically relevant information about the individual NP fractions. Therefore, for statistical reasons, HRTEM was combined with a `low-magnification’ HAADF-STEM imaging, which is indeed not able to visualize the crystal structure defects directly,
i.e., on the atomic scale, but it can recognize different NPs according to their projected shape. Furthermore, as the HAADF-STEM signal stems from the Rutherford scattering of the primary beam electrons on the atomic nuclei within the sample, the HAADF-STEM intensity measured on Au NPs depends mainly on the NP thickness [
52]. As the Au NPs contain only a single element, the dependence of the HAADF-STEM intensity on the atomic number does not play any role. Consequently, the HAADF-STEM signal from Au NPs can be used to determine their 3D form.
The information about the NP thickness and especially the information about the thickness variation within individual NPs can help to classify the Au NPs into the `flat triangular’ NPLs and `spherical’ MTPs. For a quantitative classification of the NPs, the relative variance of the HAADF-STEM intensity within individual NPs (
Figure 2a) was used:
The mean HAADF-STEM intensity,
and its variance,
were calculated from the HAADF-STEM intensity values (
) that were measured within individual NPs.
N in equations (
4) and (
5) denotes the number of pixels within the respective NP.
In order to improve the reliability of the classification procedure,
from Equation (
3) was correlated with the area-equivalent diameter of the respective NP that was calculated from the projected area
A:
The area itself was determined from the number of pixels, which were assigned to the respective NP. In
Figure 2b, the correlation between
and
is visualized using two bivariate density estimators [
53]
that are based on the Gaussian kernel functions
K.
and
are the characteristics of the
NP, and
and
are the bandwidths of the kernel function that were determined according to Scott’s rule [
54].
This classification procedure allows NPLs to be distinguished from MTPs, and the respective NP fraction to be determined. It can be seen from
Figure 2a and
Figure 2b that NPLs show a smaller variation of the HAADF-STEM intensity and that are typically larger than MTPs. The transition between NPLs and MTPs (red line in
Figure 2b) was determined using the indicator function
For MTPs,
is equal to unity, while for NPLs,
. The number of MTPs was obtained from the `weighted’ integration of the function
,
The MTP fraction was calculated as
where
is the total number of NPs. The fraction of NPLs is equal to
. The statistical analysis carried out with ∼ 2,700 Au NPs revealed that the sample under study contains approx. 90% of MTPs and about 10% of NPLs.
In conjunction with the results of HRTEM and FFT/HRTEM, the result of HAADF-STEM provides an important insight into the kinetics of the growth process of Au NPs, because the kind and the density of planar defects and consequently the morphology of the Au NPs are controlled by the reaction rate during the synthesis [
20,
25]. At a sufficiently high reaction rate, defect-free NPs with truncated octahedral (TOh) shape develop [
21,
55]. Reduction of the reaction rate leads to the stabilization of multiply twinned particles (MTPs) [
21,
50,
56]. A further reduction of the reaction rate promotes formation of NPs with plate-like morphology that contain a high density of planar defects preferentially on a single system of the
lattice planes [
17,
57,
58]. The absence of TOh NPs and the presence of ∼ 90% MTPs and ∼ 10% NPLs in the Au NPs under study indicate that the reaction conditions were moderate. However, the reduction rate was not sufficiently high to produce exclusively MTPs.
Although the statistical HAADF-STEM analysis of the Au NPs that is supported by the FFT/HRTEM investigation of few selected NPs reveals a valuable information about the kind and distribution of the planar defects, and although this technique provides an important insight into the growth conditions, further details, for instance about the stacking fault density, can only be obtained using the analytical methods operating on the atomic scale.
3.2. Classification of Differently Faceted Au Nanorods
Similarly to Au NPs, the Au nanorods (Au NRs) show also unique plasmonic properties [
59]. However, in contrast to Au NPs that possess typically one LSPR band, Au NRs produce two LSPR bands, which correspond to the plasmonic oscillations along their short and long axes [
46]. The optical properties of Au NRs depend mainly on their size and aspect ratio. For applications in catalysis, the crystallographic orientation of the Au NRs facets plays a crucial role. Zhang
et al. [
15] have shown that Au NRs terminated by high-index facets possess a higher catalytic activity than Au NRs with low-index facets.
The example from
Section 3.1 illustrated the capability of the statistical classification of Au NPs by HAADF-STEM supported by HRTEM and FFT/HRTEM. However, in that example the information obtained from HRTEM was not fully linked with the information obtained from HAADF-STEM, as the exact projected NP shape was not considered when analyzing the HAADF-STEM images. The NPLs with triangular (less truncated) and hexagonal (heavily truncated) form or the MTPs with decahedral and icosahedral shape were not classified separately (
Figure 2). Thus, the densities of planar defects were not determined statistically. The example presented in this section illustrates, how Au NRs can be classified into four categories, when their 3D shape is used.
The in-depth characterization of the Au NRs on the atomic scale was carried out using HRTEM and FFT/HRTEM on several tens of nanorods. The HRTEM images taken on reclined Au NRs (
Figure 3a) confirmed that they are single-crystalline, possess the
fcc structure, and are defect-free and elongated along one of the crystallographic directions
, as it was already reported by Zhang
et al. [
15] and Ye
et al. [
24]. The caps of the Au NRs are usually formed by the high-index facets
. The angle between these facets is 143°, as it is visible in the HRTEM image from
Figure 3a, where the facets
and
are highlighted in red. The cross sections of the Au NRs (
Figure 3b) are typically terminated by the high-index facets
[
24], but in some cases also the low-index facets of the
and
types were found. This result was confirmed by the presence of differently oriented reclined NRs (s.
Figure 3d).
According to
Figure 3c, the angles between the high-index facets
are either 143.1°,
e.g., between
and
, or 126.9° like for
and
. The angles between the corresponding crystallographic directions
, and
and
are 36.9° and 53.1°, respectively. The high-index facets are believed to develop in the final stage of a seed-mediated Au NR synthesis, as they smooth the sides of the NRs by removing the edges between the previous side facets
and
[
28]. In general, a high number of crystallographically equivalent lattice planes distributed along a specific zone axis (
in
Figure 3c) smooths the kinks between the neighboring facets. From the crystallographic point of view, the number of crystallographically equivalent lattice planes increases with increasing number of non-equal
h,
k and
l values within the Miller index
. Thus, a combination of the high-index and low-index facets also facilitates the smoothing of the NR surface. From the atomistic point of view, however, the high-angle facets produce steps on the surface of the NRs.
For subsequent statistical analysis using HAADF-STEM, the Au NRs were divided into four categories of `faceted cylinders’ capped by differently oriented lattice planes. The cross-sections of the cylinders are either
complete or
incomplete. The NRs with
complete cross-sections are terminated by the high-index facets
only. The NRs with
incomplete cross-section possess one low-index facet, typically of the type
(
Figure 3b). The caps of the NRs are either
symmetric or
asymmetric. The NRs with
symmetric caps are terminated on both ends by the same facet type,
i.e.,
and
or
and
. The NRs with
asymmetric caps are terminated by facets, which are mutually rotated by 90° at the respective end,
i.e.,
and
or
and
.
Since the majority of NRs lies horizontally on the TEM grid, the NR symmetry is typically visible directly from the 2D HAADF-STEM projection (
Figure 4a). In the semi-automatic segmentation and classification routine [
40], the ratio of the NR diameters measured at the top and at the bottom of the NRs (
Figure 4b),
is employed as the parameter quantifying the degree of the NRs symmetry. For NRs with
symmetric caps,
. For NRs with
asymmetric caps,
(
Figure 5). The completeness of the NR cross-sections is quantified by
where
is the width of the plateau in the HAADF-STEM intensity profile measured across the NR (
Figure 4c) and
D the width of the whole NR. For
complete NRs,
(cf.
Figure 5), where
is the angle between the directions
and
,
i.e., 18.4°. For
incomplete NRs,
is less than 0.42. In the example from
Figure 4c showing the presence of a single
facet,
would converge to zero.
In
Figure 5, the two-dimensional distribution density of the parameters
and
is depicted in terms of the bivariate kernel density estimator,
, which was determined in analogy to Equation (
7). Individual maxima of
,
i.e.,
,
,
and
, were fitted by 2D Gaussian functions and assigned to NRs with
incomplete cross sections and
symmetric caps, to NRs with
complete cross sections and
symmetric caps, to NRs with
incomplete cross sections and
asymmetric caps, and to NRs with
complete cross sections and
asymmetric caps, respectively.
The fitted functions
,
,
and
were used to identify the boundaries between the NR categories (red lines in
Figure 5) and to determine the amounts of NRs having the respective facet configuration. The boundaries between the NR categories obey the relationships
The amount of NRs having the particular facet configuration follows from the integration of the bivariate kernel density estimator
over the respective
region. For example, the amount of NRs with complete cross sections and symmetric caps is equal to
where
is the indicator function for NRs with the particular facet configuration. The amounts of NRs with other configurations of the facets are calculated analogously.
The statistical evaluation of the HAADF-STEM measurements summarized in
Figure 5 revealed that about 33% of the Au NRs have
complete cross sections and
symmetric caps. According to HRTEM, these NRs are fully terminated by the high-index facets
. Approximately 24% of the Au NRs still possess
complete cross sections, but have
asymmetric caps. In total, the cross sections of
% Au NRs were
incomplete. Approximately 28% of them have
symmetric and about 15%
asymmetric caps. In this context, it should be noted that the completeness of the cross sections quantified using Equation(
12) is primarily an indicator of the presence or absence of parallel opposite facets oriented perpendicular to the direction of the primary electron beam (
Figure 4c). Nevertheless, assuming an equal growth rate in the crystallographically equivalent directions
, the second maxima of
appearing at
(
Figure 5) can only be explained by the presence of other facets than
that are perpendicular to the primary beam. Still, the non-parallel opposite facets like in
Figure 3b, whose presence would lead to extremely low
values, are rare.
The stabilization of specific facets in noble metal NPs is generally attributed to the presence of certain surface capping agents [
31] and structure-directing ions [
23,
60,
61]. The formation of high-index facets
on the Au NRs observed in this study was facilitated by the Ag ions [
23,
60] stemming from the addition of AgNO
to the seed-mediated solution. Stabilization of the facets
is often attributed to the presence of bromide ions [
29,
62]. In our case, the bromide ions stem from the seeds that were synthesized using the CTAB protocol [
24,
40].
In literature, Au NRs with alternating high-index
and low-index
facets bordering their cross sections are reported more frequently than the Au NRs terminated exclusively by the high-index facets [
63,
64,
65]. Also in our sample,
% of the Au NRs were terminated by mixed high-index and low-index facets. An analogous result was also obtained for the NR caps, which were predominantly
symmetric (
Figure 5). This result is in a good agreement with previous investigations [
23,
24,
66,
67] showing that, under uninhibited growth conditions, the
symmetric caps occur more frequently than the
asymmetric caps. The
symmetric or
asymmetric arrangement of the NR caps can, however, depend on the aspect ratio of the NRs. Ye
et al. [
24] reported that the amount of NRs with
asymmetric caps increases, when the aspect ratio of the NRs increases.
3.3. Hierarchical Architecture of Multi-Core Iron Oxide Nanoparticles
In both previous examples, individual nanoparticles and nanorods were sufficiently separated from each other. Thus, they could be identified and quantified quite straightforwardly from the HAADF-STEM images using a shape-based segmentation routine [
40]. However, this routine fails, when it is applied to overlying objects,
e.g., to multi-core iron oxide nanoparticles (IONPs). An example of a multi-core IONP is shown in
Figure 6a. The multi-core nature of the IONP becomes clearly visible from the FFTs of the HRTEM images (
Figure 6b and
Figure 6c), which disclose different orientations of the cores A and B along the direction of the primary electron beam,
i.e.,
and
.
The attachment of differently oriented cores is facilitated by the agreement in the distances of the parallel lattice planes of both counterparts at their interface. In the mutual orientation relationship of the cores from
Figure 6a, the following lattice planes are almost parallel and their interplanar distances equal:
,
and
, s.
Figure 6b and
Figure 6c. Slight differences in the local orientations within the cores are visible from the rigid rotation field (
Figure 6d) that was determined using the geometric phase analysis (GPA) [
38,
39]. For GPA, the reflections
,
,
and
were used. While the GPA done on the core A did not reveal noticeable orientation variations, the GPA carried out on the core B indicated that this core consists possibly of two parts, which are separated by a low-angle boundary (dotted line in
Figure 6a). These characteristics of the multi-core IONPs obtained from HRTEM and FFT/HRTEM, in particular the detailed information about the crystallographically oriented attachment of individual cores, are extremely helpful for understanding the formation of multi-core IONPs and their magnetic properties.
The multi-core IONPs show superparamagnetic behavior with a high saturation magnetization and a good biocompatibility, which makes them to favored materials for applications in biomedicine, especially in magnetic hyperthermia. The magnetic properties of the multi-core IONPs depend mainly on their size and chemical composition [
8], but they are also affected by the disorder of the magnetic moments [
68,
69,
70] and/or by the magnetic coupling between the neighboring cores [
71]. Bulk iron ferrites are strongly ferrimagnetic, and exhibit spontaneous magnetic moment and hysteresis [
72]. When the size of the IONPs is reduced, they become superparamagnetic [
8] with nearly zero remanent magnetization and coercivity, but retain a high magnetic susceptibility. A further decrease of the IONP size leads to a decrease of the saturation magnetization of the IONPs, which is attributed to the magnetic disorder in the surface layer [
73,
74]. Another reason for the reduction of the saturation magnetization of IONPs is the presence of vacancies on the cation sites in the crystal structure of Fe
O
, which also contributes to the disorder of the magnetic moments and thus to the decrease of the saturation magnetization [
75,
76]. While the saturation magnetization of vacancy-free magnetite (Fe
O
, SG
) is 92.8 Am
/kg [
77], the saturation magnetization of maghemite (Fe
O
, SG
or
[
78]) is only 74.3 Am
/kg [
77].
In order to be able to distinguish individual parts of the multi-core IONPs from each other, the shape-based segmentation routine [
40] was extended (
Figure 7). In the first step (
Figure 7a and
Figure 7b), the size of individual multi-core IONPs and the corresponding size distribution were determined from the projected area as described in
Section 3.1. In the second step (
Figure 7c-g), the individual IONPs were separated into the individual cores, which were quantified analogously. This two-step segmentation and quantification procedure revealed the size distribution functions, which are depicted in
Figure 8a. The distribution density was calculated according to
where
is the number of particles within the respective “size category”,
i.e., having their diameter between
and
.
D in the argument of the distribution function is the mean value of
and
.
is the total number of particles. The distribution functions from
Figure 8a were compared to the distribution function of the magnetic particle size (
Figure 8b), which was determined using the Langevin fit of the magnetization curve measured by alternating gradient magnetometry on an ensemble of IONPs (inset in
Figure 8b) [
71].
In analogy with the model of the multi-core IONPs suggested on the basis of the HRTEM and HAADF-STEM results (
Figure 6 and
Figure 7), the presence of three categories of magnetic particles was assumed, corresponding to (i) the fragments of the cores (nanocrystals within the cores that are frequently terminated by low-angle boundaries), (ii) the individual cores and (iii) the multi-core IONPs. Hence, the distribution of the magnetic particle size (
Figure 8b) was composed of three log-normal functions:
In Equation (
17),
is the fraction of the magnetic particles having the mean size
;
is the width of the respective log-normal function. This distribution function was used to calculate the magnetization curve [
82,
83]:
where
is the Langevin function, which argument,
depends on the saturation magnetization (
), on the size of the magnetic particles (
), on the temperature (
T), and on the strength of the external magnetic field (
H).
is the permeability of vacuum,
the Boltzmann constant. The calculated magnetization curve was fitted to the measured one (inset in
Figure 8b). The refineable parameters were the saturation magnetization (
), the fraction of the magnetic particles (
), the medians of the magnetic particle sizes (
), and the widths of the log-normal functions (
). The starting values of
and
were taken from the HRTEM and HAADF-STEM results (
Figure 8a).
For multi-core IONPs, the size distribution functions from
Figure 8a and
Figure 8b match quite good. Both techniques,
i.e., the HAADF-STEM imaging and the determination of the size of the superparamagnetic nanoparticles from the magnetization curve [
82], revealed the same mean size of the multi-core IONPs and the same width of the size distribution, (
) nm. This means that in multi-core IONPs, the magnetic moments in individual cores are highly coupled. Thus, such multi-core IONPs resemble large magnetic particles. The coupling of the magnetic moments in neighboring cores is apparently facilitated by their specific orientation relationships, as it was shown exemplarily in
Figure 6.
Mean size of the individual cores obtained from the magnetization curve, (
) nm, also agrees with the size of the cores determined by HAADF-STEM. In this context, it should be noted that the `magnetic’ size of the cores is typically smaller than their geometric size, because the magnetization of small particles is reduced by a disordered spin layer at their surface [
74]. Still, the size distribution function obtained for the individual cores from the HAADF-STEM analysis (
Figure 8a) is much broader (and thus less intense) than the size distribution function determined from the magnetization curve (
Figure 8b). As the analysis of the HAADF-STEM images classified the IONPs according to their 2D morphology (projected area) only and as it was not complemented by the information about the crystallographic orientations of adjacent cores, it cannot distinguish between the cores with coupled and uncoupled magnetic moments. Consequently, the size distribution function obtained from the HAADF-STEM measurements includes also multi-core IONPs with elliptical projected shape that consist of magnetically uncoupled cores. Finally, the fitting of the magnetization curve confirmed presence of the core fragments. However, their fraction was very low, because the majority of the core fragments possessed distinct mutual orientation relationships that facilitated a coordinated arrangement of the magnetic moments.