Percolation Criticality of Amorphous-Amorphous Transitions in Compressed Glasses

본 연구는 대규모 분자 역학 시뮬레이션과 퍼콜레이션 이론을 사용하여 압축된 실리카 유리의 저밀도-고밀도 전이가 구조적 클러스터의 임계 퍼콜레이션에 의해 구동됨을 밝히며, 이는 강성 퍼콜레이션 메커니즘을 시사하는 임계 지수를 나타내고 결합형 및 비결합형 무정형 물질 사이의 공통적인 변형 원리를 강조한다.

핵심

유리잔의 물이나 창문 유리를 상상해 보세요. 여러분은 그것들을 단단하고 경직된 물체로 알고 있습니다. 하지만 여러분이 원자 크기로 작아져서 그 내부를 들여다본다면, 아주 작고 뒤엉킨 미세한 블록들의 혼란스러운 그물망을 보게 될 것입니다. 실리카 유리(창문의 재료)에서 이 블록들은 규소와 산소로 이루어진 피라미드 형태(사면체)를 띠고 있습니다.

이 논문은 마치 고성능 카메라로 확대하여 관찰하는 영화와 같습니다. 엄청난 압력—대기압의 무려 35만 배에 달하는 압력—으로 유리를 쥐어짤 때 어떤 일이 일어나는지를 보여줍니다. 과학자들은 유리가 녹거나 깨지지 않고 어떻게 모양을 바꾸는지, 즉 "비정질-비정질 전이(amorphous-amorphous transition)"라고 불리는 과정을 이해하고자 했습니다.

다음은 그들의 발견을 쉬운 비유를 통해 들려주는 이야기입니다:

1. 콘서트장의 인파 (구조)

유리를 콘서트장에 모인 거대한 인파라고 생각해 보세요.

• 정상 압력일 때: 모든 사람이 느슨하고 열린 대형으로 서 있습니다. 실리카 유리에서 "사람들"은 규소 원자이며, 이들은 네 명의 이웃과 손을 잡고 완벽한 피라미드 모양(사면체)을 형성합니다. 이들 사이에는 빈 공간이 많아서 구조가 "흐느적거리고" 압축하기 쉽습니다.
• 압력이 증가함에 따라: 콘서트 홀이 점점 줄어든다고 상상해 보세요. 인파가 압착됩니다. 사람들은 더 이상 완벽한 피라미드 모양을 유지할 수 없습니다. 사람들은 서로 부딪히기 시작하며, 얼마나 많은 이웃과 손을 잡을지를 바꿉니다. 어떤 이들은 5명, 그다음에는 6명과 손을 잡기 시작합니다.

2. "퍼콜레이션(Percolation)" 게임 (거대한 변화)

과학자들은 퍼콜레이션이라는 개념을 사용했습니다. 스펀지에 물을 붓는 장면을 상상해 보세요.

• 낮은 압력: 스펀지에는 구멍이 있지만, 이 구멍들은 모두 떨어져 있습니다. 물을 부으면 물은 작은 주머니 안에 갇히게 됩니다. 물이 끝까지 흐르지 못합니다. 유리 내부에서도 이 "피라미드" 모양들은 고립된 섬과 같습니다.
• 임계 압력: 더 세게 쥐어짤수록 구멍들이 서로 연결되기 시작합니다. 갑자기, 스펀지의 위에서 아래까지 이어지는 거대하고 연속적인 경로가 형성됩니다. 이제 물이 흐릅니다!
• 유리 내부에서: 과학자들은 특정 압력 지점에서 새로운 모양들(5면체 또는 6면체 블록 등)이 갑자기 서로 연결되어 유리 전체를 가로지르는 거대하고 연속적인 사슬을 형성한다는 것을 발견했습니다. 이것이 바로 "퍼콜레이션 전이"입니다. 이것은 유리가 근본적으로 더 밀도가 높은 상태로 재편되는 순간입니다.

3. 인파를 바라보는 두 가지 관점

