DUST 論文 詳細レポート

W. Dai ら “Characterization of Drive-Induced Unwanted State Transitions in Superconducting Circuits” の背景、数理、実験、結果、設計上の示唆を、超伝導量子回路の制御・キャリブレーション実装に接続して整理する。

主題
DUST / 強駆動下の不要遷移
対象系
固定周波数 3D transmon
手法
pump-probe + Floquet + EM simulation
作成日
2026-05-26

1. エグゼクティブサマリー

この論文は、超伝導量子回路でマイクロ波ドライブを強くしたときに現れる DUST、すなわち Drive-Induced Unwanted State Transitions を、単なる経験的な「危険なパワー上限」ではなく、物理機構ごとに分類・予測・回避するための実験的・理論的枠組みとして提示している。

一文で言うと: 高速ゲートや高速読み出しのために強いマイクロ波を使うと、量子ビットが計算空間内外へ意図せず遷移する。この論文は、それを A: ac-Stark shift による TLS 共鳴B: トランズモン固有の多光子共鳴C: 外部モードを介した非弾性散乱 に分解し、pump-probe spectroscopy と Floquet/EM simulation で見分ける方法を示す。

研究上の価値

DUST を「測定で見える謎のリーク」から「分類可能なスペクトル特徴」へ変換した点が重要。強駆動制御の設計段階で避けるべき周波数・振幅領域を同定できる。

実験上の価値

固定周波数 transmon でも ac-Stark shift をノブとして使い、TLS や外部モード起因の遷移を時間分解 pump-probe spectroscopy で診断できる。

実装上の価値

Qubex のような制御フレームワークでは、DUST map を readout / gate / Stark drive calibration の 安全制約 として利用できる。

2. 書誌情報と位置づけ

論文名Characterization of Drive-Induced Unwanted State Transitions in Superconducting Circuits
著者W. Dai, S. Hazra, D. K. Weiss, P. D. Kurilovich, T. Connolly, H. K. Babla, S. Singh, V. R. Joshi, A. Z. Ding, P. D. Parakh, J. Venkatraman, X. Xiao, L. Frunzio, M. H. Devoret
掲載Physical Review X 16, 011011 (2026)
DOI10.1103/zdpg-mhpc
公開2026-01-15 published; arXiv:2506.24070 は 2025-06-30 初稿、2025-12-22 v3。
ライセンス公開ページ上では CC BY 4.0 として扱われている。

arXiv 版の abstract は、マイクロ波ドライブが超伝導量子回路の制御・読み出しに不可欠である一方、ドライブ強度を上げると不要遷移が起こり、操作速度と fidelity を制限すると述べている。論文は固定周波数 qubit に対して約 9 GHz 幅のドライブ周波数を掃引し、DUST の起源を3種類に分類する。

この論文の位置づけ: MIST、ionization、readout-induced transitions、drive-induced leakage といった関連現象を、より一般的な 強駆動下の不要遷移スペクトロスコピー として捉え直している。対象は 3D transmon だが、著者らは fluxonium、multimodal circuit、SNAIL、coupler など他の超伝導非線形回路にも適用可能だと位置づけている。

3. 問題設定:なぜ DUST が本質的な制限になるのか

超伝導量子回路の操作は、マイクロ波ドライブを通じて実行される。単一量子ビットゲート、二量子ビットゲート、読み出し、リセット、ボソニックモード制御、reservoir engineering など、ほぼすべての実験的操作に drive tone が関与する。

操作時間を短くするには、Rabi rate や readout photon number を上げる必要がある。直感的には、強いドライブを使えば高速になる。しかし強いドライブは、トランズモンの非線形性、計算空間外準位、TLS、読み出し共振器、パッケージモード、配線モード、材料欠陥を同時に巻き込む。そのため、ゲートや読み出しの fidelity は単純な decoherence time だけでなく、強駆動時にどの不要遷移が開くかに制限される。

高速操作したい gate / readout / reset ドライブを強くする 周波数 × 振幅を探索 不要遷移が開く relaxation / excitation / leakage fidelity / QND 性 / reset / QEC cycle を劣化 見かけの SNR や Rabi rate だけでは安全性を判定できない
図1:DUST の工学的問題。高速化のための強ドライブが、計算空間外準位や環境自由度を巻き込み、性能指標を壊す。

