Understanding How Acoustic Vibrations Modulate the Optical

United States. § Center for Nanoscale Materials, Argonne National Laboratory, Argonne, Illinois 60439, United States. ACS Nano , Article ASAP. DO...
9 downloads 19 Views 4MB Size
Subscriber access provided by UNIVERSITY OF ADELAIDE LIBRARIES

Article

Understanding How Acoustic Vibrations Modulate the Optical Response of Plasmonic Metal Nanoparticles Aftab Ahmed, Matthew Pelton, and Jeffrey R. Guest ACS Nano, Just Accepted Manuscript • DOI: 10.1021/acsnano.7b04789 • Publication Date (Web): 17 Aug 2017 Downloaded from http://pubs.acs.org on August 18, 2017

Just Accepted “Just Accepted” manuscripts have been peer-reviewed and accepted for publication. They are posted online prior to technical editing, formatting for publication and author proofing. The American Chemical Society provides “Just Accepted” as a free service to the research community to expedite the dissemination of scientific material as soon as possible after acceptance. “Just Accepted” manuscripts appear in full in PDF format accompanied by an HTML abstract. “Just Accepted” manuscripts have been fully peer reviewed, but should not be considered the official version of record. They are accessible to all readers and citable by the Digital Object Identifier (DOI®). “Just Accepted” is an optional service offered to authors. Therefore, the “Just Accepted” Web site may not include all articles that will be published in the journal. After a manuscript is technically edited and formatted, it will be removed from the “Just Accepted” Web site and published as an ASAP article. Note that technical editing may introduce minor changes to the manuscript text and/or graphics which could affect content, and all legal disclaimers and ethical guidelines that apply to the journal pertain. ACS cannot be held responsible for errors or consequences arising from the use of information contained in these “Just Accepted” manuscripts.

ACS Nano is published by the American Chemical Society. 1155 Sixteenth Street N.W., Washington, DC 20036 Published by American Chemical Society. Copyright © American Chemical Society. However, no copyright claim is made to original U.S. Government works, or works produced by employees of any Commonwealth realm Crown government in the course of their duties.

Page 1 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

TOC figure 82x44mm (300 x 300 DPI)

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Understanding How Acoustic Vibrations Modulate the Optical Response of Plasmonic Metal Nanoparticles Aftab Ahmed §,‡,*, Matthew Pelton†,‡, Jeffrey R. Guest‡,* §

Department of Electrical Engineering, California State University Long Beach, CA 90840, United States





Department of Physics, University of Maryland, Baltimore County, MD 21250, United States

Center for Nanoscale Materials, Argonne National Laboratory, Argonne, Illinois 60439, United States

* Address correspondence to [email protected] [email protected]

KEYWORDS: nanoparticle, acoustic vibrations, plasmons, transient absorption, phonons, optical response

1

ACS Paragon Plus Environment

Page 2 of 31

Page 3 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

ABSTRACT

Measurements of acoustic vibrations in nanoparticles provide an opportunity to study mechanical phenomena at nanometer length scales and picosecond time scales. Vibrations in noble-metal nanoparticles have attracted particular attention, because they couple to plasmon resonances in the nanoparticles, leading to strong modulation of optical absorption and scattering. There are three mechanisms that transduce the mechanical oscillations into changes in the plasmon resonance: (1) changes in the nanoparticle geometry, (2) changes in electron density due to changes in the nanoparticle volume, and (3) changes in the interband transition energies due to compression / expansion of the nanoparticle (deformation potential). These mechanisms have been studied in the past to explain the origin of the experimental signals; however, a thorough quantitative connection between the coupling of phonon and plasmon modes has not yet been made, and the separate contribution of each coupling mechanism has not yet been quantified. Here, we present a numerical method to quantitatively determine the coupling between vibrational and plasmon modes in noble-metal nanoparticles of arbitrary geometries, and apply it to silver and gold spheres, shells, rods, and cubes in the context of time-resolved measurements. We separately determine the parts of the optical response that are due to shape changes, changes in electron density, and changes in deformation potential. We further show that coupling is in general strongest when the regions of largest electric field (plasmon mode) and largest displacement (phonon mode) overlap. These results clarify reported experimental results, and should help guide future experiments and potential applications.

2

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 4 of 31

Measurements of acoustic vibrations in metal nanoparticles have revealed fundamental mechanical properties of the nanoparticles1, 2 and have been used as a probe of the mechanical properties of the local fluid environment, revealing phenomena such as the viscoelastic response of simple liquids at high frequencies.3, 4 Monitoring the vibrational frequency of the nanoparticles as molecules adsorb on their surface has the potential to enable mass-based sensing in a variety of environments with high sensitivity,5,

6

complementing other metal-nanoparticle chemical

sensing techniques based on enhanced Raman scattering7 or plasmon frequency shifts.8-10 The majority of experimental measurements of nanoparticle vibrations have employed noblemetal nanoparticles, because the optical spectra of these particles exhibit strong plasmon resonances whose frequencies depend on the sizes and shapes of the nanoparticles.11 Femtosecond transient absorption (TA) spectroscopy can thus be employed as a sensitive, noncontact means of detecting the transient phonon oscillations. Energy from an ultrafast pump laser pulse is absorbed by electrons in the nanoparticle and transferred to the lattice within a few picoseconds, resulting in a rapid thermal expansion of the particle. This expansion, in turn, excites coherent acoustic oscillations, or low-frequency phonon modes, that result in a periodic variation in shape and size of the nanoparticle. The strong optical absorption of the plasmonic metal nanoparticles and the rapid heating of the lattice enable efficient excitation of the acoustic modes; high-sensitivity detection of the vibrations is enabled by the sensitivity of the plasmon resonances to changes in the size and shape of the particle, resulting in modulation of the transmission of a probe laser pulse. There are three mechanisms that transduce the mechanical oscillations into changes in the plasmon resonance frequency: (1) changes in shape during oscillation (shape effect); (2) changes in electron density (ED) due to changes in nanoparticle volume (volume effect); and (3) changes in inter-band transition energies through the deformation potential (DP). Experimental work has mainly focused on the frequency shift1, 12-14 and damping15-17 of the vibrations, while vibrational amplitudes have received less attention. Spherical nanoparticles have been treated in detail

12, 13, 18

due to the availability of analytical

solutions for these geometries. The vibrational modes of spherical particles were first computed by Lamb,19 who found that the period of the symmetric breathing mode is proportional to the radius of the particle. This has been experimentally verified for gold nanoparticles in water with sizes ranging from 8 nm to 120 nm.20 Crut et al. presented a semianalytical method to study the

3

ACS Paragon Plus Environment

Page 5 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

vibrational modes of core-shell spherical nanoparticles.21 The contribution of various acoustoplasmon coupling mechanisms, particularly DP, to Raman scattering has been studied analytically for spherical nanoparticles.18 Recently, mechanical coupling in gold nanoparticles was studied, revealing ultralow frequency Raman scattering from a dimer configuration when excited near the plasmon resonance frequency.22 Finite-element-method calculations have been used to determine the vibrational modes in particles without spherical symmetry.23

In most cases, only mechanical simulations were

performed, and plasmon-vibration coupling strengths and mechanisms were not considered. Only a limited number of studies reported both optical and mechanical calculations; however, these studies – which revealed the role of “hot spots” and focused primarily on Raman scattering - either dealt with only one geometry or only one coupling mechanism.24,

25

Raman and TA

