跳到主要内容

東京工業大学 理学院 物理学系 2019年8月実施 午前 [3]

Author​

祭音Myyura (co-authored with GPT 6 Astra)

Description​

以下の文章の空欄に当てはまる適切な数または数式を答えよ。解答には、プランク定数 hh ではなく、換算プランク定数 ℏ=h/(2π)\hbar=h/(2\pi) を用いること。

[A]​

質量 mm、角振動数 ω\omega の 1 次元調和振動子のラグランジアンは以下のように与えられる。

L=m2x˙2−mω22x2L=\frac{m}{2}\dot{x}^{2}-\frac{m\omega^{2}}{2}x^{2}

(1) 上のラグランジアンに対応するハミルトニアンは H(p,x)=p22m+mω22x2H(p,x)=\dfrac{p^{2}}{2m}+\dfrac{m\omega^{2}}{2}x^{2} である。正準量子化により xx、pp を演算子 x^\hat{x}、p^\hat{p} に置き換える。x^\hat{x}、p^\hat{p} は、正準交換関係 [x^,p^]=iℏ[\hat{x},\hat{p}]=i\hbar を満たす。生成演算子 a^†\hat{a}^{\dagger} 及び消滅演算子 a^\hat{a} を以下のように導入する。

a^†=12ℏ(mω x^−imωp^),a^=12ℏ(mω x^+imωp^)\begin{aligned} \hat{a}^{\dagger} &=\frac{1}{\sqrt{2\hbar}}\left(\sqrt{m\omega}\,\hat{x}-\frac{i}{\sqrt{m\omega}}\hat{p}\right),\\ \hat{a} &=\frac{1}{\sqrt{2\hbar}}\left(\sqrt{m\omega}\,\hat{x}+\frac{i}{\sqrt{m\omega}}\hat{p}\right) \end{aligned}

このとき [a^,a^†]=あ[\hat{a},\hat{a}^{\dagger}]=\boxed{\text{あ}} であり、H^\hat{H} を a^†\hat{a}^{\dagger}、a^\hat{a} を用いて表すと H^=い\hat{H}=\boxed{\text{い}} となる。また、a^†\hat{a}^{\dagger}、a^\hat{a} の交換関係を用いると、[H^,a^†]=う[\hat{H},\hat{a}^{\dagger}]=\boxed{\text{う}}、[H^,a^]=え[\hat{H},\hat{a}]=\boxed{\text{え}} となる。

(2) a^H†(t)\hat{a}^{\dagger}_{H}(t)、a^H(t)\hat{a}_{H}(t) をハイゼンベルク描像での生成演算子、消滅演算子とする。ここでハイゼンベルク描像での演算子 O^H(t)\hat{O}_{H}(t) とは、シュレーディンガー描像の演算子 O^S\hat{O}_{S} を、

O^H(t)=exp⁡(iℏH^t)O^Sexp⁡(−iℏH^t)\hat{O}_{H}(t) =\exp\left(\frac{i}{\hbar}\hat{H}t\right) \hat{O}_{S} \exp\left(-\frac{i}{\hbar}\hat{H}t\right)

のように時間発展させたものである。生成・消滅演算子に対するハイゼンベルク方程式は

ddta^H†(t)=お a^H†(t),ddta^H(t)=か a^H(t)\frac{d}{dt}\hat{a}^{\dagger}_{H}(t) =\boxed{\text{お}}\,\hat{a}^{\dagger}_{H}(t), \qquad \frac{d}{dt}\hat{a}_{H}(t) =\boxed{\text{か}}\,\hat{a}_{H}(t)

となる。初期条件 a^H†(0)=a^†\hat{a}^{\dagger}_{H}(0)=\hat{a}^{\dagger}、a^H(0)=a^\hat{a}_{H}(0)=\hat{a} を使うと、この方程式の解は

a^H†(t)=き,a^H(t)=く\hat{a}^{\dagger}_{H}(t)=\boxed{\text{き}}, \qquad \hat{a}_{H}(t)=\boxed{\text{く}}

となる。したがって、ハイゼンベルク描像での位置演算子 x^H(t)\hat{x}_{H}(t) は

x^H(t)=け a^†+こ a^\hat{x}_{H}(t)=\boxed{\text{け}}\,\hat{a}^{\dagger} +\boxed{\text{こ}}\,\hat{a}

