跳到主要内容

京都大学 理学研究科 物理学・宇宙物理学専攻 2020年度 問題II-3-AB

Author​

祭音Myyura (refined with GPT 6 Astra)

Description​

II-3A​

Re⁡z>0\operatorname{Re}z>0 に対し Γ(z)=∫0∞tz−1e−tdt\Gamma(z)=\int_0^\infty t^{z-1}e^{-t}dt と定義する。

  1. 自然数 NN に対して Γ(N+1)=N!\Gamma(N+1)=N! を示せ。
  2. Γ(1/2)\Gamma(1/2) を求めよ。ガウス積分を用いてよい。
  3. 大きな NN に対して log⁡N!=N(log⁡N−1)+Δ\log N!=N(\log N-1)+\Delta、Δ≪N\Delta\ll N と書けることを示し、Δ\Delta の主要項を求めよ。必要なら、最大点 t0t_0 で f′′(t0)<0f''(t_0)<0 となる滑らかな関数について、m≫1m\gg1 のとき ∫emf(t)dt≃emf(t0)−2π/[mf′′(t0)]\int e^{mf(t)}dt\simeq e^{mf(t_0)}\sqrt{-2\pi/[mf''(t_0)]} を用いてよい。
  4. α>0\alpha>0 のとき ∫0∞tα/(et−1)dt\int_0^\infty t^\alpha/(e^t-1)dt を Γ\Gamma と ζ(s)=∑n=1∞n−s\zeta(s)=\sum_{n=1}^\infty n^{-s} で表せ。

II-3B​

運動エネルギー K=5 MeVK=5\,\mathrm{MeV}、電荷 2e2e の α\alpha 粒子を金薄膜に入射する。金原子核の電荷は 79e79e、質量数は197であり、α\alpha 粒子に比べ十分重いとする。古典的クーロン散乱として扱い、図1の衝突係数を bb、散乱角を θ\theta、最接近距離を DD とする。mec2=0.50 MeVm_ec^2=0.50\,\mathrm{MeV}、古典電子半径 e2/(4πε0mec2)=2.8×10−15 me^2/(4\pi\varepsilon_0m_ec^2)=2.8\times10^{-15}\,\mathrm m を用いてよい。

原題 図1:金原子核によるアルファ粒子の散乱

  1. 正面衝突での最接近距離 Dmin⁡D_{\min} を求めよ。

  2. 単位面積・単位時間あたり nn 個の粒子を入射する。衝突係数が bb と b+dbb+db の間にある粒子数率 dndn を求めよ。

  3. cot⁡(θ/2)=2b/Dmin⁡\cot(\theta/2)=2b/D_{\min} を使い、θ\theta と θ+dθ\theta+d\theta に散乱される粒子数率と微分散乱断面積 dσ/dΩd\sigma/d\Omega を求めよ。

  4. 入射粒子束を NaN_a、薄膜の単位面積あたり原子数を NbN_b とする。多重散乱を無視し、図2のように散乱角 60∘60^\circ、距離 rr の場所に検出面積 SS、効率100%の計数器を向けたところ、単位時間に NcN_c 個を検出した。rr は薄膜の大きさより十分大きい。dσ(60∘)/dΩd\sigma(60^\circ)/d\Omega を求めよ。

    原題 図2:金薄膜と粒子計数器

  5. (4) の測定値から Dmin⁡D_{\min} を求めよ。

题目描述​

A:由伽马函数积分定义证明阶乘关系,求半整数值,以鞍点近似求 log⁡N!\log N! 中的主要修正,并用伽马函数和黎曼 zeta 函数表示玻色型积分。

B:如上方图1,5 MeV 的 α\alpha 粒子对重金核发生库仑散射,电荷分别为 2e,79e2e,79e。(1) 计算正碰最接近距离;(2) 求冲量参数环带的粒子数率;(3) 利用 cot⁡(θ/2)=2b/Dmin⁡\cot(\theta/2)=2b/D_{\min} 推导 Rutherford 微分截面;(4) 按上方图2配置,入射通量 NaN_a、靶面数密度 NbN_b,距离 rr、面积 SS、方向 60∘60^\circ 的理想探测器计数率为 NcN_c,求微分截面;(5) 反求 Dmin⁡D_{\min}。

Kai​

II-3A​

(1) 部分積分により Γ(z+1)=zΓ(z)\Gamma(z+1)=z\Gamma(z)。Γ(1)=1\Gamma(1)=1 と合わせて Γ(N+1)=N!\boxed{\Gamma(N+1)=N!}。

