[논문 리뷰] Liquid, Glass and Crystal in Two-dimensional Hard disks

Ludger Santen, Werner Krauth|arXiv (Cornell University)|2001. 07. 23.

Material Dynamics and Properties참고 문헌 1인용 수 9

한 줄 요약

이 연구는 고도로 발전된 몬테카를로 방법을 사용하여 이차원 다분산 경질 디스크의 열역학적 및 동역학적 상 행동을 조사한다. 연구는 동역학적 액체-유리 전이가 다분산성과 무관하게 밀도 ρᴳ에서 발생하며, 이는 열역학적 전이를 수반하지 않으며, 다분산성이 증가함에 따라 유도되는 결정체에서 무질서한 상으로의 전이가 연속적(코스터리츠-투스 유형)임을 입증한다. 이는 이전에 1차 전이로 가정된 바와 도전하는 결과이다.

ABSTRACT

We study the thermodynamic and dynamic phase transitions in two-dimensional polydisperse hard disks using Monte Carlo methods. A conventional local Monte Carlo algorithm allows us to observe a dynamic liquid-glass transition at a density, which depends very little on the degree of polydispersity. We furthermore apply Monte Carlo methods which sample the Boltzmann equilibrium distribution at any value of the density and polydispersity, and remain ergodic even far within the glass. We find that the dynamical transition is not accompanied by a thermodynamic transition in this two-dimensional system so that the glass is thermodynamically identical to the liquid. Moreover, we scrutinize the polydispersity-driven transition from the crystal into the disordered phase (liquid or glass). Our results indicate the presence of a continuous (Kosterlitz-Thouless type) transition upon increase of the polydispersity.

연구 동기 및 목표

• 이차원 경질 디스크에서의 동역학적 유리 전이가 열역학적 상전이를 수반하는지 여부에 대한 오랫동안 지속된 논쟁을 해결하기 위해.
• 이차원 경질 디스크에서 다분산성이 증가함에 따라 유도되는 결정체에서 무질서한 상으로의 전이의 성격을 조사하기 위해.
• 유리상 영역 깊이까지도 에르고딕성을 유지하는 비국소 몬테카를로 알고리즘을 개발하고 적용하여 평형 열역학적 성질에 접근하기 위해.
• 다양한 다분산성 수준에서 동역학적 유리 전이의 정확한 위치와 특성을 규명하기 위해.
• 최종 다분산성의 역할이 결정화를 억제하고 평형 시뮬레이션에서 유리상에 도달하는 데 기여하는지를 조사하기 위해.

제안 방법

• 유리 영역에서도 봄올츠만 분포의 에르고딕 샘플링을 가능하게 하는 비국소 몬테카를로 알고리즘을 사용하여 평형 열역학 관측량을 정확히 계산한다.
• 고정 확률 증분 모델을 사용하여 입자 반지름을 균일 누적분포에 따라 할당함으로써 시스템 전반에 일관된 다분산도 ε를 확보한다.
• 국소 몬테카를로 동역학을 사용한 [NV] 상자에서 입자 크기 및 밀도에 따른 효과적 확산계수 Dᵢ를 계산한다.
• 확산계수 Dᵢ(ρ)를 스케일링 형태 Dᵢ(ρ) ∼ (ρᴳᵢ − ρ)^α를 사용하여 0으로 외삽하여 유리 전이 밀도 ρᴳ을 정의한다.
• 결합 정렬 순서 매개변수 ⟨Ψ₆⟩를 측정하고 그 확률 분포를 분석하여 1차 전이 또는 연속 전이의 징후를 탐지한다.
• 밀도 ρ와 다분산도 ε의 범위에서 체계적인 시뮬레이션을 수행하며, 특히 최종 다분산도 한계인 ε = 0.58까지 포함하여 전체 상도표를 매핑한다.

실험 결과

연구 질문

• RQ1 이차원 다분산 경질 디스크에서의 동역학적 액체-유리 전이가 열역학적 상전이를 수반하는가?
• RQ2 유리 전이 밀도 ρᴳ은 다분산도 ε의 정도에 따라 어떻게 달라지는가?
• RQ3 다분산도가 증가함에 따라 결정체에서 무질서한 상(액체 또는 유리)으로의 전이의 성격은 무엇인가?
• RQ4 체적 또는 결합 정렬 순서 매개변수 ⟨Ψ₆⟩의 확률 분포에서 히스테리시스 또는 双극분포가 나타나는가? 이는 1차 전이를 시사한다.
• RQ5 비국소 몬테카를로 방법을 사용하여 유리상 영역 깊이까지 평형 열역학적 성질(예: 상태방정식, 압축성 등)을 신뢰성 있게 계산할 수 있는가?

주요 결과

• 동역학적 액체-유리 전이가 다분산도 ε와 거의 무관한 밀도 ρᴳ에서 발생하며, 이는 단일분산 시스템의 결정화 밀도보다 크게 높다.
• 동역학적 유리 전이와 함께 열역학적 전이가 수반되지 않으며, 유리상은 열역학적으로 액체상과 동일하다. 이는 자유 에너지의 불연속성이 없이 에르고딕성 붕괴가 발생함을 의미한다.
• 다분산도 증가에 의해 유도되는 결정체에서 무질서한 상으로의 전이가 연속적임을 확인하였으며, 부피 또는 ⟨Ψ₆⟩의 확률 분포에서 히스테리시스의 증거가 없었다.
• 평균 결합 정렬 순서 매개변수 ⟨Ψ₆⟩는 다분산도 ε 증가에 따라 부드럽게 감소하며, 그 분포는 여전히 단일극이 유지되어 1차 전이가 아닌 연속 전이임을 지지한다.
• 외삽된 유리 전이 밀도 ρᴳ는 작은 ε에서는 입자 크기에 관계없이 안정적이지만, 매우 큰 ε에서는 유한 체적 효과 및 틈새를 통한 입자 탈출로 인해 약간의 이탈을 보였다.
• 최종 다분산도 한계인 ε = 0.58은 고정 확률 증분 모델에서 도달 가능한 최대값으로 확인되었으며, 이를 초월하면 결정 구조를 유지할 수 없다.