超導量子電腦架構

上週 Transmon 展示了如何選擇 \(E_J/E_C\)，在 charge sensitivity 與 anharmonicity 之間取得平衡。 但降低 parameter sensitivity 並不等於所有 error transition 都消失。這週進一步追問：我們能否直接設計量子態、 coupling operator 與 Hamiltonian term，使環境難以引發 relaxation；而當環境難以驅動 transition 時，我們自己是否也會更難控制它？

本週驅動問題：Fluxonium 透過大電感與 Josephson junction 形成不同於 Transmon 的能譜；哪些 transition 被保護、哪些 control／readout 變得更難？如何避免把單一 coherence 紀錄誤認成整體架構優勢？

查看 T2 十問的整體推進

不從 transition-rate 公式開始背，而從「系統、環境與兩者的耦合」開始。最一般的弱耦合寫法是

\[ \hat H_{\rm tot} =\hat H_{\rm S}(\lambda_0)+\hat H_{\rm B}+\hat H_{\rm SB}, \qquad \hat H_{\rm S}(\lambda_0)\ket n=E_n\ket n. \]

\(\hat H_{\rm S}\) 決定能階與 eigenstates，\(\hat H_{\rm B}\) 描述環境；真正決定環境能否改變 qubit state 的是 \(\hat H_{\rm SB}\)。若某個 circuit parameter 受到小擾動 \(\lambda(t)=\lambda_0+\delta\lambda(t)\)，直接對 Hamiltonian 展開：

\[ \begin{aligned} \hat H_{\rm S}(\lambda_0+\delta\lambda) &\simeq \hat H_0 +\delta\lambda(t) \left.\frac{\partial\hat H_{\rm S}}{\partial\lambda}\right|_{\lambda_0} +O(\delta\lambda^2)\\ &\equiv \hat H_0+F(t)\hat O. \end{aligned} \]

因此 \(F(t)\) 不是抽象的「noise」而已，它是實際 fluctuation；\(\hat O=\partial_\lambda\hat H_{\rm S}\) 則由 circuit Hamiltonian 決定。 若環境也要量子化，寫成 \(\hat H_{\rm SB}=\sum_\alpha\hat O_\alpha\otimes\hat B_\alpha\)；在 semiclassical noise model 中， \(\hat B_\alpha\) 由 stochastic function \(F_\alpha(t)\) 取代。兩種寫法都保留同一個 system-side coupling operator。

第一步：把 coupling operator 展開到 energy eigenbasis

\[ \hat O=\sum_{m,n}O_{mn}\ket m\!\bra n, \qquad O_{mn}=\langle m|\hat O|n\rangle. \]

這一步立刻把同一個 perturbation 分成兩種物理作用：

\[ \hat V(t)=F(t) \left[ \underbrace{\sum_n O_{nn}\ket n\!\bra n}_{\text{diagonal：改變相位與 transition frequency}} + \underbrace{\sum_{m\ne n}O_{mn}\ket m\!\bra n}_{\text{off-diagonal：改變 population}} \right]. \]

• Diagonal part：每個能階得到不同的瞬時 energy shift，長時間累積成隨機 phase，造成 dephasing。
• Off-diagonal part：若 \(O_{fi}\neq0\)，perturbation 才有可能把 \(\ket i\) 連到 \(\ket f\)，造成 excitation 或 relaxation。

第二步：在 interaction picture 中推導 transition amplitude

定義 Bohr frequency \(\omega_{mn}=(E_m-E_n)/\hbar\)。在 \(\hat H_0\) 的 interaction picture 中，

\[ \hat V_I(t)=F(t)\sum_{m,n}O_{mn}e^{i\omega_{mn}t}\ket m\!\bra n. \]

Interaction-picture Schrödinger equation 與 Dyson expansion 的一階項是

\[ i\hbar\frac{d}{dt}\ket{\psi_I(t)}=\hat V_I(t)\ket{\psi_I(t)}, \qquad \hat U_I(t,0)\simeq 1-\frac{i}{\hbar}\int_0^t dt'\,\hat V_I(t'). \]

若系統起始於 \(\ket i\)，投影到 \(\ket f\) 就得到一階 \(i\to f\) amplitude：

