光学と量子エレクトロニクス 第28章 §28.2 ★★★ 縞の傾きが最大の点で励起の確率を測ると、周波数のずれが確率のずれに変わる。原子 N N N 個の測定の二項分布の揺らぎ(量子射影雑音)から、安定度 σ y ( τ ) = 1 2 π ν 0 t f t c / ( N τ ) \sigma_y(\tau)=\frac{1}{2\pi\nu_0t_{\mathrm f}}\sqrt{t_{\mathrm c}/(N\tau)} σ y ( τ ) = 2 π ν 0 t f 1 t c / ( N τ ) が出る
縞の傾きでずれを読む
第24章の Ramsey の縞を、パルスが自由な時間に比べて十分短い極限で書くと、2つ目のパルスの後に原子が励起状態にいる確率は
P e ( Δ ) = cos 2 Δ t f 2 = 1 2 ( 1 + cos Δ t f ) (28.5) P_e(\Delta)=\cos^2\frac{\Delta t_{\mathrm f}}{2}=\frac12\bigl(1+\cos\Delta t_{\mathrm f}\bigr)\tag{28.5} P e ( Δ ) = cos 2 2 Δ t f = 2 1 ( 1 + cos Δ t f ) ( 28.5 )
である。Bloch 球の上では、1つ目の π / 2 \pi/2 π /2 パルスで赤道に倒れた Bloch ベクトルが、自由な時間の間に回転座標で角 Δ t f \Delta t_{\mathrm f} Δ t f だけ回り、その角が2つ目のパルスの後の高さに変わる。したがって Ramsey 分光は、原子の重ね合わせの位相と局部発振器の位相の差を、自由な時間のあいだ積算して読む位相の測定である。
時計を動かすときは、縞の頂上(Δ = 0 \Delta=0 Δ = 0 )で測らない。頂上では d P e d Δ = 0 \frac{\mathrm{d} P_e}{\mathrm{d} \Delta}=0 d Δ d P e = 0 で、小さなずれが確率を1次で変えず、ずれの符号もわからないからである。代わりに、局部発振器の周波数を中心から ± Δ m \pm\Delta_{\mathrm m} ± Δ m だけずらして交互に測る。式 (28.5) を微分すると
d P e d Δ = − t f 2 sin Δ t f (28.6) \frac{\mathrm{d} P_e}{\mathrm{d} \Delta}=-\frac{t_{\mathrm f}}{2}\sin\Delta t_{\mathrm f}\tag{28.6} d Δ d P e = − 2 t f sin Δ t f ( 28.6 )
で、傾きの大きさは Δ m t f = π / 2 \Delta_{\mathrm m}t_{\mathrm f}=\pi/2 Δ m t f = π /2 で最大値 t f / 2 t_{\mathrm f}/2 t f /2 をとる。この点では P e = 1 / 2 P_e=1/2 P e = 1/2 で、縞の高さの半分にあたる(図 28.1 )。両側で測った P e P_e P e の差は、中心が δ \delta δ だけずれていれば 2 × ( t f / 2 ) δ = t f δ 2\times(t_{\mathrm f}/2)\delta=t_{\mathrm f}\delta 2 × ( t f /2 ) δ = t f δ になり、符号でずれの向きがわかる。この差を 0 0 0 に保つように局部発振器を制御する。
図 28.1 Ramsey の縞 P e = cos 2 ( Δ t f / 2 ) P_e=\cos^2(\Delta t_{\mathrm f}/2) P e = cos 2 ( Δ t f /2 ) (青の実線、短いパルスの極限)と、局部発振器を交互にずらして測る2つの点 Δ m t f = ± π / 2 \Delta_{\mathrm m}t_{\mathrm f}=\pm\pi/2 Δ m t f = ± π /2 (赤の点)。赤の線はその点での接線で、傾きの大きさ t f / 2 t_{\mathrm f}/2 t f /2 が最大になる
量子射影雑音
1回の測定では、N N N 個の原子を Ramsey の手順に通し、励起状態にいる原子の数 N e N_e N e を数える。2つ目のパルスの後、各原子は 1 − P e ∣ g ⟩ + P e e i ϕ ∣ e ⟩ \sqrt{1-P_e}| g \rangle +\sqrt{P_e}\mathrm{e}^{\mathrm{i}\phi}| e \rangle 1 − P e ∣ g ⟩ + P e e i ϕ ∣ e ⟩ の重ね合わせにあり、状態を測ると確率 P e P_e P e で ∣ e ⟩ | e \rangle ∣ e ⟩ に、1 − P e 1-P_e 1 − P e で ∣ g ⟩ | g \rangle ∣ g ⟩ に落ちる。原子が互いに独立なら、N e N_e N e は二項分布に従い、推定した確率 P ^ e = N e / N \hat P_e=N_e/N P ^ e = N e / N の標準偏差は
σ P = P e ( 1 − P e ) N (28.7) \sigma_{P}=\sqrt{\frac{P_e(1-P_e)}{N}}\tag{28.7} σ P = N P e ( 1 − P e ) ( 28.7 )
である。原子の数え方の技術を完全にしても残るこの揺らぎを量子射影雑音 と呼ぶ(Itano たち 1993、Phys. Rev. A 47, 3554)。測定が状態を ∣ g ⟩ | g \rangle ∣ g ⟩ か ∣ e ⟩ | e \rangle ∣ e ⟩ に射影することによる揺らぎだからである。
動作点 P e = 1 / 2 P_e=1/2 P e = 1/2 では σ P = 1 / ( 2 N ) \sigma_P=1/(2\sqrt N) σ P = 1/ ( 2 N ) で、これを傾き t f / 2 t_{\mathrm f}/2 t f /2 で割ると、1回の測定で決まる離調の不確かさは
δ Δ = σ P t f / 2 = 1 t f N (28.8) \delta\Delta=\frac{\sigma_P}{t_{\mathrm f}/2}=\frac{1}{t_{\mathrm f}\sqrt N}\tag{28.8} δ Δ = t f /2 σ P = t f N 1 ( 28.8 )
になる。原子1個あたりの位相の不確かさは 1 / N 1/\sqrt N 1/ N で、これは独立な N N N 個の量子系で位相を測るときの標準量子限界である(相対エントロピーから情報と物理を組み立てる 第3部第8章 )。動作点で P e ( 1 − P e ) P_e(1-P_e) P e ( 1 − P e ) が最大になっていても、傾きも最大なので、比 σ P / ∣ d P e / d Δ ∣ \sigma_P/\lvert \mathrm{d} P_e/\mathrm{d}\Delta\rvert σ P / ∣ d P e / d Δ ∣ は動作点で最小になる。
安定度と Allan 分散
1回の測定の周期を t c t_{\mathrm c} t c (原子の準備・冷却・Ramsey の手順・検出を合わせた時間で、t f t_{\mathrm f} t f より長い)とする。時間 τ \tau τ のあいだには τ / t c \tau/t_{\mathrm c} τ / t c 回の独立な測定ができ、平均した不確かさは t c / τ \sqrt{t_{\mathrm c}/\tau} t c / τ 倍に下がる。角周波数を周波数に直し(2 π 2\pi 2 π で割り)、ν 0 \nu_0 ν 0 で割ると、相対的な周波数の不確かさは
σ y ( τ ) = 1 2 π ν 0 t f t c N τ (28.9) \sigma_y(\tau)=\frac{1}{2\pi\nu_0t_{\mathrm f}}\sqrt{\frac{t_{\mathrm c}}{N\tau}}\tag{28.9} σ y ( τ ) = 2 π ν 0 t f 1 N τ t c ( 28.9 )
になる。ν 0 t f \nu_0t_{\mathrm f} ν 0 t f は式 (28.2) の Q 値の半分なので、σ y ( τ ) = 1 π Q t c / ( N τ ) \sigma_y(\tau)=\frac{1}{\pi Q}\sqrt{t_{\mathrm c}/(N\tau)} σ y ( τ ) = π Q 1 t c / ( N τ ) とも書ける。Q 値と原子の数の平方根が大きいほど安定な時計になる。
時計の揺らぎを表すのに、ふつうの標準偏差は使えない。発振器の周波数は時間とともにゆっくり漂い、標準偏差は測る時間を延ばすほど大きくなって定まらないことがある。代わりに、長さ τ \tau τ の区間ごとに平均した y ˉ k \bar y_k y ˉ k を作り、隣り合う区間の差の2乗の平均の半分
σ y 2 ( τ ) = 1 2 ⟨ ( y ˉ k + 1 − y ˉ k ) 2 ⟩ (28.10) \sigma_y^2(\tau)=\frac12\langle (\bar y_{k+1}-\bar y_k)^2 \rangle\tag{28.10} σ y 2 ( τ ) = 2 1 ⟨( y ˉ k + 1 − y ˉ k ) 2 ⟩ ( 28.10 )
を使う。これを Allan 分散 と呼ぶ(D. W. Allan 1966、Proc. IEEE 54, 221)。隣どうしの差を取るので、y y y の一定のずれは消え、ゆっくりした漂いも小さく抑えられる。y y y が各時刻で独立に揺らぐ雑音(白色周波数雑音)なら、Allan 分散は区間の平均の分散そのものになり、σ y ( τ ) ∝ τ − 1 / 2 \sigma_y(\tau)\propto\tau^{-1/2} σ y ( τ ) ∝ τ − 1/2 になる(28.6 節の例題)。量子射影雑音は各回の測定で独立なので、この白色周波数雑音として現れ、式 (28.9) はその Allan 標準偏差である。
図 28.2 量子射影雑音だけで決まる安定度(式 (28.9) )。青は Cs の噴水時計(t f = 0.5 s t_{\mathrm f}=0.5\,\mathrm s t f = 0.5 s 、t c = 1.5 s t_{\mathrm c}=1.5\,\mathrm s t c = 1.5 s 、N = 10 6 N=10^6 N = 1 0 6 )、赤は Sr の光格子時計(t f = 0.5 s t_{\mathrm f}=0.5\,\mathrm s t f = 0.5 s 、t c = 1 s t_{\mathrm c}=1\,\mathrm s t c = 1 s 、N = 10 3 N=10^3 N = 1 0 3 )。どちらも τ − 1 / 2 \tau^{-1/2} τ − 1/2 で下がり、原子の数が 10 3 10^3 1 0 3 分の1でも光の時計は3桁以上安定である
例題 28.3 Cs 噴水と Sr 光格子の安定度 ★★★
オリジナル 量子射影雑音だけで安定度が決まるとして、2つの時計を比べる。
(1) Cs の噴水時計:N = 10 6 N=10^6 N = 1 0 6 、t f = 0.50 s t_{\mathrm f}=0.50\,\mathrm s t f = 0.50 s 、t c = 1.5 s t_{\mathrm c}=1.5\,\mathrm s t c = 1.5 s 。σ y ( 1 s ) \sigma_y(1\,\mathrm s) σ y ( 1 s ) と、σ y = 1 × 10 − 16 \sigma_y=1\times10^{-16} σ y = 1 × 1 0 − 16 に達するまでの平均の時間を求めよ。
(2) Sr の光格子時計:N = 1000 N=1000 N = 1000 、t f = 0.50 s t_{\mathrm f}=0.50\,\mathrm s t f = 0.50 s 、t c = 1.0 s t_{\mathrm c}=1.0\,\mathrm s t c = 1.0 s 。σ y ( 1 s ) \sigma_y(1\,\mathrm s) σ y ( 1 s ) と、σ y = 1 × 10 − 18 \sigma_y=1\times10^{-18} σ y = 1 × 1 0 − 18 に達するまでの時間を求めよ。
(3) 原子の数が 10 3 10^3 1 0 3 分の1なのに Sr のほうが安定な理由を、式 (28.9) の因子に分けて述べよ。
解答 (1) 式 (28.9) に入れると
σ y ( 1 s ) = 1 2 π × 9.193 × 10 9 × 0.50 1.5 10 6 = 3.46 × 10 − 11 × 1.22 × 10 − 3 = 4.24 × 10 − 14 \sigma_y(1\,\mathrm s)=\frac{1}{2\pi\times9.193\times10^{9}\times0.50}\sqrt{\frac{1.5}{10^6}}=3.46\times10^{-11}\times1.22\times10^{-3}=4.24\times10^{-14} σ y ( 1 s ) = 2 π × 9.193 × 1 0 9 × 0.50 1 1 0 6 1.5 = 3.46 × 1 0 − 11 × 1.22 × 1 0 − 3 = 4.24 × 1 0 − 14
である。σ y ∝ τ − 1 / 2 \sigma_y\propto\tau^{-1/2} σ y ∝ τ − 1/2 なので、10 − 16 10^{-16} 1 0 − 16 に達する時間は τ = ( 4.24 × 10 − 14 / 10 − 16 ) 2 s = 1.80 × 10 5 s \tau=(4.24\times10^{-14}/10^{-16})^2\,\mathrm s=1.80\times10^{5}\,\mathrm s τ = ( 4.24 × 1 0 − 14 /1 0 − 16 ) 2 s = 1.80 × 1 0 5 s 、およそ 2.1 2.1 2.1 日である。
(2) 同じように
σ y ( 1 s ) = 1 2 π × 4.292 × 10 14 × 0.50 1.0 1000 = 7.42 × 10 − 16 × 3.16 × 10 − 2 = 2.35 × 10 − 17 \sigma_y(1\,\mathrm s)=\frac{1}{2\pi\times4.292\times10^{14}\times0.50}\sqrt{\frac{1.0}{1000}}=7.42\times10^{-16}\times3.16\times10^{-2}=2.35\times10^{-17} σ y ( 1 s ) = 2 π × 4.292 × 1 0 14 × 0.50 1 1000 1.0 = 7.42 × 1 0 − 16 × 3.16 × 1 0 − 2 = 2.35 × 1 0 − 17
で、10 − 18 10^{-18} 1 0 − 18 に達する時間は ( 2.35 × 10 − 17 / 10 − 18 ) 2 s = 550 s (2.35\times10^{-17}/10^{-18})^2\,\mathrm s=550\,\mathrm s ( 2.35 × 1 0 − 17 /1 0 − 18 ) 2 s = 550 s 、10分足らずである。
(3) 式 (28.9) の因子 1 / ( ν 0 t f ) 1/(\nu_0t_{\mathrm f}) 1/ ( ν 0 t f ) は、t f t_{\mathrm f} t f が同じなので周波数の比 4.67 × 10 4 4.67\times10^4 4.67 × 1 0 4 だけ Sr が小さい。因子 t c / N \sqrt{t_{\mathrm c}/N} t c / N は 1.0 / 1000 / 1.5 / 10 6 = 25.8 \sqrt{1.0/1000}/\sqrt{1.5/10^6}=25.8 1.0/1000 / 1.5/1 0 6 = 25.8 だけ Sr が大きい。差し引き 4.67 × 10 4 / 25.8 = 1.81 × 10 3 4.67\times10^4/25.8=1.81\times10^3 4.67 × 1 0 4 /25.8 = 1.81 × 1 0 3 倍だけ Sr が安定になる。原子の数は平方根でしか入らないのに対し、周波数は1乗で入る。光の時計の利点は、この周波数の比にある。
量子射影雑音の限界に実際に届くには、局部発振器自身の揺らぎが十分小さい必要がある。測定の周期のうち原子を準備する時間は原子が発振器を見ていないので、その間の発振器の揺らぎが測定の結果に折り返して入る(Dick 効果と呼ばれる)。光の時計では、時計のレーザーを超低膨張ガラスや単結晶シリコンで作った光共振器(第5章 )に合わせて、秒の程度の時間でまず安定にしておく。
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com