量子力学 (5回目)

2023/5/27(土)
量子力学 (5回目)
 
(Quantum mechanics)
 
今回はトンネル効果です
 
 
■ 記号
t :時間(s)
x :位置(m)
x :位置ベクトル(m)
m :質量(kg)
E :エネルギー(J)
h :プランク定数(6.62607015×10-34J・s)
ℏ :ディラック定数 [ℏ = h / (2π)]
φ(x):時間を含まない波動関数
V(x):ポテンシャルエネルギー
H :ハミルトニアン
∇:ナブラ [∇ = (∂/∂x)]
d :壁の長さ(m)
V0:壁のポテンシャル [V0≧Eとする]
 
時間を含まない1次元のシュレディンガー方程式
[{-ℏ2/(2m)}(∂2/∂x2) + V(x)]φ(x) = Eφ(x)
 
 
 
■ 長さ有限の壁の解
長さdの有限の壁(壁のポテンシャルV0≧E)
 
▼ x < 0 のとき
V(x) = 0
{-ℏ2/(2m)}(∂2/∂x2)φ(x) = Eφ(x)
 
(∂2/∂x2)φ(x) = (-2mE/ℏ2)φ(x)
k = √(2mE/ℏ2) … (E ≧ 0)
と置くと
(∂2/∂x2)φ(x) = -k2φ(x)
φ(x) = Aexp(ikx) + Bexp(-ikx)
(代入すると成り立つ)
 
▼ 0 ≦ x ≦ dのとき
{ℏ2/(2m)}(∂2/∂x2)φ(x) - V(x)φ(x) = -Eφ(x)
(∂2/∂x2)φ(x) = (2m/ℏ2){V(x) - E}φ(x)
V(x) = V0 より
 
(∂2/∂x2)φ(x) = (2m/ℏ2)(V0 - E)φ(x)
λ = √{(2m/ℏ2)(V0 - E)} … (V0 ≧ E)
と置くと
(∂2/∂x2)φ(x) = λ2φ(x)
φ(x) = Cexp(λx) + Dexp(-λx)
(代入すると成り立つ)
 
▼ x > d のとき
V(x) = 0
{-ℏ2/(2m)}(∂2/∂x2)φ(x) = Eφ(x)
x < 0のとき同様にして
φ(x) = Fexp(ikx) + Gexp(-ikx)
 
 
■ 波の進行方向
正弦波の式
y = Acos(kx - ωt)
y = Asin(kx - ωt)
は右進行波を表しているので
 
波動関数
Ψ(x,t) = A{cos(kx - ωt) + isin(kx - ωt)}
= Aexp{i(kx - ωt)}

φ(x) = A{cos(kx) + isin(kx)} = Aexp(ikx)
も右進行波を表している
 
同様に
y = Asin{-(kx - ωt)}
は左進行波を表しているので
φ(x) = Aexp(-ikx)
も左進行波を表している
 
 
 
■ トンネル効果の透過確率T
φ(x < 0  ) = Aexp(ikx) + A'exp(-ikx)
φ(0≦x≦d) = Bexp(λx) + B'exp(-λx)
φ(x > d  ) = Cexp(ikx) + C'exp(-ikx)
について
波がx < 0の位置から右進行している場合
A exp( ikx) … 壁の左側の入射波(右進行)
A'exp(-ikx) … 壁の左側の反射波(左進行)
C exp( ikx) … 壁の右側の透過波(右進行)
となる
 
存在確率 |φ(x)|2 より
透過率Tを
T = |Cexp(ikx)|2 / |Aexp(ikx)|2 = |C/A|2 
として計算してみる
 
▼ 境界条件(連続性の確保)を考える
φ'(x) = dφ(x)/dx
φ'(x < 0  ) = ik{Aexp(ikx) - A'exp(-ikx)}
φ'(0≦x≦d) = λ{Bexp(λx) - B'exp(-λx)}
φ'(x > d  ) = ik{Cexp(ikx) - C'exp(-ikx)}
 
