Open Access Article
This Open Access Article is licensed under a
Creative Commons Attribution 3.0 Unported Licence

Importance of double-resonance effects in two-photon absorption properties of Au25(SR)18

Zhongwei Hu and Lasse Jensen *
Department of Chemistry, The Pennsylvania State University, 104 Chemistry Building, University Park, Pennsylvania 16802, USA. E-mail:

Received 2nd March 2017 , Accepted 18th April 2017

First published on 19th April 2017

The two-photon absorption (TPA) cross-sections of small thiolate-protected gold clusters have been shown to be much larger than typical small organic molecules. In comparison with larger nanoparticles, their TPA cross-sections per gold atom are also found to be larger. Theoretical simulations have suggested that the large enhancement of these TPA cross-sections comes from a one-photon double-resonance mechanism. However, it remains difficult to simulate TPA cross-sections of thiolate-protected gold clusters due to their large system size and a high density of states. In this work, we report a time-dependent density functional theory (TDDFT) study of the TPA spectra of the Au25(SR)18 cluster based on a damped response theory formalism. Damped response theory enables a consistent treatment of on- and off-resonance molecular properties even for molecules with a high density of states, and thus is well-suited for studying the TPA properties of gold clusters. Our results indicate that the one- and two-photon double-resonance effect is much smaller than previously found, and thus is unlikely to be the main cause of the large TPA cross-sections found experimentally. The effect of symmetry breaking of the Au25(SR)18 cluster due to the ligands on the TPA cross-sections has been studied and was found to only slightly increase the cross-section. Furthermore, by comparing with larger nanoparticles we find that the TPA cross-section per gold atom scales linearly with the diameter of the particles, and that the Kerr non-linear response of the Au25(SR)18 cluster is on the same order as that of bulk gold films.

1 Introduction

Thiolate-protected gold nanoparticles have attracted significant interest in recent years due to their exceptional stability and applications in biomedicine, catalysis, electronics, photonics and sensing.1–6 While larger gold nanoparticles (>5 nm) are characterized by localized surface plasmon resonance, small gold clusters (<3 nm) exhibit molecular-like properties due to the quantum confinement effects.3 Although the exact boundary between molecular and plasmonic response has not been established, experiments have shown that a small nanoparticle containing only ∼330 gold atoms exhibits a plasmonic response.7 Simulations using time-dependent density functional theory (TDDFT) indicate that the boundary between molecular and plasmonic behaviors occurs for the Au144(SR)60 monolayer-protected cluster with a 1.5 nm core.8

Numerous gold clusters (Aun(SR)m) have been explored both experimentally and theoretically since the total structural determination of the nanoclusters Au102(SR)44 and Au25(SR)18 using X-ray crystallography.9–12 Au25(SR)18 is probably the most extensively studied cluster due to its extraordinary atomic packing structure and the well established structure–property relationship.10–12 The stability of the Au25(SR)18 cluster can be understood by considering it to consist of a Au135+ core surrounded by 6 anionic RS–(Au–SR)2 units. According to the “super-atom” model,13 which is commonly used to understand the stability of gold clusters, this leads to a shell-closing of 8 electrons in the core that strengthens its stability. The existence of a crystal structure of the Au25(SR)18 cluster has enabled a detailed correlation between its structure and optical properties through TDDFT simulations.12,14,15

There is a large amount of literature on the linear optical properties of small gold clusters. This is mainly because such properties are sensitive to the specific atomic arrangement, and thus can aid structural determination. Furthermore, understanding the linear properties of these small clusters provides insights into the emergence of the plasmonic response found in larger nanoparticles. The non-linear optical (NLO) properties of small gold clusters have also attracted attention due to their potential use in multiphoton imaging and optical limiting applications. Although much less work has been devoted to studying the NLO properties of thiolate-protected gold clusters in contrast to their linear counterparts, experimental two-photon absorption (TPA), non-linear transmission, hyper-Rayleigh scattering, and second- and third-harmonic generation measurements have been performed for the prototypical thiolate-protected Au25 cluster.16–21 However, few theoretical studies of the NLO properties of these small ligand-protected gold clusters are available in the literature due to a high computational burden.22–24