\[ c_f^{(1)}(t) =-\frac{i}{\hbar}O_{fi} \int_0^t dt'\,F(t')e^{i\omega_{fi}t'}. \]

所以 transition probability 是

\[ P_{i\to f}(t) =\frac{|O_{fi}|^2}{\hbar^2} \left| \int_0^t dt'\,F(t')e^{i\omega_{fi}t'} \right|^2. \]

這一行已經給出兩個必要條件：第一，\(O_{fi}\) 不能被 symmetry 強迫為零；第二，\(F(t)\) 必須在 \(|\omega_{fi}|\) 附近具有 Fourier component。 有限觀察時間下的共振寬度約為 \(1/t\)；長時間極限才逐漸逼近嚴格的 energy conservation。

第三步：同一條公式分別導出 coherent drive 與 stochastic-noise transition

對人為 drive，令 \(F(t)=F_d\cos\omega_dt\)。在 \(\omega_d\simeq|\omega_{fi}|\) 且 rotating-wave approximation 成立時，

\[ \Omega_{fi}=\frac{F_d|O_{fi}|}{\hbar}, \qquad P_{i\to f}(t)= \frac{\Omega_{fi}^2}{\Omega_{fi}^2+\Delta^2} \sin^2\!\left(\frac{t}{2}\sqrt{\Omega_{fi}^2+\Delta^2}\right), \quad \Delta=\omega_d-|\omega_{fi}|. \]

這說明 control transition 的強弱由 drive-operator matrix element 決定。對 stationary random noise，則定義 correlation function 與 two-sided spectral density：

\[ C_F(\tau)=\langle F(\tau)F(0)\rangle, \qquad S_F(\omega)=\int_{-\infty}^{\infty}d\tau\, e^{i\omega\tau}C_F(\tau). \]

把 amplitude 取 ensemble average，得到

\[ \begin{aligned} \big\langle P_{i\to f}(t)\big\rangle &=\frac{|O_{fi}|^2}{\hbar^2} \int_0^t dt_1\int_0^t dt_2\, e^{i\omega_{fi}(t_1-t_2)}C_F(t_2-t_1)\\ &\xrightarrow[t\gg\tau_c]{} t\,\frac{|O_{fi}|^2}{\hbar^2}S_F(\omega_{if}), \qquad \omega_{if}=\frac{E_i-E_f}{\hbar}. \end{aligned} \]

因此 Fermi golden rule 的 transition rate 為

\[ \boxed{ \Gamma_{i\to f} =\frac{1}{\hbar^2} \left|\langle f|\hat O|i\rangle\right|^2 S_F(\omega_{if})}. \]

這個結果使用 weak coupling、stationary bath、\(t\) 長於 bath correlation time，以及可忽略高階 back-action 的近似。 若教材把 \(2\pi\)、\(\hbar\) 或 density of states 吸收到 spectral density 的定義中，前因子會寫成 \(2\pi/\hbar\) 等形式；必須先對照 \(S_F\) 的定義，不能只比公式外觀。

第四步：由 upward／downward rates 得到 \(T_1\)

令 \(\omega_{01}=(E_1-E_0)/\hbar>0\)，同一個 noise channel 給出

\[ \Gamma_\downarrow =\frac{|O_{01}|^2}{\hbar^2}S_F(+\omega_{01}), \qquad \Gamma_\uparrow =\frac{|O_{01}|^2}{\hbar^2}S_F(-\omega_{01}), \qquad \frac{1}{T_1}=\Gamma_\downarrow+\Gamma_\uparrow. \]

對 thermal bath，detailed balance 給出 \(S_F(-\omega)=e^{-\hbar\omega/k_BT}S_F(+\omega)\)。低溫時通常 \(\Gamma_\uparrow\ll\Gamma_\downarrow\)，但低頻 qubit 或有效環境溫度偏高時， upward transition 與 thermal population 不能忽略。若有多個不相關 channel，總 rate 是各 channel rate 相加；相關 channel 則需要 cross-spectral density \(S_{\alpha\beta}(\omega)\)。 這裡的 upward／downward rate 使用 unsymmetrized quantum-noise spectrum；symmetrized spectrum 會丟掉能量流向，不能單獨區分兩者。

第五步：回到 diagonal coupling，推導 dephasing 與 sweet spot

Diagonal part 不直接交換 population，而是改變兩能階之差：

\[ \delta\omega_{01}(t) =\frac{F(t)}{\hbar}\bigl(O_{11}-O_{00}\bigr), \qquad \delta\phi(t)=\int_0^t dt'\,\delta\omega_{01}(t'). \]

若 \(F(t)=\delta\lambda(t)\) 且 \(\hat O=\partial_\lambda\hat H_0\)，Hellmann–Feynman theorem 直接連接 diagonal matrix element 與 frequency sensitivity：

\[ O_{nn}=\left\langle n\left|\frac{\partial\hat H_0}{\partial\lambda}\right|n\right\rangle =\frac{\partial E_n}{\partial\lambda}, \qquad \frac{O_{11}-O_{00}}{\hbar} =\frac{\partial\omega_{01}}{\partial\lambda}. \]

因此

\[ \delta\omega_{01}(t) \simeq \frac{\partial\omega_{01}}{\partial\lambda}\delta\lambda(t) +\frac12\frac{\partial^2\omega_{01}}{\partial\lambda^2}\delta\lambda^2(t)+\cdots. \]

Sweet spot \(\partial_\lambda\omega_{01}=0\) 只消除一階 frequency fluctuation；二階 sensitivity、其他 noise parameter 與 finite-frequency relaxation 仍然存在。 Ramsey、echo 或 dynamical-decoupling sequence 量到的 dephasing，則是低頻 \(S_\lambda(\omega)\) 經各自 filter function 加權後的結果。

第六步：把抽象的 \(F(t)\hat O\) 映射回 superconducting circuit

以含 charge、Josephson 與 inductive energy 的 Hamiltonian 為例：

\[ \hat H_0 =4E_C(\hat n-n_g)^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2. \]

小幅 parameter noise 不是額外猜一個 operator，而是直接微分：

\[ \delta\hat H(t) \simeq \delta n_g(t)\underbrace{\big[-8E_C(\hat n-n_g)\big]}_{\hat O_{n_g}} +\delta\varphi_{\rm ext}(t) \underbrace{\big[-E_L(\hat\varphi-\varphi_{\rm ext})\big]}_{\hat O_{\Phi}} +\delta E_J(t)\underbrace{\big[-\cos\hat\varphi\big]}_{\hat O_{E_J}}. \]

若使用實際 flux \(\Phi_{\rm ext}\)，還要乘上 \(\partial\varphi_{\rm ext}/\partial\Phi_{\rm ext}=2\pi/\Phi_0\)。同樣地，電容式 microwave drive 通常透過 \(\hat n\) 耦合。 所以同一組 eigenstates 必須分別計算不同 operator 的 diagonal sensitivity 與 off-diagonal matrix element。

Charge-like noise：\(\hat n\) 與 offset charge

Charge noise 可來自 oxide defect、two-level system（TLS）或 gate-voltage fluctuation，並耦合到 Cooper-pair number：

\[ \hat V_{\rm charge}(t) =-8E_C\,\delta n_g(t)(\hat n-n_g) \simeq-8E_C\,\delta n_g(t)\hat n, \qquad \delta n_g=\frac{C_g\delta V_g}{2e}. \]

低頻 charge noise 常呈現 \(1/f\) 特性。Transmon 增大 \(E_J/E_C\) 後，受到指數壓低的是能量對 offset charge \(n_g\) 的 charge dispersion；這是 frequency sensitivity 的鈍化，不能等同於 \(\langle0|\hat n|1\rangle\) 指數下降。

Flux-like noise：低頻 dephasing 與 finite-frequency transition 必須分開

Flux noise 可來自 surface-spin bath。低頻部分主要透過 frequency sensitivity 造成 dephasing：

\[ \text{dephasing weight} \propto \left(\frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}}\right)^2 S_\Phi(\omega\approx0). \]

同一個 flux channel 在 \(\omega_{01}\) 附近的 noise，則透過 off-diagonal matrix element 造成 transition：

\[ \Gamma_{1\to0}^{(\Phi)} =\frac{1}{\hbar^2} \left| \left\langle0\left| \frac{\partial\hat H_0}{\partial\Phi_{\rm ext}} \right|1\right\rangle \right|^2 S_\Phi(+\omega_{01}). \]

所以 sweet spot 與 relaxation protection 不是同一個敘述；前者控制 diagonal sensitivity，後者還要檢查 finite-frequency spectrum 與 off-diagonal matrix element。

完整判斷鏈

在 weak-coupling、Markovian exponential decay 的近似下，整條推導可壓成：

\[ \boxed{ \text{circuit 與 bias} \Rightarrow \hat H_0(\lambda) \Rightarrow \{E_n,\ket n\} \Rightarrow \hat O_\lambda=\partial_\lambda\hat H_0 \Rightarrow \begin{cases} O_{nn}\ \text{與}\ S_\lambda(0) &\Rightarrow T_\phi,\\ O_{mn}\ \text{與}\ S_\lambda(\omega_{nm}) &\Rightarrow T_1\ \text{與 excitation}, \end{cases} \Rightarrow T_2^{-1}=(2T_1)^{-1}+T_\phi^{-1}.} \]

外界擾動造成某一條 transition，必須同時滿足：環境經由明確的 \(\hat O\) 耦合、\(\langle f|\hat O|i\rangle\neq0\)、 且 \(S_F\) 在所需 Bohr frequency 有 spectral weight。Allowed 只表示不被 symmetry 禁止，不代表 transition 一定很強； forbidden 是 exact zero，suppressed 則是非零但很小。比較 coherence 時也不能只看單一 \(T_1\)：matrix-element suppression 與較低 environmental spectral density 對 device scaling 的意義不同。

問題一已得到 transition 的兩個乘因子：environmental spectrum 與 system matrix element。問題二先固定環境無法完全消失，改問能否設計 \(\ket0,\ket1\) 的 wavefunction、symmetry 與 logical encoding，使 dominant noise operator 幾乎看不見 logical transition， 同時讓另一個可控制的 operator 仍能作用。

環境造成的 relaxation 與人為施加的 control 具有相同的數學骨架：

\[ \hat V_{\rm noise}(t)=\hat O_{\rm noise}F_{\rm noise}(t), \qquad \Gamma_\downarrow= \frac{1}{\hbar^2} \left|\langle0|\hat O_{\rm noise}|1\rangle\right|^2 S_{\rm noise}(+\omega_{01}), \]

