Skip to content

平成31年度 東大院 物理学専攻 修士 物理学 解答

Raw

This content is not available in your language yet.

4時間で 4 問という構成の年度です。前半 3 問は量子力学・統計力学・力学の王道で、いずれも教科書の標準的な道具(生成消滅演算子、転送行列、ラグランジアンと保存量)を最後まで走らせる体力を問います。後半の選択 3 問は実験寄りで、第3問後半と第4問後半、第6問は数値評価まで求められるため、単位と桁の管理が得点を分けます。

問題分野主題
第1問量子力学生成消滅演算子と四次非調和項の摂動、準古典極限
第2問統計力学1次元2状態スピン鎖の分配関数・転送行列・相関長
第3問解析力学・天体物理重力2体問題と連星ブラックホールのインスパイラル
第4問光学・フーリエ解析マイケルソン干渉計によるフーリエ分光と CMB 測定
第5問素粒子・原子核反応断面積、不変質量、制動放射、運動量測定の誤差
第6問電磁気学・電気回路同軸ケーブルの伝送線路方程式と表皮効果

第1問から第3問が必答で、第4問から第6問のうち 1 問を選択する形式ですが、ここでは全 6 問の解答を載せます。

第1問 調和振動子と四次非調和項

Section titled “第1問 調和振動子と四次非調和項”

1次元調和振動子

H0=12mp2+mω22x2H_0 = \frac{1}{2m}p^2 + \frac{m\omega^2}{2}x^2

を出発点とし、[x,p]=i[x,p]=i\hbar のもとで生成消滅演算子

a=mω2(x+ipmω),a=mω2(xipmω),N=aaa = \sqrt{\frac{m\omega}{2\hbar}}\left(x+\frac{ip}{m\omega}\right),\qquad a^\dagger = \sqrt{\frac{m\omega}{2\hbar}}\left(x-\frac{ip}{m\omega}\right),\qquad N = a^\dagger a

を導入します。後半では V=λx4V=\lambda x^4λ>0\lambda>0)を加えた H=H0+VH=H_0+V を、摂動として、さらに摂動でない場合を準古典的に扱います。1次元束縛ポテンシャル系のハミルトニアンが対角化可能で縮退がないことは既知として使います。

定義式から x+ipmωx+\dfrac{ip}{m\omega}xipmωx-\dfrac{ip}{m\omega} の交換子を計算します。

[x+ipmω,xipmω]=imω[x,p]+imω[p,x]=imω(i)+imω(i)=2mω\left[x+\frac{ip}{m\omega},\,x-\frac{ip}{m\omega}\right] = -\frac{i}{m\omega}[x,p] + \frac{i}{m\omega}[p,x] = -\frac{i}{m\omega}(i\hbar) + \frac{i}{m\omega}(-i\hbar) = \frac{2\hbar}{m\omega}

前の因子 mω/(2)m\omega/(2\hbar) を掛けて

[a,a]=1.[a,a^\dagger] = 1 .

これを使えば

[N,a]=[aa,a]=a[a,a]+[a,a]a=a,[N,a]=a[a,a]+[a,a]a=a.[N,a] = [a^\dagger a,\,a] = a^\dagger[a,a] + [a^\dagger,a]a = -a,\qquad [N,a^\dagger] = a^\dagger[a,a^\dagger] + [a^\dagger,a^\dagger]a = a^\dagger .

つぎに NNx,px,p に戻します。

N=aa=mω2(xipmω)(x+ipmω)=mω2(x2+p2m2ω2+imω[x,p])=1ω(p22m+mω22x2)12N = a^\dagger a = \frac{m\omega}{2\hbar}\left(x-\frac{ip}{m\omega}\right)\left(x+\frac{ip}{m\omega}\right) = \frac{m\omega}{2\hbar}\left(x^2+\frac{p^2}{m^2\omega^2}+\frac{i}{m\omega}[x,p]\right) = \frac{1}{\hbar\omega}\left(\frac{p^2}{2m}+\frac{m\omega^2}{2}x^2\right)-\frac{1}{2}

したがって H0=ω(N+12)H_0 = \hbar\omega\left(N+\tfrac12\right) であり、答えは α=1\alpha=1β=12\beta=\tfrac12 です。

NN の固有値を ν\nu、対応する規格化された固有状態を ν|\nu\rangle とします。

(i) 非負性は内積の正定値性から直接出ます。

ν=νNν=νaaν=aν20\nu = \langle\nu|N|\nu\rangle = \langle\nu|a^\dagger a|\nu\rangle = \bigl\lVert a|\nu\rangle \bigr\rVert^2 \ge 0

等号成立は aν=0a|\nu\rangle=0 のときに限ります。

(ii) [N,a]=a[N,a]=-a より Naν=(aNa)ν=(ν1)aνNa|\nu\rangle = (aN-a)|\nu\rangle = (\nu-1)a|\nu\rangle なので、aνa|\nu\rangle はゼロでなければ固有値 ν1\nu-1 の固有状態です。これを繰り返すと akνa^k|\nu\rangle はゼロでなければ固有値 νk\nu-k の固有状態で、そのノルムは

akν2=ν(ν1)(νk+1)\bigl\lVert a^k|\nu\rangle \bigr\rVert^2 = \nu(\nu-1)\cdots(\nu-k+1)

となります(k=1k=1 が (i) の計算で、あとは aνk2=νk\lVert a\,|\nu-k\rangle\rVert^2=\nu-k を使った帰納法です)。ここで ν\nu が非負整数でないと仮定し、ν\nu を超えない最大の整数を n0n_0 とします。kn0+1k\le n_0+1 のとき νk+1νn0>0\nu-k+1 \ge \nu-n_0 > 0 なので上の積の因子はすべて正、つまり an0+1ν0a^{n_0+1}|\nu\rangle \neq 0 です。ところがこれは固有値 νn01\nu-n_0-1 の固有状態であり、1<νn01<0-1<\nu-n_0-1<0 は (i) に反します。よって ν\nu は非負整数でなければなりません。逆に ν=n\nu=n が非負整数なら an+1n2=n(n1)10=0\lVert a^{n+1}|n\rangle\rVert^2 = n(n-1)\cdots 1\cdot 0 = 0 となって連鎖がちょうど止まり、矛盾は生じません。

NN の固有値は 0,1,2,0,1,2,\dots に限られ、an2=n(N+1)n=n+1>0\lVert a^\dagger|n\rangle\rVert^2 = \langle n|(N+1)|n\rangle = n+1 > 0 から aa^\dagger を作用させれば固有値は 1 ずつ上がり続けます。したがって固有値の集合はちょうど {0,1,2,}\{0,1,2,\dots\} で、縮退がないので低い順に番号を付ければ nn 番目の固有値が nn です。H0=ω(N+12)H_0=\hbar\omega(N+\tfrac12) より

En=ω(n+12).E_n = \hbar\omega\left(n+\frac{1}{2}\right).

固有状態は ana^\dagger|n\rangle のノルムが n+1\sqrt{n+1} であることから n+1=an+1n|n+1\rangle = \dfrac{a^\dagger}{\sqrt{n+1}}|n\rangle(位相の自由度は規格化された固有状態の任意性として吸収します)。これを繰り返して

n=1n!(a)n0(n1).|n\rangle = \frac{1}{\sqrt{n!}}\left(a^\dagger\right)^n |0\rangle \qquad (n\ge 1).

(i) H=H0+VH=H_0+V の固有値問題を VV の次数で展開します。H0k=Ek(0)kH_0|k\rangle = E_k^{(0)}|k\rangleVkl=kVlV_{kl}=\langle k|V|l\rangle とし、摂動を受けた固有状態を n+n(1)+|n\rangle + |n^{(1)}\rangle + \cdots、固有値を En(0)+En(1)+E_n^{(0)}+E_n^{(1)}+\cdots と書きます。位相と規格化の自由度を使って nn(1)=0\langle n|n^{(1)}\rangle = 0 と取ります。1次の項を集めると

H0n(1)+Vn=En(0)n(1)+En(1)n.H_0|n^{(1)}\rangle + V|n\rangle = E_n^{(0)}|n^{(1)}\rangle + E_n^{(1)}|n\rangle .

左から n\langle n| を掛けると nH0=En(0)n\langle n|H_0 = E_n^{(0)}\langle n| より第1項が消えて

En(1)=Vnn.E_n^{(1)} = V_{nn}.

左から k\langle k|knk\ne n)を掛けると、n(1)=knck(1)k|n^{(1)}\rangle = \sum_{k\ne n}c_k^{(1)}|k\rangle として

ck(1)=VknEn(0)Ek(0).c_k^{(1)} = \frac{V_{kn}}{E_n^{(0)}-E_k^{(0)}} .

2次の項は

H0n(2)+Vn(1)=En(0)n(2)+En(1)n(1)+En(2)nH_0|n^{(2)}\rangle + V|n^{(1)}\rangle = E_n^{(0)}|n^{(2)}\rangle + E_n^{(1)}|n^{(1)}\rangle + E_n^{(2)}|n\rangle

で、左から n\langle n| を掛けて nn(1)=0\langle n|n^{(1)}\rangle=0 を使うと En(2)=nVn(1)E_n^{(2)} = \langle n|V|n^{(1)}\rangle、すなわち

