초냉각 스틸링거‑웨버 규소의 엔트로피 비단조 온도 의존성

초록

스틸링거‑웨버(SW) 포텐셜을 이용한 규소의 초냉각 영역에서, 고밀도 액체(HDL)와 결정상 사이의 과잉 깁스 자유에너지(G⁽ᵉ⁾)를 정밀하게 계산하였다. G⁽ᵉ⁾의 온도 기울기를 통해 HDL의 절대 엔트로피가 1065 K 이상에서는 평형값을 유지하지만, 1060 K 근처 비평형 상태로 전이하면서 급격히 증가함을 확인했다. 또한 부피 분포가 넓어지면서 압력‑부피 곡선에 반데르발스(VDW) 루프가 나타나, HDL‑LDL 전이의 열역학적 원인을 제시한다.

상세 분석

본 연구는 최근 개발된 열역학적 적분(thermodynamic integration) 방법을 NPT‑Monte Carlo 시뮬레이션에 적용하여, SW 규소의 고밀도 액체(HDL)와 결정상 사이의 과잉 깁스 자유에너지 G⁽ᵉ⁾를 온도 1060 K–1070 K 구간에서 정밀히 측정하였다. G⁽ᵉ⁾(T)의 기울기 dG⁽ᵉ⁾/dT = –S⁽ᵉ⁾ 를 이용해 절대 엔트로피 변화를 직접 계산했으며, 1065 K 이상에서는 HDL이 열역학적 평형을 이루어 엔트로피가 감소하는 일반적인 거동을 보인다. 그러나 1060 K에서 시뮬레이션이 비평형 상태에 머무르는 경우, 엔탈피가 평형값보다 낮아지면서 엔트로피가 비정상적으로 상승한다. 이는 자유에너지 장벽이 얕아 HDL이 LDL(저밀도 액체) 혹은 비정질 상으로 전이하기 전, 충분히 샘플링되지 못한 결과로 해석된다.

부피 분포 분석에서는 온도가 낮아질수록 확률밀도 함수가 비대칭적으로 넓어지고, 압력‑부피(p‑v) 곡선에 양의 기울기를 가진 구간이 나타난다. 특히 1060 K와 1065 K에서는 두 개의 접선선을 그어 만든 VDW 루프가 관찰되었으며, 이는 동일 온도·압력에서 화학 퍼텐셜이 같은 두 상이 공존할 수 있음을 의미한다. 루프 양쪽 끝의 접점에서 ΔU와 pΔV가 동일 부호를 갖는 점은 엔탈피 차이가 존재함을 보여, 전이가 단순한 잠열 방출이 아니라 엔트로피와 부피 변동을 동반한 복합적인 과정임을 시사한다.

또한 저자들은 자유에너지 곡선 F(v) 주변을 스핀odal(불안정점) 기준으로 삼아 3차 테일러 전개 F(v)=F_s –p_s(v–v_s)+ (1/3!)F’’’_s (v–v_s)³ 로 근사하였다. 이 전개는 시뮬레이션 데이터와 매우 잘 맞으며, 스핀odal 양쪽이 각각 다른 상의 방정식을 갖는다는 물리적 해석을 가능하게 한다.

열역학적 적분 경로는 세 단계(잠재적 상호작용 감소 → 외부 가우시안 포텐셜 적용 → 원래 포텐셜 복원)로 구성되며, 각 단계에서 Bennett Acceptance Ratio(BAR) 방법을 이용해 자유에너지 차이를 정확히 계산하였다. 순방향·역방향 적분 결과가 일치함을 통해 히스테리시스가 없음을 확인했고, 이는 경로가 완전 가역적임을 의미한다.

결과적으로, HDL‑LDL 전이는 1060 K 근처에서 자유에너지 장벽이 매우 얕아져 비평형 상태가 장기적으로 유지될 수 있음을 보여준다. 엔트로피가 비단조적으로 증가하는 현상은 부피 변동과 VDW 루프 형성, 그리고 스핀odal 근처의 비선형 자유에너지 구조와 직접 연관된다. 이러한 열역학적 통찰은 실험적 초냉각 규소에서 관측되는 비정질 전이 메커니즘을 이해하는 데 중요한 기반을 제공한다.