院試hub

東京大学 院試 過去問 解答例

東大 理学系研究科 物理学専攻 専門科目(物理学) 2012年度 院試 過去問 解答例・解説(全6問)

全6問。電磁気学・回路2問・熱・統計力学1問・物性物理1問。テーマタグは6件(静電場・分配関数・ハミルトン形式)。2011年度と共通のテーマは分配関数・ハミルトン形式。

最終更新:

このページで公開
解説6問・解答6問・最終答6問(全6問)
公開範囲
2012年度は解答・最終答まで全文公開
問題本文
非収録

東大 専門科目(物理学) 2012年度 院試 過去問の出題内容(全6問)

この6問の分野は電磁気学・回路2問・熱・統計力学1問・物性物理1問・量子力学1問です。

大問分野主題解説の小見出し最終答
第1問量子力学スピンと回転磁場磁場の向きを取り違えない / 相互作用表示の役割あり
第2問熱・統計力学剛体回転子とオルト・パラ水素古典分配関数の / オルト・パラは熱平衡と組成凍結を分けるあり
第3問電磁気学・回路ローレンツ振動子と屈折率分散分極ではなく電流密度から入る / 共鳴の位置が色順を決めるあり
第4問しきい値・飛行時間・寿命測定しきい値は不変量で処理する / 指数分布の寿命推定あり
第5問物性物理中性子散乱と拡散係数横軸の単位に注意 / グラフ読み取りは傾きで安定化するあり
第6問電磁気学・回路同軸ケーブルと整合減衰器単位長さあたりの量で統一する / T型減衰器は二条件で決まるあり

2012年度の出題テーマと、同じテーマを出した他大学・他年度

この年度は5問に6テーマが出ています。

前年度(2011年度)との違い

大問数
2011年度 6問 → 2012年度 6問
両年度に出たテーマ
分配関数・ハミルトン形式
2012年度で新しく出たテーマ
静電場・最尤推定・マクスウェル方程式・角運動量
2011年度のページを見る

第1問 — 量子力学:スピンと回転磁場

磁場の向きを取り違えない

この問題の一定磁場は xx 方向である。zz 方向磁場として解くと固有状態も期待値も別物になる。まずハミルトニアンの行列を磁場の向きから作るのが安全である。

相互作用表示の役割

H1=−aSz\mathcal H_1=-aS_z による速い zz 軸回転を外すと,横磁場は位相因子 e±i(ω−a)te^{\pm i(\omega-a)t} だけを持つ。共鳴でこの因子が消えることが,ラビ振動の計算を一気に簡単にする。

解答

角運動量とスピンの基本代数から順に整理する。以下,設問の後半では ℏ=1\hbar=1 とする。

  1. 軌道角運動量を Lx=ypz−zpy,Ly=zpx−xpz,Lz=xpy−ypx L_x=yp_z-zp_y,\quad L_y=zp_x-xp_z,\quad L_z=xp_y-yp_x とおく。例えば [Lx,Ly]=[ypz−zpy,  zpx−xpz]=[ypz,zpx]−[ypz,xpz]−[zpy,zpx]+[zpy,xpz]. \begin{aligned} [L_x,L_y] &=[yp_z-zp_y,\;zp_x-xp_z] \\ &= [yp_z,zp_x]-[yp_z,xp_z]-[zp_y,zp_x]+[zp_y,xp_z]. \end{aligned} 異なる成分の座標・運動量は可換なので,残る項だけを計算すると [ypz,zpx]=y[pz,z]px=−iℏypx,[zpy,xpz]=x[z,pz]py=iℏxpy. [yp_z,zp_x]=y[p_z,z]p_x=-i\hbar yp_x, \qquad [zp_y,xp_z]=x[z,p_z]p_y=i\hbar xp_y. したがって [Lx,Ly]=iℏ(xpy−ypx)=iℏLz. [L_x,L_y]=i\hbar(xp_y-yp_x)=i\hbar L_z. 巡回置換により [Li,Lj]=iℏϵijkLk [L_i,L_j]=i\hbar\epsilon_{ijk}L_k が得られる。
  2. スピン 1/21/2 の行列は同じ交換関係を満たす必要がある。与えられた Sx,SzS_x,S_z に対し Sy=ℏ2(0−ii0) S_y=\frac{\hbar}{2} \begin{pmatrix} 0&-i\\ i&0 \end{pmatrix} とすれば,[Sx,Sy]=iℏSz[S_x,S_y]=i\hbar S_z などが成り立つ。
  3. 一定磁場が xx 方向,H=(H,0,0)\boldsymbol H=(H,0,0) であるから H=−μHSx=−aSx=−a2(0110),a=μH. \mathcal H=-\mu H S_x=-aS_x =-\frac{a}{2} \begin{pmatrix} 0&1\\ 1&0 \end{pmatrix}, \qquad a=\mu H. SxS_x の固有状態を用いれば直ちに E1=−a2,∣1⟩=12(11),E2=a2,∣2⟩=12(1−1) E_1=-\frac{a}{2},\quad |1\rangle=\frac{1}{\sqrt2} \begin{pmatrix}1\\1\end{pmatrix}, \qquad E_2=\frac{a}{2},\quad |2\rangle=\frac{1}{\sqrt2} \begin{pmatrix}1\\-1\end{pmatrix} である。エネルギーは磁場方向のスピン成分で決まる。
  4. 初期状態を SzS_z 上向き ∣ϕ(0)⟩=(10) |\phi(0)\rangle= \begin{pmatrix}1\\0\end{pmatrix} とする。時間発展は ∣ϕ(t)⟩=e−iHt∣ϕ(0)⟩=eiaSxt∣ϕ(0)⟩=(cos⁡(at/2)isin⁡(at/2)). |\phi(t)\rangle =e^{-i\mathcal Ht}|\phi(0)\rangle =e^{iaS_xt}|\phi(0)\rangle = \begin{pmatrix} \cos(at/2)\\ i\sin(at/2) \end{pmatrix}. よって期待値は ⟨Sx⟩=0,⟨Sy⟩=12sin⁡at,⟨Sz⟩=12cos⁡at. \langle S_x\rangle=0,\qquad \langle S_y\rangle=\frac12\sin at,\qquad \langle S_z\rangle=\frac12\cos at. 磁場が xx 軸まわりの回転を作るので,SzS_z 成分と SyS_y 成分が入れ替わる。
  5. 回転磁場を H=(H0cos⁡ωt, −H0sin⁡ωt, H) \boldsymbol H=(H_0\cos\omega t,\,-H_0\sin\omega t,\,H) と書くと H=H1+H2,H1=−aSz,H2=−a0(Sxcos⁡ωt−Sysin⁡ωt) \mathcal H=\mathcal H_1+\mathcal H_2,\qquad \mathcal H_1=-aS_z,\qquad \mathcal H_2=-a_0(S_x\cos\omega t-S_y\sin\omega t) である。ただし a0=μH0a_0=\mu H_0。相互作用表示 ∣ψ(t)⟩=eiH1t∣ϕ(t)⟩ |\psi(t)\rangle=e^{i\mathcal H_1t}|\phi(t)\rangle では iddt∣ψ(t)⟩=eiH1tH2e−iH1t∣ψ(t)⟩ i\frac{d}{dt}|\psi(t)\rangle = e^{i\mathcal H_1t}\mathcal H_2e^{-i\mathcal H_1t} |\psi(t)\rangle となる。行列表現で計算すると eiH1tH2e−iH1t=−a02(0ei(ω−a)te−i(ω−a)t0). e^{i\mathcal H_1t}\mathcal H_2e^{-i\mathcal H_1t} = -\frac{a_0}{2} \begin{pmatrix} 0&e^{i(\omega-a)t}\\ e^{-i(\omega-a)t}&0 \end{pmatrix}. したがって,回転磁場の角周波数とゼーマン分裂 aa のずれだけが相互作用表示の時間依存性として残る。
  6. 共鳴条件 ω=a\omega=a では相互作用表示のハミルトニアンが HI=−a0Sx \mathcal H_I=-a_0S_x という時間によらない形になる。初期状態を再び SzS_z 上向きに取ると ∣ψ(t)⟩=eia0Sxt(10)=(cos⁡(a0t/2)isin⁡(a0t/2)). |\psi(t)\rangle=e^{ia_0S_xt} \begin{pmatrix}1\\0\end{pmatrix} = \begin{pmatrix} \cos(a_0t/2)\\ i\sin(a_0t/2) \end{pmatrix}. SzS_z は zz 軸まわりの戻し変換と可換なので,シュレーディンガー表示でも ⟨Sz⟩=12cos⁡(a0t) \langle S_z\rangle=\frac12\cos(a_0t) である。共鳴では弱い横磁場の大きさ H0H_0 が振動周期を決める。