En(2)=knVnkVknEn(0)Ek(0)=knVnk2En(0)Ek(0).E_n^{(2)} = \sum_{k\ne n}\frac{V_{nk}V_{kn}}{E_n^{(0)}-E_k^{(0)}} = \sum_{k\ne n}\frac{\lvert V_{nk}\rvert^2}{E_n^{(0)}-E_k^{(0)}} .

縮退がないので分母はどこでも 0 になりません。

(ii) x=(a+a)x = \ell\,(a+a^\dagger)2mω\ell \equiv \sqrt{\dfrac{\hbar}{2m\omega}} と書き、(a+a)40(a+a^\dagger)^4|0\rangle を順に作用させます。

(a+a)0=1,(a+a)1=0+22,(a+a)(0+22)=31+63,(a+a)(31+63)=30+622+264.\begin{aligned} (a+a^\dagger)|0\rangle &= |1\rangle,\\ (a+a^\dagger)|1\rangle &= |0\rangle + \sqrt{2}\,|2\rangle,\\ (a+a^\dagger)\left(|0\rangle+\sqrt{2}\,|2\rangle\right) &= 3|1\rangle + \sqrt{6}\,|3\rangle,\\ (a+a^\dagger)\left(3|1\rangle+\sqrt{6}\,|3\rangle\right) &= 3|0\rangle + 6\sqrt{2}\,|2\rangle + 2\sqrt{6}\,|4\rangle . \end{aligned}

よって x40=4(30+622+264)x^4|0\rangle = \ell^4\left(3|0\rangle+6\sqrt2\,|2\rangle+2\sqrt6\,|4\rangle\right) で、V=λx4V=\lambda x^4 の行列要素は

V00=3λ4,V20=62λ4,V40=26λ4,V_{00} = 3\lambda\ell^4,\qquad V_{20} = 6\sqrt{2}\,\lambda\ell^4,\qquad V_{40} = 2\sqrt{6}\,\lambda\ell^4,

それ以外の kk では Vk0=0V_{k0}=0 です(4=2/(4m2ω2)\ell^4 = \hbar^2/(4m^2\omega^2))。x4x^4 が偶パリティなので奇数の kk が落ち、(a+a)4(a+a^\dagger)^4 が生成演算子を 4 個までしか含まないので k5k\ge 5 も落ちます。V00=3λ4V_{00}=3\lambda\ell^4 は、幅 x20=2\langle x^2\rangle_0=\ell^2 のガウス分布に対する x4=3x22\langle x^4\rangle = 3\langle x^2\rangle^2 と一致します。

基底状態のエネルギーは E0(0)Ek(0)=kωE_0^{(0)}-E_k^{(0)} = -k\hbar\omega を使って

E0=ω2+3λ4+(62λ4)22ω+(26λ4)24ω+O(λ3)=ω2+3λ4λ28ω(722+244)+O(λ3)=ω2+324m2ω2λ2138m4ω5λ2+O(λ3).\begin{aligned} E_0 &= \frac{\hbar\omega}{2} + 3\lambda\ell^4 + \frac{(6\sqrt2\,\lambda\ell^4)^2}{-2\hbar\omega} + \frac{(2\sqrt6\,\lambda\ell^4)^2}{-4\hbar\omega} + O(\lambda^3)\\ &= \frac{\hbar\omega}{2} + 3\lambda\ell^4 - \frac{\lambda^2\ell^8}{\hbar\omega}\left(\frac{72}{2}+\frac{24}{4}\right) + O(\lambda^3)\\ &= \frac{\hbar\omega}{2} + \frac{3\hbar^2}{4m^2\omega^2}\,\lambda - \frac{21\hbar^3}{8m^4\omega^5}\,\lambda^2 + O(\lambda^3). \end{aligned}

答えは E0=ω2[1+32λm2ω3214(λm2ω3)2+]E_0 = \dfrac{\hbar\omega}{2}\left[1+\dfrac{3}{2}\dfrac{\hbar\lambda}{m^2\omega^3}-\dfrac{21}{4}\left(\dfrac{\hbar\lambda}{m^2\omega^3}\right)^2+\cdots\right] です。λ/(m2ω3)=λ4/(ω)×4\hbar\lambda/(m^2\omega^3) = \lambda\ell^4/(\hbar\omega)\times 4 の形なので括弧の中は無次元で、展開は λ4ω\lambda\ell^4 \ll \hbar\omega、つまり非調和項の典型的な大きさが準位間隔よりずっと小さいときに有効です。

nn が大きいと転回点 x0x_0 も大きくなり、λx4\lambda x^412mω2x2\tfrac12 m\omega^2x^2 の比 2λx02/(mω2)2\lambda x_0^2/(m\omega^2) は発散するので、高い準位では四次項が支配します。ボーア・ゾンマーフェルト量子化

pdx=2x0x02m(Eλx4)dx=2π(n+12),x0=(Eλ)1/4\oint p\,dx = 2\int_{-x_0}^{x_0}\sqrt{2m\left(E-\lambda x^4\right)}\,dx = 2\pi\hbar\left(n+\frac{1}{2}\right), \qquad x_0=\left(\frac{E}{\lambda}\right)^{1/4}

において x=x0ux=x_0u と置くと

pdx=22mEx0111u4du=c0mλ1/4E3/4,c0=22111u4du4.9\oint p\,dx = 2\sqrt{2mE}\,x_0\int_{-1}^{1}\sqrt{1-u^4}\,du = c_0\sqrt{m}\,\lambda^{-1/4}E^{3/4},\qquad c_0 = 2\sqrt{2}\int_{-1}^{1}\sqrt{1-u^4}\,du \simeq 4.9

となり、位相空間の面積は E3/4E^{3/4} に比例します。これを 2πn2\pi\hbar n と等置すると E3/4nE^{3/4}\propto n、すなわち

Enn4/3(En4/3λ1/3m2/3n4/3).E_n \propto n^{4/3}\qquad\left(E_n \sim \hbar^{4/3}\lambda^{1/3}m^{-2/3}\,n^{4/3}\right).

答えは nn4/34/3 乗です。一般に VxsV\propto |x|^s では同じ計算で Enn2s/(s+2)E_n\propto n^{2s/(s+2)} となり、s=2s=2n1n^1(調和振動子)、s=4s=4n4/3n^{4/3}ss\to\inftyn2n^2(無限井戸)に一致します。4/3λ1/3m2/3\hbar^{4/3}\lambda^{1/3}m^{-2/3} の次元は、[λ]=J/m4[\lambda]=\mathrm{J/m^4} を代入すればエネルギーになることが確かめられます。

LL 個のスピン変数 Si\mathbf{S}_ix^=(1,0,0)\hat x=(1,0,0)y^=(0,1,0)\hat y=(0,1,0) の 2 値をとり、開いた鎖のエネルギーが

E(S1,,SL)=Ji=1L1SiSi+1,J>0E(\mathbf{S}_1,\dots,\mathbf{S}_L) = -J\sum_{i=1}^{L-1}\mathbf{S}_i\cdot\mathbf{S}_{i+1},\qquad J>0

で与えられる古典系です。x^x^=y^y^=1\hat x\cdot\hat x=\hat y\cdot\hat y=1x^y^=0\hat x\cdot\hat y=0 なので、各ボンドは両端が同じ向きなら J-J、直交していれば 00 を寄せるだけです。逆温度は β=1/(kBT)\beta=1/(k_\mathrm{B}T)、分配関数は ZL(β)Z_L(\beta) とします。

L=2L=2 の配置は 4 通りです。

(x^,x^): E=J,(y^,y^): E=J,(x^,y^): E=0,(y^,x^): E=0.\begin{aligned} (\hat x,\hat x):&\ E=-J, & (\hat y,\hat y):&\ E=-J,\\ (\hat x,\hat y):&\ E=0, & (\hat y,\hat x):&\ E=0. \end{aligned}

したがって

Z2(β)=2eβJ+2=2(eβJ+1).Z_2(\beta) = 2e^{\beta J} + 2 = 2\left(e^{\beta J}+1\right).

L=3L=3 の配置は 8 通りで、エネルギーは隣接する平行ペアの個数で決まります。

(x^,x^,x^), (y^,y^,y^):E=2J(x^,x^,y^), (y^,x^,x^), (x^,y^,y^), (y^,y^,x^):E=J(x^,y^,x^), (y^,x^,y^):E=0\begin{aligned} (\hat x,\hat x,\hat x),\ (\hat y,\hat y,\hat y):&\quad E=-2J\\ (\hat x,\hat x,\hat y),\ (\hat y,\hat x,\hat x),\ (\hat x,\hat y,\hat y),\ (\hat y,\hat y,\hat x):&\quad E=-J\\ (\hat x,\hat y,\hat x),\ (\hat y,\hat x,\hat y):&\quad E=0 \end{aligned}

内訳は E=2JE=-2J が 2 通り、E=JE=-J が 4 通り、E=0E=0 が 2 通りです。よって

Z3(β)=2e2βJ+4eβJ+2=2(eβJ+1)2.Z_3(\beta) = 2e^{2\beta J} + 4e^{\beta J} + 2 = 2\left(e^{\beta J}+1\right)^2 .

