Coherent vibrational dynamics of Au₁₄₄(SR)₆₀ nanoclusters

Abstract

The coherent vibrational dynamics of gold nanoclusters (NCs) provides important information on the coupling between vibrations and electrons as well as their mechanical properties, which is critical for understanding the evolution from a metallic state to a molecular state with diminishing size. Coherent vibrations have been widely explored in small-sized atomically precise gold NCs, while it remains a challenge to observe them in large-sized gold NCs. In this work, we report the coherent vibrational dynamics of atomically precise Au₁₄₄(SR)₆₀ NCs via temperature-dependent femtosecond transient absorption (TA) spectroscopy. The population dynamics of Au₁₄₄(SR)₆₀ consists of three relaxation processes: internal conversion, core–shell charge transfer and relaxation to the ground state. After removing the population dynamics from the TA kinetics, fast Fourier transform analysis on the residual oscillation reveals distinct vibrational modes at 1.5 THz (50 cm⁻¹) and 2 THz (67 cm⁻¹), which arise from the wavepacket motions along the ground-state and excited-state potential energy surfaces (PES), respectively. These results are helpful for understanding the physical properties of gold nanostructures with a threshold size that lies in between those of molecular-like NCs and metallic-state nanoparticles.

---

Introduction

Gold nanostructures exhibit unique optical and electronic properties and have received wide attention from researchers owing to promising applications in optical sensing, catalysis, therapeutics, etc.^1–6 These distinct properties are highly dependent on the structure, size, and metal composition of nanocrystals. Gold nanoparticles (NPs) with a diameter larger than ∼2 nm are well known to show surface plasmon resonance (SPR) because of the collective oscillations of conduction band electrons.^2,7 However, the quantum confinement effect starts to govern the optical properties of sub-2 nm nanoparticles (also known as nanoclusters, NCs), which possess discrete energy levels and molecular-like properties such as single-electron transitions or excitons.^4,8,9 The acoustic response of these nanostructures with different sizes has been widely reported,^8–18 which originate from different mechanisms. Normally, the photoexcited coherent breathing vibration mode for gold NPs is explained as the impulsive heating of the lattice upon short pulse laser excitation that expands the volume of NPs and results in decreased electronic density, which further gives rise to lattice vibration and periodic shift of the SPR band.^11,19 In molecular-like gold NCs, however, the coherent oscillations arise from the wavepacket motions along the ground- and excited-state potential energy surface (PES) generated by ultrafast laser pulses.^20–22 The coherent vibrations in small-sized atomically precise gold NCs were well explored experimentally or theoretically, including Au₁₀, Au₁₅, Au₂₅, Au₃₀, Au₃₇, Au₃₈, Au₁₀₂, etc.,^23–29 while such phenomena in large-sized gold NCs (n > 100) have rarely been observed.³⁰ The coherent vibrations have also received wide coverage in the literature previously, including organics, inorganics, biological systems, light harvesting systems, etc.^31–44

Among different sized gold NCs, Au₁₄₄L₆₀ (L = ligands) NCs play a vital role in understanding the size-dependent electronic and optical evolution. Au₁₄₄L₆₀ has a diameter of ∼2 nm that lies in between those of ultrasmall molecular NCs and larger-size plasmonic gold nanoparticles. The crystal structure of Au₁₄₄(SR)₆₀ was successfully solved by Yan et al. in 2018.⁴⁵ They found that the inner core is a hollow icosahedral 12-gold unit (the first layer, Au₁₂, magenta in Fig. 1a), which shows a similar structure to the center-filled icosahedral 13-gold structure (Au₁₃, Fig. 1b and c). The second layer is composed of 42 gold atoms (blue in Fig. 1a), and then 60 gold atoms (third layer, light green in Fig. 1a) enclose the 54 (12 + 42) gold atoms and assemble into a 114-gold core (Au₁₁₄), on which 30 Au(SR)₂ staples (orange in Fig. 1a, non-gold atoms not shown) ride astride, giving rise to a complete Au₁₄₄(SR)₆₀ structure. In the same year, Lei et al.⁴⁶ reported the crystal structure of Au₁₄₄(C≡CAr)₆₀ (Ar = 2-FC₆H₄⁻), which showed a similar kernel structure to that of Au₁₄₄(SR)₆₀.

Fig. 1 (a) The crystal structure of an Au₁₄₄(SR)₆₀ NC:⁴⁵ the inner core Au₁₂ (first layer) is shown in magenta, the second layer of 42 gold atoms is shown in blue, the third layer of 60 gold atoms is shown in light green, and the 30 Au atoms in staples are shown in orange. The non-gold atoms are not shown for clarity. (b) The structure of the hollow Au₁₂ inner core of Au₁₄₄(SR)₆₀. (c) The structure of icosahedral Au₁₃; the center gold atom is shown in yellow. (d) The matrix-assisted laser desorption ionization (MALDI) mass spectrum of Au₁₄₄PET₆₀; note that the peak at 17.2 k is doubly charged and the 34 k peak is 1+ charged (in-source ionized). (e) The steady-state absorption spectrum of Au₁₄₄PET₆₀ dissolved in toluene.

