Jack S. Googasian,Sara E. Skrabalak
Practical Considerations for Simulating the Plasmonic Properties of Metal Nanoparticles
FDTD,DDA,plasmonic,nanoparticle,simulation,modeling,tomography
Simulating the plasmonic properties of colloidally derived metal nanoparticles with accuracy to their experimentally observed measurements is challenging due to the many structural and compositional parameters that influence their scattering and absorption properties. Correlation between single nanoparticle scattering measurements and simulated spectra emphasize these strong structural and compositional relationships, providing insight into the design of plasmonic nanoparticles. This Perspective builds from this history to highlight how the structural features of models used in simulation methods such as those based on the Finite-Difference Time-Domain (FDTD) method and Discrete Dipole Approximation (DDA) are of critical consideration for correlation with experiment and ultimately prediction of new nanoparticle properties. High-level characterizations such as electron tomography are discussed as ways to advance the accuracy of models used in such simulations, allowing the plasmonic properties of structurally complex nanoparticles to be better understood. However, we also note that the field is far from bringing experiment and simulation into agreement for plasmonic nanoparticles with complex compositions, reflecting analytical challenges that inhibit accurate model generation. Potential directions for addressing these challenges are also presented.
Practical Considerations for Simulating the Plasmonic Properties of Metal Nanoparticles

Simulating the plasmonic properties of colloidally derived metal nanoparticles with accuracy to their experimentally observed measurements is challenging due to the many structural and compositional parameters that influence their scattering and absorption properties. Correlation between single nanoparticle scattering measurements and simulated spectra emphasize these strong structural and compositional relationships, providing insight into the design of plasmonic nanoparticles. This Perspective builds from this history to highlight how the structural features of models used in simulation methods such as those based on the Finite-Difference Time-Domain (FDTD) method and Discrete Dipole Approximation (DDA) are of critical consideration for correlation with experiment and ultimately prediction of new nanoparticle properties. High-level characterizations such as electron tomography are discussed as ways to advance the accuracy of models used in such simulations, allowing the plasmonic properties of structurally complex nanoparticles to be better understood. However, we also note that the field is far from bringing experiment and simulation into agreement for plasmonic nanoparticles with complex compositions, reflecting analytical challenges that inhibit accurate model generation. Potential directions for addressing these challenges are also presented.

Nanoscale materials have properties that differ from their bulk counterparts, giving them applications in sensing, catalysis, medicine, and more. Optically active nanomaterials have been used in stained glass since ancient Rome, with the famous Lycurgus cup as a beautiful example where the apparent color of the cup given by the nanoparticles embedded within the glass depends on the location of the light source. Colloidal metal nanoparticles have attracted much interest for their widely tunable properties that are accessible from scalable, bottom-up approaches. Many colloidal syntheses produce metal nanoparticles of high symmetry; however, nanoparticles with structural and compositional complexity can also be obtained. For example, Au nanostars are useful plasmonic nanoparticles for biomedical applications but lack defined symmetry, and multimetallic nanoparticles are of increasing importance to applications such as catalysis. Moreover, colloidal syntheses of metal nanoparticles produce samples with inherent heterogeneity; i.e., slight structural variations from one nanoparticle to the next in a sample. As will be outlined herein, this complexity of nanoparticle samples often necessitates corroboration of their nanoscale properties with simulation. Nanoparticle features such as particle size, shape, composition, and architecture can be manipulated through colloidal syntheses. In turn, these features influence the localized surface plasmon resonances (LSPRs) of a metal nanoparticle. An LSPR is a collective oscillation of the conduction electrons of a nanoparticle induced by incident light. The energy of this resonance depends on the structure, composition, and surrounding environment of the nanoparticle, highlighting the attractiveness of colloidal syntheses but also the need to accurately capture these qualities in any simulation. Gustav Mie analytically solved Maxwell’s equations for the scattering properties of an electromagnetic plane wave from a homogeneous sphere in 1908. Systematic studies have subsequently shown excellent correlation between experimentally measured and analytically derived optical properties, highlighting how representative models of simple convex shapes were relatively straightforward to construct. However, the Mie solution is limited to simple structures such as spheres and cylinders. For more complex shapes, the finite-difference time-domain (FDTD) method—a numerical modeling tool for electrodynamics—is appropriate. Originally published by Yee in 1966, FDTD discretizes 3D space into compartments known as Yee cells across which Maxwell’s equations may be solved. After the simulation has been broken down into a 3D grid, the electric and magnetic field components are then sequentially and recursively updated as they propagate across the mesh in defined time steps. The interaction of the electronic and magnetic field components is then recorded for interpretation. The implementation of the FDTD method has historically been limited by the required computational resources but has become more generally available as those costs have decreased. The FDTD method is now widely used due to its simple implementation and generality to many systems. However, the meshing size of the Yee cells may result in discontinuities from staircasing curved surfaces to a Euclidean grid. This staircasing effect can be reduced by discretizing the 3D space to arbitrary precision, but the resources required then to process the simulation increase with mesh step size. In this Perspective, practical considerations for simulations of metal nanoparticles are discussed, such as this staircasing effect, with particular emphasis on bringing simulation and experiment into agreement for structurally and compositionally complex nanoparticles. We begin by highlighting select ground-breaking single-particle correlation studies that give insight into the importance of realistic nanoparticle models for FDTD and DDA simulations. Like FDTD, the DDA method is used to compute the scattering of light by particles of arbitrary size and shape, where in this case the particle is approximated as an array of dipoles that interact with the light and each other. We then discuss practical considerations when creating models for such simulations by showing the spectral changes that arise from subtle changes to a given model or simulation parameter. These considerations lay the foundation for the emerging interest in using tomographical reconstructions to create accurate models for simulation, which is proving quite useful to the field of chiroptical nanomaterials. We note that the examples presented here are not comprehensive but were selected to illustrate the main conceptual approaches toward high-quality model generation. We conclude with challenges to the field in accurately capturing nanomaterials with complex compositions for simulations.

