院試hub

筑波大学 院試 過去問 解答例

筑波大 数理物質科学研究群 物理学学位プログラム 物理 2024年度 院試 過去問 解答例・解説(全5問)

全5問。電磁気学・回路1問・熱・統計力学1問・量子力学1問。テーマタグは11件(固有値・固有ベクトル・留数定理・剛体の回転運動)。

最終更新:

このページで公開
解説5問と大問1問の途中式・最終答(全5問)
解答PDFに収録
途中式と最終答(最終答つき5問)
問題本文
非収録

筑波大 物理 2024年度 院試 過去問の出題内容(全5問)

この5問の分野は電磁気学・回路1問・熱・統計力学1問・量子力学1問・力学1問です。

大問分野主題解説の小見出し最終答
第1問—物理数学問2の直交性 / 問3の符号あり
第2問力学—平行軸と垂直軸 / 中心力の保存則あり
第3問量子力学—反エルミート行列 / 波束の運動量あり
第4問電磁気学・回路電磁気学グラフ選択の根拠 / ポテンシャルの検算あり
第5問熱・統計力学熱統計力学分配関数からエントロピーへ / 混合エントロピーの意味あり

この年度の解説には典型ミス2件が付いています。

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

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

第1問 — 物理数学

方針

問1、2、5、8、9、10は標準公式をそのまま使えるが、答案では「どの条件でその公式が使えるか」を一行添えると強い。問3、4は互いに関連しており、問3の sinc 関数の積分を問4のデルタ列として使う。

問2の直交性

有限区間の正弦基底は、端点 m=0,N+1m=0,N+1 で消える離散ラプラシアンの固有ベクトルである。異なる固有値に属する固有ベクトルが直交する、という見方をすると、余弦和を最後まで計算しなくても結果を説明できる。ただし試験答案では、少なくとも n=n′n=n' の規格化係数 (N+1)/2(N+1)/2 を明示すること。

問3の符号

経路が原点を下側によけるため、e2iz/z2e^{2iz}/z^2 は上半平面で閉じたときに原点の極を含む。一方 e−2iz/z2e^{-2iz}/z^2 は下半平面で閉じると原点を含まない。ここを逆にすると答えの符号が崩れる。2/z22/z^2 の項は留数がないだけでなく、小円弧を含めた経路積分としてゼロであることを確認しておくと安全である。

問4の検算

核 Kt(x)=sin⁡2(tx)tx2 K_t(x)=\frac{\sin^2(tx)}{t x^2} は幅が 1/t1/t、面積が π\pi の山になる。したがって極限は πδ(x)\pi\delta(x) と同じ働きをし、結果が πf(0)\pi f(0) になる。f(∞)=0f(\infty)=0 は遠方で余計な寄与が残らないための条件である。

典型ミス

∇×(∇×A⃗)\nabla\times(\nabla\times\vec A) の符号を逆に書くミスが多い。成分表示で ϵijkϵklm=δilδjm−δimδjl\epsilon_{ijk}\epsilon_{klm}=\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl} を一度書けば防げる。球座標の勾配では、θ\theta 成分に 1/r1/r、ϕ\phi 成分に 1/(rsin⁡θ)1/(r\sin\theta) が付く理由を、線素 dr⃗d\vec r から説明するのが最も確実である。

解答