\[ \hat V_{\rm drive}(t)=\hat O_{\rm drive}A\cos\omega_dt, \qquad \Omega_R=\frac{A}{\hbar} \left|\langle0|\hat O_{\rm drive}|1\rangle\right|. \]

兩式的差別不在「noise 是壞的、drive 是好的」，而在它們是否經由同一個 physical port 與同一個 operator 耦合。 若 \(\hat O_{\rm noise}=\hat O_{\rm drive}\)，壓低 qubit matrix element 會同時延長 \(T_1\) 與降低 Rabi rate； 若兩者具有不同的 symmetry、spatial profile 或 frequency dependence，才可能對 environment 弱、對 control 強。

第一步：把所有 operator 投影到 logical subspace

定義 logical projector

\[ \hat P_{\rm L}=\ket0\!\bra0+\ket1\!\bra1. \]

任一 Hermitian operator 在 logical subspace 中都可以分解成

\[ \hat P_{\rm L}\hat O\hat P_{\rm L} =\bar O\,\hat I +\frac{O_{11}-O_{00}}{2}\hat\sigma_z +\operatorname{Re}O_{01}\,\hat\sigma_x -\operatorname{Im}O_{01}\,\hat\sigma_y, \qquad \bar O=\frac{O_{00}+O_{11}}{2}. \]

因此 state protection 的物理要求可以直接讀成：

• 環境難以分辨 logical states：\(O_{11}^{\rm noise}-O_{00}^{\rm noise}\approx0\)，降低一階 dephasing。
• 環境難以交換 logical population：\(O_{01}^{\rm noise}\approx0\)，降低 relaxation／excitation。
• 人為控制不是 identity：\(O_{01}^{\rm drive}\neq0\)，保留 transverse control；需要 logical \(Z\) control 時，drive 的 diagonal difference 也必須可調。

理想情況可濃縮為

\[ \boxed{ \hat P_{\rm L}\hat O_{\rm noise}\hat P_{\rm L} \simeq c\hat I, \qquad \left|\langle0|\hat O_{\rm drive}|1\rangle\right| \neq0.} \]

這個條件必須逐一套用到 dominant noise channels。對 charge noise 成立，不代表對 flux、critical-current、dielectric-loss 或 quasiparticle channel 也成立。

第二步：exact symmetry 如何產生 selection rule？

假設存在 parity-type unitary symmetry \(\hat{\mathcal P}\)，滿足 \(\hat{\mathcal P}^2=\hat I\) 與 \([\hat H_0,\hat{\mathcal P}]=0\)。能量 eigenstates 可同時選成 symmetry eigenstates：

\[ \hat{\mathcal P}\ket n=p_n\ket n, \qquad p_n=\pm1. \]

若 coupling operator 在此 symmetry 下具有 parity \(\eta_O=\pm1\)，

\[ \hat{\mathcal P}\hat O\hat{\mathcal P}^{\dagger}=\eta_O\hat O, \]

則 matrix element 必須滿足

\[ \begin{aligned} O_{mn} &=\langle m|\hat O|n\rangle\\ &=\langle m|\hat{\mathcal P}^{\dagger} (\hat{\mathcal P}\hat O\hat{\mathcal P}^{\dagger}) \hat{\mathcal P}|n\rangle\\ &=p_mp_n\eta_O\,O_{mn}. \end{aligned} \]

所以

\[ \boxed{p_mp_n\eta_O=-1\quad\Longrightarrow\quad O_{mn}=0.} \]

• 若 \(\ket m,\ket n\) parity 相反，even operator 的 matrix element 為零，odd operator 可以連接兩態。
• 若兩態 parity 相同，odd operator 的 matrix element 為零，even operator 才可能連接兩態。
• 「可以」只表示 symmetry 沒有強迫它為零；實際大小仍取決於 wavefunction。

這正是分離 noise 與 control 的第一種方法：讓 \(\hat O_{\rm noise}\) 具有被禁止的 parity，讓 \(\hat O_{\rm drive}\) 具有被允許的 parity。 在 half-flux Fluxonium 的對稱位能中，states 具有確定 parity；但要特別注意，\(\hat n\) 與相對於對稱中心的 phase operator 都是 odd operator， 因此相反 parity 的 \(0\leftrightarrow1\) transition 是 parity-allowed。若該 matrix element 很小，物理來源是 wavefunction／tunneling suppression，而不是把它誤稱為 parity-forbidden。

第三步：exact selection rule 如何被 symmetry-breaking perturbation 打開？

若真實裝置含有小的 symmetry-breaking term，

\[ \hat H=\hat H_0+\epsilon\hat W, \qquad [\hat W,\hat{\mathcal P}]\neq0, \]

一階 perturbation theory 給出

\[ \ket{\widetilde n} \simeq\ket n+\epsilon\sum_{k\neq n} \frac{\langle k|\hat W|n\rangle}{E_n-E_k}\ket k. \]

原本嚴格為零的 \(\langle0|\hat O|1\rangle\)，會因 opposite-symmetry admixture 變成 \(O(\epsilon)\)；對應 transition rate 通常由 matrix element 的平方開始成為 \(O(\epsilon^2)\)。因此聲稱「symmetry protected」時，還必須列出 flux detuning、junction mismatch、stray capacitance、offset charge 或 control-line asymmetry 如何破壞 symmetry，而不能只展示理想 Hamiltonian。

第四步：wavefunction overlap suppression 不等於 exact zero

在廣義座標 \(x\) 的 representation 中，

\[ O_{01}=\langle0|\hat O|1\rangle =\int dx\,\psi_0^*(x)\,\hat O(x,-i\partial_x)\,\psi_1(x). \]

若 logical wavefunctions localized 在不同 wells、不同 regions，或多維 configuration space 中近乎 disjoint support，對某些 local／smooth operator 的積分可以非常小。 這種 suppression 可以隨 barrier action、有效質量或 energy-scale ratio 呈 exponential scaling，但通常不是 exact zero：finite tunneling、disorder、mode mixing 與 operator 的 derivative structure 都可能留下 residual matrix element。

所以必須把兩種敘述分開：selection-rule forbidden 是由 symmetry algebra 得到零； overlap suppressed 是由具體 eigenfunction 積分得到小量。前者主要問 symmetry-breaking perturbation，後者主要問參數 scaling 與 fabrication variation。

多自由度或 encoded state structure 還可以把 dominant fluctuation 變成 common-mode action，使 \(\hat P_{\rm L}(\sum_j c_j\hat O_j)\hat P_{\rm L}\simeq c\hat I\)。但這只保護指定的 collective channel； local disorder 或 differential noise 若投影成 \(\sigma_z\) 或 \(\sigma_{x,y}\)，仍會造成 decoherence。

第五步：sweet spot 是 diagonal protection，不是 transition selection rule

對 parameter noise \(\lambda\)，問題一的 Hellmann–Feynman relation 給出

\[ \frac{\partial\omega_{01}}{\partial\lambda} =\frac{1}{\hbar} \left( \left\langle1\left|\frac{\partial\hat H_0}{\partial\lambda}\right|1\right\rangle - \left\langle0\left|\frac{\partial\hat H_0}{\partial\lambda}\right|0\right\rangle \right). \]

Sweet spot \(\partial_\lambda\omega_{01}=0\) 表示兩個 logical states 對該 fluctuation 具有相同的一階 diagonal response，環境較難由 frequency shift 分辨它們。 它並沒有要求 \(\langle0|\partial_\lambda\hat H_0|1\rangle=0\)，所以不能自動推出較長 \(T_1\)。反過來，off-diagonal matrix element 很小也不保證 \(O_{11}-O_{00}\) 很小。

