院試hub

京都大学 院試 過去問 解答例

京大 理学研究科 物理学・宇宙物理学専攻 専門科目(物理学) 2022年度 院試 解答例・解説

京都大学 理学研究科 物理学・宇宙物理学専攻 専門科目(物理学) 2022年度の院試 過去問について、設問ごとの解法方針・部分点の置き所を解説。全12問収録の解答・解説PDFと併用できます。問題本文は含みません。

最終更新:

設問ごとの解法方針・部分点の置き所を無料で公開しています。

完全な途中式・最終答は解答・解説PDFに収録しています。問題本文は含まれません。

1 — I-1 量子力学(2 次元等方調和振動子と一様磁場)

方針 — M^,N^\hat{M},\hat{N} は「左右回り励起」の生成子

2 次元等方調和振動子は, x,yx,y 方向の独立な調和振動子として量子化するか, あるいは「右回り(角運動量 -\hbar)」「左回り(角運動量 ++\hbar)」の独立な 2 モードに分解するかで等価. 後者のモードは b^+=a^x+ia^y2,b^=a^xia^y2 \hat{b}_{+} = \frac{\hat{a}_{x}+i\hat{a}_{y}}{\sqrt{2}},\qquad \hat{b}_{-} = \frac{\hat{a}_{x}-i\hat{a}_{y}}{\sqrt{2}} で定義され, M^=2b^,  N^=2b^+\hat{M}=\sqrt{2}\,\hat{b}_{-}^{\dagger},\;\hat{N}=\sqrt{2}\,\hat{b}_{+}^{\dagger}. すなわち M^\hat{M} は左回り励起(++\hbar の角運動量), N^\hat{N} は右回り励起(-\hbar). この見方をすると (3) の交換関係が「M^\hat{M}H^\hat{H}ω\hbar\omega 上げ, L^\hat{L}\hbar 上げる」という記述で透けて見える.

なぜ H^B\hat{H}_{B} の固有値が 2 つの独立振動子の和になるか

H^B\hat{H}_{B} は, 2 次元等方調和振動子の ωΩ\omega\to\Omega 置換 ++ ωLL^-\omega_{L}\hat{L} 補正. L^\hat{L}H^\hat{H} の対角化基底でも対角なので, 2 つは同時に対角化できる. 結局 H^B=(ΩωL)N^M+(Ω+ωL)N^N+Ω. \hat{H}_{B} = \hbar(\Omega-\omega_{L})\hat{N}_{M} + \hbar(\Omega+\omega_{L})\hat{N}_{N} + \hbar\Omega. ここで N^M=b^b^,N^N=b^+b^+\hat{N}_{M}=\hat{b}_{-}^{\dagger}\hat{b}_{-},\hat{N}_{N}=\hat{b}_{+}^{\dagger}\hat{b}_{+}. 強磁場極限 ωLω\omega_{L}\gg\omega では ΩωL0\Omega-\omega_{L}\to 0 となり, N^M\hat{N}_{M} モードがランダウ準位の縮退方向(各ランダウ準位内の状態を指定する)を担うことが分かる.

ベクトルポテンシャルの選び方とゲージ

(B) は対称ゲージ A=12B×r\vec{A}=\frac{1}{2}\vec{B}\times\vec{r} に他ならず, この選択は問題設定の中心力 VV の対称性(回転対称)と相性が良い. ランダウゲージ A=(By,0,0)\vec{A}=(-By,0,0) を選ぶと並進対称性が一方向で簡単化される代わりに L^\hat{L} との関係が見えにくくなる. 「L^\hat{L} を保存量にする」目的で対称ゲージを選んだ問題と理解できる.