この論文の問いは、「実験で観測される多数の不要遷移を、どの物理機構に由来するか切り分けられるか」である。これに答えられれば、対策も変わる。TLS なら材料・界面・電場集中の問題、多光子共鳴なら周波数割当・パルス設計の問題、外部モード非弾性散乱ならパッケージ・RF hygiene・広帯域フィルタリングの問題になる。

4. 実験系と測定プロトコル

実験は、3D rectangular waveguide cavity 中の固定周波数 transmon で行われる。キャビティの fundamental mode は readout resonator として使われ、pump は専用 drive port から入る。入力 drive line は広い周波数範囲で十分な power delivery ができるように設計される一方、qubit coherence を壊さないよう cryogenic attenuation も確保される。

主要パラメータ

意味
qubit 周波数ωq/2π = 4.5285 GHz|0t⟩ ↔ |1t⟩ 遷移
非調和性αq/2π = -184.2 MHztransmon の弱い負の anharmonicity
readout resonator 周波数ωr/2π = 9.0342 GHzreadout mode; spurious/inelastic processes にも関わる
resonator linewidthκ/2π = 7.20 MHzreadout photon decay の指標
分散シフトχ/2π = -1.23 MHzmulti-state single-shot readout 用
pump 周波数範囲およそ 2.8-11.8 GHz約 9 GHz 幅の broadband spectroscopy
pump 長1 μsDUST を励起・観測する stimulation pulse
pump 後 wait400 ns残留 cavity photon を抜く
coherenceT1 ≈ 85 μs, T2E ≈ 135 μs測定期間中の代表値

pump-probe spectroscopy

基本的な実験系列は次のとおり。

  1. transmon を |0t⟩ または |1t⟩ に初期化する。
  2. 周波数 ωd と power/amplitude を変えた pump pulse を 1 μs 入れる。
  3. 400 ns 待って残留 photon を減らす。
  4. single-shot readout で最終状態 population を測る。
  5. P(0t), P(1t), P(2t), P(rest) を集計し、transition probability map を作る。
P(rest) = 1 - P(0_t) - P(1_t) - P(2_t) 初期状態 |s_t> から出た確率: P_leave,s(ω_d, A_d) = 1 - P(s_t | initial = s_t, ω_d, A_d)

ac-Stark shift による power calibration

この論文で特に重要なのは、横軸を単なる室温発生器 power ではなく、qubit に実際に誘起された ac-Stark shift で較正している点である。同じ room-temperature amplitude でも、周波数依存の transfer function、drive line、mixer、filter、readout resonator 近傍の応答によって、transmon へ届く有効強度は変わる。

低 power 領域で各 ω_d ごとに近似: Δ_q^ac(ω_d, A_d) ≈ s(ω_d) · A_d^2 ここで s(ω_d) は測定から得る周波数依存 calibration slope。

そのため DUST map は、実装上も drive frequency × AWG amplitude だけでなく、drive frequency × estimated ac-Stark shift として保存・表示するべきである。

5. DUST の数理モデル

DUST は単一の現象名だが、数理的には複数の機構の総称である。共通する骨格は、周期駆動された非線形量子系が、裸の計算基底だけでは記述できない dressed/Floquet 状態へ移ることである。

5.1 駆動トランズモン Hamiltonian

transmon の最小モデルは、Josephson cosine potential と charging energy からなる。

H_t / ℏ = 4E_C( n̂ - n_g )^2 / ℏ - E_J cos(φ̂) / ℏ 駆動項を入れると: H(t) = H_t + H_d(t) H_d(t) ≈ ε_d n̂ cos(ω_d t) E_C: charging energy E_J: Josephson energy n_g: offset charge n̂, φ̂: charge and phase operators

ドライブの存在により、transmon eigenstate |jt⟩ は単純な裸状態ではなく、ドライブで変位・混成した dressed state になる。弱駆動では ac-Stark shift として扱えるが、強駆動では多光子共鳴・branch swap・非摂動効果が重要になる。

