Decoupling of single-particle and collective dynamics in arrested phase-separating glassy mixtures

기술 요약: 정지된 상분리 유리 혼합물에서의 단일 입자 및 집단 역학의 디커플링(Decoupling)

문제 제기
본 연구는 부드러운 스타 폴리머(soft star polymers)와 단단한 콜로이드 추적체(hard colloidal tracers)로 구성된 이성분 혼합물의 구조적 및 역학적 특성을 조사하며, 특히 부드러운 스타 폴리머가 주성분을 이루는 "단백질 한계(protein limit)" 내에서의 현상을 다룬다. 기존 연구들은 경질 구형(hard-sphere) 또는 인력 상호작용을 가진 유리 형성 액체를 광범위하게 규명해 왔으나, 부드러운 매트릭스 내에서의 추적체 역학은 상대적으로 덜 탐구되었다. 부드러운-단단한 혼합물에서 유리화(vitrification)와 초기 상분리(incipient phase separation, demixing) 사이의 상호작용을 지배하는 미시적 메커니즘을 이해하는 데에는 중요한 공백이 존재한다. 구체적으로, 저자들은 단단한 콜로이드의 첨가가 스타 폴리머 유리의 에르고드성(ergodicity), 구조, 그리고 유리 용융(glass melting)과 정지된 상분리 사이의 경쟁에서 발생하는 복잡한 다중 척도 역학에 어떻게 영향을 미치는지 밝히고자 한다.

방법론
저자들은 LAMMPS 패키지를 사용하여 거친 입자(coarse-grained) 분자 역학(MD) 시뮬레이션을 수행하였다.

• 모델 시스템: 시스템은 단단한 콜로이드(가파른 WCA 포텐셜로 모델링됨)와 스타 폴리머(기능성 $f$에 의존하는 부드럽고 초연성인 포텐셜로 모델링됨)로 구성된다. 단단한 성분과 부드러운 성분 사이의 교차 상호작용은 적분 방정식 이론으로부터 도출되었다.
• 파라미터: 연구는 스타 폴리머의 기능성 $f=166$ 및 크기 비율 $R_{gyr}^S / R_H = 2.25$에 초점을 맞추었다. 스타 폴리머 밀도는 순수 스타 시스템의 유리화 선보다 약간 높은 값인 $\rho_S \sigma_S^3 = 0.4$로 고정되었다.
• 시뮬레이션 프로토콜: 결정화를 방지하면서 평형 상태의 유리 상태를 생성하기 위해, 저자들은 스타 폴리머 직경에 다분산성(polydispersity)을 도입하였으며, 기능성 $f$를 액체 상태($f=10$)에서 목표 유리 상태($f=166$)까지 점진적으로 증가시키는 "퀜칭(quenching)" 프로토콜을 활용하였다. 여기에 7가지 서로 다른 농도의 단단한 콜로이드($\rho_H \sigma_S^3 \approx 0.02$ ~ $0.08$)를 첨가하였다.
• 분석: 구조적 특성은 정적 구조 인자($S_{ij}(q)$)와 조성-조성 구조 인자($S_{cc}(q)$)를 통해 분석되었다. 역학적 특성은 평균 제곱 변위(MSD), 자가(self) 및 집단(collective) 모드에 대한 중간 산란 함수(ISF), 반 호브(van Hove) 상관 함수, 그리고 비가우시안(non-Gaussian, NG) 파라미터 $\alpha_2(t)$를 사용하여 특징지어졌다. 시뮬레이션에서 관찰된 인구 분할(population splitting)을 해석하기 위해 분석적인 이상태 전환 확산 모델(two-state switching diffusion model)이 개발되었다.

주요 기여 및 결과

1. 유리 용융 및 초기 상분리:
   단단한 콜로이드를 스타 폴리머 유리에 첨가하면 부드러운 매트릭스의 역학이 가속화되어 유리를 액체 상태로 "녹이는" 효과를 낸다. 이는 MSD 플래토(plateau)의 단축, 확산 계수의 증가, 그리고 스타 폴리머의 이완 시간 감소를 통해 입증된다. 동시에, 구조 인자 $S_{SS}(q)$와 $S_{HH}(q)$는 $q \to 0$에서 피크를 발달시키며, 이는 장범위 농도 요동과 초기 상분리(탈혼합)의 시작을 신호한다. 저자들은 유리 용융과 초기 상분리가 비교적 낮은 단단한 콜로이드 농도에서 공존함을 보여준다.

