헬름홀츠 자유에너지

편집 역사 토론
물리 계산물리 분자동역학 마지막 수정: 2026-08-15 04:26:09

1. 개요[편집]

헬름홀츠 자유에너지
Helmholtz Free Energy
기호F (IUPAC 권장 표기는 A)
정의F = U − TS
자연변수T, V, N
대응 앙상블정준(NVT) 앙상블
통계역학F = −kBT ln Z
의미등온 과정에서 뽑아낼 수 있는 최대 일

헬름홀츠 자유에너지는 온도와 부피가 고정된 계에서 자발성을 판정하는 열역학 퍼텐셜이자, 등온 과정에서 계로부터 뽑아낼 수 있는 일의 상한을 주는 양 F=U−TSF = U - TS 이다. 내부에너지에서 “온도 TT 에서는 어차피 엔트로피 형태로 묶여 있어 일로 못 바꾸는 몫” TSTS 를 미리 떼어 낸 것이라, 남은 FF 가 자유롭게 쓸 수 있는 에너지라는 뜻에서 자유에너지다.1

깁스 자유에너지 GG 와의 분업은 명확하다. 실험실 비커는 압력이 고정되어 있으므로 GG, 컴퓨터 안의 주기경계 상자는 부피가 고정되어 있으므로 FF. 그래서 시뮬레이션 하는 사람에게는 FF 가 기본값이다. NVT로 돌린 분자동역학이나 메트로폴리스 몬테카를로 방법이 표본하는 앙상블의 특성 함수가 정확히 FF 이고, 상평형·결합상수·용해도 같은 실측 가능한 양은 전부 여기서 파생된다.

범위 정리: 앙상블 자체의 정의와 에너지 요동은 정준 앙상블, 정온정압 쪽 논의는 깁스 자유에너지, ΔF\Delta F 를 실제로 계산하는 알고리즘의 세부는 자유에너지 섭동과 우산 표본추출 문서에 있다. 여기서는 FF 라는 양 자체를 다룬다.

2. 최소 원리와 최대 일[편집]

FF 가 왜 하필 저 조합인지는 제2법칙에서 곧장 나온다. 온도 TT 인 열저장조에 담긴 계가 열 δQ\delta Q 를 받고 일 δW\delta W 를 할 때, 우주 전체의 엔트로피 증가 조건 δQ≤T dS\delta Q \le T\,dS 와 제1법칙 dU=δQ−δWdU = \delta Q - \delta W 를 합치면

δW≤− dU+T dS=− d(U−TS)∣T=− dF\delta W \le -\,dU + T\,dS = -\,d(U - TS)\Big|_T = -\,dF

즉 등온에서 계가 할 수 있는 일은 아무리 잘해야 −ΔF-\Delta F 만큼이고, 등호는 준정적(가역) 과정에서만 성립한다. 여기서 헬름홀츠 본인이 붙인 독일어 이름 freie Energie, 그리고 IUPAC이 권장하는 기호 AA(Arbeit = 일)의 유래가 나온다. 일을 전혀 하지 않는 계라면 δW=0\delta W = 0 이므로

dF≤0(T,V 고정)dF \le 0 \quad (T, V \text{ 고정})

이 되어, 정온정적 조건에서 FF 는 평형에서 최소가 된다. 고립계의 “엔트로피 최대”를 계 변수만으로 다시 쓴 것이 이 최소 원리이며, 주변의 엔트로피 변화를 −ΔU/T-\Delta U/T 라는 계 쪽 양으로 대신 계산해 준다는 점이 회계상의 요점이다.

전미분과 그 부산물은 다음과 같다.

dF=−S dT−p dV+∑iμi dni  ⟹  S=−(∂F∂T)V,N,p=−(∂F∂V)T,N,μi=(∂F∂ni)T,VdF = -S\,dT - p\,dV + \sum_i \mu_i\, dn_i \;\Longrightarrow\; S = -\left(\frac{\partial F}{\partial T}\right)_{V,N},\quad p = -\left(\frac{\partial F}{\partial V}\right)_{T,N},\quad \mu_i = \left(\frac{\partial F}{\partial n_i}\right)_{T,V}

