量子力学の様々な構成方法 第2部 第5章 §2.5.5 ★★★ Cooper 対の箱を n = 0 , 1 n=0,1 n = 0 , 1 へ、transmon を最低の2準位へ射影すると共通例題②の Hamiltonian になり、Rabi 振動 ( Ω / Ω ′ ) 2 sin 2 ( Ω ′ t / 2 ) (\Omega/\Omega')^2\sin^2(\Omega't/2) ( Ω/ Ω ′ ) 2 sin 2 ( Ω ′ t /2 ) が回路の言葉で出る
Cooper 対の箱:静的な2準位系
E J ≪ E C E_{\mathrm J}\ll E_{\mathrm C} E J ≪ E C で n g n_{\mathrm g} n g が 1 2 \frac12 2 1 に近いとき、電荷の状態 ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ と ∣ 1 ⟩ | 1 \rangle ∣1 ⟩ のエネルギー 4 E C n g 2 4E_{\mathrm C}n_{\mathrm g}^2 4 E C n g 2 と 4 E C ( 1 − n g ) 2 4E_{\mathrm C}(1-n_{\mathrm g})^2 4 E C ( 1 − n g ) 2 は近く、ほかの電荷の状態は ∼ 4 E C \sim4E_{\mathrm C} ∼ 4 E C 以上離れている。(2.5.20) を { ∣ 1 ⟩ , ∣ 0 ⟩ } \{| 1 \rangle ,| 0 \rangle \} { ∣1 ⟩ , ∣0 ⟩} の張る部分空間へ射影し、∣ e ⟩ = ∣ 1 ⟩ | e \rangle =| 1 \rangle ∣ e ⟩ = ∣1 ⟩ 、∣ g ⟩ = ∣ 0 ⟩ | g \rangle =| 0 \rangle ∣ g ⟩ = ∣0 ⟩ (σ z ∣ e ⟩ = + ∣ e ⟩ \sigma_z| e \rangle =+| e \rangle σ z ∣ e ⟩ = + ∣ e ⟩ )と置くと、定数を除いて
H ^ b o x = 4 E C ( 1 − 2 n g ) 2 σ z − E J 2 σ x = − ℏ Δ 2 σ z − ℏ Ω 2 σ x , ℏ Ω = E J , ℏ Δ = − 4 E C ( 1 − 2 n g ) (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} H ^ box = 2 4 E C ( 1 − 2 n g ) σ z − 2 E J σ x = − 2 ℏΔ σ z − 2 ℏΩ σ x , ℏΩ = E J , ℏΔ = − 4 E C ( 1 − 2 n g ) ( 2.5.26 )
となる。2つの電荷の状態の差 4 E C ( 1 − 2 n g ) 4E_{\mathrm C}(1-2n_{\mathrm g}) 4 E C ( 1 − 2 n g ) は 4 E C ( 1 − n g ) 2 − 4 E C n g 2 4E_{\mathrm C}(1-n_{\mathrm g})^2-4E_{\mathrm C}n_{\mathrm g}^2 4 E C ( 1 − n g ) 2 − 4 E C n g 2 から出る。
これは共通例題②の Hamiltonian − ( ℏ Δ / 2 ) σ z − ( ℏ Ω / 2 ) σ x -(\hbar\Delta/2)\sigma_z-(\hbar\Omega/2)\sigma_x − ( ℏΔ/2 ) σ z − ( ℏΩ/2 ) σ x と同じ形で、ここでは Ω \Omega Ω が外から当てた駆動でなく接合のトンネル、Δ \Delta Δ がゲート電圧で決まる電荷のエネルギーの差である。n g n_{\mathrm g} n g を 1 2 \frac12 2 1 から遠い値(基底状態がほぼ ∣ 0 ⟩ | 0 \rangle ∣0 ⟩ )から急に n g n_{\mathrm g} n g へ切り替えて時間 t t t だけ置くと、第1部第2章 ・第2部第3章 で出した式のとおり
P ∣ 1 ⟩ ( t ) = ( Ω Ω ′ ) 2 sin 2 Ω ′ t 2 , Ω ′ = Ω 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} P ∣1 ⟩ ( t ) = ( Ω ′ Ω ) 2 sin 2 2 Ω ′ t , Ω ′ = Ω 2 + Δ 2 ( 2.5.27 )
で島に Cooper 対が1個入っている確率が振動する。縮退点 n g = 1 2 n_{\mathrm g}=\frac12 n g = 2 1 では Δ = 0 \Delta=0 Δ = 0 で、周期 h / E J h/E_{\mathrm J} h / E 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 \} { ∣0 ⟩ , ∣1 ⟩} の張る部分空間(ほかの電荷の状態までの距離 ∼ 4 E C \sim4E_{\mathrm C} ∼ 4 E C が E J E_{\mathrm J} E J より十分大きい) | 証明 条件付き(射影した2準位の固有値の誤差は2次の摂動で O ( E J 2 / E C ) \mathcal{O}(E_{\mathrm J}^2/E_{\mathrm C}) O ( E J 2 / E C ) 。図 2.5.2 左の E J / E C = 1 E_{\mathrm J}/E_{\mathrm C}=1 E J / E C = 1 では第3準位が離れているが、誤差の厳密な上界は、Josephson の項の大きさ(ノルム ≤ E J \le E_{\mathrm J} ≤ E 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 ∣ g ⟩ = ∣ k = 0 ⟩ 、∣ e ⟩ = ∣ k = 1 ⟩ | e \rangle =| k=1 \rangle ∣ e ⟩ = ∣ k = 1 ⟩ の間隔は ℏ ω q = E 01 \hbar\omega_{\mathrm q}=E_{01} ℏ ω q = E 01 である。ゲートに V g = V d c + V 0 cos ω d t V_{\mathrm g}=V_{\mathrm{dc}}+V_0\cos\omega_{\mathrm d}t V g = V dc + V 0 cos ω d t を加えると、(2.5.18) の ( n ^ − n g ) 2 (\hat n-n_{\mathrm g})^2 ( n ^ − n g ) 2 の交差項から
H ^ d ( t ) = 2 e C g C Σ V 0 cos ( ω d t ) 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} H ^ d ( t ) = 2 e C Σ C g V 0 cos ( ω d t ) n ^ ( 2.5.28 )
が出る(8 E C / ( 2 e ) = 2 e / C Σ 8E_{\mathrm C}/(2e)=2e/C_\Sigma 8 E C / ( 2 e ) = 2 e / C Σ )。n ^ \hat n n ^ の行列要素は (2.5.22) により ⟨ e ∣ n ^ ∣ g ⟩ ≈ ( i / 2 ) ( E J / 8 E C ) 1 / 4 \langle e | \hat n | g \rangle \approx(\mathrm{i}/\sqrt2)(E_{\mathrm J}/8E_{\mathrm C})^{1/4} ⟨ e ∣ n ^ ∣ g ⟩ ≈ ( i / 2 ) ( E J /8 E C ) 1/4 で、純虚数である。∣ e ⟩ | e \rangle ∣ e ⟩ の位相を選び直して行列要素を実数 n 01 = ∣ ⟨ e ∣ n ^ ∣ g ⟩ ∣ n_{01}=\lvert \langle e | \hat n | g \rangle \rvert n 01 = ∣⟨ e ∣ n ^ ∣ g ⟩∣ にし、2準位へ射影するとH ^ = ( ℏ ω q / 2 ) σ z + ℏ Ω cos ( ω d t ) σ x \hat H=(\hbar\omega_{\mathrm q}/2)\sigma_z+\hbar\Omega\cos(\omega_{\mathrm d}t)\sigma_x H ^ = ( ℏ ω q /2 ) σ z + ℏΩ cos ( ω d t ) σ x 、ℏ Ω = 2 e ( C g / C Σ ) V 0 n 01 \hbar\Omega=2e(C_{\mathrm g}/C_\Sigma)V_0n_{01} ℏΩ = 2 e ( C g / C Σ ) V 0 n 01 である(n ^ \hat n n ^ の対角要素は σ z \sigma_z σ z の係数を少しずらすだけなので落とした)。e i ω d t σ z / 2 \mathrm{e}^{\mathrm{i}\omega_{\mathrm d}t\sigma_z/2} e i ω d t σ z /2 で回る座標へ移り、e ± 2 i ω d t \mathrm{e}^{\pm2\mathrm{i}\omega_{\mathrm d}t} e ± 2 i ω d t で振動する項を捨てると(回転波近似)
H ^ r o t = − ℏ Δ 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} H ^ rot = − 2 ℏΔ σ z + 2 ℏΩ σ x , Δ = ω d − ω q ( 2.5.29 )
を得る。σ x \sigma_x σ x の係数の符号は駆動の位相 π \pi π のずれで反転でき、(2.5.27) は Ω 2 \Omega^2 Ω 2 にしか依らないので、共通例題②の正答がそのまま当てはまる。この節の行列要素を E J / E C = 50 E_{\mathrm J}/E_{\mathrm C}=50 E J / E C = 50 で対角化と比べると、n 01 = 1.088 n_{01}=1.088 n 01 = 1.088 (振動子の近似 1.118 1.118 1.118 )、∣ ⟨ 2 ∣ n ^ ∣ 1 ⟩ ∣ / n 01 = 1.370 \lvert \langle 2 | \hat n | 1 \rangle \rvert /n_{01}=1.370 ∣⟨ 2∣ n ^ ∣1 ⟩∣ / n 01 = 1.370 (振動子なら 2 = 1.414 \sqrt2=1.414 2 = 1.414 )である。
近似の裏付け
型 B(極限) | 小さいもの Ω / ω q → 0 \Omega/\omega_{\mathrm q}\to0 Ω/ ω q → 0 かつ ∣ Δ ∣ / ω q → 0 \lvert \Delta\rvert /\omega_{\mathrm q}\to0 ∣ Δ ∣ / ω q → 0 (e ± 2 i ω d t \mathrm{e}^{\pm2\mathrm{i}\omega_{\mathrm d}t} e ± 2 i ω d t の項を時間平均で捨てる) | 証明 条件付き(時間 1 / Ω 1/\Omega 1/Ω の桁の範囲で誤差 O ( Ω / ω q ) \mathcal{O}(\Omega/\omega_{\mathrm q}) O ( Ω/ ω q ) 。周期的な係数の常微分方程式に対する平均化の定理による。判定は第2部第4章の札に揃えた) | 行き先 近似の数理 第7章
第3準位への漏れ
transmon は Cooper 対の箱と違い、第3準位 ∣ 2 ⟩ | 2 \rangle ∣2 ⟩ が ∣ e ⟩ | e \rangle ∣ e ⟩ から E 01 + α E_{01}+\alpha E 01 + α の所にあって、α \alpha α は E 01 E_{01} E 01 の数 % にすぎない。駆動は ∣ e ⟩ → ∣ 2 ⟩ | e \rangle \to| 2 \rangle ∣ e ⟩ → ∣2 ⟩ の遷移も 2 Ω \sqrt2\,\Omega 2 Ω 程度の強さで起こし、その離調は α \alpha α である。∣ e ⟩ | e \rangle ∣ e ⟩ と ∣ 2 ⟩ | 2 \rangle ∣2 ⟩ だけの2準位として見積もると、漏れの確率の最大は 2 Ω 2 / ( 2 Ω 2 + α 2 ) 2\Omega^2/(2\Omega^2+\alpha^2) 2 Ω 2 / ( 2 Ω 2 + α 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 G H z \alpha/h=-0.345\,\mathrm{GHz} α / h = − 0.345 GHz )を Ω / 2 π = 20 M H z \Omega/2\pi=20\,\mathrm{MHz} Ω/2 π = 20 MHz で駆動する。
(1) 共鳴(Δ = 0 \Delta=0 Δ = 0 )で ∣ g ⟩ | g \rangle ∣ g ⟩ から ∣ e ⟩ | e \rangle ∣ e ⟩ へ完全に移す時間(π パルスの長さ)を求めよ。
(2) 駆動の振動数が 5 M H z 5\,\mathrm{MHz} 5 MHz ずれている(Δ / 2 π = 5 M H z \Delta/2\pi=5\,\mathrm{MHz} Δ/2 π = 5 MHz )と、∣ e ⟩ | e \rangle ∣ e ⟩ にいる確率の最大はいくらか。
(3) ∣ e ⟩ → ∣ 2 ⟩ | e \rangle \to| 2 \rangle ∣ e ⟩ → ∣2 ⟩ への漏れの確率の最大を、∣ e ⟩ | e \rangle ∣ e ⟩ , ∣ 2 ⟩ | 2 \rangle ∣2 ⟩ の2準位と見なして見積もれ(結合 2 Ω \sqrt2\,\Omega 2 Ω 、離調 α \alpha α )。(1) のパルスが漏れを抑えている理由を述べよ。
解答 (1) (2.5.27) で Δ = 0 \Delta=0 Δ = 0 なら P e = sin 2 ( Ω t / 2 ) P_e=\sin^2(\Omega t/2) P e = sin 2 ( Ω t /2 ) で、Ω t = π \Omega t=\pi Ω t = π のとき1。t = π / Ω = 1 / ( 2 × 20 × 10 6 ) s = 25 n s t=\pi/\Omega=1/(2\times20\times10^6)\,\mathrm s=25\,\mathrm{ns} t = π /Ω = 1/ ( 2 × 20 × 1 0 6 ) s = 25 ns 。
(2) ( Ω / Ω ′ ) 2 = 20 2 / ( 20 2 + 5 2 ) = 400 / 425 = 0.941 (\Omega/\Omega')^2=20^2/(20^2+5^2)=400/425=0.941 ( Ω/ Ω ′ ) 2 = 2 0 2 / ( 2 0 2 + 5 2 ) = 400/425 = 0.941 。Δ \Delta Δ は Rabi 振動の振幅を下げ、振動数を Ω ′ = 425 × 2 π M H z \Omega'=\sqrt{425}\times2\pi\,\mathrm{MHz} Ω ′ = 425 × 2 π MHz に上げる。
(3) 結合 2 Ω \sqrt2\,\Omega 2 Ω 、離調 α \alpha α の2準位の移りの確率の最大は 2 Ω 2 / ( 2 Ω 2 + α 2 ) = 2 × 20 2 / ( 2 × 20 2 + 345 2 ) = 6.7 × 10 − 3 2\Omega^2/(2\Omega^2+\alpha^2)=2\times20^2/(2\times20^2+345^2)=6.7\times10^{-3} 2 Ω 2 / ( 2 Ω 2 + α 2 ) = 2 × 2 0 2 / ( 2 × 2 0 2 + 34 5 2 ) = 6.7 × 1 0 − 3 (周波数は MHz 単位)。離調 ∣ α ∣ \lvert \alpha\rvert ∣ α ∣ が結合 2 Ω \sqrt2\,\Omega 2 Ω より12倍以上大きいので、∣ 2 ⟩ | 2 \rangle ∣2 ⟩ へはほとんど移らない。Ω \Omega Ω を大きくして操作を速くするほど漏れは Ω 2 / α 2 \Omega^2/\alpha^2 Ω 2 / α 2 で増えるので、非調和性 ∣ α ∣ ≈ E C \lvert \alpha\rvert \approx E_{\mathrm C} ∣ α ∣ ≈ E C が操作の速さの上限を決める。
回路の2準位系を LC 共振器と容量でつなぐと、共振器の光子と量子ビットの間の結合 ℏ g ( a ^ σ + + a ^ † σ − ) \hbar g(\hat a\sigma_++\hat a^\dagger\sigma_-) ℏ g ( a ^ σ + + a ^ † σ − ) が得られ、原子と共振器の量子電磁力学と同じ Jaynes–Cummings 模型になる(A. Blais et al. (2004), Phys. Rev. A 69 , 062320)。その物理は光学と量子エレクトロニクス 第34章 、量子計算の部品としての使い方は相対エントロピーから情報と物理を組み立てる 第6部第1章 にある。共振器どうしの結合は章末の総合問題で扱う。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com