因此必須區分三種 protection：

1. Sweet-spot protection：\(\partial\omega_{01}/\partial\lambda=0\)，主要降低低頻 parameter fluctuation 所造成的 dephasing。
2. Selection-rule protection：因 exact symmetry 與 operator parity，使特定 matrix element 嚴格為零。
3. Matrix-element suppression：transition 未被 symmetry 禁止，但因 wavefunction localization 或 overlap 很小而受到抑制。

第六步：direct transition 太弱時，control 還有哪些路徑？

保留 controllability 不等於一定要把 protected \(0\leftrightarrow1\) transition 做得很強。常見方法有：

1. 使用不同 operator／port：例如 dominant loss 透過 charge-like operator，但 control 使用 flux-like operator；或反過來利用 operator parity。
2. 使用 auxiliary level：讓 \(0\leftrightarrow a\) 與 \(1\leftrightarrow a\) 為 allowed transitions，再用 two-tone／Raman process 實現 logical rotation。
3. 暫時打破 protection：idle 時停在 protected point，gate 時施加 bias 或 parametric modulation，短暫混合 states，pulse 後回到 protected Hamiltonian。
4. 經由 virtual transitions readout：即使 \(O_{01}\) 很小，\(\ket0\)、\(\ket1\) 到 higher states 的不同 virtual coupling 仍可產生 dispersive contrast。

以 detuned auxiliary state \(\ket a\) 為例，消去 \(\ket a\) 後的有效 Raman coupling 具有

\[ \Omega_{\rm eff}\sim \frac{\Omega_{0a}\Omega_{a1}}{2\Delta_a}, \]

但 auxiliary-state population、AC Stark shift、additional decay path 與 calibration complexity 都會成為代價。暫時打破 symmetry 也不是免費的： control pulse 可能同時打開原本被禁止的 noise channel。

Dispersive readout 也可以利用 higher-state virtual transitions，而不要求直接 qubit transition 很強。在 rotating-wave、dispersive approximation 下可示意為

\[ \chi_j\simeq \sum_{kj}\frac{|g_{jk}|^2}{\omega_{kj}-\omega_r}, \qquad \chi=\chi_1-\chi_0, \qquad \omega_{ab}=\frac{E_a-E_b}{\hbar}>0\ \text{for }a>b. \]

這條路徑仍需檢查 critical photon number、Purcell decay、near-resonant higher transition 與 measurement-induced leakage。

第七步：較強 drive 不能用來掩蓋所有 controllability 問題

若 \(O_{01}^{\rm drive}\) 很小，增加 amplitude \(A\) 的確可以恢復 \(\Omega_R\)，但同一個 pulse 也會耦合到 \(\ket{k\ge2}\)。定義與 drive 最近的 leakage detuning：

\[ \Delta_{\rm leakage} =\min_{i\in\{0,1\},\,k\ge2} \left|\omega_{ki}-\omega_d\right|. \]

對 off-resonant leakage channel，最低階估計為

\[ P_{\rm leakage} \sim \sum_{i\in\{0,1\},\,k\ge2} \left| \frac{A\langle k|\hat O_{\rm drive}|i\rangle} {\hbar(\omega_{ki}-\omega_d)} \right|^2. \]

因此合理的 state structure 不只要使 \(O_{01}^{\rm noise}\) 小，還要讓 desired drive matrix element 與 undesired leakage matrix elements 形成足夠大的 hierarchy。 對直接 resonant control，通常需要 \(\Omega_R\ll\Delta_{\rm leakage}\)；實際 gate 還要加入 pulse spectrum、AC Stark shift、counter-rotating term 與 calibration error。

完整設計檢查順序

1. 指定 logical states \(\ket0,\ket1\) 與所有 dominant \(\hat O_{\rm noise}^{(\alpha)}\)。
2. 計算每個 projected operator：\(O_{11}-O_{00}\) 與 \(O_{01}\)，分別判斷 dephasing 與 transition protection。
3. 說明小 matrix element 來自 exact symmetry、wavefunction overlap，還是只來自較低 environmental \(S(\omega)\)。
4. 加入 realistic asymmetry，檢查 forbidden matrix element 對 mismatch／bias error 的 scaling。
5. 選擇 \(\hat O_{\rm drive}\)、direct 或 auxiliary transition，計算 Rabi rate。
6. 列出所有 nearby transitions，檢查 leakage、readout、reset 與 thermalization，而不是只驗證 idle coherence。

為什麼不只繼續調整 Transmon 的 \(E_J/E_C\)？

Fluxonium 可以視為對 charge／flux qubit 設計的重新安排：它不是只在原有 Transmon Hamiltonian 裡換一組數值， 而是加入非常大的電感，形成新的 global confinement、flux dependence 與 multi-well potential landscape。 與 Transmon 相比，Fluxonium 常具有較低 computational frequency、較大的 anharmonicity，並可在特定工作點得到較長 relaxation time。

Fluxonium 的三個主要元件

• Josephson junction：提供 nonlinear cosine potential 與 \(E_J\)。
• Capacitor：提供 charging energy \(E_C\)。
• Superinductor：提供很小的 inductive energy \(E_L\) 與 phase confinement。

Superinductor 必須同時具有非常大的電感、低損耗與低 parasitic capacitance。實作可能使用 Josephson-junction array、 高 kinetic-inductance 材料或其他高阻抗結構；這些實作也會帶來 array modes、製程變異與寄生耦合。

第一步：從 branch flux 寫出三個 circuit-element energies

選擇跨 small junction 的 branch flux 作為 generalized coordinate：