Besides the crystal structure determination of Au₁₄₄(SR)₆₀ NCs, the physical properties have also attracted considerable research interest. Malola et al.⁴⁷ reported a computational study of several monolayer-protected Au_n NCs by applying the ab initio method together with the atomistic method, and they found that the threshold gold-core diameter for the emergency of the localized surface plasmon effect is between 1.5 nm (Au₁₄₄(SH)₆₀) and 2 nm (Au₃₁₄(SH)₉₆). Weissker et al.⁴⁸ presented a low-temperature absorption spectrum of Au₁₄₄(SR)₆₀ with clearly visible individual peaks from the UV to the NIR region, conveying quantum confinement characters in the clusters. Yi et al.⁴⁹ provided an electronic relaxation study of [Au₁₄₄(SC₆H₁₃)₆₀]^q (q = −1, 0, +1, +2) by employing femtosecond transient absorption spectroscopy, and they thought that the pump power-dependent bleach recovery kinetics and oxidation state-dependent excited-state absorption kinetics should be assigned to metallic behavior while stimulated emission showed signatures of “superatom” character. Mustalahti et al.⁵⁰ reported the ultrafast electronic relaxation and vibrational cooling dynamics of Au₁₄₄(SC₂H₄Ph)₆₀ using transient mid-IR spectroscopy, in which two relaxation processes corresponding to fast relaxation via electron–phonon coupling and heat dissipation were found. Low-temperature absorption spectroscopy and density functional theory (DFT) calculations conducted by Negishi et al.⁵¹ revealed that the metal core of Au₁₄₄(SC₁₂)₆₀ and smaller NCs have molecular-like electronic structures, whereas the core structures of larger NCs resemble that of bulk gold. Zhou et al.⁵² performed systematic ultrafast spectroscopy studies on atomically precise gold NCs with different sizes, and they found that Au₁₄₄(SR)₆₀ and smaller NCs exhibit molecular-like behavior. Lei et al.⁴⁶ reported that the steady-state absorption spectrum of Au₁₄₄(C≡CAr)₆₀ NCs shows two bands at around 560 and 620 nm, indicating molecular-like behavior. Besides, they found that the Au–Au distances in the Au₁₂ core (2.7470–2.7587 Å) are shorter than the average Au–Au distance of bulk face-centered cubic gold (2.884 Å), which also suggests that the bonding is more molecular than metallic. Deng et al.⁵³ found that the origin of the structural stability of Au₁₄₄(SR)₆₀ and Au₁₄₄(C≡CR)₆₀ is due to the thermodynamic stability by applying DFT energy calculation and energy decomposition analysis (EDA); in addition, they also suggested that the Au₁₄₄ NCs are unique in size because the larger (Au₂₇₉) and smaller (Au₁₃₀, Au₁₃₃) gold NCs both prefer the center-filled gold polyhedrons. Recent work by Du et al.⁵⁴ reported temperature-dependent (room temperature to 2.5 K) absorption spectra, in which they observed several new features at cryogenic temperatures, which provides strong evidence of discrete energy levels within Au₁₄₄PET₆₀ (where, PET = SCH₂CH₂Ph).

Despite the great progress made in the study of the photophysical properties of Au₁₄₄L₆₀ NCs, the coherent vibrational dynamics still remains unknown, which is thus critical to pursue, and abundant information about interactions between the ambient environment and electronic relaxation may be revealed. In this work, using temperature-dependent femtosecond broadband transient absorption (TA) spectroscopy, we report the photoinduced excited-state population and vibration dynamics of Au₁₄₄PET₆₀ NCs. By tracking the kinetic traces at different probe wavelengths, the coherent oscillations with frequencies at around 1.5 THz (50 cm⁻¹) and 2 THz (67 cm⁻¹) are obtained. Furthermore, detailed analysis of coherent oscillations has identified the possible contributions of the two vibrational modes in Au₁₄₄PET₆₀.

Results and discussion

The synthesis of Au₁₄₄PET₆₀ NCs followed a previous method,⁵⁵ and the MALDI mass spectrum of Au₁₄₄PET₆₀ is displayed in Fig. 1d, suggesting the high purity of the sample. Shown in Fig. 1e is the steady-state absorption spectrum of Au₁₄₄PET₆₀ in toluene, spanning from the UV to the NIR region. There are several broad absorption peaks at around 337, 463, 517 and 700 nm, which agrees with previous reports.^50,54,55 From the continuous absorption bands in Au₁₄₄PET₆₀, it is difficult to extract the energy gap. Koivisto et al.⁵⁶ collected the optical spectra of an Au₁₄₄PET₆₀ thin film in IR and NIR regimes, and combined with linear response time-dependent DFT (TD-DFT) calculations, the optical gap of Au₁₄₄PET₆₀ was determined to be 0.19 eV. Chen et al.⁵⁷ revealed a 0.17 eV gap for Au₁₄₄(BM)₆₀ (BM is the abbreviation of benzyl mercaptan) at 298 K by means of cyclic and differential pulse voltammetry experiments. The conclusions of these two reports coincide with each other, suggesting a very small energy gap of Au₁₄₄ NCs (around 0.2 eV).

