院試hub

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

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

全4問。熱・統計力学1問・力学1問。テーマタグは7件(調和振動子・固有値・固有ベクトル・分配関数)。2022年度と共通のテーマは調和振動子・ラグランジアン・シュレディンガー方程式。

最終更新:

収録6年度分の解答PDF:大阪大学 理学研究科 物理学専攻・宇宙地球科学専攻 物理学(¥3,200・紙面見本あり)

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

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

この4問の分野は熱・統計力学1問・力学1問です。

大問分野主題解説の小見出し最終答
第1問力学連成振動符号の検算 / 試験で書くべき点あり
第2問—磁気ミラー物理的な見方 / 符号の注意あり
第3問—スピンの時間発展位相の扱い / 歳差運動の意味あり
第4問熱・統計力学調和トラップ中の統計力学長さの意味 / 凝縮条件の読み方あり

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

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

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

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

大問数
2022年度 4問 → 2023年度 4問
2022年度のページを見る

第1問 — 連成振動

方針

この問題の要点は、中央のばねのポテンシャルを 12kAB(xA−xB)2\frac12k_{AB}(x_A-x_B)^2 と書けるかどうかである。等質量・等ばね定数の前半では xA+xBx_A+x_B と xA−xBx_A-x_B がそのまま基準振動になる。後半では、単独で同じ角振動数をもつという条件により、重心座標と相対座標の交差項が消える。

符号の検算

同相運動 xA=xBx_A=x_B では中央のばねは伸びないので、角振動数は k/m\sqrt{k/m} のままでなければならない。逆相運動 xA=−xBx_A=-x_B では中央のばねの伸びが 2xA2x_A となるため、復元力が強くなり角振動数が大きくなる。この物理的検算により、ωA−B\omega_{A-B} の +2kAB+2k_{AB} の符号を確認できる。

典型ミス

重心座標を単純平均 (xA+xB)/2(x_A+x_B)/2 としてしまうと、mA≠mBm_A\ne m_B の場合に運動エネルギーが対角化されない。また、xr=xA−xBx_r=x_A-x_B と定義したなら、中央のばねのエネルギーは符号によらず 12kABxr2\frac12k_{AB}x_r^2 である。ここで力の符号を直接追うより、ラグランジアンから導く方が誤りにくい。

試験で書くべき点

係数問題では、まず xG(t)x_G(t) と xr(t)x_r(t) を初期条件から決め、最後に xA,xBx_A,x_B へ戻す流れを答案に残すと部分点が取りやすい。特に初速度がどちらも 0 なので、sin⁡\sin 項の係数がすべて 0 になる理由を一行で示しておくとよい。

解答