Analyzing the relationship between particle structure and its optical properties is fundamentally difficult due to the inherent heterogeneity of NP samples, where commonly ensemble responses are measured and an assumption is made that the average structure observed accounts for the response. Yet, there is evidence that this assumption is a poor one. For example, in the seminal 1997 paper on single-molecule and single-nanoparticle surface-enhanced Raman scattering (SERS), Nie and Emory posited that the bulk of the Raman signal originated from less than 1% of the nanoparticles in the population, suggesting that the enhancement associated with this subpopulation was 106–107 higher than the ensemble averaged value. Single-particle studies can address the effects of heterogeneity from one NP to the next. For example, correlative studies can be undertaken in which the scattering properties of single NPs are measured alongside high-level structural characterization by electron microscopy. Such approaches allow for fine structural information to be correlated to the optical properties of specific NPs. In turn, structure–property trends can be established, and NPs with outlier properties can be identified. The importance of fine structural features to the plasmonic properties of metal NPs was particularly evident for a systematic study by Van Duyne and co-workers, where the optical responses of individual Au and Ag nanocubes were measured by dark-field scattering spectroscopy and these results were correlated to structural features from high-resolution transmission electron microscopy (TEM) images. To bring the experimentally obtained scattering spectra into agreement with FDTD simulations, the models needed to accurately capture the subtle structural features of the Au and Ag nanocubes. That is, the simplest models of platonic solids have perfectly defined edges and vertices, but as the TEM image of a Ag nanocube shows (Figure 1A), atomically sharp corners and edges are not observed. Instead, rounded corners are evident. The models developed by Van Duyne and co-workers for Au and Ag nanocubes captured the fine structural features observed by TEM by including face-to-face distance (c) and the radius of curvature (r) of isolated nanocubes (Figure 1B). However, this model was insufficient to accurately bring experiment and simulation in agreement. The dielectric environment surrounding the NPs had to be accurately captured, meaning that the substrate supporting the NPs and its distance from the NPs must be included in the model. Unfortunately, measuring the distance between the nanocubes and substrate is experimentally challenging. To gain insight into this feature, the optical response of Au nanocubes was systematically studied as a function of position and substrate identity by simulation, where a substrate of thickness hl, nl refractive index, and h distance away from the nanocube was added to the simulation (Figure 1B). As shown in Figure 1C, good agreement between simulation and experiment was achieved when the structural features of the Au nanocubes were included in conjunction with the substrate (either Si3N4 or Formvar), with the Au nanocubes positioned at the surface. Van Duyne and co-workers did note a few subtle discrepancies in terms of peak width and peak dips and readers are referred to the original reference for the excellent in-depth discussion. For our purposes, however, this systematic study shows quite well how important the fine features of models are to bring simulation and experiment into agreement. These same ideas hold when moving to metal nanoparticles of other shapes. For example, Van Duyne and co-workers also studied the plasmonic properties of single Ag triangular nanoprisms by dark-field optical microscopy and spectroscopy, bringing agreement between experiment and electrodynamic modeling based on the DDA method. This agreement was achieved through careful consideration of nanoprism height, edge length, and tip sharpness. This study was particularly impactful for showing how simulation can be applied to polydisperse nanoparticle samples, especially when coupled with single-particle approaches. That is, nanoparticle samples with large polydispersity will give rise to broad LSPRs as each nanoparticle in the sample has a slightly different spectrum that contributes to the ensemble response. This feature is evident in Figure 2A, where the black curve corresponds to the ensemble spectrum of the triangular Ag nanoprisms, while the scattering spectra from single Ag nanoprisms are much narrower as denoted by the colored curves, spanning a λmax of ∼600 to 700 nm depending on the edge length (Figure 2A). Notably, the scattering spectra of the single triangular nanoprisms fall within the envelope of the ensemble spectrum. Moreover, good agreement was achieved between simulation and experiment as shown in Figure 2B, where the model for simulations based on the DDA method was created from dimensions extracted from atomic force microscopy. Many studies report correlation between nanoparticle structures with their observed optical properties through simulation. These examples highlight how simulations can address complex samples with heterogeneity, providing relevant structure–property relationships. Significantly, the structural features of nanoscale systems must be accurately captured in models for sufficient correlation with experimental data.

