# Effects of Interference between Energy Absorption Processes of Molecule and Surface Plasmons on Light Emission Induced by Scanning Tunneling Microscopy

\abstEffects of coupling between an molecular exciton and a surface plasmon (exciton-plasmon coupling) on the luminescence properties of the molecule and the surface plasmons are investigated using nonequilibrium Green’s function method. Molecular absorption and enhancement by molecular electronic and vibrational modes (molecular modes) lead to dip and peak structures in the luminescence spectra of the surface plasmons. The re-absorption by the surface plasmons leads to a dent structure in their luminescence spectra. It is found that the processes of the molecular absorption and the re-absorption by the surface plasmons interfere with each other. Enhancement and suppression of these energy-absorption processes are due to the interference. Moreover, due to the exciton-plasmon coupling, excitation channels of the molecule arise in the energy range lower than the first electronic excitation energy of the molecule. It is found that the electron transitions via these excitation channels give rise to the molecular luminescence and the vibrational excitations at the bias voltage lower than the ratio of the first electronic excitation energy of the molecule to the elementary charge. The results also indicate that the vibrational excitations assist the emission of photons, whose energy exceeds the product of the elementary charge and the bias voltage (upconverted luminescence). \kwordmolecular luminescence, surface plasmon, scanning tunneling microscopy, molecular vibration, exciton-plasmon coupling, nonequilibrium Green’s function method, interference

## 1 Introduction

The optical property of the system, where electronically excited molecules and surface plasmons interact with each other, has been a topic of great interest in various research fields. Surface plasmons localized near metal nanostrcutures generate an intense electromagnetic field in a nano scale region. Therefore, field enhancement by the surface plasmons are commonly employed in many applications, such as surface-enhanced Raman spectroscopy, surface-enhanced infrared absorption spectroscopy, high-resolution biosensors, light-emitting devices, photovoltaics, and so forth. [1, 2, 3, 4, 5, 6, 7, 8, 9, 10] These pioneering works have mainly focused on how an electromagnetic field generated by surface plasmons affects the optical properties of the target molecules. In addition to this, recent studies suggest that dynamics of the molecules including luminescence and absorption can also affect the optical property of the surface plasmons. For example, coupling of the surface plasmons to a molecular exciton (exciton-plasmon coupling) leads to structures (e.g., peak, dip, and Fano-like structures) in luminescence and absorption spectra of the surface plasmons. [11, 12, 13, 14] Microscopic understanding of effects of the exciton-plasmon coupling will help in elucidating the true nature of the energy transfer between the metal nanostructures and the molecules.

The highly localized tunneling current of a scanning tunneling microscope (STM) can be utilized as an atomic-scale source to induce light emission from the system. STM-induced light emission (STM-LE) spectroscopy has a unique advantage of having the potential to unveil optical properties with submolecular spatial resolution. Moreover, one of the most intersting things regarding STM-LE from clean and molecule-covered metal substrates is that the surface plasmons localized near the tip-sample gap region contribute to luminescence processes.

Gimzewski et al. has first reported STM-LE from a metal surface. [15] Nowadays, it is well confirmed that light emission from clean metal surfaces is due to the radiative decay of surface plasmons localized near the tip-sample gap region. [16, 17, 18, 19, 20] When a molecule is located near the gap region, there are two radiative processes, i.e., luminescence from the surface plasmons and the molecule. If the molecule is directly adsorbed on a metal substrate, charge and energy transfer between them lead to the quenching of molecular luminescence. Thus, the molecule acts as a spacer between the tip and the metal substrate, which slightly modifies the plasmon-mediated emission. [21, 22, 23] If the molecule is electrically decoupled from the metal substrate by dielectric films or molecular multilayers, intrinsic molecular luminescence, which is associated with the electronic and vibrational transitions in the molecule, can be observed. [24, 25, 26, 27]

Recent studies have reported the importance of effects of the interaction between the intense electromagnetic field generated by the surface plasmons and the transition moments for the molecular excitations and de-excitations. Due to the fact that the molecular luminescence can be in the same energy window as the plasmon-mediated emission, light emission can be induced by the two kinds of radiative processes simultaneously. Uemura et al. and Liu et al. have shown that the plasmon-mediated emission overlapping with the molecular absorption enhances molecular luminescence intensity. [28, 29, 30] Dong et al. have reported selective enhancement of luminescence intensity by spectrally tuning the surface plasmon mode to match a particular transition in the molecule. [31] Moreover, fascinating new phenomena, i.e., luminescence associated with radiative decay from excited vibrational levels for the excited electronic state of the molecule (hot luminescence) and luminescence at energies beyond the product of the elementary charge and the applied bias voltage (upconverted luminescence) have also been observed. [31]

