查看“︁德拜模型”︁的源代码

{{NoteTA|G1=物理學}}
{{roughtranslation}}
[[File:Debye100.jpg|thumb|200px|[[彼得·德拜]]]]
在[[热力学]]和[[固体物理学]]中，'''德拜模型'''（英語：'''Debye model'''）是由[[彼得·德拜]]在1912年提出的方法<ref>{{cite journal |first=Peter |last=Debye |title=Zur Theorie der spezifischen Waerme  |language=de|journal=[[Annalen der Physik]] |volume=39 |issue=4 |pages=789–839 |year=1912 |doi= 10.1002/andp.19123441404 |bibcode = 1912AnP...344..789D }}</ref>，用于估算[[声子]]对固体的[[比热]]（热容）的贡献。德拜模型把[[原子]][[晶格]]的[[振动]]（熱）当作盒中的[[聲子]]处理，而与此不同的[[爱因斯坦模型]]则将固体作为许多单独的、不相互作用的[[量子谐振子]]处理。德拜模型正确地预言了低温时固体的热容：与<math>T^3</math>成正比（德拜<math>T^3</math>定律）。与[[爱因斯坦模型]]一样，它在高温时也与[[杜隆-珀蒂定律]]相符合。但由于模型的简单假设，它在中间的温度不太准确。

==推导==
德拜模型是[[普朗克黑体辐射定律]]的固体版本，后者把[[电磁辐射]]视为[[盒中气体|盒中的光子气体]]，而德拜模型把原子的振动视为盒中的[[声子]]（盒子就是这块固体）。它们大部分的计算步骤都是相同的，因为[[声子]]气和光子气都是有线性色散关系的[[玻色气体]]。

考虑一个边长为<math>L</math>的立方体。从[[盒中粒子]]一文可知，盒中的（沿某一边长方向传播的）声波由于共振效应，其波长只能为：

:<math>\lambda_n = {2L\over n}\,,</math>

其中<math>n</math>是整数。一个声子的能量是：

:<math>E_n\ =h\nu_n\,,</math>

其中<math>h</math>是[[普朗克常数]]，<math>\nu_{n}</math>是声子的频率。近似认为频率与波长成反比，则有：

:<math>E_n=h\nu_n={hc_s\over\lambda_n}={hc_sn\over 2L}\,,</math>

其中<math>c_s</math>是固体中的声速。推广到三维的情况（沿任意方向传播的[[声子]]）：

:<math>E_n^2={p_n}^2{c_s}^2=\left({hc_s\over2L}\right)^2\left(n_x^2+n_y^2+n_z^2\right)\,.</math>

其中<math>p_n</math>为声子的动量。频率与波长成反比的近似（即声速恒定）对于低能量声子是准确的，但对于高能量声子则不然（参见[[声子]]）。这是德拜模型的局限之一，导致其对中间温度的预言不正确，而对低温和高温都很精确。

下面来计算盒中的总能量：

:<math>U = \sum_n E_n\,\bar{N}(E_n)\,,</math>

其中<math>\bar{N}(E_n)</math>是盒中能量为<math>E_n</math>的声子的数目。换句话说，总能量等于某能量的声子数目乘以其能量再对各能量求和。在三维空间中，我们有：

:<math>U = \sum_{n_x}\sum_{n_y}\sum_{n_z}E_n\,\bar{N}(E_n)\,.</math>

这是德拜模型和[[普朗克黑体辐射定律]]的不同之处。与盒中的电磁辐射不一样，由于声子不能有无限大的频率，所以盒中声子只有有限个能量状态。它的频率由它的传播介质——固体的原子晶格所约束。考虑以下的横向声子的插图：

:::::[[Image:Phonons modes.jpg|400px]]

可以合理假设声子的最小波长是原子间距的两倍，如最下一图所示。固体中有<math>N</math>个原子，我们的固体是正方体，因此每一条边有<math>\sqrt[3]{N}</math>个原子，原子间距为<math>L/\sqrt[3]{N}</math>，从而最小波长为：

:<math>\lambda_{\rm min} = {2L \over \sqrt[3]{N}}\,,</math>

最大的模数<math>n</math>（对于[[光子]]是无限大）为：

:<math>n_{\rm max} = \sqrt[3]{N}\,.</math>

这设置了三重能量求和的上限：