L1L-1 個のスピンからなる鎖の各配置に、右端に SL\mathbf{S}_L を継ぎ足すことを考えます。追加されるボルツマン因子は SL=SL1\mathbf{S}_L=\mathbf{S}_{L-1} なら eβJe^{\beta J}、そうでなければ 11 で、その和 eβJ+1e^{\beta J}+1SL1\mathbf{S}_{L-1} の向きによらない定数です。したがって

ZL(β)=(eβJ+1)ZL1(β)(L3)Z_L(\beta) = \left(e^{\beta J}+1\right)Z_{L-1}(\beta)\qquad (L\ge 3)

が成り立ち、Z2=2(eβJ+1)Z_2 = 2(e^{\beta J}+1) から

ZL(β)=2(eβJ+1)L1.Z_L(\beta) = 2\left(e^{\beta J}+1\right)^{L-1}.

L=2,3L=2,3 で設問1・2 の結果を再現します。この式は、L1L-1 本のボンドがそれぞれ独立に「平行(重み eβJe^{\beta J})か直交(重み 11)」を選び、さらに左端のスピンの向きに 2 通りの自由度がある、という数え方に対応します。

FL(β)=1βlnZL(β)=1β[ln2+(L1)ln(1+eβJ)].F_L(\beta) = -\frac{1}{\beta}\ln Z_L(\beta) = -\frac{1}{\beta}\left[\ln 2 + (L-1)\ln\left(1+e^{\beta J}\right)\right].

内部エネルギーは

UL=βlnZL=(L1)JeβJ1+eβJU_L = -\frac{\partial}{\partial\beta}\ln Z_L = -(L-1)\frac{Je^{\beta J}}{1+e^{\beta J}}

なので、LL\to\infty でスピン 1 個あたりは

u(β)=limLULL=JeβJ1+eβJ=J1+eβJ.u(\beta) = \lim_{L\to\infty}\frac{U_L}{L} = -\frac{Je^{\beta J}}{1+e^{\beta J}} = -\frac{J}{1+e^{-\beta J}} .

高温極限 βJ0\beta J\to 0 では uJ/2u\to -J/2、低温極限 βJ\beta J\to\infty では uJu\to -J です。

出現確率との関係は次のように読めます。LL\to\infty ではスピン 1 個あたりのボンド数が 1 に近づき、各ボンドは独立に、両端が平行である確率

p=eβJeβJ+1=11+eβJp_\parallel = \frac{e^{\beta J}}{e^{\beta J}+1} = \frac{1}{1+e^{-\beta J}}

をもちます。平行なボンドのエネルギーは J-J、直交なボンドは 00 なので u=Jpu = -J p_\parallel です。高温極限では 1 本のボンドの 4 通りの配置がすべて等確率に現れ、そのうち 2 通りが平行なので p=1/2p_\parallel = 1/2、したがって u=J/2u=-J/2 になります。低温極限ではボルツマン因子 eβJe^{\beta J} が効いて平行な配置しか現れず、p1p_\parallel\to 1、つまり全スピンが揃った 2 つの基底状態(エネルギー J(L1)-J(L-1))に確率が集中して uJu\to -J となります。

S1=a^\mathbf{S}_1=\hat a に固定した和を L1L-1 番目まで取ったあと、右端に SL=b^\mathbf{S}_L=\hat b を継ぎ足す操作を書き下します。ボンドの重みを w(c^,b^)=eβJw(\hat c,\hat b)=e^{\beta J}c^=b^\hat c=\hat b)、w(c^,b^)=1w(\hat c,\hat b)=1c^b^\hat c\ne\hat b)と書くと

ZLPL(a^,b^)=c^[ZL1PL1(a^,c^)]w(c^,b^)Z_L\,P_L(\hat a,\hat b) = \sum_{\hat c}\left[Z_{L-1}P_{L-1}(\hat a,\hat c)\right]w(\hat c,\hat b)

であり、設問3 の ZL1/ZL=1/(eβJ+1)Z_{L-1}/Z_L = 1/(e^{\beta J}+1) を使えば式(2)の形になります。したがって

(Q(x^,x^)Q(x^,y^)Q(y^,x^)Q(y^,y^))=1eβJ+1(eβJ11eβJ),\begin{pmatrix} Q(\hat x,\hat x) & Q(\hat x,\hat y)\\ Q(\hat y,\hat x) & Q(\hat y,\hat y)\end{pmatrix} = \frac{1}{e^{\beta J}+1}\begin{pmatrix} e^{\beta J} & 1\\ 1 & e^{\beta J}\end{pmatrix},

すなわち Q(x^,x^)=Q(y^,y^)=eβJeβJ+1Q(\hat x,\hat x)=Q(\hat y,\hat y)=\dfrac{e^{\beta J}}{e^{\beta J}+1}Q(x^,y^)=Q(y^,x^)=1eβJ+1Q(\hat x,\hat y)=Q(\hat y,\hat x)=\dfrac{1}{e^{\beta J}+1} です。各行の和が 1 になっており、Q(a^,b^)Q(\hat a,\hat b)SL1=a^\mathbf{S}_{L-1}=\hat a を与えたときに SL=b^\mathbf{S}_L=\hat b となる条件付き確率だと解釈できます。

設問6 の関係を繰り返すと、PL=P2QL2P_L = P_2\,Q^{L-2} です。L=2L=2 では設問1 から

P2=12(eβJ+1)(eβJ11eβJ)=12QP_2 = \frac{1}{2\left(e^{\beta J}+1\right)}\begin{pmatrix} e^{\beta J} & 1\\ 1 & e^{\beta J}\end{pmatrix} = \frac{1}{2}Q

なので PL=12QL1P_L = \tfrac12 Q^{L-1} となります。QQ の固有値は固有ベクトル (1,1)(1,1) に対して 11(1,1)(1,-1) に対して

qeβJ1eβJ+1=tanhβJ2q \equiv \frac{e^{\beta J}-1}{e^{\beta J}+1} = \tanh\frac{\beta J}{2}

なので、対応する射影演算子を使って

QL1=12(1+qL11qL11qL11+qL1),PL(x^,x^)=PL(y^,y^)=1+qL14.Q^{L-1} = \frac{1}{2}\begin{pmatrix} 1+q^{L-1} & 1-q^{L-1}\\ 1-q^{L-1} & 1+q^{L-1}\end{pmatrix}, \qquad P_L(\hat x,\hat x) = P_L(\hat y,\hat y) = \frac{1+q^{L-1}}{4}.

a^b^\hat a\cdot\hat b は両者が同じ向きのときだけ 1 なので

S1SL=PL(x^,x^)+PL(y^,y^)=12[1+(eβJ1eβJ+1)L1]=12+12tanhL1 ⁣βJ2.\langle \mathbf{S}_1\cdot\mathbf{S}_L\rangle = P_L(\hat x,\hat x)+P_L(\hat y,\hat y) = \frac{1}{2}\left[1+\left(\frac{e^{\beta J}-1}{e^{\beta J}+1}\right)^{L-1}\right] = \frac{1}{2}+\frac{1}{2}\tanh^{L-1}\!\frac{\beta J}{2}.

L=2L=2eβJ/(eβJ+1)=pe^{\beta J}/(e^{\beta J}+1)=p_\parallel となり設問5 の結果と、β0\beta\to 01/21/2(無相関)、β\beta\to\infty11(完全に揃う)となり物理的な期待と一致します。

式(3)と比べれば A=B=12A=B=\tfrac12、かつ exp(1/ξ)=q\exp(-1/\xi)=q、すなわち

ξ=1lntanh(βJ/2)=1lncoth(βJ/2).\xi = -\frac{1}{\ln\tanh(\beta J/2)} = \frac{1}{\ln\coth(\beta J/2)} .

高温極限 βJ1\beta J\ll 1 では tanh(βJ/2)βJ/2\tanh(\beta J/2)\simeq \beta J/2 なので

ξ1ln(2/βJ),\xi \simeq \frac{1}{\ln\left(2/\beta J\right)} ,

これは 1 サイト分より短く、隣同士でしか相関が残らないことを表します。低温極限 βJ1\beta J\gg 1 では tanh(βJ/2)=1eβJ1+eβJ12eβJ\tanh(\beta J/2) = \dfrac{1-e^{-\beta J}}{1+e^{-\beta J}} \simeq 1-2e^{-\beta J} より lntanh(βJ/2)2eβJ\ln\tanh(\beta J/2)\simeq -2e^{-\beta J} なので

ξ12eβJ.\xi \simeq \frac{1}{2}e^{\beta J} .

相関長は T0T\to 0 で指数関数的に発散します。1 本のボンドを反転させる励起エネルギーが JJ で、そのような「ドメイン壁」の平均間隔が eβJ\sim e^{\beta J} になる、という描像です。有限温度では ξ\xi が有限なので、この 1 次元系には有限温度の秩序相がありません。

第3問 重力2体問題と連星ブラックホール

Section titled “第3問 重力2体問題と連星ブラックホール”

質量 m1,m2m_1,m_2 の 2 質点が相互の重力で束縛運動する系を、相対座標 r=r2r1\mathbf{r}=\mathbf{r}_2-\mathbf{r}_1 の軌道面極座標 (r,φ)(r,\varphi) で記述します。M=m1+m2M=m_1+m_2μ=m1m2/(m1+m2)\mu=m_1m_2/(m_1+m_2) として

