tecktom.com

2準位系を切り出す:共通例題②

2026年10月17日10分で読めますAI生成教科書『虎の巻』第2部 第5章 回路の量子化:LC 回路と Josephson 接合
目次
  1. Cooper 対の箱:静的な2準位系
  2. transmon:マイクロ波で駆動する2準位系
  3. 第3準位への漏れ
  4. 木下の層

AI 生成・木下が未勉強

量子力学の様々な構成方法 第2部 第5章 §2.5.5 ★★★ Cooper 対の箱を n=0,1n=0,1 へ、transmon を最低の2準位へ射影すると共通例題②の Hamiltonian になり、Rabi 振動 (Ω/Ω′)2sin⁡2(Ω′t/2)(\Omega/\Omega')^2\sin^2(\Omega't/2) が回路の言葉で出る

Cooper 対の箱:静的な2準位系

EJ≪ECE_{\mathrm J}\ll E_{\mathrm C} で ngn_{\mathrm g} が 12\frac12 に近いとき、電荷の状態 ∣0⟩| 0 \rangle と ∣1⟩| 1 \rangle のエネルギー 4ECng24E_{\mathrm C}n_{\mathrm g}^2 と 4EC(1−ng)24E_{\mathrm C}(1-n_{\mathrm g})^2 は近く、ほかの電荷の状態は ∼4EC\sim4E_{\mathrm C} 以上離れている。(2.5.20) を {∣1⟩,∣0⟩}\{| 1 \rangle ,| 0 \rangle \} の張る部分空間へ射影し、∣e⟩=∣1⟩| e \rangle =| 1 \rangle、∣g⟩=∣0⟩| g \rangle =| 0 \rangle(σz∣e⟩=+∣e⟩\sigma_z| e \rangle =+| e \rangle)と置くと、定数を除いて

H^box=4EC(1−2ng)2σz−EJ2σx=−ℏΔ2σz−ℏΩ2σx,ℏΩ=EJ,ℏΔ=−4EC(1−2ng)(2.5.26)\hat H_{\mathrm{box}}=\frac{4E_{\mathrm C}(1-2n_{\mathrm g})}{2}\sigma_z-\frac{E_{\mathrm J}}{2}\sigma_x =-\frac{\hbar\Delta}{2}\sigma_z-\frac{\hbar\Omega}{2}\sigma_x,\qquad \hbar\Omega=E_{\mathrm J},\quad \hbar\Delta=-4E_{\mathrm C}(1-2n_{\mathrm g})\tag{2.5.26}

となる。2つの電荷の状態の差 4EC(1−2ng)4E_{\mathrm C}(1-2n_{\mathrm g}) は 4EC(1−ng)2−4ECng24E_{\mathrm C}(1-n_{\mathrm g})^2-4E_{\mathrm C}n_{\mathrm g}^2 から出る。

これは共通例題②の Hamiltonian −(ℏΔ/2)σz−(ℏΩ/2)σx-(\hbar\Delta/2)\sigma_z-(\hbar\Omega/2)\sigma_x と同じ形で、ここでは Ω\Omega が外から当てた駆動でなく接合のトンネル、Δ\Delta がゲート電圧で決まる電荷のエネルギーの差である。ngn_{\mathrm g} を 12\frac12 から遠い値(基底状態がほぼ ∣0⟩| 0 \rangle)から急に ngn_{\mathrm g} へ切り替えて時間 tt だけ置くと、第1部第2章・第2部第3章で出した式のとおり

