Chapter 02

플라즈마의 척도: 데바이·주파수·준중성

이온화된 기체가 "플라즈마"가 되려면 하전 입자가 충분히 많아서 서로를 가리고, 함께 움직이고, 전체로는 중성이어야 한다. 이 장은 그 말을 세 개의 숫자 — 데바이 길이 λ_D, 플라즈마 주파수 f_pe, 플라즈마 파라미터 N_D — 로 바꾼다. PL-300 IC 모듈의 기준 플라즈마(Ar 10 mTorr, 1.5 kW → Tₑ ≈ 2.5 eV, nₑ ≈ 4×10¹¹ cm⁻³)에서 λ_D ≈ 20 µm, f_pe ≈ 6 GHz, N_D ≈ 10⁴이다. 이 세 수가 왜 시스가 수백 µm인지, 왜 13.56 MHz 전기장이 플라즈마 안에 들어가지 못하는지, 왜 2.45 GHz 마이크로파 플라즈마는 밀도 상한이 있는지를 설명한다.

작은 전하 불균형이 만드는 전위를 포아송 식으로 계산하고, 준중성이 성립하는 조건 δ ≪ (λ_D/L)²를 유도할 수 있다.
볼츠만 전자를 넣은 포아송 식을 풀어 데바이 차폐를 보이고, 선형화가 언제 맞고 언제 틀리는지 수치 해로 확인할 수 있다.
전자·이온 플라즈마 주파수를 계산하고, 공정 RF 주파수(0.4–60 MHz, 2.45 GHz)에 전자와 이온이 각각 따라가는지 판단할 수 있다.
플라즈마 파라미터 N_D와 결합 파라미터 Γ로 공정 플라즈마가 약결합 이상 기체임을 보이고, 사하 식으로 공정 플라즈마가 열평형에서 얼마나 먼지 설명할 수 있다.
차단 밀도로 전자기파 투과 여부를 판정하고, nₑ·Tₑ·압력·크기로 플라즈마의 세 조건을 판정할 수 있다.

준중성: 전하는 왜 거의 정확히 상쇄되는가

공정 플라즈마의 몸통(벌크)에서 전자 밀도 $n_e$ 와 이온 밀도 $n_i$ 는 거의 같다. 이 성질을 준중성(quasineutrality)이라 하고 $n_e \approx n_i \approx n$ 으로 쓴다. "거의"가 얼마나 정확한지는 포아송 식이 알려 준다. 폭 $L$ 인 판 모양 영역에서 이온이 전자보다 비율 $\delta$ 만큼 많다고 하자.

\frac{d^2\phi}{dx^2} = -\frac{e(n_i - n_e)}{\varepsilon_0} = -\frac{e\,n\,\delta}{\varepsilon_0}\quad\Longrightarrow\quad \phi_{\max} = \frac{e\,n\,\delta\,L^2}{8\varepsilon_0},\qquad E_{\max} = \frac{e\,n\,\delta\,L}{2\varepsilon_0}

양쪽 벽을 접지(

\phi = 0

)로 두었다. 이 전위가 전자의 열에너지

T_e

(V 단위)를 넘으면 전자가 그 영역으로 끌려와 불균형을 메운다.

\phi_{\max} = T_e

로 두면 허용되는 최대 불균형은

\delta^* = 8\varepsilon_0 T_e/(e n L^2) = 8(\lambda_D/L)^2

이다.

PL-300 IC( $n = 4\times10^{17}$ m⁻³)에서 지름 10 cm 영역의 전하가 $\delta = 10^{-6}$ 만 어긋나도 $\phi_{\max} \approx 9$ V가 생긴다. 전자 온도 2.5 eV의 약 3.6배다. 그러니 벌크에서 $\delta$ 는 $\delta^* \approx 3\times10^{-7}$ 보다 작아야 하고, 실제로 그렇다. 준중성은 가정이 아니라 큰 전기장이 강제하는 결과다. 이 규칙이 깨지는 곳은 $L$ 이 $\lambda_D$ 와 비슷해지는 곳, 곧 벽 앞의 시스뿐이다(7장).

그림 2-1. 벽 사이 플라즈마의 밀도와 전위. 벌크에서는 두 밀도가 그림 두께 안에서 겹친다. 전자가 이온보다 훨씬 빨라 벽으로 먼저 빠지므로, 벽 앞 수십 λ_D에 이온이 남는 시스가 생기고 플라즈마 전위는 벽보다 몇 Tₑ 높아진다.

SIMULATOR

전하 불균형 δ가 만드는 전위: 준중성의 크기

밀도 n 전자 온도 Tₑ 영역 크기 L 불균형 δ = (nᵢ−nₑ)/n

φ_max—

E_max—

허용 불균형 δ* = 8(λ_D/L)²—

판정—

