ACS Publications. Most Trusted. Most Cited. Most Read
My Activity
CONTENT TYPES

Figure 1Loading Img
RETURN TO ISSUEPREVCondensed Matter, In...Condensed Matter, Interfaces, and MaterialsNEXT

Simulating Plasmon Resonances of Gold Nanoparticles with Bipyramidal Shapes by Boundary Element Methods

  • Jacopo Marcheselli
    Jacopo Marcheselli
    SISSA—Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
  • Denis Chateau
    Denis Chateau
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Frederic Lerouge
    Frederic Lerouge
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Patrice Baldeck
    Patrice Baldeck
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Chantal Andraud
    Chantal Andraud
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Stephane Parola
    Stephane Parola
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Stefano Baroni*
    Stefano Baroni
    SISSA—Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
    *Email: [email protected]
  • Stefano Corni*
    Stefano Corni
    Dipartimento di Scienze Chimiche, Università di Padova, 35131 Padova, Italy
    Istituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, Italy
    *Email: [email protected]
  • Marco Garavelli*
    Marco Garavelli
    Dipartimento di Chimica Industriale “Toso Montanari”, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
    *Email: [email protected]
  • , and 
  • Ivan Rivalta*
    Ivan Rivalta
    Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
    Dipartimento di Chimica Industriale “Toso Montanari”, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, Italy
    *Email: [email protected]
    More by Ivan Rivalta
Cite this: J. Chem. Theory Comput. 2020, 16, 6, 3807–3815
Publication Date (Web):May 7, 2020
https://doi.org/10.1021/acs.jctc.0c00269

Copyright © 2022 American Chemical Society. This publication is licensed under these Terms of Use.

  • Open Access

Article Views

2307

Altmetric

-

Citations

LEARN ABOUT THESE METRICS
PDF (4 MB)
Supporting Info (1)»

Abstract

Computational modeling and accurate simulations of localized surface plasmon resonance (LSPR) absorption properties are reported for gold nanobipyramids (GNBs), a class of metal nanoparticle that features highly tunable, geometry-dependent optical properties. GNB bicone models with spherical tips performed best in reproducing experimental LSPR spectra while the comparison with other geometrical models provided a fundamental understanding of base shapes and tip effects on the optical properties of GNBs. Our results demonstrated the importance of averaging all geometrical parameters determined from transmission electron microscopy images to build representative models of GNBs. By assessing the performances of LSPR absorption spectra simulations based on a quasi-static approximation, we provided an applicability range of this approach as a function of the nanoparticle size, paving the way to the theoretical study of the coupling between molecular electron densities and metal nanoparticles in GNB-based nanohybrid systems, with potential applications in the design of nanomaterials for bioimaging, optics and photocatalysis.

Introduction

ARTICLE SECTIONS
Jump To

Metal nanoparticles, noble metals in particular, feature a wide spectral tuning of the localized surface plasmon resonance (LSPR), (1,2) as a function of their size, shape, and dielectric environment. (3−6) Gold nanoparticles (NPs) with bipyramidal shape (gold nanobipyramids, GNBs, see Figure 1) are particularly well-known for their high tunability and optical sensitivity, (7) featuring a narrower plasmon band with respect to other nonspherical nanoparticles, like rods or dog-bone (8−10) and allowing fine-tuning of the plasmon resonance energy. (11−17) Moreover, the tips of gold nanobipyramids enhance the effects of the local electric field generated by the metal nanoparticle, which make these nano-objects particularly sensitive to local changes in the dielectric environment, suitable even for detection of single molecules (18,19) and various technological applications. In particular, given the possibility of producing GNBs with variable size and aspect ratio (AR) at high yields and in almost monodispersed samples, (20,21) these nano-objects represent ideal inorganic partners for the assembly of hybrid organic–inorganic nanoparticles with peculiar optical properties, with potential impact in technological devices for bioimaging, photocatalysis, optics, solar cells, and biotechnology. (22−24)

Figure 1

Figure 1. (a) The GNBs’ models used in this work, including, from left to right the bipyramidal with pentagonal base and spherical tip (BPS), the biconical with spherical tip (BCS) and the biconical with flat tip (BCF) shapes. (b) A schematic representation of the basic geometric parameters of GNBs, with base length (w), height (h), and 2D projected area (in yellow) being extracted from TEM images, from which the nontruncated, ideal height (h′) can be derived.

The interaction between the optically induced collective excitations of free charges confined in the metal nanoparticle surface and the excitons of organic chromophores is one of the main features of metal–organic nanohybrids. (25) Moreover, upon resonant excitation, LSPR can exceed the diffraction limit and collect light into a subwavelength region, greatly enhancing the electric field acting over optically active molecules adjacent (adsorbed or physisorbed) to the metal surface. (26−28) In general, the LSPR operates quenching or enhancing the photophysical properties of the nearby molecular systems, resulting for instance in modified absorbance and/or fluorescence, or in preserving the photochemical stability of the chromophores by opening nonreactive excited state deactivation processes. (14,29) To allow and control such intriguing synergistic optical properties and to design new nanomaterials for versatile technological applications, computational modeling could be extremely beneficial, as it could provide rationalization of experimental evidence and eventually in silico prediction of relevant properties. However, computing in details the physicochemical properties of metal–organic nanohybrids proves to be a difficult computational challenge. While the nature of the organic component requires a full quantum-mechanical (QM) treatment in order to be predictive and properly interpret molecular electronic properties, the size of the inorganic component renders the QM treatment computationally too expensive, hampering the employment of an homogeneous technique for supramolecular nanohybrid systems. For this reason, hybrid models have been developed combining a QM level treatment and a classical electromagnetism approach for molecules and the inorganic nanostructures, respectively. (30) In particular, state-of-the-art methodologies (31−35) treat the inorganic metallic component within a polarizable continuum model (PCM) approach, (36−38) i.e., as a continuous body characterized by its own frequency dependent dielectric function, while the molecular excitations can be determined by employing the time-dependent density functional theory (TDDFT). (39) Still, such hybrid computational modeling has not been reported for nanohybrids containing gold bipyramids, despite their relevant optical properties, their potential application as inorganic components in nanohybrid materials and their use as seeds for regrowing nanoparticles with various shapes. (11,20) The first theoretical challenge is to achieve an accurate and comprehensive modeling of pure GNBs, with previously reported simulations possessing limited accuracy both in describing structural details of various bipyramidal morphologies (e.g., smooth edges and tips) and in reproducing experimental LSPR absorption maxima as a function of the GNB’s AR. (12,14)
The theoretical treatment of LSPR in noble metal nanoparticles has long been established, (1,40) with methods that generally stem from classical electrodynamics Mie scattering theory, which gives analytical results only for simple spheres, such as spheres, imposing the use of numerical approaches for other shapes, like bipyramids. These methods generally include modeling of the environment around the metallic object as a continuum dielectric, usually proving to be effective in the description of LSPR spectra and near-field enhancements, in the interpretation of experimental observations (41) and in the design new materials for technological applications. (42) For GNBs, the discrete dipole approximation (DDA) (43) has been used to solve numerically the electromagnetic (EM) scattering equations (12,14) for a nano-object, where the three-dimensional solid particle is replaced by a finite array of polarizable point dipoles, with the spacing between the dipoles small compared to the wavelength of the incoming light. Each of these dipoles has an oscillating polarization in response to both an incident plane wave and the electric fields due to all of the other dipoles in the array. This numerical method allows the employment of arbitrary shapes, instead of being restricted to a small number of regular shapes like spheres, spheroids and infinite cylinders, along with the opportunity to introduce anisotropies inside the particle. However, this approach has some limitations in terms of affordability with respect to the number of dipoles employed, which does not allow for a very detailed representation of complex shapes. An alternative approach is to solve numerically the EM scattering equations exploiting the boundary element method (BEM). BEM, in general, is a numerical computational method of solving linear partial differential equations which were reformulated as integral equations. Rather than solving a numerical problem of a set of partial differential equations defined throughout the whole space, BEM maps the problem into an integral equation defined over the boundary of the system. In practice, a three-dimensional problem (like that solved in DDA method) is mapped onto a two-dimensional surface mesh, thus demanding a computational cost that scales with the size of the surface instead of the volume and then allowing for a greater control of the details of the chosen object. The solid particle is represented by a set of apparent surface charges and currents that can be used to solve EM scattering equations and compute the absorption and scattering of the dielectric object. (44−47) In particular, the metal response to the electric fields can be determined using the widespread assumption that the metal behaves as a dielectric for time dependent fields. Within this approach, the solvent can still be described as a continuum dielectric which occupies all the space free from the metal specimen.
In this work, we report the unprecedented modeling of GNBs by means of the BEM approach, comparing the resulting optical properties with DDA results and with a large set of experimental data, and we introduce a BEM-based modeling for future studies of organic–inorganic nanohybrids based on GNBs. We first focused on providing a comprehensive theoretical investigation of GNBs to link the LSPR absorption maxima to their morphological characteristics, i.e., various kinds of tips or possible base geometries, even when a distribution of GNBs (as synthetic nanoparticles usually appear in solution) is taken into account, thus determining the effect of the statistical distribution of geometrical parameters on the LSPR absorption of GNBs with respect to single GNBs. If the size of the NP is comparable to the size of the incident wave, then the full-EM equation is needed. If the size of the NP is much smaller than the incident wave, the full-EM (i.e., Maxwell) equations reduce to just a time dependent Poisson problem. This approximation is called quasi-static (or dipole) approximation (QSA), and it allows straightforwardly to couple organic chromophores described by first principle approaches with the metal NP. (30) However, within this approximation the particle size does not affect the computed LSPR absorption energy but only its intensity, making the QSA simulations of LSPR absorption energies independent from the particle sizes. Here, we employed the BEM approach within the QSA to determine the maximum size of the GNBs that can be modeled with reasonable accuracy, paving the way to the study of GNB-based nanohybrids optical properties.

