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 dipoledipole 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 dipoledipole interaction, axial distribution, and dipolar orientation of nanocrystals for closely spaced nanocrystals.
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 chargedensity 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 g
^{∝}
1/√
V_{eff}
can be achieved near the surface of metal nanostructures, since the effective mode volume
V_{eff}
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 nearfield coupling between emitters and SPPs in metal nanostructures has stimulated a number of studies on nanoplasmonics
[15

18]
.
Nearfield 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 dipoledipole 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]
where
q
stands for inside or outside of the metallic NW with
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.
Here
J_{n}
and
H_{n}
^{(1)}
are Bessel and Hankel function, respectively. The dielectric function
ε _{q}
(
ω
) is taken from the Drude Lorentz model, i.e., ∊
_{in}
(
ω
) = ∊
_{∞}
[1
ω_{p}^{2}
/
ω
(
ω
+
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
k_{z}
. 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]
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 modeexpansion coefficients
a_{q}
(
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
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
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
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 = 10^{7}/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 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
=
H
_{0}
H
_{int}
, where
H
_{0}
is the unperturbed part and the interaction Hint is given by
where
is the transition dipole moment of the
j
th NC and
is the electric field at the location of the
j
th NC due to the surface plasmon polaritons in the nanowire.
r_{mj}
and
are the distance and the directional unit vector between
m
th and
j
th NCs, respectively. The first term represents the interaction between NCs and NW plasmons, while the second term describes the dipoledipole 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}
,
z_{j}
) with transition energy
ω _{j}
to a fundamental SPP of mode with wave number
k
in metallic NW is obtained as
where
and Δ
_{kj}
= ω
_{k}
–ω
_{j}
with
being the electric field of the surface plasmon mode of wave number
k
at the location of the
j
th NC. Δ
_{kj}
describes the detuning between
j
th 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
The evenorder 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
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
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
where
J_{mj}
is the dipoledipole interaction energy between
m
th and
j
th NCs in Eq. (8). Similarly, when all NCs are not identical, the second order transition amplitude of the SE is calculated as
III. NUMERICAL EVALUATION AND DISCUSSION
The first order TP,
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
and
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
From Eq. (10), the third order TP for the SE of an NC to an SPP of mode
is obtained by
where
and
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}
,
z_{j}
) 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(
z_{n}
–
z_{m}
) in Eq. (16) is always nonnegligible since
k
Δz
~
π
even for
Δz
~ O(1
nm
). 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 NCSPP coupling strength ?
g_{j,k}
? and its phase θj,k are expressed as
and
respectively, where
α_{j}
is the orientation of the transition dipole measured from the wire axis. 
g_{jk}
 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
