화학 퍼텐셜

편집 역사 토론
계산물리 계산화학 물리 마지막 수정: 2026-08-04 04:43:30

1. 개요[편집]

화학 퍼텐셜
Chemical Potential
기호μ
정의입자 하나를 더 넣는 데 드는 자유에너지
단위J/입자 또는 J/mol (kJ/mol이 관행)
평형 조건접촉한 두 상의 μ가 같아진다
대표 측정법위도 삽입법, 열역학적 적분, CBMC
다른 얼굴페르미 준위, 전기화학 퍼텐셜

온도가 열의 평형 조건이고 압력이 부피의 평형 조건이라면, 화학 퍼텐셜은 개수의 평형 조건이다.

화학 퍼텐셜(chemical potential) μ\mu 는 다른 시강 변수를 고정한 채 계에 입자 하나를 더 넣을 때 드는 자유에너지 변화로, μ=(∂G/∂N)T,P\mu = (\partial G/\partial N)_{T,P} 로 정의된다. 물이 왜 증발하는지, 소금이 왜 녹는지, 전자가 왜 저쪽 금속으로 넘어가는지가 전부 ”μ\mu 가 높은 쪽에서 낮은 쪽으로”라는 한 문장으로 정리된다.

범위 정리: μ\mu 를 고정한 앙상블과 GCMC 알고리즘은 대정준 앙상블, 두 상태 사이 자유에너지 차이를 재는 일반론은 자유에너지 섭동, 전자의 점유 통계 자체는 페르미-디랙 분포에 있다. 이 문서는 μ\mu 라는 양 자체와 그것을 재는 방법을 다룬다.

2. 앙상블마다 다른 얼굴, 같은 양[편집]

열역학 퍼텐셜 넷을 각각 NN 으로 미분하면 전부 μ\mu 가 나온다. 뭘 고정하고 미분하느냐만 다르다.

μ=(∂U∂N)S,V=(∂F∂N)T,V=(∂G∂N)T,P=−T(∂S∂N)U,V\mu = \left(\frac{\partial U}{\partial N}\right)_{S,V} = \left(\frac{\partial F}{\partial N}\right)_{T,V} = \left(\frac{\partial G}{\partial N}\right)_{T,P} = -T\left(\frac{\partial S}{\partial N}\right)_{U,V}

마지막 표현이 은근히 중요하다. 입자를 넣으면 자유에너지가 오르는 게 아니라 저장조의 엔트로피가 준다는 관점이며, μ\mu 가 음수일 수 있는 이유가 여기 있다(희박한 기체의 μ\mu 는 크게 음수다). 시성 변수의 1차 동차성에서 오일러 관계 U=TS−PV+∑iμiNiU = TS - PV + \sum_i \mu_i N_i 가 나오고, 단성분계에서는 곧장

μ=GN\mu = \frac{G}{N}

즉 입자당 깁스 자유에너지가 된다. 여기서 다시 미분하면 깁스 자유에너지의 시강 변수들을 묶는 깁스-듀엠 관계 S dT−V dP+∑iNi dμi=0S\,dT - V\,dP + \sum_i N_i\,d\mu_i = 0 이 떨어진다. 성분이 cc 개인 계에서 독립적인 시강 변수가 c+1c+1 개뿐이라는, 상률(phase rule)의 뼈대다.

3. 이상기체 + 과잉 분해[편집]

시뮬레이션에서 μ\mu 를 다룰 때는 항상 두 조각으로 쪼갠다.

μ=kBTln⁡ ⁣(ρΛ3)⏟μid+μex⏟상호작용\mu = \underbrace{k_BT \ln\!\left(\rho \Lambda^3\right)}_{\mu^{\text{id}}} + \underbrace{\mu^{\text{ex}}}_{\text{상호작용}}