Computational Details

ARTICLE SECTIONS
Jump To

Numerical computations of LSPR absorption maxima have been performed using two different codes: Scuff-em, a free, open-source software implementation of the BEM method developed by Reid et al. (47,48) This code solves the full-EM scattering equations, taking into account explicitly the size of the GNBs; and TDPLAS, (49) that is routinely used to compute interaction between molecular electron densities and metal NPs, and presently relies on the QSA approximation.
The BEM approach implemented in Scuff-em makes use of effective (or apparent) electric (K) and magnetic (N) current densities located on the boundary of the nanoparticle and tangential to it to express the EM fields (electric and magnetic field) everywhere in space. By imposing the proper EM boundary conditions at the NP surfaces, it turns out that such current densities should satisfy the linear equation:
(1)
where the integrals are on the NP boundary from the inside (∂Ωr) and from the outside (∂Ωs) of the NP and ΓXX,t is the proper electromagnetic dyadic Green’s function (X=Electric or Magnetic) for the interior (t = r) or the exterior (t = s) of the NP evaluated at the incident EM field frequency. ΓXX,t contains the dependence on the NP dielectric function ϵ(ω) (as well as on the magnetic susceptibility μ(ω), that however we take equal to the vacuum susceptibility μ0 for all frequencies since gold is not a magnetic material); we omit dependence on frequency in the formulas for clarity. The expressions of ΓXX,t can be found in ref (50). Einc and Hinc in eq 1 are the components tangential (∥) to the NP surface of the incident electric and magnetic fields. To numerically solve eq 1, the surface integrals are discretized in finite boundary elements, and the surface current densities K, N are expanded on a finite basis set:
(2)
Converting eq 1 in a standard linear equation problem:
(3)
with
(4)
and
(5)
Z0 is the impedance of free space and |fα⟩ is a vector of basis set elements that are by construction tangent to the surface. Scuff-em makes use in particular of the RWG basis set. (51) Once the vectors K and N have been obtained, they can be used to calculate the EM field power absorbed by the NP (Pabs) via the calculation of the Poynting vector flux through the NP surface:
(6)
TDPLAS exploits the QSA where it is possible to simplify the full EM BEM eq 1. In particular, the theory can be reformulated in terms of apparent charge densities (rather than currents). An account on how to move from full EM BEM equations to QSA ones is given, e.g., in ref (52). In TDPlas implementation, the surface charge density is assumed to be constant on each boundary element i, giving rise to a total charge qi. Such charges are also obtained by solving a linear equation leading to (49)
(7)
where gi(ω) is given by −Einc·si (si is the center of the boundary element i) and:
(8)
is the NP response matrix containing the dielectric function ϵ(ω), the diagonal matrix A collects the areas of the boundary elements and the S and D matrices are related to electric field and electric potential integrals:
(9)
for out of diagonal elements (for diagonal elements see ref (31)). nj is the unit vector normal to the boundary element j and pointing outward. Since ϵ(ω) is complex, also q are complex. Once the charges q are obtained for a unit electric field set along a given direction e, a polarizability elements αfe can be obtained as αfe = ∑fiqi and the overall absorbed power at the frequency ω can be obtained in atomic units as where is the imaginary part and Tr is the trace of.
GNBs’ modeling is based on experimental geometrical parameters extracted from transmission electron microscopy (TEM) data of synthetic GNBs, as reported by Chateau et al., (12) Sánchez-Iglesias et al., (20) and Chateau et al. (21) Among the GNBs synthesized by Chateau et al. (12) we have considered the NPs for which the DDA simulated and the experimental LSPR maxima were reported and whose volumes were explicitly indicated, as determined from TEM images. For the GNBs synthesized by Sánchez-Iglesias et al., (20) we selected the samples where the NPs appeared less aggregated, allowing a more accurate image analysis of the TEM pictures (i.e., the GNBs reported in Figure S10a–c in the Supporting Information of ref (20)).
The analysis of TEM images show that the selected GNBs could have either a pentagonal (bipyramids) or a spherical (bicones) base shape. Indeed, the penta-twinned seeds used for the synthesis of GNBs have pentagonal base but GNBs’ edges are generally strongly smoothed, as observed in 3D tomographic electron microscopy at the single-particle level and in high-resolution TEM images. (11,16,53) Therefore, three GNB model shapes have been here considered: (i) bipyramid with pentagonal base and spherical tip (bipyramid spherical, BPS); (ii) bicone with spherical tip (bicone spherical, BCS); and (iii) bicone with flat tip (bicone flat, BCF). The model shapes and their meshes, built using the code GMSH, (54) are shown in Figure 1a. The smoothed edges are necessary to comply with the BEM method requirement of continuously differentiable surfaces and to minimize as much as possible the nonphysical effects arising from infinitely sharp edges, providing realistic descriptions of tip or edge effects. While for the GNBs synthesized by Chateau et al. (12) the shape (biconical), the AR (including the truncation level from ideal bipyramids), and the size of GNBs (in terms of volume) were explicitly reported, for the remaining GNBs, these parameters were extracted from the TEM images using the image analysis tools provided by ImageJ. (55) The GNBs models have been thus determined, for all shapes, from three experimental geometrical parameters for each single NP: (i) the width of the bipyramidal base (w); (ii) the AR (h/w), where w and height (h) of single particles were considered to be the Feret’s minimum and maximum diameter, respectively; (56) and (iii) the ”ideal” AR (h′/w), namely ARId, as computed by approximating the 2D projection of the GNB to a truncated rhombus (or the union of two symmetric isosceles trapezoids), as depicted in Figure 1b. By simple geometrical considerations, the ARId can be rewritten as a function of parameters that can be directly extracted from TEM images, reading
where A is the 2D projected area of the GNBs as obtained from ImageJ analysis of the TEM images. To build the GNB model of a given sample (associated with an experimental LSPR spectrum), thus, we extracted from its TEM image the w, h and A values for each NP, allowing computing the average ⟨w⟩, ⟨AR⟩, and ⟨ARId⟩ over the full set of NPs present in the sample. These three averaged parameters allow the construction of an average GNB structure for each sample. The standard deviations of these structural parameters, i.e., σ(⟨w⟩), σ(⟨AR⟩), and σ(⟨ARId⟩), were also computed, providing a way to estimate the geometrical uncertainty over the simulated LSPR values, due to the use of a single, average GNB structure as representative of each sample. A full list of the GNBs’s parameters is reported in the SI, see Table S2. Finally, for the dielectric properties of gold, we employed the analytic approach reported in ref (57and58).

Results and Discussion

ARTICLE SECTIONS
Jump To

Assessment of GNBs’ Model Shapes

Figure 2a illustrates the effects of AR and base size on the LSPR absorption maxima of GNBs particles with the same ARId and shows the comparison between experimental and theoretical values computed solving the full-EM scattering equations, i.e., without limitation in accuracy due to the specific particle sizes and light wavelengths. Previously reported experimental data and DDA simulations (12) are compared with our BEM results, computed with the Scuff-em code for the three model shapes (i.e., BPS, BCS and BCF) here considered (see Computational Details section). First, we point out that the DDA results (obtained with a biconical shape) match experimental LSPR absorption maxima only for larger ARs, while at smaller ARs the absolute errors increase and also the trend is not perfectly reproduced (see Figure 2a and Table S1 in the SI). This behavior is most likely due to the limitations of DDA approach in defining with great details the GNB shape (due to the computational cost of the 3D mapping), which become marginal only when the AR and the size of the NP increase significantly. Our BEM computations, instead, do not suffer of such limitations and indeed they provide excellent agreement with experimental data even at small ARs if the BCS model, i.e., with round base and spherical tips, is chosen. Notably, if the pentagonal base is considered instead, i.e., the BPS model, the BEM simulations reproduce properly the experimental trends but the LSPR maxima are offset by a consistent and systematic red-shift, always larger than 0.07 eV (see Table S1 in the SI). This systematic red-shift with respect to the biconical BCS shape is due to the presence of the extended lateral edges associated with the pentagonal base, an effect that is clearly out-weighting the (small) blue-shift one should expect when the overall particle volume reduces by passing from a biconical to a bipyramidal shape. In fact, the decrease of metal NPs size is generally accompanied by a blue-shift of the plasmonic resonance absorption.

Figure 2

Figure 2. (a) Effects of AR on the LSPR absorption maxima of GNBs particles. Experimental data (12) (in red) are compared with computed values solving the full-EM equations, including previously reported (12) DDA simulations (in blue), and our BEM results adopting the BPS (in yellow), the BCS (in green), and the BCF (in violet) models. (b) Experimental (normalized) absorption electronic spectra of GNBs featuring the various ARs reported in (a), including AR = 2.25 (in dark gray), AR = 2.40 (in red), AR = 2.75 (in blue), and AR = 3.10 (in green), are compared with spectra computed by our BEM approach using the BCS model.

