Hydrodynamic stability and pattern formation in hexatic epithelial layers

기술 요약: 육방정 격자 상피층의 유체역학적 안정성 및 패턴 형성

문제 제기
본 논문은 변형 가능한 세포들의 조밀한 충전(close-packing)에서 발생하는 육방정 질서(6회 회전 대칭성)의 역할을 중심으로, 상피층의 유체역학적 안정성과 패턴 형성을 다룬다. 능동적 액정 유체 역학은 극성(1회 대칭) 및 네마틱(2회 대칭) 시스템에 광범위하게 적용되어 집단 세포 이동과 결함 역학을 설명해 왔으나, 육방정 상피의 구체적인 유체역학적 거동에 대한 연구는 미진한 상태이다. 핵심 과제는 기계적 활성, 측면 접착력, 그리고 육방정 질서의 다중 스케일 특성 사이의 상호작용이 정지 상태(quiescent state)의 안정성과 집단 운동의 출현에 어떻게 영향을 미치는지 이해하는 것이다. 불안정성이 일반적으로 발생하는 표준 능동 네마틱 시스템과 달리, 육방정 시스템은 서로 다른 길이 스케일에 걸친 에너지 주입과 소산 사이의 더 복잡한 균형을 보일 수 있다.

방법론
저자들은 수동적 $p$-atic 액정 이론으로부터 유도된 연속체 유체역학 모델을 능동적 육방정 시스템($p=6$)으로 확장하여 사용한다.

• 지배 방정식: 모델은 질량 밀도, 운동량 밀도, 그리고 국소적 육방정 방향을 나타내는 계수-6 차수 텐서($Q_6$)를 결합한다. 능동 응력은 세포 수준에서 생성되는 힘을 반영하여 $\sigma^{(a)} = \alpha_6 \nabla^{\otimes 4} \odot Q_6$로 모델링된다.
• 차원 축소: 연구는 두 가지 영역에 집중한다:
  1. 채널 구속(Channel Confinement): 폭 $L$을 가진 2D 채널로, 시스템은 종방향에 대해 불변이다. 이는 육방정 각도 $\theta(y)$와 종방향 속도 $v_x(y)$에 대한 결합 방정식으로 문제를 단순화한다.
  2. 주기적 도메인: 완전한 벡터 흐름과 방향 역학을 허용하는 2D 주기적 도메인이다.
• 소산의 정규화: 표준 점성 모델이 작은 스케일에서 보이는 비물리적인 불안정성을 해결하기 위해 저자들은 하이퍼점성(hyperviscosity, 구체적으로 고차 속도 구배인 $\eta_6 \nabla^4$)을 도입한다. 이는 에너지 소산 스케일이 작은 스케일에서의 에너지 주입율($\sim q^6$)과 일치하도록 하여 발산을 방지하기 위함이다. 이로 인해 표준 전단 점성($\eta_2$)과 하이퍼점성($\eta_6$)을 모두 포함하는 스케일 의존적 소산율이 도출된다.
• 분석: 저자들은 정지 상태에 대한 선형 안정성 분석을 수행하여 분산 관계를 도출한다. 이들은 데보라 수($De$)와 교차 길이 스후$\ell_\times = (\eta_6/\eta_2)^{1/4}$를 포함한 임계 파라미터들을 식별하며, 이는 점성 소산 영역과 하이퍼점성 소산 영역 사이의 전이를 결정한다. 전체 비선형 방정식의 수치 시뮬레이션을 사용하여 해석적 예측을 검증하고 전이 후의 역학을 탐구한다.

주요 기여 및 결과

1. 스케일 의依存적 안정성 메커니즘:
   본 논문은 육방정 상피의 유체역학적 안정성이 (능동 네마틱에서처럼) 단순히 시스템 크기에 의해 결정되는 것이 아니라, 길이 스케일에 따른 에너지 주입과 소산 사이의 균형에 의해 결정됨을 입증한다.

