콘텐츠로 이동

7부 고체물리 7.6

자유전자 모형

앞 다섯 장은 이온이 주인공이었다. 이온이 어디에 놓이고(7.1·7.2), 왜 붙어 있으며(7.3), 어떻게 흔들리는지(7.4·7.5)를 다뤘고 전자는 결합을 만드는 풀로만 등장했다. 이 장부터 전자가 주인공이 된다. 그리고 첫 모형은 놀랄 만큼 거칠다 — 이온을 통째로 지운다. 금속의 원자가 전자를 빈 상자 속의 기체로 보고, 이온은 전하 중성을 맞추는 균일한 양전하 배경으로만 남긴다. 전자끼리의 쿨롱 반발도 무시한다. 남는 것은 상자 속의 자유 입자뿐이다.

이런 모형이 맞을 리 없어 보이는데 — 이온의 쿨롱 퍼텐셜은 수 eV 깊이고 전자 사이 거리는 원자 간격과 같다 — 실제로는 알칼리 금속과 귀금속의 전기전도·열전도·열용량·홀 효과·광학 반사를 정량적으로 맞힌다. 왜 맞는지는 이 장이 아니라 뒤 장들이 답한다. 주기적 이온 퍼텐셜 속의 전자는 자유전자와 같은 꼴로 움직이되 질량만 m∗m^{*} 로 바뀐다는 것(7.7 블로흐 정리, 7.9 유효질량), 전자 기체가 이온 전하를 짧은 거리에서 가려 버린다는 것(7.10 차폐), 파울리 배타가 전자끼리의 산란을 거의 막는다는 것(7.10 페르미 액체)이 그 이유다. 이 장은 그 결과를 받아 쓰는 자리가 아니라, 그 결과가 필요해지는 자리다 — 자유전자 모형이 어디까지 맞고 어디서 틀리는지를 정확히 알아야 7.7 이후가 무엇을 고치는지 보인다.

장은 두 층으로 되어 있다. 1절의 드루데 모형(1900)은 전자를 고전 기체로 본다. 옴의 법칙과 홀 효과를 한 줄로 내고 비데만–프란츠 법칙까지 맞히지만, 열용량에서 100배 틀린다. 2절부터의 좀머펠트 모형(1927)은 운동 방정식은 그대로 두고 통계만 맥스웰–볼츠만에서 페르미–디랙으로 바꾼다. 그 한 번의 교체가 열용량(4절), 평균 자유 행정(5절), 열전도(8절)의 100배 오차를 전부 없앤다 — 3절의 페르미 구가 그 교체의 기하다. 6절의 비저항과 7절의 홀 효과는 두 모형이 같은 답을 주지만, 답의 해석은 다르다.

반도체에서 이 장이 어디로 가는지 미리 적어 둔다.

  • 캐리어 통계 (11.1) — 반도체의 전자 농도 n=∫D(E)f(E) dEn=\int D(E)f(E)\,dE 는 3절의 상태밀도와 2절의 페르미–디랙 분포를 m→m∗m\to m^{*}, E→E−EcE\to E-E_c 로 바꿔 쓴 것이다. 축퇴 반도체 (11.1.7)는 이 장의 금속 그 자체다.
  • 이동도와 전도도 (11.2) — σ=ne2τ/m∗\sigma=ne^2\tau/m^{*}, μn=eτ/m∗\mu_n=e\tau/m^{*}. 금속과 반도체의 전도도가 10810^{8} 배 다른 것은 τ\tau 가 아니라 nn 의 차이다 (5절).
  • 홀 측정 (11.2.9·19부) — 캐리어의 부호와 농도를 재는 유일한 직접 수단이 7절의 RH=−1/neR_H=-1/ne 다.
  • 금속 배선 (18.22·21.7) — Cu 의 평균 자유 행정 40 nm 가 배선 폭과 같아지면서 비저항이 2~3배 오른다. 6절의 마티센 규칙과 크기 효과가 그 계산이다.
  • 광학과 플라스몬 (7.10·7.11·11.8) — 5절의 σ(ω)\sigma(\omega) 와 플라스마 주파수 ωp\omega_p 가 금속의 반사, 도핑 반도체의 자유 캐리어 흡수, 투명 전극의 설계를 정한다.
  • 열 설계 (3.7.13·18부) — 8절의 비데만–프란츠 법칙은 배선의 전기 저항에서 열전도도를 바로 읽게 해 준다. 전기가 잘 통하는 금속이 열도 잘 통하는 이유가 여기 있다.

표기는 Kittel 을 따르되 에너지와 전기장의 글자는 이 사이트의 뒤 장들(7.9의 EE–kk 다이어그램, 11부)과 맞춘다. 두 교재 모두 전기장에 EE 를 쓰고 에너지에는 ε\varepsilon (Kittel) 또는 E\mathcal{E} (Ashcroft·Mermin)를 쓴다. 여기서는 반대로 에너지가 EE, 전기장이 E\mathcal{E} 다 — 반도체 교재(Pierret·Sze)의 관례이고, 밴드 다이어그램의 세로축을 EE 로 읽는 뒤 장들과 어긋나지 않기 위해서다. Kittel 의 εF\varepsilon_F 는 이 장의 EFE_F 다.

이 장은 대소문자만 다른 짝과 같은 글자의 다른 뜻이 유난히 많다. 장 전체에 걸리므로 여기 한 번 정리한다.

기호읽는 법무엇인가단위 · 비고
nn, NN전자 농도, 전자 수n=N/Vn=N/V. 소문자는 단위 부피당, 대문자는 전체m−3\mathrm{m^{-3}}, —
ee기본 전하양수 1.602×10−19 C1.602\times10^{-19}\ \mathrm{C}. 전자의 전하는 −e-eC
mm, m∗m^{*}전자 질량, 유효질량이 장은 전부 자유전자 질량 mm. m∗m^{*} 는 7.9에서kg
EE, E\mathcal{E}에너지, 전기장서체가 다르다. EFE_F 는 페르미 에너지eV, V/m
k\mathbf{k}, kBk_B파수벡터, 볼츠만 상수kFk_F 는 페르미 파수. kBk_B 는 언제나 첨자 Bm−1\mathrm{m^{-1}}, J/K
μ\mu, μn\mu_n화학퍼텐셜, 이동도μ(T=0)=EF\mu(T{=}0)=E_F. 이동도는 11부 표기대로 첨자 nneV, m2/V s\mathrm{m^2/V\,s}
τ\tau, ll완화 시간, 평균 자유 행정l=vFτl=v_F\tau (5절). 드루데는 vthτv_{\mathrm{th}}\tau 로 잘못 잡았다s, m
σ\sigma, ρ\rho전도도, 비저항ρ=1/σ\rho=1/\sigma. 질량 밀도는 ρm\rho_m 으로 구분S/m, Ω m\Omega\,\mathrm{m}
D(E)D(E)상태밀도단위 에너지당 상태 수 (스핀 포함). 부피당은 g(E)=D/Vg(E)=D/VeV−1\mathrm{eV^{-1}}
f(E)f(E)페르미–디랙 분포에너지 EE 인 상태가 차 있을 확률—
RHR_H홀 계수Ey/(jxB)\mathcal{E}_y/(j_xB). 자유전자는 −1/ne-1/nem3/C\mathrm{m^3/C}
L0L_0, LL로렌츠 수κ/σT\kappa/\sigma T. 첨자 0 은 이론값 2.44×10−82.44\times10^{-8}. 3절의 상자 변 LL 과 무관W Ω K−2\mathrm{W\,\Omega\,K^{-2}}
Θ\Theta디바이 온도7.5의 ΘD\Theta_D. 비저항 맞춤에는 ΘR\Theta_RK

금속은 이온 격자에 담긴 전자 기체다

섹션 제목: “금속은 이온 격자에 담긴 전자 기체다”

톰슨이 전자를 발견한 지 3년 뒤인 1900년, 드루데는 금속을 기체 운동론(3.8.8기체 운동론과 평균 자유 행정)으로 다뤘다. 원자가 전자 ZZ 개가 원자를 떠나 결정 전체를 돌아다니고, 남은 이온은 제자리에 박혀 있다. 전자는 기체 분자처럼 직선으로 날다가 무엇인가에 부딪혀 방향을 바꾼다. 가정을 넷으로 적으면 이렇다 (Ashcroft·Mermin Ch.1).

  1. 자유전자 근사. 충돌과 충돌 사이에 전자는 이온이나 다른 전자로부터 아무 힘도 받지 않는다. 외부 전기장·자기장만 작용한다.
  2. 독립전자 근사. 전자끼리의 상호작용을 무시한다.
  3. 충돌은 순간적이다. 어느 전자든 다음 dtdt 동안 충돌할 확률은 dt/τdt/\tau 로, 전자의 위치나 속도와 무관하다. τ\tau 를 완화 시간(relaxation time) 또는 충돌 시간이라 한다.
  4. 충돌이 열평형을 만든다. 충돌 직후의 속도는 충돌 전 속도와 무관하게, 그 자리의 온도에 맞는 열분포에서 무작위로 뽑힌다.

가정 1과 2는 명백히 거칠다 — 그럼에도 왜 되는지가 7.7과 7.10의 주제다. 가정 3과 4는 드루데가 “이온과의 충돌”로 그렸지만, 무엇과 부딪히는지는 이 모형 안에서 정할 수 없고 정할 필요도 없다. τ\tau 를 하나의 실험 상수로 두면 나머지가 다 나온다. 완벽한 주기 격자는 전자를 산란시키지 않으며(7.7), 실제로 부딪히는 상대는 격자의 진동(포논)과 결함이라는 사실(6절)은 30년 뒤에야 드러났다. 그때도 τ\tau 의 해석만 바뀌고 식은 그대로 살아남았다.

원자가 ZZ, 원자량 AA, 질량 밀도 ρm\rho_m 인 금속의 전도 전자 농도는

n=NV=NA ZρmAn=\frac{N}{V}=N_A\,\frac{Z\rho_m}{A}

이다 (NAN_A 아보가드로 수). 값은 1022∼1023 cm−310^{22}\sim10^{23}\ \mathrm{cm^{-3}} 로, 고전 기체(상온 상압 2.7×1019 cm−32.7\times10^{19}\ \mathrm{cm^{-3}})의 천 배다. 전자 하나가 차지하는 부피를 반지름 rsr_s 인 구로 나타내면

4π3rs3=1n⟹rs=(34πn)1/3\frac{4\pi}{3}r_s^3=\frac1n \qquad\Longrightarrow\qquad r_s=\Bigl(\frac{3}{4\pi n}\Bigr)^{1/3}

이고, 보어 반지름 a0=0.529a_0=0.529 Å 으로 나눈 rs/a0r_s/a_0 가 금속마다 2~6 이다. 곧 전자 사이 거리가 원자 크기와 같다 — 이온 퍼텐셜도, 전자끼리의 반발도 무시할 크기가 아니다. 가정 1·2가 왜 되는지 묻는 것이 정당한 이유다.

금속ZZnn (1022 cm−310^{22}\ \mathrm{cm^{-3}})rs/a0r_s/a_0ρ\rho (273 K, μΩ cm\mu\Omega\,\mathrm{cm})τ\tau (273 K, 10−1410^{-14} s)
Na12.653.934.23.2
K11.404.866.14.1
Cu18.472.671.562.7
Ag15.863.021.514.0
Au15.903.012.043.0
Al318.12.072.450.80

τ\tau 열은 아래의 σ=ne2τ/m\sigma=ne^2\tau/m 을 실측 ρ\rho 로 거꾸로 푼 값이다 (Ashcroft·Mermin Table 1.1–1.3). 10−1410^{-14} s — 이 자릿수를 기억해 둔다. 5절에서 이 시간에 페르미 속도를 곱하면 평균 자유 행정이 나온다.

운동 방정식 — 마찰이 있는 뉴턴 법칙

