[논문 리뷰] Random close packing of polydisperse hard spheres

Michiel Hermes, Marjolein Dijkstra|arXiv (Cornell University)|2009. 03. 24.

Material Dynamics and Properties인용 수 44

한 줄 요약

이 연구는 이벤트 기반 분자 동역학 시뮬레이션을 사용하여 다분산 경량 구체의 봉쇄 밀도를 조사하며, 압축 속도가 감소함에 따라 평형화가 향상되어 봉쇄 밀도가 증가하는 것으로 밝혀졌다. 결과적으로 10% 다분산 구체의 경우 봉쇄 패킹은 0.638에서 0.658 사이이며, 단일분산 구체의 경우 0.635에서 0.645 사이로 나타나, 봉쇄 상태가 압축 중 정지된 유리 상태의 무한 압력 극한임을 지지한다.

ABSTRACT

We study jammed configurations of hard spheres as a function of compression speed using an event-driven molecular dynamics algorithm. We find that during the compression, the pressure follows closely the metastable liquid branch until the system gets arrested into a glass state as the relaxation time exceeds the compression speed. Further compression yields a jammed configuration that can be regarded as the infinite pressure configuration of that glass state. Consequently, we find that the density of jammed packings varies from 0.638 to 0.658 for polydisperse hard spheres and from 0.635 to 0.645 for pure hard spheres upon decreasing the compression rate. This demonstrates that the density at which the systems falls out of equilibrium determines the density at which the system jams at infinite pressure. In addition, we give accurate data for the jamming density as a function of compression rate and size polydispersity.

연구 동기 및 목표

• 경량 구체 시스템에서 봉쇄 밀도가 압축 속도와 크기 다분산성에 어떻게 의존하는지 조사하기.
• 무작위 밀접패킹이 고정된 상태인지, 운동적 경로에 따라 변하는지 명확히 하기.
• 실험적 콜로이드 시스템과 비교하기 위해 정확한, 압축 속도에 따라 달라지는 봉쇄 밀도 데이터를 제공하기.
• 봉쇄 구성이 유리 상태의 무한 압력 극한과 대응한다는 가설을 시험하기.
• 압축 속도와 다분산성 영향을 고려하여 보고된 무작위 밀접패킹 밀도의 모순을 해결하기.

제안 방법

• 입자 크기 증가 동안 온도와 다분산성을 일정하게 유지하기 위해 수정된 알고리즘을 사용한 이벤트 기반 분자 동역학 시뮬레이션.
• 압축 속도는 시간에 따라 변화하는 직경 증가율 Γ = dσ/dt를 통해 제어하며, MD 시간 단위를 사용.
• 음수 직경을 방지하고 물리적 현실성을 확보하기 위해 다분산성을 로그정규분포에서 샘플링.
• 봉쇄 근처에서 Eq. (1)을 사용해 상태방정식 데이터를 피팅하여 무한 압력으로의 외삽을 수행하고 φJ를 결정.
• 통계적 정확도를 확보하기 위해 각 조건에 대해 50회의 독립 시뮬레이션을 실시.
• 유한한 크기 효과는 2000에서 200,000개의 입자를 대상으로 비교하여 테스트하였으며, 통계 오차 범위 내에서 양호한 일치를 보였다.

실험 결과

연구 질문

• RQ1 경량 구체의 봉쇄 밀도 φJ는 압축 속도 Γ에 어떻게 의존하는가?
• RQ2 크기 다분산성이 봉쇄 밀도에 미치는 영향은 어느 정도이며, 결정화를 방지하는 데 기여하는가?
• RQ3 봉쇄 구성은 유리 상태의 무한 압력 극한으로 해석될 수 있는가?
• RQ4 왜 실험 및 시뮬레이션에서의 무작위 밀접패킹 밀도 값이 크게 다를까?
• RQ5 봉쇄 밀도는 무한히 느린 압축으로 어떻게 외삽되며, 그 한계 값은 얼마인가?

주요 결과

• 단일분산 경량 구체의 봉쇄 밀도 φJ는 Γ⁻¹이 10에서 1000으로 증가함에 따라 0.635에서 0.645 사이로 변동한다.
• 10% 다분산 경량 구체의 경우, 압축 속도가 감소함에 따라 φJ는 0.638에서 0.658로 증가하여 느린 속도에서 더 강한 평형화가 일어남을 나타낸다.
• 압축 중 압력은 체계가 비평형 상태로 빠지기 전까지는 메타안정성 유체 상태방정식을 따르며, 이후 유리 상태로 정지된다.
• 봉쇄 근처에서 압력의 역수 1/(βPσ³)는 부피 분율에 대해 선형적 의존성을 보이며, 봉쇄와 유리 상태의 무한 압력 극한 사이의 이론적 연관성을 지지한다.
• 느린 압축 속도(Γ⁻¹ > 1000)에서 단일분산 체계에서는 결정화가 발생하지만, 크기 다분산성이 이를 억제하여 더 넓은 압축 속도 범위에 접근할 수 있도록 한다.
• Schaertl 등(1994)의 실험적 값은 단일 실행에서 통계 정확도가 낮아 우리 결과와 크게 다름을 보이며, 보정에 사용하는 데 의심을 제기한다.