最終答

% Sy=ℏ2(0−ii0). S_y=\frac{\hbar}{2}\begin{pmatrix}0&-i\\ i&0\end{pmatrix}. H=(H,0,0)\boldsymbol H=(H,0,0) では H=−aSx,E=∓a2,∣E=−a/2⟩=12(11). \mathcal H=-aS_x,\quad E=\mp \frac{a}{2},\quad |E=-a/2\rangle=\frac{1}{\sqrt2}\binom{1}{1}. 初期 SzS_z 上向きなら ⟨Sx⟩=0,⟨Sy⟩=12sin⁡at,⟨Sz⟩=12cos⁡at. \langle S_x\rangle=0,\quad \langle S_y\rangle=\frac12\sin at,\quad \langle S_z\rangle=\frac12\cos at. 回転磁場の相互作用表示では HI=−a02(0ei(ω−a)te−i(ω−a)t0), \mathcal H_I=-\frac{a_0}{2} \begin{pmatrix}0&e^{i(\omega-a)t}\\ e^{-i(\omega-a)t}&0\end{pmatrix}, 共鳴 ω=a\omega=a で ⟨Sz⟩=12cos⁡(a0t)\langle S_z\rangle=\frac12\cos(a_0t)。

第2問 — 統計力学:剛体回転子とオルト・パラ水素

古典分配関数の sin⁡θ\sin\theta

pϕp_\phi のガウス積分から sin⁡θ\sin\theta が出る。これが球面の面積要素と合わさり,最終的に 4π4\pi が現れる。ここを落とすと係数だけでなく温度依存の理解も曖昧になる。

オルト・パラは熱平衡と組成凍結を分ける

急冷後の 3:13:1 は熱平衡比ではなく組成が凍結した比である。各成分の内部では低温の最低回転準位に落ちる,という二段階で考えると平均エネルギーを誤らない。

解答

