E - 東北大学電気・情報系

Download Report

Transcript E - 東北大学電気・情報系

Communications Devices
Lecture on 6/24, 7,1, 7/8, 7/15, 7/22, 7/29, 2014
Hirohito YAMADA
About lecture
1. Schedule
6/24 Basic of semiconductor photonic devices
7/1 Matter-electromagnetic wave interaction based on semi-classical theory
7/8 Electromagnetic field quantization and quantum theory
7/15 Optical transition in semiconductor, Photo diode, Laser diode
7/22 Optical amplifier, Optical modulator, Optical switch,
Optical wavelength filter, and Optical multiplexer/demultiplexer
7/29 Summary
2. Textbook written in Japanese
米津 宏雄 著、光通信素子工学 - 発光・受光素子 -、工学図書
霜田 光一 編著、量子エレクトロニクス、裳華房
山田 実著、電子・情報工学講座15 光通信工学、培風館
伊藤弘昌 編著、フォトニクス基礎、朝倉書店 第5章
3. Questions E-mail: [email protected], ECEI 2nd Bld. Room 203
4. Lecture note dounload URL: http://www5a.biglobe.ne.jp/~babe
Data traffic explosion on the Internet
Trend of data traffic processed by a domestic network node
Data rate (Gbit/s)
Rapid increase of double per 3 years
Double/3 years
double/year
M/D/Y
http://www.jpix.ad.jp/en/technical/traffic.html
Growth of internet traffic in Japan
Total download traffic in Japan was about 2T bps at the end of 2012
Daily average value
Annual growth rate: 40%
Total download traffic in Japan
Total upload traffic in Japan
2003
2004
2005
2006
2007
year
2008
2009
2010
2011
2012
Cited from: H25年度版情報通信白書
Optical fiber submarine networks
Cited from http://www1.alcatel-lucent.com/submarine/refs/index.htm
Power consumption forecast of network equipments
Total internet traffic (Tbps)
Annual power consumption of
network equipments (×1011 Wh)
Domestic internet traffic is increasing 40%/year
If increasing trend continue, by 2024, power consumption of
ICT equipments will exceed total power generation at 2007
Total power generation
at 2007
year
http://www.aist-victories.org/jp/about/outline.html
Expanding applied area of optical communication
Nowadays, application area of optical communications are spreading
from rack-to-rack of server to universal-bass-interface of PCs
Backplane of a server
(Orange color cables are optical fibers)
Universal Bass interface (Light Peak)
installed in SONY VAIO Z
Photonic devices used for optical communications
Various optical devices for use in optical networks
1. Passive optical device, Passive photonic device
- Optical waveguide, optical fiber
- Optical splitter
: Treated in our lecture
- Optical directional coupler
- Optical wavelength filter
- Wavelength multiplexer/demultiplexer (MUX/DEMUX)
- Light polarizer
- Wave plate
- Dispersion control device
- Optical attenuator
- Optical isolator
- Optical circulator
- Optical switch, Photonic switch
- Photo detector, Photo diode (PD)
2. Active optical device, Active photonic device
- Light-emitting diode (LED)
- Semiconductor laser, Laser diode (LD)
- Optical amplifier
3. Other devices(Wavelength converter, Optical coherent receiver, etc.)
Photonic devices supporting life of 21st century
Various photonic devices used for applications other than optical communication
- charge-coupled device (CCD) image sensor
- CMOS image sensor
- solar cell, photovoltaic cell
- photo-multiplier
- image pick-up tube
- CRT: cathode-ray tube, Braun tube
- liquid crystal display (LCD)
- plasma display
- organic light emitting display
- various recording materials (CD, DVD, BLD, hologram, film, bar-code)
- various lasers (gas laser, solid laser, liquid laser)
- non-linear optical devices
These devices collectively means “Photonic devices”
What is Photonic Device
Manipulating electron charge
Electronics
Magnetic Recording
voltage and current
Manipulating
Electron Tube, Diode, Tr, FET, LSI wavefunction of electron
Tunnel effect devices
Manipulating
both spin and Spintronics
charge of
GMR HDD
electrons
MRAM
Photo-electronics
(Opt-electronics)
Display
CCD, CMOS sensor Manipulating
photon and
solar cell
electron charge
LD, LED
Unexplored Photo detector, PD Manipulating
photon
Magnetics
Photonics
Opto-spinics
Energy and
number
Laser
magnet-optical
HDD
Optical disk
Electromagnetic
disk
Magnetic tape
wave
Optical fiber
Manipulating spin
amplitude and
of materials
phase
?
Basis of semiconductor photonic devices
Properties required for semiconductor used as photonic devices
105
・ Transparent at operating wavelength
・ Small nonlinear optical effect
(as distinct from nonlinear optical devices)
・ Better for small material dispersion
・ Small birefringence
(polarization independence)
Active devices(Light-emitting devices)
・ Moderately-opaque at operating
wavelength
・ High radiant transition probability
in case of light-emitting devices
(Direct transition semiconductor )
・ To be obtained pn-junction
(To be realized current injection
devices)
Optical absorption coefficient: α (cm-1)
Passive devices(Non light-emitting devices)
T = 300K
Si
104
In0.53Ga0.47As
GaAs
InP
Ge
103
102
0.4 0.6 0.8 1.0 1.2 1.4 1.6
Wavelength: λ (μm)
Optical absorption coefficient of
major semiconductor materials
Band structure of semiconductors
Semiconductor and band structure
According to Bloch theorem, wave function of an electrons in crystal
is described as a quantum number called “wave number” This
predicts existing dispersion relation between energy and wave
number of electron. This relation is called energy band (structure)
Electron
Hole
In bulk Si, Holes distribute at around
the Γ point, on the other hand,
electrons distribute at around the X
point (Indirect transition
semiconductors)
Band gap ~1.1eV
Dispersion relation of electron energy in Si
Band structure of semiconductors
Band structure of compound semiconductors
Both electrons and holes distribute at around the Γ point
(Direct transition semiconductors)
Electron
Electron
Hole
Hole
GaAs
InP
Band structure of semiconductors
Band structure of Ge
Ge which is group Ⅳ semiconductor is also
indirect transition type semiconductor, but
by adding tensile strain, it changes a direct
transition-like band structure
Recent year, Ge laser diode(RT, Pulse)
was realized by current injection
Band structure of Ge
Conduction
band
Conduction
band
Indirect transition
Direct transition
1.6%
tensile strain
Valence band
Valence band
Basis of semiconductor photonic devices
Material dispersion
Values of dielectric constant (refractive index), magnetic permeability depend
on frequency of electromagnetic wave interacting with the material
n ( ),
 ( )
Dielectric constant (refractive index) significantly
changes at the resonance frequency of materials.
In linear response, between real part and
imaginary part of frequency response function
holds Kramers-Kronig relation.
Real part
 ( ),
Photon energy (eV)
Phenomenologically-derived equation of relation
between wavelength and refractive index
n  1 
2
Ai 
2
  i
2
2
Here, λi = c/νi, c: light speed, νi: resonance
frequency of material, Ai: Constant
Imaginary part
W. Sellmeier equation
Photon energy (eV)
Calculation of dielectric function ε of Si
Basis of semiconductor photonic devices
Birefringence
In anisotropic medium, dielectric constant (refractive index), magnetic permeability
are tensor
 Dx

 Dy
D
 z
   xx
 
    yx
 
  zx
 xy
 xz   E x 
 yy
 yz
 zy
 zz

 E y

 E z




(Material equation)
Incident ray
Optical axis
Outgoing ray
extraordinary
ray
Crystal is optically anisotropic medium
When light beam enter a crystal, it splits two
beams (ordinary ray and extraordinary ray)
ordinary ray
In birefringence crystals, an incident direction
where light beam does not split is called optical
axis (correspond to c-axis of the crystal)
Birefringence of calcite
Basis of semiconductor photonic devices
Nonlinear optical effect
Values of dielectric constant (refractive index), magnetic permeability depend
on amplitude of electromagnetic field interacting with the material
 ( E ),
n ( E ),
 (H )
When strong electric field (light) applied to a material, nonlinear optical effects
emerge. Wavelength conversion devices use this effects.
When intensity of incident light is weak, linear polarization P which proportional
to the electric field E is induced.
D ( x , t )   E ( x , t )   0 r E ( x , t )
Linear polarization
  0 (1   e ) E ( x , t )   0 E ( x , t )   0  e E ( x , t )
  0 E ( x,t)  P ( x,t)
P  0 E
 : Electric susceptibility
When intensity of incident light become strong, electric susceptibility become
depended on the electric field E
P   0 (E )E
 ( E )   0 
(1 )

(2)
E
(3)
EE  
(4)
EEE    
Basis of semiconductor photonic devices
Electron transition in direct transition type
and indirect transition type semiconductors
In case of indirect transition, phonon
intervenes light emission or absorption
Band structure of (a) GaAs and (b) Si
Light emission from materials
Light absorption and emission in material
Excited state
Nucleus
All light came from atoms
Sun light‥‥Nuclear fusion of hydrogen
Fluorescence of fireflies ‥‥Chemical reaction of
organic materials
Light
Light from burning materialsの‥‥ Chemical reaction
of organic materials
Ground state
Excited state
ν
Electroluminescence from LED‥‥Electron transition ΔE
in semiconductors
Light
Ground state
ΔE = hν
Why materials emit light?
You have to learn about interaction mechanism of matter with
electromagnetic field if you want to understand this
In the field of Quantum electronics
Light emission from materials
According to classical electromagnetics, Rutherford
atom model is unstable. It predicts lifetime of atoms
are order of 10-11 S. (See the final subject in my
lecture note ElectromagneticsⅡ)
v
ω
Electron -e
m
Proton
r
+e
In order to solve this antinomy, Quantum mechanics
was proposed
N. Bohr proposed an atom model which electrons
exist as standing wave of matter-wave. The shape
of the standing wave is defined by the quantum
condition, and it is arrowed in several discreet
states. When an electron transits from one steady
state to other state, it emit / absorb photon which
energy correspond to energy difference between
the two states. (ΔE = hν)
Why photon (light) is emitted / absorbed
when electron transition between the states
Rutherford atom model
Proton
-e
+e
Bohr hydrogen atom model
Theory which describe light absorption and emission
There are three method describing interaction of matter with electromagnetic fields
1. Classical theory
2. Semi-classical theory
Energy of electrons in semiconductor is quantized (Band structure),
On the other hand, energy of electromagnetic field is treated by
classical electromagnetics
3. Quantum theory
Energy of electrons in semiconductor is quantized (Band structure),
Electromagnetic field is also quantized (Field quantization)
Electromagnetic Optical
field
absorption
Stimulated Spontaneous
emission
emission
Method
Matter
Classical theory
Classical
Classical
Possible
Impossible
Impossible
Semi-classical theory Quantum
Classical
Possible
Possible
Impossible
Quantum
Possible
Possible
Possible
Quantum theory
Quantum
Description of electromagnetic field
In order to understand interaction of matter with electromagnetic field, we need to
understand description of electromagnetic field and its fundamental characteristics
Maxwell equations
  E ( x,t)  
