위상공간 다중히스토그램을 이용한 화학퍼텐셜 계산법

초록

본 논문은 하드 입자 시스템의 자유에너지와 화학퍼텐셜을 구하기 위해 파라미터 의존 포텐셜을 도입한 위상공간 다중히스토그램(PSMH) 방법을 제안한다. 무한히 얇은 정사각형 플레이트릿을 대상으로 등압·등부피 몬테카를로 시뮬레이션을 수행하고, 히스토그램 재가중 기법을 통해 상전이 압력을 정확히 추정하였다.

상세 분석

PSMH 방법은 기존 다중-정준(multicanonical) 기법을 확장한 형태로, 하드 입자 시스템과 이상기체 사이를 매끄럽게 연결하는 파라미터 h에 의존하는 포텐셜 U_h = h n_c를 도입한다. 여기서 n_c는 입자 쌍이 겹치는 횟수이며, h→∞이면 하드 입자 모델(겹침 금지)로, h→0이면 겹침이 허용된 이상기체로 수렴한다. 이 연속적인 변환을 통해 U_h 에 대한 히스토그램 ψ_h(n_c) 를 여러 h 값에서 수집하고, 재가중(reweighting) 관계식 ψ_{h+Δh}(n_c)=C_{h+Δh} C_h^{-1} ψ_h(n_c) exp(-Δh n_c)를 이용해 상태밀도 w(n_c) 를 역산한다. 상태밀도의 최소값 w(0) 은 하드 입자 시스템의 전체 위상공간 부피와 직접 연결되며, 이를 로그 변환해 자유에너지 F 와 화학퍼텐셜 μ 를 얻는다(βF≈-ln w(0)+ln C_0, βμ=βF/N+p* N(L/D)^3).

시뮬레이션 절차는 크게 두 단계로 나뉜다. 첫 번째는 전통적인 등압(NPT) MC를 이용해 다양한 축소압력 p* 에서 시스템을 평형화하고, 압력‑화학퍼텐셜 관계를 초기 추정한다. 두 번째는 등부피(NV_h) MC를 수행하면서 h 를 점진적으로 감소시켜 ψ_h(n_c) 를 수집한다. 히스토그램이 충분히 겹치는 구간을 확보하면 다중히스토그램 기법을 적용해 w(n_c) 를 연속적으로 연결할 수 있다.

정사각형 플레이트릿(N=120) 시스템에 적용한 결과, 등방성‑방사상(이소트로픽‑넴틱) 전이가 p*≈1.13 에서 발생함을 확인했다. 화학퍼텐셜은 βμ≈4.61 이며, 이는 기존 Bates 연구(βμ≈7.3, 보정 후 ≈4.1‑4.6)와 좋은 일치를 보인다. 또한 히스테리시스가 약하게 나타나며, 입자 삽입이 가능한 하드 입자 시스템에서 PSMH가 유효함을 입증했다. 반면 구형 입자와 같이 삽입이 거의 불가능한 경우에는 히스토그램이 겹치지 않아 방법이 실패한다는 한계도 제시하였다.

이와 같이 PSMH는 파라미터 h 를 통해 이상기체와 하드 입자 사이의 연속적인 경로를 만들고, 다중히스토그램 재가중을 이용해 정확한 자유에너지와 화학퍼텐셜을 계산한다는 점에서 기존 방법보다 직관적이며 구현이 비교적 간단하다. 특히, 입자 삽입이 가능한 복합형 하드 입자(플레이트, 원통 등) 시스템의 상전이 연구에 유용한 도구가 될 것으로 기대된다.