섹션 제목: “운동 방정식 — 마찰이 있는 뉴턴 법칙”

전자 하나의 운동량 p=mv\mathbf p=m\mathbf v 가 외력 F(t)\mathbf F(t) 아래 어떻게 변하는지 보자. 시각 tt 에 운동량 p(t)\mathbf p(t) 인 전자를 잡아 dtdt 뒤를 본다. 확률 dt/τdt/\tau 로 충돌했고, 충돌한 전자의 운동량은 가정 4에 따라 방향이 무작위라 평균 0 이다. 확률 1−dt/τ1-dt/\tau 로 충돌하지 않았고, 그 전자는 p+F dt\mathbf p+\mathbf F\,dt 를 가진다. 평균 운동량은

p(t+dt)=(1−dtτ)(p(t)+F dt)+dtτ⋅0\mathbf p(t+dt)=\Bigl(1-\frac{dt}{\tau}\Bigr)\bigl(\mathbf p(t)+\mathbf F\,dt\bigr) +\frac{dt}{\tau}\cdot 0

이고, dt2dt^2 을 버리면

dpdt=F(t)−pτ\frac{d\mathbf p}{dt}=\mathbf F(t)-\frac{\mathbf p}{\tau}

가 된다. 뉴턴 법칙에 마찰항 −p/τ-\mathbf p/\tau 가 붙은 꼴이다. 1.4의 감쇠 진동 (1.4.42계 선형 ODE와 상수 계수)과 같은 방정식이고, 외력을 끄면 p∝e−t/τ\mathbf p\propto e^{-t/\tau} 로 사라진다 — 그래서 τ\tau 를 완화 시간이라 부른다. 마찰항의 정체는 점성이 아니라 충돌이 방향의 기억을 지우는 속도다.

전기장 E\mathcal{E} 를 걸면 F=−eE\mathbf F=-e\boldsymbol{\mathcal E} 다. 정상 상태에서 dp/dt=0d\mathbf p/dt=0 이므로

p=−eτE,vd=pm=−eτmE\mathbf p=-e\tau\boldsymbol{\mathcal E}, \qquad \mathbf v_d=\frac{\mathbf p}{m}=-\frac{e\tau}{m}\boldsymbol{\mathcal E}

— 이 평균 속도가 표류 속도(drift velocity)다. 전류 밀도는 전자 농도 × 전하 × 속도이므로

j=−ne vd=ne2τmE≡σE,  σ=ne2τm  \mathbf j=-ne\,\mathbf v_d=\frac{ne^2\tau}{m}\boldsymbol{\mathcal E} \equiv\sigma\boldsymbol{\mathcal E}, \qquad \boxed{\;\sigma=\frac{ne^2\tau}{m}\;}

이다. 이것이 옴의 법칙(2.6.1옴의 법칙과 전도)이고, 드루데 모형이 처음으로 준 전도도의 미시적 표현이다. 전기장에 비례하는 전류 — 비례 상수에 전자 농도, 전하의 제곱, 완화 시간이 들어가고 질량이 분모에 있다. 전하가 제곱인 이유는 두 번 쓰이기 때문이다: 한 번은 힘을 받을 때(−eE-e\mathcal E), 한 번은 전류를 나를 때(−ev-e\mathbf v). 그래서 전하의 부호는 전도도에 나타나지 않는다 — 정공도 전자와 같은 방향의 전류를 만든다.

충돌 사이의 자유 비행 — 전기장은 매 비행을 조금씩 휘어 표류를 만든다 전기장 없음 — 제자리 근처를 맴돈다 출발 전기장 𝓔 (+x) — 같은 충돌열이 왼쪽으로 밀린다 출발 𝓔 알짜 변위 ∝ vd​ t 옅은 선: 전기장 없을 때의 같은 충돌열 ○ 충돌 — 직후 속도는 이전 속도와 무관하게 다시 뽑힌다. 표류 속도 vd​ = −e𝓔τ/m 는 열속도보다 훨씬 작다.
그림 7.6-1. 충돌 사이의 자유 비행과 표류. 같은 충돌열(출발 방향과 비행 시간이 같은 열여섯 번의 비행)을 두 번 그렸다. 왼쪽은 전기장이 없을 때 — 충돌마다 방향이 초기화되어 출발점 근처를 맴돈다. 오른쪽은 +x 방향 전기장을 걸었을 때 — 전자가 받는 힘은 −e𝓔 이므로 매 비행이 −x 쪽으로 휘고, 같은 충돌열인데도 전체가 왼쪽으로 밀린다(옅은 선이 왼쪽 궤적). 표류 속도 v_d = −e𝓔τ/m 는 한 비행의 속력보다 훨씬 작지만, 방향이 있어 평균에서 살아남는다. 가속은 눈에 보이도록 과장했다.직접 작도

드루데의 평균 자유 행정 — 첫 번째 균열

섹션 제목: “드루데의 평균 자유 행정 — 첫 번째 균열”

충돌 사이에 전자가 나는 거리, 평균 자유 행정 l=vτl=v\tau 를 어림해 보자. 드루데는 전자를 고전 기체로 보았으니 속력에 등분배 정리의 열속도 vth=3kBT/m≈1.1×105 m/sv_{\mathrm{th}}=\sqrt{3k_BT/m}\approx1.1\times10^{5}\ \mathrm{m/s} (273 K) 를 넣었다. τ≈3×10−14\tau\approx3\times10^{-14} s 와 곱하면 l≈3l\approx3 nm, 곧 원자 간격의 열 배 안팎이다. “이온에 부딪히며 날아가는 전자”라는 그림과 그럴듯하게 맞았고, 드루데는 만족했다.

이것이 틀렸다는 것은 두 갈래로 드러난다. 첫째, 저온 순수 금속의 τ\tau 는 실온의 수천 배로 커져 ll 이 밀리미터에 이른다 — 이온이 산란체라면 있을 수 없는 길이다. 둘째, 5절에서 보듯 전자의 속력은 vthv_{\mathrm{th}} 가 아니라 그보다 열 배 큰 페르미 속도 vF≈1.6×106v_F\approx1.6\times10^{6} m/s 라서, 실온에서도 l≈40l\approx40 nm 다. 이 40 nm 가 반도체 배선 설계의 기준 길이가 된다(6절).

성적표 — 무엇을 맞히고 무엇을 틀리나

섹션 제목: “성적표 — 무엇을 맞히고 무엇을 틀리나”

드루데 모형은 한 세대 동안 금속의 표준 이론이었고, 지금도 반도체 소자 방정식의 뼈대다 (σ=neμ\sigma=ne\mu, 홀 계수, 플라스마 주파수는 전부 드루데 식 그대로다). 성적을 정직하게 적는다.

현상드루데의 예측실측판정
옴의 법칙σ=ne2τ/m\sigma=ne^2\tau/mτ∼10−14\tau\sim10^{-14} s 로 맞춤○ (τ\tau 는 실험 상수)
홀 계수 (7절)RH=−1/neR_H=-1/ne알칼리 금속 ○, Al·Be 는 부호 반대△
비데만–프란츠 (8절)κ/σT=\kappa/\sigma T= 상수상수 맞음, 값은 반△ (두 오차가 상쇄)
전자 열용량 (4절)32nkB\tfrac32 nk_B예측의 1 %✕
열전 계수 (8절)∼100 μ\sim100\ \muV/K∼1 μ\sim1\ \muV/K✕ (100배)
평균 자유 행정원자 간격 수준저온에서 mm✕
ρ(T)\rho(T)설명 없음∝T\propto T (고온), ∝T5\propto T^5 (저온)—

열용량과 열전 계수의 100배는 같은 뿌리다 — 전자를 고전 기체로 본 것. 2절에서 통계를 바꾸면 이 둘이 한꺼번에 고쳐지고 비데만–프란츠의 값도 맞는다. 홀 계수의 부호와 ρ(T)\rho(T) 는 통계로는 안 되고 밴드 이론(7.7·7.9)과 포논 산란(6절)이 필요하다.

전자는 페르미온이다 — 바꾸는 것은 통계 하나뿐

섹션 제목: “전자는 페르미온이다 — 바꾸는 것은 통계 하나뿐”

드루데가 놓친 것은 역학이 아니라 셈법이다. 전자는 스핀 12\tfrac12 인 페르미온이라 파울리 배타원리(5.6.5파울리 배타원리)를 따르고, 한 양자 상태에는 (스핀까지 세어) 전자 하나만 들어간다. 그러면 열평형에서 에너지 EE 인 상태가 차 있을 확률은 맥스웰–볼츠만 인자 e−E/kBTe^{-E/k_BT} 가 아니라 페르미–디랙 분포(6.4.7페르미-디랙 분포의 유도)

f(E)=1e(E−μ)/kBT+1f(E)=\frac{1}{e^{(E-\mu)/k_BT}+1}

이다. μ\mu 는 화학퍼텐셜 — 전자를 하나 더 넣을 때 드는 자유에너지이고, 실제로는 “전자 수가 NN 이 되도록” 정해지는 규격화 상수다. 절대영도에서의 값을 페르미 에너지 EF≡μ(T=0)E_F\equiv\mu(T{=}0) 라 하고, 그것을 온도로 바꾼 TF=EF/kBT_F=E_F/k_B 를 페르미 온도라 한다.

좀머펠트(1927)는 드루데의 가정 1~3을 그대로 두고 가정 4의 “열분포”를 이 분포로 바꿨다. 그것이 전부다. 운동 방정식도, τ\tau 도, σ=ne2τ/m\sigma=ne^2\tau/m 도 살아남는다. 그런데 열용량· 열전도·평균 자유 행정의 100배 오차가 전부 사라진다. 어떻게 그럴 수 있는지가 이 절과 다음 두 절이다.

원리를 보는 데는 1차원이면 충분하다. 길이 LL 인 선 위에 전자 NN 개를 놓고 양 끝을 무한 벽으로 막는다. 허용 상태는 무한 우물의 정상파(5.3.1무한 우물)

ψn(x)=Asin⁡nπxL,En=ℏ22m(nπL)2,n=1,2,3,…\psi_n(x)=A\sin\frac{n\pi x}{L}, \qquad E_n=\frac{\hbar^2}{2m}\Bigl(\frac{n\pi}{L}\Bigr)^2, \qquad n=1,2,3,\dots

이다. 절대영도에서 전자들은 파울리 배타 때문에 한 준위에 둘씩(스핀 ↑↓) 아래부터 채운다. NN 개를 다 넣으면 마지막으로 찬 준위의 양자수는 nF=N/2n_F=N/2 이고, 그 에너지가 페르미 에너지다:

EF=ℏ22m(nFπL)2=ℏ22m(Nπ2L)2E_F=\frac{\hbar^2}{2m}\Bigl(\frac{n_F\pi}{L}\Bigr)^2 =\frac{\hbar^2}{2m}\Bigl(\frac{N\pi}{2L}\Bigr)^2

여기서 이미 핵심이 보인다. EFE_F 는 온도와 무관하게 전자 밀도 N/LN/L 만으로 정해지고, 밀도가 금속 수준이면 그 값이 수 eV 다. 상온의 kBT=0.026k_BT=0.026 eV 는 그 백분의 일이다. 전자 기체는 실온에서도 사실상 절대영도에 있다.

온도가 하는 일 — 계단의 모서리만 둥글어진다

섹션 제목: “온도가 하는 일 — 계단의 모서리만 둥글어진다”

T=0T=0 에서 f(E)f(E) 는 EFE_F 아래에서 1, 위에서 0 인 계단이다. 온도를 올리면 계단이 어떻게 되는지 식에서 읽는다. E=μE=\mu 에서 f=12f=\tfrac12 — 언제나. E−μ≫kBTE-\mu\gg k_BT 면 f≈e−(E−μ)/kBTf\approx e^{-(E-\mu)/k_BT} — 볼츠만 꼬리. μ−E≫kBT\mu-E\gg k_BT 면 1−f≈e−(μ−E)/kBT1-f\approx e^{-(\mu-E)/k_BT} — 아래쪽도 대칭으로 빈다. 곧 온도는 μ\mu 를 중심으로 폭 몇 kBTk_BT 의 띠 안에서만 점유를 흐리고, 그 바깥은 T=0T=0 과 구별되지 않는다.

