超伝導量子ビットの物理・数理：全体の地図
この資料は、transmon/circuit-QED 系を軸に、単一量子ビット、1Q–1R、2Q、CR gate、キャリブレーション、制御装置までを「何を近似したか」「どの式変形をしたか」が追える形でまとめたものです。各ページの Interactive lab でパラメータを動かすと、近似式の依存性が視覚的に確認できます。
最短の概念ルート
LC 回路
→
Josephson 非線形
$-E_J\cos\varphi$
$-E_J\cos\varphi$
→
弱非調和人工原子
Duffing oscillator
Duffing oscillator
→
microwave drive
$X/Y/Z$ gate
$X/Y/Z$ gate
→
resonator 読み出し
dispersive shift $\chi$
dispersive shift $\chi$
→
2Q interaction
$XX+YY$, $ZZ$, CR $ZX$
$XX+YY$, $ZZ$, CR $ZX$
基本式の骨格：
1. Transmon: H = 4E_C(n-n_g)^2 - E_J cos φ
2. Duffing: H/ℏ = ω_q b†b + (α/2)b†b(b†b-1)
3. 1Q drive: H_rot/ℏ = -Δ Z/2 + Ω_x X/2 + Ω_y Y/2
4. JC: H/ℏ = ω_r a†a + ω_q Z/2 + g(a†σ_- + aσ_+)
5. 2Q: H/ℏ = Σ_i[ω_i n_i + α_i n_i(n_i-1)/2] + J(b_1†b_2+b_1b_2†)
6. CR: H_CR/ℏ = 1/2 Σ_{μν} ω_{μν} σ_μ⊗σ_ν, 欲しい主成分は ZX
この資料の読み方
- まず各テーマの「物理の絵」を見る。
- 次に「導出」を開き、変数定義・座標変換・摂動展開・捨てた項を確認する。
- Interactive lab のスライダーで、近似がいつ破綻するかを見る。
- 最後に calibration ページで、理論量がどの実験シーケンスで測られるかを対応づける。
「省略なし」を実務的に満たすため、各セクションに 近似ログ を置いています。完全な第一原理シミュレーションでは、ここで示す有効ハミルトニアンではなく、チップの full circuit model、寄生モード、線路歪み、開放系、測定 chain を含む数値モデルへ進みます。
記号規約
| 記号 | 意味 |
|---|---|
| $E_J$ | Josephson energy。接合の非線形インダクタンスを決める。 |
| $E_C=e^2/2C_\Sigma$ | charging energy。transmon の非調和性をほぼ決める。 |
| $\alpha$ | 非調和性。transmon では概ね $\alpha\simeq -E_C/\hbar$。 |
| $g$ | qubit–resonator coupling。 |
| $\Delta$ | detuning。文脈により $\omega_q-\omega_r$ または $\omega_c-\omega_t$。 |
| $J$ | qubit–qubit exchange coupling。 |
| $\chi$ | dispersive shift。qubit 状態依存の resonator 周波数シフト。 |
1. 1Q：transmon 単一量子ビット
transmon は「Josephson 非線形で少しだけ非等間隔になった LC 振動子」です。二準位だけを使いますが、物理的には多準位系です。そのため、高速パルス・DRAG・leakage・AC Stark shift を理解するには Duffing oscillator として扱う必要があります。
回路量子化から transmon Hamiltonian まで
導出 1：node flux から Hamiltonian
1
単一 Josephson 接合と shunt capacitance の自由度を node flux $\Phi$ で書く。位相は $\varphi = 2\pi\Phi/\Phi_0$。
Lagrangian:
L = 1/2 C_Σ \dot Φ^2 + E_J cos(2πΦ/Φ_0) + 2e n_g \dot Φ
= (ℏ^2 C_Σ/8e^2) \dot φ^2 + E_J cos φ + ℏ n_g \dot φ
2
共役運動量をとる。$Q=\partial L/\partial \dot\Phi=C_\Sigma\dot\Phi+2en_g$。Cooper 対数演算子を $n=Q/2e$ とすると $[\varphi,n]=i$。
Hamiltonian:
H = Q^2/(2C_Σ) - E_J cos φ
= 4E_C (n-n_g)^2 - E_J cos φ,
E_C = e^2/(2C_Σ)
$n_g$ は offset charge。transmon 領域 $E_J/E_C\gg 1$ では charge dispersion が指数関数的に小さくなり、通常のゲート導出では $n_g$ 依存を無視する。
導出 2：cos 展開 → harmonic oscillator + quartic perturbation
cos φ = 1 - φ^2/2 + φ^4/24 - φ^6/720 + ...
H ≈ 4E_C n^2 + (E_J/2)φ^2 - E_J - (E_J/24)φ^4 + O(E_J φ^6)
1
二次部分を harmonic oscillator に同一視するため、zero-point fluctuation を導入する。
φ = φ_zpf(a+a†), n = -i n_zpf(a-a†)
φ_zpf = (2E_C/E_J)^(1/4), n_zpf = (E_J/32E_C)^(1/4)
φ_zpf n_zpf = 1/2
2
二次 Hamiltonian は
H_2 = -E_J + √(8E_JE_C)(a†a + 1/2)
3
四次項を一次摂動で評価する。harmonic oscillator の期待値は $\langle m|(a+a^\dagger)^4|m\rangle=6m^2+6m+3$。
ΔE_m^(4) = -E_J/24 · φ_zpf^4 · (6m^2+6m+3)
= -E_C/12 · (6m^2+6m+3)
E_m ≈ -E_J + √(8E_JE_C)(m+1/2) - E_C(6m^2+6m+3)/12
ℏω_01 = E_1-E_0 ≈ √(8E_JE_C)-E_C
ℏω_12 = E_2-E_1 ≈ √(8E_JE_C)-2E_C
ℏα = ℏ(ω_12-ω_01) ≈ -E_C
導出 3：二準位射影と Bloch Hamiltonian
H_q = E_0 |0⟩⟨0| + E_1 |1⟩⟨1|
= (E_0+E_1)I/2 + (E_0-E_1)Z/2
global phase を捨てる： H_q/ℏ = -ω_01 Z/2
規約を変えて H_q/ℏ = +ω_q Z/2 と書くこともある。
符号規約は実験制御ソフトで混乱しやすい。最終的には「$|1\rangle$ が $|0\rangle$ より $\hbar\omega_q$ 高い」という物理だけが不変です。
導出 4：lab frame drive → rotating frame → RWA
H_lab/ℏ = ω_q Z/2 + Ω(t) cos(ω_d t + φ) X
1
$U=e^{+iω_d t Z/2}$ で回転座標へ。
H_rot = U H_lab U† + iℏ \dot U U†
2
$U X U^\dagger=X\cos ω_dt - Y\sin ω_dt$ を使う。
cos(ω_dt+φ)[X cosω_dt - Y sinω_dt]
= 1/2[ X cosφ - Y sinφ ] + terms oscillating at 2ω_d
3
RWA: $2ω_d$ で高速振動する項を平均ゼロとして捨てる。
H_RWA/ℏ = (ω_q-ω_d) Z/2 + Ω(t)/2 [X cosφ - Y sinφ]
= -Δ Z/2 + Ω_x X/2 + Ω_y Y/2
ただし Δ=ω_d-ω_q, Ω_x=Ω cosφ, Ω_y=-Ω sinφ
DRAG：leakage を抑える一次の考え方
三準位 RWA 模型：
H/ℏ = 0·|0⟩⟨0| + 0·|1⟩⟨1| + α |2⟩⟨2|
+ Ω_x(t)/2 [ |0⟩⟨1| + √2 |1⟩⟨2| + h.c.]
+ Ω_y(t)/2 [ -i|0⟩⟨1| - i√2 |1⟩⟨2| + h.c. ]
adiabatic elimination の一次条件：
Ω_y(t) ≈ -β_DRAG · \dot Ω_x(t)/α
直感：$Ω_x$ が急に変わると、瞬間固有基底が追従できず $|2\rangle$ 成分が出る。直交 quadrature に時間微分を足すことで、その非断熱成分を打ち消す。
Interactive lab：transmon 近似
decoherence：開放系の最小モデル
dρ/dt = -i[H,ρ]/ℏ + Γ_1 D[σ_-]ρ + (Γ_φ/2) D[Z]ρ
D[A]ρ = AρA† - 1/2 A†Aρ - 1/2 ρA†A
1/T_2 = 1/(2T_1) + 1/T_φ
$T_1$ は qubit 周波数付近の横ノイズ、$T_\phi$ は低周波の縦ノイズに支配されます。実験の Ramsey は $T_2^*$、echo は低周波ノイズを refocus した $T_2^{echo}$ を測ります。
近似ログ：1Q
| 近似 | 条件 | 捨てたもの | 破綻時の症状 |
|---|---|---|---|
| lumped element | 回路サイズ $\ll$ 波長 | 分布定数モード | spurious mode, package mode |
| cos の四次まで | $φ_{zpf}\ll1$、すなわち $E_J/E_C\gg1$ | $φ^6,φ^8,...$ | 周波数・非調和性のずれ |
| 二準位近似 | $|Ω|\ll |α|$、パルス帯域 $\ll |α|$ | $|2\rangle$ 以上 | leakage, DRAG 必要 |
| RWA | $Ω,|Δ|\llω_q$ | $2ω_d$ 成分 | Bloch-Siegert shift |
| Markov Lindblad | 環境相関時間が短い | 非 Markov 記憶 | 非指数緩和、TLS 揺らぎ |
2. 1Q–1R：circuit QED と分散読み出し
1 量子ビット + 1 読み出し resonator の最小模型です。核心は、qubit 状態に応じて resonator 周波数がずれることです。この状態依存周波数シフト $\chi$ を読み出し信号に変換します。
Jaynes–Cummings から dispersive Hamiltonian へ
導出 1：RWA 後の JC Hamiltonian
H/ℏ = ω_r a†a + ω_q Z/2 + g(a†σ_- + aσ_+)
ここで $a$ は resonator photon、$σ_-$ は qubit の下降演算子。counter-rotating 項 $aσ_- + a^†σ_+$ は $ω_q+ω_r$ で高速振動するため RWA で捨てます。
導出 2：Schrieffer–Wolff 変換
H = H_0 + V
H_0/ℏ = ω_r a†a + ω_q Z/2
V/ℏ = g(a†σ_- + aσ_+)
Δ = ω_q - ω_r
1
$V$ を一次で消す反 Hermitian 生成子 $S$ を選ぶ。
S = (g/Δ)(a†σ_- - aσ_+)
条件： [S,H_0] = -V
2
$e^S H e^{-S}=H+[S,H]+\frac12[S,[S,H]]+...$ を二次まで展開。
H_eff ≈ H_0 + 1/2[S,V]
[S,V]/ℏ = 2(g^2/Δ)(a†a + 1/2)Z + constant
H_disp/ℏ = (ω_r + χ Z)a†a + 1/2(ω_q + χ)Z
χ = g^2/Δ （二準位模型）
導出 3：transmon 多準位性を入れた $χ$
ω_{j,j+1}=ω_q+jα, g_j≈g√(j+1), Δ_j=ω_{j,j+1}-ω_r=Δ+jα
resonator shift for qubit level j:
δω_r^{(j)} = g_j^2/Δ_j - g_{j-1}^2/Δ_{j-1}, g_{-1}=0
δω_r^{(0)} = g^2/Δ
δω_r^{(1)} = 2g^2/(Δ+α) - g^2/Δ
2χ = δω_r^{(1)} - δω_r^{(0)}
= 2g^2/(Δ+α) - 2g^2/Δ
= -2g^2 α/[Δ(Δ+α)]
したがって、ここでは $χ=-g^2 α/[Δ(Δ+α)]$。文献により $H=(ω_r-χZ)a†a$ と置くため符号が反転することがあります。
読み出し pointer state の式
drive frame で、qubit 状態 s∈{g,e} に対する resonator coherent amplitude:
\dot α_s = -[κ/2 + i(ω_r^s-ω_d)]α_s - iε(t)
steady state:
α_s^{ss} = -iε / [κ/2 + i(ω_r^s-ω_d)]
SNR ≈ √(ηκT_int) |α_e-α_g|
読み出しは「$|g\rangle$ と $|e\rangle$ が作る IQ 平面上の coherent state の分離」を測る問題です。
Purcell 効果と critical photon number
Γ_Purcell ≈ κ (g/Δ)^2
n_crit ≈ Δ^2/(4g^2)
強い読み出しは SNR を上げますが、分散近似の破綻、qubit excitation、MIST、resonator 非線形性を引き起こします。
Interactive lab：分散シフト
測定の誤差要因
| 誤差 | 物理 | 対策 |
|---|---|---|
| overlap error | IQ 分布が重なる | $χ,κ,T_{int},η$ の最適化、matched filter |
| $T_1$ 中の decay | 測定中に $|1\rangle\to|0\rangle$ | 短時間化、Purcell filter、SNR 向上 |
| residual photon | 測定後に cavity photon が残る | depletion pulse、待ち時間、filter 設計 |
| MIST/leakage | 強駆動で JTC ladder の共鳴に乗る | photon 数制限、detuning 設計、leakage reset |
近似ログ：1Q–1R
| 近似 | 条件 | 捨てたもの | 破綻 |
|---|---|---|---|
| JC RWA | $g\llω_q,ω_r$ | counter-rotating | Bloch-Siegert, non-RWA transition |
| dispersive SW | $|g/Δ|\ll1$ | 高次 hybridization | 強混合、QND 性低下 |
| linear resonator | 低 photon、Kerr 小 | resonator 非線形 | bifurcation, response distortion |
| Gaussian noise readout | 増幅器雑音が支配的 | 非 Gaussian outlier | threshold fidelity 低下 |
3. 2Q：2 量子ビットの相互作用
2Q の基本は、2 つの Duffing oscillator を exchange coupling $J$ で結合した模型です。計算基底だけを見ると $XX+YY$ と $ZZ$ が現れますが、$ZZ$ と leakage は多準位性に由来します。
2 transmon Hamiltonian と static ZZ
導出 1：Duffing 2体系
H/ℏ = Σ_{i=1,2} [ω_i b_i†b_i + α_i/2 · b_i†b_i(b_i†b_i-1)]
+ J(b_1†b_2 + b_1 b_2†)
$J$ は直接容量結合、bus resonator、tunable coupler などから生じる有効交換結合です。
導出 2：二準位射影で $XX+YY$
b_i → σ_i^- , b_i† → σ_i^+
J(σ_1^+σ_2^- + σ_1^-σ_2^+)
= J/2 (X_1X_2 + Y_1Y_2)
この項は $|10\rangle\leftrightarrow|01\rangle$ を交換します。共鳴に近いと iSWAP 型、遠 detuned では仮想遷移として周波数シフトを作ります。
導出 3：static ZZ の二次摂動
定義は dressed energy で
ζ = E_{11}-E_{10}-E_{01}+E_{00}
主な寄与は $|11\rangle$ が $|20\rangle,|02\rangle$ と仮想的に混ざること。
E_{10}=ω_1, E_{01}=ω_2, E_{11}=ω_1+ω_2
E_{20}=2ω_1+α_1, E_{02}=2ω_2+α_2
Δ=ω_1-ω_2
⟨20|V|11⟩ = √2 J, E_{11}-E_{20}= -Δ-α_1
⟨02|V|11⟩ = √2 J, E_{11}-E_{02}= Δ-α_2
δE_{11}^{(2)} = 2J^2/( -Δ-α_1 ) + 2J^2/( Δ-α_2 )
δE_{10}^{(2)} + δE_{01}^{(2)} = J^2/Δ - J^2/Δ = 0
ζ ≈ 2J^2[1/(Δ-α_2) - 1/(Δ+α_1)]
= -2J^2(α_1+α_2)/[(Δ+α_1)(α_2-Δ)]
分母が小さい点は avoided crossing / frequency collision。摂動式は使えません。
CZ と iSWAP の見方
iSWAP: exp[-i θ/2 (XX+YY)] → |10⟩,|01⟩ の交換
CZ: diag(1,1,1,-1) → φ_ZZ = φ_11-φ_10-φ_01+φ_00 = π
CZ は $|11\rangle$ を $|02\rangle$ 近傍に近づけて条件付き位相を蓄積する、または tunable coupler で $ZZ$ 相互作用を動的に作る方式として実装されます。
Interactive lab：ZZ と衝突点
2Q 設計で避けるべき周波数条件
| 条件 | 意味 | 症状 |
|---|---|---|
| $ω_1≈ω_2$ | $|10\rangle$ と $|01\rangle$ が混合 | idle swap, addressability 低下 |
| $ω_1+ω_2≈2ω_2+α_2$ | $|11\rangle≈|02\rangle$ | ZZ 増大、leakage。ただし CZ に使える |
| $ω_1+ω_2≈2ω_1+α_1$ | $|11\rangle≈|20\rangle$ | 同上 |
| drive harmonics | CR/flux pulse の sideband が他遷移に当たる | spectator error, leakage |
近似ログ：2Q
| 近似 | 条件 | 破綻時 |
|---|---|---|
| Duffing truncation | 低励起、transmon 領域 | 高準位の実 spectrum がずれる |
| exchange RWA | $J\llω_i$ | counter-rotating coupling |
| 二次摂動 ZZ | $J\ll |Δ|,|Δ+α_1|,|Δ-α_2|$ | avoided crossing では full diagonalization |
| isolated pair | spectator qubits が遠 detuned | many-body ZZ, crosstalk |
4. 1Q のキャリブレーション
1Q calibration は「Hamiltonian の未知パラメータを測る」ことと「パルス誤差を増幅して補正する」ことの組み合わせです。Rabi は回転角、Ramsey は detuning、DRAG は leakage と quadrature error、RB は平均誤差を測ります。
測定シーケンスと数式
Rabi：amplitude → 回転角
共鳴駆動 Δ=0:
U(A,T)=exp[-i X/2 ∫_0^T Ω(A,t)dt]
θ(A)=∫Ω(A,t)dt
|0⟩ から開始： P_1(A)=sin^2[θ(A)/2]
線形領域なら θ(A)=cA, π pulse は cA_π=π
Ramsey：detuning を fringe 周波数に変換
X_{π/2} - wait τ - X_{π/2}
wait 中の Hamiltonian: H/ℏ = -Δ Z/2
P_1(τ)=1/2[1 + C exp(-τ/T_2^*) cos(Δτ+φ)]
fringe 周波数が意図した detuning とずれていれば、qubit 周波数または drive 周波数を更新します。
Echo：低周波 dephasing の refocus
X_{π/2} - τ/2 - X_π - τ/2 - X_{π/2}
準静的 detuning $ξZ/2$ は前半と後半で符号が反転するため相殺されます。高周波ノイズは残ります。
DRAG calibration：何を最小化するか
| 観測量 | シーケンス例 | 調整 |
|---|---|---|
| leakage | repeated $X_π$ 後の $|2\rangle$ 読み出し | DRAG 係数 |
| axis error | $X_{π/2}Y_{π/2}$ 型の誤差増幅 | IQ phase / DRAG |
| AC Stark $Z$ | interleaved Ramsey | drive detuning / virtual Z |
RB：平均 gate error への変換
Clifford 長 m の survival:
P(m)=A p^m + B
single-qubit average error:
r = (1-p)/2
$A,B$ は SPAM を吸収する係数。RB は coherent error を depolarizing parameter に圧縮するため、誤差の種類の同定には別途 tomography や error amplification が必要です。
Interactive lab：Rabi/Ramsey
1Q calibration の実務フロー
resonator spectroscopy
→
qubit spectroscopy
→
Rabi
→
Ramsey
→
DRAG
→
readout IQ
→
RB / IRB
実験ではこの順序を一度だけ行うのではなく、ドリフト監視と自動再較正でループさせます。1Q gate の安定性は、制御器の振幅・位相ドリフト、温度、LO/PLL、IQ mixer の安定性にも依存します。
5. CR：Cross Resonance gate
CR は固定周波数 transmon の代表的 2Q gate です。control qubit に target 周波数の drive を当てると、qubit 間結合と dressed basis により、target 回転が control 状態に依存します。欲しい有効項は $ZX$ です。
CR の有効 Hamiltonian
出発点：2 qubit + control drive
H_0/ℏ = ω_c Z_c/2 + ω_t Z_t/2 + J/2 (X_cX_t+Y_cY_t)
H_d/ℏ = Ω(t) cos(ω_t t+φ) X_c
drive は control line に入れるが、周波数は target の $ω_t$。control から見ると off-resonant drive、dressed basis では target への条件付き drive 成分を持ちます。
Pauli 展開
H_CR/ℏ = 1/2[ ω_IX IX + ω_IY IY + ω_IZ IZ
+ ω_ZX ZX + ω_ZY ZY + ω_ZZ ZZ + ω_ZI ZI + ...]
理想 entangler:
U_ZX(θ)=exp[-i θ ZX/2]
θ=∫ω_ZX(t)dt
θ=π/2 なら local gates と合わせて CNOT と等価
弱駆動・多準位 transmon でのスケーリング
ω_ZX ≈ - J Ω / Δ · α_c/(Δ+α_c)
ここで $Δ=ω_c-ω_t$、$α_c$ は control transmon の非調和性。厳密な係数と符号は Hamiltonian 規約、drive 位相、dressed basis 定義で変わります。
この式が伝える本質：$J=0$ なら出ない、$Ω$ にほぼ比例、$Δ$ が大きいと遅い、$Δ≈0$ や $Δ≈-α_c$ などの衝突点は危険。
echo の代数：なぜ不要項が消えるか
Echoed CR の概念：
CR(+Ω) - X_π^c - CR(-Ω) - X_π^c
| 項 | $X_π^c$ での符号 | $Ω→-Ω$ での符号 | echo 後 |
|---|---|---|---|
| $IX,IY$ | 不変 | 反転 | キャンセルしやすい |
| $ZX,ZY$ | 反転 | 反転 | 残る |
| $ZI,ZZ$ | 反転 | 主に不変 | キャンセル傾向 |
実機では finite rise/fall、drive 非線形、Stark shift、crosstalk のため完全には消えず、cancellation tone と virtual Z 補正を併用します。
Interactive lab：CR rate
CR gate の主な誤差
| 項 | 意味 | 補償 |
|---|---|---|
| $IX/IY$ | target の無条件回転。classical crosstalk でも出る。 | echo, cancellation tone |
| $ZI$ | control の Stark shift | virtual Z |
| $IZ$ | target の Stark shift | virtual Z |
| $ZZ$ | 条件付き位相、static + drive-induced | echo, phase compensation |
| leakage | $|2\rangle$ 以上への遷移 | amplitude 制限、pulse shaping、周波数設計 |
近似ログ：CR
| 近似 | 条件 | 破綻 |
|---|---|---|
| 弱駆動摂動 | $Ω/|Δ|\ll1$、高準位衝突なし | Pauli 係数が非線形、leakage |
| pair-only model | spectator が遠い | spectator CR, crowding |
| linear crosstalk | target への漏れ drive が小さい | IY/IX の非線形変化 |
| echo cancellation | 2 pulse が対称 | rise/fall, drift で残差 |
6. 2Q のキャリブレーション
2Q calibration は、有効 Hamiltonian の Pauli 係数を測り、欲しい非局所角だけを残し、局所位相・crosstalk・leakage を補償する作業です。
CR calibration の導出と測定
control 状態別 target Rabi から $IX,ZX$ を分離
H_eff/ℏ ≈ 1/2[(ω_IX I + ω_ZX Z_c)X_t + (ω_IY I + ω_ZY Z_c)Y_t + ...]
control を $|0\rangle$、$|1\rangle$ に準備すると、target の有効回転ベクトルは
control |0⟩: Ω⃗_0 = (ω_IX + ω_ZX, ω_IY + ω_ZY, ...)
control |1⟩: Ω⃗_1 = (ω_IX - ω_ZX, ω_IY - ω_ZY, ...)
ω_IX = (Ω_x^0 + Ω_x^1)/2, ω_ZX = (Ω_x^0 - Ω_x^1)/2
ω_IY = (Ω_y^0 + Ω_y^1)/2, ω_ZY = (Ω_y^0 - Ω_y^1)/2
entangling angle の合わせ込み
θ_ZX = ∫_0^T ω_ZX(t) dt
CNOT 等価ゲートでは、実装の分解に応じて $θ_{ZX}=π/2$ を目標にします。echo pulse, rotary tone, cancellation tone が入ると実効面積は単純な矩形面積ではなく、実測 fringe で合わせます。
local phase compensation
U_meas ≈ exp[-i( θ_ZX ZX + θ_ZI ZI + θ_IZ IZ + θ_ZZ ZZ + ...)/2]
$ZI,IZ$ は virtual Z で補償できます。$ZZ$ は局所ゲートでは消えないため、echo・pulse shape・追加補償が必要です。
CZ/tunable coupler calibration
φ_ZZ = φ_11 - φ_10 - φ_01 + φ_00
CZ 条件： φ_ZZ = π mod 2π, leakage P_leak small
- flux amplitude で avoided crossing を探索。
- duration を掃引して $φ_{ZZ}$ を測る。
- $|02\rangle,|20\rangle$ leakage を最小化する波形にする。
- 単一 qubit 位相を virtual Z で戻す。
- 2Q RB / XEB / cycle benchmarking で評価。
Interactive lab：Hamiltonian tomography
2Q benchmarking の違い
| 手法 | 測るもの | 注意 |
|---|---|---|
| 2Q RB | 平均 Clifford error | 漏れや coherent error の詳細は潰れる |
| interleaved RB | 特定 gate の相対誤差 | 参照 RB 依存 |
| XEB | ランダム回路の fidelity | モデル依存、SPAM 処理 |
| cycle benchmarking | 特定 cycle の Pauli error | QEC サイクルに近い評価 |
| leakage RB | 計算部分空間外への遷移 | leakage readout calibration が必要 |
2Q calibration checklist
- pair spectroscopy: $ω_c,ω_t,α_c,α_t,J,ζ$。
- static ZZ: target Ramsey with control $|0\rangle/|1\rangle$。
- CR amplitude/duration sweep: $ZX$ rate。
- Hamiltonian tomography: $IX,IY,ZX,ZY,ZZ,ZI,IZ$。
- echo/cancellation tone/virtual Z。
- leakage measurement and reset compatibility。
- 2Q RB + simultaneous RB。
7. その他：実験系・制御系・最適制御
超伝導量子ビットは Hamiltonian だけでは動きません。希釈冷凍機、マイクロ波配線、フィルタ、増幅器、AWG/ADC、FPGA、制御ソフトウェア、データ解析、最適制御が一体になって動作します。
希釈冷凍機と熱励起
thermal excited population:
P_1/P_0 = exp[-ℏω_q/(k_B T)]
5 GHz → ℏω/k_B ≈ 240 mK
T=10 mK なら exp(-24) ≈ 3.8×10^-11
実際の qubit 有効温度は、赤外線、同軸線雑音、残留 photon、非平衡 quasiparticle によって冷凍機温度より高くなることがあります。
制御スタック
compiler
→
pulse schedule
→
AWG/DAC
→
IQ/LO
→
chip
→
resonator
→
JPA/TWPA/HEMT
→
ADC/FPGA
実験 API では、qubit、readout resonator、coupler を物理チャネルに対応づけ、pulse sequence をデバイスの言葉で記述します。
量子最適制御
目的関数と制約
minimize L = 1 - F_gate
+ λ_leak P_leak
+ λ_amp ∫|Ω(t)|^2 dt
+ λ_bw ∫|dΩ/dt|^2 dt
+ λ_robust E_ξ[1-F(ξ)]
subject to:
|Ω(t)|≤Ω_max, bandwidth≤B, waveform samples compatible with AWG,
P_leak≤ε_leak, calibration time finite
GRAPE の最小導出
時間を N piece に分割：
U(T)=U_N...U_2U_1, U_k=exp[-i(H_0+Σ_j u_{j,k}H_j)Δt]
fidelity 例： F = |Tr(U_target† U(T))|^2/d^2
gradient は前向き伝搬 $U_k...U_1$ と後ろ向き伝搬 $U_{target}^†U_N...U_{k+1}$ を使って効率よく計算します。実機では model mismatch のため閉ループ最適化や calibration 実験と併用します。
Interactive lab：システム全体
実務でよく起きる「理論と実機の差」
| 差分 | 原因 | 観測方法 | 対策 |
|---|---|---|---|
| 波形歪み | ケーブル、フィルタ、IQ mixer、skin effect | scope, cryoscope, Ramsey distortion | predistortion |
| crosstalk | line coupling、package mode、制御器漏れ | simultaneous RB, Hamiltonian tomography | cancellation matrix, scheduling |
| frequency drift | TLS, flux noise, 温度、quasiparticle | tracking Ramsey | adaptive recalibration |
| readout outlier | MIST, leakage, amplifier nonlinearity | IQ clustering, leakage readout | readout power制限、reset |
近似・用語集：どこで何を捨てたか
超伝導量子ビットの理解で最も重要なのは、式そのものよりも「その式がどの条件で正しいか」です。このページは全セクションの近似を一覧化します。
近似一覧
| 近似 | 数学操作 | 小パラメータ | 捨てた項 | 破綻時に使うもの |
|---|---|---|---|---|
| cos 展開 | $\cosφ$ を Taylor 展開 | $φ_{zpf}$ | $φ^6$ 以上 | charge basis / numerical diagonalization |
| Duffing | quartic perturbation | $\sqrt{E_C/E_J}$ | 高次非線形、charge dispersion | full transmon spectrum |
| 二準位 | projector $P=|0\rangle\langle0|+|1\rangle\langle1|$ | $Ω/α$ | $|2\rangle$ 以上 | multi-level simulation |
| RWA | 高速振動項の平均化 | $Ω/ω$, $g/ω$ | $e^{±2iωt}$ | Floquet / lab-frame integration |
| SW/dispersive | unitary block diagonalization | $g/Δ$ | 高次 hybridization | exact diagonalization / dressed basis |
| ZZ 摂動 | 二次エネルギー補正 | $J/δE$ | 高次混合 | full multi-level diagonalization |
| Lindblad | Markov bath、Born 近似 | bath correlation time / system time | memory | non-Markov model, TLS model |
| RB depolarizing fit | Clifford twirl | gate-independent noise 近似 | 誤差構造 | GST, cycle benchmarking, tomography |
よく使う単位変換
1 GHz ↔ h·1GHz/k_B ≈ 48 mK
5 GHz ↔ 240 mK
angular frequency と ordinary frequency:
ω = 2π f
論文の MHz が ω/2π なのか ω なのか必ず確認。
記号の符号に関する注意
$Z=|0\rangle\langle0|-|1\rangle\langle1|$ と置くか、逆に置くかで Hamiltonian の符号が変わります。$Δ=ω_d-ω_q$ と置くか $ω_q-ω_d$ と置くかでも Ramsey 式の符号が変わります。実験上の不変量は、測定される fringe 周波数、回転方向、エネルギー差です。
資料中では、物理量の依存性を見やすくするため、符号規約の違いがあり得る式には注記を入れています。
参考文献・出典メモ
この HTML は、アップロード済み PDF 群と標準的な circuit-QED/transmon 理論に基づき、教育用に再構成したものです。
| 主題 | 参照 |
|---|---|
| 超伝導 qubit 総説、transmon の非等間隔性 | Kjaergaard et al., “Superconducting Qubits: Current State of Play”. |
| circuit QED / dispersive readout | Blais et al., “Cavity quantum electrodynamics for superconducting electrical circuits”. |
| DRAG | Motzoi et al., “Simple Pulses for Elimination of Leakage in Weakly Nonlinear Qubits”. |
| CR effective Hamiltonian | Magesan & Gambetta, “Effective Hamiltonian models of the cross-resonance gate”. |
| ZZ-free / dynamically corrected 1Q gate | Watanabe et al., “ZZ-Interaction-Free Single-Qubit-Gate Optimization in Superconducting Qubits”. |
| Measurement-induced state transitions | Khezri et al., “Measurement-induced state transitions in a superconducting qubit”. |
| dynamic dispersive shift readout | Swiadek et al., “Enhancing Dispersive Readout...”. |
| 量子最適制御 | Koch et al., “Quantum optimal control in quantum technologies”. |
| 制御ソフトウェア / QCS | Keysight Quantum Control System (QCS) data sheet. |
| 希釈冷凍機 | 福山寛「希釈冷凍機の原理としくみ」. |