5.2 ac-Stark shifted transition

ドライブ下の遷移周波数を、概念的に次のように書く。

ω̃_ij(ξ, ω_d) = [Ẽ_i(ξ, ω_d) - Ẽ_j(ξ, ω_d)] / ℏ 特に qubit 遷移では: ω̃_q(ξ, ω_d) = ω_q + Δ_q^ac(ξ, ω_d)

ここで ξ は drive strength を表す無次元量と考えればよい。実験では ξ を直接制御するのではなく、各 drive frequency ごとに測った Δqac によって drive strength を再パラメータ化する。

5.3 Floquet 表現

ドライブ Hamiltonian が周期 T = 2π/ωd を持つとき、Floquet theorem により解は次の形で表せる。

|ψ_α(t)> = exp(-i ε_α t / ℏ) |u_α(t)> |u_α(t + T)> = |u_α(t)> ε_α: quasienergy |u_α(t)>: Floquet mode

多光子共鳴は、drive amplitude を変えたときに Floquet branch が避け交差または branch swap を起こす点として現れる。論文では、この branch analysis と hybridization parameter を使って intrinsic multi-photon transitions を検出する。

5.4 hybridization parameter

強いドライブでは、単に「裸状態との overlap が小さい」だけでは不要遷移と判定できない。ac-Stark shift や coherent displacement だけでも overlap は変わるからである。そこで論文は、理想的に displaced された状態を基準にして、そこからの逸脱を共鳴の指標にする。

Θ_j(ξ, ω_d) = 1 - | < ˜j_t(ξ,ω_d) | ȷ̄_t(ξ,ω_d) > |^2 ˜j_t: 実際の Floquet mode ȷ̄_t: ac-Stark displacement だけを含む ideal-displaced state Θ_j ≈ 0 : 単なる Stark/displacement Θ_j が急増 : 非計算準位との多光子混成候補

5.5 DUST transition probability

実験的には、理論的な branch swap そのものではなく、pump 後の population 変化を測る。したがって DUST map の基礎量は population である。

初期状態 |0_t> の DUST: P_DUST,0 = 1 - P(0_t) 初期状態 |1_t> の DUST: P_DUST,1 = 1 - P(1_t) leakage 指標: P_leak ≈ P(2_t) + P(rest) ※定義は用途に応じて調整

QEC や readout validation では、単に Pleave だけでなく、Prest、Pf、再初期化不能な長寿命 population を分離して見る必要がある。

6. 3つの DUST 機構

A. TLS 共鳴 ac-Stark shift が qubit を動かす q TLS 横線状: Δ_ac のみで決まる B. 多光子共鳴 transmon Hamiltonian 固有 |0/1> |k> m photons Floquet branch swap で検出 C. 非弾性散乱 外部モードへ余剰エネルギー放出 q EM Re Z_env(ω_s) のピークで増強
図2:論文の DUST 分類。A は ac-Stark shift による TLS との共鳴、B は transmon 固有の多光子共鳴、C は外部モードへエネルギーを逃がす非弾性散乱。
機構物理像スペクトル上の特徴主な対策
A ac-Stark/TLS ドライブで qubit 周波数が Stark shift し、TLS と共鳴してエネルギー交換する。 drive frequency そのものには弱く依存し、Δqac 軸でほぼ横線になる。 材料改善、界面損失低減、電場集中低減、TLS monitoring。
B intrinsic multiphoton トランズモンの非線形 Hamiltonian により、複数 drive photon を吸収して高準位へ遷移する。 ωd と power の両方に依存する斜め線。Floquet simulation で再現される。 drive frequency allocation、Floquet screening、offset charge sensitivity 低減。
C external-mode inelastic scattering drive が transmon 遷移を誘起し、余剰エネルギーを readout/package/TLS など外部モードへ放出する。 斜め線。transmon-only Floquet では出ず、外部モードや EM simulation と対応する。 RF package design、spurious mode engineering、broadband filtering、EM co-simulation。

A. ac-Stark-shift-induced resonant exchange with TLS

最小モデルでは、駆動でシフトした qubit と TLS の transverse coupling を考える。