To probe the excited-state population and vibration dynamics, we first collected the broadband femtosecond transient absorption of Au₁₄₄PET₆₀ dissolved in toluene in the visible spectral region with a pump laser of 400 nm (3.1 eV). As shown in Fig. S1,† there is no sign of obvious coherent oscillations in the pseudo-color map, which may be due to the weak coherent vibrations of Au₁₄₄PET₆₀ in solution that can be easily perturbed by the solvent fluctuation and environmental heat. To obtain coherent oscillation dynamics, the temperature-dependent TA spectra of Au₁₄₄PET₆₀ dispersed in a PMMA thin film were recorded. As there is no absorption in the visible region for PMMA, the TA experiments of Au₁₄₄PET₆₀ cannot be affected by the polymer film.

Fig. 2a shows the TA data map of the Au₁₄₄PET₆₀ film at 77 K; the energy of the excitation pulse (3.1 eV) is significantly higher than the band gap of ∼0.2 eV.^56,57 Significant oscillations superimposed on the TA dynamics can be observed. 400 nm excitation forms an excited-state that is localized in the Au₁₁₄ core (12 + 42 + 60 = 114, as shown in Fig. 1a and b) based on a previous report,⁴⁸ which is named the core state. Immediately after photon excitation, broad excited-state absorption (ESA) that spanned over most of the probed region was observed instantaneously; and most of the excited-state population decays to the ground state within the first 10 ps (Fig. 2a and b). After global analysis on the TA data,⁵⁸ the wavelength-dependent amplitudes of characteristic transient components were extracted. As shown in Fig. 2c, the best fitting results with three decay-associated spectra (DAS) are portrayed, and temporal analysis of the three distinct components gives time constants of 1.3 ps, 2.0 ps, and 24.6 ps, respectively. One could observe that DAS1 and DAS2 are almost negatively correlated,²⁴ implying that the amplitude of DAS2 increases on a time scale of 1.3 ps. We assign the first evolution process from DAS1 to DAS2 to internal conversion (from the core state to a relaxed core state) because the 400 nm excitation gives rise to excess excited-state energy, which can also be further evidenced by the kinetic traces shown in Fig. 2d. The ESA at 740 nm was obtained immediately after photoexcitation, whereas 600 nm ESA exhibits a slow growth process in the first ∼1.7 ps, which points to an internal conversion process from the core state to a relaxed core state, i.e., the ESA at 600 nm originates from the relaxed core state to higher excited states. Subsequently, the initial excitation localized in the Au₁₁₄ core relaxes to a core–shell state because of energy redistribution. The last time-constant (24.6 ps) should be assigned to the relaxation to the ground state. Compared to those smaller gold NCs, the relatively fast decay herein results from the ultrasmall energy gap (∼0.2 eV) of Au₁₄₄PET₆₀.³⁰

Fig. 2 (a) The TA data map of Au₁₄₄PET₆₀ dispersed in a polymethyl methacrylate (PMMA) thin film under the excitation of 400 nm at 77 K. (b) The TA spectra as a function of time delay within 7 ps. (c) Decay-associated spectra (DAS) obtained from global analysis on the TA data. (d) TA kinetic traces probed at different wavelengths and the corresponding fitting results of Au₁₄₄PET₆₀.

The TA spectra at 295 K in the film phase are shown in the S2 part of the ESI.† As shown in Fig. S2a and b,† the TA spectra show some clear trends as the temperature decreases: (1) both ESA and ground-state bleach (GSB) became shaper; (2) both GSB and ESA evolved to higher energies; (3) the intensity of the entire spectra increased. A similar phenomenon was observed in rod-shaped [Au₃₇(PPh₃)₁₀(SC₂H₄Ph)₁₀X₂]⁺ (X = Cl/Br) previously,²⁹ in which the blue-shift and enhancement of ESA at low temperatures implied that the electron–phonon coupling played a vital role in the relaxation process. To understand the thermal effect on the electronic dynamics, the normalized TA kinetic traces measured at different temperatures and corresponding fitting are shown in Fig. S2c–f.† One could find that the electronic dynamics measured at 295 K is slightly faster than that measured at 77 K. This suggests that the heat transfer/dissipation are not the dominant factors during the excited-state relaxation. The blue-shift of the GSB signal is attributed to the increase of the HOMO–LUMO gap as the temperature decreases,⁵⁷ which was also observed in Au₃₅, Au₃₈, etc.,^51,54,59,60 with these phenomena often explained as the interaction between the Au core and phonons at the staple.^51,59 The detailed TA spectra and fitting results are shown in Fig. S3 and Table S1.†

