Progress in Nano-Electro Optics III-Springer (2004)
Progress in Nano-Electro Optics III-Springer (2004)
Progress in Nano-Electro Optics III-Springer (2004)
optical sciences 96
founded by H.K.V. Lotsch
Editor-in-Chief
William T. Rhodes Ferenc Krausz
Georgia Institute of Technology Vienna University of Technology
School of Electrical and Computer Engineering Photonics Institute
Atlanta, GA 30332-0250, USA Gusshausstraße 27/387
E-mail: [email protected] 1040 Wien, Austria
E-mail: [email protected]
and
Editorial Board Max-Planck-Institut für Quantenoptik
Hans-Kopfermann-Straße 1
Toshimitsu Asakura 85748 Garching, Germany
Hokkai-Gakuen University
Faculty of Engineering Bo Monemar
1-1, Minami-26, Nishi 11, Chuo-ku Department of Physics
Sapporo, Hokkaido 064-0926, Japan and Measurement Technology
E-mail: [email protected] Materials Science Division
Linköping University
58183 Linköping, Sweden
Karl-Heinz Brenner E-mail: [email protected]
Chair of Optoelectronics
University of Mannheim Herbert Venghaus
Institute of Computer Engineering Heinrich-Hertz-Institut
B6, 26 für Nachrichtentechnik Berlin GmbH
68131 Mannheim, Germany Einsteinufer 37
E-mail: [email protected] 10587 Berlin, Germany
E-mail: [email protected]
Theodor W. Hänsch
Horst Weber
Max-Planck-Institut für Quantenoptik
Hans-Kopfermann-Straße 1 Technische Universität Berlin
85748 Garching, Germany Optisches Institut
E-mail: [email protected] Straße des 17. Juni 135
10623 Berlin, Germany
E-mail: [email protected]
Takeshi Kamiya
Ministry of Education, Culture, Sports Harald Weinfurter
Science and Technology Ludwig-Maximilians-Universität München
National Institution for Academic Degrees Sektion Physik
3-29-1 Otsuka, Bunkyo-ku Schellingstraße 4/III
Tokyo 112-0012, Japan 80799 München, Germany
E-mail: [email protected] E-mail: [email protected]
Motoichi Ohtsu (Ed.)
Progress
in Nano-Electro-Optics III
Industrial Applications and Dynamics
of the Nano-Optical System
123
Professor Dr. Motoichi Ohtsu
The University of Tokyo
Department of Electronics Engineering
School of Engineering
7-3-1 Hongo, Bunkyo-ku
Tokyo 113-8656
Japan
E-mail: [email protected]
ISSN 0342-4111
ISBN 3-540-21050-4 Springer Berlin Heidelberg New York
I hope that this volume will be a valuable resource for the readers and future
specialists.
Index . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . . 221
List of Contributors
S. Mononobe
1 Introduction
Recently, scanning near-field optical microscopy (SNOM) employing an aper-
tured probe has been seen as a powerful tool for nano-optical applications
due to its high spatial resolution capability down to a few tens of nanometers
or less. The resolution is mainly determined by the aperture size of the probe
and the sample–probe distance because the SNOM is a scanning probe mi-
croscope based on the short-range electromagnetic interaction between the
probe and sample, which are much smaller than the optical wavelength. It
is apparent that fabricating such probes and controlling the sample–probe
distance have been the most important factors in the development of high-
resolution SNOM.
A near-field optical interaction system can be operated in two comple-
mentary ways called illumination mode and collection mode, where the probe
acts as a selective generator of a localized optical near-field and a sensitive
scatterer of an optical near-field generated over the sample, respectively. Fig-
ures 1a and b show schematic illustrations of i-mode SNOM and c-mode
SNOM, respectively. In the c-mode SNOM where the light is incident to
the total internal reflection, the sample–probe distance can be controlled by
using near-field optical (NFO) variation as a feedback control signal. For i-
mode SNOM, a shear-force feedback technique [1, 2] has been widely used to
regulate the sample–probe distance.
In order to realize such a scattering probe and a generating probe, forming
a dielectric taper and metallizing the taper has been used as an effective
method. Tapered probes have been fabricated by the following methods:
– Etching a quartz crystal rod [3]
– Heating and pulling a glass pipette by a micropipette puller [4]
– Meniscus etching of an optical fiber [5–13]
– Heating and pulling an optical fiber [14–18]
– Selective etching of an optical fiber [19–23]
– Heating and pulling and etching of an optical fiber [24, 25]
– Meniscus etching and selective etching of an optical fiber [26, 27]
For metallizing the tapered probes, a vacuum evaporation method [28] has
been applied, where the dielectric taper is coated with an aluminum or gold
2 S. Mononobe
Fig. 1. Schematic illustrations of (a) the i-mode, (b) the c-mode, and (c) the i-
c mode. In the c-mode SNOM, the light is incident to the total internal reflection.
The three-dimensional optical near field generated and localized on the sample
surface is scattered by the probe, and part of the scattered field is collected and
detected through the probe. The principle of operation of i-mode is similar except
that the probe works as a generator of the optical near field to illuminate the sample
surface. The scattered field by the sample is collected by conventional optics. In i-c
mode, a sample is excited by an optical near field on the probe. The light generated
on the sample is scattered and collected by the probe
film by a vacuum evaporation unit except for its apex region. Among the
metallized tapered probes, fiber-type probes have high transmission efficien-
cies due to their waveguide structure, and moreover, have been demonstrated
to have the molecular sensitivity, the nanometric spatial resolution, and the
locally spectroscopic capability in i-mode SNOM, c-mode SNOM, and hy-
brid SNOM called illumination-collection mode (i-c mode). Figure 1c shows
schematic illustration of i-c mode SNOM where the apertured probe func-
Near-Field Optical Fiber Probes and the Imaging Applications 3
tions as both an optical generator exciting the sample and a scatterer of the
excited optical near-fields.
i-mode SNOM and i-c mode SNOM are often employed for near-field
imaging of dye-doped samples and spectroscopic investigation of semicon-
ductor devices. In such spectroscopic applications where one must cope with
extremely low detected power, the probe should have high throughput in i-
mode SNOM, and be highly sensitive in c-mode SNOM. Furthermore, in the
i-mode fashion, to avoid thermal damage to the sample and the probe, the
probe should be used with as low an input power as possible. Therefore, the
resolution capability and throughput of the tapered probe have to be opti-
mized depending on SNOM applications. This optimization should be done
by varying the probe shape, i.e., the cone angle. Furthermore, to obtain a
highly resolved image, one must fabricate a metallized probe with an apex
region emerging from a metal film.
We have developed several application-oriented fiber probes called protru-
sion-type probe, double-tapered probe, pure silica fiber probe, and performed
near-field imaging of biomolecules and dye-doped samples, spectroscopic
study of semiconductors, and ultraviolet applications. In this chapter, the
near-field optical fiber probes and near-field microscopy applications are de-
scribed. Section 2 discusses the basic techniques for tapering and metallizing
fibers. The protrusion-type probe and the biological imaging are described in
Sect. 3. Section 4 describes apertured probes and their application to near-
field imaging of dye-doped samples. The double-tapered probe and pure silica
fiber probe are described in Sects. 5 and 6, respectively. Table 1 summarizes
the SNOM images appeared in this chapter.
Fig. 2. Schematic illustration of the micropipette puller used for tapering an optical
fiber with a diameter of 125 µm
Fig. 3. Schematic illustrations of (a) the vacuum evaporation unit for metallizing
the pulled fiber and (b) the metallized fiber
as small as 50 nm. In the tapered portion, strong optical leaky modes are
generated due to the varied core diameter.
For i-mode SNOM, the pulled fiber with an apex diameter of about 50 nm
must be metallized except for its apex region. To metallize the fiber, the
pulled fiber is rotated while evaporating aluminum in vacuum as shown in
Fig. 3a so that the metallized probe has a thickness profile as shown in
Fig. 3b. Here, the typical radial thickness is around 150 nm. The apex region
is aluminized with a thickness smaller than the half-radial thickness due to
the throwing of evaporated vapor. The metal thickness covering the apex
region can be reduced to a quarter of the radial thickness by inclining the
rotating fiber in vacuum.
6 S. Mononobe
Fig. 4. Schematic illustrations of meniscus etching of a fiber at (a) the start, (b)
tapering, and (c) stop
ammonium fluoride (NH4 F), the cone angle can be controlled in a wide region
from 20◦ to 180◦ for an apex diameter less than 10 nm. Furthermore, selective
etching is the most highly reproducible technique among the three tapering
techniques. This method can be applied to any single-mode fiber produced
by vapor-phase axial deposition (VAD) [40].
Figure 5a shows a cross-sectional profile of the refractive index of a silica
fiber with a GeO2 -doped core and a pure silica clad. Here, n1 and n2 are the
refractive indexes of the core and clad, respectively. r1 and r2 are the radii
of the core and clad, respectively. On immersing the fiber in a buffered HF
solution (BHF) with a volume ratio of [40%–NH4 F aqueous sol.]: [50%–HF
acid]: [H2 O]= X:1:1 at 25◦ C, the core is hollowed at X = 0 and is tapered at
X = 10. Figure 5b shows a schematic explanation of the geometrical model
for the tapering process based on selective etching. Bright shading and dark
shading in the upper part represent the cross-sectional profiles of the fiber
before and after the etching with an etching time τ , respectively. τ is the
etching time required for making the apex diameter zero, θ and d are the cone
angle and apex diameter of the probe, respectively. The lower part shows the
dissolution rates R1 and R2 of the core and clad, respectively. Here, R1 < R2 .
Assuming that the dissolution rates R1 and R2 are constant within the core
and clad regions, respectively, relations between the cone angle θ, the length
L of the tapered core, and the apex diameter d are represented by
8 S. Mononobe
θ R1
sin = , (1)
2 R2
r1 − d/2
L= , (2)
tan(θ/2)
and
2r1 (1 − T /τ ) (T < τ ) ,
d(T ) = (3)
0 (T ≥τ ) .
Here, the etching time τ that is required for making the apex diameter zero
is expressed as
r1 R 1 + R 2
τ= . (4)
R1 R2 − R1
When a fiber with a GeO2 -doped core and a pure silica clad is immersed in
BHFs with volume ratios of X:1:1, the core region is hollowed in X < 1.7,
and is tapered in X > 1.7 as shown in Fig. 5. Based on (1), the cone angle
is determined by the dissolution rate ratio of the core and clad. Figure 6a
shows variations of the dissolution rates R1 and R2 and the cone angle θ
as a function of X. Here, the fiber used was produced by vapor-phase axial
deposition (VAD) so as to have an index difference of 2.5%. The dissolution
rate ratio R1 /R2 decreases with increasing NH4 F volume ratio of X, and
Fig. 6. (a) Dependencies of the dissolution rate R1 (broken curve) of the GeO2 -
doped core, R2 (dotted line) of the pure silica clad, and the cone angle θ (solid curve)
on X. Here, the value of ∆n defined as (n21 − n22 )/2n21 is 2.5%; (b) Dependency of
θ on ∆n at X = 10
Near-Field Optical Fiber Probes and the Imaging Applications 9
Pure silica core fiber has high transmittance in a wide region from near
ultraviolet to near infrared (see Fig. 33). To make an index difference between
the pure silica core and the clad, the clad is often made of fluorine-doped
silica. We produced a pure silica core fiber with the fluorine-doped clad and
an index difference −0.7%.2 By immersing the fiber in BHFs with volume
ratios 1.7:1:1 and 10:1:1, the obtained probe has cone angles of 62◦ and 77◦ ,
respectively. The fiber could not be hollowed at any volume ratios of BHFs.
In selective etching of a pure silica fiber, varying the index difference or the
fluorine-doping ratio controls the cone angle. Among single-mode fibers, only
1500-nm pure silica fiber (PSF) with an index difference of −0.3% and core
diameter of 10 µm is commercially available.
For the shear-force feedback technique, the clad diameter of the probe is
one of the main parameters governing the resonance frequency of dithering a
probe. In our experience, tapered probes with a clad diameter of 20–30 µm
are usually needed to obtain a proper resonance frequency, i.e., 20–40 kHz.
Figure 7 schematically shows an etching method to fabricate tapered probes
2
For the following discussion, we define the relative refractive index difference ∆n
of doped glass to pure silica, which is expressed as (n22 −n21 )/2n22 and (n21 −n22 )/2n21
for a pure silica core fiber and a pure silica clad fiber, respectively.
10 S. Mononobe
with the reduced clad diameter. Here, r2 is the clad radius before etching, D
is the clad diameter of the probe, θ is the cone angle of the tapered core, d
is the apex diameter, and L is the length of the tapered core. This method
involves two steps: (A) reducing the clad thickness and (B) tapering the core.
In step A, the fiber is immersed in BHF with volume ratio of [NH4 F
aqueous solution (40%)]: [HF acid (50%)]: [H2 O]=1.7:1:1 for an etching time
TA . Then, the clad diameter is reduced to [2r2 − 2R2A TA ]. The core end is
kept flat since the dissolution rates R1A and R2A of the core and clad are
equal.
In step B, the fiber is selectively etched in XB :1:1 (where XB >1.7) for
etching time TB . The core is tapered with the cone angle θ represented by
θ R1B
sin = . (5)
2 R2B
Near-Field Optical Fiber Probes and the Imaging Applications 11
Fig. 8. SEM micrographs of (a) a shoulder-shaped probe and (b) the magnified
apex region. D = 25 µm; θ = 20◦ ; d < 10 nm
To make the apex diameter zero, the etching time TB must be longer than
the time τ , which is given by
r1 R1B + R2B
τ= . (6)
R1B R2B − R1B
The clad diameter D is represented by
Step Volume ratio of BHF Etching time (min) Dissolution rates (µm h−1 )
A 1.7:1:1 86 R1 = R2 = 30
B 10:1:1 75 R1 = 1.1, R2 = 6.5
12 S. Mononobe
The uncoated fiber probe with a nanometric apex can be employed for the
c-mode SNOM as shown in Fig. 1b. Although the principal factor governing
the resolution is the small apex size, the resolution is affected by propagat-
ing components and the scattering of low spatial frequency components of
the optical near field. To obtain highly resolved c-mode images, such un-
wanted scattered light must be suppressed. To perform high-resolution near-
field imaging, we developed a fiber probe with a protruding tip from a metal
film, which is called a protrusion-type probe [41]. Figure 9a shows the cross-
sectional schematic illustration of the protrusion-type probe. Here, θ, d, df ,
and tM represent the cone angle, the apex diameter, the foot diameter, and
the thickness of the metal film, respectively. This probe scatters frequency
components between (1/df ) and (1/d) and suppresses the generation of com-
ponents of less than (1/df ). To realize such a probe, a selective resin-coating
(SRC) method [41] was developed, which can be applied to tapered probes
fabricated by the etching method as shown in Fig. 7. Apertured probes are
also effective for suppressing unwanted light. However, it is difficult to pro-
duce apertured probes with an apex diameter of a few tens of nanometers by
a vacuum evaporation unit. If applying the vacuum evaporation method to
a fiber tip with a cone angle of 20◦ , an entirely coated fiber probe as shown
in Fig. 9b is obtained. The thickness tA of metal covering the apex could not
be made less than 40 nm for tM = 150 nm.
Figure 10 shows a schematic diagram of the SRC involving four steps: (A)
metal coating, (B) resin coating, (C) preferential etching of metal covering
the apex region, (D) removing the resin. Figure 11 shows a scanning electron
micrograph of the top region of the fabricated protrusion-type probe. The
probe has a protrusion from a gold film and a foot diameter df of less than
30 nm.
Before applying the SRC method, a tapered probe was fabricated with a
clad diameter D = 45 µm and a cone angle θ = 20◦ based on selective etching.
In step A, the probe was coated with 120-nm thick gold by a magnetron
sputtering unit. In step B, the probe is dipped in an acrylic resin solution
and then removed from the acrylic solution with a withdrawal speed of VD =
5 cm s−1 . Using an acrylic resin solution with a low viscosity coefficient of
11 cP, we succeed in making a fiber tip protruding from resin film. The density
of the acrylic solution is 0.85 g cm−3 at 25◦ C. In step C, the probe was etched
for 2 min in a solution KI-I2 -H2 O, mixed with a weight ratio of 20:1:400 and
diluted 50 times with deionized water. In step D, the acrylic film was removed
by immersing the probe in acetone.
Fig. 10. Schematic diagram of the selective resin-coating method for fabricating
a protrusion-type probe; D, clad diameter of the probe; r1 , core radius; θ, L, cone
angle and length of the tapered core, respectively; VD , withdrawal speed of the fiber
from the resin solution; LR , length of the tapered core where the resin is coated;
df , foot diameter of the protrusion
14 S. Mononobe
Fig. 12. Dependencies of the foot diameter df on (a) the clad diameter D for a
withdrawal speed of 5 cm s−1 and (b) the withdrawal speed VD for a clad diameter
of 45 µm
Figures 12a and 12b show the dependencies of the foot diameter on the
clad diameter and the withdrawal speed, respectively. The foot diameter can
be controlled by varying the viscosity of the resin solution, withdrawal speed,
and the clad diameter.
Fig. 14. Variation of the optical near-field intensity as a function of the sample–
probe distance, obtained with c-mode SNOM employing a protrusion-type probe
Sample
The samples were fabricated by fixing the separated filaments from the body
on a hydrophilic glass substrate. A single filament has a typical length of a
few micrometers and a diameter of around 25 nm as shown in Fig. 15b. The
sample substrate was mounted on the prism of the optical system as Fig. 13.
Fig. 15. (a) Schematic illustration of the fabricated protrusion-type fiber probe
with cone angle of 20◦ , a foot diameter of 30 nm, and apex diameter less than 10 nm;
(b) A transmission electron micrograph of the Salmonella flagellar filaments. The
diameter of a flagellar filament, as indicated by arrows, are 25 nm; (c) A c-mode
SNOM image of Salmonella flagellar filaments. The arrows at the top right-hand side
represent the directions of wavevector k and electric field vector E; (d) Perspective
view of the near-field image of Salmonella flagellar filaments in water. The bright
region, as indicated by arrows, is 50 nm. The arrows at the top right-hand side are
defined as the directions of wavevector k and electric field vector E
optical imaging in water was also done by c-mode SNOM under NFO in-
tensity feedback control as in air. The filament sample was fixed on a glass
substrate and was fitted with an acrylic ring. The columnar bath with an
inner diameter of 10 mm and a height of 2 mm was filled with water during
the observation. The sample–probe distance was estimated to be less than
30 nm. The bright segments correspond to the filaments. The FWHM of the
18 S. Mononobe
Zvyagin et al. [31] isolated microtubules from the brain of a pig, and observed
with c-mode SNOM using a protrusion-type probe with a foot diameter of
30 nm. For this observation, a c-mode SNOM with s-polarized incident light
was used. To regulate the sample–probe distance, a shear-force feedback sys-
tem [1, 2] was added to c-mode SNOM. In the shear-force feedback system,
the differential force and frictional force are detected by dithering the probe
at its resonance frequency at some amplitude, and the amplitude of dithering
is measured by an optical interference technique. Figure 17 shows the varia-
tion of the amplitude of dithering of the protrusion-type probe as a function
of the sample–probe distance. The inset is the resonance frequency spectrum.
Figure 18a shows the transmission electron micrograph of microtubules.
The dark region seen in the middle of a strand corresponds to microtubules.
The width of the microtubules, as indicated by arrows, is 25 nm. Figure 18b
shows a shear-force topographic image of microtubules. The sample–probe
distance was kept at less than 5 nm during the scanning. Here, the long
bright structure and the round ones correspond to a microtubule and pro-
tein aggregates, respectively. From this high-resolution imaging result, the
Near-Field Optical Fiber Probes and the Imaging Applications 19
Fig. 17. Variation of the dithering amplitude of the probe as a function of the
sample–probe distance. The inset is the resonance frequency spectrum of the dither-
ing probe
Fig. 18. (a) Transmission electron micrograph of a microtubule sample; (b) Shear-
force topographic image of the microtubule sample; (c) c-mode SNOM image in the
rectangular region in (b). The FWHM of the dark-line structure is approximately
40 nm
For this near-field imaging, other probe types were also used, which we
fabricated by pulling and etching the highly GeO2 -doped fiber. In the heating-
and-pulling process, the fiber was tapered with the highest pulling strength
by a commercial puller (Sutter, P-2000) so as to have a dip structure at the
taper end. By the etching in a BHF with a volume ratio of 10:1:1, the dip
was flattened, and then, the core is protruded from the pure silica clad. The
tapered fibers were aluminized by the vacuum evaporation method. Some of
the probes had throughputs higher than the protrusion-type probes. However,
we could not fabricate high-throughput probes with reproducibility because it
is difficult to control the sizes of the dip generated at the pulling. Figures 20a
Near-Field Optical Fiber Probes and the Imaging Applications 21
Fig. 20. Photoluminescence spectra on maintaining the protruding tip (a) 200 nm
above the quantum dots, (b) less than 10 nm above the quantum dots, (d) less than
10 nm above the SiO2 mask; (c–f ) Monochromatic photoluminescence images at the
energy regions labeled A–C, respectively, and simultaneously obtained topographic
image
Fig. 21. Scanning electron micrographs of (a) a pulled and etched probe and (d)
its magnified view
closed circles and open triangles correspond to the variations for the cone
angles of 20◦ and 50◦ , respectively. To collect the scattered light from the
apex of the probe, an objective lens with a numerical aperture of 0.4 was
inserted in front of the photodetector. At dF = 200 nm, The throughput of
the probe with a large cone angle of 50◦ is 10 times larger than that with
θ = 20◦ . Increasing the cone angle enhances the throughput. For a near-field
spectroscopic application requiring high throughput, the tapered probe with
a large cone angle, i.e., 50–90◦ should be fabricated by a selective-etching
method.
Fig. 23. (a) Ag-MgF2 -Al-coated fiber probe; (b) Variation of the optical near-field
intensity as a function of the sample–probe distance
24 S. Mononobe
a cone angle of 30◦ . We call this tapered fiber a pencil-shaped probe. The
method for fabricating the pencil-shaped probe with a metal-dielectric-metal
coating involves two steps: first sharpening of an optical fiber a into pencil-
shaped structure and then coating in a vacuum evaporation unit.
The pencil-shaped probe with a cone angle of 30◦ and an apex diameter
less than 10 nm was fabricated with almost 100% reproducibility based on
selective etching of a multistep index fiber with the GeO2 -doped silica core,
pure silica clad, and a F-doped silica support. This fiber is specially designed
and fabricated by the VAD method to produce pencil-shaped probes and
triple-tapered probe. For details of the etching process, refer to the Appendix.
The double-stranded plasmid DNA with a ring structure (pUC18 and 2868
base pairs) was used as the sample. It was diluted with distilled water to a
final concentration of 5 ng/µL. A quantity of 2 µL of this was dropped onto
Near-Field Optical Fiber Probes and the Imaging Applications 25
Figure 24a shows the image of DNA on the ultrasmooth sapphire surface
obtained using a noncontact-mode atomic force microscope [43]. The bright-
ring structures correspond to DNA molecules. The observed width of a single
strand, as indicated by open arrows, is estimated to be around 10 nm.
The experimental system of c-mode as appears in Fig. 13 was used for
near-field optical imaging of single DNA molecules. The DNA sample is
mounted on the prism, and the optical near field is sensitively picked up by
the Ag-MgF2 -Al-coated probe. Figure 24b shows the c-mode SNOM image,
obtained under near-field optical intensity feedback control. The pixel size is
5 nm × 5 nm. The arrows at the top left-hand side of the figure indicate the
propagation k and electric field E vectors of the incident light with a 680 nm
wavelength. The coiled loop and single strands of DNA are observed in this
figure. The FWHM of the intensity variation at the portion, as indicated by
the line, is around 20 nm. Figure 24c shows the magnified view of the single
strands of DNA. The pixel size is 2 nm × 2 nm. The FWHM, as indicated
by the arrows, is about 4 nm. This nanometric spatial resolution reflects the
special care taken during preparation of the sample and of the probe to avoid
unwanted scattered light.
We tried to obtain near-field optical images of DNA molecules by em-
ploying the protrusion-type probe with a foot diameter of 30 nm. However,
successful results could not be obtained. The metal-dielectric-metal-coated
probe tip with a layer of thin silver coating acts as a scattering probe much
more sensitive than a dielectric tip protruding from metal.
26 S. Mononobe
Fig. 24. (a) AFM image and (b, c) near-field optical images of a DNA sample,
obtained with c-mode SNOM employing the fiber probe in Fig. 23a
Fig. 26. (a) Cross-sectional illustration of the fabricated metallized fiber by the
vacuum evaporation method. Its top region enclosed by the dotted rectangular was
removed by rubbing with the substrate under shear-force feedback control before
the near-field imaging; (b) near-field optical image of neurons of rats. The magnified
view of (c) shear-force image and (d) i-mode SNOM images of the same area
dark arrows and the bright arrow correspond to the neural process and the
cell body, respectively. Before this imaging, the metal covering the top region
of the probe, enclosed by the dotted rectangular in Fig. 26a, was removed
by rubbing with the sample substrate under shear-force feedback control.
Figures 26c and 26d correspond to the simultaneous shear-force and near-
field images of the magnified region, as shown by the long dark arrow in the
Near-Field Optical Fiber Probes and the Imaging Applications 29
neuron process in Fig. 26b. The fringe-like structures are seen in the near-field
image of Fig. 26d, however, they do not appear in the topographic image.
This shows that the fringe-like structures labeled with toluidine blue lie just
underneath the cell membrane on the surface. The structures are identified as
microtubules, which are the main constituent elements of the neural process.
The FWHM of near-field intensity variation at the portion, as indicated by
the line in Fig. 26d, is estimated to be 26 nm, which agrees with the typical
diameter of 25 nm of a single microtubule.
Recently, another apertured probe with a large cone angle, i.e., 50–90◦
has been chosen for performing i-mode and i-c mode imaging applications
including single-dye-molecule detection because the probe has a much higher
throughput than the tapered probe as shown in Fig. 26a. Hosaka et al. [35]
fabricated a double-tapered probe [44] with a cone angle of 90◦ and an aper-
ture diameter of a few tens of nanometers as shown in Fig. 27a, and succeeded
in obtaining near-field fluorescence imaging of dye molecules with a high reso-
lution of around 10 nm. Figure 27b shows the near-field fluorescence image of
single Rhodamin 6G molecules dispersed on the quartz glass at 700 ± 20 nm.
In this experiment, the dye was excited by a He-Ne laser with a 633 nm
wavelength. The bright spots correspond to emitting single dye molecules.
The estimated resolution is around 10 nm. For details of the double-tapered
probe see Sect. 6.
Fig. 27. (a) Schematic illustration of the double-tapered probe with a cone angle
of 90◦ ; (b) i-c mode near-field fluorescence image of single Rhodamin 6G molecules
dispersed on a quartz substrate at 700 ± 20 nm
30 S. Mononobe
the wavelength size λ, and d and df are the apex and foot diameters of the
protrusion, respectively.
As discussed in Sect. 3.6, the first cone angle of the double-tapered probe
has to be increased to enhance the throughput of the probe. Since the light
entering the tapered core is strongly attenuated by the metal film in the
subwavelength cross-sectional portion, by making the base diameter the value
of λ or more, the double-tapered probe with the first cone angle of 90◦ has a
transmission coefficient as high as the single-tapered probe with a cone angle
of 90◦ . Based on FDTD simulation [45], the double-tapered probe can have a
collection efficiency higher than the single-tapered probe. The double-tapered
structure is suitable for coupling the scattered light into the core.
The method for fabricating the probe involves two steps: (A) controlling the
clad diameter and forming the first taper; (B) forming the second taper,
which are based on selective etching as described in Sect. 2.3. Figure 29a
shows the method schematically. In step A, by immersing the fiber in BHF
with XA :1:1 (where XA > 1.7), the core is protruded from the clad with the
cone angle θ1 represented by
θ1 R1A
sin = . (8)
2 R2A
In step B, immersing the fiber in BHF with 10:1:1 the core is tapered with
cone angles of θ1 and θ2 , where
θ2 R1B
sin = . (9)
2 R2B
The base diameter dB is represented by
TB
dB = 2r1 (1 − ) (TB < τ ) , (10)
τ
which is controlled by varying the etching time TC . Here, τ is defined by (4).
The clad diameter D is given by (7).
In actual conditions as shown in Table 4, a probe with a double-tapered
core with cone angles of θ1 = 90◦ and θ2 = 20◦ , and a base diameter of
dB = 500 nm was obtained. Figures 29b and c show SEM micrographs of
the tapered core region of a double-tapered probe and its magnified top
region.
32 S. Mononobe
Fig. 29. (a) Etching method for fabricating a double-tapered probe. The fiber is
consecutively etched in buffered HF with volume ratios of XA :1:1 (1.7 < XA < 10)
and 10:1:1 in steps A and B, respectively. Here, 2r1 and 2r2 are the core and
clad diameters, respectively. θ1 and θ2 are the first and second cone angles of the
double-tapered probe; SEM micrographs of (b) a double-tapered core and (c) the
magnified top region. θ1 = 90◦ ; θ2 = 20◦ ; dB = 500 nm
Fig. 30. Schematic illustrations of (a) a double-tapered probe and (b) a i-c mode
SNOM for investigating the lateral p-n junction
34 S. Mononobe
Fig. 31. (a) Shear-force topographic and (b) i-c mode SNOM images of the GaAs
lateral p-n junction
have been obtained with i-c mode SNOM employing a double-tapered probe
with an aperture diameter of 20–30 nm, and a resolution of 20–30 nm has
been achieved. For details of these studies, see Saiki [47].
Fig. 32. (a) Near-field Raman spectra of a tabular polydiacetylene single crystal;
(b) Two-dimensional mapping the relative near-field intensity at peak S to peak L
Fig. 33. Optical transmission spectrums of a GeO2 -doped core fiber with a re-
fractive index difference of 1% (dotted curve) and a pure silica core fiber (solid
curve)
Fig. 34. Pure silica core probes with (a) the conical core protruded from the flat
clad end (b) tapered clad and core, and (c) triple-taper structure
Fig. 35. Scanning electron micrograph of a tapered fiber with a pure silica core,
obtained by immersing a 1300-nm fiber with a core diameter of 10 µm for 240 min
in a BHF with a volume ratio of 10:1:1. The cone angle is 120◦
The inset in Fig. 36 shows the index-difference profile of the double-core fiber.
The diameters of sections 1–5 are 100 nm, 2.9 µm, 40 µm, 42 µm, and 125 µm,
respectively. Sections 3–5 correspond to the clad region. Sections 3 and 5
are made of low fluorine-doped silica with an index difference of −0.35%.
Section 4 is made of the high fluorine-doped silica with an index difference of
−0.94%. The core region consists of sections 1 and 2, which are made of GeO2 -
doped silica and pure silica, respectively. The index difference of section 1 to
section 2 is estimated to be around 0.7%. The diameter of section 1 is only
100 nm. Since 99.9% of the cross-sectional area of the core region is occupied
by the pure silica, the double-core fiber has transmittance as high as that of
the pure silica fiber as shown in Fig. 33.
To fabricate this double-core fiber, we fabricated a preform glass rod by a
combination of vapor-phase axial deposition and plasma-activated chemical
vapor deposition and then drew the 125-µm diameter fiber by heating the
preform. Table 5 summarizes the parameters for drawing the fiber.
As shown in Fig. 36, the double-core fiber is tapered by a selective etching
method, which involves three steps: (A) hollowing of sections 1 and 4, (B)
38 S. Mononobe
Fig. 36. Etching process to fabricate a triple-tapered probe. The inset shows the
cross-sectional profile of the relative refractive index difference ∆n of a multistep
index fiber produced to fabricate the probe
Fig. 37. Cross-sectional profile of the index difference ∆n of a multistep index fiber
produced to fabricate the probe; SEM micrographs of (a) a triple-tapered probe
and (b) its magnified apex region. θ1 = 60◦ ; θ2 = 120◦ ; θ3 = 60◦ ; d < 10 nm
cone angle θ2 , and the third cone angle θ3 are 60◦ (±1◦ ), 120◦ (±1◦ ), and
60◦ (±1◦ ), respectively. Table 6 summarizes the parameters of the etching
process for fabricating the probe of Fig. 37a. After the etching, the triple-
tapered probe was coated with a 200-nm thick aluminum film by a vacuum
evaporation method.
40 S. Mononobe
Fig. 38. Absorption (dotted curves) and photoluminescence (solid curves) spectra
of polydihexylsilane (PDHS) and n-decyl-(s)-2-methylbutyl silane (Chiral-PS), re-
spectively. These polysilanes have different backbone structures called transplanar
and 7/3 helical, respectively
Near-Field Optical Fiber Probes and the Imaging Applications 41
Fig. 39. Schematic diagram of the experimental setup for UV spectroscopic applica-
tions. Here, the photon-counting detection is performed at i-mode SNOM imaging.
The spectrum of the sample is obtained by using the CCD camera and monochro-
mator
Fig. 40. (a) Triple-tapered fiber probe with a pure silica core and (b) UV near-
fild photoluminescnece image of transplanar-type polydihexylsilane obtained with
i-mode SNOM and probe as shown in (a) and (b), respectively. The top region of
the probe, surrounded by the dotted rectangular in (b), was removed before the
imaging
42 S. Mononobe
Fig. 41. (a) Topographic image and (b) i-mode near-field photoluminescence image
of chiral-PS molecules at 325 nm
luminescence image of PDHS at 380 nm. The FWHM of the bright portion is
estimated to be around 100 nm. Such a near-field signal could not be detected
through double-tapered probes with the GeO2 -doped core due to background
originated from the luminescence of the core.
Next, by applying the above-described UV-SNOM technique to chiral-
PS, we succeeded in obtaining i-mode near-field photoluminescence images
of chiral-PS molecules dispersed on ultraflat sapphire substrate. Figure 41a
shows a topographic image of chiral-PS. The height of the observed fibril-like
structure is 7–50 nm. Since the diameter of Chiral-PS molecules has been esti-
mated to be around 1 nm, we guess the fibril-like structures are constructed
of several molecules. Figure 41b shows a UV near-field photoluminescence
image of chiral-PS molecules at 325 nm in the same area as Fig. 41a. For
the excitation of the sample, UV light with a wavelength of 300 nm was
used. The achieved spatial resolution is around 60 nm. The probe used was a
UV triple-tapered probe with a small core diameter of around 2 µm. It was
fabricated based on selective etching of a new double-core fiber, which was
produced by only vapor-phase axial deposition. Details of this fiber will be
discussed elsewhere.
Among commercial single-mode fibers, only 1.3-µm fiber (PSF) with a pure
silica core has a high transmittance at near UV. We developed a new method
of tapering a 1.3-µm PSF based on a combination of pulling and etching.
The method involves two steps: (A) heating-and-pulling the fiber by a mi-
cropipette puller and (B) etching the fiber in buffered hydrogen fluoride so-
lution as shown schematically in Fig. 43a. Figure 43b shows the magnified
top region of the tapered shape formed by step B. Here, θ and 2r1E are the
cone angle in the apex region of the fiber and the reduced core diameter at
the end of the tapered core, respectively. Figure 43c shows the cross-sectional
profile of the tapered probe obtained by increasing the etching time in step B.
Details of this profile are discussed later. The fabricated tapered fiber has the
tapered shape and the magnified top region as seen in the scanning electron
micrographs of Figs. 44a and b, respectively. Figures 44c and d correspond
to the magnified view of regions indicated by squares in b and c, respectively.
The cone angle and the apex diameter are θ = 65◦ and 10–20 nm, respec-
tively. By investigating 20 fiber samples, we obtained 80% reproducibility for
a cone angle θ = 60◦ (±5◦ ) and an apex diameter of 10–20 nm. We now
describe the method for producing this pure silica core fiber.
44 S. Mononobe
Fig. 43. (a) Method for fabricating a tapered fiber with the pure silica core. (b)
Cross-sectional profile of the top of the fiber formed by step A in (a). r1E , reduced
core diameter at the end of the tapered core; θB1 , θB2 , cone angles of the core and
clad, respectively. (c) Cross-sectional profile obtained by increasing the etching time
in step B
In step A, a PSF with core diameter of 10 µm, a clad diameter of 125 µm,
and an index difference of −0.3% was heated and pulled by a micropipette
puller (Sutter Instrument, P-2000) combined with a CO2 laser. This puller
has been mechanically adjusted to fabricate a micropipette with a diameter of
1 mm. To pull the fiber, the two puller bars were exchanged to optional ones
Near-Field Optical Fiber Probes and the Imaging Applications 45
on which a 125-µm bare fiber can be attached. In our case, to produce a probe
with a bare-portion length as small as 1 cm, a plastic-coated portion of the
fiber sample was carefully attached on one bar, and a bare portion was fixed
to another bar. The puller is adjusted with its parameters shown in Table 7.
Here, the heat parameter H (0≤H≤999) fixes the CO2 laser power. The
filament parameter F is the length of the fiber that is scanned with the CO2
laser beam. The velocity parameter V (1≤V≤255) shows the velocity of the
puller bar at the end of the heating time. The delay parameter D (0≤D≤255)
represents the delay time between the end of the heating and the beginning
of the pulling (in milliseconds). The puller is mechanically adjusted to make
the delay zero at D = 125. The pull parameter P decides the pull strength
and is controlled in a region of 0≤P≤255. If step A is repeated using the
same puller and the same fiber, some of the parameters, such as H and D,
may be changed to fabricate a pulled fiber as seen in Fig. 44a.
In step B, the fiber was etched by immersing for 10 min in a buffered
hydrogen fluoride solution (BHF) with a volume ratio of 40%-NH4 F aqueous
solution : 50%-HF acid : deionized water = 1.7:1:1. The temperature of the
BHF was 25◦ C (±0.1◦ C). The fluorine-doped clad and the pure silica core
have dissolution rates of R1 = 6.6 µm h−1 and R2 = 7.6 µm h−1 , respectively.
The pulled and etched probe with a few tens of nanometers was obtained by
stopping the etching before the core protruded from the fluoride-doped clad.
The cone angle θ is increased by increasing the etching time T . In our exper-
iments, the cone angles were 35◦ and 65◦ for etching times of 5 and 10 min,
respectively. However, we could not realize the pulled and etched probe with
a reproducibility of more than 80% due to the mechanical misalignment of
the puller.
On the other hand, once the core is exposed from the fluorine-doped silica,
the core is selectively etched due to the difference of the dissolution rate R1
of the core and the dissolution rate R2 of the clad (> R1 ). Then, the fiber
has two cone angles of θB1 and θB2 , as shown in Fig. 43c, where dB is the
base diameter of the conical core, with the angle θB1 given by
θB1 R1 θB2
sin = sin . (11)
2 R2 2
The pulled and etched probe as shown in Fig. 21a was fabricated by selective
etching and has a GeO2 tip protruding from the silica clad. In the case of
using a PSF, the cross section of the apex region was elliptical, with a longer
principal diameter of around 200 nm. Therefore, we could not employ the
selectively etched PSF as a probe. The value of dB is equal to the reduced
core diameter 2r1E of the fiber as shown in Fig. 43a. The reduced core diam-
eter is estimated to be 2 µm from a SEM image of selectively etched fibers
that were immersed in BHF with 1.7:1:1 for times longer than 10 min. It is
straightforward to obtain a minimum ratio of r1E /r1 = 1/5.
46 S. Mononobe
Fig. 44. Scanning electron micrographs of (a) the pulled and etched probe with
the pure silica core and (b) its magnified top region; (c, d) The magnified apex
regions as indicated by the squares in (b) and (c), respectively. The cone angle and
apex diameter are 65◦ and 20 nm, respectively
Fig. 45. (a) Schematic illustration of the probe fabricated by pulling and etching
a commercial fiber with the pure silica core; (b) Shear-force topographic image and
(c) i-mode near-field photoluminescence images of polydihexylsilane aggregates at
380 nm
8 Outlook
The resolution capability and throughput of SNOM have been greatly im-
proved by employing application-oriented probes such as a double-tapered
probe with an aperture diameter of a few tens of nanometers and a pure
silica fiber probe, and various imaging and spectroscopic applications have
been successfully carried out. However, the techniques for mass producing
these probes is not yet established due to the low reproducibility of the met-
allizing process for forming a nanometric aperture. Furthermore, single-mode
fiber probes should be newly developed for some imaging applications such
as polarization measurements in the ultraviolet and visible regions. If these
techniques are established, SNOM and related techniques will be further used
widely by many researchers and engineers. To realize the mass production of
near-field optical probes, the author has concentrated on the development of
a novel metallization technique based on electroless plating [51] and on the
fabrication of new single-mode fibers.
Acknowledgments
Appendix
Fabrication of Tapered Fibers Based on Hybrid Selective Etching
of a Double-Clad Fiber
The double-clad fiber [22] was developed for fabricating pencil-shaped probes
and triple-tapered probes. These probes can be fabricated with almost 100%
reproducibility based on selective etching. In this section, we describe the
structure of a triple-tapered probe and design/fabrication of the fiber.
Triple-Tapered Probe
Fig. 46. (a) Cross-sectional profiles of the single-tapered probe and (b) the triple-
tapered probe; θ, cone angle of the tapered core; λ, optical wavelength; L, length
of the tapered core with the foot diameter equal to λ; tS , skin depth of the metal;
θ1 , θ2 , θ3 , cone angles of the first, second, and third tapers, respectively; dB1 , dB2 ,
base diameters of the first and second tapers, respectively; L1 , length of the first
taper
increasing the cone angle. In the case of fabricating the tapered probe with
a large cone angle, i.e., 120◦ , one must accept the limited resolution affected
by optical leaking from the metal around the apex region of the probe. In
the triple-tapered probe, we can decrease the length L by increasing θ2 in
order to enhance the throughput. Thus, the resolution capability is increased
by simultaneously decreasing the cone angle and decreasing the first taper
length L1 to a few hundreds of nanometers, which corresponds to several
times the skin depth of metal.
Figure 47a shows the cross-sectional profile of the index difference of the
double-clad fiber. Here, the fiber involves three sections: (1) a GeO2 -doped
silica core, (2) a pure silica clad, and (3) a fluorine-doped silica support. The
index differences of sections 1 and 3 with respect to section 2 are 1.2% and
−0.7%, respectively. The radii of sections 1–3 are r1 = 0.65 µm, r2 = 13.5 µm,
and r3 = 62.5 µm, respectively. Figure 47b shows a schematic diagram of the
etching process consisting of three steps A–C. By defining the dissolution
rate of section i in step j as Rij (where i = 1, 2, 3; j=A, B, C), the relative
dissolution rates Rij /R3j is represented by Fig. 47c. As described in Sect. 2.3,
the relative dissolution rate of GeO2 -doped glass to pure silica glass depends
on the concentration of buffered HF mixed with a volume ratio of 40%NH4 F
sol. : 50%HF acid: H2 O = X:1:Y . We now use 1.7:1:1, 10:1:1, and 1.7:1:5, in
50 S. Mononobe
Fig. 47. (a) Cross-sectional profile of the index difference ∆n. r1 , radius of the
GeO2 -doped silica core; r2 , radius of the pure silica clad; r3 , radius of a fluorine-
doped silica support; (b) Etching method for fabricating application-oriented
probes; (c) Cross-sectional profile of the relative dissolution rate
We now discuss the etching process using 1.7:1:1, 10:1:1, and 1.7:1:5 in
steps A–C, respectively.
In step A, the fiber is tapered to an angle of θA represented by
If the fiber diameter is equal to [2r2 ] after step A, the etching time TA is
represented by
Near-Field Optical Fiber Probes and the Imaging Applications 51
The tapered fiber will have an apex diameter smaller than 2r1 if TA ≥τA ,
where τA is the time required to make an apex diameter of 2r1 and is
r2 − r1 R2A + R3A
τA = . (16)
R3A R3A − R2A
Thus, it is straightforward to find that the radius r2 of section 2 must be
smaller than the critical radius r2P expressed as
r3 + ξ·r1 R2A /R3A + 1
r2P = , where ξ = . (17)
1+ξ 1 − R2A /R3A
We obtain pencil-shaped probes with zero apex diameter and cone angles θA
and θB when the etching time TB is larger than τB , expressed as
r1 R1B + R3B
τB = . (19)
R1B R3B − R1B
Furthermore, to obtain a triple-tapered probe, we perform step C, when the
largest cone angle θC2 is given by
θC2 R1C
sin = (where R1C > R2C ) . (20)
2 R3C
The cone angle θC1 is increased from θB by increasing the etching time TC ,
and is equal to θC2 at TC > τC , where τC is given by
r1 R1C + R3C
τC = . (21)
R1C R3C − R1C
Therefore, we can obtain a triple-tapered probe and a pencil-shaped probe
with a cone angle θB = θC2 at 0 < TC < τC and TC > τC , respectively.
To obtain a cone angle θA as small as 62◦ , the index difference of section 3
is estimated from (15) and (19) to be as high as 0.7%. To obtain a cutoff
wavelength of around 400 nm, section 1 is tailored with an index difference
of 1.2% and a core radius of r1 = 0.65 µm. Then we obtain an estimated
value of θA = 17◦ from (13) and (19). Furthermore, when the clad radius r3
is a standard value of 62.5 µm, we obtain a critical radius of r2P = 23 µm
from (14) and (18). We make the radius r2 have a value of 13.5 µm, which is
smaller than the critical radius.
To realize the fiber as shown in Fig. 47a, we fabricated a preform glass
rod by vapor-phase axial deposition and drew the fiber. To suppress the
52 S. Mononobe
Fig. 48. (a–c) SEM micrographs of a pencil-shaped probe, its magnified tapered
core, and its magnified apex region. θA = 62◦ ; θB = 17◦ ; d < 10 nm; (d) SEM
micrograph of a triple-tapered probe. θA = 62◦ ; θC1 = 50◦ ; θC2 = 85◦ . (e) SEM
micrograph of a pencil-shaped probe. θA = 62◦ ; θB = 85◦
Fig. 49. (a) Dependencies of the taper angle θB on the volume ratio XB of NH4 F
aqueous solution in the etching solution; (b) Dependency of the taper angle θC2 on
the volume ratio of H2 O in the etching solution
magnified apex region, respectively. Here, the apex diameter is less than
10 nm. The cone angles are θA = 62◦ and θB = 17◦ , which agree with the
estimated values from (15) and (19).
Furthermore, by etching the pencil-shaped probe in a 1.7:1:5 solution
for TC = 2 min, we obtained a triple-tapered probe with a high-resolution
capability and high throughput. Figure 48d shows a SEM micrograph of
the triple-tapered probe. The probe has three cone angles of θC1 = 50◦ ,
θC2 = 85◦ , and θA = 62◦ and an apex diameter of less than 10 nm. The two
base diameters of tapers with θC1 = 50◦ and θC2 = 85◦ are dB1 = 250 nm
and dB2 = 1.28 µm, respectively. At TC = 2.75 min, we obtained another
type of pencil-shaped probe with a large angle of 85◦ near the apex region as
shown in Fig. 48e. Figures 49a and b show the dependencies of θB and θC2 ,
respectively, on the volume ratios of BHF. Here, 1.7:1:1, XB :1:1, and 1.7:1:YC
are used in steps A–C. The cone angles θB and θC2 are controlled by varying
XB and YC , respectively.
References
1. E. Betzig, P.L. Finn, J. S. Weiner: Appl. Phys. Lett. 60, 2484 (1992)
2. R. Toledo-Crow, P.C. Yang, Y. Chen, Vaez-Iravani: Appl. Phys. Lett. 60, 2957
(1992)
3. D.W. Pohl, W. Denk, M. Lanz: Appl. Phys. Lett. 44, 651 (1984)
4. E. Betzig, M. Isaacson, A. Lewis: Appl. Phys. Lett. 51, 2088 (1987)
5. D.R. Turner: US Patent, 4,469,554 (1983)
6. D. Courjon, K. Sarayeddine, M. Spajer: Opt. Commun. 71 23 (1989)
7. R.C. Reddick, R.J. Warmack, T L. Ferrel: Phys. Rev. B 39 767 (1989)
8. K.M. Takahashi: J. Colloid Interface Sci. 134, 181 (1990)
54 S. Mononobe
1 A General Method
for Preparing Surface-Bound Metal Particles
1.1 Surface-Bound Metal Particles
The usual method for stabilizing colloidal particles is to chemically tailor their
surface in such a way that the repulsive force among the spheres prevents
them from sticking together. By gradually decreasing the repulsive force,
for example, through an increase in the salt concentration, in the presence
of a high-energy surface, spheres begin to adsorb on the surface, forming a
monolayer. This is the basis of our methodology [10]. This simple scheme can
be expanded along a number of different themes.
We have purchased monodisperse polystyrene spheres, from Polysciences,
Warrington, PA, USA (Polystyrene bead Select Certified Size Standard Mi-
crosphere) or Dyno Particles, Norway (Dynospheres) and used them as re-
ceived. For most experiments we used various polystyrene (PS) labware for
life-science experiments as substrates. These were purchased from Nunc,
Nalgene etc. and used without washing. Prior to forming a monolayer of
polystyrene spheres, the PS substrate was coated with 20 nm gold in a vac-
uum evaporator (Type EBV-6 DA, ULVAC, Kanagawa, Japan). Gold, purity
of 99.95%, was purchased from Furuya Metal, Tokyo, Japan and was evapo-
rated from a tungsten boat at a base pressure of 2 × 10−6 Torr at an evapora-
tion rate of 0.1 nm s−1 ; PS substrates were held at least 30 cm away from the
tungsten boat in order to prevent melting of the surface. PS was chosen as the
substrate material because it is optically clear and flat PS pieces are readily
60 H. Takei and M. Himmelhaus
Sample Characterization
We list below physical parameters that can be readily and reproducibly con-
trolled.
Sphere Diameter
Adsorption Density
While the dense coverage of the surface with a monolayer of spheres is over
within ten minutes, by adjusting the incubation time one can control the
adsorption density of surface-bound spheres; Fig. 4 shows SEM pictures of
samples prepared with reduced incubation times. As the adsorption density
is reduced, spheres tend to form independent clusters.
Fig. 4. SEM micrograph of adsorbed spheres (sphere diameter 209 nm) at reduced
densities. Typical density at saturation is approximately 15.3 µm−2 , but by reducing
the sphere incubation period, the density can be reduced. (a) Density of 12.2 µm−2 .
(b) Density of 10.3 µm−2
Deposition Thickness
One can readily control the deposition thickness of the metal layers. There
are two layers involved, the one deposited directly on the substrate, and the
64 H. Takei and M. Himmelhaus
other layer deposited on top of the sphere layer. Varying the thickness of each
layer affects the optical property of the sample, but it is the layer deposited
on top of the sphere monolayer that has the most profound effect.
Gap Control
Fig. 6. SEM micrograph of spheres etched in oxygen plasma. The sphere size is
209 nm. (a) Etching duration of 60 s. (b) Etching duration of 50 s
Heat Treatment
To a lesser degree, we can control the shape with temperature. By raising the
temperature slightly above the onset of polymer deformation, it is possible
to deform the particle. Figure 7 is an SEM picture of a sample that has been
heated up to 120◦ C for a brief period. The sphere has become deformed and
now shows a greater contact area to the substrate. Interestingly, the presence
of the gold cap has the effect of retaining the original shape of the sphere. In
contrast, when surface-adsorbed PS spheres without gold are heated under
the same condition, spheres simply turn into shapeless globs, as shown in
Fig. 8.
66 H. Takei and M. Himmelhaus
Fig. 8. SEM micrograph showing the effect of heating bare polystyrene spheres
without gold on top. As with Fig. 7, the temperature was raised temporarily up to
120◦ C; in contrast to Fig. 7 lack of gold cap leads to complete merging of adjacent
particles. Image was taken at an angle of 45◦ from the surface normal
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 67
Patterned Particles
We have found that we can control the adsorption process of spheres with
patterns of alkanethiol molecules [11]. We prepared a gold-coated substrate
consisting of regions of different alkanethiol molecules. By controlling the
functionality of the alkanethiol molecule, it was possible to form patterned
layers of spheres. To gain control on the adsorption process a distinct interac-
tion between the spheres in suspension and selected areas of the surface needs
to be introduced. While it is possible to utilize physical phenomena, such as
changes in polarity, surface energy, or charge distribution in different regions
of the surface, we have obtained the best results for chemically patterned
surfaces in connection with the use of 1-ethyl-3-(3-(dimethylamino)propyl)
carbodiimide (EDC). Here, the role of EDC in the adsorption process of the
spheres is threefold. Due to its amino group, which is protonated in a neu-
tral pH environment, the molecule is positively charged and thus behaves as a
positively charged surfactant: it adsorbs on the negatively charged sphere sur-
face, thereby changing the total charge of the spheres, i.e., their ζ-potential,
and thus reducing interparticle repulsion. At the same time, the spheres are
converted into carbodiimide functionalized particles and thus provide poten-
tial reaction sites for chemical bonding. Furthermore, dissolved excess EDC
molecules act as macromolecules in the colloidal suspension and force aggre-
gation due to depletion forces. We have found that we have to adjust the EDC
concentration in suspension in such a way, that the body of EDC molecules
adsorbs on the spheres and the remaining dissolved molecules are too few
to exert depletion forces beyond kT . In this regime, which corresponds to
0.1–0.6 mM EDC (for suspensions with 2.6 vol.% spheres and particle sizes
of some 100 nm) adsorption turned out to be sensitive to the chemistry of
the surface. Note that this concentration is significantly lower than the value
for preparing uniform, non-patterned samples. At higher EDC concentration,
this sensitivity is screened by physisorption of the spheres due to depletion
forces. We have proved this concept by replacing EDC with a poly(ethylene
glycol) (PEG) molecule of similar size and molecular weight. PEG is non-
adsorbing and chemically inert, so that it solely exerts depletion forces on
the suspended spheres. At high concentration (above 0.6 mM), where the
depletion forces exceed the Brownian motion of the suspended particles, we
obtained similar results to those found for EDC at similar concentration.
At low concentration, sphere adsorption failed independent of the chemical
functionalities present at the surface. However, by using EDC and proper
surface functionalization, e.g., by means of carboxyl-terminated alkanethiols,
the carbodiimide functionalized spheres adsorb at low concentration due to
chemical bonding. We have combined this selected chemical reactivity with
the technique of microcontact printing to form patterned sphere layers. For
that, a nonfunctionalized thiol (e.g., octadecanethiol; C18 ) is printed onto a
native gold surface by means of a microstructured elastomer stamp (Fig. 9).
In this way, the C18 adsorbs only on the gold areas in contact with the stamp.
68 H. Takei and M. Himmelhaus
e la s to m e r s ta m p
C 1 8
=
A u
+ & P S
T G s o lu tio n C 1 8 T G C 1 8
> A u A u
P S P S
s p h e r e s u s p e n s io n
T G C 1 8
T G C 1 8
? A u A u
P S P S
In a second step, the gold substrate is immersed into a solution containing the
functionalized alkanethiol molecule (e.g., carboxyl-terminated such as thio-
glycolate; TG), filling the remaining native areas on the surface. In this way,
a chemical pattern of functionalized and nonfunctionalized areas has formed.
Exposure of this pattern to an EDC-activated sphere suspension at low EDC
concentration yields the formation of patterned sphere layers. Figure 10 is a
SEM picture of a patterned sample.
Fig. 10. SEM micrograph of patterned spheres. Spheres are selectively adsorbed on
regions covered by thioglycolate separated by regions covered by octadecanethiol
Fig. 11. Photograph of gold particles in the shape of letters. Twenty nanometers
of gold were evaporated onto a polystyrene substrate, followed by formation of a
uniform monolayer of spheres (110 nm in diameter). By using a stencil in the shape
of letters, an additional 20 nm of gold were evaporated within well-defined regions.
This is to illustrate visually the optical contrast formed by the gold particle
70 H. Takei and M. Himmelhaus
Fig. 12. Photograph of gold particles formed with PS spheres of different diam-
eters. Spheres are adsorbed within circles defined by the rim while the number
superimposed indicates the sphere diameter in nanometers. The spheres on the
left-hand half of the sample are left bare whereas twenty nanometers of gold were
evaporated on top of the spheres on the right-hand half
The deposition thickness has been found to have a completely different effect
depending on whether it is the top or bottom layer that is affected. While the
bottom layer must be present in order to enhance the absorption spectrum,
its exact thickness matters relatively little. On the other hand, the amount of
deposited gold on top of the PS sphere has a great bearing on the absorption
spectrum, as shown in Fig. 15 for 110-nm spheres. For other sphere diameters
as well, initially a broad peak begins to form at longer wavelengths, and as the
deposition thickness is increased, the peak height increases and shifts toward
shorter wavelengths. Eventually, further increases fail to produce additional
shifts. The saturation limit begins to set in at around 20 nm.
Fig. 15. Effect of the deposition thickness on the absorption spectrum. The sphere
diameter is 110 nm for all the samples
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 73
Sphere Density
Figure 16 shows spectra of samples with different sphere densities. The den-
sity corresponds to 10.3 µm−2 , 12.2 µm−2 , and 15.3 µm−2 from sparse,
moderate and dense samples, respectively, for 209-nm spheres in (a), and
Fig. 16. Effect of the adsorption density on the absorption spectrum. (a) 209-nm
spheres: the density is 10.3 µm−2 , 12.2 µm−2 , and 15.3 µm−2 from sparse, moderate
and dense samples, respectively. (b) 110-nm spheres: the density is 39.8 µm−2 ,
41.9 µm−2 , and 61.9 µm−2 from sparse, moderate and dense samples, respectively
74 H. Takei and M. Himmelhaus
39.8 µm−2 , 41.9 µm−2 , and 61.9 µm−2 from sparse, moderate and dense
samples, respectively, for 110-nm spheres in (b).
Mixed Samples
Fig. 17. Effect of mixing two sphere sizes at different mixture ratios. (a) 60.9 µm−2
for 110 nm. (b) 44.7 µm−2 for 110 nm and 5 µm−2 for 209 nm. (c) 20.4 µm−2 for
110 nm and 10.1 µm−2 for 209 nm. (d) 15.5 µm−2 for 209 nm
Material Dependence
We have concentrated our effort with gold particles simply due to its stabil-
ity in air, but on occasion we have also experimented with other metals such
as silver, copper and platinum. The spectra for silver samples are shown in
Fig. 18. Twenty nanometers of silver were evaporated on spheres with diam-
eters 55, 110, 152, 209 nm. They were quite susceptible to the atmosphere,
as has been reported before in different systems [22].
Gap Dependence
The ability to control the gap distance at will between adjacent particles is
one of the unique features of the current particle system. It is precisely this
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 75
Fig. 18. Absorption spectra of silver particles. Twenty nanometers of silver were
deposited on spheres of various diameters adsorbed on a 20-nm thick gold film
Fig. 19. Absorption spectra of samples prepared from spheres treated to oxygen
plasma etching prior to deposition of the second gold layer. It can be seen that with
increasingly longer etching periods, the absorption maximum shifts toward shorter
wavelengths
The present chapter deals mostly with a new methodology. We have presented
some examples of what can be prepared and characterized them only in terms
of absorption spectra and SEM micrographs. Most of the physical parameters
that have been mentioned are well known to affect absorption spectra of
noble metal particles; what is new here is the fact that these parameters can
be, in any combination, modified reproducibly, and resulting samples are of
sufficient size that makes macroscopic characterization an easy task. This is
important not only for practical application of nano-optics based on metal
particles, but also for rapid confirmation of theories with experimental data.
The size dependence is something that has always been part of the dis-
cussion over noble metal particles. While analysis is the most stringent for
an isolated spherical solid particle of limited size, it has been extended to
particles of different shapes such as ellipsoid, hemisphere, rod and cube [2].
In general, the fundamental mode of collective free-electron oscillation dom-
inates in the quasistatic limit where the size is significantly smaller than the
wavelength; modeling with free electrons applies well only for alkali metals,
while noble metals, particularly gold, show pronounced deviations. For shapes
other than a sphere such as the cap-shaped particle in this chapter, one would
expect more than one fundamental mode even in the quasistatic limit. As it is,
the sizes of particles prepared are beyond this limit, so that if the cap-shaped
particle were suspended in isolation, one would expect to observe more than
one peak. One might expect that if particles are oriented under an external
force, it would be possible to selectively observe these peaks, but studies on
cubic and hemispherical particles have revealed that the presence of edges
and corners would preclude homogeneous polarization even when particles
are oriented. Such an experiment is currently planned.
As the deposition thickness is increased, the diameter stays more or less
constant as determined by the cross section of the polystyrene sphere, but
the shape undergoes changes as the cap becomes thicker. In contrast to hemi-
spherical particles and core-shell particles [2,20,23], cap-shaped particles have
rarely been studied as far as we are aware. The core-shell particles, particu-
larly with a dielectric core and a gold shell, have been receiving much interest
in recent years [24–26], as they exhibit absorption peaks in the near infrared
region, and the precise peak position can be controlled by the relative thick-
ness of the shell. It can be imagined that as the deposition thickness is in-
creased for the cap-shaped particle, the characteristics would shift from that
of a core-shell particle to a hemispherical particle, if observed in isolation.
Such an investigation is currently underway.
Difficulties of analyzing the actual samples prepared stem from the fact
that particles are adsorbed on a substrate at high densities. Thus it would
be necessary to take into account image interactions due to the substrate
and interactions among adjacent particles. According to earlier studies on
the effect of the substrate [23], multipole effects are enhanced and many
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 77
secondary modes appear. In addition, with the sample under discussion here,
the underlying gold film apparently constitutes part of an optical cavity, thus
enhancing the overall spectrum [10, 27].
The effect of interactions among adjacent particles is observed in more
than one way. Varying the sphere density is one way, and the shift in the
peak position as seen in Fig. 16 is one indication (the peak height is not
a straightforward indicator of the sphere density as reduction in the density
increases the area of pure gold surface that by reflection reduces the apparent
peak height). Mixed samples are another indication. Particle systems with a
certain size distribution are quite common, and they are characterized as
mixtures of independent components. The spectra shown in Fig. 17 show,
however, that spheres of two different sizes interact to give rise to a single
new absorption peak. It should be mentioned that control over the sphere
adsorption density has received attention recently [28].
Samples where a clear gap is formed among adjacent particles provide the
clearest evidence for interactions. In recent years, there has been a revival
of interest in coupled particles, particularly in connection with the enhanced
near field in the gap region. There have been a number of theoretical re-
ports [29, 30], but experimentally it has proven difficult to create samples
consisting of coupled particles with a well-controlled gap distance; it should
be mentioned that Garcia-Santamaria et al. [31] have made an attempt with
silica-coated gold nanoparticles. To our knowledge, we are reporting here
for the first time an experimental technique that allows preparation of such
samples. For now, we limit ourselves to the fact that the gap distance can be
controlled by varying the etching duration, and the extent of the blue shift
becomes greater with longer etching durations. There is no doubt that the
etching also reduces the thickness of PS spheres so that this needs to be taken
into account for future analysis. Moreover, to facilitate comparison with the-
oretical analysis, it would be desirable to prepare samples in such a way that
the axis connecting all coupled particles would be the same because the rela-
tive angle between the axis and the polarization direction of the probing light
plays an important role in the type of excited modes. One such sample would
consist of parallel lines of stringed particles; it should be possible to prepare
such samples by exploiting the patterning technique mentioned earlier.
It is quite exciting to speculate on and eventually verify experimentally
the strength of the near field in the gap region. Theoretical work suggests [30]
that the gap of coupled spheres is an ideal location to generate an intense
“hot” spot that can be utilized for surface-enhanced Raman spectroscopy [32–
34]. This might be verifiable by a combination of Raman spectroscopy and
atomic force microscopy in the fashion reported by [35]; they attempted to
correlate between the morphology of particle clusters and the strength of
observed Raman signal but could not draw a conclusion on the shape of
the ideal “hot” particle. Their conclusion, was however, that the intrinsic
Raman enhancement factors of “hot” particle were on the order of 1014 to
78 H. Takei and M. Himmelhaus
Experimental Setup
Fig. 20. Schematic diagram of the experimental setup. Gold particles were formed
at the bottom of a polystyrene well with a volume of 10 µl. Light from a tungsten
halogen lamp is led to the well bottom via a fiber bundle, and the reflected light
is collected for real-time analysis by a spectrometer. The sample is mechanically
agitated to assure immediate and uniform mixing
80 H. Takei and M. Himmelhaus
Surface Coating
There are many methods for adsorbing biomolecules onto the sensor surface.
(1) Direct physisorption
It is well known that many biomolecules bind irreversibly to a bare gold
surface. It is thus possible to inject a buffer solution containing a biomolecule.
The selection of the buffer solution is not crucial, but we often use phosphate
buffer saline (PBS) maintained at pH 7.4. Although the concentration of
1 mg/ml quickly saturates the surface in a matter of a few minutes, the
concentration of 100 µg/ml is enough with a slightly longer incubation period.
(2) Coupling via protein A or protein G
Direct physisorption described above is rather straightforward, but there is
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 81
some concern that the antibody may partially lose its ability to bind its
antigen through denaturation. This possibility is reduced if the antibody is
adsorbed indirectly via an intermediate layer. Protein A and protein G can
serve such a purpose; they both have a multiple number of binding sites for
antibodies (four for protein A and two for protein G) so that physisorbed
protein A or protein G is likely to retain its function. As in (1), protein A or
protein G at the concentration of 100 µg/ml can be injected into the sensor
well, and after some 10 minutes, the coverage is more or less complete. Then,
a solution containing an antibody is injected. The origin of the antibody,
whether from mouse, goat, human, etc., has some effect on the affinity of
protein A or protein G toward the antibody.
(3) Coupling via biotin–avidin (streptavidin) reaction
The great affinity between biotin and avidin (or similarly for streptavidin) is
well known, and most biomolecules can be readily modified by a biotin moi-
ety; thus physisorbed avidin can serve as the binding layer for a biotinylated
biomolecule.
(4) Activated carboxyl group
Formation of a covalent bond between an activated carboxyl group and an
amino group is widely used as a linking method; most biomolecules possess
an amino group so that a surface-bound carboxyl group can serve as a start-
ing point. Gold is a convenient surface because an alkanethiol molecule is
known to readily form a self-assembled monolayer (SAM) and alkanethiol
molecules with carboxyl functional groups are commercially available. An
ethanol solution containing 10-carboxy-1-decane thiol (10-carbo.), purchased
from Dojindo, is prepared at the concentration of 1 mM. Exposure of the
sensor well to the solution for a period of 10 minutes is enough to prime the
sensor surface with a carboxyl group. A freshly prepared mixture of EDC
and N-hydroxy succinimide (NHS), 0.1 M, is added for some 10 minutes and
immediately afterwards the biomolecule to be adsorbed is injected.
Antigen–Antibody Reaction 1
in the peak wavelength. Once the binding becomes saturated, the unbound
antigen is washed off with phosphate buffer. To remove the antigen, addition
of diluted hydrochloric acid, 10 mM, is enough to break up the antigen–
antibody reaction. This is reflected in a reduction in the peak wavelength.
Addition of diluted acid, however, leaves the antibody intact so that when
the antigen is added again, the surface-bound antibody can capture it. This
process can be repeated.
Antigen–Antibody Reaction 2
Fig. 23. Dependence of the sensorgram on the sample concentration. The covered
range is from 0.25 µg/ml to 32 µg/ml. At t = 200 s, anti-streptavidin antibody
is injected onto a sensor covered by streptavidin. At t = 800 s, phosphate buffer
saline (PBS) is used to rinse off unbound antibody. After t = 900 s, the sensor is
exposed to diluted hydrochloric acid to break the bond between streptavidin and
anti-streptavidin
Protein–DNA Interactions
While there are many types of biosensors [6, 8, 41–44], one of the most im-
portant parameters for a biosensor is sensitivity. Naturally it is roughly de-
84 H. Takei and M. Himmelhaus
Fig. 24. Obtaining kinetic data. Sensorgrams from Fig. 23 are fitted to a single
exponential function, A − B exp(−kt), and the binding constant, k, was obtained
for various concentrations. All the binding constants were plotted as a function of
the concentration; linear regression was used for fitting. The slope of the fitted data
is the association time constant, kass , and the y intercept is the dissociation time
constant, kdiss . The equilibrium dissociation constant, kD , is kdiss /kass . kD of 63 nM
is fairly typical for an antigen–antibody reaction. The quality of linear regression
also indicates the reliability of our sensor
The sensor is exposed to various fluids of different refractive indices. This was
achieved by preparing mixtures of water and glycerol with different mixture
ratios as shown in Fig. 28. With refractive indices of 1.33 for water and 1.47
for glycerol, the sensorgram shows that the change in the peak wavelength
for a unit change in the refractive index is 61 nm. In contrast, the litera-
ture figure for the traditional SPR, when measured in terms of wavelength,
is typically 3000 nm per unit change in the refractive index [42, 47]. This
86 H. Takei and M. Himmelhaus
Fig. 28. Sensorgram showing sensor response toward water/glycerol mixtures. The
volume concentration of glycerol is incrementally increased by 10%; the first mixture
contains 10% glycerol, and the last mixture added at t = 1900 s has 80% glycerol.
With refractive indices, n = 1.33 for water and 1.47 for glycerol, it is roughly
estimated that the shift is 61 nm per one refractive index unit, δn = 1
would suggest that the sensitivity of our sensor is significantly lower. Fig-
ure 29, however, shows an interesting comparison between the SPR sensor
and our sensor subjected to an identical protocol; the first half consists of a
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 87
Fig. 29. The following protocol was used to illustrate the reduced “bulk” effect of
our LSPR instrument with respect to the SPR instrument; at t = 100 s, 1 mg/ml
avidin was added to the sensor to initiate physisorption. After saturation, unbound
avidin was rinsed with PBS buffer at t = 400 s. 100 µg/ml of anti-avidin antibody
was injected to start the antigen–antibody reaction; the overall shift in the wave-
length is typical. The portion of the sensorgram after t = 1200 s is identical to the
initial portion of Fig. 28. (a) Response of an SPR sensor from Texas Instruments,
Spreeta, that was chosen as a representative SPR sensor. (b) Our instrument. Note
how the shifts induced by injection of glycerol mixtures are suppressed in (b)
88 H. Takei and M. Himmelhaus
Temperature Dependence
Figure 30 shows the response of the sensor when water at various temper-
atures is injected. Initially it is at room temperature (20◦ C), and the tem-
perature of water injected is increased incrementally from 1 to 63◦ C; each
injection results in a peak because the temperature returns to room temper-
ature rather rapidly. There are two noteworthy features. One is that even
though the refractive index is higher for water at lower temperature in this
temperature range, the shift in the peak wavelength is opposite; an increase in
the density results in a decrease in the peak wavelength. This suggests that a
change in temperature has an effect on the sensor through a mechanism other
than a simple change in the refractive index, possibly a structural change in
the nanostructure of the sensor itself. The second noteworthy feature is that
there is hysteresis; the peak wavelength at room temperature, after the se-
ries of temperature jumps, is permanently changed. Again, this might reflect
a structural change in the nanostructure of the sensor itself. While this is-
sue might call for attention in the future, the fact that temperature-induced
shifts are quite small and most if not all binding experiments will not subject
the sensor to such large jumps in temperature renders this a minor nuisance
rather than a problem. In contrast, a typical SPR instrument is characterized
by a jump of 0.01 angular degree (e.g., 100 resonance units in the industry
parlance) per 1 degree Celsius [49]. In terms of the wavelength of our instru-
ment, it would loosely translate to 0.5 nm/degree Celsius.
In summary, we have developed a novel biosensor based on optical proper-
ties of gold nanoparticles. There have been numerous reports in recent years
on similar sensors also based on localized surface plasmon resonance, LSPR,
of gold nanoparticles [6, 8]. Advantages cited for the LSPR sensor over the
SPR sensor are simplicity of the instrument, low costs, ability to give a result
that can be read off visually, etc. On the other hand, there have been few
reports on the actual performance. It has been observed that the dependence
of the peak wavelength on the refractive index of the surrounding material is
significantly lower, by more than one order of magnitude. We also make sim-
ilar observations, but we emphasize that as far as detection of biomolecules
is concerned, it has a sensitivity significantly higher than what the result
of the bulk dependence leads one to believe. The reduced “bulk” effect is a
tremendous advantage because molecular-binding events can be monitored
stably in the presence of fluctuations in temperature and salt concentrations
of the buffer solution.
One outstanding characteristic of our LSPR sensor is the extraordinarily
high optical density and the narrow width of the absorption peak. Typically
the optical density of surface-adsorbed gold nanoparticles is less than 0.5,
but the figure for our sensor is more than four times greater. This facilitates
accurate determination of the peak wavelength, better than 0.02 nm. As the
shift induced upon formation of a SAM, corresponding to a mass loading of
1 ng/mm2 , is typically 5.5 nm, under the assumption that linearity holds,
this leads to the detection limit of under 5 pg/mm2 . This compares well with
the reported practical limit of 10 pg/mm2 for the SPR sensor.
3 Curious Observation
While in most cases, the response of our sensor is quite comparable to that
of the SPR sensor with the generally observed difference mentioned earlier,
some protocols have led to curious observations. One such case is shown in
90 H. Takei and M. Himmelhaus
Fig. 31. Binding of amino dextran (MW 10 kilodalton) onto the sensor coated with
10-carbo. Ostensibly, the negative charge of the thiol molecule attracts electrostat-
ically positively charged amino dextran suspended in water. Remarkably, the shift
is well over 120 nm. If the refractive-index dependence in Fig. 28 is any indication,
this would reflect change of more than two refractive index units. Obviously, this
cannot be so, thus inviting further exploration for an explanation
Fig. 31. Here, The sensor was first exposed to an ethanol solution containing
10-carbo at the concentration of 1 mM. Then, excess 10-carbo was rinsed with
pure water, resulting in reduction in the peak wavelength. Then, a solution
containing 10 kilodalton amino dextran was injected. Opposite charges on the
surface-bound carboxyl group and amino dextran are expected to lead to for-
mation of an amino dextran layer. Indeed this was what we observed, but the
shift was phenomenal at well over 120 nm. Although the exact shift depends
on the molecular weight of the amino dextran used and its concentration,
such phenomenal shifts are routinely observed with this particular protocol.
When we recall that a shift due to a SAM formation was well under 10 nm,
and a shift caused by a change in the refractive index was some 61 nm per
unit change, a shift of more than 100 nm is surprising. While we suspected
a permanent change in the sensor structure, observation of the sensor with
a scanning electron microscope did not reveal any noticeable change, not to
mention removal of the spheres.
4 Conclusion
We have shown a new method for forming surface-adsorbed metal parti-
cles; The beauty of this method is its simplicity and practicality. While this
A Novel Method for Forming Uniform Surface-Adsorbed Metal Particles 91
Acknowledgements
This work has been made possible by a grant MF-16 from the Ministry
of Health and Labor. We express much appreciation to Dr. T. Fujimura,
Dr. J. Pipper, Dr. Z. Cao, Dr. T. Okamoto, Dr. H. Sugiyama, Dr. K. Ni-
nomiya, Dr. T. Iwayanagi, and Dr. H. Taguchi in the course of this project.
Without their help, the project would not have been possible. We also thank
the Ministry of Health and Labor, and Dr. Inoue of the National Institute
of Health Science of Japan for their partial support of this project. Professor
M. Furuya and Dr. H. Kambara have been constant sources of encourage-
ment.
References
1. T. Kume, N. Nakagawa, S. Hayashi, K. Yamamoto: Solid State Commun. 93,
171 (1995)
2. Optical Properties of Metal Clusters by U. Kreibig, M. Vollmer (Springer, Berlin
1995)
3. Nanoparticles and Nanostructured Films ed. by J.H. Fendler (Wiley-VCH,
Weimheim 1998)
4. R. Micheletto, H. Fukuda, M. Ohtsu: Langmuir 11, 3333 (1995)
5. G. Kalyuzhny, M.A. Schneeweiss, A. Shanzer, A. Vaskevich, I. Rubinstein:
J. Am. Chem. Soc. 123, 3177 (2001)
6. F. Meriaudeau, T. Downey, A. Wig, A. Passian, M. Buncick, T.L. Ferrell: Sens.
Actuators B 54, 106 (1999)
7. M. Brust, D. Bethell, C.J. Kiely, D.J. Schiffrin: Langmuir 14, 5425 (1998)
8. T. Okamoto, I. Yamaguchi, T. Kobayashi: Opt. Lett. 25, 372 (2000)
9. G.A. Niklasson, H.G. Craighead: Thin Solid Films 125, 165 (1985)
10. H. Takei: J. Vac. Sci. Technol. B 17, 1906 (1999)
11. M. Himmelhaus, H. Takei: Phys. Chem. Chem. Phys. 4, 496 (2002)
12. A.S. Dimitrov, C.D. Dushkin, H. Yoshimura, K. Nagayama: Langmuir 10, 432
(1994)
13. T. Fujimura, T. Itoh, A. Imada, R. Shimada, T. Koda, N. Chiba, H. Mura-
matsu, H. Miyazaki, K. Ohtaka: J. Lumino. 954, 87 (2000)
14. Y. Masuda, M. Itoh, T. Yonezawa, K. Koumoto: Langmuir 18, 4155 (2002)
15. H. Takei, N. Shimizu: Langmuir 13, 1865 (1997)
16. C. Haginoya, M. Ishibashi, K. Koike: Appl. Phys. Lett. 71, 2934 (1997)
92 H. Takei and M. Himmelhaus
17. J.C. Riboh, A.J. Haes, A.D. McFarland, C.R. Yonzon, R.P. Van Duyne:
J. Phys. Chem. B 107, 1772 (2003)
18. M.D. Malinsky, K.L. Kelly, G.C. Schatz, R.P. Van Duyne: J. Am. Chem. Soc.
123, 1471 (2001)
19. U.C. Fischer, D.W. Pohl: Phys. Rev. Lett. 62, 458 (1989)
20. P. Barnickel, A. Wokaun: Mol. Phys. 67, 1353 (1989)
21. G.L. Hornyak, C.R. Martin: Thin Solid Films 303, 84 (1997)
22. Y.W. Cao, R. Jin, C.A. Mirkin: J. Am. Chem. Soc. 123, 7961 (2001)
23. R.J. Warmack, S.L. Humphrey: Phys. Rev. B 34, 2246 (1986)
24. H.S. Zhou, I. Homma, H. Komiyama, J.W. Haus: Phys. Rev. B 50, 12052 (1994)
25. R.D. Averitt, D. Sarkar, N.J. Halas: Phys. Rev. Lett. 78, 4217 (1997)
26. T. Pham, J.B. Jackson, N.J. Halas, T.R. Lee: Langmuir 18, 4915 (2002)
27. A. Leitner, Z. Zhao, H. Brunner, F.R. Aussenegg, A. Wokaun: Appl. Opt. 32,
102 (1993)
28. R.R. Bhat, D.A. Fischer, J. Genzer: Langmuir 18, 5640 (2002)
29. J.P. Kottmann, O.J.F. Martin: Opt. Lett. 26, 1096 (2001)
30. H. Xu, M. Käll: Phys. Rev. Lett. 89, 246802 (2002)
31. F. Garcia-Santamaria, V. Salgeuirino-Maerira, C. Lopez, L.M. Liz-Marzan:
Langmuir 18, 4519 (2002)
32. P.F. Liao, M.B. Stern: Opt. Lett. 7, 483 (1982)
33. R.G. Freeman, K.C. Grabar, K.J. Allison, R.M. Bright, J.A. Davis,
A.P. Guthrie, M.B. Hommer, M.A. Jackson, P.C. Smith, D.G. Walter,
M.J. Natan: Science 267 1629 (1995)
34. K. Kneipp, Y. Wang, H. Kneipp, L.T. Perelman, I. Itzkan, R.R. Dasari,
M.S. Feld: Phys. Rev. Lett. 78, 1667 (1997)
35. S. Nie, S.R. Emory: Science 275, 1102 (1997)
36. A. Dogariu, T. Thio, L.J. Wang, T.W. Ebbesen, H.J. Lezec: Opt. Lett. 26, 450
(2001)
37. R. Elghanian, J.J. Storhoff, R.C. Mucic, R.L. Letsinger, C.A. Mirkin: Science
277, 1078 (1997)
38. P. Bao, A.G. Frutos, C. Greef, J. Lahiri, U. Muller, T.C. Peterson, L. Warden,
X. Xie: Anal. Chem. 74, 1792 (2002)
39. H. Takei: In: SPIE Microfluidics Devices and Systems (Santa Clara, California,
1998) 278
40. M. Himmelhaus, H. Takei: Sens. Actuators B 63, 24 (2000)
41. C.R. Musil, D. Jeggle, H.W. Lehmann, L. Scandella, J. Gobrecht, M. Döbeli:
J. Vac. Sci. Technol. B 13, 2781 (1995)
42. J. Homola, S.S. Yee, G. Gauglitz: Sens. Acutators B 54, 3 (1999)
43. N.-P. Huang, J. Voros, S.M. De Paul, M. Textor, N.D. Spencer: Langmuir 18,
220 (2002)
44. A.J. Thiel, A.G. Frutos, C.E. Jordan, R.M. Corn, L.M. Smith: J. Vac. Sci. Tech-
nol. B 13, 2781 (1995)
45. A. Ulman: Chem. Rev. 96, 1533 (1996)
46. S. Maltis, I. Rubinstein: Langmuir 14, 1116 (1998)
47. L.S. Jung, C.T. Campbell, T.M. Chinowsky, M.N. Mar, S.S. Yee: Langmuir 14,
5636 (1998)
48. H. Takei, M. Himmelhaus, T. Okamoto: Opt. Lett. 27, 342 (2002)
49. Biomolecular Sensors ed. by E. Gizeli, C.R. Lowe (Taylor & Francis, London
2002)
Near-Field Optical-Head Technology for
High-Density, Near-Field Optical Recording
T. Matsumoto
1 Introduction
density cannot be increased beyond this limit. In the case of magnetic record-
ing, the recorded information is expected to decay spontaneously because of
thermal fluctuation in the magnetic domains when the recording density ap-
proaches 1 Tb/in2 . A solution is to increase the coercivity of the magnetic
recording medium. However, conventional writing heads cannot generate a
magnetic field large enough to write data on such recording medium because
of a fundamental limit on the saturation moment of the writing head material.
As a method to overcome these limits, near-field optical recording is now
drawing much attention. In this recording, the data bits are written and
read by using an optical near field generated near a nanometer-scale object.
Because the size of the optical spot is not limited by the diffraction of light,
and can be reduced to less than a few tens of nanometers [1], this method is
expected to enable a recording density of over 1 Tb/in2 .
In a near-field optical recording device, the key component is a near-field
optical head that generates the optical near field. In this chapter, after briefly
reviewing near-field optical recording, we describe technical issues regarding
the near-field optical head and introduce some solutions to resolve them.
Among these issues, improved efficiency in generating the optical near field
is the most important. We describe a highly efficient near-field optical head
that uses a metallic plate; the intensity of the optical near field is enhanced
by a plasmon generated in the metallic plate.
The recording density of conventional optical recording drives has been in-
creased by reducing the spot size of the focused laser. However, the minimum
spot size, D, is governed by the law of diffraction and is approximately
D = λ/N A , (1)
D = λ/nN A . (2)
In addition, in the case of the Weierstrass-sphere lens, the incident rays are
refracted at the sphere’s surface. This corresponds to an increased N A, and
the spot size is further reduced to
D = λ/n2 N A . (3)
This small optical spot enables increased recording density. For example, a
recording density of about 70 Gb/in2 has been realized using an SIL made of
Bi4 Ge3 O12 with a refractive index of 2.23 [4]. This type of recording method is
sometimes classified as near-field optical recording because it uses an evanes-
cent wave generated by a ray whose incident angle is larger than the critical
angle. However, the spot size is still limited by the diffraction limit of light,
and no dramatic increase of recording density is expected.
In near-field optical recording, the data bits are written and read by using
the optical near field. When light is introduced into a nanometer-scale object,
such as a subwavelength aperture or scatterer, a localized electromagnetic
field – the optical near field – is generated near the object. The distribution
of the optical near field is determined by the object’s shape and size, and
does not depend on the light wavelength. For example, when the optical
near field is generated by the subwavelength aperture, the width of the near-
field distribution is as small as the aperture diameter [1]. Therefore, a high
96 T. Matsumoto
recording density beyond the diffraction limit can be realized by writing and
reading data bits using the optical near field.
Betzig et al. [5] have demonstrated near-field optical recording using a
magneto-optical recording medium. They generated the optical near field us-
ing a fiber probe, which is widely used in scanning near-field optical mi-
croscopy (SNOM). In the writing process, a Pt/Co multilayer recording
medium was heated to the Curie temperature by the optical near field to
reverse the magnetization of the recording medium. In the reading process,
transmitted light from the recording medium was collected by the fiber probe
and the rotation of the light polarization caused by the Faraday effect was
detected. They have written and read recording marks with a diameter of
60 nm, which corresponds to a recording density of 170 Gb/in2 .
Near-field optical recording using a phase-change recording medium or a
photochromic recording medium has been reported by, respectively, Hosaka
et al. [6] and Jiang et al. [7]. In the phase-change recording, the recording
medium is heated by light to change its phase from crystalline to amorphous,
or vice versa. Hosaka et al. used GeSbTe recording medium, and wrote and
read recording marks with a diameter of 60 nm using a fiber probe. In the
photochromic recording, the recording marks are written through a photo-
chemical reaction. Jiang et al. wrote and read recording marks with a diam-
eter of 130 nm on a Langmuir–Blodgett film of photochromic material.
3 Technical Issues
Regarding the Near-Field Optical Head
In the primary experiments done on near-field optical recording, a fiber probe
was used as the near-field optical head. However, this probe was designed
Near-Field Optical-Head Technology 97
mainly for near-field optical microscopy and is not suitable for practical
recording devices. To realize a practical near-field optical head, the following
requirements must be met.
1. Precise control of the spacing between the head and the recording medium
2. Integration of peripheral components into the near-field optical head
3. Highly efficient generation of the optical near field
The optical near field is localized near the aperture or the scatterer and its
intensity decreases exponentially as the distance from the surface of the head
increases. Because the decay length is about 10–20 nm, the spacing between
the head and the recording medium must be maintained within a range of
10–20 nm.
With the fiber probe, the spacing is controlled by using the shear force
acting between the fiber probe and the recording medium [10, 11]. In this
method, the fiber probe is dithered parallel to the surface of the recording
medium. Because the dithering amplitude is changed by the shear force, a
feedback loop can be constructed by detecting the change in the dithering
amplitude. However, in a practical recording device, the recording disk rotates
at a very high speed. For example, when the recording density is 1 Tb/in2 ,
and the data transmission rate is 1 Gbps, the linear velocity of the recording
medium is about 30 m/s. If the fiber probe is placed near recording medium
rotating at such speed, the fiber probe soon crashes because the response
based on shear-force feedback is too slow.
To solve this problem, planar near-field optical heads have been proposed.
In this type of head, the structure for generating the optical near field is
mounted on the slider used for the hard disk drive. In the hard disk drive,
the spacing between the head and the recording medium is maintained at a
constant value by air pressure acting between the slider and the recording
medium. (The air pressure is determined by the design of pads formed on
the slider surface.) The head–medium spacing of hard disk drives is currently
near 10 nm, and is expected to decrease to less than 10 nm within a few years.
By applying this technology to the near-field optical head, we can precisely
control and maintain the spacing between the near-field optical head and the
recording medium.
Figure 3 shows an example of a planar near-field optical head.
Lee et al. [12] have proposed a silicon planar near-field optical head that
has an array of square apertures (Fig. 3a). The 60 × 60 nm apertures were
fabricated by an isotropic etching process with high reproducibility, and a line
pattern with linewidths of 250 nm was read out with this near-field optical
head.
98 T. Matsumoto
Fig. 3. Examples of planar near-field optical heads: (a) Head with an array of
apertures formed in a Si substrate, (b) Head with an aperture formed on a glass
protrusion, (c) Head with an array of apertures formed on Si protrusions
Isshiki et al. [13] have proposed a planar near-field optical head with an
aperture formed on a pyramidal protrusion (Fig. 3b). They fabricated the
pyramidal protrusion by cutting the glass substrate with a focused ion beam,
and fabricated the aperture by coating the protrusion with metal and slicing
the apex with a focused ion beam. They read out a 250-nm line pattern on
a disk rotating at a linear velocity of 0.38 m/s. The data transmission rate
and carrier-to-noise ratio were, respectively, 1.5 Mbps and 37 dB (bandwidth:
10 kHz).
Yatsui et al. [14] have also proposed a planar near-field optical head with
an array of apertures formed on silicon pyramidal protrusions (Fig. 3c). They
fabricated the array of silicon pyramidal protrusions through an isotropic
etching process, and formed the apertures at the apices of the protrusions.
They used a phase-change recording medium, and succeeded in writing and
reading data bits using the apertures. The recording medium was rotated at
a linear velocity of 0.43 m/s, and the minimum mark length and the data
transmission rate were, respectively, 110 nm and 2.0 MHz. Note that parallel
writing and reading using the aperture array is expected to enable a high
data transmission rate.
When a subwavelength aperture is used to generate the optical near field, the
efficiency (defined as the ratio between the power of the optical near field and
that of the incident light) is quite low. For example, the efficiency is 10−4 –
10−5 when the aperture diameter is 60 nm [5], and it decreases further as the
aperture diameter becomes smaller. This is because a cutoff diameter exists
for mode propagation in the waveguide. When the diameter of the tapered
waveguide becomes smaller than the cutoff diameter, the propagating light
suffers a huge loss and its intensity decays exponentially.
This low efficiency is a serious hindrance to realizing practical near-field
optical recording devices. For example, in the case of hybrid recording, the
temperature rise at the heated point in the recording medium is 40◦ C when
a recording medium moving at a linear velocity of 15 m/s is irradiated by
light with a spot size of 50 nm, peak power of 1 mW, and a pulsewidth of
2 ns [18]. Given these values and the output power of available laser diodes,
the required efficiency for writing is 1–10% if we assume that the temperature
rise required for writing is 200◦ C.
To raise efficiency, several methods have been proposed, such as
1. Optimization of the waveguide shape
2. Use of a surface plasmon
3. Use of a metallic waveguide
4. Use of a planar antenna
Low efficiency has also been a problem for SNOM when observing a weak
optical signal, such as photoluminescence, or a nonlinear signal. Although
most of these methods have been developed for SNOM applications, they
100 T. Matsumoto
between the second and third tapers excites the HE-plasmon mode, which
has a smaller cutoff diameter. The efficiency of the double-tapered probe is
10 times higher than that of a conventional probe, and the efficiency of the
triple-tapered probe is 1000 times higher than that of a conventional probe.
Fischer et al. [21,22] have proposed a tetrahedral probe that uses the surface
plasmon generated on a surface of a metallic film. The surface plasmon is a
collective oscillation of electrons generated on a surface of a thin metallic film
by irradiating the film with light [23]. When the surface plasmon is excited,
the electric field near the metallic field is enhanced. In the tetrahedral probe,
the surface plasmon is excited on the metallic film on the tetrahedral glass
chip, and propagates toward the apex. Thus, the electric field at the apex is
enhanced.
Thio et al. [24] have proposed a structure that consists of a subwavelength
aperture created in a metallic film and a set of concentric circular grooves
surrounding the aperture. When light is introduced into this structure, the
surface plasmon is excited on the metallic film by the circular grooves, and
the optical near-field intensity at the aperture is enhanced because the sur-
face plasmon assists the light transmitting through the aperture. In their
experiment, the transmission of the aperture with the grooves was 50 times
that of an aperture without grooves.
Fischer and Pohl [25] have also proposed a probe that uses a surface
plasmon generated at a subwavelength metallic sphere.
When light is introduced into the subwavelength metallic sphere, the
dipole moment p inside the sphere is given by
ε(ω) − 1
p = 4πε0 a3 E0 , (4)
ε(ω) + 2
where ε(ω) is the dielectric constant of the sphere, ω is the angular frequency
of the light, a is the radius of the sphere, ε0 is the dielectric constant of
vacuum, and E0 is the magnitude of the applied electric field. The field Es
at the surface of the sphere is the superposition of the field generated by the
dipole and the applied field E0 , and is given by
ε(ω)
Es = 3 E0 . (5)
ε(ω) + 2
As shown in this equation, Es is enhanced when the light frequency satisfies
the equation
Re[ε(ω)] = −2 . (6)
Fig. 5. Relation between the enhancement factor and the wavelength for a gold
sphere calculated by a dipole model
for a gold sphere. In this calculation, the experimental value given in [34] was
used as the dielectric constant of the gold. As shown in this figure, the field
at the surface became greater when the wavelength was near 500 nm.
The peak corresponds to a resonance condition of the collective oscillation
of the electrons in the metallic sphere. When the sphere is illuminated by
light, charges in the sphere move to the surface (Fig. 6). These charges give
rise to an internal field Ei , hence a restoring force acts on the charges. This
force causes a charge oscillation that has a resonance frequency satisfying
(6). If the light frequency corresponds to the resonance frequency of the
charge oscillation, light energy is strongly absorbed by the sphere and the
optical near field at the surface is enhanced. The quantum of the collective
oscillation of the charges in the small metallic particle is called a localized
surface plasmon or a localized plasmon.
A coaxial probe and a metallic pin probe have been proposed by, respectively,
Fischer [27] and Takahara et al. [28]. These probes have metallic cores, and
high efficiency is expected because there exists a mode propagating at the
surface of the metallic core, and this mode does not have a cutoff diameter.
Grober et al. [29] have proposed the use of a bow-tie antenna, which has been
used for microwave applications, as an optical near-field generator. When an
electromagnetic wave is incident upon such an antenna, current is induced in
the two arms, and charges accumulate at the apices. These charges work as
a dipole and a strong electromagnetic field is generated at the gap. To check
the feasibility of this approach, they used microwaves with a wavelength of
13.6 cm (2.2 GHz), and used a dipole antenna to observe the electromagnetic
field localized at the gap.
4.1 Principle
Fig. 7. Schematic of the near-field optical head with a wedge-shaped metallic plate
To check the feasibility of the proposed head, the distribution of the optical
near field was calculated through a finite difference time-domain (FDTD)
calculation [32, 33].
In this calculation, the incident light was assumed to be a plane wave.
Experimental values given in [34] were used as the dielectric constants of the
metal. The cell size was made large near the boundary and small near the
object to be calculated to minimize the computer memory required. When
the object was the wedge-shaped metallic plate, the cell size was designed
to have the smallest value near the apex, and the minimum cell size was
1/10 of the radius of the apex. As the boundary condition, Mur’s absorption
boundary condition was used. A parallel computer with 24 processors was
used for the calculation.
To begin with, for reference we will show the distribution generated near a
subwavelength aperture and a circular metallic plate.
Aperture
Fig. 9. Intensity distribution of the optical near field generated near an aperture
created in a metallic film. It was assumed that the metal was gold, the aperture
diameter was 30 nm, and the thickness of the metallic film was 100 nm. The dis-
tribution was calculated on a plane 5 nm from the metallic plate
Figure 11 shows the intensity distribution of the optical near field near a
circular metallic plate as shown in Fig. 10. For this calculation, we assumed
that both the diameter a and the thickness t were 30 nm, the metal was
gold, and the wavelength of the incident light was 650 nm. The polarization
of the incident light was in the direction of the arrows shown in Fig. 11.
The distance d between the observation plane and the aperture was 5 nm.
As shown in this figure, the near-field intensity became high at the edges A
Near-Field Optical-Head Technology 107
Fig. 11. Intensity distribution of the optical near field generated near a circular
metallic plate. It was assumed that the metal was gold, and the circle diameter and
the metallic film thickness were both 30 nm. The distribution was calculated on a
plane 5 nm from the metallic plate
Fig. 12. Intensity of the optical near-field generated near the circular metallic plate
as a function of the incident light wavelength. The intensity was calculated on a
plane 2 nm from the metallic plate
and A of the metallic plate. This is because the charges in the metallic plate
accumulated at these edges. The intensity at the peak was 3.5 times that of
the incident light. Figure 12 shows the dependence of the optical near-field
108 T. Matsumoto
In this subsection, we will present our simulation results for the wedge-shaped
metallic plate placed in air, and describe the fundamental properties of the
optical near field generated near the metallic plate, such as the intensity
distribution, efficiency, and resonance property of the plasmon.
Near-Field Distribution
Figure 14 shows the near-field distribution near the metallic plate in the shape
of a sector as shown in Fig. 13a. For this simulation, we assumed that the
metal was gold, the radius r of the apex was 10 nm, the length l was 150 nm,
the thickness t was 30 nm, the apex angle θ was 60◦ , and the wavelength of
the incident light was 690 nm. The distribution was calculated on a plane
5 nm from the metallic plate. The unit of intensity was the ratio between the
intensity (power density) of the optical near field and that of the incident
light. As shown in this figure, a strong optical near field was generated at
the apex, and the peak intensity was 350 times that of the incident light.
This intensity was 8750 times that of the 30-nm aperture created in the
metallic film, and 85 times that of the circular metallic plate. The distribution
measured at half-maximum was 19 nm wide in the x direction and 25 nm
wide in the y direction.
For this type of near-field optical head, the polarization must be in the
x direction as shown in Fig. 13. Figure 15 shows a near-field distribution
when the polarization of the incident light was in the y direction. As shown
Fig. 13. Simulation model for the wedge-shaped metallic plate: (a) metallic plate
in the shape of a sector, (b) metallic plate in the shape of a triangle, (c) side view
Near-Field Optical-Head Technology 109
Fig. 14. Intensity distribution of the optical near field generated near a metallic
plate in the shape of a sector. It was assumed that the radius r of the apex was
10 nm, the length l was 150 nm, the thickness t was 30 nm, the apex angle θ was
60◦ , and the incident light wavelength was 690 nm. The unit of intensity was the
ratio between the intensity of the optical near field and that of the incident light
Fig. 15. Intensity distribution of the optical near field when the polarization was
in the y direction
110 T. Matsumoto
in this figure, the optical near field became strong along the side edges of the
metallic plate, but weak at the apex. This is because the charges oscillate in
the direction parallel to the y axis; hence the charges accumulate along the
side edges instead of accumulating at the apex.
Fig. 16. Intensity distribution of the optical near field generated near a metallic
plate in the shape of a triangle. It was assumed that the metal was gold, the radius
r of apex A was 10 nm, the radius r of apices B was 30 nm, the apex angle θ was
60◦ , the thickness t was 30 nm and the incident light wavelength was 690 nm. The
distribution was calculated on a plane 5 nm from the metallic plate
Figure 16 shows the near-field distribution when the shape of the metallic
plate was a triangle. For this calculation, we assumed that the metal was
gold, the radius r of apex A was 10 nm, the radius r of apices B was 30 nm,
the apex angle θ was 60◦ , the thickness t was 30 nm and the wavelength of
the incident light was 690 nm. The distribution was calculated on a plane
5 nm from the metallic plate. As shown in this figure, a strong optical near-
field was also generated at apex A as with the metallic plate in the shape of
the sector. The distribution at the apex was the same as that of the metallic
plate in the shape of the sector. Note that the radius r of apices B should
be larger than the radius r of the apex A, otherwise charges also concentrate
at apices B, and a high optical near-field intensity is also generated there.
Near-Field Optical-Head Technology 111
Because the optical near field is localized near the metallic plate, the size of
the optical spot depends on the distance from the metallic plate d. Figure 17
shows the spot size as a function of the distance from the metallic plate for a
metallic plate in the shape of the sector, and Fig. 18 shows the distributions
for distances of 2, 5, and 10 nm. In Fig. 17, the solid line represents the
width measured at half-maximum in the x direction dx, and the broken line
represents the width measured at half maximum in the y direction dy. As
shown, the spot size decreased as the distance became smaller. For example,
the spot size was dx = 9 nm, dy = 16 nm when the distance was 2 nm, while
it was dx = 28 nm, dy = 30 nm when the distance was 10 nm.
70
60
50
Spot size [nm]
40 dy
30
20 dx
10
0
0 5 10 15 20
Distance [nm]
Fig. 17. Dependence of the spot size on the distance from the metallic plate when
the radius of the apex was 10 nm
The size of the optical spot also depends on the radius r of the apex. Figure 19
shows the size of the optical spot as a function of the radius of the apex
simulated for the metallic plate in the shape of the sector. The solid line
shows the width measured at half-maximum in the x direction dx, and the
broken line shows the width measured at half-maximum in the y direction
dy. The distribution was calculated on a plane 5 nm from the metallic plate.
As shown, the spot size decreased linearly as the radius of the apex became
smaller.
112 T. Matsumoto
Fig. 18. Intensity distribution of the optical near field generated near a metallic
plate in the shape of a sector for a distance of (a) 2 nm, (b) 5 nm, (c) 10 nm
50
40
dy
Spot size [nm]
30
Efficiency
where pnear is the power density (intensity) of the optical near field, pin is
that of the incident light, and S and S are the areas where the power density
is higher than half-maximum. The efficiency calculated for the near-field dis-
tribution in Fig. 14 was about 15% if we assume that the incident beam was
a Gaussian beam with a full width at half-maximum of 1 µm (corresponding
to an optical spot focused by a lens with an N A of 0.35). When we calculated
the near-field distribution, we assumed the incident light was a plane wave.
Therefore, in this estimation, the diameter of the incident beam was assumed
to be 1 µm so that the incident light intensity would be nearly constant at the
metallic plate. However, under such a condition, the energy loss is large be-
cause much of the light is not incident on the metallic plate. Higher efficiency
should be attainable by reducing the diameter of the incident light.
In the case of the aperture, the near-field intensity became lower, and the
efficiency decreased rapidly as the optical spot became smaller. However, for
the head using the metallic plate, the near-field intensity increased when the
spot size became smaller. As a consequence, the efficiency did not decrease
as rapidly as the aperture. The solid line in Fig. 20 shows the intensity as a
function of the apex radius simulated for the metallic plate in the shape of the
sector. The intensity represents the peak intensity in a distribution calculated
on a plane 5 nm from the metallic plate. As shown, the near-field intensity
rose as the apex radius became smaller. This is because the oscillating charges
concentrated in a smaller area, and the charge density increased when the
apex radius became smaller. The broken line in Fig. 20 shows efficiency as a
function of the apex radius. The efficiency was about 10% even if the radius
was 5 nm (corresponding to an optical spot 15 nm by 18 nm).
The wavelength of the incident light must correspond to the resonance wave-
length of the plasmon to generate a strong optical near field. However, the
available wavelength is limited if we use a semiconductor laser as the light
source. Therefore, we need to tune the resonance wavelength of the plasmon
to the wavelength of the semiconductor laser. For the wedge-shaped metallic
plate, the resonance wavelength of the plasmon can be adjusted by changing
the material or the length l of the metallic plate. In this subsection, we will
describe how the resonance wavelength depends on these parameters.
114 T. Matsumoto
600 25
500 Efficiency
20
Efficiency [%]
400
Intensity
15
300
Intensity 10
200
5
100
0 0
0 10 20
Radius [nm]
Fig. 20. Dependence of the optical near-field intensity and efficiency on the apex
radius. It was assumed that the length l, the apex angle θ, and the incident light
wavelength were 150 nm, 60◦ , and 690 nm, respectively. The intensity represents the
peak intensity in a distribution measured on a plane 2 nm from the metallic plate.
For the calculation of efficiency, the incident light was assumed to be a Gaussian
beam with a full width at half-maximum of 1 µm
Figure 21 shows the relation between the optical near-field intensity and
the incident light wavelength for metallic plates made of gold, silver, and
aluminum. For this calculation, we assumed the shape of the metallic plate
was a sector with an apex radius r of 20 nm, a length l of 100 nm, a thickness
Fig. 21. Relation between the near-field intensity and the wavelength for metallic
plates made of silver, gold, and aluminum. The shape of the metallic plate was
assumed to be a sector with an apex radius r of 20 nm, a length l of 100 nm, a
thickness t of 30 nm, and an apex angle θ of 60◦ . The closed triangles, closed circles,
and closed squares represent the values for gold, silver, and aluminum, respectively.
The intensity represents the peak intensity in the distribution calculated on a plane
2 nm from the metallic plate
Near-Field Optical-Head Technology 115
t of 30 nm, and an apex angle θ of 60◦ . The intensity represents the peak
intensity in a distribution calculated on a plane 2 nm from the metallic plate.
The closed triangles, closed circles, and closed squares represent the values for
gold, silver, and aluminum, respectively. The resonance wavelength of gold
was about 650 nm, which is near the wavelengths of a red semiconductor laser
(635 nm, 650 nm, and 670 nm). The resonance wavelength of aluminum was
370 nm, which is near the wavelength of a blue semiconductor laser (405 nm).
In the case of silver, the resonance wavelength was 550 nm. Although, there
is no semiconductor laser of such wavelength, the resonance wavelength can
be shifted toward the red or near-infrared regions by changing the length l
as we explain next.
The near-field intensity at the resonance wavelength became larger in the
order Al < Au < Ag. This indicates that the near-field intensity becomes
large when the imaginary part of the dielectric constant becomes small. This
is because the damping of the charge oscillation, which causes energy loss,
becomes small when the imaginary part of the dielectric constant is small.
To realize high efficiency, we therefore need to use a metal whose dielectric
constant has a small imaginary part.
Figure 22 shows the resonance curves when the length l was varied from
100 nm to 600 nm. For this calculation, we assumed that the shape of the
metallic plate was a sector with an apex radius r of 20 nm, a thickness t of
30 nm, and an apex angle θ of 60◦ , and that the metal was silver. The intensity
Fig. 22. Relation between the near-field intensity and the wavelength for metallic
plates of various lengths. It was assumed that the shape of the metallic plate was
a sector with an apex radius r of 20 nm, a thickness t of 30 nm, and an apex angle
θ of 60◦ , and that the metal was silver. The intensity represents the peak intensity
in the distribution calculated on a plane 2 nm from the metallic plate
116 T. Matsumoto
In the previous calculation, the metallic plate was placed in air. However, in
recording devices, the metallic plate is formed on the slider (a transparent
substrate), and placed near the recording medium. In this subsection, we will
describe the resonance property of the plasmon and the near-field distribution
in such a situation.
Simulation Model
To realize high durability, the metallic plate should be embedded in the slider
surface as shown in Fig. 23. Thus, we assumed that the metallic plate was
embedded in a glass substrate with a refractive index of 1.45. The recording
medium was assumed to be a thin film of TbFeCo or GeSbTe with a thickness
of 20 nm, and the separation between the recording medium and the metallic
plate was assumed to be 10 nm. For the metallic plate, we assumed that the
metal was gold, the apex radius r was 10 nm, the length l was 150 nm, the
thickness t was 30 nm, and the apex angle θ was 60◦ .
Fig. 23. Simulation model when the recording medium and slider were placed near
the metallic plate: (a) when there was no recording medium, (b) when there was a
recording medium
Near-Field Optical-Head Technology 117
First, we show the influence of the glass substrate on the resonance prop-
erty of the plasmon. Figure 24 shows the resonance curve when there was
no recording medium. The solid line represents the resonance curve for the
metallic plate embedded in the glass substrate and the broken line represents
the resonance curve when the metallic plate was placed in air. The intensity
represents the peak intensity in the distribution calculated on a plane 2 nm
from the metallic plate. As shown, the resonance wavelength shifted by about
50 nm toward a longer wavelength when the metallic plate was embedded in
the glass substrate. We believe this was because the dipoles in the dielec-
tric substrate interact with the charges in the metallic plate, and weaken the
restoring force working on the charges.
Figure 25 shows the resonance curves when the recording medium was
placed near the metallic plate embedded in the glass substrate. The solid line
represents the resonance curve when the TbFeCo recording medium was used,
the broken line represents the resonance curve when the GeSbTe medium
was used, and the dash-dot line represents the resonance curve when there
was no recording medium. The intensity represents the peak intensity in the
distribution calculated on a plane 2 nm from the metallic plate. As shown, the
resonance wavelength was shifted toward a longer wavelength by interaction
with the recording medium. The shift was about 50 nm for both TbFeCo and
GeSbTe.
Fig. 24. Relation between the near-field intensity and the wavelength for the metal-
lic plates embedded in a glass substrate. The solid line represents the resonance
curve for the metallic plate embedded in the glass substrate and the broken line
represents the resonance curve when the metallic plate was placed in air. The in-
tensity represents the peak intensity in the distribution calculated on a plane 2 nm
from the metallic plate
118 T. Matsumoto
Fig. 25. Relation between the near-field intensity and the wavelength when a
recording medium was placed near the metallic plate. The solid line represents the
resonance curve for a TbFeCo recording medium, the broken line represents the res-
onance curve for a GeSbTe recording medium, and the dash-dot line represents the
resonance curve when there was no recording medium. The separation between the
recording medium and the metallic plate was assumed to be 10 nm. The intensity
represents the peak intensity in the distribution calculated on a plane 2 nm from
the metallic plate
Near-Field Distribution
Figure 26 shows the intensity distribution of the optical near-field when there
was no recording medium, and Fig. 27 shows the intensity distribution when
the TbFeCo recording medium was placed near the metallic plate. The wave-
length of the incident light was assumed to be at the resonance wavelength
of the plasmon – 780 nm when there was no recording medium, and 830 nm
when there was a recording medium. Both distributions were calculated on a
plane 10 nm from the metallic plate (on the surface of the recording medium
when a recording medium was used). As shown in these figures, the opti-
cal spot was semicircular when the recording medium was used, but was
crescent-shaped when there was no recording medium. With the recording
medium, the distribution measured at half-maximum was about 30 nm wide
in both the x and y directions. If we compare this spot size with that of a 4.7-
GB, DVD-RAM (full width at half-maximum of the optical spot = 580 nm,
recording density = 3.7 Gb/in2 ), this spot size corresponds to a recording
density of approximately 1 Tb/in2 . The calculated efficiency for this distri-
bution was 20% when we assumed that the incident beam was a Gaussian
beam with a full width at half-maximum of 1 µm.
When we placed the recording medium near the metallic plate, the in-
tensity of the optical near field increased by a factor of about 2.5. This was
because the decay length of the optical near field changes when the recording
Near-Field Optical-Head Technology 119
Fig. 28. Intensity of the optical near field as a function of the distance from the
metallic plate. The solid line represents the curve when a recording medium was
used, and the broken line represents the curve when there was no recording medium
medium is placed near the metallic plate. Figure 28 shows the dependence of
the near-field intensity on the distance from the surface of the metallic plate.
The solid line represents the curve when a recording medium was used, and
the broken line represents the curve when there was no recording medium.
The near-field intensity decreased exponentially as the distance increased,
and its decay length increased when the recording medium was placed near
the metallic plate. We believe this was because image charges were induced in
the recording medium, and the intensity of the optical near field between the
head and the recording medium became greater because of the interaction of
the charges induced at both the apex of the metallic plate and the recording
medium.
When a wedge-shaped metallic plate is used, the spot size can be further
decreased by reducing the apex radius or the spacing between the head and
the recording medium. An alternative way to decrease the spot size is to use
two metallic plates (Fig. 29). When the two metallic plates are illuminated by
light, an optical near field is generated between the two apices because of the
interaction of charges concentrated at these apices. This mechanism is similar
to that of the bow-tie antenna used at microwave frequencies. However, at
optical frequencies, the intensity of the optical near field can be enhanced
when the plasmon is excited in the metallic plates by optimizing the material
and the size of the metallic plate.
Figure 30 shows the intensity distribution of the optical near field gener-
ated near the two metallic plates. For this calculation, we assumed that the
shape of the metallic plates was a sector, the metal was silver, the gap width
Near-Field Optical-Head Technology 121
Fig. 30. Intensity distribution of the optical near field generated near the double
metallic plate. It was assumed that the metallic plates had a sector shape, the metal
was silver, the gap width g was 5 nm, the apex radius r was 20 nm, the length l was
200 nm, the thickness t was 30 nm, the apex angle θ was 60◦ , and the wavelength
was 830 nm. The distribution was calculated on a plane 2 nm from the metallic
plates
g was 5 nm, the radius r of each apex was 20 nm, the length l was 200 nm,
the thickness t was 30 nm, the apex angle θ was 60◦ , and the wavelength was
830 nm. The distribution was calculated on a plane 2 nm from the metallic
plates. As shown, a strong optical near field was generated at the gap. The
peak intensity was about 2000 times that of the incident light. The distribu-
122 T. Matsumoto
Fig. 31. Relation between the near-field intensity and the wavelength for the double
metallic plates with a gap width of 24 nm, 16 nm, and 8 nm
As described in Sect. 4.1, a near-field optical head with a metallic plate can
be easily fabricated by using electron beam lithography. Figure 32 shows an
example of fabrication process. In this process, the head is fabricated using a
sacrificial substrate to embed the metallic plate in the slider. First, an array of
metallic plates is fabricated, using electron beam lithography, on a sacrificial
substrate coated with a sacrificial film (a thin metallic film) in steps (a) to
(c). Secondly, a photoresist pattern is formed on the sacrificial substrate to
make pads near the metallic plates, and a dielectric film (Al2 O3 ) is formed
on the substrate by sputtering. Thirdly, a glass substrate was glued to the
Near-Field Optical-Head Technology 123
Fig. 33. SEM image of the metallic plates: (a) one metallic plate, (b) two metallic
plates
top of the structure, and then the photoresist and the sacrificial substrate are
removed. Figure 33 shows a scanning electron microscope image of metallic
plates fabricated by the method described above. In this fabrication, gold
was used as the metal, and the thickness was 30 nm. As shown in this figure,
the radius of the apex was about 10 nm, and the gap width was 20 nm for
the two metallic plates. Based on our simulation results, a recording density
of about 1 Tb/in2 should be achievable by using the fabricated near-field
optical head.
124 T. Matsumoto
5 Summary
Near-field optical recording is a promising way to realize a recording density
of over 1 Tb/in2 . In this chapter, we focused on the near-field optical head,
which is a key device for near-field optical recording. First, we explained the
technical issues regarding the near-field optical head and introduced some
solutions to these issues. We focused on a highly efficient near-field optical
head that uses a wedge-shaped metallic plate, and described its optical prop-
erties based on a simulation using a finite-difference time-domain method.
The simulation results confirmed that a strong optical near field is generated
at the apex of the metallic plate when a plasmon is excited in the metallic
plate. When a TbFeCo recording medium was placed 10 nm from the near-
field optical head, the size of the optical spot was 30 nm, which corresponds
to an areal recording density of approximately 1 Tb/in2 . The efficiency was
20% if we assume that the incident beam was a Gaussian beam with a full
width at half-maximum of 1 µm. Furthermore, we discussed an optical head
using two metallic plates. We confirmed through our simulation that a highly
localized optical near field was generated at the gap when the plasmon was
excited in the metallic plates. The distribution was 5 nm by 5 nm when the
two apices were separated by 5 nm. These results indicate that the recording-
density limit now being encountered in conventional recording devices can be
overcome with this type of near-field optical head.
To realize such advanced recording devices, development of a recording
medium suitable for near-field optical recording is also important. For in-
stance, in the case of magnetic recording medium, a recording medium with
a small grain size and thin cover layer will be needed, and tuning of the mag-
netic and thermal characteristics will be necessary. High-precision tracking
servo technology will also be required. For example, a track pitch of less than
40 nm is required for a recording density of more than 1 Tb/in2 ; that of the
current Blu-ray disk is 320 nm. We believe that these technologies can be
realized on the basis of technology that has been developed for conventional
optical or magnetic recording devices. However, various technical challenges
will have to be overcome to realize such technology because the required size
and precision are less than 1/10 those of today’s conventional technology.
Acknowledgements
References
1. Near-field Nano/Atom Optics and Technology, ed. by M. Ohtsu (Springer,
Tokyo 1998)
2. S.M. Mansfield, G. S. Kino: Appl. Phys. Lett. 57, 2615 (1990)
3. B.D. Terris, H.J. Mamin, D. Rugar, W.R. Studenmund, G.S. Kino:
Appl. Phys. Lett. 65, 388 (1994)
4. M. Shinoda, K. Saito, T. Ishimoto, T. Kondo, A. Nakaoki, M. Furuki,
M. Takeda, M. Yamamoto: In: Technical Digest of Optical Data Storage 2003
(Canada 2003) 120
5. E. Betzig, J.K. Trautman, R. Wolfe, E.M. Gyorgy, P.L. Finn, M.H. Kryder,
C.-H. Chang: Appl. Phys. Lett. 61, 142 (1992)
6. S. Hosaka, T. Shintani, M. Miyamoto, A. Kikukawa, A. Hirtsune, M. Terao,
M. Yoshida, K. Fujita, S. Krammer: J. Appl. Phys. 79, 8082 (1996)
7. S. Jiang, J. Ichihashi, H. Monobe, M. Fujihira, M. Ohtsu: Opt. Commun. 106,
173 (1994)
8. H. Saga, H. Nemoto, H. Sukeda, M. Takahashi: Jpn. J. Appl. Phys. 38, Part
1, 1839 (1999)
9. H. Sukeda, H. Saga, H. Nemoto, Y. Ito, C. Haginoya: IEEE Trans. Magn. 37,
1234 (2001)
10. M. Vasz-Iravani, R. Toledo-Crow, Y. Chen: J. Vac. Sci. Technol. A 11, 742
(1993)
11. E. Betzig, P.L. Finn, J.S. Weiner: Appl. Phys. Lett. 60, 2484 (1995)
12. M.B. Lee, M. Kourogi, T. Yatsui, K. Tsutsui, N. Atoda, M. Ohtsu: Appl. Opt.
38, 3566 (1999)
13. F. Isshiki, K. Ito, K. Etoh, S. Hosaka: Appl. Phys. Lett. 76, 804 (2000)
14. T. Yatsui, M. Kourogi, K. Tsutsui, J. Takahashi, M. Ohtsu: Opt. Lett. 25, 1279
(2000)
15. A. Partovi, D. Peale, M. Wutting, C.A. Murray, G. Zydzik, L. Hopkins, K. Bald-
win. W.S. Hobson, J. Wynn, J. Lopata, L. Dhar, R. Chichester, J.H.-J. Yeh:
Appl. Phys. Lett. 75, 1515 (1999)
16. K. Goto: Jpn. J. Appl. Phys. 37, Part 1, 2274 (1998)
17. K. Kato, S. Ichihara, M. Oumi, H. Maeda, T. Niwa, Y. Mitsuoka, K. Nakajima,
T. Ohkubo, K. Itao: Jpn. J. Appl. Phys. 42, Part 1, 5102 (2003)
18. T.W. McDaniel, W.A. Challener: In: Technical Digest of MORIS’02 (France
2002) 117
19. T. Saiki, S. Mononobe, M. Ohtsu, N. Saito, J. Kusano: Appl. Phys. Lett. 68,
2612 (1996)
20. T. Yatsui, M. Kourogi, M. Ohtsu: Appl. Phys. Lett. 73, 2090 (1998)
21. U.C. Fischer et al.: In: Near-Field Optics ed. by D.W. Pohl, D. Courjon (Kluwer
Academic, Dordrecht 1993) 255
22. J. Koglin, U.C. Fischer, H. Fuchs: Phys. Rev. B 55, 7977 (1997)
23. H. Raether: Surface Plasmon on Smooth and Rough Surface and on Grating
(Springer-Verlag, Berlin 1988)
24. T. Thio, K.M. Pellerin, R.A. Linke, H.J. Lezec, T.W. Ebbesen: Opt. Lett. 26,
1972 (2001)
25. U.C. Fischer, D.W. Pohl: Phys. Rev. Lett. 62, 458 (1989)
26. T. Matsumoto, T. Ichimura, T. Yatsui, M. Kourogi, T. Saiki, M. Ohtsu:
Opt. Rev. 5, 369 (1998)
126 T. Matsumoto
1 Introduction
The data storage density is advancing because of the requirement of infor-
mation data growth as shown in Fig. 1. A 1 Tb/inch2 density will be re-
quired around 2010. Physical limits, however, are predicted both for magnetic
recording and for optical means. To overcome the optical limit (diffraction
limit), near-field optics is required [1].
The magnetic recording density is currently increasing at rates of up to
100 % per year. To obtain the higher densities, grain sizes of the conventional
continuous magnetic films should be reduced to maintain a necessary signal-
to-noise ratio. The small grain sizes, however, reduce the thermal stability
of the magnetization of each bit. This, referred to as superparamagnetism,
can be overcome by increasing KuV /kT , where Ku is an energy barrier to
reversal per grain volume (anisotropic energy), V is a volume per grain, k
is the Boltzmann constant, and T is temperature. High Ku materials such
signal contrasts are 1% and 10%, respectively. In future near-field storage de-
vices, signal contrast will be insufficient to achieve practical read-out speeds.
As for photochromic materials, although better signal contrast is expected,
read out without destroying the recorded information has been a major issue.
We showed that amorphous organic dye molecules formed droplet-like nm-
scale structures and were aligned by the AASA method [17]. It is possible to
inject and confine charges into a single dot using AFM [18, 19]. Fluorescence
changes in dye molecules associated with charge injection into their thin films
was observed [20]. Fluorescence measurement has such a high signal-to-noise
ratio that it can be observed even in a single molecule. Therefore, read out
using fluorescence is expected to achieve high-speed scanning. Furthermore,
electrical writing, namely charge confinement, will prevent record rupture at
read out.
A Ni master disk possessing spiral patterns with 60-, 100-, 150-, 200- and
250-nm width lands and a 400-nm width groove (110-nm depth) was pressed
to a resist film on a CoCrPt film to transfer the spiral patterns at room
temperature and at a pressure of 1000 bar. The master disk was prepared by
the photoprocess as shown in Fig. 9.
132 K. Naito et al.
Fig. 10. AFM images of spiral patterns of the Ni master (upper images) and of
the resist film (lower images). The lengths are the land widths of the master disk
Fig. 11. Phase-separation dot patterns of the diblock copolymers in grooves with
different widths: (a) 60 nm, (b) 150 nm, (c) 200 nm, (d) 250 nm
The AFM images of spiral patterns transferred to the resist film are shown
in Fig. 10. The land-groove spiral patterns were precisely transferred.
A diblock copolymer solution was cast into the obtained grooves, and then
annealed to prepare self-assembling dot structures aligned along the grooves.
Figure 11 indicates dot structures in grooves with different widths. One to
four dot lines could be obtained by changing the groove width. The groove
width should be controlled to obtain the regularly aligned dot lines along the
grooves. In the figure, the groove widths were considered rather suitable for
the two and four dot lines. Figure 12 shows the phase-separation dot patterns
of another diblock copolymer in grooves. The dot size was 15 nm and the dot
pitch was 35 nm. This small dots were obtained by using the low molecular
weight diblock copolymer.
The PMMA dots were selectively removed by the oxygen plasma treat-
ment [12]. The resulting holes were filled by spin-on-glass (SOG). The lower
magnetic films were patterned by ion milling using the SOG dots as a mask.
Nano-Optical Media for Ultrahigh-Density Storage 133
Fig. 13. (a) SEM image of the patterned magnetic media with a 40-nm diameter,
(b) the whole image of the patterned media disk prepared on a 2.5-inch HDD glass
plate
Figure 13a shows the SEM image of the patterned magnetic media with
40 nm diameter. There are fluctuations in size and in position of the magnetic
dots. The irregularities are probably caused by the diblock copolymer dot
disorder and by the nanopatterning process fluctuation. The whole image of
the patterned media disk prepared on a 2.5-inch diameter HDD glass plate
is shown in Fig. 13b. Interference colors based on the wide groove lines were
observed. Figure 14 shows the TEM image of the patterned dot with about
40 nm diameter. This image indicates that the magnetic dot is completely
separated.
134 K. Naito et al.
FePt is an ordered alloy that has a high Ku. High temperature is required
to obtain the ordered structure form the sputtered film as shown in Fig. 15.
We found that addition of Cu decreases the temperature [22].
Using an electron-beam mastering apparatus recently developed will en-
able reduction of the groove width to increase the density of magnetic cells.
It may be impossible to produce circumferential dot patterns with a 40 nm
diameter like those in the figure even by the present electron beam apparatus.
We used the AASA method to prepare the etching mask for the magnetic
films. For cheaper preparation, we prepared the Ni master disk itself by the
AASA method. The details will be shown elsewhere.
Figure 16 shows the M –H loops of the raw continuous film (CoCrPt) and
the patterned media. A coercive force and squareness ratio of the patterned
film increased compared to the continuous film. The magnetization of the
patterned media decreased due to the decrement of the material volume.
The squareness ratio was almost unity.
Figure 17 indicates the MFM images of the patterned magnetic media
after AC erase and after DC erase. Three magnetic dot lines are observed
corresponding to Fig. 11c and Fig. 13a. Bright and dark dots are randomly
Nano-Optical Media for Ultrahigh-Density Storage 135
Fig. 16. M –H loops of the raw continuous film (CoCrPt) and the patterned media
3 Organic-Dye-Patterned Media
3.1 Preparation
The organic molecule mainly used here was TTPAE or TPD, shown in Fig. 18.
These molecules give amorphous films as a hole transport layer in organic EL
devices [23]. The films were formed by vacuum evaporation onto thermally
oxidized p-Si substrates with a resistivity of 0.001–0.01 Ω cm. The thickness
of the SiO2 layer was 20 nm. The electrodes were fabricated on the back
of the Si substrates by evaporating Cr and Au. Before the evaporation of
molecules, the SiO2 /Si substrates were made hydrophobic by treatment in
an atmosphere of hexamethyldisilazane.
Figure 19 shows the contact-mode topographic image of the TTPAE film
observed using an AFM. It is clearly seen that the film is constructed of
spatially isolated droplet-like domain structures. The size and the density of
domains were varied according to the amount of deposition. The diameter,
height and density were varied from 20 to 300 nm, from 10 to 100 nm, and
from 1011 to 109 cm−2 , respectively. When the amount of deposition was
increased to the film thickness of more than 40 nm, the film structure changed
into a continuous layer with a flat surface. We verified that thin films of several
amorphous molecules, including TPD, also showed similar domain structures.
The observed uniform droplet-like domain structures suggest that the
film growth is based on the Volmer–Weber mechanism. In this mechanism,
growth characteristics are known to depend strongly on the surface energies
of the film and of the substrate. When the films were formed onto hydrophilic
SiO2 /Si substrates, it was found that the domain size was much larger than
films on the hydrophobic substrates.
The droplet size and alignment could be controlled by the AASA method.
The hole arrays were prepared on polymethylmethacrylate film by using
AFM. The pitch was 30 nm. TPD molecules were evaporated onto the ar-
rays. TPD amorphous dots were selectively formed on the holes with the same
size as shown in Fig. 20. Using a nanoimprint method should produce the
large-area hole arrays to obtain the large-area organic-dye-patterned media
as shown in Fig. 21.
Figure 22 shows a schematic diagram of the structure of the device used in this
chapter. To apply voltage to the organic films, Au films of 10 nm thickness
were fabricated on the substrates following the deposition of organic films. Au
films of 10 nm were selected as the thinnest providing good electrical contact
and sufficient fluorescence intensity from the underlying organic films.
Scanning of the applied voltage started at 0 V. Fluorescence measure-
ments were first obtained in the positive-voltage regions, and then in the
negative regions. The voltage scan speed was 0.02 Hz. All fluorescence mea-
surements were carried out in air and at room temperature.
In order to investigate the electrical properties of TTPAE films, displace-
ment-current measurements were performed. Since details of displacement-
current measurements were described in [24], only distinctive features are
shown here. Figure 23 shows a schematic diagram of the displacement-current
measurement. The measured current is expressed by the capacitance and scan
speed of the applied voltage. In the simplest case, the organic film behaves
as a dielectric and the measured current reflects the capacitance of the entire
dielectric film, that is, organic and SiO2 (CSiO2 +Org ). As for the donor or
acceptor molecules, carrier injection or ejection through metal/org interfaces
occur. If carrier injection through the metal/org interface occurs, the mea-
sured current will reflect the capacitance of the SiO2 (CSiO2 ) film. Therefore,
the displacement current is sensitive to carrier injection into the organic film.
Figure 24 shows the applied-voltage dependence of fluorescence in the
TTPAE film. These data were obtained in the first scan of voltage for positive
and negative voltages, respectively. As shown later, large hysteresis effects
were observed in fluorescence quenching. As shown in Fig. 24, the fluorescence
intensity decreases as applied voltage V increases at positive voltages. On the
other hand, at negative voltages, no variations with V are observed. Changes
in the fluorescence spectrum and excitation spectrum are shown in Fig. 25. It
was found that the observed fluorescence quenching did not follow any shift
of spectral peaks.
Figure 26 shows the displacement-current characteristic of a TTPAE film
under visible light irradiation. At positive bias voltages, hole injection and
ejection were clearly observed, and at negative bias voltages, no carrier injec-
tion was observed. These findings are consistent with the fact that TTPAE
is an electron-donor material. Figure 27 shows a comparison between fluores-
cence measurements and displacement-current measurements under UV ir-
radiation. Since both fluorescence variations and displacement-currents show
large hysteresis especially at negative bias voltages, the results for the first
cycle of voltage are shown separately from those for the second and third cy-
cles. At positive bias voltages, fluorescence and current into the TTPAE film
were observed to vary in the same manner. Both fluorescence variations and
displacement currents are reproducible at positive bias voltages. On the other
hand, at negative bias voltages, large hysteretic variations were observed in
both measurements. On the way to −25 V in the first cycle, no fluorescence
change or negative current into the TTPAE film was observed. On the way
back to 0 V, a positive current and fluorescence quenching were observed
at the same voltage. In the second or third cycle, two current peaks were
observed on the way to −25 V and one peak was observed on the way back
to 0 V. Fluorescence quenching was observed at lower voltages where current
was observed on the way to −25 V, and at the same voltage as the current
observed on the way back to 0 V.
Figure 28 shows the fluorescence-quenching efficiency for various thick-
nesses of TTPAE films. The fluorescence-quenching efficiency increased as
the film thickness decreased and had a peak value at a thickness of 15 nm.
In the case of films thicker than 15 nm (c), uniform film structures were
Nano-Optical Media for Ultrahigh-Density Storage 141
observed. On the other hand, in the case of thinner films, separated film
structures (b) or isolated dot structures (a) were observed. As for the sepa-
rated structures or the dot structures, the effective thickness of the film was
found to be thicker than the case for uniform films.
We directly observed fluorescence quenching for the TTPAE dots by a
near-field optical microscope. Charge injection is carried out through contact
electrification by a gold-coated AFM tip with a positive bias. The charges
were found to be very stable over a long period. This phenomena will be
the same as in a flash silicon memory, where electrons are stored in a float-
ing gate. Injected holes can be extracted by a negatively biased tip. Holes
were injected to diamine dots at 3 V. These holes were extracted at −2 V.
Nano-Optical Media for Ultrahigh-Density Storage 143
4 Conclusion
We have described the circumferential magnetic patterned media with a 40-
nm diameter on a 2.5-inch diameter glass plate. The media were prepared
by an artificially assisted self-assembling (AASA) method, which includes
simple nanopatterning of land-groove spirals using a nanoimprint and self-
assembling of the diblock polymer aligned in the grooves. Magnetic films were
etched by ion milling using the nanodot structures of the diblock copolymer
as a mask. Fluctuations of size and position of the patterned cell, observed for
the present media, should be reduced for the realization of patterned media.
We have also investigated organic-dye-patterned media prepared by the
AASA method. Fluorescence quenching of TTPAE dots with charge injection
applying bias voltages. Displacement-current measurements show that the
quenching efficiencies depend strongly on the polarities of injected charges
and that quenching cannot be explained by the external-field-induced exciton
dissociation. Charge injection properties show that the hole injections play
an essential role in fluorescence quenching. A 60% fluorescence-quenching
efficiency is achieved. This will be a novel optical storage medium with high
contrast and high signal-to-noise ratio.
These nm-scale dot structures should prove to be effective in the applica-
tion of near-field optics to ultrahigh-density storage media.
Acknowledgements
The authors are very grateful to T. Maeda for providing them with the sam-
ples of the magnetic films. They also acknowledge useful discussions with
H. Yoda, J. Akiyama, T. Hiraoka and Y. Yanagita.
References
1. Modified by K.N. from the roadmap (Fig. 2.3.1.1) In: Research survey on near-
field light technology, OITDA (2000)
2. S.Y. Chou, M.S. Wei, P.R. Krauss, P.B. Fisher: J. Appl. Phys. 76, 6673 (1994)
3. R.L. White, R.M.H. New, RE.F.W. Pease: IEEE Trans. Magn. 33, 990 (1997)
144 K. Naito et al.
1 Introduction
1.1 What are Optical Near Fields?
Since the insight of Synge in 1928 [1], the localization of light near a material
system has been recognized as one of the unique properties of an optical near
field, and for more than two decades a variety of theoretical and experimental
studies have been intensively conducted [2–9]. At the beginning of this chap-
ter, we explain the spatial localization property of light, using two examples,
in order to give a physical insight into the optical near-field phenomenon.
The first example is evanescent light, which is generated on a dielectric
surface, as is well known, when the incident angle θ of a plane wave of light
exceeds the critical angle θc and total internal reflection occurs as shown in
Fig. 1 [2, 10]. The penetration or attenuation depth normal to the surface
(along the z axis) is of the order of the wavelength of light, which indicates
the spatial localization of light along the z axis. This differs from ordinary
propagating light, and much attention has been paid to its characteristics [11–
27]. Note that the total internal reflection is a macroscopic phenomenon, and
that it is valid only on a scale at which the refractive index n of bulk material
has meaning.
Another example is illustrated in Fig. 2 [2,9]. Figure 2a shows a plane wave
of light incident on a plane with a small hole of radius a that is much smaller
than the wavelength λ of light. Little of the propagating far-field light passes
through the hole, but light exists as an optical near field near the aperture on
the other side. The localization range of the optical near field is much smaller
than the wavelength of light and is determined by the aperture size a. This
kind of field, which a scanning near-field optical microscope (SNOM) uses at
the tip of an optical near-field probe, has been investigated in detail [28–32].
A more general situation is depicted in Fig. 2b. By irradiating a dielectric
that is much smaller than the wavelength of the incident light, optical near
fields coiling around the object are generated together with scattering far
fields. The localization range of the optical near fields is independent of the
wavelength of the incident light and is determined by the size of the object:
it depends on the radius a of the dielectric sphere as shown in Fig. 2b. Such
spatial localization reflecting the structure of the source of the fields (i.e.,
146 A. Shojiguchi et al.
Fig. 1. (a) The infinite planar boundary between a dielectric surface and a vacuum.
Incident light from the dielectric surface is generally reflected and refracted. (b)
Evanescent light on the dielectric surface generated by total internal reflection of
the incident light
Fig. 2. Optical near fields in various boundary conditions. (a) An optical near field
leaking from a hole that is much smaller than the wavelength of light. The leaking
field attenuates as it leaves the aperture. The localization range of light is much
smaller than the wavelength of light, depending on the size of the aperture. (b) An
optical near field surrounding a small dielectric sphere and the scattered optical far
field produced by an incident optical far field. The localization range of the optical
near field is determined by the sphere radius and is independent of the wavelength
of light
e−γr
φ= . (1)
γr
The fields satisfying the spherical boundary condition can be expanded in
terms of plane waves, or can be expressed using angular spectrum repre-
sentation. From a numerical analysis, it follows that the distribution of the
angular spectrum has a peak at γ = 1/a [2]. Applying the result to (1), we see
that the localization range of the optical near fields around a sphere is of the
order of the radius a, i.e., the size of the material system. Here, it is obvious
that an optical near field should be derived consistently from microscopic-
induced polarization in matter because they are mutually related, and several
approaches in this direction are outlined in the following subsection 1.2.
The main difference between a material system coupled with optical near
fields and a material system coupled with propagating far fields, is that one
is coupled locally and the other is coupled globally, as depicted schematically
in Fig. 3.
Since the propagating field is usually expanded in terms of a plane-wave
basis, photons emitted as plane waves propagate from one site to distant
A Phenomenological Description of Optical Near Fields 149
Fig. 3. Schematic drawing of (a) a global and (b) a local coupling systems. Each
element of the system in (a) can interact directly with arbitrary elements at a
distant site, while the elements at the near site can interact only with each other
in (b)
sites in the material system and can also excite the far sites (global excita-
tion [50]), as shown in Fig. 3a. By contrast, as shown in Fig. 3b, “near-field
photons” or “localized photons” can move only to nearby sites and excite
them (individual excitation). Generally speaking, the dynamics of a locally
coupled system differ from those of a globally coupled system with respect
to the relaxation speed towards equilibrium [51]. Therefore, we expect the
dynamics and equilibrium states for a locally coupled system to differ from
those of an ordinary globally coupled system, such as a propagating field–
matter system [52]. From this perspective, it is important to investigate the
nanometric material system that is interacting with optical near fields via
individual excitations [46–48, 53].
This chapter presents a model of a material system, i.e., for a quantum-
dot system, interacting with optical near fields, and discusses the dynamics
of the electric dipole moment and the radiation properties of the system in
detail. For comparison, a material system interacting with a radiation field
is also examined. We study the propagation of an initial excitation signal
in a quantum dot that is prepared by an optical near field locally and in-
dividually, and describe it in terms of excitons in N two-level systems. In a
dilute limit, excitons are approximated as bosons, and a rigorous solution of
the Heisenberg equation shows that the dipole moments representing quan-
tum coherence between any two levels follow linear dynamics. In a nondilute
case, excitons at one site obey the fermion commutation relations, while ex-
citons at other sites satisfy the boson commutation relations; this results in
nonlinear equations of motion. We predict a coherent oscillation of all the
dipoles in the system (dipole ordering), and strong radiation from some of
the dipole-ordered states, which are close to Dicke’s superradiance [49].
150 A. Shojiguchi et al.
2 Model Hamiltonian
One of the most important features of optical near fields is the localization
property, where it is not suitable to use a broad spreading wave as a basis
function of quantization of the fields. Therefore, it is important to find a
152 A. Shojiguchi et al.
Fig. 4. Schematic diagram of a model system: N two-level QDs labelled by the site
number are closely located in a ring, and are interacting with localized photons
A Phenomenological Description of Optical Near Fields 153
QDs and its manipulation by the localized photons. The model Hamiltonian
of the system can be written as
H = Ha + Hb + Hint , (4)
where Ha represents localized photons, Hb describes excitons, and Hint rep-
resents the localized photon–exciton interaction. Each Hamiltonian can be
expressed as
N
Ha = εa†n an + V (a†n+1 an + a†n an+1 ) , (5a)
n=1
N
Hb = E b†n bn , (5b)
n=1
N
Hint = U (a†n bn + b†n an ) , (5c)
n=1
where n indicates the site number, and an (a†n ) and bn (b†n ) represent annihi-
lation (creation) operators of a localized photon and an exciton, respectively.
The periodic boundary condition requires that the (N +1)-th site corresponds
to the first site. The constant energies of the localized photons and excitons
are represented by ε = ω and E = Ω, respectively. The hopping energy
of the localized photons is represented as V = v, and U = g gives the
strength of the conventional dipolar coupling between a localized photon and
an exciton in the rotating-wave approximation.
We apply the boson commutation relations to the localized photons as
[an , a†n ] = δnn , [an , an ] = [a†n , a†n ] = 0 . (6)
The creation operator of an exciton, b†n ,
can be written in terms of the anni-
hilation operator of a valence electron of the n-th site, cn0 , and the creation
operator of a conducting electron of n-th site, c†nf , as b†n = c†nf cn0 . In the
same way, the annihilation operator of an exciton, bn , can be written as
bn = c†n0 cnf . Assuming that only one exciton is generated in a quantum dot
site, we obtain the relation as Nf + N0 = c†nf cnf + c†n0 cn0 = 1. With the help
of these expressions and the fermi commutation relations for electrons, the
commutation relations for excitons [72, 73] are derived as
[bn , b†n ] = δnn (1 − 2b†n bn ) , (7)
which shows that excitons behave as fermions at intrasite locations and as
bosons at intersite ones. It follows from (7) that excitons are approximated
as bosons in a dilute limit of the exciton density, b†n bn ≡ Nn 1.
With the boson approximation, the Heisenberg equation of motion can
be rigorously solved. First, we obtain a rigorous solution of the Heisenberg
equation for bosonic excitons to investigate the dynamical properties. Then
we investigate the dynamics of fermionic excitons, solving the Heisenberg
equation both perturbatively and numerically.
154 A. Shojiguchi et al.
It follows from this expression that when the expectation value of the number
of excitons is small as b†n bn N , the second term on the right-hand side
of (10) can be neglected, and excitons are approximated as bosons as
where the abbreviations e = ε + 2V cos k and X † (k) = (A†k , Bk† ) are used.
The matrix kernel of the Hamiltonian, M , is diagonalized as
T λ+ 0
S MS = , (17)
0 λ−
by an orthogonal matrix S as
⎛ ⎞
L+K L−K
S=⎝ 2L 2L ⎠ . (18)
L−K L+K
2L − 2L
where Pn (t) = bn (t) + b†n (t) denotes the dipole operator that describes the
electric dipole moment of each quantum dot. Moreover, two other indepen-
dent operators are defined as
Vn (t) = i(bn (t) − b†n (t)) , Wn (t) = b†n (t)bn (t) − bn (t)b†n (t) , (27)
where the latter represents the population difference of the system. The unit
operator and these three operators form a basis for any two-level system,
A Phenomenological Description of Optical Near Fields 157
Fig. 5. Time evolution of the dipole-moment distribution when the total number
of sites is eight. The parameters E = 2, ε = 1, and V = 1 are used in the
calculation. The vertical and horizontal axes represent the dipole amplitude and
the site number, respectively. The dipole moment is initially set only at the first
site
can obtain two copies of coherence at the farthest pair of sites as illustrated
in Fig. 7b.
Moreover, it might be possible to use the system as a nanophotonic device
of transporting or splitting quantum coherence.
A Phenomenological Description of Optical Near Fields 159
Fig. 6. Time evolution of the dipole moment of (a) the first site and (b) the fifth
site. When the dipoles are large in (a), corresponding dipoles in (b) are small, and
vice versa
Noticing that the Hamiltonian for localized photons Ha given by (5a) can be
written in a quadratic form as
⎛ ⎞
⎛ ⎞ a1
ε V V ⎜ ⎟
⎜V ε V ⎟⎜ ⎟
⎜ ⎟ ⎜ .. ⎟
⎜ . . . . . . ⎟ ⎜ . ⎟
Ha = (a†1 , · · · , a†N ) ⎜
⎜
⎟⎜
⎟⎜ . ⎟
⎟ (31)
⎜ . . . . . . ⎟ ⎜ .. ⎟
⎝ V ε V ⎠⎜ ⎝
⎟
⎠
V V ε
aN
†
≡ an Rnm am , (32)
nm
where the notation G0 (t) = eiL0 t is used. Substituting (40) into (38) and
using (36), we can obtain the time evolution of the exciton operator bn (t) =
G(2) (t)bn . Suppose that localized photons are initially in the vacuum and
Vn = 0, then the expectation value of bn (t) is expressed as
N
bn (t) = e−iΩt bn + g 2 Γj (t)
j=1
N
× Pnj Pmj (1 − δnm )Wn bm − 2
Pnj bn , (41)
m=1
is used, and cj (t) and dj (t) represent the real and imaginary parts of Γj (t),
respectively. The expression given by (41) and its Hermitian conjugate pro-
vide a time evolution of the expectation value of the dipole at an arbitrary
QD site n as
⎛ ⎞
Pn (t) = Pn ⎝cos Ωt − g 2 2 ⎠
cj (t)Pnj
j
+g 2 cj (t)Pnj Pmj Pm Wn . (43)
j m=n
In the next section we discuss the dynamical properties of the system, on the
basis of the solution given by (44).
Flip Hypothesis
If the sign of the initial population difference of the n-th site, Wn (0), is
inverted, then the direction of the dipole moment of the n-th site at arbitrary
time, Pn (t), is flipped.
Since the dipoles at sites 4, 5 and 6 in Fig. 8 are directed opposite to the
others, we invert the sign of Wn for n = 4, 5, 6. Figure 9 shows the result
that the direction of the dipoles at sites 4, 5, and 6 are flipped, and that all
the dipoles oscillate with the same phase but with different amplitudes. Thus
we observe that the system is transferred from a locally excited state to a
coherently oscillating state of the dipoles, in other words, to a dipole-ordered
state. In the ordered state the total dipole is N times larger than each single
dipole. Moreover, since radiation from an oscillating dipole is proportional to
164 A. Shojiguchi et al.
the square of the dipole moment, we can expect a high intensity of radiation
from the dipole-ordered state. The radiation property of the system will be
discussed in Sect. 7.1.
It is well known that the electric interaction between molecules with no dipole
moments is weak when they are separated in a macroscopic distance, and
that the interaction becomes strong when they are very close to each other.
Taking this into account, we investigate the dynamics of the system driven
by localized photons, from an initial state that alternating dipoles are set
and thus the total dipole of the system vanishes for an even number of sites
(see Fig. 10a). Such a dipole-forbidden state can be manipulated by localized
photons, not by propagating far fields.
As an initial condition we set Pn = (−1)n , and for simplicity, localized
photons are in the vacuum, Vn = 0, and Wn = 0. The system parameters
E = 2, ε = 1, and V = 1 are used as before. Figure 11 shows the time
evolution of such a dipole-forbidden state.
A Phenomenological Description of Optical Near Fields 165
√
and Wn = 1/ 2 for all n, the result is the same as shown in Fig. 11 except
for its amplitude of the oscillation.
Next, we manipulate the distribution of the population difference nonuni-
formly, so that the signs are set as opposite to those of the corresponding
dipoles:
√
Pn = (−, +, −, +, −, +, −, +)/ 2 ,
√
Wn = (+, −, +, −, +, −, +, −)/ 2 .
Figure 12 presents the result that a dipole-forbidden state as shown in
Fig. 10a is converted to a dipole-ordered state as illustrated in Fig. 10b.
Since the dipole ordering has occurred by manipulating the initial distri-
bution of the population differences, the result can be interpreted by the “flip
hypothesis” proposed in Sect. 4.2. As a collective oscillation of the dipoles oc-
curred, the system evolved from a nonradiative state to a radiative state
through the localized photon–exciton interaction.
A Phenomenological Description of Optical Near Fields 167
Here we should note that the second-order perturbative solutions break the
unitarity, and that the long-time behaviors of the dynamics given by the
solutions are not necessarily the same as the exact ones. Thus we have to
return to the Heisenberg equations, (30a) and (30b), to evaluate the time
evolution of physical observables. The Heisenberg equations are first solved
by neglecting the quantum correlations between excitons and localized pho-
tons such as Wn yn = Wn yn , and later the quantum correlations are
estimated in Sect. 8.2. In this semiclassical approximation, (30a) and (30b)
can be converted to the following coupled differential equations
⎧
⎨ ∂t Pn = −ΩVn + gWn yn ,
∂t Vn = ΩPn − gWn xn , (46)
⎩
∂t Wn = g (Vn xn − Pn yn ) ,
∂t xn = −ωyn − v (yn−1 + yn+1) − gVn ,
(47)
∂t yn = ωxn + v xn−1 + xn+1 + gPn ,
where the notations xn = an + a†n , yn = i(an − a†n ), and v = V / are used.
To check the reliability of the semiclassical approximation, we again inves-
tigate the system examined in Sect. 4.2. Figure 13 shows the time evolution
of the dipole distribution obtained from the semiclassical approach for the
system corresponding to Fig. 8.
Comparing Fig. 13 with Fig. 8, we find that both profiles are the same, and
that the amplitude of each dipole obtained from the Heisenberg equations is
less than 1, which means that the unitarity of the time evolution is conserved.
In Fig. 14, we show the time evolution of the dipole distribution obtained
from the semiclassical approach when the initial population differences of
sites 4, 5, and 6 are inverted. It follows from the figure that each dipole is
ordered as the flip hypothesis predicts, and that the flip hypothesis is still
valid in the numerical solution of the semiclassical Heisenberg equation.
Numerical results obtained from (46) and (47) with a semiclassical approx-
imation predict the dipole dynamics with qualitative similarities as shown
from the perturbative solution (44). The flip hypothesis proposed in Sect. 4.2
is also verified by using the numerical solutions of (46) and (47), or those of
the Heisenberg equations without quantum correlations. These results show
that the approach discussed in this section, maintaining an advantage to con-
serve the unitarity, qualitatively describe the dipole dynamics of the system
in a similar way as the other two approaches. The validity of the approxi-
mation employed here, which we will examine in Sect. 8.2, is reported in the
Dicke model that superradiant phenomena can be described with a semiclas-
sical approximation that neglects quantum correlations among the atoms [57],
and that the contribution from the quantum correlations is of the order of
1/N for large N [58].
168 A. Shojiguchi et al.
Fig. 13. Time evolution of the dipole distribution for the system (E = 2, ε = 1, V =
1) obtained from the Heisenberg equation with the semiclassical approximation.
Initially all of the sites except the first site are in the ground state
Fig. 14. Time evolution of the dipole distribution for the system (E = 2, ε = 1, V =
1) obtained from the Heisenberg equation with the semiclassical approximation. The
initial population differences of sites 4, 5, and 6 are inverted
where the exciton energy E is denoted as Ω. With the help of the compo-
nents of the Bloch vector (Pn , Vn , Wn ) the effective Hamiltonian (51a) can
170 A. Shojiguchi et al.
be rewritten as
1 + Wn
Heff = (Ω + ∆Ωn ) + ∆Ωnm (Pn Pm + Vn Vm ) , (52)
n
2 4 n
m=n
Using the effective Hamiltonian Heff , we can obtain the time evolution of the
dipole moment at site n in a similar way as discussed in Sect. 4.1:
Pn (t) = Pn cos[(Ω + ∆Ωn )t] + ∆Ωnm tWn Pm sin(Ωt) . (54)
m=n
∆Ωn+1,m+1 = ∆Ωn,m .
This is due to the fact that the periodic boundary condition is imposed on the
system, and that the Hamiltonian given in (4) has translational invariance.
Owing to the symmetric structure as
the dipole distribution of the system initiated by the dipole excitation at site
2 can be discussed in terms of ∆Ω1n that governs the problem driven by
the site-one excitation. In general the problem is reduced to examining the
structure of the coefficients ∆Ω1n that depend on the material parameters
Ω and λj , or E, ε, and V . Since the parameter sets crossing the resonant
points (Ω = λj ) change the pattern of the dipole distribution in the system,
we can classify quasisteady states of the dipole distribution into four groups
in Table 1.
The system with E = 2, ε = 1, and V = 1, which has been investigated
in Sect. 4, belongs to the first group. Since the coefficients ∆Ω1n in the first
group have different signs at sites 4, 5, and 6 from the others, the quasisteady
state of the dipole distribution that is predicted by the effective Hamiltonian
is the one described by the second-order perturbation solution (see (44) and
Fig. 8). The third and fourth groups are, in particular, interesting among
the four groups of the quasisteady states. The third group corresponds to a
“ferromagnetic” system, in which all the electric dipole moments are aligned
in the same direction, while the fourth group corresponds to an “antiferro-
magnetic” system, where the dipoles are aligned to have alternating signs.
Using numerical results of (46) and (47), we examine the dynamics of the
third and fourth groups whether they behave as predicted by the effective
Hamiltonian: “ferromagnetic,” or “antiferromagnetic.” The parameter set of
E = 3.01, ε = 1, and V = 1 gives
and thus the system belongs to the third group. The result is shown in Fig. 15,
which predicts the “ferromagnetic” behavior as expected.
When the parameter set of E = 1, ε = 3.01, and V = 1 is chosen, the
system corresponds to the fourth group. Figure 16 presents the time evolution
172 A. Shojiguchi et al.
Fig. 15. Time evolution of the dipole distribution of the system belonging to the
third group. The initial conditions are set as Pn = δn1 and Wn = 1 − δn1 . All
the electric dipole moments are aligned in the same direction, and a dipole-ordered
state appears
Fig. 16. Time evolution of the dipole distribution of the system belonging to the
fourth group. The initial conditions are set as Pn = δn1 and Wn = 1 − δn1 .
The dipole distribution similar to the “antiferromagnetic” state is shown, where
the dipoles alternatingly change the directions
Fig. 18. Time evolution of the dipole distribution of the system belong-
ing to the fourth group. The initial population differences are set as Wn =
(0, −, +, −, +, −, +, −). According to the flip hypothesis, the alternating dipoles
change their directions and are aligned to form a dipole-ordered state
Until now, the state of localized photons at each site was initially taken to be
a vacuum. Since localized photons in the Heisenberg equations have different
degrees of freedom, this section examines the temporal evolution given an
initial distribution of localized photons. The initial conditions take the state
A Phenomenological Description of Optical Near Fields 175
given below, in which the distribution of localized photons is given for site 1
only.
This state can be achieved by having the localized photons at site 1 take on
a coherent state in which the eigenvalue is a real number. The coherent state
is the eigenstate of the annihilation operator of a simple harmonic oscillator.
Therefore, the initial state satisfying the following condition needs to be
established at site 1:
Fig. 19. The temporal evolution of the dipole distribution of the system for the
first group in which the localized photonic states are initially in the coherent state
at the first site and the other sites are in the vacuum state and all the excitonic
states are in the totally inverted state
176 A. Shojiguchi et al.
Fig. 20. Quasisteady states arise from the initial condition in which there is lo-
calized photon excitation at a single site when the system parameters have various
values (, V, E). The results in (a), (b), and (c) correspond to the first, third, and
fourth groups in Table 1, respectively
6 Dicke’s Superradiance
In the preceding section, we have shown that electric transition dipole mo-
ments of the system are aligned to have a large value, depending on the
material parameters of the system. It is then expected that strong radiation
can be emitted due to a large total dipole moment, and it is intriguing to
examine radiation properties of the dipole-ordered states obtained from our
localized photon model. Before such an investigation, we will outline superra-
diance [49, 54–57, 60, 63] as a cooperative phenomenon, where radiation from
N two-level systems interacting with propagating light exceeds the sum of
the one individually emitted from the systems. More precisely, a radiation
pulse with a sharp width of 1/N and a high intensity of N 2 is emitted [49].
This kind of system is called the Dicke model, and was extensively studied
by the middle of the 1970s since the pioneering work of Dicke in 1954 [49].
In this section we first introduce the Dicke model and Dicke states, and
then discuss superradiant states, or superradiance, following Mandel and
Wolf [77]. In order to examine how such states are dynamically obtained, a
master equation is derived with the help of the projection operator method,
and the time evolution of radiation is studied. Numerical solutions of the
Dicke master equation are presented for a small system whose size is much
smaller than the wavelength of the radiation field. At the same time, ana-
lytic solutions are obtained within a semiclassical approximation to discuss
the importance of quantum correlation. Finally, we examine the effects of the
dipole–dipole interaction for an inhomogeneous system, as well as the effects
of reabsorption and re-emission of photons for a large system, both of which
violate the cooperative states and destroy some part of the coherence in the
system.
where Ω, b†n , and bn denote the excitation energy, creation and annihilation
operators of excitation in the n-th site, respectively. These operators can also
(n)
be defined by spin operators Ri (i = 1, 2, 3) and Bloch vectors (Pn , Vn , Wn )
as follows.
178 A. Shojiguchi et al.
1 1
(bn + b†n ) = Pn ,
(n)
R1 = (59a)
2 2
(n) i † 1
R2 = (bn − bn ) = Vn , (59b)
2 2
(n) 1 † † 1
R3 = (bn bn − bn bn ) = Wn . (59c)
2 2
(n)
The spin operators Ri satisfy the commutation relations
(n) (n ) (n)
[Ri , Rj ] = iδnn ijk Rk , (60)
with Levi-Civita’s symbol ijk . Introducing total spin operators for the N
two-level systems as
N
N
N
(n) (n) (n)
R1 = R1 , R2 = R2 , R3 = R3 , (61)
n=1 n=1 n=1
where the total spin operators also obey the commutation relations
similar to (60).
Let us consider a state in which N1 of the N two-level systems are in the
lower state |1 and N2 are in the excited state |2, where
1
N = N1 + N 2 , m= (N2 − N1 ) . (64)
2
Evidently m is a measure of the total population inversion, and it is an
eigenvalue of the spin operator R3 . Introducing the operator R2 = R12 + R22 +
R32 and using (63), we obtain the commutation relations
so that the collective operators Ri obey the same algebra as the angular mo-
mentum. Thus there exist the states |l, m that are simultaneous eigenstates
of both R3 and R2 as
1 1
R3 |l, m = m|l, m , for N ≥ m ≥ − N , (66a)
2 2
1
R2 |l, m = l(l + 1)|l, m , for N ≥ l ≥ |m| . (66b)
2
The quantum number l is called the cooperation number [49] and determines
the collectivity of cooperative phenomena. For the raising and lowering op-
erators as
A Phenomenological Description of Optical Near Fields 179
N
N
R+ = b†n = R1 + iR2 , R− = bn = R1 − iR2 , (67)
n=1 n=1
Superradiant States
In the following, let us assume that the system size is much smaller than the
wavelength of the radiation, and that the long-wave approximation is valid.
The interaction Hamiltonian between the two-level system and radiation in
the dipole coupling is given by [49]
where the electric field is evaluated at the origin owing to the long-wave
approximation, and the total electric dipole of the system is written as
N
µtot = µ (bn + b†n ) = µ(R+ + R− ) . (70)
n=1
Here we suppose a transition from an initial state |l, m|vac to a final state
|ψ, where |l, m and |vac represent a Dicke state and the vacuum of the
radiation field, respectively. Then the transition probability can be written
as
where E − (t) is the e−iωt -dependent part of the electric field E(t). By squaring
the absolute value of the transition probability and summing over all possible
final states, the probability of emission of radiation photons wemit is given
as [49, 77]
wemit = |ψ|µ · E(t)R− (t)|l, m|vac|2
ψ
= l, m|R+ (t)R− (t)|l, mvac|µ · E + (t)µ · E − (t)|vac
= R+ (t)R− (t)A = (l + m)(l − m + 1)A , (72)
where (68) is used for the final expression, and A = vac|µ · E + (t)µ ·
E − (t)|vac is identified as the Einstein A coefficient for each two-level system
like an atom or a quantum dot.
If all quantum dots are in the ground state |1, then m = −N/2, and
from (66b) we must have l = N/2. Hence, from (72) the radiation rate wemit
180 A. Shojiguchi et al.
is zero as expected. If all quantum dots are in the excited states |2, then
m = N/2 and we again have l = N/2. Hence wemit = N A, which is just
what one would expect from a group of N independently radiating quantum
dots. However, the situation is quite different if the initial state is not the
completely excited state. Let us consider a state in which half of the quantum
dots are excited and half are not, so that m = 0. Then we have
and l can have any value between 0 and N/2. The larger the value of l,
the larger is the collective rate of radiation of the system. In particular, if
l = N/2, then
1
wemit = N (N + 2)A , (74)
4
which maximizes the value of cooperative spontaneous radiation. Since the
maximum radiation probability is proportional to the square of the site
number N , it becomes very large for large N . This phenomenon is called
superradiance, and the state |l = N/2, m = 0 is called a superradiant
Dickestate. From the
definition of the total Bloch vector (P, V, W ) as
P = n Pn , V = n Vn , W = n Wn , it follows that the Bloch vec-
tor for the superradiant state is on the P –V plane and its length reaches a
maximum, as shown in Fig. 21.
In this section the Dicke master equation is derived for a small system to
discuss dynamical manifestation of superradiant states described in the pre-
ceding section.
According to Mandel and Wolf [77, 78], we outline the derivation of a gener-
alized master equation by projection operators to project out of the density
A Phenomenological Description of Optical Near Fields 181
operator for a total system that part on which our interest is focused [79–82].
The total system can be described in the interaction representation by the
Liouville–von Neumann equation for the density operator ρ(t) as
∂ρ(t) 1
= [HI (t), ρ(t)] = −iL(t)ρ(t) , (75)
∂t i
where HI (t) is the interaction Hamiltonian, and L · · · = [HI (t), · · · ]/ is the
Liouvillian operator, as before. We now introduce a time-independent pro-
jection operator P :
P2 = P , (76)
which is chosen so as to project out the most relevant part of ρ(t). Using the
complementary operator (1 − P ), we can express ρ(t) as
On multiplying both sides of (75) by P on the left and then using (77) to
substitute for ρ(t) on the right, we have
∂ρ(t)
P = −iP L(t)P ρ(t) − iP L(t)(1 − P )ρ(t) . (78)
∂t
Similarly, after multiplying (75) by (1 − P ) on the left, we obtain
∂ρ(t)
(1 − P ) = −i(1 − P )L(t)P ρ(t) − i(1 − P )L(t)(1 − P )ρ(t) . (79)
∂t
Formally integrating the first-order differential equation (79) and solving for
(1 − P )ρ(t), we find
t
(1 − P )ρ(t) = exp −i(1 − P ) L(t )dt (1 − P )ρ(0)
0
t t
−i exp −i(1 − P ) L(t )dt
0 τ
×(1 − P )L(τ )P ρ(τ )dτ . (80)
Substitution of (1 − P )ρ(t) into (78) gives
∂ρ(t)
P = −iP L(t)P ρ(t)
∂t
t
−iP L(t) exp −i(1 − P ) L(t )dt (1 − P )ρ(0)
0
t t
−P L(t) exp −i(1 − P ) L(t )dt
0 τ
×(1 − P )L(τ )P ρ(τ )dτ . (81)
182 A. Shojiguchi et al.
If we set
t
U (t, τ ) ≡ exp −i(1 − P ) L(t )dt , (82)
τ
Under these conditions, the first term on the right-hand side of (83),
P L(t)P ρ(t), can be written as
1
P L(t)P ρ(t) = |vacvac|TrF [HFint (t), |vacvac|ρA (t)]
⎛ ⎞
= |vacvac| ⎝ gkλ vac|a† |vacR− ei(ω−Ω)t ρA (t) − h.c.⎠
kλ
kλ
=0,
(89)
because
The rate of photon emission was obtained perturbatively in Sect. 6.1. In this
section, we examine the time evolution of this emission rate, using the Dicke
master equation (99). Since the rate of energy dissipation from the two-level
(atomic) system is proportional to the emission rate, −∂R3 /∂t should be
evaluated. Multiplying (99) by R3 on the right and taking the trace, one finds
∂
R3 = β{−R3 R+ R− + 2R+ R3 R− − R+ R− R3 }
∂t
−iγ{R3 R+ R− − R+ R− R3 } . (101)
In order to simplify the above equation, we use the following relations
[R+ , R3 ] = −R+ , (102)
1 †
R+ R− = R3 + N bn bm , (103)
2 n m=n
which follow from the collective operators defined in (67) and the commuta-
tion relations in (63). In addition, using (103), we can show that
A Phenomenological Description of Optical Near Fields 185
[R+ R− , R3 ] = [b†n bm , R3 ]
n m=n
= [b†n bm − b†n bm ] = 0 , (104)
n m=n
∂
R3 = 2β[R+ R3 R− − R3 R+ R− ] . (105)
∂t
With the help of (102), this leads to the final equation
∂
R3 = 2β[(R3 R+ − R+ )R− − R3 R+ R− ]
∂t
= −2βR+ R− . (106)
It follows from the above equation that the time evolution of R+ R− is
required to evaluate the rate of energy dissipation. Multiplying (99) by R+ R−
on the right and taking the trace, we then find
∂
R+ R− = −2β{R+ R− R+ R− − R+ R+ R− R− }
∂t
−iγ[R+ R− , R+ R− ]
= −4β{R+ R− − R+ R− R3 } , (107)
where we used the following relation to obtain the second line from the first
line
we then have
∂
R+ R− = −2βR+ R− (1 − R3 ) . (110)
∂t
Simultaneously solving (106) and (110), one can obtain the final result. Or
putting z = R3 , this leads to the following nonlinear equation
This allows us to evaluate the rate of energy dissipation, −∂ R3 /∂t = −ż,
and thus to obtain the rate of photon emission.
Under the following initial conditions
Fig. 22. Time evolution of the rate of photon emission that is numerically obtained
from (111) with the initial condition of a totally inverted state. The vertical axis is
represented in units of the Einstein A coefficient A = 2β for β = 0.05. The dotted,
dashed, and solid curves show the results for the site number N = 12, 15, and 20,
respectively
which correspond to l(0) = N/2 and m(0) = N/2 in terms of Dicke states,
we solve (111) numerically. The result is shown in Fig. 22 for the case of
β = 0.05 and the site number N = 12, 15, and 20. It follows from Fig. 22
that the peak value of the rate of emission becomes higher and the width
gets narrower, as the site number becomes larger. As will be shown in the
next subsection, a pulse of radiation whose height is proportional to N 2 and
width is proportional to 1/N is emitted from the two-level system prepared
in the initial condition.
Fig. 23. Time evolution of the rate of photon emission plotted according to the
analytic solution obtained with the semiclassical approximation. The vertical axis
is represented in units of the Einstein A coefficient A = 2β. The parameters are set
to β = 0.05, t0 = 3, and N = 20, corresponding to Fig. 22
where the electric field operator E is divided into the positive and negative
frequency parts E + and E − defined as
ck
E (r n ) = i
+
eλ akλ eik·r n , (123a)
20 V
kλ
E − (r n ) = (E + (r n ))† . (123b)
The creation and annihilation operators a† and akλ act on the mode of the
kλ
radiation field with wavevector k and polarization eλ , and obey the boson
commutation relation. An arbitrary quantization volume is denoted as V .
The dipole operator is similarly defined as
P (r n ) = µ(b†n + bn ) . (124)
A Phenomenological Description of Optical Near Fields 189
Using this interaction Hamiltonian, we can obtain the following master equa-
tion in the Born–Markov approximation [62, 83]
dρ
= −i Ωnm [b†n bm , ρ] − γnm [b†n bm ρ + ρb†n bm − 2bm ρb†n ] , (125)
dt nm
n=m
with
µ2 Ω 3
γnm = F (Ωrnm /c) , (126a)
16π 2 c3 0
∞
µ2 P
Ωnm = − dkk 3 F (krnm ) , (126b)
3
16π 0 0 k − Ω/c
F (krnm ) = dΩk (e · )2 eik·r nm , (126c)
Fig. 24. Arrangement of three two-level systems. Each two-level system is config-
ured at the apex of an isosceles triangle
from which it follows that R2 is generally not conserved, but is a conserved
quantity only if the site number N is equal to 2, or only if Ωnm is constant for
all suffixes [62, 69]. Therefore the state is not expressed by Dicke states, and
quantum coherence of the two-level system is destroyed by the dipole–dipole
interactions.
Now we consider how the dipole–dipole interactions affect superradiance.
Let the site number be three and be arranged as shown in Fig. 24. For sim-
plicity, we assume Ω12 = Ω13 = Ω23 and γmn = γ.
The following eight bases are employed to take matrix elements of (125)
and to solve the simultaneous differential equations.
where the sign + indicates that the upper level is occupied, while the sign −
means that the lower level is occupied. Under the initial condition of + +
+|ρ(0)| + ++ = 1, we obtain the matrix elements of the density operator
ρ(t) at time t to express the radiation intensity
I(t) = b†n bm = R+ R− . (132)
n,m
Using the Pauli matrices σµ (µ = 0, 1, 2, 3), we can write nonzero components
of the density matrices as µ xµ σµ . The vector x = (x1 , x2 , x3 ) executes fast
√
precession about a “torque” Γ = 2 2Ω12 , 0, Ω23 (see Fig. 25). When we
introduce the rotation R β, 2̂ that rotates about the 2 axis through an angle
√
2 2Ω12
tan β = , (133)
Ω23
then the vector Γ is transformed as R β, 2̂ Γ = (0, 0, Γ ) with
2
√
Γ = 2 .
2 2Ω12 + Ω23
A Phenomenological Description of Optical Near Fields 191
Fig. 25. Relation between the angles β and χ. In the case of pure superradiance
(Ω12 = Ω23 ), the torque Γ is given by Γ0
The solutions of the differential equations for the density matrices can be
obtained from a similar rotation of x with cos χ defined as [62]
1
cos (χ + β) = . (134)
3
When cos χ = 1 or χ = 0, Ω12 = Ω23 is satisfied and the Dicke superradiant
states become exact eigenstates, while cos χ = 0 gives orthogonal states to
the superradiant states. This means that the deviation of χ from zero shows
the inhomogeneity of the system, or a measure of the coherence properties.
According to Coffey and Friedberg [62], we present the radiation intensity
I(t) in Fig. 26, which shows the peak and tail of I(t) become reduced and
longer, respectively, as cos χ varies from one to zero corresponding to the
sequence of a, c, d, and b in Fig. 26.
Tokihiro et al. [69] also examined the effects of the dipole–dipole interac-
tion on superradiance, using a linear excitonic system, in which excitons can
hop from one site to its nearest neighbors via the dipole–dipole interaction.
They used the following master equation in the interaction representation
dρ
= −iJ [b†n bn+1 + b†n+1 bn , ρ]
dt n
−γ[R+ R− ρ + ρR+ R− − 2R− ρR+ ] , (135)
where the hopping energy of an exciton and the radiation decay rate are
denoted as
µ2 ( · a)2
2J = 3
1−3 , (136)
4π0 a a2
µ2 Ω 3
γ= , (137)
3πc3 0
respectively. The first term in (135) corresponds to the Hamiltonian of the
XY-model in a spin system whose exact solutions are well known [84]. Tak-
192 A. Shojiguchi et al.
Fig. 26. Normalized radiation intensity for three two-level systems with different
Ω12 /Ω23 : (a) Ω12 = Ω23 , cos χ = 1, corresponding to pure superradiance, (b)
cos χ = 0, (c) Ω12 = 8Ω23 , cos χ = 0.96, (d) Ω12 = Ω23 /8, cos χ = 0.63, (e)
I = 3e−γt , where three two-level systems radiate independently. The vertical axis
is represented in units of the Einstein A coefficient while the horizontal axis is
represented by γt
ated pulses [63, 65]. We will analyze this effect semiclassically, according to
Banifacio et al. [60, 63, 85], and comment on experimental observations.
As samples have a volume with linear dimensions much larger than the wave-
length of the radiation, some quantitative differences appear in the dynamics
of the system evolution due to the effects connected with the propagation of
the electromagnetic field along the emitting sample. According to Bonifacio
et al., we use a laser model to describe equations of motion as
where LAF denotes the Liouvillian operator, and ΛF is the dissipation term
that represents photons coming away from the cavity. They are defined as
follows:
1
LAF ρ = [HAF , ρ] , (139a)
HAF = g(aR+ + a† R− ) , (139b)
† †
ΛFX = γ{[aρ, a ] + [a, ρa ]} , (139c)
where the constants (g and γ) and the operators have their usual meanings,
as before. If we multiply (138) on the left by a, R+ , and R3 , respectively, and
if we take the trace, we obtain equations for α(t) = ia(t), r1 (t) = R+ (t),
and r3 (t) = R3 (t) as
ṙ1 = 2gαr3 , (140a)
ṙ3 = −2gαr1 , (140b)
α̇ = gr1 − γα , (140c)
where we have used the semiclassical approximation to neglect quantum cor-
relations. Introducing the polar angle θ of a Bloch vector as
N N
r1 (t) = sin θ , r3 (t) = cos θ , (141)
2 2
one finds the following differential equation that shows the damped harmonic
oscillation [60, 63]
and the equations for conservation of the magnitude of the Bloch vector R2
and energy of the system
∂t (r12 + r32 ) = 0 , (143a)
∂t (α2 + r3 ) = −2γα2 . (143b)
194 A. Shojiguchi et al.
Fig. 27. Time evolution of the rate of radiation from a large sample. The vertical
axis is represented in units of the Einstein A coefficient as before. The parameters
used are γ = 0.1, g = 1/14, and N = 20
It follows from (143b) that the dissipation energy from the system per unit
time is 2γα2 , and it corresponds to radiation energy from the system. Thus
we can write radiation intensity I(t) as
γ 2
I(t) = 2γα2 = θ̇ . (144)
2g 2
N g2
ż = − (1 − z 2 ) , (146)
γ
which is the same as (115), the Dicke master equation for a small system
obtained semiclassically.
At the end of this subsection, we briefly comment on experimental obser-
vation of superradiance. Skribanovitz et al. [64] experimentally showed large
A Phenomenological Description of Optical Near Fields 195
which indicates that there are two elements mainly contributing to the ra-
diation factor; one is the collectiveness of the system measured by the co-
operation number l, and the other is the total energy of the quantum-dot
system given by R3 = m. Since the dynamics of the system is driven in our
system by localized photons, the evolution of the radiation factor is also de-
scribed in terms of l(t) and m(t) that are developed according to the localized
photon–exciton interaction.
Numerical results of (46) and (47) are shown in Figs. 28 and 29. The
upper parts in Fig. 28 show the time evolution of the radiation factors, while
the lower parts illustrate the dipole distribution when each radiation factor
has a maximal value. The total dipole in Fig. 28a is smaller than those in
Figs. 28b and 28c, and the peak value of the radiation factor is also smaller. It
follows from the figure that the radiation factor increases as the total dipole
becomes larger. The peak values of the radiation factor in Figs. 28b and 28c,
15 and 14, correspond to the value for the Dicke’s superradiance, which is
obtained as 20 for N = 8 from (74). We thus expect that quasisteady states
shown in Figs. 28b and 28c are close to the superradiant states, and that the
total Bloch vectors for such states are on the P − V plane.
196 A. Shojiguchi et al.
Fig. 28. Time evolution of the radiation factor R+ R− (upper) and the dipole
distribution at time indicated by the arrows when the radiation factors take max-
imal values (lower) for (a) a system belonging to the first group with an initial
population difference Wn = (0 + + − − − ++), (b) a system belonging to the
third group with an initial population difference Wn = (0 + + + + + ++), and
(c) a system belonging to the fourth group with an initial population difference
Wn = (0 − + − + − +−)
Fig. 29. Time evolution of the radiation factor (left) and the polar angle of the
total Bloch vector (right) for a system belonging to the third group. The polar angle
is measured from the Wn = (− − − − − − −−) axis
In order to check whether the total Bloch vectors for the states belonging
to the third group in Table 1 (see Fig. 28b) are on the P − V plane or not, we
examine the time evolution of the polar angles of the Bloch vectors as well
as the radiation factors. As shown in Fig. 29, the polar angle of the Bloch
vector takes 90 degrees, that is, on the P − V plane when the radiation factor
has a maximal value. Therefore we conclude that the system belonging to
the third group in Table 1 is in transition to a quasisteady state close to
A Phenomenological Description of Optical Near Fields 197
Fig. 30. The dependence of the width and the peak point of superradiant pulse
on localized photon–exciton couplings g: double logarithmic plot (left) and the
definitions of the pulse width δt, the pulse peak points TA , and the bottom point
of pulse valley TB (right)
the superradiant Dicke state, judging from its large radiation factor and the
polar angle of the total Bloch vector.
In order to see the mechanism of the superradiance phenomena in our
model we shall investigate the dependence of the localized photon–exciton
couplings g on the width of superradiant pulse. The result is shown in Fig. 30.
From this plot we find a relation as follows
1
T ∝ . (148)
g
Recalling the mechanism of the oscillation of superradiant pulse in the Dicke
model explained by Bonifasio (see (142)) we know that this result means
the oscillation of superradiant pulses in our model is caused by emission and
absorption of localized photons.
Figures 28b and 28c indicate that multiple peaks appear in the radiation,
or that multiple pulses are emitted from the system. One may think, as a
possible origin, that such a phenomenon stems from the recurrence inherent
in an isolated system. However, such multiple pulses may survive even if the
system becomes dissipative, which will be examined in detail in the next
section.
where HQDeff is the effective Hamiltonian given by (51a) that describes the N
two-level quantum-dot system interacting with localized photons. The Hamil-
tonians HF and HFint describe the free-radiation field and the exciton–free
photon interaction, respectively. Explicit Hamiltonians are written as
HF = ωk a† ak,λ , (150)
k,λ
k,λ
i(ω −Ω)t
HFint (t) = gk,λ a† R− e k + h.c. , (151)
k,λ
k,λ
where creation and annihilation operators of a free photon with wavevector k,
polarization λ, and frequency ωk are denoted as a† and ak,λ , respectively.
k,λ
The coupling constant between the free photon and exciton is given as
iΩ
gk,λ = √ µ · ∗k,λ , (152)
V 2ωk 0
P · · · = |00|TrF · · · , (154)
with the vacuum of the radiation field |0. Moving to the interaction rep-
resentation, we obtain equations of motion for the density operator ρ̃(t) as
follows:
∂ ρ̃
P = −iPL2 (t)P ρ̃(t) − iPL2 (t)U (t, 0)(1 − P)ρ̃(0)
∂t
t
−PL2 (t) dτ U (t, τ )(1 − P)L2 (τ )P ρ̃(τ ) , (155)
0
where L2 (t) is the Liouville operator associated with H2 (t), and the operator
U (t, τ ) is defined as
t
U (t, τ ) = exp −i(1 − P) L2 (t )dt . (156)
τ
A Phenomenological Description of Optical Near Fields 199
where LQDeff and LFint are the Liouville operators associated with the Hamil-
tonians HQDeff and HFint , respectively. Moreover, using the Born approxima-
tion that neglects the exciton operators of higher than the second order [77],
we approximate
ρ̃A = ρn , (162)
n
and then solve (160). Noticing that the dynamics governed by the original
Hamiltonian (4) is not rigorously identical to the one described by the effec-
tive Hamiltonian (51a), we use an isolated system described by the effective
200 A. Shojiguchi et al.
Fig. 31. Time evolution of the radiation factor for dissipative systems (upper) and
isolated systems (lower). We consider the following three cases for both systems:
(a) a “ferromagnetic” case, (b) an “antiferromagnetic” case that is turned into a
dipole-ordered state after manipulating the initial population differences, and (c) a
“dipole-forbidden” case discussed in Sect. 4.3. The parameters β = 0.05 and γ = 0
are used
It is well known that superradiance in the Dicke model [49] occurs from a
state where all excited states of all sites are occupied. On the other hand,
the semiclassical approach discussed above cannot predict the occurrence of
superradiance of the system when the total dipole of the system is zero as
an initial condition [58]. This means that quantum fluctuations and corre-
lations should be properly included so as to correctly describe the radiation
properties of a system with no initial dipoles, and that the semiclassical ap-
proximation is not appropriate in this case. Thus we numerically solve the
master equation (160), taking quantum correlations into account, and we
investigate the radiation properties of the dissipative system. In particular,
we compare the results obtained from our model, i.e., the localized photon
model with those obtained from the Dicke model, for which the first term of
the right-hand side is dropped from (160). Some of such results are shown in
Figs. 32 and 33.
Figure 32a shows the time evolution of the radiation factor for the case
that all the populations are completely in the excited states and there are
no dipoles as an initial condition. The solid curves are the results for our
model, while the dashed curves represent the results for the Dicke model. It
follows from the figure that a single superradiant pulse is emitted in both
models, but that the peak value of the superradiant pulse is reduced, while
the tail is extended in our model. In Fig. 32b, we present the result obtained
from the initial√condition as a zigzag profile of the dipole distribution √ of
Pn = (−1)n / 2 and flat population differences of Wn = 1/ 2. The
system corresponds to a dipole-forbidden state as shown in Fig. 11, where the
total dipole is always zero. The radiation profiles obtained for both models are
qualitatively, same as found as in Fig. 32a. The time evolution of the radiation
factor is illustrated in Fig.
√ 32c for the case that initially the alternating dipole
is set as Pn = (−1)n / 2 and the sign of the population difference in each
QD is set opposite to that of the corresponding dipole as Wn = −Pn .
Owing to the flip hypothesis, a dipole-ordered state emerges in this case,
and a distinct difference between the two models is observed: our model
(solid curve) shows superradiance, while the Dicke model (dashed curve) does
not. Since we can infer that the difference stems from the dipole-ordering
phenomenon, we further examine this case in order to clarify the relation
between the radiation factor and the total dipole of the system.
Figure 33 presents the time evolution of the radiation factor (upper) and
the total dipole (lower) by changing the exciton–localized-photon coupling
constant g as g = 0.5, 0.8, and 1.2. It is found that the frequency of the
total dipole increases as the coupling constant becomes strong. As a result,
the oscillation frequency of the radiation factor increases, and thus multiple
pulses are emitted superradiantly from the dissipative system. Therefore the
difference between the two models, as we expected, originates from the occur-
202 A. Shojiguchi et al.
Fig. 32. Time evolution of the radiation factor obtained with quantum correlations.
The solid curves are the results for our model, while the dashed curves represent the
results for the Dicke model. The dissipative system is assumed, and the parameters
N = 4, β = 0.05, and γ = 0 are used. In addition, the following
√ initial conditions
√
are used: (a) Pn = 1 and Wn = 1, (b) Pn = (−1)n / 2 and Wn = 1/ 2
which correspond
√ to a dipole-forbidden state as shown in Fig. 11, and (c) Pn =
(−1)n / 2 = − Wn , which corresponds to a dipole-ordered state as shown in
Fig. 12
Fig. 33. Time evolution of the radiation factor (upper) and the total dipole moment
(lower) obtained with quantum correlations. The oscillating curves are the results
for our model, while the monotonically decreasing curves represent the results for
the Dicke model. The system is assumed to be dissipative and in a dipole-ordered
state evolved from a dipole-forbidden state, i.e., the same as in Fig. 32c. The initial
conditions are also the same as in Fig. 32c, except for the exciton–localized-photon
coupling (a) g = 0.5, (b) g = 0.8, and (c) g = 1.2
4.5 1.4
5.0
Radiation factor
Radiation factor
1.2
Radiation factor
Fig. 34. Radiation profile obtained from three different methods. The system and
the initial conditions to be considered are the same as in Fig. 33. The result is
obtained (a) for the isolated system by the Liouvillian dynamics of the effective
Hamiltonian with the semiclassical approximation, (b) for the dissipative system
by the master equation (160) with the semiclassical approximation, and (c) for the
dissipative system by the master equation (160) with quantum correlations. The
parameters N = 4, β = 0.05, and γ = 0 are used
that the semiclassical approach can describe qualitatively the radiation prop-
erties of both isolated and dissipative systems when the total dipole of the
system is not zero. The strong radiation coming from the dipole-ordering phe-
nomenon, or the nonlinearity and the collective phenomena of the dynamics
of the system considered in this chapter can be qualitatively predicted by the
semiclassical approach.
204 A. Shojiguchi et al.
Fig. 35. The QDs interacting with localized photons are configured in a line. At
each end of the line, we set different boundary conditions
Fig. 37. Behavior 2: the dynamics of the distribution of localized photons of an open
system with reservoirs of localized photons in the coherent states at the zeroth and
ninth site. In this case, we observe the growth of a structure similar to a standing
wave
to the right again. Consequently, it travels back and forth. In other words,
linear behaviors like transmission of a wave can be observed in the temporal
evolution of localized photons despite the nonlinear nature of the system’s
dynamics. In Fig. 37, there is a clear division between the sites where local-
ized photons accumulate (peaks) and the sites where localized photons are
rare (nodes). Localized photons continue to accumulate in the system and no
stationary state emerges. This can be viewed as formation
of a standing wave.
Examination of the total number of localized photons n an (t)† an (t) in
the system yields the results shown in Fig. 38.
Figure 38 clearly distinguishes between the two types of system when
looking at the number of photons. In the system shown in Fig. 36, in which
the transportation of localized photons follows a linear wave motion, the
total number of photons is approximately constant. This means that localized
photons in the system are flowing through the system. Therefore we call this
through flowing mode. In the system shown in Fig. 37, the total number
of localized photons increases monotonically and localized photons continue
to accumulate in the system. We call this phenomenon the photon storage
A Phenomenological Description of Optical Near Fields 207
Fig. 38. Temporal evolution of the total localized photon number n an (t)† an (t)
(feet). (a) and (b) correspond to the systems in Fig. 36 and Fig. 37, respectively
mode. The fact that the relaxation phenomenon does not occur implies that
either photons can accumulate in the system, or that the system relaxes very
slowly and our calculations are insufficient for it to reach equilibrium.
Following the discussion given in Sect. 7 and using (147), we then cal-
culated the radiation intensity for the two systems, as shown in Fig. 39. It
follows from Fig. 39 that no cooperative effect can be seen in the system
shown in Fig. 36, while a cooperative effect (superradiance) and a time-
irreversible phenomenon (dacaying profile of radiation) is observed in the
system shown in Fig. 37. Here the system parameters used in Fig. 37 are
same as those used in Fig. 28b, but it should be noted that a superradiant
radiation intensity decays time-irreversibly in the open system as assumed in
Fig. 37 while a recursive superradiance occurs in the closed system consid-
Fig. 39. Temporal evolution of the radiation intensity of the open systems. (a) and
(b) correspond to the systems in Fig. 36 and Fig. 37, respectively. Strong radiation
occurs in the system that is accumulating the localized photons and not in the
system where there is no accumulation
208 A. Shojiguchi et al.
Fig. 40. Temporal evolution of the total number of localized photons in the system
corresponding to Fig. 38b when the coupling constant g equals 0.65, 0.67, and 0.69.
We observe a drastic change of the number of localized photons between g = 0.65
and g = 0.67. The difference amounts to a single digit
Linear System
The Heisenberg equations of the system become the following when an exciton
is approximated by a boson:
⎧
⎪
⎪ ∂t Pn = −ΩVn − gyn ,
⎨
∂t Vn = ΩPn + gxn ,
(164)
⎪
⎪ ∂t xn = −ωyn − v(yn−1 + yn+1 ) − gVn ,
⎩
∂t yn = ωxn + v(xn−1 + xn+1 ) + gPn .
These are clearly linear equations. Since the equation of motion of the num-
ber of localized photons an (t)† an (t) is the same as in the case of (163) dis-
cussed above, which is obviously linear, in the boson-approximated system
210 A. Shojiguchi et al.
Fig. 41. Dependence of the total localized photon number on the coupling con-
stant g at t = 100. (a) and (b) correspond to the systems in Fig. 36 and Fig. 37,
respectively
Fig. 42. The dependence of the coupling constant g of the total localized photon
number at t = 100 for the system in which there is boson-approximation of excitons.
(a) and (b) take the same system parameters as the systems in Fig. 36 and Fig. 37,
respectively
all the equations of motion are linear. Solving these equations, Fig. 42 shows
the dependency of the total number of localized photons on the coupling
constant g. Although random behavior does not appear in this case, resonat-
ing coupling constant dependence is found. The meaning of the resonating
coupling constants is not yet clear; however, at least there is no disorderly,
unpredictable behavior. From these observations, we conclude that the dis-
orderliness of Fig. 41 is a chaotic behavior originating from the nonlinearity
of the dynamics.
Here, one must be cautious, because a system consisting of eight sites has
many modes, and it is possible that an apparently disorderly dependence has
A Phenomenological Description of Optical Near Fields 211
emerged via the competition of many modes (for example, the superposition
of many types of quasiperiodic behavior). Quasiperiodic behavior does not
have a period either, so chaos does not exist simply because of a lack of
periodicity. In the next section, the phenomenon is examined by simplifying
the problem by reducing the degrees of freedom to the minimal value, while
considering the above.
Fig. 43. The dependence of the coupling constant g on the total localized photon-
number of the two-site open system. (A) is the plot for the range from g = 0 to
4.1. (a)–(d) show the chaotic regions in (A). (a) shows the range from g = 2.40
to 2.41, (b) is from g = 2.95 to 3.02, (c) is from g = 3.09 to 3.28, and (d) is from
g = 3.76 to 3.965
Fig. 44. A double logarithmic plot of the histogram of the localized photon numbers
for the regions depicted in Fig. 43. Figures (a)–(d) correspond to chaotic regions
(a) to (d) in Fig. 43
10 Conclusions
In order to investigate the dynamics of the system, we presented a model
of N two-level quantum-dot systems interacting with optical near fields that
are explicitly expressed in terms of the localized-photonic degrees of free-
dom, which characterize the unique property of the localization of optical
near fields. At the low-density limit, excitons (N two-level systems) are ap-
proximated as bosons, and a rigorous solution of the Heisenberg equation
is obtained. Using this solution, we examined the dynamics of the excitonic
214 A. Shojiguchi et al.
system, and showed that the dipole moments are transferred linearly within
the system. Since the dipole moments in the system represent the quantum
coherence between any two energy levels, this phenomenon might be applica-
ble to a photonic device on a nanometer scale or to the transfer of quantum
information.
For fermionic excitons, the Heisenberg equation becomes nonlinear, and
the dynamics is more complicated. We obtained a perturbative solution given
by (43) within the second order with respect to the localized photon–exciton
interaction, to investigate the dipole dynamics. The study revealed that sev-
eral oscillating quasisteady states exist, depending on the material parame-
ters. Using the effective Hamiltonian obtained from renormalization of the
localized photonic degrees of freedom, we classified such quasisteady states
into several groups: one is a “ferromagnetic” state with all the dipoles aligned
in the same direction, the other is an “antiferromagnetic” state with alter-
nating dipoles, as shown in Table 1.
Note that if the sign of the population difference Wn (0) at an arbitrary
site n and time t = 0 is reversed, then the dipole of the same site at arbitrary
time t, Pn (t), also changes sign (see (44)). Using this flip hypothesis, we can
transform an arbitrary dipole distribution of the system into a dipole-ordered
state after manipulating the initial distribution of the population differences.
This hypothesis is based on the perturbative solution (44), which determines
the sign of Pn (t) according to the sign of the product of P1 (0)Wn (0).
It also originates from the fermionic property of excitons, which gives the
Heisenberg equations of motion for bn (∝ Pn ) as
[Hb−b , bn ] = ∆Ωnm Wn bm . (165)
m=n
The right-hand side of (165) is proportional to Pm Wn for fermions, while
it is proportional to Pm for bosons. Therefore, the occurrence of this kind of
nonlinearity for fermionic excitons is a possible origin of the flip hypothesis.
Moreover, it is not trivial; this property can be applicable to results obtained
from the semiclassical Heisenberg equations (46) and (47).
We examined the radiation property of our system, solving the Heisen-
berg equations (46) and (47) using the semiclassical approximation. In the
investigation, we treated the system as isolated, since the coupling between
radiation fields and the sytem is very weak. The numerical analysis showed
that dipole-ordered states, which have large total dipole moments, show a
large radiation probability comparable to Dicke’s superradiance. In partic-
ular, it predicted that multiple pulses are emitted from the system super-
radiantly. In order to verify whether such a phenomenon is inherent in an
isolated system, we solved the master equation (160) for a dissipative system
using the semiclassical approximation, and found that such multiple pulses
in the radiation profile can even survive in a dissipative system coupled to a
radiation reservoir.
A Phenomenological Description of Optical Near Fields 215
of chaos. Indeed, it was recently reported that the Dicke model is a quan-
tum chaotic system [93]. The term “quantum chaotic system” means that the
distribution of the energy-level spacing of the Hamiltonian differs from that
of a regular system [88]. Since there is no instability of orbits in a quantum
system, due to the uncertainty principle, we need new signatures of chaos for
a quantum system whose classical correspondence is also chaotic. One of the
quantum signatures of chaos is the level statistics [89], and the analysis of
these level statistics remains an intriguing open problem. Since our system
has a richer structure than the Dicke model, we expect our model to produce
more fruitful properties as a quantum chaotic system. Finally, we point out
that the open system discussed using our localized photon model is one of
the simplest nano-optical structures to use to examine the optical response
of nanostructures and to discuss the susceptibilities of nanometric material
systems.
Acknowledgements
References
1. E.A. Synge: Philos. Mag. 6, 356 (1928)
2. M. Ohtsu, H. Hori: Near-Field Nano-Optics (Kluwer Academic/Plenum pub-
lishers, New York 1999)
3. M. Ohtsu (ed.): Near-Field Nano/Atom Optics and Technology (Springer-
Verlag, Tokyo 1998)
4. D.W. Pohl, D. Courjon: Near Field Optics (Kluwer Academic Publishers, Dor-
drecht 1993)
5. J.P. Fillard: Near Field Optics and Nanoscopy (World Scientific, Singapore
1996)
6. M.A. Paesler, P.J. Moyer: Near-Field Optics: Theory, Instrumentation, and
Application (Wiley, New York 1996)
7. M. Ohtsu (ed.): Progress in Nano Electro-Optics I (Springer-Verlag, Berlin
2002)
8. M. Ohtsu (ed.): Progress in Nano Electro-Optics II (Springer-Verlag, Berlin
2003)
9. M. Ohtsu, K. Kobayashi: Optical Near Fields (Springer-Verlag, Heidelberg
2004)
10. K. Shimoda: Introduction to Laser Physics (Springer-Verlag, Berlin and Hei-
delberg 1986)
11. C. Carniglia, L. Mandel: Phys. Rev. D 3, 280 (1971)
218 A. Shojiguchi et al.