Thinning-by-spinning: shear rheology of dense chiral fluids

기술 요약: 스피닝에 의한 박화(Thinning-by-spinning): 고밀도 카이랄 유체의 전단 유변학

문제 제기
자기 회전(self-spinning) 입자로 구성된 카이랄 유체는 비평형 상태의 물질을 구현하는 다재다 рука한 플랫폼이다. 여기서 회전은 홀(odd, 횡방향) 힘을 생성하여 상 거동과 역학을 재형성한다. 자기 회전이 육방정계 질서(hexatic order)를 깨뜨리고 결함 증식을 촉진함으로써 2D 고체를 유동화한다는 사실은 이미 확립되어 있으나, 자기 회전과 국소 전단(local shear) 사이의 동등성, 그리고 이러한 카이랄 유체가 외부에서 가해진 전단에 어떻게 반응하는지에 대한 근본적인 질문이 남아 있다. 특히, 밀도가 높은 계(고체 및 유체 영역을 포괄하는)에서 고유한 카이랄성과 외부 구동 사이의 상호작용은 통일된 유변학적 프레임워크를 필요로 한다.

방법론
저자들은 2차원 입자 기반 시뮬레이션을 사용하여 모델 카이랄 유체의 선형 및 비선형 유변학을 조사한다. 시스템은 레나드-존스(Lennard-Jones) 퍼텐셜에 비보존적 횡방향 힘이 추가된 $N=4096$개의 구형 입자(디스크)로 구성된다. 역학은 다음을 포함하는 과감쇠 랑제빈 방정식(overdamped Langevin equation)에 의해 지배된다:

1. 보존적 상호작용: 레나드-존스 퍼텐셜 $U(r_{ij})$로부터 유도됨.
2. 횡방향 힘: $f_{ij} = -\Gamma u_{ij}$로 정의되며, 여기서 $u_{ij}$는 분리 벡터와 평면에 수직인 단위 벡터의 외적에 비례하며, 카이랄 파라미터 $\Omega$에 의존하는 토크를 도입한다.
3. 외부 전단: 전단율 $\dot{\gamma}$를 가진 리스-에드워즈(Lees-Edwards) 경계 조건을 통해 구현됨.

시뮬레이션은 다양한 전단율($10^{-5} \leq \dot{\gamma} \leq 1$)과 카이랄 강도($\Omega$)를 탐색하며, 양(+)의 값(전단 와도와 정렬된 방향)과 음(-)의 값(전단 와도와 반대 방향)을 모두 포함한다. 시스템은 고체 영역(고밀도 $\rho=0.7$, 저온 $k_BT=0.35$)과 액체 영역(밀도 $\rho=0.6$, $k_BT=0.47$) 모두에서 분석된다.

주요 관측량은 다음과 같다:

• 응력 텐서: 이링-커크우드(Irving-Kirkwood) 응력 텐서 $\sigma_{\alpha\beta}$로부터 대칭화된 전단 응력 $\sigma_s$와 점도 $\eta$를 도출한다.
• 구조적 질서: 상전이와 결함 증식을 모니터링하기 위해 국소 및 전역 육방정계 질서 파라미터($\psi_6$ 및 $\Psi_6$)를 사용한다.
• 유효 온도($T_{eff}$): 플럭추에이션-소산 정리(FDT)의 위반, 특히 스토크스-아인슈타인 관계식 $k_B T_{eff} = D/\mu$를 통해 계산된다. 여기서 $D$는 확산 계수이고 $\mu$는 작은 외부 힘에 대한 선형 응답으로부터 얻은 이동도이다.

주요 기여 및 결과

1. 고체의 유동화 및 전단 박화(Shear Thinning):
   고체 영역에서 시스템은 항복 응력(yield stress) 거동을 보인다. 즉, 평형 상태의 시스템($\Omega=0$)은 낮은 전단율에서 흐르지 않는다. 그러나 카이랄 파라미터 $|\Omega|$를 증가시키면 시스템이 유동화되어 제로 전단 한계에서 유한한 점도를 갖게 된다. 전단율과 카이랄성을 모두 높이면 점도가 감소하며, 시스템을 전단 박화 영역으로 몰아넣는다. 점도 스케일링은 $\eta \sim 1/\dot{\gamma}$ (모드 결합 이론과 일치) 및 $\eta \sim 1/|\Omega|$를 따르며, 이는 카이랄성이 고유한 전단 응력원 역할을 함을 나타낸다.