앞은 상호작용이 없는 이상기체 몫으로 밀도 ρ=N/V\rho = N/V 와 열적 드브로이 파장 Λ\Lambda 만 알면 해석적으로 끝난다. 실제로 계산이 필요한 건 뒤의 과잉 화학 퍼텐셜 μex\mu^{\text{ex}} 뿐이다. 이 분해가 유용한 이유는 명확하다 — Λ\Lambda 에 들어 있는 플랑크 상수는 고전 시뮬레이션에서 아무 물리도 하지 않고 두 상태를 비교할 때 상쇄되므로, 상평형·용해도 같은 문제는 전부 μex\mu^{\text{ex}} 의 경쟁으로 환원된다. 실험 화학의 활동도 계수나 퓨가시티 계수도 결국 μex\mu^{\text{ex}} 를 다르게 포장한 것이다.

4. 위도 삽입법 — 유령 입자를 던져 본다[편집]

μex\mu^{\text{ex}} 를 직접 재는 가장 우아한 방법이 위도 삽입법(Widom insertion, 1963)이다. 정준 앙상블 ZN+1/ZNZ_{N+1}/Z_N 을 정직하게 쓰면 다음이 나온다.

μex=−kBTln⁡⟨e−βΔU+⟩N\mu^{\text{ex}} = -k_BT \ln \left\langle e^{-\beta \Delta U^{+}} \right\rangle_{N}

ΔU+\Delta U^{+} 는 상자 안 무작위 위치에 입자 하나를 시험 삽입했을 때 늘어나는 퍼텐셜 에너지다. 핵심은 그 입자를 실제로 넣지 않는다는 것 — 에너지만 계산하고 버린다. 그래서 “유령 입자”라 부르며, 평형 궤적을 돌리면서 매 프레임 수백 번 던져도 궤적 자체는 전혀 교란되지 않는다. 딱딱한 구 계에서는 이 식이 특히 예쁘다. e−βΔUe^{-\beta\Delta U} 가 0 아니면 1이므로 μex=−kBTln⁡P0\mu^{\text{ex}} = -k_BT\ln P_0, 즉 겹치지 않고 들어갈 확률의 로그가 된다.

문제는 밀집계에서 그 P0P_0 가 사라진다는 것이다. 액체 밀도에서 무작위로 던진 입자는 거의 항상 기존 입자와 겹치고, 레너드-존스 퍼텐셜의 반발 코어에서 ΔU\Delta U 가 수십 kBTk_BT 로 튀면 e−βΔUe^{-\beta\Delta U} 는 10−1010^{-10} 수준이 된다. 그런데 평균값을 지배하는 것은 드물게 열리는 빈 공동이다. 표본의 99.99%가 0을 기여하고 아주 가끔 등장하는 공동이 전부를 결정하는, 자유에너지 섭동에서 익숙한 지수 평균의 저주 그 자체다. 게다가 이 편향은 조용하다 — 통계 오차만 보면 수렴한 것처럼 보이는데 값 자체가 체계적으로 치우친다.1 사슬 분자라면 사슬 전체를 한 번에 꽂아야 하니 아예 시도조차 무의미하다.

5. 그래서 실제로는 이렇게 잰다[편집]

  • 열역학적 적분. 결합 파라미터 λ\lambda 로 입자를 0에서 1까지 천천히 켜면서 μex=∫01⟨∂U/∂λ⟩λ dλ\mu^{\text{ex}} = \int_0^1 \langle \partial U/\partial\lambda\rangle_\lambda \, d\lambda 를 적분한다. 한 방에 삽입하는 대신 여러 단계로 나누는 것이라 밀집계에서도 작동한다. 단 λ→0\lambda \to 0 근처에서 LJ 코어가 특이점을 만들기 때문에 소프트코어 퍼텐셜로 코어를 뭉개는 처방이 필수다.2
  • 형태편향 몬테카를로(CBMC). 사슬 분자를 통째로 던지지 않고 원자 하나씩 성장시키되, 매 단계에서 에너지가 낮은 방향으로 편향 표본추출하고 그 편향을 로젠블루스 가중치로 정확히 되돌린다. 고분자·알케인의 μ\mu 계산을 가능하게 만든 결정적 도구.
  • 팽창앙상블. λ\lambda 를 이산 사다리로 만들어 앙상블 자체를 확장하고, 계가 사다리를 랜덤워크하도록 가중치를 조정한다. 병렬 템퍼링과 사촌지간이며, 삽입-삭제를 한 번에 하지 않고 왕복시킨다는 발상이 핵심이다.
  • 오버랩 표본추출·BAR. 정방향과 역방향 작업 분포를 함께 써서 편향을 줄이는 계열. 삽입만 하는 위도법이 한쪽 꼬리에 의존한다는 약점을 정면으로 공략한다.