\[ \Phi(t)=\int^t V(t')dt', \qquad \varphi=\frac{2\pi\Phi}{\Phi_0}, \qquad \varphi_{\rm ext}=\frac{2\pi\Phi_{\rm ext}}{\Phi_0}. \]

Capacitor、Josephson junction 與 linear inductor 分別提供

\[ T_C=\frac12C_\Sigma\dot\Phi^2, \qquad U_J=-E_J\cos\varphi, \qquad U_L=\frac{(\Phi-\Phi_{\rm ext})^2}{2L}. \]

這裡 external flux 的位置取決於 gauge 與 coordinate choice；物理內容是 fluxoid constraint 使 small-junction phase、superinductor phase 與 external flux 不再彼此獨立。

第二步：由 Lagrangian 做 Legendre transform

Lagrangian 為 kinetic energy 減 potential energy：

\[ \mathcal L =\frac12C_\Sigma\dot\Phi^2 +E_J\cos\!\left(\frac{2\pi\Phi}{\Phi_0}\right) -\frac{(\Phi-\Phi_{\rm ext})^2}{2L}. \]

與 \(\Phi\) 共軛的 canonical charge 是

\[ Q=\frac{\partial\mathcal L}{\partial\dot\Phi} =C_\Sigma\dot\Phi, \qquad n=\frac{Q}{2e}, \qquad [\hat\Phi,\hat Q]=i\hbar \Longleftrightarrow [\hat\varphi,\hat n]=i. \]

做 Legendre transform：

\[ \begin{aligned} \hat H &=\hat Q\dot{\hat\Phi}-\mathcal L\\ &=\frac{\hat Q^2}{2C_\Sigma} -E_J\cos\hat\varphi +\frac{(\hat\Phi-\Phi_{\rm ext})^2}{2L}\\ &=4E_C\hat n^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2, \end{aligned} \]

其中

\[ E_C=\frac{e^2}{2C_\Sigma}, \qquad E_L=\frac{1}{L}\left(\frac{\Phi_0}{2\pi}\right)^2. \]

所以問題三最直接的答案是：superinductor 增加 quadratic inductive-energy term \(\tfrac12E_L(\hat\varphi-\varphi_{\rm ext})^2\)。但更重要的反直覺是 \(L\) 越大，\(E_L\) 越小；superinductor 提供的是曲率很小、範圍卻遍及許多 cosine wells 的 global confinement。

第三步：兩種常見 gauge 寫法是同一個 Hamiltonian

把新 coordinate 定義成 \(\hat\theta=\hat\varphi-\varphi_{\rm ext}\)，可改寫為

\[ \hat H =4E_C\hat n^2 +\frac12E_L\hat\theta^2 -E_J\cos(\hat\theta+\varphi_{\rm ext}). \]

因此文獻中常見的兩種形式

\[ -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2 \qquad\text{與}\qquad -E_J\cos(\hat\theta+\varphi_{\rm ext}) +\frac12E_L\hat\theta^2 \]

只差 coordinate shift。推導與數值計算必須從頭到尾選定同一個 gauge，不能同時在 cosine 與 parabola 重複加入 external-flux shift。

第四步：為什麼 inductive shunt 會改變 phase coordinate 與 offset charge？

Transmon-like circuit 寫成

\[ \hat H_{\rm T}=4E_C(\hat n-n_g)^2-E_J\cos\hat\varphi. \]

沒有 inductive term 時，potential 是 \(2\pi\)-periodic，phase 是 compact coordinate，\(n_g\) 會進入 wavefunction 的 boundary condition。 加入 inductor 後，quadratic term 使 \(\varphi\) 成為定義在整條 real line 的 extended coordinate，\(\varphi\to\varphi+2\pi\) 會移到另一個 cosine well， 不再是同一個 configuration point。

若暫時保留 offset charge，

\[ \hat H=4E_C(\hat n-n_g)^2+U(\hat\varphi), \qquad \hat U_g=e^{in_g\hat\varphi}, \qquad \hat U_g^\dagger(\hat n-n_g)\hat U_g=\hat n. \]

因 \(\varphi\in\mathbb R\)，這個 unitary transformation 不會引入 compact-coordinate 的 twisted boundary condition，static \(n_g\) 可由 low-energy Hamiltonian 消去。 這是 Fluxonium「free of charge offsets」的結構性來源；它不表示所有 capacitive loss 或 \(\hat n\)-coupled finite-frequency noise 都消失。

第五步：Josephson-junction array 為什麼可以近似成 superinductor？

實作上常用 \(N\) 個較大的 Josephson junction 串聯。令每顆 array junction 的 phase drop 為 \(\theta_j\)，且 \(E_{J,a}\gg E_{C,a}\)，使每顆 junction 很少發生 quantum phase slip。對小 \(\theta_j\) 展開：

\[ -E_{J,a}\cos\theta_j \simeq-E_{J,a}+\frac12E_{J,a}\theta_j^2. \]

若 collective phase drop 為 \(\Theta=\sum_{j=1}^{N}\theta_j\)，最低能量分布近似 \(\theta_j=\Theta/N\)，因此

\[ U_{\rm array} \simeq\sum_{j=1}^{N}\frac12E_{J,a}\left(\frac{\Theta}{N}\right)^2 =\frac12\frac{E_{J,a}}{N}\Theta^2 \equiv\frac12E_L\Theta^2. \]

所以在最簡單且均勻的 array approximation 中，

\[ E_L\simeq\frac{E_{J,a}}{N}, \qquad L_{\rm array}\simeq N\frac{(\Phi_0/2\pi)^2}{E_{J,a}}. \]

串聯 junction 數增加，總電感增大而 effective \(E_L\) 下降。但 \(N\) 不能任意增加：junction capacitance、island-to-ground capacitance、parameter disorder 與 array length 會共同產生一整組 internal plasma／array modes。

第六步：新增 quadratic term 後，位能在數學上發生什麼？

Effective potential 是

\[ U(\varphi) =-E_J\cos\varphi +\frac12E_L(\varphi-\varphi_{\rm ext})^2. \]

Local minima \(\varphi_m\) 滿足

\[ E_J\sin\varphi_m +E_L(\varphi_m-\varphi_{\rm ext})=0, \]

而每個 minimum 附近的 curvature 與 small-oscillation energy scale 約為

\[ U''(\varphi_m)=E_J\cos\varphi_m+E_L, \qquad \hbar\omega_{p,m}\simeq \sqrt{8E_C\big(E_J\cos\varphi_m+E_L\big)}. \]

當 \(E_J/E_L\) 很大時，cosine 產生許多 local wells，緩慢變化的 parabola 決定各 wells 的 global energy ordering。忽略 local zero-point correction 時，

\[ U_m\approx-E_J+\frac12E_L(2\pi m-\varphi_{\rm ext})^2. \]

External flux 因而可以重新排列不同 fluxoid-like wells；在 \(\varphi_{\rm ext}=\pi\) 時，相鄰 wells 可成為對稱的 degenerate pair。 \(E_C\) 與 barrier 決定 inter-well tunneling，這正是問題四要繼續追蹤的 wavefunction engineering。

第七步：什麼條件下 array 才能只留下單一 \(E_L\)？

真實 array 並不是沒有 dynamics 的理想線圈。保留所有 array-island fluxes 時，Hamiltonian 具有 capacitance matrix 與多個 Josephson cosine：

\[ \hat H_{\rm array} =\frac12\hat{\mathbf Q}^{\mathsf T}\mathbf C^{-1}\hat{\mathbf Q} -\sum_{j=1}^{N}E_{J,a}^{(j)}\cos\hat\theta_j. \]

把它約化成單一 quadratic term，至少需要：

1. 每顆 array junction 位於 \(E_{J,a}/E_{C,a}\gg1\) 的近線性 regime，quantum phase-slip rate 足夠低。
2. 最低 internal array-mode frequency 高於 qubit、control 與 readout 的 relevant frequency range。
3. Array modes 與 qubit collective coordinate 的 hybridization 足夠弱，或已明確納入 multimode Hamiltonian。
4. Ground capacitance、junction disorder、dielectric loss 與 quasiparticle loss 不會把 superinductor 變成新的 dominant noise source。

一旦這些條件失效，「加入 superinductor」就不再只是多一個 \(\tfrac12E_L\varphi^2\)：還會多出 parasitic modes、avoided crossings、additional decay／leakage paths 與 fabrication sensitivity。

Superinductor 為什麼不是直接換成大 resistor？

Superinductance 的目標是零 DC resistance、低 dissipation，並在工作頻帶提供高於 Cooper-pair resistance quantum \(R_Q=h/(2e)^2\approx6.45\,\mathrm{k}\Omega\) 的 impedance，同時把 self-resonance 推到工作頻帶之外。 大 resistor 雖然也可能有高 impedance，卻帶來實部 admittance、Johnson–Nyquist noise 與能量耗散，不能提供同樣的 coherent Hamiltonian term。

• \(\hat n\)：與 phase \(\hat\varphi\) 共軛的 Cooper-pair number operator。
• \(\hat\varphi\)：small junction 的 phase difference。
• \(\varphi_{\rm ext}=2\pi\Phi_{\rm ext}/\Phi_0\)：external flux 所造成的 phase offset。
• \(E_C\)：charging energy，控制 charge kinetic term。
• \(E_J\)：Josephson energy，控制 cosine wells 的深度。
• \(E_L\)：inductive energy，由 superinductor 決定。

Layout-to-Hamiltonian mapping

Chip design 不是只畫電路，而是直接雕刻 Hamiltonian 的參數空間；同時，製程變異也會直接變成 \(E_J,E_C,E_L\) spread、 frequency offset 與 yield 問題。

問題三已由 circuit 得到 Fluxonium Hamiltonian。問題四不再只問「多了哪一項」，而是追蹤每個 energy scale 如何改變 eigenvalue problem：

\[ \hat H =4E_C\hat n^2 -E_J\cos\hat\varphi +\frac12E_L(\hat\varphi-\varphi_{\rm ext})^2, \qquad [\hat\varphi,\hat n]=i. \]

第一步：把 Hamiltonian 寫成 phase-basis Schrödinger equation

在 phase basis 中 \(\hat n=-i\partial_\varphi\)，因此

\[ \left[ -4E_C\frac{\partial^2}{\partial\varphi^2} +U(\varphi) \right]\psi_k(\varphi) =E_k\psi_k(\varphi), \]

其中

\[ U(\varphi)=-E_J\cos\varphi +\frac12E_L(\varphi-\varphi_{\rm ext})^2. \]

與普通一維 Schrödinger equation 比較，phase particle 的 effective mass 為

\[ m_\varphi=\frac{\hbar^2}{8E_C}. \]

所以 \(E_C\) 不只是「charging-energy 數字」：它決定 phase particle 的量子動能與有效質量；\(E_J,E_L,\varphi_{\rm ext}\) 則共同決定它所看見的 potential landscape。

三個能量尺度的競爭

Fluxonium 的物理特色不是某一個單獨參數，而是 \(E_C,E_J,E_L\) 的競爭。調整它們會同時改變 wavefunction localization、 tunneling、transition frequency、anharmonicity、operator matrix element 與 leakage spectrum。

第二步：先在單一 well 內讀出 localization 與 zero-point fluctuation

令 \(\varphi_m\) 為某一 local minimum，並定義

\[ K_m=U''(\varphi_m)=E_J\cos\varphi_m+E_L. \]

在 \(\varphi_m\) 附近做 harmonic expansion，得到

\[ \hbar\omega_{p,m}\simeq\sqrt{8E_CK_m}, \qquad \varphi_{{\rm zpf},m}=\left(\frac{2E_C}{K_m}\right)^{1/4}, \qquad n_{{\rm zpf},m}=\left(\frac{K_m}{32E_C}\right)^{1/4}. \]

這三個量直接給出參數直覺：

• 增大 \(E_C\)：phase particle 變輕，\(\psi(\varphi)\) 變寬，inter-well tunneling 增強；charging kinetic term 在能譜中的作用也更強。
• 增大 \(E_J\)：cosine wells 變深、barrier 變高，local curvature 增大，但不同 wells 間的 tunneling 通常下降。
• 減小 \(E_L\)：global parabola 變平，相鄰 fluxoid-like wells 的 energy detuning 下降，更多 wells 進入低能 spectrum。
• 改變 \(\varphi_{\rm ext}\)：不直接改變 cosine barrier，而是平移 parabola，重新排列不同 wells 的相對能量。

Deep-well regime 的 tunnel splitting 可概略寫成 \(\Delta\sim\hbar\omega_p e^{-S}\)，其中 semiclassical action \(S\) 隨 barrier height 與 effective mass 增加。 因此較大 \(E_J/E_C\) 通常使 states 更 localized、\(\Delta\) 更小；精確值仍須由完整 potential 求解。

第三步：external flux 如何選擇不同 wells？

在 \(E_J/E_L\gg1\) 時，local minima 接近 \(\varphi_m\simeq2\pi m\)，其 localized-state energy 可先估為

\[ E_m^{\rm loc} \approx -E_J +\frac12E_L(2\pi m-\varphi_{\rm ext})^2 +\frac12\hbar\omega_{p,m}. \]

對 external flux 的斜率對應 persistent current：

\[ I_m=-\frac{\partial E_m^{\rm loc}}{\partial\Phi_{\rm ext}} \approx \frac{2\pi E_L}{\Phi_0}(2\pi m-\varphi_{\rm ext}). \]

當兩個 localized wells 的 energies 接近時，tunneling 使 crossing 變成 avoided crossing。External flux 控制的是 well detuning； \(E_J/E_C\) 主要控制 tunnel gap；\(E_L\) 則同時控制 well detuning 的尺度與 persistent-current slope。

Half-flux potential、parity 與 wavefunction

在 \(\varphi_{\rm ext}=\pi\) 令 \(x=\varphi-\pi\)，Hamiltonian 變成

\[ \hat H_{1/2} =-4E_C\frac{\partial^2}{\partial x^2} +E_J\cos x+\frac12E_Lx^2, \qquad U(x)=U(-x). \]

Potential 具有 exact reflection symmetry， eigenstates 因而可具有確定 parity：ground state \(\ket0\) 為偶、first excited state \(\ket1\) 為奇。 Bonding／anti-bonding 直覺可以輔助理解，但 selection rule 的嚴格推論來自 symmetry 與 operator parity。

若只保留左右兩個 localized states \(\ket L,\ket R\)，low-energy Hamiltonian 可寫成

\[ \hat H_{\rm 2w} =-\frac12\left(\Delta\hat\sigma_x+\epsilon\hat\sigma_z\right), \qquad \epsilon\simeq2I_p\,\delta\Phi, \qquad \delta\Phi=\Phi_{\rm ext}-\frac{\Phi_0}{2}. \]

其 qubit splitting 與 half-flux eigenstates 為

\[ \hbar\omega_{01}=\sqrt{\Delta^2+\epsilon^2}, \qquad \ket0\simeq\frac{\ket L+\ket R}{\sqrt2}, \qquad \ket1\simeq\frac{\ket L-\ket R}{\sqrt2}. \]

偏離 half flux 時 \(|\epsilon|\) 增加，eigenstates 逐漸 localized；在 half flux 則形成 even／odd superpositions。 這個 two-well model 提供直覺，但 higher plasmon states、其他 wells 與真實 matrix elements仍須由完整 Hamiltonian 計算。

• \(\hat n\) 與 \(x=\varphi-\pi\) 都是 odd operator，原則上可以連接偶、奇態。
• 若 \(\langle0|\hat n|1\rangle\) 或 \(\langle0|x|1\rangle\) 很小，不能只說「parity forbidden」；必須檢查 device-specific wavefunction overlap。
• \(\hat n^2\)、\(x^2\) 等 even operator 在相反 parity states 間可由 selection rule 嚴格禁止。

Exact parity 來自 half-flux 處的 \(U(x)=U(-x)\)。\(E_J\gg E_L\) 與 \(E_J\gg E_C\) 是 deep-well／bonding 圖像常用的條件， 但不是 parity symmetry 成立的必要條件。偏離 half flux 時 symmetry 與一階 flux sensitivity 會恢復，其大小需依裝置參數計算。

第四步：從 symmetry 與 two-well model 推導 flux sweet spot

External-flux coupling operator 是

\[ \frac{\partial\hat H}{\partial\Phi_{\rm ext}} =-\frac{2\pi E_L}{\Phi_0} (\hat\varphi-\varphi_{\rm ext}). \]

在 half flux，\(x=\varphi-\pi\) 是 odd operator，而任何 parity eigenstate 都滿足 \(\langle n|x|n\rangle=0\)。由 Hellmann–Feynman theorem，

\[ \left.\frac{\partial E_n}{\partial\Phi_{\rm ext}}\right|_{\Phi_0/2}=0, \qquad \boxed{ \left.\frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}} \right|_{\Phi_0/2}=0.} \]

