図で読む DUST 論文詳細レポート

W. Dai et al., “Characterization of drive-induced unwanted state transitions in superconducting circuits” の内容を、論文中の図を挿入して、実験・数理・設計上の意味がつながるように解説するレポート。

TopicDrive-Induced Unwanted State Transitions
SystemFixed-frequency 3D transmon
Core methodPump-probe spectroscopy + Floquet
Report date2026-05-26

1. 要旨: この論文は何を解いたのか

この論文の主張は明確である。超伝導量子回路では、ゲートや読み出しを高速化するためにマイクロ波ドライブを強くしたい。しかし強ドライブは、qubit を計算空間内外の別状態へ勝手に遷移させる。著者らはこの現象を DUST: drive-induced unwanted state transitions と呼び、実験スペクトル上で似て見える多数の feature を三つの物理機構に分解する。

ATLS/ac-Stark

ac-Stark shift により qubit 遷移周波数が TLS と共鳴し、エネルギー交換や緩和が起きる。

BIntrinsic multi-photon

トランズモンの非線形 Hamiltonian そのものが、複数光子吸収による高準位遷移を許す。

CExternal-mode scattering

ドライブがトランズモン遷移を起こし、余剰エネルギーを spurious EM mode や TLS へ放出する。

このレポートの読み方: まず Fig. 1 で三機構の物理像をつかみ、Fig. 2 で実測 map の軸と色を理解する。その後、Fig. 3 が A、Fig. 4-5 が B、Fig. 6-8 が C の証拠としてどう使われているかを見る。
図の出典について: 本文の Fig. 1-8 はすべて埋め込み済み。Appendix は公開 HTML レンダリングから取得できた主要補助図を入れている。アップロード PDF 側は期限切れだったため、現行 v3 の全 18 図を PDF から直接抽出するのではなく、公開 arXiv/ar5iv で取得できた図を使用した。
Fig. 1

DUST の三つの物理機構

Fig. 1: DUST の三つの物理機構

図の要点: 3D トランズモン、ac-Stark shift による TLS 共鳴、トランズモン固有の多光子遷移、外部モードを介した非弾性散乱を一枚でまとめた概念図。右下の Re Z_env のピーク列は、読み出しモードやパッケージ・キャビティの寄生電磁モードが「散乱先」になることを示す。

読み方:
  • この図が論文全体の地図である。DUST は単一の現象名ではなく、強ドライブ下で観測される不要遷移の総称であり、著者らはそれを Mechanism A/B/C に分解して読む。
  • A は qubit 周波数そのものが ac-Stark shift で動き、近傍 TLS と共鳴してエネルギー交換する機構。B は外部モードなしでもトランズモンの非線形性だけで起こる多光子遷移。C はドライブがトランズモン遷移を起こし、余剰エネルギーを別周波数のフォトンとして環境へ逃がす非弾性散乱である。
  • 実験的には、A は ac-Stark shift 軸でほぼ水平に現れる。B は transmon-only Floquet で予測できる。C は環境インピーダンスや RF/EM モードを見ないと説明できない。

2. 実験: DUST map は何を測っているのか

著者らの基本測定は time-resolved pump-probe spectroscopy である。qubit を |0_t> または |1_t> に初期化し、周波数 ω_d と強度を掃引した stimulation/pump pulse を印加し、その後に多準位 readout で最終人口を読む。

P_leave,j(ω_d, A_d) = 1 - P(final = j | initial = j, ω_d, A_d) P_leak,j(ω_d, A_d) = P(final in non-computational states | initial = j)

論文で特に重要なのは、縦軸を発生器 power のまま扱わず、qubit に実際に誘起された ac-Stark shift に変換している点である。これにより、周波数依存の伝達関数や detuning による drive efficiency の違いをある程度吸収できる。

Fig. 2

pump-probe spectroscopy による DUST map

Fig. 2: pump-probe spectroscopy による DUST map