페르미–디랙 분포 — 온도는 EF​ 주변 몇 kB​T 폭만 건드린다 f(E) E 1 ½ EF​ T = 0 kB​T = 0.02 EF​ kB​T = 0.1 EF​ −∂f/∂E E EF​ ≈ 3.5 kB​T 이 창 안의 상태만이 열용량·전도·열전에 참여한다 미분
그림 7.6-2. 페르미–디랙 분포. 왼쪽은 f(E). T = 0 에서 E_F 에서 뚝 떨어지는 계단이고, 온도가 오르면 E_F 주변 몇 k_BT 폭만 흐려진다 — 실온의 구리(k_BT/E_F ≈ 0.004)는 그림의 0.02 곡선보다도 훨씬 계단에 가깝다. f(E_F) = ½ 은 온도와 무관하다. 오른쪽은 −∂f/∂E — E_F 를 중심으로 반높이 폭 3.53 k_BT 인 봉우리로, 넓이가 1 이다. 4절의 열용량, 5절의 전도도, 8절의 열전 계수는 전부 이 창 안의 상태만 본다: 창 바깥의 전자는 갈 곳이 없어(아래) 또는 있지도 않아(위) 아무것에도 응답하지 못한다.직접 작도

이 그림이 드루데의 100배를 설명한다. 온도를 ΔT\Delta T 만큼 올릴 때 에너지를 받을 수 있는 전자는 창 안의 전자, 곧 전체의 ∼T/TF\sim T/T_F 뿐이다. 구리는 TF=8.2×104T_F=8.2\times10^{4} K 라 실온에서 T/TF≈0.004T/T_F\approx0.004 — 전자 250개 중 하나만 온도를 “느낀다”. 나머지는 파울리 배타에 갇혀 있다: 위로 가려 해도 그 자리가 이미 차 있다. 그래서 열용량이 고전값의 1 % 다 (4절에서 정확히 센다).

한 식으로 금속과 반도체를 잇는다

섹션 제목: “한 식으로 금속과 반도체를 잇는다”

f(E)f(E) 의 볼츠만 꼬리를 다시 본다. E−μ≫kBTE-\mu\gg k_BT 인 상태의 점유는 f≈e−(E−μ)/kBTf\approx e^{-(E-\mu)/k_BT} 로, 맥스웰–볼츠만 분포와 같다. 곧 페르미–디랙 분포는 상태가 드문드문 찰 때 고전 통계로 돌아간다. 전도대 바닥이 μ\mu 보다 몇 kBTk_BT 위에 있는 보통의 반도체가 정확히 이 경우라, 11.1의 캐리어 농도 n=Nc e−(Ec−EF)/kBTn=N_c\,e^{-(E_c-E_F)/k_BT} 는 이 꼬리다. 금속은 μ\mu 가 밴드 깊숙이 들어가 계단 쪽을, 반도체는 꼬리 쪽을 쓴다 — 같은 함수의 두 끝이고, 축퇴 반도체가 그 사이를 잇는다.

한 변 LL 인 정육면체 상자(V=L3V=L^3) 속의 자유전자를 본다. 슈뢰딩거 방정식은 −ℏ22m∇2ψ=Eψ-\frac{\hbar^2}{2m}\nabla^2\psi=E\psi 이고, 해는 평면파(5.3.4자유 입자와 파속)

ψk(r)=1V eik⋅r,Ek=ℏ2k22m,p=ℏk,v=ℏkm\psi_{\mathbf k}(\mathbf r)=\frac{1}{\sqrt V}\,e^{i\mathbf k\cdot\mathbf r}, \qquad E_{\mathbf k}=\frac{\hbar^2k^2}{2m}, \qquad \mathbf p=\hbar\mathbf k,\quad \mathbf v=\frac{\hbar\mathbf k}{m}

이다. 1차원에서는 벽에서 ψ=0\psi=0 인 정상파를 썼지만, 3차원에서는 주기 경계 조건 ψ(x+L,y,z)=ψ(x,y,z)\psi(x+L,y,z)=\psi(x,y,z) 를 쓴다. 이유는 둘이다. 부피 성질은 표면이 어떻든 같아야 하므로 편한 조건을 골라도 되고, 정상파는 ⟨p⟩=0\langle\mathbf p\rangle=0 이라 전류를 나르지 못하는 반면 평면파 eik⋅re^{i\mathbf k\cdot\mathbf r} 는 운동량 ℏk\hbar\mathbf k 의 고유 상태라 5절의 전도를 다루기에 맞다. 주기 조건은 eikxL=1e^{ik_xL}=1 을 요구하므로

kx, ky, kz=0, ±2πL, ±4πL, …k_x,\,k_y,\,k_z=0,\ \pm\frac{2\pi}{L},\ \pm\frac{4\pi}{L},\ \dots

— 허용되는 k\mathbf k 는 kk 공간에서 간격 2π/L2\pi/L 인 격자점이다. 점 하나가 차지하는 kk 공간의 부피는 (2π/L)3=8π3/V(2\pi/L)^3=8\pi^3/V 이고, 한 점에 스핀 둘이 들어간다. 이 한 문장이 3절 전체다 — 상태를 세는 일이 kk 공간의 부피를 재는 기하로 바뀌었다.

페르미 구 — 전자 밀도가 모든 것을 정한다

섹션 제목: “페르미 구 — 전자 밀도가 모든 것을 정한다”

절대영도에서 전자 NN 개는 kk 공간의 원점에서 가까운 점부터 채운다. 에너지가 kk 의 크기만으로 정해지므로 채워진 영역은 구이고, 그 반지름을 페르미 파수 kFk_F, 표면을 페르미 면이라 한다. 전자 수를 세면

N=2×43πkF38π3/V=V3π2kF3⟹  kF=(3π2n)1/3  N=2\times\frac{\tfrac43\pi k_F^3}{8\pi^3/V}=\frac{V}{3\pi^2}k_F^3 \qquad\Longrightarrow\qquad \boxed{\;k_F=(3\pi^2n)^{1/3}\;}

이고, 나머지는 따라온다:

EF=ℏ22m(3π2n)2/3,vF=ℏkFm,TF=EFkBE_F=\frac{\hbar^2}{2m}(3\pi^2n)^{2/3}, \qquad v_F=\frac{\hbar k_F}{m}, \qquad T_F=\frac{E_F}{k_B}

전자 농도 하나가 넷을 전부 정한다. 온도도, 시료 크기도, 이온의 종류도 들어가지 않는다.

k 공간에서는 구, 에너지에서는 √E — 같은 상태를 두 번 센다 k 공간 (단면) kx​ ky​ kF​ 2π/L 페르미 면 ● 채워짐 (스핀 2) ○ 비어 있음 D(E) — 단위 에너지당 상태 수 E EF​ D(EF​) D(E) ∝ √E 채워진 상태 D(E) f(E) ~kB​T 폭만 흐려진다 면적 = N
그림 7.6-3. k 공간의 페르미 구와 에너지 공간의 상태밀도. 왼쪽: 주기 경계 조건이 허용하는 k 점은 간격 2π/L 의 격자이고, T = 0 에서 전자는 원점에서 가까운 점부터 스핀 둘씩 채우므로 채워진 영역은 반지름 k_F 인 구(그림은 단면)가 된다. N = 2 × (구 부피)/(점 하나의 부피)가 k_F = (3π²n)^{1/3} 을 준다. 오른쪽: 같은 상태를 에너지로 세면 D(E) ∝ √E 이고, 유한 온도에서 채워진 상태 D(E)f(E) 는 E_F 근처 k_BT 폭에서만 흐려진다. 채워진 영역의 면적이 N 이다.직접 작도
금속nn (1022 cm−310^{22}\ \mathrm{cm^{-3}})kFk_F (Å⁻¹)vFv_F (10610^{6} m/s)EFE_F (eV)TFT_F (10410^{4} K)
Na2.650.921.073.243.77
K1.400.750.862.122.46
Cu8.471.361.577.008.16
Ag5.861.201.395.496.38
Au5.901.211.405.536.42
Al18.11.752.0311.713.6

(Kittel Table 6.1.) 눈금을 익혀 둔다. kF∼1k_F\sim1 Å⁻¹ — 페르미 면 위 전자의 드브로이 파장 2π/kF2\pi/k_F 가 원자 간격과 같다. 우연이 아니라 n∼1/a3n\sim1/a^3 이기 때문이고, 그래서 페르미 면이 브릴루앙 영역 경계에 닿을 수 있다 — 7.7의 밴드 갭이 거기서 열린다. vF∼106v_F\sim10^{6} m/s — 빛의 0.5 %, 열속도의 열 배. EF∼E_F\sim 수 eV — 화학 결합 에너지와 같은 자릿수이고 kBTk_BT 의 백 배.

상태밀도 — 제곱근 곡선의 모양

섹션 제목: “상태밀도 — 제곱근 곡선의 모양”

에너지 EE 아래의 상태 수는 위의 셈에서 kF→k=2mE/ℏk_F\to k=\sqrt{2mE}/\hbar 로 바꾸면

N(E)=V3π2(2mEℏ2)3/2N(E)=\frac{V}{3\pi^2}\Bigl(\frac{2mE}{\hbar^2}\Bigr)^{3/2}

이고, 상태밀도(단위 에너지당 상태 수, 스핀 포함)는 그 미분이다:

D(E)=dNdE=V2π2(2mℏ2)3/2E1/2D(E)=\frac{dN}{dE}=\frac{V}{2\pi^2}\Bigl(\frac{2m}{\hbar^2}\Bigr)^{3/2}E^{1/2}

ln⁡N=32ln⁡E+\ln N=\tfrac32\ln E+ 상수 를 미분하면 더 짧다 — dN/N=32 dE/EdN/N=\tfrac32\,dE/E 이므로

D(E)=3N(E)2E,D(EF)=3N2EFD(E)=\frac{3N(E)}{2E}, \qquad D(E_F)=\frac{3N}{2E_F}

페르미 준위의 상태밀도는 전자 수를 페르미 에너지로 나눈 것의 32\tfrac32 배다. 이 한 줄이 4절의 열용량과 8절의 열전도를 계산 없이 준다. 형태는 E\sqrt E — 원점에서 수직으로 출발해 점점 완만해지는 곡선이다. 이유는 kk 공간 구 껍질의 넓이가 k2∝Ek^2\propto E 이고 dk/dE∝E−1/2dk/dE\propto E^{-1/2} 라서, 곱이 E1/2E^{1/2} 이다.

같은 셈을 차원만 바꿔 하면(6.5.9차원별 상태밀도) 2차원은 껍질이 원둘레 ∝k\propto k 라 D2D=mA/πℏ2D_{2D}=mA/\pi\hbar^2 로 에너지와 무관한 상수이고, 1차원은 D1D∝E−1/2D_{1D}\propto E^{-1/2} 로 바닥에서 발산한다. 3차원의 E\sqrt E, 2차원의 계단, 1차원의 뾰족한 봉우리 — 양자 우물· 나노선·양자점의 상태밀도가 이 순서로 나오고(10부), MOSFET 반전층의 부밴드(13부)가 2차원 계단이다.

절대영도의 총 에너지는 EFE_F 까지 E D(E)E\,D(E) 를 적분한 것이다. D∝E1/2D\propto E^{1/2} 이므로

U0=∫0EFE D(E) dE=3N2EF3/2∫0EFE3/2 dE=35 NEFU_0=\int_0^{E_F}E\,D(E)\,dE =\frac{3N}{2E_F^{3/2}}\int_0^{E_F}E^{3/2}\,dE =\frac35\,NE_F