典型ミス

  • (2) で零点エネルギー ω/2×2\hbar\omega/2 \times 2 を忘れる. x,yx,y 各方向に ω/2\hbar\omega/2 ずつあり合計 ω\hbar\omega.
  • (3) L^\hat{L} を昇降演算子で書く際の符号. L^=x^p^yy^p^x\hat{L}=\hat{x}\hat{p}_{y}-\hat{y}\hat{p}_{x} から始めて, 中点で a^xa^ya^ya^x\hat{a}_{x}\hat{a}_{y}^{\dagger}-\hat{a}_{y}\hat{a}_{x}^{\dagger} などが出るが, モードが独立なので a^xa^y=a^ya^x\hat{a}_{x}\hat{a}_{y}^{\dagger}=\hat{a}_{y}^{\dagger}\hat{a}_{x}. 慎重に.
  • (5) で qBL/(2m)-qB\,L/(2m) ではなく +qBL/(2m)+qB\,L/(2m) になる. ypxxpy=Ly\,p_{x}-x\,p_{y} = -L の符号を確認すること(L=xpyypxL=xp_{y}-yp_{x} の定義).
  • (7) で \ell の最大値を DD ではなく D1D-1 と勘違い. {N,,N}\ell\in\{-N,\dots,N\}N=D1N=D-1.
  • (8) で「(b) は B=0B=0 で水平」と書く. (a) のみ水平で, (b) は線形傾き n(Dn)q/(2m)-n(D-n)\hbar q/(2m) が残る. n(Dn)=0n(D-n)=0 となるのは n=0n=0n=Dn=D の場合だけ.

検算 — 角運動量保存と縮退度

(3) で [H^B=0,L^]=0[\hat{H}_{B=0},\hat{L}]=0 になることは, ハミルトニアンが回転対称(中心力ポテンシャル + ベクトルポテンシャル無し)である事実の言い換え. 一方, L^\hat{L} が保存量であることから, エネルギー準位の縮退は L^\hat{L} の固有値の数で説明できて D=N+1D=N+1.

(7) を n=Dn=D, n=0n=0 で評価すると Etotn=D=D2ΩE_{\text{tot}}|_{n=D}=D^{2}\hbar\Omega, Etotn=0=0E_{\text{tot}}|_{n=0}=0. 後者は粒子無しなのでゼロが正しい. 前者は DD 個の状態すべてのエネルギーの和で, \ell の総和がゼロ(=N,step 2N=0\sum_{\ell=-N,\,\text{step }2}^{N}\ell=0)になるので ωL=0-\omega_{L}\hbar\sum\ell=0 から DΩD=D2ΩD\cdot \hbar\Omega\cdot D=D^{2}\hbar\Omega も整合する.

背景 — Fock-Darwin スペクトルとランダウ準位

このハミルトニアンは Fock-Darwin スペクトルとして知られ, 量子ドットの低エネルギー物理の標準モデル. 強磁場極限でランダウ準位(ω0\omega\to 0 または ωLω\omega_{L}\gg\omega)に, 弱磁場極限で 2 次元調和振動子(ゼーマン分裂のみ)に連続的に補間する. (8) のグラフは, 半導体量子ドットでの磁気光学スペクトルの実測パターン(Aufbau 原理 + Hund 則の競合)の出発点.

完全な解答(途中式・最終答)はPDFに収録

2 — I-2A 電磁気学(電気双極子・四重極子)

方針 — 多重極展開を「打ち消し合いの順序」で見る

電荷分布を遠方から見るときの 1/r1/r 展開:

  • 1/r1/r(単極子): 全電荷.
  • 1/r21/r^{2}(双極子): 電気双極子モーメント p=qiri\vec{p}=\sum q_{i}\vec{r}_{i}.
  • 1/r31/r^{3}(四重極子): 電気四重極子テンソル.

本問の ϕ1\phi_{1}(2 電荷の双極子)は全電荷 00 なので単極子項が消え, 双極子項が先頭. ϕ3\phi_{3} は全電荷 00 かつ双極子モーメント 00 なので, 四重極子項が先頭. ϕ4\phi_{4} は全電荷 2q02q\neq 0 なので単極子項が先頭になる. これを見抜けると, 計算前から「先頭の落ち方」が予測できる.

(1+x)1/2(1+x)^{-1/2} 展開の使いどころ

1/r±a1/|\vec{r}\pm \vec{a}|a\vec{a} の小さい場合に展開するときは, 中身を完全平方の形 r2(1±2r^a/r+)r^{2}(1\pm 2\hat{r}\cdot\vec{a}/r+\dots) に整理してから (A) を使う. このとき注意すべきは「次数」:

  • a/r\vec{a}/r1 次 のみ取れば双極子レベル(問 (2)).
  • 2 次 まで取れば四重極子レベル(問 (4),(6)).

(1+x)1/21x/2+(3/8)x2(1+x)^{-1/2}\approx 1-x/2+(3/8)x^{2}xx そのものに ε2\varepsilon^{2} の寄与もある(問 (4) の ε2\varepsilon^{2} 項)ので, x/2-x/2 から出る ε2\varepsilon^{2} も忘れないこと.

