Quantum-Accelerated Nonlinear Seismic Simulation for Megathrust-Prone Urban Environments: A Three-Order-of-Magnitude Speedup via HHL-Embedded Newton–Raphson Integration

技術要約：量子加速による非線形地震シミュレーション

問題提起
非線形時刻歴解析（NLTHA）は、構造物の耐震性能を評価するための厳格な標準であり、ASCE 7-22やペルーのNTE E.030などの規制によって、不整形構造物や重要構造物に対して義務付けられている。しかし、その適用は計算コストによって著しく制限されている。$N$自由度（DOF）を持つシステムの場合、古典的な解法は各ニュートン・ラプソン反復において線形方程式系を解く必要があり、結果として1タイムステップあたりの時間計算量は$O(N^3)$となる。このスケーリングにより、NLTHAは単一の建物研究に限定されており、アンデス沈み込み帯のメガスラストの影響を受ける都市環境のような高ハザード地域における、リアルタイムかつ都市規模のリスク評価を阻害している。

手法
本研究では、計算時間を3桁削減するために設計された、量子強化ニュートン・ラプソン（QNR）統合フレームワークを提案する。この手法は、ハロー・ハッシディム・ロイド（HHL）量子線形ソルバーを、古典的なヒルバー・ヒューズ・テイラー（HHT-$\alpha$）タイムステップ法に統合するものである。

1. 古典的定式化: $N$-DOF非線形システムの運動方程式をHHT-$\alpha$を用いて離散化する。各タイムステップにおいて、非線形残差はニュートン・ラプソン反復によって解決され、そこでは$K_T^{(i)} \Delta u^{(i)} = R^{(i)}$（ここで$K_T$は有効接線剛性行列）を解くことが要求される。
2. 量子マッピング:
   • 状態エンコーディング: 残差力ベクトル$R$を、$n = \lceil \log_2 N \rceil$個の量子ビットを用いて振幅エンコーディングし、量子状態 $|b\rangle$ に変換する。
   • ハミルトニアン分解: 正規化された接線剛性行列 $A = K_T / \|K_T\|$ を、量子ハミルトニアンを形成するためにパウリ演算子のテンソル積の和へと分解する。
   • HHLの実行: HHLアルゴリズムは、量子位相推定と制御回転を実行し、状態 $|x\rangle \propto A^{-1}|b\rangle$ を生成することで、実質的に線形システムを解く。
   • 読み出し: 変位増分は、古典シャドウ・トモグラフィーを介して抽出され、疎な観測量に対しては$O(\log N)$の測定のみを必要とする。
3. 検証プロトコル: 本フレームワークは、64-DOFの非線形せん断型建物モデルを用いた1,000回のモンテカルロ研究を通じて検証された。地震動は、リマ（土壌型S3）の設計スペクトルに一致するようにPEER NGA-West2データベースからスケーリングされた。統計的有意性は、二標本t検定およびシャピロ・ウィルク正規性検定を用いて評価された。

主要な貢献

• QNRフレームワーク: 構造力学におけるHHT-$\alpha$ニュートン・ラプソン・ループ内にHHLを組み込んだ、数学的に厳密な定式化。
• 物理的定式化: SDOFおよびMDOFの物理原則を量子状態へマッピングすることを定式化するための、自由体図（FBD）の使用。
• 定量的ベンチマーク: 古典的な$O(N^3)$の時間計算量および$O(N^2)$の空間計算量から、量子的な$O(s^2\kappa^2 \log N)$の時間計算量および$O(\log N)$の空間計算量への移行を示すマルチメトリック比較行列。
• 統計的堅牢性: 平均加速率1,002.5倍（95%信頼区間 [980.2, 1,024.8]）および、フル精度古典リファレンスに対する平均二乗誤差（MSE）$1.2 \times 10^{-4}$を示す包括的な検証。

結果
QNRフレームワークは、古典的なHHT-$\alpha$ベースラインと比較して、平均加速係数（$\mu$）1,002.5倍を達成した。

• 統計的有意性: 結果は統計的に有意であり（$t = 2,814.5, p < 0.0001$）、95%信頼区間は「加速なし」という帰無仮説を除外している。加速率の分布は、シャピロ・ウィルク検定（$W = 0.998, p = 0.31$）により正規分布であることが確認された。
• 精度: MSE $1.2 \times 10^{-4}$ は、ドリフトに基づく構造設計の1%の工学的許容閾値を十分に下回っている。
• スケーラビリティ: 理論上、本フレームワークは、古典的な実用限界である$N \approx 10^3$に対し、$N \approx 10^6$ DOFまでのモデルをサポートする。
• 効率性: 本研究は、1解あたりの相対的なエネルギー削減が約100倍であることを推定している。

意義と主張
本論文は、このフレームワークがNLTHAを、単一の建物に対する計算負荷の高いツールから、リアルタイムの都市規模地震リスク評価のための実行可能なプラットフォームへと移行させるものであると主張している。具体的には、ナスカ・南米沈み込み帯に位置する軟弱地盤盆地上にある大都市リマにおいて、QNRフレームワークは数千の構造物の同時非線形評価を可能にする。この能力は、NTE E.030およびCENEPREDの災害リスク軽減命令への準拠を直接支援するものである。

著者らは、ペルーの文脈における3つの変革的な応用を特定している：

1. リアルタイム都市規模シミュレーション: 約10,000個の鉄筋コンクリートフレームを同時にシミュレートすることにより、緊急対応（例：CENEPRED）のための即時損傷状態マップを生成する。
2. 確率論的地震ハザード解析（PSHA）: リマの地質に特有の地盤・構造物相互作用を組み込んだ非線形脆弱性曲線を生成するために、数百万回の実行を行う。
3. 地震早期警戒システムの統合: 量子NLTHAをSWAVEシステムと結合し、破壊的なS波が到着する前に、重要インフラ（病院、橋、ダム）の非線形応答を予測する。

限界と今後の課題
著者らは、現在のNISQ（Noisy Intermediate-Scale Quantum）デバイスのゲートエラー率（$\sim 10^{-3}$）が回路深度を制限するため、現在の結果がフォールトトレラント・プロセッサの理論的なゲート数モデルに基づいていることを認めている。さらに、行列をロードするための状態準備のオーバーヘッドは、最悪のシナリオにおいて$O(N^2)$でスケールする可能性があり、メッシュ細分化に伴う条件数の増大は、プリコンディショニングなしでは量子優位性を損なう可能性がある。今後の課題として、近未来のハードウェア向けの変分量子線形ソルバー（VQLS）の実装、既存の量子プロセッサ（IBM Quantum, IonQ）でのベンチマーク、およびリマの地質プロファイルに校正された非線形地盤モデルの組み込みが提案されている。