2. 단일 입자 및 집단 역학의 디커플링(Decoupling):
   단단한 콜로이드의 자가 역학(self-dynamics)과 집단 역학(collective dynamics) 사이의 디커플링은 본 연구의 핵심 발견이다. 단단한 콜로이드는 MSD에서 확산적 거동(단순 브라운 운동을 시사)을 보이지만, 자가 중간 산란 함수($F_s$)는 중간 길이 척도에서 비지수적, 로그형 이완을 나타낸다. 반면, 집단 중간 산란 함수($F_c$)는 낮은 파수($q$)에서 훨씬 더 느리게 이완된다. 이러한 디커플링은 단단한 콜로이드 농도가 증가함에 따라 더욱 뚜렷해지며, 이는 주변의 느린 역학에 의해 집단 운동이 저해되는 크고 불균일한 영역이 형성되었음을 나타낸다.

3. 인구 분할 및 역학적 불균일성:
   정지된 상분리는 단단한 콜로이드를 콜로이드가 풍부한 영역과 콜로이드가 부족한 영역의 두 가지 뚜렷한 그룹으로 나누는 "인구 분할(population splitting)"을 야기한다.

• 공간적 불균일성: 국부 밀도 분포 분석은 콜로이드가 풍부한 영역이 콜로이드가 부족한 영역보다 현저히 더 이동성이 높음을 보여준다.
• 비가우시안 역학: 단단한 콜로이드는 강한 비가우시안 거동을 보인다. 부드러운 성분의 NG 파라미터가 전형적인 유리 역학(단일 피크)을 따르는 것과 달리, 단단한 성분의 NG 파라미터는 긴 시간대에서 두 번째 피크를 나타낸다. 이 두 번째 피크는 인구 분할에 기인한다.
• 이상태 모델: 저자들은 콜로이드가 빠른(콜로이드 풍부) 상태와 느린(콜로이드 부족) 확산 상태 사이를 교차하는 전환 확산 모델을 제안한다. 이 모델은 MSD는 확산적으로 보이지만 변위 분포는 두 인구 집단의 서로 다른 확산 계수로 인해 비가우시안으로 남는 "브라운 운동적이면서도 비가우시안인(Brownian yet non-Gaussian)" 역학의 출현을 성공적으로 설명한다.

4. 비가우시안성의 스케일링:
   가우시안 회복 특성 시간($\tau_G$)에 의해 시간을 재조정(rescaling)하면, 단단한 콜로이드의 NG 파라미터는 $t^{-0.8}$의 멱법칙(power-law) 붕괴를 가지는 보편적 곡선으로 수렴한다. 이 지수는 단순 이상태 모델이 예측하는 $t^{-1}$ 스케일링과 다르며, 이는 매트릭스의 고유한 유리성과 인구 분할 사이의 상호작용이 더 복잡한 역학적 지형을 생성함을 시사한다.

의의 및 주장
본 논문은 부드러운-단단한 혼합물에서 정지된 상분리와 유리성이 상호작용하여 표준 평균장 기술이나 단순 추적자 확산 모델로는 포착할 수 없는 복잡한 다중 척도 현상을 생성한다는 것을 입증한다고 주장한다. 구체적으로, 본 연구는 다음을 강조한다:

• 부드러운 유리 매트릭스 내의 단단한 콜로이드는 자가 역학 및 집단 역학 사이의 디커플링을 드러내는 프로브 역할을 하며, 이는 상분리 정지의 징후이다.
• 시스템은 콜로이드의 국부 환경에 따라 빠른 인구와 느린 인구로 나뉘는 "인구 분할"을 보이며, 이는 전체적인 운동이 확산적으로 보임에도 불구하고 비가우시안 역학을 유발한다.
• 유리 용융과 상분리의 공존은 구조적 불균일성(콜로이드 풍부 vs 콜로이드 부족 도메인)이 역학적 불균일성을 직접적으로 결정하는 불균일한 환경을 조성한다.

저자들은 이러한 발견이 이러한 혼합물에서 관찰되는 복잡한 유변학적 거동에 대한 미시적 설명을 제공하며, 로그형 ISF 붕괴 및 NG 파라미터의 특정 스케일링과 같은 구체적인 징후들이 향후 스타 폴리머-콜로이드 혼합물에 대한 실험을 통해 검증될 수 있을 것이라고 제언한다.