線形二原子分子の回転自由度は二つである。古典極限,量子準位,核スピン統計を分けて扱う。

  1. 古典的な剛体回転子のエネルギーは E=12I(pθ2+pϕ2sin⁡2θ) E=\frac{1}{2I}\left(p_\theta^2+\frac{p_\phi^2}{\sin^2\theta}\right) である。一分子の回転分配関数を Zcl=1h2∫0πdθ∫02πdϕ∫−∞∞dpθ∫−∞∞dpϕ e−βE Z_{\mathrm{cl}} =\frac{1}{h^2}\int_0^\pi d\theta\int_0^{2\pi}d\phi \int_{-\infty}^{\infty}dp_\theta \int_{-\infty}^{\infty}dp_\phi\, e^{-\beta E} と書く。運動量積分は ∫dpθe−βpθ2/(2I)=2πIβ,∫dpϕe−βpϕ2/(2Isin⁡2θ)=sin⁡θ2πIβ. \int dp_\theta e^{-\beta p_\theta^2/(2I)} =\sqrt{\frac{2\pi I}{\beta}}, \quad \int dp_\phi e^{-\beta p_\phi^2/(2I\sin^2\theta)} =\sin\theta\sqrt{\frac{2\pi I}{\beta}}. よって Zcl=1h22πIβ∫02πdϕ∫0πsin⁡θ dθ=8π2IkBTh2. Z_{\mathrm{cl}} =\frac{1}{h^2}\frac{2\pi I}{\beta} \int_0^{2\pi}d\phi\int_0^\pi\sin\theta\,d\theta =\frac{8\pi^2Ik_{\mathrm{B}}T}{h^2}. Θrot=ℏ22IkB \Theta_{\mathrm{rot}}=\frac{\hbar^2}{2Ik_{\mathrm{B}}} を用いれば Zcl=TΘrot. Z_{\mathrm{cl}}=\frac{T}{\Theta_{\mathrm{rot}}}. したがって U=−∂∂βlog⁡Zcl=kBT,C=∂U∂T=kB. U=-\frac{\partial}{\partial\beta}\log Z_{\mathrm{cl}} =k_{\mathrm{B}}T, \qquad C=\frac{\partial U}{\partial T}=k_{\mathrm{B}}.
  2. 量子力学では回転準位は El=ℏ22Il(l+1)=kBΘrotl(l+1),gl=2l+1 E_l=\frac{\hbar^2}{2I}l(l+1) =k_{\mathrm{B}}\Theta_{\mathrm{rot}}l(l+1), \qquad g_l=2l+1 である。よって Zq=∑l=0∞(2l+1)exp⁡[−ΘrotTl(l+1)]. Z_{\mathrm{q}}=\sum_{l=0}^{\infty}(2l+1) \exp\left[-\frac{\Theta_{\mathrm{rot}}}{T}l(l+1)\right]. 高温 T≫ΘrotT\gg\Theta_{\mathrm{rot}} では和を積分で近似でき, Zq≃∫0∞(2l+1)e−Θrotl(l+1)/T dl=TΘrot. Z_{\mathrm{q}}\simeq\int_0^\infty (2l+1) e^{-\Theta_{\mathrm{rot}}l(l+1)/T}\,dl =\frac{T}{\Theta_{\mathrm{rot}}}. したがって古典結果 U=kBT, C=kBU=k_{\mathrm{B}}T,\ C=k_{\mathrm{B}} に一致する。
  3. 低温 T≪ΘrotT\ll\Theta_{\mathrm{rot}} では l=0l=0 と l=1l=1 だけを残せばよい。 Z≃1+3e−2Θrot/T. Z\simeq 1+3e^{-2\Theta_{\mathrm{rot}}/T}. 内部エネルギーは U≃3(2kBΘrot)e−2Θrot/T1+3e−2Θrot/T≃6kBΘrote−2Θrot/T. U\simeq \frac{3(2k_{\mathrm{B}}\Theta_{\mathrm{rot}})e^{-2\Theta_{\mathrm{rot}}/T}} {1+3e^{-2\Theta_{\mathrm{rot}}/T}} \simeq 6k_{\mathrm{B}}\Theta_{\mathrm{rot}}e^{-2\Theta_{\mathrm{rot}}/T}. したがって C≃12kB(ΘrotT)2e−2Θrot/T. C\simeq 12k_{\mathrm{B}}\left(\frac{\Theta_{\mathrm{rot}}}{T}\right)^2 e^{-2\Theta_{\mathrm{rot}}/T}. 比熱は等分配則の値へ滑らかに行くのではなく,励起準位が熱的に使えるようになるまで指数的に小さい。
  4. 水素分子では二つの陽子がフェルミ粒子である。電子・振動の基底状態を対称と見れば,全波動関数の交換反対称性は 核スピン部分×回転部分 \text{核スピン部分}\times\text{回転部分} で決まる。核スピン三重項は対称であり,回転波動関数は反対称,すなわち奇数 ll を取る。これがオルト水素である。核スピン一重項は反対称であり,回転波動関数は対称,すなわち偶数 ll を取る。これがパラ水素である。 高温では偶数 ll と奇数 ll の回転状態数がほぼ等しく,核スピン縮退度だけが効くので Northo:Npara=3:1 N_{\mathrm{ortho}}:N_{\mathrm{para}}=3:1 となる。
  5. 急冷でオルト・パラ比が 3:13:1 のまま凍結されたとする。T=30 KT=30\,\mathrm{K} では Θrot=90 K\Theta_{\mathrm{rot}}=90\,\mathrm{K} より,パラ水素は最低の l=0l=0,オルト水素は最低の l=1l=1 にほぼ入る。回転エネルギーは El=0=0,El=1=2kBΘrot=180 kBK. E_{l=0}=0,\qquad E_{l=1}=2k_{\mathrm{B}}\Theta_{\mathrm{rot}}=180\,k_{\mathrm{B}}\mathrm{K}. したがって一分子あたり平均内部エネルギーは UkB≃34⋅180 K=135 K≃1.4×102 K. \frac{U}{k_{\mathrm{B}}} \simeq \frac34\cdot 180\,\mathrm{K} =135\,\mathrm{K} \simeq 1.4\times 10^2\,\mathrm{K}.

最終答

% 古典分配関数は Zcl=8π2IkBTh2=TΘrot,Θrot=ℏ22IkB, Z_{\mathrm{cl}}=\frac{8\pi^2Ik_{\mathrm{B}}T}{h^2} =\frac{T}{\Theta_{\mathrm{rot}}},\qquad \Theta_{\mathrm{rot}}=\frac{\hbar^2}{2Ik_{\mathrm{B}}}, したがって U=kBT, C=kBU=k_{\mathrm{B}}T,\ C=k_{\mathrm{B}}。量子回転子では Z=∑l=0∞(2l+1)e−Θrotl(l+1)/T. Z=\sum_{l=0}^\infty (2l+1)e^{-\Theta_{\mathrm{rot}}l(l+1)/T}. 低温では C≃12kB(Θrot/T)2e−2Θrot/T. C\simeq 12k_{\mathrm{B}}(\Theta_{\mathrm{rot}}/T)^2e^{-2\Theta_{\mathrm{rot}}/T}. 水素分子はオルトが奇数 ll,パラが偶数 ll を取り,高温比は 3:13:1。急冷後 30 K30\,\mathrm{K} の平均回転エネルギーは U/kB≃1.4×102 K. U/k_{\mathrm{B}}\simeq 1.4\times 10^2\,\mathrm{K}.