Besides the population dynamics, the coherent vibrational dynamics also provides important information of the excited-state behavior. Upon resonant pump excitation (400 nm), the electronic transition can project a wavepacket onto the excited-state potential, and a wavepacket was formed on the ground-state potential through resonant Raman transition.⁶¹ In the initial 3 picoseconds, one can observe significant coherent vibrations at almost all probed wavelengths (Fig. 2a), and the extracted kinetics at selected wavelengths are shown in Fig. S4.† By subtracting the population dynamics of the kinetic traces, the fast Fourier transform (FFT) analyses of the oscillatory components (Fig. 3a–c) give two frequencies of ∼1.5 THz (50 cm⁻¹) and ∼2 THz (67 cm⁻¹) that predominantly contribute to the oscillations of the TA spectra. After careful analysis on the coherent oscillation in Au₁₄₄PET₆₀ NCs, we found that the ∼1.5 THz oscillation dominates the GSB signal between 500 nm and 550 nm. ∼2 THz oscillation is more significant from 550 nm to 700 nm, in which the excited-state components play a vital role in this region. Of note, there is a negative peak at 660 to 680 nm that was assigned to the stimulated emission (SE) signal previously,⁴⁹ which also reflects the excited-state dynamics. Overall, we can confirm that the ∼1.5 THz mode originates from the wavepacket motions in the ground-state PES, whereas the ∼2 THz mode arises from excited-state PES.

Fig. 3 The FFT results (a–c) and analysis of oscillatory components (d–f) on the coherent vibrational dynamics of Au₁₄₄PET₆₀ (77 K) probed at 505 nm, 588 nm, and 675 nm.

The oscillatory components were fitted with a curve that resembles a cosine curve oscillating in the envelope of an exponential damping function and considering the instrumental response function. The convolution form of the damped oscillation is expressed as follows:^24,62,63

where

f _osc(t) represents the cosine oscillating component, f_damp(t) represents the damping component, and f_IRF(t) represents the instrumental response function (IRF). The form of convolution can be rewritten as

in which A denotes the amplitude of the oscillation, t_osc denotes the time constant of the oscillation, t_d denotes the time constant of the damping, t_c denotes the time zero of the IRF, and FWHM denotes the full-width-half-maximum of the IRF, φ is the relative phase of the oscillation, and .

The oscillation period (t_osc) monitored at 505 nm was fitted to be 0.64 ps, while t_osc probed at 588 nm and 675 nm was ∼0.51 ps. Analysis on the transient absorption at 295 K was also performed, but the oscillatory signals are weaker than those at 77 K. The detailed coherent-vibration dynamics of Au₁₄₄PET₆₀ obtained at 295 K are shown in Fig. S5 and S6.† Comparing the vibrational data in solution (295 K) and film (77 K and 295 K) states, we found that (1) there are no vibrational signals in the solution state (Fig. S1†), while oscillatory features can be clearly observed in the film state (Fig. 2a and S5a†); (2) at cryogenic temperature, the coherent vibrations become more significant. These can be explained by: (1) the conformational relaxation and vibration of ligands and the Au core will be suppressed in the film state; (2) the reduced environmental noise at low temperatures will enhance the coherent oscillations in the TA kinetics of Au₁₄₄PET₆₀. The frequency comparison of Au₁₄₄PET₆₀ collected at 77 K and 295 K is shown in Fig. 4 and Table 1 summarises the parameters from the fitting of oscillatory data.

Fig. 4 The oscillation frequencies as a function of probe wavelength at 77 K and 295 K. Note that the error bars were calculated based on the frequencies probed at around 505 nm, 588 nm, and 675 nm, respectively.

Table 1 The parameters of coherent oscillations obtained from the TA data

Wavelength/nm	Frequency/THz		t osc/ps		Assignment
	77 K	295 K	77 K	295 K
505	∼1.5		0.63	0.66	Ground-state PES
588	∼2		0.52	0.55	Excited-state PES
675	∼2		0.50	0.51	Excited-state PES

---