H_A / ℏ = ½ ω̃_q(ξ) σ_z + ½ ω_TLS τ_z + g_TLS (σ_+ τ_- + σ_- τ_+) 共鳴条件: ω̃_q(ξ) = ω_TLS 冷たく散逸的な TLS なら、qubit excitation が TLS へ移り、最終的に |1_t> → |0_t> の relaxation として見える。

この機構の判定基準は明快である。共鳴位置が drive frequency ではなく、誘起された ac-Stark shift の値で決まるなら TLS/ac-Stark 型候補と見る。論文では、TLS の周波数 drift、telegraph-noise-like switching、数時間スケールの出現・消失が観測されている。

B. Intrinsic multi-photon transitions

transmon の cosine potential は弱非線形なので、強い drive 下では多光子過程が開く。代表的な共鳴条件は概念的に次の形になる。

m ω_d ≈ ω̃_ij(ξ, ω_d) m: 吸収される drive photon 数 i,j: transmon の裸状態または dressed state labels

低 power では摂動論的な resonance line で見通せるが、高 power では複数の高準位、non-RWA 効果、offset charge drift、branch swap が絡むため、Floquet 解析が必要になる。

C. Inelastic scattering involving external modes

外部モードがあると、drive photon のエネルギーは transmon 遷移と外部 photon へ分配される。

概念的なエネルギー保存: ω_s ≈ | m ω_d - ω̃_ij(ξ, ω_d) | 遷移率の形: Γ_C ∝ |M_m(ξ)|^2 · Re Z_env(ω_s) ω_s: 散乱 photon の周波数 Z_env: transmon が見る環境インピーダンス M_m: m 光子 mixing matrix element

環境インピーダンスの実部にピークがある、つまり readout mode、package mode、spurious EM mode、TLS defect が存在する周波数では、この過程が増強される。論文では、EM simulation で得られる package mode が複数の DUST feature を説明することを示している。

7. 主要結果の読み解き

7.1 DUST map は「危険領域の地図」

Fig. 2 相当の DUST map では、x 軸に drive frequency、y 軸に drive power または normalized ac-Stark shift、色に transition probability を置く。すると、drive condition の平面上に線状・帯状の resonance feature が現れる。

この map から読み取るべき点は、単に「高 power が危険」ではないことである。ある周波数では比較的高 power まで安全でも、別の周波数では低 power で鋭い不要遷移が出る。したがって、readout や gate の drive 周波数は、SNR や Rabi rate だけでなく DUST map 上の safe window を考慮して選ぶ必要がある。

7.2 quasi-horizontal feature は TLS/ac-Stark 型

O, P のような特徴は、drive frequency に対してほぼ水平で、ac-Stark shift の特定値で起こる。これは qubit の有効周波数が Stark shift によって TLS と共鳴したときに発生するため、drive photon energy の整数倍条件ではなく、qubit-TLS detuning が支配的である。

7.3 Floquet simulation は intrinsic feature を選別する

transmon-only Hamiltonian の Floquet simulation で再現される feature は、外部モードを仮定しなくても起こる intrinsic multiphoton resonance と判定される。論文では hybridization parameter と branch analysis により、実験で観測された複数の feature が transmon 固有の多光子遷移として同定される。

注意: Floquet simulation の hybridization parameter は、観測される transition probability と定量的に一対一対応するわけではない。論文自身も、Θ と実測 transition probability の定量接続は今後の課題としている。したがって実装では、Floquet は「危険候補の位置予測」として扱い、実測で validation するのが安全である。

7.4 external mode を入れないと説明できない feature がある

transmon-only Floquet では現れない K, L のような feature は、readout mode や spurious EM mode を含めた拡張 Hamiltonian で説明される。この結果は、DUST が単なる transmon の非線形問題ではなく、デバイス全体の RF 環境問題でもあることを示す。

特に重要なのは、動作帯域から離れた package mode でも inelastic scattering の受け皿になり得ることである。したがって、設計時には qubit/readout 周波数近傍だけでなく、drive によって散乱 photon が出得る広い周波数範囲の EM mode を評価する必要がある。

