[論文レビュー] On the generalization of the Wigner semicircle law to real symmetric tensors

Razvan Gurău|arXiv (Cornell University)|Apr 6, 2020

Tensor decomposition and applications参考文献 24被引用数 9

ひとこと要約

本稿は、$ p \geq 3 $ の実対称テンソルに対して、Wignerの半円則を一般化する。テンソルレゾルベントを導入し、テンソル不変量の積分表現を導出する。$ N \to \infty $ の極限において、期待レゾルベントは有限なスペクトルカットを示し、その不連続性は任意のテンソル次数に一般化された半円則に従う普遍的法則に従う。スピikedテンソルモデルでは、信号対雑音比が臨界値に達する際に鋭い相転移が観測される。

ABSTRACT

We propose a simple generalization of the matrix resolvent to a resolvent for real symmetric tensors $T\in \otimes^p \mathbb{R}^N$ of order $p\ge 3$. The tensor resolvent yields an integral representation for a class of tensor invariants and its singular locus can be understood in terms of the real eigenvalues of tensors. We then consider a random Gaussian (real symmetric) tensor. We show that in the large $N$ limit the expected resolvent has a finite cut in the complex plane and that the associated "spectral density", that is the discontinuity at the cut, obeys a universal law which generalizes the Wigner semicircle law to arbitrary order. Finally, we consider a spiked tensor for $p\ge 3$, that is the sum of a fixed tensor $b\,v^{\otimes p}$ with $v\in \mathbb{R}^N$ (the signal) and a random Gaussian tensor $T$ (the noise). We show that in the large $N$ limit the expected resolvent undergoes a sharp transition at some threshold value of the signal to noise ratio $b$ which we compute analytically.

研究の動機と目的

• 実対称テンソルの次数 $ p \geq 3 $ に対して、ランダム行列理論で用いられる行列レゾルベントを拡張した一般化されたレゾルベントを定義すること。
• テンソルレゾルベントの特異点とスペクトル密度の間の関係を確立し、高次テンソルへの固有値分布の概念を一般化すること。
• ガウス確率テンソルモデルの $ N \to \infty $ 極限を解析し、レゾルベントの不連続性に対する普遍的スペクトル則を導出すること。
• 固定された信号テンソルがガウス確率テンソルに加えられたスパikedテンソルモデルを研究し、信号検出のための臨界閾値を同定すること。

提案手法

• 行列レゾルベント $ N^{-1} \operatorname{Tr}[(w - T)^{-1}] $ を高次テンソルに一般化した、対称的 $ p $-線形形式としてのテンソルレゾルベント $ \omega(w;T) $ を導入する。
• 曲線 $ \gamma $ における $ w $ の contour 積分を用いて、テンソル不変量 $ I(T) $ の積分表現を、表現関数 $ h_I(w) $ を用いて導出する。
• 分散 $ \frac{N^{p-1}}{2p} \sum (T_{a_1\dots a_p})^2 $ を持つガウス確率テンソルに対して、$ N \to \infty $ 極限におけるレゾルベントを解析し、実軸上に有限なカットへの収束を示す。
• カットを越える際のレゾルベントの不連続性としてスペクトル密度を計算し、一般の $ p \geq 3 $ に対してMeijer G関数を用いて表現する。
• 信号テンソル $ A = \frac{b}{N^{p/2 - 1}} v^{\otimes p} + T $ を含むスパikedテンソルモデルを研究し、$ v $ を単位ベクトル、$ T $ をガウス分布に従うものとして、レゾルベントの特異点を用いて信号の存在を検出する。
• 最急降下法とインスタントン解析を用いて、分配関数の不連続性を計算し、レゾルベントの解析的構造における相転移を同定する。

実験結果

リサーチクエスチョン

• RQ1 行列レゾルベントは、次数 $ p \geq 3 $ の実対称テンソルに対してどのように一般化可能か。また、この一般化されたレゾルベントの解析的構造はいかなるものか。
• RQ2 $ p \geq 3 $ のガウス確率テンソルの $ N \to \infty $ 極限における、レゾルベントの不連続性に対する普遍的スペクトル則は何か。
• RQ3 ランク1の信号を持つスパikedテンソルモデルでは、信号対雑音比 $ b $ の増加に伴い、レゾルベントの特異点集合に鋭い相転移が現れるか。
• RQ4 次数 $ p \geq 3 $ の場合に、$ N \to \infty $ 極限において信号検出の臨界閾値 $ b_t $ を解析的に計算可能か。
• RQ5 テンソルレゾルベントのスペクトル密度は、高次テンソルに対する何らかの一般化された固有値の概念と関係があるか。

主な発見

• ガウス確率テンソルの次数 $ p \geq 3 $ に対して、$ N \to \infty $ 極限における期待レゾルベントは、実軸上に有限なカット $ \left[-\sqrt{p^p / (p-1)^{p-1}}, \sqrt{p^p / (p-1)^{p-1}} \right] $ を示す。
• スペクトル密度（カットを越える際のレゾルベントの不連続性）は、任意の $ p \geq 3 $ に対してWigner半円則を一般化する普遍的法則に従い、Meijer G関数を用いて表現可能である。
• $ p=2 $ の場合、一般化された法則は標準的なWigner半円則に還元されるが、$ p \geq 3 $ の場合、法則は特異的かつ非ガウス的である。
• スパikedテンソルモデルにおいて、信号対雑音比の臨界値 $ b_t^2 = \frac{(p-1)^p}{(p-2)^{p-2}} $ で鋭い転移が発生し、レゾルベントの最大で除去できない特異点が $ \sqrt{p^p / (p-1)^{p-1}} $ から $ p^{p/2} $ へとジャンプする。
• この転移は、テンソル内にランク1の信号が存在するかどうかの検出閾値として解釈され、レゾルベントの解析的構造が信号優位性の始まりを示唆する。
• レゾルベントの不連続性はインスタントン寄与によって計算可能であり、$ \mathcal{Z}_+(y) - \mathcal{Z}_-(y) = \sum_\nu \frac{1}{\sqrt{\mathcal{S}''(\phi_\nu)}} e^{-\mathcal{S}(\phi_\nu)} $ という式で与えられる。ここで $ \phi_\nu $ は $ \mathcal{S}'(\phi_\nu) = 0 $ を満たす実インスタントンである。