measurements provide complementary information, as they involve detection of different vibrational modes. As the origins of TA signals remain relatively unexplored, we focus here on acousto-plasmon coupling in TA experiments. Specifically, we perform numerical calculations to quantitatively determine the coupling between vibrational and plasmon modes in noble-metal nanoparticles. We integrate finiteelement (FE) calculations of the mechanical response with finite-difference time-domain (FDTD) calculations of the optical response in order to determine the optical signals that are expected for different vibrational modes. We focus on four commonly investigated nanoparticle geometries: spheres,1, 13, 16, 20, 26-28 rods,16, 17, 29-33 shells,34, 35 and cubes.1, 36-38 We separately determine ,

the parts of the optical response that are due to shape changes, volume changes (changes in local ED), and changes in interband transition energies (DP). Although we present results that are relevant to time-resolved measurements, this method can be readily extended to include asymmetric vibrational modes that are Raman active.

These results illuminate reported

experimental results, and will help guide the optimization of future experiments and guide design for potential applications.

4

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

CALCULATIONS We calculate the frequencies and mode shapes of acoustic vibrations using FE simulations, using the bulk mechanical properties of silver and gold. The initial condition for the oscillations is taken to be an isotropic expansion of the particle, corresponding to uniform heating by an incident laser pulse. Deformed geometries are then determined corresponding to the maximum compression and maximum expansions for the various vibrational modes. Finally, the optical response of the deformed geometries are calculated using the FDTD method. We separately calculate the contributions to the overall shift in plasmon resonance due to shape deformation, local changes in ED, and local DP.

Calculation details regarding simulating deformed

geometries and calculation of permittivity based on local effects are given in the Methods section, with further details in the Supporting Information (SI). Experimentally, the initial expansion of the particle is proportional to the energy in the pump pulse, with greater expansions resulting in greater amplitudes of oscillation for each of the normal modes. It is expected that a larger amplitude of deformation will result in a larger plasmon shift. However, a quantitative relation showing the dependence of plasmon shift on amplitude of deformation has not been reported. Here, we investigate this relation by calculating the optical response corresponding to deformation amplitudes of 1 nm, 0.5 nm, and 0.25 nm for the cases of gold and silver nanospheres as well as gold nanorods. Our calculations show a linear relation between plasmon shift and deformation amplitude, as discussed below. This allows for the use of an arbitrary deformation amplitude in the FDTD calculations; the plasmon shift corresponding to any given vibrational amplitude can then be calculated by linear interpolation of the simulation results.

Conversely, the vibrational amplitudes in a TA

experiment can be estimated by analyzing the observed shift in plasmon frequency. In our calculations, we use an amplitude of 1 nm for the dominant vibrational mode, allowing the calculations to be performed without requiring very small mesh sizes, and thus reducing the required computational resources. It should be noted that the phonon modes are normalized to the amplitude of the dominant mode, and thus the amplitudes of the remaining modes are smaller than 1 nm.

5

ACS Paragon Plus Environment

Page 6 of 31

Page 7 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

NANOSPHERES Vibrational Modes of Nanospheres Spherical metal nanoparticles have relatively simple plasmonic and vibrational mode structures that can be solved analytically due to their high degree of symmetry, and thus serve as a useful test case to validate our numerical simulations. Vibrations of an elastic, isotropic sphere can be calculated by solving Navier’s equation. The eigenmode solutions correspond to the series of spherical harmonics; however, for the uniform initial stress that results from isotropic heating, only spheroidal modes with zero angular momentum are excited.27 The displacement vector 𝑢 is given by: 𝑢 = 𝑎! 𝑑 𝑑𝑟 [𝑗! 𝜔! 𝑟 𝑐! ], where 𝑗! is the spherical Bessel function of the first kind, 𝑎! is the radial unit vector, and cl is the longitudinal speed of sound in the particle.39, 40 The eigenfrequencies 𝜔! can be approximated analytically as 𝜔! 𝐿, 𝑛 = 𝜒! 𝑐! 𝑅, where R is the radius of the sphere and 𝜒! is the nth solution of: 𝜒! cot 𝜒! = 1 − (𝜒! 𝑐! /2𝑐! )! , where 𝑐! is the transverse speed of sound in the particle.39 The mode shapes of the fundamental breathing mode (𝑛 = 0) and the first higher-order radial mode (𝑛 = 1) of a 60-nm silver nanosphere are shown in the inset of Fig. 1a. Figure 1b compares the analytical (solid line) and numerical results (circles) for the absolute displacement as a function of position within the sphere for these two modes (fundamental breathing mode normalized to 1 nm); the results are in excellent agreement. The above equations predicts resonance frequencies of 55 GHz and 119 GHz for the two modes of a 60-nm silver sphere, in good agreement with the numerically calculated values of 56 GHz and 119 GHz, respectively. The analytical solutions thus validate the accuracy of our numerical calculations. Optical Response of Vibrating Nanospheres The optical response of a spherical metal nanoparticle is dominated by the lowest-order, dipolar plasmon mode. In the quasi-static limit, when the sphere diameter is small compared to the optical wavelength, the frequency of this mode is given by the condition11 𝜖! 𝜔 + 2𝜖! = 0,

(1)

where 𝜖! and 𝜖! are the dielectric functions of the metal and the surroundings, respectively. The frequency of the plasmon mode is independent of the size of the particle in this limit, so the

6

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

radial modes that are excited by heating are not expected to result in a shift in the plasmon frequency due to the shape effect only. We have verified this by solving Maxwell’s equations directly in the FDTD calculations, without employing the quasi-static approximation. The extinction cross-section of the particle, on the other hand, is sensitive to the particle volume;11 we see a modulation in the extinction cross-section (with negligible peak shift) as the particle becomes larger and smaller. However, the increase in cross-section from the size change is almost completely cancelled by ED and DP effects. The dominant overall effect of the nanoparticle vibration is thus a shift of the frequency of the plasmon resonance, as seen in Fig. 1c, with the largest contribution coming from DP. The relation between the optical response and deformation amplitude is investigated by calculating the plasmon shift of the fundamental breathing mode of silver and gold nanospheres for varying deformation amplitudes, considering all three coupling mechanisms. Results for a silver nanosphere are shown in Fig. 2; results for a gold nanosphere are given in Fig. S4 of the SI. A simple interpolation algorithm is discussed in the SI along with an example showing the interpolated optical response corresponding to vibrational amplitudes of 125 pm and 250 pm (see Fig. S5). The higher-order radial mode results in a small change in the volume of the sphere, reflected by a surface displacement of about one fifth of the fundamental mode (Fig. 1a), and experimental evidence for this mode has been obtained only through detailed analysis of transient-absorption data.26 Further details about the optical response of nanospheres and other geometries studied in this work are provided in the SI, where the contributions from all three coupling mechanisms are separately presented; a summary is presented in Fig. 7.

7

ACS Paragon Plus Environment

Page 8 of 31

