ボースニック量子情報のための制御量子ビットとしてのフラクソニウム
Bosonic codes in superconducting resonators are a hardware-efficient avenue for quantum error correction and benefit from the inherent bias toward relaxation errors provided by long-lived cavities compared to typical...
背景と学術的系譜
起源と学術的系譜
本論文で取り上げられる問題は、超伝導共振器に格納されるボースニックコードを利用した量子コンピューティングの勃興分野における、堅牢な量子誤り訂正の極めて重要な必要性から生じている。これらのボースニックコードは、固有の長い記憶時間と偏ったエラー特性により、耐故障性量子コンピュータへのハードウェア効率的な道を提供する一方で、その実用的な実装は重大な課題に直面している。
歴史的に、これらのボースニックモードの制御と読み出しは、補助的な制御量子ビットとの結合に依存してきた。この目的で最も一般的に使用されてきた補助量子ビットは、トランスモン量子ビットである。しかし、トランスモンは、その広範な使用にもかかわらず、ボースニックコードの固有の利点を損なういくつかの有害な影響を導入する。これらの限界には以下が含まれる。
- 有限寿命: トランスモン量子ビットは有限のコヒーレンス時間を持つため、共振器操作の期間と忠実度が直接制限される。この有限寿命は、結合されたボースニック共振器に過剰なデコヒーレンスを誘発する可能性がある。
- 望ましくない非線形性: トランスモン量子ビットは、必然的に共振器に自己カー非線形性を誘発する。この望ましくない自己相互作用は、光子数依存の周波数シフトを引き起こし、望ましくない状態進化につながり、量子操作の忠実度を制限する。
- 新たなエラーチャネル: 結合自体が、パーセル減衰や量子ビット誘起デフェージングなどの新たなエラーチャネルを導入する可能性があり、回路性能をさらに低下させる。
トランスモンを用いたこれらの問題の軽減を目的とした過去の試みには、量子ビット-共振器結合の最小化またはチューニング、あるいはパラメトリック制御方式などの戦略が含まれていた。しかし、これらのアプローチは、制御速度の低下や回路設計・操作の複雑化の増加といったトレードオフを伴うことが多かった。既存のトランスモンベースのアーキテクチャのこの根本的な限界—高品位ボースニック共振器を「台無しにする」傾向—は、研究者たちに、高速で普遍的なボースニック制御を提供しつつ、これらの継承された有害な影響を排除できる代替補助量子ビットを模索することを強いた。本論文では、フラクソニウム量子ビットを有望な代替候補として調査する。これは、その長い寿命、フラックスチューニング可能性、および柔軟なハミルトニアン設計により、望ましくない非線形性を最小限に抑えるか排除するように量子ビット-共振器相互作用を調整できることに動機づけられている。
直感的なドメイン用語
- ボースニックコード: 情報が単純なオン/オフスイッチとしてではなく、瓶の中のマーブルの正確な数、あるいはグラスの中の水の正確な高さとして格納されると想像してください。各マーブルまたは水の高さは量子状態を表し、単一の「コンテナ」により複雑な情報を格納できます。
- 超伝導共振器: これは、完全に滑らかで超効率的なエコーチェンバー、または非常に高品質のベルと考えることができます。音(この場合はマイクロ波光子、光の小さなパケットのようなもの)を中に入れると、音はエネルギーを失うことなく非常に、非常に長い時間反響し、量子情報を格納するのに理想的な場所になります。
- 補助制御量子ビット: これは、特殊なリモコンや熟練したアシスタントのようなものです。それ自体は主要な情報を格納しませんが、主要な記憶装置(共振器)に直接触れることなく、その状態を正確に操作、読み出し、準備するために使用されます。
- トランスモン量子ビット: これは広く使用されてきた一般的なタイプのリモコン(補助量子ビット)です。一般的に信頼性がありますが、いくつかの欠点があります。それは比較的早く状態を「忘れる」(有限寿命)傾向があり、記憶装置と相互作用すると、エコーチェンバーが音を歪ませるように、記憶装置が望ましくない非線形な方法で動作するように誤って動作させることがあります。
- 自己カー非線形性: これは、音波が望ましくない方法で互いに干渉し始め、歪みを引き起こす「エコーチェンバー」(共振器)の欠陥のようなものです。各音波が独立して動作するのではなく、互いに影響を与え始め、全体的な音を正確に制御することが困難になります。これは、格納された情報の純粋さを損なう望ましくない自己相互作用です。
記法表
| 記法 | 説明 |
|---|---|
問題定義と制約
中核問題の定式化とジレンマ
本論文が取り組む中心的な問題は、量子誤り訂正の固有の利点を損なうことなく、超伝導共振器内のボースニックモードの堅牢な量子制御を達成する方法である。
入力/現在の状態:
超伝導共振器に実装されたボースニックコードは、量子誤り訂正へのハードウェア効率的な道を提供する。これらの共振器は長い量子状態記憶時間を誇り、緩和に対して強く偏ったエラーを示すため、誤り訂正の要件が簡素化される。これらのボースニックモードの普遍的な制御を達成するために、補助制御量子ビットが共振器に結合される。歴史的に、トランスモン量子ビットがこの役割の標準的な選択肢であった。
望ましい終点/目標状態:
最終的な目標は、ボースニック共振器の有効な読み出しと普遍的な制御能力の両方を提供する共振器-量子ビット結合を実現することである。極めて重要なのは、この結合が共振器を「台無しにする」有害な影響を導入しないこと、すなわち、共振器の長いコヒーレンス時間と有利なエラーバイアスを維持することである。具体的には、望ましい結果は、高速制御のための大きな分散シフト ($\chi$) を可能にすると同時に、自己カー非線形性 ($K$) を排除または大幅に抑制し、追加のデコヒーレンスを最小限に抑えるシステムである。
欠落しているリンクと数学的ギャップ:
正確な欠落しているリンクは、高速で普遍的なボースニック制御のための強い分散結合 ($\chi$) を提供すると同時に、自己カー非線形性 ($K$) とデコヒーレンスを最小限に抑えるか排除できる制御量子ビットアーキテクチャである。トランスモン量子ビットは、制御を提供するものの、本質的に $\chi$ と $K$ の間に制約されたトレードオフを導入する。ここで、高い $\chi$ を達成すると、しばしば大きな $K$ が生じ、その逆も同様であるか、あるいは共振器コヒーレンス(強いハイブリダイゼーション)を損なう領域での動作が必要となる。数学的なギャップは、特に高い $\chi/K$ 比を可能にするシステムハミルトニアン(またはシステム内の調整可能なパラメータのセット)を見つけることである。理想的には $K \approx 0$ である。本論文は、フラクソニウム量子ビットを調査することでこのギャップを埋めようとする。フラクソニウムの設計柔軟性は、この達成を可能にするようにハミルトニアンを調整できると仮定されている。
ジレンマ:
過去の研究者を閉じ込めてきた中心的なジレンマは、強力で高速な量子制御の達成と、ボースニック共振器の量子コヒーレンスと線形性の維持との間の痛みを伴うトレードオフである。補助量子ビット(トランスモンのような)が制御のために共振器に結合されると、通常、2つの主要な有害な影響を導入する。
1. 過剰なデコヒーレンス: トランスモン量子ビットの有限寿命(一般的に共振器の寿命よりも短い)は、新たなエラーチャネルとして機能し、結合された共振器にデコヒーレンスを誘発し、その長い記憶時間を無効にする。
2. 望ましくない非線形性: トランスモン量子ビットの非調和性は共振器に継承され、自己カー非線形性 ($K$) として現れる。この非線形性は、光子数依存の共振器周波数シフトを引き起こし、望ましくない状態進化につながり、操作忠実度を制限する。
これらの問題を軽減するための過去の試み(最小化またはチューニング可能な量子ビット-共振器結合やパラメトリック制御方式の使用など)は、しばしば制御速度の低下または回路設計・操作の複雑化の増加につながる。トランスモンシステムは、$\chi$ と $K$ が互いにほぼ比例するという根本的な制約に直面しており、強い量子ビット-共振器ハイブリダイゼーションの領域に入ることなく高い $\chi/K$ 比を達成することが困難である。これはそれ自体が共振器コヒーレンスの維持には望ましくない。これは、制御の改善(より高い $\chi$)がしばしば非線形性 ($K$) の増加またはコヒーレンスの低下という代償を伴うことを意味し、永続的で困難なトレードオフを生み出す。
制約と失敗モード
高忠実度、低エラーのボースニック量子制御を達成するという問題は、いくつかの厳しく現実的な制約によって非常に困難になっている。
-
物理的制約:
- 量子ビット寿命とデコヒーレンス: 補助量子ビットの有限寿命は主要な制約である。比較的短い寿命を持つトランスモン量子ビットは、結合された共振器に過剰なデコヒーレンスを誘発し、その長いコヒーレンス利点を損なう。フラクソニウムは、報告されているミリ秒単位の寿命[36-38]のために調査されている。これは、量子ビット誘起共振器デコヒーレンスを最小限に抑える可能性がある。
- 共振器寿命: 理想的な量子ビットであっても、記憶共振器自体は有限のエネルギー緩和時間 ($T_{1,s}$) を持つ。平面プロトタイプデバイスでは、単一光子寿命は12 µsと測定された。この比較的短い寿命は、状態準備とトモグラフィーの忠実度に著しく影響を与える。なぜなら、完全な実験シーケンスは、この時間の相当な割合を占める可能性があるからである。これは2Dオンチップアーキテクチャの実用的な限界であるが、3D実装はより高い品質係数を提供する。
- 継承された非線形性: 補助量子ビットの非調和性は、本質的に共振器に自己カー非線形性 ($K$) を付与する。これは、光子数依存の周波数シフトと望ましくない状態進化を引き起こす量子ビット設計(例:トランスモン)の根本的な物理的特性である。目標は、この継承された非線形性が最小化または排除されるシステムを設計することである。
- 量子ビット-共振器ハイブリダイゼーション: 強い結合は制御のために望ましいが、量子ビットと共振器の状態間の強いハイブリダイゼーション(例:量子ビットと共振器の周波数が近すぎる場合)は、共振器のコヒーレンスと線形性を損なう可能性があるため、一般に望ましくない。システムは、低いハイブリダイゼーションの領域で動作する必要がある。
- フラックスチューニング可能性: 外部磁束を介してシステムパラメータをその場で調整できる能力は、重要な制御ノブである。フラクソニウムのフラックスチューニング可能性は重要な利点であり、システム特性の動的な再構成を可能にするが、正確なフラックス制御の必要性も導入する。
-
計算上の制約:
- 複雑なハミルトニアンモデリング: 結合された量子ビット-共振器システムの挙動を正確に予測するには、複雑なハミルトニアン(例:式(1))の数値対角化と、マスター方程式(式(C5))を用いたシステムダイナミクスのシミュレーションが必要である。これには、かなりの計算リソースと洗練されたモデリングツール(scQubitsやQuTiPなど)が必要となる。
- パラメータチューニングと最適化: 理論モデルと実験データの間の定量的一致を達成し、高忠実度制御ゲート(SNAPゲートのような)を設計するには、多数の回路パラメータの正確なチューニングとパルスシーケンスの数値最適化が必要である。これは計算集約的であり、初期条件に敏感である可能性がある。
-
データ駆動型制約:
- 測定選択性とコヒーレンス: フラックススイートスポットから離れると、量子ビットのコヒーレンスが低下するため、光子数選択的な測定を実行する能力が制限され、変位キャリブレーションと共振器状態の特性評価がより困難になる。これは、$\chi$ や $K$ のような抽出されたパラメータの精度に影響を与える。
- 二準位系(TLS)からの時間的変動: 超伝導回路内のTLSの存在は、システムパラメータ(例:量子ビット分散シフト $\chi$、量子ビット $T_1$ および $T_2$)の時間的変動を引き起こす可能性がある。数日間にわたる時間スケールで観測されるこれらの変動は、キャリブレーションと測定の安定性と信頼性に影響を与え、一貫した高忠実度操作を困難にする。
- 量子ビット初期化の非効率性: 補助量子ビットの基底状態への不完全な初期化(例:熱平衡時の約66%の基底状態人口)は、全体的な不忠実度に寄与する。これには特定の冷却と事後選択手順が必要であり、実験に複雑さと時間を追加する。
- 読み出し背景とクロストーク: 共振器のクロストークと測定ベース冷却後の量子ビット減衰は、読み出し信号に歪んだ背景をもたらす可能性があり、真の量子ビット信号を分離するために慎重な補正手順が必要となる。
- 実験的分解能限界: 小さなカー非線形性 ($K$) の値を分解する能力は、実験ノイズ、量子ビットコヒーレンス、およびラムゼーフリンジに必要なシフトによって制限される。例えば、本論文では最小分解可能な $K/2\pi \approx 300$ Hz と推定しており、これは現在の方法ではより小さな非線形性を信頼性をもって検出および特性評価することが困難であることを意味する。これは、$K$ の抑制を検証できる精度を制限する。
FIG. 1. Fluxonium-resonator system. (a) Schematic compar- ison of an idealized cavity quantum electrodynamics system versus a circuit implementation. Our target is a dispersive shift χ as the only interaction (left). In practice, a resonator coupled to an artificial atom implemented by a superconducting circuit also inherits (at least) a self-Kerr nonlinearity K arising from the qubit’s anharmonicity (right). With sufficient control over cir- cuit parameters, K can be tuned and potentially suppressed. (b) False-colored optical images of the fabricated device. A fluxo- nium qubit is capacitively coupled to storage and readout modes, implemented as coplanar waveguide resonators
なぜこのアプローチなのか
選択の必然性
著者らがフラクソニウムを補助量子ビットとして探求することを選択したのは、ボースニック量子情報処理に広く採用されているトランスモン量子ビットに固有の根本的な限界の直接的な結果であり、恣意的なものではなかった。本論文は明確に、「一般的に使用されるトランスモン量子ビットは…有限寿命と必然的に誘発される自己カー非線形性により、共振器操作に限界を課す」(p.1, 要旨)と述べている。この認識は、主にトランスモンベースのシステムである従来の「最先端」(SOTA)手法が、ボースニック共振器に「過剰なデコヒーレンスや望ましくない非線形性などの非常に有害な影響」(p.1, 要旨)を導入したため不十分であったことを強調している。具体的には、「トランスモンを使用すると、制御不能な非線形性は避けられない」(p.2, 序論)ため、それらを軽減しようとすると制御速度の低下や回路の複雑化の増加といった問題が生じる。中心的な問題は、これらの継承された有害な影響なしに、高速で普遍的なボースニック制御を可能にする量子ビットを見つけることだった。フラクソニウムは、その独特の特性により、これらの重大な欠点を克服するための唯一の実行可能な候補として浮上した。
比較優位性
フラクソニウム量子ビットは、このアプリケーションにおいてトランスモンよりも圧倒的に優れた、いくつかの質的および構造的な利点を提供する。第一に、フラクソニウムは「ミリ秒単位の寿命[36-38]」を実証しており、これはトランスモンにおける重大な問題である量子ビット誘起共振器デコヒーレンスを最小限に抑えるために重要である(p.2)。これは、ボースニック共振器の長い記憶時間のより良い維持に直接つながる。第二に、その「フラックスチューニング可能性により、量子ビット-共振器結合のその場での調整が可能[39,40]」(p.2)であり、トランスモンよりも制約の少ない動的なシステム特性制御を提供する。最も重要なのは、フラクソニウムの「豊かなエネルギー準位構造」と、さまざまな回路パラメータ[41,42]を通じてハミルトニアンを変更できる能力により、「望ましくない非線形性が最小化または排除されるように、実効的な量子ビット-共振器ハミルトニアンを調整できる」(p.2)ことである。これは深遠な構造的利点である。本論文は、フラクソニウム-共振器システムが「トランスモン-共振器システムで可能なことを大幅に超える $\chi/K$ 比」を達成できることを実証している(p.5, 結論)。高い $\chi$(分散シフト)は高速操作に望ましく、低い $K$(自己カー非線形性)は制御忠実度を低下させる光子数依存デフェージングを防ぐために不可欠である。本質的な $\chi/K$ 関係の制約により、トランスモンでは単純に実現できない高忠実度制御を可能にする、「半分のフラックスで大きな $\chi$ を維持しながら消滅する $K$ を達成する」能力は、主要な差別化要因である(図5)。
制約との整合性
選択されたフラクソニウムアプローチは、ボースニック量子情報処理に定められた厳格な制約と完全に一致している。主な目標は、ボースニックコードの固有の利点(長い記憶時間や偏ったエラーチャネルなど)を維持しながら、オシレータの普遍的な量子制御を達成することであり、極めて重要なのは、補助量子ビットからの有害な影響を導入しないことである。フラクソニウムのユニークな特性は、これらの要件との「結婚」を形成する。
- デコヒーレンスの最小化: 共振器コヒーレンスを維持するという制約は、フラクソニウムの「長い寿命」(ミリ秒範囲)によって満たされ、これは量子ビット誘起共振器デコヒーレンスを直接低減する(p.2)。これは、トランスモンの有限寿命の限界からの大幅な改善である。
- 非線形性の制御: 主要な制約は、光子数依存デフェージングを引き起こす自己カー効果などの望ましくない非線形性の抑制である。フラクソニウムの「ハミルトニアンの設計柔軟性」と「豊かなエネルギー準位構造」により、「望ましくない非線形性が最小化または排除されるように、実効的な量子ビット-共振器ハミルトニアンを調整できる」(p.2)。特に、大きな分散シフトを維持しながらカー非線形性を消滅させる能力は、トランスモンシステムを悩ませる制御不能な非線形性の問題を直接解決する。
- 高速で普遍的な制御: 高速で普遍的なボースニック制御の必要性は、フラクソニウムの「大きな分散結合」(p.2)の能力と、その優れた $\chi/K$ 比によってサポートされており、デフェージングを最小限に抑えながら高速ゲート操作を可能にする。フラックスチューニング可能性は、特定の制御プロトコルに最適化するために、結合のその場での調整を可能にすることで、これをさらに強化する。
代替案の却下
本論文は、高性能ボースニック制御のための主要な代替案としてトランスモン量子ビットを暗黙的かつ明示的に却下している。トランスモンはボースニックコード機能[19-26]に広く使用されてきたが、著者らの中心的な動機は、その固有の限界から生じている。本論文は、トランスモンが「有限寿命と必然的に誘発される自己カー非線形性により、共振器操作に限界を課す」(p.1, 要旨)と述べている。これらは、フラクソニウムが克服するように設計されている根本的な問題である。
トランスモンベースのシステムで言及されている他のアプローチ、「最小化またはチューニング可能な量子ビット-共振器結合[9,12,29,30]またはパラメトリック制御方式[27,30–34]」も不十分と見なされている。本論文は、これらの代替案が「制御速度の低下、または回路設計や操作の複雑化の増加につながる可能性がある」(p.2)と指摘している。これは、それらが特定の問題に対する部分的な解決策を提供する可能性がある一方で、全体的な高速、高忠実度、単純なボースニック制御の目標を損なう新しいトレードオフを導入することを示唆している。対照的に、フラクソニウムは、これらの追加の複雑さや速度低下なしに優れたパフォーマンスを達成するために、量子ビットの特性を根本的に変更することにより、よりエレガントで直接的なソリューションを提供する。本論文は、この特定の共振器QEDコンテキストにおける補助量子ビットの役割として、GANsや他の量子ビットタイプ(例:トポロジカル量子ビット)のような他の一般的な量子コンピューティングアーキテクチャについては、直接の代替案として議論していない。
FIG. 4. Bosonic control using the fluxonium. (a) Pulse sequence for the preparation and characterization of Fock states in the storage resonator. A selective number-dependent arbitrary phase gate is used to prepare specific Fock states, which are characterized using either qubit spectroscopy or Wigner tomography. (b) Fluxonium spectroscopy with the storage resonator prepared in |1⟩(top) and 1 √
数学的・論理的メカニズム
マスター方程式
システムダイナミクス(コヒーレント進化とデコヒーレンスの両方を含む)を支配する絶対的な中心方程式は、リンドブラッドマスター方程式である。この方程式は、本論文でフラクソニウム-共振器システムの性能と限界を理解する上で中心となる開量子系のモデリングに不可欠である。具体的には、本論文では付録Cの式(C5)として言及している。
$$ \frac{\partial \hat{\rho}}{\partial t} = -i[H_{\text{eff,rot}} + H_{\text{drive,rot}}, \hat{\rho}] + \sum_k [L_k \hat{\rho} L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \hat{\rho}\}] $$
項別解剖
この方程式を分解し、各成分の数学的定義、物理的/論理的役割、およびその使用の根拠を説明する。
-
$\frac{\partial \hat{\rho}}{\partial t}$
- 数学的定義: 密度演算子 $\hat{\rho}$ の偏時間微分である。
- 物理的/論理的役割: システムの量子状態の変化の瞬間的な変化率を表す。本質的に、それは時間とともにシステムの状態がどのように進化するかを示す。
- なぜ使用されるか: これは、特に環境と相互作用する開量子系を考慮する場合に、量子系の時間進化を記述するための標準的な数学的表現である。
-
$\hat{\rho}$
- 数学的定義: 密度演算子(または密度行列)。正半定値エルミート演算子であり、トレースは1である。
- 物理的/論理的役割: 密度演算子は、量子系の状態の完全な統計的記述を提供する。それは純粋な量子状態(システムが確定した量子状態にある場合)と混合状態(システムが量子状態の古典的な確率的混合にある場合)の両方を表すことができ、デコヒーレンスにさらされた現実的なシステムのモデリングに不可欠である。
- なぜ使用されるか: 開量子系の場合、環境との相互作用により、システムが純粋な状態にあることはまれである。密度演算子は、これらの混合状態とその進化を記述するための適切なツールである。
-
$-i[H_{\text{eff,rot}} + H_{\text{drive,rot}}, \hat{\rho}]$
- 数学的定義: 交換子項であり、$[A, B] = AB - BA$ である。ここで、$H_{\text{total,rot}} = H_{\text{eff,rot}} + H_{\text{drive,rot}}$ は回転系における全ハミルトニアンである。したがって、項は $-i(H_{\text{total,rot}}\hat{\rho} - \hat{\rho} H_{\text{total,rot}})$ である。
- 物理的/論理的役割: 方程式のこの部分は、量子状態のコヒーレント、ユニタリ、可逆な進化を記述する。それは、内部エネルギー(ハミルトニアン)と適用される外部制御場によるシステム状態の変化を決定する。これは、純粋状態に対するシュレーディンガー方程式に類似している。
- なぜ使用されるか: 交換子形式は、進化がユニタリであることを保証し、散逸がない場合に密度行列のトレース(全確率)と状態の純粋性を維持する。虚数単位 $i$ は、時間進化を記述する量子力学の基本である。
-
$H_{\text{eff,rot}}$
- 数学的定義: 回転系における実効ハミルトニアンであり、$H_{\text{eff,rot}}/\hbar = \chi \hat{a}^\dagger \hat{a} |e\rangle \langle e| + \frac{K}{2} \hat{a}^\dagger \hat{a}^\dagger \hat{a} \hat{a}$ で与えられる。
- 物理的/論理的役割: このハミルトニアンは、分散領域におけるフラクソニウム量子ビット(二準位系としてモデル化)と超伝導記憶共振器との間の固有の静的な相互作用を記述する。それは、分散結合と自己カー非線形性の主要な物理現象を捉える。
- なぜ使用されるか: この簡略化された形式は、強い分散結合条件の下で有効であり、量子ビット-共振器相互作用の明確で直感的な理解を可能にする。回転系への変換は、高速振動項を除去し、解析を簡略化する。
-
$\chi \hat{a}^\dagger \hat{a} |e\rangle \langle e|$
- 数学的定義: 分散シフト項。$\chi$ は分散結合強度、$\hat{a}^\dagger$ および $\hat{a}$ は共振器の生成および消滅演算子、$|e\rangle \langle e|$ は量子ビットの励起状態への射影演算子である。
- 物理的/論理的役割: この項は、量子ビットが励起状態 $|e\rangle$ にあるときに記憶共振器の共振周波数をシフトさせる。逆に、量子ビットの遷移周波数は、共振器内の光子数に依存してシフトする。この相互作用は、量子ビット読み出しとボースニック共振器の量子制御の実装の基盤である。
- なぜ使用されるか: この項は、量子ビットと共振器間の強い分散結合から自然に生じる。これは、システムへのエネルギーシフトを表すため、ハミルトニアンの加算項である。
-
$\frac{K}{2} \hat{a}^\dagger \hat{a}^\dagger \hat{a} \hat{a}$
- 数学的定義: 自己カー非線形性項。$K$ はカー係数、$\hat{a}^\dagger \hat{a}^\dagger \hat{a} \hat{a}$ は光子数演算子 $\hat{n} = \hat{a}^\dagger \hat{a}$ を用いて $\hat{n}(\hat{n}-1)$ と書き換えられる演算子である。
- 物理的/論理的役割: この項は、共振器の周波数が既に共振器内に存在する光子の数に依存する非線形性を導入する。これは、異なる光子数状態が異なる速度で位相を蓄積させる原因となり、望ましくない状態進化(例:コヒーレント状態の圧縮)を引き起こし、量子操作の忠実度を制限する。$1/2$ の係数は、カー項の従来のスケールである。
- なぜ使用されるか: この項は、非調和フラクソニウム量子ビットとの相互作用から共振器が継承する最低次の非線形性を捉える。これはシステムへの追加のエネルギー寄与を表すため、ハミルトニアンの加算項である。
-
$H_{\text{drive,rot}}$
- 数学的定義: 回転系における駆動ハミルトニアンであり、$H_{\text{drive,rot}}/\hbar = \sum_n e_n(t)e^{-i\delta_n t} |e\rangle \langle g| + \text{h.c.}$ で与えられる。
- 物理的/論理的役割: このハミルトニアンは、フラクソニウム量子ビットを能動的に制御するために適用される外部マイクロ波場を記述する。これらの駆動は、量子ビットの基底状態 ($|g\rangle$) と励起状態 ($|e\rangle$) 間の遷移を誘発するように設計されており、ラビ振動や選択的位相ゲートなどの操作を可能にする。
- なぜ使用されるか: 外部制御は量子状態を操作するために不可欠である。この項は、量子ビットを駆動する精密に整形されたマイクロ波パルスのモデリングを可能にし、それが分散結合された共振器を制御する。
-
$\sum_n e_n(t)e^{-i\delta_n t} |e\rangle \langle g| + \text{h.c.}$
- 数学的定義: これは、$n$ 番目のマイクロ波駆動の複素振幅 $e_n(t)$、回転位相因子 $e^{-i\delta_n t}$、および量子ビット昇順演算子 $|e\rangle \langle g|$ の和である。「h.c.」はエルミート共役を表し、降順演算子 $|g\rangle \langle e|$ を含む。$\delta_n = \omega_n - \Omega$ は、$n$ 番目の駆動周波数 $\omega_n$ と量子ビット遷移周波数 $\Omega$ の間のデチューニングである。
- 物理的/論理的役割: 和の各項は、量子ビットに適用される特定のマイクロ波駆動を表す。$e_n(t)$ はパルスの強度と時間形状を制御し、$e^{-i\delta_n t}$ は回転系に対する駆動の位相進化を考慮し、$|e\rangle \langle g|$ は基底状態から励起状態への遷移(および励起状態から基底状態への遷移のためのその共役)を誘発する。これにより、量子ビット状態の精密な操作が可能になる。
- なぜ使用されるか: 和は、複数の、潜在的に異なる駆動トーンの適用を可能にする。指数関数的な位相因子は、回転系への変換の結果であり、共振時($\delta_n=0$)には駆動が実質的に時間不変になる。
-
$\sum_k [L_k \hat{\rho} L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \hat{\rho}\}]$
- 数学的定義: これはリンドブラッド超演算子であり、量子システムのコヒーレントでない、散逸的な進化を表す。$\{A, B\} = AB + BA$ は反交換子である。
- 物理的/論理的役割: この項は、システムが環境と不可逆的に相互作用することによるデコヒーレンスとエネルギー緩和を考慮する。各 $L_k$ は、特定の損失メカニズム(例:光子損失、量子ビット緩和、デフェージング)に対応する崩壊演算子である。それはシステムを混合状態に向かわせ、量子コヒーレンスを低減する。
- なぜ使用されるか: リンドブラッド形式は、開量子系のマルコフ的、トレース保存、完全正の進化を記述する最も一般的な方法である。和はすべての関連する独立した損失チャネルを含み、特定の形式は密度行列が物理的に有効であることを保証する。
-
$L_k$
- 数学的定義: 崩壊演算子。本論文では以下のように指定されている。
- 記憶モードエネルギー減衰: $L_1 = \sqrt{\kappa} \hat{a}$
- フラクソニウム緩和: $L_2 = \sqrt{\Gamma_1} |g\rangle \langle e|$
- フラクソニウム励起: $L_3 = \sqrt{\Gamma_1} |e\rangle \langle g|$
- 純粋デフェージング: $L_4 = \sqrt{\Gamma_\phi/2} (|e\rangle \langle e| - |g\rangle \langle g|)$
- 物理的/論理的役割: 各 $L_k$ は、環境との特定の種類の不可逆相互作用を記述する。$\sqrt{\kappa} \hat{a}$ は共振器からの光子損失をモデル化する。$\sqrt{\Gamma_1} |g\rangle \langle e|$ は、量子ビットが励起状態から基底状態にエネルギーを失うことを表す。$\sqrt{\Gamma_1} |e\rangle \langle g|$ は、量子ビットの熱励起を考慮する。$\sqrt{\Gamma_\phi/2} (|e\rangle \langle e| - |g\rangle \langle g|)$ は、エネルギー交換なしに位相情報が失われる純粋デフェージングをモデル化する。
- なぜ使用されるか: これらの演算子は、超伝導量子回路で観測される主要な物理的損失メカニズムを正確にモデル化するために選択される。平方根係数は、これらの散逸プロセスの正しいレートを保証する。
- 数学的定義: 崩壊演算子。本論文では以下のように指定されている。
ステップバイステップの流れ
密度行列 $\hat{\rho}$ で表される単一の抽象データポイントの正確なライフサイクルが、この数学的エンジンを通過する様子を追ってみよう。
- 初期状態注入: プロセスは、時刻 $t=0$ におけるフラクソニウム-共振器システムの量子状態を表す初期密度行列 $\hat{\rho}(0)$ から始まる。これは、不完全な初期化による基底状態、コヒーレント状態、または混合状態である可能性がある。
- コヒーレント進化組立ライン:
- 時間が進むにつれて、マスター方程式のコヒーレント部分、$-i[H_{\text{eff,rot}} + H_{\text{drive,rot}}, \hat{\rho}]$ が $\hat{\rho}$ に作用する。
- 固有の相互作用: まず、静的な実効ハミルトニアン $H_{\text{eff,rot}}$ が $\hat{\rho}$ の形状を形成し始める。分散シフト項 $\chi \hat{a}^\dagger \hat{a} |e\rangle \langle e|$ は、フラクソニウム量子ビットが基底状態か励起状態かによって、共振器の周波数が「回転」するのを変える。これにより、状態依存の周波数シフトが作成される。
- 非線形位相蓄積: 同時に、自己カー非線形性 $\frac{K}{2} \hat{a}^\dagger \hat{a}^\dagger \hat{a} \hat{a}$ は、共振器内で光子数依存の位相シフトを誘発する。$\hat{\rho}$ がコヒーレント状態を表す場合、この項は異なる光子数成分が異なる速度で進化するため、位相空間で「せん断」または歪みを引き起こす。
- 外部制御入力: 次に、外部マイクロ波パルス($H_{\text{drive,rot}}$ で記述)が適用される。これらのパルスは、フラクソニウム量子ビットの遷移(例:$|g\rangle \leftrightarrow |e\rangle$)を誘発するように、正確なタイミングと形状で整形される。量子ビットを操作することにより、分散結合された共振器の状態を間接的に制御できる。例えば、特定の位相シフトを付与したり、状態転送を実行したりすることによって。
- これらのコヒーレント操作は $\hat{\rho}$ をユニタリ変換し、ヒルベルト空間内のその成分を予測可能で可逆な方法で回転させる。
- コヒーレントでない減衰とデフェージングステーション:
- コヒーレント進化と同時に、マスター方程式のコヒーレントでない部分、$\sum_k [L_k \hat{\rho} L_k^\dagger - \frac{1}{2}\{L_k^\dagger L_k, \hat{\rho}\}]$ が $\hat{\rho}$ を継続的に処理する。
- 光子損失: 共振器の崩壊演算子 $\sqrt{\kappa} \hat{a}$ は、共振器からの光子の脱出をモデル化する。これにより、共振器内の総光子数が減少し、量子状態が不可逆的にエネルギーを失う。
- 量子ビット緩和/励起: 量子ビットの崩壊演算子 $\sqrt{\Gamma_1} |g\rangle \langle e|$ および $\sqrt{\Gamma_1} |e\rangle \langle g|$ は、量子ビットが環境にエネルギーを失う(緩和)またはエネルギーを得る(熱励起)ことを表す。これらのプロセスは、量子ビットの状態が熱平衡に向かって減衰する原因となる。
- 純粋デフェージング: 純粋デフェージング演算子 $\sqrt{\Gamma_\phi/2} (|e\rangle \langle e| - |g\rangle \langle g|)$ は、エネルギー交換なしに量子ビットが量子位相情報を失う原因となる。これは重ね合わせを破壊し、量子ビットを混合状態に向かわせる。
- これらのコヒーレントでないプロセスは、時間とともに $\hat{\rho}$ をより「混合」させ、対角外要素(コヒーレンス)を低減し、システムを古典的な統計的混合に向かわせる。
- 出力状態: 一定の進化時間 $t$ の後、マスター方程式は最終密度行列 $\hat{\rho}(t)$ を生成する。この $\hat{\rho}(t)$ は、意図されたコヒーレント変換と避けられないコヒーレントでない劣化の両方を経たシステムの量子状態を表す。この最終密度行列から、量子ビットの人口確率、ウィグナー関数、または忠実度メトリックなどの観測可能な量を計算して、システムの性能を特徴付けることができる。
最適化ダイナミクス
リンドブラッドマスター方程式によって記述されるメカニズムは、典型的な機械学習の意味での「学習」アルゴリズムではなく、量子システム性能の特性評価、制御、および最適化のための強力なツールである。ここでの「最適化」は、望ましい量子状態を達成するかエラーを最小限に抑えるために、最適なシステムパラメータまたは制御パルスシーケンスを見つけることを指す。
-
パラメータ抽出と損失地形マッピング:
- このメカニズムの主な用途は、システムのダイナミクスをシミュレーションすることである。これらのシミュレーションを実験データ(例:量子ビット分光、ラムゼー干渉測定、ウィグナートモグラフィー)と比較することにより、著者らは分散シフト $\chi$、カー非線形性 $K$、およびさまざまな減衰率($\kappa$、$\Gamma_1$、$\Gamma_\phi$)などの重要な物理パラメータを抽出する。
- これにはフィッティングプロセスが含まれる。モデルパラメータを変化させることでシミュレーション結果が調整され、実験データに最もよく一致するようになる。この文脈では、「損失地形」はトレーニング損失関数ではなく、モデルの予測が実験的観測と一致する度合いの地形である。目標は、不一致を最小限に抑えるパラメータ値を見つけることであり、通常は最小二乗法または同様のフィッティングアルゴリズムが使用される。
- 本論文では、システム損失を考慮してマスター方程式を用いてラムゼーフリンジをフィッティングし、実験データから $\chi$ および $K$ を抽出することを明示的に言及している。
- フラクソニウムのユニークなフラックスチューニング可能性は、$\chi$ および $K$ の動的な制御を可能にする。これらのパラメータを外部フラックスの関数としてマッピングすることにより(図3)、著者らは、例えばカー非線形性 $K$ が抑制され、大きな分散シフト $\chi$ が維持される「スイートスポット」を特定できる。これは、高忠実度ボースニック制御の重要な最適化目標である。
-
制御パルス最適化:
- 特定の量子制御タスク(例:共振器フォック状態または重ね合わせ状態の準備、選択的数依存任意位相(SNAP)ゲートの使用)には、正確な制御パルスシーケンスが必要である。これらのパルスは、駆動ハミルトニアン $H_{\text{drive,rot}}$ における時間依存振幅 $e_n(t)$ とデチューニング $\delta_n$ によって定義される。
- 本論文では、「選択的 $\pi$ パルスの振幅と位相は、コヒーレントエラーを抑制するために数値的に最適化された」と述べている。この最適化プロセスには、異なるパルスパラメータでマスター方程式(C5)を順方向に繰り返し実行することが含まれる。各パラメータセットについて、結果として得られる量子状態の忠実度(目標状態にどれだけ近いか)が計算される。
- このシナリオにおける「損失地形」は、パルスパラメータの関数としての非忠実度(1 - 忠実度)である。数値最適化アルゴリズム(例:勾配降下法、進化アルゴリズム、または最適制御法)は、この地形をナビゲートし、望ましい量子操作の忠実度を最大化するパルスパラメータを見つけるために使用される。
- コヒーレントでないエラーのシミュレーション(図4(f))も、システム寿命($T_{1,s}$、$T_{\phi,f}$)がゲート忠実度にどのように影響するかを示すことで、この最適化を導き、将来のデバイス設計と実験プロトコルに情報を提供する。
-
反復状態更新:
- その核心において、マスター方程式自体は、微小時間ステップ $dt$ にわたる密度行列 $\hat{\rho}$ の連続的な反復更新を記述する。数値シミュレーションでは、これは通常、時間ステップ法(例:ルンゲ=クッタ法)を使用して解決される。各ステップで、$\hat{\rho}(t+dt)$ は、コヒーレントおよびコヒーレントでない進化演算子の組み合わせを適用することによって $\hat{\rho}(t)$ から計算される。この反復プロセスにより、さまざまな実験条件と制御パルス下での量子状態の完全なライフサイクルを追跡することができ、特性評価と最適化の両方の基礎を形成する。
結果、限界、結論
実験設計とベースライン
私たちの実験セットアップは、ボースニック量子情報のための制御量子ビットとしてのフラクソニウムの能力を検証するために細心の注意を払って設計された。私たちは、超伝導記憶共振器とフラクソニウム量子ビットからなる最小回路を構築し、量子ビット状態測定のための追加の読み出し共振器を付加した(付録A 3で詳述)。この概念実証作業のために、私たちは2次元(2D)オンチップアーキテクチャを選択した。これは、共振器の品質係数とのトレードオフを認識しつつも、フラクソニウムと共振器の回路パラメータに対する正確な制御を可能にした。
フラクソニウムの最も有利な動作点は、その寿命とコヒーレンスが最大化される半フラックス領域で特定された。この時点で、私たちは量子ビットの固有特性を特徴付け、デバイスAの緩和時間 ($T_1$) を123 µs、ハーンエコーコヒーレンス時間 ($T_2$) を90 µsと測定した(表I)。量子ビット初期化は、測定ベースの冷却と事後選択手順によって達成され、一貫して90%を超える基底状態人口が得られた(付録B 1)。
量子ビット-共振器相互作用を特徴付けるために、私たちは変位パルスを使用して記憶モードに振幅 $\alpha$ のコヒーレント状態を準備し、次にフラクソニウムのスペクトルを測定した(図2(a))。$\chi/2\pi = 1.0$ MHz のピーク分離を持つ、明確な数分割量子ビット共鳴の観測は、強い分散結合の決定的な証拠となった。この強い結合は、ボースニック状態制御と読み出しの基盤となるフォック状態選択的量子ビット回転の基礎である。私たちはさらに、変位後の共振器が真空状態に戻る確率を監視することにより、共振器の緩和率を測定し、単一光子寿命を12 µsと抽出した(図2(b))。この値は、主に製造された2Dデバイスの内部品質係数によって制限されていた。
記憶モードの継承された自己カー非線形性 ($K$) を抽出することが中心的な目標であった。私たちは、「キャビティラムゼー」干渉実験を使用してこれを達成した。この実験では、一連の共振器変位と数選択的測定(図2(c))を使用した。より大きな変位振幅で観測されたラムゼーフリンジのシフトは、光子数依存デチューニングの直接的な兆候であり、そこから自己カー係数 $K/2\pi = 3.6$ kHz を抽出した。
極めて重要なのは、私たちは $\chi$ と $K$ のフラックス依存性を調査した。私たちは、さまざまな外部フラックス値にわたるフラクソニウムスペクトルを測定し、抽出された回路パラメータを使用して、$\chi$ と $K$ の予想されるフラックス依存性を数値的に計算した。この特性評価は、異なる記憶共振器周波数を持つ2つの異なるデバイス(AおよびB)で実行された(図3)。
最後に、私たちは共振器フォック状態とその重ね合わせ状態を準備および特徴付けることにより、ボースニック制御を実証した。これは、共振器変位操作と選択的数依存任意位相(SNAP)ゲート(図4(a))を使用して達成された。準備された状態は、量子ビット分光法とウィグナートモグラフィー(図4(b), 4(c))の両方を使用して特徴付けられた。この作業全体における私たちの主なベースライン比較対象は、超伝導回路で広く使用されている補助量子ビットであるトランスモン量子ビットであった。私たちは、フラクソニウム-共振器システムが、トランスモンベースのアーキテクチャと比較して、優れた $\chi/K$ 比を達成できることを厳密に実証することを目指した。トランスモンベースのシステムでは、$\chi^2/|K|$ の理論的上限が $|K|/2\pi \ge (\chi/2\pi)^2 / (2.12 \text{ GHz})$ として確立されていた(付録C3)。また、私たちの発見を文脈化するために、多数の文献値(表II)も参照した。
証拠が証明すること
私たちが収集した実験的証拠は、特にトランスモンベースのシステムに固有の限界を克服する能力において、フラクソニウムの高性能ボースニック制御量子ビットとしての可能性を説得力のある証明を提供している。
第一に、$\chi/2\pi = 1.0$ MHz の分散シフトを持つ明確な数分割量子ビット共鳴の観測(図2(a))は、強い分散結合を断言している。これは、ボースニック状態制御と読み出しに必要な、正確なフォック状態選択的量子ビット回転を可能にする、私たちのボースニック制御メカニズムの基盤である。記憶共振器の測定された単一光子寿命12 µs(図2(b))は、2Dプロトタイプによって制限されていたものの、意味のある期間量子情報を格納できることを確認している。
第二に、私たちの「キャビティラムゼー」干渉実験は、共振器の自己カー非線形性の決定的で否定できない証拠を提供した。ラムゼーフリンジの振幅依存シフト(図2(d))は、カー効果を直接明らかにし、 $K/2\pi = 3.6$ kHz を抽出することを可能にした。私たちの実験データとリンドブラッドマスター方程式シミュレーション(図2(e))との定量的一致は、重要な証拠である。それは、我々が非線形性と損失の相互作用に関する堅牢で数学的な理解を持っていることを証明しており、低エラーボースニック制御のためのデバイスパラメータを信頼性をもって予測することを可能にする。
第三に、フラクソニウムのユニークなチューニング可能性は、徹底的に証明された。外部フラックスの関数としての $\chi$ と $K$ の測定(デバイスAの図3(c), 3(e);デバイスBの図3(d), 3(f))は、数値予測と非常によく一致する、驚くほど広い変動を示した。最も重要なのは、データは $K$ が符号を変え、特定のフラックスバイアスポイントでゼロを横切ることができることを示したことである(図3(e))。このチューニング可能性はゲームチェンジャーであり、標準的なトランスモン量子ビットではほとんどアクセスできない、カー誘起状態進化を動的に抑制する能力を実証している。
第四に、私たちは共振器の状態の量子制御を首尾よく実証した。共振器変位とSNAPゲートを使用することにより、単一光子フォック状態 $|1\rangle$ と重ね合わせ状態 $(|0\rangle - |1\rangle)/\sqrt{2}$ の両方を準備および特徴付けた。実験的に測定されたウィグナー関数(図4(c))は、これらの状態の準備が成功したことを明確に示している。これらの測定されたウィグナー関数とマスター方程式シミュレーション(図4(d), 4(e))との強い定量的一致は、コヒーレントでないエラーの影響を含む、システムダイナミクスとエラーモデルの理解をさらに検証している。
最後に、そしておそらく最も重要なのは、証拠はフラクソニウム-共振器システムが、重要な指標においてトランスモンと比較して優れたパフォーマンスを達成できることを証明している。私たちの分析、特に $\chi^2/|K|$ 対フラックスの図(図5(d))は、フラクソニウムがいくつかのフラックス領域でシミュレートされたトランスモン境界を超えることができることを示している。デバイスBの場合、2つの測定点($\Phi_{\text{ext}} \approx 0.1919\Phi_0$ および $\Phi_{\text{ext}} \approx 0.3358\Phi_0$)は、それぞれ $\chi/2\pi = \{-1.83, -3.31\}$ MHz および $K/2\pi = \{0.31, 1.04\}$ kHz で、トランスモン境界を明確に上回っている。これは、フラクソニウムのコアメカニズムが分散結合とカー非線形性の間のより有利なバランスを可能にするという、高忠実度ボースニック量子情報処理の重要な利点であるという、決定的な、否定できない証拠である。半フラックスで大きな $\chi$ を持ちながら消滅する $K$ を達成するフラクソニウム回路のシミュレートされた例(図5(b))は、この可能性をさらに強化している。
限界と将来の方向性
私たちの仕事は、フラクソニウムをボースニック量子情報処理の有望な候補として確立しているが、現在のプロトタイプデバイスにおけるいくつかの限界は、達成可能な忠実度を制限しており、将来の開発のためのエキサイティングな方向性を示唆している。
主要な限界は、12 µs と測定された、私たちの2D共振器の比較的低い共振器寿命である。ウィグナートモグラフィーの約3.5 µs の完全な実験シーケンスは、この緩和時間の相当な割合(>25%)を表す。その結果、フォック状態 $|1\rangle$ (79%) および重ね合わせ状態 $(|0\rangle - |1\rangle)/\sqrt{2}$ (91%) の準備忠実度は、現在、この短い共振器寿命と量子ビット初期化の非効率性によって制限されている(表III)。耐故障性量子計算に必要な高忠実度を達成するには、システムコヒーレンス時間の顕著な改善が不可欠である。
工学的な観点から、高Q共振器(3Dアーキテクチャで使用されるものなど)とフラクソニウム量子ビットを統合することは、依然として大きな課題である。私たちの2Dプラットフォームは概念実証には優れていましたが、より高Qのシステムへのスケーリングには、回路パラメータを確実にターゲットにし、高コヒーレントなフラクソニウム量子ビットを製造するための製造プロセスのさらなる進歩が必要となる。
将来の探求のもう一つの領域は、高次の非線形性に関係している。私たちの現在のボースニック制御実験は主に少数の共振器光子で動作しており、そこでは主要なカー非線形性 ($K_2$) が支配的である。しかし、ボースニックコードがより多くの光子数を利用するように進化するにつれて、より高次のカー係数($p > 2$ の $K_p$)がますます関連性を持つようになるだろう。私たちのシミュレーションは、デバイスAおよびBの場合、約10個の光子後(図14(c))にこれらの高次の補正が重要になることを示唆している。将来の研究では、堅牢な制御のためにこれらの高次の項を軽減または利用する方法に対処する必要があるだろう。
将来を見据えると、フラクソニウムの固有のフラックスチューニング可能性($\chi$ および $K$)は、システム特性の動的な再構成のための強力なノブを提供する。これは、アイドル期間中のカー誘起状態進化を抑制したり、特定の量子プロトコル用にパラメータを動的に最適化したりするために利用できる。異なる $K/\chi$ 比を工学的に設計できる能力は、高速操作と最小限の非線形性のバランスをとる新しい制御スキームの設計の可能性を開く。
最終的な目標は、「カーフリーボースニック制御」を実現することである。これは、フラクソニウムが強い分散結合領域で共振器に結合され、自己カー非線形性が実質的に排除されることを意味し、特に半フラックスで顕著である。私たちのシミュレーションはこれが達成可能であると予測しており、その実験的実現は重要な次のステップである。これには、高Q共振器統合のさらなる進歩と、予測されたハミルトニアンパラメータを達成するための洗練された製造プロセスの開発が含まれるだろう。
直接的な改善を超えて、確立された制御技術は、高Q共振器を組み込んだより高度なアーキテクチャに直接転送可能であり、カーフリーボースニック制御の実現への道を開く。さらに、量子ビット初期化のためのアクティブリセット方法の改善や、よりコヒーレントな量子ビットおよび共振器システムの開発などの現在のエラー源に対処することは、コヒーレントでない損失を軽減するために重要となるだろう。
フラクソニウムのユニークな設計自由度は、その豊かなエネルギー準位構造とフラックスチューニング可能性に由来し、従来の共振器QEDシステムではアクセスできない、量子ビット-光子相互作用を調整するための比類のない機会を提供する。これは、将来の研究が、SNAPプロトコルを超えた新しいボースニック制御スキームを探求できることを示唆しており、新しいゲート設計、誤り訂正戦略、あるいはフラクソニウムの汎用性を完全に活用する量子情報処理の全く新しいパラダイムにつながる可能性がある。既存の量子ビット技術のバイアスから解放されたこの広範な視点は、批判的思考を刺激し、ボースニック量子コンピューティングの進化を推進するだろう。
FIG. 3. Hamiltonian parameters as function of external flux. (a) and (b) Fluxonium spectra of device A (left; solid circles) and device B (right; open circles), fit to the |g⟩→|e⟩transition. The higher transitions |g⟩→|f ⟩and |g⟩→|h⟩are shown in gray. The storage resonator level (orange) crosses different higher fluxonium levels in each device. This results in distinct flux dependence of the dispersive shift χ shown in (c) and (d), and the self-Kerr nonlinearity K shown in (e) and (f). K can change sign and cross zero at specific flux bias points (star; see Appendix D 1 for raw data). Solid lines are expected parameters based on the fit to the fluxonium spectrum. All simulations are performed numer- ically and show good agreement with the measurements. Error bars are smaller than the marker size (Appendix A 5)
FIG. 5. Relation of χ and K for the fluxonium and trans- mon. (a) K (purple line) and overlap between bare and dressed states (green line) as a function of for a (hypothetical) transmon-resonator system with parameters EC/2π = 530 MHz and EJ/2π = 26.5 GHz for fixed |χ|/2π = 1 MHz. Although K reaches zero near EC ≈ , the qubit and resonator states become strongly hybridized at this point. The black star shows example transmon parameters for the bound (black dashed line) on χ/K discussed in the main text. (b) The same quanti- ties as in (a), shown for a fluxonium-resonator system at half flux with EC/2π = 1.19 GHz, EL/2π = 556 MHz, EJ/2π = 3.04 GHz. Blue star indicates a set of fluxonium parameters that beat the transmon bound, chosen with a detuning of 100 MHz from the K = 0 case, to illustrate that no unrealistic fine-tuning is required for small Kerr nonlinearity. (c) Our derived χ/K bound (dashed black line) plotted with values from the literature [9,23–25,44,45,48–53] for transmons (black circles are measured data; black diamonds have measured χ and simulated K), and devices A (blue circle) and B (open blue circle) of this work (at half-flux). (d) The figure of merit χ2/K vs flux. The ratio is plotted for measured data (blue circles) and fit to simula- tion (blue line) for device A (device B) in the top (bottom) panel. The black dashed lines are the transmon bound. For device B, two measured points fall above the transmon bound, with χ/2π = {−1.83, −3.31} MHz and K/2π = {0.31, 1.04} kHz at ?ext ≈{0.1919, 0.3358}?0 (see Appendix D 3 for details)
FIG. 14. Higher-order Kerr in the fluxonium devices. (a) Sim- ulated pth order shift Kp as a function of detuning for the same fluxonium as presented in Fig. 5 for p = 2 (purple), 3 (blue), 4 (orange), and 5 (green). At each , g is tuned such that |χ| = 1 MHz. The blue star at /2π = −2.25 MHz is the same example point as in Fig. 5. (b) At the blue star point of sub-figure a, the simulated deviation ω (orange stars) of the cavity fre- quency is plotted for different photon numbers. The deviation only due to K2/2π = 1.76 Hz is plotted for comparison (dashed purple line). (c) The simulated deviation ω for devices A and B (closed and open circles, respectively). The deviations only due to K2 are plotted as dash-dotted and dotted lines, respectively