The most intriguing NLO property of the Au25(SR)18 cluster is probably its large TPA cross-section in the communication wavelength region, which has been reported by Ramakrishna et al.16 as 2700 GM for excitation at 1290 nm. This value is much larger than that of many organic chromophores, and promotes the Au25(SR)18 cluster as a potentially well-qualified candidate for various NLO applications including biological imaging, nanolithography, and optical limiting.17 More interestingly, a very large TPA cross-section of 427[thin space (1/6-em)]000 GM was also found for this cluster for excitation at 800 nm.16 Furthermore, the TPA cross-sections per gold atom for it and a few other small gold clusters were shown to exhibit a different size-scaling as compared to larger nanoparticles.16 To understand this unusual behavior, Day et al.24 used TDDFT simulations based on the single residue of quadratic response functions to obtain the TPA cross-section of the Au25(SR)18 cluster. These calculations suggested that a one-photon double-resonance effect could lead to the large TPA cross-sections observed experimentally. Recent work has also found that the optical Kerr response in 3 nm gold films is many orders of magnitude larger than that of the bulk metal.25 These results suggest that quantum size effects could lead to a significant third-order non-linear response for small gold clusters. However, simulations of resonance non-linear properties of metal clusters pose a significant challenge due to the high number of states contributing to the spectra, which raises the possibility of resonance effects that need to be dealt with carefully to avoid unphysically large response properties. This is particularly a problem when using standard response theory to calculate the molecular properties since the response functions diverge when optical frequencies, or the sum of them, equal an excitation energy.

In this work, we report simulated TPA spectra for the Au25(SR)18 cluster using TDDFT. To avoid unphysical resonance effects, we will use a recently implemented damped cubic response formalism26 to calculate the TPA cross-sections. Damped response theory26–30 takes the broadening of electronic states into account and thus avoids the unphysical behaviors for molecular properties on resonance. In addition to TPA cross-sections, we will also characterize the resonance optical Kerr effect for the Au25(SR)18 cluster. Both of these optical processes can be described by a third-order non-linear response tensor obtained using damped response theory. Our results show that the one-photon double resonance effect is smaller than previously found. We also find that the quantum size effects for the Au25(SR)18 cluster do not lead to a significantly enhanced third-order non-linear response.

2 Theory

Considering the simultaneous absorbance process for two linearly polarized photons with identical energies,31 one can utilize the imaginary part of orientationally averaged third-order response properties to express the TPA cross-section (σTPA) as:32
image file: c7sc00968b-t1.tif(1)
where αf is the fine structure constant, ω is the incident frequency, and N is an integer value related to the experimental setup.33 In this work, N = 4 is used for all simulated TPA spectra and the unit of σTPA is given as Göppert-Mayer (1 GM = 10−50 cm4 s per photon).34 The term γ(−ω; ω, ω, −ω) is denoted as γIDRI in this work as its real part corresponds to the intensity-dependent refractive index (IDRI). The imaginary part of γIDRI includes both saturated linear absorption and two-photon absorption, and in the traditional sum-over-states (SOS) approach the two processes can be related to the negative and two-photon terms (N-terms and T-terms), respectively.35 The N-terms provide large negative contributions to the γ tensors36,37 and correspond to purely one-photon processes, hence should not be considered when describing TPA. By eliminating the N-terms, one can obtain a reduced form for the IDRI, of which the imaginary part corresponds to the pure TPA process. This reduced IDRI, termed γTPA in this work, can be given using the SOS expression as:
image file: c7sc00968b-t2.tif(2)
where ωn0 is the excitation energy of state n, Γ is the energy broadening parameter, and Sαβ is the TPA transition moment that involves the transition dipole moments between the ground and excited states (〈0|μα|m〉) as well as two excited states (image file: c7sc00968b-t3.tif), i.e.
image file: c7sc00968b-t4.tif(3)

The corresponding expression for γTPA in a damped response formalism has been previously reported by Kristensen et al.38 and Hu et al.26 The use of γTPA avoids any negative TPA intensities caused by the pure one-photon processes and allows for appropriate σTPA calculations in the presence of one- and two-photon double-resonance effects. Therefore, we adopt it for all TPA simulations in this work.

