院試hub

大阪公立大学 院試 過去問 解答例

大阪公立大 理学研究科 物理学専攻 物理 2023年度 院試 過去問 解答例・解説(全4問)

全4問。テーマタグは4件(分配関数・ラグランジアン・電磁波)。

最終更新:

収録3年度分の解答PDF:大阪公立大学 理学研究科 物理学専攻 物理(¥2,580・紙面見本あり)

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

大阪公立大 物理 2023年度 院試 過去問の出題内容(全4問)

大問主題解説の小見出し最終答
第1問二体中心力と有効ポテンシャル保存量の出どころ / 重心座標で交差項が消える理由あり
第2問同軸円電流の相互作用近似の次数 / 双極子場としての確認あり
第3問大正準分布と光子気体大分配関数は規格化定数 / ボース粒子で積に分解できる理由あり
第4問ガウス波束の時間発展期待値が実数になる理由 / 最小不確定性の見抜き方あり

この年度の解説には検算1件が付いています。

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

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

第1問 — 二体中心力と有効ポテンシャル

保存量の出どころ

ポテンシャルが ∣r1−r2∣|\mathbf r_1-\mathbf r_2| だけの関数であることが本質である。空間並進に対する不変性が全運動量保存を、空間回転に対する不変性が全角運動量保存を与える。答案では対称性の言葉だけで済ませず、設問の指示通り運動方程式から右辺を打ち消す計算を書くと減点されにくい。

重心座標で交差項が消える理由

r˙1=R˙+(m2/M)r˙\dot{\mathbf r}_1=\dot{\mathbf R}+(m_2/M)\dot{\mathbf r}, r˙2=R˙−(m1/M)r˙\dot{\mathbf r}_2=\dot{\mathbf R}-(m_1/M)\dot{\mathbf r} を代入すると、 R˙⋅r˙\dot{\mathbf R}\cdot\dot{\mathbf r} の係数は m1m2/M−m2m1/M=0m_1m_2/M-m_2m_1/M=0 になる。この一行を書いておくと、換算質量 μ\mu がどこから出たかが明確になる。

運動平面の示し方

l=μr×r˙\mathbf l=\mu\mathbf r\times\dot{\mathbf r} なので常に l⋅r=0\mathbf l\cdot\mathbf r=0 である。l\mathbf l が保存して向きも固定されるため、r\mathbf r は固定平面から出られない。単に「角運動量保存なので平面運動」と書くだけでは、どの平面かが曖昧である。

有効ポテンシャルの読み方

l2/(2μr2)l^2/(2\mu r^2) は遠心力ポテンシャルで、r→0r\to0 での振る舞いを決める。したがって l=0l=0 と l≠0l\ne0 では原点近くの許容範囲が全く違う。特に l≠0l\ne0 では r=0r=0 に到達できない。

典型的な間違い

E>0E>0 のとき二次方程式の二つの根のうち一方は負である。正根より外側が許容範囲であり、二つの正の転回点を書くのは誤りである。また E=0E=0 は極限として別に書くと符号ミスを避けやすい。

解答

