跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2014年9月実施 II-2

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

我々が観測する状態がハミルトニアンの固有状態と一致するとは限らない。例えば観測されるニュートリノの弱い相互作用の固有状態は、ニュートリノ一粒子状態のハミルトニアンの固有状態とは一致しない。このような場合に起こる量子振動現象を以下の二つの調和振動子からなる系の量子力学で考えてみよう。以下では、自然単位系 ℏ=c=1\hbar=c=1 を用いる。

以下の交換関係を満たす二組の生成消滅演算子 a^1,a^1†,a^2,a^2†\hat a_1,\hat a_1^\dagger,\hat a_2,\hat a_2^\dagger で表される調和振動子を考える。

[a^i,a^j†]=δij,[a^i,a^j]=0.(A)[\hat a_i,\hat a_j^\dagger]=\delta_{ij},\quad[\hat a_i,\hat a_j]=0.\tag{A}

ここで、i,j=1,2i,j=1,2、δij\delta_{ij} はクロネッカーのデルタである。系の真空 ∣0⟩|0\rangle は、

a^1∣0⟩=0,a^2∣0⟩=0(B)\hat a_1|0\rangle=0,\quad\hat a_2|0\rangle=0\tag{B}

で定義される。二つの調和振動子が相互作用している場合として、以下のハミルトニアンを考えよう。第2項を通じて二つの調和振動子は相互作用している。

H^=ωAa^1†a^1+ωB(a^1†a^2+a^2†a^1)+ωCa^2†a^2=(a^1†,a^2†)(ωAωBωBωC)(a^1a^2).(C)\begin{aligned} \hat H&=\omega_A\hat a_1^\dagger\hat a_1+\omega_B(\hat a_1^\dagger\hat a_2+\hat a_2^\dagger\hat a_1)+\omega_C\hat a_2^\dagger\hat a_2\\ &=(\hat a_1^\dagger,\hat a_2^\dagger) \begin{pmatrix}\omega_A&\omega_B\\\omega_B&\omega_C\end{pmatrix} \begin{pmatrix}\hat a_1\\\hat a_2\end{pmatrix}. \end{aligned}\tag{C}

ただし、ωA,ωB,ωC\omega_A,\omega_B,\omega_C は実数とする。

(1) 新しく消滅演算子を、

b^1=a^1cos⁡θ+a^2sin⁡θ,b^2=−a^1sin⁡θ+a^2cos⁡θ\begin{aligned} \hat b_1&=\hat a_1\cos\theta+\hat a_2\sin\theta,\\ \hat b_2&=-\hat a_1\sin\theta+\hat a_2\cos\theta \end{aligned}

のように定義する。このとき実パラメーター θ\theta を適当に選ぶことで、ハミルトニアンは、

H^=∑i=1,2ωib^i†b^i=(b^1†,b^2†)(ω100ω2)(b^1b^2)\hat H=\sum_{i=1,2}\omega_i\hat b_i^\dagger\hat b_i =(\hat b_1^\dagger,\hat b_2^\dagger) \begin{pmatrix}\omega_1&0\\0&\omega_2\end{pmatrix} \begin{pmatrix}\hat b_1\\\hat b_2\end{pmatrix}

の形に書ける。このとき、ω1,ω2,tan⁡2θ\omega_1,\omega_2,\tan2\theta を ωA,ωB,ωC\omega_A,\omega_B,\omega_C を用いて表せ。ただし、ω1≤ω2\omega_1\le\omega_2 とする。

(2) 時刻 t=0t=0 で与えられる初期状態が ∣α⟩=a^1†∣0⟩|\alpha\rangle=\hat a_1^\dagger|0\rangle であるとき、時刻 tt での状態ベクトルを ∣ψ(t)⟩|\psi(t)\rangle とし、時刻 tt で ∣β⟩=a^2†∣0⟩|\beta\rangle=\hat a_2^\dagger|0\rangle にある確率 P(α→β)(t)=∣⟨β∣ψ(t)⟩∣2P(\alpha\to\beta)(t)=|\langle\beta|\psi(t)\rangle|^2 を、ω1,ω2,θ\omega_1,\omega_2,\theta を用いて表せ。

(3) 前問では時刻 t=0t=0 に発生した状態 ∣α⟩|\alpha\rangle を時刻 tt に状態 ∣β⟩|\beta\rangle として観測する確率を計算したが、実際に観測される確率は、以下の時間間隔 [t−ΔT/2,t+ΔT/2][t-\Delta T/2,t+\Delta T/2] に対する時間平均

