院試hub

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

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

全4問。力学1問。テーマタグは4件(剛体の回転運動・静電場・固有振動数と振動系)。2024年度と共通のテーマはラグランジアン。

最終更新:

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

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

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

この4問の分野は力学1問です。

大問分野主題解説の小見出し最終答
第1問力学二重振り子と剛体振り子二重振り子の最初の山場 / 固有振動数の確認あり
第2問—電場・双極子・誘電体球一様帯電球のグラフ / 双極子近似の符号あり
第3問—磁束を囲む荷電粒子演算子の順序 / 磁場ゼロとベクトルポテンシャルゼロは違うあり
第4問—平均場イジング模型平均場方程式の意味 / 二重数えの補正あり

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

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

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

大問数
2024年度 4問 → 2025年度 4問
両年度に出たテーマ
ラグランジアン
2025年度で新しく出たテーマ
剛体の回転運動・静電場・固有振動数と振動系
2024年度のページを見る

第1問 — 二重振り子と剛体振り子

二重振り子の最初の山場

座標を正しく置けば、あとは v22v_2^2 の交差項を落とさないことが重要である。 cos⁡(α1−α2)\cos(\alpha_1-\alpha_2) は微小振動では 11 に近似するが、 α˙1α˙2\dot\alpha_1\dot\alpha_2 自体は二次の微小量なので残る。 この交差項を落とすと質量行列が対角になり、二重振り子らしい連成が消えてしまう。

固有振動数の確認

ω2\omega^2 は二次方程式で決まる。判別式の中身は (l1+l2)2−4m1m1+m2l1l2=(l1−l2)2+4m2m1+m2l1l2>0 (l_1+l_2)^2-\frac{4m_1}{m_1+m_2}l_1l_2 =(l_1-l_2)^2+\frac{4m_2}{m_1+m_2}l_1l_2>0 であり、固有振動数は実数である。また次元は g/lg/l 型になっている。 ここまで確認しておくと、符号ミスをかなり見つけやすい。

剛体との対応

剛体の場合、α2\alpha_2 は二本目の糸の角ではなく剛体の向きを表す角である。 ただし重心が左端から距離 hh にあるので、重心の並進運動は 二重振り子の第二質点と同じ形をもつ。違いは、剛体の回転エネルギー 12Iα˙22\frac12I\dot\alpha_2^2 が加わる点である。 この追加分のために γ\gamma が現れ、対応条件から h=l2/γh=l_2/\gamma が出る。

基準振動の物理的意味

γ=2, l1=l2\gamma=2,\ l_1=l_2 の場合、剛体は二つの等質量点が両端に集中したダンベルに対応する。 低振動数モードでは α1\alpha_1 と α2\alpha_2 が同符号で、糸と剛体がほぼ一体となって ゆっくり揺れる。高振動数モードでは両者が逆符号で、重心の移動を抑えながら剛体が大きく 回転する。答案では「同相・逆相」だけでなく、振幅比まで書くと議論が明確になる。

解答

質点の座標を (x1,y1)=(l1sin⁡α1,−l1cos⁡α1), (x_1,y_1)=(l_1\sin\alpha_1,-l_1\cos\alpha_1), (x2,y2)=(l1sin⁡α1+l2sin⁡α2,−l1cos⁡α1−l2cos⁡α2) (x_2,y_2)=(l_1\sin\alpha_1+l_2\sin\alpha_2, -l_1\cos\alpha_1-l_2\cos\alpha_2) とおく。したがって v12=l12α˙12,v22=l12α˙12+l22α˙22+2l1l2α˙1α˙2cos⁡(α1−α2) v_1^2=l_1^2\dot\alpha_1^2,\qquad v_2^2=l_1^2\dot\alpha_1^2+l_2^2\dot\alpha_2^2 +2l_1l_2\dot\alpha_1\dot\alpha_2\cos(\alpha_1-\alpha_2) である。基準面を y=0y=0 とすれば Tm=12(m1+m2)l12α˙12+12m2l22α˙22+m2l1l2α˙1α˙2cos⁡(α1−α2), T_m=\frac12(m_1+m_2)l_1^2\dot\alpha_1^2 +\frac12m_2l_2^2\dot\alpha_2^2 +m_2l_1l_2\dot\alpha_1\dot\alpha_2\cos(\alpha_1-\alpha_2), Um=−(m1+m2)gl1cos⁡α1−m2gl2cos⁡α2. U_m=-(m_1+m_2)gl_1\cos\alpha_1-m_2gl_2\cos\alpha_2 . 微小振動では、定数項を除いて Lm=12(m1+m2)l12α˙12+12m2l22α˙22+m2l1l2α˙1α˙2−12(m1+m2)gl1α12−12m2gl2α22. L_m=\frac12(m_1+m_2)l_1^2\dot\alpha_1^2 +\frac12m_2l_2^2\dot\alpha_2^2 +m_2l_1l_2\dot\alpha_1\dot\alpha_2 -\frac12(m_1+m_2)gl_1\alpha_1^2-\frac12m_2gl_2\alpha_2^2 . これより (m1+m2)l12α¨1+m2l1l2α¨2+(m1+m2)gl1α1=0, (m_1+m_2)l_1^2\ddot\alpha_1 +m_2l_1l_2\ddot\alpha_2 +(m_1+m_2)gl_1\alpha_1=0, m2l1l2α¨1+m2l22α¨2+m2gl2α2=0. m_2l_1l_2\ddot\alpha_1 +m_2l_2^2\ddot\alpha_2 +m_2gl_2\alpha_2=0 . αi=Aieiωt\alpha_i=A_i e^{i\omega t} とおくと ∣(m1+m2)(gl1−l12ω2)−m2l1l2ω2−m2l1l2ω2m2(gl2−l22ω2)∣=0. \begin{vmatrix} (m_1+m_2)(gl_1-l_1^2\omega^2)&-m_2l_1l_2\omega^2\\ -m_2l_1l_2\omega^2&m_2(gl_2-l_2^2\omega^2) \end{vmatrix}=0. s=ω2s=\omega^2 と書けば m1l1l2s2−(m1+m2)g(l1+l2)s+(m1+m2)g2=0 m_1l_1l_2s^2-(m_1+m_2)g(l_1+l_2)s+(m_1+m_2)g^2=0 であるから、 ω±2=(m1+m2)g2m1l1l2{(l1+l2)±(l1+l2)2−4m1m1+m2l1l2}. \omega_\pm^2= \frac{(m_1+m_2)g}{2m_1l_1l_2} \left\{(l_1+l_2)\pm \sqrt{(l_1+l_2)^2-\frac{4m_1}{m_1+m_2}l_1l_2}\right\}.