Simulations of single-particle scattering spectra highlight the sensitivity of their λmax to slight structural perturbations. In this section, we raise awareness of the structural considerations for models used in FDTD simulations. Toward this end, we created a structural model of a Au nanocube similar to Van Duyne and co-workers, with the corners and edges are defined by a radius of curvature r and face-to-face distance c (Figure 1A). However, the refractive index of the medium was set at 1.333 (no substrate) to mimic colloidal measurements in water; the simulation mesh was set to 2 nm. The face-to-face distance and radius of curvature were systematically changed to investigate changes in λmax. Finally, the dielectric information reported by Johnson and Christy for bulk Au is used for the nanocube material. The dielectric information was measured by ellipsometry on thin films, but this information may also be measured by electron energy loss spectroscopy (EELS) and X-ray photoelectron spectroscopy (XPS). We note, however, that measured dielectric properties may vary from theory due to inter- and intraband transitions unique from free-electron contributions. We further acknowledge that there are deviations from the measured true dielectric function when approximated by a polynomial expression for FDTD; however, there is historically good agreement between experimental and numerical optical properties suggesting that the approximation sufficiently represents the dielectric behavior at least for many simple compositions. Figure 3A shows models of the Au nanocubes used in simulations, where c = 100 nm and r was varied from 0 to 40 nm in steps of 10 nm. The resulting scattering spectra are shown in Figure 3B, where a blue shift in λmax from ∼660 nm to ∼600 nm is evident from r = 0 to 40 nm. This blue shift is consistent with the work by Van Duyne and co-workers and highlights how models for FDTD simulations should incorporate realistic facet truncations to achieve experimentally realistic positions in λmax. We observed that the edge length, the distance between vertices on the cube, decreases with increasing rounding. Thus, an additional rounded model with edge length of 100 nm was created as well to compare rounding and face-to-face distance with edge length directly. We observe that although the rounding causes broadening of the scattering spectrum, the simulated λmax of both the unrounded and 100 nm edge length cubes are roughly equal (Figure 3A,B). Regardless, simulations of this nature can be quite useful for revealing trends in λmax with systematic changes in structure, aiding in plasmonic nanocrystal design. Still, for such trends to be successfully extracted, the simulation parameters must be carefully selected to avoid undesirable staircasing effects. Selecting the mesh size is of critical importance with the FDTD method, with structures with fine structural features (e.g., sharp corners) necessitating smaller mesh sizes than large structural features. While decreasing the mesh size generally decreases the accumulated error, the computational resources required increase. To illustrate the effect of mesh size, we performed simulations on a model nanocube system (c = 150 nm and r = 60 nm) with mesh sizes of 2, 4, and 8 nm step length (Figure 3C). The scattering spectra are shown in Figure 3D, with the effect of improper meshing particularly evident in the spectra obtained with 4 and 8 nm step lengths. With the larger step lengths, multiple peaks appear in the spectra that are inconsistent with the symmetry of the nanocube as well as experimental observations. These results highlight how those conducting such simulations must strike a reasonable balance between mesh size and computational time, where ensuring that the results make physical sense is essential. Finally, one must also be cognizant of the direction of light propagation when designing simulations and its effect on the results. The scattering properties of nanoparticles are orientation dependent, so the direction of light propagation must be considered for structures of low symmetry. That is, rotation of the model in 3D space is necessary to achieve an accurate average of the scattering properties. This effect is illustrated well by considering the simulation of intrinsically chiral Au nanocrystals such as those grown from tetrahedral seeds shown in Figure 3E; these nanocrystals consist of two flattened ellipsoids stacked on top of each other but offset by 90°. For brevity, these chiral nanocrystals will be referred to in this text as tetrahedral-seeded chiral Au nanocrystals. In Figure 3F are the simulated scattering spectra of the tetrahedral-seeded chiral Au nanocrystals, where each color trace corresponds with a g-factor spectrum from a different orientation. Notably, the simulated chiroptical behavior of this nanocrystal changes significantly with orientation although the average spectrum (black trace) agrees well with the experimentally measured response for these nanocrystals in solution. If a single orientation were taken instead, the chiral scattering properties would be inaccurate and, in some cases, a nearly complete inversion.