となる。

[B]​

上の調和振動子に時刻 t=0t=0 で外部から撃力を加えたとする。対応するラグランジアンは以下のように与えられる。

L=m2x˙2−mω22x2+F(t)x,F(t)=F0δ(t)L=\frac{m}{2}\dot{x}^{2}-\frac{m\omega^{2}}{2}x^{2}+F(t)x, \qquad F(t)=F_{0}\delta(t)

ここで F0F_{0} は定数、δ(t)\delta(t) はディラックのデルタ関数である。

(3) ハイゼンベルク描像での位置演算子 x^H(t)\hat{x}_{H}(t)、運動量演算子 p^H(t)\hat{p}_{H}(t) に対するハイゼンベルク方程式から p^H(t)\hat{p}_{H}(t) を消去すると、

d2dt2x^H(t)=さ\frac{d^{2}}{dt^{2}}\hat{x}_{H}(t)=\boxed{\text{さ}}

となる。

(4) t=0t=0 を境にして、x^H(t)\hat{x}_{H}(t) は

lim⁡ϵ→+0(x^H(ϵ)−x^H(−ϵ))=し,lim⁡ϵ→+0(ddtx^H(ϵ)−ddtx^H(−ϵ))=す\begin{aligned} \lim_{\epsilon\to+0} \left(\hat{x}_{H}(\epsilon)-\hat{x}_{H}(-\epsilon)\right) &=\boxed{\text{し}},\\ \lim_{\epsilon\to+0} \left(\frac{d}{dt}\hat{x}_{H}(\epsilon) -\frac{d}{dt}\hat{x}_{H}(-\epsilon)\right) &=\boxed{\text{す}} \end{aligned}

を満たす。

(5) t<0t<0 では、ハイゼンベルク方程式の解は問 (2) で求めたように

x^H(t)=け a^†+こ a^\hat{x}_{H}(t)=\boxed{\text{け}}\,\hat{a}^{\dagger} +\boxed{\text{こ}}\,\hat{a}

と書くことができる。同様にして t>0t>0 でのハイゼンベルク方程式の解は t>0t>0 での生成・消滅演算子 b^†\hat{b}^{\dagger}、b^\hat{b} を用いて

x^H(t)=け b^†+こ b^\hat{x}_{H}(t)=\boxed{\text{け}}\,\hat{b}^{\dagger} +\boxed{\text{こ}}\,\hat{b}

と書くことができる。このとき問 (4) で求めた条件を b^†\hat{b}^{\dagger}、b^\hat{b} について解くと

b^†=せ,b^=そ\hat{b}^{\dagger}=\boxed{\text{せ}}, \qquad \hat{b}=\boxed{\text{そ}}

が得られる。

(6) t<0t<0 で考えている系が基底状態 ∣0⟩|0\rangle であったとする(即ち a^∣0⟩=0\hat{a}|0\rangle=0)。このとき t>0t>0 では

⟨0∣x^H(t)∣0⟩=た,⟨0∣(x^H(t)−た)2∣0⟩=ち\begin{aligned} \langle0|\hat{x}_{H}(t)|0\rangle&=\boxed{\text{た}},\\ \left\langle0\left| \left(\hat{x}_{H}(t)-\boxed{\text{た}}\right)^{2} \right|0\right\rangle&=\boxed{\text{ち}} \end{aligned}

となる。ここで、∣0⟩|0\rangle は ⟨0∣0⟩=1\langle0|0\rangle=1 と規格化されているものとする。

(7) t<0t<0 で考えている系が第 1 励起状態 ∣1⟩|1\rangle であったとする。このとき t>0t>0 では

⟨1∣x^H(t)∣1⟩=つ,⟨1∣(x^H(t)−つ)2∣1⟩=て\begin{aligned} \langle1|\hat{x}_{H}(t)|1\rangle&=\boxed{\text{つ}},\\ \left\langle1\left| \left(\hat{x}_{H}(t)-\boxed{\text{つ}}\right)^{2} \right|1\right\rangle&=\boxed{\text{て}} \end{aligned}

となる。ここで、∣1⟩|1\rangle は ⟨1∣1⟩=1\langle1|1\rangle=1 と規格化されているものとする。