L=12μ(drdt)2+12μr2(dφdt)2+GμMr\mathcal{L} = \frac{1}{2}\mu\left(\frac{dr}{dt}\right)^2 + \frac{1}{2}\mu r^2\left(\frac{d\varphi}{dt}\right)^2 + \frac{G\mu M}{r}

が与えられています。後半では重力波放射によるエネルギー損失を、ニュートン力学の枠内で半径 a(t)a(t) の準円軌道の断熱的な縮みとして扱い、E=GμM/(2a)E=-G\mu M/(2a) と放射率

LGW=32G45c5μ2M3a5L_{\mathrm{GW}} = \frac{32G^4}{5c^5}\frac{\mu^2M^3}{a^5}

を使います。μ,M\mu,M の時間変化は無視します。

Lr˙=μr˙\dfrac{\partial\mathcal L}{\partial \dot r} = \mu\dot rLr=μrφ˙2GμMr2\dfrac{\partial \mathcal L}{\partial r} = \mu r\dot\varphi^2 - \dfrac{G\mu M}{r^2} より

μd2rdt2=μr(dφdt)2GμMr2.\mu\frac{d^2r}{dt^2} = \mu r\left(\frac{d\varphi}{dt}\right)^2 - \frac{G\mu M}{r^2}.

φ\varphiL\mathcal L に現れないので

ddt(μr2dφdt)=0.\frac{d}{dt}\left(\mu r^2\frac{d\varphi}{dt}\right) = 0 .

L\mathcal Ltt を陽に含まないので、エネルギーは E=r˙L/r˙+φ˙L/φ˙LE=\dot r\,\partial\mathcal L/\partial\dot r + \dot\varphi\,\partial\mathcal L/\partial\dot\varphi - \mathcal L で与えられます。

E=12μ(drdt)2+12μr2(dφdt)2GμMr,J=Lφ˙=μr2dφdt.E = \frac{1}{2}\mu\left(\frac{dr}{dt}\right)^2 + \frac{1}{2}\mu r^2\left(\frac{d\varphi}{dt}\right)^2 - \frac{G\mu M}{r}, \qquad J = \frac{\partial\mathcal L}{\partial\dot\varphi} = \mu r^2\frac{d\varphi}{dt} .

JJ については、L\mathcal Lφ\varphi を含まないので

dJdt=ddtLφ˙=Lφ=0.\frac{dJ}{dt} = \frac{d}{dt}\frac{\partial\mathcal L}{\partial\dot\varphi} = \frac{\partial\mathcal L}{\partial\varphi} = 0 .

EE については直接微分して運動方程式を代入します。

dEdt=μr˙r¨+μrr˙φ˙2+μr2φ˙φ¨+GμMr2r˙\frac{dE}{dt} = \mu\dot r\ddot r + \mu r\dot r\dot\varphi^2 + \mu r^2\dot\varphi\ddot\varphi + \frac{G\mu M}{r^2}\dot r

に、設問1 の μr¨=μrφ˙2GμM/r2\mu\ddot r = \mu r\dot\varphi^2 - G\mu M/r^2 と、ddt(μr2φ˙)=0\dfrac{d}{dt}(\mu r^2\dot\varphi)=0 から従う μr2φ¨=2μrr˙φ˙\mu r^2\ddot\varphi = -2\mu r\dot r\dot\varphi を入れると

dEdt=(μrr˙φ˙2GμMr˙r2)+μrr˙φ˙22μrr˙φ˙2+GμMr˙r2=0.\frac{dE}{dt} = \left(\mu r\dot r\dot\varphi^2 - \frac{G\mu M\dot r}{r^2}\right) + \mu r\dot r\dot\varphi^2 - 2\mu r\dot r\dot\varphi^2 + \frac{G\mu M\dot r}{r^2} = 0 .

よって EEJJ はいずれも運動の定数です。前者は L\mathcal L の時間並進対称性、後者は軌道面内の回転対称性に対応します。

E=GμM2aE=-\dfrac{G\mu M}{2a} を時間で微分すると dEdt=GμM2a2dadt\dfrac{dE}{dt} = \dfrac{G\mu M}{2a^2}\dfrac{da}{dt} です。これを dE/dt=LGWdE/dt=-L_{\mathrm{GW}} に代入して

GμM2a2dadt=32G45c5μ2M3a5dadt=64G35c5μM2a3.\frac{G\mu M}{2a^2}\frac{da}{dt} = -\frac{32G^4}{5c^5}\frac{\mu^2M^3}{a^5} \quad\Longrightarrow\quad \frac{da}{dt} = -\frac{64G^3}{5c^5}\frac{\mu M^2}{a^3}.

右辺が負なので軌道は縮み続けます。a3a^3 が分母にあるため縮みは加速し、有限時間で a0a\to 0 に達します。

dPdt=A(PcP)5/3\dfrac{dP}{dt} = -A\left(\dfrac{P_\mathrm{c}}{P}\right)^{5/3} を変数分離して

P5/3dP=APc5/3dt38(P8/3P08/3)=APc5/3t.P^{5/3}\,dP = -A P_\mathrm{c}^{5/3}\,dt \quad\Longrightarrow\quad \frac{3}{8}\left(P^{8/3}-P_0^{8/3}\right) = -AP_\mathrm{c}^{5/3}\,t .

P=0P=0 となる時刻が合体時刻なので

τGW=3P08/38APc5/3=38APc(P0Pc)8/3.\tau_{\mathrm{GW}} = \frac{3\,P_0^{8/3}}{8A\,P_\mathrm{c}^{5/3}} = \frac{3}{8A}P_\mathrm{c}\left(\frac{P_0}{P_\mathrm{c}}\right)^{8/3}.

AA が無次元、PcP_\mathrm{c} が時間の次元なので右辺は時間の次元をもちます。

ケプラーの方程式 GMP2=4π2a3GMP^2=4\pi^2a^3 を微分すると 2GMPP˙=12π2a2a˙2GMP\,\dot P = 12\pi^2a^2\dot a、すなわち

dPdt=6π2a2GMPdadt.\frac{dP}{dt} = \frac{6\pi^2a^2}{GMP}\frac{da}{dt} .

設問4 の結果を入れると

dPdt=6π2a2GMP64G3μM25c5a3=384π2G2μM5c5Pa,\frac{dP}{dt} = -\frac{6\pi^2a^2}{GMP}\cdot\frac{64G^3\mu M^2}{5c^5a^3} = -\frac{384\pi^2G^2\mu M}{5c^5\,P\,a},

さらに a=(GM4π2)1/3P2/3a = \left(\dfrac{GM}{4\pi^2}\right)^{1/3}P^{2/3} を代入して

dPdt=384π2(4π2)1/35G5/3μM2/3c5P5/3=965(2π)8/3(Gμ3/5M2/5c3)5/3P5/3.\frac{dP}{dt} = -\frac{384\pi^2\left(4\pi^2\right)^{1/3}}{5}\,\frac{G^{5/3}\mu M^{2/3}}{c^5}\,P^{-5/3} = -\frac{96}{5}(2\pi)^{8/3}\left(\frac{G\mu^{3/5}M^{2/5}}{c^3}\right)^{5/3}P^{-5/3}.

これが式(6)です。比較すれば

Pc=Gμ3/5M2/5c3,α=35,β=25,γ=3,P_\mathrm{c} = \frac{G\mu^{3/5}M^{2/5}}{c^3},\qquad \alpha = \frac{3}{5},\quad \beta = \frac{2}{5},\quad \gamma = -3, A=38422/3π8/35=965(2π)8/32.6×103.A = \frac{384\cdot 2^{2/3}\pi^{8/3}}{5} = \frac{96}{5}(2\pi)^{8/3} \simeq 2.6\times 10^{3}.

PcP_\mathrm{c} の次元は、[G]=m3kg1s2[G]=\mathrm{m^3\,kg^{-1}\,s^{-2}} から m33kg1+3/5+2/5s2+3=s\mathrm{m^{3-3}\,kg^{-1+3/5+2/5}\,s^{-2+3}}=\mathrm{s} となり時間で正しく、α+β=1\alpha+\beta=1γ=3\gamma=-3 は次元解析だけからも決まります。μ3/5M2/5\mu^{3/5}M^{2/5} はチャープ質量と呼ばれる組み合わせで、重力波の周期変化から直接読み取れる質量はこれです。

設問5 の結果を PcP_\mathrm{c} について解きます。

Pc=(38AτGW)3/5P08/5P_\mathrm{c} = \left(\frac{3}{8A\tau_{\mathrm{GW}}}\right)^{3/5}P_0^{8/5}

等質量 mm の連星では M=2mM=2mμ=m/2\mu=m/2 なので μ3/5M2/5=(m/2)3/5(2m)2/5=21/5m\mu^{3/5}M^{2/5} = (m/2)^{3/5}(2m)^{2/5}=2^{-1/5}m、したがって

Pc=Gm21/5c3mM=21/5Pcc3GM.P_\mathrm{c} = \frac{Gm}{2^{1/5}c^3} \quad\Longrightarrow\quad \frac{m}{M_\odot} = 2^{1/5}\,P_\mathrm{c}\,\frac{c^3}{GM_\odot} .