次に剛体の重心を (xG,yG)=(l1sin⁡α1+hsin⁡α2,−l1cos⁡α1−hcos⁡α2) (x_G,y_G)=(l_1\sin\alpha_1+h\sin\alpha_2, -l_1\cos\alpha_1-h\cos\alpha_2) とする。重心まわりの慣性モーメントが I=(γ−1)Mh2I=(\gamma-1)Mh^2 であるので、微小振動のラグランジアンは LG=12Ml12α˙12+Ml1hα˙1α˙2+12Mγh2α˙22−12Mgl1α12−12Mghα22. L_G=\frac12M l_1^2\dot\alpha_1^2 +Ml_1h\dot\alpha_1\dot\alpha_2 +\frac12M\gamma h^2\dot\alpha_2^2 -\frac12Mgl_1\alpha_1^2-\frac12Mgh\alpha_2^2 . Lm=LGL_m=L_G が成り立つには M=m1+m2,Mh=m2l2,Mγh2=m2l22 M=m_1+m_2,\qquad Mh=m_2l_2,\qquad M\gamma h^2=m_2l_2^2 でなければならない。したがって h=l2γ,m1=(γ−1)m2,M=γm2. h=\frac{l_2}{\gamma},\qquad m_1=(\gamma-1)m_2,\qquad M=\gamma m_2 .

特に γ=2, l1=l2=l\gamma=2,\ l_1=l_2=l では m1=m2m_1=m_2, h=l/2h=l/2 であり、 ω±2=gl(2±2). \omega_\pm^2=\frac{g}{l}(2\pm\sqrt2). 対応する基準振動の振幅比は A2A1=2(g/l−ω2)ω2 \frac{A_2}{A_1}= \frac{2(g/l-\omega^2)}{\omega^2} なので、 ω−2=gl(2−2):A2=2A1,ω+2=gl(2+2):A2=−2A1. \omega_-^2=\frac{g}{l}(2-\sqrt2):\quad A_2=\sqrt2 A_1, \qquad \omega_+^2=\frac{g}{l}(2+\sqrt2):\quad A_2=-\sqrt2 A_1 .

最終答

ω±2=(m1+m2)g2m1l1l2{(l1+l2)±(l1+l2)2−4m1m1+m2l1l2},h=l2γ. \omega_\pm^2= \frac{(m_1+m_2)g}{2m_1l_1l_2} \left\{(l_1+l_2)\pm \sqrt{(l_1+l_2)^2-\frac{4m_1}{m_1+m_2}l_1l_2}\right\}, \qquad h=\frac{l_2}{\gamma}. γ=2, l1=l2=l\gamma=2,\ l_1=l_2=l では ω±2=(g/l)(2±2)\omega_\pm^2=(g/l)(2\pm\sqrt2) であり、低振動数モードは同相、 高振動数モードは逆相である。

第2問 — 電場・双極子・誘電体球

一様帯電球のグラフ

内部電場は包む電荷が r3r^3 に比例する一方、ガウス面の面積が r2r^2 に比例するため E∝rE\propto r となる。外部では全電荷 QQ が中心にある場合と同じで E∝r−2E\propto r^{-2}。 グラフでは r=ar=a で折れ曲がるが、値は連続である。