図の要点: 初期状態を準備し、周波数・パワー可変の stimulation/pump pulse を入れ、待ち時間後に最終人口を読み出す。主図は初期 |1_t> から外れた確率 1 - P(1_t) を drive frequency と normalized ac-Stark shift の平面に描いた DUST map。

読み方:
  • 横軸はドライブ周波数、縦軸はドライブ強度を ac-Stark shift |Δ_q^ac/α_q| で規格化したもの。単なる AWG amplitude や室温 power ではなく ac-Stark shift で較正することで、伝達関数や detuning の違いを吸収している。
  • 赤い線や帯が「その周波数・強度で |1_t> から出る」条件を表す。ローマ字ラベル A-P は代表的な共鳴 feature。緑/マゼンタ/青のラベル色は、著者らの分類、すなわち intrinsic multiphoton、inelastic scattering、TLS/ac-Stark 型を示す。
  • この図を見れば、強ドライブを安全に使える領域と、避けるべき危険領域が直感的に分かる。Qubex に実装する場合の最重要出力は、まさにこのような DUST map と safe mask になる。
実験 map の実用的意味: 赤い feature を避ければよい、という単純な安全地図としても使える。ただし、なぜ赤いのかを分類しないと対策が分からない。TLS なら材料・界面、intrinsic なら周波数割当と Floquet、external mode ならパッケージ/EM/RF 設計が主な対策になる。

3. Mechanism A: ac-Stark shift による TLS 共鳴

Mechanism A は、ドライブが qubit を直接遷移させるというより、ドライブで qubit 周波数が動いた結果、材料欠陥などの TLS と共鳴してエネルギー交換する機構である。DUST map 上では、drive frequency ではなく ac-Stark shift の値に支配されるため、準水平の feature として見える。

Dressed qubit frequency under drive: ω̃_01(ξ, ω_d) = ω_01 + Δ_q^ac(ξ, ω_d) + higher-order shifts Mechanism A condition: ω̃_01(ξ, ω_d) ≈ ω_TLS
Fig. 3

Mechanism A: ac-Stark shift で TLS と共鳴する

Fig. 3: Mechanism A: ac-Stark shift で TLS と共鳴する

図の要点: ac-Stark shift により qubit 遷移周波数が TLS 周波数へ近づき、TLS と resonant exchange を起こす。DUST map では drive frequency よりも Δ_q^ac に支配される quasi-horizontal feature として出る。

読み方:
  • 数式で言えば、ドライブ下の dressed qubit 周波数 ω~_01(ξ,ω_d) が TLS 周波数 ω_TLS に等しくなると、横結合 g によって励起が TLS に移る。条件は ω~_01 = ω_TLS であり、ドライブ周波数そのものの光子数整合ではない。
  • 右下の時間追跡図は TLS bath が静的でないことを示す。長時間観測すると、TLS の周波数は連続的に drift したり、telegraph noise 的に切り替わったりする。したがって TLS 由来の safe/unsafe map は時間依存で、永続 calibration として固定してはいけない。
  • 対策は制御パルスの設計だけでは限界がある。材料、界面、電場集中、TLS density の低減が本質的な mitigation になる。

4. Mechanism B: トランズモン固有の多光子遷移

Mechanism B は、外部モードなしで説明できる。トランズモンは完全な二準位系ではなく、弱い非調和性を持つ多準位非線形振動子である。強い周期ドライブの下では、n 個の drive photons が同時に吸収され、高い非計算準位への共鳴が生じる。

Driven transmon Hamiltonian: H(t) = 4E_C (n - n_g)^2 - E_J cos(φ) + E_d n cos(ω_d t) Intrinsic multiphoton resonance, schematic: ω̃_{i+m,i}(ξ, n_g) ≈ n ω_d Hybridization indicator: Θ(|j_t>) = 1 - | < Floquet_j(ξ,ω_d) | ideal-displaced_j(ξ,ω_d) > |^2
Fig. 4

Mechanism B の局所例: intrinsic multi-photon transition

Fig. 4: Mechanism B の局所例: intrinsic multi-photon transition

図の要点: feature E を例に、トランズモンが複数の drive photons を吸収して非計算準位へ遷移する様子を示す。branch analysis では Floquet branch の swap として共鳴が可視化される。