2. 스피닝에 의한 박화 메커니즘:
   저자들은 고정된 전단율에서도 카이랄 파라미터 $|\Omega|$의 크기에 따라 점도가 단조 감소하는 "스피닝에 의한 박화(thinning-by-spinning)" 메커니즘을 확인하였다. 이 효과는 거시적이지만, 순수하게 미시적인 쌍 단위 카이랄 상호작용에 의해 구동된다.

3. 유효 온도를 통한 통일된 프레임워크:
   선형 응답 영역에서 저자들은 그린-쿠보(Green-Kubo) 관계의 일반화를 제안한다. 열역학적 온도 $T$를 카이랄성에 의존하는 유효 온도 $T_{eff}(\Omega)$로 대체함으로써, 평형 그린-쿠보 적분이 비평형 전단 점도를 정량적으로 예측할 수 있음을 보여준다. 이는 카이랄성이 시스템에 에너지를 주입하는 방식이 거시적 척도에서 높아진 유효 온도로 포착될 수 있음을 시사한다.

4. 유동 곡선의 붕괴(Flow Curve Collapse):
   선형 응답을 넘어, 서로 다른 전단율과 회전율에 대한 유동 곡선은 점도를 $T_{eff}$로 스케일링하고 전단율을 $|\Omega|$로 스케일링했을 때 하나의 마스터 커브로 붕괴된다. 이 붕괴는 자기 회전이 외부 변형과 직접 경쟁하는 고유한 응력원으로서 작용하여, 전단과 스피닝의 유동화 경로를 통합함을 입증한다.

5. 손잡이 방향 비대칭성 및 층상 구조(Layering):
   강한 구동 조건에서는 스피닝과 전단 사이의 대응 관계가 깨지며 뚜렷한 손잡이 방향 비대칭성이 드러난다.

• 반대 방향 손잡이 ($\Omega < 0$): 횡방향 힘이 가해진 전단에 저항할 때, 시스템은 임계 전단율 $\dot{\gamma}^*$에서 구조적 전이를 겪는다. 유체는 실 모양의 준-1차원(quasi-1D) 질서 도메인(층상 구조)으로 자가 조직화된다. 이 구조는 흐름을 용이하게 하여 응력을 크게 완화하고 전단 박화를 강화한다.
• 정렬된 손잡이 ($\Omega > 0$): 시스템은 전단 방향을 따라 중간 정도로 길어진 육방정계 패치를 보이지만, 훨씬 높은 전단율에 도달하기 전까지는 앞서 언급한 뚜렷한 실 모양의 채널을 형성하지 못한다.
• 층상 구조로의 전이는 평균 육방정계 질서 파라미터 $|\psi_6|$의 최솟값과 구조 인자 피크(structure factor peak)의 급격한 증가에 의해 특징지어진다.

의의
본 논문은 카이랄성이 밀집된 입자계에서 유동화를 일으키는 일반적인 메커니즘임을 확립한다. 저자들은 카이랄성과 외부 구동이 서로를 강화하여 고체 및 유체 상 모두에서 전단 박화 거동을 만들어내는 통일된 유변학적 프레임워크를 제공한다. 연구 결과는 선형 응답가 유효 온도로 설명될 수 있는 반면, 비선형 영역은 외부 변형과 고유한 카이랄성 사이의 경쟁에 의해 지배되며, 상대적인 구동 방향(handedness)에 따라 층상 구조와 같은 뚜렷한 구조적 상태를 형성함을 보여준다. 이러한 결과는 밀집된 카이랄 현탁액, 능동 물질(active matter), 그리고 횡방향 상호작용이 재료의 특성을 근본적으로 변화시키는 홀 역학적(odd mechanical) 메타물질의 유변학을 이해하는 데 중요한 참조가 된다.