Two-well model 給出同一結論：

\[ \frac{\partial\omega_{01}}{\partial\Phi_{\rm ext}} =\frac{1}{\hbar} \frac{\epsilon}{\sqrt{\Delta^2+\epsilon^2}} \frac{\partial\epsilon}{\partial\Phi_{\rm ext}} \xrightarrow[\epsilon\to0]{}0. \]

在此 operating point，一階 flux sensitivity 消失，dephasing 由 higher-order term 主導。但 sweet spot 不自動保證 relaxation protection； relaxation 仍須檢查相應 operator 的 matrix element 與 \(S_F(\omega_{01})\)。

而且在 two-well approximation 中

\[ \left.\frac{\partial^2\omega_{01}}{\partial\Phi_{\rm ext}^2} \right|_{\Phi_0/2} \simeq\frac{(2I_p)^2}{\hbar\Delta}. \]

所以把 tunnel gap \(\Delta\) 做得非常小，雖可降低 \(\omega_{01}\) 與部分 transition matrix element，卻會增加 sweet spot 的二階 curvature， 並加重 thermal initialization 與慢控制問題。Protection 仍然是多目標 trade-off。

第五步：逐一檢查 noise operator，而不是籠統說「Fluxonium 被保護」

Fluxonium 的低能態可分布在不同 phase wells。若兩個 states 的 wavefunction overlap 很小，某些 operator 的 transition matrix element 也會變小，使 \(\ket1\rightarrow\ket0\) 經該 noise channel 的機率降低。這是 wavefunction engineering；它與 exact parity selection rule 必須分開命名。