By comparing the two biconical models, BCS and BCF, it is possible to define the effect of flattening the GNB’s tip, which showed to feature two main contributions: (i) the increase of volume particle passing from a spherical to a flat tip for GNBs, which blue-shifts the BCS’s LSPR absorption with respect to BCF and (ii) the effect of surface sharpening due to the presence of tips, which induces a red-shift of the plasmonic resonance. As shown in Figure 2a, when the GNB with the smallest AR is considered (i.e., AR = 2.25), which also features the smallest base size (i.e., w = 11 nm), the volume variation due to change of tip shape is significant if related to the overall, relatively small particle volume, indicating the former ”volume” effect is probably the dominant one. In fact, this smallest GNB also features a significant truncation percentage of 38% (see Table S1), defined as (1 – h/h′) in %, and thus possesses a blunt tip, this minimizing the tip effect independently on the tips shape (spherical or flat). Overall, this explains why, in this case, the BCS’s absorption is blue-shifts with respect to the BCF one. Conversely, as the AR and the base size are slightly increased (i.e., to AR = 2.40 and w = 15 nm) and the truncation percentage reduced to 27%, the BCS absorption is red-shifted with respect to the BCF model. This result indicates that, in this case, the volume modification has a smaller contribution than in the previous case (reducing the LSPR blue-shift) and the presence of a sharpened tip determines a sizable tip effect (and thus a red-shift) that is clearly larger for the sharp spherical tip than for the planar flat tip. This trend is essentially maintained if the AR and the base size are further increased and the truncation percentage further reduced, but the BCS vs BCF absorption difference tends to shrink as a consequence of the reducing size of the tip with the elongation of the GNB particle, which will make fading the importance of the choice of the tip shape as the AR increases, see Figure 2a.
The BCS model thus provides the best agreement with experimental evidence and hereafter will be used as the reference shape for all remaining simulations. In Figure 2b, the experimental (normalized) absorption spectra of the GNBs at different ARs are compared with BEM simulations using the BCS reference shape. The experimental plasmon resonances are, as expected, strongly prominent over the interband transitions of bulk gold lying at shorter wavelengths (i.e., in the UV spectral window not shown in Figure 2b) while weaker absorptions are found at wavelengths in the 500–600 nm range. Our results suggest that in this spectral region several contributions give rise to such absorptions: (i) the presence of spherical spheres (resulting as byproduct of the GNBs synthesis); (ii) the transverse plasmon resonances due to the interacting light polarized along the short axis of GNBs; and (iii) multipolar resonance absorption peaks. As shown in Figure S4 in the SI, we checked the effects of the polarization of the exciting electromagnetic field over the plasmonic resonances, and we found that, for all GNBs here considered, the transverse resonances are located at 506 nm. Moreover, for the largest GNB (AR = 3.1) an additional absorption peak at around 600 nm due to multipolar resonance can be detected in the simulated spectrum (which disregard contributions from spheres and transverse resonances), suggesting that part of the absorption detected experimentally for this large, elongated GNB, feature such of a contribution. Notably, experimental spectra of even more elongated GNBs than that with AR = 3.1 reported in ref (21) feature clearer signatures of this multipolar contribution, at correspondingly more red-shifted wavelenghts. Finally, as shown in Figure S3, the intensities of the plasmonic bands depend on the volume of the GNBs, i.e., the greater the volume of the particle the more intense is its absorption band, in line with what is known for simpler shapes.

Base Size and Truncation Effects

Given the good agreement between experimental and theoretical LSPR spectra computed using the BEM method for the BSC model, with better accuracy than previously reported DDA computations, we performed simulations on a larger set of synthetic GNBs featuring a wider range of base sizes (up to ca. 37 nm), AR and ARId (low to 1.84 and 3.16, respectively) and tip truncation percentages (up to 48%), as reported in Figure 3 and Table S2 in the SI. As reported in Table 1, the LSPR absorption maxima computed, within the full-EM framework, with the BEM approach and BCS geometrical model feature overall a good agreement with all experimental data, with energy difference discrepancies always below 0.06 eV (see Figure 3b). The simulated LSPR absorptions are computed using an average GNB particle, i.e. averaging the geometrical parameters, i.e. height ⟨h⟩, aspect ratios ⟨AR⟩ and ⟨ARId⟩, that can be extracted from the experimental TEM images, as described in the Computational Details section.
Table 1. Experimental (in nm) and Computed (in nm and eV) LSPR Absorption Maxima for All GNBs Samplesa
exp. LSPR (nm)theory LSPR (nm)theory LSPR (eV)ΔE (eV) (Exp–Comp)geom. var. ξ (eV)W (nm)
639b6221.993–0.050.0834.3
650b6601.8780.030.136.7
650b6261.981–0.070.0523.0
677c6741.839–0.0080.0518.1
698c6831.815–0.040.0720.1
714b7261.7080.030.0421.9
732c7551.6420.050.0736.3
787b7851.579–0.0040.0626.3
a

Theoretical values are computed with BEM method and BSC morphology, using a single GNB particle model based on the average of geometrical parameters (reported in Table S2 in the SI) extracted from experimental TEM images of each GNB sample. Differences between experimental and computed values are reported as ΔE, in eV. Geometrically-induced maximum variations on the computed LSPR energies (ξ, in eV) observed for the ensembles of GNBs in each sample are also reported.

b

Data from ref (21).

c

Data from ref (20).

Figure 3

Figure 3. (a) Comparison of experimental (black solid circles) LSPR absorption maxima and theoretical values (green solid circles), computed with average single-particle BCS models, for all the GNBs considered. These GNBs feature a large geometrical variability, including different base sizes and various ARIds (in colored text) and truncation percentages (in square brackets). For the GNBs TEM images analyzed in this work, the largest difference (ξ) in LSPR energies, with respect to the average single-particle computations, are reported as (±ξ) error bars of theoretical values. (b) Correlation plot of experimental vs simulated absorption maxima.

In order to assess the uncertainty of the computed LSPR associated with the geometrical variance of the synthetic GNBs, we analyzed TEM images from the samples reported in refs (20and21). We used the standard deviations (±1σ) of the three geometrical parameters ⟨w⟩, ⟨AR⟩, and ⟨ARId⟩ (see the Computational Details section) to build eight (23) GNBs’ “extreme cases” models for each sample and we computed their individual LSPR absorption maxima to have a rough estimate of their variation as a function of these parameters. In Figure 3a we report the largest difference (ξ) in LSPR energy with respect to the average GNBs as (±ξ) error bars, showing that the geometrical uncertainty derived from the extraction of the GNBs’ parameters from TEM images is significantly large at already ±1σ from the average and generally even larger than the difference between the simulated and experimental LSPR absorption maxima. This fact suggests that the sensitivity and accuracy of the BEM approach is sufficiently high that the results of the simulations are mostly affected by the choice of the NP extracted from the TEM image and from the corresponding geometrical parameters. This outcome highlights the importance of averaging all of the geometrical parameters extracted from TEM images rather than selecting one random single NP from the sample. Finally, we tested how much the result of such protocol differs from the explicit computations of all the LSPR absorption maxima of any NP present in the sample. Since these explicit ensemble computations have a computational cost significantly larger than a (averaged) single particle computation, we performed the simulations for only one representative sample. In particular, we selected the GNB sample featuring both a large geometrically induced variation on the LSPR energy (i.e., a large ξ) and one of the largest discrepancy between the experimental LSPR maximum and the single-particle theoretical value (i.e., a “large” ΔE(Exp–Comp)). We have thus computed the LSPR spectrum of the full ensemble of the GNB featuring an experimental LSPR maximum at 639 nm (i.e., 1.940 eV), featuring ξ = 0.08 eV and ΔE(Exp–Comp) = −0.05 eV, and compared it with the spectrum of the (averaged) single-particle, showing a deviation of LSPR maximum of only 0.01 eV and an almost identical LSPR spectra (with obviously a broader line shape in the full ensemble spectrum), as reported in Figures S1 and S2 in the SI.

Base Size Limits for QSA Computations

After validating BEM approach to EM scattering for GNBs particles with broad geometrical variety in terms of base size, ARId, and tip truncation (i.e., AR), we investigated how well QSA results could compare against full-EM computations. The assessment of QSA computations, in fact, is crucial to allow theoretical estimation of interactions between GNBs and molecular transition dipole moments in nanohybrids materials. Figure 4 reports the comparison of the experimental and full-EM simulations data described in section Assessment of GNBs’ model shapes with the QSA estimates of the LSPR absorption maxima using the BSC model. The outcome shows good agreement between the QSA results and the full-EM (and the experimental) data. As expected from the intrinsic limitation of the QSA approach (see the Introduction), the accordance worsens while increasing the GNBs’ particle size, which (in the reported cases) is associated with an increase of the AR.

Figure 4

Figure 4. Comparison of experimental (black solid circles) LSPR absorption maxima and theoretical values for the BCS models computed with full-EM (dark green solid circles) and QSA (light green solid circles) methods.

In order to better evaluate the error of the QSA calculations as a function of the GNBs’ size, we built a series of GNB models with increasing base size w, keeping the AR and the ARId fixed and then we compared the LSPR absorption maxima computed with the QSA and the full-EM approaches. The outcome of these computations is reported in Figure 5a and it shows how the particle size-independent of QSA estimates of the LSPR absorption maxima converge to the full-EM results for decreasing GNBs’ base sizes.