平衡位置からの変位を右向きを正として測る。ばねの伸びは、左右の壁につながるものが xA,xBx_A,x_B、中央のばねが xB−xAx_B-x_A である。

  1. mA=mB=m, kA=kB=km_A=m_B=m,\ k_A=k_B=k の場合、運動方程式は mx¨A=−(k+kAB)xA+kABxB,mx¨B=kABxA−(k+kAB)xB.\begin{align} m\ddot{x}_A&=-(k+k_{AB})x_A+k_{AB}x_B,\\ m\ddot{x}_B&=k_{AB}x_A-(k+k_{AB})x_B. \end{align}
  2. 和と差をとると md2dt2(xA+xB)=−k(xA+xB),md2dt2(xA−xB)=−(k+2kAB)(xA−xB)\begin{align} m\frac{d^2}{dt^2}(x_A+x_B)&=-k(x_A+x_B),\\ m\frac{d^2}{dt^2}(x_A-x_B)&=-(k+2k_{AB})(x_A-x_B) \end{align} となる。したがって ωA+B=km,ωA−B=k+2kABm. \omega_{A+B}=\sqrt{\frac{k}{m}},\qquad \omega_{A-B}=\sqrt{\frac{k+2k_{AB}}{m}}.
  3. 中央のばねを外したとき、同じ角振動数 ω\omega で単振動する条件は kAmA=kBmB=ω2 \frac{k_A}{m_A}=\frac{k_B}{m_B}=\omega^2 である。
  4. 一般の場合のラグランジアンは L=12mAx˙A2+12mBx˙B2−12kAxA2−12kBxB2−12kAB(xA−xB)2. L=\frac12m_A\dot{x}_A^2+\frac12m_B\dot{x}_B^2 -\frac12k_Ax_A^2-\frac12k_Bx_B^2 -\frac12k_{AB}(x_A-x_B)^2 .
  5. xG=mAxA+mBxBmA+mB,xr=xA−xB,mG=mA+mB,mr=mAmBmA+mB x_G=\frac{m_Ax_A+m_Bx_B}{m_A+m_B},\quad x_r=x_A-x_B,\quad m_G=m_A+m_B,\quad m_r=\frac{m_Am_B}{m_A+m_B} とおくと xA=xG+mBmGxr,xB=xG−mAmGxr. x_A=x_G+\frac{m_B}{m_G}x_r,\qquad x_B=x_G-\frac{m_A}{m_G}x_r . 条件 kA=mAω2, kB=mBω2k_A=m_A\omega^2,\ k_B=m_B\omega^2 を用いると交差項が消え、 L=12mGx˙G2−12mGω2xG2+12mrx˙r2−12(mrω2+kAB)xr2. L=\frac12m_G\dot{x}_G^2-\frac12m_G\omega^2x_G^2 +\frac12m_r\dot{x}_r^2 -\frac12\left(m_r\omega^2+k_{AB}\right)x_r^2 .
  6. オイラー・ラグランジュ方程式より x¨G+ω2xG=0,x¨r+(ω2+kABmr)xr=0. \ddot{x}_G+\omega^2x_G=0,\qquad \ddot{x}_r+\left(\omega^2+\frac{k_{AB}}{m_r}\right)x_r=0.
  7. 相対運動の角振動数は ωr=ω2+kABmr. \omega_r=\sqrt{\omega^2+\frac{k_{AB}}{m_r}}.
  8. 初期条件から xG(0)=mAmGx0,x˙G(0)=0,xr(0)=x0,x˙r(0)=0 x_G(0)=\frac{m_A}{m_G}x_0,\quad \dot{x}_G(0)=0,\quad x_r(0)=x_0,\quad \dot{x}_r(0)=0 である。したがって xA(t)=mAmA+mBx0cos⁡ωt+mBmA+mBx0cos⁡ωrt,xB(t)=mAmA+mBx0cos⁡ωt−mAmA+mBx0cos⁡ωrt.\begin{align} x_A(t)&=\frac{m_A}{m_A+m_B}x_0\cos\omega t +\frac{m_B}{m_A+m_B}x_0\cos\omega_rt,\\ x_B(t)&=\frac{m_A}{m_A+m_B}x_0\cos\omega t -\frac{m_A}{m_A+m_B}x_0\cos\omega_rt . \end{align} よって係数は アイウエオカキク値mAx0mA+mB0mBx0mA+mB0mAx0mA+mB0−mAx0mA+mB0 \begin{array}{c|cccccccc} &\text{ア}&\text{イ}&\text{ウ}&\text{エ}&\text{オ}&\text{カ}&\text{キ}&\text{ク}\\ \hline \text{値} &\dfrac{m_Ax_0}{m_A+m_B}&0& \dfrac{m_Bx_0}{m_A+m_B}&0& \dfrac{m_Ax_0}{m_A+m_B}&0& -\dfrac{m_Ax_0}{m_A+m_B}&0 \end{array} である。
  9. mA≫mBm_A\gg m_B では mG≃mA, mr≃mBm_G\simeq m_A,\ m_r\simeq m_B であり、 ωr≃kB+kABmB=ω1+kABkB≃ω(1+kAB2kB). \omega_r\simeq \sqrt{\frac{k_B+k_{AB}}{m_B}} =\omega\sqrt{1+\frac{k_{AB}}{k_B}} \simeq \omega\left(1+\frac{k_{AB}}{2k_B}\right). したがって xA(t)≃x0cos⁡ωt x_A(t)\simeq x_0\cos\omega t で、重い物体 AA はほぼ固有振動を保つ。一方 xB(t)≃x0{cos⁡ωt−cos⁡ωrt}=2x0sin⁡(ωr+ω)t2sin⁡(ωr−ω)t2. x_B(t)\simeq x_0\{\cos\omega t-\cos\omega_rt\} =2x_0\sin\frac{(\omega_r+\omega)t}{2} \sin\frac{(\omega_r-\omega)t}{2}. 物体 BB はほぼ角振動数 ω\omega の速い振動をしながら、差周波 ωr−ω≃ωkAB/(2kB)\omega_r-\omega\simeq \omega k_{AB}/(2k_B) によるゆっくりした振幅変調を受ける。 包絡線の最大振幅はおよそ 2x02x_0 で、節を通過すると見かけの位相が反転する。

