索末菲展開是由阿諾德·索末菲發展的一種近似計算方法,專門用於計算在凝聚體物理和統計物理中出現的一類特定的積分。在物理中,這類積分表示的是採用費米-狄拉克分佈計算的統計平均。
在熱力學beta的值較大的情況下,我們可以把以下形式的積分關於
展開為:[1][2]
![{\displaystyle \int _{-\infty }^{\infty }{\frac {H(\varepsilon )}{e^{\beta (\varepsilon -\mu )}+1}}\,\mathrm {d} \varepsilon =\int _{-\infty }^{\mu }H(\varepsilon )\,\mathrm {d} \varepsilon +{\frac {\pi ^{2}}{6}}\left({\frac {1}{\beta }}\right)^{2}H^{\prime }(\mu )+O\left({\frac {1}{\beta \mu }}\right)^{4}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f983830f3d4215cb10eb0d9cd34ba5df0b96e227)
上式即為索末菲展開的一般形式。其中
表示一個任意函數,
表示
在
處的導數;
指
的n階最小量(參見大O符號),表示此展開中所有不大於
的項。只有當
在
時趨向於零,且
時
的增長速度不快於任意多項式,我們才能運用這個展開做近似計算。
在自由電子模型中的應用
此類積分常常在計算固體的自由電子模型時出現。在這些計算中,上述積分表示的是
的期望值。通過計算這些積分,我們可以進一步確定
和化學勢
的關係。
考慮在一給定空間中相互獨立運動且處於熱力學平衡的一群全同粒子。如果這群粒子滿足鮑利不相容原理,能量為
的單粒子量子態的佔據機率遵循費米-狄拉克分佈:
[3]
假設
為單粒子的狀態密度,N為導帶電子的總數。則
![{\displaystyle N=\int _{-\infty }^{\infty }D(\varepsilon )f(\varepsilon ,\mu ,T)\mathrm {d} \varepsilon }](https://wikimedia.org/api/rest_v1/media/math/render/svg/63957bd009159d34df821ca5559e808560912fdc)
對於任意的狀態密度函數
,我們可能無法直接計算出此積分。但如果我們應用索末菲展開,我們可以得到以下近似結果:
![{\displaystyle N=\int _{-\infty }^{\mu }D(\varepsilon )\,\mathrm {d} \varepsilon +{\frac {\pi ^{2}}{6}}\left({\frac {1}{\beta }}\right)^{2}D^{\prime }(\mu )+O\left({\frac {1}{\beta \mu }}\right)^{4}}](https://wikimedia.org/api/rest_v1/media/math/render/svg/fe72a5e07ed378f656fcb7c7a1d141945703d0ef)
對於自由電子氣體,我們有
。由此經過一系列計算,我們可以得到:
[4]
其中
。
根據此近似計算所需的假設,索末菲展開常被用於對低溫系統的近似計算。
索末菲展開的推導
在此章節中,我們需要對我們研究的積分關於
作二階展開,其中
是溫度和波茲曼常數的乘積。
我們先作變數代換
:
![{\displaystyle I=\int _{-\infty }^{\infty }{\frac {H(\varepsilon )}{e^{\beta (\varepsilon -\mu )}+1}}\,\mathrm {d} \varepsilon =\tau \int _{-\infty }^{\infty }{\frac {H(\mu +\tau x)}{e^{x}+1}}\,\mathrm {d} x\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c85146a017e4bd34d2559c9058efb1cce669ee3f)
將積分範圍劃分成兩部分,
,並對
作變數代換
:
![{\displaystyle I=\underbrace {\tau \int _{-\infty }^{0}{\frac {H(\mu +\tau x)}{e^{x}+1}}\,\mathrm {d} x} _{I_{1}}+\underbrace {\tau \int _{0}^{\infty }{\frac {H(\mu +\tau x)}{e^{x}+1}}\,\mathrm {d} x} _{I_{2}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/4e1a6b9cad594b91700c3c06a1b61f3dc2bee175)
![{\displaystyle I_{1}=\tau \int _{-\infty }^{0}{\frac {H(\mu +\tau x)}{e^{x}+1}}\,\mathrm {d} x=\tau \int _{0}^{\infty }{\frac {H(\mu -\tau x)}{e^{-x}+1}}\,\mathrm {d} x\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c446d6e21bb52ecb8dfb2e6e430fd9fa9231b22d)
接下來,通過使用以下等式
![{\displaystyle {\frac {1}{e^{-x}+1}}=1-{\frac {1}{e^{x}+1}}\,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/f532aa437ec686b81dae0b0727c69bb24f05d489)
可被化為下述形式:
![{\displaystyle I_{1}=\tau \int _{0}^{\infty }H(\mu -\tau x)\,\mathrm {d} x-\tau \int _{0}^{\infty }{\frac {H(\mu -\tau x)}{e^{x}+1}}\,\mathrm {d} x\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/3eb7e3223d027e0929c1de6dd095838638a42591)
再對第一項作變數代換
將
變換回原來的變數。結合
,我們可以得到:
![{\displaystyle I=\int _{-\infty }^{\mu }H(\varepsilon )\,\mathrm {d} \varepsilon +\tau \int _{0}^{\infty }{\frac {H(\mu +\tau x)-H(\mu -\tau x)}{e^{x}+1}}\,\mathrm {d} x\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ce8d10947fbca99e9c5f9c8ba6e12f1ef598d952)
若
足夠小,
足夠平滑,第二項的分子可以被如下近似到第一階導數:
![{\displaystyle \Delta H=H(\mu +\tau x)-H(\mu -\tau x)\approx 2\tau xH'(\mu )+\cdots \,,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/a9fae98dfb0cded51e475ff09ff57d4b9dafdd59)
代入前式可得:
![{\displaystyle I=\int _{-\infty }^{\mu }H(\varepsilon )\,\mathrm {d} \varepsilon +2\tau ^{2}H'(\mu )\int _{0}^{\infty }{\frac {x\mathrm {d} x}{e^{x}+1}}\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/92607fc98eed201cf9ce556b45d6d91b8f29d908)
已知第二項定積分的值為[5] :
.
因此,
![{\displaystyle I=\int _{-\infty }^{\infty }{\frac {H(\varepsilon )}{e^{\beta (\varepsilon -\mu )}+1}}\,\mathrm {d} \varepsilon \approx \int _{-\infty }^{\mu }H(\varepsilon )\,\mathrm {d} \varepsilon +{\frac {\pi ^{2}}{6\beta ^{2}}}H'(\mu )\,}](https://wikimedia.org/api/rest_v1/media/math/render/svg/d8bb47dc5c8caa4b62a99e1cd2c3bb5610d8793b)
母函數
費米分佈的矩的母函數是:
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}e^{\tau \epsilon /2\pi }\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{\tau }}\left\{{\frac {({\frac {\tau T}{2}})}{\sin({\frac {\tau T}{2}})}}e^{\tau \mu /2\pi }-1\right\},\quad 0<\tau T/2\pi <1.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5284e2e04749c0926b29c48be39ab41d589d8964)
這裏,
,且我們通過減去一個單位階躍函數
去掉了溫度為零的情況下發散的函數值。關於
,計算其各次展開後可以得到以下結果:[6]
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}=\left({\frac {\mu }{2\pi }}\right),}](https://wikimedia.org/api/rest_v1/media/math/render/svg/b5f87d5cb601e027145fd5d96538a9f63885a2c7)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}\left({\frac {\epsilon }{2\pi }}\right)\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{2!}}\left({\frac {\mu }{2\pi }}\right)^{2}+{\frac {T^{2}}{4!}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/63889301fa874652fb1e78c1bb3da85c5b63ed97)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}{\frac {1}{2!}}\left({\frac {\epsilon }{2\pi }}\right)^{2}\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{3!}}\left({\frac {\mu }{2\pi }}\right)^{3}+\left({\frac {\mu }{2\pi }}\right){\frac {T^{2}}{4!}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/ab0c3b33b04517dff8a7f98705377e0508771a33)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}{\frac {1}{3!}}\left({\frac {\epsilon }{2\pi }}\right)^{3}\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{4!}}\left({\frac {\mu }{2\pi }}\right)^{4}+{\frac {1}{2!}}\left({\frac {\mu }{2\pi }}\right)^{2}{\frac {T^{2}}{4!}}+{\frac {7}{8}}{\frac {T^{4}}{6!}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/1ffd6b9d1a636fde5c99e13ccc7fa43f5e7495e9)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}{\frac {1}{4!}}\left({\frac {\epsilon }{2\pi }}\right)^{4}\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{5!}}\left({\frac {\mu }{2\pi }}\right)^{5}+{\frac {1}{3!}}\left({\frac {\mu }{2\pi }}\right)^{3}{\frac {T^{2}}{4!}}+\left({\frac {\mu }{2\pi }}\right){\frac {7}{8}}{\frac {T^{4}}{6!}},}](https://wikimedia.org/api/rest_v1/media/math/render/svg/c526b7a329ee5fbe3e7fac2f4f588a303bccf5d8)
![{\displaystyle \int _{-\infty }^{\infty }{\frac {d\epsilon }{2\pi }}{\frac {1}{5!}}\left({\frac {\epsilon }{2\pi }}\right)^{5}\left\{{\frac {1}{1+e^{\beta (\epsilon -\mu )}}}-\theta (-\epsilon )\right\}={\frac {1}{6!}}\left({\frac {\mu }{2\pi }}\right)^{6}+{\frac {1}{4!}}\left({\frac {\mu }{2\pi }}\right)^{4}{\frac {T^{2}}{4!}}+{\frac {1}{2!}}\left({\frac {\mu }{2\pi }}\right)^{2}{\frac {7}{8}}{\frac {T^{4}}{6!}}+{\frac {31}{24}}{\frac {T^{6}}{8!}}.}](https://wikimedia.org/api/rest_v1/media/math/render/svg/5b3dd4499f1d5bbf459edfaa3821fa28d59e4629)
對於玻色函數的奇矩(odd moment),我們有相似的母函數:
註釋
- ^ Ashcroft & Mermin 1976,第760頁 harvnb error: no target: CITEREFAshcroft_&_Mermin1976 (help).
- ^ Fabian, J. "Sommerfeld's expansion" (頁面存檔備份,存於互聯網檔案館) (PDF).
- ^ Grosso 2014,第112頁 harvnb error: no target: CITEREFGrosso2014 (help).
- ^ Grosso 2014,第116頁 harvnb error: no target: CITEREFGrosso2014 (help).
- ^ "Definite integrals containing exponential functions" (頁面存檔備份,存於互聯網檔案館).
- ^ R. Loganayagam, P. Surówka. Anomaly/Transport in an Ideal Weyl gas. JHEP. 2012, 04: 2012:97. Bibcode:2012JHEP...04..097L. arXiv:1201.2812
. doi:10.1007/JHEP04(2012)097.
參考文獻