3 Computational details

All calculations in this work were carried out through a locally modified version of the Amsterdam Density Functional (ADF) 2014 program package.39–41 The starting geometry for the Au25(SH)18 cluster is based on the crystal structure,10,12 and the initial atomic coordinates for the Au25(SH)16(SPh)2 cluster were obtained from ref. 22. Geometry optimization was performed using the Becke–Perdew (BP86)42,43 XC functional with a small frozen-core triple-ζ polarized Slater-type (TZP) basis set from the ADF library. The BP86 XC functional with a large frozen-core TZP basis set was adopted to calculate the first- and third-order response properties. Scalar relativistic effects have been accounted for by means of the zeroth-order regular approximation (ZORA).44,45 Solvent effects are not included in the simulations, but good agreement between theory and experiment has been previously12 demonstrated for the absorption spectrum of the Au25(SR)18 cluster. The finite lifetime of the electronic excited states is included phenomenologically using a damping parameter of 0.0034 a.u. (∼0.1 eV), which was previously found to be acceptable29,30 and also roughly the same as the Lorentzian fitting width used in ref. 12. The conversion factor to SI and cgs units for γ is:46 1 a.u. = 7.0423 × 10−54 m5 V−2 = 5.0367 × 10−40 esu.

4 Results and discussion

The one-photon absorption spectrum of the Au25(SH)18 cluster is characterized by three main bands found at 1.8, 2.8, and 3.1 eV, respectively.12 The lowest band is a HOMO to LUMO transition which can be characterized as an intraband transition (sp ← sp), the second band arises from mixed intraband (sp ← sp) and interband (sp ← d) transitions, while the third band arises predominantly from interband transitions (sp ← d).12,14,15,47 While initially the lowest transition at 1.8 eV was described in terms of the electronic and geometric structure of the Au13 core, recent work has shown that the optical absorption spectra are not separable into core and ligand contributions.14 A comparison between the simulated and experimental absorption spectra for the Au25(SR)18 cluster is shown in Fig. 1. The simulated spectrum is obtained in the gas phase using the –SH group as the ligand while the experimental one is measured in toluene for the gold cluster passivated by the SCH2CH2Ph ligands.12
image file: c7sc00968b-f1.tif
Fig. 1 (a) Simulated absorption spectrum for the Au25(SR)18 (R = H) cluster in the gas phase. (b) Experimental absorption spectrum for the Au25(SR)18 (R = CH2CH2Ph) cluster in toluene taken from ref. 12.

The theoretical spectrum consists of a major band at 1.43 eV and a broader band ranging from 2.2 → 3.1 eV. The former one (labeled as “a”) corresponds to the lowest band in the experimental spectrum at 1.8 eV (labeled as “a′”) and is primarily characterized by a HOMO to LUMO transition. The latter one has contributions from three groups of transitions at 2.4, 2.6, and 2.8 eV, labeled as “b”, “c” and “d”, respectively. The “b” band corresponds to the mixed intraband and interband experimental transitions at 2.8 eV (labeled as “b′”), and the “c” band corresponds to the experimental transitions at 3.1 eV (labeled as “c′”). Although the “d” band is not resolved in the experimental spectrum due to thermal broadening and its weak oscillator strength, low temperature measurements of the absorption spectrum have shown several additional bands above 3 eV.48 The red-shift of the simulated spectrum is likely a result of the neglection of solvent effects, the choice of XC functionals, and the simplified ligand used in the simulations.49 The band assignment presented here follows that of ref. 12, where a larger splitting of the “b” and “c” bands was obtained by using the SAOP XC potential. We refer to ref. 15 for a comprehensive discussion of the optical absorption of the thiolate-protected Au25 cluster.