Regarding the excitation processes of the molecule, there are two dominant mechanisms involved. One is the excitation by the injection of electrons and/or holes from the electrodes. The other is the excitation by the absorption of surface plasmons that are excited by the inelastic tunneling between the tip and the substrate. Recent experiments suggest the importance of the plasmon-mediated excitations [28, 29, 30, 31, 32]. In surface-plasmon fluorescence spectroscopy, it has reported that the molecule is excited by surface plasmons. [33] Furthermore, a similar plasmon-mediated excitation of molecules has been proposed in previous theoretical studies. In these studies, molecular dynamics have been investigated within the framework of the one-body problem, where the surface plasmon is treated as a classical electric field. [34, 35, 36, 37] Tian et al. [36, 37] used the plasmon-mediated excitation mechanism to explain some spectral features observed in Ref. \citenDong2010. The results indicate that the molecular dynamics induced by the surface plasmons play essential roles.

Direct experimental evidences of the plasmon-mediated excitations of molecules have been obtained in luminescence spectra acquired with tips with a molecule-covered shaft over clean metal surface. [13] Although no electron tunneling to the molecules takes place, the observed spectra, that can be considered as the luminescence spectra of the surface plasmons, show peak and dip structures. The position of some of these structures matches the position of the peak structures in molecular luminescence and absorption spectra. Thus, the molecular dynamics including the molecular luminescence and absorption can be expected to have an impact on luminescence-spectral profiles of the surface plasmons. To understand this from a microscopic point of view, there is a need to investigate the interplay between the dynamics of the molecule and the surface plasmons within the framework of the quantum many-body theory.

We have recently investigated the effects of coupling between a molecular exciton and the surface plasmon (exciton-plasmon coupling) on the luminescence properties with the aid of nonequilibrium Green’s function method. [38] The results have first revealed that in addition to the molecular luminescence and absorption, the re-absorption by the surface plasmons modifies the luminescence-spectral profiles of the surface plasmons. It has also been shown that the electron transitions, which give rise to the molecular luminescence and the vibrational excitations, occur at the bias voltage lower than the ratio of the energy of the molecular exciton to the elementary charge . [39]

In this paper, we provide the details on the calculation procedure and the analysis processes. In addition, since the processes of the energy absorption by the molecule and the surface plasmons are found to interfere with each other, the effects of the interference on the luminescence-spectral profiles are discussed. The mechanisms for occurrence of the electronic excitations of the molecule at are also discussed in detail.

## 2 Model

We consider the model, where the molecular exciton interacts with a vibron (a molecular-vibrational quantum) and with the surface plasmon. For simplicity, the molecular vibrations are approximated by a single harmonic oscillator. [37] In the present study, three vibrational levels (ground, first-excited and second-excited vibrational levels) are mainly populated. Therefore the harmonic oscillator model will provide valuable insight into the essential physical features. [40, 41, 42] The surface plasmon modes localized near the tip-sample gap region are approximated by a single energy mode. [37, 36, 43, 16] Plasmon mediated emissions acquired with W-tip on clean Au(111) surface and one- or two-monolayers of TPP molecules on Au(111) surface exhibit a single, broad peak in energy range considered in this study. [31] Hence higher-order modes, that result in peaks at higher energy range, are ignored. Creation (annihilation) of the molecular exciton is induced by absorption (emission) of the surface plasmon. [37, 36, 13, 33] Consequently, the Hamiltonian of the system is

(1) |

where and are creation and annihilation operators (for later convenience assumed to satisfy boson commutation rules) for the molecular exciton with energy . Operators and are creation and annihilation operators for the molecular vibrational mode with energy , and are for a surface plasmon mode with energy , and are for a phonon mode in a thermal phonon bath, and . The energy parameters , and correspond to the dipolar exciton-vibron coupling, the exciton-plasmon coupling, and the coupling between the vibron and the phonon in the thermal phonon bath. The model system Eq. (1) is a generalized one that can be utilized to investigate STM-LE from molecular layers [31] and from a clean metal surface with a molecule-covered tip. [13] Since the parameter values are expected to be of the same order in the situation where the molecules are on molecular layers or adsorbed on a tip surface, the same model system Eq. (1) is applied in both cases.