Figure 5

Figure 5. (a) Comparison between full-EM (colored solid circles) and QSA calculations (colored flat lines) of LSPR absorption maxima, as a function of the GNBs base sizes and for different ARs (color coded from 2.25 to 3.10) and fixed ARId (equal to 3.6), with empty colored circles reporting the experimental data (at specific base sizes). (b) Energy errors (in percentage) on the LSPR absorptions computed with QSA method with respect to the full-EM reference values, as a function of the GNBs’ base size, for different ARs (color coded from 2.25 to 3.10) and fixed ARId (equal to 3.6).

As previously mentioned, the full-EM computations show how the increase of the particle base size (that, for a given AR and ARId, corresponds to an increase of particle volume) is associated with a red-shift of the LSPR absorption, limiting the accuracy of the QSA computations to a finite range of base sizes, which depends on the chosen AR and ARId. Thus, in order to evaluate the extension of this range we computed the error (in percentage) in QSA absorption energies with respect to full-EM results as a function of the base size, for each AR and for a fixed ARId (equal to 3.6, see Table S1). As shown in Figure 5b, to stay within an error of 5% in LSPR absorption energy, the maximum base size should be lower than 36 nm for this class of GNBs. Considering that the typical base sizes of GNB objects are below 40 nm (as all GNBs investigated here, see Tables S1 and S2 in the SI) and that ARIds smaller than the one selected (that eventually would slightly reduce the maximum base size value) are not quite common, our result definitively represents a good estimation of the typical base size limit that one can reach one modeling GNBs LSPR absorption spectra within the QSA methodology.

Conclusions

ARTICLE SECTIONS
Jump To

In this work we proved that classical electrodynamics full-EM scattering calculations, solved numerically employing BEM approaches, can accurately describe the LSPR absorption properties of synthetic GNBs’ samples, if biconical models with spherical tips are employed. BEM simulations are sensitive enough to the GNBs’ geometries that can provide results fully consistent with experimental evidence and with higher accuracy than DDA computations. The GNBs’ models can be constructed by extracting geometrical parameters from TEM images of synthetic samples, using a single (random) representative particle, the full set of particles, or an average single-particle obtained by averaging the geometrical parameters over the ensemble of particles. Our results demonstrate that the latter option provides results very similar to the full-ensemble computations, with a gain in computational cost, and suggest that random selection of a single nanoparticle from TEM images to model a GNB could entail a large error on the computed LSPR absorption. Finally, the performances of QSA simulations of GNBs’ LSPR absorption spectra have been compared to the accurate full-EM approach, since reliable QSA computations would allow appropriate characterization of the interactions between molecular electron densities and gold nanoparticles that occur in GNB-based organic nanohybrids. Despite QSA computations of LSPR absorption energies are not sensitive to the particle sizes, they proved to be in good agreement with full-EM calculations for an extended range of GNBs’ base sizes, proving the applicability of this approach to synthetic GNBs and paving the way to the study of the optical properties of GNB-based nanohybrid materials.

Supporting Information

ARTICLE SECTIONS
Jump To

The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.0c00269.

  • Additional results as shown in Figures S1–S4, Tables S1 and S2, and refs 1 and 2 (PDF)

Terms & Conditions

Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Author Information

ARTICLE SECTIONS
Jump To

  • Corresponding Authors
    • Stefano Baroni - SISSA—Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy Email: [email protected]
    • Stefano Corni - Dipartimento di Scienze Chimiche, Università di Padova, 35131 Padova, ItalyIstituto di Nanoscienze, Consiglio Nazionale delle Ricerche CNR-NANO, 41125 Modena, ItalyOrcidhttp://orcid.org/0000-0001-6707-108X Email: [email protected]
    • Marco Garavelli - Dipartimento di Chimica Industriale “Toso Montanari”, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, ItalyOrcidhttp://orcid.org/0000-0002-0796-289X Email: [email protected]
    • Ivan Rivalta - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, FranceDipartimento di Chimica Industriale “Toso Montanari”, Università degli Studi di Bologna, Viale del Risorgimento 4, I-40136 Bologna, ItalyOrcidhttp://orcid.org/0000-0002-1208-602X Email: [email protected]
  • Authors
    • Jacopo Marcheselli - SISSA—Scuola Internazionale Superiore di Studi Avanzati, Via Bonomea 265, 34136 Trieste, Italy
    • Denis Chateau - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
    • Frederic Lerouge - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, FranceOrcidhttp://orcid.org/0000-0003-2909-527X
    • Patrice Baldeck - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
    • Chantal Andraud - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
    • Stephane Parola - Laboratoire de Chimie UMR 5182, CNRS, Université Lyon 1, Univ Lyon, Ens de Lyon, F-69342 Lyon, France
  • Notes
    The authors declare no competing financial interest.

Acknowledgments

ARTICLE SECTIONS
Jump To

I.R., M.G., F.L., P.B., C.A., and S.P. thank ENS Lyon (France) and the financial support from the LABEX iMUST (ANR-10-LABX-0064) of Université de Lyon, within the program “Investissements d’Avenir” (ANR-11-IDEX-0007) operated by the French National Research Agency (ANR). S.C. acknowledge funding from the ERC under the Grant ERC-CoG-681285 TAME-Plasmons. I.R. acknowledges the use of HPC resources of the “Pôle Scientifique de Modélisation Numérique” at the ENS Lyon, France.

References

ARTICLE SECTIONS
Jump To