题目描述​

在下文的空白处填入适当的数值或表达式。答案中须使用约化普朗克常数 ℏ=h/(2π)\hbar=h/(2\pi),而不使用普朗克常数 hh。

[A]​

质量为 mm、角频率为 ω\omega 的一维谐振子的拉格朗日量为

L=m2x˙2−mω22x2.L=\frac{m}{2}\dot x^2-\frac{m\omega^2}{2}x^2.

(1) 上述拉格朗日量对应的哈密顿量为 H(p,x)=p2/(2m)+mω2x2/2H(p,x)=p^2/(2m)+m\omega^2x^2/2。通过正则量子化,将 x,px,p 替换为算符 x^,p^\hat x,\hat p,它们满足正则对易关系 [x^,p^]=iℏ[\hat x,\hat p]=i\hbar。如下引入产生算符 a^†\hat a^\dagger 和湮灭算符 a^\hat a:

a^†=12ℏ(mωx^−imωp^),a^=12ℏ(mωx^+imωp^).\begin{aligned} \hat a^\dagger&=\frac{1}{\sqrt{2\hbar}} \left(\sqrt{m\omega}\hat x-\frac{i}{\sqrt{m\omega}}\hat p\right),\\ \hat a&=\frac{1}{\sqrt{2\hbar}} \left(\sqrt{m\omega}\hat x+\frac{i}{\sqrt{m\omega}}\hat p\right). \end{aligned}

此时 [a^,a^†]=あ[\hat a,\hat a^\dagger]=\boxed{\text{あ}};用 a^†,a^\hat a^\dagger,\hat a 表示 H^\hat H,可得 H^=い\hat H=\boxed{\text{い}}。此外,利用 a^†,a^\hat a^\dagger,\hat a 的对易关系可得

[H^,a^†]=う,[H^,a^]=え.[\hat H,\hat a^\dagger]=\boxed{\text{う}},\qquad [\hat H,\hat a]=\boxed{\text{え}}.

(2) 设 a^H†(t),a^H(t)\hat a_H^\dagger(t),\hat a_H(t) 分别为海森堡绘景中的产生算符和湮灭算符。这里,海森堡绘景中的算符 O^H(t)\hat O_H(t) 由薛定谔绘景中的算符 O^S\hat O_S 按照

O^H(t)=exp⁡(iℏH^t)O^Sexp⁡(−iℏH^t)\hat O_H(t)=\exp\left(\frac{i}{\hbar}\hat Ht\right) \hat O_S\exp\left(-\frac{i}{\hbar}\hat Ht\right)

作时间演化得到。产生、湮灭算符满足的海森堡方程为

ddta^H†(t)=おa^H†(t),ddta^H(t)=かa^H(t).\begin{aligned} \frac{d}{dt}\hat a_H^\dagger(t)&=\boxed{\text{お}}\hat a_H^\dagger(t),\\ \frac{d}{dt}\hat a_H(t)&=\boxed{\text{か}}\hat a_H(t). \end{aligned}

利用初始条件 a^H†(0)=a^†\hat a_H^\dagger(0)=\hat a^\dagger、a^H(0)=a^\hat a_H(0)=\hat a,这些方程的解为

a^H†(t)=き,a^H(t)=く.\hat a_H^\dagger(t)=\boxed{\text{き}},\qquad \hat a_H(t)=\boxed{\text{く}}.

因此,海森堡绘景中的位置算符 x^H(t)\hat x_H(t) 为

x^H(t)=けa^†+こa^.\hat x_H(t)=\boxed{\text{け}}\hat a^\dagger +\boxed{\text{こ}}\hat a.

[B]​

设上述谐振子在时刻 t=0t=0 受到外部冲击力。相应的拉格朗日量为

L=m2x˙2−mω22x2+F(t)x,F(t)=F0δ(t).\begin{aligned} L&=\frac{m}{2}\dot x^2-\frac{m\omega^2}{2}x^2+F(t)x,\\ F(t)&=F_0\delta(t). \end{aligned}

其中 F0F_0 是常数,δ(t)\delta(t) 是狄拉克 delta 函数。