The assignment of the oscillatory components is crucial for understanding the excited-state deactivation dynamics. As we discussed above, the inner core of Au₁₄₄ NCs is hollow icosahedral Au₁₂, which is structurally similar to center-filled icosahedral Au₁₃ (see Fig. 1). Previous coherent oscillatory analysis of rod-shaped [Au₃₇(PPh₃)₁₀(SC₂H₄Ph)₁₀X₂]⁺ (X = Cl/Br) indicated that a 2.1 THz (70 cm⁻¹) vibration arises from the radial breathing mode of the Au₁₃ unit.²⁹ We also note that recent work by Martinet et al.⁶⁴ reported the low-frequency Raman spectra of Au₂₅(SG)₁₈ (SG = glutathione) and Au₂₅(p-MBA)₁₈ (p-MBA = 4-mercaptobenzoic acid) NCs, in which the 1.9 THz vibration was assigned to the breathing mode of the icosahedral Au₁₃ core. Considering that there are several layers of gold atoms that enclose the Au₁₂ inner core within the Au₁₄₄PET₆₀ structure, the 2 THz oscillations of Au₁₄₄PET₆₀ should be ascribed to the breathing mode of the icosahedral Au₁₂ unit. We also observed that the TA profiles probed at 6 ps of Au₁₄₄ and Au₂₅ are similar to each other (Fig. S7†). Such an observation may suggest some connections of the icosahedral cores (Au₁₂vs. Au₁₃), while further theoretical simulations are required to understand the similar excited-state spectral profiles of Au₂₅ and Au₁₄₄.

For large counterparts of gold NCs, the continuum mechanical model can be used to predict the coherent oscillation, which gives rise to the empirical relationship ν = 50/R (R is the particle radius).^24,65 By substituting R with 0.95 nm (note that the diameter of the entire Au₁₄₄PET₆₀ is 1.9 nm),⁶⁶ we obtained 52 cm⁻¹ (1.6 THz). Since Au₁₄₄PET₆₀ is a spherical structure, one could conclude that the 1.5 THz originates from the entire Au₁₄₄PET₆₀ NC, which is symmetrically distributed on the entire Au₁₄₄PET₆₀.

In addition, we note that Au₁₄₄PET₆₀ has a core–shell spherical structure, which is similar to [Au₁₂(SCH₃)₉]⁺, [Au₂₅(SCH₃)₁₈]⁻, and Au₁₀₂(SCH₃)₄₄.²³ Maioli et al.²³ found that all three gold NCs exhibit two vibrational modes: a low-frequency quadrupolar-like mode (QM for short) and a high-frequency breathing-acoustic mode (BM for short). The BM is a highly symmetric fundamental mode that is formed from in-phase radially symmetric vibration of each atom around the equilibrium position, while the QM is composed of periodical expansion and contraction along one direction accompanied by out-of-phase vibration in the perpendicular direction. It is deduced that the 2 THz and 1.5 THz vibrations of Au₁₄₄PET₆₀ may arise from the BM and QM, respectively. The accurate assignments of coherent vibration modes of Au₁₄₄ NCs still require further experimental and theoretical evidence.

Based on the TA data analysis, the excited-state deactivation diagram and the mechanism of coherent vibrations of Au₁₄₄ are shown in Fig. 5. As discussed above, the 600 nm ESA originates from the relaxed core state, whereas the 740 nm ESA arises from the core state. With 400 nm (3.1 eV) excitation, the Au₁₁₄ core is directly populated (core state), and then the energy is redistributed over the entire Au₁₄₄PET₆₀ NC after the 1.3 ps internal conversion is followed by a 2.0 ps charge transfer (Fig. 5a). The last process is the relaxation to the ground state. Such a fast excited-state relaxation processes agrees with the very small energy gap of Au₁₄₄PET₆₀ (∼0.2 eV). The energy gap decreases as temperature increases,⁵⁷ which accelerates the excited-state relaxation at higher temperatures. The 400 nm photoexcitation leads to simultaneous generation of wavepacket motions on the excited- and ground-states⁶¹ with frequencies of 2 THz and 1.5 THz, respectively (Fig. 5b). These wavepacket motions result in coherent oscillations in the TA kinetics probed at different wavelengths.

Fig. 5 (a) Proposed excited-state deactivation mechanism of Au₁₄₄PET₆₀. (b) Schematic diagram of the coherent vibrations of Au₁₄₄PET₆₀.

Understanding the excited-state population- and vibration-dynamics of gold NCs is crucial for their applications in solar cells,⁶⁷ photoluminescence⁶⁸ and nonlinear optical properties.⁶⁹ The population dynamics pictures the relaxation of excited states, while the coherent vibrational dynamics reflects the detailed interactions between electrons and structural vibrations. The excited-state processes of gold NCs, including charge transfer, energy transfer and excited-state coherence, are detrimental to boost the applications of NCs in photo-energy storage and conversion. Therefore, unravelling the excited-state relaxation mechanisms and measuring the accurate time-constants of the excited-state processes are of great importance.

Large-sized gold nanoclusters (more than 100 gold atoms) are critical for understanding the transition from discrete energy levels to continuous bands in metal nanoparticles. They exhibit broad absorption bands and relatively weak excited state signals. Therefore, it remains a challenge to accurately probe their population- and vibration-dynamics. Among those large sized gold NCs, Au₁₄₄(SR)₆₀ is a typical example and it has a very small E_g of 0.17 eV.^56,57 Therefore, it has no gap emission so that we chose ultrafast TA to probe the population and vibration dynamics of Au₁₄₄. By applying low-temperature TA, we successfully observed two coherent vibration modes at 77 K, which were hidden at room temperature. Detailed analysis on these coherent vibrations can reveal the origin of electron–phonon coupling in different sized nanoparticles. Based on the current observations and analysis, we found that temperature dependent TA can be a new strategy to investigate the coherent vibrational dynamics of other large-sized metal NCs, which should provide new insight into research in this field.