指定は5問選択であるが、ここでは確認用に10問すべてを解く。

  1. ∫1∞dx1+x2=[arctan⁡x]1∞=π2−π4=π4. \int_1^\infty \frac{dx}{1+x^2} =\left[\arctan x\right]_1^\infty =\frac{\pi}{2}-\frac{\pi}{4} =\frac{\pi}{4}.
  2. Snn′=∑m=1Nsin⁡πnmN+1sin⁡πn′mN+1 S_{nn'}=\sum_{m=1}^N \sin\frac{\pi n m}{N+1}\sin\frac{\pi n' m}{N+1} とおく。積和公式より Snn′=12∑m=1N[cos⁡π(n−n′)mN+1−cos⁡π(n+n′)mN+1]. S_{nn'}=\frac12\sum_{m=1}^N \left[ \cos\frac{\pi(n-n')m}{N+1} -\cos\frac{\pi(n+n')m}{N+1} \right]. これは離散正弦関数の直交性であり、 Snn′=N+12δnn′ S_{nn'}=\frac{N+1}{2}\delta_{nn'} である。例えば n=n′n=n' では ∑m=1Nsin⁡2πnmN+1=N+12, \sum_{m=1}^N\sin^2\frac{\pi n m}{N+1} =\frac{N+1}{2}, n≠n′n\ne n' では上の2つの余弦和が打ち消し合う。
  3. 経路 CC は原点を下側によけるものとする。分解式の3つの積分を I+=∫Ce2izz2 dz,I−=∫Ce−2izz2 dz,I0=∫C2z2 dz I_+=\int_C\frac{e^{2iz}}{z^2}\,dz,\quad I_-=\int_C\frac{e^{-2iz}}{z^2}\,dz,\quad I_0=\int_C\frac{2}{z^2}\,dz と書く。 I+I_+ は上半平面で閉じる。大円弧の寄与は消え、原点の2位の極の留数は Res⁡z=0e2izz2=ddze2iz∣z=0=2i \operatorname{Res}_{z=0}\frac{e^{2iz}}{z^2} =\left.\frac{d}{dz}e^{2iz}\right|_{z=0}=2i なので I+=2πi⋅2i=−4π. I_+=2\pi i\cdot 2i=-4\pi. I−I_- は下半平面で閉じると大円弧の寄与が消え、原点は経路の外側にあるから I−=0. I_-=0. また 2/z22/z^2 は留数を持たず、原点下側の小半円を含めて寄与が相殺されるので I0=0. I_0=0. したがって ∫−∞∞(sin⁡xx)2dx=−14(I++I−−I0)=−14(−4π)=π. \int_{-\infty}^{\infty}\left(\frac{\sin x}{x}\right)^2dx = -\frac14(I_+ + I_- - I_0) =-\frac14(-4\pi)=\pi.
  4. 変数変換 u=txu=tx を用いると ∫−∞∞f(x)sin⁡2(tx)tx2 dx=∫−∞∞f ⁣(ut)(sin⁡uu)2du. \int_{-\infty}^{\infty}f(x)\frac{\sin^2(tx)}{t x^2}\,dx = \int_{-\infty}^{\infty}f\!\left(\frac{u}{t}\right) \left(\frac{\sin u}{u}\right)^2du. (sin⁡u/u)2\left(\sin u/u\right)^2 は可積分であり、問3の結果より積分値は π\pi である。よって lim⁡t→∞∫−∞∞f(x)sin⁡2(tx)tx2 dx=f(0)∫−∞∞(sin⁡uu)2du=πf(0). \lim_{t\to\infty} \int_{-\infty}^{\infty}f(x)\frac{\sin^2(tx)}{t x^2}\,dx = f(0)\int_{-\infty}^{\infty}\left(\frac{\sin u}{u}\right)^2du =\pi f(0).
  5. 成分表示で {∇×(∇×A⃗)}i=ϵijk∂jϵklm∂lAm=∂i∂mAm−∂j∂jAi. \{\nabla\times(\nabla\times\vec A)\}_i = \epsilon_{ijk}\partial_j\epsilon_{klm}\partial_l A_m = \partial_i\partial_m A_m-\partial_j\partial_j A_i. したがって ∇×(∇×A⃗)=∇(∇⋅A⃗)−(∇⋅∇)A⃗. \nabla\times(\nabla\times\vec A) =\nabla(\nabla\cdot\vec A)-(\nabla\cdot\nabla)\vec A.
  6. 位置ベクトルの微小変化は dr⃗=e⃗r dr+re⃗θ dθ+rsin⁡θ e⃗ϕ dϕ d\vec r=\vec e_r\,dr+r\vec e_\theta\,d\theta +r\sin\theta\,\vec e_\phi\,d\phi である。一方 df=∇f⋅dr⃗=∂f∂rdr+∂f∂θdθ+∂f∂ϕdϕ. df=\nabla f\cdot d\vec r =\frac{\partial f}{\partial r}dr +\frac{\partial f}{\partial\theta}d\theta +\frac{\partial f}{\partial\phi}d\phi. e⃗r,e⃗θ,e⃗ϕ\vec e_r,\vec e_\theta,\vec e_\phi は互いに直交する単位ベクトルなので、係数比較から ∇f=∂f∂re⃗r+1r∂f∂θe⃗θ+1rsin⁡θ∂f∂ϕe⃗ϕ. \nabla f = \frac{\partial f}{\partial r}\vec e_r +\frac1r\frac{\partial f}{\partial\theta}\vec e_\theta +\frac1{r\sin\theta}\frac{\partial f}{\partial\phi}\vec e_\phi. 左辺を直交座標成分で書けば、問題の式そのものになる。
  7. 曲面を F(x,y,z)=x2a2+y2a2+z2b2−1=0 F(x,y,z)=\frac{x^2}{a^2}+\frac{y^2}{a^2}+\frac{z^2}{b^2}-1=0 とおく。法線ベクトルは勾配で与えられるから ∇F=(2xa2,2ya2,2zb2). \nabla F = \left(\frac{2x}{a^2},\frac{2y}{a^2},\frac{2z}{b^2}\right). 点 (0,a/2,b/2)(0,a/\sqrt2,b/\sqrt2) では ∇F=(0,2a,2b). \nabla F = \left(0,\frac{\sqrt2}{a},\frac{\sqrt2}{b}\right). よって、これに比例する任意の非零ベクトル、例えば (0,1a,1b) \left(0,\frac{1}{a},\frac{1}{b}\right) が曲面に垂直である。
  8. 行列 JJ について直接計算すると J2=I. J^2=I. 指数関数を偶数次と奇数次に分ければ exp⁡(iα2J)=Icos⁡α2+iJsin⁡α2. \exp\left(\frac{i\alpha}{2}J\right) = I\cos\frac{\alpha}{2} +iJ\sin\frac{\alpha}{2}.
  9. 対称性より Vv⃗=λTv⃗,Vv⃗′=λ′Tv⃗′ V\vec v=\lambda T\vec v,\qquad V\vec v'=\lambda' T\vec v' である。第1式を左から (v⃗′)t(\vec v')^t でなく、後で比較しやすいように v⃗t\vec v^t から第2式へ掛けると v⃗tVv⃗′=λ′v⃗tTv⃗′. \vec v^tV\vec v'=\lambda'\vec v^tT\vec v'. また第1式を転置して右から v⃗′\vec v' を掛けると、T,VT,V が対称なので v⃗tVv⃗′=λv⃗tTv⃗′. \vec v^tV\vec v'=\lambda\vec v^tT\vec v'. 両式を引いて (λ−λ′)v⃗tTv⃗′=0. (\lambda-\lambda')\vec v^tT\vec v'=0. λ≠λ′\lambda\ne\lambda' より v⃗⋅Tv⃗′=0. \vec v\cdot T\vec v'=0.
  10. A=(abcd),A−1=1det⁡A(d−b−ca) A=\begin{pmatrix}a&b\\ c&d\end{pmatrix},\qquad A^{-1}=\frac1{\det A}\begin{pmatrix}d&-b\\ -c&a\end{pmatrix} である。したがって det⁡A Tr⁡(A−1dAdx)=d a′−b c′−c b′+a d′ \det A\,\operatorname{Tr}\left(A^{-1}\frac{dA}{dx}\right) = d\,a'-b\,c'-c\,b'+a\,d' である。一方 ddxdet⁡A=ddx(ad−bc)=a′d+ad′−b′c−bc′, \frac{d}{dx}\det A = \frac{d}{dx}(ad-bc) = a'd+ad'-b'c-bc', これは上式と一致する。よって ddxdet⁡A(x)=det⁡A(x)Tr⁡[A−1(x)dA(x)dx]. \frac{d}{dx}\det A(x) = \det A(x)\operatorname{Tr}\left[ A^{-1}(x)\frac{dA(x)}{dx} \right].

最終答

問1 π4,問2 N+12δnn′,問3 ∫−∞∞(sin⁡xx)2dx=π,問4 πf(0),問5 ∇×(∇×A⃗)=∇(∇⋅A⃗)−(∇⋅∇)A⃗,問6 ∇f=fre⃗r+1rfθe⃗θ+1rsin⁡θfϕe⃗ϕ,問7 (0,1a,1b) の任意の非零倍,問8 Icos⁡α2+iJsin⁡α2,問9 v⃗⋅Tv⃗′=0,問10 ddxdet⁡A=det⁡A Tr⁡(A−1dAdx). \begin{gathered} \text{問1 } \frac{\pi}{4},\quad \text{問2 } \frac{N+1}{2}\delta_{nn'},\quad \text{問3 } \displaystyle\int_{-\infty}^{\infty}\left(\frac{\sin x}{x}\right)^2dx=\pi,\\ \text{問4 } \pi f(0),\quad \text{問5 } \nabla\times(\nabla\times\vec A)=\nabla(\nabla\cdot\vec A)-(\nabla\cdot\nabla)\vec A,\\ \text{問6 } \displaystyle \nabla f=f_r\vec e_r+\frac1r f_\theta\vec e_\theta+ \frac1{r\sin\theta}f_\phi\vec e_\phi,\\ \text{問7 } \displaystyle \left(0,\frac1a,\frac1b\right)\ \text{の任意の非零倍},\quad \text{問8 } I\cos\frac{\alpha}{2}+iJ\sin\frac{\alpha}{2},\\ \text{問9 } \vec v\cdot T\vec v'=0,\quad \text{問10 } \displaystyle \frac{d}{dx}\det A=\det A\,\operatorname{Tr}\left(A^{-1}\frac{dA}{dx}\right). \end{gathered}

第2問 — 力学

方針

[A] は保存則と慣性モーメントの基礎事項である。[B] は「エネルギー保存」「円運動の拘束」「拘束力が消える条件」を組み合わせる問題で、未定乗数 λ\lambda が垂直抗力に対応していることを見抜くのが核心である。

平行軸と垂直軸

平行軸の定理では、加わる項は必ず全質量 MM と軸間距離の二乗 ℓ2 \ell^2 の積である。円板の面内軸まわりの慣性モーメントでは、円板の対称性により Ix=IyI_x=I_y とする一行が必要である。

中心力の保存則

角運動量保存は「力が中心向きなのでトルクがゼロ」という内容である。答案では r⃗˙×p⃗=0\dot{\vec r}\times\vec p=0 も明示すること。ここを省くと、r⃗×p⃗˙\vec r\times\dot{\vec p} だけを見た不完全な証明になりやすい。

離脱条件の意味

λ=0\lambda=0 は、輪が質点を押す垂直抗力がゼロになることを表す。値 y=2a3 y=\frac{2a}{3} は頂点 y=ay=a より低く、最下点より高いので物理的に妥当である。もし y=a/3y=a/3 や負の値が出た場合は、エネルギー式または拘束条件の2階微分の符号を疑うとよい。

典型ミス

未定乗数付きラグランジアンで、∂[−mgy]/∂y=−mg\partial[-mgy]/\partial y=-mg の符号を落とすと問4の λ\lambda が誤る。また、離脱時に速度をゼロと置いてはいけない。ゼロになるのは輪から受ける垂直抗力であり、速度ではない。

力学の途中式・最終答をPDFで見る

第3問 — 量子力学

方針

[A] は線形代数、波束、エネルギーと運動量の誤差伝播である。[B] はパウリ行列の指数関数とスピンの歳差運動で、状態ベクトルそのものよりも「全体位相は期待値に影響しない」という点を意識すると計算が短くなる。

反エルミート行列

エルミート行列の固有値が実数であることと対になる事実である。A†=−AA^\dagger=-A なら期待値 v⃗†Av⃗\vec v^\dagger A\vec v は純虚数になる。固有値の証明では、固有ベクトルのノルム v⃗†v⃗\vec v^\dagger\vec v が正であることを書いて割るのが大切である。

波束の運動量

ガウス包絡から出る −x/a2-x/a^2 の項は奇関数として消える。残るのは平面波因子 eikxe^{ikx} の波数 kk による ℏk\hbar k である。正規化定数を細かく計算しなくても、奇偶性と規格化だけで期待値は求まる。

不確定性の見積もり

E=p2/(2m)E=p^2/(2m) から ΔE=(p/m)Δp\Delta E=(p/m)\Delta p とするのは、揺らぎが平均値に比べて十分小さいときの一次近似である。ここでは ΔE/E=10−4\Delta E/E=10^{-4} なので妥当である。最後の Δx\Delta x は等号ではなく下限であり、「少なくともこの程度」と読む。

スピン時間発展の検算

zz 軸磁場で −x-x 方向のスピンが半周期だけ回って +x+x 方向になるため、問5で ⟨Sx⟩=+ℏ/2\langle S_x\rangle=+\hbar/2 となる。さらに yy 軸まわりに半回転させると初期方向へ戻るため、問6の期待値は初期状態と同じになる。ただし状態ベクトルには全体位相 ii が残るので、重なりは 11 ではなく −i-i である。

量子力学の途中式・最終答をPDFで見る

第4問 — 電磁気学

方針

[A] はガウスの法則とアンペールの法則から、内側と外側の rr 依存性を判定する問題である。選択肢を見比べる前に、まず r<Rr<R と r>Rr>R の式を出すと迷いにくい。[B] は電磁誘導によるエネルギー散逸を、減衰振動の運動方程式に読み替える問題である。

グラフ選択の根拠

一様帯電球の内部電場は包む電荷が r3r^3 に比例するため E∝rE\propto r である。球殻や導体球の内部電場はゼロだが、ポテンシャルはゼロではなく一定である。この「電場ゼロ」と「ポテンシャルゼロ」を混同しないことが重要である。

ポテンシャルの検算

ポテンシャルは E=−dV/drE=-dV/dr でつながる。球殻・導体球の内部で VV が一定なら内部電場はゼロになる。一様帯電球の内部ポテンシャルは中心で最大だが有限で、r=Rr=R で外部の 1/r1/r 型と滑らかにつながる。

誘導起電力の符号

磁束を Φ=BAsin⁡θ\Phi=BA\sin\theta と置くと、Faraday の法則から E=−dΦ/dt\mathcal E=-d\Phi/dt である。符号は正方向の取り方に依存するが、電力 P1=E2/RP_1=\mathcal E^2/R と減衰係数は符号に依存しない。運動方程式では、誘導電流の効果が必ず θ˙\dot\theta に反対向きの制動項になることを確認する。

減衰条件

θ¨+γθ˙+ω02θ=0 \ddot\theta+\gamma\dot\theta+\omega_0^2\theta=0 は、γ<2ω0\gamma<2\omega_0 のときだけ振動しながら減衰する。等号は臨界減衰、超えると非振動的に戻る。試験では、条件式だけでなく、失われたエネルギーが抵抗のジュール熱になることまで書くと問5の要求を満たす。

電磁気学の途中式・最終答をPDFで見る

第5問 — 熱統計力学

方針

[A] は熱力学ポテンシャルの定義を、分配関数で書き直すだけでよい。[B] は「異種気体では各成分が全体積を使えるようになる」ことと、「同種気体では粒子を区別して数えてはいけない」ことが中心である。

分配関数からエントロピーへ

S=−(∂F∂T)V=kBlog⁡Z+kBT∂log⁡Z∂T S=-\left(\frac{\partial F}{\partial T}\right)_V = k_B\log Z+k_BT\frac{\partial\log Z}{\partial T} の第2項を落とさないこと。理想気体では log⁡Z\log Z に (3N/2)log⁡T(3N/2)\log T が含まれるため、この微分から (3/2)NkB(3/2)Nk_B が出る。

混合エントロピーの意味

異なる気体では、混合後に気体1も気体2も体積 V1+V2V_1+V_2 を探索できる。そのため N1kBlog⁡(V/V1)+N2kBlog⁡(V/V2) N_1k_B\log(V/V_1)+N_2k_B\log(V/V_2) だけエントロピーが増える。この式は質量 m1,m2m_1,m_2 に依存しない。質量依存項は混合前後で変化しないためである。

不可逆性

断熱的に仕切りを外して異種気体を混ぜる過程では、外界に何も残さずに元の分離状態へ戻すことはできない。エントロピー増大 ΔS>0\Delta S>0 がその熱力学的な表現である。可逆経路を使うのは、実際の混合過程が可逆だからではなく、エントロピーが状態量だからである。

ギブスのパラドックス

同種気体を異種気体として数えると、存在しない混合エントロピー NkBlog⁡2Nk_B\log2 が出てしまう。正しくは混合後の同一粒子全体に対して N!N! で割る。この不可区別性の補正により、同じ温度・圧力の同種気体を混ぜてもエントロピーは変化しない。

熱統計力学の途中式・最終答をPDFで見る

筑波大 物理 院試 過去問の収録3年度