By applying the Lang-Firsov (LF) canonical transformation, [44] with , becomes

(2) |

where and . Here, the term including is neglected. This approximation is expected to be valid, since the population of the molecular exciton is low .

The luminescence spectra of the molecule and the surface plasmons are obtained by the relations, [45]

(3) | ||||

(4) |

Here and are the lesser projections of the Green’s function of the molecular exciton and the surface plasmon on the Keldysh contour ( and ), respectively. In STM-LE from molecular layers and from a clean metal surface with a molecule-covered tip, quenching of molecular luminescence is suppressed and relatively significant, respectively. When the quenching is suppressed, is expected to be detected experimentally. [25] When the quenching is relatively significant, the luminescence spectra detected experimentally are expected to correspond to . [46]

To investigate the excitation properties and analyze the structures in the luminescence spectra, the spectral functions of the molecule and the surface plasmons, which correspond to their absorption spectra, are calculated by the relations,

(5) | ||||

(6) |

where and are the retarded projections of and .

To obtain the spectral functions and the luminescence spectra, it is needed to calculate the LF transformed and , which are defined as

(7) | ||||

(8) |

where denotes statistical average in representation by system evolution for . The variable is the Keldysh contour time variable. The operator denotes the time ordering along the Keldysh contour. [47] The explicit appearance of the operators and in the LF transformed Green’s function of the molecular exciton (polarization propagator) Eq. (7) arise from the dipolar exciton coupling to the vibron. This is in contrast to the case of equivalent monopole electron and hole coupling to local boson studied in Ref. \citenGumhalter1984, Sec. 5C, in which these terms are absent from the polarization propagator.

To calculate the Green’s function containing the exciton-plasmon coupling , we take into account a random-phase-approximation (RPA)-type diagrammatic series. This approximation is expected to hold for . [49] The integral equation for is given by

(9) |

where

(10) | ||||

(11) | ||||

(12) |

and is the plasmon Green’s function for . The diagrammatic representation of Eq. (9) is given in Fig. 1.

Here, we employ an approximation of decoupling electronic and vibrational parts in , and calculate and . The integral equation for is given by

(13) | ||||

(14) |

where is for . The diagrammatic representation of Eq. (13) is shown in Fig. 2.

The integral equation for the plasmon Green’s function is

(15) |

As for the correlation function , it is expressed in terms of the vibrational momentum Green’s function , which is defined by

(16) |

where . For the second order cumulant expansion of , this connection is given by

(17) |

The integral equation for is given by

(18) | ||||

(19) |

Here, is the vibrational momentum Green’s function for . It contains the self-energy emerged from the coupling to the thermal phonon bath . [44]

By assuming the condition of stationary current, the distribution function of the surface plasmons excited by inelastic tunneling between the tip and the substrate is given by

(20) |

where is a coefficient related to the current amplitude due to the inelastic tunneling, [50] is the elementary charge, and is the bias voltage. The lesser and greater projections of are obtained through the relations, and , where is the retarded projection of . [47] The equations mentioned above are self-consistently solved as described in the next section.

## 3 Self-Consistent Calculation Scheme

Due to the connection between the correlation function and the vibrational momentum Green’s function , represented by Eq. (17), it is easier to express the lesser and greater projections of in terms of in the time domain. At the same time, expressions for the lesser and greater projections of , , , and are easier to be written in the energy domain. Hence, we are working in both time and energy domains, by implementing a fast Fourier transform between them. The ensuring self-consistent calculation scheme consists of the following steps:

1. We start with the retarded projection of the Green’s functions in the energy domain for

(22) |

where represents the vibrational damping due to , which is given by

(23) |

and represents the decay of the surface plasmons for . The lesser and greater projections of these Green’s functions are obtained by

(24) | ||||

(25) | ||||

(26) | ||||

(27) |

where

(28) |

Here, is the Bose-Einstein distribution function. The lesser and greater projections are then transformed to the time domain by the Fourier transformation;

(29) |

where represent the lesser and greater projections of the Green’s functions , and is the real time variable.