— 전자 하나의 평균 에너지가 35EF\tfrac35E_F, 수 eV 다. 고전 기체의 32kBT\tfrac32k_BT 와 견주면 백 배이고, 절대영도인데도 0 이 아니다. 파울리 배타가 전자를 높은 상태로 밀어 올린 결과다.

이 에너지는 부피에 의존한다(EF∝V−2/3E_F\propto V^{-2/3}). 부피를 줄이면 U0U_0 이 오르므로 기체는 바깥으로 미는 압력을 낸다:

P=−∂U0∂V=23U0V=25 nEF,B=−V∂P∂V=53P=23 nEFP=-\frac{\partial U_0}{\partial V}=\frac23\frac{U_0}{V}=\frac25\,nE_F, \qquad B=-V\frac{\partial P}{\partial V}=\frac53P=\frac23\,nE_F

BB 는 부피 탄성률(7.3.9)이다. 구리에 넣으면 P≈38P\approx38 GPa, 곧 4×1054\times10^{5} 기압이고, B≈63B\approx63 GPa 다. 실측 부피 탄성률은 Na 6.4, K 3.1, Cu 137, Al 76 GPa 이고 이 식은 Na 9.2, K 3.2, Cu 63, Al 226 GPa 를 준다. 알칼리 금속에서는 자릿수와 경향이 맞는다 — 금속을 누르기 어려운 이유의 절반은 전자 기체의 축퇴 압력이다. Cu·Al 에서 어긋나는 것은 d 전자와 이온 반발이 함께 들기 때문이고, 그것은 7.3의 결합 이론이 맡는다.

등분배 정리(6.3.5등분배 정리와 고체의 열용량)로 자유 입자 NN 개의 열용량은 C=32NkBC=\tfrac32Nk_B 다. 원자가 전자가 하나씩인 금속이라면 이온의 3NkB3Nk_B(뒬롱–프티)에 전자의 32NkB\tfrac32Nk_B 가 얹혀 실온 열용량이 4.5NkB4.5Nk_B 여야 한다. 실측은 3NkB3Nk_B 에서 거의 벗어나지 않는다 — 전자의 기여는 예측의 1 % 다. 드루데 시대 최대의 수수께끼였고, 2절의 그림이 답을 반쯤 주었다: 온도를 느끼는 전자는 EFE_F 근처 kBTk_BT 폭 안의 전자, 곧 전체의 T/TFT/T_F 뿐이다.

정성적 셈 — 답의 꼴을 먼저 본다

섹션 제목: “정성적 셈 — 답의 꼴을 먼저 본다”

창 안의 전자 수를 Neff≈N (T/TF)N_{\mathrm{eff}}\approx N\,(T/T_F) 로 잡고, 이 전자들이 온도 TT 에서 각각 ∼kBT\sim k_BT 의 열에너지를 얻는다고 하면

U(T)−U0≈NTTF kBT⟹Cel=∂U∂T≈2NkBTTFU(T)-U_0\approx N\frac{T}{T_F}\,k_BT \qquad\Longrightarrow\qquad C_{\mathrm{el}}=\frac{\partial U}{\partial T}\approx 2Nk_B\frac{T}{T_F}

이다. 계수는 믿을 것이 못 되지만 두 가지는 이미 맞다. 온도에 비례한다 — 고전의 상수가 아니라. 그리고 크기가 NkBNk_B 에 T/TF∼10−2T/T_F\sim10^{-2} 를 곱한 것이다 — 100배의 차이가 여기서 난다. 정확한 계수는 π2/2\pi^2/2 이고, 아래에서 센다.

EFE_F 를 기준으로 재면, 온도 TT 에서 전자 기체가 더 가진 에너지는

ΔU=∫0∞E D(E)f(E) dE−∫0EFE D(E) dE\Delta U=\int_0^\infty E\,D(E)f(E)\,dE-\int_0^{E_F}E\,D(E)\,dE

이다. 전자 수가 보존되므로 ∫0∞Df dE=∫0EFD dE=N\int_0^\infty D f\,dE=\int_0^{E_F}D\,dE=N 이고, 양변에 EFE_F 를 곱해 위 식에서 빼면 (Kittel §6, 식 6.30–6.36)

ΔU=∫EF∞(E−EF) D(E) f(E) dE+∫0EF(EF−E) D(E) [1−f(E)] dE\Delta U=\int_{E_F}^{\infty}(E-E_F)\,D(E)\,f(E)\,dE +\int_0^{E_F}(E_F-E)\,D(E)\,\bigl[1-f(E)\bigl]\,dE

가 된다. 두 항 모두 양수이고 뜻이 분명하다 — 첫째는 EFE_F 위로 올라간 전자들이 얻은 에너지, 둘째는 EFE_F 아래에 생긴 구멍을 만드는 데 든 에너지다. 온도로 미분하면 D(E)D(E) 는 온도와 무관하니 ff 만 미분되고, 두 적분이 합쳐져

Cel=d ΔUdT=∫0∞(E−EF) D(E) ∂f∂T dEC_{\mathrm{el}}=\frac{d\,\Delta U}{dT}=\int_0^\infty(E-E_F)\,D(E)\,\frac{\partial f}{\partial T}\,dE

이다. ∂f/∂T\partial f/\partial T 는 EFE_F 근처 kBTk_BT 폭에서만 0 이 아니므로 D(E)≈D(EF)D(E)\approx D(E_F) 로 꺼낼 수 있다. μ≈EF\mu\approx E_F 로 두고 x=(E−EF)/kBTx=(E-E_F)/k_BT 를 쓰면

∂f∂T=E−EFkBT2 ex(ex+1)2\frac{\partial f}{\partial T}=\frac{E-E_F}{k_BT^2}\,\frac{e^x}{(e^x+1)^2}

이고, 적분의 아래끝 −EF/kBT≈−300-E_F/k_BT\approx-300 을 −∞-\infty 로 보내면

Cel=D(EF) kB2T∫−∞∞x2ex(ex+1)2 dx=π23 D(EF) kB2 TC_{\mathrm{el}}=D(E_F)\,k_B^2T\int_{-\infty}^{\infty}\frac{x^2e^x}{(e^x+1)^2}\,dx =\frac{\pi^2}{3}\,D(E_F)\,k_B^2\,T

— 적분값이 π2/3\pi^2/3 이다. 3절의 D(EF)=3N/2EFD(E_F)=3N/2E_F 를 넣으면

  Cel=π22 NkB TTF  \boxed{\;C_{\mathrm{el}}=\frac{\pi^2}{2}\,Nk_B\,\frac{T}{T_F}\;}

정성적 셈의 계수 2 가 π2/2≈4.9\pi^2/2\approx4.9 로 바뀌었을 뿐 모든 것이 같다. 구리 실온에서 T/TF=0.0037T/T_F=0.0037 → Cel=0.018 NkBC_{\mathrm{el}}=0.018\,Nk_B, 고전값 1.5 NkB1.5\,Nk_B 의 1.2 %. 100배가 정확히 나왔다.

계산에서 μ≈EF\mu\approx E_F 로 둔 것이 마음에 걸리면 좀머펠트 전개로 확인한다. H(E)H(E) 가 EFE_F 근처에서 매끄러우면

∫0∞H(E)f(E) dE=∫0μH(E) dE+π26(kBT)2H′(μ)+O[(kBT)4]\int_0^\infty H(E)f(E)\,dE =\int_0^{\mu}H(E)\,dE+\frac{\pi^2}{6}(k_BT)^2H'(\mu)+O\bigl[(k_BT)^4\bigr]

이다 (Ashcroft·Mermin Ch.2, 부록 C). 첫 항은 T=0T=0 의 계단, 둘째 항이 창의 분산 × 기울기다. H=DH=D 로 놓고 NN 을 고정하면 μ\mu 가 온도에 따라 어떻게 움직이는지 나온다:

μ(T)=EF[1−π212(TTF)2]\mu(T)=E_F\Bigl[1-\frac{\pi^2}{12}\Bigl(\frac{T}{T_F}\Bigr)^2\Bigr]

3차원에서는 DD 가 오르막이라 μ\mu 가 살짝 내려간다 — 창의 위쪽에 상태가 더 많으니 같은 전자 수를 맞추려면 중심을 낮춰야 한다. 구리 실온에서 10−510^{-5} 의 보정이고, 열용량 계산에는 (T/TF)2(T/T_F)^2 차수라 들어오지 않는다. 2차원에서는 DD 가 상수라 이 보정이 0 이다 — 반전층의 μ\mu 가 온도에 둔감한 이유다.

실측과 비교 — C/T 를 T² 에 대해 그린다

섹션 제목: “실측과 비교 — C/T 를 T² 에 대해 그린다”

저온 금속의 열용량은 전자 항과 포논 항(7.5.4디바이 모형과 T³ 법칙)의 합이다:

C=γT+AT3,γ=π23D(EF)kB2=π2NkB22EFC=\gamma T+AT^3, \qquad \gamma=\frac{\pi^2}{3}D(E_F)k_B^2=\frac{\pi^2Nk_B^2}{2E_F}

TT 에 비례하는 항과 T3T^3 에 비례하는 항을 가르는 표준 요령이 양변을 TT 로 나누어 T2T^2 에 대해 그리는 것이다: C/T=γ+AT2C/T=\gamma+AT^2 는 직선이고, 절편이 γ\gamma, 기울기가 AA 다.

C/T 를 T² 에 대해 그리면 절편이 전자, 기울기가 포논이다 C/T (mJ mol⁻¹ K⁻²) T² (K²) 0.1 0.2 0.3 2.0 2.4 2.8 γ = 2.08 전자 (절편) A = 2.57 포논 (기울기) 칼륨(K): C/T = γ + AT² C (mJ mol⁻¹ K⁻¹) T (K) 0.5 1 1.5 2 10 20 T* = √(γ/A) ≈ 0.9 K 전자 γT 포논 AT³ 합
그림 7.6-4. 금속의 저온 열용량. 칼륨의 값 γ = 2.08 mJ mol⁻¹ K⁻², A = 2.57 mJ mol⁻¹ K⁻⁴ (Kittel Fig. 6.9). 왼쪽: C/T 를 T² 에 대해 그리면 직선이고, T² → 0 절편이 전자 항 γ, 기울기가 포논 항 A 다 — 두 기여가 한 그래프에서 분리된다. 오른쪽: 같은 두 항을 T 에 대해 그리면 T* = √(γ/A) ≈ 0.9 K 아래에서 전자가, 위에서 포논이 이긴다. 실온에서는 포논이 전자의 백 배라 전자 항이 보이지 않고, 1 K 아래에서는 전자만 남는다.직접 작도
금속γobs\gamma_{\mathrm{obs}} (mJ mol⁻¹ K⁻²)γfree\gamma_{\mathrm{free}}mth/m=γobs/γfreem_{\mathrm{th}}/m=\gamma_{\mathrm{obs}}/\gamma_{\mathrm{free}}
Na1.381.091.26
K2.081.671.25
Cu0.6950.5051.38
Ag0.6460.6451.00
Au0.7290.6421.14
Al1.350.9121.48

(Kittel Table 6.2.) 자유전자 값은 표의 nn 만으로 계산한 것이다. 단순 금속에서 실측이 자유전자 값의 11.5배 안에 든다 — TT 에 비례한다는 사실과 크기의 자릿수를 모두 맞혔다. 전이 금속(Fe·Mn·Ni)은 515배로 크게 벗어나는데, 좁은 d 밴드가 EFE_F 에 걸려 D(EF)D(E_F) 가 큰 탓이다(7.13).