• 일반적 안정성: 측면 접착 상호작용( $\eta_2$로 표현됨)이 강할 때, 시스템은 크기에 관계없이 일반적으로 안정적이다. 에너지가 주입되는 스케일에서 효율적으로 소산되기 때문에, 정지 상태는 미세한 섭동에 대해 안정성을 유지한다.
• 불안정성 임계치: 불안정성은 활성이 소산에 비해 충분히 높을 때만 발생하며, 이는 교차 스케일 $\ell_\times$에서 정의된 임계 데보라 수($De_\times > 16/3$)에 의해 특징지어진다.

2. 패턴 형성 불안정성의 계층 구조:
   구속된 채널 내에서 시스템은 역방향 흐름 레인(counter-flowing lanes) 형성을 유도하는 일련의 불안정성 폭포를 보인다.

• 자발적 흐름 전이: 능동 네마틱과 유사하게, 1차 불안정성은 세포들이 서로 반대 방향으로 흐르는 두 개의 레인을 형성한다.
• 레인 증식: 네마틱과 달리, 다수의 고유한 길이 스케일의 존재는 추가적인 계층적 불안정성을 허용한다. 채널 폭 $L$이 증가하거나 활성이 높아짐에 따라, 시스템은 증가하는 수의 레인($n_0 \approx 2\lfloor q_0/q_{min} \rfloor$)으로 분할된다.
• 고정된 레인 폭: 결정적으로, 이러한 레인의 폭은 시스템 크기에 비례하여 스케일링되는 것이 아니라, 교차 길이 $\ell_\times$에 비례하는 스케일 독립적인 값으로 수렴한다. 이는 레인 구조가 재료의 류학(rheology)에 따른 본질적인 특성임을 시사한다.

3. 2차원 역학 및 "능동 난류(Active Turbulence)":
   2D 주기적 도메인에서 1D 레인 구조는 와류(vortex)의 증식에 의해 불안정해진다.

• 낮은 활성: 시스템은 파장 $\sim \ell_\times$를 갖는 정지된 질서 띠(헤링본 패턴)를 형성한다.
• 높은 활성: 시스템은 시공간적 혼돈(spatiotemporal chaos) 상태로 전이된다. 이 영역은 네마틱에서 관찰되는 "능동 난류"와 유사한 특징을 공유하지만, 저자들은 뚜렷한 차이점을 언급한다. 혼돈 상태는 선형 불안정성의 "기억"을 유지하며, 특정 악성 중간엽 종양(예: 섬유육종)에서 관찰되는 헤링본 패턴과 유사한 방향성 클러스터 및 와류 구조를 보여준다.

의의 및 주장
저자들은 이 연구가 육방정 상피의 유체역학적 안정성에 대한 최초의 포괄적인 분석을 제공한다고 주장한다. 주요 의의는 육방정 질서가 다중 스케일 조직화를 도입하여 능동 유체의 안정성 기준을 근본적으로 변화시킨다는 점을 확립한 데 있다.

• 네마틱과의 대비: 본 논문은 벌크 능동 네마틱과 극명한 대조를 강조한다. 네마틱에서는 도메인이 충분히 크다면 어떤 유한한 활성이라도 시스템을 불안정하게 만들 수 있다. 반면, 육방정 상피에서는 강한 측면 접착(높은 $\eta_2$)이 시스템 크기에 관계없이 정지 상태를 영구적으로 안정화할 수 있다.
• 생물학적 관련성: 이 모델은 상피 조직의 기계적 안정성이 측면 접착(예: E-cadherin 발현)을 통해 조절될 수 있음을 시사한다. 고정된 폭의 레인과 헤링본 패턴의 출현은 복잡한 세포 수준의 규칙을 도입하지 않고도 이동하는 세포 클러스터의 공간적 조직화와 특정 병리학적 양상을 이해하기 위한 이론적 틀을 제공한다.
• 이론적 프레임워크: 본 연구는 $p$-atic 질서와 스케일 의존적 소산 사이의 상호작용이 복잡한 집단 행동의 출현을 지배한다는 연속체 프레임워크를 구축하며, 이는 극성 또는 네마틱 능동 물질에서 전형적인 단일 질서 매개변수 설명을 넘어선다.