This article references 58 other publications.

  1. 1
    Quinten, M. Optical Properties of Nanoparticle Systems: Mie and beyond; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2011.
  2. 2
    Mayer, K. M.; Hafner, J. H. Localized Surface Plasmon Resonance Sensors. Chem. Rev. 2011, 111, 38283857,  DOI: 10.1021/cr100313v
  3. 3
    Lermé, J.; Baida, H.; Bonnet, C.; Broyer, M.; Cottancin, E.; Crut, A.; Maioli, P.; Del Fatti, N.; Vallée, F.; Pellarin, M. Size Dependence of the Surface Plasmon Resonance Damping in Metal Nanospheres. J. Phys. Chem. Lett. 2010, 1, 29222928,  DOI: 10.1021/jz1009136
  4. 4
    Grillet, N.; Manchon, D.; Bertorelle, F.; Bonnet, C.; Broyer, M.; Cottancin, E.; Lermé, J.; Hillenkamp, M.; Pellarin, M. Plasmon Coupling in Silver Nanocube Dimers: Resonance Splitting Induced by Edge Rounding. ACS Nano 2011, 5, 94509462,  DOI: 10.1021/nn2041329
  5. 5
    Hao, F.; Nehl, C. L.; Hafner, J. H.; Nordlander, P. Plasmon Resonances of a Gold Nanostar. Nano Lett. 2007, 7, 729732,  DOI: 10.1021/nl062969c
  6. 6
    Guerrero-Martínez, A.; Barbosa, S.; Pastoriza-Santos, I.; Liz-Marzán, L. M. Nanostars shine bright for you: Colloidal synthesis, properties and applications of branched metallic nanoparticles. Curr. Opin. Colloid Interface Sci. 2011, 16, 118127,  DOI: 10.1016/j.cocis.2010.12.007
  7. 7
    Burda, C.; Chen, X.; Narayanan, R.; El-Sayed, M. A. Chemistry and Properties of Nanocrystals of Different Shapes. Chem. Rev. 2005, 105, 10251102,  DOI: 10.1021/cr030063a
  8. 8
    Link, S.; El-Sayed, M. A. Spectral Properties and Relaxation Dynamics of Surface Plasmon Electronic Oscillations in Gold and Silver Nanodots and Nanorods. J. Phys. Chem. B 1999, 103, 84108426,  DOI: 10.1021/jp9917648
  9. 9
    Carbó-Argibay, E.; Rodríguez-González, B.; Pastoriza-Santos, I.; Pérez-Juste, J.; Liz-Marzán, L. M. Growth of pentatwinned gold nanorods into truncated decahedra. Nanoscale 2010, 2, 23772383,  DOI: 10.1039/c0nr00239a
  10. 10
    Saute, B.; Narayanan, R. Solution-based direct readout surface enhanced Raman spectroscopic (SERS) detection of ultra-low levels of thiram with dogbone shaped gold nanoparticles. Analyst 2011, 136, 527532,  DOI: 10.1039/C0AN00594K
  11. 11
    Lee, J.-H.; Gibson, K. J.; Chen, G.; Weizmann, Y. Bipyramid-templated synthesis of monodisperse anisotropic gold nanocrystals. Nat. Commun. 2015, 6, 8671,  DOI: 10.1038/ncomms8571
  12. 12
    Chateau, D.; Liotta, A.; Vadcard, F.; Navarro, J. R. G.; Chaput, F.; Lermé, J.; Lerouge, F.; Parola, S. From gold nanobipyramids to nanojavelins for a precise tuning of the plasmon resonance to the infrared wavelengths: experimental and theoretical aspects. Nanoscale 2015, 7, 19341943,  DOI: 10.1039/C4NR06323F
  13. 13
    Kou, X.; Ni, W.; Tsung, C.-K.; Chan, K.; Lin, H.-Q.; Stucky, G.; Wang, J. Growth of Gold Bipyramids with Improved Yield and Their Curvature-Directed Oxidation. Small 2007, 3, 21032113,  DOI: 10.1002/smll.200700379
  14. 14
    Navarro, J. R. G.; Lerouge, F.; Micouin, G.; Cepraga, C.; Favier, A.; Charreyre, M. T.; Blanchard, N. P.; Lermé, J.; Chaput, F.; Focsan, M.; Kamada, K.; Baldeck, P. L.; Parola, S. Plasmonic bipyramids for fluorescence enhancement and protection against photobleaching. Nanoscale 2014, 6, 5138,  DOI: 10.1039/c3nr06425e
  15. 15
    Navarro, J. R. G.; Manchon, D.; Lerouge, F.; Blanchard, N. P.; Marotte, S.; Leverrier, Y.; Marvel, J.; Chaput, F.; Micouin, G.; Gabudean, A.-M.; Mosset, A.; Cottancin, E.; Baldeck, P. L.; Kamada, K.; Parola, S. Synthesis of PEGylated gold nanostars and bipyramids for intracellular uptake. Nanotechnology 2012, 23, 465602,  DOI: 10.1088/0957-4484/23/46/465602
  16. 16
    Navarro, J. R. G.; Manchon, D.; Lerouge, F.; Cottancin, E.; Lermé, J.; Bonnet, C.; Chaput, F.; Mosset, A.; Pellarin, M.; Parola, S. Synthesis, electron tomography and single-particle optical response of twisted gold nano-bipyramids. Nanotechnology 2012, 23, 145707,  DOI: 10.1088/0957-4484/23/14/145707
  17. 17
    Meena, S. K.; Sulpizi, M. From Gold Nanoseeds to Nanorods: The Microscopic Origin of the Anisotropic Growth. Angew. Chem., Int. Ed. 2016, 55, 1196011964,  DOI: 10.1002/anie.201604594
  18. 18
    Le Ru, E. C.; Grand, J.; Sow, I.; Somerville, W. R. C.; Etchegoin, P. G.; Treguer-Delapierre, M.; Charron, G.; Félidj, N.; Lévi, G.; Aubard, J. A Scheme for Detecting Every Single Target Molecule with Surface-Enhanced Raman Spectroscopy. Nano Lett. 2011, 11, 50135019,  DOI: 10.1021/nl2030344
  19. 19
    Rye, J.-M.; Bonnet, C.; Lerouge, F.; Pellarin, M.; Lermé, J.; Parola, S.; Cottancin, E. Single gold bipyramids on a silanized substrate as robust plasmonic sensors for liquid environments. Nanoscale 2018, 10, 1609416101,  DOI: 10.1039/C8NR03400A
  20. 20
    Sánchez-Iglesias, A.; Winckelmans, N.; Altantzis, T.; Bals, S.; Grzelczak, M.; Liz-Marzán, L. M. High-Yield Seeded Growth of Monodisperse Pentatwinned Gold Nanoparticles through Thermally Induced Seed Twinning. J. Am. Chem. Soc. 2017, 139, 107110,  DOI: 10.1021/jacs.6b12143
  21. 21
    Chateau, D.; Desert, A.; Lerouge, F.; Landaburu, G.; Santucci, S.; Parola, S. Beyond the Concentration Limitation in the Synthesis of Nanobipyramids and Other Pentatwinned Gold Nanostructures. ACS Appl. Mater. Interfaces 2019, 11, 3906839076,  DOI: 10.1021/acsami.9b12973
  22. 22
    Jiang, R.; Li, B.; Fang, C.; Wang, J. Metal/Semiconductor Hybrid Nanostructures for Plasmon-Enhanced Applications. Adv. Mater. 2014, 26, 52745309,  DOI: 10.1002/adma.201400203
  23. 23
    Navarro, J. R. Nanocarriers with ultrahigh chromophore loading for fluorescence bio-imaging and photodynamic therapy. Biomaterials 2013, 34, 83448351,  DOI: 10.1016/j.biomaterials.2013.07.032
  24. 24
    Campu, A.; Lerouge, F.; Chateau, D.; Chaput, F.; Baldeck, P.; Parola, S.; Maniu, D.; Craciun, A. M.; Vulpoi, A.; Astilean, S.; Focsan, M. Gold NanoBipyramids Performing as Highly Sensitive Dual-Modal Optical Immunosensors. Anal. Chem. 2018, 90, 85678575,  DOI: 10.1021/acs.analchem.8b01689
  25. 25
    Törmä, P.; Barnes, W. L. Strong coupling between surface plasmon polaritons and emitters: a review. Rep. Prog. Phys. 2015, 78, 013901,  DOI: 10.1088/0034-4885/78/1/013901
  26. 26
    Lakowicz, J. R. Radiative decay engineering 5: metal-enhanced fluorescence and plasmon emission. Anal. Biochem. 2005, 337, 171194,  DOI: 10.1016/j.ab.2004.11.026
  27. 27
    Geddes, C. D.; Lakowicz, J. R. Editorial: Metal-Enhanced Fluorescence. J. Fluoresc. 2002, 12, 121129,  DOI: 10.1023/A:1016875709579
  28. 28
    Pompa, P. P.; Martiradonna, L.; Torre, A. D.; Sala, F. D.; Manna, L.; De Vittorio, M.; Calabi, F.; Cingolani, R.; Rinaldi, R. Metal-enhanced fluorescence of colloidal nanocrystals with nanoscale control. Nat. Nanotechnol. 2006, 1, 126130,  DOI: 10.1038/nnano.2006.93
  29. 29
    Sau, T. K.; Rogach, A. L.; Jäckel, F.; Klar, T. A.; Feldmann, J. Properties and Applications of Colloidal Nonspherical Noble Metal Nanoparticles. Adv. Mater. 2010, 22, 18051825,  DOI: 10.1002/adma.200902557
  30. 30
    Mennucci, B.; Corni, S. Multiscale modelling of photoinduced processes in composite systems. Nat. Rev. Chem. 2019, 3, 315330,  DOI: 10.1038/s41570-019-0092-4
  31. 31
    Corni, S.; Tomasi, J. Enhanced response properties of a chromophore physisorbed on a metal particle. J. Chem. Phys. 2001, 114, 37393751,  DOI: 10.1063/1.1342241
  32. 32
    Corni, S.; Tomasi, J. Surface enhanced Raman scattering from a single molecule adsorbed on a metal particle aggregate: A theoretical study. J. Chem. Phys. 2002, 116, 11561164,  DOI: 10.1063/1.1428349
  33. 33
    Caricato, M.; Andreussi, O.; Corni, S. Semiempirical (ZINDO-PCM) Approach to Predict the Radiative and Nonradiative Decay Rates of a Molecule Close to Metal Particles. J. Phys. Chem. B 2006, 110, 1665216659,  DOI: 10.1021/jp0626418
  34. 34
    Vukovic, S.; Corni, S.; Mennucci, B. Fluorescence Enhancement of Chromophores Close to Metal Nanoparticles. Optimal Setup Revealed by the Polarizable Continuum Model. J. Phys. Chem. C 2009, 113, 121133,  DOI: 10.1021/jp808116y
  35. 35
    Andreussi, O.; Biancardi, A.; Corni, S.; Mennucci, B. Plasmon-Controlled Light-Harvesting: Design Rules for Biohybrid Devices via Multiscale Modeling. Nano Lett. 2013, 13, 44754484,  DOI: 10.1021/nl402403v
  36. 36
    Amovilli, C.; Barone, V.; Cammi, R.; Cancès, E.; Cossi, M.; Mennucci, B.; Pomelli, C. S.; Tomasi, J. In Adv. Quantum Chem.; Löwdin, P.-O., Ed.; Academic Press: New York, 1998; Vol. 32, pp 227261.
  37. 37
    Cancès, E.; Mennucci, B. New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals. J. Math. Chem. 1998, 23, 309326,  DOI: 10.1023/A:1019133611148
  38. 38
    Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 29993094,  DOI: 10.1021/cr9904009
  39. 39
    Marques, M. A. L., Maitra, N. T., Nogueira, F. M. S., Gross, E. K. U., Rubio, A., Eds.; Fundamentals of Time-Dependent Density Functional Theory; Lecture Notes in Physics; Springer-Verlag: Berlin, 2012.
  40. 40
    Maier, S. A. Plasmonics: Fundamentals and Applications; Springer US: New York, 2007.
  41. 41
    Sinzig, J.; Quinten, M. Scattering and absorption by spherical multilayer particles. Appl. Phys. A: Solids Surf. 1994, 58, 157162,  DOI: 10.1007/BF00332172
  42. 42
    Mahmoud, M. A.; El-Sayed, M. A. Gold Nanoframes: Very High Surface Plasmon Fields and Excellent Near-Infrared Sensors. J. Am. Chem. Soc. 2010, 132, 1270412710,  DOI: 10.1021/ja104532z
  43. 43
    Draine, B. T.; Flatau, P. J. Discrete-Dipole Approximation For Scattering Calculations. J. Opt. Soc. Am. A 1994, 11, 1491,  DOI: 10.1364/JOSAA.11.001491
  44. 44
    Garcia de Abajo, F. J.; Howie, A. Retarded field calculation of electron energy loss in inhomogeneous dielectrics. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 65, 115418,  DOI: 10.1103/PhysRevB.65.115418
  45. 45
    Myroshnychenko, V.; Rodriguez-Fernandez, J.; Pastoriza-Santos, I.; Funston, A. M.; Novo, C.; Mulvaney, P.; Liz-Marzan, L. M.; Garcia de Abajo, F. J. Modelling the optical response of gold nanoparticles. Chem. Soc. Rev. 2008, 37, 17921805,  DOI: 10.1039/b711486a
  46. 46
    Hohenester, U.; Trügler, A. MNPBEM - A Matlab toolbox for the simulation of plasmonic nanoparticles. Comput. Phys. Commun. 2012, 183, 370381,  DOI: 10.1016/j.cpc.2011.09.009
  47. 47
    Reid, M. T. H. scuff-em: A comprehensive and full-featured computational physics suite for boundary-element analysis , 2018.
  48. 48
    Reid, M. T. H.; Johnson, S. G. Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations. IEEE Trans. Antennas Propag. 2015, 63, 35883598,  DOI: 10.1109/TAP.2015.2438393
  49. 49
    Pipolo, S.; Corni, S. Real-Time Description of the Electronic Dynamics for a Molecule Close to a Plasmonic Nanoparticle. J. Phys. Chem. C 2016, 120, 2877428781,  DOI: 10.1021/acs.jpcc.6b11084
  50. 50
    Novotny, L.; Hecht, B. Principles of Nano-Optics; Cambridge University Press: Cambridge, U.K., 2006.
  51. 51
    Rao, S.; Wilton, D.; Glisson, A. Electromagnetic scattering by surfaces of arbitrary shape. IRE Trans. Antennas Propag. 1982, 30, 409418,  DOI: 10.1109/TAP.1982.1142818
  52. 52
    Pipolo, S.; Corni, S.; Cammi, R. The cavity electromagnetic field within the polarizable continuum model of solvation. J. Chem. Phys. 2014, 140, 164114,  DOI: 10.1063/1.4871373
  53. 53
    Canbek, Z. C.; Cortes-Huerto, R.; Testard, F.; Spalla, O.; Moldovan, S.; Ersen, O.; Wisnet, A.; Wang, G.; Goniakowski, J.; Noguera, C.; Menguy, N. Twinned Gold Nanoparticles under Growth: Bipyramids Shape Controlled by Environment. Cryst. Growth Des. 2015, 15, 36373644,  DOI: 10.1021/acs.cgd.5b00121
  54. 54
    Geuzaine, C.; Remacle, J.-F. Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Meth. Engng. 2009, 79, 13091331,  DOI: 10.1002/nme.2579
  55. 55
    Schindelin, J.; Rueden, C. T.; Hiner, M. C.; Eliceiri, K. W. The ImageJ ecosystem: An open platform for biomedical image analysis. Mol. Reprod. Dev. 2015, 82, 518529,  DOI: 10.1002/mrd.22489
  56. 56
    Merkus, H. G. Particle Size Measurements: Fundamentals, Practice, Quality; Particle Technology Series; Springer Netherlands: Amsterdam, 2009.
  57. 57
    Etchegoin, P. G.; Le Ru, E. C.; Meyer, M. An analytic model for the optical properties of gold. J. Chem. Phys. 2006, 125, 164705,  DOI: 10.1063/1.2360270
  58. 58
    Etchegoin, P. G.; Le Ru, E. C.; Meyer, M. Erratum: An analytic model for the optical properties of gold. J. Chem. Phys. 2006, 125, 164705,  DOI: 10.1063/1.2360270 .
    J. Chem. Phys. 2007, 127, 189901.

