하드 구슬 시스템의 유리 전이·체킹을 아우르는 파라미터화 방정식

초록

감마분포 기반 포텐셜 에너지 랜드스케이프 이론을 활용해, 무작위 조밀도 η_J 를 입력 변수로 하는 파라미터화된 압축인자 Z(η_J) 식을 제시한다. 이 식은 η_J = 0.62 ~ 0.66 구간의 메타안정·유리 영역을 정확히 묘사하고, 안정 유동 영역까지 자연스럽게 연결한다. 두 특수 경로(η_J = 0.64 ≈ rcp와 이상 유리 전이) 를 상세히 검증하고, 전이점에서 열용량 피크와 전도성 변화가 나타나는 것을 확인한다. 또한, Arrhenius와 과잉 엔트로피 스케일링을 적용한 전송 특성도 η ≈ 0.555 에서 급격히 변함을 보고한다.

상세 분석

본 논문은 기존의 가우시안 기반 포텐셜 에너지 랜드스케이프(PEL) 이론이 유리·체킹 상태에서 압력 발산(pole)을 제공하지 못한다는 한계를 지적한다. 이를 극복하기 위해 감마분포를 채택해 고유구조( inherent structure, IS ) 에너지 분포의 비대칭성을 반영하고, 분포의 형태 파라미터 α 와 스케일 β 를 부피 의존 함수로 설정한다. 감마분포식 f(E) ∝ E^{α‑1} exp(‑E/β) 는 α→∞ 일 때 가우시안으로, α→1 일 때 지수분포로 수렴하므로, 메타안정·유리 영역에서 α 값을 조절해 압축인자 Z 에 필요한 특이점(η_J‑pole)을 자연스럽게 도입한다.

식(12)에서 제시된 전체 자유에너지 A = A_IS + A_vib + A_jam 은 각각 IS 기여, 진동 기여, 그리고 “jammed” 기여로 구분된다. 특히 jammed 기여는 Z_jam = C · (1‑η/η_J)^{‑γ} 형태로, η → η_J 접근 시 압력이 무한히 커지는 특성을 부여한다. 여기서 γ는 감마분포의 α와 연계된 지수이며, η_J는 입력 파라미터이다.

유동 영역에서는 전통적인 Carnahan‑Starling(CS)식이 고밀도에서 비물리적 극점을 갖는 문제를 피하기 위해, 11차까지의 정확한 베릴 계수를 이용한 폐쇄 베릴식(Closed‑virial) Z = 1 + ∑_{n=2}^{11} B_n η^{n‑1} + C · (1‑η/η_cp)^{‑δ} 을 채택한다. η_cp는 결정성 고밀도(≈0.7405)이며, δ는 고밀도에서의 발산 정도를 조절한다.

두 특수 경로는 (1) η_J = 0.64(≈rcp) 경로와 (2) 이상 유리 전이 경로이다. 첫 번째 경로에서는 시뮬레이션 데이터(η = 0.55 ~ 0.66)와 비교해 압축인자, 등온 압축성, 열용량 C_p 등을 계산했으며, 전반적인 오차가 1 % 이하로 매우 우수함을 보였다. 특히 η≈0.555에서 C_p가 급격히 상승하고, 이는 IS 기여와 jammed 기여가 동등하게 작용하기 시작함을 의미한다. 두 번째 경로에서는 η_J를 임의로 0.68 이상으로 설정해 이상 유리 전이(η_g ≈ 0.58)를 재현했으며, 이 경우 Z 는 η → η_g 에서 연속적이지만 기울기가 변하는 “kink”를 보인다. 이는 전통적인 자유에너지 곡면이 아닌, 두 개의 서로 다른 자유에너지 지배 구역이 교차함을 시사한다.

전송 특성 분석에서는 Arrhenius 법칙 D ∝ exp(‑E_a/kT)와 과잉 엔트로피 스케일링 D ∝ exp(‑α S_ex) 을 적용했으며, 두 경우 모두 η≈0.555에서 기울기가 변한다. 이는 구조적 재배열이 급격히 제한되어 동적 활성화 에너지가 증가하고, 과잉 엔트로피가 감소함을 반영한다.

결론적으로, 감마분포 기반 PEL 이론은 유리·체킹 영역의 비대칭 에너지 분포를 정량적으로 포착하고, η_J 를 입력 변수로 하는 파라미터화된 방정식 Z(η_J) 는 전체 밀도 구간(η = 0 ~ 0.66)을 일관되게 기술한다. 이는 기존의 가우시안 기반 모델이 제공하지 못했던 압력 발산과 전이점에서의 열역학·동역학 변화를 동시에 설명할 수 있는 강력한 프레임워크이다.