第3問 — 電磁気学:ローレンツ振動子と屈折率分散

分極ではなく電流密度から入る

この設問は束縛電子の速度を使って電流密度を作らせている。分極 P=NerP=Ne\boldsymbol r から解いても同じ結果になるが,アンペール・マクスウェルの式に戻す流れを明示すると係数を落としにくい。

共鳴の位置が色順を決める

水晶では可視域が共鳴より低振動数側にあり,屈折率は振動数とともに増える。共鳴が可視域に入ると曲線が分断され,プリズム上の並びが単純な赤から紫の順ではなくなる。

解答

束縛電子の運動から電流密度を作り,それをマクスウェル方程式に戻す。

  1. 電子の変位を r\boldsymbol r,外部電場を E=E0cos⁡ωt\boldsymbol E=\boldsymbol E_0\cos\omega t とする。復元力を −mω02r-m\omega_0^2\boldsymbol r とすれば md2rdt2=−mω02r+eE0cos⁡ωt m\frac{d^2\boldsymbol r}{dt^2} = -m\omega_0^2\boldsymbol r +e\boldsymbol E_0\cos\omega t である。ここで e<0e<0 であるが,最終的な分散式には e2e^2 が現れる。
  2. 定常解を r=Rcos⁡ωt \boldsymbol r=\boldsymbol R\cos\omega t とおくと m(ω02−ω2)R=eE0,R=em(ω02−ω2)E0. m(\omega_0^2-\omega^2)\boldsymbol R=e\boldsymbol E_0, \qquad \boldsymbol R= \frac{e}{m(\omega_0^2-\omega^2)}\boldsymbol E_0. 速度は v=−ωRsin⁡ωt. \boldsymbol v =-\omega\boldsymbol R\sin\omega t. 単位体積あたりの束縛電子数を NN とすれば,電流密度は id=Nev=−Ne2ωm(ω02−ω2)E0sin⁡ωt. \boldsymbol i_d =Ne\boldsymbol v = -\frac{Ne^2\omega}{m(\omega_0^2-\omega^2)} \boldsymbol E_0\sin\omega t.
  3. 真空中のアンペール・マクスウェルの式を 1μ0∇×B=ϵ0∂E∂t+i \frac{1}{\mu_0}\nabla\times\boldsymbol B = \epsilon_0\frac{\partial\boldsymbol E}{\partial t} +\boldsymbol i と書く。両辺の発散を取ると左辺は常にゼロなので 0=ϵ0∂∂t(∇⋅E)+∇⋅i. 0=\epsilon_0\frac{\partial}{\partial t}(\nabla\cdot\boldsymbol E) +\nabla\cdot\boldsymbol i. ガウスの法則 ∇⋅E=ρ/ϵ0\nabla\cdot\boldsymbol E=\rho/\epsilon_0 を使えば ∂ρ∂t+∇⋅i=0 \frac{\partial\rho}{\partial t}+\nabla\cdot\boldsymbol i=0 となり,電荷保存が従う。
  4. 真空で伝導電流がない場合,E=E0cos⁡ωt\boldsymbol E=\boldsymbol E_0\cos\omega t に対して ∇×B=−μ0ϵ0ωE0sin⁡ωt. \nabla\times\boldsymbol B = -\mu_0\epsilon_0\omega\boldsymbol E_0\sin\omega t. 誘電体中では束縛電子の電流密度を加えるので ∇×B=−μ0ω[ϵ0+Ne2m(ω02−ω2)]E0sin⁡ωt. \nabla\times\boldsymbol B = -\mu_0\omega \left[ \epsilon_0+ \frac{Ne^2}{m(\omega_0^2-\omega^2)} \right]\boldsymbol E_0\sin\omega t. したがって同じ形 ∇×B=−CωE0sin⁡ωt \nabla\times\boldsymbol B=-{\cal C}\omega\boldsymbol E_0\sin\omega t で書けば C=μ0[ϵ0+Ne2m(ω02−ω2)]. {\cal C}=\mu_0 \left[ \epsilon_0+ \frac{Ne^2}{m(\omega_0^2-\omega^2)} \right].
  5. 屈折率は波の位相速度 vph=1/Cv_{\mathrm{ph}}=1/\sqrt{{\cal C}} と c=1/ϵ0μ0c=1/\sqrt{\epsilon_0\mu_0} の比から n(ω)=cvph=Cϵ0μ0=1+Ne2mϵ0(ω02−ω2). n(\omega)=\frac{c}{v_{\mathrm{ph}}} =\sqrt{\frac{{\cal C}}{\epsilon_0\mu_0}} = \sqrt{ 1+ \frac{Ne^2}{m\epsilon_0(\omega_0^2-\omega^2)} }. プラズマ振動数 ωp=Ne2mϵ0 \omega_p=\sqrt{\frac{Ne^2}{m\epsilon_0}} を用いれば n(ω)=1+ωp2ω02−ω2. n(\omega)= \sqrt{1+\frac{\omega_p^2}{\omega_0^2-\omega^2}}.

最終答

% 電子の運動方程式は mr¨+mω02r=eE0cos⁡ωt. m\ddot{\boldsymbol r}+m\omega_0^2\boldsymbol r =e\boldsymbol E_0\cos\omega t. 定常解から id=−Ne2ωm(ω02−ω2)E0sin⁡ωt. \boldsymbol i_d= -\frac{Ne^2\omega}{m(\omega_0^2-\omega^2)} \boldsymbol E_0\sin\omega t. 誘電体中では ∇×B=−μ0ω[ϵ0+Ne2m(ω02−ω2)]E0sin⁡ωt. \nabla\times\boldsymbol B =-\mu_0\omega \left[ \epsilon_0+\frac{Ne^2}{m(\omega_0^2-\omega^2)} \right]\boldsymbol E_0\sin\omega t. よって n(ω)=1+ωp2ω02−ω2,ωp2=Ne2mϵ0. n(\omega)= \sqrt{1+\frac{\omega_p^2}{\omega_0^2-\omega^2}}, \qquad \omega_p^2=\frac{Ne^2}{m\epsilon_0}. 可視域より高い共鳴なら通常分散,可視域近傍の共鳴なら異常分散で色順が乱れる。

