2차원 이방성 입자에서 형태가 결정하는 고체‑고체 상전이 동역학 경로

초록

본 연구는 2차원 볼‑스틱 다각형(오각형, 육각형, 팔각형) 시스템을 분자동역학 시뮬레이션으로 조사하여, 동일한 격자구조를 유지하는 등축 전이 동안 전단(번역) 운동은 균일 팽창만을 보이는 반면, 회전(방향) 운동은 입자 형태에 따라 전혀 다른 결함 패턴을 형성한다는 점을 밝혀냈다. 고정‑COM 시뮬레이션을 통해 오각형은 회전‑주도, 육각형은 전단‑주도, 팔각형은 두 운동이 동시 진행되는 준평형 경로를 따른다. 이러한 경로 차이는 전이 속도와 냉각 시 트랩 형성에 영향을 미치며, 물질 설계 시 원하는 동역학 특성을 조절하는 새로운 설계 원칙을 제시한다.

상세 분석

이 논문은 기존의 구형 콜로이드 모델이 다루기 어려운 ‘이방성 입자’의 고체‑고체 상전이 동역학을 이해하기 위해, 각 꼭짓점에 LJ 구와 인접 구를 연결하는 결합을 갖는 2차원 볼‑스틱 다각형 모델을 도입하였다. 오각형, 육각형, 팔각형 모두 동일한 압력·온도 조건에서 close‑packing → rotator‑crystal 로의 등축 전이를 보이며, 격자 대칭은 변하지 않는다. 전이 과정에서 Voronoi 면적이 균일하게 증가함을 통해 전단 운동이 단순 팽창임을 확인하였다. 반면, 각 입자의 몸체‑방향을 나타내는 단위벡터 필드에서 급격한 방향 변화를 보이는 ‘결함’들을 추적했을 때, 오각형은 흐릿한 스트라이프 형태, 육각형은 완전 무작위, 팔각형은 뚜렷한 스트라이프를 형성한다는 점이 드러났다.

핵심 메커니즘을 밝히기 위해 ‘고정‑COM’ 시뮬레이션을 수행하였다. 이는 전단 좌표를 고정하고 회전만 허용함으로써, 주어진 부피에서 회전 자유도가 자유에너지 최소화에 얼마나 기여하는지를 측정한다. 결과는 세 가지 전형적인 커플링 모드를 제시한다. 첫째, 오각형에서는 고정‑시뮬레이션이 원래 MD보다 약간 높은 몸체‑방향 순서를 유지하므로 회전이 전단보다 빠르게 진행되어 회전‑주도 경로를 따른다. 둘째, 육각형은 고정‑시뮬레이션이 원래 MD보다 낮은 순서를 보이며, 전단 팽창이 회전보다 앞서 진행되어 전단‑주도 경로가 지배한다. 셋째, 팔각형은 두 시뮬레이션이 거의 일치해 전단과 회전이 동시 진행되는 준평형(quasi‑equilibrium) 경로를 따른다.

이러한 차이는 입자 간 거리와 회전 제약에 기인한다. 육각형은 팽창 후 입자 간격이 커 회전이 독립적으로 일어나 무작위 결함을 만든다. 반면, 오각형과 팔각형은 초기 간격이 작아 회전이 이웃 입자에 전이되기 쉬워 스트라이프 형태의 결함 전파가 나타난다. 특히 오각형은 이미 부모상에 소수의 결함이 존재해 스트라이프 외에 고립된 결함도 관찰된다.

전이 속도 측면에서, 회전‑주도 경로는 회전 장벽이 존재하므로 전단‑주도보다 느리지만, 준평형 경로는 두 운동이 동시에 진행돼 가장 빠른 전이를 보인다. 냉각 과정에서는 동일 메커니즘이 적용되지만, 회전‑주도와 전단‑주도 경로가 혼재해 다양한 ‘동역학 트랩’이 형성되어 역전이 속도가 크게 변동한다.

이 연구는 입자 형태가 전이 동역학을 결정짓는 핵심 변수임을 실증적으로 보여주며, 회전‑전단 커플링을 설계 변수로 활용해 전이 속도와 경로를 맞춤 설계할 수 있는 새로운 설계 패러다임을 제시한다.