As the previous section illustrated, modeling and simulation of relatively simple structures such as cubes and spheres are both fairly straightforward with high-quality microscopy imaging. This ability is aided by the high symmetry of such structures. However, the previous section also showed the sensitivity of accurate simulation to fine structural features, which can become difficult to accurately capture in complex nanoparticles of lower symmetry compared to, for example, cubes or spheres. Robust methods for model generation are required. One method that has been leveraged to gain fine 3D structural information for complex nanocrystals is focused ion beam (FIB) milling for SEM imaging. FIB SEM enables acquisition of cross-sectional images that can reveal internal or sandwiched features of a nanoparticle, which can be directly incorporated or compiled together as a slice-wise structural model. For example, Lee et al. detailed the synthesis of intrinsically chiral Au nanocrystals through reduction of HAuCl4 onto Au seeds, where the chirality was induced by chiral capping molecules during growth. The structure of 432 helicoid III was probed by FIB helium ion milling coupled with SEM imaging as shown in Figure 4A. By gently ablating several nanocrystals, sequential cross-sectional SEM images were captured. From these images, the 3D nature of the chiral gaps within the nanocrystals could be more realistically understood as denoted by the red, green, and blue traces in Figures 4B–D and used to create a model for FDTD simulations (Figure 4C). Figure 4D shows the 3D interpolation of the curved surface based on the depth and curvature information obtained through FIB SEM analysis. The FDTD simulations corroborate the experimentally measured chiroptical activity of 432 helicoid III. Although the technique itself is quite sophisticated and the resulting model still idealized, FIB SEM imaging was essential to capturing the chiral gap features and bringing experiment and simulation into alignment. Unfortunately, FIB SEM is a destructive method, meaning that other analyses cannot be readily coupled with the high level of structural information captured. In contrast, electron tomography is a nondestructive quantitative method that can produce 3D reconstructions of nanoparticles through acquisition of scanning transmission electron microscopy (STEM) images at stepwise angles. Materials must be supported on a TEM grid and rotated through a wide range of angles for full 3D information. Thus, the technique is well-suited for nanoparticles that are simply dropcast on a grid but can be challenging for larger materials and those on substrates. The reconstruction with high shape accuracy is generated by importing these projections (i.e., images obtained at specific angles) into a processing algorithm, where we encourage the reader to read work by Jenkinson et al. for an in-depth discussion. An advantage of this method over FIB milling coupled with SEM imaging is the ability to generate models without introducing bias during data analysis through the destructive nature of FIB SEM itself. This advantage is particularly evident when the fine structural features of a nanocrystal cannot be inferred. Nondestructive analytical methods are particularly important to observe morphological and optical changes of single nanoparticles over time. Tomography is well-equipped to analyze polydisperse samples and complex shapes such as Au nanostars. Vanrompay et al. reported the synthesis and annealing of Au nanostars, measuring morphological changes by fast in situ electron tomography as shown in Figure 4E–H. Contrasting with traditional tomography which may take up to an hour in acquisition time, fast in situ tomography may be conducted in as few as 6 min. The small acquisition time of fast in situ tomography lowered the effect of the electron beam, which has been shown to affect the thermal stability of otherwise beam-insensitive materials. Therefore, the Au nanostars were able to be heated for a given time, cooled to room temperature, and reconstructed using electron tomography to observe gradual morphological changes due to thermal annealing. The tomographic reconstructions after annealing at 300 °C for 0, 30, and 1200 s are shown in Figure 4E–G, respectively. These reconstructions were then used as models to evaluate the effect of annealing time on the extinction cross-section of the Au nanostars by the boundary element method (BEM) of optical simulation. BEM differs from previously discussed methods as only the fields on a surface are discretized, allowing greater scalability when compared to volume element based methods. The simulated extinction spectra are shown in Figure 4H at 0 s (black trace), 30 s (red trace), and 1200 s (blue trace). As the nanostar shape became more rounded, the electric field enhancement of the nanostar greatly decreased and discrete LSPR bands coalesced into a broad extinction spectrum. Note that the magnitude of the calculated electric field enhancement is not consistent across all reconstructions in Figure 4E–G, decreasing with the thermal annealing time. This change is consistent with a decrease in tip sharpness and highlights the importance of accurately representing the radius of curvature for a nanostar’s tips in near-field simulations for analysis of SERS measurements. Electron tomography is therefore a powerful tool to correlate nuanced structural features to optical properties when coupled with simulation. A similar approach was used to shed insight into the tetrahedral-seeded chiral Au nanocrystals shown in Figure 3E; this model was created through tomographic reconstruction using high angle annular dark field (HAADF) scanning transmission electron microscopy (STEM). This level of characterization was essential as, from SEM images, we were unable to unequivocally define the intrinsic chirality that accounted for their chiroptical activity shown in Figure 4I (red trace), along with their extinction properties (black trace). However, this tomographically generated model was used in FDTD simulations, and by averaging over the various orientations as discussed in the previous section, good agreement between experiment and simulation was possible (Figure 4I versus Figure 4J). We particularly found this approach beneficial as there is also a need to untangle the chiroptical activity from near-field enhancement emanating from surface-bound chiral growth molecules. We anticipate that the application of tomographic reconstructions to create models for simulations will prove useful in many complex nanoparticle systems as shown similarly for Au nanostars, chiral Au nanorods, and chiral nanoparticle assemblies.

