Carroll - Distributed Feedback Semiconductor Lasers
Carroll - Distributed Feedback Semiconductor Lasers
Distributed feedback
semiconductor lasers
John Carroll,
James Whiteaway
& Dick Plumb
John Carroll,
James Whiteaway
& Dick Plumb
Copublished by:
SPIE—The International Society for Optical Engineering
PO Box 10, Bellingham, Washington 98227-0010
Phone: 360/676-3290
Fax: 360/647-1445
E-mail: [email protected]
www: http://www.spie.org
While the authors and the publishers believe that the information and
guidance given in this work is correct, all parties must rely upon their own
skill and judgment when making use of it. Neither the author nor the
publishers assume any liability to anyone for any loss or damage caused
by any error or omission in the work, whether such error or omission is
the result of negligence or any other cause. Any and all such liability is
disclaimed.
The moral right of the authors to be identified as authors of this work has
been asserted by him/her in accordance with the Copyright, Designs and
Patents Act 1988.
Preface xiii
Acknowledgments xv
Principal abbreviations xvii
Principal notation xix
3.5 Scaling 90
3.6 Horizontal guiding: effective-index method 90
3.7 Orthogonality of fields 92
3.8 Far fields 92
3.9 Waveguiding with quantum-well materials 93
3.10 Summary and conclusions 95
3.11 References 96
Index 405
Preface
The authors have had the pleasure of working for many years within
one of the exciting areas of research and development: sources for
optical communications, an area which has seen phenomenal growth
over three decades. The semiconductor laser diode has proved to be
an essential device and one form or another of a distributed feedback
semiconductor laser (DFB laser) is a key component for modern
optical communication systems. Technically, the term DFB implies a
particular form of laser but to the authors the term has come to mean
any laser in which there is distributed feedback or where Bragg
gratings are integrated within the device.
This is a research monograph which combines a high level of
tutorial material with a research review. It will be useful for
postgraduate courses in optoelectronics studying either systems or
devices, research workers in optical systems, lasers and optoelectronic
devices, and also be of use forfinal-year-projectstudents in BEng and
MEng degree courses. Lecturers and research workers in the field of
laser diodes will find this a useful source of ideas, references and
tutorials which reach the explanations other books do not reach.
Combined with this tutorial material are accounts of the practical
operation and design of DFB lasers as used by the communications
industry along with discussions of the operational reasons for
requiring key modifications within the laser in order to meet specific
performance targets.
The tutorial material is distributed throughout the text and the
appendices. The different orders of operation of a Bragg grating are
outlined, novel concepts of the stopband diagram are explained, and
the modelling of laser interaction is based on sound fundamentals of
energy balance so that both particle models of lasers and electro-
magnetic-field models can be tied together. There are many different
ways of estimating the spontaneous emission, and again the con-
xiv Preface
The names of all our friends and colleagues who have participated in
Laser Workshops, COST240, and key laser conferences over relevant
years are too numerous to be mentioned individually but they all have
had an influence on our thoughts about lasers and all are thanked. We
ask forgiveness of those who feel they should have been listed.
Thanks are also given to the IEE and in particular Jonathan Simpson
for support in publishing this book.
The support from our families, while apparently taken for granted,
has never been forgotten by us.
Principal abbreviations
AR antireflection
CCCH Auger process involving three conduction-band
electrons and one heavy hole
CNR carrier-to-noise ratio
DBR distributed Bragg reflector
DFB distributed feedback
EDFA erbium-doped fibre amplifier
FP Fabry-Perot
GaAs gallium arsenide
GaAsxP,_x gallium arsenide phosphide
Ga.Al^As gallium aluminium arsenide
C^rCt Tn Ac P gallium indium arsenide phosphide
GAVCF grating-assisted vertical-coupler filter
HOFC higher-order Fourier coefficients
InP indium phosphide
IVBA intervalence-band absorption
L-I light-current
MATLAB registered trade mark of a MathWorks computing
environment
MBE molecular-beam epitaxy
MOCVD metal-organic chemical-vapour deposition
MQW multiquantum well
SCH separate-confinement heterostructure
RIN relative intensity noise
SiO2 silicon dioxide
Si3N4 silicon nitride
SMSR side-mode-suppression ratio
SNR signal-to-noise ratio
SOA semiconductor optical amplifier
TDM time-division multiplexing
xviii Principal abbreviations
Sf signal bandwidth
A/ system bandwidth, gain bandwidth
A/w/ linewidth
Av/>/ spontaneous bandwidth for laser
A^/ change in frequency with change in gain
8(zl — Zo) Dirac delta function of position
Ap elemental reflection per Bragg period
AT elemental transmission coefficient per Bragg
period
8t computational step length in time
Ax uncertainty in position
A/? uncertainty in m o m e n t u m
8z space step related to St= 8z/vg
e, exp exponential
E optical electric field also as a function of zco
time, etc.: E(z), E(<o), E(t), E(y, z)
E o, •^racb A/;/(r) subscripted versions of electric field
,Er forward, reverse electric field
energy of particle
energy levels (e.g.: c, v —^conduction, valence
band)
\ > ©act,% :lad energy gaps appropriate to different mate-
rials
P quasi-Fermi levels for electrons (holes)
permittivity of free space (Chapter 4)
eA relative permittivity, same in layers 1 and 2
effective relative permittivity
&rreffr S
ri relative a n d effective, real a n d imaginary
permittivities
e gain-saturation parameter where G(N) =
G(N-Ntt)/(l + eS)
h Planck's constant
rjexl external quantum efficiency
k Boltzmann's constant
xxii Principal notation
L laser length
A, Ao free-space wavelength
A^ free-space Bragg wavelength
kelec(rm electron wavelength
A,m, Xx, A2 wavelength in material/guide, material 1, 2
etc.
\pmk wavelength of peak gain
A, A' Bragg-period first-order grating, period of
second-order grating A' = 2A
q electronic charge
Q cavity Q ratio of stored energy to (2TT • energy
lost per cycle)
t time
t(. thickness of the contact layer
T transmission coefficient at facets
T, Tn absolute temperature, noise temperature
T integer giving time steps
T j , T2 transistor labels
r time constant of spectral filter
Tp9 Tfp effective p h o t o n lifetime i n c l u d i n g t h e dis-
tributed losses
rsp lifetime for s p o n t a n e o u s radiative emission
rr spontaneous recombination time constant
including all mechanisms
rro recombination time constant when A/= Ntr
hnimband intraband relaxation time ~10~ 12 s
Principal notation xxv
y co-ordinate distance
z axial distance
Z integer giving space steps
Chapter 1
The semiconductor-diode laser
1.1 Background
The invention of semiconductor diode lasers [1-4] along with their
ability to be directly modulated with information [5,6] and the
invention of optical-fibre communication [7-9] in the 1960s-1970s
have together had far-reaching effects on telecommunications leading
to extensive optical-fibre communication networks and systems
[10-12]. This book explores one aspect of this success story. High-
performance optical-fibre systems need sources with a number of
features:
(i) stable single-frequency output;
(ii) modulation capability of gigabit/s;
(iii) stable operating lifetimes measured in 106—•lO8 hours at room
temperature; and
(iv) manufacturability.
The distributed-feedback (DFB) laser has these qualities and is now a
standard device. The design and development of these lasers has
proved to be fascinating work for those involved, and the physical and
numerical modelling which has gone into developing these lasers is
the topic of this book. Although the text concentrates on lasers with
distributed feedback, the modelling work presented is more widely
applicable and it is hoped that the reader will be as excited as the
authors by the unfolding story.
Eo=EoPleflPnghle^-^Le~^L (1.3)
This results in the 'amplitude' condition for stable oscillation (lasing)
given by
/^A^e^'-^M (1.4)
or
4 Distributed feedback semiconductor lasers
divergent
output beam
reflection
optical length 0 Pright
(N half wavelengths)
reflection n
r .,
left
Figure 1.1 Operation of a Fabry-Perot laser
Showing schematic of standing electric-field waves inside the resonator
1
(1.5)
Pleft Pright
0.5
A. 7V A A
-2000 0 2000
Frequency, GHz
Figure 1.2 Power spectrum for a Fabry-Perot laser
250 ^m Fabry-Perot laser operating at a central wavelength of around
1.55 jmm
magnitude of the round-trip gain Ground is ever so slightly less than unity
and it is the random spontaneous emission (discussed in Chapter 2)
which is amplified through gain and feedback, by the factor of
1/(1 - Gmund)~l03 to 106. Because the gain and spontaneous emission
spectra cover several nanometres change of wavelength, the phase
condition /3L=NTT (N integer) is satisfied at a number of frequencies,
or modes, within that band but the gain condition ensures that the
spontaneous emission is only amplified significantly at modes where
11/(1 - Ground) I is very large. As IGrowwrfl becomes smaller away from
Apeak, the amplification of the modes by the random excitation grows
weaker. The spectrum of the FP laser in general shows multiple modes
with frequencies determined from the phase condition giving a mode
spacing AFPf=(vg/2L) where vg is the group velocity for the optical
fields in the laser's waveguide. Figure 1.2 gives a typical spectrum for
a 250 juim Fabry-Perot laser operating at a central wavelength of
around 1.55 |xm showing these features.
With an ideal design, the laser's output power is proportional to the
drive current above threshold, as in Figure 1.3. This is useful for direct
amplitude modulation of the optical output by modulating the current
drive. This ideal linear behaviour is modified by heating which reduces
the optical gain and so increases the threshold. Heating and gain
variations can also alter the spatial variations in optical and electronic
6 Distributed feedback semiconductor lasers
T ' ^
3
optical
output / / >''
power, 2
mW
/ / ''''
1
/
/ /
0
20 40 60
input current, mA
Figure 1.3 Typical optical output-power/current curve for FP laser
and curves show the effect of increasing temperature T (see
Section 2.6, eqn. 2.44)
densities within the laser, thereby also changing the ideal linearity.
Modulation at gigabit rates creates additional phenomena such as a
change of frequency (called 'chirp'; discussed in Chapter 4).
It follows that lasing action can occur only in semiconductors which
exhibit optical gain; this normally implies a direct-band-gap semi-
conductor (see Chapter 2), and as a consequence silicon, so widely
used for electronic devices, is unsuited to make lasers. Furthermore,
the need for low threshold currents and stable optical modes requires
the complex technology of different materials for heterojunction
diodes as discussed briefly next.
0 > 0 critical for total internal reflection cross sectional side view
optical
intensity
ray
Figure 1.5 Rays and intensity for confinement of optical fields to an active layer
a Ray picture (total internal reflection)
b Field or modal picture
tion also normally decreases its refractive index. In Figure 1.5, the
outer cladding or confining layers have a slightly lower refractive index
than the active or gain layer. The layers act as a slab waveguide
confining the photons closely around the active layer through the
action of total internal reflection. The simple optical-ray picture is
shown in Figure l.ba. Total internal reflection (Appendix 1) requires
that the refractive indexes of each of the outer materials, n^ and ?%, are
smaller than the refractive index nx in the active region. A guided ray
has to propagate at a sufficiently shallow angle to the active region, and
from Snell's law, the angle of incidence 6>dcrilical. One may therefore
write:
sin 0>sin 9ailiral = n^/nx = \ - An/nY here An=nx - n% and n^ru, (1.7)
The angle <f> shown in Figure 1.5 then has to be 'glancing' such that
<f>= ( T T / 2 - 0 ) <COS~1 ( l - A n / n ^ V ^ A n / n ) (1.8)
where n is the mean index. GaAs has a refractive index of 3.6,
corresponding to nY, and Ga^^Al^As is lower by An— — 0.66x; typically
x=0A, so An is about 0.26.
The conditions in eqns. 1.7 and 1.8 ensure that optical rays which
propagate down the active region are reflected from side to side. This
ray picture, however, fails to show the evanescent fields in the cladding
layers which have to accompany total internal reflection. The more
accurate wave-guiding picture, shown schematically in Figure 1.56,
indicates that, although optical power does not propagate outwards
from the active layer, an evanescent wave penetrates to the order of a
wavelength into the cladding layers. The details of waveguiding in
multislab guides are discussed in Chapter 3 and Appendices 1 and 2.
The semiconductor-diode laser 9
-1501
emitting spot s t r j p e contact
• - 0 .5 |um
Vertical'
V ^ heterojunctions
p contact ^ P GaAIAs
semi-insulating GaAIAs
\ n GaAIAs
n contact
contact to ridge
& p contact metal
• - - • emitting spot
-- n confining layer
(& substrate)
Figure 1.8 End viezv of 'ridge-waveguide' laser
Ridge typically may be 5-10 JJLITI wide with active layer beneath, -0.2 juim
thick
r
with effective active
width w
confining layer a vert
,' to ^
1 monitor.
diode laser
(a) (b)
AR coating
output
grating lens
diode laser
Figure 1.12 Laser with external grating reflector
design allows linewidths of less than 100 kHz and tunability over tens
of nanometres to be achieved, but the mechanical assembly is too
expensive and delicate for such lasers to be used in the robust
environment of commercial communication systems, especially as the
maximum frequency for direct modulation is limited.
optical amplifier
or Q B Q O D O Q Q D O Q QQQ Of
(erbium doped fibre)
laser + modulator
M-HH photo- spectral
nun '\C
input signal dispersed and detector filter
attenuated
modulated optical amplifier
carrier demodulator
TRANSMITTER signal conditioner
HfHi mm
restored/retimed signal
UjOaQflBQQQQlLQffllD
standard fibre and laser sources becoming centred around the 1.3 |xm
and 1.55 |xm wavelengths.
Although the attenuation in optical fibres is sufficiently low to allow
links of well over 100 km using standard 1 mW FP lasers at 140 Mbit/s,
at higher bit rates there is a significant dynamic wavelength shift (or
'chirp') in the laser's emission which causes loss of information,
regardless of the laser's output power, because of the dispersion in the
fibre. Dispersion is the phenomenon whereby, because of wavelength-
dependent phenomena within a material, different wavelengths in the
fibre have different speeds of propagation. Consequently a pulse
which contains a range of wavelengths will spread as it travels.
Eventually, given a long enough fibre link, one pulse spreads into an
adjacent pulse and causes interference. At 1.55 |xm, where the
attenuation in standard silica fibre is lowest, dispersion is significant
amounting to about I7ps/kmnm. For such standard silica fibre,
dispersion is nominally zero at 1.3 |mm, but laser wavelengths are very
rarely matched exactly to the dispersion minimum of the particular
fibre in a link (due to manufacturing and installation tolerances, as
well as temperature variations) so low chirp is desirable from lasers,
The semiconductor-diode laser 17
angle of
incidence angle of
diffraction
grating period
blazed A
reflection
grating
The Bragg grating buried inside a laser may not be the final answer
to the problems of precise frequency determination and high CNR,
but it provides an excellent starting point and a fascinating study. This
chapter therefore concludes with an introduction to the Bragg grating
and how it may be buried inside a practical laser.
angle of
wave-front diffraction
4TIM
^ incidence^
wave-front
Figure 1.15 Diffraction from grating embedded in laser waveguide
For a grating embedded inside a laser, 0rrr/2 and for feedback the
diffracted wave 0d~7r/2. However, one can also have 0d~O when the
diffracted wave gives radiation out from the guide. The diffracted
wavefront is reinforced because there are 2T7"M radians per pitch where
M is an integer and is known as the order of diffraction
A sin ft A sin 6d
(1.10)
m2
where n} and n^ are the refractive indices in the central and outer
regions of the waveguide, respectively. Inserting eqn. 1.11 into eqn.
1.10 gives
20 Distributed feedback semiconductor lasers
(1.12)
I
05
C
"I
cO .Q
05
1
i radiation los
II 0
0 CO
LL CC CE
CO CO CO
>
forv
UO1
forv
1 O T3 I D
0 0
ID
0
1
II 0 0 0
y of feed
i-1
c
di "5 g1
1 0
Sum
radiation loss only in the first and second diffraction orders at angles
symmetrically disposed about the normal, but gives feedback in the
third diffraction order. In general, feedback is obtained from a grating
of order N in diffraction order N, but diffraction orders intermediate
between 0 and N result in radiation loss for such a grating.
In principle, just as the blaze or shape of the grating in a
spectrometer alters the power division between the different orders of
diffraction, so also one can change the reflection and radiation
[42-44] by tailoring the grating-tooth shape. One may also need to
distinguish between the TE and TM polarisation within the guides
[45]. In general, the requirement to mass manufacture gratings buried
within a laser gives severe constraints on the tooth shape which is
usually controlled by the chemical etching and crystal orientation.
Details of wave propagation and gain within such gratings with a
distributed reflectivity are covered at the end of Chapter 4 and
beginning of Chapter 5. A method of calculating the reflectivity per
unit length is discussed in Chapter 4 with the help of a MATLAB
program slabexec.
1.6.2 Fabrication of gratings inside lasers
The theoretical concepts of lasers with distributed feedback [46] were
developed well before the technology could deliver a reliable and
useful method of fabrication. The early semiconductor lasers with
distributed feedback operated at low temperatures and had short
operating lifetimes [47,48] but in the mid 1970s GaAs lasers with
distributed feedback could be made to operate at room temperature
[49], and then developments with new materials led to room-
temperature operation of lasers, with distributed feedback, around
1.3 fxrn and 1.55 |xm [50,51] in the early 1980s. The details of
fabrication of semiconductor lasers and the techniques of substrate
formation and subsequent epitaxial-layer growth using metal organic
chemical vapour-deposition (MOCVD) or molecular-beam epitaxy
(MBE) are outside the scope of this book and are far from trivial. The
research and development of materials-growth techniques for forming
the confinement-layer/active-layer heterojunctions of lasers have been
vital. All that can be given here is a mere glimpse of this marvellous
interaction of materials and technology that has led to reliable
commercial devices. Figure 1.16 illustrates some principal processing
steps needed to form a semiconductor laser with a uniform embedded
grating which creates the feedback as described in Section 1.6.1.
The first step is to grow, by the methods of epitaxial growth, good
laser structures with the substrate followed by the confining layers and
The semiconductor-diode laser 23
460 nm waveguide
period
grating
active
layer
stripe contact
55° p-lnP
effective
permittivity
for 0.46 |im * e min
mode
Figure 1.17 Schematic cross-section of DFB laser grating etched in epitaxial-layer
structure
passive
waveguide Bragg grating optical gain Bragg grating
layer reflector section reflector
;XX>: A, 12 \ N
ooo 12
m
n
l
X 12 X 12 X 12
m m m
reflector section centre (gain) section reflector section
Figure 1.19 Grouping of a grating into \m/2 sections where all reflections Ap from
each section add
gratinc
CO OC 00 CO
distributed feedback
forward
reverse
scattered into the optical wave in the reverse direction. Both waves
grow as they travel towards their respective facets because of the
feedback and the gain. The abrupt changes in refractive-index, with
each section of length Am/4, are purely schematic to emphasise the
main features. This is the distributed-feedback (DFB) referred to
earlier. However, it is found that the DFB laser oscillates at two possible
frequencies slightly removed from the Bragg frequency depending on
KL (Chapters 4 and 5). The reason for this is that the phase condition
for oscillation cannot be met at the Bragg frequency, as discussed
below.
Now reconsider the DBR laser. It is found that the presence of gain
enhances the reflectivity of the Bragg gratings and this increased
reflectivity dramatically reduces the required gain for the FP-gain
section. A particular utility of the argument used in Section 1.7.1 based
on the Am/2 sections above comes in the ideal design of the central
section for a DBR laser required to operate at the Bragg frequency of
the grating where there is the maximum reflection. It is possible to
ensure that the round-trip phase matching occurs, as in a Fabry-Perot
laser, if the two reflectors are identical mirror versions of each other
with a real reflection of the same sign, and then these reflectors are
separated by a central section of length NxAm/2 (Figure 1.19). The
shortest possible Fabry-Perot section is Aw/2 long (N = 1) and this then
requires the insertion of an additional Am/4, say into the lower-
refractive-index region to form a length Am/2 between two mirrored
sections of Bragg gratings (Figure 1.21a).
28 Distributed feedback semiconductor lasers
(a) DFB laser with single ^ m /4 phase shift (b) DFB laser with two J^/8 phase shifts
i'' h
IJUUMJ.JUUU
h
gam sections gain section pnase section gain section
X
mmmmmm
a can be used for low chirp
b and a for tunable lasers
c integrated modulator
gain sections
driven harder than the other, and under modulation the drives to the
ends are reversed, this being known as 'push-pull' operation [55] (see
also Chapter 8). The total gain and electron density in the laser remain
constant, but optical power is switched from end to end (Section 8.3).
Linewidths of less than 0.01 nm spread in wavelength are predicted to
be possible at 10 Gbit/s with these lasers. Structures such as that shown
in Figure 1.226, with a plain gain region in the middle, can be used to
give a wide range of tuning for the laser [56,57].
'DFB' lasers with more complex grating designs have been designed
to give wide-range tunability though not allowing high-speed modula-
tion [58,59]. Here the lasers are designed with an additional
periodicity of several wavelengths superimposed on the A/2 Bragg
wavelength. This then allows a different periodic reflection-spectrum
from each end of the laser, each of which can be shifted slightly by the
small refractive-index reduction caused by carrier injection. 'Vernier'
action, discussed in Section 8.4.3, then occurs so that the laser operates
only at a wavelength where reflection peaks from both of the gratings
coincide. This has allowed quasicontinuous tuning over greater than
100 nm centred on 1.5 fxm. Vernier action is also possible where the
waveguide for the laser is branched into a Y with the favoured mode
having a common peak amplification/reflection from each arm of the
Y (see Section 8.4.4) [60].
For applications which demand narrow stable linewidths, structures
such as that shown in Figure 1.22c are sometimes used [61,62]. The
main part of the laser is run continuously, and hence has a narrow
spectrum, and the modulator section simply switches the output on
and off. In practice, the modulator interacts with the master laser, and
some excess chirp results, but designs are improving, and many more
multisection, multicontact lasers with increased functionality may be
expected in the future.
30 Distributed feedback semiconductor lasers
Bragg contact
reflectors f
top diameter - 10 \i m
height ~ 20 |i m
1.8 Summary
This chapter has set the scene for the operation and fabrication of
semiconductor lasers, emphasising why the authors see lasers with
grating reflectors as so important. The earliest lasers were examined,
showing that a key requirement is to confine the photons and
electrons to the same physical region. A geometrical factor T, known as
the confinement factor, expresses the degree to which this has been
The semiconductor-diode laser 31
achieved. The fact that light is harder to confine closely around a given
region than it is to confine the electrons around the same region
means that F takes values typically ranging from 0.3 to 0.01, depending
in detail on the technology that is used. The simplest requirements for
lasing were put forward so that the reader understands the impor-
tance, just as in electronic feedback oscillators, of the round-trip gain
and phase in achieving a stable oscillation. Many factors influence the
net efficiency of a laser other than the facet-reflectivities (Section
1.4.1). Nonradiative recombination, waveguide absorption and scatter-
ing losses reduce the quantum efficiency for converting electron/hole
pairs to photons from the ideal of 100%. Further reductions occur
because not all charge carriers, injected at the contacts, reach the
active region. Carriers spread away from contacts, leak across blocking
layers or over heterobarriers and suffer interfacial recombination.
Typical single-facet effieciencies -10-30% dependant on structure and
processing.
The importance of single-mode lasers for modern optical commu-
nication along silica fibres over many tens of kilometres was noted. A
variety of different types of laser were outlined, and in particular it was
observed that to gain a stable single mode the favoured technology is
to incorporate some form of Bragg grating. The fundamentals of
gratings were introduced and the elements of different laser designs
with gratings buried in their structures were presented.
1.9 Bibliography
1.10 References
1 HALL, R.N., FENNER, G.E., KINGSLEY, J.D., SOLTYS, TJ. and CARL-
SON, R.O.: 'Coherent light emission form GaAs junctions', Phys. Rev. Lett,
1962, 9, pp.366-368
2 NATHAN, M.I., DUMKE, W.P., BURNS, G., DILL, EH. and LASHER, G.:
'Stimulated emission of radiation from GaAs p-n junctions', Appl. Phys.
Lett, 1962, 1, pp. 62-64
3 QUIST, T.M., REDIKER, R.H., KEYES, R.J., KRAG, W.E., LAX, B.,
MCWHORTER, A.L. and ZEIGLER, H.J.: 'Semiconductor laser of GaAs',
Appl. Phys. Lett., 1962, 1, pp. 91-92
The semiconductor-diode laser 33
2.1 Introduction
The operation of semiconductor lasers is strongly affected by the
materials from which they are made. Spontaneous emission which
initiates lasing, optical gain which is essential to achieve lasing, and
other processes involved in lasing all use quantum processes at the
level of single atoms and electrons within the lasing material. Indeed,
without simplifications, the physics and mathematics necessary to
describe such atomic systems fully is too complicated, certainly for the
level of this book. The simplification and approximations must,
however, be done in a way which is adapted to the requirements of
semiconductor lasers, and the limitations must be understood. The
results and implications of quantum physics are discussed here but
there are only illustrative outlines of any derivations. Greater depth
may be obtained from specialist material [1,2] (see Bibliography).
In the context of lasers, the atomic behaviour of semiconductors can
be summarised under three headings:
(i) the distribution of energy states;
(ii) the occupation probability of those states; and
(iii) transitions between states.
The first few sections of this chapter will cover these headings at the
atomic level, while the later sections will apply them at a macroscopic
level to optical gain, spontaneous emission and noise.
momentum (p)
a
Figure 2.1 Energy-momentum diagrams
(i) In free space
(ii) In a periodic lattice of period a with p=hk/27T and k=27r/\dmron
may move around in the crystal. Simple Newtonian mechanics forms a
starting point where the kinetic energy of a particle % - %c is measured
from the conduction band energy edge %n and given in terms of its
momentum p and effective mass mc from
(2.1)
(2.2)
heavy electron
conduction normal electron
bands
valence _
bands _»#T jL ^ • • L ^ ^T~ h e a v Y h o l e
light hole
hole split off
Figure 2.2 %-k (energy-momentum) diagram for a direct-band-gap semiconductor
Lowest energy in conduction band is directly above highest energy in
valence band. Hole-split-off band is caused by coupling of electron spins
and orbital angular momentum
conduction
band
energy
momentum
Figure 2.3 Transitions in indirect band-gap material
For an electron to move from the conduction band to the valence band
requires momentum and energy to be conserved. Only a photon has
enough energy to balance the energy loss and only the phonon has
enough momentum to balance the momentum change
conduction
conduction
band band
energy
photon momentum
momentum
h/X <« h/a
m emission absorption
Figure 2.5 Transitions between energy levels in the conduction and valence bands
As A^ changes so \k\ varies. Over a band of energy d%, there is a
corresponding band or range d& in I k\
(i) Since Auger processes involve three particles, they increase with
the cube of carrier density.
(ii) All Auger rates increase with decreasing bandgap because, in
general, the range of likely transitions increases with decreasing
bandgap.
(iii) Processes such as CCCH, where an increase in temperature
increases the numbers of electrons available to participate (as in
practical lasers), will have rates which increase with temperature.
(iv) Since numbers of dopant atoms or defects are almost constant
with temperature, Auger processes involving these increase less
rapidly with temperature than the processes of (iii) above.
As a direct consequence of (ii) above, it can be seen that 1.55 jmm,
1.3 juim and 0.9 jxm lasers are steadily less affected by Auger recombi-
nation as the wavelengths decrease, as indicated by the sensitivity of
the threshold currents in lasers to temperature with Auger processes
typically becoming significant at around 50°C for 1.5 |mm devices, 70°C
at 1.3 |xm, and 120°C + for devices operating at 0.9 |j,m or less.
N 1 states - Q . ^ ^ . ^
dS /dt
SPont
n g electrons p holes
isotropically and also over a much wider band of frequencies than the
lasing linewidth. Consequently only a very small proportion f3sp
(typically 10~6 to 10~2) couple to the lasing mode and eqn. 2.9 must be
modified to read
dS/dt=KstimS(NN^PNl - NXN2) +(3spKspNP (2.10)
Under lasing conditions when injected charge dominates over doping,
N^P, and P may be eliminated, allowing a photon-interaction
equation often used in modelling:
*= G'n(N- Nh) S+p^BN2 (2.11)
One may link the parameters in eqn. 2.11 with those in eqn. 2.10 as
follows:
p (2.12)
where Gm is the differential power gain/unit time, and Nlr is the
transparency value of electron density at which the stimulated optical
gain is cancelled by the stimulated optical absorption—the material is
transparent for the light at that frequency.
A nonzero acceptor density Na in the gain region causing P=N+Na
from electrical neutrality could, from the theory, reduce the transpar-
ency electron density. However, a high density of acceptors increases the
nonradiative recombination and in practice this technique is not
normally used deliberately, though ^dopants such as Zn may diffuse
from the />-cladding layer during growth. Here doping in the gain
region is ignored.
From eqns. 2.7 and 2.8, if gain is to exceed absorption, then one
requires
NP~ (Nj - N) (N2 - P) =NP{1 - (1 - Nx/N) (1 - N2/P)}>0 (2.13)
But iVand Pare related to Nx and N2 via the Fermi relationships of
eqns. 2.3 and 2.4:
(2.14)
N2/P={l+exp(%jp-%2)/kT} = 1/FP(%2) (2.15)
and so, for net gain (substituting the above into eqn. 2.13)
k i
t { ( « « ) ( ^ V l / ^ ( ^ ^ ^
Eqn. 2.16 confirms the earlier result that the voltage Vj applied to the
laser junction, given from the separation in quasi-Fermi levels, must be
greater than the equivalent bandgap voltage to obtain lasing action [9].
A further implication is that, if the contact resistances in the laser
Gain, loss and spontaneous emission 49
structure (see Section 5.7) could be made very low, the carrier density
in the active region would be determined solely by the applied voltage
and the optical gain could then be regarded as a direct function of
applied voltage, rather than the electron density.
Now in thermal equilibrium, the net photon-emission rate dS/dt is
zero with no voltage applied across the junction (i.e. %JN-%J^). Also,
select the situation where spontaneous and stimulated emission cover
the same frequency ranges and angles so that the coupling factor
/3sp=l. These conditions force the system to be in equilibrium at a
temperature T with black-body thermal radiation described by
Planck's formula [19] for the photon density over a range of
frequencies, say kspf, around an optical frequency/:
S=87rfn\ Aspf/c3{cxp(hf/kT) ~ 1} (2.17)
Here n is the refractive index and ngis the group index, i.e. c/ng=vgis
the group velocity or the velocity at which the energy propagates while
c/n=vp is the phase velocity. In equilibrium, dS/dt=0, and eqn. 2.10
becomes
0=KsltmSNP- K^SNP^/N- 1) (N 2 /P- 1) + KspNP (2.18)
Also in equilibrium %jN=%jTso that eqns. 2.14 and 2.15 yield
S=(Ksp/Kslim)/{exp(%i -%2)/kT- 1) (2.19)
In the simple two-level system, hf=%x —%2 so comparing eqns. 2.17 and
2.19:
(Ksp/Kslm) =8irfn\ bspf/c5 (2.20)
A relationship having been established between spontaneous and
stimulated emission coefficients, these may now be related to the
differential of the (gain/unit time) with electron density given from
Gn. From eqn. 2.12, Gm=Kstim{H^H2)) and writing (N^N*,) =NreJ then
using this as a reference density where N=P=Nref in eqn. 2.6 so that
Ksp= l/Nr(,jTsp, o n e obtains
G ^ f f f n g r ^ J ) (2.21)
Using the group velocity vg to link the spatial distance z travelled by the
light with the time Ogives z=tvg=ct/ngand the material's optical-power-
gain/absorption per unit distance may be written as:
(1/5) dS/dz=c2(iY-A^r)/(87r/2w2T,/;A,/;j0 (2.22)
showing the important link between spontaneous recombination and
stimulated gain and absorption. This is related to the quantum theory
50 Distributed feedback semiconductor lasers
The analysis above refers to just two levels of energy when in reality
there is a band of energies in both the conduction and valence band
as in Figure 2.5 where one is dealing with pairs of electrons and holes
just inside the conduction band and valence band, respectively. Using
the notation of Figure 2.5, and assuming parabolic %/k relationships
with ^-selection of the electron/hole interactions with negligible
photon momentum:
(2.24)
where \/mr=\/mc+\/mv gives a combined effective mass ('reduced
mass') from the effective masses in the conduction and valence bands.
Now for a spread of photon energies d(h/) there is a spread dk in
the magnitude of the electron momenta and a consequential
spread of their energies in the conduction band given from
d% = (f\2k/mf) dk=(mr/mr) d(h/). There then has to be a similar spread
of energies d(S&v) of holes because they have the same k in the valence
band, and so d(Mt) = (fi2k/m,) dk={mr/mv)d%.
Now consider a density of photons $ d(h/) spread over a range of
frequencies d/ around a frequency / interacting with a density of
available electrons in the energy states equivalent to Nj at an energy %x
in Section 2.3.1 but now %1=%=%(.+ 8%r and takes a spread d%. This
density of available electrons is found by first finding the density of
electron states (given by eqn. 2.2) with the range of energies d%
associated with the spectral spread d/:
emission as
(
-5) (%-~%c)l/2FN(%) d% (2.26)
h-/
Using a similar argument, the density of vacant sites or holes available
to receive the electrons for stimulated recombination is the equivalent
of N2FP(%2) in Section 2.3.1 where the energy %2=(%-hf) will be
below the valence-band edge. Bearing in mind that the spread of
energies (mc/mt) d% in the valence band corresponds to the spread d^
in the conduction band for the same ^-values, then
N2FP(%2) becomes
3/2
/9 \
47r( - J (%-%+hf)l/2FPCZ-hf)(mc/mt) &% (2.27)
Now h/along with the spread d(h/) will be held constant but the whole
range of % can be permitted, consistent with %^%(. and % —hf^%v,
i.e.
«p (2.28)
Then on summing or integrating over the whole energy range of
permitted %, the equivalent stimulated emission rate given from eqn. 2.7
now becomes
(dS'/d 0 slim misswn d (h/)
=AslmS'{j(% -%)x/2{%v~~%^f)l/2FP{%--hf)FN{%) &%} d(h/) (2.29)
(2.30)
giving the result that
52 Distributed feedback semiconductor lasers
10 _
relative
(l-FN)(l-FP)=FNFPexP{(hf-qVl)/kT} (2.31)
Hence one finds for the absorption
(dS'/dt)stimahsmpiton d(h/) =AstmS' exp{(h/- qVf)/kT}x
(2.32)
with the same range of % as in eqn. 2.28.
Now the gain spectrum is given from the net stimulated (emis-
sion — absorption) so that approximating the Fermi-Dirac distribu-
tions to unity over the limited range of energies, one can determine
the gain spectrum as approximately
(dS'/dz)/S'=2gtn*(ng/c)Astim x
{l-expj(h/-^)/kr}}[J(«-«,) l / 2 {« r --« + (h/-%)} 1/2 d < g] (2.33)
The integral may be evaluated by putting y2=(%-%c) to find a rough
but explicit expression for the gain spectrum with ^-selection as a
function of the junction voltage:
gm(hf, F/) = (7rn/ ( ,y8c)(h/-^) 2 [l-exp{(h/- 9 V,)/kT}] (2.34)
Figure 2.8 sketches this dependence of gain on the optical frequency
and the junction voltage. As the junction voltage increases, so the peak
optical gain increases in frequency. One can also see that as kT
Gain, loss and spontaneous emission 53
gain, cm
200
1.8 x 10 1 8 cm' 3
150 -
100 -
-50
-60
-70
-1
gain, cm Af ~k7)h|,
50 r~ 1.6 x 10 1 8 cm' 3
-50
-100
(a\ spontaneous
dipole solid angle for evert
radiation into free space
2 effective area A
= Q. '/«eff
FWHM emission angles Q ,0
vert ho
frequency
Figure 2.12 Schematic of the principal elements forming the spontaneous emission-
coupling coefficient
100
gain 1.33JV
50
1.33AT
(1/cm)
o
-0.01 0.01 -0.01 0 0.01 -0.01 0 0.01
log(atfQ) log(cnflB)
Figure 2.13 Schematic links between gain and refractive index changes with frequency
caused by Kramers-Kronig relationships
The frequency H is the frequency of the gain peak. The real refractive
index is given by the phase, measured in inverse distance to have the
same dimensions as the gain (imaginary refractive index). The real
index level can be set at an arbitrary level without altering the frequency
dependence; only the relative changes of phase are of concern. As the
electron density increases, the gain (imaginary refractive index)
increases and the real refractive index decreases. The parameter aH can
then be estimated
A7.16 and limiting the spectral range to values around the gain peak.
Even at the peak of the optical-gain mechanism, there is a separate
physical effect caused by a plasma of nearly free electrons which adds
a component into the relative permittivity, reducing the real refractive
index with increasing electron density (Appendix 7). As the frequency
moves away from the peak gain, this change in the (real) refractive
index is modified by the gain curvature through the Kramers-Kronig
relationship. At any particular wavelength, one can define the ratio
known as 'Henry's a factor' or aH:
aH= - {(dnr/dN)/(dnt/dN)} (2.40)
where the negative sign is inserted to make aH positive in this work.
When amplitude modulating a laser, by turning it on and off, there
have to be changes in the gain. From eqn. 2.39, the field gain/unit
distance g=coni/c = 2Trni/Afs, so that from eqn. 2.40 changes Agin the
gain give changes in the refractive index from
An = -aHAfsAg/27r (2.41)
The laser frequency is usually controlled by an optical length such as
the optical length of the Bragg period or the Fabry-Perot cavity, and
is therefore directly proportional to the average refractive index.
Consequently the magnitude of aH directly affects the change of
frequency A^/with a change in gain Agin a laser by approximately
1.6 x 10 18 c m " 3
0CL
where the bar denotes some appropriate average over the laser. These
dynamic wavelength shifts are known as (dynamic) chirp with dynamic
changes in frequency following dynamic changes in the gain. Different
laser structures using the same material will have different chirp
because of structure-related effects, and especially if the operating
frequency can be controlled in relation to the gain peak. Figure 2.13
shows a rough estimate of the changes in aH with operating frequency,
indicating that aH rises markedly at the lower frequencies and reaches
a minimum above the gain peak.
This effect can be understood by splitting the refractive-index
dependence into the component caused by free electrons (also
referred to as the plasma effect) where increasing the electron density
reduces the real refractive index over a broad range of frequencies
(Appendix 7, Figure A7.3). The localisation of the gain gives rise to an
associated ^-shaped change in the real refractive index, centred about
the gain peak. Above the frequency of this gain peak, the net real
refractive-index change with electron-density change reaches a mini-
mum. Figure 2.14 shows measured data indicating the reduction in aH
at high frequencies (short wavelength), as also shown schematically in
Figure 2.13.
64 Distributed feedback semiconductor lasers
In Fabry-Perot lasers, the laser must operate close to the gain peak
so that major variations in aH are not possible. In DFB lasers, the
choice of operating frequency is determined more by the grating, so
that the choice of the wavelength for peak gain gives an extra degree
of freedom and allows one to operate at frequencies where there is a
lower value of aH than for a Fabry-Perot laser, and so in principle one
can reduce the chirp. At 1550 nm it is common to offset the lasing
wavelength by around 20 nm down in wavelength from the gain peak
to take advantage of this lower value of chirp.
Quantum-well material allows one to engineer the gain profile to a
much greater extent than merely changing the composition in bulk
materials [8]. If the gain peak could be made narrower about its centre
frequency, then the ^-shaped change in the real refractive index
becomes stronger and gives a bigger offset for the plasma effect [33]
(see Figure A7.3). At the correct wavelength, it may be possible to
design quantum-well material to approach a 'chirpless' laser where
aH—*0. A detailed study of Henry's alpha factor has been undertaken
by several authors, with a notable discussion by Osinksi and Buus
[34].
f conduction band
% act
n^ ^ tf Q2
3 act 'Q1
valence band
a b
Figure 2.15 Energy-band/distance diagrams
a Bulk material with active region widths over 20 nm so that a
continuum of energy bands is the appropriate approximation rather
than discrete levels
b Quantum-well material with active region widths in 1-10 nm range for
which quantised levels have to be considered
f
density of states
bulk
Figure 2.16 Schematic comparison of density of states in bulk material and quantum-
tuell material
2.7.2 Gain saturation and the need for multiple quantum wells
Many of the benefits of quantum-well materials result from the lower
density of states and the lower carrier densities needed to reach
transparency; however, these features can also be detrimental. A
separate (optical) confinement structure is essential in quantum-well
lasers to ensure that the optical filament is concentrated around the
quantum well(s); the quantum well does give some positive waveguid-
ing, but it is far too weak to confine the mode. However, because of the
nanometre scale of the quantum regions, confinement factors F<0.01
are typical, in contrast to F~0.5 for bulk lasers. Hence the required
gain for quantum wells may need to be nearly two orders of magnitude
larger than the gain in similar lasers using bulk material. This higher
gain leads to a more rapid filling of the available states at a particular
energy, making any increase in gain at that desired energy (i.e. at a
desired wavelength \0) impossible. As the drive current increases
steadily, yet higher energy levels in the wells may be filled, increasing
the gain at shorter and shorter wavelengths, but this will not increase
the gain at Ao. The gain therefore saturates at a given Ao with increasing
electron density, as was pointed out in Section 2.4.3. These serious
effects of gain saturation, especially with 1.5 |xm material, are
mitigated by stacking, say, three to seven quantum wells in an active
Gain, loss and spontaneous emission 69
ih
momentum k
Figure 2.17 Schematic effects of strain
n-contact
i-contact
valence band conduction band
separate
confinement
region
quantum well
layers
separate
confinement
region
p-contact p-contact
valence band conduction band
Figure 2.18 Transport into quantum wells
Energy/depth band diagram showing the lengths dSCH of the separate
confinement zones, designed to confine the photons close to the
quantum wells. These must be kept small to allow electrons and holes to
move rapidly into the quantum wells. Graded impurity composition of
the confining regions can be used to assist in the carrier transport
into the wells. The electric field is particularly low in these separate-
confinement regions, so that carrier transport is often dominated by
diffusion with a limiting time constant rdij-^q d\CH//JL^T which is of the
order of the time taken for holes (with their lower mobility than
electrons) to diffuse across their separate confinement region. The
problem is reduced if the waveguide layers are thin (dSCHis small), and
reduced still further if a graded-composition separate confinement
heterostructure (GRIN-SCH) is used. Here either a gradient or a series
of steps in the composition of the SCH layers can change the
conduction-band edge with depth, thereby assisting the carriers to
transport from the contacts to the wells in much the same way that
varying composition in the base of a transistor can help the transport
of charge carriers through the base. The small-signal dynamics of
carrier transport are discussed in Appendix 4, along with implications
for changes in emission frequency caused by changes of the free
carrier density in the confinement regions. Once this effect is
understood and the confinement lengths are reduced appropriately,
then transport effects can often be neglected, at least for modulation
frequencies up to tens of gigahertz [45].
72 Distributed feedback semiconductor lasers
2.8 Summary
Starting from a background of degree-level knowledge about semi-
conductors and band diagrams, this chapter provides key tutorial
background for understanding photon-electron interactions in semi-
conductors so that one can appreciate how laser materials provide
optical gain and loss, spontaneous and stimulated emission, and how
these quantities change with electron density and with wavelength.
Important discussions include a two-level model for a semiconductor
which gives the key rate equations for photon and electron inter-
actions and includes the full Fermi-Dirac statistics; this progresses to
an illustrative account of determining the gain as a function of
frequency and excitation.
The various material factors which are important in laser models are
discussed. Among these effects are the coupling factor for the
spontaneous emission to the lasing mode and the changes that occur
with gain and nonuniformity (the Petermann factors); and Henry's
spectral linewidth-broadening factor aH, which describes the changes
of refractive index and gain with electron density. Important effects
such as Auger recombination, temperature dependence, strain in
material, quantum-well material and carrier transport are all touched
on to provide the required background for later chapters.
2.9 Bibliography
Some helpful books on materials:
ADACHI, S.: 'Physical properties of III-V semiconductor compounds' (Wiley
1992)
RAZEGHI, M: 'The MOCVD challenge' (Adam Hilger-IOP Publishing,
1989), vol. 1
MURARKA, S.P, and PERKERAR, M.C.: 'Electronic materials: science and
technology' (Academic Press, Boston 1989) (although this is not about
opto-electronic materials, it contains helpful accounts of a wide range of
materials technology)
PEARSALL, T.R: 'GalnAsP alloy semiconductors' (Wiley, Chichester, 1982)
2.10 References
1 SALE, T.E.: 'Gain calculations for strained InGaAs/GaAs quantum wells'
in 'Vertical cavity surface emitting lasers' (Research Studies Press,
Taunton; Wiley, New York, 1995), chap. 3
2 CORZINE, S.W., YAN, R.-H., and COLDREN, L.A.: 'Optical gain in III-V
bulk and quantum well semiconductors' in ZORY, P.S. (Ed.): 'Quantum
well lasers' (Academic Press, San Diego, 1993)
3 COLES, B.R., and CAPLIN, A.D.: 'The electronic structure of solids'
(Arnold, London, 1976), p. 88
4 SHOCKLEY, W.: 'Electrons and holes in semiconductors' (Van Nostrand,
New York, 1950)
5 KITTEL, C: 'Introduction to solid state physics' (J. Wiley, New York,
1966), lsted., pp. 316-317
6 CARROLL, J.E.: 'Physical models of semiconductor devices' (Arnold,
London, 1974), Section 6.5
7 ZORY, P.S. (Ed.): 'Quantum well lasers' (Academic Press, San Diego,
1993)
8 O'REILLY, E.P.: 'Valence band engineering in strained layer structures',
Semicond. Sd. Technol 1989, 4, pp. 121-137
9 BERNARD, M.G.A., and DURAFFOURG, G.: 'Lasers conditions in
semiconductors', Phys. Status Solid., 1961, 1, pp. 699-703
10 CARROLL, J.E.: 'Physical models of semiconductor devices' (Arnold,
London, 1974), Section 3.5
11 KANE, E.O.: 'Thomas-Fermi approach to impure semiconductor band
structure', Phys. Rev., 1963, 131, pp. 79-88
12 STERN, F: 'Calculated spectral dependence in gain in excited GaAs',/.
Appl Phys., 1976, 47, pp. 5382-5386
13 HALPERIN, B.I., and LAX, M.: 'Impurity band tails in the high density
limit. 1: Minimum counting methods', Phys. Rev., 1966, 148, pp. 722-740
14 BRINKMAN, W.F., and LEE, P.A.: 'Coulomb effects on the gain spectrum
of semiconductors', Phys. Rev. Lett, 1973, 33, pp. 237-240
15 BEATTIE, A.R., and LANDSBERG, P.T.: 'Auger effect in semiconductors',
Proc. Roy. Soc. London, 1959, 249, pp. 16-29
16 DUTTA, N.K., and NELSON, R.J.: 'The case for Auger recombination in
In^Ga, As/,.,,/. Appl. Phys., 1982, 53, pp. 74-92
17 BECK, A.H., and AHMED, H.: 'An introduction to physical electronics'
(Arnold, London, 1968), chap. 12
18 CARROLL, J.E.: 'Rate equations in semiconductor electronics' (Cam-
bridge University Press, Cambridge, 1985).
19 HECHT, E.: 'Optics' (Addison-Wesley, 1987), 2nd ed., Section 13.2
20 GOWAR, J.: 'Optical communication systems' (Prentice-Hall, Hemel
Hempstead, 1993), Section 18.2
74 Distributed feedback semiconductor lasers
3.1 Introduction
horizontal
'vertical variation
variation of intensity
of intensity
marised the concepts where one views the light, propagating along the
layers of dielectric, as plane electromagnetic waves (see Appendix 1)
with a cone of directions determining a pencil of optical rays. Each ray
experiences total internal reflection (TIR) at shallow ray angles of
incidence when reaching the outer boundaries of the guide. It is this
TIR which confines the light to one main region which has an
incrementally larger refractive index than the other layers. To keep the
principles clear, only nonmagnetic material, with permeability |xr
sufficiently close to unity, is considered. The complex refractive index
n=nr+jni and the complex relative permittivity er=srr-\-jeH of the
material are found to be related by n2=et.
For a plane wave in a uniform slab of such a material, the wave
propagates as exp{j(cot— (3z)\ exp gzwhere
(/3+jg)2=co2{fjL0sr((o)s0} (3.1)
Then the local field intensity exhibits standing waves with a periodicity
AOT/2 of the form
E*E=u(x, y)*u(x, y) {F*F+R*R+2 \FR\cos(2/30z+2(f>)} (3.5)
At first sight the periodic changes of intensity might be thought to
react with the charge carriers and set up periodic changes to the
optical gain in the active region. However, in general, these standing
waves have little effect because diffusion of the charge carriers along
the length of the guide smooths out and greatly reduces any impressed
periodicity in the electron density. There are exceptions when the
standing waves do have some effect, and these are discussed briefly at
appropriate times. Over one whole optical wavelength the product
term \FR\ cos(2/3oz+2 0) normally averages to have zero effect (but
see Section 6.4.3).
The modal pattern u (x, y) may be normalised to have a peak value
I u(x, y) \pmk=l permitting definition of an effective area for the mode
in the guide, see Figure 3.2, to be given from
\u(x,y)\2dxdy (3.6)
• / /
The optical gain of the material may now be introduced from the
imaginary part n{ (x, y) of the refractive index. The gain is confined to
the active region which occupies an area F si where the electrons and
holes recombine. The density of photons S interacts with the material
gain per unit distance gm only within the active region and F is known
as the confinement factor, given from
(3.8)
A specific value for the confinement factor F is determined by taking
gm to be the peak value of {co wf(x, y)/c] and taking &iSas determined
from eqn. 3.7.
A slab waveguide, which is used to confine plane electromagnetic
waves in one dimension, provides an instructive starting point for
optical waveguiding to determine u(y) and erejj in one plane for a
particular mode. MATLAB programs are used to discover the main
effects of guiding and then an effective refractive-index model is used
for estimating u(x, y) in two dimensions.
plane wave
TM O
angle to the Oz axis so that Oz', the direction of launch of the plane
wave, is not exactly parallel to the Oz axis. This means that the Hy field
is not exactly parallel to the Oy axis and there is therefore a small
component of longitudinal or Hz field (see also Figure Al.l). Provided
that the angle of launching into the slab guide is shallow enough, then
the plane wave is reflected through TIR at the boundaries. Here any
one optical ray can be thought of as performing a slalom down the
guide with components of Ex, H, and Hz defined by the right-handed
set of axes shown in Figure 3.3. In general, there is a distribution, or
pencil, of rays at a number of different shallow angles at any one
position z. Appendix Al discusses plane waves further.
The TM mode is very similar to the TE mode except that it is now
the Hx field which is purely transverse to the junction (i.e. in the
horizontal direction as in Figure 3.36). The angle of launching of the
plane wave means that there has to be an E2 field as well as an E^ field
so that there is a longitudinal component of the E field.
with HTE = (j/ co/UL0jULr) VxErE and JEVM = ( -j/ (os0er) V x HTM
With the axes as shown in the figures, Ex is the key field component for
the TE mode while Hx is the key field component for the TM mode.
The boundary conditions are determined by seeking the field
components which are tangential to the dielectric interfaces, as these
84 Distributed feedback semiconductor lasers
1 dE
Ex and Hz= are continuous for TE fields
j dy
1 r)//
Hx and Ez= - - are continuous for TM fields
dy
One notices that, for the TE fields with nonmagnetic material, where
/zr=l everywhere, the Ex field and its derivative are continuous across
the boundaries. For the TM modes there will be slight discontinuities
in Ey [= -/3HJ(cos0£t)] at boundaries between layers of different
permittivity whereas dEy/dy [ =j/3EJ is continuous.
An important caveat is required. Although interchanging E and H
can be helpful to keep a similar mathematical formalism for the
calculations of field profiles for TE and TMfields,care is needed when
considering energy or power interactions. Optical gain, in either TE or
TM waves, relies on the electricfieldinteracting with those oscillating
electronic dipoles formed by the recombining holes and electrons
which can radiate into the waveguide in the semiconductor material.
When it comes to the exchange of power and energy between
electrons and photons, it is therefore essential to consider the
interaction of the appropriate transverse E fields taken here as Ex for
the TE mode and Ey=(/3Hx/(oe) for the TM mode in the important
gain region.
an effective permittivity seff as in eqns. 3.2 and 3.3 with the fields
propagating over the whole vertical and horizontal extent of the guide
as in eqn. 3.4. The concept of an effective permittivity also can be
approached through a weighted mean relative permittivity. By the
definition of 'guiding', the optical fields have to vanish far enough
away from the main guiding layer. Now consider a solution Ex(y, z) to
the wave equation
{d%(y, z)/dz2} + {d%(y, z)/df) + kler(y)Ex(y, z)=0 (3.10)
2
with kl = co jULoso.
>z) = (3.11)
Then one may define
Using
[d%(y,z)/dy2)dy=0 (3.13)
If the mode is asymmetric such that J^ Ex(y9 z) dy=O, then one has to
consider different limits between, say, yx and y2, values of y for which
dEx(y, z)/dy vanishes so that with the new limits eqn. 3.13 still holds.
The value of eqn. 3.15 lies in showing that the effective refractive index
is a weighted mean, rather than any utility in computation.
For guides with no loss or gain, the real effective refractive index
n(Jjj~verejj always lies between the maximum refractive index of all the
layers and the maximum refractive index of the two outermost layers.
This gives a clear limited area of search for appropriate numerical
solutions. It is instructive to consider the physical reasons for these
limits. Vertical variations exp( — yyy) or exp(yyy) occur with
jy-koV {srej]~ € !) when erejj>sr and such fields are said to be evanescent
To confine the optical power, the fields in both outer layers have to
decay as the magnitude of 3; increases (i.e. be evanescent). However,
when srejj<er the vertical variations are of the form cos (/3yy+(f)) where
/3v = &o v (s.,. — ereff) and can be said to give vertical propagation. If the
fields were to propagate vertically across the whole set of layers, the
power would dissipate away from the guiding region. Equally, the
modes cannot evanesce everywhere. There has to be at least one layer
where vertical 'propagation' occurs with both outer layers giving an
evanescent vertical variation with negligible energy at the far edges of
the guide. These requirements then force the effective real relative
permittivity to lie between the extreme actual values of the real relative
permittivity.
evanescent
transmission
decays away
from active area
confining layer (bottom)
The example shown in Figure 3.7 makes the layer sufficiently wide
so that, at 1.55 |xm wavelength, it supports three modes. The zero-
order mode has no nulls in the Efield,while the first and second order
modes have one and two null intensities, respectively. The E-field
intensities in the TM mode are very similar in this instance to those for
the TE mode but, because of the different boundary conditions, there
are small discontinuities in the field intensities at the dielectric
interfaces. Only the second-order TM mode has been illustrated
because the other orders show hardly any significant difference from
the TE mode in their field profiles for this geometry. The modal
indices are very similar for both TE and TM modes. One notices,
particularly for the higher-order modes, how the fields evanesce most
strongly in the layer of lowest refractive index so that the optical
intensities do not penetrate as far into the bottom layer (index 3.2) as
Principles of modelling guided waves 89
' / Y ^Ss-^
3.1994 eft. index
200/cm 0
3.277
. • 3.1989
// \ \ V"
. 300/cm / TE field \
TE field
intensity
1intensities \
gain/cm 9air\
refractive
index n=3.4
profile
3.1950 p - - - - jprofil^^
^ 1000/cm 1 \g \
I J_ gain '0' 3.2
they penetrate into the top layer (index 3.25). With loss/gain in the
guide, the fields do not have strict nulls although the strong minima/
maxima remain.
guide with a 0.3 [xm-thick gain region would require a much larger
value of gain than for the 2 fxm guide because of the high diffraction
losses. Trying out numerical examples on a multipurpose slab-layer
program can rapidly give a sense of what can and cannot be achieved
in terms of guiding.
3.5 Scaling
Classic early work [9] on three-slab guides showed how to formulate
the guide in terms of normalised parameters which can be useful once
experience has been gained in the required orders of magnitude for
these normalised parameters. Here the required scaling laws are
discussed briefly to show how to change a successful design of a guide
to another wavelength.
Taking some appropriate vertical layer thickness d, a normalised
parameter can be defined:
V(y)=k,dV{sr(y)-srFff} (3.16)
This enables eqn. 3.10 to be written in the form
{d%(y, z)/d(y/d)*HV(y)%(y, z)=0 (3.17)
For a symmetrical three-slab guide where the central thickness is wx = d
with relative permittivity srX, then a V number can be defined from
koav(srl — er0) and this permits scaling of the guide with straightfor-
ward rules. Typical V numbers lie within the range 1-3 to ensure a
clear single 'zero'-order mode. With multilayers, one can see that,
provided that V(y) is maintained with y scaled with d, then the shape
of the field profile has the same shape, at least on the scale of d.
Narrowing the guide by a factor of two requires a doubling of
[\/{sr(y) — £r(,jjM which is determined by the differences of the relative
permittivity.
1 for m= n I
un{x,y)*um(x,y)&x&y=8nnsin where Smn= f (3.21)
0 for m^n I
The value of sin gives the modal equivalent area. The conjugation is
significant because it leads to what Marcuse refers to as being 'energy
orthogonal' [17]. The energy storage can be split into separate modes
with all the consequences of equipartition of modal energy for
quantum statistics which is significant in mode counting for sponta-
neous emission (Appendix 9). There is an orthogonality (without the
conjugation) for guides with strong loss or gain [1] but these are not
true modes when it comes to counting modes for spontaneous
emission, as is highlighted in [18].
1
polar obliquity
intensity factor
n = 3.5
eff TM
-1 1
outer layer
contact
(eqn. 3.21) with respect to any other modal pattern. The MATLAB
programs provided permit the user to explore the validity of this for
any particular slab-guide design by simply changing the material gain
within the guide.
3.11 References
1 COLON, R.E.: 'Field theory of guided waves' (McGrawHill, 1960)
2 ADAMS, M.J.: 'An introduction to optical waveguides' (Wiley, 1981)
3 SNYDER, A.W., and LOVE, J.D.: 'Optical waveguide theory' (Chapman
and Hall, 1983)
4 CHIANG, K.S.: 'Review of numerical and approximate methods for the
modal-analysis of general optical dielectric wave-guides', Opt. Quantum
Electron., 1994, 26, pp. sll3-sl34
5 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (Wiley,
Chichester, 1980), Section 6.4
6 HESS, O., KOCH, S.W., and MOLONEY, J.V.: 'Filamentation and beam-
propagation in broad-area semiconductor-lasers', IEEE J. Quantum
Electron., 1995, 31, pp. 35-43
7 HAUSS, H.A.: 'Waves and fields in optoelectronics' (Prentice Hall,
Englewood Cliffs, 1984), p. 175
8 ROZZI, T.E., and IN'TVELD, G.H.: 'Field and network analysis of
interacting step discontinuities in planar dielectric waveguides', IEEE
Trans., 1979, MTT-27, pp. 303-309
9 KOGELNIK, H., and RAMASWAMY, V.: 'Scaling rules for thin fibre optical
waveguides', AppL Opt, 1974, 8, p. 1857
10 MUNOWITZ, M., and VEZZETTI, D.J.: 'Beam-propagation computations
in one and 2 transverse dimensions', Opt. Comm., 1993, 100, pp. 43-47
11 LEE, J.S., and SHIN, S.Y: 'On the validity of the effective-index method
for rectangular dielectric wave-guides',/. Lightwave Technol, 1993, 11, pp.
1320-1324
12 CHIANG, K. S.: 'Performance of the effective-index method for the
analysis of dielectric wave-guides', Opt Lett, 1991, 16, pp. 714-716
13 SILVESTER, P.P., and FERRARI, R.L.: 'Finite elements for electrical
engineers' (Cambridge University Press, 1996), 3rd ed.
14 KOSHIBA, M., and INOUE, K: 'Simple and efficient finite-element
analysis of microwave and optical wave-guides', IEEE Trans., 1992, MTT-40,
pp. 371-377
15 YOUNG, T.P.: 'CAD tools for optoelectronic subsystems', GECJ. Res., 1994,
11, pp.110-121
16 HOEKSTRA, H.J.W.M.: 'On beam propagation methods for modelling in
integrated optics', Opt. Quantum Electron., 1997, 29, pp. 157-171 (over 70
references)
17 MARCUSE, D.: ' Quantum mechanical explanation of spontaneous
emission K-factor', Electron. Lett, 1982, 18, pp. 820-822
18 RAMO, S., WHINNERY, J.R., and VAN DUZER, X: Fields and waves in
communication electronics' (Wiley, New York, 1994), 3rd. ed., Section
12.4
Chapter 4
Optical energy exchange in guides
4.LI Introduction
There are two basic classical methods of modelling mathematically the
operation of semiconductor lasers. The first method applies the
concepts of photon/electron particle exchange, outlined in Chapter
2, where one discusses the rate of absorption and emission of photons
along with the rate of recombination of holes and electrons, ensuring
at each stage that there is a detailed balance between photon
generation and electron/hole recombination leading to particle
conservation and energy conservation. This is the standard rate-
equation approach which is robust and well researched but can be
difficult to apply when there are strong nonuniformities, and even
more difficult when the phase of the electromagnetic field is
important. For distributed-feedback lasers, both the phase of the field
and nonuniformities are important and so one has to abandon the
photon-rate equation in favour of an approach based on interactions
between electromagnetic fields and the electric dipoles in an active
optical medium. This chapter presents both methods and shows
(mainly in Appendix 5) that they are, indeed, just different aspects of
the same physics of energy conservation and are wholly compatible
with one another.
The electromagnetic-field analysis is essential when the refractive
index/permittivity changes periodically inside the laser. The chapter
then concludes with an analysis of the coupled-mode equations which
determine how a fraction of the forward-travelling field is coupled into
the reverse-travelling field with a medium which has a periodic
permittivity, i.e. the waveguide contains a Bragg grating.
98 Distributed feedback semiconductor lasers
'vertical' /
/ <
; 'horizontal'
A
distance J / y distance
/ ' ' electron .. p^"} photon
N
^ 7 ^ density ^- -,'--
^ t J ^ density S
w^z——^
axi
arear A ^ ^st a n c e ^
•—;
_. /_ y '
\
1A
______L^
* /
/ .
(4.1)
at
Here the concepts outlined in Section 2.3 are followed and that
Optical energy exchange in guides 99
section should be read first. The first term on the right-hand side in
eqn. 4.1 gives the rate of stimulated emission of photons by the
electrons within the laser's active region, as determined by the
material's optical power-gain coefficient Gm(N) per unit time. The
photon density is taken to be uniform so that any standing-wave effects
have been neglected, as in eqn. 3.7. Only those photons within the
region of electronic gain, of volume FT, experience gain so the
confinement factor F appears here. The second term is the optical-
power absorption per unit time as determined by a material/guide
parameter Am. This term accounts for absorption of photons outside
the gain region. [There can be an additional loss term ( — Tb N SY)
caused by IVBA, as discussed in Section 2.7.3, but it is omitted in this
first analysis.]
The photon lifetime rp gives an estimate for the rate of escape of
photons out of the cavity and can be determined in a very rough-and-
ready manner by assuming that a proportion T of the photons get
transmitted at each facet on each pass down the lasing cavity so that a
photon has an approximate effective lifetime rp~L/(Tvg) within the
cavity. More accurate estimates require the use of field calculations
(see eqn. 4.10). In the final term, N/rsp gives the spontaneous emission
caused by the electrons spontaneously recombining with holes, giving
out a random emission of light over a range of frequencies and all
directions. Only those frequencies and directions which can couple to
the lasing mode can participate in the stimulated interaction (see
Section 2.4.1) so that typically there is a coupling factor /3^~10~5<^l.
When the device is lasing strongly, one can often ignore this
spontaneous contribution but it is essential around threshold to excite
lasing action.
Obviously in eqn. 4.1 the volume T cancels but its inclusion at this
stage makes it more clear as to where and why the confinement factor
F occurs. The electron-hole recombination provides the stimulated
emission over the smaller volume F T than the volume Y occupied by
the photons as a whole. The confinement factor is then a convenient
'single-parameter' method of modelling complicated transverse dis-
tributions of photons and electrons on the assumption that these
transverse patterns are fixed by the guiding mechanisms and do not
change significantly with optical intensity or carrier density. That
assumption is not always true and then serious modifications have to
be made. It is usually a good approximation when the photons are
guided through well designed changes of the refractive index in the
material structure rather than being guided by the electronically
100 Distributed feedback semiconductor lasers
dt r q
where the first term on the right-hand side is the counterpart of the
similarly placed term in eqn. 4.1 and gives the rate of removal of
electrons from the interaction region F T caused by stimulated
recombination associated with the material optical power gain. In a
unit volume, N/rr gives the net spontaneous rate of recombination,
representing this for simplicity by a single (carrier-dependent) rate
constant l/r r . The difference between the time constants rr and rsp is
that the latter relates only to radiative recombination giving out
photons that can interact with the lasing field while the former covers
all recombinations including nonradiative recombinations.
This recombination in detail is made up from the following main
factors:
N/ rr=AN+B(NP) + C (N2F) + C"(NP2) (4.3)
The first term on the right-hand side is recombination proportional to
iVand is usually nonradiative; the second term is the bipolar radiative
recombination of electrons, density N, and holes, density P. The last
pair of terms in eqn. 4.3 represent Auger recombination where there
is no radiation given out, as discussed in Chapter 2. The fact that,
Optical energy exchange in guides 101
Link between power gain/loss per unit time and field gain/loss per unit distance
G=2gug Gm=2gmvg Am=2amvm
where vg is group (energy) velocity in material; the factor of 2 arises because
t = GS=vgdS/dz=2vgE dE/dz= (2vgg)E2
under operational conditions, N~P within the laser means that one
can approximately simplify these terms into three major terms:
N/rr~AN+BN2+CN3 (4.4)
The analysis will show shortly that once the laser is lasing, iV remains
roughly constant or clamped and a linear approximation (N/rr) with
r r in the nanosecond range is a useful first approximation.
The final term in eqn. 4.2 is determined by the driving current /
which pushes electrons at a rate I/q into the conduction band within
the interaction region. This drive current is the source of electrons
and, under equilibrium conditions, must balance the net rate of loss of
102 Distributed feedback semicondudor lasers
(4.5)
rp
dt
104 Distributed feedback semiconductor lasers
~) /\^-rGm(N0)\ (4.8)
W/ K J
with
(4.9)
T
P
With N0=Nth9 the net steady-state photon output rate O0=T S0/rp is also
determined from
(4U)
dS,
(4.12)
Optical energy exchange in guides 107
where Gmi = G^A^ with GJn giving the differential material (power)
gain/unit time (Table 4.1). This pair of coupled equations becomes a
classic second-order equation
where
-Wso+i and
Note that the rates of change with time refer to modulation rates
(typically microwave rates up to several gigahertz) and not to rates of
change at optical frequencies (terahertz).
It is now helpful to list the basic results augmented by calculations
from a large-signal analysis outlined at the end of Appendix 4.
(i) Eqn. 4.14 shows that there is a damped resonant response of the
light output to the changes in the electron drive. Figure 4.3
demonstrates that gain saturation at the higher optical powers
[putting Gm—*Gm/ (1 + 65)] has a strong dynamic effect on the damping
of the resonance, and indeed without it one obtains an unrealistically
low damping of the transient. Typically, the damping rate is measured
on the nanosecond scale. The resonance frequency, referred to as the
photon—electron resonant frequency, increases as the square root of the
steady-state light output and takes values measured in the gigahertz
range. This resonance, which provides an upper limit to the useful
bandwidth of direct modulation, can be pushed up into the tens of
gigahertz by using material with high optical gain.
Horizontal carrier diffusion in the plane of the layers, referred to in
Section 1.3.3 but not considered in detail in this work, adds an
additional rate of loss of electrons and changes the details of the
damping of the electron-photon resonance. In general, these carrier
losses increase the damping and reduce the peak magnification in the
AM response (shown in Figure 4.5) with the maximum effect when the
width of the optical mode is approximately the same as the diffusion
length [11,12].
108 Distributed feedback semiconductor lasers
D = 2 . 5 ; a = 0 . 3 ; Q = 1000
normalised light
1
normalised carrier
w density
0
0 500 1000 1500
time / photon lifetime
Figure 4.3 Normalised dynamic light/current characteristics
The electron-photon resonance is shown from a large-signal dynamic
model with and without gain saturation: Q=Tr/TfJ; D is normalised drive;
a=Nlr/Nlh (Appendix 4). The drive current is switched on at 7^=0 . The
parameters are normalised as in Appendix 4. The normalising power
might be -10 mW with the photon lifetime measured in picoseconds.
Because the light has to grow from statistically varying spontaneous
emission, this emission will alter the turn-on delay of the stimulated
output
(ii) The light output (or rate of photon emission) can increase rapidly
only if the gain rises significantly above the threshold value. Similarly,
the light output can decrease only if the gain falls below the threshold
value so that dynamic changes of the light must be accompanied by
changes in the gain. For small signals, these changes of gain and light
are in quadrature (eqns. 4.12 and 4.13 and Figure 4.4) where the peak
gain (peak electron density) gives the greatest rate of increase of light
while the minimum gain gives the greatest rate of decrease of light.
Figure 4.4 also holds the key to the explanation of significant
changes in frequency and spectral-line broadening when one turns
lasers on and off by changing the gain. Recall eqn. 2.42, repeated here
for convenience, relating the changes in the frequency Agf with
changes in that gain:
Optical energy exchange in guides 109
20 r
relative
response (dB)
-20
-1 0 1
log (modulation frequency /FT)
43A Introduction
Whenever the phase of the optical field becomes an important effect
in the physics, straightforward rate equations have to be augmented.
For example, the modes of a Fabry-Perot laser have to be determined
from arguments about the field and the round-trip phase. In a DFB
laser the phase of the field, which is reflected back from the grating in
relation to the phase of the forward-travelling field, is a crucial factor.
A field analysis is therefore essential deriving, from Maxwell's
equations, the propagation of the forward field combined with some
distributed reflection, caused by the grating, of the reverse field into a
forward component. The propagation of the reverse field and the
coupling of the forward field will follow. It is also important to
understand that energy conservation as exhibited by the photon rate
equation of Section 4.1 is also demonstrable from Maxwell's classical
equations. Appendix 5 does this necessary detailed work developing:
(i) the role of electronic polarisation in the optical medium as the
key method of energy exchange between electrons and photons;
(ii) the connection between 'particle' exchange and electromagnetic
energy conservation; and
(iii) the detailed calculations for forward and reverse waves in a
guided medium with gain and periodic refractive indices.
This second part of the chapter presents a more tutorial account
referring, for detail, to Appendix 5.
where the propagation coefficient /30 is identical over the whole (x, y)
cross-section of the guide. Within the permitted propagation ranges
for a single mode, there is a unique /30 for each a)0. The transverse
variation which gives the guided field pattern is given from a
distribution u(x, y) (as in Figure 4.6) normalised so that the modulus
of the peak value is unity to define an effective area si as in eqn. 3.6.
112 Distributed feedback semiconductor lasers
contact
horizontal
'vertical variation
variation of intensity
of intensity
z
Figure 4.6 TE mode
Schematic field patterns guided by a laser structure with the E field
mainly in the Ox direction. In many semiconductor diode lasers, the field
pattern in the x-direction is broader than that in the ^direction.
= (Tgm-aJF (4.19)
dz
Optical energy exchange in guides 113
f l II II ii i i
where the first two terms ensure that, with no gain or loss, the fields
propagate with a fixed wave-packet envelope at the group velocity
which is the velocity at which energy propagates and is slightly lower
than the phase velocity in semiconductor materials. The terms on the
right-hand side of eqn. 4.19 give the net optical field gain using the
notation of Table 4.1. Here the field gain gm is a function of the
electron density N and may also saturate as the photon density S
increases. This field gain is reduced by the confinement factor F and
also by any absorption of the fields outside the gain region denoted by
am. The term (Tgm- am) appears regularly and, as shown in Table 4.1,
will be replaced with the net field gain per unit distance g.
Next one has to include the excitation of the fields by spontaneous
emission so that the advection equation for the forward wave becomes,
for example:
dF 1 dF ^ /A v
—+ =gF+isi)f (4.20)
u
dz vgdt
The detailed evaluation of the spontaneous excitation, written here as
ispfT is discussed in Chapters 2 and 5 and Appendix 9. This emission is
emitted over a broad range of frequencies and angles so that only a
small fraction /3sp couples to the lasing mode in the narrow waveguide,
114 Distributed feedback semiconductor lasers
(4.22)
vg
with Ef=F(t, z) u(x, y) exp{j(a)ot- /3hz)}. This has now allowed the fields
to have a central frequency co0 and at the same time to appear to have
a central propagation coefficient (3h as if the wave operated at the
Bragg frequency.
Because the spontaneous emission is random, the phase factor
exp(— j8z) multiplying this emission is found to be of no con-
sequence. The key random properties of spontaneous emission are
Optical energy exchange in guides 115
(4.23)
(4-24)
paired with eqn. 4.22. Similarly, one can associate a photon density per
unit guide length Sr with this reverse-travelling field and obtain the
equation comparable to eqn. 4.23:
Pnr)=ha>vgSf{r)d (4.26)
It follows that, apart from the dimensions and calculation of the
spontaneous excitation, the advection equations for Pj{r) are the same
asforS /(r) .
emission, and for simplicity ignore any dispersion, i.e. take the
permittivity to be independent of frequency so that phase and group
velocity are the same and take the parameter 5=0 so that (31=0% S^SQ/ULQ,
to give:
d2E d2E
{ { ) (2/3)}=0 (4.31)
Insert the fields of eqn. 4.30 into eqn. 4.31, assuming slowly varying
changes relative to the propagation coefficient or frequency so that
only first-order derivatives in time and space are retained (Appendix 5
shows how this can be done in more detail). The 'slow' variations in
the forward wave F are placed on the left-hand side to give the first-
order (i.e. ignoring terms such as srfl dF/df) rates of change of field:
dF dF
2j(3bu(x, y) -—+2j(o0srrs0fi0u(x, y) —
oz ot
=« - (o20srfl(x, y) cos(2/3hz)s0iLL0Eexp( -ja)ot+jpbz) (4.32)
The terms on the right-hand side couple strongly to the terms on the
left-hand side only if both sides vary with the same spatial frequency.
This is called the phase-matching condition. Any higher-order Fourier
components (HOFCs) such as £OJULOE exp(-ja)ot+fef3bz) change so
rapidly in space that any terms they excite over the first half of a
wavelength will be excited in the opposite phase over the next half
wavelength and the excitation effectively cancels. It is found that
elimination of HOFCs varying as~exp(/l\|/3/,z), N = ±l; ±2 etc., is avery
good approximation, when integrating eqn. 4.31 over a few wave-
lengths. Taking E as in eqn. 4.29 and writing
cos(2(3bz) =j{exp(j2(3bz) +exp( ~~j2/3bz)} (4.33)
the elimination process discussed above selects the slowly varying
envelope of the reverse fields R(z, t) and allows the variations
exp(ja)0t+j/3bz) to be cancelled to give
a>o8rf\(x, y) cos(2/3bz)s0/ui0Eexp( -j(o0t+jfibz)
- 4 o)lsrJl(x, y)eo/uLoR(t, z) u(x, y) (4.34)
Eqn. 4.32 then simplifies into
dF 1 dF !
()~—=--jj(oj0/vg){srjl(x,y)/en}R(t9z)u(x,y) (4.35)
Optical energy exchange in guides 119
where
I u(x, y) 12 dxdy
• / /
(4 38)
'
120 Distribu ted feedback semiconductor lasers
n
gain and guidance region
XIII'llllllllll
Substrate////
1111111 ilia i
The equivalent equation for the reverse wave is found by reversing the
sign of z:
dR 1 dR
(4.39)
With the phase of the grating chosen relative to the spatial origin as in
eqns. 4.29 and 4.31, Kjr-Krj=K. Chapter 5 revisits these equations, gives
an alternative perspective on their physics and highlights also the
meaning of the factory in front of the coupling coefficient K.
effective
t titt I
A =A, m
dielectric
constant
distributed
feedback
into
laser
Figure 4.10 Feedback and scattering properties of first- and second-order gratings
a First-order grating
b Second-order grating
Symmetrical radiation is assumed for second-order gratings here but
reflective layers would modify this assumption
The modified eqn. 4.32 now has an extra term including Emd:
dF dF
2j/3bu(y) •r-+2j(o0erre0iuL0u(y) —
OZ 01
where
and
u(x, y) dxdy
JJ
• / /
with
ll{8rf[(y)/EJu(x,y)v(y)Emddxdy (4.49)
124 Distribu ted feedback semiconductor lasers
t x,y)
— exp( -j/3hy) dy dx/si (4.50)
4.5 Summary
Starting with the classic rate equations which approximate the physics
of a uniform Fabry-Perot laser, several key features about lasing were
discovered:
(i) The round-trip complex gain including the feedback has to be
unity so that both the phase and magnitude of the field gain play
a significant role in the mode selection.
(ii) There is a resonance between the electron and photon energy
exchange which limits the capabilities for modulating a laser
directly.
(iii) There are increases/decreases in the electron density which
decrease/increase the refractive index and increase/decrease
the lasing frequency as the laser is turned on/off. The
magnitude of this significant effect is determined by aH, the
linewidth-broadening factor.
(iv) In a Fabry-Perot laser there are very few mechanisms to keep
the laser lasing in just one single mode, which is the most
desirable condition for modern communication systems.
(v) The necessity to consider both the phase and the gain suggests
that one should consider travelling-wave field equations for a
laser, and the necessary mathematics was set up to permit wave
packets to be considered with independent frequency and
propagation coefficients centred about relatively arbitrary
choices so as to prepare the ground for the final main section.
(vi) The introduction of a periodic grating into the waveguide which
forms the laser introduces a frequency-selective element that is
126 Distributed feedback semiconductor lasers
the main feature of this book—the Bragg grating inside the laser
diode.
(vii) The wave packets are centred about a frequency close to, but not
necessarily identical to, the peak reflection frequency of the
grating, but the central propagation for the wave packet is taken
to be governed by the periodicity of the Bragg grating.
(viii) The field equations coupling the forward and reverse waves are
discussed and methods for finding, from Fourier analysis, the
correct coupling coefficients given the waveguide structure are
briefly outlined and a MATLAB program has been made
available for the reader to explore this aspect further.
Chapters 5 and 6 continue with the physics and design of DFB lasers.
The detailed numerical modelling is left to Chapter 7.
4.6 References
1 STATZ, H., and DE MARS, G.: 'Transients and oscillation pulses in
masers', in TOWNES, C.H. (Ed.): 'Quantum electronics' (Columbia
University Press, New York), 1960, pp. 530-537
2 IKEGAMI, T , and SUEMATSU, Y: 'Direct modulation of semiconductor
junction laser', Electron. Commun. Japan, 1968, 51(3), pp. 51-58
3 PAOLI, T , and RIPPER, J.E.: 'Direct modulation of semiconductor lasers',
Proc. IEE, 1970, 58, pp. 1457-1465
4 LANG, R., and KOBAYASHI, K.: 'Suppression of the relaxation oscillation
in the modulated light output of semiconductor lasers', IEEEJ. Quantum
Electron., 1976, 12, pp. 194-199
5 CARROLL, J.E.: 'Rate equations in semiconductor electronics' (Cam-
bridge University Press, 1990), chap. 7
6 OSINSKI, M., and ADAMS, M.J.: 'Transient time-averaged spectra of
rapidly modulated semiconductor lasers', lEEProc.J, 1985, 128, pp. 34-7
7 KAWAGUCHI, H.: 'Instabilities and non-linearities in laser diodes'
(Artech House, Boston) Section 2.1.7 and Chapter 5
8 ADAMS, M.J., and OSINSKI, M.: 'Longitudinal mode competition in
semiconductor lasers—rate-equations revisited', IEEProc.-I, 1982, 129, pp.
271-274
9 YABLONOVITCH, E.: 'Photonic band-gap crystals', / Phys. Condens.
Matter, 1993, 5, pp. 2443-2460
10 DEMARTINI, F., MARROCCO, M., MATALONI, P., MURRA, D., and
LOUDON, R.: 'Spontaneous and stimulated-emission in the thresholdless
microlaser',/. Opt. Soc. Am. B-Opt. Phys., 1993, 10, pp. 360-381
11 FURUYA, K., SUEMATSU, Y, and HONG, T: 'Reduction of resonance
like peak in direct modulation due to carrier diffusion in injection laser',
Applied Optics, 1978, 17, pp. 1949-1952
12 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (Wiley,
1980) Section 7.3.4, 'Diffusion of damping laser oscillations', pp. 424—
433
Optical energy exchange in guides 127
5.1 Introduction
distance Z
Figure 5.1 Choice of phase position for equal Krj and Kfr
terms (#cry, Kfr) at the centre of symmetry for this type of grating is
observed always to be positive and remains the same if the reference
point moves.
Now consider a forward field F propagating in a waveguide at an
angular frequency co0 with an axial complex propagation coefficient /30.
The complex amplitude of F at %=0 is transformed to amplitude
Fexp( —j/30A+gA) at z=A and similarly a backward wave of amplitude
R incident at z=A is transformed to amplitude R exp( -j/30A+gA) at
z=0. The gain gis the net field gain per unit length making allowance
for losses and confinement factor. The value of the propagation
coefficient may be written as
Po=P>+8 (5.1)
where f3b gives the Bragg wavelength for maximum reflectivity and
therefore, as in Section 4.3.3, 5 is a measure of the detuning of the
optical wave from the central Bragg condition. Allowing for gain,
phase shifts and coupling, one may write, for one Bragg period,
(5.2)
(5.3)
F and R are defined at the reference planes one grating period apart,
while gA and <5A are assumed to be sufficiently small compared with
unity so that the approximations in eqns. 5.2 and 5.3 hold, which for
practical gratings and gains is usually correct. The 'finite-difference'
scheme of eqns. 5.2 and 5.3 can then be rearranged with only a slight
change of notation into the differential form in the steady state:
dF
+—=(g-j8)F+JKR (5.4)
~=(g-j8)R+j*F (5.5)
represents a grating where the gain maxima (or loss minima) coincide
with the maxima in real index (and vice versa if Kgain were negative).
With Kindex>0 and Kgain>0 one may refer to a gain grating 'in phase with'
the index grating (or a loss grating 'in antiphase with' the index
grating). With Kindex>0 and Kgain<0 one refers to a gain grating 'in
antiphase with' the index grating (or a loss grating 'in phase with' the
index grating).
As before, the phase of the reference point can change the phase of
the imaginary component so that it is easy to be confused between loss
and gain in a grating, and it is safest to stick to writing the equations
with a consistent reference point. A pure index grating had a positive
Krf Kfr product and, using similar arguments, a pure gain grating has a
negative K^ Kfr product regardless of the point of phase reference.
n
8= ^-^ (5.10)
c
The two eigenmodes then become
A{exp( ~j/3hz) +5+ exp(j/3bz)} exp( -jftz) (5.11)
B{exp(jM +5_ exp( -jfa)} exp(jftz) (5.12)
where /3, is a propagation coefficient offset from the Bragg value, in the
presence of coupling', /3e is to be found for each 5 and g The solution to
any problem with a uniform grating can then be obtained by matching
the eigenmodes to the boundary conditions. The solutions for s+/_ are
found by substituting eqns. 5.11 and 5.12, and into eqns. 5.4 and 5.5 to
obtain
s _[(8+jg)±{(8+jg)KfrKrfn
«
«rf
Basic principles of lasers with distributed feedback 133
s jv+mwtfK**^ (514)
Kfr
0.5 . g/K
g/K 0 0 8/K 8/K
Figure 5.2 Normalised dependence of real and imaginary parts of fie on 8 and the
field gain gfor a pure index guide
Results are symmetric in 5, and gain has been taken as positive
real (P e / K)
2.5 v
imag (Pe
g/K
Figure 53 Normalised dependence of real and imaginary parts of (3e on 8 and the
field gain gfor a pure-gain grating (real K)
Results are symmetric in 8, and gain has been taken as positive
8/K
phase (radians)
Figure 5.4 Dependence of the real part of fie on detuning and phase (arg Kgmn/Kindex)
of K with Kr>0
/3 can take positive or negative values, and there is symmetry in 8
its ends, are now discussed. In Appendix 6 it is shown how the coupled
equations may be solved analytically and directly with a matrix
formulation where the forward and reverse fields are arranged as a
vector with
M ss]
The fields at any position z can be put in terms of the fields at z=0
(5.16)
through M(z):
M(z) = cos (
1 0 0 ~j
(g-jS) — K (5.17)
0 J 0
1
reflectivity
1.5-.
0.5-
Figure 5.5 Reflectivity of a uniform grating of length KL=2as a function ofgain and
frequency deviation
The reflection is symmetrical in 8 so only one-half of the diagram is drawn
to avoid too much detail
laser. Other devices introduce, for example, two Aw/8 phase shifts
appropriately displaced about the centre of the laser. For such devices,
the lasing line lies within a band of amplified spontaneous emission
which has significantly lower amplitude where the reflection sculp-
tures the emission spectrum. It is this band of low spontaneous
emission which is commonly referred to as the 'stopband' and can be
found by numerical calculations as discussed later.
Figure 5.6 shows the dependence of this extended definition of a
normalised * stopband width' (measured as a change in free-space
wavelength 2AAv/0/;) on the KL product of the grating for:
(i) a uniform grating laser,
(ii) a symmetric Am/4-phase-shifted structure, and
(iii) a DFB laser with 2 x Am/8 phase shifts positioned symmetrically in
the cavity and spaced by 70% of the cavity length (see Section
5.6.5).
The curves are plotted as dimensionless groups so as to generalise
them:
27T 1
^normalised (5.25)
10
8--
A
normalised
6
stop band -
• 2x
•••
width /8
normalised
•«•
A uniform
KL
Figure 5.6 Normalised 'stopband width' against KL
The stopband width is defined here as the separation between the two
main peaks of spontaneous emission at high gain and so varies for lasers
with different structures as discussed in Section 5.5
function of KL means that, for low values of KL, AAV/0/, « 1/L, and for
high KL products, AAv/0/; <* /c, i.e. for very short grating lengths the
'stopband width'-—•infinity because there is no frequency selectivity,
while for high KL products the 'stopband width' <* K as one might
expect [6].
ex
PI (5.26)
140 Distributed feedback semiconductor lasers
The confinement factor V has now been inserted and modifies the
material gain; T had not been introduced at the start of Chapter 1.
This general principle that the round-trip complex gain is unity
applies to a DFB laser as to any laser. The facet reflectivities pkfl and plighl
are replaced by equivalent grating reflectivities with appropriate
phases of reflection, and the actual cavity length has to be changed to
an equivalent cavity length. This is the essential physics which is hidden
in the initial lasing equation (eqn. 5.23).
One might think that the strong reflection at the Bragg frequency
fbmgg would ensure oscillation a t / ^ ^ but eqn. 5.23 showed no solution
at this central frequency for real K but instead indicated modes towards
the edges of the 'stopband'. This paradox arises because, although the
reflection and gain are strong enough at fbmgg, the feedback is in
antiphase which prevents lasing action. This phasing can be under-
stood by considering a microcavity formed by one grating 'tooth'
(taken to be AOT/4 long with a low index value) exactly at the centre of
the device, for example with the index stepping to the higher value at
each end of this tooth. All grating reflections at the left-hand end of
this tooth are in phase with each other and can be lumped into one
real reflection (p,oolh ^ <* K) at the step to the adjacent higher-index
tooth. Likewise, all the reflections on the right-hand end can be
collected into a reflection (p(oothrigh) at the step to the adjacent right-
hand higher-index tooth. From symmetry p(oofhrighl=p\oolh ¥=pml gives the
net reflection at the ends of this central tooth of the pure index
grating. Each single pass is Am/4 long (a phase shift of TT/2) giving the
net round-trip feedback with one pass there and one pass back as
{ionrtexp(-;i7r)}2<0 (5.27)
•t
relative position
L
of 2.5 • 10.735
jump
(a.u.)
--P^- 10.725 • = - ! -
intensity
distance along laser distance along laser distance along laser
400 u m W
spontaneous emission spectrum
schematic structure
Side modes
V4
phase shift
- 400 u m -
laser length
1 565 nm P574
Figure 5.7 Simulated performance ofhigh-KL KJ4 phase-shifted DFB laser operating
at low output power
Output power at each facet 0.37 mW, /cL=3, zero facet reflectivity,
2000 A/cm2 drive current, effective recombination time 4 ns, transpar-
ency density 10 cm"3
18
F(z) =
I ^ UMHA^HMJ I M (5.31)
If 0=0, then the oscillation condition (eqn. 5.23) for the uniform DFB
laser may be recovered. A Am/4-phase-shifted DFB laser is modelled
with 0= (TT/2) or exp(j<f>) -j which is found to yield
Basic principles of lasers with distributed feedback 143
where Mx and M3 are given by eqn. 5.17 for the appropriate section
lengths. An analytic solution for eqn. 5.34 is not straightforward and it
is preferable to compute the results. However, one can see the utility of
transfer matrices for studying complex structures.
Figure 5.8 shows the simulated field and electron-density profiles
along the cavity of a 2xA w /8 laser analysed as a concatenation of
uniform grating sections. The significant point about the 2xA w /8-
phase-shift DFB laser is that it provides a more uniform field and
electron density along its length over a wider range of current inputs
than the single Am/4-phase-shift structure for appropriate KL products
- 2 , provided that the phase shifts are positioned correctly at
appropriate values of AL from the facets as shown. With a flat field
profile, the gain is more effectively utilised, there are narrower line
widths because of small variations in the refractive index along the
structure and there is less variation in the output wavelength with
power level than for many other types of Bragg laser.
While considering transfer matrices, it is worthwhile making a short
digression to find the transfer matrix for N pairs of teeth of a Bragg
grating, at for example its central Bragg frequency. A transfer-matrix
formulation gives
M= (Mw/4M+(i)ApMw/2M_ ( j )Ap M, /4 ) N (5.35)
where a reflection + 2A/0 from a step change in refractive index at one
position is joined to a reflection of - |Ap by a transfer matrix M^ as in
eqn. 5.30 but with a phase shift of 4>= IT/2 corresponding to the optical
length of Am/4. This unit is then completed with transfer matrices of
phase shift TT/4 on either side to complete the half-wavelength-long
unit which concatenates with N other such Bragg elements. The
elemental transfer matrix for a step change of refractive index,
assuming plane-wave fronts, is M+>A/9 given from
144 Distributed feedback semiconductor lasers
intensity
10.676
distance along laser distance along laser distance along laser
Spontaneous emission spectrum t~ L = 300 u. m
schematic structure
AL AL
10
relative " -»'« * •
emission
Intensity
10 C X /8 X /8
phase shift phase shift
nm
1560 1570 250 u r n
separation
|
(5.36)
-Ap/2 1 J
where T=V{1 — (1A/9)2} and §Ap= (n, - n^/{nx + n^). Then looking
ahead to a future section on Bragg reflectors for VCSELs (Section 8.5),
an estimate for the transfer matrix of 20 pairs of teeth with ^=3.5
and ng=3 (very approximately the values for GaAs and AlAs) is given
from
\20
20
-1.0119 -J0.1548
+70.1548 -1.0119
Basic principles of lasers with distributed feedback 145
io 9349
R
i
-
-jlO.8891
> i o - 8 8 9 i ir^i
10.9349 J | ^ J
(5>37)
so that, with zero input on the right-hand side giving RL=0,
\ Ro/Fo\2=\ M2l/ M22\2=0.9916 or 99% power reflectivity. This is an
optimistic calculation, neglecting loss and assuming that the optical
energy is confined fully to the regions containing the dielectric steps,
giving the maximum reflection possible. This type of estimate will be
used later in Section 8.5.
the carrier density N(z) and photon density S(z) along the z direction.
The changes are relatively slow, so that diffusion does not significantly
moderate the effect and diffusion is ignored. The coupling of electron
density and photon density is found from eqn. 4.2 divided by the
volume T. The drive-current density replaces the drive current
because one now has to allow for nonuniformities in the electron
density N(z). The gain Gm(N) is replaced with the differential gain
expanding in N about the transparency density. The linear-to-cubic
recombination terms of eqn. 4.4 are all considered to give
frequency frequency
Bragg frequency.
0.5 S Uniform
DFB
0.06 0.02 0.006 0.010.02 0.04 0.08
0.1 ••
1 2 K L 3
Figure 5.11 Longitudinal-intensity-variance contours for 2 x A J8 DFB lasers
Parameters give variance as measured by k\E\2/\E\2mnx and ordinate
AL/L gives position of phase shifts: note that the values of AL/L give
two distinct modal points; see Figure 5.13
Uniform
DFB
58 62 66 70
W AAAAAAAAA/ W i
i
| I J^^<^T ,A 1/1=0.5: n A uniform DFB
\J
* I = 400 u m >jS?
/ I =0:^/4 DFB O g off Bragg'
U
^ 1 • I
^ , ^ ^ ^ • *
A/=—7 (5.39)
Pch
where p( is the resistivity and tc is the thickness of the contact layer. For
jfr-InP-doped material, the resistivity pr( = l/qPl-ihoU) ~8x 10~4 Urn for a
hole density P~lx 1018 cm" 3 and mobility /^o/,~75 cm 2 /Vs [14].
Hence for a contact layer of such material with t(.= 2 |xm and
A/~1.5kA/cm 2 we find that AV^~25meV. The variations of the
injected carrier density along the laser are determined by the
variations in the stimulated photon rate. To avoid significant spatial-
hole burning in such a laser, then, it is necessary for these variations to
be much less than 1.5 kA/cm2. Typical threshold-current densities are
of this order of magnitude, so that spatial-hole burning is likely when
one is operating at two or more times threshold.
If the device was fabricated on a p-type substrate, it would be
appropriate to consider the voltage drop arising from injection
through n-type material. A typical value of the resistivity of n-lnP
doped at lxl0 1 8 cm~ 3 is about 3xl0~ 5 dm [15] which would imply
that longitudinal hole burning could be significantly reduced as a
result of the low series resistance for such a device.
A more detailed analysis involves relating the junction voltage to the
injected carrier density along the laser, taking into account the
resistivity and thickness of the intervening layers, and the voltage
difference between the contact and junction. The relationship
between junction voltage and carrier density could be obtained using
the conventional Fermi-Dirac functions for the occupation probabil-
ity, but this approach requires an integration to obtain the total carrier
density. An alternative is to use the approximate but usefully accurate
analytic relationship between Fermi-level energy and carrier density
for electrons and holes, respectively [16], defining contact parameters
Acm and Bcon from
^-^(1+0.154J -^~
V l\
and
154 Distributed feedback semiconductor lasers
V - V.
con
R
J
\7
(one segment of laser) (one segment of laser)
where %g is the bandgap energy. The carrier-rate (eqn. 5.38), now with
dA//d£=0 in the steady state and an external applied contact voltage of
Vam, is then modified to become
_/ = (V(:on VJ)J)
= I-N(r)S(z) (5.43)
qd pt(qd
F F F
o M
M+l M+2 P-1 p
R p.
fields fields
section numbers
Figure 5.15 Subdivided laser for transfer-matrix methods
Electron density iVis uniform in each subsection
(5.44)
At the facets, the appropriate field reflectivities pright and pkfl give
conditions RL=prighlFL a n d F0=pleftR().
The method is iterative and, if there are large numbers of significant
modes, then the iteration takes longer. It is therefore suited to DFB
lasers where there are perhaps only three or four significant spectral
modes which need to be tracked. The iteration starts by ignoring
stimulated emission and calculating the carrier density in each
subsection, allowing for the appropriate carrier-recombination terms
(as in eqn. 4.2 but for static conditions). The effective refractive index
of each subsection is then calculated allowing for the plasma effect and
Henry's alpha factor (Chapter 2) [18] which relates the change in the
real and imaginary refractive index with carrier density:
Basic principles of lasers with distributed feedback 157
aH =—
(5.45)
(5.46)
2TT
Once the carrier density and refractive index have been found in each
subsection, the round-trip gain and phase, taking some reference
point near the centre of the laser, are computed as a function of
frequency using eqn. 5.44. From these transfer-matrix operations the
forward and reverse fields at each frequency can be computed along
the laser and from these fields the mean photon density at any
frequency in each section may be found.
The peaks in the gain spectrum are the potential lasing modes. All
the potential modes then have to be examined carefully to see which
of these modes are at or above their own threshold gain for lasing. For
each mode that is above its threshold gain, it will subsequently increase
in power and therefore its associated forward and backward field
amplitudes, at the modal frequency, in each subsection are stored
along with the total photon density over all (say, ML) such modes. The
frequencies of these modes around threshold are also stored along
with the photon densities in each section. Using this model, as
discussed later in Section 6.2, the linewidth is 'Lorentzian' and this fact
permits one to fit the solutions from a few frequencies around a
potential mode so as to find the net optical power, thereby reducing
the calculation time. Finding this photon density is an iterative process
in itself, adjusting the carrier density in each subsection to allow for
the mean stimulated emission rate for a self-consistent solution. Note
that the refractive index is not altered at this stage because this would
change the frequency of the mode significantly. The change in gain
introduces only a small frequency shift.
158 Distributed feedback semiconductor lasers
0.8 i
F>ower output, mW waveicmgth shift, nm
measured
measured
— O -
4
" O
0.6
/ / computed - 400 urn
2 c/ ' computed
A/4 phase shifted DFE
0 / ' \ 0.4 I
C) drive current, mA 10
° (5 power output, mW 5
a h
Figure 5.16 Static characteristics
For high KL ( = 3) A/4 phase-shifted DFB showing curvature of the
light/current characteristics, along with the proportional shift of the
emission wavelength across the stopband
The next step is to find the electron density in each subsection from
the steady-state carrier-rate equation, eqn. 4.2 (see also Appendix 4)
modified to take into account the different recombination mecha-
nisms and the total sum of all the photon densities SM in all the ML
modes with different differential gains GmM and transparency densities
NlrM. The electron rate equation in the steady state now reads as
L G; M {MZ)-^M}5M(Z) (5.47)
M=l
1 2 M N
Figure 5.17 Spontaneous fields in Mth section (reflection and transmission)
(5.50)
For the output from the left-hand facet, the subscripts / and r are
interchanged.
It can be seen that, if there is zero spontaneous emission, then a
finite output requires 11 — pMrpM/l2=0. The fact that (PMTPM/) is a
regular complex function of frequency means that this is zero only at
discrete real frequencies which define the permissible lasing fre-
quencies. With finite spontaneous emission, one discovers, around the
lasing frequency a)0, that II -pMrpMf\ = l - GMr-j 2(co- (o0) YMr where
GMr and YMr are parameters associated with the round-trip gain (see
Basic principles of lasers with distributed feedback 161
IAAA
1/10 1/10 I I I j I
1575 158O 1575 1580
wavelength, nm wavelength, nm
also Section 6.2) as seen from the centre of the Mth section and (o0 is
the same 'lasing' frequency for all values of M. This form of the
denominator, as a function of frequency, ensures that the line shape of
any semiconductor laser is effectively a Lorentzian line shape within
the limitations of the model; and fitting a few frequencies to this shape
permits the prediction of the output around the lasing frequency,
without the need to compute the precise maximum output frequency.
Moreover, it is found that, if 11 - pMrpMj I2=O for the Mth cell, then it is
also zero for the M'th cell and the calculation of the lasing frequency
is independent of which cell is chosen as the 'central' cell but is
determined by the total effect of all the sections. However, the
closeness of the denominator to zero in eqn. 5.50 sometimes leads to
errors unless the local gain is iterated appropriately. Figure 5.18 shows
typical emission spectra calculated by this method, which has also been
formalised into a technique referred to as the power-matrix method
[24]. Figure 5.166 shows the spectral relative shifts within the stop
band for a AM/4 phase shifted DFB.
5.9 Summary
The chapter started with a discussion of the physics of the grating-
coupling process, a section which should be read in conjunction with
work in Sections 1.6 and 4.4. This section gave alternative insights into
162 Distributed feedback semiconductor lasers
5.10 References
1 KOGELNIK, H., and SHANK, C.V.: 'Coupled-wave theory of distributed
feedback lasers',/ Appl Phys., 1972, 43, pp. 2327-2335
2 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (John
Wiley 8c Sons, 1980), pp. 480-481
3 WANG, S.: 'Principles of distributed feedback and distributed Bragg-
reflector lasers', IEEEJ. Quantum Electron., 1974, 10, pp. 413-427
4 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (John
Wiley & Sons, 1980), pp. 501-502
Basic principles of lasers with distributed feedback 163
6.1 Introduction
6.2 Linewidth
6.2.1 General
The spectral channels in wireless-communication systems are deter-
mined by precise radio or microwave carrier frequencies with the
spacing between these channels determined by the linewidth of the
modulated carrier. At radio frequencies, the modulation bandwidth
dominates this linewidth so that one can say (for the amplitude
modulation format at a maximum frequency of fM and allowing for the
upper and lower sidebands) that the modulated carrier linewidth is
close to 2/Jvf. The principle remains with other modulation tech-
niques—the modulated linewidth is determined by the modulation
and its format. Channel spacing follows with appropriate spectral
'guard bands' to ensure that negligible interference from neighbour-
ing channels can statistically occur. The DFB laser, with its ability to
determine the optical frequency, can define channels in an optical
communication system so that its linewidth and the accuracy of
specifying its frequency become important system parameters in
determining minimum optical channel spacing for a WDM system.
Unfortunately, even when the modulation frequencies fM~0 the
linewidth of a DFB laser is significant and cannot be ignored. But
worse, as the modulation frequency increases then the carrier's
linewidth broadens much more than the value of 2fM that would be
considered for the classical radio system mentioned above. The laser
exhibits frequency shifts called dynamic chirp because of changes in
refractive index caused by changes of gain as the laser is modulated, as
outlined in Sections 2.5 and 4.2 (Figure 4.4). On switching the laser on
and off there are changes in gain, modulated by the photon-electron
resonance; hence the laser is frequency-modulated at the photon-
electron resonance frequencies by an amount which depends on the
More advanced distributed feedback laser design 167
sidemode
a m p.i,ude
Figure 6.1 Schematic spectrum of modulated microwave oscillator and diode laser
where Pnetsp = 7T(o01 TEsptot I 2 /2Q gives the net spontaneous emission over
a spontaneous emission-noise bandwidth -TTTCOQ/2Q in this model.
Hence the linewidth may also be written as
(6.8)
oul) (6.12)
The generality of eqns. 6.10 and 6.11, which go back to the roots of
electromagnetic theory, forms a basis to estimate or measure the
effective value of aH for arbitrary laser structures. Consider a change
in output intensity which requires a small change 8Ground in the round-
trip gain of the laser; then this has to be accompanied by a phase
change to give a total round-trip change of complex gain of
8Ground(l+jaH). The oscillation/lasing condition (neglecting changes
in the spontaneous emission) is therefore effectively
E()U((co) [(1 - Ground - 5G r a w w J+72Q{(w- OJ0) / co0}-j8GroundaH]
-~(l+jaH)TEspiol=0 (6.13)
From eqn. 6.13, the change in intensity causes a change in the central
frequency of Sfo/fo = — aH(8Grmnd/2 Q). Further, because the intrinsic
3 dB linewidth is (A/VK/) = &>0(l — Ground)/27rQJ this also has to change by
8(AIwf) = — a)QSGmind/27TQ. The effective value of aH therefore may be
determined from 8f0 and <5(A/W/) to give
This value usually will not be the same as the material's intrinsic
linewidth-enhancement factor given in eqn. 6.8 but is modified by the
field distributions and waveguiding within the laser.
To use eqn. 6.14 to find numerically the effective linewidth-
enhancement factor, first calculate the steady-state solution (for
example as in Figure 6.3a for a AOT/4 phase-shifted laser with KL>2).
The mode intensity everywhere is, say, enhanced by ±1% and the new
charge carrier densities N+ (z) and N_ (z) are recalculated from, say,
eqn. 5.38. These then change the local refractive index n+ (z) and
n_ (z) which in turn leads to new central frequencies fo+ or fQ_ where
there is unity complex round-trip gain. The two linewidths for these
two conditions are calculated so that the net change in the central
frequency, given from 8fo=fo+ ~fo- , along with the net change in the
full bandwidth S(ALWf) =A, w / + — &Lwf- then give the effective value of
&Heijfromeqn. 6.14.
This approach has been used successfully to predict the effective
linewidth-enhancement factor for 2xA w /8 DFB laser structures with
varying KL product and varying phase-shift position AL/L. Figure 6.4
plots the effective linewidth-enhancement factor for a uniform
(AL/L=4) and Aw/4-phase-shifted (AL/L=0) DFB laser as a function
of output power and the intrinsic linewidth-enhancement factor of the
material. The DFB structure significantly increases the effective aH at
174 Distributed feedback semiconductor lasers
amplified con
spontaneous
intensity co_
Aco
frequency
b
super-linear increase
in power at facets
low power levels where a change in the power level creates more
change in spatial-hole burning than at high power levels. With the
phase-shifted DFB the differences between the effective and intrinsic
values of aH are small at both high and low powers, with larger
differences mainly at intermediate powers.
Another striking example of structural effects on the effective value
of aH is given by a DFB laser when there is a complex coupling
coefficient Kind(,x+JKgain. For such a grating with the notation used in
this book, the peaks in the gain and the refractive-index periodicity
coincide when Kindex and Kgain are both positive ('in phase'). With aH>\
an increase in the carrier density increases the strength of the gain
grating, but weakens the index grating by a greater amount and
reduces the round-trip feedback or round-trip gain. The net differ-
ential round-trip gain with respect to carrier density is reduced thereby
increasing the effective linewidth-enhancement factor. Alternatively,
with Kindrx>0 but K^^ <0 then both the gain and index gratings increase
in strength giving a lower cavity loss. The net differential mode gain is
enhanced thereby reducing the linewidth-enhancement factor and
hence the linewidth.
A more sophisticated analysis of the result given in eqn. 6.14 has
been given for a Fabry-Perot laser [17]. Provided that one uses the
effective aH, eqn. 6.14 still holds regardless of the laser being a DFB or
FP. However, care is needed if one attempts to estimate the linewidth
More advanced distributed feedback laser design 175
effective value of a
H
----- uniform DFB
14 A/4 phase shifted DFB
12
-\ ^ material value of, a TT
10 - '* **^ ^ n
S
** //^j; _- 10
8 "~ - p
iiij r Tr
6 "" ^ -^7r.
!l"^!U:-^z_=--=— e
4
4
2 2 .-
0 "i i i I I i i l I 11 I l l l I I I I I I I I i i i i i 11 I i j I
0 1 2 3 4
power output/facet, mW
Figure 6.4 Estimated effective linewidth-enhancement values for uniform and k/4
phase-shifted DFB lasers for different values of material linewidth-
enhancement values
Laser is 400 |xm long /with K=75 cm"1. Results at 'zero' power agree with
Amann [16] who assumed no spatial-hole burning
for a DFB laser on a similar basis to the calculation for a FP laser: the
DFB can have half an order of magnitude larger linewidths. As
outlined in Section 2.4.2, when using standard rate equations, it was
necessary to modify the spontaneous-emission terms by the Petermann
factors, Klr to account for enhanced spontaneous emission trans-
versely, and Kz to account for enhanced effects of spontaneous
emission because of longitudinal distributions of optical intensity. The
net effect is a spontaneous-emission factor Ktof = KlrxKz>l. The
reasons for a significant Kz in certain types of DFB lie in these having
end sections acting more as amplifiers rather than as feedback
elements. The factor Klnl is then very approximately a multiplier of the
spontaneous emission that would be obtained in an ideal uniform FP
laser of similar material and length. A similar factor, Kir, which applies
to gain guiding that can occur horizontally under the injecting contact
of a laser structure [18] where also there is an increased spontaneous-
emission rate into the lasing mode leading to a larger linewidth-power
product.
176 Distributed feedback semiconductor lasers
unity, do not clamp the gain and, in such modes, the amplitude
fluctuations induced by the spontaneous emission give them a
relatively broad linewidth. At high enough powers, the negatively
cross-correlated noise in the main and side modes, combined with a
difference in their nonlinear gain coefficients, give rise to rebroaden-
ing of the main lasing mode. A detailed analysis [19] then predicts an
extra linewidth contribution ALW+f which varies with the cube of the
sidemode intensity PSM, and as the square of the difference Ae in the
nonlinear gain coefficients e for the main and side modes:
P
^ = [M,] [MJ [M3] [M+] [M,] I ^ I (6.16)
+M22
P (6.17)
relative frequency
facet 1: p h a s e 8 0
(degrees) 2.9 3
stop band width (nm)
coated front facet with a random phase with respect to the grating
(corresponding to most manufacturing conditions); allowed for
longitudinal-mode spatial-hole burning and examined the main mode
and principal side mode. Values of KL-0.75 optimised the single-mode
yield for high power [31] and the rear-facet phase had a strong
influence on the maximum output power of this single mode. Later
work [32] examined the effects of a small reflectivity for the front
facet, showing that in general this reduced the single-mode yield for
K L ~ 0 . 7 5 , but that this yield was nearly independent of KL provided
that the front-facet reflectivity <5%. Threshold current and linewidth-
power product are also influenced by facet reflections [33].
Using the previous range of facet reflectivities, with phases stepped
at 45° intervals, the performance of a 2 x Am/8 DFB laser was examined
when modulated with a 40 mA peak-to-peak current drive at 2.5 Gbit/s
and with a bias level selected to give an extinction ratio of about 10:1,
values appropriate to real systems. The simulation considered the
main mode and principal side mode together with spacial hole
burning. Figure 6.6a shows the distribution of threshold currents
obtained and Figure 6.66 shows the relative frequencies, for given
values of sidemode-suppression ratio, as the facet phase is stepped.
The simulations above took into account submount parasitic
impedances to produce realistic current drives into the laser at this
frequency. This device was modelled with parameters given in Table
A11.2 and had significant speed limitations because of the time taken
182 Distributed feedback semiconductor lasers
relative
frequency
10.
10
6.4,1 General
The concept that deep periodic etching into the active-gain region can
give a gain grating has already been met. In terms of Fresnel plane-
wave reflections at each step, as outlined in Section 1.7.1, with a
complex refractive index change between adjacent teeth of a grating,
the reflectivity tip per Bragg period (two reflections per period) in
Section 1.7.1 becomes:
)
)
However, this gives a significant overestimate per step because the
deep etch does not give an abrupt junction penetrating fully into the
optical fields but rather a modulation of the thickness of an active
region or an absorbing layer, as discussed next. The numerical
techniques such as those outlined in the program slabexec in
directory slab have to be used, in general.
g
A wrA r
15
K L
gain
,0.06
1
index
A relative intensity
10 3 T K L * 0.2+/2
10
1568 15701572
wavelength, nm
A relative intensity relative intensity |
10
optical carriers over the access network with cable television is another
clear application for lasers which give out more power than the
conventional 'communication' source. In WDM applications, power
combiners for the different wavelengths usually lose significant
amounts of power, and again higher power is required.
There is, then, a strong commercial interest in high-speed low-
wavelength chirp sources which are composed of narrow linewidth CW
DFB-laser sources coupled monolithically, or in a hybrid manner, to
external modulators. The modelling of devices with external mod-
ulators is left until Chapter 8 and for the present it is noted that this is
a tried and trusted system designer's choice for controlled low-chirp
modulation, though it adds complexity and cost. The attenuation in
such external modulators can be significant, and in hybrid arrange-
ments there is the additional loss associated with the coupling between
the laser and modulator, and consequently the source needs perhaps
10 dB more power than previously.
One way of increasing the power output from a laser is to change
from using a symmetric device, which usually wastes power by radiating
equally from both facets. A high front-to-back emission ratio can
concentrate the power where required but this leads to surprising and
significant drawbacks.
laser drive
booster amplifier drive
^anti-reflection
grating quantum wells
coated facet
barriers
Figure 6.10 Schematic monolithically integrated DFB laser and flared amplifier
Slightly modified from the drawing by MJ. Fice, Nortel [57] with
permission
0.2 _
0.1 -
output power, mW
5 "
wavelength, nm
1551
1551
wavelength, nm
1552J
1552.
1.E0-J
mode power, mW
1.E-2-
time, ns
output power
5 _
mW
time, ns
energy, pj
10 -3dB: 0.0997nms .
1 -10dB: 0.296nm
0.1 -20dB: 0.382nm
0.01 50dB SMSR
-30dB: 0.439nm
0.001
0.0001 long wavelength
sidemode^
r 1 1 i i 1
1551.2 1551.6 1552.0 1552.4 1552.8
emission wavelength, nm
(simulated spectrometer resolution = 0.1 n m )
Figure 6.12 Simulated performance of 2x\m/8 DFB laser with significant delay
(-70 ps) in injecting carriers into active regions [79]
(a) Total intensity against time (b) Wavelength of each mode against
time (c) Power in each mode against time (d) Unfiltered transmitted eye
diagram (e) Time-averaged spectrum including simulated spectrometer
resolution
More advanced distributed feedback laser design 199
6.9
100 200 300
distance, \im
Figure 6.13 Longitudinal hole burning along laser's length changing with time
As illustrated by plot of active-region carrier-density N profile against
time. Device is switched from a 'zero' at 1.1 of threshold current to a
'one' level with 40 mA more current
5 -
1 time, ns
wavelength, nm
1552.6 J
lasing mode
1552.4-
time, ns b -
wavelength, nm
Jong wavelength
1553.4 -
side mode
1553.2-
1 2 time, ns
output power, mW
J\ J\
1
time, ns
-3dB 0.0921 nm
-10dB 0.169 nm
-20dB 0..243 nm
-30dB 0..299 nm 51dBSMSR
long wavelength
6.7 Summary
This chapter began by looking at how to estimate the linewidth starting
from a theory which applied to most electromagnetic oscillators but
was then found to be inadequate for semiconductor lasers because of
the plasma and Kramers-Kronig effects where changes in the gain
(changes in electron density) give changes in the refractive index
through aH, Henry's linewidth-broadening factor. This increases both
the static linewidth and the dynamic linewidth and in the latter case
has detrimental consequences in increasing the dispersion penalty for
high-bit-rate communication over long-haul optical fibres. The effec-
tive value of the linewidth enhancement as measured by the output
spectra from a device was noted often to be larger than the intrinsic
linewidth-enhancement factor of the material and this was explained
More advanced distributed feedback laser design 203
by considering the averaging effects over the laser. The dynamic effects
of linewidth enhancement have been noted in Chapter 4.
Perfectly uniform DFB lasers need modifications to be acceptable to
ensure that the DFB operates only in the longer-wavelength mode of
its two main modes. It is only this mode which is found to have an
acceptable performance within an optical-communication system,
taking into account SMSR, eye-opening etc. One of these mode-
selection methods was seen in Chapter 5 to be the introduction of
2x Am/8 phase shifts into the grating, but even so the studies discussed
in this chapter show that operation is critically affected by the strength
and phase of the facet reflections, and reflections should be reduced
to below 0.1% for a reliable commercial product; and by the same
arguments one has to protect the laser carefully from reflections
coming from the system if the best performance is to be obtained.
An alternative change to the uniform DFB laser discussed in this
chapter is to introduce some gain coupling into the grating. This latter
can be accomplished, for example, by etching the grating through to
the gain region thereby giving some periodic gain as well as periodic
changes in refractive index. Provided that the gain coupling was 'in
phase' with the index coupling, then the desired longer-wavelength
mode could be selected. There is evidence that such complex grating
lasers may be more tolerant to reflections if the imaginary coupling
coefficient is large enough.
There is a need for more optical output power, especially for
pumping EDFAs and for use in WDM systems. An increase of at least
10 dB can be obtained by using integrated amplifiers and, when
combined with an integrated modulator, this combination is often a
'preferred' source for designers requiring high-bit-rate modulation of
DFB lasers with the capability of being used in a WDM system.
If perfectly uniform lasers were to operate consistently in the desired
long-wavelength mode, then studies discussed in this chapter show that
the phase of the facet relative to the periodic structure of the grating
has to be controlled to a fraction of a radian and, with normal
manufacturing tolerances, this would give unacceptable yields. The
2x Am/8-phase-shifted DFB laser is much better in selecting the
preferred longer-wavelength mode when appropriately designed. Part
of this design when the lasers are to be directly modulated at gigabit-
per-second rates comes in making thin enough separate confinement
heterojunction layers. The final section on modelling showed the
importance of reducing the time taken to transport the carriers across
these layers to well below the period of the modulation.
204 Distributed feedback semiconductor lasers
6.8 References
1 WHITEAWAY, J.E.A., THOMPSON, G.H.B., COLLAR, A.J., and ARMIS-
TEAD, C.J.: T h e design and assessment of A/4 phase shifted DFB laser
structures', IEEEJ. Quantum Electron., 1989, 25, pp. 1261-1279
2 GHAFOURI SHIRAZ, H., and LO, B.S.K.: 'Structural dependence of
three phase shift distributed feedback semiconductor laser diodes at
threshold using the transfer matrix method (TMM)', Semicond. Sci.
TechnoL, 1994, 9, pp. 1126-1132
3 ZHOU, P., and LEE, G.S.: 'Mode selection and spatial hole burning
suppression of a chirped grating distributed feedback laser', Appl. Phys.
Lett, 1990, 56, pp. 1400-1402
4 OKAI, M., TSUCHIYA, T., UOMI, K., CHINONE, N., and HARADA, T.:
'Corrugation-pitch-modulated MQW DFB laser with narrow linewidth
(170 kHz), IEEE Photonics TechnoL Lett, 1990, 2, pp. 529-530
5 WHITEAWAY, J.E.A., GARRETT, B., THOMPSON, G.H.B., COLLAR, A.J.,
ARMISTEAD, C.J., and FICE, M.J.: T h e static and dynamic characteristics
of single and multiple phase-shifted DFB laser structures', IEEEJ. Quantum
Electron., 1992, 28, pp. 1277-1293
6 KOGELNIK, H., and SHANK, C.V.: 'Coupled-wave theory of distributed
feedback lasers',/. Appl. Phys., 1972, 43, pp. 2327-2335
7 FLANIGAN, B.J., and CARROLL, J.E.: 'Mode selection in complex-
coupled DFB lasers', Electron. Lett, 1995, 31, pp. 977-979
8 CHAN, Y.C., PREMARATNE, M., and LOWERY, A.J.: 'Semiconductor laser
linewidth from the transmission-line laser model', ZEE Proc.—Optoelec-
tronics, 1997, 144,(4), pp. 246-252
9 OHTSU, M.: 'Highly coherent semiconductor lasers' (Artech House,
Boston, 1992)
10 HENRY, C.H.: 'Theory of the linewidth of semiconductor lasers', IEEEJ.
Quantum Electron., 1982, 18, pp. 259-264
11 HENRY, C.H.: 'Theory of phase noise and power spectrum in a single
mode injection laser', IEEE J Quantum Electron., 1983, 19, pp. 1391-1397
12 HENRY, C.H.: 'Phase noise in semiconductor lasers',/. Lightwave TechnoL,
1986,4, pp. 289-311
13 PROBERT, P.J., and CARROLL, J.E.: 'Lumped circuit model prediction of
linewidth of Fabry-Perot and DFB lasers, including external cavity
devices', Proc. IEE.-J., 1989, 134, pp. 22-32
14 SLATER, J.C.: 'Microwave electronics' (Van Nostrand, 1950; Dover, New
York, 1969)
15 SCHAWLOW, A.L., and TOWNES, C.H.: 'Infra-red and optical masers',
Phys. Rev., 1959, 112, pp. 940-949
16 AMANN, M.C.: 'Linewidth enhancement in distributed-feedback semi-
conductor lasers', Electron. Lett, 1990, 26, pp. 569-571
17 LU, H., BLAAUW, C , BENYON, B., LI, G.P., and MAKINO, T: 'High
power and high speed performance of 1.3 |xm strained MQW gain
coupled DFB lasers', IEEEJ. Quantum Electron., 1995, 31, pp. 375-381
18 PETERMANN, K.: 'Calculated spontaneous emission factor for double-
heterostructure injection lasers with gain-induced waveguiding', ZEEE/.
Quantum Electron., 1979, 15, pp. 566-570
More advanced distributed feedback laser design 205
37 TSANG, W.T., CHOA, F.S., WU, M.C., CHEN, Y.K., LOGAN, R.A.,
TANBUN-EK, T., CHU, S.N.G., SERGENT, A.M., MAGILL, P., REICH-
MANN, K., and BURRUS, C.A.: 'Gain-coupled long wavelength
InGaAsP/InP distributed feedback lasers with quantum well gratings
grown by chemical beam epitaxy'. Proceedings of 13th IEEE International
Semiconductor Laser Conference, 1992, paper B-l, pp. 12-13
38 ZOZ, J., and BORCHERT, B.: 'Dynamic behaviour of complex-coupled
DFB lasers with in-phase absorptive grating', Electron. Lett., 1994, 30, pp.
39-40
39 LU, H., MAKINO T., and LI, G.P.: 'Dynamic properties of partly gain-
coupled 1.55 |xm strained DFB lasers', IEEEf. Quantum Electron., 1995, 31,
pp. 1443-1450
40 MORTHIER, G., VANKWIKELBERGE, P., DAVID, K., and BAETS, R.:
'Improved performance of AR-coated DFB lasers by the introduction of
gain coupling', IEEE Photonics Technol Lett, 1990, 2, pp. 170-172
41 DAVID, K., BUUS, J., and BAETS, R.G.: 'Basic analysis of AR-coated, partly
gain-coupled DFB lasers: the standing wave effect', IEEE f. Quantum
Electron., 1992, 28, pp. 427-433
42 MORTHIER, G., and VANKWIKELBERGE, P.: 'Handbook of distributed
feedback laser diodes' (Artech House Boston, 1997), chap. 3
43 LOWERY, A.J.: 'Large-signal effective a factor of complex-coupled DFB
semiconductor lasers', Electron. Lett, 1992, 28, pp. 2295-2297
44 LOWERY, A.J., and NOVAK, D.: 'Enhanced maximum intrinsic modula-
tion bandwidth of complex-coupled DFB semiconductor lasers', Electron.
Lett, 1993, 29, pp. 461-463
45 PAN, X., TROMBORG, B., OLESEN, H., and LASSEN, H.E.: 'Effective
linewidth enhancement factor and spontaneous emission rate of DFB
lasers with gain coupling', IEEE Photonics Technol. Lett, 1992, 4, pp.
1213-1215
46 LI, G.P., MAKINO, T , MOORE, R., and PUETZ, N.: '1.55 |jim index/gain
coupled DFB lasers with strained layer multiquanturn-well active grating',
Electron. Lett, 1992, 28, pp. 1726-1727
47 LI, G.P., MAKINO, T, MOORE, R., PUETZ, N., LEONG, K-W., and LU,
H.: 'Partly gain-coupled 1.55 |xm strained-layer multiquantum-well DFB
lasers', IEEEf. Quantum Electron., 1993, 29, pp. 1736-1742
48 MAKINO, T , and GLINSKI, J.: 'Effects of radiation loss on the
performance of second-order DFB semiconductor lasers', IEEE]. Quantum
Electron., 1988, 24, pp. 73-82
49 BAETS, R.G., DAVID, K, and MORTHIER, G.: 'On the distinctive features
of gain coupled DFB lasers and DFB lasers with second order grating',
IEEEJ. Quantum Electron., 1993, 29, pp. 1792-1798
50 DAVID, K, MORTHIER, G., VANKWIKELBERGE, P., and BAETS, R.:
'Yield analysis of non-AR-coated DFB lasers with combined index and gain
coupling', Electron. Lett, 1990, 26, pp. 238-239
51 ADAMS, D.M., CASSIDY, D.T., and BRUCE, D.M.: 'Scanning photo-
luminescence technique to determine the phase of the grating at the
facets of gain-coupled DFB's', IEEE f. Quantum Electron., 1996, 32, pp.
1237-1242
More advanced distributed feedback laser design 207
7.1 Introduction
A commonly used starting point for modelling lasers is a small-signal
analysis of the rate equations of the form discussed in Section 4.1
where perturbations from the steady state are considered [1-3]. These
solutions are often explicitly analytic and/or they can be rapidly
computed. Appendix 4 gives an outline of such an analysis including
carrier transport [4], from contacts to the radiative recombination
region, which is not included in this chapter. Small-signal methods
help to elucidate the physics of modulation and noise [5-8], especially
around steady-state values which can be computed more readily than
large-signal dynamic states. For Fabry-Perot lasers, coupling of
electron equations and photon equations has been done in a variety of
ways well reviewed by Buus [9] but the power of computers has moved
far in the last decade. DFB lasers have more complex structures and
have led to new methods specifically to aid in this understanding.
Transfer-matrix techniques, mentioned in Chapter 5, are also known as
transmission matrices [10,11] and are used for the analysis and design
of multisection and nonuniform lasers by tracking their performance
around specific frequencies. However, random spontaneous inputs to
a laser give randomly varying outputs which then need averaging. The
power-matrix method [12] is a transfer-matrix method which was
specially developed to compute the ramn-square values of optical
power even though the laser is excited by a stochastic 'spontaneous'
input. As has already been seen from Chapter 5, transfer-matrix
techniques operate well for 'single-mode' lasers with only one or two
clear weak sidemodes which can be all tracked as the laser changes
frequency, but particular care is needed with many modes or where
spectral mixing occurs. The techniques can be specially helpful for
long or complicated structures where time-domain analysis requires
excessive numerical storage.
210 Distributed feedback semiconductor lasers
This chapter aims to provide tutorial material for the reader with
limited experience in numerical modelling. The philosophy is to use
straightforward low-order techniques where the synergy between the
physics and numerical methods is reasonably good. Speed of computa-
tion is gained from modern high clock rates and cheap random-access
memory rather than the sophistication of high-order computational
techniques with large step sizes. The reader is referred to texts and
handbooks [21-24] rather than the research literature. Sections on
ordinary differential equations and hyperbolic partial differential
equations can then be selected. To help build up the reader's
confidence, the work starts with first-order ordinary differential
equations before moving in graded steps to the coupled travelling-
wave equations of the DFB. A series of tutorial MATLAB programs,
able to run on the student version of MATLAB [25], is provided via the
'net'. The reader with a full copy of MATLAB 4.0 or better can
enhance these programs, interfacing with C + + to gain computational
speed if required. In laser diodes, short enough pulses of optical
energy can disperse and change their shape through nonlinear
physical interactions, and it is important that the programs do not
confuse numerically induced gain or distortion with similar physical
effects. The good numerical analyst would design algorithms with step
lengths determined by error control [26] built into the algorithm, but
the ready availability of high-speed computation has seduced the
authors into more limited programs, checking that there are no
singularities at small step lengths, and then comparing results with
different step lengths to determine that a satisfactory value was used.
12)
Z Z+l
Figure 7.1 Forward, central, and reverse differences (3, 2, 2, respectively)
dF/dz—{^(Zs+ss) -F(Zs-k)}/s
(7.5)
The half step may suggest that one is introducing an additional layer
of points but this is simply an initial convenient notation to ensure
symmetry between the +ss displacement and the — £S displacement to
give a net step of s. The symmetry with Lax averaging proves to be
important in the (spatial or temporal) frequency domain because a
finite-difference solution 'samples' the electromagnetic fields leading
to the so-called 'Nyquist' limitations on frequency [28] which occur in
any discrete signal. This will be discussed later. The half-space intervals
may be 'removed' through shifting the origin by writing Z+2—>Z' + 1
and Z - 5—>Z' and then dropping the prime to give
l
(7.6)
7.2.2 Accuracy
The accuracy of these finite-difference schemes in the end depends on
the mathematical result [29] that for all values of the complex 'gain'
y=(g+j8)L
(7.7)
Lim^00[(l+ir/N)N/(l-lr/N)N]=expr (7.8)
Programs stepr and stepj in directory diff are available to help the
reader gain a 'feel' for the problem of convergence as the number of
steps N is changed. The programs integrate dF/dz=zF (Figure 7.2) and
dF/dz=jzF, respectively, from z=0-3 starting with F=l. The test
examples of either exponential growth or pure imaginary rate of
change in F(z) have been chosen because these (or mixtures) are
common in laser physics. It will be found in stepj that eqn. 7.5 is
slightly modified so as to ensure d(F*F)/dz=0 demonstrating the
considerable benefits of appropriate central-difference schemes which
can avoid numerically induced growth when there is no physical
growth.
214 Distributed feedback semiconductor lasers
df/dz=fz
with step length
f=l at z = 0.1
50
central
40
forward
30
20 analytic
f=exp( 1 / 2 z 2 )
10
0
0 1 2
Figure 7.2 Comparison of difference schemes
7.3.1 Introduction
The two equations
1 dF dF 1 dR dR
+—= 0 —=0 (7.9)
vg dt dz Vg dt dz
are the first-order wave equations for waves propagating in either the
forward or reverse directions with group velocities ±vg. In the
literature of partial differential equations, these are often referred to
as the advection equations. Their physical solutions are straightfor-
wardly given from
F=f(t-z/vg) and R=r(t+z/v) (7.10)
where f(t) and r(t) are arbitrary functions of time. It is possible to
consider immediately the equations with gain and phase changes
varying as
1 dF dF 1 dR dR
_ +— (7.11)
Vg dt dz v dt dz
because, on putting
Numerical modelling for DFB lasers 215
+—= 0 - —= 0 (7.13)
vg dt dz vg dt dz
(7.14)
Low values of a imply that the pulse travels more slowly while larger
values of a make the pulse travel faster. Notice the use of 'Lax
averaging' in space where
F(T4,Z)^{F(T+lz4)+F(T+l Z-\)\ (7.15)
with similar Lax averaging in time. Different arrangements of Lax
averaging can give related finite-difference schemes that can then be
explored for stability and accuracy. Eqn. 7.14 can be rearranged with
future time values of F given in terms of previous values:
|,Z-|) (7.16)
The half-space values may be 'removed' by a shift of the time and space
216 Distributed feedback semiconductor lasers
b
Figure 7.3 Plots of pulse propagation
Figure 7.3 shows how pulses change their shape in space as they
propagate. The time of arrival at each point in space depends on the
group velocity, and, as the group velocity changes, the numerical
methods introduce some distortion. With differences of group velocity of
a few per cent less than the normalised value of 1, the distortion is usually
acceptable. Here by having a 50% variation in the normalised group
velocity one can see significant distortion induced at the pulse edges; this
could be reduced by having shorter steps giving more steps but less net
distortion over the same total distance
(7.21)
where again when a = 1 and the half spaces are 'removed', as explained
after eqn. 7.16, to give
i?(T+l,Z)=/?(T,Z+l) (7.22)
The requirement of a constant group velocity is not too serious a
limitation. It can be relaxed easily if the group velocity does not
change with time or frequency but changes only with space. Uniform
time steps 8t are then maintained but the spatial steps may be adjusted
to correspond to vg{z) 8t at each local space point Only if vg(t, z)
varies with time as well as space is it necessary to consider a more
complicated system, and this is not done in this book.
1 dF dF 1 dR dR
J J
vg dt dz vg dt dz
(7.29)
f,Z-|)}
(7.30)
1 -J\KS\ \F\(T+\),
i 1 JL |
1 JZKS \\ r\(J
I 72ACS 1 \ \ R ( T
Write
cos jd
(7.32)
with tan \6=\ /cs and eqn. 7.32 then gives a matrix U{6) such that
iflid) U(0) = l wheref implies conjugation and transposing, f/is said to
be a unitary matrix. Changing the origin of time and space as before to
'remove' the half steps in eqn. 7.31 one can get back almost to the
conventional advection relationships linking ^{(T+l), (Z+l)} and
J?{(T+1), Z} with F{T, Z} and R{T, (Z+l)}:
where sinh 0r = Kgain s/(l — 4 K2ain s2) the gain-coupling coefficient. Now
one finds thatF*Fz+1 — i?*i?z is correctly conserved, again supporting
the strength of the central-difference scheme combined with Lax
averaging in obtaining a useful approximation to the physics.
(7 37)
'
(7 38)
-
1 dF dF 1 dR dR
_„+___ = y#c/j — =jKp (739)
J J ;
vg dt dz vg dt dz
(7.40)
where sin 0=KS/(1+4~K S ). Frequently in a DFB, one finds that there
2 2
are long sections with a uniform K SO that a series of space steps can all
have the same values of ft
The stability criterion should be revisited because now, at each space
step, the forward field changes by exp(j$Z) and with /3~0 an increase
at each time step is given from exp{(g—jS)s}. This growth at each step
is the correct physics and does not indicate numerical instability.
Numerical instability would arise if the rate of increase in time became
unbounded or unphysical as the value of 0 either increased or took
critical values. Here, following the work leading to eqn. 7.36, the
dispersion relationship is
cos(d)~sS-jsg)=cos 0cos 0 (7.41)
The gain remains bounded for the high spatial frequencies and one
may expect convergence and adequately low distortion for sufficiently
small steps of length s.
13.2 MATLABcode
A follow-up to eqn. 7.40 must show how to concatenate the spatial steps
of length s. Initially, the field numbers are considered from 1 to N f 1
where there are N sections making the laser's length L=N>. For the
moment, the laser is considered to be uniform. The 'present' fields
(forward and reverse), labelled as ff(n) and fr(n) where n refers to the
space-step numbering, are stored as vectors ff and fr. The 'new' fields,
created at the next time step, are stored as ffn and frn .The 'complex
magnification' per step is
mL=exp{(gL-j*dL)/N}-(2*N + gL-j # dL)/(2*N-gL+rdL) (7.42)
Here gL — j*dL is the code for (g~~ jS)L with L the laser length. The last
term is the central-difference approximation to the exponential. This
magnification is stored as an array mL(n) defining gain and phase
changes at each spatial step and is updated at each time step as
required. From eqn. 7.33 , the Bragg grating introduces a reflection
term sin 8 { = KS/(1+4K2S2) ~KS} and a transmission term
cos 0~ v (1 — K2S2). Numerical 'energy' conservation requires
2 2
cos 6 4-sin 6=1 so, although rounding sin 6 to ACS only gives a small
224 Distributed feedback semiconductor lasers
(7.49)
R(32) = 0
10 20 30
distance
Figure 7.4 Numerical solution for a uniform-DFB-structure amplifier
KL=2; gL=0.5; 8L=4 with a 31-section laser (32 field points) calculation.
On the scale of this drawing one cannot distinguish the difference
between the analytic and the numerical solutions from dfbamp
a(z)+d(z) 1
(7.50)
ib(L)/[a(L)-b(L)}
dF 1 dF
-r+-—=gF+ispf(t,z) (7.51)
Hl
dz vg dt
Excitation of the fields is through the random spontaneous emission,
represented by the term isp(t, z) whose magnitude is studied in
Nu men cal modelling for DFB lasers 227
power spectrum
coherent + coherent +
white noise field
filtered noise field
imaginary
(a) random walk
(c)
-100 frequency 100
filtered power spectrum
20 •
dB
'instantaneous (d)
instantaneous field
-20
-40
-60 real
-100 frequency 100
Figure 7.5 Complex coherentfieldwith added noise
A program spontd provides numerical demonstrations of adding sponta-
neous emission (e.g. white noise) to a coherent signal (i.e. a constant
complex field) at the central frequency. The demonstration shows
(a) A display of the sampled instantaneous fields with time
(b) The white-frequency-domain spectrum apart from the coherent
spike
(c) The effect in the time domain of filtering in the spectral domain
(discussed in Section 7.9) to give
(d) A coherent signal plus narrowband noise where the coherent vector
performs a random walk in the phase plane as in (c)
F1 F2 F3 F4
t
o
N1 ! N2 ! A/3
1 2 3 4
• distance
Figure 7.6 Field and spontaneous excitations
At different space points, for one time
4
temporal spectral
amplitude amplitude
Nyquist limits
—nvo +Jcv
g
spectral
temporal amplitude d
amplitude
l
ill,
time angular frequency
step st/Vg step 2n Vg /N
Figure 7.7 Discretisation: schematic relation between time step and spectral limits
A main period
temporal (a)
amplitude
l1,
time-
main period
4 1
: ' (c)
spectral
amplitude
, 1l Jjili I Ilk 11,
angular frequen
aliasing aliasing
Figure 7.8 Discretisation: periodicity and aliasing
= fl/n(+i)+a/n(-i) (7.60)
Writing K=(2r vg/s), and changing to unit time shifts, this filter
rearranges to give
Numerical modelling for DFB lasers 235
(7.61)
This program is given in filti under the directory filter. The theory can
be tested out analytically by noting that, in the normalised frequency
domain, there is a phase shift of exp(jd)M) per whole time step (i.e.
(bM-(x)Ms/vg) and then use of eqn. 7.60 with half steps gives the
spectral form
BOu(=A[n cos(j(bM)/{co${jd)M)+jKsin(Td)M)} (7.62)
As already noted, Lax averaging replaces aIn with [aIn (+2) + % ( - 2M/2
as in eqn. 7.60 and filters the signal as cos(2 coM) with the Nyquist
limitations — 7T<coM< TT. The parameter fallows the curvature of the
gain with frequency in the central part of the filter to be changed so
that, if JK=1, there is a spectral power filter varying as I cos(2 (x>M) '2- If
0 < K< 1 the spectrum is flattened with respect to the situation with
K=\ while K> 1 narrows the gain bandwidth, increasing the gain
curvature. The filter is stable provided that K> 0.
It is frequently more useful to have the spectral filter offset from the
central laser frequency by an amount co0gsel so that (oM—*coM — (o0jj-sel:
I B(hu 12= \Aln I V [1 + {(coM - (o<gM) r ) 2 ] (7.63)
The implementation of this offset frequency in the conventional
frequency domain replaces exp(ja)Mt) with exp[j(a)M — (*)Ojfsel)i\. In the
normalised frequency domain with one positive space step, there is a
phase shift of exp( — j(oojpe(). This phase shift can be incorporated into
the digital filter by writing @ = Gxp(jd)0^sel) and then putting
(a) -5 00
-5
-10 -10
-15 -15
-5
-10
-15
the noise power falls to half of its central value. Figure 7.9c shows the
output from the same filter as given in Figure 1.9b but now
implemented with a white-noise input (using the program filtin).
On examining the result of Figure 1.9 c and considering this applied
to gain, the reader may be disturbed at seeing the 'output' rise again
at the Nyquist-frequency edges. This problem, if it is a problem, is
removed by an alternative filtering demonstrated in the programs filt2.
In this alternative program, the gain always falls to zero at the Nyquist-
band edges but the gain profile is not symmetrical about the gain peak.
Numerical modelling for DFB lasers 237
7.9.1 General
Care is needed in implementing the gain filter of Section 7.8 because
the field gain over a distance z is given from Gjidd =exp(j gdz) while the
power gain over the same distance is given from Gpoxm.=exp(2 \ gdz).
The integrals are not straightforward when the gain varies with time.
Now the physics requires the gain g per unit length to be some
function of frequency so that the temporal filter must affect g, as in
Section 7.8, rather than act on GfieU or Gpower, where there would be a
more complicated spectral action. Note that, even if g((o) = g(co0), a
real value at co=(o0f then on moving away from this central frequency,
the Kramers-Kronig relationships (Appendix 7) force the gain to take
complex values g((o)—+grmi(co)+jgimag(a))9 where the magnitude of
gimng((o) depends on the local curvature of grml(a>). However, it is only
gmd that changes the power gain per step; the imaginary component
merely changes the phase change per step. It was for this reason that,
in Figure 7.9a and b, the bandwidth was related only to the real part of
the gain.
A further complication comes from the fact that, as the gain
increases with an increase in the electron density, so, at the same time,
there is an associated reduction in the permittivity of the material
leading to a complex differential gain g(l+jaH) which changes with
frequency and electron density. Here aH is known as Henry's
linewidth-broadening factor [36, 37], taking values typically between 1
and 7 depending on the material and structure of the lasing waveguide
[38]. However, it has to be recognised that the net complex refractive
index is made up of several contributions with each contribution
varying in distinctive ways. Appendix 7 addresses this issue with three
238 Distributed feedback semiconductor lasers
components: gain, loss and the plasma effect. At the central gain peak
g(o)0) at angular frequency OJ0 it is changes in the plasma effect (or
oscillation of a plasma of nearly free electrons) to which is attributed
changes in the phase term jg((o0)aH, with the value at the trans-
parency density taken as the reference level. Over the frequencies
where the gain is strong, this phase contribution is assumed to remain
approximately constant. The Lorentzian filter is then applied only to
the net-gain/loss term and not to the term in aH:
-ii}0)} (7.65)
with
grml(a))=g(co0)/{l + i*(w-<*0)2} (7.66)
and
^(ai--a) o ) 2 } (7.67)
For each spatial step of length s, the complex magnification mL of eqn.
7.42 is given by exp{g(&>) (l+jaH)s} and is now modelled as
mL=((1+Tgs)/(1-ygs))*exp(j*ahs) (7.68)
where the MATLAB vector ahs is the array of values of {g(<o0)aHs} that
are allowed to vary at each step dependent on the central gain
associated with the electron density at that step, but are taken to be
independent of the frequency of lasing. The real part of this
magnification at the central frequency is given by exp{g(a)) s} per step
represented in eqn. 7.68 by its central difference approximation of
( l + | g s ) / ( l - | g s ) as in eqn. 7.42 (with the substitution that L/N = s)
but now as the frequency changes, g(co) s, or equivalently the term gs,
will be modelled in accordance with the filter principles just outlined.
At the time of writing, the definitive paper in this area of complex gain
modelling with frequency has perhaps not yet been written, and more
needs to be done numerically, experimentally and physically, but [35]
provides important ideas for transmission-line laser modelling.
where
Bo(n + 1) = {ff(n + 1)-Co(n)} with Ao(n) = {ff(n + 1) + Co(n)}
with
Co(n) = exp{j ahs(n)}{cost(n)ffo(n)+j sint(n)fro(n + 1)}
Eqn. 7.71 is then unravelled, retaining the assumption that the gain
and phase changes caused by the electron density vary sufficiently
slowly, to give the result
ffn(n + 1) = ma(n)C(n>-mb(n)Co(n) + mc(n)ff(n + 1) (7.72)
where
ma = [{K+1+Tg(n)s}/{K+1-yg(n)s}] (7.73)
mb=Ep[{K-1-yg(n)s}/{K+1-Tg(n)s}] (7.74)
mc=Ep[{K-1+yg(n)s}/{K+1-ig(n)s}] (7.75)
The updated reverse field frn(n) is found in a similar fashion and
gives
frn(nr) = ma(nr)D(nr)-mb(nr)Do(nr) + mc(nr)fr(nr) (7.76)
where
D(nr) = exp{j ahs(nr)}{cost(nr)fr(nr + 1)+j sint(nr)ff(nr)}
and
Do(nr) = exp{j ahs(nr)}{cost(nr)fro(nr+1)+j sint(nr)ffo(nr)}
Equations 7.72-7.76 achieve the objectives of giving the new fields
{ffn(n) frn(n)} at each time point in terms of the present fields {ff(n)
fr(n)} and old fields {ffo(n) fro(n)}. The gain filter now is a function of
modulation frequency and has been designed to give zero gain at the
Nyquist-band edges, adjusted for any frequency offset. This sig-
nificantly reduces the risks of aliasing errors or unphysical
oscillations.
20 frequ. offset 1
-10
gain gives strong attenuation at the Nyquist limits unless there are
exceptional frequency shifts. Further shaping of the digital filter could
be accomplished by making K complex as, for example, when using
the program filt2 in the directory filter.
L =j{KL) R+ [n {gmL)
\ V Yt+fz) ~j
-j(8L)\-(aL)]R+(Lisl,r) (7.78)
(F*F+R*R\ I
)7 (7.79)
vg ) qL
where
N ^
+ gM
r vr/
vgrr
(7.8!)
vgrm [q(NfrTML)}
•601
0 100 200 300 400 500 600 700 800 900 1000 1100 1200 1300 10 15 20 25 30
field section number
a b c
normalised power spectrum v. frequency. GHz normalised electron density v. dtetance photon/s (solid sum, dotted forward, dashed reverse)
1.5
Figure 7.11 Typical output for dfbH (a DFB with one central KJ4 phase shift)
a initial run; b continuation; c field profile with spontaneous emission;
d spectrum; e electron density-distance; / field profile without
spontaneous emission
7.11 Summary
The chapter started with discussions on the numerical solution of an
ordinary differential equation with complex rates of growth, high-
lighting the importance of the straightforward central-difference
method. The advection equations or travelling-wave equations were
considered, focusing on a central-difference technique with appro-
priate steps in space and time. Gain and phase changes were
introduced along with coupled reflections. At each stage, numerical
stability was considered and model test programs given. A uniform
Bragg-grating amplifier was considered and the numerical model of a
forward travelling wave inside a Bragg grating was implemented in
MATLAB and checked against an analytic model.
The processes of spontaneous emission are modelled by uncorre-
lated Gaussian excitations injected into each section in the laser, and
tutorial models of the spontaneous-emission processes were made
available in MATLAB code.
The process of discretisation leads to limitations, and these were
briefly discussed, particularly in relation to the fast Fourier transform.
Strategies for Lorentzian spectral filters were outlined and then
implemented to ensure that the gain inside the laser varied with
frequency in an appropriate manner. The use of stochastic excitation
to investigate broadband behaviour of a filter was demonstrated.
Finally, a basic program for a DFB laser operating around the peak
gain with spontaneous emission was discussed and implemented in
MATLAB code.
It is a pleasure to acknowledge that many of the ideas in this chapter
owe a considerable debt to research students Chi Tsang [39] who
participated in the European COST-240 exercise [40] , Mark Nowell
[41] and Barry Flanigan [42], Si Fung Yu [43] and David Jones [44]
and especially to research associates Li Ming Zhang [43,45] (see also
[12]) and Dominique Marcenac [45-47]. Their work has been pulled
together with a more consistent use of central-difference methods for
laser-diode DFB analysis than appears to have been previously
Numerical modelling for DFB lasers 249
7.12 References
1 AGRAWAL, G.R, and DUTTA, N.K.: 'Semiconductor lasers' (Van Nos-
trand Reinhold, New York, 1993), 2nd ed., Section 6.6.1
2 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (J. Wiley,
Chichester, 1980), Section 7.3.2
3 CARROLL, J.E.: 'Rate equations in semiconductor electronics' (Cam-
bridge University Press, Cambridge, 1985), Section 6.4.2
4 NAGARAJAN, R., ISHIKAWA, M., FUKUSHIMA, T., GEELS, R.S., BOW-
ERS, J.E., and COLD REN, LA.: 'High speed quantum well lasers and
carrier transport effects', IEEE J. Quantum Electron., 1992, 28, pp.
1990-2007
5 OHTSU, M.: 'Highly coherent semiconductor lasers' (Artech House,
Boston, 1991), chap. 2
6 PETERMANN, K.: 'Laser diode modulation and noise' (Kluwer Academic
Publishers, Dordrecht, 1991)
7 VANKWILEBERGE, P., MORTHIER, G., and BAETS, R. 'CLADISS—a
longitudinal multimode model for the analysis of static and stochastic
behaviour of diode lasers with distributed feedback', IEEE J. Quantum
Electron., 1988, 24, pp. 2160-2169
8 MORTHIER, G., and VANKWIKELBERGE, P.: 'Handbook of distributed
feedback lasers' (Artech, Boston, 1997)
9 BUUS, J.: 'Principles of semiconductor laser modelling', IEE Proc.—-J.,
1985, 132, pp. 42-51
10 COLDREN, LA., and CORZINE, S.W.: 'Diode lasers and photonic
integrated circuits (j. Wiley, New York, 1995), Section 3.5.2
11 BJORK, G., and NILSSON, O.: 'A new and efficient numerical matrix
theory of complicated laser structures. Properties of asymmetric phase
shifted DFB-laser'/ Lightwave Technol, 1987, 5, pp. 140-146
12 ZHANG, L.M., and CARROLL, J.E.: 'Large signal dynamic model of the
DFB laser', IEEEJ. Quantum Electron., 1992, 28, pp. 604-611
13 LOWERY, A.J.: Transmission line modelling of semiconductor lasers: the
transmission-line laser model', Int.J. Numer. Modell, 1989, 2, pp. 249-265
14 LOWERY, A.J.: 'A qualitative comparison between two semiconductor
laser amplifier equivalent circuit models', IEEE.J. Quantum Electron., 1990,
26, pp.1369-1375
15 LOWERY, A.J., GURNEY, P.C.R., WANG, X-H., NGUYEN, L.V.T., CHAN,
Y.C., and PREMARATNE, M.: 'Time-domain simulation of photonic
devices, circuits and systems'. Proceedings of Physics and Simulation of
Optoelectronic Devices, SPIE Photonics West, 1996, 2693, pp. 624-635
250 Distributed feedback semiconductor lasers
8.1 Introduction
Computer-aided design (CAD) has been a vital initial ingredient in the
realisation of very large-scale integrated (VLSI) electronic circuits.
Now in optoelectronics projects, mathematical modelling, although at
present lacking the sophistication of CAD in VLSI, is being recognised
as a key component in the design of devices and systems [1]. One
thrust of the book has been to lay out the physical- and mathematical-
modelling techniques for the electromagnetic and electronic
interactions within distributed feedback lasers to give better explana-
tions and to facilitate new designs. This chapter outlines further areas
in optoelectronics where such modelling will be at the forefront in the
design of more effective prototypes. Optical systems where arrays of
devices may be interconnected are one such important area and are
discussed in Section 8.2. Novel concepts such as that of the push-pull
laser (Section 8.3), tunable lasers with Bragg gratings (Section 8.4),
and surface-emitting lasers (Section 8.5) are all candidate areas for
applications of the modelling principles given in this text. The future
for mathematical modelling, coupled with a sound physical under-
standing, is bright, extensive and assured.
(a) optical spectra and waveforms with selected time windows and, if
required, effects of dispersion would be included;
(b) baseband spectra and filtered waveforms using RC, Butterworth or
Bessel filters, typically with 1-7 poles;
(c) eye diagrams giving a qualitative indication of the system
performance;
(d) instantaneous frequency against time (from rate of change of
phase of complex field);
256 Distributed feedback semiconductor lasers
(e) instantaneous and mean powers over selected time windows; and
(/) data-transfer techniques for concatenating simulations and trans-
fer to other computer systems and output hardware.
To find wide acceptance within any research/development organisa-
tion, such model libraries need to ensure that the software is very 'user
friendly' and does not need significant prior knowledge of program-
ming languages or complex physics. The difficult problem of the user
interface is not addressed here but programs like MATLAB 5 are
coming to the rescue with more user-friendly methods for data input
and output.
In this section, results are presented from:
(i) the modelling of SOA post-amplification of output from high-
speed digitally modulated optical sources;
(ii) DFB laser integrated with electroabsorption modulator allowing
for residual reflections from modulator-output facet;
(iii) cross-gain and four-wave-mixing wavelength conversion in a
semiconductor optical amplifier; and
(iv) cross-phase wavelength conversion in a Mach-Zehnder inter-
ferometer structure incorporating two semiconductor optical
amplifiers.
10-
0.0005-
5 -
0
-0.05 0.05 0.10 0.15 -0.05 0 0.05 O10 O15
relative time, ns relative time, ns
a
b d
Figure 8.1 Simulated amplifier performance for 1 mW peak optical input power
4a
20.
-40.
filter detector
relative time, ns
Figure 8.2 Simulated amplifier performance for 5 mW peak optical input power
T 5 4
relative time, ns relative time, ns
10.
5.
1 2 3 4 5 1 \ I Y 5
relative time, ns relative time, ns
Figure 83 Simulation of effect of reflections from output facet of an integrated DFB
laser and electroabsorption modulator on dynamic performance of source
Time step for the computation is 0.0956 ps while the curves are plotted
every 4.11 ps
1 -
power, mW power, mW
10(h
0.2 1559 nm
50-
0.1-
1557 nm
1 2 3 4 3 4
time, ns time, ns
power, mW
power, mW
1556 nm
1558 nm
0.15- I
il
i—. i—
J 1
50- 0.1-
0.05-
^—I
1 2 3 4 1 2 3 4 5
time, ns time, ns
ideal load
semiconductor
optical
amplifier output
ideal optical
transmitter waveform
semiconductor
CW optical waveform optical
at 1557nm amplifier
ideal optical
transmitter
Mach Zehnder with SOA
in each arm ! Modulated input optical
waveform at 1558 nm
Figure 8.7 Schematic diagram of cross-phase ivavelength-conversion system
power, mW
radians
amplified probe waveform wavelength offset 1557 nm
5 -
10- 0—
b
0.02-
0.01-
0.00
-0.1
relative time, ns relative time, ns
-V
Figure 8.9 Emitter-coupled amplifier giving push-pull output
The output from T, which drives T2 is the emitter E. Consequently, T, is
driven as an emitter follower, while T2 is driven as a grounded-base
transistor. The (near) constant current I ensures that increment 81 in the
current flowing in T, is a decrement 81 in the current flowing in T2.
Hence outputs 0, and O2 are in antiphase or push-pull. The configura-
tion is also known as a 'long-tail pair'
Figure 8.10 Stopband and energy/distance diagrams for push-pull DFB operation
a Right output high; b Right output low
Consider a symmetrical split-contact DFB laser as in Figure 8.10. The
modulation of this laser is achieved by increasing IR and simultaneously
decreasing IL so that the stored energy moves from the left to the right.
The digital output changes from low (zero) to high (one). To change
back to the low state, one decreases IR and simultaneously increases IL.
From the symmetry, one can see that an output train of information Q
at the right-hand facet is matched by a complementary output Q at the
left-hand facet. Again, from the symmetry it follows that the steady
frequency in the 'zero' state must be identical to that in the 'one' state.
The dynamic chirp or change of frequency is limited to the transient
region where the energy actually moves, and it will be seen shortly that
there are ways of optimising this transient change.
To understand how the push-pull laser works without delving into
deep mathematics, use can be made of the stopband/distance
diagrams introduced in Section 5.6.2 and one can consider how these
diagrams change as the modulation changes. Figure 8.10 shows the
laser structure and the optical energy/distance and stopband/
distance diagrams schematically. When the greatest current is driven
into the right-hand contact as in Figure 8.10a, the gain is greatest in
the right-hand section and light generated in the middle of this section
propagates to the left and right but, because of the Bragg grating, the
light evanesces on the left leaving very little light emerging from the
left-hand facet. More light emerges from the nearer right-hand facet.
The laser grating is designed so that the laser naturally prefers to lase
in the — 1 mode (the lower-frequency mode in the stopband diagram).
268 Distributed feedback semiconductor lasers
k ,d A/
AJUUUUUUU1 output lUUUULAJUUU^^ output
current drive into facet facet
current drives into split contacts
single contact
2.6 |-
relative
relative wavelength shift,
wavelength shift, Angstroms
Angstroms
Figure 8.11 Output intensity and time-resolved chirp from commercial single-contact
MQWDFB laser contrasted with unoptimised 3-contact push-pull bulk
DFB laser
Pattern is 010 at 2.5 Gbit/s. The total bias level for single contact is
60 mA compared with 90 mA for push-pull and the modulation depth
is 2 VH, (data after Flanigan [27])
IL Id.c. IR
luUUlAAJUUlAiULaAAAJUUUUUUUI
iT Ax
>*
10.0'
+20mA
8.0 0mA
(symmetrical)
-10mA
6.0
-20mA
-30mA
4.0 -40mA
(single ended)
in-phase
2.0
waveleijgth shift,
Angst oms
0
time, ns 1.0
4.8
+60 ps
4.0
+40 ps
32
wavelength +20 ps
shift,
Angstroms 2.4 Ops
(symmetric)
-20 ps
1.6
-40 ps
0.8
-60 ps
0.
time, ns 1
Figure 8:14 Modelled effect of temporal asymmetry on the time-resolved chirp for a
three-contact MQW device push-pull modulated at 1.5 Gbit/s with a
010 pattern.
The output was taken from the right facet, and a positive temporal
asymmetry means that the modulation switch on the left contact leads
the switch on the right contact. The DC bias levels on the end contacts
were 22.5 mA, the centre contact bias was 50 mA. The modulation
depth on the end contacts was 40 mA.The asymmetry is as defined in
Figure 8.12 (after Flanigan [27]).
400km
motion
effective
gain (cm" ) | Af~k7?h
100
bulk
-100 in
1/1650 1/1600 1/1550
frequency (in reciprocal nm)
conduction band
-i
/ WWWWV\AAA\A/WVW
s vwwwvw
i VWSAMAA/V
1 2 3
2 contact DFB 2 contacts: gain +DFB reflector 3 contacts: gain + phase +DFB
When contacts A and B are driven identically, the device behaves like
a conventional DFB laser. One way of making a tunable laser drives one
contact (say A) alone until the laser reaches threshold and then, on
increasing the current through contact B, tuning occurs. The
important feature is that the asymmetry along the laser's length in the
drive leads to an associated asymmetry in the effective linewidth-
broadening factor aH and also asymmetry in the gain saturation
enabling the gain and phase tuning to be partly decoupled to different
sections of the laser. In this device it is helpful to use quantum-well
active regions so that, on the side which is close to gain saturation, the
aH factor is increased, thereby increasing the change of frequency for
a given change in the carrier density. At 1500 nm, reductions in
wavelength of 2-3 nm are possible by increasing the current in
contact B before heating causes the refractive index, and thus the
wavelength, to increase again. The device is easy to make but the
tuning range is limited if relatively predictable.
Laser 2 may be regarded as a Fabry-Perot cavity (contact C) with
feedback determined by one single-facet mirror, and one reflecting
stack or Bragg grating (contact D) which is primary in determining the
FP mode that is selected. The net round-trip gain from the reflector
and the FP section must be appropriate but the important thing is that,
at approximately the lasing wavelength, increasing current into D directly
reduces the refractive index of the waveguide via the plasma effect,
and hence reduces the wavelength of the peak reflection from the
Bragg grating, thereby tuning the laser. To avoid the difficulties of
interaction between gain and frequency in the Bragg section, this
section can be fabricated in a wide-bandgap waveguide which neither
absorbs nor gives gain to the lasing wavelength. Any spontaneous
emission from this guide lies outside the frequencies of operational
significance. In general with only two contacts, sufficient change in
wavelength will lead to a jump to an adjacent Fabry-Perot mode
('mode hopping') so as to maintain the required round-trip phase
matching.
This drawback of mode hopping is avoided in laser 3, where the
Fabry-Perot section of laser 2 is now divided into two distinct regions.
The region with contact F is a section of plain waveguide made from
material with a bandgap which is too large to absorb or give gain at the
lasing wavelength but can have its refractive index changed through
the plasma effect and so provide the correct round-trip phase as the
lasing frequency changes without creating a mode hop. The region
with contact E provides the required round-trip gain (with some effect
Future devices, modelling and systems analysis 279
on the phase) as the Bragg reflector is tuned via the plasma effect as
in laser 2.
Variations can be created, for example, by having the region, with
contact F, formed from an intrinsic layer sandwiched within a P-N
junction (forming a P-I-A/junction). Then with a variable reverse bias
with strong enough electric fields (>10V/|mm) the conduction and
valence-band edges within the I-region shift so as to increase the
bandgap energy. This is known as the Franz-Keldysh effect [35] in
bulk materials, or the quantum-confined Stark effect [36,37] in
quantum-well materials, and is primarily employed in electroabsorp-
tion modulators for high-speed low-chirp applications. The change in
absorption also implies (from the Kramers-Kronig relationships) that
there are changes in the refractive index which may be employed for
tuning a laser. However, with those structures and voltages which are
compatible with laser operation, the total refractive-index change
available is typically less than 0.1% with a fractional tuning range of the
same order. Four-section devices exist, similar to laser 3 but with
tunable Bragg reflectors at both ends and phase- and gain-compensat-
ing regions.
L.
mm*
-1 AQ ^+1
tuning TV
anti-reflection coatings
. Contact A y Contact B
IIP v///)>/////////,
contact C y
modes coincide
T M11 in 11 I I I combA + B
^ I
mode I I I I l l lI I I I I I combA + C
amplitude
frequency
contact A
contact B
active I;
passive tuning
guide layer
lasing
filament
Contact C
In summary, while there are many candidates for the role of the
tunable semiconductor laser, the really broadband devices are all
difficult to characterise and require complex arrangements for the
control of several drive currents. These features combined with the
thermal limitations are a recurring problem of such tunable lasers and
become particularly acute if the laser degrades with time. The control
problem will need to be solved before such lasers become suitable for
mass applications.
o
optical output optical output
\ \ gain section
2nd order grating reflector & output coupler
a
optical output
optical output
optional anti-reflector / \
gain
0.2
0.15
99% reflectivity
8n/Sn
0.1
0.05
0 10 20 30 40 50 60 70 80
approximate number of layer pairs
lasers, and may well be even more serious with VCSELs unless these
lasers can operate with lower threshold-current densities. The strong
temperature dependence of the nonradiative Auger recombination
will also be of concern.
Making allowance for such changes, analysis and experiment does
confirm that there is an overall increase in the speed of response [87]
for the following reasons:
(i) VCSELS tend to use high-gain material, at high carrier densities
(~5 X 1018 cm" 3 ), which then reduces the stimulated and sponta-
neous carrier lifetime, though it restricts the maximum
temperature of operation;
(ii) the small dimensions of VCSELs give small electrical capacitances
which are easier to charge and discharge; and
(iii) the resonance frequency increases as V(/— I(h) (Appendix 4) and
the low threshold currents in VCSELs means that one can drive
the laser at, say, 10 times threshold whereas edge emitters are
rarely driven at more than three times threshold.
If one uses time-domain modelling for VCSELs, some changes from
the techniques used for edge-emitting DFB/DBR lasers are advised.
The significant refractive-index steps in the Bragg reflectors, com-
bined with their length, suggests that the appropriate space step for
the computation is that of a single layer or perhaps layer pair within
the stack compared with the tens of grating periods used previously.
Using such small step lengths gives excellent results for computing the
reflection spectra of the stack mirrors, but also requires correspond-
ingly smaller time steps ~10~15s, rather than ~10~13s for
conventional DFBs.
8.5.3 Construction of VCSELs
VCSELs are fabricated using virtually identical technologies to those
developed for edge emitters but one critical process (facet cleaving) is
replaced by a different critical process: reflector stack fabrication.
Figure 8.28 indicates some varieties. Testing is potentially much easier
than with edge emitters since it may be done without dicing chips from
their wafer; this is a very considerable advantage.
Formation of dielectric-mirror stacks by vacuum evaporation of
materials such as SiO2, Mg2F, TiO2 etc. is a well established and reliable
technology. In-situ monitoring of reflectivity during deposition is often
employed, so that drifts in deposition rates, optical densities etc. can
292 Distributed feedback semiconductor lasers
active
region
oxidised proton
annulu bombarded
8.6 Summary
Four areas have been covered in this final chapter. The first has been
systems modelling where advanced numerical models can predict the
performance of assemblages of devices all of which are modelled
separately but can be brought together. It is believed that, as device
modelling improves in accuracy, this area will grow extensively and it
will reduce the costs of trying out new systems and speed up their
Future devices, modelling and systems analysis 295
8.7 References
1 LOWERY, A.J.: 'Computer-aided photonics design', IEEE Spectrum, 1997,
34, pp. 26-31
2 LOWERY, A. J.: 'Semiconductor device and lightwave systems perform-
ance modelling'. Technical digest of Optical Fiber Communication, 97
(OFC'97), OSA Technical Digest Series, 1997, 6, paper ThG6, pp. 267-268
3 LOWERY, A.J.: 'Transmission line modelling of semiconductor lasers: the
transmission line laser model', Int. J. Numer. Model, 1989, 2, pp. 249-265
4 LOWERY, A.J.: 'Transmission-line laser modelling of semiconductor laser
amplified optical communications systems', TEE Proc. J, 1992, 139, pp.
180-188
5 LOWERY, A.J., GURNEY, P.C.R., WANG, X-H., NGUYEN, L.V.T., CHAN,
Y.C., and PREMARATNE, M.: 'Time-domain simulation of photonic
devices, circuits and systems', Proceedings of Physics and Simulation of
Optoelectronic Devices, at SPIE Photonics West, 1996, 2693, pp. 624-635
6 YOUNG, T.P.: 'CAD tools for optoelectronic subsystems', GEC J. Res.,
1994, 11, pp. 110-121
296 Distributed feedback semiconductor lasers
63 PARK, C.Y, KIM, D.B., YOON, T.H., KIM, J.S., OH, K.R., LEE, S.W., LEE,
S.M., AHN, J.H., KIM, H.M., and PYUN, K.E.: 'Fabrication of wavelength-
tunable InGaAsP/InP grating-assisted codirectional coupler filter with
very narrow bandwidth', Electron. Lett, 1997, 33, pp. 773-774
64 ASGHARI, M., ZHU, B., WHITE, I.H., SELTZER, C.P., NICE, C ,
HENNING, I.D., BURNESS, A.L., and THOMPSON, G.H.B.: 'Demon-
stration of an integrated multichannel grating cavity laser for WDM
applications', Electron. Lett., 1994, 30, pp. 1674-1675
65 MILLER, B.L, DREYER, K, BEHRINGER, R.E., KOREN, M.G., CHIEN,
M., RAYBON, G., and CAPIK, R.J.: 'Low chirp wavelength selectable 1 x
6 arrays suitable for WDM applications', 15th IEEE International Semi-
conductor Laser Conference, 1996, paper Wl.3, pp. 129-130
66 SATO, K, and MURAKAMI, M. 'Experimental investigation of thermal
crosstalk in a distributed feedback laser array', IEEE Photonics Technol
Lett, 1991, 3, pp. 501-503
67 AMMAN, M.C.: 'Recent progress on wavelength tunable laser diodes',
Inst. Phys. Conf Sen, 1994, 136, pp. 257-264
68 MOLLENAUER, L.F., and WHITE,J.C. (Eds.): 'Tunable lasers' (Springer
Verlag, 1987)
69 OKUDA, H., SODA, H., MORIKI, K., MOTEGI, Y, and IGA, K.:
'GalnAsP-InP surface emitting injection laser with buried hetero-
structures',Jm./. Appl. Phys., 1981, 20, pp. L563-L566
70 MOTEGI, Y, SODA, H., and IGA, K: 'Surface emitting GalnAsP-InP
injection laser with short cavity length', Electron. Lett, 1982, 18, pp.
461-463
71 IGA, K, ISHIKAWA, S., OHKOUCHI, S., and NISHIMURA, T: 'Room
temperature pulsed operation of GaAlAs-GaAs surface emitting injec-
tion laser', Appl. Phys. Lett, 1984, 45, pp. 348-350
72 OU, S.S., YANG, J.J., andJANSEN, M.: '635 nm GalnP/GaAlInP surface
emitting laser diodes', Appl Phys. Lett, 1993, 63, pp. 3262-3264
73 TAKAMORI, T, COLDREN, L.A., and MERZ, J.L.: 'Folded cavity
transverse junction stripe surface emitting laser', Appl. Phys. Lett, 1989,
55, pp. 1053-1055
74 SAITO, H., and KONDO, Y: ' 4 x 4 surface emitting 1.55>m InGaAsP/
InP laser arrays with microcoated reflectors fabricated by reactive ion
etching', fpn.J. Appl. Phys. Lett, 1991, 30, pp. L599-L601
75 GOODHUE, W.D., DONNELLY, J.P., WANG, C.A., LINCOLN, G.A.,
RAUSCHENBACH, K., BAILEY, R.J., and JOHNSON, G.D.: 'Monolithic
2-dimensional InGaAs/AlGaAs and AHnGaAs/AlGaAs diode laser arrays
with over 50% differential quantum efficiencies', Appl. Phys. Lett, 1991,
59, pp. 632-634
76 STEGMULLER, B., WESTERMEIER, H., THULKE, W., FRANZ, G., and
SACHER, D.: 'Surface emitting InGaAsP/InP distributed feedback laser
diode with monolithic integrated microlens', IEEE Photonics Technol. Lett,
1991,3, pp. 776-778
77 EVANS, G.A., CARLSON, N.W., BOUR, D.P., LURIE, M., DEFREEZ,
R.K., and BOSSERT, D.J.: '14-W peak power grating surface-emitting
laser array', Electron. Lett, 1990, 26, pp. 1380-1381
Future devices, modelling and systems analysis 301
78 YU, S.F., PLUMB, R.G.S., and CARROLL, J.E.: 'Spatial active optical
switching by using grating coupled surface emitting DFB lasers', Electron.
Lett., 1993, 29, pp. 1147-1148
79 JEWELL, J.L., HUANG, K.E, TAI, K., LEE, Y.H., FISCHER, R.J.,
MCCALL, S.L., and CHO, A.Y.: 'Vertical cavity single quantum well laser',
Appl. Phys. Lett., 1989, 55, pp. 424-426
80 JEWELL, J.L., HARBISON, J.P., SCHERER, A., LEE, Y.H., and FLOREZ,
L.T.: 'Vertical-cavity surface-emitting lasers—design, growth, fabrication,
characterization', IEEEJ. Quantum Electron., 1991, 27, pp. 1332-1346
81 IGA, K., ISHIKAWA, S., OHKOUCHI, S., and NISHIMURA, T.: 'Room
temperature pulsed operation of GaALAs-GaAs surface emitting injec-
tion laser', Appl. Phys. Lett., 1984, 45, pp. 348-350
82 SALES, T.E., 'Vertical cavity surface emitting lasers', (Wiley, New York,
1995)
83 SHTENGEL, G., TEMKIN, H., BRUSENBACH, P., UCHIDA, T., KIM, M.,
PARSONS, C , QUINN, W.E., and SWIRHUN, S.E.: 'High-speed vertical-
cavity surface-emitting laser', IEEE Photonics Technol. Lett., 1993, 5, pp.
1359-1362
84 BABA, T., HAMANO, T., KOYAMA, E, and IGA, K: 'Spontaneous
emission factor of a microcavity DBR surface emitting laser', IEEE J.
Quantum Electron., 1991, 27, pp. 1347-1358
85 BABA, T., HAMANO, T, KOYAMA, E, and IGA, K.: 'Spontaneous
emission factor of a microcavity DBR surface emitting laser, 2, effects of
electron quantum confinements', IEEEJ. Quantum Electron., 1992, 28, pp.
1310-1319
86 HAMANO, T., IGA, K., and BABA, T.: 'Spontaneous emission behavior
and its injection level dependence in 3-dimensional microcavity surface-
emitting lasers', Electron. Commun. Jpn. 2, Electron., 1996, 79, pp. 46-54
87 SANDUSKY, J.V., and BRUECK, S.R.J.: 'Observation of spontaneous
emission microcavity effects in an external-cavity surface-emitting laser
structure', Appl Phys. Lett, 1996, 69, pp. 3993-3995
88 WATANABE, I., KOYAMA, E, and IGA, K.: 'GalnAsP-InP CBH surface
emitting laser with a dielectric multilayer reflector', Jpn. J. Appl Phys. 1,
Regul Pap. Short Notes, 1987, 26, pp. 1598-1599
89 CHAILERTVANITKUL, A., IGA, K., and MORIKI, K.: 'GalnAsP-InP
surface emitting laser (A=1.4[xm, 77 K) with heteromultilayer Bragg
reflector', Electron. Lett., 1985, 21, pp. 303-304
90 HOUNG, Y.M., TAN, M.R.T., LIANG, B.W., WANG, S.Y, YANG, L., and
MARS, D.E.: 'InGaAs(0.98 jxm)/GaAs vertical cavity surface emitting
laser grown by gas-source molecular beam epitaxy', / Crystal Growth,
1994, 136, pp. 216-220
91 SUGIMOTO, M., KOSAKA, H., KURIHARA, K., OGURA, I., NUMAI, T.,
and KASAHARA, K.: 'Very low threshold current density in vertical cavity
surface emitting laser diodes with periodically doped distributed Bragg
reflectors', Electron. Lett., 1992, 28, pp. 385-387
92 ZEEB, E., MOLLER, B., REINER, G., RIES, M., HACKBARTH, T., and
EBELING, K. J.: 'Planar proton implanted VCSELs and fiber-coupled 2-D
VCSEL arrays', IEEE J. Select. Topics Quantum Electron., 1995, 1 , pp.
302 Distributed feedback semiconductor lasers
106 SCOTT, J.W., GEELS, R.S., CORZINE, S.W., and COLDREN, L.A.:
'Modelling temperature effects and spatial hole burning to optimise
vertical cavity surface emitting laser performance', IEEE J. Quantum
Electron., 1993, 29, pp. 1295-1308
107 WILSON, G.C., KUCHTA, D.M., WALKER, J.D., and SMITH, J.S.: 'Spatial
hole-burning and self focusing in vertical cavity surface emitting laser
diodes', Appl Phys. Lett, 1994, 64, pp. 542-544
108 WIPIEJEWSKI, T, PETERS, M.G., THIBEAULT, B.J., YOUNG, D.B., and
COLDREN, L.A.: 'Size dependent output power saturation of vertical
cavity surface emitting laser diodes', IEEE Photonics Technol. Lett., 1996, 8,
pp. 10-12
109 FIEDLER, U., REINER, G., SCHNITZER, P., and EBELING, K.J.: 'Top
surface-emitting vertical cavity laser diodes for 10 Gbit/s data transmis-
sion', IEEE Photonics Technol. Lett, 1996, 8, pp. 746-748
110 NUMAI, T, KURIHARA, K., KUHN, K., KOSAKA, H., OGURA, I.,
KAJITA, M., SAITO, H., and KASAHARA, K.: 'Control of light output
polarization for surface emitting laser type device by strained active layer
grown on misoriented substrate', IEEE]. Quantum Electron., 1995, 31, pp.
636-642
111 TAKAHASHI, M., VACCARO, P., FUJITA, K., WATANABE, T, MUKAI-
HARA, T., KOYAMA, E, and IGA, K.: 'InGaAs-GaAs vertical cavity surface
emitting laser grown on GaAs (311) A substrate having low threshold and
stable polarization', IEEE Photonics Technol. Lett., 1996, 8, pp. 737-739
Appendix 1
Maxwell, plane waves and reflections
with c= (l//i o eo)1/2 = 3x 108 m/s, the velocity of light in vacuum. The
phase velocity of light within material is cw=3x 108/ nr m/s which varies
slightly with the optical frequency. In general, with optical gain or loss
where n{ is nonzero, (3 is complex.
The //fields are given from eqn. Al.l by putting d/dx=d/dy=0 with
d/dt-^jo) so that only a component Hy is found. For forward-
propagating fields:
(curl H)x= - dHy/dz=jPHy=ere0dEx/dt=ja>ereoEx (A1.6)
Hy=(8rs0/firfjL0)]/2Ex-^Ex/Hy=(377/n) O (A1.7)
In isotropic material, the E field could lie along any direction but the
orientation offieldswith respect to the vector direction of propagation
/?= ((3X, (3y, f3z) must be maintained. For plane linearly polarised waves,
the forward E field lies along one direction with Ep i^and /^forming
a right-handed set for the forward wave as shown in Figure Al.l and
similarly for the reverse wave Er9 Hr and /3r form a right-handed set.
I R
showing the choice of axes in Figure A1.2. Because all tangential fields
have to be continuous, the periodicity of these three waves, resolved
along the interface must be identical (Figure A1.3) to one another.
This periodicity is determined from /3zi, the axial component of the
propagation vector(0, fiyi, /3zi) of the incident wave, and similarly
replacing the subscript i with r and t for the reflected and transmitted
waves, respectively. Note that (3xi/r/l=0 shows that the light rays lie in the
yz plane. Hence, from the equal periodicity condition,
/3zl=(3zr=f3zt
The ray angles of the incident transmitted and reflected directions are
Appendix 1 307
given from:
tan 6l=pu/pyi tan 0=Pzr/Pyr tan 6t=Pzt/Py( (A1.9)
From eqn. A1.5, one has for the i and rwaves in region 1 and t wave in
region 2, respectively,
Hence from eqn. A1.8 one recovers Snell's law for reflection and
refraction in lossless material with real refractive indices
sin di=sin ^ ( r ^ / n j s i n 0t (Al.ll)
The 'critical ray' has ^=90° when sin fli=(n2/n1)<l. When sin dt>\
it is not possible to find any real solution to Snell's law of eqn. A l . l l
but explicit expressions can still be found for the vertical 'propagation'
in the region 2:
i i \
T 2
Figure A1.5 Reflection/transmission amplitudes at surfaces
F=£forTE fields
F=J/forTM fields
giving
S77(co/c)H=-/3yEx (A1.15)
The continuity of the tangential-field components Ex and Hz requires
Exi+Exr=Exl ^yiE^liyrEx=^JExt (A1.16)
Noting that fiyi/r=± (wn 1 /c)cos 0,-;fiyt=- (co n<>/c)cos 0t and 0,-= 0r
nx cos O^— nYcos 0rExr=n2 cos 0tExt (Al.l7)
Eliminating appropriately from eqns. A1.15 and Al.l7 leads, after a
little work, to reflection and transmission coefficients:
cos
p=Exr/Exi={(nl cos 0,-^2 &t)/ (ni c o s 0,-+w2 cos 0,)} (Al.l8)
T=Ext/Exi={2nx cos 0J(nx cos 0,-+n2 cos 0,)} (A1.19)
Appendix 1 309
(A1.25)
Noting again that fiyi/r=±{a)nl/c) cos 6-, /3yi/r= - (con^/c) cos 0t and still
(A1.30)
giving
(iv) Brewster's angle: A new special case arises in that it is now possible
to have pTM-Q at a critical angle: the 'Brewster' angle. This happens
when
0,)}=O (A1.32)
Using SnelFs law and writing n^ (1 — sin2 #,) = ?%2 (1 — sin2 6t)
gives sin ^ =/22/V(n2+72|) or tan d^n^/rix (A1.33)
This is the Brewster angle of incidence where there is no reflection of
the TM wave and only a TE wave is reflected. This is the principle on
which Polaroid glasses help to cut out some of the glare from
reflection. This angle is of special significance in gas lasers, and has
also been used in semiconductor lasers to cut out facet reflections (at
one polarisation) at the semiconductor/air interface. However, the
steepness of the Brewster angle (#i~70°) for a typical semiconductor/
air interface makes this unattractive compared with the present
technology for low-reflectivity surfaces where multiple dielectric thin
films are deposited with precision over the surface and can keep
reflections below 0.1% for near normal incidence.
A1.9 References
1 STRATTON, J.A.: 'Electromagnetic theory' (McGraw Hill, New York,
1941)
2 RAMO, S., WHINNERY, J.R. and VAN DUZER, T.: 'Fields and waves in
communication electronics', (J.Wiley, New York, 1993), 3rd Ed.
3 LONGAIR, M.S.: 'Theoretical concepts in physics' (Cambridge University
Press, 1984) (case study 2 gives an excellent historical perspective)
4 HERZINGER, CM., LU, C.C., DETEMPLE, T.A., and CHEW, W.C.: T h e
semiconductor waveguide facet problem', IEEEJ. Quantum Electron., 1993,
29, pp. 2273-81
5 YU, B.M., and LIU, J.M.: 'Gain margin analysis of distributed feedback
lasers for both transverse electric and magnetic-modes', IEEEJ. Quantum
Electron., 1992, 28, pp. 822-832
6 SALE, T.E.: 'Vertical cavity surface emitting lasers' (Research Studies Press,
Taunton; Wiley, New York, 1995)
Appendix 2
Algorithms for the multilayer
slab guide
exp(+ja>0 and the magnetic fields relate to the electric fields from
jmodHz/dt=dEx/dy that with the wave equation give
7\
iclA 7,-77.
(A2.3)
(A2.4)
Then at the 2/1 interface one has the fields at the start of region 2
given from
Then at the 3/2 interface one has the fields at the end of region 2
given from
71727374 v
I tx^ J 4 3
316 Distributed feedback semiconductor lasers
It is again important to get the sign of y4 right because, just as for the
top layer T, the bottom layer (here layer 4) has to have an evanescent
wave with only Ex4_ evanescing away from the guiding layer. One also
notices that £M+ is zero and consequently we have the situation that
7i727s74
Z=EJHX= -jyr/ere0co= -j
2 o =( / i 0 /e 0 ) 1 / 2 (A2.12)
As before, one starts with the top layer Twith a single evanescent wave
varying as exp ( — yry) with a wave impedance of — j7^{yje,.#b)- The
layer 1, of thickness du has a 'forward' wave exp( — yxy) with a wave
impedance of —jl^^yjefc) and a 'reverse' wave varying as exp(yxy)
with a wave impedance of jZ^{yr/erko). As in the TE case, r+ and r—
indicate, respectively, the 'forward' and 'reverse' types of evanescent
wave in the layer r. Then at the interface we have matching H fields
and E fields, respectively:
(Hxl++HxX _ ) n = (HxT++Hxl^_) n
-H^),, (A2.13)
1= (ST7I) ~
l - Jn
Jn \{sryi)~ (exyT) {sTyx) + (exyT)
(A2.14)
Moving a distance — dx (i.e. back from the interface) to reach the
interface 2 / 1 , one accounts for the propagation inside the layer using
Then at the 2/1 interface one has the fields at the start of region 2
given from
Then at the 3/2 interface one has the fields at the end of region 2
given from
(A2.17)
12
As for the TE mode, this process is repeated until the bottom layer is
reached. For example, with the five layers as illustrated, define
(A2.18)
M= D^sD^H^H^ r (A2.19)
Jx4+ l
=M\
nx4x4-_ 143
(A2.20)
leading to the requirement for a solution, similar to that for the TE
mode, that
l/9l = IM n /M 12 l=0 (A2.21)
The solution technique is followed as for the TE mode.
The word of warning given elsewhere is repeated: it is the E fields
which interact with the electronic dipoles in the material caused by
electrons and holes recombining to give gain. Consequently it is the E
fields which are important in determining the strength of the
interaction and in determining the confinement factor in a laser.
When discussing the TM modes, any computing programs constructed
may well find it easiest initially to compute the H field but the final
interaction should compute the total E fields to indicate the strength
of the interaction with the electronic gain along the guide and to
estimate the confinement factor using the E fields.
Appendix 2 319
E(y)
forward radiation
+cos9]
Vz [1 + COS6 ]
i
/2 p^ [1 - cos e ]
refractive index n
-cose]
refractive Index n
Figure A 2.5 Radiation with H field normal to plane of emission
The radiating aperture (white) with its known fields are placed just in
front of the laser facet (shaded) to show more clearly that the reverse
reradiated light is reflected from the facet
Figure A2.6 Replacing etched/refilled thick layer with many thin layers for calculation
of a slab wavguide which includes a Bragg grating
Permittivity is then averaged along the length of the guide to find the
field patterns
using the form of eqn. 4.37. However, the reader is advised to become
familiar with running the slab-waveguide program without gratings
before requesting the calculation for a Bragg grating.
A2.5 References
1 SNYDER, A.W., and LOVE, J.D.: 'Optical waveguide theory' (Chapman and
Hall, London, 1983)
2 ADAMS, M.J.: 'Introduction to optical guided waves' (Wiley, Chichester,
1981)
3 ROZZI, T.E., and IN'TVELD, G.H.: 'Fields and network analysis of
interacting step discontinuities in planar dielectric waveguides', IEEE
Trans., 1979, MTT-27, pp. 303-309
4 RAMO, S., WHINNERY, J.R., and VAN DUZER, X: 'Fields and waves in
communication electronics' (Wiley, 1994), 3rd ed., pp. 618-627
5 HOCKHAM, G.A.: 'Radiation from a solid state laser', Electron. Lett, 1973,
9, pp. 389-391
6 LEWIN, L.: 'Obliquity factor for radiation from solid state laser', Electron.
Lett., 1974, 10, pp. 134-135
7 LAU, T., and BALLANTYNE, J.M.: 'Two-dimensional analysis of a dielectric
waveguide mirror',/ Lightwave Technol, 1997, 15, pp. 551-558
Appendix 3
Group refractive index of laser
waveguides
The group velocity is the velocity of a wave packet (i.e. the velocity of
energy) that is centred on a central carrier frequency f = (0/ 2TT. The
group velocity is different from the phase velocity for two reasons: (i)
the waveguide changes the propagation coefficient as a function of
frequency; and (ii) the material permittivity changes with frequency so
that the propagation coefficient in the material changes with
frequency. This appendix illustrates these two roles using a symmetric
three-layer waveguide so that one can appreciate the physics through
putting numbers into an analytic solution.
The analytic solution for the propagation coefficient as a function of
frequency is well documented for the TE mode in a symmetrical three-
slab guide [1-4] . The guide is illustrated in Figure 3.5, where the
effect of weak transverse guiding is neglected and only the strong
lateral guiding is considered. The thickness of the central layer is taken
here to be d and it has a relative permittivity erl and is surrounded by
layers with relative permittivity er2. The axial-propagation coefficient /3Z
determines, through the wave equation with an angular frequency co,
that the central region has a lateral propagation coefficient f3':
where /^=V (co2 e0 /x0) = (co/c) =2TT/ A with A the free-space wavelength.
In the outer regions there is an evanescent wave decaying spatially in
y with a rate
y)=Pl-kler2 (A3.2)
There are several different notations used in the literature giving
equivalent results, but the use of symmetry of the guide gives a
Appendix 3 325
straightforward result [3,4] that the field profile with continuity of the
E fields must be of the form:
EX=EO cos((3yy) [0<y<d/2] (A3.3)
(A3.5)
JeA.JjL.*^ (A3.6)
where Ser= (srl — s^). A normalised thickness for such a slab guide may
be given in terms of the 'Vparameter' where V= (co d/c) V(Ser) and
then, defining the optical phase change across half the width of the
guide ® = (fiy d/2), the equation which determines the guide's main
mode of operation is:
&
cos 0 = 2 - (A3.7)
(odV(Ser) 0
\^IL 2 ~ (A3.8)
c =V= cos 0
326 Distributed feedback semicondu ctor lasers
c/n
0.29
c/n V
0 1 2
normalised frequency, V
Figure A3.1 Phase velocity against normalised guide frequency
The outera nlayers of the slab guide have the lower refractive index
n.2 = V £,2 d with thin enough guides the energy resides in the outer
layers giving the higher phase velocity. As the optical thickness
increases, V increases and the energy is stored in the central part of the
guide. At high enough V, the phase velocity is c/n} where nx= Ve r] but
at such high Vnumbers the guide is overmoded
cod
=—=c/V{erl - &r(cos (A3.10)
J1
(3d
From the equations, observe that when cos 0 - 0 then vp~c/v(srl).
Physically, this arises because the optical width of the guide is so large
that all the stored optical energy is, effectively, in the central region.
However, when cos @ ~ 1 where the optical width is small then all the
energy is effectively in the cladding and so vp~c/v(er2)- The
fundamental mode only extends usefully up to lvalues of around IT
when the first higher-order mode starts to propagate so that @ is only
considered up to values around TT/4 rather than the full range. Figure
A3.1 plots the phase velocity as a function of the normalised frequency
V for materials given in Table A3.1; these are typical materials for
heterostructure DFB lasers. Typical V numbers for d=0.1 to 0.2 jxm-
thick guiding layers at wavelengths around 1.55 fxm are in the range
0.6-1.3.
Appendix 3 327
(A3.12)
+ 0 tan©
approxfrhate
4.4 useful range - " " 'group index:
guide + material
4
index
value
3.6 • group index: guide only
- V-V-----
3.2 phase index
0 1 2 3 4
normalised frequency, V
Figure A3.2 Effective group and phase indices against normalised frequency
Group and phase indices as appropriate for a guide formed from a
cladding layer InP with a refractive index of 3.168 and a core layer of
In^Ga^ASyP^, with a refractive index 3.531. For a core-layer thickness of
0.2 |xm, one finds V~ 1.25 for A ~ 1.55 |xm
A3.1 References
1 HAUSS, H.A.: 'Waves and fields in optoelectronics' (Prentice Hall,
Englewood Cliffs, 1984), p. 175
2 KOGELNIK, H., and RAMASWAMY, V.: 'Scaling rules for thin-fiber optical
waveguides', Appl. Opt, 1974, 8, p. 1857
3 COLDREN, L.A., and CORZEINE, S.W.: 'Diode lasers and photonic
integrated circuits' (Wiley, 1995)
4 ANDERSON, W.W.: 'Mode confinement and gain injunction lasers', IEEE
J. Quantum Electron., 1965, 6, pp. 228-236
Appendix 4
Small-signal analysis
of single-mode laser
The effective photon lifetime r'p has been used as in Chapter 4. The
intervalence-band-absorption term bNS is added here for complete-
ness and the spontaneous emission is recognised as being formed from
bipolar recombination so that, when the holes and electrons are of
equal density, the radiative spontaneous emission is proportional to
N2. The gain is linearised about the transparency density Ntr but gain
saturation with optical density is included. The spontaneous-emission
coupling /3sp is now assumed to have absorbed the confinement factor.
The electronic eqn. 4.2 is also written more fully as
330 Distributed feedback semiconductor lasers
(A4.5)
With So zero, the electron density increases with the current input
until N=N(h, when
^AN^BNl+CNl (A4.8)
qd
As / continues to increase, the optical output then increases almost
linearly with drive current above threshold, provided that the gain
saturation is negligible:
qdGm(Nth-Ntr)
The dynamic first-order equations are found by retaining only the first-
order terms in eqns. A4.1 and A4.2 using eqns. A4.4. In particular, the
small-signal version of eqn. A4.1 may be rewritten after some
manipulation as
EjiV, (A4.10)
V W
where
where
1 9
-=A+2BN0+3CN20 (A4.12)
fH2Si==~ (A4.13)
T
r/2/
where
and B2=Gm(N0-Nlr)
r
1
ram(N0~Ntr) rl-=r6iv0+- 7 (A4.14)
7
0 Tp
Rearranging to obtain the change in photon intensity Sx as driven by
the changes in the current density J{,
(-co2M+ja>My+a>l)S^^ (A4.15)
where
T
rf2T./l \T^2 T
d\/
Neglecting the IVBA term b gives the damping and resonance from
Appendix 4 333
Eqn. A4.15 shows a classic damped resonant response with the overall
damping rate determined principally by the effective recombination
time but with increased damping at high output-power levels and with
gain saturation. The resonance frequency increases as JSQ [or
alternatively as ^(Io — Ilh) when the steady current increases above the
threshold value Ith, assuming that there is negligible gain saturation].
Notice that G'm in eqn. A4.16 is the differential power gain per unit
time so that knowledge of the resonance frequency as a function of
intensity 50 can yield an estimate of this differential gain. The most
important consequence of gain saturation is that it can both lead to a
marked reduction in the resonant frequency cor at high power levels
and also give additional damping of the optical response to changes in
current. Therefore at high enough optical-power levels, gain satura-
tion slows down the dynamic response of a laser to direct modulation
of the current drive.
At this point a little care is needed to note that S} already represents
an optical power so that the — 3 dB bandwidth of the optical-power
modulation for the laser occurs when
\~a>M2+jcoMy+cD2r\=2aj2r (A4.17)
2
Usually a>My<(o ., allowing the approximate — 3 dB bandwidth to be
defined as
*>M-3<*ir(V3)<», (A4.18)
Measurements of the small-signal AM response are often normalised to
the response at low frequency. The magnification at resonance is given
by
Magnification \ramam = - r = ^ Tp
' <- (A4.19)
y
)
which has a maximum when the photon density is given by
1
(A4.20)
334 Distributed feedback semiconductor lasers
are
(A4.21)
(Xt mAs.ni * e*r * cap
UN0<po-N,)S\-- '
qdqw(l+jo)MTmp)
(A4.24)
^-^^ANqu-BNlv~- l {
CNsqw- Gm(NHqw-Ntl)S (A4.25)
dt r 'hap
r r, *
which limits the amplitude response rather than the capture rate of
charge carriers into the quantum wells. On the other hand, if it takes
too long for the wells to capture the charge carriers, there is a
significant first-order time constant which can limit the response well
before the photon-electron resonance plays a role and this pheno-
menon exhibits itself principally by an additional roll-off in the small-
signal AM response.
(A4.29)
Appendix 4 337
dn dn
r (A4.30)
diV wg
'""dN
The detailed results are straightforward if 'messy' and shed little light
on the physics other than to once again indicate that the key effect is
to add an additional delay so that the form of the FM response has the
general form
N(/
Wl (A4.31)
(J\ / qdqj (1 +J(OMTMP) ( - (o2Ma+ja)My+ o)2r)
where a, y and cor2 are similar to the parameters that previously would
have been found for the AM response and £ is again proportional to
j(oM as a)M increases. The important point is that transport adds an
additional roll-off term for the FM response just as it does for the AM
response.
Figure A4.1 summarises, in its upper sketch, the schematic features
of AM responses as found from the classic electron-photon rate
equations using logarithmic scales. The damping factor has been
chosen arbitrarily to give a peak magnification of 3 dB. It may be seen
that the transport time constant needs to be such that a>r rcnp< 1 for
transport effects to have negligible effects on the frequency response.
The lower sketch in Figure A4.1 shows the FM response. If there were
no nonlinear gain, the FM response would tend to zero as the
frequency fell. Nonlinear gain flattens the response. Note that the
338 Distributed feedback semiconductor lasers
10 r
relative output dB
0
3 dB
-10
-20
AM Response
-30
-40
-2 -1 coxyr/co 0
M r
1
relative output dB
10
-30
-1 co^/co
(A4.32)
]
dt T\dt Tp/ rr qW
^=A^+F7^7 (A4.34)
and this, along with the effective photon lifetime, is used to normalise
parameters:
a^ «3 xy
r.4 eS^ (A 4.35)
'p th 'p >C
d^=(A^)^__
A J aT H
dTL \ —a -
dLn
7+Ln )-NH(l- frp) +D (A4.38)
^dTn \dT,
—a '
1
—=-[-Ln-(M-(Ln/0}(l-^)+Z)] (A4.41)
U 1n £y
Note how there are then four normalised parameters ($sp, a, Q and es
which determine the normalised light/time characteristics for a given
normalised drive D. Note also that, in general, the light changes at a
much faster rate than the electrons because of the factor l/Qin the
normalised rate equation for the electrons. The electrons vary on a
time scale governed mainly by the electron-recombination time, which
is longer by a factor of Q (typically ~ 1000) than the photon lifetime
which governs the rate of change of the photons. This is an important
point in modelling because it is the photon lifetime which determines
the temporal step that is required.
There are two MATLAB programs, fpstat and fpdyna, in directory
fabpero. These allow exploration of the qualitative effects of sponta-
neous-emission parameters and gain saturation.
In summary, the rate-equation approach has been extremely
successful in dealing with several major characteristics of Fabry—Perot
lasers where the fields are sufficiently uniform and phase effects of
fields are not important. However, considerable care is required in
applying the results from such equations to distributed-feedback lasers
where phasing, nonuniform fields, and the distribution of the
feedback change the laser's behaviour in important ways.
A4.6 Reference
1 VANKWIKELBERGE, P., BUYTAERT, E, FRANCHOIS, A., BAETS R.,
KUINDERSMA, PL, and FREDRIKSZ, C.W.: 'Analysis of the carrier
induced FM response of DFB lasers', IEEE J. Quantum Electron., 1989, 25,
pp. 2239-2254
Appendix 5
Electromagnetic energy exchange
p =
s~\ equilibrium
\Q } charge
gives
(A5.4)
The Poynting vector 2?= (Ex H) gives the instantaneous power per unit
area and, to give the power from the laser, must be integrated over the
xy plane of the appropriate laser facet of effective area si:
To see how these results can be reconciled with the particle model of
Section 4.1, one must recognise that the electrical polarisation has
three distinct components:
(i) A real susceptibility ^ r with P=xrs0E
In this first component, the dipole current {dP(t)/dt} is in quadrature
with the electric field and from eqn. A5.4 there is no mean real power
exchange. There is reactive power where the power flows back and
forth between the electric and magnetic fields, and so the effects of Xr
are reflected in the stored electric energy within the material which at
any instant is given from
U=UeoenE.E)={k^XrE.E)+UeoE.E)} (A5.6)
using the fact that srr= l+xr-
There is an important theorem in electromagnetism which states
that, at resonance or in the normal mode of a guide, the mean electric-
energy density is identical to the mean magnetic-energy density [1, 2].
The laser is both resonant and uses a guide for the optical fields so
that, at the lasing frequencies, there is an equivalence of the mean
values of magnetic and electric stored energies where, using the
phasor notation:
U=}(sosn.E.E*) = Um=j([jL0H.H*) (A5.7)
Here the fields refer to the peak values and the additional half has
appeared because this is the mean of the instantaneous energy which
peaks every half cycle.
(ii) An imaginary susceptibility Xt w ith P=jXi ®oE
This second component gives [dP(t) / dt\ directly in phase or in
antiphase with the fields and results in a direct exchange of power.
Here sign conventions are important. If one assumes, as here, a time
344 Distributed feedback semiconductor lasers
(A5.ll)
Using the vector identity for the divergence of the Poynting vector,
which gives the electromagnetic-power flux,
On substituting for 3H/dt and dE/dt from eqns. A5.9 and A5.10, the
instantaneous rate of energy-density exchange is given from
(ExH)zdxdy
f f ff
I I Ufdxdy= I I sn.SQEj.Ejdxdy=h.(x)QSjSi (A5.15)
(EfxHf)zdxdy=fi(o0Sfs&vg (A5.16)
• / /
where vg is the group velocity. Here the important new point has been
made which differs from the rate equations of Section 4.1, namely that
power flow and energy stored can both be associated with the forward
(and also the reverse) directions along the guide. Eqn. A5.14 may be
recast to give the mean forward-power density and associated energy-
density exchange as in eqn. 4.23 from
1 8Sr dSr
~~+—- - 2g(iV) Sf > = (fico0si) sspf (A5.17)
The factor of 2 for 2 g(N) is again here because g(N) is used as the
field gain per unit distance as function of electron density //and 2g(N)
refers to the power gain per unit distance. Here this net optical power
Appendix 5 347
co0sris0Ef.Efdxdy (A5.18)
/ /
The spontaneous emission combined with its coupling factor /3sp into
the forward-propagating field is determined from the mean temporal
value of
=V
Equations A5.17 and A5.20 combine to give S= (Sj+Sr). The averaging
over one cycle, as in eqn. A5.15, is kept in future equations but not
shown explicitly. Integrating over the laser's length with Y = s&L gives:
U dx dy dz=
The power flowing out at the 'right' and 'left' facet ends gives the rate
of escaping photons from
348 Distributed feedback semiconductor lasers
(ExH)zdxdy\
The mean rate of net stimulated emission over the whole volume is
given from
o0sris0E.Edxdydz (A5.23)
Note that, because the calculation is concerned with the fields, it has
by default included the confinement factor and any loss which was put
in separately in Section 4.1. The net output of the radiative
spontaneous emission (N/TSP) to the lasing mode taking into account
both directions of propagation of energy is given from
(A5.24)
^ j > ) (A5.25)
dz
as the predominant components. The magnetic-field components may
be eliminated from Maxwell's equations to obtain the wave equation
for Ex, retaining a frequency-dependent real permittivity for the
material but ignoring the frequency dependence of the much smaller
imaginary permittivity:
djsponlx
y) }8
* °
(A5.26)
A single-mode solution to this equation is solved for steady-state
amplitudes with negligible loss and gain at one particular frequency o>
and a corresponding value of /3. These values of a) and yS are the same
over the whole cross-section of the guide, so that
Ex=Ef(z, t) u(x,y) where Ef{z, t)=Eoexp (](a)t-/3z)}
For this steady-state solution Eo cancels with
After solving this equation for the real Vertical' variation u(x, y) with
an index-guided wave (as, for example, in Chapter 3) relate co and /?
giving an effective permittivity srreg-:
(A5.28)
Now suppose that the mode no longer has a steady-state amplitude but
that
Ef (z, t) =F(z, t) exp{j(cobt- ftz)} (A5.29)
where F(z, t) varies slowly in space and time with respect to
exp{j(cob t— f3hz)}. In using phasor notation, one can 'equate' co with
(—jd/dt) but now one needs to separate the low frequencies, varying
with rates noted by (d/dt), from the optical frequency ca/r This is done
by writing <0—>a)h—j(d/dt) and a)2—>cob2 — 2 cobj {d/dt) to a first order
in the slow changes denoted by {d/dt). It follows from Taylor's
theorem that the frequency-dependent permittivity becomes a tempo-
rally changing permittivity
350 Distributed feedback semiconductor lasers
The result is more useful if one uses the effective modal area:
y)
u(x,yydxdy=gsl/pb (A5.36)
J J 2e^
dR 1 (dR\
(A5.39)
vg
A5.4 References
1 COLON, R.E.: 'Field theory of guided waves' (McGraw Hill, 1960)
2 SLATER, J.C.: 'Microwave electronics' (Van Nostrand, 1950)
Appendix 6
Pauli equations
This appendix gives the detailed calculations for finding the steady-
state fields in a uniform DFB laser with uniform gain, and thereby
finding the threshold conditions. The results are essential if one wishes
to make comparisons with the numerical algorithms to estimate the
accuracy of these algorithms. The appendix also provides a tutorial on
the use of Pauli matrices for coupled differential equations which
enables one to extend the concept of exp(yz) to exp(Mz) where Mis
a matrix: a useful concept for linear differential equations.
Take the coupled-wave equations
1 r) /? r) /?
*ft(>'"t ( A 6
-2)
Pauli matrices are widely used in the quantum theory of spin and
polarisation of light and the way in which the fields are coupled, but
they can also be used to provide neat solutions for other physical
problems (see, for example, [1]) and in particular solutions for
coupled linear equations. Define then the vector F and the Pauli
matrices
cr9 =
+a
^Jt "Tz~JK(riF~ ^S-jS)F=O (A6.4)
If, for example, there were no gain, phase or space variations, then
- — ~JKaxF=0 (A6.5)
vgdt
The reader who has not met the concept of matrix exponentials will
rapidly appreciate that one may, using a Taylor series, formally write
exp(jatf)=l+/irtf- (1/2!) ( ^ ) 2 - 7 ' ( l / 3 ! ) (criy)*+> • • (A6.6)
Using crf=l (more strictly this is the identity matrix), one recovers the
identity
expiry) =cos(y) +jax sin(y) (A6.7)
A solution then for eqn. A6.5 using Pauli matrices yields
F( t) = cos (Kvgt) +jcrx sin (Kvgt) (A6.8)
The object now is to integrate in space with a constant-frequency input
with the offset frequency determined by 5 so that we rearrange eqn.
A6.4 by multiplying through by cr3 to obtain
dF 1 dF
—+or 3 -a3(g~j8)F+Kcr2F=0 (A6.9)
dz Vg dt
Neglecting the time variation for the moment, the formal solution of
this equation is
F(z) =exp{ -j(o-3D- a2JK)z}F(0) (A6.10)
where, for a shorthand in this appendix only, D=8+j g.
However, because cr3 a2=~~ a2 a^ it may be shown by using Taylor
series as in eqn. A6.6 that exp{ — j(a3 D— cr2 JK)Z} cannot be the same
as expK-703 D)z] exp{(- CT2K)Z). The problem, stated briefly, is that
the matrices <r3 and a2 do not commute and so the matrices formed
from exp{( —703 D) z] and exp{( — <T2K) Z) cannot commute. One has to
consider
(o"$D~ CF2JK) (cr3D— CT2JK) = (D2 — K2) — (asa2+a2as)jDK= (D2 — K2)
(A6.ll)
showing that
- a2JK)2n=(D2-K2) n
(A6.12)
354 Distributed feedback semiconductor lasers
Apply eqn. A6.14 to find the solutions for eqn. A6.9 for a length z:
S 3
F(z) =cos(&z)F(0) + ^ ^ {(g-j8)crs-K(T2}F(0) (A6.15)
F(z) a+ d jb F(0)
(A6.16)
R(z) —jb a—d R(0)
Ja+d * ] =1
This then permits eqn. A6.16 to be arranged more usefully so that the
two outputs are in terms of the two inputs:
KL=2
KL=2
-5 +5 -5 5L +5
cos(PeL)={(g-j8)/P,) si (A6.18)
One can explore this threshold equation, for example with program
dfbthr in directory grating from:
At threshold with zero power output, the total intensity along the laser
can be shown as being proportional to \a(z) — d(z) 12+ \b(z) 12;
however, this is not of much interest at operating power levels because
spatial-hole burning, with variations in gain and photon density, play
key roles. The field profiles given in Chapter 6 or from computation
using the techniques of Chapter 7 are then more relevant.
356 Distributed feedback semiconductor lasers
0.2
error KL = 2 1
/ / \
0
[/gL =1.8 N
\r \
-2 -1 mode number +1 +2
+5
Figure A 6,2 Error
Error is zero where oscillation modes appear. This is shown for KL=2
and one can see here four modes. Two modes with gain gL~ 1 with ±1
modes while, if the gain is increased to gL~ 1.8, these modes closest to
the stopband grow rapidly with time and the ±2 modes further away
from the stopband just reach threshold
A6.1 Reference
1 GOLDSTEIN, H.: 'Classical mechanics' (Addison-Wesley, Reading, 1980),
2nd ed., p. 156 and exercises on pp. 185-186
Appendix 7
Kramers-Krdnig relationships
A7.1 Causality
The Kramers-Kronig relationships provide fundamental rules for the
relationships between the real and imaginary parts of the spectrum of
any real physically realisable quantity when expressed as a complex
function of frequency. One cannot, for example, design the optical
gain spectrum to be any desired function of frequency without
discovering that the phase spectrum is then closely prescribed.
Frequently, such connections appear to be abstract and mathemat-
ically based. This appendix looks at three different ways of discovering
these relationships which should help the reader to understand the
fundamental nature and physics of the Kramers-Kronig relationships.
The appendix includes at the end a collection of approximations for
the real refractive indices of relevant laser materials.
The first way relies on a real 'causal' system where no real output
o{ t) can occur from any system until after a. real input i(t) has occurred.
An elegant way to proceed splits this real output o(t) from such a system
into even and odd parts, as in Figure A7.1, noting that
oft) oo(t)
Now assume that the Fourier transform of o(t) is Or((o) •\-jOi(<a) in the
spectral domain so that the Fourier relationships give:
«.(')= I \Or(
(A7.2)
(A7.3)
Using eqn. A7.1 links the real and imaginary parts of the spectrum [1]
from
Or(a>) =
lit 0,-(<w') sin(6>' I/I)}
do/
2TT
cos(wif)
I
df (A7.4)
d(of
0,-(«) = cos(a)ft)} sin(a>UI) d^ (A7.5)
2TT J
The fact that cosine and sine are Hilbert-transform pairs permits one
to show that ( ^ - ^ ( O , ) and Oz=3€ (Or) where it corresponds to
taking the Hilbert transform in frequency [2].
To apply this to optics, consider a short length h of optical medium
with a real gain/loss g((o)h and imaginary phase shift S(f>(a))h as
functions of frequency, and an input I(a)) (plane waves with frequencies
co) with a corresponding output 0(a)) related to I((o) by
0(co) ==/(a)) exp h{g(a)) +j8(f>(co)} or
H I +g((o) h+jS<f>((o) h}I(w) (A7.6)
Splitting into real and imaginary parts:
Or(co)={l+g(w)h}I(a)) and Of.(a>) = &|>(<»)A/(a>) (A7.7)
If one allows a delta function or real impulse of field as an input in the
time domain with /(a») = l, there cannot be any real output o(t) for
Appendix 7 359
t<0. This therefore means that the results of eqns. A7.4 and A7.5 have
to hold again and force an intimate connection between the gain and
phase as a function of frequency. These relationships hold for any
physically realisable outputs, of which an optical material's complex
gain/phase is but one such quantity.
The next step is to take the real and imaginary parts so, with the
contribution around the semicircles tending to zero, one has
360 Distributed feedback semiconductor lasers
complex
CO co'+./a
-R, plane
x(a>r)
. da)' = O
XX<o)=-- if
TT - . ( « ' - (o)
d(o' (A7.9)
(A7.10)
Provided that one uses a physical model for a system, then one finds
that the above difficult mathematics can be circumvented and one
always satisfies the Kramers-Kronig relationships. For example, if each
electron were 'bound' to a 'centre' with some damping and restoring
force, then a physical equation of damped harmonic motion of an
electronic dipole (see Appendix 5) is developed:
Appendix 7 361
d x d^c
7-2 +lyl —+n 2 x= ? £(<)/m* (A7.ll)
Of (Xt
x loss (©j
-2.5
real susc.
loss: imag. susc
0.02
where
•n=ifZi/V&0(%\-%'§) (A7.18)
and
« = 1.2398/ A (A7.19)
with % in electron volts and A in micrometres. %0 and %d are parameters
dependent on the composition of the alloys. This expression is valid
provided that the energy of the incident radiation is less than the
bandgap energy; in other words, the material is transparent. Hence, it
cannot be used for the active region material in a laser which, in any
case, will also be influenced by the injected-carrier concentration. We
can now detail the specific parameter values applicable to the
following alloy material systems:
(i) In^Ga^ASjP^ lattice-matched to InP [7,8]
^=1.35-0.72)H-0.12/ (A7.20)
*=0.1894y (0.4184 - 0.013);) (A7.21)
where
(A7.22)
364 Distributed feedback semiconductor lasers
A7.5 References
1 LEE, Y.W.: 'Statistical theory of communication' (John Wiley, 1960)
2 DEUTSCH, R.: 'System analysis techniques' (Prentice Hall, Englewood
Cliffs, 1969), p. 108
3 ERDELYI, A. (Ed.): 'Tables of integral transforms, Vol. 2' (McGraw-Hill,
1954) Use partial fractions and entries 6 and 7, Chap. XV, p. 244
4 ASHCROFT, N.W., and MERMIN, N.D.: 'Solid state physics' (Saunders
College, Philadelphia, 1976), pp. 539-544
5 GREENHOW, M. 'High-frequency and low-frequency asymptotic conse-
quences of the Kramers-Kronig relations' /. Eng. Math., 1986, 20, pp.
293-306
* GREENE, P.D.: Private communication.
Appendix 7 365
where Pamier and Pnohe give the photon power for the optical carrier
(which is being digitally modulated) and optical noise, respectively.
The current is matched into a resistor so that the ratio of equivalent
electrical carrier power to noise power is
^^elechicalequhmlmr (^ramVr) V (Pnoise)^ (A8.2)
2
where (P) indicates a temporal average. The electrical relative-
intensity noise is the inverse, so that one writes
Appendix 8 367
probability
P
false decision
probability
p(1, given 0)
time
Figure A8.1 Importance of relative-intensity noise
(a) Shows an ideal digital signal
(b) Photon detection within a decision window varies with a Gaussian
distribution
(c) Shows carrier power plus noise power. Decisions are made as to
whether one or zero has been transmitted according to a detection of
photons above or below the decision level within the decision window of
time
(d) Noise peaks give a probability of a false decision being made so that
a 1 instead of 0 is detected. Scale of both noise and probability are
enhanced manyfold to illustrate the problem
(A8.3)
detected as a 'one' because the signal is above the threshold level for
long enough within the decision window to return a 'one' state. One
can increase the decision level until there is negligible probability of
this happening; however, one then has to consider the probability of
the detector returning a 'zero' when a 'one' has been transmitted. The
decision level has to be chosen to minimise both probabilities and
need not be exactly at the 'midpoint' as indicated here.
The probability of a false decision giving a 'one' is indicated by the
small shaded area under the probability curve in Figure A8.1d
indicating the integral over the relevant section of the probability
distribution of the signal. A useful model is to take this probability
distribution to be Gaussian (Figure A8.1&) with a variance which is
determined by the mean-square deviation of the photon count within
the decision window given by Seount:
~KJ~2=cr2 (ASA)
The Gaussian probability function p(Smml) for finding a photon level
Scoum given a mean level Scount0 is then
p(SC(mn() dS^d/Vv) exp{- (S^-SM
(A8.5)
Taking for simplicity Smunl0 = 0 for a zero, and take the decision level to
be SdedsionJ then the probability of the signal rising above the decision
level when a zero has been transmitted is, say, pom^wn7jno
In general, the decision level is adjusted to give the minimum value for
pnf>[. The analysis is simplified here by assuming that the noise at both
Appendix 8 369
the 'one' and 'zero' levels is the same and the decision level is adjusted
to the midlevel so that Qj=Sone/ (a 2v2) where Sone is the maximum
output-photon count. For most digital systems, a maximum value of
this false probability is usually taken as 10"9 or even smaller, so that the
probability of a false decision is then
iaP!
mIF ill
-165
H InWISH
-170
i -170
-175
WILJ
-20 -10 " f r 0
frequency, GHz
f
r 10 20
-175
1 ! ! i Iff
-20 -10 "*r0 f
r 10 20
frequency,GHz
A8.1 References
1 COLDREN, LA., and CORZINE, S.W.: 'Diode lasers and photonic
integrated circuits' (Wiley, 1995), pp. 230-241
2 GREEN, P.E.: 'Fibre optic networks' (Prentice Hall, Englewood Cliffs,
1993), chap. 9 gives an overview about the relative importance of different
forms of noise on different systems
3 GOLUBEV, Y.M., and SOLOKEV, I.V.: 'Photon anti-bunching in a coherent
light source and suppression of the photorecording noise', Sov. Phys.-JETP,
1984, 60, pp. 234-238
4 MACHIDA, S., YAMAMOTO, Y, and ITAYA, Y: 'Observation of amplitude
squeezing in a constant current driven semiconductor laser', Phys. Rev.
Lett., 1987, 58, pp. 1000-1003
5 RICHARDSON, W.H., MACHIDA, S., and YAMAMOTO, Y: 'Squeezed
photon numbtr noise and sub-Poissonian electrical partition noise in a
semiconductor laser', Phys. Rev. Lett., 1991, 66, pp. 2867-2870
6 LI, Y-Q., EDWARDS, P.J., LYNAM, P., and CHEUNG, W.N.: 'Quantum-
correlated light from transverse junction stripe laser diodes', Int. J.
Optoelectron., 1996, 10, pp. 417-421
Appendix 9
Thermal, quantum and numerical
noise
A9.1 Introduction
Now recall that detection takes the signal down to baseband where the
noise bandwidth runs from — 8f to + 5 / and so the net bandwidth
A/=2 8f so that eqn. A9.6 confirms the result of eqn. A9.3.
With coherent optical signals, it is possible to recover a factor of 2 by
using homodyne detection which requires that the phase of the output
is known. This additional information about the phase reduces the
minimum uncertainty to i A/A/rather than the value of A/A/obtained
by power (photon) detection. The minimum uncertainty in measure-
ment of the optical-field power then directly relates to the best possible
signal-to-noise-power ratio, given that one uses ideal measurement
methods.
1 dE dE
_ + —
i (t,z) (A9.7)
Vff dt dz
Here ispf is proportional to the random dipole 'current' as a function
of time that excites the spontaneous emission for the forwardfieldand
g is the field gain, at present taken to be a constant. To simplify the
analysis, redefine the field variables:
Appendix 9 375
(A9.9)
[<
F=Fin exp(gz) + exp(gz) I ispf(t, zr) exp( — gzf) dzf (A9.11)
Jo
Squaring up on either side and ignoring terms like F*isfwhich average
to zero
F*F=Fin*Finexp(2gz) +
, * f'^f- , , ,- ,
p y n ( y fry \ I fl 7 I 1 (t 7 i ^ PVTII — (T7 I 1 ( t 7 \ PYT1 ( — GT7 i fl7
Jo Jo
(A9.12)
Spontaneous excitations are assumed to be spatially uncorrelated,
such that h,,f(t, f)% (t, z')=r (t)28(z' - z") (A9.13)
f spJ
time gives the mean-square noise power per unit length of amplifier.
Integration over the length L gives the output at the end of the
amplifier region from
For future comparison, note that the rate of change of the mean power
flow given by eqn. A9.15 can also be written at a position z as a sum of
stimulated growth plus addition of spontaneous emission:
dPf/dz = 2gPf+fspf(t)2 (A9.16)
Writing GL=exp(2gL) as the power gain along the guide gives
(A9..7,
(A9.18)
Now consider the interesting situation. At the input there was a power
Pin but with a quantum uncertainty of the power (using direct
detection) given from h/5/over a bandwidth of Sf. However, this
uncertainty is not power that can get amplified but is merely an
uncertainty in the measurement of the fields. At the end of an ideal
quantum-amplification process with a sufficiently large gain GL> the
quantum uncertainty can be negligible compared with output power
levels, but now it appears that there is an added available noise power
of fsp/(t)2/2gproduced by amplifying the available thermal noise. The
measurements at the input and the output are ideal measurement
processes and, if the amplification is ideal quantum amplification with
no absorptive loss of photons but only quantum spontaneous emission
along with ideal stimulated gain, it is argued that the signal-to-noise
ratio at the output must equal the signal-to-noise ratio Pfin/ (hfSf) at
the input. Quantifying this result gives
2gPfm/fspf(tf=P/in/hf8f or fspf(t)2=2gh/Sf (A9.19)
Here then is a statistical measure of the mean-square spontaneous
Appendix 9 377
PfMU=APfin+ (1 - A) - ^ - (A9.20)
- ^ - = h / A , / / { e x p ( h / / k r j -1} (A9.22)
The transition from attenuation to gain can be made via the concept
of a negative reciprocal temperature where ( 1 / T J is large near
absolute zero but becomes small by pumping sufficient electrons and
holes into the interaction region. This pumping initially effectively
increases their equivalent temperature with the attenuation a tending
to zero until transparency is reached when a- 0. Then with still harder
pumping, a changes sign to become gain g( = — a) with
(1/ Tw)—• — 1(1/ 7JI and with yet harder pumping still to give large
enough 'population inversion' (1/TJ—• —oo with exp(h//kjTJ—»0
and one recovers eqn. A9.17. The conclusion that one arrives at is that
378 Distributed feedback semiconductor lasers
for the material. The connection between jR^, and KsHm was found from
eqn. 2.19 considering the equilibrium conditions and black-body-
radiation density to give eqn. 2.20 written as
lvg) (A9.28)
On rearranging eqn. A9.27 with the result of eqn. A9.28, explicit
values for the spontaneous excitation and spontaneous-coupling
coefficients are obtained as
(A9.29)
effective
erreff aperture
optical guide
(A9.33)
where spontaneous excitation ispf(Ts/vg, Zs) is inserted at each of the
space points [T, Z] which are separated by s so that Z and T are integers
defining the step numbers. The lack of correlation between each space
Appendix 9 381
and time point requires for distinct Z and Z' and distinct T and T':
sispf(Ts/vg, Z's)sispf(Ts/vg9 Zs) =0 (A9.34)
sispf(Ts/vr Zs)sispf(Ts/vg, Zs)=0 (A9.35)
If one wishes to restrict the noise bandwidth below the Nyquist limit,
one should introduce temporal correlation, as will be discussed
shortly.
Now similar results can be formed for the reverse wave: there is
nothing special about the forward direction other than the sponta-
neous excitation for the forward wave is uncorrelated with the
excitation for the reverse wave. Hence one has
sispr(Ts/vg, Ts)sispr(Ts/vg, Zs)=0 (A9.36)
sispr(Ts/vgf Zs)stapr(Js/vg9 Zs)=0 (A9.37)
with
sispr(Ts/vg, Z's)sispf(Ts/vg9 Zs)=0 (A9.38)
but now for any T, T, Z and Z\ even though these may be equal. This
ensures that there is no coupling between the forward and reverse
emission when there is no scattering from the forward field into the
reverse field.
When space and time are equal in both excitation terms
sisp/T/vg, Zs)sispJ(Ts/vg, Zs) =srspf(t)2 (A9.39)
2
sp (T/vg, Zs)sispr(Ts/vg, Z s ) = s V ( 0
siispr (A9.40)
with
(A9.41)
As the length of the laser waveguide L=N s increases (retaining g=0
for this calculation), the mean spontaneous forward power increases
as
l (A9.42)
These requirements are met for the forward field by inserting, at each
space step and at every time step, a random excitation given by
sispf(Ts/vg, Zs)={l(a/47r) (Sf/AJ) (rsps&sBPNhf)\l/2gi(A9A3)
where S?l=9?+7% with % and °}J normally distributed random real
numbers such that the mean-square value 2ft 2 = 1. The randomness of
the complex 91 gives the correct zero-mean correlation in amplitude
and in phase from point to point and time to time, as required for
382 Distributed feedback semiconductor lasers
spontaneous power
(a) I
frequency
spnum
spontaneous power
(b)
frequency
spontaneous power
(c)
frequency
[
physical ] numerical! 1111II
spontaneous
L
power
frequency
(A9.48)
and s is sufficiently short, the numerical-noise spectrum, may extend
beyond the physical-noise spectrum, as in Figure A9.3, but could
approximate the significant central region adequately.
A9.9 References
1 HENRY, C.H.: 'Theory of spontaneous emission noise in open resonators
and its application to lasers and optical amplifiers',/. Lightwave TechnoL,
1986, 4, pp. 288-297
2 MARCENAC, D.D., and CARROLL, J.E.: 'Quantum mechanical model for
realistic Fabry-Perot lasers', IEEProc.-J., 1993, 140, pp. 157-171
3 THOMPSON, G.H.B.: 'Physics of semiconductor laser devices' (Wiley,
1980)
4 FREEMAN, J.J.: 'Principles of noise' (Wiley, 1958)
5 KITTEL, C, and KROEMER, H.: Thermal physics' (W.H. Freeman and
Co., 1980), 2nd ed.
6 MARCUSE, D.: 'Principles of quantum electronics' (Academic Press,
1980)
7 GHAFOURI-SHIRAZ, H.: 'Fundamentals of laser diode amplifiers' (Wiley,
Chichester, 1995)
8 RAMO, S., WHINNERY, J.R., and VAN DUZER, T: 'Fields and waves in
communication electronics' (Wiley, 1994), 3rd ed.
Appendix 10
Laser packaging
A10.1 Introduction
The extra requirements of packaging optoelectronic rather than just
electronic devices are challenging [1], and real packages involve
numerous compromises between different aspects of performance, as
well as manufacturing costs. To make the best of a well designed
device, DFB-laser chips must be packaged so as
(i) to protect them from the environment;
(ii)to enable monitoring of the optical output;
(iii)
to provide thermal heat sinking;
(iv)to provide a rigid and stable optical coupling platform to a fibre
(or other output system); and then
(v) to allow electrical interfaces between the appropriate 'DC and
modulated RF current supplies and the laser.
V V,n
driver"1' ' • ' • i . • - • • • •# . -*•— '•
ensure that more power emerges from the front facet). Power for the
fibre is extracted from the front facet [often antireflection (AR)
coated]. In Fabry-Perot lasers, the front and back powers track each
other very well, but in DFB lasers spatial-hole burning and longitudi-
nal shifts of the lasing power can result in considerable divergences
(see Figure 5.18). In addition, efficient coupling of optical power from
laser to fibre requires submicron tolerances, and is likely to vary with
temperature, whereas power coupled to the monitor photodiode is
stable due to the relatively large sensitive area. Hence monitoring of
power in the system fibre is likely only to be accurate to ± 1-2 dB, and
long-term 'creep' of the package may increase this figure. Caution is
therefore needed in analogue or other high-linearity applications in
using the monitor diode to control the laser's precise performance.
A10.6 References
1 MATTHEWS, M.R., MACDONALD, R.M., and PRESTON, K.R.: 'Optical
components—the new challenge in packaging', IEEE Trans. Compon.,
Hybrids Manuf. TechnoL, 1990, 13, pp. 798-806
2 XENG, J-X.: 'A 3D electromagnetic simulator for high frequency
applications', IEEE Trans. Compon., Packag. Manuf. TechnoL B, Adv. Package
1995, 18, pp.578-595
3 TUCKER, R.S., and POPE, D.J.: 'Microwave circuit models of semi-
conductor injection lasers', IEEE Trans., 1983, MTT—31, pp. 289-294
4 DELPIANO, R, PAOLETTI, R., AUDAGNOTTO, P., and PULEO, M.:
'High-frequency modeling and characterization of high-performance
DFB laser modules', IEEE Trans. Compon. Packag. Manuf. TechnoL B, Adv.
Packag., 1994, 17, pp. 412-417
5 TIWARJ, S.: 'Transmission line delay limitations of laser bandwidths', IEE
Proc. Optoelectron., 1994, 141, pp. 163-166
6 NEIFELD, M.A., and CHOU, W-C: 'Electrical Packaging Impact on
Source Components in Optical Interconnects', IEEE Trans. Compon.,
Packag. Manuf. TechnoL, B, Adv. Packag., 1995, 18, pp. 578-595
7 HELMS, J., KURTZKE, C , and PETERMANN, K.: 'External feedback
requirements for coherent optical communication systems', J. Lightwave
TechnoL, 1992, 10, pp. 1137-1141
8 ROBERTSON, M.J., LEALMAN, I.F., and COLLINS, J.V.: 'The expanded
mode laser - a route to low cost optoelectronics', IEICE Trans. Electron.,
1997, E80C, pp. 17-23
9 GREEN, P.E.: 'Fibre optic networks', (Prentice Hall, Englewood Cliffs,
1993), p. 87
10 GOWAR, J.: 'Optical communication systems' (Prentice Hall, 1993), 2nd.
edn., p. 438
11 HALL, S.A., LANE, R., WANG, H.C., and GARERI, A.: 'Assembly of laser-
fibre arrays',/. Lightwave TechnoL, 1994, 12, pp. 1820-1826
12 PEALL, R.G., SHAW, B.J., AYLIFFE, P.J., PRIDDLE, H.F.M., BRICHENO,
T, and GURTON, P.: '1 x8, 8 Gbit/s transmitter module for optical space
switch applications', Electron. Lett, 1997, 33, pp. 1250-1252
13 CANN, R., HARRISON, P., and SPEAR, D.: 'Use of silicon vee groove
technology in the design and volume manufacture of optical devices',
Proc. SPIE, 1997, 3004, pp. 170-173
Appendix 11
Tables of device parameters
and simulated performance for
DFB laser structures
Table A 11.1 Summary of material and device parameters used for large-signal
dynamic modelling of uniform-grating DFB laser
* The static modelling used similar parameters but where they differ from those used
for the dynamic modelling they are shown in brackets
Table All.2 Material and device parameters used for modelling unstrained
(inGa)As-(InGa) (AsP) 6QWy 2x \J8DFB laser structure
s)
source matching bond
R.
\
laser
bond-wire C X i
laser T '
laser
Grating parameters
2 xA m /8 grating, #r=60cnrr 1 , KL=2A, L=350 jxm
Am/8 phase shifts positioned symmetrically in cavity and spaced 245 |jim apart,
i.e. AL=122.5 fxm
and inserts, as in Table A12.1, the entries for interrogating the new
subdirectories. The punctuation marks at the start and ends of the
lines need to be followed in exact detail. MATLAB will then
automatically interrogate these directories/subdirectories on start up
Appendix 12 397
'C:\MATLAB\toolbox\matlab\demos;',...
'C:\MATLAB\laser\dfb;1,...
'C:\MATLAB\laser\diff;1,...
'C:\MATLAB\laser\fabpero;1,...
etc.
"C:\MATLAB\laser\spontan;1,...
above 100 MHz, and it is worth remembering that enthusiasm can die
when results take too long to appear.
Random-access memory (RAM) is important. The MATLAB figures
produced by the code appear on our machines to use very roughly 1-2
Mbyte of memory, so that it is easy to find 16 Mbyte of memory being
used simply in storing a relatively small number of figures, and it seems
that 32 Mbyte of RAM is highly desirable for easy computation. For any
serious future user of this code, 64 Mbyte and 166 MHz clock rate is
recommended, remembering that even at the time of writing this is
already out of date! The M-file delfig, which is in several sub-
directories, permits the rapid deletion of figures to avoid clogging
memory.
Each machine, in our experience, can have slight differences in its
preferences about the directories in which the MATLAB codes are
stored. It appears to be safest to store the working programs actually
within a subdirectory within the MATLAB main directory. However,
above all, store the core programs in a back-up version on your own
machine to avoid continual reliance on the Net. Modify only copied
versions of this back-up and then back-up the modified versions!
Remember that programs are stored with the extension .m (e.g.
slabexec.m) but to run the program slabexec.m, for example, one
only types slabexec at the command-window prompt. The MATLAB
command window must, of course, be open and address the right
directory or it will simply return:
? Undefined function or variable slabexec
Ensuring that one has the right window is best accomplished by
clicking on
File , Run M-file, Browse
and then switching into the required directory and program as per
menu.
This appendix lists, under the different directories, the programs
which have been made available along with a very brief description of
them, and where sample outputs can be found within the book.
Opening the directory will give access to a README file about
programs in the directory. Opening the program as if to modify the M-
file will give access to a short account of the program, along with
'running' comments with the line header given by the symbol % .
Within the programs, there are default settings so that one usually
presses <return> or types return when asked to make an entry from
the keyboard. The default setting is then automatically returned and
400 Distributed feedback semiconductor lasers
the program should run. In this way it is hoped that the reader will at
least get started, but the real exercise for the reader is to both run and
modify the programs.
The contents of the various directories are briefly recounted below,
in approximately the order the programs appear in the book.
slab
•4
3.4 •0.0002467
-4.5
-5 0.1
2 3 4 -40 -20 0 20 40 60
T E E fields (mod squared) TE far fields - angular distribution
Figure A12.1 Sequence of plots obtained from the default mode for running slabexec
grating
/ \^3'4 / I
-1
-1.5 \ 3 , /0.000246A / /
^0.0002467 '
-2 V
-2.5 3.6 0.002467
depth
-3
Figure A12.2 Modified output from slabexec when inserting a Bragg grating
etched into the waveguide with an infill of different material
Note that guide is too wide and the overlap with the grating is
inadequate, leading to a low K
diff
(i) Programs stepr and stepj in directory diff integrate AF/ dz= zFand
dF/ dz = j zF, respectively, showing how step length alters the accuracy.
Figure 7.2 gives a sample of the output expected.
(ii) Program advec provides a basic test program for the advection
equations where the forward field is stored at all space points but for
only two time steps. Figure 7.3 provides a sample of the output
expected.
spontan
(i) Program spont in directory spontan gives the reader an indication
of the modelling process for spontaneous emission. It provides a
dynamic output of the static display, shown in Figure 7.5.
Appendix 12 403
(ii) Programs amp1, amp2 and amp3 provide three short demonstra-
tions of generation and amplification based on an 'ideal' uniform laser
amplifier, demonstrating that spontaneous emission and gain are
linked. Program ampif is similar to amp 1 but has a filtered noise input
rather than relying on white spontaneous noise.
fftexamp
filter
dfb
fabpero
There are two MATLAB programs fpstat and fpdyna in directory
fabpero which are able to give the static and dynamic characteristics of
a Fabry-Perot-type laser using normalised rate equations discussed at
the end of Appendix 4. The outputs are similar to those in Figures 4.2
and 4.3.
A12.3 Reference
1 'MATLAB 1995, Student Edition' (Prentice Hall and The MathWorks Inc,
1995)
Index