비율 γobs/γfree\gamma_{\mathrm{obs}}/\gamma_{\mathrm{free}} 를 열 유효질량 mth/mm_{\mathrm{th}}/m 이라 부른다. γ∝D(EF)∝m\gamma\propto D(E_F)\propto m 이므로(EF∝1/mE_F\propto1/m), 실측 γ\gamma 가 큰 것을 “전자가 무거운 것”으로 읽은 것이다. 벗어남의 원인은 셋이다.

  • 밴드 구조 — 이온 퍼텐셜이 E(k)E(\mathbf k) 를 자유전자 포물선에서 휘어 EFE_F 에서의 기울기, 곧 D(EF)D(E_F) 를 바꾼다. 7.9의 유효질량 m∗m^{*} 가 이것이다.
  • 전자–포논 상호작용 — 전자가 격자를 끌고 다니면서 무거워진다. m∗→m∗(1+λ)m^{*}\to m^{*}(1+\lambda) 로, λ\lambda 는 7.14 초전도의 결합 상수와 같은 양이다. Pb 같은 강결합 금속에서 1.5 에 이른다.
  • 전자–전자 상호작용 — 7.10.8의 페르미 액체 보정.

이 셋이 겹쳐 CeAl₃·UPt₃ 같은 무거운 페르미온 화합물에서는 mth/m∼1000m_{\mathrm{th}}/m\sim1000 이 된다 — 자유전자 모형이 가장 극적으로 실패하는 자리이자, γ\gamma 하나로 그 실패를 정량화할 수 있다는 점에서 가장 쓸모 있는 자리다.

반고전 운동 방정식 — 페르미 구가 통째로 움직인다

섹션 제목: “반고전 운동 방정식 — 페르미 구가 통째로 움직인다”

좀머펠트 전자의 운동량은 ℏk\hbar\mathbf k 다. 드루데의 운동 방정식을 p=ℏk\mathbf p=\hbar\mathbf k 로 다시 쓰면, 외력 아래 전자 하나의 k\mathbf k 는

ℏdkdt=F=−e(E+v×B)\hbar\frac{d\mathbf k}{dt}=\mathbf F=-e\bigl(\boldsymbol{\mathcal E}+\mathbf v\times\mathbf B\bigr)

를 따른다 (로렌츠 힘, 2.4.1로렌츠 힘 법칙과 전류; 이 식이 결정 속에서도 성립한다는 것이 7.9.2의 주제다). 자기장이 없고 충돌이 없으면 모든 전자의 k\mathbf k 가 같은 속도로 옮겨간다: δk=−eEt/ℏ\delta\mathbf k=-e\boldsymbol{\mathcal E}t/\hbar. 페르미 구 전체가 kk 공간에서 통째로 밀린다. 충돌은 이 이동을 되돌린다 — 평균 τ\tau 마다 초기화되므로 정상 상태에서 구는

δk=−eτℏE\delta\mathbf k=-\frac{e\tau}{\hbar}\boldsymbol{\mathcal E}

만큼 옮겨진 자리에 멈춘다. 전자 전부가 속도 δv=ℏ δk/m\delta\mathbf v=\hbar\,\delta\mathbf k/m 를 더 가지므로 전류 밀도는

j=−ne δv=−ne ℏm(−eτℏE)=ne2τmE\mathbf j=-ne\,\delta\mathbf v=-ne\,\frac{\hbar}{m}\Bigl(-\frac{e\tau}{\hbar}\boldsymbol{\mathcal E}\Bigr) =\frac{ne^2\tau}{m}\boldsymbol{\mathcal E}

— 드루데와 완전히 같은 σ=ne2τ/m\sigma=ne^2\tau/m 이다. 통계를 바꿨는데 전도도가 그대로인 이유는, 옴의 법칙이 전자 하나하나의 속도 분포가 아니라 전체 운동량의 평균만 묻기 때문이다. 전자 NN 개가 각각 ℏ δk\hbar\,\delta\mathbf k 씩 더 가지면 총 운동량은 통계와 무관하게 Nℏ δkN\hbar\,\delta\mathbf k 다.

전기장은 페르미 구를 통째로 밀고, 산란은 앞쪽 초승달을 뒤쪽으로 돌려보낸다 kx​ ky​ δk = −e𝓔τ/ħ 산란 k → k′ 새로 채워진 상태 비워진 상태 평형 (점선) 정상 상태 𝓔 (+x) 구 안쪽 전자도 전부 δk 만큼 옮겨졌지만, 채워짐이 달라진 곳은 두 초승달뿐이다. 알짜 속도 ≠ 0 인 전자 수 ≈ D(EF​) × (평균 에너지 변화) ∝ n · δk / kF​ 전류 j = −n e ⟨v⟩ = n e² τ 𝓔 / m 산란은 페르미 면 위에서만 일어난다 — 안쪽은 갈 빈자리가 없다(파울리). 평균 자유 행정 l = vF​ τ 열속도가 아니라 페르미 속도가 들어간다.
그림 7.6-5. 전기장 아래의 페르미 구. 점선이 평형 페르미 구(단면), 실선이 +x 방향 전기장 아래 정상 상태의 구다. 모든 전자의 k 가 같은 δk = −e𝓔τ/ħ 만큼 옮겨졌지만, 채워짐이 실제로 달라진 곳은 두 초승달뿐이다 — 왼쪽(−k_x 쪽)에 새로 채워진 상태, 오른쪽에 비워진 상태. 정상 상태를 유지하는 것은 앞쪽 표면의 전자가 뒤쪽 빈 상태로 산란되는 과정이고, 파울리 배타 때문에 구 안쪽 전자는 갈 빈자리가 없어 이 과정에 참여할 수 없다. 그래서 산란은 페르미 면 위의 전자, 곧 속력 v_F 인 전자만 겪고, 평균 자유 행정은 v_F τ 다. δk 는 눈에 보이도록 과장했다 — 배선 수준의 전류 밀도에서 실제 비율 δk/k_F 는 10⁻¹⁰ 안팎이다.직접 작도

무엇이 달라졌나 — 누가 산란되는가

섹션 제목: “무엇이 달라졌나 — 누가 산란되는가”

식은 같지만 그림은 다르다. 드루데에서는 모든 전자가 열속도로 날다가 이온에 부딪혔다. 좀머펠트에서는 다음 넷이 바뀐다.

  1. 산란은 페르미 면 위에서만 일어난다. 산란은 전자를 k\mathbf k 에서 k′\mathbf k' 로 옮기는 일인데, 구 안쪽의 k′\mathbf k' 는 전부 차 있다. 앞쪽 초승달의 전자가 뒤쪽 초승달의 빈자리로 가는 것만 허용된다. 온도가 τ\tau 를 정할 때 관여하는 전자는 이들뿐이다.
  2. 속력은 vFv_F 다. 산란되는 전자가 페르미 면 위에 있으니 평균 자유 행정은 l=vFτl=v_F\tau 이고, vFv_F 는 열속도의 열 배다. 구리 273 K: vF=1.57×106v_F=1.57\times10^{6} m/s, τ=2.7×10−14\tau=2.7\times10^{-14} s, l=42l=42 nm — 원자 간격의 150배다. 77 K 에서는 ρ\rho 가 1/8 이라 l≈300l\approx300 nm, 고순도 단결정 4 K 에서는 mm 급이다.
  3. δk\delta k 는 아주 작다. 배선의 전형적인 전류 밀도 j=1 A/mm2j=1\ \mathrm{A/mm^2} 에서 표류 속도는 vd=j/ne≈0.07v_d=j/ne\approx0.07 mm/s, 곧 δk=mvd/ℏ≈0.6 m−1\delta k=mv_d/\hbar\approx0.6\ \mathrm{m^{-1}} — kF=1.4×1010 m−1k_F=1.4\times10^{10}\ \mathrm{m^{-1}} 의 10−1010^{-10} 이다. 전기장으로 말하면 E=j/σ≈0.2\mathcal E=j/\sigma\approx0.2 mV/cm 에 지나지 않는다. 옴의 법칙이 선형인 이유가 이것이다: 구가 아주 조금만 옮겨져서 응답이 1차항으로 끝난다. 반도체에서는 nn 이 작아 같은 전류에 훨씬 큰 vdv_d 가 필요하고, 전기장이 10410^{4} V/cm 를 넘으면 이 조건이 깨져 속도가 포화한다(11.7).
  4. 이온 격자는 산란체가 아니다. ll 이 원자 간격의 백 배가 넘는다는 것은 전자가 이온 수백 개를 지나면서 부딪히지 않는다는 뜻이다. 완벽한 주기 퍼텐셜은 블로흐 파를 산란시키지 않고(7.7), 산란하는 것은 주기성의 깨짐 — 포논과 결함이다(6절).

ρ=m/ne2τ\rho=m/ne^2\tau 에 τ=l/vF\tau=l/v_F 와 vF=ℏkF/mv_F=\hbar k_F/m 을 넣으면 질량이 사라지고

ρ l=ℏkFne2\rho\,l=\frac{\hbar k_F}{ne^2}

가 남는다. 우변은 전자 농도만의 함수라 금속마다 정해진 상수다. 구리는 6.6×10−16 Ω m26.6\times10^{-16}\ \Omega\,\mathrm{m^2} 이고, 실온 ρ=1.7 μΩ cm\rho=1.7\ \mu\Omega\,\mathrm{cm} 로 나누면 l=39l=39 nm 가 바로 나온다. 이 곱 ρl\rho l 이 6절에서 좁은 배선 금속을 고르는 성적표가 된다.

금속ρ\rho (273 K, μΩ cm\mu\Omega\,\mathrm{cm})τ\tau (10−1410^{-14} s)l=vFτl=v_F\tau (nm)μn=eτ/m\mu_n=e\tau/m (cm2/V s\mathrm{cm^2/V\,s})
Cu1.562.74247
Ag1.514.05670
Au2.043.04253
Al2.450.801614
Na4.23.23456

마지막 열이 뜻밖일 수 있다. 구리의 이동도는 47 cm2/V s47\ \mathrm{cm^2/V\,s} 로, 실리콘 전자의 14001400 보다 30배 작다. 실리콘은 τ\tau 가 비슷하거나 조금 길고 m∗=0.26 mm^{*}=0.26\,m 이라 eτ/m∗e\tau/m^{*} 가 크다. 그런데도 구리의 전도도가 실리콘(1015 cm−310^{15}\ \mathrm{cm^{-3}} 도핑)보다 10810^{8} 배 큰 것은 오로지 nn 이 102210^{22} 대 101510^{15} 이기 때문이다. 1절의 “손잡이 셋” 가운데 금속과 반도체를 가르는 것은 nn 뿐이다.

교류 전도도와 플라스마 주파수

섹션 제목: “교류 전도도와 플라스마 주파수”

전기장이 E(t)=E0e−iωt\boldsymbol{\mathcal E}(t)=\boldsymbol{\mathcal E}_0e^{-i\omega t} 로 흔들리면 운동 방정식 m dv/dt=−eE−mv/τm\,d\mathbf v/dt=-e\boldsymbol{\mathcal E}-m\mathbf v/\tau 의 정상 해는 v=v0e−iωt\mathbf v=\mathbf v_0e^{-i\omega t} 이고, 대입하면 (−iω+1/τ)mv0=−eE0(-i\omega+1/\tau)m\mathbf v_0=-e\boldsymbol{\mathcal E}_0 이다. j=−nev\mathbf j=-ne\mathbf v 로

σ(ω)=σ01−iωτ,σ0=ne2τm\sigma(\omega)=\frac{\sigma_0}{1-i\omega\tau}, \qquad \sigma_0=\frac{ne^2\tau}{m}

이다. ωτ≪1\omega\tau\ll1 이면 직류와 같고, ωτ≫1\omega\tau\gg1 이면 σ≈i ne2/mω\sigma\approx i\,ne^2/m\omega 로 순허수 — 전자가 마찰 없이 관성만으로 흔들리는 상태다. 실수부(흡수)와 허수부(관성)가 ωτ=1\omega\tau=1 에서 교차한다. 구리는 1/τ≈4×10131/\tau\approx4\times10^{13} rad/s, 곧 적외선(파장 ∼50 μ\sim50\ \mum)이 경계다. 그 아래(마이크로파·THz 까지)에서 금속은 직류와 같은 저항체이고, 그 위(가시광)에서는 관성체다.