B ( x , t )
  B ( x , t )  
  E ( x,t) 
t
E ( x , t )
t
 ie 
1 E ( x , t )
v
2
t
  ie
e

  B ( x,t)  0
Electric field E and magnetic field B can be also described as follows with
electromagnetic potential A(x, t) and ϕ(x, t)
E 
A
t
  ,
B  A
Therefore, Maxwell equations can be replaced to equations with A and ϕ as
values of electromagnetic field, instead of E and B
Description of electromagnetic field
Electromagnetic potentials can be described as follows with arbitrary scalar
function χ(x, t).
A L ( x , t )  A ( x , t )  grad  ( x , t ),
L ( x,t)   ( x,t) 
 ( x , t )
t
This function χ(x, t) is called a “gauge function”, and selecting these new
electromagnetic potential AL and ϕL is called “gauge transformation”.
When χ(x, t) was selected as AL and ϕL satisfying the following relation,
div A L ( x , t )  
 L ( x , t )
t
0
the gauge is called “Lorenz gauge”. In this case, basic equations that describe
electromagnetic phenomena is reduced to two simple equations regarding AL
and ϕL as follows.
2


   
2

t


 A L ( x , t )    i e ( x , t ),


2


   
2

t


1
 L ( x , t )    e ( x , t )



Description of electromagnetic field
Other than Lorenz gauge, when A was selected as satisfying   A  0 condition,
the gauge is called “Coulomb gauge”.
In this case, basic equations that describe electromagnetic phenomena is as follows
  


2


    2

t



 A  
   i


t

In free space where both electric charge ρe and electric current ie do not exist,
  0,
2


   
2

t


A  0


Therefore, E and B is derived from A with the following relations
E 
A
t
B  A
Interaction of charged particles with electromagnetic fields
Hamiltonian of the particle when single charged particle (electron) is in
electromagnetic field is described as follows
H 
1
( p  e A )  eV
2
2m
Here, p is the momentum operator, m is electron mass, V is potential of
electron, e is elementary charge, and A is vector potential. Here, H0 is an
Unperturbed Hamiltonian which is for an electron in space without
electromagnetic field, and Hint is an Interaction Hamiltonian which is
originated from interaction between an electron and electromagnetic field.
H 
1
2m
H0 
( p  e A )  eV  H 0  H int
2
1
p  eV
2
2m
H int  
e
2m
( p  A  A  p) 
e
2
2m
A
2
となる。
The last term in Hint is proportional to A2, and it reveals higher order effects
(nonlinear optical effect). Here, we ignore it because the contribution is small.
Electrical-dipole approximation of the interaction
Position of a charged particle is described as
r  r0 e
i t
 r0 e
*
 i t
Momentum of a charged particle is described as
p  m r  m ( i  r0 e
i t
 i  r0 e
*
 i t
)
Vector potential is described as
A  A0 e
i t
 A0 e
*
Polarization
r of electron
cloud
E
Polarization of atom
 i t
+e
e+iωt
E
Electric field is
E 
A
t
  i A0 e
i t
 i A0 e
*
r
 i t
e−iωt
−e
Therefore,
Electrical dipole
H int  
e
( p  A  A  p)  
2m
e
m
m ( i  r0  A 0  i  r0  A 0 )
*
 e r  E
 R  E
Here, R = er, and it called electric dipole moment
*
ei2ωt や e−i2ωt の項は、
時間積分すると消える
Electrical-dipole approximation of the interaction
Interaction Hamiltonian
H int  
e
( p  A  A  p)
2m
には、電場による影響(クーロン力)と磁場による影響(ローレンツ力)の両方が含まれ
ているはずであるが、上で導いた荷電粒子との相互作用の式では、電場による影響
RE しか含まれていない。これは、荷電粒子の運動の形態を、 r  r0 e i t  r0* e  i t
のように、ある限られた場所での振動と仮定したためである。もし荷電粒子に平行運
動のような運動形態を仮定したならば、磁場による影響も含まれてくるはずである。
このように相互作用を、電場のみに影響すると仮定して RE のように近似することを、
電気双極子近似という。近似による変位の部分を電場で展開すると、電気四重極子
や電気八重極子のような多重極の分極が現れることもある。
工学で扱う物理現象は非常に複雑なものが多い。従って、全ての物理現象を取り入
れた完璧な理論を構築することは不可能である。
良きエンジニアとは、それら複雑な物理現象の中で、何が本質的に重要かを見極め、
近似をうまく使い、無視できる物理現象は思い切って無視し、シンプルな理論を構築
できる人である。ただし、どんな近似を使ったのかは決して忘れてはいけない。
Electrical dipole in semiconductor
In semiconductor, an electron and hole pair forms a electrical dipole
Conduction
band
−e
Electron
+
Electron
Electrical dipole
E
Electrical dipole
+e
Hole
Hole
Valence band
Electrical dipole in semiconductor
Quantum statistics and density matrix
電流値など、多くの粒子(電子)が関与する物理量は、多くの粒子による統計的な
平均値となる。また単一事象において、複数回観測(測定)を行った場合の期待値
なども統計的な扱いが必要となる。そのような多数の粒子或いは多数回の測定に
ついての統計的性質について扱う。まず、ν 番目の粒子または ν 回目の測定にお
ける状態を、

( )

C
( )
n
と書くこととする。
n
n
ここで n は、一つの粒子のみが存在する場合のエネルギー固有状態である。完備
性が仮定されているので、 n の1次結合によって任意の  ( ) をつくれるはずであ
り、 n には添え字 (ν) はいらない。
If operator for a physical quantity is assumed as A , ν 番目の粒子における期待値は、

( )
A
( )

C
( )*
n
( )
Cm
n A m
となる。
n ,m
次に、粒子の集団全体における期待値の平均(ensemble average)を求める。
ν 番目の粒子が寄与する割合(粒子の抽出確率)を P(ν) とおき、規格化しておく。


P
( )
1
Quantum statistics and density matrix
Statistical average of expected value is described as
 A 

P
( )

( )
A
( )



P
( )
( )*
Cn
( )
Cm
n Am
 ,n ,m
ここで、以下のような書き換えを行う。
A   A ,
A nm  n A m ,
 mn 
P
( )
( )*
Cn
( )
Cm
 C nC m
*

  11

 21

 



ρmn を要素にもつ行列 ρ をis density matrix
密度行列を用いると、
A 

 mn A nm 
n ,m


m
n ,m


m   n n  A m 
 n

I 


m  n n A m

n n

m A m
 12

 22












と書き改められる。
m
は恒等演算子(identity operator)
n

m A m
は行列 ρA の対角要素を全て足し合わせたもの、即ちTrace であり、
m
A 

m
m  A m  Tr (  A )
と表わすことができる。また、 Tr (  )  1 である。
密度行列の運動方程式
( )
密度行列 ρ の性質が分かれば、集団内の個々の粒子についての状態 
や抽出
確率 P(ν) を知らなくても、統計性を含めた期待値 A を知ることができる。そこで、密度
行列 ρ を表現する方程式、つまり ρ が従うべき方程式を求めてみる。
まず、密度行列の行列要素の定義式を、以下のように書き直す。
 mn 
P
( )
( )*
Cn
( )
Cm




P
( )

( )
m 
n
( )

m 
( )
P
( )

( )

従って密度行列は、  



n  m  
 
( )
P
( )

( )
( )
P
( )

( )

 n

と書くことができる。

この式の両辺を時間 t で微分すると、
d
dt



 d  ( )
( )
P


 dt
( )
 
( )
P
( )
d 
dt
( )



ただし、抽出確率 P(ν) の時間依存性はないとしている。
となる。
密度行列の運動方程式
Schroedinger equation i 
 i
d

( )
 
( )
d

( )
 H 
( )
dt
and its Hermitian conjugate
より、
H
dt
d

dt



1
i
H


 i
( )
P
( )
H   H  

1
i
( )
 
H ,  
( )
P
( )

( )
H 

i 
となる。
これは、密度行列の時間発展を表す式であり、密度行列の運動方程式或いは、
量子リウヴィル(Liouville)方程式もしくはリウヴィル-フォン・ノイマン方程式とも呼
ばれる。
双極子との相互作用がある場合の密度行列
電子系の主ハミルトニアンを H0 とし、電気双極子能率を R とする。前に述べたように、
電場が存在するときの相互作用ハミルトニアン Hint は、−RE となり、電子系全体のハ
ミルトニアン H は、
H  H 0  H int  H 0  RE
となる。
これを、密度行列の運動方程式に代入すると、
d

dt

1
i
1
i
H   H  
1
i
H 0   H 0  
H 0   H 0  
1
i
  R  R  E
1
i
  RE
 RE  
となる。
ここで、最後の ≈ では、1個の電子が存在する領域が電磁波の波長に比べて十分に
小さく、その範囲内で電場の分布は一定と見なせることを仮定している。実際、気体
原子に束縛されている電子の存在範囲はせいぜい数Å程度であり、また半導体中の
電子に至ってもせいぜい数十Å程度である。それに対して、相互作用する光の波長
は数千Åもあるので、この仮定は妥当である。
双極子との相互作用がある場合の密度行列
次に、密度行列の行列要素に対する方程式を導出する。エネルギー固有状態には
時間依存性が無いので、
d  mn

dt
d  mn
d
m  n  m
d
dt

dt
1
i
と表すことができ、従って、
n
dt
m H 0 n 
1
i
m H 0 n 
1
i
m R n E 
1
i

1




m H 0  l l  n 
m    l l H 0 n
i
i
 l

 l


1




m    l l R n E 
m R  l l  n E
i
i
 l

 l


1
1
1