読み方:
  • この例では、おおまかに ω~_15 ≈ 2ω_d という二光子条件で |1_t> から高準位へ遷移する。縦軸を強度にすると、ac-Stark shift によりエネルギー差が動くため、共鳴線は斜めに曲がる。
  • 下段の branch analysis は、Floquet モードをドライブ強度に沿って追跡し、どの裸状態ラベルの枝が入れ替わるかを見る手法。branch swap は、駆動下固有状態が強く混成し、不要遷移が起こりやすい条件である。
  • 強いドライブ領域では高準位の charge dispersion が大きくなり、offset charge のゆらぎによって共鳴条件が smeared/fuzzy になる。これは実測スペクトルで線が太く不鮮明になる理由の一つである。
Fig. 5

Floquet simulation で intrinsic feature を同定する

Fig. 5: Floquet simulation で intrinsic feature を同定する

図の要点: 実測 DUST map と transmon-only Floquet simulation を並べ、トランズモン固有の多光子遷移だけを拾い出す。simulation は hybridization parameter Θ を用いる。

読み方:
  • 著者らは、裸トランズモン状態と Floquet mode の単純な overlap ではなく、ac-Stark displacement による自明な変形を差し引いた ideal-displaced state を基準にする。その上で Θ(|j_t>) = 1 - |<Floquet_j | ideal-displaced_j>|^2 のような量を計算する。
  • Θ が小さい領域は、強ドライブでも単に状態が dressed/displaced されているだけ。Θ が線状に大きくなる領域は、非計算準位との多光子共鳴を示す。
  • 重要なのは、Fig. 2 の全 feature が Fig. 5 で再現されるわけではない点である。再現されるものは Mechanism B、再現されないものは TLS または外部モードを疑う、という分類手順が立つ。
分類上のポイント: Fig. 5 の transmon-only Floquet simulation と実験 map の feature が重なれば Mechanism B と判定できる。重ならない feature は、TLS または外部モードを疑う。

5. Mechanism C: 外部モードを介した非弾性散乱

Mechanism C は、実験デバイスらしさが最も強く出る機構である。ドライブはトランズモンを励起または緩和させ、その余剰エネルギーを別周波数の photon として環境へ放出する。環境には読み出しモード、キャビティ高次モード、パッケージモード、ケーブル/構造共振、TLS が含まれる。

Schematic energy conservation for inelastic scattering: ω̃_{i+m,i}(ξ) + p ω_s ≈ n ω_d Rate trend, schematic: Γ_C ∝ |matrix element|^2 × Re Z_env(ω_s)
Fig. 6

Mechanism C: spurious EM mode を介した非弾性散乱

Fig. 6: Mechanism C: spurious EM mode を介した非弾性散乱

図の要点: トランズモン遷移と同時に、余剰エネルギーがスプリアス電磁モードへ放出される例。excitation 型と relaxation 型で DUST map 上の slope が逆向きになる。

読み方:
  • 非弾性散乱では、n 個の drive photons が吸収され、トランズモン状態が変わり、余剰エネルギーが ω_s の環境モードへ放出される。概念的には ω~_{i+m,i}(ξ) + pω_s = nω_d のようなエネルギー保存条件で表せる。
  • この過程の強さは、放出周波数で qubit が見る環境インピーダンス Re Z_env(ω_s) に敏感である。パッケージ、読み出しキャビティ、高次モード、ケーブルや構造共振がピークを持つと、DUST が増強される。
  • 設計上は、qubit/readout 周波数だけを見ればよいわけではない。数十 GHz 近くまでの RF/EM モードが、強ドライブ下では間接的に制御 fidelity を制限し得る。
Fig. 7

transmon-only では見えない高次 inelastic process

Fig. 7: transmon-only では見えない高次 inelastic process

図の要点: feature K, L は実験に現れるが transmon-only Floquet には出ない。読み出し resonator を含めた二モード Floquet simulation では再現される。