관성체의 유전 응답을 보자. 전도 전류를 분극으로 읽어 ε(ω)=1+iσ/ε0ω\varepsilon(\omega)=1+i\sigma/\varepsilon_0\omega 로 쓰고 ωτ≫1\omega\tau\gg1 을 넣으면

ε(ω)=1−ωp2ω2,  ωp2=ne2ε0m  \varepsilon(\omega)=1-\frac{\omega_p^2}{\omega^2}, \qquad \boxed{\;\omega_p^2=\frac{ne^2}{\varepsilon_0m}\;}

— ωp\omega_p 가 플라스마 주파수다. ω<ωp\omega<\omega_p 에서 ε<0\varepsilon<0 이라 굴절률 ε\sqrt\varepsilon 이 허수이고, 파는 매질에 들어가지 못하고 반사된다 — 금속의 광택이다. ω>ωp\omega>\omega_p 에서는 ε>0\varepsilon>0 이라 투과한다. ℏωp\hbar\omega_p 는 Na 5.9 eV(파장 210 nm), Al 15.8 eV, Cu 10.8 eV 로, 단순 금속은 자외선에서 투명해진다. 구리·금이 색을 띠는 것은 dd 밴드에서 ss 밴드로의 밴드간 흡수가 가시광에 걸려서다(7.11).

ωτ = 1 이 저항과 관성의 경계, ω = ωp​ 가 반사와 투과의 경계 σ(ω) / σ0​ ωτ 1 2 3 4 1 ½ Re σ — 흡수(저항) ∝ 1/ω² 로 사라진다 Im σ — 관성(유도성) ∝ 1/ω 로 남는다 ε(ω) = 1 − ωp​²/ω² ω -3 -2 -1 1 ωp​ ε < 0 반사 (금속 광택) ε > 0 투과 (자외선 투명) → 1
그림 7.6-6. 드루데 교류 전도도와 플라스마 주파수. 왼쪽: σ(ω) = σ₀/(1 − iωτ) 의 실수부와 허수부. ωτ ≪ 1 이면 직류와 같고, ωτ = 1 에서 둘이 ½ 에서 만나며, ωτ ≫ 1 이면 실수부(흡수)는 1/ω² 로 사라지고 허수부(관성)만 1/ω 로 남는다. 오른쪽: 충돌을 무시한 유전함수 ε(ω) = 1 − ω_p²/ω². ω < ω_p 에서 음수 — 파가 들어가지 못하고 반사된다(금속 광택). ω > ω_p 에서 양수 — 투명하다. 알칼리 금속은 자외선(Na: ~210 nm)에서 투명해진다.직접 작도

마티센 규칙 — 산란율은 더해진다

섹션 제목: “마티센 규칙 — 산란율은 더해진다”

전자를 산란시키는 상대가 둘 이상이면 — 포논과 불순물 — 각각의 산란이 독립이라는 가정 아래 산란율이 더해진다:

1τ=1τph+1τi⟹ρ=mne2(1τph+1τi)=ρph(T)+ρ0\frac1\tau=\frac{1}{\tau_{\mathrm{ph}}}+\frac{1}{\tau_{\mathrm i}} \qquad\Longrightarrow\qquad \rho=\frac{m}{ne^2}\Bigl(\frac{1}{\tau_{\mathrm{ph}}}+\frac{1}{\tau_{\mathrm i}}\Bigr) =\rho_{\mathrm{ph}}(T)+\rho_0

이것이 마티센 규칙(Matthiessen’s rule)이다. 확률의 덧셈이라 τ\tau 가 아니라 1/τ1/\tau 가 더해지고, 그래서 비저항이 더해진다 — 전도도가 아니라. 포논 항 ρph\rho_{\mathrm{ph}} 는 온도에 의존하고 T→0T\to0 에서 사라지지만, 불순물·결함 항 ρ0\rho_0 는 온도와 무관하다. 그래서 순수하지 않은 금속의 비저항은 절대영도에서 0 이 아니라 잔류 비저항 ρ0\rho_0 에 멈추고, ρ(300 K)/ρ(4 K)\rho(300\ \mathrm K)/\rho(4\ \mathrm K) — 잔류 저항비(RRR) — 가 시료 순도의 표준 척도가 된다. 상용 구리 배선은 수십, 고순도 단결정은 수천이다.

잔류 비저항은 곡선을 위로 올릴 뿐 모양을 바꾸지 않는다 — 마티센 규칙 ρ / ρ_ph(Θ) T / Θ 0.25 0.5 0.75 1 1.25 0.5 1.0 1.5 고온 점근 ∝ T (점선) 합금 · 결함 많음 보통 순도 이상적 순금속 ρ0​ (잔류) ∝ T⁵ 기울기가 같다 — 1/τ = 1/τph​ + 1/τi​ 의 합 RRR = ρ(300 K)/ρ(4 K) 는 순도의 척도 — 상용 Cu 배선 ~ 수십, 고순도 단결정 ~ 수천.
그림 7.6-7. 금속 비저항의 온도 의존성 — 마티센 규칙 ρ = ρ₀ + ρ_ph(T). 포논 항은 블로흐–그뤼나이젠 식 ρ_ph ∝ (T/Θ)⁵ J₅(Θ/T) 로 계산했다. T ≳ Θ/3 에서 T 에 비례하고, T ≪ Θ 에서 T⁵ 로 사라진다. 잔류 비저항 ρ₀ 는 곡선을 통째로 위로 올릴 뿐 모양을 바꾸지 않는다 — 세 곡선의 기울기가 같다. 실온 구리(Θ_R ≈ 320 K, T/Θ ≈ 0.9)는 직선 구간 안에 있어 온도 계수가 0.4 %/K 다.직접 작도

포논 산란 — 왜 고온에서 T, 저온에서 T⁵ 인가

섹션 제목: “포논 산란 — 왜 고온에서 T, 저온에서 T⁵ 인가”

포논은 격자 주기성의 흔들림이고, 전자는 그 흔들림에 산란된다. 산란율은 포논 수에 비례한다. 고온(T≫ΘT\gg\Theta)에서 각 모드의 포논 수는 등분배로 nph≈kBT/ℏωn_{\mathrm{ph}}\approx k_BT/\hbar\omega (7.5.1포논 열용량과 플랑크 분포) 라 TT 에 비례하고, 따라서

ρph∝T(T≳Θ/3)\rho_{\mathrm{ph}}\propto T \qquad(T\gtrsim\Theta/3)

이다. 실온 구리의 온도 계수 dρ/dT=0.0068 μΩ cm/Kd\rho/dT=0.0068\ \mu\Omega\,\mathrm{cm/K}, 곧 0.4 %/K 가 이것이다.

저온(T≪ΘT\ll\Theta)에서는 두 가지가 겹쳐 지수가 5 로 뛴다. ① 열적으로 들뜬 포논은 ℏω≲kBT\hbar\omega\lesssim k_BT 인 것뿐이고 그 수는 T3T^3 에 비례한다(디바이 T3T^3 법칙과 같은 셈). ② 그런 포논의 파수는 q≲kBT/ℏvs≪kFq\lesssim k_BT/\hbar v_s\ll k_F 라, 전자를 각도 θ≈q/kF∝T\theta\approx q/k_F\propto T 만큼만 꺾는다. 5절의 주의에서 말한 대로 운동량 완화에는 (1−cos⁡θ)≈θ2/2∝T2(1-\cos\theta)\approx\theta^2/2\propto T^2 의 무게가 붙는다. 곱하면

ρph∝T3×T2=T5(T≪Θ)\rho_{\mathrm{ph}}\propto T^3\times T^2=T^5 \qquad(T\ll\Theta)

— 블로흐 T5T^5 법칙이다. 둘을 잇는 보간식이 블로흐–그뤼나이젠 식

ρph(T)=A(TΘR)5∫0ΘR/Tx5 dx(ex−1)(1−e−x)\rho_{\mathrm{ph}}(T)=A\Bigl(\frac{T}{\Theta_R}\Bigr)^5\int_0^{\Theta_R/T}\frac{x^5\,dx}{(e^x-1)(1-e^{-x})}

이고, 그림의 곡선이 이것이다. ΘR\Theta_R 은 비저항에서 맞춘 디바이 온도로 열용량의 ΘD\Theta_D 와 비슷하다(Cu 320 K, Al 395 K). 실제 금속에서는 움클랍 과정(7.5.8움클랍 과정)이 저온의 지수를 조금 바꾸고, 페르미 면이 브릴루앙 영역에 닿는 금속에서는 더 복잡해지지만, “고온 TT, 저온 T5T^5, 사이는 블로흐–그뤼나이젠”이 표준 그림이다.

불순물 원자 하나는 산란 단면적을 가진 점 산란체다. 농도 xx (원자 분율)가 낮으면 ρi∝x\rho_{\mathrm i}\propto x 이고, 이원 합금 전체로는 노르트하임 규칙 ρi∝x(1−x)\rho_{\mathrm i}\propto x(1-x) — 두 성분이 반반일 때 최대다. 산란 세기는 불순물이 모체와 얼마나 다른가에 달렸고, 같은 주기의 원소를 구리에 1 at % 씩 넣으면 Δρ\Delta\rho 가 Zn 0.3, Ga 1.4, Ge 3.7, As 6.7 μΩ cm\mu\Omega\,\mathrm{cm} 안팎으로 원자가 차이의 제곱에 비례한다(린데 규칙). 여분의 전하가 만드는 스크린된 쿨롱 퍼텐셜(7.10.3)의 산란 단면적이 (ΔZ)2(\Delta Z)^2 이기 때문이다. 공공·전위·입계도 각각 ρ0\rho_0 에 더해진다 — 급랭한 금속의 공공, 냉간 가공한 금속의 전위, 미세 결정립의 입계가 어닐링으로 줄어드는 것을 ρ0\rho_0 로 추적한다.

한 가지 예외를 적어 둔다. Fe·Mn 같은 자성 불순물이 든 Cu·Au 는 저온에서 비저항이 내려가다 다시 오른다 — ρ\rho 에 −ln⁡T-\ln T 항이 붙어 극소가 생긴다(곤도 효과). 불순물의 스핀이 전도 전자 스핀과 교환 상호작용하면서 산란 진폭이 온도에 의존하는 결과로, 마티센 규칙의 “온도 무관 ρ0\rho_0” 가 깨지는 드문 경우다.

크기 효과 — 배선 폭이 평균 자유 행정과 같아질 때

섹션 제목: “크기 효과 — 배선 폭이 평균 자유 행정과 같아질 때”

5절의 l≈40l\approx40 nm 는 이제 벌크 상수가 아니라 설계 길이다. 박막의 두께 dd 나 배선의 폭 ww 가 ll 과 비슷해지면 표면이 새로운 산란체가 된다. 표면에서 전자가 거울처럼 반사되면 (specular, 확률 pp) 운동량이 보존되어 비저항이 오르지 않지만, 난반사되면(diffuse, 1−p1-p) 표면이 충돌 한 번과 같다. 푸흐스–존트하이머(Fuchs–Sondheimer) 이론의 근사식은

ρρbulk≈1+38(1−p)ld(d≳l)\frac{\rho}{\rho_{\mathrm{bulk}}}\approx1+\frac38(1-p)\frac{l}{d} \qquad(d\gtrsim l)

이고, 다결정 배선에서는 결정립계가 전자를 반사(확률 RR)해 마야다스–샤츠케스 (Mayadas–Shatzkes) 항이 따로 더해진다. 구리 배선에서 p≈0p\approx0, R≈0.3∼0.5R\approx0.3\sim0.5 로 잡히며, 폭 20 nm 배선의 비저항은 벌크의 2~3배다. 마티센 규칙의 항이 하나 더 늘어난 것이다: ρ=ρph+ρ0+ρsurf+ρgb\rho=\rho_{\mathrm{ph}}+\rho_0+\rho_{\mathrm{surf}}+\rho_{\mathrm{gb}}.