典型ミス

  • (1) の概略図で原点でゼロにし忘れる. ϕ1(0,0,0)=q/(4πε0)(2/d2/d)=0\phi_{1}(0,0,0)=q/(4\pi\varepsilon_{0})\cdot(2/d - 2/d)=0.
  • (2) で 1/r1/r の項を残す. 全電荷 00 なので 1/r1/r は完全に打ち消し合う. 計算結果が A/rA/r の形になるなら符号ミス.
  • (3) で電場を p/r3-\vec{p}/r^{3} だけにしてしまう. 微分の連鎖則で (pr)(1/r3)(\vec{p}\cdot\vec{r})\nabla(1/r^{3}) から出る 3(pr)r/r53(\vec{p}\cdot\vec{r})\vec{r}/r^{5} の項を忘れない.
  • (4) で ε2\varepsilon^{2} 展開の係数 3/83/81/21/2 にする. (1+x)1/2(1+x)^{-1/2} の 2 次係数は (1/2)(3/2)2=38\frac{(-1/2)(-3/2)}{2}=\frac{3}{8}.
  • (5) の EθE_{\theta}rr の冪を間違える. ϕ31/r3\phi_{3}\propto 1/r^{3} なのに Eθ=(1/r)ϕ/θE_{\theta}=-(1/r)\partial\phi/\partial\theta なので 1/r41/r^{4} になる(1/r31/r=1/r41/r^{3}\cdot 1/r=1/r^{4}).
  • (6) で ϕ4\phi_{4} の先頭を 00 と書く. 全電荷 2q02q\neq 0 なので単極子項 2q/r2q/r が先頭.

検算 — 双極子モーメントのチェック

(2) の結果は, 双極子モーメント定義から導いた一般公式 ϕ(r)=pr/(4πε0r3)\phi(\vec{r})=\vec{p}\cdot\vec{r}/(4\pi\varepsilon_{0}r^{3}) そのもの. 双極子の方向 p^=e^z\hat{p}=\hat{e}_{z} で確認すると pr=qdz\vec{p}\cdot\vec{r}=qd\,z, よって ϕ=qdz/(4πε0r3)\phi=qd\,z/(4\pi\varepsilon_{0}r^{3}) が我々の式と一致.

(4) の極限を確認. θ=0\theta=0(極軸方向)で 3cos2θ1=23\cos^{2}\theta-1=2, ϕ3=2qd2/(4πε0r3)>0\phi_{3}=2qd^{2}/(4\pi\varepsilon_{0}r^{3})>0. θ=π/2\theta=\pi/2(赤道面)で 3cos2θ1=13\cos^{2}\theta-1=-1, ϕ3=qd2/(4πε0r3)<0\phi_{3}=-qd^{2}/(4\pi\varepsilon_{0}r^{3})<0. 物理的には, 極軸近傍は +q+q に近く, 赤道面は 2q-2q の影響が相対的に大きいので符号は妥当.

背景 — Legendre 多項式と多重極子

1/ra1/|\vec{r}-\vec{a}| の展開係数は本質的にルジャンドル多項式 P(cosθ)P_{\ell}(\cos\theta) で, 一般に 1ra==0ar+1P(cosθ)(r>a). \frac{1}{|\vec{r}-\vec{a}|} = \sum_{\ell=0}^{\infty}\frac{a^{\ell}}{r^{\ell+1}}P_{\ell}(\cos\theta)\quad(r>a). 四重極子の角度依存 3cos2θ1=2P2(cosθ)3\cos^{2}\theta-1=2P_{2}(\cos\theta) はこの展開の =2\ell=2 成分. (5) の EθE_{\theta} に現れる sin(2θ)\sin(2\theta) も, P2(cosθ)cosθP_{2}'(\cos\theta)\propto\cos\theta から自然に出てくる結果.

完全な解答(途中式・最終答)はPDFに収録

3 — I-2B 物理数学(1 次元波動方程式の初期値問題)

方針 — d'Alembert 解の物理的解釈

1 次元波動方程式 (A) の特性曲線(波が進む直線)は x±ct=constx\pm ct=\text{const} で, この線に沿って位相情報が保存される. 変数変換 ξ=x+ct,η=xct\xi=x+ct,\eta=x-ct はこの特性曲線を新しい座標軸に取ったもので, 方程式が uξη=0u_{\xi\eta}=0 という最も単純な形になる. 解は「左進行波 ψ(ξ)=ψ(x+ct)\psi(\xi)=\psi(x+ct)」と「右進行波 ϕ(η)=ϕ(xct)\phi(\eta)=\phi(x-ct)」の重ね合わせ.