第4問 — 原子核・素粒子:しきい値・飛行時間・寿命測定

しきい値は不変量で処理する

実験室系では生成粒子が全て静止することはできない。重心系で相対運動が消えるという条件を ss に入れるのが,しきい値問題の標準手順である。

指数分布の寿命推定

崩壊時間が指数分布に従う場合,寿命の最尤推定量は標本平均である。個々の長い崩壊時間を外れ値として捨てるのではなく,指数分布の長い尾として含める。

解答

相対論的不変量,検出器の偶然同時計数,指数分布の最尤推定をまとめて扱う。

  1. 標的陽子が静止している実験室系で p+p→p+p+π0 p+p\to p+p+\pi^0 を考える。入射陽子の全エネルギーを ElabE_{\mathrm{lab}} とすると s=(p1+p2)2=2mp2c4+2mpc2Elab. s=(p_1+p_2)^2 =2m_p^2c^4+2m_pc^2E_{\mathrm{lab}}. しきい値では重心系で生成粒子が相対運動を持たないため sthr=(2mp+mπ)2c4. s_{\mathrm{thr}}=(2m_p+m_\pi)^2c^4. したがって入射陽子の運動エネルギーは Kthr=(2mp+mπ)2−2mp22mpc2−mpc2=(2mπ+mπ22mp)c2. \begin{aligned} K_{\mathrm{thr}} &= \frac{(2m_p+m_\pi)^2-2m_p^2}{2m_p}c^2 -m_pc^2 \\ &= \left(2m_\pi+\frac{m_\pi^2}{2m_p}\right)c^2. \end{aligned} mpc2=938 MeVm_pc^2=938\,\mathrm{MeV}, mπc2=135 MeVm_\pi c^2=135\,\mathrm{MeV} を使うと Kthr≃270 MeV+13522⋅938 MeV≃2.8×102 MeV. K_{\mathrm{thr}}\simeq 270\,\mathrm{MeV} +\frac{135^2}{2\cdot 938}\,\mathrm{MeV} \simeq 2.8\times 10^2\,\mathrm{MeV}.
  2. π0\pi^0 の崩壊で生じる高エネルギー光子は,物質中で主に以下の過程を起こす。 光電効果,コンプトン散乱,電子陽電子対生成. \text{光電効果},\qquad \text{コンプトン散乱},\qquad \text{電子陽電子対生成}. 数十 MeV 以上の領域では対生成が特に重要であり,その後の電磁シャワーを検出することで光子を捉える。
  3. 飛行距離と飛行時間から v=ℓt=1.200 m5.000 ns=2.40×108 m/s=0.800c. v=\frac{\ell}{t} = \frac{1.200\,\mathrm{m}}{5.000\,\mathrm{ns}} = 2.40\times 10^8\,\mathrm{m/s} =0.800c. 相対誤差は δvv=(0.0121.200)2+(0.1205.000)2≃2.6×10−2. \frac{\delta v}{v} = \sqrt{ \left(\frac{0.012}{1.200}\right)^2+ \left(\frac{0.120}{5.000}\right)^2 } \simeq 2.6\times 10^{-2}. したがって v=(2.40±0.06)×108 m/s=(0.800±0.021)c. v=(2.40\pm0.06)\times10^8\,\mathrm{m/s} =(0.800\pm0.021)c.
  4. 運動量 p=4.0 GeV/cp=4.0\,\mathrm{GeV}/c の粒子について,長さ ℓ\ell の飛行時間は t=ℓβc=ℓc1+(mcp)2. t=\frac{\ell}{\beta c} = \frac{\ell}{c} \sqrt{1+\left(\frac{mc}{p}\right)^2}. p≫mcp\gg mc なので t≃ℓc[1+12(mcp)2]. t\simeq \frac{\ell}{c} \left[ 1+\frac12\left(\frac{mc}{p}\right)^2 \right]. したがって KK 中間子と π\pi 中間子の時間差は Δt≃ℓ2cmK2−mπ2p2. \Delta t \simeq \frac{\ell}{2c} \frac{m_K^2-m_\pi^2}{p^2}. mKc2≃0.49 GeV,mπc2≃0.14 GeV m_Kc^2\simeq0.49\,\mathrm{GeV},\qquad m_\pi c^2\simeq0.14\,\mathrm{GeV} を用いると mK2−mπ2p2=0.492−0.1424.02≃1.4×10−2. \frac{m_K^2-m_\pi^2}{p^2} = \frac{0.49^2-0.14^2}{4.0^2} \simeq 1.4\times10^{-2}. Δt≥300 ps\Delta t\ge 300\,\mathrm{ps} とするには ℓ≳2cΔt1.4×10−2≃1.3×101 m. \ell \gtrsim \frac{2c\Delta t}{1.4\times10^{-2}} \simeq 1.3\times10^1\,\mathrm{m}. よって必要な飛行距離は約 13 m13\,\mathrm{m} である。
  5. 二つの計数管のノイズがそれぞれ独立に r=1.0 kHzr=1.0\,\mathrm{kHz} で入り,パルス幅を τ=50 ns\tau=50\,\mathrm{ns} とする。二つの矩形パルスが少しでも重なる条件は開始時刻差が ∣Δt∣<τ|\Delta t|<\tau である。偶然同時計数率は Racc≃2r2τ=2(103 s−1)2(50×10−9 s)=1.0×10−1 s−1. R_{\mathrm{acc}}\simeq 2r^2\tau = 2(10^3\,\mathrm{s^{-1}})^2(50\times10^{-9}\,\mathrm{s}) =1.0\times10^{-1}\,\mathrm{s^{-1}}. すなわち約 0.10 Hz0.10\,\mathrm{Hz} である。
  6. 寿命分布を f(t;τ)=1τe−t/τ(t≥0) f(t;\tau)=\frac{1}{\tau}e^{-t/\tau}\qquad(t\ge0) とする。独立な測定値 tit_i に対する対数尤度は log⁡L=−nlog⁡τ−1τ∑i=1nti. \log L=-n\log\tau-\frac{1}{\tau}\sum_{i=1}^n t_i. ∂log⁡L/∂τ=0\partial\log L/\partial\tau=0 から τ^=1n∑i=1nti. \hat\tau=\frac{1}{n}\sum_{i=1}^n t_i. 与えられた九つの値の和は 0.34+0.74+4.76+1.81+1.35+9.80+0.65+0.62+2.43=22.50 μs. 0.34+0.74+4.76+1.81+1.35+9.80+0.65+0.62+2.43 =22.50\,\mu\mathrm{s}. したがって τ^=22.509 μs=2.50 μs. \hat\tau=\frac{22.50}{9}\,\mu\mathrm{s} =2.50\,\mu\mathrm{s}.