読み方:
  • この図は論文の説得力を高める重要図である。実験 map には K, L が見えるが、トランズモンだけの Hamiltonian では現れない。したがって、これは「トランズモン固有の多光子遷移」ではない。
  • 読み出し resonator を明示的に Hilbert space に入れると、K/L が再現される。これは、読み出しモードが単なる readout channel ではなく、強ドライブ下では DUST を媒介する動的自由度になり得ることを示す。
  • ただし全ての外部モードを Floquet に入れるのはスケーラブルではない。現実的には、少数モード Floquet、EM simulation、実測 DUST map を組み合わせる必要がある。
Fig. 8

TLS-assisted inelastic scattering: 線が切れる理由

Fig. 8: TLS-assisted inelastic scattering: 線が切れる理由

図の要点: 外部モードが安定した EM mode ではなく、周波数を switching する TLS の場合、DUST feature は連続線ではなく broken line として現れる。

読み方:
  • Fig. 3 の TLS は qubit 近傍の周波数で resonant exchange を起こしたが、Fig. 8 では TLS が遠方周波数の散乱先として働いている。つまり TLS は Mechanism A だけでなく Mechanism C にも関与する。
  • TLS 周波数が二値的に切り替わると、エネルギー保存条件を満たす drive frequency/power の組も切り替わる。そのため、実測スペクトルでは線が途中で途切れたり、二本の近い線になったりする。
  • Qubex の自動 feature extraction では、こうした broken/telegraphic feature を単純な連続線検出だけで見落とさない工夫が必要になる。

6. Appendix 図で補強される実験・解析の詳細

Appendix の図は、本文の主張を支える実験技術と追加検証を示す。DUST の実装や再現を考える場合は、むしろ appendix 図の方が具体的な設計情報を与える。

Fig. A1

多準位 single-shot readout の分離

Fig. A1: 多準位 single-shot readout の分離

図の要点: |0_t> から |4_t> までの readout clusters。DUST 後にどの準位へ人口が移ったかを判定する基礎になる。

読み方:
  • DUST map は、単に「読めた/読めない」の図ではない。最終状態の人口 P(0_t), P(1_t), P(2_t), P(rest) を分類できる readout が必要である。
  • readout 自体が DUST を起こすと pump-probe の解釈が壊れるため、著者らは readout power を制限し、刺激ドライブ由来の遷移だけを測るようにしている。
Fig. A3

初期 |0_t>, |1_t> の full DUST landscape

Fig. A3: 初期 |0_t>, |1_t> の full DUST landscape

図の要点: 実験 map を single-quantum transitions と higher-state leakage に分け、さらに transmon-only Floquet simulation と比較した appendix の総覧図。

読み方:
  • Fig. 2 は主に |1_t> 初期化の代表 map だが、Appendix では |0_t> と |1_t> の両方を interleaved に測っている。左右の対応を見ると、同じ外部モードや同じ非線形過程が、初期状態の違いに応じて少し異なる周波数に現れる。
  • 上段の single-quantum transitions は Mechanism A/C が主に寄与し、下段の higher-state leakage には Mechanism B/C が強く現れる。下の simulation と照合して、green labels は intrinsic と判定される。
Fig. A4

Floquet model に使う lumped circuit

Fig. A4: Floquet model に使う lumped circuit

図の要点: transmon-only model と readout resonator を含む model の境界を示す回路図。外部モードを含めると Hilbert space が増えるが、Mechanism C の再現には必要になる場合がある。

読み方:
  • トランズモン単体の Hamiltonian では Mechanism B を扱える。一方、readout mode や spurious mode を明示的に含むと、Mechanism C の一部を Floquet simulation で予測できる。
  • Qubex 実装では、最初から全ての環境モードを入れるより、transmon-only を標準、selected-mode simulation を optional にする設計が現実的である。
Fig. A5

readout mode を含む Floquet branch analysis

Fig. A5: readout mode を含む Floquet branch analysis

図の要点: joint transmon-resonator Hilbert space での hybridization map と branch analysis。Fig. 7 の K/L のような feature を説明するための appendix 解析。