i
m H0 l l  n 
1
m  l l R n E
l

m R n E
i

l
1

i
m  l l H0 n
l
1
i

m R l l  n E
l
と書くことができる。なお。ここで用いている固有状態は、主ハミルトニアン H0 の固有
状態である。
双極子との相互作用がある場合の密度行列
従って、状態 n は、エネルギー固有値 Wn をとり、状態間には直交性があるので、
d  mn

dt

1
i
1
i
Wm m  n 
1
i
 mn W m  W n  
 i  mn  nm 
1
i
m  n Wn 
 
1
i
ml

1
i
 ml R ln E 
l
R ln  R ml  ln  E

m  l l R n E
l
1
i
R
ml
1
i

m R l l  n E
l
 ln E
l
となる。
l
ただし ωmn は、順位 m と順位 n とのエネルギー差に対応する角周波数であり、
W n  W m    nm
で与えられる。
2準位系での相互作用
話を簡単にするためにまず、電子の取り得るエネルギー準位が2つしかない2準位系
において、電磁場との相互作用を考える。
b
Wb
上側(励起)準位を b 、下側(基底)準位を a とすると、
d  bb
dt
 i  bb  bb 
1
i
  ba R ab
  bb R bb  R ba  ab  R bb  bb  E
  ba
a
Wa
となる。ここで、  bb  W b  W b  /   0 であり、双極子能率 R の行列要素の対角要素
Rbb も、一般的な媒質では 0 となる。従って、
d  bb

dt
d  aa
dt
d  ba
dt
d  ab
dt

1
  ba R ab
 R ba  ab  E
  ab R ba
 R ab  ba  E
i
1
i
 i  ba  ab 
 i  ab  ba 
が得られる。
1
i
1
i
となる。同様に、
  ba R aa
  bb R ba  R ba  aa  R bb  ba  E   i  ba  ba 
  aa R ab
  ab R bb  R aa  ab  R ab  bb  E  i  ab  ba 
1
i
1
i
  bb   aa R ba E
  aa
  bb  R ab E
2準位系での相互作用
ここで  bb ,  aa は各々、上側および下側準位の電子分布であり、 Tr (  )   aa   bb  1
である。一方  ab ,  ba は2つの状態間での量子相関(量子コヒーレンス)を表している。
これらの式は、それらの時間的変化を記述しており、第1式は上側準位の電子分布
の時間的変化を、第2式は下側準位の電子分布の時間的変化を、そして第3式と4
式は、それら準位間での量子コヒーレンスの時間的変化を記述している。
これらの式は、光の吸収(absorption)や誘導放出(stimulated emission)を記述して
いるが、実際にはこれら方程式に含まれていない以下の現象が存在する。
上側準位の電子は誘導放出により下側準位に遷移し、やがて熱平衡状態での電子
分布になったら正味の発光は起きなくなる。そこで、連続して発光させるためには、
下側準位の電子を上側準位に汲み上げてやる必要がある。これをポンピングと言
い、半導体レーザーや発光ダイオードなどでは、pn接合に電流を注入することによ
りポンピングを行っている。
上側準位の電子は、電磁場が存在しなくてもある一定の割合で光を放出して下側準
位に遷移する。これを自然放出と言う。自然放出は電磁場を量子化した時の不確定
性に起因するもので、電磁場を古典論で扱う場合(半古典論)には導出できない。
これまでの密度行列の導出においては、粒子同士の衝突やエネルギー準位の揺ら
ぎなどは考慮していなかった。しかし実際は、粒子の衝突などによって双極子の調和
振動が乱され、双極子振動は減衰していく。これを電子緩和効果という。
2準位系での相互作用
このような効果を現象論的ではあるが付加してやる必要がある。
ポンピングによって上側準位の電子密度が増加(その分下側準位は減少するが、)
する割合を Λp、自然放出による電子寿命時間(縦緩和時間とも言う)を τs、粒子同士
の衝突などにより双極子振動が減衰していくが、その寿命時間(デコヒーレンス時
間)(横緩和時間とも言う)を τd とする。ポンピングと自然放出は、密度行列の対角要
素に付加され、双極子の寿命時間は非対角要素の変化として導入される。従って、
これら現象論的な補正をした方程式として、
d  bb

1
dt
i
d  aa
1
dt
d  ba
dt
d  ab
dt

i
  ba R ab
 R ba  ab  E 
  ab R ba
 R ab  ba  E 
  i  ba  ba 
 i  ab  ba 
が得られる。
1
i
1
i
  bb
  aa
 bb
s
 bb
s
 Λp
 Λp
  aa  R ba E 
  bb  R ab E 
 ba
d
 ab
d
 (1)
 (2)
Wb
自然放出
τs
  ba
 (3)
 (4)
ρbb
Wa
ρaa
b
ポンピング
Λp
a
2準位系での分極率
次に、微分方程式
d  ab
 i  ab  ba 
dt
1
i
  aa
  bb  R ab E 
 ab
d
 ( 4 ) の解を求めてみる。
電場が存在しない場合、上の方程式は同次方程式
d  ab
dt
 i  ab  ba 
 ab
d

1 
  ab
  i  ba 
 d 


1 
となり、その解は  ab  C exp  i  ba   t
d 

の形となる。ただし C は任意定数。そこで、電場が存在する場合については U(t) を
変数として、

 ab  U ( t ) exp  i  ba 

1 
t
 d 
と表わせると考え、これを式(4)に代入すると、



1 
1 
1 
 t   i  ba 
U exp  i  ba 
t
exp  i  ba 




dt
d  
d 
 d 


dU

1 
1
  ab 
  aa   bb  R ab E
  i  ba 


i

d 


1 

U ( t ) exp
  i  ba 
 d 


1 
1
 i  ba 
t 
  aa   bb  R ab E


d 
i

となる。
2準位系での分極率

1 
1
t 
  aa   bb R ab E
exp  i  ba 

dt
d 
i

dU
従って、
の関係が得られる。
電場 E を E  E  e i t  E * e  i t と表し、上式を時間に関して積分すると、
U 

1
i
1
i
0   aa
t
 
t
0
aa

1  
  Ed 
  bb  R ab exp    i  ba 
 d  

  bb  R ab


1  

  E  e i  t  E * e  i  t d 
exp    i  ba 


d

 

1 
t
 ab  U exp  i  ba 
 d 



  aa
  bb  R ab
i
  aa
  bb  R ab
i

となる。
積分時間内 0 ~ t で (ρaa – ρbb)Rab は一定とみなすと、

1  t
 t exp
exp  i  ba 
 d  0


1  

  E  e i   E * e  i  d 

i



ba
 d  


1 

 i    b a 


1   t
d 

t  E  e
exp  i  ba 
d 
 0

d  


次のスライドに続く、、、、




t
0
*
E e

1 
  i    b a 

d 



d 

2準位系での分極率



  aa
  bb  R ab
i
  aa
  bb  R ab
i
  aa
  bb  R ab
i
e
e

1
 i ba 

d


1
 i ba 

d


 

1 
t 
 i    b a 

 
E
d 



e

1
  i    ba  
d
 


t




1

 i b a 


E
 i t
d


e  e
 i      1 

ba

d
  aa
  bb  R ab
i



1 
*


  i    b a 

E
d 

 
e
1

  i     
ba


d
0











0
t




1 
1 
*






 i    b a 
t
  i    b a 
t
E
E
d 
d 


e
e
 1 
 1 

1 
1
  i     


 i    ba  




ba


d
d
ここで、t >> τd の定常状態を考えると e
 ab 
t


t


*


E

  i t

e

e
1
  i    ba  

d
t
d
0
であるから、

i t
*  i t

E e
E e


1
1
 i    ba  
 i    ba  

d
d





となる。

1
 i ba 

d


t






 

2準位系での分極率
さらに ω ≈ ωbaであるから、{ }内の第1項の分母の虚部はほぼ 0 となるが、第2項の分
母の虚部はかなり大きい。従って、 { }内の第2項は無視できる。そうすると、
 ab  
  bb
  aa  R ab
i
E e
i t
i    ba  
1
となる。
d
同様にして ρba も求められるが、ρba は ρab の複素共役であり、 ρba = ρ*ab である。
 ba 
  bb
  aa  R ba
i
*
E e
 i t
i  ba    
となる。
1
d
ところで、古典電磁気学において分極 P は、
P   0  ( ) E  e
従って、  ( ) 
j t
  ( ) E  e
1
 0 E  Δt
*
*

t
Pe
t  Δt
 j t
 j t
dt ,

で与えられた。
 ( ) 
*
1
 0 E  Δt
*

t
Pe
t  Δt
ただし、Δt は 1/ω の数倍程度の時間間隔をとることとする。
j t
dt
2準位系での分極率
ところで、双極子能率 R は、R = er と表わせたが、これを量子力学的な演算子とし、
電子密度を Nt とすると分極 P は、
P  N t R  N t Tr   R   N t   mn R nm
となる。
m ,n
従って、
 ( ) 
Nt
 0 E  Δt


t
t  Δt
  aa R aa
N t   bb   aa 
i  0
  ab R ba   bb R bb   ba R ab e
R ab
 j t
dt
2
i    ba  
1
となる。
d
ただしここでは、前にも述べたように Raa = Rbb = 0 と見なしている。
光増幅利得
分極率 χ の虚部は、それが正の場合は光を増幅する割合を示す光学利得定数を、
負の場合は光減衰係数を表わしている。つまり、光学利得定数 g は、
g 



 0  Im  



N t   bb   aa  R ab

 0  r  N t   bb   aa  R ab
 0 r

磁性体でなければ、  r  1,
2


i
Im 
1
 i    ba  

d





1
2
d

  ba 
2
 1
 
d
 r  nr




2

 0  N t R ab
0
nr 
  bb
2

  aa 
  ba 
2
1
d
 1
 
d




となる。
準位 b と a に存在する電子数の合計は単位体積当たり Nt である。準位 b の電子密
度を Nb、準位 a の電子密度を Na とする。ρbb と ρaa とは電子分布の比率を示すので、
N t  bb  N b ,
g 
N t  aa  N a ,
 0  R ab
0
nr 
2
Nb  Na  Nt
 N b  N a 1  d 
   ba 2  1  d 2