最終答

% Kthr=(2mπ+mπ22mp)c2≃2.8×102 MeV. K_{\mathrm{thr}}= \left(2m_\pi+\frac{m_\pi^2}{2m_p}\right)c^2 \simeq 2.8\times10^2\,\mathrm{MeV}. 高エネルギー光子の主過程は光電効果・コンプトン散乱・対生成で,高エネルギーでは対生成が重要。 v=(2.40±0.06)×108 m/s. v=(2.40\pm0.06)\times10^8\,\mathrm{m/s}. 4.0 GeV/c4.0\,\mathrm{GeV}/c の π/K\pi/K 識別で 300 ps300\,\mathrm{ps} を得るには ℓ≃13 m. \ell\simeq 13\,\mathrm{m}. 偶然同時計数率は 0.10 Hz0.10\,\mathrm{Hz}。寿命の最尤推定値は τ^=2.50 μs. \hat\tau=2.50\,\mu\mathrm{s}.

第5問 — 物性実験:中性子散乱と拡散係数

横軸の単位に注意

式のローレンツ幅は角周波数の幅 Γω=DQ2\Gamma_\omega=DQ^2 である。一方,実験図は ℏω\hbar\omega を横軸にしているため,読み取った幅は ΓE=ℏDQ2\Gamma_E=\hbar DQ^2 である。

グラフ読み取りは傾きで安定化する

各ピークの半値半幅を個別に使うより,ΓE\Gamma_E と Q2Q^2 の直線性を使って傾きを読む方が誤差に強い。三点がほぼ原点を通る直線に乗ることも,ブラウン拡散モデルの検算になる。

解答

飛行時間法の運動学,ローレンツ型スペクトル,アレニウス解析を一つずつ処理する。数値は有効数字二桁で示す。

  1. グラファイトの面間隔を d=3.4 A˚d=3.4\,\text{\AA},一次反射を用いる。ブラッグ条件 λ=2dsin⁡θB \lambda=2d\sin\theta_{\mathrm{B}} であり,入射角・反射角は 45∘45^\circ なので λ=2(3.4 A˚)sin⁡45∘≃4.8 A˚. \lambda=2(3.4\,\text{\AA})\sin45^\circ \simeq 4.8\,\text{\AA}. 弾性散乱では ∣ki∣=∣kf∣=2πλ. |\boldsymbol k_i|=|\boldsymbol k_f|=\frac{2\pi}{\lambda}. 散乱角が 2θ2\theta であるから Q=∣ki−kf∣=22πλsin⁡θ=4πλsin⁡θ. Q=|\boldsymbol k_i-\boldsymbol k_f| = 2\frac{2\pi}{\lambda}\sin\theta = \frac{4\pi}{\lambda}\sin\theta.
  2. 入射中性子のエネルギーは Ei=pi22m=h22mλ2. E_i=\frac{p_i^2}{2m} =\frac{h^2}{2m\lambda^2}. 散乱後,中性子が試料から検出器まで距離 LL を時間 Δt\Delta t で飛ぶなら Ef=12m(LΔt)2. E_f=\frac12m\left(\frac{L}{\Delta t}\right)^2. 試料へのエネルギー移行を ΔE=Ei−Ef\Delta E=E_i-E_f と取ると ΔE=h22mλ2−mL22(Δt)2. \Delta E= \frac{h^2}{2m\lambda^2} -\frac{mL^2}{2(\Delta t)^2}.
  3. 時間相関関数が C(Q,t)=ae−DQ2∣t∣ C(Q,t)=ae^{-DQ^2|t|} である。偶関数性を使うと I(Q,ω)=12π∫−∞∞ae−DQ2∣t∣eiωt dt=aπ∫0∞e−DQ2tcos⁡ωt dt=aπDQ2(DQ2)2+ω2. \begin{aligned} I(Q,\omega) &= \frac{1}{2\pi}\int_{-\infty}^{\infty} ae^{-DQ^2|t|}e^{i\omega t}\,dt \\ &= \frac{a}{\pi}\int_0^\infty e^{-DQ^2t}\cos\omega t\,dt \\ &= \frac{a}{\pi} \frac{DQ^2}{(DQ^2)^2+\omega^2}. \end{aligned} よって準弾性ピークはローレンツ型になる。
  4. 最大値は ω=0\omega=0 で I(Q,0)=aπDQ2. I(Q,0)=\frac{a}{\pi DQ^2}. 半値条件は aπDQ2(DQ2)2+ω2=12aπDQ2. \frac{a}{\pi}\frac{DQ^2}{(DQ^2)^2+\omega^2} = \frac12\frac{a}{\pi DQ^2}. したがって ω2=(DQ2)2,Γω=DQ2. \omega^2=(DQ^2)^2, \qquad \Gamma_\omega=DQ^2. 横軸をエネルギー ℏω\hbar\omega で読むときの半値半幅は ΓE=ℏDQ2 \Gamma_E=\hbar DQ^2 である。
  5. アレニウス式 D=D0exp⁡(−ΔEaRT) D=D_0\exp\left(-\frac{\Delta E_a}{RT}\right) の対数を取ると ln⁡D=ln⁡D0−ΔEaR1T. \ln D=\ln D_0-\frac{\Delta E_a}{R}\frac{1}{T}. 端の二点を用いると ln⁡1.80×10−51.10×10−5≃ln⁡(1.64)≃0.49, \ln\frac{1.80\times10^{-5}}{1.10\times10^{-5}} \simeq \ln(1.64) \simeq 0.49, Δ(1T)=3.37×10−3−3.64×10−3=−2.7×10−4 K−1. \Delta\left(\frac1T\right) = 3.37\times10^{-3}-3.64\times10^{-3} = -2.7\times10^{-4}\,\mathrm{K^{-1}}. したがって直線の傾きは約 −1.8×103 K -1.8\times10^3\,\mathrm{K} であり ΔEa≃R(1.8×103)=8.3×1.8×103 J/mol≃1.5×104 J/mol=15 kJ/mol. \Delta E_a \simeq R(1.8\times10^3) = 8.3\times1.8\times10^3\,\mathrm{J/mol} \simeq 1.5\times10^4\,\mathrm{J/mol} = 15\,\mathrm{kJ/mol}. この活性化エネルギーは,水分子が拡散する際に水素結合ネットワークを組み替えるためのエネルギーと主に関係する。