2. By utilizing the Langreth theorem [51] and using in lesser and greater projections of Eq.(17), the lesser and greater correlation functions for can be obtained with [44]

(30) | ||||

(31) |

In the first iteration step, is applied.

3. The obtained are convoluted with to obtain according to the relations

(32) | ||||

(33) |

The lesser and greater projections of are then transformed to the energy domain by the inverse Fourier transformation. They are used for calculating the retarded projection of .

4. The self-energies , obtained in the previous step, are used for calculating (retarded, lesser and greater projections) in the energy domain. The retarded projection of are

(34) |

Similarly, by utilizing the Langreth theorem and the inverse Fourier transformation, the lesser and greater projections of are obtained in the energy domain. Using and in lesser and greater projections of Eq.(10) and the Langreth theorem, the lesser and greater projections of can be obtained as

(35) | ||||

(36) |

These functions are then transformed to the energy domain and used for calculating the retarded projection of .

5. The obtained and are used for calculating , according to Eq.(9). The retarded projection of is

(37) |

Utilizing the Langreth theorem and the inverse Fourier transformation, the lesser and greater projections of become

(38) | ||||

(39) |

The obtained are used for calculating , according to Eq. (15). The retarded, lesser and greater projections of are

(40) |

(41) | ||||

(42) |

In this step, and are transformed to the time domain. They are substituted into the lesser and greater projections of Eq.(19) to obtain . They are then transformed to the energy domain and used for calculating .

6. Then is updated by Eq. (18). The retarded, lesser and greater projections of are

(43) |

(44) | ||||

(45) |

They are used again in step 2. Hence, , , , , , and are updated in the same way. The above procedure is repeated until the proportion of the absolute change of the population of the molecular exciton in two subsequent iterations to the latest is less than a predefined threshold of . Here, is given by

(46) |

After convergence, the spectral functions and the luminescence spectra are calculated. To analyze the luminescence spectra, the lesser projections of the right in Eqs. (9) and (15) are divided into four terms using the Langreth theorem.

(47) |

(48) |

The parameters given correspond to the experiment on the STM-LE from tetraphenylporphyrin(TPP) molecules on gold surface: [31] , and . The statistical average is taken for temperature . [31] The square of is reported to be 0.61 on the basis of first-principles calculations. [37] The parameter is given so that the molecular vibrational lifetime due to the coupling to the thermal phonon bath is 13 ps. [31] A Markovian decay is assumed for the surface plasmon to give a plasmon lifetime for of 4.7 fs. [31, 37] The coefficient is set to , where a tunneling current is and an excitation probability of the surface plasmons per electron tunneling event is considered to be .

## 4 Results and Discussion

### 4.1 Plasmon luminescence

Figure 3 shows the luminescence spectra of the surface plasmons for various bias voltages .

At , shows a smooth curve. At , complicated peak and dip structures appear.

The position of the peak structures in the energy range (near 1.71 eV and 1.55 eV) correspond to that of the peak structures in the luminescence spectra of the molecule [Fig. 4(b)]. The position of the dip structures in the energy range (near 1.88 eV, 2.03 eV and 2.19 eV) correspond to that of the peak structures in the spectral function of the molecule [Fig. 4(c)]. A dent structure appears near the energy of the surface plasmon mode (2.05 eV), where the spectral function of the surface plasmons has its maximum intensity [Fig. 4(d)].

To analyze the origin of these structures in , contributions from each term in Eq. (48) are shown in Fig. 5.

The contribution from the first term on the right side of Eq. (48) (green dashed line) represents the plasmon luminescence where the exciton-plasmon coupling is ignored, i.e., . The third term (magenta dotted line) expresses the contribution of the energy transfer from the molecular exciton. The luminescence spectra of the surface plasmon (red line) shows peak structures near 1.71 eV and 1.55 eV. These peak structures are, therefore, attributed to enhancement by the molecular electronic and vibrational modes (molecular modes). The second and fourth terms (magenta dotted and cyan dashed-dotted lines) give the processes where the energy of the surface plasmons are absorbed by the molecule and the surface plasmons (molecular absorption and re-absorption by the surface plasmons, respectively). For the contribution from the second and fourth terms, the dip structures appear near 1.87 eV, 2.03 eV and 2.19 eV in . Since the spectral intensity in the spectral function of the surface plasmons is high in the energy range from 1.85 eV to 2.2 eV, the re-absorption by the surface plasmons leads to the dent structure in this energy range.