双極子近似の符号

正電荷側を θ=0\theta=0 とすれば、軸上正方向では ϕd>0\phi_d>0 になる。 この符号確認をしておくと、距離の展開で R+R_+ と R−R_- を取り違えにくい。 電場の θ\theta 成分は Eθ=−1r∂ϕ∂θ E_\theta=-\frac1r\frac{\partial\phi}{\partial\theta} であるため、sin⁡θ\sin\theta の符号もここで決まる。

反電場の意味

誘電体球の内部で分極電荷が作る電場は、分極 P\bm P と逆向きである。 この電場は反電場と呼ばれ、球では係数が 1/31/3 になる。 平板や細長い針状物体では係数が変わるので、ここでは「球だから 1/31/3」と明記するのが大切である。

外部電場の書き方

外部電場を求めるとき、E0\bm E_0 を忘れて分極双極子の場だけを書いてしまう誤答が多い。 問題が「分極のみ」と指定している箇所では双極子の場だけを描き、全外部電場を求める箇所では 一様電場と双極子場の和を書く。答案ではこの二つを明確に分けると減点されにくい。

電場・双極子・誘電体球の途中式・最終答をPDFで見る

第3問 — 磁束を囲む荷電粒子

演算子の順序

(p−qA)2(\bm p-q\bm A)^2 では、p\bm p が右側の関数にも A\bm A にも作用する。 そのため (∇⋅A)ψ(\nabla\cdot\bm A)\psi の項が一度現れる。 最終的にこの問題の領域では ∇⋅A=0\nabla\cdot\bm A=0 だが、最初から消すと展開の根拠が曖昧になる。

磁場ゼロとベクトルポテンシャルゼロは違う

粒子のいる領域では ∇×A=0\nabla\times\bm A=0 である。しかし ∮A⋅dl=πr02B=Φ \oint \bm A\cdot d\bm l=\pi r_0^2B=\Phi はゼロとは限らない。環状領域は穴をもつため、A\bm A を単純に全域で勾配として消去できない。 ここが古典的なローレンツ力の議論との違いである。

量子数の役割

nn は狭い幅 Δr\Delta r 方向の箱型量子化を表し、mm は角度方向の周期性から来る量子数である。 Δr≪r0\Delta r\ll r_0 では径方向エネルギー ℏ2n2π2/(2MΔr2)\hbar^2n^2\pi^2/(2M\Delta r^2) が大きく、磁束の効果は主に角運動量項のずれとして現れる。

答案での考察

考察問題では「磁場がないのに力を受ける」と書くと不正確である。 定常状態のエネルギーが変わる理由は、波動関数の位相が閉曲線に沿って exp⁡(iqΦ/ℏ)\exp(iq\Phi/\hbar) の影響を受け、許される角度方向の波数がずれるためである。 局所的な力ではなく、領域のトポロジーと磁束による位相効果として説明する。

磁束を囲む荷電粒子の途中式・最終答をPDFで見る

第4問 — 平均場イジング模型

平均場方程式の意味

平均場近似では、各スピンは周囲のスピンそのものではなく平均値 mm が作る有効場を感じる。 ただし、その有効場によって計算した平均値が、最初に仮定した mm と一致しなければならない。 この自己無撞着条件が m=tanh⁡(zJm/kBT)m=\tanh(zJm/k_BT) である。

二重数えの補正

全エネルギーの係数 1/21/2 は重要である。一つのスピンに隣接スピンが zz 個あるので 単純には NzNz 個の結合を数えるが、結合 (i,j)(i,j) と (j,i)(j,i) を二回数えている。 そのため結合数は Nz/2Nz/2 になる。

エントロピーは組合せから出す

平均磁化 mm を固定すると、上向きと下向きの個数が決まる。 あとは NN 個の場所から上向きの場所を選ぶだけなので W(m)=N!N+!N−! W(m)=\frac{N!}{N_+!N_-!} となる。答案では N±=N(1±m)/2N_\pm=N(1\pm m)/2 を書いてからスターリング近似に入ると、 式の由来が明確になる。

安定性の判定

自己無撞着方程式の解が複数あるとき、全てが熱平衡状態とは限らない。 平衡状態は自由エネルギーの最小で決まる。T<TcT<T_c で m=0m=0 は方程式の解として残るが、 d2F/dm2<0d^2F/dm^2<0 なので最大であり、安定な解ではない。

臨界指数

転移点近傍で m0∝(Tc−T)1/2 m_0\propto (T_c-T)^{1/2} となる。平均場近似における秩序パラメータの臨界指数は 1/21/2 である。 展開では m0≠0m_0\ne0 の低温側解を扱っているので、途中で m0m_0 で割る理由を一言添えるとよい。

平均場イジング模型の途中式・最終答をPDFで見る

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

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

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

  • 2024年度(全4問)

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

  • 2023年度(全4問)

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