Cooperative Spontaneous Emission from Nanocrystals to a Surface Plasmon Polariton in a Metallic Nanowire

  • cc icon
  • ABSTRACT

    We analyze the cooperative spontaneous emission of optically excited nanocrystals into surface plasmon polaritons propagating on the surface of a cylindrical metallic nanowire. The spontaneous emission probability of the nanocrystals is obtained by perturbative expansions with and without dipole-dipole interaction among nanocrystals in order to see the cooperative effects. The spontaneous emission probability depends on the radial and axial distributions, as well as on the dipolar orientation of nanocrystals. It is shown that the spontaneous emission probability is strongly influenced by dipole-dipole interaction, axial distribution, and dipolar orientation of nanocrystals for closely spaced nanocrystals.


  • KEYWORD

    Plasmonic nanowire , Nanocrystals , Cooperative spontaneous emission

  • I. INTRODUCTION

    The optical properties of composite nanostructures including metal are strongly affected by surface plasmon polaritons’ (SPPs) collective oscillations of conduction electrons with respect to the positive ion background [1-5]. SPPs are the propagating excitations of charge-density waves, and their associated intense electromagnetic fields are generated in the vicinity of the surface of a conductor. They are bound to and propagate along the metal/dielectric interface. Recently, great attention has been paid to the study of the strong coupling between individual emitters and SPPs in metallic nanowires [6-9]. The emission properties of nanoscale optical emitters can be significantly modified by the proximity of a nanowire that supports surface plasmons because an SPP is generally very intense and highly localized [5,7,10]. This subwavelength localization is accompanied by a dramatic condensation of optical fields. A substantial increase in the coupling strength g1/√Veff can be achieved near the surface of metal nanostructures, since the effective mode volume Veff of the SPPs can be considerably reduced. For an optical emitter placed within the evanescent SPP mode tail, the spontaneous emission probability of the excited energy of an emitter into the surface plasmons is tremendously enhanced [11-14]. Thus, the controllable strong near-field coupling between emitters and SPPs in metal nanostructures has stimulated a number of studies on nano-plasmonics [15-18].

    Near-field coupling of excited colloidal nanocrystals (NCs) to SPPs in a metallic nanowire (NW) is a challenging problem since nanocrystals are randomly dispersed relative to the metallic nanowire [17]. In this paper, we investigate the spontaneous emission (SE) of the excited energy of one NC in NC aggregation to the SPPs in a nearby metallic nanowire for some specified configurations and for specific transition dipolar orientation of NCs. The radius of the NW is to be so small that it is sufficient to consider only the fundamental TE SPP mode in our calculations. In this work, we analyze the cooperative effect in the spontaneous decay of an excited NC in an NC aggregate into a SPP mode with and without the dipole-dipole interaction between NCs [8]. The spontaneous emission rate into surface plasmons can be suppressed or enhanced for specially arranged NCs, depending on their relative positions and transition dipolar orientations due to their quantum interference between the interactions. This variable emission rate can be characterized by cooperativeness, which for thin wires and closely spaced NCs is predicted to be large. In the next section, a theoretical modeling is employed to calculate the transition amplitude of the spontaneous emission under the perturbative expansion scheme. Results of the numerical evaluation for the transition probability (TP) of the spontaneous emission are presented in Sec. III followed by a discussion on its characteristics. A brief conclusion is given in Sec. IV.

    II. THEORETICAL MODEL

    The system of interest in this paper is schematically shown in Fig. 1. NCs are supposed to be aggregated around the silver nanowire along the axis, where only one of the NCs is excited. First, we briefly review the properties of SPPs propagating along the surface of metallic NW [19,20], which have been observed experimentally [21-24]. We are especially interested in azimuthally independent SPPs, which are supported under the condition that the radius of the wire is small enough such that only the fundamental TM surface plasmon polariton mode exists while all other higher order modes are cut off [7]. In this case the electric and magnetic fields of the fundamental TM surface plasmon polariton mode are given by [19]

    image
    image

    where q stands for inside or outside of the metallic NW with

    image
    image
    image

    Here Jn and Hn(1) are Bessel and Hankel function, respectively. The dielectric function ε q(ω) is taken from the Drude- Lorentz model, i.e., ?in(ω) = ?[1-ωp2/ ω (ω+i γ)], where γ is the relaxation rate due to ohmic loss. Since the system does not have azimuthal dependence, the wave number of the SPP has only z dependence. From now on we drop the subscript z for the wave number kz. The dispersion relation (k, ωk) of the SPP mode is obtained by applying the boundary condition of the TM mode fields on the surface of the metallic NW [8, 11, 19, 20]

    image

    where a is the radius of the NW. The dispersion relation for the fundamental mode of SPP on metallic NW is shown in Fig. 2. The mode-expansion coefficients aq(k) in Eqs. (1) and (2) are determined by matching the boundary condition of the fields on the metal nanowire/dielectric interface and also by normalizing the electromagnetic fields to their vacuum fluctuation energies [1,8,11]. After performing these operations, the resulting electric field outside the metal becomes

    image

    where b(k) is again the mode expansion coefficient but normalized by electric permittivity of the vacuum and nanowire length L, which is averaged over the azimuthal angle for the wire. The radial component of the field is (almost) real and the axial component of the field is imaginary, so that they are out of phase by π /2. The absolute value of the

    electric field outside the nanowire as a function of wave number is shown in Fig. 3(a); solid and dashed lines are for the radial and axial components of the field, respectively. Fig. 3(b) is the electric field as a function of radial distance from the wire center. Dramatic concentration of optical fields is well manifested.

    The Hamiltonian of this composite system is written by H = H0 Hint, where H0 is the unperturbed part and the interaction Hint is given by

    image

    where

    image

    is the transition dipole moment of the jth NC and

    image

    is the electric field at the location of the jth NC due to the surface plasmon polaritons in the nanowire. rmj and

    image

    are the distance and the directional unit vector between mth and jth NCs, respectively. The first term represents the interaction between NCs and NW plasmons, while the second term describes the dipole-dipole interaction (DDI) between NCs. The transition energy of an NC is assumed to be detuned from the surface plasmon polarion mode energy of the metallic nanowire.

    The transition amplitude of the SE is calculated from perturbation theory. Within the dipole approximation, the first order transition amplitude of the spontaneous emission from an excited j th NC at (ρj, zj) with transition energy ω j to a fundamental SPP of mode with wave number k in metallic NW is obtained as

    image

    where

    image

    and Δ kj = ω k?ωj with

    image

    being the electric field of the surface plasmon mode of wave number k at the location of the jth NC. Δkj describes the detuning between jth NC’s transition energy and SPP mode energy with wave number k. If all NCs are identical with ωm = ωj and the direct DDI between NCs is ignorable, the third order perturbative transition amplitude of the SE is calculated as

    image

    The even-order perturbation terms are null, since the coupling is linear to the dipole moment. Since the NC’s transition energy is detuned from the SPPs in the metallic nanowire and sin(Δ kk' t/2)/Δ kk' is a very strongly peaked function of Δkk', the last term dominates. Thus Eq. (10) can be approximated by

    image

    When all NCs are not identical, the third order perturbative transition amplitude of the SE is similarly calculated, but in this case the detuning parameter Δmj between NCs is, of course, explicitly included in the expression

    image

    On the other hand, if the direct DDI between NCs cannot be ignored, the second order transition amplitude of the SE term survives. From the calculation carried out for the case when all NCs are identical with only the nearest neighbor interaction, the transition amplitude of the SE is

    image

    where Jmj is the dipole-dipole interaction energy between mth and jth NCs in Eq. (8). Similarly, when all NCs are not identical, the second order transition amplitude of the SE is calculated as

    image

    III. NUMERICAL EVALUATION AND DISCUSSION

    The first order TP,

    image

    for the SE of an excited NC to a SPP of mode k does not exhibit any collective effect, since it only involves single NC coupling to SPP mode regardless of the existence of other NCs. In order to see the collective effect, the higher order perturbative terms must be examined. We first consider the case without DDI. When we evaluate the total TP composed of

    image

    and

    image

    the TP would have the cross product contributions between the first and third order transition amplitudes of the SEs. However, neither do they show any substantial collective effect on the same footing as the first order contribution alone. Then, the most important contribution with respect to the cooperative effect comes from the square of

    image

    From Eq. (10), the third order TP for the SE of an NC to an SPP of mode

    image

    is obtained by

    image

    where

    image

    and

    image

    Here, f (k, k', k", t) is the absolute square of the temporal part in Eq. (10) and it is a function of Δk0, Δ k'0, Δ k"0, Δ kk', and Δ kk", so that it represents the consequence of detuning effect. As seen from Eq. (15), characteristic collective behavior of this composed system in the third order TP is actually described by the phase factor exp(iΦ n, m, k', k").

    The cooperativeness in the third order TP is manifested by the phase Φn, m, k', k", which describes the correlation among the whole NC aggregate. Since the mode of SPPs does not depend on the azimuthal angle, the phase Φ n, m, k', k" of Eq. (16) depends only on the radial and axial locations of NCs. Thus, the cooperativeness is determined by the radial and axial arrangement of NCs, i.e., the locations of NCs ( ρj, zj) and by the dipolar orientation. In order to analyze the collective effect, we consider the dependence of the phase factor on the radial and axial arrangement as well as on the dipolar orientation of NCs. First, the phase difference caused by the axial distribution of NCs in the form of k(zn?  zm) in Eq. (16) is always non-negligible since kΔz ~ π even for Δz ~ O(1nm). Thus the axial distribution of NCs along the wire can be an important factor to induce the collective effect in the SE. Next, for the analysis of the phase caused by the radial arrangement and the dipolar orientation, which is represented by the term of θj,k in Eq. (16), we consider, for simplicity, the case when the transition dipole vector is in the plane of (ρ, z),. In this case, the absolute magnitude of the NC-SPP coupling strength ?gj,k? and its phase θj,k are expressed as

    image

    and

    image

    respectively, where αj is the orientation of the transition dipole measured from the wire axis. |gjk| and θjk as functions of k and αj are shown in Fig. 4(a) and 4(b), respectively.

    Now we consider the collective effect of the SE in the third order TP. As was pointed out in the previous section, the terms with (sin Δ kk' t/2)/Δ kk' and (sin Δ ;kk" t/2)/Δ kk" in f(k,k',k",t) of Eq. (15) dominate, since NC’s transition energy is detuned from the SPPs in the metallic NW. The amplitude of the function (sin Δkk' t/2)/Δkk', which is very highly peaked at Δ kk' = 0, increases as t and decreases to zero when Δkk' = 2π/t. When time t ~ zn (dk/dω) is small enough as ~ fs or less, non-resonant coupling between SPP modes is allowed in the third order process. In this situation, the phase Φof Eq. (16) is not negligible. That is, the phase difference with the form of zn(k?k') in Eq. (16) is order of π for the time being ~ fs or less in which Δ z is assumed to be ±10nm. On the other hand, in the longtime limit, Δkk' ? Δkk" ?0 and for the most part only the resonant coupling between SPP modes is allowed in the third order process. In this circumstance, there are many cancellations among the phases in Eq. (16), so that no appreciable collective effect is expected. In this resonant coupling approximation, the third order transition probability of the SE can be approximated as

    image

    which is plotted as a function of k in Fig. 5(a). Since the transition energy of an NC is assumed to be detuned from the SPP mode energy of the NW and to be above the mode energy, the TP of the SE will increase as the size of NCs increases. That is, when the size of NCs increases, the transition energy of NCs decreases and then Δkj in Eq. (19) will decrease. In the case of non-identical NCs, it can be further approximated as Pk(3)(t)(1 Δ jm / Δ kj)-2, which shows reduced transition probability of the SE for increased disparity in the size of NCs. As a reference, the first order TP of the SE to the k mode of SPP Pk(1)(t) is given by

    image

    and is shown in Fig. 5(b), for comparison. No noticeable

    change for TP is observed with regard to including multiple interactions between NC and SPP.

    Now we analyze the cooperative effect due to direct DDI. When all NCs are identical, the second order TP of the SE to the k mode of SPP, Pk(2)(t) is obtained as

    image

    where

    image

    The cooperativeness of NCs is involved in the phase Φas the case without DDI. As described above, the phase difference caused by the axial distance between NCs, is order of π for Δz ~ 10 nm. On the other hand, the phase difference, Δθ = θ m,kθn,k is induced by the radial distribution and the dipolar orientation of NCs, as seen from Eq. (18). Since the radial component of the electric field is (almost) real and the axial component of the electric field is imaginary, there is no interference effect coming from the NC-SSP coupling difference as seen from Eq. (18), if the transition dipole is arranged only along the wire or only perpendicular to the wire axis. Thus the phase difference Δ θ in Eq. (22) depends strongly on the dipolar orientation. This character is plotted as a function of α and k in Fig. 6(a). The figure shows that, for a given transition dipolar orientation of NCs, the SE of an excited NC into several specified modes of the fundamental SPP in the metallic nanowire is totally inhibited.

    When the radial distances of NCs from the wire center are almost the same among the NC aggregates (?gm,k???gk?) and the orientations of NC’s transition dipoles are parallel to each other with DDI coupling J, the transition probability of the SE is approximated as

    image

    In this case, the phase difference Δθ of Eq. (22) almost vanishes and the phase is determined by the difference of axial distance between NCs, which is plotted as a function of Δz and k in Fig. 6(c). The effect of interference on the SE due to the second order TP caused by the direct DDI is explicitly shown in the figure. For a given configuration of NCs, the SE to a certain specified SPP modes is totally

    suppressed, as expected. And when all NCs are arranged along the NW, the TP of the SE will decrease as the distance between the NCs and the NW increases as seen from Fig. 3(b).

    In practice, the phase difference Δθ caused by the radial distribution between NCs is appreciable, however the effect of the radial distribution of NC on the cooperativeness in the second order TP is negligible due to the huge difference of the field |E(ρn)|-|E(ρm)|, which causes low visibility for the interference pattern of Eq. (21). Since the phase difference Δθ depends strongly on the dipolar orientation, the effect of dipolar orientation on the cooperativeness in the second order TP is considered by dealing with two simple cases. For instance, when there are three NCs with similar radial distances ρjρ and if

    image
    image

    the second order TP as a function of α and k is obtained as shown in Fig. 7(a). On the other hand, if the first and the third dipoles point along the radial direction, i.e.,

    image

    the transition probability of the SE as a function of α and k behaves like Fig. 7(b). Here, we have analyzed the TP of the SE when only one NC is

    excited. Meanwhile, if two or more NCs are excited, then the TP of the SE will strongly depend on the initially prepared state. For the case without the interaction between excitons, if the NCs are excited in the form of a superradiative (subradiative) state, the emission rate will be enhanced (suppressed) [25].

    IV. CONCLUSION

    In summary, we have analyzed the characteristics of the cooperative feature of the spontaneous emission of one excited NC among NC aggregate into a mode of SPP in a metallic nanowire. In order to see the collective effect of the SE, higher order perturbation expansions of the transition probability of the SE have been derived. In particular, the influence of the dipolar orientation and the axial and radial arrangement of NCs has been considered. The first-order transition probability of the SE does not show collective behavior. Also the collective effect disappears in the thirdorder process of the spontaneous emission when direct DDI is not involved under the resonant coupling approximation between SPPs. However, the second-order transition probability of the SE yields a cooperative effect of the SE. In particular, the spontaneous emission of the excited NC energy is suppressed or enhanced depending on the relative orientation of the dipoles and the arrangement of NCs.

  • 1. Maier S. A. 2007 Plasmonics: Fundamentals and Applications google
  • 2. Wang H., Brandl D. W., Nordlander P., Halas N. J. (2007) “Plasmonic nanostructures: artificial molecules” [Acc. Chem. Res.] Vol.40 P.53-62 google doi
  • 3. Fernandez-Dominguez A. I., Martin-Moreno L., Garcia- Vidal F. J., Andrews S. R., Maier S. A. (2008) “Spoof surface plasmon polariton modes propagating along periodically corrugated wires” [IEEE J. Quantum Electron] Vol.14 P.1515-1521 google doi
  • 4. Yang J., Cao Q., Zhou C. (2010) “Analytical recurrence formula for the zeroth-order metal wire plasmon of terahertz waves” [J. Opt. Soc. Am. A] Vol.27 P.1608-1612 google doi
  • 5. Wei H., Ratchford D., Li X., Xu H., Shih C.-K. (2009) “Propagating surface plasmon induced photon emission from quantum dots” [Nano Lett.] Vol.9 P.4168-4171 google doi
  • 6. Ditlbacher H., Hohenau A., Wagner D., Kreibig U., Rogers M., Hofer F., Aussenegg F. R., Krenn J. R. (2005) “Silver nanowires as surface plasmon resonators” [Phys. Rev. Lett.] Vol.95 P.257403 google doi
  • 7. Chang D. E., Sørensen A. S., Hemmer P. R., Lukin M. D. (2007) “Strong coupling of single emitters to surface plasmons” [Phys. Rev. B] Vol.76 P.035420 google doi
  • 8. Cheng M.-T., Luo Y.-Q., Song Y.-Y., Zhao G.-X. (2010) “Plasmonic waveguides mediated energy transfer between two distant quantum dots” [J. Mod. Opt.] Vol.57 P.2177-2181 google doi
  • 9. Kien F. L., Gupta S. D., Nayak K. P., Hakuta K. (2005) “Nanofiber-mediated radiative transfer between two distant atoms” [Phys. Rev. B] Vol.72 P.063815 google
  • 10. Chang D. E., Sørensen A. S., Demler E. A., Lukin M. D. (2007) “A single-photon transistor using nanoscale surface plasmons” [Nat. Phys.] Vol.3 P.807-812 google doi
  • 11. Chen Y. N., Chen G. Y., Chuu D. S., Brandes T. (2009) “Quantum-dot exciton dynamics with a surface plasmon: band-edge quantum optics” [Phys. Rev. A] Vol.79 P.033815 google doi
  • 12. Akimov A. V., Mukherjee, A., Yu C. L., Chang D. E., Zibrov A. S., Hemmer P. R., Park H., Lukin M. D. (2007) “Generation of single optical plasmons in metallic nanowires coupled to quantum dots” [Nature] Vol.450 P.402-406 google doi
  • 13. Chen G.-Y., Chen Y.-N., Chuu D.-S. (2008) “Spontaneous emission of quantum dot excitons into surface plasmons in a nanowire” [Opt. Lett.] Vol.33 P.2212-2214 google doi
  • 14. Rukhlenko I. D., Handapangoda D., Premaratne M., Fedorov A. V., Baranov A. V., Jagadish C. (2009) “Spontaneous emission of guided polaritons by quantum dot coupled to metallic nanowire: beyond the dipole approximation” [Opt. Express] Vol.17 P.17570-17581 google doi
  • 15. Witthaut D., Sorensen A. S. (2010) “Photon scattering by a three-level emitter in a one-dimensional waveguide” [New J. Phys.] Vol.12 P.043052 google doi
  • 16. Lei, D. Y., Aubry A., Maier S. A., Pendry J. B. (2010) “Broadband nano-focusing of light using kissing nanowires” [New J. Phys.] Vol.12 P.093030 google doi
  • 17. O'Carroll D. M., Hofmann C. E., Atwater H. A. (2010) “Conjugated polymer/metal nanowire heterostructure plasmonic antennas” [Adv. Mater.] Vol.22 P.1223-1227 google doi
  • 18. John S., Quang T. (1994) “Spontaneous emission near the edge of a photonic band gap” [Phys. Rev. A] Vol.50 P.1764-1769 google doi
  • 19. Stratton J. A. 1941 Electromagnetic Theory google
  • 20. Pfeiffer C. A., Economou E. N., Ngai K. L. (1974) “Surface polaritons in a circularly cylindrical interface: surface plasmons” [Phys. Rev. B] Vol.10 P.3038-3051 google doi
  • 21. Schider G., Krenn J. R., Hohenau A., Ditlbacher H., Leitner A., Aussenegg F. R., Schaich W. L., Puscasu I., Monacelli B., Boreman G. (2003) “Plasmon dispersion relation of Au and Ag nanowires” [Phys. Rev. B] Vol.68 P.155427 google doi
  • 22. Vesseur E. J. R., de Waele R., Kuttge M., Polman A. (2007) “Direct observation of plasmonic modes in Au nanowires using high-resolution cathodoluminescence spectroscopy” [Nano Lett.] Vol.9 P.2843-2846 google
  • 23. Schmidt M. A., Russell P. St. J. (2008) “Long-range spiralling surface plasmon modes on metallic nanowires” [Opt. Express] Vol.16 P.13617-13623 google doi
  • 24. Sanders A. W., Routenberg D. A., Wiley B. J., Xia Y., Dufresne E. R., Reed M. A. (2006) “Observation of plasmon propagation, redirection, and fan-out in silver nanowires” [Nano Lett.] Vol.6 P.1822-1826 google doi
  • 25. Hong S.-K., Nam S. W., Yeon K.-H. (2007) “Coherent exciton states of excitonic nanocrystal-molecules” [Phys. Rev. B] Vol.76 P.115330 google doi
  • [FIG. 1.] Schematic view of the system. NCs are located in the vicinity of silver nanowire. Cylindrical coordinates are used with z axis as the wire axis.
    Schematic view of the system. NCs are located in the vicinity of silver nanowire. Cylindrical coordinates are used with z axis as the wire axis.
  • [FIG. 2.] Dispersion relation of the fundamental SPP mode on lossless silver, where ? ∞ = 9.6 for silver, ? ∞ = 5.3 for dielectric (GaN), ω p = 3.76 eV/h, γ = 0, and the radius of the NW is set to be a = 5 nm
    Dispersion relation of the fundamental SPP mode on lossless silver, where ?
∞
 = 9.6 for silver, ?
∞
= 5.3 for dielectric (GaN), ω
p = 3.76 eV/h, γ
 = 0, and the radius of the NW is set to be a = 5 nm
  • [FIG. 3.] Electric field outside the wire. (a) as a function of wave number for ρ = 1.2a and (b) as a function of radial distance from the wire center for k = 107/cm. Arbitrary unit is used for vertical axis. The solid and dashed lines are for the radial and axial components of the field, respectively. The radial component is real and the axial component is imaginary.
    Electric field outside the wire. (a) as a function of wave number for ρ = 1.2a and (b) as a function of radial distance from the wire center for k = 107/cm. Arbitrary unit is used for vertical axis. The solid and dashed lines are for the radial and axial components of the field, respectively. The radial component is real and the axial component is imaginary.
  • [FIG. 4.] NC-SPP coupling constant gj,k. (a) The coupling strength ?gj,k? as a function of the dipolar orientation and the wave number for ρ = 1.2a. (b) The phase of the coupling constant θj,k as a function of the dipolar orientation αand the wave number k.
    NC-SPP coupling constant gj,k. (a) The coupling strength ?gj,k? as a function of the dipolar orientation and the wave number for ρ = 1.2a. (b) The phase of the coupling constant θj,k as a function of the dipolar orientation αand the wave number k.
  • [FIG. 5.] Transition probability of the spontaneous emission to the k mode of SPP without temporal part as a function of time and wave number for ρ =1.2a. (a) The third order transition probability of the SE under resonant approximation. (b) The first order transition probability of the SE.
    Transition probability of the spontaneous emission to the k mode of SPP without temporal part as a function of time and wave number for ρ
 =1.2a. (a) The third order transition probability of the SE under resonant approximation. (b) The first order transition probability of the SE.
  • [FIG. 6.] The second order transition probability of the SE to the k mode of SPP when all DDI coupling is the same and dipolar directions are all the same, where Δz = zj 1- zj = α , and f(t) is set to be unity. (a) TP as a function of k and α . (b) TP as a function of for α = 0. Solid and dashed lines are for μm = ρ and μ m = z , respectively. (c) TP as a function of k and Δ z for α=π /4. The range of the inter-NC distance Δz is from a to 6α .
    The second order transition probability of the SE to the k mode of SPP when all DDI coupling is the same and dipolar directions are all the same, where Δz
= zj 1- zj = α
, and f(t) is set to be unity. (a) TP as a function of k and α
. (b) TP as a function of for α = 0. Solid and dashed lines are for μm = ρ and μ
m = z , respectively. (c) TP as a function of k and Δ
z for α=π
/4. The range of the inter-NC distance Δz is from a to 6α
.
  • [FIG. 7.] The second order transition probability of the SE to the k mode of SPP when ρj ?ρ with r mj= z, μ 1 = ρ, μ2 ?z = cosθand Δz = a. (a) TP as a function of k and α when μ 3=z . (b) TP as a function of k and αwhen μ3 = ρ.
    The second order transition probability of the SE to the k mode of SPP when ρj ?ρ with r mj= z, μ
1 = ρ, μ2 ?z = cosθand Δz = a. (a) TP as a function of k and α when μ
3=z . (b) TP as a function of k and αwhen μ3 = ρ.