이 자리에서 5절의 ρl=ℏkF/ne2\rho l=\hbar k_F/ne^2 가 성적표가 된다. 배선 폭 w≪lw\ll l 인 극한에서 비저항은 ρbulk l/w\rho_{\mathrm{bulk}}\,l/w 에 비례하므로, 좁은 배선의 저항은 ρ\rho 가 아니라 ρl\rho l 곱으로 정해진다. Co·Ru·Mo 는 벌크 비저항이 Cu 의 3~4배지만 ll 이 10 nm 안팎으로 짧아 ρl\rho l 이 Cu 와 비슷하거나 작다 — 폭 10 nm 대에서 Cu 를 대체하는 후보가 되는 이유이고, 18.22와 10.4.8에서 이어진다.

자기장 속의 운동 방정식과 사이클로트론 주파수

섹션 제목: “자기장 속의 운동 방정식과 사이클로트론 주파수”

로렌츠 힘을 넣은 드루데 방정식은

mdvdt=−e(E+v×B)−mvτm\frac{d\mathbf v}{dt}=-e\bigl(\boldsymbol{\mathcal E}+\mathbf v\times\mathbf B\bigr)-\frac{m\mathbf v}{\tau}

이다. 전기장과 충돌을 잠시 끄면 전자는 자기장에 수직인 면에서 원운동을 하고, 그 각진동수가 사이클로트론 주파수

ωc=eBm(fc=ωc/2π=28 GHz/T)\omega_c=\frac{eB}{m} \qquad(f_c=\omega_c/2\pi=28\ \mathrm{GHz/T})

다. 충돌이 있으면 전자는 한 바퀴를 채 돌기 전에 산란된다 — 기준은 ωcτ\omega_c\tau 다. ωcτ≪1\omega_c\tau\ll1 이면 자기장은 궤도를 조금 휘는 섭동이고, ωcτ≫1\omega_c\tau\gg1 이면 전자가 여러 바퀴를 돈다. 실온 구리에서 1 T 는 ωcτ≈5×10−3\omega_c\tau\approx5\times10^{-3} 로 약한 장이고, 4 K 고순도 구리(RRR 1000)에서는 τ\tau 가 천 배라 같은 1 T 가 ωcτ≈5\omega_c\tau\approx5 의 강한 장이다. 강한 장에서 전자가 그리는 kk 공간의 궤도가 페르미 면의 단면이고, 그것이 7.8의 드하스–판알펜 효과와 궤도 양자화로 간다.

전류가 xx 방향, 자기장이 zz 방향인 막대를 본다. 전자는 −x-x 로 움직이고 로렌츠 힘 −ev×B-e\mathbf v\times\mathbf B 를 −y-y 방향으로 받아 아래쪽 가장자리로 밀린다. 그러나 막대는 폐회로가 아니므로 yy 방향 전류는 흐를 수 없다 — 전자가 아래 가장자리에 쌓이고 위는 양전하가 남아, 로렌츠 힘을 정확히 상쇄하는 횡전기장 Ey\mathcal E_y 가 생긴다. 이것이 홀 전기장이다.

정상 상태(dv/dt=0d\mathbf v/dt=0)의 두 성분 방정식을 적으면

0=−eEx−evyB−mvxτ,0=−eEy+evxB−mvyτ\begin{aligned} 0&=-e\mathcal E_x-ev_yB-\frac{m v_x}{\tau},\\ 0&=-e\mathcal E_y+ev_xB-\frac{m v_y}{\tau} \end{aligned}

이고, jy=0j_y=0 이므로 vy=0v_y=0 이다. 첫 식에서 vx=−eτEx/mv_x=-e\tau\mathcal E_x/m (드루데 그대로), 둘째 식에서

Ey=vxB=−eτmB Ex=−ωcτ Ex\mathcal E_y=v_xB=-\frac{e\tau}{m}B\,\mathcal E_x=-\omega_c\tau\,\mathcal E_x

이다. jx=−nevxj_x=-nev_x 로 vxv_x 를 지우면 Ey=−jxB/ne\mathcal E_y=-j_xB/ne 가 되고, 홀 계수를 정의하면

  RH≡EyjxB=−1ne  \boxed{\;R_H\equiv\frac{\mathcal E_y}{j_xB}=-\frac{1}{ne}\;}

τ\tau 가 사라졌다. 홀 계수는 산란과 무관하게 캐리어의 부호와 농도만 담는다 — 부호는 전하의 부호, 크기는 농도의 역수. 이것이 홀 효과가 반도체 계측의 중심에 있는 이유다.

전류는 +x, 자기장은 지면 밖 — 전자는 아래로 밀리고 홀 전기장이 그것을 막는다 B ⊙ (지면 밖, +z) jx​ jx​ −e v F = −e v × B + + + + + + + + + + − − − − − − − − − − 𝓔H​ (−y) VH​ = 𝓔H​ w − (아래가 낮다) x y z (지면 밖) 정상 상태: −e 𝓔H​ = e v B → 𝓔H​ = −vB = −jx​B/(ne) RH​ = 𝓔H​/(jx​B) = −1/(ne) — 부호는 캐리어의 부호, 크기는 1/농도. τ 가 사라진다.
그림 7.6-8. 홀 효과의 기하. 전류 j_x 는 +x, 자기장 B 는 지면 밖(+z). 전자는 −x 로 움직이며 로렌츠 힘 F = −e v × B 를 −y 방향으로 받아 아래 가장자리에 쌓이고, 위 가장자리에는 양전하가 남는다. 그 결과 아래를 향하는 홀 전기장 𝓔_H 가 생겨 로렌츠 힘과 맞서고, 위아래 사이의 전압이 홀 전압 V_H = 𝓔_H w 다. 정상 상태의 힘의 균형 e𝓔_H = e v B 에서 𝓔_H = −j_x B/(ne) — 완화 시간이 사라지고 캐리어의 부호와 농도만 남는다. 캐리어가 정공이면 같은 전류에 대해 위아래가 뒤집힌다.직접 작도

같은 계산에서 두 가지가 더 나온다. 홀 저항률 ρxy=Ey/jx=B/ne\rho_{xy}=\mathcal E_y/j_x=B/ne 는 BB 에 정비례하고, 세로 저항률 ρxx=Ex/jx=m/ne2τ\rho_{xx}=\mathcal E_x/j_x=m/ne^2\tau 는 BB 와 무관하다 — 드루데 모형의 자기저항은 0 이다. 실제 금속은 대부분 B2B^2 에 비례하는 양의 자기저항을 보이는데, 그것은 페르미 면이 구가 아니거나 캐리어가 두 종류일 때 생기고 7.8·7.12.7에서 다룬다.

RHR_H 를 자유전자 값으로 나눈 비 −1/(RHne)-1/(R_Hne) 를 보면, 원자가 전자 수로 계산한 nn 이 얼마나 맞는지가 보인다 (Ashcroft·Mermin Table 1.4, 강한 장):

금속LiNaKRbCsCuAgAuBeMgAl
−1/(RHne)-1/(R_Hne)0.81.21.11.00.91.51.31.5−0.2-0.2−0.4-0.4−0.3-0.3

알칼리 금속은 1 에 가깝다 — 자유전자 모형의 가장 깨끗한 성공이다. 귀금속은 1.3~1.5 로 어긋나고, Be·Mg·Al 은 부호가 반대다: 전류를 나르는 것이 양전하처럼 보인다. 알루미늄의 −0.3-0.3 은 원자당 정공 하나(nh=n/3n_h=n/3)로 읽으면 정확히 맞는다. 자유전자 모형 안에는 양전하 캐리어가 없으므로, 이것은 통계로도 산란으로도 고칠 수 없는 균열이다. 답은 7.9.3의 정공 — 거의 찬 밴드의 빈자리가 양전하·양질량의 입자처럼 행동한다는 것 — 이고, 그 이전에 7.7의 밴드 갭이 있어야 “거의 찬 밴드”가 정의된다. 1절의 성적표에서 홀 계수를 △로 둔 이유다.

반도체에서 — 캐리어를 세는 자

섹션 제목: “반도체에서 — 캐리어를 세는 자”

반도체에서 홀 효과는 자유전자 모형이 그대로 쓰이는 자리다. 한 종류의 캐리어만 있으면

RH=−1ne (n형),RH=+1pe (p형)R_H=-\frac{1}{ne}\ (\text{n형}), \qquad R_H=+\frac{1}{pe}\ (\text{p형})

이고, 홀 전압의 부호가 캐리어의 종류, 크기가 농도를 준다. 전도도와 합쳐 홀 이동도 μH=∣RH∣σ\mu_H=|R_H|\sigma 를 얻는데, 산란 기구에 따라 실제 이동도와 rH=1.1∼1.2r_H=1.1\sim1.2 의 홀 인자만큼 다르다(포논 산란 3π/83\pi/8, 이온화 불순물 산란 315π/512315\pi/512 — 7.12의 완화시간 평균에서 나온다). 두 캐리어가 함께 있으면

RH=pμp2−nμn2e (pμp+nμn)2R_H=\frac{p\mu_p^2-n\mu_n^2}{e\,(p\mu_p+n\mu_n)^2}

으로 이동도의 제곱이 무게가 되어, 진성에 가까운 시료는 수가 적어도 빠른 전자가 부호를 정한다. 11.2.9의 홀 측정과 반 데르 파우 방법, 19부의 계측이 전부 이 식 위에 있다.

수를 넣어 보면 금속과 반도체의 자리가 보인다. 두께 tt 인 판에 전류 II 가 흐를 때 홀 전압은 VH=IB/netV_H=IB/net 이다. 구리 판(t=100 μt=100\ \mum, I=1I=1 A, B=1B=1 T)은 VH≈0.7 μV_H\approx0.7\ \muV — 잡음 수준이다. n=1015 cm−3n=10^{15}\ \mathrm{cm^{-3}} 실리콘(t=500 μt=500\ \mum, I=1I=1 mA, B=1B=1 T)은 VH≈12V_H\approx12 mV — 만 배가 넘는다. 홀 센서(17부)를 반도체로 만드는 이유가 RH∝1/nR_H\propto1/n 이고, 같은 이유로 금속의 RHR_H 측정은 정밀 실험이다. 강한 장·저온·2차원 전자 기체에서 ρxy\rho_{xy} 가 h/νe2h/\nu e^2 의 계단으로 양자화되는 양자 홀 효과는 10.4.7에서 다룬다.

기체 운동론의 열전도도 — 이번에는 전자로

섹션 제목: “기체 운동론의 열전도도 — 이번에는 전자로”

기체의 열전도도는 3.7·7.5에서 쓴 공식(7.5.7열전도와 포논 기체의 열저항)

κ=13 CV v l\kappa=\frac13\,C_V\,v\,l

이다 — 단위 부피당 열용량 × 속력 × 평균 자유 행정. 포논 기체에 썼던 식을 전자 기체에 그대로 쓴다. 그런데 이제 세 인자의 값을 안다. 부피당 열용량은 4절의 CV=π22nkB(T/TF)C_V=\frac{\pi^2}{2}nk_B(T/T_F), 속력은 5절의 vFv_F, 자유 행정은 l=vFτl=v_F\tau 다. 넣으면

κel=13⋅π22 nkBkBTEF⋅vF2τ=π23 nkB2Tτm\kappa_{\mathrm{el}}=\frac13\cdot\frac{\pi^2}{2}\,nk_B\frac{k_BT}{E_F}\cdot v_F^2\tau =\frac{\pi^2}{3}\,\frac{nk_B^2T\tau}{m}