However, the origin of the structure near 1.81 eV in (red line in Fig. 5), which looks like a peak structure, cannot be explained by each of these processes. There is a need to consider a superposition of the processes of the molecular absorption and the re-absorption by the surface plasmons. To analyze this, the real and imaginary parts of and are shown in Fig. 6.

The product of and gives the fourth term on the right side of Eq. (48). As is well known, the lesser projection is a pure imaginary number in the energy domain. [47] The function gives the spectral function of the surface plasmons for , and gives the spectral function of the molecule. These functions correspond to their absorption spectra.

One of the possible mechanisms for the occurrence of the structure near 1.81 eV in is an interference between the processes of the molecular absorption and the re-absorption by the surface plasmons. To investigate the interference effects, the real parts of and , the product of which leads to the interference between these processes, are analyzed. Since the real and imaginary parts of satisfy the Kramers-Kronig relations, the real part of is given by , where and is the principal value. Therefore, the real part of is negative for , and positive for [red line in Fig. 6(a)].

Regarding , the processes of the molecular absorption are classified into two kinds: the first one is accompanied by the vibrational transitions while the other is not (inelastic and elastic processes, respectively). From the spectral function of the molecule , it is confirmed that for , the elastic process is dominant over the inelastic processes. [52] Therefore the real part of is negative when the energy is lower than the energy for the elastic process (1.87 eV), where the shift of the peak position from is scaled as . [38] The real part, , tends to be positive in the energy range [red line in Fig. 6(b)].

Therefore, the product gives a positive contribution to the molecular luminescence intensity when is lower than both and . Thus, the interference between the processes of the molecular absorption and the re-absorption by the surface plasmons suppresses these energy absorption processes in this energy range. The product tends to be negative when is between and . Hence, the interference tends to enhance the energy absorption processes in this energy range. When the energy is higher than both and , the interference tends to suppress the energy absorption processes. This trend can be seen even when the energy of the surface plasmon mode is varied. Thus, the structure near 1.81 eV in (red line in Fig. 5), which looks like a peak structure, can be attributed to the interference between these absorption processes.

The peak and dip structures in can be seen in the recent experimental result. [13] Moreover, our results show that the re-absorption by the surface plasmons leads to the dent structure in . Thus, in addition to the molecular luminescence and absorption, the re-absorption by the surface plasmons plays an important role in determining the luminescence-spectral profiles of the surface plasmons. Also, it is found that the effects of the interference between the processes of the molecular absorption and the re-absorption by the surface plasmons can be seen in . We expect that these results can be verified by experiments, for example by a comparison between the luminescence spectra acquired with a molecule-covers tip and with a clean metallic tip on a clean metal surface, which correspond to the red solid and green dashed lines in Fig. 5, respectively.

### 4.2 Molecular luminescence at the high bias voltage

To analyze the luminescence-spectral profile of the molecule, contributions from each term in Eq. (47) are shown in Fig. 7.

The contribution from the first term on the right side of Eq. (47) (green dashed line) leads to peak structures near 1.87 eV, 1.71 eV, 1.55 eV and so on. These structure correspond to the molecular luminescence associated with the transitions from the ground vibrational level for the excited electronic state to the vibrational levels for the ground electronic state of the molecule. The physical interpretation of the processes is that since the electronic and vibrational part in are decoupled, the creation (annihilation) of the molecular exciton induced by the absorption (emission) of the surface plasmons is not accompanied by the vibrational motions. Hence the first term expresses the molecular luminescence associated with the transitions from the ground vibrational level where the molecular vibration is initially distributed.

The third term (magenta dotted line) leads to peak structures near 2.05 eV, 2.21eV, and 2.37 eV. The position of these structures correspond to the energy of transitions from the excited vibrational levels for the excited electronic state to the vibrational levels for the ground electronic state of the molecule. The third term describes the processes, where the motion of the molecular exciton is accompanied by the vibrational motions. The molecule is excited to the vibrational levels for the excited electronic state by the absorption of the surface plasmons, and then de-excited to the ground vibrational level for the ground electronic state by subsequent emission of the surface plasmons. Thus, the third term leads to the luminescence, which are attributed to the radiative decay from the vibrational excited levels for the excited electronic state (hot luminescence). [31, 13]