연구진은 마치 서로 다른 카메라 각도에서 영화를 보는 것처럼, 두 가지 방식으로 이 인파를 관찰했습니다.

• "결합된" 관점 (악수): 그들은 누가 직접 손을 잡고 있는지(화학 결в)를 보았습니다. 그들은 피라미드 모양들이 어떻게 악수 방식을 바꾸는지 관찰했습니다.
• "비결합된" 관점 (개인 공간): 그들은 악수를 무시하고, 서로 손을 잡고 있든 아니든 상관없이 누가 서로 가까이 있는지만을 보았습니다. 이것은 사람들이 손을 잡고 있지는 않지만, 그냥 서로 근처에 서 있는 군중을 보는 것과 같습니다.

놀라운 점: 두 카메라 모두 정확히 같은 이야기를 보여주었습니다! "악수" 관점과 "개인 공간" 관점 모두에서, 유리가 동일한 순서로 변형된다는 것이 나타났습니다: 먼저 느슨한 모양들이 연결되고, 그다음 밀도가 높은 모양들이 주도권을 잡습니다. 이는 유리가 어떻게 변하는지를 결정하는 규칙이, 원자들이 "손을 잡고 있든"(실리카처럼) 혹은 그냥 "서로 부딪히고 있든"(얼어붙은 물/얼음처럼) 관계없이 보편적임을 시사합니다.

4. "마법의 숫자"와 게임의 규칙

과학자들은 이 변형이 표준적인 규칙(확률 게임과 같은)을 따르는지, 아니면 자신만의 특별한 규칙을 가지고 있는지 알고 싶어 했습니다.

• 사면체 (4면체 모양): 원래의 피라미드 모양(4개의 손을 잡는 모양)이 무너질 때, 그것은 정확히 무작위 확률 게임처럼 무너졌습니다. 이는 "표준적인" 행동이었습니다.
• 고차 구조 (5, 6 또는 그 이상의 손): 새로운 밀도가 높은 모양들이 형성되고 연결될 때, 그것들은 표준 규칙을 깨뜨렸습니다. 그것들은 더 복잡하고 다른 규칙을 따랐습니다. 과학자들은 이를 "강성 퍼콜레이션(rigidity percolation)"이라고 부릅니다. 이는 단순히 군중이 무작위로 연결된 것이 아니라, 전체 구조를 갑자기 훨씬 더 딱딱하고 견고하게 만드는 방식으로 연결되었음을 의미합니다.

5. 핵심 요약

이 논문은 유리를 쥐어짤 때 단순히 크기가 줄어드는 것이 아니라, 내부 구조가 새로운 밀도가 높은 "상태"로 재편되는 극적인 상전이와 같은 사건이 일어난다고 결론짓습니다.

• 이 전이는 특정 "임계" 압력에서 발생합니다.
• 새로운 구조가 연결되는 방식은 무작위적인 확률(기존 모양의 경우)과 더 견고하고 구조적인 규칙(새로운 밀도 높은 모양의 경우)의 혼합입니다.
• 이러한 행동은 실리카 유리와 비정질 얼음에서 유사하게 나타나며, 이는 자연이 압력 하에서 서로 다른 종류의 유리질 물질을 재배열하기 위해 유사한 "설계도"를 사용함을 시사합니다.

요약하자면, 이 논문은 유리를 쥐어짤 때 미세한 "골격"이 어떻게 부서지고, 이동하고, 스스로를 재건하는지를 정확하게 그려내며, 느슨하고 흐느적거리는 유리에서 밀도가 높고 견고한 유리로 변하는 과정이 특정한 예측 가능한 임계점을 통해 일어난다는 것을 밝혀냈습니다.

기술 요약

기술 요약: 압축된 유리 내 비정질-비정질 전이의 퍼콜레이션 임계성

