東大 工学系研究科 物理工学専攻 物理学 2013年度 院試 過去問 解答例・解説

1　円柱Bの二つの摩擦力の向きと作用点

3本の円柱の中心を正三角形に置く。下側右の円柱を $B$、上側を $C$ とし、$B$ が $C$ から受ける法線力を $N$、接線方向摩擦力の大きさを $F$、床から受ける摩擦力の大きさを $F'$ とする。

$B$ について、$C$ との接触点の位置ベクトルは $$\left(-\frac{a}{2},\,\frac{\sqrt{3}a}{2}\right)$$ である。静止するには、$C$ から受ける摩擦力は接線に沿って左下向き、床から受ける摩擦力も左向きでなければならない。中心まわりのモーメントをとると、法線力と重力、床の垂直抗力は寄与しないので $$aF-aF'=0$$ となる。したがって $$F=F'$$ である。水平方向のつり合いから

作用点と向きは次の単独のBの図の通りである。円柱間の摩擦は接点の接線に沿う。

2　静止を支える二種類の摩擦係数

上側の円柱 $C$ について左右対称性を用いると、左右の接触で法線力の大きさはともに $N$、摩擦力の大きさはともに $F$ である。鉛直方向のつり合いは $$\sqrt{3}N+F=mg.$$ よって $$F=\frac{mg}{2(2+\sqrt{3})}, \qquad N=\frac{mg}{2}.$$ $B$ が床から受ける垂直抗力を $R$ とすると $$R-mg-\frac{\sqrt{3}}{2}N-\frac{1}{2}F=0$$ であり、上の値を代入して $$R=\frac{3mg}{2}.$$ したがって静止に必要な最小摩擦係数は

3　円柱の慣性モーメントの積分

一様円柱を同軸の薄い円筒殻に分ける。密度$\rho=m/(\pi a^2l)$、半径$r$の殻の質量は$dm=2\pi\rho l r\,dr$なので

4　右向き運動の停止距離と回転

接触点は初め右向きに$v_0+a\omega_0$で滑るので、動摩擦は左向きである。$\mu''$を原題の動摩擦係数とすれば $$v(t)=v_0-\mu''gt,\qquad \omega(t)=\omega_0-2\mu''gt/a.$$ 向きを変える時刻$t_T=v_0/(\mu''g)$で

5　滑りが止まる時の左向きの速さ

純粋転がりが始まる条件は、接触点の相対速度 $$u=v+a\omega$$ が 0 になることである。したがって $$v_0+a\omega_0-3\mu''gt_R=0,\qquad t_R=\frac{v_0+a\omega_0}{3\mu''g}.$$ この時の重心速度は $$v_R=v_0-\mu''gt_R =\frac{2v_0-a\omega_0}{3}.$$ したがって左向きに戻る場合の速さは

6　Sへ戻った時に無滑りである条件

純粋転がりが始まる位置は $$x_R=v_0t_R-\frac{1}{2}\mu''gt_R^2 =\frac{(v_0+a\omega_0)(5v_0-a\omega_0)}{18\mu''g}.$$ 転がり開始後は摩擦力を必要とせず等速で進む。したがって点 $S$ を転がりながら通過するには、 $$v_R<0,\qquad x_R\ge0$$ が必要十分である。よって

である。上限の等号では、Sに戻った瞬間に接触点の滑りがちょうど0となるため、この条件に含める。

陰極を $x=0$、陽極を $x=L$ とし、電位を $\phi(x)$ と書く。電子の電荷の大きさを $e$、質量を $m$、電子密度を $n(x)$、速度を $v(x)$ とする。

1. 陰極から初速 0 で出た電子について、静電場がした仕事が運動エネルギーになる。陰極での電位を 0 とすれば $$\frac{1}{2}mv^2=e\phi.$$ よって $$v(\phi)=\sqrt{\frac{2e\phi}{m}}.$$

2. 定常一次元流では、任意の断面を単位時間に通過する電荷量が等しい。電流密度を $$i(x)=-e n(x)v(x)$$ と書くと、電荷保存より $$\frac{di}{dx}=0.$$ したがって $i(x)=-i_0$ は $x$ に依存しない。