P‾(α→β)(t)=1ΔT∫t−ΔT/2t+ΔT/2dt′ P(α→β)(t′)(D)\overline P(\alpha\to\beta)(t)=\frac1{\Delta T}\int_{t-\Delta T/2}^{t+\Delta T/2}dt'\,P(\alpha\to\beta)(t')\tag{D}

であるとする。ただし、t>ΔTt>\Delta T とする。このとき P‾(α→β)(t)\overline P(\alpha\to\beta)(t) の表式を求めよ。

(4) 問 (3) で求めた P‾(α→β)(t)\overline P(\alpha\to\beta)(t) を ω2>ω1\omega_2>\omega_1 のときに ΔT\Delta T の関数として図示し、これの ΔT=0\Delta T=0 での値、及び ΔT≫1/(ω2−ω1)\Delta T\gg1/(\omega_2-\omega_1) での極限値を図に明記せよ。

题目描述​

我们观测到的状态不一定与哈密顿量的本征态一致。例如,观测到的中微子弱相互作用本征态与中微子单粒子哈密顿量的本征态并不相同。下面通过由两个谐振子组成的量子力学系统,考虑这种情况下发生的量子振荡现象。以下采用自然单位制 ℏ=c=1\hbar=c=1。

考虑由两组产生、湮灭算符 a^1,a^1†,a^2,a^2†\hat a_1,\hat a_1^\dagger,\hat a_2,\hat a_2^\dagger 表示的谐振子,算符满足对易关系

[a^i,a^j†]=δij,[a^i,a^j]=0.(A)[\hat a_i,\hat a_j^\dagger]=\delta_{ij},\quad[\hat a_i,\hat a_j]=0.\tag{A}

其中 i,j=1,2i,j=1,2,δij\delta_{ij} 是克罗内克符号。系统的真空态 ∣0⟩|0\rangle 定义为

a^1∣0⟩=0,a^2∣0⟩=0.(B)\hat a_1|0\rangle=0,\quad\hat a_2|0\rangle=0.\tag{B}

考虑两个谐振子存在相互作用的情况,哈密顿量为

H^=ωAa^1†a^1+ωB(a^1†a^2+a^2†a^1)+ωCa^2†a^2=(a^1†,a^2†)(ωAωBωBωC)(a^1a^2).(C)\begin{aligned} \hat H&=\omega_A\hat a_1^\dagger\hat a_1+\omega_B(\hat a_1^\dagger\hat a_2+\hat a_2^\dagger\hat a_1)+\omega_C\hat a_2^\dagger\hat a_2\\ &=(\hat a_1^\dagger,\hat a_2^\dagger) \begin{pmatrix}\omega_A&\omega_B\\\omega_B&\omega_C\end{pmatrix} \begin{pmatrix}\hat a_1\\\hat a_2\end{pmatrix}. \end{aligned}\tag{C}

两个谐振子通过第二项相互作用,且 ωA,ωB,ωC\omega_A,\omega_B,\omega_C 均为实数。

(1) 定义新的湮灭算符

b^1=a^1cos⁡θ+a^2sin⁡θ,b^2=−a^1sin⁡θ+a^2cos⁡θ.\begin{aligned} \hat b_1&=\hat a_1\cos\theta+\hat a_2\sin\theta,\\ \hat b_2&=-\hat a_1\sin\theta+\hat a_2\cos\theta. \end{aligned}

适当选择实参数 θ\theta,哈密顿量可以写成

H^=∑i=1,2ωib^i†b^i=(b^1†,b^2†)(ω100ω2)(b^1b^2).\hat H=\sum_{i=1,2}\omega_i\hat b_i^\dagger\hat b_i =(\hat b_1^\dagger,\hat b_2^\dagger) \begin{pmatrix}\omega_1&0\\0&\omega_2\end{pmatrix} \begin{pmatrix}\hat b_1\\\hat b_2\end{pmatrix}.

用 ωA,ωB,ωC\omega_A,\omega_B,\omega_C 表示 ω1,ω2,tan⁡2θ\omega_1,\omega_2,\tan2\theta,其中 ω1≤ω2\omega_1\le\omega_2。

(2) 时刻 t=0t=0 的初态为 ∣α⟩=a^1†∣0⟩|\alpha\rangle=\hat a_1^\dagger|0\rangle。设时刻 tt 的状态矢量为 ∣ψ(t)⟩|\psi(t)\rangle,用 ω1,ω2,θ\omega_1,\omega_2,\theta 表示时刻 tt 处于 ∣β⟩=a^2†∣0⟩|\beta\rangle=\hat a_2^\dagger|0\rangle 的概率 P(α→β)(t)=∣⟨β∣ψ(t)⟩∣2P(\alpha\to\beta)(t)=|\langle\beta|\psi(t)\rangle|^2。