In Fig. 2 we plot the simulated TPA spectrum for the Au25(SH)18 cluster as a function of the one-photon energy. The spectrum is dominated by a broad band at 1.4 eV with a shoulder at 1.3 eV and a weaker band at 1.1 eV. The weak low-energy band corresponds to two-photon excitations into a set of weaker states that are found as a shoulder to the “b” band in the one-photon absorption spectrum. The stronger shoulder at 1.3 eV in the two-photon spectrum corresponds to excitation into the strong “c” band in the one-photon absorption spectrum. Finally, the largest two-photon absorption cross-section is found around 753 GM for excitation into the “d” band at 1.4 eV. This σTPA value is comparable to that of large organic TPA chromophores,50 mainly due to a double-resonance effect where the two-photon transition into the “d” band is enhanced by a one-photon resonance with the “a” band. We demonstrate this double-resonance effect by plotting a schematic energy diagram with the most important transitions in Fig. 3. From eqn (3), we see that the two-photon transition moments S0n become large for photon energies near the one-photon transitions.

image file: c7sc00968b-f2.tif
Fig. 2 Simulated TPA spectrum for the Au25(SH)18 clusters in the gas phase.

image file: c7sc00968b-f3.tif
Fig. 3 Schematic energy diagram showing the most important transitions for the Au25(SH)18 cluster in the gas phase.

Experimentally, a huge TPA cross-section of 427[thin space (1/6-em)]000 GM was found for a photon energy of 1.55 eV.16 Since this corresponds to excitation into the “a” band in the one-photon absorption spectrum, the huge TPA cross-section likely results from enhancement by one-photon processes, either a double-resonance effect or other effects such as excited state absorption. Previous simulations have found very large TPA cross-sections for small monolayer protected gold and silver clusters due to double-resonance effects.24,51,52 A TPA cross-section of 620[thin space (1/6-em)]000 GM at a photon energy of 1.58 eV was reported for the Au25(SH)18 cluster using B3LYP and a SDD-DZ basis set, with similar values found for the other functionals tested.24 The large TPA cross-section is in good agreement with the experimental observation and was attributed to resonance enhancement from the lowest excited state. However, in contrast to these huge TPA cross-sections found in the previous simulations, our results are several orders of magnitude smaller with the largest value found around 753 GM. Previous simulations used traditional quadratic response theory within a TDDFT formalism to simulate the TPA cross-section, which can produce unphysical large TPA cross-sections due to double-resonance effects.26,38,53 A major advantage of using damped response theory is that the two-photon transition moments remain finite even if in the vicinity of one-photon resonances, and thus can correctly describe this double-resonance effect.24

Experimentally, it was also found that the TPA cross-sections per gold atom for the small gold clusters were significantly larger than those for larger nanoparticles, following a different size dependence.16 This was attributed to quantum size effects and indicated a transition from small clusters characterized in terms of discrete transitions to larger nanoparticles characterized in terms of plasmon resonances. For the Au25(SH)18 cluster, a TPA cross-section per gold atom of 17[thin space (1/6-em)]080 GM was reported16 while we find a much smaller value of 30 GM per gold atom. However, our results compare well with the expected value based on the size-scaling observed for the TPA cross-sections of the larger nanoparticles. This is illustrated in the ESI, where we plot the experimental σTPA/gold as a function of the size of the Au976 (3.0 nm) and the Au2406 (4.0 nm) nanoparticles16 compared with the simulated results for the Au25(SH)18 cluster.

The experiments also reported a large TPA cross-section of 2700 GM for a photon energy of 0.96 eV, corresponding to a two-photon excitation into the “a′” band.16 As this is the lowest one-photon band, the possible double-resonance effects can be ruled out. Considering the fact that the orientation of the ligands with respect to the Au25S18 core makes the entire cluster roughly centrosymmetric,10–12 one should expect the lowest excited state to be one-photon allowed but two-photon forbidden. This is consistent with our simulations, where no significant TPA cross-section is found for excitation into the “a” band. This is also in agreement with the results reported by Day et al.24 The large TPA cross-section at 0.96 eV was speculated by Day et al.24 as the result of the tail of the very large TPA cross-section for excitation into the “d” band at 1.4 eV. However, that is not seen in our simulations.