문제 정의
압축된 유리, 특히 사면체 기반 물질인 실리카($a$-SiO$_2$) 및 비정질 얼음($a$-H$_2$O)에서 발생하는 구조적 변형을 지배하는 메커니즘은 여전히 중요한 과제로 남아 있다. 결정질 상전이와 같은 날카로운 전이와 달리, 이러한 변환은 국부적인 구조적 모티프가 서로 다른 다형체(polyamorphs)로 점진적으로 진화하는 과정을 포함한다. 기존 연구들은 압력 하에서 고배위 다면체(high-coordination polyhedra)의 거대하고 퍼콜레이션되는 클러스터가 출현함을 확인했으나, 유리에서의 이러한 임계적 성질은 알려진 바가 없었다. 특히, 이러한 구조적 전이가 표준(무작위) 상관관계 없는 퍼콜레이션의 보편성 클래스(universality class)에 속하는지, 아니면 강성 퍼콜레이션(rigidity percolation)과 같은 별개의 임계 거동을 보이는지는 불분명하다. Ab initio 접근 방식에 내재된 작은 시스템 크기로 인해 기존의 임계 지수 추정 시도는 한계가 있었다.

방법론
이러한 한계를 해결하기 위해 저자들은 대규모 고전 분자 동역학(MD) 시뮬레이션을 활용하였다.

• 시스템 및 포텐셜: 실리카 유리에 대해서는 LAMMPS 패키지 내에서 SHIK 쌍 포텐셜을 사용하여, 비정질 얼음에 대해서는 Gromacs 내에서 TIP4P 힘의 장(force field)을 사용하여 시뮬레이션을 수행하였다. 시스템 크기는 실리카의 경우 약 1,000에서 $10^6$개 원자, 얼음의 경우 최대 49,000개 분량으로 구성되었다.
• 프로토콜: 샘플을 평형화한 후, 실리카의 경우 300 K, 얼음의 경우 124 K의 상온으로 냉각(quench)하였으며, 실리카는 35 GPa까지, 얼음은 15 kbar까지 준정적으로 압축하였다.
• 구조적 기술자(Descriptors): 두 가지 상호 보완적인 접근법이 사용되었다:
  1. 결합 모델(Bonded Model): 공유 결합된 Si-O 연결성을 기반으로 하여, SiO$_Z$ 다면체($Z$는 배위수)와 그 연결성(꼭짓점 공유, 모서리 공유, 면 공유)을 분석하였다.
  2. 비결합 모델(Non-bonded Model): 차단 반경(cutoff radius) 내에서의 중심 원자 배치(실리카의 경우 Si-Si, 얼음의 경우 O-O)를 기반으로 하며, SiSi$_Z$ 및 OO$_Z$ 다면체를 정의한다. 이 접근법은 상호작용을 비결합 힘으로 취급하여, 결합형 유리(실리카)와 비결합형 유리(얼음) 간의 직접적인 비교를 가능하게 한다.
• 퍼콜레이션 분석: 클러스터를 식별하기 위해 Union-Find 알고리즘을 구현한 전용 Python 코드인 Nexus-CAT가 사용되었다. 퍼콜레이션 특성인 상관 길이($\xi$), 질서 매개변수($P_\infty$), 평균 클러스터 크기($\langle S \rangle$), 최대 클러스터 크기($S_{max}$)가 계산되었다.
• 임계 지수: 임계 지수($\nu, \beta, \gamma$)와 프랙탈 차원($D_f, D_f^*$)을 구하기 위해 유한 크기 스케일링(finite-size scaling) 방법(데이터 붕괴 기법 사용)이 적용되었다.

주요 결과

1. 퍼콜레이션 전이의 순서: 압축 과정에서 연구는 뚜렷한 퍼콜레이션 이벤트의 순서를 확인하였다. $a$-SiO$_2$에서 저밀도(LD) 사면체 네트워크($Z=4$)는 고밀도(HD) 클러스터($Z=5, 6$)가 출현하고 퍼콜레이션함에 따라 탈퍼콜레이션(depercolate)된다. 구체적으로:

   • SiO$_5$ 클러스터는 $\approx$12 GPa에서 퍼콜레이션된다.
   • SiO$_6$ 클러스터는 $\approx$23 GPa에서 퍼콜레이션된다.
   • "스티쇼바이트 유사(Stishovite-like)" 클러스터(두 개의 모서리 공유 연결을 가진 SiO$_6$)는 $\approx$30 GPa에서 퍼콜레이션된다.
   • $a$-H$_2$O에서도 유사한 순서(LD, HD, VHD 클러스터)가 관찰되었으며, 이는 실리카보다 약 20배 낮은 압력에서 발생하였다. 이는 물질 특이적 상호작용에 의해 스케일링된 보편적 메커니즘을 시사한다.

