想像してみてください。あなたは、極小の、目に見えないダンサーに、完璧なピルエット（回転）を教えようとしています。そのダンサーは肉や骨でできているのではなく、純粋なエネルギーと確率でできており、物理法則が風呂場の石鹸のように滑りやすい世界に生きています。これが量子コンピューティングの世界です。この領域では、「ダンサー」は量子ビット（qubit）と呼ばれ、「ピルエット」は論理ゲート、つまりコンピュータに思考させるための基本的な命令です。しかし、ここには落とし穴があります。これらのダンサーは非常に繊細なのです。もし強く押しすぎたり、速すぎたり、あるいはリズムを間違えたりすると、彼らはつまずき、バランスを崩し、パフォーマンス全体が混沌とした惨状に変わってしまいます。だからこそ、科学者には「量子最適制御（Quantum Optimal Control）」が必要です。それは、量子ビットを転ばせることなく、正確に目的地へと導くための、目に見えない「押し（パルス）」の完璧なシーケンスを設計する技術なのです。目標は、これらの動きを非常に精密かつ信頼できるものにし、今日のスーパーコンピュータが解くのに数千年かかるような問題を解決できるコンピュータを構築することです。
ここで、ドイツとスウェーデンの研究チームが登場します。彼らは、これらの量子ダンサーのために、新しいダンスインストラクターを試してみることにしました。彼らは、複雑な機械に歩き方や飛び方を教えるロボット工学の世界ですでにスター的存在である「反復線形二次レギュレータ（iLQR）」という手法を取り入れ、それを量子ゲートの振付を教える方法へと昇華させました。単に推測して試すのではなく、彼らの新しいアプローチは、この問題を「熱いか冷たいか」のゲームを逆再生して行うようなものとして扱います。まずゴール地点からスタートし、そこに到達するために「あるべきだった完璧な動き」は何かを導き出し、そこからステップごとに逆算して制御を調整していくのです。しかし、彼らは単にロボットの手法をコピー＆ペーストしたわけではありません。彼らは、その「ダンスの動き」をより滑らかにするために改良を加えました。現実の世界では、パルスを生成する装置（任意波形発生器）は、速度の急激でぎこちない変化を処理することができません。そのため、研究者たちは急激なジャンプにペナルティを与えるルールを追加し、アルゴリズムが滝のように崩れ落ちるのではなく、川の流れのようにスムーズな経路を見つけ出すように強制したのです。
彼らがシミュレーションで見出した結果は、非常に有望なものです。この新しい、滑らかにされたバージョンのiL_QRを使用することで、彼らは極めて正確な1量子ビットおよび2量子ビットゲートのパルスシーケンスを設計することができました。1量子ビットによる「Xゲート」（基本的な反転）において、彼らが達成した精度は、エラーがほとんど目に見えないほど小さく、10億回の試行につきわずか4回程度でした。2つのダンサーを結合させる、より複雑な2量子ビットの「クロスレゾナンス」ゲートにおいても、エラー率はわずかに高いものの、10万回につき約6回という高い忠実度の解を見出しました。最も素晴らしい点は、事前に「正解」を知る必要がなかったことです。優れた初期値を与えないと行き詰まってしまう他の手法とは異なり、この新しい手法は、たとえランダムで乱れた動きから始まったとしても、完璧な経路を見つけ出すことができるようです。
しかし、この物語には「ただし書き」があります。これらの結果は、現在はあくまでシミュレーション、つまりコンピュータ上でのデジタルなリハーサルであり、実際の量子マシン上でのパフォーマンスではありません。研究者たちは、彼らが見つけたパルスはまだ少し長く、完了までに最大240ナノ秒を要しており、これは可能なものと比較するとスローモーションのダンスのようなものであると認めています。また、彼らの手法はハードウェアに優しい非常に滑らかで穏やかなパルスを生み出しますが、既存の他の手法（L-BFGS-B GRAPEなど）は、非常に具体的で適切に選ばれた初期値を与えれば、時にはさらに精密な解を見つけることができるとも指摘しています。しかし、もし完璧な初期値が分からなくても、この新しいiLQRの手法は、ゼロからルーチンを組み立てることができる信頼できるコーチのような存在です。
チームの研究は、この「ロボット工学スタイル」のアプローチが、量子コンピューティングの未来における強力な新ツールになり得ることを示唆しています。制御問題を軌道最適化タスクとして扱い、エッジを滑らかにすることで、解を見るための水晶玉を持たずとも、完璧なゲートに非常に近づけることを彼らは証明しました。彼らはまだすべての問題を解決したわけではありません。現実世界のハードウェアは混沌としており、シミュレーションでは予測できない予期せぬ事態が常に起こります。しかし、彼らは新しい扉を開きました。彼らの次のステップは、このデジタルな振付を実際の量子チップ上でテストし、シミュレーションの中で見せたような美しいダンスを、現実の世界でも披露できるかどうかを確認することです。
技術要約：量子最適制御のための反復線形二次レギュレータ（iLQR）
問題提起
量子最適制御（QOC）は、量子システムに対して高忠実度のゲート操作を実現するために、精密な外部電場パルスシーケンスを設計することを目指しており、具体的にはエラー訂正の閾値を下回る忠実度を目標としている。GRAPE（Gradient Ascent Pulse Engineering）やKrotov法といった確立された手法が広く用いられているが、これらはしばしば特定の初期化戦略に依存したり、収束保証や複雑な制約の処理に関する制限があったりする。著者らは、ロボティクスにおける最先端の軌道最適化手法である反復線形二次レギュレータ（iLQR）を、量子ゲート合成の領域に適用する際の課題を特定している。課題は、iLQRを量子系の非線形ダイナミクス（シュレディンガー方程式によって支配される）に適応させつつ、滑らかなパルスエンベロープや任意波形発生器（AWG）に求められるゼロ境界条件といった、ハードウェア固有の制約を組み込むことにある。
手法
本論文では、特に固定周波数型トランモン量子ビットを対象とした、量子ゲート最適化のためのiLLRアルゴリズムの適応を提案している。その手法は、以下の主要な技術的コンポーネントで構成されている。
- 同型表現（Isomorphic Representation）: 量子発展の複素数値ユニタリ行列を扱うために、著者らは状態空間をベクトル化している。複素行列 U∈Cd×d は、実数値のベクトルおよび 2d2 次元の行列に写像され、標準的な実数値最適化手法の使用を可能にしている。
- パデ近似（Padé Approximation）: 高次系において計算コストが高くなる可能性がある正確な行列指数関数の代わりに、著者らはパデ近似（具体的には8次の冪級数）を利用して、前方伝播に必要な行列指数関数 exp(−iH(uk)Δt) を近似している。
- iLQR定式化: QOC問題は、実行コストと最終コストからなるコスト汎関数 J を最小化する軌道最適化問題として定式化される。このアルゴリズムは、値関数の微分とフィードバックゲイン（κk,Kk）を計算するための後方パス（backward pass）にガウス・ニュートン近似（ダイナミクスの二階微分を無視する）を用い、続いて制御軌道を更新するための前方パス（forward pass）を採用している。ヘッセ行列近似の可逆性を確保するために、レーベンバーグ・マルカートタ正規化が使用されている。
- 滑らかな制御制約: 急激な振幅変化が望ましくないハードウェアの制限に対処するため、著者らは制御変数を修正している。パルス振幅 uk を直接最適化するのではなく、パルスエンベロープの微分（u˙k）を最適化している。実際の制御は uk+1=uk+u˙kΔt を通じて更新される。この定式化により、コスト関数で変化率、パルス振幅、および最終値を罰することが可能となり、パルスがゼロから始まりゼロで終わることを保証している。
- システムモデル: 本手法は以下のモデルでテストされている：
- 単一量子ビット系: 2準位および3準位のエネルギーレベルを持つダフィング振動子としてモデル化。
- 2量子ビット系: Schrieffer-Wolf変換を用いて導出された有効ハミルトニアンを用いた、分散領域における結合トランモンとしてモデル化。2準位および3準位の両モデルがシミュレーションされている。
主な貢献
- アルゴリズムの適応: 量子ゲート最適化へのiLQRの適応に成功し、ロボティクスの軌道最適化とQOCの間の関連性を確立した。
- 滑らかなパルス合成: 微分ベースの制御定式化を導入することで、複雑な不等式制約を必要とせずに、滑らかなパルス形状とゼロ境界条件を本質的に強制する手法を提示した。
- 包括的なシミュレーション: 固定周波数型トランモンにおける X ゲート（単一量子ビット）およびクロス共鳴（CR）ゲート（2量子ビット）の高忠実度シミュレーション結果。これには、2準位および3準位（リーク）モデルの両方が含まれる。
- 比較分析: L-BFGS-B GRAPEとの直接比較を行い、収束挙動と初期パルス形状への感度の違いを明らかにした。
結果
著者らは、IBM Eagleプロセッサ（ibm_brisace）のパラメータを用いてゲート合成をシミュレートした：
- 単一量子ビット (2準位): X ゲートにおいて、微分制御を用いたiLQR法は、滑らかな放物線状のパルス形状で 4.0×10−9 の不完全度（infidelity）を達成した。これは、ガウス型の初期推定値を用いたL-BFGS-B GRAPEの解による 1.9×10−13 の不完全度と比較されている。
- 単一量子ビット (3準位): ∣2⟩ 状態へのリークを含めた場合、iLQRは 2.1×10−7 の不完全度を達成した。得られたパルスは、位相の遅れた成分が同相成分の微分に比例するというDRAG（Derivative Removal by Adiabatic Gate）の手順に似た特性を自然に示した。
- 2量子ビット (2準位): クロス共鳴ゲート（期間240 ns）において、iLQRは 1.1×10−8 の不完全度を達成した。
- 2量子ビット (3準位): リークを考慮した場合、不完全度は 5.9×10−5 であった。
- ハイパーパラメータの感度: 本研究では、二次コスト行列（Qf,Rd,Rc,Rf）の設計が、初期パルス形状よりも最終的な忠実度と滑らかさに大きな影響を与えることが判明した。
意義と主張
本論文は、事前の解析解が未知であるシナリオにおいて、iLQRが既存のQOC手法に対する堅牢な代替案となることを主張している。
- 初期化からの独立性: 初期パルス形状（例：ガウス型）に近い解を生成する傾向があるL-BFGS-B GRAPEとは異なり、iLQRのアプローチは、コスト関数が適切に調整されていれば、ランダムな初期化に関わらず滑らかで高忠実度な解を見出す。
- ハードウェア適合性: 微分を最適化することにより、本手法はAWGの制約（滑らかさ、開始/終了時のゼロ）を自然に遵守するパルスを生成する。これは実験的な実装において極めて重要である。
- 今後の方向性: 著者らは、現在のゲートは高忠実度ではあるものの、比較的長い（例：240 ns）ことを控えめに述べている。今後の課題は、時間依存のコスト関数を用いてより短いパルスを最適化すること、および超伝導システムに固有のモデルの不正確さを軽減するためにハードウェアフィードバックを統合することである。本手法は、可変周波数型トランモンを含む様々な超伝導アーキテクチャに適用可能であると提示されている。
毎週最高の quantum physics 論文をお届け。
スタンフォード、ケンブリッジ、フランス科学アカデミーの研究者に信頼されています。
受信トレイを確認して登録を完了してください。
問題が発生しました。もう一度お試しください。
スパムなし、いつでも解除可能。
週刊ダイジェスト — 最新の研究をわかりやすく。登録