数値を入れます。38AτGW=38×2500×0.15=1.0×103\dfrac{3}{8A\tau_{\mathrm{GW}}} = \dfrac{3}{8\times 2500\times 0.15} = 1.0\times 10^{-3} なのでその 3/53/5 乗は 101.810^{-1.8}、また P08/5=(6×102)8/5=663/5×103.26×2.9×103.2=17.4×103.2P_0^{8/5} = (6\times 10^{-2})^{8/5} = 6\cdot 6^{3/5}\times 10^{-3.2} \simeq 6\times 2.9\times 10^{-3.2}=17.4\times 10^{-3.2} です。よって

Pc17.4×101.83.2=1.7×104s.P_\mathrm{c} \simeq 17.4\times 10^{-1.8-3.2} = 1.7\times 10^{-4}\,\mathrm{s}.

2GM/c3105s2GM_\odot/c^3\approx 10^{-5}\,\mathrm{s} すなわち GM/c35×106sGM_\odot/c^3 \approx 5\times 10^{-6}\,\mathrm{s} を使って

mM1.1×1.74×1045×10638.\frac{m}{M_\odot} \simeq \frac{1.1\times 1.74\times 10^{-4}}{5\times 10^{-6}} \simeq 38 .

答えは m4×101Mm \approx 4\times 10^{1}\,M_\odot、太陽質量の 40 倍程度です。GW150914 の実際の推定値(約 30 から 36 太陽質量)と同じ桁に収まっています。

μ,M\mu,M を一定として、a0a_0 から amina_{\min} まで縮む間に系が失うエネルギー Erad=E(a0)E(amin)E_{\mathrm{rad}} = E(a_0)-E(a_{\min})

Erad=GμM2aminGμM2a0GμM2aminE_{\mathrm{rad}} = \frac{G\mu M}{2a_{\min}} - \frac{G\mu M}{2a_0} \simeq \frac{G\mu M}{2a_{\min}}

です(a0amina_0\gg a_{\min} を使いました)。aaamina_{\min} に達したあとは放射されないと仮定するので、これが放射された重力波の全エネルギーで、それが mm の約 1 割に相当する静止エネルギー 0.1mc20.1\,mc^2 に等しいとします。等質量なので μM=m22m=m2\mu M = \dfrac{m}{2}\cdot 2m = m^2

Gm22amin=mc210amin=5Gmc2.\frac{Gm^2}{2a_{\min}} = \frac{mc^2}{10} \quad\Longrightarrow\quad a_{\min} = \frac{5Gm}{c^2}.

これが求める関係式です。1 つのブラックホールのシュワルツシルト半径 rs=2Gm/c2r_\mathrm{s}=2Gm/c^2 を使えば amin=2.5rsa_{\min}=2.5\,r_\mathrm{s} で、2 つの地平面(合わせて 2rs2r_\mathrm{s})がちょうど触れ合う直前のスケールになっており、ニュートン的な 2 質点描像が破綻する距離として妥当です。設問7 の m38Mm\approx 38M_\odot を入れると GM/c2=c(GM/c3)1.5kmGM_\odot/c^2 = c\cdot(GM_\odot/c^3)\approx 1.5\,\mathrm{km} より amin5×38×1.5km3×102kma_{\min}\approx 5\times 38\times 1.5\,\mathrm{km} \approx 3\times 10^{2}\,\mathrm{km} です。

第4問 マイケルソン干渉計によるフーリエ分光

Section titled “第4問 マイケルソン干渉計によるフーリエ分光”

光源からの光がビームスプリッタ(BS)で分割され、反射光は鏡 A、透過光は鏡 B へ向かい、それぞれ戻って BS で再び合流します。検出器へ向かう光は「A から戻って BS を透過した成分」と「B から戻って BS で反射された成分」の和、光源へ戻る光はその逆の組み合わせです。複素電場を Eei2πνtEe^{-i2\pi\nu t} と書くので、往復時間 Δt\Delta t の伝播は複素振幅に e+i2πνΔte^{+i2\pi\nu\Delta t} を掛けます。BS の振幅透過・反射係数を p,rp,r、透過による位相差を ϕp\phi_\mathrm{p}、BS と両鏡での反射による位相差をすべて ϕr\phi_\mathrm{r}、鏡の反射率は 100 パーセントとします。後半では r=p=1/2r=p=1/\sqrt2ϕr=π\phi_\mathrm{r}=\piϕp=π/2\phi_\mathrm{p}=\pi/2、BS と鏡 A・B の距離を L+lL+lLlL-l とします。h=6.6×1034Jsh=6.6\times10^{-34}\,\mathrm{J\,s}c=3.0×108m/sc=3.0\times10^{8}\,\mathrm{m/s}kB=1.4×1023J/Kk_\mathrm{B}=1.4\times10^{-23}\,\mathrm{J/K} です。

鏡 A の側へ行く光は、BS で 1 回反射(係数 rr、位相 ϕr\phi_\mathrm{r})し、鏡 A で 1 回反射(係数 1、位相 ϕr\phi_\mathrm{r})して戻ります。往復の伝播で ei2πνΔtAe^{i2\pi\nu\Delta t_\mathrm{A}} が掛かるので

EA=rEiei(2ϕr+2πνΔtA).E_\mathrm{A} = r\,E_\mathrm{i}\,e^{i\left(2\phi_\mathrm{r}+2\pi\nu\Delta t_\mathrm{A}\right)} .

与えられた EB=pEiei(ϕp+ϕr+2πνΔtB)E_\mathrm{B}=pE_\mathrm{i}e^{i(\phi_\mathrm{p}+\phi_\mathrm{r}+2\pi\nu\Delta t_\mathrm{B})} と比べると、prp\leftrightarrow rϕpϕr\phi_\mathrm{p}\leftrightarrow\phi_\mathrm{r} が入れ替わった形になっています。

検出器へ向かう成分は、A からの光が BS を透過(ppϕp\phi_\mathrm{p})したものと、B からの光が BS で反射(rrϕr\phi_\mathrm{r})したものの和です。

Eo=peiϕpEA+reiϕrEB=prEiei(ϕp+2ϕr)[ei2πνΔtA+ei2πνΔtB].E_\mathrm{o} = p\,e^{i\phi_\mathrm{p}}E_\mathrm{A} + r\,e^{i\phi_\mathrm{r}}E_\mathrm{B} = p\,r\,E_\mathrm{i}\,e^{i\left(\phi_\mathrm{p}+2\phi_\mathrm{r}\right)}\left[e^{i2\pi\nu\Delta t_\mathrm{A}} + e^{i2\pi\nu\Delta t_\mathrm{B}}\right].

光源へ戻る成分は逆の組み合わせで、A からの光が BS で反射、B からの光が BS を透過したものの和です。

Es=reiϕrEA+peiϕpEB=Eieiϕr[r2ei(2ϕr+2πνΔtA)+p2ei(2ϕp+2πνΔtB)].E_\mathrm{s} = r\,e^{i\phi_\mathrm{r}}E_\mathrm{A} + p\,e^{i\phi_\mathrm{p}}E_\mathrm{B} = E_\mathrm{i}\,e^{i\phi_\mathrm{r}}\left[r^2 e^{i\left(2\phi_\mathrm{r}+2\pi\nu\Delta t_\mathrm{A}\right)} + p^2 e^{i\left(2\phi_\mathrm{p}+2\pi\nu\Delta t_\mathrm{B}\right)}\right].

EoE_\mathrm{o} では 2 つの経路がどちらも透過 1 回・反射 1 回を経るので係数が prpr で揃い、EsE_\mathrm{s} では r2r^2p2p^2 に分かれるのが違いです。

往復距離から ΔtA=2(L+l)c\Delta t_\mathrm{A} = \dfrac{2(L+l)}{c}ΔtB=2(Ll)c\Delta t_\mathrm{B} = \dfrac{2(L-l)}{c} なので

2πν(ΔtAΔtB)=2πν4lc=8πνlc.2\pi\nu\left(\Delta t_\mathrm{A}-\Delta t_\mathrm{B}\right) = 2\pi\nu\cdot\frac{4l}{c} = \frac{8\pi\nu l}{c}.

EoE_\mathrm{o} では ϕp,ϕr\phi_\mathrm{p},\phi_\mathrm{r} が全体の位相にしか現れないので絶対値をとると落ちます。eiθ1+eiθ22=4cos2θ1θ22\lvert e^{i\theta_1}+e^{i\theta_2}\rvert^2 = 4\cos^2\dfrac{\theta_1-\theta_2}{2}pr=1/2pr=1/2 を使って

Eo2=(pr)2Ei24cos2(4πνlc)=Ei2cos2(4πνlc)=12Ei2[cos(8πνlc)+1],\lvert E_\mathrm{o}\rvert^2 = (pr)^2\lvert E_\mathrm{i}\rvert^2\cdot 4\cos^2\left(\frac{4\pi\nu l}{c}\right) = \lvert E_\mathrm{i}\rvert^2\cos^2\left(\frac{4\pi\nu l}{c}\right) = \frac{1}{2}\lvert E_\mathrm{i}\rvert^2\left[\cos\left(\frac{8\pi\nu l}{c}\right)+1\right],

最後は cos2θ=(1+cos2θ)/2\cos^2\theta = (1+\cos2\theta)/2 を使いました。これが式(1)です。

EsE_\mathrm{s} では r2=p2=1/2r^2=p^2=1/22ϕr=2π2\phi_\mathrm{r}=2\pi2ϕp=π2\phi_\mathrm{p}=\pi なので