p=−(∂F/∂V)Tp = -(\partial F/\partial V)_T 는 시뮬레이션에서 상태방정식을 뽑는 정공법이고, 교차 2계 미분이 같다는 사실에서 맥스웰 관계식 (∂S/∂V)T=(∂p/∂T)V(\partial S/\partial V)_T = (\partial p/\partial T)_V 가 떨어진다. 내부에너지를 되찾고 싶으면 기브스-헬름홀츠 형태를 쓴다.

U=F+TS=−T2[∂∂T(FT)]VU = F + TS = -T^2\left[\frac{\partial}{\partial T}\left(\frac{F}{T}\right)\right]_V

르장드르 변환으로 이웃 퍼텐셜들과 이어지는 관계도 한 줄이다. U(S,V)→F(T,V)→G(T,p)U(S,V) \to F(T,V) \to G(T,p), 즉 F=U−TSF = U - TS, G=F+pVG = F + pV. 자세한 계보는 르장드르 변환 참고.

퍼텐셜자연변수최소가 되는 조건앙상블
UUS,V,NS, V, N등엔트로피·등적미시정준
FFT,V,NT, V, N등온·등적정준
GGT,p,NT, p, N등온·등압등온등압
Ω\OmegaT,V,μT, V, \mu등온·등적·등화학퍼텐셜대정준

3. 분배함수와의 다리[편집]

통계역학이 열역학에 연결되는 지점이 바로 FF 다. 정준 앙상블의 분배함수 Z=∑ie−βEiZ = \sum_i e^{-\beta E_i} 에 대해

F=−kBTln⁡ZF = -k_B T \ln Z

이 한 줄이 다리 전체다. 여기에 S=−(∂F/∂T)VS = -(\partial F/\partial T)_V 를 적용하면 엔트로피의 깁스 표현 S=−kB∑ipiln⁡piS = -k_B\sum_i p_i \ln p_i 가 나오고, ∂/∂V\partial/\partial V 를 하면 압력이, ∂/∂β\partial/\partial\beta 를 하면 평균 에너지가 나온다. ZZ 하나만 알면 그 계의 열역학 전체를 안다는 말의 실체가 이것이다.

로그가 붙어 있다는 점이 중요하다. 독립된 두 부분계의 분배함수는 곱해지므로(Z=Z1Z2Z = Z_1 Z_2) 자유에너지는 더해진다(F=F1+F2F = F_1 + F_2). 즉 FF 는 크기 성질(extensive)이고, 이 성질 덕분에 FF 를 입자 수로 나눈 값이 열역학적 극한에서 잘 정의된다. 상호작용이 없는 고전 이상기체라면 적분이 다 떨어져

F=−NkBT[ln⁡VNΛ3+1],Λ=h2πmkBTF = -N k_B T\left[\ln\frac{V}{N\Lambda^3} + 1\right], \qquad \Lambda = \frac{h}{\sqrt{2\pi m k_B T}}

가 되고(Λ\Lambda 는 열적 드브로이 파장), 여기서 p=−(∂F/∂V)T=NkBT/Vp = -(\partial F/\partial V)_T = Nk_BT/V 라는 이상기체 법칙이 자동으로 나온다. N!N! 을 빼먹으면 FF 가 크기 성질을 잃고 깁스 역설이 생기는데, 그 유명한 사고가 일어나는 자리도 정확히 이 식이다.

4. 왜 시뮬레이션은 F 를 직접 못 재는가[편집]

여기가 계산 쪽 사람들이 제일 자주 부딪히는 벽이다. FF 는 앙상블 평균이 아니다.

시뮬레이션이 잘하는 일은 ⟨A⟩=∑iAie−βEi/Z\langle A \rangle = \sum_i A_i e^{-\beta E_i}/Z 를 재는 것이다. 이 식에서 ZZ 는 분자와 분모에 함께 있어 약분되므로, 표본을 볼츠만 분포에서 뽑기만 하면 ZZ 를 몰라도 평균이 나온다. 그런데 F=−kBTln⁡ZF = -k_BT\ln Z 에는 ZZ 가 통째로 들어간다. 약분될 상대가 없다. 어떤 물리량 A(x)A(\mathbf{x}) 를 골라 평균해도 FF 가 되지 않는다는 뜻이고, 같은 이유로 엔트로피도 힘이나 압력처럼 궤적에서 곧바로 읽히는 “역학적 관측량”이 아니다.2

