光学と量子エレクトロニクス 第29章 §29.3 ★★★ 一様な重力の中で、Mach–Zehnder 型の位相差は Δ ϕ = ( k e f f ⋅ g − α c h ) T 2 \Delta\phi=(\bm{k}_{\mathrm{eff}}\cdot\bm{g}-\alpha_{\mathrm{ch}})T^2 Δ ϕ = ( k eff ⋅ g − α ch ) T 2 。経路積分の位相に分けると、伝搬の位相は2本で打ち消し合い、光の位相だけが残る
光の位相から読む
重力の中で、原子の位置は r ( t ) = r 0 + v 0 t + 1 2 g t 2 \bm{r}(t)=\bm{r}_0+\bm{v}_0t+\frac12\bm{g}t^2 r ( t ) = r 0 + v 0 t + 2 1 g t 2 である。式 (29.5) を式 (29.6) の Δ ϕ \Delta\phi Δ ϕ に入れると、2階差分は r ( 0 ) − 2 r ( T ) + r ( 2 T ) = g T 2 \bm{r}(0)-2\bm{r}(T)+\bm{r}(2T)=\bm{g}T^2 r ( 0 ) − 2 r ( T ) + r ( 2 T ) = g T 2 なので、光の位相の寄与は k e f f ⋅ g T 2 \bm{k}_{\mathrm{eff}}\cdot\bm{g}T^2 k eff ⋅ g T 2 である。落ちる原子にとって、2光子の共鳴の周波数は Doppler のずれ k e f f ⋅ g t \bm{k}_{\mathrm{eff}}\cdot\bm{g}t k eff ⋅ g t で時間とともに変わる。各パルスで共鳴を保つため、光の周波数の差を ω 1 − ω 2 = 一定 + α c h t \omega_1-\omega_2=\text{一定}+\alpha_{\mathrm{ch}}t ω 1 − ω 2 = 一定 + α ch t のように一定の速さ α c h \alpha_{\mathrm{ch}} α ch で掃引すると、Φ ( t ) = α c h t 2 / 2 \Phi(t)=\alpha_{\mathrm{ch}}t^2/2 Φ ( t ) = α ch t 2 /2 の2階差分 α c h T 2 \alpha_{\mathrm{ch}}T^2 α ch T 2 が引かれる。合わせて
Δ ϕ = ( k e f f ⋅ g − α c h ) T 2 (29.7) \Delta\phi=\bigl(\bm{k}_{\mathrm{eff}}\cdot\bm{g}-\alpha_{\mathrm{ch}}\bigr)T^2\tag{29.7} Δ ϕ = ( k eff ⋅ g − α ch ) T 2 ( 29.7 )
を得る。掃引の速さ α c h \alpha_{\mathrm{ch}} α ch を変えながら出口の確率を測ると、α c h = k e f f ⋅ g \alpha_{\mathrm{ch}}=\bm{k}_{\mathrm{eff}}\cdot\bm{g} α ch = k eff ⋅ g で Δ ϕ = 0 \Delta\phi=0 Δ ϕ = 0 になり、そこで P s P_s P s が最小になる。この点は T T T によらないので、いくつかの T T T で縞を取ると、すべての縞が共通に暗くなる点として k e f f ⋅ g \bm{k}_{\mathrm{eff}}\cdot\bm{g} k eff ⋅ g が決まる(図 29.3 )。重力加速度は、光の周波数の掃引の速さと光の波数という、どちらも周波数の標準で測れる2つの量の比 g = ∣ α c h ∣ / k e f f g=\lvert \alpha_{\mathrm{ch}}\rvert /k_{\mathrm{eff}} g = ∣ α ch ∣ / k eff として読み出される。
図 29.3 掃引の速さ α c h \alpha_{\mathrm{ch}} α ch を変えたときの出口の確率 P s = ( 1 − cos Δ ϕ ) / 2 P_s=(1-\cos\Delta\phi)/2 P s = ( 1 − cos Δ ϕ ) /2 (式 (29.7) 、コントラスト 1)。赤の実線は T = 10 m s T=10\,\mathrm{ms} T = 10 ms 、青の点線は T = 20 m s T=20\,\mathrm{ms} T = 20 ms 。縞の間隔は 1 / T 2 1/T^2 1/ T 2 に比例して狭くなり、α c h = k e f f ⋅ g \alpha_{\mathrm{ch}}=\bm{k}_{\mathrm{eff}}\cdot\bm{g} α ch = k eff ⋅ g の点だけはどの T T T でも暗い
経路積分の位相の3つの寄与
上の計算は、光の位相を原子の古典的な軌道の上で評価しただけで、2本の経路を進む間に波束が受け取る位相を数えていない。この点を経路積分で確かめる。解析力学から様々な理論分野を組み立てる 第4部第8章(解析力学から様々な理論分野を組み立てる 第4部第8章 )で示したとおり、Lagrangian が位置と速度の2次以下なら、伝播の核は古典の経路の作用 S c l S_{\mathrm{cl}} S cl を使って e i S c l / ℏ \mathrm{e}^{\mathrm{i} S_{\mathrm{cl}}/\hbar} e i S cl /ℏ に前因子を掛けた形に厳密に書ける。重力の中の自由な原子の Lagrangian L = 1 2 m r ˙ 2 + m g ⋅ r L=\frac12m\dot{\bm{r}}^2+m\bm{g}\cdot\bm{r} L = 2 1 m r ˙ 2 + m g ⋅ r はこの形である。
経路積分を2本の経路に当てはめた結果を、証明を省いて引く1 1 P. Storey and C. Cohen-Tannoudji, “The Feynman path integral approach to atomic interferometry. A tutorial, J. Phys. II France 4 , 1999 (1994). 。
命題 29.1(干渉計の位相の3つの寄与)
Lagrangian が位置と速度の2次以下で、パルスが瞬間的なとき、Mach–Zehnder 型の干渉計の位相差は
Δ ϕ = 1 ℏ ( ∫ 上 L d t − ∫ 下 L d t ) ⏟ 伝搬 + ( ϕ 1 − ϕ 2 上 + ϕ 3 上 ) − ϕ 2 下 ⏟ 光 + p ˉ ⋅ ( r 下 − r 上 ) ℏ ⏟ 分離 (29.8) \Delta\phi=\underbrace{\frac{1}{\hbar}\Bigl(\int_{\text{上}}L\,\mathrm{d} t-\int_{\text{下}}L\,\mathrm{d} t\Bigr)}_{\text{伝搬}}
+\underbrace{\bigl(\phi_1-\phi_2^{\text{上}}+\phi_3^{\text{上}}\bigr)-\phi_2^{\text{下}}}_{\text{光}}
+\underbrace{\frac{\bar{\bm{p}}\cdot(\bm{r}_{\text{下}}-\bm{r}_{\text{上}})}{\hbar}}_{\text{分離}}\tag{29.8} Δ ϕ = 伝搬 ℏ 1 ( ∫ 上 L d t − ∫ 下 L d t ) + 光 ( ϕ 1 − ϕ 2 上 + ϕ 3 上 ) − ϕ 2 下 + 分離 ℏ p ˉ ⋅ ( r 下 − r 上 ) ( 29.8 )
の3つの和になる。伝搬の項は各経路の古典の作用(内部状態のエネルギー − ℏ ω s g -\hbar\omega_{sg} − ℏ ω s g を ∣ s ⟩ | s \rangle ∣ s ⟩ にいる間だけ L L L に足す)、光の項は各経路が受けたパルスの位相をその経路の位置で評価したもの、分離の項は最後のパルスで2本の経路の端点がずれていたときの補正(p ˉ \bar{\bm{p}} p ˉ は2本の運動量の平均)である。
一様な重力のとき、3つを1つずつ数える。上の経路は時刻 0 0 0 に速度を ℏ k e f f / m \hbar\bm{k}_{\mathrm{eff}}/m ℏ k eff / m だけ増やし、時刻 T T T に元へ戻す。下の経路は時刻 T T T に同じだけ増やす。図 29.2 のとおり、2本は時刻 2 T 2T 2 T に同じ点で閉じるので、分離の項は 0 0 0 である。伝搬の項は、前半で上の経路が速い分の運動エネルギーと位置エネルギーの差と、後半で下の経路が速い分の差が、ちょうど打ち消し合って 0 0 0 になる(29.3 の例題の (3) で確かめる)。内部状態のエネルギーも、2本がそれぞれ時間 T T T ずつ ∣ s ⟩ | s \rangle ∣ s ⟩ にいるので打ち消し合う。光の項は、各経路の位置で光の位相を評価すると、上の経路の位置が時刻 T T T に下の経路より ℏ k e f f T / m \hbar k_{\mathrm{eff}}T/m ℏ k eff T / m だけずれていることを含めても k e f f ⋅ g T 2 \bm{k}_{\mathrm{eff}}\cdot\bm{g}T^2 k eff ⋅ g T 2 になる。したがって式 (29.7) は、光の位相を古典の軌道で評価した近似でなく、一様な重力では厳密な結果である。
この打ち消し合いは、位相の別の書き方からも分かる。2本の経路の高さの差を Δ z ( t ) \Delta z(t) Δ z ( t ) とすると、位置エネルギーの差 m g Δ z ( t ) mg\Delta z(t) m g Δ z ( t ) が時間 d t \mathrm{d} t d t の間に位相 m g Δ z d t / ℏ mg\Delta z\,\mathrm{d} t/\hbar m g Δ z d t /ℏ を与える。2本の経路が ( z , t ) (z,t) ( z , t ) 面で囲む面積は ∫ Δ z d t = ℏ k e f f T 2 / m \int\Delta z\,\mathrm{d} t=\hbar k_{\mathrm{eff}}T^2/m ∫ Δ z d t = ℏ k eff T 2 / m なので
∣ Δ ϕ ∣ = m g ℏ ∫ Δ z d t = m g ℏ ⋅ ℏ k e f f T 2 m = k e f f g T 2 (29.9) \lvert \Delta\phi\rvert =\frac{mg}{\hbar}\int\Delta z\,\mathrm{d} t=\frac{mg}{\hbar}\cdot\frac{\hbar k_{\mathrm{eff}}T^2}{m}=k_{\mathrm{eff}}gT^2\tag{29.9} ∣ Δ ϕ ∣ = ℏ m g ∫ Δ z d t = ℏ m g ⋅ m ℏ k eff T 2 = k eff g T 2 ( 29.9 )
である。原子干渉計の重力の位相は、2本の経路が時空で囲む面積に比例する。運動エネルギーの差の寄与と光の位相の寄与の間の入れ替わりは書き方の違いで、合計はどちらで数えても同じになる。
重力の勾配
地球の重力は高さとともに弱くなる。鉛直上向きに z z z をとり、z ¨ = − g + γ g z \ddot z=-g+\gamma_{\mathrm g}z z ¨ = − g + γ g z (γ g = 2 g / R E ≈ 3.08 × 10 − 6 s − 2 \gamma_{\mathrm g}=2g/R_{\mathrm E}\approx3.08\times10^{-6}\,\mathrm s^{-2} γ g = 2 g / R E ≈ 3.08 × 1 0 − 6 s − 2 、R E R_{\mathrm E} R E は地球の半径)とすると、γ g \gamma_{\mathrm g} γ g の1次までで
z ( 0 ) − 2 z ( T ) + z ( 2 T ) = − g T 2 + γ g ( z 0 T 2 + v 0 T 3 − 7 12 g T 4 ) (29.10) z(0)-2z(T)+z(2T)=-gT^2+\gamma_{\mathrm g}\Bigl(z_0T^2+v_0T^3-\frac{7}{12}gT^4\Bigr)\tag{29.10} z ( 0 ) − 2 z ( T ) + z ( 2 T ) = − g T 2 + γ g ( z 0 T 2 + v 0 T 3 − 12 7 g T 4 ) ( 29.10 )
になる。z 0 z_0 z 0 , v 0 v_0 v 0 ははじめの位置と速度(上向きが正)である。最後の項は、原子が落ちて重力の強い所へ入るほど加速が増えることの効果で、静止から落とした原子では位相を相対的に 7 12 γ g T 2 \frac{7}{12}\gamma_{\mathrm g}T^2 12 7 γ g T 2 だけ増やす。勾配があると Lagrangian に 1 2 m γ g z 2 \frac12m\gamma_{\mathrm g}z^2 2 1 m γ g z 2 が加わる。これも2次なので式 (29.8) は使えるが、2本の経路は時刻 2 T 2T 2 T に完全には閉じず、分離の項が残る。その大きさは ℏ k e f f / m \hbar k_{\mathrm{eff}}/m ℏ k eff / m に比例する小さい補正で、ここでは扱わない。
例題 29.3 ルビジウムの原子重力計 ★★★
オリジナル 静止から落とした 87 R b {}^{87}\mathrm{Rb} 87 Rb の雲で、k e f f \bm{k}_{\mathrm{eff}} k eff を鉛直に向けた Mach–Zehnder 型の干渉計を T = 0.1 s T=0.1\,\mathrm s T = 0.1 s で動かす。
(1) 掃引をしないときの位相 k e f f g T 2 k_{\mathrm{eff}}gT^2 k eff g T 2 、縞を暗く保つ掃引の速さ ∣ α c h ∣ / 2 π \lvert \alpha_{\mathrm{ch}}\rvert /2\pi ∣ α ch ∣ /2 π 、時刻 2 T 2T 2 T までに雲が落ちる距離を求めよ。
(2) 1回の測定で原子 N = 10 6 N=10^6 N = 1 0 6 個を数え、位相の不確かさが δ ϕ = 1 / N \delta\phi=1/\sqrt N δ ϕ = 1/ N のとき、g g g の相対的な不確かさ δ g / g \delta g/g δ g / g を求めよ。それは地表の高さにして何 m m \mathrm{mm} mm の変化に当たるか(g ∝ 1 / r 2 g\propto1/r^2 g ∝ 1/ r 2 、地球の半径 6.371 × 10 6 m 6.371\times10^{6}\,\mathrm m 6.371 × 1 0 6 m )。
(3) 一様な重力のとき、2本の経路の伝搬の位相の差 Δ S / ℏ \Delta S/\hbar Δ S /ℏ が 0 0 0 になることを示せ。
(4) 重力の勾配による補正の位相を求め、(2) の不確かさと比べよ。
解答 (1) k e f f g T 2 = ( 1.6106 × 10 7 ) ( 9.81 ) ( 0.1 ) 2 = 1.58 × 10 6 r a d k_{\mathrm{eff}}gT^2=(1.6106\times10^{7})(9.81)(0.1)^2=1.58\times10^{6}\,\mathrm{rad} k eff g T 2 = ( 1.6106 × 1 0 7 ) ( 9.81 ) ( 0.1 ) 2 = 1.58 × 1 0 6 rad である。掃引の速さの大きさは k e f f g = 1.58 × 10 8 r a d s − 2 k_{\mathrm{eff}}g=1.58\times10^{8}\,\mathrm{rad\,s^{-2}} k eff g = 1.58 × 1 0 8 rad s − 2 で、2 π 2\pi 2 π で割って 25.1 M H z s − 1 25.1\,\mathrm{MHz\,s^{-1}} 25.1 MHz s − 1 になる。落ちる距離は g ( 2 T ) 2 / 2 = 0.196 m g(2T)^2/2=0.196\,\mathrm m g ( 2 T ) 2 /2 = 0.196 m で、装置の高さは 20 c m 20\,\mathrm{cm} 20 cm ほどあればよい。
(2) δ g / g = δ ϕ / ( k e f f g T 2 ) = 10 − 3 / ( 1.58 × 10 6 ) = 6.33 × 10 − 10 \delta g/g=\delta\phi/(k_{\mathrm{eff}}gT^2)=10^{-3}/(1.58\times10^{6})=6.33\times10^{-10} δ g / g = δ ϕ / ( k eff g T 2 ) = 1 0 − 3 / ( 1.58 × 1 0 6 ) = 6.33 × 1 0 − 10 、すなわち δ g = 6.2 × 10 − 9 m s − 2 \delta g=6.2\times10^{-9}\,\mathrm{m\,s^{-2}} δ g = 6.2 × 1 0 − 9 m s − 2 である。g ∝ 1 / r 2 g\propto1/r^2 g ∝ 1/ r 2 から δ g / g = 2 δ r / R E \delta g/g=2\delta r/R_{\mathrm E} δ g / g = 2 δ r / R E なので、δ r = ( δ g / g ) R E / 2 = 2.02 m m \delta r=(\delta g/g)R_{\mathrm E}/2=2.02\,\mathrm{mm} δ r = ( δ g / g ) R E /2 = 2.02 mm である。1回の測定で、机の上と床の上の重力の違いの 1 / 400 1/400 1/400 程度の差が見分けられる。
(3) k e f f \bm{k}_{\mathrm{eff}} k eff を上向きにとり、上の経路を z 上 ( t ) = z 0 + ( v 0 + u ) t − g t 2 / 2 z_{\text{上}}(t)=z_0+(v_0+u)t-gt^2/2 z 上 ( t ) = z 0 + ( v 0 + u ) t − g t 2 /2 (0 ≤ t ≤ T 0\le t\le T 0 ≤ t ≤ T )、z 0 + v 0 t − g t 2 / 2 + u T z_0+v_0t-gt^2/2+uT z 0 + v 0 t − g t 2 /2 + u T (T ≤ t ≤ 2 T T\le t\le2T T ≤ t ≤ 2 T )、下の経路を z 0 + v 0 t − g t 2 / 2 z_0+v_0t-gt^2/2 z 0 + v 0 t − g t 2 /2 (0 ≤ t ≤ T 0\le t\le T 0 ≤ t ≤ T )、z 0 + v 0 t − g t 2 / 2 + u ( t − T ) z_0+v_0t-gt^2/2+u(t-T) z 0 + v 0 t − g t 2 /2 + u ( t − T ) (T ≤ t ≤ 2 T T\le t\le2T T ≤ t ≤ 2 T )と書く(u = ℏ k e f f / m u=\hbar k_{\mathrm{eff}}/m u = ℏ k eff / m )。前半の Lagrangian の差は L 上 − L 下 = m u ( v 0 − g t ) + 1 2 m u 2 − m g u t L_{\text{上}}-L_{\text{下}}=mu(v_0-gt)+\frac12mu^2-mgut L 上 − L 下 = m u ( v 0 − g t ) + 2 1 m u 2 − m g u t で、0 0 0 から T T T まで積分すると m u ( v 0 T + 1 2 u T − g T 2 ) mu\bigl(v_0T+\frac12uT-gT^2\bigr) m u ( v 0 T + 2 1 u T − g T 2 ) である。後半の差は − m u ( v 0 − g t ) − 1 2 m u 2 − m g u ( 2 T − t ) -mu(v_0-gt)-\frac12mu^2-mgu(2T-t) − m u ( v 0 − g t ) − 2 1 m u 2 − m g u ( 2 T − t ) で、T T T から 2 T 2T 2 T まで積分すると − m u ( v 0 T + 1 2 u T − 3 2 g T 2 ) − 1 2 m g u T 2 -mu\bigl(v_0T+\frac12uT-\frac32gT^2\bigr)-\frac12mguT^2 − m u ( v 0 T + 2 1 u T − 2 3 g T 2 ) − 2 1 m g u T 2 である。2つを足すと m u ( − g T 2 + 3 2 g T 2 − 1 2 g T 2 ) = 0 mu(-gT^2+\frac32gT^2-\frac12gT^2)=0 m u ( − g T 2 + 2 3 g T 2 − 2 1 g T 2 ) = 0 になる。
(4) 式 (29.10) で z 0 = v 0 = 0 z_0=v_0=0 z 0 = v 0 = 0 として、補正の相対的な大きさは 7 12 γ g T 2 = 7 12 ( 3.08 × 10 − 6 ) ( 0.01 ) = 1.80 × 10 − 8 \frac{7}{12}\gamma_{\mathrm g}T^2=\frac{7}{12}(3.08\times10^{-6})(0.01)=1.80\times10^{-8} 12 7 γ g T 2 = 12 7 ( 3.08 × 1 0 − 6 ) ( 0.01 ) = 1.80 × 1 0 − 8 、位相にして 1.58 × 10 6 × 1.80 × 10 − 8 = 0.0284 r a d 1.58\times10^{6}\times1.80\times10^{-8}=0.0284\,\mathrm{rad} 1.58 × 1 0 6 × 1.80 × 1 0 − 8 = 0.0284 rad である。これは (2) の不確かさ 10 − 3 r a d 10^{-3}\,\mathrm{rad} 1 0 − 3 rad の約 28 28 28 倍なので、補正せずに g g g を求めると誤る。原子を上へ打ち上げる配置(噴水。頂点で v 0 = g T v_0=gT v 0 = g T )では v 0 T 3 v_0T^3 v 0 T 3 の項が加わり、補正は − 5 12 γ g T 2 -\frac{5}{12}\gamma_{\mathrm g}T^2 − 12 5 γ g T 2 と符号が変わる。
1 P. Storey and C. Cohen-Tannoudji, “The Feynman path integral approach to atomic interferometry. A tutorial, J. Phys. II France 4 , 1999 (1994).
AI(Claude)が書いた教材を、木下が勉強しながら公開しています。式と数値は検算していますが、誤りが残っているかもしれません。誤りのご指摘・削除のご依頼をいただければ直します。 連絡先:support@tecktom.com