(3) 从海森堡绘景中位置算符 x^H(t)\hat x_H(t)、动量算符 p^H(t)\hat p_H(t) 所满足的海森堡方程中消去 p^H(t)\hat p_H(t),可得

d2dt2x^H(t)=さ.\frac{d^2}{dt^2}\hat x_H(t)=\boxed{\text{さ}}.

(4) 在 t=0t=0 两侧,x^H(t)\hat x_H(t) 满足

lim⁡ϵ→+0(x^H(ϵ)−x^H(−ϵ))=し,lim⁡ϵ→+0(ddtx^H(ϵ)−ddtx^H(−ϵ))=す.\begin{aligned} \lim_{\epsilon\to+0}\bigl(\hat x_H(\epsilon)-\hat x_H(-\epsilon)\bigr) &=\boxed{\text{し}},\\ \lim_{\epsilon\to+0}\left(\frac{d}{dt}\hat x_H(\epsilon) -\frac{d}{dt}\hat x_H(-\epsilon)\right) &=\boxed{\text{す}}. \end{aligned}

(5) 当 t<0t<0 时,如第 (2) 问所求,海森堡方程的解可写为

x^H(t)=けa^†+こa^.\hat x_H(t)=\boxed{\text{け}}\hat a^\dagger +\boxed{\text{こ}}\hat a.

同样,t>0t>0 时海森堡方程的解可用 t>0t>0 的产生、湮灭算符 b^†,b^\hat b^\dagger,\hat b 写为

x^H(t)=けb^†+こb^.\hat x_H(t)=\boxed{\text{け}}\hat b^\dagger +\boxed{\text{こ}}\hat b.

此时,将第 (4) 问得到的条件对 b^†,b^\hat b^\dagger,\hat b 求解,得到

b^†=せ,b^=そ.\hat b^\dagger=\boxed{\text{せ}},\qquad \hat b=\boxed{\text{そ}}.

(6) 设 t<0t<0 时所考虑的系统处于基态 ∣0⟩|0\rangle(即 a^∣0⟩=0\hat a|0\rangle=0)。那么在 t>0t>0 时,

⟨0∣x^H(t)∣0⟩=た,⟨0∣(x^H(t)−た)2∣0⟩=ち.\begin{aligned} \langle0|\hat x_H(t)|0\rangle&=\boxed{\text{た}},\\ \langle0|\bigl(\hat x_H(t)-\boxed{\text{た}}\bigr)^2|0\rangle &=\boxed{\text{ち}}. \end{aligned}

其中 ∣0⟩|0\rangle 已归一化,即 ⟨0∣0⟩=1\langle0|0\rangle=1。

(7) 设 t<0t<0 时所考虑的系统处于第一激发态 ∣1⟩|1\rangle。那么在 t>0t>0 时,

⟨1∣x^H(t)∣1⟩=つ,⟨1∣(x^H(t)−つ)2∣1⟩=て.\begin{aligned} \langle1|\hat x_H(t)|1\rangle&=\boxed{\text{つ}},\\ \langle1|\bigl(\hat x_H(t)-\boxed{\text{つ}}\bigr)^2|1\rangle &=\boxed{\text{て}}. \end{aligned}

其中 ∣1⟩|1\rangle 已归一化,即 ⟨1∣1⟩=1\langle1|1\rangle=1。

Kai​

(1)​

正準交換関係 [x^,p^]=iℏ[\hat x,\hat p]=i\hbar より、

[a^,a^†]=12ℏ(−i[x^,p^]+i[p^,x^])=1,x^=ℏ2mω(a^†+a^),p^=imℏω2(a^†−a^).\begin{aligned} [\hat a,\hat a^\dagger] &=\frac{1}{2\hbar}\left(-i[\hat x,\hat p]+i[\hat p,\hat x]\right)=1,\\ \hat x&=\sqrt{\frac{\hbar}{2m\omega}}(\hat a^\dagger+\hat a),\\ \hat p&=i\sqrt{\frac{m\hbar\omega}{2}}(\hat a^\dagger-\hat a). \end{aligned}

これらをハミルトニアンに代入すると、

H^=ℏω2(a^†a^+a^a^†)=ℏω(a^†a^+12).\hat H=\frac{\hbar\omega}{2}(\hat a^\dagger\hat a+\hat a\hat a^\dagger) =\hbar\omega\left(\hat a^\dagger\hat a+\frac12\right).