이다 (EF=12mvF2E_F=\tfrac12mv_F^2 를 썼다). 고전 기체와 다른 점이 명확하다 — 열을 나르는 것은 T/TFT/T_F 의 소수 전자뿐이지만, 그 전자들이 vFv_F 로 날아 ll 이 길다. 두 효과가 곱해진다.

비데만–프란츠 법칙과 로렌츠 수

섹션 제목: “비데만–프란츠 법칙과 로렌츠 수”

이 κ\kappa 를 같은 τ\tau 의 σ=ne2τ/m\sigma=ne^2\tau/m 으로 나누면 nn, τ\tau, mm 이 전부 사라진다:

  κσT=π23(kBe)2≡L0=2.44×10−8 W Ω K−2  \boxed{\;\frac{\kappa}{\sigma T}=\frac{\pi^2}{3}\Bigl(\frac{k_B}{e}\Bigr)^2\equiv L_0 =2.44\times10^{-8}\ \mathrm{W\,\Omega\,K^{-2}}\;}

— 비데만–프란츠 법칙(1853)이고, L0L_0 이 로렌츠 수다. 금속의 종류, 순도, 온도에 관계없이 열전도도와 전기전도도의 비가 온도에만 비례한다. 우변에 기본 상수 둘뿐이라 L0L_0 은 보편 상수다. 실측 L=κ/σTL=\kappa/\sigma T 를 적는다 (10−8 W Ω K−210^{-8}\ \mathrm{W\,\Omega\,K^{-2}}, Kittel Table 6.5):

금속AgAuCuPbPtWZn
0 °C2.312.352.232.472.513.042.31
100 °C2.372.402.332.562.603.202.33

Cu 와 Ag 는 5 % 안에서 맞고, 나머지도 30 % 안이다. 구리 실온의 κ=L0Tσ=2.44×10−8×273×6.4×107≈430\kappa=L_0T\sigma =2.44\times10^{-8}\times273\times6.4\times10^{7}\approx430 W/m·K 는 실측 400 과 같은 값이다. 금속의 열전도는 거의 전부 전자가 나른다 — 포논의 몫은 1~2 % 로, 절연체(Si 150, 다이아몬드 2000 W/m·K)에서 포논만으로 내는 열전도와 자릿수는 같지만 금속에서는 전자가 그것을 압도한다.

드루데 모형으로 같은 계산을 하면 CV=32nkBC_V=\tfrac32nk_B, 12mv2=32kBT\tfrac12mv^2=\tfrac32k_BT 로

κσT=32(kBe)2=1.11×10−8 W Ω K−2\frac{\kappa}{\sigma T}=\frac32\Bigl(\frac{k_B}{e}\Bigr)^2=1.11\times10^{-8}\ \mathrm{W\,\Omega\,K^{-2}}

— 옳은 값의 절반이다. 그런데 드루데의 1900년 논문은 계산 실수로 이것의 두 배를 얻어 2.2×10−82.2\times10^{-8} 을 보고했고, 실측과 “놀랍게” 맞았다. 잘못된 열용량(100배 큼)과 잘못된 속도 제곱(100배 작음)이 서로 지웠고 그 위에 산술 오류가 겹쳐 맞은 것이다. 이 우연이 드루데 모형의 신뢰를 한 세대 동안 지탱했다. 좀머펠트가 두 100배를 각각 바로잡자 비는 그대로 상수로 남았고 값이 π2/3\pi^2/3 으로 정확해졌다 — 상쇄가 우연이 아니라 같은 전자가 두 흐름을 나른다는 구조였음이 드러난 셈이다.

어디서 깨지나 — 탄성 산란과 비탄성 산란

섹션 제목: “어디서 깨지나 — 탄성 산란과 비탄성 산란”

법칙의 조건은 “전하와 열이 같은 τ\tau 를 겪는다”는 것이었다. 불순물 산란은 탄성이다 — 전자의 에너지는 그대로이고 방향만 바뀐다. 이 산란은 전하의 흐름과 열의 흐름을 같은 비율로 깎으므로 L=L0L=L_0 이 정확히 성립한다. 고온의 포논 산란도 ℏωph≪kBT\hbar\omega_{\mathrm{ph}}\ll k_BT 라 거의 탄성이고, 역시 L≈L0L\approx L_0 이다.

중간 온도(T∼Θ/10T\sim\Theta/10)의 포논 산란은 다르다. 포논 에너지 ℏω∼kBT\hbar\omega\sim k_BT 가 전자 에너지의 창 폭과 같아, 한 번의 산란이 전자의 에너지를 창 밖으로 옮겨 열류를 지운다 — 방향은 작은 각도만 바뀌어 전류는 거의 그대로인데. 열류가 전류보다 더 잘 깎이니 κ\kappa 가 상대적으로 작고 L<L0L<L_0 이다. 순수한 금속에서 L/L0L/L_0 이 Θ/10\Theta/10 근처에서 0.3~0.5 까지 내려가는 이유이고, 불순물이 많으면 탄성 산란이 지배해 L0L_0 으로 돌아간다.

열전도도는 봉우리를 만든다 — 저온은 불순물이, 고온은 포논이 정한다 κ / κ(T ≫ Θ) T / Θ 0.5 1 1 2 3 4 ρ0​ = 0.05 ρ0​ = 0.01 (더 순수) ∝ T (불순물) → 일정 (포논, ρ ∝ T) L / L0​ T / Θ 0.5 1 0.5 1.0 포논만 (순수 극한, 점선) 불순물이 L 을 L0​ 로 되돌린다 탄성 산란 → L = L0​
그림 7.6-9. 금속의 열전도도와 로렌츠 수. 전기 비저항은 ρ = ρ₀ + ρ_ph(T), 열저항은 W = ρ₀/(L₀T) + W_ph(T) 로 두고 κ = 1/W 를 그렸다(Wilson 식). 왼쪽: κ(T) 는 저온에서 T 에 비례해 오르고(불순물 산란, C ∝ T), 봉우리를 만든 뒤 고온에서 일정해진다(포논 산란, ρ ∝ T 가 C ∝ T 를 지운다). 더 순수한 시료일수록 봉우리가 높고 낮은 온도에 있다 — 고순도 구리는 20 K 근처에서 실온의 열 배가 넘는 κ 를 낸다. 오른쪽: L/L₀ 는 저온(불순물)과 고온(고온 포논)에서 1 이고, 중간 온도에서 비탄성 포논 산란 때문에 1 아래로 내려간다. 순수 극한(점선)에서는 0 으로 간다.직접 작도

열전 효과 — 창의 기울기가 만드는 전압

섹션 제목: “열전 효과 — 창의 기울기가 만드는 전압”

같은 창에서 하나가 더 나온다. 온도 기울기가 있으면 뜨거운 쪽에서 오는 전자와 차가운 쪽에서 오는 전자의 에너지 분포가 다르고, D(E)D(E) 나 τ(E)\tau(E) 가 에너지에 따라 변하면 두 흐름이 상쇄되지 않아 알짜 전류 — 열린 회로에서는 전압 — 가 생긴다. 제벡 계수 SS 는 전류가 0 일 때의 E=S ∇T\mathcal E=S\,\nabla T 로 정의되고, 좀머펠트 전개로 (모트 공식)

S=−π23kB2Te dln⁡σ(E)dE∣EFS=-\frac{\pi^2}{3}\frac{k_B^2T}{e}\,\frac{d\ln\sigma(E)}{dE}\Bigl|_{E_F}

이다 — σ(E)\sigma(E) 는 “에너지 EE 인 전자만 있을 때의 전도도”이고, 그 로그 미분이 답이다. 자유전자에서 σ(E)∝E3/2\sigma(E)\propto E^{3/2} (전자 수 ∝E3/2\propto E^{3/2}, τ\tau 일정) 라

S=−π22kBeTTF≈−1.6 μV/K(Cu, 300 K)S=-\frac{\pi^2}{2}\frac{k_B}{e}\frac{T}{T_F}\approx-1.6\ \mu\mathrm{V/K} \qquad(\text{Cu, 300 K})

— 마이크로볼트다. 드루데의 고전 계산은 T/TFT/T_F 인자가 없어 −32kB/e≈−130 μ-\frac32k_B/e\approx-130\ \muV/K, 100배 크다. 1절 성적표의 마지막 100배가 여기서 닫힌다. 구리의 실측은 +1.8 μ+1.8\ \muV/K 로 크기는 맞고 부호가 반대다 — 홀 계수와 같은 종류의 실패이고, EFE_F 에서 σ(E)\sigma(E) 의 기울기가 밴드 구조 때문에 음수인 탓이다. 반도체의 SS 가 수백 μ\muV/K 인 것은 EFE_F 가 밴드 가장자리 근처에 있어 dln⁡σ/dE∼1/kBTd\ln\sigma/dE\sim1/k_BT 로 크기 때문이고, 그것이 열전 소자와 열전대의 원리다(7.12.7, 17부).

반도체에서 — 열은 배선을 따라 나간다

섹션 제목: “반도체에서 — 열은 배선을 따라 나간다”

칩의 열은 트랜지스터에서 나지만 실리콘과 절연막(SiO₂ 1.4 W/m·K)을 지나 밖으로 나가는 길은 좁다. 금속 배선과 비아, 범프, 히트 스프레더가 열의 고속도로이고, 그 열전도도는 재지 않아도 전기 비저항에서 L0T/ρL_0T/\rho 로 읽는다. TCAD 의 전열 해석(20부)이 금속 영역의 κ\kappa 를 이렇게 넣고, 배선의 줄 발열 ρj2\rho j^2 이 같은 ρ\rho 로 커지면서 κ\kappa 는 1/ρ1/\rho 로 줄어드는 결합이 자기 발열의 양성 되먹임을 만든다(3.7.13·18부). 좁은 배선에서 6절의 크기 효과로 ρ\rho 가 세 배 오르면 κ\kappa 는 그만큼 내려간다 — 두 문제가 하나의 τ\tau 다.

이 장이 반도체에서 다시 나오는 곳

섹션 제목: “이 장이 반도체에서 다시 나오는 곳”

한 문장으로 줄이면 — 전자 기체의 성질은 EFE_F 근처 kBTk_BT 폭의 창이 정하고, 그 창의 상태 수 D(EF)D(E_F)·속력 vFv_F·기억 시간 τ\tau 가 열용량·전도·열전도의 전부다. 목적지를 적어 둔다.

  • 캐리어 통계 (11.1) — D(E)∝E−EcD(E)\propto\sqrt{E-E_c}, f(E)f(E), NcN_c, 축퇴 반도체의 페르미–디랙 적분(3절·2절). 2차원 계단 상태밀도와 부밴드(13부).
  • 이동도와 산란 (11.2·7.12) — σ=neμn\sigma=ne\mu_n, μn=eτ/m∗\mu_n=e\tau/m^{*}, 마티센 규칙, 포논·이온화 불순물 산란의 온도 의존성(5·6절). 볼츠만 방정식이 τ\tau 를 정의한다.
  • 홀 측정 (11.2.9·19부) — RHR_H, 홀 이동도, 두 캐리어 공식, 반 데르 파우(7절). 양자 홀 효과(10.4.7).
  • 배선 (18.22·21.7·10.4.8) — l=40l=40 nm 와 크기 효과, ρl\rho l 성적표, 온도 계수(6절), 비데만–프란츠로 읽는 열전도(8절).
  • 광학 (7.10·7.11·11.8·16부) — σ(ω)\sigma(\omega), ωp\omega_p, 자유 캐리어 흡수, 투명 전극, 플라스몬(5절).
  • 다음 장들. 이온을 돌려놓는다 — 블로흐 정리와 밴드 갭(7.7), 유효질량과 정공(7.9), 전자 기체의 차폐와 페르미 액체(7.10). 자유전자 모형이 틀린 세 곳(홀 부호·열전 부호· ρ(T)\rho(T))이 거기서 고쳐진다.