跳到主要内容

京都大学 情報学研究科 数理工学専攻 2022年度 専門科目 問題5(物理統計学)

Author​

祭音Myyura (Based on Miyake's answer refined with GPT 6 Astra)

Description​

エネルギーレベルが

En=hν(12+n),n=0,1,2,3,…E_n=h\nu\left(\frac12+n\right),\qquad n=0,1,2,3,\ldots

なる振動数 ν>0\nu>0 の振動子系を考える。ここで h>0h>0 は定数であり、エネルギーレベルの縮退は無く、同系の分配関数 ZZ は

Z=∑n=0∞exp⁡(−EnkT)Z=\sum_{n=0}^{\infty}\exp\left(-\frac{E_n}{kT}\right)

で与えられるとする。ただし、k>0k>0 をボルツマン定数、TT を絶対温度とする。以下の問いに答えよ。

(i) 分配関数 ZZ を計算せよ。

(ii) エネルギー EE の期待値 ⟨E⟩\langle E\rangle を求めよ。

(iii) 比熱 C=d⟨E⟩/dTC=d\langle E\rangle/dT を求めよ。

(iv) 比熱 CC の低温極限 T→0T\to0 を求めよ。

(v) 比熱 CC の高温極限 T→∞T\to\infty を求めよ。

题目描述​

频率为 ν>0\nu>0 的振子具有无简并能级

En=hν(12+n),n=0,1,2,…,E_n=h\nu\left(\frac12+n\right),\qquad n=0,1,2,\ldots,

其中 h>0h>0 为常数。其配分函数为

Z=∑n=0∞e−En/(kT),Z=\sum_{n=0}^{\infty}e^{-E_n/(kT)},

k>0k>0 为玻尔兹曼常数,TT 为绝对温度。

(i) 求配分函数 ZZ。

(ii) 求能量的期望值 ⟨E⟩\langle E\rangle。

(iii) 求比热 C=d⟨E⟩/dTC=d\langle E\rangle/dT。

(iv) 求比热在 T→0T\to0 时的低温极限。

(v) 求比热在 T→∞T\to\infty 时的高温极限。

Kai​

以下では β=1/(kT)\beta=1/(kT)、ε=hν\varepsilon=h\nu とおく。

(i)​

T>0T>0 において 0<e−βε<10<e^{-\beta\varepsilon}<1 なので、等比級数を用いると

Z=e−βε/2∑n=0∞e−nβε=e−βε/21−e−βε=12sinh⁡(βε/2).Z=e^{-\beta\varepsilon/2}\sum_{n=0}^{\infty}e^{-n\beta\varepsilon} =\boxed{\frac{e^{-\beta\varepsilon/2}}{1-e^{-\beta\varepsilon}} =\frac{1}{2\sinh(\beta\varepsilon/2)}}.

(ii)​

正準分布の関係式から

⟨E⟩=−∂log⁡Z∂β=ε2+εeβε−1=hν2coth⁡(hν2kT).\langle E\rangle=-\frac{\partial\log Z}{\partial\beta} =\boxed{\frac{\varepsilon}{2}+\frac{\varepsilon}{e^{\beta\varepsilon}-1}} =\frac{h\nu}{2}\coth\left(\frac{h\nu}{2kT}\right).

(iii)​

ε\varepsilon は温度に依存しないため、(ii) を TT で微分して

C=(hν)2kT2ehν/(kT)(ehν/(kT)−1)2.\boxed{C=\frac{(h\nu)^2}{kT^2} \frac{e^{h\nu/(kT)}}{\left(e^{h\nu/(kT)}-1\right)^2}}.

(iv)​

x=hν/(kT)x=h\nu/(kT) とおくと

C=kx2ex(ex−1)2.C=kx^2\frac{e^x}{(e^x-1)^2}.

T→0+T\to0^+ では x→∞x\to\infty であり、C∼kx2e−xC\sim kx^2e^{-x} なので

lim⁡T→0+C=0.\boxed{\lim_{T\to0^+}C=0}.

(v)​

T→∞T\to\infty では x→0x\to0 である。(ex−1)/x→1(e^x-1)/x\to1 を用いて

lim⁡T→∞C=klim⁡x→0ex((ex−1)/x)2=k.\boxed{\lim_{T\to\infty}C =k\lim_{x\to0}\frac{e^x}{\left((e^x-1)/x\right)^2}=k}.

Reference​