2
と書かれる。従って光学利得定数は、
となる。
反転分布
ポンピングを行わない熱平衡 (thermal equilibrium)状態では、電子分布はMaxwellBoltzmann分布に従い、Nb < Na であり、従って光学利得定数 g の値は負となる。つ
まり、媒質によって光が吸収される。これを光物性学では、基礎吸収と呼んでいる。
一方、ポンピングによって Nb > Na が実現されると、利得定数が正となり、光が増幅
されて出てくる。これが誘導放出である。 Nb > Na の状態を反転分布(population
inversion)という。また、Nb = Na の時、媒質は透明となる。
エネルギー E
E
b
Nb Wb
b
a
Na
Wa
a
粒子数 N
熱平衡状態 (Nb < Na)
Maxwell-Boltzmann分布
P(E )  e

E
kT
k: ボルツマン定数
T: 媒質の温度
Nb
Wb
Na Wa
粒子数 N
反転分布 (Nb > Na)
光増幅のメカニズム
Re χ
分極 P に光電場 Ei が入射する場合、
P   0  E i   0  Re   i Im  E i
と表わせ、Im χ が負の場合、分極
は入射光 Ei より 0 ~ π/2 位相が遅
れて振動する。この分極によって
放射される光電場 Er は、 Δz を分
極が存在する範囲として、
Er  
i

2

Re χ
ωba
ω
Im χ
Im χ
 P Δz
一方、Im χ が正の場合、分極は入
射光 Ei より 0 ~ π/2 位相が進んで
振動する。この分極によって放射さ
れる光電場 Er は分極振動に対し
てπ/2 位相が遅れて放射される。
従って、出力光電場は増幅される。
ω
ωba
と表わされ、分極振動に対してπ/2
位相が遅れて放射される。従って、
出力光電場は Ei とEr との合成とな
るため、出射光は減衰する。(右図)
ω
ωba
ρbb < ρaa
ρbb > ρaa
0 ~ π/2
Ei
Er
P
0 ~ π/2
Eout
P
減衰
ω
ωba
Ei
Er
増幅
Eout
レーザー
レーザとは、光の発振器
正帰還回路
光の正帰還回路
光増幅媒体
+ Amp.
鏡
電気の発振器
レーザー
光増幅媒体とはどのようなものか?
物質(原子系)と光との相互作用 以下の3つの課程が同時に起きている
電子など
E2
減衰
入射光
出射光
増幅
入射光
発光
出射光
E1
二準位系
(原子など)
光の吸収
誘導放出
自然放出
熱平衡状態
Maxwell-Boltzmann分布
E
E2
k: ボルツマン定数
T: 媒質の温度
P(E )  e

n2: 励起状態の原子数
E
kT
誘導放出
正味では減衰
E1
吸収
P(E)
熱平衡状態では、励起準位の原子数
は基底準位の原子数よりも少ない n1> n2
吸収
吸収
n1: 基底状態の原子数
自然放出の起きる確率 = An2
A: アインシュタインのA係数
吸収の起きる確率 = Bn1 I
B: アインシュタインのB係数
誘導放出の起きる確率 = Bn2 I
I: 入射光の強度
Bn1 I > Bn2 I
熱平衡状態では、吸収の確率 > 誘導放出の確率となり、入射光は減衰して出てくる
反転分布
反転分布
Tが負(負温度状態)
E
n2: 励起状態の原子数
E2
P(E )  e

E
kT
誘導放出 誘導放出
E1
P(E)
励起準位の原子数が基底準位の原
子数よりも多い状態を反転分布という
n1< n2
誘導放出
吸収
正味では増幅
n1: 基底状態の原子数
Bn1 I < Bn2 I
反転分布では、誘導放出の確率 > 吸収の確率となり、入射光は増幅されて出てくる
レーザーとは、何らかの方法で反転分布を作り出し、放射の誘導放出(Stimulated
emission)を用いて光を増幅する装置
レーザー
レーザー(laser)の語源
light amplification by stimulated emission of radiation
光増幅
誘導放出
放射
メーザー(maser) microwave amplification by stimulated emission of radiation
レーザーの種類
A. 固体レーザー
a. 結晶体レーザー
ルビー(Cr3+: Al2O3)レーザー 1960年世界で最初に発振したレーザー。
光励起で発振。λ = 694.3nm(赤)
Nd3+: YAGレーザー λ = 1064nm(近赤外)
b. ガラス・レーザー 光学ガラスにNdを
ドープ。高出力
レーザー
B. 液体レーザー
有機色素(ローダミンなど)を光励起して発振。広範囲に波長可変可能。
C. 気体レーザー
a. 中性気体レーザー
He-Neレーザー He-Neの混合ガスで
グロー放電を起こして発振。
λ = 632.8nm(赤)
b. イオンレーザー
Ar+レーザー, Kr+レーザー 中性原子のガスでアーク放電を起こして発振。
Ar+レーザー λ = 488nm(青), 514.5nm(緑)
Kr+レーザー λ = 647nm.1, 676.4nm(赤),
c. 分子レーザー
CO2レーザー 大出力のため、金属加工用や医療用レーザーメスなどに使用
λ = 10.6μm(赤外)
D. 半導体レーザー コンパクトで低消費電力のため、光通信用途を始め各種
民生機器で広く使用。
レーザーとコヒーレント光
レーザー光は、コヒーレント(coherent)な光
コヒーレントとは、波の位相が揃った状態。高スペクトル純度、良好な収束性を有する
時間的コヒーレンス 空間的コヒーレンス
光
の
電
界
t
光
の
強
度
f 又は λ
光
の
強
度
光
の
電
界
コヒーレント光
t
f 又は λ
インコヒーレント光
(コヒーレントでない)
自然界に存在する光は全てインコヒーレント光
例: 太陽光、炎から出る光、蛍の光、白熱電球、蛍光灯、LED
コヒーレントな光を人工的に発生させる装置がレーザー
場の量子論基礎
電磁場の量子化
誘電率 ε の物質中の電磁場のベクトルポテンシャル A は、以下のように与えられる。
A
1
Q F ( x, y, z ) x

ここで Q は、時間変化を表わす量子力学的演算子であり、x は電磁場(ベクトル
ポテンシャル)の偏り(偏波)の方向を表わす単位ベクトルである。F(x, y, z) は電磁
場の空間分布を表わす実数関数で、古典的なMaxwell方程式から求められる関
数と同じであり、波動方程式
  
0

2
 F ( x, y, z )  0
を満たす。
ω は、古典的電磁場の角周波数であり、上式においては、演算子 Δ において、
電磁界分布 F(x, y, z) を定める固有値として導入されている。
なお電磁界分布関数 F(x, y, z) は、定義された空間全体 V で以下の様に規格化
しておく。

V
F ( x , y , z ) dV  1
2
場の量子論基礎
扱う物質の運動速度が光速度に比べて遥かに小さい場合、相対論は考慮しなく
てもよい。非相対論的な量子力学では、時間に依存する物理量のみが量子力学
的演算子として扱われる(ハイゼンベルク描像) 。つまり Q は、電磁場の振幅の
時間変動を表しているので演算子であり、 F(x, y, z) は電磁場の空間分布を表わ
すが時間には依存しないので、量子力学的演算子ではなく、古典的なMaxwell方
程式から得られる分布関数と考えてよい。
電場 E と磁場 B は、真電荷が存在しない場合即ちスカラーポテンシャルが全
空間で一定の場合、ベクトルポテンシャル A のみから、
E 
A
t
,
B  A
で与えられる。
演算子 Q の時間微分に対応した新しい演算子 P を定義する。
(量子力学では、古典力学のように物理量を直接的に「微分」できないことに注意。)
P 
Q
t
場の量子論基礎
従って、電場 E と磁場 B は、
E 
A
t
1 Q

B  A 
 t
1

F ( x, y, z ) x  
Q   F ( x , y , z ) x 
1

P F ( x, y, z ) x
と書ける。
ここで、量子力学における要請より、演算子 Q と P との間には以下の不確定性関係
が成立しなければならない。
Q , P   QP
 PQ  i 
ここで、演算子 Q と P とからなる以下の新たな演算子 a と a† を定義する。
a
†
a 
1
2
 Q  iP 
1
2
 Q  iP 
これまでの議論より、演算子 Q と P は古典的には cosωt と ‒ωsinωt に対応し、a と
a† は、オイラーの公式を用いて、ω2 =1 と置いてやれば e‒iωt と ‒ eiωt に対応している
ことが推測できる。
場の量子論基礎
a と a† の間には、[a, a†] = 1 の交換関係が成り立つ。何故なら、
a , a   aa
†
†
1

2
1

2
1
†
a a

2
2
 Q  iP  Q  iP    Q  iP  Q  iP 
Q  i  QP  i  PQ  P   Q  i  QP  i  PQ  P
2
2 i   PQ  QP  
2
i

2
2
2

i
P , Q    i    1

また、n を整数値とし、 n を整数値 n (ただし、n ≥ 0)に対応する固有状態とすると、
a は整数値を一つ下げ、 a† は逆に整数値を一つ上げる演算子である。
a n 
†
a n 
n n 1
n 1 n 1
従って、 a を消滅演算子(annihilation operator)、a† を生成演算子(creation
operator)と呼ぶ。
場の量子論基礎
また、 a† と a の積をとった演算子 N = a†aは、
†
N n a a n a
†
†
n n 1 
n a n 1 
n n n n
n
から分かるように、状態 n に作用すると、整数値 n を固有値として返す演算子で
あり、数演算子(number operator)と呼ばれている。
1
†
aa 
2
1
†
aa 
2
 Q  iP  
 Q  iP  
1
2
1
2
 Q  iP  
 Q  iP   i
2

Q
2

P
より、 Q 
より、

†
2
P  i
a  a 

2
a  a 
従って a と a† を用いると、電場と磁場は、
1
E 
B 

1

P F ( x, y, z ) x  i
Q   F ( x , y , z ) x  

2
a  a F ( x , y , z ) x
†

2 
a  a   F ( x , y , z ) x 
†
と表せる。
†
場の量子論基礎
これらから、電磁場のハミルトニアンを求めると、
H
f

1
2

V

1 2
1
1
 †

2
  E  B  dV     a a       N  

2
2




となり、電磁場エネルギーの期待値が、
H
f
1
1


n    N   n    n   n
2
2


