Chirality is invisible to displacement data: a general symmetry-group criterion and taxonomy for the identifiability of active tissue deformations

技術要約：変位データにおいてカイラリティ（対掌性）は不可視である

問題提起
生物学的組織は、細胞プロセス（増殖、伸長、再配向）によって駆動される能動的な変形を起こすが、これらには回転またはカイラルな成分が含まれることが多い。モルフォエラスティシティ（形態弾性学）の文脈において、変形勾配 $F$ は能動テンソル $F_a$ と弾性テンソル $F_e$ に乗法的に分解される。ここで、ある重要な逆問題が生じる。すなわち、結果として得られる変位場 $u$ の観測から、能動テンソル $F_a$ を一意に特定できるか、という問いである。

先行研究（著者による同論文の随伴論文）では、二次元（2D）の識別問題における二分性が示された。等方性、対称な異方性、および対称なせん断の各クラスにおける $F_a$ は、ほぼマシン精度で識別可能であったが、歪対称（skew-symmetric）なクラスは、根本的かつ不可避な符号の曖昧性（すなわち、$F_a(s)$ と $F_a(-s)$ が同一の変位を生じさせること）を示した。しかし、これらの知見は、ケースごとの数値探索と特定の代数的一致に基づいたものであり、この曖昧性が特定のパラメータ化によるアーティファクトなのか、あるいは前方問題の一般的な構造的特性なのかという点は未解決のままであった。さらに、テストされていない $F_a$ のクラスに対して識別可能性を予測するための一般的な基準も存在していなかった。

手法
本論文は、識別問題を特定のパラメータ化の探索としてではなく、構成エネルギー密度 $W$ の対称群に関する問いとして再定式化する。

1. 前方演算子と等価性： 著者は、二つの能動テンソル $F_a$ と $F'_a$ が同一の変位場を与える場合、それらは観測的に等価であるという等価関係を定義する。これらすべての等価なテンソルの集合が、前方演算子の「ファイバー（fiber）」を形成する。
2. 対称群の解析： 本研究は、前方演算子のファイバーが、材料の保存エネルギー密度の材料対称群の下での $F_a$ の左軌道（left orbit）を含むことを特定している。具体的には、$Q$ が材料対称群の要素であり、$W(F_e Q) = W(F_e)$ を満たすとき、$F_a$ と $Q F_a$ は区別不能となる。
3. 構成則の導出： 圧縮性ネオ・フック型エネルギー（随伴研究で使用された構成則）について、著者は、材料対称群が正確に正回転群 $SO(d)$ であることを解析的に証明している。これは、能動テンソルの剛体的な事前回転が、変位データからは不可視であることを意味する。
4. 軌道交差基準： 著者は、クラス限定の識別可能性に関する一般的な基準を導出した。特定のパラメータ化されたクラス $\mathcal{C}$ 内の能動テンソル $F_a$ は、その $SO(d)$ 軌道が $F_a$ 自身以外にクラス $\mathcal{C}$ へ再進入しない場合に限り、識別可能である。
5. 検証： 理論的結果は以下の方法で検証された：
   • 記号代数チェック（対称群と回転恒等式の確認）。
   • ランダムなテンソルクラスを用いたモンテカルロ数値スキャン。
   • 回転対称性を破る境界条件下で $F_a$ と $Q F_a$ の前方問題を解き、同一の変位場がマシン精度で一致することを確認する直接的な非線形有限要素シミュレーション。

主な貢献

• 一般的な対称性基準： 本論文は、単一の、構成則に基づいた、幾何学に依存しない識別可能性の基準を確立した。識別可能性は、純粋に「クラスのメンバーの $SO(d)$ 軌道が、自身以外にそのクラスと交差するかどうか」によって決定される。
• 曖昧性の統一的説明： 基準は、歪対称クラスで見られた符号の曖昧性が、そのクラスが（再パラメータ化を除いて）$SO(2)$ 回転の下で閉じていることの帰結であることを説明する。逆に、対称なクラスは、テンソルが物理的に許容（正定値）である限り、その軌道が（単位元を除いて）クラス内に再進入しないため、識別可能であることを証明している。
• 予測能力： 基準は、これまでテストされていなかった二パラメータのオフダイアゴナル（非対角）クラスに対して適用された。これにより、非識別性はこの大きなクラス内における $p = -q$（歪対称のライン）に厳密に限定されていることが正しく予測され、数値スキャンによって確認された。これは、事前の知識なしに非識別的な領域を特定できる本ツールの能力を示している。
• 識別可能性の分類学： 本結果は、個別のパラメータ化ではなく、構成的対称性に基づき、散在する識別結果を統一された分類学へと整理している。

結果

• 理論的結果： 圧縮性ネオ・フック型材料において、材料対称群は正確に $SO(d)$ である。したがって、$F_a \sim Q F_a$ （ただし $Q \in SO(d)$）となる。
• 2D特有の詳細：
  • 識別可能なクラス： 等方性、一般的な対称な異方性、および対称なせん断のクラスは、正定性を維持したまま非自明な回転がクラスのメンバーをクラス内に写すことがないため、識別可能である。
  • 非識別なクラス： 歪対称クラスは識別不可能である。なぜなら、特定の回転 $Q(s) \in SO(2)$ がクラス内で一方を他方へと写すため、$F_a(s)$ と $F_a(-s)$ は観測的に等価となる。
  • 二パラメータ・クラス： 一般的なオフダイアゴナル・クラス $F_a(p,q)$ において、非識別性は $p+q=0$（歪対称のライン）のときにのみ発生する。それ以外のすべての点は識別可能である。
• 数値的結果： 有限要素シミュレーションにより、$F_a$ と $Q F_a$ に対して前方問題を解くと、境界条件が回転対称性を破る場合であっても、変位場がマシン精度（相対誤差 $10^{-14}$ から $10^{-15}$）で一致することが確認された。

意義と主張
本論文は、ケースごとの数値的な知見と、一般的な構造的理解との間の溝を埋めるものである。「曖昧性を探す」ことから「対称軌道の交差をチェックする」ことへと焦点を移すことで、著者は、識別可能性を a priori（事前）に判断するための信頼できる、検証可能なツールを提供している。

その意義は以下の点にある：

1. ケースごとの分析からの脱却： この基準により、研究者は徹底的な数値探索を行うことなく、新しいパラメータ化されたクラスに対して識別可能性を判断できる。
2. カイラリティの性質の解明： 変位データのみに依存する場合、カイラルな能動変形は、その大きさまでは識別できるが、その「手（向き）」については決して識別できないという、精密な力学的言明を提供している。
3. 既存文献への補完： 著者は、数値的な再構成を行う前に、まず解析的なステップが必要であると位置づけている。つまり、数値的な反転を試みる前に、構成的な非識別性（対称性によるもの）を排除すべきであると主張している。
4. 実用的な有用性： 未テストの二パラメータ・クラスにおける非識別的な箇所を予測することに成功したこと（当初の手計算による誤りを修正した結果）は、この基準が単なる抽象的な存在証明ではなく、堅牢で実用的なツールであることの具体的な証拠となっている。

著者は、結果が圧縮性ネオ・フック型エネルギーに特有であること、および本基準が非一意性のための「十分条件」（対称性によるもの）を提供するものであり、幾何学固有の他の非一意性の原因を排除するものではないことを述べ、範囲について謙虚な姿勢を保っている。また、2Dの場合は離散的な曖昧性が生じるが、3Dの場合（$SO(3)$）は、ジェネリックに連続的な非識別性につながる可能性があるとし、これを今後の課題としている。