3. ポアソン方程式は $$\frac{d^2\phi}{dx^2}=-\frac{\rho}{\varepsilon_0}$$ である。電子の電荷密度は $\rho=-en$ なので $$\frac{d^2\phi}{dx^2}=\frac{en}{\varepsilon_0}.$$ 一方、$-en v=-i_0$ より $en=i_0/v$ であるから $$\frac{d^2\phi}{dx^2} =\frac{i_0}{\varepsilon_0} \sqrt{\frac{m}{2e}}\phi^{-1/2}.$$ したがって $$A=\frac{i_0}{\varepsilon_0}\sqrt{\frac{m}{2e}}, \qquad \alpha=-\frac{1}{2}.$$

4. 電場が陰極で 0 という条件は $\phi'(0)=0$ である。 $$\phi''=A\phi^{-1/2}$$ に $\phi'$ を掛けると $$\phi''\phi'=A\phi^{-1/2}\phi'.$$ 積分して $$\frac{1}{2}(\phi')^2=2A\phi^{1/2}$$ となる。境界条件により積分定数は 0 である。よって $$\frac{d\phi}{dx}=2\sqrt{A}\,\phi^{1/4}.$$ さらに積分して $$\frac{4}{3}\phi^{3/4}=2\sqrt{A}\,x,$$ すなわち $$\phi(x)=\left(\frac{3}{2}\sqrt{A}\,x\right)^{4/3}.$$

5. $\phi(L)=V$ を代入すると $$V=\left(\frac{3}{2}\sqrt{A}\,L\right)^{4/3}, \qquad A=\frac{4}{9}\frac{V^{3/2}}{L^2}.$$ したがって電流密度の大きさは $$i_0=\frac{4}{9}\varepsilon_0 \sqrt{\frac{2e}{m}} \frac{V^{3/2}}{L^2}.$$

1　無摂動の双極子選択則

電気双極子近似では、遷移行列要素は $$\langle n'|X|n\rangle$$ に比例する。$X$ は $a+a^\dagger$ に比例するので $$\langle n'|X|n\rangle\neq 0$$ となるのは

のときだけである。

2　四次の非調和項の一次エネルギー補正

非調和項を $$H_{\mathrm{ah}}=\lambda X^4$$ とする。一次の摂動エネルギーは $$\Delta E_n^{(1)}=\lambda\langle n|X^4|n\rangle$$ である。$Y=a+a^\dagger$と置くと、$Y^2|n\rangle=\sqrt{n(n-1)}|n-2\rangle+(2n+1)|n\rangle+\sqrt{(n+1)(n+2)}|n+2\rangle$である。このノルムの二乗を取れば $$\langle n|(a+a^\dagger)^4|n\rangle =6n^2+6n+3.$$ したがって

3　両側の状態補正と相殺される遷移

$x_0=\sqrt{\hbar/(2m\omega_0)}$、$Y=a+a^\dagger$と置く。一次補正した固有状態では $$|n'\rangle=|n\rangle+\frac{\lambda x_0^4}{\hbar\omega_0} \sum_{j\ne n}\frac{\langle j|Y^4|n\rangle}{n-j}|j\rangle+O(\lambda^2).$$ 双極子行列要素では、初期状態と終状態の両方の補正を加えるので、係数は $$\langle l'|X|n'\rangle=x_0Y_{ln} +\frac{\lambda x_0^5}{\hbar\omega_0} \left[\sum_{j\ne n}\frac{Y_{lj}(Y^4)_{jn}}{n-j} +\sum_{j\ne l}\frac{(Y^4)_{lj}Y_{jn}}{l-j}\right]+O(\lambda^2).$$ $Y^4$の$\Delta n=\pm2,\pm4$と$Y$の$\pm1$から、一見$\pm5$も候補になる。しかし$l=n+5$では二項の係数が $$\sqrt{(n+1)(n+2)(n+3)(n+4)(n+5)}\left(-\frac14+\frac14\right)=0$$ となり、$l=n-5$も同様に相殺する。$\Delta n=\pm3$は相殺せず、$n=8$からの新しい行列要素は $$\langle5'|X|8'\rangle=\frac{4\sqrt{21}\lambda x_0^5}{\hbar\omega_0},\qquad \langle11'|X|8'\rangle=\frac{3\sqrt{110}\lambda x_0^5}{\hbar\omega_0}$$ である。従って一次までで許される遷移先は

低エネルギー側への放出に限れば$5,7$で、$3,13$はこの次数では許されない。

4　ハイゼンベルク表示の位置演算子

ハイゼンベルク方程式から $$\dot X=\frac{P}{m},\qquad \dot P=-m\omega_0^2X$$ である。したがって

5　数状態と振動する重ね合わせの比較

数状態では$\langle n|X(0)|n\rangle=\langle n|P(0)|n\rangle=0$なので $$\langle n|X(t)|n\rangle=0.$$ 有限振幅の例として$|\psi\rangle=(|0\rangle+|1\rangle)/\sqrt2$を取ると、$\langle X(0)\rangle=x_0$、$\langle P(0)\rangle=0$より

これは振幅が0でない調和振動である。

一粒子エネルギーの基準を基底状態で $E=0$ とする。温度 $T$、化学ポテンシャル $\mu$ に対する平均占有数は $$f(E)=\frac{1}{\exp\{(E-\mu)/(k_{\mathrm B}T)\}-1}$$ である。分母が正であるためには $$\mu\leq 0$$ でなければならない。

1. 固定した $E>0$ では、$\mu$ が大きいほど分母が小さくなり、占有数は増える。凝縮点に近づくと $\mu\to0$ となる。

2. 体積 $V$ の三次元箱で、スピン縮退を考えない一粒子状態数を運動量空間で数えると $$D(E)dE =\frac{V}{h^3}4\pi p^2\,dp, \qquad E=\frac{p^2}{2m}.$$ よって $$D(E)=\frac{V}{4\pi^2} \left(\frac{2m}{\hbar^2}\right)^{3/2}E^{1/2}.$$

3. 励起状態に入る粒子数は $$N_{\mathrm{ex}}(T,\mu) =\int_0^\infty D(E) \frac{dE}{\exp\{(E-\mu)/(k_{\mathrm B}T)\}-1}$$ である。最大値は $\mu=0$ で $$N_{\mathrm{ex}}^{\max}(T) =V\left(\frac{mk_{\mathrm B}T}{2\pi\hbar^2}\right)^{3/2} \zeta\!\left(\frac{3}{2}\right).$$ $T\to0$ ではこの値は 0 に近づく。

4. 全粒子数 $N$ が固定されているとき、凝縮温度 $T_c$ は $$N =V\left(\frac{mk_{\mathrm B}T_c}{2\pi\hbar^2}\right)^{3/2} \zeta\!\left(\frac{3}{2}\right)$$ で定まる。したがって $$T_c=\frac{2\pi\hbar^2}{mk_{\mathrm B}} \left\{\frac{N}{V\zeta(3/2)}\right\}^{2/3}.$$

5. $T<T_c$ では化学ポテンシャルは $\mu=0$ に張り付き、励起状態には $$N_{\mathrm{ex}}(T)=N\left(\frac{T}{T_c}\right)^{3/2}$$ だけ入る。残り $$N_0=N\left\{1-\left(\frac{T}{T_c}\right)^{3/2}\right\}$$ が基底状態に巨視的に占有される。

6. $T<T_c$ の内部エネルギーは、基底状態のエネルギーを 0 とすれば励起粒子だけが寄与する。 $$U =\int_0^\infty D(E) \frac{E\,dE}{e^{E/(k_{\mathrm B}T)}-1}.$$ 積分公式 $$\int_0^\infty \frac{x^{s-1}}{e^x-1}\,dx =\Gamma(s)\zeta(s)$$ を用いて $$U=\frac{3}{2} \frac{\zeta(5/2)}{\zeta(3/2)} Nk_{\mathrm B}T\left(\frac{T}{T_c}\right)^{3/2}.$$ よって体積一定の熱容量は $$C_V=\left(\frac{\partial U}{\partial T}\right)_V =\frac{15}{4} \frac{\zeta(5/2)}{\zeta(3/2)} Nk_{\mathrm B}\left(\frac{T}{T_c}\right)^{3/2}.$$

7. 液体ヘリウムを理想ボース気体として粗く見積もる。質量密度を $\rho=1.5\times10^2\,\mathrm{kg/m^3}$、原子質量を $m=6.7\times10^{-27}\,\mathrm{kg}$ とすれば $$\frac{N}{V}=\frac{\rho}{m}\simeq2.2\times10^{28}\,\mathrm{m^{-3}}.$$ $\zeta(3/2)\simeq2.61$、$\hbar\simeq1.1\times10^{-34}\,\mathrm{J\,s}$、$k_{\mathrm B}\simeq1.4\times10^{-23}\,\mathrm{J/K}$ を用いると $$T_c\simeq 3\,\mathrm{K}$$ の程度になる。実在の液体ヘリウムの転移温度からずれる主な理由は、粒子間相互作用が強く、理想気体近似が粗いことである。

角振動数を $\omega$、真空中波数を $k_0=\omega/c$ とする。屈折率 $n_2$ のコアを $|x|<d$、屈折率 $n_1$ のクラッドを $|x|>d$ とし、$n_2>n_1$ とする。伝搬方向は $z$ 方向で、場の $z$ 依存性を $e^{i(\omega t-\beta z)}$ と置く。

1. コア内の平面波成分の波数の大きさは $k_0n_2$ であり、その $z$ 成分が $\beta$ である。クラッドで指数減衰する導波モードには $$k_0n_1<\beta<k_0n_2$$ が必要である。これは幾何光学的には全反射条件に対応する。

2. TE モードでは $$\boldsymbol E=(0,E_y(x),0)e^{i(\omega t-\beta z)}$$ と書ける。ファラデーの法則から $$H_x=-\frac{\beta}{\omega\mu_0}E_y,\qquad H_z=\frac{i}{\omega\mu_0}\frac{dE_y}{dx},\qquad H_y=0$$ である。アンペールの法則を合わせると、各領域で $$\frac{d^2E_y}{dx^2} +\left(k_0^2n_j^2-\beta^2\right)E_y=0$$ が成り立つ。

3. $H_z$ は $E_y$ の $x$ 微分に比例する。もし $H_z=0$ なら $E_y$ は各領域で定数となる。しかしクラッドで減衰し、かつ境界で連続な非自明解は作れない。したがって導波するTEモードでは一般に $$H_z\neq0$$ である。

4. コア内では $$h=\sqrt{k_0^2n_2^2-\beta^2}$$ とおくと $$E_y''+h^2E_y=0.$$ クラッドでは $$\gamma=\sqrt{\beta^2-k_0^2n_1^2}$$ とおくと $$E_y''-\gamma^2E_y=0.$$

5. 偶モードを考えると $$E_y(x)= \begin{cases} A\cos hx, & |x|<d,\\ B e^{-\gamma(x-d)}, & x>d,\\ B e^{\gamma(x+d)}, & x<-d. \end{cases}$$ $x=d$ で $E_y$ と $H_z$、すなわち $dE_y/dx$ が連続であるから $$A\cos hd=B,\qquad -Ah\sin hd=-\gamma B.$$ これより $$h\tan hd=\gamma, \qquad \tan\!\left(d\sqrt{k_0^2n_2^2-\beta^2}\right) = \frac{\sqrt{\beta^2-k_0^2n_1^2}} {\sqrt{k_0^2n_2^2-\beta^2}}.$$

6. $$V=d k_0\sqrt{n_2^2-n_1^2}$$ とおく。偶モードの分散式は $$d h-\tan^{-1}\left(\frac{\gamma}{h}\right)=m\pi \qquad(m=0,1,2,\ldots)$$ と同値である。左辺は $\beta=k_0n_1$ で $V$、$\beta\to k_0n_2$ で $-\pi/2$ に単調に下がる。したがって少なくとも $m=0$ の偶モードは存在する。偶モードが一つだけである条件は $$V\leq \pi, \qquad d k_0\sqrt{n_2^2-n_1^2}\leq \pi$$ である。

質量 $m$ の原子が交互にばね定数 $C_1>C_2$ のばねでつながれ、単位胞長を $a$ とする。単位胞 $n$ の二つの変位を $u_n,v_n$ と書く。

1. 運動方程式は $$m\ddot u_n=-C_1(u_n-v_n)-C_2(u_n-v_{n-1}),$$ $$m\ddot v_n=-C_1(v_n-u_n)-C_2(v_n-u_{n+1})$$ である。平面波 $$u_n=Ue^{i(kna-\omega t)},\qquad v_n=Ve^{i(kna-\omega t)}$$ を代入すると $$m\omega^2 \begin{pmatrix}U\\V\end{pmatrix} = \begin{pmatrix} C_1+C_2 & -(C_1+C_2e^{-ika})\\ -(C_1+C_2e^{ika}) & C_1+C_2 \end{pmatrix} \begin{pmatrix}U\\V\end{pmatrix}.$$ 固有値は $$m\omega^2=C_1+C_2 \pm\sqrt{C_1^2+C_2^2+2C_1C_2\cos ka}.$$ 符号 $-$ が音響分枝、符号 $+$ が光学分枝である。

2. $k=0$ では $$\omega_{\mathrm a}(0)=0,\qquad \omega_{\mathrm o}(0)=\sqrt{\frac{2(C_1+C_2)}{m}}.$$ ゾーン端 $k=\pi/a$ では $$\omega_{\mathrm a}\!\left(\frac{\pi}{a}\right)=\sqrt{\frac{2C_2}{m}}, \qquad \omega_{\mathrm o}\!\left(\frac{\pi}{a}\right)=\sqrt{\frac{2C_1}{m}}.$$

3. 長波長で $$\cos ka\simeq1-\frac{(ka)^2}{2}$$ と展開すると、音響分枝は $$\omega_{\mathrm a}\simeq a|k|\sqrt{\frac{C_1C_2}{2m(C_1+C_2)}}.$$ よって音速は $$v_{\mathrm a} =a\sqrt{\frac{C_1C_2}{2m(C_1+C_2)}}.$$ 長波長の音響モードでは二つの原子が同位相に動き、光学モードでは逆位相に動く。

4. 光の入射波数ベクトルを $\boldsymbol K_i$、散乱後を $\boldsymbol K_f$、フォノン波数を $\boldsymbol q$ とする。フォノンを生成するストークス過程では $$\omega_i-\omega_f=\omega_{\mathrm B},\qquad \boldsymbol K_i-\boldsymbol K_f=\boldsymbol q$$ が成り立つ。光の振動数変化は小さいので $|\boldsymbol K_i|\simeq|\boldsymbol K_f|=n\omega_i/c$ と近似できる。散乱角を $\theta$ とすると $$q=|\boldsymbol K_i-\boldsymbol K_f| =2\frac{n\omega_i}{c}\sin\frac{\theta}{2}.$$

5. 長波長音響フォノンでは $$\omega_{\mathrm B}=v_{\mathrm a}q =2v_{\mathrm a}\frac{n\omega_i}{c}\sin\frac{\theta}{2}.$$

6. 後方散乱では $\theta=\pi$ である。入射光の波長を $\lambda_i=680\,\mathrm{nm}$、屈折率を $n=2.4$、ブリルアンシフトの波数表示を $1/\lambda_{\mathrm B}=3.0\,\mathrm{cm^{-1}}$ とすると $$\frac{1}{\lambda_{\mathrm B}} =\frac{2nv_{\mathrm a}}{c\lambda_i}.$$ よって $$v_{\mathrm a} =\frac{c\lambda_i}{2n\lambda_{\mathrm B}} \simeq 1.3\times10^4\,\mathrm{m/s}.$$

7. 中性子の質量を $M_n$、プランク定数を $h=2\pi\hbar$ とする。入射波長 $\lambda_i$、散乱後波長 $\lambda_f$ に対してフォノン生成のエネルギーは $$\hbar\omega =\frac{h^2}{2M_n} \left(\frac{1}{\lambda_i^2}-\frac{1}{\lambda_f^2}\right).$$ 二つの入射波長 $\lambda_{i1}<\lambda_{i2}$ がゾーン端 $q=\pi/a$ の二分枝に対応すると、短い入射波長の方が大きいエネルギーを与えるので光学分枝に対応する。したがって $$\omega_1= \frac{h^2}{2M_n\hbar} \left(\frac{1}{\lambda_{i1}^2}-\frac{1}{\lambda_f^2}\right), \qquad \omega_2= \frac{h^2}{2M_n\hbar} \left(\frac{1}{\lambda_{i2}^2}-\frac{1}{\lambda_f^2}\right)$$ と置くと $$C_1=\frac{m\omega_1^2}{2},\qquad C_2=\frac{m\omega_2^2}{2}.$$

8. すでに得た音速 $$v_{\mathrm a} =a\sqrt{\frac{C_1C_2}{2m(C_1+C_2)}}$$ を用いると $$a=v_{\mathrm a} \sqrt{\frac{2m(C_1+C_2)}{C_1C_2}}.$$ $\omega_1,\omega_2$ で書けば $$a=2v_{\mathrm a} \frac{\sqrt{\omega_1^2+\omega_2^2}}{\omega_1\omega_2}.$$