Es=12Eieiϕr[ei2πνΔtAei2πνΔtB],Es2=14Ei24sin2(4πνlc)=Ei2sin2(4πνlc).E_\mathrm{s} = \frac{1}{2}E_\mathrm{i}e^{i\phi_\mathrm{r}}\left[e^{i2\pi\nu\Delta t_\mathrm{A}} - e^{i2\pi\nu\Delta t_\mathrm{B}}\right], \qquad \lvert E_\mathrm{s}\rvert^2 = \frac{1}{4}\lvert E_\mathrm{i}\rvert^2\cdot 4\sin^2\left(\frac{4\pi\nu l}{c}\right) = \lvert E_\mathrm{i}\rvert^2\sin^2\left(\frac{4\pi\nu l}{c}\right).

したがって

Eo2+Es2=Ei2[cos2(4πνlc)+sin2(4πνlc)]=Ei2\lvert E_\mathrm{o}\rvert^2 + \lvert E_\mathrm{s}\rvert^2 = \lvert E_\mathrm{i}\rvert^2\left[\cos^2\left(\frac{4\pi\nu l}{c}\right)+\sin^2\left(\frac{4\pi\nu l}{c}\right)\right] = \lvert E_\mathrm{i}\rvert^2

が成り立ちます。検出器側で暗くなる分はちょうど光源側へ戻り、エネルギーは保存しています。この相補性は ϕr\phi_\mathrm{r}ϕp\phi_\mathrm{p} の差が π/2\pi/2 になっていることで保証されます。

I(ν)=I1δ(νν1)+I2δ(νν2)I(\nu)=I_1\delta(\nu-\nu_1)+I_2\delta(\nu-\nu_2) を式(3)に代入するだけです。

Pok=12I1[cos(2πν1kΔlc)+1]+12I2[cos(2πν2kΔlc)+1].P_\mathrm{o}^k = \frac{1}{2}I_1\left[\cos\left(\frac{2\pi\nu_1 k\Delta l}{c}\right)+1\right] + \frac{1}{2}I_2\left[\cos\left(\frac{2\pi\nu_2 k\Delta l}{c}\right)+1\right].

α2πΔlc\alpha \equiv \dfrac{2\pi\Delta l}{c} と置きます。(I1,I2)=(I0,0)(I_1,I_2)=(I_0,0) の場合は Pok=I02[cos(kα(νˉ+Δν))+1]P_\mathrm{o}^k=\dfrac{I_0}{2}\left[\cos\left(k\alpha(\bar\nu+\Delta\nu)\right)+1\right](0,I0)(0,I_0) の場合は νˉΔν\bar\nu-\Delta\nu に置き換わったものです。両者の差は和積の公式で

cos(kα(νˉ+Δν))cos(kα(νˉΔν))=2sin(kανˉ)sin(kαΔν)\cos\left(k\alpha(\bar\nu+\Delta\nu)\right) - \cos\left(k\alpha(\bar\nu-\Delta\nu)\right) = -2\sin\left(k\alpha\bar\nu\right)\sin\left(k\alpha\Delta\nu\right)

となります。これが任意の Δν\Delta\nu とすべての k=0,1,2,k=0,1,2,\dots でゼロになる条件は sin(ανˉ)=0\sin(\alpha\bar\nu)=0、すなわち

2πΔlcνˉ=πN(N=1,2,)νˉ=Nc2Δl\frac{2\pi\Delta l}{c}\bar\nu = \pi N\quad (N=1,2,\dots) \qquad\Longleftrightarrow\qquad \bar\nu = \frac{Nc}{2\Delta l}

です(このとき sin(kανˉ)=sin(kπN)=0\sin(k\alpha\bar\nu)=\sin(k\pi N)=0 ですべての kk で消えます)。最小のものが求めるナイキスト周波数で

νn=c2Δl.\nu_\mathrm{n} = \frac{c}{2\Delta l}.

式(3)の位相 2πνkΔl/c2\pi\nu k\Delta l/c は、実効的なサンプリング間隔 Δl/c\Delta l/c で時間信号を離散化したのと同じ形なので、その半分の逆数がナイキスト周波数になる、という標準的な結果と一致します。鏡の移動幅が Δl/4\Delta l/4 で光路差の変化が Δl\Delta l になることに注意します。

xhν/(kBT)x\equiv h\nu/(k_\mathrm{B}T) とします。hνkBTh\nu\ll k_\mathrm{B}T では ex1xe^{x}-1\simeq x なので

I(ν;T)2hν3c2kBThν=2ν2kBTc2I(\nu;T) \simeq \frac{2h\nu^3}{c^2}\cdot\frac{k_\mathrm{B}T}{h\nu} = \frac{2\nu^2 k_\mathrm{B}T}{c^2}

がレイリー・ジーンズの放射法則です。hνkBTh\nu\gg k_\mathrm{B}T では ex1exe^{x}-1\simeq e^{x} なので

I(ν;T)2hν3c2exp(hνkBT)I(\nu;T) \simeq \frac{2h\nu^3}{c^2}\exp\left(-\frac{h\nu}{k_\mathrm{B}T}\right)

がウィーンの放射法則です。前者に hh が残らないのが古典極限であることの印で、後者の指数関数的減衰が高周波側を強く抑えます。

I(ν;T)I(\nu;T)x=hν/(kBT)x=h\nu/(k_\mathrm{B}T) の関数 f(x)=x3/(ex1)f(x)=x^3/(e^x-1) に定数を掛けた形なので、最大は x2.8x\approx 2.8 で起こります。

νmax=2.8kBTh=2.8×1.4×1023×2.76.6×1034=1.06×10226.6×1034=1.6×1011Hz\nu_{\max} = \frac{2.8\,k_\mathrm{B}T}{h} = \frac{2.8\times 1.4\times 10^{-23}\times 2.7}{6.6\times 10^{-34}} = \frac{1.06\times 10^{-22}}{6.6\times 10^{-34}} = 1.6\times 10^{11}\,\mathrm{Hz}

すなわち約 160GHz160\,\mathrm{GHz}(波長 c/νmax1.9mmc/\nu_{\max}\approx 1.9\,\mathrm{mm})です。次に νn=10νmax=1.6×1012Hz\nu_\mathrm{n}=10\nu_{\max}=1.6\times 10^{12}\,\mathrm{Hz} となるように Δl\Delta l を決めます。設問5 の νn=c/(2Δl)\nu_\mathrm{n}=c/(2\Delta l) より

Δl=c2νn=3.0×1082×1.6×1012=9.4×105m\Delta l = \frac{c}{2\nu_\mathrm{n}} = \frac{3.0\times 10^{8}}{2\times 1.6\times 10^{12}} = 9.4\times 10^{-5}\,\mathrm{m}

です。答えは Δl9.4×105m\Delta l \approx 9.4\times 10^{-5}\,\mathrm{m}(約 94μm94\,\mu\mathrm{m})で、鏡の移動ステップ Δl/4\Delta l/4 は約 23μm23\,\mu\mathrm{m} になります。Δl=λmax/20\Delta l = \lambda_{\max}/20 という関係で覚えておけば桁の確認が速いです。

中性子と原子核の反応断面積、中性 π\pi 中間子の 2 光子崩壊からの質量再構成、制動放射の次元解析、磁場中の運動量測定という 4 つの独立な話題を扱います。真空中の光速を cc とします。

中性子は 940MeV/c2940\,\mathrm{MeV}/c^2 で e、陽子は 938MeV/c2938\,\mathrm{MeV}/c^2 で d、中性 π\pi 中間子は 135MeV/c2135\,\mathrm{MeV}/c^2 で b です。残りは a が電子(0.5110.511)、c が KK 中間子(494494)、f が Δ\Delta 粒子(12321232)に対応します。中性子が陽子よりわずかに重いことが、中性子のベータ崩壊が起こり陽子が安定である理由です。

np+e+νˉen \to p + e^- + \bar\nu_e

生じる粒子は陽子、電子、反電子ニュートリノの 3 つです。

ΔI=nσIΔx\Delta I = -n\sigma I\Delta x を微分方程式にすると dIdx=nσI\dfrac{dI}{dx} = -n\sigma I で、n,σn,\sigma が一様なら

I(x)=I0enσx,I1=I(d)=I0enσd.I(x) = I_0 e^{-n\sigma x},\qquad I_1 = I(d) = I_0 e^{-n\sigma d}.

これを σ\sigma について解いて

σ=1ndlnI0I1.\sigma = \frac{1}{nd}\ln\frac{I_0}{I_1}.

[n]=m3[n]=\mathrm{m^{-3}}[d]=m[d]=\mathrm{m} なので σ\sigma は面積の次元をもち、断面積として正しい次元です。

π0γγ\pi^0\to\gamma\gamma の 4 元運動量保存から、不変質量を計算します。光子は pi=Ei/c\lvert \vec p_i\rvert = E_i/c なので

m2c4=(E1+E2)2p1+p22c2=(E1+E2)2(E12+E22+2E1E2cosϕ)=2E1E2(1cosϕ).m^2c^4 = \left(E_1+E_2\right)^2 - \left\lvert \vec p_1+\vec p_2\right\rvert^2c^2 = \left(E_1+E_2\right)^2 - \left(E_1^2+E_2^2+2E_1E_2\cos\phi\right) = 2E_1E_2\left(1-\cos\phi\right).