(2) t=x2t=x^2 と置けば Γ(1/2)=2∫0∞e−x2dx=π\Gamma(1/2)=2\int_0^\infty e^{-x^2}dx=\boxed{\sqrt\pi}。

(3) Γ(N+1)=∫0∞e−t+Nlog⁡tdt\Gamma(N+1)=\int_0^\infty e^{-t+N\log t}dt の指数は t=Nt=N で最大となる。二次まで展開して

−t+Nlog⁡t=Nlog⁡N−N−(t−N)22N+⋯ .-t+N\log t=N\log N-N-\frac{(t-N)^2}{2N}+\cdots.

従って N!≃NNe−N2πNN!\simeq N^Ne^{-N}\sqrt{2\pi N} であり

Δ=12log⁡(2πN)+O(N−1).\boxed{\Delta=\frac12\log(2\pi N)+O(N^{-1}).}

これは NN に比べ十分小さい。

(4) 正の項の級数 1/(et−1)=∑n=1∞e−nt1/(e^t-1)=\sum_{n=1}^\infty e^{-nt} を積分して

∫0∞tαet−1dt=∑n=1∞n−α−1∫0∞xαe−xdx=Γ(α+1)ζ(α+1).\int_0^\infty\frac{t^\alpha}{e^t-1}dt =\sum_{n=1}^\infty n^{-\alpha-1}\int_0^\infty x^\alpha e^{-x}dx =\boxed{\Gamma(\alpha+1)\zeta(\alpha+1)}.

α>0\alpha>0 により原点近傍の積分も収束する。

II-3B​

(1) エネルギー保存から K=(2⋅79)e2/(4πε0Dmin⁡)K=(2\cdot79)e^2/(4\pi\varepsilon_0D_{\min})。よって

Dmin⁡=158mec2Kre=1580.505(2.8×10−15)≃4.4×10−14 m.\boxed{D_{\min}=158\frac{m_ec^2}{K}r_e =158\frac{0.50}{5}(2.8\times10^{-15}) \simeq4.4\times10^{-14}\,\mathrm m.}

(2) 環帯の面積を掛けて dn=2πnb db\boxed{dn=2\pi n b\,db}。

(3) b=(Dmin⁡/2)cot⁡(θ/2)b=(D_{\min}/2)\cot(\theta/2)、∣db/dθ∣=(Dmin⁡/4)csc⁡2(θ/2)|db/d\theta|=(D_{\min}/4)\csc^2(\theta/2)。従って

dn(θ)=πnDmin⁡24cot⁡θ2csc⁡2θ2 dθ,\boxed{dn(\theta)=\frac{\pi nD_{\min}^2}{4}\cot\frac\theta2\csc^2\frac\theta2\,d\theta,}
dσdΩ=dnn 2πsin⁡θ dθ=Dmin⁡216sin⁡4(θ/2).\boxed{\frac{d\sigma}{d\Omega}=\frac{dn}{n\,2\pi\sin\theta\,d\theta} =\frac{D_{\min}^2}{16\sin^4(\theta/2)}.}

(4) 検出器の立体角は S/r2S/r^2。照射面積を AirrA_{\rm irr} とする。NaN_a が単位面積・単位時間あたりの入射粒子数であるため、全入射粒子数率は NaAirrN_aA_{\rm irr} であり

Nc=NaAirrNbdσdΩ∣60∘Sr2.N_c=N_aA_{\rm irr}N_b\left.\frac{d\sigma}{d\Omega}\right|_{60^\circ}\frac S{r^2}.

従って

dσdΩ∣60∘=Ncr2NaAirrNbS.\boxed{\left.\frac{d\sigma}{d\Omega}\right|_{60^\circ} =\frac{N_cr^2}{N_aA_{\rm irr}N_bS}.}

NcN_c を総計数率として用いる場合、照射面積の情報が必要である。一方、NaN_a を総入射粒子数率として定義する、または NcN_c を単位照射面積あたりの検出計数率として扱う規約なら、同じ式は Ncr2/(NaNbS)N_cr^2/(N_aN_bS) と書ける。

(5) 16sin⁡430∘=116\sin^4 30^\circ=1 より

Dmin⁡=rNcNaAirrNbS.\boxed{D_{\min}=r\sqrt{\frac{N_c}{N_aA_{\rm irr}N_bS}}.}

(4) の面積を吸収した規約では Dmin⁡=rNc/(NaNbS)D_{\min}=r\sqrt{N_c/(N_aN_bS)} となる。

Reference​