Alternatively, symmetry breaking of the ligand shell could lead to an increased TPA cross-section. Previous work has demonstrated that such symmetry breaking leads to the observation of hyper-Rayleigh scattering of the Au25(SR)18 cluster which otherwise would be symmetry forbidden due the centrosymmetry.19,22 To investigate the symmetry breaking effects on the TPA cross-section, we also considered the Au25(SH)16(SPh)2 cluster where two phenyl groups have been substituted. This structure is taken from previous work investigating the symmetry breaking effects on the first hyperpolarizabilities.22 The TPA spectrum for the Au25(SH)16(SPh)2 cluster is shown in Fig. 4. In comparison to the Au25(SH)18 cluster, the addition of the two phenyl groups leads to slightly larger TPA cross-sections with a maximum of 905 GM. However, no significant TPA cross-section is found for excitation into the “a” band, and thus symmetry breaking is unlikely to be the reason for the large TPA cross-section observed experimentally.

image file: c7sc00968b-f4.tif
Fig. 4 Simulated TPA spectrum for the Au25(SH)16(SPh)2 cluster in the gas phase.

To investigate the importance of other one-photon resonance enhancements in the third-order non-linearity of Au25(SH)18, we simulated the optical Kerr response corresponding to the γ(−ω; ω, ω, −ω) tensor with all one-photon terms included. The real part of the optical Kerr response is related to the IDRI and its imaginary part corresponds to two-photon and saturated one-photon absorption processes. In Fig. 5 we plot both Re[γIDRI] and Im[γIDRI] for the Au25(SH)18 cluster as a function of the one-photon energy, together with the contribution arising only from two-photon absorption, i.e. Im[γTPA]. The Re[γIDRI] curve is characterized by a large positive band at 1.43 eV and a small negative band at 1.31 eV. In contrast, a large negative band at 1.37 eV with a small positive band at 1.46 eV mainly constitutes the Im[γIDRI] curve. Im[γTPA] is much larger than Im[γIDRI] demonstrating a strong destructive interference between the one- and two-photon processes. This seems to indicate that the one-photon processes are more dominant compared to the two-photon processes, and could be the reason that only relatively modest TPA cross-sections are found for the Au25(SH)18 cluster. The magnitude of the optical Kerr response determined here is comparable to that of the pure electronic response in organic molecules.54 Although there are no direct measurements of the optical Kerr response of small gold clusters, Qian et al.25 reported the optical Kerr response of 3 nm gold films. They found a large Kerr susceptibility of 2.06 × 10−15 m2 V−2, which is about four orders of magnitude larger than that determined for bulk gold films. The large non-linear response observed was attributed to quantum size effects in the thin gold film. Here we calculated |χ3| to be 3.26 × 10−19 m2 V−2 for the Au25(SH)18 cluster by assuming a diameter of ∼1.1 nm. Our calculated value is comparable to the bulk value of 9.11 × 10−19 m2 V−2 for a 15 nm thick gold film reported in the same work, but much smaller than the thin film result.25 Thus, in general we do not find that quantum size effects lead to a significantly enhanced non-linear response in the small monolayer protected clusters.

image file: c7sc00968b-f5.tif
Fig. 5 Calculated Re[γIDRI] (blue solid line), Im[γTPA] (green dashed line), and Im[γIDRI] (magenta dash-dotted line) for the Au25(SH)18 cluster in the gas phase. The black dashed line shows a value of 0 esu.

Although our simulations show that the TPA cross-sections for the Au25(SH)18 cluster are enhanced by a double-resonance effect, they are much smaller than what has been indicated by previous simulations or found experimentally. In general, our results seem to suggest that quantum size effects in these small Au25(SR)18 clusters do not lead to extremely large TPA cross-sections. However, it is important to note that the results presented here are sensitive to the exact value used for the energy broadening parameter (Γ). In this work, Γ was chosen to match the experimental absorption spectrum. Using a smaller Γ value would lead to larger TPA cross-sections, but one would need unphysical small values to match the experimental TPA cross-sections. Finally, the TDDFT simulations presented here all employ an adiabatic approximation for the exchange–correlation kernels. Previously, it was shown that such a typical approximation could cause spurious pole effects near one-photon resonances.26 For small molecules, these spurious poles lead to significantly larger TPA cross-sections, yet the behavior is still unknown for systems like the Au25(SR)18 clusters that have a high density of states.