7.5 TLS は qubit 近傍だけでなく高周波 mixing にも関わる

Sec. VIII では、遠く detune した TLS が inelastic scattering の外部モードとして関与し、broken line や telegraphic switching のような特徴を作る例が示される。これは、TLS が単に qubit 周波数近傍にある場合だけでなく、多光子 mixing を通じて広い周波数範囲で DUST に寄与し得ることを示している。

8. 論文から抽出できる DUST 分類アルゴリズム

論文の判断手順は、実験フレームワークに落とし込むと次のようになる。

pump-probe DUST map feature extraction 線状特徴・閾値超え・人口移動パターンを検出 Δ_ac のみで決まるか? quasi-horizontal feature? Yes A TLS/ac-Stark No transmon-only Floquet で出るか? hybridization / branch swap Yes B intrinsic No C external mode
図3:論文の判定ロジックを実装向けに整理したもの。A → B → C の順に切り分ける。

擬似コード

for feature in detected_features(dust_map):
    if feature.position_depends_only_on_ac_stark_shift():
        feature.mechanism = "A: TLS/ac-Stark resonant exchange"
    elif transmon_only_floquet.predicts(feature):
        feature.mechanism = "B: intrinsic multiphoton transition"
    else:
        feature.mechanism = "C: inelastic scattering via external mode"
        feature.candidate_modes = match_to_em_modes(feature, em_mode_list)

この順序が重要である。A は drive frequency matching ではなく ac-Stark shift matching で起こるため、まず横線状 feature として分離できる。B は理想 transmon Hamiltonian に含まれるため、transmon-only Floquet で検出できる。残りが C だが、C の中にも stable EM mode と switching TLS defect が混在する。

9. 設計・実装への示唆

9.1 frequency allocation は DUST map を見て行うべき

readout frequency、CR drive frequency、Stark tone、parametric drive、reset pulse などの周波数は、従来は desired interaction rate、detuning、crosstalk、SNR などを見て選ばれることが多い。しかし DUST の観点では、それらに加えて DUST-free window を選ぶ必要がある。

実験指針: 「最も速い drive 条件」ではなく、「不要遷移確率を閾値以下に保ったまま最も速い drive 条件」を選ぶ。これは gate/readout optimizer の目的関数を、性能最大化から制約付き最適化へ変えることを意味する。

9.2 readout は assignment fidelity だけでは不十分

高 power readout では、IQ separation や assignment fidelity が良く見えても、測定後に transmon が |2⟩ や higher/rest states へ漏れている可能性がある。したがって readout calibration には、測定中・測定後の DUST/MIST leakage validation を含めるべきである。

9.3 package/RF 設計は広帯域で評価する

機構 C は、散乱 photon 周波数における Re Zenv に敏感である。したがって、qubit 周波数と readout 周波数だけでなく、mixing により到達可能な広い周波数帯に対して spurious mode を確認する必要がある。論文では複数の spurious electromagnetic modes が feature を説明し、遠い RF mode も DUST に寄与し得ることが示されている。

9.4 TLS は「材料問題」と「診断信号」の両方である

TLS は DUST の原因としては厄介だが、論文の pump-probe method は固定周波数 transmon を使った TLS spectroscopy としても機能する。flux-tunable qubit なしに、ac-Stark shift を使って TLS bath の drift や switching を観測できる点は、材料評価ツールとしても意味がある。

10. Qubex 実装への接続

本レポートの主対象は論文そのものだが、これまでの議論に沿って Qubex へ組み込む場合の見取り図も示す。DUST は gate primitive ではなく、強ドライブ条件の診断・制約生成モジュールとして実装するのが自然である。

10.1 主要出力

出力内容用途
DUST heatmapPleaved, Ad) または Pleaved, Δac)危険領域の可視化
population mapsPg, Pe, Pf, Prestrelaxation, excitation, leakage の分離
feature tableresonance feature の位置、強度、候補機構物理原因の分類と再測定対象の特定
safe mask閾値以下の drive condition を True にする boolean mapcalibration optimizer の探索制約
recommended windows安全な周波数範囲と最大 amplitudereadout/gate/Stark tone 設計
guardrailDUST score を返す callable constraint自動 calibration に組み込む