初期速度ゼロ条件と「半分ずつ分裂」

初期速度がゼロのとき, ψ=ϕ\psi'=\phi' から ψ=ϕ\psi=\phi (++ 定数)となり, 初期波形が左右にちょうど半分ずつ分裂する形で進む. これが (6) の式 u(x,t)=12[u0(x+ct)+u0(xct)]u(x,t)=\frac{1}{2}[u_{0}(x+ct)+u_{0}(x-ct)](u0u_{0} は初期波形)の意味. 一般の(初期速度 0\neq 0)場合は d'Alembert の公式 u(x,t)=12[u0(x+ct)+u0(xct)]+12cxctx+ctv0(ξ)dξ u(x,t) = \frac{1}{2}\bigl[u_{0}(x+ct)+u_{0}(x-ct)\bigr] + \frac{1}{2c}\int_{x-ct}^{x+ct}v_{0}(\xi)\,d\xi で表され, 第 2 項が初期速度 v0v_{0} の寄与を担う.

典型ミス

  • (2) で /t=c/ξ+c/η\partial/\partial t = c\partial/\partial\xi + c\partial/\partial\eta と符号を取り違える. ξ=x+ct\xi=x+ct なので ξ/t=+c\partial\xi/\partial t = +c, η=xct\eta=x-ct なので η/t=c\partial\eta/\partial t = -c. したがって /t=c/ξc/η\partial/\partial t=c\partial/\partial\xi - c\partial/\partial\eta.
  • (4) で u/t\partial u/\partial tψ(x+ct)\psi'(x+ct) ではなく ψ(x+ct)c\psi(x+ct)\cdot c と書く. 連鎖則で外側微分 cψc\psi'.
  • (5)(6) の積分定数を最終解に残す. C1,C2C_{1},C_{2}C1+C2=0C_{1}+C_{2}=0 かつ uuψ+ϕ\psi+\phi の形でしか効かないので消える. 任意定数を残すと初期条件と矛盾(無限大に発散しないため exp\exp で抑えても定数のみ残るのは不合理).
  • (7) で t=3/ct=3/c のピーク位置を ±c3\pm c\cdot 3 などにする. x=±ctx=\pm ctt=3/ct=3/c なら x=±3x=\pm 3. 単位の整合性に注意.

検算 — 波動方程式の直接代入

(6) の解が (A) を満たすかを確認: 2ut2=12[(2c)2(2(x+ct)21)exp((x+ct)2)+], \frac{\partial^{2}u}{\partial t^{2}} = \frac{1}{2}\Bigl[(-2c)^{2}\bigl(2(x+ct)^{2}-1\bigr)\exp(-(x+ct)^{2}) + \cdots\Bigr], このように c22u/x2c^{2}\partial^{2}u/\partial x^{2} と一致する. 簡単にチェックするには「左進行と右進行のそれぞれが波動方程式を満たす」ことだけ確かめれば線形性で十分.

背景 — Cauchy 問題の well-posedness

本問の初期値問題は, u(x,0)u(x,0)tu(x,0)\partial_{t}u(x,0) の 2 つの初期データから一意な解を構成する Cauchy 問題. 1 次元波動方程式は典型的な双曲型方程式で, 任意の連続な初期値に対して d'Alembert 解が古典解を与える. これが熱方程式や波動方程式 + 散逸など他の PDE の解析の出発点.

完全な解答(途中式・最終答)はPDFに収録

4 — I-3A 力学(回転系・コリオリ力・ラーモア定理)

方針 — 「ベクトルの時間微分公式」を 2 度使う

(A) はベクトル一般に対する公式で, 回転系で見た時間微分 d/dtd'/dt と慣性系の d/dtd/dt の関係. (1) でやるべきことは「r\vec{r} への適用 → dr/dt=dr/dt+ω×rd\vec{r}/dt = d'\vec{r}/dt + \vec{\omega}\times\vec{r}」を出発点に, さらに「vdr/dt\vec{v}\equiv d'\vec{r}/dt という ベクトル に同じ公式を適用する」こと. 公式を 2 段階で用いる感覚さえ掴めば, 加速度が「真の加速度 ++ 2 倍コリオリ ++ 遠心」の 3 項に分解されるのは機械的.