對 Hamiltonian \(\hat H=4E_C\hat n^2+U(\hat\varphi)\)，commutator 提供一條非常有用的 exact relation：

\[ [\hat H,\hat\varphi]=-8iE_C\hat n, \qquad \left|\langle0|\hat n|1\rangle\right| =\frac{\hbar\omega_{01}}{8E_C} \left|\langle0|\hat\varphi|1\rangle\right|. \]

這說明 half-flux 的 \(0\leftrightarrow1\) transition 即使具有不小的 phase dipole，低 \(\omega_{01}\) 仍會壓低 charge matrix element。 因此 capacitive／dielectric relaxation 可以受到抑制，而 flux-like operator 的 relaxation 不會自動得到同樣的 suppression；仍須乘上各自的 \(S_{\rm noise}(\omega_{01})\)。

第六步：protection 必須與 control、readout、leakage 一起驗證

對 capacitive 與 flux-like drive，直接 Rabi rates 分別具有

\[ \Omega_R^{(n)}\propto A_n\left|\langle0|\hat n|1\rangle\right|, \qquad \Omega_R^{(\varphi)}\propto A_\varphi\left|\langle0|x|1\rangle\right|. \]

Charge matrix element 小會讓 capacitive direct gate 變慢；Fluxonium 的 higher transition，例如 \(1\leftrightarrow2\)，可能具有大得多的 charge matrix element， 因此 control／readout 可利用 operator hierarchy、auxiliary level 或不同 coupling port。代價是 pulse sequence、frequency allocation 與 leakage analysis 更複雜。

至少要同時計算

\[ \alpha=\omega_{12}-\omega_{01}, \qquad \Delta_{\rm leakage} =\min_{i\in\{0,1\},\,k\ge2} |\omega_{ki}-\omega_d|, \qquad \Omega_R\ll\Delta_{\rm leakage}. \]

Fluxonium 常有較大 anharmonicity，但同時具有 rich multilevel spectrum；「anharmonicity 大」不能取代逐條檢查 drive matrix element、near resonance、 array modes 與 multiphoton leakage。

Protected spectrum 的直接代價：低頻與 thermal population

Fluxonium 的 protected regime 常伴隨較低 computational frequency。以 \(f_{01}=500\,\mathrm{MHz}\) 為例， \(hf/k_B\approx24\,\mathrm{mK}\)：在 \(T=50\,\mathrm{mK}\) 時，\(p_1/p_0\approx e^{-24/50}\approx0.62\)， \(p_1\approx0.38\)；在 \(T=20\,\mathrm{mK}\) 時，\(p_1/p_0\approx0.30\)，\(p_1\approx0.23\)。

另一個常用尺度是 \(f_{01}=800\,\mathrm{MHz}\)，此時 \(hf/k_B\approx38\,\mathrm{mK}\)，在 \(20\,\mathrm{mK}\) 時 \(p_1\approx0.13\)（只作 thermal-order estimate；實際值仍受非平衡 photons 與有效溫度影響）。低頻 transition 使 active reset、 thermalization 與額外 readout transition 成為設計問題。

完整驗證流程：由參數走到 protection claim

1. 指定 \(E_C,E_J,E_L,\Phi_{\rm ext}\)，求出 \(U(\varphi)\)、local minima 與 barrier。
2. 對完整 Hamiltonian diagonalize，畫出 \(E_k(\Phi)\) 與 \(\psi_k(\varphi)\)，不能只依賴 two-well 圖像。
3. 計算 \(\partial_\lambda\omega_{01}\) 與二階 curvature，確認 sweet spot 對哪個 parameter 成立。
4. 對每個 dominant noise channel 計算 \(\langle0|\hat O_{\rm noise}|1\rangle\) 與對應 \(S_{\rm noise}(\omega_{01})\)。
5. 計算 drive／readout operator 的 matrix elements、Rabi rate、dispersive path 與 \(\Delta_{\rm leakage}\)。
6. 把 \(\omega_{01}\) 與 refrigerator／effective temperature 比較，加入 thermal population、reset 與 bandwidth constraint。

前四問已用 Fluxonium 建立從 noise source、operator、Hamiltonian 到 eigenstate／matrix element 的完整鏈。第五問把 Fluxonium 視為第一個完整案例，並將方向反過來： 若先指定 dominant error operator，我們需要改變參數、增加自由度、創造新的 symmetry，還是改造 Josephson element 本身？

Protected-qubit 的四個不可省略條件

這四項不是四個可任選的優點，而是一個 protected-qubit proposal 必須同時回答的 design specification：

1. Dephasing：哪一個 fluctuation parameter \(\lambda\) 被鈍化？只在單一 sweet spot，還是在一段 parameter range？
2. Relaxation：哪一個 \(\hat O_{\rm noise}\) 的 matrix element 被壓低？是 exact symmetry 還是 approximate overlap suppression？
3. Control：若 noise coupling 被壓低，使用哪一個不同的 \(\hat O_{\rm drive}\)、哪一條 transition 或哪一個 auxiliary level 控制？
4. Leakage：保護 logical states 時，有沒有同時製造鄰近 parasitic states、dense spectrum 或新的 leakage path？

從 error target 反推 Hamiltonian 的六個步驟