원리적으로 ZZ 를 재려면 위상공간 부피 전체를 훑어야 하는데, 3N3N 차원 공간에서 볼츠만 인자가 유의미한 영역은 상상할 수 없이 작은 조각이다. 그래서 현실의 처방은 하나다 — 절대값을 포기하고 차이만 잰다.

  • 열역학 적분(TI): 두 상태를 λ\lambda 로 잇고 ΔF=∫01⟨∂U/∂λ⟩λ dλ\Delta F = \int_0^1 \langle \partial U/\partial\lambda\rangle_\lambda\, d\lambda. 피적분함수는 멀쩡한 앙상블 평균이므로 잴 수 있다.
  • 자유에너지 섭동(FEP): 츠방치히 항등식 ΔF=−kBTln⁡⟨e−βΔU⟩A\Delta F = -k_BT\ln\langle e^{-\beta\Delta U}\rangle_A. 엄밀하지만 지수 평균이라 꼬리 표본이 지배한다. BAR·MBAR이 같은 데이터에서 분산을 줄여 준다.
  • 우산 표본추출: 물리적 반응좌표를 따라가는 F(ξ)F(\xi) 프로파일. 역시 얻는 것은 곡선의 모양이지 절대 높이가 아니다.
  • 자바르진스키 등식: 평형이 아닌 유한 속도 스위칭에서 한 일 WW 에 대해 ⟨e−βW⟩=e−βΔF\langle e^{-\beta W}\rangle = e^{-\beta \Delta F}. 비가역 과정으로도 평형 자유에너지 차이를 뽑을 수 있다는 반직관적 결과이며, 광집게 단분자 실험에서 실제로 쓰인다.
  • 위덤 삽입법: 가상 입자를 무작위 위치에 꽂아 보고 μex=−kBTln⁡⟨e−βΔU⟩\mu^{\text{ex}} = -k_BT\ln\langle e^{-\beta\Delta U}\rangle 로 초과 화학퍼텐셜을 뽑는다. 밀도가 높으면 삽입이 죄다 다른 입자와 겹쳐 유효 표본이 0에 수렴하므로, 조밀한 액체나 고체에서는 사실상 못 쓴다.

기준점이 필요하면 자유에너지를 아는 계에 연결하는 수밖에 없다. 기체는 이상기체로, 결정은 아인슈타인 결정으로 이어 붙여 적분하는 프렌켈-래드 방식이 표준이며, 이렇게 얻은 것이 이른바 절대 자유에너지다. 결국 “절대”라는 말도 어딘가에 해석적으로 풀리는 기준계를 놓았다는 뜻이지 마법이 아니다.

5. 고체의 F — 조화 근사와 준조화 근사[편집]

고체는 사정이 낫다. 원자가 평형 위치 근처에서만 떨기 때문에 퍼텐셜을 2차까지 전개하는 조화 근사가 잘 먹히고, 그러면 FF 가 해석적으로 떨어진다. 힘상수 행렬을 대각화해 얻은 포논 진동수 ωqν\omega_{\mathbf{q}\nu} 에 대해

F=U0+∑qν[ℏωqν2+kBTln⁡(1−e−βℏωqν)]F = U_0 + \sum_{\mathbf{q}\nu}\left[\frac{\hbar\omega_{\mathbf{q}\nu}}{2} + k_B T \ln\left(1 - e^{-\beta\hbar\omega_{\mathbf{q}\nu}}\right)\right]

첫 항이 정적 격자 에너지, 둘째가 영점 진동, 셋째가 열적 기여다. 고온 고전 극한에서는 F→U0+kBT∑ln⁡(βℏω)F \to U_0 + k_BT\sum \ln(\beta\hbar\omega) 로 단순해진다. 표본추출이 전혀 필요 없다는 것이 결정적인 장점 — 진동수만 구하면 온도 의존성이 공짜로 나온다.

문제는 조화 근사 안에서는 진동수가 부피와 무관하다는 점이다. 그러면 열팽창이 0이고, 정적 비열과 정압 비열이 같아지고, 열전도도가 무한대가 된다. 현실과 안 맞는다. 준조화 근사(quasi-harmonic approximation, QHA)는 최소한의 수선으로 이 문제를 푼다. 진동수가 부피에 의존한다고 보고 ωqν(V)\omega_{\mathbf{q}\nu}(V) 를 여러 부피에서 계산해 F(V,T)F(V,T) 곡면을 만든 뒤, 각 온도에서 F+pVF + pV 를 최소화하는 부피를 찾는 것이다. 그렇게 하면 열팽창계수·상태방정식·그뤼나이젠 파라미터 γ=− dln⁡ω/dln⁡V\gamma = -\,d\ln\omega/d\ln V 가 전부 따라 나온다. 지구 내부 광물의 고압 상평형이나 밀도범함수이론 기반 재료 계산에서 표준 도구로 쓰인다.