Page 9 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Figure 1. Calculation results for a 60-nm silver sphere. (a) Absolute surface displacement as a function of frequency. Isotropic radial vibrational modes (fundamental and higherorder breathing modes) are shown in the inset. The mode shape shows the absolute local displacement of the deformed nanoparticle from the equilibrium position. The vibrational amplitudes are normalized to reduce computational resources, and are large compared to experimentally accessible values (see text), however, linear acoustoplasmon coupling allows for the estimation of plasmon shift in the limit of small deformations. (b) Absolute normalized displacement as a function of position within the sphere for the two modes. Solid lines represent the theoretical results, and circles and squares are the simulation results for the fundamental and higher-order modes respectively. (c) and (d) Net change in extinction cross-section, including combined contributions from shape, electron-density (ED), and deformation-potential (DP) mechanisms, for the fundamental (c) and higher order breathing modes (d). The relative contribution of each coupling mechanism is presented in the Supporting Information. Extinction spectra for expanded, non-deformed, and contracted nanospheres are shown in the lower panels of (c, d) by dashed-dot blue, solid black, and dashed red curves, respectively. The change in extinction cross-section (𝛿𝜎!"# ) is plotted in the upper panels. The acoustic modes in the insets of (c) and (d) use the same color scales as in (a). A TA experiment measures the difference between the spectrum of the deformed nanoparticle and the un-deformed nanoparticle as a function of time as the particle undergoes acoustic vibrations. A TA measurement on a vibrating silver nanosphere is thus predicted to show strong modulation when the probe wavelength is detuned half a linewidth away from the initial plasmon 8

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

resonance in either direction, as shown in Fig. 1c. The modulation on the blue side of the resonance is 180 degrees out of phase with the modulation on the red side of the resonance, and the modulation passes through zero near the plasmon resonance frequency. Changes in linewidth result from the different optical absorptivity of silver at the plasmon resonance as the plasmon resonance frequency oscillates, and are predicted to be small. All of these numerical predictions are consistent with experimental results.27

Figure 2. Calculation results for a 60-nm silver sphere showing plasmon resonance shift for varying phonon vibrational amplitudes. (a) Changes in extinction cross-section for vibrational amplitude of 1 nm (green dots), 0.5 nm (dashed-dot red), 0.25 nm (dashed blue) and 0 nm (solid black). (b) Plasmon resonance wavelength as a function of deformation amplitude. Blue squares represent the calculated values and the dashed red line is a linear fit. The mechanical properties of silver and gold are very similar, resulting in nearly identical vibrational modes, as shown in Fig. 3a. On the other hand, gold has stronger optical absorption than silver at wavelengths around 500 nm, corresponding to the blue edge of the sphericalnanoparticle plasmon resonance. The result is that the extinction cross-section of the particle undergoes significant modulation as the plasmon resonance passes into and out of wavelength regions where gold is strongly absorbing. The calculated change in extinction is significant on the red side of the plasmon resonance, where optical absorption of gold is relatively weak, but is 9

ACS Paragon Plus Environment

Page 10 of 31

Page 11 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

much smaller on the blue side of the resonance, where the absorption of gold is stronger; this is shown in Fig. 3b and is consistent with experimental observations.41 The higher order mode again does not induce large changes in the extinction cross section (Fig. 3c).

Figure 3. Calculation results for a 60-nm gold sphere. (a) Displacement of a point on the surface as a function of frequency. Isotropic radial modes (fundamental and higher order breathing modes) are shown in the inset. The mode shape shows the absolute local displacement of the deformed nanoparticle from the equilibrium position. (b-c) Changes in extinction cross-section for the two vibrational modes.

10

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

NANORODS Vibrational Modes of Nanorods While the high symmetry of nanospheres means that the vibrational and plasmon mode structures of the particles and the coupling between vibrations and plasmons are relatively straightforward, cylindrically-symmetric nanorods manifest a more complex coupling.42 Specifically, the two lowest-order vibrational modes of a nanorod are the fundamental extensional mode, which primarily affects the aspect ratio (length / width) of the rod, and a higher-frequency mode that primarily affects the nanorod volume. Fig. 4a illustrates these two modes for a 55 nm long and 26 nm wide gold nanorod. The calculated mode frequencies of 19 GHz and 102 GHz agree well with the theoretical predictions of Ref. 29, which suggests 18.4 GHz and 96 GHz for the extensional and breathing modes, respectively. While our calculations suggest that the lowest-order mode is indeed primarily an extensional mode, the shape of the higher frequency mode does not correspond exclusively to a breathing mode; instead it is a combination of the fundamental (n = 0) breathing and the higher order (n = 2) extensional mode, due to the proximity of their resonant frequencies.29 For the small aspect ratio of 2.1 shown in Fig. 4a, the high-frequency acoustic mode is excited with an amplitude comparable to that of the extensional mode. However, as shown in Fig. 4d, this high-frequency mode gets weaker as the aspect ratio increases, consistent with experimental findings.29

11

ACS Paragon Plus Environment

Page 12 of 31

Page 13 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Figure 4. Calculation results for a 55-nm long and 26-nm wide gold nanorod (left panels) and 61-nm long and 22-nm wide nanorod (right panels). (a, d) Displacement of the rod end (solid green) and of the mid-point (dashed magenta) of the nanorod as a function of frequency. Mode shapes are also shown along with a circle and arrow indicating the point on the rod’s surface where displacement is measured. (b, e) Changes in extinction cross-section in response to the extensional vibrational mode. (c, f) Changes in extinction cross-section in response to the higher order vibrational mode. Optical Response of Vibrating Nanorods Similar to the vibrational response, the plasmonic response of the nanorod splits into a lowerfrequency (long-wavelength) longitudinal mode and a higher-frequency (short-wavelength) transverse mode. For the longitudinal plasmon, the spatial mode overlaps well with the regions of largest displacement for the extensional acoustic mode, resulting in efficient modulation of the plasmon by this vibration, as shown in Fig. 4b and 4e. This interplay is primarily due to changes in the aspect ratio of the nanorod. The longitudinal plasmon mode can be thought of as a Fabry12

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 14 of 31

Perot (FP) resonance for surface-plasmon polaritons propagating along a metal nanowire and reflected from the two ends of the rod. The resulting standing-wave plasmon mode, with a node at the center of the rod, is sensitive to changes in length; the extensional vibrational mode thus produces a significant shift in the longitudinal plasmon resonance. As for gold and silver nanospheres, the plasmon shift as a function of the amplitude of the extensional mode of the gold nanorod shows a linear relation, as shown in the SI (Fig. S6). The effect of the oscillating aspect ratio is a periodic modulation of the plasmon resonance frequency, as it is for the breathing modes of nanospheres. However, the longitudinal plasmon frequency occurs at wavelengths where optical absorption in the bulk metal is weak; the optical response is thus predicted to be symmetric (but 180 degrees out of phase) for the red and blue sides of the plasmon resonance, as observed experimentally.29

Since the metal absorption

changes little as the plasmon resonance changes, the linewidth of the plasmon is predicted to be virtually unaffected by the acoustic oscillations; this is again consistent with experiment.30 For the higher-frequency acoustic mode, expansion of the nanorod reduces its aspect ratio, resulting in a blue shift of the plasmon resonance when only shape effects are considered. However, the expansion also significantly increases the volume of the nanorod, modifying the ED and the DP, and thereby shifting the plasmon resonance towards the red (Fig. S9, S10 and 4 c,f). The overall optical signal is a balance between these competing effects, and is thus sensitive to the size and aspect ratio of the nanoparticles. A non-negligible modulation of the longitudinal plasmon resonance due to the breathing mode is predicted for gold nanorods with lengths of 55 nm and widths of 26 nm.29 Such a modulation has been observed, together with the larger modulation due to the extensional mode, in single-particle transient-absorption measurements on gold nanorods.16,

17, 31

The corresponding signal was absent in ensemble

measurements,30 but this may have been a result of rapid dephasing of the breathing-mode signal due to inhomogeneities in the sample. Both the extensional and breathing acoustic modes also modulate the transverse plasmon mode of the gold nanorod, as shown in Fig. S9 and S10. The effect is negligible, however, for the case of gold, due to strong metal absorption at the transverse plasmon wavelength. In the case of silver, coupling of the extensional vibrational mode to the transverse plasmon mode is calculated

13

ACS Paragon Plus Environment

Page 15 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

to be an order of magnitude weaker than its coupling to the longitudinal plasmon mode. The breathing mode of the silver rod couples to the transverse plasmon mode with about 1/3 of the strength of its coupling to the longitudinal plasmon mode. Optical measurements to date have therefore focused nearly exclusively on the longitudinal plasmon mode. For comparison, Fig. S11 shows vibrational as well as optical results for a silver nanorod; because of the large spectral gap between the inter-band transition frequencies and plasmon resonance of silver nanorods, DP effects are negligible. NANOSHELLS Vibrational Modes of Nanoshells The acoustic vibrational modes of a nanoshell can be understood in terms of coupling between motion of the inner and outer surfaces. The vibrational frequencies thus depend on the outer radius, Rout, and on the radius ratio (Rin/Rout), where Rin is the inner radius. Fig. 5 shows calculation results for the two lowest-order vibrational modes for two cases from Ref. 35: (i) Rout = 40 nm, radius ratio of 0.75 (Fig. 5a), (ii) Rout = 25 nm, radius ratio of 0.6 (Fig. 5d). In both cases, the initial isotropic stress in the nanoparticle predominantly excites the lowest-order, symmetric mode, where the inner and outer surfaces of the shell are displaced in phase. The higher frequency, antisymmetric vibrational mode is excited to a lesser degree, corresponding to out-of-phase displacement of the two surfaces. The period of the symmetric vibrational mode is longer than that of a solid sphere of the same size and increases with decreasing shell thickness.34 Optical Response of Vibrating Nanoshells The plasmon modes on the inner and outer surfaces of the shell couple electrostatically, and the resulting optical response is dominated by the symmetric coupled dipolar mode.43 The plasmon resonance frequency shifts as the thickness of the shell varies, providing a widely tunable plasmonic response. The response of the plasmon resonances in nanoshells to the vibrational modes is shown in Fig. 5. In both cases, the fundamental breathing (symmetric) mode produces a strong oscillation in plasmon resonance frequency.

However, the antisymmetric vibrational mode produces a 14

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

significant plasmon oscillation only in the case of the smaller shell (Fig. 5f); this is consistent with experimental results.35 This occurs for two reasons. Firstly, the amplitude of displacement is larger for the smaller nanoshell; secondly, in the case of the smaller nanoshell, all three coupling mechanisms produce plasmon shifts towards the red, resulting in a stronger modulation of the probe beam.

Figure 5. Calculation results for a gold nanoshell with inner radius of 30 nm and outer radius of 40 nm (left panels), and a nanoshell with inner radius of 15 nm and outer radius of 25 nm (right panels). (a, d) Radial displacement of the inner surface (solid green) and outer surface (dashed magenta) of the nanoshell as a function of frequency. The isotropic radial phonon modes are also shown along with a circle and arrow indicating the points where displacement is measured. (b, e) Changes in extinction cross-section in response to the fundamental and (c, f) higher order breathing modes.

15

ACS Paragon Plus Environment

Page 16 of 31

Page 17 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

As the hollow shell expands and contracts, the inner and outer surfaces move by different amounts, changing the radius ratio of the shell. For the smaller nanoshell, the solid shell occupies a larger fraction of the interior of the particle, so that the fractional change in volume of the solid portion is larger than for the larger nanoshell. This means that the smaller nanoshell experiences a larger change in ED and also larger effects due to DP, leading to an additional shift of the plasmon resonance frequency. The fact that shape changes dominate the optical response is similar to the case of nanorods, and is in contrast to the case of solid nanospheres, where the dominant contribution to the optical signal comes from changes in ED and DP. A similarity between the response of the solid sphere and the smaller nanoshell is that the metal absorption is strong on the blue side of the plasmon resonance, so that a more significant transient signal is expected on the red side of the resonance (Fig. 5b and 5e), consistent with experimental observations.35 For the higher-order radial mode, the relative displacements of the inner and outer surfaces is sensitive to radius ratio. For the larger particle of the two that we simulated, this mode happens to result in nearly equal displacements of the two surfaces. This means a minimal change in the shell thickness, ED, and DP, resulting in a small change in the plasmon resonance (Fig. 5c). For the smaller particle, the relative motion of the two surfaces is larger, leading to a larger optical signal (Fig. 5f). This is again consistent with the experimental results.35 NANOCUBES Vibrational Modes of Nanocubes Nanocubes have lower symmetry than either nanorods or nanoshells, and exhibit correspondingly more complex optical and mechanical behavior. Of the many acoustic modes that are possible, two have been predicted to dominate the mechanical response.38 Fig. 6a shows the calculated deformations corresponding to these two lowest-order acoustic modes for a silver nanocube with 35-nm edge length. These results show that the lowest-order mode primarily involves deformation of the eight corners of the cube, whereas the higher-order mode primarily involves displacement of the centers of the six faces of the cube.

16

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Optical Response of Vibrating Nanocubes Just as only two vibrational modes dominate the mechanical response of a nanocube, only six dominant plasmon modes account for more than 90% of the optical absorption within the visible spectral range.44, 45 These dominant plasmon modes are sensitive primarily to the shape of the corners of the cubes, where the electric fields are strongly concentrated.45 They are therefore expected to couple strongly to the fundamental vibrational mode, which involves significant deformation of the corners, and only weakly to the higher-order mode, which does not. This is confirmed by the calculated changes in optical spectra shown in Fig. 6 b-c. Even though the mechanical amplitude of the higher-order mode is nearly as large as that of the fundamental mode, it results in a much smaller shift in the plasmon peak. This is consistent with TA measurements on single nanocubes, which have shown only a single vibrational frequency for the majority of particles studied.36 A similar importance of field localization has previously been recognized in the context of Raman scattering from metal nanoparticles.24, 25

17

ACS Paragon Plus Environment

Page 18 of 31

Page 19 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Figure 6. Calculation results for a 35-nm silver nanocube. (a) Displacement as a function of frequency at the corner (solid green) and face-center (dashed magenta) of the nanocube. The fundamental and overtone breathing modes are also shown. (b-c) Changes in extinction cross-section for the two modes. The signal produced by the higher-order mode is further suppressed by the opposing effects of shape and volume changes. For the fundamental vibrational mode, expansion of the cube produces an outward motion of the corners that results in a red shift of the plasmon mode; this adds to the red shift produced by the decreasing ED and a further slight red shift by DP, resulting in a larger overall signal.

For the higher-order mode, by contrast, expansion of the cube

produces an inward motion of the corners and outward motion of the faces, causing a blue shift of the plasmon mode; this works against the red shift produced by the decreasing ED and DP, resulting in a smaller overall signal.

18

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

The magnitude and phase of the signal produced by the higher-order mode will depend on the balance between the shape and volume effects, which will, in turn, be determined by the size of the cube. Fig. S15 shows results for a larger cube, with edge length of 85 nm. Strong optical scattering from the 85 nm nanocube significantly modifies the extinction spectra as compared to the smaller nanocube, for which scattering is relatively weak. In this case, the optical signal due to the higher-order mode is significantly stronger than for the smaller cube, and is comparable to the optical signal due to the fundamental mode. This is consistent with experimental results on ensembles of nanocubes, where smaller cubes have shown signals due only to the fundamental vibrational frequency, whereas cubes with edge lengths larger than 50 nm have shown signals that correspond to both modes.38 DISCUSSION Table 1 presents a summary of our numerical results along with a comparison to theoretical predictions and to experimentally determined values where available. Differences between model and experiment are largely due to the effects of the surrounding medium, which are not included in our mechanical calculations. For most of the particles, the effect of the surroundings is primarily to damp the vibrations.16 This results in a relatively small shift in resonance frequency, so that our calculations are in reasonable agreement with experiment.3 The nanoshells are the exception: motion of the inner surface of real colloidal nanoshells is restricted by the nearly incompressible fluid filling the shell, resulting in significantly higher vibrational frequencies. Accounting for this effect quantitatively will require coupling the finite-element treatment of the nanoparticle vibration, based on the Navier equation, to a finite-element modeling of the enclosed liquid, based on the Navier-Stokes equation; such a model is beyond the scope of this paper. Nonetheless, the full fluid-dynamic treatment is not expected to have a significant effect on the geometry of the vibrational modes, so the coupling between vibrational and plasmonic modes should be mostly unchanged by these effects.

19

ACS Paragon Plus Environment

Page 20 of 31

Page 21 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Table 1. Summary of the numerical results for different geometries investigated in this work.

Shape

Material

Dimensions (nm)

Silver

Radius = 30

Phonon Modes (GHz) Expt. Sim. (this work) 55

Gold

Gold

Silver

Length = 55 Width = 26

50.7

54

--

113

16.1

19

91

102

15.3

17

[29]

Length = 61 Width = 22

[29, 46]

Rin = 30 Rout = 40

87.7

117

15

28 [35]

Rin = 15 Rout = 25

--

135

37.4

45 [35]

81

171

52.6

57

--

80

16.1

23

Edge = 35

[38]

Cube Silver

122 [13, 20]

Shell Gold

-Radius = 30

Rod Gold

56 [27]

Sphere Gold

Expt. Ref.

Edge = 85

[38] 23.8

20

ACS Paragon Plus Environment

33

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Figure 7. Plasmon resonance shift corresponding to normalized vibrational modes for metal nanoparticles with different geometries (amplitude of the dominant vibrational mode is 1 nm). Contributions to the overall shift as the particle expands are indicated from shape only (red circle), shape + electron density (blue square), and shape + electron density + deformation potential (green diamond). L and H represent the low-frequency and high-frequency vibrational modes, respectively. The dashed vertical line represents the plasmon resonance of the undeformed nanoparticle. Figure 7 summarizes the calculated shifts in plasmon resonance frequency for each vibrational mode of the various nanoparticles due to shape changes, ED, and DP. In the case of spherical nanoparticles, the fundamental breathing mode is the dominant mode, with surface displacement about 5 times larger than that of the next higher order breathing mode (see Fig. 1 and 3), in

21

ACS Paragon Plus Environment

Page 22 of 31

Page 23 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

agreement with experimental results.16, 20, 27, 47, 48 The higher-order breathing mode is excited 180 degrees out of phase with respect to the fundamental mode, which would suggest a blue shift of the plasmon resonance; however, the local volume change involves both expansion and compression, as seen in the inset of Fig. 1a, modulating the local permittivity in both directions from the equilibrium position. The weighted average permittivity is thus slightly modified from the equilibrium position, producing a negligible red shift in plasmon resonance for a 60 nm sphere.

In terms of coupling mechanism, DP plays a significant role in the case of gold

nanoparticles as compared to silver because of the spectral location of the inter-band frequencies. The magnitude of expansion in experiments can be estimated using the absorbed pump energy and the corresponding change in temperature, ∆𝑇 = 𝑄/𝐶𝑉, where Q is the absorbed energy, and C and V are the heat capacity and volume of the nanoparticle, respectively. Q is calculated as the product of the energy density in the pump pulse and the extinction cross-section of the nanoparticle. Finally, the change in volume is calculated using the volumetric thermal expansion of the particle. Using typical values of pump pulse energy (100 nJ) and the calculated extinction cross-section of a 60 nm gold sphere (1.4 x 10-14 m2), we estimate a 56 pm change in the radius of the particle. Considering the demonstrated linearity of the plasmon shift (Fig. S4b), this corresponds to a plasmon shift of about 1.23 nm, using the results of Fig. 3 and Fig. 7. In the case of nanorods, we can perform the reverse calculation, and estimate the actual deformations of nanoparticles by comparing transient absorption measurements to the numerically calculated change in plasmon resonance wavelength (δλr).

Ruijgrok et al.

performed pump-probe measurements on single gold nanorods with an aspect ratio of 2.4 and length of 60 nm.16 From the pump-probe data and using the results of Fig. 4e, we estimate a deformation of about 0.3 nm. We can also estimate the expected change in the length of the nanorod using ∆𝐿 = 𝐿𝛼∆𝑇, where 𝛼 is the thermal expansion coefficient and L is the initial length of the rod. For the same experiments, this results in Δ𝐿 = 0.299 nm, in good agreement with our spectroscopic determination. In the case of nanoshells, a plasmon shift of approximately 2.5 nm has been reported in response to the fundamental breathing mode of a nanoshell with an outer radius of 40 nm.35 Comparing the experimental plasmon resonance shift to our numerical results, the experimental plasmon shift

22

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 24 of 31

corresponds to a 0.25 nm change in inner radius of the shell and a 0.168 nm change in the outer radius of the shell. On the other hand, using the absorbed pump energy and thermal expansion calculations an approximate change of 1.4 nm in the inner radius and 0.97 nm in the outer radius is predicted. As this calculation ignores the heat transferred from the metal shell to the liquid within the shell, it is likely that the temperature rise is lower than predicted and the smaller amplitudes of deformation determined spectroscopically based on our calculations are more accurate. In the case of nanocubes, our results predict the detection of a single vibrational mode (the fundamental breathing mode) for smaller nanocubes with an edge length of 35 nm. The higherorder breathing mode produces negligible plasmon shift, because the small shift in plasmon resonance due to shape modulation is almost completely cancelled by ED effects, as seen in Fig. S14b. On the other hand, for larger nanocubes with an edge length of 85 nm, optical scattering becomes significant and produces a broad feature on the red side of the resonance. In this case, our results predict the detection of two vibrational modes, because the two modes produce comparable plasmons shifts 180 degree out of phase with one another, as shown in Fig. S15. These results are again in agreement with experimental findings.38 CONCLUSIONS We have developed a numerical approach that predicts the optical response due to the mechanical vibrations of plasmonic metal nanoparticles. The method is applicable to arbitrary geometries, and provides insight into the mechanisms responsible for coupling mechanical vibrations to plasmon resonances. This allows for studying additional variables in experimental measurements of the nanoparticle vibrations: instead of being limited to analysis of only vibrational frequencies and damping rates, we are now able to quantitatively relate optical signals to mode shapes and vibrational amplitudes.

Since the calculations predict the full

spectral response due to the vibrating nanoparticles, they can be used to optimize experimental measurements by predicting the probe wavelengths that will be most sensitive to specific vibrational modes. Moreover, they can be used to optimize nanoparticle geometries in order to maximize the optical signal due to a given vibrational amplitude, with potential application to ultrasensitive mass detection.6

23

ACS Paragon Plus Environment

Page 25 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

Future work will be directed towards applying the method to additional nanoparticle geometries with lower symmetry that have been studied experimentally, such as bipyramids,3, 16 nanodisks,49 crosses,50 and core-shell nanoparticles.5, 51 In addition, the proposed method can be applied to study vibrational modes that are accessed using Raman spectroscopy and compare the results with experimental observations.18, 24, 25 Further connection with experiment can be achieved by coupling the Navier-equation treatment of nanoparticle vibrations to a Navier-Stokes treatment of surrounding liquid in order to quantify the effects of the nanoparticle environment on the signal transduction. Our method can be extended to separately and quantitatively study the effects of mechanical vibrations on the real and imaginary parts of the dielectric function of a single nanoparticle, enabling comparison to time-resolved interferometric measurements.28, 52 METHODS The vibrational modes of the nanoparticles are calculated by solving the Navier equation with continuity of stress and displacement at the boundary of the nanoparticle. Impulsively excited modes are calculated in the frequency domain using the finite-element (FE) solver COMSOL Multiphysics.50,

53

The metal nanoparticles are assumed to be elastically isotropic and

homogeneous, bulk elastic constants are used, and the mechanical response is modeled using continuum mechanics.54, 55

The COMSOL built-in material database is used to represent gold

and silver using their bulk elastic constants and sound velocities. We ignore the mechanical properties of the surrounding medium and perform the FE calculations for metal nanoparticles in vacuum.

In solution-based experiments, mechanical interaction with the surrounding

environment can slightly shift the vibrational frequencies and damp the oscillations,56 though in many cases the frequency shifts are negligible;17 an isotropic loss factor has been introduced in the FE calculations to avoid arbitrarily large deformations. Small changes in the vibration frequency and damping due to the environment will not affect conclusions for most of the cases studied in this work, including spheres, rods and cubes; for nanoshells, however, the calculated mode frequencies and deformations may be affected by the core material. In addition, we do not take into account changes in the plasmon resonance due to heating of the nanoparticle or its environment. Provided that the change in temperature of the particle and its vibrational amplitude are both small – a condition that holds in most experiments – then the effect of heating on the plasmon resonance will be nearly independent of the effect of vibration, 24

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

Page 26 of 31

and the two can be separated by appropriate fitting of experimental data.16 We thus focus here exclusively on the coupling between mechanical vibrations and plasmon resonances. Laser-induced heating is modeled by a uniform initial raised temperature, which is accurate for particles larger than a few nanometers, so that the initial heat transfer from the electrons to the lattice is much faster than the vibrational period of the nanoparticle. The dielectric response of a non-deformed nanoparticle is modeled using the Lorentz-Drude (LD) model: !

𝜀! 𝜔 = 𝜀! ∞ + 𝜔!!

𝜔! !!! !,!

𝑓!  ,                                                                                         − 𝜔 ! − 𝑗𝜔Γ!

where 𝜔! is the plasma frequency, 𝑓! is the oscillator strength, 𝜔!,! is the inter-band transition frequency, and Γ! is the linewidth of the mth Lorentzian term. These parameters are obtained for gold and silver by fitting LD model to the experimental data of Refs. 57, 58, and are given in Table 1 of the SI. The bulk plasma frequency 𝜔! =

𝑁𝑒 !

𝜀! 𝑚    ,   is modified by changes in the free-electron

density, 𝑁, of a deformed nanoparticle: 𝑁 = 𝑁! (1 + ∆𝑉/𝑉! ), where e and m are the charge and effective mass of the electron, respectively, 𝑁! and 𝑉! are the free electron density and volume of the non-deformed nanoparticle, respectively, and ∆𝑉 is the change in volume. The local volume change is obtained from the finite-element calculations. Aside from modifying the ED, the deformation of the particle also modifies the inter-band transition frequencies via the DP. The shifts in 𝜔!,! are calculated using the DPs for gold and silver: 𝛿𝐸!,! = 𝜉! 𝑑𝑉 𝑉 , where 𝜉! and 𝛿𝐸!,! are the DPs and change in energy of the mth transition, and the inter-band transition frequency and energy are related by 𝐸!,! = ℏ𝜔!,! (Refs. 59, 60). Details of the calculation of the dielectric constant are provided in the SI. Equations S3 and S4 of the SI describe two possible ways (local and weighted average respectively) of calculating the optical response of a deformed nanoparticle, based on local or weighted-average dielectric constants. Fig. S2 and S3 of the SI shows that both methods

25

ACS Paragon Plus Environment

Page 27 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

produce almost identical results, and thus Equation S4 can be employed to reduce the required computational resources as well as the simulation time. Extinction cross-sections of the particles are calculated using Lumerical FDTD Solutions.61 The deformed geometry is transferred from COMSOL to Lumerical for calculations of optical properties in the time domain; details are provided in the SI. Time-domain results are Fourier transformed to produce plasmon resonance spectra as a function of frequency. The FDTD simulation domain is terminated by a perfectly matched layer. A mesh override region is defined for accurate representation of the nanoparticle, with the size of the mesh determined according to the magnitude of deformations that need to be modeled. A minimum mesh size of 0.15 nm was used, and conformal variant 1 mesh was used. The polarization of the incident light determines which plasmon modes are excited. In the case where the particle has spherical symmetry, the polarization is inconsequential. For the case of rods and cubes, we employed a source with electric field at 45 degrees with respect to the long axis of the rod or the edge of the cube. ACKNOWLEDGMENTS Use of the Center for Nanoscale Materials, an Office of Science user facility, was supported by the U. S. Department of Energy, Office of Science, Office of Basic Energy Sciences, under Contract No. DE-AC02-06CH11357. This material is based upon work supported by Laboratory Directed Research and Development (LDRD) funding from Argonne National Laboratory, provided by the Director, Office of Science, of the U.S. Department of Energy under Contract No. DE-AC02-06CH11357. This material is based in part on work supported by the National Science Foundation under Grant No. DMR-1554895. The authors thank Prof. John Sader and Dr. Debadi Chakraborty for helpful discussions. SUPPORTING INFORMATION AVAILABLE: This document includes details of calculations of the dielectric response of the deformed nanoparticles, local vs. averaged effects, methods used to transform deformed geometries from COMSOL to Lumerical, interpolation of numerical results, and calculated extinction crosssections for each nanoparticle showing the contribution of each coupling mechanism separately. This material is available free of charge via the Internet at http://pubs.acs.org.

26

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

REFERENCES 1. Hartland, G. V. Optical Studies of Dynamics in Noble Metal Nanostructures. Chem. Rev. 2011, 111, 3858-3887. 2. Tang, Y.; Ouyang, M. Tailoring Properties and Functionalities of Metal Nanoparticles through Crystallinity Engineering. Nat. Mater. 2007, 6, 754-759. 3. Pelton, M.; Chakraborty, D.; Malachosky, E.; Guyot-Sionnest, P.; Sader, J. E. Viscoelastic Flows in Simple Liquids Generated by Vibrating Nanostructures. Phys. Rev. Lett. 2013, 111, 244502. 4. Yu, K.; Major, T. A.; Chakraborty, D.; Devadas, M. S.; Sader, J. E.; Hartland, G. V. Compressible Viscoelastic Liquid Effects Generated by the Breathing Modes of Isolated Metal Nanowires. Nano Lett. 2015, 15, 3964-3970. 5. Dacosta Fernandes, B.; Spuch-Calvar, M.; Baida, H.; Tréguer-Delapierre, M.; Oberlé, J.; Langot, P.; Burgin, J. Acoustic Vibrations of Au Nano-Bipyramids and Their Modification under Ag Deposition: A Perspective for the Development of Nanobalances. ACS Nano 2013, 7, 76307639. 6. Arlett, J. L.; Myers, E. B.; Roukes, M. L. Comparative Advantages of Mechanical Biosensors. Nat. Nanotechnol. 2011, 6, 203-215. 7. Stiles, P. L.; Dieringer, J. A.; Shah, N. C.; Van Duyne, R. P. Surface-Enhanced Raman Spectroscopy. Annu. Rev. Anal. Chem. 2008, 1, 601-626. 8. Homola, J.; Yee, S. S.; Gauglitz, G. Surface Plasmon Resonance Sensors: Review. Sens. Actuators, B 1999, 54, 3-15. 9. Anker, J. N.; Hall, W. P.; Lyandres, O.; Shah, N. C.; Zhao, J.; Van Duyne, R. P. Biosensing with Plasmonic Nanosensors. Nat. Mater. 2008, 7, 442-453. 10. Willets, K. A.; Van Duyne, R. P. Localized Surface Plasmon Resonance Spectroscopy and Sensing. Annu. Rev. Phys. Chem. 2007, 58, 267-297. 11. Pelton, M.; Bryant, G. W. Introduction to Metal-Nanoparticle Plasmonics. Wiley: New Jersey, 2013. 12. Hartland, G. V. Measurements of the Material Properties of Metal Nanoparticles by Time-Resolved Spectroscopy. Phys. Chem. Chem. Phys. 2004, 6, 5263-5274. 13. Hartland, G. V. Coherent Excitation of Vibrational Modes in Metallic Nanoparticles. Annu. Rev. Phys. Chem. 2006, 57, 403-430. 14. Burgin, J.; Langot, P.; Arbouet, A.; Margueritat, J.; Gonzalo, J.; Afonso, C. N.; Vallée, F.; Mlayah, A.; Rossell, M. D.; Van Tendeloo, G. Acoustic Vibration Modes and Electron– Lattice Coupling in Self-Assembled Silver Nanocolumns. Nano Lett. 2008, 8, 1296-1302. 15. Pelton, M.; Sader, J. E.; Burgin, J.; Liu, M.; Guyot-Sionnest, P.; Gosztola, D. Damping of Acoustic Vibrations in Gold Nanoparticles. Nat. Nanotechnol. 2009, 4, 492-495. 16. Ruijgrok, P. V.; Zijlstra, P.; Tchebotareva, A. L.; Orrit, M. Damping of Acoustic Vibrations of Single Gold Nanoparticles Optically Trapped in Water. Nano Lett. 2012, 12, 10631069. 17. Yu, K.; Zijlstra, P.; Sader, J. E.; Xu, Q.-H.; Orrit, M. Damping of Acoustic Vibrations of Immobilized Single Gold Nanorods in Different Environments. Nano Lett. 2013, 13, 2710-2716. 18. Bachelier, G.; Mlayah, A. Surface Plasmon Mediated Raman Scattering in Metal Nanoparticles. Phys. Rev. B 2004, 69, 205408. 19. Lamb, H. On the Vibrations of an Elastic Sphere. Proc. London Math. Soc. 1881, s1-13, 189-212.

27

ACS Paragon Plus Environment

Page 28 of 31

Page 29 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

20. Hodak, J. H.; Henglein, A.; Hartland, G. V. Size Dependent Properties of Au Particles: Coherent Excitation and Dephasing of Acoustic Vibrational Modes. J. Chem. Phys. 1999, 111, 8613-8621. 21. Crut, A.; Juvé, V.; Mongin, D.; Maioli, P.; Del Fatti, N.; Vallée, F. Vibrations of Spherical Core-Shell Nanoparticles. Phys. Rev. B 2011, 83, 205430. 22. Girard, A.; Gehan, H.; Crut, A.; Mermet, A.; Saviot, L.; Margueritat, J. Mechanical Coupling in Gold Nanoparticles Supermolecules Revealed by Plasmon-Enhanced Ultralow Frequency Raman Spectroscopy. Nano Lett. 2016, 16, 3843-3849. 23. Crut, A.; Maioli, P.; Del Fatti, N.; Vallée, F. Acoustic Vibrations of Metal Nano-Objects: Time-Domain Investigations. Phys. Rep. 2015, 549, 1-43. 24. Large, N.; Saviot, L.; Margueritat, J.; Gonzalo, J.; Afonso, C. N.; Arbouet, A.; Langot, P.; Mlayah, A.; Aizpurua, J. Acousto-Plasmonic Hot Spots in Metallic Nano-Objects. Nano Lett. 2009, 9, 3732-3738. 25. Tripathy, S.; Marty, R.; Lin, V. K.; Teo, S. L.; Ye, E.; Arbouet, A.; Saviot, L.; Girard, C.; Han, M. Y.; Mlayah, A. Acousto-Plasmonic and Surface-Enhanced Raman Scattering Properties of Coupled Gold Nanospheres/Nanodisk Trimers. Nano Lett. 2011, 11, 431-437. 26. Nelet, A.; Crut, A.; Arbouet, A.; Del Fatti, N.; Vallée, F.; Portalès, H.; Saviot, L.; Duval, E. Acoustic Vibrations of Metal Nanoparticles: High Order Radial Mode Detection. Appl. Surf. Sci. 2004, 226, 209-215. 27. Del Fatti, N.; Voisin, C.; Chevy, F.; Vallée, F.; Flytzanis, C. Coherent Acoustic Mode Oscillation and Damping in Silver Nanoparticles. J. Chem. Phys. 1999, 110, 11484-11487. 28. van Dijk, M. A.; Lippitz, M.; Orrit, M. Detection of Acoustic Oscillations of Single Gold Nanospheres by Time-Resolved Interferometry. Phys. Rev. Lett. 2005, 95, 267406. 29. Hu, M.; Wang, X.; Hartland, G. V.; Mulvaney, P.; Juste, J. P.; Sader, J. E. Vibrational Response of Nanorods to Ultrafast Laser Induced Heating. Theoretical and Experimental Analysis. J. Am. Chem. Soc. 2003, 125, 14925-14933. 30. Park, S.; Pelton, M.; Liu, M.; Guyot-Sionnest, P.; Scherer, N. F. Ultrafast Resonant Dynamics of Surface Plasmons in Gold Nanorods. J. Phys. Chem. C 2007, 111, 116-123. 31. Zijlstra, P.; Tchebotareva, A. L.; Chon, J. W. M.; Gu, M.; Orrit, M. Acoustic Oscillations and Elastic Moduli of Single Gold Nanorods. Nano Lett. 2008, 8, 3493-3497. 32. Baida, H.; Mongin, D.; Christofilos, D.; Bachelier, G.; Crut, A.; Maioli, P.; Del Fatti, N.; Vallée, F. Ultrafast Nonlinear Optical Response of a Single Gold Nanorod near Its Surface Plasmon Resonance. Phys. Rev. Lett. 2011, 107, 057402. 33. Soavi, G.; Tempra, I.; Pantano, M. F.; Cattoni, A.; Collin, S.; Biagioni, P.; Pugno, N. M.; Cerullo, G. Ultrasensitive Characterization of Mechanical Oscillations and Plasmon Energy Shift in Gold Nanorods. ACS Nano 2016, 10, 2251–2258. 34. Guillon, C.; Langot, P.; Del Fatti, N.; Vallée, F.; Kirakosyan, A. S.; Shahbazyan, T. V.; Cardinal, T.; Treguer, M. Coherent Acoustic Vibration of Metal Nanoshells. Nano Lett. 2007, 7, 138-142. 35. Dowgiallo, A.-M.; Schwartzberg, A. M.; Knappenberger, K. L. Structure-Dependent Coherent Acoustic Vibrations of Hollow Gold Nanospheres. Nano Lett. 2011, 11, 3258-3262. 36. Staleva, H.; Hartland, G. V. Vibrational Dynamics of Silver Nanocubes and Nanowires Studied by Single-Particle Transient Absorption Spectroscopy. Adv. Funct. Mater. 2008, 18, 3809-3817. 37. Staleva, H.; Hartland, G. V. Transient Absorption Studies of Single Silver Nanocubes. J. Phys. Chem. C 2008, 112, 7535-7539.

28

ACS Paragon Plus Environment

ACS Nano

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

38. Petrova, H.; Lin, C.-H.; de Liejer, S.; Hu, M.; McLellan, J. M.; Siekkinen, A. R.; Wiley, B. J.; Marquez, M.; Xia, Y.; Sader, J. E.; Hartland, G. V. Time-Resolved Spectroscopy of Silver Nanocubes: Observation and Assignment of Coherently Excited Vibrational Modes. J. Chem. Phys. 2007, 126, 094709. 39. Galstyan, V.; Pak, O. S.; Stone, H. A. A Note on the Breathing Mode of an Elastic Sphere in Newtonian and Complex Fluids. Phys. Fluids 2015, 27, 032001. 40. Eringen, A. C.; Şuhubi, E. S. Linear Theory. Academic Press: New York, 1975; pp 717862. 41. Hartland, G. V. Coherent Vibrational Motion in Metal Particles: Determination of the Vibrational Amplitude and Excitation Mechanism. J. Chem. Phys. 2002, 116, 8048-8055. 42. Crut, A.; Maioli, P.; Fatti, N. D.; Vallee, F. Anisotropy Effects on the Time-Resolved Spectroscopy of the Acoustic Vibrations of Nanoobjects. Phys. Chem. Chem. Phys. 2009, 11, 5882-5888. 43. Prodan, E.; Radloff, C.; Halas, N. J.; Nordlander, P. A Hybridization Model for the Plasmon Response of Complex Nanostructures. Science 2003, 302, 419-422. 44. Fuchs, R. Theory of the Optical Properties of Ionic Crystal Cubes. Phys. Rev. B 1975, 11, 1732-1740. 45. Zhang, S.; Bao, K.; Halas, N. J.; Xu, H.; Nordlander, P. Substrate-Induced Fano Resonances of a Plasmonic Nanocube: A Route to Increased-Sensitivity Localized Surface Plasmon Resonance Sensors Revealed. Nano Lett. 2011, 11, 1657-1663. 46. Hu, M.; Hillyard, P.; Hartland, G. V.; Kosel, T.; Perez-Juste, J.; Mulvaney, P. Determination of the Elastic Constants of Gold Nanorods Produced by Seed Mediated Growth. Nano Lett. 2004, 4, 2493-2497. 47. Arbouet, A.; Del Fatti, N.; Vallee, F. Optical Control of the Coherent Acoustic Vibration of Metal Nanoparticles. J. Chem. Phys. 2006, 124, 144701. 48. Del Fatti, N.; Voisin, C.; Christofilos, D.; Vallée, F.; Flytzanis, C. Acoustic Vibration of Metal Films and Nanoparticles. J. Phys. Chem. A 2000, 104, 4321-4326. 49. Chang, W.-S.; Wen, F.; Chakraborty, D.; Su, M.-N.; Zhang, Y.; Shuang, B.; Nordlander, P.; Sader, J. E.; Halas, N. J.; Link, S. Tuning the Acoustic Frequency of a Gold Nanodisk through Its Adhesion Layer. Nat. Commun. 2015, 6, 7022. 50. O’Brien, K.; Lanzillotti-Kimura, N. D.; Rho, J.; Suchowski, H.; Yin, X.; Zhang, X. Ultrafast Acousto-Plasmonic Control and Sensing in Complex Nanostructures. Nat. Commun. 2014, 5, 4042. 51. Cardinal, M. F.; Mongin, D.; Crut, A.; Maioli, P.; Rodríguez-González, B.; Pérez-Juste, J.; Liz-Marzán, L. M.; Del Fatti, N.; Vallée, F. Acoustic Vibrations in Bimetallic Au@Pd Core– Shell Nanorods. J. Phys. Chem. Lett. 2012, 3, 613-619. 52. Masia, F.; Langbein, W.; Borri, P. Measurement of the Dynamics of Plasmons inside Individual Gold Nanoparticles Using a Femtosecond Phase-Resolved Microscope. Phys. Rev. B 2012, 85, 235403. 53. Ayedh, H. M.; Wacker, A. Acoustic Phonons in Nanowires with Embedded Heterostructures. J. Nanomater. 2011, 2011, 7. 54. Juvé, V.; Crut, A.; Maioli, P.; Pellarin, M.; Broyer, M.; Del Fatti, N.; Vallée, F. Probing Elasticity at the Nanoscale: Terahertz Acoustic Vibration of Small Metal Nanoparticles. Nano Lett. 2010, 10, 1853-1858. 55. Combe, N.; Saviot, L. Acoustic Modes in Metallic Nanoparticles: Atomistic Versus Elasticity Modeling. Phys. Rev. B 2009, 80, 035411.

29

ACS Paragon Plus Environment

Page 30 of 31

Page 31 of 31

1 2 3 4 5 6 7 8 9 10 11 12 13 14 15 16 17 18 19 20 21 22 23 24 25 26 27 28 29 30 31 32 33 34 35 36 37 38 39 40 41 42 43 44 45 46 47 48 49 50 51 52 53 54 55 56 57 58 59 60

ACS Nano

56. Saviot, L.; Netting, C. H.; Murray, D. B. Damping by Bulk and Shear Viscosity of Confined Acoustic Phonons for Nanostructures in Aqueous Solution. J. Phys. Chem. B 2007, 111, 7457-7461. 57. Yang, H. U.; D'Archangel, J.; Sundheimer, M. L.; Tucker, E.; Boreman, G. D.; Raschke, M. B. Optical Dielectric Function of Silver. Phys. Rev. B 2015, 91, 235137. 58. Olmon, R. L.; Slovick, B.; Johnson, T. W.; Shelton, D.; Oh, S.-H.; Boreman, G. D.; Raschke, M. B. Optical Dielectric Function of Gold. Phys. Rev. B 2012, 86, 235147. 59. Szczepanek, P.; Glosser, R. Piezo-Optical Constants of Gold. Solid State Commun. 1974, 15, 1425-1429. 60. Tups, H.; Otto, A.; Syassen, K. Pressure Dependence of the D-Band to Fermi-Level Excitation Threshold in Silver. Phys. Rev. B 1984, 29, 5458-5461. 61. Lee, A.; Andrade, G. F. S.; Ahmed, A.; Souza, M. L.; Coombs, N.; Tumarkin, E.; Liu, K.; Gordon, R.; Brolo, A. G.; Kumacheva, E. Probing Dynamic Generation of Hot-Spots in SelfAssembled Chains of Gold Nanorods by Surface-Enhanced Raman Scattering. J. Am. Chem. Soc. 2011, 133, 7563-7570.

30

ACS Paragon Plus Environment