5 Conclusions

In summary, we have reported first-principles simulations of TPA spectra for two thiolate-protected Au25 clusters based on a damped cubic response formalism within TDDFT. We find that the calculated TPA cross-sections are much smaller than their experimental counterparts, which indicates that the previously suggested one- and two-photon double resonance effect is unlikely to be the only cause for the large TPA intensities reported experimentally. The calculation of TPA cross-sections on a per gold atom basis, as well as the Kerr non-linear responses, is in-line with those expected from larger nanoparticles. Symmetry breaking was shown to only lead to small enhancements of the TPA cross-sections. Overall, this work represents the first cubic response theory approach to TPA simulations of the Au25(SR)18 clusters, and shows that quantum size effects do not lead to a significantly enhanced third-order non-linear response. This is in agreement with the molecular origin of the electronic transitions in the small gold clusters.


This work was supported by the NSF award CHE-1362825. We acknowledge support received from the Research Computing and Cyber infrastructure, a unit of the Information Technology Services at Penn State.


  1. M.-C. Daniel and D. Astruc, Chem. Rev., 2004, 104, 293 CrossRef CAS PubMed.
  2. H. Häkkinen, Chem. Soc. Rev., 2008, 37, 1847–1859 RSC.
  3. R. Jin, Nanoscale, 2010, 2, 343 RSC.
  4. E. C. Dreaden, A. M. Alkilany, X. Huang, C. J. Murphy and M. A. El-Sayed, Chem. Soc. Rev., 2012, 41, 2740–2779 RSC.
  5. A. Fernando, K. L. D. M. Weerawardene, N. V. Karimova and C. M. Aikens, Chem. Rev., 2015, 115, 6112–6216 CrossRef CAS PubMed.
  6. R. Jin, C. Zeng, M. Zhou and Y. Chen, Chem. Rev., 2016, 116, 10346–10413 CrossRef CAS PubMed.
  7. H. Qian, Y. Zhu and R. Jin, Proc. Natl. Acad. Sci. U. S. A., 2012, 109, 696–700 CrossRef CAS PubMed.
  8. S. Malola, L. Lehtovaara, J. Enkovaara and H. Häkkinen, ACS Nano, 2013, 7, 10263–10270 CrossRef CAS PubMed.
  9. P. D. Jadzinsky, G. Calero, C. J. Ackerson, D. A. Bushnell and R. D. Kornberg, Science, 2007, 318, 430–433 CrossRef CAS PubMed.
  10. M. W. Heaven, A. Dass, P. S. White, K. M. Holt and R. W. Murray, J. Am. Chem. Soc., 2008, 130, 3754 CrossRef CAS PubMed.
  11. J. Akola, M. Walter, R. L. Whetten, H. Häkkinen and H. Grönbeck, J. Am. Chem. Soc., 2008, 130, 3756 CrossRef CAS PubMed.
  12. M. Zhu, C. M. Aikens, F. J. Hollander, G. C. Schatz and R. Jin, J. Am. Chem. Soc., 2008, 130, 5883 CrossRef CAS PubMed.
  13. M. Walter, J. Akola, O. Lopez-Acevedo, P. D. Jadzinsky, G. Calero, C. J. Ackerson, R. L. Whetten, H. Grönbeck and H. Häkkinen, Proc. Natl. Acad. Sci. U. S. A., 2008, 105, 9157 CrossRef CAS PubMed.
  14. C. M. Aikens, J. Phys. Chem. C, 2008, 112, 19797 CAS.
  15. C. M. Aikens, J. Phys. Chem. A, 2009, 113, 10811 CrossRef CAS PubMed.
  16. G. Ramakrishna, O. Varnavskit, J. Kim, D. Lee and T. Goodson, J. Am. Chem. Soc., 2008, 130, 5032 CrossRef CAS PubMed.
  17. S. H. Yau, O. Varnavski and I. Theodore Goodson, Acc. Chem. Res., 2013, 46, 1506–1516 CrossRef CAS PubMed.
  18. R. Philip, P. Chantharasupawong, H. Qian, R. Jin and J. Thomas, Nano Lett., 2012, 12, 4661–4667 CrossRef CAS PubMed.
  19. I. Russier-Antoine, F. Bertorelle, M. Vojkovic, D. Rayane, E. Salmon, C. Jonin, P. Dugourd, R. Antoine and P.-F. Brevet, Nanoscale, 2014, 6, 13572–13578 RSC.
  20. S. Knoppe, M. Vanbel, S. van Cleuvenbergen, L. Vanpraet, T. Bürgi and T. Verbiest, J. Phys. Chem. C, 2015, 119, 6221–6226 CAS.
  21. J. Olesiak-Banska, M. Waszkielewicz, K. Matczyszyn and M. Samoc, RSC Adv., 2016, 6, 98748–98752 RSC.
  22. S. Knoppe, H. Häkkinen and T. Verbiest, J. Phys. Chem. C, 2015, 119, 27676–27682 CAS.
  23. N. Van Steerteghem, S. Van Cleuvenbergen, S. Deckers, C. Kumara, A. Dass, H. Häkkinen, K. Clays, T. Verbiest and S. Knoppe, Nanoscale, 2016, 8, 12123–12127 RSC.
  24. P. N. Day, K. A. Nguyen and R. Pachter, J. Chem. Theory Comput., 2010, 6, 2809 CrossRef CAS PubMed.
  25. H. Qian, Y. Xiao and Z. Liu, Nat. Commun., 2016, 7, 13153 CrossRef CAS PubMed.
  26. Z. Hu, J. Autschbach and L. Jensen, J. Chem. Theory Comput., 2016, 12, 1294–1304 CrossRef CAS PubMed.
  27. P. Norman, D. M. Bishop, H. J. A. Jensen and J. Oddershede, J. Chem. Phys., 2001, 115, 10323–10334 CrossRef CAS.
  28. P. Norman, D. M. Bishop, H. J. A. Jensen and J. Oddershede, J. Chem. Phys., 2005, 123, 194103 CrossRef PubMed.
  29. L. Jensen, J. Autschbach and G. C. Schatz, J. Chem. Phys., 2005, 122, 224115 CrossRef CAS PubMed.
  30. Z. Hu, J. Autschbach and L. Jensen, J. Chem. Phys., 2014, 141, 124305 CrossRef PubMed.
  31. W. M. McClain, J. Chem. Phys., 1971, 55, 2789 CrossRef.
  32. D. W. Silverstein and L. Jensen, J. Chem. Phys., 2012, 136, 064111 CrossRef PubMed.
  33. M. T. P. Beerepoot, D. H. Friese, N. H. List, J. Kongsted and K. Ruud, Phys. Chem. Chem. Phys., 2015, 17, 19306–19314 RSC.
  34. M. Göppert-Mayer, Ann. Phys., 1931, 401, 273 CrossRef.
  35. T. R. Ensley, H. Hu, M. Reichert, M. R. Ferdinandus, D. Peceli, J. M. Hales, J. W. Perry, Z. Li, S.-H. Jang, A. K.-Y. Jen, S. R. Marder, D. J. Hagan and E. W. V. Stryland, J. Opt. Soc. Am. B, 2016, 33, 780–796 CrossRef CAS.
  36. J. Pérez-Moreno and M. G. Kuzyk, J. Chem. Phys., 2005, 123, 194101 CrossRef PubMed.
  37. J. Pérez-Moreno, K. Clays and M. G. Kuzyk, J. Chem. Phys., 2008, 128, 084109 CrossRef PubMed.
  38. K. Kristensen, J. Kauczor, A. J. Thorvaldsen, P. Jørgensen, T. Kjærgaard and A. Rizzo, J. Chem. Phys., 2011, 134, 214104 CrossRef PubMed.
  39. G. te Velde, F. M. Bickelhaupt, E. J. Baerends, C. Fonseca Guerra, S. J. A. van Gisbergen, J. G. Snijders and T. Ziegler, J. Comput. Chem., 2001, 22, 931–967 CrossRef CAS.
  40. C. Fonseca Guerra, J. G. Snijders, G. te Velde and E. J. Baerends, Theor. Chem. Acc., 1998, 99, 391–403 Search PubMed.
  41. E. Baerends, T. Ziegler, J. Autschbach, D. Bashford, A. Bérces, F. Bickelhaupt, C. Bo, P. Boerrigter, L. Cavallo, D. Chong, L. Deng, R. Dickson, D. Ellis, M. van Faassen, L. Fan, T. Fischer, C. F. Guerra, M. Franchini, A. Ghysels, A. Giammona, S. van Gisbergen, A. Götz, J. Groeneveld, O. Gritsenko, M. Grüning, S. Gusarov, F. Harris, P. van den Hoek, C. Jacob, H. Jacobsen, L. Jensen, J. Kaminski, G. van Kessel, F. Kootstra, A. Kovalenko, M. Krykunov, E. van Lenthe, D. McCormack, A. Michalak, M. Mitoraj, S. Morton, J. Neugebauer, V. Nicu, L. Noodleman, V. Osinga, S. Patchkovskii, M. Pavanello, P. Philipsen, D. Post, C. Pye, W. Ravenek, J. Rodríguez, P. Ros, P. Schipper, H. van Schoot, G. Schreckenbach, J. Seldenthuis, M. Seth, J. Snijders, M. Solà, M. Swart, D. Swerhone, G. te Velde, P. Vernooijs, L. Versluis, L. Visscher, O. Visser, F. Wang, T. Wesolowski, E. van Wezenbeek, G. Wiesenekker, S. Wolff, T. Woo and A. Yakovlev, Amsterdam Density Functional, 2014 Search PubMed.
  42. A. D. Becke, Phys. Rev. A, 1988, 38, 3098–3100 CrossRef CAS.
  43. J. P. Perdew, Phys. Rev. B: Condens. Matter Mater. Phys., 1986, 33, 8822–8824 CrossRef.
  44. E. van Lenthe, E. J. Baerends and J. G. Snijders, J. Chem. Phys., 1993, 99, 4597–4610 CrossRef CAS.
  45. E. van Lenthe, E. J. Baerends and J. G. Snijders, J. Chem. Phys., 1994, 101, 9783–9792 CrossRef CAS.
  46. D. P. Shelton and J. E. Rice, Chem. Rev., 1994, 94, 3–29 CrossRef CAS.
  47. C. M. Aikens, J. Phys. Chem. Lett., 2011, 2, 99–104 CrossRef CAS PubMed.
  48. M. S. Devadas, S. Bairu, H. Qian, E. Sinn, R. Jin and G. Ramakrishna, J. Phys. Chem. Lett., 2011, 2, 2752–2758 CrossRef CAS.
  49. H.-F. Zhang, M. Stender, R. Zhang, C. Wang, J. Li and L.-S. Wang, J. Phys. Chem. B, 2004, 108, 12259–12263 CrossRef CAS.
  50. F. Terenziani, C. Katan, E. Badaeva, S. Tretiak and M. Blanchard-Desce, Adv. Mater., 2008, 20, 4641–4678 CrossRef CAS.
  51. P. N. Day, R. Pachter, K. A. Nguyen and T. P. Bigioni, J. Phys. Chem. A, 2016, 120, 507–518 CrossRef CAS PubMed.
  52. Ž. Sanader, M. Krstić, I. Russier-Antoine, F. Bertorelle, P. Dugourd, P.-F. Brevet, R. Antoine and V. Bonačić-Koutecký, Phys. Chem. Chem. Phys., 2016, 18, 12404–12408 RSC.
  53. S. M. Parker, S. Roy and F. Furche, J. Chem. Phys., 2016, 145, 134105 CrossRef PubMed.
  54. S. R. Vigil and M. G. Kuzyk, J. Opt. Soc. Am. B, 2001, 18, 679–691 CrossRef CAS.


Electronic supplementary information (ESI) available. See DOI: 10.1039/c7sc00968b

This journal is © The Royal Society of Chemistry 2017