Electro-absorption Modulators EE232, Spring 2001 Pei-Cheng Ku April 24 & 26, 2001
Dec 20, 2015
Electro-absorption Modulators
EE232, Spring 2001Pei-Cheng Ku
April 24 & 26, 2001
Optical Modulation
Direct modulation on semiconductor lasers: Output frequency drifts
carrier induced (chirp) temperature variation due to carrier modulation
Limited modulation depth (don’t want to turn off laser)
External modulation Electro-optical modulation (low efficiency) Electroabsorption (EA) modulation (smaller
modulation bandwidth)
Electroabsorption (EA) Modulator
EA modulator is a semiconductor device with the same structure as the laser diode.
In laser diodes, we inject large enough current to achieve stimulated emission. While in EA modulator, we apply electric field (reverse bias) to change the absorption spectrum. No carriers are injected into the active region. However, carriers are generated due to absorption of light.
Advantages of EA modulators
Zero biasing voltage Low driving voltage Low/negative chirp high bandwidth Integrated with DFB
Schematics of an EA modulator
Waveguide type is more popular.
From textbook, figure 13.10.
Outline of our discussion
Know how absorption spectrum changes with applied electric field (effective-mass approximation)
See how extinction ratio of modulation can be enhanced in the quantum well (quantum confined Stark effect). Excitons will be introduced.
We will follow the textbook (chapter 13). A recommended reference book is “Introduction to Semiconductor Optics” by N. Peyghambarian et al. (1993).
Physics behind EA Modulators
How absorption spectrum in semiconductors can be changed? Physical model: effective-mass equation
Single-particle representation Two-particel representation
Coulomb interaction between electrons and holes: Excitons Electric field effect: Franz-Keldysh effect Coulomb+Electric field: DC Stark effect Coulomb+Electric field+QW: QCSE
Absorption Spectrum Change under Applied Electric Field
Franz-Keldysh Effect Neglect Coulomb interaction between electrons
and holes. DC Stark Effect
Franz-Keldysh effect plus Coulomb interaction between electrons and holes (excitons).
Quantum confined Start effect (QCSE) DC Stark Effect in quantum wells Excitons been confined in quantum well. Stark
effect enhanced.
Effective mass approximation Electron moving in semiconductor crystal with periodic potential V (r)
under the influence of a potential U (r).
Effective mass approximation: if the potential is slowly varying compared to the period of the lattice, the above equation can be replaced by effective-mass equation as follows (section 4.4.1 and appendix B):
that is the potential causes the envelope of the fast-varying Bloch function to change. Without potential U (r), the envelope (r) reduces to plane wave
)()()()(2 crystal
2
0
2
rrrr
EUV
m
)( )()( )()0()()(2
2*
2
rrrrkrr kk nn uEEUm
*
22
2 and )exp(
m
kEi
rk
Wavefunction for electrons moving in periodic potential
For example, at the valence bandedge. The wavefunction is p-like (L = 3/2 in sp3 orbitals).
)(rkvu
)exp()( rkr i
)()()( rrr u
Form of U (r)
Franz-Keldysh Effect
DC Stark Effect (slow-varying approximation only holds for Wannier excitons)
Quantum confined Start effect (QCSE)
rFr eU )(
)(4)(
2
heee rrrFr
r
eeU
)()(4
)(2
ehe
ee rrr
rFr er
Ve
eU
Effective mass approximation in Two-Particle Representation
Since U (r) is in general a function of both electron and hole, it’s more convenient to combine two single-particle equations of motion for electron and hole together.
Effective-mass equation for electrons:
Effective-mass equation for holes:
Add them together
)()0()()(2
2*
2
eee rkrr eeneeeee
EEUm
)()0()()(2
2*
2
hhh rkrr hhnhhhhh
EEUm
)()()()()(),( is ion wavefunctoriginal
and )()()(),( where
),(),()()(22
,,,,
,,
2*
22
*
2
hkekhkk,k
ekehe
hkk,k
ekehe
hehehe
rrrrk,krr
rrk,krr
rrrrrr
heh
he
e
h
he
e
vcheh
ehheh
ghehh
ee
uuB
EEEB
EEUUmm
Effective mass approximation in Two-Particle Representation (cont.)
In the path of generalizing single-particle effective mass equation to two-particle effective mass equation, we didn’t make new approximations. The new envelope wavefunction (which is a linear combination of the single-particle envelope functions for electrons and holes) still has to satisfy the slowly-varying approximation. This will be valid if the potential seen by electrons can be smoothly transferred to the potential profile seen by holes.
V
e
V
eA
kk
B
ii
he
he
hh
he
ee
h
he
e
rk
k,k
rk
hehe
hkk,k
ekhehe
k,krr
rrk,krr
)(),(
follows. as and of modes waveplane of expansion thein writtenbe also can
)()()(),( ,,
Validity of Slowly-Varying Approximationin two-particle representation
Franz-Keldysh Effect
For 5V applied voltage, there is 500kV/cm electric field across 0.1m active region. The variation of electric field across a lattice site is 500kV/cm x 5.6Å ~ 0.028V which gives an energy variation of 0.028eV for an electron.
For an electron moving in periodic atomic potential, it sees a Coulombic attraction from the nucleus.
The potential energy variation across a lattice site is around V(a) ~ 2.569eV >> 0.028eV!
|| )(),( eeeU hehe rrFrr
nucleus.) screened"" seesIt atom. each from available is electron valenceone (Also,
used.) been has of instead (Note, 4
)( 00
2
rr
erV
Excitons
The effective-mass equation for an electron-hole pair is:
We have to sum over all possible interactions between electrons and holes. But if the carrier density is low enough, we can neglect the contribution from the rest of the electrons and holes and consider only bare Coulomb interaction, that is:
n)interactiopair of because 1/2factor a put to(remember
||4
1
2
1)()( where
),(),()()(22
, ,
2
2*
22
*
2
ji jhe,irhe
he
ghheehh
ee
e
NNUU
EErUrUmm
rrrr
rrrr
he
hehe
||4)()(
2
hehe rrrr
r
he
eUU
Effective-mass equation for Excitons
With the bare Coulomb interaction between electron and hole, we have the effective-mass equation for an exciton as follows.
This equation is formally identical to a hydrogen atom problem except the masses for the electron and proton are replaced by the effective masses for the electrons and holes, respectively.
Any central-force problem (i.e. the potential is a function of relative coordinate between two particles) can be separated into center-of-mass motion and relative motion.
),(),(||422
22
*
22
*
2
hehehe
rrrrrr
gr
hh
ee
EEe
mm
),(),(
:mass ofCenter :coordinate Relative**
**
rRrr
rrRrrr
he
hehe
he
he
mm
mm
Center-of-mass motion and Relative motion
Effective-mass equation in two-particle representation for any central-force problem can be rewritten as follows.
Center-of-mass motion is a free-running solution:
),()(),()(22
22
22
rRrRrrR
g
r
EEUmM
M
KEE
Um
V
e
g
r
i
2
)()()(2
)(),(
22
22
rrr
rrR
r
RK
Solutions for Excitons Same as hydrogen atoms (3D or 2D) Exciton Bohr Radius and Rydberg energy
Energy levels: bound states and continuum
Ionization energy for ground state exciton in 2D case is 4 times larger than it is in 3D case. Therefore in QW, it’s easier to form bound exciton states.
22
4
2
2
)4(2
4
r
ry
r
rB
emR
ema
2
2
)2/1( :(2D) states boundFor
:(3D) states boundFor
n
RE
n
RE
yn
yn
Energy Diagram of ExcitonsElectron-hole energy E
n=1
n=2
n=3
K
In GaAs,
E1=4.2meV (3D) or
E1=16.8meV (2D)
continuum
Wavefunction for Excitons (3D)
Coulomb potential is spherical symmetric with respect to the center of mass.
Three quantum numbers: n, l, m.
“Accidental degeneracy” for same quantum number n. Orbitals designated by 1s, 2s, 2p, 3s, 3p, 3d, …
),(Y)()( lmnl rr
llllm
nl
n
,1 , 0, , ,1 ,
1 , 0,
3, 2, 1,
Radial Solution (3D)
.polynomial Laguerreassociated theis where
2exp
])!1[(2
)!1(22)(
12
12
2/1
3
3
qL
na
rL
na
r
nn
ln
nana
rr
lln
B
lln
BB
l
Bnl
Bar
B
s ea
r
31
1)( e.g.
2 4 6 8 10
0.002
0.004
0.006
0.008
0.01
0.012probability
r/aB
2p
2s
1s
Spherical Harmonics Ylm (3D)
|| 00Y || 20Y|| 11Y
|| 10Y || 12Y
|| 22Y
Wavefunction Representation for Excitons (3D)
s-like electron (L=1/2) bloch function
p-like hole (L=3/2) bloch function
Sitting at holereference frame.
S-like exciton envelope function (e.g. n=1)
)(rP-like exciton envelope function (e.g. n=2)
)(r
Wavefunction for Excitons (2D)
Coulomb potential is cylindrical symmetric with respect to the center of mass.
Two quantum numbers: n, m.
“Accidental degeneracy” for same quantum number n.
2)()(
im
nm
er
1 , 0, , ,1
3, 2, 1,
nnm
n
Radial Solution (2D)
.polynomial Laguerreassociated theis where
)2/1(
2
)2/1(exp
)!1||()12(
)!1||(
)2/1(
24)(
||21||
||21||
2/1
3
||
qL
an
rL
an
r
mnn
mn
an
r
ar
mmn
B
mmn
B
m
BBnm
Bar
B
ea
r2
10
4)( e.g.
Ba
r
B
s ea
r
31
1)( tocompared
2D excitons are more tightly bounded.
Validity of Approximations used in the Calculation of Excitons
Slowly-Varying Potential Profile Approximation:GaAs/InGaAs/AlGaAs lattice density ~ 5.6Å
Bare Coulomb interaction approximation:Exciton radius aB for GaAs ~ 140ÅExciton density:
3D: 8.7e16/cm-3 2D: 8.1e10/cm-2
For EA modulators, no carriers are injected into active region. Carrier density ~ intrinsic doping level ~ 1e16/cm-3 (bulk) or 1e10/cm-2 (100Å QW). Bare coulomb interaction is a good approximation.
But for laser applications, carrier density is high: 5e18/cm -3 (bulk) or 5e12/cm-2 (100Å QW). Coulomb interaction between electrons and holes is screened.
potential) s(Yukawa' )exp(
4)(
2
r
rerU
s
Optical Absorption without U (r)
From Fermi’s golden rule, the absorption coefficient can be derived as in equation (9.1.23) between two energy bands:
where matrix element is calculated as in section 9.3.1 in the case of no external electric field and neglecting the Coulomb interaction between electrons and holes. The result is as follows (equation 9.3.11).
))((22
)(2'
cvvccvk k
ffEEHV
v c
))((ˆ2
)(
ˆ2
2
0
,0
0'
cvvck
cv
kkcvcv
ffEEpeV
C
pem
eAH
vc
'cvH
Optical Absorption with U (r)
unctioninit wavef|ˆ|ion wavefunctfinal2
:becomesnow element matrix interband The
0
0'cv
icv pee
m
eAH op rk
)()()()()(),( is ion wavefunct that theRecall ,,,, hkekhkk,k
ekehe rrrrk,krrheh
he
e vcheh uuB
V
e
V
eAB
ii
he
hh
he
ee
h
he
e
rk
k,k
rk
hehkk,k
ekhehe k,krrk,krr
)()()()(),( function envelope thewhere ,,
.|)0(||)0(| problem, force central In
)(|)0(||ˆ|E2)( spectrum absorption
)0(ˆ)(eˆ
dˆ)(
22
statespossible all
220
'
rr
r
rk,k
rk,k
e
hh
he
ee
kee
RK
rkrk
k,k
rk
he
EpedC
VpeApe
eV
e
V
epeAH
cv
cvi
cv
iii
cvcvop
Absorption Spectrum Calculation Summary
Absorption is proportional to the probability that we find electron and hole at the same “position” (i.e. r=0).
For central force problem, we only need to evaluate the envelope wavefunction for the relative motion and take the square of its absolute value.
Note the factor of 2 in the formula of the absorption spectrum comes from the fact that we consider two particles at one time.
The delta function takes care of the (joint) density of states available.
We have to sum over all possible states to get the total absorption.
Remarks on the first assertion
Same “position” actually means electron and hole are in the same lattice site. What happens in the unit cell is governed by the oscillation strength pcv which depends on the details of the Bloch functions inside the lattice site.
This is because is an “envelope” function. The true wavefunction is the product of the envelope function with the Bloch function.
is the probability of finding the electron and hole in the same lattice site.
2)0( r
),( he rr
2)0( r
Density of States function g (E)
Density of states means total number of degenerate states for a given energy level. For example in bulk materials, since there’s no restriction on the motion of the electrons, for a given energy E, the number of states that satisfy
is the density of states at that energy.
)(2
222*
2
zyx kkkm
E
)()( d)(d
statespossible all
EgEfEEfE
g (E) as a function of dimensionality
EEEm
kkk
Em
E
kEkkkk
EEE
m
kE
EE
kEk
EEEk
EkV
r
r
r
n
g(E) d2
2d4d :3n
constantg(E) d2
d
dd2d2d :2n
/1g(E) d1
2dd
dd
ddd :1n
)(g(E) d)(d :0n
g(E)d d)2(
22
3/2
223
22
2
n
n
k
Diagrammatic representation for g (E)
3D
1D
E
g(E)
2D
0D
For a given energy range, the number of carriers necessary to fill out these density of states: 3D>2D>1D>0D.
Franz-Keldysh Effect Franz-Keldysh effect is also a central-force problem.
2/33/1
2/33/1
2/33/1
2/33/1
2
2
3/1
2
2
22
3
2
3
2
3)()(
3
2
3
2
3)()(
0for
0)(d
)(d
2 variableof change
)()(2
ZJZJZ
ZAiZ
ZIZIZ
ZAiZ
Z
ZZZ
Z
eF
Ez
eFmZ
zEzeFzdz
d
m
zr
zr
|| )(),( eerreFrrU hehe
Airy Function Ai (Z)
Z>0: electron-hole energy+Eg < electric field potential
Z<0: electron-hole energy+Eg > electric field potential, i.e. above bandgap oscillation wavefunction
-10 -7.5 -5 -2.5 2.5 5
-0.4
-0.2
0.2
0.4
Z
Ai(Z)Smaller period
Absorption Spectrum of Franz-Keldysh Effect
Absorption spectrum reduce to the familiar square root dependence of energy when F0.
Exciton Absorption Spectrum
3D Excitons
2D Excitons
factor.t enhancemen Sommerfeld is 1
/2)(
where
114
2)(
/23
132332
0
xD
y
g
n y
gD
y
g
By
e
xxS
R
E
R
ES
nR
E
naR
A
factor.t enhancemen Sommerfeld is 1
2)(
where
)2/1(
1
)2/1(
14
2)(
/22
12232
0
xD
n y
gD
y
g
By
exS
R
ES
nR
E
naR
A
Exciton Absorption Spectra (schematics)
2D and 3D exciton absorption spectra with zero/finite linewidth. The spectra is very sensitive to the temperature.
With dampingWithout damping
Excitonic effect at different temperature
T Bulk GaAs
Quantum well (quasi 2D)
Absorption peak blueshifts because the bandgap becomes larger at lower temperature.
Sommerfeld Enhancement
Even without excitonic peaks, bandedge absorption is enhanced due to Coulombic interaction between electrons and holes.
Lots of bound states nearthe onset of continnum sumtogether to give Sommerfeldenhancement.
Quantum Confined Stark Effect*
h2
*
2
e2
*
2
2
||2
||2
where
),()(),(4
Vem
HVem
H
EEe
HH
hh
h
he
e
e
hghhr
he
rFrF
rrrrrr
e
eee
Quantum well potential
Because the potential of the quantum well breaks the spherical symmetry, it’s not a central force problem anymore! This symmetry breaking, unfortunately, forces us to solve this problem numerically. However, we’ll try to factorize this equation and get more physical insights.
*D. A. B. Miller et al., Phys. Rev. B., 32, 1043 (1985)
Factorization of Schrödinger equation for QCSE
),,(2
),,()()(4d
d
2d
d
22
follows. as rewritten be can equationer Schroeding original Now the
.coordinate mass ofcenter on dependt doesn' function potential thebecause
),,()(
as factorized be also can )( function envelope The
. and ˆˆ where
and
:planexy in scoordinate relative and mass ofcenter Define
22
22
2
2
2
*
2
2
2
*
22
2
**
**
het
ghehhee
rhheer
he
i
he
he
zzGM
KEEzzGzVzV
z
e
zmzmm
zzGA
e,
,
mmMyx
M
mm
ρρ
ρrr
rr
yxρ
ρρRρρρ
tt RK
he
he
hethe
Solutions to QCSE
Consider solutions inside the quantum well. For given z, we can see that the problem reduces to 2D-exciton problem. For given , this is a quantum well problem. Therefore we can expect the solution to this equation can be expanded in terms of the product of 2D-exciton envelope function and the quantum well wavefunction. We then expect that this expansion should only need a few terms to give a good approximation.
Physically, because of the carrier confinement due to the quantum well, we should expect most of the excitons have their electron and hole either both in the quantum well or both outside the quantum well.
As we have seen before, 2D exciton has 4 times larger ionization energy and therefore is much more stable than excitons outside the quantum well.
Absorption Spectrum for QCSE
Exciton absorption peak red-shifts with increasing electric field strength.
Absorption peak shifts quadratically with applied electric field
If the quantum well “only perturbs” the solution a little bit, the absorption peak position shift can be calculated from time-independent perturbation theory.
The energy shift, for example for the ground state exciton in infinite quantum well, due to the electric field is:
Absorption peak shifts quadratically with applied electric field.
2422**1
11
11
)0()0(
2221
111
/)()0()(
0 eh1||eh1 order termfirst The
|keh,||jeh,|eh1||eh1)0()(
effheeh
eh
kj jk
eh
eh
LFemmCEFE
z
EE
zFezeFEFE
Potential profile has reflection symmetry along z-direction. Thus exciton wavefunction is either evenor odd with respect to z=0.
Design considerations for EA modulators
Operation principle Contrast ratio Insertion loss Modulation efficiency Chirp considerations and optimization Packaging and Integration
EA Modulators use QCSE
By applying electric field parallel to the quantum well growth direction, the absorption can be changed dramatically at the desired wavelength.
Contrast Ratio
The larger the contrast ratio, the better it is. Contrast ratio can be made as large as possible
by increasing the length of the modulator. But propagation loss then becomes an issue.
LVRdBR
e
e
VVP
VPR
offonoffon
LV
L
offout
onoutoffon
)0()(343.4)log(10)(
)(
)0(
//
)(
)0(
/
Insertion Loss
The longer the modulator is, the larger the insertion loss. Therefore, contrast ratio should not be increased by lengthen the device too much.
Waveguide structure also needs careful design to increase the coupling efficiency from the single-mode-fiber output.
L
in
outin eP
VPP )0(1)0(
Loss
Modulation Efficiency Modulation efficiency quantifies how much voltage do we
need to modulate the optical signal.
Smaller detuning will increase the modulation efficiency. However, it also results in a larger insertion loss.
FV
LV
V
R offon
343.4
)0()(343.4/
Chirp
Usually, the light to the input of the EA modulator comes directly from the output of a laser diode (e.g. DFB). This light signal is almost single frequency but will have a finite linewidth due to the frequency noise inside the laser diodes. This variation of frequency is called “chirp”.
Chirp results from spontaneous emission and carrier induced phase fluctuation from the coupling between imaginary and real parts of the refractive index (Kramers-Kronig relation).
spp
eFWHM
e
RN
f
dNdg
dNdn
dNn
dNn
)1(
/
/4
/]dIm[
/]dRe[
2
'
Fundamental linewidth
Chirp generated from Direct Modulation
Recall that the generalized rate equation for laser diodes we derived in the class:
The amount of chirp generated from direct modulation depends on e and the signal waveform (note that S is the photon number inside the cavity and the output photon number is S (1-R) where R is the reflectivity of the output facet).
Large-signal modulation: blue-shift at the beginning of the pulse and red-shift at the tail.
t
S
Svg
t
gSvN
qd
J
t
N
RSgvt
S
eeg
g
spg
d
d
2 )(
2
1
d
d
d
d
)(d
d
Neglect spontaneous emission
chirp
Disadvantage of Chirp
In optical fiber, the dispersion is anomalous at 1.55m, that is the refractive index of the fiber increases with the wavelength (cf. the regular optical glass in which the blue color sees higher refractive index than the red color does). The speed of light in a media is inverse proportional to the refractive index. Therefore in fiber, longer wavelength signal travels in a slower speed.
Because of this property, if the signal pulse has higher frequency component (positive chirp) at the beginning of the pulse and lower frequency (negative chirp) at the tail, the pulse will experience broadening during the propagation through the fiber.
If we can have opposite frequency chirp for the signal pulse, we can compensate the dispersion in the fiber.
Chirp generated in EA Modulator
We have loss instead of gain. With the increase of carrier density inside the active region due to the absorption of signal, we expect the linewidth enhancement factor to become negative if generated carriers can be removed from the cavity fast enough (i.e. no saturation of absorption or equivalently speaking no gain is generated by the small amount of inversion carriers).
This will give us positive chirp when signal has negative slope (corresponding to the tail).
Aside: In semiconductor optical amplifier, the chirp has opposite sign to laser diodes if gain is saturated.
t
S
Se
d
d
2
N
Nn
Ng
Nn
dNn
dNne d/)(d
d/d4
d/d
d/d4
/]dIm[
/]dRe[ ''
Negative Chirp in EA Modulator
Chirp is a strong function of wavelength and bias
40GHz Modulator
Long device length for easy packagingshort modulator section to reduce capacitancelow driving voltage of 2.8Vlow insertion loss of 8dB
Integrated DFB-EA Transmitter
10Gb/s module, Ith = 20mA, Pmax = 4mW @80mA , extinction ratio = 15dB for -2.5V.