量子力学の様々な構成方法 第2部 第5章 §2.5.4 ★★★ 接合のエネルギー − E J cos ( 2 π Φ / Φ 0 ) -E_{\mathrm J}\cos(2\pi\Phi/\Phi_0) − E J cos ( 2 π Φ/ Φ 0 ) を足すと H ^ = 4 E C ( n ^ − n g ) 2 − E J cos φ ^ \hat H=4E_{\mathrm C}(\hat n-n_{\mathrm g})^2-E_{\mathrm J}\cos\hat\varphi H ^ = 4 E C ( n ^ − n g ) 2 − E J cos φ ^ 。電荷が Cooper 対の単位で離散になり、位相が角度になる。E J ≫ E C E_{\mathrm J}\gg E_{\mathrm C} E J ≫ E C で非調和な振動子(transmon)になる
Josephson の2つの関係
2枚の超伝導体を薄い絶縁膜で隔てた接合を考える。Josephson は、超伝導の秩序変数の位相の差 φ \varphi φ について、接合を流れる電流と両端の電圧が
I = I c sin φ , V = ℏ 2 e φ ˙ (2.5.15) I=I_{\mathrm c}\sin\varphi,\qquad V=\frac{\hbar}{2e}\dot\varphi\tag{2.5.15} I = I c sin φ , V = 2 e ℏ φ ˙ ( 2.5.15 )
で結ばれることを予言した(B. D. Josephson (1962), Phys. Lett. 1 , 251)。I c I_{\mathrm c} I c は臨界電流である。第2式の 2 e 2e 2 e は、トンネルして電荷を運ぶ Cooper 対(電子2個)の電荷である。第2式と (2.5.1) の V = Φ ˙ V=\dot\Phi V = Φ ˙ を比べると、位相の差と節点の磁束は
φ = 2 e ℏ Φ = 2 π Φ Φ 0 , Φ 0 = h 2 e = 2.068 × 10 − 15 W b (2.5.16) \varphi=\frac{2e}{\hbar}\Phi=2\pi\frac{\Phi}{\Phi_0},\qquad \Phi_0=\frac{h}{2e}=2.068\times10^{-15}\,\mathrm{Wb}\tag{2.5.16} φ = ℏ 2 e Φ = 2 π Φ 0 Φ , Φ 0 = 2 e h = 2.068 × 1 0 − 15 Wb ( 2.5.16 )
で結ばれる(定数の差は原点の選び方に吸収する)。Φ 0 \Phi_0 Φ 0 は超伝導の磁束量子である。
接合にたまるエネルギーは、電源が接合にした仕事 ∫ I V d t \int IV\,\mathrm{d} t ∫ I V d t で、
I V = I c sin φ ⋅ ℏ 2 e φ ˙ = d d t ( − E J cos φ ) , E J = ℏ I c 2 e = I c Φ 0 2 π (2.5.17) IV=I_{\mathrm c}\sin\varphi\cdot\frac{\hbar}{2e}\dot\varphi=\frac{\mathrm{d} }{\mathrm{d} t}\Bigl(-E_{\mathrm J}\cos\varphi\Bigr),\qquad E_{\mathrm J}=\frac{\hbar I_{\mathrm c}}{2e}=\frac{I_{\mathrm c}\Phi_0}{2\pi}\tag{2.5.17} I V = I c sin φ ⋅ 2 e ℏ φ ˙ = d t d ( − E J cos φ ) , E J = 2 e ℏ I c = 2 π I c Φ 0 ( 2.5.17 )
だから、接合のポテンシャルは U J ( Φ ) = − E J cos ( 2 π Φ / Φ 0 ) U_{\mathrm J}(\Phi)=-E_{\mathrm J}\cos(2\pi\Phi/\Phi_0) U J ( Φ ) = − E J cos ( 2 π Φ/ Φ 0 ) である。φ \varphi φ が小さいときは U J ≈ − E J + Φ 2 / ( 2 L J ) U_{\mathrm J}\approx-E_{\mathrm J}+\Phi^2/(2L_{\mathrm J}) U J ≈ − E J + Φ 2 / ( 2 L J ) 、L J = ℏ / ( 2 e I c ) L_{\mathrm J}=\hbar/(2eI_{\mathrm c}) L J = ℏ/ ( 2 e I c ) で、接合はインダクタンス L J L_{\mathrm J} L J のインダクタとして振る舞う。φ \varphi φ が大きくなると放物線から外れるので、接合は「非線形のインダクタ」である。
接合の Hamiltonian
図 2.5.1 (b) のように、接合と並列に静電容量 C C C を置き、島をゲート容量 C g C_{\mathrm g} C g を通して電圧 V g V_{\mathrm g} V g につなぐ。Lagrangian は
L = C 2 Φ ˙ 2 + C g 2 ( Φ ˙ − V g ) 2 + E J cos 2 π Φ Φ 0 L=\frac C2\dot\Phi^2+\frac{C_{\mathrm g}}{2}(\dot\Phi-V_{\mathrm g})^2+E_{\mathrm J}\cos\frac{2\pi\Phi}{\Phi_0} L = 2 C Φ ˙ 2 + 2 C g ( Φ ˙ − V g ) 2 + E J cos Φ 0 2 π Φ
で、運動量は Q = C Σ Φ ˙ − C g V g Q=C_\Sigma\dot\Phi-C_{\mathrm g}V_{\mathrm g} Q = C Σ Φ ˙ − C g V g (C Σ = C + C g C_\Sigma=C+C_{\mathrm g} C Σ = C + C g )、Legendre 変換はH = ( Q + C g V g ) 2 / ( 2 C Σ ) − E J cos ( 2 π Φ / Φ 0 ) H=(Q+C_{\mathrm g}V_{\mathrm g})^2/(2C_\Sigma)-E_{\mathrm J}\cos(2\pi\Phi/\Phi_0) H = ( Q + C g V g ) 2 / ( 2 C Σ ) − E J cos ( 2 π Φ/ Φ 0 ) (定数 − C g V g 2 / 2 -C_{\mathrm g}V_{\mathrm g}^2/2 − C g V g 2 /2 を落とした)を与える。電荷を Cooper 対の数で、磁束を位相で測り直すと
H ^ = 4 E C ( n ^ − n g ) 2 − E J cos φ ^ , n ^ = Q ^ 2 e , φ ^ = 2 e ℏ Φ ^ , E C = e 2 2 C Σ , n g = − C g V g 2 e (2.5.18) \hat H=4E_{\mathrm C}(\hat n-n_{\mathrm g})^2-E_{\mathrm J}\cos\hat\varphi,\qquad
\hat n=\frac{\hat Q}{2e},\quad \hat\varphi=\frac{2e}{\hbar}\hat\Phi,\quad E_{\mathrm C}=\frac{e^2}{2C_\Sigma},\quad n_{\mathrm g}=-\frac{C_{\mathrm g}V_{\mathrm g}}{2e}\tag{2.5.18} H ^ = 4 E C ( n ^ − n g ) 2 − E J cos φ ^ , n ^ = 2 e Q ^ , φ ^ = ℏ 2 e Φ ^ , E C = 2 C Σ e 2 , n g = − 2 e C g V g ( 2.5.18 )
となる(4 E C n 2 = 4 e 2 n 2 / ( 2 C Σ ) = Q 2 / ( 2 C Σ ) 4E_{\mathrm C}n^2=4e^2n^2/(2C_\Sigma)=Q^2/(2C_\Sigma) 4 E C n 2 = 4 e 2 n 2 / ( 2 C Σ ) = Q 2 / ( 2 C Σ ) )。n g n_{\mathrm g} n g の符号は電荷の向きの取り方で変わるが、後で示すようにスペクトルは n g n_{\mathrm g} n g の偶関数なので答えに関わらない。交換関係 (2.5.6) は [ φ ^ , n ^ ] = ( 2 e / ℏ ) ( 1 / 2 e ) i ℏ = i \left[\hat\varphi,\,\hat n\right]=(2e/\hbar)(1/2e)\mathrm{i}\hbar=\mathrm{i} [ φ ^ , n ^ ] = ( 2 e /ℏ ) ( 1/2 e ) i ℏ = i になる。E C E_{\mathrm C} E C は電子1個を島に載せる静電エネルギーで、n g n_{\mathrm g} n g はゲートが誘起した電荷を Cooper 対の数で測ったもの(連続の値をとる)である。
電荷は離散、位相は角度
島が接合と容量だけで外とつながっているとき、島の電荷が変わるのは接合を Cooper 対がトンネルするときだけである。そのため島の過剰な Cooper 対の数 n n n は整数で、n ^ \hat n n ^ の固有値は n ∈ Z n\in\mathbb{Z} n ∈ Z に限られる(固有状態を ∣ n ⟩ | n \rangle ∣ n ⟩ と書く)。n ^ \hat n n ^ のスペクトルが整数なら、それと共役な φ \varphi φ は 2 π 2\pi 2 π を周期とする角度で、φ \varphi φ と φ + 2 π \varphi+2\pi φ + 2 π は同じ状態を表す。これは平面上の回転子の角運動量 L ^ z \hat L_z L ^ z と方位角の関係と同じである。
このとき φ ^ \hat\varphi φ ^ そのものは自己共役な演算子として定義できず(周期関数の掛け算で閉じない)、使えるのは周期関数 e ± i φ ^ \mathrm{e}^{\pm\mathrm{i}\hat\varphi} e ± i φ ^ だけである。交換関係は
[ n ^ , e i φ ^ ] = e i φ ^ , e i φ ^ ∣ n ⟩ = ∣ n + 1 ⟩ , cos φ ^ = 1 2 ∑ n ∈ Z ( ∣ n + 1 ⟩ ⟨ n ∣ + ∣ n ⟩ ⟨ n + 1 ∣ ) (2.5.19) \left[\hat n,\,\mathrm{e}^{\mathrm{i}\hat\varphi}\right]=\mathrm{e}^{\mathrm{i}\hat\varphi},\qquad \mathrm{e}^{\mathrm{i}\hat\varphi}| n \rangle =| n+1 \rangle ,\qquad
\cos\hat\varphi=\frac12\sum_{n\in\mathbb{Z}}\bigl(| n+1 \rangle \langle n |+| n \rangle \langle n+1 |\bigr)\tag{2.5.19} [ n ^ , e i φ ^ ] = e i φ ^ , e i φ ^ ∣ n ⟩ = ∣ n + 1 ⟩ , cos φ ^ = 2 1 n ∈ Z ∑ ( ∣ n + 1 ⟩ ⟨ n ∣ + ∣ n ⟩ ⟨ n + 1∣ ) ( 2.5.19 )
と置き換わる。e i φ ^ \mathrm{e}^{\mathrm{i}\hat\varphi} e i φ ^ は Cooper 対を1個島へ運び込む演算子で、− E J cos φ ^ -E_{\mathrm J}\cos\hat\varphi − E J cos φ ^ は Cooper 対が1個ずつトンネルする項である。電荷の基底 { ∣ n ⟩ } \{| n \rangle \} { ∣ n ⟩} で書くと (2.5.18) は3重対角行列
⟨ n ′ ∣ H ^ ∣ n ⟩ = 4 E C ( n − n g ) 2 δ n ′ n − E J 2 ( δ n ′ , n + 1 + δ n ′ , n − 1 ) (2.5.20) \langle n' | \hat H | n \rangle =4E_{\mathrm C}(n-n_{\mathrm g})^2\delta_{n'n}-\frac{E_{\mathrm J}}{2}\bigl(\delta_{n',n+1}+\delta_{n',n-1}\bigr)\tag{2.5.20} ⟨ n ′ ∣ H ^ ∣ n ⟩ = 4 E C ( n − n g ) 2 δ n ′ n − 2 E J ( δ n ′ , n + 1 + δ n ′ , n − 1 ) ( 2.5.20 )
になり、∣ n ∣ ≤ 40 \lvert n\rvert \le40 ∣ n ∣ ≤ 40 程度で打ち切って対角化すれば低い準位を十分な精度で得ることができる。この章の数値はすべてこの対角化による。
命題 2.5.1(n g n_{\mathrm g} n g についての周期性と対称性)
(2.5.18) の固有値 E m ( n g ) E_m(n_{\mathrm g}) E m ( n g ) は n g n_{\mathrm g} n g の周期1の偶関数である:E m ( n g + 1 ) = E m ( n g ) E_m(n_{\mathrm g}+1)=E_m(n_{\mathrm g}) E m ( n g + 1 ) = E m ( n g ) 、E m ( − n g ) = E m ( n g ) E_m(-n_{\mathrm g})=E_m(n_{\mathrm g}) E m ( − n g ) = E m ( n g ) 。
証明 S ^ = e i φ ^ \hat S=\mathrm{e}^{\mathrm{i}\hat\varphi} S ^ = e i φ ^ (S ^ ∣ n ⟩ = ∣ n + 1 ⟩ \hat S| n \rangle =| n+1 \rangle S ^ ∣ n ⟩ = ∣ n + 1 ⟩ )はユニタリで、S ^ † n ^ S ^ = n ^ + 1 \hat S^\dagger\hat n\hat S=\hat n+1 S ^ † n ^ S ^ = n ^ + 1 、S ^ † cos φ ^ S ^ = cos φ ^ \hat S^\dagger\cos\hat\varphi\,\hat S=\cos\hat\varphi S ^ † cos φ ^ S ^ = cos φ ^ を満たす。したがって S ^ † H ^ ( n g + 1 ) S ^ = 4 E C ( n ^ + 1 − n g − 1 ) 2 − E J cos φ ^ = H ^ ( n g ) \hat S^\dagger\hat H(n_{\mathrm g}+1)\hat S=4E_{\mathrm C}(\hat n+1-n_{\mathrm g}-1)^2-E_{\mathrm J}\cos\hat\varphi=\hat H(n_{\mathrm g}) S ^ † H ^ ( n g + 1 ) S ^ = 4 E C ( n ^ + 1 − n g − 1 ) 2 − E J cos φ ^ = H ^ ( n g ) で、両者はユニタリ同値だから固有値が等しい。偶関数であることは、P ^ ∣ n ⟩ = ∣ − n ⟩ \hat P| n \rangle =| -n \rangle P ^ ∣ n ⟩ = ∣ − n ⟩ で定めるユニタリ P ^ \hat P P ^ が P ^ n ^ P ^ = − n ^ \hat P\hat n\hat P=-\hat n P ^ n ^ P ^ = − n ^ 、P ^ cos φ ^ P ^ = cos φ ^ \hat P\cos\hat\varphi\hat P=\cos\hat\varphi P ^ cos φ ^ P ^ = cos φ ^ を満たし、P ^ H ^ ( n g ) P ^ = H ^ ( − n g ) \hat P\hat H(n_{\mathrm g})\hat P=\hat H(-n_{\mathrm g}) P ^ H ^ ( n g ) P ^ = H ^ ( − n g ) となることから従う。
2つの領域:Cooper 対の箱と transmon
比 E J / E C E_{\mathrm J}/E_{\mathrm C} E J / E C で、系のふるまいは2つの領域に分かれる(図 2.5.2 )。E J ≲ E C E_{\mathrm J}\lesssim E_{\mathrm C} E J ≲ E C では静電エネルギーが支配的で、準位は n g n_{\mathrm g} n g によって大きく動く。電荷の基底 ∣ n ⟩ | n \rangle ∣ n ⟩ がほぼ固有状態で、トンネルの項はそれらを混ぜる小さな摂動である。これを Cooper 対の箱と呼ぶ(V. Bouchiat et al. (1998), Phys. Scr. T76 , 165)。準位が n g n_{\mathrm g} n g に敏感なので、島の周りの電荷のゆらぎがそのまま準位のゆらぎになる。
E J ≫ E C E_{\mathrm J}\gg E_{\mathrm C} E J ≫ E C では位相のポテンシャル − E J cos φ -E_{\mathrm J}\cos\varphi − E J cos φ が深い井戸を作り、波動関数は φ = 0 \varphi=0 φ = 0 の近くに局在する。準位は n g n_{\mathrm g} n g にほとんど依らなくなる。Koch らはこの領域を使って電荷のゆらぎに強い量子ビットを設計し、transmon と名づけた(J. Koch et al. (2007), Phys. Rev. A 76 , 042319)。
図 2.5.2 (2.5.18) の低い3準位の n g n_{\mathrm g} n g 依存(電荷の基底の対角化による厳密な値なので点線)。縦軸は n g = 1 2 n_{\mathrm g}=\frac12 n g = 2 1 での E 1 − E 0 E_1-E_0 E 1 − E 0 を単位にした。左の Cooper 対の箱(E J / E C = 1 E_{\mathrm J}/E_{\mathrm C}=1 E J / E C = 1 )では準位が n g n_{\mathrm g} n g で大きく動き、n g = ± 1 2 n_{\mathrm g}=\pm\frac12 n g = ± 2 1 で m = 0 m=0 m = 0 , 1 1 1 が避け合う(2準位に射影した近似では間隔 E J E_{\mathrm J} E J )。右の E J / E C = 10 E_{\mathrm J}/E_{\mathrm C}=10 E J / E C = 10 では準位がほぼ平らになり、間隔が m m m とともに狭まる(非調和性)。どちらも周期1・偶関数である(命題 2.5.1 )。
transmon の準位:非調和な振動子
E J ≫ E C E_{\mathrm J}\gg E_{\mathrm C} E J ≫ E C のとき、φ \varphi φ のゆらぎは小さいので cos φ = 1 − φ 2 / 2 + φ 4 / 24 − ⋯ \cos\varphi=1-\varphi^2/2+\varphi^4/24-\cdots cos φ = 1 − φ 2 /2 + φ 4 /24 − ⋯ と展開し、φ \varphi φ を角度でなく実数直線の座標として扱う。2次までで止めると H ^ ≈ − E J + 4 E C n ^ 2 + E J φ ^ 2 / 2 \hat H\approx-E_{\mathrm J}+4E_{\mathrm C}\hat n^2+E_{\mathrm J}\hat\varphi^2/2 H ^ ≈ − E J + 4 E C n ^ 2 + E J φ ^ 2 /2 で、これは「質量」1 / ( 8 E C ) 1/(8E_{\mathrm C}) 1/ ( 8 E C ) (ℏ \hbar ℏ を含めた単位で)、「ばね定数」E J E_{\mathrm J} E J の振動子である。振動の量子は
ℏ ω p = 8 E J E C (2.5.21) \hbar\omega_{\mathrm p}=\sqrt{8E_{\mathrm J}E_{\mathrm C}}\tag{2.5.21} ℏ ω p = 8 E J E C ( 2.5.21 )
で、接合のプラズマ振動数と呼ぶ(ω p = 1 / L J C Σ \omega_{\mathrm p}=1/\sqrt{L_{\mathrm J}C_\Sigma} ω p = 1/ L J C Σ と同じ)。この振動子の生成消滅演算子 b ^ \hat b b ^ , b ^ † \hat b^\dagger b ^ † を使うと
φ ^ = ( 2 E C E J ) 1 / 4 ( b ^ + b ^ † ) , n ^ = i 2 ( E J 2 E C ) 1 / 4 ( b ^ † − b ^ ) (2.5.22) \hat\varphi=\Bigl(\frac{2E_{\mathrm C}}{E_{\mathrm J}}\Bigr)^{1/4}(\hat b+\hat b^\dagger),\qquad
\hat n=\frac{\mathrm{i}}2\Bigl(\frac{E_{\mathrm J}}{2E_{\mathrm C}}\Bigr)^{1/4}(\hat b^\dagger-\hat b)\tag{2.5.22} φ ^ = ( E J 2 E C ) 1/4 ( b ^ + b ^ † ) , n ^ = 2 i ( 2 E C E J ) 1/4 ( b ^ † − b ^ ) ( 2.5.22 )
である(φ z p f n z p f = 1 2 \varphi_{\mathrm{zpf}}n_{\mathrm{zpf}}=\frac12 φ zpf n zpf = 2 1 が [ φ ^ , n ^ ] = i \left[\hat\varphi,\,\hat n\right]=\mathrm{i} [ φ ^ , n ^ ] = i に対応する)。E J / E C E_{\mathrm J}/E_{\mathrm C} E J / E C が大きいほど φ \varphi φ のゆらぎ ( 2 E C / E J ) 1 / 4 (2E_{\mathrm C}/E_{\mathrm J})^{1/4} ( 2 E C / E J ) 1/4 は小さく、これは (2.5.11) の「Z ≪ R K Z\ll R_{\mathrm K} Z ≪ R K なら磁束のゆらぎは Φ 0 \Phi_0 Φ 0 より小さい」と同じことを接合の言葉で言っている。
4次の項 − E J φ 4 / 24 = − ( E C / 12 ) ( b ^ + b ^ † ) 4 -E_{\mathrm J}\varphi^4/24=-(E_{\mathrm C}/12)(\hat b+\hat b^\dagger)^4 − E J φ 4 /24 = − ( E C /12 ) ( b ^ + b ^ † ) 4 を1次の摂動として扱う。⟨ k ∣ ( b ^ + b ^ † ) 4 ∣ k ⟩ = 6 k 2 + 6 k + 3 \langle k | (\hat b+\hat b^\dagger)^4 | k \rangle =6k^2+6k+3 ⟨ k ∣ ( b ^ + b ^ † ) 4 ∣ k ⟩ = 6 k 2 + 6 k + 3 (第2部第4章 の行列要素 k \sqrt k k を4回掛けて数える)から
E k ≈ − E J + 8 E J E C ( k + 1 2 ) − E C 12 ( 6 k 2 + 6 k + 3 ) , (2.5.23) E_k\approx-E_{\mathrm J}+\sqrt{8E_{\mathrm J}E_{\mathrm C}}\Bigl(k+\frac12\Bigr)-\frac{E_{\mathrm C}}{12}(6k^2+6k+3),\tag{2.5.23} E k ≈ − E J + 8 E J E C ( k + 2 1 ) − 12 E C ( 6 k 2 + 6 k + 3 ) , ( 2.5.23 )
E 01 ≡ E 1 − E 0 ≈ 8 E J E C − E C , α ≡ ( E 2 − E 1 ) − ( E 1 − E 0 ) ≈ − E C (2.5.24) E_{01}\equiv E_1-E_0\approx\sqrt{8E_{\mathrm J}E_{\mathrm C}}-E_{\mathrm C},\qquad
\alpha\equiv(E_2-E_1)-(E_1-E_0)\approx-E_{\mathrm C}\tag{2.5.24} E 01 ≡ E 1 − E 0 ≈ 8 E J E C − E C , α ≡ ( E 2 − E 1 ) − ( E 1 − E 0 ) ≈ − E C ( 2.5.24 )
を得る。準位の間隔は上へ行くほど E C E_{\mathrm C} E C ずつ狭まる。この非調和性 α \alpha α が、最低の2準位だけを選んで操作することを可能にする(2.5.5 節)。解析力学から様々な理論分野を組み立てる 第1部第15章 は同じ間隔の狭まりを古典の作用変数 H ( J ) H(J) H ( J ) の曲がりから出した。量子の摂動で出した (2.5.24) はその量子版である。
近似の裏付け
型 B(極限) | 小さいもの E C / E J → 0 E_{\mathrm C}/E_{\mathrm J}\to0 E C / E J → 0 (位相のゆらぎ ( 2 E C / E J ) 1 / 4 (2E_{\mathrm C}/E_{\mathrm J})^{1/4} ( 2 E C / E J ) 1/4 が小さい) | 証明 証明済み(非退化な極小のまわりの調和近似で、低い固有値が E C / E J \sqrt{E_{\mathrm C}/E_{\mathrm J}} E C / E J の漸近展開を持つ。B. Simon (1983), Ann. Inst. H. Poincaré A 38 , 295) | 行き先 近似の数理 第7章
図 2.5.3 は E J / E C = 10 E_{\mathrm J}/E_{\mathrm C}=10 E J / E C = 10 の井戸と準位を、調和近似と並べたものである。上の準位ほど調和近似から下へずれ、間隔が詰まる。
図 2.5.3 E J / E C = 10 E_{\mathrm J}/E_{\mathrm C}=10 E J / E C = 10 、n g = 0 n_{\mathrm g}=0 n g = 0 の位相の井戸。点線(黒)が厳密なポテンシャル − E J cos φ -E_{\mathrm J}\cos\varphi − E J cos φ 、実線(灰)が調和近似 − E J + E J φ 2 / 2 -E_{\mathrm J}+E_{\mathrm J}\varphi^2/2 − E J + E J φ 2 /2 。青の点線が厳密な準位 E k / E C = − 5.80 , 2.10 , 7.45 , 16.65 E_k/E_{\mathrm C}=-5.80,\ 2.10,\ 7.45,\ 16.65 E k / E C = − 5.80 , 2.10 , 7.45 , 16.65 、赤の実線が調和近似の準位 − 5.53 , 3.42 , 12.36 -5.53,\ 3.42,\ 12.36 − 5.53 , 3.42 , 12.36 である。厳密な準位の間隔は 7.90 , 5.35 7.90,\ 5.35 7.90 , 5.35 と狭まる。
電荷分散:周期ポテンシャルのバンドとして
展開で φ \varphi φ を実数直線の座標として扱ったとき、n g n_{\mathrm g} n g は Hamiltonian から消えた(n ^ → n ^ − n g \hat n\to\hat n-n_{\mathrm g} n ^ → n ^ − n g は平行移動で消せる)。n g n_{\mathrm g} n g への依存が戻ってくるのは、φ \varphi φ が角度であること、すなわち隣の井戸 φ + 2 π \varphi+2\pi φ + 2 π が同じ点であることを思い出したときである。φ \varphi φ の表示で (2.5.18) を見ると、これは周期 2 π 2\pi 2 π のポテンシャル − E J cos φ -E_{\mathrm J}\cos\varphi − E J cos φ の中の粒子で、n g n_{\mathrm g} n g は Bloch の擬運動量の役をする(2.5.6 節で境界条件の形にする)。準位の n g n_{\mathrm g} n g 依存はバンドの分散であり、その幅は隣の井戸へのトンネルの振幅で決まる。Koch らは Mathieu 方程式の漸近形から、n g = 1 2 n_{\mathrm g}=\frac12 n g = 2 1 と 0 0 0 での準位の差 ϵ m \epsilon_m ϵ m を
ϵ m ≈ ( − 1 ) m E C 2 4 m + 5 m ! 2 π ( E J 2 E C ) m 2 + 3 4 e − 8 E J / E C (2.5.25) \epsilon_m\approx(-1)^mE_{\mathrm C}\,\frac{2^{4m+5}}{m!}\sqrt{\frac2\pi}\Bigl(\frac{E_{\mathrm J}}{2E_{\mathrm C}}\Bigr)^{\frac m2+\frac34}\mathrm{e}^{-\sqrt{8E_{\mathrm J}/E_{\mathrm C}}}\tag{2.5.25} ϵ m ≈ ( − 1 ) m E C m ! 2 4 m + 5 π 2 ( 2 E C E J ) 2 m + 4 3 e − 8 E J / E C ( 2.5.25 )
と与えた(同論文の式 (2.5))。指数 e − 8 E J / E C \mathrm{e}^{-\sqrt{8E_{\mathrm J}/E_{\mathrm C}}} e − 8 E J / E C のために、E J / E C E_{\mathrm J}/E_{\mathrm C} E J / E C を増やすと電荷分散は指数関数的に小さくなり、非調和性 α ≈ − E C \alpha\approx-E_{\mathrm C} α ≈ − E C は冪でしか減らない。これが transmon の設計の要点である。電荷の基底の対角化と比べると、E J / E C = 50 E_{\mathrm J}/E_{\mathrm C}=50 E J / E C = 50 で (2.5.25) に対する数値の比は m = 0 m=0 m = 0 で 0.96 0.96 0.96 、m = 1 m=1 m = 1 で 0.83 0.83 0.83 で、E J / E C = 20 E_{\mathrm J}/E_{\mathrm C}=20 E J / E C = 20 (0.93 0.93 0.93 、0.73 0.73 0.73 )より1に近い。前の因子は漸近形の主要項なので、比は E J / E C → ∞ E_{\mathrm J}/E_{\mathrm C}\to\infty E J / E C → ∞ で1に近づく。
近似の裏付け
型 B(極限) | 小さいもの トンネルの作用 8 E J / E C → ∞ \sqrt{8E_{\mathrm J}/E_{\mathrm C}}\to\infty 8 E J / E C → ∞ (隣の井戸 φ ± 2 π \varphi\pm2\pi φ ± 2 π が遠い) | 証明 証明済み(周期作用素の孤立したバンドの幅が井戸間の Agmon 距離で指数的に小さいこと:B. Simon (1984)、Helffer–Sjöstrand (1984)。(2.5.25) の前の因子は Mathieu 方程式の漸近形による) | 行き先 近似の数理 第12章
例題 2.5.4 transmon の設計値を読む ★★★
オリジナル E C / h = 0.3 G H z E_{\mathrm C}/h=0.3\,\mathrm{GHz} E C / h = 0.3 GHz 、E J / h = 15 G H z E_{\mathrm J}/h=15\,\mathrm{GHz} E J / h = 15 GHz (E J / E C = 50 E_{\mathrm J}/E_{\mathrm C}=50 E J / E C = 50 )の transmon を考える。
(1) 全静電容量 C Σ C_\Sigma C Σ と臨界電流 I c I_{\mathrm c} I c を求めよ。
(2) (2.5.24) で f 01 = E 01 / h f_{01}=E_{01}/h f 01 = E 01 / h と α / h \alpha/h α / h を見積もれ。電荷の基底の対角化では f 01 = 5.683 G H z f_{01}=5.683\,\mathrm{GHz} f 01 = 5.683 GHz 、α / h = − 0.345 G H z \alpha/h=-0.345\,\mathrm{GHz} α / h = − 0.345 GHz になる。見積もりとの差は何の次数か。
(3) 第1励起準位の電荷分散 ∣ ϵ 1 ∣ / h \lvert \epsilon_1\rvert /h ∣ ϵ 1 ∣ / h は対角化で約 12 k H z 12\,\mathrm{kHz} 12 kHz である。島の電荷が n g n_{\mathrm g} n g の全範囲でゆらいでも、f 01 f_{01} f 01 は相対的にどれだけしか動かないか。
解答 (1) C Σ = e 2 / ( 2 E C ) = ( 1.602 × 10 − 19 ) 2 / ( 2 × 6.626 × 10 − 34 × 0.3 × 10 9 ) = 64.6 f F C_\Sigma=e^2/(2E_{\mathrm C})=(1.602\times10^{-19})^2/(2\times6.626\times10^{-34}\times0.3\times10^9)=64.6\,\mathrm{fF} C Σ = e 2 / ( 2 E C ) = ( 1.602 × 1 0 − 19 ) 2 / ( 2 × 6.626 × 1 0 − 34 × 0.3 × 1 0 9 ) = 64.6 fF 。I c = 2 e E J / ℏ = 2 e ⋅ 2 π × 15 × 10 9 s − 1 = 30.2 n A I_{\mathrm c}=2eE_{\mathrm J}/\hbar=2e\cdot2\pi\times15\times10^9\,\mathrm{s^{-1}}=30.2\,\mathrm{nA} I c = 2 e E J /ℏ = 2 e ⋅ 2 π × 15 × 1 0 9 s − 1 = 30.2 nA 。
(2) 8 E J E C / h = 8 × 15 × 0.3 G H z = 6.0 G H z \sqrt{8E_{\mathrm J}E_{\mathrm C}}/h=\sqrt{8\times15\times0.3}\,\mathrm{GHz}=6.0\,\mathrm{GHz} 8 E J E C / h = 8 × 15 × 0.3 GHz = 6.0 GHz なので f 01 ≈ 6.0 − 0.3 = 5.7 G H z f_{01}\approx6.0-0.3=5.7\,\mathrm{GHz} f 01 ≈ 6.0 − 0.3 = 5.7 GHz 、α / h ≈ − 0.3 G H z \alpha/h\approx-0.3\,\mathrm{GHz} α / h ≈ − 0.3 GHz 。対角化との差は f 01 f_{01} f 01 で 0.017 G H z 0.017\,\mathrm{GHz} 0.017 GHz (0.06 E C 0.06E_{\mathrm C} 0.06 E C )、α \alpha α で 0.045 G H z 0.045\,\mathrm{GHz} 0.045 GHz (0.15 E C 0.15E_{\mathrm C} 0.15 E C )である。展開 (2.5.22) の小さいパラメータは ( E C / E J ) 1 / 2 (E_{\mathrm C}/E_{\mathrm J})^{1/2} ( E C / E J ) 1/2 で、次の補正は E C × O ( ( E C / E J ) 1 / 2 ) E_{\mathrm C}\times\mathcal{O}((E_{\mathrm C}/E_{\mathrm J})^{1/2}) E C × O (( E C / E J ) 1/2 ) の大きさ(E C × 0.14 E_{\mathrm C}\times0.14 E C × 0.14 の桁)なので、差はこの次数にある。
(3) f 01 f_{01} f 01 の変化は ∣ ϵ 1 − ϵ 0 ∣ / h ≈ ∣ ϵ 1 ∣ / h ≈ 12 k H z \lvert \epsilon_1-\epsilon_0\rvert /h\approx\lvert \epsilon_1\rvert /h\approx12\,\mathrm{kHz} ∣ ϵ 1 − ϵ 0 ∣ / h ≈ ∣ ϵ 1 ∣ / h ≈ 12 kHz 程度で、5.7 G H z 5.7\,\mathrm{GHz} 5.7 GHz の 2 × 10 − 6 2\times10^{-6} 2 × 1 0 − 6 にすぎない。E J / E C = 1 E_{\mathrm J}/E_{\mathrm C}=1 E J / E C = 1 の Cooper 対の箱では、準位の動きは遷移エネルギーそのものと同じ桁だった(図 2.5.2 左)。
量子化された位相の実験
接合の位相 φ \varphi φ は、超伝導体の中の膨大な数の Cooper 対が共有する1つの集団の座標である。これに量子力学がそのまま当てはまるかは、1980年代に実験で確かめられた。Martinis・Devoret・Clarke は電流を流した接合の位相の井戸にマイクロ波を当て、準位が離散的であることを測った(J. M. Martinis, M. H. Devoret and J. Clarke (1985), Phys. Rev. Lett. 55 , 1543)。測った準位の位置は、古典の領域で測った接合の値を入れた量子の計算と一致した。この一連の研究に 2025 年のノーベル物理学賞が与えられた(電気回路での巨視的な量子トンネルとエネルギーの量子化の発見)。この章の言葉で言えば、正準量子化の処方を回路の磁束と電荷に当ててよいことの実験的な裏付けである。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com