한계도 분명하다. 융점의 대략 2/3 이상에서는 진짜 비조화 효과(포논-포논 산란, 결함 생성)가 커져 QHA가 무너지고, 부피를 키우면 특정 모드의 진동수가 허수가 되어 계산 자체가 정의되지 않는다. 그 위로는 다시 표본추출 기반 기법(자기무모순 포논, 열역학 적분)으로 돌아가야 한다.

6. 오해하기 쉬운 것들[편집]

  • FF 는 “저장된 에너지”가 아니다. FF 가 줄어드는 과정에서 에너지가 사라지는 게 아니라, UU 는 그대로인데 엔트로피가 늘어 −TS-TS 가 내려가는 경우가 흔하다. 고무줄을 당기면 온도가 오르는 현상(굴-줄 효과)이 순전히 엔트로피 항으로 굴러가는 대표 사례다.
  • FF 의 절대값에는 임의의 상수가 있다. 에너지 원점, 힘장의 결합항 정의, Λ\Lambda 에 들어간 플랑크 상수 처리에 따라 값이 통째로 이동한다. 그래서 논문에 찍힌 FF 단일 값은 거의 항상 의미가 없고, 같은 코드·같은 설정에서 잰 차이만 의미가 있다.
  • F(V)F(V) 는 볼록해야 한다. 열역학적 극한에서 FF 는 VV 에 대해 아래로 볼록하며, 판데르발스 방정식 등에서 보이는 비볼록 고리는 평균장 근사의 인공물이다. 이 고리를 공통접선으로 잘라 내는 조작이 곧 맥스웰 구성이고, 1차 상전이의 공존 영역이 여기서 정의된다.
  • FF 와 AA 는 같은 것이다. 그런데 옛 문헌, 특히 화학 쪽에서 FF 로 깁스 자유에너지를 가리킨 사례가 있어서, 오래된 논문의 FF 는 어느 쪽인지 자연변수를 보고 확인해야 한다.3
  • 유한계에서는 앙상블마다 값이 다르다. FF 와 GG 의 관계는 열역학적 극한에서 르장드르 변환으로 정확하지만, 입자 수백 개짜리 상자에서는 ln⁡N\ln N 규모의 차이가 남는다. 작은 클러스터를 다룰 때 이걸 무시하면 결론이 뒤집히기도 한다.

7. 관련 문서[편집]

8. Footnotes[편집]

  1. 헤르만 폰 헬름홀츠(1821~1894)는 이 개념을 1882년 화학 반응의 “결합력” 논쟁 속에서 정리했다. 당시 통설(톰센-베르텔로 원리)은 “반응의 구동력은 발열량”이라는 것이었는데, 흡열인데도 저절로 일어나는 반응들이 자꾸 나와 곤란한 상황이었다. 답은 엔트로피 항이었고, 그 뒤로 인류는 “발열이면 자발적”이라는 직관을 100년 넘게 못 버리고 있다. ↩

  2. 그래서 “이 시뮬레이션의 엔트로피를 출력해 줘”라는 요청은 “이 시뮬레이션의 압력을 출력해 줘”와 난이도가 다르다. 압력은 비리얼로 매 스텝 계산되지만 엔트로피는 위상공간 부피에 대한 질문이라, 궤적 하나를 아무리 들여다봐도 안 나온다. 이 차이를 모르고 “왜 엔트로피 출력 옵션이 없냐”고 항의하는 것이 MD 입문자의 통과 의례. ↩

  3. 기호 전쟁의 흔적은 지금도 남아 있다. 물리 쪽 교과서는 대체로 FF, 화학 쪽과 IUPAC은 AA. 어느 쪽을 쓰든 논문 첫 등장에서 정의를 한 번 적어 주는 것이 예의이고, 안 적으면 리뷰어가 반드시 물어본다. ↩