Conclusions

In summary, we report the coherent vibrational dynamics of Au₁₄₄PET₆₀ NCs by temperature-dependent transient absorption analysis. Upon photoexcitation, the excited-state relaxation pathways of Au₁₄₄PET₆₀ can be extracted with three processes, corresponding to the internal conversion, core–shell state and decay to the ground state, respectively. Significant coherent oscillations superimposed on the population dynamics can be observed in the first 3 ps at 77 K and 295 K. We demonstrate that the high-frequency 2 THz oscillation originates from wavepacket motions along the excited-state PES, whereas the low-frequency 1.5 THz oscillation arises from wavepacket motions along the ground-state PES. Based on previous results and the continuum mechanical model, we have proposed two possible assignments of the vibrational modes: (1) the 2 THz mode is localized on the 12-gold-atom inner core while the 1.5 THz mode is symmetrically distributed on the entire Au₁₄₄PET₆₀; (2) the 2 THz and 1.5 THz modes may also arise from the breathing mode and quadrupolar-like mode of Au₁₄₄PET₆₀, respectively. We also proposed that low temperature transient absorption should be a new strategy to probe the vibration dynamics of large-sized metal NCs. These results are very helpful for elucidating the electron-vibration interactions of gold NCs with size around 2 nm, which should stimulate broad research interest in this topic.

Data availability

Materials and methods, additional TA data and coherent vibrations have been uploaded as the ESI.† The data that support the findings of this study are available from the corresponding author(s) upon reasonable request.

Author contributions

The manuscript was written through contributions of all authors. All authors have given approval to the final version of the manuscript.

Conflicts of interest

There are no conflicts to declare.

Acknowledgements

M. Z. acknowledges the support from the Chinese Academy of Sciences (YSBR-007).

References