Larmor 定理の意味

(3) の結論は ラーモアの定理: 弱い磁場のもとで, ラーモア角速度 ωL=qB/(2m)\vec{\omega}_{L}=-q\vec{B}/(2m) で回転する系を取ると, 磁場が O(B2)O(B^{2}) までは消え去って見える. これは原子物理での「ゼーマン効果は磁場 1 次に比例」「反磁性は磁場 2 次」という階層構造の力学的根拠. 第 1 問 (5)(6) の Fock-Darwin スペクトルでも, ωLL-\omega_{L}L の項としてこの定理が顔を出している.

典型ミス

  • (1) の係数 2 を忘れる. コリオリ項に 2ω×v2\vec{\omega}\times\vec{v}\,' という 22 が現れるのは, dr/dtd'\vec{r}/dt の時間微分から 1 つ, さらに (A) を 2 度目に当てる際にもう 1 つ, 同じ ω×v\vec{\omega}\times\vec{v}\,' が出るため.
  • (2) で v=v\vec{v}=\vec{v}\,' と取り違える. Lorentz 力は 慣性系での 速度 v\vec{v} に依存. 回転系の速度 v\vec{v}\,' ではない.
  • (3) で ω=qB/mk^\vec{\omega}=qB/m\hat{k} と書く(2 で割らない). これはサイクロトロン振動数で, Larmor 振動数の倍. Larmor 振動数は qB/(2m)qB/(2m).
  • (3) で符号を逆にする. q>0q>0, B\vec{B}+k^+\hat{k} 方向のとき, ラーモア角速度は k^-\hat{k} 方向(粒子の磁気力による旋回方向と逆向きの 半分). 物理的に, 磁場による旋回の半分速度で逆方向に回ると, 系から見て磁場効果が消える.

検算 — 単位と次元

ωL=qB/(2m)\vec{\omega}_{L}=-q\vec{B}/(2m) の次元: [q][B]/[m]=(C)(T)/(kg)=(Ckg/(Cs))/kg=1/s[q][B]/[m] = (\text{C})(\text{T})/(\text{kg}) = (\text{C}\cdot\text{kg/(C}\cdot\text{s)})/\text{kg} = 1/\text{s}. 確かに角速度の次元(rad/s).

(1) の関係式 ainertial=a+2ω×v+ω×(ω×r)\vec{a}_{\text{inertial}} = \vec{a}\,' + 2\vec{\omega}\times\vec{v}\,' + \vec{\omega}\times(\vec{\omega}\times\vec{r}) を地球の自転(東京で北を見て)に適用すると, 地球表面で右(東向き)に投げた物体は北半球で右側にずれる. これがコリオリ力の有名な台風や貿易風の説明.

背景 — Larmor 定理と磁気力学

ラーモア定理は弱磁場での近似的な等価性で, 中心力ポテンシャル下での電子の運動を, 磁場無し問題に「ラーモア回転」を加えただけのもので近似できる. これが原子の磁気モーメントの計算で, 「軌道角運動量と磁場が結合して μB-\vec{\mu}\cdot\vec{B} のゼーマン項を生む」ことの古典的根拠. 反磁性(O(B2)O(B^{2}))を生む遠心力ポテンシャルの ω×(ω×r)\vec{\omega}\times(\vec{\omega}\times\vec{r}) 項からは, 磁場 2 次の補正として Langevin の反磁性公式が導かれる.

完全な解答(途中式・最終答)はPDFに収録

5 — I-3B 熱力学(熱力学的考察によるシュテファン・ボルツマン則)

方針 — 熱力学だけで温度依存性を絞り込む

光子気体の状態方程式について必要な情報は次の 2 点だけ.

  1. 圧力とエネルギー密度の比例関係 p=u~/3p=\tilde{u}/3(電磁場の Maxwell 応力テンソルから).
  2. u~=u~(T)\tilde{u}=\tilde{u}(T) (VV に依存しない. 強度パラメータが温度のみ).

これと熱力学関係式 (U/V)T=T(p/T)Vp(\partial U/\partial V)_{T}=T(\partial p/\partial T)_{V}-p を組み合わせると, u~\tilde{u} の温度依存性が「変数分離型 ODE」 4u~=Tdu~/dT4\tilde{u} = T\,d\tilde{u}/dT で決まり, 解 u~T4\tilde{u}\propto T^{4} が出る. 微視的物理(光子の 2 偏光自由度 ++ ボーズ統計 ++ 分散 ω=ck\omega=ck)無しで T4T^{4} 依存性が決まるのが熱力学の威力.