最終答

ωA+B=k/m\omega_{A+B}=\sqrt{k/m}、ωA−B=(k+2kAB)/m\omega_{A-B}=\sqrt{(k+2k_{AB})/m}。 一般の場合は重心運動と相対運動に分離し、 ωr=ω2+kAB/mr\omega_r=\sqrt{\omega^2+k_{AB}/m_r}。 初期条件付きの解は xA=x0{mAcos⁡ωt+mBcos⁡ωrt}/(mA+mB)x_A=x_0\{m_A\cos\omega t+m_B\cos\omega_rt\}/(m_A+m_B)、 xB=x0mA{cos⁡ωt−cos⁡ωrt}/(mA+mB)x_B=x_0m_A\{\cos\omega t-\cos\omega_rt\}/(m_A+m_B)。

第2問 — 磁気ミラー

方針

この問題は、ローレンツ力、∇⋅B=0\nabla\cdot\boldsymbol{B}=0、エネルギー保存を順につなぐ構成である。最初に一様磁場で回転の向きと符号を確定し、その符号を非一様磁場の計算にも引き継ぐことが重要である。

物理的な見方

dBz/dz>0dB_z/dz>0 では、進行方向に磁場が強くなる。保存量 vϕ2/Bzv_\phi^2/B_z のため、横方向速度は大きくなる。その分、全運動エネルギー保存により vzv_z が小さくなり、十分強い磁場に達すると vz=0v_z=0 となって反射する。これは磁気ミラーの基本機構である。

符号の注意

正電荷で Bz>0B_z>0 の場合、円運動は +z+z 方向から見ると時計回りである。したがって、問題で定められた正の ϕ\phi 方向とは逆向きになり vϕ<0v_\phi<0 となる。ここを落とすと mzm_z と FzF_z の符号が連鎖的に誤る。

検算

半径 r∝Bz−1/2r\propto B_z^{-1/2} は、強い磁場ほどラーマー半径が小さくなるという直観と一致する。また Φ=πr2Bz\Phi=\pi r^2B_z が一定であることは、断熱不変量が保たれていることの別表現であり、最後の設問の強い検算になる。

磁気ミラーの途中式・最終答をPDFで見る

第3問 — スピンの時間発展

方針

ハミルトニアンは n⋅σ\boldsymbol{n}\cdot\boldsymbol{\sigma} だけで決まる二準位系である。固有値はパウリ行列の代数から一行で決まり、時間発展はエネルギー固有状態に位相因子を付けるだけでよい。磁場の向きが時間変化する後半では、基底そのものの時間微分が a1,a2a_1,a_2 を混合する。

位相の扱い