10.2 DUST score の例

S_DUST(ω_d, A_d) = w1 · P_leave + w2 · P_leak + w3 · P_rest + w4 · proximity_to_detected_feature S_DUST < 0.01 safe 0.01-0.05 caution 0.05-0.20 risky > 0.20 forbidden
DUST score は絶対的な物理量ではなく、optimizer 用の工学的 risk metric として設計する。

10.3 API イメージ

dust = qx.contrib.dust.scan(
    exp,
    target="Q03",
    drive_frequencies=np.linspace(3.0e9, 8.0e9, 401),
    drive_amplitudes=np.linspace(0.0, 0.35, 121),
    initial_states=("g", "e"),
    pump_duration_ns=1000,
    post_pump_wait_ns=400,
    n_shots=1024,
)

dust.plot_transition_map(initial="e")
dust.plot_leakage_map(initial="g")
safe = dust.safe_mask(max_p_leave=0.01, max_p_leakage=0.002)
windows = dust.recommended_windows()

guardrail = dust.as_guardrail(max_p_leave=0.01, max_p_leakage=0.002)

10.4 実装方針

  1. Phase 1: pump-probe 2D scan と DUST heatmap を実装する。
  2. Phase 2: drive amplitude を ac-Stark shift 軸へ較正する。
  3. Phase 3: feature extraction と簡易ラベル付けを行う。
  4. Phase 4: safe mask / guardrail を readout/gate calibration に接続する。
  5. Phase 5: optional dependency として Floquet prediction を追加する。
  6. Phase 6: 外部 EM mode list を読み込み、機構 C の候補 mode matching を行う。

11. 限界と未解決課題

定量予測の限界

Floquet hybridization parameter は feature の位置予測には有用だが、transition probability を直接与えるわけではない。実測 validation が必要。

多モード Floquet の計算コスト

外部モードをすべて入れると Hilbert space が急増し、計算が現実的でなくなる。効率的な reduced model が必要。

TLS の予測不能性

TLS は個体差、時間 drift、telegraphic switching を示す。設計段階で完全予測するのは難しく、定期的な spectroscopy が必要。

quasiparticle など未探索要因

論文は strong drive による quasiparticle generation や quasiparticle-assisted transitions を未探索の重要要因として挙げている。

実験実装時の注意: DUST scan は強ドライブを意図的にかける測定である。装置の mixer saturation、room-temperature chain の非線形、spurious harmonics、ADC saturation、qubit heating、長寿命 leakage を監視し、abort 条件を必ず設けるべきである。

12. 読み方ガイド

  1. Abstract と Fig. 1: 3つの機構 A/B/C の全体像を把握する。
  2. Sec. II と Fig. 2: pump-probe DUST map の読み方を理解する。
  3. Sec. III: ac-Stark shift をノブにした TLS spectroscopy と quasi-horizontal feature を確認する。
  4. Sec. IV-V: intrinsic multiphoton resonance、Floquet branch analysis、hybridization parameter を読む。
  5. Sec. VI-VIII: external mode、readout mode、TLS-assisted inelastic scattering の例を追う。
  6. Sec. IX: 設計上の mitigation strategy と未解決課題を確認する。
  7. Appendix: 実装する場合は power calibration、readout threshold、state labeling、EM mode identification を重点的に読む。

13. 出典・参照情報

  1. APS / Physical Review X: W. Dai et al., Characterization of Drive-Induced Unwanted State Transitions in Superconducting Circuits, Phys. Rev. X 16, 011011 (2026), DOI: 10.1103/zdpg-mhpc.
  2. arXiv abstract page: arXiv:2506.24070, Characterization of drive-induced unwanted state transitions in superconducting circuits. Submission history: v1 2025-06-30, v2 2025-08-02, v3 2025-12-22.
  3. arXiv HTML full text: Spectroscopy/Characterization of drive-induced unwanted state transitions in superconducting circuits, sections I-IX and appendices.
  4. Associated code/data: paper reference to Zenodo data/code package, DOI 10.5281/zenodo.18011745, and open-source floquet package for reproducing simulation results.