J. Electromagnetic Analysis & Applications, 2010, 2, 672-676
doi:10.4236/jemaa.2010.212088 Published Online December 2010 (http://www.SciRP.org/journal/jemaa)
Copyright © 2010 SciRes. JEMAA
Theory of Carbon Nanotubes as Optical Nano
Waveguides
Afshin Moradi1,2
1Department of Nano Science, Kermanshah University of Technology, Kermanshah, Iran; 2Department of Nano Science, Institute for
Studies in Theoretical Physics and Mathematics (IPM), Tehran, Iran
E-mail: a.moradi@nano.ipm.ac .ir
Received November 5th, 2010; revised November 25th, 2010; accepted December 1st, 2010.
ABSTRACT
The propagation of surface plasmon waves in metallic single-walled carbon nanotubes is analyzed within the frame-
work of the classical electrodynamics. The conduction electrons of the system are modelled by an innitesimally thin
layer of free-electron gas which is described by means of the semiclassical kinetic theory of the electron dynamics. The
effects of the energy band structure is taken into accoun t and a more accura te dispersion relation fo r surface plasmon
oscillations in the zig-zag and armchair nanotubes of metallic cha ra c ter is obtained.
Keywords: Carbon Nanotubes, Dispersion Relation , Plasmon Oscillations
1. Introduction
With the discovery by Iijima [1] of carbon nanotubes
(CNTs) structures, a new class of materials with a re-
duced dimensionality has been introduced. Metallic
CNTs are considered suitable candidates in the eld of
plasmonics as new plasmonics waveguides [2-12]. These
new plasmonic waveguides can be built by some simple
and well-known methods such as CVD [12].
By using the classical electrodynamics and a semi-
classical kinetic theory, Slepyan et al [2], derived the
dispersion relation of surface waves in single-walled
carbon nanotubes (SWCNTs) and for the first time,
found that CNTs can be used as a nano waveguide for
controlling electro magnetic wave propagation in speci-
fied frequency ranges (for instance, infrared and optical).
In particular, by solving Maxwell and hydrodynamic
equations, the propagation of electromagnetic waves in
SWCNTs is studied in Reference [4,5] and it has been
shown that dispersion behaviors of the plasma waves
with TM and TE modes are quite similar. However, the
hydrodynamic theory, without any effects of the energy
band structure to be taken into account, can not be valid
enough for the investigations of plasmon waves propaga-
tion in CNTs.
SWCNTs are quasi one-dimensional material, which
could be regarded as a rolled-up graphene layer (i.e., a
mono-atomic layer of graphite) in the cylindrical form. It
has a radius of a few nanometers and lengths up to cen-
timeters. A graphene layer is a semi-metallic material.
Nevertheless, when a graphene layer is rolled up it may
become either metallic or semiconducting, depending
on its geometry. Figure 1 shows a graphene layer. The
geometric structure of a SWCNT is uniquely deter-
mined by the chiral vector
12
R,manam n ,
where m and n are integers, and a1 and a2 are the ele-
mentary vectors of the dimensional graphene lattice
(see Figure 1). The tube radius of the CNTs is given by
22
0,
2
c
a
rmmn

n
where 0 is the lattice constant of the graphite
sheet and 0
0
1.42b
=3abA
is the distance between the near-
est-neighboring carbon atoms. A SWCNTs is metallic if
3mn q
, where 0,1,2,q
q
Thus, armchair nano-
tubes are always metallic, whereas zig-zag nanotubes are
metallic only if 3m
with 1, 2,q
In this Letter, we study the energy band effects on the
dispersion relation of the surface plasmon waves in
SWCNTs of metallic character, by using the semiclassi-
cal kinetic theory of the electron dynamics. In compari-
son with previous investigatio ns [4,5] that focus on plas-
mon wave oscillations in a cylindrical electron gas as a
simple model of metallic tubes, present work stresses on
more exact analysis of geometrical effects, including the
radius and the chiral angle of the nanotube.
Theory of Carbon Nanotubes as Optical Nano Waveguides 673
Figure 1. Graphene sheet used in forming CNTs (The dots
illustrate the carbon atoms positions). The lattice basis vec-
tors are a1 and a2. The chiral vector
, where m and n are integers. The nanotubes
with are called zig-zag, those withare called
armchair and those with are called chiral.
12
Rma na
mn

,mn
n0mn
2. Formulation of the Problem
Let us consider both zig-zag and armchair
nanotubes as infinitesimally thin and infinitely
long cylindrical shells of radius c with its axis along
the -direction and regard the CNT to consist of
(,0)m
r
(,)mm
z
-electrons superimposed with equilibrium densities
(per unit area) 0. We assume that in equilibrium the n
-electron fluid has no velocity and is the perturbed
density (per unit area) of fluid, produced by the
n
-electrons themselves under the action of the electric
field generated by the fixed positive ions of the lattice.
Hydrodynamic theory describes electronic motion in
terms of two dynamical variables, namely the electron-
density fluctuation, and
(,)nxt
uxt
,=

,
z
uu
The
basic equations in this linearized hydrody- namic model
are the equations of motion, the equations of continuity,
 
0
,=
nxt nuxt
t

,
,
(1)
   
0
00
,=, ,
eff
uxt n
nn xteExtnu xt
tm

 

,,
ˆ
(2)
where

ˆ
,=zz
ExtEe Ee
eff
m
is the tangential com-
ponent of the electromagnetic field, e is the element
charge, is the effective mass of the
-electrons and
1
=(/)(/ )
z
ezae

differentiates only tangen-
tially to the nanotubes surface. In the right-hand side of
Equation (2), the first term arise from the internal
interaction force in the fluid with that is the
square of the speed of propagation of density
disturbances in a uniform 2D homogeneous electron
fluid. The second term is the force on
2
=/2
F

-electron fluid
due to the tangential component of the electric field,
evaluated at the nanotube surface and the last
term represents the effects of the scattering of the
electrons with the positive-charge background, where
=c
rr
being the friction coefficient. We note that this term has
been neglected in the previous works [4,5].
Let us deal with the surface plasmon waves with TM
modes only in this paper, the similar result for TE mode
can be obtained. The electric field vector can be
expanded in the following Fourier forms (,Ex
im qz

)t
 
=
(,)=,.
m
m
ExtdqE rqexpt
 




(3)
For the TM modes, the field components can be
expressed in terms of
z
m
E, and it is readily shown that
this satisfy
2
2=0
zm
m
E
r



2
1,
zm
dd
rE
rdr dr




2
/c
(4)
where and c is the light speed. The
parameter
222
=q

is a real quantity in the region
<.c
q
(5)
This means that we deal with the slow transverse
magnetic waves. By eliminating the velocity field
,uxt, one can obtain the following equation from
Equations (1 ) and (2)
 
,.
t
22
2,= ,nx
tnxtE x
t





 
0
eff
n
e
tm
qz
(6)
Upon solving Equation. (6) by means of the space-time
Fourier transforms for the induced density on
the cylindrical surface, we find ),( txn
 
=
,= ,
m
m
n xtdqNqexpit
 




m (7)
where
0
=,
z
r


m
meff
n
ie
N E
m

c
m
qE
(8)
and
2
m
22
=/
mc
iq

r. Now, we use the
appropriate boundary condition, we have
 
><
0
||=
e,
rm crrrm c rrm
cc
Er Er N
 (9)
Copyright © 2010 SciRes. JEMAA
Theory of Carbon Nanotubes as Optical Nano Waveguides
674
where 0
is the permittivity of free space and the radial
component rm and the azimuthal component Em
E
of
the electric field, are given by
 
2
=
zm
rm
dE r
q
Er idr
, (10)
and
 
2
=
mz
mq
Er Er
r
.
m
,
>,
(11)
On the other hand, the relevant solution of Equation (4)
is



0
=
zmz mcmc
Er EKrIrrr

(12)
and



0
=
zmz mcmc
Er EIrKrrr

(13)
where

m
I
x and

m
K
x are the modified Bessel
functions. Substituting Equations (12-13) into boundary
condition Equation (9), by using Equation (8), for
/cq
, one can obtain the dispersion equation as
below:


22
2222
22
=,
pcmc mc
cc
mm
irIrK
rr





r
(14)
where eff
is the eigen-frequency of
the
1/2
200
=/
pc
en rm

-electron gas layer in metallic SWCNTs. The
solutions of Equation (14) yield complex frequencies
=ri
i

. It may be observed that the imaginary part
i
will be given simply by /2
. In fact, by writing
=/2
ri
, the solutions for finite damping will be
of the form

1/2
22 2
2222
22
=
42
pcmc mc
cc
mm
rIrKr
rr .i
 


 





(15)
The friction coefficient is the inverse of the electron
relaxation time
. For CNTs it is taken as =
12
310
s
=0
[13], so we have Hz. Formally
speaking, the dispersion characteristics of the surface
waves in the system is dependent on the nanotube
geometry (including the radius and the chiral angle of the
nanotube), the wave number, the angular momentum,
and the friction coefficient. However, it is easy to find
that by increasing friction coefficient, the dispersion
curves shift to lower frequencies, so in the following we
set
12
=1/3 10
.
At this stage, from Equation (15), one can see for
investigation the dispersion characteristics of the zig-zag
and armchair SWCNTs, we have to give the values of
0. The parameter 0 takes into account the
influence of the atomic crystal field. By using the
semiclassical model of the
/eff
nm /eff
nm
-electron dynamics, Miano
and Villone [14], obtained the following estimation:
02
21,
F
eff c
n
mr
(16)
where
F
is the velocity of the electrons at the Fermi
level
00
3
=.
2
F
b
0
is a characteristic energy of the
graphene lattice and is the Planck constant; it results
that . The Equation (16) holds
for zig-zag nanotubes with , for armchair
nanotubes with and for chiral nanotubes with

02.7 3eV
=3 <60mq
6
10
<50
(0.9 1)/
Fms

m
=3q2nm
. In the range of validity of Equation (16),
the parameter decreases as the nanotube radius increases.
To see clearly the energy band effects on the
dispersion relation of the surface waves in SWCNTs, in
the following we consider long and short wavelength
limits of the Equation (15). For c, by using the
well-known asymptotic expressions [15],
r

 
==(
2
2
xx
mm
e),
I
xandKxewith thefinite m
x
x
the dispersion relation can be written approximately as
2
22 2
0
=
F
c
e
r
 
. (17)
One can see unlike the case of Equation (3) in Ref. 5,
where the dispersion relation is independent of the
geometrical effects of the tube, the right-hand side of
Equation (17) in the present work, depends strongly on
the radius of the nanotube. It is easy to find that as the
values of the nanotube radius increases the values of
c
r
decrease.
In the opposite limit 0
c
r
, where the phase
velocity of the surface plasmon is comparable to the
velocity of light, surface plasmon oscillations couple
with the electromagnetic wave and retardation effects are
present. Retardation effects on low-dimensional plasmons
are investigated in details in Reference [16]. If we
neglect the retardation effects, by using the well-known
expressions of Bessel functions,
 
1
=(
12
m
m),
x
Ix m

0
2ln1.12
()=()for0and=for=0
2
m
m
m
Kxm Kxm
xx
3
then we may obtain for =0m

1/2
2
2
0
41.123
=0, 0ln,
F
c
e
mr










(18)
Copyright © 2010 SciRes. JEMAA
Theory of Carbon Nanotubes as Optical Nano Waveguides 675
that is a quasiacoustic mode and for , we get 0m
2
22
22
0
=
F
cc
e
m
rr

2
.m
(19)
which is also quite sensitive to the geometric of the
nanotube. Comparing the long-wavelength and
short-wavelength limits, it can be seen that the energy
band structure play an important role in the dispersion
relation, for all values of wavelength.
To better understand the energy band effects on the
dispersion relation of the plasmon waves in the system,
we illustrate in Figure 2, the dependence of the
frequency
on the variable for different nanotube
geometries with and 1. One can see that the
dispersion curves
q
=0m
for the zig-zag (27,0) with radius
and armchair (15,15) with radius
are largely similar. This means that behaviors
of the plasmon waves are not sensitive to the types of
metallic nanotubes with same radius.
=1.056
c
r
1.017 nm nm =
c
r
In order to study the mechanism of exciting plasmon
waves on the SWCNTs for future waveguide usage, we
plot the speed lines of three electron beams in Figure
2(a), by considering the expression =q
/s
. As seen,
when the electron beam velocities locate in the range
(i.e., the velocity of
the electron beam can be equal to the phase velocity of
the surface plasmon modes), the electron beam is in
synchronization with the surface waves, and they interact
with each other and instability occurs between them.
Thus, we conclude that surface waves in the system can
only be excited by applying some relativistic electron
beams with the speed of about .
6
0.9310/<< 410/ms msv

6
10
6
m
Finally, let us look at the spatial extension of the
electromagnetic field associated with the surface-
plasmon polariton as shown in Figure 3, for a nanotube
with radii c. The attenuation length is determined at
long- wavelength s (where retardation effects are present)
by means of the penetration depth. However in the
nonretarded surface plasmon condition, i.e.,
r
/qc
q,
the penetration depth in the system is thereby
leading to a strong concentration of the electromagnetic
surface-plasmon field near the interface. This means that
the dispersion curves in Figures 2 and 3, can show the
relation between frequency and the inversion penetration
depth. From Figures 2(a) and 2(b), one can see that, as
increasing the wave frequency, the radial penetration
depth of the TM surface modes decrease
1/
3. Conclusions
In summary, a theoretical model based on the classical
electrodynamics and linearized hydrodynamic theory is
employed to describe the plasmon wave propagation on
the surface of the metallic SWCNTs, where the effects of
the energy band structure is taken into account. It has
been found that the nanotube geometry play an important
(a)
(b)
Figure 2. The dispersion curves from the surface waves for
different nanotube geometries, for (a) m = 0 and (b) m = 1,
respectively. The straight dotted lines correspond to differ-
ent v for = vq.
Figure 3. Surfac e TM modes of a nanotube as a function of
the radial coordinate r.
Copyright © 2010 SciRes. JEMAA
Theory of Carbon Nanotubes as Optical Nano Waveguides
Copyright © 2010 SciRes. JEMAA
676
role in the dispersion relation of the surface waves, for
all value of wavelength. Also, numerical results show
that behaviors of the plasmon waves are not sensitive to
the types of metallic nanotubes with same radius. In
addition, the results obtained make us believe that the
hydrodynamic theory in conjunction with semiclassical
model is available and appropriates for studies of the
plasmon wave oscillations in CNTs, especially for dif-
ferent nanotube geometries.
REFERENCES
[1] S. Iijima, “Helical Microtubules of Graphitic Carbon,”
Nature, Vol. 354, No. 6348, 1991, pp. 56-58.
[2] G. Ya. Slepyan, S. A. Maksimenko, A. Lakhtakia, O. M.
Yevtushenko and A. V. Gusakov, “Electronic and
Electromagnetic Properties of Nanotubes,” Physical
Review B, Vol. 57, 1998, pp. 9485-9497.
[3] M. V. Shuba, S. A. Maksimenko and A. Lakhtakia,
“Electromagnetic Wave Propagation in an Almost
Circular Bundle of Closely Packed Metallic Carbon
Nanotubes,” Physics Review B, Vol. 76, No. 15, 2007, pp.
1-9.
[4] L. Wei and Y. N. Wang, “Electromagnetic Wave
Propagation in Single-Wall Carbon Nanotubes,” Physics
Letters A, Vol. 333, 2004, pp. 303-309.
[5] H. Khosravi and A. Moradi, “Comment on:
Electromagnetic Wave Propagation in Single-Wall
Carbon Nanotubes,” Physical Letters A, Vol. 364, 2007,
pp. 515-516.
[6] A. Moradi, “Guided Dispersion Characteristics of
Metallic Single-Walled Carbon Nanotubes in the
Presence of Dielectric Media,” Optics Communications,
Vol. 283, No. 1, 2010, pp. 160-163.
[7] D. J. Mowbray, Z. L. Miskovic, F. O. Goodman and Y.-N.
Wang, “Interactions of Fast Ions with Carbon Nanotubes:
Two-Fluid Model,” Physical Review B, Vol. 70, No. 19,
2004, pp. 1-7.
[8] D. J. Mowbray, S. Segui, J. Gervasoni, Z. L. Miskovi,
and N. R. Arista, “Plasmon Excitations on a Single-Wall
Carbon Nanotube by External Charges: Two-
Dimensional Two-Fluid Hydrodynamic Model,” Physical
Review B, Vol. 82, No. 3, 2010, pp. 1-14.
[9] C. Javaherian and B. Shokri, “Guided Dispersion
Characteristics of Metallic Single-Wall Carbon
Nanotubes,” Journal of Physics D, Vol. 42, No. 5, 2009,
pp. 1-6.
[10] L. Liu, Z. Han and S. He, “Novel Surface Plasmon
Waveguide for High Integration,” Optics Express, Vol.
13, No. 17, 2005, pp. 6645-6650.
[11] E. Ozbay, “Plasmonics: Merging Photonics and
Electronics at Nanoscale Dimensions,” Science, Vol. 311,
No. 5758, 2006, pp. 189-193.
[12] M. Su, B. Zheng and J. Liu, “A Scalable CVD Method
for the Synthesis of Single-Walled Carbon Nanotubes
with High Catalyst Productivity,” Chemical Physics
Letters, Vol. 322, No. 5, 2000, pp. 321-326.
[13] Y. Miyamoto, S. G. Louie and M. L. Cohen, “Chiral
Conductivities of Nanotubes,” Physical Review Letteras,
Vol. 76, No. 12, 1996, pp. 2121-2124.
[14] G. Miano and F. Villone, “An Integral Formulation for
the Electrodynamics of Metallic Carbon Nanotubes
Based on a Fluid Model,” IEEE Transaction on Antennas
and Propagation, Vol. 54, No. 10, 2006, pp. 2713-2724.
[15] M. Abramowitz and I. A. Stegum, Handbook of
Mathematical Funcyions, Dover, New York, 1965.
[16] L. Wendler and T. Kraft, “Retardation Effects on Intra-
and Intersubband Plasmons in Quantum Wells and Their
Manifestations in Grating-Coupler-Assisted Optical
Transmission,” Physical Review B, Vol. 60, No. 24, 1999,
pp. 16603-16610.