エネルギーが −ℏω-\hbar\omega の状態には e+iωte^{+i\omega t}、+ℏω+\hbar\omega の状態には e−iωte^{-i\omega t} が掛かる。符号を逆にしても確率だけを見る設問では気づきにくいが、(a),(b) の虚数符号で誤りが現れる。

歳差運動の意味

⟨σz⟩=1−2sin⁡22θsin⁡2ωt\langle\sigma_z\rangle=1-2\sin^22\theta\sin^2\omega t は、初期スピンが有効磁場方向とずれているときだけ振動する。nx=0n_x=0 なら θ=0\theta=0 で振動せず、nz=0n_z=0 なら振幅が最大になる。これは古典的な歳差運動の量子版として理解できる。

断熱極限と急速極限

ゆっくり回す場合は、状態が瞬間固有状態に追随するため ∣↑⟩|\uparrow\rangle から ∣↓⟩|\downarrow\rangle に移る。速く回す場合は、状態ベクトルが実質的に動かないので ∣↑⟩|\uparrow\rangle のままである。選択肢問題では、全体位相は物理的に区別できないため無視してよい。

スピンの時間発展の途中式・最終答をPDFで見る

第4問 — 調和トラップ中の統計力学

方針

調和トラップでは「体積」に相当する量が普通の箱の体積ではなく、熱的な広がり lth3l_{\mathrm{th}}^3 で決まる。古典分配関数を Z1=(lth/λth)3Z_1=(l_{\mathrm{th}}/\lambda_{\mathrm{th}})^3 と見直すと、自由エネルギーが ρλth3\rho\lambda_{\mathrm{th}}^3 だけで書ける理由が明確になる。

長さの意味

lthl_{\mathrm{th}} はポテンシャルエネルギーが熱エネルギーと同程度になる空間的広がりであり、温度が高いほど大きく、トラップが強いほど小さい。一方 λth\lambda_{\mathrm{th}} は量子波束の熱的な波長で、温度が低いほど大きい。量子統計性は、この二つの比を通じて現れる。

凝縮条件の読み方

通常相では μ<eg\mu<e_g で、A=β(eg−μ)>0A=\beta(e_g-\mu)>0 である。温度を下げると N(βℏω)3N(\beta\hbar\omega)^3 が増加するが、励起状態が受け入れられる粒子数には I(0)I(0) による上限がある。その上限に達した後、化学ポテンシャルは ege_g を超えられず、基底状態の占有数だけが巨視的に増える。

典型ミス

N!N! を入れ忘れると古典自由エネルギーの 1−ln⁡(ρλth3)1-\ln(\rho\lambda_{\mathrm{th}}^3) が得られない。また、状態密度の係数 CC は縮退度 gs≃s2/2g_s\simeq s^2/2 から決まるため、単に3次元だから s2s^2 としてしまうと係数が2倍ずれる。

検算

C=1/2C=1/2 を用いると、臨界温度は Tc∝ωN1/3T_c\propto \omega N^{1/3} となる。粒子数を増やすほど、またトラップを強くするほど凝縮温度が上がるという物理的傾向と一致している。

調和トラップ中の統計力学の途中式・最終答をPDFで見る

阪大 物理学 院試 過去問の収録6年度

  • 2026年度(全4問)

    小振動と基準振動 / 電子ビームと電磁波 / 一次元量子井戸と散乱

  • 2025年度(全4問)

    解析力学・制限三体問題 / 電磁波・表皮効果 / 量子調和振動子

  • 2024年度(全4問)

    拘束条件付き力学 / 電磁場のエネルギー保存 / 環上粒子と摂動

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

    連成振動 / 磁気ミラー / スピンの時間発展

  • 2022年度(全4問)

    剛体振り子と連成振動 / 円環電荷と磁気双極子 / デルタ関数ポテンシャルの散乱

  • 2021年度(全4問)

    中心力・散乱・円軌道の安定性 / 回転帯電球殻の磁気双極子場 / 調和振動子・角運動量・磁場中の固有値