Single-particle correlations between experimentally obtained scattering spectra and simulation have highlighted the value of high-resolution models and have been used to extract many meaningful structure–property trends for monometallic nanoparticles. As structurally complex systems are now being handled with advanced characterization, the next challenge for the field will be accurately modeling nanoscale systems comprised of multiple elements. Multielement nanoparticles are selected to obtain synergistic effects and new properties altogether. However, the spatial arrangement of the different elements impacts the resultant properties. For example, Woessner et al. showed that the spatial distribution of Pd within/on branched Au nanoparticles with Oh symmetry impacts their sensitivity to changes in refractive index. In this study, Pd-tipped octopods were compared to core@shell Au@Pd octopods and alloy AuPd octopods. Interestingly, capturing the alloy composition accurately in FDTD simulations was challenging. Alloys are difficult to model due to a lack of data on the dielectric function of the material as well as insufficient analytical tools to adequately capture the distribution with atomic level precision. These challenges are captured well in a report by Schaak and co-workers in their study of plasmonic alloy Au1–xCux (x = 0 to 0.5) nanoparticles of quasi-spherical shape. The composition (Au:Cu ratio) of their colloidally prepared nanoparticles was extracted from the position of X-ray diffraction (XRD) reflections in accordance with Vegard’s law. Energy dispersive X-ray spectroscopy (EDS) was also used to determine the Au:Cu ratio. With this information as input, linear mixing was used to approximate the dielectric function of the alloy nanoparticles for simulation of the plasmon resonance frequencies by Mie theory. These results were compared to each sample’s spectrum determined by UV–visible spectroscopy. Discrepancies between the theoretically predicted plasmon resonance frequencies and those observed through experiment demonstrated that a mere linear combination of dielectric functions does not accurately reflect the composite dielectric function of these alloys. The origin of the discrepancies was not fully elucidated and has been reported in other systems as well. Work by our group in collaboration with the Dragnea group highlights additional challenges to accurately capturing nanoparticles with complex compositions. Specifically, a single-particle correlation study was undertaken in which stellated Au–Pd nanocrystals with Oh symmetry were analyzed by optical dark-field spectroscopy and SEM. The optical dark-field spectroscopy provided the λmax of each nanocrystal analyzed and SEM was used to obtain structural information (face diagonal and tip width of the stellated nanocrystals). Contour plots of LSPR λmax as a function of tip width and face diagonal as well as resonance width and face diagonal were obtained and compared to simulated data for stellated Au-only nanocrystals with the same structural features. Differences were evident from this comparison, revealing the impact of the complex Au–Pd distribution in these nanocrystals on their plasmonic properties. In fact, this analysis was consistent with the Au–Pd nanocrystals having Pd-rich tips, which was supported by single-particle elemental mapping by STEM-EDS. Still, we caution that EDS mapping measurements represent the 3D elemental composition as a 2D projection, which means that assumptions may be made with regards to elemental distribution when trying to use such data to create models for accurate modeling. Toward accurate spatial distributions of complex compositions, atom probe tomography (APT) is state-of-the-art. APT features resolution of subnanometer composition regardless of atomic number, although distinction between elements of similar atomic mass is difficult. In APT, a sharp sample is put under high voltage, causing a large electrostatic field at the tip. The atoms at the tip surface are then ionized and desorbed toward the ion detector. Position-sensitive detectors may be used to measure the image projections of the tip which change with composition. Field evaporation simulations calculate the trajectories of ions to determine the 2D projection on the detector. These simulations can account for aberrations such as changes in tip shape that impact the surface electric field often by combining high level simulation approaches. These simulations are then used to reconstruct the tip from the projection data. The accurate elemental reconstructions highlight the utility of APT to gain 3D compositional information; however, this method requires the samples to be sharply pointed and is restricted to a 150 nm field of view. Although nanoparticles may easily fit in this size constraint, reconstruction artifacts generated by differing evaporation fields of materials generally produce undesirable structural and compositional distortions leading to poorer spatial and morphological resolution when compared to other tomographic methods. Therefore, APT is often coupled with other methods to characterize the structure separately from the elemental composition. An example of a reconstruction from APT is shown in Figure 5A,B for complementary metal-oxide-semiconductor (CMOS) transistors, which are prominently used to manufacture integrated circuit technologies. 3D geometries of CMOS transistors display lower current leakage through tunneling effects and are therefore more desirable than their 2D planar counterparts. The example shown is in gate-all-around (Figure 5A) and trigate (Figure 5B) configurations, in which the silicon nanowire is surrounded by a gate material on all sides or three sides, respectively. The efficacy and electronic behavior of these structures are greatly influenced by the spatial elemental distribution, corresponding to dopant concentration and surface roughness. APT reconstructions performed on the materials accurately confirmed the elemental composition of both transistor types although the morphological distortions are present. Although APT can theoretically achieve atomic-scale elemental resolution, morphological distortions limit its utility for structure-dependent plasmonic simulation. Therefore, measurement of composition while also maintaining the structural information is advantageous. Toward this end, 3D tomograms with high spatial accuracy may be obtained by HAADF-STEM. HAADF-STEM imaging may then be coupled with either Z-contrast or EDS information to provide 3D elemental composition. This capability is growing with the advent of new detection capabilities for acquiring quantitative EDS data. For example, core@shell Au@Ag nanorods synthesized by Goris et al. shown in Figure 5C are an attractive system for sensing due to the ability to readily change their optical properties through the size and aspect ratio of the Au nanorod core. Accurately measuring the 3D structure and composition of these nanoparticles plays a vital role in understanding their optical properties. HAADF-STEM imaging was used to obtain a tomographical reconstruction. The images were processed using an advanced compressive sensing algorithm to lessen the presence of noise and scattering artifacts, enabling robust atomic resolution. Z-contrast was then used to distinguish Au and Ag atoms and obtain a 3D elemental reconstruction of the heterostructure as seen in Figure 5D. Spatial compositional resolution is particularly important for systems in which elements may readily alloy and diffuse. As another example of 3D elemental information, reconstructions of individual AuAg nanoparticles obtained via HAADF-STEM coupled with EDS are shown in Figure 5E,F. These AuAg nanoparticles were prepared by galvanic replacement between Ag nanoparticles and AuCl4–, giving rise to hollow structures. The compositional analysis in Figure 5E shows the surface confinement of Au, with a low level of replacement. However, greater addition of the Au precursor inverts the elemental distribution with Au concentrated to the interior and surface enrichment of Ag as shown in Figure 5F. The similar lattice parameters of Au and Ag facilitate interdiffusion, allowing for the inversion of surface composition. This transformation is difficult to measure by ordinary STEM EDS mapping because it compresses 3D compositional information into a 2D projection. However, by coupling HAADF-STEM tomography with EDS, the elemental composition information is embedded into the 3D structural tomograms, readily displaying the 3D arrangement of Au and Ag. Translating these approaches to simulation models is straightforward in systems with known dielectric constants and well-defined interfaces. In a simplified core@shell architecture assuming perfect interfacing between elements, the simulation mesh may be created by assigning proper mesh priority where the core is effectively cut out from the shell. FDTD simulations may then be used in tandem with these reconstructions to further approximate their plasmonic properties due to improved spatial elemental construction. In more complex systems such as the Au1–xCux alloys discussed previously, the dielectric information may not be directly inferred, leading to complications during analysis. Ellipsometry, XPS, and valence EELS have proven useful for determining the dielectric information in bulk monometallic systems, but direct analysis becomes difficult for 3D nanoparticles where spatial arrangement of elements may impact the measured dielectric function. Single-particle EELS measurements have been used to determine the dielectric functions of single-element systems, but may prove useful to reveal the complex dielectric functions of systems with complex elemental compositions. Therefore, these atomistic elemental characterizations highlight the challenge of accurately measuring and modeling composition for optical simulation.