2. 모델의 상호 보완성: 결합 모델(SiO$_Z$)과 비결합 모델(SiSi$_Z$/OO$_Z$) 접근법은 일관되고 상호 보완적인 통찰을 제공하였다. 비결합 모델은 구조적 변형을 성공적으로 포착하였으며, 결합형 유리(실리카)와 비결합형 유리(얼음)를 비교할 수 있는 통합된 프레임워크를 제공하였다.

3. 임계 지수 및 보편성:

   • 상관 길이 $\xi$는 시스템 크기 $L$에 대해 $\xi \propto L^a$ ($a \approx 1$)로 스케일링되어 임계 거동을 확인하였다.
   • 사면체의 탈퍼콜레이션: 저밀도 사면체 클러스터($Z=4$)의 탈퍼콜레이션은 표준 3D 무작위 퍼콜레이션의 임계 지수(예: 피셔 지수 $\tau \approx 2.12-2.14$)와 매우 유사한 결과를 나타냈다.
   • 고배위 클러스터의 퍼콜레이션: 이와 대조적으로, 고배위 클러스터($Z=5, 6$) 및 VHD 구조의 퍼콜레이션은 표준 퍼콜레이션으로부터 편차를 보였다. 프랙탈 차원($D_f$)은 표준 값인 2.53보다 체계적으로 낮았으며, 지수들은 다른 보편성 클래스를 시사했다.

4. 프랙탈 차원: 다양한 배위수에 대한 유한 클러스터(lattice animals)의 프래탈 차원 $D_f^*$는 $\approx 2.29$에서 $2.46$ 사이로 계산되었으며, 이는 일반적으로 표준 퍼콜레이션 값보다 낮았다.

의의 및 주장
본 논문은 장거리 기술자(long-range descriptors)의 도입과 대규모 시뮬레이션의 사용이 압축된 유리 내 퍼콜레이션 임계 지수의 최초의 견고한 추정을 가능하게 했다고 주장한다. 주요 의의는 비정질-비정질 전이가 균일하지 않은 임계 거동을 보인다는 발견에 있다:

• 본래의 사면체 네트워크의 붕괴는 표준 무작위 퍼콜레이션을 따르는 것으로 보인다.
• 고배위도의 강성 다면체 네트워크(밀집화와 관련된)의 출현 및 퍼콜레이션은 다른 보편성 클래스를 따를 가능성이 높으며, 이는 강성 퍼콜레이션 메커니즘을 시사한다. 이는 낮은 연결성을 가진 매질 내에서 높은 연결성을 가진 클러스터가 퍼콜레이션될 때 발생하는 위상학적 및 탄성적 불균질성에 기인한다.

본 연구는 이러한 퍼콜레이션 전이와 압축성 피크 및 소성(plasticity)의 시작과 같은 유리의 기계적 성질 사이의 연결 고리를 확립한다. 또한, 퍼콜레이션에 의한 구조적 변화를 추적하는 것이 결정 다형체(예: coesite IV, coesite V, stishovite)와 비정질 상태 사이의 유추를 위한 기초를 제공한다고 제안하며, 이는 유리 상태가 이러한 결정 구조와 연관된 에너지 메가베이슨(energy megabasins)에 대응할 수 있음을 시사한다. 이 연구는 결합형 및 비결합형 유리 시스템 모두에서 압력 유도 변형에 대한 이해를 통합하는 데 있어 비결합 접근법의 유용성을 강조한다.