1. 指定 error：先寫出 dominant fluctuation \(F(t)\)、coupling operator \(\hat O_{\rm noise}\) 與 relevant frequency range；低頻 dephasing 與 \(\omega_{01}\) 附近的 relaxation 要分開。
2. 指定 state structure：決定要利用 sweet spot、exact symmetry、wavefunction localization、near degeneracy，或多自由度的 disjoint support。
3. 合成 candidate Hamiltonian：選擇要改變的 parameter regime、增加的 energy term、額外自由度、Josephson Fourier component 或 material-dependent element。
4. 驗證 protection：求 spectrum／eigenstates，計算 \(\partial\omega_{01}/\partial\lambda\) 與 \(\langle0|\hat O_{\rm noise}|1\rangle\)，不能只靠位能圖宣稱「受到保護」。
5. 驗證可操作性：檢查 \(\hat O_{\rm drive}\)、Rabi rate、readout path、auxiliary transition、\(\Delta_{\rm leakage}\) 與 reset。
6. 映射回硬體：把每個 Hamiltonian term 映射到 capacitor、junction、superinductor、loop 或材料，再列出新引入的 loss、disorder、parasitic mode 與 fabrication sensitivity。

進階 Topic 名詞巡禮：Hamiltonian engineering 的三條延伸路線

這裡的教學目的就是建立名詞與概念入口：學生不必在本週完成每種架構的完整推導，但必須看見不同路線分別在 Hamiltonian 的哪個層次工作， 日後讀到新 qubit proposal 時，能用同一組四項判準自行拆解。

1. \(0-\pi\) qubit：增加自由度與 disconnected state support

Fluxonium 主要仍可視為有效單自由度問題 \(H=H(\varphi,n)\)；\(0-\pi\) circuit 把問題提升為多模態 Hamiltonian， 例如 \(H=H(\phi,\theta)\)。實際 circuit 使用兩個 Josephson junction、兩個大電感與兩個大電容，目的不是再找一個普通 sweet spot， 而是讓 logical states 在 generalized-coordinate space 中具有近乎不相交的支撐，同時使 logical energies 接近簡併：

• Dephasing：near degeneracy 降低部分 parameter sensitivity。
• Relaxation：logical-state support 分離可壓低 local-operator matrix elements。
• Control：必須設計能跨越 protected subspaces 的 drive path，常需 auxiliary states 或 symmetry-breaking control。
• Leakage：multimode spectrum 與寄生模態必須一起檢查。
• 工程挑戰：large superinductance、low stray capacitance、元件匹配與 disorder-induced mode coupling。

自學問題：為什麼增加一個自由度，可能做到單一 phase coordinate 難以同時做到的保護？

2. \(\cos(2\varphi)\) qubit：直接工程化 Josephson potential

一般 Josephson junction 提供 \(-E_J\cos\varphi\)。更一般地，可把 current-phase relation 的 Fourier components 寫成：

\(\cos\varphi\) 對應 single-Cooper-pair tunneling；有效 \(\cos(2\varphi)\) 可對應 two-Cooper-pair process，建立不同週期性與 Cooper-pair parity sectors，使某些 charge-induced transitions 受到 symmetry 抑制。這把 Hamiltonian engineering 從「加入 \(L\)、\(C\)」 提升成「直接選擇 Josephson potential 的 Fourier content」。

• Dephasing：檢查 logical splitting 對 offset charge、flux 與元件 asymmetry 的敏感度。
• Relaxation：利用 parity sectors 壓低會改變受保護 parity 的 transition。
• Control：控制脈衝若必須暫時打破 parity，要檢查保護與 gate speed 的交換。
• Leakage：新的 periodicity 會重排 spectrum；必須確認非計算態與 unwanted Fourier components。
• 工程挑戰：如何實現占優勢的 \(E_2\)，並抑制殘餘 \(E_1\cos\varphi\) 與 fabrication asymmetry。

自學問題：我們是否一定只能接受自然 junction 給出的 \(-E_J\cos\varphi\)，還是可以工程化它的 Fourier components？

3. Ferro-transmon：改造材料，也是在改造 Hamiltonian

Ferromagnetic Josephson-junction 路線不是單純增加 linear circuit element，而是讓 Josephson energy 或更一般的 current-phase relation 依賴磁性狀態。Hybrid ferro-transmon 可利用 hysteretic behavior，使 qubit frequency 透過 magnetic pulse 進行近似 digital tuning； 但磁性自由度也會同時帶來 magnetic fluctuation 與 domain dynamics 等新 noise source。

• Dephasing：digital／hysteretic operating point 可能降低某些連續 bias sensitivity，但必須量測 magnetic-state stability。
• Relaxation：不會僅因使用 ferroic material 自動得到 matrix-element protection；仍須指出 \(\hat O_{\rm noise}\)。
• Control：磁性狀態提供新的 nonvolatile tuning knob。
• Leakage：需檢查額外 material excitations 與 modified current-phase relation 所引入的 levels。
• 工程挑戰：magnetic fluctuation、hysteresis reproducibility、材料介面與製程整合。

核心提醒：Hamiltonian engineering 不一定等於 protection engineering；新的 control knob 也可能同時帶來新的 noise operator。

四項判準 × 三條路線

Fluxonium 自身的設計取捨：protected spectrum 不是免費的

以四項判準回頭檢查 Fluxonium，必須同時列出工作頻率、anharmonicity、sweet spot、dominant noise operator、gate primitive、 readout contrast、reset、fabrication complexity、coupler compatibility 與 frequency allocation。只比較最佳 \(T_1\) 會忽略 control stack 與 yield。

課堂整合與 Exit ticket

1. 本週最重要的一條 equation／model 是什麼？它的適用近似為何？
2. 哪一個看似局部的參數，實際會跨越 device、control 與 system 三層？
3. 若只能做一個新實驗，什麼結果最能改變目前的 architecture decision？
4. 從 \(0-\pi\)、\(\cos2\varphi\)、ferro-transmon 中選一個名詞，寫下你下一步需要查清楚的 Hamiltonian term、symmetry 與 dominant noise。

跨層驗收：重畫課程開始時的「抽象模型 → 真實硬體」訊息流圖。兩張圖的差異就是本週學習的證據。 完成條件是留下可檢查的推導、關係圖、比較表或 provisional decision，並能說明其物理假設與工程代價。

下週預告：W4 將從 driven two-level system 出發，推到 rotating frame 與 RWA。學生要帶著本週建立的 qubit Hamiltonian、 matrix element 與 leakage-gap 直覺進入 W4，才能理解 pulse control 為何受裝置能譜限制。

核心與延伸閱讀

• Fluxonium: Single Cooper-Pair Circuit Free of Charge OffsetsV. E. Manucharyan et al., Science 326, 113–116 (2009).Fluxonium 與 superinductor 的原始論文。
• High Fidelity Two-Qubit Gates on Fluxoniums Using a Tunable CouplerI. N. Moskalenko et al., npj Quantum Information 8, 130 (2022).Fluxonium 控制、reset、耦合與多量子位元工程。
• Introduction to Quantum Electromagnetic CircuitsU. Vool and M. H. Devoret, International Journal of Circuit Theory and Applications 45, 897–934 (2017).從電路圖、Lagrangian 到量子 Hamiltonian 的核心讀物。
• Engineering High-Coherence Superconducting QubitsI. Siddiqi, Nature Reviews Materials 6, 875–891 (2021).不同 qubit architecture、材料損耗與 coherence 的設計 trade-off。