r=∣r1−r2∣r=|\mathbf r_1-\mathbf r_2| とおく。ラグランジアンは L=12m1r˙1 2+12m2r˙2 2−U(∣r1−r2∣) L=\frac12 m_1\dot{\mathbf r}_1^{\,2} +\frac12 m_2\dot{\mathbf r}_2^{\,2} -U(|\mathbf r_1-\mathbf r_2|) である。したがって ∂U∂r1=U′(r)r1−r2r,∂U∂r2=−U′(r)r1−r2r \frac{\partial U}{\partial \mathbf r_1} =U'(r)\frac{\mathbf r_1-\mathbf r_2}{r},\qquad \frac{\partial U}{\partial \mathbf r_2} =-U'(r)\frac{\mathbf r_1-\mathbf r_2}{r} より、運動方程式は m1r¨1=−U′(r)r1−r2r,m2r¨2=U′(r)r1−r2r. m_1\ddot{\mathbf r}_1=-U'(r)\frac{\mathbf r_1-\mathbf r_2}{r},\qquad m_2\ddot{\mathbf r}_2= U'(r)\frac{\mathbf r_1-\mathbf r_2}{r}. 右辺は互いに逆向きで大きさが等しいので作用反作用の法則が成り立つ。よって P˙=m1r¨1+m2r¨2=0 \dot{\mathbf P} =m_1\ddot{\mathbf r}_1+m_2\ddot{\mathbf r}_2=\mathbf 0 であり、全運動量 P\mathbf P は保存する。また力は相対位置 r=r1−r2\mathbf r=\mathbf r_1-\mathbf r_2 に平行であるから l˙=r1×m1r¨1+r2×m2r¨2=(r1−r2)×[−U′(r)r1−r2r]=0. \begin{aligned} \dot{\mathbf l} &=\mathbf r_1\times m_1\ddot{\mathbf r}_1 +\mathbf r_2\times m_2\ddot{\mathbf r}_2 \\ &=(\mathbf r_1-\mathbf r_2)\times \left[-U'(r)\frac{\mathbf r_1-\mathbf r_2}{r}\right] =\mathbf 0 . \end{aligned} したがって全角運動量 l\mathbf l も保存する。

全質量 M=m1+m2M=m_1+m_2、換算質量 μ=m1m2/M\mu=m_1m_2/M、重心 R=m1r1+m2r2M,r=r1−r2 \mathbf R=\frac{m_1\mathbf r_1+m_2\mathbf r_2}{M},\qquad \mathbf r=\mathbf r_1-\mathbf r_2 を用いると r1=R+m2Mr,r2=R−m1Mr \mathbf r_1=\mathbf R+\frac{m_2}{M}\mathbf r,\qquad \mathbf r_2=\mathbf R-\frac{m_1}{M}\mathbf r である。交差項は消え、 L=12MR˙ 2+12μr˙ 2−U(r) L=\frac12 M\dot{\mathbf R}^{\,2} +\frac12\mu\dot{\mathbf r}^{\,2}-U(r) となる。以後 R˙=0\dot{\mathbf R}=0 の系で考えれば l=r×μr˙. \mathbf l=\mathbf r\times\mu\dot{\mathbf r}. l\mathbf l が一定なので、l⋅r=0\mathbf l\cdot\mathbf r=0 が任意の時刻で成り立つ。したがって r\mathbf r は、原点を通り法線が l\mathbf l である固定平面内にとどまる。

球座標では r˙ 2=r˙2+r2θ˙2+r2sin⁡2θ ϕ˙2 \dot{\mathbf r}^{\,2} =\dot r^2+r^2\dot\theta^2+r^2\sin^2\theta\,\dot\phi^2 だから L=12μ(r˙2+r2θ˙2+r2sin⁡2θ ϕ˙2)−U(r). L=\frac12\mu\left( \dot r^2+r^2\dot\theta^2+r^2\sin^2\theta\,\dot\phi^2 \right)-U(r). Euler--Lagrange 方程式から μr¨=μrθ˙2+μrsin⁡2θ ϕ˙2−U′(r), \mu\ddot r =\mu r\dot\theta^2+\mu r\sin^2\theta\,\dot\phi^2-U'(r), ddt(μr2θ˙)=μr2sin⁡θcos⁡θ ϕ˙2, \frac{d}{dt}\left(\mu r^2\dot\theta\right) =\mu r^2\sin\theta\cos\theta\,\dot\phi^2, ddt(μr2sin⁡2θ ϕ˙)=0 \frac{d}{dt}\left(\mu r^2\sin^2\theta\,\dot\phi\right)=0 を得る。

運動平面を θ=π/2\theta=\pi/2 に取ると l=μr2∣ϕ˙∣,∣ϕ˙∣=lμr2. l=\mu r^2|\dot\phi|,\qquad |\dot\phi|=\frac{l}{\mu r^2}. エネルギーは E=12μr˙2+12μr2ϕ˙2+U(r)=12μr˙2+V(r) E=\frac12\mu\dot r^2+\frac12\mu r^2\dot\phi^2+U(r) =\frac12\mu\dot r^2+V(r) であり、 V(r)=U(r)+l22μr2. V(r)=U(r)+\frac{l^2}{2\mu r^2}.

U(r)=−α/r  (α>0)U(r)=-\alpha/r\;(\alpha>0) のとき V(r)=l22μr2−αr. V(r)=\frac{l^2}{2\mu r^2}-\frac{\alpha}{r}. l=0l=0 では V(r)=−α/rV(r)=-\alpha/r で、r→0r\to0 で −∞-\infty、r→∞r\to\infty で 0−0- に近づく。 l≠0l\ne0 では r→0r\to0 で +∞+\infty、r→∞r\to\infty で 0−0- に近づき、 r0=l2μα,Vmin⁡=−μα22l2 r_0=\frac{l^2}{\mu\alpha},\qquad V_{\min}=-\frac{\mu\alpha^2}{2l^2} で最小値をとる。

運動可能条件は E≥V(r),r>0 E\ge V(r),\qquad r>0 である。

l=0l=0 の場合は E+α/r≥0E+\alpha/r\ge0 であるから {0<r≤α/(−E),E<0,0<r<∞,E≥0. \begin{cases} 0<r\le \alpha/(-E), & E<0,\\ 0<r<\infty, & E\ge0. \end{cases} l≠0l\ne0 の場合、条件は Er2+αr−l22μ≥0 Er^2+\alpha r-\frac{l^2}{2\mu}\ge0 である。最小エネルギー Emin⁡=−μα22l2 E_{\min}=-\frac{\mu\alpha^2}{2l^2} を用いると {運動不可,E<Emin⁡,r=r0=l2μα,E=Emin⁡,rmin⁡≤r≤rmax⁡,Emin⁡<E<0,r≥l22μα,E=0,r≥r+,E>0, \begin{cases} \text{運動不可}, & E<E_{\min},\\[2mm] r=r_0=\dfrac{l^2}{\mu\alpha}, & E=E_{\min},\\[3mm] r_{\min} \le r\le r_{\max}, & E_{\min}<E<0,\\[2mm] r\ge\dfrac{l^2}{2\mu\alpha}, & E=0,\\[3mm] r\ge r_+, & E>0, \end{cases} ただし rmin⁡=α−α2+2El2/μ−2E,rmax⁡=α+α2+2El2/μ−2E(E<0), r_{\min}=\frac{\alpha- \sqrt{\alpha^2+2El^2/\mu}}{-2E},\qquad r_{\max}=\frac{\alpha+ \sqrt{\alpha^2+2El^2/\mu}}{-2E}\qquad(E<0), また E>0E>0 の唯一の正根を r+=−α+α2+2El2/μ2E r_+=\frac{-\alpha+\sqrt{\alpha^2+2El^2/\mu}}{2E} と書いた。

最終答

L=12MR˙ 2+12μr˙ 2−U(r),l=μr×r˙,V(r)=U(r)+l22μr2. L=\frac12 M\dot{\mathbf R}^{\,2} +\frac12\mu\dot{\mathbf r}^{\,2}-U(r),\qquad \mathbf l=\mu\mathbf r\times\dot{\mathbf r},\qquad V(r)=U(r)+\frac{l^2}{2\mu r^2}. U=−α/rU=-\alpha/r では上の運動可能範囲が答である。

第2問 — 同軸円電流の相互作用

近似の次数

ΔB\Delta\mathbf B の前にすでに a1a_1 が掛かっている。全体として最初に残る非零項は a12a_1^2 であるため、Bx,ByB_x,B_y では R−3R^{-3} の一次補正が必要になる。一方、BzB_z ではベクトル成分中の +a1+a_1 と −A-A の一次補正の積が同じ次数で効く。

双極子場としての確認

小さい円電流は磁気双極子 m=Iπa2ez\mathbf m=I\pi a^2\mathbf e_z とみなせる。磁気双極子場 B=μ04πr5{3(m⋅r)r−r2m} \mathbf B=\frac{\mu_0}{4\pi r^5} \{3(\mathbf m\cdot\mathbf r)\mathbf r-r^2\mathbf m\} に代入すると、解答の B\mathbf B と一致する。係数の π\pi が消えて μ0Ia2/4\mu_0Ia^2/4 になる点はよく確認しておきたい。

力の向き

x,yx,y 成分は周回積分で打ち消される。残る zz 成分は t2×B2\mathbf t_2\times\mathbf B_2 のうち横磁場との外積から生じる。符号を取り違えると引力・斥力の結論が逆になるので、最後に「C2C_2 は +z+z 側にある。負の FzF_z は C1C_1 へ向かう」と言葉で確認する。

答案上の注意

この問題は途中式のベクトル積と三角関数積分で点を取りやすい。特に t1=(−sin⁡φ,cos⁡φ,0)\mathbf t_1=(-\sin\varphi,\cos\varphi,0) と t2=(−sin⁡θ,cos⁡θ,0)\mathbf t_2=(-\sin\theta,\cos\theta,0) を最初に明記すれば、以降の符号を追跡しやすい。

同軸円電流の相互作用の途中式・最終答をPDFで見る

第3問 — 大正準分布と光子気体

大分配関数は規格化定数

大正準分布の指数因子だけでは確率にならない。全状態について足し合わせて 1 になるように割る量が Ξ\Xi である。設問 (1) はこの規格化条件を一行で書けるかを見る問題である。

ボース粒子で積に分解できる理由

理想気体では相互作用がないため、エネルギーが各一粒子状態の占有数の和になる。そのため指数の中も ∑α\sum_\alpha に分かれ、占有数ごとの和が積に分解する。フェルミ粒子なら nα=0,1n_\alpha=0,1 で止まるが、ボース粒子では 00 から ∞\infty まで足す点が重要である。

光子の化学ポテンシャル

光子数は保存量ではないので、熱平衡では μ=0\mu=0 である。このため Ω=F−μN\Omega=F-\mu N が FF と一致する。ここを落とすと、グランドポテンシャルからヘルムホルツ自由エネルギーへ移る理由が不明確になる。

積分の検算

FF は T4VT^4V に比例し、符号は負である。エネルギー密度は E/V∝T4E/V\propto T^4 で、これは Stefan--Boltzmann 則と同じ温度依存性である。圧力が P=E/(3V)P=E/(3V) になることも、電磁波の輻射圧として自然な結果である。

典型的な間違い

状態密度には偏光の自由度 2 がすでに含まれている。ここでさらに 2 を掛けるとすべての係数が二倍になる。また ℏc\hbar c の三乗を ℏ3c\hbar^3c や ℏc3\hbar c^3 と誤記しやすいので、最後に次元を確認する。

大正準分布と光子気体の途中式・最終答をPDFで見る

第4問 — ガウス波束の時間発展

期待値が実数になる理由

位置演算子と運動量演算子が自己共役であることを、位置表示で確認する設問である。運動量では部分積分の境界項が消えることが要点であり、ここを書かないと形式的な複素共役計算だけになってしまう。

最小不確定性の見抜き方

初期波束は実ガウス包絡に平面波 eik0xe^{ik_0x} を掛けた形である。平面波因子は平均運動量を ℏk0\hbar k_0 だけずらすが、幅 Δp\Delta p はガウス包絡の幅で決まる。したがって k0k_0 が不確定性積に現れないことは自然である。

Fourier 変換の規格化

この問題では ψ\psi と ϕ\phi の両方に 1/2π1/\sqrt{2\pi} を入れる対称規約を使っている。規約を取り違えると ϕ(k,0)\phi(k,0) の前係数がずれる。指数部だけでなく前係数まで合わせるのが、規格化を保った答案として重要である。

波束の中心と広がり

xˉ(t)=(ℏk0/m)t\bar x(t)=(\hbar k_0/m)t は群速度 vg=dω/dk=ℏk0/mv_g=d\omega/dk=\hbar k_0/m で進むことを表す。一方、Δxt\Delta x_t は時間とともに増大する。これは自由粒子の分散関係 ω=ℏk2/(2m)\omega=\hbar k^2/(2m) が線形でないため、異なる kk 成分が異なる速度で進むからである。

検算

t=0t=0 を代入すると C(0)=(2πσ2)−1/4,γ(0)=σ−2,xˉ(0)=0 C(0)=(2\pi\sigma^2)^{-1/4},\qquad \gamma(0)=\sigma^{-2},\qquad \bar x(0)=0 となり、与えられた初期波動関数に戻る。また Δxt≥σ\Delta x_t\ge\sigma であり、自由発展で波束が狭くなることはない。

ガウス波束の時間発展の途中式・最終答をPDFで見る

大阪公立大 物理 院試 過去問の収録3年度

  • 2025年度(全4問)

    二重振り子と剛体振り子 / 電場・双極子・誘電体球 / 磁束を囲む荷電粒子

  • 2024年度(全4問)

    連成振動と連成振り子 / 電磁波と定常電流の磁場 / 井戸型ポテンシャルとデルタ関数摂動

  • 2023年度(このページ・全4問)

    二体中心力と有効ポテンシャル / 同軸円電流の相互作用 / 大正準分布と光子気体