레너드‑존스 초냉각 액체와 이상 유리의 구조 변화를 밝히는 새로운 적분 방정식

초록

본 연구는 복제된 HMSA 이론과 최적 분할(ODS) 스킴을 결합해 레너드‑존스(LJ) 유체를 매우 낮은 온도 영역까지 확장한다. 복제된 시스템의 상호작용을 통해 겹침(order parameter Q)을 계산하고, 초냉각 액체와 이상 유리 단계에서 방사형 분포함수(g(r))에 새로운 피크가 나타나는 온도 T₀를 규명한다. 이는 장거리 질서는 없지만 근거리에서 fcc‑유사 구조가 형성됨을 시사한다.

상세 분석

본 논문은 전통적인 액체 상태 이론을 초냉각 및 이상 유리 영역에 적용하기 위해 두 가지 핵심 기술을 도입한다. 첫째, 복제된 열역학적 자기일관성 하이브리드 평균구면 근사(HMSA)를 사용해 두 개의 약하게 결합된 복제(클론) 시스템을 대칭적인 이원 혼합물로 모델링한다. 복제 사이의 약한 짧은 거리 인력 u₁₂(r)은 ε₁₂→0 한계에서 사라지도록 설계되어, 복제 간 겹침 Q와 g₁₂(0)을 직접 계산할 수 있게 한다. 둘째, 레너드‑존스 포텐셜을 최적 분할(Optimized Division Scheme, ODS)으로 나누어 단거리(짧은 거리)와 장거리(긴 거리) 부분을 정의한다. ODS는 기존의 WCA 분할이 핵심 영역에서 브리지 함수와 캐비티 함수를 부정확하게 재현하는 문제를 해결하고, 특히 r≈0에서의 브리지 함수 값을 정확히 맞추도록 p=2와 다항식 계수를 조정한다. 이러한 분할은 HMSA 폐쇄식의 α 파라미터를 열역학적 자기일관성(압축성 및 비리얼 방정식) 조건에 맞게 최적화할 수 있게 하며, 결과적으로 g₁₁(r)과 g₁₂(r)의 정밀한 계산을 가능하게 한다.

복제 방법론에서는 겹침 Q=8πρ∫₀^∞ g₁₂(r) w(r) r² dr 를 정의하고, w(r)은 복제 사이 인력 형태에 따라 결정된다. 온도와 밀도가 주어지면 R‑HMSA 방정식을 수치적으로 풀어 g₁₁(r)과 g₁₂(r)를 얻는다. 고온에서는 ε₁₂를 끄면 복제가 완전히 해리되어 Q는 무작위값 Q_random, g₁₂(0)=1을 보이며, 이는 일반적인 단일 성분 HMSA 결과와 일치한다. 온도를 낮추면 Q와 g₁₂(0)가 급격히 증가하고, 특정 밀도‑온도 조합에서 복제 간 겹침이 비정상적으로 크게 나타나며, 이는 동적 전이 온도 T_D와 자유에너지 계산을 통해 얻은 임계 온도 T_cr 사이에 존재한다.

가장 중요한 발견은 T₀<T_cr 이하에서 g₁₁(r)의 첫 번째 피크와 두 번째 피크 사이에 새로운 중간 피크가 나타난다는 점이다. 이 피크는 기존의 LJ 액체에서 관찰되지 않았으며, 구조적으로는 fcc 결정의 근거리 배치와 유사한 짧은 거리 질서를 암시한다. 그러나 장거리 질서는 여전히 무질서한 상태를 유지하므로, 이상 유리의 특성을 만족한다. 이러한 구조적 전이는 밀도 의존적인 온도 T₀가 존재함을 보여주며, 이는 복제 이론이 예측한 “전구 전이(precursor transition)”와도 일맥상통한다.

결과적으로, R‑HMSA + ODS 조합은 기존 Monte‑Carlo 혹은 분자동역학 시뮬레이션이 접근하기 어려운 초저온 영역까지 정확한 구조 정보를 제공한다. 이는 복제 기반 정적 이론이 동적 전이와 구조적 변화를 동시에 포착할 수 있음을 입증하며, 이상 유리의 존재와 그 미세구조에 대한 새로운 물리적 통찰을 제공한다.