最終答

% λ=4.8 A˚,Q=4πλsin⁡θ. \lambda=4.8\,\text{\AA},\qquad Q=\frac{4\pi}{\lambda}\sin\theta. ΔE=h22mλ2−mL22(Δt)2. \Delta E= \frac{h^2}{2m\lambda^2} -\frac{mL^2}{2(\Delta t)^2}. I(Q,ω)=aπDQ2(DQ2)2+ω2,Γω=DQ2,ΓE=ℏDQ2. I(Q,\omega)= \frac{a}{\pi}\frac{DQ^2}{(DQ^2)^2+\omega^2}, \qquad \Gamma_\omega=DQ^2,\quad \Gamma_E=\hbar DQ^2. 図の読み取りから D≃7.6×10−5 cm2/s. D\simeq 7.6\times10^{-5}\,\mathrm{cm^2/s}. 水の温度依存性から ΔEa≃15 kJ/mol, \Delta E_a\simeq 15\,\mathrm{kJ/mol}, 主に水素結合の組み替えエネルギーを反映する。

第6問 — 電磁気学:同軸ケーブルと整合減衰器

単位長さあたりの量で統一する

ここでの L,CL,C は総インダクタンス・総容量ではなく単位長さあたりの線路定数である。伝送速度 1/LC1/\sqrt{LC} がケーブル長に依存しないのはこのためである。

T型減衰器は二条件で決まる

減衰だけなら抵抗比はいくらでも作れるが,反射させない条件が入力インピーダンスを ZZ に固定する。減衰比 1/21/2 と整合条件を同時に解くと R1=Z/3, R2=4Z/3R_1=Z/3,\ R_2=4Z/3 が一意に出る。

解答

同軸ケーブルの断面内の静電場・静磁場から線路定数を求め,最後に整合したT型減衰器を設計する。

  1. 中心導体半径を aa,外部導体内半径を bb とする。単位長さあたりの電荷を λ\lambda とすれば,ガウスの法則から Er(r)=λ2πϵr(a<r<b) E_r(r)=\frac{\lambda}{2\pi\epsilon r} \qquad (a<r<b) である。電位差は V=∫abEr(r) dr=λ2πϵln⁡ba. V=\int_a^b E_r(r)\,dr = \frac{\lambda}{2\pi\epsilon}\ln\frac{b}{a}. したがって単位長さあたりの電気容量は C=λV=2πϵln⁡(b/a). C=\frac{\lambda}{V} = \frac{2\pi\epsilon}{\ln(b/a)}. 電流 II が流れると,アンペールの法則から Bφ(r)=μI2πr. B_\varphi(r)=\frac{\mu I}{2\pi r}. 単位長さあたりの磁場エネルギーは ∫abBφ22μ 2πr dr=μI24πln⁡ba. \int_a^b \frac{B_\varphi^2}{2\mu}\,2\pi r\,dr = \frac{\mu I^2}{4\pi}\ln\frac{b}{a}. これを (1/2)LI2(1/2)LI^2 と比べると L=μ2πln⁡ba. L=\frac{\mu}{2\pi}\ln\frac{b}{a}.
  2. 長さ Δx\Delta x の微小区間について,直列インダクタンスは LΔxL\Delta x,並列容量は CΔxC\Delta x である。直列インダクタンスによる電圧降下は ΔV=V(x+Δx,t)−V(x,t)=−LΔx ∂I∂t. \Delta V = V(x+\Delta x,t)-V(x,t) = -L\Delta x\,\frac{\partial I}{\partial t}. Δx→0\Delta x\to0 で ∂V∂x=−L∂I∂t. \frac{\partial V}{\partial x} = -L\frac{\partial I}{\partial t}. また,微小容量に流れ込む電流は CΔx ∂V∂t C\Delta x\,\frac{\partial V}{\partial t} なので,電流の減少分は ΔI=I(x+Δx,t)−I(x,t)=−CΔx ∂V∂t. \Delta I = I(x+\Delta x,t)-I(x,t) = -C\Delta x\,\frac{\partial V}{\partial t}. したがって ∂I∂x=−C∂V∂t. \frac{\partial I}{\partial x} = -C\frac{\partial V}{\partial t}.
  3. 二つの伝送線路方程式をもう一度微分して組み合わせる。 ∂2V∂x2=−L∂∂t(∂I∂x)=LC∂2V∂t2. \frac{\partial^2 V}{\partial x^2} = -L\frac{\partial}{\partial t} \left(\frac{\partial I}{\partial x}\right) = LC\frac{\partial^2V}{\partial t^2}. 同様に ∂2I∂x2=LC∂2I∂t2. \frac{\partial^2 I}{\partial x^2} = LC\frac{\partial^2I}{\partial t^2}. 標準形 ∂2V∂x2=1v2∂2V∂t2 \frac{\partial^2 V}{\partial x^2} = \frac{1}{v^2}\frac{\partial^2V}{\partial t^2} と比べて v=1LC. v=\frac{1}{\sqrt{LC}}.
  4. 絶縁体の比誘電率が 22,比透磁率が 11 なら ϵ=2ϵ0,μ=μ0. \epsilon=2\epsilon_0,\qquad \mu=\mu_0. 上で求めた L,CL,C を使うと v=1μϵ=12μ0ϵ0=c2≃0.7c. v=\frac{1}{\sqrt{\mu\epsilon}} = \frac{1}{\sqrt{2\mu_0\epsilon_0}} = \frac{c}{\sqrt2} \simeq 0.7c. したがって光速の約 70%70\% である。
  5. 進行波 V=V0eiω(t±x/v),I=I0eiω(t±x/v) V=V_0e^{i\omega(t\pm x/v)},\qquad I=I_0e^{i\omega(t\pm x/v)} を伝送線路方程式に代入する。例えば +x+x 方向波では ∂V∂x=−iωvV=−L(iω)I. \frac{\partial V}{\partial x} = -\frac{i\omega}{v}V = -L(i\omega)I. したがって VI=Lv=L1LC=LC. \frac{V}{I}=Lv=L\frac{1}{\sqrt{LC}}=\sqrt{\frac{L}{C}}. これが特性インピーダンスである。断面から求めた L,CL,C を代入すると Z=LC=12πμϵln⁡ba. Z=\sqrt{\frac{L}{C}} = \frac{1}{2\pi}\sqrt{\frac{\mu}{\epsilon}}\ln\frac{b}{a}. ケーブル長ではなく,断面形状と媒質だけで決まる。

