헬름홀츠 자유에너지

오류 제보
물리 계산물리 분자동역학 마지막 수정: 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. 어느 쪽을 쓰든 논문 첫 등장에서 정의를 한 번 적어 주는 것이 예의이고, 안 적으면 리뷰어가 반드시 물어본다. ↩

이 문서가 문제를 푸는 데 도움이 됐나요? 한 번의 클릭이 다음 검수 순서를 정합니다.

의견을 반영할게요. 감사합니다.

다음으로 읽기 분자동역학 약 5분 더 탐색 →