1cosϕ=2sin2(ϕ/2)1-\cos\phi = 2\sin^2(\phi/2) を使って

m=2E1E2(1cosϕ)c2=2E1E2c2sinϕ2.m = \frac{\sqrt{2E_1E_2\left(1-\cos\phi\right)}}{c^2} = \frac{2\sqrt{E_1E_2}}{c^2}\sin\frac{\phi}{2}.

ϕ0\phi\to 0m0m\to 0(同方向に飛ぶ 2 光子の不変質量はゼロ)、ϕ=π\phi=\pi で静止系の対称崩壊 m=2E/c2m=2E/c^2 になり、いずれも妥当です。

MKSA 単位系で、質量 M\mathsf{M}、長さ L\mathsf{L}、時間 T\mathsf{T}、電流 I\mathsf{I} を使って各量の次元を書きます。

[P]=ML2T3,[a]=LT2,[e]=IT,[ε0]=M1L3T4I2,[c]=LT1[P] = \mathsf{M}\mathsf{L}^2\mathsf{T}^{-3},\quad [a] = \mathsf{L}\mathsf{T}^{-2},\quad [e] = \mathsf{I}\mathsf{T},\quad [\varepsilon_0] = \mathsf{M}^{-1}\mathsf{L}^{-3}\mathsf{T}^{4}\mathsf{I}^{2},\quad [c] = \mathsf{L}\mathsf{T}^{-1}

[ε0][\varepsilon_0] はクーロンの法則 F=q1q2/(4πε0r2)F=q_1q_2/(4\pi\varepsilon_0r^2) から出しました)。P=Aafegε0hciP=Aa^fe^g\varepsilon_0^hc^i の両辺で次元を合わせると

M:1=hI:0=g+2hL:2=f3h+iT:3=2f+g+4hi\begin{aligned} \mathsf{M}:&\quad 1 = -h\\ \mathsf{I}:&\quad 0 = g+2h\\ \mathsf{L}:&\quad 2 = f-3h+i\\ \mathsf{T}:&\quad -3 = -2f+g+4h-i \end{aligned}

第1式から h=1h=-1、第2式から g=2g=2。第3式は f+i=1f+i=-1、第4式は 2fi=1-2f-i=-1 となり、両者から f=2f=2i=3i=-3 です。したがって

f=2,g=2,h=1,i=3,P=Ae2a2ε0c3.f=2,\quad g=2,\quad h=-1,\quad i=-3,\qquad P = A\,\frac{e^2a^2}{\varepsilon_0c^3}.

これは A=1/(6π)A=1/(6\pi) とすればラーモアの公式そのものです。

アルミニウムと鉛の比較では、鉛のほうが PP が大きくなります。理由は次のとおりです。Pa2P\propto a^2 で、aa は原子核のクーロン場による加速度なので、電荷 ZeZe の原子核から距離 rr の点では a=Ze24πε0mer2a = \dfrac{Ze^2}{4\pi\varepsilon_0 m_e r^2} となり ZZ に比例します。したがって同じ距離を通る電子に対して PZ2P\propto Z^2 です。原子番号は ZAl=13Z_{\mathrm{Al}}=13ZPb=82Z_{\mathrm{Pb}}=82 なので、鉛では (82/13)240(82/13)^2\approx 40 倍ほど放射が強くなります。電磁カロリメータの吸収体に鉛のような高 ZZ 物質を使うのは、この Z2Z^2 依存によって放射長が短くなるからです。

磁場中では磁気力が仕事をしないので v\lvert\vec v\rvertp\lvert\vec p\rvert は一定で、運動は磁場に垂直な面内の円になります。運動方程式は相対論的にも

dpdt=qv×B,dpdt=qvB\frac{d\vec p}{dt} = q\,\vec v\times\vec B,\qquad \left\lvert\frac{d\vec p}{dt}\right\rvert = qvB

で、半径 RR の円運動では運動量ベクトルが角速度 v/Rv/R で向きを変えるだけなので dp/dt=pv/R\left\lvert d\vec p/dt\right\rvert = p\,v/R です。両者を等置して vv を消せば

pvR=qvBp=qBR.\frac{pv}{R} = qvB \qquad\Longrightarrow\qquad p = qBR .

つぎに幾何です。軌道の接線が xx 軸(入射面・出射面の法線)となす角を θ\theta、弧長を ss とすると、dxds=cosθ\dfrac{dx}{ds}=\cos\theta と、曲率一定であることから dθds=1R\dfrac{d\theta}{ds}=\dfrac{1}{R} が成り立ちます。図の配置では BB が紙面の裏向き、q>0q>0+x+x 方向に入射するので力は qv×Bq\vec v\times\vec B により図の上向きで、θ\thetass とともに増加します。この 2 式から dx=Rcosθdθdx = R\cos\theta\,d\theta で、入射面から出射面まで積分すると

L=θ1θ2Rcosθdθ=R(sinθ2sinθ1)R=Lsinθ2sinθ1.L = \int_{\theta_1}^{\theta_2}R\cos\theta\,d\theta = R\left(\sin\theta_2-\sin\theta_1\right) \qquad\Longrightarrow\qquad R = \frac{L}{\sin\theta_2-\sin\theta_1}.

磁場領域の xx 方向の厚みだけで RR が決まり、粒子が磁場領域のどの高さを通るかに依らないことがこの手法の要点です。

κ=sinθ2sinθ1Lθ2θ1L\kappa = \dfrac{\sin\theta_2-\sin\theta_1}{L} \simeq \dfrac{\theta_2-\theta_1}{L} で、θ1\theta_1θ2\theta_2 の誤差に相関がなく LL に誤差がないので、誤差の伝播は偏微分の二乗和になります。

(Δκ)2=(κθ1)2(Δθ1)2+(κθ2)2(Δθ2)2=(Δθ1)2+(Δθ2)2L2\left(\Delta\kappa\right)^2 = \left(\frac{\partial\kappa}{\partial\theta_1}\right)^2\left(\Delta\theta_1\right)^2 + \left(\frac{\partial\kappa}{\partial\theta_2}\right)^2\left(\Delta\theta_2\right)^2 = \frac{\left(\Delta\theta_1\right)^2+\left(\Delta\theta_2\right)^2}{L^2}

したがって

Δκ=(Δθ1)2+(Δθ2)2L.\Delta\kappa = \frac{\sqrt{\left(\Delta\theta_1\right)^2+\left(\Delta\theta_2\right)^2}}{L}.

p=qB/κp=qB/\kappa なので相対誤差は Δpp=Δκκ=(Δθ1)2+(Δθ2)2θ2θ1\dfrac{\Delta p}{p} = \dfrac{\Delta\kappa}{\kappa} = \dfrac{\sqrt{(\Delta\theta_1)^2+(\Delta\theta_2)^2}}{\theta_2-\theta_1} となり、高運動量(曲がりが小さい)ほど測定が苦しくなることが読めます。

同軸ケーブルの中心軸を xx 軸に取り、単位長さあたりの電気容量を CC、自己インダクタンスを LL とすると、内部導体の電流 I(x,t)I(x,t) と両導体間の電圧 V(x,t)V(x,t) は伝送線路方程式

Ix=CVt,Vx=LIt\frac{\partial I}{\partial x} = -C\frac{\partial V}{\partial t},\qquad \frac{\partial V}{\partial x} = -L\frac{\partial I}{\partial t}

に従います。内部導体半径 aa、外部導体内半径 bb、絶縁体の誘電率 ε\varepsilon、透磁率は真空と同じ μ0\mu_0 として C=2πεln(b/a)C=\dfrac{2\pi\varepsilon}{\ln(b/a)}L=μ0ln(b/a)2πL=\dfrac{\mu_0\ln(b/a)}{2\pi} が与えられています。c=3.0×108m/sc=3.0\times10^{8}\,\mathrm{m/s}μ0=4π×107H/m\mu_0=4\pi\times10^{-7}\,\mathrm{H/m} を使います。

第2式を xx で微分し、第1式を代入します。

2Vx2=LtIx=Lt(CVt)=LC2Vt2\frac{\partial^2V}{\partial x^2} = -L\frac{\partial}{\partial t}\frac{\partial I}{\partial x} = -L\frac{\partial}{\partial t}\left(-C\frac{\partial V}{\partial t}\right) = LC\frac{\partial^2V}{\partial t^2}

すなわち 2Vx2=LC2Vt2\dfrac{\partial^2V}{\partial x^2} = LC\dfrac{\partial^2 V}{\partial t^2} で、II についても順序を入れ替えて同じ式が出ます。これを 2Vx2=1v22Vt2\dfrac{\partial^2 V}{\partial x^2}=\dfrac{1}{v^2}\dfrac{\partial^2V}{\partial t^2} と比べて

v=1LC.v = \frac{1}{\sqrt{LC}} .

与えられた CCLL の積では ln(b/a)\ln(b/a) が消えて

LC=μ0ln(b/a)2π2πεln(b/a)=μ0ε,v=1μ0ε=cεrLC = \frac{\mu_0\ln(b/a)}{2\pi}\cdot\frac{2\pi\varepsilon}{\ln(b/a)} = \mu_0\varepsilon, \qquad v = \frac{1}{\sqrt{\mu_0\varepsilon}} = \frac{c}{\sqrt{\varepsilon_\mathrm{r}}}

