Efficient Optimization of Tensor Rings with Low-Rank Environments

技術要約：低ランク環境を用いたテンソルリングの効率的な最適化

問題提起
テンソルリング（TR）分解は、周期的な量子系の自然な表現を提供し、周期境界条件（PBC）において、同等の精度を得るために必要なボンド次元（$\chi_r$）がテンソル・トレイン（TT）よりも大幅に小さくなることがよくあります。しかし、その閉じた幾何学的構造のために、TRの最適化は計算量的に困難です。開端境界を持つテンソル・トレインとは異なり、TRはグローバルな標準形（canonical form）を持ちません。そのため、ノルムおよびハミルトニアンの環境（environment）の操作は非自明となります。これらの環境を直接操作することは、次元が $\chi^2$ となる二重の仮想空間上の演算を伴い、結果として $O(\chi^6)$ の計算コストを招きます。さらに、標準的な特異値分解（SVD）はグローバルな計量（metric）を考慮しないため、局所的な切り捨てステップ（DMRGスイープに不可欠なもの）が非最適になります。既存の周期DMRG法（Ring-DMRG）は、環境を低ランク表現に圧縮することで、これらの問題を緩和し、コストを $O(p\chi^3)$ に削減しています（ここで $p$ は保持される環境ランク）。しかし、従来の実装には堅牢な2サイト更新スキームが欠けており、より不安定な反復ソルバーに依存しており、また、Ring-DMRGがテンソル・トレイン手法に対して真の優位性を持つレジームを明確に特定できていませんでした。

手法
著者らは、安定性、効率性、および切り捨て精度に対処するために、いくつかの技術的改善を統合した、強化されたRing-DMRGアルゴリズムを提案しています。

1. 低ランク環境の構築: この手法は、有効なノルムおよびハミルトニアン環境を近似するために、ランダム化特異値分解（rSVD）を利用します。フル環境行列を明示的に形成する代わりに（これは $O(\chi^6)$ のコストがかかります）、アルゴリズムはテンソルネットワーク形式で行列ベクトル積を計算し、支配的な $p$ 個の特異方向のみを保持します。これにより、環境の適用コストが $O(p\chi^3)$ に減少します。
2. 一般化Davidsonソルバー: 局所的な一般化固有値問題 $H^{(\ell)}x = \lambda N^{(\ell)}x$ を解くために、著者らは標準的なLanczos法に代わり、テンソルネイティブなDavidsonソルバーを採用しています。これにより、入れ子状のKrylov構造や、大きなボンド次元において悪条件となるノルム行列 $N^{(\ell)}$ の明示的な逆行列計算を回避できます。
3. ゲージ変換とフレーム変更: 局所的な固有値問題の条件数を改善するために、著者らは特定のゲージ変換を導入しています。彼らは、左側および右側の仮想境界にある簡約化されたノルム行列を用いてノルム環境を反復的にホワイトニング（whitening）することにより、「バランスド・フレーム（balanced frame）」を構築します。これは、従来の「ドミナント・フレーム（dominant frame）」アプローチと比較して、有効ノルム環境の条件数を大幅に減少させます。
4. 2サイト切り捨てスキーム: 中心的な貢献は、堅牢な2サイト更新です。2サイトテンソル $\Pi^{(\ell, \ell+1)}$ の切り捨ては、非自明な計量 $N^{(\ell, \ell+1)}$ によって定義されるフル状態のノルムにおける誤差を最小化する、重み付き低ランク近似問題として定式化されています。
   • 分離可能な計量: 環境が十分に分離可能である場合、フレームの変更により、標準的な切り捨てSVDによる切り捨てが可能になります。
   • フル計量による切り捨て: 計量が分離可能でない場合、交互最小二乗法（ALS）が採用されます。これは、環境を固定したまま、2つのテンソルを交互に最適化し、ガラーキン部分空間法と、悪条件を扱うためのデフレーション付き直交化を用いて線形最小二乗問題を解きます。