Cited By

This article is cited by 13 publications.

  1. Yibin Jiang, Abhishek Sharma, Leroy Cronin. An Accelerated Method for Investigating Spectral Properties of Dynamically Evolving Nanostructures. The Journal of Physical Chemistry Letters 2023, 14 (16) , 3929-3938. https://doi.org/10.1021/acs.jpclett.3c00395
  2. Tommaso Giovannini, Luca Bonatti, Piero Lafiosca, Luca Nicoli, Matteo Castagnola, Pablo Grobas Illobre, Stefano Corni, Chiara Cappelli. Do We Really Need Quantum Mechanics to Describe Plasmonic Properties of Metal Nanostructures?. ACS Photonics 2022, 9 (9) , 3025-3034. https://doi.org/10.1021/acsphotonics.2c00761
  3. Antonio Carone, Pablo Mariani, Anthony Désert, Marco Romanelli, Jacopo Marcheselli, Marco Garavelli, Stefano Corni, Ivan Rivalta, Stephane Parola. Insight on Chirality Encoding from Small Thiolated Molecule to Plasmonic Au@Ag and Au@Au Nanoparticles. ACS Nano 2022, 16 (1) , 1089-1101. https://doi.org/10.1021/acsnano.1c08824
  4. Laura Fabris. Gold Nanostars in Biology and Medicine: Understanding Physicochemical Properties to Broaden Applicability. The Journal of Physical Chemistry C 2020, 124 (49) , 26540-26553. https://doi.org/10.1021/acs.jpcc.0c08460
  5. Marta Rosa, Giulia Dall’Osto, Roberto Cammi, Stefano Corni. Quantum optimal control theory for a molecule interacting with a plasmonic nanoparticle. Theoretical Chemistry Accounts 2023, 142 (8) https://doi.org/10.1007/s00214-023-03025-5
  6. Sibanisankar Sahoo, Gour Mohan Das, Venkata Ramanaiah Dantham. Experimental and Numerical Investigation on Plasmonic Properties of Single-Sprouted Potato-Shaped Au-Ag Bimetallic Nanoparticles. Plasmonics 2023, 13 https://doi.org/10.1007/s11468-023-01931-9
  7. Jieli Lyu, Claire Goldmann, Cyrille Hamon, Doru Constantin. Extracting the morphology of gold bipyramids from small-angle X-ray scattering experiments via form factor modelling. Journal of Applied Crystallography 2023, 56 (1) , 214-221. https://doi.org/10.1107/S1600576722011669
  8. Yibin Jiang, Daniel Salley, Abhishek Sharma, Graham Keenan, Margaret Mullin, Leroy Cronin. An artificial intelligence enabled chemical synthesis robot for exploration and optimization of nanomaterials. Science Advances 2022, 8 (40) https://doi.org/10.1126/sciadv.abo2626
  9. Go Kawamura, Atsunori Matsuda. Nanomaterials for Localized Surface Plasmon Resonance-Related Optical Functionalities. 2022, 37-70. https://doi.org/10.1007/978-3-031-16518-4_2
  10. Si Yin Tee, Enyi Ye, Choon Peng Teng, Yuki Tanaka, Karen Yuanting Tang, Khin Yin Win, Ming-Yong Han. Advances in photothermal nanomaterials for biomedical, environmental and energy applications. Nanoscale 2021, 13 (34) , 14268-14286. https://doi.org/10.1039/D1NR04197E
  11. Amirmostafa Amirjani, S. K. Sadrnezhaad. Computational electromagnetics in plasmonic nanostructures. Journal of Materials Chemistry C 2021, 9 (31) , 9791-9819. https://doi.org/10.1039/D1TC01742J
  12. E. Coccia, J. Fregoni, C. A. Guido, M. Marsili, S. Pipolo, S. Corni. Hybrid theoretical models for molecular nanoplasmonics. The Journal of Chemical Physics 2020, 153 (20) https://doi.org/10.1063/5.0027935
  13. Giulia Dall’Osto, Gabriel Gil, Silvio Pipolo, Stefano Corni. Real-time dynamics of plasmonic resonances in nanoparticles described by a boundary element method with generic dielectric function. The Journal of Chemical Physics 2020, 153 (18) https://doi.org/10.1063/5.0022329
  • Abstract

    Figure 1

    Figure 1. (a) The GNBs’ models used in this work, including, from left to right the bipyramidal with pentagonal base and spherical tip (BPS), the biconical with spherical tip (BCS) and the biconical with flat tip (BCF) shapes. (b) A schematic representation of the basic geometric parameters of GNBs, with base length (w), height (h), and 2D projected area (in yellow) being extracted from TEM images, from which the nontruncated, ideal height (h′) can be derived.

    Figure 2

    Figure 2. (a) Effects of AR on the LSPR absorption maxima of GNBs particles. Experimental data (12) (in red) are compared with computed values solving the full-EM equations, including previously reported (12) DDA simulations (in blue), and our BEM results adopting the BPS (in yellow), the BCS (in green), and the BCF (in violet) models. (b) Experimental (normalized) absorption electronic spectra of GNBs featuring the various ARs reported in (a), including AR = 2.25 (in dark gray), AR = 2.40 (in red), AR = 2.75 (in blue), and AR = 3.10 (in green), are compared with spectra computed by our BEM approach using the BCS model.

    Figure 3

    Figure 3. (a) Comparison of experimental (black solid circles) LSPR absorption maxima and theoretical values (green solid circles), computed with average single-particle BCS models, for all the GNBs considered. These GNBs feature a large geometrical variability, including different base sizes and various ARIds (in colored text) and truncation percentages (in square brackets). For the GNBs TEM images analyzed in this work, the largest difference (ξ) in LSPR energies, with respect to the average single-particle computations, are reported as (±ξ) error bars of theoretical values. (b) Correlation plot of experimental vs simulated absorption maxima.

    Figure 4

    Figure 4. Comparison of experimental (black solid circles) LSPR absorption maxima and theoretical values for the BCS models computed with full-EM (dark green solid circles) and QSA (light green solid circles) methods.

    Figure 5

    Figure 5. (a) Comparison between full-EM (colored solid circles) and QSA calculations (colored flat lines) of LSPR absorption maxima, as a function of the GNBs base sizes and for different ARs (color coded from 2.25 to 3.10) and fixed ARId (equal to 3.6), with empty colored circles reporting the experimental data (at specific base sizes). (b) Energy errors (in percentage) on the LSPR absorptions computed with QSA method with respect to the full-EM reference values, as a function of the GNBs’ base size, for different ARs (color coded from 2.25 to 3.10) and fixed ARId (equal to 3.6).

  • References

    ARTICLE SECTIONS
    Jump To

    This article references 58 other publications.

    1. 1
      Quinten, M. Optical Properties of Nanoparticle Systems: Mie and beyond; Wiley-VCH Verlag GmbH & Co. KGaA: Weinheim, Germany, 2011.
    2. 2
      Mayer, K. M.; Hafner, J. H. Localized Surface Plasmon Resonance Sensors. Chem. Rev. 2011, 111, 38283857,  DOI: 10.1021/cr100313v
    3. 3
      Lermé, J.; Baida, H.; Bonnet, C.; Broyer, M.; Cottancin, E.; Crut, A.; Maioli, P.; Del Fatti, N.; Vallée, F.; Pellarin, M. Size Dependence of the Surface Plasmon Resonance Damping in Metal Nanospheres. J. Phys. Chem. Lett. 2010, 1, 29222928,  DOI: 10.1021/jz1009136
    4. 4
      Grillet, N.; Manchon, D.; Bertorelle, F.; Bonnet, C.; Broyer, M.; Cottancin, E.; Lermé, J.; Hillenkamp, M.; Pellarin, M. Plasmon Coupling in Silver Nanocube Dimers: Resonance Splitting Induced by Edge Rounding. ACS Nano 2011, 5, 94509462,  DOI: 10.1021/nn2041329
    5. 5
      Hao, F.; Nehl, C. L.; Hafner, J. H.; Nordlander, P. Plasmon Resonances of a Gold Nanostar. Nano Lett. 2007, 7, 729732,  DOI: 10.1021/nl062969c
    6. 6
      Guerrero-Martínez, A.; Barbosa, S.; Pastoriza-Santos, I.; Liz-Marzán, L. M. Nanostars shine bright for you: Colloidal synthesis, properties and applications of branched metallic nanoparticles. Curr. Opin. Colloid Interface Sci. 2011, 16, 118127,  DOI: 10.1016/j.cocis.2010.12.007
    7. 7
      Burda, C.; Chen, X.; Narayanan, R.; El-Sayed, M. A. Chemistry and Properties of Nanocrystals of Different Shapes. Chem. Rev. 2005, 105, 10251102,  DOI: 10.1021/cr030063a
    8. 8
      Link, S.; El-Sayed, M. A. Spectral Properties and Relaxation Dynamics of Surface Plasmon Electronic Oscillations in Gold and Silver Nanodots and Nanorods. J. Phys. Chem. B 1999, 103, 84108426,  DOI: 10.1021/jp9917648
    9. 9
      Carbó-Argibay, E.; Rodríguez-González, B.; Pastoriza-Santos, I.; Pérez-Juste, J.; Liz-Marzán, L. M. Growth of pentatwinned gold nanorods into truncated decahedra. Nanoscale 2010, 2, 23772383,  DOI: 10.1039/c0nr00239a
    10. 10
      Saute, B.; Narayanan, R. Solution-based direct readout surface enhanced Raman spectroscopic (SERS) detection of ultra-low levels of thiram with dogbone shaped gold nanoparticles. Analyst 2011, 136, 527532,  DOI: 10.1039/C0AN00594K
    11. 11
      Lee, J.-H.; Gibson, K. J.; Chen, G.; Weizmann, Y. Bipyramid-templated synthesis of monodisperse anisotropic gold nanocrystals. Nat. Commun. 2015, 6, 8671,  DOI: 10.1038/ncomms8571
    12. 12
      Chateau, D.; Liotta, A.; Vadcard, F.; Navarro, J. R. G.; Chaput, F.; Lermé, J.; Lerouge, F.; Parola, S. From gold nanobipyramids to nanojavelins for a precise tuning of the plasmon resonance to the infrared wavelengths: experimental and theoretical aspects. Nanoscale 2015, 7, 19341943,  DOI: 10.1039/C4NR06323F
    13. 13
      Kou, X.; Ni, W.; Tsung, C.-K.; Chan, K.; Lin, H.-Q.; Stucky, G.; Wang, J. Growth of Gold Bipyramids with Improved Yield and Their Curvature-Directed Oxidation. Small 2007, 3, 21032113,  DOI: 10.1002/smll.200700379
    14. 14
      Navarro, J. R. G.; Lerouge, F.; Micouin, G.; Cepraga, C.; Favier, A.; Charreyre, M. T.; Blanchard, N. P.; Lermé, J.; Chaput, F.; Focsan, M.; Kamada, K.; Baldeck, P. L.; Parola, S. Plasmonic bipyramids for fluorescence enhancement and protection against photobleaching. Nanoscale 2014, 6, 5138,  DOI: 10.1039/c3nr06425e
    15. 15
      Navarro, J. R. G.; Manchon, D.; Lerouge, F.; Blanchard, N. P.; Marotte, S.; Leverrier, Y.; Marvel, J.; Chaput, F.; Micouin, G.; Gabudean, A.-M.; Mosset, A.; Cottancin, E.; Baldeck, P. L.; Kamada, K.; Parola, S. Synthesis of PEGylated gold nanostars and bipyramids for intracellular uptake. Nanotechnology 2012, 23, 465602,  DOI: 10.1088/0957-4484/23/46/465602
    16. 16
      Navarro, J. R. G.; Manchon, D.; Lerouge, F.; Cottancin, E.; Lermé, J.; Bonnet, C.; Chaput, F.; Mosset, A.; Pellarin, M.; Parola, S. Synthesis, electron tomography and single-particle optical response of twisted gold nano-bipyramids. Nanotechnology 2012, 23, 145707,  DOI: 10.1088/0957-4484/23/14/145707
    17. 17
      Meena, S. K.; Sulpizi, M. From Gold Nanoseeds to Nanorods: The Microscopic Origin of the Anisotropic Growth. Angew. Chem., Int. Ed. 2016, 55, 1196011964,  DOI: 10.1002/anie.201604594
    18. 18
      Le Ru, E. C.; Grand, J.; Sow, I.; Somerville, W. R. C.; Etchegoin, P. G.; Treguer-Delapierre, M.; Charron, G.; Félidj, N.; Lévi, G.; Aubard, J. A Scheme for Detecting Every Single Target Molecule with Surface-Enhanced Raman Spectroscopy. Nano Lett. 2011, 11, 50135019,  DOI: 10.1021/nl2030344
    19. 19
      Rye, J.-M.; Bonnet, C.; Lerouge, F.; Pellarin, M.; Lermé, J.; Parola, S.; Cottancin, E. Single gold bipyramids on a silanized substrate as robust plasmonic sensors for liquid environments. Nanoscale 2018, 10, 1609416101,  DOI: 10.1039/C8NR03400A
    20. 20
      Sánchez-Iglesias, A.; Winckelmans, N.; Altantzis, T.; Bals, S.; Grzelczak, M.; Liz-Marzán, L. M. High-Yield Seeded Growth of Monodisperse Pentatwinned Gold Nanoparticles through Thermally Induced Seed Twinning. J. Am. Chem. Soc. 2017, 139, 107110,  DOI: 10.1021/jacs.6b12143
    21. 21
      Chateau, D.; Desert, A.; Lerouge, F.; Landaburu, G.; Santucci, S.; Parola, S. Beyond the Concentration Limitation in the Synthesis of Nanobipyramids and Other Pentatwinned Gold Nanostructures. ACS Appl. Mater. Interfaces 2019, 11, 3906839076,  DOI: 10.1021/acsami.9b12973
    22. 22
      Jiang, R.; Li, B.; Fang, C.; Wang, J. Metal/Semiconductor Hybrid Nanostructures for Plasmon-Enhanced Applications. Adv. Mater. 2014, 26, 52745309,  DOI: 10.1002/adma.201400203
    23. 23
      Navarro, J. R. Nanocarriers with ultrahigh chromophore loading for fluorescence bio-imaging and photodynamic therapy. Biomaterials 2013, 34, 83448351,  DOI: 10.1016/j.biomaterials.2013.07.032
    24. 24
      Campu, A.; Lerouge, F.; Chateau, D.; Chaput, F.; Baldeck, P.; Parola, S.; Maniu, D.; Craciun, A. M.; Vulpoi, A.; Astilean, S.; Focsan, M. Gold NanoBipyramids Performing as Highly Sensitive Dual-Modal Optical Immunosensors. Anal. Chem. 2018, 90, 85678575,  DOI: 10.1021/acs.analchem.8b01689
    25. 25
      Törmä, P.; Barnes, W. L. Strong coupling between surface plasmon polaritons and emitters: a review. Rep. Prog. Phys. 2015, 78, 013901,  DOI: 10.1088/0034-4885/78/1/013901
    26. 26
      Lakowicz, J. R. Radiative decay engineering 5: metal-enhanced fluorescence and plasmon emission. Anal. Biochem. 2005, 337, 171194,  DOI: 10.1016/j.ab.2004.11.026
    27. 27
      Geddes, C. D.; Lakowicz, J. R. Editorial: Metal-Enhanced Fluorescence. J. Fluoresc. 2002, 12, 121129,  DOI: 10.1023/A:1016875709579
    28. 28
      Pompa, P. P.; Martiradonna, L.; Torre, A. D.; Sala, F. D.; Manna, L.; De Vittorio, M.; Calabi, F.; Cingolani, R.; Rinaldi, R. Metal-enhanced fluorescence of colloidal nanocrystals with nanoscale control. Nat. Nanotechnol. 2006, 1, 126130,  DOI: 10.1038/nnano.2006.93
    29. 29
      Sau, T. K.; Rogach, A. L.; Jäckel, F.; Klar, T. A.; Feldmann, J. Properties and Applications of Colloidal Nonspherical Noble Metal Nanoparticles. Adv. Mater. 2010, 22, 18051825,  DOI: 10.1002/adma.200902557
    30. 30
      Mennucci, B.; Corni, S. Multiscale modelling of photoinduced processes in composite systems. Nat. Rev. Chem. 2019, 3, 315330,  DOI: 10.1038/s41570-019-0092-4
    31. 31
      Corni, S.; Tomasi, J. Enhanced response properties of a chromophore physisorbed on a metal particle. J. Chem. Phys. 2001, 114, 37393751,  DOI: 10.1063/1.1342241
    32. 32
      Corni, S.; Tomasi, J. Surface enhanced Raman scattering from a single molecule adsorbed on a metal particle aggregate: A theoretical study. J. Chem. Phys. 2002, 116, 11561164,  DOI: 10.1063/1.1428349
    33. 33
      Caricato, M.; Andreussi, O.; Corni, S. Semiempirical (ZINDO-PCM) Approach to Predict the Radiative and Nonradiative Decay Rates of a Molecule Close to Metal Particles. J. Phys. Chem. B 2006, 110, 1665216659,  DOI: 10.1021/jp0626418
    34. 34
      Vukovic, S.; Corni, S.; Mennucci, B. Fluorescence Enhancement of Chromophores Close to Metal Nanoparticles. Optimal Setup Revealed by the Polarizable Continuum Model. J. Phys. Chem. C 2009, 113, 121133,  DOI: 10.1021/jp808116y
    35. 35
      Andreussi, O.; Biancardi, A.; Corni, S.; Mennucci, B. Plasmon-Controlled Light-Harvesting: Design Rules for Biohybrid Devices via Multiscale Modeling. Nano Lett. 2013, 13, 44754484,  DOI: 10.1021/nl402403v
    36. 36
      Amovilli, C.; Barone, V.; Cammi, R.; Cancès, E.; Cossi, M.; Mennucci, B.; Pomelli, C. S.; Tomasi, J. In Adv. Quantum Chem.; Löwdin, P.-O., Ed.; Academic Press: New York, 1998; Vol. 32, pp 227261.
    37. 37
      Cancès, E.; Mennucci, B. New applications of integral equations methods for solvation continuum models: ionic solutions and liquid crystals. J. Math. Chem. 1998, 23, 309326,  DOI: 10.1023/A:1019133611148
    38. 38
      Tomasi, J.; Mennucci, B.; Cammi, R. Quantum Mechanical Continuum Solvation Models. Chem. Rev. 2005, 105, 29993094,  DOI: 10.1021/cr9904009
    39. 39
      Marques, M. A. L., Maitra, N. T., Nogueira, F. M. S., Gross, E. K. U., Rubio, A., Eds.; Fundamentals of Time-Dependent Density Functional Theory; Lecture Notes in Physics; Springer-Verlag: Berlin, 2012.
    40. 40
      Maier, S. A. Plasmonics: Fundamentals and Applications; Springer US: New York, 2007.
    41. 41
      Sinzig, J.; Quinten, M. Scattering and absorption by spherical multilayer particles. Appl. Phys. A: Solids Surf. 1994, 58, 157162,  DOI: 10.1007/BF00332172
    42. 42
      Mahmoud, M. A.; El-Sayed, M. A. Gold Nanoframes: Very High Surface Plasmon Fields and Excellent Near-Infrared Sensors. J. Am. Chem. Soc. 2010, 132, 1270412710,  DOI: 10.1021/ja104532z
    43. 43
      Draine, B. T.; Flatau, P. J. Discrete-Dipole Approximation For Scattering Calculations. J. Opt. Soc. Am. A 1994, 11, 1491,  DOI: 10.1364/JOSAA.11.001491
    44. 44
      Garcia de Abajo, F. J.; Howie, A. Retarded field calculation of electron energy loss in inhomogeneous dielectrics. Phys. Rev. B: Condens. Matter Mater. Phys. 2002, 65, 115418,  DOI: 10.1103/PhysRevB.65.115418
    45. 45
      Myroshnychenko, V.; Rodriguez-Fernandez, J.; Pastoriza-Santos, I.; Funston, A. M.; Novo, C.; Mulvaney, P.; Liz-Marzan, L. M.; Garcia de Abajo, F. J. Modelling the optical response of gold nanoparticles. Chem. Soc. Rev. 2008, 37, 17921805,  DOI: 10.1039/b711486a
    46. 46
      Hohenester, U.; Trügler, A. MNPBEM - A Matlab toolbox for the simulation of plasmonic nanoparticles. Comput. Phys. Commun. 2012, 183, 370381,  DOI: 10.1016/j.cpc.2011.09.009
    47. 47
      Reid, M. T. H. scuff-em: A comprehensive and full-featured computational physics suite for boundary-element analysis , 2018.
    48. 48
      Reid, M. T. H.; Johnson, S. G. Efficient Computation of Power, Force, and Torque in BEM Scattering Calculations. IEEE Trans. Antennas Propag. 2015, 63, 35883598,  DOI: 10.1109/TAP.2015.2438393
    49. 49
      Pipolo, S.; Corni, S. Real-Time Description of the Electronic Dynamics for a Molecule Close to a Plasmonic Nanoparticle. J. Phys. Chem. C 2016, 120, 2877428781,  DOI: 10.1021/acs.jpcc.6b11084
    50. 50
      Novotny, L.; Hecht, B. Principles of Nano-Optics; Cambridge University Press: Cambridge, U.K., 2006.
    51. 51
      Rao, S.; Wilton, D.; Glisson, A. Electromagnetic scattering by surfaces of arbitrary shape. IRE Trans. Antennas Propag. 1982, 30, 409418,  DOI: 10.1109/TAP.1982.1142818
    52. 52
      Pipolo, S.; Corni, S.; Cammi, R. The cavity electromagnetic field within the polarizable continuum model of solvation. J. Chem. Phys. 2014, 140, 164114,  DOI: 10.1063/1.4871373
    53. 53
      Canbek, Z. C.; Cortes-Huerto, R.; Testard, F.; Spalla, O.; Moldovan, S.; Ersen, O.; Wisnet, A.; Wang, G.; Goniakowski, J.; Noguera, C.; Menguy, N. Twinned Gold Nanoparticles under Growth: Bipyramids Shape Controlled by Environment. Cryst. Growth Des. 2015, 15, 36373644,  DOI: 10.1021/acs.cgd.5b00121
    54. 54
      Geuzaine, C.; Remacle, J.-F. Gmsh: A 3-D finite element mesh generator with built-in pre- and post-processing facilities. Int. J. Numer. Meth. Engng. 2009, 79, 13091331,  DOI: 10.1002/nme.2579
    55. 55
      Schindelin, J.; Rueden, C. T.; Hiner, M. C.; Eliceiri, K. W. The ImageJ ecosystem: An open platform for biomedical image analysis. Mol. Reprod. Dev. 2015, 82, 518529,  DOI: 10.1002/mrd.22489
    56. 56
      Merkus, H. G. Particle Size Measurements: Fundamentals, Practice, Quality; Particle Technology Series; Springer Netherlands: Amsterdam, 2009.
    57. 57
      Etchegoin, P. G.; Le Ru, E. C.; Meyer, M. An analytic model for the optical properties of gold. J. Chem. Phys. 2006, 125, 164705,  DOI: 10.1063/1.2360270
    58. 58
      Etchegoin, P. G.; Le Ru, E. C.; Meyer, M. Erratum: An analytic model for the optical properties of gold. J. Chem. Phys. 2006, 125, 164705,  DOI: 10.1063/1.2360270 .
      J. Chem. Phys. 2007, 127, 189901.
  • Supporting Information

    Supporting Information

    ARTICLE SECTIONS
    Jump To

    The Supporting Information is available free of charge at https://pubs.acs.org/doi/10.1021/acs.jctc.0c00269.

    • Additional results as shown in Figures S1–S4, Tables S1 and S2, and refs 1 and 2 (PDF)


    Terms & Conditions

    Most electronic Supporting Information files are available without a subscription to ACS Web Editions. Such files may be downloaded by article for research use (if there is a public use license linked to the relevant article, that license may permit other uses). Permission may be obtained from ACS for other uses through requests via the RightsLink permission system: http://pubs.acs.org/page/copyright/permissions.html.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

Pair your accounts.

Export articles to Mendeley

Get article recommendations from ACS based on references in your Mendeley library.

You’ve supercharged your research process with ACS and Mendeley!

STEP 1:
Click to create an ACS ID

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

Please note: If you switch to a different device, you may be asked to login again with only your ACS ID.

MENDELEY PAIRING EXPIRED
Your Mendeley pairing has expired. Please reconnect