最終答

% C=2πϵln⁡(b/a),L=μ2πln⁡ba. C=\frac{2\pi\epsilon}{\ln(b/a)},\qquad L=\frac{\mu}{2\pi}\ln\frac{b}{a}. 伝送線路方程式は ∂V∂x=−L∂I∂t,∂I∂x=−C∂V∂t, \frac{\partial V}{\partial x}=-L\frac{\partial I}{\partial t}, \qquad \frac{\partial I}{\partial x}=-C\frac{\partial V}{\partial t}, したがって v=1LC=1μϵ. v=\frac{1}{\sqrt{LC}}=\frac{1}{\sqrt{\mu\epsilon}}. ϵr=2,μr=1\epsilon_r=2,\mu_r=1 では v≃0.7cv\simeq0.7c。特性インピーダンスは Z=LC=12πμϵln⁡ba. Z=\sqrt{\frac{L}{C}} =\frac{1}{2\pi}\sqrt{\frac{\mu}{\epsilon}}\ln\frac{b}{a}. 電圧を 1/21/2 にする整合T型減衰器は R1=Z3,R2=4Z3. R_1=\frac{Z}{3},\qquad R_2=\frac{4Z}{3}.

東大 専門科目(物理学) 院試 過去問の収録15年度

  • 2026年度(全4問)

    量子力学:周期ポテンシャルとバンドギャップ / 統計力学:1次元高分子とイジング模型 / 力学・電磁気:回転管内の荷電質点系

  • 2025年度(全4問)

    量子力学:振動摂動と和則 / 統計力学:調和トラップ中のボース気体 / 電磁気学:ポアソン方程式と鏡像法

  • 2024年度(全4問)

    量子力学:反転調和振動子とスクイーズド状態 / 統計力学:理想気体と吸着 / 電磁気:レイリー散乱と光ピンセット

  • 2023年度(全4問)

    量子力学:2次元調和振動子と角運動量 / 統計力学:ワイル型分散と低温比熱 / 電磁気学:多重極展開

  • 2022年度(全4問)

    量子力学:不確定性関係と最小不確定波束 / 統計力学:理想ボース気体とボース・アインシュタイン凝縮 / 電磁気学:平面電磁波の反射と透過

  • 2020年度(全6問)

    二状態系と相関 / 理想気体とフェルミ気体 / 同軸ケーブルのTEM波

  • 2019年度(全6問)

    量子力学:調和振動子と4次摂動 / 統計力学:2状態スピン鎖と転送行列 / 力学・重力:2体問題と重力波による軌道減衰

  • 2018年度(全6問)

    量子力学:2準位系と断熱遷移 / 統計力学:古典理想気体と低温量子気体 / 電磁気:振動双極子と2波源干渉

  • 2017年度(全6問)

    量子力学:2粒子振動子と急変 / 統計力学:反強磁性イジング模型の平均場 / 力学:二重振り子の固有振動と断熱追随

  • 2016年度(全6問)解答・最終答まで全文公開

    量子力学:局在磁束とゼロモード / 統計力学:調和振動子の古典極限と量子補正 / 力学:ロケット方程式と一定固有加速度

  • 2015年度(全6問)解答・最終答まで全文公開

    量子力学:円周上のデルタ関数ポテンシャル / 統計力学:高分子鎖の配位数と揺らぎ応答 / 電磁気学:運動する電荷列と場のローレンツ変換

  • 2014年度(全6問)解答・最終答まで全文公開

    量子力学:クーロン束縛状態と一様磁場中の荷電粒子 / 統計力学:二次元フェルミ粒子とギャップをもつ二帯模型 / 電磁気学:誘電体と良導体の境界における平面波

  • 2013年度(全6問)解答・最終答まで全文公開

    量子力学:調和振動子とコヒーレント状態 / 統計力学:理想ボース気体と断熱磁化 / 電磁気学:ポテンシャル、ゲージ変換、荷電粒子の束縛運動

  • 2012年度(このページ・全6問)解答・最終答まで全文公開

    量子力学:スピンと回転磁場 / 統計力学:剛体回転子とオルト・パラ水素 / 電磁気学:ローレンツ振動子と屈折率分散

  • 2011年度(全6問)解答・最終答まで全文公開

    量子力学:分離型ポテンシャルの束縛状態 / 統計力学:三角形上のスピン / 電磁気学:荷電粒子のドリフト