跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2017年度 問題I-3-ABC

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

I-3A(物理数学)​

拡散方程式

∂g∂t=∂2g∂x2,−∞<x<∞,t>0\frac{\partial g}{\partial t}=\frac{\partial^2g}{\partial x^2},\qquad -\infty<x<\infty,\quad t>0

を、g(x,0)=δ(x)g(x,0)=\delta(x) の条件で解け。必要なら

∫0∞e−ax2cos⁡(bx) dx=12πa e−b2/(4a),a>0\int_0^\infty e^{-ax^2}\cos(bx)\,dx=\frac12\sqrt{\frac\pi a}\,e^{-b^2/(4a)},\qquad a>0

を用いてよい。

I-3B(統計力学)​

非縮退のエネルギー準位 0,Δ0,\Delta を持つ独立な粒子 NN 個が、温度 TT で熱平衡にある。ボルツマン定数を kk とする。

  1. 分配関数 Z(T)Z(T) を書け。
  2. Z(T)Z(T) からヘルムホルツ自由エネルギー F(T)F(T) を求めよ。
  3. F(T)F(T) からエントロピー S(T)S(T) を求めよ。
  4. T→∞T\to\infty での SS とその物理的意味を、式を含めて簡潔に説明せよ。
  5. SS を T,FT,F、内部エネルギー UU だけで表せ。
  6. U(T)U(T) を求めよ。

I-3C(力学)​

質量 mm の質点が球面 x2+y2+z2=R2x^2+y^2+z^2=R^2 上に拘束され、球面を離れず滑らかに運動する。zz 軸を鉛直上向きとし、重力加速度は −gez-g\mathbf e_z とする。

  1. 南極 (0,0,−R)(0,0,-R) 付近で y=0y=0 を保つ微小振動の周期を求めよ。
  2. 高さ z=−z0z=-z_0(0<z0<R0<z_0<R)の水平面で円運動する速さ v0v_0 を求めよ。
  3. (R,0,0)(R,0,0) で yy 方向の初速度 v1v_1 を与えたとき、その後の高さの最小値 z1z_1 を求めよ。

题目描述​

I-3A: 在全实轴上解扩散方程 gt=gxxg_t=g_{xx},初值为 g(x,0)=δ(x)g(x,0)=\delta(x)。可使用题给高斯余弦积分。

I-3B: NN 个独立粒子各有两个非简并能级 0,Δ0,\Delta,温度为 TT,玻尔兹曼常数为 kk。依次求:(1) 分配函数;(2) 自由能;(3) 熵;(4) 高温熵极限及物理含义;(5) 用 T,F,UT,F,U 表示熵;(6) 内能。

I-3C: 质量 mm 的质点被约束在半径 RR 的球面上光滑运动且不脱离球面,zz 轴向上,重力向下。(1) 求保持 y=0y=0、在南极附近微振动的周期;(2) 求在 z=−z0z=-z_0(0<z0<R0<z_0<R)水平面作圆周运动的速率;(3) 从 (R,0,0)(R,0,0) 沿 yy 方向以速度 v1v_1 出发,求最低高度 z1z_1。

Kai​

I-3A​

g~(k,t)=∫g(x,t)e−ikx dx\widetilde g(k,t)=\int g(x,t)e^{-ikx}\,dx とすると

∂tg~=−k2g~,g~(k,0)=1,g~=e−k2t.\partial_t\widetilde g=-k^2\widetilde g,\quad \widetilde g(k,0)=1, \quad \widetilde g=e^{-k^2t}.

逆変換により

g(x,t)=14πtexp⁡ ⁣(−x24t),t>0.\boxed{g(x,t)=\frac1{\sqrt{4\pi t}}\exp\!\left(-\frac{x^2}{4t}\right),\qquad t>0.}

積分は常に 11 であり、t↓0t\downarrow0 でデルタ関数に分布の意味で収束する。

I-3B​

(1)​

粒子一個の分配関数は 1+e−Δ/(kT)1+e^{-\Delta/(kT)} だから

Z=(1+e−Δ/(kT))N.\boxed{Z=(1+e^{-\Delta/(kT)})^N.}

(2)​

F=−kTln⁡Z=−NkTln⁡(1+e−Δ/(kT)).\boxed{F=-kT\ln Z=-NkT\ln(1+e^{-\Delta/(kT)}).}

(3)​

S=−∂F∂T=Nkln⁡(1+e−Δ/(kT))+NΔ/T1+eΔ/(kT).\boxed{S=-\frac{\partial F}{\partial T} =Nk\ln(1+e^{-\Delta/(kT)})+\frac{N\Delta/T}{1+e^{\Delta/(kT)}}.}

(4)​

S⟶Nkln⁡2=kln⁡(2N).\boxed{S\longrightarrow Nk\ln2=k\ln(2^N).}

二準位の占有確率がともに 1/21/2 に近づき、全 2N2^N 状態が等確率になる。

(5)​

F=U−TSF=U-TS より S=(U−F)/T\boxed{S=(U-F)/T}。

(6)​

U=F+TS=NΔ1+eΔ/(kT).\boxed{U=F+TS=\frac{N\Delta}{1+e^{\Delta/(kT)}}.}

これは励起準位の占有確率に NΔN\Delta を掛けた結果とも一致する。

I-3C​

(1)​

鉛直下向きからの角度を θ\theta とすれば

L=12mR2θ˙2+mgRcos⁡θ,θ¨+gRθ≃0.L=\frac12mR^2\dot\theta^2+mgR\cos\theta, \qquad \ddot\theta+\frac gR\theta\simeq0.

よって T=2πR/g\boxed{T=2\pi\sqrt{R/g}}。

(2)​

軌道半径は ρ=R2−z02\rho=\sqrt{R^2-z_0^2}。拘束力の鉛直成分が重力を支え、水平成分が向心力となるので

mv02/ρmg=ρz0,v0=g(R2−z02)z0.\frac{mv_0^2/\rho}{mg}=\frac\rho{z_0},\qquad \boxed{v_0=\sqrt{\frac{g(R^2-z_0^2)}{z_0}}.}

(3)​

エネルギーと zz 軸周りの角運動量が保存し、

12mv2+mgz=12mv12,mρvφ=mRv1.\frac12m v^2+mgz=\frac12mv_1^2,\qquad m\rho v_\varphi=mRv_1.

最低点では z˙=0\dot z=0 なので速度は方位角方向のみであり、ρ2=R2−z12\rho^2=R^2-z_1^2。したがって

v12−2gz1=R2v12R2−z12,z1(2gz12−v12z1−2gR2)=0.v_1^2-2gz_1=\frac{R^2v_1^2}{R^2-z_1^2},\qquad z_1(2gz_1^2-v_1^2z_1-2gR^2)=0.

負の折返し点を選び、

z1=v12−v14+16g2R24g.\boxed{z_1=\frac{v_1^2-\sqrt{v_1^4+16g^2R^2}}{4g}.}

v1=0v_1=0 の場合も極限として z1=−Rz_1=-R となる。

Reference​