なぜ p=u~/3p=\tilde{u}/3 となるか

これは「電磁波が真空中で速さ cc で進む」「3 次元等方的に運動量を運ぶ」ことと相対論から導かれる. ざっくり言うと, 光子気体の応力テンソル TijT^{ij} の対角成分が u~/3\tilde{u}/3(i=ji=j, 各方向に等しく等方圧). 物質粒子(非相対論的理想気体)では p=23u~p=\frac{2}{3}\tilde{u}(粒子 1 個あたりエネルギー 12mv2\frac{1}{2}m v^{2}, 圧力 13nmv2\frac{1}{3}n m \langle v^{2}\rangle から p=23n12mv2p=\frac{2}{3}n\langle\frac{1}{2}m v^{2}\rangle). 一般に超相対論極限で p=u~/3p=\tilde{u}/3 になる.

典型ミス

  • (1) で dUdU に化学ポテンシャル項 μdN\mu\,dN を含めない. 光子気体では μ=0\mu=0(光子は数が保存しない)から本問では問題にしないが, 一般の系では含める必要がある.
  • (2) で Maxwell 関係を間違える. FF から導く式は (S/V)T=(p/T)V(\partial S/\partial V)_{T}=(\partial p/\partial T)_{V}. GGHH から導く別の Maxwell 式と取り違えない.
  • (3) で u~=u~(T,V)\tilde{u}=\tilde{u}(T,V) と書いて VV 依存性を残す. 光子気体の特殊性として u~\tilde{u} は強度量で, 体積を変えても変わらない. これは「化学ポテンシャル μ=0\mu=0 かつ 1 種類の素粒子のみ」という条件から従う.
  • 微分方程式の符号. 4u~=Tdu~/dT4\tilde{u}=T\,d\tilde{u}/dT から du~/u~=4dT/Td\tilde{u}/\tilde{u}=4\,dT/TT4T^{4}. 「3 と 4」を混同しないこと.

検算 — 古典理想気体での同じ手続き

古典理想気体は pV=NkBTpV=Nk_{B}T, すなわち p=NkBT/Vp=Nk_{B}T/V. (p/T)V=NkB/V=p/T(\partial p/\partial T)_{V}=Nk_{B}/V=p/T. (UV)T=TpTp=0. \Bigl(\frac{\partial U}{\partial V}\Bigr)_{T} = T\cdot\frac{p}{T} - p = 0. 理想気体の内部エネルギーが VV に依存せず温度のみで決まる, という有名な事実が同じ Maxwell 関係から従う. (3) と同じ手順を「光子気体に対する状態方程式 p=u~/3p=\tilde{u}/3 で動かす」のが本問の構図.

背景 — Stefan-Boltzmann 定数とランプの色温度

u~=aT4\tilde{u}=aT^{4} から「単位面積あたりの放射エネルギー流」 j=σT4j=\sigma T^{4} (σ=ac/4\sigma=ac/4) が出る. これは恒星のスペクトルや黒体放射, 太陽電池の効率(Shockley-Queisser 限界)の計算で繰り返し使われる基本式. (1)-(3) のように 1 ページで導けるのが熱力学の力.

完全な解答(途中式・最終答)はPDFに収録

6 — I-3C 物理数学(ベクトル解析・連立 ODE・無限積分)

方針 — 「テンソル成分」と「ベクトル形」を行き来する

(1) は典型的な「レビ・チビタ記号 εijk\varepsilon_{ijk} + 縮約公式 εε=δδδδ\varepsilon\varepsilon=\delta\delta-\delta\delta」のパターン. ベクトル恒等式は丸暗記より, 成分で 3 行で導けるようにしておくと圧倒的にロバスト. 中央の縮約 εijkεklm=δilδjmδimδjl\varepsilon_{ijk}\varepsilon_{klm} = \delta_{il}\delta_{jm}-\delta_{im}\delta_{jl} さえ覚えればよい.

連立 ODE は固有値問題に帰着