:<math>U = \sum_{n_x}^{\sqrt[3]{N}}\sum_{n_y}^{\sqrt[3]{N}}\sum_{n_z}^{\sqrt[3]{N}}E_n\,\bar{N}(E_n)\,.</math>

对于缓慢变化的、表现良好的函数，求和可以用积分来代替（又称为[[托马斯-费米近似]]）：

:<math>U \approx\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}} E(n)\,\bar{N}\left(E(n)\right)\,dn_x\, dn_y\, dn_z\,.</math>

至此，我们还没有提到<math>\bar{N}(E)</math>：能量为<math>E\,</math>的声子的数目。声子服从[[玻色-爱因斯坦统计]]。它们的分布由著名的玻色-爱因斯坦公式给出：

:<math>\langle N\rangle_{BE} = {1\over e^{E/kT}-1}\,.</math>

由于一个声子有三个可能的偏振态（一个[[纵波|纵向]]、两个[[横波|横向]]，大致携带相同的能量），需把以上公式乘以3：

:<math>\bar{N}(E) = {3\over e^{E/kT}-1}\,.</math>

实际上我们可以定义''等效声速''<math>c_s:=c_{{\rm eff}}</math>，也就是说，德拜温度<math>T_d</math>（见下文）与<math>c_{{\rm eff}}</math>成正比。更加精确地讲，<math>T_D^{-3}\propto c_{{\rm eff}}^{-3}:=(1/3)c_{{\rm long}}^{-3}+(2/3)c_{{\rm trans}}^{-3}</math>，其中我们区分了纵向和横向的声波速度（贡献分别为1/3和2/3）。德拜温度或有效声速是晶体的硬度的一种衡量。

把此式代入能量积分，得：

:<math>U  = \int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}}\int_0^{\sqrt[3]{N}} E(n)\,{3\over e^{E(n)/kT}-1}\,dn_x\, dn_y\, dn_z\,.</math>

这比光子的情况容易得多，因为（在半经典式的处理中）光子的频率（也即这个三重积分）是无上界的。但如上面图所示，这一点对于[[声子]]不成立。为了近似计算这个三重积分，德拜使用了[[球坐标系]]：

:<math>\ (n_x,n_y,n_z)=(n\cos \theta \cos \phi,n\cos \theta \sin \phi,n\sin \theta )</math>

并用八分之一球（球在一个卦限内的部分）来近似代替立方体：

:<math>U \approx\int_0^{\pi/2}\int_0^{\pi/2}\int_0^R E(n)\,{3\over e^{E(n)/kT}-1}n^2 \sin\theta\, dn\, d\theta\, d\phi\,,</math>

其中<math>R</math>是球的半径，通过保持立方体和八分之一球中的粒子数目相同来得出。立方体的体积是<math>N</math>： 

:<math>N = {1\over8}{4\over3}\pi R^3\,,</math>

因此：

:<math>R = \sqrt[3]{6N\over\pi}\,.</math>

用球的积分来代替立方体积分这一近似是德拜模型不准确性的另一来源。

先积分掉角向部分：

:<math>U = {3\pi\over2}\int_0^R \,{hc_sn\over 2L}{n^2\over e^{hc_sn/2LkT}-1} \,dn</math>

再利用变量代换<math>x = {hc_sn\over 2LkT}</math>：

:<math>U = {3\pi\over2} kT \left({2LkT\over hc_s}\right)^3\int_0^{hc_sR/2LkT} {x^3\over e^x-1}\, dx</math>

为简化表达式，可以定义'''德拜温度'''<math>T_D</math>——它的量纲与温度相同，因物质而异：

:<math>T_D\ \stackrel{\mathrm{def}}{=}\  {hc_sR\over2Lk} = {hc_s\over2Lk}\sqrt[3]{6N\over\pi} = {hc_s\over2k}\sqrt[3]{{6\over\pi}{N\over V}}</math>

于是比内能为：

:<math>\frac{U}{Nk} = 9T \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^3\over e^x-1}\, dx = 3T D_3 \left({T_D\over T}\right)\,,</math>

其中<math>D_3(x)</math>是（第三）[[德拜函数]]。

对<math>T</math>微分，我们便得到无量纲热容：

:<math> \frac{C_V}{Nk} = 9 \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^4 e^x\over\left(e^x-1\right)^2}\, dx\,.</math>