곡선은 영역 크기 L, L/10, L/100에서 δ에 따른 φ_max이고, 수평 점선은 Tₑ다. 점선 아래(φ_max < Tₑ)에서만 그 불균형이 버틸 수 있다. 해볼 것: ① 기본값(10 cm, δ = 10⁻⁶)에서 φ_max ≈ 9 V > Tₑ. ② L을 20 µm(≈ λ_D)로 줄이면 δ ~ 1이어도 된다 — 데바이 길이보다 작은 영역에서는 준중성이 깨진다. ③ 밀도를 높일수록 허용 δ*가 작아진다(δ* ∝ 1/n).

데바이 차폐: 포아송-볼츠만 식

플라즈마 속에 전극을 넣어 전위 $\phi_0$ 를 걸면, 전위는 전극에서 멀어지며 얼마나 빨리 0으로 돌아올까? 진공이라면 판 전극의 전위는 거리에 따라 선형으로만 변한다. 플라즈마에서는 전자가 움직여 전극을 감싸고 전위를 가린다. 이것이 데바이 차폐(Debye shielding)다. 전자는 전위 속에서 볼츠만 분포를 따르고, 무거운 이온은 이 짧은 시간 동안 움직이지 않는다고 보면( $n_i = n_0$ ) 다음 식을 얻는다.

\varepsilon_0\frac{d^2\phi}{dx^2} = e\,n_0\left[\exp\!\left(\frac{\phi}{T_e}\right) - 1\right]\quad\xrightarrow{\ \psi = \phi/T_e,\ \xi = x/\lambda_D\ }\quad \frac{d^2\psi}{d\xi^2} = e^{\psi} - 1

T_e

는 V 단위,

n_e = n_0 e^{\phi/T_e}

(볼츠만 관계). 길이를

\lambda_D = \sqrt{\varepsilon_0T_e/en_0}

로 재면 식에서 모든 상수가 사라진다. 데바이 길이는 이 식이 고른 자연 단위다.