と書かれる。
このように、電磁場のエネルギーは量子化すると離散的となり、電磁場(光)のエネ
ルギー量子のことを光子(photon)と言い、n を光子数(photon number)と呼ぶ。こ
こで、光子が無い、即ち n = 0 であっても、ħω/2 の電磁場エネルギーが残ることと
なる。これは所謂ゼロ点振動エネルギーであり、取り出すことはできないが、実は
後で述べるように、これが光の自然放出にとって非常に重要な役割を果たす。即
ち、自然放出はこのゼロ点振動があるため起こり、ゼロ点振動(真空場の揺らぎ)
がないと自然放出は起こらない。
全量子論による物質と電磁場との相互作用
電磁場中を1個の電子が運動している時、系全体のハミルトニアンは、
H 
1
( p  e A )  eV  H
2
2m
f
と表わせる。ここで m は電子の質量、 e は電子の素電荷、 p は電子の運動量演算
子、A はベクトルポテンシャル、V は電子のポテンシャル(極性は電子に対して定義)、
Hf は電磁場のハミルトニアンである。電磁場のベクトルポテンシャル A は、
B  A 
A

2 

2 
a  a   F ( x , y , z ) x 
†
a  a F ( x , y , z ) x
†
であることから、
と表わせる。
ここでベクトルポテンシャルは、クーロン・ゲージならば   A  0 の関係にあるので、
運動量演算子 p (   i   ) との交換関係は 0 となる。  p , A   p  A  A  p  0
従って、 p  A  A  p としてよい。(ただし、クーローン・ゲージの場合)
補) 系全体のハミルトニアンの導出
ところで、質量 m、電荷 q の荷電粒子が、電磁場 A 内に存在する場合の系全体の
ハミルトニアンは以下のようにして求められる。ただしこの場合、粒子のスピンは考
慮に入れていない。
電磁気学で、荷電粒子に働く力はローレンツ力
F 
dp
 q(E  v  B)
で与えられた。
dt
これは、荷電粒子が U  q (  v  A )

F   U
というポテンシャルの中に存在することを意味する。ここで、ϕ はスカラーポテンシャル。
従って、この系のラグランジュアン L は、
L  T U 
p
2
 q (  v  A )
2m
で与えられる。
これをラグランジュの方程式で、x, ẋ をそれぞれ共役変数とすると、
d  L  L
0


dt   x   x
から系の運動方程式が求められる。
補) 系全体のハミルトニアンの導出
また、 L 
p
2
 q (  v  A )
2m
 L
p
 p  qA
 x
 p
v
より、x に対する共役な変数    は、
 m
である。
従ってハミルトニアンは、
2
p

H  p  x  L  ( p  q A )  v 
 q (  v  A )
2m
 ( p  qA) 

2
p
p
m

p
2
 q ( 
2m
p
 A)
m
 q
2m
従って、電磁場内の荷電粒子のハミルトニアンは、
H 
1
( p  qA)  q
2
2m
となる。
これと、電磁場自身のハミルトニアン Hf を加えて、系全体のハミルトニアン H は
H 
1
2m
( p  qA)  q  H
2
f
となる。
全量子論による物質と電磁場との相互作用
1
系全体のハミルトニアン H 
( p  e A )  eV  H
2
2m
f
を、電子と電磁場が互いに
独立に存在している場合の主ハミルトニアン H0 と、電子と電磁場との相互作用に
基づくハミルトニアン Hint とに分解する。
H0 
1
 †
2
p  eV     a a  
2m
2

1
H int  
e
pA
m
e
2
A
2
2m
主ハミルトニアンの固有状態を u 、これに対応するエネルギー固有値を u とする。
これらは、電子の固有状態  とエネルギー固有値 W、および電磁場の固有状態
n と光子数 n により、以下のように表わせる。
H0 u  u u
1

u  W    n  
2

u  
n
全量子論による物質と電磁場との相互作用
相互作用ハミルトニアンが存在しなければ、状態は初期状態 u から変化しない。
相互作用ハミルトニアンが存在すると、状態が遷移し、その遷移確率は、
2
u f H int u i
に比例する。ここで添字 i は変化前の状態、f は変化後の状態を示す。
これをフェルミの黄金律(黄金則) Fermi’s Golden Ruleと呼び、量子系のあるエネル
ギー固有状態から別のエネルギー固有状態への単位時間あたりの遷移確率を、摂
動法の最低次数の近似によって計算する方法である。
相互作用ハミルトニアン Hint は、電子に対する演算子 p および、電磁場に対する演
算子 a と a† を含んでいるが、p は電子状態  だけに作用し、a と a† は電磁場の状
態 n にだけ作用する。従って、相互作用ハミルトニアン の行列要素は、
u f H int u i  
e

m
2 

F ( x, y, z ) n f a  a
†
n
i

f
pi
となる。
Hint のうちで A2 は、二つ以上の光子が同時に作用する場合にのみ意味を持ち、通
常その確率は小さいので、ここでは無視することにした。(つまり、非線形光学効果
は無視した。)
全量子論による物質と電磁場との相互作用
消滅および生成演算子には以下の性質があったので、
a n 
n n 1
†
a n 
n f a  a
†
n
n 1 n 1
は、以下の二つの場合にだけ値を持つ
i

nf a  a
†
n

 ni  1


 ni
i
( n f  n i  1 のとき )
( n f  n i  1 のとき )
このことから、状態 ui から状態 uf への遷移確率 Pif が、
Pif 
e
2
 m
2

f
pi
2
n  1
2
F ( x, y, z )  
n
(ただし W f  W i    )
(ただし W f  W i    )
となる。このように、電子が上側準位から下側準位に遷移して、光を放出する時の
遷移確率は n + 1 に比例し、逆に光を吸収して、電子が下側準位から上側準位に
遷移する時の遷移確率は n に比例する。
全量子論による物質と電磁場との相互作用
物質が光を放出するか、或いは吸収するかは、統計的な性質で決まる。物質内の
ある二つの電子準位に着目し、上側準位を b、下側準位を a とする。半導体中の電
子などは、フェルミ・ディラックの分布関数に従って分布しており、上側準位での分
布を fb、下側準位での分布を fa とする。電子が上側準位から下側準位へ遷移する
確率は、上側準位に電子が存在し、下側準位に電子が遷移できる空きの状態が存
在する確率に比例するので、fb(1 − fa)(n + 1) に比例する。逆に、下側準位から上側
準位へ遷移する確率は、 fa(1 − fb)n に比例する。両式の差は、(fb − fa)n + fb(1 − fa)
であり、準位 b から準位 a への正味の遷移確率 Pba は、
e 
2
Pba 
f
pi
 m
2
F ( x, y, z )
2
2
( f b 
f a ) n  f b (1  f a )
と書かれる。{ }内の (fb − fa)n に比例する項は、fb > fa のとき誘導放出を表わし、
fb < fa のときは吸収を表わす。 fb(1 − fa) に比例する項は、光子数 n が 0 でも光を
放出し、自然放出と呼ばれる。
b
a
fb
fa
b
n+1
a
光放出
fb
fa
n
光吸収
全量子論による物質と電磁場との相互作用
始状態 ui から終状態 uf への遷移確率は、先に述べたフェルミの黄金律によって与
えられるが、そこには遷移を可能にする状態がなければならない。ここで言う状態と
は、物質のエネルギー状態(基底状態や励起状態)のみならず、物質と相互作用し
ようとする光(フォトン)の状態である。
例えば簡単のために、一辺の長さが L の立方体の箱を考え、その中に光と相互作
用する物質があり、エネルギー状態間遷移によって光を吸収または放出するとする。
ただし L は、相互作用をする光の波長に比べて十分長いものとする。この立方体は
一種の光共振器であり、例えば箱の壁面で電場がゼロとなるような境界条件を満た
さなければならないとすると、この箱の中での電磁界分布は、px, py, pz を正の整数
として、
L
  py 
8
  px 
  pz 
F ( x, y, z ) 
sin 
x  sin 
y  sin 
z
L
3

L
L




L

L
と書かれる。
  
 
2
0

1
 0
2
 F ( x, y, z )  0
k
2
x
 ky  kz
2
2





L
より、各周波数 ω は、



 0 L 
物質
2
p
2
x
 py  pz
2
2

光共振器
を満足しなければならない。正の整数値 px, py, pz をモード番号と言う。
全量子論による物質と電磁場との相互作用
モード番号が作る空間を、図のように書く。個々の座標点が電磁場のモードを示す。
原点を中心とした半径 r の球を描くと、

 



r ,
 0 L 

r 
px  py  pz
2
2
2
であり、この球の体積の1/8が角周波数 0 から ω までのモードの総数である。
従って、角周波数 ω ~ ω + dω の間に含まれる電磁場モードの数は、(球の表面
積)×dr である。なお、電磁場は一つの進行方向に対して直交する 2 つの偏波
状態が存在できる。従って、 dω 当たりのモード数は、
3
pz
3 2

2  4  


 0 L 
1 (  0 )
  2d  
3

8



L  d
3
2
r
となる。上式は、最初に定義した空間の体積 L3 に比例して
いる。従って、単位体積当たりで dω に入れるモード数は、
D ( ) d  
(  0 )

3
py
3 2
 d
2
px
となり、これが電磁場の状態密度を与える。
実際に遷移が起こる確率は、この状態密度を掛けたものとなる。
表面積 4πr2
半導体中での電子遷移過程
バンド間遷移による光吸収と発光
発光過程
光吸収過程
バンド内遷移光吸収
伝導帯
光
価電子帯
(重い正孔)
価電子帯
(軽い正孔)
スピン分裂
価電子帯
非発光遷移課程
自然放出
誘導放出
価電子帯間吸収
InterValence Band Absorption
オージェ効果
2光子吸収
Two-photon Absorption
半導体pn接合
n型
p型
n型
p型
E
空乏層
電子
正孔
I
Eg: バンドギャップ
電界E 伝導帯
電
子
n
の
エ
EF
ネ
ル
Eg
ギ
p
ー
価電子帯
(a) 熱平衡状態
Vb
n
p
電
EFc
子
の
エ 再結合
eVb
ネ
ル
ギ
ー
Vb = (EFc - EFv) / e
(b) 順バイアス印加時
EFv
半導体内の電子分布
バルク半導体において、伝導帯および価電子帯の電子の
状態密度関数 ρ(E) は、
E
ρc(E)
3
1  2me  2
c (E ) 


2
2   
E  Ec
( ただし , E  E c )
Ec
3
v (E ) 
1  2mh 


