[논문 리뷰] Variational formulation of active nematics: theory and simulation

Waleed Mirza, Alejandro Torres-Sànchez외 2명|arXiv (Cornell University)|2023. 06. 02.

Micro and Nano Robotics인용 수 5

한 줄 요약

이 논문은 온사거의 형식론을 기반으로 압축 가능한 2차원 활성 난류에 대한 변분 프레임워크를 제안하며, 자유 에너지, 소산, 활성성을 통합하여 제어 방정식과 호환성 있는 유한요소 방법을 유도한다. 이는 밀도에 따라 변하는 활성 난류 역학의 정확한 시뮬레이션을 가능하게 하며, 수렴성은 뛰어나고 상처 치유와 제한된 세포 집단에서의 결함을 통한 자가조직화를 드러낸다.

ABSTRACT

The structure and dynamics of important biological quasi-two-dimensional systems, ranging from cytoskeletal gels to tissues, are controlled by nematic order, flow, defects and activity. Continuum hydrodynamic descriptions combined with numerical simulations have been used to understand such complex systems. The development of thermodynamically consistent theories and numerical methods to model active nemato-hydrodynamics is eased by mathematical formalisms enabling systematic derivations and structured-preserving algorithms. Alternative to classical nonequilibrium thermodynamics and bracket formalisms, here we develop a theoretical and computational framework for active nematics based on Onsager's variational formalism to irreversible thermodynamics, according to which the dynamics result from the minimization of a Rayleighian functional capturing the competition between free-energy release, dissipation and activity. We show that two standard incompressible models of active nemato-hydrodynamics can be framed in the variational formalism, and develop a new compressible model for density-dependent active nemato-hydrodynamics relevant to model actomyosin gels. We show that the variational principle enables a direct and transparent derivation not only of the governing equations, but also of the finite element numerical scheme. We exercise this model in two representative examples of active nemato-hydrodynamics relevant to the actin cytoskeleton during wound healing and to the dynamics of confined colonies of elongated cells.

연구 동기 및 목표

• 밀도 변화와 전환을 고려한 압축 가능한 2차원 활성 난류 시스템을 위한 열역학적으로 일관된 프레임워크를 개발하기 위해.
• 비가역 열역학에 기반한 변분 원리에서 유도된 활성 난류 모델의 물리적 해석을 명확히 하기 위해.
• 제어 편미분 방정식과 호환성 있는 유한요소 수치적 스킴을 체계적으로 유도하기 위해.
• 액틴 세척모계의 상처 치유 동안의 역학과 제한된 긴장된 세포 집단에서의 결함 진화와 같은 생물학적으로 관련된 시나리오를 시뮬레이션하기 위해.
• 메쉬 크기와 시간 간격에 대한 수렴성 분석을 통해 수치 방법을 검증하기 위해.

제안 방법

• 시간 진동이 자유 에너지 변화, 소산 및 활성성을 조합한 기능을 최소화하는 방식으로 온사거의 변분 원리를 사용하여 동역학을 기술한다.
• 밀도 ρ, 난류 순서 매개변수 S, 그리고 지향장 n을 포함한 상태 변수와 관련된 흐름 및 변화율을 도입한다.
• 자유 에너지, 소산, 활성 일 항을 포함하는 라그랑지안에 변분 원리를 적용하여 제어 방정식을 유도한다.
• 변분 구조를 그대로 유지하는 유한요소법(FEM)을 구성하여 에너지 안정성과 일관된 공간 이산화를 확보한다.
• 시간에 암시적이고 갈레르킨 기반의 FEM 스킴을 사용하며, 메쉬 크기와 시간 간격에 대한 수렴성 분석을 통해 검증한다.
• 이 프레임워크를 두 가지 생물학적 시스템에 적용한다: 큰 알세포에서의 상처 치유와 제한된 세포 집단에서의 결함 역학.

실험 결과

연구 질문

• RQ1 밀도에 따라 변하는 효과를 고려한 압축 가능한 2차원 활성 난류에 대해 통합적이고 열역학적으로 일관된 프레임워크를 어떻게 개발할 수 있는가?
• RQ2 온사거의 변분 형식론을 사용할 경우 활성 난류에서 운동 방정식과 수치 스킴을 동시에 유도하는 데 어떤 의미가 있는가?
• RQ3 압축성과 전환을 포함할 경우 활성 시스템에서 난류 구조의 자가조직화에 어떤 영향을 미치는가?
• RQ4 위상 결함과 활성 유동은 큰 알세포의 상처 치유 역학에서 어떤 역할을 하는가?
• RQ5 활동 강도와 이방성은 제한된 활성 난류 시스템에서 결함 역학과 유동 패턴에 어떻게 영향을 미치는가?

주요 결과

• 변분 프레임워크는 기존의 비압축성 활성 난류 모델을 성공적으로 복원하면서도, 밀도에 따라 변하는 시스템으로 확장한다.
• 유한요소 스킴은 최적의 공간 수렴성(비율 ∼h²)과 시간 수렴성(비율 ∼Δt)을 보이며, 수치적 안정성을 확인한다.
• 상처 치유 모델에서는 활성 난류 유동이 결함의 이동과 재조직화를 이끌어내며, 큰 세포에서 세척모계 재조직화를 촉진한다.
• 제한된 집단에서는 활동이 증가할수록 활성 난류가 발생하며, 이때 결함 역학은 ℓₐ/ℓ₀ 비율에 의해 지배되는데, 여기서 ℓₐ는 활성 난류 길이 척도이다.
• 이방성 활동(κ = λ_aniso/λ로 정량화)은 결함 운동과 유동 패턴을 크게 변화시키며, κ가 클수록 유도된 유동과 결함 쌍이 형성된다.
• 모델은 난류 질서, 활동성, 유체역학적 유동의 상호작용에 의해 아스터, 타코이드, 난류 벨트의 자가조직화에 의해 발생하는 현상을 캐릭터라이즈한다.