Non-linear visco-elasto-plastic rheology of a viscous vertex model

1. 연구 배경 및 문제 제기 (Problem)

• 생물학적 조직의 역학: 발생 생물학, 조직 공학, 암 전이 등 복잡한 생물학적 과정에서 조직의 형태 변화는 세포 집단의 역학적 상호작용에 의해 결정됩니다.
• 기존 연구의 한계:
  • 기존 점 모델 연구들은 주로 선형 유체역학 (linear rheology) 이나 탄성 특성에 집중했습니다.
  • 많은 모델이 외부 기질 (substrate) 에 대한 마찰을 가정하거나, 내부 소산 (dissipation) 과정을 명시적으로 다루지 않았습니다.
  • 조직 내부의 소산 과정은 병진 (translation) 및 회전 (rotation) 갈릴레이 불변성 (Galilean invariance) 을 만족해야 하지만, 이를 동시에 충족하는 모델은 드뭅니다.
  • 비선형 탄성, 점소성 (visco-plastic) 거동, 그리고 능동적 스트레스 (active stress) 가 결합된 비선형 유체역학적 관계를 체계적으로 규명한 연구는 부족했습니다.

2. 방법론 (Methodology)

저자들은 다음과 같은 방법론을 사용하여 모델을 구축하고 검증했습니다.

• 점 모델 (Vertex Model) 확장:
  • 2D 조직을 다각형 셀의 타일링으로 표현하며, 에너지 함수는 면적 탄성과 둘레 탄성을 포함합니다.
  • 내부 마찰 (Internal Friction): 외부 기질 마찰 대신, 세포 면적, 세포 둘레, 세포 간 인터페이스 길이에 대한 점성 마찰 항을 명시적으로 도입했습니다. 이는 병진 및 회전 갈릴레이 불변성을 만족하도록 설계되었습니다.
  • 능동적 힘 (Active Forces): 두 가지 다른 세포 기작을 모델링했습니다.
    1. 이방성 인터페이스 장력 (Anisotropic interface tension): 방향에 따라 장력이 변하는 힘.
    2. 능동적 기어오르기 (Active crawling): 세포 내 특정 꼭짓점에 작용하는 인장력 (Xenopus 배아에서 관찰된 수렴 - 신장 기작 모사).
• 변형률 분해 (Strain Decomposition):
  • 전체 전단 변형률 ($\tilde{v}$) 을 **가역적인 세포 모양 변화 ($\dot{Q}$)**와 **비가역적인 세포 재배열 (T1 전이, $R$)**로 분해하는 프레임워크를 적용했습니다.
  • 이를 통해 조직의 전단 응력 ($\tilde{\sigma}$) 을 세포 모양 ($Q$) 과 전단율 ($\tilde{v}$) 의 함수로 표현하는 구성 방정식을 유도했습니다.
• 시뮬레이션:
  • 작은 진폭의 진동 전단 (linear regime) 과 일정한 전단율, 그리고 큰 진폭의 진동 전단 (large-amplitude oscillatory shear) 시뮬레이션을 수행하여 모델의 예측력을 검증했습니다.

3. 주요 기여 및 결과 (Key Contributions & Results)

A. 비선형 전단 응력 구성 방정식 ($\tilde{\sigma}(Q, \tilde{v})$)

• 선형 영역: 작은 변형에서 저장 탄성률 ($G_0$) 과 점도 ($\eta_0$) 가 일정함을 확인했습니다.
• 비선형 탄성 영역: 큰 변형에서 전단 응력은 세포 모양 $Q$에 대해 3 차 및 5 차 다항식 형태로 비선형적으로 증가함을 발견했습니다.
• 점성 거동의 재발견:
  • 조직이 세포 모양 변화만 일으킬 때와 T1 전이 (세포 재배열) 가 발생할 때의 점도가 크게 다름을 발견했습니다.
  • T1 전이가 발생하는 경우 (소성 유동), 유효 점도는 탄성 변형만 일어날 때보다 약 3 배 더 큽니다. 이는 T1 전이 직후 세포 모양의 이완 과정에서 추가적인 전단 응력이 발생하기 때문으로 해석됩니다.
• 능동적 힘의 영향: 능동적 장력이나 기어오르기 힘은 전단 응력 곡선을 선형적으로 이동시키는 효과를 가지며, 이는 각각의 활동성 계수 ($\rho_{ait}, \rho_{crawl}$) 로 정량화되었습니다.

B. 재배열 함수 ($R(Q, \tilde{v})$) 및 항복 거동

• 비가역 분율: 세포 재배열에 의한 전단율 기여도 ($R/\tilde{v}$) 는 세포 모양 $Q$에 대해 지수 함수 형태로 증가함을 발견했습니다.
• 활동성의 영향: 능동적 힘은 재배열이 시작되는 임계 세포 모양 ($Q_0$) 을 변화시키지만, 지수 함수의 기울기 ($b$) 에는 큰 영향을 미치지 않았습니다.
• 예측 모델: 구해진 $R(Q, \tilde{v})$ 함수를 사용하여, 일정한 전단율 하에서의 정상 상태 세포 모양 ($Q_s$) 을 매우 정확하게 예측할 수 있었습니다.

C. 큰 진폭 진동 전단 (LAOS) 예측

• 개발된 평균장 모델을 사용하여 큰 진폭의 진동 전단 시뮬레이션 결과를 예측했습니다.
• 모델은 세포 모양 ($Q$) 과 전단 응력 ($\tilde{\sigma}$) 의 시간적 진화를 실험 (시뮬레이션) 데이터와 잘 일치시켰습니다.
• 전향 방향이 바뀌는 순간 (yielding) 에 약간의 편차가 관찰되었으나, 이는 과거 변형 이력에 대한 기억 효과 (memory effect) 가 평균장 모델의 단순한 $Q$ 변수로 완전히 포착되지 않았기 때문으로 분석되었습니다.

4. 의의 및 결론 (Significance)

• 이론적 통합: 세포 수준의 기계적 상호작용 (마찰, 능동력) 과 조직 수준의 비선형 유체역학적 거동을 연결하는 최초의 체계적인 평균장 모델을 제시했습니다.
• 내부 소산의 중요성: 외부 기질 마찰이 아닌 내부 점성 마찰을 고려함으로써, 실제 3D 조직이나 기질이 없는 환경에서의 조직 역학을 더 정확하게 묘사할 수 있음을 보였습니다.
• 유체역학적 복잡성 규명: 조직의 점도가 단순히 물질 상수가 아니라, 세포 재배열 (T1 전이) 의 유무에 따라 동적으로 변화한다는 점을 규명했습니다. 이는 조직이 탄성체에서 소성체로 전이되는 과정을 이해하는 데 중요한 통찰을 제공합니다.
• 일반성: 이 접근법은 다양한 세포 기반 조직 모델 (Vertex model variants) 에 적용 가능한 일반적인 구성 방정식을 구축하는 틀을 제공합니다.

이 연구는 생물학적 조직이 어떻게 복잡한 비선형, 비선형 점탄소성 거동을 보이는지에 대한 물리적 메커니즘을 규명하고, 이를 정량적으로 예측할 수 있는 강력한 수학적 도구를 제공했다는 점에서 의의가 큽니다.