~
z_{n}
(
dk
/
d
ω) is small enough as ~
fs
or less, nonresonant 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
z_{n}
(
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
NCSPP coupling constant g_{j,k.} (a) The coupling strength ?g_{j,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.
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 nonidentical NCs, it can be further approximated as
P_{k}
^{(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
P_{k}
^{(1)}
(
t
) is given by
and is shown in
Fig. 5
(b), for comparison. No noticeable
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.
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,
P_{k}
^{(2)}
(
t
) is obtained as
where
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 NCSSP 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 (∣
g_{m,k}
∣≈∣
g_{k}
∣) 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
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
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= z_{j 1} z_{j} = α, 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 interNC distance Δz is from a to 6α.
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
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.,
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
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} = ρ.
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 firstorder 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 secondorder 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.
Acknowledgements
This work has been supported by a Korea University Grant.
Maier S. A.
2007
Plasmonics: Fundamentals and Applications
Springer
New York, USA
Chapter 1.
Wang H.
,
Brandl D. W.
,
Nordlander P.
,
Halas N. J.
(2007)
“Plasmonic nanostructures: artificial molecules”
Acc. Chem. Res.
40
53 
62
DOI : 10.1021/ar0401045
FernandezDominguez A. I.
,
MartinMoreno 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
14
1515 
1521
DOI : 10.1109/JSTQE.2008.918107
Yang J.
,
Cao Q.
,
Zhou C.
(2010)
“Analytical recurrence formula for the zerothorder metal wire plasmon of terahertz waves”
J. Opt. Soc. Am. A
27
1608 
1612
DOI : 10.1364/JOSAA.27.001608
Wei H.
,
Ratchford D.
,
Li X.
,
Xu H.
,
Shih C.K.
(2009)
“Propagating surface plasmon induced photon emission from quantum dots”
Nano Lett.
9
4168 
4171
DOI : 10.1021/nl9023897
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.
95
257403 
DOI : 10.1103/PhysRevLett.95.257403
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
76
035420 
DOI : 10.1103/PhysRevB.76.035420
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.
57
2177 
2181
DOI : 10.1080/09500340.2010.532573
Kien F. L.
,
Gupta S. D.
,
Nayak K. P.
,
Hakuta K.
(2005)
“Nanofibermediated radiative transfer between two distant atoms”
Phys. Rev. B
72
063815 
Chang D. E.
,
Sørensen A. S.
,
Demler E. A.
,
Lukin M. D.
(2007)
“A singlephoton transistor using nanoscale surface plasmons”
Nat. Phys.
3
807 
812
DOI : 10.1038/nphys708
Chen Y. N.
,
Chen G. Y.
,
Chuu D. S.
,
Brandes T.
(2009)
“Quantumdot exciton dynamics with a surface plasmon: bandedge quantum optics”
Phys. Rev. A
79
033815 
DOI : 10.1103/PhysRevA.79.033815
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
450
402 
406
DOI : 10.1038/nature06230
Chen G.Y.
,
Chen Y.N.
,
Chuu D.S.
(2008)
“Spontaneous emission of quantum dot excitons into surface plasmons in a nanowire”
Opt. Lett.
33
2212 
2214
DOI : 10.1364/OL.33.002212
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
17
17570 
17581
DOI : 10.1364/OE.17.017570
Lei, D. Y.
,
Aubry A.
,
Maier S. A.
,
Pendry J. B.
(2010)
“Broadband nanofocusing of light using kissing nanowires”
New J. Phys.
12
093030 
DOI : 10.1088/13672630/12/9/093030
O'Carroll D. M.
,
Hofmann C. E.
,
Atwater H. A.
(2010)
“Conjugated polymer/metal nanowire heterostructure plasmonic antennas”
Adv. Mater.
22
1223 
1227
DOI : 10.1002/adma.200902024
Stratton J. A.
1941
Electromagnetic Theory
McGrawHill
New York, USA
Chapter 9.
Pfeiffer C. A.
,
Economou E. N.
,
Ngai K. L.
(1974)
“Surface polaritons in a circularly cylindrical interface: surface plasmons”
Phys. Rev. B
10
3038 
3051
DOI : 10.1103/PhysRevB.10.3038
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
68
155427 
DOI : 10.1103/PhysRevB.68.155427
Vesseur E. J. R.
,
de Waele R.
,
Kuttge M.
,
Polman A.
(2007)
“Direct observation of plasmonic modes in Au nanowires using highresolution cathodoluminescence spectroscopy”
Nano Lett.
9
2843 
2846
Schmidt M. A.
,
Russell P. St. J.
(2008)
“Longrange spiralling surface plasmon modes on metallic nanowires”
Opt. Express
16
13617 
13623
DOI : 10.1364/OE.16.013617
Sanders A. W.
,
Routenberg D. A.
,
Wiley B. J.
,
Xia Y.
,
Dufresne E. R.
,
Reed M. A.
(2006)
“Observation of plasmon propagation, redirection, and fanout in silver nanowires”
Nano Lett.
6
1822 
1826
DOI : 10.1021/nl052471v
Hong S.K.
,
Nam S. W.
,
Yeon K.H.
(2007)
“Coherent exciton states of excitonic nanocrystalmolecules”
Phys. Rev. B
76
115330 
DOI : 10.1103/PhysRevB.76.115330