확산은 원자의 무작위 열운동과 결함을 통한 점프에서 시작합니다. 그러나 우리가 관찰하는 알짜 흐름은 무작위가 아니라 계의 자유에너지를 낮추는 방향을 가집니다.
우리 시편의 전체 조성은 35 wt% B지만 상태도는 공존하는 액체가 32 wt% B, α 고체가 43 wt% B여야 한다고 말합니다. 새 고체가 생길 때마다 B 원자는 계면 주변에서 다시 분배되어야 합니다.
단단해지는 시편 안에서 B 원자는 어떻게 제자리를 찾아갈까요?
이번 장의 답원자는 공공·침입 자리를 통해 점프하고, 화학퍼텐셜 기울기가 만든 알짜 플럭스가 조성 분배를 진행시킵니다.
01고체에서 원자는 어떻게 움직이는가
결정 속 원자는 평형 위치 주변에서 진동하다가 충분한 열에너지를 얻으면 이웃 자리로 점프합니다. 치환형 원자는 보통 빈 격자점인 공공(vacancy)과 자리를 바꾸고, C·N·H처럼 작은 원자는 격자 사이의 침입 자리(interstitial site) 사이를 이동합니다.
- 공공 기구 — 공공의 형성 확률과 원자–공공 교환의 이동 장벽이 모두 확산속도에 영향을 줍니다.
- 침입형 기구 — 이동 가능한 빈 침입 자리가 많아 대체로 치환형 확산보다 빠릅니다.
- 교환·링 기구 — 여러 원자가 협동해 이동하는 경로도 가능하지만 보통 높은 에너지가 필요합니다.
개별 점프의 방향은 열적 요동 때문에 무작위입니다. 1차원 무작위 걸음에서 평균 변위는 0이지만 평균제곱변위는 시간에 비례합니다.
$$ \langle x\rangle=0,\qquad \langle x^2\rangle=2Dt $$
이 식의 $D$는 추적 원자 하나의 무작위 이동을 나타내는 자기·추적 확산계수와 연결됩니다. 조성 기울기가 있는 합금의 상호확산은 열역학적 상호작용과 두 성분의 이동성을 함께 포함하므로 같은 개념으로 무조건 치환하면 안 됩니다.
02진짜 구동력은 화학퍼텐셜 기울기다
“높은 농도에서 낮은 농도로 간다”는 설명은 이상적인 묽은 용액에서는 맞지만 일반적인 법칙은 아닙니다. 비이상 합금, 외부 전기장이나 응력장이 있는 계에서 원자는 농도보다화학퍼텐셜 $\mu$을 낮추는 방향으로 이동합니다.
$$ \mathbf{J}=-L\,\nabla\mu $$
$L$은 이동성을 포함한 양의 수송계수입니다. 이상적인 묽은 용액에서는 화학퍼텐셜 기울기가 농도 기울기에 비례해 익숙한 Fick 형태로 줄어듭니다. 따라서 역확산처럼 겉보기에는 낮은 농도에서 높은 농도로 움직이는 경우도 열역학적으로 가능하며, 이때도 화학퍼텐셜은 낮아지는 방향입니다.
03픽의 1법칙 — 농도 기울기와 순간 플럭스
Mehrer 교재의 연속체 설명처럼, 등방성 매질에서 가장 단순한 1차원 Fick의 1법칙은 다음과 같습니다.
$$ J_x=-D\frac{\partial C}{\partial x} $$
$J_x$는 단위 면적·단위 시간당 지나가는 입자 또는 몰의 양, $C$는 농도, $D$는 확산계수로 단위는 $\mathrm{m^2/s}$입니다. 음의 부호는 흐름과 농도 증가 방향이 반대라는 뜻입니다.
1법칙은 정상상태에서 특히 간단히 활용되지만 정상상태에만 한정된 식은 아닙니다. 비정상 확산에서도 각 순간의 국소 플럭스와 농도 기울기를 연결하며, 시간 변화는 이어지는 보존식과 2법칙으로 구합니다.
이 단순식은 등방성, 단일한 확산종, 농도 기울기가 주된 구동력이라는 가정을 담고 있습니다. 결정이 이방성이면 $D$는 방향에 따라 달라지는 텐서가 되고, 반응·전기장·응력·빠른 확산경로가 있으면 추가 항이나 더 일반적인 화학퍼텐셜 표현이 필요합니다.
04질량보존과 픽의 2법칙
작은 부피에 들어오는 양과 나가는 양의 차이는 그 안의 농도 변화가 됩니다. 내부 생성·소멸이 없는 보존식 $\partial C/\partial t=-\nabla\cdot\mathbf{J}$에 Fick의 1법칙을 넣으면 일반적인 확산 방정식을 얻습니다.
$$ \frac{\partial C}{\partial t}=\nabla\cdot(D\nabla C) $$
1차원에서 $D$가 조성·위치에 따라 달라질 수 있으면 다음처럼 써야 합니다.
$$ \frac{\partial C}{\partial t}=\frac{\partial}{\partial x}\!\left(D(C,T)\frac{\partial C}{\partial x}\right) $$
오직 $D$가 공간과 농도에 무관한 상수일 때만$\partial C/\partial t=D\,\partial^2C/\partial x^2$로 단순화됩니다. Shewmon이 강조하듯 실제 확산계수는 조성과 시간에 따라 변할 수 있어, 상수 $D$ 해를 적용하기 전 가정을 확인해야 합니다.
05초기·경계조건이 해를 결정한다
같은 방정식도 처음 농도와 표면 조건이 다르면 해가 달라집니다. 반무한 고체의 초기 농도가$C_0$이고 $t>0$부터 표면 농도를 $C_s$로 일정하게 유지하며 $D$를 상수로 보면, 대표적인 오차함수 해는 다음과 같습니다.
$$ \frac{C(x,t)-C_s}{C_0-C_s}=\operatorname{erf}\!\left(\frac{x}{2\sqrt{Dt}}\right) $$
농도 분포는 어느 지점에서 갑자기 끊기지 않고 매끄럽게 변합니다. 자연스러운 길이 척도는$\sqrt{Dt}$이며, 1차원 무작위 걸음의 RMS 거리 $\sqrt{2Dt}$와는 상수배만 다릅니다.
- 형상은 무한·반무한·박막 중 무엇인가?
- 표면은 일정 농도인가, 일정 플럭스인가?
- $D$를 상수로 볼 수 있는가?
- 내부 반응이나 이동 계면이 없는가?
06온도 의존성 — Arrhenius 관계
원자 점프에는 활성화 자유에너지가 필요합니다. 제한된 온도 범위에서 엔탈피·엔트로피 항을 상수로 묶으면 확산계수는 흔히 다음 Arrhenius 형태를 따릅니다.
$$ D=D_0\exp\!\left(-\frac{Q_d}{RT}\right) $$
공공 기구의 자기확산에서는 유효 활성화 에너지 $Q_d$에 공공 형성과 이동의 에너지가 함께 들어갈 수 있습니다. 침입형 원자는 이미 이용 가능한 빈 침입 자리가 많아 주로 이동 장벽의 영향을 받습니다. $\ln D$ 대 $1/T$가 직선이면 기울기 $-Q_d/R$에서 활성화 에너지를 구할 수 있지만, 상변화나 기구 전환이 있으면 직선이 꺾일 수 있습니다.
07실제 고체: 상호확산, Kirkendall 효과와 빠른 경로
A–B 확산쌍에서 A와 B의 고유 확산계수가 다르면 두 원자의 플럭스가 같지 않습니다. 격자에 고정된 불활성 표지자는 시편 끝에 대해 이동하고, 반대 방향으로 공공의 순플럭스가 생깁니다. 이것이Kirkendall 효과이며, 고체 확산이 원자의 직접 자리교환만으로 일어난다는 관점을 무너뜨리고 공공 기구를 뒷받침한 핵심 증거입니다.
또 완전한 결정 내부의 격자확산만 존재하는 것은 아닙니다. 입계·전위·표면은 원자 배치가 느슨해 활성화 장벽이 낮은 단락 확산(short-circuit diffusion) 경로가 될 수 있습니다. 낮은 온도나 미세립 재료에서는 전체 수송량이 작더라도 이런 경로의 상대적 기여가 커집니다.
08계산 예제 — 침탄의 길이 척도
1000 °C에서 오스테나이트 내 탄소의 유효 확산계수를 설명용 값$D=2\times10^{-11}\,\mathrm{m^2/s}$로 두고 1시간 처리한다고 합시다. 실제 $D$는 탄소 농도와 사용한 데이터 평가식에 따라 달라지므로 공정 설계에서는 Ågren 같은 평가 자료를 사용해야 합니다.
$$ \sqrt{Dt}=\sqrt{(2\times10^{-11})(3600)}\approx2.7\times10^{-4}\,\mathrm{m}=0.27\,\mathrm{mm} $$
0.27 mm는 탄소가 정확히 여기까지 들어갔다는 날카로운 침탄 깊이가 아니라 농도 분포의 특성 길이입니다. 원하는 임계 농도에서의 깊이는 오차함수 해에 $C(x,t)$를 넣어 따로 계산해야 합니다. 시간을 네 배 늘리면 이 길이 척도는 두 배가 되고, 온도를 조금 올리면 $D$가 지수적으로 커져 훨씬 큰 효과가 날 수 있습니다.
09이 글의 용어 사전
본문을 읽다가 “그래서 이 말이 정확히 뭐지?”라는 질문이 생기면 여기에서 다시 확인할 수 있습니다.
- 플럭스flux
- 단위 면적을 단위 시간 동안 통과하는 원자 또는 물질의 양입니다. 확산의 ‘이동 속도’를 면적 기준으로 표현합니다.
- 기울기gradient
- 위치에 따라 어떤 양이 얼마나 빠르게 변하는지를 나타냅니다. 농도 기울기와 화학퍼텐셜 기울기가 확산 방향과 연결됩니다.
- 화학퍼텐셜chemical potential
- 성분 하나가 이동할 때 계의 자유에너지가 얼마나 변하는지를 나타내는 양입니다. 확산의 일반적인 열역학적 구동력입니다.
- 확산계수diffusivity
- 원자나 성분이 얼마나 빠르게 퍼지는지를 나타내는 계수 D입니다. 단위는 m²/s이며 온도·조성·확산 기구에 따라 달라집니다.
- 공공vacancy
- 원래 원자가 있어야 할 결정 격자점이 비어 있는 점결함입니다. 치환형 원자가 이동할 자리를 제공합니다.
- 침입형 원자interstitial atom
- 정규 격자점 사이의 빈 공간에 들어간 작은 원자입니다. 탄소·질소·수소가 대표적입니다.
- 경계조건boundary condition
- 표면 농도나 플럭스처럼 해석 영역의 경계에서 정하는 조건입니다. 같은 확산 방정식도 경계조건에 따라 해가 달라집니다.
- 활성화 에너지activation energy
- 원자가 한 자리에서 다음 자리로 점프하기 위해 넘어야 하는 에너지 장벽입니다.
- Kirkendall 효과Kirkendall effect
- 두 성분의 고유 확산속도가 달라 불활성 표지자와 격자가 이동하는 현상입니다. 공공 확산 기구의 핵심 증거입니다.
- 단락 확산short-circuit diffusion
- 입계·전위·표면처럼 격자 내부보다 빠른 경로를 통한 확산입니다.
10참고자료와 더 읽을거리
- H. Mehrer, Diffusion in Solids — Ch. 2(Fick 법칙), Ch. 10(상호확산과 Kirkendall 효과).
- P. Shewmon, Diffusion in Solids, 2nd ed. — Ch. 1(확산 방정식), Ch. 4(Darken 분석), Ch. 6(빠른 확산경로).
- Balluffi, Allen & Carter, Kinetics of Materials — 원자 이동과 연속체 수송의 열역학적 서술.
- IUPAC, “Definitions and conventions of terms relating to diffusion in the solid state” (1999) — 고체 확산 용어의 권고 정의.
- A. Fick, “Ueber Diffusion” (1855) — Fick 법칙의 원 논문.
- J. Ågren, “A revised expression for the diffusivity of carbon in binary Fe–C austenite” (1986) — 오스테나이트 내 탄소 확산계수 평가.