読み方:
  • readout resonator を入れると状態ラベルは |transmon, resonator> の積状態に近いものになる。ただし高励起では強く混成し、明確な bare label が付けられない dressed state も出る。
  • この図は、外部モードを入れれば説明力が増す一方で、状態数とモード数の増加により計算が重くなることも示している。
Fig. A6

パッケージの spurious EM modes と実験検証

Fig. A6: パッケージの spurious EM modes と実験検証

図の要点: HFSS などの EM simulation で得た environmental impedance と、Ramsey 型 spectroscopy による EM1/EM2 の実験確認。

読み方:
  • 理想対称な 3D cavity では transmon と結合しにくいモードも、チップ位置ずれや coupling pin による対称性破れで結合する。したがって layout simulation では完全対称モデルだけを見ると危険である。
  • DUST feature A, D, G, I, J, N などは、こうした simulated/observed EM modes と対応付けられる。DUST 対策には「広帯域での microwave hygiene」が必要になる。
Fig. A7

TLS との resonant exchange の時間ダイナミクス

Fig. A7: TLS との resonant exchange の時間ダイナミクス

図の要点: ac-Stark tone で qubit を TLS に park し、保持時間を掃引した時間領域測定。強く散逸する TLS では指数減衰、弱散逸 TLS では減衰振動が現れる。

読み方:
  • Mechanism A は静的スペクトル上の horizontal feature として見るだけでなく、時間領域で resonant exchange rate や TLS dissipation を調べられる。
  • この測定は、DUST map の feature が本当に TLS exchange かどうかを確認する診断にもなる。
Fig. A8

intrinsic multi-photon transition の時間ダイナミクス

Fig. A8: intrinsic multi-photon transition の時間ダイナミクス

図の要点: feature E′ に drive 条件を合わせ、|0_t> と高準位間の人口振動を観測する。Floquet branch swap と対応する。

読み方:
  • Mechanism B は単なる損失ではなく、条件がきれいに合えば coherent population oscillation として見える。これは多光子遷移が駆動 Hamiltonian 内の実在の avoided crossing であることを示す。
  • ただし高パワーでは複数の高準位や offset charge 依存が絡み、単一の Rabi-like rate として取り出すのは難しくなる。
Fig. A9

switching TLS による extrinsic excitation の確認

Fig. A9: switching TLS による extrinsic excitation の確認

図の要点: feature B/B′ の拡大と single-shot time trace。TLS が on/off 的に切り替わると、同じ drive 条件でも excitation が出る時間窓と出ない時間窓が生じる。

読み方:
  • Fig. 8 の broken line を時間領域で裏付ける appendix 図である。single-shot trace を見ると、transition が統計平均ではなく時間的にスイッチングしていることが分かる。
  • DUST 診断では、feature の平均強度だけでなく、時間変動性も保存しておく価値がある。Qubex なら raw shots や time stamps を metadata として残す設計が望ましい。

7. 数理モデルのまとめ

7.1 測定量

実験の raw output は、初期状態 i、drive frequency ω_d、drive amplitude A_d に対する最終人口 P(k | i; ω_d,A_d) である。DUST 解析ではこれを transition probability、leakage probability、state-resolved population map に変換する。

P_leave,i = 1 - P(i | i; ω_d,A_d) P_relax,e→g = P(g | e; ω_d,A_d) P_excite,g→e/f/rest = 1 - P(g | g; ω_d,A_d) P_rest = 1 - P(g) - P(e) - P(f) # 実験上の分類不能または高準位人口

7.2 ac-Stark 軸

drive amplitude だけでは、周波数ごとの伝達関数差を吸収できない。そこで、低パワーで測った qubit resonance shift から A_d → Δ_q^ac を較正し、全 map を |Δ_q^ac/α_q| の軸で表示する。

A_d^2 at room temperature ─calibration→ Δ_q^ac(ω_d, A_d) normalized drive strength: s = |Δ_q^ac / α_q|

7.3 分類規則

