Quantum computing-based solver for interacting power grids

技術要約：相互作用する電力網のための量子コンピューティング・ベースのソルバー

問題提起
マルチターミナル送電網におけるパワーエレクトロニクスおよび再生可能エネルギー源の統合が進むにつれ、複雑な動的不安定性と高調波歪みが生じている。共振モード解析（RMA）は、システムのアドミタンス行列を分解することで臨界共振モードを特定するための標準的な手法であるが、大規模ネットワークに適用する場合、深刻な計算上のボトルネックに直面する。これらのシステムにおけるアドミタンス行列は、能動成分と無効成分の両方の存在により非エルミート（non-Hermitian）であり、かつ高次元である。非エルミート・ランチョス法のような古典的な反復手法は、膨大なクリロフ部分空間ベクトルを格納する必要があり、$n \times n$ 行列に対して $O(n^3)$ でスケールする対角化を行う。グリッドの規模が数万バスへと拡大するにつれ、これらのメモリおよび処理要件は、古典的なハイパフォーマンス・コンピューティング（HPC）の限界を超える。さらに、変分量子固有値ソルバー（VQE）のような標準的な量子固有値ソルバーは、エルミート演算子用に設計されており、電力系統の安定性解析に内在する非エルミート固有値問題を直接解くことはできない。

手法
著者らは、古典的なメモリ制約と標準的な量子アルゴリズムの両方の限界を克服するために、実分散ベースの変分量子固有値ソルバー（RVVQE）を利用した、新しい量子・古典ハイブリッド・フレームワークを提案している。その手法は以下の通りである：

1. システム定式化： グリッドは周波数スキャン解析を用いてモデル化され、節点電圧と注入電流の関係はアドミタンス行列 $Y_{B,f}$ によって定義される。共振は、$Y_{B,f}$ の固有値がゼロに近づき、モーダルインピーダンスがピークに達したときに特定される。
2. 行列分解： 非エルミートであるアドミタンス行列 $Y_{B,f}$ を量子ハードウェアと互換性を持たせるため、数学的に2つの純粋なエルミート行列 $H_f$ と $K_f$ に分解する。すなわち、$Y_{B,f} = H_f + iK_f$ である。
3. 量子マッピング： これらのエルミート成分は、$N$ 量子ビットの量子プロセッシングユニット（QPU）での実行のためにパウリ文字列へとマッピングされる。ここで $N = \lceil \log_2 n \rceil$ である。この対数的なスケーリングにより、数千のバスを持つシステムを、比較的少ない数の量子ビット（例：70,000バスに対して17量子ビット）で表現することが可能となる。
4. RVVQE アルゴリズム： パラメータ化された量子回路（アンザッツ）が試行状態 $|\psi(\theta)\rangle$ を準備する。エネルギーを最小化する代わりに、アルゴリズムは以下の実数値分散コスト関数を最小化する：
   $$C_{var}(\theta) = \langle H_f^2 \rangle_\theta - \langle H_f \rangle_\theta^2 + \langle K_f^2 \rangle_\theta - \langle K_f \rangle_\theta^2$$
   古典最適化器は、このコスト関数がゼロに向かうようにパラメータ $\theta$ を反復的に更新する。
5. 固有値抽出： 収束後、元の非エルミートシステムの複素固有値は、エルミート成分の期待値の和によって再構成される：$\lambda = \langle H_f \rangle + i\langle K_f \rangle$。その後、臨界共振モーダルインピーダンス $Z_{cm,f} = 1/\lambda$ が計算される。
6. デフレーション： 複数の固有値を抽出するために、直交性を確保するために、以前に見出された固有状態とのオーバーラップに基づいたペナルティをコスト関数に加える逐次デフレーション技術が用いられる。

主な貢献

• アルゴリズムの適応： 本論文は、非エルミート演算子用に設計された RVVQE アルゴリズムを、電力系統の共振解析という特定の文脈に適応させ、散逸的な電気グリッドと開放量子系の間の溝を埋めている。
• スケーラブルなアーキテクチャ： グリッドの状態を量子メモリに対数的にエンコードすることで、提案手法は古典的なRAMの $O(n^2)$ メモリ制限を回避し、大陸規模のネットワークを分析するための理論的にスケーラブルな経路を提供する。
• ハイブリッド・ワークフロー： 周波数スイープ、行列分解、および安定性評価をハイブリッド量子・古典ループを用いて処理する、完全な「グリッド・トゥ・量子」ワークフロー（アルゴリズム1）を定式化した。

結果
本フレームワークは、標準的な5バス送電システムを用いて検証された。非エルミート・アドミタンス行列は、100 Hz から 3000 Hz の周波数範囲にわたって生成された。5x5 行列は、3量子ビットの量子システムに適合させるために 8x8 にパディングされた。

• 精度： RVVQE アルゴリズムによって計算された臨界共振モーダルインピーダンス（$|Z_{cm}|$）は、全周波数スペクトルにわたって、標準的な行列対角化（SciPy ライブラリ）から得られた厳密な古典解とほぼ完全に一致した。
• ハードウェア制約： シミュレーションは量子誤差緩和（QEM）なしで行われた。誤差訂正がないにもかかわらず、結果は古典的な参照値に近い値を維持しており、現在のノイズの多いハードウェアに対するアルゴリズムの堅牢性を示唆している。
• 安定性解析： 本手法は共振周波数と臨界モードを正常に特定し、小信号安定性と参加係数を評価する能力を実証した。

意義と主張
本論文は、この研究が、現在古典的な計算境界を超えている大規模なネットワークにおける共振不安定性の診断が可能な、「極めてスケーラブルなアーキテクチャ」を提供すると主張している。著者らは、現在のハードウェアには制限があるものの、量子ビットの対数的なスケーリングが、将来の大規模な電力系統解析への実行可能な道筋を提供することを強調している。非エルミート行列に対する RVVQE の適用成功は、複雑な電力系統の安定性研究（将来のより大きな IEEE テストネットワーク（例：118バス、300バス）や合成グリッド ACTIVSg70k への応用を含む）に量子コンピューティングを適用するための基礎的な一歩として提示されている。本研究は、産業規模の問題を実際のハードウェア上で解決したと主張するものではなく、ベンチマークシステムにおけるアプローチの数学的妥当性と精度を実証したものである。