2
2   
2
Eg
Ev  E
( ただし , E  E v )
me: 伝導帯における電子の有効質量
mh: 価電子帯における電子(正孔)の有効質量
Ev
放物線状
ρv(E)
0
状態密度関数
ρ(E)
E
伝導帯の電子および価電子帯の正孔の存在確率はフェル
ミ-ディラック分布関数で与えられる。
f (E ) 
1
1  exp
 E  E F  / kT 
k: ボルツマン定数
T: 媒質の温度
EF: フェルミ・レベル
フェルミ・レベル EF は、電子(正孔)の存在確率が1/2と
なるエネルギーレベル。
EF
f(E)
0
1/2
1
フェルミ-ディラック分布関数
熱平衡状態での電子分布
従って、半導体内において熱平衡状態での伝導帯の電子および価電子帯の正孔
密度のエネルギー分布は、各々の状態密度関数にフェルミ-ディラック分布関数を
かけたものとなる。
伝導帯電子分布
E
E
E
ρc(E)
価電子帯
正孔分布
Ec
Ec
EF
Eg
Ev
価電子帯
電子分布
Ev
ρv(E)
ρ(E)
0
状態密度関数
0
1/2
1
f(E)
フェルミ-ディラック分布関数
0
N
伝導帯および価電子帯の
電子分布
電流注入時の電子分布
半導体pn接合に順バイアスを印加して電流注入を行うと、もはや熱平衡状態ではな
くなり、伝導帯電子および価電子帯正孔の存在確率を一つのフェルミ-ディラック分布
関数で表わすことはできなくなり、伝導帯電子および価電子帯正孔に対して各々
別々の分布関数を定義する必要がある。そこで、伝導帯電子に対するフェルミ・レベ
ルを EFc、価電子帯正孔に対するフェルミ・レベルを EFv とし、これらを擬フェルミ・レ
ベル(quasi-Fermi level)と呼ぶ。
伝導帯電子に
伝導帯電子分布
対する分布関数
E
E
E
価電子帯電子に
対する分布関数
ρc(E)
EFc
Ec
価電子帯
正孔分布
Ec
価電子帯
電子分布
Eg
Ev
ρv(E)
ρ(E)
0
状態密度関数
Ev
EFv
0
1/2
1
f(E)
フェルミ-ディラック分布関数
0
N
伝導帯および価電子帯の
電子分布
pn接合界面における電子密度分布
J点
電
子
の
エ
ネ
ル
ギ
ー
E
n
Ec
EF
p
Eg
Ev
E
ρc(E)
正孔
ρv(E)
ρ(E)
0
熱平衡状態のバンド図
E
0
1
f(E)
電子の分布関数
pn接合(J点)での状態密度関数
N
0
電子密度の分布
(a) 熱平衡状態
電
子
の
エ
ネ
ル
ギ
ー
n
J点
p
E
E
E
伝導帯
ρc(E)
EFc
Ec
EFv
Vb = (EFc − EFv) / e
順バイアス時のバンド図
Eg
Ev
0
価電子帯
反転分布
ρv(E)
正孔
ρ(E)
pn接合(J点)での状態密度関数
(b) 反転分布
0
1
電子の分布関数
f(E)
0
N
電子密度の分布
半導体における発光遷移
E
伝導帯電子分布
バンド間遷移による発光スペクトルは、
伝導帯電子および価電子帯正孔の分
布を反映する
バンドテイル
E2
自
然
放
出
光
強
度
Ec
hν = E2− E1
hν
hν
Ev
E1
価電子帯正孔分布
N
0
バンド間遷移による発光
Eg
光子エネルギー(hν)
自然放出光スペクトル
練習問題
問1. 真空中の電子は光と相互作用可能か?
それを考えるには、真空中の電子のエネルギー分散関係と光(フォトン)の
エネルギー分散関係を比較してみること。
まず、真空中の電子のエネルギー分散関係( E-k 曲線)を図にしてみよう。
そこに、フォトンのエネルギー分散関係を重ねて比較してみよう。
真空中の電子のエネルギー分散関係 (E – k 関係)は、
E k  
2
 k
2m
2
m: 真空中の電子の質量(半導体中の場合は、有効質量で
考えればよい。)
フォトンのエネルギー E = hν, フォトンの運動量 p = hν/c
半導体材料の光吸収
105
T = 300K
伝導帯電子分布
Si
光吸収係数: α (cm-1)
E
E2
Ec
hν
hν = E2− E1
104
In0.53Ga0.47As
GaAs
InP
Ge
103
Ev
E1
102
価電子帯正孔分布
N
0
バンド間遷移による光吸収
0.4 0.6 0.8 1.0 1.2 1.4 1.6
波長: λ (μm)
主な半導体材料の光吸収係数
半導体pn接合への光照射
p型領域
電子
正孔
n型領域
−
+
p型領域
−
+
光
光
光起電力: Vph
光
Vph
n型領域
Iph: 光電流
Vb
伝導帯
フェルミ
レベル
価電子帯
Vph: 光照射によって生じる電圧(光起電力)
光
Vb: 逆バイアス電圧
伝導帯
フェルミ
レベル
価電子帯
光照射による光起電力
逆バイアスを印加して光照射
太陽電池(solar cell, photovoltaic cell)
光検出器(photo diode, APD)
光検出器
PINフォトダイオード(PIN-PD)
逆バイアスされた半導体pin接合の
i層に光を照射し、光電流を取り出す
p電極
光
SiNx 反射防止膜
Zn拡散による
p+-InP領域
p+-InP層
−
+
光 Iph: 光電流
n−-InP
ウインドウ層
Vb
i-InGaAs光吸収層
光
n+-InP層
電子
正孔
i-InGaAs −
n -InP層
光吸収層
Vb: 逆バイアス電圧
伝導帯
フェルミ
レベル
n+-InP基板
n電極
価電子帯
InGaAs系 PIN-PDの構造
PIN-PDの動作原理
光検出器
逆バイアスされたpn接合に光が照射されると強度に比例した光電流が取り出せる
逆バイアス
電極
i
10-4
p+
出力信号
n+
光電流
電極
逆バイアス
暗電流: Id (A), 光電流: Iph (A)
光信号
光電流 Iph
10-6
10-8
暗電流 Id
+75℃
+50℃
+25℃
0℃
-25℃
10-10
10-12
光
出力信号
光電流
受光回路
0
20
40
逆方向印加電圧 (V)
PIN-PDの電圧-電流特性
60
光検出器
アバランシェ フォトダイオード(APD)
p電極
光
SiNx 反射防止膜
電子
正孔
InP増倍層 i-InGaAs
光吸収層 p+-InP層
n -InP層
アバランシェ効果により、光電流を増倍
するしくみを有している (高感度)
−
+
p+-InP受光領域
光 Iph: 光電流
n-InP
キャップ層
InP増倍領域
Vb
Vb: 逆バイアス電圧
i-InGaAs光吸収層
n+-InP層
n+-InP基板
n電極
InGaAs系 APDの構造
光
APDの動作原理
伝導帯
フェルミ
レベル
価電子帯
光検出器
アバランシェ フォトダイオード(APD)
10
GB積~40GHz
10-2
10-4
M=1
遮断周波数: fc (GHz)
暗電流: Id (A), 光電流: Iph (A)
アバランシェ降伏(Avalanche breakdown)または雪崩降伏とは、自由電子が電界で
加速され衝突電離を引き起こす過程が繰り返し発生することで、大電流が流れる現象
光電流 Iph
10-6
10-8
暗電流 Id
10-10
降伏電圧 VB
10-12
0
20
40
60
逆方向印加電圧 (V)
APDの電圧-電流特性
80
1
0.1
1
5
10
増倍率 M
APDの高周波特性
50
100
半導体発光素子
電界印加によって得られる発光をエレクトロルミネセンス(Electroluminescence:EL)
と呼び、電流(電子と正孔)を注入して発光させる注入型と、電界によって加速した電
子が何らかの発光中心に衝突し、その発光中心が励起されて発光する真性とに区
別される。注入型としては発光ダイオード(Light Emitting Diode: LED)がある。
発光ダイオード(Light Emitting Diode: LED)
無機半導体
波長帯(色)
近赤外
赤色
橙色
黄色
緑色
青色
紫色
紫外
白色
用途
通信
表示
表示
表示
表示
表示
表示
特殊
照明
材料
InGaAsP, AlGaAs
AlGaAs, GaAsP
GaAsP, AlGaInP, InGaN, GaN , AlGaN
GaAsP, AlGaInP, InGaN , GaN , AlGaN
AlGaInP, InGaN , GaN , AlGaN, GaP
InGaN , GaN , AlGaN, ZnSe, ZnO
InGaN , GaN , AlGaN, ZnSe, ZnO
InGaN , GaN , AlGaN, ZnO, ダイヤモンド (C)
青または紫外のLEDと白色蛍光体との組み合わせ
有機半導体
有機ELなどとして、ディスプレイなどに使われ始めている
発光ダイオード(LED)
発光ダイオード(Light Emitting Diode: LED)
n
電子の拡散長
EFc
発光再結合
eVb
ホールの拡散長
青色LEDの外観
p
EFv
LEDの構造
半導体レーザー
半導体レーザー (Laser Diode: LD)
光共振器を有しており、光の誘導放出を利用して発振動作する
特徴: ・ コンパクト (チップ本体は0.3mm角程度)
・ 取り扱い容易 (乾電池2本程度で動作可能)
・ 直接変調で数Gbpsの高速変調が可能
・ 高信頼性 (通信用のInGaAsPレーザは100万時間以上の寿命に)
・ 安価 (FTTH用LDはチップコストで数百円、CD用LDは数十円に)
Blu-ray用
青紫色LD
(三洋)
へき開面(鏡面)
LDチップの構造
出展: www.phlab.ecl.ntt.co.jp/master/04_module/002.html
光通信用バタフライモジュール
二重ヘテロ接合活性層
pn接合ダイオードに順方向に
電流注入すると、電子は p型
領域に、正孔は n型領域に
各々拡散していくが、それら
キャリヤの拡散長は数 μmあ
り、pn接合界面を挟んで数
μm程度の広い領域で発光再
結合が起こる。従ってホモ接
合では、キャリヤ密度が低い
ために効果的な発光再結合
は得にくい。
一般的なLDの活性層は、二
重ヘテロ(Double-hetero)構
造によってキャリヤ(電子と正
孔)および光を活性層に効果
的に閉じ込め、高いキャリヤ
密度により高い発光再結合
効率を実現している。
n型
電
子
の
エ
ネ
ル
ギ
ー
電
子
の
エ
ネ
ル
ギ
ー
(a) 構造
p型
活性層
n
EF
(b) 熱平衡状態
p
n
EFc
p
発光再結合
(c) 順バイアス時
EFv
光の強度分布
屈
折
率
(d) 屈折率分布
半導体量子井戸
活性層の厚さを電子のド・ブロイ波長(数nm)程度まで
第三準位
薄くすると、キャリヤが狭い領域に閉じこめられるため、
状態密度が離散化される。これを量子サイズ効果と呼
第二準位
ぶ。電子および正孔の波動関数およびエネルギー準
位も量子化される。このような構造を半導体量子井戸
第一準位
(Quantum Well: QW)構造と呼ぶ
各量子準位の波動関数
InGaAsP
量子井戸
p-InP
n-InP
InGaAsP
n-InP 量子井戸
p-InP
電子
階段状
ρ(E)
正孔
量子井戸構造
バンド構造
状態密度
キャリヤ閉じ込めの次元と状態密度との関係
p-InP
p-InP
p-InP
bulk InGaAsP
QW
n-InP
バルク活性層
E
ρc(E)
p-InP
QWire
QD
n-InP
n-InP
n-InP
量子井戸活性層
量子細線活性層
量子ドット活性層
E
ρc(E)
E
E
ρc(E)
ρc(E)
階段状
放物線状
ρv(E)
0
ρv(E)
ρv(E)
ρ(E) 0
線状
鋸歯状
ρ(E) 0
状態密度関数
ρv(E)
ρ(E) 0
ρ(E)
多重量子井戸(MQW)活性層
p-InP
MQW活性層
活性層
MQW
n-InP
MQW活性層構造
n-InP
InGaAs/
InGaAsP
MQW
MQW活性層の断面SEM写真
E
p-InP
電子
ρ(E)
EF
軽い正孔 重い正孔
バンド構造
状態密度
歪量子井戸とは
通常エピタキシャル成長では、基板の格子定数に成長させる物質の格子定数を合わ
せて結晶成長を行う。例) (100)InP基板 (格子定数: Λ=5.869Å)にIn0.53Ga0.47Asなど。
エピタキシャル
成長層 (In0.53Ga0.47As)
InP 基板
Λ
歪量子井戸とは、基板とは僅かに異なる格子
定数の材料を成長させて、格子歪を導入
In0.53Ga0.47As
エピタキシャル
成長層 (In0.6Ga0.4As)
InP 基板
Λ
Ⅲ-Ⅴ族化合物半導体の格子定数
格子歪による半導体バンド構造の変化
S. L. Chuang, Phys. Rev. B 43, 9649 (1991)
z
y
x
格子歪による In1-xGaxAs バンド構造の変化
歪量子井戸活性層
量子井戸活性層に格子歪を導入することにより、(特に価電子帯の)状態密度関数を
望ましい形に設計することが可能。
1. 面内(2軸性)圧縮歪を導入した(x < 0.468)場合
価電子帯では、Γ点での重い正孔と軽い正孔のエネルギー縮退が解けて分裂す
る。その結果、伝導帯と重い正孔のバンド間での遷移が支配的となる。エネル
ギー縮退が解けることにより、バンド混合の効果が緩和されて、 Γ点付近での重
い正孔の有効質量が小さくなる。その結果、僅かに正孔を注入すれば反転分布
が得られるようになる。従って、発振閾値の低いLDを実現するには好都合。
2. 面内(2軸性)引っ張り歪を導入した(x > 0.468)場合
価電子帯では、 Γ点での重い正孔と軽い正孔のエネルギー縮退が解けて分裂す
る。その結果、伝導帯と軽い正孔のバンド間での遷移が支配的となる。伝導帯と
軽い正孔のバンド間での遷移は、活性層に垂直な電場成分を有する(TM-like
モードの)光と相互作用をするため、TM-likeモードの光に対して光学利得を有す
る。TEとTMの両方の光に対して増幅を行いたい光増幅器などに用いられる。
このように、本来の材料の有する物性(バンド構造)を造り替えて、好ましい状態密度
関数の形に自由にデザインする技術をバンドエンジニアリングという。量子井戸、超
格子、歪量子井戸などはその典型である。
バンドエンジニアリングは、現代の科学的な錬金術である !!
半導体レーザの発振条件
利得条件
光共振器を光が往復した時、元の光強度と同じになる条件
2  g  L
1
つまり、 R f R r e
位相条件(共振条件)
光共振器を光が往復した時、光電場の位相が一致する条件
つまり、 2  L   f   r  2 q 
反射率: Rr
光出力
1−Rr
q は自然数, β は半導体内での光の伝搬定数
ϕf , ϕr は端面反射による位相回転
光増幅利得: g cm-1
光出力
1−Rf
光吸収係数: α
共振器長: L
cm-1
反射率: Rf
半導体レーザの光共振器構造
Fabry-Perot (FP)共振器レーザー
2枚の平行に向き合った鏡によるFP型光共振器
によって正帰還が得られ発振するレーザー
発振波長間隔
0
2
 