$|\psi| \ll 1$ 이면 $e^\psi - 1 \approx \psi$ 로 선형화해 $\psi = \psi_0 e^{-\xi}$ , 곧 $\phi = \phi_0 e^{-x/\lambda_D}$ 를 얻는다. 판 대칭에서는 비선형 식을 한 번 적분할 수 있다. 양변에 $\psi'$ 를 곱해 적분하면 $\tfrac12(\psi')^2 = e^{\psi} - 1 - \psi$ (멀리서 $\psi, \psi' \to 0$ )이고, 이 1계 식 $\psi' = \mp\sqrt{2(e^\psi - 1 - \psi)}$ 를 적분하면 기준 해를 얻는다. 구 대칭(점전하·작은 탐침)에서는 선형해가 유카와 꼴 $\phi = \phi_0(a/r)e^{-(r-a)/\lambda_D}$ 가 된다. 아래 시뮬레이터는 비선형 식을 격자 600개에서 뉴턴법으로 직접 풀어 두 근사와 비교한다.

\frac{\psi_{i+1} - 2\psi_i + \psi_{i-1}}{h^2} + \frac{g}{\xi_i}\frac{\psi_{i+1} - \psi_{i-1}}{2h} - \left(e^{\psi_i} - 1\right) = 0,\qquad J_{ii} = -\frac{2}{h^2} - e^{\psi_i}

g = 0

(판) 또는 2(구). 경계는

\psi(\xi_0) = \psi_0

\psi(\xi_{\max}) = 0

. 야코비 행렬이 삼중 대각이라 토머스 알고리즘으로 한 번에 풀고, 한 번에 바뀌는 양을 2 이하로 묶어(감쇠 뉴턴) 수렴시킨다.

SIMULATOR

1D 포아송-볼츠만 풀이: 비선형 차폐 vs 선형 지수

전극 전위 eφ₀/Tₑ

기하

전자 온도 Tₑ 밀도 n₀

λ_D—

φ가 1/e로 줄어드는 거리—

뉴턴 반복—

기준 해와 최대 차이—

위: 전위 ψ = eφ/Tₑ (실선: 뉴턴 수치 해, 점선: 선형 근사, 점: 판 1차 적분 기준 해). 아래: 알짜 전하 (nᵢ − nₑ)/n₀ — 양수는 전자가 밀려난 이온 과잉, 음수는 전자가 몰린 곳. 해볼 것: ① ψ₀ = −0.5처럼 작으면 세 곡선이 겹친다. ② ψ₀ = −10이면 전자가 벽 앞 수 λ_D에서 사실상 0이 된다(알짜 전하 → 1). 비선형 해는 선형보다 멀리 뻗는다 — 시스의 씨앗이다. ③ ψ₀ = +5로 전자를 끌어당기면 전자가 벽에 몰려 0.1–0.3 λ_D 안에서 대부분 가려지고 선형 근사는 크게 틀린다. ④ 구 기하에서는 1/r 감쇠가 겹쳐 더 빨리 가려진다.

그림 2-2. 데바이 차폐. 양전하 주위에 전자가 몰려(이온은 약간 밀려나) 전위가

1/r

보다 훨씬 빨리 줄어든다. 차폐에 참여하는 것은 열운동하는 전자이므로, 전자가 뜨거울수록(

T_e

↑) 잘 모이지 않아

\lambda_D

가 길어지고, 밀도가 높을수록 짧아진다.

그림의 직관을 숫자로 보자. PL-300 IC 벌크에서 $\lambda_D \approx 7430\sqrt{2.5/4\times10^{17}} \approx 19\ \mu\text{m}$ 다. 챔버 반지름(0.25 m)은 데바이 길이의 약 13,000배이므로 벌크는 완벽하게 준중성이다. 반면 랭뮤어 탐침(반지름 0.1 mm)은 $\lambda_D$ 의 다섯 배밖에 안 되어 탐침 주위 시스가 탐침 크기와 비슷하다. 탐침 이론이 까다로운 이유다(15장).

플라즈마 주파수: 전자는 얼마나 빨리 출렁이는가

준중성을 강제하는 전기장은 복원력이다. 전자 한 층을 이온 배경에 대해 $\xi$ 만큼 밀면, 양쪽에 면전하 $\pm e n \xi$ 가 생겨 평행판 축전기처럼 전기장 $E = e n\xi/\varepsilon_0$ 가 걸리고 전자를 되돌린다. 용수철과 같은 꼴이다.

m_e\frac{d^2\xi}{dt^2} = -eE = -\frac{n e^2}{\varepsilon_0}\xi\quad\Longrightarrow\quad \omega_{pe} = \sqrt{\frac{n e^2}{\varepsilon_0 m_e}},\qquad f_{pe} \approx 8.98\sqrt{n[\text{m}^{-3}]}\ \text{Hz}

PL-300 IC(

n = 4\times10^{17}

m⁻³):

f_{pe} \approx 5.7

GHz, 주기 약 0.18 ns. 이온도 같은 식을 따르되 질량이

M

이므로

\omega_{pi} = \omega_{pe}\sqrt{m_e/M}

— Ar에서

f_{pi} \approx 21

MHz다. 전자에 압력(열운동)이 더해지면 파수

k

인 진동은 보옴-그로스 분산

\omega^2 = \omega_{pe}^2 + 3k^2 T_e e/m_e

를 따르고, 기체 충돌(빈도

\nu_m

)은 진동을

e^{-\nu_m t/2}

로 감쇠시킨다.

그림 2-3. 전자 층을 ξ만큼 밀면 평행판 축전기가 생긴다. 복원력이 ξ에 비례하므로 전자는 ω_pe로 진동한다. 진동수가 층 두께나 ξ와 무관하고 밀도만의 함수라는 점이 핵심이다.

SIMULATOR

전자 층 진동: 플라즈마 주파수, 보옴-그로스, 충돌 감쇠

밀도 n 압력 (Ar) 파장 / λ_D 전자 온도 Tₑ

f_pe—

보옴-그로스 f—

충돌 빈도 ν_m / ω_pe—

진폭 1/e 감쇠—

위: 이온 배경(옅은 띠) 위 전자 층 24장의 위치와 알짜 전하(붉은 쪽 이온 과잉, 푸른 쪽 전자 과잉). 아래: 한 층의 변위 ξ(t). 화면은 플라즈마 주기 하나를 약 1.5초로 늦춰 보여 준다. 각 층은 RK4로

\ddot\xi = -\omega^2\xi - \nu_m\dot\xi

를 적분한다. 해볼 것: ① 10 mTorr에서 ν_m/ω_pe ≈ 10⁻³ — 진동이 수백 주기 넘게 산다. ② 10 Torr로 올리면 ν_m이 ω_pe에 가까워져 몇 주기 안에 죽는다(고압에서 플라즈마 조건 ③이 위태로운 이유). ③ 파장을 λ_D의 몇 배로 줄이면 열운동 압력이 복원력에 더해져 주파수가 오른다.

공정 RF(0.4–60 MHz)는 $f_{pe}$ 보다 두세 자릿수 느리다. 그래서 전자는 RF 전기장에 거의 즉시 따라가 플라즈마 안의 전기장을 가린다. RF 에너지는 시스(CCP)나 표피층(ICP)에서만 들어간다. 이온은 사정이 다르다. 2 MHz는 $f_{pi}$ 보다 느려 이온이 시스 전압의 순간값을 따라가고, 60 MHz는 빨라 평균값만 느낀다. 이중 주파수 CCP의 바이어스 주파수가 낮고 소스 주파수가 높은 이유다(9장, 이온 에너지 분포 폭은 7장).

SIMULATOR

주파수 사다리: 누가 RF를 따라가는가

밀도 n 압력 (Ar, 300 K) 전자 온도 Tₑ

밀도 위치

f_pi (이온)—

f_pe (전자)—

ν_en/2π (전자-중성)—

ν_in/2π (이온-중성)—

로그 주파수 축 위에 플라즈마 고유 주파수(색 막대)와 공정 RF 주파수(세로선)를 놓았다. 아래 표는 RF 주파수마다 이온과 전자가 따라가는지(f_RF < f_p)를 판정한다. 충돌 빈도는

\nu_{en} = n_gK_{el}(T_e)

\nu_{in} = n_g\sigma_i\bar v_i

(σᵢ = 10⁻¹⁸ m², 300 K)로 계산한다. 해볼 것: ① 시스 경계 밀도로 바꾸면 f_pi가 약 9 MHz로 떨어져 13.56 MHz도 "이온이 못 따라가는" 쪽에 들어간다. ② 압력을 올리면 ν_in이 f_pi를 넘는다 — 이온은 충돌에 막혀 RF 순간값을 따라가지 못한다.

플라즈마 파라미터와 결합 파라미터

데바이 차폐는 통계적 현상이다. 전하를 가리는 전자 구름이 "구름"이려면 반지름 $\lambda_D$ 인 구 안에 전자가 아주 많아야 한다. 그 수를 플라즈마 파라미터(plasma parameter) $N_D$ 라 한다. 같은 내용을 입자 사이 평균 거리 $a$ 에서의 쿨롱 에너지와 열에너지의 비, 결합 파라미터(coupling parameter) $\Gamma$ 로도 쓸 수 있다.

N_D = \frac{4\pi}{3}n\lambda_D^3 \approx 1.72\times10^{12}\,\frac{T_e[\text{eV}]^{3/2}}{\sqrt{n[\text{m}^{-3}]}},\qquad \Gamma = \frac{e^2/(4\pi\varepsilon_0 a)}{e\,T_e},\quad a = \left(\frac{3}{4\pi n}\right)^{1/3},\qquad \Gamma \approx \frac{1}{3}N_D^{-2/3}

N_D \gg 1

(곧

\Gamma \ll 1

)이면 약결합 플라즈마로, 입자 하나하나의 쿨롱 충돌보다 수많은 먼 입자의 집단 전기장이 운동을 지배한다. 이온-전자 쿨롱 충돌 빈도와 플라즈마 주파수의 비도 대략

\nu_{ei}/\omega_{pe} \sim \ln\Lambda/N_D

로 작다. PL-300 IC:

N_D \approx 1.1\times10^4

\Gamma \approx 7\times10^{-4}

그림 2-4. 약결합과 강결합. 공정 플라즈마는 데바이 구 안에 전자가 10³–10⁵개 있는 전형적인 약결합 플라즈마다. 예외로, 기상에서 자란 µm 크기 먼지 입자는 전하가 수천 e라서 Γ ≫ 1인 "먼지 플라즈마 결정"을 이룰 수 있다(13장).

SIMULATOR

데바이 구 안의 전자 수와 결합 파라미터

밀도 n 전자 온도 Tₑ

N_D—

평균 거리 a / λ_D—

결합 Γ—

판정—

왼쪽: λ_D 구의 단면. 구 안 전자 N_D개를 그리되 600개가 넘으면 600개만 그린다(점 하나가 여러 개를 대표). 오른쪽: (n, Tₑ) 평면에서 N_D = 1, 10³, 10⁶ 선과 Γ = 1 선. 해볼 것: ① PL-300 기본값에서 N_D ≈ 10⁴. ② 밀도를 10²⁴ cm⁻³ 쪽으로, 온도를 0.1 eV 아래로 보내면 N_D < 1 — 구 안이 텅 빈다. ③ Γ = 1 선과 N_D ≈ 0.2 선이 거의 겹친다(Γ ≈ N_D^(−2/3)/3).

전리도와 사하 식: 공정 플라즈마는 평형에서 얼마나 먼가

전리도(ionization degree) $x = n_e/(n_e + n_g)$ 는 기체 가운데 이온화된 비율이다. PL-300 IC는 $x \approx 4\times10^{17}/3.2\times10^{20} \approx 1.3\times10^{-3}$ 이다. 기체가 한 온도 $T$ 의 열평형에 있다면 전리도는 이온화와 재결합의 균형인 사하 식(Saha equation)으로 정해진다.

\frac{n_e n_i}{n_0} = \frac{2g_i}{g_0}\left(\frac{2\pi m_e kT}{h^2}\right)^{3/2}\exp\!\left(-\frac{U_{iz}}{kT}\right)\approx \frac{2g_i}{g_0}\,2.415\times10^{21}\,T^{3/2}e^{-U_{iz}/kT}\ \ [\text{m}^{-3}]

U_{iz}

는 이온화 에너지(Ar 15.76 eV),

g

는 통계 무게(Ar에서

2g_i/g_0 \approx 12

, 근사),

T

는 K. 일정 압력

p

에서

n_0 + n_i = p/kT

로 두고

x^2/(1-x) = S(T)

를 푼다. 지수 인자 때문에

kT \ll U_{iz}

이면

x

는 천문학적으로 작다.

공정 플라즈마에 사하 식을 적용하면 모순이 드러난다. 기체 온도 300 K를 넣으면 $x \sim 10^{-130}$ , 전자 온도 29,000 K를 넣으면 $x \approx 1$ (완전 전리)이다. 실제 $10^{-3}$ 은 어느 쪽과도 맞지 않는다. 공정 플라즈마의 전리도는 열평형이 아니라 이온화 속도와 벽 손실의 균형(입자 균형, 8장)으로 정해지기 때문이다. 전자는 벽으로 빠져나가 재결합할 틈이 없고, 기체는 차가운 채로 남는다.

그림 2-5. 이온화는 전자 에너지 분포의 꼬리가 한다. Tₑ = 2.5 eV에서 평균 에너지(3.75 eV)는 문턱의 1/4이지만, 15.76 eV를 넘는 전자가 약 0.6 % 있어 이온화를 유지한다. 꼬리의 모양이 반응 속도를 얼마나 바꾸는지는 4장, 속도 계수 계산은 3장.

SIMULATOR

사하 평형 전리도 vs 공정 플라즈마

기체

압력 실제 nₑ 실제 Tₑ

실제 전리도 (T_g = 300 K)—

사하 x at 300 K—

사하 x at Tₑ—

실제 x를 주는 평형 온도—

곡선은 일정 압력에서 사하 식의 평형 전리도 x(T), 수평 점선은 실제 전리도, 세로선은 실제 Tₑ다. 해볼 것: ① Ar 10 mTorr, 4×10¹¹ cm⁻³에서 실제 x ≈ 1.2×10⁻³을 평형으로 얻으려면 약 5,400 K가 필요하다. 사하 곡선은 Tₑ(29,000 K)에서 이미 1에 붙어 있다. ② Cs나 Na처럼 이온화 에너지가 작은 알칼리는 2,000 K 대기압 불꽃에서도 x ~ 10⁻⁶–10⁻⁴이다. 불꽃 플라즈마의 전자는 대부분 이런 불순물에서 나온다. ③ 압력을 대기압으로 올리면 같은 온도에서 x가 줄어든다(재결합 상대가 많아져서).

IC 모듈의 전리도와 "온도"Ar 10 mTorr(

n_g = 3.2\times10^{20}

m⁻³),

n_e = 4\times10^{17}

m⁻³ →

x \approx 1.2\times10^{-3}

. 이온은 10⁶개의 중성 원자 가운데 약 1,200개, 플라즈마에 쓰인 전력 1.5 kW 가운데 기체 가열로 가는 몫은 일부다. 그래서 기체 온도는 300–600 K에 머무르고, 이온 온도는 이온-중성 전하 교환 때문에 기체와 비슷하다(약 0.03–0.05 eV). 사하 식으로 이 전리도를 "설명"하려면 5,400 K가 필요하지만, 그런 온도는 챔버 어디에도 없다. 이 플라즈마는 세 온도(Tₑ ≈ 29,000 K, Tᵢ ≈ T_g ≈ 수백 K, 평형 등가 5,400 K)가 공존하는 비평형 계다.

전자기파 차단: 어떤 주파수가 플라즈마를 뚫는가

전자가 전기장을 가리는 성질은 전자기파에도 똑같이 적용된다. 파동의 주파수 $\omega$ 가 $\omega_{pe}$ 보다 낮으면 전자가 따라 움직여 파동을 반사하고, 높으면 전자가 따라가지 못해 파동이 통과한다. 시간 의존을 $e^{-i\omega t}$ 로 두면 충돌이 드문 플라즈마의 유전율과 분산 관계는 다음과 같다.

\varepsilon_p(\omega) = 1 - \frac{\omega_{pe}^2}{\omega(\omega + i\nu_m)},\qquad \omega^2 = \omega_{pe}^2 + c^2k^2,\qquad n_c = \frac{\varepsilon_0 m_e\omega^2}{e^2} \approx 1.24\times10^{10}\left(\frac{f}{1\,\text{GHz}}\right)^2\ \text{cm}^{-3}

n_e \gt n_c

인 과밀(overdense) 플라즈마에서는

k

가 허수가 되어 파동이 표피 깊이

\delta = c/\sqrt{\omega_{pe}^2 - \omega^2}

(무충돌,

\approx c/\omega_{pe}

) 안에서 감쇠한다. 2.45 GHz의 차단 밀도는

7.4\times10^{10}

cm⁻³, 13.56 MHz는

2.3\times10^{6}

cm⁻³다.

그림 2-6. 플라즈마 속 전자기파의 분산 관계(왼쪽)와 차단 밀도(오른쪽). 공정 RF는 모든 공정 플라즈마에서 과밀이라 플라즈마 속으로 전파하지 못한다. 마이크로파 간섭계는 차단 밀도보다 훨씬 높은 주파수를 써서 파동이 통과하며 겪는 위상 변화로 밀도를 잰다(15장).

아래 시뮬레이터는 밀도가 경사를 이루며 올라가는 플라즈마에 왼쪽에서 평면파를 쏘아, 1D 파동 방정식 $E'' + k_0^2\varepsilon_p(x)E = 0$ 을 복소수 RK4로 오른쪽 끝에서 거꾸로 적분한다. 오른쪽 끝에서는 "나가는 파동 또는 감쇠하는 파동"만 있다는 조건을 주고, 왼쪽 진공 영역에서 입사파 $A e^{ik_0x}$ 와 반사파 $Be^{-ik_0x}$ 로 나눠 반사율 $|B/A|^2$ 을 구한다.

SIMULATOR

밀도 경사로 들어가는 전자기파: 차단과 표피 깊이

파동 주파수 f 최대 밀도 n_max 압력 (충돌)

주파수 바로 가기

차단 밀도 n_c—

n_max / n_c—

반사 / 투과 / 흡수—

표피 깊이 δ—

위: 밀도 분포(1–4 cm에서 0 → n_max, 그 뒤 평탄)와 차단 밀도(점선). 아래: |E|/|E_입사|(굵은 선)와 순간 전기장(가는 선). 충돌 빈도는 ν_m = n_g K_el(3 eV). 해볼 것: ① 2.45 GHz, n_max = 3×10¹⁰ cm⁻³(n_c의 40 %)는 거의 다 통과한다. ② 밀도를 10¹¹로 올리면 차단층에서 전반사되고 그 너머는 수 mm 표피 깊이로 감쇠한다 — 2.45 GHz 표면파 소스가 과밀 플라즈마를 유지하는 원리는 10장. ③ 13.56 MHz는 어떤 공정 밀도에서도 차단이다. 압력을 올려 충돌을 늘리면 흡수가 생긴다(충돌 표피 가열).

조건	n_max / n_c	결과	공정에서의 의미
2.45 GHz, nₑ = 3×10¹⁰ cm⁻³	0.4	대부분 투과, 약간 반사	점화 직후: 마이크로파가 플라즈마를 통과해 챔버 전체를 데움
2.45 GHz, nₑ = 4×10¹¹ cm⁻³	5.4	전반사, 표피 깊이 약 9 mm	유지 상태: 전력은 창 아래 표피층과 표면파로 흡수(10장)
헤어핀 탐침 1.5–6 GHz	공진 주파수 $f_r^2 = f_0^2 + f_{pe}^2$	nₑ 약 10⁹–4×10¹¹ cm⁻³ 측정	PL-300 진단(15장). 6 GHz면 IC 기준 밀도까지 잴 수 있음
13.56 MHz, 모든 공정 밀도	10⁴–10⁶	차단, 표피 깊이 수 mm–cm	ICP: 코일 전력은 창 아래 표피층에서만 흡수. CCP: 전력은 시스로 들어감
94 GHz 간섭계, nₑ ~10¹¹–10¹²	~10⁻³–10⁻²	투과, 위상 지연	선적분 밀도 측정(15장)

플라즈마의 세 조건과 척도 계산기

이 장의 숫자들을 모으면 이온화 기체가 플라즈마로 행동하기 위한 세 조건이 된다. 셋째 조건은 진동이 충돌로 지워지지 않아야 한다는 뜻인데, 수 Torr 방전에서는 위태롭다. 그래도 그런 방전을 플라즈마라고 부르는 것은 ①과 ②가 충분히 성립해 준중성·차폐·시스 같은 플라즈마 구조가 그대로 나타나기 때문이다.

\text{①}\ \lambda_D \ll L,\qquad \text{②}\ N_D \gg 1,\qquad \text{③}\ \omega_{pe}\tau_{en} = \frac{\omega_{pe}}{\nu_{en}} \gt 1

L

은 플라즈마 크기,

\nu_{en} = n_gK_{el}(T_e)

는 전자-중성 운동량 전달 충돌 빈도. 계산기는 여기에 평균 자유 행로

\lambda_e = \bar v_e/\nu_{en}

\lambda_i = 1/(n_g\sigma_i)

와 전리도를 함께 보여 준다(기체 온도 300 K, Ar).

그림 2-7. 플라즈마의 세 조건. ①은 공간, ②는 입자 수, ③은 시간에 관한 조건이다. 공정 플라즈마는 ①과 ②를 수천–수만 배 여유로 만족하고, ③은 수 mTorr에서는 1,000배 여유가 있지만 수 Torr PECVD에서는 1 안팎이다.

SIMULATOR

척도 계산기: nₑ, Tₑ, p, L → 플라즈마인가?

전자 밀도 nₑ 전자 온도 Tₑ 압력 (Ar, 300 K) 플라즈마 크기 L

프리셋

λ_D—

f_pe—

N_D—

결합 Γ—

λₑ / λᵢ—

전리도 x—

① L/λ_D —

② N_D —

③ ω_pe/ν_en —

막대는 로그 길이 축 위의 척도들이다: 입자 간 거리 a, λ_D, 부유 시스 두께(≈ 5λ_D, 어림), 이온·전자 평균 자유 행로, 플라즈마 크기 L. 불꽃 프리셋은 2,000 K 대기압과 같은 기체 밀도(300 K로 환산해 약 114 Torr)를 넣었다. 해볼 것: ① PL-300 IC에서 L/λ_D ≈ 1.3×10⁴, N_D ≈ 10⁴, ω_pe/ν_en ≈ 1,100 — 세 조건 모두 큰 여유. ② PE PECVD(3 Torr)는 ③이 약 0.8로 경계다. ③ 불꽃은 ③이 0.5 아래로, 충돌이 지배하는 "약한" 플라즈마다. ④ 전리층은 크기가 km라 ①이 압도적이다.

그림 2-8. PL-300 IC 기준 플라즈마의 길이 사다리. 데바이 길이 → 시스 → 이온 자유 행로 → 전자 자유 행로 → 챔버 순으로 각각 한 자릿수쯤씩 커진다. 시스가 λᵢ보다 얇아 이온이 무충돌로 웨이퍼에 닿고(7장), λₑ가 챔버와 비슷해 저압 ICP의 전자 가열이 비국소적이다(4장).

플라즈마 (조건)	λ_D	f_pe	N_D	ω_pe/ν_en	전리도	판정
PL-300 IC: Ar 10 mTorr, 2.5 eV, 4×10¹¹ cm⁻³	19 µm	5.7 GHz	1.1×10⁴	~1,100	1.2×10⁻³	세 조건 모두 큰 여유
PL-300 CP: 30 mTorr, 3 eV, 5×10¹⁰	58 µm	2.0 GHz	4×10⁴	~100	5×10⁻⁵	간격 25 mm/λ_D ≈ 430
PL-300 PE: 3 Torr, 2 eV, 10¹⁰	0.1 mm	0.9 GHz	5×10⁴	~0.8	10⁻⁷	③ 경계, 간격 10 mm/λ_D ≈ 95
마그네트론: 3 mTorr, 3 eV, 10¹²	13 µm	9 GHz	9×10³	~4,600	10⁻²	큰 여유, 전리도 높음
불꽃: 대기압 2,000 K, 0.2 eV, 10¹⁰	33 µm	0.9 GHz	1.5×10³	~0.5	~10⁻⁹	③ 불만족, 충돌 지배

세 척도로 본 IC 모듈 설계창-웨이퍼 간격 10 cm는 λ_D의 5,000배, λᵢ의 30배, λₑ의 3배다. 그래서 벌크는 준중성(①), 이온은 벌크에서 여러 번 충돌하며 양극성 확산(6장)하지만 시스(0.1–0.7 mm)는 무충돌로 건넌다. 전자는 간격을 한 번 건너는 동안 몇 번밖에 충돌하지 않고 탄성 충돌로는 에너지를 거의 잃지 않으므로 창 아래 표피층(수 cm)에서 받은 에너지를 챔버 전체로 나른다. 헤어핀 탐침의 상한 6 GHz는

n_e = 4.5\times10^{11}

cm⁻³ 플라즈마의 f_pe와 비슷하다. 공진 주파수는

\sqrt{f_0^2 + f_{pe}^2}

라 탐침 설계

f_0

를 함께 고려해야 기준 밀도를 잴 수 있다.

핵심 정리

전하 불균형 δ는 전위 $\phi \approx e n\delta L^2/8\varepsilon_0$ 을 만든다. $\phi \lt T_e$ 이려면 $\delta \lt 8(\lambda_D/L)^2$ — 벌크는 $10^{-7}$ 수준까지 준중성이다.
볼츠만 전자를 넣은 포아송 식은 무차원 꼴 $\psi'' = e^\psi - 1$ 이고 길이 단위가 $\lambda_D = \sqrt{\varepsilon_0T_e/en}$ 다. $|\psi| \ll 1$ 이면 전위가 $e^{-x/\lambda_D}$ 로 가려지고, 큰 음전위에서는 전자가 고갈된 시스가 생긴다.
$\omega_{pe} = \sqrt{ne^2/\varepsilon_0m_e}$ , $f_{pe} \approx 8.98\sqrt{n}$ Hz. PL-300 IC: 5.7 GHz. 이온은 $f_{pi} = f_{pe}\sqrt{m_e/M}$ (Ar 약 21 MHz). 공정 RF는 전자에게는 느리고, 이온에게는 주파수에 따라 느리거나(2 MHz) 빠르다(60 MHz).
$N_D = \tfrac{4\pi}{3}n\lambda_D^3 \approx 10^4$ , $\Gamma \approx 7\times10^{-4}$ : 공정 플라즈마는 약결합 이상 기체다.
사하 식은 공정 플라즈마의 전리도( $10^{-3}$ )를 설명하지 못한다(300 K → $10^{-130}$ , Tₑ → 1). 전리도는 이온화-벽 손실 균형이 정하는 비평형 값이다.
차단 밀도 $n_c \approx 1.24\times10^{10}(f/\text{GHz})^2$ cm⁻³. 과밀 플라즈마에서 파동은 표피 깊이 $\approx c/\omega_{pe}$ 안에서 감쇠한다. 13.56 MHz는 늘 과밀, 2.45 GHz는 $7.4\times10^{10}$ cm⁻³에서 차단된다.
플라즈마의 세 조건: $\lambda_D \ll L$ , $N_D \gg 1$ , $\omega_{pe}/\nu_{en} \gt 1$ . PL-300 IC 길이 사다리는 λ_D 19 µm → 시스 0.1–0.7 mm → λᵢ 3 mm → λₑ 3.4 cm → 챔버 10–25 cm.

확인 퀴즈

1. 전자 온도는 그대로 두고 전자 밀도를 4배로 올리면 데바이 길이는?

\lambda_D \propto \sqrt{T_e/n}

이므로 n이 4배면 λ_D는 1/2이다. 차폐에 참여하는 전자가 많을수록 더 짧은 거리에서 전위를 가린다.

2. PL-300 IC 플라즈마(nₑ = 4×10¹¹ cm⁻³)의 전자 플라즈마 주파수에 가장 가까운 것은?

f_{pe} \approx 8.98\sqrt{4\times10^{17}} \approx 5.7\times10^9

Hz. 13.56 MHz보다 약 400배 빠르므로 전자는 RF 전기장을 거의 즉시 가린다.

3. 포아송-볼츠만 식을 선형화한 해 φ₀e^(−x/λ_D)가 가장 크게 틀리는 경우는?

선형화는

e^\psi - 1 \approx \psi

이므로 |ψ| ≪ 1에서만 맞는다. ψ = +5이면 e⁵ ≈ 148배로 전자가 몰려 훨씬 짧은 거리에서 가려진다. 밀도는 λ_D 단위로 재면 식에서 사라지므로 선형화의 정확도와 무관하다.

4. 이중 주파수 CCP에서 2 MHz 바이어스가 이온 에너지를 크게 흔들고 60 MHz 소스는 그렇지 않은 이유로 가장 알맞은 것은?

이온은

f_{pi} = f_{pe}\sqrt{m_e/M}

보다 느린 전기장만 따라간다. 시스 영역의 밀도에서 f_pi는 수 MHz–10 MHz 정도라 2 MHz 바이어스의 순간 전압을 따라가 넓은 두 봉우리 IEDF가 되고, 60 MHz에서는 시간 평균 전압만 느낀다(7장, 9장).

5. 기체 온도 300 K인 Ar 10 mTorr 플라즈마의 전리도가 10⁻³인 것을 사하 식으로 해석할 때 옳은 것은?

300 K 사하는 10⁻¹³⁰, Tₑ 사하는 거의 1이다. 저압 플라즈마에서는 전자-이온이 부피에서 재결합하기 전에 벽으로 빠져나가므로 평형이 성립하지 않는다. 실제 전리도는 8장의 입자·전력 균형으로 구한다.

6. nₑ = 1×10¹¹ cm⁻³인 플라즈마에 2.45 GHz 마이크로파를 수직으로 쏘았다(차단 밀도 7.4×10¹⁰ cm⁻³). 무엇이 일어나는가?

n > n_c이므로 플라즈마 안에서 k가 허수가 된다. 파동은

n = n_c

인 곳에서 반사되고 그 너머로는

\delta \approx c/\sqrt{\omega_{pe}^2 - \omega^2}

(수 mm–cm)만큼 스며든다. 과밀 마이크로파 플라즈마는 표면파로 전력을 넣는다.

7. 다음 중 플라즈마의 셋째 조건(ω_pe/ν_en > 1)을 가장 만족하기 어려운 것은?

ν_en은 기체 밀도에 비례한다. 3 Torr에서 ν_en은 10 mTorr의 300배이고 밀도(→ ω_pe)는 오히려 낮아 비가 1 안팎이다. 그래도 λ_D ≪ L, N_D ≫ 1은 충분히 만족해 시스와 준중성 같은 플라즈마 구조는 그대로 나타난다.