The second and fourth terms give the processes, where the energy of the molecular exciton is absorbed by the surface plasmons and re-absorbed by the molecule. Effects of interference between the processes of the absorption by the surface plasmons and the re-absorption by the molecule can be seen, as described in the previous subsection.

### 4.3 Molecular luminescence at the low bias voltage

Figure 8 shows the luminescence spectra of the molecule at various bias voltages lower than .

Despite the fact that is lower than the energy of the molecular exciton , the molecular luminescence is found to occur. The results indicate that the electron transitions in the molecule occur at this bias voltage. To analyze the mechanism for the occurrence of the electron transitions at , the spectral function of the molecule is calculated and shown in Fig. 9.

Owing to the exciton-plasmon coupling , the position and the width of the peaks in are shifted and broadened, respectively. The spectral intensities are found in the energy range lower than . It means that the excitation channels of the molecule arise in this energy range. Thus, the electron transitions in the molecule occur via the excitation channels resulting from the exciton-plasmon coupling, and give rise to the molecular luminescence at . Corresponding luminescence spectra can be seen in a recent experiment. [53] The observed luminescence spectra shows peak structures near the energies of transitions in the molecule. The observed spectra also show the bias voltage dependence similar to Fig. 8. Thus, our results first reveal that the molecular luminescence at are attributed to the electron transitions in the molecule via the excitation channels, that result from the exciton-plasmon coupling.

Figure 10 shows the luminescence spectra of the molecule at .

A peak structure is found in the energy range . Luminescence with the energy higher than are called upconverted luminescence. One of the possible mechanisms for the occurrence of the upconverted luminescence is as follows: the electronic transitions of the molecule induced by the absorption and emission of the surface plasmons are accompanied by excitations of the molecular vibration, and the vibrational excitations assist the occurrence of the upconverted luminescence. The contribution of the vibrational excitations can be found in comparison with the vibrational state in thermal equilibrium, where the molecular vibration with the energy is distributed according to the Bose-Einstein distribution function at and is therefore almost in the ground state.

The dependence of luminescence spectra on and is also shown in Fig. 10. The luminescence intensity in the energy range is proportional to , and the intensity of the upconverted luminescence is proportional to the square of . As the energy of the surface plasmon mode is shifted to the low-energy side, the luminescence intensity increase. This is attributed to the fact that since the energy of the surface plasmon mode is lower than , the electron transitions in the molecule in the energy range are enhanced by the surface plasmons.

Figure 11 shows the dependence of the vibrational occupation number [red solid lines in Fig. 11] and the population of the molecular exciton [blue dotted lines in Fig. 11].

It is confirmed that the vibrational excitations occur at . Hence the vibrational excitations are found to assist the occurrence of the upconverted luminescence. The slope of changes at of approximately 1.85 eV for [Figs. 11(a) and 11(b)] and at of approximately 1.90 eV for [Figs. 11(c) and 11(d)]. At this bias voltage, the excitation channels of the molecule increase. The slope of the difference between and the vibrational occupational number for the vibrational state in thermal equilibrium changes at the same bias voltage. Therefore, the electron transitions in the molecule via the excitation channels resulting from the exciton-plasmon coupling give rise to the vibrational excitations at . It is found that is scaled by , where represents the vibrational damping for . Thus, even at , the vibrational excitations occur with the electron transitions in the molecule via the excitation channels resulting from the exciton-plasmon coupling and assist the occurrence of the upconverted luminescence.

## 5 Conclusion

The luminescence properties of the molecule and the surface plasmons are greatly influenced by the interplay between the dynamics of the molecule and the surface plasmons. The peak structures in arise due to the enhancement by the molecular modes and the dip structures due to the molecular absorption. The re-absorption by the surface plasmons leads to the dent structure in . The interference between the processes of the molecular absorption and the re-absorption by the surface plasmons is found to enhance or suppress the luminescence intensity of the surface plasmons, depending on the energy relative to and . The effects of the interference can be seen in . Moreover, due to the exciton-plasmon coupling, the excitation channels of the molecule arise even in the energy range lower than . The electron transitions via these excitation channels give rise to the molecular luminescence at . It is also found that even at , the vibrational excitation occur with the electron transitions via these excitation channels and assist the occurrence of the upconverted luminescence.