P∣1⟩(t)=(ΩΩ′)2sin⁡2Ω′t2,Ω′=Ω2+Δ2(2.5.27)P_{| 1 \rangle }(t)=\Bigl(\frac{\Omega}{\Omega'}\Bigr)^2\sin^2\frac{\Omega't}{2},\qquad \Omega'=\sqrt{\Omega^2+\Delta^2}\tag{2.5.27}

で島に Cooper 対が1個入っている確率が振動する。縮退点 ng=12n_{\mathrm g}=\frac12 では Δ=0\Delta=0 で、周期 h/EJh/E_{\mathrm J} で完全に入れ替わる。Nakamura・Pashkin・Tsai は、ゲートに短い電圧パルスを加えてこの振動を観測し、固体の電子回路で2準位の量子状態を時間領域で制御した最初の実験とした(Y. Nakamura, Yu. A. Pashkin and J. S. Tsai (1999), Nature 398, 786)。

近似の裏付け

型 A(射影) | 小さいもの 電荷の状態 {∣0⟩,∣1⟩}\{| 0 \rangle ,| 1 \rangle \} の張る部分空間(ほかの電荷の状態までの距離 ∼4EC\sim4E_{\mathrm C} が EJE_{\mathrm J} より十分大きい) | 証明 条件付き(射影した2準位の固有値の誤差は2次の摂動で O(EJ2/EC)\mathcal{O}(E_{\mathrm J}^2/E_{\mathrm C})。図 2.5.2 左の EJ/EC=1E_{\mathrm J}/E_{\mathrm C}=1 では第3準位が離れているが、誤差の厳密な上界は、Josephson の項の大きさ(ノルム ≤EJ\le E_{\mathrm J})が2準位とほかの電荷の状態の隔たりの半分より小さい範囲で、Schrieffer–Wolff 変換が定義できて有効ハミルトニアンの摂動級数が収束することから得られる:S. Bravyi, D. P. DiVincenzo, D. Loss (2011), Ann. Phys. 326, 2793–2826 の第3節) | 行き先 近似の数理 第4章

transmon:マイクロ波で駆動する2準位系

transmon の最低の2準位 ∣g⟩=∣k=0⟩| g \rangle =| k=0 \rangle、∣e⟩=∣k=1⟩| e \rangle =| k=1 \rangle の間隔は ℏωq=E01\hbar\omega_{\mathrm q}=E_{01} である。ゲートに Vg=Vdc+V0cos⁡ωdtV_{\mathrm g}=V_{\mathrm{dc}}+V_0\cos\omega_{\mathrm d}t を加えると、(2.5.18) の (n^−ng)2(\hat n-n_{\mathrm g})^2 の交差項から

H^d(t)=2e CgCΣ V0cos⁡(ωdt) n^(2.5.28)\hat H_{\mathrm d}(t)=2e\,\frac{C_{\mathrm g}}{C_\Sigma}\,V_0\cos(\omega_{\mathrm d}t)\,\hat n\tag{2.5.28}

が出る(8EC/(2e)=2e/CΣ8E_{\mathrm C}/(2e)=2e/C_\Sigma)。n^\hat n の行列要素は (2.5.22) により ⟨e∣n^∣g⟩≈(i/2)(EJ/8EC)1/4\langle e | \hat n | g \rangle \approx(\mathrm{i}/\sqrt2)(E_{\mathrm J}/8E_{\mathrm C})^{1/4} で、純虚数である。∣e⟩| e \rangle の位相を選び直して行列要素を実数 n01=∣⟨e∣n^∣g⟩∣n_{01}=\lvert \langle e | \hat n | g \rangle \rvert にし、2準位へ射影するとH^=(ℏωq/2)σz+ℏΩcos⁡(ωdt)σx\hat H=(\hbar\omega_{\mathrm q}/2)\sigma_z+\hbar\Omega\cos(\omega_{\mathrm d}t)\sigma_x、ℏΩ=2e(Cg/CΣ)V0n01\hbar\Omega=2e(C_{\mathrm g}/C_\Sigma)V_0n_{01} である(n^\hat n の対角要素は σz\sigma_z の係数を少しずらすだけなので落とした)。eiωdtσz/2\mathrm{e}^{\mathrm{i}\omega_{\mathrm d}t\sigma_z/2} で回る座標へ移り、e±2iωdt\mathrm{e}^{\pm2\mathrm{i}\omega_{\mathrm d}t} で振動する項を捨てると(回転波近似)

H^rot=−ℏΔ2σz+ℏΩ2σx,Δ=ωd−ωq(2.5.29)\hat H_{\mathrm{rot}}=-\frac{\hbar\Delta}{2}\sigma_z+\frac{\hbar\Omega}{2}\sigma_x,\qquad \Delta=\omega_{\mathrm d}-\omega_{\mathrm q}\tag{2.5.29}

を得る。σx\sigma_x の係数の符号は駆動の位相 π\pi のずれで反転でき、(2.5.27) は Ω2\Omega^2 にしか依らないので、共通例題②の正答がそのまま当てはまる。この節の行列要素を EJ/EC=50E_{\mathrm J}/E_{\mathrm C}=50 で対角化と比べると、n01=1.088n_{01}=1.088(振動子の近似 1.1181.118)、∣⟨2∣n^∣1⟩∣/n01=1.370\lvert \langle 2 | \hat n | 1 \rangle \rvert /n_{01}=1.370(振動子なら 2=1.414\sqrt2=1.414)である。

近似の裏付け

型 B(極限) | 小さいもの Ω/ωq→0\Omega/\omega_{\mathrm q}\to0 かつ ∣Δ∣/ωq→0\lvert \Delta\rvert /\omega_{\mathrm q}\to0(e±2iωdt\mathrm{e}^{\pm2\mathrm{i}\omega_{\mathrm d}t} の項を時間平均で捨てる) | 証明 条件付き(時間 1/Ω1/\Omega の桁の範囲で誤差 O(Ω/ωq)\mathcal{O}(\Omega/\omega_{\mathrm q})。周期的な係数の常微分方程式に対する平均化の定理による。判定は第2部第4章の札に揃えた) | 行き先 近似の数理 第7章

第3準位への漏れ

transmon は Cooper 対の箱と違い、第3準位 ∣2⟩| 2 \rangle が ∣e⟩| e \rangle から E01+αE_{01}+\alpha の所にあって、α\alpha は E01E_{01} の数 % にすぎない。駆動は ∣e⟩→∣2⟩| e \rangle \to| 2 \rangle の遷移も 2 Ω\sqrt2\,\Omega 程度の強さで起こし、その離調は α\alpha である。∣e⟩| e \rangle と ∣2⟩| 2 \rangle だけの2準位として見積もると、漏れの確率の最大は 2Ω2/(2Ω2+α2)2\Omega^2/(2\Omega^2+\alpha^2) で、Ω\Omega を ∣α∣\lvert \alpha\rvert より十分小さくすれば抑えることができる。そのため transmon の操作のパルスは数十 ns の長さにとり、さらに漏れを打ち消すようにパルスの形を整える(Motzoi らの DRAG)。2準位への射影は、ℏΩ≪∣α∣\hbar\Omega\ll\lvert \alpha\rvert の条件でだけ使える近似である。

例題 2.5.5 transmon の π パルスと漏れ ★★

オリジナル例題 2.5.4 の transmon(α/h=−0.345 GHz\alpha/h=-0.345\,\mathrm{GHz})を Ω/2π=20 MHz\Omega/2\pi=20\,\mathrm{MHz} で駆動する。

  1. (1) 共鳴(Δ=0\Delta=0)で ∣g⟩| g \rangle から ∣e⟩| e \rangle へ完全に移す時間(π パルスの長さ)を求めよ。
  2. (2) 駆動の振動数が 5 MHz5\,\mathrm{MHz} ずれている(Δ/2π=5 MHz\Delta/2\pi=5\,\mathrm{MHz})と、∣e⟩| e \rangle にいる確率の最大はいくらか。
  3. (3) ∣e⟩→∣2⟩| e \rangle \to| 2 \rangle への漏れの確率の最大を、∣e⟩| e \rangle, ∣2⟩| 2 \rangle の2準位と見なして見積もれ(結合 2 Ω\sqrt2\,\Omega、離調 α\alpha)。(1) のパルスが漏れを抑えている理由を述べよ。
解答

(1) (2.5.27) で Δ=0\Delta=0 なら Pe=sin⁡2(Ωt/2)P_e=\sin^2(\Omega t/2) で、Ωt=π\Omega t=\pi のとき1。t=π/Ω=1/(2×20×106) s=25 nst=\pi/\Omega=1/(2\times20\times10^6)\,\mathrm s=25\,\mathrm{ns}。

(2) (Ω/Ω′)2=202/(202+52)=400/425=0.941(\Omega/\Omega')^2=20^2/(20^2+5^2)=400/425=0.941。Δ\Delta は Rabi 振動の振幅を下げ、振動数を Ω′=425×2π MHz\Omega'=\sqrt{425}\times2\pi\,\mathrm{MHz} に上げる。

(3) 結合 2 Ω\sqrt2\,\Omega、離調 α\alpha の2準位の移りの確率の最大は 2Ω2/(2Ω2+α2)=2×202/(2×202+3452)=6.7×10−32\Omega^2/(2\Omega^2+\alpha^2)=2\times20^2/(2\times20^2+345^2)=6.7\times10^{-3}(周波数は MHz 単位)。離調 ∣α∣\lvert \alpha\rvert が結合 2 Ω\sqrt2\,\Omega より12倍以上大きいので、∣2⟩| 2 \rangle へはほとんど移らない。Ω\Omega を大きくして操作を速くするほど漏れは Ω2/α2\Omega^2/\alpha^2 で増えるので、非調和性 ∣α∣≈EC\lvert \alpha\rvert \approx E_{\mathrm C} が操作の速さの上限を決める。

回路の2準位系を LC 共振器と容量でつなぐと、共振器の光子と量子ビットの間の結合 ℏg(a^σ++a^†σ−)\hbar g(\hat a\sigma_++\hat a^\dagger\sigma_-) が得られ、原子と共振器の量子電磁力学と同じ Jaynes–Cummings 模型になる(A. Blais et al. (2004), Phys. Rev. A 69, 062320)。その物理は光学と量子エレクトロニクス 第34章、量子計算の部品としての使い方は相対エントロピーから情報と物理を組み立てる 第6部第1章にある。共振器どうしの結合は章末の総合問題で扱う。

木下の勉強ステータス

木下はまだこの記事を勉強していません。

AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com

← AI生成教科書『虎の巻』の記事一覧へ