[论文解读] Variational formulation of active nematics: theory and simulation

Waleed Mirza, Alejandro Torres-Sànchez+2 more|arXiv (Cornell University)|Jun 2, 2023

Micro and Nano Robotics被引用 5

一句话总结

本文基于昂萨格形式理论，提出了一种用于可压缩二维活性向列系统的变分框架，统一了自由能、耗散与活性，推导出控制方程及相容的有限元格式。该方法可精确模拟密度依赖的活性向列动力学，展现出良好的收敛性，并揭示了伤口愈合与受限细胞群落中的缺陷介导自组织现象。

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.

研究动机与目标

• 开发一种适用于具有密度变化和周转的可压缩二维活性向列系统、热力学一致的理论框架。
• 通过基于非平衡热力学的变分原理推导，阐明活性向列模型的物理解释。
• 系统推导控制偏微分方程及相容的有限元数值格式。
• 模拟生物相关情景，包括伤口愈合过程中肌动蛋白细胞骨架的动力学，以及受限伸长细胞群落中的缺陷演化。
• 通过空间与时间离散化的收敛性分析，验证数值方法的可靠性。

提出的方法

• 采用昂萨格变分原理描述动力学，其中时间演化通过最小化包含自由能变化、耗散与活性的泛函实现。
• 引入状态变量：密度 ρ、向列序参量 S 和指向场 n，及其对应的通量与变化率。
• 通过在包含自由能、耗散与活性功项的拉格朗日函数上应用变分原理，推导出控制方程。
• 构建继承变分结构的有限元方法（FEM），确保能量稳定性与一致的空间离散化。
• 采用隐式时间格式与伽辽金有限元格式，结合自适应时间步长，并通过网格尺寸与时间步长的收敛性分析进行验证。
• 将该框架应用于两类生物系统：大型卵细胞中的伤口愈合过程，以及受限细胞群落中的缺陷动力学。

实验结果

研究问题

• RQ1 如何为具有密度依赖效应的可压缩二维活性向列系统，构建一个统一且热力学一致的理论框架？
• RQ2 使用昂萨格变分形式理论，对推导活性向列系统的运动方程与数值格式有何影响？
• RQ3 压缩性与周转的引入如何影响活性系统中向列结构的自组织行为？
• RQ4 拓扑缺陷与活性流在大型卵细胞伤口愈合动力学中扮演何种角色？
• RQ5 活性强度与各向异性如何影响受限活性向列系统中的缺陷动力学与流场模式？

主要发现

• 该变分框架在成功恢复标准不可压缩活性向列模型的同时，将其推广至可压缩、密度依赖的系统。
• 有限元格式展现出最优的空间收敛性（收敛率 ∼h²）与时间收敛性（收敛率 ∼Δt），证实了数值方法的鲁棒性。
• 在伤口愈合模型中，活性向列流驱动缺陷迁移与重排，促进大型细胞中细胞骨架的再组织。
• 在受限群落中，增加活性导致活性湍流的出现，缺陷动力学由比例 ℓₐ/ℓ₀ 决定，其中 ℓₐ 为活性向列长度尺度。
• 各向异性活性（以 κ = λ_aniso/λ 表征）显著改变缺陷运动与流场模式，κ 值越高，越促进定向流动与缺陷成对形成。
• 该模型通过向列序、活性与流体动力学相互作用，成功捕捉了星形、 tactoids 与向列束等结构的自组织涌现。