공통 교훈은 하나다. μ\mu 는 평균값이 아니라 자유에너지라서 궤적에서 그냥 읽히지 않는다. 온도나 압력처럼 순간 추정량이 있는 양이 아니고, 항상 별도의 계산 설계를 요구한다.

6. 페르미 준위는 사실 화학 퍼텐셜이다[편집]

물리·화학·전자공학이 서로 다른 이름으로 부르던 양이 같은 것이었다는 사실은 볼 때마다 유쾌하다.

전하를 띤 입자에는 정전 퍼텐셜 ϕ\phi 가 얹히므로 실제 평형을 결정하는 것은 전기화학 퍼텐셜 μ~=μ+zeϕ\tilde{\mu} = \mu + ze\phi 다. 두 상이 접하면 μ~\tilde\mu 가 같아질 때까지 전하가 이동하고, 그 과정에서 생긴 ϕ\phi 차이가 접촉 전위로 측정된다. 반도체에서 전자의 μ~\tilde\mu 를 부르는 이름이 바로 페르미 준위다. 페르미-디랙 분포의 분모에 들어 있는 μ\mu 가 그것이고, 서로 다른 도핑의 두 반도체를 붙였을 때 페르미 준위가 일직선으로 정렬되며 밴드가 휘는 그림이 곧 ”μ\mu 를 맞추는 과정”이다. 비평형에서 전자와 정공에 준페르미 준위를 따로 주는 반도체 소자 시뮬레이션의 관행도 같은 언어다.

배터리도 같은 이야기다. 리튬 이온 전지의 개방회로 전압은 양극과 음극에서 리튬의 화학 퍼텐셜 차이를 전자 전하로 나눈 값이다.

VOCV=−μLi양극−μLi음극eV_{\text{OCV}} = -\frac{\mu_{\text{Li}}^{\text{양극}} - \mu_{\text{Li}}^{\text{음극}}}{e}

밀도범함수이론으로 삽입 화합물의 전 에너지를 계산해 전압을 예측하는 계산재료 분야가 정확히 이 식 하나에 서 있다.3 통계역학의 시험 삽입, 고체물리의 페르미 준위, 전기화학의 셀 전압이 전부 “입자 하나를 더 넣는 값”이라는 같은 뿌리에서 자란 셈이다.4

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 위도 삽입법의 오차 막대는 “이 값이 얼마나 흔들리는가”만 말해 준다. 표본이 진짜 봉우리를 한 번도 못 밟았다면 흔들림도 작다. 오차 막대가 예쁘게 작은 위도법 결과가 가장 위험한 결과인 이유다. ↩

  2. 소프트코어를 안 쓰고 λ\lambda 를 그냥 선형으로 켜면 적분값이 λ→0\lambda\to 0 에서 발산하는 것을 수치적으로 목격하게 된다. 이걸 ”λ\lambda 창을 더 촘촘히 나누면 되겠지”로 대응하다가 며칠을 태우는 것이 이 바닥의 통과의례. ↩

  3. 예측 정확도는 대개 0.1~0.3 V 수준으로, 후보 물질을 걸러내기에는 충분하고 셀 설계 스펙으로 쓰기에는 모자라다. 그래도 물질 하나 합성하는 데 드는 시간과 비교하면 압도적으로 싸다. ↩

  4. 그래서 화학과·물리과·전자과 학생이 같은 테이블에 앉으면 30분쯤 지나서야 “아 그거 우리가 페르미 준위라고 부르는 거네”가 나온다. 학제 간 협업의 절반은 용어 사전을 맞추는 일이다. ↩