さらに [AB,C]=A[B,C]+[A,C]B[AB,C]=A[B,C]+[A,C]B を用いて、

[H^,a^†]=ℏωa^†,[H^,a^]=−ℏωa^.[\hat H,\hat a^\dagger]=\hbar\omega\hat a^\dagger, \qquad [\hat H,\hat a]=-\hbar\omega\hat a.

したがって、

(あ)=1,(い)=ℏω(a^†a^+12),(う)=ℏωa^†,(え)=−ℏωa^.\boxed{\begin{aligned} \text{(あ)}&=1,\\ \text{(い)}&=\hbar\omega\left(\hat a^\dagger\hat a+\frac12\right),\\ \text{(う)}&=\hbar\omega\hat a^\dagger,\\ \text{(え)}&=-\hbar\omega\hat a. \end{aligned}}

(2)​

ハイゼンベルク方程式 O^˙H=(i/ℏ)[H^,O^H]\dot{\hat O}_H=(i/\hbar)[\hat H,\hat O_H] と (1) より、

da^H†dt=iωa^H†,da^Hdt=−iωa^H.\frac{d\hat a_H^\dagger}{dt}=i\omega\hat a_H^\dagger, \qquad \frac{d\hat a_H}{dt}=-i\omega\hat a_H.

初期条件を満たす解は

a^H†(t)=eiωta^†,a^H(t)=e−iωta^\hat a_H^\dagger(t)=e^{i\omega t}\hat a^\dagger, \qquad \hat a_H(t)=e^{-i\omega t}\hat a

なので、

x^H(t)=ℏ2mω(eiωta^†+e−iωta^).\hat x_H(t)=\sqrt{\frac{\hbar}{2m\omega}} \left(e^{i\omega t}\hat a^\dagger+e^{-i\omega t}\hat a\right).

よって、

(お)=iω,(か)=−iω,(き)=eiωta^†,(く)=e−iωta^,(け)=ℏ2mωeiωt,(こ)=ℏ2mωe−iωt.\boxed{\begin{aligned} \text{(お)}&=i\omega,&\text{(か)}&=-i\omega,\\ \text{(き)}&=e^{i\omega t}\hat a^\dagger,& \text{(く)}&=e^{-i\omega t}\hat a,\\ \text{(け)}&=\sqrt{\frac{\hbar}{2m\omega}}e^{i\omega t},& \text{(こ)}&=\sqrt{\frac{\hbar}{2m\omega}}e^{-i\omega t}. \end{aligned}}

(3)​

p=∂L/∂x˙=mx˙p=\partial L/\partial\dot x=m\dot x より、外力を含むハミルトニアンは

H^′(t)=p^22m+mω22x^2−F(t)x^.\hat H'(t)=\frac{\hat p^2}{2m}+\frac{m\omega^2}{2}\hat x^2-F(t)\hat x.

これに対応するハイゼンベルク方程式は

x^˙H=p^Hm,p^˙H=−mω2x^H+F(t).\dot{\hat x}_H=\frac{\hat p_H}{m}, \qquad \dot{\hat p}_H=-m\omega^2\hat x_H+F(t).

p^H\hat p_H を消去して、

(さ)=−ω2x^H(t)+F0mδ(t).\boxed{\text{(さ)}=-\omega^2\hat x_H(t)+\frac{F_0}{m}\delta(t)}.

(4)​

t=0t=0 における位置と速度の跳びを、それぞれ J=x^H(0+)−x^H(0−)J=\hat x_H(0^+)-\hat x_H(0^-)、K=x^˙H(0+)−x^˙H(0−)K=\dot{\hat x}_H(0^+)-\dot{\hat x}_H(0^-) とおく。分布微分では x^¨H\ddot{\hat x}_H に Jδ′(t)+Kδ(t)J\delta'(t)+K\delta(t) が現れる。一方、(3) の右辺には δ′\delta' の項がなく、δ\delta の係数は F0/mF_0/m である。係数比較から

(し)=0,(す)=F0m.\boxed{\text{(し)}=0,\qquad\text{(す)}=\frac{F_0}{m}}.