这个公式给出了任何温度下德拜模型的结果。下面更简单的公式给出了低温和高温极限下的渐近表现。前面已经提到，这两个渐进表现是精确的，而中间温度的表现则不然。其根本原因在于，低温和高温近似下，德拜模型都给出了其一，低频下精确的{{zh-em|[[色散关系]]}}<math>E(\nu )</math>，其二，精确的{{zh-em|[[状态密度]]}}（单位频率间隔内的声子数目）<math>(\int g(\nu ) \, {\rm d\nu}\equiv 3N)</math>。

==德拜的推导==

实际上，德拜用不同且更简单的方法推出了这个公式。利用[[连续介质力学|连续介质]]的[[固体力学]]，他发现频率小于某个特定值的振动状态的数目趋近于：

:<math> n \sim {1 \over 3} \nu^3 V F\,, </math>

其中<math> V </math>是体积，<math> F </math>是一个从[[弹性 (物理学)|弹性系数]]和密度计算出的因子。把这与温度T的量子谐振子的期望能量（已经由爱因斯坦在他的模型中使用）结合，便给出能量：

:<math>U = \int_0^\infty \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,</math>

如果振动频率趋于无穷大。这个形式给出了<math>T^3</math>的表现，它在低温时是正确的。但德拜意识到N个原子不可能有超过<math>3N</math>个振动状态。他假设在原子固体中，振动状态的频谱将继续遵循以上的规则，直到一个最大频率<math>\nu_m</math>使得总的状态数目为<math>3N</math>为止：

:<math> 3N = {1 \over 3} \nu_m^3 V F \,.</math>

德拜知道这个假设不是真正正确的（较高的频率比假设要更加密集），但它保证了高温时的正确表现（[[杜隆-珀蒂定律]]）。于是，能量由以下给出：

:<math>U = \int_0^{\nu_m} \,{h\nu^3 V F\over e^{h\nu/kT}-1}\, d\nu\,,</math>

::<math> = V F kT (kT/h)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,</math>

其中<math>T_D</math>是<math>h\nu_m/k</math>。

::<math> = 9 N k T (T/T_D)^3 \int_0^{T_D/T} \,{x^3 \over e^x-1}\, dx\,,</math>

::<math> = 3 N k T D_3(T_D/T)\,,</math>

其中<math>D_3</math>是一个函数，后来命名为三阶[[德拜函数]]。

==低温极限==
[[Image:DebyeVSEinstein.jpg|thumb|比較德拜、愛因斯坦分別對於熱容與溫度之間關係的預測。|200px]]
德拜固体的温度称为低的，如果<math>T \ll T_D</math>，在这个情况下：

:<math> \frac{C_V}{Nk} \sim 9 \left({T\over T_D}\right)^3\int_0^{\infty} {x^4 e^x\over \left(e^x-1\right)^2}\, dx</math>

这个定积分可以精确计算：

:<math> \frac{C_V}{Nk} \sim {12\pi^4\over5} \left({T\over T_D}\right)^3</math>

在低温极限中，德拜模型的局限不适用，它给出了（声子）热容、温度、弹性系数，以及每个原子的体积（后面的数量是包含在德拜温度之中的）之间的正确关系。

==高温极限==

德拜固体的温度称为高的，如果<math>T \gg T_D</math>。<math>e^x - 1\approx  x</math>如果<math>|x|<<1</math>，得出：

:<math> \frac{C_V}{Nk} \sim 9 \left({T\over T_D}\right)^3\int_0^{T_D/T} {x^4 \over x^2}\, dx </math>

:<math>\frac{C_V}{Nk} \sim 3\,.</math>

这就是[[杜隆-珀蒂定律]]，它是相当准确的，虽然它没有考虑非谐性，这造成了热容进一步上升。如果固体是[[导体]]或[[半导体]]，那么它的总热容还可能含有电子的不可忽略的贡献。

==参见==
*[[玻色气体]]
*[[盒中气体]]

==参考文献==
<references />

==外部链接==
*[http://www.thphys.de/F14_Debye_Experiment.html 用低温恒温器来决定石英的比热、热导率和电导率] {{Wayback|url=http://www.thphys.de/F14_Debye_Experiment.html |date=20180824135105 }}

{{固体物理学}}
{{统计力学}}

{{DEFAULTSORT:Debye model}}
[[Category:热力学]]
[[Category:统计力学]]
[[Category:固体物理学]]
[[Category:荷兰发明]]