5. 信念伝播（Belief Propagation）との類似性: rSVDによる低ランク環境の構築は、信念伝播（BP）の一般化として解釈されます。標準的なBPは単一のメッセージ（ランク1）を保持しますが、rSVDアプローチは行列値メッセージ（ランク $r$）を伝播させ、環境の近似を系統的に改善します。

主要な結果

• 数値的安定性と速度: バランスド・フレームとDavidsonソルバーの組み合わせにより、収束性が大幅に向上しました。HeisenbergモデルおよびXXZモデルにおいて、バランスド・フレームは、ゲージなしの表現と比較して、ノルム環境の条件数を3桁減少させました。Davidsonソルバーは、$\chi=70$ において、Lanczos法と比較してスイープ時間を約5分の1に短縮しました。
• 切り捨て精度: 2サイト更新において、フル計量を用いたALSによる切り捨ては、複数の環境モードが大きく寄与するXXZモデルにおいて、最も小さな切り捨て誤差をもたらします。これは、標準的なSVD切り捨てやシングルフレーム近似と比較して、より正確な基底状態エネルギーと実時間ダイナミクスをもたらします。
• スケーリングと優位性のレジーム: 本論文では、必要な環境ランク $p$ がシステムサイズ $N$ およびボンド次元 $\chi$ に対してどのようにスケールするかを分析しています。
  • ギャップのある系: ギャップのある相では、環境ランク $p$ は小さく $N$ に依存せず、テンソルリングのボンド次元 $\chi_r$ はテンソル・トレインの等価な $\chi_t$ よりも大幅に小さくなります（$\chi_t \propto \chi_r^{1.5}$ とスケール）。ここで、Ring-DMRGは実質的な計算上の優位性を提供します。
  • 臨界系: 熱力学的極限に近づく臨界系について、著者らは「有限サイズ–有限エンタングルメント」レジームを特定しています。周期的なグローバル幾何学への感度を維持するためには、有限エンタングルメント相関長 $\xi_\chi \sim N$ となるようにボンド次元をスケールさせる必要があります。この軌跡に沿って、保持される環境ランク $p$ は一定に保たれます。その結果、Ring-DMRGは、ボンド次元に対して $O(\chi^3)$ のスケーリングを維持し（ただし $\chi$ はより小さい）、周期的なグローバル特性に敏感な観測量（例：長距離相関、中心電荷の推定）に対して、テンソル・トレインよりも計算量的に好ましいものとなります。
• 中心電荷の推定: スピン1/2 Heisenberg鎖の中心電荷を推定するテストケースにおいて、テンソルリング表現は $\chi_{TR}=150$ で $c_{eff} \approx 0.993$ の精度を達成しましたが、テンソル・トレインでは $\chi_{TT} \approx 1500$ を必要としました。これは、このレジームにおいて、リング表現が8倍程度のスピードアップの可能性を持つことを示唆しています。

意義と主張
本論文は、周期的なテンソルネットワーク最適化における歴史的なボトルネックを克服する、堅牢かつ効率的なRing-DMRGの定式化を提供することを主張しています。環境を考慮した切り捨てを伴う2サイト更新の導入と、ゲージ変換によるソルバーの安定化により、著者らは、Ring-DMRGが数値的に安定しているだけでなく、特定のクラスの問題に対してテンソル・トレイン手法よりも計算量的に優れていることを示しています。

その意義は、特に短い相関長を持つ周期境界条件や、有限サイズレジームにおける臨界系を含む、テンソルリングのアンザッツが真の優位性を持つ具体的なレジームを特定したことにあります。本研究は、テンソルネットワークの最適化と信念伝播の間の架け橋となり、有限次元の環境モードの部分空間を保持すること（rSVD経由）が、標準的なBPの系統的な一般化であることを示唆しています。著者らは、現在の実装が（壁時計時間で競合するために）ブロック疎演算などのさらなる最適化を必要とするものの、理論的枠組みは、Ring-DMRGが周期的な量子多体系に対する実行可能かつしばしばより優れた代替手段であることを確立していると結論付けています。