(3) 上一问计算了时刻 t=0t=0 产生的状态 ∣α⟩|\alpha\rangle 在时刻 tt 被观测为 ∣β⟩|\beta\rangle 的概率。现设实际观测的概率为时间区间 [t−ΔT/2,t+ΔT/2][t-\Delta T/2,t+\Delta T/2] 内的时间平均:

P‾(α→β)(t)=1ΔT∫t−ΔT/2t+ΔT/2dt′ P(α→β)(t′).(D)\overline P(\alpha\to\beta)(t)=\frac1{\Delta T}\int_{t-\Delta T/2}^{t+\Delta T/2}dt'\,P(\alpha\to\beta)(t').\tag{D}

其中 t>ΔTt>\Delta T。求 P‾(α→β)(t)\overline P(\alpha\to\beta)(t) 的表达式。

(4) 在 ω2>ω1\omega_2>\omega_1 时,将问 (3) 所得 P‾(α→β)(t)\overline P(\alpha\to\beta)(t) 作为 ΔT\Delta T 的函数画图,并在图中明确标出 ΔT=0\Delta T=0 时的值和 ΔT≫1/(ω2−ω1)\Delta T\gg1/(\omega_2-\omega_1) 时的极限值。

Kai​

(1)​

D=(ωA−ωC)2+4ωB2D=\sqrt{(\omega_A-\omega_C)^2+4\omega_B^2} とおく。特性方程式より

ω1,2=ωA+ωC∓D2.\boxed{\omega_{1,2}=\frac{\omega_A+\omega_C\mp D}{2}}.

変換後の非対角項は ωBcos⁡2θ+(ωC−ωA)sin⁡2θ/2\omega_B\cos2\theta+(\omega_C-\omega_A)\sin2\theta/2 なので

tan⁡2θ=2ωBωA−ωC.\boxed{\tan2\theta=\frac{2\omega_B}{\omega_A-\omega_C}}.

固有値の順序を指定するには、D>0D>0 に対して

cos⁡2θ=ωC−ωAD,sin⁡2θ=−2ωBD\cos2\theta=\frac{\omega_C-\omega_A}{D},\qquad \sin2\theta=-\frac{2\omega_B}{D}

を選べばよい。D=0D=0 では角度は任意である。

(2)​

逆変換 a1†=b1†cos⁡θ−b2†sin⁡θa_1^\dagger=b_1^\dagger\cos\theta-b_2^\dagger\sin\theta により

∣ψ(t)⟩=(e−iω1tcos⁡θ b1†−e−iω2tsin⁡θ b2†)∣0⟩.|\psi(t)\rangle=(e^{-i\omega_1t}\cos\theta\,b_1^\dagger -e^{-i\omega_2t}\sin\theta\,b_2^\dagger)|0\rangle.

したがって

⟨β∣ψ(t)⟩=sin⁡θcos⁡θ(e−iω1t−e−iω2t),P(t)=sin⁡22θ sin⁡2Dt2.\begin{gathered} \langle\beta|\psi(t)\rangle=\sin\theta\cos\theta(e^{-i\omega_1t}-e^{-i\omega_2t}),\\ \boxed{P(t)=\sin^22\theta\,\sin^2\frac{Dt}{2}}. \end{gathered}

(3)​

A=12sin⁡22θA=\frac12\sin^22\theta と書くと P(t)=A(1−cos⁡Dt)P(t)=A(1-\cos Dt)。積分して

P‾(t)=A[1−cos⁡(Dt)sin⁡(DΔT/2)DΔT/2].\boxed{\overline P(t)=A\left[1-\cos(Dt)\frac{\sin(D\Delta T/2)}{D\Delta T/2}\right]}.

D=0D=0 では遷移確率はゼロ。

(4)​

s=DΔT/2s=D\Delta T/2 に対する曲線は A[1−cos⁡(Dt)sinc⁡s]A[1-\cos(Dt)\operatorname{sinc}s] である。

P‾∣ΔT=0=P(t),P‾⟶A=12sin⁡22θ.\boxed{\overline P|_{\Delta T=0}=P(t),\qquad \overline P\longrightarrow A=\frac12\sin^22\theta}.

s=nπs=n\pi(n=1,2,…n=1,2,\ldots)で高さ AA を通り、AA の上下で振幅が 1/ΔT1/\Delta T に比例して減衰する。cos⁡Dt>0\cos Dt>0 なら AA より下から出発し、負なら上下を反転した形、ゼロなら水平線となる。次の図に代表的な三つの場合を示す。

時間平均遷移確率

t>ΔTt>\Delta T の制限内で大きい時間窓を考えるには Dt≫1Dt\gg1 が必要である。

Reference​