1. S. Yamazoe, K. Koyasu and T. Tsukuda, Acc. Chem. Res., 2014, 47, 816–824.
2. K. M. Mayer and J. H. Hafner, Chem. Rev., 2011, 111, 3828–3857.
3. G. Li and R. Jin, Acc. Chem. Res., 2013, 46, 1749–1758.
4. R. Jin, C. Zeng, M. Zhou and Y. Chen, Chem. Rev., 2016, 116, 10346–10413.
5. S. S. Kumar, K. Kwak and D. Lee, Electroanalysis, 2011, 23, 2116–2124.
6. R. García-Álvarez, L. Chen, A. Nedilko, A. Sánchez-Iglesias, A. Rix, W. Lederle, V. Pathak, T. Lammers, G. von Plessen, K. Kostarelos, L. M. Liz-Marzán, A. J. C. Kuehne and D. N. Chigrin, ACS Photonics, 2020, 7, 646–652.
7. G. V. Hartland, Chem. Rev., 2011, 111, 3858–3887.
8. M. Zhou, X. Du, H. Wang and R. Jin, ACS Nano, 2021, 15, 13980–13992.
9. M. Zhou and R. Jin, Annu. Rev. Phys. Chem., 2021, 72, 121–142.
10. B. T. Diroll, M. S. Kirschner, P. Guo and R. D. Schaller, Annu. Rev. Phys. Chem., 2019, 70, 353–377.
11. G. V. Hartland, Annu. Rev. Phys. Chem., 2006, 57, 403–430.
12. M. Iwamura, K. Kimoto, K. Nozaki, H. Kuramochi, S. Takeuchi and T. Tahara, J. Phys. Chem. Lett., 2018, 9, 7085–7089.
13. S. H. Yau, O. Varnavski and T. Goodson, Acc. Chem. Res., 2013, 46, 1506–1516.
14. M. Zhou, S. Vdović, S. Long, M. Zhu, L. Yan, Y. Wang, Y. Niu, X. Wang, Q. Guo, R. Jin and A. Xia, J. Phys. Chem. A, 2013, 117, 10294–10303.
15. R. J. Newhouse, H. Wang, J. K. Hensel, D. A. Wheeler, S. Zou and J. Z. Zhang, J. Phys. Chem. Lett., 2011, 2, 228–235.
16. R. Marty, A. Arbouet, C. Girard, A. Mlayah, V. Paillard, V. K. Lin, S. L. Teo and S. Tripathy, Nano Lett., 2011, 11, 3301–3306.
17. Y. U. Staechelin, D. Hoeing, F. Schulz and H. Lange, ACS Photonics, 2021, 8, 752–757.
18. M. S. Kirschner, C. M. Lethiec, X.-M. Lin, G. C. Schatz, L. X. Chen and R. D. Schaller, ACS Photonics, 2016, 3, 758–763.
19. M. Pelton, J. E. Sader, J. Burgin, M. Liu, P. Guyot-Sionnest and D. Gosztola, Nat. Nanotechnol., 2009, 4, 492–495.
20. C. J. Bardeen, Q. Wang and C. V. Shank, Phys. Rev. Lett., 1995, 75, 3410–3413.
21. G. D. Scholes, G. R. Fleming, L. X. Chen, A. Aspuru-Guzik, A. Buchleitner, D. F. Coker, G. S. Engel, R. van Grondelle, A. Ishizaki, D. M. Jonas, J. S. Lundeen, J. K. McCusker, S. Mukamel, J. P. Ogilvie, A. Olaya-Castro, M. A. Ratner, F. C. Spano, K. B. Whaley and X. Zhu, Nature, 2017, 543, 647–656.
22. S. G. Rafiq S, J. Phys. Chem. A, 2016, 120, 6297–6299.
23. P. Maioli, T. Stoll, H. E. Sauceda, I. Valencia, A. Demessence, F. Bertorelle, A. Crut, F. Vallée, I. L. Garzón, G. Cerullo and N. Del Fatti, Nano Lett., 2018, 18, 6842–6849.
24. M. Y. Sfeir, H. Qian, K. Nobusada and R. Jin, J. Phys. Chem. C, 2011, 115, 6200–6207.
25. M. Zhou, T. Higaki, G. Hu, M. Y. Sfeir, Y. Chen, D.-e. Jiang and R. Jin, Science, 2019, 364, 279–282.
26. R. D. Senanayake, E. B. Guidez, A. J. Neukirch, O. V. Prezhdo and C. M. Aikens, J. Phys. Chem. C, 2018, 122, 16380–16388.
27. H. Qian, M. Y. Sfeir and R. Jin, J. Phys. Chem. C, 2010, 114, 19935–19940.
28. M. Zhou, S. Tian, C. Zeng, M. Y. Sfeir, Z. Wu and R. Jin, J. Phys. Chem. C, 2017, 121, 10686–10693.
29. M. Zhou, R. Jin, M. Y. Sfeir, Y. Chen, Y. Song and R. Jin, Proc. Natl. Acad. Sci., 2017, 114, E4697–E4705.
30. M. Zhou, C. Zeng, Y. Song, J. W. Padelford, G. Wang, M. Y. Sfeir, T. Higaki and R. Jin, Angew. Chem., Int. Ed., 2017, 56, 16257–16261.
31. L. Bolzonello, F. Fassioli and E. Collini, J. Phys. Chem. Lett., 2016, 7, 4996–5001.
32. S. Rafiq, B. Fu, B. Kudisch and G. D. Scholes, Nat. Chem., 2021, 13, 70–76.
33. M. Park, A. J. Neukirch, S. E. Reyes-Lillo, M. Lai, S. R. Ellis, D. Dietze, J. B. Neaton, P. Yang, S. Tretiak and R. A. Mathies, Nat. Commun., 2018, 9, 2525.
34. A. Chenu and G. D. Scholes, Annu. Rev. Phys. Chem., 2015, 66, 69–96.
35. N. Christensson, H. F. Kauffmann, T. Pullerits and T. Mančal, J. Phys. Chem. B, 2012, 116, 7449–7454.
36. J. R. Caram, H. Zheng, P. D. Dahlberg, B. S. Rolczynski, G. B. Griffin, A. F. Fidler, D. S. Dolzhnikov, D. V. Talapin and G. S. Engel, J. Phys. Chem. Lett., 2014, 5, 196–204.
37. T. Pullerits, D. Zigmantas and V. Sundström, Proc. Natl. Acad. Sci., 2013, 110, 1148–1149.
38. A. S. Rury, S. A. Sorenson and J. M. Dawlaty, J. Phys. Chem. Lett., 2017, 8, 181–187.
39. T. Debnath, D. Sarker, H. Huang, Z.-K. Han, A. Dey, L. Polavarapu, S. V. Levchenko and J. Feldmann, Nat. Commun., 2021, 12, 2629.
40. G. D. Scholes, J. Phys. Chem. Lett., 2018, 9, 1568–1572.
41. A. W. Chin, J. Prior, R. Rosenbach, F. Caycedo-Soler, S. F. Huelga and M. B. Plenio, Nat. Phys., 2013, 9, 113–118.
42. P. J. M. Johnson, A. Halpin, T. Morizumi, V. I. Prokhorenko, O. P. Ernst and R. J. D. Miller, Nat. Chem., 2015, 7, 980–986.
43. S. M. Falke, C. A. Rozzi, D. Brida, M. Maiuri, M. Amato, E. Sommer, A. D. Sio, A. Rubio, G. Cerullo, E. Molinari and C. Lienau, Science, 2014, 344, 1001–1005.
44. E. Thyrhaug, S. A. Bogh, M. R. Carro-Temboury, C. S. Madsen, T. Vosch and D. Zigmantas, Nat. Commun., 2017, 8, 15577.
45. N. Yan, N. Xia, L. Liao, M. Zhu, F. Jin, R. Jin and Z. Wu, Sci. Adv., 2018, 4, eaat7259.
46. Z. Lei, J.-J. Li, X.-K. Wan, W.-H. Zhang and Q.-M. Wang, Angew. Chem., Int. Ed., 2018, 57, 8639–8643.
47. S. Malola, L. Lehtovaara, J. Enkovaara and H. Häkkinen, ACS Nano, 2013, 7, 10263–10270.
48. H. C. Weissker, H. B. Escobar, V. D. Thanthirige, K. Kwak, D. Lee, G. Ramakrishna, R. L. Whetten and X. López-Lozano, Nat. Commun., 2014, 5, 3785.
49. C. Yi, M. A. Tofanelli, C. J. Ackerson and K. L. Knappenberger, J. Am. Chem. Soc., 2013, 135, 18222–18228.
50. S. Mustalahti, P. Myllyperkiö, T. Lahtinen, K. Salorinne, S. Malola, J. Koivisto, H. Häkkinen and M. Pettersson, J. Phys. Chem. C, 2014, 118, 18233–18239.
51. Y. Negishi, T. Nakazaki, S. Malola, S. Takano, Y. Niihori, W. Kurashige, S. Yamazoe, T. Tsukuda and H. Häkkinen, J. Am. Chem. Soc., 2015, 137, 1206–1212.
52. M. Zhou, C. Zeng, Y. Chen, S. Zhao, M. Y. Sfeir, M. Zhu and R. Jin, Nat. Commun., 2016, 7, 13240.
53. S. Deng, J. Li, P. Wang and Y. Pei, Nanoscale, 2021, 13, 18134–18139.
54. X. Du, Z. Liu, T. Higaki, M. Zhou and R. Jin, Chem. Sci., 2022, 13, 1925–1932.
55. H. Qian and R. Jin, Chem. Mater., 2011, 23, 2209–2217.
56. J. Koivisto, S. Malola, C. Kumara, A. Dass, H. Häkkinen and M. Pettersson, J. Phys. Chem. Lett., 2012, 3, 3076–3080.
57. S. Chen, T. Higaki, H. Ma, M. Zhu, R. Jin and G. Wang, ACS Nano, 2020, 14, 16781–16790.
58. J. Snellenburg, S. Laptenok, R. Seger, K. Mullen and I. Van Stokkum, J. Stat. Software, 2012, 49, 1–22.
59. M. S. Devadas, S. Bairu, H. Qian, E. Sinn, R. Jin and G. Ramakrishna, J. Phys. Chem. Lett., 2011, 2, 2752–2758.
60. V. D. Thanthirige, M. Kim, W. Choi, K. Kwak, D. Lee and G. Ramakrishna, J. Phys. Chem. C, 2016, 120, 23180–23188.
61. C. Fitzpatrick, J. H. Odhner and R. J. Levis, J. Phys. Chem. A, 2020, 124, 6856–6866.
62. O. Varnavski, G. Ramakrishna, J. Kim, D. Lee and T. Goodson, ACS Nano, 2010, 4, 3406–3412.
63. https://webhome.phy.duke.edu/∼rgb/Class/phy51/phy51/node24https://webhome.phy.duke.edu/%7Ergb/Class/phy51/phy51/node24.html.
64. Q. Martinet, A. Berthelot, A. Girard, B. Donoeva, C. Comby-Zerbino, E. Romeo, F. Bertorelle, M. van der Linden, N. Tarrat, N. Combe and J. Margueritat, J. Phys. Chem. C, 2020, 124, 19324–19332.
65. J. H. Hodak, A. Henglein and G. V. Hartland, J. Phys. Chem. B, 2000, 104, 9954–9965.
66. J. Zhang, Z. Li, J. Huang, C. Liu, F. Hong, K. Zheng and G. Li, Nanoscale, 2017, 9, 16879–16886.
67. K. G. Stamplecoskie and A. Swint, J. Mater. Chem. A, 2016, 4, 2075–2081.
68. M. R. Narouz, S. Takano, P. A. Lummis, T. I. Levchenko, A. Nazemi, S. Kaappa, S. Malola, G. Yousefalizadeh, L. A. Calhoun, K. G. Stamplecoskie, H. Häkkinen, T. Tsukuda and C. M. Crudden, J. Am. Chem. Soc., 2019, 141, 14997–15002.
69. J. Olesiak-Banska, M. Waszkielewicz, P. Obstarczyk and M. Samoc, Chem. Soc. Rev., 2019, 48, 4087–4117.

---

Footnote

† Electronic supplementary information (ESI) available. See https://doi.org/10.1039/d2sc02246j

---