縦多モード発振
Δλ
2 n eff L
λ0 : 発振波長の中心値
neff : 実効屈折率
L : 素子長
λ0
発振スペクトル

へき開面(鏡面)
FPレーザーの構造
回折格子
分布帰還(DFB)型レーザー
回折格子によるBragg反射により、光の分布帰還
が得られ、 Bragg波長近傍の単一波長で発振
発振波長
単一縦モード発振
  2 n eff 
Λ : 回折格子の周期
neff : 実効屈折率
発振スペクトル

DFBレーザーの構造
出展: www.matsuoka-lab.imr.tohoku.ac.jp/purposes.html
半導体レーザの光共振器構造
半導体レーザの各種光共振器構造
へき開面
へき開面
反射防止膜
反射防止膜
活性層
光出力
光出力
光出力
活性層
光出力
回折格子
FP型
λ/4位相シフト
DFB型
光出力
反射防止膜
DBR反射鏡
反射防止膜
活性層
活性領域
光出力
光出力
回折格子
VCSEL型
DBR型
半導体レーザの発振特性
0 10 20 40
60 80
100℃
15
光出力 (mW)
120℃
10
5
0
0
50
100
150
注入電流 (mA)
1.3μm InGaAsP MQW FPレーザの電流-光出力(L-I)特性
半導体レーザの発振特性
発振状態
1302
光強度 (10dB/div)
光強度 (10dB/div)
1.3μm InGaAsP MQWレーザの発振スペクトル
1312
波長 (nm)
FPレーザ
1322
発振状態
(SMSR= 50dB)
1298
1308
波長 (nm)
DFBレーザ
1318
半導体レーザの発振特性
光出射特性
y
x
z
レーザビーム
基板平行方向光強度分布
y
-θ
0
+θ
-θ
基板垂直方向
0
基板平行方向
(a) FFP
x
+θ
基板垂直方向光強度分布
(b) NFP
光変調特性
光信号
光
出
力
Ith
Ib
電流
変調信号(電気)
規格化変調光強度 (3dB/div)
半導体レーザの駆動電流に変調信号を重畳し、直接光出力を変調する方式を直接変調
と呼ぶ。
Ib=1.04 Ith1.10 Ith 1.23 I
th
1.37 Ith
0
2
4
6
1.64 Ith
8
変調周波数: f (GHz)
小信号変調応答周波数特性
半導体レーザの電流-光出力(I-L)特性
10
光変調特性
大信号変調特性
I
注
入
電
流
Ib + Ip
Ith
Ib
t
(a)
活
性
層
内
キ
ャ
リ
ヤ
密
度
n
キャリヤ寿命時間 τs での減衰
nth
nb
t
(b)
S
光
子
密
度
振動周期はほぼ 1/fr
Sop
0
t
(c)
通信用途以外での半導体レーザ
世界初、発振波長530nm帯で100mW以上の光出力を有する純緑色
半導体レーザーを開発
~ 従来の窒化ガリウム系緑色レーザーに比べて、約2倍の高輝度を実現 ~
住友電気工業株式会社
ソニー株式会社
2012年6月21日プレスリリース
http://www.sony.co.jp/SonyInfo/News/Press/201206/12-0621/
光増幅器
半導体光増幅器(Semiconductor optical amplifier: SOA)
半導体レーザーの両端面に無反射膜
を形成するなどして、光共振器をなく
したもの (光の正帰還がかからなくな
るのでレーザー発振はしない)
無反射加工
無反射加工
半導体レーザーチップ
特徴
・ 小型で高効率
・ 電流注入により光増幅
・ 光増幅利得 ~ 20dB
・ 増幅波長帯域 30 ~ 40nm
・ 自然放出光雑音大
・ 偏光依存性有り
通信用SOAモジュール
光増幅器
光ファイバー増幅器(Optical fiber amplifier)
Er添加光ファイバー増幅器
コアに、エルビウム(Er3+)などの希土類を添加
波長980 nmや1480 nmなどの光で励起すると
波長1.54 μm付近に光増幅利得発生
Er3+の準位
特徴
・ 光増幅利得 ~ 30dB
・ 増幅波長帯域 ~ 30nm
・ 自然放出光雑音小
・ 偏光依存性無
・ 光励起
・ 低効率
光ファイバー増幅器の構成
市販の光ファイバー増幅器
光増幅器
ラマン増幅器(Raman Amplifier)
光ファイバーに非常に強い励起光を入射すると、石英ガラスの分子振動エネルギー
に対応して、誘導ラマン散乱(Stimulated Raman Scattering: SRS)が生じ、励起光
波長より100 nm程度長い波長域に光学利得が得られる。
信号光 増幅
中間状態
(仮想準位)
反転分布
励起光
特徴
・ 励起光波長を選べば、任意の波長で
増幅が可能
・ 偏光依存性無
・ 光励起、低効率
励起状態
物質固有の振動エネルギー
基底状態
SRSによる光増幅原理
光ファイバーラマン増幅器の構成
波長可変レーザ
波長多重伝送やフォトニックネットワークでは、発振波長を自由に変えられる波長可変
レーザが望まれる。波長可変レーザは、光増幅器と波長選択素子とで構成できる。
コリメート レンズ
回折格子
半導体光増幅器
回折格子による機械式波長可変レーザの構成
マイクロヒーター
SOA
リング共振器
Si細線導波路
スポットサイズ
変換
リング共振器による熱光学式波長可変レーザ
光変調器
光通信において、半導体レーザーなどは直接変調が可能であるが、直接変調を行え
ない光源(ガスレーザーや固体レーザーなど)には、外部で光変調を行うための光変
調器がある。
光変調には、
光の強度を変調する光強度変調器
1
アナログ変調
0
1
0
デジタル (バイナリ)変調
光の位相を変調する光位相変調器
0
アナログ変調
1
0
デジタル (バイナリ)変調
光変調器
電界吸収(EA)型光強度変調器
化合物半導体などのpn接合に逆バイアスを印加すると、
印加電界によって光吸収特性が変化(フランツ・ケル
ディッシュ効果、量子閉じ込めシュタルク効果)し、これを
利用して光の強度変調を行うもの。
40GbpsEA変調器(沖電気)
Franz–Keldysh effectの説明
quantum-confined Stark effect(QCSE)の説明
光変調器
ニオブ酸リチウム(lithium niobate: LN) によるMZI型光強度変調器
ニオブ酸リチウム(LiNbO3)は、電圧を印加すると屈折率が変化する電気光学
(E-O)効果(Pockels effect)を有している。従って、 LNによる光導波路では、伝
搬光の位相変調が可能。LNによる光導波路によってMach-Zehnder(MZ)型の
光干渉計を構成し、光の位相変化を強度変化に変換し光強度変調も行える。
光の位相を変調
LN光導波路
電極
強度変調された出射光
M-Z干渉計
入射光
LN光変調器の原理
LiNbO3(LN)光変調器
(住友大阪セメント)
光変調器
LN光変調器の構造
住友大阪セメント(株)の資料より
光分岐(光スプリッター)
Y分岐
0次モード
1次偶モードのみ
選択的に励振される
出射光
入射光
断熱的モード変換部
出射光
MMI(multi-mode interference)カプラー型分岐
多モード干渉導波路
入射光
出射光
光分岐(光スプリッター)
光方向性結合器
0次モード
1次偶モードのみ
選択的に励振される
出射光
上に比べてπ/2だけ
位相が遅れて出射
入射光
方向性結合器
完全結合長
出射光
P=
P×
完全結合長
分岐比が50:50になる
P=
結合導波路長: L
奇モード 偶モード
方向性結合器の動作
P×
結合導波路長: L
光合分波器
波長多重光伝送においては、異なる波長の光搬送波を1本の光ファイバーを用いて
伝送する。従って、異なる波長の光を束ねる(合波器)と、光を波長によって分ける (分
光器あるいは分波器)が不可欠。
Arrayed Waveguide Grating
クラッド
この一本一本が
コア
このような光導
波路からなる
0.5 m
0.5 m
Si 基板
石英光導波路
50 mm
Arrayed Waveguide Grating (AWG)
スラブ
導波路   
12
N
AWGの動作原理
導波路の波長分散を利用
光合分波器
波長によって導波路を伝搬する光の速度が異なる
λ21
光分波特性
NTTエレクトロニクス社製品
光スイッチ
電気制御-光スイッチ (光の経路を切り換えるが、ON-OFFの制御は電気で行う)
スイッチング機構
特 徴
出力ファイバー
メカニカル
(MEMS)
Port1
Port2
入力ファイバー
入力1
熱光学(T-O)効果
電気光学(E-O)効果
出力1
mS~Sオーダーの切換え速度
比較的安価
ヒーター
入力2
+
mSオーダーの遅い切換え速度
安価
出力2
-
電界印加
nSオーダーの高速切換え
高価
その他に、磁気光学(M-O)型、音響光学(A-O)型などもある
光制御-光スイッチ (光-光スイッチ or All光スイッチ)
ON-OFF制御も光でやる
現在研究開発中 将来のフォトニックネットワークや全光信号処理システムに
使われるかも?
Si細線導波路による熱光学(T-O)光スイッチ
マイクロヒーター
チップサイズは僅か1.4 mm×2 mm
電極
1.4 mm
Port1
Port2
Si細線導波路
T ran sm issio n (d B )
MZ型Si細線導波路光スイッチ素子
2 mm
-2 0
-3 0
-4 0
250 m
P o rt 1
P o rt 2
0
50
100
150
H eatin g p o w er (m W )
スイッチング特性
Port8
200
1×8光スイッチの写真
T. Chu et al., Optics Express 13, 10109 (2005)
T. Chu et al., Proc. SPIE 6477 (2007)
光アイソレーター
通信用デバイスでは、戻り光が様々な悪影響を及ぼすことがある。従って、それを
除去する光アイソレーターが必要。磁気光学素子と偏光子により、光を一方向にの
み通過させる素子を光相反素子または光アイソレーターと呼ぶ。
光アイソレーターの原理
各種光アイソレーター
住友金属鉱山(株)HPより
分散補償素子
光ファイバーや通信用光デバイスの中では、波長によって光の伝搬速度が異なる
ために、光信号波形が乱されることがある。波長によって光の伝搬速度がことなる
現象を chromatic dispersion と呼ぶ。これを補償する素子が分散補償 素子である。
分散補償素子
光ファイバーを用いたもの
分散補償光ファイバー(dispersion compensation fiber: DCF)
通常の光ファイバーとは逆の分散特性を有する光ファイバー
AWGとSLM(空間位相変調器)を
用いたもの
SLMとしては液晶による光位相
変調器LCOSを使用
LCOS: Liquid Crystal on Silicon
大庭他、NTT技術ジャーナル2010.2
光ヘテロダイン検波
コヒーレント光伝送方式では、光ヘテロダイン/ホモダイン検波などの光検波方式が
用いられる
合波器(BS)
信号光
fS
受光器(PD)
IFアンプ/
フィルター
fIF=fS-fLO
ベースバンド
復調器
局部発振光
fLO
周波数弁別器
LD
光ヘテロダイン検波回路のブロック図
信号光は局部発振光と混合(合波)されて共に受光器(PD)に入る。従って信号光の
電界の位相情報がPDの検波出力電流に保存される。
通常、LDの発振波長(周波数)は揺らいでいるため、LDからの局部発振光と信号
光との周波数差(中間周波数)をモニターし、局部発振LDの駆動電流にフィードバッ
クをかけることにより、局部発振光LDの発振波長を信号光にロックする光PLLなど
の手法が用いられる
デジタルコヒーレント光伝送
デジタルコヒーレント方式の概要
出典: 信学会誌H24年12月、号 総合報告
フォトニックネットワーク
光リンクを流れる光信号を一旦電気に変換することなく、光のまま交換する
次世代の光ネットワーク
- Optical Add/Drop Multiplexer(OADM)
OADM
- Optical Cross Connect(OXC)
OADM
OXC
WDM ring NW
WDM ring NW
OADM
OADM
リング型ネットワーク
- λ-MPLS, G-MPLS
(MPLS: Multi-protocol label switching)
OPS
Mesh NW
- Optical Burst Switching(OBS)
OPS
OPS
- Optical Packet Switching(OPS)
OPS
メッシュ型ネットワーク
Optical Add/Drop Multiplexer(OADM)
波長多重化された信号の中から、ある特定波長(チャンネル)の信号のみ取り出したり、
加えたりするもの
OADM
1 ‥‥ n
OADM
OADM
WDMリングNW
OADM
OADM
WDM信号
i
i
R-OADM (Reconfigurable OADM) Add/Dropする波長を任意に設定できるもの
1 ‥‥ n
OADM
OADM
WDM信号
OADM
WDMリングNW
OADM
i
i
OADM
波長合分波器(AWG)
N×N AWG 1個で N×N波長ルーターが構成可能
SiO2 クラッド
SiO2 コア
0.5 m
λ1, λ2, λ3, …, λN
λ1, λ2, λ3, …, λN
λ1, λ2, λ3, …, λN
λ1, λ2, λ3, …, λN
λ2, λ3, λ4, …, λ1
λ3, λ4, λ5, …, λ2
λ1, λ2, λ3, …, λN
λN, λ1, λ2, …, λN-1
0.5 m
Si 基板
このような石英光導
波路からなるAWG
シリコン光導波路を用いれば、
極微小なAWGが実現可能
50 m
大きさ1/1000
50 mm
Arrayed Waveguide Grating (AWG)
Si細線光導波路によるAWG
四光波混合による40Gb/s波長変換
Y.-H. Kuo et al., Optics Express 14, 11721 (2006)
逆バイアスされた 8cm長 SOI pin リブ導波路による40G波長変換
波長変換効率: -8.6dB
Pump power: 450 mW
40Gbps NRZ波長変換実験系
波長変換出力スペクトル
40Gbps Eye Diagram
左: 入力信号
右: 波長変換出力信号
光による符合ラベル処理
グレーティング光ファイバーを用いる方式
1 4 3 2
2 1 3 4
データ
データ
データ
t
波長ラベル
波長ラベル
2
3
4
1
サキュレータ
グレーティング光ファイバー
グレーティングパターンとラベルが一致した場合
一致しないと
データ
データ
データ
t
光バッファ
1. 光遅延線路と光スイッチによる
光遅延線路
光遅延線路
光遅延線路
光スイッチ
光スイッチ
光スイッチ
2. Slow Lightによる
電磁誘導透過 EIT: Electromagnetically Induced Transparency
0.9μK(約−273℃) ナトリウム
300,000km/s → 28m/s
70~90K(-203~-183℃)
ルビジウムRb
300,000km/s → 1km/s
|3>
|1>
|2>
レポート課題
私の講義(全5回分)の内容に基づいて、通信用光デバイスに関して学んだこと
を、各自でテーマを決めて、A4レポート用紙10頁以内(図なども含む)で述べよ。
レポートには、選んだ課題のテーマ名を必ず明記すること。
例) 半導体レーザーのモード競合雑音について、、、、など。
評価のポイントは、選んだテーマとは関係無く、自分なりの考えや主張に基づ
いた起承転結のある論旨の明確なレポートを高く評価します。Webなどからの
切り貼りや、何を主張したいのか分からない、結論が無い論文やレポートにつ
いては評価は低くなります。また、講義に関するご意見や感想も歓迎します。
提出〆切: 8月8日(木)17:00
提出先: 私のメールボックスまたは研究室(青葉山 電気系2号館2階202号室)