Modeling and simulation of nanoscale systems allows for understanding of the underlying phenomena. As computational costs decrease, there is greater access to advanced simulation techniques. Such methods are now being widely implemented to numerically calculate the optical properties of metal nanostructures. This ability can provide insight into the origin of experimental observations, establish trends with regards to structure–property relationships, and give a pathway to addressing polydisperse samples. Yet, this ability relies heavily on the quality of the models used in the simulations, where methods like tomography are capturing fine details of nanomaterials. When moving to the use of simulation for the prediction of nanocrystals with new or enhanced properties, we encourage using insights from these studies to include such fine features from the start (e.g., corner and edge truncations, defects). Such details are already known to greatly impact near-field effects central to surface-enhanced Raman spectroscopy. That is, the electric field enhancement is greatly impacted by fine structural features as discussed previously for Au nanostars, as well as for many other nanostructures. Simulations must accurately represent these features (and an appropriate mesh size must be selected) because small volume changes can lead to large changes in the electric field enhancement, ultimately affecting the use of a given nanostructure as a SERS substrate. These constraints again highlight the importance of accurate modeling for systems with fine structural features. We are further encouraged by advances to couple tomographic information with composition as a pathway toward modeling nanoparticles of multiple material domains. The spatial arrangement of elements throughout nanostructures can greatly impact their properties and these advances in characterization that enable accurate compositional information in 3D has the potential to transform understanding of these materials when coupled with single-particle property measurements and even application of in situ stimuli. This advancement should shed light on when the assumption of linear mixing of dielectric functions is sufficient for compositions such as alloys, moving us further to the use of simulations for prediction as well.

FDTD simulations were performed using Lumerical FDTD Solutions software. The refractive index was set to 1.333. The dielectric functions for the Au models were fit to optical data collected by Johnson and Christy. Unless otherwise noted, the simulation mesh size was set to 2 nm. Nanocube models were constructed of 8 spheres of radius r at the corners, 12 cylinders of radius r, and height (c – 2r) on the edges, and 3 orthogonal rectangular prisms of lengths c, (c – 2r), and (c – 2r) for the cube faces. Rotation of the models was performed in 15 deg increments from 0 to 75 deg in both the x- and y-axes.