x = 0のとき、φ(0)、φ'(0)より
A + A' = B + B' … ①
ik(A - A') = λ(B - B') … ②
x = dのとき、φ(d)、φ'(d)より
Bexp(λd) + B'exp(-λd) = Cexp(ikd) … ③
λ{Bexp(λd) - B'exp(-λd)} = ikCexp(ikd) … ④
 
③±④/λの1/2倍より
B = C{(1 + ik/λ)/2}exp(ikd)/exp(λd)
= C{(λ + ik)/(2λ)}exp(ikd)/exp(λd)
B'= C{(1 - ik/λ)/2}exp(ikd)/exp(-λd)
= C{(λ - ik)/(2λ)}exp(ikd)/exp(-λd)
 
B,B'を①に代入して分子分母をik倍する
A + A' = C{exp(ikd)/(2λ)}
{(λ + ik)/exp(λd) + (λ - ik)/exp(-λd)
= C{exp(ikd)/(2λ)}
{(λ + ik)exp(-λd) + (λ - ik)exp(λd)}
= C{exp(ikd)/(2ikλ)}
{(ikλ - k2)exp(-λd) + (ikλ + k2)exp(λd)}
 
B,B'を②/ikに代入する
A - A' = (λ/ik)C{exp(ikd)/(2λ)}
{(λ + ik)exp(-λd) - (λ - ik)exp(λd)}
= C{exp(ikd)/(2ikλ)}
{(λ2 + ikλ)exp(-λd) - (λ2 - ikλ)exp(λd)}
 
{A + A' + (A - A')}/2 = A
= C{exp(ikd)/(4ikλ)}
{(ikλ - k2)exp(-λd) + (ikλ + k2)exp(λd)}
+ C{exp(ikd)/(4ikλ)}
{(λ2 + ikλ)exp(-λd) - (λ2 - ikλ)exp(λd)}
= C{exp(ikd)/(4ikλ)}
{(ikλ - k2)exp(-λd) + (ikλ + k2)exp(λd)
+ (λ2 + ikλ)exp(-λd) - (λ2 - ikλ)exp(λd)}
= C{exp(ikd)/(4ikλ)}
{(k2 + 2ikλ - λ2)exp(λd)
-(k2 - 2ikλ - λ2)exp(-λd)}
= C{exp(ikd)/(4ikλ)}
{(k + iλ)2exp(λd) - (k - iλ)2exp(-λd)}
 
C/A = 4ikλexp(-ikd)/
{(k + iλ)2exp(λd) - (k - iλ)2exp(-λd)}
 
 
 
▼ 透過率T
sinh(x) = (ex - e-x)/2
cosh(x) = (ex + e-x)/2
ex  = cosh(x) + sinh(x)
e-x = cosh(x) - sinh(x)
eix  = cos(x) + isin(x)
e-ix  = cos(-x) + isin(-x) = cos(x) - isin(x)
|e-ix|2 = √[{cos(x)}2 - {isin(x)}2]2 
= cos2(x) + sin2(x) = 1
sinh2(x) - cosh2(x) = 1
|z|2 = zz(
)
|a+ib|2 = (a+ib)(a-ib)
 
 
(A/C){4ikλexp(-ikd)}
= (k + iλ)2exp(λd) - (k - iλ)2exp(-λd)
= (k2 + 2ikλ - λ2)exp(λd)
+ (-k2 + 2ikλ + λ2)exp(-λd)
= (k22){exp(λd) - exp(-λd)}
+ 2ikλ{exp(λd) + exp(-λd)}
= 2(k22)sinh(λd) + 4ikλcosh(λd)
 
C/A
= {2ikλexp(-ikd)}
/{(k22)sinh(λd) + 2ikλcosh(λd)}
 
|C/A|2 
= |2ikλexp(-ikd)|2 
/|(k22)sinh(λd) + 2ikλcosh(λd)|2 
= |-4k2λ2|
/|(k22)sinh(λd) + 2ikλcosh(λd)|2 
= (4k2λ2)
/√[{(k22)sinh(λd)}2 - {2ikλcosh(λd)}2]2 
= (4k2λ2)
/{(k22)2sinh2(λd) + 4k2λ2cosh2(λd)}
= (4k2λ2)
/{(k4-2k2λ24)sinh2(λd) + 4k2λ2cosh2(λd)}
= (4k2λ2) / {(k4+2k2λ24)sinh2(λd)
+ 4k2λ2cosh2(λd) - 4k2λ2sinh2(λd)}
 
T = (4k2λ2)/{(k2 + λ2)2sinh2(λd) + 4k2λ2}
 
 
▼ k,λを元に戻す
k = √(2mE/ℏ2) … (E ≧ 0)
λ = √{(2m/ℏ2)(V0 - E)} … (V0 ≧ E)
 
4k2λ2 
= 4(2mE/ℏ2)(2m/ℏ2)(V0 - E)
= (16m2/ℏ4)E(V0 - E)
 
T = (4k2λ2)/{(k2 + λ2)2sinh2(λd) + 4k2λ2}
= (16m2/ℏ4)E(V0 - E)
/ [{(2mE/ℏ2) + (2m/ℏ2)(V0 - E)}2 
sinh2[d√{(2m/ℏ2)(V0 - E)}] + (16m2/ℏ4)E(V0 - E)
= (16m2/ℏ4)E(V0 - E) / (4m2/ℏ4){E + (V0 - E)}2 
sinh2[(d/ℏ)√{2m(V0 - E)}) + (16m2/ℏ4)E(V0 - E)
= (16m2/ℏ4)E(V0 - E) / (4m2/ℏ4)(E + V0 - E)2 
sinh2[(d/ℏ)√{2m(V0 - E)}) + (16m2/ℏ4)E(V0 - E)
 
T
= 4E(V0 - E)
/ [4E(V0 - E) + V02sinh2[(d/ℏ)√{2m(V0 - E)}]]
 
 
▼ E/V0の形にする
T
= 4E(V0 - E)
/ [4E(V0 - E) + V02sinh2[(d/ℏ)√{2m(V0 - E)}]]
= 4E/V02(V0 - E)
/ [4E/V02(V0 - E) + sinh2[(d/ℏ)√{2mV0(1 - E/V0)}]]
= 4E/V0(1 - E/V0)
/ [4E/V0(1 - E/V0) + sinh2[(d/ℏ)√{2mV0(1 - E/V0)}]]
 

▼ 結果(透過率T) 
r = E/V0 
T = 4r(1 - r)
/ [4r(1 - r) + sinh2[(d/ℏ)√{2mV0(1 - r)}]]
 
 
▼ トンネル効果の反射率R
反射率Rは|A'/A|2 または
R = 1 - T
で求められる
 
 

)

このブログの人気の投稿

NEWS

N88-BASICでゲーム (1回目)