4/19 No. 1
Advanced Plasma and Laser Science プラズマ・レーザー特論E
Quick review of quantum mechanics
量子力学の復習
Takeshi Sato and Kenichi Ishikawa
http://ishiken.free.fr/english/lecture.html [email protected]
[email protected]
4/19 No. 2
Hydrogen atom 水素原子の波動関数 Atomic unit 原子単位
Rabi oscillation ラビ振動
4/19 No. 3
Hydrogen-like atom
水素原子の波動関数
4/19 No. 4
Schrödinger equation シュレーディンガー方程式
ポテンシャル
V(r)
中の質量m
の電子€
i ∂ψ
∂ t = −
22 m ∇
2ψ (r,t ) + V (r ) ψ (r,t )
steady state
定常状態€
−
22 m ∇
2ϕ (r) + V (r) ϕ (r) = εϕ ( r)
: Energy eigenvalue
エネルギー固有値(エネルギー準位): Eigen function
固有波動関数:
Wave function
Eigenvalue problem
固有値問題€
ψ (r,t ) = ϕ (r)e
−iωt€
ε = ω
€
ψ (r ,t ) Particle of mass m moving in a potential V(r)
:波動関数
€
ϕ (r )
4/19 No. 5
Hydrogen-like atom 水素様原子
原子核のクーロンポテンシャル
€
V (r) = V (r) = − Ze
24 πε
0r
€
−
22m ∇
2ϕ (r) − Ze
24 πε
0r ϕ (r) = εϕ (r)
係数が煩雑
Introduction of atomic unit (a.u.)
原子単位の導入€
− 1
2 ∇
2ϕ (r) − Z
r ϕ (r) = εϕ (r)
(Time-independent Schrödinger equation)
シュレーディンガー方程式Bare Coulomb potential from the nucleus
cumbersome coefficients
4/19 No. 6
Atomic unit 原子単位
€
= m = e = e
24 πε
0= 1
となるような単位系長さ
€
a
0=
2m e
24 πε
0$
% & ' ( )
= 4πε
0
2me
2= 5.292 ×10
−11m
ボーア半径
時間
エネルギー
€
e
24 πε
0a
0= 27.21 eV
€
1 eV = 1.602 × 10
−19J
€
3m e
24 πε
0$
% & ' ( )
2
= a
0α c = 0.0242 fs
€
α = e
24 πε
0c = 7.297 ×10
−3= 1 137.0
微細構造 定数
速度
€
a
0÷ a
0αc = α c
Electron
電子Unit system in which Length
Energy
Time
Velocity
Bohr radius
2 (ionization potential of H)
fine structure constant
Atomic scale of length,
energy, and time
4/19 No. 7
Atomic unit is closely related to Bohr hydrogen atom
Dimension Expression Value Meaning
length 5.29 10 -11 m Bohr radius
energy 27.2 eV
Coulomb potential energy at the Bohr radius
velocity 2.19 10 6 m/s electron orbital
velocity
time 24.2 attoseconds
time during which the electron
proceeds 1 radian
electric field 5.14 10 11 V/m field at the Bohr
radius laser
intensity 3.51 10 16 W/cm 2
laser field = electric field at the Bohr radius
a 0 = 4 0 2 /me 2 E h = me 4
(4 0 ) 2 = e 2 4 0 a 0 v = e 2
4 0 = c E h = a 0
v F = e
4 0 a 2 0 1
2 c 0 F 2
4/19 No. 8
Hydrogen-like atom 水素様原子
原子核のクーロンポテンシャル
€
V (r) = V (r ) = − Ze
24 πε
0r = − Z r
€
−
22m ∇
2ϕ (r) − Ze
24 πε
0r ϕ (r) = εϕ (r) Polar coordinate
極座標系€
r = (r, θ , φ )
固有波動関数
€
ϕ (r) = R
nl(r)Y
lm( θ , φ ) Bound state
束縛状態€
ε < 0
エネルギー固有値
€
n = 1,2,3
動径波動関数
Spherical harmonics
球面調和関数€
0 ≤ l ≤ n − 1
€
− l ≤ n ≤ l Bare Coulomb potential from the nucleus
(Time-independent Schrödinger equation)
シュレーディンガー方程式€
− 1
2 ∇
2ϕ (r) − Z
r ϕ (r) = εϕ (r)
€
ε
n= − Z
2me
44 πε
0( )
22
21
n
2= − Z
22n
2Energy eigenvalue
Eigen function
Radial wave function
4/19 No. 9
Bound states 束縛状態
エネルギー固有値
€
ε
n= − Z
2me
44πε
0( )
22
21
n
2= − Z
22n
2€
n = 1,2,3
€
ε
1= − me
44 πε
0( )
22
2= −13.6 eV
Ground state
基底状態r in a 0 (Bohr radius
ボーア半径)
€
a
0= 4 πε
0
2me
2= 5.3 ×10
−11m = 0.053 nm
€
ϕ (r) = R
nl(r )Y
lm( θ , φ )
€
0 ≤ l ≤ n −1
€
−l ≤ n ≤ l
1s
2s, 2p
3s, 3p, 3d
En e rg y ( e V)
Coulomb potential
Energy eigenvalue
4/19 No. 10
エネルギー固有値
€
ε
n= − Z
2n
2€
n = 1,2,3 Energy eigenvalue
Balmer series
Lyman series
4/19 No. 11
Radial wave function and spherical harmonics 動径波動関数と球面調和関数
€
R
1s= 1 a
0"
# $ %
&
'
3 / 2
2e
−r/a0€
R
2s= 1 a
0"
# $
%
&
'
3 / 2
1
2 e
−r/ 2a01 − r 2a
0"
# $
%
&
'
€
R
2p= 1 a
0"
# $ %
&
'
3 / 2
1
2 6 e
−r/ 2a0r a
0€
R
3s= 1 a
0"
# $
%
&
'
3 / 2
2
3 3 e
−r/ 3a01 − 2 3
r
a
0+ 2 27
r a
0"
# $
%
&
' )
2* + +
, - . .
€
R
nl∗( r)R
n l#(r)
0
∫
∞r
2dr = δ
nn #Orthonormality
規格直交性Z = 1
の場合€
Y
00= 1 4 π
€
Y
1,0= 3
4π cos θ
€
Y
1,±1= 3
8 π sin θ e
±iφ€
Y
2,0= 5
16π ( 3cos
2θ − 1 )
€
Y
2,±1= 15
8π sin θ cos θ e
±iφ€
Y
2,±2= 15 32 π sin
2
θ e
±2iφOrthonormality
規格直交性€
Y
lm∗( θ , φ )
∫ Y
l & m &( θ , φ )sin θ d θ d φ = δ
ll &δ
mm &€
ϕ
nlm∗∫ ϕ
n % l % m %r
2sin θ drd θ d φ = δ
nn %δ
ll %δ
mm %4/19 No. 12
Probability density
存在確率密度Radial wave function
動径波動関数r (atomic unit)
1s 2s
2p
3s
3p
3d
4/19 No. 13
Continuum states 自由状態、連続状態
€
ε > 0
€
ϕ (r) = R
εl(r)Y
lm( θ , φ )
€
ε > 0
€
l ≥ 0
€
−l ≤ n ≤ l
Arbitrary positive number
任意の正の実数Necessary when ionization is considered
イオン化を考えるときに必要€
R
εl(r) = 2 Z
1 − e
−2πn %s
2+ n %
2s=1 l
∏ (2l (2 kr) + 1)!
le
−ikrF (i n % + l +1,2l + 2,2ikr )
€
k = 2 mE / = 2 E
€
"
n = Z k
Radial wave function
動径波動関数→ Coulomb wave function
クーロン波動関数合流型超幾何関数
€
R
ε∗l(r )R
ε $ l(r )
0
∫
∞r
2dr = 0
€
ε ≠ $ ε
€
R
εl( r)
20
∫
∞r
2dr > 0 Density of states
状態密度confluent hypergeometric function
wave number
波数4/19 No. 14
Radial wave function 動径波動関数
Continuum states
自由状態(連続状態)Bound states
束縛状態r
の単位はa 0 (
ボーア半径)
4/19 No. 15
Coulomb wave function vs. electron in a free space V ( r )=0 クーロン波動関数と自由空 間の電子波動関数のとの比較
€
− 1
2 ∇
2ϕ (r) = εϕ (r)
€
− 1 2
d
2dr
2+ 2
r d
dr − l (l +1) r
2#
$ % &
' ( R(r ) = ε R(r)
€
V (r) = 0
In a free space
€
R
El( r) = 2 k
π j
l( kr) %
r→∞% % → 2 π k
1
r cos kr − π
2 (l +1) '
( ) * + , Spherical Bessel function
Coulomb wave function
€
R
El(r) $
r→∞$ $ → 2 π k
1
r cos kr + Z
k log 2kr − π
2 (l +1) − σ
l(
) * +
, -
Phase shift
位相シフト(位相のずれ)€
σ
l= arg Γ(l +1 + iZ / k )
10 20 30 40 50
0.5 0.5
r
€
rR
El(r)
€
l = 1 (p-wave)
E = 13.6 eV
Couomb
V(r)=0
4/19 No. 16
Short-range potential V ( r )=0 at r > r 0
短距離ポテンシャル
€
− 1
2 ∇
2ϕ (r) = εϕ (r)
€
− 1 2
d
2dr
2+ 2
r d
dr − l (l +1) r
2#
$ % &
' ( R(r ) = ε R(r)
€
V (r ) = 0
j
l(kr ) "
r→0"" → (kr )
l(2l +1)!! "
r→∞"" → 1
kr cos kr − π
2 (l + 1)
%
&'
( )*
y
l(kr ) "
r→0"" − → (2l − 1)!!
(kr )
l+1"
r→∞"" → 1
kr sin kr − π
2 (l + 1)
%
&'
( )*
Spherical Bessel function
Phase shift
位相シフト(位相のずれ)r > r
010.48 Graphs 263
Figure 10.48.1: jn(x), n= 0(1)4,0 x 12. Figure 10.48.2: yn(x), n= 0(1)4,0< x 12.
Figure 10.48.3: j5(x),y5(x), j25(x) +y25(x), 0 x 12. Figure 10.48.4: j5(x), y5(x), ⇥
j52(x) +y52(x), 0 x 12.
Figure 10.48.5: i(1)0 (x), i(2)0 (x),k0(x), 0 x 4. Figure 10.48.6: i(1)1 (x),i(2)1 (x),k1(x), 0 x 4.
Figure 10.48.7: i(1)5 (x), i(2)5 (x), k5(x), 0 x 8.
R
El(r ) = 2k
π ( c
jj
l(kr ) + c
yy
l(kr ) )
4/19 No. 17
Temporal evolution by an external field 外場との相互作用による時間発展
€
i ∂ψ
∂t = − 1
2 ∇
2ψ (r,t ) − Z
r ψ (r,t ) + V
I(r,t ) ψ (r,t )
相互作用
Interaction
Without the external field
相互作用項がない場合€
ψ
n(r,t ) = ϕ
n(r)e
−iωnt€
ω
n= ε
n Eigen state
固有状態€
i ∂ψ
∂ t = ( H
0+ H
I) ψ (r,t )
€
H
0= − 1
2 ∇
2− Z r
€
H
I= V
I(r,t )
€
H
0ϕ
n(r) = ε
nϕ
n(r)
With the external field
相互作用項がある場合€
ψ ( r,t ) = c
nϕ
n(r)e
−iωntn
∑
€
c
n= e
iωnt∫ ϕ
n*(r) ψ ( r,t )dV = e
iωntn ψ
€
H
0n = ω
nn (atomic unit)
4/19 No. 18
€
i ∂
∂ t ψ = ( H
0+ H
I) ψ
€
i ∂
∂t n ψ = n H
0+ H
Iψ = n H
0ψ + n H
Iψ = ω
nn ψ + n H
Iψ
€
n ψ = c
ne
−iωnt€
i c ˙
n= n H
Iψ e
iωnt€
m m
m
∑ = I Identity operator
単位演算子can be inserted anywhere
€
i c ˙
n= n H
Im m ψ e
iωntm
∑ = n H
Im c
me
i(ωn−ωm)tm
∑
€
i c ˙
n= n H
Im c
me
i(ωn−ωm)tm
∑
€
n H
Im Transition matrix element
遷移行列要素Image
イメージTransition from m to n due to the interaction H I
状態
m
が相互作用H I
によって状態n
に遷移するThe interaction H I couples m to n.
4/19 No. 19
Important example: Rabi oscillation 重要な例:ラビ振動
€
ω
0€
ε
2€
ε
1Resonance frequency
遷移振動数(共鳴振動数)€
ω
0= ε
2− ε
1 2準位系Two-level atom
2準位系€
ψ ( r,t ) = C
1(t ) ψ
1( r,t ) + C
2(t ) ψ
2(r,t )
€
C
2 2€
C
1 2 光の振動数がω 0
に近いときは、放射過程に関与するのは選ばれた二 つの原子状態のみ。
If the laser frequency ω is close to ω 0 ,
only the two levels are relevant.
4/19 No. 20
€
i ∂ψ
∂t = − 1
2 ∇
2ψ (r,t ) − Z
r ψ (r,t ) + V
I(r,t ) ψ (r,t )
€
ψ (r,t ) = C
1(t ) ψ
1(r,t ) + C
2(t ) ψ
2(r,t )
€
ψ (r,t )
2d
3r
∫ = C
1(t )
2+ C
2(t )
2= 1
€
V
I( C
1ψ
1+ C
2ψ
2) = i ∂ C
1∂ t ψ
1+ ∂ C
2∂ t ψ
2$
% & ' ( )
€
ψ
1∗ を左からかけて空間積分€
i ∂ C
1∂ t = C
1V
11+ C
2V
12e
−iω0tSimilarly
同様に
€
V
ij= i V
Ij = ∫ ϕ
i∗V
Iϕ
jd
3r
€
i ∂ C
2∂t = C
1e
iω0tV
21+ C
2V
22€
ω
0€
ε
2€
ε
1€
C
2 2€
C
12multiply with from the left and take a volume integral
€
ψ
1∗4/19 No. 21
Interaction Hamiltonian 相互作用ハミルトニアン
Complete Hamiltonian for the interaction of an atom with an electromagnetic field is rather complicated. 電磁場と原子 の間の相互作用に対するハミルトニアンの完全な形は複雑
x
y z
€
E
0cos( kx − ω t )
€
H
0cos( kx − ω t ) Ze k
r
€
k = 2 π
波数λ
Wavelength
波長€
x << λ
€
kx << 1
€
E
0cos(kx − ω t )
€
E
0cos ω t
€
V
I= zE
0cos ω t Dipole approximation
電気双極子近似Dipole approximation is often sufficient.
レーザーに関しては、多くの場合、電気双極子近似で十分
(原子単位)
Wave number
4/19 No. 22
€
V
I= zE
0cos ω t
€
i ∂ C
1∂ t = C
1V
11+ C
2V
12e
−iω0t€
V
ij= i V
Ij = ∫ ϕ
i∗V
Iϕ
jd
3r = cos ω t zE ∫
0ϕ
i∗ϕ
jd
3r = X
ijcos ω t
€
i ∂ C
2∂ t = C
1e
iω0t
V
21+ C
2V
22€
X
11= X
22= 0
€
i ∂ C
1∂ t = 2 γ C
2e
−iω0tcos ω t
€
i ∂ C
2∂ t = 2 γ C
1e
iω0tcos ω t
€
i ∂ C
1∂ t = γ C
2[ e
i(
ω−ω0)
t+ e
−i(
ω+ω0)
t]
€
i ∂ C
2∂t = γ C
1[ e
i(
ω+ω0)
t+ e
−i(
ω−ω0)
t]
€
X
12= X
21= 2 γ
(Real
実数)€
V
ij= i V
Ij = ∫ ϕ
i∗V
Iϕ
jd
3r
How V I couples the two levels.
「
V I
のおかげでj → i
に遷移する」割合4/19 No. 23
Rabi oscillation ラビ振動
回転波近似
€
i ∂ C
1∂ t = γ C
2[ e
i(
ω−ω0)
t+ e
−i(
ω+ω0)
t]
€
i ∂ C
2∂t = γ C
1[ e
i(
ω+ω0)
t+ e
−i(
ω−ω0)
t]
€
i ∂ C
1∂t = γ e
i(
ω−ω0)
tC
2€
i ∂ C
2∂t = γ e
−i(
ω−ω0)
tC
1 初期条件€
C
1= 1, C
2= 0
€
C
1(t ) = cos Ω t − i ( ω − ω
0)
2 Ω sin Ω t
%
&
' (
) * exp i
2 ( ω − ω
0) t
+ , -
. / 0
€
C
2(t ) = − i γ
Ω sin Ωt exp − i
2 ( ω − ω
0) t
&
' ( )
* +
€
Ω = γ
2+ ( ω − ω
0)
24
€
ω
0€
ε
2€
ε
1€
C
2 2€
C
12Rotating wave approximation
Initial condition
4/19 No. 24
Rabi oscillation ラビ振動
€
Ω = γ
2+ ( ω − ω
0)
2€ 4
C
2(t )
2= γ
2Ω
2sin
2Ωt
€
C
1(t )
2= 1 − C
2(t )
2Population
ポピュレーション€
ω = ω
0€
ω − ω
0= 0.92 γ
€
γ t
€
γ t
€
γ t
€
C
1(t )
2€
C
2(t )
2€
ω − ω
0= 3.5 γ
吸収 放出 吸収 放出