すなわち、位置は連続であり、運動量は p^H(0+)−p^H(0−)=F0\hat p_H(0^+)-\hat p_H(0^-)=F_0 だけ増加する。

(5)​

c=ℏ/(2mω)c=\sqrt{\hbar/(2m\omega)} とおく。(4) の位置の連続性と速度の跳びより、

b^†+b^=a^†+a^,iωc[(b^†−b^)−(a^†−a^)]=F0m.\begin{aligned} \hat b^\dagger+\hat b&=\hat a^\dagger+\hat a,\\ i\omega c\left[(\hat b^\dagger-\hat b)-(\hat a^\dagger-\hat a)\right] &=\frac{F_0}{m}. \end{aligned}

これらを連立して、

(せ)=a^†−iF02mℏω,(そ)=a^+iF02mℏω.\boxed{\begin{aligned} \text{(せ)}&=\hat a^\dagger-\frac{iF_0}{\sqrt{2m\hbar\omega}},\\ \text{(そ)}&=\hat a+\frac{iF_0}{\sqrt{2m\hbar\omega}}. \end{aligned}}

定数のずれなので [b^,b^†]=1[\hat b,\hat b^\dagger]=1 も保たれる。これを位置演算子に代入すると、t>0t>0 では

x^H(t)=X^(t)+xcl(t),X^(t)=c(eiωta^†+e−iωta^),xcl(t)=F0mωsin⁡ωt.\begin{aligned} \hat x_H(t)&=\hat X(t)+x_{\mathrm{cl}}(t),\\ \hat X(t)&=c\left(e^{i\omega t}\hat a^\dagger+e^{-i\omega t}\hat a\right),\\ x_{\mathrm{cl}}(t)&=\frac{F_0}{m\omega}\sin\omega t. \end{aligned}

(6)​

a^∣0⟩=0\hat a|0\rangle=0 とそのエルミート共役より ⟨0∣X^(t)∣0⟩=0\langle0|\hat X(t)|0\rangle=0 なので、

(た)=F0mωsin⁡ωt.\boxed{\text{(た)}=\frac{F_0}{m\omega}\sin\omega t}.

平均値を引くと定数項 xclx_{\mathrm{cl}} が消える。また、

X^(t)2=c2(e2iωt(a^†)2+e−2iωta^2+a^†a^+a^a^†).\begin{aligned} \hat X(t)^2=c^2\bigl(&e^{2i\omega t}(\hat a^\dagger)^2 +e^{-2i\omega t}\hat a^2\\ &+\hat a^\dagger\hat a+\hat a\hat a^\dagger\bigr). \end{aligned}

基底状態では最後の項の期待値だけが 11 となるため、

(ち)=⟨0∣X^(t)2∣0⟩=ℏ2mω.\boxed{\text{(ち)}=\langle0|\hat X(t)^2|0\rangle =\frac{\hbar}{2m\omega}}.

(7)​

a^∣1⟩=∣0⟩\hat a|1\rangle=|0\rangle、a^†∣1⟩=2∣2⟩\hat a^\dagger|1\rangle=\sqrt2|2\rangle および異なる数状態の直交性から、⟨1∣X^(t)∣1⟩=0\langle1|\hat X(t)|1\rangle=0。したがって、

(つ)=F0mωsin⁡ωt.\boxed{\text{(つ)}=\frac{F_0}{m\omega}\sin\omega t}.

さらに

⟨1∣a^2∣1⟩=⟨1∣(a^†)2∣1⟩=0,⟨1∣a^†a^∣1⟩=1,⟨1∣a^a^†∣1⟩=2\begin{aligned} \langle1|\hat a^2|1\rangle &=\langle1|(\hat a^\dagger)^2|1\rangle=0,\\ \langle1|\hat a^\dagger\hat a|1\rangle&=1,\\ \langle1|\hat a\hat a^\dagger|1\rangle&=2 \end{aligned}

を (6) の式に用いると、

(て)=⟨1∣X^(t)2∣1⟩=3ℏ2mω.\boxed{\text{(て)}=\langle1|\hat X(t)^2|1\rangle =\frac{3\hbar}{2m\omega}}.

撃力は位置の平均値に古典的振動 xcl(t)x_{\mathrm{cl}}(t) を加えるが、これらの状態における位置の分散は変えない。

Reference​