[논문 리뷰] Crystallization of the Wahnström Binary Lennard-Jones Liquid

Ulf R. Pedersen, Nicholas P. Bailey외 2명|ArXiv.org|2007. 06. 06.

Liquid Crystal Research Advancements참고 문헌 1인용 수 3

한 줄 요약

이 연구는 마이크로초 시간스케일에서 분자역학 시뮬레이션을 통해 와른스트롬 이성분 레너드-존스 액체의 자발적 결정화를 보고하며, 안정한 결정 구조로 MgZn₂ (C14) 레이브스 상이 확인되었다. 결정화 메커니즘은 nucleation과 성장에 의해 지배되며, 최대 결정화 속도는 T = 0.60 ± 0.02에서 발생하고, 밀도 ρ = 0.82일 때 융해 온도는 Tₘ = 0.93 ± 0.02로 추정된다.

ABSTRACT

We report observation of crystallization of the glass-forming binary Lennard-Jones liquid first used by Wahnström [G. Wahnström, Phys. Rev. A 44, 3752 (1991)]. Molecular dynamics simulations of the metastable liquid on a timescale of microseconds were performed. The liquid crystallized spontaneously. The crystal structure was identified as MgZn_2. Formation of transient crystallites is observed in the liquid. The crystallization is investigate at different temperatures and compositions. At high temperature the rate of crystallite formation is the limiting factor, while at low temperature the limiting factor is growth rate. The melting temperature of the crystal is estimated to be T_m=0.93 at rho=0.82. The maximum crystallization rate of the A_2B composition is T=0.60+/-0.02.

연구 동기 및 목표

• 장기간 시간스케일에서 불안정한 유리형성 와른스트롬 이성분 레너드-존스 액체의 결정화 거동을 조사하기 위해.
• 결정화 과정 중 형성된 결정 구조를 규명하고, 구조 분석을 통해 그 정체성을 검증하기 위해.
• 다양한 온도에서 결정화를 지배하는 동역학적 제약 조건—핵형성과 성장 중 어느 것이 한계를 이루는지 규명하기 위해.
• 결정 상의 융해 온도 Tₘ과 최대 결정화 속도를 갖는 온도를 추정하기 위해.
• 액체, 초냉각 액체, 유리상 및 결정상 상태를 모두 포함하는 시스템에서 구조적 회복, 일시적 결정립 형성 및 결정 성장 간의 상호작용을 탐색하기 위해.

제안 방법

• 시간 단위 t* = 0.01을 사용한 Gromacs 소프트웨어 패키지를 이용한 비평형 분자역학 시뮬레이션을 수행하였다.
• 레너드-존스 포텐셜을 사용하였으며, 파라미터는 σ_BB = 1.2σ_AA, m_B = 2m_A, 그리고 모든 조합에 대해 ε_ij = 1로 설정하였다.
• 일정한 온도를 유지하기 위해 Nosé-Hoover 열역학장치를 사용하였고, 주기적 경계 조건을 적용하였다.
• 2.0σ에서 3.5σ로 이격된 스위치 함수를 적용하여 3.5σ에서 포텐셜을 부드럽게 잘라내었다.
• N_A = N_B = 512 및 N_A = 1024, N_B = 512로 설정된 밀도 ρ = 0.75와 ρ = 0.82σ⁻³_AA에서 시뮬레이션을 수행하였다.
• VMD를 사용해 시각화하고, 공통 이웃 분석을 통해 결정립을 식별하며 결정 구조를 확인하였다.

실험 결과

연구 질문

• RQ1 와른스트롬 이성분 레너드-존스 액체는 마이크로초 시간스케일에서 자발적으로 결정화되는가?
• RQ2 결정화 과정 중 형성된 결정 구조는 무엇이며, A₂B 조성과 어떻게 관련되어 있는가?
• RQ3 다양한 온도에서 결정화를 지배하는 동역학적 제약 조건—핵형성 또는 성장 중 어느 것이 한계를 이루는가?
• RQ4 주어진 밀도와 조성 조건에서 MgZn₂ 상의 융해 온도 Tₘ는 얼마인가?
• RQ5 A₂B 조성에서 결정화 속도가 최대가 되는 온도는 얼마인가?

주요 결과

• 와른스트롬 이성분 레너드-존스 액체는 마이크로초 시간스케일에서 시뮬레이션을 통해 자발적으로 결정화되며, 장기적으로는 불안정한 액체 상태가 진정으로 안정하지 않음을 보여준다.
• 형성된 결정 구조는 A₂B 스토이히오메트리와 입자 크기 비율에 부합하는 P6₃/mmc 대칭을 가지는 MgZn₂ (C14) 레이브스 상으로 확인되었다.
• 초냉각 액체 상태에서 일시적인 결정립이 형성되며, 공통 이웃 분석을 통해 확인되었고, 전체 결정화 이전에 일시적인 에너지 감소로 관찰되었다.
• 고온에서는 (T = 0.624) 핵형성이 속도 제한 단계이며, 저온에서는 (T = 0.583) 결정 성장이 제한 요소이다.
• MgZn₂ 결정의 융해 온도는 ρ = 0.82에서 Tₘ = 0.93 ± 0.02로 추정되었고, A₂B 조성 조건을 만족한다.
• 최대 결정화 속도는 T = 0.60 ± 0.02에서 발생하여, 결정화 동역학에 최적의 온도 창이 존재함을 시사한다.