$\mathrm M$-Eigenpairs of Partially Symmetric Tensors: Exact Reformulation and Perturbation Bounds

技術要約：部分対称テンソルのM-固有対：厳密な再定式化と摂動境界

問題提起
本論文は、弾性理論において自然に発生する弾性テンソルである、4次の部分対称テンソルにおけるすべての実M-固有対を見つけるという計算上の課題に取り組んでいる。テンソル $A \in \mathbb{R}^{m \times n \times m \times n}$ のM-固有対 $(\lambda, x, y)$ は、独立した球面制約（$x^\top x = 1, y^\top y = 1$）に従う2つのベクトル変数を含む結合非線形系によって定義される。既存の文献は、包含境界、強弾性性に関する十分条件、および極値（最大または最小）M-固有値のための数値解法に焦点を当ててきたが、M-固有対の完全な集合を計算することは依然として困難である。標準的なZ-固有値問題へと変換する既存のリフティング戦略は、特定のユークリッド正規化に依存しており、より一般的な正規化フレームワークが厳密なスペクトル対応を保持できるかという問いが残されていた。

手法
著者らは、M-固有値問題を一般化されたテンソル固有値問題として再定式化する、パラメータ化されたリフティング・フレームワークを提案している。その手法は以下のステップで進む：

1. テンソル・リフティング： 4次の部分対称テンソル $A$ を、関連するリフトされたテンソル $B_A$ を構築することにより、高次元空間 $\mathbb{R}^{(m+n) \times (m+n) \times (m+n) \times (m+n)}$ に埋め込む。この構成は、結合されたM-固有値問題の式を、ベクトル $w = (x^\top, y^\top)^\top$ を用いた単一のモード1テンソル方程式へと写像する。
2. 対称化： 構築されたテンソル $B_A$ を対称化して $\bar{B}_A$ を作成し、すべてのインデックス置換を満たすようにする。
3. パラメータ化された正規化： 重み付き正規化行列 $B_{\alpha,\beta} = \text{diag}(2\beta^2 I_m, 2\alpha^2 I_n)$ を導入する（ここで $\alpha, \beta > 0$）。これは以下の一般化固有値問題を定義する：
   $$\bar{B}_A w^3 = \lambda B_{\alpha,\beta} w, \quad w^\top B_{\alpha,\beta} w = 1$$
   この定式化は、標準的なZ-固有値問題（$\alpha = \beta = 1/\sqrt{2}$ の場合に回収される）を一般化している。
4. 厳密な対応と回収： 本論文は、元のテンソル $A$ のM-固有対と、リフトされたテンソル $\bar{B}_A$ の$B_{\alpha,\beta}$-固有対との間の厳密かつ正確な対応関係を確立している。リフトされた空間の解を元のM-固有対へと写すための具体的な回収公式が導出されており、ゼロ固有値の扱いについても特定されている。
5. 摂動解析： 加法的摂動下での最大M-固有値に対する摂動境界を導出する。著者らは、リフティング変換が（スケーリング因子を除いて）これらの境界を正確に保持し、追加の緩和を導入しないことを示している。

主な貢献
本論文の主な貢献は以下の3点である：

• パラメータ化された再定式化： 標準的なZ-固有値再定式化を、より一般的な重み付きフレームワークへと拡張する、パラメータ化された$B_{\alpha,\beta}$-正規化を伴うリフトされたテンソル構成を提案した。
• 厳密な対応： 元のテンソルのM-固有対と、リフトされたテンソルの$B_{\alpha,\beta}$-固有対との間の一対一の対応関係を確立し、すべての実M-固有対を回収するための明示的な手順を提供した。
• 摂動境界の保持： 元の問題とリフトされた形式の両方において最大M-固有値の摂動境界を導出し、リフティング変換がスペクトル情報と摂動推定値を、追加の緩和なしに正確に保持することを証明した。

数値結果
提案手法の有効性は、MATLABで実装された5つの数値実験を通じて検証されている：

• 全スペクトルの回収： 例4.1および4.2（$2\times2$ テンソルおよびドロマイトの$3\times3$ 弾性係数テンソルを含む）において、アルゴリズムは既知の参照結果と一致する、すべての異なる実M-固有値（それぞれ6個および22個）を正常に回収した。
• パラメータ感度： 例4.3は、さまざまなパラメータ $(\alpha, \beta)$ の選択が、高い数値精度（残差 $\approx 10^{-15}$）で同一のM-固有値集合をもたらすことを示した。一部の非標準的なパラメータ選択では、標準的なZ-固有値の場合と比較して計算時間がわずかに短縮された。
• スケーラビリティ： 例4.4では、様々な次元のランダムテンソルに対して、パラメータ化された $B_{\alpha,\beta}$ 定式化が標準的な定式化と同じ数の固有対を一貫して計算しつつ、より低い平均CPU時間を達成した（加速比は1.08から3.01の範囲）。
• 摂動境界： 例4.5は、計算された摂動後の最大M-固有値が、常に理論的に導出された下限および上限の範囲内に存在することを確認し、摂動レベルが減少するにつれて区間の幅が急速に減少することを示した。

意義と主張
本論文は、提案するフレームワークが、極値固有値の計算よりも注目が少ないタスクである「完全な」実M-固有対を計算するための、堅牢で柔軟な手法を提供すると主張している。パラメータ化された正規化を導入することで、著者らは、標準的なZ-固有値アプローチが彼らのより広範なフレームワークの特殊なケースに過ぎないことを示している。極めて重要な点として、本論文は、この一般化が理論的保証を損なわないことを強調している。具体的には、リフティングプロセスは元の問題の厳密な摂動境界を保持する。数値実験は、この柔軟性を活用することで、精度を犠牲にすることなく計算効率を向上させられる可能性を示唆している。