となります(ε=εrε0\varepsilon=\varepsilon_\mathrm{r}\varepsilon_0c=1/μ0ε0c=1/\sqrt{\mu_0\varepsilon_0})。速度は導体の寸法によらず絶縁体の比誘電率だけで決まるので、この設問では ln(b/a)=1.8\ln(b/a)=1.8 は使いません。εr=2.2\varepsilon_\mathrm{r}=2.2 から 2.2=1.48\sqrt{2.2}=1.48

v=3.0×1081.48=2.0×108m/s,λ=vf=2.0×1085.0×1080.4m.v = \frac{3.0\times 10^{8}}{1.48} = 2.0\times 10^{8}\,\mathrm{m/s}, \qquad \lambda = \frac{v}{f} = \frac{2.0\times 10^{8}}{5.0\times 10^{8}} \approx 0.4\,\mathrm{m}.

答えは λ0.4m\lambda \approx 0.4\,\mathrm{m} です。ケーブルの長さがこの程度になると分布定数回路として扱う必要が出てきます。

xx 軸正方向に進む波は V=V(xvt)V=V(x-vt)I=I(xvt)I=I(x-vt) と書けます。これを Vx=LIt\dfrac{\partial V}{\partial x}=-L\dfrac{\partial I}{\partial t} に入れると V=LvIV' = Lv\,I'、積分定数をゼロとして V=LvIV = Lv\,I、つまり

Z0=Lv=LLC=LC.Z_0 = Lv = \frac{L}{\sqrt{LC}} = \sqrt{\frac{L}{C}} .

与えられた C,LC,L を代入すると

Z0=μ0ln(b/a)/2π2πε/ln(b/a)=ln(b/a)2πμ0ε=ln(b/a)2πμ0cεr.Z_0 = \sqrt{\frac{\mu_0\ln(b/a)/2\pi}{2\pi\varepsilon/\ln(b/a)}} = \frac{\ln(b/a)}{2\pi}\sqrt{\frac{\mu_0}{\varepsilon}} = \frac{\ln(b/a)}{2\pi}\cdot\frac{\mu_0c}{\sqrt{\varepsilon_\mathrm{r}}} .

μ0c=4π×107×3.0×108=120π3.8×102Ω\mu_0c = 4\pi\times10^{-7}\times 3.0\times10^{8} = 120\pi \approx 3.8\times10^{2}\,\Omega なので

Z0=1.82π120π1.48=1.8×1202×1.487×101ΩZ_0 = \frac{1.8}{2\pi}\cdot\frac{120\pi}{1.48} = \frac{1.8\times 120}{2\times 1.48} \approx 7\times 10^{1}\,\Omega

(詳しくは 73Ω73\,\Omega)です。テレビ用同軸ケーブルの規格値 75Ω75\,\Omega とよく合っています。

入射パルスと反射パルスの重ね合わせが終端で抵抗器のオームの法則を満たす、という条件で決まります。+x+x 方向へ進む波では V=Z0IV=Z_0Ix-x 方向へ進む波では設問3 と同じ計算で V=Z0IV=-Z_0I になるので、終端における全電圧と全電流は

V1+V2,V1Z0V2Z0V_1+V_2,\qquad \frac{V_1}{Z_0}-\frac{V_2}{Z_0}

です。これらが V=rIV=rI を満たすので

V1+V2=rV1V2Z0V2=rZ0r+Z0V1.V_1+V_2 = r\,\frac{V_1-V_2}{Z_0} \qquad\Longrightarrow\qquad V_2 = \frac{r-Z_0}{r+Z_0}V_1 .

r=Z0r=Z_0V2=0V_2=0(無反射)、開放端 rr\to\inftyV2=V1V_2=V_1(同符号で全反射)、短絡端 r=0r=0V2=V1V_2=-V_1(符号反転して全反射)となり、いずれも既知の振る舞いです。

V,Iei(kxωt)V,I \propto e^{i(kx-\omega t)}kk は複素数)と置いて、Vx=LItRI\dfrac{\partial V}{\partial x}=-L\dfrac{\partial I}{\partial t}-RIIx=CVt\dfrac{\partial I}{\partial x}=-C\dfrac{\partial V}{\partial t} に代入します。

ikV~=(iωLR)I~,ikI~=iωCV~ik\tilde V = \left(i\omega L - R\right)\tilde I,\qquad ik\tilde I = i\omega C\tilde V

から I~\tilde I を消すと

k2=(iωLR)ωCi=ω2LC(1+iRωL).k^2 = \frac{\left(i\omega L-R\right)\omega C}{i} = \omega^2LC\left(1+\frac{iR}{\omega L}\right).

RLωR\ll L\omega として平方根を展開すると

kωLC(1+iR2ωL)=ωv+i2RCL=ωv+iR2Z0.k \simeq \omega\sqrt{LC}\left(1+\frac{iR}{2\omega L}\right) = \frac{\omega}{v} + \frac{i}{2}R\sqrt{\frac{C}{L}} = \frac{\omega}{v}+\frac{iR}{2Z_0}.

虚部が正なので eikxe^{ikx}+x+x 方向へ進みながら減衰します。右端は Z0Z_0 で終端されていて反射がないため、ケーブル内は純粋な進行波で、長さ dd の右端での振幅は

V2=V1exp(Rd2Z0)=V1exp(Rd2CL).V_2 = V_1\exp\left(-\frac{Rd}{2Z_0}\right) = V_1\exp\left(-\frac{Rd}{2}\sqrt{\frac{C}{L}}\right).

この近似では減衰量が ω\omega に依らないことに注意します。周波数依存性は、次の設問で見る表皮効果によって RRω\omega に依存することから生じます。

軸方向の電流密度 j=J(r)eiωtj=J(r)e^{-i\omega t}、円周方向の磁束密度 BφB_\varphi、軸方向の誘導電場 EE について、良導体で変位電流を無視したアンペールの法則とファラデーの法則は

1r(rBφ)r=μj,Er=Bφt\frac{1}{r}\frac{\partial\left(rB_\varphi\right)}{\partial r} = \mu j, \qquad \frac{\partial E}{\partial r} = \frac{\partial B_\varphi}{\partial t}

です(後者は ×E=B/t\nabla\times\vec E = -\partial\vec B/\partial tφ\varphi 成分)。j=σEj=\sigma E と時間依存性 eiωte^{-i\omega t} を入れ、Bφ=B(r)eiωtB_\varphi = \mathcal B(r)e^{-i\omega t} とすると

1σdJdr=iωBiωσ1rddr(rdJdr)=μJ.\frac{1}{\sigma}\frac{dJ}{dr} = -i\omega\mathcal B \qquad\Longrightarrow\qquad \frac{i}{\omega\sigma}\frac{1}{r}\frac{d}{dr}\left(r\frac{dJ}{dr}\right) = \mu J .

δa\delta\ll a、すなわち JJ が表面近傍で厚さ δ\delta の範囲で変化するとき、曲率の効果 1rdJdr\dfrac{1}{r}\dfrac{dJ}{dr}d2Jdr2\dfrac{d^2J}{dr^2} に対して δ/a\delta/a の比で小さいので落とせて

d2Jdr2=iμωσJ.\frac{d^2J}{dr^2} = -i\mu\omega\sigma\,J .

JeκrJ\propto e^{\kappa r} と置くと κ2=iμωσ\kappa^2 = -i\mu\omega\sigma、つまり κ=eiπ/4μωσ=1i2μωσ\kappa = e^{-i\pi/4}\sqrt{\mu\omega\sigma} = \dfrac{1-i}{\sqrt2}\sqrt{\mu\omega\sigma}(内部へ向かって減衰する、すなわち Reκ>0\mathrm{Re}\,\kappa>0 の根を選びました)。よって

J(r)=J(a)exp[(Reκ)(ra)],1δ=Reκ=μωσ2\lvert J(r)\rvert = \lvert J(a)\rvert \exp\left[\left(\mathrm{Re}\,\kappa\right)\left(r-a\right)\right], \qquad \frac{1}{\delta} = \mathrm{Re}\,\kappa = \sqrt{\frac{\mu\omega\sigma}{2}}

すなわち

δ=2μωσ.\delta = \sqrt{\frac{2}{\mu\omega\sigma}} .

ω\omega が大きいほど、また導電率が高いほど電流は表面に集まります。実効断面積が 2πaδ2\pi a\delta に減るため単位長さの抵抗は R1/δωR\propto 1/\delta \propto \sqrt{\omega} で増え、設問5 の減衰は高周波で exp(constωd)\exp\left(-\mathrm{const}\cdot\sqrt{\omega}\,d\right) の形で強くなります。次元は [μωσ]=Hm1s1Ω1m1=m2[\mu\omega\sigma] = \mathrm{H\,m^{-1}\cdot s^{-1}\cdot \Omega^{-1}m^{-1}} = \mathrm{m^{-2}} なので δ\delta は長さで正しくなっています。

出典: 東京大学大学院理学系研究科 物理学専攻 平成31年度 修士課程 入学試験問題 物理学。問題文は要約して引用しています。

Report an error in this article ・Operated by: Mugen Giken LLCPricingTermsLegal notice

© 2026 夢現技研合同会社 ・Feeding the text to an LLM is welcome. Code samples are MIT licensed.