線形定数係数の連立 ODE は, 行列 MM を作ってその固有値・固有ベクトルを求めれば終わり. 解は Cieλixvi\sum C_{i}\,e^{\lambda_{i}x}\vec{v}_{i}. 本問は 2×22\times 2 行列で固有値が異なる(対角化可能)ケース. 重根や複素固有値があれば Jordan 標準形や sin/cos\sin/\cos への変換が必要だが, 本問は素直.

典型ミス

  • (1) で縮約 εε\varepsilon\varepsilon の符号. εijkεklm=δilδjmδimδjl\varepsilon_{ijk}\varepsilon_{klm}=\delta_{il}\delta_{jm}-\delta_{im}\delta_{jl} で, 添字の対応は「1 番目同士・2 番目同士」.
  • (2) で div(A×B)=ArotBBrotA\mathrm{div}(\vec{A}\times\vec{B}) = \vec{A}\cdot\mathrm{rot}\vec{B} - \vec{B}\cdot\mathrm{rot}\vec{A} と符号を逆にする. 正しくは BrotAArotB\vec{B}\cdot\mathrm{rot}\vec{A} - \vec{A}\cdot\mathrm{rot}\vec{B}. これは多くの教科書で間違えやすい点.
  • (3) で MM の取り方を (z,y)=(z',y')=\dots にして固有値・固有ベクトルを誤る. y,zy',z' の順を一貫させること.
  • (4) で項別積分の正当化を忘れる. テイラー級数 log(1x)=xn/n\log(1-x)=-\sum x^{n}/nx<1|x|<1 でしか収束しないが, x=1x=1 までの広義積分が収束(被積分関数 log(1x)/x\log(1-x)/xx1x\to 1 で可積分: log(1x)\log(1-x) が積分可能な特異点)するので Abel の定理(あるいは単調収束)で項別積分が正当化される.

検算 — Spence 関数(ジロサインドラー積分)

Li2(z)=0zlog(1t)tdt=n=1znn2\mathrm{Li}_{2}(z) = -\int_{0}^{z}\frac{\log(1-t)}{t}dt = \sum_{n=1}^{\infty}\frac{z^{n}}{n^{2}} は二重対数関数(dilogarithm). 本問の積分は Li2(1)=π2/6-\mathrm{Li}_{2}(1) = -\pi^{2}/6 そのもの. これは Euler によって 1735 年に解決された有名な「バーゼル問題」 1/n2=π2/6\sum 1/n^{2}=\pi^{2}/6 の積分表示.

背景 — Helmholtz 分解と (B) 式

(B) は流体力学の Bernoulli 則の出発点. 非粘性流体の Euler 方程式 vt+(v)v=pρ \frac{\partial\vec{v}}{\partial t} + (\vec{v}\cdot\nabla)\vec{v} = -\frac{\nabla p}{\rho} で, (v)v=(v2/2)v×(×v)(\vec{v}\cdot\nabla)\vec{v} = \nabla(\vec{v}^{2}/2) - \vec{v}\times(\nabla\times\vec{v}) と (B) の符号反転を使うと「定常で渦無しの流れでは p+ρv2/2=p+\rho v^{2}/2 = 一定」が出る. (1) は単なる代数的恒等式に見えて, 流体力学・磁場の MHD まで波及する基本ツール.

完全な解答(途中式・最終答)はPDFに収録

7 — I-3D 量子力学(角運動量の合成)

方針 — 「最高ウェイト + 降下」で Clebsch-Gordan 展開を作る

角運動量合成の Clebsch-Gordan 係数を導くアルゴリズム:

  1. 合成系で最高の Sz=SS_{z}=S を持つ状態は 1 通りに決まる(両側で最高ウェイトを取った積状態).
  2. これに S^\hat{S}^{-} を順に作用させて, 同じ SS 多重項の他の状態を得る.
  3. 異なる S<SS'<S の最高ウェイト状態は, 与えられた SzS_{z} の部分空間内で既に得た S,Sz|S,S_{z}\rangle に直交するように構成する.
  4. ステップ 2 と 3 を繰り返して, すべての S,Sz|S,S_{z}\rangle を生成.

本問 (1)–(3) はまさにこの手順. (4) では (i) のステップで, スピン 11/21\otimes 1/2 の最高 S=3/2,Sz=3/2S=3/2,S_{z}=3/2 から S^\hat{S}^{-} を 1 回適用して 3/2,1/2|3/2,1/2\rangle を出している.

なぜ 0,0|0,0\rangle は反対称か

