[논문 리뷰] Topological Constraints and Rigidity of Network Glasses from Molecular Dynamics Simulations

Mathieu Bauchy|arXiv (Cornell University)|2015. 06. 22.

Complex Network Analysis Techniques인용 수 4

한 줄 요약

이 논문은 분자의 운동역학(MD) 시뮬레이션에 대해 위상적 제약 이론을 엄밀하게 적용할 수 있는 방법을 제안한다. 원자 척도의 구속력과 구부림 제약를 원자 운동의 통계적 분석을 통해 정량화함으로써, 각도 표준편차를 사용하여 손상되지 않은(낮은 각도 편차) 및 손상된(높은 각도 편차) 제약을 구분할 수 있음을 보여주며, 이는 정확한 강성 예측과 실제 유리 형성 체계에 대한 이론의 적용 가능성을 검증한다.

ABSTRACT

Due to its non-crystalline nature, the glassy state has remained one the most exciting scientific challenges. To study such materials, Molecular Dynamics (MD) simulations have been extensively used because they provide a direct view into its microscopic structure. MD is therefore used not only to reproduce real system properties but also benefits from detailed atomic scale analysis. Unfortunately, MD shows inherent limitations because of the limited computational power. For instance, only the simulations of small systems are currently permitted, which prevents from studying small compositional changes, although it is well known that they can dramatically alter system properties. At this stage, it is tempting to follow topological constraint theory, which aims at describing macroscopic properties of the glass relying only on the connectivity of individual atoms, thus considering the complicated glass network as simple mechanical trusses. Thanks to only basic hand calculations, this theory has been successful in predicting complex composition and temperature behavior of glass properties, such as the glass forming ability, the viscosity of the corresponding liquid or the elasticity. The purpose of my PhD work is to connect these successful approaches based only on the topology of the underlying low-temperature network with properties that can be obtained from MD calculations. It should thus allow for an increased applicability of rigidity theory.

연구 동기 및 목표

• 위상적 제약 이론과 분자의 운동역학 시뮬레이션을 연결하기 위해 제약 수를 원자 척도의 동역학에 기반하여 정립한다.
• 특히 비사면체 또는 비팔전자 규칙을 만족하지 않는 원자에 대해 전통적인 제약 수 계산에서 발생하는 모순을 해결하기 위해 MD에서 유도된 통계적 측정치를 사용한다.
• 열적 또는 압력에 의한 응력 하에서 기계적으로 관련 있는(손상되지 않은) 제약와 손상된 제약를 식별한다.
• 온도 및 압력에 따른 손상되지 않은 구부림 제약의 비율을 정량화하여 유리에서의 강성 전이 존재를 검증한다.
• 복잡한 산화물 및 켈레오젠화물 유리에 강성 이론을 적용하기 위한 물리적으로 타당하고 데이터 기반의 방법을 제공한다.

제안 방법

• 다양한 온도와 압력에서 네트워크 유리의 원자 구성을 MD 시뮬레이션을 통해 생성한다.
• 각 결합 각도(BO)의 각도 표준편차(σθ)를 시간에 따라 계산하여 구부림 제약(BB)의 안정성을 평가한다.
• 이중 가우시안 피팅을 각도 분포 함수에 적용하여 손상되지 않은(낮은 σθ) 및 손상된(높은 σθ) 제약를 분리한다.
• 이중 분포의 특징을 바탕으로 손상되지 않은 것과 손상된 것 사이의 임계값(σθ ≈ 17°)을 정의한다.
• 이중 가우시안 모델에 대한 피팅을 통해 온도 함수로 손상되지 않은 BB 제약의 비율 q(T)를 추적한다.
• 첫 번째 衍射 피크의 반폭(full width at half maximum, FWHM)과 같은 구조적 서명과 강성 전이 창문의 존재를 관련시킨다.

실험 결과

연구 질문

• RQ1 분자의 운동역학 시뮬레이션을 사용하여 네트워크 유리에서 기계적으로 손상되지 않은 것과 손상된 구부림 제약를 식별할 수 있는가?
• RQ2 손상되지 않은 구부림 제약의 비율은 온도에 따라 어떻게 변화하며, 이는 강성 전이와 관련이 있는가?
• RQ3 각도 편차(σθ)의 이중분포가 손상되지 않은 것과 손상된 것 사이의 임계값을 정의하는 데 사용될 수 있는가?
• RQ4 최대 상관 길이(L)를 가진 밀도 창이 강성 또는 구조적 질서가 향상된 영역과 일치하는가?
• RQ5 압력과 조성의 조합이 제약 통계에 의해 드러나는 네트워크 유리의 강성에 어떻게 영향을 미치는가?

주요 결과

• 300 K에서 각도 분포는 낮은 σθ에서 뚜렷한 피크를 보이며, 이는 모든 구부림 제약가 손상되지 않았음을 나타낸다.
• 4500 K에서 각도 분포는 크게 넓어지며, 이는 열 활성화로 인해 모든 구부림 제약가 손상되었음을 나타낸다.
• 유리 전이 온도(Tg) 근처에서 이중분포가 나타나며, 하나의 피크는 손상되지 않은 제약를, 다른 하나는 손상된 제약를 나타낸다.
• 손상되지 않은 구부림 제약의 비율 q(T)는 온도가 증가함에 따라 감소하며, 이 데이터는 Mauro-Gupta 공식에 잘 맞는다.
• 이중분포의 최소값을 바탕으로 σθ ≈ 17°가 손상되지 않은 것과 손상된 것 사이의 경계로 식별된다.
• ρ = 2.1에서 3.5 g/cm³ 사이의 밀도 창은 상관 길이(L = 2π / FWHM)의 최댓값을 보이며, 이는 점도의 최소값과 네트워크 형성 원자의 확산 최댓값과 일치하여 구조 전이를 나타낸다.