パルス最適化
パルス最適化は、ゲートの不忠実度 C = 1 − F に対する勾配降下法によって制御波形を数値的に彫刻することです — GRAPEの解析的勾配(2005)から、勾配がノイズそのものを貫いて流れる現代の微分可能な開放系オプティマイザまで。
どういう意味か
ゲートを制御包絡線ε(t)として理解すれば、最良のゲートを見つけることは最適化問題になります: 実現された発展と目標ユニタリの間の忠実度Fに対し、コストC = 1 − Fを最小化するのです。勾配降下法はパルスパラメータθを勾配の逆方向に更新します: θ_{k+1} = θ_k − η∇_θC(θ_k)。GRAPE(Khanejaら、2005)がこの問題に解析的勾配を導入し、自動微分(autograd)がそれを一般化します — 時間発展全体(例: トロッター分割されたリンドブラッド伝播器)をPyTorchの演算で構成すれば、1回の逆伝播ですべてのパラメータの正確な勾配が同時に得られます。有限差分ならパラメータごとに1回のシミュレーションが必要です。現代のパルスは数百〜数千のパラメータを持つため、これが決定的です。最前線は開放系最適化です: ノイズのないシュレーディンガー・シミュレーションで最適化してから点数を事後割引(e^{−T/T1}を掛ける)する代わりに、GradPulseのような微分可能オプティマイザはリンドブラッド・マスター方程式を通じてT1緩和・純位相緩和(Tφ)・漏れをforward passの中に入れます — オプティマイザは本物のノイズ世界のスコアボードの上を歩くのです。代償は大きい — 密度行列はd×dなので、27次元の明示的カプラモデルは729個の要素を持ちます — そして伝播器を書き直すことは、独立ソルバー(正確な行列指数、適応型ODEのQuTiP)との約10⁻¹⁴レベルの相互検証を要求し、それはCIで強制されます。日常のたとえ
よくある誤解
- 勾配なし(gradient-free)の最適化はスケールしません — パルスパラメータが数百〜数千個になると、Nelder-Meadのような手法は次元とともに急激に遅くなりますが、逆伝播は巻き戻し1回ですべての勾配を与えます。微分可能性は好みではなく、規模における生存条件です。
- 2つのシミュレータが一致しても、コードを共有していたり同じ近似系列を使っていたりすれば何も証明しません — 共有されたトロッター型の時間離散化は、両者の比較では見えない共通モード誤差を生みます。独立した再実装(正確な行列指数 vs トロッター分割 vs 適応型ODE)だけが、一致を証拠にします。
要点
- GRAPE(2005)がパルス設計に解析的勾配を導入し、autogradがそれを任意の微分可能シミュレーションへ一般化します。
- 逆伝播は1回のbackward passですべてのパルスパラメータの勾配を与えます — 有限差分ならパラメータごとに1回のシミュレーションが必要です。
- 開放系のコスト: 密度行列はd×dなので、27次元の明示的カプラモデルは729個の要素を持ちます — 正直なノイズモデリングは高くつきます。
- 微分可能な開放系オプティマイザは、事後のe^{−T/T1}割引の代わりにT1・Tφ・漏れをforward passの中に入れます — 勾配がノイズを貫いて流れます。
- 三重ソルバー相互検証(PyTorchトロッター vs 正確なNumPy行列指数 vs QuTiP適応型ODE)の約10⁻¹⁴の一致が転写バグを捕らえ、どれか2つのソルバーがずれればCIゲートがビルドを失敗させます。
理解度チェック
規模が大きくなるとき、微分可能性が(単なる好みではなく)生存条件である理由は何ですか?
- A.アナログハードウェアが微分可能な波形しか受け付けないから
- B.パラメータが数百〜数千個あるとき、逆伝播の1回のbackward passだけが勾配計算を現実的にするから
- C.勾配なしの手法ではガウス包絡線を表現できないから
- D.微分可能なコードはGPUで動き、GPUが必須だから
答えを見る
答え: B. パラメータが数百〜数千個あるとき、逆伝播の1回のbackward passだけが勾配計算を現実的にするから
理由: 有限差分はパラメータごとに1回のシミュレーションが必要で、勾配なしの探索は次元とともに急激に遅くなります。逆伝播は巻き戻し1回ですべてのパラメータの正確な勾配を与えるので、数百〜数千のパルスセグメントの規模では唯一の現実的な道です。
微分可能な開放系オプティマイザで「勾配がノイズを貫いて流れる」とは正確には何を意味しますか?
- A.局所最小値から脱出するためにオプティマイザが勾配にランダムノイズを加える
- B.ノイズのあるハードウェア実行を多数回行って勾配を平均する
- C.損失の計算経路にリンドブラッドのノイズ項(T1、Tφ、漏れ)がautogradグラフの中に含まれ、∂C/∂θがノイズの効果まで反映する
- D.ノイズはランダムなので勾配から相殺されて消える
答えを見る
答え: C. 損失の計算経路にリンドブラッドのノイズ項(T1、Tφ、漏れ)がautogradグラフの中に含まれ、∂C/∂θがノイズの効果まで反映する
理由: 緩和・位相緩和・漏れの項を含むリンドブラッド発展全体が微分可能な演算で構成されているため、計算された勾配は各パラメータがノイズ世界で忠実度をどう変えるかをすでに反映しています — 事後のe^{−T/T1}割引との本質的な違いです。
前提となる概念
GRAPE gradient core is established (Khaneja 2005); the differentiable open-system tooling landscape incl. GradPulse case study (PureStateLabs/gradpulse, existence+claims verified on GitHub 2026-07-03) is current-state, hence timebound.