(3) の 0,0|0,0\rangle が「2 粒子交換について反対称」なのは, スピン==00 の合成系は シングレット(回転対称な唯一の状態)であり, スピン部分の対称性は粒子交換について反対称になる. これは 2 電子の場合「スピン反対称 → 軌道対称(空間に対称)」となり, 化学結合の covalent 結合(分子軌道のシグマ)と整合する.

逆に S=1S=1 の トリプレット 1,m|1,m\rangle (m=1,0,1m=-1,0,1) は対称な部分空間, 軌道反対称となるので Hund 則(同じスピンで電子が並ぶ方が安定)を理解する基礎.

典型ミス

  • (2) で規格化定数 1/21/\sqrt{2}1/21/2 にする. 2 つの規格化された直交状態の和の規格化は 1/21/\sqrt{2}. 「2 つあるから 1/2 で割る」と勘違いしないこと.
  • (3) で 0,0|0,0\rangle の符号を逆にする. 規約的には 1/2,1/2;1/2,1/20,0\langle 1/2,1/2;1/2,-1/2|0,0\rangle を正にとる(Condon–Shortley 規約)のが標準. 全体に絶対値の符号は物理に効かないが, テキスト間で食い違いが出やすい.
  • (4) で S^S,Sz\hat{S}^{-}|S,S_{z}\rangle の係数 (S+Sz)(SSz+1)\sqrt{(S+S_{z})(S-S_{z}+1)} を取り違える. 上昇演算子は (SSz)(S+Sz+1)\sqrt{(S-S_{z})(S+S_{z}+1)} で対比的. 公式自体を覚えるか, S^2S^z2S^z\hat{S}^{2}-\hat{S}_{z}^{2}-\hbar\hat{S}_{z} の対角化から導けるようにしておく.
  • (4) で全状態 3/2,1/2|3/2,1/2\rangle1,01/2,1/2|1,0\rangle|1/2,1/2\rangle だけで終わらせる. 1,11/2,1/2|1,1\rangle|1/2,-1/2\rangle も足し合わせるのが正しい(両者とも Sz=1/2S_{z}=1/2 の部分空間に居る).

検算 — 規格化と直交性

(2) と (3) を内積で確認: 1,01,0=(1/2)(1+1)=1\langle 1,0|1,0\rangle = (1/2)(1+1) = 1 ✓. 1,00,0=(1/2)(11)=0\langle 1,0|0,0\rangle = (1/2)(1-1) = 0 ✓ (符号差で打ち消し). 0,00,0=(1/2)(1+1)=1\langle 0,0|0,0\rangle = (1/2)(1+1) = 1 ✓.

(4) も同様に 3/2,1/22=2/3+1/3=1|\,|3/2,1/2\rangle\,|^{2} = 2/3 + 1/3 = 1 で規格化されている.

背景 — Clebsch-Gordan 表と既約分解

スピン j1j2j_{1}\otimes j_{2} の合成は次のように既約分解される: j1j2=J=j1j2j1+j2J. \mathbf{j_{1}}\otimes\mathbf{j_{2}} = \bigoplus_{J=|j_{1}-j_{2}|}^{j_{1}+j_{2}}\mathbf{J}. 本問の (3) は 1/21/2=011/2\otimes 1/2 = 0\oplus 1 の構成, (4) は 11/2=1/23/21\otimes 1/2 = 1/2\oplus 3/2 の構成. Clebsch-Gordan 係数表(3j3j-symbol や 6j6j-symbol)はこれを系統的にまとめたもの. 多体系のスピンや原子の項記号(例 2S+1LJ^{2S+1}L_{J})を扱う場合に必須の道具.

完全な解答(途中式・最終答)はPDFに収録

8

この設問は図表を含むため、解説はPDF版でご確認いただけます。

完全な解答(途中式・最終答)はPDFに収録

9

この設問は図表を含むため、解説はPDF版でご確認いただけます。

完全な解答(途中式・最終答)はPDFに収録

10

この設問は図表を含むため、解説はPDF版でご確認いただけます。

完全な解答(途中式・最終答)はPDFに収録

11

この設問は図表を含むため、解説はPDF版でご確認いただけます。

完全な解答(途中式・最終答)はPDFに収録

12

この設問は図表を含むため、解説はPDF版でご確認いただけます。

完全な解答(途中式・最終答)はPDFに収録

京都大学 専門科目(物理学) — 他の年度