観測上の特徴候補機構追加確認主な対策
ac-Stark shift 軸でほぼ水平。時間的に drift/switching することがある。Mechanism A: TLS/ac-StarkTLS parking / time-domain exchange材料・界面改善、電場集中低減、safe amplitude band 回避
transmon-only Floquet の Θ map に出る。Mechanism B: intrinsic multi-photonbranch analysis、時間領域 population oscillationdrive frequency allocation、charge dispersion 低減、パルス上限設定
transmon-only では出ないが、EM mode/readout mode/TLS を含めると説明できる。Mechanism C: external-mode inelastic scatteringEM simulation、Ramsey-based mode spectroscopy、selected-mode Floquetpackage/EM 設計、広帯域フィルタ、spurious mode detuning

8. 診断 workflow と Qubex への落とし込み

Qubex に組み込むなら、DUST は「強ドライブ安全性診断」として実装するのが自然である。出力は論文の Fig. 2/A3 に対応する map と、その map から得られる feature table、safe mask、recommended windows である。

取得するデータ

initial_state × drive_frequency × drive_amplitude の sweep を行い、state-resolved population を保存する。

result.populations[initial][final][f_idx, a_idx] result.metrics.P_leave[initial] result.metrics.P_leakage[initial]

返すべき出力

DUST heatmap、population maps、feature table、safe/unsafe mask、recommended windows、calibration guardrail。

safe_mask = (P_leave < threshold_leave) AND (P_leak < threshold_leak) feature = line/patch in map with severity, mechanism_candidate, action

擬似 API

dust = qx.contrib.dust.scan( exp, target="Q03", drive_frequencies=np.linspace(3.0e9, 11.8e9, 401), drive_amplitudes=np.linspace(0.0, 0.35, 121), initial_states=("g", "e"), pump_duration_ns=1000, post_pump_wait_ns=400, classify_states=("g", "e", "f", "rest"), ) dust.plot_transition_map(initial="e") dust.plot_population_map(initial="e", final="f") dust.safe_mask(max_p_leave=0.01, max_p_leakage=0.002) dust.recommended_windows() dust.features.to_dataframe()
実験価値: readout amplitude、CR/Stark drive、reset pulse、qudit pulse の calibration は、性能最大化だけでなく「DUST が出ない範囲で性能最大化」へ変えられる。

9. 限界と注意点

  • TLS は時間変動する。 Fig. 3 と Fig. A9 が示すように、TLS feature は drift/switching する。DUST map は定期更新が必要である。
  • Floquet の Θ は rate そのものではない。 Θ は共鳴条件を見つける強力な指標だが、実験 transition probability との定量対応は論文でも今後の課題として残る。
  • 全外部モード Floquet は重い。 モード数が増えると Hilbert space が爆発する。実用上は EM simulation で重要モードを絞り、selected-mode simulation を行う。
  • 広帯域強ドライブ scan は装置リスクがある。 mixer saturation、ADC saturation、加熱、quasiparticle generation、長寿命 leakage を監視し、abort 条件を必ず入れるべきである。

10. 出典・ライセンス

画像は ar5iv による arXiv:2506.24070 の HTML レンダリングから取得した図を、この解説レポート内に埋め込んだもの。本文 Fig. 1-8 はすべて挿入し、Appendix は公開レンダリング上で取得できた主要補助図を v3 の番号に合わせて A1, A3-A9 として掲載している。現行 arXiv v3 は 18 figures だが、アップロード PDF 側が期限切れで PDF からの全図直接抽出ができなかったため、取得可能な図を優先した。本文説明は論文内容に基づく要約・解説であり、キャプションは原文の逐語転載ではなく説明用に再構成している。

論文: W. Dai et al., “Characterization of drive-induced unwanted state transitions in superconducting circuits,” arXiv:2506.24070 / Physical Review X. arXiv ページ上の license link は Creative Commons Attribution 4.0 International (CC BY 4.0) を指している。図を再利用する場合は、著者・論文・arXiv/DOI・CC BY 4.0 への attribution を保持すること。

この HTML は自己完結型で、図は base64 として埋め込まれている。