[論文レビュー] Numerical Methods for Biomembranes: conforming subdivision methods versus non-conforming PL methods

Jingmin Chen, Thomas Yu+4 more|arXiv (Cornell University)|Jan 28, 2019

Advanced Numerical Analysis Techniques参考文献 64被引用数 4

ひとこと要約

本稿は、バイオ膜モデリングにおける幾何的変分問題、特にウィルモアエネルギー最小化の文脈で、適合型サブディビジョンサーフェス（SS）法と非適合型ピecewise線形（PL）法を比較している。PL法に基づく離散ウィルモアエネルギーの最小化は、滑らかでないウィルモア最小化子への収束に失敗することを証明している一方、SS法による最小化は、高次精度と適合性のおかげで成功する。これにより、コンformal幾何学を用いたPL法の正則化が提案された。

ABSTRACT

The Canham-Helfrich-Evans models of biomembranes consist of a family of geometric constrained variational problems. In this article, we compare two classes of numerical methods for these variational problems based on piecewise linear (PL) and subdivision surfaces (SS). Since SS methods are based on spline approximation and can be viewed as higher order versions of PL methods, one may expect that the only difference between the two methods is in the accuracy order. In this paper, we prove that a numerical method based on minimizing any one of the `PL Willmore energies' proposed in the literature would fail to converge to a solution of the continuous problem, whereas a method based on minimization of the bona fide Willmore energy, well-defined for SS but not PL surfaces, succeeds. Motivated by this analysis, we propose also a regularization method for the PL method based on techniques from conformal geometry. We address a number of implementation issues crucial for the efficiency of our solver. A software package called Wmincon accompanies this article, provides parallel implementations of all the relevant geometric functionals. When combined with a standard constrained optimization solver, the geometric variational problems can then be solved numerically. To this end, we realize that some of the available optimization algorithms/solvers are capable of preserving symmetry, while others manage to break symmetry; we explore the consequences of this observation.

研究の動機と目的

• バイオ膜の幾何的変分問題を解く際のピースワイズ線形（PL）法とサブディビジョンサーフェス（SS）法の収束特性を分析すること。
• 理論的に一貫性があるにもかかわらず、PL法に基づくウィルモアエネルギー最小化が滑らかな解に収束しない理由を特定すること。
• コンformal幾何学を用いた正則化技術をPL法に提案し、収束特性を改善すること。
• PLおよびSSサーフェス上で幾何的関数を計算するための並列ソフトウェアパッケージ Wmincon を開発・公開すること。
• ウィルモア型問題に適用した最適化ソルバーにおける対称性の保存と破れを調査すること。

提案手法

• ウィルモアエネルギーを連続的な幾何的関数として定式化し、PLおよびSSサーフェス用の離散類似物を定義する。
• 一貫性のあるPLウィルモアエネルギー最小化が、適合性の欠如と不十分な正則性のため、滑らかなウィルモア最小化子に収束しないことを証明する。
• 離散的コンフォーマル幾何学に基づくコンフォーマルパラメトリゼーションを用いたPL法の正則化を導入し、収束を安定化する。
• Wmincon と呼ばれる並列 C++ ソフトウェアパッケージに、ウィルモア、面積、体積、平均曲率などの幾何的関数を実装する。
• Wmincon を用いて標準的な制約付き最適化ソルバーを用いて、Canham-Helfrich-Evans 問題を数値的に解く。
• 最適化ソルバーの対称性の挙動を分析し、対称性を保存するか、対称性を破するかのアルゴリズムを区別する。

実験結果

リサーチクエスチョン

• RQ1 理論的に一貫性があるにもかかわらず、なぜPL法に基づく離散ウィルモアエネルギー最小化は滑らかなウィルモア最小化子に収束しないのか？
• RQ2 コンフォーマル幾何学に基づく正則化によって、PL法におけるウィルモアエネルギー最小化の収束性を回復させることができるか？
• RQ3 異なる最適化ソルバーは、幾何的変分問題における解の対称性にどのように影響を与えるか？
• RQ4 適合性と高次近似が滑らかな解への収束を保証するために果たす役割は何か？
• RQ5 提案された正則化技術は、幾何PDEにおける他の非適合有限要素法へ一般化可能か？

主な発見

• 任意の一貫性のあるPLサーフェス上での離散ウィルモアエネルギーは、命題 3.9 および 3.12 で示されるように、滑らかなウィルモア最小化子に収束しない。
• サブディビジョンサーフェス法は、高次近似と適合性のおかげで、PL法とは異なり、滑らかなウィルモア最小化子に成功して収束する。
• コンフォーマル幾何学に基づくPL法の正則化により、収束性が回復され、数値的安定性が向上する。
• ソフトウェアパッケージ Wmincon を用いることで、PLおよびSSサーフェス上で幾何的関数を効率的に並列計算でき、変分問題の数値的解法が促進される。
• 一部の最適化ソルバーは対称性を保存し、対称な最小化子に収束するが、他のソルバーは対称性を破り、非対称な解に収束するため、数値結果の信頼性に影響を与える。
• 理論的分析により、PLサーフェスの空間は、 genus 0 であっても、滑らかなウィルモア最小化子への一貫性のある収束をサポートするにはあまりに大きく不規則であることが示唆されている。