院試hub

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

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

東京大学 理学系研究科 物理学専攻 専門科目(物理学) 2012年度の院試 過去問について、設問ごとの解法方針・部分点の置き所を解説。全6問収録の解答・解説PDFと併用できます。問題本文は含みません。

最終更新:

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

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

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

相互作用表示の役割

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

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

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

  1. 軌道角運動量を Lx=ypzzpy,Ly=zpxxpz,Lz=xpyypx L_x=yp_z-zp_y,\quad L_y=zp_x-xp_z,\quad L_z=xp_y-yp_x とおく。例えば [Lx,Ly]=[ypzzpy,  zpxxpz]=[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=iypx,[zpy,xpz]=x[z,pz]py=ixpy. [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(xpyypx)=iLz. [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(0ii0) S_y=\frac{\hbar}{2} \begin{pmatrix} 0&-i\\ i&0 \end{pmatrix} とすれば,[Sx,Sy]=iSz[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(11) 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)=eiHtϕ(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=12sinat,Sz=12cosat. \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ωtSysinω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)=eiH1tH2eiH1tψ(t) i\frac{d}{dt}|\psi(t)\rangle = e^{i\mathcal H_1t}\mathcal H_2e^{-i\mathcal H_1t} |\psi(t)\rangle となる。行列表現で計算すると eiH1tH2eiH1t=a02(0ei(ωa)tei(ω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_zzz 軸まわりの戻し変換と可換なので,シュレーディンガー表示でも Sz=12cos(a0t) \langle S_z\rangle=\frac12\cos(a_0t) である。共鳴では弱い横磁場の大きさ H0H_0 が振動周期を決める。

最終答

% Sy=2(0ii0). 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=12sinat,Sz=12cosat. \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)tei(ω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=aSz=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 は熱平衡比ではなく組成が凍結した比である。各成分の内部では低温の最低回転準位に落ちる,という二段階で考えると平均エネルギーを誤らない。

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

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

  1. 古典的な剛体回転子のエネルギーは E=12I(pθ2+pϕ2sin2θ) E=\frac{1}{2I}\left(p_\theta^2+\frac{p_\phi^2}{\sin^2\theta}\right) である。一分子の回転分配関数を Zcl=1h20π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/(2Isin2θ)=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=βlogZcl=kBT,C=UT=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}} では和を積分で近似でき, Zq0(2l+1)eΘrotl(l+1)/Tdl=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=0l=1l=1 だけを残せばよい。 Z1+3e2Θrot/T. Z\simeq 1+3e^{-2\Theta_{\mathrm{rot}}/T}. 内部エネルギーは U3(2kBΘrot)e2Θrot/T1+3e2Θrot/T6kBΘrote2Θ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}. したがって C12kB(ΘrotT)2e2Θ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=30KT=30\,\mathrm{K} では Θrot=90K\Theta_{\mathrm{rot}}=90\,\mathrm{K} より,パラ水素は最低の l=0l=0,オルト水素は最低の l=1l=1 にほぼ入る。回転エネルギーは El=0=0,El=1=2kBΘrot=180kBK. E_{l=0}=0,\qquad E_{l=1}=2k_{\mathrm{B}}\Theta_{\mathrm{rot}}=180\,k_{\mathrm{B}}\mathrm{K}. したがって一分子あたり平均内部エネルギーは UkB34180K=135K1.4×102K. \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}. 低温では C12kB(Θrot/T)2e2Θrot/T. C\simeq 12k_{\mathrm{B}}(\Theta_{\mathrm{rot}}/T)^2e^{-2\Theta_{\mathrm{rot}}/T}. 水素分子はオルトが奇数 ll,パラが偶数 ll を取り,高温比は 3:13:1。急冷後 30K30\,\mathrm{K} の平均回転エネルギーは U/kB1.4×102K. U/k_{\mathrm{B}}\simeq 1.4\times 10^2\,\mathrm{K}.

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

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

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

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

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

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

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

  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=ϵ0Et+i \frac{1}{\mu_0}\nabla\times\boldsymbol B = \epsilon_0\frac{\partial\boldsymbol E}{\partial t} +\boldsymbol i と書く。両辺の発散を取ると左辺は常にゼロなので 0=ϵ0t(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 に入れるのが,しきい値問題の標準手順である。

指数分布の寿命推定

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

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

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

  1. 標的陽子が静止している実験室系で p+pp+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π)22mp22mpc2mpc2=(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=938MeVm_pc^2=938\,\mathrm{MeV}, mπc2=135MeVm_\pi c^2=135\,\mathrm{MeV} を使うと Kthr270MeV+13522938MeV2.8×102MeV. 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.200m5.000ns=2.40×108m/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)22.6×102. \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)×108m/s=(0.800±0.021)c. v=(2.40\pm0.06)\times10^8\,\mathrm{m/s} =(0.800\pm0.021)c.
  4. 運動量 p=4.0GeV/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}. pmcp\gg mc なので tc[1+12(mcp)2]. t\simeq \frac{\ell}{c} \left[ 1+\frac12\left(\frac{mc}{p}\right)^2 \right]. したがって KK 中間子と π\pi 中間子の時間差は Δt2cmK2mπ2p2. \Delta t \simeq \frac{\ell}{2c} \frac{m_K^2-m_\pi^2}{p^2}. mKc20.49GeV,mπc20.14GeV m_Kc^2\simeq0.49\,\mathrm{GeV},\qquad m_\pi c^2\simeq0.14\,\mathrm{GeV} を用いると mK2mπ2p2=0.4920.1424.021.4×102. \frac{m_K^2-m_\pi^2}{p^2} = \frac{0.49^2-0.14^2}{4.0^2} \simeq 1.4\times10^{-2}. Δt300ps\Delta t\ge 300\,\mathrm{ps} とするには 2cΔt1.4×1021.3×101m. \ell \gtrsim \frac{2c\Delta t}{1.4\times10^{-2}} \simeq 1.3\times10^1\,\mathrm{m}. よって必要な飛行距離は約 13m13\,\mathrm{m} である。
  5. 二つの計数管のノイズがそれぞれ独立に r=1.0kHzr=1.0\,\mathrm{kHz} で入り,パルス幅を τ=50ns\tau=50\,\mathrm{ns} とする。二つの矩形パルスが少しでも重なる条件は開始時刻差が Δt<τ|\Delta t|<\tau である。偶然同時計数率は Racc2r2τ=2(103s1)2(50×109s)=1.0×101s1. 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.10Hz0.10\,\mathrm{Hz} である。
  6. 寿命分布を f(t;τ)=1τet/τ(t0) f(t;\tau)=\frac{1}{\tau}e^{-t/\tau}\qquad(t\ge0) とする。独立な測定値 tit_i に対する対数尤度は logL=nlogτ1τi=1nti. \log L=-n\log\tau-\frac{1}{\tau}\sum_{i=1}^n t_i. logL/τ=0\partial\log L/\partial\tau=0 から τ^=1ni=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)c22.8×102MeV. 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)×108m/s. v=(2.40\pm0.06)\times10^8\,\mathrm{m/s}. 4.0GeV/c4.0\,\mathrm{GeV}/cπ/K\pi/K 識別で 300ps300\,\mathrm{ps} を得るには 13m. \ell\simeq 13\,\mathrm{m}. 偶然同時計数率は 0.10Hz0.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_EQ2Q^2 の直線性を使って傾きを読む方が誤差に強い。三点がほぼ原点を通る直線に乗ることも,ブラウン拡散モデルの検算になる。

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

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

  1. グラファイトの面間隔を d=3.4A˚d=3.4\,\text{\AA},一次反射を用いる。ブラッグ条件 λ=2dsinθB \lambda=2d\sin\theta_{\mathrm{B}} であり,入射角・反射角は 4545^\circ なので λ=2(3.4A˚)sin454.8A˚. \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=kikf=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=EiEf\Delta E=E_i-E_f と取ると ΔE=h22mλ2mL22(Δt)2. \Delta E= \frac{h^2}{2m\lambda^2} -\frac{mL^2}{2(\Delta t)^2}.
  3. 時間相関関数が C(Q,t)=aeDQ2t C(Q,t)=ae^{-DQ^2|t|} である。偶関数性を使うと I(Q,ω)=12πaeDQ2teiωtdt=aπ0eDQ2tcosωtdt=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=0I(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) の対数を取ると lnD=lnD0ΔEaR1T. \ln D=\ln D_0-\frac{\Delta E_a}{R}\frac{1}{T}. 端の二点を用いると ln1.80×1051.10×105ln(1.64)0.49, \ln\frac{1.80\times10^{-5}}{1.10\times10^{-5}} \simeq \ln(1.64) \simeq 0.49, Δ(1T)=3.37×1033.64×103=2.7×104K1. \Delta\left(\frac1T\right) = 3.37\times10^{-3}-3.64\times10^{-3} = -2.7\times10^{-4}\,\mathrm{K^{-1}}. したがって直線の傾きは約 1.8×103K -1.8\times10^3\,\mathrm{K} であり ΔEaR(1.8×103)=8.3×1.8×103J/mol1.5×104J/mol=15kJ/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.8A˚,Q=4πλsinθ. \lambda=4.8\,\text{\AA},\qquad Q=\frac{4\pi}{\lambda}\sin\theta. ΔE=h22mλ2mL22(Δ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. 図の読み取りから D7.6×105cm2/s. D\simeq 7.6\times10^{-5}\,\mathrm{cm^2/s}. 水の温度依存性から ΔEa15kJ/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 が一意に出る。

2012年度は公開から時間が経過しているため、解答・最終答まで全文公開しています

同軸ケーブルの断面内の静電場・静磁場から線路定数を求め,最後に整合した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πϵlnba. 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πrdr=μI24πlnba. \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πlnba. 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ΔxIt. \Delta V = V(x+\Delta x,t)-V(x,t) = -L\Delta x\,\frac{\partial I}{\partial t}. Δx0\Delta x\to0Vx=LIt. \frac{\partial V}{\partial x} = -L\frac{\partial I}{\partial t}. また,微小容量に流れ込む電流は CΔxVt C\Delta x\,\frac{\partial V}{\partial t} なので,電流の減少分は ΔI=I(x+Δx,t)I(x,t)=CΔxVt. \Delta I = I(x+\Delta x,t)-I(x,t) = -C\Delta x\,\frac{\partial V}{\partial t}. したがって Ix=CVt. \frac{\partial I}{\partial x} = -C\frac{\partial V}{\partial t}.
  3. 二つの伝送線路方程式をもう一度微分して組み合わせる。 2Vx2=Lt(Ix)=LC2Vt2. \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}. 同様に 2Ix2=LC2It2. \frac{\partial^2 I}{\partial x^2} = LC\frac{\partial^2I}{\partial t^2}. 標準形 2Vx2=1v22Vt2 \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=c20.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 方向波では Vx=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πμϵlnba. Z=\sqrt{\frac{L}{C}} = \frac{1}{2\pi}\sqrt{\frac{\mu}{\epsilon}}\ln\frac{b}{a}. ケーブル長ではなく,断面形状と媒質だけで決まる。

最終答

% C=2πϵln(b/a),L=μ2πlnba. C=\frac{2\pi\epsilon}{\ln(b/a)},\qquad L=\frac{\mu}{2\pi}\ln\frac{b}{a}. 伝送線路方程式は Vx=LIt,Ix=CVt, \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 では v0.7cv\simeq0.7c。特性インピーダンスは Z=LC=12πμϵlnba. 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}.

東京大学 専門科目(物理学) — 他の年度