Our results will pave the way for understanding the dynamics in the system where multiple bosonic particles (e.g. molecular exciton, vibron, surface plasmon, etc.) are coupled to each other. The analysis of the molecular and plasmon luminescence may contribute to the development of plasmon enhanced organic light-emitting diode (OLED) devise and high-resolution spectroscopy and sensing. Especially, the upconverted luminescence may open up new opportunities to increase the luminescence intensity in the high-energy region and explore highly efficient blue light-emitting materials for OLED devices.

This work is supported in part by MEXT through the G-COE program, Grant-in-Aid for Scientific Research on Innovative Areas Program (2203-22104008) and Scientific Research (c) Program (22510107). It was also supported in part by JST through ALCA Program and Strategic Japanese-Croatian Cooperative Program on Materials Science. Some of the calculations presented here were performed using the ISSP Super Computer Center, University of Tokyo. The authors are deeply grateful to Professor Wilson Agerico Diño and Professor Hiroshi Nakanishi in Osaka University for useful discussions.

### References

- F. Tam, G. P. Goodrich, B. R. Johnson, and N. J. Halas: Nano Lett. 7 (2007) 496.
- M. Moskovits: Rev. Mod. Phys. 57 (1985) 783.
- S. Nie and S. R. Emory: Science 275 (1997) 1102.
- B. Pettinger, B. Ren, G. Picardi, R. Schuster, and G. Ertl: Phys. Rev. Lett. 92 (2004) 096101.
- D. Enders and A. Pucci: Appl. Phys. Lett. 88 (2006) 184104.
- J. Kundu, F. Le, P. Nordlander, and N. J. Halas: Chem. Phys. Lett. 452 (2008) 115.
- J. R. Lakowicz: Anal. Biochem. 298 (2001) 1.
- A. Fujiki, T. Uemura, N. Zettsu, M. Akai-Kasaya, A. Saito, and Y. Kuwahara: Appl. Phys. Lett. 96 (2010) 043307.
- H. A. Atwater and A. Polman: Nat. Materials 9 (2010) 205.
- S. Pillai and M. A. Green: Sol. Energy Mater. Sol. Cells 94 (2010) 1481.
- H. Ishihara, A. Nobuhiro, M. Nakatani, and Y. Mizumoto: J. Photochem. Photobio. A: Chemistry 221 (2011) 148.
- A. J. White, B. D. Fainberg, and M. Galperin: J. Phys. Chem. Lett. 3 (2012) 2738.
- N. L. Schneider and R. Berndt: Phys. Rev. B 86 (2012) 035445.
- E. J. Osley, C. G. Biris, P. G. Thompson, R. R. F. Jahromi, P. A. Warburton, and N. C. Panoiu: Phys. Rev. Lett. 110 (2013) 087402.
- J. K. Gimzewski, B. Reihl, J. H. Coombs, and R. R. Schlittler: Z. Phys. B 72 (1988) 497.
- R. Berndt, J. K. Gimzewski, and P. Johansson: Phys. Rev. Lett. 67 (1991) 3796.
- Y. Uehara, T. Fujita, and S. Ushioda: Phys. Rev. Lett. 83 (1999) 2445.
- G. Hoffmann, J. Kliewer, and R. Berndt: Phys. Rev. Lett. 87 (2001) 176803.
- M. Sakurai, C. Thirstrup, and M. Aono: Applied Physics A 80 (2005) 1153.
- F. Rossel, M. Pivetta, and W. D. Schneider: Surf. Sci. Rep. 65 (2010) 129.
- G. Hoffmann, L. Libioulle, and R. Berndt: Phys. Rev. B 65 (2002) 212107.
- K. Amemiya: Phys. Rev. B 67 (2003) 075409.
- X. Tao, Z. C. Dong, J. L. Yang, Y. Luo, J. G. Hou, and J. Aizpurua: J. Chem. Phys. 130 (2009) 084706.
- X. H. Qiu, G. V. Nazin, and W. Ho: Science 299 (2003) 542.
- Z. C. Dong, X. L. Guo, A. S. Trifonov, P. S. Dorozhkin, K. Miki, K. Kimura, S. Yokoyama, and S. Mashiko: Phys. Rev. Lett. 92 (2004) 086801.
- D. Ino, T. Yamada, and M. Kawai: J. Chem. Phys. 129 (2008) 014701.
- C. Chen, P. Chu, C. A. Bobisch, D. L. Mills, and W. Ho: Phys. Rev. Lett. 105 (2010) 217402.
- T. Uemura, M. Furumoto, T. Nakano, M. Akai-Kasaya, A. Saito, M. Aono, and Y. Kuwahara: Chem. Phys. Lett. 448 (2007) 232.
- H. W. Liu, Y. Ie, R. Nishitani, Y. Aso, and H. Iwasaki: Phys. Rev. B 75 (2007) 115429.
- H. W. Liu, R. Nishitani, T. Z. Han, Y. Ie, Y. Aso, and H. Iwasaki: Phys. Rev. B 79 (2009) 125415.
- Z. C. Dong, X. L. Zhang, H. Y. Gao, Y. Luo, C. Zhang, L. G. Chen, R. Zhang, X. Tao, Y. Zhang, J. L. Yang, and J. G. Hou: Nat. Photonics 4 (2010) 50.
- C. Zhang, R. Zhang, S. Jiang, L. Zhang, H. Y. Gao, X. L. Zhang, L. G. Chen, Y. Liao, and Z. C. Dong: Appl. Phys. Lett. 100 (2012) 073111.
- A. Ishida and T. Majima: Nanotechnology 10 (1999) 308.
- H. Xu, X.-H. Wang, M. P. Persson, H. Q. Xu, M. Käll, and P. Johansson: Phys. Rev. Lett. 93 (2004) 243002.
- P. Johansson, H. Xu, and M. Käll: Phys. Rev. B 72 (2005) 035427.
- G. Tian, J. C. Liu, and Y. Luo: Phys. Rev. Lett. 106 (2011) 177401.
- G. Tian and Y. Luo: Phys. Rev. B 84 (2011) 205419.
- K. Miwa, M. Sakaue, and H. Kasai: to be published in J. Phys. Soc. Jpn.
- K. Miwa, M. Sakaue, and H. Kasai: Nanoscale Research Letters 8 (2013) 204.
- It should be noted that a more general model for the molecular vibrations, e.g., the one taking into account multiple vibrational modes and/or an intramolecular anharmonicity, can be found in the literature, and/or the intramolecular anharmonicity will be investigated in future works.
- H. Kasai and A. Okiji: Progress in Surface Science 44 (1993) 101.
- H. Kasai, W. A. Diño, and A. Okiji: Surface Science Reports 43 (2001) 1.
- B. Pettinger, K. F. Domke, D. Zhang, R. Schuster, and G. Ertl: Phys. Rev. B 76 (2007) 113409.
- M. Galperin, A. Nitzan, and M. A. Ratner: Phys. Rev. B 73 (2006) 045314.
- M. Galperin, M. A. Ratner, and A. Nitzan: J. Chem. Phys. 130 (2009) 144109.
- For example, when the molecule is directly adsorbed on the tip surface, the molecular exciton has a finite lifetime due to the energy transfer to the metal tip. This effect can be taken into account by introducing a relaxation constant . The lifetime of the molecular exciton for is proportional to . When is large, the luminescence intensities of the molecule are expected to be much lower than the surface plasmons. In this case, the photon emission rate is approximately proportional to . Although has not been taken into account here, we have confirmed that leads only to variation in the luminescence intensities and the peak width in . Thus, the luminescence spectra detected experimentally are expected to correspond to .
- L. V. Keldysh: Sov. Phys. JETP 20 (1965) 1018.
- B. Gumhalter: Prog. Surf. Sci. 15 (1984) 1.
- S. Maier, T. L. Schmidt, and A. Komnik: Phys. Rev. B 83 (2011) 085401.
- B. N. J. Persson and A. Baratoff: Phys. Rev. Lett. 68 (1992) 3224.
- H. Haug and A.-P. Jauho: Quantum Kinetics in Transport and Optics of Semiconductors (Springer, Berlin, 2007) 2nd edition ed., p. 362.
- G. D. Mahan: Many-Particle Physics (Plenum Pub Corp, New York, 2000) 3rd edition ed., p. 785.
- A. Fujiki, Y. Miyake, Y. Oshikane, M. Akai-Kasaya, A. Saito, and Y. Kuwahara: Nanoscale Research Letters 6 (2011) 347.